Confidence Intervals for Population Allele Frequencies: The General Case of Sampling from a Finite Diploid Population of Any Size

The estimation of population allele frequencies using sample data forms a central component of studies in population genetics. These estimates can be used to test hypotheses on the evolutionary processes governing changes in genetic variation among populations. However, existing studies frequently do not account for sampling uncertainty in these estimates, thus compromising their utility. Incorporation of this uncertainty has been hindered by the lack of a method for constructing confidence intervals containing the population allele frequencies, for the general case of sampling from a finite diploid population of any size. In this study, we address this important knowledge gap by presenting a rigorous mathematical method to construct such confidence intervals. For a range of scenarios, the method is used to demonstrate that for a particular allele, in order to obtain accurate estimates within 0.05 of the population allele frequency with high probability (%), a sample size of is often required. This analysis is augmented by an application of the method to empirical sample allele frequency data for two populations of the checkerspot butterfly (Melitaea cinxia L.), occupying meadows in Finland. For each population, the method is used to derive % confidence intervals for the population frequencies of three alleles. These intervals are then used to construct two joint % confidence regions, one for the set of three frequencies for each population. These regions are then used to derive a % confidence interval for Jost''s D, a measure of genetic differentiation between the two populations. Overall, the results demonstrate the practical utility of the method with respect to informing sampling design and accounting for sampling uncertainty in studies of population genetics, important for scientific hypothesis-testing and also for risk-based natural resource management.     

Killer cell immunoglobulin-like receptor (KIR) gene content variation in the HGDP-CEPH populations

In the present study, we investigate patterns of variation in the KIR cluster in a large and well-characterized sample of worldwide human populations in the Human Genome Diversity Project-Centre d'Etude du Polymorphisme Humain (HGDP-CEPH) panel in order to better understand the patterns of diversity in the region. Comparison of KIR data with that from other genomic regions allows control for strictly demographic factors; over 500,000 additional genomic markers have been typed in this panel by other investigators and the data made publicly available. Presence/absence frequencies and haplotypic associations for the KIR region are analyzed in the 52 populations comprising the panel and in accordance with major world regions (Africa, Middle East, Central Asia, East Asia, Europe, Americas, and Oceania). These data represent the first overview of KIR population genetics in the well-documented HGDP-CEPH panel and suggest different evolutionary histories and recent selection in the KIR gene cluster.     

Assessing the Effect of Sequencing Depth and Sample Size in Population Genetics Inferences

Next-Generation Sequencing (NGS) technologies have dramatically revolutionised research in many fields of genetics. The ability to sequence many individuals from one or multiple populations at a genomic scale has greatly enhanced population genetics studies and made it a data-driven discipline. Recently, researchers have proposed statistical modelling to address genotyping uncertainty associated with NGS data. However, an ongoing debate is whether it is more beneficial to increase the number of sequenced individuals or the per-sample sequencing depth for estimating genetic variation. Through extensive simulations, I assessed the accuracy of estimating nucleotide diversity, detecting polymorphic sites, and predicting population structure under different experimental scenarios. Results show that the greatest accuracy for estimating population genetics parameters is achieved by employing a large sample size, despite single individuals being sequenced at low depth. Under some circumstances, the minimum sequencing depth for obtaining accurate estimates of allele frequencies and to identify polymorphic sites is , where both alleles are more likely to have been sequenced. On the other hand, inferences of population structure are more accurate at very large sample sizes, even with extremely low sequencing depth. This all points to the conclusion that under various experimental scenarios, in cost-limited population genetics studies, large sample sizes at low sequencing depth are desirable to achieve high accuracy. These findings will help researchers design their experimental set-ups and guide further investigation on the effect of protocol design for genetic research.     

The influence of artificial cranial deformation on discontinuous morphological traits

