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.

Preliminary population studies of the grassland swallowtail butterfly Euryades corethrus (Lepidoptera, Papilionidae)

Euryades corethrus is a Troidini butterfly (Papilionidae, Papilioninae), endemic to grasslands in southern Brazil, Uruguay, Argentina and Paraguay. Formerly abundant, nowadays it is in the Red list of endangered species for those areas. During its larval stage, it feeds on Aristolochia spp, commonly found in southern grasslands. These native grassland areas are diminishing, being converted to crops and pastures, causing habitat loss for Aristolochia and E. corethrus. This study aimed to assess the genetic diversity, population structure and demographic history of E. corethrus. We sampled eight populations from Rio Grande do Sul, Brazil and based on Cytochrome Oxidase subunit I (COI) molecular marker, our results suggest a low genetic variability between populations, presence of gene flow and, consequently, lack of population structure. A single maternally inherited-genetic marker is insufficient for population-level decisions, but barcoding is a useful tool during early stages of population investigation, bringing out genomic diversity patterns within the target species. Those populations likely faced a bottleneck followed by a rapid expansion during the last glaciation and subsequent stabilization in effective population size. Habitat loss is a threat, which might cause isolation, loss of genetic variability and, ultimately, extinction of E. corethrus if no habitat conservation policy is adopted.

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.

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

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

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.

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

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

A review of the occurrence and diversity of the sphragis in butterflies (Lepidoptera, Papilionoidea)

Males of many butterfly species secrete long-lasting mating plugs to prevent their mates from copulating with other males, thus ensuring their sperm will fertilize all future eggs laid. Certain species have further developed a greatly enlarged, often spectacular, externalized plug, termed a sphragis. This distinctive structure results from complex adaptations in both male and female genitalia and is qualitatively distinct from the amorphous, internal mating plugs of other species. Intermediate conditions between internal plug and external sphragis are rare. The term sphragis has often been misunderstood in recent years, hence we provide a formal definition based on accepted usage throughout most of the last century. Despite it being a highly apparent trait, neither the incidence nor diversity of the sphragis has been systematically documented. We record a sphragis or related structure in 273 butterfly species, representing 72 species of Papilionidae in 13 genera, and 201 species of Nymphalidae in 9 genera. These figures represent respectively, 13% of Papilionidae, 3% of Nymphalidae, and 1% of known butterfly species. A well-formed sphragis evolved independently in at least five butterfly subfamilies, with a rudimentary structure also occurring in an additional subfamily. The sphragis is probably the plesiomorphic condition in groups such as Parnassius (Papilionidae: Parnassiinae) and many Acraeini (Nymphalidae: Heliconiinae). Some butterflies, such as those belonging to the Parnassius simo group, have apparently lost the structure secondarily. The material cost of producing the sphragis is considerable. It is typically offset by production of a smaller spermatophore, thus reducing the amount of male-derived nutrients donated to the female during mating for use in oogenesis and/or somatic maintenance. The sphragis potentially represents one of the clearest examples of mate conflict known. Investigating its biology should yield testable hypotheses to further our understanding of the selective processes at play in an 'arms race' between the sexes. This paper provides an overview, which will inform future study.

How the Aridification of Australia Structured the Biogeography and Influenced the Diversification of a Large Lineage of Australian Cicadas

