np2QTL: networking phenotypic plasticity quantitative trait loci across heterogeneous environments

Despite its critical importance to our understanding of plant growth and adaptation, the question of how environment‐induced plastic response is affected genetically remains elusive. Previous studies have shown that the reaction norm of an organism across environmental index obeys the allometrical scaling law of part‐whole relationships. The implementation of this phenomenon into functional mapping can characterize how quantitative trait loci (QTLs) modulate the phenotypic plasticity of complex traits to heterogeneous environments. Here, we assemble functional mapping and allometry theory through Lokta?Volterra ordinary differential equations (LVODE) into an R‐based computing platform, np²QTL, aimed to map and visualize phenotypic plasticity QTLs. Based on LVODE parameters, np²QTL constructs a bidirectional, signed and weighted network of QTL?QTL epistasis, whose emergent properties reflect the ecological mechanisms for genotype?environment interactions over any range of environmental change. The utility of np²QTL was validated by comprehending the genetic architecture of phenotypic plasticity via the reanalysis of published plant height data involving 3502 recombinant inbred lines of maize planted in multiple discrete environments. np²QTL also provides a tool for constructing a predictive model of phenotypic responses in extreme environments relative to the median environment.     

Mendelian Factors Underlying Quantitative Traits in Tomato: Comparison across Species, Generations, and Environments

As part of ongoing studies regarding the genetic basis of quantitative variation in phenotype, we have determined the chromosomal locations of quantitative trait loci (QTLs) affecting fruit size, soluble solids concentration, and pH, in a cross between the domestic tomato (Lycopersicon esculentum Mill.) and a closely-related wild species, L. cheesmanii. Using a RFLP map of the tomato genome, we compared the inheritance patterns of polymorphisms in 350 F2 individuals with phenotypes scored in three different ways: (1) from the F2 progeny themselves, grown near Davis, California; (2) from F3 families obtained by selfing each F2 individual, grown near Gilroy, California (F3-CA); and (3) from equivalent F3 families grown near Rehovot, Israel (F3-IS). Maximum likelihood methods were used to estimate the approximate chromosomal locations, phenotypic effects (both additive effects and dominance deviations), and gene action of QTLs underlying phenotypic variation in each of these three environments. A total of 29 putative QTLs were detected in the three environments. These QTLs were distributed over 11 of the 12 chromosomes, accounted for 4.7-42.0% of the phenotypic variance in a trait, and showed different types of gene action. Among these 29 QTLs, 4 were detected in all three environments, 10 in two environments, and 15 in only a single environment. The two California environments were most similar, sharing 11/25 (44%) QTLs, while the Israel environment was quite different, sharing 7/20 (35%) and 5/26 (19%) QTLs with the respective California environments. One major goal of QTL mapping is to predict, with maximum accuracy, which individuals will produce progeny showing particular phenotypes. Traditionally, the phenotype of an individual alone has been used to predict the phenotype of its progeny. Our results suggested that, for a trait with low heritability (soluble solids), the phenotype of F3 progeny could be predicted more accurately from the genotype of the F2 parent at QTLs than from the phenotype of the F2 individual. For a trait with intermediate heritability (fruit pH), QTL genotype and observed phenotype were about equally effective at predicting progeny phenotype. For a trait with high heritability (mass per fruit), knowing the QTL genotype of an individual added little if any predictive value, to simply knowing the phenotype. The QTLs mapped in the L. esculentum X L. cheesmanii F2 appear to be at similar locations to many of those mapped in a previous cross with a different wild tomato (L. chmielewskii).(ABSTRACT TRUNCATED AT 400 WORDS)     

A functional-structural model of rice linking quantitative genetic information with morphological development and physiological processes

Background and Aims

Although quantitative trait loci (QTL) analysis of yield-related traits for rice has developed rapidly, crop models using genotype information have been proposed only relatively recently. As a first step towards a generic genotype–phenotype model, we present here a three-dimensional functional–structural plant model (FSPM) of rice, in which some model parameters are controlled by functions describing the effect of main-effect and epistatic QTLs.

