Effective Population Size,Extended Linkage Disequilibrium and Signatures of Selection in the Rare Dog Breed Lundehund

The Lundehund is an old dog breed with remarkable anatomical features including polydactyly in all four limbs and extraordinary flexibility of the spine. We genotyped 28 Lundehund using the canine Illumina high density beadchip to estimate the effective population size (N_e) and inbreeding coefficients as well as to identify potential regions of positive selection. The decay of linkage disequilibrium was slow with r² = 0.95 in 50 kb distance. The last 7-200 generations ago, Ne was at 10-13. An increase of N_e was noted in the very recent generations with a peak value of 19 for N_e at generation 4. The FROH estimated for 50-, 65- and 358-SNP windows were 0.87, 087 and 0.81, respectively. The most likely estimates for F_ROH after removing identical-by-state segments due to linkage disequilibria were at 0.80-0.81. The extreme loss of heterozygosity has been accumulated through continued inbreeding over 200 generations within a probably closed population with a small effective population size. The mean inbreeding coefficient based on pedigree data for the last 11 generations (F_Ped = 0.10) was strongly biased downwards due to the unknown coancestry of the founders in this pedigree data. The long-range haplotype test identified regions with genes involved in processes of immunity, olfaction, woundhealing and neuronal development as potential targets of selection. The genes QSOX2, BMPR1B and PRRX2 as well as MYOM1 are candidates for selection on the Lundehund characteristics small body size, increased number of digits per paw and extraordinary mobility, respectively.     

Effective population size in the captive breeding program of the Lake Victoria cichlid Paralabidochromis chilotes

Microsatellite DNA markers were used to monitor levels of genetic variation in 3 generations of the American Zoo and Aquarium Association [AZA] Species Survival Plan [SSP] captive breeding program for the Lake Victoria cichlid Paralabidochromis chilotes. Temporal changes in the frequency of 15 alleles, across four polymorphic loci, were used to estimate effective population size (N_e). The upper limit of the 95% confidence interval for N_e never exceeded eight individuals, with all of the corresponding N_e/N ratios falling below 0.15. A test of the proportion of expected heterozygous individuals between the F₁ and F₃ generations indicated a significant decline in expected heterozygosity of 5% per generation. Alternative husbandry protocols, including subdividing the captive population, are addressed to reduce the further loss of genetic variation. Zoo Biol 18:215–222, 1999. © 1999 Wiley‐Liss, Inc.     

Genetic diversity and structure of an endemic and critically endangered stream river salamander (Caudata: Ambystoma leorae) in Mexico

Small or isolated populations are highly susceptible to stochastic events. They are prone and vulnerable to random demographic or environmental fluctuations that could lead to extinction due to the loss of alleles through genetic drift and increased inbreeding. We studied Ambystoma leorae an endemic and critically threatened species. We analyzed the genetic diversity and structure, effective population size, presence of bottlenecks and inbreeding coefficient of 96 individuals based on nine microsatellite loci. We found high levels of genetic diversity expressed as heterozygosity (H_o = 0.804, H_e = 0.613, H_e* = 0.626 and H_Nei = 0.622). The population presents few alleles (4–9 per locus) and genotypes (3–14 per locus) compared with other mole salamanders species. We identified three genetically differentiated subpopulations with a significant level of genetic structure (F_ST = 0.021, R_ST = 0.044 y D_est = 0.010, 95 % CI). We also detected a reduction signal in population size and evidence of a genetic bottleneck (M = 0.367). The effective population size is small (Ne = 45.2), but similar to another mole salamanders with restricted distributions or with recently fragmented habitat. The inbreeding coefficient levels detected are low (F_IS = ?0.619–0.102) as is gene flow. Despite, high levels of genetic diversity A. leorae is critically endangered because it is a small isolated population.     

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.     

Genetic diversity and molecular differentiation of Chinese toad based on microsatellite markers

