Integrated Analysis of Small RNA, Transcriptome, and Degradome Sequencing Reveals the MiR156, MiR5488 and MiR399 are Involved in the Regulation of Male Sterility in PTGMS Rice

Profiling of Known and Novel microRNAs in an Oleaginous Crop Native to the Amazon Basin, Sacha Inchi (Plukenetia volubilis), Through smallRNA-Seq

BACKGROUND

MicroRNAs (miRNAs) play crucial roles in regulating tissue-specific gene expression and plant development. This study explores the identification and functional characterization of miRNAs in Plukenetia volubilis (sacha inchi), an economically and nutritionally significant crop native to the Amazon basin, across three organs: root, stem, and leaf.

METHODS

Small RNA libraries were sequenced on the Illumina Novaseq 6000 platform, yielding high-quality reads that facilitated the discovery of known and novel miRNAs using miRDeep-P.

RESULTS

A total of 277 miRNAs were identified, comprising 71 conserved and 206 novel miRNAs, across root, stem, and leaf tissues. In addition, differential expression analysis using DESeq2 identified distinct miRNAs exhibiting tissue-specific regulation. Notably, novel miRNAs like novel_1, novel_88, and novel_189 showed significant roles in processes such as auxin signaling, lignin biosynthesis, and stress response. Functional enrichment analysis of miRNA target genes revealed pathways related to hormonal regulation, structural reinforcement, and environmental adaptation, highlighting tissue-specific functions. The Principal Component Analysis and PERMANOVA confirmed clear segregation of miRNA expression profiles among tissues, underlining organ-specific regulation. Differential expression patterns emphasized unique regulatory roles in each organ: roots prioritized stress response and nutrient uptake, leaves focused on photosynthesis and UV protection, and stems contributed to structural integrity and nutrient transport, suggesting evolutionary adaptations in P. volubilis.

CONCLUSIONS

This study identified novel miRNA-mediated networks that regulate developmental and adaptive processes in P. volubilis, underscoring its molecular adaptations for resilience and productivity. By characterizing both conserved and novel miRNAs, the findings lay a foundation for genetic improvement and molecular breeding strategies aimed at enhancing agronomic traits, stress tolerance, and the production of bioactive compounds.

Integrated analysis of transcriptome, sRNAome, and degradome involved in the drought-response of maize Zhengdan958

Drought is a major abiotic stress in restricting the growth, development, and yield of maize. As a significant epigenetic regulator, small RNA also functions in connecting the transcriptional and post-transcriptional regulatory network. Further to help comprehending the molecular mechanisms underlying drought adaptability and tolerance of maize, an integrated multi-omics analysis of transcriptome, sRNAome, and degradome was performed on the seedling roots of an elite hybrid Zhengdan958 under drought stress. In this study, 2,911 genes, 32 conserved miRNAs, and 12 novel miRNAs showed a significantly differential expression under drought stress. Moreover, 6,340 target genes of 445 miRNAs were validated using degradome sequencing, forming 281 miRNA-mRNA pairs in control (CK) and drought-stressed (DS) library. These target genes were mainly involved in the plant hormone signal transduction and phenylpropanoid biosynthesis pathways. The integrated multi-omics analysis revealed that five DEmiRNA-mRNA pairs displayed negatively correlated expression patterns, which were also verified by qRT-PCR. Tissue-specific expression profile and regulatory network analysis revealed that miR528a/b-Zm00001d021850, miR408a/b-Zm00001d020794, and miR164e-Zm00001d003414 might be essential in root-specific drought stress response of maize Zhengdan958 seedlings. These worthwhile will promote the functional characterization of miRNA-mRNA modules response to drought stress, and potentially contribute to drought-resistance breeding of maize.

Insights into the Early Steps of the Symbiotic Interaction between Soybean (Glycine max) and Sinorhizobium fredii Symbiosis Using Transcriptome, Small RNA, and Degradome Sequencing

