Phenotypic plasticity in tropical butterflies is linked to climatic seasonality on a macroevolutionary scale.

Phenotypic plasticity can be adaptive in fluctuating environments by providing rapid environment-phenotype matching and this applies particularly in seasonal environments. African Bicyclus butterflies have repeatedly colonized seasonal savannahs from ancestral forests around the Late Miocene and many species now exhibit seasonal polyphenism. On a macroevolutionary scale, it can be expected that savannah species will exhibit higher plasticity due to them experiencing stronger environmental seasonality than forest species. We quantified seasonality using environmental niche modelling, and surveyed the degree of plasticity in a key wing pattern element (eyespot size) using museum specimens. We show that species occurring in highly seasonal environments display strong plasticity, while species in less seasonal or aseasonal environments exhibit surprisingly variable degrees of plasticity, including strong to no plasticity. Furthermore, eyespot size plasticity has a moderate phylogenetic signal and the ancestral Bicyclus likely exhibited some degree of plasticity. We propose hypotheses to explain the range of plasticity patterns seen in less seasonal environments, and generate testable predictions for the evolution of plasticity in Bicyclus. Our study provides one of the most compelling cases showing links between seasonality and phenotypic plasticity on a macroevolutionary scale and the potential role of plasticity in facilitating the colonization of novel environments.

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.].

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

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

Seasonal environments drive convergent evolution of a faster pace-of-life in tropical butterflies.

New ecological niches that may arise due to climate change can trigger diversification, but their colonisation often requires adaptations in a suite of life-history traits. We test this hypothesis in species-rich Mycalesina butterflies that have undergone parallel radiations in Africa, Asia, and Madagascar. First, our ancestral state reconstruction of habitat preference, using c. 85% of extant species, revealed that early forest-linked lineages began to invade seasonal savannahs during the late Miocene-Pliocene. Second, rearing replicate pairs of forest and savannah species from the African and Malagasy radiation in a common garden experiment, and utilising published data from the Asian radiation, demonstrated that savannah species consistently develop faster, have smaller bodies, higher fecundity with an earlier investment in reproduction, and reduced longevity, compared to forest species across all three radiations. We argue that time-constraints for reproduction favoured the evolution of a faster pace-of-life in savannah species that facilitated their persistence in seasonal habitats.

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.

Whole-chromosome hitchhiking driven by a male-killing endosymbiont.

Neo-sex chromosomes are found in many taxa, but the forces driving their emergence and spread are poorly understood. The female-specific neo-W chromosome of the African monarch (or queen) butterfly Danaus chrysippus presents an intriguing case study because it is restricted to a single 'contact zone' population, involves a putative colour patterning supergene, and co-occurs with infection by the male-killing endosymbiont Spiroplasma. We investigated the origin and evolution of this system using whole genome sequencing. We first identify the 'BC supergene', a broad region of suppressed recombination across nearly half a chromosome, which links two colour patterning loci. Association analysis suggests that the genes yellow and arrow in this region control the forewing colour pattern differences between D. chrysippus subspecies. We then show that the same chromosome has recently formed a neo-W that has spread through the contact zone within approximately 2,200 years. We also assembled the genome of the male-killing Spiroplasma, and find that it shows perfect genealogical congruence with the neo-W, suggesting that the neo-W has hitchhiked to high frequency as the male-killer has spread through the population. The complete absence of female crossing-over in the Lepidoptera causes whole-chromosome hitchhiking of a single neo-W haplotype, carrying a single allele of the BC supergene and dragging multiple non-synonymous mutations to high frequency. This has created a population of infected females that all carry the same recessive colour patterning allele, making the phenotypes of each successive generation highly dependent on uninfected male immigrants. Our findings show how hitchhiking can occur between the physically unlinked genomes of host and endosymbiont, with dramatic consequences.

To mate, or not to mate: The evolution of reproductive diapause facilitates insect radiation into African savannahs in the Late Miocene.

Many tropical environments experience cyclical seasonal changes, frequently with pronounced wet and dry seasons, leading to a highly uneven temporal distribution of resources. Short-lived animals inhabiting such environments often show season-specific adaptations to cope with alternating selection pressures. African Bicyclus butterflies show strong seasonal polyphenism in a suite of phenotypic and life-history traits, and their adults are thought to undergo reproductive diapause associated with the lack of available larval host plants during the dry season. Using 3 years of longitudinal field data for three species in Malawi, dissections demonstrated that one forest species reproduces continuously, whereas two savannah species undergo reproductive diapause in the dry season, either with or without pre-diapause mating. Using additional data from field-collected and museum samples, we then documented the same three mating strategies for a further 37 species. Phylogenetic analyses indicated that the ancestral state was a non-diapausing forest species, and that habitat preference and mating strategy evolved in a correlated fashion. Bicyclus butterflies underwent rapid diversification during the Late Miocene, coinciding with expansions into more open savannah habitat. We conclude that the ability to undergo reproductive diapause was a key trait that facilitated colonization and eventual radiation into savannahs in the Late Miocene.

