Population-Specific Links Between Heterozygosity and the Rate Human Microsatellite Evolution

Microsatellites form an abundant class of DNA sequences used widely as genetic markers. Surprisingly, the length of human microsatellites varies highly predictably with distance from Africa, apparently following the linear decline in variability that arose as we colonised the world. Such patterns have been used to argue that heterozygosity modulates the rate of microsatellite evolution. Here I test the ensuing prediction that variation in demographic history will cause individual populations predictably either to lead or to lag any given trend in length. I find that they do: larger populations with locally higher heterozygosity have microsatellites that are longer when a locus is expanding and shorter when a locus is contracting. These patterns remain even after controlling for the stepwise way in which heterozygosity and allele lengths decline across the world. This analysis provides support for a strongly discontinuous model for how human genetic variability is distributed and shows how individual populations differ in the average rate their microsatellites are evolving. Such patterns have the potential to provide a new window onto historical demography.     

Growth rate correlates to individual heterozygosity in the European eel, Anguilla anguilla L

Heterozygosity-fitness correlations (HFCs) have been reported in populations of many species. We provide evidence for a positive correlation between genetic variability and growth rate at 12 allozyme loci in a catadromous marine fish species, the European eel (Anguilla anguilla L.). More heterozygous individuals show a significantly higher length and weight increase and an above average condition index in comparison with more homozygous individuals. To a lesser extent, six microsatellite loci show a similar pattern, with positive but not significant correlations between heterozygosity and growth rate. The HFCs observed could be explained by an effect of either direct allozyme over-dominance or associative overdominance. Selection affecting some of the allozyme loci would explain the greater strength of the HFCs found at allozymes in comparison with microsatellites and the lack of correlation between MLH at allozymes and MLH at microsatellites. Associative overdominance (where allozyme loci are merely acting as neutral markers of closely linked fitness loci) might provide an explanation for the HFCs if we consider that allozyme loci have a higher chance than microsatellites to be in linkage disequilibrium with fitness loci.     

Comparing RADseq and microsatellites for estimating genetic diversity and relatedness — Implications for brown trout conservation

The conservation and management of endangered species requires information on their genetic diversity, relatedness and population structure. The main genetic markers applied for these questions are microsatellites and single nucleotide polymorphisms (SNPs), the latter of which remain the more resource demanding approach in most cases. Here, we compare the performance of two approaches, SNPs obtained by restriction‐site‐associated DNA sequencing (RADseq) and 16 DNA microsatellite loci, for estimating genetic diversity, relatedness and genetic differentiation of three, small, geographically close wild brown trout (Salmo trutta) populations and a regionally used hatchery strain. The genetic differentiation, quantified as F_ST, was similar when measured using 16 microsatellites and 4,876 SNPs. Based on both marker types, each brown trout population represented a distinct gene pool with a low level of interbreeding. Analysis of SNPs identified half‐ and full‐siblings with a higher probability than the analysis based on microsatellites, and SNPs outperformed microsatellites in estimating individual‐level multilocus heterozygosity. Overall, the results indicated that moderately polymorphic microsatellites and SNPs from RADseq agreed on estimates of population genetic structure in moderately diverged, small populations, but RADseq outperformed microsatellites for applications that required individual‐level genotype information, such as quantifying relatedness and individual‐level heterozygosity. The results can be applied to other small populations with low or moderate levels of genetic diversity.     

Patterns, sources and ecological implications of clonal diversity in apomictic Ranunculus carpaticola (Ranunculus auricomus complex, Ranunculaceae)