Symbiotic nitrogen fixation carried out by the soybean-rhizobia symbiosis increases soybean yield and reduces the amount of nitrogen fertilizer that has been applied. MicroRNAs (miRNAs) are crucial in plant growth and development, prompting an investigation into their role in the symbiotic interaction of soybean with partner rhizobia. Through integrated small RNA, transcriptome, and degradome sequencing analysis, 1215 known miRNAs, 314 of them conserved, and 187 novel miRNAs were identified, with 44 differentially expressed miRNAs in soybean roots inoculated with Sinorhizobium fredii HH103 and a ttsI mutant. The study unveiled that the known miRNA gma-MIR398a-p5 was downregulated in the presence of the ttsI mutation, while the target gene of gma-MIR398a-p5, Glyma.06G007500, associated with nitrogen metabolism, was upregulated. The results of this study offer insights for breeding high-efficiency nitrogen-fixing soybean varieties, enhancing crop yield and quality.

Omics approaches to understand the MADS-box gene family in common bean (Phaseolus vulgaris L.) against drought stress

MADS-box genes are known to play important roles in diverse aspects of growth/devolopment and stress response in several plant species. However, no study has yet examined about MADS-box genes in P. vulgaris. In this study, a total of 79 PvMADS genes were identified and classified as type I and type II according to the phylogenetic analysis. While both type I and type II PvMADS classes were found to contain the MADS domain, the K domain was found to be present only in type II PvMADS proteins, in agreement with the literature. All chromosomes of the common bean were discovered to contain PvMADS genes and 17 paralogous gene pairs were identified. Only two of them were tandemly duplicated gene pairs (PvMADS-19/PvMADS-23 and PvMADS-20/PvMADS-24), and the remaining 15 paralogous gene pairs were segmentally duplicated genes. These duplications were found to play an important role in the expansion of type II PvMADS genes. Moreover, the RNAseq and RT-qPCR analyses showed the importance of PvMADS genes in response to drought stress in P. vulgaris.

The roles of non-coding RNAs in male reproductive development and abiotic stress responses during this unique process in flowering plants

Successful male reproductive development is the guarantee for sexual reproduction of flowering plants. Male reproductive development is a complicated and multi-stage process that integrates physiological processes and adaptation and tolerance to a myriad of environmental stresses. This well-coordinated process is governed by genetic and epigenetic machineries. Non-coding RNAs (ncRNAs) play pleiotropic roles in the plant growth and development. The identification, characterization and functional analysis of ncRNAs and their target genes have opened a new avenue for comprehensively revealing the regulatory network of male reproductive development and its response to environmental stresses in plants. This review briefly addresses the types, origin, biogenesis and mechanisms of ncRNAs in plants, highlights important updates on the roles of ncRNAs in regulating male reproductive development and emphasizes the contribution of ncRNAs, especially miRNAs and lncRNAs, in responses to abiotic stresses during this unique process in flowering plants.

Identification and Analysis of the MIR399 Gene Family in Grapevine Reveal Their Potential Functions in Abiotic Stress

MiR399 plays an important role in plant growth and development. The objective of the present study was to elucidate the evolutionary characteristics of the MIR399 gene family in grapevine and investigate its role in stress response. To comprehensively investigate the functions of miR399 in grapevine, nine members of the Vvi-MIR399 family were identified based on the genome, using a miRBase database search, located on four chromosomes (Chr 2, Chr 10, Chr 15, and Chr 16). The lengths of the Vvi-miR399 precursor sequences ranged from 82 to 122 nt and they formed stable stem-loop structures, indicating that they could produce microRNAs (miRNAs). Furthermore, our results suggested that the 2 to 20 nt region of miR399 mature sequences were relatively conserved among family members. Phylogenetic analysis revealed that the Vvi-MIR399 members of dicots (Arabidopsis, tomato, and sweet orange) and monocots (rice and grapevine) could be divided into three clades, and most of the Vvi-MIR399s were closely related to sweet orange in dicots. Promoter analysis of Vvi-MIR399s showed that the majority of the predicted cis-elements were related to stress response. A total of 66.7% (6/9) of the Vvi-MIR399 promoters harbored drought, GA, and SA response elements, and 44.4% (4/9) of the Vvi-MIRR399 promoters also presented elements involved in ABA and MeJA response. The expression trend of Vvi-MIR399s was consistent in different tissues, with the lowest expression level in mature and young fruits and the highest expression level in stems and young leaves. However, nine Vvi-MIR399s and four target genes showed different expression patterns when exposed to low light, high light, heat, cold, drought, and salt stress. Interestingly, a putative target of Vvi-MIR399 targeted multiple genes; for example, seven Vvi-MIR399s simultaneously targeted VIT_213s0067g03280.1. Furthermore, overexpression of Vvi_MIR399e and Vvi_MIR399f in Arabidopsis enhanced tolerance to drought compared with wild-type (WT). In contrast, the survival rate of Vvi_MIR399d-overexpressed plants were zero after drought stress. In conclusion, Vvi-MIR399e and Vvi-MIR399f, which are related to drought tolerance in grapevine, provide candidate genes for future drought resistance breeding.

