Hitchhiking effect of a beneficial mutation spreading in a subdivided population

A central problem in population genetics is to detect and analyze positive natural selection by which beneficial mutations are driven to fixation. The hitchhiking effect of a rapidly spreading beneficial mutation, which results in local removal of standing genetic variation, allows such an analysis using DNA sequence polymorphism. However, the current mathematical theory that predicts the pattern of genetic hitchhiking relies on the assumption that a beneficial mutation increases to a high frequency in a single random-mating population, which is certainly violated in reality. Individuals in natural populations are distributed over a geographic space. The spread of a beneficial allele can be delayed by limited migration of individuals over the space and its hitchhiking effect can also be affected. To study this effect of geographic structure on genetic hitchhiking, we analyze a simple model of directional selection in a subdivided population. In contrast to previous studies on hitchhiking in subdivided populations, we mainly investigate the range of sufficiently high migration rates that would homogenize genetic variation at neutral loci. We provide a heuristic mathematical analysis that describes how the genealogical structure at a neutral locus linked to the locus under selection is expected to change in a population divided into two demes. Our results indicate that the overall strength of genetic hitchhiking--the degree to which expected heterozygosity decreases--is diminished by population subdivision, mainly because opportunity for the breakdown of hitchhiking by recombination increases as the spread of the beneficial mutation across demes is delayed when migration rate is much smaller than the strength of selection. Furthermore, the amount of genetic variation after a selective sweep is expected to be unequal over demes: a greater reduction in expected heterozygosity occurs in the subpopulation from which the beneficial mutation originates than in its neighboring subpopulations. This raises a possibility of detecting a "hidden" geographic structure of population by carefully analyzing the pattern of a selective sweep.     

DNA Polymorphism in a Subdivided Population: The Expected Number of Segregating Sites in the Two-Subpopulation Model

Using the two subpopulation model, the expected numbers of segregating sites in a number of DNA sequences randomly sampled from a subdivided population were examined for several types of population subdivisions. It is shown that, in the case where the pattern of migration is symmetrical such as the finite island model, the expected number of segregating sites is independent of the migration rate when two or three DNA sequences are randomly sampled from the same subpopulation, but depends on the migration rate when more than three DNA sequences are sampled. It is also shown that the population subdivision can increase the amount of DNA polymorphism even in a subpopulation in some cases.     

Genetic Divergence and Signatures of Natural Selection in Marginal Populations of a Keystone,Long-Lived Conifer,Eastern White Pine (Pinus strobus) from Northern Ontario

Marginal populations are expected to provide the frontiers for adaptation, evolution and range shifts of plant species under the anticipated climate change conditions. Marginal populations are predicted to show genetic divergence from central populations due to their isolation, and divergent natural selection and genetic drift operating therein. Marginal populations are also expected to have lower genetic diversity and effective population size (N _e) and higher genetic differentiation than central populations. We tested these hypotheses using eastern white pine (Pinus strobus) as a model for keystone, long-lived widely-distributed plants. All 614 eastern white pine trees, in a complete census of two populations each of marginal old-growth, central old-growth, and central second-growth, were genotyped at 11 microsatellite loci. The central populations had significantly higher allelic and genotypic diversity, latent genetic potential (LGP) and N _e than the marginal populations. However, heterozygosity and fixation index were similar between them. The marginal populations were genetically diverged from the central populations. Model testing suggested predominant north to south gene flow in the study area with curtailed gene flow to northern marginal populations. Signatures of natural selection were detected at three loci in the marginal populations; two showing divergent selection with directional change in allele frequencies, and one balancing selection. Contrary to the general belief, no significant differences were observed in genetic diversity, differentiation, LGP, and N _e between old-growth and second-growth populations. Our study provides information on the dynamics of migration, genetic drift and selection in central versus marginal populations of a keystone long-lived plant species and has broad evolutionary, conservation and adaptation significance.     

Patterns of Gene Variation in Central and Marginal Populations of DROSOPHILA ROBUSTA

