Candidate Genes and Quantitative Trait Loci for Grain Yield and Seed Size in Durum Wheat

A bird's-eye view: exploration of the flavin-containing monooxygenase superfamily in common wheat

The Flavin Monooxygenase (FMO) gene superfamily in plants is involved in various processes most widely documented for its involvement in auxin biosynthesis, specialized metabolite biosynthesis, and plant microbial defense signaling. The roles of FMOs in defense signaling and disease resistance have recently come into focus as they may present opportunities to increase immune responses in plants including leading to systemic acquired resistance, but are not well characterized. We present a comprehensive catalogue of FMOs found in genomes across vascular plants and explore, in depth, 170 wheat TaFMO genes for sequence architecture, cis-acting regulatory elements, and changes due to Transposable Element insertions. A molecular phylogeny separates TaFMOs into three clades (A, B, and C) for which we further report gene duplication patterns, and differential rates of homoeologue expansion and retention among TaFMO subclades. We discuss Clade B TaFMOs where gene expansion is similarly seen in other cereal genomes. Transcriptome data from various studies point towards involvement of subclade B2 TaFMOs in disease responses against both biotrophic and necrotrophic pathogens, substantiated by promoter element analysis. We hypothesize that certain TaFMOs are responsive to both abiotic and biotic stresses, providing potential targets for enhancing disease resistance, plant yield and other important agronomic traits. Altogether, FMOs in wheat and other crop plants present an untapped resource to be exploited for improving the quality of crops.

A novel variation of TaGW2-6B increases grain weight without penalty in grain protein content in wheat (Triticum aestivum L.)

Yield and quality are two crucial breeding objects of wheat therein grain weight and grain protein content (GPC) are two key relevant factors correspondingly. Investigations of their genetic mechanisms represent special significance for breeding. In this study, 199 F2 plants and corresponding F2:3 families derived from Nongda3753 (ND3753) and its EMS-generated mutant 564 (M564) were used to investigate the genetic basis of larger grain and higher GPC of M564. QTL analysis identified a total of 33 environmentally stable QTLs related to thousand grain weight (TGW), grain area (GA), grain circle (GC), grain length (GL), grain width (GW), and GPC on chromosomes 1B, 2A, 2B, 4D, 6B, and 7D, respectively, among which QGw.cau-6B.1, QTgw.cau-6B.1, QGa.cau-6B.1, and QGc.cau-6B.1 shared overlap confidence interval on chromosome 6B. This interval contained the TaGW2 gene playing the same role as the QTLs, so TaGW2-6B was cloned and sequenced. Sequence alignment revealed two G/A SNPs between two parents, among which the SNP in the seventh exon led to a premature termination in M564. A KASP marker was developed based on the SNP, and single-marker analysis on biparental populations showed that the mutant allele could significantly increase GW and TGW, but had no effect on GPC. Distribution detection of the mutant allele through KASP marker genotyping and sequence alignment against databases ascertained that no materials harbored this allele within natural populations. This allele was subsequently introduced into three different varieties through molecular marker-assisted backcrossing, and it was revealed that the allele had a significant effect on simultaneously increasing GW, TGW, and even GPC in all of three backgrounds. Summing up the above, it could be concluded that a novel elite allele of TaGW2-6B was artificially created and might play an important role in wheat breeding for high yield and quality.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01455-y.

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.

QTL mapping for kernel-related traits in a durum wheat x T. dicoccum segregating population

Durum wheat breeding relies on grain yield improvement to meet its upcoming demand while coping with climate change. Kernel size and shape are the determinants of thousand kernel weight (TKW), which is a key component of grain yield, and the understanding of the genetic control behind these traits supports the progress in yield potential. The present study aimed to dissect the genetic network responsible for kernel size components (length, width, perimeter, and area) and kernel shape traits (width-to-length ratio and formcoefficient) as well as their relationships with kernel weight, plant height, and heading date in durum wheat. Quantitative Trait Locus (QTL) mapping was performed on a segregating population of 110 recombinant inbred lines, derived from a cross between the domesticated emmer wheat accession MG5323 and the durum wheat cv. Latino, evaluated in four different environments. A total of 24 QTLs stable across environments were found and further grouped in nine clusters on chromosomes 2A, 2B, 3A, 3B, 4B, 6B, and 7A. Among them, a QTL cluster on chromosome 4B was associated with kernel size traits and TKW, where the parental MG5323 contributed the favorable alleles, highlighting its potential to improve durum wheat germplasm. The physical positions of the clusters, defined by the projection on the T. durum reference genome, overlapped with already known genes (i.e., BIG GRAIN PROTEIN 1 on chromosome 4B). These results might provide genome-based guidance for the efficient exploitation of emmer wheat diversity in wheat breeding, possibly through yield-related molecular markers.

