Optical maps refine the bread wheat Triticum aestivum cv. Chinese Spring genome assembly

Wheat Omics: Advancements and Opportunities

Plant omics, which includes genomics, transcriptomics, metabolomics and proteomics, has played a remarkable role in the discovery of new genes and biomolecules that can be deployed for crop improvement. In wheat, great insights have been gleaned from the utilization of diverse omics approaches for both qualitative and quantitative traits. Especially, a combination of omics approaches has led to significant advances in gene discovery and pathway investigations and in deciphering the essential components of stress responses and yields. Recently, a Wheat Omics database has been developed for wheat which could be used by scientists for further accelerating functional genomics studies. In this review, we have discussed various omics technologies and platforms that have been used in wheat to enhance the understanding of the stress biology of the crop and the molecular mechanisms underlying stress tolerance.

Multiomics Based Association Mapping in Wheat Reveals Genetic Architecture of Quality and Allergenic Related Proteins

Wheat is an important staple crop since its proteins contribute to human and animal nutrition and are important for its end-use quality. However, wheat proteins can also cause adverse human reactions for a large number of people. We performed a genome wide association study (GWAS) on 114 proteins quantified by LC-MS-based proteomics and expressed in an environmentally stable manner in 148 wheat cultivars with a heritability > 0.6. For 54 proteins, we detected quantitative trait loci (QTL) that exceeded the Bonferroni-corrected significance threshold and explained 17.3−84.5% of the genotypic variance. Proteins in the same family often clustered at a very close chromosomal position or the potential homeolog. Major QTLs were found for four well-known glutenin and gliadin subunits, and the QTL segregation pattern in the protein encoding the high molecular weight glutenin subunit Dx5 could be confirmed by SDS gel-electrophoresis. For nine potential allergenic proteins, large QTLs could be identified, and their measured allele frequencies open the possibility to select for low protein abundance by markers as long as their relevance for human health has been conclusively demonstrated. A potential allergen was introduced in the beginning of 1980s that may be linked to the cluster of resistance genes introgressed on chromosome 2AS from Triticum ventricosum. The reported sequence information for the 54 major QTLs can be used to design efficient markers for future wheat breeding.

Genome-Wide Association Study of Arabinoxylan Content from a 562 Hexaploid Wheat Collection

The selection of wheat varieties with high arabinoxylan (AX) levels could effectively improve the daily consumption of dietary fiber. However, studies on the selection of markers for AX levels are scarce. This study analyzed AX levels in 562 wheat genotypes collected from 46 countries using a GWAS with the BLINK model in the GAPIT3. Wheat genotypes were classified into eight subpopulations that exhibited high genetic differentiation based on 31,926 SNP loci. Eight candidate genes were identified, among which those encoding F-box domain-containing proteins, disease resistance protein RPM1, and bZIP transcription factor 29 highly correlated with AX levels. The AX level was higher in the adenine allele than in the guanine alleles of these genes in the wheat collection. In addition, the AX level was approximately 10% higher in 3 adenine combinations than 2 guanine, 1 adenine, and 3 guanine combinations in genotypes of three genes (F-box domain-containing proteins, RPM1, and bZIP transcription factor 29). The adenine allele, present in 97.46% of AX-95086356 SNP, exhibited a high correlation with AX levels following classification by country. Notably, the East Asian wheat genotypes contain high adenine alleles in three genes. These results highlight the potential of these three SNPs to serve as selectable markers for high AX content.

RNA-Seq Data Processing in Plant-Pathogen Interaction System: A Case Study

In RNA-seq data processing, short reads are usually aligned from one species against its own genome sequence; however, in plant-pathogen interaction systems, reads from both host and pathogen samples are blended together. In contrast with single-genome analyses, both pathogen and host reference genomes are involved in the alignment process. In such circumstances, the order in which the alignment is carried out, whether the host or pathogen is aligned first, or if both genomes are aligned simultaneously, influences the read counts of certain genes. This is a problem, especially at advanced infection stages. It is crucial to have an appropriate strategy for aligning the reads to their respective genomes, yet the existing strategies of either sequential or parallel alignment become problematic when mapping mixed reads to their corresponding reference genomes. The challenge lies in the determination of which reads belong to which species, especially when homology exists between the host and pathogen genomes. This chapter proposes a combo-genome alignment strategy, which was compared with existing alignment scenarios. Simulation results demonstrated that the degree of discrepancy in the results is correlated with phylogenetic distance of the two species in the mixture which was attributable to the extent of homology between the two genomes involved. This correlation was also found in the analysis using two real RNA-seq datasets of Fusarium-challenged wheat plants. Comparisons of the three RNA-seq processing strategies on three simulation datasets and two real Fusarium-infected wheat datasets showed that an alignment to a combo-genome, consisting of both host and pathogen genomes, improves mapping quality as compared to sequential alignment procedures.

Agronomic Trait Analysis and Genetic Mapping of a New Wheat Semidwarf Gene Rht-SN33d

Plant height is a key agronomic trait that is closely to the plant morphology and lodging resistance in wheat. However, at present, the few dwarf genes widely used in wheat breeding have narrowed wheat genetic diversity. In this study, we selected a semi-dwarf wheat mutant dwarf33 that exhibits decreased plant height with little serious negative impact on other agronomic traits. Genetic analysis and mutant gene mapping indicated that dwarf33 contains a new recessive semi-dwarf gene Rht-SN33d, which was mapped into ~1.3 Mb interval on the 3DL chromosome. The gibberellin metabolism-related gene TraesCS3D02G542800, which encodes gibberellin 2-beta-dioxygenase, is considered a potential candidate gene of Rht-SN33d. Rht-SN33d reduced plant height by approximately 22.4% in mutant dwarf33. Further study revealed that shorter stem cell length may be the main factor causing plant height decrease. In addition, the coleoptile length of dwarf33 was just 9.3% shorter than that of wild-type Shaannong33. These results will help to expand our understanding of new mechanisms of wheat height regulation, and obtain new germplasm for wheat improvement.

A community resource to mass explore the wheat grain proteome and its application to the late-maturity alpha-amylase (LMA) problem

BACKGROUND

Late-maturity alpha-amylase (LMA) is a wheat genetic defect causing the synthesis of high isoelectric point alpha-amylase following a temperature shock during mid-grain development or prolonged cold throughout grain development, both leading to starch degradation. While the physiology is well understood, the biochemical mechanisms involved in grain LMA response remain unclear. We have applied high-throughput proteomics to 4,061 wheat flours displaying a range of LMA activities. Using an array of statistical analyses to select LMA-responsive biomarkers, we have mined them using a suite of tools applicable to wheat proteins.

RESULTS

We observed that LMA-affected grains activated their primary metabolisms such as glycolysis and gluconeogenesis; TCA cycle, along with DNA- and RNA- binding mechanisms; and protein translation. This logically transitioned to protein folding activities driven by chaperones and protein disulfide isomerase, as well as protein assembly via dimerisation and complexing. The secondary metabolism was also mobilized with the upregulation of phytohormones and chemical and defence responses. LMA further invoked cellular structures, including ribosomes, microtubules, and chromatin. Finally, and unsurprisingly, LMA expression greatly impacted grain storage proteins, as well as starch and other carbohydrates, with the upregulation of alpha-gliadins and starch metabolism, whereas LMW glutenin, stachyose, sucrose, UDP-galactose, and UDP-glucose were downregulated.

CONCLUSIONS

To our knowledge, this is not only the first proteomics study tackling the wheat LMA issue but also the largest plant-based proteomics study published to date. Logistics, technicalities, requirements, and bottlenecks of such an ambitious large-scale high-throughput proteomics experiment along with the challenges associated with big data analyses are discussed.

Chromosome-Scale Genome Assembly and Annotation of Allotetraploid Annual Bluegrass (Poa annua L.)

Poa annua L. is a globally distributed grass with economic and horticultural significance as a weed and as a turfgrass. This dual significance, and its phenotypic plasticity and ecological adaptation, have made P. annua an intriguing plant for genetic and evolutionary studies. Because of the lack of genomic resources and its allotetraploid (2n = 4x = 28) nature, a reference genome sequence would be a valuable asset to better understand the significance and polyploid origin of P. annua. Here we report a genome assembly with scaffolds representing the 14 haploid chromosomes that are 1.78 Gb in length with an N50 of 112 Mb and 96.7% of BUSCO orthologs. Seventy percent of the genome was identified as repetitive elements, 91.0% of which were Copia- or Gypsy-like long-terminal repeats. The genome was annotated with 76,420 genes spanning 13.3% of the 14 chromosomes. The two subgenomes originating from Poa infirma (Knuth) and Poa supina (Schrad) were sufficiently divergent to be distinguishable but syntenic in sequence and annotation with repetitive elements contributing to the expansion of the P. infirma subgenome.

