Effect of epistasis and environment on flowering time in barley reveals a novel flowering-delaying QTL allele.

Flowering time is a complex trait and has a key role in crop yield and adaptation to environmental stressors such as heat and drought. This study aimed to better understand the interconnected dynamics of epistasis and environment and look for novel regulators. We investigated 534 spring barley MAGIC DH lines for flowering time at various environments. Analysis of quantitative trait loci (QTLs), epistatic interactions, QTL × environment (Q×E) interactions, and epistasis × environment (E×E) interactions were performed with single SNP and haplotype approaches. In total, 18 QTLs and 2420 epistatic interactions were detected, including intervals harboring major genes such as Ppd-H1, Vrn-H1, Vrn-H3, and denso/sdw1. Epistatic interactions found in field and semi-controlled conditions were distinctive. Q×E and E×E interactions revealed that temperature influenced flowering time by triggering different interactions between known and newly detected regulators. A novel flowering-delaying QTL allele was identified on chromosome 1H (named 'HvHeading') and was shown to be engaged in epistatic and environment interactions. Results suggest that investigating epistasis, environment, and their interactions, rather than only single QTLs, is an effective approach for detecting novel regulators. We assume that barley can adapt flowering time to the environment via alternative routes within the pathway.

Adaptive selection of founder segments and epistatic control of plant height in the MAGIC winter wheat population WM-800.

BACKGROUND: Multi-parent advanced generation intercross (MAGIC) populations are a newly established tool to dissect quantitative traits. We developed the high resolution MAGIC wheat population WM-800, consisting of 910 F4:6 lines derived from intercrossing eight recently released European winter wheat cultivars. RESULTS: Genotyping WM-800 with 7849 SNPs revealed a low mean genetic similarity of 59.7% between MAGIC lines. WM-800 harbours distinct genomic regions exposed to segregation distortion. These are mainly located on chromosomes 2 to 6 of the wheat B genome where founder specific DNA segments were positively or negatively selected. This suggests adaptive selection of individual founder alleles during population development. The application of a genome-wide association study identified 14 quantitative trait loci (QTL) controlling plant height in WM-800, including the known semi-dwarf genes Rht-B1 and Rht-D1 and a potentially novel QTL on chromosome 5A. Additionally, epistatic effects controlled plant height. For example, two loci on chromosomes 2B and 7B gave rise to an additive epistatic effect of 13.7 cm. CONCLUSION: The present study demonstrates that plant height in the MAGIC-WHEAT population WM-800 is mainly determined by large-effect QTL and di-genic epistatic interactions. As a proof of concept, our study confirms that WM-800 is a valuable tool to dissect the genetic architecture of important agronomic traits.

Genome-Wide Association Study in Wheat Identifies Resistance to the Cereal Cyst Nematode Heterodera filipjevi.

The cyst nematode Heterodera filipjevi is a plant parasite causing substantial yield loss in wheat. Resistant cultivars are the preferred method of controlling cyst nematodes. Association mapping is a powerful approach to detect associations between phenotypic variation and genetic polymorphisms; in this way favorable traits such as resistance to pathogens can be located. Therefore, a genome-wide association study of 161 winter wheat accessions was performed with a 90K iSelect single nucleotide polymorphism (SNP) chip. Population structure analysis grouped into two major subgroups and first principal component accounted 6.16% for phenotypic diversity. The genome-wide linkage disequilibrium across wheat was 3 cM. Eleven quantitative trait loci (QTLs) on chromosomes 1AL, 2AS, 2BL, 3AL, 3BL, 4AS, 4AL, 5BL, and 7BL were identified using a mixed linear model false discovery rate of P < 0.01 that explained 43% of total genetic variation. This is the first report of QTLs conferring resistance to H. filipjevi in wheat. Eight QTLs on chromosomes 1AL, 2AS, 2BL, 3AL, 4AL, and 5BL were linked to putative genes known to be involved in plant-pathogen interactions. Two other QTLs on 3BL and one QTL on 7BL linked to putative genes known to be involved in abiotic stress.

Identifying developmental QTL alleles with favorable effect on grain yield components under late-terminal drought in spring barley MAGIC population.

Barley is the fourth most cultivated cereal worldwide, and drought is a major cause of its yield loss by negatively affecting its development. Hence, better understanding developmental mechanisms that control complex polygenic yield-related traits under drought is essential to uncover favorable yield regulators. This study evaluated seven above-ground yield-related traits under well-watered (WW) and late-terminal drought (TD) treatment using 534 spring barley multiparent advanced generation intercross double haploid (DH) lines. The analysis of quantitative trait loci (QTL) for WW, TD, marker by treatment interaction, and drought stress tolerance identified 69, 64, 25, and 25 loci, respectively, for seven traits from which 15 loci were common for at least three traits and 17 were shared by TD and drought stress tolerance. Evaluation of allelic effects for a QTL revealed varying effect of parental alleles. Results showed prominent QTL located on major flowering time gene Ppd-H1 with favorable effects for grain weight under TD when flowering time was not significantly affected, suggesting that this gene might be linked with increasing grain weight by ways other than timing of flowering under late-terminal drought stress. Furthermore, a desirable novel QTL allele was identified on chromosome 5H for grain number under TD nearby sucrose transporter gene HvSUT2. The findings indicated that spring barley multiparent advanced generation intercross population can provide insights to improve yield under complex condition of drought.

