Pleiotropic QTL influencing spikelet number and heading date in common wheat (Triticum aestivum L.)

Identification and map-based cloning of an EMS-induced mutation in wheat gene TaSP1 related to spike architecture

KEY MESSAGE

A candidate gene TaSP1 related to spike shape was cloned, and the gene-specific marker was developed to efficiently track the superior haplotype in common wheat. Spike shape, an important factor that affects wheat grain yield, is mainly defined by spike length (SPL), spikelet number (SPN), and compactness. Zhoumai32 mutant 1160 (ZM1160), a mutant obtained from ethyl methane sulfonate (EMS) treatment of hexaploid wheat variety Zhoumai32, was used to identify and clone the candidate gene that conditioned the spike shape. Genetic analysis of an F2 population derived from a cross of ZM1160 and Bainong207 suggested that the compact spike shape in ZM1160 was controlled by a single recessive gene, and therefore, the mutated gene was designated as Tasp1. With polymorphic markers identified through bulked segregant analysis (BSA), the gene was mapped to a 2.65-cM interval flanked by markers YZU0852 and MIS46239 on chromosome 7D, corresponding to a 0.42-Mb physical interval of Chinese spring (CS) reference sequences (RefSeq v1.0). To fine map TaSP1, 15 and seven recombinants were, respectively, screened from 1599 and 1903 F3 plants derived from the heterozygous F2 plants. Finally, TaSP1 was delimited to a 21.9 Kb (4,870,562 to 4,892,493 bp) Xmis48123-Xmis48104 interval. Only one high-confidence gene TraesCS7D02G010200 was annotated in this region, which encodes an unknown protein with a putative vWA domain. Quantitative reverse transcription PCR (qRT-PCR) analysis showed that TraesCS7D02G010200 was mainly expressed in the spike. Haplotype analysis of 655 wheat cultivars using the candidate gene-specific marker Xg010200p2 identified a superior haplotype TaSP1b with longer spike and more spikelet number. TaSP1 is beneficial to the improvement in wheat spike shape.

Genetic Analysis of Adaptive Traits in Spring Wheat in Northeast China

The dissection of the genetic architecture and the detection of the loci for adaptive traits are important for marker-assisted selection (MAS) for breeding. A spring wheat diversity panel with 251 cultivars, mainly from China, was obtained to conduct a genome-wide association study (GWAS) to detect the new loci, including the heading date (HD), maturating date (MD), plant height (PH), and lodging resistance (LR). In total, 41 loci existing in all 21 chromosomes, except for 4A and 6B, were identified, and each explained 4.3-18.9% of the phenotypic variations existing in two or more environments. Of these, 13 loci are overlapped with the known genes or quantitative trait loci (QTLs), whereas the other 28 are likely to be novel. The 1A locus (296.9-297.7 Mb) is a multi-effect locus for LR and PH, whereas the locus on chromosome 6D (464.5-471.0 Mb) affects both the HD and MD. Furthermore, four candidate genes for adaptive traits were identified, involved in cell division, signal transduction, and plant development. Additionally, two competitive, allele-specific PCR (KASP) markers, Kasp_2D_PH for PH and Kasp_6D_HD for HD, were developed and validated in another 162 spring wheat accessions. Our study uncovered the genetic basis of adaptive traits and provided the associated SNPs and varieties with more favorable alleles for wheat MAS breeding.

A High-Resolution Linkage Map Construction and QTL Analysis for Morphological Traits in Anthurium (Anthurium andraeanum Linden)

Anthurium andraeanum Linden is a prominent ornamental plant belonging to the family Araceae and is cultivated worldwide. The morphology characteristics are crucial because they significantly impact ornamental values, commercial properties, and the efficiency of space utilization in production. However, only a few related investigations have been conducted in anthurium to date. In this study, an F1 genetic segregation population containing 160 progenies was generated through hybridization between potted and cut anthurium varieties. Fifteen morphological traits were assessed and revealed substantial levels of genetic variation and widespread positive correlation. Based on specific length amplified fragment (SLAF) sequencing technology, 8171 single nucleotide polymorphism (SNP) markers were developed, and the high-density linkage map of 2202.27 cM in length distributed on 15 linkage groups was constructed successfully, with an average distance of 0.30 cM. Using the inclusive composite interval mapping (ICIM) method, 59 QTLs related to 15 key morphological traits were successfully identified, which explained phenotypic variance (PVE) ranging from 6.21% to 17.74%. Thirty-three of those associated with 13 traits were designated as major QTLs with PVE > 10%. These findings offer valuable insights into the genetic basis of quantitative traits and are beneficial for molecular marker-assisted selection (MAS) in anthurium breeding.

Genome-wide association analysis of time to heading and maturity in bread wheat using 55K microarrays