Following the elucidation by geneticists of the nature of minor skeletal variants in the mouse, anthropologists have stressed the potential of these traits for tracing the affinities and movements of extinct human populations. Earlier Sullivan observed that discrete traits could be particularly valuable where artificial cranial deformation limits the use of craniometry. Twenty-eight minor variants were studied in bifronto-occipitally deformed and undeformed skulls of a sample of 78 from a single Hopewell mound. The pattern of frequency differences between deformed and undeformed with respect to traits at the back of the vault and in the frontal region, interpreted in developmental terms, reveals a hypostotic effect in these regions in the deformed skull; while, in contrast, traits of the lateral vault, facial skeleton and cranial base point to a general hyperostotic effect in these regions. Each of three emissaria tends to be more constant in the deformed. That minor cranial variants manifest a plastic response to this type of environmentally-imposed stress is consistent with the nature of such variants as elucidated by genetics research in mice. The findings suggest that deformed crania should be excluded from population studies in which genetic divergence between groups is estimated in terms of cranial trait frequencies.     

Hitchhiking and recombination in birds: evidence from Mhc-linked and unlinked loci in Red-winged Blackbirds (Agelaius phoeniceus)

Hitchhiking phenomena and genetic recombination have important consequences for a variety of fields for which birds are model species, yet we know virtually nothing about naturally occurring rates of recombination or the extent of linkage disequilibrium in birds. We took advantage of a previously sequenced cosmid clone from Red-winged Blackbirds (Agelaius phoeniceus) bearing a highly polymorphic Mhc class II gene, Agph-DABI, to measure the extent of linkage disequilibrium across approximately 40 kb of genomic DNA and to determine whether non-coding nucleotide diversity was elevated as a result of physical proximity to a target of balancing selection. Application of coalescent theory predicts that the hitchhiking effect is enhanced by the larger effective population size of blackbirds compared with humans, despite the presumably higher rates of recombination in birds. We surveyed sequence polymorphism at three Mhc-linked loci occurring 1.5-40 kb away from Agph-DAB1 and found that nucleotide diversity was indistinguishable from that found at three presumably unlinked, non-coding introns (beta-actin intron 2, beta-fibrinogen intron 7 and rhodopsin intron 2). Linkage disequilibrium as measured by Lewontin's D' was found only across a few hundred base pairs within any given locus, and was not detectable among any Mhc-linked loci. Estimated rates of the per site recombination rate p derived from three different analytical methods suggest that the amounts of recombination in blackbirds are up to two orders of magnitude higher than in humans, a discrepancy that cannot be explained entirely by the higher effective population size of blackbirds relative to humans. In addition, the ratio of the number of estimated recombination events per mutation frequently exceeds 1, as in Drosophila, again much higher than estimates in humans. Although the confidence limits of the blackbird estimates themselves span an order of magnitude, these data suggest that in blackbirds the hitchhiking effect for this region is negligible and may imply that the per site per individual recombination rate is high, resembling those of Drosophila more than those of humans.     

A population genetic analysis of migration: the case of the noctule bat (Nyctalus noctula)

Although rarely assessed, the population genetics of hibernating colonies can help to understand some aspects of population structure, even when samples from nursery or mating colonies are not available, or in studies of migration when both types of samples are available and can be compared. Here we illustrate both points in a survey of mitochondrial DNA (mtDNA) control region sequences used to study the population genetics of hibernating colonies of a migrating species, the noctule bat (Nyctalus noctula). Lacking samples from Scandinavian nursery colonies, we use a North European hibernacula to suggest that Scandinavian populations are isolated from Central and East European colonies. Then, we compare genetic diversities of nursery and hibernating colonies. We find a significantly higher haplotype diversity in hibernacula, confirming that they consist of individuals from different nursery colonies. Finally, we show that pairwise comparisons of the haplotype frequencies of nursery and hibernating colonies contain some information on the migration direction of the noctule bat.     

Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants

A compilation was made of 307 studies using nuclear DNA markers for evaluating among- and within-population diversity in wild angiosperms and gymnosperms. Estimates derived by the dominantly inherited markers (RAPD, AFLP, ISSR) are very similar and may be directly comparable. STMS analysis yields almost three times higher values for within-population diversity whereas among-population diversity estimates are similar to those derived by the dominantly inherited markers. Number of sampled plants per population and number of scored microsatellite DNA alleles are correlated with some of the population genetics parameters. In addition, maximum geographical distance between sampled populations has a strong positive effect on among-population diversity. As previously verified with allozyme data, RAPD- and STMS-based analyses show that long-lived, outcrossing, late successional taxa retain most of their genetic variability within populations. By contrast, annual, selfing and/or early successional taxa allocate most of the genetic variability among populations. Estimates for among- and within-population diversity, respectively, were negatively correlated. The only major discrepancy between allozymes and STMS on the one hand, and RAPD on the other hand, concerns geographical range; within-population diversity was strongly affected when the former methods were used but not so in the RAPD-based studies. Direct comparisons between the different methods, when applied to the same plant material, indicate large similarities between the dominant markers and somewhat lower similarity with the STMS-based data, presumably due to insufficient number of analysed microsatellite DNA loci in many studies.     

