Genetic diversity and linkage disequilibrium in Chinese bread wheat (Triticum aestivum L.) revealed by SSR markers

DIW1 encoding a clade I PP2C phosphatase negatively regulates drought tolerance by de-phosphorylating TaSnRK1.1 in wheat

Drought seriously impacts wheat production (Triticum aestivum L.), while the exploitation and utilization of genes for drought tolerance are insufficient. Leaf wilting is a direct reflection of drought tolerance in plants. Clade A PP2Cs are abscisic acid (ABA) co-receptors playing vital roles in the ABA signaling pathway, regulating drought response. However, the roles of other clade PP2Cs in drought tolerance, especially in wheat, remain largely unknown. Here, we identified a gain-of-function drought-induced wilting 1 (DIW1) gene from the wheat Aikang 58 mutant library by map-based cloning, which encodes a clade I protein phosphatase 2C (TaPP2C158) with enhanced protein phosphatase activity. Phenotypic analysis of overexpression and CRISPR/Cas9 mutant lines demonstrated that DIW1/TaPP2C158 is a negative regulator responsible for drought resistance. We found that TaPP2C158 directly interacts with TaSnRK1.1 and de-phosphorylates it, thus inactivating the TaSnRK1.1-TaAREB3 pathway. TaPP2C158 protein phosphatase activity is negatively correlated with ABA signaling. Association analysis suggested that C-terminal variation of TaPP2C158 changing protein phosphatase activity is highly correlated with the canopy temperature, and seedling survival rate under drought stress. Our data suggest that the favorable allele with lower phosphatase activity of TaPP2C158 has been positively selected in Chinese breeding history. This work benefits us in understanding the molecular mechanism of wheat drought tolerance, and provides elite genetic resources and molecular markers for improving wheat drought tolerance.

Genetic Diversity and Population Structure in Türkiye Bread Wheat Genotypes Revealed by Simple Sequence Repeats (SSR) Markers

Wheat genotypes should be improved through available germplasm genetic diversity to ensure food security. This study investigated the molecular diversity and population structure of a set of Türkiye bread wheat genotypes using 120 microsatellite markers. Based on the results, 651 polymorphic alleles were evaluated to determine genetic diversity and population structure. The number of alleles ranged from 2 to 19, with an average of 5.44 alleles per locus. Polymorphic information content (PIC) ranged from 0.031 to 0.915 with a mean of 0.43. In addition, the gene diversity index ranged from 0.03 to 0.92 with an average of 0.46. The expected heterozygosity ranged from 0.00 to 0.359 with a mean of 0.124. The unbiased expected heterozygosity ranged from 0.00 to 0.319 with an average of 0.112. The mean values of the number of effective alleles (Ne), genetic diversity of Nei (H) and Shannon's information index (I) were estimated at 1.190, 1.049 and 0.168, respectively. The highest genetic diversity (GD) was estimated between genotypes G1 and G27. In the UPGMA dendrogram, the 63 genotypes were grouped into three clusters. The three main coordinates were able to explain 12.64, 6.38 and 4.90% of genetic diversity, respectively. AMOVA revealed diversity within populations at 78% and between populations at 22%. The current populations were found to be highly structured. Model-based cluster analyses classified the 63 genotypes studied into three subpopulations. The values of F-statistic (Fst) for the identified subpopulations were 0.253, 0.330 and 0.244, respectively. In addition, the expected values of heterozygosity (He) for these sub-populations were recorded as 0.45, 0.46 and 0.44, respectively. Therefore, SSR markers can be useful not only in genetic diversity and association analysis of wheat but also in its germplasm for various agronomic traits or mechanisms of tolerance to environmental stresses.

Wheat NAC-A18 regulates grain starch and storage proteins synthesis and affects grain weight

KEY MESSAGE

Wheat NAC-A18 regulates both starch and storage protein synthesis in the grain, and a haplotype with positive effects on grain weight showed increased frequency during wheat breeding in China. Starch and seed storage protein (SSP) directly affect the processing quality of wheat grain. The synthesis of starch and SSP are also regulated at the transcriptional level. However, only a few starch and SSP regulators have been identified in wheat. In this study, we discovered a NAC transcription factor, designated as NAC-A18, which acts as a regulator of both starch and SSP synthesis. NAC-A18, is predominately expressed in wheat developing grains, encodes a transcription factor localized in the nucleus, with both activation and repression domains. Ectopic expression of wheat NAC-A18 in rice significantly decreased starch accumulation and increased SSP accumulation and grain size and weight. Dual-luciferase reporter assays indicated that NAC-A18 could reduce the expression of TaGBSSI-A1 and TaGBSSI-A2, and enhance the expression of TaLMW-D6 and TaLMW-D1. A yeast one hybrid assay demonstrated that NAC-A18 bound directly to the cis-element "ACGCAA" in the promoters of TaLMW-D6 and TaLMW-D1. Further analysis indicated that two haplotypes were formed at NAC-A18, and that NAC-A18_h1 was a favorable haplotype correlated with higher thousand grain weight. Based on limited population data, NAC-A18_h1 underwent positive selection during Chinese wheat breeding. Our study demonstrates that wheat NAC-A18 regulates starch and SSP accumulation and grain size. A molecular marker was developed for the favorable allele for breeding applications.

Tiller Number1 encodes an ankyrin repeat protein that controls tillering in bread wheat

Wheat (Triticum aestivum L.) is a major staple food for more than one-third of the world's population. Tiller number is an important agronomic trait in wheat, but only few related genes have been cloned. Here, we isolate a wheat mutant, tiller number1 (tn1), with much fewer tillers. We clone the TN1 gene via map-based cloning: TN1 encodes an ankyrin repeat protein with a transmembrane domain (ANK-TM). We show that a single amino acid substitution in the third conserved ankyrin repeat domain causes the decreased tiller number of tn1 mutant plants. Resequencing and haplotype analysis indicate that TN1 is conserved in wheat landraces and modern cultivars. Further, we reveal that the expression level of the abscisic acid (ABA) biosynthetic gene TaNCED3 and ABA content are significantly increased in the shoot base and tiller bud of the tn1 mutants; TN1 but not tn1 could inhibit the binding of TaPYL to TaPP2C via direct interaction with TaPYL. Taken together, we clone a key wheat tiller number regulatory gene TN1, which promotes tiller bud outgrowth probably through inhibiting ABA biosynthesis and signaling.

Genetic Diversity and Population Structure Analysis of Castanopsis hystrix and Construction of a Core Collection Using Phenotypic Traits and Molecular Markers

Castanopsis hystrix is a valuable native, broad-leaved, and fast-growing tree in South China. In this study, 15 phenotypic traits and 32 simple sequence repeat (SSR) markers were used to assess the genetic diversity and population structure of a natural population of C. hystrix and to construct a core germplasm collection by a set of 232 accessions. The results showed that the original population of C. hystrix had relatively high genetic diversity, with the number of alleles (Na), effective number of alleles (Ne), observed heterozygosity (Ho), expected heterozygosity (He), Shannon's information index (I), and polymorphism information content (PIC) averaging at 26.188, 11.565, 0.863, 0.897, 2.660, and 0.889, respectively. Three sub-populations were identified based on a STRUCTURE analysis, indicating a strong genetic structure. The results from the phylogenetic and population structures showed a high level of agreement, with 232 germplasms being classified into three main groups. The analysis of molecular variance (AMOVA) test indicated that 96% of the total variance was derived from within populations, which revealed a low differentiation among populations. A core collection composed of 157 germplasms was firstly constructed thereafter, of which the diversity parameters non-significantly differed from the original population. These results revealed the genetic diversity and population structure of C. hystrix germplasms, which have implications for germplasm management and genome-wide association studies on C. hystrix, as well as for core collection establishment applications in other wood-producing hardwood species.

TaGSNE, a WRKY transcription factor, overcomes the trade-off between grain size and grain number in common wheat and is associated with root development

Wheat is one of the world's major staple food crops, and breeding for improvement of grain yield is a priority under the scenarios of climate change and population growth. WRKY transcription factors are multifaceted regulators in plant growth, development, and responses to environmental stimuli. In this study, we identify the WRKY gene TaGSNE (Grain Size and Number Enhancer) in common wheat, and find that it has relatively high expression in leaves and roots, and is induced by multiple abiotic stresses. Eleven single-nucleotide polymorphisms were identified in TaGSNE, forming two haplotypes in multiple germplasm collections, named as TaGSNE-Hap-1 and TaGSNE-Hap-2. In a range of different environments, TaGSNE-Hap-2 was significantly associated with increases in thousand-grain weight (TGW; 3.0%) and spikelet number per spike (4.1%), as well as with deeper roots (10.1%) and increased root dry weight (8.3%) at the mid-grain-filling stage, and these were confirmed in backcross introgression populations. Furthermore, transgenic rice lines overexpressing TaGSNE had larger panicles, more grains, increased grain size, and increased grain yield relative to the wild-type control. Analysis of geographic and temporal distributions revealed that TaGSNE-Hap-2 is positively selected in China and Pakistan, and TaGSNE-Hap-1 in Europe. Our findings demonstrate that TaGSNE overcomes the trade-off between TGW/grain size and grain number, leading us to conclude that these elite haplotypes and their functional markers could be utilized in marker-assisted selection for breeding high-yielding varieties.

