Meta-analysis of Correlated Traits via Summary Statistics from GWASs with an Application in Hypertension

Genome-wide association studies (GWASs) have identified many genetic variants underlying complex traits. Many detected genetic loci harbor variants that associate with multiple—even distinct—traits. Most current analysis approaches focus on single traits, even though the final results from multiple traits are evaluated together. Such approaches miss the opportunity to systemically integrate the phenome-wide data available for genetic association analysis. In this study, we propose a general approach that can integrate association evidence from summary statistics of multiple traits, either correlated, independent, continuous, or binary traits, which might come from the same or different studies. We allow for trait heterogeneity effects. Population structure and cryptic relatedness can also be controlled. Our simulations suggest that the proposed method has improved statistical power over single-trait analysis in most of the cases we studied. We applied our method to the Continental Origins and Genetic Epidemiology Network (COGENT) African ancestry samples for three blood pressure traits and identified four loci (CHIC2, HOXA-EVX1, IGFBP1/IGFBP3, and CDH17; p < 5.0 × 10⁻⁸) associated with hypertension-related traits that were missed by a single-trait analysis in the original report. Six additional loci with suggestive association evidence (p < 5.0 × 10⁻⁷) were also observed, including CACNA1D and WNT3. Our study strongly suggests that analyzing multiple phenotypes can improve statistical power and that such analysis can be executed with the summary statistics from GWASs. Our method also provides a way to study a cross phenotype (CP) association by using summary statistics from GWASs of multiple phenotypes.     

Harvesting the Promising Fruits of Genomics: Applying Genome Sequencing Technologies to Crop Breeding

Next generation sequencing (NGS) technologies are being used to generate whole genome sequences for a wide range of crop species. When combined with precise phenotyping methods, these technologies provide a powerful and rapid tool for identifying the genetic basis of agriculturally important traits and for predicting the breeding value of individuals in a plant breeding population. Here we summarize current trends and future prospects for utilizing NGS-based technologies to develop crops with improved trait performance and increase the efficiency of modern plant breeding. It is our hope that the application of NGS technologies to plant breeding will help us to meet the challenge of feeding a growing world population.

This article is part of the PLOS Biology Collection “The Promise of Plant Translational Research.”

     

Classical phenotyping and deep learning concur on genetic control of stomatal density and area in sorghum

Stomatal density (SD) and stomatal complex area (SCA) are important traits that regulate gas exchange and abiotic stress response in plants. Despite sorghum (Sorghum bicolor) adaptation to arid conditions, the genetic potential of stomata-related traits remains unexplored due to challenges in available phenotyping methods. Hence, identifying loci that control stomatal traits is fundamental to designing strategies to breed sorghum with optimized stomatal regulation. We implemented both classical and deep learning methods to characterize genetic diversity in 311 grain sorghum accessions for stomatal traits at two different field environments. Nearly 12,000 images collected from abaxial (Ab) and adaxial (Ad) leaf surfaces revealed substantial variation in stomatal traits. Our study demonstrated significant accuracy between manual and deep learning methods in predicting SD and SCA. In sorghum, SD was 32%–39% greater on the Ab versus the Ad surface, while SCA on the Ab surface was 2%–5% smaller than on the Ad surface. Genome-Wide Association Study identified 71 genetic loci (38 were environment-specific) with significant genotype to phenotype associations for stomatal traits. Putative causal genes underlying the phenotypic variation were identified. Accessions with similar SCA but carrying contrasting haplotypes for SD were tested for stomatal conductance and carbon assimilation under field conditions. Our findings provide a foundation for further studies on the genetic and molecular mechanisms controlling stomata patterning and regulation in sorghum. An integrated physiological, deep learning, and genomic approach allowed us to unravel the genetic control of natural variation in stomata traits in sorghum, which can be applied to other plants.

High-throughput phenotyping using deep learning tools integrated with genome-wide association studies revealed genes that control SD and area in grain sorghum.     

Distinct genetic architectures for male and female inflorescence traits of maize