To investigate the genetic mechanisms underlying the reproductive traits (time to flowering and maturity) in wheat and identify candidate genes associated, a phenotypic analysis was conducted on 239 wheat accessions (lines) from around the world. A genome-wide association study (GWAS) of wheat heading and maturity phases was performed using the MLM (Q+K) model in the TASSLE software, combined with the Wheat 55K SNP array. The results revealed significant phenotypic variation in heading and maturity among the wheat accessions across different years, with coefficients of variation ranging from 0.96% to 1.97%. The phenotypic data from different years exhibited excellent correlation, with a genome-wide linkage disequilibrium (LD) attenuation distance of 3 Mb. Population structure analysis, evolutionary tree analysis, and principal component analysis indicated that the 239 wheat accessions formed a relatively homogeneous natural population, which could be divided into three subgroups. The GWAS results identified a total of 293 SNP marker loci that were significantly associated with wheat heading and maturity stages (P ≤ 0.001) in different environments. Among them, nine stable SNP marker loci were consistently detected in multiple environments. These marker loci were distributed on wheat chromosomes 1A、1B、2D、3A、5B、6D and 7A. Each individual locus explained 4.03%-16.06% of the phenotypic variation. Furthermore, through careful analysis of the associated loci with large phenotypic effect values and stable inheritance, a total of nine candidate genes related to wheat heading and maturity stages were identified. These findings have implications for molecular marker-assisted selection breeding programs targeting specific wheat traits at the heading and maturity stages. In summary, this study conducted a comprehensive GWAS of wheat heading and maturity phases, revealing significant associations between genetic markers and key developmental stages in wheat. The identification of candidate genes and marker loci provides valuable information for further studies on wheat breeding and genetic improvement targeted at enhancing heading and maturity traits.

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.

Deciphering spike architecture formation towards yield improvement in wheat

Wheat is the most widely grown crop globally, providing 20% of the daily consumed calories and protein content around the world. With the growing global population and frequent occurrence of extreme weather caused by climate change, ensuring adequate wheat production is essential for food security. The architecture of the inflorescence plays a crucial role in determining the grain number and size, which is a key trait for improving yield. Recent advances in wheat genomics and gene cloning techniques have improved our understanding of wheat spike development and its applications in breeding practices. Here, we summarize the genetic regulation network governing wheat spike formation, the strategies used for identifying and studying the key factors affecting spike architecture, and the progress made in breeding applications. Additionally, we highlight future directions that will aid in the regulatory mechanistic study of wheat spike determination and targeted breeding for grain yield improvement.

Two QTL regions for spike length showing pleiotropic effects on Fusarium head blight resistance and thousand-grain weight in bread wheat (Triticum aestivum L.)

Spike length (SL) plays an important role in the yield improvement of wheat and is significantly associated with other traits. Here, we used a recombinant inbred line (RIL) population derived from a cross between Yangmai 12 (YM12) and Yanzhan 1 (YZ1) to construct a genetic linkage map and identify quantitative trait loci (QTL) for SL. A total of 5 QTL were identified for SL, among which QSl.yaas-3A and QSl.yaas-5B are two novel QTL for SL. The YZ1 alleles at QSl.yaas-2D and QSl.yaas-5A, and the YM12 alleles at QSl.yaas-2A, QSl.yaas-3A, and QSl.yaas-5B conferred increasing SL effects. Two major QTL QSl.yaas-5A and QSl.yaas-5B explained 9.11-15.85% and 9.01-12.85% of the phenotypic variations, respectively. Moreover, the positive alleles of QSl.yaas-5A and QSl.yaas-5B could significantly increase Fusarium head blight (FHB) resistance (soil surface inoculation and spray inoculation were used) and thousand-grain weight (TGW) in the RIL population. Kompetitive allele-specific PCR (KASP) markers for QSl.yaas-5A and QSl.yaas-5B were developed and validated in an additional panel of 180 wheat cultivars/lines. The cultivars/lines harboring both the positive alleles of QSl.yaas-5A and QSl.yaas-5B accounted for only 28.33% of the validation populations and had the longest SL, best FHB resistance (using spray inoculation), and highest TGW. A total of 358 and 200 high-confidence annotated genes in QSl.yaas-5A and QSl.yaas-5B were identified, respectively. Some of the genes in these two regions were involved in cell development, disease resistance, and so on. The results of this study will provide a basis for directional breeding of longer SL, higher TGW, and better FHB resistance varieties and a solid foundation for fine-mapping QSl.yaas-5A and QSl.yaas-5B in future.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01427-8.

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.

An intron-located single nucleotide variation of TaGS5-3D is related to wheat grain size through accumulating intron retention transcripts

KEY MESSAGE

Thirty-three stable QTL for 13 yield-related traits across ten environments were identified in the PD34/MY47 RIL population, and a candidate gene TaGS5-3D in Qmt.nwafu.3D was preliminarily identified to affect grain-related traits through accumulation of specific transcripts. Dissecting the genetic basis of yield-related traits is pivotal for improvement of wheat yield potential. In this study, a recombinant inbred line (RIL) population genotyped by SNP markers was used to detect quantitative trait loci (QTL) related to yield-related traits in ten environments. A total of 102 QTL were detected, including 33 environmentally stable QTL and 69 putative QTL. Among them, Qmt.nwafu.3D was identified as a pleiotropic QTL in the physical interval of 149.77-154.11 Mb containing a potential candidate gene TaGS5-3D. An SNP (T > C) was detected in its ninth intron, and TaGS5-3D-C was validated as a superior allele associated with larger grains using a CAPS marker. Interestingly, we found that TaGS5-3D-C was closely related to significantly up-regulated expression of intron-retained transcript (TaGS5-3D-PD34.1), while TaGS5-3D-T was related to dominant expression of normal splicing transcript (TaGS5-3D-MY47.1). Our results indicated that alternative splicing associated with the SNP T/C could be involved in the regulation of grain-related traits, laying a foundation for the functional analysis of TaGS5-3D and its greater potential application in high-yield wheat breeding.

Identification and validation of new quantitative trait loci for spike-related traits in two RIL populations

