QTL mapping for grain yield and three yield components in a population derived from two high-yielding spring wheat cultivars

Genetic factors of grain cadmium concentration in Polish wheat (Triticum polonicum L.)

Wheat (Triticum aestivum L.) is one of the most important crops worldwide and a major source of human cadmium (Cd) intake. Limiting grain Cd concentration (Gr_Cd_Conc) in wheat is necessary to ensure food safety. However, the genetic factors associated with Cd uptake, translocation and distribution and Gr_Cd_Conc in wheat are poorly understood. Here, we mapped quantitative trait loci (QTLs) for Gr_Cd_Conc and its related transport pathway using a recombinant inbred line (RIL) population derived from 2 Polish wheat varieties (RIL_DT; dwarf Polish wheat [DPW] and tall Polish wheat [TPW]). We identified 29 novel major QTLs for grain and tissue Cd concentration; 14 novel major QTLs for Cd uptake, translocation, and distribution; and 27 major QTLs for agronomic traits. We also analyzed the pleiotropy of these QTLs. Six novel QTLs (QGr_Cd_Conc-1A, QGr_Cd_Conc-3A, QGr_Cd_Conc-4B, QGr_Cd_Conc-5B, QGr_Cd_Conc-6A, and QGr_Cd_Conc-7A) for Gr_Cd_Conc explained 8.16% to 17.02% of the phenotypic variation. QGr_Cd_Conc-3A, QGr_Cd_Conc-6A, and QGr_Cd_Conc-7A pleiotropically regulated Cd transport; 3 other QTLs were organ-specific for Gr_Cd_Conc. We fine-mapped the locus of QGr_Cd_Conc-4B and identified the candidate gene as Cation/Ca exchanger 2 (TpCCX2-4B), which was differentially expressed in DPW and TPW. It encodes an endoplasmic reticulum membrane/plasma membrane-localized Cd efflux transporter in yeast. Overexpression of TpCCX2-4B reduced Gr_Cd_Conc in rice. The average Gr_Cd_Conc was significantly lower in TpCCX2-4BDPW genotypes than in TpCCX2-4BTPW genotypes of the RIL_DT population and 2 other natural populations, based on a Kompetitive allele-specific PCR marker derived from the different promoter sequences between TpCCX2-4BDPW and TpCCX2-4BTPW. Our study reveals the genetic mechanism of Cd accumulation in wheat and provides valuable resources for genetic improvement of low-Cd-accumulating wheat cultivars.

Identification of Quantitative Trait Loci Associated with Plant Adaptation Traits Using Nested Association Mapping Population

This study evaluated 290 recombinant inbred lines (RILs) of the nested association mapping (NAM) population from the UK. The population derived from 24 families, where a common parent was "Paragon," one of the UK's spring wheat cultivar standards. All genotypes were tested in two regions of Kazakhstan at the Kazakh Research Institute of Agriculture and Plant Industry (KRIAPI, Almaty region, Southeast Kazakhstan, 2019-2022 years) and Alexandr Barayev Scientific-Production Center for Grain Farming (SPCGF, Shortandy, Akmola region, Northern Kazakhstan, 2019-2022 years). The studied traits consisted of plant adaptation-related traits, including heading date (HD, days), seed maturation date (SMD, days), plant height (PH, cm), and peduncle length (PL, cm). In addition, the yield per m2 was analyzed in both regions. Based on a field evaluation of the population in northern and southeastern Kazakhstan and using 10,448 polymorphic SNP (single-nucleotide polymorphism) markers, the genome-wide association study (GWAS) allowed for detecting 74 QTLs in four studied agronomic traits (HD, SMD, PH, and PL). The literature survey suggested that 16 of the 74 QTLs identified in our study had also been detected in previous QTL mapping studies and GWASs for all studied traits. The results will be used for further studies related to the adaptation and productivity of wheat in breeding projects for higher grain productivity.

Quantitative trait loci detection for three tiller-related traits and the effects on wheat (Triticum aestivum L.) yields

KEY MESSAGE

A total of 38 putative additive QTLs and 55 pairwise putative epistatic QTLs for tiller-related traits were reported, and the candidate genes underlying qMtn-KJ-5D, a novel major and stable QTL for maximum tiller number, were characterized. Tiller-related traits play an important role in determining the yield potential of wheat. Therefore, it is important to elucidate the genetic basis for tiller number when attempting to use genetic improvement as a tool for enhancing wheat yields. In this study, a quantitative trait locus (QTL) analysis of three tiller-related traits was performed on the recombinant inbred lines (RILs) of a mapping population, referred to as KJ-RILs, that was derived from a cross between the Kenong 9204 (KN9204) and Jing 411 (J411) lines. A total of 38 putative additive QTLs and 55 pairwise putative epistatic QTLs for spike number per plant (SNPP), maximum tiller number (MTN), and ear-bearing tiller rate (EBTR) were detected in eight different environments. Among these QTLs with additive effects, three major and stable QTLs were first documented herein. Almost all but two pairwise epistatic QTLs showed minor interaction effects accounting for no more than 3.0% of the phenotypic variance. The genetic effects of two colocated major and stable QTLs, i.e., qSnpp-KJ-5D.1 and qMtn-KJ-5D, for yield-related traits were characterized. The breeding selection effect of the beneficial allele for the two QTLs was characterized, and its genetic effects on yield-related traits were evaluated. The candidate genes underlying qMtn-KJ-5D were predicted based on multi-omics data, and TraesKN5D01HG00080 was identified as a likely candidate gene. Overall, our results will help elucidate the genetic architecture of tiller-related traits and can be used to develop novel wheat varieties with high yields.

