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

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.

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 alleles1. 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 filling2-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.

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.

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.

WheatCENet: A Database for Comparative Co-expression Networks Analysis of Allohexaploid Wheat and Its Progenitors

Genetic and epigenetic changes after polyploidization events could result in variable gene expression and modified regulatory networks. Here, using large-scale transcriptome data, we constructed co-expression networks for diploid, tetraploid, and hexaploid wheat species, and built a platform for comparing co-expression networks of allohexaploid wheat and its progenitors, named WheatCENet. WheatCENet is a platform for searching and comparing specific functional co-expression networks, as well as identifying the related functions of the genes clustered therein. Functional annotations like pathways, gene families, protein-protein interactions, microRNAs (miRNAs), and several lines of epigenome data are integrated into this platform, and Gene Ontology (GO) annotation, gene set enrichment analysis (GSEA), motif identification, and other useful tools are also included. Using WheatCENet, we found that the network of WHEAT ABERRANT PANICLE ORGANIZATION 1 (WAPO1) has more co-expressed genes related to spike development in hexaploid wheat than its progenitors. We also found a novel motif of CCWWWWWWGG (CArG) specifically in the promoter region of WAPO-A1, suggesting that neofunctionalization of the WAPO-A1 gene affects spikelet development in hexaploid wheat. WheatCENet is useful for investigating co-expression networks and conducting other analyses, and thus facilitates comparative and functional genomic studies in wheat. WheatCENet is freely available at http://bioinformatics.cpolar.cn/WheatCENet and http://bioinformatics.cau.edu.cn/WheatCENet.

Fine mapping of reduced height locus RHT26 in common wheat

We fine mapped RHT26 for plant height in wheat, confirmed its genetic effects in a panel of wheat cultivars and predicted candidate genes. Development of wheat cultivars with appropriate plant height (PH) is an important goal in breeding. Utilization of semi-dwarfing genes Rht-B1b and Rht-D1b triggered wheat Green Resolution in the 1960s. Since these genes also bring unfavorable features, such as reduced coleoptile length and grain weight, it is necessary to identify alternative reduced height genes without yield penalty. Here we constructed a high-density genetic map of a recombinant inbred line population derived from the cross of Zhongmai175 and Lunxuan987 and detected a stable genetic locus for PH, designated RHT26, on chromosome arm 3DL in all of six environments, accounting for 6.8-14.0% of the phenotypic variances. RHT26 was delimited to an approximate 1.4 Mb physical interval (517.1-518.5 Mb) using secondary mapping populations derived from 22 heterozygous recombinant plants and 24 kompetitive allele-specific PCR markers. Eleven high-confidence genes were annotated in the physical interval according to the Chinese Spring reference genome, and four of them were predicted as candidates for RHT26 based on genome and transcriptome sequencing analyses. We also confirmed that RHT26 had significant effects on PH, but not grain yield in a panel of wheat cultivars; its dwarfing allele has been frequently used in wheat breeding. These findings lay a sound foundation for map-based cloning of RHT26 and provide a breeding-applicable tool for marker-assisted selection.

The TaTCP4/10-B1 cascade regulates awn elongation in wheat (Triticum aestivum L.)

Awns are important morphological markers for wheat and exert a strong physiological effect on wheat yield. The awn elongation suppressor B1 has recently been cloned through association and linkage analysis in wheat. However, the mechanism of awn inhibition centered around B1 remains to be clarified. Here, we identified an allelic variant in the coding region of B1 through analysis of re-sequencing data; this variant causes an amino acid substitution and premature termination, resulting in a long-awn phenotype. Transcriptome analysis indicated that B1 inhibited awn elongation by impeding cytokinin- and auxin-promoted cell division. Moreover, B1 directly repressed the expression of TaRAE2 and TaLks2, whose orthologs have been reported to promote awn development in rice or barley. More importantly, we found that TaTCP4 and TaTCP10 synergistically inhibited the expression of B1, and a G-to-A mutation in the B1 promoter attenuated its inhibition by TaTCP4/10. Taken together, our results reveal novel mechanisms of awn development and provide genetic resources for trait improvement in wheat.

A Survey of the Transcriptomic Resources in Durum Wheat: Stress Responses, Data Integration and Exploitation

Durum wheat (Triticum turgidum subsp. durum (Desf.) Husn.) is an allotetraploid cereal crop of worldwide importance, given its use for making pasta, couscous, and bulgur. Under climate change scenarios, abiotic (e.g., high and low temperatures, salinity, drought) and biotic (mainly exemplified by fungal pathogens) stresses represent a significant limit for durum cultivation because they can severely affect yield and grain quality. The advent of next-generation sequencing technologies has brought a huge development in transcriptomic resources with many relevant datasets now available for durum wheat, at various anatomical levels, also focusing on phenological phases and environmental conditions. In this review, we cover all the transcriptomic resources generated on durum wheat to date and focus on the corresponding scientific insights gained into abiotic and biotic stress responses. We describe relevant databases, tools and approaches, including connections with other "omics" that could assist data integration for candidate gene discovery for bio-agronomical traits. The biological knowledge summarized here will ultimately help in accelerating durum wheat breeding.

Positional cloning and characterization reveal the role of a miRNA precursor gene ZmLRT in the regulation of lateral root number and drought tolerance in maize

Lateral roots play essential roles in drought tolerance in maize (Zea mays L.). However, the genetic basis for the variation in the number of lateral roots in maize remains elusive. Here, we identified a major quantitative trait locus (QTL), qLRT5-1, controlling lateral root number using a recombinant inbred population from a cross between the maize lines Zong3 (with many lateral roots) and 87-1 (with few lateral roots). Fine-mapping and functional analysis determined that the candidate gene for qLRT5-1, ZmLRT, expresses the primary transcript for the microRNA miR166a. ZmLRT was highly expressed in root tips and lateral root primordia, and knockout and overexpression of ZmLRT increased and decreased lateral root number, respectively. Compared with 87-1, the ZmLRT gene model of Zong3 lacked the second and third exons and contained a 14 bp deletion at the junction between the first exon and intron, which altered the splicing site. In addition, ZmLRT expression was significantly lower in Zong3 than in 87-1, which might be attributed to the insertions of a transposon and over large DNA fragments in the Zong3 ZmLRT promoter region. These mutations decreased the abundance of mature miR166a in Zong3, resulting in increased lateral roots at the seedling stage. Furthermore, miR166a post-transcriptionally repressed five development-related class-III homeodomain-leucine zipper genes. Moreover, knockout of ZmLRT enhanced drought tolerance of maize seedlings. Our study furthers our understanding of the genetic basis of lateral root number variation in maize and highlights ZmLRT as a target for improving drought tolerance in maize.

Genome-Wide Investigation of Apyrase (APY) Genes in Peanut (Arachis hypogaea L.) and Functional Characterization of a Pod-Abundant Expression Promoter AhAPY2-1p

Peanut (Arachis hypogaea L.) is an important food and feed crop worldwide and is affected by various biotic and abiotic stresses. The cellular ATP levels decrease significantly during stress as ATP molecules move to extracellular spaces, resulting in increased ROS production and cell apoptosis. Apyrases (APYs) are the nucleoside phosphatase (NPTs) superfamily members and play an important role in regulating cellular ATP levels under stress. We identified 17 APY homologs in A. hypogaea (AhAPYs), and their phylogenetic relationships, conserved motifs, putative miRNAs targeting different AhAPYs, cis-regulatory elements, etc., were studied in detail. The transcriptome expression data were used to observe the expression patterns in different tissues and under stress conditions. We found that the AhAPY2-1 gene showed abundant expression in the pericarp. As the pericarp is a key defense organ against environmental stress and promoters are the key elements regulating gene expression, we functionally characterized the AhAPY2-1 promoter for its possible use in future breeding programs. The functional characterization of AhAPY2-1P in transgenic Arabidopsis plants showed that it effectively regulated GUS gene expression in the pericarp. GUS expression was also detected in flowers of transgenic Arabidopsis plants. Overall, these results strongly suggest that APYs are an important future research subject for peanut and other crops, and AhPAY2-1P can be used to drive the resistance-related genes in a pericarp-specific manner to enhance the defensive abilities of the pericarp.

