Priors and Posteriors in Bayesian Timing of Divergence Analyses: The Age of Butterflies Revisited

Deciphering the patterns and timing of diversification of the genus Melanastera (Hemiptera: Psylloidea: Liviidae) in the Neotropics

Even after decades of research on diversification in the Neotropics, our understanding of the evolutionary processes that shaped Neotropical clades is still incomplete. In the present study, we used different divergence times and likelihood-based methods to investigate the influence of biogeography and host plant associations on the diversification of the most species-rich psyllid genus Melanastera (Liviidae) from the Neotropics as a model group of herbivorous insects. We used molecular phylogenetic data from seven gene fragments (four mitochondrial and three nuclear). The putatively monophyletic group of Neotropical Melanastera species has an estimated crown node age of 20.2 Ma (ML, CI 20.2-30.6) or 23.2 Ma (BI, 95 % HPD 16.6-32.6), with diversification occurring mainly in the Upper Miocene, although some species groups diversified in the Pliocene or Pleistocene. Biogeographic analysis suggests that the Neotropical Melanastera originated from the Pacific region of South and Central America. We detected a shift in diversification rates that likely occurred either at the time of origin of Melanastera or during the main colonisation of the Atlantic and Amazon Forests, followed by a subsequent slowdown in speciation rates. State-dependent speciation and extinction models revealed a significant relationship between this diversification shift and the shift of Melanastera to the plant families Melastomataceae and Annonaceae, reflecting the impact of host switching on speciation rates in this group. This period also coincides with several independent dispersal events from the Atlantic and Amazon Forests to other parts of the Neotropics. Taken together, the results of the current study suggest that diversification of Melanastera was facilitated by shifts to new host families, which may have promoted the dispersal of Melanastera into new adaptive zones with subsequent processes of local speciation.

Characterization of the Complete Mitochondrial Genome of Three Satyrid Butterfly Species (Satyrinae:Amathusiini) and Reconstructed Phylogeny of Satyrinae

Satyrinae, one of the most species-rich groups within the Nymphalidae family, has traditionally relied on morphological characteristics for classification. However, this approach encounters challenges due to issues such as cryptic species and paraphyletic groups. Recent molecular phylogenetic studies have revealed the complex evolutionary history of Satyrinae, leading to the reclassification of the originally polyphyletic Satyrini into multiple independent tribes and confirming the monophyletic status of groups such as Amathusiini. Nevertheless, the phylogenetic relationships and divergence times of certain tribes remain contentious. This study focuses on three species of the Amathusiini tribe (Faunis aerope, Stichophthalma howqua, and Aemona lena), constructing a phylogenetic tree by sequencing the complete mitochondrial genome and integrating 13 protein-coding genes, including COI and ND5. The results indicate that the mitogenome lengths for the three satyrid species are 15,512 bp for Faunis aerope, 13,914 bp for Stichophthalma howqua, and 15,288 bp for Aemona lena. The genetic composition and sequencing of the newly obtained mitogenomes exhibit high conservation and are distinctive to this group of butterflies. Each of the three mitogenomes contains a characteristic collection of 37 genes along with an AT-rich region. Notably, the tRNA genes across these mitogenomes display a conventional cloverleaf configuration; however, the tRNASer stem (AGN) lacks the dihydrouridine (DHU) arm. The three species exhibit varying lengths of AT-rich regions, resulting in differences in their mitochondrial genome sizes. Finally, the phylogenetic analysis supports the relationships among the four tribes of Satyrinae as: (Satyrini + (Amathusiini + Elymniini)) + Melanitini.

Rapid Evolution of Host Repertoire and Geographic Range in a Young and Diverse Genus of Montane Butterflies

Evolutionary changes in geographic distribution and larval host plants may promote the rapid diversification of montane insects, but this scenario has been rarely investigated. We studied the rapid radiation of the butterfly genus Colias, which has diversified in mountain ecosystems in Eurasia, Africa, and the Americas. Based on a data set of 150 nuclear protein-coding genetic loci and mitochondrial genomes, we constructed a time-calibrated phylogenetic tree of Colias species with broad taxon sampling. We then inferred their ancestral geographic ranges, historical diversification rates, and the evolution of host use. We found that the most recent common ancestor of Colias was likely geographically widespread and originated ~3.5 Ma. The group subsequently diversified in different regions across the world, often in tandem with geographic expansion events. No aspect of elevation was found to have a direct effect on diversification. The genus underwent a burst of diversification soon after the divergence of the Neotropical lineage, followed by an exponential decline in diversification rate toward the present. The ancestral host repertoire included the legume genera Astragalus and Trifolium but later expanded to include a wide range of Fabaceae genera and plants in more distantly related families, punctuated with periods of host range expansion and contraction. We suggest that the widespread distribution of the ancestor of all extant Colias lineages set the stage for diversification by isolation of populations that locally adapted to the various different environments they encountered, including different host plants. In this scenario, elevation is not the main driver but might have accelerated diversification by isolating populations.

Tempo and mode of winter diapause evolution in butterflies

Quantifying the tempo and mode via modern phylogenetic comparative methods can provide key insights into how selection and constraints shape trait evolution on a macroevolutionary time scale. Here, we elucidate the evolution of hibernation (winter) diapause, a complex and defining life-history trait that allows temporal escape from harsh winters in temperate regions for many insects, including our model system, butterflies. Butterflies can diapause in all major life stages, and the availability of global-scale phylogenies makes them an ideal model system for studying diapause evolution. First, using a thorough literature survey, we scored the developmental stage of hibernation diapause (egg, larva, pupa, adult) vs. absence of diapause. We find that larval diapause is most common, while pupal, egg, and adult diapause are relatively rare. Next, we determined that the loss of diapause occurred at a much higher rate and that gains primarily occurred from the non-diapause state. While ancestral state estimation at deeper nodes remained uncertain, we found consistent patterns for some families and strong evidence for extensive convergence in diapause evolution. Contrary to expectations, we find no support for increased gain of diapause during the Eocene-Oligocene glaciation (~35 million years ago). Overall, the evolution of diapause in butterflies has a complex history, has evolved convergently, and has likely predated the major glaciation event consistent with the deep history of diapause evolution in insects. This study advances our understanding of the evolution of a complex and important life-history trait and establishes a macroevolutionary foundation for future studies on the ultimate and proximate basis of diapause evolution.

Genome assembly of an endemic butterfly (Minois Aurata) shed light on the genetic mechanisms underlying ecological adaptation to arid valley habitat

BACKGROUND

The Hengduan Mountains, one of the global biodiversity hotspots with exceptional species richness and high endemism, contains numerous arid valleys that create a distinctive geographical and ecological landscape. However, the adaptive evolutionary mechanisms of organism in the arid valley remain poorly understood. Minois aurata, an endemic butterfly species found exclusively in the arid valley of the upper Minjiang River, represents an attractive model system for studying adaptive evolutionary mechanisms to arid valley environments.

RESULTS

Here, we present the first chromosome-level genome assembly for Minois aurata, with a total size of approximately 609.17 Mb, and a scaffold N50 size of 23.88 Mb. These scaffolds were further clustered and anchored onto 29 chromosomes based on Hi-C data. A total of 16,163 protein-coding genes were predicted, of which 91.83% were functionally annotated. The expansion of transposable elements (TEs) accounts for the relatively large genome size of M. aurata, potentially aiding its adaptation to environmental conditions. Phylogenomic analyses based on 3,785 single-copy genes revealed that M. aurata is most closely related to Hipparchia semele. Further mitochondrial genome analysis of four Minois species placed M. aurata in a basal position within the genus, supporting it as an independent species. A total of 185 rapidly evolving and 232 specific gene families were identified in M. aurata. Functional enrichment analysis indicated that these gene families were mainly associated with ultraviolet radiation, heat and hypoxia responses. We also identified 234 positive selected genes in M. aurata, some of which are related to compound eye photoreceptor development, osmotic stress, and light stimulus response. Demographic analysis indicated that the effective population size of M. aurata decreased around 0.4 and 0.04 million years ago, respectively, coinciding with the localized sub-glaciation.

