A complete enumeration and classification of two-locus disease models

There are 512 two-locus, two-allele, two-phenotype, fully penetrant disease models. Using the permutation between two alleles, between two loci, and between being affected and unaffected, one model can be considered to be equivalent to another model under the corresponding permutation. These permutations greatly reduce the number of two-locus models in the analysis of complex diseases. This paper determines the number of nonredundant two-locus models (which can be 102, 100, 96, 51, 50, or 58, depending on which permutations are used, and depending on whether zero-locus and single-locus models are excluded). Whenever possible, these nonredundant two-locus models are classified by their property. Besides the familiar features of multiplicative models (logical AND), heterogeneity models (logical OR), and threshold models, new classifications are added or expanded: modifying-effect models, logical XOR models, interference and negative interference models (neither dominant nor recessive), conditionally dominant/recessive models, missing lethal genotype models, and highly symmetric models. The following aspects of two-locus models are studied: the marginal penetrance tables at both loci, the expected joint identity-by-descent (IBD) probabilities, and the correlation between marginal IBD probabilities at the two loci. These studies are useful for linkage analyses using single-locus models while the underlying disease model is two-locus, and for correlation analyses using the linkage signals at different locations obtained by a single-locus model.     

Evolution of coadaptation in a two-locus epistatic system

Although recent advances in genome biology have dramatically increased our understanding of the contribution of gene interactions to the development of complex phenotypes, we still lack general agreement on the process and mechanisms responsible for the evolution of epistatic systems. Even if genes in a species are indeed integrated into coadapted complexes of interacting components, simple additive evolution may eventually result in epistatic differentiation of populations. Consequently, the prevalence of epistatic gene action does not tell us anything about the role of epistatic selection in the history of population divergence. To elucidate the contribution of epistatic selection in the evolution of coadaptation, we investigate the fixation process of two mutations that interact synergistically to enhance fitness. We show by diffusion analysis and simulations that epistatic selection on cosegregating variants does not by itself promote the evolution of epistatic systems; rather, accumulation of neutral mutations may play a crucial role, creating an appropriate genetic milieu for adaptive evolution in the future generations.     

Some epistatic two-locus models of disease. II. The confounding of linkage and association

This paper continues to examine the model discussed in the preceding paper. Specifically, it will be shown how a linkage analysis performed in the presence of a disease-marker association can give rise to erroneous and misleading results.     

Some epistatic two-locus models of disease. I. Relative risks and identity-by-descent distributions in affected sib pairs

A two-locus disease model is presented in which a marker locus interacts epistatically with another unlinked trait to cause the disease. Such a model can lead to disease-marker associations and distortions in the sharing of marker types among affected family members. These effects are quantified. In the case of HLA-disease associations, this model is presented as an alternative to the “hitchhiking” theory of tight linkage leading to linkage disequilibrium.     

A simple method for testing two-locus models of inheritance

A graphic method for testing simple two-locus models of inheritance is developed. The model assumes two alleles at each locus where both loci exhibit dominant, both exhibit recessive, or one locus exhibits dominant and one locus exhibits recessive inheritance. Examples of applying the graphs using published data on three diseases are given.     

Two-stage two-locus models in genome-wide association

Studies in model organisms suggest that epistasis may play an important role in the etiology of complex diseases and traits in humans. With the era of large-scale genome-wide association studies fast approaching, it is important to quantify whether it will be possible to detect interacting loci using realistic sample sizes in humans and to what extent undetected epistasis will adversely affect power to detect association when single-locus approaches are employed. We therefore investigated the power to detect association for an extensive range of two-locus quantitative trait models that incorporated varying degrees of epistasis. We compared the power to detect association using a single-locus model that ignored interaction effects, a full two-locus model that allowed for interactions, and, most important, two two-stage strategies whereby a subset of loci initially identified using single-locus tests were analyzed using the full two-locus model. Despite the penalty introduced by multiple testing, fitting the full two-locus model performed better than single-locus tests for many of the situations considered, particularly when compared with attempts to detect both individual loci. Using a two-stage strategy reduced the computational burden associated with performing an exhaustive two-locus search across the genome but was not as powerful as the exhaustive search when loci interacted. Two-stage approaches also increased the risk of missing interacting loci that contributed little effect at the margins. Based on our extensive simulations, our results suggest that an exhaustive search involving all pairwise combinations of markers across the genome might provide a useful complement to single-locus scans in identifying interacting loci that contribute to moderate proportions of the phenotypic variance.     