The central and marginal populations of D. robusta differ greatly in the level of inversion polymorphism; the marginal populations are monomorphic or nearly so and the central populations are highly polymorphic. This paper presents the frequencies of alleles at forty gene loci in various populations of D. robusta, studied by electrophoresis of proteins and enzymes. Population samples were obtained from eight widely separated populations of D. robusta which included the central, the extreme marginal and the intervening populations between the center and the margins. We find that the proportion of polymorphic loci and average heterozygosity per individual is slightly higher in the marginal populations than the central populations. In D. robusta on an average, 39% of the loci are polymorphic and the average proportion of loci heterozygous per individual is 11%. A breakdown of loci in three categories, viz, hydrolytic enzymes and some other enzymes, larval proteins and glycolytic and Kreb's cycle enzymes, shows that in all populations the level of polymorphism is highest in the hydrolytic enzymes, intermediate in larval proteins and least in the glycolytic and Kreb's cycle enzymes. On the average, the proportion of loci heterozygous per individual for three groups of loci is: hydrolytic enzymes and others (.164), larval proteins (.115) and glycolytic and Kreb's cycle enzymes (.037). We also observe that in all populations the level of polymorphism on the X chromosome is far less than the expected 38%; in salivary gland cells the euchromatic length of the X chromosome is 38% of the entire genome. Lower levels of polymorphism for the X chromosome loci are explained due to low probability of balanced polymorphisms for the X-linked loci since the conditions for establishment of balanced polymorphism for X-linked loci are more restrictive than for the autosomal loci.-The polymorphic loci can be grouped according to pattern of allele frequencies in different populations as follows: (1) The allele frequencies are similar in all populations at the XDH, Pep-1 and Hex-1 loci. (2) The alleles at the Est-1, Est-2, Amy loci and the AP-4(1.0) and the LAP-1(.90) alleles show north south clinal change in frequency. (3) There is north south and east west differentiation at the Pt-5, Pt-8 and Pt-9 loci and the allele AP-4(.81). (4) Polymorphism at loci such as Fum, B.Ox, Hex-8, Pep-2 and Pep-3 are restricted to only one or two of the populations. (5) Allele frequencies at the MDH and ODH loci fluctuate between populations. (6) Allele frequencies at many polymorphic loci such as Est-1, Est-2, LAP-1, AP-4, Pt-5, Pt-8, Pt-9, Pt-16, MDH, Fum change clinally within a gene arrangement. The pattern of gene variation in D. robusta is very complex and cannot be easily explained due to migration of neutral alleles between once-isolated populations or to semi-isolation of neutral alleles. The observations of the pattern of allele variation in different populations, high levels of polymorphism in the marginal populations which have small population size and low levels of polymorphism of the X chromosome loci all support the argument in favor of balancing selection as the main mechanism for the maintenance of these polymorphisms. Environmental factors must play a role in the maintenance of a great deal of these polymorphisms, since we observe clinal allele frequency changes even within a given inversion type.     

A representative model of the hierarchic structure of human populations]

The paper deals with the concept of a hierarchy population. The population is hierarchial, if it can be divided into a certain amount of subpopulations in such a way that the set may naturally break into classes (levels). A migration may exist from each sub-population or into a higher one. Such population structures are frequently encountered among the human populations. It is shown that in a hierarchial population a selection is more efficient than in a non-subdivided one. In this respect the role of small populations in evolution is suggested. It was possible to demonstrate that in a hierarchial population there was a higher degree of polymorphism, a higher velocity of the evolutional process, that in a non-subdivided population of the same size. Besides, the level of polymorphism in the hierarchial population increases monotomously with the growth of the amount of generations during rather a long period of time.     

Population genetic structure and differentiation of Anthoxanthum alpinum in the subalpine-alpine ecocline of Swiss Alps