Methods

The model simulates the growth and development of rice based on selected ecophysiological processes, such as photosynthesis (source process) and organ formation, growth and extension (sink processes). It was devised using GroIMP, an interactive modelling platform based on the Relational Growth Grammar formalism (RGG). RGG rules describe the course of organ initiation and extension resulting in final morphology. The link between the phenotype (as represented by the simulated rice plant) and the QTL genotype was implemented via a data interface between the rice FSPM and the QTLNetwork software, which computes predictions of QTLs from map data and measured trait data.

Key Results

Using plant height and grain yield, it is shown how QTL information for a given trait can be used in an FSPM, computing and visualizing the phenotypes of different lines of a mapping population. Furthermore, we demonstrate how modification of a particular trait feeds back on the entire plant phenotype via the physiological processes considered.

Conclusions

We linked a rice FSPM to a quantitative genetic model, thereby employing QTL information to refine model parameters and visualizing the dynamics of development of the entire phenotype as a result of ecophysiological processes, including the trait(s) for which genetic information is available. Possibilities for further extension of the model, for example for the purposes of ideotype breeding, are discussed.Key words: Functional–structural plant model, ecophysiology, QTL analysis, plant modelling, quantitative genetics     

A multiparental cross population for mapping QTL for agronomic traits in durum wheat (Triticum turgidum ssp. durum)

Multiparental cross designs for mapping quantitative trait loci (QTL) provide an efficient alternative to biparental populations because of their broader genetic basis and potentially higher mapping resolution. We describe the development and deployment of a recombinant inbred line (RIL) population in durum wheat (Triticum turgidum ssp. durum) obtained by crossing four elite cultivars. A linkage map spanning 2664 cM and including 7594 single nucleotide polymorphisms (SNPs) was produced by genotyping 338 RILs. QTL analysis was carried out by both interval mapping on founder haplotype probabilities and SNP bi‐allelic tests for heading date and maturity date, plant height and grain yield from four field experiments. Sixteen QTL were identified across environments and detection methods, including two yield QTL on chromosomes 2BL and 7AS, with the former mapped independently from the photoperiod response gene Ppd‐B1, while the latter overlapped with the vernalization locus VRN‐A3. Additionally, 21 QTL with environment‐specific effects were found. Our results indicated a prevalence of environment‐specific QTL with relatively small effect on the control of grain yield. For all traits, functionally different QTL alleles in terms of direction and size of genetic effect were distributed among parents. We showed that QTL results based on founder haplotypes closely matched functional alleles at known heading date loci. Despite the four founders, only 2.1 different functional haplotypes were estimated per QTL, on average. This durum wheat population provides a mapping resource for detailed genetic dissection of agronomic traits in an elite background typical of breeding programmes.     

Extreme QTL mapping of germination speed in Arabidopsis thaliana

Seed germination is a key life history transition for annual plants and partly determines lifetime performance and fitness. Germination speed, the elapsed time for a nondormant seed to germinate, is a poorly understood trait important for plants’ competitiveness and fitness in fluctuating environments. Germination speed varied by 30% among 18 Arabidopsis thaliana populations measured, and exhibited weak negative correlation with flowering time and seed weight, with significant genotype effect (P < 0.005). To dissect the genetic architecture of germination speed, we developed the extreme QTL (X‐QTL) mapping method in A. thaliana. The method has been shown in yeast to increase QTL mapping power by integrating selective screening and bulk‐segregant analysis in a very large mapping population. By pooled genotyping of top 5% of rapid germinants from ~100 000 F₃ individuals, three X‐QTL regions were identified on chromosomes 1, 3 and 4. All regions were confirmed as QTL regions by sequencing 192 rapid germinants from an independent F₃ selection experiment. Positional overlaps were found between X‐QTLs and previously identified seed, life history and fitness QTLs. Our method provides a rapid mapping platform in A. thaliana with potentially greater power. One can also relate identified X‐QTLs to the A. thaliana physical map, facilitating candidate gene identification.     