Genome and population dynamics under selection and neutrality: an example of S-allele diversity in wild cherry (Prunus avium L.)

Self-incompatibility, a common attribute of plant development, forms a classical paradigm of balancing selection in natural populations, in particular negative frequency-dependent selection. Under negative frequency-dependent selection population genetics theory predicts that the S-locus, being in command of self-incompatibility, keeps numerous alleles in equal frequencies demonstrating a wide allelic range. Moreover, while natural populations exhibit a higher within population genetic diversity, a reduction of population differentiation and increase of effective migration rate is expected in comparison to neutral loci. Allelic frequencies were investigated in terms of distribution and genetic structure at the gametophytic self-incompatibility locus in five wild cherry (Prunus avium L.) populations. Comparisons were also made between the differentiation at the S-locus and at the SSR loci. Theoretical expectations under balancing selection were congruent to the results observed. The S-locus showed broad multiplicity (16 S-alleles), high genetic diversity, and allelic isoplethy. Genetic structure at the self-incompatibility locus was almost four times lower than at 11 nSSR loci. Analysis of molecular variance revealed that only 5?% of the total genetic variation concerns differentiation among populations. In conclusion, the wealth of S-allelic diversity found in natural wild cherry populations in Greece is useful not only in advancing basic population genetics research of self-incompatibility systems in wild cherry but also in the development of breeding programs.     

Population genetics and the evolution of geographic range limits in an annual plant

Abstract Theoretical models of species' geographic range limits have identified both demographic and evolutionary mechanisms that prevent range expansion. Stable range limits have been paradoxical for evolutionary biologists because they represent locations where populations chronically fail to respond to selection. Distinguishing among the proposed causes of species' range limits requires insight into both current and historical population dynamics. The tools of molecular population genetics provide a window into the stability of range limits, historical demography, and rates of gene flow. Here we evaluate alternative range limit models using a multilocus data set based on DNA sequences and microsatellites along with field demographic data from the annual plant Clarkia xantiana ssp. xantiana. Our data suggest that central and peripheral populations have very large historical and current effective population sizes and that there is little evidence for population size changes or bottlenecks associated with colonization in peripheral populations. Whereas range limit populations appear to have been stable, central populations exhibit a signature of population expansion and have contributed asymmetrically to the genetic diversity of peripheral populations via migration. Overall, our results discount strictly demographic models of range limits and more strongly support evolutionary genetic models of range limits, where adaptation is prevented by a lack of genetic variation or maladaptive gene flow.     

Post-colonization temporal genetic variation of an introduced fly, Rhagoletis completa

Evolutionary biologists have been puzzled by the success of introduced species: despite founder effects that reduce genetic variability, invasive species are still successful at colonizing new environments. It is possible that the evolutionary processes during the post-colonization period may increase the genetic diversity and gene flow among invasive populations over time, facilitating their long-term success. Therefore, genetic diversity and population structure would be expected to show greater temporal variation for successful introduced populations than for native populations. We studied the population genetics of the walnut husk fly, Rhagoletis completa, which was introduced into California from the Midwestern US in the early 1900s. We used microsatellites and allozymes to genotype current and historic fly populations, providing a rare perspective on temporal variability in population genetic parameters. We found that introduced populations showed greater temporal fluctuations in allele frequencies than native populations. Some introduced populations also showed an increase in genetic diversity over time, indicating multiple introductions had occurred. Population genetic structure decreased in both native and introduced populations over time. Our study demonstrates that introduced species are not at equilibrium and post-colonization processes may be important in ameliorating the loss of genetic diversity associated with biological invasions.     