Neutral two-locus multiple allele models with recombination

General formulae for the homozygosity and variance of linkage disequilibrium are derived for neutral, stationary, two-locus multiple allele models where there is a symmetric type of mutation at each locus. Particular cases examined are K allele models, the infinite alleles model, and the stepwise mutation model. The two-locus infinite allele model is examined at the molecular level and a joint probability generating function is found for the number of heterozygous sites at each locus in two randomly chosen gametes.     

Detecting purely epistatic multi-locus interactions by an omnibus permutation test on ensembles of two-locus analyses

Background  

Purely epistatic multi-locus interactions cannot generally be detected via single-locus analysis in case-control studies of complex diseases. Recently, many two-locus and multi-locus analysis techniques have been shown to be promising for the epistasis detection. However, exhaustive multi-locus analysis requires prohibitively large computational efforts when problems involve large-scale or genome-wide data. Furthermore, there is no explicit proof that a combination of multiple two-locus analyses can lead to the correct identification of multi-locus interactions.     

A semi-symmetric two-locus model

The two-locus symmetric viability model characterized by its invariance with respect to the exchange of alleles at each locus, is a well-studied model of classical two-locus theory. The symmetric model introduced by Lewontin and Kojima is among the few multi-locus models with epistatic interactions between loci for which a polymorphism with linkage equilibrium can be stable and this happens when recombination is sufficiently large. We show that an analogous property holds true for a different model, in which symmetry need exist at only one locus. The properties of this new semi-symmetric model are compared with those of the classical symmetric model. For tight linkage, two classes of polymorphisms are possible, depending on the magnitude of additive epistasis. The recombination rate above which linkage equilibrium becomes stable is derived analytically. As in the symmetric model, intervals of recombination in which no polymorphism is stable are possible, and stable polymorphisms can coexist with stable fixations.     

Four simultaneously stable polymorphic equilibria in two-locus two-allele models

The existence of four simultaneously stable equilibria with both loci polymorphic is shown for the Lewontin-Kojima version of the two-locus two-allele symmetric viability model, using bifurcation theory. This exceeds the previously claimed bound of two stable polymorphisms. Biological implications of the result are discussed.     

Simulation studies of segregation analysis: application to two-locus models

We tested the power of a segregation analysis method (first proposed by Elandt-Johnson) to distinguish between single-locus and two-locus models, with and without environmentally caused reduced penetrance. We also looked at the effect of ascertainment probability on the analysis and at the proband-conditioned ascertainment correction proposed by Cannings and Thompson. We found that: (1) the segregation analysis has sufficient power to distinguish between the fully-penetrant double-recessive (RR) model and the fully-penetrant single-locus dominant and recessive models; (2) the method can also distinguish fairly well between the dominant-recessive (DR) and RR models, even when one does not take into account the population prevalence; (3) the method has much less power to distinguish between the fully-penetrant RR model and the single-locus models with reduced penetrance; (4) when environmental penetrance is taken account of in the analysis, the power of the method to distinguish between the one- and two-locus models improved substantially; (5) the estimates of ascertainment probability, pi, were robust, regardless of the model under which the data were generated; and (6) the Cannings-Thompson approach to ascertainment correction worked well only when the pi used to generate the data was less than .1.     

Rates of decay of linkage disequilibrium under two-locus models of selection

The effect of selection and linkage on the decay of linkage disequilibrium, D, is investigated for a hierarchy of two-locus models. The method of analysis rests upon a qualitative classification of the dynamic of D under selection relative to the neutral dynamic. To eliminate the confounding effects of gene frequency change, the behavior of D is first studied with gene frequencies fixed at their invariant values. Second, the results are extended to certain special situations where gene frequencies are changing simultaneously.A wide variety of selection regimes can cause an acceleration of the rate of decay of D relative to the neutral rate. Specifically, the asymptotic rate of decay is always faster than the neutral rate in the neighborhood of a stable equilibrium point, when viabilities are additive or only one locus is selected. This is not necessarily the case for models in which there is nonzero additive epistasis. With multiplicative viabilities, decay is always accelerated near a stable boundary equilibrium, but decay is only faster near the stable central equilibrium (with = 0) if linkage is sufficiently loose. In the symmetric viability model, decay may even be retarded near a stable boundary equilibrium. Decay is only accelerated near a stable corner equilibrium when the double homozygote is more fit than the double heterozygotes. Decay near a stable edge equilibrium may be retarded if there is loose linkage. With symmetric viabilities there is usually an acceleration of the decay process for gene frequencies near 1/2 when the central equilibrium (with = 0) is stable. This is always the case when the sign of the epistasis is negative or zero.Conversely, the decay ofD is retarded in the neighborhood of a stable equilibrium in the multiplicative and symmetric viability models if any of the conditions above are violated. Near an unstable equilibrium of any of the models considered,D may either increase or decay at a rate slower than, equal to, or faster than the neutral rate. These analytic results are supplemented by numerical studies of the symmetric viability model.     

