Genomic Scans of Zygotic Disequilibrium and Epistatic SNPs in HapMap Phase III Populations

Previous theory indicates that zygotic linkage disequilibrium (LD) is more informative than gametic or composite digenic LD in revealing natural population history. Further, the difference between the composite digenic and maximum zygotic LDs can be used to detect epistatic selection for fitness. Here we corroborate the theory by investigating genome-wide zygotic LDs in HapMap phase III human populations. Results show that non-Africa populations have much more significant zygotic LDs than do Africa populations. Africa populations (ASW, LWK, MKK, and YRI) possess more significant zygotic LDs for the double-homozygotes (D_AABB) than any other significant zygotic LDs (D_AABb, D_AaBB, and D_AaBb), while non-Africa populations generally have more significant D_AaBb’s than any other significant zygotic LDs (D_AABB, D_AABb, and D_AaBB). Average r-squares for any significant zygotic LDs increase generally in an order of populations YRI, MKK, CEU, CHB, LWK, JPT, CHD, TSI, GIH, ASW, and MEX. Average r-squares are greater for D_AABB and D_AaBb than for D_AaBB and D_AABb in each population. YRI and MKK can be separated from LWK and ASW in terms of the pattern of average r-squares. All population divergences in zygotic LDs can be interpreted with the model of Out of Africa for modern human origins. We have also detected 19735-95921 SNP pairs exhibiting strong signals of epistatic selection in different populations. Gene-gene interactions for some epistatic SNP pairs are evident from empirical findings, but many more epistatic SNP pairs await evidence. Common epistatic SNP pairs rarely exist among all populations, but exist in distinct regions (Africa, Europe, and East Asia), which helps to understand geographical genomic medicine.     

The effect of positive assortative mating at one locus on a second linked locus

Summary Considerations proceed from a model of positive assortative mating based on genotype at one locus, with an arbitrary number of alleles, assuming no selection, mutation, or migration, hypothetically infinite population size, and discrete non-overlapping generations. From these conditions, inferences are made about the genotypic structure at a linked locus, as well as about the corresponding 2-locus gametic structure.The following main results are presented: in the course of the generations, the genotypic structure at the second locus and the 2-locus gametic structure always tend to a limit responsive to the initial conditions concerning the joint genotypic structure at the two loci and the degree of assortativity and linkage. A complete, analytical representation of the limits is given. In particular, if assortative mating is only partial and at the same time linkage is not complete, a population is not able to maintain a permanent deviation of the gametic structure from linkage equilibrium, and thus the genotypic structure at the second locus tends to Hardy-Weinberg proportions. On the other hand, if initial linkage disequilibrium is combined with partial assortative mating and complete linkage (or with complete assortative mating and unlinked loci) the population maintains this disequilibrium and thus the genotypic structure at the second locus need not tend to Hardy-Weinberg proportions. It turns out that the conditions not only of complete linkage, but also of unlinked loci together with complete assortativity, imply no change in gametic structure from the initial structure.In order to demonstrate the influence of several parameters on the speed of convergence to and the magnitude of the respective limits, several graphs are included.     

Recombination disequilibrium in ideal and natural populations

Following Hardy-Weinberg disequilibrium (HWD) occurring at a single locus and linkage disequilibrium (LD) between two loci in generations, we here proposed the third genetic disequilibrium in a population: recombination disequilibrium (RD). RD is a measurement of crossover interference among multiple loci in a random mating population. In natural populations besides recombination interference, RD may also be due to selection, mutation, gene conversion, drift and/or migration. Therefore, similarly to LD, RD will also reflect the history of natural selection and mutation. In breeding populations, RD purely results from recombination interference and hence can be used to build or evaluate and correct a linkage map. Practical examples from F₂, testcross and human populations indeed demonstrate that RD is useful for measuring recombination interference between two short intervals and evaluating linkage maps. As with LD, RD will be important for studying genetic mapping, association of haplotypes with disease, plant breading and population history.     

The population genetics of anisogamy

This paper analyses the population genetics of anisogamy controlled by a single locus, in both the haploid and diploid cases. The conclusions of Parker et al. (1972), based on computer calculations, are confirmed analytically. The effects of the existence of two mating types on the evolution of anisogamy are examined. Close linkage between a mating type locus and the gamete size locus may produce non-random associations of alleles, leading to disassortative fusion with respect to gamete size. With loose linkage, there is random association of alleles, but selection favours closer linkage.     

