Whole‐genome resequencing reveals the pleistocene temporal dynamics of Branchiostoma belcheri and Branchiostoma floridae populations

Global climatic fluctuations governed the ancestral demographic histories of species and contributed to place the current population status into a more extensive ecological and evolutionary context. Genetic variations will leave unambiguous signatures in the patterns of intraspecific genetic variation in extant species since the genome of each individual is an imperfect mosaic of the ancestral genomes. Here, we report the genome sequences of 20 Branchiostoma individuals by whole‐genome resequencing strategy. We detected over 140 million genomic variations for each Branchiostoma individual. In particular, we applied the pairwise sequentially Markovian coalescent (PSMC) method to estimate the trajectories of changes in the effective population size (N_e) of Branchiostoma population during the Pleistocene. We evaluated the threshold of sequencing depth for proper inference of demographic histories using PSMC was ≥25×. The PSMC results highlight the role of historical global climatic fluctuations in the long‐term population dynamics of Branchiostoma. The inferred ancestral N_e of the Branchiostoma belcheri populations from Zhanjiang and Xiamen (China) seawaters was different in amplitude before the first (mutation rate = 3 × 10^?9) or third glaciation (mutation rate = 9 × 10^?9) of the Pleistocene, indicating that the two populations most probably started to evolve in isolation in their respective seas after the first or third glaciation of the Pleistocene. A pronounced population bottleneck coinciding with the last glacial maximum was observed in all Branchiostoma individuals, followed by a population expansion occurred during the late Pleistocene. Species that have experienced long‐term declines may be especially vulnerable to recent anthropogenic activities. Recently, the industrial pollution and the exploitation of sea sand have destroyed the harmonious living environment of amphioxus species. In the future, we need to protect the habitat of Branchiostoma and make full use of these detected genetic variations to facilitate the functional study of Branchiostoma for adaptation to local environments.     

Effect of Yangtze River on population genetic structure of the relict plant Parrotia subaequalis in eastern China

Parrotia subaequalis (Hamamelidaceae) is a Tertiary relic species endemic in eastern China. We used inter‐simple sequence repeat (ISSR) markers to access genetic diversity and population genetic structure in natural five populations of P. subaequalis. The levels of genetic diversity were higher at species level (H = 0.2031) but lower at population level (H = 0.1096). The higher genetic diversity at species levels might be attributed to the accumulation of distinctive genotypes which adapted to the different habitats after Quaternary glaciations. Meanwhile, founder effects on the early stage, and subsequent bottleneck of population regeneration due to its biological characteristics, environmental features, and human activities, seemed to explain the low population levels of genetic diversity. The hierarchical AMOVA revealed high levels (42.60%) of among‐population genetic differentiation, which was in congruence with the high levels of Nei's genetic differentiation index (G_ST = 0.4629) and limited gene flow (N_m = 0.5801) among the studied populations. Mantel test showed a significant isolation‐by‐distance, indicating that geographic isolation has a significant effect on genetic structure in this species. Unweighted pair‐group method with arithmetic average clustering, PCoA, and Bayesian analyses uniformly recovered groups that matched the geographical distribution of this species. In particular, our results suggest that Yangtze River has served as a natural barrier to gene flow between populations occurred on both riversides. Concerning the management of P. subaequalis, the high genetic differentiation among populations indicates that preserving all five natural populations in situ and collecting enough individuals from these populations for ex situ conservation are necessary.     

Are habitat fragmentation,local adaptation and isolation‐by‐distance driving population divergence in wild rice Oryza rufipogon?

