Microsatellite variation, population structure, and bottlenecks in the threatened copperbelly water snake

Habitat loss and isolation is pervasive in the Midwest U.S. Wetlands are experiencing particularly dramatic declines, yet there is a paucity of information on the genetic impacts of these losses to obligate wetland vertebrates. We quantified the genetic variation of extant populations of a shallow wetland specialist and evaluated potential reductions in population size (i.e. bottlenecks) using seven polymorphic microsatellite DNA markers. We analyzed 228 copperbelly water snakes (Nerodia erythrogaster neglecta), representing populations from three states. Moderate genetic differentiation exists among all three regions (F _ST = 0.12, P < 0.001), with evidence for low levels of differentiation within the federally protected Ohio region (F _ST = 0.025, P = 0.007), and moderate to strong differentiation within the Indiana region (F _ST = 0.23, P < 0.001). Furthermore, Bayesian clustering (i.e. STRUCTURE) supports the separation of the Indiana sites, both from each other and from all other sampling sites. However, it does not support the separation of the Ohio sites from the Kentucky sites. Differentiation among sampling sites did not appear to be related to geographic distance, but rather depended on the quality of terrestrial corridors used for dispersal. Mode shifts in allele frequencies and excess heterozygosity tests were negative, while M-ratio tests were nearly all positive, indicating the likelihood of historical rather than contemporary population bottlenecks. However, potential subspecific intergradation in the Kentucky region may have artificially lowered the M-ratio, and we suggest caution when using the M-ratio approach if intergradation is suspected. Our results have conservation implications for wetland management and management of the copperbelly populations, and emphasizes the importance of protecting wetland complexes.     

Assessing population genetic structure via the maximisation of genetic distance

Background

The inference of the hidden structure of a population is an essential issue in population genetics. Recently, several methods have been proposed to infer population structure in population genetics.

Methods

In this study, a new method to infer the number of clusters and to assign individuals to the inferred populations is proposed. This approach does not make any assumption on Hardy-Weinberg and linkage equilibrium. The implemented criterion is the maximisation (via a simulated annealing algorithm) of the averaged genetic distance between a predefined number of clusters. The performance of this method is compared with two Bayesian approaches: STRUCTURE and BAPS, using simulated data and also a real human data set.

Results

The simulations show that with a reduced number of markers, BAPS overestimates the number of clusters and presents a reduced proportion of correct groupings. The accuracy of the new method is approximately the same as for STRUCTURE. Also, in Hardy-Weinberg and linkage disequilibrium cases, BAPS performs incorrectly. In these situations, STRUCTURE and the new method show an equivalent behaviour with respect to the number of inferred clusters, although the proportion of correct groupings is slightly better with the new method. Re-establishing equilibrium with the randomisation procedures improves the precision of the Bayesian approaches. All methods have a good precision for F_ST ≥ 0.03, but only STRUCTURE estimates the correct number of clusters for F_ST as low as 0.01. In situations with a high number of clusters or a more complex population structure, MGD performs better than STRUCTURE and BAPS. The results for a human data set analysed with the new method are congruent with the geographical regions previously found.

Conclusion

This new method used to infer the hidden structure in a population, based on the maximisation of the genetic distance and not taking into consideration any assumption about Hardy-Weinberg and linkage equilibrium, performs well under different simulated scenarios and with real data. Therefore, it could be a useful tool to determine genetically homogeneous groups, especially in those situations where the number of clusters is high, with complex population structure and where Hardy-Weinberg and/or linkage equilibrium are present.     

Genetic differentiation of <Emphasis Type="Italic">Sorbus torminalis</Emphasis> in Eastern Europe as determined by microsatellite markers

The genetic variation in fourteen Sorbus torminalis (L.) Crantz. populations distributed over the eastern and south-eastern part of its range was studied using seven nuclear microsatellite loci. The differentiation level was relatively high (F _ST = 0.228), as expected for a species with a fragmented range. The distance-based approach to the analysis of differentiation patterns (neighbour-joining tree based on pairwise coefficients of differentiation) did not reveal a clear geographical structure. On the other hand, model-based Bayesian methods (BAPS and STRUCTURE) gave geographically continuous clusters of populations. The occurrence of populations deviating strongly from the general pattern is attributed to founder effect. In spite of a generally high differentiation, a significant isolation-by-distance pattern was found, which might be a consequence of postglacial migration and gene flow among descendants of different refugia.     

Genetic and phenotypic variation in central and northern European populations of Aedes (Aedimorphus) vexans (Meigen, 1830) (Diptera,Culicidae)