Folate (vitamin B9) content analysis in bread wheat (Triticum aestivum L.)

Vitamin B9, particularly folic acid, is an essential molecule for human health. Wheat flour is one of the major sources of calorie intake by humans. The selection of folate-rich genotypes in wheat breeding can enhance the natural folate value in the daily diet. This study used a precise, high-performance liquid chromatography (HPLC) assay to analyze folate content in a 262-accession Chinese wheat mini-core collection (MCC) grown under three environments. Four folate derivatives in grains including tetrahydrofolate (THF), 5-methyltetrahydrofolate (5-CH₃-THF), 5-formyltetrahydrofolate (5-CHO-THF), and 5,10-methenyltetrahydrofolate (5,10-CH⁺THF) were considered. An association analysis of water regimes, accession types, released years, geographical origin, and agronomic traits with folate content was conducted for the first time. There was a large amount of variation in folate content in the analyzed accessions, with genotype identified as the main influencing factor. Total folate content was significantly correlated with the content of the four MCC derivatives under the three environments. 5-CH₃-THF and 5-CHO-THF were the most abundant among the four folate derivatives and were positively correlated with high folate content. The 12 accessions with the highest folate content showed an average of more than 80 μg/100 g. The analysis demonstrated that this Chinese wheat had not undergone extensive selection for folate content during breeding, which is unrelated to the geographical origin, accession types, winter/spring types, and grain colors of wheat. The content of THF, 5-CH₃-THF, and 5,10-CH⁺THF was significantly negatively correlated with grain width, grain thickness, and thousand kernel weight. A relatively weak negative relationship manifested between folate contents and flowering date, whereas no significant correlation with tiller number, grain number per spike, maturity date, height, and spike length was detected. The investigation benefits wheat breeders for folate enhancement.

Building a Calibration Set for Genomic Prediction, Characteristics to Be Considered, and Optimization Approaches

The efficiency of genomic selection strongly depends on the prediction accuracy of the genetic merit of candidates. Numerous papers have shown that the composition of the calibration set is a key contributor to prediction accuracy. A poorly defined calibration set can result in low accuracies, whereas an optimized one can considerably increase accuracy compared to random sampling, for a same size. Alternatively, optimizing the calibration set can be a way of decreasing the costs of phenotyping by enabling similar levels of accuracy compared to random sampling but with fewer phenotypic units. We present here the different factors that have to be considered when designing a calibration set, and review the different criteria proposed in the literature. We classified these criteria into two groups: model-free criteria based on relatedness, and criteria derived from the linear mixed model. We introduce criteria targeting specific prediction objectives including the prediction of highly diverse panels, biparental families, or hybrids. We also review different ways of updating the calibration set, and different procedures for optimizing phenotyping experimental designs.

Genome-Wide Association Study of Grain Number in Common Wheat From Shanxi Under Different Water Regimes

Water availability is a crucial environmental factor on grain number in wheat, which is one of the important yield-related traits. In this study, a diverse panel of 282 wheat accessions were phenotyped for grain number per spike (GNS), spikelet number (SN), basal sterile spikelet number (BSSN), and apical sterile spikelet number (ASSN) under different water regimes across two growing seasons. Correlation analysis showed that GNS is significantly correlated with both SN and BSSN under two water regimes. A total of 9,793 single nucleotide polymorphism (SNP) markers from the 15 K wheat array were employed for genome-wide association study (GWAS). A total of 77 significant marker-trait associations (MTAs) for investigated traits as well as 8 MTAs for drought tolerance coefficient (DTC) were identified using the mixed linear model. Favored alleles for breeding were inferred according to their estimated effects on GNS, based on the mean difference of varieties. Frequency changes in favored alleles associated with GNS in modern varieties indicate there is still considerable genetic potential for their use as markers for genome selection of GNS in wheat breeding.

Association analysis for agronomic traits in wheat under terminal heat stress

•

Terminal heat stress leads to irreversible damage in wheat.

•

Marker assisted selection and gene pyramiding for portrayal of heat tolerance.

•

Allelic frequency and polymorphic information showed significant variability.

•

Markers xcfa2147 and xwmc671 could be potentail for heat stress tolerance.

Terminal heat stress causes irreversible damage to wheat crop productivity. It reduces the vegetative growth and flowering period that consequently declines the efficiency to capture available stem reserves (carbohydrates) in grains. Markers associated with thermotolerant traits ease in marker assisted selection (MAS) for crop improvement. It identifies the genomic regions associated with thermotolerant traits in wheat, but the scarcity of markers is the major hindrance in crop improvement. Therefore, 158 wheat genotypes were subjected to genotyping with 165 simple sequence repeat markers dispersed on three genomes (A, B and D). Allelic frequency and polymorphic information content values were highest on genome A (5.34 (14% greater than the lowest value at genome D) and 0.715 (3% greater than the lowest value at genome D)), chromosome 4 (5.40 (16% greater than the lowest value at chromosome 2) and 0.725 (5% greater than the lowest value at chromosome 6)) and marker xgwm44 (13.0 (84% greater than the lowest value at marker xbarc148) and 0.916 (46% greater than the lowest value at marker xbarc148)). Bayesian based population structure discriminated the wheat genotypes into seven groups based on genetic similarity indicating their ancestral origin and geographical ecotype. Linkage disequilibrium pattern had highest significant (P < 0.001) linked loci pairs 732 on genome A at r² > 0.1 whereas, 58 on genome B at r² > 0.5. Linkage disequilibrium decay (P < 0.01 and r² > 0.1) had larger LD block (5–10 cM) on genome A. Highly significant MTAs (P < 0.000061) under heat stress conditions were identified for flag leaf area (xwmc336), spikelet per spike (xwmc553), grains per spike (cxfa2147, xwmc418 and xwmc121), biomass (xbarc7) and grain yield (xcfa2147 and xwmc671). The identified markers in this study could facilitate in MAS and gene pyramiding against heat stress in wheat.

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.

Association Analysis Revealed That TaPYL4 Genes Are Linked to Plant Growth Related Traits in Multiple Environment

Abscisic acid (ABA), one of phytohormones, plays an important regulatory role in plant growth and development. ABA receptor PYL4 (pyrabactin resistance 1-like 4) was previously detected to be involved in plant response to a variety of stresses. TaPYL4 overexpression could enhance wheat (Triticum aestivum) drought resistance. In order to further investigate TaPYL4's role in regulating development of other major agronomic traits in wheat, genes of TaPYL4-2A, TaPYL4-2B, and TaPYL4-2D were cloned from wheat, respectively. Polymorphism analysis on TaPYL4 sequences revealed that encoding regions of the three genes were highly conserved, without any SNP (single nucleotide polymorphism) presence. However, nine SNPs and four SNPs were identified in the promoter regions of TaPYL4-2A and TaPYL4-2B, respectively. Functional molecular markers were developed based on these polymorphisms, which were then used to scan a natural population of 323 common wheat accessions for correlation analysis between genotype and the target phenotypic traits. Both TaPYL4-2A and TaPYL4-2B markers were significantly correlated with plant growth-related traits under multiple environments (well-watered, drought and heat stress treatments). The additive effects of TaPYL4-2A and TaPYL4-2B were verified by the combinational haplotype (Hap-AB1∼Hap-AB4) effects determined from field data. Cis-acting elements were analyzed in the promoters of TaPYL4-2A and TaPYL4-2B, showing that a TGA-element bound by ARFs (auxin response factors) existed only in Hap-2A-1 of TaPYL4-2A. Gene expression assays indicated that TaPYL4-2A was constitutively expressed in various tissues, with higher expression in Hap-2A-1 genotypes than in Hap-2A-2 materials. Notably, TaARF4 could act as TaPYL4-2A transcription activator in Hap-2A-1 materials, but not in Hap-2A-2 genotypes. Analysis of geographic distribution and temporal frequency of haplotypes indicated that Hap-AB1 was positively selected in wheat breeding in China. Therefore, TaPYL4-2A and TaPYL4-2B could be a valuable target gene in wheat genetic improvement to develop the ideal plant architecture.

Association of TaD14-4D, a Gene Involved in Strigolactone Signaling, with Yield Contributing Traits in Wheat

