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 favorable-warmer and wetter. Our results suggest changing Neogene climate, uplift of eastern African orogens, and biotic interactions have had different effects on the various subclades of Bicyclus, producing one of the most spectacular butterfly radiations in Africa. [Afrotropics; biodiversity; biome; biotic interactions; Court Jester; extinction; grasslands; paleoclimates; Red Queen; refugia forests; dependent-diversification; speciation.].

A release from developmental bias accelerates morphological diversification in butterfly eyespots.

Development can bias the independent evolution of traits sharing ontogenetic pathways, making certain evolutionary changes less likely. The eyespots commonly found on butterfly wings each have concentric rings of differing colors, and these serially repeated pattern elements have been a focus for evo-devo research. In the butterfly family Nymphalidae, eyespots have been shown to function in startling or deflecting predators and to be involved in sexual selection. Previous work on a model species of Mycalesina butterfly, Bicyclus anynana, has provided insights into the developmental control of the size and color composition of individual eyespots. Experimental evolution has also shown that the relative size of a pair of eyespots on the same wing surface is highly flexible, whereas they are resistant to diverging in color composition, presumably due to the underlying shared developmental process. This fixed color composition has been considered as a prime example of developmental bias with significant consequences for wing pattern evolution. Here, we test this proposal by surveying eyespots across the whole subtribe of Mycalesina butterflies and demonstrate that developmental bias shapes evolutionary diversification except in the genus Heteropsis which has gained independent control of eyespot color composition. Experimental manipulations of pupal wings reveal that the bias has been released through a novel regional response of the wing tissue to a conserved patterning signal. Our study demonstrates that development can bias the evolutionary independence of traits, but it also shows how bias can be released through developmental innovations, thus, allowing rapid morphological change, facilitating evolutionary diversification.

Out of Asia: Intercontinental dispersals after the Eocene-Oligocene transition shaped the zoogeography of Limenitidinae butterflies (Lepidoptera: Nymphalidae).

Most members of the nymphalid subfamily Limenitidinae are distributed in tropical regions of Africa, Asia, and the Americas. Previous studies have inferred their higher-level phylogeny and found that Southeast Asia seems to be the center of origin, with numerous dispersal events to other continents. However, the complete biogeographic history of Limenitidinae butterflies is still largely unknown. We sampled 181 taxa from 164 species and used a metagenomic method to obtain 40 genes (mitogenomes and three nuclear ribosomal loci) for inferring the historical biogeography of the group. We find that Limenitidinae originated in eastern Asia during the early Eocene (ca. 52 Ma) and started to diversify and disperse into Africa before the end of Eocene. Intercontinental exchanges between Africa and eastern Asia continued in the early Miocene: Asian Adoliadini and Asian endemic taxa Bhagadatta had African origins in the Oligocene, whereas African Neptini dispersed in the opposite direction from Asia in the early Miocene. In addition, ancestors of the tribes Limenitidini and Adoliadini dispersed into the Neotropics and Australasia multiple times during the early-to-middle Miocene. Eastern Asia is the center of origin of the tribe Limenitidini, with several taxa disjunctly distributed in eastern Asia and the Americas. Our work provides a robust phylogenetic hypothesis of relationships in the tribe Limenitidini and suggests that the alala-species group of Adelpha should be placed in the genus Limenitis. Renamed taxa comb. nov. based on our findings are listed in the text.

Is Sexual Conflict a Driver of Speciation? A Case Study With a Tribe of Brush-footed Butterflies.

Understanding the evolutionary mechanisms governing the uneven distribution of species richness across the tree of life is a great challenge in biology. Scientists have long argued that sexual conflict is a key driver of speciation. This hypothesis, however, has been highly debated in light of empirical evidence. Recent advances in the study of macroevolution make it possible to test this hypothesis with more data and increased accuracy. In the present study, we use phylogenomics combined with four different diversification rate analytical approaches to test whether sexual conflict is a driver of speciation in brush-footed butterflies of the tribe Acraeini. The presence of a sphragis, an external mating plug found in most species among Acraeini, was used as a proxy for sexual conflict. Diversification analyses statistically rejected the hypothesis that sexual conflict is associated with shifts in diversification rates in Acraeini. This result contrasts with earlier studies and suggests that the underlying mechanisms driving diversification are more complex than previously considered. In the case of butterflies, natural history traits acting in concert with abiotic factors possibly play a stronger role in triggering speciation than does sexual conflict. [Acraeini butterflies; arms race; exon capture phylogenomics; Lepidoptera macroevolution; sexual selection; sphragis.].

Quantifying the effects of species traits on predation risk in nature: A comparative study of butterfly wing damage.

