Landscape fragmentation overturns classical metapopulation thinking

Habitat loss and isolation caused by landscape fragmentation represent a growing threat to global biodiversity. Existing theory suggests that the process will lead to a decline in metapopulation viability. However, since most metapopulation models are restricted to simple networks of discrete habitat patches, the effects of real landscape fragmentation, particularly in stochastic environments, are not well understood. To close this major gap in ecological theory, we developed a spatially explicit, individual-based model applicable to realistic landscape structures, bridging metapopulation ecology and landscape ecology. This model reproduced classical metapopulation dynamics under conventional model assumptions, but on fragmented landscapes, it uncovered general dynamics that are in stark contradiction to the prevailing views in the ecological and conservation literature. Notably, fragmentation can give rise to a series of dualities: a) positive and negative responses to environmental noise, b) relative slowdown and acceleration in density decline, and c) synchronization and desynchronization of local population dynamics. Furthermore, counter to common intuition, species that interact locally ("residents") were often more resilient to fragmentation than long-ranging "migrants." This set of findings signals a need to fundamentally reconsider our approach to ecosystem management in a noisy and fragmented world.

Dynamics of a host-parasitoid interaction clarified by modelling and DNA sequencing

It has been hypothesised that the 2-year oscillations in abundance of Xestia moths are mediated by interactions with 1-year Ophion parasitoid wasps. We tested this hypothesis by modelling a 35-year time series of Xestia and Ophion from Northern Finland. Additionally, we used DNA barcoding to ascertain the species diversity of Ophion and targeted amplicon sequencing of their gut contents to confirm their larval hosts. Modelling of the time-series data strongly supported the hypothesised host-parasitoid dynamics and that periodic occurrence of Xestia moths is mediated by Ophion. DNA barcodes revealed that Ophion included five species rather than just one while targeted amplicon sequencing verified that Ophion does parasitise Xestia. At least one Ophion species employs 1-year Syngrapha interrogationis as an alternate host, but it did not detectably affect Xestia-Ophion dynamics. We also demonstrate the previously unrecognised complexity of this system due to cryptic parasitoid diversity.

Scaling up the effects of inbreeding depression from individuals to metapopulations

Inbreeding is common in nature, and many laboratory studies have documented that inbreeding depression can reduce the fitness of individuals. Demonstrating the consequences of inbreeding depression on the growth and persistence of populations is more challenging because populations are often regulated by density- or frequency-dependent selection and influenced by demographic and environmental stochasticity. A few empirical studies have shown that inbreeding depression can increase extinction risk of local populations. The importance of inbreeding depression at the metapopulation level has been conjectured based on population-level studies but has not been evaluated. We quantified the impact of inbreeding depression affecting the fitness of individuals on metapopulation persistence in heterogeneous habitat networks of different sizes and habitat configuration in a context of natural butterfly metapopulations. We developed a spatial individual-based simulation model of metapopulations with explicit genetics. We used Approximate Bayesian Computation to fit the model to extensive demographic, genetic and life-history data available for the well-studied Glanville fritillary butterfly (Melitaea cinxia) metapopulations in the Åland islands in SW Finland. We compared 18 semi-independent habitat networks differing in size and fragmentation. The results show that inbreeding is more frequent in small habitat networks, and consequently, inbreeding depression elevates extinction risks in small metapopulations. Metapopulation persistence and neutral genetic diversity maintained in the metapopulations increase with the total habitat amount in and mean patch size of habitat networks. Dispersal and mating behaviour interact with landscape structure to determine how likely it is to encounter kin while looking for mates. Inbreeding depression can decrease the viability of small metapopulations even when they are strongly influenced by stochastic extinction-colonization dynamics and density-dependent selection. The findings from this study support that genetic factors, in addition to demographic factors, can contribute to extinctions of small local populations and also of metapopulations.

Responses of generalist and specialist species to fragmented landscapes