Tillering is a crucial agronomic trait of wheat; it determines yield and plant architecture. Strigolactones (SLs) have been reported to inhibit plant branching. D14, a receptor of SLs, has been described to affect tillering in rice, yet it has seldomly been studied in wheat. In this study, three TaD14 homoeologous genes, TaD14-4A, TaD14-4B, and TaD14-4D, were identified. TaD14-4A, TaD14-4B, and TaD14-4D were constitutively expressed, and TaD14-4D had a higher expression level in most tissues. TaD14 proteins were localized in both cytoplasm and nucleus. An SNP and a 22 bp insertion/deletion (Indel) at the exon regions of TaD14-4D were detected, forming three haplotypes, namely 4D-HapI, 4D-HapII, and 4D-HapIII. Due to the frameshift mutation in the coding region of 4D-HapII, the interaction of 4D-HapII with TaMAX2 and TaD53 was blocked, which led to the blocking of SL signal transduction. Based on the two variation sites, two molecular markers, namely dCAPS-250 and Indel-747, were developed. Association analysis suggested that haplotypes of TaD14-4D were associated with effective tillering number (ETN) and thousand kernel weight (TKW) simultaneously in four environments. The favorable haplotype 4D-HapIII underwent positive selection in global wheat breeding. This study provides insights into understanding the function of natural variations of TaD14-4D and develops two useful molecular markers for wheat breeding.

Variety Discrimination Power: An Appraisal Index for Loci Combination Screening Applied to Plant Variety Discrimination

Molecular marker technology is used widely in plant variety discrimination, molecular breeding, and other fields. To lower the cost of testing and improve the efficiency of data analysis, molecular marker screening is very important. Screening usually involves two phases: the first to control loci quality and the second to reduce loci quantity. To reduce loci quantity, an appraisal index that is very sensitive to a specific scenario is necessary to select loci combinations. In this study, we focused on loci combination screening for plant variety discrimination. A loci combination appraisal index, variety discrimination power (VDP), is proposed, and three statistical methods, probability-based VDP (P-VDP), comparison-based VDP (C-VDP), and ratio-based VDP (R-VDP), are described and compared. The results using the simulated data showed that VDP was sensitive to statistical populations with convergence toward the same variety, and the total probability of discrimination power (TDP) method was effective only for partial populations. R-VDP was more sensitive to statistical populations with convergence toward various varieties than P-VDP and C-VDP, which both had the same sensitivity; TDP was not sensitive at all. With the real data, R-VDP values for sorghum, wheat, maize and rice data begin to show downward tendency when the number of loci is 20, 7, 100, 100 respectively, while in the case of P-VDP and C-VDP (which have the same results), the number is 6, 4, 9, 19 respectively and in the case of TDP, the number is 6, 4, 4, 11 respectively. For the variety threshold setting, R-VDP values of loci combinations with different numbers of loci responded evenly to different thresholds. C-VDP values responded unevenly to different thresholds, and the extent of the response increased as the number of loci decreased. All the methods gave underestimations when data were missing, with systematic errors for TDP, C-VDP, and R-VDP going from smallest to biggest. We concluded that VDP was a better loci combination appraisal index than TDP for plant variety discrimination and the three VDP methods have different applications. We developed the software called VDPtools, which can calculate the values of TDP, P-VDP, C-VDP, and R-VDP. VDPtools is publicly available at https://github.com/caurwx1/VDPtools.git.

A thorough screening based on QTLs controlling zinc and copper accumulation in the grain of different wheat genotypes

Excess trace metals may cause damage to human health due to the consumption of food grain grown in contaminated soils. This study was designed to understand the genetic mechanisms of copper (Cu) and zinc (Zn) accumulation in wheat grain under stressed environments. The differences of Cu/Zn contents in the grain among 246 wheat varieties were analyzed, and the wheat varieties with low or high accumulation of Cu and Zn in the safe range were also screened out. The accumulation of Cu and Zn in grains of "Chushanbao" was lowest, which could be used as a novel germplasm for wheat breeding under heavy metal stress. We found that Cu contents of wheat grain were significantly and positively correlated with Zn. The quantitative trait loci (QTLs) for grain Cu content (GCuC) and grain Zn content (GZnC) were detected by genome-wide association study (GWAS). Twenty-three loci affecting GCuC were identified on chromosomes 1A, 1D, 2A, 2B, 2D, 3A, 3B, 3D, 4A, 4B 4D, 5A, 6D, 7A, and 7B, explaining 2.6-5.8% of the phenotypic variation. Sixteen loci associated with the GZnC on 11 different chromosomes 1B, 2B, 2D, 3A, 3D, 4A, 4B, 5A, 5D, 6B, and 7D were detected, which could explain 2.7~6.6% of phenotypic variance. We also determined five associated loci on chromosomes 2B, 2D, 3A, 4B, and 5A were in pleiotropic regions affecting both GCuC and GZnC. This study would help in better understanding the molecular basis of Cu/Zn accumulation in wheat grain, and the associated markers may be useful for marker-assisted selection (MAS) breeding program.

A large-scale genomic association analysis identifies the candidate causal genes conferring stripe rust resistance under multiple field environments

The incorporation of resistance genes into wheat commercial varieties is the ideal strategy to combat stripe or yellow rust (YR). In a search for novel resistance genes, we performed a large-scale genomic association analysis with high-density 660K single nucleotide polymorphism (SNP) arrays to determine the genetic components of YR resistance in 411 spring wheat lines. Following quality control, 371 972 SNPs were screened, covering over 50% of the high-confidence annotated gene space. Nineteen stable genomic regions harbouring 292 significant SNPs were associated with adult-plant YR resistance across nine environments. Of these, 14 SNPs were localized in the proximity of known loci widely used in breeding. Obvious candidate SNP variants were identified in certain confidence intervals, such as the cloned gene Yr18 and the major locus on chromosome 2BL, despite a large extent of linkage disequilibrium. The number of causal SNP variants was refined using an independent validation panel and consideration of the estimated functional importance of each nucleotide polymorphism. Interestingly, four natural polymorphisms causing amino acid changes in the gene TraesCS2B01G513100 that encodes a serine/threonine protein kinase (STPK) were significantly involved in YR responses. Gene expression and mutation analysis confirmed that STPK played an important role in YR resistance. PCR markers were developed to identify the favourable TraesCS2B01G513100 haplotype for marker-assisted breeding. These results demonstrate that high-resolution SNP-based GWAS enables the rapid identification of putative resistance genes and can be used to improve the efficiency of marker-assisted selection in wheat disease resistance breeding.

Transcription Factor TaWRKY51 Is a Positive Regulator in Root Architecture and Grain Yield Contributing Traits

Wheat is one of the staple food crops. The utilization of elite genetic resources to develop resource-efficient wheat varieties is an effective approach to deal with the challenges of climate change and population growth. WRKY transcription factors (TFs) are multifaceted regulators of plant growth and development and response to environmental stress. The previous studies have shown that TaWRKY51 positively regulates the development of lateral roots, while its roles in agronomic trait development are not clear, and there is no functional marker for molecular breeding. To bridge the gap, we cloned the three members of TaWRKY51 and found they were highly expressed in the roots and flag leaves at the flowering stage and were induced by the multiple abiotic stresses and phytohormones. The highest expression level was observed in TaWRKY51-2D, followed by TaWRKY51-2A and -2B. The two haplotypes/alleles for each member were identified in the natural populations, and functional markers were developed accordingly. The association assays revealed that Hap-2A-I was an elite haplotype for the large spike, Hap-2B-II and allele-G were favorable haplotypes/alleles for long root. However, only Hap-2A-I was selected for wheat breeding in China. The results of transgenic experiments showed that the rice lines overexpressing TaWRKY51 had large panicle, high thousand-grain-weight, and more crown and lateral roots, which further confirmed the results of association analysis. In short, TaWRKY51 is a positive regulator of the root architecture and grain yield (GY) contributing traits. The elite gene resources and functional markers may be utilized in the marker-assisted selection for high-yield breeding in wheat.

Population structure of Nepali spring wheat (Triticum aestivum L.) germplasm

BACKGROUND

Appropriate information about genetic diversity and population structure of germplasm improves the efficiency of plant breeding. The low productivity of Nepali bread wheat (Triticum aestivum L.) is a major concern particularly since Nepal is ranked the 4th most vulnerable nation globally to climate change. The genetic diversity and population structure of Nepali spring wheat have not been reported. This study aims to improve the exploitation of more diverse and under-utilized genetic resources to contribute to current and future breeding efforts for global food security.

RESULTS