A unified DNA sequence and non‐DNA sequence mapping model of complex traits

Increasing evidence shows that quantitative inheritance is based on both DNA sequence and non‐DNA sequence variants. However, how to simultaneously detect these variants from a mapping study has been unexplored, hampering our effort to illustrate the detailed genetic architecture of complex traits. We address this issue by developing a unified model of quantitative trait locus (QTL) mapping based on an open‐pollinated design composed of randomly sampling maternal plants from a natural population and their half‐sib seeds. This design forms a two‐level hierarchical platform for a joint linkage‐linkage disequilibrium analysis of population structure. The EM algorithm was implemented to estimate and test DNA sequence‐based effects and non‐DNA sequence‐based effects of QTLs. We applied this model to analyze genetic mapping data from the OP design of a gymnosperm coniferous species, Torreya grandis, identifying 25 significant DNA sequence and non‐DNA sequence QTLs for seedling height and diameter growth in different years. Results from computer simulation show that the unified model has good statistical properties and is powerful for QTL detection. Our model enables the tests of how a complex trait is affected differently by DNA‐based effects and non‐DNA sequence‐based transgenerational effects, thus allowing a more comprehensive picture of genetic architecture to be charted and quantified.     

Combining quantitative trait Loci analysis and an ecophysiological model to analyze the genetic variability of the responses of maize leaf growth to temperature and water deficit

Ecophysiological models predict quantitative traits of one genotype in any environment, whereas quantitative trait locus (QTL) models predict the contribution of alleles to quantitative traits under a limited number of environments. We have combined both approaches by dissecting into effects of QTLs the parameters of a model of maize (Zea mays) leaf elongation rate (LER; H. Ben Haj Salah, F. Tardieu [1997] Plant Physiol 114: 893-900). Response curves of LER to meristem temperature, water vapor pressure difference, and soil water status were established in 100 recombinant inbred lines (RILs) of maize in six experiments carried out in the field or in the greenhouse. All responses were linear and common to different experiments, consistent with the model. A QTL analysis was carried out on the slopes of these responses by composite interval mapping confirmed by bootstrap analysis. Most QTLs were specific of one response only. QTLs of abscisic acid concentration in the xylem sap colocalized with QTLs of response to soil water deficit and conferred a low response. Each parameter of the ecophysiological model was computed as the sum of QTL effects, allowing calculation of parameters for 11 new RILs and two parental lines. LERs were simulated and compared with measurements in a growth chamber experiment. The combined model accounted for 74% of the variability of LER, suggesting that it has a general value for any RIL under any environment.     

Mapping of quantitative trait loci based on growth models

An approach called growth model-based mapping (GMM) of quantitative trait loci (QTLs) is proposed in this paper. The principle of the approach is to fit the growth curve of each individual or line with a theoretical or empirical growth model at first and then map QTLs based on the estimated growth parameters with the method of multiple-trait composite interval mapping. In comparison with previously proposed approaches of QTL mapping based on growth data, GMM has several advantages: (1) it can greatly reduce the amount of phenotypic data for QTL analysis and thus alleviate the burden of computation, particularly when permutation tests or simulation are performed to estimate significance thresholds; (2) it can efficiently analyze unbalanced phenotype data because both balanced and unbalanced data can be used for fitting growth models; and (3) it may potentially help us to better understand the genetic basis of quantitative trait development because the parameters in a theoretical growth model may often have clear biological meanings. A practical example of rice leaf-age development is presented to demonstrate the utility of GMM.     

Computational identification of genes modulating stem height–diameter allometry

