Phylogenetic and Comparative Genomics Study of Papilionidae Based on Mitochondrial Genomes

Most species of Papilionidae are large and beautiful ornamental butterflies. They are recognized as model organisms in ecology, evolutionary biology, genetics, and conservation biology but present numerous unresolved phylogenetic problems. Complete mitochondrial genomes (mitogenomes) have been widely used in phylogenetic studies of butterflies, but mitogenome knowledge within the family Papilionidae is limited, and its phylogeny is far from resolved. In this study, we first report the mitogenome of Byasa confusa from the subfamily Papilioninae of Papilionidae. The mitogenome of B. confusa is 15,135 bp in length and contains 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and an AT-rich control region (CR), closely mirroring the genomic structure observed in related butterfly species. Comparative analysis of 77 Papilionidae mitogenomes shows gene composition and order to be identical to that of an ancestral insect, and the AT bias, Ka/Ks, and relative synonymous codon usage (RSCU) are all consistent with that of other reported butterfly mitogenomes. We conducted phylogenetic analyses using maximum-likelihood (ML) and Bayesian-inference (BI) methods, with 77 Papilionidae species as ingroups and two species of Nymphalidae and Lycaenidae as outgroups. The phylogenetic analysis indicated that B. confusa were clustered within Byasa. The phylogenetic trees show the monophyly of the subfamily Papilioninae and the tribes Leptocircini, Papilionini, and Troidini. The data supported the following relationships in tribe level on Papilioninae: (((Troidini + Papilionini) + Teinopalpini) + Leptocircini). The divergence time analysis suggests that Papilionidae originated in the late Creataceous. Overall, utilizing the largest number of Papilionidae mitogenomes sequenced to date, with the current first exploration in a phylogenetic analysis on Papilionidae (including four subfamilies), this study comprehensively reveals the mitogenome characteristics and mitogenome-based phylogeny, providing information for further studies on the mitogenome, phylogeny, evolution, and taxonomic revision of the Papilionidae family.

Male and female contributions to diversity among birdwing butterfly images

Machine learning (ML) newly enables tests for higher inter-species diversity in visible phenotype (disparity) among males versus females, predictions made from Darwinian sexual selection versus Wallacean natural selection, respectively. Here, we use ML to quantify variation across a sample of > 16,000 dorsal and ventral photographs of the sexually dimorphic birdwing butterflies (Lepidoptera: Papilionidae). Validation of image embedding distances, learnt by a triplet-trained, deep convolutional neural network, shows ML can be used for automated reconstruction of phenotypic evolution achieving measures of phylogenetic congruence to genetic species trees within a range sampled among genetic trees themselves. Quantification of sexual disparity difference (male versus female embedding distance), shows sexually and phylogenetically variable inter-species disparity. Ornithoptera exemplify high embedded male image disparity, diversification of selective optima in fitted multi-peak OU models and accelerated divergence, with cases of extreme divergence in allopatry and sympatry. However, genus Troides shows inverted patterns, including comparatively static male embedded phenotype, and higher female than male disparity - though within an inferred selective regime common to these females. Birdwing shapes and colour patterns that are most phenotypically distinctive in ML similarity are generally those of males. However, either sex can contribute majoritively to observed phenotypic diversity among species.

Complete mitochondrial genome data and phylogenetic analysis of Papilio macilentus Janson, 1877 (Lepidoptera: Papilionoidea: Papilionidae)

In the present study, the complete mitochondrial genome (mitogenome) of the Papilio macilentus (Lepidoptera: Papilionoidea: Papilionidae) was sequenced by next-generation sequencing method. The mitochondrial genome is a circular DNA molecule of 15,264 bp in size with 80.7% AT content, including 37 genes (13 protein-coding genes, 2 rRNA genes, and 22 tRNA genes), and a long non-coding region (Control region). All protein-coding genes are initiated by ATN codons, and terminated with TAA, TAG, or single T. All tRNAs can be folded into common clover leaf secondary structure, except trn-S1. Phylogenetic analyses based on 13 protein-coding genes and 2 rRNA genes using maximum likelihood and Bayesian inference confirmed that P. macilentus and Papilio memnon are clustered into a clade, and revealed the relationships between Papilionini, Troidini, Teinopaippini and Leptocircini.

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.

Air temperature drives the evolution of mid-infrared optical properties of butterfly wings

This study uncovers a correlation between the mid-infrared emissivity of butterfly wings and the average air temperature of their habitats across the world. Butterflies from cooler climates have a lower mid-infrared emissivity, which limits heat losses to surroundings, and butterflies from warmer climates have a higher mid-infrared emissivity, which enhances radiative cooling. The mid-infrared emissivity showed no correlation with other investigated climatic factors. Phylogenetic independent contrasts analysis indicates the microstructures of butterfly wings may have evolved in part to regulate mid-infrared emissivity as an adaptation to climate, rather than as phylogenetic inertia. Our findings offer new insights into the role of microstructures in thermoregulation and suggest both evolutionary and physical constraints to butterflies' abilities to adapt to climate change.

