Genetic Analysis of Central Carbon Metabolism Unveils an Amino Acid Substitution That Alters Maize NAD-Dependent Isocitrate Dehydrogenase Activity

Background

Central carbon metabolism (CCM) is a fundamental component of life. The participating genes and enzymes are thought to be structurally and functionally conserved across and within species. Association mapping utilizes a rich history of mutation and recombination to achieve high resolution mapping. Therefore, applying association mapping in maize (Zea mays ssp. mays), the most diverse model crop species, to study the genetics of CCM is a particularly attractive system.

Methodology/Principal Findings

We used a maize diversity panel to test the CCM functional conservation. We found heritable variation in enzyme activity for every enzyme tested. One of these enzymes was the NAD-dependent isocitrate dehydrogenase (IDH, E.C. 1.1.1.41), in which we identified a novel amino-acid substitution in a phylogenetically conserved site. Using candidate gene association mapping, we identified that this non-synonymous polymorphism was associated with IDH activity variation. The proposed mechanism for the IDH activity variation includes additional components regulating protein level. With the comparison of sequences from maize and teosinte (Zea mays ssp. Parviglumis), the maize wild ancestor, we found that some CCM genes had also been targeted for selection during maize domestication.

Conclusions/Significance

Our results demonstrate the efficacy of association mapping for dissecting natural variation in primary metabolic pathways. The considerable genetic diversity observed in maize CCM genes underlies heritable phenotypic variation in enzyme activities and can be useful to identify putative functional sites.     

Joint QTL linkage mapping for multiple-cross mating design sharing one common parent

Background

Nested association mapping (NAM) is a novel genetic mating design that combines the advantages of linkage analysis and association mapping. This design provides opportunities to study the inheritance of complex traits, but also requires more advanced statistical methods. In this paper, we present the detailed algorithm of a QTL linkage mapping method suitable for genetic populations derived from NAM designs. This method is called joint inclusive composite interval mapping (JICIM). Simulations were designed on the detected QTL in a maize NAM population and an Arabidopsis NAM population so as to evaluate the efficiency of the NAM design and the JICIM method.

Principal Findings

Fifty-two QTL were identified in the maize population, explaining 89% of the phenotypic variance of days to silking, and nine QTL were identified in the Arabidopsis population, explaining 83% of the phenotypic variance of flowering time. Simulations indicated that the detection power of these identified QTL was consistently high, especially for large-effect QTL. For rare QTL having significant effects in only one family, the power of correct detection within the 5 cM support interval was around 80% for 1-day effect QTL in the maize population, and for 3-day effect QTL in the Arabidopsis population. For smaller-effect QTL, the power diminished, e.g., it was around 50% for maize QTL with an effect of 0.5 day. When QTL were linked at a distance of 5 cM, the likelihood of mapping them as two distinct QTL was about 70% in the maize population. When the linkage distance was 1 cM, they were more likely mapped as one single QTL at an intermediary position.

Conclusions

Because it takes advantage of the large genetic variation among parental lines and the large population size, NAM is a powerful multiple-cross design for complex trait dissection. JICIM is an efficient and specialty method for the joint QTL linkage mapping of genetic populations derived from the NAM design.     

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.     

Genetic dynamics underlying phenotypic development of biomass yield in triticale

Background

The nature of dynamic traits with their phenotypic plasticity suggests that they are under the control of a dynamic genetic regulation. We employed a precision phenotyping platform to non-invasively assess biomass yield in a large mapping population of triticale at three developmental stages.

Results

Using multiple-line cross QTL mapping we identified QTL for each of these developmental stages which explained a considerable proportion of the genotypic variance. Some QTL were identified at each developmental stage and thus contribute to biomass yield throughout the studied developmental phases. Interestingly, we also observed QTL that were only identified for one or two of the developmental stages illustrating a temporal contribution of these QTL to the trait. In addition, epistatic QTL were detected and the epistatic interaction landscape was shown to dynamically change with developmental progression.

Conclusions

In summary, our results reveal the temporal dynamics of the genetic architecture underlying biomass accumulation in triticale and emphasize the need for a temporal assessment of dynamic traits.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-458) contains supplementary material, which is available to authorized users.     

