Determinants of Genetic Structure in a Nonequilibrium Metapopulation of the Plant Silene latifolia

Population genetic differentiation will be influenced by the demographic history of populations, opportunities for migration among neighboring demes and founder effects associated with repeated extinction and recolonization. In natural populations, these factors are expected to interact with each other and their magnitudes will vary depending on the spatial distribution and age structure of local demes. Although each of these effects has been individually identified as important in structuring genetic variance, their relative magnitude is seldom estimated in nature. We conducted a population genetic analysis in a metapopulation of the angiosperm, Silene latifolia, from which we had more than 20 years of data on the spatial distribution, demographic history, and extinction and colonization of demes. We used hierarchical Bayesian methods to disentangle which features of the populations contributed to among population variation in allele frequencies, including the magnitude and direction of their effects. We show that population age, long-term size and degree of connectivity all combine to affect the distribution of genetic variance; small, recently-founded, isolated populations contributed most to increase F _ST in the metapopulation. However, the effects of population size and population age are best understood as being modulated through the effects of connectivity to other extant populations, i.e. F _ST diminishes as populations age, but at a rate that depends how isolated the population is. These spatial and temporal correlates of population structure give insight into how migration, founder effect and within-deme genetic drift have combined to enhance and restrict genetic divergence in a natural metapopulation.     

PHASE III OF WRIGHT'S SHIFTING BALANCE PROCESS AND THE VARIANCE AMONG DEMES IN MIGRATION RATE

Interdemic selection by the differential migration of individuals out from demes of high fitness and into demes of low fitness (Phase III) is one of the most controversial aspects of Wright's Shifting Balance Theory. I derive a relationship between Phase III migration and the interdemic selection differential, S, and show its potential effect on F_ST. The relationship reveals a diversifying effect of interdemic selection by Phase III migration on the genetic structure of a metapopulation. Using experimental metapopulations, I explored the effect of Phase III migration on F_ST by comparing the genetic variance among demes for two different patterns of migration: (1) island model migration and (2) Wright's Phase III migration. Although mean migration rates were the same, I found that the variance among demes in migration rate was significantly higher with Phase III than with island model migration. As a result, F_ST for the frequency of a neutral marker locus was higher with Phase III than it was with island model migration. By increasing F_ST, Phase III enhanced the genetic differentiation among demes for traits not subject to interdemic selection. This feature makes Wright's process different from individual selection which, by reducing effective population size, decreases the genetic variance within demes for all other traits. I discussed this finding in relation to the efficacy of Phase III and random migration for effecting peak shifts, and the contribution of genes with indirect effects to among‐deme variation.     

Genetic load,inbreeding depression and heterosis in an age‐structured metapopulation

In a metapopulation, the process of recurrent local extinction and recolonization gives rise to an age structure among demes. Recently established demes will tend to differ from older demes in terms of the levels of genetic diversity found within them and the way this diversity is distributed among demes in the same and different ages. The effects of population turnover on average levels of genetic diversity among demes in a metapopulation have been the focus of much attention, both for neutral and nonneutral loci, but much less is known about the distribution of nonneutral genetic diversity among demes of different ages. In this paper, we used computer simulations to study the distribution of genetic load, inbreeding depression and heterosis in an age‐structured metapopulation. We found that, for mildly deleterious mutations, within‐deme inbreeding depression increased, whereas heterosis and genetic load decreased with deme age following severe colonization bottlenecks. In contrast, recessive lethal alleles tended to be purged during colonization, with older populations showing higher genetic load and higher within‐deme inbreeding depression. Heterosis caused by recessive lethal alleles and resulting from gene flow among different demes tended to be greatest for young demes, because the mutations responsible tended to be purged in the first few generations after colonization, but its effects increased again as populations grow older as a result of immigration. Our results point to a need for estimates of genetic diversity, genetic load, within‐deme inbreeding depression and heterosis in demes of different age classes separately.     

Coalescence in a metapopulation with recurrent local extinction and recolonization