Genome-wide macroevolutionary signatures of key innovations in butterflies colonizing new host plants

The mega-diversity of herbivorous insects is attributed to their co-evolutionary associations with plants. Despite abundant studies on insect-plant interactions, we do not know whether host-plant shifts have impacted both genomic adaptation and species diversification over geological times. We show that the antagonistic insect-plant interaction between swallowtail butterflies and the highly toxic birthworts began 55 million years ago in Beringia, followed by several major ancient host-plant shifts. This evolutionary framework provides a valuable opportunity for repeated tests of genomic signatures of macroevolutionary changes and estimation of diversification rates across their phylogeny. We find that host-plant shifts in butterflies are associated with both genome-wide adaptive molecular evolution (more genes under positive selection) and repeated bursts of speciation rates, contributing to an increase in global diversification through time. Our study links ecological changes, genome-wide adaptations and macroevolutionary consequences, lending support to the importance of ecological interactions as evolutionary drivers over long time periods.

Tick mitochondrial genomes: structural characteristics and phylogenetic implications

Ticks are obligate blood-sucking arachnid ectoparasites from the order Acarina, and many are notorious as vectors of a wide variety of zoonotic pathogens. However, the systematics of ticks in several genera is still controversial. The mitochondrial genome (mt-genome) has been widely used in arthropod phylogeny, molecular evolution and population genetics. With the development of sequencing technologies, an increasing number of tick mt-genomes have been sequenced and annotated. To date, 63 complete tick mt-genomes are available in the NCBI database, and these genomes have become an increasingly important genetic resource and source of molecular markers in phylogenetic studies of ticks in recent years. The present review summarizes all available complete mt-genomes of ticks in the NCBI database and analyses their characteristics, including structure, base composition and gene arrangement. Furthermore, a phylogenetic tree was constructed using mitochondrial protein-coding genes (PCGs) and ribosomal RNA (rRNA) genes from ticks. The results will provide important clues for deciphering new tick mt-genomes and establish a foundation for subsequent taxonomic research.

Sequence capture across large phylogenetic scales by using pooled PCR-generated baits: A case study of Lepidoptera

Sequence capture across large phylogenetic scales is not easy because hybridization capture is only effective when the genetic distance between the bait and target is small. Here, we propose a simple but effective strategy to tackle this issue: pooling DNA from a number of selected representative species of different clades to prepare PCR-generated baits to minimize the genetic distance between the bait and target. To demonstrate the utility of this strategy, we newly developed a set of universal nuclear markers (including 94 nuclear protein-coding genes) for Lepidoptera, a superdiverse insect group. We used a DNA pool from six lepidopteran species (representing six superfamilies) to prepare PCR baits for the 94 markers. These homemade PCR baits were used to capture sequence data from 43 species of 17 lepidopteran families, and 94% of the target loci were recovered. We constructed two data sets from the obtained data (one containing ~90 kb target coding sequences and the other containing ~120 kb target + flanking coding sequences). Both data sets yielded highly similar and well-resolved trees with 90% of nodes having >95% bootstrap support. Our capture experiment indicated that using DNA mixtures pooled from different clade-representative species of Lepidoptera to prepare PCR baits can reliably capture a large number of targeted nuclear markers across different Lepidoptera lineages. We hope that this newly developed nuclear marker set will serve as a new phylogenetic tool for Lepidoptera phylogenetics, and the PCR bait preparation strategy can facilitate the application of sequence capture techniques by researchers to accelerate data collection.

Whole Genome Shotgun Phylogenomics Resolves the Pattern and Timing of Swallowtail Butterfly Evolution

