BACKGROUND
drought is a severe abiotic stress that affects the production and productivity of rice.
drought stress at the reproductive stage is the most devastating  <cit> .
because of the ongoing process of climate change, the rainfall pattern has become more irregular in the cropping season, causing widespread drought in rice-growing areas, which results in severe yield losses  <cit> .
the development of drought-tolerant varieties that maintain good yield under drought is a priority area of rice research for sustainable rice production.
marker-assisted mapping and introgression of major-effect qtl for grain yield under drought could be an efficient and fast-track approach for breeding drought-tolerant rice varieties  <cit> .
however, the successful use of qtl in marker-assisted selection depends on their effect and consistency across genetic backgrounds and environments.
most of the qtl for grain yield under drought have been mapped against a single genetic background in early-segregating generations  evaluated in a limited number of environments.
such qtl may not provide a consistent effect because of variation in the genetic background and environment.
additionally, the qtl may not be transferrable to other backgrounds because of unfavorable epistatic interactions resulting in reduced or even no effects in a new genetic background  <cit> .
considering all these facts, it is difficult to predict the usefulness of qtl for mas based only on their performance in an individual genetic background in any particular study.
a more efficient way to select qtl for mas is to compare the identified qtl with earlier reported studies for their consistency of location and effect across genetic backgrounds and environments.
consistently identified qtl at the same chromosomal location, explaining high phenotypic variance and having a major effect on a trait, can be effectively used in mas  <cit> .
qtl meta-analysis is an approach to identify consensus qtl across studies, to validate qtl effects across environments/genetic backgrounds, and also to refine qtl positions on the consensus map  <cit> .
qtl meta-analysis requires independent qtl for the same trait obtained from different populations, different locations, or different environmental conditions  <cit> .
the consistent qtl identified by meta-analysis for a set of qtl at a confidence interval  of 95% are called meta-qtl .
the meta-qtl with the smallest ci and having a consistent and large effect on a trait are useful in mas.
in plants, the concept of meta-analysis has been applied to the analysis of qtl/genes for blast resistance  <cit> , root traits and drought tolerance in rice  <cit> , lint fiber length in cotton  <cit> , cyst nematode resistance in soybean  <cit> , fusarium head blight resistance in wheat  <cit> , flowering time  <cit> , drought tolerance in maize  <cit> , and disease resistance in cocoa  <cit> .
qtl validation is another approach to confirm the effect of qtl across different genetic backgrounds.
qtl regions harbor many genes; among them, a few key genes could be more important in the regulation of a complex trait.
meta-qtl regions with refined positions are more accurate for short-listing of candidate genes.
the common candidate genes short-listed across the meta-qtl are more likely candidates that regulate yield  <cit> .
in this study, qtl meta-analysis was carried out for yield qtl under drought to develop a consensus map and to identify consensus yield qtl under drought with the objective to provide markers of mqtl with high effects and small confidence intervals for possible use in mas or for fine-mapping qtl for gene discovery.
also, markers linked to  <dig> major qtl for grain yield were validated on a set of random drought-tolerant lines, including landraces and improved drought breeding lines developed at irri, to know the frequency of their universal presence.
further, a comparative genomics approach was used to identify the homologous regions of mqtl in other cereal crops such as maize, sorghum, wheat, and barley .
materials and methods
meta-qtl analysis
three steps were employed for the identification of a consensus qtl for grain yield under drought.
first, in a bibliographic review, reliable data on qtl for yield per plant were compiled.
second, a consensus map was created and on this map the qtl of individual studies was projected.
in the third step, a meta-analysis was performed on qtl clusters to identify the consensus mqtl.
bibliographic review and synthesis of yield qtl data
qtl information was collected from published reports involving mapping of qtl for grain yield under drought.
there were  <dig> reports of a qtl mapping for grain yield under drought.
the details of the parents used in developing the mapping population, size of the mapping population, markers used, and yield qtl identified are given in table  <dig>  in all,  <dig> qtl were reported for yield.
bsa = bulk segregant analysis; wg = whole genotyping; aflp = amplified fragment length polymorphism; rflp =restricted fragment length polymorphism, ssr = simple sequence repeats, bac = bacterial artificial chromosome
development of a consensus map
a consensus genetic map was constructed and meta-analysis was performed using biomercator v <dig>  .
the rice genetic linkage map of temnykh et al.
<cit>  was used as a reference map, on which the markers of  <dig> studies were projected to develop a consensus map.
chromosomes connected with fewer than two common markers to the reference map were excluded before the creation of the consensus map.
inversions of marker sequences were filtered out by discarding inconsistent loci with the exception of very closely linked markers.
after the integration of all maps, the consensus map contained  <dig> markers, including ssr, rflp, aflp markers, and genes.
the consensus map covered a total length of  <dig> cm, with an average distance of  <dig>  cm between markers.
qtl projections
for all studies, the 95% confidence intervals of initial qtl on their original maps were estimated using the approach described by darvasi and soller  <cit> :
where n is the population size and r <dig> the proportion of the phenotypic variance explained by the qtl.
the ci was re-estimated to control the heterogeneity of ci calculation methods across studies.
projection of qtl positions was performed by using a simple scaling rule between the original qtl flanking marker interval and the corresponding interval on the consensus chromosome.
for a given qtl position, the new ci on the consensus linkage group was approximated with a gaussian distribution around the most likely qtl position.
all projections of qtl onto the consensus map were performed using the biomercator  .
meta-analysis
meta-analysis was performed on the qtl clusters on each chromosome using biomercator  .
the akaike information criterion  was used to select the qtl model on each chromosome  <cit> .
according to this, the qtl model with the lowest aic value is considered a significant model indicating the number of meta-qtl.
qtl meta-analysis requires independent qtl for the same trait obtained from different plant populations, different locations, or different environmental conditions  <cit> .
qtl validation
genotyping
all molecular marker work was conducted in the gene array and molecular marker analysis  laboratory, plant breeding, genetics and biotechnology  division, irri.
for dna extraction, freeze-dried samples were used.
freeze-dried leaf samples were cut in eppendorf tubes and ground through a geno grinder.
extraction was carried out by the modified ctab method.
dna samples were stored in 2-ml deep-well plates .
dna samples were quantified on  <dig> % agarose gel and concentration adjusted to approximately  <dig> ng μl- <dig>  pcr amplification was done with a 15-μl reaction mixture having  <dig> ng dna,  <dig> × pcr buffer,  <dig> μm dntps,  <dig> μm primers, and  <dig> unit taq polymerase enzyme.
the pcr profiles started with an initial denaturation of dna at 94°c for  <dig> minutes, followed by  <dig> amplification cycles of denaturation at 94°c for  <dig> minute, annealing temperatures varied from 55°c to 58°c for  <dig> seconds based on the primer, extension at 72°c for  <dig> minute and final extension at 72°c for  <dig> minutes.
the pcr products were resolved on 8% non-denaturing polyacrilamide gels .
the gels were scored taking respective qtl donor alleles as reference band and scores were used for qtl validation.
the details of the peak markers of the  <dig> major effect qtl are given in additional file  <dig>
twelve major effect drought grain yield qtl were validated on a panel of  <dig> drought tolerant lines consisting of traditional drought tolerant donors, drought tolerant breeding lines developed through conventional breeding approaches and random high yielding lines under drought from qtl mapping populations.
the peak marker of all the twelve major effect qtl were amplified on the drought panel lines.
the lines were scored taking qtl donor allele as a base.
the list of lines is given in the additional file  <dig>
gene content analysis
the  <dig> meta-qtl were analyzed for gene content to know the presence of genes and gene clusters responsible for drought.
a comparative genomics approach was followed to analyze the genes present in meta-qtl.
gene content was noted based on annotated data of homologous regions in nipponbare using rap, build <dig> .
it is assumed that the genes identified in nipponbare regions are homologous and collinear to those underlying the yield qtl under drought mapped in different studies involving different donors and recipients.
comparative genomics to identify homologous regions in cereals
a comparative genomics approach was followed to identify homologous regions between rice and maize using the genomic databases .
homologous regions identified were checked for the presence of drought grain yield qtl of maize .
in sorghum, wheat, and barley, grain yield qtl reported were collected from a literature survey and these were compared with the meta-qtl using the comparative maps available in the gramene database .
RESULTS
overview of qtl and development of a consensus map
in the  <dig> populations of rice screened for drought tolerance to map qtl, population size ranged from  <dig>  <cit>  to  <dig> lines  <cit> .
the number of markers used ranged from  <dig> to  <dig>  <cit> .
the number of locations for phenotyping varied from  <dig> to  <dig>  from the  <dig> studies,  <dig> yield qtl were reported, which were distributed on all the chromosomes except chromosome  <dig> .
the number of qtl per population ranged from  <dig> to  <dig>  the proportion of qtl per chromosome ranged from one qtl each on chromosomes  <dig> and  <dig> to  <dig> yield qtl on chromosome  <dig>  the distribution of yield qtl on different chromosomes showed that chromosomes  <dig>   <dig>  and  <dig> have the highest number,  <dig>   <dig>  and  <dig> qtl, respectively .
the phenotypic variance of the initial qtl varied from  <dig> % to 40% and the confidence interval of the markers varied from  <dig> to  <dig> cm.
the rice genetic map of temnykh et al.
<cit>  was used as a reference map to develop a consensus map as this is a widely used genetic map of rice and it contained most of the markers used in the different studies.
the consensus map consisted of  <dig> markers with a total map length of  <dig> cm.
the average distance between the markers was  <dig>  cm, thus enabling the identification of a precise location of qtl.
there were very few marker inversions in the consensus map, which were discarded from the final map and further analysis.
meta-analysis and qtl validation
it is widely believed that qtl are accurate and can be positioned onto chromosomal locations by molecular mapping  <cit> .
however, their complex nature and context dependency in different genetic backgrounds and environments are constraints in identifying their precise location.
the identification of the most accurate and precise major-effect qtl across genetic backgrounds and environments is a prerequisite for the successful use of qtl in mas.
meta-analysis of qtl identified in different studies helps to identify the most precise and concise qtl, which can be further pursued for mas or the identification of candidate genes.
in our study, we attempted to identify the meta-qtl for grain yield under drought by genome-wide meta-analysis.
from a bibliographic survey, a total of  <dig> qtl were short-listed for grain yield under drought from  <dig> studies.
all  <dig> qtl were projected on a consensus map.
the chromosomal regions with only one qtl were not considered for meta-analysis since meta-analysis by definition involves more than one qtl.
thus,  <dig> qtl were used for meta-analysis and meta-qtl were short-listed based on the akaike information criterion .
accordingly, the qtl model with the lowest aic value was considered a significant model indicating the number of meta-qtl.
the number of meta-qtl along with their aic values and confidence intervals are given in table  <dig>  in total,  <dig> independent meta-qtl were identified at a confidence interval of 95% on seven chromosomes, and meta-analysis successfully reduced the total qtl by 63% .
the meta-qtl identified on each chromosome varied from  <dig> to  <dig>  there were four meta-qtl on chromosome 1; two on chromosomes  <dig>   <dig>   <dig> and 10; and one each on chromosomes  <dig> and  <dig>  the phenotypic variance of the meta-qtl varied from 4% to 28%.
at  <dig> of the  <dig> meta-qtl, the mean phenotypic variance was more than 10%.
in general, the confidence intervals at most of the meta-qtl were narrower than their respective original qtl.
at nine loci on chromosomes  <dig>   <dig>   <dig>   <dig>   <dig>  and  <dig>  meta-qtl were narrower than the mean of their initial qtl.
however, at five loci, the meta-qtl were broader than the mean of the initial qtl.
the confidence intervals of the meta-qtl varied from  <dig>  cm between the marker intervals rg <dig> and rm <dig> on chromosome  <dig> to  <dig>  cm between the marker intervals rm <dig> and rm <dig> on chromosome  <dig>  at two regions, meta-qtl <dig>   and mqtl <dig> , the ci declined to around  <dig> cm.
the physical intervals of the meta-qtl varied from  <dig>  mb to  <dig>  mb.
three meta-qtl were less than  <dig> kb.
the meta-qtl regions with small genetic and physical intervals are useful in mas.
it is significant to see that seven qtl that had less than  <dig>  mb intervals also had a genetic interval of around  <dig> cm with a phenotypic variance of more than 10% .
three of these meta-qtl were on chromosome  <dig> and one each on chromosomes  <dig>   <dig> and  <dig>  the physical intervals of mqtl <dig> , mqtl <dig> , and mqtl <dig>  were less than  <dig> kb, that of mqtl <dig>  was  <dig> kb, and those of mqtl <dig> , mqtl <dig> , and mqtl <dig>  were 1mb,  <dig>  mb, and  <dig>  mb, respectively.
mqtl <dig> , mqtl <dig> , and mqtl <dig>  had phenotypic variance of more than 20%.
the seven mqtl regions with small genetic and physical intervals are important regions for mas, fine mapping, candidate gene identification, and functional analysis.
these qtl can be introgressed in popular rice mega-varieties to develop drought-tolerant and high-yielding lines.
in addition to meta-analysis of qtl, the markers linked to the  <dig> major-effect qtl for grain yield were also validated on a panel of drought-tolerant lines to confirm their presence in larger set of lines.
it is notable that major-effect qtl dty <dig>  was present in 85% of the lines.
dty <dig> , dty <dig>  dty <dig> , dty <dig> , and dty <dig>  were present in more than 50% of the lines .
the amplification of the rm <dig> and rm <dig> peak markers of dty <dig>  and dty <dig>  in a set of  <dig> drought tolerant panel lines is presented in additional file  <dig>  the result indicates the presence of at least one of the major-effect grain yield qtl in the drought panel lines.
in general, the major-effect qtl identified for grain yield under drought have a genetic gain of 10% to 30%, with a yield advantage of around  <dig> to  <dig> kg/ha over recipient parents.
however, considering practical benefit to farmers, the development of drought-tolerant rice varieties with a yield advantage of at least  <dig> ton/ha could be the desired target for rice breeders.
the marker-aided qtl pyramiding of the major-effect mqtl identified in this study can be considered as an option for achieving this target.
aic = akaike information criterion, ci = confidence interval, cm = centimorgan, mqtl = meta-qtl; kb = kilobase, mb = megabase
a comparison was made between the meta-qtl identified in this study with the meta-qtl identified for root traits in two earlier studies  <cit> .
it is very interesting to note that mqtl <dig> , mqtl <dig> , mqtl <dig> , mqtl <dig> , and mqtl <dig>  coincided with qtl clusters for root and leaf morphology traits associated with drought tolerance/avoidance in rice  <cit> .
all the  <dig> independent meta-qtl coincided with at least one meta-qtl identified for root traits under drought  <cit> .
earlier studies on meta-analysis of qtl for root traits  <cit>  and blast resistance in rice  <cit> , fusarium head blight resistance in wheat  <cit> , flowering time in maize  <cit> , nematode resistance in soybean  <cit> , and lint fiber length in cotton  <cit>  identified precise and concise meta-qtl.
meta-qtl were also used to deduce candidate genes and were recommended for mas in some of these studies.
comparative genomics of mqtl
the existence of an evolutionary relationship among the grass families is a well-known fact.
the syntenic relationship can be used to identify the homologous regions among these species, which in turn is useful in defining their role in plant growth, development, and adaptation across species.
we compared meta-qtl regions for synteny in other cereal crops.
the major-effect mqtl <dig>  was also found in maize on chromosome  <dig> near marker msu <dig>  in wheat on chromosome 4b near marker rht-b <dig>  and in barley on chromosome 6h near marker bmac <dig>  while major-effect mqtl <dig>  was also found in maize on chromosome  <dig> near marker umc107a .
all these markers were linked to grain yield under drought in their respective crops.
the largest parts of chromosomes  <dig> and  <dig> of rice have a syntenic relation with chromosomes  <dig> and  <dig> of maize, so their respective homologous qtl were also found on the corresponding chromosomes.
an interesting observation is that, near the sd <dig> locus on chromosome  <dig> of rice, qtl for grain yield under drought were identified most frequently.
sd <dig> is a major locus responsible for semidwarf plant stature in rice and its corresponding locus in wheat is rht-b <dig> on chromosome 4b.
mqtl <dig>  is near the sd <dig> locus and also on its corresponding locus rht-b <dig> in wheat, major qtl for grain yield under drought were detected.
gene content analysis and identification of candidate genes
meta qtl with precise and narrow confidence intervals are useful in short listing the candidate genes.
using the annotated gene information available in the rice database, the genes present in the  <dig> meta-qtl regions were analyzed by comparative genomics approach and candidate genes were shortlisted.
the short-listed candidate genes can be further confirmed by transgenic approaches by loss or gain of function studies.
most of the genes present in the mqtl were genes for hypothetical and expressed proteins, pseudo genes, genes for signal transduction, and transposable elements.
however, there were many annotated genes/gene families that were common across the mqtl regions; these are probable candidate genes for yield under drought.
it was found that three kinds of genes frequently occurred together in these regions.
the genes/gene families were stress-inducible genes, growth and development-related genes, and sugar transport-related genes.
table  <dig> lists the important genes underlying mqtl for grain yield under drought.
in six mqtl with less than a  <dig> mb region, lrr kinase, leucine zipper, cell division-controlling proteins, sugar transport protein-like genes, no apical meristem , pentatricopeptide repeat proteins, cytokinin oxidase, f-box proteins, ap2-domain containing proteins, and zinc-finger transcription factors were present.
the candidacy of these genes for yield and yield traits has already been proved in rice and other crops.
cytochrome p <dig> has a role in bassinosteroid homeostasis and had an influence on leaf angle leading to increased yield in rice  <cit> .
pentatricopeptide repeats are present in the promoter region of rf genes, which restore fertility and also play a role in embryogenesis in arabidopsis  <cit> .
zinc-finger -like proteins are known to be involved in stress tolerance.
zinc-finger protein in rice are induced after different types of stresses, namely, cold, desiccation, salt, submergence, heavy metals, and mechanical injury.
over expression of the zinc-finger gene in transgenic tobacco conferred tolerance of cold, dehydration, and salt stress at the seed germination/seedling stage  <cit> .
f-box proteins play an important role in floral development and stress tolerance.
in addition, f-box proteins appear to serve as the key components of the machinery involved in regulating plant growth and development throughout the plant's life cycle and their expression is influenced by light and abiotic stresses  <cit> .
leucine zippers are a class of transcription factor involved in aba-independent stress tolerance.
over expression of osbzip <dig> in rice triggered clusters of genes regulating stress adaptations  <cit> .
the no apical meristem gene  plays an important role in the growth and development of meristematic tissue.
the root-specific expression of this gene resulted in enhanced root growth and improved drought tolerance in rice  <cit> .
the other important genes that harbored the meta-qtl were the erecta and dreb genes.
erecta is a leucine-rich repeat receptor-like kinase gene known for its influence on inflorescence development, stomatal density, epidermal cell expansion, and mesophyll cell proliferation.
this gene is mainly involved in transpiration efficiency and enhanced drought response  <cit> .
dreb is a well-known transcription factor that is induced by drought and it activates many down stream stress-responsive genes to ultimately improve the drought and chilling tolerance of rice  <cit> .
some of these short-listed genes can be considered as positional candidate genes that determine grain yield under drought.
however, it is also well known that yield and adaptability to stress are complex in nature and highly negatively correlated.
the qtl/genes for these two are often co-located.
even though individual genes have been proved to regulate yield under controlled drought experiments, a well-coordinated response of many genes is essential for drought tolerance under field conditions.
this is evident from the presence of three different groups of gene clusters in most of the meta-qtl regions.
CONCLUSIONS
meta-analysis of grain yield qtl is an effective approach in identifying concise and precise consensus qtl.
the seven meta-qtl identified with small genetic and physical intervals could be useful in mas/pyramiding.
validation of the major-effect qtl confirmed the consistency of the major-effect grain yield qtl under drought in different drought-tolerant panel lines.
the comparative genomics approach to identify the consistency of drought grain yield qtl across species revealed the conservation of some of the loci, indicating their evolutionary significance.
the presence of gene clusters in the meta-qtl indicates that a well-coordinated response of many genes is essential to achieve drought tolerance under field conditions.
authors' contributions
ak conceived the idea of meta- analysis, qtl validation and comparative genomics of grain yield qtl under drought.
bpms compiled and analyzed the data, carried out the qtl validation and comparative genomics.
sd helped in data compilation and analysis.
bpms, pv, sd, hua and ak were responsible for drafting and editing the manuscript.
all authors have read and approved the final manuscript.
supplementary material
additional file 1
details of the markers used for qtl validation.
this file contains the list of major effect qtls for grain yield under drought and peak markers of the qtls.
primer sequence, product size of the markers and annealing temperatures  used for amplifying the markers.
click here for file
additional file 2
drought panel lines for qtl validation.
this table shows the list of drought panel lines and type of the breeding material.
these lines were used for validating the major effect qtls for grain yield under drought.
click here for file
additional file 3
amplification of rm <dig> and rm <dig> peak markers of qtl <dig>  and qtl <dig>  in a set of  <dig> drought tolerant panel lines.
the gel picture shows the amplication of rm <dig> and rm <dig> peak markers of qtl <dig>  and qtl <dig>  in a set of  <dig> drought tolerant panel lines.
click here for file
acknowledgements
financial support to this study was provided by the bill & melinda gates foundation, usa.

 drought is a severe abiotic stress that affects the production and productivity of rice. drought stress at the reproductive stage is the most devastating. because of the ongoing process of climate change, the rainfall pattern has become more irregular in the cropping season, causing widespread drought in rice-growing areas, which results in severe yield losses. the development of drought-tolerant varieties that maintain good yield under drought is a priority area of rice research for sustainable rice production. marker-assisted mapping and introgression of major-effect qtl for grain yield under drought could be an efficient and fast-track approach for breeding drought-tolerant rice varieties. however, the successful use of qtl in marker-assisted selection depends on their effect and consistency across genetic backgrounds and environments. most of the qtl for grain yield under drought have been mapped against a single genetic background in early-segregating generations evaluated in a limited number of environments. such qtl may not provide a consistent effect because of variation in the genetic background and environment. additionally, the qtl may not be transferrable to other backgrounds because of unfavorable epistatic interactions resulting in reduced or even no effects in a new genetic background. considering all these facts, it is difficult to predict the usefulness of qtl for mas based only on their performance in an individual genetic background in any particular study. a more efficient way to select qtl for mas is to compare the identified qtl with earlier reported studies for their consistency of location and effect across genetic backgrounds and environments. consistently identified qtl at the same chromosomal location, explaining high phenotypic variance and having a major effect on a trait, can be effectively used in mas. qtl meta-analysis is an approach to identify consensus qtl across studies, to validate qtl effects across environments/genetic backgrounds, and also to refine qtl positions on the consensus map. qtl meta-analysis requires independent qtl for the same trait obtained from different populations, different locations, or different environmental conditions. the consistent qtl identified by meta-analysis for a set of qtl at a confidence interval of 95% are called meta-qtl. the meta-qtl with the smallest ci and having a consistent and large effect on a trait are useful in mas. in plants, the concept of meta-analysis has been applied to the analysis of qtl/genes for blast resistance, root traits and drought tolerance in rice, lint fiber length in cotton, cyst nematode resistance in soybean, fusarium head blight resistance in wheat, flowering time, drought tolerance in maize, and disease resistance in cocoa. qtl validation is another approach to confirm the effect of qtl across different genetic backgrounds. qtl regions harbor many genes; among them, a few key genes could be more important in the regulation of a complex trait. meta-qtl regions with refined positions are more accurate for short-listing of candidate genes. the common candidate genes short-listed across the meta-qtl are more likely candidates that regulate yield. in this study, qtl meta-analysis was carried out for yield qtl under drought to develop a consensus map and to identify consensus yield qtl under drought with the objective to provide markers of mqtl with high effects and small confidence intervals for possible use in mas or for fine-mapping qtl for gene discovery. also, markers linked to major qtl for grain yield were validated on a set of random drought-tolerant lines, including landraces and improved drought breeding lines developed at irri, to know the frequency of their universal presence. further, a comparative genomics approach was used to identify the homologous regions of mqtl in other cereal crops such as maize, sorghum, wheat, and barley. materials and methods meta-qtl analysis three steps were employed for the identification of a consensus qtl for grain yield under drought. first, in a bibliographic review, reliable data on qtl for yield per plant were compiled. second, a consensus map was created and on this map the qtl of individual studies was projected. in the third step, a meta-analysis was performed on qtl clusters to identify the consensus mqtl. bibliographic review and synthesis of yield qtl data qtl information was collected from published reports involving mapping of qtl for grain yield under drought. there were reports of a qtl mapping for grain yield under drought. the details of the parents used in developing the mapping population, size of the mapping population, markers used, and yield qtl identified are given in table in all, qtl were reported for yield. bsa = bulk segregant analysis; wg = whole genotyping; aflp = amplified fragment length polymorphism; rflp =restricted fragment length polymorphism, ssr = simple sequence repeats, bac = bacterial artificial chromosome development of a consensus map a consensus genetic map was constructed and meta-analysis was performed using biomercator v. the rice genetic linkage map of temnykh et al. was used as a reference map, on which the markers of studies were projected to develop a consensus map. chromosomes connected with fewer than two common markers to the reference map were excluded before the creation of the consensus map. inversions of marker sequences were filtered out by discarding inconsistent loci with the exception of very closely linked markers. after the integration of all maps, the consensus map contained markers, including ssr, rflp, aflp markers, and genes. the consensus map covered a total length of cm, with an average distance of cm between markers. qtl projections for all studies, the 95% confidence intervals of initial qtl on their original maps were estimated using the approach described by darvasi and soller: where n is the population size and r the proportion of the phenotypic variance explained by the qtl. the ci was re-estimated to control the heterogeneity of ci calculation methods across studies. projection of qtl positions was performed by using a simple scaling rule between the original qtl flanking marker interval and the corresponding interval on the consensus chromosome. for a given qtl position, the new ci on the consensus linkage group was approximated with a gaussian distribution around the most likely qtl position. all projections of qtl onto the consensus map were performed using the biomercator. meta-analysis meta-analysis was performed on the qtl clusters on each chromosome using biomercator. the akaike information criterion was used to select the qtl model on each chromosome. according to this, the qtl model with the lowest aic value is considered a significant model indicating the number of meta-qtl. qtl meta-analysis requires independent qtl for the same trait obtained from different plant populations, different locations, or different environmental conditions. qtl validation genotyping all molecular marker work was conducted in the gene array and molecular marker analysis laboratory, plant breeding, genetics and biotechnology division, irri. for dna extraction, freeze-dried samples were used. freeze-dried leaf samples were cut in eppendorf tubes and ground through a geno grinder. extraction was carried out by the modified ctab method. dna samples were stored in 2-ml deep-well plates. dna samples were quantified on % agarose gel and concentration adjusted to approximately ng μl- pcr amplification was done with a 15-μl reaction mixture having ng dna, × pcr buffer, μm dntps, μm primers, and unit taq polymerase enzyme. the pcr profiles started with an initial denaturation of dna at 94°c for minutes, followed by amplification cycles of denaturation at 94°c for minute, annealing temperatures varied from 55°c to 58°c for seconds based on the primer, extension at 72°c for minute and final extension at 72°c for minutes. the pcr products were resolved on 8% non-denaturing polyacrilamide gels. the gels were scored taking respective qtl donor alleles as reference band and scores were used for qtl validation. the details of the peak markers of the major effect qtl are given in additional file twelve major effect drought grain yield qtl were validated on a panel of drought tolerant lines consisting of traditional drought tolerant donors, drought tolerant breeding lines developed through conventional breeding approaches and random high yielding lines under drought from qtl mapping populations. the peak marker of all the twelve major effect qtl were amplified on the drought panel lines. the lines were scored taking qtl donor allele as a base. the list of lines is given in the additional file gene content analysis the meta-qtl were analyzed for gene content to know the presence of genes and gene clusters responsible for drought. a comparative genomics approach was followed to analyze the genes present in meta-qtl. gene content was noted based on annotated data of homologous regions in nipponbare using rap, build. it is assumed that the genes identified in nipponbare regions are homologous and collinear to those underlying the yield qtl under drought mapped in different studies involving different donors and recipients. comparative genomics to identify homologous regions in cereals a comparative genomics approach was followed to identify homologous regions between rice and maize using the genomic databases. homologous regions identified were checked for the presence of drought grain yield qtl of maize. in sorghum, wheat, and barley, grain yield qtl reported were collected from a literature survey and these were compared with the meta-qtl using the comparative maps available in the gramene database.  overview of qtl and development of a consensus map in the populations of rice screened for drought tolerance to map qtl, population size ranged from to lines. the number of markers used ranged from to. the number of locations for phenotyping varied from to from the studies, yield qtl were reported, which were distributed on all the chromosomes except chromosome. the number of qtl per population ranged from to the proportion of qtl per chromosome ranged from one qtl each on chromosomes and to yield qtl on chromosome the distribution of yield qtl on different chromosomes showed that chromosomes and have the highest number, and qtl, respectively. the phenotypic variance of the initial qtl varied from % to 40% and the confidence interval of the markers varied from to cm. the rice genetic map of temnykh et al. was used as a reference map to develop a consensus map as this is a widely used genetic map of rice and it contained most of the markers used in the different studies. the consensus map consisted of markers with a total map length of cm. the average distance between the markers was cm, thus enabling the identification of a precise location of qtl. there were very few marker inversions in the consensus map, which were discarded from the final map and further analysis. meta-analysis and qtl validation it is widely believed that qtl are accurate and can be positioned onto chromosomal locations by molecular mapping. however, their complex nature and context dependency in different genetic backgrounds and environments are constraints in identifying their precise location. the identification of the most accurate and precise major-effect qtl across genetic backgrounds and environments is a prerequisite for the successful use of qtl in mas. meta-analysis of qtl identified in different studies helps to identify the most precise and concise qtl, which can be further pursued for mas or the identification of candidate genes. in our study, we attempted to identify the meta-qtl for grain yield under drought by genome-wide meta-analysis. from a bibliographic survey, a total of qtl were short-listed for grain yield under drought from studies. all qtl were projected on a consensus map. the chromosomal regions with only one qtl were not considered for meta-analysis since meta-analysis by definition involves more than one qtl. thus, qtl were used for meta-analysis and meta-qtl were short-listed based on the akaike information criterion. accordingly, the qtl model with the lowest aic value was considered a significant model indicating the number of meta-qtl. the number of meta-qtl along with their aic values and confidence intervals are given in table in total, independent meta-qtl were identified at a confidence interval of 95% on seven chromosomes, and meta-analysis successfully reduced the total qtl by 63%. the meta-qtl identified on each chromosome varied from to there were four meta-qtl on chromosome 1; two on chromosomes and 10; and one each on chromosomes and the phenotypic variance of the meta-qtl varied from 4% to 28%. at of the meta-qtl, the mean phenotypic variance was more than 10%. in general, the confidence intervals at most of the meta-qtl were narrower than their respective original qtl. at nine loci on chromosomes and meta-qtl were narrower than the mean of their initial qtl. however, at five loci, the meta-qtl were broader than the mean of the initial qtl. the confidence intervals of the meta-qtl varied from cm between the marker intervals rg and rm on chromosome to cm between the marker intervals rm and rm on chromosome at two regions, meta-qtl and mqtl, the ci declined to around cm. the physical intervals of the meta-qtl varied from mb to mb. three meta-qtl were less than kb. the meta-qtl regions with small genetic and physical intervals are useful in mas. it is significant to see that seven qtl that had less than mb intervals also had a genetic interval of around cm with a phenotypic variance of more than 10%. three of these meta-qtl were on chromosome and one each on chromosomes and the physical intervals of mqtl, mqtl, and mqtl were less than kb, that of mqtl was kb, and those of mqtl, mqtl, and mqtl were 1mb, mb, and mb, respectively. mqtl, mqtl, and mqtl had phenotypic variance of more than 20%. the seven mqtl regions with small genetic and physical intervals are important regions for mas, fine mapping, candidate gene identification, and functional analysis. these qtl can be introgressed in popular rice mega-varieties to develop drought-tolerant and high-yielding lines. in addition to meta-analysis of qtl, the markers linked to the major-effect qtl for grain yield were also validated on a panel of drought-tolerant lines to confirm their presence in larger set of lines. it is notable that major-effect qtl dty was present in 85% of the lines. dty, dty dty, dty, and dty were present in more than 50% of the lines. the amplification of the rm and rm peak markers of dty and dty in a set of drought tolerant panel lines is presented in additional file the result indicates the presence of at least one of the major-effect grain yield qtl in the drought panel lines. in general, the major-effect qtl identified for grain yield under drought have a genetic gain of 10% to 30%, with a yield advantage of around to kg/ha over recipient parents. however, considering practical benefit to farmers, the development of drought-tolerant rice varieties with a yield advantage of at least ton/ha could be the desired target for rice breeders. the marker-aided qtl pyramiding of the major-effect mqtl identified in this study can be considered as an option for achieving this target. aic = akaike information criterion, ci = confidence interval, cm = centimorgan, mqtl = meta-qtl; kb = kilobase, mb = megabase a comparison was made between the meta-qtl identified in this study with the meta-qtl identified for root traits in two earlier studies. it is very interesting to note that mqtl, mqtl, mqtl, mqtl, and mqtl coincided with qtl clusters for root and leaf morphology traits associated with drought tolerance/avoidance in rice. all the independent meta-qtl coincided with at least one meta-qtl identified for root traits under drought. earlier studies on meta-analysis of qtl for root traits and blast resistance in rice, fusarium head blight resistance in wheat, flowering time in maize, nematode resistance in soybean, and lint fiber length in cotton identified precise and concise meta-qtl. meta-qtl were also used to deduce candidate genes and were recommended for mas in some of these studies. comparative genomics of mqtl the existence of an evolutionary relationship among the grass families is a well-known fact. the syntenic relationship can be used to identify the homologous regions among these species, which in turn is useful in defining their role in plant growth, development, and adaptation across species. we compared meta-qtl regions for synteny in other cereal crops. the major-effect mqtl was also found in maize on chromosome near marker msu in wheat on chromosome 4b near marker rht-b and in barley on chromosome 6h near marker bmac while major-effect mqtl was also found in maize on chromosome near marker umc107a. all these markers were linked to grain yield under drought in their respective crops. the largest parts of chromosomes and of rice have a syntenic relation with chromosomes and of maize, so their respective homologous qtl were also found on the corresponding chromosomes. an interesting observation is that, near the sd locus on chromosome of rice, qtl for grain yield under drought were identified most frequently. sd is a major locus responsible for semidwarf plant stature in rice and its corresponding locus in wheat is rht-b on chromosome 4b. mqtl is near the sd locus and also on its corresponding locus rht-b in wheat, major qtl for grain yield under drought were detected. gene content analysis and identification of candidate genes meta qtl with precise and narrow confidence intervals are useful in short listing the candidate genes. using the annotated gene information available in the rice database, the genes present in the meta-qtl regions were analyzed by comparative genomics approach and candidate genes were shortlisted. the short-listed candidate genes can be further confirmed by transgenic approaches by loss or gain of function studies. most of the genes present in the mqtl were genes for hypothetical and expressed proteins, pseudo genes, genes for signal transduction, and transposable elements. however, there were many annotated genes/gene families that were common across the mqtl regions; these are probable candidate genes for yield under drought. it was found that three kinds of genes frequently occurred together in these regions. the genes/gene families were stress-inducible genes, growth and development-related genes, and sugar transport-related genes. table lists the important genes underlying mqtl for grain yield under drought. in six mqtl with less than a mb region, lrr kinase, leucine zipper, cell division-controlling proteins, sugar transport protein-like genes, no apical meristem, pentatricopeptide repeat proteins, cytokinin oxidase, f-box proteins, ap2-domain containing proteins, and zinc-finger transcription factors were present. the candidacy of these genes for yield and yield traits has already been proved in rice and other crops. cytochrome p has a role in bassinosteroid homeostasis and had an influence on leaf angle leading to increased yield in rice. pentatricopeptide repeats are present in the promoter region of rf genes, which restore fertility and also play a role in embryogenesis in arabidopsis. zinc-finger -like proteins are known to be involved in stress tolerance. zinc-finger protein in rice are induced after different types of stresses, namely, cold, desiccation, salt, submergence, heavy metals, and mechanical injury. over expression of the zinc-finger gene in transgenic tobacco conferred tolerance of cold, dehydration, and salt stress at the seed germination/seedling stage. f-box proteins play an important role in floral development and stress tolerance. in addition, f-box proteins appear to serve as the key components of the machinery involved in regulating plant growth and development throughout the plant's life cycle and their expression is influenced by light and abiotic stresses. leucine zippers are a class of transcription factor involved in aba-independent stress tolerance. over expression of osbzip in rice triggered clusters of genes regulating stress adaptations. the no apical meristem gene plays an important role in the growth and development of meristematic tissue. the root-specific expression of this gene resulted in enhanced root growth and improved drought tolerance in rice. the other important genes that harbored the meta-qtl were the erecta and dreb genes. erecta is a leucine-rich repeat receptor-like kinase gene known for its influence on inflorescence development, stomatal density, epidermal cell expansion, and mesophyll cell proliferation. this gene is mainly involved in transpiration efficiency and enhanced drought response. dreb is a well-known transcription factor that is induced by drought and it activates many down stream stress-responsive genes to ultimately improve the drought and chilling tolerance of rice. some of these short-listed genes can be considered as positional candidate genes that determine grain yield under drought. however, it is also well known that yield and adaptability to stress are complex in nature and highly negatively correlated. the qtl/genes for these two are often co-located. even though individual genes have been proved to regulate yield under controlled drought experiments, a well-coordinated response of many genes is essential for drought tolerance under field conditions. this is evident from the presence of three different groups of gene clusters in most of the meta-qtl regions.  meta-analysis of grain yield qtl is an effective approach in identifying concise and precise consensus qtl. the seven meta-qtl identified with small genetic and physical intervals could be useful in mas/pyramiding. validation of the major-effect qtl confirmed the consistency of the major-effect grain yield qtl under drought in different drought-tolerant panel lines. the comparative genomics approach to identify the consistency of drought grain yield qtl across species revealed the conservation of some of the loci, indicating their evolutionary significance. the presence of gene clusters in the meta-qtl indicates that a well-coordinated response of many genes is essential to achieve drought tolerance under field conditions. authors' contributions ak conceived the idea of meta- analysis, qtl validation and comparative genomics of grain yield qtl under drought. bpms compiled and analyzed the data, carried out the qtl validation and comparative genomics. sd helped in data compilation and analysis. bpms, pv, sd, hua and ak were responsible for drafting and editing the manuscript. all authors have read and approved the final manuscript. supplementary material additional file 1 details of the markers used for qtl validation. this file contains the list of major effect qtls for grain yield under drought and peak markers of the qtls. primer sequence, product size of the markers and annealing temperatures used for amplifying the markers. click here for file additional file 2 drought panel lines for qtl validation. this table shows the list of drought panel lines and type of the breeding material. these lines were used for validating the major effect qtls for grain yield under drought. click here for file additional file 3 amplification of rm and rm peak markers of qtl and qtl in a set of drought tolerant panel lines. the gel picture shows the amplication of rm and rm peak markers of qtl and qtl in a set of drought tolerant panel lines. click here for file acknowledgements financial support to this study was provided by the bill & melinda gates foundation, usa.