Demographic history and distinct selection signatures of two domestication genes in mungbean

Domestication is the long and complex process underlying the evolution of crops, in which artificial directional selection transformed wild progenitors into the desired form, affecting genomic variation and leaving traces of selection at targeted loci. However, whether genes controlling important domestication traits follow the same evolutionary pattern expected under the standard selective sweep model remains unclear. With whole-genome resequencing of mungbean (Vigna radiata), we investigated this issue by resolving its global demographic history and targeted dissection of the molecular footprints of genes underlying 2 key traits representing different stages of domestication. Mungbean originated in Asia, and the Southeast Asian wild population migrated to Australia about 50 thousand generations ago. Later in Asia, the cultivated form diverged from the wild progenitor. We identified the gene associated with the pod shattering resistance trait, VrMYB26a, with lower expression across cultivars and reduced polymorphism in the promoter region, reflecting a hard selective sweep. On the other hand, the stem determinacy trait was associated with VrDet1. We found that 2 ancient haplotypes of this gene have lower gene expression and exhibited intermediate frequencies in cultivars, consistent with selection favoring independent haplotypes in a soft selective sweep. In mungbean, contrasting signatures of selection were identified from the detailed dissection of 2 important domestication traits. The results suggest complex genetic architecture underlying the seemingly simple process of directional artificial selection and highlight the limitations of genome-scan methods relying on hard selective sweeps.

Identification of stable QTL controlling multiple yield components in a durum × cultivated emmer wheat population under field and greenhouse conditions

Crop yield gains are needed to keep pace with a growing global population and decreasing resources to produce food. Cultivated emmer wheat is a progenitor of durum wheat and a useful source of genetic variation for trait improvement in durum. Here, we evaluated a recombinant inbred line population derived from a cross between the North Dakota durum wheat variety Divide and the cultivated emmer wheat accession PI 272527 consisting of 219 lines. The population was evaluated in 3 field environments and 2 greenhouse experiments to identify quantitative trait locus associated with 11 yield-related traits that were expressed in a consistent manner over multiple environments. We identified 27 quantitative trait locus expressed in at least 2 field environments, 17 of which were also expressed under greenhouse conditions. Seven quantitative trait locus regions on chromosomes 1B, 2A, 2B, 3A, 3B, 6A, and 7B had pleiotropic effects on multiple yield-related traits. The previously cloned genes Q and FT-B1, which are known to be associated with development and morphology, were found to consistently be associated with multiple traits across environments. PI 272527 contributed beneficial alleles for quantitative trait locus associated with multiple traits, especially for seed morphology quantitative trait locus on chromosomes 1B, 2B, and 6A. Three recombinant inbred lines with increased grain size and weight compared to Divide were identified and demonstrated the potential for improvement of durum wheat through deployment of beneficial alleles from the cultivated emmer parent. The findings from this study provide knowledge regarding stable and robust quantitative trait locus that breeders can use for improving yield in durum wheat.

Genome-wide association studies and genomic selection assays made in a large sample of cacao (Theobroma cacao L.) germplasm reveal significant marker-trait associations and good predictive value for improving yield potential