CONCLUSION

The chromosome-level genome offers a comprehensive genomic basis for understanding the evolutionary and adaptive strategies of Minois aurata in the unique arid valley environment of the Hengduan Mountains, while also providing valuable insights into the broader mechanisms of organism adaptation to such habitats.

Whole genomes from the extinct Xerces Blue butterfly can help identify declining insect species

The Xerces Blue (Glaucopsyche xerces) is considered to be the first butterfly to become extinct in historical times. It was notable for its chalky lavender wings with conspicuous white spots on the ventral wings. The last individuals were collected in their restricted habitat, in the dunes near the Presidio military base in San Francisco, in 1941. We sequenced the genomes of four 80- to 100-year-old Xerces Blue, and seven historical and one modern specimens of its closest relative, the Silvery Blue (Glaucopsyche lygdamus). We compared these to a novel annotated genome of the Green-Underside Blue (Glaucopsyche alexis). Phylogenetic relationships inferred from complete mitochondrial genomes indicate that Xerces Blue was a distinct species that diverged from the Silvery Blue lineage at least 850,000 years ago. Using nuclear genomes, both species experienced population growth during the Eemian interglacial period, but the Xerces Blue decreased to a very low effective population size subsequently, a trend opposite to that observed in the Silvery Blue. Runs of homozygosity and deleterious load in the former were significantly greater than in the later, suggesting a higher incidence of inbreeding. These signals of population decline observed in Xerces Blue could be used to identify and monitor other insects threatened by human activities, whose extinction patterns are still not well known.

Genome-Wide Gene Birth-Death Dynamics Are Associated with Diet Breadth Variation in Lepidoptera

Comparative analyses of gene birth-death dynamics have the potential to reveal gene families that played an important role in the evolution of morphological, behavioral, or physiological variation. Here, we used whole genomes of 30 species of butterflies and moths to identify gene birth-death dynamics among the Lepidoptera that are associated with specialist or generalist feeding strategies. Our work advances this field using a uniform set of annotated proteins for all genomes, investigating associations while correcting for phylogeny, and assessing all gene families rather than a priori subsets. We discovered that the sizes of several important gene families (e.g. those associated with pesticide resistance, xenobiotic detoxification, and/or protein digestion) are significantly correlated with diet breadth. We also found 22 gene families showing significant shifts in gene birth-death dynamics at the butterfly (Papilionoidea) crown node, the most notable of which was a family of pheromone receptors that underwent a contraction potentially linked with a shift to visual-based mate recognition. Our findings highlight the importance of uniform annotations, phylogenetic corrections, and unbiased gene family analyses in generating a list of candidate genes that warrant further exploration.

Phylogenomic analysis of Stylops reveals the evolutionary history of a Holarctic Strepsiptera radiation parasitizing wild bees

Holarctic Stylops is the largest genus of the enigmatic insect order Strepsiptera, twisted winged parasites. Members of Stylops are obligate endoparasites of Andrena mining bees and exhibit extreme sexual dimorphism typical of Strepsiptera. So far, molecular studies on Stylops have focused on questions on species delimitation. Here, we utilize the power of whole genome sequencing to infer the phylogeny of this morphologically challenging genus from thousands of loci. We use a species tree method, concatenated maximum likelihood analysis and Bayesian analysis with a relaxed clock model to reconstruct the phylogeny of 46 Stylops species, estimate divergence times, evaluate topological consistency across methods and infer the root position. Furthermore, the biogeographical history and coevolutionary patterns with host species are assessed. All methods recovered a well resolved topology with close to all nodes maximally supported and only a handful of minor topological variations. Based on the result, we find that included species can be divided into 12 species groups, seven of them including only Palaearctic species, three Nearctic and two were geographically mixed. We find a strongly supported root position between a clade formed by the spreta, thwaitesi and gwynanae species groups and the remaining species and that the sister group of Stylops is Eurystylops or Eurystylops + Kinzelbachus. Our results indicate that Stylops originated in the Western Palaearctic or Western Palaearctic and Nearctic in the early Neogene or late Paleogene, with four independent dispersal events to the Nearctic. Cophylogenetic analyses indicate that the diversification of Stylops has been shaped by both significant coevolution with the mining bee hosts and host-shifting. The well resolved and strongly supported phylogeny will provide a valuable phylogenetic basis for further studies into the fascinating world of Strepsipterans.

Two Notorious Nodes: A Critical Examination of Relaxed Molecular Clock Age Estimates of the Bilaterian Animals and Placental Mammals

The popularity of relaxed clock Bayesian inference of clade origin timings has generated several recent publications with focal results considerably older than the fossils of the clades in question. Here, we critically examine two such clades: the animals (with a focus on the bilaterians) and the mammals (with a focus on the placentals). Each example displays a set of characteristic pathologies which, although much commented on, are rarely corrected for. We conclude that in neither case does the molecular clock analysis provide any evidence for an origin of the clade deeper than what is suggested by the fossil record. In addition, both these clades have other features (including, in the case of the placental mammals, proximity to a large mass extinction) that allow us to generate precise expectations of the timings of their origins. Thus, in these instances, the fossil record can provide a powerful test of molecular clock methodology, and why it goes astray, and we have every reason to think these problems are general. [Cambrian explosion; mammalian evolution; molecular clocks.].

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.

Wolbachia infection status and molecular diversity in the species of tribe Tagiadini Mabille, 1878 (Lepidoptera: Hesperiidae) collected in China

Wolbachia, one of the most ubiquitous heritable symbionts in lepidopteran insects, can cause mitochondrial introgression in related host species. We recently found mito-nuclear discordance in the Lepidopteran tribe Tagiadini Mabille 1878 from which Wolbachia has not been reported. In this study, we found that 13 of the 46 species of Tagiadini species tested were positive for Wolbachia. Overall, 14% (15/110) of Tagiadini specimens were infected with Wolbachia and nine new STs were found from 15 isolates. A co-phylogenetic comparison, divergence time estimation and Wolbachia recombination analysis revealed that mito-nuclear discordance in Tagiadini species is not mediated by Wolbachia, but Wolbachia acquisition in Tagiadini appears to have occurred mainly through horizontal transmission rather than codivergence.

The angiosperm radiation played a dual role in the diversification of insects and insect pollinators

Interactions with angiosperms have been hypothesised to play a crucial role in driving diversification among insects, with a particular emphasis on pollinator insects. However, support for coevolutionary diversification in insect-plant interactions is weak. Macroevolutionary studies of insect and plant diversities support the hypothesis that angiosperms diversified after a peak in insect diversity in the Early Cretaceous. Here, we used the family-level fossil record of insects as a whole, and insect pollinator families in particular, to estimate diversification rates and the role of angiosperms on insect macroevolutionary history using a Bayesian process-based approach. We found that angiosperms played a dual role that changed through time, mitigating insect extinction in the Cretaceous and promoting insect origination in the Cenozoic, which is also recovered for insect pollinator families only. Although insects pollinated gymnosperms before the angiosperm radiation, a radiation of new pollinator lineages began as angiosperm lineages increased, particularly significant after 50 Ma. We also found that global temperature, increases in insect diversity, and spore plants were strongly correlated with origination and extinction rates, suggesting that multiple drivers influenced insect diversification and arguing for the investigation of different explanatory variables in further studies.

Genomics of the relict species Baronia brevicornis sheds light on its demographic history and genome size evolution across swallowtail butterflies

