WheatExp: an RNA-seq expression database for polyploid wheat

Functional divergences of natural variations of TaNAM-A1 in controlling leaf senescence during wheat grain filling

Leaf senescence is an essential physiological process related to grain yield potential and nutritional quality. Green leaf duration (GLD) after anthesis directly reflects the leaf senescence process and exhibits large genotypic differences in common wheat; however, the underlying gene regulatory mechanism is still lacking. Here, we identified TaNAM-A1 as the causal gene of the major loci qGLD-6A for GLD during grain filling by map-based cloning. Transgenic assays and TILLING mutant analyses demonstrated that TaNAM-A1 played a critical role in regulating leaf senescence, and also affected spike length and grain size. Furthermore, the functional divergences among the three haplotypes of TaNAM-A1 were systematically evaluated. Wheat varieties with TaNAM-A1d (containing two mutations in the coding DNA sequence of TaNAM-A1) exhibited a longer GLD and superior yield-related traits compared to those with the wild type TaNAM-A1a. All three haplotypes were functional in activating the expression of genes involved in macromolecule degradation and mineral nutrient remobilization, with TaNAM-A1a showing the strongest activity and TaNAM-A1d the weakest. TaNAM-A1 also modulated the expression of the senescence-related transcription factors TaNAC-S-7A and TaNAC016-3A. TaNAC016-3A enhanced the transcriptional activation ability of TaNAM-A1a by protein-protein interaction, thereby promoting the senescence process. Our study offers new insights into the fine-tuning of the leaf functional period and grain yield formation for wheat breeding under various geographical climatic conditions.

