Strong isolation by distance among local populations of an endangered butterfly species (Euphydryas aurinia)

The marsh fritillary (Euphydryas aurinia) is a critically endangered butterfly species in Denmark known to be particularly vulnerable to habitat fragmentation due to its poor dispersal capacity. We identified and genotyped 318 novel SNP loci across 273 individuals obtained from 10 small and fragmented populations in Denmark using a genotyping‐by‐sequencing (GBS) approach to investigate its population genetic structure. Our results showed clear genetic substructuring and highly significant population differentiation based on genetic divergence (F _ST) among the 10 populations. The populations clustered in three overall clusters, and due to further substructuring among these, it was possible to clearly distinguish six clusters in total. We found highly significant deviations from Hardy–Weinberg equilibrium due to heterozygote deficiency within every population investigated, which indicates substructuring and/or inbreeding (due to mating among closely related individuals). The stringent filtering procedure that we have applied to our genotype quality could have overestimated the heterozygote deficiency and the degree of substructuring of our clusters but is allowing relative comparisons of the genetic parameters among clusters. Genetic divergence increased significantly with geographic distance, suggesting limited gene flow at spatial scales comparable to the dispersal distance of individual butterflies and strong isolation by distance. Altogether, our results clearly indicate that the marsh fritillary populations are genetically isolated. Further, our results highlight that the relevant spatial scale for conservation of rare, low mobile species may be smaller than previously anticipated.     

Experimental test of genetic rescue in isolated populations of brook trout

Genetic rescue is an increasingly considered conservation measure to address genetic erosion associated with habitat loss and fragmentation. The resulting gene flow from facilitating migration may improve fitness and adaptive potential, but is not without risks (e.g., outbreeding depression). Here, we conducted a test of genetic rescue by translocating ten (five of each sex) brook trout (Salvelinus fontinalis) from a single source to four nearby and isolated stream populations. To control for the demographic contribution of translocated individuals, ten resident individuals (five of each sex) were removed from each recipient population. Prior to the introduction of translocated individuals, the two smallest above‐barrier populations had substantially lower genetic diversity, and all populations had reduced effective number of breeders relative to adjacent below‐barrier populations. In the first reproductive bout following translocation, 31 of 40 (78%) translocated individuals reproduced successfully. Translocated individuals contributed to more families than expected under random mating and generally produced larger full‐sibling families. We observed relatively high (>20%) introgression in three of the four recipient populations. The translocations increased genetic diversity of recipient populations by 45% in allelic richness and 25% in expected heterozygosity. Additionally, strong evidence of hybrid vigour was observed through significantly larger body sizes of hybrid offspring relative to resident offspring in all recipient populations. Continued monitoring of these populations will test for negative fitness effects beyond the first generation. However, these results provide much‐needed experimental data to inform the potential effectiveness of genetic rescue‐motivated translocations.     

Dispersal and gene flow in the rare, parasitic Large Blue butterfly Maculinea arion

Dispersal is crucial for gene flow and often determines the long‐term stability of meta‐populations, particularly in rare species with specialized life cycles. Such species are often foci of conservation efforts because they suffer disproportionally from degradation and fragmentation of their habitat. However, detailed knowledge of effective gene flow through dispersal is often missing, so that conservation strategies have to be based on mark–recapture observations that are suspected to be poor predictors of long‐distance dispersal. These constraints have been especially severe in the study of butterfly populations, where microsatellite markers have been difficult to develop. We used eight microsatellite markers to analyse genetic population structure of the Large Blue butterfly Maculinea arion in Sweden. During recent decades, this species has become an icon of insect conservation after massive decline throughout Europe and extinction in Britain followed by reintroduction of a seed population from the Swedish island of Öland. We find that populations are highly structured genetically, but that gene flow occurs over distances 15 times longer than the maximum distance recorded from mark–recapture studies, which can only be explained by maximum dispersal distances at least twice as large as previously accepted. However, we also find evidence that gaps between sites with suitable habitat exceeding ～20 km induce genetic erosion that can be detected from bottleneck analyses. Although further work is needed, our results suggest that M. arion can maintain fully functional metapopulations when they consist of optimal habitat patches that are no further apart than ～10 km.     

