Heterozygosity-fitness correlations among wild populations of European tree frog (Hyla arborea) detect fixation load

Importance of Genetic-Fitness Correlations for the Conservation of Amphibians

Endangered animals suffer from isolation of their habitats. Isolation leads to a reduction in population size as well as a decrease in genetic diversity and a concomitant increase in the risk of extinction. Amphibians are the most endangered vertebrate class. Besides habitat loss, fragmentation and isolation, amphibians are threatened by emerging diseases e.g., chytrid fungus or Ranavirus. By employing experiments, researchers investigate whether changes in genetic diversity within or among isolated populations affect amphibian fitness. While genetic diversity estimates are based on molecular markers, typically microsatellites, fitness is mostly measured as tadpole performance in rearing experiments often under varying environmental conditions. Tadpole performances (e.g., body mass, growth rate and survival) have been found to be negatively affected by low genetic diversity, as several studies have found a positive association between genetic diversity and these fitness traits. Moreover, infection with pathogens also seems to be more likely in individuals or populations with lower genetic diversity. Overall, these genetic-fitness correlations seem to be more pronounced or detectable in smaller, declining populations but not in larger populations. Genomic studies, which sample a larger fraction of the genome, are still scarce in the conservation genetic literature on amphibians. These are likely to increase in upcoming years and may reveal adaptive variants that protect against dangerous pathogens or environmental changes. Altogether, genetic-fitness correlation studies should be a priority in order to develop effective management plans for the genetic rescue of isolated, imperilled amphibian populations.

Southern introgression increases adaptive immune gene variability in northern range margin populations of Fire-bellied toad

Northern range margin populations of the European fire-bellied toad (Bombina bombina) have rapidly declined during recent decades. Extensive agricultural land use has fragmented the landscape, leading to habitat disruption and loss, as well as eutrophication of ponds. In Northern Germany (Schleswig-Holstein) and Southern Sweden (Skåne), this population decline resulted in decreased gene flow from surrounding populations, low genetic diversity, and a putative reduction in adaptive potential, leaving populations vulnerable to future environmental and climatic changes. Previous studies using mitochondrial control region and nuclear transcriptome-wide SNP data detected introgressive hybridization in multiple northern B. bombina populations after unreported release of toads from Austria. Here, we determine the impact of this introgression by comparing the body conditions (proxy for fitness) of introgressed and nonintrogressed populations and the genetic consequences in two candidate genes for putative local adaptation (the MHC II gene as part of the adaptive immune system and the stress response gene HSP70 kDa). We detected regional differences in body condition and observed significantly elevated levels of within individual MHC allele counts in introgressed Swedish populations, associated with a tendency toward higher body weight, relative to regional nonintrogressed populations. These differences were not observed among introgressed and nonintrogressed German populations. Genetic diversity in both MHC and HSP was generally lower in northern than Austrian populations. Our study sheds light on the potential benefits of translocations of more distantly related conspecifics as a means to increase adaptive genetic variability and fitness of genetically depauperate range margin populations without distortion of local adaptation.

Evolutionary principles guiding amphibian conservation

The Anthropocene has witnessed catastrophic amphibian declines across the globe. A multitude of new, primarily human-induced drivers of decline may lead to extinction, but can also push species onto novel evolutionary trajectories. If these are recognized by amphibian biologists, they can be engaged in conservation actions. Here, we summarize how principles stemming from evolutionary concepts have been applied for conservation purposes, and address emerging ideas at the vanguard of amphibian conservation science. In particular, we examine the consequences of increased drift and inbreeding in small populations and their implications for practical conservation. We then review studies of connectivity between populations at the landscape level, which have emphasized the limiting influence of anthropogenic structures and degraded habitat on genetic cohesion. The rapid pace of environmental changes leads to the central question of whether amphibian populations can cope either by adapting to new conditions or by shifting their ranges. We gloomily conclude that extinction seems far more likely than adaptation or range shifts for most species. That said, conservation strategies employing evolutionary principles, such as selective breeding, introduction of adaptive variants through translocations, ecosystem interventions aimed at decreasing phenotype-environment mismatch, or genetic engineering, may effectively counter amphibian decline in some areas or for some species. The spread of invasive species and infectious diseases has often had disastrous consequences, but has also provided some premier examples of rapid evolution with conservation implications. Much can be done in terms of setting aside valuable amphibian habitat that should encompass both natural and agricultural areas, as well as designing protected areas to maximize the phylogenetic and functional diversity of the amphibian community. We conclude that an explicit consideration and application of evolutionary principles, although certainly not a silver bullet, should increase effectiveness of amphibian conservation in both the short and long term.