Wheat (Triticum aestivum L.) is one of the most important cereal crops for ensuring food security worldwide. Identification of major quantitative trait loci (QTL) for spike-related traits is important for improvement of yield potential in wheat breeding. In this study, by using the wheat 55K single nucleotide polymorphism (SNP) array and diversity array technology (DArT), two recombinant inbred line populations derived from crosses avocet/chilero and avocet/huites were used to map QTL for kernel number per spike (KNS), total spikelet number per spike (TSS), fertile spikelet number per spike (FSS), and spike compactness (SC). Forty-two QTLs were identified on chromosomes 2A (4), 2B (3), 3A (2), 3B (7), 5A (11), 6A (4), 6B, and 7A (10), explaining 3.13-21.80% of the phenotypic variances. Twelve QTLs were detected in multi-environments on chromosomes 2A, 3B (2), 5A (4), 6A (3), 6B, and 7A, while four QTL clusters were detected on chromosomes 3A, 3B, 5A, and 7A. Two stable and new QTL clusters, QKns/Tss/Fss/SC.haust-5A and QKns/Tss/Fss.haust-7A, were detected in the physical intervals of 547.49-590.46 Mb and 511.54-516.15 Mb, accounting for 7.53-14.78% and 7.01-20.66% of the phenotypic variances, respectively. High-confidence annotated genes for QKns/Tss/Fss/SC.haust-5A and QKns/Tss/Fss.haust-7A were more highly expressed in spike development. The results provide new QTL and molecular markers for marker-assisted breeding in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01401-4.

Identification of novel putative alleles related to important agronomic traits of wheat using robust strategies in GWAS

Principal component analysis (PCA) is widely used in various genetics studies. In this study, the role of classical PCA (cPCA) and robust PCA (rPCA) was evaluated explicitly in genome-wide association studies (GWAS). We evaluated 294 wheat genotypes under well-watered and rain-fed, focusing on spike traits. First, we showed that some phenotypic and genotypic observations could be outliers based on cPCA and different rPCA algorithms (Proj, Grid, Hubert, and Locantore). Hubert's method provided a better approach to identifying outliers, which helped to understand the nature of these samples. These outliers led to the deviation of the heritability of traits from the actual value. Then, we performed GWAS with 36,000 single nucleotide polymorphisms (SNPs) based on the traditional approach and two robust strategies. In the conventional approach and using the first three components of cPCA as population structure, 184 and 139 marker-trait associations (MTAs) were identified for five traits in well-watered and rain-fed environments, respectively. In the first robust strategy and when rPCA was used as population structure in GWAS, we observed that the Hubert and Grid methods identified new MTAs, especially for yield and spike weight on chromosomes 7A and 6B. In the second strategy, we followed the classical and robust principal component-based GWAS, where the first two PCs obtained from phenotypic variables were used instead of traits. In the recent strategy, despite the similarity between the methods, some new MTAs were identified that can be considered pleiotropic. Hubert's method provided a better linear combination of traits because it had the most MTAs in common with the traditional approach. Newly identified SNPs, including rs19833 (5B) and rs48316 (2B), were annotated with important genes with vital biological processes and molecular functions. The approaches presented in this study can reduce the misleading GWAS results caused by the adverse effect of outlier observations.

Association mapping identifies loci and candidate genes for grain-related traits in spring wheat in response to heat stress

Heat stress is a limiting factor in wheat production along with global warming. Development of heat-tolerant wheat varieties and generation of suitable pre-breeding materials are the major goals in current wheat breeding programs. Our understanding on the genetic basis of thermotolerance remains sparse. In this study, we genotyped a collection of 211 core spring wheat accessions and conducted field trials to evaluate the grain-related traits under heat stress and non-stress conditions in two different locations for three consecutive years. Based on SNP datasets and grain-related traits, we performed genome-wide association study (GWAS) to detect stable loci related to thermotolerance. Thirty-three quantitative trait loci (QTL) were identified, nine of them are the same loci as previous studies, and 24 are potentially novel loci. Functional candidate genes at these QTL are predicted and proved to be relevant to heat stress and grain-related traits such as TaELF3-A1 (1A) for earliness per se (Eps), TaHSFA1-B1 (5B) influencing heat tolerance and TaVIN2-A1 (6A) for grain size. Functional markers of TaELF3-A1 were detected and converted to KASP markers, with their function and genetic diversity being analyzed in the natural populations. In addition, our results unveiled favor alleles controlling agronomic traits and/or heat stress tolerance. In summary, we provide insights into heritable correlation between yield and heat stress tolerance, which will accelerate the development of new cultivars with high and stable yield of wheat in the future.

Fine mapping and genetic analysis identified a C2H2-type zinc finger as a candidate gene for heading date regulation in wheat

KEY MESSAGE