Characteristics and Expression Analysis of Invertase Gene Family in Common Wheat (Triticum aestivum L.)

Invertase (INV) irreversibly catalyzes the conversion of sucrose into glucose and fructose, playing important role in plant development and stress tolerance. However, the functions of INV genes in wheat have been less studied. In this study, a total of 126 TaINV genes were identified using a genome-wide search method, which could be classified into five classes (TaCWI-α, TaCWI-β, TaCI-α, TaCI-β, and TaVI) based on phylogenetic relationship. A total of 101 TaINVs were collinear with their ancestors in the synteny analysis, and we speculated that polyploidy events were the main force in the expansion of the TaINV gene family. Compared with TaCI, TaCWI and TaVI are more similar in gene structure and protein properties. Transcriptome sequencing analysis showed that TaINVs expressed in multiple tissues with different expression levels. Among 19 tissue-specific expressed TaINVs, 12 TaINVs showed grain-specific expression pattern and might play an important role in wheat grain development. In addition, qRT-PCR results further confirmed that TaCWI50 and TaVI27 show different expression in grain weight NILs. Our results demonstrated that the high expression of TaCWI50 and TaVI27 may be associated with a larger TGW phenotype. This work provides the foundations for understanding the grain development mechanism.

Gene Expression Datasets for Two Versions of the Saccharum spontaneum AP85-441 Genome

Sugarcane is a species of tall grass with high biomass and sucrose production, and the world’s largest crop by production quantity. Its evolutionary environment adaptation and anthropogenic breeding response have resulted in a complex autopolyploid genome. Few efforts have been reported in the literature to document this organism’s gene co-expression and annotation, and, when available, use different gene identifiers that cannot be easily associated across studies. This data descriptor paper presents a dataset that consolidates expression matrices of two Saccharum spontaneum AP85-441 genome versions and an algorithm implemented in Python to mechanically obtain this dataset. The data are processed from the allele-level information of the two sources, with BLASTn used bidirectionally to suggest feasible mappings between the two sets of alleles, and a graph-matching optimization algorithm to maximize global identity and uniqueness of genes. Association tables are used to consolidate the expression values from alleles to genes. The contributed expression matrices comprise 96 experiments and 109,050 and 35,516 from the two genome versions. They can represent significant computational cost reduction for further research on, e.g., sugarcane co-expression network generation, functional annotation prediction, and stress-specific gene identification.

Genome-wide identification of cystathionine beta synthase genes in wheat and its relationship with anther male sterility under heat stress

Cystathionine beta synthase (CBS) domains containing proteins (CDCPs) plays an important role in plant development through regulation of the thioredoxin system, as well as its ability to respond to biotic and abiotic stress conditions. Despite this, no systematic study has examined the wheat CBS gene family and its relation to high temperature-induced male sterility. In this study, 66 CBS family members were identified in the wheat genome, and their gene or protein sequences were used for subsequent analysis. The TaCBS gene family was found to be unevenly distributed on 21 chromosomes, and they were classified into four subgroups according to their gene structure and phylogeny. The results of collinearity analysis showed that there were 25 shared orthologous genes between wheat, rice and Brachypodium distachyon, and one shared orthologous gene between wheat, millet and barley. The cis-regulatory elements of the TaCBS were related to JA, IAA, MYB, etc. GO and KEGG pathway analysis identified these TaCBS genes to be associated with pollination, reproduction, and signaling and cellular processes, respectively. A heatmap of wheat plants based on transcriptome data showed that TaCBS genes were expressed to a higher extent in spikelets relative to other tissues. In addition, 29 putative tae-miRNAs were identified, targeting 41 TaCBS genes. Moreover, qRT-PCR validation of six TaCBS genes indicated their critical role in anther development, as five of them were expressed at lower levels in heat-stressed male sterile anthers than in Normal anthers. Together with anther phenotypes, paraffin sections, starch potassium iodide staining, and qRT-PCR data, we hypothesized that the TaCBS gene has a very important connection with the heat-stressed sterility process in wheat, and these data provide a basis for further insight into their relationship.

Wheat transcriptomic responses to extended feeding by wheat curl mites

The economic importance of wheat and its contribution to human and livestock diets has been already demonstrated. However, wheat production is impacted by pests that induce yield reductions. Among these pests, wheat curl mite (WCM, Aceria tosichella Keifer) impacts wheat all around the world. WCM are tiny pests that feed within the whorl of developing leaves, and their feeding causes leaf curling by preventing them from unfurling. The curling of the leaves provides a protective niche for the WCM. Additionally, WCM are also the vector of serious viruses in wheat. Little is known regarding the impact of the WCM on wheat transcriptome, and to date, only one article has been published describing the wheat transcriptomic changes after 1 day of WCM feeding. To better understand the wheat transcriptome variation after extended feeding by WCM [10 days post infestation (dpi)], we used an RNA-seq approach. We collected WCM-infested and uninfested leaves from two wheat cultivars: Byrd (WCM resistant) and Settler CL (WCM susceptible) at 10 dpi. Our transcriptomic analysis revealed the common and specific transcriptomic variations in WCM resistant and susceptible wheat cultivars, chromosome 3D specific location of the differentially expressed genes with functions involved in defense and stress response, and also identified the gene functions related to lipid signaling and membrane integrity, and phytohormone pathways potentially contributing to WCM resistance. Collectively, our study provides important insights on wheat defense mechanisms against WCM after extended feeding.

Automated Quantitative Measurement of Yellow Halos Suggests Activity of Necrotrophic Effectors in Septoria tritici Blotch

Many necrotrophic plant pathogens utilize host-selective toxins or necrotrophic effectors during the infection process. We hypothesized that the chlorotic yellow halos frequently observed around necrotic lesions caused by the wheat pathogen Zymoseptoria tritici could result from the activity of necrotrophic effectors interacting with the products of toxin sensitivity genes. As an initial step toward testing this hypothesis, we developed an automated image analysis (AIA) workflow that could quantify the degree of yellow halo formation occurring in wheat leaves naturally infected by a highly diverse pathogen population under field conditions. This AIA based on statistical learning was applied to more than 10,000 naturally infected leaves collected from 335 wheat cultivars grown in a replicated field experiment. We estimated a high heritability (h² = 0.71) for the degree of yellow halo formation, suggesting that this quantitative trait has a significant genetic component. Using genome-wide association mapping, we identified six chromosome segments significantly associated with the yellow halo phenotype. Most of these segments contained candidate genes associated with targets of necrotrophic effectors in other necrotrophic pathogens. Our findings conform with the hypothesis that toxin sensitivity genes could account for a significant fraction of the observed variation in quantitative resistance to Septoria tritici blotch. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Cytogenetic and genomic characterization of a novel tall wheatgrass-derived Fhb7 allele integrated into wheat B genome

We identified and integrated the novel FHB-resistant Fhb7^The2 allele into wheat B genome and made it usable in both common and durum wheat breeding programs without yellow flour linkage drag. A novel tall wheatgrass-derived (Thinopyrum elongatum, genome EE) Fhb7 allele, designated Fhb7^The2, was identified and integrated into the wheat B genome through a small 7B-7E translocation (7BS·7BL-7EL) involving the terminal regions of the long arms. Fhb7^The2 conditions significant Type II resistance to Fusarium head blight (FHB) in wheat. Integration of Fhb7^The2 into the wheat B genome makes this wild species-derived FHB resistance gene usable for breeding in both common and durum wheat. By contrast, other Fhb7 introgression lines involving wheat chromosome 7D can be utilized only in common wheat breeding programs, not in durum wheat. Additionally, we found that Fhb7^The2 does not have the linkage drag of the yellow flour pigment gene that is tightly linked to the decaploid Th. ponticum-derived Fhb7 allele Fhb7^Thp. This will further improve the utility of Fhb7^The2 in wheat breeding. DNA sequence analysis identified 12 single nucleotide polymorphisms (SNPs) in Fhb7^The2, Fhb7^Thp, and another Th. elongatum-derived Fhb7 allele Fhb7^The1, which led to seven amino acid conversions in Fhb7^The2, Fhb7^Thp, and Fhb7^The1, respectively. However, no significant variation was observed in their predicted protein configuration as a glutathione transferase. Diagnostic DNA markers were developed specifically for Fhb7^The2. The 7EL segment containing Fhb7^The2 in the translocation chromosome 7BS·7BL-7EL exhibited a monogenic inheritance pattern in the wheat genetic background. This will enhance the efficacy of marker-assisted selection for Fhb7^The2 introgression, pyramiding, and deployment in wheat germplasm and varieties.

Mapping and validation of major and stable QTL for flag leaf size from tetraploid wheat