Sources and implications of genetic diversity in agamic complexes are still under debate. Population studies (amplified fragment length polymorphisms, microsatellites) and karyological methods (Feulgen DNA image densitometry and flow cytometry) were employed for characterization of genetic diversity and ploidy levels of 10 populations of Ranunculus carpaticola in central Slovakia. Whereas two diploid populations showed high levels of genetic diversity, as expected for sexual reproduction, eight populations are hexaploid and harbour lower degrees of genotypic variation, but maintain high levels of heterozygosity at many loci, as is typical for apomicts. Polyploid populations consist either of a single AFLP genotype or of one dominant and a few deviating genotypes. genotype/genodive and character incompatibility analyses suggest that genotypic variation within apomictic populations is caused by mutations, but in one population probably also by recombination. This local facultative sexuality may have a great impact on regional genotypic diversity. Two microsatellite loci discriminated genotypes separated by the accumulation of few mutations ('clone mates') within each AFLP clone. Genetic diversity is partitioned mainly among apomictic populations and is not geographically structured, which may be due to facultative sexuality and/or multiple colonizations of sites by different clones. Habitat differentiation and a tendency to inhabit artificial meadows is more pronounced in apomictic than in sexual populations. We hypothesize that maintenance of genetic diversity and superior colonizing abilities of apomicts in temporally and spatially heterogeneous environments are important for their distributional success.     

A Study on a Nearly Neutral Mutation Model in Finite Populations

As a nearly neutral mutation model, the house-of-cards model is studied in finite populations using computer simulations. The distribution of the mutant effect is assumed to be normal. The behavior is mainly determined by the product of the population size, N, and the standard deviation, sigma, of the distribution of the mutant effect. If 4N sigma is large compared to one, a few advantageous mutants are quickly fixed in early generations. Then most mutation becomes deleterious and very slow increase of the average selection coefficient follows. It takes very long for the population to reach the equilibrium state. Substitutions of alleles occur very infrequently in the later stage. If 4N sigma is the order of one or less, the behavior is qualitatively similar to that of the strict neutral case. Gradual increase of the average selection coefficient occurs and in generations of several times the inverse of the mutation rate the population almost reaches the equilibrium state. Both advantageous and neutral (including slightly deleterious) mutations are fixed. Except in the early stage, an increase of the standard deviation of the distribution of the mutant effect decreases the average heterozygosity. The substitution rate is reduced as 4N sigma is increased. Three tests of neutrality, one using the relationship between the average and the variance of heterozygosity, another using the relationship between the average heterozygosity and the average number of substitutions and Watterson's homozygosity test are applied to the consequences of the present model. It is found that deviation from the neutral expectation becomes apparent only when 4N sigma is more than two. Also a simple approximation for the model is developed which works well when the mutation rate is very small.     

Population genetic inferences using immune gene SNPs mirror patterns inferred by microsatellites

Single nucleotide polymorphisms (SNPs) are replacing microsatellites for population genetic analyses, but it is not apparent how many SNPs are needed or how well SNPs correlate with microsatellites. We used data from the gopher tortoise, Gopherus polyphemus—a species with small populations, to compare SNPs and microsatellites to estimate population genetic parameters. Specifically, we compared one SNP data set (16 tortoises from four populations sequenced at 17 901 SNPs) to two microsatellite data sets, a full data set of 101 tortoises and a partial data set of 16 tortoises previously genotyped at 10 microsatellites. For the full microsatellite data set, observed heterozygosity, expected heterozygosity and F_ST were correlated between SNPs and microsatellites; however, allelic richness was not. The same was true for the partial microsatellite data set, except that allelic richness, but not observed heterozygosity, was correlated. The number of clusters estimated by structure differed for each data set (SNPs = 2; partial microsatellite = 3; full microsatellite = 4). Principle component analyses (PCA) showed four clusters for all data sets. More than 800 SNPs were needed to correlate with allelic richness, observed heterozygosity and expected heterozygosity, but only 100 were needed for F_ST. The number of SNPs typically obtained from next‐generation sequencing (NGS) far exceeds the number needed to correlate with microsatellite parameter estimates. Our study illustrates that diversity, F_ST and PCA results from microsatellites can mirror those obtained with SNPs. These results may be generally applicable to small populations, a defining feature of endangered and threatened species, because theory predicts that genetic drift will tend to outweigh selection in small populations.     