Evolutionary relationships have remained unresolved in many well-studied groups, even though advances in next-generation sequencing and analysis, using approaches such as transcriptomics, anchored hybrid enrichment, or ultraconserved elements, have brought systematics to the brink of whole genome phylogenomics. Recently, it has become possible to sequence the entire genomes of numerous nonbiological models in parallel at reasonable cost, particularly with shotgun sequencing. Here, we identify orthologous coding sequences from whole-genome shotgun sequences, which we then use to investigate the relevance and power of phylogenomic relationship inference and time-calibrated tree estimation. We study an iconic group of butterflies—swallowtails of the family Papilionidae—that has remained phylogenetically unresolved, with continued debate about the timing of their diversification. Low-coverage whole genomes were obtained using Illumina shotgun sequencing for all genera. Genome assembly coupled to BLAST-based orthology searches allowed extraction of 6621 orthologous protein-coding genes for 45 Papilionidae species and 16 outgroup species (with 32% missing data after cleaning phases). Supermatrix phylogenomic analyses were performed with both maximum-likelihood (IQ-TREE) and Bayesian mixture models (PhyloBayes) for amino acid sequences, which produced a fully resolved phylogeny providing new insights into controversial relationships. Species tree reconstruction from gene trees was performed with ASTRAL and SuperTriplets and recovered the same phylogeny. We estimated gene site concordant factors to complement traditional node-support measures, which strengthens the robustness of inferred phylogenies. Bayesian estimates of divergence times based on a reduced data set (760 orthologs and 12% missing data) indicate a mid-Cretaceous origin of Papilionoidea around 99.2 Ma (95% credibility interval: 68.6–142.7 Ma) and Papilionidae around 71.4 Ma (49.8–103.6 Ma), with subsequent diversification of modern lineages well after the Cretaceous-Paleogene event. These results show that shotgun sequencing of whole genomes, even when highly fragmented, represents a powerful approach to phylogenomics and molecular dating in a group that has previously been refractory to resolution.

Climate is a strong predictor of near-infrared reflectance but a poor predictor of colour in butterflies

Colour variation across climatic gradients is a common ecogeographical pattern; yet there is long-standing contention over underlying causes, particularly selection for thermal benefits. We tested the evolutionary association between climate gradients and reflectance of near-infrared (NIR) wavelengths, which influence heat gain but are not visible to animals. We measured ultraviolet (UVA), visible (Vis) and NIR reflectance from calibrated images of 372 butterfly specimens from 60 populations (49 species, five families) spanning the Australian continent. Consistent with selection for thermal benefits, the association between climate and reflectance was stronger for NIR than UVA-Vis wavelengths. Furthermore, climate predicted reflectance of the thorax and basal wing, which are critical to thermoregulation; but it did not predict reflectance of the entire wing, which has a variable role in thermoregulation depending on basking behaviour. These results provide evidence that selection for thermal benefits has shaped the reflectance properties of butterflies.

Phylogenetics of moth-like butterflies (Papilionoidea: Hedylidae) based on a new 13-locus target capture probe set

The Neotropical moth-like butterflies (Hedylidae) are perhaps the most unusual butterfly family. In addition to being species-poor, this family is predominantly nocturnal and has anti-bat ultrasound hearing organs. Evolutionary relationships among the 36 described species are largely unexplored. A new, target capture, anchored hybrid enrichment probe set ('BUTTERFLY2.0') was developed to infer relationships of hedylids and some of their butterfly relatives. The probe set includes 13 genes that have historically been used in butterfly phylogenetics. Our dataset comprised of up to 10,898 aligned base pairs from 22 hedylid species and 19 outgroups. Eleven of the thirteen loci were successfully captured from all samples, and the remaining loci were captured from ≥94% of samples. The inferred phylogeny was consistent with recent molecular studies by placing Hedylidae sister to Hesperiidae, and the tree had robust support for 80% of nodes. Our results are also consistent with morphological studies, with Macrosoma tipulata as the sister species to all remaining hedylids, followed by M. semiermis sister to the remaining species in the genus. We tested the hypothesis that nocturnality evolved once from diurnality in Hedylidae, and demonstrate that the ancestral condition was likely diurnal, with a shift to nocturnality early in the diversification of this family. The BUTTERFLY2.0 probe set includes standard butterfly phylogenetics markers, captures sequences from decades-old museum specimens, and is a cost-effective technique to infer phylogenetic relationships of the butterfly tree of life.

Jamaica's Critically Endangered Butterfly: A Review of the Biology and Conservation Status of the Homerus Swallowtail (Papilio (Pterourus) homerus Fabricius)

The Homerus swallowtail, Papilio (Pterourus) homerus Fabricius, is listed as an endangered species and is endemic to the Caribbean island of Jamaica. The largest butterfly in the Western Hemisphere, P. homerus once inhabited seven of Jamaica's 14 parishes and consisted of at least three populations; however, now only two stronghold populations remain, a western population in the rugged Cockpit Country and an eastern population in the Blue and John Crow Mountains. Despite numerous studies of its life history, much about the population biology, including estimates of total numbers of individuals in each population, remains unknown. In addition, a breeding program is needed to establish an experimental population, which could be used to augment wild populations and ensure the continued survival of the species. Here, we present a review of the biology of P. homerus and recommendations for a conservation plan.

Phylogeny and Evolution of Lepidoptera

Until recently, deep-level phylogeny in Lepidoptera, the largest single radiation of plant-feeding insects, was very poorly understood. Over the past two decades, building on a preceding era of morphological cladistic studies, molecular data have yielded robust initial estimates of relationships both within and among the ∼43 superfamilies, with unsolved problems now yielding to much larger data sets from high-throughput sequencing. Here we summarize progress on lepidopteran phylogeny since 1975, emphasizing the superfamily level, and discuss some resulting advances in our understanding of lepidopteran evolution.