The flag leaf is an important photosynthetic organ of wheat (Triticum aestivum L.). Appropriate flag leaf size can effectively increase grain yield. In this study, a tetraploid wheat population of recombinant inbred lines (RILs) and a genetic map constructed based on a wheat 55K single-nucleotide polymorphism (SNP) array were used to identify quantitative trait loci (QTL) for flag leaf length (FLL), flag leaf width (FLW), flag leaf area (FLA), and the flag leaf length/width ratio (FLR). A novel and major interval flanked by markers AX-111633224 and AX-109317229 was identified. This interval includes QTL for FLL (QFll.sau-AM-4B.2), for FLW (QFlw.sau-AM-4B.4), for FLA (QFla.sau-AM-4B), and for FLR (QFlr.sau-AM-4B). Based on the genotypes of the closely linked KASP (Kompetitive allele-specific polymerase chain reaction [PCR]) marker (KASP-AX-108756198), QFlw.sau-AM-4B.4 and QFla.sau-AM-4B were successfully verified in two F₃ populations with different genetic backgrounds. Genetic associations between flag leaf-related traits and other agronomic traits were detected and analyzed. Four genes in this interval were likely involved in the growth and development of the flag leaf size. In conclusion, this study provides clues for excavating genes related to flag leaf size and breeding variety with ideal plant structure.

Nested association mapping-based GWAS for grain yield and related traits in wheat grown under diverse Australian environments

Utilising a nested association mapping (NAM) population-based GWAS, 98 stable marker-trait associations with 127 alleles unique to the exotic parents were detected for grain yield and related traits in wheat. Grain yield, thousand-grain weight, screenings and hectolitre weight are important wheat yield traits. An understanding of their genetic basis is crucial for improving grain yield in breeding programmes. Nested association mapping (NAM) populations are useful resources for the dissection of the genetic basis of complex traits such as grain yield and related traits in wheat. Coupled with phenotypic data collected from multiple environments, NAM populations have the power to detect quantitative trait loci and their multiple alleles, providing germplasm that can be incorporated into breeding programmes. In this study, we evaluated a large-scale wheat NAM population with two recurrent parents in unbalanced trials in nine diverse Australian field environments over three years. By applying a single-stage factor analytical linear mixed model (FALMM) to the NAM multi-environment trials (MET) data and conducting a genome-wide association study (GWAS), we detected 98 stable marker-trait associations (MTAs) with their multiple alleles. 74 MTAs had 127 alleles that were derived from the exotic parents and were absent in either of the two recurrent parents. The exotic alleles had favourable effects on 46 MTAs of the 74 MTAs, for grain yield, thousand-grain weight, screenings and hectolitre weight. Two NAM RILs with consistently high yield in multiple environments were also identified, highlighting the potential of the NAM population in supporting plant breeding through provision of germplasm that can be readily incorporated into breeding programmes. The identified beneficial exotic alleles introgressed into the NAM population provide potential target alleles for the genetic improvement of wheat and further studies aimed at pinpointing the underlying genes.

Separation of the effects of two reduced height (Rht) genes and genomic background to select for less Fusarium head blight of short-strawed winter wheat (Triticum aestivum L.) varieties

FHB resistance shared pleiotropic loci with plant height and anther retention. Genomic prediction allows to select for genomic background reducing FHB susceptibility in the presence of the dwarfing allele Rht-D1b. With the high interest for semi-dwarf cultivars in wheat, finding locally adapted resistance sources against Fusarium head blight (FHB) and FHB-neutral reduced height (Rht) genes is of utmost relevance. In this study, 401 genotypes of European origin without/with dwarfing alleles of Rht-D1 and/or Rht24 were analysed across five environments on FHB severity and the morphological traits such as plant height (PH), anther retention (AR), number of spikelets per ear, ear length and ear density. Data were analysed by combined correlation and path analyses, association mapping and coupling single- and multi-trait genome-wide association studies (ST-GWAS and MT-GWAS, respectively) and genomic prediction (GP). All FHB data were corrected for flowering date or heading stage. High genotypic correlation (r_g = 0.74) and direct path effect (0.57) were detected between FHB severity and anther retention (AR). Moderate correlation (r_g = - 0.55) was found between FHB severity and plant height (PH) with a high indirect path via AR (- 0.31). Indirect selection for FHB resistance should concentrate on AR and PH. ST-GWAS identified 25 quantitative trait loci (QTL) for FHB severity, PH and AR, while MT-GWAS detected six QTL across chromosomes 2A, 4D, 5A, 6B and 7B conveying pleiotropic effects on the traits. Rht-D1b was associated with high AR and FHB susceptibility. Our study identified a promising positively acting pleiotropic QTL on chromosome 7B which can be utilized to improve FHB resistance while reducing PH and AR. Rht-D1b genotypes having a high resistance genomic background exhibited lower FHB severity and AR. The use of GP for estimating the genomic background was more effective than selection of GWAS-detected markers. We demonstrated that GP has a great potential and should be exploited by selecting for semi-dwarf winter wheat genotypes with higher FHB resistance due to their genomic background.

A draft chromosome-scale genome assembly of a commercial sugarcane

Sugarcane accounts for a large portion of the worlds sugar production. Modern commercial cultivars are complex hybrids of S. officinarum, S. spontaneum, and several other Saccharum species, resulting in an auto-allopolyploid with 8-12 copies of each chromosome. The current genome assembly gold standard is to generate a long read assembly followed by chromatin conformation capture sequencing to scaffold. We used the PacBio RSII and chromatin conformation capture sequencing to sequence and assemble the genome of a South East Asian commercial sugarcane cultivar, known as Khon Kaen 3. The Khon Kaen 3 genome assembled into 104,477 contigs totalling 7 Gb, which scaffolded into 56 pseudochromosomes containing 5.2 Gb of sequence. Genome annotation produced 242,406 genes from 30,927 orthogroups. Aligning the Khon Kaen 3 genome sequence to S. officinarum and S. spontaneum revealed a high level of apparent recombination, indicating a chimeric assembly. This assembly error is explained by high nucleotide identity between S. officinarum and S. spontaneum, where 91.8% of S. spontaneum aligns to S. officinarum at 94% identity. Thus, the subgenomes of commercial sugarcane are so similar that using short reads to correct long PacBio reads produced chimeric long reads. Future attempts to sequence sugarcane must take this information into account.

Genome-wide association study for resistance in bread wheat (Triticum aestivum L.) to stripe rust (Puccinia striiformis f. sp. tritici) races in Argentina

BACKGROUND

Wheat stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most devastating diseases of the wheat crop. It causes significant reductions in both grain yield and grain quality. In recent years, new and more virulent races have overcome many of the known resistance genes in Argentinian germplasm. In order to identify loci conferring resistance to the local races of Pst for effective utilization in future breeding programs, a genome-wide association study (GWAS) was performed using a collection of 245 bread wheat lines genotyped with 90 K SNPs.

RESULTS

To search for adult plant resistance (APR) the panel was evaluated for disease severity (DS) and area under disease progress curve (AUDPC) in field trials during two years under natural infection conditions. To look for seedling or all-stage resistance (ASR) the panel was evaluated to determine infection type (IT) under greenhouse conditions against two prevalent races in Argentina. The phenotypic data showed that the panel possessed enough genetic variability for searching for sources of resistance to Pst. Significant correlations between years were observed for Pst response in the field and high heritability values were found for DS (H² = 0.89) and AUDPC (H² = 0.93). Based on GWAS, eight markers associated with Pst resistance (FDR < 0.01) were identified, of these, five were associated with ASR (on chromosomes 1B, 2A, 3A and 5B) and three with APR (on chromosomes 3B and 7A). These markers explained between 2% and 32.62% of the phenotypic variation. Five of the markers corresponded with previously reported Yr genes/QTL, while the other three (QYr.Bce.1B.sd.1, QYr.Bce.3A.sd and QYr.Bce.3B.APR.2) might be novel resistance loci.

CONCLUSION

Our results revealed high genetic variation for resistance to Argentinian stripe rust races in the germplasm used here. It constitutes a very promising step towards the improvement of Pst resistance of bread wheat in Argentina. Also, the identification of new resistance loci would represent a substantial advance for diversifying the current set of resistance genes and to advance in the improvement of the durable resistance to the disease.

Leaf rolling in bread wheat (Triticum aestivum L.) is controlled by the upregulation of a pair of closely linked/duplicate zinc finger homeodomain class transcription factors during moisture stress conditions