Empirical studies have shown that, unlike species with specialized resource requirements, generalist species may benefit from habitat destruction. We use a family of models to probe the causes of the contrasting responses of these two types of species to habitat destruction. Our approach allows a number of mechanisms to be switched on and off, thereby making it possible to study their marginal and joint effects. Unlike many previous models, we do not assume any intrinsic competitive asymmetry between the species, and we assume pre-emptive rather than displacement competition. Under these assumptions, in the mean-field model the prevalences of all species decrease monotonically with decreasing habitat availability, independently of the degree of specialization. However, in the stochastic and spatial individual-based simulations of the same model, the specialists dominate in landscapes of high quality, whereas generalists thrive in landscapes of intermediate quality; no species persist in very poor landscapes. The same pattern also occurs in a non-spatial stochastic model but not in a deterministic spatial model, showing that demographic stochasticity plays a key role in shaping the outcome of competitive interactions.

Genetic effects on life-history traits in the Glanville fritillary butterfly

BACKGROUND

Adaptation to local habitat conditions may lead to the natural divergence of populations in life-history traits such as body size, time of reproduction, mate signaling or dispersal capacity. Given enough time and strong enough selection pressures, populations may experience local genetic differentiation. The genetic basis of many life-history traits, and their evolution according to different environmental conditions remain however poorly understood.

METHODS

We conducted an association study on the Glanville fritillary butterfly, using material from five populations along a latitudinal gradient within the Baltic Sea region, which show different degrees of habitat fragmentation. We investigated variation in 10 principal components, cofounding in total 21 life-history traits, according to two environmental types, and 33 genetic SNP markers from 15 candidate genes.

RESULTS

We found that nine SNPs from five genes showed strong trend for trait associations (p-values under 0.001 before correction). These associations, yet non-significant after multiple test corrections, with a total number of 1,086 tests, were consistent across the study populations. Additionally, these nine genes also showed an allele frequency difference between the populations from the northern fragmented versus the southern continuous landscape.

DISCUSSION

Our study provides further support for previously described trait associations within the Glanville fritillary butterfly species across different spatial scales. Although our results alone are inconclusive, they are concordant with previous studies that identified these associations to be related to climatic changes or habitat fragmentation within the Åland population.

Significant disparities in allergy prevalence and microbiota between the young people in Finnish and Russian Karelia

BACKGROUND

Atopic allergy has been more common among schoolchildren in Finland, as compared to Russian Karelia. These adjacent regions show one of the most contrasting socio-economical differences in the world.

OBJECTIVE

We explored changes in allergy from school age to young adulthood from 2003 to 2010/2012 in these two areas. The skin and nasal microbiota were also compared.

METHODS

Randomly selected children from Finnish (n = 98) and Russian Karelia (n = 82) were examined in 2003, when the children were 7-11 years of age, and again in 2010 (Finnish Karelia) and 2012 (Russian Karelia). We analysed self-reported allergy symptoms and sensitization to common allergens by serum sIgE values. The skin (volar forearm) and nasal mucosa microbiota, collected in 2012 (aged 15-20 years), identified from DNA samples, were compared with multivariate methods.

RESULTS

Asthma, hay fever, atopic eczema, self-reported rhinitis, as well as atopic sensitization, were threefold to 10-fold more common in Finland, as compared to Russian Karelia. Hay fever and peanut sensitization were almost non-existent in Russia. These patterns remained throughout the 10-year follow-up. Skin microbiota, as well as bacterial and fungal communities in nasal mucosa, was contrastingly different between the populations, best characterized by the diversity and abundance of genus Acinetobacter; more abundant and diverse in Russia. Overall, diversity was significantly higher among Russian subjects (P_skin < 0.0001, P_{nasal-bacteria} < 0.0001 and P_nasal-fungi < 0.01). Allergic diseases were not associated with microbial diversity in Finnish subjects.

CONCLUSIONS AND CLINICAL RELEVANCE

Differences in allergic phenotype, developed in early life, remain between populations. A parallel difference in the composition of skin and nasal microbiota suggests a potential underlying mechanism. Our results also suggest that high abundance and diversity of Acinetobacter might contribute to the low allergy prevalence in Russia. Implications of early-life exposure to Acinetobacter should be further investigated.

Ecological and genetic basis of metapopulation persistence of the Glanville fritillary butterfly in fragmented landscapes

