A Collinearity-Incorporating Homology Inference Strategy for Connecting Emerging Assemblies in the Triticeae Tribe as a Pilot Practice in the Plant Pangenomic Era

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.

Natural variant of Rht27, a dwarfing gene, enhances yield potential in wheat

KEY MESSAGE

Discovery of Rht27, a dwarf gene in wheat, showed potential in enhancing grain yield by reducing plant height. Plant height plays a crucial role in crop architecture and grain yield, and semi-dwarf Reduced Height (Rht) alleles contribute to lodging resistance and were important in "Green Revolution." However, the use of these alleles is associated with some negative side effects in some environments, such as reduced coleoptile length, low nitrogen use efficiency, and reduced yield. Therefore, novel dwarf gene resources are needed to pave an alternative route to overcome these side effects. In this study, a super-dwarf mutant rht27 was obtained by the mutagenesis of G1812 (Triticum urartu, the progenitor of the A sub-genome of common wheat). Genetic analysis revealed that the dwarf phenotype was regulated by a single recessive genetic factor. The candidate region for Rht27 was narrowed to a 1.55 Mb region on chromosome 3, within which we found two potential candidate genes that showed polymorphisms between the mutant and non-mutagenized G1812. Furthermore, the natural variants and elite haplotypes of the two candidates were investigated in a natural population of common wheat. The results showed that the natural variants affect grain yield components, and the dwarf haplotypes show the potential in improving agronomic traits and grain yield. Although the mutation in Rht27 results in severe dwarf phenotype in T. urartu, the natural variants in common wheat showed desirable phenotype, which suggests that Rht27 has the potential to improve wheat yield by utilizing its weak allelic mutation or fine-tuning its expression level.

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.

Identification of candidate genes for adult plant stripe rust resistance transferred from Aegilops ventricosa 2NvS into wheat via fine mapping and transcriptome analysis

KEY MESSAGE

An adult plant gene for resistance to stripe rust was narrowed down to the proximal one-third of the 2NvS segment translocated from Aegilops ventricosa to wheat chromosome arm 2AS, and based on the gene expression analysis, two candidate genes were identified showing a stronger response at the adult plant stage compared to the seedling stage. The 2NvS translocation from Aegilops ventricosa, known for its resistance to various diseases, has been pivotal in global wheat breeding for more than three decades. Here, we identified an adult plant resistance (APR) gene in the 2NvS segment in wheat line K13-868. Through fine mapping in a segregating near-isogenic line (NIL) derived population of 6389 plants, the candidate region for the APR gene was narrowed down to between 19.36 Mb and 33 Mb in the Jagger reference genome. Transcriptome analysis in NILs strongly suggested that this APR gene conferred resistance to stripe rust by triggering plant innate immune responses. Based on the gene expression analysis, two disease resistance-associated genes within the candidate region, TraesJAG2A03G00588940 and TraesJAG2A03G00590140, exhibited a stronger response to Puccinia striiformis f. sp. tritici (Pst) infection at the adult plant stage than at the seedling stage, indicating that they could be potential candidates for the resistance gene. Additionally, we developed a co-dominant InDel marker, InDel_31.05, for detecting this APR gene. Applying this marker showed that over one-half of the wheat varieties approved in 2021 and 2022 in Sichuan province, China, carry this gene. Agronomic trait evaluation of NILs indicated that the 2NvS segment effectively mitigated the negative effects of stripe rust on yield without affecting other important agronomic traits. This study provided valuable insights for cloning and breeding through the utilization of the APR gene present in the 2NvS segment.

Identification of the succinate-CoA ligase protein gene family reveals that TaSUCL1-1 positively regulate cadmium resistance in wheat

The Succinate-CoA ligase (SUCL1) gene family is involved in energy metabolism, phytohormone signaling, and plant growth, development, and tolerance to stress. This is the first study to analyze the SUCL1 gene family in wheat (Triticum aestivum). 17 TaSUCL1 genes were identified in the complete genome sequence and classified into five subfamilies based on related genes found in three other species. The 17 TaSUCL1 genes were unevenly distributed across 11 chromosomes, and the collinearity of these genes was further investigated. Through using real-time qPCR (RT-qPCR) analysis, we identified the expression patterns of the TaSUCL1 genes under various tissues and different heavy metal stress conditions. The functions of selected TaSUCL1-1 gene were investigated by RNA interference (RNAi). This study provided a comprehensive analysis of the TaSUCL1 gene family. Within the TaSUCL1 genes, the exon-intron structure and motif composition exhibited significant similarity among members of the same evolutionary branch. Homology analysis and phylogenetic comparison of the SUCL1 genes in different plants offered valuable insights for studying the evolutionary characteristics of the SUCL1 genes. The expression levels of the TaSUCL1 genes in different tissues and under various metal stress conditions reveal its important role in plant growth and development. Gene function analysis demonstrated that TaSUCL1-1 silenced wheat plants exhibited a decrease in the total cadmium (Cd) concentrations and gene expression levels compared to the wild type (WT). Additionally, TaSUCL1-1 belonging to class c physically interacts with the β-amylase protein TaBMY1 as verified by yeast two-hybridization. This research provides a useful resource for further study of the function and molecular genetic mechanism of the SUCL1 gene family members.

Epigenetic control on transcription of vernalization genes and whole-genome gene expression profile induced by vernalization in common wheat

Vernalization is necessary for winter wheat to flower. However, it is unclear whether vernalization is also required for spring wheat, which is frequently sown in fall, and what molecular mechanisms underlie the vernalization response in wheat varieties. In this study, we examined the molecular mechanisms that regulate vernalization response in winter and spring wheat varieties. For this purpose, we determined how major vernalization genes (VRN1, VRN2, and VRN3) respond to vernalization in these varieties and whether modifications to histones play a role in changes in gene expression. We also identified genes that are differentially regulated in response to vernalization in winter and spring wheat varieties. We found that in winter wheat, but not in spring wheat, VRN1 expression decreases when returned to warm temperature following vernalization. This finding may be associated with differences between spring and winter wheat in the levels of tri-methylation of lysine 27 on histone H3 (H3K27me3) and tri-methylation of lysine 4 on histone H3 (H3K4me3) at the VRN1 gene. Analysis of winter wheat transcriptomes before and after vernalization revealed that vernalization influences the expression of several genes, including those involved in leucine catabolism, cysteine biosynthesis, and flavonoid biosynthesis. These findings provide new candidates for further study on the mechanism of vernalization regulation in wheat.

Using UAV-Based Temporal Spectral Indices to Dissect Changes in the Stay-Green Trait in Wheat

Stay-green (SG) in wheat is a beneficial trait that increases yield and stress tolerance. However, conventional phenotyping techniques limited the understanding of its genetic basis. Spectral indices (SIs) as non-destructive tools to evaluate crop temporal senescence provide an alternative strategy. Here, we applied SIs to monitor the senescence dynamics of 565 diverse wheat accessions from anthesis to maturation stages over 2 field seasons. Four SIs (normalized difference vegetation index, green normalized difference vegetation index, normalized difference red edge index, and optimized soil-adjusted vegetation index) were normalized to develop relative stay-green scores (RSGS) as the SG indicators. An RSGS-based genome-wide association study identified 47 high-confidence quantitative trait loci (QTL) harboring 3,079 single-nucleotide polymorphisms associated with SG and 1,085 corresponding candidate genes. Among them, 15 QTL overlapped or were adjacent to known SG-related QTL/genes, while the remaining QTL were novel. Notably, a set of favorable haplotypes of SG-related candidate genes such as TraesCS2A03G1081100, TracesCS6B03G0356400, and TracesCS2B03G1299500 are increasing following the Green Revolution, further validating the feasibility of the pipeline. This study provided a valuable reference for further quantitative SG and genetic research in diverse wheat panels.