Zinc finger-homeodomain (ZF-HDs) class IV transcriptional factors (TFs) is a plant-specific transcription factor and play a key role in stress responses, plant growth, development, and hormonal signaling. In this study, two new leaf rolling TFs genes, namely TaZHD1 and TaZHD10, were identified in wheat using comparative genomic analysis of the target region that carried a major QTL for leaf rolling identified through multi-environment phenotyping and high throughput genotyping of a RIL population. Structural and functional annotation of the candidate ZHD genes with its closest rice orthologs reflects the species-specific evolution and, undoubtedly, validates the notions of remote-distance homology concept. Meanwhile, the morphological analysis resulted in contrasting difference for leaf rolling in extreme RILs between parental lines HD2012 and NI5439 at booting and heading stages. Transcriptome-wide expression profiling revealed that TaZHD10 transcripts showed significantly higher expression levels than TaZHD1 in all leaf tissues upon drought stress. The relative expression of these genes was further validated by qRT-PCR analysis, which also showed consistent results across the studied genotypes at the booting and anthesis stage. The contrasting modulation of these genes under drought conditions and the available evidenced for its epigenetic behavior that might involve the regulation of metabolic and gene regulatory networks. Prediction of miRNAs resulted in five Tae-miRs that could be associated with RNAi mediated control of TaZHD1 and TaZHD10 putatively involved in the metabolic pathway controlling rolled leaf phenotype. Gene interaction network analysis indicated that TaZHD1 and TaZHD10 showed pleiotropic effects and might also involve other functions in wheat in addition to leaf rolling. Overall, the results increase our understanding of TaZHD genes and provide valuable information as robust candidate genes for future functional genomics research aiming for the breeding of wheat varieties tolerant to leaf rolling.

Heat Stress Tolerance 2 confers basal heat stress tolerance in allohexaploid wheat (Triticum aestivum L.)

Heat stress substantially reduces the yield potential of wheat (Triticum aestivum L.), one of the most widely cultivated staple crops, and greatly threatens global food security in the context of global warming. However, few studies have explored the heat stress tolerance (HST)-related genetic resources in wheat. Here, we identified and fine-mapped a wheat HST locus, TaHST2, which is indispensable for HST in both the vegetative and reproductive stages of the wheat life cycle. The studied pair of near isogenic lines (NILs) exhibited diverse morphologies under heat stress, based on which we mapped TaHST2 to a 485 kb interval on chromosome arm 4DS. Under heat stress, TaHST2 confers a superior conversion rate from soluble sugars to starch in wheat grains, resulting in faster grain filling and a higher yield potential. A further exploration of genetic resources indicated that TaHST2 underwent strong artificial selection during wheat domestication, suggesting it is an essential locus for basal HST in wheat. Our findings provide deeper insights into the genetic basis of wheat HST and might be useful for global efforts to breed heat-stress-tolerant cultivars.

Identification of Disease Resistance Parents and Genome-Wide Association Mapping of Resistance in Spring Wheat

The likelihood of success in developing modern cultivars depend on multiple factors, including the identification of suitable parents to initiate new crosses, and characterizations of genomic regions associated with target traits. The objectives of the present study were to (a) determine the best economic weights of four major wheat diseases (leaf spot, common bunt, leaf rust, and stripe rust) and grain yield for multi-trait restrictive linear phenotypic selection index (RLPSI), (b) select the top 10% cultivars and lines (hereafter referred as genotypes) with better resistance to combinations of the four diseases and acceptable grain yield as potential parents, and (c) map genomic regions associated with resistance to each disease using genome-wide association study (GWAS). A diversity panel of 196 spring wheat genotypes was evaluated for their reaction to stripe rust at eight environments, leaf rust at four environments, leaf spot at three environments, common bunt at two environments, and grain yield at five environments. The panel was genotyped with the Wheat 90K SNP array and a few KASP SNPs of which we used 23,342 markers for statistical analyses. The RLPSI analysis performed by restricting the expected genetic gain for yield displayed significant (p < 0.05) differences among the 3125 economic weights. Using the best four economic weights, a subset of 22 of the 196 genotypes were selected as potential parents with resistance to the four diseases and acceptable grain yield. GWAS identified 37 genomic regions, which included 12 for common bunt, 13 for leaf rust, 5 for stripe rust, and 7 for leaf spot. Each genomic region explained from 6.6 to 16.9% and together accounted for 39.4% of the stripe rust, 49.1% of the leaf spot, 94.0% of the leaf rust, and 97.9% of the common bunt phenotypic variance combined across all environments. Results from this study provide valuable information for wheat breeders selecting parental combinations for new crosses to develop improved germplasm with enhanced resistance to the four diseases as well as the physical positions of genomic regions that confer resistance, which facilitates direct comparisons for independent mapping studies in the future.

Computational genomics insights into cold acclimation in wheat

Development of cold acclimation in crops involves transcriptomic reprograming, metabolic shift, and physiological changes. Cold responses in transcriptome and lipid metabolism has been examined in separate studies for various crops. In this study, integrated computational approaches was employed to investigate the transcriptomics and lipidomics data associated with cold acclimation and vernalization in four wheat genotypes of distinct cold tolerance. Differential expression was investigated between cold treated and control samples and between the winter-habit and spring-habit wheat genotypes. Collectively, 12,676 differentially expressed genes (DEGs) were identified. Principal component analysis of these DEGs indicated that the first, second, and third principal components (PC1, PC2, and PC3) explained the variance in cold treatment, vernalization and cold hardiness, respectively. Differential expression feature extraction (DEFE) analysis revealed that the winter-habit wheat genotype Norstar had high number of unique DEGs (1884 up and 672 down) and 63 winter-habit genes, which were clearly distinctive from the 64 spring-habit genes based on PC1, PC2 and PC3. Correlation analysis revealed 64 cold hardy genes and 39 anti-hardy genes. Cold acclimation encompasses a wide spectrum of biological processes and the involved genes work cohesively as revealed through network propagation and collective association strength of local subnetworks. Integration of transcriptomics and lipidomics data revealed that the winter-habit genes, such as COR413-TM1, CIPKs and MYB20, together with the phosphatidylglycerol lipids, PG(34:3) and PG(36:6), played a pivotal role in cold acclimation and coordinated cohesively associated subnetworks to confer cold tolerance.

Identification of the major QTL QPm.cas-7D for adult plant resistance to wheat powdery mildew

Developing effective and durable host plant resistance is crucial for controlling powdery mildew, a devastating disease caused by Blumeria graminis f. sp. tritici (Bgt). In the present study, we dissected the genetic basis of the adult plant resistance to powdery mildew using a recombinant inbred line (RIL) composed of 176 F₉ RILs population derived from a cross between PuBing 3228 (P3228) and susceptible cultivar Gao 8901. P3228 exhibits stable adult-plant resistance to powdery mildew in the field over consecutive years. We identified two QTLs on chromosomes 7DS (QPm.cas-7D) and 1AL (QPm.cas-1A) contributed by P3228, and one QTL on 3DS (QPm.cas-3D) contributed by Gao 8901, which could explain 65.44%, 3.45%, and 2.18% of the phenotypic variances, respectively. By analyzing the annotated genes in the 1.168 Mb physical interval of the major QTL QPm.cas-7D, we locked a previously cloned adult-plant resistance gene Pm38 that was most probably the candidate gene of QPm.cas-7D. Sequence alignment analysis revealed that the candidate gene of QPm.cas-7D in P3228 was identical to the reported Pm38 sequence. Two haplotypes QPm-7D-R and QPm-7D-S were identified in the whole Pm38 genomic regions between P3228 and Gao 8901. To apply QPm.cas-7D in wheat breeding, we developed a kompetitive allele-specific PCR (KASP) marker Kasp5249 that is closely linked with these haplotypes. It is worth mentioning that the QPm-7D-R haplotype significantly decreased TKW and underwent negative selection for higher yields in China wheat breeding. In this study, we identified a major QTL QPm.cas-7D and revealed the relationship between its resistance and yield, which could be beneficial for further applications in wheat disease resistance and high-yield breeding.

Identification of the genes at S and Z reveals the molecular basis and evolution of grass self-incompatibility

Self-incompatibility (SI) is a feature of many flowering plants, whereby self-pollen is recognized and rejected by the stigma. In grasses (Poaceae), the genes controlling this phenomenon have not been fully elucidated. Grasses have a unique two-locus system, in which two independent genetic loci (S and Z) control self-recognition. S and Z are thought to have arisen from an ancient duplication, common to all grasses. With new chromosome-scale genome data, we examined the genes present at S- and Z-loci, firstly in ryegrass (Lolium perenne), and subsequently in ~20 other grass species. We found that two DUF247 genes and a short unstructured protein (SP/ZP) were present at both S- and Z- in all SI species, while in self-compatible species these genes were often lost or mutated. Expression data suggested that DUF247 genes acted as the male components and SP/ZP were the female components. Consistent with their role in distinguishing self- from non-self, all genes were hypervariable, although key secondary structure features were conserved, including the predicted N-terminal cleavage site of SP/ZP. The evolutionary history of these genes was probed, revealing that specificity groups at the Z-locus arose before the advent of various grass subfamilies/species, while specificity groups at the S-locus arose after the split of Panicoideae, Chloridoideae, Oryzoideae and Pooideae. Finally, we propose a model explaining how the proteins encoded at the S and Z loci might function to specify self-incompatibility.

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.