Allozyme variation in 6 enzyme systems coding 10 loci was estimated for 18 subpopulations of Anthoxanthum alpinum from three altitudinal transects in two localities of the Swiss Alps. Mean proportions of polymorphic loci (95% criterion), average number of alleles per locus, and mean expected heterozygosity were 64.9%, 2.37 and 0. 252, respectively. Mean genetic distance among populations was 0.011, and 79% of the genetic variation resided within populations. Based on allozyme analysis, marginal subpopulations appeared to have similar level of genetic variability to central subpopulations. Relatively high genetic differentiation, low gene flow values and small neighbourhood sizes suggested that inbreeding followed by genetic drift was possible causes of lowgenetic variability in Arpette A. alpinum populations.     

阿尔卑斯山高山-亚高山过渡区高山黄花茅的群体遗传结构和分化研究

应用等位酶分析,在瑞士阿尔卑斯山的阿尔拜特(Arpette)和拜阿尔普(Belalp),沿三个不同的海拔梯度,研究了高山黄花茅3个自然居群的遗传变异和分化。研究结果表明,平均的多态性位点比例为64．9％,每个位点平均等位基因数为2．37,平均期望杂合性为0．252。亚居群间平均的遗传距离为0．011,发现79％的遗传变异存在于居群内。基于等位酶分析,边缘亚居群与中心亚居群似乎有类似的遗传变异性。相对比较高的遗传分化、低的亚居群间基因流和小的邻居大小值暗示,近交和随后的遗传漂变可能是导致阿尔拜特黄花茅居群遗传变异性较低的主要原因。     

Subdivision and haplotype structure in natural populations of Arabidopsis lyrata

High levels of inbreeding are expected to cause a strong reduction in levels of genetic variability, effective recombination rates and in adaptation compared with related outcrossing populations. We examined patterns of DNA polymorphism at five nuclear loci and one chloroplast locus within and between four populations of the outcrossing plant Arabidopsis lyrata, a close relative of the highly self-fertilizing model species A. thaliana. The observed patterns are compared with species-wide polymorphism at orthologous loci, as well as within- and between-population patterns at other studied loci in A. thaliana. In addition to evidence for much higher average within-population diversity, species-wide levels of silent polymorphism are generally higher in A. lyrata than in A. thaliana, unlike the results from a previous study of the ADH locus. However, polymorphism is also low in the North American A. lyrata subspecies lyrata compared with the European subspecies petraea, suggesting either a population bottleneck in North American populations or recent admixture involving diverged European populations. Differentiation between the two subspecies is strong, although there are few fixed differences, suggesting that their isolation is recent. Estimates of intralocus recombination rates and analysis of haplotype structure in European A. lyrata populations indicate lower recombination than predicted based on the variability together with physical recombination rates estimated from A. thaliana. This may be due to strong population subdivision, or to recent departures from demographic equilibrium such as a bottleneck or population admixture. Alternatively, there may be consistently lower recombination rates in the outcrossing species. In contrast, estimates of recombination rates from species-wide samples of A. thaliana are close to the values expected assuming a high rate of self-fertilization. Complex population histories in both A. thaliana and A. lyrata complicate theoretical predictions and empirical tests of the effects of inbreeding on polymorphism and molecular evolution.     

High levels of genetic diversity in marginal populations of the marsh orchid Dactylorhiza majalis subsp. majalis

Because of harsh conditions, suboptimal habitat quality and poor connectivity to other populations, plant populations at the margin of a distribution are expected to be less genetically diverse, but to be more divergent from each other than populations in the centre of a distribution. In northern Europe, northern marginal populations may also be younger than populations further to the south, and may have had less time to accumulate genetic diversity by mutation and gene flow. However, orchids have very small seeds, which are easily dispersed long distances by wind, and orchids are therefore expected to show less differentiation between marginal and central populations than other groups of seed plants. Here, we analysed whether Scandinavian populations of the tetraploid marsh orchid Dactylorhiza majalis subsp. majalis differ from central European populations in genetic diversity patterns. A total of 220 plants from eight central European and ten Scandinavian populations was examined for variation at five nuclear microsatellite loci, nuclear ITS and 13 polymorphic sites in noncoding regions of the plastid genome. The total genetic diversity was slightly lower in Scandinavia than in central Europe, both in plastid and nuclear markers, but the differences were small. Also, the Scandinavian populations were less diverse and somewhat more strongly differentiated from each other than the central European ones. Dactylorhiza majalis subsp. majalis has apparently colonized Scandinavia on multiple independent occasions and from different source areas in the south. Seed flow between Scandinavian populations has still not fully erased the patterns imprinted by early colonization. Our results suggest that marginal populations of orchids may be as important as central ones in preserving genetic diversity through Pleistocene glacial cycles. We also predict that orchids with their light seeds are better adapted than many other plants to respond to future climate changes by dispersing into new suitable areas.     

