The Genetics of DROSOPHILA SUBOBSCURA Populations. Xvii. Further Genic Heterogeneity within Electromorphs by Urea Denaturation and the Effect of the Increased Genic Variability on Linkage Disequilibrium Studies

Urea denaturation of allozymes was used to provide finer resolution of allelic states within classes of different electrophoretic mobility. This method gives perfectly repeatable results. About 170 isogenic strains for the O chromosome of Drosophila subobscura, derived from two natural populations, were constructed. Their gene arrangements were studied, as well as eight polymorphic genes located on the O chromosome (Est-5, Odh, Ao, ME, Xdh, Lap, Pept-1 and Acph). Crosses performed indicate that differences in urea sensitivity are genetically controlled by the same genes that control electrophoretic mobility. Twice as many alleles have been detected in comparison to the usual electrophoretic method. However, the effective number of alleles did not increase considerably.Studies of linkage disequilibria, by taking into account the finer resolution of allelic states, gave results nearly identical with those obtained in studies where the usual electrophoretic method was used. Although the power of the test is diminished, the absence of genic associations seems to indicate that there are no hidden linkage disequilibria in electrophoretic studies (because of consolidation effects of real alleles into few electromorph classes). The paucity of linkage disequilibria would indicate that there are no epistatic interactions such as those suggested in the model of Franklin and Lewontin (1970).     

Disequilibrium among several linked neutral genes in finite population 1. mean changes in disequilibrium

Formulae are developed for computing changes in expected values in a finite population of linkage disequilibrium among neutral genes from more than two loci, although the exact analysis is taken up to only six loci. An essentially haploid model is used. As with two loci, the three-locus disequilibrium declines exponentially at all generations, but for m > 3 loci a matrix has to be constructed to give joint changes in the m-locus disequilibrium and products of disequilibria with fewer loci, for example of two $\text{m}\text{2}\text{-}\text{locus}$ disequilibria. The asymptotic rates of change in multilocus disequilibria depend on the arrangement of genes on the chromosome as well as its total length, but the initial rate of breakdown of disequilibrium from a line cross base is less dependent on the arrangement. With equally spaced loci the asymptotic rate of breakdown of m locus disequilibrium is roughly proportional to m. Although mutation and interference are excluded from the main analysis, it is shown how they can be incorporated.     

The Genetical Response to Natural Selection by Varied Environments. IV. Gametic Disequilibrium in Spatially Varied Environments

Gametic disequilibria between allozyme loci were related to spatial variation of the environment in caged populations of Drosophila melanogaster . Two experiments, one with flies collected at "Chateau Tahbilk," South Australia, and the other with flies from "Groningen," The Netherlands, were sampled at generations 16 and 32. Spatial variation of the environment was stimulated using three food media. Eight polymorphic allozyme loci were used to estimate gametic disequilibria from digenic combinations of allotypes. All populations were duplicated within an environment and maintained at about 2500 adults. Standardized gametic disequilibria were compared by a weighted least squares analysis of the z-transformed statistical correlation of allele frequencies. Gametic disequilibria were strongly dependent upon food niche and food-niche interactions. The effects also varied with sampling time and were similar in duplicate populations. Gametic disequilibria were most often detected in the "Groningen"-derived populations and their strength was not strongly associated with recombination fraction. Many of the disequilibria concerned unlinked loci. The strength of selection was probably considerable and populations were evolving genetic architectures which reflected niche selection by the different foods without marked genetic isolation between foods; gene frequencies did not vary between niches within a population cage.     

The Genetic Structure of Natural Populations of DROSOPHILA MELANOGASTER. Xii. Linkage Disequilibrium in a Large Local Population