Properties of the major gene index for three-allele and two-locus models

The expected value of the MGI(2) statistic was evaluated for several three-allele and two-locus models for which we have specified dominance relationships among phenotype classes based on heterozygote and homozygote genotype groupings. In a parallel simulation study we investigated the nature of MGI(alpha) for alpha = 1/2, 1, 2 under more elaborate continuous trait expressions incurred from superposition of a background distribution upon each of the phenotypic mean effects.     

Cytonuclear models of epistatic mating with backcrossing and selection

We develop hybrid zone models that explore the combined effects of mating system and either backcrossing or viability selection on the disequilibria between nuclear and cytoplasmic genes. In the epistatic mating plus backcrossing model, we find patterns of permanent cytonuclear disequilibria like those found when epistatic mating is the only factor, as well as a novel combination of significant cytonuclear disequilibria sign patterns. The second group of models evaluates the potential of epistatic mating and postzygotic viability selection to maintain cytonuclear disequilibria. Simulations are used to evaluate nine patterns of selection, each of which represent differing forms of selection against hybrids, and show that while all disequilibria usually decay to zero, under certain circumstances a number of different patterns of significant cytonuclear disequilibria are possible at equilibrium. The results from these models are compared to the observed cytonuclear disequilibria previously found in a hybrid population of Hyla treefrogs.     

A two-locus model of speciation.

Speciation is considered as the evolution of partial or complete cross-incompatibility between the carriers of genes (at a locus called "object locus") that distinguish the prospective species populations. The mating relations at the object locus are modified by the alleles at a second mating modifier locus. Based on a widely applicable concept of fitness and mating preference, it is shown that heterozygote disadvantage in fitness at the object locus is necessary for speciation, which corroborates Wallace's hypothesis. It is pointed out that the difference between sympatric and parapatric speciation essentially lies in the mechanisms stabilizing the polymorphism required at the object locus as a prerequisite for speciation. In the presence of recombination between the object and mating modifier locus speciation may be prevented by forces maintaining gametic phase imbalance between these loci such as can result from unidirectional gene flow between parapatric populations.     

A complete solution for dissecting pure main and epistatic effects of QTL in triple testcross design

Epistasis plays an important role in genetics, evolution and crop breeding. To detect the epistasis, triple test cross (TTC) design had been developed several decades ago. Classical procedures for the TTC design use only linear transformations Z(1), Z(2) and Z(3), calculated from the TTC family means of quantitative trait, to infer the nature of the collective additive, dominance and epistatic effects of all the genes. Although several quantitative trait loci (QTL) mapping approaches in the TTC design have been developed, these approaches do not provide a complete solution for dissecting pure main and epistatic effects. In this study, therefore, we developed a two-step approach to estimate all pure main and epistatic effects in the F(2)-based TTC design under the F(2) and F(∞) metric models. In the first step, with Z(1) and Z(2) the augmented main and epistatic effects in the full genetic model that simultaneously considered all putative QTL on the whole genome were estimated using empirical Bayes approach, and with Z(3) three pure epistatic effects were obtained using two-dimensional genome scans. In the second step, the three pure epistatic effects obtained in the first step were integrated with the augmented epistatic and main effects for the further estimation of all other pure effects. A series of Monte Carlo simulation experiments has been carried out to confirm the proposed method. The results from simulation experiments show that: 1) the newly defined genetic parameters could be rightly identified with satisfactory statistical power and precision; 2) the F(2)-based TTC design was superior to the F(2) and F(2:3) designs; 3) with Z(1) and Z(2) the statistical powers for the detection of augmented epistatic effects were substantively affected by the signs of pure epistatic effects; and 4) with Z(3) the estimation of pure epistatic effects required large sample size and family replication number. The extension of the proposed method in this study to other base populations was further discussed.