Multi-omics analysis of small RNA, transcriptome, and degradome to identify putative miRNAs linked to MeJA regulated and oridonin biosynthesis in Isodon rubescens

Isodon rubescens has garnered much attention due to its anti-tumor or anti-cancer properties. However, little is known about the molecular mechanism of oridonin biosynthesis leveraging the regulatory network between small RNAs and mRNAs. In this study, the regulatory networks of miRNAs and targets were examined by combining mRNA, miRNA, and degradome. A total of 348 miRNAs, including 287 known miRNAs and 61 novel miRNAs, were identified. Among them, 51 miRNAs were significantly expressed, and 36 miRNAs responded to MeJA. A total of 3066 target genes were associated with 228 miRNAs via degradome sequencing. Multi-omics analysis demonstrated that 27 miRNA-mRNA pairs were speculated to be involved in MeJA regulation, and 36 miRNA-mRNA pairs were hypothesized to be involved in the genotype-dependence of I. rubescens. Furthermore, 151 and 7 miRNA-mRNA modules were likely engaged in oridonin biosynthesis as identified by psRNATarget and degradome sequencing, respectively. Some miRNA-mRNA modules were confirmed via RT-qPCR. Moreover, miRNAs targeting plant hormone signal transduction pathway genes were identified, such as miR156, miR167, miR393, and PC-3p-19822_242. Collectively, our results demonstrate for the first time that miRNAs are identified in I. rubescens, and laid a solid foundation for further research on the molecular mechanism of oridonin biosynthesis mediated by miRNA.

Heat-responsive microRNAs participate in regulating the pollen fertility stability of CMS-D2 restorer line under high-temperature stress

Anther development and pollen fertility of cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) restorer lines are susceptible to continuous high-temperature (HT) stress in summer, which seriously hinders the large-scale application of "three-line" hybrids in production. Here, integrated small RNA, transcriptome, degradome, and hormone profiling was performed to explore the roles of microRNAs (miRNAs) in regulating fertility stability in mature pollens of isonuclear alloplasmic near-isogenic restorer lines NH and SH under HT stress at two environments. A total of 211 known and 248 novel miRNAs were identified, of which 159 were differentially expressed miRNAs (DEMs). Additionally, 45 DEMs in 39 miRNA clusters (PmCs) were also identified, and most highly expressed miRNAs were significantly induced in SH under extreme HT, especially four MIR482 and six MIR6300 family miRNAs. PmC28 was located in the fine-mapped interval of the Rf1 gene and contained two DEMs, gra-miR482_L-2R + 2 and gma-miR2118a-3p_R + 1_1ss18TG. Transcriptome sequencing identified 6281 differentially expressed genes, of which heat shock protein (HSP)-related genes, such as HSP70, HSP22, HSP18.5-C, HSP18.2 and HSP17.3-B, presented significantly reduced expression levels in SH under HT stress. Through integrating multi-omics data, we constructed a comprehensive molecular network of miRNA-mRNA-gene-KEGG containing 35 pairs of miRNA/target genes involved in regulating the pollen development in response to HT, among which the mtr-miR167a_R + 1, tcc-miR167c and ghr-miR390a, tcc-miR396c_L-1 and ghr-MIR169b-p3_1ss6AG regulated the pollen fertility by influencing ARF8 responsible for the auxin signal transduction, ascorbate and aldarate metabolism, and the sugar and lipid metabolism and transport pathways, respectively. Further combination with hormone analysis revealed that HT-induced jasmonic acid signaling could activate the expression of downstream auxin synthesis-related genes and cause excessive auxin accumulation, followed by a cascade of auxin signal transduction, ultimately resulting in pollen abortion. The results provide a new understanding of how heat-responsive miRNAs regulate the stability of fertility restoration for CMS-D2 cotton under heat stress.