Bulked segregant RNA-seq reveals complex resistance expression profile to powdery mildew in wild emmer wheat W762

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is one of the most destructive fungal diseases threatening global wheat production. Exploring powdery mildew resistance (Pm) gene(s) and dissecting the molecular mechanism of the host resistance are critical to effectively and reasonably control this disease. Durum wheat (Triticum turgidum L. var. durumDesf.) is an important gene donor for wheat improvement against powdery mildew. In this study, a resistant durum wheat accession W762 was used to investigate its potential resistance component(s) and profile its expression pattern in responding to Bgt invasion using bulked segregant RNA-Seq (BSR-Seq) and further qRT-PCR verification. Genetic analysis showed that the powdery mildew resistance in W762 did not meet monogenic inheritance and complex genetic model might exist within the population of W762 × Langdon (susceptible durum wheat). After BSR-Seq, 6,196 consistently different single nucleotide polymorphisms (SNPs) were called between resistant and susceptible parents and bulks, and among them, 763 SNPs were assigned to the chromosome arm 7B. Subsequently, 3,653 differentially expressed genes (DEGs) between resistant and susceptible parents and bulks were annotated and analyzed by Gene Ontology (GO), Cluster of Orthologous Groups (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment. The potential regulated genes were selected and analyzed their temporal expression patterns following Bgt inoculation. As a result, nine disease-related genes showed distinctive expression profile after Bgt invasion and might serve as potential targets to regulate the resistance against powdery mildew in W762. Our study could lay a foundation for analysis of the molecular mechanism and also provide potential targets for the improvement of durable resistance against powdery mildew.

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.

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.

Novel genomic regions on chromosome 5B controlling wheat powdery mildew seedling resistance under Egyptian conditions

Wheat powdery mildew (PM) causes significant yield losses worldwide. None of the Egyptian wheat cultivars was detected to be highly resistant to such a severe disease. Therefore, a diverse spring wheat panel was evaluated for PM seedling resistance using different Bgt conidiospores collected from Egyptian fields in two growing seasons. The evaluation was done in two separate experiments. Highly significant differences were found between the two experiments suggesting the presence of different isolates populations. Highly significant differences were found among the tested genotypes confirming the ability to improve PM resistance using the recent panel. Genome-wide association study (GWAS) was done for each experiment separately and a total of 71 significant markers located within 36 gene models were identified. The majority of these markers are located on chromosome 5B. Haplotype block analysis identified seven blocks containing the significant markers on chromosome 5B. Five gene models were identified on the short arm of the chromosome. Gene enrichment analysis identified five and seven pathways based on the biological process and molecular functions respectively for the detected gene models. All these pathways are associated with disease resistance in wheat. The genomic regions on 5B seem to be novel regions that are associated with PM resistance under Egyptian conditions. Selection of superior genotypes was done and Grecian genotypes seem to be a good source for improving PM resistance under Egyptian conditions.

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.

BarleyExpDB: an integrative gene expression database for barley

BACKGROUND

RNA-sequencing (RNA-seq) has been widely used to study the dynamic expression patterns of transcribed genes, which can lead to new biological insights. However, processing and analyzing these huge amounts of histological data remains a great challenge for wet labs and field researchers who lack bioinformatics experience and computational resources.

RESULTS

We present BarleyExpDB, an easy-to-operate, free, and web-accessible database that integrates transcriptional profiles of barley at different growth and developmental stages, tissues, and stress conditions, as well as differential expression of mutants and populations to build a platform for barley expression and visualization. The expression of a gene of interest can be easily queried by searching by known gene ID or sequence similarity. Expression data can be displayed as a heat map, along with functional descriptions as well as Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Proteins Families Database, and Simple Modular Architecture Research Tool annotations.

CONCLUSIONS

BarleyExpDB will serve as a valuable resource for the barley research community to leverage the vast publicly available RNA-seq datasets for functional genomics research and crop molecular breeding.

Machine Learning-Assisted Approaches in Modernized Plant Breeding Programs

In the face of a growing global population, plant breeding is being used as a sustainable tool for increasing food security. A wide range of high-throughput omics technologies have been developed and used in plant breeding to accelerate crop improvement and develop new varieties with higher yield performance and greater resilience to climate changes, pests, and diseases. With the use of these new advanced technologies, large amounts of data have been generated on the genetic architecture of plants, which can be exploited for manipulating the key characteristics of plants that are important for crop improvement. Therefore, plant breeders have relied on high-performance computing, bioinformatics tools, and artificial intelligence (AI), such as machine-learning (ML) methods, to efficiently analyze this vast amount of complex data. The use of bigdata coupled with ML in plant breeding has the potential to revolutionize the field and increase food security. In this review, some of the challenges of this method along with some of the opportunities it can create will be discussed. In particular, we provide information about the basis of bigdata, AI, ML, and their related sub-groups. In addition, the bases and functions of some learning algorithms that are commonly used in plant breeding, three common data integration strategies for the better integration of different breeding datasets using appropriate learning algorithms, and future prospects for the application of novel algorithms in plant breeding will be discussed. The use of ML algorithms in plant breeding will equip breeders with efficient and effective tools to accelerate the development of new plant varieties and improve the efficiency of the breeding process, which are important for tackling some of the challenges facing agriculture in the era of climate change.

Chromosome groups 5, 6 and 7 harbor major quantitative trait loci controlling root traits in bread wheat (Triticum aestivum L.)

Identifying genomic regions for root traits in bread wheat can help breeders develop climate-resilient and high-yielding wheat varieties with desirable root traits. This study used the recombinant inbred line (RIL) population of Synthetic W7984 × Opata M85 to identify quantitative trait loci (QTL) for different root traits such as rooting depth (RD), root dry mass (RM), total root length (RL), root diameter (Rdia) and root surface areas (RSA1 for coarse roots and RSA2 for fine roots) under controlled conditions in a semi-hydroponic system. We detected 14 QTL for eight root traits on nine wheat chromosomes; we discovered three QTL each for RD and RSA1, two QTL each for RM and RSA2, and one QTL each for RL, Rdia, specific root length and nodal root number per plant. The detected QTL were concentrated on chromosome groups 5, 6 and 7. The QTL for shallow RD (Q.rd.uwa.7BL: Xbarc50) and high RM (Q.rm.uwa.6AS: Xgwm334) were validated in two independent F₂ populations of Synthetic W7984 × Chara and Opata M85 × Cascade, respectively. Genotypes containing negative alleles for Q.rd.uwa.7BL had 52% shallower RD than other Synthetic W7984 × Chara population lines. Genotypes with the positive alleles for Q.rm.uwa.6AS had 31.58% higher RM than other Opata M85 × Cascade population lines. Further, we identified 21 putative candidate genes for RD (Q.rd.uwa.7BL) and 13 for RM (Q.rm.uwa.6AS); TraesCS6A01G020400, TraesCS6A01G024400 and TraesCS6A01G021000 identified from Q.rm.uwa.6AS, and TraesCS7B01G404000, TraesCS7B01G254900 and TraesCS7B01G446200 identified from Q.rd.uwa.7BL encoded important proteins for root traits. We found germin-like protein encoding genes in both Q.rd.uwa.7BL and Q.rm.uwa.6AS regions. These genes may play an important role in RM and RD improvement. The identified QTL, especially the validated QTL and putative candidate genes are valuable genetic resources for future root 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.

CRISPR/Cas9 editing of wheat Ppd-1 gene homoeologs alters spike architecture and grain morphometric traits

Mutations in Photoperiod-1 (Ppd-1) gene are known to modify flowering time and yield in wheat. We cloned TaPpd-1 from wheat and found high similarity among the three homoeologs of TaPpd-1. To clarify the characteristics of TaPpd-1 homoeologs in different photoperiod conditions for inflorescence architecture and yield, we used CRISPR/Cas9 system to generate Tappd-1 mutant plants by simultaneous modification of the three homoeologs of wheat Ppd-1. Tappd-1 mutant plants showed no off-target mutations. Four T₀-edited lines under short-day length and three lines under long-day length conditions with the mutation frequency of 25% and 21%, respectively. These putative transgenic plants of all the lines were self-fertilized and generated T₁ and T₂ progenies and were evaluated by phenotypic and expression analysis. Results demonstrated that simultaneously edited TaPpd-1- A1, B1, and D1 homoeologs gene copies in T2_SDL-8-4, T2_SDL-4-5, T2_SDL-3-9, and T2_{_}LDL-10-9 showed similar spike inflorescence, flowering time, and significantly increase in 1000-grain weight, grain area, grain width, grain length, plant height, and spikelets per spike due to mutation in both alleles of Ppd-B1 and Ppd-D1 homoeologs but only spike length was decreased in T2_SDL-8-4, T2_SDL-4-5, and T2_{_}LDL-13-3 mutant lines due to mutation in both alleles of Ppd-A1 homoeolog under both conditions. Our results indicate that all TaPpd1 gene homoeologs influence wheat spike development by affecting both late flowering and earlier flowering but single mutant TaPpd-A1 homoeolog affect lowest as compared to the combination with double mutants of TaPpd-B1 and TaPpd-D1, TaPpd-A1 and TaPpd-B1, and TaPpd-A1 and TaPpd-D1 homoeologs for yield enhancement. Our findings further raised the idea that the relative expression of the various genomic copies of TaPpd-1 homoeologs may have an impact on the spike inflorescence architecture and grain morphometric features in wheat cultivars.

In-Depth Understanding of the Genetic Control of Stripe Rust Resistance (Puccinia striiformis f. sp. tritici) Induced in Wheat (Triticum aestivum) by Trichoderma asperellum T34

Wheat stripe rust (caused by Puccinia striiformis f. tritici Erikss.) causes severe yield losses worldwide. Due to the continuous appearance of new stripe rust races, resistance has been broken in most of the highly resistant genotypes in Egypt and worldwide. Therefore, looking for new ways to resist such a severe disease is urgently needed. Trichoderma asperellum strain T34 has been known as an effective bioagent against many crop diseases. It exists naturally in Egyptian fields. Therefore, in our study, the effectiveness of strain T34 was tested as a bioagent against wheat stripe rust. For this purpose, 198 spring wheat genotypes were tested for their resistance against two different P. striiformis f. tritici populations collected from the Egyptian fields. The most highly aggressive P. striiformis f. tritici population was used to test the effectiveness of strain T34. Highly significant differences were found between strain T34 and stripe rust, suggesting the effectiveness of strain T34 in stripe rust resistance. A genome-wide association study identified 48 gene models controlling resistance under normal conditions and 46 gene models controlling strain T34-induced resistance. Of these gene models, only one common gene model was found, suggesting the presence of two different genetic systems controlling resistance under each condition. The pathways of the biological processes were investigated under both conditions. This study provided in-depth understanding of genetic control and, hence, will accelerate the future of wheat breeding programs for stripe rust resistance.

Genome-wide characterization of FK506-binding proteins, parvulins and phospho-tyrosyl phosphatase activators in wheat and their regulation by heat stress

Peptidyl-prolyl cis-trans isomerases (PPIases) are ubiquitous proteins which are essential for cis-trans isomerisation of peptide bonds preceding the proline residue. PPIases are categorized into four sub-families viz., cyclophilins, FK506-binding proteins (FKBPs), parvulins and protein phosphatase 2A phosphatase activators (PTPAs). Apart from catalysing the cis-trans isomerization, these proteins have also been implicated in diverse cellular functions. Though PPIases have been identified in several important crop plants, information on these proteins, except cyclophilins, is scanty in wheat. In order to understand the role of these genes in wheat, we carried out genome-wide identification using computational approaches. The present study resulted in identification of 71 FKBP (TaFKBP) 12 parvulin (TaPar) and 3 PTPA (TaPTPA) genes in hexaploid wheat genome, which are distributed on different chromosomes with uneven gene densities. The TaFKBP and TaPar proteins, besides PPIase domain, also contain additional domains, indicating functional diversification. In silico prediction also revealed that TaFKBPs are localized to ER, nucleus, chloroplast and cytoplasm, while the TaPars are confined to cytoplasm and nucleus. The TaPTPAs, on the contrary, appear to be present only in the cytoplasm. Evolutionary studies predicted that most of the TaFKBP, TaPar and TaPTPA genes in hexaploid wheat have been derived from their progenitor species, with some events of loss or gain. Syntenic analysis revealed the presence of many collinear blocks of TaFKBP genes in wheat and its sub-genome donors. qRT-PCR analysis demonstrated that expression of TaFKBP and TaPar genes is regulated differentially by heat stress, suggesting their likely involvement in thermotolerance. The findings of this study will provide basis for further functional characterization of these genes and their likely applications in crop improvement.

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.

Transcriptome and Proteome Co-Profiling Offers an Understanding of Pre-Harvest Sprouting (PHS) Molecular Mechanisms in Wheat (Triticum aestivum)

While wheat (Triticum aestivum L.) is a widely grown and enjoyed crop, the diverse and complex global situation and climate are exacerbating the instability of its supply. In particular, pre-harvest sprouting (PHS) is one of the major abiotic stresses that frequently occurs due to irregular climate conditions, causing serious damage to wheat and its quality. In this study, transcriptomic analysis with RNA-seq and proteomic analysis with LC-MS/MS were performed in PHS-treated spikes from two wheat cultivars presenting PHS sensitivity and tolerance, respectively. A total of 13,154 differentially expressed genes (DEGs) and 706 differentially expressed proteins (DEPs) were identified in four comparison groups between the susceptible/tolerant cultivars. Gene function and correlation analysis were performed to determine the co-profiled genes and proteins affected by PHS treatment. In the functional annotation of each comparative group, similar functions were confirmed in each cultivar under PHS treatment; however, in Keumgang PHS+7 (K7) vs. Woori PHS+7 (W7), functional annotations presented clear differences in the "spliceosome" and "proteasome" pathways. In addition, our results indicate that alternative splicing and ubiquitin-proteasome support the regulation of germination and seed dormancy. This study provides an advanced understanding of the functions involved in transcription and translation related to PHS mechanisms, thus enabling specific proposals for the further analysis of germination and seed dormancy mechanisms and pathways in wheat.

Genome-wide identification, new classification, expression analysis and screening of drought & heat resistance related candidates in the RING zinc finger gene family of bread wheat (Triticum aestivum L.)

Background

RING (Really Interesting New Gene) zinc finger (RING-zf) proteins belong to an important subclass of zinc fingers superfamily, which play versatile roles during various developmental stages and in abiotic stress responses. Based on the conserved cysteine and histidine residues, the RING-zf domains are classified into RING-HC (C3HC4), RING-H2 (C3H2C3), RING-v, RING-D, RING-S/T, RING-G, and RING-C2. However, little is known about the function of the RING-zfs of wheat.

Results

In this study, 129 (93.5%) of 138 members were found in nucleus, indicating TaRING-zf were primarily engaged in the degradation of transcription factors and other nuclear-localized proteins. 138 TaRING-zf domains can be divided into four canonical or modified types (RING-H2, RING-HC, RING-D, and RING-M). The RING-M was newly identified in T. aestivum, and might represent the intermediate other states between RING-zf domain and other modified domains. The consensus sequence of the RING-M domain can be described as M-X₂-R-X₁₄-Cys-X₁-H-X₂-Cys-X₂-Cys-X₁₀-Cys-X₂-Cys. Further interspecies collinearity analyses showed that TaRING-zfs were more closely related to the genes in Poaceae. According to the public transcriptome data, most of the TaRING-zfs were expressed at different 15 stages of plant growth, development, and some of them exhibited specific responses to drought/heat stress. Moreover, 4 RING-HC (TraesCS2A02G526800.1, TraesCS4A02G290600.1, TraesCS4B02G023600.1 and TraesCS4D02G021200.1) and 2 RING-H2 (TraesCS3A02G288900.1 and TraesCS4A02G174600.1) were significantly expressed at different development stages and under drought stress. These findings provide valuable reference data for further study of their physiological functions in wheat varieties.

Conclusions

Taken together, the characterization and classifications of the TaRING-zf family were extensively studied and some new features about it were revealed. This study could provide some valuable targets for further studies on their functions in growth and development, and abiotic stress responses in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08905-x.

Consensus genomic regions associated with grain protein content in hexaploid and tetraploid wheat

A meta-analysis of QTLs associated with grain protein content (GPC) was conducted in hexaploid and tetraploid wheat to identify robust and stable meta-QTLs (MQTLs). For this purpose, as many as 459 GPC-related QTLs retrieved from 48 linkage-based QTL mapping studies were projected onto the newly developed wheat consensus map. The analysis resulted in the prediction of 57 MQTLs and 7 QTL hotspots located on all wheat chromosomes (except chromosomes 1D and 4D) and the average confidence interval reduced 2.71-fold in the MQTLs and QTL hotspots compared to the initial QTLs. The physical regions occupied by the MQTLs ranged from 140 bp to 224.02 Mb with an average of 15.2 Mb, whereas the physical regions occupied by QTL hotspots ranged from 1.81 Mb to 36.03 Mb with a mean of 8.82 Mb. Nineteen MQTLs and two QTL hotspots were also found to be co-localized with 45 significant SNPs identified in 16 previously published genome-wide association studies in wheat. Candidate gene (CG) investigation within some selected MQTLs led to the identification of 705 gene models which also included 96 high-confidence CGs showing significant expressions in different grain-related tissues and having probable roles in GPC regulation. These significantly expressed CGs mainly involved the genes/gene families encoding for the following proteins: aminotransferases, early nodulin 93, glutamine synthetases, invertase/pectin methylesterase inhibitors, protein BIG GRAIN 1-like, cytochrome P450, glycosyl transferases, hexokinases, small GTPases, UDP-glucuronosyl/UDP-glucosyltransferases, and EamA, SANT/Myb, GNAT, thioredoxin, phytocyanin, and homeobox domains containing proteins. Further, eight genes including GPC-B1, Glu-B1-1b, Glu-1By9, TaBiP1, GSr, TaNAC019-A, TaNAC019-D, and bZIP-TF SPA already known to be associated with GPC were also detected within some of the MQTL regions confirming the efficacy of MQTLs predicted during the current study.

Comprehensive In Silico Analysis of RNA Silencing-Related Genes and Their Regulatory Elements in Wheat (Triticum aestivum L.)

Dicer-like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) are known as the three major gene families that act as the critical components of RNA interference or silencing mechanisms through the noncoding small RNA molecules (miRNA and siRNA) to regulate the expressions of protein-coding genes in eukaryotic organisms. However, most of their characteristics including structures, chromosomal location, subcellular locations, regulatory elements, and gene networking were not rigorously studied. Our analysis identified 7 TaDCL, 39 TaAGO, and 16 TaRDR genes as RNA interference (RNAi) genes from the wheat genome. Phylogenetic analysis of predicted RNAi proteins with the RNAi proteins of Arabidopsis and rice showed that the predicted proteins of TaDCL, TaAGO, and TaRDR groups are clustered into four, eight, and four subgroups, respectively. Domain, 3D protein structure, motif, and exon-intron structure analyses showed that these proteins conserve identical characteristics within groups and maintain differences between groups. The nonsynonymous/synonymous mutation ratio (Ka/Ks) < 1 suggested that these protein sequences conserve some purifying functions. RNAi genes networking with TFs revealed that ERF, MIKC-MADS, C2H2, BBR-BPC, MYB, and Dof are the key transcriptional regulators of the predicted RNAi-related genes. The cis-regulatory element (CREs) analysis detected some important CREs of RNAi genes that are significantly associated with light, stress, and hormone responses. Expression analysis based on an online database exhibited that almost all of the predicted RNAi genes are expressed in different tissues and organs. A case-control study from the gene expression level showed that some RNAi genes significantly responded to the drought and heat stresses. Overall results would therefore provide an excellent basis for in-depth molecular investigation of these genes and their regulatory elements for wheat crop improvement against different stressors.