A Theoretical Model for Estimating Linkage in F2 Populations with Distorted Single Gene Segregation

A model is presented which allows estimation of linkage from dihybrid F₂ populations with distorted single gene segregation by applying the maximum-likelihood method. For different selection processes operating on one locus at either the gametic or the zygotic level, it can be demonstrated that, if the deficit is previously taken into account, testing for free recombination can be carried out without prior knowledge of the causes of this deficit. In the presence of linkage, the expected frequencies of two phenotypic classes depend on whether gametic or zygotic selection is operating. The remaining two classes can be utilized for the estimation of linkage as their frequency ratio is independent of these selection types. The application of this procedure to situations with coupling, incomplete penetrance, gametic and zygotic selection is discussed.     

Effects of epistasis on tests for linkage in self-pollinated species

Summary Tests for linkage based on covariances among relatives in self-pollinated species are usually based upon an assumption that epistasis is not important. This study was conducted to determine the impact of epistasis on, and to investigate the sensitivity of, such tests. Thirty covariances were calculated for each of ten non-epistatic and ten epistatic genetic models with varying probabilities of recombination between two coupling or repulsion loci. Each set of covariances was tested for linkage by comparing covariances calculated for the model with those expected for an additive-dominance model with no linkage. Results showed that the test for linkage is quite insensitive to the effects of linkage due to the disproportionate influence of inbreeding. Repulsion linkages should be easier to detect than coupling linkages for all models. Epistasis was found to mimic or counteract the effects of linkage. Tests for linkage based on covariances within a hierarchical mating design appear to be insensitive to linkage and may confuse the effects of linkage and epistasis.     

A slow selection analysis of two locus, two allele traits.

A deterministic, continuous time model for the dynamics of two locus, two allele Mendelian traits in a large randomly mating diploid population is derived. The model allows for frequency and time dependent birth and death rates. It is analyzed under the assumption that the selective forces acting in the population are small. Slow selection approximations to the system's solution are then constructed. Two particular cases are studied. First, when linkage between loci is tight, the population is shown to rapidly approach Hardy-Weinberg proportions, which then may vary on a (slow) time scale determined by differential fitness. In the case of constant birth and death rates, a measure of the population's fitness is shown to increase on the slow time scale after an initial rapid adjustment. The second case considered is for loose linkage; a population near linkage equilibrium is studied. It is shown that the epistatic parameters cancel and that the results agree with the tight linkage case to leading order. The linkage disequlibrium is described in both cases.     

Geographical Variation in a Quantitative Character

A model for the evolution of the local averages of a quantitative character under migration, selection, and random genetic drift in a subdivided population is formulated and investigated. Generations are discrete and nonoverlapping; the monoecious, diploid population mates at random in each deme. All three evolutionary forces are weak, but the migration pattern and the local population numbers are otherwise arbitrary. The character is determined by purely additive gene action and a stochastically independent environment; its distribution is Gaussian with a constant variance; and it is under Gaussian stabilizing selection with the same parameters in every deme. Linkage disequilibrium is neglected. Most of the results concern the covariances of the local averages. For a finite number of demes, explicit formulas are derived for (i) the asymptotic rate and pattern of convergence to equilibrium, (ii) the variance of a suitably weighted average of the local averages, and (iii) the equilibrium covariances when selection and random drift are much weaker than migration. Essentially complete analyses of equilibrium and convergence are presented for random outbreeding and site homing, the Levene and island models, the circular habitat and the unbounded linear stepping-stone model in the diffusion approximation, and the exact unbounded stepping-stone model in one and two dimensions.     

Is Fisher's model necessary for the theory of population improvement?

Summary It is shown here that genetic advance in one cycle of recurrent selection can be formulated directly in terms of covariances between relatives by application of the general statistical principle of linear prediction. For practical use of such formulae it is necessary to estimate the corresponding covariance between relatives from the mating design used. With General Combining Ability selection such estimation is direct. For other types of selection, it is necessary to derive associated covariances from other types of covariances but it is not necessary to use classical results of covariances between relatives in terms of genetic effects. Indeed, covariances can be derived without factorial decomposition of the genetic effects at one locus, i.e., without the concept of additivity and dominance. This approach allows a simple derivation of the genetic advance after n cycles of selection, followed by m generations of intercrossing, with a minimum of assumptions.     