High-quality genome assembly enables prediction of allele-specific gene expression in hybrid poplar

Poplar (Populus) is a well-established model system for tree genomics and molecular breeding, and hybrid poplar is widely used in forest plantations. However, distinguishing its diploid homologous chromosomes is difficult, complicating advanced functional studies on specific alleles. In this study, we applied a trio-binning design and PacBio high-fidelity long-read sequencing to obtain haplotype-phased telomere-to-telomere genome assemblies for the 2 parents of the well-studied F1 hybrid "84K" (Populus alba × Populus tremula var. glandulosa). Almost all chromosomes, including the telomeres and centromeres, were completely assembled for each haplotype subgenome apart from 2 small gaps on one chromosome. By incorporating information from these haplotype assemblies and extensive RNA-seq data, we analyzed gene expression patterns between the 2 subgenomes and alleles. Transcription bias at the subgenome level was not uncovered, but extensive-expression differences were detected between alleles. We developed machine-learning (ML) models to predict allele-specific expression (ASE) with high accuracy and identified underlying genome features most highly influencing ASE. One of our models with 15 predictor variables achieved 77% accuracy on the training set and 74% accuracy on the testing set. ML models identified gene body CHG methylation, sequence divergence, and transposon occupancy both upstream and downstream of alleles as important factors for ASE. Our haplotype-phased genome assemblies and ML strategy highlight an avenue for functional studies in Populus and provide additional tools for studying ASE and heterosis in hybrids.

Dosage effect genes modulate grain development in synthesized Triticum durum-Haynaldia villosa allohexaploid

Polyploidization in plants often leads to increased cell size and grain size, which may be affected by the increased genome dosage and transcription abundance. The synthesized Triticum durum (AABB)-Haynaldia villosa (VV) amphiploid (AABBVV) has significantly increased grain size, especially grain length, than the tetraploid and diploid parents. To investigate how polyploidization affects grain development at the transcriptional level, we perform transcriptome analysis using the immature seeds of T. durum, H. villosa, and the amphiploid. The dosage effect genes are contributed more by differentially expressed genes from genome V of H. villosa. The dosage effect genes overrepresent grain development-related genes. Interestingly, the vernalization gene TaVRN1 is among the positive dosage effect genes in the T. durum‒H. villosa and T. turgidum‒Ae. tauschii amphiploids. The expression levels of TaVRN1 homologs are positively correlated with the grain size and weight. The TaVRN1-B1 or TaVRN1-D1 mutation shows delayed florescence, decreased cell size, grain size, and grain yield. These data indicate that dosage effect genes could be one of the important explanations for increased grain size by regulating grain development. The identification and functional validation of dosage effect genes may facilitate the finding of valuable genes for improving wheat yield.

The era of panomics-driven gene discovery in plants

Panomics is an approach to integrate multiple 'omics' datasets, generated using different individuals or natural variations. Considering their diverse phenotypic spectrum, the phenome is inherently associated with panomics-based science, which is further combined with genomics, transcriptomics, metabolomics, and other omics techniques, either independently or collectively. Panomics has been accelerated through recent technological advancements in the field of genomics that enable the detection of population-wide structural variations (SVs) and hence offer unprecedented insights into the genetic variations contributing to important agronomic traits. The present review provides the recent trends of panomics-driven gene discovery toward various traits related to plant development, stress tolerance, accumulation of specialized metabolites, and domestication/dedomestication. In addition, the success stories are highlighted in the broader context of enhancing crop productivity.

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.

Wheat powdery mildew resistance gene Pm13 encodes a mixed lineage kinase domain-like protein

Wheat powdery mildew is one of the most destructive diseases threatening global wheat production. The wild relatives of wheat constitute rich sources of diversity for powdery mildew resistance. Here, we report the map-based cloning of the powdery mildew resistance gene Pm13 from the wild wheat species Aegilops longissima. Pm13 encodes a mixed lineage kinase domain-like (MLKL) protein that contains an N-terminal-domain of MLKL (MLKL_NTD) domain in its N-terminus and a C-terminal serine/threonine kinase (STK) domain. The resistance function of Pm13 is validated by mutagenesis, gene silencing, transgenic assay, and allelic association analyses. The development of introgression lines with significantly reduced chromosome segments of Ae. longissima encompassing Pm13 enables widespread deployment of this gene into wheat cultivars. The cloning of Pm13 may provide valuable insights into the molecular mechanisms underlying Pm13-mediated powdery mildew resistance and highlight the important roles of kinase fusion proteins (KFPs) in wheat immunity.

Genetic dissection and identification of stripe rust resistance genes in the wheat cultivar Lanhangxuan 121, a cultivar selected from a space mutation population

Stripe rust is a devastating disease of wheat worldwide. Chinese wheat cultivar Lanhangxuan 121 (LHX121), selected from an advanced line L92-47 population that had been subjected to space mutation breeding displayed a consistently higher level of resistance to stipe rust than its parent in multiple field environments. The aim of this research was to establish the number and types of resistance genes in parental lines L92-47 and LHX121 using separate segregating populations. The first population developed from a cross between LHX121 and susceptible cultivar Xinong 822 comprised 278 F2:3 lines. The second validation population comprised 301 F2:3 lines from a cross between L92-47 and susceptible cultivar Xinong 979. Lines of two population were evaluated for stripe rust response at three sites during the 2018-2020 cropping season. Affymetrix 660 K SNP arrays were used to genotype the lines and parents. Inclusive composite interval mapping detected QTL QYrLHX.nwafu-2BS, QYrLHX.nwafu-3BS, and QYrLHX.nwafu-5BS for resistance in all three environments. Based on previous studies and pedigree information, QYrLHX.nwafu-2BS and QYrLHX.nwafu-3BS were likely to be Yr27 and Yr30 that are present in the L92-47 parent. QYrLHX.nwafu-5BS (YrL121) detected only in LHX121 was mapped to a 7.60 cM interval and explained 10.67-22.57% of the phenotypic variation. Compared to stripe rust resistance genes previously mapped to chromosome 5B, YrL121 might be a new adult plant resistance QTL. Furthermore, there were a number of variations signals using 35 K SNP array and differentially expressed genes using RNA-seq between L92-47 and LHX121 in the YrL121 region, indicating that they probably impair the presence and/or function of YrL121.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01461-0.

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.

Two haplotype-resolved genome assemblies for AAB allotriploid bananas provide insights into banana subgenome asymmetric evolution and Fusarium wilt control

Bananas (Musa spp.) are one of the world's most important fruit crops and play a vital role in food security for many developing countries. Most banana cultivars are triploids derived from inter- and intraspecific hybridizations between the wild diploid ancestor species Musa acuminate (AA) and M. balbisiana (BB). We report two haplotype-resolved genome assemblies of the representative AAB-cultivated types, Plantain and Silk, and precisely characterize ancestral contributions by examining ancestry mosaics across the genome. Widespread asymmetric evolution is observed in their subgenomes, which can be linked to frequent homologous exchange events. We reveal the genetic makeup of triploid banana cultivars and verify that subgenome B is a rich source of disease resistance genes. Only 58.5% and 59.4% of Plantain and Silk genes, respectively, are present in all three haplotypes, with >50% of genes being differentially expressed alleles in different subgenomes. We observed that the number of upregulated genes in Plantain is significantly higher than that in Silk at one-week post-inoculation with Fusarium wilt tropical race 4 (Foc TR4), which confirms that Plantain can initiate defense responses faster than Silk. Additionally, we compared genomic and transcriptomic differences among the genes related to carotenoid synthesis and starch metabolism between Plantain and Silk. Our study provides resources for better understanding the genomic architecture of cultivated bananas and has important implications for Musa genetics and breeding.