Ecologists are challenged to construct models of the biological consequences of habitat loss and fragmentation. Here, we use a metapopulation model to predict the distribution of the Glanville fritillary butterfly during 22 years across a large heterogeneous landscape with 4,415 small dry meadows. The majority (74%) of the 125 networks into which the meadows were clustered are below the extinction threshold for long-term persistence. Among the 33 networks above the threshold, spatial configuration and habitat quality rather than the pooled habitat area predict metapopulation size and persistence, but additionally allelic variation in a SNP in the gene Phosphoglucose isomerase (Pgi) explains 30% of variation in metapopulation size. The Pgi genotypes are associated with dispersal rate and hence with colonizations and extinctions. Associations between Pgi genotypes, population turnover and metapopulation size reflect eco-evolutionary dynamics, which may be a common feature in species inhabiting patch networks with unstable local dynamics.

" Nature step" to prevent noncommunicable inflammatory disease

The prevention of many diseases has significantly improved by intervening in known risk factors. However, the causes of the increase in allergy and type 1 diabetes are unknown. These diseases are often associated with a low-grade inflammation and immunological imbalance. The lifestyle and environment of urbanized populations have changed causing imbalance in the human normal flora and affecting immune regulation. We discuss everyday factors affecting immune regulation, using allergy as an example. Health may be promoted through the "nature step", by supporting the connection between humans and nature.

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.

Effects of ambient and preceding temperatures and metabolic genes on flight metabolism in the Glanville fritillary butterfly

Flight is essential for foraging, mate searching and dispersal in many insects, but flight metabolism in ectotherms is strongly constrained by temperature. Thermal conditions vary greatly in natural populations and may hence restrict fitness-related activities. Working on the Glanville fritillary butterfly (Melitaea cinxia), we studied the effects of temperature experienced during the first 2 days of adult life on flight metabolism, genetic associations between flight metabolic rate and variation in candidate metabolic genes, and genotype-temperature interactions. The maximal flight performance was reduced by 17% by 2 days of low ambient temperature (15 °C) prior to the flight trial, mimicking conditions that butterflies commonly encounter in nature. A SNP in phosphoglucose isomerase (Pgi) had a significant association on flight metabolic rate in males and a SNP in triosephosphate isomerase (Tpi) was significantly associated with flight metabolic rate in females. In the Pgi SNP, AC heterozygotes had higher flight metabolic rate than AA homozygotes following low preceding temperature, but the trend was reversed following high preceding temperature, consistent with previous results on genotype-temperature interaction for this SNP. We suggest that these results on 2-day old butterflies reflect thermal effect on the maturation of flight muscles. These results highlight the consequences of variation in thermal conditions on the time scale of days, and they contribute to a better understanding of the complex dynamics of flight metabolism and flight-related activities under conditions that are relevant for natural populations living under variable thermal conditions.

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.

Wolbachia Infection in a Natural Parasitoid Wasp Population

The maternally transmitted bacterium Wolbachia pipientis is well known for spreading and persisting in insect populations through manipulation of the fitness of its host. Here, we identify three new Wolbachia pipientis strains, wHho, wHho2 and wHho3, infecting Hyposoter horticola, a specialist wasp parasitoid of the Glanville fritillary butterfly. The wHho strain (ST435) infects about 50% of the individuals in the Åland islands in Finland, with a different infection rate in the two mitochondrial (COI) haplotypes of the wasp. The vertical transmission rate of Wolbachia is imperfect, and lower in the haplotype with lower infection rate, suggesting a fitness trade-off. We found no association of the wHho infection with fecundity, longevity or dispersal ability of the parasitoid host. However, preliminary results convey spatial associations between Wolbachia infection, host mitochondrial haplotype and parasitism of H. horticola by its hyperparasitoid, Mesochorus cf. stigmaticus. We discuss the possibility that Wolbachia infection protects H. horticola against hyperparasitism.

Large-Scale Habitat Corridors for Biodiversity Conservation: A Forest Corridor in Madagascar

