Deep transcriptome sequencing provides new insights into the structural and functional organization of the wheat genome

Genome-wide analysis of salicylic acid and jasmonic acid signalling marker gene families in wheat

Jasmonic acid (JA) and salicylic acid (SA) phytohormone pathways are important regulators of stress tolerance. Knowledge regarding the diversity, phylogeny and functionality of wheat genes involved in JA and SA response is limited. Using Arabidopsis, rice and wheat genomic and wheat disease transcriptomic data, we deduced the size, phylogenetic diversity and pathogen-responsiveness of seven hormone-responsive gene families, and thus selected 14 candidates as potential hormone responsive gene markers. Gene-specific expression studies assessed the impact of exogenous JA and SA on their transcriptional activation in leaves of two distinct wheat cultivars. RNAseq data were interrogated to assess their disease responsiveness and tissue-specific expression. This study elucidated the number, phylogeny and pathogen-responsiveness of wheat genes from seven families, including 12 TaAOS, 6 TaJAMyb, 256 TaWRKY group III, 85 TaPR1, 205 TaPR2, 76 TaPR3 and 124 TaPR5. This included the first description of the wheat AOS, JAMyb, PR2, PR3 and PR5 gene families. Gene expression studies delineated TaAOS1-5B and TaJAMyb-4A as JA-responsive in leaves, but not significantly responsive to SA treatment, while TaWRKY45-B was a SA- but not a JA-responsive marker. Other candidate genes were either unresponsive or non-specific to SA or JA. Our findings highlight that all seven gene families are greatly expanded in wheat as compared to other plants (up to 7.6-fold expansion), and demonstrate disparity in the response to biotic stress between some homoeologous and paralogous sequences within these families. The SA- and JA-responsive marker genes identified herein will prove useful tools to monitor these signalling pathways in wheat.

Investigation of two-pore K+ (TPK) channels in Triticum aestivum L. suggests their role in stress response

Two-pore K+ (TPK) channels are voltage-independent and involved in stress response in plants. Herein, we identified 12 TaTPK genes located on nine chromosomes in the Triticum aestivum genome. The majority of TaTPK genes comprised two exons. Each TaTPK channel comprised four transmembrane (TM) helices, N- and C-terminal ion-channel domains, two EF-hand domains and one 14-3-3 binding site. Additionally, highly conserved 'GYGD' motif responsible for K+ ion specificity, was found in between the TMs in both the ion-channel domains. Nine TaTPK channels were predicted to be localised at the plasma membrane, while three were vacuolar. The protein-protein and protein-chemical interactions indicated the coordinated functioning of the TaTPK channels with the other K+ transporters and their possible interaction with the Ca2+-signaling pathway. Expression studies suggested their importance in both vegetative and reproductive tissues development. Significantly modulated expression of various TaTPK genes during heat, drought, combined heat and drought and salt stresses, and after fungal infestation, depicted their function in stress responses. The miRNAs and transcription factors interaction analyses suggested their role in the hormone, light, growth and development-related, and stress-responsive signaling cascades. The current study suggested vital functions of various TaTPK genes, especially in stress response, and would provide an opportunity for their detailed characterization in future studies.

Genome wide characterization and expression analysis of CrRLK1L gene family in wheat unravels their roles in development and stress-specific responses

Catharanthus roseus receptor-like kinase 1-like (CrRLK1L) genes encode a subfamily of receptor-like kinases (RLK) that regulate diverse processes during plant growth, development, and stress responses. The first CrRLK1L was identified from the Catharanthus roseus, commonly known as Madagascar periwinkle. Subsequently, CrRLK1L gene families have been characterized in many plants. The genome of T. aestivum encodes 15 CrRLK1L genes with 43 paralogous copies, with three homeologs each, except for -2-D and -7-A, which are absent. Chromosomal localization analysis revealed a markedly uneven distribution of CrRLK1L genes across seven different chromosomes, with chromosome 4 housing the highest number of genes, while chromosome 6 lacked any CrRLK1L genes. Tissue-specific gene expression analysis revealed distinct expression patterns among the gene family members, with certain members exhibiting increased expression in reproductive tissues. Gene expression analysis in response to various abiotic and biotic stress conditions unveiled differential regulation of gene family members. Cold stress induces CrRLK1Ls -4-B and -15-A while downregulating -3-A and -7B. Drought stress upregulates -9D, contrasting with the downregulation of -7D. CrRLK1L-15-B and -15-D were highly induced in response to 1 hr of heat, and combined drought and heat stress, whereas -10-B is downregulated. Similarly, in response to NaCl stress, only CrRLK1L1 homeologs were induced. Fusarium graminearum and Claviceps purpurea inoculation induces homeologs of CrRLK1L-6 and -7. The analysis of cis-acting elements in the promoter regions identified elements crucial for plant growth and developmental processes. This comprehensive genome-wide analysis and expression study provides valuable insights into the essential functions of CrRLK1L members in wheat.

Modulation in gene expression and enzyme activity suggested the roles of monodehydroascorbate reductase in development and stress response in bread wheat

Monodehydroascorbate reductase (MDHAR) is a crucial enzymatic antioxidant of the ascorbate-glutathione pathway involved in reactive oxygen species scavenging. Herein, we identified 15 TaMDHAR genes in bread wheat. Phylogenetic analysis revealed their clustering into three groups, which are also related to the subcellular localization in the peroxisome matrix, peroxisome membrane, and chloroplast. Each TaMDHAR protein consisted of two conserved domains; Pyr_redox and Pyr_redox_2 of the pyridine nucleotide disulfide oxidoreductase family. The occurrence of diverse groups of cis-regulatory elements in the promoter region and their interaction with numerous transcription factors suggest assorted functions of TaMDHARs in growth and development and in light, phytohormones, and stress responses. Expression analysis in various tissues further revealed their importance in vegetative and reproductive development. In addition, the differential gene expression and enhanced enzyme activity during drought, heat, and salt treatments exposed their role in abiotic stress response. Interaction of MDHARs with various antioxidant enzymes and biochemicals related to the ascorbate-glutathione cycle exposed their synchronized functioning. Interaction with auxin indicated the probability of cross-talk between antioxidants and auxin signaling. The miR168a, miR169, miR172 and others interaction with various TaMDHARs further directed their association with developmental processes and stress responses. The current study provides extensive information about the importance of TaMDHARs, moreover, the precise role of each gene needs to be established in future studies.

Plant Tolerance to Drought Stress with Emphasis on Wheat

Environmental stresses, such as drought, have negative effects on crop yield. Drought is a stress whose impact tends to increase in some critical regions. However, the worldwide population is continuously increasing and climate change may affect its food supply in the upcoming years. Therefore, there is an ongoing effort to understand the molecular processes that may contribute to improving drought tolerance of strategic crops. These investigations should contribute to delivering drought-tolerant cultivars by selective breeding. For this reason, it is worthwhile to review regularly the literature concerning the molecular mechanisms and technologies that could facilitate gene pyramiding for drought tolerance. This review summarizes achievements obtained using QTL mapping, genomics, synteny, epigenetics, and transgenics for the selective breeding of drought-tolerant wheat cultivars. Synthetic apomixis combined with the msh1 mutation opens the way to induce and stabilize epigenomes in crops, which offers the potential of accelerating selective breeding for drought tolerance in arid and semi-arid regions.

Transcriptome meta-analysis of abiotic stresses-responsive genes and identification of candidate transcription factors for broad stress tolerance in wheat

