Population genetic structure of the sea bass (Lateolabrax japonicus) in Korea based on multiplex PCR assays with 12 polymorphic microsatellite markers

Population genetics has been recognized as a key component of policy development for fisheries and conservation management. In this study, natural sea bass (Lateolabrax japonicus) populations in three ocean basins in Korea were assessed using multiplex assays with 12 highly polymorphic microsatellite loci; 203 alleles and similarly high levels of genetic diversity [mean number of alleles (N_A) = 14.43, mean expected heterozygosity (He) = 0.84] were detected. All populations showed significant heterozygote deficiency at four loci, which could be explained by the presence of null alleles. The genetic population subdivision was low and was significantly different according to F-statistics (overall F _ST = 0.003, R _ST = 0.005). However, this substructure was not supported by an analysis of molecular variance test, analyses of isolation by distance or Bayesian analysis. The passive dispersal of eggs/larvae via the main currents appears to facilitate gene flow. The possibility of a recent genetic bottleneck was observed in all three populations of L. japonicus, indicating that overfishing and degradation of the environment in recent years has led to a decline in the sea bass populations in Korea. Our study demonstrates that sea bass in Korea do not appear to be genetically partitioned and should be managed as a single unit; however, the potential for a rapid loss of genetic diversity remains. Information regarding the genetic characteristics of Korean sea bass populations has important implications for fishery management and conservation efforts and will aid in the sustainable exploitation of fishing resources and the preservation of biodiversity.     

Combining multiple analytical approaches for the identification of population structure and genetic delineation of two subspecies of the endemic Arabian burnet moth Reissita simonyi (Zygaenidae; Lepidoptera)

Habitat fragmentation and landscape topology may influence the genetic structure and connectivity between natural populations. Six microsatellite loci were used to infer the population structure of 35 populations (N = 788) of the alpine Arabian burnet moth Reissita simonyi (Lepidoptera, Zygaenidae) in Yemen and Oman. Due to the patchy distribution of larval food plants, R. simonyi is not continuously distributed throughout the studied area and the two recognized subspecies of this endemic species (Reissita s. simonyi/R. s. yemenicola) are apparently discretely distributed. All microsatellites showed prevalence of null alleles and therefore a thorough investigation of the impact of null alleles on different population genetic parameters (F _ST, inbreeding coefficients, and Population Graph topologies) is given. In general, null alleles reduced genetic covariance and independence of allele frequencies resulting in a more connected genetic topology in Population Graphs and an overestimation of pairwise F _ST values and inbreeding coefficients. Despite the presence of null alleles, Population Graphs also showed a much higher genetic connectivity within subspecies (and lower genetic differentiation (via F _ST)) than between; supporting existing taxonomic distinction. Partial Mantel tests showed that both geographical distance and altitude were highly correlated with the observed distribution of genetic structure within R. simonyi. In conclusion, we identified geographical and altitudinal distances in R. simonyi as well as an intervening desert area to be the main factors for spatial genetic structure in this species and show that the taxonomic division into two subspecies is confirmed by genetic analysis.     

Estimation and adjustment of microsatellite null alleles in nonequilibrium populations

Nonamplified (null) alleles are a common feature of microsatellite genotyping and can bias estimates of allele and genotype frequencies, thereby hindering population genetic analyses. The frequency of microsatellite null alleles in diploid populations can be estimated for populations that are in Hardy–Weinberg equilibrium. However, many microsatellite data sets are from nonequilibrium populations, often with known inbreeding coefficients (F) or fixation indices (F_IS or F_ST). Here, we propose a novel null allele estimator that can be used to estimate the null allele frequency and adjust visible allele frequencies in populations for which independent estimates of F, F_IS or F_ST are available. The algorithm is currently available as an Excel macro that can be downloaded at no cost from http://www.microchecker.hull.ac.uk/ and will be incorporated into the software micro ‐checker .     

Population structure of the malaria vector Anopheles sinensis (Diptera: Culicidae) in China: two gene pools inferred by microsatellites

Background

Anopheles sinensis is a competent malaria vector in China. An understanding of vector population structure is important to the vector-based malaria control programs. However, there is no adequate data of A. sinensis population genetics available yet.

Methodology/Principal Findings

