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.

Macroevolutionary shifts of WntA function potentiate butterfly wing-pattern diversity.

Butterfly wing patterns provide a rich comparative framework to study how morphological complexity develops and evolves. Here we used CRISPR/Cas9 somatic mutagenesis to test a patterning role for WntA, a signaling ligand gene previously identified as a hotspot of shape-tuning alleles involved in wing mimicry. We show that WntA loss-of-function causes multiple modifications of pattern elements in seven nymphalid butterfly species. In three butterflies with a conserved wing-pattern arrangement, WntA is necessary for the induction of stripe-like patterns known as symmetry systems and acquired a novel eyespot activator role specific to Vanessa forewings. In two Heliconius species, WntA specifies the boundaries between melanic fields and the light-color patterns that they contour. In the passionvine butterfly Agraulis, WntA removal shows opposite effects on adjacent pattern elements, revealing a dual role across the wing field. Finally, WntA acquired a divergent role in the patterning of interveinous patterns in the monarch, a basal nymphalid butterfly that lacks stripe-like symmetry systems. These results identify WntA as an instructive signal for the prepatterning of a biological system of exuberant diversity and illustrate how shifts in the deployment and effects of a single developmental gene underlie morphological change.

Flight-induced transgenerational maternal effects influence butterfly offspring performance during times of drought.

Maternal condition can generate resource-related maternal effects through differential egg provisioning that can negatively affect offspring performance especially when offspring growth occurs in stressful or sub-optimal environments. Using the Speckled Wood butterfly, Pararge aegeria (L.) we tested the hypothesis that repeated periods of intensive flight during female oviposition affects egg provisioning and reduces offspring performance when larval development occurs under stressful conditions on drought stressed host plants. We investigated whether (after controlling for egg size) maternal age and flight treatment resulted in changes in egg provisioning and whether this contributed to variation in offspring traits across life stages. Age-related changes in maternal condition were found to generate resource-related maternal effects that influenced offspring traits across all life stages. Flight-induced changes in maternal egg provisioning were found to have direct consequences for offspring development in the egg and larval stages. There were significant interactive effects between maternal age and flight on larval development and growth. Compared to offspring from forced flight mothers, offspring from control (no forced flight) mothers that hatched from eggs laid early in the oviposition period (i.e. by younger mothers) had shorter larval development times and heavier pupal masses, suggesting that offspring from mothers in relatively good condition may be able to buffer some of the costs associated with growth on drought stressed host plants. Our multi-factor study demonstrates the importance of considering the various, and often interacting, mechanisms by which maternal effects may influence offspring performance in stressful environments.

Vertically transmitted rhabdoviruses are found across three insect families and have dynamic interactions with their hosts.

A small number of free-living viruses have been found to be obligately vertically transmitted, but it remains uncertain how widespread vertically transmitted viruses are and how quickly they can spread through host populations. Recent metagenomic studies have found several insects to be infected with sigma viruses (Rhabdoviridae). Here, we report that sigma viruses that infect Mediterranean fruit flies (Ceratitis capitata), Drosophila immigrans, and speckled wood butterflies (Pararge aegeria) are all vertically transmitted. We find patterns of vertical transmission that are consistent with those seen in Drosophila sigma viruses, with high rates of maternal transmission, and lower rates of paternal transmission. This mode of transmission allows them to spread rapidly in populations, and using viral sequence data we found the viruses in D. immigrans and C. capitata had both recently swept through host populations. The viruses were common in nature, with mean prevalences of 12% in C. capitata, 38% in D. immigrans and 74% in P. aegeria We conclude that vertically transmitted rhabdoviruses may be widespread in a broad range of insect taxa, and that these viruses can have dynamic interactions with their hosts.

Ancient expansion of the hox cluster in lepidoptera generated four homeobox genes implicated in extra-embryonic tissue formation.

Gene duplications within the conserved Hox cluster are rare in animal evolution, but in Lepidoptera an array of divergent Hox-related genes (Shx genes) has been reported between pb and zen. Here, we use genome sequencing of five lepidopteran species (Polygonia c-album, Pararge aegeria, Callimorpha dominula, Cameraria ohridella, Hepialus sylvina) plus a caddisfly outgroup (Glyphotaelius pellucidus) to trace the evolution of the lepidopteran Shx genes. We demonstrate that Shx genes originated by tandem duplication of zen early in the evolution of large clade Ditrysia; Shx are not found in a caddisfly and a member of the basally diverging Hepialidae (swift moths). Four distinct Shx genes were generated early in ditrysian evolution, and were stably retained in all descendent Lepidoptera except the silkmoth which has additional duplications. Despite extensive sequence divergence, molecular modelling indicates that all four Shx genes have the potential to encode stable homeodomains. The four Shx genes have distinct spatiotemporal expression patterns in early development of the Speckled Wood butterfly (Pararge aegeria), with ShxC demarcating the future sites of extraembryonic tissue formation via strikingly localised maternal RNA in the oocyte. All four genes are also expressed in presumptive serosal cells, prior to the onset of zen expression. Lepidopteran Shx genes represent an unusual example of Hox cluster expansion and integration of novel genes into ancient developmental regulatory networks.