TaMADS29 interacts with TaNF-YB1 to synergistically regulate early grain development in bread wheat

Grain development is a crucial determinant of yield and quality in bread wheat (Triticum aestivum L.). However, the regulatory mechanisms underlying wheat grain development remain elusive. Here we report how TaMADS29 interacts with TaNF-YB1 to synergistically regulate early grain development in bread wheat. The tamads29 mutants generated by CRISPR/Cas9 exhibited severe grain filling deficiency, coupled with excessive accumulation of reactive oxygen species (ROS) and abnormal programmed cell death that occurred in early developing grains, while overexpression of TaMADS29 increased grain width and 1,000-kernel weight. Further analysis revealed that TaMADS29 interacted directly with TaNF-YB1; null mutation in TaNF-YB1 caused grain developmental deficiency similar to tamads29 mutants. The regulatory complex composed of TaMADS29 and TaNF-YB1 exercises its possible function that inhibits the excessive accumulation of ROS by regulating the genes involved in chloroplast development and photosynthesis in early developing wheat grains and prevents nucellar projection degradation and endosperm cell death, facilitating transportation of nutrients into the endosperm and wholly filling of developing grains. Collectively, our work not only discloses the molecular mechanism of MADS-box and NF-Y TFs in facilitating bread wheat grain development, but also indicates that caryopsis chloroplast might be a central regulator of grain development rather than merely a photosynthesis organelle. More importantly, our work offers an innovative way to breed high-yield wheat cultivars by controlling the ROS level in developing grains.

A chromosome-scale genome assembly of Dasypyrum villosum provides insights into its application as a broad-spectrum disease resistance resource for wheat improvement

Dasypyrum villosum is one of the most valuable gene resources in wheat improvement, especially for disease resistance. The mining of favorable genes from D. villosum is frustrated by the lack of a whole genome sequence. In this study, we generated a doubled-haploid line, 91C43^DH, using microspore culture and obtained a 4.05-GB high-quality, chromosome-scale genome assembly for D. villosum. The assembly contains39 727 high-confidence genes, and 85.31% of the sequences are repetitive. Two reciprocal translocation events were detected, and 7VS-4VL is a unique translocation in D. villosum. The prolamin seed storage protein-coding genes were found to be duplicated; in particular, the genes encoding low-molecular-weight glutenin at the Glu-V3 locus were significantly expanded. RNA sequencing (RNA-seq) analysis indicated that, after Blumeria graminearum f.sp tritici (Bgt) inoculation, there were more upregulated genes involved in the pattern-triggered immunity and effector-triggered immunity defense pathways in D. villosum than in Triticum urartu. MNase hypersensitive sequencing (MH-seq) identified two Bgt-inducible MH sites (MHSs), one in the promoter and one in the 3' terminal region of the powdery mildew resistance (Pm) gene NLR1-V. Each site had two subpeaks and they were termed MHS1 (MHS1.1/1.2) and MHS2 (MHS2.1/2.2). Bgt-inducible MHS2.2 was uniquely present in D. villosum, and MHS1.1 was more inducible in D. villosum than in wheat, suggesting that MHSs may be critical for regulation of NLR1-V expression and plant defense. In summary, this study provides a valuable genome resource for functional genomics studies and wheat-D. villosum introgression breeding. The identified regulatory mechanisms may also be exploited to develop new strategies for enhancing Pm resistance by optimizing gene expression in wheat.

Comparative gene retention analysis in barley, wild emmer, and bread wheat pangenome lines reveals factors affecting gene retention following gene duplication

BACKGROUND

Gene duplication is a prevalent phenomenon and a major driving force underlying genome evolution. The process leading to the fixation of gene duplicates following duplication is critical to understand how genome evolves but remains fragmentally understood. Most previous studies on gene retention are based on gene duplicate analyses in single reference genome. No population-based comparative gene retention analysis has been performed to date.

RESULTS

Taking advantage of recently published genomic data in Triticeae, we dissected a divergent homogentisate phytyltransferase (HPT2) lineage caught in the middle stage of gene fixation following duplication. The presence/absence of HPT2 in barley (diploid), wild emmer (tetraploid), and bread wheat (hexaploid) pangenome lines appears to be associated with gene dosage constraint and environmental adaption. Based on these observations, we adopted a phylogeny-based orthology inference approach and performed comparative gene retention analyses across barley, wild emmer, and bread wheat. This led to the identification of 326 HPT2-pattern-like genes at whole genome scale, representing a pool of gene duplicates in the middle stage of gene fixation. Majority of these HPT2-pattern-like genes were identified as small-scale duplicates, such as dispersed, tandem, and proximal duplications. Natural selection analyses showed that HPT2-pattern-like genes have experienced relaxed selection pressure, which is generally accompanied with partial positive selection and transcriptional divergence. Functional enrichment analyses showed that HPT2-pattern-like genes are over-represented with molecular-binding and defense response functions, supporting the potential role of environmental adaption during gene retention. We also observed that gene duplicates from larger gene family are more likely to be lost, implying a gene dosage constraint effect. Further comparative gene retention analysis in barley and bread wheat pangenome lines revealed combined effects of species-specific selection and gene dosage constraint.

CONCLUSIONS

Comparative gene retention analyses at the population level support gene dosage constraint, environmental adaption, and species-specific selection as three factors that may affect gene retention following gene duplication. Our findings shed light on the evolutionary process leading to the retention of newly formed gene duplicates and will greatly improve our understanding on genome evolution via duplication.

A wheat integrative regulatory network from large-scale complementary functional datasets enables trait-associated gene discovery for crop improvement