Seven hundred and three second chromosomes were extracted from a Raleigh, North Carolina population of Drosophila melanogaster in 1970. Additionally, four hundred and eighty-nine third chromosomes were extracted from a large cage population founded from the flies in the 1970 Raleigh collection. The alpha glycerol-3-phosphate dehydrogenase-1, malate dehydrogenase-1, alcohol dehydrogenase, and alpha amylase loci were studied from the second chromosomes, and the esterase-6, esterase-C, and octanol dehydrogenase loci were analyzed from the third chromosomes. Inversions, relative viability and fecundity were studied for both classes of chromosomes. The following significant findings were obtained: (1) All loci examined were polymorphic or had at least two alleles at appreciable frequencies. Analysis of the combined data from this experiment with that of Mukai, Mettler and Chigusa (1971) revealed that the frequencies of the genes in the second chromosomes collected in early August were approximately the same over three years. (2) Linkage disequilibria between and among isozyme genes inter se were not detected except in a few cases which can be considered due to non-random sampling. (3) Linkage disequilibria between isozyme genes and polymorphic inversions were detected when the recombination values between the breakage points of the inversions and the genes in question were small. In only a few cases, were second and third order linkage disequilibria including polymorphic inversions detected. (4) Evidence for either variation among genotypes within loci or cumulative effects of heterozygosity was found for viability and fecundity. As a result of these findings, it was tentatively concluded that although selection might be perceptibly operating on some polymorphic isozyme loci, most of the polymorphic isozyme genes are selectively neutral or near-neutral in the populations studied.     

Disequilibrium among several linked neutral genes in finite population I. Mean changes in disequilibrium

Formulae are developed for computing changes in expected values in a finite population of linkage disequilibrium among neutral genes from more than two loci, although the exact analysis is taken up to only six loci. An essentially haploid model is used. As with two loci, the three-locus disequilibrium declines exponentially at all generations, but for m > 3 loci a matrix has to be constructed to give joint changes in the m-locus disequilibrium and products of disequilibria with fewer loci, for example of two m/2-locus disequilibria. The asymptotic rates of change in multilocus disequilibria depend on the arrangement of genes on the chromosome as well as its total length, but the initial rate of breakdown of disequilibrium from a line cross base is less dependent on the arrangement. With equally spaced loci the asymptotic rate of breakdown of m locus disequilibrium is roughly proportional to m. Although mutation and interference are excluded from the main analysis, it is shown how they can be incorporated.     

Gametic associations between inversion and allozyme polymorphisms in Drosophila buzzatii.

Gametic disequilibria between second chromosome polymorphic arrangements and seven linked allozyme loci were estimated in seven populations of Drosophila buzzatii from Argentina. Significant and consistent associations across populations were detected for Est-1, Est-2, Aldox, and XDH: Phenograms based on Nei's genetic distance showed that chromosomes carrying the 2ST arrangement were more similar to each other, irrespective of the population from which they were extracted, than to chromosomes carrying the derived 2J and 2JZ3.Restriction of recombination in heterokaryotypes seems to be the best explanation for the significant linkage disequilibria between inversions and the loci located inside the rearranged segments, for example, Est-1 and Aldox, or close to the break points, for example, Est-2. However, epistatic interactions between Xdh, which is outside the inversions and not near the break points, and loci tightly linked to the inversions, is the most likely explanation for the association between Xdh and chromosomal arrangements. Some of the associations detected in endemic Argentinean populations are coincident with data obtained in colonizing populations of the Old World and Australia. Thus historical processes that took place in the original area of the species' distribution can account for these linkage disequilibria in colonized populations of D. buzzatii.     

Genetic structure of five susceptibility gene regions for coronary artery disease: disequilibria within and among regions

We analyzed two-locus disequilibria for 16 polymorphic loci of seven susceptibility genes for coronary artery disease located in five chromosomal regions distributed across four chromosomes. Included were the genes coding for apolipoprotein B (ApoB, chromosome 2, four marker loci), lipoprotein lipase (LPL, chromosome 8, three marker loci), apolipoproteins AI, CIII, AIV (ApoAI–CIII–AIV, chromosome 11, three marker loci), apolipoprotein E (ApoE, chromosome 19, two marker loci), and the low density lipoprotein receptor (LDLR, chromosome 19, four marker loci). Our sample included 540 unrelated individuals from the Rochester, Minn. population. There were no statistically significant deviations of single-locus genotypes from Hardy-Weinberg equilibrium. The strongest associations within genes were for composite diallelic disequilibria; 17/19 were significant (13 at Pr <0.001, 1 at Pr <0.01, 3 at Pr <0.05). These observations suggest marker alleles within genes have a shared evolutionary history reflected by disequilibria that have not been dissipated by recombination. Disequilibrium was not generally concordant with the physical orderings of markers. Only two significant higher-order disequilibria were observed although 12 triallelic disequilibria were at maximum possible values. We observed 19 statistically significant disequilibria (Pr <0.05; 4 composite diallelic, 13 triallelic, and 2 quadriallelic) between 101 pairs of marker loci, where each locus in a pair was from a different unlinked region. These unexpected results are most likely explained by recent historical factors, including worldwide population expansion and amalgamation with continuous admixture, that influence the genetic structure (organization of alleles and non-alleles into genotypes) of a population. We conclude that disequilibria between loci from unlinked regions may be more extensive than is commonly assumed. Our findings also suggest that it is, on average, at least 15 times more likely to not detect significant disequilibrium among unlinked loci when it is really present than to make a false positive inference. Disequilibria between functional loci within or between regions will impact estimates of genetic variance associated with particular functional mutations within a susceptibility gene region. Received: 15 January 1998 / Accepted: 24 June 1998     

