Influence of Gene Interaction on Complex Trait Variation with Multilocus Models

Although research effort is being expended into determining the importance of epistasis and epistatic variance for complex traits, there is considerable controversy about their importance. Here we undertake an analysis for quantitative traits utilizing a range of multilocus quantitative genetic models and gene frequency distributions, focusing on the potential magnitude of the epistatic variance. All the epistatic terms involving a particular locus appear in its average effect, with the number of two-locus interaction terms increasing in proportion to the square of the number of loci and that of third order as the cube and so on. Hence multilocus epistasis makes substantial contributions to the additive variance and does not, per se, lead to large increases in the nonadditive part of the genotypic variance. Even though this proportion can be high where epistasis is antagonistic to direct effects, it reduces with multiple loci. As the magnitude of the epistatic variance depends critically on the heterozygosity, for models where frequencies are widely dispersed, such as for selectively neutral mutations, contributions of epistatic variance are always small. Epistasis may be important in understanding the genetic architecture, for example, of function or human disease, but that does not imply that loci exhibiting it will contribute much genetic variance. Overall we conclude that theoretical predictions and experimental observations of low amounts of epistatic variance in outbred populations are concordant. It is not a likely source of missing heritability, for example, or major influence on predictions of rates of evolution.     

Quantitative genetic models for the balance between migration and stabilizing selection

The evolution of a quantitative trait subject to stabilizing selection and immigration, with the immigrants deviating from the local optimum, is considered under a number of different models of the underlying genetic basis of the trait. By comparing exact predictions under the infinitesimal model obtained using numerical methods with predictions of a simplified approximate model based on ignoring linkage disequilibrium, the increase in the expressed genetic variance as a result of linkage disequilibrium generated by migration is shown to be relatively small and negligible, provided that the genetic variance relative to the squared deviation of immigrants from the local optimum is sufficiently large or selection and migration is sufficiently weak. Deviation from normality is shown to be less important by comparing predictions of the infinitesimal model with a model presupposing normality. For a more realistic symmetric model, involving a finite number of loci only, no linkage and equal effects and frequencies across loci, additional changes in the genetic variance arise as a result of changes in underlying allele frequencies. Again, provided that the genetic variance relative to the squared deviation of the immigrants from the local optimum is small, the difference between the predictions of infinitesimal and the symmetric model are small unless the number of loci is very small. However, if the genetic variance relative to the squared deviation of the immigrants from the local optimum is large, or if selection and migration are strong, both linkage disequilibrium and changes in the genetic variance as a result of changes in underlying allele frequencies become important.     

WILL POPULATION BOTTLENECKS AND MULTILOCUS EPISTASIS INCREASE ADDITIVE GENETIC VARIANCE?

We apply new analytical methods to understand the consequences of population bottlenecks for expected additive genetic variance. We analyze essentially all models for multilocus epistasis that have been numerically simulated to demonstrate increased additive variance. We conclude that for biologically plausible models, large increases in expected additive variance--attributable to epistasis rather than dominance--are unlikely. Naciri-Graven and Goudet (2003) found that as the number of epistatically interacting loci increases, additive variance tends to be inflated more after a bottleneck. We argue that this result reflects biologically unrealistic aspects of their models. Specifically, as the number of loci increases, higher-order epistatic interactions become increasingly important in these models, with an increasing fraction of the genetic variance becoming nonadditive, contrary to empirical observations. As shown by Barton and Turelli (2004), without dominance, conversion of nonadditive to additive variance depends only on the variance components and not on the number of loci per se. Numerical results indicating that more inbreeding is needed to produce maximal release of additive variance with more loci follow directly from our analytical results, which show that high levels of inbreeding (F > 0.5) are needed for significant conversion of higher-order components. We discuss alternative approaches to modeling multilocus epistasis and understanding its consequences.     

Evaluating the utility of microsatellites for investigations of autopolyploid taxa