Polyploidy before and after domestication of crop species

Recent advances in the genomics of polyploid species answer some of the long-standing questions about the role of polyploidy in crop species. Here, we summarize the current literature to reexamine scenarios in which polyploidy played a role both before and after domestication. The prevalence of polyploidy can help to explain environmental robustness in agroecosystems. This review also clarifies the molecular basis of some agriculturally advantageous traits of polyploid crops, including yield increments in polyploid cotton via subfunctionalization, modification of a separated sexuality to selfing in polyploid persimmon via neofunctionalization, and transition to a selfing system via nonfunctionalization combined with epistatic interaction between duplicated S-loci. The rapid progress in genomics and genetics is discussed along with how this will facilitate functional studies of understudied polyploid crop species.

Genome-wide association study and genomic prediction of resistance to stripe rust in current Central and Northern European winter wheat germplasm

We found two loci on chromosomes 2BS and 6AL that significantly contribute to stripe rust resistance in current European winter wheat germplasm. Stripe or yellow rust, caused by the fungus Puccinia striiformis Westend f. sp. tritici, is one of the most destructive wheat diseases. Sustainable management of wheat stripe rust can be achieved through the deployment of rust resistant cultivars. To detect effective resistance loci for use in breeding programs, an association mapping panel of 230 winter wheat cultivars and breeding lines from Northern and Central Europe was employed. Genotyping with the Illumina® iSelect® 25 K Infinium® single nucleotide polymorphism (SNP) genotyping array yielded 8812 polymorphic markers. Structure analysis revealed two subpopulations with 92 Austrian breeding lines and cultivars, which were separated from the other 138 genotypes from Germany, Norway, Sweden, Denmark, Poland, and Switzerland. Genome-wide association study for adult plant stripe rust resistance identified 12 SNP markers on six wheat chromosomes which showed consistent effects over several testing environments. Among these, two marker loci on chromosomes 2BS (RAC875_c1226_652) and 6AL (Tdurum_contig29607_413) were highly predictive in three independent validation populations of 1065, 1001, and 175 breeding lines. Lines with the resistant haplotype at both loci were nearly free of stipe rust symptoms. By using mixed linear models with those markers as fixed effects, we could increase predictive ability in the three populations by 0.13-0.46 compared to a standard genomic best linear unbiased prediction approach. The obtained results facilitate an efficient selection for stripe rust resistance against the current pathogen population in the Northern and Central European winter wheat gene pool.

Chromosome painting reveals inter-chromosomal rearrangements and evolution of subgenome D of wheat

Aegilops species represent the most important gene pool for breeding bread wheat (Triticum aestivum). Thus, understanding the genome evolution, including chromosomal structural rearrangements and syntenic relationships among Aegilops species or between Aegilops and wheat, is important for both basic genome research and practical breeding applications. In the present study, we attempted to develop subgenome D-specific fluorescence in situ hybridization (FISH) probes by selecting D-specific oligonucleotides based on the reference genome of Chinese Spring. The oligo-based chromosome painting probes consisted of approximately 26 000 oligos per chromosome and their specificity was confirmed in both diploid and polyploid species containing the D subgenome. Two previously reported translocations involving two D chromosomes have been confirmed in wheat varieties and their derived lines. We demonstrate that the oligo painting probes can be used not only to identify the translocations involving D subgenome chromosomes, but also to determine the precise positions of chromosomal breakpoints. Chromosome painting of 56 accessions of Ae. tauschii from different origins led us to identify two novel translocations: a reciprocal 3D-7D translocation in two accessions and a complex 4D-5D-7D translocation in one accession. Painting probes were also used to analyze chromosomes from more diverse Aegilops species. These probes produced FISH signals in four different genomes. Chromosome rearrangements were identified in Aegilops umbellulata, Aegilops markgrafii, and Aegilops uniaristata, thus providing syntenic information that will be valuable for the application of these wild species in wheat breeding.

A GBS-based genome-wide association study reveals the genetic basis of salinity tolerance at the seedling stage in bread wheat (Triticum aestivum L.)

High salinity levels affect 20% of the cultivated area and 9%–34% of the irrigated agricultural land worldwide, ultimately leading to yield losses of crops. The current study evaluated seven salt tolerance-related traits at the seedling stage in a set of 138 pre-breeding lines (PBLs) and identified 63 highly significant marker-trait associations (MTAs) linked to salt tolerance. Different candidate genes were identified in in silico analysis, many of which were involved in various stress conditions in plants, including glycine-rich cell wall structural protein 1-like, metacaspase-1, glyceraldehyde-3-phosphate dehydrogenase GAPA1, and plastidial GAPA1. Some of these genes coded for structural protein and participated in cell wall structure, some were linked to programmed cell death, and others were reported to show abiotic stress response roles in wheat and other plants. In addition, using the Multi-Trait Genotype-Ideotype Distance Index (MGIDI) protocol, the best-performing lines under salt stress were identified. The SNPs identified in this study and the genotypes with favorable alleles provide an excellent source to impart salt tolerance in wheat.

Patterns of genetic variation and QTLs controlling grain traits in a collection of global wheat germplasm revealed by high-quality SNP markers

BACKGROUND

Establish a molecular breeding program involved assembling a diverse germplasm collection and generating accurate genotypes to characterize their genetic potential and associate them with agronomic traits. In this study, we acquired over eight hundred wheat accessions from international gene banks and assessed their genetic relatedness using high-quality SNP genotypes. Understanding the scope of genomic variation in this collection allows the breeders to utilize the genetic resources efficiently while improving wheat yield and quality.

RESULTS

A wheat diversity panel comprising 39 durum wheat, 60 spelt wheat, and 765 bread wheat accessions was genotyped on iSelect 90 K wheat SNP arrays. A total of 57,398 SNP markers were mapped to IWGSC RefSeq v2.1 assembly, over 30,000 polymorphic SNPs in the A, B, D genomes were used to analyze population structure and diversity, the results revealed the separation of the three species and the differentiation of CIMMYT improved breeding lines and landraces or widely grown cultivars. In addition, several chromosomal regions under selection were detected. A subset of 280 bread wheat accessions was evaluated for grain traits, including grain length, width, surface area, and color. Genome-wide association studies (GWAS) revealed that several chromosomal regions were significantly linked to known quantitative trait loci (QTL) controlling grain-related traits. One of the SNP peaks at the end of chromosome 7A was in strong linkage disequilibrium (LD) with WAPO-A1, a gene that governs yield components.

CONCLUSIONS

Here, the most updated and accurate physical positions of SNPs on 90 K genotyping array are provided for the first time. The diverse germplasm collection and associated genotypes are available for the wheat researchers to use in their molecular breeding program. We expect these resources to broaden the genetic basis of original breeding and pre-breeding materials and ultimately identify molecular markers associated with important agronomic traits which are evaluated in diverse environmental conditions.

The genome and gene editing system of sea barleygrass provide a novel platform for cereal domestication and stress tolerance studies

The tribe Triticeae provides important staple cereal crops and contains elite wild species with wide genetic diversity and high tolerance to abiotic stresses. Sea barleygrass (Hordeum marinum Huds.), a wild Triticeae species, thrives in saline marshlands and is well known for its high tolerance to salinity and waterlogging. Here, a 3.82-Gb high-quality reference genome of sea barleygrass is assembled de novo, with 3.69 Gb (96.8%) of its sequences anchored onto seven chromosomes. In total, 41 045 high-confidence (HC) genes are annotated by homology, de novo prediction, and transcriptome analysis. Phylogenetics, non-synonymous/synonymous mutation ratios (Ka/Ks), and transcriptomic and functional analyses provide genetic evidence for the divergence in morphology and salt tolerance among sea barleygrass, barley, and wheat. The large variation in post-domestication genes (e.g. IPA1 and MOC1) may cause interspecies differences in plant morphology. The extremely high salt tolerance of sea barleygrass is mainly attributed to low Na⁺ uptake and root-to-shoot translocation, which are mainly controlled by SOS1, HKT, and NHX transporters. Agrobacterium-mediated transformation and CRISPR/Cas9-mediated gene editing systems were developed for sea barleygrass to promote its utilization for exploration and functional studies of hub genes and for the genetic improvement of cereal crops.

Identification and validation of a locus for wheat maximum root length independent of parental reproductive environment