Cytological observation and RNA-seq analysis reveal novel miRNAs high expression associated with the pollen fertility of neo-tetraploid rice

BACKGROUND

Neo-tetraploid rice lines exhibit high fertility and strong heterosis and harbor novel specific alleles, which are useful germplasm for polyploid rice breeding. However, the mechanism of the fertility associated with miRNAs remains unknown. In this study, a neo-tetraploid rice line, termed Huaduo21 (H21), was used. Cytological observation and RNA-sequencing were employed to identify the fertility-related miRNAs in neo-tetraploid rice.

RESULTS

H21 showed high pollen fertility (88.08%), a lower percentage of the pollen mother cell (PMC) abnormalities, and lower abnormalities during double fertilization and embryogenesis compared with autotetraploid rice. A total of 166 non-additive miRNAs and 3108 non-additive genes were detected between H21 and its parents. GO and KEGG analysis of non-additive genes revealed significant enrichments in the DNA replication, Chromosome and associated proteins, and Replication and repair pathways. Comprehensive multi-omics analysis identified 32 pairs of miRNA/target that were associated with the fertility in H21. Of these, osa-miR408-3p and osa-miR528-5p displayed high expression patterns, targeted the phytocyanin genes, and were associated with high pollen fertility. Suppression of osa-miR528-5p in Huaduo1 resulted in a low seed set and a decrease in the number of grains. Moreover, transgenic analysis implied that osa-MIR397b-p3, osa-miR5492, and osa-MIR5495-p5 might participate in the fertility of H21.

CONCLUSION

Taken together, the regulation network of fertility-related miRNAs-targets pairs might contribute to the high seed setting in neo-tetraploid rice. These findings enhance our understanding of the regulatory mechanisms of pollen fertility associated with miRNAs in neo-tetraploid rice.

The role of microRNAs in responses to drought and heat stress in peanut (Arachis hypogaea)

MicroRNAs (miRNAs) are 21-24 nt small RNAs (sRNAs) that negatively regulate protein-coding genes and/or trigger phased small-interfering RNA (phasiRNA) production. Two thousand nine hundred miRNA families, of which ∼40 are deeply conserved, have been identified in ∼80 different plant species genomes. miRNA functions in response to abiotic stresses is less understood than their roles in development. Only seven peanut MIRNA families are documented in miRBase, yet a reference genome assembly is now published and over 480 plant-like MIRNA loci were predicted in the diploid peanut progenitor Arachis duranensis genome. We explored by computational analysis of a leaf sRNA library and publicly available sRNA, degradome, and transcriptome datasets the miRNA and phasiRNA space associated with drought and heat stresses in peanut. We characterized 33 novel candidate and 33 ancient conserved families of MIRNAs and present degradome evidence for their cleavage activities on mRNA targets, including several noncanonical targets and novel phasiRNA-producing noncoding and mRNA loci with validated novel targets such as miR1509 targeting serine/threonine-protein phosphatase7 and miRc20 and ahy-miR3514 targeting penta-tricopeptide repeats (PPRs), in contradistinction to other claims of miR1509/173/7122 superfamily miRNAs indirectly targeting PPRs via TAS-like noncoding RNA loci. We characterized the inverse correlations of significantly differentially expressed drought- and heat-regulated miRNAs, assayed by sRNA blots or transcriptome datasets, with target mRNA expressions in the same datasets. Meta-analysis of an expression atlas and over representation of miRNA target genes in co-expression networks suggest that miRNAs have functions in unique aspects of peanut gynophore development. Genome-wide MIRNA annotation of the published allopolyploid peanut genome can facilitate molecular breeding of value-added traits.

Plant non-coding RNAs function in pollen development and male sterility