A genome-wide association study (GWAS) was undertaken to unravel marker-trait associations (MTAs) between SNP markers and phenotypic traits. It involved a subset of 421 cacao accessions from the large and diverse collection conserved ex situ at the International Cocoa Genebank Trinidad. A Mixed Linear Model (MLM) in TASSEL was used for the GWAS and followed by confirmatory analyses using GAPIT FarmCPU. An average linkage disequilibrium (r2) of 0.10 at 5.2 Mb was found across several chromosomes. Seventeen significant (P ≤ 8.17 × 10-5 (-log10 (p) = 4.088)) MTAs of interest, including six that pertained to yield-related traits, were identified using TASSEL MLM. The latter accounted for 5 to 17% of the phenotypic variation expressed. The highly significant association (P ≤ 8.17 × 10-5) between seed length to width ratio and TcSNP 733 on chromosome 5 was verified with FarmCPU (P ≤ 1.12 × 10-8). Fourteen MTAs were common to both the TASSEL and FarmCPU models at P ≤ 0.003. The most significant yield-related MTAs involved seed number and seed length on chromosome 7 (P ≤ 1.15 × 10-14 and P ≤ 6.75 × 10-05, respectively) and seed number on chromosome 1 (P ≤ 2.38 × 10-05), based on the TASSEL MLM. It was noteworthy that seed length, seed length to width ratio and seed number were associated with markers at different loci, indicating their polygenic nature. Approximately 40 candidate genes that encode embryo and seed development, protein synthesis, carbohydrate transport and lipid biosynthesis and transport were identified in the flanking regions of the significantly associated SNPs and in linkage disequilibrium with them. A significant association of fruit surface anthocyanin intensity co-localised with MYB-related protein 308 on chromosome 4. Testing of a genomic selection approach revealed good predictive value (genomic estimated breeding values (GEBV)) for economic traits such as seed number (GEBV = 0.611), seed length (0.6199), seed width (0.5435), seed length to width ratio (0.5503), seed/cotyledon mass (0.6014) and ovule number (0.6325). The findings of this study could facilitate genomic selection and marker-assisted breeding of cacao thereby expediting improvement in the yield potential of cacao planting material.

Genome-Wide Association Study in Bread Wheat Identifies Genomic Regions Associated with Grain Yield and Quality under Contrasting Water Availability

The objectives of this study were to identify genetic loci in the bread wheat genome that would influence yield stability and quality under water stress, and to identify accessions that can be recommended for cultivation in dry and hot regions. We performed a genome-wide association study (GWAS) using a panel of 232 wheat accessions spanning diverse ecogeographic regions. Plants were evaluated in the Israeli Northern Negev, under two environments: water-limited (D; 250 mm) and well-watered (W; 450 mm) conditions; they were genotyped with ~71,500 SNPs derived from exome capture sequencing. Of the 14 phenotypic traits evaluated, 12 had significantly lower values under D compared to W conditions, while the values for two traits were higher under D. High heritability (H² = 0.5–0.9) was observed for grain yield, spike weight, number of grains per spike, peduncle length, and plant height. Days to heading and grain yield could be partitioned based on accession origins. GWAS identified 154 marker-trait associations (MTAs) for yield and quality-related traits, 82 under D and 72 under W, and identified potential candidate genes. We identified 24 accessions showing high and/or stable yields under D conditions that can be recommended for cultivation in regions under the threat of global climate change.

Quantitative trait loci identification and breeding value estimation of grain weight-related traits based on a new wheat 50K single nucleotide polymorphism array-derived genetic map

Mining novel and less utilized thousand grain weight (TGW) related genes are useful for improving wheat yield. In this study, a recombinant inbred line population from a cross between Zhongkemai 138 (ZKM138, high TGW) and Chuanmai 44 (CM44, low TGW) was used to construct a new Wheat 50K SNP array-derived genetic map that spanned 1,936.59 cM and contained 4, 139 markers. Based on this map, ninety-one quantitative trait loci (QTL) were detected for eight grain-related traits in six environments. Among 58 QTLs, whose superior alleles were contributed by ZKM138, QTgw.cib-6A was a noticeable major stable QTL and was also highlighted by bulked segregant analysis with RNA sequencing (BSR-Seq). It had a pyramiding effect on TGW enhancement but no significant trade-off effect on grain number per spike or tiller number, with two other QTLs (QTgw.cib-2A.2 and QTgw.cib-6D), possibly explaining the excellent grain performance of ZKM138. After comparison with known loci, QTgw.cib-6A was deduced to be a novel locus that differed from nearby TaGW2 and TaBT1. Seven simple sequence repeat (SSR) and thirty-nine kompetitive allele-specific PCR (KASP) markers were finally developed to narrow the candidate interval of QTgw.cib-6A to 4.1 Mb. Only six genes in this interval were regarded as the most likely candidate genes. QTgw.cib-6A was further validated in different genetic backgrounds and presented 88.6% transmissibility of the ZKM138-genotype and a 16.4% increase of TGW in ZKM138 derivatives. And the geographic pattern of this locus revealed that its superior allele is present in only 6.47% of 433 Chinese modern wheat varieties, indicating its potential contribution to further high-yield breeding.