Limits on the reproducibility of marker associations with southern leaf blight resistance in the maize nested association mapping population

Background

A previous study reported a comprehensive quantitative trait locus (QTL) and genome wide association study (GWAS) of southern leaf blight (SLB) resistance in the maize Nested Association Mapping (NAM) panel. Since that time, the genomic resources available for such analyses have improved substantially. An updated NAM genetic linkage map has a nearly six-fold greater marker density than the previous map and the combined SNPs and read-depth variants (RDVs) from maize HapMaps 1 and 2 provided 28.5 M genomic variants for association analysis, 17 fold more than HapMap 1. In addition, phenotypic values of the NAM RILs were re-estimated to account for environment-specific flowering time covariates and a small proportion of lines were dropped due to genotypic data quality problems. Comparisons of original and updated QTL and GWAS results confound the effects of linkage map density, GWAS marker density, population sample size, and phenotype estimates. Therefore, we evaluated the effects of changing each of these parameters individually and in combination to determine their relative impact on marker-trait associations in original and updated analyses.

Results

Of the four parameters varied, map density caused the largest changes in QTL and GWAS results. The updated QTL model had better cross-validation prediction accuracy than the previous model. Whereas joint linkage QTL positions were relatively stable to input changes, the residual values derived from those QTL models (used as inputs to GWAS) were more sensitive, resulting in substantial differences between GWAS results. The updated NAM GWAS identified several candidate genes consistent with previous QTL fine-mapping results.

Conclusions

The highly polygenic nature of resistance to SLB complicates the identification of causal genes. Joint linkage QTL are relatively stable to perturbations of data inputs, but their resolution is generally on the order of tens or more Mbp. GWAS associations have higher resolution, but lower power due to stringent thresholds designed to minimize false positive associations, resulting in variability of detection across studies. The updated higher density linkage map improves QTL estimation and, along with a much denser SNP HapMap, greatly increases the likelihood of detecting SNPs in linkage with causal variants. We recommend use of the updated genetic resources and results but emphasize the limited repeatability of small-effect associations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1068) contains supplementary material, which is available to authorized users.     

Genome-wide QTL mapping of nine body composition and bone mineral density traits in pigs

Background

Since the pig is one of the most important livestock animals worldwide, mapping loci that are associated with economically important traits and/or traits that influence animal welfare is extremely relevant for efficient future pig breeding. Therefore, the purpose of this study was a genome-wide mapping of quantitative trait loci (QTL) associated with nine body composition and bone mineral traits: absolute (Fat, Lean) and percentage (FatPC, LeanPC) fat and lean mass, live weight (Weight), soft tissue X-ray attenuation coefficient (R), absolute (BMC) and percentage (BMCPC) bone mineral content and bone mineral density (BMD).

Methods

Data on the nine traits investigated were obtained by Dual-energy X-ray absorptiometry for 551 pigs that were between 160 and 200 days old. In addition, all pigs were genotyped using Illumina’s PorcineSNP60 Genotyping BeadChip. Based on these data, a genome-wide combined linkage and linkage disequilibrium analysis was conducted. Thus, we used 44 611 sliding windows that each consisted of 20 adjacent single nucleotide polymorphisms (SNPs). For the middle of each sliding window a variance component analysis was carried out using ASReml. The underlying mixed linear model included random QTL and polygenic effects, with fixed effects of sex, housing, season and age.

Results

Using a Bonferroni-corrected genome-wide significance threshold of P < 0.001, significant peaks were identified for all traits except BMCPC. Overall, we identified 72 QTL on 16 chromosomes, of which 24 were significantly associated with one trait only and the remaining with more than one trait. For example, a QTL on chromosome 2 included the highest peak across the genome for four traits (Fat, FatPC, LeanPC and R). The nearby gene, ZNF608, is known to be associated with body mass index in humans and involved in starvation in Drosophila, which makes it an extremely good candidate gene for this QTL.

Conclusions

Our QTL mapping approach identified 72 QTL, some of which confirmed results of previous studies in pigs. However, we also detected significant associations that have not been published before and were able to identify a number of new and promising candidate genes, such as ZNF608.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0068-2) contains supplementary material, which is available to authorized users.     