The developmental variation in stem height with respect to stem diameter is related to a broad range of ecological and evolutionary phenomena in trees, but the underlying genetic basis of this variation remains elusive. We implement a dynamic statistical model, functional mapping, to formulate a general procedure for the computational identification of quantitative trait loci (QTLs) that control stem height–diameter allometry during development. Functional mapping integrates the biological principles underlying trait formation and development into the association analysis of DNA genotype and endpoint phenotype, thus providing an incentive for understanding the mechanistic interplay between genes and development. Built on the basic tenet of functional mapping, we explore two core ecological scenarios of how stem height and stem diameter covary in response to environmental stimuli: (i) trees pioneer sunlit space by allocating more growth to stem height than diameter and (ii) trees maintain their competitive advantage through an inverse pattern. The model is equipped to characterize ‘pioneering’ QTLs (piQTLs) and ‘maintaining’ QTLs (miQTLs) which modulate these two ecological scenarios, respectively. In a practical application to a mapping population of full‐sib hybrids derived from two Populus species, the model has well proven its versatility by identifying several piQTLs that promote height growth at a cost of diameter growth and several miQTLs that benefit radial growth at a cost of height growth. Judicious application of functional mapping may lead to improved strategies for studying the genetic control of the formation mechanisms underlying trade‐offs among quantities of assimilates allocated to different growth parts.     

QTL and QTL x environment effects on agronomic and nitrogen acquisition traits in rice

Agricultural environments deteriorate due to excess nitrogen application.Breeding for low nitrogen responsive genotypes can reduce soil nitrogen input.Rice genotypes respond variably to soil available nitrogen.The present study attempted quantification of genotype x nitrogen level interaction and mapping of quantitative trait loci (QTLs) associated with nitrogen use efficiency (NUE) and other associated agronomic traits.Twelve parameters were observed across a set of 82 double haploid (DH) lines derived from IR64/Azucena.Three nitrogen regimes namely,native (0 kg/ha; no nitrogen applied),optimum (100 kg/ha) and high (200 kg/ha) replicated thrice were the environments.The parents and DH lines were significantly varying for all traits under different nitrogen regimes.All traits except plant height recorded significant genotype x environment interaction.Individual plant yield was positively correlated with nitrogen use efficiency and nitrogen uptake.Sixteen QTLs were detected by composite interval mapping.Eleven QTLs showed significant QTL x environment interactions.On chromosome 3,seven QTLs were detected associated with nitrogen use,plant yield and associated traits.A QTL region between markers RZ678,RZ574 and RZ284 was associated with nitrogen use and yield.This chromosomal region was enriched with expressed gene sequences of known key nitrogen assimilation genes.     

A method for identification of the expression mode and mapping of QTL underlying embryo-specific characters

Embryos of crop seeds are one of the major sources of the plant protein and lipid for human nutrition. The genetic expression for embryo-specific characters in crop seeds can be controlled exclusively by the embryo or the maternal genotypes and sometimes by both simultaneously. However, current methods for mapping quantitative trait loci (QTLs) underlying characters of maternal plants have not been effective in dealing with the QTL analysis of embryo characters. On the basis of the expression feature of embryo, a statistical method was proposed for the identification of expression mode and mapping of QTL controlling embryo traits. The maximum likelihood method implemented via the expectation maximization algorithm was used to estimate parameters of a putative embryo-specific QTL. The QTL expression mode was identified by the likelihood ratio test statistic. Statistical power and other properties of the proposed method were investigated under a variety of scenarios through simulation studies. The results showed that the mapping method neglecting the effects of embryo genotype or maternal effects could neither identify the expression mode of QTL nor estimate its genetic effects accurately, whereas the proposed method could effectively map the embryo-specific QTL of various expression modes.     

Quantitative trait loci for growth trajectories in Populus