Molecular Characterization, Evolutionary Analysis, and Expression Profiling of BOR Genes in Important Cereals

Boron (B) is an essential micronutrient of plants. Plants grapple with a narrow range of B between its toxicity and deficiency. B homeostasis mechanism is required to rescue plants from such a quagmire. B transporters are specialized proteins involved in the homeostasis of B. In the present study, a total of 29 BOR genes were identified in five major cereals, including three BORs in each Brachypodium distachyon and Sorghum bicolor, four in Oryza sativa, six in Zea mays, and 13 in Triticum aestivum. Multiple sequence alignments, domain structure analyses, and phylogenetic analysis indicated the conserved nature of the BOR protein family. Duplication events and Ka/Ks analysis of TaBORs showed the role of segmental duplication events and purifying selection in the expansion of the BOR family in T. aestivum. Furthermore, in silico expression and co-expression analyses under biotic and abiotic stress conditions depicted their involvement in combating such conditions. Moreover, qRT-PCR of TaBORs in B treatment suggested the roles of BOR genes in B stress management. The present study hints at the conserved nature of BOR proteins and their different aspects. The study will lay down a way for several crop improvement programs.

Genome-Wide Identification and Functional Characterization of the Chloride Channel TaCLC Gene Family in Wheat (Triticum aestivum L.)

