Genome-wide transcriptome study in wheat identified candidate genes related to processing quality, majority of them showing interaction (quality x development) and having temporal and spatial distributions

Temporal Gene Expression Profiles From Pollination to Seed Maturity in Sorghum Provide Core Candidates for Engineering Seed Traits

Sorghum (Sorghum bicolor (L.) Moench) is a highly nutritional multipurpose millet crop. However, the genetic and molecular regulatory mechanisms governing sorghum grain development and the associated agronomic traits remain unexplored. In this study, we performed a comprehensive transcriptomic analysis of pistils collected 1-2 days before pollination, and developing seeds collected -2, 10, 20 and 30 days after pollination of S. bicolor variety M35-1. Out of 31 337 genes expressed in these stages, 12 804 were differentially expressed in the consecutive stages of seed development. These exhibited 10 dominant expression patterns correlated with the distinct pathways and gene functions. Functional analysis, based on the pathway mapping, transcription factor enrichment and orthology, delineated the key patterns associated with pollination, fertilization, early seed development, grain filling and seed maturation. Furthermore, colocalization with previously reported quantitative trait loci (QTLs) for grain weight/size revealed 48 differentially expressed genes mapping to these QTL regions. Comprehensive literature mining integrated with QTL mapping and expression data shortlisted 25, 17 and 8 core candidates for engineering grain size, starch and protein content, respectively.

RNA-Seq transcriptome profiling of immature grain wheat is a technique for understanding comparative modeling of baking quality

Improving the baking quality is a primary challenge in the wheat flour production value chain, as baking quality represents a crucial factor in determining its overall value. In the present study, we conducted a comparative RNA-Seq analysis on the high baking quality mutant "O-64.1.10" genotype and its low baking quality wild type "Omid" cultivar to recognize potential genes associated with bread quality. The cDNA libraries were constructed from immature grains that were 15 days post-anthesis, with an average of 16.24 and 18.97 million paired-end short-read sequences in the mutant and wild-type, respectively. A total number of 733 transcripts with differential expression were identified, 585 genes up-regulated and 188 genes down-regulated in the "O-64.1.10" genotype compared to the "Omid". In addition, the families of HSF, bZIP, C2C2-Dof, B3-ARF, BES1, C3H, GRF, HB-HD-ZIP, PLATZ, MADS-MIKC, GARP-G2-like, NAC, OFP and TUB were appeared as the key transcription factors with specific expression in the "O-64.1.10" genotype. At the same time, pathways related to baking quality were identified through Kyoto Encyclopedia of Genes and Genomes. Collectively, we found that the endoplasmic network, metabolic pathways, secondary metabolite biosynthesis, hormone signaling pathway, B group vitamins, protein pathways, pathways associated with carbohydrate and fat metabolism, as well as the biosynthesis and metabolism of various amino acids, have a great deal of potential to play a significant role in the baking quality. Ultimately, the RNA-seq results were confirmed using quantitative Reverse Transcription PCR for some hub genes such as alpha-gliadin, low molecular weight glutenin subunit and terpene synthase (gibberellin) and as a resource for future study, 127 EST-SSR primers were generated using RNA-seq data.

Comparative transcriptome analyses revealed key genes involved in high amylopectin biosynthesis in wheat

High amylopectin starch is an important modified starch for food processing industries. Despite a thorough understanding of starch biosynthesis pathway, the regulatory mechanism responsible for amylopectin biosynthesis is not well explored. The present study utilized transcriptome sequencing approach to understand the molecular basis of high amylopectin content in three high amylopectin mutant wheat lines ('TAC 6', 'TAC 358', and 'TAC 846') along with parent variety 'C 306'. Differential scanning calorimetry (DSC) of high amylopectin starch identified a high thermal transition temperature and scanning electron microscopy (SEM) revealed more spherical starch granules in mutant lines compared to parent variety. A set of 4455 differentially expressed genes (DEGs) were identified at two-fold compared to the parent variety in high amylopectin wheat mutants. At ten-fold, 279 genes, including two starch branching genes (SBEIIa and SBEIIb), were up-regulated and only 30 genes, including the starch debranching enzyme (DBE), were down-regulated. Among the genes, different isoforms of sucrose non-fermenting-1-related protein kinase-1 (TaSnRK1α2-3B and TaSnRK1α2-3D) and its regulatory subunit, sucrose non-fermenting-4 (SNF-4-2A, SNF-4-2B, and SNF-4-5D), were found to be highly up-regulated. Further, expression of the DEGs related to starch biosynthesis pathway and TaSnRK1α2 and SNF-4 was performed using qRT-PCR. High expression of TaSnRK1α2, SNF-4, and SBEII isoforms suggests their probable role in high amylopectin starch biosynthesis in grain endosperm.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03364-3.