The floodwater mosquito Aedes vexans can be a massive nuisance in the flood plain areas of mainland Europe, and is the vector of Tahyna virus and a potential vector of Dirofilaria immitis. This epidemiologically important species forms three subspecies worldwide, of which Aedes vexans arabiensis has a wide distribution in Europe and Africa. We quantified the genetic and phenotypic variation in Ae. vexans arabiensis in populations from Sweden (northern Europe), Hungary, and Serbia (central Europe). A landscape genetics approach (F_ST, STRUCTURE, BAPS, GENELAND) revealed significant differentiation between northern and southern populations. Similar to genetic data, wing geometric morphometrics revealed two different clusters, one made by Swedish populations, while another included Hungarian and Serbian populations. Moreover, integrated genetic and morphometric data from the spatial analysis suggested groupings of populations into three clusters, one of which was from Swedish and Hungarian populations. Data on spatial analysis regarding an intermediate status of the Hungarian population was supported by observed Isolation‐by‐Distance patterns. Furthermore, a low proportion of interpopulation vs intrapopulation variance revealed by AMOVA and low‐to‐moderate F_ST values on a broader geographical scale indicate a continuous between‐population exchange of individuals, including considerable gene flow on the regional scale, are likely to be responsible for the maintenance of the observed population similarity in Aе. vexans. We discussed data considering population structure in the light of vector control strategies of the mosquito from public health importance.     

Low levels of quantitative and molecular genetic differentiation among natural populations of Medicago ciliaris Kroch. (Fabaceae) of different Tunisian eco-geographical origin

In this paper, we analyze the genetic variability in four Tunisian natural populations of Medicago ciliaris using 19 quantitative traits and six polymorphic microsatellite loci. We investigated the amplification transferability of 30 microsatellites developed in the model legume M. truncatula to M. ciliaris. Results revealed that about 56.66% of analyzed markers are valuable genetic markers for M. ciliaris. The most genetic diversity at quantitative traits and microsatellite loci was found to occur within populations (>80%). Low differentiations among populations at quantitative traits Q _ST  = 0.146 and molecular markers F _ST  = 0.18 were found. The majority of measured traits exhibited no significant difference in the level of Q _ST and F _ST . Furthermore, significant correlations established between these traits and eco-geographical factors suggested that natural selection should be invoked to explain the level of phenotypic divergence among populations rather than drift. There was no significant correlation between population differentiation at quantitative traits and molecular markers. Significant spatial genetic structure consistent with models of isolation by distance was detected within all studied populations. The site-of-origin environmental factors explain about 9.07% of total phenotypic genetic variation among populations. The eco-geographical factors that influence more the variation of measured traits among populations are the soil texture and altitude. Nevertheless, there were no consistent pattern of associations between gene diversity (He) and environmental factors.     

Morphological and microsatellite diversity associated with ecological factors in natural populations of Medicago laciniata Mill. (Fabaceae)

Genetic variability in 10 natural Tunisian populations of Medicago laciniata were analysed using 19 quantitative traits and 12 polymorphic microsatellite loci. A large degree of genetic variability within-populations and among-populations was detected for both quantitative characters and molecular markers. High genetic differentiation among populations for quantitative traits was seen, with Q _ST = 0.47, and F _ST = 0.47 for microsatellite markers. Several quantitative traits displayed no statistical difference in the levels of Q _ST and F _ST . Further, significant correlations between quantitative traits and eco-geographical factors suggest that divergence in the traits among populations may track environmental differences. There was no significant correlation between genetic variability at quantitative traits and microsatellite markers within populations. The site-of-origin of eco-geographical factors explain between 18.13% and 23.40% of genetic variance among populations at quantitative traits and microsatellite markers, respectively. The environmental factors that most influence variation in measured traits among populations are assimilated phosphorus (P₂0₅) and mean annual rainfall, followed by climate and soil texture, altitude and organic matter. Significant associations between eco-geographical factors and gene diversity, H _e , were established in five-microsatellite loci suggesting that these simple sequence repeats (SSRs) are not necessarily biologically neutral.     

On the Apportionment of Population Structure