A transposable element insertion is associated with an alternative life history strategy.

Tradeoffs affect resource allocation during development and result in fitness consequences that drive the evolution of life history strategies. Yet despite their importance, we know little about the mechanisms underlying life history tradeoffs. Many species of Colias butterflies exhibit an alternative life history strategy (ALHS) where females divert resources from wing pigment synthesis to reproductive and somatic development. Due to this reallocation, a wing color polymorphism is associated with the ALHS: either yellow/orange or white. Here we map the locus associated with this ALHS in Colias crocea to a transposable element insertion located downstream of the Colias homolog of BarH-1, a homeobox transcription factor. Using CRISPR/Cas9 gene editing, antibody staining, and electron microscopy we find white-specific expression of BarH-1 suppresses the formation of pigment granules in wing scales and gives rise to white wing color. Lipid and transcriptome analyses reveal physiological differences associated with the ALHS. Together, these findings characterize a mechanism for a female-limited ALHS.

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.

Geographic origin and migration phenology of European red admirals (Vanessa atalanta) as revealed by stable isotopes.

BACKGROUND: Long-distance migration has evolved multiple times in different animal taxa. For insect migrants, the complete annual migration cycle covering several thousand kilometres, may be performed by several generations, each migrating part of the distance and reproducing. Different life-cycle stages and preferred orientation may thus, be found along the migration route. For migrating red admirals (Vanessa atalanta) it has been questioned if they reproduce in the most northern part of the range. Here we present migration phenology data from a two-year time series of migrating red admirals captured at Rybachy, Kaliningrad, in the northern part of Europe investigating time for migration, life-history stage (migration, reproduction) as well as site of origin in individual butterflies. METHODS: Red admirals were captured daily at a coastal site during spring, summer and autumn in 2004 and 2005. For the sampled individuals, reproductive status and fuel content were estimated by visual inspection, and hydrogen isotopes (δ 2H) were analysed in wing samples. δ 2H values was compared with samples from two nearby reference sites in Estonia and Poland. RESULTS: Analysis of hydrogen isotopes (δ 2H) in red admiral wings showed that the spring cohort were of a southerly origin, while those caught in August or later in the autumn were from the local region or areas further to the north. All females caught during spring had developing eggs in their abdomen, but no eggs were found in late summer/autumn. There was a male-biased sex ratio during autumn and a difference in lipid content between years. When comparing the isotopic data with inland nearby locations, it was clear that the range of δ 2H values (- 181 to - 78) was wider at Rybachy as compared to the two reference sites in Estonia and Poland (- 174 to - 100). CONCLUSIONS: During spring, migratory female red admirals arrived from the south and were ready to reproduce, while the autumn passage mainly engaged local and more northern individuals carrying large fuel deposits in preparation for long-distance migration. The phenology data suggest that individuals select to migrate in favourable weather conditions and that numbers may differ between years. Future studies should focus on individual sampling at a wide range of sites to reveal differential migration strategies and timing of migration between sexes and populations of migrating butterflies.

Wolbachia in the Genus Bicyclus: a Forgotten Player.

Bicyclus butterflies are key species for studies of wing pattern development, phenotypic plasticity, speciation and the genetics of Lepidoptera. One of the key endosymbionts in butterflies, the alpha-Proteobacterium Wolbachia pipientis, is affecting many of these biological processes; however, Bicyclus butterflies have not been investigated systematically as hosts to Wolbachia. In this study, we screen for Wolbachia infection in several Bicyclus species from natural populations across Africa as well as two laboratory populations. Out of the 24 species tested, 19 were found to be infected, and no double infection was found, but both A- and B-supergroup strains colonise this butterfly group. We also show that many of the Wolbachia strains identified in Bicyclus butterflies belong to the ST19 clonal complex. We discuss the importance of our results in regard to routinely screening for Wolbachia when using Bicyclus butterflies as the study organism of research in eco-evolutionary biology.

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.

Two new species of Bebearia Hemming, 1960, as further evidence of centre of endemism of butterflies in Western Nigeria (Lepidoptera: Nymphalidae: Limenitinae).