Maximum root length (MRL) plays an important role in the uptake of nutrients and resisting abiotic stresses. Understanding the genetic mechanism of root development is of great significance for genetic improvement of wheat. Previous studies have confirmed that parental reproductive environment (PRE) has a significant impact on growth and development of the next generation in the whole life cycle of a given plant. In this study, a recombinant inbred line population genotyped using the Wheat55K SNP array, was used to map quantitative trait loci (QTL) for wheat seedling MRL based on the harvested seeds from five different PREs. A total of 5 QTL located on chromosomes 3D and 7A were identified. Among them, QMrl.sicau-2SY-3D.2 located in a 4.0 cM interval on chromosome 3D was likely independent of PREs. QMrl.sicau-2SY-7A.2 was detected in two tests and probably influenced by PREs. The effect of QMrl.sicau-2SY-3D.2 was further validated using the tightly linked kompetitive allele specific PCR (KASP) marker, KASP-AX-111589572, in populations with different genetic backgrounds. Lines with a combination of positive alleles from QMrl.sicau-2SY-3D.2 and QMrl.sicau-2SY-7A.2 have significantly longer MRL. Furthermore, four genes (TraesCS3D03G0612000, TraesCS3D03G0608400, TraesCS3D03G0613600, and TraesCS3D03G0602400) mainly expressed in wheat root were predicted to be associated with root growth. Taken together, this study reports on a major QTL independent of PREs and lays a foundation for understanding the regulation mechanism of wheat MRL at the seedling stage.

GENESPACE tracks regions of interest and gene copy number variation across multiple genomes

The development of multiple chromosome-scale reference genome sequences in many taxonomic groups has yielded a high-resolution view of the patterns and processes of molecular evolution. Nonetheless, leveraging information across multiple genomes remains a significant challenge in nearly all eukaryotic systems. These challenges range from studying the evolution of chromosome structure, to finding candidate genes for quantitative trait loci, to testing hypotheses about speciation and adaptation. Here, we present GENESPACE, which addresses these challenges by integrating conserved gene order and orthology to define the expected physical position of all genes across multiple genomes. We demonstrate this utility by dissecting presence–absence, copy-number, and structural variation at three levels of biological organization: spanning 300 million years of vertebrate sex chromosome evolution, across the diversity of the Poaceae (grass) plant family, and among 26 maize cultivars. The methods to build and visualize syntenic orthology in the GENESPACE R package offer a significant addition to existing gene family and synteny programs, especially in polyploid, outbred, and other complex genomes.

The genome is the complete DNA sequence of an individual. It is a crucial foundation for many studies in medicine, agriculture, and conservation biology. Advances in genetics have made it possible to rapidly sequence, or read out, the genome of many organisms. For closely related species, scientists can then do detailed comparisons, revealing similar genes with a shared past or a common role, but comparing more distantly related organisms remains difficult.

One major challenge is that genes are often lost or duplicated over evolutionary time. One way to be more confident is to look at ‘synteny’, or how genes are organized or ordered within the genome. In some groups of species, synteny persists across millions of years of evolution. Combining sequence similarity with gene order could make comparisons between distantly related species more robust.

To do this, Lovell et al. developed GENESPACE, a software that links similarities between DNA sequences to the order of genes in a genome. This allows researchers to visualize and explore related DNA sequences and determine whether genes have been lost or duplicated. To demonstrate the value of GENESPACE, Lovell et al. explored evolution in vertebrates and flowering plants. The software was able to highlight the shared sequences between unique sex chromosomes in birds and mammals, and it was able to track the positions of genes important in the evolution of grass crops including maize, wheat, and rice.

Exploring the genetic code in this way could lead to a better understanding of the evolution of important sections of the genome. It might also allow scientists to find target genes for applications like crop improvement. Lovell et al. have designed the GENESPACE software to be easy for other scientists to use, allowing them to make graphics and perform analyses with few programming skills.

Identification, expression, and association analysis of calcineurin B-like protein–interacting protein kinase genes in peanut

Plants usually respond to the external environment by initiating a series of signal transduction processes mediated by protein kinases, especially calcineurin B-like protein–interacting protein kinases (CIPKs). In this study, 54 CIPKs were identified in the peanut genome, of which 26 were from cultivated species (named AhCIPKs) and 28 from two diploid progenitors (Arachis duranensis—AdCIPKs and Arachis ipaensis—AiCIPKs). Evolution analysis revealed that the 54 CIPKs were composed of two different evolutionary branches. The CIPK members were unevenly distributed at different chromosomes. Synteny analysis strongly indicated that whole-genome duplication (allopolyploidization) contributed to the expansion of CIPK. Comparative genomics analysis showed that there was only one common collinear CIPK pairs among peanut, Arabidopsis, rice, grape, and soybean. The prediction results of cis-acting elements showed that AhCIPKs, AdCIPKs, and AiCIPKs contained different proportions of transcription factor binding motifs involved in regulating plant growth, abiotic stress, plant hormones, and light response elements. Spatial expression profiles revealed that almost all AhCIPKs had tissue-specific expression patterns. Furthermore, association analysis identified one polymorphic site in AdCIPK12 (AhCIPK11), which was significantly associated with pod length, seed length, hundred seed weight, and shoot root ratio. Our results provide valuable information of CIPKs in peanut and facilitate better understanding of their biological functions.

Genome-wide association mapping identifies common bunt (Tilletia caries) resistance loci in bread wheat (Triticum aestivum) accessions of the USDA National Small Grains Collection

Association mapping and phenotypic analysis of a diversity panel of 238 bread wheat accessions highlights differences in resistance against common vs. dwarf bunt and identifies genotypes valuable for bi-parental crosses. Common bunt caused by Tilletia caries and T. laevis was successfully controlled by seed dressings with systemic fungicides for decades, but has become a renewed threat to wheat yield and quality in organic agriculture where such treatments are forbidden. As the most efficient way to address this problem is the use of resistant cultivars, this study aims to broaden the spectrum of resistance sources available for breeders by identifying resistance loci against common bunt in bread wheat accessions of the USDA National Small Grains Collection. We conducted three years of artificially inoculated field trials to assess common bunt infection levels in a diversity panel comprising 238 wheat accessions for which data on resistance against the closely related pathogen Tilletia controversa causing dwarf bunt was already available. Resistance levels against common bunt were higher compared to dwarf bunt with 99 accessions showing [Formula: see text] 1% incidence. Genome-wide association mapping identified six markers significantly associated with common bunt incidence in regions already known to confer resistance on chromosomes 1A and 1B and novel loci on 2B and 7A. Our results show that resistance against common and dwarf bunt is not necessarily controlled by the same loci but we identified twenty accessions with high resistance against both diseases. These represent valuable new resources for research and breeding programs since several bunt races have already been reported to overcome known resistance genes.

Genetic dissection of lutein content in common wheat via association and linkage mapping

Genetic architecture controlling grain lutein content of common wheat was investigated through an integration of genome-wide association study (GWAS) and linkage analysis. Putative candidate genes involved in carotenoid metabolism and regulation were identified, which provide a basis for gene cloning and development of nutrient-enriched wheat varieties through molecular breeding. Lutein, known as 'the eye vitamin', is an important component of wheat nutritional and end-use quality. However, the genetic manipulation of grain lutein content (LUC) in common wheat has not previously been well studied. Here, quantitative trait loci (QTL) associated with the LUC measured by high performance liquid chromatography (HPLC) were first identified by integrating a genome-wide association study (GWAS) and linkage mapping. A Chinese wheat mini-core collection (MCC) of 262 accessions and a doubled haploid (DH) population derived from Jinchun 7 and L1219 were genotyped using the 90K SNP array. A total of 124 significant marker-trait associations (MTAs) on all 21 wheat chromosomes except for 1A, 4D, and 5B that formed 58 QTL were detected. Among them, six stable QTL were identified on chromosomes 2AL, 2DS, 3BL, 3DL, 7AL, and 7BS. Meanwhile, three of the ten QTL identified in the DH population, QLuc.5A.1 and QLuc.5A.2 on chromosome 5AL and QLuc.6A.2 on 6AS, were stable and independently explained 5.58-10.86% of the phenotypic variation. The QLuc.6A.2 region colocalized with two MTAs identified by GWAS. Moreover, 71 carotenoid metabolism-related candidate genes were identified, and the allelic effects were analyzed in the MCC panel based on the 90K array. Results revealed that the genes CYP97A3 (Chr. 6B) and CCD1 (Chr. 5A) were significantly associated with LUC. Additionally, the gene PSY3 (QLuc.5A.1) and several candidate genes involved in the methylerythritol 4-phosphate (MEP) pathways colocalized with stable QTL regions. The present study provides potential targets for future functional gene exploration and molecular breeding in common wheat.

Whole-genome analysis of hard winter wheat germplasm identifies genomic regions associated with spike and kernel traits

