Search

Wheat is a staple food worldwide and provides 40% of the calories in the diet. Climate change and global warming pose a threat to wheat production, however, and demand a deeper understanding of how heat stress might impact wheat production and wheat biology. However, it is difficult to identify novel heat stress associated genes when the genomic information is not available. Wheat has a very large and complex genome that is about 37 times the size of the rice genome. The present study sequenced the whole transcriptome of the wheat cv. HD2329 at the flowering stage, under control (22° ± 3°C) and heat stress (42°C, 2 h) conditions using Illumina HiSeq and Roche GS-FLX 454 platforms. We assembled more than 26.3 and 25.6 million high-quality reads from the control and HS-treated tissues transcriptome sequences respectively. About 76,556 (control) and 54,033 (HS-treated) contigs were assembled and annotated de novo using different assemblers and a total of 21,529 unigenes were obtained. Gene expression profile showed significant differential expression of 1525 transcripts under heat stress, of which 27 transcripts showed very high (>10) fold upregulation. Cellular processes such as metabolic processes, protein phosphorylation, oxidations-reductions, among others were highly influenced by heat stress. In summary, these observations significantly enrich the transcript dataset of wheat available on public domain and show a de novo approach to discover the heat-responsive transcripts of wheat, which can accelerate the progress of wheat stress-genomics as well as the course of wheat breeding programs in the era of climate change.

Calcium acts as a signaling molecule and plays very important role in plants-tolerance to the abiotic stresses including the heat stress. We studied the effect of exogenous Ca2+ (10 mM) prior to the heat stress (42 °C, 2 h) on different molecular and biochemical parameters associated with thermotolerance in wheat (Triticum aestivum). We could observe 64 and 52 unique protein spots in HD2967 (thermotolerant) and HD2329 (thermosusceptible) cultivars under normal and heat stress condition. Similarly, 13 and 9 (HD2967) and 19 and 9 proteins (HD2329) were observed to be up-regulated and down-regulated in response to HS. MALDI-TOF-TOF/MS characterization identified the differentially expressed protein (DEP) spots as calcium dependent protein kinase (CDPK), oxygen evolving enhancer protein, HSP17, HSP70, Rubisco activase etc. Transcript profiling of identified stress-associated genes (SAGs) showed very high expression of CDPK, HSFA4a, HSP17, SOD and APX in response to Ca2+ + HS in HD2967 compared to HD2329 cultivars of wheat. Similarly accumulation of signaling molecules (H2O2 and CDPK) as well as osmolyte was observed maximum in response to Ca2+ + HS in HD2967 compared to HD2329 cultivars. Very high activities of guaiacol-peroxidase and ascorbate peroxidase were observed in Ca-treated HD2967 compared to HD2329 when exposed to HS. It is the abundance of these chaperones and antioxidant enzymes in thermotolerant cultivar which limit the accumulation of H2O2 in response to Ca2+ and HS. Exogenous Ca2+ application showed negative correlation with lipid peroxidation and positive correlation with total antioxidant capacity of the cell system under the elevated temperature. Alteration in Ca2+ efflux triggers the activities of kinases and peroxidases more in tolerant compared to susceptible cultivars which regulate the accumulation of ROS inside the cell and attenuate the lipid peroxidation process. Ultimately, it enhances the total antioxidant capacity and thermotolerance of the plants under the heat stress.

Global warming is a major threat for agriculture and food security, and in many cases the negative impacts are already apparent. Wheat is one of the most important staple food crops and is highly sensitive to the heat stress (HS) during reproductive and grain-filling stages. Here, whole transcriptome analysis of thermotolerant wheat cv. HD2985 was carried out at the post-anthesis stage under control (22 ± 3 °C) and HS-treated (42 °C, 2 h) conditions using Illumina Hiseq and Roche GS-FLX 454 platforms. We assembled ~24 million (control) and ~23 million (HS-treated) high-quality trimmed reads using different assemblers with optimal parameters. De novo assembly yielded 52,567 (control) and 59,658 (HS-treated) unigenes. We observed 785 transcripts to be upregulated and 431 transcripts to be downregulated under HS; 78 transcripts showed >10-fold upregulation such as HSPs, metabolic pathway-related genes, etc. Maximum number of upregulated genes was observed to be associated with processes such as HS-response, protein-folding, oxidation-reduction and photosynthesis. We identified 2008 and 2483 simple sequence repeats (SSRs) markers from control and HS-treated samples; 243 SSRs were observed to be overlying on stress-associated genes. Polymorphic study validated four SSRs to be heat-responsive in nature. Expression analysis of identified differentially expressed transcripts (DETs) showed very high fold increase in the expression of catalytic chaperones (HSP26, HSP17, and Rca) in contrasting wheat cvs. HD2985 and HD2329 under HS. We observed positive correlation between RNA-seq and qRT-PCR expression data. The present study culminated in greater understanding of the heat-response of tolerant genotype and has provided good candidate genes for the marker development and screening of wheat germplasm for thermotolerance.