Relict species, like coelacanth, gingko, tuatara, are the remnants of formerly more ecologically and taxonomically diverse lineages. It raises the questions of why they are currently species-poor, have restrained ecology, and are often vulnerable to extinction. Estimating heterozygosity level and demographic history can guide our understanding of the evolutionary history and conservation status of relict species. However, few studies have focused on relict invertebrates compared to vertebrates. We sequenced the genome of Baronia brevicornis (Lepidoptera: Papilionidae), which is an endangered species, the sister species of all swallowtail butterflies, and is the oldest lineage of all extant butterflies. From a dried specimen, we were able to generate both long-read and short-read data and assembled a genome of 406 Mb for Baronia. We found a fairly high level of heterozygosity (0.58%) compared to other swallowtail butterflies, which contrasts with its endangered and relict status. Taking into account the high ratio of recombination over mutation, demographic analyses indicated a sharp decline of the effective population size initiated in the last million years. Moreover, the Baronia genome was used to study genome size variation in Papilionidae. Genome sizes are mostly explained by transposable elements activities, suggesting that large genomes appear to be a derived feature in swallowtail butterflies as transposable elements activity is recent and involves different transposable elements classes among species. This first Baronia genome provides a resource for assisting conservation in a flagship and relict insect species as well as for understanding swallowtail genome evolution.

Unraveling phylogenetic relationships and species boundaries in the arid adapted Gerbillus rodents (Muridae: Gerbillinae) by RAD-seq data

Gerbillus is one of the most speciose genera among rodents, with ca. 51 recognized species. Previous attempts to reconstruct the evolutionary history of Gerbillus mainly relied on the mitochondrial cyt-b marker as a source of phylogenetic information. In this study, we utilize RAD-seq genomic data from 37 specimens representing 11 species to reconstruct the phylogenetic tree for Gerbillus, applying concatenation and coalescence methods. We identified four highly supported clades corresponding to the traditionally recognized subgenera: Dipodillus, Gerbillus, Hendecapleura and Monodia. Only two uncertain branches were detected in the resulting trees, with one leading to diversification of the main lineages in the genus, recognized by quartet sampling analysis as uncertain due to possible introgression. We also examined species boundaries for four pairs of sister taxa, including potentially new species from Morocco, using SNAPP. The results strongly supported a speciation model in which all taxa are treated as separate species. The dating analyses confirmed the Plio-Pleistocene diversification of the genus, with the uncertain branch coinciding with the beginning of aridification of the Sahara at the the Plio-Pleistocene boundary. This study aligns well with the earlier analyses based on the cyt-b marker, reaffirming its suitability as an adequate marker for estimating genetic diversity in Gerbillus.

Niche breadth explains the range size of European-centred butterflies, but dispersal ability does not

Aim

The breadth of ecological niches and dispersal abilities have long been discussed as important determinants of species' range sizes. However, studies directly comparing the relative effects of both factors are rare, taxonomically biased and revealed inconsistent results.

Location

Europe.

Time Period

Cenozoic.

Major Taxa

Butterflies, Lepidoptera.

Methods

We relate climate, diet and habitat niche breadth and two indicators of dispersal ability, wingspan and a dispersal tendency index, to the global range size of 369 European-centred butterfly species. The relative effects of these five predictors and their variation across the butterfly phylogeny were assessed by means of phylogenetic generalized least squares models and phylogenetically weighted regressions respectively.

Results

Climate niche breadth was the most important single predictor, followed by habitat and diet niche breadth, while dispersal tendency and wingspan showed no relation to species' range size. All predictors together explained 59% of the variation in butterfly range size. However, the effects of each predictor varied considerably across families and genera.

Main Conclusions

Range sizes of European-centred butterflies are strongly correlated with ecological niche breadth but apparently independent of dispersal ability. The magnitude of range size-niche breadth relationships is not stationary across the phylogeny and is often negatively correlated across the different dimensions of the ecological niche. This variation limits the generalizability of range size-trait relationships across broad taxonomic groups.

Lepidopteran Synteny Units reveal deep chromosomal conservation in butterflies and moths

DNA is compacted into individual particles or chromosomes that form the basic units of inheritance. However, different animals and plants have widely different numbers of chromosomes. This means that we cannot readily tell which chromosomes are related to which. Here, we describe a simple technique that looks at the similarity of genes on each chromosome and thus gives us a true picture of their homology or similarity through evolutionary time. We use this new system to look at the chromosomes of butterflies and moths or Lepidoptera. We term the associated synteny units, Lepidopteran Synteny Units (LSUs). Using a sample of butterfly and moth genomes from across evolutionary time, we show that LSUs form a simple and reliable method of tracing chromosomal homology back through time. Surprisingly, this technique reveals that butterfly and moth chromosomes show conserved blocks dating back to their sister group the Trichoptera. As Lepidoptera have holocentric chromosomes, it will be interesting to see if similar levels of synteny are shown in groups of animals with monocentric chromosomes. The ability to define homology via LSU analysis makes it considerably easier to approach many questions in chromosomal evolution.

Warning Coloration, Body Size, and the Evolution of Gregarious Behavior in Butterfly Larvae

AbstractMany species gain antipredator benefits by combining gregarious behavior with warning coloration, yet there is debate over which trait evolves first and which is the secondary adaptive enhancement. Body size can also influence how predators receive aposematic signals and potentially constrain the evolution of gregarious behavior. To our knowledge, the causative links between the evolution of gregariousness, aposematism, and larger body sizes have not been fully resolved. Here, using the most recently resolved butterfly phylogeny and an extensive new dataset of larval traits, we reveal the evolutionary interactions between important traits linked to larval gregariousness. We show that larval gregariousness has arisen many times across butterflies, and aposematism is a likely prerequisite for gregariousness to evolve. We also find that body size may be an important factor for determining the coloration of solitary, but not gregarious, larvae. Additionally, by exposing artificial larvae to wild avian predation, we show that undefended, cryptic larvae are heavily predated when aggregated but benefit from solitariness, whereas the reverse is true for aposematic prey. Our data reinforce the importance of aposematism for gregarious larval survival while identifying new questions about the roles of body size and toxicity in the evolution of grouping behavior.

Pattern variation is linked to anti-predator coloration in butterfly larvae

Prey animals typically try to avoid being detected and/or advertise to would-be predators that they should be avoided. Both anti-predator strategies primarily rely on colour to succeed, but the specific patterning used is also important. While the role of patterning in camouflage is relatively clear, the design features of aposematic patterns are less well understood. Here, we use a comparative approach to investigate how pattern use varies across a phylogeny of 268 species of cryptic and aposematic butterfly larvae, which also vary in social behaviour. We find that longitudinal stripes are used more frequently by cryptic larvae, and that patterns putatively linked to crypsis are more likely to be used by solitary larvae. By contrast, aposematic larvae are more likely to use horizontal bands and spots, but we find no differences in the use of individual pattern elements between solitary and gregarious aposematic species. However, solitary aposematic larvae are more likely to display multiple pattern elements, whereas those with no pattern are more likely to be gregarious. Our study advances our understanding of how pattern variation, coloration and social behaviour covary across lepidopteran larvae, and highlights new questions about how patterning affects larval detectability and predator responses to aposematic prey.

Molecular phylogeny, historical biogeography, and classification of Pseudocoladenia butterflies (Lepidoptera: Hesperiidae)

