QTL Mapping for Wheat Flour Color with Additive, Epistatic, and QTL×Environmental Interaction Effects

Mapping of QTLs and meta-QTLs for Heterodera avenae Woll. resistance in common wheat (Triticum aestivum L.)

BACKGROUND

In hexaploid wheat, quantitative trait loci (QTL) and meta-QTL (MQTL) analyses were conducted to identify genomic regions controlling resistance to cereal cyst nematode (CCN), Heterodera avenae. A mapping population comprising 149 RILs derived from the cross HUW 468 × C 306 was used for composite interval mapping (CIM) and inclusive composite interval mapping (ICIM).

RESULTS

Eight main effect QTLs on three chromosomes (1B, 2A and 3A) were identified using two repeat experiments. One of these QTLs was co-localized with a previously reported wheat gene Cre5 for resistance to CCN. Seven important digenic epistatic interactions (PVE = 5% or more) were also identified, each involving one main effect QTL and another novel E-QTL. Using QTLs earlier reported in literature, two meta-QTLs were also identified, which were also used for identification of 57 candidate genes (CGs). Out of these, 29 CGs have high expression in roots and encoded the following proteins having a role in resistance to plant parasitic nematodes (PPNs): (i) NB-ARC,P-loop containing NTP hydrolase, (ii) Protein Kinase, (iii) serine-threonine/tyrosine-PK, (iv) protein with leucine-rich repeat, (v) virus X resistance protein-like, (vi) zinc finger protein, (vii) RING/FYVE/PHD-type, (viii) glycosyl transferase, family 8 (GT8), (ix) rubisco protein with small subunit domain, (x) protein with SANT/Myb domain and (xi) a protein with a homeobox.

CONCLUSION

Identification and selection of resistance loci with additive and epistatic effect along with two MQTL and associated CGs, identified in the present study may prove useful for understanding the molecular basis of resistance against H. avenae in wheat and for marker-assisted selection (MAS) for breeding CCN resistant wheat cultivars.

Genetics of spot blotch resistance in bread wheat (Triticum aestivum L.) using five models for GWAS

Genetic architecture of resistance to spot blotch in wheat was examined using a Genome-Wide Association Study (GWAS) involving an association panel comprising 303 diverse genotypes. The association panel was evaluated at two different locations in India including Banaras Hindu University (BHU), Varanasi (Uttar Pradesh), and Borlaug Institute for South Asia (BISA), Pusa, Samastipur (Bihar) for two consecutive years (2017-2018 and 2018-2019), thus making four environments (E1, BHU 2017-18; E2, BHU 2018-19; E3, PUSA, 2017-18; E4, PUSA, 2018-19). The panel was genotyped for 12,196 SNPs based on DArT-seq (outsourced to DArT Ltd by CIMMYT); these SNPs included 5,400 SNPs, which could not be assigned to individual chromosomes and were therefore, described as unassigned by the vendor. Phenotypic data was recorded on the following three disease-related traits: (i) Area Under Disease Progress Curve (AUDPC), (ii) Incubation Period (IP), and (iii) Lesion Number (LN). GWAS was conducted using each of five different models, which included two single-locus models (CMLM and SUPER) and three multi-locus models (MLMM, FarmCPU, and BLINK). This exercise gave 306 MTAs, but only 89 MTAs (33 for AUDPC, 30 for IP and 26 for LN) including a solitary MTA detected using all the five models and 88 identified using four of the five models (barring SUPER) were considered to be important. These were used for further analysis, which included identification of candidate genes (CGs) and their annotation. A majority of these MTAs were novel. Only 70 of the 89 MTAs were assigned to individual chromosomes; the remaining 19 MTAs belonged to unassigned SNPs, for which chromosomes were not known. Seven MTAs were selected on the basis of minimum P value, number of models, number of environments and location on chromosomes with respect to QTLs reported earlier. These 7 MTAs, which included five main effect MTAs and two for epistatic interactions, were considered to be important for marker-assisted selection (MAS). The present study thus improved our understanding of the genetics of resistance against spot blotch in wheat and provided seven MTAs, which may be used for MAS after due validation.

High-Resolution Mapping in Two RIL Populations Refines Major "QTL Hotspot" Regions for Seed Size and Shape in Soybean (Glycine max L.)

Seed size and shape are important traits determining yield and quality in soybean. However, the genetic mechanism and genes underlying these traits remain largely unexplored. In this regard, this study used two related recombinant inbred line (RIL) populations (ZY and K3N) evaluated in multiple environments to identify main and epistatic-effect quantitative trait loci (QTLs) for six seed size and shape traits in soybean. A total of 88 and 48 QTLs were detected through composite interval mapping (CIM) and mixed-model-based composite interval mapping (MCIM), respectively, and 15 QTLs were common among both methods; two of them were major (R2 > 10%) and novel QTLs (viz., qSW-1-1ZN and qSLT-20-1K3N). Additionally, 51 and 27 QTLs were identified for the first time through CIM and MCIM methods, respectively. Colocalization of QTLs occurred in four major QTL hotspots/clusters, viz., "QTL Hotspot A", "QTL Hotspot B", "QTL Hotspot C", and "QTL Hotspot D" located on Chr06, Chr10, Chr13, and Chr20, respectively. Based on gene annotation, gene ontology (GO) enrichment, and RNA-Seq analysis, 23 genes within four "QTL Hotspots" were predicted as possible candidates, regulating soybean seed size and shape. Network analyses demonstrated that 15 QTLs showed significant additive x environment (AE) effects, and 16 pairs of QTLs showing epistatic effects were also detected. However, except three epistatic QTLs, viz., qSL-13-3ZY, qSL-13-4ZY, and qSW-13-4ZY, all the remaining QTLs depicted no main effects. Hence, the present study is a detailed and comprehensive investigation uncovering the genetic basis of seed size and shape in soybeans. The use of a high-density map identified new genomic regions providing valuable information and could be the primary target for further fine mapping, candidate gene identification, and marker-assisted breeding (MAB).