Under field conditions, wheat is subjected to single or multiple stress conditions. The elucidation of the molecular mechanism of stress response is a key step to identify candidate genes for stress resistance in plants. In this study, RNA-seq data analysis identified 17.324, 10.562, 5.510, and 8.653 differentially expressed genes (DEGs) under salt, drought, heat, and cold stress, respectively. Moreover, the comparison of DEGs from each stress revealed 2374 shared genes from which 40% showed highly conserved expression patterns. Moreover, co-expression network analysis and GO enrichment revealed co-expression modules enriched with genes involved in transcription regulation, protein kinase pathway, and genes responding to phytohormones or modulating hormone levels. The expression of 15 selected transcription factor encoding genes was analyzed under abiotic stresses and ABA treatment in durum wheat. The identified transcription factor genes are excellent candidates for genetic engineering of stress tolerance in wheat.

Fine mapping of KLW1 that conditions kernel weight mainly through regulating kernel length in wheat (Triticum aestivum L.)

KLW1 was localized to a 0.6 cM interval near the centromere of chromosome 4B and found to be dominant in conditioning longer kernels and higher kernel weight. Kernel weight is a major wheat yield component and affected by kernel dimensions, filling process and kernel density. Because of this complexity, the mechanism underlying kernel weight is still far from clear. Qtgw.nau-4B or KLW1 was a major kernel weight QTL identified in the Nanda2419 × Wangshuibai population. We showed that introduction of the Nanda2419 allele into elite cultivar Wenmai6 resulted in longer kernels as well as higher kernel weight, without affecting other traits such as spike number per plant, plant height, spike length, spikelet number per spike, and kernel number per spike. KLW1 was dominant in conditioning higher kernel weight and functioned mainly through affecting kernel length. Using F₂ plants derived from KLW1 NIL, a high-density genetic map covering the QTL was constructed. KLW1 was consequently confined to the 0.6 cM Xwgrc4219-Xwgrc4067 interval by evaluating the recombinant lines in three field trials. KLW1 is complementary to KT1, the QTL on chromosome 5A of Nanda2419 for thicker and heavier kernels, in producing larger kernels with higher commercial value, augmenting its usefulness in wheat breeding.

Exploration of Piezo Channels in Bread Wheat (Triticum aestivum L.)

Piezo channels belong to an important class of cell membrane-bound, Ca2+-permeable, mechanosensitive channels consisting of a pore and multiple transmembrane helices. In plants, the functional aspects of Piezo channels have been less studied than other mechanosensitive channels. However, a few studies that have been carried out indicate the involvement of Piezo channels in stress response and developmental processes. In our analysis, we identified a total of three Piezo genes in the Triticum aestivum genome. The phylogenetic analysis revealed the monocot and dicot-specific clustering of Piezo proteins. The gene and protein structure analysis indicated their conserved architecture. The promoter region of each of the three Piezo genes contained light-, growth-and development-, hormone-, and stress-responsive cis-regulatory elements. Moreover, the differential expression of Piezo genes in tissue developmental stages and under abiotic and biotic stress conditions indicated their probable role in plant growth and development and various stresses. The quantitative real-time polymerase chain reaction (qRT-PCR) analysis suggested that TaPiezo1-D might be involved in Ca2+ homeostasis. In addition, protein–protein interaction indicated their precise role in glucose, hormone and stress responses. The miRNA interaction analysis further suggested their participation in signaling cascades and biological processes. The present study will extend our understanding about Piezo channels in Ca2+ mediated signaling in plants under various stresses and provide a path for the functional validation of TaPiezo genes in future research.

A computational study on the structure-function relationships of plant caleosins

Plant cells store energy in oil bodies constructed by structural proteins such as oleosins and caleosins. Although oil bodies usually accumulate in the seed and pollen of plants, caleosins are present in various organs and organelles. This issue, coupled with the diverse activities of caleosins, complicates the description of these oleo-proteins. Therefore, the current article proposes a new classification based on the bioinformatics analysis of the transmembrane topology of caleosins. Accordingly, the non-membrane class are the most abundant and diverse caleosins, especially in lower plants. Comparing the results with other reports suggests a stress response capacity for these caleosins. However, other classes play a more specific role in germination and pollination. A phylogenetic study also revealed two main clades that were significantly different in terms of caleosin type, expression profile, molecular weight, and isoelectric point (P < 0.01). In addition to the biochemical significance of the findings, predicting the structure of caleosins is necessary for constructing oil bodies used in the food and pharmaceutical industries.

Genome-wide identification and expression analysis of calmodulin and calmodulin-like genes in wheat (Triticum aestivum L.)

Calmodulin (CaM) and calmodulin-like (CML) genes are widely involved in plant growth and development and mediating plant stress tolerance. However, the whole genome scale studies about CaM and CML gene families have not been done in wheat, and the possible functions of most wheat CaM/CML gene members are still unknown. In this study, a total of 18 TaCaM and 230 TaCML gene members were identified in wheat genome. Among these genes, 28 TaCaM/CML gene members have 74 duplicated copies, while 21 genes have 48 transcript variants, resulting in 321 putative TaCaM/CML transcripts totally. Phylogenetic tree analysis showed that they can be classified into 7 subfamilies. Similar gene structures and protein domains can be found in members of the same gene cluster. The TaCaM/CML genes were spread among all 21 chromosomes with unbalanced distributions, while most of the gene clusters contained 3 homoeologous genes located in the same homoeologous chromosome group. Synteny analysis showed that most of TaCaM/CMLs gene members can be found with 1-4 paralogous genes in T. turgidum and Ae. Tauschii. High numbers of cis-acting elements related to plant hormones and stress responses can be observed in the promoters of TaCaM/CMLs. The spatiotemporal expression patterns showed that most of the TaCaM/TaCML genes can be detected in at least one tissue. The expression levels of TaCML17, 21, 30, 50, 59 and 75 in the root or shoot can be up-regulated by abiotic stresses, suggesting that TaCML17, 21, 30, 50, 59 and 75 may be related with responses to abiotic stresses in wheat. The spatiotemporal expression patterns of TaCaM/CML genes indicated they may be involved widely in wheat growth and development. Our results provide important clues for exploring functions of TaCaMs/CMLs in growth and development as well as responses to abiotic stresses in wheat in the future.

Wheat transcriptomic responses to extended feeding by wheat curl mites

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

OSCA Genes in Bread Wheat: Molecular Characterization, Expression Profiling, and Interaction Analyses Indicated Their Diverse Roles during Development and Stress Response

The hyperosmolality-gated calcium-permeable channels (OSCA) are pore-forming transmembrane proteins that function as osmosensors during various plant developmental processes and stress responses. In our analysis, through in silico approaches, a total of 42 OSCA genes are identified in the Triticum aestivum genome. A phylogenetic analysis reveals the close clustering of the OSCA proteins of Arabidopsis thaliana, Oryza sativa, and T. aestivum in all the clades, suggesting their origin before the divergence of dicots and monocots. Furthermore, evolutionary analyses suggest the role of segmental and tandem duplication events (Des) and purifying selection pressure in the expansion of the OSCA gene family in T. aestivum. Expression profiling in various tissue developmental stages and under abiotic and biotic stress treatments reveals the probable functioning of OSCA genes in plant development and the stress response in T. aestivum. In addition, protein-protein and protein-chemical interactions reveal that OSCA proteins might play a putative role in Ca²⁺-mediated developmental processes and adaptive responses. The miRNA interaction analysis strengthens the evidence for their functioning in various biological processes and stress-induced signaling cascades. The current study could provide a foundation for the functional characterization of TaOSCA genes in future studies.

Location and Identification on Chromosome 3B of Bread Wheat of Genes Affecting Chiasma Number