Abstract Many species exist as metapopulations in balance between local population extinction and recolonization. The effect of these processes on average population differentiation, within-deme diversity, and specieswide diversity has been considered previously. In this paper, coalescent simulations of Slatkin's propagule-pool and migrant-pool models are used to characterize the distribution of neutral genetic diversity within demes (π_s), diversity in the metapopulation a whole (TT_T), the ratio F _ST= (π_t–π_S)/π_T, Tajima's D statistic, and several ratios of gene-tree branch lengths. Using these distributions, power to detect differences in key metapopulation parameter values is determined under contrasting sampling regimes. The results indicate that it will be difficult to use sequence data from a single locus to detect a history of extinctions and recolonizations in a metapopulation because of high genealogical variance, the loss of diversity due to reductions in effective population size, and the fact that a genealogy of lineages from different demes under Slatkin's model differs from a neutral coalescent only in its time scale. Genetic indices of gene-tree shape that capture the effects of extinction/recolonization on both external branches and the length of the genealogy as a whole will provide the best indication of metapopulation dynamics if several lineages are sampled from several different demes.     

Parasite‐mediated selection in a natural metapopulation of Daphnia magna

Parasite‐mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease‐related traits. However, nonadaptive processes like migration and extinction‐(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life‐history and disease‐related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite‐mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life‐history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST–F_ST‐like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction‐(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.     

Effects of metapopulation processes on measures of genetic diversity

Many species persist as a metapopulation under a balance between the local extinction of subpopulations or demes and their recolonization through dispersal from occupied patches. Here we review the growing body of literature dealing with the genetic consequences of such population turnover. We focus our attention principally on theoretical studies of a classical metapopulation with a 'finite-island' model of population structure, rather than on 'continent-island' models or 'source-sink' models. In particular, we concern ourselves with the subset of geographically subdivided population models in which it is assumed that all demes are liable to extinction from time to time and that all demes receive immigrants. Early studies of the genetic effects of population turnover focused on population differentiation, such as measured by F(ST). A key advantage of F(ST) over absolute measures of diversity is its relative independence of the mutation process, so that different genes in the same species may be compared. Another advantage is that F(ST) will usually equilibrate more quickly following perturbations than will absolute levels of diversity. However, because F(ST) is a ratio of between-population differentiation to total diversity, the genetic effects of metapopulation processes may be difficult to interpret in terms of F(ST) on its own, so that the analysis of absolute measures of diversity in addition is likely to be informative. While population turnover may either increase or decrease F(ST), depending on the mode of colonization, recurrent extinction and recolonization is expected always to reduce levels of both within-population and species-wide diversity (piS and piT, respectively). One corollary of this is that piS cannot be used as an unbiased estimate of the scaled mutation rate, theta, as it can, with some assumptions about the migration process, in species whose demes do not fluctuate in size. The reduction of piT in response to population turnover reflects shortened mean coalescent times, although the distribution of coalescence times under extinction colonization equilibrium is not yet known. Finally, we review current understanding of the effect of metapopulation dynamics on the effective population size.     

Low genetic diversity and strong but shallow population differentiation suggests genetic homogenization by metapopulation dynamics in a social spider

Mating systems and population dynamics influence genetic diversity and structure. Species that experience inbreeding and limited gene flow are expected to evolve isolated, divergent genetic lineages. Metapopulation dynamics with frequent extinctions and colonizations may, on the other hand, deplete and homogenize genetic variation, if extinction rate is sufficiently high compared to the effect of drift in local demes. We investigated these theoretical predictions empirically in social spiders that are highly inbred. Social spiders show intranest mating, female‐biased sex ratio, and frequent extinction and colonization events, factors that deplete genetic diversity within nests and populations and limit gene flow. We characterized population genetic structure in Stegodyphus sarasinorum, a social spider distributed across the Indian subcontinent. Species‐wide genetic diversity was estimated over approximately 2800 km from Sri Lanka to Himalayas, by sequencing 16 protein‐coding nuclear loci. We found 13 SNPs in 6592 bp (π = 0.00045) indicating low species‐wide nucleotide diversity. Three genetic lineages were strongly differentiated; however, only one fixed difference among them suggests recent divergence. This is consistent with a scenario of metapopulation dynamics that homogenizes genetic diversity across the species' range. Ultimately, low standing genetic variation may hamper a species' ability to track environmental change and render social inbreeding spiders ‘evolutionary dead‐ends’.     

Human‐facilitated metapopulation dynamics in an emerging pest species,Cimex lectularius

The number and demographic history of colonists can have dramatic consequences for the way in which genetic diversity is distributed and maintained in a metapopulation. The bed bug (Cimex lectularius) is a re‐emerging pest species whose close association with humans has led to frequent local extinction and colonization, that is, to metapopulation dynamics. Pest control limits the lifespan of subpopulations, causing frequent local extinctions, and human‐facilitated dispersal allows the colonization of empty patches. Founder events often result in drastic reductions in diversity and an increased influence of genetic drift. Coupled with restricted migration, this can lead to rapid population differentiation. We therefore predicted strong population structuring. Here, using 21 newly characterized microsatellite markers and approximate Bayesian computation (ABC), we investigate simplified versions of two classical models of metapopulation dynamics, in a coalescent framework, to estimate the number and genetic composition of founders in the common bed bug. We found very limited diversity within infestations but high degrees of structuring across the city of London, with extreme levels of genetic differentiation between infestations (F_ST = 0.59). ABC results suggest a common origin of all founders of a given subpopulation and that the numbers of colonists were low, implying that even a single mated female is enough to found a new infestation successfully. These patterns of colonization are close to the predictions of the propagule pool model, where all founders originate from the same parental infestation. These results show that aspects of metapopulation dynamics can be captured in simple models and provide insights that are valuable for the future targeted control of bed bug infestations.     

TEMPORAL FLUCTUATIONS IN DEMOGRAPHIC PARAMETERS AND THE GENETIC VARIANCE AMONG POPULATIONS

Contrary to assumptions commonly made in the study of population genetics, the demographic properties of many populations are not always constant. Important characteristics of populations such as migration rate and population size may vary in time and space. Moreover, local populations often come and go; the rate of extinction and the properties of colonization may also vary. In this paper, the approach to equilibrium following a disturbance in the genetic variance among populations is described. The rate of migration is shown to be critical in determining the extent to which extinction and recolonization affects genetic differentiation. Perturbations and variations through time and space in demographic parameters such as population size and migration rate are shown to be important in determining the partitioning of genetic variance. Equations are given to predict the average through time of genetic differentiation among populations in the event of a single disturbance or in constant fluctuations in the pertinent demographic parameters. In general, these fluctuations increase the F_ST of a species. Spatial demographic variation affects F_STmuch more than temporal variation. These demographic properties make some species unsuitable for the empirical analysis of migration with indirect genetic measures. Demographic instability may play a large role in the evolution of genetic variation.     

Retrospective coalescent methods and the reconstruction of metapopulation histories in the sea

Phylogeographic analyses are a key interface between ecological and evolutionary ways of knowing because such analyses integrate the cumulative effects of demographic (ecological) processes over geological (evolutionary) time scales. Newly developed coalescent methods allow evolutionary ecologists to overcome some limitations associated with inferring population history from classic methods such as Wright’s F _ST. Here we briefly contrast classic and coalescent methods for looking backward in time through a population genetic lens, focusing on the key advantages of the isolation-with-migration (IM) class of coalescent methods for distinguishing ancient connectedness from actual recurrent contemporary gene flow as causes of genetic similarity or differentiation among populations. Making this critical distinction can lead to the discovery of otherwise obscured histories underlying conventional patterns of spatial variation. We illustrate the importance of these insights using analyses of Pacific fishes, snails, and sea stars in which population sizes and divergence times are more important than rates of contemporary gene flow as determinants of population genetic differentiation. We then extend the IM method to genetic data from two model metapopulation species (California abalone, Australian damselfish). The analyses show the potential use of non-equilibrium IM methods for differentiating among metapopulation models that make different predictions about population parameters and have different implications for the design of marine protected areas and other conservation goals. At face value, the results largely rule out classic metapopulation dynamics (dominated by extinction and colonization rather than connectivity via ongoing recurrent gene flow) but, at the same time, do not strongly support a modern marine metapopulation dynamic (ecologically significant connectivity between demes). However, the results also highlight the need for much more data (i.e., loci) sampled on different spatial scales in order to determine whether metapopulation dynamics might exist on smaller scales than are typically sampled by most phylogeographers and landscape geneticists.     

Is local selection so widespread in river organisms? Fractal geometry of river networks leads to high bias in outlier detection

Identifying local adaptation is crucial in conservation biology to define ecotypes and establish management guidelines. Local adaptation is often inferred from the detection of loci showing a high differentiation between populations, the so‐called F_ST outliers. Methods of detection of loci under selection are reputed to be robust in most spatial population models. However, using simulations we showed that F_ST outlier tests provided a high rate of false‐positives (up to 60%) in fractal environments such as river networks. Surprisingly, the number of sampled demes was correlated with parameters of population genetic structure, such as the variance of F_STs, and hence strongly influenced the rate of outliers. This unappreciated property of river networks therefore needs to be accounted for in genetic studies on adaptation and conservation of river organisms.     

Population genetic structure and diversity of the endangered Cantabrian capercaillie

The Cantabrian capercaillie (Tetrao urogallus cantabricus) occupies the southwestern edge of the grouse family distribution range in Eurasia. It is endemic to the Cantabrian Mountains in northwestern Spain and is geographically isolated and separated from the neighboring population in the Pyrenees by a distance of 300 km. Over the last decades, the population has undergone a dramatic decline and is now threatened with extinction. This study presents the genetic analysis of the Cantabrian capercaillie population using non-invasive samples. We performed genotyping of 45 individuals using 20 microsatellites and a sex marker. The data highlight the need for using a large number of markers when considering fragmented small populations. Genetic diversity (H_E = 0.50) and average number of alleles (3.40) in the population were low. The population is fragmented into 2 clusters (F_ST = 0.113) that fit with areas on both sides of the transportation ways that divide its range. Both clusters exhibited additional heterozygote deficits. Geographical distance was negatively correlated with genetic relatedness (r = −0.44, P ≤ 0.001). The data show a recent decline in effective population size that can be related to an ongoing process of population reduction and fragmentation. Conservation actions should focus on the protection of local demes by maintaining a dense network of suitable patches to maximize reproductive output and the number of potential dispersers to reconnect the 2 subpopulations. © 2012 The Wildlife Society.     

Novel tetranucleotide microsatellite DNA markers for the wood frog,Rana sylvatica

Fifteen tetranucleotide microsatellite loci were identified and characterized for wood frogs (Rana sylvatica) collected from three vernal pools in the southeastern US. These markers revealed a high degree of genetic diversity (nine to 34 alleles per locus), heterozygosity (30.6–92.3%) and allelic heterogeneity (69% of comparisons were statistically significant). Considerable differentiation among populations was observed as genetic distances (chord) ranged between 0.40 and 0.55 and all F_ST values (0.02–0.05) were statistically significant. Genotypic assignment tests correctly classified 103 of 113 individuals to their respective collection. These markers should prove useful for investigating fine‐scale population structure and metapopulation dynamics.     

THE ROLE OF DEME SIZE,REPRODUCTIVE PATTERNS,AND DISPERSAL IN THE DYNAMICS OF t-LETHAL HAPLOTYPES

The t-lethal haplotypes (t) found in house mouse (Mus musculus) populations are recessive lethals favored by gametic selection whereby male heterozygotes exhibit a non-Mendelian transmission ratio of about 95% t. The expected equilibrium frequency is 0.385; however, empirical values are lower, averaging close to 0.13. We examined the hypothesis that interdemic selection is the cause of the low empirical values by using a deme-structured simulation model that included overlapping generations, a realistic breeding system, differential deme productivity, and a large total population. We found that under some conditions interdemic selection could lower t frequency below 0.13 in the face of immigration rates up to 5%. Low frequencies were correlated with effective deme size (n_e), regardless of whether n_e was changed through changing deme size (n) or through changing the proportion of breeding adults. Earlier workers showed how the first two phases of interdemic selection (random genetic differentiation and mass selection) interacted to reduce the haplotype frequency, but here we show the importance of the third phase (differential productivity of demes) once demes are linked by dispersal. The effect of this phase is not due to the (negative) covariation between deme productivity and haplotype frequency, but occurs when differential deme productivity generates a difference in t frequency between the population of juveniles recruited into their natal deme and the population of juvenile dispersers. This difference was maximized when the average productivity of demes was low, either because few adult females bred at any one time and/or because fecundity was low. Contrary to an earlier prediction, male-biased dispersal also reduced haplotype frequency, and this probably stems from the relative excess of wild-type genotypes among dispersers compared to the deme residents. Another unexpected finding was that the randomly generated excess of heterozygotes (F_IS < 0) found in small demes favored t haplotypes; however, the effect was only seen when the more powerful influence of the third phase of interdemic selection was removed. Simulations of neutral polymorphisms showed that a deme structure giving F_ST ≤ 0.6 is inconsistent with a haplotype frequency below 0.13. Based on current empirical estimates of F_ST (about 0.2), we concluded that immigration rates in the field are too high for interdemic selection alone to cause the observed deficit of lethal haplotypes. One factor that could combine with population structure effects is the observation that the transmission ratio is lowered to around 0.6 in litters produced from postpartum estrus (PPE). Incorporating this factor, we showed that interdemic selection could be effective in lowering the frequency of t below 0.13 when F_ST was above 0.43 even when migration rates were up to 10%. These results suggest that if empirical haplotype and F_ST estimates are accurate, then additional factors such as a lowered fitness of heterozygotes may be involved.     

Isolated in an ocean of grass: low levels of gene flow between termite subpopulations

Habitat fragmentation is one of the most important causes of biodiversity loss, but many species are distributed in naturally patchy habitats. Such species are often organized in highly dynamic metapopulations or in patchy populations with high gene flow between subpopulations. Yet, there are also species that exist in stable patchy habitats with small subpopulations and presumably low dispersal rates. Here, we present population genetic data for the ‘magnetic’ termite Amitermes meridionalis, which show that short distances between subpopulations do not hinder exceptionally strong genetic differentiation (F_ST: 0.339; R_ST: 0.636). Despite the strong genetic differentiation between subpopulations, we did not find evidence for genetic impoverishment. We propose that loss of genetic diversity might be counteracted by a long colony life with low colony turnover. Indeed, we found evidence for the inheritance of colonies by so‐called ‘replacement reproductives’. Inhabiting a mound for several generations might result in loss of gene diversity within a colony but maintenance of gene diversity at the subpopulation level.     

Isolation by Time During an Arctic Phytoplankton Spring Bloom

The arctic phytoplankton spring bloom, which is often diatom‐dominated, is a key event that provides the high latitude communities with a fundamental flux of organic carbon. During a bloom, phytoplankton may increase its biomass by orders of magnitude within days. Yet, very little is known about phytoplankton bloom dynamics, including for example how blooming affects genetic composition and diversity of a population. Here, we quantified the genetic composition and temporal changes of the diatom Fragilariopsis cylindrus, which is one of the most important primary producers in the Arctic, during the spring bloom in western Greenland, using 13 novel microsatellite markers developed for this study. We found that genetic differentiation (quantified using sample‐specific F_ST) decreased between time points as the bloom progressed, with the most drastic changes in F_ST occurring at the start of the bloom; thus the genetic structure of the bloom is characterized by isolation by time. There was little temporal variation in genetic diversity throughout the bloom (mean H_E = 0.57), despite marked fluctuations in F. cylindrus cell concentrations and the temporal change in sample‐specific F_ST. On the basis of this novel pattern of genetic differentiation, we suggest that blooming behavior may promote genetic diversity of a phytoplankton population.     

Background selection and FST: Consequences for detecting local adaptation

Background selection is a process whereby recurrent deleterious mutations cause a decrease in the effective population size and genetic diversity at linked loci. Several authors have suggested that variation in the intensity of background selection could cause variation in F_ST across the genome, which could confound signals of local adaptation in genome scans. We performed realistic simulations of DNA sequences, using recombination maps from humans and sticklebacks, to investigate how variation in the intensity of background selection affects F_ST and other statistics of population differentiation in sexual, outcrossing species. We show that, in populations connected by gene flow, Weir and Cockerham's (1984; Evolution, 38 , 1358) estimator of F_ST is largely insensitive to locus‐to‐locus variation in the intensity of background selection. Unlike F_ST, however, d_XY is negatively correlated with background selection. Moreover, background selection does not greatly affect the false‐positive rate in F_ST outlier studies in populations connected by gene flow. Overall, our study indicates that background selection will not greatly interfere with finding the variants responsible for local adaptation.     

Population structure and colonization of Bactrocera dorsalis (Diptera: Tephritidae) in China,inferred from mtDNA COI sequences

The oriental fruit fly, Bactrocera dorsalis, is a serious pest of fruits and vegetables in South‐east Asia, and, because of quarantine restrictions, impedes international trade and economic development in the region. Revealing genetic variation in oriental fruit fly populations will provide a better understanding of the colonization process and facilitate the quarantine and management of this species. The genetic structure in 15 populations of oriental fruit fly from southern China, Laos and Myanmar in South‐east Asia was examined with a 640‐bp sequence of the mitochondrial cytochrome oxidase subunit I (COI) gene. The highest levels of genetic diversity were found in Laos and Myanmar. Low to medium levels of genetic differentiation (F_ST ≤ 0.134) were observed among populations. Pooled populations from mainland China differed from those in Laos and Myanmar (F_ST = 0.024). Genetic structure across the region did not follow the isolation‐by‐distance model. The high genetic diversity observed in Laos and Myanmar supports the South‐east Asian origin of B. dorsalis. High genetic diversity and significant differentiation between some populations within mainland China indicate B. dorsalis populations have been established in the region for an extended period of time. High levels of genetic diversity observed among the five populations from Hainan Island and similarity between the Island and Chinese mainland populations indicate that B. dorsalis was introduced to Hainan from the mainland and has been on the island for many years. High genetic diversity in the recently established population in Shanghai (Pudong) suggests multiple introductions or a larger number of founders.     

THE EFFECTIVE SIZE OF A HIERARCHICALLY STRUCTURED POPULATION

A knowledge of the effective size of a population (N_e) is important in understanding its current and future evolutionary potential. Unfortunately, the effective size of a hierarchically structured population is not, in general, equal to the sum of its parts. In particular, the inbreeding structure has a major influence on N_e. Here I link N_e to Wright's hierarchical measures of inbreeding, F_IS and F_ST, for an island-structured population (or metapopulation) of size N_T. The influence of F_ST depends strongly on the degree to which island productivity is regulated. In the absence of local regulation (the interdemic model), interdemic genetic drift reduces N_e. When such drift is combined with local inbreeding under otherwise ideal conditions, the effects of F_IS and F_ST are identical: increasing inbreeding either within or between islands reduces N_e, with N_e = N_T/[(1 + F_IS)(1 + F_ST) ? 2F_ISF_ST]. However, if islands are all equally productive because of local density regulation (the traditional island model), then N_e = N_T/[(1 + F_IS)(1 –F_ST)] and the effect of F_ST is reversed. Under the interdemic model, random variation in the habitat quality (and hence productivity) of islands act to markedly decrease N_e. This variation has no effect under the island model because, by definition, all islands are equally productive. Even when no permanent island structure exists, spatial differences in habitat quality can significantly increase the overall variance in reproductive success of both males and females and hence lower N_e. Each of these basic results holds when other nonideal factors are added to the model. These factors, deviations from a 1:1 sex ratio, greater than Poisson variance in female reproductive success, and variation in male mating success due to polygynous mating systems, all act to lower N_e. The effects of male and female variance on N_e have important differences because only females affect island productivity. Finally, it is noted that to use these relationships, F_IS and F_ST must be estimated according to Wright's definition (and corrected to have a zero expectation under the null model). A commonly used partitioning (θ, θ_g) can be biased if either island size or the number of islands is small.     

Contemporary habitat loss reduces genetic diversity in an ecologically specialized butterfly

Aim This study investigated the influence of contemporary habitat loss on the genetic diversity and structure of animal species using a common, but ecologically specialized, butterfly, Theclinesthes albocincta (Lepidoptera: Lycaenidae), as a model. Location South Australia. Methods We used amplified fragment length polymorphism (AFLP) and allozyme datasets to investigate the genetic structure and genetic diversity among populations of T. albocincta in a fragmented landscape and compared this diversity and structure with that of populations in two nearby landscapes that have more continuous distributions of butterflies and their habitat. Butterflies were sampled from 15 sites and genotyped, first using 363 informative AFLP bands and then using 17 polymorphic allozyme loci (n = 248 and 254, respectively). We complemented these analyses with phylogeographic information based on mitochondrial DNA (mtDNA) haplotype information derived from a previous study in the same landscapes. Results Both datasets indicated a relatively high level of genetic structuring across the sampling range (AFLP, F_ST = 0.34; allozyme, F_ST = 0.13): structure was greatest among populations in the fragmented landscape (AFLP, F_ST = 0.15; allozyme, F_ST = 0.13). Populations in the fragmented landscape also had significantly lower genetic diversity than populations in the other two landscapes: there were no detectable differences in genetic diversity between the two continuous landscapes. There was also evidence (r² = 0.33) of an isolation by distance effect across the sampled range of the species. Main conclusions The multiple lines of evidence, presented within a phylogeographic context, support the hypothesis that contemporary habitat fragmentation has been a major driver of genetic erosion and differentiation in this species. Theclinesthes albocincta populations in the fragmented landscape are thus likely to be at greater risk of extinction because of reduced genetic diversity, their isolation from conspecific subpopulations in other landscapes, and other extrinsic forces acting on their small population sizes. Our study provides compelling evidence that habitat loss and fragmentation have significant rapid impacts on the genetic diversity and structure of butterfly populations, especially specialist species with particular habitat preferences and poor dispersal abilities.