Gene regulation is central to all aspects of organism growth, and understanding it using large-scale functional datasets can provide a whole view of biological processes controlling complex phenotypic traits in crops. However, the connection between massive functional datasets and trait-associated gene discovery for crop improvement is still lacking. In this study, we constructed a wheat integrative gene regulatory network (wGRN) by combining an updated genome annotation and diverse complementary functional datasets, including gene expression, sequence motif, transcription factor (TF) binding, chromatin accessibility, and evolutionarily conserved regulation. wGRN contains 7.2 million genome-wide interactions covering 5947 TFs and 127 439 target genes, which were further verified using known regulatory relationships, condition-specific expression, gene functional information, and experiments. We used wGRN to assign genome-wide genes to 3891 specific biological pathways and accurately prioritize candidate genes associated with complex phenotypic traits in genome-wide association studies. In addition, wGRN was used to enhance the interpretation of a spike temporal transcriptome dataset to construct high-resolution networks. We further unveiled novel regulators that enhance the power of spike phenotypic trait prediction using machine learning and contribute to the spike phenotypic differences among modern wheat accessions. Finally, we developed an interactive webserver, wGRN (http://wheat.cau.edu.cn/wGRN), for the community to explore gene regulation and discover trait-associated genes. Collectively, this community resource establishes the foundation for using large-scale functional datasets to guide trait-associated gene discovery for crop improvement.

Plant pan-genomics and its applications

Plant genomes are so highly diverse that a substantial proportion of genomic sequences are not shared among individuals. The variable DNA sequences, along with the conserved core sequences, compose the more sophisticated pan-genome that represents the collection of all non-redundant DNA in a species. With rapid progress in genome sequencing technologies, pan-genome research in plants is now accelerating. Here we review recent advances in plant pan-genomics, including major driving forces of structural variations that constitute the variable sequences, methodological innovations for representing the pan-genome, and major successes in constructing plant pan-genomes. We also summarize recent efforts toward decoding the remaining dark matter in telomere-to-telomere or gapless plant genomes. These new genome resources, which have remarkable advantages over numerous previously assembled less-than-perfect genomes, are expected to become new references for genetic studies and plant breeding.

Gapless genome assembly of azalea and multi-omics investigation into divergence between two species with distinct flower color

The genus Rhododendron (Ericaceae), with more than 1000 species highly diverse in flower color, is providing distinct ornamental values and a model system for flower color studies. Here, we investigated the divergence between two parental species with different flower color widely used for azalea breeding. Gapless genome assembly was generated for the yellow-flowered azalea, Rhododendron molle. Comparative genomics found recent proliferation of long terminal repeat retrotransposons (LTR-RTs), especially Gypsy, has resulted in a 125 Mb (19%) genome size increase in species-specific regions, and a significant amount of dispersed gene duplicates (13 402) and pseudogenes (17 437). Metabolomic assessment revealed that yellow flower coloration is attributed to the dynamic changes of carotenoids/flavonols biosynthesis and chlorophyll degradation. Time-ordered gene co-expression networks (TO-GCNs) and the comparison confirmed the metabolome and uncovered the specific gene regulatory changes underpinning the distinct flower pigmentation. B3 and ERF TFs were found dominating the gene regulation of carotenoids/flavonols characterized pigmentation in R. molle, while WRKY, ERF, WD40, C2H2, and NAC TFs collectively regulated the anthocyanins characterized pigmentation in the red-flowered R simsii. This study employed a multi-omics strategy in disentangling the complex divergence between two important azaleas and provided references for further functional genetics and molecular breeding.

Positional cloning identified HvTUBULIN8 as the candidate gene for round lateral spikelet (RLS) in barley (Hordeum vulgare L.)

Map-based cloning, subcellular localization, virus-induced-gene-silencing and transcriptomic analysis reveal HvTUB8 as a candidate gene with pleiotropic effects on barley spike and leaf development via ethylene and chlorophyll metabolism. Barley lateral spikelet morphology and grain shape play key roles in grain physical quality and yield. Several genes and QTLs for these traits have been cloned or fine mapped previously. Here, we report the phenotypic and genotypic analysis of a barley mutant with round lateral spikelet (rls) from cv. Edamai 934. rls had round lateral spikelet, short but round grain, shortened awn, thick glume and dark green leaves. Histocytologic and ultrastructural analysis revealed that the difference of grain shape of rls was caused by change of cell arrangement in glume, and the dark leaf color resulted from enlarged chloroplast. HvTUBULIN8 (HvTUB8) was identified as the candidate gene for rls by combination of RNA-Seq, map-based-cloning, virus-induced-gene-silencing (VIGS) and protein subcellular location. A single G-A substitution at the third exon of HvTUB8 resulted in change of Cysteine 354 to tyrosine. Furthermore, the mutant isoform Hvtub8 could be detected in both nucleus and cytoplasm, whereas the wild-type protein was only in cytoplasm and granular organelles of wheat protoplasts. Being consistent with the rare phenotype, the "A" allele of HvTUB8 was only detected in rls, but not in a worldwide barley germplasm panel with 400 accessions. VIGS confirmed that HvTUB8 was essential to maintain spike integrity. RNA-Seq results suggested that HvTUB8 may control spike morphogenesis via ethylene homeostasis and signaling, and control leaf color through chlorophyll metabolism. Collectively, our results support HvTUB8 as a candidate gene for barley spike and leaf morphology and provide insight of a novel mechanism of it in barley development.

Genetic dissection of quantitative trait loci for flag leaf size in bread wheat (Triticum aestivum L.)

Flag leaf size is a crucial trait influencing plant architecture and yield potential in wheat. A recombinant inbred line (RIL) population derived from the cross of W7268 and Chuanyu 12 was employed to identify quantitative trait loci (QTL) controlling flag leaf length (FLL), flag leaf width (FLW), and flag leaf area (FLA) in six environments and the best linear unbiased estimator (BLUE) datasets. Using a 55 K SNP-based genetic map, six major and stable QTL were detected with 6.33-53.12% of explained phenotypic variation. Except for QFlw.cib-4B.3, the other five major QTL were co-located within two intervals on chromosomes 2B and 2D, namely QFll/Fla.cib-2B and QFll/Flw/Fla.cib-2D, respectively. Their interactions and effects on the corresponding traits and yield-related traits were also assessed based on flanking markers. QFll/Fla.cib-2B showed pleiotropic effects on spikelet number per spike (SNS). QFlw.cib-4B.3 and QFll/Flw/Fla.cib-2D had effects on grain number per spike (GNS) and thousand-grain weight (TGW). Comparison analysis suggested that QFll/Fla.cib-2B was likely a new locus. Two candidate genes, TraesCS2B03G0222800 and TraesCS2B03G0230000, associated with leaf development within the interval of QFll/Fla.cib-2B were identified based on expression-pattern analysis, gene annotation, ortholog analysis, and sequence variation. The major QTL and markers reported here provide valuable information for understanding the genetic mechanism underlying flag leaf size as well as breeding utilization in wheat.

Pinpointing genomic loci for drought-induced proline and hydrogen peroxide accumulation in bread wheat under field conditions

BACKGROUND

Proline (Pro) and hydrogen peroxide (H₂O₂) play a critical role in plants during drought adaptation. Genetic mapping for drought-induced Pro and H₂O₂ production under field conditions is very limited in crop plants since their phenotyping with large populations is labor-intensive. A genome-wide association study (GWAS) of a diversity panel comprised of 184 bread wheat cultivars grown in natural field (control) and rain-out shelter (drought) environments was performed to identify candidate loci and genes regulating Pro and H₂O₂ accumulation induced by drought.

RESULTS

The GWAS identified top significant marker-trait associations (MTAs) on 1A and 2A chromosomes, respectively for Pro and H₂O₂ in response to drought. Similarly, MTAs for stress tolerance index (STI) of Pro and H₂O₂ were identified on 5B and 1B chromosomes, respectively. Total 143 significant MTAs were identified including 36 and 71 were linked to drought and 2 and 34 were linked to STI for Pro and H₂O_2, respectively. Next, linkage disequilibrium analysis revealed minor alleles of significant single-markers and haplotypes were associated with higher Pro and H₂O₂ accumulation under drought. Several putative candidate genes for Pro and H₂O₂ content encode proteins with kinase, transporter or protein-binding activities.

CONCLUSIONS

The identified genetic factors associated with Pro and H₂O₂ biosynthesis underlying drought adaptation lay a fundamental basis for functional studies and future marker-assisted breeding programs.

An established protocol for generating transgenic wheat for wheat functional genomics via particle bombardment

Wheat is one of the most important food crops in the world and is considered one of the top targets in crop biotechnology. With the high-quality reference genomes of wheat and its relative species and the recent burst of genomic resources in Triticeae, demands to perform gene functional studies in wheat and genetic improvement have been rapidly increasing, requiring that production of transgenic wheat should become a routine technique. While established for more than 20 years, the particle bombardment-mediated wheat transformation has not become routine yet, with only a handful of labs being proficient in this technique. This could be due to, at least partly, the low transformation efficiency and the technical difficulties. Here, we describe the current version of this method through adaptation and optimization. We report the detailed protocol of producing transgenic wheat by the particle gun, including several critical steps, from the selection of appropriate explants (i.e., immature scutella), the preparation of DNA-coated gold particles, and several established strategies of tissue culture. More importantly, with over 20 years of experience in wheat transformation in our lab, we share the many technical details and recommendations and emphasize that the particle bombardment-mediated approach has fewer limitations in genotype dependency and vector construction when compared with the Agrobacterium-mediated methods. The particle bombardment-mediated method has been successful for over 30 wheat genotypes, from the tetraploid durum wheat to the hexaploid common wheat, from modern elite varieties to landraces. In conclusion, the particle bombardment-mediated wheat transformation has demonstrated its potential and wide applications, and the full set of protocol, experience, and successful reports in many wheat genotypes described here will further its impacts, making it a routine and robust technique in crop research labs worldwide.

Genome-wide association study reveals structural chromosome variations with phenotypic effects in wheat (Triticum aestivum L.)

Structural chromosome variations (SCVs) are large-scale genomic variations that can be detected by fluorescence in situ hybridization (FISH). SCVs have played important roles in the genome evolution of wheat (Triticum aestivum L.), but little is known about their genetic effects. In this study, a total of 543 wheat accessions from the Chinese wheat mini-core collection and the Shanxi Province wheat collection were used for chromosome analysis using oligonucleotide probe multiplex FISH. A total of 139 SCVs including translocations, pericentric inversions, presence/absence variations (PAVs), and copy number variations (CNVs) in heterochromatin were identified at 230 loci. The distribution frequency of SCVs varied between ecological regions and between landraces and modern cultivars. Structural analysis using SCVs as markers clearly divided the landraces and modern cultivars into different groups. There are very clear instances illustrating alien introgression and wide application of foreign germplasms improved the chromosome diversity of Chinese modern wheat cultivars. A genome-wide association study (GWAS) identified 29 SCVs associated with 12 phenotypic traits, and five (RT4AS•4AL-1DS/1DL•1DS-4AL, Mg2A-3, Mr3B-10, Mr7B-13, and Mr4A-7) of them were further validated using a doubled haploid population and advanced sib-lines, implying the potential value of these SCVs. Importantly, the number of favored SCVs that were associated with agronomic trait improvement was significantly higher in modern cultivars compared to landraces, indicating positive selection in wheat breeding. This study demonstrates the significant effects of SCVs during wheat breeding and provides an efficient method of mining favored SCVs in wheat and other crops.

Genome-Wide Identification and Analysis of FKBP Gene Family in Wheat (Triticum asetivum)

FK506-binding protein (FKBP) genes have been found to play vital roles in plant development and abiotic stress responses. However, limited information is available about this gene family in wheat (Triticum aestivum L.). In this study, a total of 64 FKBP genes were identified in wheat via a genome-wide analysis involving a homologous search of the latest wheat genome data, which was unevenly distributed in 21 chromosomes, encoded 152 to 649 amino acids with molecular weights ranging from 16 kDa to 72 kDa, and was localized in the chloroplast, cytoplasm, nucleus, mitochondria, peroxisome and endoplasmic reticulum. Based on sequence alignment and phylogenetic analysis, 64 TaFKBPs were divided into four different groups or subfamilies, providing evidence of an evolutionary relationship with Aegilops tauschii, Brachypodium distachyon, Triticum dicoccoides, Arabidopsis thaliana and Oryza sativa. Hormone-related, abiotic stress-related and development-related cis-elements were preferentially presented in promoters of TaFKBPs. The expression levels of TaFKBP genes were investigated using transcriptome data from the WheatExp database, which exhibited tissue-specific expression patterns. Moreover, TaFKBPs responded to drought and heat stress, and nine of them were randomly selected for validation by qRT-PCR. Yeast cells expressing TaFKBP19-2B-2 or TaFKBP18-6B showed increased influence on drought stress, indicating their negative roles in drought tolerance. Collectively, our results provide valuable information about the FKBP gene family in wheat and contribute to further characterization of FKBPs during plant development and abiotic stress responses, especially in drought stress.

Genome-Wide Analysis of SQUAMOSA-Promoter-Binding Protein-like Family in Flowering Pleioblastus pygmaeus

SQUAMOSA Promoter-Binding Protein-Like (SPL) family is well-known for playing an important role in plant growth and development, specifically in the reproductive process. Bamboo plants have special reproductive characteristics with a prolonged vegetative phase and uncertain flowering time. However, the underlying functions of SPL genes in reproductive growth are undisclosed in bamboo plants. In the study, a total of 28 SPLs were screened from an ornamental dwarf bamboo species, Pleioblastus pygmaeus. Phylogenetic analysis indicates that 183 SPLs from eight plant species can be classified into nine subfamilies, and the 28 PpSPLs are distributed among eight subfamilies. Homologous analysis shows that as many as 32 pairs of homologous genes were found between P. pygmaeus and rice, and 83 pairs were found between P. pygmaeus and Moso bamboo, whose Ka/Ks values are all <1. MiRNA target prediction reveals that 13 out of the 28 PpSPLs have recognition sites complementary to miRNA156. To screen the SPLs involved in the reproductive growth of bamboo plants, the mRNA abundance of the 28 PpSPLs was profiled in the different tissues of flowering P. pygmaeus and non-flowering plants by RNA-Seq. Moreover, the relative expression level of eight PpSPLs is significantly higher in flowering P. pygmaeus than that in non-flowering plants, which was also validated by RT-qPCR. Combined with phylogenetic analysis and homologous analysis, the eight significant, differentially expressed PpSPLs were identified to be associated with the reproductive process and flower organ development. Among them, there are four potential miRNA156-targeting PpSPLs involved in the flowering process. Of significant interest in the study is the identification of 28 SPLs and the exploration of four key flowering-related SPLs from P. pygmaeus, which provides a theoretic basis for revealing the underlying functions of SPLs in the reproductive growth of bamboo plants.

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

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

Alteration of synonymous codon usage bias accompanies polyploidization in wheat

The diploidization of polyploid genomes is accompanied by genomic variation, including synonymous nucleotide substitutions that may lead to synonymous codon usage bias (SCUB). SCUB can mirror the evolutionary specialization of plants, but its effect on the formation of polyploidies is not well documented. We explored this issue here with hexaploid wheat and its progenitors. Synonymous codons (SCs) ending in either cytosine (NNC) or guanidine (NNG) were more frequent than those ending in either adenosine (NNA) or thymine (NNT), and the preference for NNC/G codons followed the increase in genome ploidy. The ratios between NNC/G and NNA/T codons gradually decreased in genes with more introns, and the difference in these ratios between wheat and its progenitors diminished with increasing ploidy. SCUB frequencies were heterogeneous among exons, and the bias preferred to NNA/T in more internal exons, especially for genes with more exons; while the preference did not appear to associate with ploidy. The SCUB alteration of the progenitors was different during the formation of hexaploid wheat, so that SCUB was the homogeneous among A, B and D subgenomes. DNA methylation-mediated conversion from cytosine to thymine weakened following the increase of genome ploidy, coinciding with the stronger bias for NNC/G SCs in the genome as a function of ploidy, suggesting that SCUB contribute to the epigenetic variation in hexaploid wheat. The patterns in SCUB mirrored the formation of hexaploid wheat, which provides new insight into genome shock-induced genetic variation during polyploidization. SCs representing non-neutral synonymous mutations can be used for genetic dissection and improvement of agricultural traits of wheat and other polyploidies.

Histone acetyltransferase TaHAG1 interacts with TaPLATZ5 to activate TaPAD4 expression and positively contributes to powdery mildew resistance in wheat

Plants have evolved a two-branched innate immune system to detect and cope with pathogen attack, which are initiated by cell-surface and intracellular immune receptors leading to pattern-triggered immunity (PTI) and effector-triggered immunity (ETI), respectively. A core transducer including PAD4-EDS1 node is proposed as the convergence point for a two-tiered immune system in conferring pathogen immunity. However, the transcriptional regulatory mechanisms controlling expression of these key transducers remain largely unknown. Here, we identified histone acetyltransferase TaHAG1 as a positive regulator of powdery mildew resistance in wheat. TaHAG1 regulates expression of key transducer gene TaPAD4 and promotes SA and reactive oxygen species accumulation to accomplish resistance to Bgt infection. Moreover, overexpression and CRISPR-mediated knockout of TaPAD4 validate its role in wheat powdery mildew resistance. Furthermore, TaHAG1 physically interacts with TaPLATZ5, a plant-specific zinc-binding protein. TaPLATZ5 directly binds to promoter of TaPAD4 and together with TaHAG1 to potentiate the expression of TaPAD4 by increasing the levels of H3 acetylation. Our study revealed a key transcription regulatory node in which TaHAG1 acts as an epigenetic modulator and interacts with TaPLATZ5 that confers powdery mildew resistance in wheat through activating a convergence point gene between PTI and ETI, which could be effective for genetic improvement of disease resistance in wheat and other crops.

Small secreted peptides encoded on the wheat (triticum aestivum L.) genome and their potential roles in stress responses

Small secreted peptides (SSPs) are important signals for cell-to-cell communication in plant, involved in a variety of growth and developmental processes, as well as responses to stresses. While a large number of SSPs have been identified and characterized in various plant species, little is known about SSPs in wheat, one of the most important cereal crops. In this study, 4,981 putative SSPs were identified on the wheat genome, among which 1,790 TaSSPs were grouped into 38 known SSP families. The result also suggested that a large number of the putaitive wheat SSPs, Cys-rich peptides in particular, remained to be characterized. Several TaSSP genes were found to encode multiple SSP domains, including CLE, HEVEIN and HAIRPININ domains, and two potentially novel TaSSP family DYY and CRP8CI were identified manually among unpredicted TaSSPs. Analysis on the transcriptomic data showed that a great proportion of TaSSPs were expressed in response to abiotic stresses. Exogenous application of the TaCEPID peptide encoded by TraesCS1D02G130700 enhanced the tolerance of wheat plants to drought and salinity, suggesting porential roles of SSPs in regulating stress responses in wheat.

Wheat genomic study for genetic improvement of traits in China

Bread wheat (Triticum aestivum L.) is a major crop that feeds 40% of the world's population. Over the past several decades, advances in genomics have led to tremendous achievements in understanding the origin and domestication of wheat, and the genetic basis of agronomically important traits, which promote the breeding of elite varieties. In this review, we focus on progress that has been made in genomic research and genetic improvement of traits such as grain yield, end-use traits, flowering regulation, nutrient use efficiency, and biotic and abiotic stress responses, and various breeding strategies that contributed mainly by Chinese scientists. Functional genomic research in wheat is entering a new era with the availability of multiple reference wheat genome assemblies and the development of cutting-edge technologies such as precise genome editing tools, high-throughput phenotyping platforms, sequencing-based cloning strategies, high-efficiency genetic transformation systems, and speed-breeding facilities. These insights will further extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture in China and throughout the world.

Reduction of Rhizoctonia cerealis Infection on Wheat Through Host- and Spray-Induced Gene Silencing of an Orphan Secreted Gene

Rhizoctonia cerealis is a soilborne fungus that can cause sharp eyespot in wheat, resulting in massive yield losses found in many countries. Due to the lack of resistant cultivars, fungicides have been widely used to control this pathogen. However, chemical control is not environmentally friendly and is costly. Meanwhile, the lack of genetic transformation tools has hindered the functional characterization of virulence genes. In this study, we attempted to characterize the function of virulence genes by two transient methods, host-induced gene silencing (HIGS) and spray-induced gene silencing (SIGS), which use RNA interference to suppress the pathogenic development. We identified ten secretory orphan genes from the genome. After silencing these ten genes, only the RcOSP1 knocked-down plant significantly inhibited the growth of R. cerealis. We then described RcOSP1 as an effector that could impair wheat biological processes and suppress pathogen-associated molecular pattern-triggered immunity in the infection process. These findings confirm that HIGS and SIGS can be practical tools for researching R. cerealis virulence genes. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Fine mapping and characterization of a major QTL for grain weight on wheat chromosome arm 5DL

We fine mapped QTL QTKW.caas-5DL for thousand kernel weight in wheat, predicted candidate genes and developed a breeding-applicable marker. Thousand kernel weight (TKW) is an important yield component trait in wheat, and identification of the underlying genetic loci is helpful for yield improvement. We previously identified a stable quantitative trait locus (QTL) QTKW.caas-5DL for TKW in a Doumai/Shi4185 recombinant inbred line (RIL) population. Here we performed fine mapping of QTKW.caas-5DL using secondary populations derived from 15 heterozygous recombinants and delimited the QTL to an approximate 3.9 Mb physical interval from 409.9 to 413.8 Mb according to the Chinese Spring (CS) reference genome. Analysis of genomic synteny showed that annotated genes in the physical interval had high collinearity among CS and eight other wheat genomes. Seven genes with sequence variation and/or differential expression between parents were predicted as candidates for QTKW.caas-5DL based on whole-genome resequencing and transcriptome assays. A kompetitive allele-specific PCR (KASP) marker for QTKW.caas-5DL was developed, and genotyping confirmed a significant association with TKW but not with other yield component traits in a panel of elite wheat cultivars. The superior allele of QTKW.caas-5DL was frequent in a panel of cultivars, suggesting that it had undergone positive selection. These findings not only lay a foundation for map-based cloning of QTKW.caas-5DL but also provide an efficient tool for marker-assisted selection.

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

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

Identification of ankyrin-transmembrane-type subfamily genes in Triticeae species reveals TaANKTM2A-5 regulates powdery mildew resistance in wheat

The ankyrin-transmembrane (ANKTM) subfamily is the most abundant subgroup of the ANK superfamily, with critical roles in pathogen defense. However, the function of ANKTM proteins in wheat immunity remains largely unexplored. Here, a total of 381 ANKTMs were identified from five Triticeae species and Arabidopsis, constituting five classes. Among them, class a only contains proteins from Triticeae species and the number of ANKTM in class a of wheat is significantly larger than expected, even after consideration of the ploidy level. Tandem duplication analysis of ANKTM indicates that Triticum urartu, Triticum dicoccoides and wheat all had experienced tandem duplication events which in wheat-produced ANKTM genes all clustered in class a. The above suggests that not only did the genome polyploidization result in the increase of ANKTM gene number, but that tandem duplication is also a mechanism for the expansion of this subfamily. Micro-collinearity analysis of Triticeae ANKTMs indicates that some ANKTM type genes evolved into other types of ANKs in the evolution process. Public RNA-seq data showed that most of the genes in class d and class e are expressed, and some of them show differential responses to biotic stresses. Furthermore, qRT-PCR results showed that some ANKTMs in class d and class e responded to powdery mildew. Silencing of TaANKTM2A-5 by barley stripe mosaic virus-induced gene silencing compromised powdery mildew resistance in common wheat Bainongaikang58. Findings in this study not only help to understand the evolutionary process of ANKTM genes, but also form the basis for exploring disease resistance genes in the ANKTM gene family.

Utility of Triti-Map for bulk-segregated mapping of causal genes and regulatory elements in Triticeae

Triticeae species, including wheat, barley, and rye, are critical for global food security. Mapping agronomically important genes is crucial for elucidating molecular mechanisms and improving crops. However, Triticeae includes many wild relatives with desirable agronomic traits, and frequent introgressions occurred during Triticeae evolution and domestication. Thus, Triticeae genomes are generally large and complex, making the localization of genes or functional elements that control agronomic traits challenging. Here, we developed Triti-Map, which contains a suite of user-friendly computational packages specifically designed and optimized to overcome the obstacles of gene mapping in Triticeae, as well as a web interface integrating multi-omics data from Triticeae for the efficient mining of genes or functional elements that control particular traits. The Triti-Map pipeline accepts both DNA and RNA bulk-segregated sequencing data as well as traditional QTL data as inputs for locating genes and elucidating their functions. We illustrate the usage of Triti-Map with a combination of bulk-segregated ChIP-seq data to detect a wheat disease-resistance gene with its promoter sequence that is absent from the reference genome and clarify its evolutionary process. We hope that Triti-Map will facilitate gene isolation and accelerate Triticeae breeding.

The Calcium-Dependent Protein Kinase TaCDPK27 Positively Regulates Salt Tolerance in Wheat

As essential calcium ion (Ca2+) sensors in plants, calcium-dependent protein kinases (CDPKs) function in regulating the environmental adaptation of plants. However, the response mechanism of CDPKs to salt stress is not well understood. In the current study, the wheat salt-responsive gene TaCDPK27 was identified. The open reading frame (ORF) of TaCDPK27 was 1875 bp, coding 624 amino acids. The predicted molecular weight and isoelectric point were 68.905 kDa and 5.6, respectively. TaCDPK27 has the closest relationship with subgroup III members of the CDPK family of rice. Increased expression of TaCDPK27 in wheat seedling roots and leaves was triggered by 150 mM NaCl treatment. TaCDPK27 was mainly located in the cytoplasm. After NaCl treatment, some of this protein was transferred to the membrane. The inhibitory effect of TaCDPK27 silencing on the growth of wheat seedlings was slight. After exposure to 150 mM NaCl for 6 days, the NaCl stress tolerance of TaCDPK27-silenced wheat seedlings was reduced, with shorter lengths of both roots and leaves compared with those of the control seedlings. Moreover, silencing of TaCDPK27 further promoted the generation of reactive oxygen species (ROS); reduced the activities of superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT); aggravated the injury to photosystem II (PS II); and increased programmed cell death (PCD) in wheat leaves under NaCl treatment, confirming that the TaCDPK27-silenced seedlings exhibited more NaCl injury than control seedlings. Taken together, the decrease in NaCl tolerance in TaCDPK27-silenced seedlings was due to excessive ROS accumulation and subsequent aggravation of the NaCl-induced PCD. TaCDPK27 may be essential for positively regulating salt tolerance in wheat seedlings.