In plants, chloride channels (CLC) are involved in a series of specific functions, such as regulation of nutrient transport and stress tolerance. Members of the wheat Triticum aestivum L. CLC (TaCLC) gene family have been proposed to encode anion channels/transporters that may be related to nitrogen transportation. To better understand their roles, TaCLC family was screened and 23 TaCLC gene sequences were identified using a Hidden Markov Model in conjunction with wheat genome database. Gene structure, chromosome location, conserved motif, and expression pattern of the resulting family members were then analyzed. Phylogenetic analysis showed that the TaCLC family can be divided into two subclasses (I and II) and seven clusters (-a, -c1, -c2, -e, -f1, -f2, and -g2). Using a wheat RNA-seq database, the expression pattern of TaCLC family members was determined to be an inducible expression type. In addition, seven genes from seven different clusters were selected for quantitative real-time PCR (qRT-PCR) analysis under low nitrogen stress or salt stress conditions, respectively. The results indicated that the gene expression levels of this family were up-regulated under low nitrogen stress and salt stress, except the genes of TaCLC-c2 cluster which were from subfamily -c. The yeast complementary experiments illustrated that TaCLC-a-6AS-1, TaCLC-c1-3AS, and TaCLC-e-3AL all had anion transport functions for NO₃⁻ or Cl⁻, and compensated the hypersensitivity of yeast GEF1 mutant strain YJR040w (Δgef1) in restoring anion-sensitive phenotype. This study establishes a theoretical foundation for further functional characterization of TaCLC genes and provides an initial reference for better understanding nitrate nitrogen transportation in wheat.

Gene Co-expression Network Analysis

Gene co-expression analysis is a data analysis technique that helps identify groups of genes with similar expression patterns across several different conditions. By means of these techniques, different groups have been able to assign putative metabolic pathways and functions to understudied genes and to identify novel metabolic regulation networks for different metabolites. Some groups have even used network comparative studies to understand the evolution of these networks from green algae to land plants. In this chapter, we will go over the basic definitions required to understand network topology and gene module identification. Additionally, we offer the reader a walk-through a standard analysis pipeline as implemented in the package WGCNA that takes as input raw fastq files and obtains co-expressed gene clusters and representative gene expression patterns from each module for downstream applications.

Proton-pumping pyrophosphatase homeolog expression is a dynamic trait in bread wheat (Triticum aestivum)

Proton-pumping pyrophosphatases (H+-PPases) have been shown to enhance biomass and yield. However, to date, there has been little work towards identify genes encoding H+-PPases in bread wheat (Triticum aestivum) (TaVPs) and limited knowledge on how the expression of these genes varies across different growth stages and tissue types. In this study, the IWGSC database was used to identify two novel TaVP genes, TaVP4 and TaVP5, and elucidate the complete homeolog sequences of the three known TaVP genes, bringing the total number of bread wheat TaVPs from 9 to 15. Gene expression levels of each TaVP homeolog were assessed using quantitative real-time PCR (qRT-PCR) in four diverse wheat varieties in terms of phenotypic traits related to high vacuolar pyrophosphatase expression. Homeolog expression was analyzed across multiple tissue types and developmental stages. Expression levels of the TaVP homeologs were found to vary significantly between varieties, tissues and plant developmental stages. During early development (Z10 and Z13), expressions of TaVP1 and TaVP2 homeologs were higher in shoot tissue than root tissue, with both shoot and root expression increasing in later developmental stages (Z22). TaVP2-D was expressed in all varieties and tissue types and was the most highly expressed homeolog at all developmental stages. Expression of the TaVP3 homeologs was restricted to developing grain (Z75), while TaVP4 homeolog expression was higher at Z22 than earlier developmental stages. Variation in TaVP4B was detected among varieties at Z22 and Z75, with Buck Atlantico (high biomass) and Scout (elite Australian cultivar) having the highest levels of expression. These findings offer a comprehensive overview of the bread wheat H+-PPase family and identify variation in TaVP homeolog expression that will be of use to improve the growth, yield, and abiotic stress tolerance of bread wheat.

Genome-wide identification and characterization of the Lateral Organ Boundaries Domain (LBD) gene family in polyploid wheat and related species

BACKGROUND

Wheat (Triticum aestivum) originated from three different diploid ancestral grass species and experienced two rounds of polyploidization. Exploring how certain wheat gene subfamilies have expanded during the evolutionary process is of great importance. The Lateral Organ Boundaries Domain (LBD) gene family encodes plant-specific transcription factors that share a highly conserved LOB domain and are prime candidates for this, as they are involved in plant growth, development, secondary metabolism and stress in various species.

METHODS

Using a genome-wide analysis of high-quality polyploid wheat and related species genome sequences, a total of 228 LBD members from five Triticeae species were identified, and phylogenetic relationship analysis of LBD members classified them into two main classes (classes I and II) and seven subgroups (classes I a-e, II a and II b).

RESULTS

The gene structure and motif composition analyses revealed that genes that had a closer phylogenetic relationship in the same subgroup also had a similar gene structure. Macrocollinearity and microcollinearity analyses of Triticeae species suggested that some LBD genes from wheat produced gene pairs across subgenomes of chromosomes 4A and 5A and that the complex evolutionary history of TaLBD4B-9 homologs was a combined result of chromosome translocation, polyploidization, gene loss and duplication events. Public RNA-seq data were used to analyze the expression patterns of wheat LBD genes in various tissues, different developmental stages and following abiotic and biotic stresses. Furthermore, qRT-PCR results suggested that some TaLBDs in class II responded to powdery mildew, regulated reproductive growth and were involved in embryo sac development in common wheat.

Using precision phenotyping to inform de novo domestication

An update on the use of precision phenotyping to assess the potential of lesser cultivated species as candidates for de novo domestication or similar development for future agriculture.

Reduced free asparagine in wheat grain resulting from a natural deletion of TaASN-B2: investigating and exploiting diversity in the asparagine synthetase gene family to improve wheat quality

BACKGROUND

Understanding the determinants of free asparagine concentration in wheat grain is necessary to reduce levels of the processing contaminant acrylamide in baked and toasted wheat products. Although crop management strategies can help reduce asparagine concentrations, breeders have limited options to select for genetic variation underlying this trait. Asparagine synthetase enzymes catalyse a critical step in asparagine biosynthesis in plants and, in wheat, are encoded by five homeologous gene triads that exhibit distinct expression profiles. Within this family, TaASN2 genes are highly expressed during grain development but TaASN-B2 is absent in some varieties.

RESULTS

Natural genetic diversity in the asparagine synthetase gene family was assessed in different wheat varieties revealing instances of presence/absence variation and other polymorphisms, including some predicted to affect the function of the encoded protein. The presence and absence of TaASN-B2 was determined across a range of UK and global common wheat varieties and related species, showing that the deletion encompassing this gene was already present in some wild emmer wheat genotypes. Expression profiling confirmed that TaASN2 transcripts were only detectable in the grain, while TaASN3.1 genes were highly expressed during the early stages of grain development. TaASN-A2 was the most highly expressed TaASN2 homeologue in most assayed wheat varieties. TaASN-B2 and TaASN-D2 were expressed at similar, lower levels in varieties possessing TaASN-B2. Expression of TaASN-A2 and TaASN-D2 did not increase to compensate for the absence of TaASN-B2, so total TaASN2 expression was lower in varieties lacking TaASN-B2. Consequently, free asparagine concentrations in field-produced grain were, on average, lower in varieties lacking TaASN-B2, although the effect was lost when free asparagine accumulated to very high concentrations as a result of sulphur deficiency.

CONCLUSIONS

Selecting wheat genotypes lacking the TaASN-B2 gene may be a simple and rapid way for breeders to reduce free asparagine concentrations in commercial wheat grain.

Designing Future Crops: Genomics-Assisted Breeding Comes of Age

