Small RNA and Transcriptome Sequencing Reveals miRNA Regulation of Floral Thermogenesis in Nelumbo nucifera

Gene expression and metabolite levels converge in the thermogenic spadix of skunk cabbage

The inflorescence (spadix) of skunk cabbage (Symplocarpus renifolius) is strongly thermogenic and can regulate its temperature at around 23 °C even when the ambient temperature drops below freezing. To elucidate the mechanisms underlying developmentally controlled thermogenesis and thermoregulation in skunk cabbage, we conducted a comprehensive transcriptome and metabolome analysis across 3 developmental stages of spadix development. Our RNA-seq analysis revealed distinct groups of expressed genes, with selenium-binding protein 1/methanethiol oxidase (SBP1/MTO) exhibiting the highest levels in thermogenic florets. Notably, the expression of alternative oxidase (AOX) was consistently high from the prethermogenic stage through the thermogenic stage in the florets. Metabolome analysis showed that alterations in nucleotide levels correspond with the developmentally controlled and tissue-specific thermogenesis of skunk cabbage, evident by a substantial increase in AMP levels in thermogenic florets. Our study also reveals that hydrogen sulfide, a product of SBP1/MTO, inhibits cytochrome c oxidase (COX)-mediated mitochondrial respiration, while AOX-mediated respiration remains relatively unaffected. Specifically, at lower temperatures, the inhibitory effect of hydrogen sulfide on COX-mediated respiration increases, promoting a shift toward the dominance of AOX-mediated respiration. Finally, despite the differential regulation of genes and metabolites throughout spadix development, we observed a convergence of gene expression and metabolite accumulation patterns during thermogenesis. This synchrony may play a key role in developmentally regulated thermogenesis. Moreover, such convergence during the thermogenic stage in the spadix may provide a solid molecular basis for thermoregulation in skunk cabbage.

The potential mechanism of response to light intensity in energy metabolism mediated by miRNA in Isatis indigotica

Light is the main source of energy for plant growth. Studies have shown that I. indigotica is a light-demanding plant and its yield and various active components are positively correlated with light intensity, but no studies of light intensity affecting energy metabolism in I. indigotica have been reported. Mitochondria are the main site of energy metabolism, and miRNAs are important factors in regulating gene expression, this experiment attempts to study the effects of different light intensities on energy metabolism from the perspective of mitochondria and miRNAs. The results show that the biomass、mitochondrial structural integrity and energy metabolism in I. indigotica were found to be positively correlated with light intensity. Small RNA and transcriptome sequencing identified 241 miRNAs and 36,372 mRNAs, and degradomic technology identified 72 miRNAs targeting 106 mRNAs, among which 12 pairs of miRNA-mRNAs were annotated on mitochondria. Combined with RT-qPCR validation, it was concluded that miR167a-5p positively regulates LETM1 and affects mitochondrial structure, miR400-5p and mIR169m-p3_1ss15CT negatively regulate GRXS15 and CMC4, respectively, affecting SDH and CCO activities, and miR395a-APS4 may affect the utilization of ATP and sulfate assimilation. In summary, the results of this study complement and enrich knowledge of light effects on mitochondria from the perspective of miRNA, while providing guidance for the cultivation of I. indigotica.

Identification of long non-coding RNA-microRNA-mRNA regulatory modules and their potential roles in drought stress response in wheat (Triticum aestivum L.)

Drought is one of the most severe abiotic stresses that influence wheat production across the globe. Understanding the molecular regulatory network of wheat in response to drought is of great importance in molecular breeding. Noncoding RNAs influence plant development and resistance to abiotic stresses by regulating gene expression. In this study, whole-transcriptome sequencing was performed on the seedlings of two wheat varieties with contrasting levels of drought tolerance under drought and control conditions to identify long noncoding RNAs (lncRNAs), micro RNAs (miRNAs), and mRNAs related to drought stress and explore the potential lncRNA-miRNA-mRNA regulatory modules in controlling wheat drought stress response. A total of 1515 differentially expressed lncRNAs (DELs), 209 differentially expressed miRNAs (DEMs), and 20462 differentially expressed genes (DEGs) were identified. Of the 20462 DEGs, 1025 were identified as potential wheat drought resistance-related DEGs. Based on the regulatory relationship and expression patterns of DELs, DEMs, and DEGs, 10 DEL-DEM-DEG regulatory modules related to wheat drought stress response were screened, and preliminary expression verification of two important candidate modules was performed. Our results revealed the possible roles of lncRNA-miRNA-mRNA modules in regulatory networks related to drought tolerance and provided useful information as valuable genomic resources in molecular breeding of wheat.

The Genome-Wide Identification of Long Non-Coding RNAs Involved in Floral Thermogenesis in Nelumbo nucifera Gaertn

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.

Integrated Analysis of miRNAs Associated With Sugarcane Responses to Low-Potassium Stress