Genetic dissection of yield component traits including spike and kernel characteristics is essential for the continuous improvement in wheat yield. Genome-wide association studies (GWAS) have been frequently used to identify genetic determinants for spike and kernel-related traits in wheat, though none have been employed in hard winter wheat (HWW) which represents a major class in US wheat acreage. Further, most of these studies relied on assembled diversity panels instead of adapted breeding lines, limiting the transferability of results to practical wheat breeding. Here we assembled a population of advanced/elite breeding lines and well-adapted cultivars and evaluated over four environments for phenotypic analysis of spike and kernel traits. GWAS identified 17 significant multi-environment marker-trait associations (MTAs) for various traits, representing 12 putative quantitative trait loci (QTLs), with five QTLs affecting multiple traits. Four of these QTLs mapped on three chromosomes 1A, 5B, and 7A for spike length, number of spikelets per spike (NSPS), and kernel length are likely novel. Further, a highly significant QTL was detected on chromosome 7AS that has not been previously associated with NSPS and putative candidate genes were identified in this region. The allelic frequencies of important quantitative trait nucleotides (QTNs) were deduced in a larger set of 1,124 accessions which revealed the importance of identified MTAs in the US HWW breeding programs. The results from this study could be directly used by the breeders to select the lines with favorable alleles for making crosses, and reported markers will facilitate marker-assisted selection of stable QTLs for yield components in wheat breeding.

New insight into the sucrose biosynthesis pathway from genome-wide identification, gene expression analysis, and subcellular localization in hexaploid wheat (Triticum aestivum L.)

Sucrose, the main synthesized product and transported form of photoassimilates, moves from leaves to support plant growth and storage substance biosynthesis occurring in the heterotrophic sink organs. Enhancing sucrose biosynthesis efficiency is a top priority for crop yield breeding programs. However, the molecular mechanism of sucrose biosynthesis is still elusive, especially in wheat. We performed transcriptome sequencing, subcellular localization, and bioinformatics analysis to identify key isoforms and metabolic branches associated with sucrose biosynthesis in wheat. Our analysis identified 45 genes from 13 families that exhibited high expression in young leaves with an evident diurnal change. The carbon flux from photoassimilates to sucrose was divided into two pathways. In the cytoplasm, assimilates initiating at phosphotrioses (TPs) exported by TaTPT1 from chloroplasts flowed along the TP-Sucrose branch formed by TaALD6, TaFBP5, TacPGI, TacPGM, TaUGP1, TaSPS5, and TaSPP1. Intermediates either from the Calvin cycle or TP-Sucrose branch were converted into ADPGlc to synthesize the simple starch, which was transiently degraded by a series of enzymes, including TaBAM4 and TaSEX4 in the chloroplast. Similar to the starch-biosynthesis branch in endosperms, the TP-Sucrose branch is the most prominent in leaves because each reaction can be catalyzed by at least one highly expressed isoform with expected cytosolic localization. The key isoforms and major branches highlighted in the wheat sucrose biosynthesis pathway expand our molecular understanding of crop sucrose biosynthesis and provide clues to increase wheat yield by enhancing the sucrose synthesis efficiency of leaves.

A catalogue of resistance gene homologs and a chromosome-scale reference sequence support resistance gene mapping in winter wheat

A resistance gene atlas is an integral component of the breeder's arsenal in the fight against evolving pathogens. Thanks to high-throughput sequencing, catalogues of resistance genes can be assembled even in crop species with large and polyploid genomes. Here, we report on capture sequencing and assembly of resistance gene homologs in a diversity panel of 907 winter wheat genotypes comprising ex situ genebank accessions and current elite cultivars. In addition, we use accurate long-read sequencing and chromosome conformation capture sequencing to construct a chromosome-scale genome sequence assembly of cv. Attraktion, an elite variety representative of European winter wheat. We illustrate the value of our resource for breeders and geneticists by (i) comparing the resistance gene complements in plant genetic resources and elite varieties and (ii) conducting genome-wide associations scans (GWAS) for the fungal diseases yellow rust and leaf rust using reference-based and reference-free GWAS approaches. The gene content under GWAS peaks was scrutinized in the assembly of cv. Attraktion.

CRISPR-Cas9 Gene Editing of the Sal1 Gene Family in Wheat

The highly conserved Sal1 encodes a bifunctional enzyme with inositol polyphosphate-1-phosphatase and 3′ (2′), 5′-bisphosphate nucleotidase activity and has been shown to alter abiotic stress tolerance in plants when disrupted. Precise gene editing techniques were used to generate Sal1 mutants in hexaploid bread wheat. The CRISPR (Clustered Regulatory Interspaced Short Palindromic Repeats) Cas9 system with three guide RNAs (gRNAs) was used to inactivate six Sal1 homologous genes within the Bobwhite wheat genome. The resulting mutant wheat plants with all their Sal1 genes disabled had slimmer stems, had a modest reduction in biomass and senesced more slowly in water limiting conditions, but did not exhibit improved yield under drought conditions. Our results show that multiplexed gRNAs enabled effective targeted gene editing of the Sal1 gene family in hexaploid wheat. These Sal1 mutant wheat plants will be a resource for further research studying the function of this gene family in wheat.

Genome-wide identification and expression analysis of the TaRRA gene family in wheat (Triticum aestivum L.)

Cytokinin is an important endogenous hormone in plants performing a wide spectrum of biological roles. The type-A response regulators (RRAs) are primary cytokinin response genes, which are important components of the cytokinin signaling pathway and are involved in the regulation of plant growth and development. By analysis of the whole genome sequence of wheat, we identified 20 genes encoding RRAs which were clustered into eight homologous groups. The gene structure, conserved motifs, chromosomal location, and cis-acting regulatory elements of the TaRRAs were analyzed. Quantitative real-time polymerase chain reaction (qRT-PCR) results showed that the expression levels of most of the TaRRAs increased rapidly on exogenous cytokinin application. Moreover, the TaRRA family members displayed different expression profiles under the stress treatments of drought, salt, cold, and heat. This study provides valuable insights into the RRA gene family in wheat and promotes the potential application of these genes in wheat genetic improvement.

Variations of wheat (Triticum aestivum L.) chromosomes caused by the 5A chromosomes with complex cytological structure

To study the effects of structural alterations of chromosomes caused by tandem repeats on the meiotic recombination, the wheat (Triticum aestivum L.) 5A chromosomes with different structure from ten wheat cultivars were used to investigate their meiotic recombination using non-denaturing fluorescence in situ hybridization (ND-FISH) technology. Fifteen cross combinations were carried out and they were divided into seven F₁ categories. The structural difference between the intercalary regions of the long arms of the two 5A chromosomes (5AL) in the F₁ categories III, VI, and VII was greater than that in the categories I and II, subsequently, the recombination frequencies in the distal regions of the 5AL arm in the progenies from the three categories were significantly lower than that from the categories I and II. For the two 5A chromosomes in the F₁ categories VI and VII, the structural differences in the distal regions of both of the two arms were greater than that in the categories IV and V. So, the recombination frequencies in the intercalary region of the 5AL arm in the progeny from the categories IV and V were higher than that in the progeny from the categories VI and VII. The breakage of 5A chromosome together with the 5A translocations and the breakage of some other chromosomes were observed in the progeny from the F₁ categories V, VI, and VII. These chromosomal variations were not observed in the progenies from the other four F₁ categories. In conclusion, the smaller structural difference between the 5A chromosomes in distal regions of the two arms resulted in a higher recombination frequency in interstitial region and vice versa. The 5A chromosome with complex cytological structure can be used to induce genetic variations of wheat genome.

A major quantitative trait locus for wheat total root length associated with precipitation distribution

Optimizing root system architecture (RSA) allows crops to better capture water and nutrients and adapt to harsh environment. Parental reproductive environment (PRE) has been reported to significantly affect growth and development throughout the life cycle of the next generation. In this study, 10 RSA-related traits were evaluated in seedling stage from five independent hydroponic tests using seeds harvested from five different PREs. Based on the Wheat55K SNP array-based genetic map, quantitative trait loci (QTL) for these traits were detected in a recombinant inbred line population. Twenty-eight putative QTL for RSA-related traits were detected, covering thirteen chromosomal regions. A major QTL, QTrl.sicau-2SY-4D for total root length (TRL), which was likely independent of PREs, explained 15.81-38.48% of phenotypic variations and was located at 14.96-19.59 Mb on chromosome arm 4DS. Interestingly, it showed pleiotropic effects on TRL, root area, root volume, root forks, root dry weight, and shoot dry weight. The functional marker KASP-Rht-D1 for Rht-D1 was used to genotype 2SY population and remapping QTL for TRL showed that QTrl.sicau-2SY-4D was not linked to Rht-D1. The kompetitive allele-specific PCR (KASP) marker, KASP-AX-110527441 linked to this major QTL, was developed and used to successfully validate its effect in three different genetic populations. Further analysis suggested that the positive allele at QTrl.sicau-2SY-4D was mainly utilized in wheat breeding of northwest China where precipitation was significantly lower, indicating that wheat requires longer TRL to capture water and nutrients in arid or semi-arid regions due to deficient precipitation. Additionally, four genes (TraesCS4D03G0059800, TraesCS4D03G0057800, TraesCS4D03G0064000, and TraesCS4D03G0064400) possibly related to root development were predicted in physical interval of QTrl.sicau-2SY-4D. Taken together, these results enrich our understanding on the genetic basis of RSA and provide a potentially valuable TRL QTL for wheat breeding.