Sampling for Microsatellite-Based Population Genetic Studies: 25 to 30 Individuals per Population Is Enough to Accurately Estimate Allele Frequencies

One of the most common questions asked before starting a new population genetic study using microsatellite allele frequencies is “how many individuals do I need to sample from each population?” This question has previously been answered by addressing how many individuals are needed to detect all of the alleles present in a population (i.e. rarefaction based analyses). However, we argue that obtaining accurate allele frequencies and accurate estimates of diversity are much more important than detecting all of the alleles, given that very rare alleles (i.e. new mutations) are not very informative for assessing genetic diversity within a population or genetic structure among populations. Here we present a comparison of allele frequencies, expected heterozygosities and genetic distances between real and simulated populations by randomly subsampling 5–100 individuals from four empirical microsatellite genotype datasets (Formica lugubris, Sciurus vulgaris, Thalassarche melanophris, and Himantopus novaezelandia) to create 100 replicate datasets at each sample size. Despite differences in taxon (two birds, one mammal, one insect), population size, number of loci and polymorphism across loci, the degree of differences between simulated and empirical dataset allele frequencies, expected heterozygosities and pairwise F_ST values were almost identical among the four datasets at each sample size. Variability in allele frequency and expected heterozygosity among replicates decreased with increasing sample size, but these decreases were minimal above sample sizes of 25 to 30. Therefore, there appears to be little benefit in sampling more than 25 to 30 individuals per population for population genetic studies based on microsatellite allele frequencies.     

Conservation genetics and phylogeography of southern dunlins Calidris alpina schinzii

Breeding populations of southern dunlin Calidris alpina schinzii in South Fennoscandia and the Baltic are severely fragmented and declining dramatically. Information on the genetic structure and diversity is therefore of importance for the conservation and management of these populations. Here we present the results of comparative genetic analyses of these populations with other populations of the schinzii , alpina and arctica subspecies in northern Europe. We sequenced the mitochondrial DNA control region and the Z-chromosome intron VLDLR-9, and analyzed microsatellites and AFLPs, for analyses of within-population genetic diversity. We also extended previous analyses of the phylogeographic structure of dunlins in northern Europe with a larger sample of individuals and populations. Our results revealed no evidence of reduced genetic diversity or increased levels of inbreeding in the small and fragmented populations around the Baltic Sea as compared to the more vital and larger populations elsewhere. Nevertheless, their small population sizes and presumably high degree of isolation may lead to local extinctions, indicating that demographic and ecological factors may pose a greater threat to the survival of these populations than purely genetic factors. Phylogeographically, the schinzii populations in Scandinavia and the Baltic do not form a separate genetic clade, but are part of larger cline of genetic variation encompassing several recognized subspecies of dunlins in the western Palearctic region. Only the Icelandic population showed some distinctiveness in genetic structure and might therefore be considered a separate management unit. Our study highlights the general problem of lack of genetic support for subspecies in avian taxonomy and conservation genetics.     

Molecular evolution and population genetics of Greater Caribbean green turtles (Chelonia mydas) as inferred from mitochondrial DNA control region sequences

The molecular evolution and population genetics of migratory green turtles (Chelonia mydas) in the Greater Caribbean were examined with mitochondrial DNA (mtDNA) control region I sequences. A total of 488 base positions (bp) per individual were aligned for 44 individuals from four nesting populations in Florida, Costa Rica, Aves Island (Venezuela), and Surinam. Twelve sequence polymorphisms were detected, representing ten transitions, one transversion, and one 10-bp repeat. Sequence analyses of within- and between-population diversity revealed a deep divergence between western and eastern Caribbean nesting colonies and an inverse relationship between reproductive female population size and mtDNA diversity. In small populations, genetic admixture was important to maintaining high diversity, whereas larger populations appear to have experienced historical bottlenecks or resulted from founder effects. Mitochondrial DNA sequences of the control region offer an order of magnitude greater resolution than restriction site data for addressing questions about mtDNA variation, both within and between populations of green turtles.     

Genetic diversity and condition factor: a significant relationship in Flemish but not in German populations of the European bullhead (Cottus gobio L.)

