Tagging large CNV blocks in wheat boosts digitalization of germplasm resources by ultra-low-coverage sequencing.

BACKGROUND: The massive structural variations and frequent introgression highly contribute to the genetic diversity of wheat, while the huge and complex genome of polyploid wheat hinders efficient genotyping of abundant varieties towards accurate identification, management, and exploitation of germplasm resources. RESULTS: We develop a novel workflow that identifies 1240 high-quality large copy number variation blocks (CNVb) in wheat at the pan-genome level, demonstrating that CNVb can serve as an ideal DNA fingerprinting marker for discriminating massive varieties, with the accuracy validated by PCR assay. We then construct a digitalized genotyping CNVb map across 1599 global wheat accessions. Key CNVb markers are linked with trait-associated introgressions, such as the 1RS·1BL translocation and 2NvS translocation, and the beneficial alleles, such as the end-use quality allele Glu-D1d (Dx5 + Dy10) and the semi-dwarf r-e-z allele. Furthermore, we demonstrate that these tagged CNVb markers promote a stable and cost-effective strategy for evaluating wheat germplasm resources with ultra-low-coverage sequencing data, competing with SNP array for applications such as evaluating new varieties, efficient management of collections in gene banks, and describing wheat germplasm resources in a digitalized manner. We also develop a user-friendly interactive platform, WheatCNVb ( http://wheat.cau.edu.cn/WheatCNVb/ ), for exploring the CNVb profiles over ever-increasing wheat accessions, and also propose a QR-code-like representation of individual digital CNVb fingerprint. This platform also allows uploading new CNVb profiles for comparison with stored varieties. CONCLUSIONS: The CNVb-based approach provides a low-cost and high-throughput genotyping strategy for enabling digitalized wheat germplasm management and modern breeding with precise and practical decision-making.

Ploidy variation induces butterfly effect on chromatin topology in wheat.

Polyploidy is a prominent driver of plant diversification, accompanied with dramatic chromosomal rearrangement and epigenetic changes that affect gene expression. How chromatin interactions within and between subgenomes adapt to ploidy transition remains poorly understood. We generate open chromatin interaction maps for natural hexaploid wheat (AABBDD), extracted tetraploid wheat (AABB), diploid wheat progenitor Aegilops tauschii (DD) and resynthesized hexaploid wheat (RHW, AABBDD). Thousands of intra- and interchromosomal loops are de novo established or disappeared in AB subgenomes after separation of D subgenome, in which 37-95% of novel loops are lost again in RHW after merger of D genome. Interestingly, more than half of novel loops are formed by cascade reactions that are triggered by disruption of chromatin interaction between AB and D subgenomes. The interaction repressed genes in RHW relative to DD are expression suppressed, resulting in more balanced expression of the three homoeologs in RHW. The interaction levels of cascade anchors are decreased step-by-step. Leading single nucleotide polymorphisms of yield- and plant architecture-related quantitative trait locus are significantly enriched in cascade anchors. The expression of 116 genes interacted with these anchors are significantly correlated with the corresponding traits. Our findings reveal trans-regulation of intrachromosomal loops by interchromosomal interactions during genome merger and separation in polyploid species.

Molecular characterization of QTL for grain zinc and iron concentrations in wheat landrace Chinese Spring.

KEY MESSAGE: Three stable QTL for grain zinc concentration were identified in wheat landrace Chinese Spring. Favorable alleles were more frequent in landraces than in modern wheat cultivars. Wheat is a major source of dietary energy for the growing world population. Developing cultivars with enriched zinc and iron can potentially alleviate human micronutrient deficiency. In this study, a recombinant inbred line (RIL) population with 245 lines derived from cross Zhou 8425B/Chinese Spring was used to detect quantitative trait loci (QTL) for grain zinc concentration (GZnC) and grain iron concentration (GFeC) across four environments. Three stable QTL for GZnC with all favorable alleles from Chinese Spring were identified on chromosomes 3BL, 5AL, and 5BL. These QTL explaining maxima of 8.7%, 5.8%, and 7.1% of phenotypic variances were validated in 125 resequenced wheat accessions encompassing both landraces and modern cultivars using six kompetitive allele specific PCR (KASP) assays. The frequencies of favorable alleles for QGZnCzc.caas-3BL, QGZnCzc.caas-5AL and QGZnCzc.caas-5BL were higher in landraces (90.4%, 68.0%, and 100.0%, respectively) compared to modern cultivars (45.9%, 35.4%, and 40.9%), suggesting they were not selected in breeding programs. Candidate gene association studies on GZnC in the cultivar panel further delimited the QTL into 8.5 Mb, 4.1 Mb, and 47.8 Mb regions containing 46, 4, and 199 candidate genes, respectively. The 5BL QTL located in a region where recombination was suppressed. Two stable and three less stable QTL for GFeC with favorable alleles also from Chinese Spring were identified on chromosomes 4BS (Rht-B1a), 4DS (Rht-D1a), 1DS, 3AS, and 6DS. This study sheds light on the genetic basis of GZnC and GFeC in Chinese Spring and provides useful molecular markers for wheat biofortification.

Fine-mapping and validation of the major quantitative trait locus QFlANG-4B for flag leaf angle in wheat.

KEY MESSAGE: This study precisely mapped and validated a quantitative trait locus (QTL) located on chromosome 4B for flag leaf angle in wheat. Flag leaf angle (FLANG) is closely related to crop architecture and yield. We previously identified the quantitative trait locus (QTL) QFLANG-4B for FLANG on chromosome 4B, located within a 14-cM interval flanked by the markers Xbarc20 and Xzyh357, using a mapping population of recombinant inbred lines (RILs) derived from a cross between Nongda3331 (ND3331) and Zang1817. In this study, we fine-mapped QFLANG-4B and validated its associated genetic effect. We developed a BC3F3 population using ND3331 as the recurrent parent through marker-assisted selection, as well as near-isogenic lines (NILs) by selfing BC3F3 plants carrying different heterozygous segments for the QFLANG-4B region. We obtained eight recombinant types for QFLANG-4B, narrowing its location down to a 5.3-Mb region. This region contained 76 predicted genes, 7 of which we considered to be likely candidate genes for QFLANG-4B. Marker and phenotypic analyses of individual plants from the secondary mapping populations and their progeny revealed that the FLANG of the ND3331 allele is significantly higher than that of the Zang1817 allele in multiple environments. These results not only provide a basis for the map-based cloning of QFLANG-4B, but also indicate that QFLANG-4B has great potential for marker-assisted selection in wheat breeding programs designed to improve plant architecture and yield.

A k-mer-based pangenome approach for cataloging seed-storage-protein genes in wheat to facilitate genotype-to-phenotype prediction and improvement of end-use quality.

Wheat is a staple food for more than 35% of the world's population, with wheat flour used to make hundreds of baked goods. Superior end-use quality is a major breeding target; however, improving it is especially time-consuming and expensive. Furthermore, genes encoding seed-storage proteins (SSPs) form multi-gene families and are repetitive, with gaps commonplace in several genome assemblies. To overcome these barriers and efficiently identify superior wheat SSP alleles, we developed "PanSK" (Pan-SSP k-mer) for genotype-to-phenotype prediction based on an SSP-based pangenome resource. PanSK uses 29-mer sequences that represent each SSP gene at the pangenomic level to reveal untapped diversity across landraces and modern cultivars. Genome-wide association studies with k-mers identified 23 SSP genes associated with end-use quality that represent novel targets for improvement. We evaluated the effect of rye secalin genes on end-use quality and found that removal of ω-secalins from 1BL/1RS wheat translocation lines is associated with enhanced end-use quality. Finally, using machine-learning-based prediction inspired by PanSK, we predicted the quality phenotypes with high accuracy from genotypes alone. This study provides an effective approach for genome design based on SSP genes, enabling the breeding of wheat varieties with superior processing capabilities and improved end-use quality.

Positional cloning and characterization reveal the role of TaSRN-3D and TaBSR1 in the regulation of seminal root number in wheat.

Seminal roots play a critical role in water and nutrient absorption, particularly in the early developmental stages of wheat. However, the genes responsible for controlling SRN in wheat remain largely unknown. Genetic mapping and functional analyses identified a candidate gene (TraesCS3D01G137200, TaSRN-3D) encoding a Ser/Thr kinase glycogen synthase kinase 3 (STKc_GSK3) that regulated SRN in wheat. Additionally, experiments involving hormone treatment, nitrate absorption and protein interaction were conducted to explore the regulatory mechanism of TaSRN-3D. Results showed that the TaSRN-3D4332 allele inhibited seminal roots initiation and development, while loss-of-function mutants showed significantly higher seminal root number (SRN). Exogenous application of epi-brassinolide could increase the SRN in a HS2-allelic background. Furthermore, chlorate sensitivity and 15N uptake assays revealed that a higher number of seminal roots promoted nitrate accumulation. TaBSR1 (BIN2-related SRN Regulator 1, orthologous to OsGRF4/GL2 in rice) acts as an interactor of TaSRN-3D and promotes TaBSR1 degradation to reduce SRN. This study provides valuable insights into understanding the genetic basis and regulatory network of SRN in wheat, highlighting their roles as potential targets for root-based improvement in wheat breeding.

The potassium transporter TaNHX2 interacts with TaGAD1 to promote drought tolerance via modulating stomatal aperture in wheat.

Drought is a major global challenge in agriculture that decreases crop production. γ-Aminobutyric acid (GABA) interfaces with drought stress in plants; however, a mechanistic understanding of the interaction between GABA accumulation and drought response remains to be established. Here we showed the potassium/proton exchanger TaNHX2 functions as a positive regulator in drought resistance in wheat by mediating cross-talk between the stomatal aperture and GABA accumulation. TaNHX2 interacted with glutamate decarboxylase TaGAD1, a key enzyme that synthesizes GABA from glutamate. Furthermore, TaNHX2 targeted the C-terminal auto-inhibitory domain of TaGAD1, enhanced its activity, and promoted GABA accumulation under drought stress. Consistent with this, the tanhx2 and tagad1 mutants showed reduced drought tolerance, and transgenic wheat with enhanced TaNHX2 expression had a yield advantage under water deficit without growth penalty. These results shed light on the plant stomatal movement mechanism under drought stress and the TaNHX2-TaGAD1 module may be harnessed for amelioration of negative environmental effects in wheat as well as other crops.

Innovative computational tools provide new insights into the polyploid wheat genome.

Bread wheat (Triticum aestivum) is an important crop and serves as a significant source of protein and calories for humans, worldwide. Nevertheless, its large and allopolyploid genome poses constraints on genetic improvement. The complex reticulate evolutionary history and the intricacy of genomic resources make the deciphering of the functional genome considerably more challenging. Recently, we have developed a comprehensive list of versatile computational tools with the integration of statistical models for dissecting the polyploid wheat genome. Here, we summarize the methodological innovations and applications of these tools and databases. A series of step-by-step examples illustrates how these tools can be utilized for dissecting wheat germplasm resources and unveiling functional genes associated with important agronomic traits. Furthermore, we outline future perspectives on new advanced tools and databases, taking into consideration the unique features of bread wheat, to accelerate genomic-assisted wheat breeding.

Genomic insights into the origin and evolution of spelt (Triticum spelta L.) as a valuable gene pool for modern wheat breeding.

Spelt (Triticum aestivum ssp. spelta) is an important wheat subspecies mainly cultivated in Europe before the 20th century that has contributed to modern wheat breeding as a valuable genetic resource. However, relatively little is known about the origins and maintenance of spelt populations. Here, using resequencing data from 416 worldwide wheat accessions, including representative spelt wheat, we demonstrate that European spelt emerged when primitive hexaploid wheat spread to the west and hybridized with pre-settled domesticated emmer, the putative maternal donor. Genomic introgression regions from domesticated emmer confer spelt's primitive morphological characters used for species taxonomy, such as tenacious glumes and later flowering. We propose a haplotype-based "spelt index" to identify spelt-type wheat varieties and to quantify utilization of the spelt gene pool in modern wheat cultivars. This study reveals the genetic basis for the establishment of the spelt wheat subspecies in a specific ecological niche and the vital role of the spelt gene pool as a unique germplasm resource in modern wheat breeding.

TaNAM-6A is essential for nitrogen remobilisation and regulates grain protein content in wheat (Triticum aestivum L.).

Grain protein content (GPC) is a crucial quality trait in bread wheat, which is influenced by the key transcription factor TaNAM. However, the regulatory mechanisms of TaNAM have remained largely elusive. In this study, a new role of TaNAM was unveiled in regulating nitrogen remobilisation which impacts GPC. The TaNAM knockout mutants generated by clustered regularly interspaced short palindromic repeats/Cas9 exhibited significantly delayed senescence and lower GPC, while overexpression of TaNAM-6A resulted in premature senility and much higher GPC. Further analysis revealed that TaNAM directly activates the genes TaNRT1.1 and TaNPF5.5s, which are involved in nitrogen remobilisation. This activity aids in the transfer of nitrogen from leaves to grains for protein synthesis. In addition, an elite allele of TaNAM-6A, associated with high GPC, was identified as a candidate gene for breeding high-quality wheat. Overall, our work not only elucidates the potential mechanism of TaNAM-6A affecting bread wheat GPC, but also highlights the significance of nitrogen remobilisation from senescent leaves to grains for protein accumulation. Moreover, our research provides a new target and approach for improving the quality traits of wheat, particularly the GPC.

Natural variation of STKc_GSK3 kinase TaSG-D1 contributes to heat stress tolerance in Indian dwarf wheat.

Heat stress threatens global wheat (Triticum aestivum) production, causing dramatic yield losses worldwide. Identifying heat tolerance genes and comprehending molecular mechanisms are essential. Here, we identify a heat tolerance gene, TaSG-D1E286K, in Indian dwarf wheat (Triticum sphaerococcum), which encodes an STKc_GSK3 kinase. TaSG-D1E286K improves heat tolerance compared to TaSG-D1 by enhancing phosphorylation and stability of downstream target TaPIF4 under heat stress condition. Additionally, we reveal evolutionary footprints of TaPIF4 during wheat selective breeding in China, that is, InDels predominantly occur in the TaPIF4 promoter of Chinese modern wheat cultivars and result in decreased expression level of TaPIF4 in response to heat stress. These sequence variations with negative effect on heat tolerance are mainly introduced from European germplasm. Our study provides insight into heat stress response mechanisms and proposes a potential strategy to improve wheat heat tolerance in future.

Methyltransferase TaSAMT1 mediates wheat freezing tolerance by integrating brassinosteroid and salicylic acid signaling.

Cold injury is a major environmental stress affecting the growth and yield of crops. Brassinosteroids (BRs) and salicylic acid (SA) play important roles in plant cold tolerance. However, whether or how BR signaling interacts with the SA signaling pathway in response to cold stress is still unknown. Here, we identified an SA methyltransferase, TaSAMT1 that converts SA to methyl SA (MeSA) and confers freezing tolerance in wheat (Triticum aestivum). TaSAMT1 overexpression greatly enhanced wheat freezing tolerance, with plants accumulating more MeSA and less SA, whereas Tasamt1 knockout lines were sensitive to freezing stress and accumulated less MeSA and more SA. Spraying plants with MeSA conferred freezing tolerance to Tasamt1 mutants, but SA did not. We revealed that BRASSINAZOLE-RESISTANT 1 (TaBZR1) directly binds to the TaSAMT1 promoter and induces its transcription. Moreover, TaBZR1 interacts with the histone acetyltransferase TaHAG1, which potentiates TaSAMT1 expression via increased histone acetylation and modulates the SA pathway during freezing stress. Additionally, overexpression of TaBZR1 or TaHAG1 altered TaSAMT1 expression and improved freezing tolerance. Our results demonstrate a key regulatory node that connects the BR and SA pathways in the plant cold stress response. The regulatory factors or genes identified could be effective targets for the genetic improvement of freezing tolerance in crops.

Reshaped DNA methylation cooperating with homoeolog-divergent expression promotes improved root traits in synthesized tetraploid wheat.

Polyploidization is a major event driving plant evolution and domestication. However, how reshaped epigenetic modifications coordinate gene transcription to generate phenotypic variations during wheat polyploidization is currently elusive. Here, we profiled transcriptomes and DNA methylomes of two diploid wheat accessions (SlSl and AA) and their synthetic allotetraploid wheat line (SlSlAA), which displayed elongated root hair and improved root capability for nitrate uptake and assimilation after tetraploidization. Globally decreased DNA methylation levels with a reduced difference between subgenomes were observed in the roots of SlSlAA. DNA methylation changes in first exon showed strong connections with altered transcription during tetraploidization. Homoeolog-specific transcription was associated with biased DNA methylation as shaped by homoeologous sequence variation. The hypomethylated promoters showed significantly enriched binding sites for MYB, which may affect gene transcription in response to root hair growth. Two master regulators in root hair elongation pathway, AlCPC and TuRSL4, exhibited upregulated transcription levels accompanied by hypomethylation in promoter, which may contribute to the elongated root hair. The upregulated nitrate transporter genes, including NPFs and NRTs, also are significantly associated with hypomethylation, indicating an epigenetic-incorporated regulation manner in improving nitrogen use efficiency. Collectively, these results provided new insights into epigenetic changes in response to crop polyploidization and underscored the importance of epigenetic regulation in improving crop traits.

Natural variations of HvSRN1 modulate the spike rachis node number in barley.

Grain number, one of the major determinants of yield in Triticeae crops, is largely determined by spikelet number and spike rachis node number (SRN). Here, we identified three quantitative trait loci (QTLs) for SRN using 145 recombinant inbred lines derived from a barley R90/1815D cross. qSRN1, the major-effect QTL, was mapped to chromosome 2H and explained up to 38.77% of SRN variation. Map-based cloning revealed that qSRN1 encodes the RAWUL domain-containing protein HvSRN1. Further analysis revealed that two key SNPs in the HvSRN1 promoter region (∼2 kb upstream of the transcription start site) affect the transcript level of HvSRN1 and contribute to variation in SRN. Similar to its orthologous proteins OsLAX2 and ZmBA2, HvSRN1 showed protein-protein interactions with HvLAX1, suggesting that the LAX2-LAX1 model for spike morphology regulation may be conserved in Poaceae crops. CRISPR-Cas9-induced HvSRN1 mutants showed reduced SRN but increased grain size and weight, demonstrating a trade-off effect. Our results shed light on the role of HvSRN1 variation in regulating the balance between grain number and weight in barley.

Type I MADS-box transcription factor TaMADS-GS regulates grain size by stabilizing cytokinin signalling during endosperm cellularization in wheat.

Grain size is one of the important traits in wheat breeding programs aimed at improving yield, and cytokinins, mainly involved in cell division, have a positive impact on grain size. Here, we identified a novel wheat gene TaMADS-GS encoding type I MADS-box transcription factor, which regulates the cytokinins signalling pathway during early stages of grain development to modulate grain size and weight in wheat. TaMADS-GS is exclusively expressed in grains at early stage of seed development and its knockout leads to delayed endosperm cellularization, smaller grain size and lower grain weight. TaMADS-GS protein interacts with the Polycomb Repressive Complex 2 (PRC2) and leads to repression of genes encoding cytokinin oxidase/dehydrogenases (CKXs) stimulating cytokinins inactivation by mediating accumulation of the histone H3 trimethylation at lysine 27 (H3K27me3). Through the screening of a large wheat germplasm collection, an elite allele of the TaMADS-GS exhibits higher ability to repress expression of genes inactivating cytokinins and a positive correlation with grain size and weight, thus representing a novel marker for breeding programs in wheat. Overall, these findings support the relevance of TaMADS-GS as a key regulator of wheat grain size and weight.

Mutation of a highly conserved amino acid in RPM1 causes leaf yellowing and premature senescence in wheat.

KEY MESSAGE: A point mutation of RPM1 triggers persistent immune response that induces leaf premature senescence in wheat, providing novel information of immune responses and leaf senescence. Leaf premature senescence in wheat (Triticum aestivum L.) is one of the most common factors affecting the plant's development and yield. In this study, we identified a novel wheat mutant, yellow leaf and premature senescence (ylp), which exhibits yellow leaves and premature senescence at the heading and flowering stages. Consistent with the yellow leaves phenotype, ylp had damaged and collapsed chloroplasts. Map-based cloning revealed that the phenotype of ylp was caused by a point mutation from Arg to His at amino acid 790 in a plasma membrane-localized protein resistance to Pseudomonas syringae pv. maculicola 1 (RPM1). The point mutation triggered excessive immune responses and the upregulation of senescence- and autophagy-associated genes. This work provided the information for understanding the molecular regulatory mechanism of leaf senescence, and the results would be important to analyze which mutations of RPM1 could enable plants to obtain immune activation without negative effects on plant growth.

Paternally imprinted LATE-FLOWERING2 transcription factor contributes to paternal-excess interploidy hybridization barriers in wheat.

Interploidy hybridization between hexaploid and tetraploid genotypes occurred repeatedly during genomic introgression events throughout wheat evolution, and is commonly employed in wheat breeding programs. Hexaploid wheat usually serves as maternal parent because the reciprocal cross generates progeny with severe defects and poor seed germination, but the underlying mechanism is poorly understood. Here, we performed detailed analysis of phenotypic variation in endosperm between two interploidy reciprocal crosses arising from tetraploid (Triticum durum, AABB) and hexaploid wheat (Triticum aestivum, AABBDD). In the paternal- versus the maternal-excess cross, the timing of endosperm cellularization was delayed and starch granule accumulation in the endosperm was repressed, causing reduced germination percentage. The expression profiles of genes involved in nutrient metabolism differed strongly between these endosperm types. Furthermore, expression patterns of parental alleles were dramatically disturbed in interploidy versus intraploidy crosses, leading to increased number of imprinted genes. The endosperm-specific TaLFL2 showed a paternally imprinted expression pattern in interploidy crosses partially due to allele-specific DNA methylation. Paternal TaLFL2 binds to and represses a nutrient accumulation regulator TaNAC019, leading to reduced storage protein and starch accumulation during endosperm development in paternal-excess cross, as confirmed by interploidy crosses between tetraploid wild-type and clustered regularly interspaced palindromic repeats (CRISPR) - CRISPR-associated protein 9 generated hexaploid mutants. These findings reveal a contribution of genomic imprinting to paternal-excess interploidy hybridization barriers during wheat evolution history and explains why experienced breeders preferentially exploit maternal-excess interploidy crosses in wheat breeding programs.

Phosphorylation of KAT-2B by WKS1/Yr36 redirects the lipid flux to jasmonates to enhance resistance against wheat stripe rust.

Wheat (Triticum aestivum) is one of the most essential human energy and protein sources. However, wheat production is threatened by devastating fungal diseases such as stripe rust, caused by Puccinia striiformis Westend. f. sp. tritici (Pst). Here, we reveal that the alternations in chloroplast lipid profiles and the accumulation of jasmonate (JA) in the necrosis region activate JA signaling and trigger the host defense. The collapse of chloroplasts in the necrosis region results in accumulations of polyunsaturated membrane lipids and the lipid-derived phytohormone JA in transgenic lines of Yr36 that encodes Wheat Kinase START 1 (WKS1), a high-temperature-dependent adult plant resistance protein. WKS1.1, a protein encoded by a full-length splicing variant of WKS1, phosphorylates and enhances the activity of keto-acyl thiolase (KAT-2B), a critical enzyme catalyzing the ß-oxidation reaction in JA biosynthesis. The premature stop mutant, kat-2b, accumulates less JA and shows defects in the host defense against Pst. Conversely, overexpression of KAT-2B results in a higher level of JA and limits the growth of Pst. Moreover, JA inhibits the growth and reduces pustule densities of Pst. This study illustrates the WKS1.1‒KAT-2B‒JA pathway for enhancing wheat defense against fungal pathogens to attenuate yield loss.