The Genetics of DROSOPHILA SUBOBSCURA Populations. IX. Studies on Linkage Disequilibrium in Four Natural Populations

Gametic frequencies were obtained in four natural populations of D. sub-obscura by extracting wild chromosomes and subsequently analyzing them for inversions and allozymes. The genes Lap and Pept-1, both located within the same inversions of chromosome O, were found in striking nonrandom associations with them of the same kind and degree in all populations studied. On the contrary, the gene Acph, also located within the previously mentioned inversions, was found in linkage disequilibrium with them only in two populations and of opposite directions. This is also the case for the genes Est-9 and Hk, both located within chromosome E inversions. While the gene Est-9 was in strong linkage disequilibrium with the inversions, of the same kind and degree in all populations studied, Hk was found to be in linkage equilibrium. Allele frequencies for the 29 genes studied do not show geographical variation except for the genes Lap, Pept-1 and Est-9, the ones found in linkage disequilibria with the geographically varying gene arrangements. Although mechanical or historical explanations for these equilibria cannot be ruled out, these data cannot be explained satisfactorily by the "middle gene explanation," which states that loci displaying such linkage disequilibria are the ones located near the break points of inversions, while the ones displaying linkage equilibria with them are located in the middle of them. There is no evidence for consistent linkage disequilibria between pairs of loci, except for the closely linked genes of the complex locus, Est-9. This would imply, if it is not a peculiarity of the Est-9 complex, that the linkage disequilibria are found only between very closely linked loci or that, for less closely linked genes, the associations are too weak to be detected by the usual samples sizes.     

The Genetic Structure of Natural Populations of DROSOPHILA MELANOGASTER Xiii. Further Studies on Linkage Disequilibrium

The Raleigh, North Carolina, population of Drosophila melanogaster was examined for linkage disequilibrium in 1974, several years after previous analyses in 1968, 1969, and 1970. alphaglycerol-3-phosphate dehydrogenase-1 (alphaGpdh-1), malate dehydrogenase-1 (Mdh-1), alcohol dehydrogenase (Adh), and hexokinase-C (Hex-C, tentative name, F. M. Johnson, unpublished; position determined by the present authors to be 2-74.5) were assayed for 617 second chromosomes, and esterase-C (Est-C) and octanol dehydrogenase (Odh) were assayed for 526 third chromosomes. In addition, two polymorphic inversions in the second chromosomes [In(2L)t and In(2R)NS] were examined, and the following findings were obtained: (1) No linkage disequilibrium between isozyme genes was detected. Significant linkage disequilibria were found only between the polymorphic inversions and isozyme genes [In(2L)t vs. Adh, and In(2R)NS vs. Hex-C]. Significant disequilibrium was not detected between In(2L)t and alphaGpdh-1, which is included in the inversion, but a tendency toward disequilibrium was consistently found from 1968 to 1974. The frequency of two-strand double crossovers within inversion In(2L)t involving a single crossover on each side of alphaGpdh-1 was estimated to be 0.00022. Thus, the consistent but not significant linkage disequilibrium between the two factors can be explained by recombination after the inversion occurred. (2) Previously existing linkage disequilibrium between Adh and In(2R)NS (the distance is about 30 cM, but the effective recombination value is about 1.75%) was found to have disappeared. (3) No higher-order linkage disequilibrium was detected. (4) Linkage disequilibrium between Odh and Est-C (the distance of which was estimated to be 0.0058 +/- 0.002) could not be detected (chi(2) (df=1) = 0.9).-From the above results, it was concluded that linkage disequilibria among isozyme genes are very rare in D. melanogaster, so that the Franklin-Lewontin model (Franklin and Lewontin 1970) is not applicable to these genes. The linkage disequilibria between some isozyme genes and polymorphic inversions may be explained by founder effect.     