Two new endemic butterfly species from the genus Bebearia: B. oshogbo sp. nov. and B. wojtusiaki sp. nov., are described from western Nigeria; B. oshogbo is most closely related to the Guineo-Congolian B. tentyris (Hewitson) and the Upper Guinean B. osyris (Schultze), whereas B. wojtusiaki constitutes a morphological and biogeographic link between the Central African B. plistonax (Hewitson) and the Upper Guinean endemic B. arcadius (Fabricius). The finding of these new species gives further strong evidence that western Nigeria should be recognized as a distinct biogeographic sub-region of West Africa, as the area hosts a substantial number of endemic taxa (listed in the discussion).

Differentiation in putative male sex pheromone components across and within populations of the African butterfly Bicyclus anynana as a potential driver of reproductive isolation.

Sexual traits are often the most divergent characters among closely related species, suggesting an important role of sexual traits in speciation. However, to prove this, we need to show that sexual trait differences accumulate before or during the speciation process, rather than being a consequence of it. Here, we contrast patterns of divergence among putative male sex pheromone (pMSP) composition and the genetic structure inferred from variation in the mitochondrial cytochrome oxidase 1 and nuclear CAD loci in the African butterfly Bicyclus anynana (Butler, 1879) to determine whether the evolution of "pheromonal dialects" occurs before or after the differentiation process. We observed differences in abundance of some shared pMSP components as well as differences in the composition of the pMSP among B. anynana populations. In addition, B. anynana individuals from Kenya displayed differences in the pMSP composition within a single population that appeared not associated with genetic differences. These differences in pMSP composition both between and within B. anynana populations were as large as those found between different Bicyclus species. Our results suggest that "pheromonal dialects" evolved within and among populations of B. anynana and may therefore act as precursors of an ongoing speciation process.

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.

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.

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.

On the fate of seasonally plastic traits in a rainforest butterfly under relaxed selection.

Many organisms display phenotypic plasticity as adaptation to seasonal environmental fluctuations. Often, such seasonal responses entails plasticity of a whole suite of morphological and life-history traits that together contribute to the adaptive phenotypes in the alternative environments. While phenotypic plasticity in general is a well-studied phenomenon, little is known about the evolutionary fate of plastic responses if natural selection on plasticity is relaxed. Here, we study whether the presumed ancestral seasonal plasticity of the rainforest butterfly Bicyclus sanaos (Fabricius, 1793) is still retained despite the fact that this species inhabits an environmentally stable habitat. Being exposed to an atypical range of temperatures in the laboratory revealed hidden reaction norms for several traits, including wing pattern. In contrast, reproductive body allocation has lost the plastic response. In the savannah butterfly, B. anynana (Butler, 1879), these traits show strong developmental plasticity as an adaptation to the contrasting environments of its seasonal habitat and they are coordinated via a common developmental hormonal system. Our results for B. sanaos indicate that such integration of plastic traits - as a result of past selection on expressing a coordinated environmental response - can be broken when the optimal reaction norms for those traits diverge in a new environment.

Identification and biosynthesis of novel male specific esters in the wings of the tropical butterfly, Bicyclus martius sanaos.

Representatives of the highly speciose tropical butterfly genus Bicyclus (Lepidoptera: Nymphalidae) are characterized by morphological differences in the male androconia, a set of scales and hair pencils located on the surface of the wings. These androconia are assumed to be associated with the release of courtship pheromones. In the present study, we report the identification and biosynthetic pathways of several novel esters from the wings of male B. martius sanaos. We found that the volatile compounds in this male butterfly were similar to female-produced moth sex pheromones. Components associated with the male wing androconial areas were identified as ethyl, isobutyl and 2-phenylethyl hexadecanoates and (11Z)-11-hexadecenoates, among which the latter are novel natural products. By topical application of deuterium-labelled fatty acid and amino acid precursors, we found these pheromone candidates to be produced in patches located on the forewings of the males. Deuterium labels from hexadecanoic acid were incorporated into (11Z)-11-hexadecenoic acid, providing experimental evidence of a Δ11-desaturase being active in butterflies. This unusual desaturase was found previously to be involved in the biosynthesis of female-produced sex pheromones of moths. In the male butterflies, both hexadecanoic acid and (11Z)-11-hexadecenoic acid were then enzymatically esterified to form the ethyl, isobutyl and 2-phenylethyl esters, incorporating ethanol, isobutanol, and 2-phenylethanol, derived from the corresponding amino acids L-alanine, L-valine, and L-phenylalanine.