The population genetics of the Antarctic neritic krill species Euphausia crystallorophias was examined by nucleotide sequence variation in its mitochondrial DNA. A 616 base pair region of the cytochrome c oxidase subunit I (COI) gene was screened for mutations by single-strand conformational polymorphism (SSCP) combined with restriction digestion. E. crystallorophias caught in three different regions of the Antarctic coastline were used--two samples from the Mertz Glacier Polynya and one sample each from the western side of the Antarctic Peninsula and from the Davis Sea. Significant genetic differences between krill samples were identified. However, the extent of these differences did not correlate with the degree of geographic separation between the sampling sites. This suggests that the genetic structuring may be the result of small-scale differentiation rather than differentiation between resident populations in separate parts of the Southern Ocean. The possibility that genetic differences between samples within a region are as important as differences between regions has implications for other studies of krill population genetics.     

Analysis of genetic diversity in red clover (Trifolium pratense L.) breeding populations as revealed by RAPD genetic markers.

Red clover is an important forage legume species for temperate regions and very little is known about the genetic organization of its breeding populations. We used random amplified polymorphic DNA (RAPD) genetic markers to address the genetic diversity and the distribution of variation in 20 breeding populations and cultivars from Chile, Argentina, Uruguay, and Switzerland. Genetic distances were calculated for all possible pairwise combinations. A high level of polymorphism was found and the proportion of polymorphic loci across populations was 74.2%. A population derived from a non-certified seedlot displayed a higher proportion of polymorphic loci than its respective certified seedlot. Gene diversity values and population genetics parameters suggest that the populations analyzed are diverse. An analysis of molecular variance (AMOVA) revealed that the largest proportion of variation (80.4%) resides at the within population level. RAPD markers are a useful tool for red clover breeding programs. A dendrogram based on genetic distances divided the breeding populations analyzed into three distinct groups. The amount and partition of diversity observed can be of value in identifying the populations that parents of synthetic cultivars are derived from and to exploit the variation available in the populations analyzed.     

Genome-enabled hitchhiking mapping identifies QTLs for stress resistance in natural Drosophila

Identification of genes underlying complex traits is an important problem. Quantitative trait loci (QTL) are mapped using marker-trait co-segregation in large panels of recombinant genotypes. Most frequently, recombinant inbred lines derived from two isogenic parents are used. Segregation patterns are also studied in pedigrees from multiple families. Great advances have been made through creative use of these techniques, but narrow sampling and inadequate power represent strong limitations. Here, we propose an approach combining the strengths of both techniques. We established a mapping population from a sample of natural genotypes, and applied artificial selection for a complex character. Selection changed the frequencies of alleles in QTLs contributing to the selection response. We infer QTLs with dense genotyping microarrays by identifying blocks of linked markers undergoing selective changes in allele frequency. We demonstrated this approach with an experimental population composed from 20 isogenic strains. Selection for starvation survival was executed in three replicated populations with three control non-selected populations. Three individuals per population were genotyped using Affymetrix GeneChips. Two regions of the genome, one each on the left arms of the second and third chromosomes, showed significant divergence between control and selected populations. For the former region, we inferred allele frequencies in selected and control populations by pyrosequencing. We conclude that the allele frequency difference, averaging approximately 40% between selected and control lines, contributed to selection response. Our approach can contribute to the fine scale decomposition of the genetics of direct and indirect selection responses, and genotype by environment interactions.     

Genetic diversity and population structure of Aeromonas hydrophila, Aer. bestiarum, Aer. salmonicida and Aer. popoffii by multilocus enzyme electrophoresis (MLEE)

Genetic diversity, genetic relationship, identification and population structure of 120 Aeromonas strains (including Aer. hydrophila, Aer. bestiarum, Aer. salmonicida and Aer. popoffii) isolated from various sources were studied by analysis of 15 genetic loci by multilocus enzyme electrophoresis (MLEE). All 15 loci were polymorphic, with an average of 9.4 alleles per locus and a mean genetic diversity (H) of 0.64. Cluster analysis defined at H < or = 0.7 differentiated most of the taxa analysed except the Aer. popoffii and Aer. bestiarum strains, which showed a close genetic relationship. Allelic frequencies of five loci (EST1, HEX, IDH, LDH1 and MDH) identified 94% of the strains. The index of association (IA) for the total sample was 2.38 and IA values calculated for the different populations were always significantly different from zero. These results suggest that the population structure of this Aeromonas sample is strongly clonal, confirm the taxonomic status of the analysed species in population genetics terms, and show the usefulness of MLEE for identifying Aeromonas species.     