Understanding meiotic crossover (CO) variation in crops like bread wheat (Triticum aestivum L.) is necessary as COs are essential to create new, original and powerful combinations of genes for traits of agronomical interest. We cytogenetically characterized a set of wheat aneuploid lines missing part or all of chromosome 3B to identify the most influential regions for chiasma formation located on this chromosome. We showed that deletion of the short arm did not change the total number of chiasmata genome-wide, whereas this latter was reduced by ~35% while deleting the long arm. Contrary to what was hypothesized in a previous study, deletion of the long arm does not disturb the initiation of the synaptonemal complex (SC) in early meiotic stages. However, progression of the SC is abnormal, and we never observed its completion when the long arm is deleted. By studying six different deletion lines (missing different parts of the long arm), we revealed that at least two genes located in both the proximal (C-3BL2-0.22) and distal (3BL7-0.63-1.00) deletion bins are involved in the control of chiasmata, each deletion reducing the number of chiasmata by ~15%. We combined sequence analyses of deletion bins with RNA-Seq data derived from meiotic tissues and identified a set of genes for which at least the homoeologous copy on chromosome 3B is expressed and which are involved in DNA processing. Among these genes, eight (CAP-E1/E2, DUO1, MLH1, MPK4, MUS81, RTEL1, SYN4, ZIP4) are known to be involved in the recombination pathway.

Insight into the Roles of Proline-Rich Extensin-like Receptor Protein Kinases of Bread Wheat (Triticum aestivum L.)

Proline-rich extensin-like receptor protein kinases (PERKs) are known for their roles in the developmental processes and stress responses of many plants. We have identified 30 TaPERK genes in the genome of T. aestivum, exploring their evolutionary and syntenic relationship and analyzing their gene and protein structures, various cis-regulatory elements, expression profiling, and interacting miRNAs. The TaPERK genes formed 12 homeologous groups and clustered into four phylogenetic clades. All the proteins exhibited a typical domain organization of PERK and consisted of conserved proline residue repeats and serine-proline and proline-serine repeats. Further, the tyrosine-x-tyrosine (YXY) motif was also found conserved in thirteen TaPERKs. The cis-regulatory elements and expression profiling under tissue developmental stages suggested their role in plant growth processes. Further, the differential expression of certain TaPERK genes under biotic and abiotic stress conditions suggested their involvement in defense responses as well. The interaction of TaPERK genes with different miRNAs further strengthened evidence for their diverse biological roles. In this study, a comprehensive analysis of obtained TaPERK genes was performed, enriching our knowledge of TaPERK genes and providing a foundation for further possible functional analyses in future studies.

Characterization and Expression of the Pirin Gene Family in Triticum aestivum

Pirins are nuclear bicupin proteins, encoded by genes that are one of several gene families that comprise the Cupin superfamily in plants. Pirin genes have been implicated in stress response pathways studied in Arabidopsis and At-Pirin1 has been shown to interact with the heterotrimeric G-protein alpha subunit (GPA1). The aim of this study was to identify the members of the Pirin gene family in Triticum aestivum, to correct their annotations in the whole genome and gain an insight into their tissue-specific expression as well as their response to abiotic and biotic stresses. The Pirin gene family in T. aestivum is comprised of 18 genes that represent six paralogous gene copies, each having an A, B and D homeolog. Expression analysis of the Pirin genes in T. aestivum Illumina RNA-seq libraries, which included sampling from differing tissue types as well as abiotic and biotic stresses, indicates that the members of the Pirin gene family have specialized expression and play a role in stress responses. Pirin gene families are also identified in other monocots including Aegilops tauschii, Hordeum vulgare, Brachypodium distachyon, Oryza sativa, Zea mays, Sorghum bicolor and the dicot Arabidopsis thaliana.

The Use and Limitations of Exome Capture to Detect Novel Variation in the Hexaploid Wheat Genome

The bread wheat (Triticum aestivum) pangenome is a patchwork of variable regions, including translocations and introgressions from progenitors and wild relatives. Although a large number of these have been documented, it is likely that many more remain unknown. To map these variable regions and make them more traceable in breeding programs, wheat accessions need to be genotyped or sequenced. The wheat genome is large and complex and consequently, sequencing efforts are often targeted through exome capture. In this study, we employed exome capture prior to sequencing 12 wheat varieties; 10 elite T. aestivum cultivars and two T. aestivum landrace accessions. Sequence coverage across chromosomes was greater toward distal regions of chromosome arms and lower in centromeric regions, reflecting the capture probe distribution which itself is determined by the known telomere to centromere gene gradient. Superimposed on this general pattern, numerous drops in sequence coverage were observed. Several of these corresponded with reported introgressions. Other drops in coverage could not be readily explained and may point to introgressions that have not, to date, been documented.

Characterization of the heterotrimeric G protein gene families in Triticum aestivum and related species

This study characterizes the heterotrimeric G protein gene families in Triticum aestivum, their tissue-specific expression patterns during development and in response to biotic and abiotic stress conditions. There are three Gα genes, three Gβ and 12 Gγ genes, totaling 18 genes encoding heterotrimeric G proteins in the hexaploid wheat genome. Each haploid genome of the hexaploid T. aestivum has a single gene encoding the α subunit of the heterotrimeric G protein complex, GA1, a single Gβ and four Gγ genes. Each gene has three homeologous copies in the A, B and D genomes. The physical interaction between the Gβ (Gpb) and two Gγ subunits, Gpg1 and Gpg2, was shown through bimolecular fluorescence complementation (BiFC). The gene expression in response to biotic and abiotic stresses showed both up-regulation and down-regulation of members of the gene families. Gγ2-B and Gγ2-D are significantly upregulated during heat stress, GA1-D is upregulated by cold stress and Gγ1-A and Gγ1-D were upregulated by Fusarium graminearum inoculation in a F. graminearum resistant cultivar. This suggests that these members may play roles in biotic and abiotic signaling pathways and the roles of these genes within these pathways need further investigation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03156-9.

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.

Exploration of glutathione reductase for abiotic stress response in bread wheat (Triticum aestivum L.)

A total of seven glutathione reductase (GR) genes were identified in Triticum aestivum, which were used for comparative structural characterization, phylogenetic analysis and expression profiling with the GR genes of other cereal plants. The modulated gene expression and enzyme activity revealed the role of GRs in abiotic stress response in T. aestivum. Glutathione reductase (GR) is an enzymatic antioxidant that converts oxidized glutathione (GSSG) into reduced glutathione (GSH) through the ascorbate-glutathione cycle. In this study, a total of seven GR genes forming two homeologous groups were identified in the allohexaploid genome of bread wheat (Triticum aestivum). Besides, we identified three GR genes in each Aegilops tauschii, Brachypodium distachyon, Triticum urartu and Sorghum bicolor, which were used for comparative characterization. Phylogenetic analysis revealed the clustering of GR proteins into two groups; class I and class II, which were predicted to be localized in cytoplasm and chloroplast, respectively. The exon-intron and conserved motif patterns were almost conserved in each group, in which a maximum of 10 and 17 exons were present in chloroplastic and cytoplasmic GRs, respectively. The protein structure analysis confirmed the occurrence of conserved pyridine nucleotide disulfide oxidoreductase (Pyr_redox) and pyridine nucleotide disulfide oxidoreductase dimerization (Pyr_redox_dim) domains in each GR. The active site of GR proteins consisted of two conserved cysteine residues separated by four amino acid residues. Promoter analysis revealed the occurrence of growth and stress-related cis-active elements. Tissue-specific expression profiling suggested the involvement of GRs in both vegetative and reproductive tissue development in various plants. The differential expression of TaGR genes and enhanced GR enzyme activity suggested their roles under drought, heat, salt and arsenic stress. Interaction of GRs with other proteins and chemical compounds of the ascorbate-glutathione cycle revealed their coordinated functioning. The current study will provide a foundation for the validation of the precise role of each GR gene in future studies.

Investigation of Roles of TaTALE Genes during Development and Stress Response in Bread Wheat