Measures of population differentiation, such as F_ST, are traditionally derived from the partition of diversity within and between populations. However, the emergence of population clusters from multilocus analysis is a function of genetic structure (departures from panmixia) rather than of diversity. If the populations are close to panmixia, slight differences between the mean pairwise distance within and between populations (low F_ST) can manifest as strong separation between the populations, thus population clusters are often evident even when the vast majority of diversity is partitioned within populations rather than between them. For any given F_ST value, clusters can be tighter (more panmictic) or looser (more stratified), and in this respect higher F_ST does not always imply stronger differentiation. In this study we propose a measure for the partition of structure, denoted E_ST, which is more consistent with results from clustering schemes. Crucially, our measure is based on a statistic of the data that is a good measure of internal structure, mimicking the information extracted by unsupervised clustering or dimensionality reduction schemes. To assess the utility of our metric, we ranked various human (HGDP) population pairs based on F_ST and E_ST and found substantial differences in ranking order. E_ST ranking seems more consistent with population clustering and classification and possibly with geographic distance between populations. Thus, E_ST may at times outperform F_ST in identifying evolutionary significant differentiation.     

Genomic diversity and differentiation of a managed island wild boar population

The evolution of island populations in natural systems is driven by local adaptation and genetic drift. However, evolutionary pathways may be altered by humans in several ways. The wild boar (WB) (Sus scrofa) is an iconic game species occurring in several islands, where it has been strongly managed since prehistoric times. We examined genomic diversity at 49 803 single-nucleotide polymorphisms in 99 Sardinian WBs and compared them with 196 wild specimens from mainland Europe and 105 domestic pigs (DP; 11 breeds). High levels of genetic variation were observed in Sardinia (80.9% of the total number of polymorphisms), which can be only in part associated to recent genetic introgression. Both Principal Component Analysis and Bayesian clustering approach revealed that the Sardinian WB population is highly differentiated from the other European populations (F_ST=0.126–0.138), and from DP (F_ST=0.169). Such evidences were mostly unaffected by an uneven sample size, although clustering results in reference populations changed when the number of individuals was standardized. Runs of homozygosity (ROHs) pattern and distribution in Sardinian WB are consistent with a past expansion following a bottleneck (small ROHs) and recent population substructuring (highly homozygous individuals). The observed effect of a non-random selection of Sardinian individuals on diversity, F_ST and ROH estimates, stressed the importance of sampling design in the study of structured or introgressed populations. Our results support the heterogeneity and distinctiveness of the Sardinian population and prompt further investigations on its origins and conservation status.     

Genetic diversity and gene flow in a Caribbean tree <Emphasis Type="Italic">Pterocarpus officinalis</Emphasis> Jacq.: a study based on chloroplast and nuclear microsatellites

We analysed the molecular diversity of Pterocarpus officinalis, a tree species distributed in Caribbean islands, South and Central America to quantify the genetic variation within island, to assess the pattern of differentiation and infer levels of gene flow; with the overall goal of defining a strategy of conservation. Two hundred two individuals of 9 populations were analysed using three chloroplast and six nuclear microsatellite markers. The observed heterozygosity varied markedly among the populations for nuclear (H _Onuc = 0.20–0.50) and chloroplast microsatellites (H _cp = 0.22–0.68). The continental population from French Guyana showed a higher value of H _Onuc than island populations, and the differences were significant in some cases. The fixation index F _IS ranged from −0.043 to 0.368; a significant heterozygote deficit was detected in 7 populations. The heterozygosity excess method suggested that two populations in Guadeloupe have undergone a recent bottleneck. Global and pairwise F _ST were high for both nuclear (F _STnuc = 0.29) and chloroplast microsatellites (F _STcp = 0.58). The neighbour-joining tree based on both markers, presented a differentiation pattern that can be explained by the seed dispersal by flotation and marine stream. The comparison of Bayesian approach and the method based on allelic frequency demonstrate a very limited number of migrants between populations.     

基于SSR标记探讨三种金花茶植物的遗传多样性和遗传结构

薄叶金花茶、小花金花茶和小瓣金花茶是三种濒危金花茶植物,为了解珍稀濒危植物遗传多样性和遗传结构,该研究利用微卫星标记对他们的7个种群共184个个体进行了遗传多样性和遗传结构分析。结果表明:11个位点共检测到等位基因92个。在物种水平上,小瓣金花茶平均等位基因数(N_A)为3.9、有效等位基因数(N_E)为2.328、观测杂合度(H_o)为0.520、期望杂合度(H_e)为0.501,高于薄叶金花茶和小花金花茶。在种群水平上,有效等位基因数(N_E)在1.788～2.466之间,期望杂合度(H_e)在0.379～0.543之间;种群间遗传分化系数(FST)在0.143 7～0.453 3之间,种群间基因流(N_m)在0.301 5～1.488 9之间。AMOVA分子变异分析显示65.72%的变异存在于种群内。三种金花茶具有较低水平的遗传多样性和高水平的种群间遗传分化。STRUCTURE和PCoA种群遗传结构分析结果将取样种群分为2组,即薄叶金花茶和小花金花茶大部分个体分为一组,小瓣金花茶大部分个体分为一组。现存所有种群应根据实际情况尽快采取就地保护或迁地保护措施。     