Microsatellite markers indicate genetic differences between cultivated and natural populations of endangered Taxus yunnanensis

Taxus yunnanensis (Taxaceae) is an endangered plant known for its effective natural anti‐cancer metabolite, taxol. To assess the effectiveness of T. yunnanensis ex situ conservation in China, the genetic diversity and genetic structure of nine cultivated and 14 natural populations were compared using 11 polymorphic microsatellite markers. The purpose of this comparison was to determine whether the cultivated populations were genetic representatives of natural types. Among the 14 natural populations studied, three were from Tibet and 11 were from the Hengduan Mountains. There were fewer genetic variations among cultivated populations than among natural populations, but this difference was not significant. However, 23 unique alleles were detected in cultivated populations and 14 were detected in natural populations. Striking differences in genetic structure (G_ST = 0.365) was observed between them. STRUCTURE, AMOVA and cluster analyses also consistently showed a clear genetic separation between the 14 natural and seven cultivated populations. Only five natural populations from the south‐central Hengduan Mountains were genetically represented in two cultivated collections. The reason for the occurrence of distinct genetic differences between cultivated and natural populations and the loss of natural population‐specific alleles in the cultivated populations can be attributed to uneven selection during establishment. These results suggest that the genetic structure of T. yunnanensis may not have been conserved effectively. Based on our study, practical conservation strategies for ex situ preservation of T. yunnanensis have been developed. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 177 , 450–461.     

Low genetic variation reduces cross-compatibility and offspring fitness in populations of a narrow endemic plant with a self-incompatibility system

Genetic variation was shown earlier to bereduced in smaller populations of the narrowendemic putatively self-incompatible Cochlearia bavarica. To test whether thisnegatively affects plant fitness by reducedavailability of compatible mates and byinbreeding depression, we studied effects ofpopulation size and pollination treatments oncross-compatibility and offspring fitness in 16isolated populations of this plant. After openpollination, compatibility of crosses (i.e.,whether at least one fruit developed per markedflower), fruit set of compatible crosses, andcumulative fitness (number of plants permaternal ovule) after 14 months in a commongarden were lower for plants from smallerpopulations. Throughout the study, cumulativefitness was lower after hand pollination withpollen of one donor than after open pollination(finally 73.4% lower), suggesting that severalpollen donors or single pollen donors of higherquality are involved in open pollination.Moreover, cumulative fitness was lower afterhand selfing than after hand outcrossing(finally 69.4% lower), indicating bothinbreeding depression and reduced compatibilityafter selfing. High self-compatibility(40.6%), dry stigmas, and differences in thecompatibility of 11 of 33 experimentalreciprocal crosses between plant pairsconfirmed that C. bavarica has asporophytic self-incompatibility system, as iscommon in the Brassicaceae. Our studydemonstrates, that plants in smallerpopulations of species with a sporophyticself-incompatibility system can experiencetwofold fitness reductions associated withreduced genetic variability, i.e., twofoldgenetic Allee effects: via reducedcross-compatibility and via reduced offspringfitness.     

Rapid loss of genetic variation in an endangered possum

The endangered mountain pygmy possum is the only Australian marsupial that hibernates under snow cover. Most of its alpine habitat was burnt by a rare fire in 2003, and habitat loss and disturbance have also occurred owing to ski resort development. Here we show that there has been a rapid loss of genetic variation following habitat loss associated with resort development, but no detectable loss of alleles or decrease in heterozygosity following the fire.     

Dynamic landscapes shape post‐wildfire recolonisation and genetic structure of the endangered Hermes copper (Lycaena hermes) butterfly