We used genotyping-by-sequencing (GBS) to characterize a panel of 318 spring wheat accessions from Nepal including 166 landraces, 115 CIMMYT advanced lines, and 34 Nepali released varieties. We identified 95 K high-quality SNPs. The greatest genetic diversity was observed among the landraces, followed by CIMMYT lines, and released varieties. Though we expected only 3 groupings corresponding to these 3 seed origins, the population structure revealed two large, distinct subpopulations along with two smaller and scattered subpopulations in between, with significant admixture. This result was confirmed by principal component analysis (PCA) and UPGMA distance-based clustering. The pattern of LD decay differed between subpopulations, ranging from 60 to 150 Kb. We discuss the possibility that germplasm explorations during the 1970s-1990s may have mistakenly collected exotic germplasm instead of local landraces and/or collected materials that had already cross-hybridized since exotic germplasm was introduced starting in the 1950s.

CONCLUSION

We suggest that only a subset of wheat "landraces" in Nepal are authentic which this study has identified. Targeting these authentic landraces may accelerate local breeding programs to improve the food security of this climate-vulnerable nation. Overall, this study provides a novel understanding of the genetic diversity of wheat in Nepal and this may contribute to global wheat breeding initiatives.

The MYB transcription factor TaPHR3-A1 is involved in phosphate signaling and governs yield-related traits in bread wheat

Improved inorganic phosphate (Pi) use efficiency in crops will be important for sustainable agriculture. Exploring molecular mechanisms that regulate Pi uptake could provide useful information for breeding wheat with improved Pi use efficiency. Here, a TaPHR3-A1 (Gene ID: TraesCS7A02G415800) ortholog of rice OsPHR3 that functions in transcriptional regulation of Pi signaling was cloned from wheat chromosome 7A. Ectopic expression of TaPHR3-A1 in Arabidopsis and rice produced enhanced vegetative growth and more seeds. Overexpression in transgenic rice led to increased biomass, grain number, and primary panicle branching by 61.23, 42.12, and 36.34% compared with the wild type. Transgenic wheat lines with down-regulation of TaPHR3-A1 exhibited retarded growth and root hair development at the seedling stage, and showed yield-related effects at the adult stage when grown in both low- and sufficient Pi conditions, indicating that TaPHR3-A1 positively regulated tolerance to low Pi. Introgression lines further confirmed the effect of TaPHR3-A1 in improving grain number. The Chinese wheat mini core collection and a recombinant inbred line analysis demonstrated that the favorable allele TaPHR3-A1-A associated with higher grain number was positively selected in breeding. A TaPHR3-A1-derived cleaved amplified polymorphic sequence marker effectively identified haplotype TaPHR3-A1-A. Our results suggested that TaPHR3-A1 was a functional regulatory factor for Pi uptake and provided useful information for marker-assisted selection for high yield in wheat.

Starch Metabolism in Wheat: Gene Variation and Association Analysis Reveal Additive Effects on Kernel Weight

Kernel weight is a key determinant of yield in wheat (Triticum aestivum L.). Starch consists of amylose and amylopectin and is the major constituent of mature grain. Therefore, starch metabolism in the endosperm during grain filling can influence kernel weight. In this study, we sequenced 87 genes involved in starch metabolism from 300 wheat accessions and detected 8,141 polymorphic sites. We also characterized yield-related traits across different years in these accessions. Although the starch contents fluctuated, thousand kernel weight (TKW) showed little variation. Polymorphisms in six genes were significantly associated with TKW. These genes were located on chromosomes 2A, 2B, 4A, and 7A; none were associated with starch content or amylose content. Variations of 15 genes on chromosomes 1A and 7A formed haplotype blocks in 26 accessions. Notably, accessions with higher TKWs had more of the favorable haplotypes. We thus conclude that these haplotypes contribute additive effects to TKW.

Genome-wide association analysis of stem water-soluble carbohydrate content in bread wheat

GWAS identified 36 potentially new loci for wheat stem water-soluble carbohydrate (WSC) contents and 13 pleiotropic loci affecting WSC and thousand-kernel weight. Five KASP markers were developed and validated. Water-soluble carbohydrates (WSC) reserved in stems contribute significantly to grain yield (GY) in wheat. However, knowledge of the genetic architecture underlying stem WSC content (SWSCC) is limited. In the present study, 166 diverse wheat accessions from the Yellow and Huai Valleys Winter Wheat Zone of China and five other countries were grown in four well-watered environments. SWSCC at 10 days post-anthesis (10DPA), 20DPA and 30DPA, referred as WSC10, WSC20 and WSC30, respectively, and thousand-kernel weight (TKW) were assessed. Correlation analysis showed that TKW was significantly and positively correlated with WSC10 and WSC20. Genome-wide association study was performed on SWSCC and TKW with 373,106 markers from the wheat 660 K and 90 K SNP arrays. Totally, 62 stable loci were detected for SWSCC, with 36, 24 and 19 loci for WSC10, WSC20 and WSC30, respectively; among these, 36 are potentially new, 16 affected SWSCC at two or three time-points, and 13 showed pleiotropic effects on both SWSCC and TKW. Linear regression showed clear cumulative effects of favorable alleles for increasing SWSCC and TKW. Genetic gain analyses indicated that pyramiding favorable alleles of SWSCC had simultaneously improved TKW. Kompetitive allele-specific PCR markers for five pleiotropic loci associated with both SWSCC and TKW were developed and validated. This study provided a genome-wide landscape of the genetic architecture of SWSCC, gave a perspective for understanding the relationship between WSC and GY and explored the theoretical basis for co-improvement of WSC and GY. It also provided valuable loci and markers for future breeding.

Genetic Diversity and Population Structure of Asian and European Common Wheat Accessions Based on Genotyping-By-Sequencing

Obtaining information on the genetic diversity and population structure of germplasm facilitates its use in wheat breeding programs. Recently, with the development of next-generation sequencing technology, genotyping-by-sequencing (GBS) has been used as a high-throughput and cost-effective molecular tool for examination of the genetic diversity of wheat breeding lines. In this study, GBS was used to characterize a population of 180 accessions of common wheat originating from Asia and Europe between the latitudes 30° and 45°N. In total, 24,767 high-quality single-nucleotide polymorphism (SNP) markers were used for analysis of genetic diversity and population structure. The B genome contained the highest number of SNPs, followed by the A and D genomes. The polymorphism information content was in the range of 0.1 to 0.4, with a mean of 0.26. The distribution of SNPs markers on the 21 chromosomes ranged from 243 on chromosome 4D to 2,337 on chromosome 3B. Structure and cluster analyses divided the panel of accessions into two subgroups (G1 and G2). G1 principally consisted of European and partial Asian accessions, and G2 comprised mainly accessions from the Middle East and partial Asia. Molecular analysis of variance showed that the genetic variation was greater within groups (99%) than between groups (1%). Comparison of the two subgroups indicated that G1 and G2 contained a high level of genetic diversity. The genetic diversity of G2 was slightly higher as indicated by the observed heterozygosity (H o) = 0.23, and unbiased diversity index (uh) = 0.34. The present results will not only help breeders to understand the genetic diversity of wheat germplasm on the Eurasian continent between the latitudes of 30° and 45°N, but also provide valuable information for wheat genetic improvement through introgression of novel genetic variation in this region.

RING finger ubiquitin E3 ligase gene TaSDIR1-4A contributes to determination of grain size in common wheat

Salt and drought-induced RING finger1 (SDIR1) is a RING-type E3 ubiquitin ligase that plays a key role in ABA-mediated responses to salinity and drought stress via the ubiquitination pathway in some plant species. However, its function in wheat (Triticum aestivum) is unknown. Here, we isolated a SDIR1 member in wheat, TaSDIR1-4A, and characterized its E3 ubiquitin ligase activity. DNA polymorphism assays showed the presence of two nucleotide variation sites in the promoter region of TaSDIR1-4A, leading to the detection of the haplotypes Hap-4A-1 and Hap-4A-2 in wheat populations. Association analysis showed that TaSDIR1-4A haplotypes were associated with 1000-grain weight (TGW) across a variety of different environments, including well-watered and heat-stress conditions. Genotypes with Hap-4A-2 had higher TGW than those with Hap-4A-1. Phenotypes in both gene-silenced wheat and transgenic Arabidopsis showed that TaSDIR1-4A was a negative regulator of grain size. Gene expression assays indicated that TaSDIR1-4A was most highly expressed in flag leaves, and expression was higher in Hap-4A-1 accessions than in Hap-4A-2 accessions. The difference might be attributable to the fact that TaERF3 (ethylene response factor) can act as a transcriptional repressor of TaSDIR1-4A in Hap-4A-2 but not in Hap-4A-1. Examination of modern wheat varieties shows that the favorable haplotype has been positively selected in breeding programs in China. The functional marker for TaSDIR1-4A developed in this study should be helpful for future wheat breeding.

Selection criteria for high-yielding and early-flowering bread wheat hybrids under heat stress