Habitat fragmentation weakens the connection between populations and is accompanied with isolation by distance (IBD) and local adaptation (isolation by adaptation, IBA), both leading to genetic divergence between populations. To understand the evolutionary potential of a population and to formulate proper conservation strategies, information on the roles of IBD and IBA in driving population divergence is critical. The putative ancestor of Asian cultivated rice (Oryza sativa) is endangered in China due to habitat loss and fragmentation. We investigated the genetic variation in 11 Chinese Oryza rufipogon populations using 79 microsatellite loci to infer the effects of habitat fragmentation, IBD and IBA on genetic structure. Historical and current gene flows were found to be rare (m_h = 0.0002–0.0013, m_c = 0.007–0.029), indicating IBD and resulting in a high level of population divergence (F_ST = 0.343). High within‐population genetic variation (H_E = 0.377–0.515), relatively large effective population sizes (N_e = 96–158), absence of bottlenecks and limited gene flow were found, demonstrating little impact of recent habitat fragmentation on these populations. Eleven gene‐linked microsatellite loci were identified as outliers, indicating local adaptation. Hierarchical AMOVA and partial Mantel tests indicated that population divergence of Chinese O. rufipogon was significantly correlated with environmental factors, especially habitat temperature. Common garden trials detected a significant adaptive population divergence associated with latitude. Collectively, these findings imply that IBD due to historical rather than recent fragmentation, followed by local adaptation, has driven population divergence in O. rufipogon.     

Genetic diversity and population structure of three Indian horse breeds

The genetic relationships of three Indian horse breeds—Marwari, Spiti, and Kathiawari were studied by genotyping 96 individuals with 20 polymorphic microsatellite markers. A total of 157 alleles were detected across 20 polymorphic loci. The Marwari population showed the highest allelic diversity (A = 5.7 and Ar = 5.14), followed by Spiti (A = 4.9 and Ar = 4.74) and Kathiawari (A = 4.1 and Ar = 3.82). The gene diversity was highest in the Spiti population (He = 0.67), followed by Marwari (He = 0.66) and Kathiawari (He = 0.59). Within population inbreeding estimates (f) in Marwari, Spiti and Kathiawari breeds were 0.18, 0.08, and 0.07, respectively, suggesting high level of inbreeding in these breeds. Analysis of bottleneck revealed evidence of recent bottleneck in Spiti and Kathiawari populations. Pair-wise Fst analysis, AMOVA and assignment tests demonstrated high genetic differentiation and low gene flow between populations. The information about genetic diversity and population structure will be useful for the future development of effective breeding management in order to preserve these Indian horse breeds.     

Trapped within the city: integrating demography,time since isolation and population‐specific traits to assess the genetic effects of urbanization

Urbanization is a severe form of habitat fragmentation that can cause many species to be locally extirpated and many others to become trapped and isolated within an urban matrix. The role of drift in reducing genetic diversity and increasing genetic differentiation is well recognized in urban populations. However, explicit incorporation and analysis of the demographic and temporal factors promoting drift in urban environments are poorly studied. Here, we genotyped 15 microsatellites in 320 fire salamanders from the historical city of Oviedo (Est. 8th century) to assess the effects of time since isolation, demographic history (historical effective population size; N_e) and patch size on genetic diversity, population structure and contemporary N_e. Our results indicate that urban populations of fire salamanders are highly differentiated, most likely due to the recent N_e declines, as calculated in coalescence analyses, concomitant with the urban development of Oviedo. However, urbanization only caused a small loss of genetic diversity. Regression modelling showed that patch size was positively associated with contemporary N_e, while we found only moderate support for the effects of demographic history when excluding populations with unresolved history. This highlights the interplay between different factors in determining current genetic diversity and structure. Overall, the results of our study on urban populations of fire salamanders provide some of the very first insights into the mechanisms affecting changes in genetic diversity and population differentiation via drift in urban environments, a crucial subject in a world where increasing urbanization is forecasted.     

Rise and fall of a wolf population: genetic diversity and structure during recovery,rapid expansion and drastic decline

The grey wolves (Canis lupus) of Finland have had a varied history, with a period of rapid population expansion after the mid‐1990s followed by a decline with a current census size of about 140 wolves. Here, we investigate the impact of unstable population size and connectivity on genetic diversity and structure in a long‐term genetic study of 298 Finnish wolves born in 1995–2009 and genotyped for 17 microsatellite loci. During the initial recovery and prior to population expansion, genetic diversity was high (1995–1997: LD‐N_e = 67.2; H_o = 0.749; H_e = 0.709) despite a small census size and low number of breeders (N_c < 100; N_b < 10) likely reflecting the status of the Russian source population. Surprisingly, observed heterozygosity decreased significantly during the study period (t = ?2.643, P = 0.021) despite population expansion, likely a result of an increase in inbreeding (F_IS = 0.108 in 2007–2009) owing to a low degree of connectivity with adjacent Russian wolf population (m = 0.016–0.090; F_ST = 0.086, P < 0.001) and population crash after 2006. However, population growth had a temporary positive impact on N_e and number of family lines. This study shows that even strong population growth alone might not be adequate to retain genetic diversity, especially when accompanied with low amount of subsequent gene flow and population decline.     