Dissecting the genetic basis of Fusarium crown rot resistance in wheat by genome wide association study

KEY MESSAGE

A locus conferring Fusarium crown rot resistance was identified on chromosome arm 3DL through genome wide association study and further validated in two recombinant inbred lines populations. Fusarium crown rot (FCR) is a severe soil borne disease in many wheat growing regions of the world. In this study, we attempted to detect loci conferring FCR resistance through a new seedling inoculation assay. A total of 223 wheat accessions from different geography origins were used to assemble an association panel for GWAS analysis. Four genotypes including Heng 4332, Luwanmai, Pingan 998 and Yannong 24 showed stable resistance to FCR. A total of 54 SNPs associated with FCR resistance were identified. Among the 10 putative QTLs represented by these SNPs, seven QTLs on chromosome 2B, 3A, 3D, 4A, 7A and 7B were novel and were consistently detected in at least two of the three trials conducted. Qfcr.cau.3D-3, which was targeted by 38 SNPs clustered within a genomic region of approximately 5.57 Mb (609.12-614.69 Mb) on chromosome arm 3DL, was consistently detected in all the three trials. The effects of Qfcr.cau.3D-3 were further validated in two recombinant inbred line populations. The presence of this locus reduced FCR severity up to 21.55%. Interestingly, the collinear positions of sequences containing the four SNPs associated with two FCR loci (Qfcr.cau.3A and Qfcr.cau.3B) were within the regions of Qfcr.cau.3D-3, suggesting that genes underlying these three loci may be homologous. Our results provide useful information for improving FCR resistance in wheat.

The mechanisms behind the contrasting responses to waterlogging in black-grass (Alopecurus myosuroides) and wheat (Triticum aestivum)

Black-grass (Alopecurus myosuroides ) is one of the most problematic agricultural weeds of Western Europe, causing significant yield losses in winter wheat (Triticum aestivum ) and other crops through competition for space and resources. Previous studies link black-grass patches to water-retaining soils, yet its specific adaptations to these conditions remain unclear. We designed pot-based waterlogging experiments to compare 13 biotypes of black-grass and six cultivars of wheat. These showed that wheat roots induced aerenchyma when waterlogged whereas aerenchyma-like structures were constitutively present in black-grass. Aerial biomass of waterlogged wheat was smaller, whereas waterlogged black-grass was similar or larger. Variability in waterlogging responses within and between these species was correlated with transcriptomic and metabolomic changes in leaves of control or waterlogged plants. In wheat, transcripts associated with regulation and utilisation of phosphate compounds were upregulated and sugars and amino acids concentrations were increased. Black-grass biotypes showed limited molecular responses to waterlogging. Some black-grass amino acids were decreased and one transcript commonly upregulated was previously identified in screens for genes underpinning metabolism-based resistance to herbicides. Our findings provide insights into the different waterlogging tolerances of these species and may help to explain the previously observed patchiness of this weed's distribution in wheat fields.

Transcriptomic insights into shared responses to Fusarium crown rot infection and drought stresses in bread wheat (Triticum aestivum L.)

KEY MESSAGE

Shared changes in transcriptomes caused by Fusarium crown rot infection and drought stress were investigated based on a single pair of near-isogenic lines developed for a major locus conferring tolerance to both stresses. Fusarium crown rot (FCR) is a devastating disease in many areas of cereal production worldwide. It is well-known that drought stress enhances FCR severity but possible molecular relationship between these two stresses remains unclear. To investigate their relationships, we generated several pairs of near isogenic lines (NILs) targeting a locus conferring FCR resistance on chromosome 2D in bread wheat. One pair of these NILs showing significant differences between the two isolines for both FCR resistance and drought tolerance was used to investigate transcriptomic changes in responsive to these two stresses. Our results showed that the two isolines likely deployed different strategies in dealing with the stresses, and significant differences in expressed gene networks exist between the two time points of drought stresses evaluated in this study. Nevertheless, results from analysing Gene Ontology terms and transcription factors revealed that similar regulatory frameworks were activated in coping with these two stresses. Based on the position of the targeted locus, changes in expression following FCR infection and drought stresses, and the presence of non-synonymous variants between the two isolines, several candidate genes conferring resistance or tolerance to these two types of stresses were identified. The NILs generated, the large number of DEGs with single-nucleotide polymorphisms detected between the two isolines, and the candidate genes identified would be invaluable in fine mapping and cloning the gene(s) underlying the targeted locus.

A high-performance computational workflow to accelerate GATK SNP detection across a 25-genome dataset

BACKGROUND

Single-nucleotide polymorphisms (SNPs) are the most widely used form of molecular genetic variation studies. As reference genomes and resequencing data sets expand exponentially, tools must be in place to call SNPs at a similar pace. The genome analysis toolkit (GATK) is one of the most widely used SNP calling software tools publicly available, but unfortunately, high-performance computing versions of this tool have yet to become widely available and affordable.

RESULTS

Here we report an open-source high-performance computing genome variant calling workflow (HPC-GVCW) for GATK that can run on multiple computing platforms from supercomputers to desktop machines. We benchmarked HPC-GVCW on multiple crop species for performance and accuracy with comparable results with previously published reports (using GATK alone). Finally, we used HPC-GVCW in production mode to call SNPs on a "subpopulation aware" 16-genome rice reference panel with ~ 3000 resequenced rice accessions. The entire process took ~ 16 weeks and resulted in the identification of an average of 27.3 M SNPs/genome and the discovery of ~ 2.3 million novel SNPs that were not present in the flagship reference genome for rice (i.e., IRGSP RefSeq).

CONCLUSIONS

This study developed an open-source pipeline (HPC-GVCW) to run GATK on HPC platforms, which significantly improved the speed at which SNPs can be called. The workflow is widely applicable as demonstrated successfully for four major crop species with genomes ranging in size from 400 Mb to 2.4 Gb. Using HPC-GVCW in production mode to call SNPs on a 25 multi-crop-reference genome data set produced over 1.1 billion SNPs that were publicly released for functional and breeding studies. For rice, many novel SNPs were identified and were found to reside within genes and open chromatin regions that are predicted to have functional consequences. Combined, our results demonstrate the usefulness of combining a high-performance SNP calling architecture solution with a subpopulation-aware reference genome panel for rapid SNP discovery and public deployment.

New insights into molecular features of the genome-wide AOX family and their responses to various stresses in common wheat (Triticum aestivum L.)

