Natural selection opposing artificial selection: a two-locus deterministic simulation

A two-locus, two-allele metric trait was submitted to artificial truncation selection and to three types of opposing natural selection (two-locus extensions of directional selection, overdominance and underdominance) by numerical simulation in a large random-mating population. Limits to selection were generally reached by generation 100. Intermediate selection plateaus were found, with minor genes, for all three modes of opposing natural selection, but they were least frequent with underdominance. Multiple outcomes were common. In particular, fixation of the genotype favored by artificial selection was often associated with fixation of another genotype and/or with a central equilibrium; the end point actually reached depended on the genetic starting point of the simulation. In general, when the alleles favored by truncation selection were combined (positive linkage disequilibrium) in the base population, or when the trait was determined by major genes, artificial selection would prevail. Limitations inherent to this type of work are discussed, and possible avenues for further work on the antagonism between artificial and natural selection are proposed.     

Mutation-Selection Balance with Stochastic Selection

Diffusion theory has been used to analyze a model of mutation-selection balance in which the selection process is assumed to be stochastic in time. The limiting outcome of the mutation-stochastic selection process is determined qualitatively by the geometric mean fitnesses of the genotypes, and the conditions for fixation or polymorphism are similar to those that determine the outcome of the mutation-selection process when selection is constant. However, in the case of a completely recessive allele, detailed numerical study of the polymorphism associated with stochastic selection has shown that the average allele frequency maintained is greater than the equilibrium frequency expected when selection is constant, even when the geometric mean fitness of the recessive homozygotes is identical in the stochastic and deterministic models. Thus, allele frequencies in natural populations that are too high to be plausibly explained by a balance between mutation and constant selection can be accounted for if selection is stochastic.     

不同QTL增效基因初始频率下标记辅助选择的效果

采用随机模拟方法模拟了在一个闭锁群体内连续对单个性状选择10个世代的情形。在假定选择性状受一个位于常染色体上的QTL和多基因共同控制的情况下,采用动物模型标记辅助最佳线性无偏预测方法估计个体育种值并据此进行种畜的选留,并在此基础上系统地比较了QTL增效基因初始频率对标记辅助选择效果的影响。结果表明：当群体中QTL增效基因的初始频率较低时,选择所获得的QTL基因型值的进展会更大,标记辅助选择在单位时间内可获得较大的遗传进展;此时,尽管QTL增效基因在群体中固定所需的世代数会更长一些,但其频率上升的速度却更快。而QTL增效基因初始频率的高低对群体近交增量的影响不是很大。     

Reinforcement and the genetics of nonrandom mating

Abstract.— The occurrence of reinforcement is compared when premating isolation is caused by the spread of a gene causing females to prefer to mate with males carrying a population-specific trait (a "preference" model) and by a gene that causes females to prefer to mate with males that share their own trait phenotype (an "assortative mating" model). Both two-island models, which have symmetric gene flow, and continent-island models, which have one-way gene flow, are explored. Reinforcement is found to occur much more easily in a two-island assortative mating model than in any of the other three models. This is due primarily to the fact that in this model the assortative mating allele will automatically become genetically associated in each population with the trait allele that is favored by natural selection on that island. In contrast, natural selection on the trait both favors and opposes the evolution of premating isolation in the two-island preference model, depending on the particular population. These results imply that species recognition in the context of mating may evolve particularly easily when it targets cues that are favored by natural selection in each population. In the continent-island models, reinforcement is found to occur more often under the preference model than the assortative mating model, thus reversing the trend from the two-island models. Patterns of population subdivision may therefore play a role in determining what types of premating isolation may evolve.     

Neofunctionalization of duplicated genes under the pressure of gene conversion