Genetic diversity and population structure of 9 populations of Bufo gargarizans with total 111 samples in China were assessed using seven microsatellite loci. The analysed microsatellite markers produced 161 alleles, varied from 9 to 38 alleles each locus. The number of alleles per population per locus ranged from 4.43 to 10.29. Polymorphic information content showed that all seven loci were highly informative (mean = 0.810 ± 0.071). The average observed heterozygosity was less than the expected (0.353 ± 0.051 and 0.828 ± 0.067, respectively). All tested populations gave significant departures from Hardy–Weinberg equilibrium. Genetic differentiation among the populations was considerably high with the overall and pairwise F _ST values (mean = 0.160 ± 0.039), and showed fairly high level of inbreeding (indicated by a mean F _IS value of 0.504 ± 0.051) and global heterozygote deficit. In comparison to other amphibian studies; however, our results suggested that the level of genetic structuring in B. gargarizans was relatively low in the geographical scale of the study area. Interestingly, the speculated population bottleneck was found to be absent and the analyses provide only weak evidence for a recent contraction in size even though there was severe inbreeding (indicated by the F _IS value) in the Chinese toad populations.     

Effects of fragmentation phenomena on the genetic structure and gene flow in Centaurea cineraria group (Asteraceae) in the Mediterranean Basin

The complex history of the Mediterranean region illustrates how ancient and recent phenomena are closely associated with species distribution and the creation of phylogeographic divisions within Mediterranean flora. A good model to explore the genetic consequences of fragmentation can be found in Centaurea cineraria and its close relatives. We applied simple sequence repeat molecular markers to a dense population sampling throughout the distribution area of all C. cineraria taxa to study how fragmentation has altered the genetic structure and distribution of C. cineraria. The average gene diversity (H_e) was 0.286, and the average allelic richness (A_r) was 3.65 and ranged from 2.15 (C. gymnocarpa) to 5.25 (C. busambarensis). The F_IS averaged a relatively high 0.223, ranging from ? 0.724 in C. aeolica subsp. aeolica to 0.589 in C. leucadea. Our results indicate that habitat fragmentation over several generations reduced heterozygosity due to random genetic drift in populations of C. cineraria. This heterozygosity erosion becomes more severe when the inbreeding coefficient is positive and the outcrossing rates show a significant increase. The results observed for outcrossing rates and inbreeding coefficient could also indirectly support the possibility of disrupted gene flow or mating pattern changes in fragmented C. cineraria populations.     

Population dynamics of Walia ibex (Capra walie) at Simien Mountains National Park,Ethiopia

Intensive total direct counts of Walia ibex (Capra walie) population were performed at Simien Mountains National Park (SMNP) in 2009. Historical data were collected from SMNP and literature reviews. Different models were suited to determine population growth rates and intrinsic rate of increase. The population size estimated was 745 animals. The correlation between the two repeated counts was significant (r = 0.99 and P < 0.01). Mean instantaneous growth rate (r), growth rate per capita (λ) and population annual growth rate (Λ) were 2.6 ± 2.6, 0.03 ± 0.18 and 19.5 ± 50.4, respectively. Instantaneous growth rate and growth rate per capita were positively correlated (r = 0.958, P < 0.01). Average growth rate (rΛ) and intrinsic rate of increase (rr) under ideal (r = 0.950, P < 0.01) and random environments (r = 0.810, P < 0.01) were positively correlated. The population grows by 2.5% under ideal environments with an intrinsic increase of 0.04 (0.006%) and by 0.13% under random environments with intrinsic rate of decrease of ?0.184 or ?0.025% per year, respectively. The mean rank of the flock structure of whole population was 3.13, 3.88, 2.00 and 1.00 for males, females, juveniles and unidentified, respectively.     

Pedigree analysis of the Turkish Arab horse population: structure,inbreeding and genetic variability