Approximate Bayesian estimation of extinction rate in the Finnish Daphnia magna metapopulation

Approximate Bayesian computation (ABC) is useful for parameterizing complex models in population genetics. In this study, ABC was applied to simultaneously estimate parameter values for a model of metapopulation coalescence and test two alternatives to a strict metapopulation model in the well‐studied network of Daphnia magna populations in Finland. The models shared four free parameters: the subpopulation genetic diversity (θ_S), the rate of gene flow among patches (4Nm), the founding population size (N₀) and the metapopulation extinction rate (e) but differed in the distribution of extinction rates across habitat patches in the system. The three models had either a constant extinction rate in all populations (strict metapopulation), one population that was protected from local extinction (i.e. a persistent source), or habitat‐specific extinction rates drawn from a distribution with specified mean and variance. Our model selection analysis favoured the model including a persistent source population over the two alternative models. Of the closest 750 000 data sets in Euclidean space, 78% were simulated under the persistent source model (estimated posterior probability = 0.769). This fraction increased to more than 85% when only the closest 150 000 data sets were considered (estimated posterior probability = 0.774). Approximate Bayesian computation was then used to estimate parameter values that might produce the observed set of summary statistics. Our analysis provided posterior distributions for e that included the point estimate obtained from previous data from the Finnish D. magna metapopulation. Our results support the use of ABC and population genetic data for testing the strict metapopulation model and parameterizing complex models of demography.     

Contrasting demographic history and gene flow patterns of two mangrove species on either side of the Central American Isthmus

Comparative phylogeography offers a unique opportunity to understand the interplay between past environmental events and life‐history traits on diversification of unrelated but co‐distributed species. Here, we examined the effects of the quaternary climate fluctuations and palaeomarine currents and present‐day marine currents on the extant patterns of genetic diversity in the two most conspicuous mangrove species of the Neotropics. The black (Avicennia germinans, Avicenniaceae) and the red (Rhizophora mangle, Rhizophoraceae) mangroves have similar geographic ranges but are very distantly related and show striking differences on their life‐history traits. We sampled 18 Atlantic and 26 Pacific locations for A. germinans (N = 292) and R. mangle (N = 422). We performed coalescence simulations using microsatellite diversity to test for evidence of population change associated with quaternary climate fluctuations. In addition, we examined whether patterns of genetic variation were consistent with the directions of major marine (historical and present day) currents in the region. Our demographic analysis was grounded within a phylogeographic framework provided by the sequence analysis of two chloroplasts and one flanking microsatellite region in a subsample of individuals. The two mangrove species shared similar biogeographic histories including: (1) strong genetic breaks between Atlantic and Pacific ocean basins associated with the final closure of the Central American Isthmus (CAI), (2) evidence for simultaneous population declines between the mid‐Pleistocene and early Holocene, (3) asymmetric historical migration with higher gene flow from the Atlantic to the Pacific oceans following the direction of the palaeomarine current, and (4) contemporary gene flow between West Africa and South America following the major Atlantic Ocean currents. Despite the remarkable differences in life‐history traits of mangrove species, which should have had a strong influence on seed dispersal capability and, thus, population connectivity, we found that vicariant events, climate fluctuations and marine currents have shaped the distribution of genetic diversity in strikingly similar ways.     

Influence of translocation strategy and mating system on the genetic structure of a newly established population of island ptarmigan