Hybrid performance during wheat breeding can be improved by analyzing genetic distance (GD) among wheat genotypes and determining its correlation with heterosis. This study evaluated the GD between 16 wheat genotypes by using 60 simple sequence repeat (SSR) markers to classify them according to their relationships and select those with greater genetic diversity, evaluate the correlation of the SSR marker distance with heterotic performance and specific combining ability (SCA) for heat stress tolerance, and identify traits that most influence grain yield (GY). Eight parental genotypes with greater genetic diversity and their 28 F1 hybrids generated using diallel crossing were evaluated for 12 measured traits in two seasons. The GD varied from 0.235 to 0.911 across the 16 genotypes. Cluster analysis based on the GD estimated using SSRs classified the genotypes into three major groups and six sub-groups, almost consistent with the results of principal coordinate analysis. The combined data indicated that five hybrids showed 20% greater yield than mid-parent or better-parent. Two hybrids (P2 × P4) and (P2 × P5), which showed the highest performance of days to heading (DH), grain filling duration (GFD), and GY, and had large genetic diversity among themselves (0.883 and 0.911, respectively), were deemed as promising heat-tolerant hybrids. They showed the best mid-parent heterosis and better-parent heterosis (BPH) for DH (-11.57 and -7.65%; -13.39 and -8.36%, respectively), GFD (12.74 and 12.17%; 12.09 and 10.59%, respectively), and GY (36.04 and 20.04%; 44.06 and 37.73%, respectively). Correlation between GD and each of BPH and SCA effects based on SSR markers was significantly positive for GFD, hundred kernel weight, number of kernels per spike, harvest index, GY, and grain filling rate and was significantly negative for DH. These correlations indicate that the performance of wheat hybrids with high GY and earliness could be predicted by determining the GD of the parents by using SSR markers. Multivariate analysis (stepwise regression and path coefficient) suggested that GFD, hundred kernel weight, days to maturity, and number of kernels per spike had the highest influence on GY.

Genome-Wide Association Mapping of Adult-Plant Resistance to Stripe Rust in Common Wheat (Triticum aestivum)

Stripe rust, caused by Puccinia striiformis f. sp. tritici, is a globally devastating disease of common wheat (Triticum aestivum L.), resulting in substantial economic losses. To identify effective resistance genes, a genome-wide association study was conducted on 120 common wheat lines from different wheat-growing regions of China using the wheat 90K iSelect SNP array. Seventeen loci were identified, explaining 9.5 to 21.8% of the phenotypic variation. Most of these genes were detected in the A (seven) and B (seven) genomes, with only three in the D genome. Among them, 11 loci were colocated with known resistance genes or quantitative trait loci reported previously, whereas the other six are likely new resistance loci. Annotation of flanking sequences of significantly associated SNPs indicated the presence of three important candidate genes, including E3 ubiquitin-protein ligase, F-box repeat protein, and disease resistance RPP13-like protein. This study increased our knowledge in understanding the genetic architecture for stripe rust resistance and identified wheat varieties with multiple resistance alleles, which are useful for improvement of stripe rust resistance in breeding.

A Genome-Wide Association Study of Coleoptile Length in Different Chinese Wheat Landraces

From the perspective of wheat yield improvement, the coleoptile is vital for successful crop establishment, and long coleoptile lengths (CLs) are preferred in wheat-growing regions where deep planting is practiced. To determine the genetic basis of CL, we performed a genome-wide association study on a set of 707 Chinese wheat landraces using 18,594 single-nucleotide polymorphisms and 38,678 diversity array technology sequencing markers. We accordingly detected a total of 29 significant markers [-log₁₀ ^(P) > 4.76] distributed on chromosomes 2B, 2D, 3A, 4A, 5A, 6A, 6B, 6D, and 7B. Based on linkage disequilibrium decay distance, we identified a total of 17 quantitative trait loci associated with CL, among which QCl.sicau-6B.2, located at 508.17-509.26 Mb on chromosome 6B, was recognized as a novel major locus. We subsequently developed a high-resolution melt marker for QCl.sicau-6B.2, which was validated in an F _{2 : 3} population. Our findings provide important insights into the genetic mechanisms underlying coleoptile growth and could be applied to marker-assisted wheat selection.

Association Analysis Identifies New Loci for Resistance to Chinese Yr26-Virulent Races of the Stripe Rust Pathogen in a Diverse Panel of Wheat Germplasm

Stripe rust caused by Puccinia striiformis f. sp. tritici (Pst) is one of the most destructive fungal diseases of wheat worldwide. The expanding Yr26-virulent Pst race (V26) group overcomes almost all currently deployed resistance genes in China and has continued to accumulate new virulence. Investigating the genetic architecture of stripe rust resistance in common wheat is an important basis for a successful utilization of resistance in breeding programs. A panel of 410 exotic wheat germplasms was used for characterizing new stripe rust resistance loci. This panel was genotyped using high-density wheat 660K single-nucleotide polymorphism (SNP) array, and phenotypic evaluation of seedlings for stripe rust resistance was performed using multiple Pst races. Thirty-five loci conferring resistance were identified through genome-wide association mapping, and explained phenotypic variances ranged from 53 to 75%. Of these, 14 were colocated in the proximity of the known loci, including cataloged Yr genes Yr9, Yr10, Yr26, Yr33, Yr47, Yr56, Yr57, Yr64, Yr67, Yr72, and Yr81 and three temporarily designated as YrCen, YrNP63, and YrRC detected in our quantitative trait locus (QTL) mapping studies. Seven of them (Yr9, Yr10, Yr24/26, Yr81, YrCEN, YrNP63, and YrRC) were confirmed by molecular detection or genetic analysis. New loci that were identified to be different from reported Yr genes need further confirmation. Nine QTL with significantly large phenotypic effect on resistance to all tested races were considered as major loci for effective resistance. The identified loci enrich our stripe rust resistance gene pool, and the linked SNPs should be useful for marker-assisted selection in breeding programs.

The Circadian Clock Gene, TaPRR1, Is Associated With Yield-Related Traits in Wheat (Triticum aestivum L.)

Timing of flowering is crucial for the transformation from vegetative to reproductive growth in the important food crop, wheat (Triticum aestivum L.). The circadian clock is a central part of photoperiod regulation, with Pseudo-Response Regulators (PRRs) representing key components of circadian networks. However, little is known about the effects of PRR family members on yield-related traits in crop plants. In this study, we identified polymorphisms and haplotypes of TaPRR1, demonstrating that natural variations in TaPRR1 are associated with significant differences in yield-related traits including heading date, plant height and thousand grain weight. TaPRR1-6A-Hapla showed an earlier heading date, advanced by 0.9 to 1.7%. TaPRR1-6B-Hapla and TaPRR1-6D-Hapla displayed reduced plant height and increased thousand grain weight of up to 13.3 to 26.4% and 6.3 to 17.3%, respectively. Subcellular localization and transcriptional activity analysis showed that TaPRR1 is a nuclear localization protein with transcriptional activity controlled by an IR domain. The expression profiles of TaPRR1 genes over a 48-h period were characterized by circadian rhythms, which had two peaks under both short- and long- day conditions. In addition, geographical distribution analysis indicated higher distribution frequencies of TaPRR1-6A-Hapla, TaPRR1-6B-Haplb, and TaPRR1-6D-Haplb in different agro-ecological production zones. Furthermore, analysis of molecular variance of the distribution frequency of TaPRR1 haplotypes suggested significant differences in haplotype distribution frequency between landraces and modern cultivars. Our study provides a basis for in-depth understanding of TaPRR1 function on yield-related traits in wheat, as well as establishing theoretical guidance for wheat molecular marker-assisted breeding.

Disentangling group specific QTL allele effects from genetic background epistasis using admixed individuals in GWAS: An application to maize flowering

When handling a structured population in association mapping, group-specific allele effects may be observed at quantitative trait loci (QTLs) for several reasons: (i) a different linkage disequilibrium (LD) between SNPs and QTLs across groups, (ii) group-specific genetic mutations in QTL regions, and/or (iii) epistatic interactions between QTLs and other loci that have differentiated allele frequencies between groups. We present here a new genome-wide association (GWAS) approach to identify QTLs exhibiting such group-specific allele effects. We developed genetic materials including admixed progeny from different genetic groups with known genome-wide ancestries (local admixture). A dedicated statistical methodology was developed to analyze pure and admixed individuals jointly, allowing one to disentangle the factors causing the heterogeneity of allele effects across groups. This approach was applied to maize by developing an inbred "Flint-Dent" panel including admixed individuals that was evaluated for flowering time. Several associations were detected revealing a wide range of configurations of allele effects, both at known flowering QTLs (Vgt1, Vgt2 and Vgt3) and new loci. We found several QTLs whose effect depended on the group ancestry of alleles while others interacted with the genetic background. Our GWAS approach provides useful information on the stability of QTL effects across genetic groups and can be applied to a wide range of species.