We compared the genetic architecture of thirteen maize morphological traits in a large population of recombinant inbred lines. Four traits from the male inflorescence (tassel) and three traits from the female inflorescence (ear) were measured and studied using linkage and genome-wide association analyses and compared to three flowering and three leaf traits previously studied in the same population. Inflorescence loci have larger effects than flowering and leaf loci, and ear effects are larger than tassel effects. Ear trait models also have lower predictive ability than tassel, flowering, or leaf trait models. Pleiotropic loci were identified that control elongation of ear and tassel, consistent with their common developmental origin. For these pleiotropic loci, the ear effects are larger than tassel effects even though the same causal polymorphisms are likely involved. This implies that the observed differences in genetic architecture are not due to distinct features of the underlying polymorphisms. Our results support the hypothesis that genetic architecture is a function of trait stability over evolutionary time, since the traits that changed most during the relatively recent domestication of maize have the largest effects.     

Optical topometry and machine learning to rapidly phenotype stomatal patterning traits for maize QTL mapping

Stomata are adjustable pores on leaf surfaces that regulate the tradeoff of CO₂ uptake with water vapor loss, thus having critical roles in controlling photosynthetic carbon gain and plant water use. The lack of easy, rapid methods for phenotyping epidermal cell traits have limited discoveries about the genetic basis of stomatal patterning. A high-throughput epidermal cell phenotyping pipeline is presented here and used for quantitative trait loci (QTL) mapping in field-grown maize (Zea mays). The locations and sizes of stomatal complexes and pavement cells on images acquired by an optical topometer from mature leaves were automatically determined. Computer estimated stomatal complex density (SCD; R² = 0.97) and stomatal complex area (SCA; R² = 0.71) were strongly correlated with human measurements. Leaf gas exchange traits were genetically correlated with the dimensions and proportions of stomatal complexes (r_g = 0.39–0.71) but did not correlate with SCD. Heritability of epidermal traits was moderate to high (h² = 0.42–0.82) across two field seasons. Thirty-six QTL were consistently identified for a given trait in both years. Twenty-four clusters of overlapping QTL for multiple traits were identified, with univariate versus multivariate single marker analysis providing evidence consistent with pleiotropy in multiple cases. Putative orthologs of genes known to regulate stomatal patterning in Arabidopsis (Arabidopsis thaliana) were located within some, but not all, of these regions. This study demonstrates how discovery of the genetic basis for stomatal patterning can be accelerated in maize, a C₄ model species where these processes are poorly understood.

Optical topometry and machine learning tools allow assessment of epidermal cell patterning and analysis of its genetic architecture alongside leaf photosynthetic gas exchange in maize.     

Plasticity Regulators Modulate Specific Root Traits in Discrete Nitrogen Environments

Plant development is remarkably plastic but how precisely can the plant customize its form to specific environments? When the plant adjusts its development to different environments, related traits can change in a coordinated fashion, such that two traits co-vary across many genotypes. Alternatively, traits can vary independently, such that a change in one trait has little predictive value for the change in a second trait. To characterize such “tunability” in developmental plasticity, we carried out a detailed phenotypic characterization of complex root traits among 96 accessions of the model Arabidopsis thaliana in two nitrogen environments. The results revealed a surprising level of independence in the control of traits to environment – a highly tunable form of plasticity. We mapped genetic architecture of plasticity using genome-wide association studies and further used gene expression analysis to narrow down gene candidates in mapped regions. Mutants in genes implicated by association and expression analysis showed precise defects in the predicted traits in the predicted environment, corroborating the independent control of plasticity traits. The overall results suggest that there is a pool of genetic variability in plants that controls traits in specific environments, with opportunity to tune crop plants to a given environment.     

Four linked genes participate in controlling sporulation efficiency in budding yeast