The Metacaspase TaMCA-Id Negatively Regulates Salt-Induced Programmed Cell Death and Functionally Links With Autophagy in Wheat

Metacaspases (MCAs), a family of caspase-like proteins, are important regulators of programmed cell death (PCD) in plant defense response. Autophagy is an important regulator of PCD. This study explored the underlying mechanism of the interaction among PCD, MCAs, and autophagy and their impact on wheat response to salt stress. In this study, the wheat salt-responsive gene TaMCA-Id was identified. The open reading frame (ORF) of TaMCA-Id was 1,071 bp, coding 356 amino acids. The predicted molecular weight and isoelectric point were 38,337.03 Da and 8.45, respectively. TaMCA-Id had classic characteristics of type I MCAs domains, a typical N-terminal pro-domain rich in proline. TaMCA-Id was mainly localized in the chloroplast and exhibited nucleocytoplasmictrafficking under NaCl treatment. Increased expression of TaMCA-Id in wheat seedling roots and leaves was triggered by 150 mM NaCl treatment. Silencing of TaMCA-Id enhanced sensitivity of wheat seedlings to NaCl stress. Under NaCl stress, TaMCA-Id-silenced seedlings exhibited a reduction in activities of superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT), higher accumulation of H₂O₂ and O 2 . - , more serious injury to photosystem II (PSII), increase in PCD level, and autophagy activity in leaves of wheat seedlings. These results indicated that TaMCA-Id functioned in PCD through interacting with autophagy under NaCl stress, which could be used to improve the salt tolerance of crop plants.