In biodiversity conservation, habitat corridors are assumed to increase landscape-level connectivity and to enhance the viability of otherwise isolated populations. While the role of corridors is supported by empirical evidence, studies have typically been conducted at small spatial scales. Here, we assess the quality and the functionality of a large 95-km long forest corridor connecting two large national parks (416 and 311 km²) in the southeastern escarpment of Madagascar. We analyze the occurrence of 300 species in 5 taxonomic groups in the parks and in the corridor, and combine high-resolution forest cover data with a simulation model to examine various scenarios of corridor destruction. At present, the corridor contains essentially the same communities as the national parks, reflecting its breadth which on average matches that of the parks. In the simulation model, we consider three types of dispersers: passive dispersers, which settle randomly around the source population; active dispersers, which settle only in favorable habitat; and gap-avoiding active dispersers, which avoid dispersing across non-habitat. Our results suggest that long-distance passive dispersers are most sensitive to ongoing degradation of the corridor, because increasing numbers of propagules are lost outside the forest habitat. For a wide range of dispersal parameters, the national parks are large enough to sustain stable populations until the corridor becomes severely broken, which will happen around 2065 if the current rate of forest loss continues. A significant decrease in gene flow along the corridor is expected after 2040, and this will exacerbate the adverse consequences of isolation. Our results demonstrate that simulation studies assessing the role of habitat corridors should pay close attention to the mode of dispersal and the effects of regional stochasticity.

Helsinki alert of biodiversity and health

Urban living in built environments, combined with the use of processed water and food, may not provide the microbial stimulation necessary for a balanced development of immune function. Many chronic inflammatory disorders, including allergic, autoimmune, metabolic, and even some behavioural disorders, are linked to alteration in the human commensal microbiota. Sedentary lifestyle is associated with reduced exposure to a broad spectrum of environmental micro-organisms and surplus energy balance, both risk factors of chronic inflammatory disorders. According to the Biodiversity Hypothesis, an environment with diverse macrobiota and microbiota modifies and enriches the human microbiota, which in turn is crucial in the development and maintenance of appropriate immune function. These issues were discussed in the symposium 'Chronic Inflammation, Lifestyle and Environment', held in Helsinki, 20-22 August 2014, under the sponsorship of the Yrjö Jahnsson Foundation. This paper briefly outlines the recent findings in the context of the environment, lifestyle, and health; discusses the forces that undermine immune tolerance in urban environments; and highlights the possibilities to restore broken immune tolerance among urban dwellers, summarizing the main messages in four statements and calling for actions to combat major public health threats.

Green areas around homes reduce atopic sensitization in children

Background

Western lifestyle is associated with high prevalence of allergy, asthma and other chronic inflammatory disorders. To explain this association, we tested the ‘biodiversity hypothesis’, which posits that reduced contact of children with environmental biodiversity, including environmental microbiota in natural habitats, has adverse consequences on the assembly of human commensal microbiota and its contribution to immune tolerance.

Methods

We analysed four study cohorts from Finland and Estonia (n = 1044) comprising children and adolescents aged 0.5–20 years. The prevalence of atopic sensitization was assessed by measuring serum IgE specific to inhalant allergens. We calculated the proportion of five land-use types – forest, agricultural land, built areas, wetlands and water bodies – in the landscape around the homes using the CORINE2006 classification.

Results

The cover of forest and agricultural land within 2–5 km from the home was inversely and significantly associated with atopic sensitization. This relationship was observed for children 6 years of age and older. Land-use pattern explained 20% of the variation in the relative abundance of Proteobacteria on the skin of healthy individuals, supporting the hypothesis of a strong environmental effect on the commensal microbiota.

Conclusions

The amount of green environment (forest and agricultural land) around homes was inversely associated with the risk of atopic sensitization in children. The results indicate that early-life exposure to green environments is especially important. The environmental effect may be mediated via the effect of environmental microbiota on the commensal microbiota influencing immunotolerance.

Acinetobacter species in the skin microbiota protect against allergic sensitization and inflammation

BACKGROUND

The human commensal microbiota interacts in a complex manner with the immune system, and the outcome of these interactions might depend on the immune status of the subject.

OBJECTIVE

Previous studies have suggested a strong allergy-protective effect for Gammaproteobacteria. Here we analyze the skin microbiota, allergic sensitization (atopy), and immune function in a cohort of adolescents, as well as the influence of Acinetobacter species on immune responses in vitro and in vivo.

METHODS

The skin microbiota of the study subjects was identified by using 16S rRNA sequencing. PBMCs were analyzed for baseline and allergen-stimulated mRNA expression. In in vitro assays human monocyte-derived dendritic cells and primary keratinocytes were incubated with Acinetobacter lwoffii. Finally, in in vivo experiments mice were injected intradermally with A lwoffii during the sensitization phase of the asthma protocol, followed by readout of inflammatory parameters.

RESULTS