Interlineage Mytilus galloprovincialis Lmk. 1819 hybridization yields inconsistent genetic outcomes in the Southern hemisphere

A Southern hemisphere lineage of the blue mussel Mytilus galloprovincialis has been diverging in allopatry from Northern hemisphere conspecifics for 0.84–1.2 million years. Secondary contact between Southern and Northern hemisphere mussels in Chile, New Zealand and Australia provides an opportunity to better understand the extent and consequences of extensive range expansion. Non-native M. galloprovincialis and hybrids, as detected from RFLP assays of nuclear and mitochondrial DNA, are present in all three countries and significant cytonuclear disequilibria exist for native homozygotes in Chile and New Zealand, non-native homozygotes in Chile and non-native heterozygotes in New Zealand. Introductions into Australia are rare events given that no pure non-native mussels were detected. Immigration from one or both taxa into the hybrid zone may underlie disequilibria in New Zealand, whilst gender-directional crossing with limited ongoing hybridization contributes to disequilibria in Chile. Hybridization dynamics do not pose a threat to the Southern lineage in Chile and Australia, but in New Zealand, introgression, continued immigration and slight hybridization gender bias towards non-native maternal parents could lead to the regional extirpation of the native lineage.     

Comparative Effects of Pollen and Seed Migration on the Cytonuclear Structure of Plant Populations. I. Maternal Cytoplasmic Inheritance

We explicitly solve and analyze a series of deterministic continent-island models to delimit the effects of pollen and seed migration on cytonuclear frequencies and disequilibria in random-mating, mixed-mating and self-fertilized populations. Given the critical assumption of maternal cytoplasmic inheritance, five major findings are (i) nonzero cytonuclear disequilibria will be maintained in the island population if and only if at least some migration occurs each generation through seeds with nonrandom cytonuclear associations; (ii) immigrant seeds with no cytonuclear disequilibria can strongly affect the genetic structure of the island population by generating significant and long-lasting transient associations; (iii) with all else being equal, substantially greater admixture disequilibria are generally found with higher rates of seed migration into, or higher levels of self-fertilization within, the island population (with the possible exception of the heterozygote disequilibrium); (iv) pollen migration can either enhance or reduce the cytonuclear disequilibria caused by seed migration, or that due to mixed-mating in the absence of seed migration, but the effect is usually small and appears primarily to make a noticeable difference in predominantly outcrossing populations; and (v) pollen migration alone cannot generate even transient disequilibria de novo in populations with completely random associations. This same basic behavior is exhibited as long as there is some random outcrossing in the island population. Self-fertilized populations represent a special case, however, in that they are necessarily closed to pollen migration, and nonzero disequilibria can be maintained even in the absence of seed migration. All of these general results hold whether the population is censused as adults or as seeds, but the ability to detect nonrandom cytonuclear associations can depend strongly on the life stage censused in populations with a significant level of random outcrossing. We suggest how these models might be used for the estimation of seed and pollen migration.     

Exploring a Local Genetic Interaction Network Using Evolutionary Replay Experiments

Understanding how genes interact is a central challenge in biology. Experimental evolution provides a useful, but underutilized, tool for identifying genetic interactions, particularly those that involve non-loss-of-function mutations or mutations in essential genes. We previously identified a strong positive genetic interaction between specific mutations in KEL1 (P344T) and HSL7 (A695fs) that arose in an experimentally evolved Saccharomyces cerevisiae population. Because this genetic interaction is not phenocopied by gene deletion, it was previously unknown. Using “evolutionary replay” experiments, we identified additional mutations that have positive genetic interactions with the kel1-P344T mutation. We replayed the evolution of this population 672 times from six timepoints. We identified 30 populations where the kel1-P344T mutation reached high frequency. We performed whole-genome sequencing on these populations to identify genes in which mutations arose specifically in the kel1-P344T background. We reconstructed mutations in the ancestral and kel1-P344T backgrounds to validate positive genetic interactions. We identify several genetic interactors with KEL1, we validate these interactions by reconstruction experiments, and we show these interactions are not recapitulated by loss-of-function mutations. Our results demonstrate the power of experimental evolution to identify genetic interactions that are positive, allele specific, and not readily detected by other methods, shedding light on an underexplored region of the yeast genetic interaction network.     