Island populations and populations established by reintroductions are prone to extinction, in part because they are vulnerable to deterministic and stochastic phenomena associated with geographic isolation and small population size. As population size declines, reduced genetic diversity can result in decreased fitness and reduced adaptive potential, which may hinder short- or long-term population viability. We used 32 microsatellite markers to investigate the conservation genetics of a newly established population of Evermann’s Rock Ptarmigan (Lagopus muta evermanni) at Agattu Island, in the western Aleutian Archipelago, Alaska. We found low genetic diversity (observed heterozygosity = 0.41, allelic richness = 2.2) and a small effective population size (N _e  = 28.6), but a relatively large N _e /N ratio = 0.55, which was attributed to multiple paternity in 80% of the broods and low reproductive skew among males (λ = 0.29). Moreover, successful breeding pairs were less related to each other than random male–female pairs. For conservation efforts based on reintroductions, a mating system with high rates of multiple paternity may facilitate retention of genetic diversity, thereby reducing the potential for inbreeding in small or isolated populations. Our results underscore the importance of quantifying genetic diversity and understanding the breeding behavior of translocated populations.     

Using genomic information for management planning of an endangered perennial,Viola uliginosa

Species occupying habitats subjected to frequent natural and/or anthropogenic changes are a challenge for conservation management. We studied one such species, Viola uliginosa, an endangered perennial wetland species typically inhabiting sporadically flooded meadows alongside rivers/lakes. In order to estimate genomic diversity, population structure, and history, we sampled five sites in Finland, three in Estonia, and one each in Slovenia, Belarus, and Poland using genomic SNP data with double‐digest restriction site‐associated DNA sequencing (ddRAD‐seq). We found monophyletic populations, high levels of inbreeding (mean population F_SNP = 0.407–0.945), low effective population sizes (N_e = 0.8–50.9), indications of past demographic expansion, and rare long‐distance dispersal. Our results are important in implementing conservation strategies for V. uliginosa, which should include founding of seed banks, ex situ cultivations, and reintroductions with individuals of proper origin, combined with continuous population monitoring and habitat management.     

Does origin always matter? Evaluating the influence of nonlocal seed provenances for ecological restoration purposes in a widespread and outcrossing plant species

For restoration purposes, nature conservation generally enforces the use of local seed material based on the “local‐is‐best” (LIB) approach. However, in some cases recommendations to refrain from this approach have been made. Here we test if a common widespread species with no obvious signs of local adaptation may be a candidate species for abandoning LIB during restoration. Using 10 microsatellite markers we compared population genetic patterns of the generalist species Daucus carota in indigenous and formerly restored sites (nonlocal seed provenances). Gene diversity overall ranged between H_e = 0.67 and 0.86 and showed no significant differences between the two groups. Hierarchical AMOVA and principal component analysis revealed very high genetic population admixture and negligible differentiation between indigenous and restored sites (F_CT = 0.002). Moreover, differentiation between groups was caused by only one outlier population, where inbreeding effects are presumed. We therefore conclude that the introduction of nonlocal seed provenances in the course of landscape restoration did not jeopardize regional species persistence by contributing to inbreeding or outbreeding depressions, or any measurable adverse population genetic effect. On the basis of these results, we see no obvious objections to the current practice to use the 10‐fold cheaper, nonlocal seed material of D. carota for restoration projects.     

Low genetic diversity and local adaptive divergence of Dracaena cambodiana (Liliaceae) populations associated with historical population bottlenecks and natural selection: an endangered long‐lived tree endemic to Hainan Island,China

Historical population bottlenecks and natural selection have important effects on the current genetic diversity and structure of long‐lived trees. Dracaena cambodiana is an endangered, long‐lived tree endemic to Hainan Island, China. Our field investigations showed that only 10 populations remain on Hainan Island and that almost all have been seriously isolated and grow in distinct habitats. A considerable amount of genetic variation at the species level, but little variation at the population level, and a high level of genetic differentiation among the populations with limited gene flow in D. cambodiana were detected using inter‐simple sequence repeat (ISSR) and random amplified polymorphic DNA (RAPD) analyses. No significant correlation was found between genetic diversity and actual population size, as the genetic diversities were similar regardless of population size. The Mantel test revealed that there was no correlation between genetic and geographic distances among the 10 populations. The UPGMA, PCoA and Bayesian analyses showed that local adaptive divergence has occurred among the D. cambodiana populations, which was further supported by habitat‐private fragments. We suggest that the current genetic diversity and population differentiation of D. cambodiana resulted from historical population bottlenecks and natural selection followed by historical isolation. However, the lack of natural regeneration of D. cambodiana indicates that former local adaptations with low genetic diversity may have been genetically weak and are unable to adapt to the current ecological environments.     