A minor-effect QTL, Qhd.2AS, that affects heading date in wheat was mapped to a genomic interval of 1.70-Mb on 2AS, and gene analysis indicated that the C2H2-type zinc finger protein gene TraesCS2A02G181200 is the best candidate for Qhd.2AS. Heading date (HD) is a complex quantitative trait that determines the regional adaptability of cereal crops, and identifying the underlying genetic elements with minor effects on HD is important for improving wheat production in diverse environments. In this study, a minor QTL for HD that we named Qhd.2AS was detected on the short arm of chromosome 2A by Bulked Segregant Analysis and validated in a recombinant inbred population. Using a segregating population of 4894 individuals, Qhd.2AS was further delimited to an interval of 0.41 cM, corresponding to a genomic region spanning 1.70 Mb (from 138.87 to 140.57 Mb) that contains 16 high-confidence genes based on IWGSC RefSeq v1.0. Analyses of sequence variations and gene transcription indicated that TraesCS2A02G181200, which encodes a C2H2-type zinc finger protein, is the best candidate gene for Qhd.2AS that influences HD. Screening a TILLING mutant library identified two mutants with premature stop codons in TraesCS2A02G181200, both of which exhibited a delay in HD of 2-4 days. Additionally, variations in its putative regulatory sites were widely present in natural accession, and we also identified the allele which was positively selected during wheat breeding. Epistatic analyses indicated that Qhd.2AS-mediated HD variation is independent of VRN-B1 and environmental factors. Phenotypic investigation of homozygous recombinant inbred lines (RILs) and F2:3 families showed that Qhd.2AS has no negative effect on yield-related traits. These results provide important cues for refining HD and therefore improving yield in wheat breeding programs and will deepen our understanding of the genetic regulation of HD in cereal plants.

Contemplation on wheat vernalization

Vernalization is a period of low non-freezing temperatures, which provides the competence to flower. This mechanism ensures that plants sown before winter develop reproductive organs in more favourable conditions during spring. Such an evolutionary mechanism has evolved in both monocot and eudicot plants. Studies in monocots, represented by temperate cereals like wheat and barley, have identified and proposed the VERNALIZATION1 (VRN1) gene as a key player in the vernalization response. VRN1 belongs to MADS-box transcription factors and is expressed in the leaves and the apical meristem, where it subsequently promotes flowering. Despite substantial research advancement in the last two decades, there are still gaps in our understanding of the vernalization mechanism. Here we summarise the present knowledge of wheat vernalization. We discuss VRN1 allelic variation, review vernalization models, talk VRN1 copy number variation and devernalization phenomenon. Finally, we suggest possible future directions of the vernalization research in wheat.

Identification of quantitative trait loci of agronomic traits in bread wheat using a Pamyati Azieva × Paragon mapping population harvested in three regions of Kazakhstan

Background

Although genome-wide association studies (GWAS) are an increasingly informative tool in the mining of new quantitative trait loci (QTLs), a classical biparental mapping approach is still a powerful, widely used method to search the unique genetic factors associated with important agronomic traits in bread wheat.

Methods

In this study, a newly constructed mapping population of Pamyati Azieva (Russian Federation) × Paragon (UK), consisting of 94 recombinant inbred lines (RILs), was tested in three different regions of Kazakhstan with the purpose of QTL identification for key agronomic traits. The RILs were tested in 11 environments of two northern breeding stations (Petropavlovsk, North Kazakhstan region, and Shortandy, Aqmola region) and one southeastern station (Almalybak, Almaty region). The following eight agronomic traits were studied: heading days, seed maturation days, plant height, spike length, number of productive spikes, number of kernels per spike, thousand kernel weight, and yield per square meter. The 94 RILs of the PAxP cross were genotyped using Illumina's iSelect 20K single nucleotide polymorphism (SNP) array and resulted in the identification of 4595 polymorphic SNP markers.

Results

The application of the QTL Cartographer statistical package allowed the identification of 53 stable QTLs for the studied traits. A survey of published studies related to common wheat QTL identification suggested that 28 of those 53 QTLs were presumably novel genetic factors. The SNP markers for the identified QTLs of the analyzed agronomic traits of common wheat can be efficiently applied in ongoing breeding activities in the wheat breeding community using a marker-assisted selection approach.

Fine mapping of the reduced height gene Rht22 in tetraploid wheat landrace Jianyangailanmai (Triticum turgidum L.)

Rht22 was fine mapped in the interval of 0.53-1.48 Mb on 7AS, which reduces cell number of internode to cause semi-dwarfism in Jianyangailanmai. As a valuable germplasm resource for wheat genetic improvement, tetraploid wheat has several reduced height (Rht) and enhanced harvest index genes. Rht22, discovered in Jianyangailanmai (JAM, Triticum turgidum L., 2n = 4x = 28, AABB), significantly increases the spikelet number per spike, but its accurate chromosomal position is still unknown. In this study, a high-density genetic map was constructed using specific-length amplified fragment sequencing in an F₇ RIL_DJ population, which was derived from a cross between dwarf Polish wheat (T. polonicum L., 2n = 4x = 28, AABB) and JAM. Two plant height loci, Qph.sicau-4B and Qph.sicau-7A, were mapped on chromosomes 4BS and 7AS, respectively. Qph.sicau-7A was mapped to the 0.33-4.46 Mb interval on 7AS and likely represents the candidate region of Rht22. Fine mapping confirmed and narrowed Rht22 on chromosome arm 7AS between Xbag295.s53 and Xb295.191 in three different populations. The physical region ranged from 0.53 to 1.48 Mb and included 18 candidate genes. Transcriptome analysis of two pairs of near-isogenic lines revealed that 135 differentially expressed genes (DEGs) were associated with semi-dwarfism. Of these, the expression of 83 annotated DEGs involved in hormones synthesis and signal transduction, cell wall composition, DNA replication, microtubule and phragmoplast arrays was significantly down-regulated in the semi-dwarf line. Therefore, Rht22 causes semi-dwarfism in JAM by disrupting these cellular processes, which impairs cell proliferation and reduces internode cell number.

Overexpression of TaLBD16-4D alters plant architecture and heading date in transgenic wheat