1. The coastal sage scrub vegetation community experiences frequent fires, so the long‐term survival of species depends on the rate of recolonisations exceeding the rate of local extinctions. Recolonisation of these post‐wildfire habitats probably requires long‐distance dispersal events. These movements can also counter detrimental impacts associated with inbreeding. 2. The Hermes copper (Lycaena hermes) is an extremely rare butterfly inhabiting coastal sage scrub adjacent to San Diego, California, USA. Habitat loss due to urbanisation and impacts of recent wildfires has greatly restricted its range, prompting the United States Fish and Wildlife Service to list the Hermes copper as a Candidate Species in 2011. 3. Surveys for Hermes copper butterflies in 2010–2013 documented only two recolonisation events following wildfires in 2003 and 2007. Larger populations were readily detected each year, but detection of smaller populations was inconsistent. 4. Amplified fragment length polymorphism was used to identify potential genetic discontinuities within this species across the landscape. Results indicated that movements across much of the landscape were possible historically. However, individuals from three peripheral populations exhibited a higher degree of differentiation, suggesting more restricted dispersal in these regions. 5. From the results, it can be concluded that historically Hermes copper butterflies were able to move among habitat patches prior to recent changes in the landscape. However, low post‐fire recolonisation rates suggest limited contemporary dispersal, probably due to recent habitat fragmentation. This fragmentation is a relatively new event, as the human population in San Diego County experienced substantial growth in the late 20th Century.     

Post-bottleneck genetic diversity of elephant populations in South Africa, revealed using microsatellite analysis

Widespread hunting had fragmented and severely reduced elephant populations in South Africa by 1900. Elephant numbers increased during the 1900s, although rates of recovery of individual populations varied. The Kruger National Park elephant population increased rapidly, to more than 6000 by 1967, with recruitment boosted by immigration from Mozambique. The Addo Elephant National Park population was reduced to 11 elephants in 1931 and remains relatively small (n = 325). Loss of genetic variation is expected to occur whenever a population goes through a bottleneck, especially when post-bottleneck recovery is slow. Variation at nine polymorphic microsatellite loci was analysed for Kruger and Addo elephants, as well as museum specimens of Addo elephants shot prior to the population bottleneck. Significantly reduced genetic variation and heterozygosity were observed in Addo in comparison to Kruger (mean alleles/locus and H(E): Addo 1.89, 0.18; Kruger 3.89, 0.44). Two alleles not present in the current Addo population were observed in the museum specimens. Addo elephants represent a genetic subset of the Kruger population, with high levels of genetic differentiation resulting from rapid genetic drift. The Kruger population is low in genetic diversity in comparison to East African elephants, confirming this population also suffered an appreciable bottleneck.     

Survival of the endangered butterfly Lycaena helle in a fragmented environment: Genetic analyses over 15 years

Temporal changes in allele frequencies are often assumed in studies addressing the history of populations affected by different anthropogenic and natural impacts at different time scales. Yet, few studies directly compare the genetic composition of populations over time spans of more than 10 years. Therefore, to test the genetic effects of 15 years of population isolation in the butterfly Lycaena helle, we analysed 472 individuals from 27 samples, of which nine were collected in 1991 and 18 in 2006. Sampling was performed in five mountain regions (Pyrenees, Massif Central, Jura, Vosges and Ardennes). Genetic analyses were performed using five polymorphic microsatellites. Old and new samples of identical or neighbouring populations revealed similar genetic differentiations among these five mountain regions. A comparatively strong genetic differentiation among populations combined with a high amount of private alleles for each mountain area was detected, but mountain‐specific alleles were in most cases identical in 1991 and 2006. Nevertheless, the obtained data also indicate moderate changes between 1991 and 2006 in the species’ genetic structure – genetic differentiation among local populations increased marginally and allele frequencies showed corresponding modifications. A significant decline in genetic diversity was not detectable, and nine private alleles exclusive to a single mountain region were only detected in samples from the year 1991, whereas eleven were only observed in the individuals collected in 2006. These observations might indicate the results of genetic drift within isolated populations.     

Microsatellite analysis of population structure and genetic differentiation within and between populations of the root vole, Microtus oeconomus in the Netherlands