Quick mapping and characterization of a co-located kernel length and thousand-kernel weight-related QTL in wheat

A co-located KL and TKW-related QTL with no negative effect on PH and AD was rapidly identified using BSA and wheat 660 K SNP array. Its effect was validated in a panel of 218 wheat accessions. Kernel length (KL) and thousand-kernel weight (TKW) of wheat (Triticum aestivum L.) contribute significantly to kernel yield. In the present study, a recombinant inbred line (RIL) population derived from the cross between the wheat line S849-8 with larger kernels and more spikelets per spike and the line SY95-71 was developed. Further, of both the bulked segregant analysis (BSA) and the wheat 660 K single nucleotide polymorphism (SNP) array were used to rapidly identify genomic regions for kernel-related traits from this RIL population. Kompetitive Allele Specific PCR markers were further developed in the SNP-enriched region on the 2D chromosome to construct a genetic map. Both QKL.sicau-SSY-2D for KL and QTKW.sicau-SSY-2D for TKW were identified at multiple environments on chromosome arm 2DL. These two QTLs explained 9.68-23.02% and 6.73-18.32% of the phenotypic variation, respectively. The effects of this co-located QTL were successfully verified in a natural population consisting of 218 Sichuan wheat accessions. Interestingly, the major QTL was significantly and positively correlated with spike length, but did not negatively affect spikelet number per spike (SNS), plant height, or anthesis date. These results indicated that it is possible to synchronously improve kernel weight and SNS by using this QTL. Additionally, several genes associated with kernel development and filling rate were predicted and sequenced in the QTL-containing physical intervals of reference genomes of 'Chinese spring' and Aegilops tauschii. Collectively, these results provide a QTL with great breeding potential and its linked markers which should be helpful for fine mapping and molecular breeding.

A genome-wide association study unveils key chromosome regions involved in determining sodium accumulation in wheat under conditions of low potassium supply

Improving nutrient use efficiency is an important objective in modern breeding programs. In this work, we examined potassium utilization efficiency (KUtE) and traits potentially related to it in a formerly genotyped, geographically diverse population of bread wheat (Triticum aestivum) under low potassium supply conditions. Our results unveil the existence of a large variation within the population for the traits examined. A genome-wide association study, based on a single-locus model, identified 15 markers associated with some of those traits. No marker-trait association was found using that tool for KUtE, but the use of a multi-locus approach suggested that additional marker-trait associations may be present, including whole-plant KUtE. Besides, the existence of a significant correlation between KUtE and sodium accumulation in shoots suggests the possibility of pyramiding traits associated with sodium homeostasis to improve this efficiency. In this regard, two discrete regions mapped on the long arm of chromosome 1B (1BLA and 1BLB) were associated with variation in sodium accumulation as detected with the single and multi-locus models used. Further exploration of the potential function of the genes placed in these regions, and their expression patterns, suggested likely candidates for this trait. Among the candidates placed in 1BLA region, we found TraesCS1B02G370500, TraesCS1B02G370600, and TraesCS1B02G370900, coding for putative Calcineurin B like proteins. Region 1BLB contain TraesCS1B02G388900 coding for a kinase and other genes including TraesCS1B02G389700, TraesCS1B02G389800 and TraesCS1B02G389900 coding for Ethylene-responsive transcription factors. The information here provided can be useful in breeding programs aimed to manipulate sodium accumulation through marker-assisted selection.

Haplotype variants of Sr46 in Aegilops tauschii, the diploid D genome progenitor of wheat

Stem rust resistance genes, SrRL5271 and Sr672.1 as well as SrCPI110651, from Aegilops tauschii, the diploid D genome progenitor of wheat, are sequence variants of Sr46 differing by 1-2 nucleotides leading to non-synonymous amino acid substitutions. The Aegilops tauschii (wheat D-genome progenitor) accessions RL 5271 and CPI110672 were identified as resistant to multiple races (including the Ug99) of the wheat stem rust pathogen Puccinia graminis f. sp. tritici (Pgt). This study was conducted to identify the stem rust resistance (Sr) gene(s) in both accessions. Genetic analysis of the resistance in RL 5271 identified a single dominant allele (SrRL5271) controlling resistance, whereas resistance segregated at two loci (SR672.1 and SR672.2) for a cross of CPI110672. Bulked segregant analysis placed SrRL5271 and Sr672.1 in a region on chromosome arm 2DS that encodes Sr46. Molecular marker screening, mapping and genomic sequence analysis demonstrated SrRL5271 and Sr672.1 are alleles of Sr46. The amino acid sequence of SrRL5271 and Sr672.1 is identical but differs from Sr46 (hereafter referred to as Sr46_h1 by following the gene nomenclature in wheat) by a single amino acid (N763K) and is thus designated Sr46_h2. Screening of a panel of Ae. tauschii accessions identified an additional allelic variant that differed from Sr46_h2 by a different amino acid (A648V) and was designated Sr46_h3. By contrast, the protein encoded by the susceptible allele of Ae. tauschii accession AL8/78 differed from these resistance proteins by 54 amino acid substitutions (94% nucleotide sequence gene identity). Cloning and complementation tests of the three resistance haplotypes confirmed their resistance to Pgt race 98-1,2,3,5,6 and partial resistance to Pgt race TTRTF in bread wheat. The three Sr46 haplotypes, with no virulent races detected yet, represent a valuable source for improving stem resistance in wheat.

Identification of QTLs for wheat heading time across multiple-environments

The genetic response to changing climatic factors selects consistent across the tested environments and location-specific thermo-sensitive and photoperiod susceptible alleles in lower and higher altitudes, respectively, for starting flowering in winter wheat. Wheat breeders select heading date to match the most favorable conditions for their target environments and this is favored by the extensive genetic variation for this trait that has the potential to be further explored. In this study, we used a germplasm with broad geographic distribution and tested it in multi-location field trials across Germany over three years. The genotypic response to the variation in the climatic parameters depending on location and year uncovered the effect of photoperiod and spring temperatures in accelerating heading date in higher and lower latitudes, respectively. Spring temperature dominates other factors in inducing heading, whereas the higher amount of solar radiation delays it. A genome-wide scan of marker-trait associations with heading date detected two QTL: an adapted allele at locus TaHd102 on chromosome 5A that has a consistent effect on HD in German cultivars in multiple environments and a non-adapted allele at locus TaHd044 on chromosome 3A that accelerates flowering by 5.6 days. TaHd102 and TaHd044 explain 13.8% and 33% of the genetic variance, respectively. The interplay of the climatic variables led to the detection of environment specific association responding to temperature in lower latitudes and photoperiod in higher ones. Another locus TaHd098 on chromosome 5A showed epistatic interactions with 15 known regulators of flowering time when non-adapted cultivars from outside Germany were included in the analysis.

The brassinosteroid biosynthesis gene TaD11-2A controls grain size and its elite haplotype improves wheat grain yields

TaD11-2A affects grain size and root length and its natural variations are associated with significant differences in yield-related traits in wheat. Brassinosteroids (BRs) control many important agronomic traits and therefore the manipulation of BR components could improve crop productivity and performance. However, the potential effects of BR-related genes on yield-related traits and stress tolerance in wheat (Triticum aestivum L.) remain poorly understood. Here, we identified TaD11 genes in wheat (rice D11 orthologs) that encoded enzymes involved in BR biosynthesis. TaD11 genes were highly expressed in roots (Zadoks scale: Z11) and grains (Z75), while expression was significantly suppressed by exogenous BR (24-epiBL). Ectopic expression of TaD11-2A rescued the abnormal panicle structure and plant height (PH) of the clustered primary branch 1 (cpb1) mutant, and also increased endogenous BR levels, resulting in improved grain yields and grain quality in rice. The tad11-2a-1 mutant displayed dwarfism, smaller grains, sensitivity to 24-epiBL, and reduced endogenous BR contents. Natural variations in TaD11-2A were associated with significant differences in yield-related traits, including PH, grain width, 1000-grain weight, and grain yield per plant, and its favorable haplotype, TaD11-2A-HapI was subjected to positive selection during wheat breeding. Additionally, TaD11-2A influenced root length and salt tolerance in rice and wheat at seedling stages. These results indicated the important role of BR TaD11 biosynthetic genes in controlling grain size and root length, and also highlighted their potential in the molecular biological analysis of wheat.