Ten genes and two topologies: an exploration of higher relationships in skipper butterflies (Hesperiidae)

Despite multiple attempts to infer the higher-level phylogenetic relationships of skipper butterflies (Family Hesperiidae), uncertainties in the deep clade relationships persist. The most recent phylogenetic analysis included fewer than 30% of known genera and data from three gene markers. Here we reconstruct the higher-level relationships with a rich sampling of ten nuclear and mitochondrial markers (7,726 bp) from 270 genera and find two distinct but equally plausible topologies among subfamilies at the base of the tree. In one set of analyses, the nuclear markers suggest two contrasting topologies, one of which is supported by the mitochondrial dataset. However, another set of analyses suggests mito-nuclear conflict as the reason for topological incongruence. Neither topology is strongly supported, and we conclude that there is insufficient phylogenetic evidence in the molecular dataset to resolve these relationships. Nevertheless, taking morphological characters into consideration, we suggest that one of the topologies is more likely.

Nutrition shapes life-history evolution across species

Nutrition is a key component of life-history theory, yet we know little about how diet quality shapes life-history evolution across species. Here, we test whether quantitative measures of nutrition are linked to life-history evolution across 96 species of butterflies representing over 50 independent diet shifts. We find that butterflies feeding on high nitrogen host plants as larvae are more fecund, but their eggs are smaller relative to their body size. Nitrogen and sodium content of host plants are also both positively related to eye size. Some of these relationships show pronounced lineage-specific effects. Testis size is not related to nutrition. Additionally, the evolutionary timing of diet shifts is not important, suggesting that nutrition affects life histories regardless of the length of time a species has been adapting to its diet. Our results suggest that, at least for some lineages, species with higher nutrient diets can invest in a range of fitness-related traits like fecundity and eye size while allocating less to each egg as offspring have access to a richer diet. These results have important implications for the evolution of life histories in the face of anthropogenic changes in nutrient availability.

Tsukiyamaia, a new genus of the tribe Baorini (Lepidoptera, Hesperiidae, Hesperiinae)

Skippers of the tribe Baorini are evidently a monophyletic group in the subfamily Hesperiinae. In this study, a new Baorini member Tsukiyamaia albimacula gen. n. et sp. n. is described from north Myanmar, southwest China and north Vietnam. Despite its peculiar and striking wing-pattern, this new genus has some important characters of Baorini, such as a broad and bifid uncus and a well-developed gnathos. Based on an analysis of male genitalia and the molecular phylogenies inferred from both mitochondrial and nuclear genes (28 taxa, total aligned length: 2968 bp), it is proposed that the genus Tsukiyamaia is closely related to the genus Polytremis, which has high species diversity in China. This study not only describes a new skipper but also highlights that Tsukiyamaia is important in clarifying phylogenetic relationship of Polytremis and its allies.

Repeated Reticulate Evolution in North American Papilio machaon Group Swallowtail Butterflies

Hybridization between distinct populations or species is increasingly recognized as an important process for generating biodiversity. However, the interaction between hybridization and speciation is complex, and the diverse evolutionary outcomes of hybridization are difficult to differentiate. Here we characterize potential hybridization in a species group of swallowtail butterflies using microsatellites, DNA sequences, and morphology, and assess whether adaptive introgression or homoploid hybrid speciation was the primary process leading to each putative hybrid lineage. Four geographically separated hybrid populations were identified in the Papilio machaon species group. One distinct mitochondrial DNA clade from P. machaon was fixed in three hybrid taxa (P. brevicauda, P. joanae, and P. m. kahli), while one hybrid swarm (P. zelicaon x machaon) exhibited this hybrid mtDNA clade as well as widespread parental mtDNA haplotypes from both parental species. Microsatellite markers and morphology showed variable admixture and intermediacy, ranging from signatures of prolonged differential introgression from the paternal species (P. polyxenes/P. zelicaon) to current gene flow with both parental species. Divergences of the hybrid lineages dated to early- to mid-Pleistocene, suggesting that repeated glaciations and subsequent range shifts of parental species, particularly P. machaon hudsonianus, facilitated initial hybridization. Although each lineage is distinct, P. joanae is the only taxon with sufficient evidence (ecological separation from parental species) to define it as a homoploid hybrid species. The repetition of hybridization in this group provides a valuable foundation for future research on hybridization, and these results emphasize the potential for hybridization to drive speciation in diverse ways.

Phylogeny and Historical Biogeography of Asian Pterourus Butterflies (Lepidoptera: Papilionidae): A Case of Intercontinental Dispersal from North America to East Asia