AFLP markers suggest low population genetic differentiation of the black rockfish Sebastes schlegelii

AFLP markers were used to investigate the population genetic differentiation among populations of Sebastes schlegelii from China and Japan. A total of 405 bands were amplified for 180 individuals from 9 populations by 3 pairs of selective primers. S. schlegelii exhibited high Nei's gene diversity with an average value of 0.29 for all populations. No significant genealogical branches or clusters corresponding to sampling localities were detected by UPGMA tree. The results of AMOVA analysis and pairwise F_ST values showed no genetic divergence among different geographic populations. The result of STRUCTURE demonstrated all populations of S. schlegelii examined in the study shared one single gene pool.     

Patterns of nuclear and mitochondrial DNA variation in Iberian populations of <Emphasis Type="Italic">Emys orbicularis</Emphasis> (Emydidae): conservation implications

The European pond turtle (Emys orbicularis) is threatened and in decline in several regions of its natural range, due to habitat loss combined with population fragmentation. In this work, we have focused our efforts on studying the genetic diversity and structure of Iberian populations with a fine-scale sampling (254 turtles in 10 populations) and a representation from North Africa and Balearic island populations. Using both nuclear and mitochondrial markers (seven microsatellites, ∼1048 bp nDNA and ∼1500 bp mtDNA) we have carried out phylogenetic and demographic analyses. Our results show low values of genetic diversity at the mitochondrial level although our microsatellite dataset revealed relatively high levels of genetic variability with a latitudinal genetic trend decreasing from southern to northern populations. A moderate degree of genetic differentiation was estimated for Iberian populations (genetic distances, F _ST values and clusters in the Bayesian analysis). The results in this study combining mtDNA and nDNA, provide the most comprehensive population genetic data for E. orbicularis in the Iberian Peninsula. Our results suggest that Iberian populations within the Iberian–Moroccan lineage should be considered as a single subspecies with five management units, and emphasize the importance of habitat management rather than population reinforcement (i.e. captive breeding and reintroduction) in this long-lived species.     

Comparisons of microsatellite variability and population genetic structure of two endangered wild rice species, <Emphasis Type="Italic">Oryza rufipogon</Emphasis> and <Emphasis Type="Italic">O. officinalis</Emphasis>, and their conservation implications

Conserving endangered wild rice species requires a thorough understanding of their population genetic structure and appropriate approaches. We applied six and seven microsatellite loci to study the genetic structure of six populations throughout the range of Chinese Oryza rufipogon and Oryza officinalis, respectively. The results showed that O. rufipogon possesses higher levels of genetic diversity but lower differentiation (R_S = 3.2713, P = 100.0%, H_O = 0.1401, H_S = 0.5800, F_ST = 0.271) than O. officinalis (R_S = 2.0545, P = 57.14%, H_O = 0.0470, H_S = 0.2830, F_ST = 0.554). Mean population F_IS was slightly larger for O. officinalis (F_IS = 0.844) than that for O. rufipogon (F_IS = 0.755), indicating that O. officinalis has slightly higher departures from Hardy–Weinberg expectations and heterozygosity deficits than O. rufipogon. In addition to different origins and evolutionary histories, O. officinalis has restricted gene flow, high inbreeding, isolated small populations and fewer opportunities of hybridization with other taxa, which may determine major differences in population genetic structure from O. rufipogon. Our results suggest the adoption of a plan of involving fewer populations but more individuals within populations for O. rufipogon, while both the number of populations and the individuals for a sampled population should be almost equally considered for O. officinalis. The known high degree of inbreeding in the populations of both species implies that conservation and restoration genetics should particularly focus on the maintenance of historically significant processes such as high levels of outbreeding, gene flow and large effective population sizes. We finally proposed to further estimate the role of rice gene flow in the conservation of O. rufipogon, and to perform detailed analysis of mating systems in both species for better conservation perspectives of their ecological and evolutionary processes.     

Molecular and quantitative trait variation within and among small fragmented populations of the endangered plant species Psilopeganum sinense

Background and Aims