Alternative oxidase (AOX) is an important terminal oxidase involved in the alternative oxidation pathway in plants, which is closely related to various biotic and abiotic stress responses. However, a comprehensive research on AOX gene family of wheat is still lacking. In this study, the members of wheat AOX (TaAOX) family were identified, and their molecular characteristics and gene expression patterns were systematically investigated. Seventeen TaAOX genes were identified from Chinese Spring (CS) genome, which were mapped on 7 chromosomes and mainly clustered on the long arm's distal end of the second homologous groups. Phylogenetic analysis showed that TaAOX genes were classified into four subgroups (Ia, Ib, Ic, and Id), and the Ia subgroup possessed the most members. Tandem duplication and segmental duplication events were found during the evolution of TaAOX genes and they were affected by purifying selection demonstrated by Ka/Ks analysis. The exon numbers of this family gene varied greatly from 1 to 9. Except for Ta3BSAOX14, all the proteins encoded by the other 16 TaAOX genes contained the amino acid residues of the key active sites in the AOX domain (cd01053). The expression patterns of TaAOX genes in various tissues and under abiotic and biotic stresses were analyzed using public transcriptome data, furthermore, qRT-PCR analysis was performed for some selected TaAOX genes, and the results suggested that most members of this gene family play an important role in response to different stresses in common wheat. Our results provide basic information and valuable reference for further exploring the gene function of TaAOX family by using gene editing, RNAi, VIGS, and other technologies.

RNAseq of diverse spring wheat cultivars released during last 110 years

Here, we performed RNA-seq based expression analysis of root and leaf tissues of a set of 24 historical spring wheat cultivars representing 110 years of temporal genetic variations. This huge 130 tissues RNAseq dataset was initially used to study expression pattern of 97 genes regulating root growth and development in wheat. Root system architecture (RSA) is an important target for breeding stress-resilient and high-yielding wheat cultivars under climatic fluctuations. However, root transcriptome analysis is usually obscured due to challenges in root research due to their below ground presence. We also validated the dataset by performing correlation analysis between expression of RSA related genes in roots and leaves with 25 root traits analyzed under varying moisture conditions and 10 yield-related traits. The Pearson's correlation coefficients between root phenotypes and expression of root-specific genes varied from -0.72 to 0.78, and strong correlations with genes such as DRO1, TaMOR, ARF4, PIN1 was observed. The presented datasets have multiple uses such as a) studying the change in expression pattern of genes during time, b) differential expression of genes in two very important tissues of wheat i.e., leaf and roots, and c) studying customized expression of genes associated with important phenotypes in diverse wheat cultivars. The initial findings presented here provided key insights into understanding the transcriptomic basis of phenotypic variability of RSA in wheat cultivars.

Identification of Whirly transcription factors in Triticeae species and functional analysis of TaWHY1-7D in response to osmotic stress

Osmotic stress poses a threat to the production and quality of crops. Whirly transcription factors have been investigated to enhance stress tolerance. In this study, a total of 18 Whirly genes were identified from six Triticeae species, which were classified into Whirly1 and Whirly2. The exon-intron structure, conserved motif, chromosomal location, collinearity, and regulatory network of Whirly genes were also analyzed. Real-time PCR results indicated that TaWHY1 genes exhibited higher expression levels in leaf sheaths and leaves during the seedling stage, while TaWHY2 genes were predominantly expressed in roots. Under PEG stress, the expression levels of TaWHY1-7A, TaWHY2-6A, TaWHY2-6B, and TaWHY2-6D were increased, TaWHY1-7D was reduced, and TaWHY1-4A had no significant change. All TaWHY genes were significantly up-regulated in response to NaCl stress treatment. In addition, TaWHY1-7A and TaWHY1-7D mainly enhanced the tolerance to oxidative stress in yeast cells. TaWHY2s mainly improved NaCl stress tolerance and were sensitive to oxidative stress in yeast cells. All TaWHYs slightly improved the yeast tolerance to d-sorbitol stress. The heterologous expression of TaWHY1-7D greatly improved drought and salt tolerance in transgenic Arabidopsis. In conclusion, these results provide the foundation for further functional study of Whirly genes aimed at improving osmotic stress tolerance in wheat.

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.

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.

Rht12b, a widely used ancient allele of TaGA2oxA13, reduces plant height and enhances yield potential in wheat

KEY MESSAGE

We identified a new wheat dwarfing allele Rht12b conferring reduced height and higher grain yield, pinpointed its causal variations, developed a breeding-applicable marker, and traced its origin and worldwide distribution. Plant height control is essential to optimize lodging resistance and yield gain in crops. RHT12 is a reduced height (Rht) locus that is identified in a mutationally induced dwarfing mutant and encodes a gibberellin 2-oxidase TaGA2oxA13. However, the artificial dwarfing allele is not used in wheat breeding due to excessive height reduction. Here, we confirmed a stable Rht locus, overlapping with RHT12, in a panel of wheat cultivars and its dwarfing allele reduced plant height by 5.4-8.2 cm, equivalent to Rht12b, a new allele of RHT12. We validated the effect of Rht12b on plant height in a bi-parent mapping population. Importantly, wheat cultivars carrying Rht12b had higher grain yield than those with the contrasting Rht12a allele. Rht12b conferred higher expression level of TaGA2oxA13. Transient activation assays defined SNP-390(C/A) in the promoter of TaGA2oxA13 as the causal variation. An efficient kompetitive allele-specific PCR marker was developed to diagnose Rht12b. Conjoint analysis showed that Rht12b plus the widely used Rht-D1b, Rht8 and Rht24b was the predominant Rht combination and conferred a moderate plant height in tested wheat cultivars. Evolutionary tracking uncovered that RHT12 locus arose from a tandem duplication event with Rht12b firstly appearing in wild emmer. The frequency of Rht12b was approximately 70% (700/1005) in a worldwide wheat panel and comparable to or higher than those of other widely used Rht genes, suggesting it had been subjected to positive selection. These findings not only identify a valuable Rht gene for wheat improvement but also develop a functionally diagnostic tool for marker-assisted breeding.

Deciphering spike architecture formation towards yield improvement in wheat

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

Deciphering the evolution and complexity of wheat germplasm from a genomic perspective

Bread wheat provides an essential fraction of the daily calorific intake for humanity. Due to its huge and complex genome, progress in studying on the wheat genome is substantially trailed behind those of the other two major crops, rice and maize, for at least a decade. With rapid advances in genome assembling and reduced cost of high-throughput sequencing, emerging de novo genome assemblies of wheat and whole-genome sequencing data are leading to a paradigm shift in wheat research. Here, we review recent progress in dissecting the complex genome and germplasm evolution of wheat since the release of the first high-quality wheat genome. New insights have been gained in the evolution of wheat germplasm during domestication and modern breeding progress, genomic variations at multiple scales contributing to the diversity of wheat germplasm, and complex transcriptional and epigenetic regulations of functional genes in polyploid wheat. Genomics databases and bioinformatics tools meeting the urgent needs of wheat genomics research are also summarized. The ever-increasing omics data, along with advanced tools and well-structured databases, are expected to accelerate deciphering the germplasm and gene resources in wheat for future breeding advances.

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

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

Genome-wide identification of ETHYLENE INSENSITIVE 2 in Triticeae species reveals that TaEIN2-4D.1 regulates cadmium tolerance in Triticum aestivum