Evading predators is a fundamental aspect of the ecology and evolution of all prey animals. In studying the influence of prey traits on predation risk, previous researchers have shown that crypsis reduces attack rates on resting prey, predation risk increases with increased prey activity, and rapid locomotion reduces attack rates and increases chances of surviving predator attacks. However, evidence for these conclusions is nearly always based on observations of selected species under artificial conditions. In nature, it remains unclear how defensive traits such as crypsis, activity levels and speed influence realized predation risk across species in a community. Whereas direct observations of predator-prey interactions in nature are rare, insight can be gained by quantifying bodily damage caused by failed predator attacks. We quantified how butterfly species traits affect predation risk in nature by determining how defensive traits correlate with wing damage caused by failed predation attempts, thereby providing the first robust multi-species comparative analysis of predator-induced bodily damage in wild animals. For 34 species of fruit-feeding butterflies in an African forest, we recorded wing damage and quantified crypsis, activity levels and flight speed. We then tested for correlations between damage parameters and species traits using comparative methods that account for measurement error. We detected considerable differences in the extent, location and symmetry of wing surface loss among species, with smaller differences between sexes. We found that males (but not females) of species that flew faster had substantially less wing surface loss. However, we found no correlation between cryptic coloration and symmetrical wing surface loss across species. In species in which males appeared to be more active than females, males had a lower proportion of symmetrical wing surface loss than females. Our results provide evidence that activity greatly influences the probability of attacks and that flying rapidly is effective for escaping pursuing predators in the wild, but we did not find evidence that cryptic species are less likely to be attacked while at rest.

Evolution of Hypolimnas butterflies (Nymphalidae): Out-of-Africa origin and Wolbachia-mediated introgression.

Hypolimnas butterflies (Nymphalidae), commonly known as eggflies, are a popular model system for studying a wide range of ecological questions including mimicry, polymorphism, wing pattern evolution, and Wolbachia-host interactions. The lack of a time-calibrated phylogeny for this group has precluded understanding its evolutionary history. We reconstruct a species-level phylogeny using a nine gene dataset and estimate species divergence times. Based on the resulting tree, we investigate the taxon's historical biogeography, examine the evolution of host plant preferences, and test the hypothesis that the endosymbiotic bacterium Wolbachia mediates gene transfer between species. Our analyses indicate that the species are grouped within three strongly supported, deeply divergent clades. However, relationships among these three clades are uncertain. In addition, many Hypolimnas species are not monophyletic or monophyletic with weak support, suggesting widespread incomplete lineage sorting and/or introgression. Biogeographic analysis strongly indicates that the genus diverged from its ancestor in Africa and subsequently dispersed to Asia; the strength of this result is not affected by topological uncertainties. While the larvae of African species feed almost exclusively on Urticaceae, larvae of species found further east often feed on several additional families. Interestingly, we found an identical mitochondrial haplotype in two Hypolimnas species, H. bolina and H. alimena, and a strong association between this mitotype and the Wolbachia strain wBol1a. Future investigations should explore the plausibility of Wolbachia-mediated introgression between species.

Conserved patterns of integrated developmental plasticity in a group of polyphenic tropical butterflies.

BACKGROUND: Developmental plasticity is thought to have profound macro-evolutionary effects, for example, by increasing the probability of establishment in new environments and subsequent divergence into independently evolving lineages. In contrast to plasticity optimized for individual traits, phenotypic integration, which enables a concerted response of plastic traits to environmental variability, may affect the rate of local adaptation by constraining independent responses of traits to selection. Using a comparative framework, this study explores the evolution of reaction norms for a variety of life history and morphological traits across five related species of mycalesine butterflies from the Old World tropics. RESULTS: Our data indicate that an integrated response of a suite of key traits is shared amongst these species. Interestingly, the traits that make up the functional suite are all known to be regulated by ecdysteroid signalling in Bicyclus anynana, one of the species included in this study, suggesting the same underlying hormonal regulator may be conserved within this group of polyphenic butterflies. We also detect developmental thresholds for the expression of alternative morphs. CONCLUSIONS: The phenotypic plasticity of a broad suite of morphological and life history traits is integrated and shared among species from three geographically independent lineages of mycalesine butterflies, despite considerable periods of independent evolution and exposure to disparate environments. At the same time, we have detected examples of evolutionary change where independent traits show different patterns of reaction norms. We argue that the expression of more robust phenotypes may occur by shifting developmental thresholds beyond the boundaries of the typical environmental variation.

Revision of the Bicyclus sciathis species group (Lepidoptera: Nymphalidae) with descriptions of four new species and corrected distributional records.