Neofunctionalization occurs when a neofunctionalized allele is fixed in one of duplicated genes. This is a simple fixation process if duplicated genes accumulate mutations independently. However, the process is very complicated when duplicated genes undergo concerted evolution by gene conversion. Our simulations demonstrate that the process could be described with three distinct stages. First, a newly arisen neofunctionalized allele increases in frequency by selection, but gene conversion prevents its complete fixation. These two factors (selection and gene conversion) that work in opposite directions create an equilibrium, and the time during which the frequency of the neofunctionalized allele drifts around the equilibrium value is called the temporal equilibrium stage. During this temporal equilibrium stage, it is possible that gene conversion is inactivated by mutations, which allow the complete fixation of the neofunctionalized allele. And then, permanent neofunctionalization is achieved. This article develops basic population genetics theories on the process to permanent neofunctionalization under the pressure of gene conversion. We obtain the probability and time that the frequency of a newly arisen neofunctionalized allele reaches the equilibrium value. It is also found that during the temporal equilibrium stage, selection exhibits strong signature in the divergence in the DNA sequences between the duplicated genes. The spatial distribution of the divergence likely has a peak around the site targeted by selection. We provide an analytical expression of the pattern of divergence and apply it to the human red- and green-opsin genes. The theoretical prediction well fits the data when we assume that selection is operating for the two amino acid differences in exon 5, which are believed to account for the major part of the functional difference between the red and green opsins.     

Apparent Selection of Enzyme Alleles in Laboratory Populations of Drosophila

Twelve laboratory populations of recently collected Drosophila willistoni were begun with different frequencies of alleles at three enzyme loci, six populations at 25 degrees and six at 19 degrees . Periodic sampling of the populations allowed monitoring of the frequency changes in allozymes over time.-At Lap-5 (a locus coding for leucine amino peptidase), three alleles converged to the same frequencies in all populations at both temperatures. The apparent equilibrium frequency of the major allele was about.75; this is different from the frequency (.57) found in the natural population from which the experimental populations were begun. Allele frequency changes at the esterase-5 locus (Est-5) were slower but consistent in all cages. It is difficult to determine if an equilibrium has been reached. However, the frequency of the rare allele in all cages is about the same as in wild populations, 5%. Alleles at both Lap-5 and Est-5 are non-randomly associated with inversions in the chromosomes onto which they map. Because of these associations, it is impossible to unambiguously attribute the change in allele frequencies to selection at the loci being observed.-After one year, no significant gene frequency changes were detected at Est-7, the third locus studied.     

Simulating natural selection in landscape genetics