Screening of Induced Mutants Led to the Identification of Starch Biosynthetic Genes Associated with Improved Resistant Starch in Wheat

Several health benefits are obtained from resistant starch, also known as healthy starch. Enhancing resistant starch with genetic modification has huge commercial importance. The variation of resistant starch content is narrow in wheat, in relation to which limited improvement has been attained. Hence, there is a need to produce a wheat population that has a wide range of variations in resistant starch content. In the present study, stable mutants were screened that showed significant variation in the resistant starch content. A megazyme kit was used for measuring the resistant starch content, digestible starch, and total starch. The analysis of variance showed a significant difference in the mutant population for resistant starch. Furthermore, four diverse mutant lines for resistant starch content were used to study the quantitative expression patterns of 21 starch metabolic pathway genes; and to evaluate the candidate genes for resistant starch biosynthesis. The expression pattern of 21 starch metabolic pathway genes in two diverse mutant lines showed a higher expression of key genes regulating resistant starch biosynthesis (GBSSI and their isoforms) in the high resistant starch mutant lines, in comparison to the parent variety (J411). The expression of SBEs genes was higher in the low resistant starch mutants. The other three candidate genes showed overexpression (BMY, Pho1, Pho2) and four had reduced (SSIII, SBEI, SBEIII, ISA3) expression in high resistant starch mutants. The overexpression of AMY and ISA1 in the high resistant starch mutant line JE0146 may be due to missense mutations in these genes. Similarly, there was a stop_gained mutation for PHO2; it also showed overexpression. In addition, the gene expression analysis of 21 starch metabolizing genes in four different mutants (low and high resistant starch mutants) shows that in addition to the important genes, several other genes (phosphorylase, isoamylases) may be involved and contribute to the biosynthesis of resistant starch. There is a need to do further study about these new genes, which are responsible for the fluctuation of resistant starch in the mutants.

Effects of the different waxy proteins on starch biosynthesis, starch physicochemical properties and Chinese noodle quality in wheat

Noodles are an important food in Asia. Wheat starch is the most important component in Chinese noodles. Loss of the waxy genes leads to lower activity of starch synthesis enzymes and decreased amylose content that further affects starch properties and noodle quality. To study the effects of different waxy (Wx) protein subunits on starch biosynthesis and processing quality, the high-yielding wheat cultivar Jimai 22 was treated with the mutagen ethyl methane sulfonate (EMS) to produce a population of Wx lines and chosen 7 Wx protein combinations. The amylose content increased but swelling power decreased as the number of Wx proteins increased. Both GBSS activity and gene expression were the lowest for the waxy mutant, followed by the mutants with 1 Wx protein. The combinations of these mutant alleles lead to reductions in both RNA expression and protein levels. Noodles made from materials with 2 Wx protein subunits had the highest score, which agreed with peak viscosity. The influence of the Wx-B1 protein on amylose synthesis and noodle quality was the highest, whereas the influence of Wx-A1 protein was the lowest. Mutants with lower amylose content caused by the absence of 1 subunit, especially the Wx-B1 subunit, had superior noodle quality. Additionally, the identified mutant lines can be used as intermediate materials to improve wheat quality.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01292-x.

Goyal A, Chaudhury A. Molecular and Transcriptional Regulation of Seed Development in Cereals: Present Status and Future Prospects

Cereals are a rich source of vitamins, minerals, carbohydrates, fats, oils and protein, making them the world’s most important source of nutrition. The influence of rising global population, as well as the emergence and spread of disease, has the major impact on cereal production. To meet the demand, there is a pressing need to increase cereal production. Optimal seed development is a key agronomical trait that contributes to crop yield. The seed development and maturation is a complex process that includes not only embryo and endosperm development, but also accompanied by huge physiological, biochemical, metabolic, molecular and transcriptional changes. This chapter discusses the growth of cereal seed and highlights the novel biological insights, with a focus on transgenic and new molecular breeding, as well as biotechnological intervention strategies that have improved crop yield in two major cereal crops, primarily wheat and rice, over the last 21 years (2000–2021).