Autopolyploid taxa present numerous challenges for population genetic analyses due to difficulties determining allele dosage. Dosage ambiguity hinders accurate assessment of allele frequencies, multilocus genotypes (MLGTs), as well as levels and patterns of clonality. The pervasiveness of polyploidy in the evolutionary history of plant taxa makes this a recurring problem. Whereas diploidization of loci may occur over time, duplication of at least some loci is still frequently evident. Fortunately, with high-quality allozyme gels, it is possible to accurately infer allele dosage and, thus, determine exact MLGTs. However, accurately assessing dosage of microsatellite peaks is nearly impossible when studying wild populations with a large number of alleles per locus. Even if precise knowledge of genotypes is not required, for comparable numbers of alleles per locus and loci, the number of "phenotypes" is always lower with microsatellites than allozymes due to the inability to assess allele dosage. Microsatellite loci typically have more alleles per locus relative to allozymes although fewer loci are generally employed. Here, we present a mathematical model for comparing the relative utility of simple sequence repeat (SSR) versus allozyme markers to discriminate MLGTs. For example, the average plant allozyme study (2.6 alleles per locus, 10 polymorphic loci) has better discriminating power than SSR markers with 10 alleles at each of 3 loci, 9 alleles at 4 loci, 6 alleles at 5 loci, 5 alleles at 6 loci, and 4 alleles at 8 loci, demonstrating the value of assessing the relative discriminating power of these markers.     

Estimation of additive, dominance and epistatic variance components using finite locus models implemented with a single-site Gibbs and a descent graph sampler

In a previous contribution, we implemented a finite locus model (FLM) for estimating additive and dominance genetic variances via a Bayesian method and a single-site Gibbs sampler. We observed a dependency of dominance variance estimates on locus number in the analysis FLM. Here, we extended the FLM to include two-locus epistasis, and implemented the analysis with two genotype samplers (Gibbs and descent graph) and three different priors for genetic effects (uniform and variable across loci, uniform and constant across loci, and normal). Phenotypic data were simulated for two pedigrees with 6300 and 12,300 individuals in closed populations, using several different, non-additive genetic models. Replications of these data were analysed with FLMs differing in the number of loci. Simulation results indicate that the dependency of non-additive genetic variance estimates on locus number persisted in all implementation strategies we investigated. However, this dependency was considerably diminished with normal priors for genetic effects as compared with uniform priors (constant or variable across loci). Descent graph sampling of genotypes modestly improved variance components estimation compared with Gibbs sampling. Moreover, a larger pedigree produced considerably better variance components estimation, suggesting this dependency might originate from data insufficiency. As the FLM represents an appealing alternative to the infinitesimal model for genetic parameter estimation and for inclusion of polygenic background variation in QTL mapping analyses, further improvements are warranted and might be achieved via improvement of the sampler or treatment of the number of loci as an unknown.     

Detection of putative quantitative trait loci in line crosses under infinitesimal genetic models

Quantitative trait locus (QTL) mapping studies often employ segregating generations derived from a cross between genetically divergent inbred lines. In the analysis of such data it is customary to fit a single QTL and use a null hypothesis which assumes that the genomic region under study contributes no genetic variance. To explore the situation in which multiple linked genes contribute to the genetic variance, we simulated an F₂-mapping experiment in which the genetic difference between the two original inbred strains was caused by a large number of loci, each having equal effect on the quantitative trait. QTLs were either in coupling, dispersion or repulsion phase in the base population of inbred lines, with the expected F₂ genetic variance explained by the QTLs being equivalent in the three models. Where QTLs were in coupling phase, one inbred line was fixed for all plus alleles, and the other line was fixed for minus alleles. Where QTLs were in dispersion phase, they were assumed to be randomly fixed for one or other allele (as if the inbred lines had evolved from a common ancestor by random drift). Where QTLs were in repulsion phase alleles within an inbred line were alternating plus and minus at adjacent loci, and alternative alleles were fixed in the two inbred lines. In all these genetic models a standard interval mapping test statistic used to determine whether there is a QTL of large effect segregating in the population was inflated on average. Furthermore, the use of a threshold for QTL detection derived under the assumption that no QTLs were segregating would often lead to spurious conclusions regards the presence of genes of large effects (i.e. type I errors). The employment of an alternative model for the analysis, including linked markers as cofactors in the analysis of a single interval, reduced the problem of type I error rate, although test statistics were still inflated relative to the case of no QTLs. It is argued that in practice one should take into account the difference between the strains or the genetic variance in the F₂ population when setting significance thresholds. In addition, tests designed to probe the adequacy of a single-QTL model or of an alternative infinitesimal coupling model are described. Such tests should be applied in QTL mapping studies to help dissect the true nature of genetic variation.     