The aim of this study was to evaluate genetic variability in the Turkish Arab horse population using pedigree information. This study is the first detailed pedigree analysis of the breed in Turkey. Pedigree data were collected from the National Studbook. The pedigree data for 23 668 horses, born between 1904 and 2014, were used in the analysis. From this data set, a reference population (RP) of 14 838 animals symbolising the last generation was defined. Demographic parameters, the inbreeding level (F), the average relatedness (AR), the effective population size (N_e), the effective number of founders (f_e), the effective number of ancestors (f_a) and the number of founder genome equivalents (f_g) were calculated for the population. The average generation interval for the RP was 12.2±4.6 years, whereas the calculated pedigree completeness levels were 98.2%, 96.6% and 95.0% for the first, second and third known generations. The mean equivalent generations (t), the average complete generations and the mean maximum generations for the RP were 7.8, 5.4 and 12.2, respectively, whereas the meanFand AR were 4.6% and 9.5% for the RP. The rate of inbred animals was 94.2% for the RP, whereas the number of founders, the number of ancestors and thef_e,f_aandf_gwere 342, 223, 40, 22 and 9.6 for the RP. The large differences observed betweenf_e, and the number of founders demonstrates that genetic diversity decreased between the founder and the RP. Contribution of the 14 most influential founder to the RP was 50.0%, whereas just eight ancestral horses can account for 50% of the genetic variability.N_eestimated via an individual increase in inbreeding per generation (${\overline{N}}_{\text{e}}$), and paired increase in coancestry$\left(\text{}{\overline{N}}_{eC}\right)$, were 74.4±3.9 and 73.5±0.58, respectively. The inbreeding increases with the pedigree knowledge. In addition, the decrease in inbreeding in last years is more noticeable.     

Effective population size is strongly correlated with breeding pond size in the endangered California tiger salamander, <Emphasis Type="Italic">Ambystoma californiense</Emphasis>

Maintaining genetic diversity and population viability in endangered and threatened species is a primary concern of conservation biology. Genetic diversity depends on population connectivity and effective population size (N_e), both of which are often compromised in endangered taxa. While the importance of population connectivity and gene flow has been well studied, investigating effective population sizes in natural systems has received far less attention. However, N_e plays a prominent role in the maintenance of genetic diversity, the prevention of inbreeding depression, and in determining the probability of population persistence. In this study, we examined the relationship between breeding pond characteristics and N_e in the endangered California tiger salamander, Ambystoma californiense. We sampled 203 individuals from 10 breeding ponds on a local landscape, and used 11 polymorphic microsatellite loci to quantify genetic structure, gene flow, and effective population sizes. We also measured the areas of each pond using satellite imagery and classified ponds as either hydrologically-modified perennial ponds or naturally occurring vernal pools, the latter of which constitute the natural breeding habitat for A. californiense. We found no correlation between pond area and heterozygosity or allelic diversity, but we identified a strong positive relationship between breeding pond area and N_e, particularly for vernal pools. Our results provide some of the first empirical evidence that variation in breeding habitat can be associated with differences in N_e and suggest that a more complete understanding of the environmental features that influence N_e is an important component of conservation genetics and management.     

Estimating effective population size using RADseq: Effects of SNP selection and sample size

Effective population size (N_e) is a key parameter of population genetics. However, N_e remains challenging to estimate for natural populations as several factors are likely to bias estimates. These factors include sampling design, sequencing method, and data filtering. One issue inherent to the restriction site‐associated DNA sequencing (RADseq) protocol is missing data and SNP selection criteria (e.g., minimum minor allele frequency, number of SNPs). To evaluate the potential impact of SNP selection criteria on N_e estimates (Linkage Disequilibrium method) we used RADseq data for a nonmodel species, the thornback ray. In this data set, the inbreeding coefficient F_IS was positively correlated with the amount of missing data, implying data were missing nonrandomly. The precision of N_eestimates decreased with the number of SNPs. Mean N_e estimates (averaged across 50 random data sets with2000 SNPs) ranged between 237 and 1784. Increasing the percentage of missing data from 25% to 50% increased N_e estimates between 82% and 120%, while increasing the minor allele frequency (MAF) threshold from 0.01 to 0.1 decreased estimates between 71% and 75%. Considering these effects is important when interpreting RADseq data‐derived estimates of effective population size in empirical studies.     