Decoding the genome of superior chapatti quality Indian wheat variety 'C 306' unravelled novel genomic variants for chapatti and nutrition quality related genes

An Indian wheat variety, 'C 306' has good chapatti quality, which is controlled by multiple genes that have not been explored. We report the high quality de novo assembled genome of 'C 306' by combining short and long read sequencing data. The hybrid assembly covered 93% of gene space and identified about 142 K coding genes, 34% repetitive DNA and ~ 501 K SSR motifs. The phylogenetic analysis of about 83 K orthologous protein groups suggested the closest relationship with T. turgidum, T. aestivum and Ae. tauschii. Genome wide analysis annotated 69,217,536 genomic variants. Out of them, 1423 missense and 117 deleterious variants identified in processing, nutrition, and chapatti quality related genes such as alpha- and beta-gliadin, SSI, SSIII, SUT1, SBEI, CHS, YSL, DMAS, and NAS encoded proteins. These variants may affect quality genes. The genomic data will be potential genomic resources in wheat breeding programs for quality improvement.

A null allele of granule bound starch synthase (Wx-B1) may be one of the major genes controlling chapatti softness

Chapatti (unleavened flatbread) is a staple food in northern India and neighboring countries but the genetics behind its processing quality are poorly understood. To understand the genes determining chapatti quality, differentially expressed genes were selected from microarray data of contrasting chapatti cultivars. From the gene and trait association studies, a null allele of granule bound starch synthase (GBSS; Wx-B1) was found to be associated with low amylose content and good chapatti quality. For validation, near-isogenic lines (NILs) of this allele were created by marker assisted backcross (MAB) breeding. Background screening indicated 88.2 to 96.7% background recovery in 16 selected BC3F5 NILs. Processing quality and sensory evaluation of selected NILs indicated improvement in chapatti making quality. Traits that showed improvement were mouthfeel, tearing strength and softness indicating that the Wx-B1 may be one of the major genes controlling chapatti softness.

Rising Atmospheric Temperature Impact on Wheat and Thermotolerance Strategies

Temperature across the globe is increasing continuously at the rate of 0.15–0.17 °C per decade since the industrial revolution. It is influencing agricultural crop productivity. Therefore, thermotolerance strategies are needed to have sustainability in crop yield under higher temperature. However, improving thermotolerance in the crop is a challenging task for crop scientists. Therefore, this review work was conducted with the aim of providing information on the wheat response in three research areas, i.e., physiology, breeding, and advances in genetics, which could assist the researchers in improving thermotolerance. The optimum temperature for wheat growth at the heading, anthesis, and grain filling duration is 16 ± 2.3 °C, 23 ± 1.75 °C, and 26 ± 1.53 °C, respectively. The high temperature adversely influences the crop phenology, growth, and development. The pre-anthesis high temperature retards the pollen viability, seed formation, and embryo development. The post-anthesis high temperature declines the starch granules accumulation, stem reserve carbohydrates, and translocation of photosynthates into grains. A high temperature above 40 °C inhibits the photosynthesis by damaging the photosystem-II, electron transport chain, and photosystem-I. Our review work highlighted that genotypes which can maintain a higher accumulation of proline, glycine betaine, expression of heat shock proteins, stay green and antioxidant enzymes activity viz., catalase, peroxidase, super oxide dismutase, and glutathione reductase can tolerate high temperature efficiently through sustaining cellular physiology. Similarly, the pre-anthesis acclimation with heat treatment, inorganic fertilizer such as nitrogen, potassium nitrate and potassium chloride, mulches with rice husk, early sowing, presoaking of a 6.6 mM solution of thiourea, foliar application of 50 ppm dithiothreitol, 10 mg per kg of silicon at heading and zinc ameliorate the crop against the high temperature. Finally, it has been suggested that modern genomics and omics techniques should be used to develop thermotolerance in wheat.

Early Root Transcriptomic Changes in Wheat Seedlings Colonized by Trichoderma harzianum Under Different Inorganic Nitrogen Supplies