Growth trajectories are a biological process important to plant and animal breeding, and to evolutionary genetic studies. In this article, we report the detection of quantitative trait loci (QTLs) responsible for growth trajectories in poplars that are used as a model system for the study of forest biology. These QTLs were localized on a genetic linkage map of polymorphic markers using a statistical mapping method incorporating growth-curve models. The effects of the QTLs on growth are described as a function of age, so that age-specific changes in QTL effects can be readily projected throughout the entire growth process. The QTLs identified display increased effects on growth when trees age, yet the timing of QTL activation is earlier for stem height than diameter, which is consistent with the ecological viewpoint of canopy competition. The implications of the results for breeding and silviculture are discussed.     

Assessing the importance of genotype × environment interaction for root traits in rice using a mapping population II: conventional QTL analysis

Modifying plant root systems is considered a means of crop improvement targeted to low-resource environments, particularly low nutrient and drought-prone agriculture. The identification of quantitative trait loci (QTLs) for root traits has stimulated marker-assisted breeding to this end, but different QTLs have been detected in different populations of the same species, and importantly, in the same population when grown in different experimental environments. The presence of QTL × environment interaction is implicated, and this must be characterised if the utility of the target QTLs is to be realised. Previous attempts to do this suffer from a lack of control over replicate environments and inadequate statistical rigour. The Bala × Azucena mapping population was grown in two replicate experiments of four treatment environments, a control, a low light, a low soil nitrogen and a low soil water treatment. After a 4 weeks growth, maximum root length, maximum root thickness, root mass below 50 cm, total plant dry mass, % root mass and shoot length were measured. A summary of the overall results is presented in an accompanying paper. Here, QTL analysis by composite interval mapping is presented. A total of 145 QTLs were detected, mapping to 37 discrete loci on all chromosomes. Superficial evidence of QTL × E (great difference in LOD score) was tested by single-marker analysis which confirmed QTL × E for five loci representing only five individual trait-loci interactions. Some loci appeared to be stable across environments. Some QTLs were clearly more or less active under low light, low nitrogen or drought. A few notable loci on chromosomes 1, 2, 3, 5, 7 and 9 are briefly discussed. Also discussed are some remaining statistical shortcomings that will be addressed in another companion paper.     

Two reproductive traits show contrasting genetic architectures in Plantago lanceolata