The three amino acid loop extension (TALE) genes of the homeobox superfamily are responsible for numerous biological functions in plants. Herein, we identified a total of 72 TaTALE genes in the allohexaploid genome of bread wheat (Triticum aestivum L.) and performed a comprehensive investigation for gene and protein structural properties, phylogeny, expression patterns, and multilevel gene regulations. The identified TaTALE proteins were further classified into two groups, TaBLHs and TaKNOXs, which were tightly clustered into the phylogeny. The negative Ka/Ks ratio of duplicated genes suggested purifying selection pressure with confined functional divergence. Various signature domains and motifs were found conserved in both groups of proteins. The occurrence of diverse cis-regulatory elements and modulated expression during various developmental stages and in the presence of abiotic (heat, drought, salt) and two different fungal stresses suggested their roles in development and stress response, as well. The interaction of TaTALEs with the miRNAs and other development-related homeobox proteins also suggested their roles in growth and development and stress response. The present study revealed several important aspects of TaTALEs that will be useful in further functional validation of these genes in future studies.

Uncovering genomic regions controlling plant architectural traits in hexaploid wheat using different GWAS models

Wheat is a major food crop worldwide. The plant architecture is a complex trait mostly influenced by plant height, tiller number, and leaf morphology. Plant height plays a crucial role in lodging and thus affects yield and grain quality. In this study, a wheat population was genotyped by using Illumina iSelect 90K single nucleotide polymorphism (SNP) assay and finally 22,905 high-quality SNPs were used to perform a genome-wide association study (GWAS) for plant architectural traits employing four multi-locus GWAS (ML-GWAS) and three single-locus GWAS (SL-GWAS) models. As a result, 174 and 97 significant SNPs controlling plant architectural traits were detected by ML-GWAS and SL-GWAS methods, respectively. Among these SNP makers, 43 SNPs were consistently detected, including seven across multiple environments and 36 across multiple methods. Interestingly, five SNPs (Kukri_c34553_89, RAC875_c8121_1490, wsnp_Ex_rep_c66315_64480362, Ku_c5191_340, and tplb0049a09_1302) consistently detected across multiple environments and methods, played a role in modulating both plant height and flag leaf length. Furthermore, candidate SNPs (BS00068592_51, Kukri_c4750_452 and BS00022127_51) constantly repeated in different years and methods associated with flag leaf width and number of tillers. We also detected several SNPs (Jagger_c6772_80, RAC875_c8121_1490, BS00089954_51, Excalibur_01167_1207, and Ku_c5191_340) having common associations with more than one trait across multiple environments. By further appraising these GWAS methods, the pLARmEB and FarmCPU models outperformed in SNP detection compared to the other ML-GWAS and SL-GWAS methods, respectively. Totally, 152 candidate genes were found to be likely involved in plant growth and development. These finding will be helpful for better understanding of the genetic mechanism of architectural traits in wheat.

New insights into homoeologous copy number variations in the hexaploid wheat genome

Bread wheat is an allohexaploid species originating from two successive and recent rounds of hybridization between three diploid species that were very similar in terms of chromosome number, genome size, TE content, gene content and synteny. As a result, it has long been considered that most of the genes were in three pairs of homoeologous copies. However, these so-called triads represent only one half of wheat genes, while the remaining half belong to homoeologous groups with various number of copies across subgenomes. In this study, we examined and compared the distribution, conservation, function, expression and epigenetic profiles of triads with homoeologous groups having undergone a deletion (dyads) or a duplication (tetrads) in one subgenome. We show that dyads and tetrads are mostly located in distal regions and have lower expression level and breadth than triads. Moreover, they are enriched in functions related to adaptation and more associated with the repressive H3K27me3 modification. Altogether, these results suggest that triads mainly correspond to housekeeping genes and are part of the core genome, while dyads and tetrads belong to the Triticeae dispensable genome. In addition, by comparing the different categories of dyads and tetrads, we hypothesize that, unlike most of the allopolyploid species, subgenome dominance and biased fractionation are absent in hexaploid wheat. Differences observed between the three subgenomes are more likely related to two successive and ongoing waves of post-polyploid diploidization, that had impacted A and B more significantly than D, as a result of the evolutionary history of hexaploid wheat.

Molecular characterization revealed the role of catalases under abiotic and arsenic stress in bread wheat (Triticum aestivum L.)

Catalases are crucial antioxidant enzymes that reduce the excessive level of H₂O₂ caused by various environmental stresses and metal toxicity and hence protect the plant cells. In this study, a total of ten TaCAT genes, forming three homeologous groups, were identified in the genome of bread wheat (Triticum aestivum L.) and named as per the wheat gene symbolization guidelines. The identified catalases were characterized for various structural and physicochemical features. The proximal active-site (F(D/A)RERIPERVVHAKGASA) and heme-ligand (R(I/V)F(S/A)Y(A/S)DTQ) signature motifs, catalytic residues and peroxisomal targeting peptides were found conserved. Phylogenetic analysis clustered TaCATs into three classes, which showed conserved functional specialization based on their tissue specific expression. Modulated spatio-temporal expression of various TaCAT genes and alteration in total catalase enzyme activity during heat, drought, salt and arsenic (As^III and As^V) treatment suggested their roles in abiotic stress response and arsenic tolerance. Molecular cloning and overexpression of TaCAT3-B gene in Escherichia coli showed tolerance against heat, drought, salt and varied concentrations of As^III and As^V treatments. The results further confirmed their role in stress tolerance and recommended that these genes can be used in future stress management strategies for the development of abiotic and arsenic stress resistant transgenic crops.

Ph2 encodes the mismatch repair protein MSH7-3D that inhibits wheat homoeologous recombination

Meiotic recombination is a critical process for plant breeding, as it creates novel allele combinations that can be exploited for crop improvement. In wheat, a complex allohexaploid that has a diploid-like behaviour, meiotic recombination between homoeologous or alien chromosomes is suppressed through the action of several loci. Here, we report positional cloning of Pairing homoeologous 2 (Ph2) and functional validation of the wheat DNA mismatch repair protein MSH7-3D as a key inhibitor of homoeologous recombination, thus solving a half-century-old question. Similar to ph2 mutant phenotype, we show that mutating MSH7-3D induces a substantial increase in homoeologous recombination (up to 5.5 fold) in wheat-wild relative hybrids, which is also associated with a reduction in homologous recombination. These data reveal a role for MSH7-3D in meiotic stabilisation of allopolyploidy and provides an opportunity to improve wheat's genetic diversity through alien gene introgression, a major bottleneck facing crop improvement.

Molecular Characterization Revealed the Role of Thaumatin-Like Proteins of Bread Wheat in Stress Response

Thaumatin-like proteins (TLPs) are related to pathogenesis-related-5 (PR-5) family and involved in stress response. Herein, a total of 93 TLP genes were identified in the genome of Triticum aestivum. Further, we identified 26, 27, 39, and 37 TLP genes in the Brachypodium distachyon, Oryza sativa, Sorghum bicolor, and Zea mays genomes for comparative characterization, respectively. They could be grouped into small and long TLPs with conserved thaumatin signature motif. Tightly clustered genes exhibited conserved gene and protein structure. The physicochemical analyses suggested significant differences between small and long TLPs. Evolutionary analyses suggested the role of duplication events and purifying selection in the expansion of the TLP gene family. Expression analyses revealed the possible roles of TLPs in plant development and abiotic and fungal stress response. Recombinant expression of TaTLP2-B in Saccharomyces cerevisiae provided significant tolerance against cold, heat, osmotic, and salt stresses. The results depicted the importance of TLPs in cereal crops that would be highly useful in future crop improvement programs.

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.

Genome-wide identification, characterization, and expression patterns analysis of the SBP-box gene family in wheat (Triticum aestivum L.)

