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.

Glutinous rice has very low amylose content and is a good material for determining the structure and physicochemical properties of amylopectin. We selected 29 glutinous rice varieties and determined the amylopectin structure by high-performance anion exchange chromatography with the pulsed amperometric detection method. We also determined the correlation between amylopectin structure and the physicochemical properties of starch extracted from these varieties. The results showed that the amylopectin chain ratio Σdegree of polymerization (DP) ≤ 11/ΣDP ≤ 24 of 29 glutinous rice varieties was greater than 0.26, signifying that these varieties contained type II amylopectin. The results of the correlation analysis with gelatinization temperature showed that ΣDP 6–11 was significantly negatively correlated with the onset gelatinization temperature (GT) (T_O), peak GT (T_P), and conclusion GT (T_C). Among the thermodynamic properties, ΣDP 12–24 was significantly positively correlated with To, Tp, and Tc, ΣDP 25–36 was significantly negatively correlated with To, Tp, and Tc, and ΣDP ≥ 37 had no correlation with the thermodynamic properties. The results of correlation analysis with RVA spectrum characteristic values showed that ΣDP 6–11 was significantly negatively correlated with hot paste viscosity (HPV), cool paste viscosity (CPV), consistency viscosity (CSV), peak time (PeT), and pasting temperature (PaT) among the Rapid Visco Analyzer (RVA) profile characteristics, ΣDP 12–24 was significantly positively correlated with HPV, CPV, CSV, PeT, and PaT, and ΣDP ≥ 25 had no correlation with the viscosity characteristics. Therefore, we concluded that the amylopectin structure had a greater effect on the T_O, T_P, T_C, ΔH and peak viscosity, HPV, CPV, CSV, PeT, and PaT. The glutinous rice varieties with a higher distribution of short chains and a lower distribution of medium and long chains in the amylopectin structure resulted in lower GT and RVA spectrum characteristic values.

Salt stress is one of the most important factors limiting rice growth and yield increase. Salt tolerance of rice at the bud burst (STB) stage determines whether germinated seeds can grow normally under salt stress, which is very important for direct seeding. However, reports on quantitative trait loci (QTLs) and candidate genes for STB in rice are very limited. In this study, a natural population of 130 indica and 81 japonica rice accessions was used to identify STB-related QTLs and candidate genes using a genome-wide association study (GWAS). Nine QTLs, including five for relative shoot length (RSL), two for relative root length (RRL), and two for relative root number (RRN), were identified. Five of these STB-related QTLs are located at the same site as the characterized salt tolerance genes, such as OsMDH1, OsSRFP1, and OsCDPK7. However, an important QTL related to RSL, qRSL1-2, has not been previously identified and was detected on chromosome 1. The candidate region for qRSL1-2 was identified by linkage disequilibrium analysis, 18 genes were found to have altered expression levels under salt stress through the RNA-seq database, and 10 of them were found to be highly expressed in the shoot. It was also found that, eight candidate genes (LOC_Os01g62980, LOC_Os01g63190, LOC_Os01g63230, LOC_Os01g63280, LOC_Os01g63400, LOC_Os01g63460, and LOC_Os01g63580) for qRSL1-2 carry different haplotypes between indica and japonica rice, which exactly corresponds to the significant difference in RSL values between indica and japonica rice in this study. Most of the accessions with elite haplotypes were indica rice, which had higher RSL values. These genes with indica-japonica specific haplotypes were identified as candidate genes. Rice accessions with elite haplotypes could be used as important resources for direct seeding. This study also provides new insights into the genetic mechanism of STB.

Highlight: Abrupt drought-flood alteration is a frequent meteorological disaster that occurs during summer in southern China and the Yangtze river basin, which often causes a large area reduction of rice yield. We previously reported abrupt drought-flood alteration effects on yield and its components, physiological characteristics, matter accumulation and translocation, rice quality of rice. However, the molecular mechanism of rice yield reduction caused by abrupt drought-flood alternation has not been reported. In this study, four treatments were provided, no drought and no floods (control), drought without floods (duration of drought 10 d), no drought with floods (duration of floods 8 d), and abrupt drought-flood alteration (duration of drought 10 d and floods 8 d). The quantitative analysis of spike metabolites was proceeded by LC-MS (liquid chromatograph-mass spectrometry) firstly. Then the Heat-map, PCA, PLS-DA, OPLS-DA and response ranking test of OPLS-DA model methods were used to analysis the function of differential metabolites (DMs) during the rice panicle differentiation stage under abrupt drought-flood alteration. In addition, relative quantitative analysis of spike total proteins under the treatment was conducted iTRAQ (isobaric tags for relative and absolute quantification) and LC-MS. In this study, 5708 proteins were identified and 4803 proteins were quantified. The identification and analysis of DEPs function suggested that abrupt drought-flood alteration treatment can promote carbohydrate metabolic, stress response, oxidation-reduction, defense response, and energy reserve metabolic process, etc, during panicle differentiation stage. In this study relative quantitative proteomics, metabolomics and physiology data (soluble protein content, superoxide dismutase activity, hydrogen peroxidase activity, peroxidase activity, malondialdehyde content, free proline content, soluble sugar content and net photosynthetic rate) analysis were applied to explicit the response mechanism of rice panicle differentiation stage under abrupt drought-flood alteration and provides a theoretical basis for the disaster prevention and mitigation. Abrupt drought-flood alternation is a meteorological disaster that frequently occurs during summer in southern China and the Yangtze river basin, often causing a significant loss of rice production. In this study, a quantitative analysis of spike metabolites was conducted via liquid chromatograph-mass spectrometry (LC-MS), and Heat-map, PCA, PLS-DA, OPLS-DA, and a response ranking test of OPLS-DA model methods were used to analyze functions of differential metabolites (DMs) during the rice panicle differentiation stage under abrupt drought-flood alternation. The results showed that 102 DMs were identified from the rice spike between T1 (abrupt drought-flood alternation) and CK0 (control) treatment, 104 DMs were identified between T1 and CK1 (drought) treatment and 116 DMs were identified between T1 and CK2 (flood) treatment. In addition, a relative quantitative analysis of spike total proteins was conducted using isobaric tags for relative and absolute quantification (iTRAQ) and LC-MS. The identification and analysis of DEPs functions indicates that abrupt drought-flood alternation treatment can promote carbohydrate metabolic, stress response, oxidation-reduction, defense response, and energy reserve metabolic process during the panicle differentiation stage. In this study, relative quantitative metabolomics and proteomics analyses were applied to explore the response mechanism of rice panicle differentiation in response to abrupt drought-flood alternation. Abbreviations: CK0 no drought and no floods CK1 drought without floods CK2 no drought with floods T1 abrupt drought-flood alteration LC-MS liquid chromatograph-mass spectrometry PCA principle component analysis (O)PLS-DA (orthogonal) partial least-squares-discriminant analysis DMs differential metabolites iTRAQ isobaric tags for relative and absolute quantification DEPs differentially expressed proteins KEGG kyoto encyclopedia of genes and genomes GO gene ontology SOD superoxide dismutase CAT hydrogen peroxidase POD peroxidase MDA malondialdehyde Pn net photosynthetic rate ROS reactive oxygen species VIP variable importance in the projection FC fold change