The phylogenetic status of the well-known Asian butterflies often known as Agehana (a species group, often treated as a genus or a subgenus, within Papilio sensu lato) has long remained unresolved. Only two species are included, and one of them especially, Papilio maraho, is not only rare but near-threatened, being monophagous on its vulnerable hostplant, Sassafras randaiense (Lauraceae). Although the natural history and population conservation of "Agehana" has received much attention, the biogeographic origin of this group still remains enigmatic. To clarify these two questions, a total of 86 species representatives within Papilionidae were sampled, and four genes (concatenated length 3842 bp) were used to reconstruct their phylogenetic relationships and historical scenarios. Surprisingly, "Agehana" fell within the American Papilio subgenus Pterourus and not as previously suggested, phylogenetically close to the Asian Papilio subgenus Chilasa. We therefore formally synonymize Agehana with Pterourus. Dating and biogeographic analysis allow us to infer an intercontinental dispersal of an American ancestor of Asian Pterourus in the early Miocene, which was coincident with historical paleo-land bridge connections, resulting in the present "East Asia-America" disjunction distribution. We emphasize that species exchange between East Asia and America seems to be a quite frequent occurrence in butterflies during the Oligocene to Miocene climatic optima.

Spectrally tuned structural and pigmentary coloration of birdwing butterfly wing scales

The colourful wing patterns of butterflies play an important role for enhancing fitness; for instance, by providing camouflage, for interspecific mate recognition, or for aposematic display. Closely related butterfly species can have dramatically different wing patterns. The phenomenon is assumed to be caused by ecological processes with changing conditions, e.g. in the environment, and also by sexual selection. Here, we investigate the birdwing butterflies, Ornithoptera, the largest butterflies of the world, together forming a small genus in the butterfly family Papilionidae. The wings of these butterflies are marked by strongly coloured patches. The colours are caused by specially structured wing scales, which act as a chirped multilayer reflector, but the scales also contain papiliochrome pigments, which act as a spectral filter. The combined structural and pigmentary effects tune the coloration of the wing scales. The tuned colours are presumably important for mate recognition and signalling. By applying electron microscopy, (micro-)spectrophotometry and scatterometry we found that the various mechanisms of scale coloration of the different birdwing species strongly correlate with the taxonomical distribution of Ornithoptera species.

Deciphering the evolution of birdwing butterflies 150 years after Alfred Russel Wallace

One hundred and fifty years after Alfred Wallace studied the geographical variation and species diversity of butterflies in the Indomalayan-Australasian Archipelago, the processes responsible for their biogeographical pattern remain equivocal. We analysed the macroevolutionary mechanisms accounting for the temporal and geographical diversification of the charismatic birdwing butterflies (Papilionidae), a major focus of Wallace's pioneering work. Bayesian phylogenetics and dating analyses of the birdwings were conducted using mitochondrial and nuclear genes. The combination of maximum likelihood analyses to estimate biogeographical history and diversification rates reveals that diversity-dependence processes drove the radiation of birdwings, and that speciation was often associated with founder-events colonizing new islands, especially in Wallacea. Palaeo-environment diversification models also suggest that high extinction rates occurred during periods of elevated sea level and global warming. We demonstrated a pattern of spatio-temporal habitat dynamics that continuously created or erased habitats suitable for birdwing biodiversity. Since birdwings were extinction-prone during the Miocene (warmer temperatures and elevated sea levels), the cooling period after the mid-Miocene climatic optimum fostered birdwing diversification due to the release of extinction. This also suggests that current global changes may represent a serious conservation threat to this flagship group.

A new Heraclides swallowtail (Lepidoptera, Papilionidae) from North America is recognized by the pattern on its neck

Heraclidesrumiko Shiraiwa & Grishin, sp. n. is described from southwestern United States, Mexico, and Central America (type locality: USA, Texas, Duval County). It is closely allied to Heraclidescresphontes (Cramer, 1777) and the two species are sympatric in central Texas. The new species is diagnosed by male genitalia and exhibits a nearly 3% difference from Heraclidescresphontes in the COI DNA barcode sequence of mitochondrial DNA. The two Heraclides species can usually be told apart by the shape and size of yellow spots on the neck, by the wing shape, and the details of wing patterns. "Western Giant Swallowtail" is proposed as the English name for Heraclidesrumiko. To stabilize nomenclature, neotype for Papiliocresphontes Cramer, 1777, an eastern United States species, is designated from Brooklyn, New York, USA; and lectotype for Papiliothoas Linnaeus, 1771 is designated from Suriname. We sequenced DNA barcodes and ID tags of nearly 400 Papilionini specimens completing coverage of all Heraclides species. Comparative analyses of DNA barcodes, genitalia, and facies suggest that Heraclidesoviedo (Gundlach, 1866), reinstated status, is a species-level taxon rather than a subspecies of Heraclidesthoas (Linnaeus, 1771); and Heraclidespallas (G. Gray, [1853]), reinstated status, with its subspecies HeraclidesPapiliobajaensis (J. Brown & Faulkner, 1992), comb. n., and Heraclidesanchicayaensis Constantino, Le Crom & Salazar, 2002, stat. n., are not conspecific with Heraclidesastyalus (Godart, 1819).