SQUAMOSA promoter-binding protein (SBP)-box genes encode a family of plant-specific transcription factors that play roles in plant growth and development. The characteristics of SBP-box genes in rice (Oryza sativa) and Arabidopsis have been reported, but their potential roles in wheat (Triticum aestivum) are not fully understood. In this study, 48 SBP-box genes (TaSBPs) were identified; they were located in all wheat chromosomes except for 4B and 4D. Six TaSBPs were identified as tandem duplication genes that formed three tandem duplication pairs, while 22 were segmentally duplicated genes that formed 16 segmental duplication pairs. Subcellular localization prediction showed TaSBPs were located in nucleus. Among the 48 TaSBPs, 24 were predicted to be putative targets of TamiR156. Phylogenetic analysis showed that TaSBPs, AtSBPs, and OsSBPs that shared similar functions were clustered into the same subgroups. The phylogenetic relationships between the TaSBPs were supported by the identification of highly conserved motifs and gene structures. Four types of cis-elements––transcription-related, development-related, hormone-related, and abiotic stress-related elements––were found in the TaSBP promoters. Expression profiles indicated most TaSBPs participate in flower development and abiotic stress responses. This study establishes a foundation for further investigation of TaSBP genes and provides novel insights into their biological functions.

Genome-Wide Development and Validation of Cost-Effective KASP Marker Assays for Genetic Dissection of Heat Stress Tolerance in Maize

Maize is the third most important cereal crop worldwide. However, its production is vulnerable to heat stress, which is expected to become more and more severe in coming years. Germplasm resilient to heat stress has been identified, but its underlying genetic basis remains poorly understood. Genomic mapping technologies can fill the void, provided robust markers are available to tease apart the genotype-phenotype relationship. In the present investigation, we used data from an RNA-seq experiment to identify single nucleotide polymorphisms (SNPs) between two contrasting lines, LM11 and CML25, sensitive and tolerant to heat stress, respectively. The libraries for RNA-seq were made following heat stress treatment from three separate tissues/organs, comprising the top leaf, ovule, and pollen, all of which are highly vulnerable to damage by heat stress. The single nucleotide variants (SNVs) calling used STAR mapper and GATK caller pipelines in a combined approach to identify highly accurate SNPs between the two lines. A total of 554,423, 410,698, and 596,868 SNVs were discovered between LM11 and CML25 after comparing the transcript sequence reads from the leaf, pollen, and ovule libraries, respectively. Hundreds of these SNPs were then selected to develop into genome-wide Kompetitive Allele-Specific PCR (KASP) markers, which were validated to be robust with a successful SNP conversion rate of 71%. Subsequently, these KASP markers were used to effectively genotype an F₂ mapping population derived from a cross of LM11 and CML25. Being highly cost-effective, these KASP markers provide a reliable molecular marker toolkit to not only facilitate the genetic dissection of the trait of heat stress tolerance but also to accelerate the breeding of heat-resilient maize by marker-assisted selection (MAS).

Genome-wide characterization and expression analysis suggested diverse functions of the mechanosensitive channel of small conductance-like (MSL) genes in cereal crops

Mechanosensitive ion channels are pore-forming transmembrane proteins that allow ions to move down their electrochemical gradient in response to mechanical stimuli. They participate in many plant developmental processes including the maintenance of plastid shape, pollen tube growth, etc. Herein, a total of 11, 10, 6, 30, 9, and 8 MSL genes were identified in Aegilops tauschii, Hordeum vulgare, Sorghum bicolor, Triticum aestivum, Triticum urartu, and Zea mays, respectively. These genes were located on various chromosomes of their respective cereal, while MSLs of T. urartu were found on scaffolds. The phylogenetic analysis, subcellular localization, and sequence homology suggested clustering of MSLs into two classes. These genes consisted of cis-regulatory elements related to growth and development, responsive to light, hormone, and stress. Differential expression of various MSL genes in tissue developmental stages and stress conditions revealed their precise role in development and stress responses. Altered expression during CaCl2 stress suggested their role in Ca2+ homeostasis and signaling. The co-expression analysis suggested their interactions with other genes involved in growth, defense responses etc. A comparative expression profiling of paralogous genes revealed either retention of function or pseudo-functionalization. The present study unfolded various characteristics of MSLs in cereals, which will facilitate their in-depth functional characterization in future studies.

Structural Variations Affecting Genes and Transposable Elements of Chromosome 3B in Wheats

Structural variations (SVs) such as copy number and presence–absence variations are polymorphisms that are known to impact genome composition at the species level and are associated with phenotypic variations. In the absence of a reference genome sequence, their study has long been hampered in wheat. The recent production of new wheat genomic resources has led to a paradigm shift, making possible to investigate the extent of SVs among cultivated and wild accessions. We assessed SVs affecting genes and transposable elements (TEs) in a Triticeae diversity panel of 45 accessions from seven tetraploid and hexaploid species using high-coverage shotgun sequencing of sorted chromosome 3B DNA and dedicated bioinformatics approaches. We showed that 23% of the genes are variable within this panel, and we also identified 330 genes absent from the reference accession Chinese Spring. In addition, 60% of the TE-derived reference markers were absent in at least one accession, revealing a high level of intraspecific and interspecific variability affecting the TE space. Chromosome extremities are the regions where we observed most of the variability, confirming previous hypotheses made when comparing wheat with the other grasses. This study provides deeper insights into the genomic variability affecting the complex Triticeae genomes at the intraspecific and interspecific levels and suggests a phylogeny with independent hybridization events leading to different hexaploid species.

Molecular characterization of ascorbate peroxidase (APX) and APX-related (APX-R) genes in Triticum aestivum L

Ascorbate peroxidases (APXs) are heme-dependent H₂O₂ scavenging enzymes involved in myriad biological processes. Herein, a total of 21 TaAPX and six TaAPX-R genes were identified from the A, B and D sub-genomes of Triticum aestivum L. The occurrence of three paralogous gene pairs with unequal evolutionary rate suggested functional divergence. The phylogenetic analysis formed four distinct clades having conserved gene and protein architecture, and sub-cellular localization. The tertiary structure analysis revealed the presence of helices and coils and residues involved in ligand binding. Transcriptional profiling of each TaAPX and TaAPX-R gene suggested their specific role during development and stress response. Modulated transcript expression and APX enzyme activity during various stress conditions indicated their role in stress response. Interaction analyses suggested their association with other genes, miRNAs and various legends. The present study reported numerous features of these genes, and may provide a platform for their detailed functional characterization in future studies.

Thaumatin-like protein kinases: Molecular characterization and transcriptional profiling in five cereal crops

Thaumatin-like protein kinases (TLPKs) are defense related proteins having antimicrobial property. Herein, we identified two TLPKs in the genome of Brachypodium distachyon and Oryza sativa, four in Hordeum vulgare and Sorghum bicolor, and 16 in Triticum aestivum. All the TLPKs were located at only one chromosome in each plant except T. aestivum, where they were located on chromosome 2 and chromosome 3. Paralogous analysis suggested the occurrence of one duplication event (DE) in each B. distachyon and O. sativa, two in H. vulgare while four DEs in T. aestivum genome during the evolution of TLPKs. The majority of TLPKs were intron less, while a few contains one or two introns. The introns were found in each 0, 1 and 2 phase. Protein structure analysis suggested the occurrence of a thaumatin and a kinase domain with a transmembrane (TM) helix in each TLPK. Further, a thaumatin family signature motif "GX[GF]XCXT[GA]DCX(1,2)GX(2,3)C", a "REDDD" motif and 16 cysteine residues were found conserved in the majority of TLPKs. Expression analysis indicated variable expression of TLPKs in various tissues of different cereal crops. They were high expressing in reproductive tissues in B. distachyon, while in leaves in T. aestivum. Modulated expression of TaTLPKs in the presence of fungal pathogen, and heat, drought and salt stress in T. aestivum suggested their roles in stress response. Co-expression analysis showed interaction of TLPKs with various development and stress related genes. The results indicated diverse roles of TLPKs, which can be utilized for the development of eco-friendly pest resistant crops in future.