A multilocus-multiallele analysis of frequency-dependent selection induced by intraspecific competition

The study of the mechanisms that maintain genetic variation has a long history in population genetics. We analyze a multilocus-multiallele model of frequency- and density-dependent selection in a large randomly mating population. The number of loci and the number of alleles per locus are arbitrary. The n loci are assumed to contribute additively to a quantitative character under stabilizing or directional selection as well as under frequency-dependent selection caused by intraspecific competition. We assume the strength of stabilizing selection to be weak, whereas the strength of frequency dependence may be arbitrary. Density-dependence is induced by population regulation. Our main result is a characterization of the equilibrium structure and its stability properties in terms of all parameters. It turns out that no equilibrium exists with more than two alleles segregating per locus. We give necessary and sufficient conditions on the strength of frequency dependence to ensure the maintenance of multilocus polymorphism. We also give explicit formulas on the number of polymorphic loci maintained at equilibrium. These results are based on the assumption that selection is sufficiently weak compared with recombination, so that linkage equilibrium can be assumed. If additionally the population size is assumed to be constant, we prove that the dynamics of the model form a generalized gradient system. For the model in its general form we are able to derive necessary and sufficient conditions for the stability of the monomorphic equilibria. Furthermore, we briefly analyze a special symmetric two-locus two-allele model for a constant population size but allowing for linkage disequilibrium. Finally, we analyze a single diallelic locus with dominance to illustrate the complications that can occur if the assumption of additivity is relaxed.     

Genotype‐by‐environment interaction and the Dobzhansky–Muller model of postzygotic isolation

The Dobzhansky–Muller (D–M) model of reproductive isolation (RI) posits that hybrid sterility and inviability result from negative epistatic interactions between alleles at a minimum of two genes. This standard model makes several implicit assumptions, including a lack of environmental effects and genotype‐by‐environment interactions (GEI) involving hybrid sterility and hybrid inviability loci. Here we relax this assumption of the standard D–M model. By doing so, several patterns of the genetic architecture of RI change. First, a novel single‐locus model of postzygotic RI emerges. Several indirect lines of evidence are discussed in support of the model, but we conclude that this new single‐locus model is currently no more supported than previous ones. Second, when multilocus D–M models incorporating GEI are considered, we find that the number of potential negative epistatic interactions increases dramatically over the number predicted by the standard D–M model, even when only the most simple case of two‐allele interactions are considered. Third, these multilocus models suggest that some previous generalizations about the evolutionary genetics of postzygotic RI may not necessarily hold. Our findings also suggest that the evolution of postzygotic RI may be more likely when the expression of traits driving speciation is affected by the environment, since there appears to be a greater spectrum of potential hybrid incompatibilities under the D–M model incorporating GEI.     

Clines in polygenic traits

This article outlines theoretical models of clines in additive polygenic traits, which are maintained by stabilizing selection towards a spatially varying optimum. Clines in the trait mean can be accurately predicted, given knowledge of the genetic variance. However, predicting the variance is difficult, because it depends on genetic details. Changes in genetic variance arise from changes in allele frequency, and in linkage disequilibria. Allele frequency changes dominate when selection is weak relative to recombination, and when there are a moderate number of loci. With a continuum of alleles, gene flow inflates the genetic variance in the same way as a source of mutations of small effect. The variance can be approximated by assuming a Gaussian distribution of allelic effects; with a sufficiently steep cline, this is accurate even when mutation and selection alone are better described by the 'House of Cards' approximation. With just two alleles at each locus, the phenotype changes in a similar way: the mean remains close to the optimum, while the variance changes more slowly, and over a wider region. However, there may be substantial cryptic divergence at the underlying loci. With strong selection and many loci, linkage disequilibria are the main cause of changes in genetic variance. Even for strong selection, the infinitesimal model can be closely approximated by assuming a Gaussian distribution of breeding values. Linkage disequilibria can generate a substantial increase in genetic variance, which is concentrated at sharp gradients in trait means.     