The Persistence of Facultative Parthenogenesis in Drosophila albomicans

Parthenogenesis has evolved independently in more than 10 Drosophila species. Most cases are tychoparthenogenesis, which is occasional or accidental parthenogenesis in normally bisexual species with a low hatching rate of eggs produced by virgin females; this form is presumed to be an early stage of parthenogenesis. To address how parthenogenesis and sexual reproduction coexist in Drosophila populations, we investigated several reproductive traits, including the fertility, parthenogenetic capability, diploidization mechanisms, and mating propensity of parthenogenetic D. albomicans. The fertility of mated parthenogenetic females was significantly higher than that of virgin females. The mated females could still produce parthenogenetic offspring but predominantly produced offspring by sexual reproduction. Both mated parthenogenetic females and their parthenogenetic-sexual descendants were capable of parthenogenesis. The alleles responsible for parthenogenesis can be propagated through both parthenogenesis and sexual reproduction. As diploidy is restored predominantly by gamete duplication, heterozygosity would be very low in parthenogenetic individuals. Hence, genetic variation in parthenogenetic genomes would result from sexual reproduction. The mating propensity of females after more than 20 years of isolation from males was decreased. If mutations reducing mating propensities could occur under male-limited conditions in natural populations, decreased mating propensity might accelerate tychoparthenogenesis through a positive feedback mechanism. This process provides an opportunity for the evolution of obligate parthenogenesis. Therefore, the persistence of facultative parthenogenesis may be an adaptive reproductive strategy in Drosophila when a few founders colonize a new niche or when small populations are distributed at the edge of a species'' range, consistent with models of geographical parthenogenesis.     

Size homoplasy and mutational processes of interrupted microsatellites in two bee species, Apis mellifera and Bombus terrestris (Apidae)

Similar microsatellite electromorphs (PCR products of the same size) can arise from independent mutational events. Such alleles are not identical by descent. This phenomenon, termed size homoplasy, was studied by sequencing electromorphs of two microsatellite loci in which the stretch of basic repeats is interrupted by different short (1-2 bp) DNA motifs. The number and position of these interruptions were established for electromorphs from closely and distantly related populations of honeybees and bumblebees. No sequence difference was found when electromorphs came from the same subspecies or from closely related subspecies, suggesting that they were probably identical by descent. In contrast, sequence differences were often detected in distantly related subspecies, showing that size homoplasy frequently occurs at this level of population differentiation. Size homoplasy is increased by limits to free length variation of alleles, a phenomenon that seems to act on interrupted microsatellites when comparing distantly related taxa, that is, honeybee subspecies from different evolutionary lineages. Electromorph sequences suggest that, within the scope of these limits, large mutation events have occurred frequently at both interrupted loci studied. In good agreement with the molecular data, computations based on the observed heterozygosity and number of electromorphs and simulation studies showed that neither locus fits the one-step stepwise mutant model (SMM). We speculate that interrupted microsatellites in general could be characterized by a higher variance in repeat number and consequently a lower homoplasy rate than pure ones. Hence, interrupted microsatellites should be most appropriate for investigating population differentiation and evolutionary relationship between relatively distant populations.      

Genetic variation in BoLA microsatellite loci in Portuguese cattle breeds

Major histocompatibility complex (MHC) typing based on microsatellites can be a valuable approach to understanding the selective processes occurring at linked or physically close MHC genes and can provide important information on variability and relationships of populations. Using microsatellites within or in close proximity with bovine lymphocyte antigen (BoLA) genes, we investigated the polymorphisms in the bovine MHC, known as the BoLA, in eight Portuguese cattle breeds. Additional data from non-BoLA microsatellite loci were also used to compare the variability between these regions. Diversity was higher in BoLA than in non-BoLA microsatellites, as could be observed by the number of alleles, allelic richness and observed heterozygosity. Brava de Lide, a breed selected for aggressiveness and nobility, presented the lowest values of observed heterozygosity and allelic richness in both markers. Results from neutrality tests showed few statistically significant differences between the observed Hardy–Weinberg homozygosity ( F ) and the expected homozygosity ( F _E), indicating the apparent neutrality of the BoLA microsatellites within the analysed breeds. Nevertheless, we detected a trend of lower values of observed homozygosity compared with the expected one. We also detected some differences in the levels of allelic variability among the four BoLA microsatellites. Our data showed a higher number of alleles at the BoLA-DRB3 locus than at the BoLA-DRBP1 locus. These differences could be related to their physical position in the chromosome and may reflect functional requirements for diversity.     