Genetically inherited tolerance may unveil trait dominance patterns in an amphibian model

Chemical contamination may cause genetic erosion in natural populations by wiping out the most sensitive genotypes. This is of upmost concern if the loss of genetic variability is irreversible due to contaminant-driven elimination of alleles, which may happen if tolerance is a recessive or incompletely dominant trait - the recessive tolerance inheritance (working-) hypothesis. Accordingly, this work investigated the tolerance inheritance to lethal levels of a metal-rich acid mine drainage (AMD) and to copper sulphate in a population of Pelophylax perezi. Time-to-death for each egg, after being exposed to 60% of a sample of acid mine drainage and to 9 mg/L Cu, was registered, and, for each egg mass, the median lethal time (LT₅₀) and respective quartiles (LT₂₅ and LT₇₅) were computed. Results suggested that genetically determined tolerance could be probably driven by incomplete dominance (with possible maternal effect influence), preliminarily supporting the initial hypothesis.

Local effects drive heterozygosity-fitness correlations in an outcrossing long-lived tree

Heterozygosity-fitness correlations (HFCs) have been used to understand the complex interactions between inbreeding, genetic diversity and evolution. Although frequently reported for decades, evidence for HFCs was often based on underpowered studies or inappropriate methods, and hence their underlying mechanisms are still under debate. Here, we used 6100 genome-wide single nucleotide polymorphisms (SNPs) to test for general and local effect HFCs in maritime pine (Pinus pinaster Ait.), an iconic Mediterranean forest tree. Survival was used as a fitness proxy, and HFCs were assessed at a four-site common garden under contrasting environmental conditions (total of 16 288 trees). We found no significant correlations between genome-wide heterozygosity and fitness at any location, despite variation in inbreeding explaining a substantial proportion of the total variance for survival. However, four SNPs (including two non-synonymous mutations) were involved in significant associations with survival, in particular in the common gardens with higher environmental stress, as shown by a novel heterozygosity-fitness association test at the species-wide level. Fitness effects of SNPs involved in significant HFCs were stable across maritime pine gene pools naturally growing in distinct environments. These results led us to dismiss the general effect hypothesis and suggested a significant role of heterozygosity in specific candidate genes for increasing fitness in maritime pine. Our study highlights the importance of considering the species evolutionary and demographic history and different spatial scales and testing environments when assessing and interpreting HFCs.

Predictable allele frequency changes due to habitat fragmentation in the Glanville fritillary butterfly

Describing the evolutionary dynamics of now extinct populations is challenging, as their genetic composition before extinction is generally unknown. The Glanville fritillary butterfly has a large extant metapopulation in the Åland Islands in Finland, but declined to extinction in the nearby fragmented southwestern (SW) Finnish archipelago in the 20th century. We genotyped museum samples for 222 SNPs across the genome, including SNPs from candidate genes and neutral regions. SW Finnish populations had significantly reduced genetic diversity before extinction, and their allele frequencies gradually diverged from those in contemporary Åland populations over 80 y. We identified 15 outlier loci among candidate SNPs, mostly related to flight, in which allele frequencies have changed more than the neutral expectation. At outlier loci, allele frequencies in SW Finland shifted in the same direction as newly established populations deviated from old local populations in contemporary Åland. Moreover, outlier allele frequencies in SW Finland resemble those in fragmented landscapes as opposed to continuous landscapes in the Baltic region. These results indicate selection for genotypes associated with good colonization capacity in the highly fragmented landscape before the extinction of the populations. Evolutionary response to habitat fragmentation may have enhanced the viability of the populations, but it did not save the species from regional extinction in the face of severe habitat loss and fragmentation. These results highlight a potentially common situation in changing environments: evolutionary changes are not strong enough to fully compensate for the direct adverse effects of environmental change and thereby rescue populations from extinction.

