Improved chromosome-level genome assembly of the Glanville fritillary butterfly (Melitaea cinxia) integrating Pacific Biosciences long reads and a high-density linkage map.

BACKGROUND: The Glanville fritillary (Melitaea cinxia) butterfly is a model system for metapopulation dynamics research in fragmented landscapes. Here, we provide a chromosome-level assembly of the butterfly's genome produced from Pacific Biosciences sequencing of a pool of males, combined with a linkage map from population crosses. RESULTS: The final assembly size of 484 Mb is an increase of 94 Mb on the previously published genome. Estimation of the completeness of the genome with BUSCO indicates that the genome contains 92-94% of the BUSCO genes in complete and single copies. We predicted 14,810 genes using the MAKER pipeline and manually curated 1,232 of these gene models. CONCLUSIONS: The genome and its annotated gene models are a valuable resource for future comparative genomics, molecular biology, transcriptome, and genetics studies on this species.

Multidimensional plasticity in the Glanville fritillary butterfly: larval performance is temperature, host and family specific.

Variation in environmental conditions during development can lead to changes in life-history traits with long-lasting effects. Here, we study how variation in temperature and host plant (i.e. the consequences of potential maternal oviposition choices) affects a suite of life-history traits in pre-diapause larvae of the Glanville fritillary butterfly. We focus on offspring survival, larval growth rates and relative fat reserves, and pay specific attention to intraspecific variation in the responses (G × E × E). Globally, thermal performance and survival curves varied between diets of two host plants, suggesting that host modifies the temperature impact, or vice versa. Additionally, we show that the relative fat content has a host-dependent, discontinuous response to developmental temperature. This implies that a potential switch in resource allocation, from more investment in growth at lower temperatures to storage at higher temperatures, is dependent on the larval diet. Interestingly, a large proportion of the variance in larval performance is explained by differences among families, or interactions with this variable. Finally, we demonstrate that these family-specific responses to the host plant remain largely consistent across thermal environments. Together, the results of our study underscore the importance of paying attention to intraspecific trait variation in the field of evolutionary ecology.

Spatial and temporal genetic structure at the fourth trophic level in a fragmented landscape.

A fragmented habitat becomes increasingly fragmented for species at higher trophic levels, such as parasitoids. To persist, these species are expected to possess life-history traits, such as high dispersal, that facilitate their ability to use resources that become scarce in fragmented landscapes. If a specialized parasitoid disperses widely to take advantage of a sparse host, then the parasitoid population should have lower genetic structure than the host. We investigated the temporal and spatial genetic structure of a hyperparasitoid (fourth trophic level) in a fragmented landscape over 50 × 70 km, using microsatellite markers, and compared it with the known structures of its host parasitoid, and the butterfly host which lives as a classic metapopulation. We found that population genetic structure decreases with increasing trophic level. The hyperparasitoid has fewer genetic clusters (K = 4), than its host parasitoid (K = 15), which in turn is less structured than the host butterfly (K = 27). The genetic structure of the hyperparasitoid also shows temporal variation, with genetic differentiation increasing due to reduction of the population size, which reduces the effective population size. Overall, our study confirms the idea that specialized species must be dispersive to use a fragmented host resource, but that this adaptation has limits.

Oxygen and energy availability interact to determine flight performance in the Glanville fritillary butterfly.