ETHYLENE INSENSITIVE 2 (EIN2), as the core component of the ethylene signaling pathway, can widely regulate plant growth, development, and stress responses. However, the comprehensive study and function of EIN2 in wheat Cadmium (Cd) stress remain largely unexplored. Here, we identified 33 EIN2 genes and designated as TaEIN2-2B to TaEIN2-Un.3 in Triticum aestivum. The analysis of cis-regulatory elements in promoter regions and RNA-Seq showed that TaEIN2s were functionally related to plant growth and development, as well as the response to biotic and abiotic stress. qRT-PCR analysis of TaEIN2s indicated their sensitivity to Cd stress. Compared with WT plants, TaEIN2-4D.1-RNAi transgenic wheat lines showed enhanced shoot and root elongation, dry weight and chlorophyll accumulation, together with a reduced accumulation of Cd in wheat grain. In addition, TaEIN2-4D.1-RNAi transgenic wheat lines showed enhanced Reactive Oxygen Species (ROS) scavenging capacity compared with WT plants. In conclusion, our research indicates that TaEIN2 plays a key role in response to cadmium stress in wheat, which provides valuable information for crop improvement.

Thermosensitive SUMOylation of TaHsfA1 defines a dynamic ON/OFF molecular switch for the heat stress response in wheat

Dissecting genetic components in crop plants associated with heat stress (HS) sensing and adaptation will facilitate the design of modern crop varieties with improved thermotolerance. However, the molecular mechanisms underlying the ON/OFF switch controlling HS responses (HSRs) in wheat (Triticum aestivum) remain largely unknown. In this study, we focused on the molecular action of TaHsfA1, a class A heat shock transcription factor, in sensing dynamically changing HS signals and regulating HSRs. We show that the TaHsfA1 protein is modified by small ubiquitin-related modifier (SUMO) and that this modification is essential for the full transcriptional activation activity of TaHsfA1 in triggering downstream gene expression. During sustained heat exposure, the SUMOylation of TaHsfA1 is suppressed, which partially reduces TaHsfA1 protein activity, thereby reducing the intensity of downstream HSRs. In addition, we demonstrate that TaHsfA1 interacts with the histone acetyltransferase TaHAG1 in a thermosensitive manner. Together, our findings emphasize the importance of TaHsfA1 in thermotolerance in wheat. In addition, they define a highly dynamic SUMOylation-dependent "ON/OFF" molecular switch that senses temperature signals and contributes to thermotolerance in crops.

Multi-target genome editing reduces polyphenol oxidase activity in wheat (Triticum aestivum L.) grains

Introduction

Polyphenol oxidases (PPO) are dual activity metalloenzymes that catalyse the production of quinones. In plants, PPO activity may contribute to biotic stress resistance and secondary metabolism but is undesirable for food producers because it causes the discolouration and changes in flavour profiles of products during post-harvest processing. In wheat (Triticum aestivum L.), PPO released from the aleurone layer of the grain during milling results in the discolouration of flour, dough, and end-use products, reducing their value. Loss-of-function mutations in the PPO1 and PPO2 paralogous genes on homoeologous group 2 chromosomes confer reduced PPO activity in the wheat grain. However, limited natural variation and the proximity of these genes complicates the selection of extremely low-PPO wheat varieties by recombination. The goal of the current study was to edit all copies of PPO1 and PPO2 to drive extreme reductions in PPO grain activity in elite wheat varieties.

Results

A CRISPR/Cas9 construct with one single guide RNA (sgRNA) targeting a conserved copper binding domain was used to edit all seven PPO1 and PPO2 genes in the spring wheat cultivar 'Fielder'. Five of the seven edited T1 lines exhibited significant reductions in PPO activity, and T2 lines had PPO activity up to 86.7% lower than wild-type. The same construct was transformed into the elite winter wheat cultivars 'Guardian' and 'Steamboat', which have five PPO1 and PPO2 genes. In these varieties PPO activity was reduced by >90% in both T1 and T2 lines. In all three varieties, dough samples from edited lines exhibited reduced browning.

Discussion

This study demonstrates that multi-target editing at late stages of variety development could complement selection for beneficial alleles in crop breeding programs by inducing novel variation in loci inaccessible to recombination.

Integrative gene duplication and genome-wide analysis as an approach to facilitate wheat reverse genetics: An example in the TaCIPK family

INTRODUCTION

Reverse genetic studies conducted in the plant with a complex or polyploidy genome enriched with large gene families (like wheat) often meet challenges in identifying the key candidate genes related to important traits and prioritizing the genes for functional experiments.

OBJECTIVE

To overcome the above-mentioned challenges of reverse genetics, this work aims to establish an efficient multi-species strategy for genome-wide gene identification and prioritization of the key candidate genes.

METHODS

We established the integrative gene duplication and genome-wide analysis (iGG analysis) as a strategy for pinpointing key candidate genes deserving functional research. The iGG captures the evolution, and the expansion/contraction of large gene families across phylogeny-related species and integrates spatial-temporal expression information for gene function inference. Transgenic approaches were also employed to functional validation.

RESULTS

As a proof-of-concept for the iGG analysis, we took the wheat calcineurin B-like protein-interacting protein kinases (CIPKs) family as an example. We identified CIPKs from seven monocot species, established the orthologous relationship of CIPKs between rice and wheat, and characterized Triticeae-specific CIPK duplicates (e.g., CIPK4 and CIPK17). Integrated with our analysis of CBLs and CBL-CIPK interaction, we revealed that divergent expressions of TaCBLs and TaCIPKs could play an important role in keeping the stoichiometric balance of CBL-CIPK. Furthermore, we validated the function of TaCIPK17-A2 in the regulation of drought tolerance by using transgenic approaches. Overexpression of TaCIPK17 enhanced antioxidant capacity and improved drought tolerance in wheat.

CONCLUSION

The iGG analysis leverages evolutionary and comparative genomics of crops with large genomes to rapidly highlight the duplicated genes potentially associated with speciation, domestication and/or particular traits that deserve reverse-genetic functional studies. Through the identification of Triticeae-specific TaCIPK17 duplicates and functional validation, we demonstrated the effectiveness of the iGG analysis and provided a new target gene for improving drought tolerance in wheat.

A pathogen-induced putative NAC transcription factor mediates leaf rust resistance in barley

Leaf rust, caused by Puccinia hordei, is one of the most widespread and damaging foliar diseases affecting barley. The barley leaf rust resistance locus Rph7 has been shown to have unusually high sequence and haplotype divergence. In this study, we isolate the Rph7 gene using a fine mapping and RNA-Seq approach that is confirmed by mutational analysis and transgenic complementation. Rph7 is a pathogen-induced, non-canonical resistance gene encoding a protein that is distinct from other known plant disease resistance proteins in the Triticeae. Structural analysis using an AlphaFold2 protein model suggests that Rph7 encodes a putative NAC transcription factor with a zinc-finger BED domain with structural similarity to the N-terminal DNA-binding domain of the NAC transcription factor (ANAC019) from Arabidopsis. A global gene expression analysis suggests Rph7 mediates the activation and strength of the basal defence response. The isolation of Rph7 highlights the diversification of resistance mechanisms available for engineering disease control in crops.

Chromosome translocation affects multiple phenotypes, causes genome-wide dysregulation of gene expression, and remodels metabolome in hexaploid wheat