Over the past decade, genomics-assisted breeding (GAB) has been instrumental in harnessing the potential of modern genome resources and characterizing and exploiting allelic variation for germplasm enhancement and cultivar development. Sustaining GAB in the future (GAB 2.0) will rely upon a suite of new approaches that fast-track targeted manipulation of allelic variation for creating novel diversity and facilitate their rapid and efficient incorporation in crop improvement programs. Genomic breeding strategies that optimize crop genomes with accumulation of beneficial alleles and purging of deleterious alleles will be indispensable for designing future crops. In coming decades, GAB 2.0 is expected to play a crucial role in breeding more climate-smart crop cultivars with higher nutritional value in a cost-effective and timely manner.

Expansion and Molecular Characterization of AP2/ERF Gene Family in Wheat (Triticum aestivum L.)

The AP2/ERF is a large protein family of transcription factors, playing an important role in signal transduction, plant growth, development, and response to various stresses. AP2/ERF super-family is identified and functionalized in a different plant but no comprehensive and systematic analysis in wheat (Triticum aestivum L.) has been reported. However, a genome-wide and functional analysis was performed and identified 322 TaAP2/ERF putative genes from the wheat genome. According to the phylogenetic and structural analysis, TaAP2/ERF genes were divided into 12 subfamilies (Ia, Ib, Ic, IIa, IIb, IIc, IIIa, IIIb, IIIc, IVa, IVb, and IVc). Furthermore, conserved motifs and introns/exons analysis revealed may lead to functional divergence within clades. Cis-Acting analysis indicated that many elements were involved in stress-related and plant development. Chromosomal location showed that 320 AP2/ERF genes were distributed among 21 chromosomes and 2 genes were present in a scaffold. Interspecies microsynteny analysis revealed that maximum orthologous between Arabidopsis, rice followed by wheat. Segment duplication events have contributed to the expansion of the AP2/ERF family and made this family larger than rice and Arabidopsis. Additionally, AP2/ERF genes were differentially expressed in wheat seedlings under the stress treatments of heat, salt, and drought, and expression profiles were verified by qRT-PCR. Remarkably, the RNA-seq data exposed that AP2/ERF gene family might play a vital role in stress-related. Taken together, our findings provided useful and helpful information to understand the molecular mechanism and evolution of the AP2/ERF gene family in wheat.

Extensive structural variation in the Bowman-Birk inhibitor family in common wheat (Triticum aestivum L.)

BACKGROUND

Bowman-Birk inhibitors (BBI) are a family of serine-type protease inhibitors that modulate endogenous plant proteolytic activities during different phases of development. They also inhibit exogenous proteases as a component of plant defense mechanisms, and their overexpression can confer resistance to phytophagous herbivores and multiple fungal and bacterial pathogens. Dicot BBIs are multifunctional, with a "double-headed" structure containing two separate inhibitory loops that can bind and inhibit trypsin and chymotrypsin proteases simultaneously. By contrast, monocot BBIs have a non-functional chymotrypsin inhibitory loop, although they have undergone internal duplication events giving rise to proteins with multiple BBI domains.

RESULTS

We used a Hidden Markov Model (HMM) profile-based search to identify 57 BBI genes in the common wheat (Triticum aestivum L.) genome. The BBI genes are unevenly distributed, with large gene clusters in the telomeric regions of homoeologous group 1 and 3 chromosomes that likely arose through a series of tandem gene duplication events. The genomes of wheat progenitors also contain contiguous clusters of BBI genes, suggesting this family underwent expansion before the domestication of common wheat. However, the BBI gene family varied in size among different cultivars, showing this family remains dynamic. Because of these expansions, the BBI gene family is larger in wheat than other monocots such as maize, rice and Brachypodium. We found BBI proteins in common wheat with intragenic homologous duplications of cysteine-rich functional domains, including one protein with four functional BBI domains. This diversification may expand the spectrum of target substrates. Expression profiling suggests that some wheat BBI proteins may be involved in regulating endogenous proteases during grain development, while others were induced in response to biotic and abiotic stresses, suggesting a role in plant defense.

CONCLUSIONS

Genome-wide characterization reveals that the BBI gene family in wheat is subject to a high rate of homologous tandem duplication and deletion events, giving rise to a diverse set of encoded proteins. This information will facilitate the functional characterization of individual wheat BBI genes to determine their role in wheat development and stress responses, and their potential application in breeding.

Identification of individual components of a commercial wheat germ acid phosphatase preparation

Wheat germ acid phosphatase (WGAP) is a commercial preparation of partially purified protein commonly used in laboratory settings for non-specific enzymatic dephosphorylation. It is known that these preparations contain multiple phosphatase isozymes and are still relatively crude. This study therefore aimed to identify the protein components of a commercial preparation of wheat germ acid phosphatase using mass spectroscopy and comparative genomics. After one post-purchase purification step, the most prevalent fifteen proteins in the mixture included heat shock proteins, beta-amylases, glucoseribitol dehydrogenases, enolases, and an aminopeptidase. While not among the most abundant components, eight unique dephosphorylation enzymes were also present including three purple acid phosphatases. Furthermore, it is shown that some of these correspond to previously isolated isozymes; one of which has been also previously shown by transcriptome data to be overexpressed in wheat seeds. In summary, this study identified the major components of WGAP including phosphatases and hypothesizes the most active components towards a better understanding of this commonly used laboratory tool.

A wheat cysteine-rich receptor-like kinase confers broad-spectrum resistance against Septoria tritici blotch

The poverty of disease resistance gene reservoirs limits the breeding of crops for durable resistance against evolutionary dynamic pathogens. Zymoseptoria tritici which causes Septoria tritici blotch (STB), represents one of the most genetically diverse and devastating wheat pathogens worldwide. No fully virulent Z. tritici isolates against synthetic wheats carrying the major resistant gene Stb16q have been identified. Here, we use comparative genomics, mutagenesis and complementation to identify Stb16q, which confers broad-spectrum resistance against Z. tritici. The Stb16q gene encodes a plasma membrane cysteine-rich receptor-like kinase that was recently introduced into cultivated wheat and which considerably slows penetration and intercellular growth of the pathogen.

Bulked Segregant RNA-Seq Provides Distinctive Expression Profile Against Powdery Mildew in the Wheat Genotype YD588