Over the last 30 million years, Australia's landscape has undergone dramatic cooling and drying due to the establishment of the Antarctic Circumpolar Current and change in global CO$_{2}$ levels. Studies have shown that many Australian organisms went extinct during these major cooling events, while others experienced adaptive radiations and increases in diversification rates as a result of exploiting new niches in the arid zone. Despite the many studies on diversification and biogeography in Australia, few have been continent-wide and none have focused on a group of organisms adapted to feeding on plants. We studied 162 species of cicadas in the Australian Pauropsalta complex, a large generic lineage within the tribe Cicadettini. We asked whether there were changes in the diversification rate of Pauropsalta over time and if so: 1) which clades were associated with the rate change? 2) did timing of rate shifts correspond to known periods of dramatic historical climate change, 3) did increases in diversification rate along select lineages correspond to adaptive radiations with movement into the arid zone? To address these questions, we estimated a molecular phylogeny of the Pauropsalta complex using ${\sim}$5300 bp of nucleotide sequence data distributed among five loci (one mtDNA locus and four nDNA loci). We found that this large group of cicadas did not diversify at a constant rate as they spread through Australia; instead the signature of decreasing diversification rate changed roughly around the time of the expansion of the east Antarctic ice sheets ${\sim}$16 Ma and the glaciation of the northern hemisphere ${\sim}$3 Ma. Unlike other Australian taxa, the Pauropsalta complex did not explosively radiate in response to an early invasion of the arid zone. Instead multiple groups invaded the arid zone and experienced rates of diversification similar to mesic-distributed taxa. We found evidence for relictual groups, located in pre-Mesozoic habitat, that have not diversified and continue to reside on mesic hosts in isolated "habitat islands". Future work should focus on groups of similar ages with similar distribution patterns to determine whether this tempo and pattern of diversification and biogeography is consistent with evidence from other phytophagous insects.

Tectonic blocks and molecular clocks

Evolutionary timescales have mainly used fossils for calibrating molecular clocks, though fossils only really provide minimum clade age constraints. In their place, phylogenetic trees can be calibrated by precisely dated geological events that have shaped biogeography. However, tectonic episodes are protracted, their role in vicariance is rarely justified, the biogeography of living clades and their antecedents may differ, and the impact of such events is contingent on ecology. Biogeographic calibrations are no panacea for the shortcomings of fossil calibrations, but their associated uncertainties can be accommodated. We provide examples of how biogeographic calibrations based on geological data can be established for the fragmentation of the Pangaean supercontinent: (i) for the uplift of the Isthmus of Panama, (ii) the separation of New Zealand from Gondwana, and (iii) for the opening of the Atlantic Ocean. Biogeographic and fossil calibrations are complementary, not competing, approaches to constraining molecular clock analyses, providing alternative constraints on the age of clades that are vital to avoiding circularity in investigating the role of biogeographic mechanisms in shaping modern biodiversity.This article is part of the themed issue 'Dating species divergences using rocks and clocks'.

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.

Changes to the Fossil Record of Insects through Fifteen Years of Discovery

The first and last occurrences of hexapod families in the fossil record are compiled from publications up to end-2009. The major features of these data are compared with those of previous datasets (1993 and 1994). About a third of families (>400) are new to the fossil record since 1994, over half of the earlier, existing families have experienced changes in their known stratigraphic range and only about ten percent have unchanged ranges. Despite these significant additions to knowledge, the broad pattern of described richness through time remains similar, with described richness increasing steadily through geological history and a shift in dominant taxa, from Palaeoptera and Polyneoptera to Paraneoptera and Holometabola, after the Palaeozoic. However, after detrending, described richness is not well correlated with the earlier datasets, indicating significant changes in shorter-term patterns. There is reduced Palaeozoic richness, peaking at a different time, and a less pronounced Permian decline. A pronounced Triassic peak and decline is shown, and the plateau from the mid Early Cretaceous to the end of the period remains, albeit at substantially higher richness compared to earlier datasets. Origination and extinction rates are broadly similar to before, with a broad decline in both through time but episodic peaks, including end-Permian turnover. Origination more consistently exceeds extinction compared to previous datasets and exceptions are mainly in the Palaeozoic. These changes suggest that some inferences about causal mechanisms in insect macroevolution are likely to differ as well.

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.

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.

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.

Nymphalid butterflies diversify following near demise at the Cretaceous/Tertiary boundary