Eight microsatellite markers for the root vole (Microtus oeconomus) were developed to assess the amount of genetic variation for nine Dutch root vole populations from four different regions, and to evaluate the degree of differentiation and isolation. All eight microsatellite loci were found to be highly variable with observed heterozygosity values ranging from 0.61 to 0.82. These values are similar to those observed for more distant populations from Norway, Finland and Germany. Therefore, the populations seem not particularly depauperate of genetic variation at the microsatellite level. Genetically, the Dutch populations were found to have diverged considerably. Pairwise comparisons of all populations studied revealed FST values significantly greater than zero for most comparisons. However, the magnitude of these values considerably depends on the compared population pair. The level of differentiation between local populations within Dutch regions is generally significantly lower than the differentiation between Dutch regions. The level of differentiation between Dutch regions, however, is not significantly different from that between populations of larger geographical distance. This implies that the regional Dutch populations are both isolated from each other and from other European populations. The observation that even local populations show low but significant genetic differentiation may be indicative for progressive isolation of these populations.     

Microsatellite analysis of genetic variation among and within Alpine marmot populations in the French Alps

The genetic structure of the Alpine marmot, Marmota marmota, was studied by an analysis of five polymorphic microsatellite loci. Eight locations were sampled in the French Alps, one from Les Ecrins valley (n = 160), another from La Sassière valley (n = 289) and the six others from the Maurienne valley (n = 139). Information on social group structure was available for both Les Ecrins and La Sassière but not for the other samples. The high levels of genetic diversity observed are at odds with the results obtained using microsatellites, minisatellites and allozymes on Alpine marmots from Germany, Austria and Switzerland. Strong deficits in heterozygotes were found in Les Ecrins and La Sassière. They are caused by a Wahlund effect due to the family structure (i.e. differentiation between the family groups). The family groups exhibit excess of heterozygotes rather than deficits. This may be caused by outbreeding and this is compatible with recent results from the genetics of related social species when information on the social structure is taken into account. The observed outbreeding could be the result of females mating with transient males or males coming from neighbouring colonies. Both indicate that the species may not be as monogamous as is usually believed. The results are also compatible with a male-biased dispersal but do not allow us to exclude some female migration. We also found a significant correlation between geographical and genetic distance indicating that isolation by distance could be an issue in marmots. This study is the first that analysed populations of marmots taking into account the social structure within populations and assessing inbreeding at different levels (region, valley, population, and family groups). Our study clearly demonstrated that the sampling strategy and behavioural information can have dramatic effects on both the results and interpretation of the genetic data.     

Insights into recently fragmented vole populations from combined genetic and demographic data

We combined demographic and genetic data to evaluate the effects of habitat fragmentation on the population structure of the California red-backed vole (Clethrionomys californicus). We analysed variation in the mitochondrial DNA (mtDNA) control region and five nuclear microsatellite loci in small samples collected from two forest fragments and an unfragmented control site in 1990-91. We intensively sampled the same forest fragments and two different control sites in 1998 and 1999. Vole abundances fluctuated greatly at sizes below 50 individuals per fragment. Fragment populations had significantly lower mtDNA allelic diversity than controls, but not nuclear heterozygosity or numbers of alleles. The use of only trapping and/or mtDNA marker data would imply that fragment populations are at least partially isolated and vulnerable to inbreeding depression. In contrast, the abundance estimates combined with microsatellite data show that small fragment populations must be linked to nearby forests by high rates of migration. These results provide evidence for the usefulness of combining genetic and demographic data to understand nonequilibrium population structure in recently fragmented habitats.     

Microsatellite variation and population structure of a recovering Tree frog (Hyla arborea L.) metapopulation