Genome-Wide Analysis of the Peptidase M24 Superfamily in Triticum aestivum Demonstrates That TaM24-9 Is Involved in Abiotic Stress Response

The peptidase M24 (Metallopeptidase 24, M24) superfamily is essential for plant growth, stress response, and pathogen defense. At present, there are few systematic reports on the identification and classification of members of the peptidase M24 proteins superfamily in wheat. In this work, we identified 53 putative candidate TaM24 genes. According to the protein sequences characteristics, these members can be roughly divided into three subfamilies: I, II, III. Most TaM24 genes are complex with multiple exons, and the motifs are relatively conserved in each sub-group. Through chromosome mapping analysis, we found that the 53 genes were unevenly distributed on 19 wheat chromosomes (except 3A and 3D), of which 68% were in triads. Analysis of gene duplication events showed that 62% of TaM24 genes in wheat came from fragment duplication events, and there were no tandem duplication events to amplify genes. Analysis of the promoter sequences of TaM24 genes revealed that cis-acting elements were rich in response elements to drought, osmotic stress, ABA, and MeJA. We also studied the expression of TaM24 in wheat tissues at developmental stages and abiotic stress. Then we selected TaM24-9 as the target for further analysis. The results showed that TaM24-9 genes strengthened the drought and salt tolerance of plants. Overall, our analysis showed that members of the peptidase M24 genes may participate in the abiotic stress response and provided potential gene resources for improving wheat resistance.