Simultaneous improvement of grain yield and grain protein concentration in durum wheat by using association tests and weighted GBLUP

KEY MESSAGE

Simultaneous improvement for GY and GPC by using GWAS and GBLUP suggested a significant application in durum wheat breeding. Despite the importance of grain protein concentration (GPC) in determining wheat quality, its negative correlation with grain yield (GY) is still one of the major challenges for breeders. Here, a durum wheat panel of 200 genotypes was evaluated for GY, GPC, and their derived indices (GPD and GYD), under eight different agronomic conditions. The plant material was genotyped with the Illumina 25 k iSelect array, and a genome-wide association study was performed. Two statistical models revealed dozens of marker-trait associations (MTAs), each explaining up to 30%. phenotypic variance. Two markers on chromosomes 2A and 6B were consistently identified by both models and were found to be significantly associated with GY and GPC. MTAs identified for phenological traits co-mapped to well-known genes (i.e., Ppd-1, Vrn-1). The significance values (p-values) that measure the strength of the association of each single nucleotide polymorphism marker with the target traits were used to perform genomic prediction by using a weighted genomic best linear unbiased prediction model. The trained models were ultimately used to predict the agronomic performances of an independent durum wheat panel, confirming the utility of genomic prediction, although environmental conditions and genetic backgrounds may still be a challenge to overcome. The results generated through our study confirmed the utility of GPD and GYD to mitigate the inverse GY and GPC relationship in wheat, provided novel markers for marker-assisted selection and opened new ways to develop cultivars through genomic prediction approaches.

Dissection and validation of a promising QTL controlling spikelet number on 5B in bread wheat

KEY MESSAGE

Five environmentally stable QTLs for spikelet number per spike and days to heading were identified using a high-genetic map containing 95,444 SNPs, among which QSns.ucas-5B was validated using residual heterozygous line at multiple environments. Spikelet number per spike (SNS) and days to heading (DTH) play pivotal roles in the improvement of wheat yield. In this study, a high-density genetic map for a recombinant inbred lines (RILs) population derived from Zhengnong 17 (ZN17) and Yangbaimai (YBM) was constructed using 95,444 single-nucleotide polymorphism (SNP) markers from the Wheat660K SNP array. Our study identified a total of five environmentally stable QTLs for SNS and DTH, one of which was named QSns.ucas-5B, with a physical interval of approximately 545.4-552.1 Mb on the 5BL chromosome arm. Importantly, the elite haplotype within QSns.ucas-5B showed a consistent and positive effect on SNS, grain number and weight per spike, without extending the days to heading. These findings provide a foundation for future efforts to map and clone the gene(s) responsible for QSns.ucas-5B and further indicate the potential application of the developed and validated InDel marker of QSns.ucas-5B for molecular breeding purposes, aimed at improving wheat grain yield.

A major QTL simultaneously increases the number of spikelets per spike and thousand-kernel weight in a wheat line

KEY MESSAGE

A novel and stably expressed QTL QSNS.sicau-SSY-7A for spikelet number per spike in wheat without negative effects on thousand-kernel weight was identified and validated in different genetic backgrounds. Spikelet number per spike (SNS) is an important determinant of yield in wheat. In the present study, we combined bulked segregant analysis (BSA) and the wheat 660 K single-nucleotide polymorphism (SNP) array to rapidly identify genomic regions associated with SNS from a recombinant inbred line (RIL) population derived from a cross between the wheat lines S849-8 and SY95-71. A genetic map was constructed using Kompetitive Allele Specific PCR markers in the SNP-enriched region on the long arm of chromosome 7A. A major and stably expressed QTL, QSNS.sicau-SSY-7A, was detected in multiple environments. It was located in a 1.6 cM interval on chromosome arm 7AL flanked by the markers AX-109983514 and AX-109820548. This QTL explained 6.86-15.72% of the phenotypic variance, with LOD values ranging from 3.66 to 8.66. Several genes associated with plant growth and development were identified in the interval where QSNS.sicau-SSY-7A was located on the 'Chinese Spring' wheat and wild emmer reference genomes. Furthermore, the effects of QSNS.sicau-SSY-7A and WHEAT ORTHOLOG OFAPO1(WAPO1) on SNS were analyzed. Interestingly, QSNS.sicau-SSY-7A significantly increased SNS without negative effects on thousand-kernel weight, anthesis date and plant height, demonstrating its great potential for breeding aimed at improving grain yield. Taken together, these results indicate that QSNS.sicau-SSY-7A is a promising locus for yield improvement, and its linkage markers are helpful for fine mapping and molecular breeding.