Genetic Architecture of Ear Fasciation in Maize (Zea mays) under QTL Scrutiny

Maize ear fasciation

Knowledge of the genes affecting maize ear inflorescence may lead to better grain yield modeling. Maize ear fasciation, defined as abnormal flattened ears with high kernel row number, is a quantitative trait widely present in Portuguese maize landraces.

Material and Methods

Using a segregating population derived from an ear fasciation contrasting cross (consisting of 149 F_2:3 families) we established a two location field trial using a complete randomized block design. Correlations and heritabilities for several ear fasciation-related traits and yield were determined. Quantitative Trait Loci (QTL) involved in the inheritance of those traits were identified and candidate genes for these QTL proposed.

Results and Discussion

Ear fasciation broad-sense heritability was 0.73. Highly significant correlations were found between ear fasciation and some ear and cob diameters and row number traits. For the 23 yield and ear fasciation-related traits, 65 QTL were identified, out of which 11 were detected in both environments, while for the three principal components, five to six QTL were detected per environment. Detected QTL were distributed across 17 genomic regions and explained individually, 8.7% to 22.4% of the individual traits or principal components phenotypic variance. Several candidate genes for these QTL regions were proposed, such as bearded-ear1, branched silkless1, compact plant1, ramosa2, ramosa3, tasselseed4 and terminal ear1. However, many QTL mapped to regions without known candidate genes, indicating potential chromosomal regions not yet targeted for maize ear traits selection.

Conclusions

Portuguese maize germplasm represents a valuable source of genes or allelic variants for yield improvement and elucidation of the genetic basis of ear fasciation traits. Future studies should focus on fine mapping of the identified genomic regions with the aim of map-based cloning.     

Linkage disequilibrium and genome-wide association analysis for anthocyanin pigmentation and fruit color in eggplant

Background

The genome-wide association (GWA) approach represents an alternative to biparental linkage mapping for determining the genetic basis of trait variation. Both approaches rely on recombination to re-arrange the genome, and seek to establish correlations between phenotype and genotype. The major advantages of GWA lie in being able to sample a much wider range of the phenotypic and genotypic variation present, in being able to exploit multiple rounds of historical recombination in many different lineages and to include multiple accessions of direct relevance to crop improvement.

Results

A 191 accessions eggplant (Solanum melongena L.) association panel, comprising a mixture of breeding lines, old varieties and landrace selections originating from Asia and the Mediterranean Basin, was SNP genotyped and scored for anthocyanin pigmentation and fruit color at two locations over two years. The panel formed two major clusters, reflecting geographical provenance and fruit type. The global level of linkage disequilibrium was 3.4 cM. A mixed linear model appeared to be the most appropriate for GWA. A set of 56 SNP locus/phenotype associations was identified and the genomic regions harboring these loci were distributed over nine of the 12 eggplant chromosomes. The associations were compared with the location of known QTL for the same traits.

Conclusion

The GWA mapping approach was effective in validating a number of established QTL and, thanks to the wide diversity captured by the panel, was able to detect a series of novel marker/trait associations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-896) contains supplementary material, which is available to authorized users.     

Do large effect QTL fractionate? A case study at the maize domestication QTL teosinte branched1

Quantitative trait loci (QTL) mapping is a valuable tool for studying the genetic architecture of trait variation. Despite the large number of QTL studies reported in the literature, the identified QTL are rarely mapped to the underlying genes and it is usually unclear whether a QTL corresponds to one or multiple linked genes. Similarly, when QTL for several traits colocalize, it is usually unclear whether this is due to the pleiotropic action of a single gene or multiple linked genes, each affecting one trait. The domestication gene teosinte branched1 (tb1) was previously identified as a major domestication QTL with large effects on the differences in plant and ear architecture between maize and teosinte. Here we present the results of two experiments that were performed to determine whether the single gene tb1 explains all trait variation for its genomic region or whether the domestication QTL at tb1 fractionates into multiple linked QTL. For traits measuring plant architecture, we detected only one QTL per trait and these QTL all mapped to tb1. These results indicate that tb1 is the sole gene for plant architecture traits that segregates in our QTL mapping populations. For most traits related to ear morphology, we detected multiple QTL per trait in the tb1 genomic region, including a large effect QTL at tb1 itself plus one or two additional linked QTL. tb1 is epistatic to two of these additional QTL for ear traits. Overall, these results provide examples for both a major QTL that maps to a single gene, as well as a case in which a QTL fractionates into multiple linked QTL.     