The butterfly family Nymphalidae contains some of the most important non-drosophilid insect model systems for evolutionary and ecological studies, yet the evolutionary history of the group has remained shrouded in mystery. We have inferred a robust phylogenetic hypothesis based on sequences of 10 genes and 235 morphological characters for exemplars of 400 of the 540 valid nymphalid genera representing all major lineages of the family. By dating the branching events, we infer that Nymphalidae originated in the Cretaceous at 90 Ma, but that the ancestors of 10-12 lineages survived the end-Cretaceous catastrophe in the Neotropical and Oriental regions. Patterns of diversification suggest extinction of lineages at the Cretaceous/Tertiary boundary (65 Ma) and subsequent elevated speciation rates in the Tertiary.

Impact of duplicate gene copies on phylogenetic analysis and divergence time estimates in butterflies

Background

The increase in availability of genomic sequences for a wide range of organisms has revealed gene duplication to be a relatively common event. Encounters with duplicate gene copies have consequently become almost inevitable in the context of collecting gene sequences for inferring species trees. Here we examine the effect of incorporating duplicate gene copies evolving at different rates on tree reconstruction and time estimation of recent and deep divergences in butterflies.

Results

Sequences from ultraviolet-sensitive (UVRh), blue-sensitive (BRh), and long-wavelength sensitive (LWRh) opsins,EF-1α and COI were obtained from 27 taxa representing the five major butterfly families (5535 bp total). Both BRh and LWRh are present in multiple copies in some butterfly lineages and the different copies evolve at different rates. Regardless of the phylogenetic reconstruction method used, we found that analyses of combined data sets using either slower or faster evolving copies of duplicate genes resulted in a single topology in agreement with our current understanding of butterfly family relationships based on morphology and molecules. Interestingly, individual analyses of BRh and LWRh sequences also recovered these family-level relationships. Two different relaxed clock methods resulted in similar divergence time estimates at the shallower nodes in the tree, regardless of whether faster or slower evolving copies were used, with larger discrepancies observed at deeper nodes in the phylogeny. The time of divergence between the monarch butterfly Danaus plexippus and the queen D. gilippus (15.3–35.6 Mya) was found to be much older than the time of divergence between monarch co-mimic Limenitis archippus and red-spotted purple L. arthemis (4.7–13.6 Mya), and overlapping with the time of divergence of the co-mimetic passionflower butterflies Heliconius erato and H. melpomene (13.5–26.1 Mya). Our family-level results are congruent with recent estimates found in the literature and indicate an age of 84–113 million years for the divergence of all butterfly families.

Conclusion

These results are consistent with diversification of the butterfly families following the radiation of angiosperms and suggest that some classes of opsin genes may be usefully employed for both phylogenetic reconstruction and divergence time estimation.

Biogeography of butterflies in the Australian monsoon tropics

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

The genetic basis of a plant-insect coevolutionary key innovation

Ehrlich and Raven formally introduced the concept of stepwise coevolution using butterfly and angiosperm interactions in an attempt to account for the impressive biological diversity of these groups. However, many biologists currently envision butterflies evolving 50 to 30 million years (Myr) after the major angiosperm radiation and thus reject coevolutionary origins of butterfly biodiversity. The unresolved central tenet of Ehrlich and Raven's theory is that evolution of plant chemical defenses is followed closely by biochemical adaptation in insect herbivores, and that newly evolved detoxification mechanisms result in adaptive radiation of herbivore lineages. Using one of their original butterfly-host plant systems, the Pieridae, we identify a pierid glucosinolate detoxification mechanism, nitrile-specifier protein (NSP), as a key innovation. Larval NSP activity matches the distribution of glucosinolate in their host plants. Moreover, by using five different temporal estimates, NSP seems to have evolved shortly after the evolution of the host plant group (Brassicales) ( approximately 10 Myr). An adaptive radiation of these glucosinolate-feeding Pierinae followed, resulting in significantly elevated species numbers compared with related clades. Mechanistic understanding in its proper historical context documents more ancient and dynamic plant-insect interactions than previously envisioned. Moreover, these mechanistic insights provide the tools for detailed molecular studies of coevolution from both the plant and insect perspectives.