Chromosomal rearrangements (CRs) may occur in newly formed polyploids due to compromised meiotic fidelity. Moreover, CRs can be more readily tolerated in polyploids allowing their longer-term retention and hence potential spreading/fixation within a lineage. The direct functional consequences of CRs in plant polyploids remain unexplored. Here, we identified a heterozygous individual from a synthetic allohexaploid wheat in which the terminal parts of the long-arms of chromosomes 2D (approximately 193 Mb) and 4A (approximately 167 Mb) were reciprocally translocated. Five homogeneous translocation lines including both unbalanced and balanced types were developed by selfing fertilization of the founder mutant (RT [2DL; 4AL]-ter/1, reciprocal translocation). We investigated impacts of these translocations on phenotype, genome-wide gene expression and metabolome. We find that, compared with sibling wild-type, CRs in the form of both unbalanced and balanced translocations induced substantial changes of gene expression primarily via trans-regulation in the nascent allopolyploid wheat. The CRs also manifested clear phenotypic and metabolic consequences. In particular, the genetically balanced, stable reciprocal translocations lines showed immediate enhanced reproductive fitness relative to wild type. Our results underscore the profound impact of CRs on gene expression in nascent allopolyploids with wide-ranging phenotypic and metabolic consequences, suggesting CRs are an important source of genetic variation that can be exploited for crop breeding.

Pangenome-based trajectories of intracellular gene transfers in Poaceae unveil high cumulation in Triticeae

Intracellular gene transfers (IGTs) between the nucleus and organelles, including plastids and mitochondria, constantly reshape the nuclear genome during evolution. Despite the substantial contribution of IGTs to genome variation, the dynamic trajectories of IGTs at the pangenomic level remain elusive. Here, we developed an approach, IGTminer, that maps the evolutionary trajectories of IGTs using collinearity and gene reannotation across multiple genome assemblies. We applied IGTminer to create a nuclear organellar gene (NOG) map across 67 genomes covering 15 Poaceae species, including important crops. The resulting NOGs were verified by experiments and sequencing data sets. Our analysis revealed that most NOGs were recently transferred and lineage specific and that Triticeae species tended to have more NOGs than other Poaceae species. Wheat (Triticum aestivum) had a higher retention rate of NOGs than maize (Zea mays) and rice (Oryza sativa), and the retained NOGs were likely involved in photosynthesis and translation pathways. Large numbers of NOG clusters were aggregated in hexaploid wheat during 2 rounds of polyploidization, contributing to the genetic diversity among modern wheat accessions. We implemented an interactive web server to facilitate the exploration of NOGs in Poaceae. In summary, this study provides resources and insights into the roles of IGTs in shaping interspecies and intraspecies genome variation and driving plant genome evolution.

Identification of acyl-CoA-binding protein gene in Triticeae species reveals that TaACBP4A-1 and TaACBP4A-2 positively regulate powdery mildew resistance in wheat

Plant acyl-CoA-binding proteins (ACBPs), which contain the conserved ACB domain, participate in multiple biological processes, however, there are few reports on wheat ACBPs. In this study, the ACBP genes from nine different species were identified comprehensively. The expression patterns of TaACBP genes in multiple tissues and under various biotic stresses were determined by qRT-PCR. The function of selected TaACBP genes was studied by virus-induced gene silencing. A total of 67 ACBPs were identified from five monocotyledonous and four dicotyledonous species and divided into four classes. Tandem duplication analysis of the ACBPs suggested that tandem duplication events occurred in Triticum dicoccoides, but there was no tandem duplication event in wheat ACBP genes. Evolutionary analysis suggested that the TdACBPs may have experienced gene introgression during tetraploid evolution, while TaACBP gene loss events occurred during hexaploid wheat evolution. The expression pattern showed that all the TaACBP genes were expressed, and most of them were responsive to induction by Blumeria graminis f. sp. tritici or Fusarium graminearum. Silencing of TaACBP4A-1 and TaACBP4A-2 increased powdery mildew susceptibility in the common wheat BainongAK58. Furthermore, TaACBP4A-1, which belonged to class III, physically interacted with autophagy-related ubiquitin-like protein TaATG8g in yeast cells. This study provided a valuable reference for further investigations into the functional and molecular mechanisms of the ACBP gene family.

Centromere Plasticity With Evolutionary Conservation and Divergence Uncovered by Wheat 10+ Genomes

Centromeres (CEN) are the chromosomal regions that play a crucial role in maintaining genomic stability. The underlying highly repetitive DNA sequences can evolve quickly in most eukaryotes, and promote karyotype evolution. Despite their variability, it is not fully understood how these widely variable sequences ensure the homeostasis of centromere function. In this study, we investigated the genetics and epigenetics of CEN in a population of wheat lines from global breeding programs. We captured a high degree of sequences, positioning, and epigenetic variations in the large and complex wheat CEN. We found that most CENH3-associated repeats are Cereba element of retrotransposons and exhibit phylogenetic homogenization across different wheat lines, but the less-associated repeat sequences diverge on their own way in each wheat line, implying specific mechanisms for selecting certain repeat types as functional core CEN. Furthermore, we observed that CENH3 nucleosome structures display looser wrapping of DNA termini on complex centromeric repeats, including the repositioned CEN. We also found that strict CENH3 nucleosome positioning and intrinsic DNA features play a role in determining centromere identity among different lines. Specific non-B form DNAs were substantially associated with CENH3 nucleosomes for the repositioned centromeres. These findings suggest that multiple mechanisms were involved in the adaptation of CENH3 nucleosomes that can stabilize CEN. Ultimately, we proposed a remarkable epigenetic plasticity of centromere chromatin within the diverse genomic context, and the high robustness is crucial for maintaining centromere function and genome stability in wheat 10+ lines as a result of past breeding selections.

Analysis of Raffinose Synthase Gene Family in Bread Wheat and Identification of Drought Resistance and Salt Tolerance Function of TaRS15-3B

Raffinose synthase (RS) plays a crucial role in plant growth and development, as well as in responses to biotic stresses and abiotic stresses, yet few studies have been conducted on its role in bread wheat. Therefore, in this study we screened and identified a family of bread wheat raffinose synthase genes based on bread wheat genome information and analyzed their physicochemical properties, phylogenetic evolutionary relationships, conserved structural domains, promoter cis-acting elements, and expression patterns. The BSMV-induced silencing of TaRS15-3B resulted in the bread wheat seedlings being susceptible to drought and salt stress and reduced the expression levels of stress-related and ROS-scavenging genes in bread wheat plants. This further affected the ability of bread wheat to cope with drought and salt stress. In conclusion, this study revealed that the RS gene family in bread wheat plays an important role in plant response to abiotic stresses and that the TaRS15-3B gene can improve the tolerance of transgenic bread wheat to drought and salt stresses, provide directions for the study of other RS gene families in bread wheat, and supply candidate genes for use in molecular breeding of bread wheat for stress resistance.

Fine mapping and characterization of a major QTL for plant height on chromosome 5A in wheat

KEY MESSAGE

We precisely mapped QPH.caas-5AL for plant height in wheat, predicted candidate genes and confirmed genetic effects in a panel of wheat cultivars. Plant height is an important agronomic trait, and appropriately reduced height can improve yield potential and stability in wheat, usually combined with sufficient water and fertilizer. We previously detected a stable major-effect quantitative trait locus QPH.caas-5AL for plant height on chromosome 5A in a recombinant inbred line population of the cross 'Doumai × Shi 4185' using the wheat 90 K SNP assay. Here , QPH.caas-5AL was confirmed using new phenotypic data in additional environment and new-developed markers. We identified nine heterozygous recombinant plants for fine mapping of QPH.caas-5AL and developed 14 breeder-friendly kompetitive allele-specific PCR markers in the region of QPH.caas-5AL based on the genome re-sequencing data of parents. Phenotyping and genotyping analyses of secondary populations derived from the self-pollinated heterozygous recombinant plants delimited QPH.caas-5AL into an approximate 3.0 Mb physical region (521.0-524.0 Mb) according to the Chinese Spring reference genome. This region contains 45 annotated genes, and six of them were predicted as the candidates of QPH.caas-5AL based on genome and transcriptome sequencing analyses. We further validated that QPH.caas-5AL has significant effects on plant height but not yield component traits in a diverse panel of wheat cultivars; its dwarfing allele is frequently used in modern wheat cultivars. These findings lay a solid foundation for the map-based cloning of QPH.caas-5AL and also provide a breeding-applicable tool for its marker-assisted selection. Keymessage We precisely mapped QPH.caas-5AL for plant height in wheat, predicted candidate genes and confirmed genetic effects in a panel of wheat cultivars.