Pigmentary and photonic coloration mechanisms reveal taxonomic relationships of the Cattlehearts (Lepidoptera: Papilionidae: Parides)

BACKGROUND

The colorful wing patterns of butterflies, a prime example of biodiversity, can change dramatically within closely related species. Wing pattern diversity is specifically present among papilionid butterflies. Whether a correlation between color and the evolution of these butterflies exists so far remained unsolved.

RESULTS

We here investigate the Cattlehearts, Parides, a small Neotropical genus of papilionid butterflies with 36 members, the wings of which are marked by distinctly colored patches. By applying various physical techniques, we investigate the coloration toolkit of the wing scales. The wing scales contain two different, wavelength-selective absorbing pigments, causing pigmentary colorations. Scale ridges with multilayered lamellae, lumen multilayers or gyroid photonic crystals in the scale lumen create structural colors that are variously combined with these pigmentary colors.

CONCLUSIONS

The pigmentary and structural traits strongly correlate with the taxonomical distribution of Parides species. The experimental findings add crucial insight into the evolution of butterfly wing scales and show the importance of morphological parameter mapping for butterfly phylogenetics.

Phylogeography of Heliconius cydno and its closest relatives: disentangling their origin and diversification

The origins of the extraordinary diversity within the Neotropics have long fascinated biologists and naturalists. Yet, the underlying factors that have given rise to this diversity remain controversial. To test the relative importance of Quaternary climatic change and Neogene tectonic and paleogeographic reorganizations in the generation of biodiversity, we examine intraspecific variation across the Heliconius cydno radiation and compare this variation to that within the closely related Heliconius melpomene and Heliconius timareta radiations. Our data, which consist of both mtDNA and genome-scan data from nearly 2250 amplified fragment length polymorphism (AFLP) loci, reveal a complex history of differentiation and admixture at different geographic scales. Both mtDNA and AFLP phylogenies suggest that H. timareta and H. cydno are probably geographic extremes of the same radiation that probably diverged from H. melpomene prior to the Pliocene-Pleistocene boundary, consistent with hypotheses of diversification that rely on geological events in the Pliocene. The mtDNA suggests that this radiation originated in Central America or the northwestern region of South America, with a subsequent colonization of the eastern and western slopes of the Andes. Our genome-scan data indicate significant admixture among sympatric H. cydno/H. timareta and H. melpomene populations across the extensive geographic ranges of the two radiations. Within H. cydno, both mtDNA and AFLP data indicate significant population structure at local scales, with strong genetic differences even among adjacent H. cydno colour pattern races. These genetic patterns highlight the importance of past geoclimatic events, intraspecific gene flow, and local population differentiation in the origin and establishment of new adaptive forms.

Towards a mitogenomic phylogeny of Lepidoptera

The backbone phylogeny of Lepidoptera remains unresolved, despite strenuous recent morphological and molecular efforts. Molecular studies have focused on nuclear protein coding genes, sometimes adding a single mitochondrial gene. Recent advances in sequencing technology have, however, made acquisition of entire mitochondrial genomes both practical and economically viable. Prior phylogenetic studies utilised just eight of 43 currently recognised lepidopteran superfamilies. Here, we add 23 full and six partial mitochondrial genomes (comprising 22 superfamilies of which 16 are newly represented) to those publically available for a total of 24 superfamilies and ask whether such a sample can resolve deeper lepidopteran phylogeny. Using recoded datasets we obtain topologies that are highly congruent with prior nuclear and/or morphological studies. Our study shows support for an expanded Obtectomera including Gelechioidea, Thyridoidea, plume moths (Alucitoidea and Pterophoroidea; possibly along with Epermenioidea), Papilionoidea, Pyraloidea, Mimallonoidea and Macroheterocera. Regarding other controversially positioned higher taxa, Doidae is supported within the new concept of Drepanoidea and Mimallonidae sister to (or part of) Macroheterocera, while among Nymphalidae butterflies, Danainae and not Libytheinae are sister to the remainder of the family. At the deepest level, we suggest that a tRNA rearrangement occurred at a node between Adeloidea and Ditrysia+Palaephatidae+Tischeriidae.

Recurrent horizontal transfers of Chapaev transposons in diverse invertebrate and vertebrate animals