Genome-wide association mapping in elite winter wheat breeding for yield improvement

Increased grain yield (GY) is the primary breeding target of wheat breeders. We performed the genome-wide association study (GWAS) on 168 elite winter wheat lines from an ongoing breeding program to identify the main determinants of grain yield. Sequencing of Diversity Array Technology fragments (DArTseq) resulted in 19,350 single-nucleotide polymorphism (SNP) and presence-absence variation (PAV) markers. We identified 15 main genomic regions located in ten wheat chromosomes (1B, 2B, 2D, 3A, 3D, 5A, 5B, 6A, 6B, and 7B) that explained from 7.9 to 20.3% of the variation in grain yield and 13.3% of the yield stability. Loci identified in the reduced genepool are important for wheat improvement using marker-assisted selection. We found marker-trait associations between three genes involved in starch biosynthesis and grain yield. Two starch synthase genes (TraesCS2B03G1238800 and TraesCS2D03G1048800) and a sucrose synthase gene (TraesCS3D03G0024300) were found in regions of QGy.rut-2B.2, QGy.rut-2D.1, and QGy.rut-3D, respectively. These loci and other significantly associated SNP markers found in this study can be used for pyramiding favorable alleles in high-yielding varieties or to improve the accuracy of prediction in genomic selection.

Development and Identification of a Wheat-Rye Breeding Line for Harmonious Improvement Between Powdery Mildew Resistance and High Yield Potential

Powdery mildew, caused by Blumeria graminis f. sp. tritici, is a devastating disease that seriously threatens wheat yield and quality. To control this disease, host resistance is the preferred measure. However, wheat breeding is a complex process with elusive exchange and recombination of the traits from their parents. Increased resistance often leads to a decline in other key traits, such as yield and quality. Developing breakthrough germplasms with harmonious powdery mildew resistance and other key breeding traits is attractive in wheat breeding. In this study, we developed an ideal wheat breeding line AL46 that pyramided its hexaploid triticale parent-derived desirable yield traits and its wheat parent-derived powdery mildew resistance gene Pm2. Sequential genomic in situ hybridization (GISH), multicolor GISH, multicolor fluorescence in situ hybridization, and molecular marker analyses revealed that AL46 was a wheat-rye T1RS·1BL translocation line. Genetic analysis combined with function marker detection and sequence alignment were used to confirm that AL46 carried the Pm2 gene. Then, we evaluated the powdery mildew resistance and comprehensive traits of AL46, and just as we designed, AL46 showed harmonious powdery mildew resistance with some key breeding traits. This study not only developed an ideal wheat germplasm resource but also provided a successful example for pyramiding breeding, which could be a promising direction for wheat improvement in the future.

Characterization and identification of the powdery mildew resistance gene in wheat breeding line ShiCG15-009

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is a serious fungal disease that critically threatens the yield and quality of wheat. Utilization of host resistance is the most effective and economical method to control this disease. In our study, a wheat breeding line ShiCG15-009, released from Hebei Province, was highly resistant to powdery mildew at all stages. To dissect its genetic basis, ShiCG15-009 was crossed with the susceptible cultivar Yannong 21 to produce F₁, F₂ and F_2:3 progenies. After genetic analysis, a single dominant gene, tentatively designated PmCG15-009, was proved to confer resistance to Bgt isolate E09. Further molecular markers analysis showed that PmCG15-009 was located on chromosome 2BL and flanked by markers XCINAU130 and XCINAU143 with the genetic distances 0.2 and 0.4 cM, respectively, corresponding to a physic interval of 705.14-723.48 Mb referred to the Chinese Spring reference genome sequence v2.1. PmCG15-009 was most likely a new gene differed from the documented Pm genes on chromosome 2BL since its different origin, genetic diversity, and physical position. To analyze and identify the candidate genes, six genes associated with disease resistance in the candidate interval were confirmed to be associated with PmCG15-009 via qRT-PCR analysis using the parents ShiCG15-009 and Yannong 21 and time-course analysis post-inoculation with Bgt isolate E09. To accelerate the transfer of PmCG15-009 using marker-assisted selection (MAS), 18 closely or co-segregated markers were evaluated and confirmed to be suitable for tracing PmCG15-009, when it was transferred into different wheat cultivars.

Identification of a Wheat Powdery Mildew Dominant Resistance Gene in the Pm5 Locus for High-Throughput Marker-Assisted Selection

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), poses a severe threat to wheat yield and quality worldwide. Rapid identification and the accurate transference of effective resistance genes are important to the development of resistant cultivars and the sustainable control of this disease. In the present study, the wheat line AL11 exhibited high levels of resistance to powdery mildew at both the seedling and adult plant stages. Genetic analysis of the AL11 × 'Shixin 733' mapping population revealed that its resistance was controlled by a single dominant gene, tentatively designated PmAL11. Using bulked segregant RNA-Seq and molecular marker analysis, PmAL11 was mapped to the Pm5 locus on chromosome 7B where it cosegregated with the functional marker Pm5e-KASP. Sequence alignment analysis revealed that the Pm5e-homologous sequence in AL11 was identical to the reported recessive gene Pm5e in wheat landrace 'Fuzhuang 30'. It appears that PmAL11 was most probably Pm5e, but it was mediated by a dominant inheritance pattern, so it should provide a valuable resistance resource for both genetic study and wheat breeding. To efficiently use and trace PmAL11 in breeding, a new kompetitive allele-specific PCR marker AL11-K2488 that cosegregated with this gene was developed and confirmed to be applicable in the different wheat backgrounds, thus promoting its use in the marker-assisted selection of PmAL11.