Thaumatin-like protein kinases: Molecular characterization and transcriptional profiling in five cereal crops

Thaumatin-like protein kinases (TLPKs) are defense related proteins having antimicrobial property. Herein, we identified two TLPKs in the genome of Brachypodium distachyon and Oryza sativa, four in Hordeum vulgare and Sorghum bicolor, and 16 in Triticum aestivum. All the TLPKs were located at only one chromosome in each plant except T. aestivum, where they were located on chromosome 2 and chromosome 3. Paralogous analysis suggested the occurrence of one duplication event (DE) in each B. distachyon and O. sativa, two in H. vulgare while four DEs in T. aestivum genome during the evolution of TLPKs. The majority of TLPKs were intron less, while a few contains one or two introns. The introns were found in each 0, 1 and 2 phase. Protein structure analysis suggested the occurrence of a thaumatin and a kinase domain with a transmembrane (TM) helix in each TLPK. Further, a thaumatin family signature motif "GX[GF]XCXT[GA]DCX(1,2)GX(2,3)C", a "REDDD" motif and 16 cysteine residues were found conserved in the majority of TLPKs. Expression analysis indicated variable expression of TLPKs in various tissues of different cereal crops. They were high expressing in reproductive tissues in B. distachyon, while in leaves in T. aestivum. Modulated expression of TaTLPKs in the presence of fungal pathogen, and heat, drought and salt stress in T. aestivum suggested their roles in stress response. Co-expression analysis showed interaction of TLPKs with various development and stress related genes. The results indicated diverse roles of TLPKs, which can be utilized for the development of eco-friendly pest resistant crops in future.

De Novo Transcriptome Analysis of Durum Wheat Flag Leaves Provides New Insights Into the Regulatory Response to Elevated CO2 and High Temperature

Global warming is becoming a significant problem for food security, particularly in the Mediterranean basin. The use of molecular techniques to study gene-level responses to environmental changes in non-model organisms is increasing and may help to improve the mechanistic understanding of durum wheat response to elevated CO2 and high temperature. With this purpose, we performed transcriptome RNA sequencing (RNA-Seq) analyses combined with physiological and biochemical studies in the flag leaf of plants grown in field chambers at ear emergence. Enhanced photosynthesis by elevated CO2 was accompanied by an increase in biomass and starch and fructan content, and a decrease in N compounds, as chlorophyll, soluble proteins, and Rubisco content, in association with a decline of nitrate reductase and initial and total Rubisco activities. While high temperature led to a decline of chlorophyll, Rubisco activity, and protein content, the glucose content increased and starch decreased. Furthermore, elevated CO2 induced several genes involved in mitochondrial electron transport, a few genes for photosynthesis and fructan synthesis, and most of the genes involved in secondary metabolism and gibberellin and jasmonate metabolism, whereas those related to light harvesting, N assimilation, and other hormone pathways were repressed. High temperature repressed genes for C, energy, N, lipid, secondary, and hormone metabolisms. Under the combined increases in atmospheric CO2 and temperature, the transcript profile resembled that previously reported for high temperature, although elevated CO2 partly alleviated the downregulation of primary and secondary metabolism genes. The results suggest that there was a reprogramming of primary and secondary metabolism under the future climatic scenario, leading to coordinated regulation of C-N metabolism towards C-rich metabolites at elevated CO2 and a shift away from C-rich secondary metabolites at high temperature. Several candidate genes differentially expressed were identified, including protein kinases, receptor kinases, and transcription factors.

Characterizing the Role of TaWRKY13 in Salt Tolerance

The WRKY transcription factor superfamily is known to participate in plant growth and stress response. However, the role of this family in wheat (Triticum aestivum L.) is largely unknown. Here, a salt-induced gene TaWRKY13 was identified in an RNA-Seq data set from salt-treated wheat. The results of RT-qPCR analysis showed that TaWRKY13 was significantly induced in NaCl-treated wheat and reached an expression level of about 22-fold of the untreated wheat. Then, a further functional identification was performed in both Arabidopsis thaliana and Oryza sativa L. Subcellular localization analysis indicated that TaWRKY13 is a nuclear-localized protein. Moreover, various stress-related regulatory elements were predicted in the promoter. Expression pattern analysis revealed that TaWRKY13 can also be induced by polyethylene glycol (PEG), exogenous abscisic acid (ABA), and cold stress. After NaCl treatment, overexpressed Arabidopsis lines of TaWRKY13 have a longer root and a larger root surface area than the control (Columbia-0). Furthermore, TaWRKY13 overexpression rice lines exhibited salt tolerance compared with the control, as evidenced by increased proline (Pro) and decreased malondialdehyde (MDA) contents under salt treatment. The roots of overexpression lines were also more developed. These results demonstrate that TaWRKY13 plays a positive role in salt stress.

Engineering meiotic recombination pathways in rice

In the last 15 years, outstanding progress has been made in understanding the function of meiotic genes in the model dicot and monocot plants Arabidopsis and rice (Oryza sativa L.), respectively. This knowledge allowed to modulate meiotic recombination in Arabidopsis and, more recently, in rice. For instance, the overall frequency of crossovers (COs) has been stimulated 2.3- and 3.2-fold through the inactivation of the rice FANCM and RECQ4 DNA helicases, respectively, two genes involved in the repair of DNA double-strand breaks (DSBs) as noncrossovers (NCOs) of the Class II crossover pathway. Differently, the programmed induction of DSBs and COs at desired sites is currently explored by guiding the SPO11-1 topoisomerase-like transesterase, initiating meiotic recombination in all eukaryotes, to specific target regions of the rice genome. Furthermore, the inactivation of 3 meiosis-specific genes, namely PAIR1, OsREC8 and OsOSD1, in the Mitosis instead of Meiosis (MiMe) mutant turned rice meiosis into mitosis, thereby abolishing recombination and achieving the first component of apomixis, apomeiosis. The successful translation of Arabidopsis results into a crop further allowed the implementation of two breakthrough strategies that triggered parthenogenesis from the MiMe unreduced clonal egg cell and completed the second component of diplosporous apomixis. Here, we review the most recent advances in and future prospects of the manipulation of meiotic recombination in rice and potentially other major crops, all essential for global food security.

A wheat NAC interacts with an orphan protein and enhances resistance to Fusarium head blight disease

Taxonomically-restricted orphan genes play an important role in environmental adaptation, as recently demonstrated by the fact that the Pooideae-specific orphan TaFROG (Triticum aestivum Fusarium Resistance Orphan Gene) enhanced wheat resistance to the economically devastating Fusarium head blight (FHB) disease. Like most orphan genes, little is known about the cellular function of the encoded protein TaFROG, other than it interacts with the central stress regulator TaSnRK1α. Here, we functionally characterized a wheat (T. aestivum) NAC-like transcription factor TaNACL-D1 that interacts with TaFROG and investigated its' role in FHB using studies to assess motif analyses, yeast transactivation, protein-protein interaction, gene expression and the disease response of wheat lines overexpressing TaNACL-D1. TaNACL-D1 is a Poaceae-divergent NAC transcription factor that encodes a Triticeae-specific protein C-terminal region with transcriptional activity and a nuclear localisation signal. The TaNACL-D1/TaFROG interaction was detected in yeast and confirmed in planta, within the nucleus. Analysis of multi-protein interactions indicated that TaFROG could form simultaneously distinct protein complexes with TaNACL-D1 and TaSnRK1α in planta. TaNACL-D1 and TaFROG are co-expressed as an early response to both the causal fungal agent of FHB, Fusarium graminearum and its virulence factor deoxynivalenol (DON). Wheat lines overexpressing TaNACL-D1 were more resistant to FHB disease than wild type plants. Thus, we conclude that the orphan protein TaFROG interacts with TaNACL-D1, a NAC transcription factor that forms part of the disease response evolved within the Triticeae.