Quantitative traits are conditioned by several genetic determinants. Since such genes influence many important complex traits in various organisms, the identification of quantitative trait loci (QTLs) is of major interest, but still encounters serious difficulties. We detected four linked genes within one QTL, which participate in controlling sporulation efficiency in Saccharomyces cerevisiae. Following the identification of single nucleotide polymorphisms by comparing the sequences of 145 genes between the parental strains SK1 and S288c, we analyzed the segregating progeny of the cross between them. Through reciprocal hemizygosity analysis, four genes, RAS2, PMS1, SWS2, and FKH2, located in a region of 60 kilobases on Chromosome 14, were found to be associated with sporulation efficiency. Three of the four “high” sporulation alleles are derived from the “low” sporulating strain. Two of these sporulation-related genes were verified through allele replacements. For RAS2, the causative variation was suggested to be a single nucleotide difference in the upstream region of the gene. This quantitative trait nucleotide accounts for sporulation variability among a set of ten closely related winery yeast strains. Our results provide a detailed view of genetic complexity in one “QTL region” that controls a quantitative trait and reports a single nucleotide polymorphism-trait association in wild strains. Moreover, these findings have implications on QTL identification in higher eukaryotes.     

The Influence of Genetic Drift and Selection on Quantitative Traits in a Plant Pathogenic Fungus

Genetic drift and selection are ubiquitous evolutionary forces acting to shape genetic variation in populations. While their relative importance has been well studied in plants and animals, less is known about their relative importance in fungal pathogens. Because agro-ecosystems are more homogeneous environments than natural ecosystems, stabilizing selection may play a stronger role than genetic drift or diversifying selection in shaping genetic variation among populations of fungal pathogens in agro-ecosystems. We tested this hypothesis by conducting a Q _ST/F _ST analysis using agricultural populations of the barley pathogen Rhynchosporium commune. Population divergence for eight quantitative traits (Q _ST) was compared with divergence at eight neutral microsatellite loci (F _ST) for 126 pathogen strains originating from nine globally distributed field populations to infer the effects of genetic drift and types of selection acting on each trait. Our analyses indicated that five of the eight traits had Q _ST values significantly lower than F _ST, consistent with stabilizing selection, whereas one trait, growth under heat stress (22°C), showed evidence of diversifying selection and local adaptation (Q _ST>F _ST). Estimates of heritability were high for all traits (means ranging between 0.55–0.84), and average heritability across traits was negatively correlated with microsatellite gene diversity. Some trait pairs were genetically correlated and there was significant evidence for a trade-off between spore size and spore number, and between melanization and growth under benign temperature. Our findings indicate that many ecologically and agriculturally important traits are under stabilizing selection in R. commune and that high within-population genetic variation is maintained for these traits.     

High-Throughput Phenotyping to Detect Drought Tolerance QTL in Wild Barley Introgression Lines

Drought is one of the most severe stresses, endangering crop yields worldwide. In order to select drought tolerant genotypes, access to exotic germplasm and efficient phenotyping protocols are needed. In this study the high-throughput phenotyping platform “The Plant Accelerator”, Adelaide, Australia, was used to screen a set of 47 juvenile (six week old) wild barley introgression lines (S42ILs) for drought stress responses. The kinetics of growth development was evaluated under early drought stress and well watered treatments. High correlation (r = 0.98) between image based biomass estimates and actual biomass was demonstrated, and the suitability of the system to accurately and non-destructively estimate biomass was validated. Subsequently, quantitative trait loci (QTL) were located, which contributed to the genetic control of growth under drought stress. In total, 44 QTL for eleven out of 14 investigated traits were mapped, which for example controlled growth rate and water use efficiency. The correspondence of those QTL with QTL previously identified in field trials is shown. For instance, six out of eight QTL controlling plant height were also found in previous field and glasshouse studies with the same introgression lines. This indicates that phenotyping juvenile plants may assist in predicting adult plant performance. In addition, favorable wild barley alleles for growth and biomass parameters were detected, for instance, a QTL that increased biomass by approximately 36%. In particular, introgression line S42IL-121 revealed improved growth under drought stress compared to the control Scarlett. The introgression line showed a similar behavior in previous field experiments, indicating that S42IL-121 may be an attractive donor for breeding of drought tolerant barley cultivars.     

Evaluation of Linkage Disequilibrium Pattern and Association Study on Seed Oil Content in Brassica napus Using ddRAD Sequencing

