A Candidate Gene in an Ecological Model Species: Phosphoglucose Isomerase (Pgi) in the Glanville Fritillary Butterfly (Melitaea cinxia)

Parity-specific differences in spatial genetics and dispersal in the common lizard

Dispersal is a key demographic parameter that plays an important role in determining spatial population dynamics and genetic structure. Linking differences in dispersal patterns to life-history traits is often confounded by inconsistent environmental pressures experienced by different populations. To explore the relationship between dispersal and life history, we focus on a site where oviparous and viviparous lineages of the common lizard (Zootoca vivipara) are found adjacent to each other. We take advantage of this shared environment to investigate parity-specific dispersal patterns using high-resolution, individual-level spatial-genetic autocorrelation and population genomic approaches (11,726 single nucleotide polymorphisms; 293 oviparous and 310 viviparous individuals). We found isolation-by-distance patterns to be present in both the oviparous and viviparous populations. Density was 2.5 times higher in the oviparous population than the viviparous one, though heterozygosity and genetic diversity measures were similar in the two populations. We found marked differences in the extent of genetic neighbourhoods between the lineages, with the viviparous population showing both dispersal (σ) and spatial-genetic autocorrelation (Moran's I) at 2-fold greater geographic distances than the oviparous population. We found clear evidence of male-biased dispersal from genetic estimates in the viviparous population. In the oviparous population, evidence of male-biased dispersal was weak or absent. These differences are likely to be closely linked to specific requirements of the alternative reproductive strategies and may be the demographic consequences of mother-offspring interactions. Fine-scale geographic and individual-level measures are essential to understanding parity mode differences at microevolutionary scales and to better identifying their ecological and evolutionary impacts.

Traits Underlying Experimentally Evolved Dispersal Behavior in Tribolium castaneum

Dispersal is an important behavior in many animals, with profound effects on individual fitness and the evolutionary trajectories of populations. This is especially true within taxa with particular life-history strategies, for example those that exploit ephemeral habitat. Further, dispersal is commonly seen to be part of behavioral syndromes - suites of traits that covary across behavioral contexts. The red flour beetle, Tribolium castaneum (Coleoptera, Tenebrionidae), is a major post-harvest crop pest responsible for large losses through the infestation of stored grain. In this system dispersal is known to have a strong genetic basis and differential artificial selection on dispersal traits produces strong phenotypic divergence. However, it is unknown which traits are able to rapidly evolve to produce these results, or which behavioral components underlie differences in dispersal. Using replicate lines of T. castaneum previously selected for divergent dispersal behavior, we test for correlated activity and movement patterns, morphology and substrate surface use. We find robustly repeatable associations between the dispersal phenotype and higher activity, straighter paths, larger body size (but not relative leg length) and increased tendency to remain at the surface of fodder. Together our results suggest that dispersal is part of a syndrome of traits in T. castaneum, and must be treated as such when considering the evolution of dispersal in this system, and in attempting to predict and control its spread.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10905-024-09862-x.

The genetic basis of dispersal in a vertebrate metapopulation

Dispersal affects evolutionary processes by changing population size and genetic composition, influencing the viability and persistence of populations. Investigating which mechanisms underlie variation in dispersal phenotypes and whether populations harbour adaptive potential for dispersal is crucial to understanding the eco-evolutionary dynamics of this important trait. Here, we investigate the genetic architecture of dispersal among successfully recruited individuals in an insular metapopulation of house sparrows. We use an extensive long-term individual-based ecological data set and high-density single-nucleotide polymorphism (SNP) genotypes for over 2500 individuals. We conducted a genome-wide association study (GWAS), and found a relationship between dispersal probability and a SNP located near genes known to regulate circadian rhythm, glycogenesis and exercise performance, among other functions. However, this SNP only explained 3.8% of variance, suggesting that dispersal is a polygenic trait. We then used an animal model to estimate heritable genetic variation (σA 2 ), which composes 10% of the total variation in dispersal probability. Finally, we investigated differences in σA 2 across populations occupying ecologically relevant habitat types (farm vs. non-farm) using a genetic groups animal model. We found different adaptive potentials across habitats, with higher mean breeding value, σA 2 , and heritability for the habitat presenting lower dispersal rates, suggesting also different roles of environmental variation. Our results suggest a complex genetic architecture of dispersal and demonstrate that adaptive potential may be environment dependent in key eco-evolutionary traits. The eco-evolutionary implications of such environment dependence and consequent spatial variation are likely to become ever more important with the increased fragmentation and loss of suitable habitats for many natural populations.

Genetic architecture of dispersal behaviour in the post-harvest pest and model organism Tribolium castaneum

Dispersal behaviour is an important aspect of the life-history of animals. However, the genetic architecture of dispersal-related traits is often obscure or unknown, even in well studied species. Tribolium castaneum is a globally significant post-harvest pest and established model organism, yet studies of its dispersal have shown ambiguous results and the genetic basis of this behaviour remains unresolved. We combine experimental evolution and agent-based modelling to investigate the number of loci underlying dispersal in T. castaneum, and whether the trait is sex-linked. Our findings demonstrate rapid evolution of dispersal behaviour under selection. We find no evidence of sex-biases in the dispersal behaviour of the offspring of crosses, supporting an autosomal genetic basis of the trait. Moreover, simulated data approximates experimental data under simulated scenarios where the dispersal trait is controlled by one or few loci, but not many loci. Levels of dispersal in experimentally inbred lines, compared with simulations, indicate that a single locus model is not well supported. Taken together, these lines of evidence support an oligogenic architecture underlying dispersal in Tribolium castaneum. These results have implications for applied pest management and for our understanding of the evolution of dispersal in the coleoptera, the world's most species-rich order.

Investigating the genetic basis of vertebrate dispersal combining RNA-seq, RAD-seq and quantitative genetics

Although animal dispersal is known to play key roles in ecological and evolutionary processes such as colonization, population extinction and local adaptation, little is known about its genetic basis, particularly in vertebrates. Untapping the genetic basis of dispersal should deepen our understanding of how dispersal behaviour evolves, the molecular mechanisms that regulate it and link it to other phenotypic aspects in order to form the so-called dispersal syndromes. Here, we comprehensively combined quantitative genetics, genome-wide sequencing and transcriptome sequencing to investigate the genetic basis of natal dispersal in a known ecological and evolutionary model of vertebrate dispersal: the common lizard, Zootoca vivipara. Our study supports the heritability of dispersal in semi-natural populations, with less variation attributable to maternal and natal environment effects. In addition, we found an association between natal dispersal and both variation in the carbonic anhydrase (CA10) gene, and in the expression of several genes (TGFB2, SLC6A4, NOS1) involved in central nervous system functioning. These findings suggest that neurotransmitters (serotonin and nitric oxide) are involved in the regulation of dispersal and shaping dispersal syndromes. Several genes from the circadian clock (CRY2, KCTD21) were also differentially expressed between disperser and resident lizards, supporting that the circadian rhythm, known to be involved in long-distance migration in other taxa, might affect dispersal as well. Since neuronal and circadian pathways are relatively well conserved across vertebrates, our results are likely to be generalisable, and we therefore encourage future studies to further investigate the role of these pathways in shaping dispersal in vertebrates.

Alternative developmental and transcriptomic responses to host plant water limitation in a butterfly metapopulation

Predicting how climate change affects biotic interactions poses a challenge. Plant-insect herbivore interactions are particularly sensitive to climate change, as climate-induced changes in plant quality cascade into the performance of insect herbivores. Whereas the immediate survival of herbivore individuals depends on plastic responses to climate change-induced nutritional stress, long-term population persistence via evolutionary adaptation requires genetic variation for these responses. To assess the prospects for population persistence under climate change, it is therefore crucial to characterize response mechanisms to climate change-induced stressors, and quantify their variability in natural populations. Here, we test developmental and transcriptomic responses to water limitation-induced host plant quality change in a Glanville fritillary butterfly (Melitaea cinxia) metapopulation. We combine nuclear magnetic resonance spectroscopy on the plant metabolome, larval developmental assays and an RNA sequencing analysis of the larval transcriptome. We observed that responses to feeding on water-limited plants, in which amino acids and aromatic compounds are enriched, showed marked variation within the metapopulation, with individuals of some families performing better on control and others on water-limited plants. The transcriptomic responses were concordant with the developmental responses: families exhibiting opposite developmental responses also produced opposite transcriptomic responses (e.g. in growth-associated transcripts). The divergent responses in both larval development and transcriptome are associated with differences between families in amino acid catabolism and storage protein production. The results reveal intrapopulation variability in plasticity, suggesting that the Finnish M. cinxia metapopulation harbours potential for buffering against drought-induced changes in host plant quality.

Effects of environment and genotype on dispersal differ across departure, transfer and settlement in a butterfly metapopulation

Active dispersal is driven by extrinsic and intrinsic factors at the three stages of departure, transfer and settlement. Most empirical studies capture only one stage of this complex process, and knowledge of how much can be generalized from one stage to another remains unknown. Here we use genetic assignment tests to reconstruct dispersal across 5 years and 232 habitat patches of a Glanville fritillary butterfly (Melitaea cinxia) metapopulation. We link individual dispersal events to weather, landscape structure, size and quality of habitat patches, and individual genotype to identify the factors that influence the three stages of dispersal and post-settlement survival. We found that nearly all tested factors strongly affected departure probabilities, but that the same factors explained very little variation in realized dispersal distances. Surprisingly, we found no effect of dispersal distance on post-settlement survival. Rather, survival was influenced by weather conditions, quality of the natal habitat patch, and a strong interaction between genotype and occupancy status of the settled habitat patch, with more mobile genotypes having higher survival as colonists rather than as immigrants. Our work highlights the multi-causality of dispersal and that some dispersal costs can only be understood by considering extrinsic and intrinsic factors and their interaction across the entire dispersal process.

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 effect of summer drought on the predictability of local extinctions in a butterfly metapopulation

The ecological impacts of extreme climatic events on population dynamics and community composition are profound and predominantly negative. Using extensive data of an ecological model system, we tested whether predictions from ecological models remain robust when environmental conditions are outside the bounds of observation. We observed a 10-fold demographic decline of the Glanville fritillary butterfly (Melitaea cinxia) metapopulation on the Åland islands, Finland in the summer of 2018 and used climatic and satellite data to demonstrate that this year was an anomaly with low climatic water balance values and low vegetation productivity indices across Åland. Population growth rates were strongly associated with spatiotemporal variation in climatic water balance. Covariates shown previously to affect the extinction probability of local populations in this metapopulation were less informative when populations were exposed to severe drought during the summer months. Our results highlight the unpredictable responses of natural populations to extreme climatic events.

The March of the Beetles: Epistatic Components Dominate Divergence in Dispersal Tendency in Tribolium castaneum

The genetic underpinnings of traits are rarely simple. Most traits of interest are instead the product of multiple genes acting in concert to determine the phenotype. This is particularly true for behavioral traits, like dispersal. Our investigation focuses on the genetic architecture of dispersal tendency in the red flour beetle, Tribolium castaneum. We used artificial selection to generate lines with either high or low dispersal tendency. Our populations responded quickly in the first generations of selection and almost all replicates had higher dispersal tendency in males than in females. These selection lines were used to create a total of 6 additional lines: F1 and reciprocal F1, as well as 4 types of backcrosses. We estimated the composite genetic effects that contribute to divergence in dispersal tendency among lines using line cross-analysis. We found variation in the dispersal tendency of our lines was best explained by autosomal additive and 3 epistatic components. Our results indicate that dispersal tendency is heritable, but much of the divergence in our selection lines was due to epistatic effects. These results are consistent with other life-history traits that are predicted to maintain more epistatic variance than additive variance and highlight the potential for epistatic variation to act as an adaptive reserve that may become visible to selection when a population is subdivided.

Alleles in metabolic and oxygen-sensing genes are associated with antagonistic pleiotropic effects on life history traits and population fitness in an ecological model insect

Genes with opposing effects on fitness at different life stages are the mechanistic basis for evolutionary theories of aging and life history. Examples come from studies of mutations in model organisms, but there is little knowledge of genetic bases of life history tradeoffs in natural populations. Here, we test the hypothesis that alleles affecting oxygen sensing in Glanville fritillary butterflies have opposing effects on larval versus adult fitness-related traits. Intermediate-frequency alleles in Succinate dehydrogenase d, and to a lesser extent Hypoxia inducible factor 1α, are associated in larvae with variation in metabolic rate and activation of the hypoxia inducible factor (HIF) pathway, which affects tracheal development and delivery of oxygen to adult flight muscles. A dominant Sdhd allele is likely to cause antagonistic pleiotropy for fitness through its opposing effects on larval metabolic and growth rate versus adult flight and dispersal, and may have additional effects arising from sensitivity to low-iron host plants. Prior results in Glanville fritillaries indicate that fitness of alleles in Sdhd and another antagonistically pleiotropic metabolic gene, Phosphoglucose isomerase, depend strongly on the size and distribution of host plant patches. Hence, these intermediate-frequency alleles are involved in ecoevolutionary dynamics involving life history tradeoffs.

The physiology of movement

Movement, from foraging to migration, is known to be under the influence of the environment. The translation of environmental cues to individual movement decision making is determined by an individual's internal state and anticipated to balance costs and benefits. General body condition, metabolic and hormonal physiology mechanistically underpin this internal state. These physiological determinants are tightly, and often genetically linked with each other and hence central to a mechanistic understanding of movement. We here synthesise the available evidence of the physiological drivers and signatures of movement and review (1) how physiological state as measured in its most coarse way by body condition correlates with movement decisions during foraging, migration and dispersal, (2) how hormonal changes underlie changes in these movement strategies and (3) how these can be linked to molecular pathways. We reveale that a high body condition facilitates the efficiency of routine foraging, dispersal and migration. Dispersal decision making is, however, in some cases stimulated by a decreased individual condition. Many of the biotic and abiotic stressors that induce movement initiate a physiological cascade in vertebrates through the production of stress hormones. Movement is therefore associated with hormone levels in vertebrates but also insects, often in interaction with factors related to body or social condition. The underlying molecular and physiological mechanisms are currently studied in few model species, and show -in congruence with our insights on the role of body condition- a central role of energy metabolism during glycolysis, and the coupling with timing processes during migration. Molecular insights into the physiological basis of movement remain, however, highly refractory. We finalise this review with a critical reflection on the importance of these physiological feedbacks for a better mechanistic understanding of movement and its effects on ecological dynamics at all levels of biological organization.

Landscape permeability and individual variation in a dispersal-linked gene jointly determine genetic structure in the Glanville fritillary butterfly

There is now clear evidence that species across a broad range of taxa harbor extensive heritable variation in dispersal. While studies suggest that this variation can facilitate demographic outcomes such as range expansion and invasions, few have considered the consequences of intraspecific variation in dispersal for the maintenance and distribution of genetic variation across fragmented landscapes. Here, we examine how landscape characteristics and individual variation in dispersal combine to predict genetic structure using genomic and spatial data from the Glanville fritillary butterfly. We used linear and latent factor mixed models to identify the landscape features that best predict spatial sorting of alleles in the dispersal-related gene phosphoglucose isomerase (Pgi). We next used structural equation modeling to test if variation in Pgi mediated gene flow as measured by Fst at putatively neutral loci. In a year when the population was recovering following a large decline, individuals with a genotype associated with greater dispersal ability were found at significantly higher frequencies in populations isolated by water and forest, and these populations showed lower levels of genetic differentiation at neutral loci. These relationships disappeared in the next year when metapopulation density was high, suggesting that the effects of individual variation are context dependent. Together our results highlight that (1) more complex aspects of landscape structure beyond just the configuration of habitat can be important for maintaining spatial variation in dispersal traits and (2) that individual variation in dispersal plays a key role in maintaining genetic variation across fragmented landscapes.

Frontiers in Metapopulation Biology: The Legacy of Ilkka Hanski

This review of metapopulation biology has a special focus on Professor Ilkka Hanski's (1953–2016) research. Hanski made seminal contributions to both empirical and theoretical metapopulation biology throughout his scientific career. Hanski's early research focused on ecological aspects of metapopulation biology, in particular how the spatial structure of a landscape influences extinction thresholds and how habitat loss and fragmentation can result in extinction debt. Hanski then used the Glanville fritillary system as a natural laboratory within which he studied genetic and evolutionary processes, such as the influence of inbreeding on extinction risk and variation in selection for dispersal traits generated by landscape variation. During the last years of his career, Hanski's work was in the forefront of the rapidly developing field of eco-evolutionary dynamics. Hanski was a pioneer in showing how molecular-level underpinnings of trait variation can explain why evolutionary change can occur rapidly in natural populations and how these changes can subsequently influence ecological dynamics.

Environmental and genetic control of cold tolerance in the Glanville fritillary butterfly

Thermal tolerance has a major effect on individual fitness and species distributions and can be determined by genetic variation and phenotypic plasticity. We investigate the effects of developmental and adult thermal conditions on cold tolerance, measured as chill coma recovery (CCR) time, during the early and late adult stage in the Glanville fritillary butterfly. We also investigate the genetic basis of cold tolerance by associating CCR variation with polymorphisms in candidate genes that have a known role in insect physiology. Our results demonstrate that a cooler developmental temperature leads to reduced cold tolerance in the early adult stage, whereas cooler conditions during the adult stage lead to increased cold tolerance. This suggests that adult acclimation, but not developmental plasticity, of adult cold tolerance is adaptive. This could be explained by the ecological conditions the Glanville fritillary experiences in the field, where temperature during early summer, but not spring, is predictive of thermal conditions during the butterfly's flight season. In addition, an amino acid polymorphism (Ala-Glu) in the gene flightin, which has a known function in insect flight and locomotion, was associated with CCR. These amino acids have distinct biochemical properties and may thus affect protein function and/or structure. To our knowledge, our study is the first to link genetic variation in flightin to cold tolerance, or thermal adaptation in general.

Enzyme polymorphism, oxygen and injury: a lipidomic analysis of flight-induced oxidative damage in a succinate dehydrogenase d (Sdhd)-polymorphic insect

When active tissues receive insufficient oxygen to meet metabolic demand, succinate accumulates and has two fundamental effects: it causes ischemia-reperfusion injury while also activating the hypoxia-inducible factor pathway (HIF). The Glanville fritillary butterfly (Melitaea cinxia) possesses a balanced polymorphism in Sdhd, shown previously to affect HIF pathway activation and tracheal morphology and used here to experimentally test the hypothesis that variation in succinate dehydrogenase affects oxidative injury. We stimulated butterflies to fly continuously in a respirometer (3 min duration), which typically caused episodes of exhaustion and recovery, suggesting a potential for cellular injury from hypoxia and reoxygenation in flight muscles. Indeed, flight muscle from butterflies flown on consecutive days had lipidome profiles similar to those of rested paraquat-injected butterflies, but distinct from those of rested untreated butterflies. Many butterflies showed a decline in flight metabolic rate (FMR) on day 2, and there was a strong inverse relationship between the ratio of day 2 to day 1 FMR and the abundance of sodiated adducts of phosphatidylcholines and co-enzyme Q (CoQ). This result is consistent with elevation of sodiated lipids caused by disrupted intracellular ion homeostasis in mammalian tissues after hypoxia-reperfusion. Butterflies carrying the Sdhd M allele had a higher abundance of lipid markers of cellular damage, but the association was reversed in field-collected butterflies, where focal individuals typically flew for seconds at a time rather than continuously. These results indicate that Glanville fritillary flight muscles can be injured by episodes of high exertion, but injury severity appears to be determined by an interaction between SDH genotype and behavior (prolonged versus intermittent flight).

Genetics of dispersal

Dispersal is a process of central importance for the ecological and evolutionary dynamics of populations and communities, because of its diverse consequences for gene flow and demography. It is subject to evolutionary change, which begs the question, what is the genetic basis of this potentially complex trait? To address this question, we (i) review the empirical literature on the genetic basis of dispersal, (ii) explore how theoretical investigations of the evolution of dispersal have represented the genetics of dispersal, and (iii) discuss how the genetic basis of dispersal influences theoretical predictions of the evolution of dispersal and potential consequences.

Dispersal has a detectable genetic basis in many organisms, from bacteria to plants and animals. Generally, there is evidence for significant genetic variation for dispersal or dispersal‐related phenotypes or evidence for the micro‐evolution of dispersal in natural populations. Dispersal is typically the outcome of several interacting traits, and this complexity is reflected in its genetic architecture: while some genes of moderate to large effect can influence certain aspects of dispersal, dispersal traits are typically polygenic. Correlations among dispersal traits as well as between dispersal traits and other traits under selection are common, and the genetic basis of dispersal can be highly environment‐dependent.

By contrast, models have historically considered a highly simplified genetic architecture of dispersal. It is only recently that models have started to consider multiple loci influencing dispersal, as well as non‐additive effects such as dominance and epistasis, showing that the genetic basis of dispersal can influence evolutionary rates and outcomes, especially under non‐equilibrium conditions. For example, the number of loci controlling dispersal can influence projected rates of dispersal evolution during range shifts and corresponding demographic impacts. Incorporating more realism in the genetic architecture of dispersal is thus necessary to enable models to move beyond the purely theoretical towards making more useful predictions of evolutionary and ecological dynamics under current and future environmental conditions. To inform these advances, empirical studies need to answer outstanding questions concerning whether specific genes underlie dispersal variation, the genetic architecture of context‐dependent dispersal phenotypes and behaviours, and correlations among dispersal and other traits.