Lateral organ boundaries domain (LBD) proteins, a class of plant-specific transcription factors with a special domain of lateral organ boundaries (LOB), play essential roles in plant growth and development. However, there is little known about the functions of these genes in wheat to date. Our previous study demonstrated that TaLBD16-4D is conducive to increasing lateral root number in wheat. In the present work, we further examined important agronomical traits of the aerial part of transgenic wheat overexpressing TaLBD16-4D. Interestingly, it was revealed that overexpressing TaLBD16-4D could lead to early heading and multiple alterations of plant architecture, including decreased plant height, increased flag leaf size and stem diameter, reduced spike length and tillering number, improved spike density and grain width, and decreased grain length. Moreover, auxin-responsive experiments demonstrated that the expression of TaLBD16-4D in wild-type (WT) wheat plants showed a significant upregulation through 2,4-D treatment. TaLBD16-4D-overexpression lines displayed a hyposensitivity to 2,4-D treatment and reduced shoot gravitropic response. The expressions of a set of auxin-responsive genes were markedly different between WT and transgenic plants. In addition, overexpressing TaLBD16-4D affected the transcript levels of flowering-related genes (TaGI, TaCO1, TaHd1, TaVRN1, TaVRN2, and TaFT1). Notably, the expression of TaGI, TaCO1, TaHd1, TaVRN1, and TaFT1 displayed significant upregulation under IAA treatment. Collectively, our observations indicated that overexpressing TaLBD16-4D could affect aerial architecture and heading time possibly though participating in the auxin pathway.

Wheat genomic study for genetic improvement of traits in China

Bread wheat (Triticum aestivum L.) is a major crop that feeds 40% of the world's population. Over the past several decades, advances in genomics have led to tremendous achievements in understanding the origin and domestication of wheat, and the genetic basis of agronomically important traits, which promote the breeding of elite varieties. In this review, we focus on progress that has been made in genomic research and genetic improvement of traits such as grain yield, end-use traits, flowering regulation, nutrient use efficiency, and biotic and abiotic stress responses, and various breeding strategies that contributed mainly by Chinese scientists. Functional genomic research in wheat is entering a new era with the availability of multiple reference wheat genome assemblies and the development of cutting-edge technologies such as precise genome editing tools, high-throughput phenotyping platforms, sequencing-based cloning strategies, high-efficiency genetic transformation systems, and speed-breeding facilities. These insights will further extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture in China and throughout the world.

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 validation of a QTL for spikelet number on chromosome arm 6BL of common wheat (Triticum aestivum L.)

To provide food security for a growing world population, it will be necessary to increase yields of staple crops such as wheat (Triticum aestivum L.). Yield is a complex, polygenic trait influenced by grain weight and number, which are negatively correlated with one another. Spikelet number is an important determinant of grain number, but allelic variants impacting its expression are often associated with heading date, constraining their use in wheat germplasm that must be adapted for specific environments. Identification and characterization of genetic variants affecting spikelet number will increase selection efficiency through their deployment in breeding programs. In this study, a quantitative trait locus (QTL) on chromosome arm 6BL for spikelet number was identified and validated using an association mapping panel, a recombinant inbred line population, and seven derived heterogeneous inbred families. The superior allele, QSn.csu-6Bb, was associated with an increase of 0.248 to 0.808 spikelets per spike across multiple environments that varied for mean spikelet number. Despite epistatic interactions between QSn.csu-6B and three other loci (WAPO-A1, VRN-D3, and PPD-B1), genotypes with a greater number of superior alleles at these loci consistently exhibit higher spikelet number. The frequency of superior alleles at these loci varies among winter wheat varieties adapted to different latitudes of the US Great Plains, revealing opportunities for breeders to select for increased spikelet number using simple molecular markers. This work lays the foundation for the positional cloning of the genetic variant underlying the QSn.csu-6B QTL to strengthen our understanding of spikelet number determination in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01288-7.

Identification and characterization of QTL for spike morphological traits, plant height and heading date derived from the D genome of natural and resynthetic allohexaploid wheat

QHd.cau-7D.1 for heading date was delimited into the physical interval of approximately 17.38 Mb harboring three CONSTANS-like zinc finger genes. Spike morphological traits, plant height and heading date play important roles in yield improvement of wheat. To reveal the genetic factors that controlling spike morphological traits, plant height and heading date on the D genome, we conducted analysis of quantitative traits locus (QTL) using 198 F_7:8 recombinant inbred lines (RILs) derived from a cross between the common wheat TAA10 and resynthesized allohexaploid wheat XX329 with similar AABB genomes. A total of 23 environmentally stable QTL on the D sub-genome for spike length (SL), fertile spikelet number per spike (FSN), sterile spikelet number per spike (SSN), total spikelet number per spike (TSN), spike compactness (SC), plant height (PHT) and heading date (HD) were detected, among which eight appeared to be novel QTL. Furthermore, QHd.cau-7D.1 and QPht.cau-7D.2 shared identical confidence interval and were delimited into the physical interval of approximately 17.38 Mb with 145 annotated genes, including three CONSTANS-like zinc finger genes (TraesCS7D02G209000, TraesCS7D02G213000 and TraesCS7D02G220300). This study will help elucidate the molecular mechanism of the seven traits (SL, FSN, SSN, TSN, SC, PHT and HD) and provide a potentially valuable resource for genetic improvement.

Combinations of Ghd7, Ghd8, and Hd1 determine strong heterosis of commercial rice hybrids in diverse ecological regions