In healthy subjects, but not in atopic ones, the relative abundance of Acinetobacter species was associated with the expression of anti-inflammatory molecules by PBMCs. Moreover, healthy subjects exhibited a robust balance between anti-inflammatory and TH1/TH2 gene expression, which was related to the composition of the skin microbiota. In cell assays and in a mouse model, Acinetobacter species induced strong TH1 and anti-inflammatory responses by immune cells and skin cells and protected against allergic sensitization and lung inflammation through the skin.

CONCLUSION

These results support the hypothesis that skin commensals play an important role in tuning the balance of TH1, TH2, and anti-inflammatory responses to environmental allergens.

The Glanville fritillary genome retains an ancient karyotype and reveals selective chromosomal fusions in Lepidoptera

Previous studies have reported that chromosome synteny in Lepidoptera has been well conserved, yet the number of haploid chromosomes varies widely from 5 to 223. Here we report the genome (393 Mb) of the Glanville fritillary butterfly (Melitaea cinxia; Nymphalidae), a widely recognized model species in metapopulation biology and eco-evolutionary research, which has the putative ancestral karyotype of n=31. Using a phylogenetic analyses of Nymphalidae and of other Lepidoptera, combined with orthologue-level comparisons of chromosomes, we conclude that the ancestral lepidopteran karyotype has been n=31 for at least 140 My. We show that fusion chromosomes have retained the ancestral chromosome segments and very few rearrangements have occurred across the fusion sites. The same, shortest ancestral chromosomes have independently participated in fusion events in species with smaller karyotypes. The short chromosomes have higher rearrangement rate than long ones. These characteristics highlight distinctive features of the evolutionary dynamics of butterflies and moths.

Butterflies and moths (Lepidoptera) vary in chromosome number. Here, the authors sequence the genome of the Glanville fritillary butterfly, Melitaea cinxia, show it has the ancestral lepidopteran karyotype and provide insight into how chromosomal fusions have shaped karyotype evolution in butterflies and moths.

Heritability of flight and resting metabolic rates in the Glanville fritillary butterfly

Dispersal capacity is a key life-history trait especially in species inhabiting fragmented landscapes. Evolutionary models predict that, given sufficient heritable variation, dispersal rate responds to natural selection imposed by habitat loss and fragmentation. Here, we estimate phenotypic variance components and heritability of flight and resting metabolic rates (RMRs) in an ecological model species, the Glanville fritillary butterfly, in which flight metabolic rate (FMR) is known to correlate strongly with dispersal rate. We modelled a two-generation pedigree with the animal model to distinguish additive genetic variance from maternal and common environmental effects. The results show that FMR is significantly heritable, with additive genetic variance accounting for about 40% of total phenotypic variance; thus, FMR has the potential to respond to selection on dispersal capacity. Maternal influences on flight metabolism were negligible. Heritability of flight metabolism was context dependent, as in stressful thermal conditions, environmentally induced variation dominated over additive genetic effects. There was no heritability in RMR, which was instead strongly influenced by maternal effects. This study contributes to a mechanistic understanding of the evolution of dispersal-related traits, a pressing question in view of the challenges posed to many species by changing climate and fragmentation of natural habitats.

Differences in the aerobic capacity of flight muscles between butterfly populations and species with dissimilar flight abilities

Habitat loss and climate change are rapidly converting natural habitats and thereby increasing the significance of dispersal capacity for vulnerable species. Flight is necessary for dispersal in many insects, and differences in dispersal capacity may reflect dissimilarities in flight muscle aerobic capacity. In a large metapopulation of the Glanville fritillary butterfly in the Åland Islands in Finland, adults disperse frequently between small local populations. Individuals found in newly established populations have higher flight metabolic rates and field-measured dispersal distances than butterflies in old populations. To assess possible differences in flight muscle aerobic capacity among Glanville fritillary populations, enzyme activities and tissue concentrations of the mitochondrial protein Cytochrome-c Oxidase (CytOx) were measured and compared with four other species of Nymphalid butterflies. Flight muscle structure and mitochondrial density were also examined in the Glanville fritillary and a long-distance migrant, the red admiral. Glanville fritillaries from new populations had significantly higher aerobic capacities than individuals from old populations. Comparing the different species, strong-flying butterfly species had higher flight muscle CytOx content and enzymatic activity than short-distance fliers, and mitochondria were larger and more numerous in the flight muscle of the red admiral than the Glanville fritillary. These results suggest that superior dispersal capacity of butterflies in new populations of the Glanville fritillary is due in part to greater aerobic capacity, though this species has a low aerobic capacity in general when compared with known strong fliers. Low aerobic capacity may limit dispersal ability of the Glanville fritillary.

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.