Horizontal transfer (HT) of a transposable element (TE) into a new genome is regarded as an important force to drive genome variation and biological innovation. In addition, HT also plays an important role in the persistence of TEs in eukaryotic genomes. Here, we provide the first documented example for the repeated HT of three families of Chapaev transposons in a wide range of animal species, including mammals, reptiles, jawed fishes, lampreys, insects, and in an insect bracovirus. Multiple alignments of the Chapaev transposons identified in these species revealed extremely high levels of nucleotide sequence identity (79-99%), which are inconsistent with vertical evolution given the deep divergence time separating these host species. Rather, the discontinuous distribution amongst species and lack of purifying selection acting on these transposons strongly suggest that they were independently and horizontally transferred into these species lineages. The detection of Chapaev transposons in an insect bracovirus indicated that these viruses might act as a possible vector for the horizontal spread of Chapaev transposons. One of the Chapaev families was also shared by lampreys and some of their common hosts (such as sturgeon and paddlefish), which suggested that parasite-host interaction might facilitate HTs.

Climate-Driven Reshuffling of Species and Genes: Potential Conservation Roles for Species Translocations and Recombinant Hybrid Genotypes

Comprising 50%-75% of the world's fauna, insects are a prominent part of biodiversity in communities and ecosystems globally. Biodiversity across all levels of biological classifications is fundamentally based on genetic diversity. However, the integration of genomics and phylogenetics into conservation management may not be as rapid as climate change. The genetics of hybrid introgression as a source of novel variation for ecological divergence and evolutionary speciation (and resilience) may generate adaptive potential and diversity fast enough to respond to locally-altered environmental conditions. Major plant and herbivore hybrid zones with associated communities deserve conservation consideration. This review addresses functional genetics across multi-trophic-level interactions including "invasive species" in various ecosystems as they may become disrupted in different ways by rapid climate change. "Invasive genes" (into new species and populations) need to be recognized for their positive creative potential and addressed in conservation programs. "Genetic rescue" via hybrid translocations may provide needed adaptive flexibility for rapid adaptation to environmental change. While concerns persist for some conservationists, this review emphasizes the positive aspects of hybrids and hybridization. Specific implications of natural genetic introgression are addressed with a few examples from butterflies, including transgressive phenotypes and climate-driven homoploid recombinant hybrid speciation. Some specific examples illustrate these points using the swallowtail butterflies (Papilionidae) with their long-term historical data base (phylogeographical diversity changes) and recent (3-decade) climate-driven temporal and genetic divergence in recombinant homoploid hybrids and relatively recent hybrid speciation of Papilio appalachiensis in North America. Climate-induced "reshuffling" (recombinations) of species composition, genotypes, and genomes may become increasingly ecologically and evolutionarily predictable, but future conservation management programs are more likely to remain constrained by human behavior than by lack of academic knowledge.

Butterfly morphology in a molecular age -- does it still matter in butterfly systematics?

We review morphological characters considered important for understanding butterfly phylogeny and evolution in the light of recent large-scale molecular phylogenies of the group. A number of the most important morphological works from the past half century are reviewed and morphological character evolution is reassessed based on the most recent phylogenetic results. In particular, higher level butterfly morphology is evaluated based on a very recent study combining an elaborate morphological dataset with a similar molecular one. Special attention is also given to the families Papilionidae, Nymphalidae and Hesperiidae which have all seen morphological and molecular efforts come together in large, combined works in recent years. In all of the examined cases the synergistic effect of combining elaborate morphological datasets with ditto molecular clearly outweigh the merits of either data type analysed on its own (even for 'genome size' molecular datasets). It is evident that morphology, far from being obsolete or arcane, still has an immensely important role to play in butterfly (and insect) phylogenetics. Not least because understanding morphology is essential for understanding and evaluating the evolutionary scenarios phylogenetic trees are supposed to illustrate.

Biogeographic and diversification patterns of Neotropical Troidini butterflies (Papilionidae) support a museum model of diversity dynamics for Amazonia

BACKGROUND

The temporal and geographical diversification of Neotropical insects remains poorly understood because of the complex changes in geological and climatic conditions that occurred during the Cenozoic. To better understand extant patterns in Neotropical biodiversity, we investigated the evolutionary history of three Neotropical swallowtail Troidini genera (Papilionidae). First, DNA-based species delimitation analyses were conducted to assess species boundaries within Neotropical Troidini using an enlarged fragment of the standard barcode gene. Molecularly delineated species were then used to infer a time-calibrated species-level phylogeny based on a three-gene dataset and Bayesian dating analyses. The corresponding chronogram was used to explore their temporal and geographical diversification through distinct likelihood-based methods.

RESULTS

The phylogeny for Neotropical Troidini was well resolved and strongly supported. Molecular dating and biogeographic analyses indicate that the extant lineages of Neotropical Troidini have a late Eocene (33-42 Ma) origin in North America. Two independent lineages (Battus and Euryades+Parides) reached South America via the GAARlandia temporary connection, and later became extinct in North America. They only began substantive diversification during the early Miocene in Amazonia. Macroevolutionary analysis supports the "museum model" of diversification, rather than Pleistocene refugia, as the best explanation for the diversification of these lineages.