Wheat is one of the most important crops worldwide. The use of plant growth promoting microorganisms, such as those of the genus Trichoderma, constitutes an alternative to chemical fertilizers, since they are cheaper and are not detrimental to the environment. However, the interaction between Trichoderma and wheat plants has been scarcely studied, at least at a molecular level. In the present work, a microarray approach was used to study the early transcriptomic changes induced in wheat roots by Trichoderma harzianum, applied alone or in combination with different concentrations of calcium nitrate [Ca(NO₃)₂], which was last used as nitrogen (N) source. Our results show that T. harzianum causes larger transcriptomic changes than Ca(NO₃)₂ in wheat roots, and such changes are different when plants are challenged with Trichoderma alone or treated with a combination of T. harzianum and Ca(NO₃)₂. Overall, T. harzianum activates the expression of defense-related genes at early stages of the interaction with the roots, while this fungus reduces the expression of genes related to plant growth and development. Moreover, the current study in wheat roots, subjected to the different T. harzianum and Ca(NO₃)₂ combinations, reveals that the number of transcriptomic changes was higher when compared against those caused by the different Ca(NO₃)₂ concentrations than when it was compared against those caused by T. harzianum. N metabolism gene expression changes were in agreement with the levels of nitrate reductase activity measured in plants from Trichoderma plus Ca(NO₃)₂ conditions. Results were also concordant with plant phenotypes, which showed reduced growth at early interaction stages when inoculated with T. harzianum or with its combination with Ca(NO₃)₂ at the lowest dosage. These results were in a good agreement with the recognized role of Trichoderma as an inducer of plant defense.

Overlapping transcriptional expression response of wheat zinc-induced facilitator-like transporters emphasize important role during Fe and Zn stress

BACKGROUND

Hexaploid wheat is an important cereal crop that has been targeted to enhance grain micronutrient content including zinc (Zn) and iron (Fe). In this direction, modulating the expression of plant transporters involved in Fe and Zn homeostasis has proven to be one of the promising approaches. The present work was undertaken to identify wheat zinc-induced facilitator-like (ZIFL) family of transporters. The wheat ZIFL genes were characterized for their transcriptional expression response during micronutrient fluctuations and exposure to multiple heavy metals.

RESULTS

The genome-wide analyses resulted in identification of fifteen putative TaZIFL-like genes, which were distributed only on Chromosome 3, 4 and 5. Wheat ZIFL proteins subjected to the phylogenetic analysis showed the uniform distribution along with rice, Arabidopsis and maize. In-silico analysis of the promoters of the wheat ZIFL genes demonstrated the presence of multiple metal binding sites including those which are involved in Fe and heavy metal homeostasis. Quantitative real-time PCR analysis of wheat ZIFL genes suggested the differential regulation of the transcripts in both roots and shoots under Zn surplus and also during Fe deficiency. Specifically, in roots, TaZIFL2.3, TaZIFL4.1, TaZIFL4.2, TaZIFL5, TaZIFL6.1 and TaZIFL6.2 were significantly up-regulated by both Zn and Fe. This suggested that ZIFL could possibly be regulated by both the nutrient stress in a tissue specific manner. When exposed to heavy metals, TaZIFL4.2 and TaZIFL7.1 show significant up-regulation, whereas TaZIFL5 and TaZIFL6.2 remained almost unaffected.

CONCLUSION

This is the first report for detailed analysis of wheat ZIFL genes. ZIFL genes also encode for transporter of mugineic acid (TOM) proteins, that are involved in the release of phytosiderophores to enhance Fe/Zn uptake. The detailed expression analysis suggests the varying expression patterns during development of wheat seedlings and also against abiotic/biotic stresses. Overall, this study will lay foundation to prioritize functional assessment of the candidate ZIFL as a putative TOM protein in wheat.

"Expression of genes encoding protein disulfide isomerase (PDI) in cultivars and lines of common wheat with different baking quality of flour"

BACKGROUND

The subject of this research was to investigate the level of expression of genes encoding protein disulfide isomerase (PDI) in cultivars and lines of wheat with different baking value of flour. PDI plays a key role in the formation of disulfide bonds in newly formed proteins. Each of cultivars and lines had a specific set of high molecular weight glutenin subunits (HMW-GS). Based on the presence of individual subunits, the potential baking value is predicted. Sometimes this value is not confirmed during technological analysis. Since there are cases where flour has a better or worse value than expected on the basis of the genotype, the expression of PDI genes was considered as a potential cause for discrepancies mentioned.