Development of microsatellite markers and analysis of intraspecific genetic variability in chickpea (Cicer arietinum L.)

Paucity of polymorphic molecular markers in chickpea (Cicer arietinum L.) has been a major limitation in the improvement of this important legume. Hence, in an attempt to develop sequence-tagged microsatellite sites (STMS) markers from chickpea, a microsatellite enriched library from the C. arietinum cv. Pusa362 nuclear genome was constructed for the identification of (CA/GT) _n and (CT/GA) _n microsatellite motifs. A total of 92 new microsatellites were identified, of which 74 functional STMS primer pairs were developed. These markers were validated using 9 chickpea and one C. reticulatum accession. Of the STMS markers developed, 25 polymorphic markers were used to analyze the intraspecific genetic diversity within 36 geographically diverse chickpea accessions. The 25 primer pairs amplified single loci producing a minimum of 2 and maximum of 11 alleles. A total of 159 alleles were detected with an average of 6.4 alleles per locus. The observed and expected heterozygosity values averaged 0.32 (0.08–0.91) and 0.74 (0.23–0.89) respectively. The UPGMA based dendrogram was able to distinguish all the accessions except two accessions from Afghanistan establishing that microsatellites could successfully detect intraspecific genetic diversity in chickpea. Further, cloning and sequencing of size variant alleles at two microsatellite loci revealed that the variable numbers of AG repeats in different alleles were the major source of polymorphism. Point mutations were found to occur both within and immediately upstream of the long tracts of perfect repeats, thereby bringing about a conversion of perfect motifs into imperfect or compound motifs. Such events possibly occurred in order to limit the expansion of microsatellites and also lead to the birth of new microsatellites. The microsatellite markers developed in this study will be useful for genetic diversity analysis, linkage map construction as well as for depicting intraspecific microsatellite evolution.     

The rate of adaptation in asexuals

I study the population genetics of adaptation in asexuals. I show that the rate of adaptive substitution in an asexual species or nonrecombining chromosome region is a bell-shaped function of the mutation rate: at some point, increasing the mutation rate decreases the rate of substitution. Curiously, the mutation rate that maximizes the rate of adaptation depends solely on the strength of selection against deleterious mutations. In particular, adaptation is fastest when the genomic rate of mutation, U, equals the harmonic mean of selection coefficients against deleterious mutations, where we assume that selection for favorable alleles is milder than that against deleterious ones. This simple result is independent of the shape of the distribution of effects among favorable and deleterious mutations, population size, and the action of clonal interference. In the course of this work, I derive an approximation to the probability of fixation of a favorable mutation in an asexual genome or nonrecombining chromosome region in which both favorable and deleterious mutations occur.     

Comparative analysis of the within-population genetic structure in wild cherry (Prunus avium L.) at the self-incompatibility locus and nuclear microsatellites