Consequences of genetic linkage for the maintenance of sexually antagonistic polymorphism in hermaphrodites

When selection differs between males and females, pleiotropic effects among genes expressed by both sexes can result in sexually antagonistic selection (SA), where beneficial alleles for one sex are deleterious for the other. For hermaphrodites, alleles with opposing fitness effects through each sex function represent analogous genetic constraints on fitness. Recent theory based on single‐locus models predicts that the maintenance of SA genetic variation should be greatly reduced in partially selfing populations. However, selfing also reduces the effective rate of recombination, which should facilitate selection on linked allelic combinations and expand opportunities for balancing selection in a multilocus context. Here, I develop a two‐locus model of SA selection for simultaneous hermaphrodites, and explore the joint influence of linkage, self‐fertilization, and dominance on the maintainance of SA polymorphism. I find that the effective reduction in recombination caused by selfing significantly expands the parameter space where SA polymorphism can be maintained relative to single‐locus models. In particular, linkage facilitates the invasion of male‐beneficial alleles, partially compensating for the “female‐bias” in the net direction of selection created by selfing. I discuss the implications of accounting for linkage among SA loci for the maintenance of SA genetic variation and mixed mating systems in hermaphrodites.     

Estimation of effective population size and migration rate from one- and two-locus identity measures

Standard methods for inferring demographic parameters from genetic data are based mainly on one-locus theory. However, the association of genes at different loci (e.g., two-locus identity disequilibrium) may also contain some information about demographic parameters of populations. In this article, we define one- and two-locus parameters of population structure as functions of one- and two-locus probabilities for the identity in state of genes. Since these parameters are known functions of demographic parameters in an infinite island model, we develop moment-based estimators of effective population size and immigration rate from one- and two-locus parameters. We evaluate this method through simulation. Although variance and bias may be quite large, increasing the number of loci on which the estimates are derived improves the method. We simulate an infinite allele model and a K allele model of mutation. Bias and variance are smaller with increasing numbers of alleles per locus. This is, to our knowledge, the first attempt of a joint estimation of local effective population size and immigration rate.     

Dependence of expected heterozygosity on locus number with stabilizing selection and drift

I determine expected levels of heterozygosity in two allele multilocus models with mutation, stabilizing selection and drift. In the range 2 to 32 loci, the per locus heterozygosity can depend on the locus number. The per locus heterozygosity for ten loci can be as low as three fourths of the per locus heterozygosity in the limit, as the number of loci gets large. Simulations indicate that this dependence on locus number is not due to the population approaching equilibria at which the mean differs from the optimum, but is due to changes in the substitution rate as a function of the number of loci.     

Cost‐effective development of highly polymorphic microsatellite in Japanese quail facilitated by next‐generation sequencing

Next‐generation sequencing technologies permit rapid and cost‐effective identification of numerous putative microsatellite loci. Here, from the genome sequences of Japanese quail, we developed microsatellite markers containing dinucleotide repeats and employed these for characterisation of genetic diversity and population structure. A total of 385 individuals from 12 experimental and one wild‐derived Japanese quail lines were genotyped with newly developed autosomal markers. The maximum number of alleles, expected heterozygosity and polymorphic information content (PIC) per locus were 10, 0.80 and 0.77 respectively. Approximately half of the markers were highly informative (PIC ≥ 0.50). The mean number of alleles per locus and observed heterozygosity within a line were in the range of 1.3–4.1 and 0.11–0.53 respectively. Compared with the wild‐derived line, genetic diversity levels were low in the experimental lines. Genetic differentiation (F_ST) between all pairs of the lines ranged from 0.13 to 0.83. Genetic clustering analyses based on multilocus genotypes of individuals showed that most individuals formed clearly defined clusters corresponding to the origins of the lines. These results suggest that Japanese quail experimental lines are highly structured. Microsatellite markers developed in this study may be effective for future genetic studies of Japanese quail.     