The range of the butterfly genus Pseudocoladenia includes several biodiversity hotspots, such as the Himalayas, mountains of Southwest China, and Sundaland. However, the status of some of its species/subspecies remain controversial, and no previous phylogenetic or biogeographic analyses have been conducted. Herein, we determined the systematic relationships and biogeographic history of this genus by reconstructing its phylogeny based on six genes and 76 specimens as representatives of all known species/subspecies. Two species delimitation methods (Bayes Poisson Tree Processes and Bayesian Phylogenetics and Phylogeography) were also employed to assess the status of each taxon. Based on these results and morphological evidence, we identified 12 species and three subspecies in the genus and subsequently classified these into three species groups: P. fatih, P. dea, and P. dan. Five taxa, P. sadakoe (Sonan & Mitono, 1936) stat. nov., P. celebica (Fruhstorfer, 1909) stat. nov., P. fulvescens (Elwes & Edwarda, 1897) stat. nov., P. eacus (Latreille, 1823) stat. nov., and P. fabia (Evans, 1949) stat. nov. were all recognized as independent species. Additionally, two taxa: P. eacus sumatrana (Fruhstorfer, 1909) comb. nov. and P. eacus dhyana (Fruhstorfer, 1909) comb. nov., were placed under P. eacus (Latreille, 1823) stat. nov. as subspecies. Another new species distributed in N. Yunnan, Pseudocoladenia yunnana Fan, Cao & Hou sp. nov., was also discovered and described. Divergence time and ancestral range estimation indicated that the most recent common ancestor of Pseudocoladenia was distributed in the Himalayas-Hengduan Mountain region and Indochina and diverged approximately 14.00 Ma. Continuous and episodic dispersal, vicariance, and extinction were used to determine the current geographic distribution of the genus. The P. fatih group had a prominently disjunct distribution between the Himalaya-Hengduan Mountain and Taiwan. Meanwhile, the P. dan group was first derived in Indochina and subsequently dispersed into the Southeast Asian archipelagoes. This study provides a reference for the evolutionary route of transoceanic distributed species in Asia and elaborates on the causes of biodiversity.

A global phylogeny of butterflies reveals their evolutionary history, ancestral hosts and biogeographic origins

Butterflies are a diverse and charismatic insect group that are thought to have evolved with plants and dispersed throughout the world in response to key geological events. However, these hypotheses have not been extensively tested because a comprehensive phylogenetic framework and datasets for butterfly larval hosts and global distributions are lacking. We sequenced 391 genes from nearly 2,300 butterfly species, sampled from 90 countries and 28 specimen collections, to reconstruct a new phylogenomic tree of butterflies representing 92% of all genera. Our phylogeny has strong support for nearly all nodes and demonstrates that at least 36 butterfly tribes require reclassification. Divergence time analyses imply an origin ~100 million years ago for butterflies and indicate that all but one family were present before the K/Pg extinction event. We aggregated larval host datasets and global distribution records and found that butterflies are likely to have first fed on Fabaceae and originated in what is now the Americas. Soon after the Cretaceous Thermal Maximum, butterflies crossed Beringia and diversified in the Palaeotropics. Our results also reveal that most butterfly species are specialists that feed on only one larval host plant family. However, generalist butterflies that consume two or more plant families usually feed on closely related plants.

Rapid radiation of ant parasitic butterflies during the Miocene aridification of Africa

Africa has undergone a progressive aridification during the last 20 My that presumably impacted organisms and fostered the evolution of life history adaptations. We test the hypothesis that shift to living in ant nests and feeding on ant brood by larvae of phyto-predaceous Lepidochrysops butterflies was an adaptive response to the aridification of Africa that facilitated the subsequent radiation of butterflies in this genus. Using anchored hybrid enrichment we constructed a time-calibrated phylogeny for Lepidochrysops and its closest, non-parasitic relatives in the Euchrysops section (Poloyommatini). We estimated ancestral areas across the phylogeny with process-based biogeographical models and diversification rates relying on time-variable and clade-heterogeneous birth-death models. The Euchrysops section originated with the emerging Miombo woodlands about 22 million years ago (Mya) and spread to drier biomes as they became available in the late Miocene. The diversification of the non-parasitic lineages decreased as aridification intensified around 10 Mya, culminating in diversity decline. In contrast, the diversification of the phyto-predaceous Lepidochrysops lineage proceeded rapidly from about 6.5 Mya when this unusual life history likely first evolved. The Miombo woodlands were the cradle for diversification of the Euchrysops section, and our findings are consistent with the hypothesis that aridification during the Miocene selected for a phyto-predaceous life history in species of Lepidochrysops, with ant nests likely providing caterpillars a safe refuge from fire and a source of food when vegetation was scarce.

Mass production of unvouchered records fails to represent global biodiversity patterns

The ever-increasing human footprint even in very remote places on Earth has inspired efforts to document biodiversity vigorously in case organisms go extinct. However, the data commonly gathered come from either primary voucher specimens in a natural history collection or from direct field observations that are not traceable to tangible material in a museum or herbarium. Although both datasets are crucial for assessing how anthropogenic drivers affect biodiversity, they have widespread coverage gaps and biases that may render them inefficient in representing patterns of biodiversity. Using a large global dataset of around 1.9 billion occurrence records of terrestrial plants, butterflies, amphibians, birds, reptiles and mammals, we quantify coverage and biases of expected biodiversity patterns by voucher and observation records. We show that the mass production of observation records does not lead to higher coverage of expected biodiversity patterns but is disproportionately biased toward certain regions, clades, functional traits and time periods. Such coverage patterns are driven by the ease of accessibility to air and ground transportation, level of security and extent of human modification at each sampling site. Conversely, voucher records are vastly infrequent in occurrence data but in the few places where they are sampled, showed relative congruence with expected biodiversity patterns for all dimensions. The differences in coverage and bias by voucher and observation records have important implications on the utility of these records for research in ecology, evolution and conservation research.

Towards an integrative approach to understanding collective behaviour in caterpillars

To evolve, and remain adaptive, collective behaviours must have a positive impact on overall individual fitness. However, these adaptive benefits may not be immediately apparent owing to an array of interactions with other ecological traits, which can depend on a lineage's evolutionary past and the mechanisms controlling group behaviour. A coherent understanding of how these behaviours evolve, are exhibited, and are coordinated across individuals, therefore requires an integrative approach spanning traditional disciplines in behavioural biology. Here, we argue that lepidopteran larvae are well placed to serve as study systems for investigating the integrative biology of collective behaviour. Lepidopteran larvae display a striking diversity in social behaviour, which illustrates critical interactions between ecological, morphological and behavioural traits. While previous, often classic, work has provided an understanding of how and why collective behaviours evolve in Lepidoptera, much less is known about the developmental and mechanistic basis of these traits. Recent advances in the quantification of behaviour, and the availability of genomic resources and manipulative tools, allied with the exploitation of the behavioural diversity of tractable lepidopteran clades, will change this. In doing so, we will be able to address previously intractable questions that can reveal the interplay between levels of biological variation. This article is part of a discussion meeting issue 'Collective behaviour through time'.

The genetic basis of wing spots in Pieris canidia butterflies

Spots in pierid butterflies and eyespots in nymphalid butterflies are likely non-homologous wing colour pattern elements, yet they share a few features in common. Both develop black scales that depend on the function of the gene spalt, and both might have central signalling cells. This suggests that both pattern elements may be sharing common genetic circuitry. Hundreds of genes have already been associated with the development of nymphalid butterfly eyespot patterns, but the genetic basis of the simpler spot patterns on the wings of pierid butterflies has not been investigated. To facilitate studies of pierid wing patterns, we report a high-quality draft genome assembly for Pieris canidia, the Indian cabbage white. We then conducted transcriptomic analyses of pupal wing tissues sampled from the spot and non-spot regions of P. canidia at 3-6 h post-pupation. A total of 1352 genes were differentially regulated between wing tissues with and without the black spot, including spalt, Krüppel-like factor 10, genes from the Toll, Notch, TGF-β, and FGFR signalling pathways, and several genes involved in the melanin biosynthetic pathway. We identified 14 genes that are up-regulated in both pierid spots and nymphalid eyespots and propose that spots and eyespots share regulatory modules despite their likely independent origins.