This study used 5 microsatellite loci to estimate population genetic diversity, genetic differentiation and demographic history of A. sinensis from 14 representative localities in China. All 5 microsatellite loci were highly polymorphic across populations, with high allelic richness and heterozygosity. Hardy–Weinberg disequilibrium was found in 12 populations associated with heterozygote deficits, which was likely caused by the presence of null allele and the Wahlund effect. Bayesian clustering analysis revealed two gene pools, grouping samples into two population clusters; one includes six and the other includes eight populations. Out of 14 samples, six samples were mixed with individuals from both gene pools, indicating the coexistence of two genetic units in the areas sampled. The overall differentiation between two genetic pools was moderate (F _ST = 0.156). Pairwise differentiation between populations were lower within clusters (F _ST = 0.008–0.028 in cluster I and F _ST = 0.004–0.048 in cluster II) than between clusters (F _ST = 0.120–0.201). A reduced gene flow (Nm = 1–1.7) was detected between clusters. No evidence of isolation by distance was detected among populations neither within nor between the two clusters. There are differences in effective population size (Ne = 14.3-infinite) across sampled populations.

Conclusions/Significance

Two genetic pools with moderate genetic differentiation were identified in the A. sinensis populations in China. The population divergence was not correlated with geographic distance or barrier in the range. Variable effective population size and other demographic effects of historical population perturbations could be the factors affecting the population differentiation. The structured populations may limit the migration of genes under pressures/selections, such as insecticides and immune genes against malaria.     

Signs of selective pressure on genetic variants affecting human height

Many decades of scientific investigation have proved the role of selective pressure in Homo Sapiens at least at the level of individual genes or loci. Nevertheless, there are examples of polygenic traits that are bound to be under selection, but studies devoted to apply population genetics methods to unveil such occurrence are still lacking. Stature provides a relevant example of well-studied polygenic trait for which is now available a genome-wide association study which has identified the genes involved in this trait, and which is known to be under selection. We studied the behavior of F_ST in a simulated toy model to detect population differentiation on a generic polygenic phenotype under selection. The simulations showed that the set of alleles involved in the trait has a higher mean F_ST value than those undergoing genetic drift only. In view of this we looked for an increase in the mean F_ST value of the 180 variants associated to human height. For this set of alleles we found F_ST to be significantly higher than the genomic background (p = 0.0356). On the basis of a statistical analysis we excluded that the increase was just due to the presence of outliers. We also proved as marginal the role played by local adaptation phenomena, even on different phenotypes in linkage disequilibrium with genetic variants involved in height. The increase of F_ST for the set of alleles involved in a polygenic trait seems to provide an example of symmetry breaking phenomenon concerning the population differentiation. The splitting in the allele frequencies would be driven by the initial conditions in the population dynamics which are stochastically modified by events like drift, bottlenecks, etc, and other stochastic events like the born of new mutations.     

Differential introgression from a sister species explains high FST outlier loci within a mussel species

Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. F_ST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at F_ST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (F_ST = 0.03), but seven loci (2%) were strong F_ST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.     

Population genetic structure of three major river Populations of rohu,Labeo rohita (Cyprinidae: Cypriniformes) using microsatellite DNA markers

Determining the genetic structure is essential for developing conservation and stock improvement plans. Four dinucleotide microsatellite loci were analysed to reveal population genetic structure of the Indian major carp,Labeo rohita collected from three major rivers namely the Halda, the Jamuna, and the Padma in Bangladesh. The four loci were polymorphic (P ₉₅) in all the populations. The populations varied in the number and frequencies of alleles as well as heterozygosities in the loci analyzed. Population differentiation (F _ST) value between the Halda and the Jamuna population was significant (P<0.05). Relatively high level of gene flow and low level ofF _ST values were found between the Padma and the Jamuna population. The unweighted pair group method with averages (UPGMA) dendrogram based on genetic distance resulted in two clusters: the Halda population was alone in one cluster whereas the Jamuna and the Padma made another cluster. The results revealed a relatively low level of genetic variability in the river populations ofL. rohita in Bangladesh.     

Population estimators or progeny tests: what is the best method to assess null allele frequencies at SSR loci?