Microevolution due to pollution in amphibians: A review on the genetic erosion hypothesis

The loss of genetic diversity, due to exposure to chemical contamination (genetic erosion), is a major threat to population viability. Genetic erosion is the loss of genetic variation: the loss of alleles determining the value of a specific trait or set of traits. Almost a third of the known amphibian species is considered to be endangered and a decrease of genetic variability can push them to the verge of extinction. This review indicates that loss of genetic variation due to chemical contamination has effects on: 1) fitness, 2) environmental plasticity, 3) co-tolerance mechanisms, 4) trade-off mechanisms, and 5) tolerance to pathogens in amphibian populations.

Effect of biogeographic history on population vulnerability in European amphibians

The genetic diversity of populations, which contributes greatly to their adaptive potential, is negatively affected by anthropogenic habitat fragmentation and destruction. However, continental-scale losses of genetic diversity also resulted from the population expansions that followed the end of the last glaciation, an element that is rarely considered in a conservation context. We addressed this issue in a meta-analysis in which we compared the spatial patterns of vulnerability of 18 widespread European amphibians in light of phylogeographic histories (glacial refugia and postglacial routes) and anthropogenic disturbances. Conservation statuses significantly worsened with distances from refugia, particularly in the context of industrial agriculture; human population density also had a negative effect. These findings suggest that features associated with the loss of genetic diversity in post-glacial amphibian populations (such as enhanced fixation load or depressed adaptive potential) may increase their susceptibility to current threats (e.g., habitat fragmentation and pesticide use). We propose that the phylogeographic status of populations (i.e., refugial vs. post-glacial) should be considered in conservation assessments for regional and national red lists.

First-generation linkage map for the European tree frog (Hyla arborea) with utility in congeneric species

Background

Western Palearctic tree frogs (Hyla arborea group) represent a strong potential for evolutionary and conservation genetic research, so far underexploited due to limited molecular resources. New microsatellite markers have recently been developed for Hyla arborea, with high cross-species utility across the entire circum-Mediterranean radiation. Here we conduct sibship analyses to map available markers for use in future population genetic applications.

Findings

We characterized eight linkage groups, including one sex-linked, all showing drastically reduced recombination in males compared to females, as previously documented in this species. Mapping of the new 15 markers to the ~200 My diverged Xenopus tropicalis genome suggests a generally conserved synteny with only one confirmed major chromosome rearrangement.

Conclusions

The new microsatellites are representative of several chromosomes of H. arborea that are likely to be conserved across closely-related species. Our linkage map provides an important resource for genetic research in European Hylids, notably for studies of speciation, genome evolution and conservation.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-850) contains supplementary material, which is available to authorized users.

Stronger transferability but lower variability in transcriptomic- than in anonymous microsatellites: evidence from Hylid frogs