RESULTS

Analysis focused on three stages of grain development. The expression level of PDI genes was compared between wheat cultivars and lines with different genotype-phenotype combinations, which means diversified sets of HMW-GS combined with diversified qualitative classification. The highest expression level of PDI was noticed at early stage of grain development, which is consistent with the function of PDI. The expression level was evaluated by the real-time PCR technique.

CONCLUSION

Results obtained in this work did not allow for a clear statement of decisive significance of PDI in the context of shaping the final baking value. The results of this work contribute to an ever more in-depth understanding of the mechanisms governing baking value, and thus to the progress of the selection of new varieties with more beneficial properties.

Wheat seed transcriptome reveals genes controlling key traits for human preference and crop adaptation

Analysis of the transcriptome of the developing wheat grain has associated expression of genes with traits involving production (e.g. yield) and quality (e.g. bread quality). Photosynthesis in the grain may be important in retaining carbon that would be lost in respiration during grain filling and may contribute to yield in the late stages of seed formation under warm and dry environments. A small number of genes have been identified as having been selected by humans to optimize the performance of wheat for foods such as bread. Genes determining flour yield in milling have been discovered. Hardness is explained by variations in expression of pin genes. Knowledge of these genes should dramatically improve the efficiency of breeding better climate adapted wheat genotypes.

Comparative analysis of transcriptome in two wheat genotypes with contrasting levels of drought tolerance

Drought tolerance is a complex trait that is governed by multiple genes. The study presents differential transcriptome analysis between drought-tolerant (Triticum aestivum Cv. C306) and drought-sensitive (Triticum aestivum Cv. WL711) genotypes, using Affymetrix GeneChip® Wheat Genome Array. Both genotypes exhibited diverse global transcriptional responses under control and drought conditions. Pathway analysis suggested significant induction or repression of genes involved in secondary metabolism, nucleic acid synthesis, protein synthesis, and transport in C306, as compared to WL711. Significant up- and downregulation of transcripts for enzymes, hormone metabolism, and stress response pathways were observed in C306 under drought. The elevated expression of plasma membrane intrinsic protein 1 and downregulation of late embryogenesis abundant in the leaf tissues could play an important role in delayed wilting in C306. The other regulatory genes such as MT, FT, AP2, SKP1, ABA2, ARF6, WRKY6, AOS, and LOX2 are involved in defense response in C306 genotype. Additionally, transcripts with unknown functions were identified as differentially expressed, which could participate in drought responses.

Development of EMS-induced mutation population for amylose and resistant starch variation in bread wheat (Triticum aestivum) and identification of candidate genes responsible for amylose variation

BACKGROUND

Starch is a major part of cereal grain. It comprises two glucose polymer fractions, amylose (AM) and amylopectin (AP), that make up about 25 and 75 % of total starch, respectively. The ratio of the two affects processing quality and digestibility of starch-based food products. Digestibility determines nutritional quality, as high amylose starch is considered a resistant or healthy starch (RS type 2) and is highly preferred for preventive measures against obesity and related health conditions. The topic of nutrition security is currently receiving much attention and consumer demand for food products with improved nutritional qualities has increased. In bread wheat (Triticum aestivum L.), variation in amylose content is narrow, hence its limited improvement. Therefore, it is necessary to produce wheat lines or populations showing wide variation in amylose/resistant starch content. In this study, a set of EMS-induced M4 mutant lines showing dynamic variation in amylose/resistant starch content were produced. Furthermore, two diverse mutant lines for amylose content were used to study quantitative expression patterns of 20 starch metabolic pathway genes and to identify candidate genes for amylose biosynthesis.

RESULTS