CONCLUSIONS

This study demonstrates that: (i) current Neotropical biodiversity may have originated ex situ; (ii) the GAARlandia bridge was important in facilitating invasions of South America; (iii) colonization of Amazonia initiated the crown diversification of these swallowtails; and (iv) Amazonia is not only a species-rich region but also acted as a sanctuary for the dynamics of this diversity. In particular, Amazonia probably allowed the persistence of old lineages and contributed to the steady accumulation of diversity over time with constant net diversification rates, a result that contrasts with previous studies on other South American butterflies.

Phylogenetic codivergence supports coevolution of mimetic Heliconius butterflies

The unpalatable and warning-patterned butterflies Heliconius erato and Heliconius melpomene provide the best studied example of mutualistic Müllerian mimicry, thought-but rarely demonstrated-to promote coevolution. Some of the strongest available evidence for coevolution comes from phylogenetic codivergence, the parallel divergence of ecologically associated lineages. Early evolutionary reconstructions suggested codivergence between mimetic populations of H. erato and H. melpomene, and this was initially hailed as one of the most striking known cases of coevolution. However, subsequent molecular phylogenetic analyses found discrepancies in phylogenetic branching patterns and timing (topological and temporal incongruence) that argued against codivergence. We present the first explicit cophylogenetic test of codivergence between mimetic populations of H. erato and H. melpomene, and re-examine the timing of these radiations. We find statistically significant topological congruence between multilocus coalescent population phylogenies of H. erato and H. melpomene. Cophylogenetic historical reconstructions support repeated codivergence of mimetic populations, from the base of the sampled radiations. Pairwise distance correlation tests, based on our coalescent analyses plus recently published AFLP and wing colour pattern gene data, also suggest that the phylogenies of H. erato and H. melpomene show significant topological congruence. Divergence time estimates, based on a Bayesian coalescent model, suggest that the evolutionary radiations of H. erato and H. melpomene occurred over the same time period, and are compatible with a series of temporally congruent codivergence events. Our results suggest that differences in within-species genetic divergence are the result of a greater overall effective population size for H. erato relative to H. melpomene and do not imply incongruence in the timing of their phylogenetic radiations. Repeated codivergence between Müllerian co-mimics, predicted to exert mutual selection pressures, strongly suggests coevolution. Our results therefore support a history of reciprocal coevolution between Müllerian co-mimics characterised by phylogenetic codivergence and parallel phenotypic change.

Cretaceous origin and repeated tertiary diversification of the redefined butterflies

Although the taxonomy of the ca 18 000 species of butterflies and skippers is well known, the family-level relationships are still debated. Here, we present, to our knowledge, the most comprehensive phylogenetic analysis of the superfamilies Papilionoidea, Hesperioidea and Hedyloidea to date based on morphological and molecular data. We reconstructed their phylogenetic relationships using parsimony and Bayesian approaches. We estimated times and rates of diversification along lineages in order to reconstruct their evolutionary history. Our results suggest that the butterflies, as traditionally understood, are paraphyletic, with Papilionidae being the sister-group to Hesperioidea, Hedyloidea and all other butterflies. Hence, the families in the current three superfamilies should be placed in a single superfamily Papilionoidea. In addition, we find that Hedylidae is sister to Hesperiidae, and this novel relationship is supported by two morphological characters. The families diverged in the Early Cretaceous but diversified after the Cretaceous-Palaeogene event. The diversification of butterflies is characterized by a slow speciation rate in the lineage leading to Baronia brevicornis, a period of stasis by the skippers after divergence and a burst of diversification in the lineages leading to Nymphalidae, Riodinidae and Lycaenidae.

What causes latitudinal gradients in species diversity? Evolutionary processes and ecological constraints on swallowtail biodiversity

The latitudinal diversity gradient (LDG) is one of the most striking ecological patterns on our planet. Determining the evolutionary causes of this pattern remains a challenging task. To address this issue, previous LDG studies have usually relied on correlations between environmental variables and species richness, only considering evolutionary processes indirectly. Instead, we use a phylogenetically integrated approach to investigate the ecological and evolutionary processes responsible for the global LDG observed in swallowtail butterflies (Papilionidae). We find evidence for the 'diversification rate hypothesis' with different diversification rates between two similarly aged tropical and temperate clades. We conclude that the LDG is caused by (1) climatically driven changes in both clades based on evidence of responses to cooling and warming events, and (2) distinct biogeographical histories constrained by tropical niche conservatism and niche evolution. This multidisciplinary approach provides new findings that allow better understanding of the factors that shape LDGs.