RNA-seq reveals differentially expressed genes of rice (Oryza sativa) spikelet in response to temperature interacting with nitrogen at meiosis stage

Background: Rice (Oryza sativa) is one of the most important cereal crops, providing food for more than half of the world's population. However, grain yields are challenged by various abiotic stresses such as drought, fertilizer, heat, and their interaction. Rice at reproductive stage is much more sensitive to environmental temperatures, and little is known about molecular mechanisms of rice spikelet in response to high temperature interacting with nitrogen (N). Results: Here we reported the transcriptional profiling analysis of rice spikelet at meiosis stage using RNA sequencing (RNA-seq) as an attempt to gain insights into molecular events associated with temperature and nitrogen. This study received four treatments: 1) NN: normal nitrogen level (165 kg ha(-1)) with natural temperature (30 degrees C); 2) HH: high nitrogen level (330 kg ha-1) with high temperature (37 degrees C); 3) NH: normal nitrogen level and high temperature; and 4) HN: high nitrogen level and natural temperature, respectively. The de novo assembly generated 52,553,536 clean reads aligned with 72,667 unigenes. About 10 M reads were identified from each treatment. In these differentially expressed genes (DEGs), we found 151 and 323 temperature-responsive DEGs in NN-vs-NH and HN-vs-HH, and 114 DEGs were co-expressed. Meanwhile, 203 and 144 nitrogen-responsive DEGs were focused in NN-vs-HN and NH-vs-HH, and 111 DEGs were co-expressed. The temperature-responsive genes were principally associated with calcium-dependent protein, cytochrome, flavonoid, heat shock protein, peroxidase, ubiquitin, and transcription factor while the nitrogen-responsive genes were mainly involved in glutamine synthetase, transcription factor, anthocyanin, amino acid transporter, leucine zipper protein, and hormone. It is noted that, rice spikelet fertility was significantly decreased under high temperature, but it was more reduced under higher nitrogen. Accordingly, numerous spikelet genes involved in pollen development, pollen tube growth, pollen germination, especially sporopollenin biosynthetic process, and pollen exine formation were mainly down-regulated under high temperature. Moreover, the expression levels of co-expressed DEGs including 5 sporopollenin biosynthetic process and 7 pollen exine formation genes of NN-vs-NH were lower than that of HN-vs-HH. Therefore, these spikelet genes may play important roles in response to high temperature with high nitrogen and may be good candidates for crop improvement. Conclusions: This RNA-seq study will help elucidate the molecular mechanisms of rice spikelet defense response to high temperature interacting with high nitrogen level.