A population comprising 101 EMS-induced mutation lines (M4 generation) was produced in a bread wheat (Triticum aestivum) variety. Two methods of amylose measurement in grain starch showed variation in amylose content ranging from ~3 to 76 % in the population. The method of in vitro digestion showed variation in resistant starch content from 1 to 41 %. One-way ANOVA analysis showed significant variation (p < 0.05) in amylose and resistant starch content within the population. A multiple comparison test (Dunnett's test) showed that significant variation in amylose and resistant starch content, with respect to the parent, was observed in about 89 and 38 % of the mutant lines, respectively. Expression pattern analysis of 20 starch metabolic pathway genes in two diverse mutant lines (low and high amylose mutants) showed higher expression of key genes of amylose biosynthesis (GBSSI and their isoforms) in the high amylose mutant line, in comparison to the parent. Higher expression of amylopectin biosynthesis (SBE) was observed in the low amylose mutant lines. An additional six candidate genes showed over-expression (BMY, SPA) and reduced-expression (SSIII, SBEI, SBEIII, ISA3) in the high amylose mutant line, indicating that other starch metabolic genes may also contribute to amylose biosynthesis.

CONCLUSION

In this study a set of 101 EMS-induced mutant lines (M4 generation) showing variation in amylose and resistant starch content in seed were produced. This population serves as useful germplasm or pre-breeding material for genome-wide study and improvement of starch-based processing and nutrition quality in wheat. It is also useful for the study of the genetic and molecular basis of amylose/resistant starch variation in wheat. Furthermore, gene expression analysis of 20 starch metabolic genes in the two diverse mutant lines (low and high amylose mutants) indicates that in addition to key genes, several other genes (such as phosphorylases, isoamylases, and pullulanases) may also be involved in contributing to amylose/amylopectin biosynthesis.

Functional cereals for production in new and variable climates

Adaptation of cereal crops to variable or changing climates requires that essential quality attributes are maintained to deliver food that will be acceptable to human consumers. Advances in cereal genomics are delivering insights into the molecular basis of nutritional and functional quality traits in cereals and defining new genetic resources. Understanding the influence of the environment on expression of these traits will support the retention of these essential functional properties during climate adaptation. New cereals for use as whole grain or ground to flour for other food products may be based upon the traditional species such as rice and wheat currently used in these food applications but may also include new options exploiting genomics tools to allow accelerated domestication of new species.

Comparative Analysis of Phenolic Compound Characterization and Their Biosynthesis Genes between Two Diverse Bread Wheat (Triticum aestivum) Varieties Differing for Chapatti (Unleavened Flat Bread) Quality

Phenolic compounds (PCs) affect the bread quality and can also affect the other types of end-use food products such as chapatti (unleavened flat bread), now globally recognized wheat-based food product. The detailed analysis of PCs and their biosynthesis genes in diverse bread wheat (Triticum aestivum) varieties differing for chapatti quality have not been studied. In this study, the identification and quantification of PCs using UPLC-QTOF-MS and/or MS/MS and functional genomics techniques such as microarrays and qRT-PCR of their biosynthesis genes have been studied in a good chapatti variety, "C 306" and a poor chapatti variety, "Sonalika." About 80% (69/87) of plant phenolic compounds were tentatively identified in these varieties. Nine PCs (hinokinin, coutaric acid, fertaric acid, p-coumaroylqunic acid, kaempferide, isorhamnetin, epigallocatechin gallate, methyl isoorientin-2'-O-rhamnoside, and cyanidin-3-rutinoside) were identified only in the good chapatti variety and four PCs (tricin, apigenindin, quercetin-3-O-glucuronide, and myricetin-3-glucoside) in the poor chapatti variety. Therefore, about 20% of the identified PCs are unique to each other and may be "variety or genotype" specific PCs. Fourteen PCs used for quantification showed high variation between the varieties. The microarray data of 44 phenolic compound biosynthesis genes and 17 of them on qRT-PCR showed variation in expression level during seed development and majority of them showed low expression in the good chapatti variety. The expression pattern in the good chapatti variety was largely in agreement with that of phenolic compounds. The level of variation of 12 genes was high between the good and poor chapatti quality varieties and has potential in development of markers. The information generated in this study can be extended onto a larger germplasm set for development of molecular markers using QTL and/or association mapping approaches for their application in wheat breeding.

Transcriptome profiling of wheat glumes in wild emmer, hulled landraces and modern cultivars

Background

Wheat domestication is considered as one of the most important events in the development of human civilization. Wheat spikelets have undergone significant changes during evolution under domestication, resulting in soft glumes and larger kernels that are released easily upon threshing. Our main goal was to explore changes in transcriptome expression in glumes that accompanied wheat evolution under domestication.

Methods

A total of six tetraploid wheat accessions were selected for transcriptome profiling based on their rachis brittleness and glumes toughness. RNA pools from glumes of the central spikelet at heading time were used to construct cDNA libraries for sequencing. The trimmed reads from each library were separately aligned to the reference sub-genomes A and B, which were extracted from wheat survey sequence. Differentially expression analysis and functional annotation were performed between wild and domesticated wheat, to identity candidate genes associated with evolution under domestication. Selected candidate genes were validated using real time PCR.

Results

Transcriptome profiles of wild emmer wheat, wheat landraces, and wheat cultivars were compared using next generation sequencing (RNA-seq). We have found a total of 194,893 transcripts, of which 73,150 were shared between wild, landraces, and cultivars. From 781 differentially expressed genes (DEGs), 336 were down-regulated and 445 were up-regulated in the domesticated compared to wild wheat genotypes. Gene Ontology (GO) annotation assigned 293 DEGs (37.5 %) to GO term groups, of which 134 (17.1 %) were down-regulated and 159 (20.4 %) up-regulated in the domesticated wheat. Some of the down-regulated DEGs in domesticated wheat are related to the biosynthetic pathways that eventually define the mechanical strength of the glumes, such as cell wall, lignin, pectin and wax biosynthesis. The reduction in gene expression of such genes, may explain the softness of the glumes in the domesticated forms. In addition, we have identified genes involved in nutrient remobilization that may affect grain size and other agronomic traits evolved under domestication.

Conclusions

The comparison of RNA-seq profiles between glumes of wheat groups differing in glumes toughness and rachis brittleness revealed a few DEGs that may be involved in glumes toughness and nutrient remobilization. These genes may be involved in processes of wheat improvement under domestication.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1996-0) contains supplementary material, which is available to authorized users.

Genome-wide identification and expression characterization of ABCC-MRP transporters in hexaploid wheat

The ABCC multidrug resistance associated proteins (ABCC-MRP), a subclass of ABC transporters are involved in multiple physiological processes that include cellular homeostasis, metal detoxification, and transport of glutathione-conjugates. Although they are well-studied in humans, yeast, and Arabidopsis, limited efforts have been made to address their possible role in crop like wheat. In the present work, 18 wheat ABCC-MRP proteins were identified that showed the uniform distribution with sub-families from rice and Arabidopsis. Organ-specific quantitative expression analysis of wheat ABCC genes indicated significantly higher accumulation in roots (TaABCC2, TaABCC3, and TaABCC11 and TaABCC12), stem (TaABCC1), leaves (TaABCC16 and TaABCC17), flag leaf (TaABCC14 and TaABCC15), and seeds (TaABCC6, TaABCC8, TaABCC12, TaABCC13, and TaABCC17) implicating their role in the respective tissues. Differential transcript expression patterns were observed for TaABCC genes during grain maturation speculating their role during seed development. Hormone treatment experiments indicated that some of the ABCC genes could be transcriptionally regulated during seed development. In the presence of Cd or hydrogen peroxide, distinct molecular expression of wheat ABCC genes was observed in the wheat seedlings, suggesting their possible role during heavy metal generated oxidative stress. Functional characterization of the wheat transporter, TaABCC13 a homolog of maize LPA1 confirms its role in glutathione-mediated detoxification pathway and is able to utilize adenine biosynthetic intermediates as a substrate. This is the first comprehensive inventory of wheat ABCC-MRP gene subfamily.

Comparative transcriptional profiling of two wheat genotypes, with contrasting levels of minerals in grains, shows expression differences during grain filling

Wheat is one of the most important cereal crops in the world. To identify the candidate genes for mineral accumulation, it is important to examine differential transcriptome between wheat genotypes, with contrasting levels of minerals in grains. A transcriptional comparison of developing grains was carried out between two wheat genotypes- Triticum aestivum Cv. WL711 (low grain mineral), and T. aestivum L. IITR26 (high grain mineral), using Affymetrix GeneChip Wheat Genome Array. The study identified a total of 580 probe sets as differentially expressed (with log2 fold change of ≥2 at p≤0.01) between the two genotypes, during grain filling. Transcripts with significant differences in induction or repression between the two genotypes included genes related to metal homeostasis, metal tolerance, lignin and flavonoid biosynthesis, amino acid and protein transport, vacuolar-sorting receptor, aquaporins, and stress responses. Meta-analysis revealed spatial and temporal signatures of a majority of the differentially regulated transcripts.