Marker-trait association analyses revealed major novel QTLs for grain yield and related traits in durum wheat

The growing global demand for wheat for food is rising due to the influence of population growth and climate change. The dissection of complex traits by employing a genome-wide association study (GWAS) allows the identification of DNA markers associated with complex traits to improve the productivity of crops. We used GWAS with 10,045 single nucleotide polymorphism (SNP) markers to search for genomic regions associated with grain yield and related traits based on diverse panels of Ethiopian durum wheat. In Ethiopia, multi-environment trials of the genotypes were carried out at five locations. The genotyping was conducted using the 25k Illumina Wheat SNP array to explore population structure, linkage disequilibrium (LD), and marker-trait associations (MTAs). For GWAS, the multi-locus Fixed and Random Model Circulating Probability Unification (FarmCPU) model was applied. Broad-sense heritability estimates were high, ranging from 0.63 (for grain yield) to 0.97 (for thousand-kernel weight). The population structure based on principal component analysis, and model-based cluster analysis revealed two genetically distinct clusters with limited admixtures. The LD among SNPs declined within the range of 2.02-10.04 Mbp with an average of 4.28 Mbp. The GWAS scan based on the mean performance of the genotypes across the environments identified 44 significant MTAs across the chromosomes. Twenty-six of these MTAs are novel, whereas the remaining 18 were previously reported and confirmed in this study. We also identified candidate genes for the novel loci potentially regulating the traits. Hence, this study highlights the significance of the Ethiopian durum wheat gene pool for improving durum wheat globally. Furthermore, a breeding strategy focusing on accumulating favorable alleles at these loci could improve durum wheat production in the East African highlands and elsewhere.

Mapping QTL for Phenological and Grain-Related Traits in a Mapping Population Derived from High-Zinc-Biofortified Wheat

Genomic regions governing days to heading (DH), days to maturity (DM), plant height (PH), thousand-kernel weight (TKW), and test weight (TW) were investigated in a set of 190 RILs derived from a cross between a widely cultivated wheat-variety, Kachu (DPW-621-50), and a high-zinc variety, Zinc-Shakti. The RIL population was genotyped using 909 DArTseq markers and phenotyped in three environments. The constructed genetic map had a total genetic length of 4665 cM, with an average marker density of 5.13 cM. A total of thirty-seven novel quantitative trait loci (QTL), including twelve for PH, six for DH, five for DM, eight for TKW and six for TW were identified. A set of 20 stable QTLs associated with the expression of DH, DM, PH, TKW, and TW were identified in two or more environments. Three novel pleiotropic genomic-regions harboring co-localized QTLs governing two or more traits were also identified. In silico analysis revealed that the DArTseq markers were located on important putative candidate genes such as MLO-like protein, Phytochrome, Zinc finger and RING-type, Cytochrome P450 and pentatricopeptide repeat, involved in the regulation of pollen maturity, the photoperiodic modulation of flowering-time, abiotic-stress tolerance, grain-filling duration, thousand-kernel weight, seed morphology, and plant growth and development. The identified novel QTLs, particularly stable and co-localized QTLs, will be validated to estimate their effects in different genetic backgrounds for subsequent use in marker-assisted selection (MAS).

Amand P, Bai G, Baker J, Baker S, Liu S. QTL mapping of yield components and kernel traits in wheat cultivars TAM 112 and Duster

In the Southern Great Plains, wheat cultivars have been selected for a combination of outstanding yield and drought tolerance as a long-term breeding goal. To understand the underlying genetic mechanisms, this study aimed to dissect the quantitative trait loci (QTL) associated with yield components and kernel traits in two wheat cultivars `TAM 112' and `Duster' under both irrigated and dryland environments. A set of 182 recombined inbred lines (RIL) derived from the cross of TAM 112/Duster were planted in 13 diverse environments for evaluation of 18 yield and kernel related traits. High-density genetic linkage map was constructed using 5,081 single nucleotide polymorphisms (SNPs) from genotyping-by-sequencing (GBS). QTL mapping analysis detected 134 QTL regions on all 21 wheat chromosomes, including 30 pleiotropic QTL regions and 21 consistent QTL regions, with 10 QTL regions in common. Three major pleiotropic QTL on the short arms of chromosomes 2B (57.5 - 61.6 Mbps), 2D (37.1 - 38.7 Mbps), and 7D (66.0 - 69.2 Mbps) colocalized with genes Ppd-B1, Ppd-D1, and FT-D1, respectively. And four consistent QTL associated with kernel length (KLEN), thousand kernel weight (TKW), plot grain yield (YLD), and kernel spike^-1 (KPS) (Qklen.tamu.1A.325, Qtkw.tamu.2B.137, Qyld.tamu.2D.3, and Qkps.tamu.6A.113) explained more than 5% of the phenotypic variation. QTL Qklen.tamu.1A.325 is a novel QTL with consistent effects under all tested environments. Marker haplotype analysis indicated the QTL combinations significantly increased yield and kernel traits. QTL and the linked markers identified in this study will facilitate future marker-assisted selection (MAS) for pyramiding the favorable alleles and QTL map-based cloning.