The use of agent-based modelling of genetics in conservation genetics studies

Animal Landscape and Man Simulation System a genetically explicit agent-based model was used to obtain measures for the genetic and demographic status of simulated populations. This investigation aimed to test the applicability of this approach for assessing the effect of environmental perturbations on populations’ temporal and spatial dynamics. This was achieved by assessing how three simple scenarios with increasing degree of environmental disturbance, simulated by populations bottlenecks repeated at different intervals, affected the genetic and demographic characteristics of the simulated population. Model outputs from a simplified landscape scenario concurred with theoretical expectations validating the model in a qualitative way. Differences in medians, means and coefficient of variation of the observed (H_o) and expected heterozygosity (H_e), population census size (N), effective population size (N_e), inbreeding coefficient (F) and N_e/N ratio were observed for simulated populations. Impacts occurred rapidly after simulated bottleneck events and genetic estimates were less variable, and therefore more reliable, than demographic estimates. Precise genetic consequences of the bottlenecks repeated at different intervals, and resulting population perturbations, are a complex balance between effects on population sub-structure, size and founding events. Agent-based models are appropriate tools to simulate these interactions, being sufficiently flexible to mimic real population processes under a range of environmental conditions. Such models incorporating explicit genetics provide a promising new approach to evaluate the impact of environmental changes on genetic composition of populations.     

Genetic diversity in adult and seedling populations of <Emphasis Type="Italic">Primula vulgaris</Emphasis> in a fragmented agricultural landscape

Habitat fragmentation is known to generally reduce the size of plant populations and increase their isolation, leading to genetic erosion and increased between-population genetic differentiation. In Flanders (northern Belgium) Primula vulgaris is very rare and declining. Populations have incurred strong fragmentation for the last decades and are now restricted to a few highly fragmented areas in an intensively used agricultural landscape. Previous studies showed that small populations of this long-lived perennial herb still maintained high levels of genetic variation and low genetic differentiation. This pattern can either indicate recent gene flow or represent historical variation. Therefore, we used polymorphic microsatellite loci to investigate genetic variation and structure in adult (which may still reflect historical variation) and seedling (recent generation, thus affected by current processes) life stages. The recent generation (seedlings) showed a significant loss of observed heterozygosity (H _o) together with lower expected heterozygosity (H _e), a trend for higher inbreeding levels (F _IS) and higher differentiation (F _ST) between populations compared to the adult generation. This might result from (1) a reduction in effective population size, (2) higher inbreeding levels in the seedlings, (3) a higher survival of heterozygotes over time due to a higher fitness of heterozygotes (heterosis) and/or a lower fitness of homozygotes (inbreeding depression), (4) overlapping generations in the adult life stage, or (5) a lack of establishment of new (inbred) adults from seedlings due to degraded habitat conditions. Combining restoration of both habitat quality and gene flow between populations may be indispensable to ensure a sustainable conservation of fragmented populations.     

Genetic characterization and population bottleneck in the Hucul horse based on microsatellite and mitochondrial data

The aim of this work was to gather information about the origin and genetic characterization of the Central European Hucul horse based on 71 horses using 17 microsatellites and the D‐loop region of mtDNA. Their genetic relationship to the Polish Konik (N = 7), German (N = 4) and Hungarian wild Przewalski horses (N = 4) and 200 horse sequences from GenBank was also analysed. Both microsatellite and mtDNA analysis showed a high genetic variation in the Hucul. A total of 130 alleles were detected, the mean number of observed alleles per microsatellite was 7.647, and the number of effective alleles was 4.401. The average observed and expected heterozygosity were 0.706 and 0.747, respectively. The high heterozygosity values and Wright's fixation index (F_IS) (?0.128) indicated a low level of inbreeding, low or no selection pressure, and large number of alleles. mtDNA analysis revealed 18 haplotypes for the Hucul population with a total of 23 variable sites. Haplotype and nucleotide diversities were 0.935 ± 0.011 and 0.022 ± 0.012, respectively. Neutrality tests (Tajima's D and Fu's Fs) were non‐significant, and mismatch distribution was ragged, indicating that the Hucul population is in genetic equilibrium. The most frequent mtDNA D‐loop region belonged to haplogroup A (48%), which was also present in Przewalski Wild horse samples, while Polish Konik samples belonged to three haplotypes and C, F, and G haplogroups. Large and significant pairwise Φ_ST values along with a small number of common haplotypes indicated a low level of gene flow and lack of genetic structure among the three studied breeds (Hucul, Konik, and Przewalski Wild horse). The present work contributes to our knowledge of the genetic diversity of the Hucul horse and helps to define its genetic conservation. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 54–65.     