Male sterility is classified as either cytoplasmic male sterility (CMS) or genic male sterility (GMS). Generally, CMS involves mitochondrial genomes interacting with the nuclear genome, while GMS is caused by nuclear genes alone. Male sterility is regulated by multilevel mechanisms in which non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and phased small interfering RNAs (phasiRNAs), which have been proven to be critical elements. The development of high-throughput sequencing technology offers new opportunities to evaluate the genetic mechanism of ncRNAs in plant male sterility. In this review, we summarize the critical ncRNAs that regulate gene expression in ways dependent on or independent of hormones, which involve the differentiation of the stamen primordia, degradation of the tapetum, formation of microspores, and the release of pollen. In addition, the key mechanisms of the miRNA-lncRNA-mRNA interaction networks mediating male sterility in plants are elaborated. We present a different perspective on exploring the ncRNA-mediated regulatory pathways that control CMS in plants and create male-sterile lines through hormones or genome editing. A refined understanding of the ncRNA regulatory mechanisms in plant male sterility for the development of new sterile lines would be conducive to improve hybridization breeding.

Novel insights on genes and pathways involved in Pinus elliottii response to resinosis

Pinus elliottii, an important coniferous timber species, has recently become one of the most popular sources of resin in China. Resinosis is a common disease that may negatively affect pine tree growth and production. In this study, we used single-molecule real-time sequencing and Illumina RNA sequencing to generate an accurate transcriptome for P. elliottii. The transcriptome included 90,026 transcripts, 5160 long non-coding RNAs and 7710 transcription factors. We then analyzed RNA-sequencing, small RNA-sequencing and degradome data to identify genes, miRNAs and key miRNA-target pairs involved in response to resinosis in P. elliottii. We identified 1305 genes and 1151 miRNAs exhibiting significant differential expression in response to resinosis. According to the degradome sequencing analysis, 318 differentially expressed transcripts were targets of 14 differentially expressed miRNAs. Our study has provided resources for further functional characterization of genes and miRNAs involved in resinosis in P. elliottii, which should aid the future disease-resistance breeding of this species.

Integrated analysis of miRNAome transcriptome and degradome reveals miRNA-target modules governing floral florescence development and senescence across early- and late-flowering genotypes in tree peony

As a candidate national flower of China, tree peony has extremely high ornamental, medicinal and oil value. However, the short florescence and rarity of early-flowering and late-flowering varieties restrict further improvement of the economic value of tree peony. Specific miRNAs and their target genes engaged in tree peony floral florescence, development and senescence remain unknown. This report presents the integrated analysis of the miRNAome, transcriptome and degradome of tree peony petals collected from blooming, initial flowering, full blooming and decay stages in early-flowering variety Paeonia ostii 'Fengdan', an early-flowering mutant line of Paeonia ostii 'Fengdan' and late-flowering variety Paeonia suffruticosa 'Lianhe'. Transcriptome analysis revealed a transcript ('psu.G.00014095') which was annotated as a xyloglucan endotransglycosylase/hydrolase precursor XTH-25 and found to be differentially expressed across flower developmental stages in Paeonia ostii 'Fengdan' and Paeonia suffruticosa 'Lianhe'. The miRNA-mRNA modules were presented significant enrichment in various pathways such as plant hormone signal transduction, indole alkaloid biosynthesis, arachidonic acid metabolism, folate biosynthesis, fatty acid elongation, and the MAPK signaling pathway. Multiple miRNA-mRNA-TF modules demonstrated the potential functions of MYB-related, bHLH, Trihelix, NAC, GRAS and HD-ZIP TF families in floral florescence, development, and senescence of tree peony. Comparative spatio-temporal expression investigation of eight floral-favored miRNA-target modules suggested that transcript 'psu.T.00024044' and microRNA mtr-miR166g-5p are involved in the floral florescence, development and senescence associated agronomic traits of tree peony. The results might accelerate the understanding of the potential regulation mechanism in regards to floral florescence, development and abscission, and supply guidance for tree peony breeding of varieties with later and longer florescence characteristics.

Roles of microRNAs in abiotic stress response and characteristics regulation of plant