Nuclear SSRs are notorious for having relatively high frequencies of null alleles, i.e. alleles that fail to amplify and are thus recessive and undetected in heterozygotes. In this paper, we compare two kinds of approaches for estimating null allele frequencies at seven nuclear microsatellite markers in three French Fagus sylvatica populations: (1) maximum likelihood methods that compare observed and expected homozygote frequencies in the population under the assumption of Hardy-Weinberg equilibrium and (2) direct null allele frequency estimates from progeny where parent genotypes are known. We show that null allele frequencies are high in F. sylvatica (7.0% on average with the population method, 5.1% with the progeny method), and that estimates are consistent between the two approaches, especially when the number of sampled maternal half-sib progeny arrays is large. With null allele frequencies ranging between 5% and 8% on average across loci, population genetic parameters such as genetic differentiation (F _ST) may be mostly unbiased. However, using markers with such average prevalence of null alleles (up to 15% for some loci) can be seriously misleading in fine scale population studies and parentage analysis.     

Population genetics of the common long-armed octopus Octopus minor (Sasaki, 1920) (Cephalopoda: Octopoda) in Chinese waters based on microsatellite analysis

Octopus minor (Sasaki, 1920) is a commercially important cephalopod in Chinese waters. To provide a theoretical basis for resource protection and sustainable management, we investigated genetic structure of ten O. minor populations in Chinese waters using microsatellite DNA markers. Eight microsatellite loci revealed high allelic diversity with 11–26 alleles per locus. Observed and expected heterozygosity varied from 0.412 to 0.900 and from 0.337 to 0.845, respectively. The overall F_ST value was 0.198, indicating great genetic differentiation among populations. The F_ST value between Yilan and other populations reached more than 0.3 that may be indicative of subspecies rank. Mantel test showed significant correlations between genetic and geographic distance (R = 0.383, P = 0.004) indicating that genetic differentiation of O. minor conformed to a pattern of isolation-by-distance. Using the Neighbor-joining method, cluster analysis divided nine populations into three groups and divided ten populations into two groups wherein Yilan was distinguished from the other populations. Analysis based on F_ST, D_c values and clustering highlighted the heterogeneity of Yilan and the relative homogeneity between Yilan and Ganyu. The significant population genetic structure of O. minor is related to the combined effects of geographical barriers, current features and life history characteristics.     

Microsatellite loci in the house fly,Musca domestica L. (Diptera: Muscidae)

Genomic libraries from house flies enriched for (CA)₁₅ and (CAG)₁₀ repeats were constructed by using biotinylated probes. Twenty‐five loci were isolated and evaluated for polymorphisms in wild flies representing two geographically diverse populations. Fourteen of 19 dinucleotide loci, and one of six trinucleotide loci were polymorphic. One hundred and twenty‐seven alleles were detected, 39 of which were private. Average number of alleles per polymorphic locus was 8.4 ± 2.5 and average heterozygosity was 72 ± 4%. F_ST by the private allele method was 0.73. Three of 15 loci showed significant heterozygote deficiencies, attributed to null alleles. Five of 15 loci were amplified in the face fly, Musca autumnalis.     

Genetic variation and differentiation in Juniperus excelsa M. Bieb. populations in Turkey

Genetic variation within and among six populations of Juniperus excelsa M. Bieb., in a common garden in Lakes District of Turkey, was analyzed using four nuclear microsatellite primer pairs originally developed for J. communis. A total of five loci were observed as Jc037 amplified two distinct size ranges. The number of alleles observed for the species varied from 2 to 13, with an average of 4.5 alleles per locus. The mean expected heterozygosity (H _e) of populations was 0.584, after correction for null alleles. The mean F _IS value (?0.014) was close to zero showing no significant deviation from Hardy–Weinberg equilibrium. A low level of genetic differentiation was observed among populations (F _ST = 0.028; p < 0.001) and Nei’s genetic distance ranged from 0.014 to 0.120 between population pairs. Furthermore, there was no significant correlation between genetic distances on the one hand and geographic distances and trait differentiation on the other hand. However, the eastern populations Bey?ehir and Sorgun showed very similar genotypic structures and were differentiated from all other populations. A continuous monitoring of phenotypic traits and the association between nucleotide variation in functional genes and adaptive traits such as drought and frost tolerance of J. excelsa populations in common gardens will be useful to design effective conservation strategies in the future.     

Genetic diversity,population structure and phenotypic variation in European Salix viminalis L. (Salicaceae)