Small effective population sizes of two remnant ocelot populations (<Emphasis Type="Italic">Leopardus pardalis albescens</Emphasis>) in the United States

Threatened populations are vulnerable to the effects of genetic drift and inbreeding, particularly when gene flow is low and the effective population size is small. Estimates of effective population size (N _e ) provide important information on the status of endangered populations that have experienced severe fragmentation and serve as indicators of genetic viability. Genetic data from microsatellite loci were used to estimate N _e for the 2 remaining populations of the endangered ocelot (Leopardus pardalis albescens) occurring in the United States. Several methods were used to calculate N _e , resulting in estimates ranging from N _e  = 8.0 (95% CI: 3.2–23.1) to 13.9 (95% CI: 7.7–25.1) for the population located at the Laguna Atascosa Wildlife Refuge in Cameron County, Texas. The ocelot population in Willacy County, Texas, had N _e estimates of 2.9 (95% CI: 1.7–5.6) and 3.1 (95% CI: 1.9–13.5), respectively. Estimates of N _e in both populations were below the critical value recommended for short-term viability.     

Inbreeding and extinction: Effects of rate of inbreeding

Deleterious alleles may be removed (purged) bynatural selection in populations undergoinginbreeding. However, there is controversyregarding the effectiveness of selection inreducing the risk of extinction due toinbreeding, especially in relation to the rateof inbreeding. We evaluated the effect of therate of inbreeding on reducing extinction risk,in populations of Drosophila melanogastermaintained using full-sib mating (160replicates), or at effective population sizes(N _e) of 10 (80) or 20 (80).Extinction rates in the populations maintainedusing full-sib mating occurred at lower levelsof inbreeding than in the larger populations,whereas the two larger populations did notdiffer significantly from each other.Inbreeding coefficients at 50% extinction were0.62, 0.79 and 0.77 for the full-sib (N _e = 2.6), N _e = 10 and N _e = 20 treatments, respectively. Populations of N _e = 20 that remained extant after 60 generations, showed inbreeding depression, with the mean fitness of these populations being only 45% of the outbredcontrols. There was considerable variationamong the 31 inbred populations in fitness, butnone of the N _e = 20 populations hadfitness that was higher than the outbredcontrol. We conclude that purging may slow therate of extinction slightly, but it cannot berelied on to eliminate the deleterious effectsof inbreeding.     

Genetic variation in low elevation Douglas-fir of British Columbia and its relevance to gene conservation

Patterns of variation were studied at 20 isozyme loci in 49 coastal, low-elevation Douglas-fir populations in SW British Columbia and NW Washington State. Several components of variation were estimated for each population including the number of alleles per locus N _a, number of alleles per polymorphic locus N _a(95), inbreeding F, heterozygosity H, and population divergence D. F was near zero indicating nearly complete outcrossing within populations. H was quite high (16%) and in aecord with previous studies of Douglas-fir. D values were low (equivalent to Wrights F _ST of 0.08) indicating levels of gene flow sufficient to largely homogenize populations. The parameters of diversity N _a, N _a(95), H, and D showed little intercorrelation across populations. A homogenous pattern of genetic relationship among populations was shown by the clustering of populations based on their inferred relationship, and by the principal components of the matrix of inferred genetic relationship. Because of the complex nature of gene diversity and the continuous nature of population differentiation in Douglas-fir, it is difficult with isozyme markers to identify specific populations of value for genetic conservation in this species.This paper is dedicated to J.C. Heaman on the occasion of his retirement.     

