Long non-coding RNAs (lncRNAs) have been found to play a vital role in several gene regulatory networks involved in the various biological processes in plants related to stress response. However, systematic analyses of lncRNAs expressed in rice Cadmium (Cd) stress are seldom studied. Thus, we presented the characterization and expression of lncRNAs in rice root development at an early stage in response to Cd stress. The lncRNA deep sequencing revealed differentially expressed lncRNAs among Cd stress and normal condition. In the Cd stress group, 69 lncRNAs were up-regulated and 75 lncRNAs were down-regulated. Furthermore, 386 matched lncRNA-mRNA pairs were detected for 120 differentially expressed lncRNAs and 362 differentially expressed genes in cis, and target gene-related pathway analyses exhibited significant variations in cysteine and methionine metabolism pathway-related genes. For the genes in trans, overall, 28,276 interaction relationships for 144 lncRNAs and differentially expressed protein-coding genes were detected. The pathway analyses found that secondary metabolites, such as phenylpropanoids and phenylalanine, and photosynthesis pathway-related genes were significantly altered by Cd stress. All of these results indicate that lncRNAs may regulate genes of cysteine-rich peptide metabolism in cis, as well as secondary metabolites and photosynthesis in trans, to activate various physiological and biochemical reactions to respond to excessive Cd. The present study could provide a valuable resource for lncRNA studies in response to Cd treatment in rice. It also expands our knowledge about lncRNA biological function and contributes to the annotation of the rice genome.

Despite the great contributions of utilizing heterosis to crop productivity worldwide, the molecular mechanism of heterosis remains largely unexplored. Thus, the present research is focused on the grain number heterosis of a widely used late-cropping indica super hybrid rice combination in China using a high-throughput next-generation RNA-seq strategy. Here, we obtained 872 million clean reads, and at least one read could maps 27,917 transcripts out of 35,679 annotations. Transcript differential expression analysis revealed a total of 5910 differentially expressed genes (DGHP) between super-hybrid rice Wufengyou T025 (WFYT025) and its parents were identified in the young panicles. Out of the 5910 DGHP, 63.1% had a genetic action mode of over-dominance, 17.3% had a complete-dominance action, 15.6% had a partial-dominance action and 4.0% had an additive action. DGHP were significantly enriched in carotenoid biosynthesis, diterpenoid biosynthesis and plant hormone signal transduction pathways, with the key genes involved in the three pathways being up-regulated in the hybrid. By comparing the DGHP enriched in the KEGG pathway with QTLs associated with grain number, several DGHP were located on the same chromosomal segment with some of these grain number QTLs. Through young panicle development transcriptome analysis, we conclude that the over-dominant effect is probably the major contributor to the grain number heterosis of WFYT025. The DGHP sharing the same location with grain number QTLs could be considered a candidate gene and provide valuable targets for the cloning and functional analysis of these grain number QTLs.

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). 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 °C); 2) HH: high nitrogen level (330 kg ha−1) with high temperature (37 °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. This RNA-seq study will help elucidate the molecular mechanisms of rice spikelet defense response to high temperature interacting with high nitrogen level.

Oryza longistaminata is an African wild rice species that possesses special traits for breeding applications. Self-incompatibility is the main cause of sterility in O. longistaminata, but here we demonstrated that its pollen vitality are normal. Lipid and carbohydrate metabolism were active throughout pollen development. In this study, we used I₂-KI staining and TTC staining to investigate pollen viability. Aniline-blue-stained semithin sections were used to investigate important stages of pollen development. Tandem mass tags (TMT)-based quantitative analysis was used to investigate the profiles of proteins related to lipid and carbohydrate metabolism in 4-, 6-, and 8.5-mm O. longistaminata spikelets before flowering. Pollen was found to germinate normally in vitro and in vivo. We documented cytological changes throughout important stages of anther development, including changes in reproductive cells as they formed mature pollen grains through meiosis and mitosis. A total of 31,987 RNA transcripts and 8,753 proteins were identified, and 6,842 of the proteins could be quantified. RNA-seq and proteome association analysis indicated that fatty acids were converted to sucrose after the 6-mm spikelet stage, based on the abundance of most key enzymes of the glyoxylate cycle and gluconeogenesis. The abundance of proteins involved in pollen energy metabolism was further confirmed by combining quantitative real-time PCR with parallel reaction monitoring (PRM) analyses. In conclusion, our study provides novel insights into the pollen viability of O. longistaminata at the proteome level, which can be used to improve the efficiency of male parent pollination in hybrid rice breeding applications.

The insect resistance of Bacillus thuringiensis (Bt) transgenic rice (Oryza sativa L.) is mainly determined by the transcription of the CRY1C gene and the translation process of the Cry1C protein. With this in mind, we analyzed CRY1C expression and Cry1C protein content in a transgenic Bt insect‐resistant restorer line and its F1 hybrids, evaluated the resistance to Chilo suppressalis of various lines of rice, and determined the 50% lethal concentration (LC₅₀) of Cry1C for C. suppressalis. In four transgenic rice lines, the relative expression of CRY1C was highest at the heading stage in most tissues. Among different tissues from the same developmental stage, CRY1C expression was highest in leaves, followed by stems and panicles. The relative expression of CRY1C was higher in the parent restorer line than in the F1 hybrids. The LC₅₀ of the Cry1C protein for C. suppressalis was 4.016 μg g⁻¹, and the Cry1C protein expression level of transgenic insect‐resistant rice exceeded this threshold at the heading stages and in stems. We next analyzed the lethality of each strain toward C. suppressalis. After 48 h, the mortality rate of second‐instar C. suppressalis larvae feeding on transgenic stem tissue was higher at the heading stage than at the tiller stage and varied from 66.7 to 91.7%. By documenting the temporal and spatial expression of CRY1C and evaluating C. suppressalis resistance to Bt transgenic rice, the risk of pests and diseases during the rice production process can be greatly reduced.