MicroRNAs (miRNAs) are a class of non-coding endogenous small RNAs (long 20-24 nucleotides) that negatively regulate eukaryotes gene expression at post-transcriptional level via cleavage or/and translational inhibition of targeting mRNA. Based on the diverse roles of miRNA in regulating eukaryotes gene expression, research on the identification of miRNA target genes has been carried out, and a growing body of research has demonstrated that miRNAs act on target genes and are involved in various biological functions of plants. It has an important influence on plant growth and development, morphogenesis, and stress response. Recent case studies indicate that miRNA-mediated regulation pattern may improve agronomic properties and confer abiotic stress resistance of plants, so as to ensure sustainable agricultural production. In this regard, we focus on the recent updates on miRNAs and their targets involved in responding to abiotic stress including low temperature, high temperature, drought, soil salinity, and heavy metals, as well as plant-growing development. In particular, this review highlights the diverse functions of miRNAs on achieving the desirable agronomic traits in important crops. Herein, the main research strategies of miRNAs involved in abiotic stress resistance and crop traits improvement were summarized. Furthermore, the miRNA-related challenges and future perspectives of plants have been discussed. miRNA-based research lays the foundation for exploring miRNA regulatory mechanism, which aims to provide insights into a potential form of crop improvement and stress resistance breeding.

OsSPLs Regulate Male Fertility in Response to Different Temperatures by Flavonoid Biosynthesis and Tapetum PCD in PTGMS Rice

Photoperiod and thermo-sensitive genic male sterile (PTGMS) rice is an important resource for two line hybrid rice production. The SQUAMOSA–promoter binding, such as the (SPL) gene family, encode the plant specific transcription factors that regulate development and defense responses in plants. However, the reports about SPLs participating in male fertility regulation are limited. Here, we identified 19 OsSPL family members and investigated their involvement in the fertility regulation of the PTGMS rice lines, PA2364S and PA2864S, with different fertility transition temperatures. The results demonstrated that OsSPL2, OsSPL4, OsSPL16 and OsSPL17 affect male fertility in response to temperature changes through the MiR156-SPL module. WGCNA (weighted gene co-expression network analysis) revealed that CHI and APX1 were co-expressed with OsSPL17. Targeted metabolite and flavonoid biosynthetic gene expression analysis revealed that OsSPL17 regulates the expression of flavonoid biosynthesis genes CHI, and the up regulation of flavanones (eriodictvol and naringenin) and flavones (apigenin and luteolin) content contributed to plant fertility. Meanwhile, OsSPL17 negatively regulates APX1 to affect APX (ascorbate peroxidase) activity, thereby regulating ROS (reactive oxygen species) content in the tapetum, controlling the PCD (programmed cell death) process and regulating male fertility in rice. Overall, this report highlights the potential role of OsSPL for the regulation of male fertility in rice and provides a new insight for the further understanding of fertility molecular mechanisms in PTGMS rice.

An Integrated Analysis of Transcriptome and miRNA Sequencing Provides Insights into the Dynamic Regulations during Flower Morphogenesis in Petunia

Published genome sequences can facilitate multiple genome sequencing studies of flower development, which can serve as the basis for later analysis of variation in flower phenotypes. To identify potential regulators related to flower morphology, we captured dynamic expression patterns under five different developmental stages of petunia flowers, a popular bedding plant, using transcriptome and miRNA sequencing. The significant transcription factor (TF) families, including MYB, MADS, and bHLH, were elucidated. MADS-box genes exhibited co-expression patterns with BBR-BPC, GATA, and Dof genes in different modules according to a weighted gene co-expression network analysis. Through miRNA sequencing, a total of 45 conserved and 26 novel miRNAs were identified. According to GO and KEGG enrichment analysis, the carbohydrate metabolic process, photosynthesis, and phenylalanine metabolism were significant at the transcriptomic level, while the response to hormone pathways was significantly enriched by DEmiR-targeted genes. Finally, an miRNA–RNA network was constructed, which suggested the possibility of novel miRNA-mediated regulation pathways being activated during flower development. Overall, the expression data in the present study provide novel insights into the developmental gene regulatory network facilitated by TFs, miRNA, and their target genes.