Evidence for a single, ancient origin of a genus-wide alternative life history strategy

Understanding the evolutionary origins and factors maintaining alternative life history strategies (ALHS) within species is a major goal of evolutionary research. While alternative alleles causing discrete ALHS are expected to purge or fix over time, one-third of the ~90 species of Colias butterflies are polymorphic for a female-limited ALHS called Alba. Whether Alba arose once, evolved in parallel, or has been exchanged among taxa is currently unknown. Using comparative genome-wide association study (GWAS) and population genomic analyses, we placed the genetic basis of Alba in time-calibrated phylogenomic framework, revealing that Alba evolved once near the base of the genus and has been subsequently maintained via introgression and balancing selection. CRISPR-Cas9 mutagenesis was then used to verify a putative cis-regulatory region of Alba, which we identified using phylogenetic foot printing. We hypothesize that this cis-regulatory region acts as a modular enhancer for the induction of the Alba ALHS, which has likely facilitated its long evolutionary persistence.

Melanism and color saturation of butterfly assemblages: A comparison between a tropical rainforest and a xeric white forest

The study of butterfly coloration has helped to identify the ecological pressures involved in the evolution of animal coloration. However, almost all studies that addressed this issue have focused on species that inhabit more temperate environments, leaving the species and ecological factors of tropical regions mostly understudied. Here, our purpose was to evaluate whether butterfly assemblages from two distinct Neotropical biomes (i.e., tropical rainforest and xeric white forest) differ regarding their melanism and/or color saturation. Our hypotheses were that (1) tropical rainforest butterflies should be more melanic and color saturated, and that (2) butterflies from more open/arid tropical environments should be more melanic on their dorsal wing surfaces than on their ventral wings. Therefore, we quantified melanism and color saturation from dorsal and ventral surfaces of 121 different butterfly species. Comparisons show that rainforest butterflies, when contrasted to white forest butterflies, have more melanic dorsal wing surfaces, which might be seen as a form of protection against parasites. Our data also show that rainforest butterflies, but not white forest species, have darker dorsal wing surfaces, when compared to their own ventral surfaces, a trend that was also found for species inhabiting both biomes, which might be associated to thermoregulatory advantages. At last, our results also point that butterflies' dorsal wing sides present a higher variance between species (regardless of Biome), when compared to their own ventral wing side, an indication that some ecological factor (e.g., predation avoidance) might be exerting a strong homogenizing force on ventral wing coloration.

The Chromosome-Level Genome of Hestina assimilis (Lepidoptera: Nymphalidae) Reveals the Evolution of Saprophagy-Related Genes in Brush-Footed Butterflies

Most butterflies feed on nectar, while some saprophagous butterflies forage on various non-nectar foods. To date, little is known about the genomic and molecular shifts associated with the evolution of the saprophagous feeding strategy. Here, we assembled the high-quality chromosome-level genome of Hestina assimilis to explore its saprophagous molecular and genetic mechanisms. This chromosome-level genome of H. assimilis is 412.82 Mb, with a scaffold N50 of 15.70 Mb. In total, 98.11% of contigs were anchored to 30 chromosomes. Compared with H. assimilis and other Nymphalidae butterflies, the genes of metabolism and detoxification experienced expansions. We annotated 80 cytochrome P450 (CYP) genes in the H. assimilis genome, among which genes belonging to the CYP4 subfamily were significantly expanded (p < 0.01). These P450 genes were unevenly distributed and mainly concentrated on chromosomes 6-9. We identified 33 olfactory receptor (OR), 20 odorant-binding protein (OBP), and six gustatory receptor (GR) genes in the H. assimilis genome, which were fewer than in the nectarivorous Danaus plexippus. A decreased number of OBP, OR, and GR genes implied that H. assimilis should resort less to olfaction and gustation than their nectarivorous counterparts, which need highly specialized olfactory and gustatory functions. Moreover, we found one site under positive selection occurred in residue 996 (phenylalanine) of GR genes exclusive to H. assimilis, which is conservative in most lineages. Our study provides support for the adaptive evolution of feeding habits in butterflies.

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.

Assessing the relative performance of fast molecular dating methods for phylogenomic data

Advances in genome sequencing techniques produced a significant growth of phylogenomic datasets. This massive amount of data represents a computational challenge for molecular dating with Bayesian approaches. Rapid molecular dating methods have been proposed over the last few decades to overcome these issues. However, a comparative evaluation of their relative performance on empirical data sets is lacking. We analyzed 23 empirical phylogenomic datasets to investigate the performance of two commonly employed fast dating methodologies: penalized likelihood (PL), implemented in treePL, and the relative rate framework (RRF), implemented in RelTime. They were compared to Bayesian analyses using the closest possible substitution models and calibration settings. We found that RRF was computationally faster and generally provided node age estimates statistically equivalent to Bayesian divergence times. PL time estimates consistently exhibited low levels of uncertainty. Overall, to approximate Bayesian approaches, RelTime is an efficient method with significantly lower computational demand, being more than 100 times faster than treePL. Thus, to alleviate the computational burden of Bayesian divergence time inference in the era of massive genomic data, molecular dating can be facilitated using the RRF, allowing evolutionary hypotheses to be tested more quickly and efficiently.

Considerations on the Systematics of Neotropical Pierina, with the Description of Two New Species of Phulia Herrich-Schäffer from the Peruvian Andes (Lepidoptera: Pieridae, Pierinae, Pierini)

A comparative analysis of high-Andean Pierina was carried out, including a total of 25 species. Based on morphological evidence, with an emphasis on venation and genitalia and molecular data, using three genetic markers, we confirm the recent subjective synonymy of the generic names Tatochila Butler, 1870, Piercolias, Staudinger, 1894, Hypsochila Ureta, 1955, Infraphulia Field, 1958, Pierphulia Field, 1958, and Theochila Field, 1958 with Phulia Herrich-Schäffer, 1867. Two new species are described, namely Phulia stoddardi Pyrcz & Cerdeña n. sp., from the Andes of Central Peru, which occurs at an unusually high altitude of close to 5000 m a.s.l. in dry puna habitat, and Phulia phantasma Lamas, Willmott & Boyer n. sp., from dry montane forests in northern Peru and southern Ecuador. An overview of high-elevation butterflies is presented, with some discussion on adaptations to this environment.

Early diversifications of angiosperms and their insect pollinators: were they unlinked?

The present-day ubiquity of angiosperm-insect pollination has led to the hypothesis that these two groups coevolved early in their evolutionary history. However, recent fossil discoveries and fossil-calibrated molecular dating analyses challenge the notion that early diversifications of angiosperms and insects were inextricably linked. In this article, we examine (i) the discrepancies between dates of emergence for angiosperms and major clades of insects; (ii) the long history of gymnosperm-insect pollination modes, which likely shaped early angiosperm-insect pollination mutualisms; and (iii) how the K-Pg (Cretaceous-Paleogene) mass extinction event was vital in propelling modern angiosperm-insect mutualisms. We posit that the early diversifications of angiosperms and their insect pollinators were largely decoupled until the end of the Cretaceous.

A Revised Phylogenetic Classification for Viola (Violaceae)

