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 (Pskin < 0.0001, Pnasal-bacteria < 0.0001 and Pnasal-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.

Larval growth rate is associated with the composition of the gut microbiota in the Glanville fritillary butterfly.

The rapidly increasing body of literature on commensal microbiota has revealed a large phylotypic and functional diversity of microbes associated with vertebrates and invertebrates. In insects, the gut microbiota plays a role in digestion and metabolism of the host as well as protects the host against pathogens. In the study reported here, we sampled gut microbiota of the larvae of the Glanville fritillary butterfly (Melitaea cinxia). The larvae were collected from the field or reared in the laboratory. This butterfly has two host plant species, Plantago lanceolata and Veronica spicata, and the host plant species is known from previous studies to influence larval growth rate. However, our results demonstrate that about 50 % of the variation in larval growth rate can be attributed to the effect of the gut microbial composition plus the joint effect of microbiota and the host plant species, while host plant species alone makes no significant contribution. Our results support previous studies showing that diet influences the gut microbiota but, more unexpectedly, that the composition of the gut microbiota significantly influences larval growth rate. We suggest that host plant effects on larval growth and development observed in many previous studies may be mediated via the gut microbiota. While we measured the growth rate only in laboratory-reared larvae, the similarity of the gut microbial composition between samples from field-collected and laboratory-reared larvae suggests that the results can be generalized to natural 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.

Hunt for the origin of allergy - comparing the Finnish and Russian Karelia.

The Finnish and Russian Karelia are adjacent areas in northern Europe, socio-economically distinct but geoclimatically similar. The Karelia Allergy Study was commenced in 1998 to characterize the allergy profiles in the two areas. Allergy prevalence had increased in Finland since the early 1960s, but the situation in Russia was unknown. The key finding was that allergic symptoms and diseases were systematically more common in Finnish children and adults than in their Russian counterparts. For example, in the early 2000s, hay fever in school children was almost non-existent in Russian Karelia, and only 2% were sensitized to birch pollen compared with 27% in Finnish Karelia. Adult birth cohorts showed that among those born in the 1940s, the sensitization to pollens and pets was at the same low level in both countries, but among younger generation born in the late 1970s, the difference was already manifold. Seropositivity to some pathogens, microbial content in house dust and drinking water seemed to confer allergy protection in Russia. In subsequent studies, it became apparent that on the Finnish side, healthy children had a more biodiverse living environment as well as greater diversity of certain bacterial classes on their skin than atopic children. Abundance of skin commensals, especially Acinetobacter (gammaproteobacteria), associated with anti-inflammatory gene expression in blood leucocytes. In vivo experiments with the mouse model demonstrated that intradermally applied Acinetobacter protected against atopic sensitization and lung inflammation. These observations support the notion that the epidemic of allergy and asthma results from reduced exposure to natural environments with rich microbiota, changed diet and sedentary lifestyle. Genetic studies have confirmed strong influence of lifestyle and environment. With our results from the Karelia study, a 10-year National Allergy Programme was started in 2008 to combat the epidemic in Finland.

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.

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.

Cytochrome P450 gene CYP337 and heritability of fitness traits in the Glanville fritillary butterfly.

Fitness-related life history traits often show substantial heritable genetic variation in natural populations, but knowledge of the genetic architecture of these traits is limited. In the Glanville fritillary butterfly, we measured the heritability of key life history traits in a large outdoor population cage during 2 years and generations and combined this experiment with an association study of a set of candidate genes. The genes were selected on the basis of previous genomic and transcriptomic studies and have been linked to the physiology and life history of this or other arthropod species. Heritability was high and significant for two traits, post-diapause larval development time (h(2) = 0.37) and lifetime egg (and larval) production (h(2) = 0.62); the latter is closely related to lifetime reproductive success and therefore fitness. We discovered a strong association between genetic polymorphism in the cytochrome P450 gene CYP337 and lifetime egg production, which accounted for 14% of the additive variance in egg production. This gene belongs to a group of cytochrome P450 genes that have a well-documented role in host plant adaptations in Lepidoptera and other insects and is likely to play an important role in the ecology and microevolution of the Glanville fritillary. This study provides a prime example of a gene associated with heritable fitness variation, measured under semi-natural ecological conditions.

Size differentiation in Finnish house sparrows follows Bergmann's rule with evidence of local adaptation.

Bergmann's rule predicts that individuals are larger in more poleward populations and that this size gradient has an adaptive basis. Hence, phenotypic divergence in size traits between populations (PST ) is expected to exceed the level of divergence by drift alone (FST ). We measured 16 skeletal traits, body mass and wing length in 409 male and 296 female house sparrows Passer domesticus sampled in 12 populations throughout Finland, where the species has its northernmost European distributional margin. Morphometric differentiation across populations (PST ) was compared with differentiation in 13 microsatellites (FST ). We find that twelve traits phenotypically diverged more than FST in both sexes, and an additional two traits diverged in males. The phenotypic divergence exceeded FST in several traits to such a degree that findings were robust also to strong between-population environmental effects. Divergence was particularly strong in dimensions of the bill, making it a strong candidate for the study of adaptive molecular genetic divergence. Divergent traits increased in size in more northern populations. We conclude that house sparrows show evidence of an adaptive latitudinal size gradient consistent with Bergmann's rule on the modest spatial scale of ca. 600 km.

Low genetic differentiation in a sedentary bird: house sparrow population genetics in a contiguous landscape.

The house sparrow Passer domesticus has been declining in abundance in many localities, including Finland. We studied the genetic diversity and differentiation of the house sparrow populations across Finland in the 1980s, at the onset of the species' decline in abundance. We genotyped 472 adult males (the less dispersive sex) from 13 locations in Finland (covering a range of 400 × 800 km) and one in Sweden (Stockholm) for 13 polymorphic microsatellite markers. Our analysis of Finnish ringing records showed that natal dispersal distances are limited (90% <16 km), which confirmed earlier finding from other countries. The Finnish populations were panmictic, and genetically very homogeneous and the limited dispersal was sufficiently large to maintain their connectivity. However, all Finnish populations differed significantly from the Stockholm population, even though direct geographical distance to it was often smaller than among Finnish populations. Hence, the open sea between Finland and Sweden appears to form a dispersal barrier for this species, whereas dispersal is much less constrained across the Finnish mainland (which lacks geographical barriers). Our findings provide a benchmark for conservation biologists and emphasize the influence of landscape structure on gene flow.

Heritability of and strong single gene (Pgi) effects on life-history traits in the Glanville fritillary butterfly.

We estimated broad-sense heritabilities (H(2)) of 13 female and seven male life-history traits of the Glanville fritillary butterfly (Melitaea cinxia) under semi-natural conditions in a large outdoor population cage. The analysis was based on full-sib families collected as young larvae in the field and reared under common garden conditions. We found significant genetic variance in female lifespan, fecundity, number of matings and host-plant preference as well as in male body mass and mobility. Apart from host-plant preference, female traits that were more strongly correlated with lifetime reproductive success (LRS; measured as total number of eggs laid) had higher H(2). LRS itself exhibited significant heritability. Host-plant preference had very high H(2), consistent with a previously reported genetically determined geographical cline in host-plant preference in the study area. Lifespan and egg hatching rate were significantly associated with a SNP in the coding region of the Pgi gene, for which there is previous evidence for balancing selection. Selection on Pgi, which furthermore shows spatial and temporal variation, may maintain genetic variance in fitness-related life-history traits. In contrast, we found no strong evidence for life-history trade-offs.

Spatiotemporal structure of host-pathogen interactions in a metapopulation.

The ecological and evolutionary dynamics of species are influenced by spatiotemporal variation in population size. Unfortunately, we are usually limited in our ability to investigate the numerical dynamics of natural populations across large spatial scales and over long periods of time. Here we combine mechanistic and statistical approaches to reconstruct continuous-time infection dynamics of an obligate fungal pathogen on the basis of discrete-time occurrence data. The pathogen, Podosphaera plantaginis, infects its host plant, Plantago lanceolata, in a metapopulation setting where the presence of the pathogen has been recorded annually for 6 years in approximately 4,000 host populations across an area of 50 km x 70 km in Finland. The dynamics are driven by strong seasonality, with a high extinction rate during winter and epidemic expansion in summer for local pathogen populations. We are able to identify with our model the regions in the study area where overwintering has been most successful. These overwintering sites represent foci that initiate local epidemics during the growing season. There is striking heterogeneity at the regional scale in both the overwintering success of the pathogen and the encounter intensity between the host and the pathogen. Such heterogeneity has profound implications for the coevolutionary dynamics of the interaction.

Fitness differences associated with Pgi SNP genotypes in the Glanville fritillary butterfly (Melitaea cinxia).

Allozyme variation at the phosphoglucose isomerase (PGI) locus in the Glanville fritillary butterfly (Melitaea cinxia) is associated with variation in flight metabolic rate, dispersal rate, fecundity and local population growth rate. To map allozyme to DNA variation and to survey putative functional variation in genomic DNA, we cloned the coding sequence of Pgi and identified nonsynonymous variable sites that determine the most common allozyme alleles. We show that these single-nucleotide polymorphisms (SNPs) exhibit significant excess of heterozygotes in field-collected population samples as well as in laboratory crosses. This is in contrast to previous results for the same species in which other allozymes and SNPs were in Hardy-Weinberg equilibrium or exhibited an excess of homozygotes. Our results suggest that viability selection favours Pgi heterozygotes. Although this is consistent with direct overdominance at Pgi, we cannot exclude the possibility that heterozygote advantage is caused by the presence of one or more deleterious alleles at linked loci.

Genetic spatial structure in a butterfly metapopulation correlates better with past than present demographic structure.

The Glanville fritillary butterfly (Melitaea cinxia) has been studied in the Aland Islands in Finland since 1991, where it occurs as a classic metapopulation in a large network of 4000 dry meadows. Much ecological work has been conducted on this species, but population genetic studies have been hampered by paucity of suitable genetic markers. Here, using single nucleotide polymorphisms and microsatellites developed for the Glanville fritillary, we examine the correspondence between the demographic and genetic spatial structures. Given the dynamic nature of the metapopulation, the current genetic spatial structure may bear a signal of past changes in population sizes and past patterns of gene flow rather than reflect the current demographic structure or landscape structure. We analyse this question with demographic data for 10 years, using the Rand index to assess the similarity between the genetic, demographic, and landscape spatial structures. Our results show that the current genetic spatial structure is better explained by the past rather than by the current demographic spatial structure or by the spatial configuration of the habitat in the landscape. Furthermore, current genetic diversity is significantly explained by past metapopulation sizes. The time lag between major demographic events and change in the genetic spatial structure and diversity has implications for the study of spatial dynamics.

Spatially realistic theory of metapopulation ecology.

Much of spatial ecology since the late 1960s has been dominated by two theories, the dynamic theory of island biogeography and the classical metapopulation theory. The latter theory largely replaced the former one in the 1980s, especially in conservation applications. It is only recently that ecologists have fully realized that a relatively simple general theory can be readily constructed that makes some of the simplifying assumptions of the two earlier theories unnecessary. The spatially realistic metapopulation theory thereby provides a more unified framework for spatial ecology than the island theory or the classical metapopulation theory. This article describes the application of the spatially realistic metapopulation theory to real metapopulations living in highly fragmented landscapes. I discuss the principal messages for population ecology and conservation biology, and I also place this theory into a broader context of other approaches to spatial ecology.

Patterns of 2-year population cycles in spatially extended host-parasitoid systems.

Coexisting but temporally separated cohorts of insects with a multiannual life cycle may have dissimilar average abundance, resulting in periodically fluctuating population density. In the case of the boreal moth genus Xestia with a 2-year life cycle and a distinct abundance difference between the two coexisting cohorts, empirical results and a simple model suggest that the oscillatory dynamics are maintained by interaction with a parasitoid wasp. Here we report theoretical results on a spatially extended version of the basic model and relate the modeling results to empirical observations. A spatially extended model may have domains oscillating in different phases as is the case between western and eastern Finnish Lapland. Spatial heterogeneity tends to fix the location of phase boundaries. In contrast, spatially homogeneous temporal fluctuations tend to synchronize populations in large regions.

Evolution of migration rate in a spatially realistic metapopulation model.

We use an individual-based, spatially realistic metapopulation model to study the evolution of migration rate. We first explore the consequences of habitat change in hypothetical patch networks on a regular lattice. If the primary consequence of habitat change is an increase in local extinction risk as a result of decreased local population sizes, migration rate increases. A nonmonotonic response, with migration rate decreasing at high extinction rate, was obtained only by assuming very frequent catastrophes. If the quality of the matrix habitat deteriorates, leading to increased mortality during migration, the evolutionary response is more complex. As long as habitat patch occupancy does not decrease markedly with increased migration mortality, reduced migration rate evolves. However, once mortality becomes so high that empty patches remain uncolonized for a long time, evolution tends to increase migration rate, which may lead to an "evolutionary rescue" in a fragmented landscape. Kin competition has a quantitative effect on the evolution of migration rate in our model, but these patterns in the evolution of migration rate appear to be primarily caused by spatiotemporal variation in fitness and mortality during migration. We apply the model to real habitat patch networks occupied by two checkerspot butterfly (Melitaea) species, for which sufficient data are available to estimate rigorously most of the model parameters. The model-predicted migration rate is not significantly different from the empirically observed one. Regional variation in patch areas and connectivities leads to regional variation in the optimal migration rate, predictions that can be tested empirically.

Extinction-colonization dynamics and host-plant choice in butterfly metapopulations.

Species living in highly fragmented landscapes often occur as metapopulations with frequent population turnover. Turnover rate is known to depend on ecological factors, such as population size and connectivity, but it may also be influenced by the phenotypic and genotypic composition of populations. The Glanville fritillary butterfly (Melitaea cinxia) in Finland uses two host-plant species that vary in their relative abundances among distinct habitat patches (dry meadows) in a large network of approximately 1,700 patches. We found no effect of host species use on local extinction. In contrast, population establishment was strongly influenced by the match between the host species composition of an empty habitat patch and the relative host use by larvae in previous years in the habitat patches that were well connected to the target patch. This "colonization effect" could be due to spatially variable plant acceptability or resistance or to spatially variable insect oviposition preference or larval performance. We show that spatial variation in adult oviposition preference occurs at the relevant spatial scale and that the other possible causes of the colonization effect can be discounted. We conclude that the colonization effect is generated by host preference influencing the movement patterns of ovipositing females. Migrant females with dissimilar host preferences have different perceptions of relative patch quality, which influences their likelihood of colonizing patches with particular host composition.

Spatially structured metapopulation models: global and local assessment of metapopulation capacity.

We model metapopulation dynamics in finite networks of discrete habitat patches with given areas and spatial locations. We define and analyze two simple and ecologically intuitive measures of the capacity of the habitat patch network to support a viable metapopulation. Metapopulation persistence capacity lambda(M) defines the threshold condition for long-term metapopulation persistence as lambda(M)>delta, where delta is defined by the extinction and colonization rate parameters of the focal species. Metapopulation invasion capacity lambda(I) sets the condition for successful invasion of an empty network from one small local population as lambda(I)>delta. The metapopulation capacities lambda(M) and lambda(I) are defined as the leading eigenvalue or a comparable quantity of an appropriate "landscape" matrix. Based on these definitions, we present a classification of a very general class of deterministic, continuous-time and discrete-time metapopulation models. Two specific models are analyzed in greater detail: a spatially realistic version of the continuous-time Levins model and the discrete-time incidence function model with propagule size-dependent colonization rate and a rescue effect. In both models we assume that the extinction rate increases with decreasing patch area and that the colonization rate increases with patch connectivity. In the spatially realistic Levins model, the two types of metapopulation capacities coincide, whereas the incidence function model possesses a strong Allee effect characterized by lambda(I)=0. For these two models, we show that the metapopulation capacities can be considered as simple sums of contributions from individual habitat patches, given by the elements of the leading eigenvector or comparable quantities. We may therefore assess the significance of particular habitat patches, including new patches that might be added to the network, for the metapopulation capacities of the network as a whole. We derive useful approximations for both the threshold conditions and the equilibrium states in the two models. The metapopulation capacities and the measures of the dynamic significance of particular patches can be calculated for real patch networks for applications in metapopulation ecology, landscape ecology, and conservation biology.

The metapopulation capacity of a fragmented landscape.

Ecologists and conservation biologists have used many measures of landscape structure to predict the population dynamic consequences of habitat loss and fragmentation, but these measures are not well justified by population dynamic theory. Here we introduce a new measure for highly fragmented landscapes, termed the metapopulation capacity, which is rigorously derived from metapopulation theory and can easily be applied to real networks of habitat fragments with known areas and connectivities. Technically, metapopulation capacity is the leading eigenvalue of an appropriate 'landscape' matrix. A species is predicted to persist in a landscape if the metapopulation capacity of that landscape is greater than a threshold value determined by the properties of the species. Therefore, metapopulation capacity can conveniently be used to rank different landscapes in terms of their capacity to support viable metapopulations. We present an empirical example on multiple networks occupied by an endangered species of butterfly. Using this theory, we may also calculate how the metapopulation capacity is changed by removing habitat fragments from or adding new ones into specific spatial locations, or by changing their areas. The metapopulation capacity should find many applications in metapopulation ecology, landscape ecology and conservation biology.