Interspecific Genetic Variability of the Oyster Mushroom Pleurotus ostreatusAs Revealed by Allozyme Gene Analysis

Allozyme variation in natural populations of basidiomycete fungus Pleurotus ostreatus (88 individuals) from three regions of central Russia was studied. The species was shown to have 92.86% of polymorphic allozyme loci and expected heterozygosity H _e = 0.49. The mean number of alleles per locus was 3.5. The genetic differences among populations were supported by F-statistics (F _ST = 0.750). The low level of inbreeding (F _IS = 0.018) suggests that the P. ostreatus populations are panmictic, and the main reproduction mode involves basidiospores dispersing at long distances. Using cluster analysis, geographically isolated populations and intersterile groups were differentiated within the complex P. ostreatus species.     

INCREASED PROBABILITY OF EXTINCTION DUE TO DECREASED GENETIC EFFECTIVE POPULATION SIZE: EXPERIMENTAL POPULATIONS OF CLARKIA PULCHELLA

We established replicated experimental populations of the annual plant Clarkia pulchella to evaluate the existence of a causal relationship between loss of genetic variation and population survival probability. Two treatments differing in the relatedness of the founders, and thus in the genetic effective population size (N_e), were maintained as isolated populations in a natural environment. After three generations, the low N_e treatment had significantly lower germination and survival rates than did the high N_e treatment. These lower germination and survival rates led to decreased mean fitness in the low N_e populations: estimated mean fitness in the low N_e populations was only 21% of the estimated mean fitness in the high N_e populations. This inbreeding depression led to a reduction in population survival: at the conclusion of the experiment, 75% of the high N_e populations were still extant, whereas only 31% of the low N_e populations had survived. Decreased genetic effective population size, which leads to both inbreeding and the loss of alleles by genetic drift, increased the probability of population extinction over that expected from demographic and environmental stochasticity alone. This demonstrates that the genetic effective population size can strongly affect the probability of population persistence.     

RECENT EVOLUTION AND DIVERGENCE AMONG POPULATIONS OF A RARE MEXICAN ENDEMIC,CHIHUAHUA SPRUCE,FOLLOWING HOLOCENE CLIMATIC WARMING

Fragmentation and reduction in population size are expected to reduce genetic diversity. However, examples from natural populations of forest trees are scarce. The range of Chihuahua spruce retreated northward and fragmented coincident with the warming climate that marked the early Holocene. The isolated populations vary from 15 to 2441 trees, which provided an opportunity to test whether census number is a good predictor of genetic diversity. Mean expected heterozygosity, H_e, based on 24 loci in 16 enzyme systems, was 0.093 for 10 sampled populations, which is within the range reported for conifers. However, estimates varied more than twofold among populations and H_e was closely related to the logarithm of the number of mature trees in the population (r_He,N = 0.93). Diversity among populations, F_ST, was 24.8% of the total diversity, which is higher than that observed in almost all conifer species studied. Nei's genetic distance, D, was not related to geographic distance between populations, and D? was 0.033, which is higher than estimates for most wide-ranging species. Most populations had excess homozygosity and the fixation index, F_IS, was higher than that reported for all but one species of conifer. Nm, the number of migrants per generation, was 0.43 to 0.76, depending on estimation procedure, and is the smallest observed in conifers. The data suggest that populations of Chihuahua spruce have differentiated by drift and that they are effectively isolated. The results illustrate how a combination of paleontological observation and molecular markers can be used to illuminate recent evolutionary events. Multilocus estimates of outcrossing for two small populations were zero (complete selfing) and 0.153, respectively, which are in striking contrast to the near complete outcrossing observed in most conifers. The high fixation index and a high proportion of empty seeds (45%) suggest that inbreeding may be a serious problem for conservation of Chihuahua spruce.     

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.