Dissecting the Genetic Complexity of Fusarium Crown Rot Resistance in Wheat

Fusarium crown rot (FCR) is one of the most important diseases of wheat (Triticum aestivum L.). FCR is mainly caused by the fungal pathogens Fusarium culmorum and F. pseudograminearum. In order to identify new sources of resistance to FCR and to dissect the complexity of FCR resistance, a panel of 161 wheat accessions was phenotyped under growth room (GR) and greenhouse conditions (GH). Analysis of variance showed significant differences in crown rot development among wheat accessions and high heritability of genotype-environment interactions for GR (0.96) and GH (0.91). Mixed linear model analysis revealed seven novel quantitative trait loci (QTLs) linked to F. culmorum on chromosomes 2AL, 3AS, 4BS, 5BS, 5DS, 5DL and 6DS for GR and eight QTLs on chromosomes on 3AS, 3BS, 3DL, 4BS (2), 5BS, 6BS and 6BL for GH. Total phenotypic variances (R²) explained by the QTLs linked to GR and GH were 48% and 59%, respectively. In addition, five favorable epistasis interactions among the QTLs were detected for both GR and GH with and without main effects. Epistatic interaction contributed additional variation up to 21% under GR and 7% under GH indicating strong effects of environment on the expression of QTLs. Our results revealed FCR resistance responses in wheat to be complex and controlled by multiple QTLs.

Multi-Locus Genome Wide Association Mapping for Yield and Its Contributing Traits in Hexaploid Wheat under Different Water Regimes

Multi-locus genome wide association study was undertaken using a set of 320 diverse spring wheat accessions, which were each genotyped for 9,626 SNPs. The association panel was grown in replicated trials in four environments [two each in irrigated (IR) and rainfed (RF) environments], and phenotypic data were recorded for five traits including days to heading, days to maturity, plant height, thousand grain weight and grain yield. Forty-six significant marker-trait associations (MTAs) were identified for five traits. These included 20 MTAs in IR and 19 MTAs in RF environments; seven additional MTAs were common to both the environments. Five of these MTAs were co-localized with previously known QTL/MTAs and the remaining MTAs were novel and add to the existing knowledge. Three desirable haplotypes for agronomic traits, one for improvement in RF environment and two for improvement in IR environment were identified. Eighteen (18) promising candidate genes (CGs) involved in seven different biological activities were also identified. The expression profiles of four (Trehalose-6-Phosphate, APETALA2/Ethylene-responsive factor, DNA-binding One Zinc Finger and Gibberellin-dioxygenases) of the 18 genes showed that they were induced by drought stress in the wheat seedlings. The MTAs, haplotypes and CG-based markers may be used in marker-assisted breeding for drought tolerance in wheat.

Genome-wide association study of drought tolerance and biomass allocation in wheat

Genome wide association studies (GWAS) are important in discerning the genetic architecture of complex traits such as biomass allocation for improving drought tolerance and carbon sequestration potential of wheat. The objectives of this study were to deduce the population structure and marker-trait association for biomass traits in wheat under drought-stressed and non-stressed conditions. A 100-wheat (Triticum aestivum L.) genotype panel was phenotyped for days to heading (DTH), days to maturity (DTM), shoot biomass (SB), root biomass (RB), root to shoot ratio (RS) and grain yield (GY). The panel was sequenced using 15,600 single nucleotide polymorphism (SNPs) markers and subjected to genetic analysis using the compressed mixed linear model (CMLM) at false discovery rate (FDR < 0.05). Population structure analysis revealed six sub-clusters with high membership ancestry coefficient of ≤0.65 to their assigned sub-clusters. A total of 75 significant marker-trait associations (MTAs) were identified with a linkage disequilibrium threshold of 0.38 at 5cM. Thirty-seven of the MTAs were detected under drought-stressed condition and 48% were on the B genome, where most quantitative trait loci (QTLs) for RB, SB and GY were previously identified. There were seven pleiotropic markers for RB and SB that may facilitate simultaneous selection. Thirty-seven putative candidate genes were mined by gene annotation on the IWGSC RefSeq 1.1. The significant MTAs observed in this study will be useful in devising strategies for marker-assisted breeding for simultaneous improvement of drought tolerance and to enhance C sequestration capacity of wheat.

Breeding for dual-purpose wheat varieties using marker-trait associations for biomass yield and quality traits

This study demonstrates that an active breeding nursery with rotation can be used to identify marker-trait associations for biomass yield and quality parameters that are important for biorefinery purposes. Wheat straw is a valuable feedstock for bioethanol production, but due to the recalcitrant nature of lignocellulose, its efficient use in biorefineries is limited by its low digestibility and difficult conversion of structural carbohydrates into free sugars. A genome-wide association study (GWAS) was conducted to search for significant SNP markers that could be used in a breeding programme to improve the value of wheat straw in a biorefinery setting. As part of a 3-year breeding programme (2013-2016), 190 winter wheat lines were phenotyped for traits that affect the yield and quality of the harvested biomass. These traits included straw yield, plant height, lodging at three growth stages and Septoria tritici blotch (STB) susceptibility. Release of glucose, xylose and arabinose was determined after hydrothermal pretreatment and enzymatic hydrolysis of the straw. The lines were genotyped using 15 K SNP markers and 5552 SNP markers could be used after filtering. Heritability for all traits ranged from 0.02 to 0.74. GWASs were conducted using CMLM, SUPER and FarmCPU algorithms, to analyse which algorithm could detect the highest number of marker-trait associations (MTAs). Comparable tendencies were obtained from CMLM and FarmCPU, but FarmCPU produced the most significant results. MTAs were obtained for lodging, harvest index, plant height, STB, glucose, xylose and arabinose at a significance level of p < 9.01 × 10^-6. MTAs in chromosome 6A were observed for glucose, xylose and arabinose, and could be of importance for increasing sugar release for bioethanol production.

Genome-wide variation patterns between landraces and cultivars uncover divergent selection during modern wheat breeding

Genetic diversity, population structure, LD decay, and selective sweeps in 687 wheat accessions were analyzed, providing relevant guidelines to facilitate the use of the germplasm in wheat breeding. Common wheat (Triticum aestivum L.) is one of the most widely grown crops in the world. Landraces were subjected to strong human-mediated selection in developing high-yielding, good quality, and widely adapted cultivars. To investigate the genome-wide patterns of allelic variation, population structure and patterns of selective sweeps during modern wheat breeding, we tested 687 wheat accessions, including landraces (148) and cultivars (539) mainly from China and Pakistan in a wheat 90 K single nucleotide polymorphism array. Population structure analysis revealed that cultivars and landraces from China and Pakistan comprised three relatively independent genetic clusters. Cultivars displayed lower nucleotide diversity and a wider average LD decay across whole genome, indicating allelic erosion and a diversity bottleneck due to the modern breeding. Analysis of genetic differentiation between landraces and cultivars from China and Pakistan identified allelic variants subjected to selection during modern breeding. In total, 477 unique genome regions showed signatures of selection, where 109 were identified in both China and Pakistan germplasm. The majority of genomic regions were located in the B genome (225), followed by the A genome (175), and only 77 regions were located in the D genome. EigenGWAS was further used to identify key selection loci in modern wheat cultivars from China and Pakistan by comparing with global winter wheat and spring wheat diversity panels, respectively. A few known functional genes or loci found within these genome regions corresponded to known phenotypes for disease resistance, vernalization, quality, adaptability and yield-related traits. This study uncovered molecular footprints of modern wheat breeding and explained the genetic basis of polygenic adaptation in wheat. The results will be useful for understanding targets of modern wheat breeding, and in devising future breeding strategies to target beneficial alleles currently not pursued.

From landraces to improved cultivars: Assessment of genetic diversity and population structure of Mediterranean wheat using SNP markers

Assessment of genetic diversity and population structure in crops is essential for breeding and germplasm conservation. A collection of 354 bread wheat genotypes, including Mediterranean landraces and modern cultivars representative of the ones most widely grown in the Mediterranean Basin, were characterized with 11196 single nucleotide polymorphism (SNP) markers. Total genetic diversity (H_T) and polymorphic information content (PIC) were 0.36 and 0.30 respectively for both landraces and modern cultivars. Linkage disequilibrium for the modern cultivars was higher than for the landraces (0.18 and 0.12, respectively). Analysis of the genetic structure showed a clear geographical pattern for the landraces, which were clustered into three subpopulations (SPs) representing the western, northern and eastern Mediterranean, whereas the modern cultivars were structured according to the breeding programmes that developed them: CIMMYT/ICARDA, France/Italy, and Balkan/eastern European countries. The modern cultivars showed higher genetic differentiation (G_ST) and lower gene flow (0.1673 and 2.49, respectively) than the landraces (0.1198 and 3.67, respectively), indicating a better distinction between subpopulations. The maximum gene flow was observed between landraces from the northern Mediterranean SPs and the modern cultivars released mainly by French and Italian breeding programmes.