In this paper we present a thorough revision of the sciathis species group of the butterfly genus Bicyclus (Kirby). Type materials are discussed and in several cases lectotypes are assigned to specimens from original type series. Four new, and morphologically distinct, species are described (B. elishiae Brattström sp.n., B. heathi Brattström sp.n., B. sigiussidorum Brattström sp.n. and B. subtilisurae Brattström sp.n.), along with a comprehensive molecular phylogeny that includes exemplar taxa of all currently recognized species. We also investigate the types of all previously synonymized taxa and in the process invalidate the name B. ewondo Libert. This was done after finding the previously missing holotype of B. makomensis (Strand), which clearly belongs to the same species and thereby gives the older name priority. The phylogeny showed that some distinctly different species were surprisingly closely related, suggesting a high rate of morphological evolution in parts of the sciathis group. The distributional records for the group are updated after investigating over 1700 specimens kept in a range of museum collections. Many species previously thought to be broadly sympatric were found to have much more restricted ranges, with the previous overestimations probably based on misidentified specimens. The higher level of allopatry now established will make identification of many morphologically similar species easier. The fact that species often have smaller ranges than previously known, meaning that the level of endemism for African butterflies is likely to be higher than current estimates, has important implications for conservation management. An identification key for males of all 13 currently recognized species in the species group is included. This published work has been registered in ZooBank, http://zoobank.org/urn:lsid:zoobank.org:pub:837A9D4C-779A-4497-8176-7151D409DFA5.

Systematics and historical biogeography of the old world butterfly subtribe Mycalesina (Lepidoptera: Nymphalidae: Satyrinae).

BACKGROUND: Butterflies of the subtribe Mycalesina have radiated successfully in almost all habitat types in Africa, Madagascar, the Indian subcontinent, Indo-China and Australasia. Studies aimed at understanding the reasons behind the evolutionary success of this spectacular Old World butterfly radiation have been hampered by the lack of a stable phylogeny for the group. Here, we have reconstructed a robust phylogenetic framework for the subtribe using 10 genes from 195 exemplar taxa. RESULTS: We recovered seven well supported clades within the subtribe corresponding to the five traditional genera (Lohora, Heteropsis, Hallelesis, Bicyclus, Mycalesis), one as recently revised (Mydosama) and one newly revised genus (Culapa). The phylogenetic relationships of these mycalesine genera have been robustly established for the first time. Within the proposed phylogenetic framework, we estimated the crown age of the subtribe to be 40 Million years ago (Mya) and inferred its ultimate origin to be in Asia. Our results reveal both vicariance and dispersal as factors responsible for the current widespread distribution of the group in the Old World tropics. We inferred that the African continent has been colonized at least twice by Asian mycalesines within the last 26 and 23 Mya. In one possible scenario, an Asian ancestor gave rise to Heteropsis on continental Africa, which later dispersed into Madagascar and most likely back colonised Asia. The second colonization of Africa by Asian ancestors resulted in Hallelesis and Bicyclus on continental Africa, the descendants of which did not colonise other regions but rather diversified only in continental Africa. The genera Lohora and Mydosama are derivatives of ancestors from continental Asia. CONCLUSION: Our proposed time-calibrated phylogeny now provides a solid framework within which we can implement mechanistic studies aimed at unravelling the ecological and evolutionary processes that culminated in the spectacular radiation of mycalesines in the Old World tropics.

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.

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.

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.

Out-of-Africa again: a phylogenetic hypothesis of the genus Charaxes (Lepidoptera: Nymphalidae) based on five gene regions.

Despite the long popularity of Charaxes among collectors and researchers, their evolutionary history is largely unknown. The current and accepted species groupings and relationships within the genus are based exclusively on adult morphology and life histories. Here, we examine the monophyly and evolutionary affinities of the species-groups within the genus Charaxes and explore how they relate to members of their closest genera (Euxanthe, Polyura and Palla) using 4167bp of sequence data from five (1 mitochondrial and 4 nuclear) gene regions. Within the proposed phylogenetic framework, we estimate ages of divergence within the genus and also reconstruct their historical biogeography. We included representatives of all known species-groups in Africa and Asia, all known species of Euxanthe and Palla and two exemplar species of Polyura. We found the genus Charaxes to be a paraphyletic group with regard to the genera Polyura and Euxanthe, contrary to the earlier assumption of monophyly. We found that 13 out of 16 morphologically defined species-groups with more than one species were strongly supported monophyletic clades. Charaxes nichetes is the sister group to all the other Charaxes. Polyura grouped with the Zoolina and Pleione species-groups as a well-supported clade, and Euxanthe grouped with the Lycurgus species-group. Our results indicated that the common ancestor of Charaxes diverged from the common ancestor of Palla in the mid Eocene (45 million years ago) in (Central) Africa and began diversifying to its extant members 15 million years later. Most of the major diversifications within the genus occurred between the late Oligocene and Miocene when the global climates were putatively undergoing drastic fluctuations. A considerable number of extant species diverged from sister species during the Pliocene. A dispersal-vicariance analysis suggests that many dispersal rather than vicariance events resulted in the distribution of the extant species. The genus Polyura and the Indo-Australian Charaxes are most likely the results of three independent colonizations of Asia by African Charaxes in the Miocene. We synonymize the genera Polyura (syn. nov.) and Euxanthe (syn. nov.) with Charaxes, with the currently circumscribed Charaxes subdivided into five subgenera to reflect its phylogeny.