Identification and validation of plant height, spike length and spike compactness loci in common wheat (Triticum aestivum L.)

BACKGROUND

Plant height (PH), spike length (SL) and spike compactness (SCN) are important agronomic traits in wheat due to their strong correlations with lodging and yield. Thus, dissection of their genetic basis is essential for the improvement of plant architecture and yield potential in wheat breeding. The objective of this study was to map quantitative trait loci (QTL) for PH, SL and SCN in a recombinant inbred line (RIL) population derived from the cross 'PuBing3228 × Gao8901' (PG-RIL) and to evaluate the potential values of these QTL to improve yield.

RESULTS

In the current study, Five, six and ten stable QTL for PH, SL, and SCN, respectively, were identified in at least two individual environments. Five major QTL QPh.cas-5A.3, QPh.cas-6A, QSl.cas-6B.2, QScn.cas-2B.2 and QScn.cas-6B explained 5.58-25.68% of the phenotypic variation. Notably, two, three and three novel stable QTL for PH, SL and SCN were identified in this study, which could provide further insights into the genetic factors that shape PH and spike morphology in wheat. Conditional QTL analysis revealed that QTL for SCN were mainly affected by SL. Moreover, a Kompetitive Allele Specific PCR (KASP) marker tightly linked to stable major QTL QPh.cas-5A.3 was developed and verified using the PG-RIL population and a natural population.

CONCLUSIONS

Twenty-one stable QTL related to PH, SL, and SCN were identified. These stable QTL and the user-friendly marker KASP8750 will facilitate future studies involving positional cloning and marker-assisted selection in breeding.

High-Density Linkage Mapping of Agronomic Trait QTLs in Wheat under Water Deficit Condition using Genotyping by Sequencing (GBS)

Improvement of grain yield is the ultimate goal for wheat breeding under water-limited environments. In the present study, a high-density linkage map was developed by using genotyping-by-sequencing (GBS) of a recombinant inbred line (RIL) population derived from the cross between Iranian landrace #49 and cultivar Yecora Rojo. The population was evaluated in three locations in Iran during two years under irrigated and water deficit conditions for the agronomic traits grain yield (GY), plant height (PH), spike number per square meter (SM), 1000 kernel weight (TKW), grain number per spike (GNS), spike length (SL), biomass (BIO) and harvest index (HI). A linkage map was constructed using 5831 SNPs assigned to 21 chromosomes, spanning 3642.14 cM of the hexaploid wheat genome with an average marker density of 0.62 (markers/cM). In total, 85 QTLs were identified on 19 chromosomes (all except 5D and 6D) explaining 6.06–19.25% of the traits phenotypic variance. We could identify 20 novel QTLs explaining 8.87–19.18% of phenotypic variance on chromosomes 1A, 1B, 1D, 2B, 3A, 3B, 6A, 6B and 7A. For 35 out of 85 mapped QTLs functionally annotated genes were identified which could be related to a potential role in drought stress.

Identification of the powdery mildew resistance gene in wheat breeding line Yannong 99102-06188 via bulked segregant exome capture sequencing

Powdery mildew of wheat (Triticum aestivum), caused by Blumeria graminis f.sp. tritici (Bgt), is a destructive disease that seriously threatens the yield and quality of its host. Identifying resistance genes is the most attractive and effective strategy for developing disease-resistant cultivars and controlling this disease. In this study, a wheat breeding line Yannong 99102-06188 (YN99102), an elite derivative line from the same breeding process as the famous wheat cultivar Yannong 999, showed high resistance to powdery mildew at the whole growth stages. Genetic analysis was carried out using Bgt isolate E09 and a population of YN99102 crossed with a susceptible parent Jinhe 13-205 (JH13-205). The result indicated that a single recessive gene, tentatively designated pmYN99102, conferred seedling resistance to the Bgt isolate E09. Using bulked segregant exome capture sequencing (BSE-Seq), pmYN99102 was physically located to a ~33.7 Mb (691.0-724.7 Mb) interval on the chromosome arm 2BL, and this interval was further locked in a 1.5 cM genetic interval using molecular markers, which was aligned to a 9.0 Mb physical interval (699.2-708.2 Mb). Based on the analysis of physical location, origin, resistant spectrum, and inherited pattern, pmYN99102 differed from those of the reported powdery mildew (Pm) resistance genes on 2BL, suggesting pmYN99102 is most likely a new Pm gene/allele in the targeted interval. To transfer pmYN99102 to different genetic backgrounds using marker-assisted selection (MAS), 18 closely linked markers were tested for their availability in different genetic backgrounds for MAS, and all markers expect for YTU103-97 can be used in MAS for tracking pmYN99102 when it transferred into those susceptible cultivars.