In many species, temperature‐sensitive phenotypic plasticity (i.e., an individual's phenotypic response to temperature) displays a positive correlation with latitude, a pattern presumed to reflect local adaptation. This geographical pattern raises two general questions: (a) Do a few large‐effect genes contribute to latitudinal variation in a trait? (b) Is the thermal plasticity of different traits regulated pleiotropically? To address the questions, we crossed individuals of Plantago lanceolata derived from northern and southern European populations. Individuals naturally exhibited high and low thermal plasticity in floral reflectance and flowering time. We grew parents and offspring in controlled cool‐ and warm‐temperature environments, mimicking what plants would encounter in nature. We obtained genetic markers via genotype‐by‐sequencing, produced the first recombination map for this ecologically important nonmodel species, and performed quantitative trait locus (QTL) mapping of thermal plasticity and single‐environment values for both traits. We identified a large‐effect QTL that largely explained the reflectance plasticity differences between northern and southern populations. We identified multiple smaller‐effect QTLs affecting aspects of flowering time, one of which affected flowering time plasticity. The results indicate that the genetic architecture of thermal plasticity in flowering is more complex than for reflectance. One flowering time QTL showed strong cytonuclear interactions under cool temperatures. Reflectance and flowering plasticity QTLs did not colocalize, suggesting little pleiotropic genetic control and freedom for independent trait evolution. Such genetic information about the architecture of plasticity is environmentally important because it informs us about the potential for plasticity to offset negative effects of climate change.     

Genetic analysis of resistance to yellow rust in hexaploid wheat using a mixture model for multiple crosses

DNA-based molecular markers have been used in numerous studies for tagging specific genes in wheat for subsequent use in marker-assisted selection. Usually in plant breeding, procedures for mapping genes are based on analysis of a single segregating population. However, breeding programmes routinely evaluate large numbers of progeny derived from multiple-related crosses with some parental lines shared. In most such related crosses, the number of progeny is quite small. Thus, statistical techniques for detecting quantitative trait loci (QTLs) using data from conventional multi-cross breeding programmes are interesting. The objective of this study is to present a mixture model for QTL mapping in crosses of multiple inbred varieties with non-normal phenotype distributions and to use this model to map QTLs for yellow rust resistance in elite wheat breeding material. Three doubled haploid populations consisting of 41, 42 and 55 lines, respectively, originating from four parental varieties were studied. Multi-cross QTL analysis with three specific pathogen isolates of Puccinia striiformis f. sp. tritici and a mixture of the isolates revealed QTLs for resistance at four different genomic locations. These QTLs were found on chromosome 2AL, 2AS, 2BL and 6BL and explained between 21 and 41% of the phenotypic variation. Two of these QTLs, one on the long arm of chromosome 2A and one on the short arm of chromosome 2A were identical to the known yellow rust resistance genes Yr32 and Yr17, respectively, whereas the QTLs located on the long arms of chromosomes 2B and 6B may reflect types of resistance to yellow rust, which have not previously been mapped.     

Potential of a tomato MAGIC population to decipher the genetic control of quantitative traits and detect causal variants in the resequencing era

Identification of the polymorphisms controlling quantitative traits remains a challenge for plant geneticists. Multiparent advanced generation intercross (MAGIC) populations offer an alternative to traditional linkage or association mapping populations by increasing the precision of quantitative trait loci (QTL) mapping. Here, we present the first tomato MAGIC population and highlight its potential for the valorization of intraspecific variation, QTL mapping and causal polymorphism identification. The population was developed by crossing eight founder lines, selected to include a wide range of genetic diversity, whose genomes have been previously resequenced. We selected 1536 SNPs among the 4 million available to enhance haplotype prediction and recombination detection in the population. The linkage map obtained showed an 87% increase in recombination frequencies compared to biparental populations. The prediction of the haplotype origin was possible for 89% of the MAGIC line genomes, allowing QTL detection at the haplotype level. We grew the population in two greenhouse trials and detected QTLs for fruit weight. We mapped three stable QTLs and six specific of a location. Finally, we showed the potential of the MAGIC population when coupled with whole genome sequencing of founder lines to detect candidate SNPs underlying the QTLs. For a previously cloned QTL on chromosome 3, we used the predicted allelic effect of each founder and their genome sequences to select putative causal polymorphisms in the supporting interval. The number of candidate polymorphisms was reduced from 12 284 (in 800 genes) to 96 (in 54 genes), including the actual causal polymorphism. This population represents a new permanent resource for the tomato genetics community.     

QTL sequencing strategy to map genomic regions associated with resistance to ascochyta blight in chickpea

Whole‐genome sequencing‐based bulked segregant analysis (BSA) for mapping quantitative trait loci (QTL) provides an efficient alternative approach to conventional QTL analysis as it significantly reduces the scale and cost of analysis with comparable power to QTL detection using full mapping population. We tested the application of next‐generation sequencing (NGS)‐based BSA approach for mapping QTLs for ascochyta blight resistance in chickpea using two recombinant inbred line populations CPR‐01 and CPR‐02. Eleven QTLs in CPR‐01 and six QTLs in CPR‐02 populations were mapped on chromosomes Ca1, Ca2, Ca4, Ca6 and Ca7. The QTLs identified in CPR‐01 using conventional biparental mapping approach were used to compare the efficiency of NGS‐based BSA in detecting QTLs for ascochyta blight resistance. The QTLs on chromosomes Ca1, Ca4, Ca6 and Ca7 overlapped with the QTLs previously detected in CPR‐01 using conventional QTL mapping method. The QTLs on chromosome Ca4 were detected in both populations and overlapped with the previously reported QTLs indicating conserved region for ascochyta blight resistance across different chickpea genotypes. Six candidate genes in the QTL regions identified using NGS‐based BSA on chromosomes Ca2 and Ca4 were validated for their association with ascochyta blight resistance in the CPR‐02 population. This study demonstrated the efficiency of NGS‐based BSA as a rapid and cost‐effective method to identify QTLs associated with ascochyta blight in chickpea.     

Combined genetic and physiological analysis of a locus contributing to quantitative variation

The natural variation of many traits is controlled by multiple genes, individually referred to as quantitative trait loci (QTL), that interact with the environment to determine the ultimate phenotype of any individual. A QTL has yet to be described molecularly, in part because strategies to systematically identify them are underdeveloped and because the subtle nature of QTLs prevents the application of standard methods of gene identification. Therefore, it will be necessary to develop a systematic approach(es) for the identification of QTLs based upon the numerous positional data now being accumulated through molecular marker analyses. We have characterized a QTL by the following three-step approach: (1) identification of a QTL in complex populations, (2) isolation and genetic mapping of this QTL in near-isogenic lines, and (3) identification of a candidate gene using map position and physiological criteria. Using this approach we have characterized a plant height QTL in maize that maps to chromosome 9 near the centromere. Both map position and physiological criteria suggest the gibberillin biosynthesis gene dwarf3 as a candidate gene for this QTL.     

QTL mapping analysis of plant height and ear height of maize (Zea mays L.)

Genetic map containing 103 microsatellite loci obtained on 200 F2 plants derived from the cross R15 x 478 was used for quantitative trait loci (QTL) mapping in maize. QTL were characterized in a population of 200 F2:4 lines, derived from selfing the F2 plants, and were evaluated with two replications in two environments. QTL determinations were made from the mean of these two environments. Plant height (PH) and ear height (EH) were measured. Using composite interval mapping (CIM) method, a total of 14 distinct QTLs were identified: nine for PH and five for EH. Additive, partial dominance, dominance, and overdominance actions existed among all detected QTL affecting plant height and ear height. The QTL explained 78.27% of the phenotypic variance of PH and 41.50% of EH. The 14 QTLs displayed mostly dominance or partial dominance gene action and mapped to chromosomes 2, 3, 4, 8 and 9.     

Identification of quantitative trait loci for agronomically important traits and their association with genic-microsatellite markers in sorghum

The identification of quantitative trait loci (QTLs) affecting agronomically important traits enable to understand their underlying genetic mechanisms and genetic basis of their complex interactions. The aim of the present study was to detect QTLs for 12 agronomic traits related to staygreen, plant early development, grain yield and its components, and some growth characters by analyzing replicated phenotypic datasets from three crop seasons, using the population of 168 F₇ RILs of the cross 296B × IS18551. In addition, we report mapping of a subset of genic-microsatellite markers. A linkage map was constructed with 152 marker loci comprising 149 microsatellites (100 genomic- and 49 genic-microsatellites) and three morphological markers. QTL analysis was performed by using MQM approach. Forty-nine QTLs were detected, across environments or in individual environments, with 1–9 QTLs for each trait. Individual QTL accounted for 5.2–50.4% of phenotypic variance. Several genomic regions affected multiple traits, suggesting the phenomenon of pleiotropy or tight linkage. Stable QTLs were identified for studied traits across different environments, and genetic backgrounds by comparing the QTLs in the study with previously reported QTLs in sorghum. Of the 49 mapped genic-markers, 18 were detected associating either closely or exactly as the QTL positions of agronomic traits. EST marker Dsenhsbm19, coding for a key regulator (EIL-1) of ethylene biosynthesis, was identified co-located with the QTLs for plant early development and staygreen trait, a probable candidate gene for these traits. Similarly, such exact co-locations between EST markers and QTLs were observed in four other instances. Collectively, the QTLs/markers identified in the study are likely candidates for improving the sorghum performance through MAS and map-based gene isolations.