An ultra-high density bin-map for rapid QTL mapping for tassel and ear architecture in a large F2 maize population

Background

Understanding genetic control of tassel and ear architecture in maize (Zea mays L. ssp. mays) is important due to their relationship with grain yield. High resolution QTL mapping is critical for understanding the underlying molecular basis of phenotypic variation. Advanced populations, such as recombinant inbred lines, have been broadly adopted for QTL mapping; however, construction of large advanced generation crop populations is time-consuming and costly. The rapidly declining cost of genotyping due to recent advances in next-generation sequencing technologies has generated new possibilities for QTL mapping using large early generation populations.

Results

A set of 708 F₂ progeny derived from inbreds Chang7-2 and 787 were generated and genotyped by whole genome low-coverage genotyping-by-sequencing method (average 0.04×). A genetic map containing 6,533 bin-markers was constructed based on the parental SNPs and a sliding-window method, spanning a total genetic distance of 1,396 cM. The high quality and accuracy of this map was validated by the identification of two well-studied genes, r1, a qualitative trait locus for color of silk (chromosome 10) and ba1 for tassel branch number (chromosome 3). Three traits of tassel and ear architecture were evaluated in this population, a total of 10 QTL were detected using a permutation-based-significance threshold, seven of which overlapped with reported QTL. Three genes (GRMZM2G316366, GRMZM2G492156 and GRMZM5G805008) encoding MADS-box domain proteins and a BTB/POZ domain protein were located in the small intervals of qTBN5 and qTBN7 (~800 Kb and 1.6 Mb in length, respectively) and may be involved in patterning of tassel architecture. The small physical intervals of most QTL indicate high-resolution mapping is obtainable with this method.

Conclusions

We constructed an ultra-high-dentisy linkage map for the large early generation population in maize. Our study provides an efficient approach for fast detection of quantitative loci responsible for complex trait variation with high accuracy, thus helping to dissect the underlying molecular basis of phenotypic variation and accelerate improvement of crop breeding in a cost-effective fashion.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-433) contains supplementary material, which is available to authorized users.     

Mapping of quantitative trait loci for flesh colour and growth traits in Atlantic salmon (Salmo salar)

Background

Flesh colour and growth related traits in salmonids are both commercially important and of great interest from a physiological and evolutionary perspective. The aim of this study was to identify quantitative trait loci (QTL) affecting flesh colour and growth related traits in an F2 population derived from an isolated, landlocked wild population in Norway (Byglands Bleke) and a commercial production population.

Methods

One hundred and twenty-eight informative microsatellite loci distributed across all 29 linkage groups in Atlantic salmon were genotyped in individuals from four F2 families that were selected from the ends of the flesh colour distribution. Genotyping of 23 additional loci and two additional families was performed on a number of linkage groups harbouring putative QTL. QTL analysis was performed using a line-cross model assuming fixation of alternate QTL alleles and a half-sib model with no assumptions about the number and frequency of QTL alleles in the founder populations.

Results

A moderate to strong phenotypic correlation was found between colour, length and weight traits. In total, 13 genome-wide significant QTL were detected for all traits using the line-cross model, including three genome-wide significant QTL for flesh colour (Chr 6, Chr 26 and Chr 4). In addition, 32 suggestive QTL were detected (chromosome-wide P < 0.05). Using the half-sib model, six genome-wide significant QTL were detected for all traits, including two for flesh colour (Chr 26 and Chr 4) and 41 suggestive QTL were detected (chromosome-wide P < 0.05). Based on the half-sib analysis, these two genome-wide significant QTL for flesh colour explained 24% of the phenotypic variance for this trait.

Conclusions