Lep-MAP: fast and accurate linkage map construction for large SNP datasets

MOTIVATION

Current high-throughput sequencing technologies allow cost-efficient genotyping of millions of single nucleotide polymorphisms (SNPs) for hundreds of samples. However, the tools that are currently available for constructing linkage maps are not well suited for large datasets. Linkage maps of large datasets would be helpful in de novo genome assembly by facilitating comprehensive genome validation and refinement by enabling chimeric scaffold detection, as well as in family-based linkage and association studies, quantitative trait locus mapping, analysis of genome synteny and other complex genomic data analyses.

RESULTS

We describe a novel tool, called Lepidoptera-MAP (Lep-MAP), for constructing accurate linkage maps with ultradense genome-wide SNP data. Lep-MAP is fast and memory efficient and largely automated, requiring minimal user interaction. It uses simultaneously data on multiple outbred families and can increase linkage map accuracy by taking into account achiasmatic meiosis, a special feature of Lepidoptera and some other taxa with no recombination in one sex (no recombination in females in Lepidoptera). We demonstrate that Lep-MAP outperforms other methods on real and simulated data. We construct a genome-wide linkage map of the Glanville fritillary butterfly (Melitaea cinxia) with over 40 000 SNPs. The data were generated with a novel in-house SOLiD restriction site-associated DNA tag sequencing protocol, which is described in the online supplementary material.

AVAILABILITY AND IMPLEMENTATION

Java source code under GNU general public license with the compiled classes and the datasets are available from http://sourceforge.net/users/lep-map.

Evolution of camouflage drives rapid ecological change in an insect community

BACKGROUND

Evolutionary change in individual species has been hypothesized to have far-reaching consequences for entire ecological communities, and such coupling of ecological and evolutionary dynamics ("eco-evolutionary dynamics") has been demonstrated for a variety systems. However, the general importance of evolutionary dynamics for ecological dynamics remains unclear. Here, we investigate how spatial patterns of local adaptation in the stick insect Timema cristinae, driven by the interaction between multiple evolutionary processes, structure metapopulations, communities, and multitrophic interactions.

RESULTS

Observations of a wild T. cristinae metapopulation show that locally imperfect camouflage reduces population size and that the effect of such maladaptation is comparable to the effects of more traditional ecological factors, including habitat patch size and host-plant species identity. Field manipulations of local adaptation and bird predation support the hypothesis that maladaptation reduces population size through an increase in bird predation. Furthermore, these field experiments show that maladaptation in T. cristinae and consequent increase in bird predation reduce the pooled abundance and species richness of the co-occurring arthropod community, and ultimately cascade to decrease herbivory on host plants. An eco-evolutionary model of the observational data demonstrates that the demographic cost of maladaptation decreases habitat patch occupancy by T. cristinae but enhances metapopulation-level adaptation.

CONCLUSIONS

The results demonstrate a pervasive effect of ongoing evolution in a spatial context on population and community dynamics. The eco-evolutionary model makes testable predictions about the influence of the spatial configuration of the patch network on metapopulation size and the spatial scale of adaptation.

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.

Species-area relationships and extinctions caused by habitat loss and fragmentation

The species-area relationship (SAR) has been used to predict the numbers of species going extinct due to habitat loss, but other researchers have maintained that SARs overestimate extinctions and instead one should use the endemics-area relationship (EAR) to predict extinctions. Here, we employ spatially explicit simulations of large numbers of species in spatially heterogeneous landscapes to investigate SARs and extinctions in a dynamic context. The EAR gives the number of species going extinct immediately after habitat loss, but typically many other species have unviable populations in the remaining habitat and go extinct soon afterwards. We conclude that the EAR underestimates extinctions due to habitat loss, the continental SAR (with slope ~0.1 or somewhat less) gives a good approximation of short-term extinctions, while the island SAR calculated for discrete fragments of habitat (with slope ~0.25) predicts the long-term extinctions. However, when the remaining area of land-covering habitat such as forest is roughly less than 20% of the total landscape and the habitat is highly fragmented, all current SARs underestimate extinction rate. We show how the 'fragmentation effect' can be incorporated into a predictive SAR model. When the remaining habitat is highly fragmented, an effective way to combat the fragmentation effect is to aggregate habitat fragments into clusters rather than to place them randomly across the landscape.