Numbers and sizes of populations of the European tree frog in The Netherlands have dramatically decreased in the second half of the last century due to extensive habitat destruction and fragmentation. We have studied the genetic structure of a slowly recovering meta-population. Strong genetic differentiation, estimated at eight microsatellite loci, was found between clusters of populations (F _st-values above 0.2). Within clusters, consisting of ponds within a radius of about 5 km, European tree frog populations were less differentiated (F _st<0.08) and exact tests showed that most of the ponds within clusters were not significantly differentiated. Although local population sizes have been increasing since 1989, and some new ponds have been colonised in the direct vicinity of ponds that have been populated continuously, little evidence for gene flow between clusters of ponds was found (only one exception). Furthermore, levels of genetic diversity were low compared to populations in comparable areas elsewhere in Europe. Therefore, a continuous conservation effort is needed to prevent any further loss of genetic diversity. The alleviation of important barriers to dispersal between the clusters should be given a high priority for the restoration of the meta-population as a whole.     

Fine-scale spatial genetic structure and gene flow in a small, fragmented population of Sinojackia rehderiana (Styracaceae), an endangered tree species endemic to China

Populations of Sinojackia rehderiana are highly threatened and have small and scattered distribution due to habitat fragmentation and human activities. Understanding changes in genetic diversity, the fine-scale spatial genetic structure (SGS) at different life stages and gene flow of S. rehderiana is critical for developing successful conservation strategies for fragmented populations of this endangered species. In this study, 208 adults, 114 juveniles and 136 seedlings in a 50 × 100-m transect within an old-growth forest were mapped and genotyped using eight microsatellite makers to investigate the genetic diversity and SGS of this species. No significant differences in genetic diversity among different life-history stages were found. However, a significant heterozygote deficiency in adults and seedlings may result from substantial biparental inbreeding. Significant fine-scale spatial structure was found in different life-history stages within 19 m, suggesting that seed dispersal mainly occurred near a mother tree. Both historical and contemporary estimates of gene flow (13.06 and 16.77 m) indicated short-distance gene dispersal in isolated populations of S. rehderiana. The consistent spatial structure revealed in different life stages is most likely the result of limited gene flow. Our results have important implications for conservation of extant populations of S. rehderiana. Measures for promoting pollen flow should be taken for in situ conservation. The presence of a SGS in fragmented populations implies that seeds for ex situ conservation should be collected from trees at least 19-m apart to reduce genetic similarity between neighbouring individuals.     

Genetic differentiation within metapopulations of Euphydryas aurinia and Melitaea phoebe in China

We analyzed genetic differentiation within metapopulations of two species of checkerspot butterfly, Euphydryas aurinia and Melitaea phoebe, in China. To generate genetic information, we used a new molecular technique, DALP – direct amplified length polymorphism. AMOVA results showed that most of the variation occurred among individuals within local populations of both E. aurinia and M. phoebe. However, while there was differentiation among local population in E. aurinia (P < 0.001), there was no subdivision in metapopulation of M. phoebe (P = 0.210). This is consistent with the behavior of M. phoebe adults being more dispersive than E. aurinia. Within the M. phoebe metapopulation, three neighboring patches were always occupied during the observation period (1998–2000). In addition, the number of individuals in these three populations accounted for the majority of M. phoebe larvae, and hence we conclude that the M. phoebe metapopulation might exist as a source-sink metapopulation. On the other hand, the E. aurinia metapopulation is an example of a classical metapopulation. Therefore, the conservation management of these two species should reflect these differences.     

Fragmented habitat compensates for the adverse effects of genetic bottleneck

1. Download : Download high-res image (256KB)
2. Download : Download full-size image

     

Demographic fluctuations lead to rapid and cyclic shifts in genetic structure among populations of an alpine butterfly,Parnassius smintheus

