Pentatricopeptide-repeat proteins (PPRs) are characterized by tandem arrays of a degenerate 35-amino-acid (PPR motifs), which can bind RNA strands and participate in post-transcription. PPR proteins family is one of the largest families in land plants and play important roles in organelle RNA metabolism and plant development. However, the functions of PPR genes involved in biotic and abiotic stresses of rice (Oryza sativa L.) remain largely unknown. In the present study, a comprehensive genome-wide analysis of PPR genes was performed. A total of 491 PPR genes were found in the rice genome, of which 246 PPR genes belong to the P subfamily, and 245 genes belong to the PLS subfamily. Gene structure analysis showed that most PPR genes lack intron. Chromosomal location analysis indicated that PPR genes were widely distributed in all 12 rice chromosomes. Phylogenetic relationship analysis revealed the distinct difference between the P and PLS subfamilies. Many PPR proteins are predicted to target chloroplasts or mitochondria, and a PPR protein (LOC_Os10g34310) was verified to localize in mitochondria. Furthermore, three PPR genes (LOC_Os03g17634,LOC_Os07g40820,LOC_Os04g51350) were verified as corresponding miRNA targets. The expression pattern analysis showed that many PPR genes could be induced under biotic and abiotic stresses. Finally, seven PPR genes were confirmed with their expression patterns under salinity or drought stress. We found 491 PPR genes in the rice genome, and our genes structure analysis and syntenic analysis indicated that PPR genes might be derived from amplification by retro-transposition. The expression pattern present here suggested that PPR proteins have crucial roles in response to different abiotic stresses in rice. Taken together, our study provides a comprehensive analysis of the PPR gene family and will facilitate further studies on their roles in rice growth and development.

Circular RNAs (circRNAs) are known to play an important role in the regulation of gene expression in eukaryotes. Photo-thermosensitive genic male sterile (PTGMS) is a very important germplasm resource in two-line hybrid rice breeding. Although many circRNAs have been identified in rice (Oryza sativa L.), little is known about the biological roles of circRNAs in the fertility transition of the PTGMS rice line. In the present study, RNA-sequencing libraries were constructed from the young panicles of the Wuxiang S sterile line rice (WXS (S)) and its fertile line rice (WXS (F)) at three development stages with three biological replicates. A total of 9994 circRNAs were obtained in WXS rice based on high-throughput strand-specific RNA sequencing and bioinformatic approaches, of which 5305 were known circRNAs and 4689 were novel in rice. And 14 of 16 randomly selected circRNAs were experimentally validated with divergent primers. Our results showed that 186 circRNAs were significantly differentially expressed in WXS (F) compared with WXS (S), of which 97, 87 and 60 circRNAs were differentially expressed at the pollen mother cell (PMC) formation stage (P2), the meiosis stage (P3) and the microspore formation stage (P4), respectively. Fertility specific expression patterns of eight circRNAs were analysis by qRT-PCR. Gene ontology (GO) and KEGG pathway analysis of the parental genes of differentially expressed circRNAs (DECs) revealed that they mainly participated in various biological processes such as development, response to stimulation, hormonal regulation, and reproduction. Furthermore, 15 DECs were found to act as putative miRNA sponges to involved in fertility transition in PTGMS rice line. In the present study, the abundance and characteristics of circRNAs were investigated in the PTGMS rice line using bioinformatic approaches. Moreover, the expression patterns of circRNAs were different between WXS (F) and WXS (S). Our findings primarily revealed that circRNAs might be endogenous noncoding regulators of flower and pollen development, and were involved in the fertility transition in the PTGMS rice line, and guide the production and application of two-line hybrid rice.

The sacred lotus (Nelumbo nucifera) can maintain a stable floral chamber temperature between 30 and 35 °C when blooming despite fluctuations in ambient temperatures between about 8 and 45 °C, but the regulatory mechanism of floral thermogenesis remains unclear. Here, we obtained comprehensive protein profiles from receptacle tissue at five developmental stages using data-independent acquisition (DIA)-based quantitative proteomics technology to reveal the molecular basis of floral thermogenesis of N. nucifera. A total of 6913 proteins were identified and quantified, of which 3513 differentially abundant proteins (DAPs) were screened. Among them, 640 highly abundant proteins during the thermogenic stages were mainly involved in carbon metabolism processes such as the tricarboxylic acid (TCA) cycle. Citrate synthase was identified as the most connected protein in the protein-protein interaction (PPI) network. Next, the content of alternative oxidase (AOX) and plant uncoupling protein (pUCP) in different tissues indicated that AOX was specifically abundant in the receptacles. Subsequently, a protein module highly related to the thermogenic phenotype was identified by the weighted gene co-expression network analysis (WGCNA). In summary, the regulation mechanism of floral thermogenesis in N. nucifera involves complex regulatory networks, including TCA cycle metabolism, starch and sucrose metabolism, fatty acid degradation, and ubiquinone synthesis, etc.