The genus Viola (Violaceae) is among the 40-50 largest genera among angiosperms, yet its taxonomy has not been revised for nearly a century. In the most recent revision, by Wilhelm Becker in 1925, the then-known 400 species were distributed among 14 sections and numerous unranked groups. Here, we provide an updated, comprehensive classification of the genus, based on data from phylogeny, morphology, chromosome counts, and ploidy, and based on modern principles of monophyly. The revision is presented as an annotated global checklist of accepted species of Viola, an updated multigene phylogenetic network and an ITS phylogeny with denser taxon sampling, a brief summary of the taxonomic changes from Becker's classification and their justification, a morphological binary key to the accepted subgenera, sections and subsections, and an account of each infrageneric subdivision with justifications for delimitation and rank including a description, a list of apomorphies, molecular phylogenies where possible or relevant, a distribution map, and a list of included species. We distribute the 664 species accepted by us into 2 subgenera, 31 sections, and 20 subsections. We erect one new subgenus of Viola (subg. Neoandinium, a replacement name for the illegitimate subg. Andinium), six new sections (sect. Abyssinium, sect. Himalayum, sect. Melvio, sect. Nematocaulon, sect. Spathulidium, sect. Xanthidium), and seven new subsections (subsect. Australasiaticae, subsect. Bulbosae, subsect. Clausenianae, subsect. Cleistogamae, subsect. Dispares, subsect. Formosanae, subsect. Pseudorupestres). Evolution within the genus is discussed in light of biogeography, the fossil record, morphology, and particular traits. Viola is among very few temperate and widespread genera that originated in South America. The biggest identified knowledge gaps for Viola concern the South American taxa, for which basic knowledge from phylogeny, chromosome counts, and fossil data is virtually absent. Viola has also never been subject to comprehensive anatomical study. Studies into seed anatomy and morphology are required to understand the fossil record of the genus.

The evolution and diversification of oakleaf butterflies

Oakleaf butterflies in the genus Kallima have a polymorphic wing phenotype, enabling these insects to masquerade as dead leaves. This iconic example of protective resemblance provides an interesting evolutionary paradigm that can be employed to study biodiversity. We integrated multi-omic data analyses and functional validation to infer the evolutionary history of Kallima species and investigate the genetic basis of their variable leaf wing patterns. We find that Kallima butterflies diversified in the eastern Himalayas and dispersed to East and Southeast Asia. Moreover, we find that leaf wing polymorphism is controlled by the wing patterning gene cortex, which has been maintained in Kallima by long-term balancing selection. Our results provide macroevolutionary and microevolutionary insights into a model species originating from a mountain ecosystem.

Stepwise evolution of a butterfly supergene via duplication and inversion

Supergenes maintain adaptive clusters of alleles in the face of genetic mixing. Although usually attributed to inversions, supergenes can be complex, and reconstructing the precise processes that led to recombination suppression and their timing is challenging. We investigated the origin of the BC supergene, which controls variation in warning coloration in the African monarch butterfly, Danaus chrysippus. By generating chromosome-scale assemblies for all three alleles, we identified multiple structural differences. Most strikingly, we find that a region of more than 1 million bp underwent several segmental duplications at least 7.5 Ma. The resulting duplicated fragments appear to have triggered four inversions in surrounding parts of the chromosome, resulting in stepwise growth of the region of suppressed recombination. Phylogenies for the inversions are incongruent with the species tree and suggest that structural polymorphisms have persisted for at least 4.1 Myr. In addition to the role of duplications in triggering inversions, our results suggest a previously undescribed mechanism of recombination suppression through independent losses of divergent duplicated tracts. Overall, our findings add support for a stepwise model of supergene evolution involving a variety of structural changes.

This article is part of the theme issue ‘Genomic architecture of supergenes: causes and evolutionary consequences’.

Genomics-based higher classification of the species-rich Hairstreaks (Lepidoptera: Lycaenidae: Eumaeini)

We propose a higher classification of the lycaenid hairstreak tribe Eumaeini - one of the youngest and most species-rich butterfly tribes - based on autosome, Lepidopteran Z sex chromosome, and mitochondrial protein-coding genes. The subtribe Neolycaenina Korb is a synonym of Callophryidina Tutt, and subtribe Tmolusina Bálint is a synonym of Strephonotina K. Johnson, Austin, Le Crom, & Salazar. Proposed names are Rhammina Prieto & Busby, new subtribe; Timaetina Busby & Prieto, new subtribe; Atlidina Martins & Duarte, new subtribe; Evenina Faynel & Grishin, new subtribe; Jantheclina Robbins & Faynel, new subtribe; Paiwarriina Lamas & Robbins, new subtribe; Cupatheclina Lamas & Grishin, new subtribe; Parrhasiina Busby & Robbins, new subtribe; Ipideclina Martins & Grishin, new subtribe; and Trichonidina Duarte & Faynel, new subtribe. Phylogenetic results from the autosome and Z sex chromosome analyses are similar. Future analyses of datasets with hundreds of terminal taxa may be more practical time-wise by focussing on the smaller number of sex chromosome sequences (2.6% of nuclear protein-coding sequences). The phylogenetic classification and biological summaries for each subtribe suggest that a variety of factors affected Eumaeini diversification. About a dozen kinds of male secondary sexual organs with frequent evolutionary gains and losses occur in Atlidina, Evenina, and Jantheclina (141 species combined). Females have been shown to use these organs to discriminate between conspecific and non-conspecific males, facilitating sympatry among close relatives. Eumaeina, Rhammina, and Timaetina (140 species combined) are overwhelmingly montane with some evidence for a higher incidence of sympatric diversification. Seven Neotropical lineages in five subtribes invaded the temperate parts of the Nearctic Region with a diversification increase in the Callophryidina (262 species). North American Satyrium and Callophrys then invaded the Palearctic at least once each, with a major species-richness increase in Satyrium. The evolution of litter feeding detritivores within Calycopidina (172 species) resulted in an increase in diversification rate compared with its flower-feeding sister lineage. Atlidina, Strephonotina, Parrhasiina, and Strymonina (562 species combined) each contain a mixture of genera that specialize on one or two caterpillar food plant families and genera that are polyphagous. These would be appropriate subtribes to assess how the breadth of caterpillar food plants and the frequency of host shifts affected diversification.

Phylogeny and Biogeographic History of Parnassius Butterflies (Papilionidae: Parnassiinae) Reveal Their Origin and Deep Diversification in West China

Simple Summary

Butterflies of the genus Parnassius are distributed in the mountains across the Northern Hemisphere. Studies have shown that this genus originated in the regions of West China–Central Asia. To further explore the spatiotemporal pattern and driving mechanism of Parnassius diversification, we reconstructed the phylogeny and biogeographic history of Parnassius based on 45 species. Ancestral area reconstruction obtained by using the statistical dispersal–extinction cladogenesis model indicated that Parnassius originated in West China (Qinghai–Tibet Plateau and Xinjiang) during the Middle Miocene. Paleoenvironment changes and host plants were probably influenced by the dispersal of Parnassius butterflies from West China to East Asia, Europe, and North America. Furthermore, ancient gene introgression might have contributed to the spread of Parnassius butterflies from the high mountains of the Qinghai–Tibet Plateau to the low-altitude areas of Central East China. This study will provide an understanding of the phylogeny and biogeographic history of the genus Parnassius.

Abstract

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.

Macrostructural Evolution of the Mitogenome of Butterflies (Lepidoptera, Papilionoidea)

Simple Summary

Papilionoidea is a superfamily of Lepidoptera encompassing about 19,000 species. In the present work, we study the evolution of the structure of the mitogenome of these lepidopterans. The mechanisms generating the eight arrangements known for Papilionoidea were investigated analysing the movements of different mitochondrial genes. Five newly sequenced/assembled mitogenomes were included in our analysis involving more than 600 genomes. We provide new findings that help to understand the evolution of the gene orders MIQGO, IMQGO, 2S1GO, ES1GO and S1NGO in different butterflies. We demonstrate that the evolution of the 2S1GO in Lycaenidae followed a complicated pathway with multiple events of duplication and loss of trnS1 and changes in anticodon. We describe two new gene orders 2FFGO and 4QGO for Ampittia subvittatus (Hesperiidae) and Bhutanitis thaidina (Papilionidae).

Abstract