Unscrambling butterfly oogenesis.

BACKGROUND: Butterflies are popular model organisms to study physiological mechanisms underlying variability in oogenesis and egg provisioning in response to environmental conditions. Nothing is known, however, about; the developmental mechanisms governing butterfly oogenesis, how polarity in the oocyte is established, or which particular maternal effect genes regulate early embryogenesis. To gain insights into these developmental mechanisms and to identify the conserved and divergent aspects of butterfly oogenesis, we analysed a de novo ovarian transcriptome of the Speckled Wood butterfly Pararge aegeria (L.), and compared the results with known model organisms such as Drosophila melanogaster and Bombyx mori. RESULTS: A total of 17306 contigs were annotated, with 30% possibly novel or highly divergent sequences observed. Pararge aegeria females expressed 74.5% of the genes that are known to be essential for D. melanogaster oogenesis. We discuss the genes involved in all aspects of oogenesis, including vitellogenesis and choriogenesis, plus those implicated in hormonal control of oogenesis and transgenerational hormonal effects in great detail. Compared to other insects, a number of significant differences were observed in; the genes involved in stem cell maintenance and differentiation in the germarium, establishment of oocyte polarity, and in several aspects of maternal regulation of zygotic development. CONCLUSIONS: This study provides valuable resources to investigate a number of divergent aspects of butterfly oogenesis requiring further research. In order to fully unscramble butterfly oogenesis, we also now also have the resources to investigate expression patterns of oogenesis genes under a range of environmental conditions, and to establish their function.

Exploring sub-lethal effects of exposure to a nucleopolyhedrovirus in the speckled wood (Pararge aegeria) butterfly.

This study investigated the sub-lethal effects of larval exposure to baculovirus on host life history and wing morphological traits using a model system, the speckled wood butterfly Pararge aegeria (L.) and the virus Autographa californica nucleopolyhedrovirus. Males and females showed similar responses to the viral infection. Infection significantly reduced larval growth rate, whilst an increase in development time allowed the critical mass for pupation to be attained. There was no direct effect of viral infection on the wing morphological traits examined. There was, however, an indirect effect of resisting infection; larvae that took longer to develop had reduced resource investment in adult flight muscle mass.

Wing shape and its influence on the outcome of territorial contests in the damselfly Calopteryx virgo.

Male mating success is often determined by territory ownership and traits associated with successful territory defense. Empirical studies have shown that the territory owner wins the majority of fights with challenging males. Several physical and physiological traits have been found to correlate with resource holding potential. In addition, in aerial insects, wing design may also have a strong influence on resource holding potential, since it determines efficiency and precision during flight. However, this possibility has not yet been thoroughly evaluated using the modern technique of geometric morphometrics to analyze shape. Therefore, this study examined whether wing shape affects the outcome of male-male contests in the territorial damselfly, Calopteryx virgo (L.) (Odonata: Calopterygidae). Wing shape and also traditional flight-related morphological measures were compared between 27 pairs of winners and losers from experimental territorial contests. Contrary to expectations, there were no differences between winners and losers in all studied wing traits (shape, length, width, total surface, aspect ratio, and wing loading). However, highly significant differences in wing shape and size were detected between the fore- and hindwing. It is currently not known how these differences relate to flight performance, since previous biomechanical studies in damselflies assumed fore- and hindwings to have an identical planform.

Maternal effects, flight versus fecundity trade-offs, and offspring immune defence in the speckled wood butterfly, Pararge aegeria.