Identification, characterization and expression analysis of lineage-specific genes within Triticeae

Lineage-specific genes (LSGs) are a set of genes in a given taxon without significant sequence similarity to genes and intergenic sequences of other taxa and are functional. The tribe Triticeae mainly includes species of different ploidy levels, such as staple food crops wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.). This study is aimed at mining and characterizing the Triticeae-specific genes (TSGs) using expressed sequence data of wheat. A total of 3812 TSGs was identified and they were generally characterized by smaller size, fewer exons, shorter open reading frames and lower expression levels. Most TSGs were expressed with tissue preference and many of them were predominantly expressed in reproduction related tissues, especially in young stamen. Nearly one third of the TSGs were stress-responsive and inducible under abiotic and/or biotic stresses. A co-expression-based annotation supported the relevance of some TSGs with reproduction and stress responses, indicating their potential economic importance.

Genomic dissection and transcriptional profiling of Cysteine-rich receptor-like kinases in five cereals and functional characterization of TaCRK68-A

Cysteine-rich receptor-like kinases (CRK) constitute one of the largest subfamily of receptor-like kinases, which play crucial roles in plant development and stress response. In total, 43, 37, 36, 38 and 170 CRK genes including duplicated genes were identified in the genome of Brachypodium distachyon, Hordeum vulgare, Oryza sativa, Sorghum bicolor and Triticum aestivum, respectively. These CRK proteins were tightly clustered into four phylogenetic groups and exhibited close syntenic relationship among orthologous genes. Majority of CRK proteins contain a transmembrane domain for plasma membrane localization. The organization of exon/intron, domains and motifs were variably conserved. Tissue-specific expression suggested the involvement of certain CRK genes in plant development. Modulated expression revealed their specific stress-responsive functions. Co-expression and interaction analysis indicated their role in signaling. Ks value and divergence time analysis suggested duplication of TaCRK genes before the hybridization of T. aestivum sub-genomes. Expression comparison of duplicated TaCRK genes revealed functional retention, neofunctionalization or pseudo-functionalization. Recombinant expression of a stress-responsive gene TaCRK68-A in Escherichia coli and Saccharomyces cerevisiae displayed enhanced tolerance against heat, drought, cold and salinity stresses. The study suggested vital functions of CRKs during development and stresses, and provides the basis for functional characterization of each gene in future studies.

Genomic tools for durum wheat breeding: de novo assembly of Svevo transcriptome and SNP discovery in elite germplasm

BACKGROUND

The tetraploid durum wheat (Triticum turgidum L. ssp. durum Desf. Husnot) is an important crop which provides the raw material for pasta production and a valuable source of genetic diversity for breeding hexaploid wheat (Triticum aestivum L.). Future breeding efforts to enhance yield potential and climate resilience will increasingly rely on genomics-based approaches to identify and select beneficial alleles. A deeper characterisation of the molecular and functional diversity of the durum wheat transcriptome will be instrumental to more effectively harness its genetic diversity.

RESULTS

We report on the de novo transcriptome assembly of durum wheat cultivar 'Svevo'. The transcriptome of four tissues/organs (shoots and roots at the seedling stage, reproductive organs and developing grains) was assembled de novo, yielding 180,108 contigs, with a N50 length of 1121 bp and mean contig length of 883 bp. Alignment against the transcriptome of nine plant species identified 43% of transcripts with homology to at least one reference transcriptome. The functional annotation was completed by means of a combination of complementary software. The presence of differential expression between the A- and B-homoeolog copies of the durum wheat tetraploid genome was ascertained by phase reconstruction of polymorphic sites based on the T. urartu transcripts and inferring homoeolog-specific sequences. We observed greater expression divergence between A and B homoeologs in grains rather than in leaves and roots. The transcriptomes of 13 durum wheat cultivars spanning the breeding period from 1969 to 2005 were analysed for SNP diversity, leading to 95,358 non-rare, hemi-SNPs shared among two or more cultivars and 33,747 locus-specific (diploid inheritance) SNPs.

CONCLUSIONS

Our study updates and expands the de novo transcriptome reference assembly available for durum wheat. Out of 180,108 assembled transcripts, 13,636 were specific to the Svevo cultivar as compared to the only other reference transcriptome available for durum, thus contributing to the identification of the tetraploid wheat pan-transcriptome. Additionally, the analysis of 13 historically relevant hallmark varieties produced a SNP dataset that could successfully validate the genotyping in tetraploid wheat and provide a valuable resource for genomics-assisted breeding of both tetraploid and hexaploid wheats.

Identification, characterization and expression profiling of cation-proton antiporter superfamily in Triticum aestivum L. and functional analysis of TaNHX4-B

The monovalent cation proton antiporter (CPA) superfamily comprises Na⁺/H⁺ exchanger (NHX), K⁺ efflux antiporter (KEA), and cation/H⁺ exchanger (CHX) family proteins, which play vital functions in plants. A total of 107 TaCPA proteins were identified in Triticum aestivum, and phylogenetically classified into 35 TaNHX, 24 TaKEA and 48 TaCHX proteins. These families had representatives derived from all three sub-genomes. TaKEA genes consisted of higher number of exons, followed by TaNHXs and TaCHXs. The occurrence of about 10 transmembrane regions and higher composition of helices and coils support their membrane-bound and hydrophobic nature. Diverse expression in various tissues and modulated expression under stress conditions suggested their role in development and in response to stress. Co-expression analyses revealed their complex interaction networks. Expression of TaNHX4-B.1 and TaNHX4-B.4 facilitated differential abiotic stress tolerance to Escherichia coli. Our study provides comprehensive information about CPA genes, which would be useful in their future functional characterization.

Transcriptome sequencing assisted discovery and computational analysis of novel SNPs associated with flowering in Raphanus sativus in-bred lines for marker-assisted backcross breeding

The sequencing of radish genome aids in the better understanding and tailoring of traits associated with economic importance. In order to accelerate the genomics assisted breeding and genetic selection, transcriptomes of 33 radish inbred lines with diverse traits were sequenced for the development of single nucleotide polymorphic (SNP) markers. The sequence reads ranged from 2,560,543,741 bp to 20,039,688,139 bp with the GC (%) of 47.80-49.34 and phred quality score (Q30) of 96.47-97.54%. A total of 4951 polymorphic SNPs were identified among the accessions after stringent filtering and 298 SNPs with efficient marker assisted backcross breeding (MAB) markers were generated from the polymorphic SNPs. Further, functional annotations of SNPs revealed the effects and importance of the SNPs identified in the flowering process. The SNPs were predominantly associated with the four major flowering related transcription factors such as MYB, MADS box (AG), AP2/EREB, and bHLH. In addition, SNPs in the vital flowering integrator gene (FT) and floral repressors (EMBRYONIC FLOWER 1, 2, and FRIGIDA) were identified among the radish inbred lines. Further, 50 SNPs were randomly selected from 298 SNPs and validated using Kompetitive Allele Specific PCR genotyping system (KASP) in 102 radish inbred lines. The homozygosity of the inbred lines varied from 56 to 96% and the phylogenetic analysis resulted in the clustering of inbred lines into three subgroups. Taken together, the SNP markers identified in the present study can be utilized for the discrimination, seed purity test, and adjusting parental combinations for breeding in radish.

Characterization of the Esi3/RCI2/PMP3 gene family in the Triticeae

Background