Flying insects have the highest known mass-specific demand for oxygen, which makes it likely that reduced availability of oxygen might limit sustained flight, either instead of or in addition to the limitation due to metabolite resources. The Glanville fritillary butterfly (Melitaea cinxia) occurs as a large metapopulation in which adult butterflies frequently disperse between small local populations. Here, we examine how the interaction between oxygen availability and fuel use affects flight performance in the Glanville fritillary. Individuals were flown under either normoxic (21 kPa O2) or hypoxic (10 kPa O2) conditions and their flight metabolism was measured. To determine resource use, levels of circulating glucose, trehalose and whole-body triglyceride were recorded after flight. Flight performance was significantly reduced in hypoxic conditions. When flown under normoxic conditions, we observed a positive correlation among individuals between post-flight circulating trehalose levels and flight metabolic rate, suggesting that low levels of circulating trehalose constrains flight metabolism. To test this hypothesis experimentally, we measured the flight metabolic rate of individuals injected with a trehalase inhibitor. In support of the hypothesis, experimental butterflies showed significantly reduced flight metabolic rate, but not resting metabolic rate, in comparison to control individuals. By contrast, under hypoxia there was no relationship between trehalose and flight metabolic rate. Additionally, in this case, flight metabolic rate was reduced in spite of circulating trehalose levels that were high enough to support high flight metabolic rate under normoxic conditions. These results demonstrate a significant interaction between oxygen and energy availability for the control of flight performance.

Transcriptome analysis reveals signature of adaptation to landscape fragmentation.

We characterize allelic and gene expression variation between populations of the Glanville fritillary butterfly (Melitaea cinxia) from two fragmented and two continuous landscapes in northern Europe. The populations exhibit significant differences in their life history traits, e.g. butterflies from fragmented landscapes have higher flight metabolic rate and dispersal rate in the field, and higher larval growth rate, than butterflies from continuous landscapes. In fragmented landscapes, local populations are small and have a high risk of local extinction, and hence the long-term persistence at the landscape level is based on frequent re-colonization of vacant habitat patches, which is predicted to select for increased dispersal rate. Using RNA-seq data and a common garden experiment, we found that a large number of genes (1,841) were differentially expressed between the landscape types. Hexamerin genes, the expression of which has previously been shown to have high heritability and which correlate strongly with larval development time in the Glanville fritillary, had higher expression in fragmented than continuous landscapes. Genes that were more highly expressed in butterflies from newly-established than old local populations within a fragmented landscape were also more highly expressed, at the landscape level, in fragmented than continuous landscapes. This result suggests that recurrent extinctions and re-colonizations in fragmented landscapes select a for specific expression profile. Genes that were significantly up-regulated following an experimental flight treatment had higher basal expression in fragmented landscapes, indicating that these butterflies are genetically primed for frequent flight. Active flight causes oxidative stress, but butterflies from fragmented landscapes were more tolerant of hypoxia. We conclude that differences in gene expression between the landscape types reflect genomic adaptations to landscape fragmentation.

Plastic larval development in a butterfly has complex environmental and genetic causes and consequences for population dynamics.

1. In insects, the length of larval development time typically influences adult body size and individual fitness, and hence development time can be expected to respond in an adaptive manner to variation in environmental conditions. In the wild, larval growth may be influenced by individual condition, which can be affected by population-level parameters such as population density and abundance and quality of resources. 2. We sampled larvae of the Glanville fritillary butterfly (Melitaea cinxia) from 514 local populations across a large metapopulation before the winter diapause and reared the larvae in common garden conditions after diapause. Here, we report that small post-diapause larvae prolonged their development via an extra larval instar, apparently to compensate for their 'bad start' after diapause. The number of instars was additionally a plastic response to environmental conditions, as the frequency of the extra instar increased under cooler thermal conditions. 3. The benefit of the extra instar is clear, as it allows individuals to develop into larger adults, but the cost is delayed adult eclosion, which is likely to select against the extra instar especially in males, in which early eclosion is critical for mating success. In support of this, the frequency of the extra instar was significantly lower in males (7%) than in females (42%). 4. Polymorphisms in three genes, serpin-1, vitellin-degrading protease precursor and phosphoglucose isomerase, which are known to influence development in insects, were associated with the occurrence of the extra instar. 5. At the level of local populations, the frequency of the extra instar was higher in newly established populations than that in old local ones, possibly reflecting maternal effects, as new populations are often established by females with heavy investment in dispersal. The frequency of the extra instar in turn correlated with the change in population size over 1 year and the risk of local extinction in the natural metapopulation of the Glanville fritillary. 6. Our results highlight the importance of the physiological condition of individuals in shaping subsequent life-history events and even population dynamics.