Identification of a novel major QTL from Chinese wheat cultivar Ji5265 for Fusarium head blight resistance in greenhouse

A novel major QTL for FHB resistance was mapped to a 6.8 Mb region on chromosome 2D in a Chinese wheat cultivar Ji5265, and diagnostic KASP markers were developed for detecting it in a worldwide wheat collection. Fusarium head blight (FHB) is a serious disease in wheat (Triticum aestivum L.) and causes significant reductions in grain yield and quality worldwide. Breeding for FHB resistance is the most effective strategy to minimize the losses caused by FHB; therefore, identification of major quantitative trait loci (QTLs) conferring FHB resistance and development of diagnostic markers for the QTLs are prerequisites for marker-assisted selection (MAS). Ji5265 is a Chinese wheat cultivar resistant to FHB in multiple environments. An F₆ population of 179 recombinant inbred lines (RILs) was developed from Ji5265 × Wheaton. The population was genotyped by genotyping-by-sequencing (GBS) and phenotyped for FHB Type II resistance in greenhouses. A major QTL, designated as QFhb-2DL, was mapped in a 6.8 Mb region between the markers GBS10238 and GBS12056 on the long arm of chromosome 2D in Ji5265 and explained ~ 30% of the phenotypic variation for FHB resistance. The effect of QFhb-2DL on FHB resistance was validated using near-isogenic lines (NILs) derived from residual heterozygotes from an F₆ RIL of Ji5265 × Wheaton. The two flanking markers were converted into Kompetitive allele-specific PCR (KASP) markers (KASP10238 and KASP12056) and validated to be diagnostic in a collection of 2,065 wheat accessions. These results indicate that QFhb-2DL is a novel major QTL for resistance to FHB spread within a spike (Type II) and the two KASP markers can be used for MAS to improve wheat FHB resistance in wheat breeding programs.