Linking landscape effects to key evolutionary processes through individual organism movement and natural selection is essential to provide a foundation for evolutionary landscape genetics. Of particular importance is determining how spatially-explicit, individual-based models differ from classic population genetics and evolutionary ecology models based on ideal panmictic populations in an allopatric setting in their predictions of population structure and frequency of fixation of adaptive alleles. We explore initial applications of a spatially-explicit, individual-based evolutionary landscape genetics program that incorporates all factors--mutation, gene flow, genetic drift and selection--that affect the frequency of an allele in a population. We incorporate natural selection by imposing differential survival rates defined by local relative fitness values on a landscape. Selection coefficients thus can vary not only for genotypes, but also in space as functions of local environmental variability. This simulator enables coupling of gene flow (governed by resistance surfaces), with natural selection (governed by selection surfaces). We validate the individual-based simulations under Wright-Fisher assumptions. We show that under isolation-by-distance processes, there are deviations in the rate of change and equilibrium values of allele frequency. The program provides a valuable tool (cdpop v1.0; http://cel.dbs.umt.edu/software/CDPOP/) for the study of evolutionary landscape genetics that allows explicit evaluation of the interactions between gene flow and selection in complex landscapes.     

杂交群体动物基因聚合

基因聚合是通过优化设计杂交方案,选择利用目标基因或与其紧密连锁的分子标记,通过世代选择实现将来源于多个不同群体的优势目标基因或基因型聚合到同一个理想个体中,进而达到生产出超级经济性状个体的目的。针对聚合不同目标基因个数,设计4类杂交方案——两群体、三群体、四群体级联、四群体对称。在相同的杂交方案中,比较基因型选择和表型选择策略,分析不同杂交组合、性状遗传力、初始基因频率、基础群体规模对聚合设计的影响,并筛选出最佳的聚合方案。研究结果表明,在较大的基础群体规模和较高初始群体基因频率下,获得聚合多个目标基因的理想个体的可能性较大。在四群体杂交方案中,亲本的杂交次序对于级联杂交比对称杂交的影响较大。模拟结果表明,运用基因型选择进行聚合育种优于表型选择。文章所开发设计的聚合模拟育种的统计分析方法和相应软件为指导杂交育种方案和选择策略的设计提供理论参考,同时,为进一步设计开发聚合设计模拟育种平台奠定基础。     

Interchromosomal Biased Gene Conversion, Mutation and Selection in a Multigene Family

A mathematical model of the effects of interchromosomal biased gene conversion, mutation and natural selection on a multigene family is developed and analyzed. The model assumes two allelic states at each of n loci. The effects of genetic drift are ignored. The model is developed under the assumption of no recombination, but the analysis shows that, at equilibrium, there is no linkage disequilibrium, which implies that the conclusions are valid for arbitrary recombination among loci. At equilibrium, the balance between mutation, gene conversion and selection depends on the ratio of the mutation rates to the quantity [s + g(2α - 1)/ n], where s is the increment or decrement in relative fitness with each additional copy of one of the alleles, g is the conversion rate, and α is a measure of the bias in favor of one of the alleles. When this quantity is large relative to the mutation rates, the allele that has the net advantage, combining the effects of selection and conversion, will be nearly fixed in the multigene family. A comparison of these results with those from a comparable model of intrachromosomal biased conversion shows that biased interchromosomal conversion leads to approximately the same equilibrium copy number as does intrachromosomal conversion of the same strength. Interchromosomal conversion is much more effective in causing the substitution of one allele by another. The relative frequencies of interchromosomal and intrachromosomal conversion is indicated by the extent of the linkage disequilibrium among the loci in a multigene family.     

Genetic Divergence under Uniform Selection. II. Different Responses to Selection for Knockdown Resistance to Ethanol among DROSOPHILA MELANOGASTER Populations and Their Replicate Lines

We have tested the hypothesis that genetic differences among conspecific populations may result in diverse responses to selection, using natural populations of Drosophila melanogaster. Selection for ethanol tolerance in a tube measuring knockdown resistance was imposed on five West Coast populations. In 24 generations the selected lines increased their mean knockdown times, on average, by a factor of 2.7. An initially weak latitudinal cline was steepened by selection. The two southernmost populations showed the same increases in the selected character, but differed consistently in their correlated responses in characters related to ethanol tolerance. This result indicates that the populations responded to selection by different genetic changes. Selection decreased female body weight and increased resistance to acetone, suggesting components of the response unrelated to ethanol metabolism. The Adhs allele was favored by selection in all populations at the onset, but increased in frequency only in the selected lines of the southernmost population. There was a correlation between latitude and Adh frequency changes, suggesting that fitnesses of the Adh alleles were dependent on the genetic background. Genetic background also had a large effect on the loss of fitness due to selection. Genetic drift between replicate lines caused more variation in selection response than initial genetic differences between populations. This result demonstrates the importance of genetic drift in divergence among natural populations undergoing uniform selection, since the effective population sizes approached those of small natural populations. Drift caused greater divergence between selected replicates than control replicates. Implications of this result for the genetic model of selection response are discussed.     

Influence of defined gene blocks on the competitive ability of yeast genotypes

Competition experiments were carried out under varying exogenic and endogenic conditions. The genotypes were marked by combinations of two esterase loci, each with two alleles. When genotypes of the line W7 were used, there was no demonstrable influence of the gene blocks marked by the Est-1 locus on the competitive ability at temperatures of 21 and 29 C. However, genotypes carrying the fast allele of the Est-2 locus were favored. At 38 C, the outcome of the competition was reversed. The defined gene blocks showed different effects when interacting with different genetic backgrounds (line M7). Genotypes marked by the slow allele of the Est-2 locus were now favored (21 and 29 C), and even the gene blocks marked by the alleles of the Est-1 locus influenced the genotypes' competitive abilities. Again, the results were partly reversed at 38 C. The results are discussed with regard to the importance of enzyme variants for the genotypic selection value.     

An analysis of life history evolution in terms of the density-dependent Lefkovitch matrix model.

The evolution of demographic characteristics is considered in terms of the density-dependent Lefkovitch matrix model, which describes a species' population dynamics with a stage-specific pattern of reproduction and mortality. We obtain the invadability condition of a mutant-type into the wild-type population at the equilibrium state. The condition depends on the left and right eigenvectors at the equilibrium state. The condition depends on the left and right eigenvectors at the equilibrium state and the difference, between wild-type and mutant-type populations, of the values of elements in the Lefkovitch matrix at the equilibrium state. It is also shown that if elements of the density-dependent Lefkovitch matrix are decreasing functions of population density, then the equilibrium population density increases in the process of natural selection; that is, K-selection acts even on the stage-structured population. The evolution of life history in perennial plants is discussed through two models as an application of the above results. The evolution of perennial plants with no vegetative reproduction is analyzed in the first example. It is shown that whether monocarpic perennials (which reproduce once and die) or polycarpic perennial plants (which reproduce more than once) are favored depends on the cost of a produced seed. The second example concerns perennial plants that reproduce vegetatively. It is shown that whether monocarpic or polycarpic perennial plants are favored depends on the cost of a seed and that where vegetative reproduction is common, polycarpic perennials with no seed reproduction are favored.     

The Evolutionary Advantage of Recombination. II. Individual Selection for Recombination

Based on the Fisher-Muller theory of the evolution of recombination, an argument can be constructed predicting that a recessive allele favoring recombination will be favored, if there are either favorable or deleterious mutants occurring at other loci. In this case there is no clear distinction between individual and group selection. Computer simulation of populations segregating for recessive or dominant recombination alleles showed selection favoring recombination, except in the case of a dominant recombination allele with deleterious background mutants. The relationship of this work to parallel investigations by Williams and by Strobeck, Maynard Smith, and Charlesworth is explored. All seem to rely on the same phenomenon. There seems no reason to assume that the evolution of recombination must have occurred by group selection.     

家养动物选择信号研究进展

家养动物在人类生活中占有重要地位,它们都经历驯化而来,在自然和人工选择下,适应了当地环境和人类需要,形成了各类品种。驯化、自然和人工选择都会在基因组上留下选择信号。对这些选择信号研究,可以直接挖掘到功能基因,是目前最重要的功能基因筛选策略之一。当前已经对猪、鸡、牛、羊、狗等家养动物开展了选择信号研究,并挖掘了大量功能基因。本文主要概述了选择信号的种类和检测方法,简述其在家养动物中的研究,并讨论了选择信号分析的关键问题及其研究前景。     

Selection with Partial Selfing. I. Mass Selection

The expected responses to mass selection carried out before or after reproduction in a population whose members all have a fixed probability of self pollination (s) are formulated using covariances of relatives and their component quadratic functions for a model with arbitrary additive and dominance effects. The response measured in the first generation offspring after selection (immediate gain) can differ from that retained when the population has regained equilibrium (permanent gain). The population mean behaves in a predictable manner during the return to equilibrium, and its value at any time can be predicted from earlier generations. The permanent gain from selection after reproduction is always (1 + s)/2 times as large as that from selection before reproduction, but the relationship of the immediate gains depends on the genetic model assumed. Numerical analysis applied to a model with two alleles per locus and varying allele frequencies, dominance ratios and numbers of loci showed that the proportion of the immediate gain retained at equilibrium was reduced with the large inbreeding depression associated with increasing dominance levels and numbers of loci and was generally lower for selection after reproduction than before. In the absence of information as to the magnitude of genetic variances and inbreeding depression in species reproducing by partial selfing, the importance of this phenomenon is unknown.     

Regions of Stable Equilibria for Models of Differential Selection in the Two Sexes under Random Mating

The equilibrium structure of models of differential selection in the sexes is investigated. It is shown that opposing additive selection leads to stable polymorphic equilibria for only a restricted set of selection intensities, and that for weak selection the selection intensities must be of approximately the same magnitude in the sexes. General models of opposing directional selection, with arbitrary dominance, are investigated by considering simultaneously the stability properties of the trivial equilibria and the curve along which multiple roots appear. Numerical calculations lead us to infer that the average degree of dominance determines the equilibrium characteristics of models of opposing selection. It appears that if the favored alleles are, on the average, recessive, there may be multiple polymorphic equilibria, whereas only a single polymorphic equilibrium can occur when the favored alleles are, on the average, dominant. The principle that the average degree of dominance controls equilibrium behavior is then extended to models allowing directional selection in one sex with overdominance in the other sex, by showing that polymorphism is maintained if and only if the average fitness in heterozygotes exceeds one.     

Red Queen Processes Drive Positive Selection on Major Histocompatibility Complex (MHC) Genes

Major Histocompatibility Complex (MHC) genes code for proteins involved in the incitation of the adaptive immune response in vertebrates, which is achieved through binding oligopeptides (antigens) of pathogenic origin. Across vertebrate species, substitutions of amino acids at sites responsible for the specificity of antigen binding (ABS) are positively selected. This is attributed to pathogen-driven balancing selection, which is also thought to maintain the high polymorphism of MHC genes, and to cause the sharing of allelic lineages between species. However, the nature of this selection remains controversial. We used individual-based computer simulations to investigate the roles of two phenomena capable of maintaining MHC polymorphism: heterozygote advantage and host-pathogen arms race (Red Queen process). Our simulations revealed that levels of MHC polymorphism were high and driven mostly by the Red Queen process at a high pathogen mutation rate, but were low and driven mostly by heterozygote advantage when the pathogen mutation rate was low. We found that novel mutations at ABSs are strongly favored by the Red Queen process, but not by heterozygote advantage, regardless of the pathogen mutation rate. However, while the strong advantage of novel alleles increased the allele turnover rate, under a high pathogen mutation rate, allelic lineages persisted for a comparable length of time under Red Queen and under heterozygote advantage. Thus, when pathogens evolve quickly, the Red Queen is capable of explaining both positive selection and long coalescence times, but the tension between the novel allele advantage and persistence of alleles deserves further investigation.     

Drift or Selection: A Statistical Test of Gene Frequency Variation over Generations

The method used by Fisher and Ford (1947) to study the spread of a gene in a natural population has been modified to analyze the variation in allele frequencies from generation to generation in a common experimental procedure. A further analysis has been developed that is more sensitive to directional trends in the allele frequency over generations, and its use in detecting the action of directional selection on gene frequency at a locus is discussed. The power of each of these statistical tests is calculated for a number of cases, and the tests are applied to sets of isozyme data from Drosophila pseudoobscura and Zea mays.     

Conditional gene genealogies under strong purifying selection

The ancestral selection graph, conditioned on the allelic types in the sample, is used to obtain a limiting gene genealogical process under strong selection. In an equilibrium, two-allele system with strong selection, neutral gene genealogies are predicted for random samples and for samples containing at most one unfavorable allele. Samples containing more than one unfavorable allele have gene genealogies that differ greatly from neutral predictions. However, they are related to neutral gene genealogies via the well-known Ewens sampling formula. Simulations show rapid convergence to limiting analytical predictions as the strength of selection increases. These results extend the idea of a soft selective sweep to deleterious alleles and have implications for the interpretation of polymorphism among disease-causing alleles in humans.     

Genes Selection Comparative Study in Microarray Data Analysis

In response to the rapid development of DNA Microarray Technologies, many differentially expressed genes selection algorithms have been developed, and different comparison studies of these algorithms have been done. However, it is not clear how these methods compare with each other, especially when we used different developments tools. Here, we considered three commonly used differentially expressed genes selection approaches, namely: Fold Change, T-test and SAM, using Bioinformatics Matlab Toolbox and R/BioConductor. We used two datasets, issued from the affymetrix technology, to present results of used methods and software''s in gene selection process. The results, in terms of sensitivity and specificity, indicate that the behavior of SAM is better compared to Fold Change and T-test using R/BioConductor. While, no practical differences were observed between the three gene selection methods when using Bioinformatics Matlab Toolbox. In face of our result, the ROC curve shows that: on the one hand R/BioConductor using SAM is favored for microarray selection compared to the other methods. And, on the other hand, results of the three studied gene selection methods using Bioinformatics Matlab Toolbox are still comparable for the two datasets used.