The mitogenome of the species belonging to the Papilionodea (Lepidoptera) is a double stranded circular molecule containing the 37 genes shared by Metazoa. Eight mitochondrial gene orders are known in the Papilionoidea. MIQGO is the plesiomorphic gene order for this superfamily, while other mitochondrial arrangements have a very limited distribution. 2S1GO gene order is an exception and is present in several Lycaenidae and one species of Hesperiidae. We studied the macrostructural changes generating the gene orders of butterflies by analysing a large data set (611 taxa) containing 5 new mitochondrial sequences/assemblies and 87 de novo annotated mitogenomes. Our analysis supports a possible origin of the intergenic spacer trnQ-nad2, characterising MIQGO, from trnM. We showed that the homoplasious gene order IMQGO, shared by butterflies, species of ants, beetles and aphids, evolved through different transformational pathways. We identify a complicated evolutionary scenario for 2S1GO in Lycaenidae, characterised by multiple events of duplication/loss and change in anticodon of trnS1. We show that the gene orders ES1GO and S1NGO originated through a tandem duplication random loss mechanism. We describe two novel gene orders. Ampittia subvittatus (Hesperiidae) exhibits the gene order 2FFGO, characterised by two copies of trnF, one located in the canonical position and a second placed in the opposite strand between trnR and trnN. Bhutanitis thaidina (Papilionidae) exhibits the gene order 4QGO, characterised by the quadruplication of trnQ.

Differentiation of the Chestnut Tiger Butterfly Parantica sita (Lepidoptera: Nymphalidae: Danainae) in China

The chestnut tiger butterfly, Parantica sita (Kollar) (Lepidoptera: Nymphalidae: Danainae), occurs in Asia, along the Himalayas, and into the Malayan region. Previous studies found three types of mitogenomes with substantial genetic divergence in samples from China. To clarify the level of differentiation within P. sita, we investigated both molecular data and morphological features in 429 individuals from China. Upon examination, mitochondrial cytochrome oxidase subunit I (COI) sequences showed three substantially diverged haplotype groups. Based on microsatellite genotypes, the samples divided into three clusters that were consistent with the COI haplotype groups. With that genetic data, we named three distinguishable P. sita lineages: PS-A, PS-B, and PS-C. We also found obvious morphological differences in wing color, male sex brand, and genitalia structures among the three lineages. According to the published structure of male genitalia, that of PS-A is identical to that of P. s. sita, and that of PS-B is identical to that of P. pedonga. Based on all the results, we tentatively propose dividing P. sita into three species: PS-A (the former P. s. sita) is the typical Parantica sita [Kollar, (1844)], mainly distributed in southwestern China; PS-C (the former P. s. niphonica) is elevated to full species as Parantica niphonica (Moore, 1883), distributed in Taiwan Island and Japan; and PS-B will be Parantica pedongaFujioka, 1970, mainly distributed in Tibet and western Sichuan. Divergence time estimates showed that PS-A separated from the PS-B + PS-C clade about 8.79 million years ago (Ma), when the Hengduan Mountains underwent an appreciable elevation increase, isolating the Tibet population from the others. PS-B and PS-C diverged about 4.87 Ma, in accord with the formation of Taiwan Island mountains. The founder effect may explain why PS-C’s genetic diversity is lower than that of the other clades.

The Angiosperm Terrestrial Revolution and the origins of modern biodiversity

Biodiversity today has the unusual property that 85% of plant and animal species live on land rather than in the sea, and half of these live in tropical rainforests. An explosive boost to terrestrial diversity occurred from c. 100-50 million years ago, the Late Cretaceous and early Palaeogene. During this interval, the Earth-life system on land was reset, and the biosphere expanded to a new level of productivity, enhancing the capacity and species diversity of terrestrial environments. This boost in terrestrial biodiversity coincided with innovations in flowering plant biology and evolutionary ecology, including their flowers and efficiencies in reproduction; coevolution with animals, especially pollinators and herbivores; photosynthetic capacities; adaptability; and ability to modify habitats. The rise of angiosperms triggered a macroecological revolution on land and drove modern biodiversity in a secular, prolonged shift to new, high levels, a series of processes we name here the Angiosperm Terrestrial Revolution.

A supertree of Northern European macromoths

Ecological and life-history data on the Northern European macromoth (Lepidoptera: Macroheterocera) fauna is widely available and ideal for use in answering phylogeny-based research questions: for example, in comparative biology. However, phylogenetic information for such studies lags behind. Here, as a synthesis of all currently available phylogenetic information on the group, we produce a supertree of 114 Northern European macromoth genera (in four superfamilies, with Geometroidea considered separately), providing the most complete phylogenetic picture of this fauna available to date. In doing so, we assess those parts of the phylogeny that are well resolved and those that are uncertain. Furthermore, we identify those genera for which phylogenetic information is currently too poor to include in such a supertree, or entirely absent, as targets for future work. As an aid to studies involving these genera, we provide information on their likely positions within the macromoth tree. With phylogenies playing an ever more important role in the field, this supertree should be useful in informing future ecological and evolutionary studies.

Limitations of phylogenomic data can drive inferred speciation rate shifts

Biodiversity analyses of phylogenomic timetrees have produced many high-profile examples of shifts in the rate of speciation across the tree of life. Temporally correlated events in ecology, climate, and biogeography are frequently invoked to explain these rate shifts. In a re-examination of 15 genomic timetrees and 25 major published studies of the pattern of speciation through time, we observed an unexpected correlation between the timing of reported rate shifts and the information content of sequence alignments. Here, we show that the paucity of sequence variation and insufficient species sampling in phylogenomic datasets are the likely drivers of many inferred speciation rate shifts, rather than the proposed biological explanations. Therefore, data limitations can produce predictable but spurious signals of rate shifts even when speciation rates may be similar across taxa and time. Our results suggest that the reliable detection of speciation rate shifts requires the acquisition and assembly of long phylogenomic alignments with near-complete species sampling and accurate estimates of species richness for the clades of study.

Conserved ancestral tropical niche but different continental histories explain the latitudinal diversity gradient in brush-footed butterflies

The global increase in species richness toward the tropics across continents and taxonomic groups, referred to as the latitudinal diversity gradient, stimulated the formulation of many hypotheses to explain the underlying mechanisms of this pattern. We evaluate several of these hypotheses to explain spatial diversity patterns in a butterfly family, the Nymphalidae, by assessing the contributions of speciation, extinction, and dispersal, and also the extent to which these processes differ among regions at the same latitude. We generate a time-calibrated phylogeny containing 2,866 nymphalid species (~45% of extant diversity). Neither speciation nor extinction rate variations consistently explain the latitudinal diversity gradient among regions because temporal diversification dynamics differ greatly across longitude. The Neotropical diversity results from low extinction rates, not high speciation rates, and biotic interchanges with other regions are rare. Southeast Asia is also characterized by a low speciation rate but, unlike the Neotropics, is the main source of dispersal events through time. Our results suggest that global climate change throughout the Cenozoic, combined with tropical niche conservatism, played a major role in generating the modern latitudinal diversity gradient of nymphalid butterflies.

The first record of the monospecific genus Rhinopalpa (Lepidoptera: Nymphalidae) from China

Background

The family Nymphalidae is the largest group of butterflies with high species richeness. Rhinopalpapolynice (Cramer, [1779]), a forest species, was discovered in the mid-stream of the Yuanjiang-Red River Valley of Yunnan Province for the first time, which represents the first record of the genus Rhinopalpa in China.

New information

The species R. polynice (Cramer, [1779]) is the first record of the genus Rhinopalpa from China. The specimen was collected in the mid-stream of the Yuanjiang-Red River Valley of Yunnan Province. The female genitalia are described for the first time.

Horizontal transmission and recombination of Wolbachia in the butterfly tribe Aeromachini Tutt, 1906 (Lepidoptera: Hesperiidae)