Microsatellite assessment of multiple paternity in natural populations of a live-bearing fish,Gambusia holbrooki

Three polymorphic microsatellite loci were isolated and employed to examine paternity patterns in two natural populations of live-bearing mosquitofish, Gambusia holbrooki. Each locus displayed four to five alleles per population in samples of presumably unrelated adults. Nearly 900 embryos from a total of 50 pregnant females were assayed individually, and paternal alleles in each embryo were identified. Counts of paternal alleles, Mendelian segregation patterns, multilocus allelic associations and genetic relatedness coefficients were employed to estimate the minimum and effective numbers of fathers per brood. At least 90% of the assayed broods were shown to have been fathered by multiple males, a figure substantially higher than previous estimates based on less polymorphic genetic loci. However, the genetic data yield a face-value estimate of only about 2.2 fathers per brood, a number that seems perhaps surprisingly low based on frequencies of attempted copulations by males. Both biological and sampling factors that might bias mean sire counts downward are considered. Although higher sire counts per brood might be obtained from loci with even greater numbers of alleles, little statistical room remains for higher frequency estimates of multiple paternity in Gambusia.     

Multilocus Structure of Natural Populations of HORDEUM SPONTANEUM

The association of alleles among different loci was studied in natural populations of Hordeum spontaneum, the evolutionary progenitor of cultivated barley. The variance of the number of heterozygous loci in two randomly chosen gametes affords a useful measure of such association. The behavior of this statistic in several particular models is described. Generally, linkage (gametic phase) disequilibrium tends to increase the variance above the value expected under complete independence. This increase is greatest when disequilibria are such as to maximize the sum of squares of the two-locus gametic frequencies.—When data on several loci per individual are available, the observed variance may be tested for its agreement with that expected under the hypothesis of complete interlocus independence, using the sampling theory of this model. When applied to allozyme data from 26 polymorphic populations of wild barley, this test demonstrated the presence of geographically widespread multilocus organization. On average, the variance was 80% higher than expected under random association. Gametic frequencies for four esterase loci in both of these populations of wild barley and two composite crosses of cultivated barley were analyzed. Most generations of the composites showed less multilocus structure, as measured by the indices of association, than the wild populations.     

Multilocus selection in subdivided populations II. Maintenance of polymorphism under weak or strong migration

The potential of maintaining multilocus polymorphism by migration-selection balance is studied. A large population of diploid individuals is distributed over finitely many demes connected by migration. Generations are discrete and nonoverlapping, selection may vary across demes, and loci are multiallelic. It is shown that if migration and recombination are strong relative to selection, then with weak or no epistasis and intermediate dominance at every locus and in every deme, arbitrarily many alleles can be maintained at arbitrarily many loci at a stable equilibrium. If migration is weak relative to selection and recombination, then with weak or no epistasis and intermediate dominance at every locus and in every deme, as many alleles as there are demes can be maintained at arbitrarily many loci at equilibrium. In both cases open sets of such parameter combinations are constructed, thus the results are robust with respect to small, but arbitrary, perturbations in the parameters. For weak migration, the number of demes is, in fact, a generic upper bound to the number of alleles that can be maintained at any locus. Thus, several scenarios are identified under which multilocus polymorphism can be maintained by migration-selection balance when this is impossible in a panmictic population.      

凡纳滨对虾繁殖中不同亲本对子代遗传贡献率的差异

利用5个含有稀有等位基因的高度多态性微卫星位点比较了凡纳滨对虾繁殖中不同亲本对子代遗传贡献率的差异。通过稀有等位基因的5个微卫星位点能够对亲代和子代的谱系进行明确的鉴别。10个亲代个体中有8个个体对子代群体的基因库有贡献,不同个体之间的贡献率存在差别,最高为54.28%,最低为8.57%。在亲代和子代群体遗传结构的分析中,子代等位基因的数目与亲代相比降低了11.11%。子代的平均期望杂合度(He)、平均观测杂合度(Ho)和平均多态性信息含量(PIC)等指标均低于亲代。实验结果表明:亲本对子代基因库的贡献率的差异也是造成子代群体遗传变异水平降低的原因之一;微卫星标记可作为一种有效的工具用于对虾系谱的确认、人工繁育群体遗传多样性水平的监测等方面     