To investigate the potential of association genetics for willow breeding, Salix viminalis germplasm was assembled from UK and Swedish collections (comprising accessions from several European countries) and new samples collected from nature. A subset of the germplasm was planted at two sites (UK and Sweden), genotyped using 38 SSR markers and assessed for phenological and biomass traits. Population structure, genetic differentiation (F _ST ) and quantitative trait differentiation (Q _ST ) were investigated. The extent and patterns of trait adaptation were assessed by comparing F _ST ?and Q _ST parameters. Of the 505 genotyped diploid accessions, 27 % were not unique. Genetic diversity was high: 471 alleles was amplified; the mean number of alleles per locus was 13.46, mean observed heterozygosity was 0.55 and mean expected heterozygosity was 0.62. Bayesian clustering identified four subpopulations which generally corresponded to Western Russia, Western Europe, Eastern Europe and Sweden. All pairwise F _ST values were highly significant (p?F _ST ???=???0.12), and the smallest between the Swedish and Eastern European populations (F _ST ??=??0.04). The Swedish population also had the highest number of identical accessions, supporting the view that S. viminalis was introduced into this country and has been heavily influenced by humans. Q _ST values were high for growth cessation and leaf senescence, and to some extent stem diameter, but low for bud burst time and shoot number. Overall negative clines between longitudinal coordinates and leaf senescence, bud burst and stem diameter were also found.     

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.     

A spectral theory for Wright’s inbreeding coefficients and related quantities

Wright’s inbreeding coefficient, F_ST, is a fundamental measure in population genetics. Assuming a predefined population subdivision, this statistic is classically used to evaluate population structure at a given genomic locus. With large numbers of loci, unsupervised approaches such as principal component analysis (PCA) have, however, become prominent in recent analyses of population structure. In this study, we describe the relationships between Wright’s inbreeding coefficients and PCA for a model of K discrete populations. Our theory provides an equivalent definition of F_ST based on the decomposition of the genotype matrix into between and within-population matrices. The average value of Wright’s F_ST over all loci included in the genotype matrix can be obtained from the PCA of the between-population matrix. Assuming that a separation condition is fulfilled and for reasonably large data sets, this value of F_ST approximates the proportion of genetic variation explained by the first (K − 1) principal components accurately. The new definition of F_ST is useful for computing inbreeding coefficients from surrogate genotypes, for example, obtained after correction of experimental artifacts or after removing adaptive genetic variation associated with environmental variables. The relationships between inbreeding coefficients and the spectrum of the genotype matrix not only allow interpretations of PCA results in terms of population genetic concepts but extend those concepts to population genetic analyses accounting for temporal, geographical and environmental contexts.     

Genetic diversity among alfalfa (Medicago sativa) cultivars coming from a breeding program, using SSR markers

Alfalfa (Medicago sativa) is an autotetraploid, allogamous and heterozygous species whose cultivars are synthetic populations. The breeders apply selection pressure for some agronomic traits within a breeding pool to increase the frequency of favorable individuals. The objective of this study was to investigate the differentiation level among seven cultivars originating from one breeding program, and between these cultivars and the breeding pool, with eight SSR markers. These highly polymorphic and codominant markers, together with recent population genetic statistics extended to autotetraploids, offer tools to analyse genetic diversity in alfalfa. The number of alleles per locus varied between 3 and 24. All loci were at a panmictic equilibrium in the cultivars, except one, probably because of null alleles. With seven SSR loci, each cultivar was at panmictic equilibrium. The mean gene diversity was high, ranging from 0.665 to 0.717 in the cultivars. The parameter F _ST indicated a low but significant diversity among cultivars. Among 21 pairs of cultivars, 15 were significantly different. The breeding pool also had a high diversity, and was significantly different from each cultivar except the most recent one. Considering the characteristics of the breeding program and the mode of cultivar elaboration, we found that they were unable to generate a large variety differentiation. Estimation of population genetics parameters at SSR loci can be applied for assessing the differences between cultivars or populations, either for variety distinction or the management of genetic resources.     

F(ST) and Q(ST) under neutrality