Wheat powdery mildew, caused by the fungal pathogen Blumeria graminis f. sp. tritici (Bgt), is a destructive disease leading to huge yield losses in production. Host resistance can greatly contribute to the control of the disease. To explore potential genes related to the powdery mildew (Pm) resistance, in this study, we used a resistant genotype YD588 to investigate the potential resistance components and profiled its expression in response to powdery mildew infection. Genetic analysis showed that a single dominant gene, tentatively designated PmYD588, conferred resistance to powdery mildew in YD588. Using bulked segregant RNA-Seq (BSR-Seq) and single nucleotide polymorphism (SNP) association analysis, two high-confidence candidate regions were detected in the chromosome arm 2B, spanning 453,752,054-506,356,791 and 584,117,809-664,221,850 bp, respectively. To confirm the candidate region, molecular markers were developed using the BSR-Seq data and mapped PmYD588 to an interval of 4.2 cM by using the markers YTU588-004 and YTU588-008. The physical position was subsequently locked into the interval of 647.1-656.0 Mb, which was different from those of Pm6, Pm33, Pm51, Pm52, Pm63, Pm64, PmQ, PmKN0816, MlZec1, and MlAB10 on the same chromosome arm in its position, suggesting that it is most likely a new Pm gene. To explore the potential regulatory genes of the R gene, 2,973 differentially expressed genes (DEGs) between the parents and bulks were analyzed using gene ontology (GO), clusters of orthologous group (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Based on the data, we selected 23 potential regulated genes in the enriched pathway of plant-pathogen interaction and detected their temporal expression patterns using an additional set of wheat samples and time-course analysis postinoculation with Bgt. As a result, six disease-related genes showed distinctive expression profiles after Bgt invasion and can serve as key candidates for the dissection of resistance mechanisms and improvement of durable resistance to wheat powdery mildew.

Identification of Rf9, a Gene Contributing to the Genetic Complexity of Fertility Restoration in Hybrid Wheat

Wheat (Triticum aestivum L.) is a self-pollinating crop whose hybrids offer the potential to provide a major boost in yield. Male sterility induced by the cytoplasm of Triticum timopheevii is a powerful method for hybrid seed production. Hybrids produced by this method are often partially sterile, and full fertility restoration is crucial for wheat production using hybrid cultivars. To identify the genetic loci controlling fertility restoration in wheat, we produced two cytoplasmic male-sterile (CMS) backcross (BC1) mapping populations. The restorer lines Gerek 79 and 71R1203 were used to pollinate the male-sterile winter wheat line CMS-Sperber. Seed set and numbers of sterile spikelets per spike were evaluated in 340 and 206 individuals of the populations derived from Gerek 79 and 71R1203, respectively. Genetic maps were constructed using 930 and 994 single nucleotide polymorphism (SNP) markers, spanning 2,160 and 2,328 cM over 21 linkage groups in the two populations, respectively. Twelve quantitative trait loci (QTL) controlled fertility restoration in both BC1 populations, including a novel restorer-of-fertility (Rf) locus flanked by the SNP markers IWB72413 and IWB1550 on chromosome 6AS. The locus was mapped as a qualitative trait in the BC1 Gerek 79 population and was designated Rf9. One hundred-nineteen putative candidate genes were predicted within the QTL region on chromosome 6AS. Among them were genes encoding mitochondrial transcription termination factor and pentatricopeptide repeat-containing proteins that are known to be associated with fertility restoration. This finding is a promising step to better understand the functions of genes for improving fertility restoration in hybrid wheat.

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

Plant genome sequencing has dramatically increased, and some species even have multiple high-quality reference versions. Demands for clade-specific homology inference and analysis have increased in the pangenomic era. Here we present a novel method, GeneTribe (https://chenym1.github.io/genetribe/), for homology inference among genetically similar genomes that incorporates gene collinearity and shows better performance than traditional sequence-similarity-based methods in terms of accuracy and scalability. The Triticeae tribe is a typical allopolyploid-rich clade with complex species relationships that includes many important crops, such as wheat, barley, and rye. We built Triticeae-GeneTribe (http://wheat.cau.edu.cn/TGT/), a homology database, by integrating 12 Triticeae genomes and 3 outgroup model genomes and implemented versatile analysis and visualization functions. With macrocollinearity analysis, we were able to construct a refined model illustrating the structural rearrangements of the 4A-5A-7B chromosomes in wheat as two major translocation events. With collinearity analysis at both the macro- and microscale, we illustrated the complex evolutionary history of homologs of the wheat vernalization gene Vrn2, which evolved as a combined result of genome translocation, duplication, and polyploidization and gene loss events. Our work provides a useful practice for connecting emerging genome assemblies, with awareness of the extensive polyploidy in plants, and will help researchers efficiently exploit genome sequence resources.

Genome-wide characterization and expression and co-expression analysis suggested diverse functions of WOX genes in bread wheat

WUSCHEL-related homeobox (WOX) genes belong to the homeobox superfamily, are plant-specific and play vital functions in the growth and development. Herein, we identified a total of 43 TaWOX genes in the allohexaploid (AABBDD) genome of Triticum aestivum L. These genes were distributed on the various chromosomes of each subgenome (A, B and D). The phylogenetic analysis showed the clustering of TaWOXs into three clades: ancient, intermediate and modern or WUS. The gene and protein structures including exon/intron organization, intron phases, and domain and motif distribution were found to be conserved in each phylogenetic clade. The subcellular localization was predicted as nuclear. The Ka/Ks analyses suggested the purifying selection of paralogous genes. The differential expression profiling of various TaWOXs in numerous tissue developmental stages and different layers of grains suggested their role in growth and development. Moreover, a few genes exhibited modulated expression during abiotic and biotic stress conditions, which revealed their roles in stress response. The occurrence of various cis-acting regulatory elements further confirmed their role in plant development and stress tolerance. The co-expression analyses suggested the interactions of these genes with other genes, involved in various processes including plant development, signalling and stress responses. The present study reported several characteristic features of TaWOXs genes that can be useful for further characterization in future studies.

The NAC transcription factor NAC019-A1 is a negative regulator of starch synthesis in wheat developing endosperm

Starch is a major component of wheat (Triticum aestivum L.) endosperm and is an important part of the human diet. The functions of many starch synthesis genes have been elucidated. However, little is known about their regulatory mechanisms in wheat. Here, we identified a novel NAC transcription factor, TaNAC019-A1 (TraesCS3A02G077900), that negatively regulates starch synthesis in wheat and rice (Oryza sativa L.) endosperms. TaNAC019-A1 was highly expressed in the endosperm of developing grains and encoded a nucleus-localized transcriptional repressor. Overexpression of TaNAC019-A1 in rice and wheat led to significantly reduced starch content, kernel weight, and kernel width. The TaNAC019-A1-overexpression wheat lines had smaller A-type starch granules and fewer B-type starch granules than wild-type. Moreover, TaNAC019-A1 could directly bind to the 'ACGCAG' motif in the promoter regions of ADP-glucose pyrophosphorylase small subunit 1 (TaAGPS1-A1, TraesCS7A02G287400) and TaAGPS1-B1 (TraesCS7B02G183300) and repress their expression, thereby inhibiting starch synthesis in wheat endosperm. One haplotype of TaNAC019-B1 (TaNAC019-B1-Hap2, TraesCS3B02G092800) was positively associated with thousand-kernel weight and underwent positive selection during the Chinese wheat breeding process. Our data demonstrate that TaNAC019-A1 is a negative regulator of starch synthesis in wheat endosperm and provide novel insight into wheat yield improvement.

Genome Wide Identification, Characterization, and Expression Analysis of YABBY-Gene Family in WHEAT (Triticum aestivum L.)

The small YABBY plant-specific transcription factor has a prominent role in regulating plant growth and developmental activities. However, little information is available about YABBY gene family in Triticum aestivum L. Herein, we identified 21 TaYABBY genes in the Wheat genome database. Then, we performed the conserved motif and domain analysis of TaYABBY proteins. The phylogeny of the TaYABBY was further sub-divided into 6 subfamilies (YABBY1/YABBY3, YABB2, YABBY5, CRC and INO) based on the structural similarities and functional diversities. The GO (Gene ontology) analysis of TaYABBY proteins showed that they are involved in numerous developmental processes and showed response against environmental stresses. The analysis of all identified genes in RNA-seq data showed that they are expressed in different tissues of wheat. Differential expression patterns were observed in not only control samples but also in stressed samples such as biotic stress (i.e., Fusarium graminearum (F.g), septoria tritici (STB), Stripe rust (Sr) and Powdery mildew (Pm), and abiotic stress (i.e., drought, heat, combined drought and heat and phosphorus deficiency), especially at different grain development stages. All identified TaYABBY-genes were localized in the nucleus which implies their participation in the regulatory mechanisms of various biological and cellular processes. In light of the above-mentioned outcomes, it has been deduced that TaYABBY-genes in the wheat genome play an important role in mediating various development, growth, and resistance mechanism, which could provide significant clues for future functional studies.

Bulked Segregant RNA-Seq Reveals Distinct Expression Profiling in Chinese Wheat Cultivar Jimai 23 Responding to Powdery Mildew

Wheat powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is one of the most destructive fungal diseases threatening global wheat production. Host resistance is well known to be the most efficient method to control this disease. However, the molecular mechanism of wheat powdery mildew resistance (Pm) is still unclear. To analyze the molecular mechanism of Pm, we used the resistant wheat cultivar Jimai 23 to investigate its potential resistance components and profiled its expression in response to powdery mildew infection using bulked segregant RNA-Seq (BSR-Seq). We showed that the Pm of Jimai 23 was provided by a single dominant gene, tentatively designated PmJM23, and assigned it to the documented Pm2 region of chromosome 5DS. 3,816 consistently different SNPs were called between resistant and susceptible parents and the bulked pools derived from the combinations between the resistant parent Jimai23 and the susceptible parent Tainong18. 58 of the SNPs were assigned to the candidate region of PmJM23. Subsequently, 3,803 differentially expressed genes (DEGs) between parents and bulks were analyzed by GO, COG and KEGG pathway enrichment. The temporal expression patterns of associated genes following Bgt inoculation were further determined by RT-qPCR. Expression of six disease-related genes was investigated during Bgt infection and might serve as valuable genetic resources for the improvement of durable resistance to Bgt.

Genome-Wide Identification and Analysis of the Growth-Regulating Factor (GRF) Gene Family and GRF-Interacting Factor Family in Triticum aestivum L

Growth-regulating factors (GRFs) are unique transcription factors in plants. GRFs can interact with SNH (SYT N-terminal homology) domains in GRF-interacting factor (GIF) proteins via the N-terminal QLQ (Gln, Leu, Gln) domain to form functional complexes and participate in the regulation of downstream gene expression. In this study, we systematically identified the GRF gene family and GIF gene family in wheat and its relatives comprising Triticum urartu, Triticum dicoccoides, and Aegilops tauschii. Thirty GRF gene members are present in wheat, which are distributed on 12 chromosomes and they have 2-5 protein-coding regions. They all contain QLQ and WRC (Trp, Arg, Cys) conserved domains. Wheat possesses only eight members of the GIF gene family, which are distributed on six chromosomes. All wheat GIF (TaGIF) proteins have highly conserved SNH and QG (Gln, Gly) domains. The wheat GRF (TaGRF) gene family has 13 pairs of segmental duplication genes and no tandem duplication genes; the TaGIF gene family has two pairs of segmental duplication genes and no tandem duplication genes. It is speculated that segmental duplication events may be the main reason for the amplification of TaGRF gene family and TaGIF gene family. Based on published transcriptome data and qRT-PCR results of 8 TaGRF genes and 4 TaGIF genes, all of the genes responded strongly to osmotic stress, and the expression levels of TaGRF21 and TaGIF5 were also significantly upregulated under drought and cold stress conditions. The results obtained in this study may facilitate further investigations of the functions of TaGRF genes and TaGIF genes in order to identify candidate genes for use in stress-resistant wheat breeding programs.

A roadmap for gene functional characterisation in crops with large genomes: Lessons from polyploid wheat

Understanding the function of genes within staple crops will accelerate crop improvement by allowing targeted breeding approaches. Despite their importance, a lack of genomic information and resources has hindered the functional characterisation of genes in major crops. The recent release of high-quality reference sequences for these crops underpins a suite of genetic and genomic resources that support basic research and breeding. For wheat, these include gene model annotations, expression atlases and gene networks that provide information about putative function. Sequenced mutant populations, improved transformation protocols and structured natural populations provide rapid methods to study gene function directly. We highlight a case study exemplifying how to integrate these resources. This review provides a helpful guide for plant scientists, especially those expanding into crop research, to capitalise on the discoveries made in Arabidopsis and other plants. This will accelerate the improvement of crops of vital importance for food and nutrition security.

Genome-wide and SNP network analyses reveal genetic control of spikelet sterility and yield-related traits in wheat

Revealing the genetic factors underlying yield and agronomic traits in wheat are an imperative need for covering the global food demand. Yield boosting requires a deep understanding of the genetic basis of grain yield-related traits (e.g., spikelet fertility and sterility). Here, we have detected much natural variation among ancient hexaploid wheat accessions in twenty-two agronomic traits collected over eight years of field experiments. A genome-wide association study (GWAS) using 15 K single nucleotide polymorphisms (SNPs) was applied to detect the genetic basis of studied traits. Subsequently, the GWAS output was reinforced via other statistical and bioinformatics analyses to detect putative candidate genes. Applying the genome-wide SNP-phenotype network defined the most decisive SNPs underlying the traits. Six pivotal SNPs, co-located physically within the genes encoding enzymes, hormone response, metal ion transport, and response to oxidative stress have been identified. Of these, metal ion transport and Gibberellin 2-oxidases (GA2oxs) genes showed strong involvement in controlling the spikelet sterility, which had not been reported previously in wheat. SNP-gene haplotype analysis confirmed that these SNPs influence spikelet sterility, especially the SNP co-located on the exon of the GA2ox gene. Interestingly, these genes were highly expressed in the grain and spike, demonstrating their pivotal role in controlling the trait. The integrative analysis strategy applied in this study, including GWAS, SNP-phenotype network, SNP-gene haplotype, expression analysis, and genome-wide prediction (GP), empower the identification of functional SNPs and causal genes. GP outputs obtained in this study are encouraging for the implementation of the traits to accelerate yield improvement by making an early prediction of complex yield-related traits in wheat. Our findings demonstrate the usefulness of the ancient wheat material as a valuable resource for yield-boosting. This is the first comprehensive genome-wide analysis for spikelet sterility in wheat, and the results provide insights into yield improvement.

Molecular Characterization of Auxin Efflux Carrier- ABCB1 in hexaploid wheat

Auxin is an important phytohormone that regulates response, differentiation, and development of plant cell, tissue, and organs. Along with its local production, long-distance transport coordinated by the efflux/influx membrane transporters is instrumental in plant development and architecture. In the present study, we cloned and characterized a wheat (Triticum aestivum) auxin efflux carrier ABCB1. The TaABCB1 was physically localized to the proximal 15% of the short arm of wheat homoeologous group 7 chromosomes. Size of the Chinese spring (CS) homoeologs genomic copies ranged from 5.3–6.2 kb with the 7A copy being the largest due to novel insertions in its third intron. The three homoeologous copies share 95–97% sequence similarity at the nucleotide, 98–99% amino acid, and overall Q-score of 0.98 at 3-D structure level. Though detected in all analyzed tissues, TaABCB1 predominantly expressed in the meristematic tissues likely due to the presence of meristem-specific activation regulatory element identified in the promoter region. RNAi plants of TaABCB1 gene resulted in reduced plant height and increased seed width. Promoter analysis revealed several responsive elements detected in the promoter region including that for different hormones as auxin, gibberellic acid, jasmonic acid and abscisic acid, light, and circadian regulated elements.

Development of mlo-based resistance in tetraploid wheat against wheat powdery mildew

Powdery mildew is a severe disease in wheat. In barley, durable resistance exists, based on non-functionality of the Mlo gene. As a model to analyse the effects of mutagenesis in the homoeologous Mlo genes of wheat, we developed mlo-based powdery mildew resistance in tetraploid durum wheat. To obtain Mlo mutations, we screened a TILLING population developed in tetraploid wheat "Kronos" for which the captured exome sequence of  > 1500 lines is available. This resulted in 23 mutants for Mlo-A1 and 26 non-redundant mutants for Mlo-B1. Two Mlo-A1 and four Mlo-B1 mutants were crossed to obtain eight F₂ mutant lines that showed a range of phenotypes from susceptibility to full resistance. Pot experiments under semi-field conditions confirmed the resistance levels for six of the mutants without any signs of adverse pleiotropic effects. Resistance ranking was similar across six powdery mildew isolates, indicating no isolate specificity of the mlo-based resistance. The effect of mutations in the Mlo-B1 gene was stronger than in the Mlo-A1 gene, probably reflecting differences in wild-type Mlo gene expression levels. Strong partial resistance effects were observed with single mlo-B1 mutations hence, revealing a dosage effect of mlo mutant alleles. Two of the four mlo-B1 mutations (W163* and P335L) were very strong; however, the highest combined effect was observed with the MloA-P335S/MloB-P335L combination, suggesting that non-functional, but full-length Mlo proteins might have the strongest effect compared with nonsense mutations. Our results show that mlo-based resistance might offer possibilities to introduce durable protection in tetraploid wheat against powdery mildew.

Probing the effect of a plus 1bp frameshift mutation in protein-DNA interface of domestication gene, NAMB1, in wheat

Transcription factor NAM-B1 has a major role in the process of senescence, which results in higher Fe and Zn concentrations in grains of wild wheat (T. durum; Td). The absence of the wild type NAMB1 in T. aestivum (Ta), one of the cardinal crops essential for more than 1/3^rd of the global population, affects Fe and Zn remobilisation to the maturing grain from the flag leaf resulting in lesser micronutrient bioavailability. The cardinal difference in the NAMB1 gene between the two species is the absence of +1 bp allele in Ta. Insilico studies using NAMB1 from Td and Ta was performed to explore the variation in the interaction with the conserved cis-element DNA motif (CATGTG) as both the proteins share the same domain, but there are no in silico studies reported of these proteins. The secondary structure, 3D-modelling of the proteins, DNA-protein docking and dynamics have computed by Schrodinger Prime Suite. Predicted secondary structures were energy minimised using Macromodel and docking was performed based on binding energy and hydrogen bonds. Molecular dynamics simulation of NAMB1-Ta and NAMB1-Td individually and with the cis-element motif, performed for 100 ns, revealed significant variations in the protein-DNA interaction in Ta. This work provides the modelled 3D-interaction profile caused by a single bp frameshift mutation in understanding the difference in function between NAMB1 orthologs due to lack of NAC domain. The overall computational analysis reveals that NAMB1-Ta and NAMB1-Td proteins display a good amount of dissimilarity in their structure, dynamics and DNA-binding characteristics.Communicated by Ramaswamy H. Sarma.

Genome-Wide Identification and Analysis of the AP2 Transcription Factor Gene Family in Wheat (Triticum aestivum L.)

The AP2 transcription factors play important roles in regulating plant growth and development. However, limited data are available on the contributions of AP2 transcription factors in wheat (Triticum aestivum L.). In the present study, a total of 62 AP2 genes were identified in wheat from a genome-wide search against the latest wheat genome data. Phylogenetic and sequence alignment analyses divided the wheat AP2 genes into 3 clusters, euAP2, euANT, and basalANT. Chromosomal distribution, gene structure and duplication, and motif composition were subsequently investigated. The 62 TaAP2 genes were unevenly distributed on 21 chromosomes. Twenty-four homologous gene sets among A, B, and D sub-genomes were detected, which contributed to the expansion of the wheat AP2 gene family. The expression levels of TaAP2 genes were examined using the WheatExp database; most detected genes exhibited tissue-specific expression patterns. The transcript levels of 9 randomly selected TaAP2 genes were validated through qPCR analyses. Overexpression of TaAP2-10-5D, the most likely homolog of Arabidopsis ANT gene, increased organ sizes in Arabidopsis. Our results extend our knowledge of the AP2 gene family in wheat, and contribute to further functional characterization of AP2s during wheat development with the ultimate goal of improving crop production.

Reduced stomatal density in bread wheat leads to increased water-use efficiency

Wheat is a staple crop, frequently cultivated in water-restricted environments. Improving crop water-use efficiency would be desirable if grain yield can be maintained. We investigated whether a decrease in wheat stomatal density via the manipulation of epidermal patterning factor (EPF) gene expression could improve water-use efficiency. Our results show that severe reductions in stomatal density in EPF-overexpressing wheat plants have a detrimental outcome on yields. However, wheat plants with a more moderate reduction in stomatal density (i.e. <50% reduction in stomatal density on leaves prior to tillering) had yields indistinguishable from controls, coupled with an increase in intrinsic water-use efficiency. Yields of these moderately reduced stomatal density plants were also comparable with those of control plants under conditions of drought and elevated CO2. Our data demonstrate that EPF-mediated control of wheat stomatal development follows that observed in other grasses, and we identify the potential of stomatal density as a tool for breeding wheat plants that are better able to withstand water-restricted environments without yield loss.

The genome-wide impact of cadmium on microRNA and mRNA expression in contrasting Cd responsive wheat genotypes

BACKGROUND

Heavy metal ATPases (HMAs) are responsible for Cd translocation and play a primary role in Cd detoxification in various plant species. However, the characteristics of HMAs and the regulatory mechanisms between HMAs and microRNAs in wheat (Triticum aestivum L) remain unknown.

RESULTS

By comparative microRNA and transcriptome analysis, a total three known and 19 novel differentially expressed microRNAs (DEMs) and 1561 differentially expressed genes (DEGs) were found in L17 after Cd treatment. In H17, by contrast, 12 known and 57 novel DEMs, and only 297 Cd-induced DEGs were found. Functional enrichments of DEMs and DEGs indicate how genotype-specific biological processes responded to Cd stress. Processes found to be involved in microRNAs-associated Cd response include: ubiquitin mediated proteolysis, tyrosine metabolism, and carbon fixation pathways and thiamine metabolism. For the mRNA response, categories including terpenoid backbone biosynthesis and phenylalanine metabolism, and photosynthesis - antenna proteins and ABC transporters were enriched. Moreover, we identified 32 TaHMA genes in wheat. Phylogenetic trees, chromosomal locations, conserved motifs and expression levels in different tissues and roots under Cd stress are presented. Finally, we infer a microRNA-TaHMAs expression network, indicating that miRNAs can regulate TaHMAs.

CONCLUSION

Our findings suggest that microRNAs play important role in wheat under Cd stress through regulation of targets such as TaHMA2;1. Identification of these targets will be useful for screening and breeding low-Cd accumulation wheat lines.

Genome-wide identification, phylogenetic and expression analysis of the heat shock transcription factor family in bread wheat (Triticum aestivum L.)

Background

Environmental toxicity from non-essential heavy metals such as cadmium (Cd), which is released from human activities and other environmental causes, is rapidly increasing. Wheat can accumulate high levels of Cd in edible tissues, which poses a major hazard to human health. It has been reported that heat shock transcription factor A 4a (HsfA4a) of wheat and rice conferred Cd tolerance by upregulating metallothionein gene expression. However, genome-wide identification, classification, and comparative analysis of the Hsf family in wheat is lacking. Further, because of the promising role of Hsf genes in Cd tolerance, there is need for an understanding of the expression of this family and their functions on wheat under Cd stress. Therefore, here we identify the wheat TaHsf family and to begin to understand the molecular mechanisms mediated by the Hsf family under Cd stress.

Results

We first identified 78 putative Hsf homologs using the latest available wheat genome information, of which 38 belonged to class A, 16 to class B and 24 to class C subfamily. Then, we determined chromosome localizations, gene structures, conserved protein motifs, and phylogenetic relationships of these TaHsfs. Using RNA sequencing data over the course of development, we surveyed expression profiles of these TaHsfs during development and under different abiotic stresses to characterise the regulatory network of this family. Finally, we selected 13 TaHsf genes for expression level verification under Cd stress using qRT-PCR.

Conclusions

To our knowledge, this is the first report of the genome organization, evolutionary features and expression profiles of the wheat Hsf gene family. This work therefore lays the foundation for targeted functional analysis of wheat Hsf genes, and contributes to a better understanding of the roles and regulatory mechanism of wheat Hsfs under Cd stress.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5876-x) contains supplementary material, which is available to authorized users.