Life history of the Glanville fritillary butterfly in fragmented versus continuous landscapes.

Habitat loss and fragmentation threaten the long-term viability of innumerable species of plants and animals. At the same time, habitat fragmentation may impose strong natural selection and lead to evolution of life histories with possible consequences for demographic dynamics. The Baltic populations of the Glanville fritillary butterfly (Melitaea cinxia) inhabit regions with highly fragmented habitat (networks of small dry meadows) as well as regions with extensive continuous habitat (calcareous alvar grasslands). Here, we report the results of common garden studies on butterflies originating from two highly fragmented landscapes (FL) in Finland and Sweden and from two continuous landscapes (CL) in Sweden and Estonia, conducted in a large outdoor cage (32 by 26 m) and in the laboratory. We investigated a comprehensive set of 51 life-history traits, including measures of larval growth and development, flight performance, and adult reproductive behavior. Seventeen of the 51 traits showed a significant difference between fragmented versus CL. Most notably, the growth rate of postdiapause larvae and several measures of flight capacity, including flight metabolic rate, were higher in butterflies from fragmented than CL. Females from CL had shorter intervals between consecutive egg clutches and somewhat higher life-time egg production, but shorter longevity, than females from FL. These results are likely to reflect the constant opportunities for oviposition in females living in continuous habitats, while the more dispersive females from FL allocate more resources to dispersal capacity at the cost of egg maturation rate. This study supports theoretical predictions about small population sizes and high rate of population turnover in fragmented habitats selecting for increased rate of dispersal, but the results also indicate that many other life-history traits apart from dispersal are affected by the degree of habitat fragmentation.

Diet quality can play a critical role in defense efficacy against parasitoids and pathogens in the Glanville fritillary (Melitaea cinxia).

Numerous herbivorous insect species sequester noxious chemicals from host plants that effectively defend against predators, and against parasitoids and pathogens. Sequestration of these chemicals may be expensive and involve a trade off with other fitness traits. Here, we tested this hypothesis. We reared Glanville fritillary butterfly (Melitaea cinxia L.) larvae on plant diets containing low- and high-levels of iridoid glycosides (IGs) (mainly aucubin and catalpol) and tested: 1) whether IGs affect the herbivore's defense against parasitoids (measured as encapsulation rate) and bacterial pathogens (measured as herbivore survival); 2) whether parasitoid and bacterial defenses interact; and 3) whether sequestration of the plant's defense chemicals incurs any life history costs. Encapsulation rates were stronger when there were higher percentages of catalpol in the diet. Implanted individuals had greater amounts of IGs in their bodies as adults. This suggests that parasitized individuals may sequester more IGs, increase their feeding rate after parasitism, or that there is a trade off between detoxification efficiency and encapsulation rate. Larval survival after bacterial infection was influenced by diet, but probably not by diet IG content, as changes in survival did not correlate linearly with the levels of IGs in the diet. However, M. cinxia larvae with good encapsulation abilities were better defended against bacteria. We did not find any life history costs of diet IG concentration for larvae. These results suggest that the sequestering of plant defense chemicals can help herbivorous insects to defend against parasitoids.

Significant effects of Pgi genotype and body reserves on lifespan in the Glanville fritillary butterfly.

Individuals with a particular variant of the gene phosphoglucose isomerase (Pgi) have been shown to have superior dispersal capacity and fecundity in the Glanville fritillary butterfly (Melitaea cinxia), raising questions about the mechanisms that maintain polymorphism in this gene in the field. Here, we investigate how variation in the Pgi genotype affects female and male life history under controlled conditions. The most striking effect is the longer lifespan of genotypes with high dispersal capacity, especially in non-reproducing females. Butterflies use body reserves for somatic maintenance and reproduction, but different resources (in thorax versus abdomen) are used under dissimilar conditions, with some interactions with the Pgi genotype. These results indicate life-history trade-offs that involve resource allocation and genotypexenvironment interactions, and these trade-offs are likely to contribute to the maintenance of Pgi polymorphism in the natural populations.