Out-of-Africa origin and dispersal-mediated diversification of the butterfly genus Junonia (Nymphalidae: Nymphalinae)

The relative importance of dispersal and vicariance in the diversification of taxa has been much debated. Within butterflies, a few studies published so far have demonstrated vicariant patterns at the global level. We studied the historical biogeography of the genus Junonia (Nymphalidae: Nymphalinae) at the intercontinental level based on a molecular phylogeny. The genus is distributed over all major biogeographical regions of the world except the Palaearctic. We found dispersal to be the dominant process in the diversification of the genus. The genus originated and started diversifying in Africa about 20 Ma and soon after dispersed into Asia possibly through the Arabian Peninsula. From Asia, there were dispersals into Africa and Australasia, all around 5 Ma. The origin of the New World species is ambiguous; the ancestral may have dispersed from Asia via the Beringian Strait or from Africa over the Atlantic, about 3 Ma. We found no evidence for vicariance at the intercontinental scale. We argue that dispersal is as important as vicariance, if not more, in the global diversification of butterflies.

Evolution of larval host plant associations and adaptive radiation in pierid butterflies

Butterflies in the family Pieridae (Lepidoptera: Papilionoidea) feed as larvae on plants belonging primarily to three distantly related angiosperm orders: Fabales (legumes and allied plants), Brassicales (crucifers and related plants containing mustard oil glucosides), and Santalales ('mistletoes'). However, some utilize plants from 13 other families in a further eight orders. We investigated the evolutionary history of host plant use of the Pieridae in the context of a recent phylogenetic hypothesis of the family, using simple character optimization. Although there is a close association between host plant and butterfly higher classification, we find no evidence for cospeciation but a pattern of repeated colonization and specialization. The ancestral host of the family appears to be Fabaceae or Fabales, with multiple independent shifts to other orders, including three to Santalales. The shift to Brassicales, which contain secondary compounds (glucosinolates), promoted diversification and adaptive radiation within the subfamily Pierinae. Subsequent shifts from crucifers to mistletoes (aerial-stem hemiparasites) facilitated further diversification, and more recent shifts from mistletoes to mistletoe host trees led to exploitation of novel host plants outside the conventional three orders. Possible mechanisms underlying these host shifts are briefly discussed. In the Pierinae, a striking association between host plant, larval and adult behaviour, adult phenotype, and mimicry calls for further research into possible relationships between host specialization, plant chemistry and butterfly palatability.

That Awkward Age for Butterflies: Insights from the Age of the Butterfly Subfamily Nymphalinae (Lepidoptera: Nymphalidae)

The study of the historical biogeography of butterflies has been hampered by a lack of well-resolved phylogenies and a good estimate of the temporal span over which butterflies have evolved. Recently there has been surge of phylogenetic hypotheses for various butterfly groups, but estimating ages of divergence is still in its infancy for this group of insects. The main problem has been the sparse fossil record for butterflies. In this study I have used a surprisingly good fossil record for the subfamily Nymphalinae (Lepidoptera: Nymphalidae) to estimate the ages of diversification of major lineages using Bayesian relaxed clock methods. I have investigated the effects of varying priors on posterior estimates in the analyses. For this data set, it is clear that the prior of the rate of molecular evolution at the ingroup node had the largest effect on the results. Taking this into account, I have been able to arrive at a plausible history of lineage splits, which appears to be correlated with known paleogeological events. The subfamily appears to have diversified soon after the K/T event about 65 million years ago. Several splits are coincident with major paleogeological events, such as the connection of the African and Asian continents about 21 million years ago and the presence of a peninsula of land connecting the current Greater Antilles to the South American continent 35 to 33 million years ago. My results suggest that the age of Nymphalidae is older than the 70 million years speculated to be the age of butterflies as a whole.