Many populations, especially in insects, fluctuate in size, and periods of particularly low population size can have strong effects on genetic variation. Effects of demographic bottlenecks on genetic diversity of single populations are widely documented. Effects of bottlenecks on genetic structure among multiple interconnected populations are less studied, as are genetic changes across multiple cycles of demographic collapse and recovery. We take advantage of a long‐term data set comprising demographic, genetic and movement data from a network of populations of the butterfly, Parnassius smintheus, to examine the effects of fluctuating population size on spatial genetic structure. We build on a previous study that documented increased genetic differentiation and loss of spatial genetic patterns (isolation by distance and by intervening forest cover) after a network‐wide bottleneck event. Here, we show that genetic differentiation was reduced again and spatial patterns returned to the system extremely rapidly, within three years (i.e. generations). We also show that a second bottleneck had similar effects to the first, increasing differentiation and erasing spatial patterns. Thus, bottlenecks consistently drive random divergence of allele frequencies among populations in this system, but these effects are rapidly countered by gene flow during demographic recovery. Our results reveal a system in which the relative influence of genetic drift and gene flow continually shift as populations fluctuate in size, leading to cyclic changes in genetic structure. Our results also suggest caution in the interpretation of patterns of spatial genetic structure, and its association with landscape variables, when measured at only a single point in time.     

Microsatellite variation and structure of 28 populations of the common wetland plant, Lychnis flos-cuculi L., in a fragmented landscape

Habitat fragmentation is known to cause genetic differentiation between small populations of rare species and decrease genetic variation within such populations. However, common species with recently fragmented populations have rarely been studied in this context. We investigated genetic variation and its relationship to population size and geographical isolation of populations of the common plant species, Lychnis flos-cuculi L., in fragmented fen grasslands. We analysed 467 plants from 28 L. flos-cuculi populations of different sizes (60 000-54 000 flowering individuals) in northeastern Switzerland using seven polymorphic microsatellite loci. Genetic differentiation between populations is small (F(ST) = 0.022; amova; P < 0.001), suggesting that gene flow among populations is still high or that habitat fragmentation is too recent to result in pronounced differentiation. Observed heterozygosity (H(O) = 0.44) significantly deviates from Hardy-Weinberg equilibrium, and within-population inbreeding coefficient F(IS) is high (0.30-0.59), indicating a mixed mating breeding system with substantial inbreeding in L. flos-cuculi. Gene diversity is the only measure of genetic variation which decreased with decreasing population size (R = 0.42; P < 0.05). While our results do not indicate pronounced effects of habitat fragmentation on genetic variation in the still common L. flos-cuculi, the lower gene diversity of smaller populations suggests that the species is not entirely unaffected.     

Differentiation and levels of genetic variation in northern European lynx (Lynx lynx) populations revealed by microsatellites and mitochondrial DNA analysis

The European lynx (Lynx lynx) hasexperienced significant decline in populationnumbers over large parts of its formerdistribution area in central and northernEurope. In Scandinavia (Sweden and Norway), thespecies has been subject to intense hunting and inthe early 20th century the population size mayhave been as low as about 100 animals. Duringthe rest of the century there have beenalternating periods of restricted hunting andtotal protection. Future management of theScandinavian lynx population will requireinsight into what effects demographicbottlenecks may have had on genetic variabilityand structure. For this purpose, 276 lynxesfrom Sweden, Norway, Finland, Estonia andLatvia were analysed for polymorphism at 11feline microsatellite loci and at themitochondrial DNA (mtDNA) control region.Scandinavian lynxes were found to be fixed fora single mtDNA haplotype, while this and threeadditional haplotypes were seen in Finland andthe Baltic States (Estonia and Latvia); thehaplotypes were all very similar, onlydiffering at 1–4 sites within a 700 bp regionsequenced. Microsatellite variability wasmoderate (H_e = 0.51–0.62) with lowerheterozygosity and fewer alleles in Scandinaviathan in Finland and the Baltic States together,though significant so only for the latter.Heterozygosity data in Scandinavia wereconsistent with a recent population bottleneck.Various analyses (e.g. F_st, individual-basedtree, assignment test) revealed distinctgenetic differentiation between Scandinavianlynxes and animals from Finland and the BalticStates. Some structure was evident withinScandinavia as well, suggesting an isolation bydistance. The observed partition of geneticvariability between Scandinavia and the easterncountries thereof indicates that lynxpopulations from the two regions may need to beseparately managed. We discuss what factors canhave contributed to the population geneticstructure seen in northern European lynxpopulations of today.