Natural selection and genetic drift are important evolutionary forces in determining genetic and phenotypic differentiation in plant populations. The extent to which these two distinct evolutionary forces affect locally adaptive quantitative traits has been well studied in common plant and animal species. However, we know less about how quantitative traits respond to selection pressures and drift in endangered species that have small population sizes and fragmented distributions. To address this question, this study assessed the relative strengths of selection and genetic drift in shaping population differentiation of phenotypic traits in Psilopeganum sinense, a naturally rare and recently endangered plant species.

Methods

Population differentiation at five quantitative traits (Q_ST) obtained from a common garden experiment was compared with differentiation at putatively neutral microsatellite markers (F_ST) in seven populations of P. sinense. Q_ST estimates were derived using a Bayesian hierarchical variance component method.

Key Results

Trait-specific Q_ST values were equal to or lower than F_ST. Neutral genetic diversity was not correlated with quantitative genetic variation within the populations of P. sinense.

Conclusions

Despite the prevalent empirical evidence for Q_ST > F_ST, the results instead suggest a definitive role of stabilizing selection and drift leading to phenotypic differentiation among small populations. Three traits exhibited a significantly lower Q_ST relative to F_ST, suggesting that populations of P. sinense might have experienced stabilizing selection for the same optimal phenotypes despite large geographical distances between populations and habitat fragmentation. For the other two traits, Q_ST estimates were of the same magnitude as F_ST, indicating that divergence in these traits could have been achieved by genetic drift alone. The lack of correlation between molecular marker and quantitative genetic variation suggests that sophisticated considerations are required for the inference of conservation measures of P. sinense from neutral genetic markers.     

Genotype frequency and F ST analysis of polymorphisms in immunoregulatory genes in Chinese and Caucasian populations

Selection and genetic drift can create genetic differences between populations. Cytokines and chemokines play an important role in both hematopoietic development and the inflammatory response. We compared the genotype frequencies of 45 SNPs in 30 cytokine and chemokine genes in two healthy Chinese populations and one Caucasian population. Several SNPs in IL4 had substantial genetic differentiation between the Chinese and Caucasian populations (F _ST ~0.40), and displayed a strikingly different haplotype distribution. To further characterize common genetic variation in worldwide populations at the IL4 locus, we genotyped 9 SNPs at the IL4 gene in the Human Diversity Panel’s (N = 1056) individuals from 52 world geographic regions. We observed low haplotype diversity, yet strikingly different haplotype frequencies between non-African populations, which may indicate different selective pressures on the IL4 gene in different parts of the world. SNPs in CSF2, IL6, IL10, CTLA4, and CX3CR1 showed moderate genetic differentiation between the Chinese and Caucasian populations (0.15 < F _ST < 0.25). These results suggest that there is substantial genetic diversity in immune genes and exploration of SNP associations with immune-related diseases that vary in incidence across these two populations may be warranted.     

Genetic diversity and structure among Iranian tall fescue populations based on genomic-SSR and EST-SSR marker analysis

Tall fescue (Festuca arundinacea Schreb. subsp. arundinacea) is one of the most economically important forage grasses in cold and temperate regions of the world. In this study, we have assessed the genetic diversity and structure of wild Iranian tall fescue populations. Thirty-seven individuals from nine natural populations from northern, western, and southern Iranian mountain ranges were analyzed using six genomic-SSRs and five EST-SSRs primer pairs. Our analysis has demonstrated that transcribed EST-SSR regions showed levels of polymorphism similar to genomic-SSR regions. UPGMA, repeated bisection, and model-based Bayesian STRUCTURE clustering methods coupled with neighbor-net network were used to establish six divergent groups of individuals. F _ST estimates among clusters showed moderate to low genetic structure. Within-group genetic diversity estimate H and partial correlations between genetic and geographic distances among populations suggested that western Zagros population was related to the rest of the Iranian individuals. The isolation-by-distance hypothesis was not supported by SSR data and the present geographical sampling.     

Comparative population structure of Chinese sumac aphid <Emphasis Type="Italic">Schlechtendalia chinensis</Emphasis> and its primary host-plant <Emphasis Type="Italic">Rhus chinensis</Emphasis>