A comprehensive meta-QTL analysis for yield-related traits of durum wheat (Triticum turgidum L. var. durum) grown under different water regimes

Abiotic stress strongly affects yield-related traits in durum wheat, in particular drought is one of the main environmental factors that have effect on grain yield and plant architecture. In order to obtain new genotypes well adapted to stress conditions, the highest number of desirable traits needs to be combined in the same genotype. In this context, hundreds of quantitative trait loci (QTL) have been identified for yield-related traits in different genetic backgrounds and environments. Meta-QTL (MQTL) analysis is a useful approach to combine data sets and for creating consensus positions for the QTL detected in independent studies for the reliability of their location and effects. MQTL analysis is a useful method to dissect the genetic architecture of complex traits, which provide an extensive allelic coverage, a higher mapping resolution and allow the identification of putative molecular markers useful for marker-assisted selection (MAS). In the present study, a complete and comprehensive MQTL analysis was carried out to identify genomic regions associated with grain-yield related traits in durum wheat under different water regimes. A total of 724 QTL on all 14 chromosomes (genomes A and B) were collected for the 19 yield-related traits selected, of which 468 were reported under rainfed conditions, and 256 under irrigated conditions. Out of the 590 QTL projected on the consensus map, 421 were grouped into 76 MQTL associated with yield components under both irrigated and rainfed conditions, 12 genomic regions containing stable MQTL on all chromosomes except 1A, 4A, 5A, and 6B. Candidate genes associated to MQTL were identified and an in-silico expression analysis was carried out for 15 genes selected among those that were differentially expressed under drought. These results can be used to increase durum wheat grain yields under different water regimes and to obtain new genotypes adapted to climate change.

Genetic dissection of grain architecture-related traits in a winter wheat population

BACKGROUND

The future productivity of wheat (T. aestivum L.) as the most grown crop worldwide is of utmost importance for global food security. Thousand kernel weight (TKW) in wheat is closely associated with grain architecture-related traits, e.g. kernel length (KL), kernel width (KW), kernel area (KA), kernel diameter ratio (KDR), and factor form density (FFD). Discovering the genetic architecture of natural variation in these traits, identifying QTL and candidate genes are the main aims of this study. Therefore, grain architecture-related traits in 261 worldwide winter accessions over three field-year experiments were evaluated.

RESULTS

Genome-wide association analysis using 90K SNP array in FarmCPU model revealed several interesting genomic regions including 17 significant SNPs passing false discovery rate threshold and strongly associated with the studied traits. Four of associated SNPs were physically located inside candidate genes within LD interval e.g. BobWhite_c5872_589 (602,710,399 bp) found to be inside TraesCS6A01G383800 (602,699,767-602,711,726 bp). Further analysis reveals the four novel candidate genes potentially involved in more than one grain architecture-related traits with a pleiotropic effects e.g. TraesCS6A01G383800 gene on 6A encoding oxidoreductase activity was associated with TKW and KA. The allelic variation at the associated SNPs showed significant differences betweeen the accessions carying the wild and mutated alleles e.g. accessions carying C allele of BobWhite_c5872_589, TraesCS6A01G383800 had significantly higher TKW than the accessions carying T allele. Interestingly, these genes were highly expressed in the grain-tissues, demonstrating their pivotal role in controlling the grain architecture.

CONCLUSIONS

These results are valuable for identifying regions associated with kernel weight and dimensions and potentially help breeders in improving kernel weight and architecture-related traits in order to increase wheat yield potential and end-use quality.