Heterosis of grain yield is closely associated with heading date in crops. Gene combinations of the major heading date genes Ghd7, Ghd8, and Hd1 play important roles in enhancing grain yield and adaptation to ecological regions in rice. However, the predominant three-gene combinations for a specific ecological region remain unclear in both three-line and two-line hybrids. In this study, we sequenced these three genes of 50 cytoplasmic male sterile/maintainer lines, 31 photo-thermo-sensitive genic male sterile lines, and 109 restorer lines. Sequence analysis showed that hybrids carrying strong functional alleles of Ghd7 and Hd1 and non-functional Ghd8 are predominant in three-line hybrids and are recommended for rice production in the subtropics around 30°N/S. Hybrids carrying strong functional Ghd7 and Ghd8 and non-functional Hd1 are predominant in two-line hybrids and are recommended for low latitude areas around 23.5°N/S rich in photothermal resources. Hybrids carrying strong functional Ghd7 and Ghd8 and functional Hd1 were not identified in commercial hybrids in the middle and lower reaches of the Yangtze River, but they have high yield potential in tropical regions because they have the strongest photoperiod sensitivity. Based on these findings, two genic sterile lines, Xiangling 628S and C815S, whose hybrids often head very late, were diagnosed with these three genes, and Hd1 was targeted to be knocked out in Xiangling 628S and replaced with hd1 in C815S. The hybrids developed from both modified sterile lines in turn had appropriate heading dates and significantly improved grain yield. This study provides new insights for breeding design to develop hybrids for various regions.

Genome-wide association study for morphological, phenological, quality, and yield traits in einkorn (Triticum monococcum L. subsp. monococcum)

Einkorn (Triticum monococcum L. subsp. monococcum, 2n = 2× = 14, AmAm) is a diploid wheat whose cultivation was widespread in the Mediterranean and European area till the Bronze Age, before it was replaced by the more productive durum and bread wheats. Although scarcely cultivated nowadays, it has gained renewed interest due to its relevant nutritional properties and as source of genetic diversity for crop breeding. However, the molecular basis of many traits of interest in einkorn remain still unknown. A panel of 160 einkorn landraces, from different parts of the distribution area, was characterized for several phenotypic traits related to morphology, phenology, quality, and yield for 4 years in two locations. An approach based on co-linearity with the A genome of bread wheat, supported also by that with Triticum urartu genome, was exploited to perform association mapping, even without an einkorn anchored genome. The association mapping approach uncovered numerous marker-trait associations; for 37 of these, a physical position was inferred by homology with the bread wheat genome. Moreover, numerous associated regions were also assigned to the available T. monococcum contigs. Among the intervals detected in this work, three overlapped with regions previously described as involved in the same trait, while four other regions were localized in proximity of loci previously described and presumably refer to the same gene/QTL. The remaining associated regions identified in this work could represent a novel and useful starting point for breeding approaches to improve the investigated traits in this neglected species.

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

Four genomic regions on chromosomes 4A, 6A, 7B, and 7D were discovered, each with multiple tightly linked QTL (QTL clusters) associated with two to three yield components. The 7D QTL cluster was associated with grain yield, fertile spikelet number per spike, thousand kernel weight, and heading date. It was located in the flanking region of FT-D1, a homolog gene of Arabidopsis FLOWERING LOCUS T, a major gene that regulates wheat flowering. Genetic manipulation of yield components is an important approach to increase grain yield in wheat (Triticum aestivum). The present study used a mapping population comprised of 181 doubled haploid lines derived from two high-yielding spring wheat cultivars, UI Platinum and LCS Star. The two cultivars and the derived population were assessed for six traits in eight field trials primarily in Idaho in the USA. The six traits were grain yield, fertile spikelet number per spike, productive tiller number per unit area, thousand kernel weight, heading date, and plant height. Quantitative Trait Locus (QTL) analysis of the six traits was conducted using 14,236 single-nucleotide polymorphism (SNP) markers generated from the wheat 90 K SNP and the exome and promoter capture arrays. Of the 19 QTL detected, 14 were clustered in four chromosomal regions on 4A, 6A, 7B and 7D. Each of the four QTL clusters was associated with multiple yield component traits, and these traits were often negatively correlated with one another. As a result, additional QTL dissection studies are needed to optimize trade-offs among yield component traits for specific production environments. Kompetitive allele-specific PCR markers for the four QTL clusters were developed and assessed in an elite spring wheat panel of 170 lines, and eight of the 14 QTL were validated. The two parents contain complementary alleles for the four QTL clusters, suggesting the possibility of improving grain yield via genetic recombination of yield component loci.

Geographical distribution and adaptive variation of VRN-A3 alleles in worldwide polyploid wheat (Triticum spp.) species collection