A simple way to quickly optimize microsatellites in nonmodel organisms is to reuse loci available in closely related taxa; however, this approach can be limited by the stochastic and low cross-amplification success experienced in some groups (e.g. amphibians). An efficient alternative is to develop loci from transcriptome sequences. Transcriptomic microsatellites have been found to vary in their levels of cross-species amplification and variability, but this has to date never been tested in amphibians. Here, we compare the patterns of cross-amplification and levels of polymorphism of 18 published anonymous microsatellites isolated from genomic DNA vs. 17 loci derived from a transcriptome, across nine species of tree frogs (Hyla arborea and Hyla cinerea group). We established a clear negative relationship between divergence time and amplification success, which was much steeper for anonymous than transcriptomic markers, with half-lives (time at which 50% of the markers still amplify) of 1.1 and 37 My, respectively. Transcriptomic markers are significantly less polymorphic than anonymous loci, but remain variable across diverged taxa. We conclude that the exploitation of amphibian transcriptomes for developing microsatellites seems an optimal approach for multispecies surveys (e.g. analyses of hybrid zones, comparative linkage mapping), whereas anonymous microsatellites may be more informative for fine-scale analyses of intraspecific variation. Moreover, our results confirm the pattern that microsatellite cross-amplification is greatly variable among amphibians and should be assessed independently within target lineages. Finally, we provide a bank of microsatellites for Palaearctic tree frogs (so far only available for H. arborea), which will be useful for conservation and evolutionary studies in this radiation.

Conservation phylogeography: does historical diversity contribute to regional vulnerability in European tree frogs (Hyla arborea)?

Documenting and preserving the genetic diversity of populations, which conditions their long-term survival, have become a major issue in conservation biology. The loss of diversity often documented in declining populations is usually assumed to result from human disturbances; however, historical biogeographic events, otherwise known to strongly impact diversity, are rarely considered in this context. We apply a multilocus phylogeographic study to investigate the late-Quaternary history of a tree frog (Hyla arborea) with declining populations in the northern and western part of its distribution range. Mitochondrial and nuclear polymorphisms reveal high genetic diversity in the Balkan Peninsula, with a spatial structure moulded by the last glaciations. While two of the main refugial lineages remained limited to the Balkans (Adriatic coast, southern Balkans), a third one expanded to recolonize Northern and Western Europe, loosing much of its diversity in the process. Our findings show that mobile and a priori homogeneous taxa may also display substructure within glacial refugia ('refugia within refugia') and emphasize the importance of the Balkans as a major European biodiversity centre. Moreover, the distribution of diversity roughly coincides with regional conservation situations, consistent with the idea that historically impoverished genetic diversity may interact with anthropogenic disturbances, and increase the vulnerability of populations. Phylogeographic models seem important to fully appreciate the risks of local declines and inform conservation strategies.

Within- and among-population impact of genetic erosion on adult fitness-related traits in the European tree frog Hyla arborea

Assessing in wild populations how fitness is impacted by inbreeding and genetic drift is a major goal for conservation biology. An approach to measure the detrimental effects of inbreeding on fitness is to estimate correlations between molecular variation and phenotypic performances within and among populations. Our study investigated the effect of individual multilocus heterozygosity on body size, body condition and reproductive investment of males (that is, chorus attendance) and females (that is, clutch mass and egg size) in both small fragmented and large non-fragmented populations of European tree frog (Hyla arborea). Because adult size and/or condition and reproductive investment are usually related, genetic erosion may have detrimental effects directly on reproductive investment, and also on individual body size and condition that in turn may affect reproductive investment. We confirmed that the reproductive investment was highly size-dependent for both sexes. Larger females invested more in offspring production, and larger males attended the chorus in the pond more often. Our results did not provide evidence for a decline in body size, condition and reproductive effort with decreased multilocus heterozygosity both within and among populations. We showed that the lack of heterozygosity-fitness correlations within populations probably resulted from low inbreeding levels (inferior to ca. 20% full-sib mating rate), even in the small fragmented populations. The detrimental effects of fixation load were either low in adults or hidden by environmental variation among populations. These findings will be useful to design specific management actions to improve population persistence.

Heterozygosity-fitness correlations and inbreeding depression in two critically endangered mammals