The sacred lotus (Nelumbo nucifera Gaertn.) can produce heat autonomously and maintain a relatively stable floral chamber temperature for several days when blooming. Floral thermogenesis is critical for flower organ development and reproductive success. However, the regulatory role of microRNA (miRNA) underlying floral thermogenesis in N. nucifera remains unclear. To comprehensively understand the miRNA regulatory mechanism of thermogenesis, we performed small RNA sequencing and transcriptome sequencing on receptacles from five different developmental stages. In the present study, a total of 172 known miRNAs belonging to 39 miRNA families and 126 novel miRNAs were identified. Twenty-nine thermogenesis-related miRNAs and 3024 thermogenesis-related mRNAs were screened based on their expression patterns. Of those, seventeen differentially expressed miRNAs (DEMs) and 1765 differentially expressed genes (DEGs) had higher expression during thermogenic stages. The upregulated genes in the thermogenic stages were mainly associated with mitochondrial function, oxidoreductase activity, and the energy metabolism process. Further analysis showed that miR156_2, miR395a_5, miR481d, and miR319p may play an important role in heat-producing activity by regulating cellular respiration-related genes. This study provides comprehensive miRNA and mRNA expression profile of receptacle during thermogenesis in N. nucifera, which advances our understanding on the regulation of floral thermogenesis mediated by miRNA.

The two-line hybrid rice is an important factor of a global crop, but its fertility transition mechanism is unclear. Here, a comparative proteomics and transcriptomics analysis was completed on the two-line hybrid rice line Wuxiang S (WXS) to explore its molecular mechanism and protein regulation during fertility transition. A total of 340 differentially abundant proteins (DAPs) were identified using iTRAQ between the pollen mother cell formation stage (P2) and the meiosis stage (P3). There were 3541 and 4247 differentially expressed genes (DEGs) in P2 and P3 between WXS (Sterile, S)-WXS(S) and WXS (Fertile, F)-WXS(F), respectively, of which 92 and 71 DEGs had corresponding DAPs. Among the DAPs and DEGs, 65 (SP2 vs. FP2) and 55 (SP3 vs. FP3) corresponding DEGs and DAPs (cor-DEGs-DAPs) showed the same expression trend, indicating the cor-DEGs-DAPs genes might play vital roles in WXS fertility transition. Further analysis indicated that cor-DEGs-DAPs proteins were related to energy metabolism-related proteins in anther development and were accompanied by the activation of the stress response pathway and modifications to the cell wall, which ultimately affected the fertility transition of the PTGMS rice line WXS.

The sacred lotus (Nelumbo nucifera Gaertn.) can maintain a stable floral chamber temperature when blooming, despite ambient temperature fluctuations; however, the long non-coding RNAs (lncRNAs) involved in floral thermogenesis remain unclear. In the present study, we obtain comprehensive lncRNAs expression profiles from receptacles at five developmental stages by strand-specific RNA sequencing to reveal the lncRNAs regulatory mechanism of the floral thermogenesis of N. nucifera. A total of 22,693 transcripts were identified as lncRNAs, of which approximately 44.78% had stage-specific expression patterns. Subsequently, we identified 2579 differential expressed lncRNAs (DELs) regulating 2367 protein-coding genes mainly involved in receptacle development and reproductive process. Then, lncRNAs with floral thermogenesis identified by weighted gene co-expression network analysis (WGCNA) were mainly related to sulfur metabolism and mitochondrial electron transport chains. Meanwhile, 70 lncRNAs were predicted to act as endogenous target mimics (eTMs) for 29 miRNAs and participate in the regulation of 16 floral thermogenesis-related genes. Our dual luciferase reporter assays indicated that lncRNA LTCONS_00068702 acted as eTMs for miR164a_4 to regulate the expression of TrxL2 gene. These results deepen our understanding of the regulation mechanism of floral thermogenesis by lncRNAs and accumulate data for further research.