BACKGROUND: Maternal condition can generate resource-related maternal effects through differential egg provisioning, and can greatly affect offspring performance. In the present study, the speckled wood butterfly Pararge aegeria (L.) was used to investigate whether (after controlling for egg size) maternal age, and increased flight during the oviposition period, resulted in changes in egg provisioning and whether this contributed to variation in offspring performance, i) early in development (egg stage and early post-hatching development), and ii) later in larval development after being exposed to the model viral pathogen system; the baculovirus Autographa californica multinucleocapsid nucleopolyhedrovirus (AcMNPV). RESULTS: Age-related changes in maternal egg provisioning were observed to influence egg stage development only. Flight-induced changes in maternal egg provisioning had direct consequences for offspring growth and survival across each life stage from egg to adulthood; offspring from forced flight mothers had lower larval masses and longer development times. Offspring with lower larval masses also had reduced survival after exposure to the viral pathogen. CONCLUSION: The present study demonstrates that a change in maternal provisioning as a result of increased flight during the oviposition period has the potential to exert non-genetic cross-generational fitness effects in P. aegeria. This could have important consequences for population dynamics, particularly in fragmented anthropogenic landscapes.

Male-male competition, female mate choice and their interaction: determining total sexual selection.

Empirical studies of sexual selection typically focus on one of the two mechanisms of sexual selection without integrating these into a description of total sexual selection, or study total sexual selection without quantifying the contributions of all of the mechanisms of sexual selection. However, this can provide an incomplete or misleading view of how sexually selected traits evolve if the mechanisms of sexual selection are opposing or differ in form. Here, we take a two-fold approach to advocate a direction for future studies of sexual selection. We first show how a quantitative partitioning and examination of sexual selection mechanisms can inform by identifying illustrative studies that describe both male-male competition and female mate choice acting on the same trait. In our sample, the most common trait where this occurred was body size, and selection was typically linear. We found that male-male competition and female mate choice can be reinforcing or opposing, although the former is most common in the literature. The mechanisms of sexual selection can occur simultaneously or sequentially, and we found they were more likely to be opposing when the mechanisms operated sequentially. The degree and timing that these mechanisms interact have important implications for the operation of sexual selection and needs to be considered in designing studies. Our examples highlight where empirical data are needed. We especially lack standardized measures of the form and strength of selection imposed by each mechanism of sexual selection and how they combine to determine total sexual selection. Secondly, using quantitative genetic principles, we outline how the selection imposed by individual mechanisms can be measured and combined to estimate the total strength and form of sexual selection. We discuss the evolutionary consequences of combining the mechanisms of sexual selection and interpreting total sexual selection. We suggest how this approach may result in empirical progress in the field of sexual selection.

Genetic differentiation between resistance phenotypes in the phytophagous flea beetle, Phyllotreta nemorum.

The flea beetle Phyllotreta nemorum L. (Coleoptera: Chrysomelidae) is genetically polymorphic for resistance against the defences of one of its host plants, Barbarea vulgaris R.Br. (Brassicales: Brassicaceae). Whereas resistant flea beetles are able to use B. vulgaris as well as other cruciferous plants as food, non-resistant beetles cannot survive on B. vulgaris. This limitation to host plant use of non-resistant beetles could potentially lead to asymmetric gene flow and some degree of genetic isolation between the different resistance-genotypes. Therefore, we studied the extent of genetic differentiation at neutral allozyme loci between samples of flea beetles that were collected at different locations and first tested for resistance phenotype. Since earlier work has shown a weak, but significant, effect of geographical distance between the samples on their genetic differentiation, in the present study variation at the neutral allozyme loci in P. nemorum was partitioned between geographical distance and resistance-phenotype. Both sources independently contributed statistically significantly to population differentiation. Thus, there appears to be a limitation to genetic exchange between the resistant and non-resistant flea beetles when corrections are made for their geographic differentiation. This is consistent with the presence of some degree of host race formation in this flea beetle.

The effects of insecticides on butterflies - A review.

Pesticides, in particular insecticides, can be very beneficial but have also been found to have harmful side effects on non-target insects. Butterflies play an important role in ecosystems, are well monitored and are recognised as good indicators of environmental health. The amount of information already known about butterfly ecology and the increased availability of genomes make them a very valuable model for the study of non-target effects of pesticide usage. The effects of pesticides are not simply linear, but complex through their interactions with a large variety of biotic and abiotic factors. Furthermore, these effects manifest themselves at a variety of levels, from the molecular to metapopulation level. Research should therefore aim to dissect these complex effects at a number of levels, but as we discuss in this review, this is seldom if ever done in butterflies. We suggest that in order dissect the complex effects of pesticides on butterflies we need to integrate detailed molecular studies, including characterising sequence variability of relevant target genes, with more classical evolutionary ecology; from direct toxicity tests on individual larvae in the laboratory to field studies that consider the potentiation of pesticides by ecologically relevant environmental biotic and abiotic stressors. Such integration would better inform population-level responses across broad geographical scales and provide more in-depth information about the non-target impacts of pesticides.