Members of the Early Salt Induced 3 (Esi3/RCI2/PMP3) gene family in plants have been shown to be induced in response to both biotic and abiotic stresses and to enhance stress tolerance in both transgenic plants and Saccharomyces cerevisiae. Esi3 was first identified as a salt stress induced gene in the salt tolerant wild wheat grass, Lophopyrum elongatum, and subsequently homologous genes in many other species were found to be members of the gene family. These include Arabidopsis thaliana and Oryza sativa where they are referred to as Rare Cold Inducible 2 (RCI2), and Zea mays where they are referred to as Plasma Membrane Protein 3 (PMP3). This study characterizes the Esi3 family members in Triticum aestivum and explores the tissue specific expression patterns of the gene family members as well as their response to a variety of environmental stresses.

Results

The Esi3 gene family was found to have a total of 29 family members comprised of ten paralogous groups in the hexaploid T. aestivum. Each paralogous group contains three homeologous copies, one in each of the A, B and D genomes with the exception of Esi3–2 which is missing the B copy. The genes of the Esi3 gene family were also identified in four other monocot species, Aegilops tauschii, Hordeum vulgare, Secale cereale and Sorghum bicolor, and were confirmed or corrected for Brachypodium distachyon, Oryza sativa and Zea mays, as well as the dicot Arabidopsis thaliana. Gene expression of the Esi3s was analyzed using tissue-specific, abiotic and biotic stress RNA-Seq 454 sequence libraries and Affymetrix microarray data for T. aestivum.

Conclusions

Members of nearly all paralogous groups of the Esi3 genes in T. aestivum have altered gene expression in response to abiotic or biotic stress conditions. In addition, there are modest differences in gene expression among homeologous members of the gene family. This suggests that the Esi3 gene family plays an important role in the plants response to the stresses presented in this study. The Esi3–9 in T. aestivum has a unique N terminal extension placing it into Group III, a new group for the Esi3/RCI2/PMP3 gene family.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5311-8) contains supplementary material, which is available to authorized users.

Comparative functional genomics analysis of bHLH gene family in rice, maize and wheat

BACKGROUND

The basic helix-loop-helix transcription factors play important roles in diverse cellular and molecular processes. Comparative functional genomics can provide powerful approaches to draw inferences about gene function and evolution among species. The comprehensive comparison of bHLH gene family in different gramineous plants has not yet been reported.

RESULTS

In this study, a total of 183, 231 and 571 bHLHs were identified in rice, maize and wheat genomes respectively, and 1154 bHLH genes from the three species and Arabidopsis were classified into 36 subfamilies. Of the identified genes, 110 OsbHLHs, 188 ZmbHLHs and 209 TabHLHs with relatively high mRNA abundances were detected in one or more tissues during development, and some of them exhibited tissue-specific expression such as TabHLH454-459, ZmbHLH099-101 and OsbHLH037 in root, TabHLH559-562, - 046, - 047 and ZmbHLH010, - 072, - 226 in leaf, TabHLH216-221, - 333, - 335, - 340 and OsbHLH005, - 141 in inflorescence, TabHLH081, ZmbHLH139 and OsbHLH144 in seed. Forty five, twenty nine and thirty one differentially expressed bHLHs were respectively detected in wheat, maize and rice under drought stresses using RNA-seq technology. Among them, the expressions of TabHLH046, - 047, ZmbHLH097, - 098, OsbHLH006 and - 185 were strongly induced, whereas TabHLH303, - 562, ZmbHLH155, - 154, OsbHLH152 and - 113 showed significant down-regulation. Twenty two TabHLHs were induced after stripe rust infection at 24 h and nine of them were suppressed at 72 hpi, whereas 28 and 6 TabHLHs exhibited obviously down- and up-regulation after powdery mildew attack respectively. Forty one ZmbHLHs were differentially expressed in response to F. verticillioides infection. Twenty two co-expression modules were identified by the WGCNA, some of which were associated with particular tissue types. And GO enrichment analysis for the modules showed that some TabHLHs were involved in the control of several biological processes, such as tapetal PCD, lipid metabolism, iron absorption, stress responses and signal regulation.

CONCLUSION

The present study identifies the bHLH family in rice, maize and wheat genomes, and detailedly discusses the evolutionary relationships, expression and function of bHLHs. This study provides some novel and detail information about bHLHs, and may facilitate understanding the molecular basis of the plant growth, development and stress physiology.

De novo assembly of wheat root transcriptomes and transcriptional signature of longitudinal differentiation

Hidden underground, root systems constitute an important part of the plant for its development, nourishment and sensing the soil environment around it, but we know very little about its genetic regulation in crop plants like wheat. In the present study, we de novo assembled the root transcriptomes in reference cultivar Chinese Spring from RNA-seq reads generated by the 454-GS-FLX and HiSeq platforms. The FLX reads were assembled into 24,986 transcripts with completeness of 54.84%, and the HiSeq reads were assembled into 91,543 high-confidence protein-coding transcripts, 2,404 low-confidence protein-coding transcripts, and 13,181 non-coding transcripts with the completeness of >90%. Combining the FLX and HiSeq assemblies, we assembled a root transcriptome of 92,335 ORF-containing transcripts. Approximately 7% of the coding transcripts and ~2% non-coding transcripts are not present in the current wheat genome assembly. Functional annotation of both assemblies showed similar gene ontology patterns and that ~7% coding and >5% non-coding transcripts are root-specific. Transcription quantification identified 1,728 differentially expressed transcripts between root tips and maturation zone, and functional annotation of these transcripts captured a transcriptional signature of longitudinal development of wheat root. With the transcriptomic resources developed, this study provided the first view of wheat root transcriptome under different developmental zones and laid a foundation for molecular studies of wheat root development and growth using a reverse genetic approach.

Optical and physical mapping with local finishing enables megabase-scale resolution of agronomically important regions in the wheat genome

BACKGROUND

Numerous scaffold-level sequences for wheat are now being released and, in this context, we report on a strategy for improving the overall assembly to a level comparable to that of the human genome.

RESULTS

Using chromosome 7A of wheat as a model, sequence-finished megabase-scale sections of this chromosome were established by combining a new independent assembly using a bacterial artificial chromosome (BAC)-based physical map, BAC pool paired-end sequencing, chromosome-arm-specific mate-pair sequencing and Bionano optical mapping with the International Wheat Genome Sequencing Consortium RefSeq v1.0 sequence and its underlying raw data. The combined assembly results in 18 super-scaffolds across the chromosome. The value of finished genome regions is demonstrated for two approximately 2.5 Mb regions associated with yield and the grain quality phenotype of fructan carbohydrate grain levels. In addition, the 50 Mb centromere region analysis incorporates cytological data highlighting the importance of non-sequence data in the assembly of this complex genome region.

CONCLUSIONS

Sufficient genome sequence information is shown to now be available for the wheat community to produce sequence-finished releases of each chromosome of the reference genome. The high-level completion identified that an array of seven fructosyl transferase genes underpins grain quality and that yield attributes are affected by five F-box-only-protein-ubiquitin ligase domain and four root-specific lipid transfer domain genes. The completed sequence also includes the centromere.

Shifting the limits in wheat research and breeding using a fully annotated reference genome

An annotated reference sequence representing the hexaploid bread wheat genome in 21 pseudomolecules has been analyzed to identify the distribution and genomic context of coding and noncoding elements across the A, B, and D subgenomes. With an estimated coverage of 94% of the genome and containing 107,891 high-confidence gene models, this assembly enabled the discovery of tissue- and developmental stage-related coexpression networks by providing a transcriptome atlas representing major stages of wheat development. Dynamics of complex gene families involved in environmental adaptation and end-use quality were revealed at subgenome resolution and contextualized to known agronomic single-gene or quantitative trait loci. This community resource establishes the foundation for accelerating wheat research and application through improved understanding of wheat biology and genomics-assisted breeding.