The strength of the selection barrier between populations

Genetic differences among populations exposed to selection form barriers against genetic exchange by mortality among hybrids. The strength of such a selection barrier, with which one (recipient) population reacts against immigration from another (donor) population, may be measured as the cumulative mean fitness of hybrids and their descendants relative to the fitness of the recipient population. Previous work analysed a case of weak selection with pairwise epistatic interactions by assuming small genetic distance between two populations in contact. The present study allows large genetic difference between the donor and recipient populations and considers weak multilocus selection with arbitrary epistatic interactions between two or more linked loci. An approximate analytical expression for the barrier strength is obtained as an expansion in which the strength of selection plays the role of a small parameter. It is shown that allele frequencies and gametic linkage disequilibria contribute in different ways to the strength of the selection barrier.     

Effect of gametic disequilibrium on selection in an autotetraploid population

Summary The effects of a gametic disequilibrium (DSE) in an autotetraploid population on response to selection as measured by the covariance of selection were investigated. The theoretical responses were calculated for mass selection [Mass (1)] and half-sib progeny test selection (HSPT) in a two-allele (B and b), single locus, autotetraploid population. The complexity of calculations precluded analytical expressions for the covariances so numerical analysis was used assuming the following genetic models: monoplex dominance, partial monoplex dominance, duplex dominance, partial duplex dominance, and additive gene action.The results indicated the DSE could greatly affect the covariance of selection. For a constant allele frequency the DSE might double the covariance expected with selection in a population at random mating equilibrium (RME) of gametes, but in other instances approach zero. For all genetic models and the two breeding methods investigated the covariance of selection was always increased when the frequency of BB gamete exceeded p² (where p is frequency of allele B) and decreased when the frequency of BB gamete was less than p². The possible incorporation of this information into a long term breeding program and some other ramifications were briefly discussed.With the DSE the covariances of selection with HSPT and Mass (1) had a proportionality of 1:2, respectively, with the additive genetic model, but this relationship rarely occurred for other genetic models. The deviations from this ratio were not large in comparison to differences between selection in populations in DSE and RME.Cooperative investigations of the Alfalfa Production Research Unit, United State Department of Agriculture, Agricultural Research Service, and the Nevada Agricultural Experiment Station, Reno, Nevada. Paper No. 512. Scientific Journal Series, Nevada Agricultural Experiment Station     

Mating system and the evolution of quantitative traits: an experimental study of Mimulus guttatus

Abstract The mating system of a population profoundly influences its evolution. Inbreeding alters the balance of evolutionary forces that determine the amount of genetic variation within a population. It redistributes that variation among individuals, altering heritabilities and genetic correlations. Inbreeding even changes the basic relationships between these genetic statistics and response to selection. If populations differing only in mating system are exposed to the same selection pressures, will they respond in qualitatively different ways? Here, we address this question by imposing selection on an index of two negatively correlated traits (flower size and development rate) within experimental populations that reproduce entirely by outcrossing, entirely by self‐fertilizing, or by a mixture of outcrossing and selfing. Entirely selfing populations responded mainly by evolving larger flowers whereas outcrossing populations also evolved more rapid development. Divergence occurred despite an equivalent selection regime and no direct effect of mating system on fitness. The study provides an experimental demonstration of how the interaction of selection, genetic drift, and mating system can produce dramatic short‐term changes in trait means, variances, and covariances.     

Adaptive dynamics in diploid, sexual populations and the evolution of reproductive isolation

Evolutionary branching is the process whereby an initially monomorphic population evolves to a point where it undergoes disruptive selection and splits up into two phenotypically diverging lineages. We studied evolutionary branching in three models that are ecologically identical but that have different genetic systems. The first model is clonal, the second is sexual diploid with additive genetics on a single locus and the third is like the second but with an additional locus for mate choice. Evolutionary branching occurred under exactly the same ecological circumstances in all three models. After branching the evolutionary dynamics may be qualitatively different. In particular, in the diploid, sexual models there can be multiple evolutionary outcomes whereas in the corresponding clonal model there is only one. We showed that evolutionary branching favours the evolution of (partial) assortative mating and that this in turn effectively restores the results from the clonal model by rendering the alternative outcomes unreachable except for the one that also occurs in the clonal model. The evolution of assortative mating during evolutionary branching can be interpreted as the initial phase of sympatric speciation with phenotypic divergence and partial reproductive isolation.     