Water Stress Affects Development Time but Not Takeoff Performance in the Butterfly Pararge aegeria.

Most organisms are limited in the amount and type of resources they are able to extract from the environment. The juvenile environment is particularly important in this regard, as conditions over ontogeny can influence the adult phenotype. Whole-organism performance traits, such as locomotion, are susceptible to such environmental effects, yet the specific biotic and abiotic factors driving performance plasticity have received little attention. We tested whether speckled wood Pararge aegeria L. butterflies reared under conditions of water stress exhibited poorer flight morphology and performance than control individuals. Despite large differences in mortality between treatments, we found no effects of water stress treatment on takeoff performance and only minor treatment effects on flight morphology. However, butterflies reared on water-stressed diets exhibited both significantly greater mortality and longer development times than did control individuals. Pararge aegeria larvae may compensate for this stress by prolonging development, resulting in similar realized performance capacities at least in takeoff performance in surviving adult butterflies; other measures of flight performance remain to be considered. Alternatively, the adult phenotype may be insulated from environmental effects at the larval stage in these insects.

Studying Oogenesis in a Non-model Organism Using Transcriptomics: Assembling, Annotating, and Analyzing Your Data.

This chapter provides a guide to processing and analyzing RNA-Seq data in a non-model organism. This approach was implemented for studying oogenesis in the Speckled Wood Butterfly Pararge aegeria. We focus in particular on how to perform a more informative primary annotation of your non-model organism by implementing our multi-BLAST annotation strategy. We also provide a general guide to other essential steps in the next-generation sequencing analysis workflow. Before undertaking these methods, we recommend you familiarize yourself with command line usage and fundamental concepts of database handling. Most of the operations in the primary annotation pipeline can be performed in Galaxy (or equivalent standalone versions of the tools) and through the use of common database operations (e.g. to remove duplicates) but other equivalent programs and/or custom scripts can be implemented for further automation.

AFLP markers for the R-gene in the flea beetle, Phyllotreta nemorum, conferring resistance to defenses in Barbarea vulgaris.

A so-called R-gene renders the yellow-striped flea beetle Phyllotreta nemorum L. (Coleoptera: Chrysomelidae: Alticinae) resistant to the defenses of the yellow rocket Barbarea vulgaris R.Br. (Brassicacea) and enables it to use it as a host plant in Denmark. In this study, genetic markers for an autosomal R-gene, inherited as a single, dominant locus in flea beetles from the Danish locality "Kvaerkeby" are described, and a genetic linkage map around this particular R-gene is constructed, using the technique of AFLP (Amplified Fragment Length Polymorphism).

Divergent RNA Localisation Patterns of Maternal Genes Regulating Embryonic Patterning in the Butterfly Pararge aegeria.

The maternal effect genes responsible for patterning the embryo along the antero-posterior (AP) axis are broadly conserved in insects. The precise function of these maternal effect genes is the result of the localisation of their mRNA in the oocyte. The main developmental mechanisms involved have been elucidated in Drosophila melanogaster, but recent studies have shown that other insect orders often diverge in RNA localisation patterns. A recent study has shown that in the butterfly Pararge aegeria the distinction between blastodermal embryonic (i.e. germ band) and extra-embryonic tissue (i.e. serosa) is already specified in the oocyte during oogenesis in the ovariole, long before blastoderm cellularisation. To examine the extent by which a female butterfly specifies and patterns the AP axis within the region fated to be the germ band, and whether she specifies a germ plasm, we performed in situ hybridisation experiments on oocytes in P. aegeria ovarioles and on early embryos. RNA localisation of the following key maternal effect genes were investigated: caudal (cad), orthodenticle (otd), hunchback (hb) and four nanos (nos) paralogs, as well as TDRD7 a gene containing a key functional domain (OST-HTH/LOTUS) shared with oskar. TDRD7 was mainly confined to the follicle cells, whilst hb was exclusively zygotically transcribed. RNA of some of the nos paralogs, otd and cad revealed complex localisation patterns within the cortical region prefiguring the germ band (i.e. germ cortex). Rather interestingly, otd was localised within and outside the anterior of the germ cortex. Transcripts of nos-O formed a distinct granular ring in the middle of the germ cortex possibly prefiguring the region where germline stem cells form. These butterfly RNA localisation patterns are highly divergent with respect to other insects, highlighting the diverse ways in which different insect orders maternally regulate early embryogenesis of their offspring.