Relative contribution of neutral and deterministic processes in shaping fruit-feeding butterfly assemblages in Afrotropical forests.

The unified neutral theory of biodiversity and biogeography has gained the status of a quantitative null model for explaining patterns in ecological (meta)communities. The theory assumes that individuals of trophically similar species are functionally equivalent. We empirically evaluate the relative contribution of neutral and deterministic processes in shaping fruit-feeding butterfly assemblages in three tropical forests in Africa, using both direct (confronting the neutral model with species abundance data) and indirect approaches (testing the predictions of neutral theory using data other than species abundance distributions). Abundance data were obtained by sampling butterflies using banana baited traps set at the forest canopy and understorey strata. Our results indicate a clear consistency in the kind of species or species groups observed at either the canopy or understorey in the three studied communities. Furthermore, we found significant correlation between some flight-related morphological traits and species abundance at the forest canopy, but not at the understorey. Neutral theory's contribution to explaining our data lies largely in identifying dispersal limitation as a key process regulating fruit-feeding butterfly community structure. Our study illustrates that using species abundance data alone in evaluating neutral theory can be informative, but is insufficient. Species-level information such as habitat preference, host plants, geographical distribution, and phylogeny is essential in elucidating the processes that regulate biodiversity community structures and patterns.

A Comprehensive and Dated Phylogenomic Analysis of Butterflies.

Butterflies (Papilionoidea), with over 18,000 described species [1], have captivated naturalists and scientists for centuries. They play a central role in the study of speciation, community ecology, biogeography, climate change, and plant-insect interactions and include many model organisms and pest species [2, 3]. However, a robust higher-level phylogenetic framework is lacking. To fill this gap, we inferred a dated phylogeny by analyzing the first phylogenomic dataset, including 352 loci (> 150,000 bp) from 207 species representing 98% of tribes, a 35-fold increase in gene sampling and 3-fold increase in taxon sampling over previous studies [4]. Most data were generated with a new anchored hybrid enrichment (AHE) [5] gene kit (BUTTERFLY1.0) that includes both new and frequently used (e.g., [6]) informative loci, enabling direct comparison and future dataset merging with previous studies. Butterflies originated around 119 million years ago (mya) in the late Cretaceous, but most extant lineages diverged after the Cretaceous-Paleogene (K-Pg) mass-extinction 65 mya. Our analyses support swallowtails (Papilionidae) as sister to all other butterflies, followed by skippers (Hesperiidae) + the nocturnal butterflies (Hedylidae) as sister to the remainder, indicating a secondary reversal from diurnality to nocturnality. The whites (Pieridae) were strongly supported as sister to brush-footed butterflies (Nymphalidae) and blues + metalmarks (Lycaenidae and Riodinidae). Ant association independently evolved once in Lycaenidae and twice in Riodinidae. This study overturns prior notions of the taxon's evolutionary history, as many long-recognized subfamilies and tribes are para- or polyphyletic. It also provides a much-needed backbone for a revised classification of butterflies and for future comparative studies including genome evolution and ecology.

Revision of the Bicyclus ignobilis species-group (Lepidoptera: Nymphalidae: Satyrinae) with descriptions of two new species.

The ignobilis-group of the genus Bicyclus Kirby 1871 is revised. The species-group contains six species with a distinct wing pattern, but limited intraspecific variation, distributed across tropical African rainforest. We investigate a set of more than 1000 specimens from a range of museum collections, including some type material, and thoroughly update the biogeographical knowledge for the group. We also describe two new species as members of the group. The included species are: Bicyclus ignobilis (Butler 1870) stat. rev., B. rileyi Condamin 1961, B. maesseni Condamin 1971, B. brakefieldi Brattström 2012, B. ottossoni sp. nov. and B. vandeweghei sp. nov. Due to observing a gradual morphological cline within B. ignobilis without any sharp transitions we suppress the previously identified subspecies B. ignobilis eurini Condamin & Fox 1963 syn. nov. and B. ignobilis acutus Condamin 1965 syn. nov.