Gametophytic self-incompatibility (SI) systems in plants exhibit high polymorphism at the SI controlling S-locus because individuals with rare alleles have a higher probability to successfully pollinate other plants than individuals with more frequent alleles. This process, referred to as frequency-dependent selection, is expected to shape number, frequency distribution, and spatial distribution of self-incompatibility alleles in natural populations. We investigated the genetic diversity and the spatial genetic structure within a Prunus avium population at two contrasting gene loci: nuclear microsatellites and the S-locus. The S-locus revealed a higher diversity (15 alleles) than the eight microsatellites (4-12 alleles). Although the frequency distribution of S-alleles differed significantly from the expected equal distribution, the S-locus showed a higher evenness than the microsatellites (Shannon's evenness index for the S-locus: E = 0.91; for the microsatellites: E = 0.48-0.83). Also, highly significant deviations from neutrality were found for the S-locus whereas only minor deviations were found for two of eight microsatellites. A comparison of the frequency distribution of S-alleles in three age-cohorts revealed no significant differences, suggesting that different levels of selection acting on the S-locus or on S-linked sites might also affect the distribution and dynamics of S-alleles. Autocorrelation analysis revealed a weak but significant spatial genetic structure for the multilocus average of the microsatellites and for the S-locus, but could not ascertain differences in the extent of spatial genetic structure between these locus types. An indirect estimate of gene dispersal, which was obtained to explain this spatial genetic pattern, indicated high levels of gene dispersal within our population (sigma(g) = 106 m). This high gene dispersal, which may be partly due to the self-incompatibility system itself, aids the effective gene flow of the microsatellites, thereby decreasing the contrast between the neutral microsatellites and the S-locus.     

High genetic diversity and low structure in an endemic long-lived tree,Yucca capensis (Asparagaceae)

• Endemic species distributed in fragmented habitats are highly vulnerable to extinction because they may have low genetic diversity. However, some life-history traits can mitigate the effect of genetic drift on populations. We analysed the level and distribution of genetic variation and ancestral population size of Yucca capensis, a long-lived endemic plant of the Baja California Peninsula, Mexico. Its populations are scattered across a habitat that is suffering accelerated transformation.
• We used six nuclear microsatellites to genotype 224 individuals from 17 locations across the entire species' geographic range. We estimated polymorphisms, heterozygosity and genetic structure. We also evaluated the ancestral and recent effective size and time since the population started to change.
• We found high heterozygosity, high polymorphism and low differentiation among locations, suggesting a panmictic population across the range. We also detected a large ancestral effective population size, which suffered a strong reduction in the Mid-Holocene.
• Despite changes in environmental conditions caused by habitat modification, the high diversity and low differentiation in Y. capensis may result from its large ancestral effective size and life-history traits, such as plant longevity, clonal growth and mating system, which reduce the rate of loss of genetic variation. However, the dependence on a specialist pollinator that displays short flight range can reduce gene flow among the plant populations and could, shortly, lead them into an extinction vortex.

     

Genetic structure and distinctness of Apis mellifera L. populations from the Canary Islands

The genetic structure of Apis mellifera populations from the Canary Islands has been assessed by mitochondrial (restriction fragment length polymorphisms of the intergenic transfer RNAleu-COII region) and nuclear (microsatellites) studies. These populations show a low level of genetic variation in terms of average number of alleles and degree of heterozygosity. Significant differences in the distribution of alleles were found in both data sets, confirming the genetic differentiation among some of the islands but not within them. Two mitochondrial haplotypes characteristic of the Canary Islands are found at high frequencies, although populations are introgressed by imported honeybees of eastern European C lineage. This introgression is rather high on Tenerife and El Hierro and low on Gran Canaria and La Gomera, whereas on La Palma it has not been recorded. The finding of microsatellite alleles characteristic of the eastern European lineage corroborates the genetic introgression. Phylogenetic analyses indicate that the Canarian honeybees are differentiated from other lineages and provide genetic evidence of their African origin.     

Effects of the shape of distribution of mutant effect in nearly neutral mutation models

Models of the theory of nearly neutral mutation incorporate a continuous distribution of mutation effects in contrast to the theory of purely neutral mutation which allows no mutations with intermediate effects. Previous studies of one such model, namely the house-of-cards mutation model, assumed normal distribution of mutation effect. Here I study the house-of-cards mutation model in random-mating finite populations using the weak-mutation approximation, paying attention to the effects of the distribution of mutant effects. The average selection coefficient, substitution rate and average heterozygosity in the equilibrium and transient states were studied mainly by computer simulation. The main findings are: (i) Very rapid decrease of the substitution rate and very slow approach to equilibrium as selection becomes stronger are characteristics of assuming normal distribution of mutant effect. If the right tail of the mutation distribution decays more rapidly than that of the normal distribution, the decrease of substitution rate becomes slower and equilibrium is achieved more quickly. (ii) The dispersion index becomes smaller or larger than 1 depending on the time and the intensity of selection, (iii) LetN be the population size. When selection is strong the ratio of 4N times the substitution rate to the average heterozygosity, which is expected to be 1 under neutrality, is larger than 1 in earlier generations but becomes less than 1 in later generations. These findings show the importance of the distribution of mutant effect and time in determination of the behaviour of various statistics frequently used in the study of molecular evolution.     

Fluctuating asymmetry and genetic variability in the roe deer (Capreolus capreolus): a test of the developmental stability hypothesis in mammals using neutral molecular markers

Fluctuating asymmetry (FA), used as an indicator of developmental stability, has long been hypothesized to be negatively correlated with genetic variability as a consequence of more variable organisms being better suited to buffer developmental pathways against environmental stress. However, it is still a matter of debate if this is due to metabolic properties of enzymes encoded by certain key loci or rather to overall genomic heterozygosity. Previous analyses suggest that there might be a general difference between homeo- and poikilotherms in that only the latter tend to exhibit the negative correlation predicted by theory. In the present study, we addressed these questions by analysing roe deer (Capreolus capreolus) from five German populations with regard to FA in metric and non-metric skull and mandible traits as well as variability at eight microsatellite loci. Genetic variability was quantified by heterozygosity and mean d2 parameters, and although the latter did not show any relationship with FA, we found for the first time a statistically significant negative correlation of microsatellite heterozygosity and non-metric FA among populations. Because microsatellites are non-coding markers, this may be interpreted as evidence for the role of overall genomic heterozygosity in determining developmental stability. To test if the threshold character of non-metric traits is responsible for the metric vs non-metric difference we also carried out calculations where we treated our metric traits as threshold values. This, however, did not yield significant correlations between FA and genetic variability either.     

Estimating genome-wide heterozygosity: effects of demographic history and marker type

Heterozygosity–fitness correlations (HFCs) are often used to link individual genetic variation to differences in fitness. However, most studies examining HFCs find weak or no correlations. Here, we derive broad theoretical predictions about how many loci are needed to adequately measure genomic heterozygosity assuming different levels of identity disequilibrium (ID), a proxy for inbreeding. We then evaluate the expected ability to detect HFCs using an empirical data set of 200 microsatellites and 412 single nucleotide polymorphisms (SNPs) genotyped in two populations of bighorn sheep (Ovis canadensis), with different demographic histories. In both populations, heterozygosity was significantly correlated across marker types, although the strength of the correlation was weaker in a native population compared with one founded via translocation and later supplemented with additional individuals. Despite being bi-allelic, SNPs had similar correlations to genome-wide heterozygosity as microsatellites in both populations. For both marker types, this association became stronger and less variable as more markers were considered. Both populations had significant levels of ID; however, estimates were an order of magnitude lower in the native population. As with heterozygosity, SNPs performed similarly to microsatellites, and precision and accuracy of the estimates of ID increased as more loci were considered. Although dependent on the demographic history of the population considered, these results illustrate that genome-wide heterozygosity, and therefore HFCs, are best measured by a large number of markers, a feat now more realistically accomplished with SNPs than microsatellites.     

Soft selection reduces loss of heterozygosity in asexual reproduction

The adaptive value of sexual reproduction is still debated in evolutionary theory. It has been proposed that the advantage of sexual reproduction over asexual reproduction is to promote genetic diversity, to prevent the accumulation of harmful mutations or to preserve heterozygosity. Since these hypothetical advantages depend on the type of asexual reproduction, understanding how selection affects the taxonomic distribution of each type could help us discriminate between existing hypotheses. Here, I argue that soft selection, competition among embryos or offspring in selection arenas prior to the hard selection of the adult phase, reduces loss of heterozygosity in certain types of asexual reproduction. Since loss of heterozygosity leads to the unmasking of recessive deleterious mutations in the progeny of asexual individuals, soft selection facilitates the evolution of these types of asexual reproduction. Using a population genetics model, I calculate how loss of heterozygosity affects fitness for different types of apomixis and automixis, and I show that soft selection significantly reduces loss of heterozygosity, hence increases fitness, in apomixis with suppression of the first meiotic division and in automixis with central fusion, the most common types of asexual reproduction. Therefore, if sexual reproduction evolved to preserve heterozygosity, soft selection should be associated with these types of asexual reproduction. I discuss the evidence for this prediction and how this and other observations on the distribution of different types of asexual reproduction in nature is consistent with the heterozygosity hypothesis.     

我国6个地方绵羊品种微卫星DNA多态性研究

利用聚丙烯酰胺凝胶电泳技术研究了我国蒙古羊、乌珠穆沁羊、哈萨克羊、阿勒泰羊、滩羊和藏绵羊 6个地方绵羊品种 17个微卫星标记的多态性 ,以探讨其遗传多样性、起源分化及群体间的遗传亲缘关系。结果表明微卫星标记不同位点间遗传多样性差异极显著 (P <0 0 1) ,群体间多态信息含量 (PIC)、近交程度 (Fis)和观察杂合度 (Obs .Het)差异不显著 ,但基因多样性 (genediversity)和期望杂合度 (Exp .Het)差异显著 (P <0 0 5 )。所研究的我国 6个地方绵羊品种与欧洲品种具有相似的遗传多样性 ,但具有较高的近交系数。个体和群体的聚类分析结果提示我国地方绵羊品种可能起源于两类祖先。群体间的聚类分析结果还表明 ,蒙古羊与乌珠穆沁羊分化不明显且具有较近的遗传亲缘关系 ,蒙古羊与藏绵羊间分化明显且具有较远的遗传亲缘关系。滩羊、阿勒泰羊以及藏绵羊间也具有较近的遗传亲缘关系。所研究的我国 6个地方绵羊品种的遗传分化 (Fst)与西班牙绵羊品种接近 ,但明显小于欧洲其他绵羊品种     

A comparison of biallelic markers and microsatellites for the estimation of population and conservation genetic parameters in Atlantic salmon (Salmo salar)

Biallelic markers such as single nucleotide polymorphisms (SNPs) and insertion/deletion polymorphisms have become increasingly popular markers for various population genetics applications. However, the effort required to develop biallelic markers in nonmodel organisms is still substantial. In this study, we compared the estimation of various population genetic parameters (genetic divergence and structuring, isolation-by-distance, genetic diversity) using a limited number of biallelic markers (in total 7 loci) to those estimated with 14 microsatellite loci in 21 Atlantic salmon (Salmo salar) populations from northern Europe. Pairwise FST values were significantly correlated between biallelic loci and microsatellite datasets, as was overall heterozygosity when both anadromous and nonanadromous populations were analyzed together. However, when the anadromous and nonanadromous samples were analyzed separately, only genetic divergence correlations remained significant. Biallelic markers alone were not sufficient for reliable neighbor-joining clustering of populations but gave highly similar isolation-by-distance signals when compared with microsatellites. Finally, although several population prioritization measures for conservation exhibited significant correlation between different marker types, the specific populations highlighted as being most valuable for conservation purposes varied depending on the marker type and conservation criteria applied. This study demonstrates that a relatively small set of biallelic markers can be sufficient for obtaining concordant results in most of the analyses compared with microsatellites, although estimates of genetic distance are generally more concordant than estimates of genetic diversity. This suggests that a relatively small number of biallelic markers can provide useful information for various population genetic applications. However, we emphasize that the use of much higher number of loci is preferable, especially when the genetic differences between populations are subtle or individual multilocus genotype-based analyses are to be performed.