High-resolution detection of quantitative trait loci for seven important yield-related traits in wheat (Triticum aestivum L.) using a high-density SLAF-seq genetic map

BACKGROUND

Yield-related traits including thousand grain weight (TGW), grain number per spike (GNS), grain width (GW), grain length (GL), plant height (PH), spike length (SL), and spikelet number per spike (SNS) are greatly associated with grain yield of wheat (Triticum aestivum L.). To detect quantitative trait loci (QTL) associated with them, 193 recombinant inbred lines derived from two elite winter wheat varieties Chuanmai42 and Chuanmai39 were employed to perform QTL mapping in six/eight environments.

RESULTS

A total of 30 QTLs on chromosomes 1A, 1B, 1D, 2A, 2B, 2D, 3A, 4A, 5A, 5B, 6A, 6D, 7A, 7B and 7D were identified. Among them, six major QTLs QTgw.cib-6A.1, QTgw.cib-6A.2, QGw.cib-6A, QGl.cib-3A, QGl.cib-6A, and QSl.cib-2D explaining 5.96-23.75% of the phenotypic variance were detected in multi-environments and showed strong and stable effects on corresponding traits. Three QTL clusters on chromosomes 2D and 6A containing 10 QTLs were also detected, which showed significant pleiotropic effects on multiple traits. Additionally, three Kompetitive Allele Specific PCR (KASP) markers linked with five of these major QTLs were developed. Candidate genes of QTgw.cib-6A.1/QGl.cib-6A and QGl.cib-3A were analyzed based on the spatiotemporal expression patterns, gene annotation, and orthologous search.

CONCLUSIONS

Six major QTLs for TGW, GL, GW and SL were detected. Three KASP markers linked with five of these major QTLs were developed. These QTLs and KASP markers will be useful for elucidating the genetic architecture of grain yield and developing new wheat varieties with high and stable yield in wheat.

The fine mapping of dwarf gene Rht5 in bread wheat and its effects on plant height and main agronomic traits

Rht5 was narrowed to an approximately 1 Mb interval and had pleiotropic effects on plant height, spike length and grain size. TraesCS3B02G025600 was predicted as the possible candidate gene. Plant height is an important component related to plant architecture, lodging resistance, and yield performance. The utilization of dwarf genes has made great contributions to wheat breeding and production. In this study, two F₂ populations derived from the crosses of Jinmai47 and Ningchun45 with Marfed M were employed to identify the genetic region of reduce plant height 5 (Rht5), and their derived lines were used to evaluate its effects on plant height and main agronomic traits. Rht5 was fine-mapped between markers Kasp-25 and Kasp-23, in approximately 1 Mb region on chromosome 3BS, which harbored 17 high-confidence annotated genes based on the reference genome of Chinese Spring (IWGSC RefSeq v1.1). TraesCS3B02G025600 were predicted as the possible candidate gene based on its differential expression and sequence variation between dwarf and tall lines and parents. The results of phenotypic evaluation showed that Rht5 had pleiotropic effects on plant height, spike length, culm diameter, grain size and grain yield. The plant height of Rht5 dwarf lines was reduced by an average of 32.67% (32.53 cm) and 27.84% (33.62 cm) in the Jinmai47 and Ningchun45 population, respectively. While Rht5 showed significant and negative pleiotropic effects on culm diameter, aboveground biomass, grain yield, spike length, spikelet number, grain number per spike, grain size, grain weight and filling degree of basal second internode. The culm lodging resistance index (CLRI) of dwarf lines was significantly higher than that of tall lines in the two population. In conclusion, these results lay a foundation for understanding the dwarfing mechanism of Rht5.

Comprehensive Genome-Wide Identification, Characterization, and Expression Analysis of CCHC-Type Zinc Finger Gene Family in Wheat (Triticum aestivum L.)

The CCHC-type zinc finger proteins (CCHC-ZFPs) play versatile roles in plant growth, development and adaptation to the environment. However, little is known about functions of CCHC-ZFP gene family memebers in Triticum aestivum. In the present study, we identified a total of 50 TaCCHC-ZFP genes from the 21 wheat chromosomes, which were phylogenetically classified into eight groups based on their specific motifs and gene structures. The 43 segmentally duplicated TaCCHC-ZFP genes were retrieved, which formed 36 segmental duplication gene pairs. The collinearity analyses among wheat and other eight mono/dicots revealed that no gene pairs were found between wheat and the three dicots. The promoter analyses of the TaCCHC-ZFP genes showed that 636 environmental stress-responsive and phytohormone-responsive cis-elements. The gene ontology enrichment analysis indicated that all the TaCCHC-ZFP genes were annotated under nucleic acid binding and metal ion binding. A total of 91 MicroRNA (miRNA) binding sites were identified in 34 TaCCHC-ZFP genes according to the miRNA target analysis. Based on the public transcriptome data, the 38 TaCCHC-ZFP genes were identified as differentially expressed gene. The expression profiles of 15 TaCCHC-ZFP genes were verified by the quantitative real-time PCR assays, and the results showed that these genes were responsive to drought or heat treatments. Our work systematically investigated the gene structures, evolutionary features, and potential functions of TaCCHC-ZFP genes. It lays a foundation for further research and application of TaCCHC-ZFP genes in genetic improvement of T. aestivum.

Exome Sequencing from Bulked Segregant Analysis Identifies a Gene for All-Stage Resistance to Stripe Rust on Chromosome 1AL in Chinese Wheat Landrace 'Xiaohemai'

Stripe rust caused by Puccinia striiformis f. sp. tritici is one of the most destructive diseases of wheat. Identifying novel resistance genes applicable for developing disease-resistant cultivars is important for the sustainable control of wheat stripe rust. Chinese wheat landrace 'Xiaohemai' ('XHM') is an elite germplasm line with all-stage resistance (ASR) effective against predominant Chinese P. striiformis f. sp. tritici races. In this study, we performed a bulked segregant analysis coupled with exome capture sequencing (BSE-seq) to identify a candidate genomic region strongly associated with stripe rust resistance on chromosome 1AL in 173 F_2:3 lines derived from the cross 'XHM' × 'Avocet S'. The gene, designated as YrXH-1AL, was validated by a conventional quantitative trait locus analysis using newly developed Kompetitive allele-specific PCR (KASP) markers, explaining up to 48.50% of the phenotypic variance. By testing a secondary mapping population comprising 144 lines from the same cross at the seedling stage with prevalent P. striiformis f. sp. tritici race CYR34, YrXH-1AL was identified as a single Mendelian factor in a 1.5-cM interval flanked by KASP markers KP1A_484.33 and KP1A_490.09. This region corresponded to a 5.76-Mb genomic interval on 'Chinese Spring' chromosome 1AL. Furthermore, two cosegregating KASP markers showed high polymorphisms among 130 Chinese wheat cultivars and could be used for marker-assisted selection. Because no other Yr genes for ASR that originated from common wheat have been detected on chromosome 1AL, YrXH-1AL is likely a novel gene that can be incorporated into modern breeding materials to develop wheat cultivars with enhanced stripe rust resistance.

Genome sequences of five Sitopsis species of Aegilops and the origin of polyploid wheat B subgenome