DNA sequence variation in a 3.7-kb noncoding sequence 5' of the CYP1A2 gene: implications for human population history and natural selection

CYP1A2 is a cytochrome P450 gene that is involved in human physiological responses to a variety of drugs and toxins. To investigate the role of population history and natural selection in shaping genetic diversity in CYP1A2, we sequenced a 3.7-kb region 5' from CYP1A2 in a diverse collection of 113 individuals from three major continental regions of the Old World (Africa, Asia, and Europe). We also examined sequences in the 90-member National Institutes of Health DNA Polymorphism Discovery Resource (PDR). Eighteen single-nucleotide polymorphisms (SNPs) were found. Most of the high-frequency SNPs found in the Old World sample were also found in the PDR sample. However, six SNPs were detected in the Old World sample but not in the PDR sample, and two SNPs found in the PDR sample were not found in the Old World sample. Most pairs of SNPs were in complete linkage disequilibrium with one another, and there was no indication of a decline of disequilibrium with physical distance in this region. The average +/- SD nucleotide diversity in the Old World sample was 0.00043+/-0.00026. The African population had the highest level of nucleotide diversity and the lowest level of linkage disequilibrium. Two distinct haplotype clusters with broadly overlapping geographical distributions were present. Of the 17 haplotypes found in the Old World sample, 12 were found in the African sample, 8 were found in Indians, 5 were found in non-Indian Asians, and 5 were found in Europeans. Haplotypes found outside Africa were mostly a subset of those found within Africa. These patterns are all consistent with an African origin of modern humans. Seven SNPs were singletons, and the site-frequency spectrum showed a significant departure from neutral expectations, suggesting population expansion and/or natural selection. Comparison with outgroup species showed that four derived SNPs have achieved high (>0.90) frequencies in human populations, a trend consistent with the action of positive natural selection. These patterns have a number of implications for disease-association studies in CYP1A2 and other genes.     

遗传多样性的取样策略

合理取样是生物多样性有效保护、利用和研究所面临的最基本问题 ,它在很大程度上受到植物自身的生物学特性、环境条件和取样目的的影响。遗传多样性的取样策略是指对一定地理分布范围内的生物个体取样时 ,使样本具有代表性和包含尽可能多的遗传变异的最佳取样方法 ,包括了取样数目 (一个给定区域的居群数和一个居群的个体数 )以及取样方式。包括“哈迪 温伯格平衡 (Hardy WeinbergEquilibrium)”定律在内的居群遗传学基本原理是研究取样策略的理论基础 ,在此基础上可以对居群内的取样个体数及应获取的居群数进行理论计算 ,同时还可以根据物种居群的遗传结构特点和环境条件的异质性来决定取样的方式。因此 ,应该依据研究对象本身的特点和取样的目的来确定某一特定区域的居群取样数 ,以及某一居群内的样本数及取样方式。     

A comparison of single‐sample effective size estimators using empirical toad (Bufo calamita) population data: genetic compensation and population size‐genetic diversity correlations

The accuracy and precision of four single‐sample estimators of effective population size, N_e (heterozygote excess, linkage disequilibrium, Bayesian partial likelihood and sibship analysis) were compared using empirical data (microsatellite genotypes) from multiple natterjack toad (Bufo calamita) populations in Britain (n = 16) and elsewhere in Europe (n = 10). Census size data were available for the British populations. Because toads have overlapping generations, all of these methods estimated the number of effective breeders N_b rather than N_e. The heterozygote excess method only provided results, without confidence limits, for nine of the British populations. Linkage disequilibrium gave estimates for 10 British populations, but only six had finite confidence limits. The Bayesian and sibship methods both produced estimates with finite confidence limits for all the populations. Although the Bayesian method was the most precise, on most criteria (insensitivity to locus number, correlation with other effective and census size estimates and correlation with genetic diversity) the sibship method performed best. The results also provided evidence of genetic compensation in natterjack toads, and highlighted how the relationship between effective size and genetic diversity can vary as a function of geographical scale.