Fine mapping and characterization of a major QTL for grain weight on wheat chromosome arm 5DL

We fine mapped QTL QTKW.caas-5DL for thousand kernel weight in wheat, predicted candidate genes and developed a breeding-applicable marker. Thousand kernel weight (TKW) is an important yield component trait in wheat, and identification of the underlying genetic loci is helpful for yield improvement. We previously identified a stable quantitative trait locus (QTL) QTKW.caas-5DL for TKW in a Doumai/Shi4185 recombinant inbred line (RIL) population. Here we performed fine mapping of QTKW.caas-5DL using secondary populations derived from 15 heterozygous recombinants and delimited the QTL to an approximate 3.9 Mb physical interval from 409.9 to 413.8 Mb according to the Chinese Spring (CS) reference genome. Analysis of genomic synteny showed that annotated genes in the physical interval had high collinearity among CS and eight other wheat genomes. Seven genes with sequence variation and/or differential expression between parents were predicted as candidates for QTKW.caas-5DL based on whole-genome resequencing and transcriptome assays. A kompetitive allele-specific PCR (KASP) marker for QTKW.caas-5DL was developed, and genotyping confirmed a significant association with TKW but not with other yield component traits in a panel of elite wheat cultivars. The superior allele of QTKW.caas-5DL was frequent in a panel of cultivars, suggesting that it had undergone positive selection. These findings not only lay a foundation for map-based cloning of QTKW.caas-5DL but also provide an efficient tool for marker-assisted selection.

Genome-wide association study for grain yield and component traits in bread wheat (Triticum aestivum L.)

Genomic regions governing days to heading (DH), grain filling duration (GFD), grain number per spike (GNPS), grain weight per spike (GWPS), plant height (PH), and grain yield (GY) were investigated in a set of 280 diverse bread wheat genotypes. The genome-wide association studies (GWAS) panel was genotyped using a 35K Axiom Array and phenotyped in five environments. The GWAS analysis showed a total of 27 Bonferroni-corrected marker-trait associations (MTAs) on 15 chromosomes representing all three wheat subgenomes. The GFD showed the highest MTAs (8), followed by GWPS (7), GY (4), GNPS (3), PH (3), and DH (2). Furthermore, 20 MTAs were identified with more than 10% phenotypic variation. A total of five stable MTAs (AX-95024590, AX-94425015, AX-95210025 AX-94539354, and AX-94978133) were identified in more than one environment and associated with the expression of DH, GFD, GNPS, and GY. Similarly, two novel pleiotropic genomic regions with associated MTAs i.e. AX-94978133 (4D) and AX-94539354 (6A) harboring co-localized QTLs governing two or more traits were also identified. In silico analysis revealed that the SNPs were located on important putative candidate genes such as F-box-like domain superfamily, Lateral organ boundaries, LOB, Thioredoxin-like superfamily Glutathione S-transferase, RNA-binding domain superfamily, UDP-glycosyltransferase family, Serine/threonine-protein kinase, Expansin, Patatin, Exocyst complex component Exo70, DUF1618 domain, Protein kinase domain involved in the regulation of grain size, grain number, growth and development, grain filling duration, and abiotic stress tolerance. The identified novel MTAs will be validated to estimate their effects in different genetic backgrounds for subsequent use in marker-assisted selection (MAS).

High-resolution detection of quantitative trait loci for seven important yield-related traits in wheat (Triticum aestivum L.) using a high-density SLAF-seq genetic map

BACKGROUND

Yield-related traits including thousand grain weight (TGW), grain number per spike (GNS), grain width (GW), grain length (GL), plant height (PH), spike length (SL), and spikelet number per spike (SNS) are greatly associated with grain yield of wheat (Triticum aestivum L.). To detect quantitative trait loci (QTL) associated with them, 193 recombinant inbred lines derived from two elite winter wheat varieties Chuanmai42 and Chuanmai39 were employed to perform QTL mapping in six/eight environments.

RESULTS

A total of 30 QTLs on chromosomes 1A, 1B, 1D, 2A, 2B, 2D, 3A, 4A, 5A, 5B, 6A, 6D, 7A, 7B and 7D were identified. Among them, six major QTLs QTgw.cib-6A.1, QTgw.cib-6A.2, QGw.cib-6A, QGl.cib-3A, QGl.cib-6A, and QSl.cib-2D explaining 5.96-23.75% of the phenotypic variance were detected in multi-environments and showed strong and stable effects on corresponding traits. Three QTL clusters on chromosomes 2D and 6A containing 10 QTLs were also detected, which showed significant pleiotropic effects on multiple traits. Additionally, three Kompetitive Allele Specific PCR (KASP) markers linked with five of these major QTLs were developed. Candidate genes of QTgw.cib-6A.1/QGl.cib-6A and QGl.cib-3A were analyzed based on the spatiotemporal expression patterns, gene annotation, and orthologous search.