Effective population size,genetic diversity,and coalescence time in subdivided populations

A formula for the effective population size for the finite island model of subdivided populations is derived. The formula indicates that the effective size can be substantially greater than the actual number of individuals in the entire population when the migration rate among subpopulations is small. It is shown that the mean nucleotide diversity, coalescence time, and heterozygosity for genes sampled from the entire population can be predicted fairly well from the theory for randomly mating populations if the effective population size for the finite island model is used.     

Effect of population structure on the amount of polymorphism and the fixation probability under overdominant selection

Under overdominant selection, mutants substantially contribute to increase the amount of polymorphism. It is also known that under neutrality as the migration rates among demes decrease in a subdivided population, the amount of polymorphism increases along with the increase of the effective population size, N(e). In this study, under overdominant selection the effect of population subdivision on the amount of polymorphism was investigated using the diffusion approximation and the low migration approximation. It was shown that if selection is medium or strong (e.g., N(T)s > 1, where N(T) is the population size and s is the selective advantage of heterozygotes), the nucleotide diversity, pi, decreases along with the decrease of Nm against the increase of N(e), where N is the size of demes and m is the migration rate per deme. In addition, the ratio of the nucleotide diversity to the evolutionary rate also decreases along with the decrease of Nm. In some cases the ratio becomes smaller than that expected under neutrality as Nm decreases.     

Genetic Load and Viability Distribution in Central and Marginal Populations of DROSOPHILA SUBOBSCURA

Relative viabilities of individuals homozygous or randomly heterozygous for wild O chromosomes derived from a marginal (Norwegian) and a central (Greek) population of D. subobscura were obtained by means of a newly prepared marker strain. In the central and marginal populations 20.8 and 28.8 percent of all chromosomes proved lethal or semilethal in homozygous condition. Mean viability was higher for +/+ random heterozygotes than for +/+ homozygotes. This remained the case for the marginal, but not for the central populations, after exclusion of the detrimental chromosomes from the calculations. The variances of viabilities were higher for homozygotes than for heterozygotes, but the test crosses with chromosomes from the marginal population had generally higher variances than those with chromosomes from the central population. No correlation was found in either populations between the action of a chromosome in homozygous condition with its action in heterozygous condition. This is interpreted as complete recessiveness of genetic load. The results are discussed in terms of the observed reduction of the inversion polymorphism which is not paralleled by a reduction in enzyme and, as shown here, by reduction in viability variation. It is thought that the heterotic effect of inversions is due to their homeostatic action, which depends less on structural genes than on higher orders of organization due to gene interaction or regulation. Whatever the causes, it is very likely that marginal populations differ from central populations with respect to their genetic system.     

Sex-biased evolutionary forces shape genomic patterns of human diversity

Comparisons of levels of variability on the autosomes and X chromosome can be used to test hypotheses about factors influencing patterns of genomic variation. While a tremendous amount of nucleotide sequence data from across the genome is now available for multiple human populations, there has been no systematic effort to examine relative levels of neutral polymorphism on the X chromosome versus autosomes. We analyzed ~210 kb of DNA sequencing data representing 40 independent noncoding regions on the autosomes and X chromosome from each of 90 humans from six geographically diverse populations. We correct for differences in mutation rates between males and females by considering the ratio of within-human diversity to human-orangutan divergence. We find that relative levels of genetic variation are higher than expected on the X chromosome in all six human populations. We test a number of alternative hypotheses to explain the excess polymorphism on the X chromosome, including models of background selection, changes in population size, and sex-specific migration in a structured population. While each of these processes may have a small effect on the relative ratio of X-linked to autosomal diversity, our results point to a systematic difference between the sexes in the variance in reproductive success; namely, the widespread effects of polygyny in human populations. We conclude that factors leading to a lower male versus female effective population size must be considered as important demographic variables in efforts to construct models of human demographic history and for understanding the forces shaping patterns of human genomic variability.     