High-density genetic markers are the prerequisite for understanding linkage disequilibrium (LD) and genome-wide association studies (GWASs) of complex traits in crops. To evaluate the LD pattern in oilseed rape, we sequenced a previous association panel containing 189 B. napus inbred lines using double-digested restriction-site associated DNA (ddRAD) and genotyped 19,327 RAD tags. A total of 15,921 RAD tags were assigned to a published genetic linkage map and the majority (71.1%) of these tags was uniquely mapped to the draft reference genome “Darmor-bzh.” The distance of LD decay was 1,214 kb across the genome at the background level (r² = 0.26), with the distances of LD decay being 405 kb and 2,111 kb in the A and C subgenomes, respectively. A total of 361 haplotype blocks with length > 100 kb were identified in the entire genome. The association panel could be classified into two groups, P1 and P2, which are essentially consistent with the geographical origins of varieties. A large number of group-specific haplotypes were identified, reflecting that varieties in the P1 and P2 groups experienced distinct selection in breeding programs to adapt their different growth habitats. GWAS repeatedly detected two loci significantly associated with oil content of seeds based on the developed SNPs, suggesting that the high-density SNPs were useful for understanding the genetic determinants of complex traits in GWAS.     

Multiple modes of selection can influence the role of phenotypic plasticity in species' invasions: Evidence from a manipulative field experiment

In exploring the roles of phenotypic plasticity in the establishment and early evolution of invading species, little empirical attention has been given to the importance of correlational selection acting upon suites of functionally related plastic traits in nature. We illustrate how this lack of attention has limited our ability to evaluate plasticity''s role during invasion and also, the costs and benefits of plasticity. We addressed these issues by transplanting clones of European‐derived Plantago lanceolata L. genotypes into two temporally variable habitats in the species'' introduced range in North America. Phenotypic selection analyses were performed for each habitat to estimate linear, quadratic, and correlational selection on phenotypic trait values and plasticities in the reproductive traits: flowering onset and spike and scape lengths. Also, we measured pairwise genetic correlations for our “colonists.” Results showed that (a) correlational selection acted on trait plasticity after transplantation, (b) selection favored certain combinations of genetically correlated and uncorrelated trait values and plasticities, and (c) using signed, instead of absolute, values of plasticity in analyses facilitated the detection of correlational selection on trait value‐plasticity combinations and their adaptive value. Based on our results, we urge future studies on species invasions to (a) measure correlational selection and (b) retain signed values of plasticity in order to better discriminate between adaptive and maladaptive plasticity.     

High-throughput root phenotyping screens identify genetic loci associated with root architectural traits in Brassica napus under contrasting phosphate availabilities

Background and Aims

Phosphate (Pi) deficiency in soils is a major limiting factor for crop growth worldwide. Plant growth under low Pi conditions correlates with root architectural traits and it may therefore be possible to select these traits for crop improvement. The aim of this study was to characterize root architectural traits, and to test quantitative trait loci (QTL) associated with these traits, under low Pi (LP) and high Pi (HP) availability in Brassica napus.

Methods

Root architectural traits were characterized in seedlings of a double haploid (DH) mapping population (n = 190) of B. napus [‘Tapidor’ × ‘Ningyou 7’ (TNDH)] using high-throughput phenotyping methods. Primary root length (PRL), lateral root length (LRL), lateral root number (LRN), lateral root density (LRD) and biomass traits were measured 12 d post-germination in agar at LP and HP.

Key Results

In general, root and biomass traits were highly correlated under LP and HP conditions. ‘Ningyou 7’ had greater LRL, LRN and LRD than ‘Tapidor’, at both LP and HP availability, but smaller PRL. A cluster of highly significant QTL for LRN, LRD and biomass traits at LP availability were identified on chromosome A03; QTL for PRL were identified on chromosomes A07 and C06.

Conclusions

High-throughput phenotyping of Brassica can be used to identify root architectural traits which correlate with shoot biomass. It is feasible that these traits could be used in crop improvement strategies. The identification of QTL linked to root traits under LP and HP conditions provides further insights on the genetic basis of plant tolerance to P deficiency, and these QTL warrant further dissection.     