Most of our current understanding of comparative population structure has been come from studies of parasite–host systems, whereas the genetic comparison of gallnut-aphids and their host-plants remain poorly documented. Here, we examined the population genetic structure of the Chinese sumac aphid Schlechtendalia chinensis and its unique primary host-plant Rhus chinensis in a mountainous province in western China using inter-simple sequence repeat (ISSR) markers. Despite being sampled from a mountainous geographic range, analysis of molecular variance (AMOVA) showed that the majority of genetic variation occurred among individuals within populations of both the aphid and its host. The aphid populations were found to be structured similarly to their primary host populations (F _ST values were 0.239 for the aphid and 0.209 for its host), suggesting that there are similar patterns of gene flow between the populations of the aphid and between populations of its host-plant. The genetic distances (F _ST/1 − F _ST) between the aphid populations and between its host-plant populations were uncorrelated, indicating that sites with genetically similar host-plant populations may not always have genetically similar aphid populations. The lack of relationships between genetic and geographical distance matrices suggested that isolation by distance (IBD) played a negligible role at this level. This may be mainly attributed to the founder effect, genetic drift and the relative small spatial scale between populations. Zhumei Ren and Bin Zhu contributed equally to this work.     

Genetic structure of Pyrenophora teres net and spot populations as revealed by microsatellite analysis

The population structure of the fungal pathogen Pyrenophora teres, collected mainly from different regions of the Czech and Slovak Republics, was examined using a microsatellite analyses (SSR). Among 305 P. teres f. teres (PTT) and 82 P. teres f. maculata (PTM) isolates that were collected, the overall gene diversity was similar (? = 0.12 and ? = 0.13, respectively). A high level of genetic differentiation (F_ST = 0.46; P < 0.001) indicated the existence of population structure. Nine clusters that were found using a Bayesian approach represent the genetic structure of the studied P. teres populations. Two clusters consisted of PTM populations; PTT populations formed another seven clusters. An exact test of population differentiation confirmed the results that were generated by Structure. There was no difference between naturally infected populations over time, and genetic distance did not correlate with geographical distance. The facts that all individuals had unique multilocus genotypes and that the hypothesis of random mating could not be rejected in several populations or subpopulations serve as evidence that a mixed mating system plays a role in the P. teres life cycle. Despite the fact that the genetic differentiation value between PTT and PTM (F_ST = 0.30; P < 0.001) is lower than it is between the populations within each form (F_ST = 0.40 (PTT); F_ST = 0.35 (PTM); P < 0.001) and that individuals with mixed PTT and PTM genomes were found, the two forms of P. teres form genetically separate populations. Therefore, it can be assumed that these populations have most likely undergone speciation.     

Microsatellite variation and differentiation among local populations of Castanopsis species in Japan

Microsatellite variations in Castanopsis species in Japan were examined to clarify the genetic relationships among 25 local populations according to the difference in the number of layers of adaxial epidermis in the leaves. Six microsatellite loci were assayed for 629 seedlings from the populations, and these seedlings were classified into five types according to the state of the leaf epidermis. Remarkable differences in the allele frequency of the six microsatellite loci were observed among these local populations. The coefficients of genetic differentiation, R_ST, of each locus ranged from 0.209 to 0.388. An unweighted pair-group method (UPGMA) phenogram constructed on the population pairwise R_ST over the loci revealed three clusters (A–C), and six sub-clusters. These clusters reflected the differences in the occurrence frequency of seedlings in each epidermis type within a population. Our findings suggest that clusters A and C are the local populations dominated by Castanopsis sieboldii and Castanopsis cuspidata, respectively, while local populations of cluster B are composed of the two Castanopsis species and/or include many individuals derived by hybridization. The six sub-clusters were found to reflect the geographic relationship among the populations, suggesting a different process for geographic population dynamics during the postglacial period.     

Analysis of F ST outliers at allozyme loci in Pacific salmon: implications for natural selection

Natural selection has been invoked to explain the observed geographic distribution of allozyme allele frequencies for a number of teleost species. The effects of selection on allozyme loci in three species of Pacific salmon were tested. A simulation-based approach to estimate the null distribution of population differentiation (F _ST) and test for F _ST outliers was used. This approach showed that a majority of allozyme loci conform to neutral expectations predicted by the simulation model, with relatively few F _ST outliers found. No consistent F _ST outlier loci were found across species. Analysis of population sub-groups based on geography and genetic identity reduced the number of outlier loci for some species, indicating that large geographic groups may include genetically divergent populations and/or that there is geographic heterogeneity in selection pressure upon allozyme loci. Two outlier allozyme loci found in this analysis, lactate dehydrogenase-B and malic enzyme, have been shown to be influenced by selection in other teleost species. This approach is also useful in identifying allozyme loci (or other genetic markers) that meet assumptions for population genetic study.