Quantifying the Elevation of Mitochondrial DNA Evolutionary Substitution Rates Over Nuclear Rates in the Intertidal Copepod Tigriopus californicus

Mitochondrial DNA (mtDNA) genomes generally evolve rapidly in animals, but considerable variation in the rates of evolution of mtDNA occurs among taxa. Higher levels of mutation will tend to increase the amount of polymorphism, which should also scale with population size, but there are mixed signals from previous studies on the evolutionary outcomes of the interactions of these processes. The copepod Tigriopus californicus provides an interesting model in which to study the evolution of mtDNA because it has high levels of divergence among populations and there is the suggestion that this divergence could be involved in reproductive isolation. This species also appears to have an elevated mtDNA substitution rate, but previous studies did not provide an accurate measurement. This article examines the rate of mtDNA substitution versus nuclear substitution in T. californicus and finds that the mtDNA rate for synonymous sites averages 55-fold higher, a level that exceeds the rates found in most other invertebrates. Levels of polymorphism are also examined in both mtDNA and nuclear genes, and it is shown that the effective population size of mtDNA genes is much lower than that of nuclear genes. In addition, no correlation between polymorphism in mtDNA and nuclear genes is found across populations, which suggest factors other than demography may shape polymorphism in this species. The results from this study suggest that mtDNA is evolving at a very rapid rate in this copepod species, and this could increase the likelihood that mtDNA evolution is involved in the generation of reproductive isolation.     

Seeing ghosts: the effect of unsampled populations on migration rates estimated for sampled populations

In 2004, the term 'ghost population' was introduced to summarize the effect of unsampled subpopulations that exchange migrants with other subpopulations that have been sampled. Estimated long-term migration rates among populations sampled will be affected by ghost populations. Although it would be convenient to be able to define an apparent migration matrix among sampled populations that incorporate the exchange of migrants with ghost populations, no such matrix can be defined in a way that predicts all features of the coalescent process for the true migration matrix. This paper shows that if the underlying migration matrix is symmetric, it is possible to define an apparent migration matrix among sampled subpopulations that predicts the same within-population and between-population homozygosities among sampled populations as is predicted by the true migration matrix. Application of this method shows that there is no simple relationship between true and apparent migration rates, nor is there a way to place an upper bound on the effect of ghost populations. In general, ghost populations can create the appearance of migration between subpopulations that do not actually exchange migrants. Comparison with published results from the application of the program, MIGRATE, shows that the apparent migration rates inferred with that program in a three-subpopulation model differ from those based on pairwise homozygosities. The apparent migration matrix determined by the method described in this paper probably represents the upper bound on the effect of ghost populations.     

Patterns of parapatric speciation

Abstract. Geographic variation may ultimately lead to the splitting of a subdivided population into reproductively isolated units in spite of migration. Here, we consider how the waiting time until the first split and its location depend on different evolutionary factors including mutation, migration, random genetic drift, genetic architecture, and the geometric structure of the habitat. We perform large-scale, individual-based simulations using a simple model of reproductive isolation based on a classical view that reproductive isolation evolves as a by-product of genetic divergence. We show that rapid parapatric speciation on the time scale of a few hundred to a few thousand generations is plausible even when neighboring subpopulations exchange several individuals each generation. Divergent selection for local adaptation is not required for rapid speciation. Our results substantiates the claims that species with smaller range sizes (which are characterized by smaller local densities and reduced dispersal ability) should have higher speciation rates. If mutation rate is small, local abundances are low, or substantial genetic changes are required for reproductive isolation, then central populations should be the place where most splits take place. With high mutation rates, high local densities, or with moderate genetic changes sufficient for reproductive isolation, speciation events are expected to involve mainly peripheral populations.     