The distribution of early flowering alleles of VRN-A3 was found to be biased to low latitudes, and these alleles may contribute to environmental adaptability to low latitudes in cultivated emmer wheat. In wheat (Triticum spp.), the flowering time is an important trait for successful seed production and yield by adapting to the regional environment. An early flowering allele of VRN-A3 with 7- and 25-bp insertions in the promoter region (Vrn-A3a-h1) has recently been reported from the analysis of an emmer wheat (Triticum turgidum L. ssp. dicoccum) accession, TN26. This early flowering allele of VRN-A3 might be associated with the regional adaptation of wheat. In this study, we elucidated its geographic distribution to assess the importance of the early flowering allele of VRN-A3 in worldwide wheat collection. From sequence analysis, we identified six VRN-A3 alleles with the 7- and 25-bp insertions, namely, Vrn-A3a-h2, Vrn-A3a-h3, Vrn-A3a-h4, Vrn-A3a-h5, Vrn-A3a-h6, and Vrn-A3c-h2 from wild emmer wheat, while we identified two VRN-A3 alleles with these insertions, Vrn-A3a-h2 and Vrn-A3c-h1 from cultivated tetraploid and hexaploid wheat species in addition to Vrn-A3a-h1. Among VRN-A3 alleles distributed in cultivated wheat, we found that Vrn-A3a-h2 promoted early heading, whereas Vrn-A3c-h1 did not affect heading time. Our analysis showed that the distribution of early flowering alleles of VRN-A3 dominated in cultivated emmer wheat in Ethiopia and India, which actually showed an early flowering phenotype. This implied that the early flowering alleles of VRN-A3 contribute to adaptability to a low-latitude environment in cultivated emmer wheat. We could not find durum (T. turgidum L. ssp. durum) and bread wheat (T. aestivum L. ssp. aestivum) accessions with these early flowering alleles. Our findings indicated that Vrn-A3a-h1 and Vrn-A3a-h2 were useful for breeding of early flowering cultivars in durum and bread wheat varieties.

Characterizing the oligogenic architecture of plant growth phenotypes informs genomic selection approaches in a common wheat population

BACKGROUND

Genetic variation in growth over the course of the season is a major source of grain yield variation in wheat, and for this reason variants controlling heading date and plant height are among the best-characterized in wheat genetics. While the major variants for these traits have been cloned, the importance of these variants in contributing to genetic variation for plant growth over time is not fully understood. Here we develop a biparental population segregating for major variants for both plant height and flowering time to characterize the genetic architecture of the traits and identify additional novel QTL.

RESULTS

We find that additive genetic variation for both traits is almost entirely associated with major and moderate-effect QTL, including four novel heading date QTL and four novel plant height QTL. FT2 and Vrn-A3 are proposed as candidate genes underlying QTL on chromosomes 3A and 7A, while Rht8 is mapped to chromosome 2D. These mapped QTL also underlie genetic variation in a longitudinal analysis of plant growth over time. The oligogenic architecture of these traits is further demonstrated by the superior trait prediction accuracy of QTL-based prediction models compared to polygenic genomic selection models.

CONCLUSIONS

In a population constructed from two modern wheat cultivars adapted to the southeast U.S., almost all additive genetic variation in plant growth traits is associated with known major variants or novel moderate-effect QTL. Major transgressive segregation was observed in this population despite the similar plant height and heading date characters of the parental lines. This segregation is being driven primarily by a small number of mapped QTL, instead of by many small-effect, undetected QTL. As most breeding populations in the southeast U.S. segregate for known QTL for these traits, genetic variation in plant height and heading date in these populations likely emerges from similar combinations of major and moderate effect QTL. We can make more accurate and cost-effective prediction models by targeted genotyping of key SNPs.

Natural variation in rice ascorbate peroxidase gene APX9 is associated with a yield-enhancing QTL cluster

We previously identified a cluster of yield-related quantitative trait loci (QTLs) including plant height in CR4379, a near-isogenic line from a cross between Oryza sativa spp. japonica cultivar 'Hwaseong' and the wild relative Oryza rufipogon. Map-based cloning and transgenic approaches revealed that APX9, which encodes an l-ascorbate peroxidase 4, is associated with this cluster. A 3 bp InDel was observed leading to the addition of a valine in Hwaseong compared with O. rufipogon. APX9-overexpressing transgenic plants in the Hwaseong background were taller than Hwaseong. Consistent with these results, APX9 T-DNA insertion mutants in the japonica cultivar Dongjin were shorter. These results confirm that APX9 is the causal gene for the QTL cluster. Sequence analysis of APX9 from 303 rice accessions revealed that the 3 bp InDel clearly differentiates japonica (APX9HS) and O. rufipogon (APX9OR) alleles. indica accessions shared both alleles, suggesting that APX9HS was introgressed into indica followed by crossing. The finding that O. rufipogon accessions with different origins carry APX9OR suggests that the 3 bp insertion was specifically selected in japonica during its domestication. Our findings demonstrate that APX9 acts as a major regulator of plant development by controlling a valuable suite of agronomically important traits in rice.

De Novo Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

Identification and Validation of a Novel Locus Controlling Spikelet Number in Bread Wheat (Triticum aestivum L.)

Spikelet number is an important target trait for wheat yield improvement. Thus, the identification and verification of novel quantitative trait locus (QTL)/genes controlling spikelet number are essential for dissecting the underlying molecular mechanisms and hence for improving grain yield. In the present study, we constructed a high-density genetic map for the Kechengmai1/Chuanmai42 doubled haploid (DH) population using 13,068 single-nucleotide polymorphism (SNP) markers from the Wheat 55K SNP array. A comparison between the genetic and physical maps indicated high consistence of the marker orders. Based on this genetic map, a total of 27 QTLs associated with total spikelet number per spike (TSN) and fertile spikelet number per spike (FSN) were detected on chromosomes 1B, 1D, 2B, 2D, 3D, 4A, 4D, 5A, 5B, 5D, 6A, 6B, and 7D in five environments. Among them, five QTLs on chromosome 2D, 3D, 5A, and 7D were detected in multiple environments and combined QTL analysis, explaining the phenotypic variance ranging from 3.64% to 23.28%. Particularly, QTsn/Fsn.cib-3D for TSN and FSN [phenotypic variation explained (PVE) = 5.97-23.28%, limit of detection (LOD) = 3.73-18.51] is probably a novel locus and located in a 4.5-cM interval on chromosome arm 3DL flanking by the markers AX-110914105 and AX-109429351. This QTL was further validated in other two populations with different genetic backgrounds using the closely linked Kompetitive Allele-Specific PCR (KASP) marker KASP_AX-110914105. The results indicated that QTsn/Fsn.cib-3D significantly increased the TSN (5.56-7.96%) and FSN (5.13-9.35%), which were significantly correlated with grain number per spike (GNS). We also preliminary analyzed the candidate genes within this locus by sequence similarity, spatial expression patterns, and collinearity analysis. These results provide solid foundation for future fine mapping and cloning of QTsn/Fsn.cib-3D. The developed and validated KASP markers could be utilized in molecular breeding aiming to increase the grain yield in wheat.