Genome-wide characterization of the VQ genes in Triticeae and their functionalization driven by polyploidization and gene duplication events in wheat

Valine-glutamine motif-containing (VQ) proteins are transcriptional cofactors widely involved in plant growth, development, and response to various stresses. Although the VQ family has been genome-wide identified in some species, but the knowledge regarding duplication-driven functionalization of VQ genes among evolutionarily related species is still lacking. Here, 952 VQ genes have been identified from 16 species, emphasizing seven Triticeae species including the bread wheat. Comprehensive phylogenetic and syntenic analyses allow us to establish the orthologous relationship of VQ genes from rice (Oryza sativa) to bread wheat (Triticum aestivum). The evolutionary analysis revealed that whole-genome duplication (WGD) drives the expansion of OsVQs, while TaVQs expansion is associated with a recent burst of gene duplication (RBGD). We also analyzed the motif composition and molecular properties of TaVQ proteins, enriched biological functions, and expression patterns of TaVQs. We demonstrate that WGD-derived TaVQs have become divergent in both protein motif composition and expression pattern, while RBGD-derived TaVQs tend to adopt specific expression patterns, suggesting their functionalization in certain biological processes or in response to specific stresses. Furthermore, some RBGD-derived TaVQs are found to be associated with salt tolerance. Several of the identified salt-related TaVQ proteins were located in the cytoplasm and nucleus and their salt-responsive expression patterns were validated by qPCR analysis. Yeast-based functional experiments confirmed that TaVQ27 may be a new regulator to salt response and regulation. Overall, this study lays the foundation for further functional validation of VQ family members within the Triticeae species.

Identification of candidate gene for the defective kernel phenotype using bulked segregant RNA and exome capture sequencing methods in wheat

Wheat is a significant source of protein and starch worldwide. The defective kernel (Dek) mutant AK-3537, displaying a large hollow area in the endosperm and shrunken grain, was obtained through ethyl methane sulfonate (EMS) treatment of the wheat cultivar Aikang 58 (AK58). The mode of inheritance of the AK-3537 grain Dek phenotype was determined to be recessive with a specific statistical significance level. We used bulked segregant RNA-seq (BSR-seq), BSA-based exome capture sequencing (BSE-seq), and the ΔSNP-index algorithm to identify candidate regions for the grain Dek phenotype. Two major candidate regions, DCR1 (Dek candidate region 1) and DCR2, were identified on chromosome 7A between 279.98 and 287.93 Mb and 565.34 and 568.59 Mb, respectively. Based on transcriptome analysis and previous reports, we designed KASP genotyping assays based on SNP variations in the candidate regions and speculated that the candidate gene is TraesCS7A03G0625900 (HMGS-7A), which encodes a 3-hydroxy-3-methylglutaryl-CoA synthase. One SNP variation located at position 1,049 in the coding sequence (G>A) causes an amino acid change from Gly to Asp. The research suggests that functional changes in HMGS-7A may affect the expression of key enzyme genes involved in wheat starch syntheses, such as GBSSII and SSIIIa.

The translational landscape of bread wheat during grain development

The dynamics of gene expression in crop grains has typically been investigated at the transcriptional level. However, this approach neglects translational regulation, a widespread mechanism that rapidly modulates gene expression to increase the plasticity of organisms. Here, we performed ribosome profiling and polysome profiling to obtain a comprehensive translatome data set of developing bread wheat (Triticum aestivum) grains. We further investigated the genome-wide translational dynamics during grain development, revealing that the translation of many functional genes is modulated in a stage-specific manner. The unbalanced translation between subgenomes is pervasive, which increases the expression flexibility of allohexaploid wheat. In addition, we uncovered widespread previously unannotated translation events, including upstream open reading frames (uORFs), downstream open reading frames (dORFs), and open reading frames (ORFs) in long noncoding RNAs, and characterized the temporal expression dynamics of small ORFs. We demonstrated that uORFs act as cis-regulatory elements that can repress or even enhance the translation of mRNAs. Gene translation may be combinatorially modulated by uORFs, dORFs, and microRNAs. In summary, our study presents a translatomic resource that provides a comprehensive and detailed overview of the translational regulation in developing bread wheat grains. This resource will facilitate future crop improvements for optimal yield and quality.

Genome-Wide Analysis of MYB Transcription Factors in the Wheat Genome and Their Roles in Salt Stress Response

Large and rapidly increasing areas of salt-affected soils are posing major challenges for the agricultural sector. Most fields used for the important food crop Triticum aestivum (wheat) are expected to be salt-affected within 50 years. To counter the associated problems, it is essential to understand the molecular mechanisms involved in salt stress responses and tolerance, thereby enabling their exploitation in the development of salt-tolerant varieties. The myeloblastosis (MYB) family of transcription factors are key regulators of responses to both biotic and abiotic stress, including salt stress. Thus, we used the Chinese spring wheat genome assembled by the International Wheat Genome Sequencing Consortium to identify putative MYB proteins (719 in total). Protein families (PFAM) analysis of the MYB sequences identified 28 combinations of 16 domains in the encoded proteins. The most common consisted of MYB_DNA-binding and MYB-DNA-bind_6 domains, and five highly conserved tryptophans were located in the aligned MYB protein sequence. Interestingly, we found and characterized a novel 5R-MYB group in the wheat genome. In silico studies showed that MYB transcription factors MYB3, MYB4, MYB13 and MYB59 are involved in salt stress responses. qPCR analysis confirmed upregulation of the expression of all these MYBs in both roots and shoots of the wheat variety BARI Gom-25 (except MYB4, which was downregulated in roots) under salt stress. Moreover, we identified nine target genes involved in salt stress that are regulated by the four MYB proteins, most of which have cellular locations and are involved in catalytic and binding activities associated with various cellular and metabolic processes.

Low impact of polyploidization on the transcriptome of synthetic allohexaploid wheat

BACKGROUND

Bread wheat is a recent allohexaploid (genomic constitution AABBDD) that emerged through a hybridization between tetraploid Triticum turgidum (AABB) and diploid Aegilops tauschii (DD) less than 10,000 years ago. The hexaploidization can be re-created artificially, producing synthetic wheat that has been used to study immediate genomic responses to polyploidization. The scale of the consequences of polyploidization, and their mechanism of establishment, remain uncertain.

RESULTS

Here we sampled several synthetic wheats from alternative parental genotypes and reciprocal crosses, and examined transcriptomes from two different tissues and successive generations. We did not detect any massive reprogramming in gene expression, with only around 1% of expressed genes showing significant differences compared to their lower-ploidy parents. Most of this differential expression is located on the D subgenome, without consistency in the direction of the expression change. Homoeolog expression bias in synthetic wheat is similar to the pattern observed in the parents. Both differential expression and homoeolog bias are tissue-specific. While up to three families of transposable elements became upregulated in wheat synthetics, their position and distance are not significantly associated with expression changes in proximal genes.