TaBT1, affecting starch synthesis and thousand kernel weight, underwent strong selection during wheat improvement

BRITTLE1 (BT1), responsible for unidirectional transmembrane transport of ADP-glucose, plays a pivotal role in starch synthesis of cereal grain. In this study, we isolated three TaBT1 homoeologous genes located on chromosomes 6A, 6B, and 6D in common wheat. TaBT1 was mainly expressed in developing grains, and knockdown of TaBT1 in common wheat produced a decrease in grain size, thousand kernel weight (TKW), and grain total starch content. High diversity was detected at the TaBT1-6B locus, with 24 polymorphic sites forming three haplotypes (Hap1, Hap2, and Hap3). Association analysis revealed that Hap1 and Hap2 were preferred haplotypes in modern breeding, for their significant correlations with higher TKW. Furthermore, β-glucuronidase (GUS) staining and enzyme activity assays in developing grains of transgenic rice with exogenous promoters indicated that the promoters of Hap1 and Hap2 showed stronger driving activity than that of Hap3. Evolutionary analysis revealed that BT1 underwent strong selection during wheat polyploidization. In addition, the frequency distribution of TaBT1-6B haplotypes revealed that Hap1 and Hap2 were preferred in global modern wheat cultivars. Our findings suggest that TaBT1 has an important effect on starch synthesis and TKW, and provide two valuable molecular markers for marker assisted selection (MAS) in wheat high-yield breeding.

TaZIM-A1 negatively regulates flowering time in common wheat (Triticum aestivum L.)

Flowering time is a critical determinant of regional adaptation for crops and has strong effects on crop yields. Here, we report that TaZIM-A1, an atypical GATA-like transcription factor, is a negative regulator of flowering in wheat. TaZIM-A1 possessed weak transcriptional repression activity, with its CCT domain functioning as the major inhibitory region. TaZIM-A1 expression exhibited a typical circadian oscillation pattern under various light regimes. Overexpression of TaZIM-A1 caused a delay in flowering time and a decrease in thousand-kernel weight (TKW) in wheat under long-day conditions. Moreover, TaZIM-A1 directly bound to the promoters of TaCO-1 and TaFT-1 and downregulated their expression. Sequence analysis of a collection of common wheat cultivars identified three and two haplotypes for TaZIM-A1 and TaZIM-B1, respectively. Association analysis revealed that TaZIM-A1-HapI/-HapIII and TaZIM-B1-HapI have undergone strong positive selection during modern wheat breeding, likely due to their association with earlier heading and higher TKW. Diagnostic markers were developed for these haplotypes that can be used for wheat cultivar improvement, via marker-assisted breeding.

DAWN: a resource for yielding insights into the diversity among wheat genomes

BACKGROUND

Democratising the growing body of whole genome sequencing data available for Triticum aestivum (bread wheat) has been impeded by the lack of a genome reference and the large computational requirements for analysing these data sets.

RESULTS

DAWN (Diversity Among Wheat geNomes) integrates data from the T. aestivum Chinese Spring (CS) IWGSC RefSeq v1.0 genome with public WGS and exome data from 17 and 62 accessions respectively, enabling researchers and breeders alike to investigate genotypic differences between wheat accessions at the level of whole chromosomes down to individual genes.

CONCLUSIONS

Using DAWN we show that it is possible to visualise small and large chromosomal deletions, identify haplotypes at a glance and spot the consequences of selective breeding. DAWN allows us to detect the break points of alien introgression segments brought into an accession when transferring desired genes. Furthermore, we can find possible explanations for reduced recombination in parts of a chromosome, we can predict regions with linkage drag, and also look at diversity in centromeric regions.

Genome-Wide Association Mapping of Seedling Heat Tolerance in Winter Wheat

Heat stress during the seedling stage of early-planted winter wheat (Triticum aestivum L.) is one of the most abiotic stresses of the crop restricting forage and grain production in the Southern Plains of the United States. To map quantitative trait loci (QTLs) and identify single-nucleotide polymorphism (SNP) markers associated with seedling heat tolerance, a genome-wide association mapping study (GWAS) was conducted using 200 diverse representative lines of the hard red winter wheat association mapping panel, which was established by the Triticeae Coordinated Agricultural Project (TCAP) and genotyped with the wheat iSelect 90K SNP array. The plants were initially planted under optimal temperature conditions in two growth chambers. At the three-leaf stage, one chamber was set to 40/35°C day/night as heat stress treatment, while the other chamber was kept at optimal temperature (25/20°C day/night) as control for 14 days. Data were collected on leaf chlorophyll content, shoot length, number of leaves per seedling, and seedling recovery after removal of heat stress treatment. Phenotypic variability for seedling heat tolerance among wheat lines was observed in this study. Using the mixed linear model (MLM), we detected multiple significant QTLs for seedling heat tolerance on different chromosomes. Some of the QTLs were detected on chromosomes that were previously reported to harbor QTLs for heat tolerance during the flowering stage of wheat. These results suggest that some heat tolerance QTLs are effective from the seedling to reproductive stages in wheat. However, new QTLs that have never been reported at the reproductive stage were found responding to seedling heat stress in the present study. Candidate gene analysis revealed high sequence similarities of some significant loci with candidate genes involved in plant stress responses including heat, drought, and salt stress. This study provides valuable information about the genetic basis of seedling heat tolerance in wheat. To the best of our knowledge, this is the first GWAS to map QTLs associated with seedling heat tolerance targeting early planting of dual-purpose winter wheat. The SNP markers identified in this study will be used for marker-assisted selection (MAS) of seedling heat tolerance during dual-purpose wheat breeding.

Genetic Diversity and Molecular Characterization of Iranian Durum Wheat Landraces (Triticum turgidum durum (Desf.) Husn.) Using DArT Markers

Durum wheat is grown mainly in rain-fed regions of Iran and the Mediterranean district under stressful conditions. Different environmental conditions and agricultural practices among ancient communities have led to the development of locally adapted genotypes known as landraces. Landraces are a valued source of genetic variety and show definite adaptation to local environmental conditions according to their home of origin. This study aimed to explore linkage disequilibrium (LD) analysis and the population structure and genetic diversity of Iranian durum wheat landraces. In this study, population structure and genome-wide LD were investigated in 129 durum landrace accessions using 1500 DArT markers. Both structure and discriminant analysis of principal components obviously subdivided the sample collection into seven distinct groups centered on key ancestors and regions of origin of the germplasm. Genetic diversity among the populations was primarily within population (68 vs. 32%). Mean LD values across the entire population sample decayed below r² of 0.11 after 1 cM. LD decay of genomes A and B of Iranian durum wheat landrace is approximately 2-3 cM (r² = 0.11) and approximately 0.5 cM (r² = 0.12), respectively. Altogether, low LD decay, a high number of subpopulations, and the high existence of genetic diversity among and within populations were characteristics of the Iranian durum landrace collection. Hence, the existing genetic diversity within the population can be associated with the very long evolutionary history of plants in Iran. The populations we studied are hence presented as a valuable resource that can be used in basic and applied research in durum wheat breeding.

Genome-Wide Association Study for Spot Blotch Resistance in Hard Winter Wheat

Spot blotch (SB) caused by Cochliobolus sativus (anamorph: Bipolaris sorokiniana) is an economically important disease of wheat worldwide. Under a severe epidemic condition, the disease can cause yield losses up to 70%. Previous approaches like bi-parental mapping for identifying SB resistant genes/QTLs exploited only a limited portion of the available genetic diversity with a lower capacity to detect polygenic traits, and had a lower marker density. In this study, we performed genome-wide association study (GWAS) for SB resistance in hard winter wheat association mapping panel (HWWAMP) of 294 genotypes. The HWWAMP was evaluated for response to B. sorokiniana (isolate SD40), and a range of reactions was observed with 10 resistant, 38 moderately resistant, 120 moderately resistant- moderately susceptible, 111 moderately susceptible, and 15 susceptible genotypes. GWAS using 15,590 high-quality SNPs and 294 genotypes we identified six QTLs (p = <0.001) on chromosomes 2D, 3A, 4A, 4B, 5A, and 7B that collectively explained 30% of the total variation for SB resistance. Highly associated SNPs were identified for all six QTLs, QSb.sdsu-2D.1 (SNP: Kukri_c31121_1460, R2 = 4%), QSb.sdsu-3A.1 (SNP: Excalibur_c46082_440, R2 = 4%), QSb.sdsu-4A.1 (SNP: IWA8475, R2 = 5.5%), QSb.sdsu-4B.1 (SNP: Excalibur_rep_c79414_306, R2 = 4%), QSb.sdsu-5A.1 (SNP: Kukri_rep_c104877_2166, R2 = 6%), and QSb.sdsu-7B.1 (SNP: TA005844-0160, R2 = 6%). Our study not only validates three (2D, 5A, and 7B) genomic regions identified in previous studies but also provides highly associated SNP markers for marker assisted selection. In addition, we identified three novel QTLs (QSb.sdsu-3A.1, QSb.sdsu-4A.1, and QSb.sdsu-4B.1) for SB resistance in wheat. Gene annotation analysis of the candidate regions identified nine NBS-LRR and 38 other plant defense-related protein families across multiple QTLs, and these could be used for fine mapping and further characterization of SB resistance in wheat. Comparative analysis with barley indicated the SB resistance locus on wheat chromosomes 2D, 3A, 5A, and 7B identified in our study are syntenic to the previously identified SB resistance locus on chromosomes 2H, 3H, 5H, and 7H in barley. The 10 highly resistant genotypes and SNP markers identified in our study could be very useful resources for breeding of SB resistance in wheat.