A Diversity of Conserved and Novel Ovarian MicroRNAs in the Speckled Wood (Pararge aegeria).

microRNAs (miRNAs) are important regulators of animal development and other processes, and impart robustness to living systems through post-transcriptional regulation of specific mRNA transcripts. It is postulated that newly emergent miRNAs are generally expressed at low levels and with spatiotemporally restricted expression domains, thus minimising effects of spurious targeting on animal transcriptomes. Here we present ovarian miRNA transcriptome data for two geographically distinct populations of the Speckled Wood butterfly (Pararge aegeria). A total of 74 miRNAs were identified, including 11 newly discovered and evolutionarily-young miRNAs, bringing the total of miRNA genes known from P. aegeria up to 150. We find a positive correlation between miRNA age and expression level. A common set of 55 miRNAs are expressed in both populations. From this set, we identify seven that are consistently either ovary-specific or highly upregulated in ovaries relative to other tissues. This 'ovary set' includes miRNAs with known contributions to ovarian function in other insect species with similar ovaries and mode of oogenesis, including miR-989 and miR-2763, plus new candidates for ovarian function. We also note that conserved miRNAs are overrepresented in the ovary relative to the whole body.

Why Small Is Beautiful: Wing Colour Is Free from Thermoregulatory Constraint in the Small Lycaenid Butterfly, Polyommatus icarus.

We examined the roles of wing melanisation, weight, and basking posture in thermoregulation in Polyommatus Icarus, a phenotypically variable and protandrous member of the diverse Polyommatinae (Lycaenidae). Under controlled experimental conditions, approximating to marginal environmental conditions for activity in the field (= infrequent flight, long duration basking periods), warming rates are maximised with fully open wings and maximum body temperatures are dependent on weight. Variation in wing melanisation within and between sexes has no effect on warming rates; males and females which differ in melanisation had similar warming rates. Posture also affected cooling rates, consistent with cooling being dependent on convective heat loss. We hypothesise that for this small sized butterfly, melanisation has little or no effect on thermoregulation. This may be a factor contributing to the diversity of wing colours in the Polyommatinae. Because of the importance of size for thermoregulation in this small butterfly, requirements for attaining a suitable size to confer thermal stability in adults may also be a factor influencing larval feeding rates, development time and patterns of voltinism. Our findings indicate that commonly accepted views of the importance of melanisation, posture and size to thermoregulation, developed using medium and large sized butterflies, are not necessarily applicable to small sized butterflies.

Female choice depends on size but not symmetry of dorsal eyespots in the butterfly Bicyclus anynana.

The eyespots on the ventral wings of Bicyclus anynana butterflies are exposed when at rest and interact with predators. Those on the dorsal surface are not exposed in this way, and may be involved in courtship and mate choice. In this study, we examined whether the size and fluctuating asymmetry (FA) of dorsal eyespots are reliable signals of male quality. High developmental stability is considered to result in low FA, and to be associated with high quality. Individuals of high quality are predicted to produce sexually selected traits that are large and symmetrical, at a relatively low cost. In this study, we manipulated eyespot development to uncouple eyespot size and FA in order to examine their independent roles in signalling to the female. Individual females in cages were given the choice between two or three males differing in eyespot traits. The results indicate that although size per se of the eyespots is used as a signal, FA and wing size are not. We discuss the use of FA in studies of sexual selection and aspects of sexual selection on dorsal eyespot size.

Cooperative breeding favors maternal investment in size over number of eggs in spiders.

The transition to cooperative breeding may alter maternal investment strategies depending on density of breeders, extent of reproductive skew, and allo-maternal care. Change in optimal investment from solitary to cooperative breeding can be investigated by comparing social species with nonsocial congeners. We tested two hypotheses in a mainly semelparous system: that social, cooperative breeders, compared to subsocial, solitarily breeding congeners, (1) lay fewer and larger eggs because larger offspring compete better for limited resources and become reproducers; (2) induce egg size variation within clutches as a bet-hedging strategy to ensure that some offspring become reproducers. Within two spider genera, Anelosimus and Stegodyphus, we compared species from similar habitats and augmented the results with a mini-meta-analysis of egg numbers depicted in phylogenies. We found that social species indeed laid fewer, larger eggs than subsocials, while egg size variation was low overall, giving no support for bet-hedging. We propose that the transition to cooperative breeding selects for producing few, large offspring because reproductive skew and high density of breeders and young create competition for resources and reproduction. Convergent evolution has shaped maternal strategies similarly in phylogenetically distant species and directed cooperatively breeding spiders to invest in quality rather than quantity of offspring.