One of the biggest challenges for genetic studies on natural or unstructured populations is the unbalanced datasets where individuals are measured at different times and environments. This problem is also common in crop and animal breeding where many individuals are only evaluated for a single year and large but unbalanced datasets can be generated over multiple years. Many wheat breeding programs have focused on increasing bread wheat (Triticum aestivum L.) yield, but processing and end-use quality are critical components when considering its use in feeding the rising population of the next century. The challenges with end-use quality trait improvements are high cost and seed amounts for testing, the latter making selection in early breeding populations impossible. Here we describe a novel approach to identify marker-trait associations within a breeding program using a meta-genome wide association study (GWAS), which combines GWAS analysis from multi-year unbalanced breeding nurseries, in a manner reflecting meta-GWAS in humans. This method facilitated mapping of processing and end-use quality phenotypes from advanced breeding lines (n = 4,095) of the CIMMYT bread wheat breeding program from 2009 to 2014. Using the meta-GWAS we identified marker-trait associations, allele effects, candidate genes, and can select using markers generated in this process. Finally, the scope of this approach can be broadly applied in ‘breeding-assisted genomics’ across many crops to greatly increase our functional understanding of plant genomes.

Tocochromanols are recognized for nutritional content, plant stress response, and seed longevity. Here we present a systems biological approach to characterize and develop predictive assays for genes affecting tocochromanol variation in barley. Major QTL, detected in three regions of a SNP linkage map, affected multiple tocochromanol forms. Candidate genes were identified through barley/rice orthology and sequenced in genotypes with disparate tocochromanol profiles. Gene-specific markers, designed based on observed polymorphism, mapped to the originating QTL, increasing R² values at the respective loci. Polymorphism within promoter regions corresponded to motifs known to influence gene expression. Quantitative PCR analysis revealed a trend of increased expression in tissues grown at cold temperatures. These results demonstrate utility of a novel method for rapid gene identification and characterization, and provide a resource for efficient development of barley lines with improved tocochromanol profiles.