Butterfly survival on an isolated island by improved grip

On small isolated islands, natural selection is expected to reduce the dispersal capacity of organisms, as short distances do not require a high rate of dispersal, which might lead to accidental emigration from the population. In addition, individuals foregoing the high cost of maintaining flight capacity may instead allocate resources to other functions. However, in butterflies and many other insects, flight is necessary not only for dispersal but also for most other activities. A weakly flying individual would probably do worse and have an elevated rather than reduced probability of accidental emigration. Here, we report results consistent with the hypothesis that a butterfly population on an isolated island, instead of having lost its flight capacity, has evolved better grip to resist the force of wind and to avoid being blown off the island. Our study suggests that local adaptation has occurred in this population in spite of its very small size (Ne ∼ 100), complete isolation, low genetic variation and high genetic load.

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.

A long life in the fast lane: positive association between peak metabolic rate and lifespan in a butterfly

High peak metabolic rates may provide a performance advantage, but it may also entail a physiological cost. A long-held assumption is that high mass-specific energy expenditure is associated with short lifespan. To examine the relationship between energy expenditure and lifespan we asked two questions. First, do individuals have a consistent rate of metabolism throughout their life? Second, is metabolic rate correlated with lifespan? We analysed the repeatability of measurements of resting (RMR) and peak flight metabolic rate (MR(peak)) throughout the life of the Glanville fritillary butterfly (Melitaea cinxia). Measurements of MR(peak) showed significant repeatability. Senescence occurred only shortly before death. RMR showed a U-shaped relationship with age and very low repeatability. Intraspecific association between metabolic rates and lifespan was tested under three conditions: in the laboratory, under field conditions and in a laboratory experiment with repeated flight treatments. There was a significant correlation between MR(peak) and lifespan in all three experiments, but the correlation was positive, not negative. RMR was not correlated with lifespan. Both MR(peak) and lifespan may reflect physiological condition and therefore be positively correlated. Individuals with a large resource pool may be able to invest in mechanisms that slow down ageing. Individuals with high metabolic capacity may also possess adaptations against ageing. Molecular polymorphism in the gene phosphoglucose isomerase (Pgi) was significantly associated with both MR(peak) and lifespan, and may have coevolved with defence mechanisms against senescence. Generalisations such as 'live fast, die young' may be too simple to explain the complex processes affecting ageing and lifespan.

Variation in a Host-Parasitoid Interaction across Independent Populations

Antagonistic relationships between parasitoids and their insect hosts involve multiple traits and are shaped by their ecological and evolutionary context. The parasitoid wasp Cotesia melitaearum and its host butterfly Melitaea cinxia occur in several locations around the Baltic sea, with differences in landscape structure, population sizes and the histories of the populations. We compared the virulence of the parasitoid and the susceptibility of the host from five populations in a reciprocal transplant-style experiment using the progeny of five independent host and parasitoid individuals from each population. The host populations showed significant differences in the rate of encapsulation and parasitoid development rate. The parasitoid populations differed in brood size, development rate, pupal size and adult longevity. Some trait differences depended on specific host-parasitoid combinations, but neither species performed systematically better or worse in experiments involving local versus non-local populations of the other species. Furthermore, individuals from host populations with the most recent common ancestry did not perform alike, and there was no negative effect due to a history of inbreeding in the parasitoid. The complex pattern of variation in the traits related to the vulnerability of the host and the ability of the parasitoid to exploit the host may reflect multiple functions of the traits that would hinder simple local adaptation.

High genetic load in an old isolated butterfly population