A large number of significant and suggestive QTL for flesh colour and growth traits were found in an F2 population of Atlantic salmon. Chr 26 and Chr 4 presented the strongest evidence for significant QTL affecting flesh colour, while Chr 10, Chr 5, and Chr 4 presented the strongest evidence for significant QTL affecting growth traits (length and weight). These QTL could be strong candidates for use in marker-assisted selection and provide a starting point for further characterisation of the genetic components underlying flesh colour and growth.     

Genetic variation in the pleiotropic association between physical activity and body weight in mice

Background

A sedentary lifestyle is often assumed to lead to increases in body weight and potentially obesity and related diseases but in fact little is known about the genetic association between physical activity and body weight. We tested for such an association between body weight and the distance, duration, and speed voluntarily run by 310 mice from the F₂ generation produced from an intercross of two inbred lines that differed dramatically in their physical activity levels.

Methods

We used a conventional interval mapping approach with SNP markers to search for QTLs that affected both body weight and activity traits. We also conducted a genome scan to search for relationship QTLs (relQTLs), or chromosomal regions that affected an activity trait variably depending on the phenotypic value of body weight.

Results

We uncovered seven quantitative trait loci (QTLs) affecting body weight, but only one co-localized with another QTL previously found for activity traits. We discovered 19 relQTLs that provided evidence for a genetic (pleiotropic) association of physical activity and body weight. The three genotypes at each of these loci typically exhibited a combination of negative, zero, and positive regressions of the activity traits on body weight, the net effect of which was to produce overall independence of body weight from physical activity. We also demonstrated that the relQTLs produced these varying associations through differential epistatic interactions with a number of other epistatic QTLs throughout the genome.

Conclusion

It was concluded that individuals with specific combinations of genotypes at the relQTLs and epiQTLs might account for some of the variation typically seen in plots of the association of physical activity with body weight.     

Combining two Meishan F2 crosses improves the detection of QTL on pig chromosomes 2, 4 and 6

Background

In pig, a number of experiments have been set up to identify QTL and a multitude of chromosomal regions harbouring genes influencing traits of interest have been identified. However, the mapping resolution remains limited in most cases and the detected QTL are rather inaccurately located. Mapping accuracy can be improved by increasing the number of phenotyped and genotyped individuals and/or the number of informative markers. An alternative approach to overcome the limited power of individual studies is to combine data from two or more independent designs.

Methods

In the present study we report a combined analysis of two independent design (a French and a Dutch F2 experimental designs), with 2000 F2 individuals. The purpose was to further map QTL for growth and fatness on pig chromosomes 2, 4 and 6. Using QTL-map software, uni- and multiple-QTL detection analyses were applied separately on the two pedigrees and then on the combination of the two pedigrees.

Results

Joint analyses of the combined pedigree provided (1) greater significance of shared QTL, (2) exclusion of false suggestive QTL and (3) greater mapping precision for shared QTL.

Conclusions

Combining two Meishan x European breeds F2 pedigrees improved the mapping of QTL compared to analysing pedigrees separately. Our work was facilitated by the access to raw phenotypic data and DNA of animals from both pedigrees and the combination of the two designs with the addition of new markers allowed us to fine map QTL without phenotyping additional animals.     

A genome-wide association study of production traits in a commercial population of Large White pigs: evidence of haplotypes affecting meat quality

Background

Numerous quantitative trait loci (QTL) have been detected in pigs over the past 20 years using microsatellite markers. However, due to the low density of these markers, the accuracy of QTL location has generally been poor. Since 2009, the dense genome coverage provided by the Illumina PorcineSNP60 BeadChip has made it possible to more accurately map QTL using genome-wide association studies (GWAS). Our objective was to perform high-density GWAS in order to identify genomic regions and corresponding haplotypes associated with production traits in a French Large White population of pigs.

Methods

Animals (385 Large White pigs from 106 sires) were genotyped using the PorcineSNP60 BeadChip and evaluated for 19 traits related to feed intake, growth, carcass composition and meat quality. Of the 64 432 SNPs on the chip, 44 412 were used for GWAS with an animal mixed model that included a regression coefficient for the tested SNPs and a genomic kinship matrix. SNP haplotype effects in QTL regions were then tested for association with phenotypes following phase reconstruction based on the Sscrofa10.2 pig genome assembly.

Results