A commonly used test for natural selection has been to compare population differentiation for neutral molecular loci estimated by F_ST and for the additive genetic component of quantitative traits estimated by Q_ST. Past analytical and empirical studies have led to the conclusion that when averaged over replicate evolutionary histories, Q_ST = F_ST under neutrality. We used analytical and simulation techniques to study the impact of stochastic fluctuation among replicate outcomes of an evolutionary process, or the evolutionary variance, of Q_ST and F_ST for a neutral quantitative trait determined by n unlinked diallelic loci with additive gene action. We studied analytical models of two scenarios. In one, a pair of demes has recently been formed through subdivision of a panmictic population; in the other, a pair of demes has been evolving in allopatry for a long time. A rigorous analysis of these two models showed that in general, it is not necessarily true that mean Q_ST = F_ST (across evolutionary replicates) for a neutral, additive quantitative trait. In addition, we used finite-island model simulations to show there is a strong positive correlation between Q_ST and the difference Q_ST − F_ST because the evolutionary variance of Q_ST is much larger than that of F_ST. If traits with relatively large Q_ST values are preferentially sampled for study, the difference between Q_ST and F_ST will also be large and positive because of this correlation. Many recent studies have used tests of the null hypothesis Q_ST = F_ST to identify diversifying or uniform selection among subpopulations for quantitative traits. Our findings suggest that the distributions of Q_ST and F_ST under the null hypothesis of neutrality will depend on species-specific biology such as the number of subpopulations and the history of subpopulation divergence. In addition, the manner in which researchers select quantitative traits for study may introduce bias into the tests. As a result, researchers must be cautious before concluding that selection is occurring when Q_ST ≠ F_ST.     

Heterozygote deficiencies in parasite populations: an evaluation of interrelated hypotheses in the raccoon tick,Ixodes texanus

Population genetics is increasingly being used to study the biology of parasites at the scales of both the host (infrapopulation, IP) and host population (component population, CP). In this study we tested three mechanistic hypotheses that could explain deviations from Hardy–Weinberg equilibrium (HWE) expectations due to heterozygote deficits (HDs) at the CP scale in raccoon ticks (Ixodes texanus; n=718) collected from raccoons (Procyon lotor; n=91) and genotyped at 11 microsatellite loci. These hypotheses were presence of technical issues (for example, null alleles), hierarchical structure (for example, host demography) and cryptic structure (for example, kin structure). Although statistical support for null alleles existed, their presence would also be expected to lead to an underestimation in levels of relatedness, and thus kin structure. However, we found the opposite pattern: significant HD at the IP scale being more likely in CPs with significant vs non-significant levels of kin structure. Our analyses revealed that pooling of kin groups could lead to highly variable levels of F_IS among loci, a pattern usually suggestive of null alleles. We used Monte–Carlo (MC) simulations to show that the existence of subdivided breeding groups and high variance in individual reproductive success could adequately explain deviations from HWE in I. texanus. Thus, our results indicate that biological factors can lead to patterns that have usually been interpreted as technical issues (for example, null alleles), and that it is important to take such factors into consideration because loci deviating from HWE likely reflect the effects of real biological processes.     

Temporal genetic variation in a population of the pocket gopher,Geomys bursarius

Genetic variation at eight polymorphic loci was examined in a population of a subterranean rodent, the plains pocket gopher (Geomys bursarius), sampled over a 10-yr period. Two loci exhibited relatively minor changes in gene frequencies, while the remaining loci demonstrated major shifts in predominant alleles, loss of minor alleles, and addition of alleles due to migration. Significant deviations from Hardy-Weinberg expectations, concomitant to heterozygote deficiencies, were observed for several loci. Temporal heterogeneity, as measured by F_ST, was high and comparable to that exhibited by local populations sampled over relatively short periods of time. The high degree of temporal genetic variation is consistent with observations that fossorial rodents occur in locally isolated populations with small effective population sizes that are subject to genetic drift, bottlenecking, and inbreeding.     

The joint effects of selection and dominance on the QST - FST contrast

Q_ST measures the differentiation of quantitative traits between populations. It is often compared to F_ST, which measures population differentiation at neutral marker loci due to drift, migration, and mutation. When Q_ST is different from F_ST, it is usually taken as evidence that selection has either restrained or accelerated the differentiation of the quantitative trait relative to neutral markers. However, a number of other factors such as inbreeding, dominance, and epistasis may also affect the Q_ST − F_ST contrast. In this study, we examine the effects of dominance, selection, and inbreeding on Q_ST − F_ST. We compare Q_ST with F_ST at selected and neutral loci for populations at equilibrium between selection, drift, mutation, and migration using both analytic and simulation approaches. Interestingly, when divergent selection is acting on a locus, inbreeding and dominance generally inflate Q_ST relative to F_ST when they are both measured at the quantitative locus at equilibrium. As a consequence, dominance is unlikely to hide the signature of divergent selection on the Q_ST − F_ST contrast. However, although in theory dominance and inbreeding affect the expectation for Q_ST − F_ST, of most concern is the very large variance in both Q_ST and F_ST, suggesting that we should be cautious in attributing small differences between Q_ST and F_ST to selection.