Genetic variation in nuclear and mitochondrial markers supports a large sex difference in lifetime reproductive skew in a lekking species

Sex differences in skews of vertebrate lifetime reproductive success are difficult to measure directly. Evolutionary histories of differential skew should be detectable in the genome. For example, male‐biased skew should reduce variation in the biparentally inherited genome relative to the maternally inherited genome. We tested this approach in lek‐breeding ruff (Class Aves, Philomachus pugnax) by comparing genetic variation of nuclear microsatellites (θ_n; biparental) versus mitochondrial D‐loop sequences (θ_m; maternal), and conversion to comparable nuclear (N_e) and female (N_ef) effective population size using published ranges of mutation rates for each marker (μ). We provide a Bayesian method to calculate N_e (θ_n = 4N_eμ_n) and N_ef (θ_m = 2N_efμ_m) using 95% credible intervals (CI) of θ_n and θ_m as informative priors, and accounting for uncertainty in μ. In 96 male ruffs from one population, N_e was 97% (79–100%) lower than expected under random mating in an ideal population, where N_e:N_ef = 2. This substantially lower autosomal variation represents the first genomic support of strong male reproductive skew in a lekking species.     

Using genetic evaluation to guide conservation of remnant Juniperus communis (Cupressaceae) populations

• Many critically endangered plant species exist in small, genetically depauperate or inbred populations, making assisted gene flow interventions necessary for long‐term population viability. However, before such interventions are implemented, conservation practitioners must consider the genetic and demographic status of extant populations, which are strongly affected by species’ life‐history traits. In northwestern Europe, Juniperus communis, a dioecious, wind‐pollinated and bird‐dispersed gymnosperm, has been declining for the past century and largely exists in small, isolated and senescent populations.
• To provide useful recommendations for a recovery plan involving translocation of plants, we investigated genetic diversity and structure of populations in Belgium using four microsatellite and five plastid single‐nucleotide polymorphism (SNP) markers.
• We detected no clonality in the populations, suggesting predominantly sexual reproduction. Populations exhibited high genetic diversity (H_e = 0.367–0.563) and low to moderate genetic differentiation (F_ST ≤ 0.133), with no clear geographic structure. Highly positive inbreeding coefficients (F_IS = 0.221–0.507) were explained by null alleles, population substructuring and biparental inbreeding. No isolation by distance was observed among distant populations, but isolation at close geographic proximity was found. Patterns were consistent with high historical gene flow through pollen and seed dispersal at both short and long distances. We also tested four pre‐germination treatments among populations to improve germination rates; however, germination rates remained low and only cold‐stratification treatments induced germination in some populations.
• To bolster population regeneration, introductions of cuttings from several source populations are recommended, in combination with in situ management practices that improve seedling survival and with ex situ propagation.

     

Using genetic variation to infer associations with climate in the common frog,Rana temporaria

Recent and historical species' associations with climate can be inferred using molecular markers. This knowledge of population and species‐level responses to climatic variables can then be used to predict the potential consequences of ongoing climate change. The aim of this study was to predict responses of Rana temporaria to environmental change in Scotland by inferring historical and contemporary patterns of gene flow in relation to current variation in local thermal conditions. We first inferred colonization patterns within Europe following the last glacial maximum by combining new and previously published mitochondrial DNA sequences. We found that sequences from our Scottish samples were identical to (92%), or clustered with, the common haplotype previously identified from Western Europe. This clade showed very low mitochondrial variation, which did not allow inference of historical colonization routes but did allow interpretation of patterns of current fine‐scale population structure without consideration of confounding historical variation. Second, we assessed fine‐scale microsatellite‐based patterns of genetic variation in relation to current altitudinal temperature gradients. No population structure was found within altitudinal gradients (average F_ST = 0.02), despite a mean annual temperature difference of 4.5 °C between low‐ and high‐altitude sites. Levels of genetic diversity were considerable and did not vary between sites. The panmictic population structure observed, even along temperature gradients, is a potentially positive sign for R. temporaria persistence in Scotland in the face of a changing climate. This study demonstrates that within taxonomic groups, thought to be at high risk from environmental change, levels of vulnerability can vary, even within species.     