Genetic improvement of production while maintaining fitness

Selection for production tends to decrease fitness, in particular, major components such as reproductive performance. Under an infinitesimal genetic model restricted index selection can maintain reproductive performance while improving production. However, reproductive traits are thought to be controlled by a finite number of recessive alleles at low frequency. Culling for low reproduction may weed out the negative homozygous genotypes for reproduction in any generation, thus controlling the frequencies of alleles negative for reproduction. Restricted index selection, culling for low reproduction and a new method called empirical restricted index selection were compared for their efficiency in improving production while maintaining reproduction. Empirical restricted index selection selects animals that have on average the highest estimated breeding values for production and on average the same estimated breeding values for reproduction as the base population. An infinitesimal genetic model and models with a finite number of loci for reproduction with rare deleterious recessive alleles, which have additive, dominant or no pleiotropic effects on production, were considered. When reproduction was controlled by a finite number of loci with rare recessive alleles, restricted index selection could not maintain reproduction. The culling of 20% of the animals on reproduction maintained reproduction with all genetic models, except for the model where loci for reproduction had additive effects on production. Empirical restricted selection maintained reproduction with all models and yielded higher production responses than culling on reproduction, except when there were dominant pleiotropic effects on production.     

Least squares interval mapping of quantitative trait loci under the infinitesimal genetic model in outbred populations.

Genetic marker and phenotypic data for a quantitative trait were simulated on 20 paternal half-sib families with 100 progeny to investigate properties of within-family-regression interval mapping of a postulated single quantitative trait locus (QTL) in a marker interval under the infinitesimal genetic model, which has been the basis of the application of quantitative genetics to genetic improvement programs, and to investigate use of the infinitesimal model as null hypothesis in testing for presence of a major QTL. Genetic effects on the marked chromosome were generated based on a major gene model, which simulated a central biallelic QTL, or based on 101 biallelic QTL of equal effect, which approximated the infinitesimal model. The marked chromosome contained 0, 3.3%, 13.3%, or 33.3% of genetic variance and heritability was 0.25 or 0.70. Under the polygenic model with 3.3% of genetic variance on the marked chromosome, which corresponds to the infinitesimal model for the bovine, significant QTL effects were found for individual families. Correlations between estimates of QTL effects and true chromosome substitution effects were 0.29 and 0.47 for heritabilities of 0.25 and 0.70 but up to 0.85 with 33.3% of polygenic variance on the marked chromosome. These results illustrate the potential of marker-assisted selection even under the infinitesimal genetic model. Power of tests for presence of QTL was substantially reduced when the polygenic model with 3.3% of genetic variance on the chromosome was used as a null hypothesis. The ability to determine whether genetic variance on a chromosome was contributed by a single QTL of major effect or a large number of QTL with minor effects, corresponding to the infinitesimal model, was limited.     

Change in Genetic Variance under Selection in a Self-Fertilizing Population

In this study we show how the genetic variance of a quantitative trait changes in a self-fertilizing population under repeated cycles of truncation selection, with the analysis based on the infinitesimal model in which it is assumed that the trait is determined by an infinite number of unlinked loci without epistasis. The genetic variance is reduced not as a consequence of the genotypic frequency change but due to the build-up of linkage disequilibrium under truncation selection in this model. We assume that the order of the genotypic contribution from each locus is n(-1/2), where n is the number of loci involved, and investigate the change in linkage disequilibrium resulting from selection and self-fertilization using genotypic frequency dynamics in order to analyze the change in the genetic variance. Our analysis gives recurrence relations of genetic variance among the succeeding generations for the three cases of gene action, i.e., purely additive action, pure dominance without additive effect and the presence of both additive effect and dominance, respectively. Numerical examples are also given as a check on the recurrence formulas.