Wolbachia is arguably one of the most ubiquitous heritable symbionts among insects and understanding its transmission dynamics is crucial for understanding why it is so common. While previous research has studied the transmission pathways of Wolbachia in several insect lineages including Lepidoptera, this study takes advantage of data collected from the lepidopteran tribe Aeromachini in an effort to assess patterns of transmission. Twenty-one of the 46 species of Aeromachini species were infected with Wolbachia. Overall, 25% (31/125) of Aeromachini specimens tested were Wolbachia positive. All Wolbachia strains were species-specific except for the wJho strain which appeared to be shared by three host species with a sympatric distribution based on a cophylogenetic comparison between Wolbachia and the Aeromachini species. Two tests of phylogenetic congruence did not find any evidence for cospeciation between Wolbachia strains and their butterfly hosts. The cophylogenetic comparison, divergence time estimation, and Wolbachia recombination analysis revealed that Wolbachia acquisition in Aeromachini appears to have mainly occurred mainly through horizontal transmission rather than codivergence.

The hidden side of diversity: Effects of imperfect detection on multiple dimensions of biodiversity

Studies on ecological communities often address patterns of species distribution and abundance, but few consider uncertainty in counts of both species and individuals when computing diversity measures.We evaluated the extent to which imperfect detection may influence patterns of taxonomic, functional, and phylogenetic diversity in ecological communities.We estimated the true abundance of fruit-feeding butterflies sampled in canopy and understory strata in a subtropical forest. We compared the diversity values calculated by observed and estimated abundance data through the hidden diversity framework. This framework evaluates the deviation of observed diversity when compared with diversities derived from estimated true abundances and whether such deviation represents a bias or a noise in the observed diversity pattern.The hidden diversity values differed between strata for all diversity measures, except for functional richness. The taxonomic measure was the only one where we observed an inversion of the most diverse stratum when imperfect detection was included. Regarding phylogenetic and functional measures, the strata showed distinct responses to imperfect detection, despite the tendency to overestimate observed diversity. While the understory showed noise for the phylogenetic measure, since the observed pattern was maintained, the canopy had biased diversity for the functional metric. This bias occurred since no significant differences were found between strata for observed diversity, but rather for estimated diversity, with the canopy being more clustered.We demonstrate that ignore imperfect detection may lead to unrealistic estimates of diversity and hence to erroneous interpretations of patterns and processes that structure biological communities. For fruit-feeding butterflies, according to their phylogenetic position or functional traits, the undetected individuals triggered different responses in the relationship of the diversity measures to the environmental factor. This highlights the importance to evaluate and include the uncertainty in species detectability before calculating biodiversity measures to describe communities.

Phylogeny, species delimitation and biogeography of the endemic Palaearctic tribe Tomarini (Lepidoptera: Lycaenidae)

The tribe Tomarini is represented by the sole genus Tomares, comprising about eight species distributed from the western Mediterranean to Central Asia. We carried out a multilocus phylogenetic and a biogeographical analysis to test the taxonomy of the genus by several molecular species delimitation methods and reveal patterns shaping the current distribution of Tomares. The phylogenetic analysis based on four molecular markers recovered the monophyly of the genus and recovered two deep-branching lineages: an African clade and an Asian clade. Species delimitation analyses suggested six or ten putative species depending on the method applied. The haplotype network analysis of the Tomares nogelii clade revealed no phylogeographical and taxonomic structure. We consider the taxon Tomares nesimachus (syn. nov.) a synonym of T. nogelii and reinstate Tomares callimachus dentata stat. rev. for populations from south-eastern Turkey. Tomares originated between the early Oligocene and the early Miocene, most probably in south-west Asia. The split of the most recent common ancestor of Tomares occurred between the middle-late Miocene and middle-late Pliocene, probably as a response to increasing aridification and habitat fragmentation. Differentiation of the Asian clade took place in south-west Asia during the Pliocene and Pleistocene and coincided temporally with the evolution of Tomares host plants of the genus Astragalus (Fabaceae).

Miocene Climate and Habitat Change Drove Diversification in Bicyclus, Africa's Largest Radiation of Satyrine Butterflies

Compared to other regions, the drivers of diversification in Africa are poorly understood. We studied a radiation of insects with over 100 species occurring in a wide range of habitats across the Afrotropics to investigate the fundamental evolutionary processes and geological events that generate and maintain patterns of species richness on the continent. By investigating the evolutionary history of Bicyclus butterflies within a phylogenetic framework, we inferred the group's origin at the Oligo-Miocene boundary from ancestors in the Congolian rainforests of central Africa. Abrupt climatic fluctuations during the Miocene (ca. 19-17 Ma) likely fragmented ancestral populations, resulting in at least eight early-divergent lineages. Only one of these lineages appears to have diversified during the drastic climate and biome changes of the early Miocene, radiating into the largest group of extant species. The other seven lineages diversified in forest ecosystems during the late Miocene and Pleistocene when climatic conditions were more favourable-warmer and wetter. Our results suggest changing Neogene climate, uplift of eastern African orogens, and biotic interactions might have had different effects on the various subclades of Bicyclus, producing one of the most spectacular butterfly radiations in Africa.

Predictability of temporal variation in climate and the evolution of seasonal polyphenism in tropical butterfly communities

Phenotypic plasticity in heterogeneous environments can provide tight environment-phenotype matching. However, the prerequisite is a reliable environmental cue(s) that enables organisms to use current environmental information to induce the development of a phenotype with high fitness in a forthcoming environment. Here, we quantify predictability in the timing of precipitation and temperature change to examine how this is associated with seasonal polyphenism in tropical Mycalesina butterflies. Seasonal precipitation in the tropics typically results in distinct selective environments, the wet and dry seasons, and changes in temperature can be a major environmental cue. We sampled communities of Mycalesina butterflies from two seasonal locations and one aseasonal location. Quantifying environmental predictability using wavelet analysis and Colwell's indices confirmed a strong periodicity of precipitation over a 12-month period at both seasonal locations compared to the aseasonal one. However, temperature seasonality and periodicity differed between the two seasonal locations. We further show that: (a) most females from both seasonal locations synchronize their reproduction with the seasons by breeding in the wet season but arresting reproduction in the dry season. In contrast, all species breed throughout the year in the aseasonal location and (b) species from the seasonal locations, but not those from the aseasonal location, exhibited polyphenism in wing pattern traits (eyespot size). We conclude that seasonal precipitation and its predictability are primary factors shaping the evolution of polyphenism in Mycalesina butterflies, and populations or species secondarily evolve local adaptations for cue use that depend on the local variation in the environment.

Rampant Genome-Wide Admixture across the Heliconius Radiation

How frequent is gene flow between species? The pattern of evolution is typically portrayed as a phylogenetic tree, yet gene flow between good species may be an important mechanism in diversification, spreading adaptive traits and leading to a complex pattern of phylogenetic incongruence. This process has thus far been studied mainly among a few closely related species, or in geographically restricted areas such as islands, but not on the scale of a continental radiation. Using a genomic representation of 40 out of 47 species in the genus, we demonstrate that admixture has played a role throughout the evolution of the charismatic Neotropical butterflies Heliconius. Modeling of phylogenetic networks based on the exome uncovers up to 13 instances of interspecific gene flow. Admixture is detected among the relatives of Heliconius erato, as well as between the ancient lineages leading to modern clades. Interspecific gene flow played a role throughout the evolution of the genus, although the process has been most frequent in the clade of Heliconius melpomene and relatives. We identify Heliconius hecalesia and relatives as putative hybrids, including new evidence for introgression at the loci controlling the mimetic wing patterns. Models accounting for interspecific gene flow yield a more complete picture of the radiation as a network, which will improve our ability to study trait evolution in a realistic comparative framework.