CONCLUSIONS

Six major QTLs for TGW, GL, GW and SL were detected. Three KASP markers linked with five of these major QTLs were developed. These QTLs and KASP markers will be useful for elucidating the genetic architecture of grain yield and developing new wheat varieties with high and stable yield in wheat.

A single nucleotide deletion in the third exon of FT-D1 increases the spikelet number and delays heading date in wheat (Triticum aestivum L.)

The spikelet number and heading date are two crucial and correlated traits for yield in wheat. Here, a quantitative trait locus (QTL) analysis was conducted in F₈ recombinant inbred lines (RILs) derived from crossing two common wheats with different spikelet numbers. A total of 15 stable QTL influencing total spikelet number (TSN) and heading date (HD) were detected. Notably, FT-D1, a well-known flowering time gene in wheat, was located within the finely mapped interval of a major QTL on 7DS (QTsn/Hd.cau-7D). A causal indel of one G in the third exon of FT-D1 was significantly associated with total spikelet number and heading date. Consistently, CRISPR/Cas9 mutant lines with homozygous mutations in FT-D1 displayed an increase in total spikelet number and heading date when compared with wild type. Moreover, one simple and robust marker developed according to the polymorphic site of FT-D1 revealed that this one G indel had been preferentially selected to adapt to different environments. Collectively, these data provide further insights into the genetic basis of spikelet number and heading date, and the diagnostic marker of FT-D1 will be useful for marker-assisted pyramiding in wheat breeding.

Identification and characterization of a natural polymorphism in FT-A2 associated with increased number of grains per spike in wheat

We discovered a natural FT-A2 allele that increases grain number per spike in both pasta and bread wheat with limited effect on heading time. Increases in wheat grain yield are necessary to meet future global food demands. A previous study showed that loss-of-function mutations in FLOWERING LOCUS T2 (FT2) increase spikelet number per spike (SNS), an important grain yield component. However, these mutations were also associated with reduced fertility, offsetting the beneficial effect of the increases in SNS on grain number. Here, we report a natural mutation resulting in an aspartic acid to alanine change at position 10 (D10A) associated with significant increases in SNS and no negative effects on fertility. Using a high-density genetic map, we delimited the SNS candidate region to a 5.2-Mb region on chromosome 3AS including 28 genes. Among them, only FT-A2 showed a non-synonymous polymorphism (D10A) present in two different populations segregating for the SNS QTL on chromosome arm 3AS. These results, together with the known effect of the ft-A2 mutations on SNS, suggest that variation in FT-A2 is the most likely cause of the observed differences in SNS. We validated the positive effects of the A10 allele on SNS, grain number, and grain yield per spike in near-isogenic tetraploid wheat lines and in an hexaploid winter wheat population. The A10 allele is present at very low frequency in durum wheat and at much higher frequency in hexaploid wheat, particularly in winter and fall-planted spring varieties. These results suggest that the FT-A2 A10 allele may be particularly useful for improving grain yield in durum wheat and fall-planted common wheat varieties.

Analysis of Genetic Regions Related to Field Grain Number per Spike From Chinese Wheat Founder Parent Linfen 5064

Wheat founder parents have been important in the development of new wheat cultivars. Understanding the effects of specific genome regions on yield-related traits in founder variety derivatives can enable more efficient use of these genetic resources through molecular breeding. In this study, the genetic regions related to field grain number per spike (GNS) from the founder parent Linfen 5064 were analyzed using a doubled haploid (DH) population developed from a cross between Linfen 5064 and Nongda 3338. Quantitative trait loci (QTL) for five spike-related traits over nine experimental locations/years were identified, namely, total spikelet number per spike (TSS), base sterile spikelet number per spike (BSSS), top sterile spikelet number per spike (TSSS), fertile spikelet number per spike (FSS), and GNS. A total of 13 stable QTL explaining 3.91-19.51% of the phenotypic variation were found. The effect of six of these QTL, Qtss.saw-2B.1, Qtss.saw-2B.2, Qtss.saw-3B, Qfss.saw-2B.2, Qbsss.saw-5A.1, and Qgns.saw-1A, were verified by another DH population (Linfen 5064/Jinmai 47), which showed extreme significance (P < 0.05) in more than three environments. No homologs of reported grain number-related from grass species were found in the physical regions of Qtss.saw-2B.1 and Qtss.saw-3B, that indicating both of them are novel QTL, or possess novel-related genes. The positive alleles of Qtss.saw-2B.2 from Linfen 5064 have the larger effect on TSS (3.30%, 0.62) and have 66.89% in Chinese cultivars under long-term artificial selection. This study revealed three key regions for GNS in Linfen 5064 and provides insights into molecular marker-assisted breeding.

Spatial and Temporal Transcriptomic Heredity and Asymmetry in an Artificially Constructed Allotetraploid Wheat (AADD)