Combined analysis of transcriptome and metabolome reveals that sugar, lipid, and phenylpropane metabolism are essential for male fertility in temperature-induced male sterile rice

Photoperiod- and thermosensitive genic male sterility (PTGMS) rice is a vital germplasm resource consisting of two-line hybrid rice in which light and temperature strictly control their fertility changes. Variable environmental conditions present huge risks to the two-lines hybrid seed production. Explaining the regulatory mechanism of male fertility in rice PTGMS lines is an essential prerequisite to ensuring food security production. A group of near-isogenic lines (NILs) of a rice PTGMS line unique to this research group was used for this study. These lines have the same genetic background and regulate male fertility by responding to different temperature changes. Transcriptomic analysis revealed that 315 upregulated genes and 391 regulated genes regulated male fertility in response to temperature changes, and differentially expressed genes (DEGs) were mainly characterized in enrichment analysis as having roles in the metabolic pathways of sugar, lipid and phenylpropanoid. Electron microscopy analysis revealed that a lack of starch accumulation in sterile pollen grains induced by high temperature, with an abnormal exine development and a lack of inner pollen grains. Defective processes for sporopollenin synthesis, sporopollenin transport and pollen wall formation in sterile anthers were verified using qPCR. Targeted metabolomics analysis revealed that most lipids (phospholipids, sphingolipids and fatty acids) and flavonoids (flavones and flavanones) were upregulated in fertile anthers and involved in pollen wall development and male fertility formation, while lignin G units and C-type lignin were the major contributors to pollen wall development. The coding genes for trehalose 6-phosphate phosphatase, beta-1,3-glucanase, phospholipase D and 4-coumarate-CoA ligase are considered essential regulators in the process of male fertility formation. In conclusion, our results indicated that the expression of critical genes and accumulation of metabolites in the metabolism of sugar, lipid, and phenylpropanoid are essential for male fertility formation. The results provide new insights for addressing the negative effects of environmental variation on two-line hybrid rice production.

Jasmonic Acid Plays a Pivotal Role in Pollen Development and Fertility Regulation in Different Types of P(T)GMS Rice Lines

Two-line hybrid rice systems represent a new technical approach to utilizing the advantages of rice hybrids. However, the mechanism underlying the male sterile-line fertility transition in rice remains unclear. Peiai 64S (PA64S) is a photoperiod- and thermo-sensitive genic male sterile (PTGMS) line in which male sterility manifests at an average temperature above 23.5 °C under long-day (LD) conditions. Nongken 58S (NK58S) is a LD-sensitive genic male sterile (PGMS) rice that is sterile under LD conditions (above 13.75 h-day). In contrast, D52S is a short-day (SD)-PGMS line that manifests male sterility under SD conditions (below 13.5 h-day). In this study, we obtained fertile and sterile plants from all three lines and performed transcriptome analyses on the anthers of the plants. Gene ontology (GO) analysis suggested that the differentially expressed genes identified were significantly enriched in common terms involved in the response to jasmonic acid (JA) and in JA biosynthesis. On the basis of the biochemical and molecular validation of dynamic, tissue-specific changes in JA, indole-3-acetic acid (IAA) levels, gibberellin (GA) levels, and JA biosynthetic enzyme activities and expression, we proposed that JA could play a pivotal role in viable pollen production through its initial upregulation, constant fluctuation and leaf-spikelet signaling under certain fertility-inducing conditions. Furthermore, we also sprayed methyl jasmonate (MEJA) and salicylhydroxamic acid (SHAM) on the plants, thereby achieving fertility reversal in the PGMS lines NK58S and D52S, with 12.91–63.53% pollen fertility changes. Through qPCR and enzyme activity analyses, we identified two key enzymes—allene oxide synthase (AOS) and allene oxide cyclase (AOC)—that were produced and upregulated by 20–500-fold in PGMS in response to spraying; the activities of these enzymes reversed pollen fertility by influencing the JA biosynthetic pathway. These results provide a new understanding of hormone interactions and networks in male-sterile rice based on the role of JA that will help us to better understand the potential regulatory mechanisms of fertility development in rice in the future.