Do plant populations on distinct inselbergs talk to each other? A case study of genetic connectivity of a bromeliad species in an Ocbil landscape

Here, we explore the historical and contemporaneous patterns of connectivity among Encholirium horridum populations located on granitic inselbergs in an Ocbil landscape within the Brazilian Atlantic Forest, using both nuclear and chloroplast microsatellite markers. Beyond to assess the E. horridum population genetic structure, we built species distribution models across four periods (current conditions, mid‐Holocene, Last Glacial Maximum [LGM], and Last Interglacial) and inferred putative dispersal corridors using a least‐cost path analysis to elucidate biogeographic patterns. Overall, high and significant genetic divergence was estimated among populations for both nuclear and plastid DNA (Φ_ST(n) = 0.463 and Φ_ST(plastid) = 0.961, respectively, p < .001). For nuclear genome, almost total absence of genetic admixture among populations and very low migration rates were evident, corroborating with the very low estimates of immigration and emigration rates observed among E. horridum populations. Based on the cpDNA results, putative dispersal routes in Sugar Loaf Land across cycles of climatic fluctuations in the Quaternary period revealed that the populations’ connectivity changed little during those events. Genetic analyses highlighted the low genetic connectivity and long‐term persistence of populations, and the founder effect and genetic drift seemed to have been very important processes that shaped the current diversity and genetic structure observed in both genomes. The genetic singularity of each population clearly shows the need for in situ conservation of all of them.     

High genetic diversity in a small population: the case of Chilean blue whales

It is generally assumed that species with low population sizes have lower genetic diversities than larger populations and vice versa. However, this would not be the case for long‐lived species with long generation times, and which populations have declined due to anthropogenic effects, such as the blue whale (Balaenoptera musculus). This species was intensively decimated globally to near extinction during the 20th century. Along the Chilean coast, it is estimated that at least 4288 blue whales were hunted from an apparently pre‐exploitation population size (k) of a maximum of 6200 individuals (Southeastern Pacific). Thus, here, we describe the mtDNA (control region) and nDNA (microsatellites) diversities of the Chilean blue whale aggregation site in order to verify the expectation of low genetic diversity in small populations. We then compare our findings with other blue whale aggregations in the Southern Hemisphere. Interestingly, although the estimated population size is small compared with the pre‐whaling era, there is still considerable genetic diversity, even after the population crash, both in mitochondrial (N = 46) and nuclear (N = 52) markers (H_d = 0.890 and H_o = 0.692, respectively). Our results suggest that this diversity could be a consequence of the long generation times and the relatively short period of time elapsed since the end of whaling, which has been observed in other heavily‐exploited whale populations. The genetic variability of blue whales on their southern Chile feeding grounds was similar to that found in other Southern Hemisphere blue whale feeding grounds. Our phylogenetic analysis of mtDNA haplotypes does not show extensive differentiation of populations among Southern Hemisphere blue whale feeding grounds. The present study suggests that although levels of genetic diversity are frequently used as estimators of population health, these parameters depend on the biology of the species and should be taken into account in a monitoring framework study to obtain a more complete picture of the conservation status of a population.     

Genetic substructure and admixture as important factors in linkage disequilibrium‐based estimation of effective number of breeders in recovering wildlife populations