Identifying genomic regions determining shoot and root traits related to nitrogen uptake efficiency in a multiparent advanced generation intercross (MAGIC) winter wheat population in a high-throughput phenotyping facility.

In the context of a continuously increasing human population that needs to be fed, with environmental protection in mind, nitrogen use efficiency (NUE) improvement is becoming very important. To understand the natural variation of traits linked to nitrogen uptake efficiency (UPE), one component of NUE, the multiparent advanced generation intercross (MAGIC) winter wheat population WM-800 was phenotyped under two contrasting nitrogen (N) levels in a high-throughput phenotyping facility for six weeks. Three biomass-related, three root-related, and two reflectance-related traits were measured weekly under each treatment. Subsequently, the population was genetically analysed using a total of 13,060 polymorphic haplotypes and singular SNPs for a genome-wide association study (GWAS). In total, we detected 543 quantitative trait loci (QTL) across all time points and traits, which were pooled into 42 stable QTL (sQTL; present in at least three of the six weeks). Besides Rht-B1 and Rht-D1, candidate genes playing a role in gibberellic acid-regulated growth and nitrate transporter genes from the NPF gene family, like NRT 1.1, were linked to sQTL. Two novel sQTL on chromosomes 5 A and 6D showed pleiotropic effects on several traits. The high number of N-specific sQTL indicates that selection for UPE is useful specifically under N-limited conditions.

Unraveling Genomic Regions Controlling Root Traits as a Function of Nitrogen Availability in the MAGIC Wheat Population WM-800.

An ever-growing world population demands to be fed in the future and environmental protection and climate change need to be taken into account. An important factor here is nitrogen uptake efficiency (NUpE), which is influenced by the root system (the interface between plant and soil). To understand the natural variation of root system architecture (RSA) as a function of nitrogen (N) availability, a subset of the multiparent advanced generation intercross (MAGIC) winter wheat population WM-800 was phenotyped under two contrasting N treatments in a high-throughput phenotyping system at the seedling stage. Fourteen root and shoot traits were measured. Subsequently, these traits were genetically analyzed using 13,060 polymorphic haplotypes and SNPs in a genome-wide association study (GWAS). In total, 64 quantitative trait loci (QTL) were detected; 60 of them were N treatment specific. Candidate genes for the detected QTL included NRT1.1 and genes involved in stress signaling under N-, whereas candidate genes under N+ were more associated with general growth, such as mei2 and TaWOX11b. This finding may indicate (i) a disparity of the genetic control of root development under low and high N supply and, furthermore, (ii) the need for an N specific selection of genes and genotypes in breeding new wheat cultivars with improved NUpE.

A Haplotype-Based GWAS Identified Trait-Improving QTL Alleles Controlling Agronomic Traits under Contrasting Nitrogen Fertilization Treatments in the MAGIC Wheat Population WM-800.

The multi-parent-advanced-generation-intercross (MAGIC) population WM-800 was developed by intercrossing eight modern winter wheat cultivars to enhance the genetic diversity present in breeding populations. We cultivated WM-800 during two seasons in seven environments under two contrasting nitrogen fertilization treatments. WM-800 lines exhibited highly significant differences between treatments, as well as high heritabilities among the seven agronomic traits studied. The highest-yielding WM-line achieved an average yield increase of 4.40 dt/ha (5.2%) compared to the best founder cultivar Tobak. The subsequent genome-wide-association-study (GWAS), which was based on haplotypes, located QTL for seven agronomic traits including grain yield. In total, 40, 51, and 46 QTL were detected under low, high, and across nitrogen treatments, respectively. For example, the effect of QYLD_3A could be associated with the haplotype allele of cultivar Julius increasing yield by an average of 4.47 dt/ha (5.2%). A novel QTL on chromosome 2B exhibited pleiotropic effects, acting simultaneously on three-grain yield components (ears-per-square-meter, grains-per-ear, and thousand-grain-weight) and plant-height. These effects may be explained by a member of the nitrate-transporter-1 (NRT1)/peptide-family, TaNPF5.34, located 1.05 Mb apart. The WM-800 lines and favorable QTL haplotypes, associated with yield improvements, are currently implemented in wheat breeding programs to develop advanced nitrogen-use efficient wheat cultivars.

Detection of Epistasis for Flowering Time Using Bayesian Multilocus Estimation in a Barley MAGIC Population.

Gene-by-gene interactions, also known as epistasis, regulate many complex traits in different species. With the availability of low-cost genotyping it is now possible to study epistasis on a genome-wide scale. However, identifying genome-wide epistasis is a high-dimensional multiple regression problem and needs the application of dimensionality reduction techniques. Flowering Time (FT) in crops is a complex trait that is known to be influenced by many interacting genes and pathways in various crops. In this study, we successfully apply Sure Independence Screening (SIS) for dimensionality reduction to identify two-way and three-way epistasis for the FT trait in a Multiparent Advanced Generation Inter-Cross (MAGIC) barley population using the Bayesian multilocus model. The MAGIC barley population was generated from intercrossing among eight parental lines and thus, offered greater genetic diversity to detect higher-order epistatic interactions. Our results suggest that SIS is an efficient dimensionality reduction approach to detect high-order interactions in a Bayesian multilocus model. We also observe that many of our findings (genomic regions with main or higher-order epistatic effects) overlap with known candidate genes that have been already reported in barley and closely related species for the FT trait.