Prediction of C. elegans Longevity Genes by Human and Worm Longevity Networks

Intricate and interconnected pathways modulate longevity, but screens to identify the components of these pathways have not been saturating. Because biological processes are often executed by protein complexes and fine-tuned by regulatory factors, the first-order protein-protein interactors of known longevity genes are likely to participate in the regulation of longevity. Data-rich maps of protein interactions have been established for many cardinal organisms such as yeast, worms, and humans. We propose that these interaction maps could be mined for the identification of new putative regulators of longevity. For this purpose, we have constructed longevity networks in both humans and worms. We reasoned that the essential first-order interactors of known longevity-associated genes in these networks are more likely to have longevity phenotypes than randomly chosen genes. We have used C. elegans to determine whether post-developmental inactivation of these essential genes modulates lifespan. Our results suggest that the worm and human longevity networks are functionally relevant and possess a high predictive power for identifying new longevity regulators.     

Disequilibrium among several linked neutral genes in finite population. II. Variances and covariances of disequilibria

A method is derived for computing the variances and covariances of linkage disequilibria between neutral genes in finite populations, which is based on a linear transformation of results given previously for the mean values of disequilibria. The formulae obtained are limited to moments of sixth order or less, such as the variance of the three-locus disequilibrium. It is shown that there is no covariance between any pair of disequilibria in populations starting equilibrium. The pattern of change with time in variance of the three-locus disequilibrium from populations initially in equilibrium is similar to that for two loci, except that the highest values are achieved rather earlier and are smaller.     

Linkage Disequilibrium in Natural and Experimental Populations of Drosophila Melanogaster

We have studied linkage disequilibrium in Drosophila melanogaster in two samples from a wild population and in four large laboratory populations derived from the wild samples. We have assayed four polymorphic enzyme loci, fairly closely linked in the third chromosome: Sod Est-6, Pgm, and Odh. The assay method used allows us to identify the allele associations separately in each of the two homologous chromosomes from each male sampled. We have detected significant linkage disequilibrium between two loci in 16.7% of the cases in the wild samples and in 27.8% of the cases in the experimental populations, considerably more than would be expected by chance alone. We have also found three-locus disequilibria in more instances than would be expected by chance. Some disequilibria present in the wild samples disappear in the experimental populations derived from them, but new ones appear over the generations. The effective population sizes required to generate the observed disequilibria by randomness range from 40 to more than 60,000 individuals in the natural population, depending on which locus pair is considered, and from 100 to more than 60,000 in the experimental populations. These population sizes are unrealistic; the fact that different locus-pairs yield disparate estimates within the same population argues against the likelihood that the disequilibria may have arisen as a consequence of population bottlenecks. Migration, or population mixing, cannot be excluded as the process generating the disequilibria in the wild samples, but can in the experimental populations. We conclude that linkage disequilibrium in these populations is most likely due to natural selection acting on the allozymes, or on loci very tightly linked to them.     

The Effects of Assortative Mating and Migration on Cytonuclear Associations in Hybrid Zones

We examine the influence of nonrandom mating and immigration on the evolutionary dynamics of cytonuclear associations in hybrid zones. Recursion equations for allelic and genotypic cytonuclear disequilibria were generated under models of (1) migration alone, assuming hybrid zone matings are random with respect to cytonuclear genotype; and (2) migration in conjunction with refined epistatic mating, in which females of the pure parental species preferentially mate with conspecific males. Major results are as follows: (a) even the slightest migration removes the dependency of the final outcome on initial conditions, producing a unique equilibrium in which both pure parental genotypes are maintained in the hybrid zone; (b) in contrast to nuclear genes, the dynamics of cytoplasmic allele frequencies appear robust to changes in the assumed mating system, yet are particularly sensitive to gene flow; (c) continued immigration can generate permanent cytonuclear disequilibria, whether mating is random or assortative; and (d) the order of population censusing (before versus after reproduction by immigrants) can have a dramatic effect on the magnitude but not the pattern of cytonuclear disequilibria. Using the maximum likelihood method, the parameter space of migration rates and assortative mating rates was examined for best fit to observed cytonuclear disequilibria data in a hybrid population of Hyla tree frogs. An epistatic mating model with a total immigration rate of about 32% per generation produces equilibrium gene frequencies and cytonuclear disequilibria consistent with the empirical observations.     