The number of effective breeders (N_b) and effective population size (N_e) are population parameters reflective of evolutionary potential, susceptibility to stochasticity, and viability. We have estimated these parameters using the linkage disequilibrium‐based approach with LDNE through the latest phase of population recovery of the brown bears (Ursus arctos) in Finland (1993–2010; N = 621). This phase of the recovery was recently documented to be associated with major changes in genetic composition. In particular, differentiation between the northern and the southern genetic cluster declined rapidly within 1.5 generations. Based on this, we have studied effects of the changing genetic structure on N_b and N_e, by comparing estimates for whole Finland with the estimates for the two genetic clusters. We expected a potentially strong relationship between estimate sizes and genetic differentiation, which should disappear as the population recovers and clusters merge. Consistent with this, our estimates for whole Finland were lower than the sum of the estimates of the two genetic clusters and both approaches produced similar estimates in the end. Notably, we also found that admixed genotypes strongly increased the estimates. In all analyses, our estimates for N_e were larger than N_b and likely reflective for brown bears of the larger region of Finland and northwestern Russia. Conclusively, we find that neglecting genetic substructure may lead to a massive underestimation of N_b and N_e. Our results also suggest the need for further empirical analysis focusing on individuals with admixed genotypes and their potential high influence on N_b and N_e.     

Temporal changes in genetic diversity of isolated populations of perch and roach

Genetic drift, together with natural selection and gene flow, affects genetic variation and is the major source of changes in allele frequencies in small and isolated populations. Temporal shifts in allele frequencies at five polymorphic loci were used to estimate the amount of genetic drift in an isolated population of perch (Perca fluviatilis L.) and roach (Rutilus rutilus L.). Here, I used the populations from the Biotest basin at Forsmark, Sweden, to investigate genetic diversity between 1977 and 2000, during which time the population can be considered to be totally isolated from other populations. Microsatellite data reveal stable levels of gene diversity over time for both species. Estimates of genetic differentiation (F _ST) showed a significant divergence between 1977 and 2000 for both perch and roach. A positive correlation between genetic distance and time was found (Mantel test, perch: r = 0.724, P = 0.0112; roach: r = 0.59, P = 0.036). Estimates of effective population size (N _e) differed with a factor six between two different estimators (NeEstimator and TempoFS) applying the temporal method. Ratios of N _e/N ranged between 10⁻² and 10⁻³, values normally found in marine species. Despite low N _e the populations have not lost their evolutionary potential due to drift. But two decades of isolation have lead to isolation by time for populations of perch and roach, respectively.     

Genetic variation and effective population size in isolated populations of coastal cutthroat trout

Following glacial recession in southeast Alaska, waterfalls created by isostatic rebound have isolated numerous replicate populations of coastal cutthroat trout (Oncorhynchus clarkii clarkii) in short coastal streams. These replicate isolated populations offer an unusual opportunity to examine factors associated with the maintenance of genetic diversity. We used eight microsatellites to examine genetic variation within and differentiation among 12 population pairs sampled from above and below these natural migration barriers. Geological evidence indicated that the above-barrier populations have been isolated for 8,000–12,500 years. Genetic differentiation among below-barrier populations (F _ST = 0.10, 95% C.I. 0.08–0.12) was similar to a previous study of more southern populations of this species. Above-barrier populations were highly differentiated from adjacent below-barrier populations (mean pairwise F _ST = 0.28; SD 0.18) and multiple lines of evidence were consistent with asymmetric downstream gene flow that varied among streams. Each above-barrier population had reduced within-population genetic variation when compared to the adjacent below-barrier population. Within-population genetic diversity was significantly correlated with the amount of available habitat in above-barrier sites. Increased genetic differentiation of above-barrier populations with lower genetic diversity suggests that genetic drift has been the primary cause of genetic divergence. Long-term estimates of N _e based on loss of heterozygosity over the time since isolation were large (3,170; range 1,077–7,606) and established an upper limit for N _e if drift were the only evolutionary process responsible for loss of genetic diversity. However, it is likely that a combination of mutation, selection, and gene flow have also contributed to the genetic diversity of above-barrier populations. Contemporary above-barrier N _e estimates were much smaller than long-term N _e estimates, not correlated with within-population genetic diversity, and not consistent with the amount of genetic variation retained, given the approximate 10,000-year period of isolation. The populations isolated by waterfalls in this study that occur in larger stream networks have retained substantial genetic variation, which suggests that the amount of habitat in headwater streams is an important consideration for maintaining the evolutionary potential of isolated populations.