Common wheat (Triticum aestivum, BBAADD) is a major staple food crop worldwide. The diploid progenitors of the A and D subgenomes have been unequivocally identified; that of B, however, remains ambiguous and controversial but is suspected to be related to species of Aegilops, section Sitopsis. Here, we report the assembly of chromosome-level genome sequences of all five Sitopsis species, namely Aegilops bicornis, Ae. longissima, Ae. searsii, Ae. sharonensis, and Ae. speltoides, as well as the partial assembly of the Amblyopyrum muticum (synonym Aegilops mutica) genome for phylogenetic analysis. Our results reveal that the donor of the common wheat B subgenome is a distinct, and most probably extinct, diploid species that diverged from an ancestral progenitor of the B lineage to which the still extant Ae. speltoides and Am. muticum belong. In addition, we identified interspecific genetic introgressions throughout the evolution of the Triticum/Aegilops species complex. The five Sitopsis species have various assembled genome sizes (4.11-5.89 Gb) with high proportions of repetitive sequences (85.99%-89.81%); nonetheless, they retain high collinearity with other genomes or subgenomes of species in the Triticum/Aegilops complex. Differences in genome size were primarily due to independent post-speciation amplification of transposons. We also identified a set of Sitopsis genes pertinent to important agronomic traits that can be harnessed for wheat breeding. These newly assembled genome resources provide a new roadmap for evolutionary and genetic studies of the Triticum/Aegilops complex, as well as for wheat improvement.

Genome-Wide Association Study Reveals the Genetic Basis of Five Quality Traits in Chinese Wheat

Bread wheat is a highly adaptable food crop grown extensively around the world and its quality genetic improvement has received wide attention. In this study, the genetic loci associated with five quality traits including protein content (PC), gluten content (GC), baking value (BV), grain hardness (HA), and sedimentation value (SV) in a population of 253 Chinese wheat grown in Inner Mongolia were investigated through genome wide association mapping. A total of 103 QTL containing 556 SNPs were significantly related to the five quality traits based on the phenotypic data collected from three environments and BLUP data. Of these QTL, 32 QTL were continuously detected under at least two experiments. Some QTL such as qBV3D.2/qHA3D.2 on 3D, qPC5A.3/qGC5A on 5A, qBV5D/qHA5D on 5D, qBV6B.2/qHA6B.3 on 6B, and qBV6D/qHA6D.1 on 6D were associated with multiple traits. In addition, distribution of favorable alleles of the stable QTL in the association panel and their effects on five quality traits were validated. Analysis of existing transcriptome data revealed that 34 genes were specifically highly expressed in grains during reproductive growth stages. The functions of these genes will be characterized in future experiments. This study provides novel insights into the genetic basis of quality traits in wheat.

Fine mapping of QPm.caas-3BS, a stable QTL for adult-plant resistance to powdery mildew in wheat (Triticum aestivum L.)

A stable QTL QPm.caas-3BS for adult-plant resistance to powdery mildew was mapped in an interval of 431 kb, and candidate genes were predicted based on gene sequences and expression profiles. Powdery mildew is a devastating foliar disease occurring in most wheat-growing areas. Characterization and fine mapping of genes for powdery mildew resistance can benefit marker-assisted breeding. We previously identified a stable quantitative trait locus (QTL) QPm.caas-3BS for adult-plant resistance to powdery mildew in a recombinant inbred line population of Zhou8425B/Chinese Spring by phenotyping across four environments. Using 11 heterozygous recombinants and high-density molecular markers, QPm.caas-3BS was delimited in a physical interval of approximately 3.91 Mb. Based on re-sequenced data and expression profiles, three genes TraesCS3B02G014800, TraesCS3B02G016800 and TraesCS3B02G019900 were associated with the powdery mildew resistance locus. Three gene-specific kompetitive allele-specific PCR (KASP) markers were developed from these genes and validated in the Zhou8425B derivatives and Zhou8425B/Chinese Spring population in which the resistance gene was mapped to a 0.3 cM interval flanked by KASP14800 and snp_50465, corresponding to a 431 kb region at the distal end of chromosome 3BS. Within the interval, TraesCS3B02G014800 was the most likely candidate gene for QPm.caas-3BS, but TraesCS3B02G016300 and TraesCS3B02G016400 were less likely candidates based on gene annotations and sequence variation between the parents. These results not only offer high-throughput KASP markers for improvement of powdery mildew resistance but also pave the way to map-based cloning of the resistance gene.

Genome-Wide Analysis of HSP70s in Hexaploid Wheat: Tandem Duplication, Heat Response, and Regulation

HSP70s play crucial roles in plant growth and development, as well as in stress response. Knowledge of the distribution and heat response of HSP70s is important to understand heat adaptation and facilitate thermotolerance improvement in wheat. In this study, we comprehensively analyzed the distribution of HSP70s in hexaploid wheat (TaHSP70s) and its relatives, and we found an obvious expansion of TaHSP70s in the D genome of hexaploid wheat. Meanwhile, a large portion of tandem duplication events occurred in hexaploid wheat. Among the 84 identified TaHSP70s, more than 64% were present as homeologs. The expression profiles of TaHSP70s in triads tended to be expressed more in non-stressful and heat stress conditions. Intriguingly, many TaHSP70s were especially heat responsive. Tandem duplicated TaHSP70s also participated in heat response and growth development. Further HSE analysis revealed divergent distribution of HSEs in the promoter regions of TaHSP70 homeologs, which suggested a distinct heat regulatory mechanism. Our results indicated that the heat response of TaHSP70s may experience a different regulation, and this regulation, together with the expression of tandem duplicated TaHSP70s, may help hexaploid wheat to adapt to heat conditions.

Genome-Wide Analysis of Serine Hydroxymethyltransferase Genes in Triticeae Species Reveals That TaSHMT3A-1 Regulates Fusarium Head Blight Resistance in Wheat

Serine hydroxymethyltransferase (SHMT) plays a pivotal role in cellular one-carbon, photorespiration pathways and it influences the resistance to biotic and abiotic stresses. However, the function of SHMT proteins in wheat remains largely unexplored. In the present study, SHMT genes in five Triticeae species, Oryza sativa, and four dicotyledon species were identified based on whole genome information. The origin history of the target gene was traced by micro-collinearity analysis. Gene expression patterns of TaSHMTs in different tissues, various biotic stresses, exogenous hormones, and two biotic stresses were determined by Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR). The function of the selected TaSHMT3A-1 was studied by barley stripe mosaic virus-induced gene silencing in common wheat Bainong207. A total of 64 SHMT members were identified and further classified into two main classes based on the structure of SHMT proteins. The gene structure and motif composition analyses revealed that SHMTs kept relatively conserved within the same subclasses. Interestingly, there was a gene, TdSHMT7B-1, on chromosome 7B of Triticum dicoccoides, but there was no SHMT gene on chromosome 7 of other analyzed Triticeae species; TdSHMT7B-1 had fewer exons and conserved motifs than the genes in the same subclass, suggesting that the gene of TdSHMT7B-1 has a notable evolutionary progress. The micro-collinearity relationship showed that no homologs of TaSHMT3A-1 and its two neighboring genes were found in the collinearity region of Triticum urartu, and there were 27 genes inserted into the collinearity region of T. urartu. Furthermore, qRT-PCR results showed that TaSHMT3A-1 was responsive to abiotic stresses (NaCl and cold), abscisic acid, methyl jasmonate, and hydrogen peroxide. Significantly, upon Fusarium graminearum infection, the expression of TaSHMT3A-1 was highly upregulated in resistant cultivar Sumai3. More importantly, silencing of TaSHMT3A-1 compromises Fusarium head blight resistance in common wheat Bainong207. Our new findings suggest that the TaSHMT3A-1 gene in wheat plays an important role in resistance to Fusarium head blight. This provides a valuable reference for further study on the function of this gene family.