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.

DNA barcoding employs short, standardized gene regions (5’ segment of mitochondrial cytochrome oxidase subunit I for animals) as an internal tag to enable species identification. Prior studies have indicated that it performs this task well, because interspecific variation at cytochrome oxidase subunit I is typically much greater than intraspecific variation. However, most previous studies have focused on local faunas only, and critics have suggested two reasons why barcoding should be less effective in species identification when the geographical coverage is expanded. They suggested that many recently diverged taxa will be excluded from local analyses because they are allopatric. Second, intraspecific variation may be seriously underestimated by local studies, because geographical variation in the barcode region is not considered. In this paper, we analyse how adding a geographical dimension affects barcode resolution, examining 353 butterfly species from Central Asia. Despite predictions, we found that geographically separated and recently diverged allopatric species did not show, on average, less sequence differentiation than recently diverged sympatric taxa. Although expanded geographical coverage did substantially increase intraspecific variation reducing the barcoding gap between species, this did not decrease species identification using neighbour‐joining clustering. The inclusion of additional populations increased the number of paraphyletic entities, but did not impede species‐level identification, because paraphyletic species were separated from their monophyletic relatives by substantial sequence divergence. Thus, this study demonstrates that DNA barcoding remains an effective identification tool even when taxa are sampled from a large geographical area.

Uncovering cryptic biodiversity is essential for understanding evolutionary processes and patterns of ecosystem functioning, as well as for nature conservation. As European butterflies are arguably the best-studied group of invertebrates in the world, the discovery of a cryptic species, twenty years ago, within the common wood white Leptidea sinapis was a significant event, and these butterflies have become a model to study speciation. Here we show that the so-called 'sibling' Leptidea actually consist of three species. The new species can be discriminated on the basis of either DNA or karyological data. Such an unexpected discovery challenges our current knowledge on biodiversity, exemplifying how a widespread species can remain unnoticed even within an intensely studied natural model system for speciation.

Butterflies in the large Palearctic genus Agrodiaetus (Lepidoptera: Lycaenidae) are extremely uniform and exhibit few distinguishing morphological characters. However, these insects are distinctive in one respect: as a group they possess among the greatest interspecific karyotype diversity in the animal kingdom, with chromosome numbers (n) ranging from 10 to 125. The monophyly of Agrodiaetus and its systematic position relative to other groups within the section Polyommatus have been controversial. Characters from the mitochondrial genes for cytochrome oxidases I and II and from the nuclear gene for elongation factor 1α were used to reconstruct the phylogeny of Agrodiaetus using maximum parsimony and Bayesian phylogenetic methods. Ninety-one individuals, encompassing most of the taxonomic diversity of Agrodiaetus, and representatives of 14 related genera were included in this analysis. Our data indicate that Agrodiaetus is monophyletic. Representatives of the genus Polyommatus (sensu stricto) are the closest relatives. The sequences of the Agrodiaetus taxa in this analysis are tentatively arranged into 12 clades, only 1 of which corresponds to a species group traditionally recognized in Agrodiaetus. Heterogeneous substitution rates across a recovered topology were homogenized with a nonparametric rate-smoothing algorithm before the application of a molecular clock. Two published estimates of substitution rates dated the origin of Agrodiaetus between 2.51 and 3.85 million years ago. During this time, there was heterogeneity in the rate and direction of karyotype evolution among lineages within the genus. Karyotype instability has evolved independently three times in the section Polyommatus, within the lineages Agrodiaetus, Lysandra, and Plebicula. Rapid karyotype diversification may have played a significant role in the radiation of the genus Agrodiaetus.