We investigated inbreeding depression and genetic load in a small (N(e) ∼ 100) population of the Glanville fritillary butterfly (Melitaea cinxia), which has been completely isolated on a small island [Pikku Tytärsaari (PT)] in the Baltic Sea for at least 75 y. As a reference, we studied conspecific populations from the well-studied metapopulation in the Åland Islands (ÅL), 400 km away. A large population in Saaremaa, Estonia, was used as a reference for estimating genetic diversity and N(e). We investigated 58 traits related to behavior, development, morphology, reproductive performance, and metabolism. The PT population exhibited high genetic load (L = 1 - W(PT)/W(ÅL)) in a range of fitness-related traits including adult weight (L = 0.12), flight metabolic rate (L = 0.53), egg viability (L = 0.37), and lifetime production of eggs in an outdoor population cage (L = 0.70). These results imply extensive fixation of deleterious recessive mutations, supported by greatly reduced diversity in microsatellite markers and immediate recovery (heterosis) of egg viability and flight metabolic rate in crosses with other populations. There was no significant inbreeding depression in most traits due to one generation of full-sib mating. Resting metabolic rate was significantly elevated in PT males, which may be related to their short lifespan (L = 0.25). The demographic history and the effective size of the PT population place it in the part of the parameter space in which models predict mutation accumulation. This population exemplifies the increasingly common situation in fragmented landscapes, in which small and completely isolated populations are vulnerable to extinction due to high genetic load.

Environmental biodiversity, human microbiota, and allergy are interrelated

Rapidly declining biodiversity may be a contributing factor to another global megatrend--the rapidly increasing prevalence of allergies and other chronic inflammatory diseases among urban populations worldwide. According to the "biodiversity hypothesis," reduced contact of people with natural environmental features and biodiversity may adversely affect the human commensal microbiota and its immunomodulatory capacity. Analyzing atopic sensitization (i.e., allergic disposition) in a random sample of adolescents living in a heterogeneous region of 100 × 150 km, we show that environmental biodiversity in the surroundings of the study subjects' homes influenced the composition of the bacterial classes on their skin. Compared with healthy individuals, atopic individuals had lower environmental biodiversity in the surroundings of their homes and significantly lower generic diversity of gammaproteobacteria on their skin. The functional role of the gram-negative gammaproteobacteria is supported by in vitro measurements of expression of IL-10, a key anti-inflammatory cytokine in immunologic tolerance, in peripheral blood mononuclear cells. In healthy, but not in atopic, individuals, IL-10 expression was positively correlated with the abundance of the gammaproteobacterial genus Acinetobacter on the skin. These results raise fundamental questions about the consequences of biodiversity loss for both allergic conditions and public health in general.

Eco-evolutionary dynamics of dispersal in spatially heterogeneous environments

Evolutionary changes in natural populations are often so fast that the evolutionary dynamics may influence ecological population dynamics and vice versa. Here we construct an eco-evolutionary model for dispersal by combining a stochastic patch occupancy metapopulation model with a model for changes in the frequency of fast-dispersing individuals in local populations. We test the model using data on allelic variation in the gene phosphoglucose isomerase (Pgi), which is strongly associated with dispersal rate in the Glanville fritillary butterfly. Population-specific measures of immigration and extinction rates and the frequency of fast-dispersing individuals among the immigrants explained 40% of spatial variation in Pgi allele frequency among 97 local populations. The model clarifies the roles of founder events and gene flow in dispersal evolution and resolves a controversy in the literature about the consequences of habitat loss and fragmentation on the evolution of dispersal.

Origin and Diversification of Dung Beetles in Madagascar

Madagascar has a rich fauna of dung beetles (Scarabaeinae and Aphodiinae) withalmost 300 species described to date. Like most other taxa in Madagascar, dung beetles exhibit an exceptionally high level of endemism (96% of the species). Here,we review the current knowledge of the origin and diversification of Malagasy dung beetles. Based on molecular phylogenies, the extant dung beetles originate from eight colonizations, of which four have given rise to extensive radiations. These radiations have occurred in wet forests, while the few extant species in the less successfulradiations occur in open and semi-open habitats. We discuss the likely mechanisms of speciation and the ecological characteristics of the extant communities, emphasizing the role of adaptation along environmental gradients and allopatric speciation in generating the exceptionally high beta diversity in Malagasy dung beetles. Phylogeographic analyses of selected species reveal complex patterns with evidence for genetic introgression between old taxa. The introduction of cattle to Madagascar 1500 years ago created a new abundant resource, onto which a few species haveshifted and thereby been able to greatly expand their geographical ranges.