Impact of Measurement Error on Testing Genetic Association with Quantitative Traits

Measurement error of a phenotypic trait reduces the power to detect genetic associations. We examined the impact of sample size, allele frequency and effect size in presence of measurement error for quantitative traits. The statistical power to detect genetic association with phenotype mean and variability was investigated analytically. The non-centrality parameter for a non-central F distribution was derived and verified using computer simulations. We obtained equivalent formulas for the cost of phenotype measurement error. Effects of differences in measurements were examined in a genome-wide association study (GWAS) of two grading scales for cataract and a replication study of genetic variants influencing blood pressure. The mean absolute difference between the analytic power and simulation power for comparison of phenotypic means and variances was less than 0.005, and the absolute difference did not exceed 0.02. To maintain the same power, a one standard deviation (SD) in measurement error of a standard normal distributed trait required a one-fold increase in sample size for comparison of means, and a three-fold increase in sample size for comparison of variances. GWAS results revealed almost no overlap in the significant SNPs (p<10⁻⁵) for the two cataract grading scales while replication results in genetic variants of blood pressure displayed no significant differences between averaged blood pressure measurements and single blood pressure measurements. We have developed a framework for researchers to quantify power in the presence of measurement error, which will be applicable to studies of phenotypes in which the measurement is highly variable.     

RootAnalyzer: A Cross-Section Image Analysis Tool for Automated Characterization of Root Cells and Tissues

The morphology of plant root anatomical features is a key factor in effective water and nutrient uptake. Existing techniques for phenotyping root anatomical traits are often based on manual or semi-automatic segmentation and annotation of microscopic images of root cross sections. In this article, we propose a fully automated tool, hereinafter referred to as RootAnalyzer, for efficiently extracting and analyzing anatomical traits from root-cross section images. Using a range of image processing techniques such as local thresholding and nearest neighbor identification, RootAnalyzer segments the plant root from the image’s background, classifies and characterizes the cortex, stele, endodermis and epidermis, and subsequently produces statistics about the morphological properties of the root cells and tissues. We use RootAnalyzer to analyze 15 images of wheat plants and one maize plant image and evaluate its performance against manually-obtained ground truth data. The comparison shows that RootAnalyzer can fully characterize most root tissue regions with over 90% accuracy.     

Genetic Diversity and Association Studies in US Hispanic/Latino Populations: Applications in the Hispanic Community Health Study/Study of Latinos

US Hispanic/Latino individuals are diverse in genetic ancestry, culture, and environmental exposures. Here, we characterized and controlled for this diversity in genome-wide association studies (GWASs) for the Hispanic Community Health Study/Study of Latinos (HCHS/SOL). We simultaneously estimated population-structure principal components (PCs) robust to familial relatedness and pairwise kinship coefficients (KCs) robust to population structure, admixture, and Hardy-Weinberg departures. The PCs revealed substantial genetic differentiation within and among six self-identified background groups (Cuban, Dominican, Puerto Rican, Mexican, and Central and South American). To control for variation among groups, we developed a multi-dimensional clustering method to define a “genetic-analysis group” variable that retains many properties of self-identified background while achieving substantially greater genetic homogeneity within groups and including participants with non-specific self-identification. In GWASs of 22 biomedical traits, we used a linear mixed model (LMM) including pairwise empirical KCs to account for familial relatedness, PCs for ancestry, and genetic-analysis groups for additional group-associated effects. Including the genetic-analysis group as a covariate accounted for significant trait variation in 8 of 22 traits, even after we fit 20 PCs. Additionally, genetic-analysis groups had significant heterogeneity of residual variance for 20 of 22 traits, and modeling this heteroscedasticity within the LMM reduced genomic inflation for 19 traits. Furthermore, fitting an LMM that utilized a genetic-analysis group rather than a self-identified background group achieved higher power to detect previously reported associations. We expect that the methods applied here will be useful in other studies with multiple ethnic groups, admixture, and relatedness.     