The effect of positive assortative mating at one locus on a second linked locus

Summary A model for positive assortative mating based on genotype for one locus is employed to investigate the effect of this mating system on the genotypic structure of a second linked locus as well as on the joint genotypic structure of these two loci. It is shown that the second locus does not attain a precise positive assortative mating structure, but yet it shares a property that is characteristic of positive assortative mating, namely an increase in the frequency of homozygotes over that typically found in panmictic structures. Given any arbitrary genotypic structure for the parental population, the resulting offspring generation possesses a structure at the second locus that does not depend on the recombination frequency, while the joint structure of course does. In case assortative mating as well as linkage are not complete, there exists a unique joint equilibrium state for the two loci, which is characterized by complete stochastic independence between the two loci as well as by Hardy-Weinberg proportions at the second locus. For the second locus alone, Hardy-Weinberg equilibrium is realized if and only if gametic linkage equilibrium and an additionally specified condition are realized.     

Effects of genetic drift on variance components under a general model of epistasis

We analyze the changes in the mean and variance components of a quantitative trait caused by changes in allele frequencies, concentrating on the effects of genetic drift. We use a general representation of epistasis and dominance that allows an arbitrary relation between genotype and phenotype for any number of diallelic loci. We assume initial and final Hardy-Weinberg and linkage equilibrium in our analyses of drift-induced changes. Random drift generates transient linkage disequilibria that cause correlations between allele frequency fluctuations at different loci. However, we show that these have negligible effects, at least for interactions among small numbers of loci. Our analyses are based on diffusion approximations that summarize the effects of drift in terms of F, the inbreeding coefficient, interpreted as the expected proportional decrease in heterozygosity at each locus. For haploids, the variance of the trait mean after a population bottleneck is var(delta(z)) = sigma(n)k=1 FkV(A(k)), where n is the number of loci contributing to the trait variance, V(A(1)) = V(A) is the additive genetic variance, and V(A(k)) is the kth-order additive epistatic variance. The expected additive genetic variance after the bottleneck, denoted (V*(A)), is closely related to var(delta(z)); (V*(A)) = (1 - F) sigma(n)k=1 kFk-1V(A(k)). Thus, epistasis inflates the expected additive variance above V(A)(1 - F), the expectation under additivity. For haploids (and diploids without dominance), the expected value of every variance component is inflated by the existence of higher order interactions (e.g., third-order epistasis inflates (V*(AA. This is not true in general with diploidy, because dominance alone can reduce (V*(A)) below V(A)(1 - F) (e.g., when dominant alleles are rare). Without dominance, diploidy produces simple expressions: var(delta(z)) = sigma(n)k=1 (2F)kV(A(k)) and (V(A)) = (1 - F) sigma(n)k=1 k(2F)k-1V(A(k)). With dominance (and even without epistasis), var(delta(z)) and (V*(A)) no longer depend solely on the variance components in the base population. For small F, the expected additive variance simplifies to (V*(A)) approximately equal to (1 - F)V(A) + 4FV(AA) + 2FV(D) + 2FC(AD), where C(AD) is a sum of two terms describing covariances between additive effects and dominance and additive X dominance interactions. Whether population bottlenecks lead to expected increases in additive variance depends primarily on the ratio of nonadditive to additive genetic variance in the base population, but dominance precludes simple predictions based solely on variance components. We illustrate these results using a model in which genotypic values are drawn at random, allowing extreme and erratic epistatic interactions. Although our analyses clarify the conditions under which drift is expected to increase V(A), we question the evolutionary importance of such increases.     

Genetic evaluation by best linear unbiased prediction using marker and trait information in a multibreed population.

Genetic evaluation by best linear unbiased prediction (BLUP) requires modeling genetic means, variances, and covariances. This paper presents theory to model means, variances, and covariances in a multibreed population, given marker and breed information, in the presence of gametic disequilibrium between the marker locus (ML) and linked quantitative trait locus (MQTL). Theory and algorithms are presented to construct the matrix of conditional covariances between relatives (Gv) for the MQTL effects in a multibreed population and to obtain the inverse of Gv efficiently. Theory presented here accounts for heterogeneity of variances among pure breeds and for segregation variances between pure breeds. A numerical example was used to illustrate how the theory and algorithms can be used for genetic evaluation by BLUP using marker and trait information in a multibreed population.     

The Genetic Covariance between Characters Maintained by Pleiotropic Mutations

A statistical genetic model of a multivariate phenotype is derived to investigate the covariation of pleiotropic mutations with additive effects under the combined action of phenotypic selection, linkage and the mating system. Equilibrium formulas for large, randomly mating populations demonstrate that, when selection on polygenic variation is much smaller than twice the harmonic mean recombination rate between loci with interacting fitnesses, linkage disequilibrium is negligible and pleiotropy is the main cause of genetic correlations between characters. Under these conditions, approximate expressions for the dynamics of the genetic covariances due to pleiotropic mutations are obtained. Patterns of genetic covariance between characters and their evolution are discussed with reference to data on polygenic mutation, chromosomal organization and morphological integration.     

Genetic structure and linkage disequilibrium in landrace populations of barley in Sardinia

Multilocus digenic linkage disequilibria (LD) and their population structure were investigated in eleven landrace populations of barley (Hordeum vulgare ssp. vulgare L.) in Sardinia, using 134 dominant simple-sequence amplified polymorphism markers. The analysis of molecular variance for these markers indicated that the populations were partially differentiated (F(ST) = 0.18), and clustered into three geographic areas. Consistent with this population pattern, STRUCTURE analysis allocated individuals from a bulk of all populations into four genetic groups, and these groups also showed geographic patterns. In agreement with other molecular studies in barley, the general level of LD was low (13% of locus pairs, with P < 0.01) in the bulk of 337 lines, and decayed steeply with map distance between markers. The partitioning of multilocus associations into various components indicated that genetic drift and founder effects played a major role in determining the overall genetic makeup of the diversity in these landrace populations, but that epistatic homogenising or diversifying selection was also present. Notably, the variance of the disequilibrium component was relatively high, which implies caution in the pooling of barley lines for association studies. Finally, we compared the analyses of multilocus structure in barley landrace populations with parallel analyses in both composite crosses of barley on the one hand and in natural populations of wild barley on the other. Neither of these serves as suitable mimics of landraces in barley, which require their own study. Overall, the results suggest that these populations can be exploited for LD mapping if population structure is controlled.     

ON EVOLUTION UNDER SEXUAL AND VIABILITY SELECTION: A TWO-LOCUS DIPLOID MODEL

A two-locus diploid model of sexual selection is presented in which the two loci govern, respectively, a trait limited in expression in one sex (generally male) and the mating preferences of the other sex (generally female). The viability of a male depends on its genotype at the trait locus. In contrast, all females are equally viable and all individuals are equally fertile with respect to the two loci. Near fixation at both loci, evolution at the mating locus is neutral and hence a new mating preference allele will increase only through random genetic drift or through a correlated response to the increase of a new advantageous trait allele. If, however, a polymorphism is already maintained at the trait locus through overdominance in fitness then the increase of a rare preference allele depends only on the recombination rate between the loci and not on the new preference scheme.     

Population dynamics of a gametophytic factor controlling selective fertilization

The population behavior of a gametophytic factor (Ga) which involves gametic selection due to failure ofga pollen onGa Ga orGa ga styles in competition withGa pollen, was investigated by computer simulation. A constant versus randomly varying gametic selection parameter (k) and four different schemes of zygotic selection were introduced in this model for analyzing conditions favorable for the maintenance of locusGa polymorphic in a large, mixed selfing and random mating population. Stable polymorphism was obtained only with rather substantial heterozygote advantage at locusGa whereas the opposing pressures of gametic and zygotic selection yielded fixation of alleleGa orga around a critical value of k instead of a range of k-values allowing nontrivial equilibria. With weak selection and stochastic k, however, very slow rates of change in the genotypic proportions allowed transient polymorphism. In these cases, the rate of outcrossing (t) and the initial frequency ofGa were critical in determining the rate of allelic substitution. Moreover, low values of t allowed the replacement of alleleGa byga even with rather weak zygotic selection. These findings on the balance between gametic and zygotic selection and a markedly frequency-dependent process are briefly discussed in relation to the dynamics of similar factors involving the mating system.