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.

The genome sequence of the Glanville fritillary, Melitaea cinxia (Linnaeus, 1758)

We present a genome assembly from an individual male Melitaea cinxia (the Glanville fritillary; Arthropoda; Insecta; Lepidoptera; Nymphalidae). The genome sequence is 499 megabases in span. The complete assembly is scaffolded into 31 chromosomal pseudomolecules, with the Z sex chromosome assembled. Gene annotation of this assembly on Ensembl has identified 13,666 protein coding genes.

Increasing frequency of low summer precipitation synchronizes dynamics and compromises metapopulation stability in the Glanville fritillary butterfly

Climate change is known to shift species' geographical ranges, phenologies and abundances, but less is known about other population dynamic consequences. Here, we analyse spatio-temporal dynamics of the Glanville fritillary butterfly (Melitaea cinxia) in a network of 4000 dry meadows during 21 years. The results demonstrate two strong, related patterns: the amplitude of year-to-year fluctuations in the size of the metapopulation as a whole has increased, though there is no long-term trend in average abundance; and there is a highly significant increase in the level of spatial synchrony in population dynamics. The increased synchrony cannot be explained by increasing within-year spatial correlation in precipitation, the key environmental driver of population change, or in per capita growth rate. On the other hand, the frequency of drought during a critical life-history stage (early larval instars) has increased over the years, which is sufficient to explain the increasing amplitude and the expanding spatial synchrony in metapopulation dynamics. Increased spatial synchrony has the general effect of reducing long-term metapopulation viability even if there is no change in average metapopulation size. This study demonstrates how temporal changes in weather conditions can lead to striking changes in spatio-temporal population dynamics.

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.

Temperature- and sex-related effects of serine protease alleles on larval development in the Glanville fritillary butterfly

The body reserves of adult Lepidoptera are accumulated during larval development. In the Glanville fritillary butterfly, larger body size increases female fecundity, but in males fast larval development and early eclosion, rather than large body size, increase mating success and hence fitness. Larval growth rate is highly heritable, but genetic variation associated with larval development is largely unknown. By comparing the Glanville fritillary population living in the Åland Islands in northern Europe with a population in Nantaizi in China, within the source of the post-glacial range expansion, we identified candidate genes with reduced variation in Åland, potentially affected by selection under cooler climatic conditions than in Nantaizi. We conducted an association study of larval growth traits by genotyping the extremes of phenotypic trait distributions for 23 SNPs in 10 genes. Three genes in clip-domain serine protease family were associated with larval growth rate, development time and pupal weight. Additive effects of two SNPs in the prophenoloxidase-activating proteinase-3 (ProPO3) gene, related to melanization, showed elevated growth rate in high temperature but reduced growth rate in moderate temperature. The allelic effects of the vitellin-degrading protease precursor gene on development time were opposite in the two sexes, one genotype being associated with long development time and heavy larvae in females but short development time in males. Sexually antagonistic selection is here evident in spite of sexual size dimorphism.

Long-term metapopulation study of the Glanville fritillary butterfly (Melitaea cinxia): survey methods, data management, and long-term population trends

Long-term observational studies conducted at large (regional) spatial scales contribute to better understanding of landscape effects on population and evolutionary dynamics, including the conditions that affect long-term viability of species, but large-scale studies are expensive and logistically challenging to keep running for a long time. Here, we describe the long-term metapopulation study of the Glanville fritillary butterfly (Melitaea cinxia) that has been conducted since 1991 in a large network of 4000 habitat patches (dry meadows) within a study area of 50 by 70 km in the Åland Islands in Finland. We explain how the landscape structure has been described, including definition, delimitation, and mapping of the habitat patches; methods of field survey, including the logistics, cost, and reliability of the survey; and data management using the EarthCape biodiversity platform. We describe the long-term metapopulation dynamics of the Glanville fritillary based on the survey. There has been no long-term change in the overall size of the metapopulation, but the level of spatial synchrony and hence the amplitude of fluctuations in year-to-year metapopulation dynamics have increased over the years, possibly due to increasing frequency of exceptional weather conditions. We discuss the added value of large-scale and long-term population studies, but also emphasize the need to integrate more targeted experimental studies in the context of long-term observational studies. For instance, in the case of the Glanville fritillary project, the long-term study has produced an opportunity to sample individuals for experiments from local populations with a known demographic history. These studies have demonstrated striking differences in dispersal rate and other life-history traits of individuals from newly established local populations (the offspring of colonizers) versus individuals from old, established local populations. The long-term observational study has stimulated the development of metapopulation models and provided an opportunity to test model predictions. This combination of empirical studies and modeling has facilitated the study of key phenomena in spatial dynamics, such as extinction threshold and extinction debt.