Sugarcane is among the most important global crops and a key bioenergy source. Sugarcane production is restricted by limited levels of available soil potassium (K⁺). The ability of plants to respond to stressors can be regulated by a range of microRNAs (miRNAs). However, there have been few studies regarding the roles of miRNAs in the regulation of sugarcane responses to K⁺-deficiency. To understand how these non-coding RNAs may influence sugarcane responses to low-K⁺ stress, we conducted expression profiling of miRNAs in sugarcane roots under low-K⁺ conditions via high-throughput sequencing. This approach led to the identification of 324 and 42 known and novel miRNAs, respectively, of which 36 were found to be differentially expressed miRNAs (DEMs) under low-K⁺ conditions. These results also suggested that miR156-x/z and miR171-x are involved in these responses as potential regulators of lateral root formation and the ethylene signaling pathway, respectively. A total of 705 putative targets of these DEMs were further identified through bioinformatics predictions and degradome analyses, and GO and KEGG enrichment analyses revealed these target mRNAs to be enriched for catalytic activity, binding functions, metabolic processes, plant hormone signal transduction, and mitogen-activated protein kinase (MAPK) signaling. In summary, these data provide an overview of the roles of miRNAs in the regulation of sugarcane response to low-K⁺ conditions.

Characterization of Genomic Variation from Lotus (Nelumbo Adans.) Mutants with Wide and Narrow Tepals

Compared with rose, chrysanthemum, and water lily, the absence of short-wide and long-narrow tepals of ornamental lotus (Nelumbo Adans.) limits the commercial value of flowers. In this study, the genomes of two groups of lotus mutants with wide-short and narrow-long tepals were resequenced to uncover the genomic variation and candidate genes associated with tepal shape. In group NL (short for N. lutea, containing two mutants and one control of N. lutea), 716,656 single nucleotide polymorphisms (SNPs) and 221,688 insertion-deletion mutations (Indels) were obtained, while 639,953 SNPs and 134,6118 Indels were obtained in group WSH (short for ‘Weishan Hong’, containing one mutant and two controls of N. nucifera ‘Weishan Hong’). Only a small proportion of these SNPs and Indels was mapped to exonic regions of genome: 1.92% and 0.47%, respectively, in the NL group, and 1.66% and 0.48%, respectively, in the WSH group. Gene Ontology (GO) analysis showed that out of 4890 (NL group) and 1272 (WSH group) annotated variant genes, 125 and 62 genes were enriched (Q < 0.05), respectively. Additionally, in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, 104 genes (NL group) and 35 genes (WSH group) were selected (p < 0.05). Finally, there were 306 candidate genes that were sieved to determine the development of tepal shape in lotus plants. It will be an essential reference for future identification of tepal-shaped control genes in lotus plants. This is the first comprehensive report of genomic variation controlling tepal shape in lotus, and the mutants in this study are promising materials for breeding novel lotus cultivars with special tepals.

DIA-Based Quantitative Proteomics Reveals the Protein Regulatory Networks of Floral Thermogenesis in Nelumbo nucifera

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.

Integration of Transcriptional and Post-transcriptional Analysis Revealed the Early Response Mechanism of Sugarcane to Cold Stress

Cold stress causes major losses to sugarcane production, yet the precise molecular mechanisms that cause losses due to cold stress are not well-understood. To survey miRNAs and genes involved in cold tolerance, RNA-seq, miRNA-seq, and integration analyses were performed on Saccharum spontaneum. Results showed that a total of 118,015 genes and 6,034 of these differentially expressed genes (DEGs) were screened. Protein–protein interaction (PPI) analyses revealed that ABA signaling via protein phosphatase 2Cs was the most important signal transduction pathway and late embryogenesis abundant protein was the hub protein associated with adaptation to cold stress. Furthermore, a total of 856 miRNAs were identified in this study and 109 of them were differentially expressed in sugarcane responding to cold stress. Most importantly, the miRNA–gene regulatory networks suggested the complex post-transcriptional regulation in sugarcane under cold stress, including 10 miRNAs−42 genes, 16 miRNAs−70 genes, and three miRNAs−18 genes in CT vs. LT0.5, CT vs. LT1, and CT0.5 vs. LT1, respectively. Specifically, key regulators from 16 genes encoding laccase were targeted by novel-Chr4C_47059 and Novel-Chr4A_40498, while five LRR-RLK genes were targeted by Novel-Chr6B_65233 and Novel-Chr5D_60023, 19 PPR repeat proteins by Novel-Chr5C_57213 and Novel-Chr5D_58065. Our findings suggested that these miRNAs and cell wall-related genes played vital regulatory roles in the responses of sugarcane to cold stress. Overall, the results of this study provide insights into the transcriptional and post-transcriptional regulatory network underlying the responses of sugarcane to cold stress.