DISCUSSION

While only a few genes change their expression pattern after polyploidization, they can be involved in agronomically important pathways. Alternative parental combinations can lead to opposite changes on the same subset of D-located genes, which is relevant for harnessing new diversity in wheat breeding. Tissue specificity of the polyploidization-triggered expression changes indicates the remodelling of transcriptomes in synthetic wheat is plastic and likely caused by regulome interactions rather than permanent changes. We discuss the pitfalls of transcriptomic comparisons across ploidy levels that can inflate the de-regulation signal.

CONCLUSIONS

Transcriptomic response to polyploidization in synthetic AABBDD wheat is modest and much lower than some previous estimates. Homoeolog expression bias in wheat allohexaploids is mostly attributed to parental legacy, with polyploidy having a mild balancing effect.

Phylogenomic discovery of deleterious mutations facilitates hybrid potato breeding

Hybrid potato breeding will transform the crop from a clonally propagated tetraploid to a seed-reproducing diploid. Historical accumulation of deleterious mutations in potato genomes has hindered the development of elite inbred lines and hybrids. Utilizing a whole-genome phylogeny of 92 Solanaceae and its sister clade species, we employ an evolutionary strategy to identify deleterious mutations. The deep phylogeny reveals the genome-wide landscape of highly constrained sites, comprising ∼2.4% of the genome. Based on a diploid potato diversity panel, we infer 367,499 deleterious variants, of which 50% occur at non-coding and 15% at synonymous sites. Counterintuitively, diploid lines with relatively high homozygous deleterious burden can be better starting material for inbred-line development, despite showing less vigorous growth. Inclusion of inferred deleterious mutations increases genomic-prediction accuracy for yield by 24.7%. Our study generates insights into the genome-wide incidence and properties of deleterious mutations and their far-reaching consequences for breeding.

Reducing brassinosteroid signalling enhances grain yield in semi-dwarf wheat

Modern green revolution varieties of wheat (Triticum aestivum L.) confer semi-dwarf and lodging-resistant plant architecture owing to the Reduced height-B1b (Rht-B1b) and Rht-D1b alleles¹. However, both Rht-B1b and Rht-D1b are gain-of-function mutant alleles encoding gibberellin signalling repressors that stably repress plant growth and negatively affect nitrogen-use efficiency and grain filling^2-5. Therefore, the green revolution varieties of wheat harbouring Rht-B1b or Rht-D1b usually produce smaller grain and require higher nitrogen fertilizer inputs to maintain their grain yields. Here we describe a strategy to design semi-dwarf wheat varieties without the need for Rht-B1b or Rht-D1b alleles. We discovered that absence of Rht-B1 and ZnF-B (encoding a RING-type E3 ligase) through a natural deletion of a haploblock of about 500 kilobases shaped semi-dwarf plants with more compact plant architecture and substantially improved grain yield (up to 15.2%) in field trials. Further genetic analysis confirmed that the deletion of ZnF-B induced the semi-dwarf trait in the absence of the Rht-B1b and Rht-D1b alleles through attenuating brassinosteroid (BR) perception. ZnF acts as a BR signalling activator to facilitate proteasomal destruction of the BR signalling repressor BRI1 kinase inhibitor 1 (TaBKI1), and loss of ZnF stabilizes TaBKI1 to block BR signalling transduction. Our findings not only identified a pivotal BR signalling modulator but also provided a creative strategy to design high-yield semi-dwarf wheat varieties by manipulating the BR signal pathway to sustain wheat production.

Global transcriptome analysis of allopolyploidization reveals large-scale repression of the D-subgenome in synthetic hexaploid wheat

Synthetic hexaploid wheat (SHW) lines are created as pre-breeding germplasm to diversify the D subgenome of hexaploid wheat and capitalize upon the untapped genetic diversity of the Aegilops tauschii gene pool. However, the phenotypes observed in the Ae. tauschii parents are not always recovered in the SHW lines, possibly due to inter-subgenome interactions. To elucidate this post-polyploidization genome reprogramming phenomenon, we performed RNA-seq of four SHW lines and their corresponding tetraploid and diploid parents, across ten tissues and three biological replicates. Homoeologue expression bias (HEB) analysis using more than 18,000 triads suggests massive suppression of homoeoalleles of the D subgenome in SHWs. Comparative transcriptome analysis of the whole-genome gene set further corroborated this finding. Alternative splicing analysis of the high-confidence genes indicates an additional layer of complexity where all five splice events are identified, and retained intron is predominant. Homoeologue expression upon resynthesis of hexaploid wheat has implications to the usage and handling of this germplasm in breeding as it relates to capturing the effects of epistatic interaction across subgenomes upon polyploidization. Special considerations must be given to this germplasm in pre-breeding activities to consider the extent of the inter-subgenome interactions on gene expression and their impact on traits for crop improvement.

Transcriptome Analysis of Roots from Wheat (Triticum aestivum L.) Varieties in Response to Drought Stress

Under climate change, drought is one of the most limiting factors that influences wheat (Triticum aestivum L.) production. Exploring stress-related genes is vital for wheat breeding. To identify genes related to the drought tolerance response, two common wheat cultivars, Zhengmai 366 (ZM366) and Chuanmai 42 (CM42), were selected based on their obvious difference in root length under 15% PEG-6000 treatment. The root length of the ZM366 cultivar was significantly longer than that of CM42. Stress-related genes were identified by RNA-seq in samples treated with 15% PEG-6000 for 7 days. In total, 11,083 differentially expressed genes (DEGs) and numerous single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) were identified. GO enrichment analysis revealed that the upregulated genes were mainly related to the response to water, acidic chemicals, oxygen-containing compounds, inorganic substances, and abiotic stimuli. Among the DEGs, the expression levels of 16 genes in ZM366 were higher than those in CM42 after the 15% PEG-6000 treatment based on RT-qPCR. Furthermore, EMS-induced mutants in Kronos (T. turgidum L.) of 4 representative DEGs possessed longer roots than the WT after the 15% PEG-6000 treatment. Altogether, the drought stress genes identified in this study represent useful gene resources for wheat breeding.

Comparative analysis of Fusarium crown rot resistance in synthetic hexaploid wheats and their parental genotypes

BACKGROUND

Fusarium crown rot (FCR) is a chronic disease of cereals worldwide. Compared with tetraploid wheat, hexaploid wheat is more resistant to FCR infection. The underlying reasons for the differences are still not clear. In this study, we compared FCR responses of 10 synthetic hexaploid wheats (SHWs) and their tetraploid and diploid parents. We then performed transcriptome analysis to uncover the molecular mechanism of FCR on these SHWs and their parents.

RESULTS

We observed higher levels of FCR resistance in the SHWs compared with their tetraploid parents. The transcriptome analysis suggested that multiple defense pathways responsive to FCR infection were upregulated in the SHWs. Notably, phenylalanine ammonia lyase (PAL) genes, involved in lignin and salicylic acid (SA) biosynthesis, exhibited a higher level of expression to FCR infection in the SHWs. Physiological and biochemical analysis validated that PAL activity and SA and lignin contents of the stem bases were higher in SHWs than in their tetraploid parents.

CONCLUSION

Overall, these findings imply that improved FCR resistance in SHWs compared with their tetraploid parents is probably related to higher levels of response on PAL-mediated lignin and SA biosynthesis pathways.