the meta-qtl for grain yield under drought coincided with at least one of the meta-qtl identified for root and leaf morphology traits under drought in earlier reports. the integration of maps resulted in a consensus map with markers and a total map length of cm. the genes in the meta-qtl regions were analyzed by a comparative genomics approach and candidate genes were deduced for grain yield under drought. comparative genomics of meta-qtl with other cereals revealed that the homologous regions of mqtl and mqtl had qtl for grain yield under drought in maize, wheat, and barley respectively. meta-qtl with small genetic and physical intervals could be useful in marker-assisted selection individually and in combinations. validation and comparative genomics of the major-effect qtl confirmed their consistency within and across the species. however, identification of most precise and consistent qtl across the environments and genetics backgrounds is essential for their successful use in marker-assisted selection. in this study, an attempt was made to locate consistent qtl regions associated with yield increase under drought by applying a genome-wide qtl meta-analysis approach. mqtl, mqtl, mqtl, and mqtl were around kb and corresponded to a reasonably small genetic distance of to cm and they are suitable for use in marker-assisted selection. validation of major-effect qtl on a panel of random drought-tolerant lines revealed the presence of at least one major qtl in each line. in the last few years, efforts have been made to identify large effect qtl for grain yield under drought in rice. the shortlisted candidate genes can be cloned to unravel the molecular mechanism regulating grain yield under drought. three groups of genes such as stress-inducible genes, growth and development-related genes, and sugar transport-related genes were found in clusters in most of the meta-qtl. fifty-three yield qtl reported in studies were projected on a consensus map and meta-analysis was performed.
 in the last few years, efforts have been made to identify large effect qtl for grain yield under drought in rice. however, identification of most precise and consistent qtl across the environments and genetics backgrounds is essential for their successful use in marker-assisted selection. in this study, an attempt was made to locate consistent qtl regions associated with yield increase under drought by applying a genome-wide qtl meta-analysis approach.  the integration of maps resulted in a consensus map with markers and a total map length of cm. fifty-three yield qtl reported in studies were projected on a consensus map and meta-analysis was performed. fourteen meta-qtl were obtained on seven chromosomes. mqtl, mqtl, mqtl, and mqtl were around kb and corresponded to a reasonably small genetic distance of to cm and they are suitable for use in marker-assisted selection. the meta-qtl for grain yield under drought coincided with at least one of the meta-qtl identified for root and leaf morphology traits under drought in earlier reports. validation of major-effect qtl on a panel of random drought-tolerant lines revealed the presence of at least one major qtl in each line. dty was present in 85% of the lines, followed by dty in 79% and dty in 64% of the lines. comparative genomics of meta-qtl with other cereals revealed that the homologous regions of mqtl and mqtl had qtl for grain yield under drought in maize, wheat, and barley respectively. the genes in the meta-qtl regions were analyzed by a comparative genomics approach and candidate genes were deduced for grain yield under drought. three groups of genes such as stress-inducible genes, growth and development-related genes, and sugar transport-related genes were found in clusters in most of the meta-qtl.  meta-qtl with small genetic and physical intervals could be useful in marker-assisted selection individually and in combinations. validation and comparative genomics of the major-effect qtl confirmed their consistency within and across the species. the shortlisted candidate genes can be cloned to unravel the molecular mechanism regulating grain yield under drought. 