The relation among inbreeding, heterozygosity, and fitness has been studied primarily among outbred populations, and little is known about these phenomena in endangered populations. Most researchers conclude that the relation between coefficient of inbreeding estimated from pedigrees and fitness traits (inbreeding-fitness correlations) better reflects inbreeding depression than the relation between marker heterozygosity and fitness traits (heterozygosity-fitness correlations). However, it has been suggested recently that heterozygosity-fitness correlations should only be expected when inbreeding generates extensive identity disequilibrium (correlations in heterozygosity and homozygosity across loci throughout the genome). We tested this hypothesis in Mohor gazelle (Gazella dama mhorr) and Iberian lynx (Lynx pardinus). For Mohor gazelle, we calculated the inbreeding coefficient and measured heterozygosity at 17 microsatellite loci. For Iberian lynx, we measured heterozygosity at 36 microsatellite loci. In both species we estimated semen quality, a phenotypic trait directly related to fitness that is controlled by many loci and is affected by inbreeding depression. Both species showed evidence of extensive identity disequilibrium, and in both species heterozygosity was associated with semen quality. In the Iberian lynx the low proportion of normal sperm associated with low levels of heterozygosity was so extreme that it is likely to limit the fertility of males. In Mohor gazelle, although heterozygosity was associated with semen quality, inbreeding coefficient was not. This result suggests that when coefficient of inbreeding is calculated on the basis of a genealogy that begins after a long history of inbreeding, the coefficient of inbreeding fails to capture previous demographic information because it is a poor estimator of accumulated individual inbreeding. We conclude that among highly endangered species with extensive identity disequilibrium, examination of heterozygosity-fitness correlations may be an effective way to detect inbreeding depression, whereas inbreeding-fitness correlations may be poor indicators of inbreeding depression if the pedigree does not accurately reflect the history of inbreeding.

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.

Variability of individual genetic load: consequences for the detection of inbreeding depression

Inbreeding depression is a key factor affecting the persistence of natural populations, particularly when they are fragmented. In species with mixed mating systems, inbreeding depression can be estimated at the population level by regressing the average progeny fitness by the selfing rate of their mothers. We applied this method using simulated populations to investigate how population genetic parameters can affect the detection power of inbreeding depression. We simulated individual selfing rates and genetic loads from which we computed fitness values. The regression method yielded high statistical power, inbreeding depression being detected as significant (5 % level) in 92 % of the simulations. High individual variation in selfing rate and high mean genetic load led to better detection of inbreeding depression while high among-individual variation in genetic load made it more difficult to detect inbreeding depression. For a constant sampling effort, increasing the number of progenies while decreasing the number of individuals per progeny enhanced the detection power of inbreeding depression. We discuss the implication of among-mother variability of genetic load and selfing rate on inbreeding depression studies.

Genetic erosion in wild populations makes resistance to a pathogen more costly

Populations that have suffered from genetic erosion are expected to exhibit reduced average trait values or decreased variation in adaptive traits when experiencing periodic or emergent stressors such as infectious disease. Genetic erosion may consequentially modify the ability of a potential host population to cope with infectious disease emergence. We experimentally investigate this relationship between genetic variability and host response to exposure to an infectious agent both in terms of susceptibility to infection and indirect parasite-mediated responses that also impact fitness. We hypothesized that the deleterious consequences of exposure to the pathogen (Batrachochytrium dendrobatidis) would be more severe for tadpoles descended from European treefrog (Hyla arborea) populations lacking genetic variability. Although all exposed tadpoles lacked detectable infection, we detected this relationship for some indirect host responses, predominantly in genetically depleted animals, as well as an interaction between genetic variability and pathogen dose on life span during the postmetamorphic period. Lack of infection and a decreased mass and postmetamorphic life span in low genetic diversity tadpoles lead us to conclude that genetic erosion, while not affecting the ability to mount effective resistance strategies, also erodes the capacity to invest in resistance, increased tadpole growth rate, and metamorphosis relatively simultaneously.