Multivariate Analysis of Gametic Disequilibrium in the Yanomama

The gametic disequilibria between all possible pairs of loci were examined for a set of eight codominant loci in each of fifty Yanomama villages, using a multivariate correlation analysis which reduces the results to a single measure of departure from multiple-locus-gametic equilibrium. Thirty-two of the fifty villages departed significantly from multiple-locus gametic equilibrium. The largest contributions to the departure from multiple-locus equilibrium were due to the disequilibria between MN and Ss and between Rh(Cc) and Rh(Ee), indicating the effects of tight linkage. After removing the effects of these obvious sources of disequilibrium, sixteen of the fifty villages still remained significantly out of equilibrium. The disequilibrium between any particular pair of loci was highly erratic from village to village, and (with the exception of the MN-Ss and Cc-Ee disequilibria) averaged out very close to zero overall, suggesting a lack of systematic forces (epistatic selection). The departure from equilibrium in any one village is in excess of that expected from random sampling alone, and is attributed primarily to the fission-fusion mode of village formation operative in the Yanomama and the fact that a single village consists of a few extended lineages. Village allele frequencies are highly correlated across loci, and most of the non-independence is accounted for by large correlations in the average allelic frequencies of different loci for related villages. It is suggested that these correlations also are due to territorial expansion and population growth. For the tribe as a whole, all but the tightly linked markers of the MNSs and Rh complexes are approximately uncorrelated, and large departures from multiple-locus Hardy-Weinberg expectation are primarily due to substantial Wahlund variance within the tribe. There is no need to postulate a role for selection in these disequilibria.     

The effects of admixture and population subdivision on cytonuclear disequilibria

We examine the generation of cytonuclear disequilibria by admixture and continued gene flow. General formulas analogous to the nuclear case are first derived showing that the allelic and genotypic disequilibria from admixture or population subdivision equal their expected value across the contributing (sub) populations plus the covariance across these sources between the cytoplasmic gene frequency and the relevant nuclear frequency. A detailed study is then presented of the cytonuclear dynamics, in a random-mating population under two different migration scenarios. In both cases closed-form solutions are given for all variables as a function of the initial conditions and relevant migration parameters. The dynamics of the gene frequencies and allelic disequilibria, which dominate each system, are the same as those involving two unlinked nuclear loci, while the dynamics of the genotypic disequilibria and cytonuclear frequencies have no nuclear counterpart. The continent-island formulation focuses on a population receiving continued immigration from a large source of constant composition. A major discovery is that cytonuclear disequilibria can transiently build up on the "island" to levels far exceeding those found at equilibrium. In contrast, the admixture formulation focuses on the dynamics within two populations undergoing continued intermigration. Although in this case all cytonuclear associations must ultimately decay to zero, long-term transient disequilibria can develop which are many times their initial admixture values. For both migration scenarios it is shown that the time of population censusing relative to migration and reproduction dramatically affects both the amount and pattern of the nonrandom associations produced. The empirical relevance of these models is discussed in light of nuclear-mitochondrial data from a hybrid zone between European and North American eels and from a zone of racial admixture in humans.     

The genetic architecture of complex traits in teosinte (Zea mays ssp. parviglumis): new evidence from association mapping

Previous association analyses showed that variation at major regulatory genes contributes to standing variation for complex traits in Balsas teosinte, the progenitor of maize. This study expands our previous association mapping effort in teosinte by testing 123 markers in 52 candidate genes for association with 31 traits in a population of 817 individuals. Thirty-three significant associations for markers from 15 candidate genes and 10 traits survive correction for multiple testing. Our analyses suggest several new putative causative relationships between specific genes and trait variation in teosinte. For example, two ramosa genes (ra1 and ra2) associate with ear structure, and the MADS-box gene, zagl1, associates with ear shattering. Since zagl1 was previously shown to be a target of selection during maize domestication, we suggest that this gene was under selection for its effect on the loss of ear shattering, a key domestication trait. All observed effects were relatively small in terms of the percentage of phenotypic variation explained (<10%). We also detected several epistatic interactions between markers in the same gene that associate with the same trait. Candidate-gene-based association mapping appears to be a promising method for investigating the inheritance of complex traits in teosinte.