Twenty-three QTL regions were identified on autosomes and their effects ranged from 0.25 to 0.75 phenotypic standard deviation units for feed intake and feed efficiency (four QTL), carcass (12 QTL) and meat quality traits (seven QTL). The 10 most significant QTL regions had effects on carcass (chromosomes 7, 10, 16, 17 and 18) and meat quality traits (two regions on chromosome 1 and one region on chromosomes 8, 9 and 13). Thirteen of the 23 QTL regions had not been previously described. A haplotype block of 183 kb on chromosome 1 (six SNPs) was identified and displayed three distinct haplotypes with significant (0.0001 < P < 0.03) associations with all evaluated meat quality traits.

Conclusions

GWAS analyses with the PorcineSNP60 BeadChip enabled the detection of 23 QTL regions that affect feed consumption, carcass and meat quality traits in a LW population, of which 13 were novel QTL. The proportionally larger number of QTL found for meat quality traits suggests a specific opportunity for improving these traits in the pig by genomic selection.     

Whole genome association study identifies regions of the bovine genome and biological pathways involved in carcass trait performance in Holstein-Friesian cattle

Background

Four traits related to carcass performance have been identified as economically important in beef production: carcass weight, carcass fat, carcass conformation of progeny and cull cow carcass weight. Although Holstein-Friesian cattle are primarily utilized for milk production, they are also an important source of meat for beef production and export. Because of this, there is great interest in understanding the underlying genomic structure influencing these traits. Several genome-wide association studies have identified regions of the bovine genome associated with growth or carcass traits, however, little is known about the mechanisms or underlying biological pathways involved. This study aims to detect regions of the bovine genome associated with carcass performance traits (employing a panel of 54,001 SNPs) using measures of genetic merit (as predicted transmitting abilities) for 5,705 Irish Holstein-Friesian animals. Candidate genes and biological pathways were then identified for each trait under investigation.

Results

Following adjustment for false discovery (q-value < 0.05), 479 quantitative trait loci (QTL) were associated with at least one of the four carcass traits using a single SNP regression approach. Using a Bayesian approach, 46 QTL were associated (posterior probability > 0.5) with at least one of the four traits. In total, 557 unique bovine genes, which mapped to 426 human orthologs, were within 500kbs of QTL found associated with a trait using the Bayesian approach. Using this information, 24 significantly over-represented pathways were identified across all traits. The most significantly over-represented biological pathway was the peroxisome proliferator-activated receptor (PPAR) signaling pathway.

Conclusions

A large number of genomic regions putatively associated with bovine carcass traits were detected using two different statistical approaches. Notably, several significant associations were detected in close proximity to genes with a known role in animal growth such as glucagon and leptin. Several biological pathways, including PPAR signaling, were shown to be involved in various aspects of bovine carcass performance. These core genes and biological processes may form the foundation for further investigation to identify causative mutations involved in each trait. Results reported here support previous findings suggesting conservation of key biological processes involved in growth and metabolism.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-837) contains supplementary material, which is available to authorized users.     

Major regulatory genes in maize contribute to standing variation in teosinte (Zea mays ssp. parviglumis)

In plants, many major regulatory genes that control plant growth and development have been identified and characterized. Despite a detailed knowledge of the function of these genes little is known about how they contribute to the natural variation for complex traits. To determine whether major regulatory genes of maize contribute to standing variation in Balsas teosinte we conducted association mapping in 584 Balsas teosinte individuals. We tested 48 markers from nine candidate regulatory genes against 13 traits for plant and inflorescence architecture. We identified significant associations using a mixed linear model that controls for multiple levels of relatedness. Ten associations involving five candidate genes were significant after correction for multiple testing, and two survive the conservative Bonferroni correction. zfl2, the maize homolog of FLORICAULA of Antirrhinum, was associated with plant height. zap1, the maize homolog of APETALA1 of Arabidopsis, was associated with inflorescence branching. Five SNPs in the maize domestication gene, teosinte branched1, were significantly associated with either plant or inflorescence architecture. Our data suggest that major regulatory genes in maize do play a role in the natural variation for complex traits in teosinte and that some of the minor variants we identified may have been targets of selection during domestication.