Linkage mapping identifies a non-synonymous mutation in FLOWERING LOCUS T (FT-B1) increasing spikelet number per spike

Total spikelet number per spike (TSN) is a major component of spike architecture in wheat (Triticumaestivum L.). A major and consistent quantitative trait locus (QTL) was discovered for TSN in a doubled haploid spring wheat population grown in the field over 4 years. The QTL on chromosome 7B explained up to 20.5% of phenotypic variance. In its physical interval (7B: 6.37–21.67 Mb), the gene FLOWERINGLOCUST (FT-B1) emerged as candidate for the observed effect. In one of the parental lines, FT-B1 carried a non-synonymous substitution on position 19 of the coding sequence. This mutation modifying an aspartic acid (D) into a histidine (H) occurred in a highly conserved position. The mutation was observed with a frequency of ca. 68% in a set of 135 hexaploid wheat varieties and landraces, while it was not found in other plant species. FT-B1 only showed a minor effect on heading and flowering time (FT) which were dominated by a major QTL on chromosome 5A caused by segregation of the vernalization gene VRN-A1. Individuals carrying the FT-B1 allele with amino acid histidine had, on average, a higher number of spikelets (15.1) than individuals with the aspartic acid allele (14.3) independent of their VRN-A1 allele. We show that the effect of TSN is not mainly related to flowering time; however, the duration of pre-anthesis phases may play a major role.

The 55K SNP-Based Exploration of QTLs for Spikelet Number Per Spike in a Tetraploid Wheat (Triticum turgidum L.) Population: Chinese Landrace "Ailanmai" × Wild Emmer

Spikelet number per spike (SNS) is the primary factor that determines wheat yield. Common wheat breeding reduces the genetic diversity among elite germplasm resources, leading to a detrimental effect on future wheat production. It is, therefore, necessary to explore new genetic resources for SNS to increase wheat yield. A tetraploid landrace "Ailanmai" × wild emmer wheat recombinant inbred line (RIL) population was used to construct a genetic map using a wheat 55K single- nucleotide polymorphism (SNP) array. The linkage map containing 1,150 bin markers with a total genetic distance of 2,411.8 cm was obtained. Based on the phenotypic data from the eight environments and best linear unbiased prediction (BLUP) values, five quantitative trait loci (QTLs) for SNS were identified, explaining 6.71-29.40% of the phenotypic variation. Two of them, QSns.sau-AM-2B.2 and QSns.sau-AM-3B.2, were detected as a major and novel QTL. Their effects were further validated in two additional F₂ populations using tightly linked kompetitive allele-specific PCR (KASP) markers. Potential candidate genes within the physical intervals of the corresponding QTLs were predicted to participate in inflorescence development and spikelet formation. Genetic associations between SNS and other agronomic traits were also detected and analyzed. This study demonstrates the feasibility of the wheat 55K SNP array developed for common wheat in the genetic mapping of tetraploid population and shows the potential application of wheat-related species in wheat improvement programs.

Genetic Characterization of Multiple Components Contributing to Fusarium Head Blight Resistance of FL62R1, a Canadian Bread Wheat Developed Using Systemic Breeding

Fusarium head blight (FHB) is a devastating fungal disease of small-grain cereals that results in severe yield and quality losses. FHB resistance is controlled by resistance components including incidence, field severity, visual rating index, Fusarium damaged kernels (FDKs), and the accumulation of the mycotoxin deoxynivalenol (DON). Resistance conferred by each of these components is partial and must be combined to achieve resistance sufficient to protect wheat from yield losses. In this study, two biparental mapping populations were analyzed in Canadian FHB nurseries and quantitative trait loci (QTL) mapped for the traits listed above. Nine genomic loci, on 2AS, 2BS, 3BS, 4AS, 4AL, 4BS, 5AS, 5AL, and 5BL, were enriched for the majority of the QTL controlling FHB resistance. The previously validated FHB resistance QTL on 3BS and 5AS affected resistance to severity, FDK, and DON in these populations. The remaining seven genomic loci colocalize with flowering time and/or plant height QTL. The QTL on 4B was a major contributor to all field resistance traits and plant height in the field. QTL on 4AL showed contrasting effects for FHB resistance between Eastern and Western Canada, indicating a local adapted resistance to FHB. In addition, we also found that the 2AS QTL contributed a major effect for DON, and the 2BS for FDK, while the 5AL conferred mainly effect for both FDK/DON. Results presented here provide insight into the genetic architecture underlying these resistant components and insight into how FHB resistance in wheat is controlled by a complex network of interactions between genes controlling flowering time, plant height, local adaption, and FHB resistance components.