Heat stress is one of the major problems in agriculturally important cereal crops, especially wheat. Here, we have constructed a subtracted cDNA library from the endosperm of HS-treated (42°C for 2 h) wheat cv. HD2985 by suppression subtractive hybridization (SSH). We identified ~550 recombinant clones ranging from 200 to 500 bp with an average size of 300 bp. Sanger’s sequencing was performed with 205 positive clones to generate the differentially expressed sequence tags (ESTs). Most of the ESTs were observed to be localized on the long arm of chromosome 2A and associated with heat stress tolerance and metabolic pathways. Identified ESTs were BLAST search using Ensemble, TriFLD and TIGR databases and the predicted CDS were translated and aligned with the protein sequences available in pfam and InterProScan 5 databases to predict the differentially expressed proteins (DEPs). We observed eight different types of post-translational modifications (PTMs) in the DEPs corresponds to the cloned ESTs—147 sites with phosphorylation, 21 sites with sumoylation, 237 with palmitoylation, 96 sites with S-nitrosylation, 3066 calpain cleavage sites, and 103 tyrosine nitration sites, predicted to sense the heat stress and regulate the expression of stress genes. Twelve DEPs were observed to have transmembrane helixes (TMH) in their structure, predicted to play the role of sensors of HS. Quantitative Real-Time PCR of randomly selected ESTs showed very high relative expression of HSP17 under HS; up-regulation was observed more in wheat cv. HD2985 (thermotolerant), as compared to HD2329 (thermosusceptible) during grain-filling. The abundance of transcripts was further validated through northern blot analysis. The ESTs and their corresponding DEPs can be used as molecular marker for screening or targeted precision breeding program. PTMs identified in the DEPs can be used to elucidate the thermotolerance mechanism of wheat – a novel step towards the development of ‘climate-smart’ wheat.

RuBisCo activase (Rca) is a catalytic chaperone involved in modulating the activity of RuBisCo (key enzyme of photosynthetic pathway). Here, we identified eight novel transcripts from wheat through data mining predicted to be Rca and cloned a transcript of 1.4 kb from cv. HD2985, named as TaRca1 (GenBank acc. no. KC776912). Single copy number of TaRca1 was observed in wheat genome. Expression analysis in diverse wheat genotypes (HD2985, Halna, PBW621 and HD2329) showed very high relative expression of TaRca1 in Halna under control and HS-treated, as compared to other cultivars at different stages of growth. TaRca1 protein was predicted to be chloroplast-localized with numerous potential phosphorylation sites. Nothern blot analysis showed maximum accumulation of TaRca1 transcript in thermotolerant cv. during mealy-ripe stage, as compared to thermosusceptible. Decrease in the photosynthetic parameters was observed in all the cultivars, except PBW621 in response to HS. We observed significant increase in the Rca activity in all the cultivars under HS at different stages of growth. HS causes decrease in the RuBisCo activity; maximum reduction was observed during pollination stage in thermosusceptible cvs. as validated through immunoblotting. We observed uniform carbon distribution in different tissues of thermotolerant cvs., as compared to thermosusceptible. Similarly, tolerance level of leaf was observed maximum in Halna having high Rca activity under HS. A positive correlation was observed between the transcript and activity of TaRca1 in HS-treated Halna. Similarly, TaRca1 enzyme showed positive correlation with the activity of RuBisCo. There is, however, need to manipulate the thermal stability of TaRca1 enzyme through protein engineering for sustaining the photosynthetic rate under HS – a novel approach towards development of ‘climate-smart’ crop.

MicroRNAs (miRNAs) are small endogenous RNAs of ~22 nucleotides that have been shown to play regulatory role by negatively affecting the expression of genes at the post-transcriptional level. Information of miRNAs on some important crops like soybean, Arabidopsis, and rice, etc. are available, but no study on heat-responsive novel miRNAs has yet been reported in wheat (Triticum aestivum L.). In the present investigation, a popular wheat cultivar HD2985 was used in small RNA library construction and Illumina HiSeq 2000 was used to perform high-throughput sequencing of the library after cluster generation; 110,896,604 and 87,743,861 reads were generated in the control (22 °C) and heat-treated (42 °C for 2 h) samples, respectively. Forty-four precursor and mature miRNAs were found in T. aestivum from miRBase v 19. The frequencies of the miRNA families varied from 2 (tae-miR1117) to 60,672 (tae-miR159b). We identify 1052 and 902 mature miRNA sequences in HD2985 control and HS-treated samples by mapping on reference draft genome of T. aestivum. Maximum identified miRNAs were located on IWGSC_CSS_3B_scaff (chromosome 3B). We could identify 53 and 46 mature miRNA in the control and HS samples and more than 516 target genes by mapping on the reference genome of Oryza sativa, Zea mays, and Sorghum bicolor. Using different pipelines and plant-specific criteria, 37 novel miRNAs were identified in the control and treated samples. Six novel miRNA were validated using qRT-PCR to be heat-responsive. A negative correlation was, however, observed between the expression of novel miRNAs and their targets. Target prediction and pathway analysis revealed their involvement in the heat stress tolerance. These novel miRNAs are new additions to miRNA database of wheat, and the regulatory network will be made use of in deciphering the mechanism of thermotolerance in wheat.