Relationship between DNA Polymorphism and Fixation Time

When there is no recombination among nucleotide sites in DNA sequences, DNA polymorphism and fixation of mutants at nucleotide sites are mutually related. Using the method of gene genealogy, the relationship between the DNA polymorphism and the fixation of mutant nucleotide was quantitatively investigated under the assumption that mutants are selectively neutral, that there is no recombination among nucleotide sites, and that the population is a random mating population with N diploid individuals. The results obtained indicate that the expected number of nucleotide differences between two DNA sequences randomly sampled from the population is 42% less when a mutant at a particular nucleotide site reaches fixation than at a random time, and that heterozygosity is also expected to be less when fixation takes place than at a random time, but the amount of reduction depends on the value of 4Nv in this case, where v is the mutation rate per DNA sequence per generation. The formula for obtaining the expected number of nucleotide differences between the two DNA sequences for a given fixation time is also derived, and indicates that, even when it takes a large number of generations for a mutant to reach fixation, this number is 33% less than at a random time. The computer simulation conducted suggests that the expected number of nucleotide differences between the two DNA sequences at the time when an advantageous mutant becomes fixed is essentially the same as that of neutral mutant if the fixation time is the same. The effect of recombination on the amount of DNA polymorphism was also investigated by using computer simulation.     

Some consideration on diversifying selection

The diversifying selection due to genotype-environment interaction can increase the genetic variation in natural populations. It is known, however, that the conditions for stable genetic polymorphism or marginal overdominance are quite restricted in this selection model. In this paper a simple model of diversifying selection was examined, and the following results were obtained: (1) Even when the conditions for marginal overdominance are not satisfied, if the diversifying selection is operating, the frequency of mutants can be higher than that in the case of simple mutation-selection balance. (2) This selection model causes a large amount of genetic load (environment load), even when the conditions for marginal overdominance are not satisfied, namely even when the equilibrium frequency of mutant is very low. From these results it can be concluded that the number of loci on which this type of diversifying selection is operating is very small, if any.     

Gene admixture in the silk road region of China: evidence from mtDNA and melanocortin 1 receptor polymorphism

Mitochondrial DNA control region segment I sequences and melanocortin 1 receptor (MC1R) gene polymorphism were examined in ethnic populations in the silk road region of China. Both the frequencies of the MC1R variants and the results of mtDNA data in this region presented intermediate values between those of Europe and East and Southeast Asia, which suggested extensive gene admixture in this area and was in general agreement with previous studies. Phylogenetic analysis of the ethnic populations in the Silk Road region that based on mtDNA data didn't show expected cluster pattern according to their ethnogenesis. We suspect that a high migration rate in female among these closely related populations and other three demographic events might account for it.     

Low isozyme diversity in Nordic and central European populations of Vicia pisiformis and V. dumetorum (Fabaceae)

Vicia pisiformis and V. dumetorum are closely related perennial species occurring in small, disjunct populations throughout Europe. In Sweden, they are classified as "vulnerable" and "care-demanding", respectively. They possess low variation in RAPD (Random Amplified Polymorphic DNA) markers but relatively high morphological variation. We estimated the amount and distribution of isozyme diversity within and among marginal (Sweden, Norway) and central (Czech Republic) populations. Vicia dumetorum was monomorphic for all 14 isozyme loci in the 22 populations investigated, whereas V. pisiformis had 2.94 % polymorphic loci and a mean of 1.02 alleles per loci. The differentiation (F_sr) of populations of V pisiformis were 0.848, but no clear geographical pattern was found. Preliminary results from artificial pollination experiments suggests self-compatibility in V. dumetorum . The low overall polymorphism in Swedish and Czech localities in both species is explained in terms of common postglacial colonization events, recent anthropogenic activity, selfing and clonal reproduction. Finally, the implications of the data for conservation genetics is discussed.