Alternative splicing results in a lack of starch synthase IIa-D in Chinese wheat landrace

We evaluated the SGP-1 protein composition of 368 Chinese wheat landraces using SDS-PAGE. The SGP-D1 null type was identified in three accessions (Xiaoqingmang, Pushanbamai, and P119). An 18-bp deletion and 9-bp variation were found at the junction region of the 7th intron and 8th exon, leading to deletion of the intron-exon junction recognition site AG when aligned the 8261-bp DNA sequence of TaSSIIa-D in Pushanbamai with that of Chinese Spring. Four cDNA types with mis-spliced isoforms were subsequently detected through amplification of TaSSIIa-D cDNAs. Among these, nine type II cDNAs with a 16-bp deletion in the 8th exon were detected, indicating that the major transcriptional pattern of TaSSIIa in Pushanbamai is type II. In the type IV cDNA, a 97-bp sequence remains undeleted in the end of the 5th exon. The amylose content in Pushanbamai was significantly higher than that in all control lines under field conditions, which suggested that deletion of SGP-D1 has an efficient impact on amylose content. As the TaSSIIa gene plays an important role in regulating the content of amylose, it is anticipated that these natural variants of TaSSIIa-D will provide useful resources for quality improvement in wheat.

Genome-Wide Association Study for Spot Blotch Resistance in Hard Winter Wheat

Spot blotch (SB) caused by Cochliobolus sativus (anamorph: Bipolaris sorokiniana) is an economically important disease of wheat worldwide. Under a severe epidemic condition, the disease can cause yield losses up to 70%. Previous approaches like bi-parental mapping for identifying SB resistant genes/QTLs exploited only a limited portion of the available genetic diversity with a lower capacity to detect polygenic traits, and had a lower marker density. In this study, we performed genome-wide association study (GWAS) for SB resistance in hard winter wheat association mapping panel (HWWAMP) of 294 genotypes. The HWWAMP was evaluated for response to B. sorokiniana (isolate SD40), and a range of reactions was observed with 10 resistant, 38 moderately resistant, 120 moderately resistant- moderately susceptible, 111 moderately susceptible, and 15 susceptible genotypes. GWAS using 15,590 high-quality SNPs and 294 genotypes we identified six QTLs (p = <0.001) on chromosomes 2D, 3A, 4A, 4B, 5A, and 7B that collectively explained 30% of the total variation for SB resistance. Highly associated SNPs were identified for all six QTLs, QSb.sdsu-2D.1 (SNP: Kukri_c31121_1460, R² = 4%), QSb.sdsu-3A.1 (SNP: Excalibur_c46082_440, R² = 4%), QSb.sdsu-4A.1 (SNP: IWA8475, R² = 5.5%), QSb.sdsu-4B.1 (SNP: Excalibur_rep_c79414_306, R² = 4%), QSb.sdsu-5A.1 (SNP: Kukri_rep_c104877_2166, R² = 6%), and QSb.sdsu-7B.1 (SNP: TA005844-0160, R² = 6%). Our study not only validates three (2D, 5A, and 7B) genomic regions identified in previous studies but also provides highly associated SNP markers for marker assisted selection. In addition, we identified three novel QTLs (QSb.sdsu-3A.1, QSb.sdsu-4A.1, and QSb.sdsu-4B.1) for SB resistance in wheat. Gene annotation analysis of the candidate regions identified nine NBS-LRR and 38 other plant defense-related protein families across multiple QTLs, and these could be used for fine mapping and further characterization of SB resistance in wheat. Comparative analysis with barley indicated the SB resistance locus on wheat chromosomes 2D, 3A, 5A, and 7B identified in our study are syntenic to the previously identified SB resistance locus on chromosomes 2H, 3H, 5H, and 7H in barley. The 10 highly resistant genotypes and SNP markers identified in our study could be very useful resources for breeding of SB resistance in wheat.

Increased genomic prediction accuracy in wheat breeding using a large Australian panel

KEY MESSAGE

Genomic prediction accuracy within a large panel was found to be substantially higher than that previously observed in smaller populations, and also higher than QTL-based prediction. In recent years, genomic selection for wheat breeding has been widely studied, but this has typically been restricted to population sizes under 1000 individuals. To assess its efficacy in germplasm representative of commercial breeding programmes, we used a panel of 10,375 Australian wheat breeding lines to investigate the accuracy of genomic prediction for grain yield, physical grain quality and other physiological traits. To achieve this, the complete panel was phenotyped in a dedicated field trial and genotyped using a custom AxiomTM Affymetrix SNP array. A high-quality consensus map was also constructed, allowing the linkage disequilibrium present in the germplasm to be investigated. Using the complete SNP array, genomic prediction accuracies were found to be substantially higher than those previously observed in smaller populations and also more accurate compared to prediction approaches using a finite number of selected quantitative trait loci. Multi-trait genetic correlations were also assessed at an additive and residual genetic level, identifying a negative genetic correlation between grain yield and protein as well as a positive genetic correlation between grain size and test weight.

ADP-glucose pyrophosphorylase genes, associated with kernel weight, underwent selection during wheat domestication and breeding

ADP-glucose pyrophosphorylase, comprising two small subunits and two large subunits, is considered a key enzyme in the endosperm starch synthesis pathway in wheat (Triticum aestivum L.). Two genes, TaAGP-S1-7A and TaAGP-L-1B, were investigated in this study. Haplotypes of these genes were associated with thousand kernel weight (TKW) in different populations. Mean TKWs of favoured haplotypes were significantly higher than those of nonfavoured ones. Two molecular markers developed to distinguish these haplotypes could be used in molecular breeding. Frequencies of favoured haplotypes were dramatically increased in cultivars released in China after the 1940s. These favoured haplotypes were also positively selected in six major wheat production regions globally. Selection of AGP-S1 and AGP-L-1B in wheat mainly occurred during and after hexaploidization. Strong additive effects of the favoured haplotypes of with other genes for starch synthesis were also detected in different populations.

Genome-wide association mapping of black point reaction in common wheat (Triticum aestivum L.)

BACKGROUND

Black point is a serious threat to wheat production and can be managed by host resistance. Marker-assisted selection (MAS) has the potential to accelerate genetic improvement of black point resistance in wheat breeding. We performed a genome-wide association study (GWAS) using the high-density wheat 90 K and 660 K single nucleotide polymorphism (SNP) assays to better understand the genetic basis of black point resistance and identify associated molecular markers.

RESULTS

Black point reactions were evaluated in 166 elite wheat cultivars in five environments. Twenty-five unique loci were identified on chromosomes 2A, 2B, 3A, 3B (2), 3D, 4B (2), 5A (3), 5B (3), 6A, 6B, 6D, 7A (5), 7B and 7D (2), respectively, explaining phenotypic variation ranging from 7.9 to 18.0%. The highest number of loci was detected in the A genome (11), followed by the B (10) and D (4) genomes. Among these, 13 were identified in two or more environments. Seven loci coincided with known genes or quantitative trait locus (QTL), whereas the other 18 were potentially novel loci. Linear regression showed a clear dependence of black point scores on the number of favorable alleles, suggesting that QTL pyramiding will be an effective approach to increase resistance. In silico analysis of sequences of resistance-associated SNPs identified 6 genes possibly involved in oxidase, signal transduction and stress resistance as candidate genes involved in black point reaction.

CONCLUSION

SNP markers significantly associated with black point resistance and accessions with a larger number of resistance alleles can be used to further enhance black point resistance in breeding. This study provides new insights into the genetic architecture of black point reaction.