Quercus species divergence is driven by natural selection on evolutionarily less integrated traits

Functional traits are organismal attributes that can respond to environmental cues, thereby providing important ecological functions. In addition, an organism’s potential for adaptation is defined by the patterns of covariation among groups of functionally related traits. Whether an organism is evolutionarily constrained or has the potential for adaptation is based on the phenotypic integration or modularity of these traits. Here, we revisited leaf morphology in two European sympatric white oaks (Quercus petraea (Matt.) Liebl. and Quercus robur L.), sampling 2098 individuals, across much of their geographical distribution ranges. At the phenotypic level, leaf morphology traditionally encompasses discriminant attributes among different oak species. Here, we estimated in situ heritability, genetic correlation, and integration across such attributes. Also, we performed Selection Response Decomposition to test these traits for potential differences in oak species’ evolutionary responses. Based on the uncovered functional units of traits (modules) in our study, the morphological module “leaf size gradient” was highlighted among functionally integrated traits. Equally, this module was defined in both oaks as being under “global regulation” in vegetative bud establishment and development. Lamina basal shape and intercalary veins’ number were not, or, less integrated within the initially defined leaf functional unit, suggesting more than one module within the leaf traits’ ensemble. Since these traits generally show the greatest species discriminatory power, they potentially underwent effective differential response to selection among oaks. Indeed, the selection of these traits could have driven the ecological preferences between the two sympatric oaks growing under different microclimates.Subject terms: Quantitative trait, Evolutionary biology     

Genetic and genomic analysis of a fat mass trait with complex inheritance reveals marked sex specificity

The integration of expression profiling with linkage analysis has increasingly been used to identify genes underlying complex phenotypes. The effects of gender on the regulation of many physiological traits are well documented; however, “genetical genomic” analyses have not yet addressed the degree to which their conclusions are affected by sex. We constructed and densely genotyped a large F2 intercross derived from the inbred mouse strains C57BL/6J and C3H/HeJ on an apolipoprotein E null (ApoE^−/−) background. This BXH.ApoE^−/− population recapitulates several “metabolic syndrome” phenotypes. The cross consists of 334 animals of both sexes, allowing us to specifically test for the dependence of linkage on sex. We detected several thousand liver gene expression quantitative trait loci, a significant proportion of which are sex-biased. We used these analyses to dissect the genetics of gonadal fat mass, a complex trait with sex-specific regulation. We present evidence for a remarkably high degree of sex-dependence on both the cis and trans regulation of gene expression. We demonstrate how these analyses can be applied to the study of the genetics underlying gonadal fat mass, a complex trait showing significantly female-biased heritability. These data have implications on the potential effects of sex on the genetic regulation of other complex traits.     

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.     

Morphological traits defining species differences in wild relatives of maize are controlled by multiple quantitative trait loci

We analyzed the genetic basis of morphological differences between two wild species of teosinte (Zea diploperennis and Z. mays ssp. parviglumis), which are relatives of maize. These two species differ in a number of taxonomically important traits including the structure of the tassel (male inflorescence), which is the focus of this report. To investigate the genetic inheritance of six tassel traits, quantitative trait locus (QTL) mapping with 95 RFLP markers was employed on a population of 425 F2 plants. Each trait was analyzed by interval mapping (IM) and composite interval mapping (CIM) to identify and characterize the QTL controlling the differences in tassel morphology. We detected two to eight QTL for each trait. In total, 30 QTL with IM and 33 QTL with CIM were found for tassel morphology. QTL for several of the traits mapped near each other, suggesting pleiotropy and/or linkage of QTL. The QTL showed small to moderate magnitudes of effect. No QTL of exceptionally large effect were found as seen under domestication and in the case of some other natural species. Thus, the model involving major QTL of large effect seems not to apply to the traits and species analyzed. A mixture of QTL with positive and negative allelic effects was found for most tassel traits and may suggest a history of periodic changes in the direction of selection during the divergence of Z. diploperennis and Z. mays ssp. parviglumis or fixation of QTL alleles by random genetic drift rather than selection.