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.

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.

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.