Polyploidy, or whole-genome duplication (WGD), often induces dramatic changes in gene expression due to "transcriptome shock. " However, questions remain about how allopolyploidy (the merging of multiple nuclear genomes in the same nucleus) affects gene expression within and across multiple tissues and developmental stages during the initial foundation of allopolyploid plants. Here, we systematically investigated the immediate effect of allopolyploidy on gene expression variation in an artificial allopolyploidy system consisting of a constructed allotetraploid wheat (AADD genome, accession AT2) and its diploid progenitors Triticum urartu and Aegilops tauschii. We performed comprehensive RNA sequencing of 81 samples from different genotypes, tissues, and developmental stages. First, we found that intrinsic interspecific differences between the diploid parents played a major role in establishing the expression architecture of the allopolyploid. Nonetheless, allopolyploidy per se also induced dramatic and asymmetric patterns of differential gene expression between the subgenomes, and genes from the D subgenome exhibited a more drastic response. Second, analysis of homoeolog expression bias (HEB) revealed that the D subgenome exhibited significant expression bias and that de novo-generated HEB was attributed mainly to asymmetrical differential gene expression. Homoeolog-specific expression (HSE) analyses showed that the cis-only regulatory pattern was predominant in AT2, reflecting significant divergence between the parents. Co-expression network analysis revealed that homoeolog expression connectivity (HEC) was significantly correlated with sequence divergence in cis elements between subgenomes. Interestingly, allopolyploidy-induced reconstruction of network modules was also associated with different HSE patterns. Finally, a transcriptome atlas of spike development demonstrated that the phenotypic similarity of AT2 to T. urartu may be attributed to the combination of relatively stable expression of A-subgenome genes and drastic downregulation of their D-subgenome homoeologs. These findings provide a broad, multidimensional characterization of allopolyploidy-induced transcriptomic responses and suggest that allopolyploidy can have immediate and complex regulatory effects on the expression of nuclear genes.

Novel Quantitative Trait Loci for Grain Cadmium Content Identified in Hard White Spring Wheat

Cadmium (Cd) is a heavy metal that can cause a variety of adverse effects on human health, including cancer. Wheat comprises approximately 20% of the human diet worldwide; therefore, reducing the concentrations of Cd in wheat grain will have significant impacts on the intake of Cd in food products. The tests for measuring the Cd content in grain are costly, and the content is affected significantly by soil pH. To facilitate breeding for low Cd content, this study sought to identify quantitative trait loci (QTL) and associated molecular markers that can be used in molecular breeding. One spring wheat population of 181 doubled haploid lines (DHLs), which was derived from a cross between two hard white spring wheat cultivars "UI Platinum" (UIP) and "LCS Star" (LCS), was assessed for the Cd content in grain in multiple field trials in Southeast Idaho, United States. Three major QTL regions, namely, QCd.uia2-5B, QCd.uia2-7B, and QCd.uia2-7D, were identified on chromosomes 5B, 7B, and 7D, respectively. All genes in these three QTL regions were identified from the NCBI database. However, three genes related to the uptake and transport of Cd were used in the candidate gene analysis. The sequences of TraesCS5B02G388000 (TaHMA3) in the QCd.uia2-5B region and TraesCS7B02G320900 (TaHMA2) and TraesCS7B02G322900 (TaMSRMK3) in the QCd.uia2-7B region were compared between UIP and LCS. TaHMA2 on 7B is proposed for the first time as a candidate gene for grain Cd content in wheat. A KASP marker associated with this gene was developed and it will be further validated in near-isogenic lines via a gene-editing system in future studies.

Physical mapping of QTL associated with agronomic and end-use quality traits in spring wheat under conventional and organic management systems

Using phenotypic data of four biparental spring wheat populations evaluated at multiple environments under two management systems, we discovered 152 QTL and 22 QTL hotspots, of which two QTL accounted for up to 37% and 58% of the phenotypic variance, consistently detected in all environments, and fell within genomic regions harboring known genes. Identification of the physical positions of quantitative trait loci (QTL) would be highly useful for developing functional markers and comparing QTL results across multiple independent studies. The objectives of the present study were to map and characterize QTL associated with nine agronomic and end-use quality traits (tillering ability, plant height, lodging, grain yield, grain protein content, thousand kernel weight, test weight, sedimentation volume, and falling number) in hard red spring wheat recombinant inbred lines (RILs) using the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.0 physical map. We evaluated a total of 698 RILs from four populations derived from crosses involving seven parents at 3-8 conventionally (high N) and organically (low N) managed field environments. Using the phenotypic data combined across all environments per management, and the physical map between 1058 and 6526 markers per population, we identified 152 QTL associated with the nine traits, of which 29 had moderate and 2 with major effects. Forty-nine of the 152 QTL mapped across 22 QTL hotspot regions with each region coincident to 2-6 traits. Some of the QTL hotspots were physically located close to known genes. QSv.dms-1A and QPht.dms-4B.1 individually explained up to 37% and 58% of the variation in sedimentation volume and plant height, respectively, and had very large LOD scores that varied from 19.0 to 35.7 and from 16.7 to 55.9, respectively. We consistently detected both QTL in the combined and all individual environments, laying solid ground for further characterization and possibly for cloning.