Detection of stable QTLs for grain protein content in rice (Oryza sativa L.) employing high throughput phenotyping and genotyping platforms

Lack of appropriate donors, non-utilization of high throughput phenotyping and genotyping platforms with high genotype × environment interaction restrained identification of robust QTLs for grain protein content (GPC) in rice. In the present investigation a  BC₃F₄ mapping population was developed using grain protein donor, ARC10075 and high-yielding cultivar Naveen and 190 lines were genotyped using 40 K Affimetrix custom SNP array with the objective to identify stable QTLs for protein content. Three of the identified QTLs, one for GPC (qGPC1.1) and the other two for single grain protein content (qSGPC2.1, qSGPC7.1) were stable over the environments explaining  13%, 14% and 7.8% of the phenotypic variances, respectively. Stability and repeatability of these additive QTLs were supported by the synergistic additive effects of multi-environmental-QTLs. One epistatic-QTL, independent of  the  main effect QTL was detected over the environment for SGPC. A few functional genes governing seed storage protein were hypothesised inside these identified QTLs. The qGPC1.1 was validated by NIR Spectroscopy-based high throughput phenotyping in BC₃F₅ population. Higher glutelin content was estimated in high-protein lines with the introgression of qGPC1.1 in telomeric region of short arm of chromosome 1. This was supported by the postulation of probable candidate gene inside this QTL region encoding glutelin family proteins.

Genome-Wide Detection of Major and Epistatic Effect QTLs for Seed Protein and Oil Content in Soybean Under Multiple Environments Using High-Density Bin Map

Seed protein and oil content are the two important traits determining the quality and value of soybean. Development of improved cultivars requires detailed understanding of the genetic basis underlying the trait of interest. However, it is prerequisite to have a high-density linkage map for precisely mapping genomic regions, and therefore the present study used high-density genetic map containing 2267 recombination bin markers distributed on 20 chromosomes and spanned 2453.79 cM with an average distance of 1.08 cM between markers using restriction-site-associated DNA sequencing (RAD-seq) approach. A recombinant inbred line (RIL) population of 104 lines derived from a cross between Linhefenqingdou and Meng 8206 cultivars was evaluated in six different environments to identify main- and epistatic-effect quantitative trait loci (QTLs)as well as their interaction with environments. A total of 44 main-effect QTLs for protein and oil content were found to be distributed on 17 chromosomes, and 15 novel QTL were identified for the first time. Out of these QTLs, four were major and stable QTLs, viz., qPro-7-1, qOil-8-3, qOil-10-2 and qOil-10-4, detected in at least two environments plus combined environment with R² values >10%. Within the physical intervals of these four QTLs, 111 candidate genes were screened for their direct or indirect involvement in seed protein and oil biosynthesis/metabolism processes based on gene ontology and annotation information. Based on RNA sequencing (RNA-seq) data analysis, 15 of the 111 genes were highly expressed during seed development stage and root nodules that might be considered as the potential candidate genes. Seven QTLs associated with protein and oil content exhibited significant additive and additive × environment interaction effects, and environment-independent QTLs revealed higher additive effects. Moreover, three digenic epistatic QTLs pairs were identified, and no main-effect QTLs showed epistasis. In conclusion, the use of a high-density map identified closely linked flanking markers, provided better understanding of genetic architecture and candidate gene information, and revealed the scope available for improvement of soybean quality through marker assisted selection (MAS).

Genome-wide linkage mapping of flour color-related traits and polyphenol oxidase activity in common wheat

KEY MESSAGE

Fifty-six QTL for flour color-related traits and polyphenol oxidase activity were identified using a genome-wide linkage mapping of data from a RIL population derived from a Gaocheng 8901/Zhoumai 16 cross.

ABSTRACT

Flour color-related traits, including L*, a*, b*, yellow pigment content (YPC), and polyphenol oxidase (PPO) activity are important parameters influencing the quality of wheat end-use products. Mapping quantitative trait loci (QTL) for these traits and characterization of candidate genes are important for improving wheat quality. The aims of this study were to identify QTL for flour color-related traits and PPO activity and to characterize candidate genes using a high-density genetic linkage map in a common wheat recombinant inbred line (RIL) population derived from a cross between Gaocheng 8901 and Zhoumai 16. A linkage map was constructed by genotyping the RILs with the wheat 90 K iSelect array. Fifty-six QTL were mapped on 35 chromosome regions on homoeologous groups 1, 2, 5 and 7 chromosomes, and chromosomes 3B, 4A, 4B and 6B. Four QTL were for PPO activity, and the others were for flour color-related traits. Compared with previous studies, five QTL for a*, two for b*, one for L*, one for YPC and one for PPO activity were new. The new QTL on chromosome 2DL was involved in both a* and YPC, and another on chromosome 7DS affected both a* and L*. The scan for SNP sequences tightly linked to QTL for flour color-related traits against the wheat and/or related cereals genomes identified six candidate genes significantly related to these traits, and five of them were associated with the terpenoid backbone biosynthesis pathway. The high-density genetic linkage map of Gaocheng 8901/Zhoumai 16 represents a useful tool to identify QTL for important quality traits and candidate genes.