The miR399-CsUBC24 Module Regulates Reproductive Development and Male Fertility in Citrus

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.

Genome-wide identification of the C2H2-Zinc finger gene family and functional validation of CsZFP7 in citrus nucellar embryogenesis

KEY MESSAGE

Genome-wide identification of C2H2-ZF gene family in the poly- and mono-embryonic citrus species and validation of the positive role of CsZFP7 in sporophytic apomixis. The C2H2 zinc finger (C2H2-ZF) gene family is involved in plant vegetative and reproductive development. Although a large number of C2H2 zinc-finger proteins (C2H2-ZFPs) have been well characterized in some horticultural plants, little is known about the C2H2-ZFPs and their function in citrus. In this work, we performed a genome-wide sequence analysis and identified 97 and 101 putative C2H2-ZF gene family members in the genomes of sweet orange (C. sinensis, poly-embryonic) and pummelo (C. grandis, mono-embryonic), respectively. Phylogenetic analysis categorized citrus C2H2-ZF gene family into four clades, and their possible functions were inferred. According to the numerous regulatory elements on promoter, citrus C2H2-ZFPs can be divided into five different regulatory function types that indicate functional differentiation. RNA-seq data revealed 20 differentially expressed C2H2-ZF genes between poly-embryonic and mono-embryonic ovules at two stages of citrus nucellar embryogenesis, among them CsZFP52 specifically expressed in mono-embryonic pummelo ovules, while CsZFP7, 37, 44, 45, 67 and 68 specifically expressed in poly-embryonic sweet orange ovules. RT-qPCR further validated that CsZFP7 specifically expressed at higher levels in poly-embryonic ovules, and down-regulation of CsZFP7 in the poly-embryonic mini citrus (Fortunella hindsii) increased rate of mono-embryonic seeds compared with the wild type, indicating the regulatory potential of CsZFP7 in nucellar embryogenesis of citrus. This work provided a comprehensive analysis of C2H2-ZF gene family in citrus, including genome organization and gene structure, phylogenetic relationships, gene duplications, possible cis-elements on promoter regions and expression profiles, especially in the poly- and mono-embryogenic ovules, and suggested that CsZFP7 is involved in nucellar embryogenesis.

Advances in the Study of the Transcriptional Regulation Mechanism of Plant miRNAs

MicroRNAs (miRNA) are a class of endogenous, non-coding, small RNAs with about 22 nucleotides (nt), that are widespread in plants and are involved in various biological processes, such as development, flowering phase transition, hormone signal transduction, and stress response. The transcriptional regulation of miRNAs is an important process of miRNA gene regulation, and it is essential for miRNA biosynthesis and function. Like mRNAs, miRNAs are transcribed by RNA polymerase II, and these transcription processes are regulated by various transcription factors and other proteins. Consequently, the upstream genes regulating miRNA transcription, their specific expression, and the regulating mechanism were reviewed to provide more information for further research on the miRNA regulatory mechanism and help to further understand the regulatory networks of plant miRNAs.

Spatiotemporal profiles of gene activity in stamen delineate nucleo-cytoplasmic interaction in a male-sterile somatic cybrid citrus

Cytoplasmic male sterility (CMS) has long been used to produce seedless fruits in perennial woody crops like citrus. A male-sterile somatic cybrid citrus (G1 + HBP) was generated by protoplast fusion between a CMS callus parent 'Guoqing No. 1' Satsuma mandarin (Citrus unshiu, G1) and a fertile mesophyll parent Hirado Buntan pummelo (Citrus grandis, HBP). To uncover the male-sterile mechanism of G1 + HBP, we compared the transcriptome profiles of stamen organ and cell types at five stages between G1 + HBP and HBP, including the initial stamen primordia, enlarged stamen primordia, pollen mother cells, tetrads, and microspores captured by laser microdissection. The stamen organ and cell types showed distinct gene expression profiles. A majority of genes involved in stamen development were differentially expressed, especially CgAP3.2, which was downregulated in enlarged stamen primordia and upregulated in tetrads of G1 + HBP compared with HBP. Jasmonic acid- and auxin-related biological processes were enriched among the differentially expressed genes of stamen primordia, and the content of jasmonic acid biosynthesis metabolites was higher in flower buds and anthers of G1 + HBP. In contrast, the content of auxin biosynthesis metabolites was lower in G1 + HBP. The mitochondrial tricarboxylic acid cycle and oxidative phosphorylation processes were enriched among the differentially expressed genes in stamen primordia, meiocytes, and microspores, indicating the dysfunction of mitochondria in stamen organ and cell types of G1 + HBP. Taken together, the results indicate that malfunction of mitochondria-nuclear interaction might cause disorder in stamen development, and thus lead to male sterility in the citrus cybrid.

Medicago truncatula PHO2 genes have distinct roles in phosphorus homeostasis and symbiotic nitrogen fixation

Three PHO2-like genes encoding putative ubiquitin-conjugating E2 enzymes of Medicago truncatula were characterized for potential roles in phosphorous (P) homeostasis and symbiotic nitrogen fixation (SNF). All three genes, MtPHO2A, B and C, contain miR399-binding sites characteristic of PHO2 genes in other plant species. Distinct spatiotemporal expression patterns and responsiveness of gene expression to P- and N-deprivation in roots and shoots indicated potential roles, especially for MtPHO2B, in P and N homeostasis. Phenotypic analysis of pho2 mutants revealed that MtPHO2B is integral to Pi homeostasis, affecting Pi allocation during plant growth under nutrient-replete conditions, while MtPHO2C had a limited role in controlling Pi homeostasis. Genetic analysis also revealed a connection between Pi allocation, plant growth and SNF performance. Under N-limited, SNF conditions, Pi allocation to different organs was dependent on MtPHO2B and, to a lesser extent, MtPHO2C and MtPHO2A. MtPHO2A also affected Pi homeostasis associated with nodule formation. Thus, MtPHO2 genes play roles in systemic and localized, i.e., nodule, P homeostasis affecting SNF.

Genomic survey of high-throughput RNA-Seq data implicates involvement of long intergenic non-coding RNAs (lincRNAs) in cytoplasmic male-sterility and fertility restoration in pigeon pea

BACKGROUND

Long-intergenic non-coding RNAs (lincRNAs) originate from intergenic regions and have no coding potential. LincRNAs have emerged as key players in the regulation of various biological processes in plant development. Cytoplasmic male-sterility (CMS) in association with restorer-of-fertility (Rf) systems makes it a highly reliable tool for exploring heterosis for producing commercial hybrid seeds. To date, there have been no reports of lincRNAs during pollen development in CMS and fertility restorer lines in pigeon pea.

OBJECTIVE

Identification of lincRNAs in the floral buds of cytoplasmic male-sterile (AKCMS11) and fertility restorer (AKPR303) pigeon pea lines.

METHODS

We employed a computational approach to identify lincRNAs in the floral buds of cytoplasmic male-sterile (AKCMS11) and fertility restorer (AKPR303) pigeon pea lines using RNA-Seq data.

RESULTS

We predicted a total of 2145 potential lincRNAs of which 966 were observed to be differentially expressed between the sterile and fertile pollen. We identified, 927 cis-regulated and 383 trans-regulated target genes of the lincRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the target genes revealed that these genes were specifically enriched in pathways like pollen and pollen tube development, oxidative phosphorylation, etc. We detected 23 lincRNAs that were co-expressed with 17 pollen-related genes with known functions. Fifty-nine lincRNAs were predicted to be endogenous target mimics (eTMs) for 25 miRNAs, and found to be associated with pollen development. The, lincRNA regulatory networks revealed that different lincRNA-miRNA-mRNA networks might be associated with CMS and fertility restoration.

CONCLUSION

Thus, this study provides valuable information by highlighting the functions of lincRNAs as regulators during pollen development in pigeon pea and utilization in hybrid seed production.

Plant phosphate status influences root biotic interactions

Phosphorus (P) deficiency stress in combination with biotic stress(es) severely impacts crop yield. Plant responses to P deficiency overlapping with that of other stresses exhibit a high degree of complexity involving different signaling pathways. On the one hand, plants engage with rhizosphere microbiome/arbuscular mycorrhizal fungi for improved phosphate (Pi) acquisition and plant stress response upon Pi deficiency; on the other hand, this association is gets disturbed under Pi sufficiency. This nutrient-dependent response is highly regulated by the phosphate starvation response (PSR) mediated by the master regulator, PHR1, and its homolog, PHL. It is interesting to note that Pi status (deficiency/sufficiency) has a varying response (positive/negative) to different biotic encounters (beneficial microbes/opportunistic pathogens/insect herbivory) through a coupled PSR-PHR1 immune system. This also involves crosstalk among multiple players including transcription factors, defense hormones, miRNAs, and Pi transporters, among others influencing the plant-biotic-phosphate interactions. We provide a comprehensive view of these key players involved in maintaining a delicate balance between Pi homeostasis and plant immunity. Finally, we propose strategies to utilize this information to improve crop resilience to Pi deficiency in combination with biotic stresses.

Integrated transcriptome and proteome analysis provides insights into the mechanism of cytoplasmic male sterility (CMS) in tobacco (Nicotiana tabacum L.)

Cytoplasmic male sterility (CMS) is critical in maximizing crop yield and quality by utilizing tobacco heterosis. However, the mechanism of tobacco CMS formation remains unknown. Using paraffin section observation, transcriptome sequencing, and TMT proteomic analysis, this study describes the differences in expression profiles in morphology, transcription, and translation between the sua-CMS tobacco line (MSYY87) and its corresponding maintainer line (YY87). According to the microspore morphology, MSYY87 began to exhibit abnormal microspore development during the early stages of germination and differentiation (androgynous primordium differentiation stage). According to transcriptomic and proteomic analyses, 17 genes/proteins involved in lipid transport/binding and phenylpropane metabolism were significantly down-regulated at both the mRNA and protein levels. Through further analysis, we identified some key genes that may be involved in tobacco male sterility, including β-GLU related to energy metabolism, 4CL and bHLHs related to anther wall formation, nsLTPs related to pollen germination and anther cuticle, and bHLHs related to pollen tapetum degradation. We speculate that the down-regulation of these genes affects the normal physiological metabolism, making tobacco plants show male sterility. SIGNIFICANCE: Cytoplasmic male sterility (CMS) plays a vital role in utilizing tobacco heterosis and enhancing crop yield and quality. We observed paraffin sections and conducted transcriptome sequencing and mitochondrial proteomics to examine the tobacco CMS line Yunyan 87 (MSYY87) and its maintainer line Yunyan 87 (YY87). The down-regulation expression of β-GLU resulted in insufficient ATP supply, which resulted in disordered energy metabolism. The down-regulation expression of 4CL, nsLTPs and bHLHs may affect the formation of anther wall and anther cuticle, pollen germination, as well as the degradation of pollen tapetum. These various abnormal physiological processes, the male sterility of tobacco is finally caused. The findings shed light on the molecular mechanisms of tobacco CMS and serve as a model for fertility research in other flowering plants.

Differential Expression of miRNAs Involved in Response to Candidatus Liberibacter asiaticus Infection in Mexican Lime at Early and Late Stages of Huanglongbing Disease

Huanglongbing (HLB) is one of the most destructive diseases threatening citriculture worldwide. This disease has been associated with α-proteobacteria species, namely Candidatus Liberibacter. Due to the unculturable nature of the causal agent, it has been difficult to mitigate the disease, and nowadays a cure is not available. MicroRNAs (miRNAs) are key regulators of gene expression, playing an essential role in abiotic and biotic stress in plants including antibacterial responses. However, knowledge derived from non-model systems including Candidatus Liberibacter asiaticus (CLas)-citrus pathosystem remains largely unknown. In this study, small RNA profiles from Mexican lime (Citrus aurantifolia) plants infected with CLas at asymptomatic and symptomatic stages were generated by sRNA-Seq, and miRNAs were obtained with ShortStack software. A total of 46 miRNAs, including 29 known miRNAs and 17 novel miRNAs, were identified in Mexican lime. Among them, six miRNAs were deregulated in the asymptomatic stage, highlighting the up regulation of two new miRNAs. Meanwhile, eight miRNAs were differentially expressed in the symptomatic stage of the disease. The target genes of miRNAs were related to protein modification, transcription factors, and enzyme-coding genes. Our results provide new insights into miRNA-mediated regulation in C. aurantifolia in response to CLas infection. This information will be useful to understand molecular mechanisms behind the defense and pathogenesis of HLB.

Overview of Repressive miRNA Regulation by Short Tandem Target Mimic (STTM): Applications and Impact on Plant Biology

The application of miRNA mimic technology for silencing mature miRNA began in 2007. This technique originated from the discovery of the INDUCED BY PHOSPHATE STARVATION 1 (IPS1) gene, which was found to be a competitive mimic that prevents the cleavage of the targeted mRNA by miRNA inhibition at the post-transcriptional level. To date, various studies have been conducted to understand the molecular mimic mechanism and to improve the efficiency of this technology. As a result, several mimic tools have been developed: target mimicry (TM), short tandem target mimic (STTM), and molecular sponges (SPs). STTM is the most-developed tool due to its stability and effectiveness in decoying miRNA. This review discusses the application of STTM technology on the loss-of-function studies of miRNA and members from diverse plant species. A modified STTM approach for studying the function of miRNA with spatial-temporal expression under the control of specific promoters is further explored. STTM technology will enhance our understanding of the miRNA activity in plant-tissue-specific development and stress responses for applications in improving plant traits via miRNA regulation.

miR156 regulates somatic embryogenesis by modulating starch accumulation in citrus

Somatic embryogenesis (SE) is a major regeneration approach for in vitro cultured tissues of plants, including citrus. However, SE capability is difficult to maintain, and recalcitrance to SE has become a major obstacle to plant biotechnology. We previously reported that miR156-SPL modules regulate SE in citrus callus. However, the downstream regulatory pathway of the miR156-SPL module in SE remains unclear. In this study, we found that transcription factors CsAGL15 and CsFUS3 bind to the CsMIR156A promoter and activate its expression. Suppression of csi-miR156a function leads to up-regulation of four target genes, SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (CsSPL) genes, and reduction of SE efficiency. In the short tandem target mimic (STTM)-miR156a overexpression callus (MIM156), the number of amyloplasts and starch content were significantly reduced, and genes involved in starch synthesis and transport were down-regulated. csi-miR172d was down-regulated, whereas the target genes, CsTOE1.1 and CsTOE1.2, which inhibit the expression of starch biosynthesis genes, were up-regulated. In our working model, CsAGL15 and CsFUS3 activate csi-miR156a, which represses CsSPLs and further regulates csi-miR172d and CsTOEs, thus altering starch accumulation in callus cells and regulating SE in citrus. This study elucidates the pathway of miR156-SPLs and miR172-TOEs-mediated regulation of SE, and provides new insights into enhancing SE capability in citrus.

Integrated mRNA and microRNA expression analysis of root response to phosphate deficiency in Medicago sativa

The deficiency of available phosphate significantly limits plant growth and development. This study sought to investigate how alfalfa (Medicago sativa), a high-yielding and high-quality forage widely cultivated worldwide, responds to phosphate deficiency stress by integrating transcriptional and post-transcriptional data. In this study, 6,041 differentially expressed genes (DEGs) were identified in alfalfa roots under phosphate deficiency conditions. Furthermore, psRNATarget, RNAhybrid, and TargetFinder were used to predict the target genes of 137 differentially expressed miRNAs (DEMs) in the root. In total, 3,912 DEGs were predicted as target genes. Pearson correlation analysis revealed 423 pairs of miRNA-mRNA regulatory relationships. MiRNA negatively regulates mRNA involved in regulatory pathways of phosphate deficiency responses in alfalfa. miR156e targeted squamosa promoter-binding-like protein 13A (SPL13), miR160c targeted auxin response factor 18 (ARF18), and miR2587a controlled glycolysis and citrate cycle via Phosphoenolpyruvate carboxykinase (ATP) (PCKA). Novel-miR27 regulated SPX domain-containing protein that controls phosphate transport in alfalfa root, novel-miR3-targeted sulfoquinovosyl transferase SQD2 controlled sulfolipid synthesis and glutathione S-transferase (GST; mediated by miR169j/k and novel-miR159) regulated glutathione metabolism. miR399l regulated auxin-responsive protein SAUR72 involved in IAA signal transduction, while abscisic acid receptor PYL4 (regulated by novel-miR205 and novel-miR83) participated in ABA signal transduction. Combined miRNA-mRNA enrichment analysis showed that most miRNAs regulate the phosphate starvation response of alfalfa by modulating target genes involved in carbohydrate metabolism, sulfolipid metabolism, glutathione metabolism, and hormone signal transduction. Therefore, this study provides new insights into the post-transcriptional regulation mechanism of phosphate deficiency responses and new perspectives on phosphate assimilation pathways in alfalfa and other legumes.

Whole-Transcriptome Analysis Reveals Long Noncoding RNAs Involved in Female Floral Development of Hickory (Carya cathayensis Sarg.)

Hickory, an endemic woody oil and fruit tree species in China, is of great economic value. However, hickory has a long juvenile period and an inconsistent flowering of males and females, thus influencing the bearing rates and further limiting fruits yield. Currently, it is reported that long noncoding RNAs (lncRNAs) play critical regulatory roles in biological processes. However, the role of lncRNAs in the development of hickory female flowers remains unclear. In this study, a total of 6,862 putative lncRNAs were identified from the female flower transcriptomes in three different growth stages of hickory. We proposed that lncRNAs might play an important role in phytohormone signaling processes for flower formation, especially in the abscisic acid and jasmonic acid pathways, according to the results of our Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment. Moreover, we predicted the interactions among four microRNAs (miRNAs), three lncRNAs, and four genes. We proposed that facing the changing environment, LNC_002115 competes with PHOSPHATE2 (PHO2) for the binding sites on cca-miR399f, and protects PHO2 from suppression. In addition, cis-acting LNC_002115 regulates the expression of the SHORT VEGETATIVE PHASE (SVP) by influencing ABRE-binding factor (ABF). In brief, LNC_002115 regulates hickory female floral development by impacting both PHO2 and SVP. This study was the first to identify lncRNAs involved in hickory female floral development, and provided new insight to elucidate how lncRNAs and their targets play a role in female floral development in hickory, thus unfolding the opportunities for functional characterization of blossom-related lncRNAs in further studies.

Regulatory role of microRNAs (miRNAs) in the recent development of abiotic stress tolerance of plants

MicroRNAs (miRNAs) are a distinct groups of single-stranded non-coding, tiny regulatory RNAs approximately 20-24 nucleotides in length. miRNAs negatively influence gene expression at the post-transcriptional level and have evolved considerably in the development of abiotic stress tolerance in a number of model plants and economically important crop species. The present review aims to deliver the information on miRNA-mediated regulation of the expression of major genes or Transcription Factors (TFs), as well as genetic and regulatory pathways. Also, the information on adaptive mechanisms involved in plant abiotic stress responses, prediction, and validation of targets, computational tools, and databases available for plant miRNAs, specifically focus on their exploration for engineering abiotic stress tolerance in plants. The regulatory function of miRNAs in plant growth, development, and abiotic stresses consider in this review, which uses high-throughput sequencing (HTS) technologies to generate large-scale libraries of small RNAs (sRNAs) for conventional screening of known and novel abiotic stress-responsive miRNAs adds complexity to regulatory networks in plants. The discoveries of miRNA-mediated tolerance to multiple abiotic stresses, including salinity, drought, cold, heat stress, nutritional deficiency, UV-radiation, oxidative stress, hypoxia, and heavy metal toxicity, are highlighted and discussed in this review.

Advances in the regulation of plant salt-stress tolerance by miRNA

Salt stress significantly affects the growth, development, yield, and quality of plants. MicroRNAs (miRNAs) are involved in various stress responses via target gene regulation. Their role in regulating salt stress has also received significant attention from researchers. Various transcription factor families are the common target genes of plant miRNAs. Thus, regulating the expression of miRNAs is a novel method for developing salt-tolerant crops. This review summarizes plant miRNAs that mediate salt tolerance, specifically miRNAs that have been utilized in genetic engineering to modify plant salinity tolerance. The molecular mechanism by which miRNAs mediate salt stress tolerance merits elucidation, and this knowledge will promote the development of miRNA-mediated salt-tolerant crops and provide new strategies against increasingly severe soil salinization.

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.

MicroRNA Mediated Plant Responses to Nutrient Stress

To complete their life cycles, plants require several minerals that are found in soil. Plant growth and development can be affected by nutrient shortages or high nutrient availability. Several adaptations and evolutionary changes have enabled plants to cope with inappropriate growth conditions and low or high nutrient levels. MicroRNAs (miRNAs) have been recognized for transcript cleavage and translational reduction, and can be used for post-transcriptional regulation. Aside from regulating plant growth and development, miRNAs play a crucial role in regulating plant’s adaptations to adverse environmental conditions. Additionally, miRNAs are involved in plants’ sensory functions, nutrient uptake, long-distance root transport, and physiological functions related to nutrients. It may be possible to develop crops that can be cultivated in soils that are either deficient in nutrients or have extreme nutrient supplies by understanding how plant miRNAs are associated with nutrient stress. In this review, an overview is presented regarding recent advances in the understanding of plants’ responses to nitrogen, phosphorus, potassium, sulfur, copper, iron, boron, magnesium, manganese, zinc, and calcium deficiencies via miRNA regulation. We conclude with future research directions emphasizing the modification of crops for improving future food security.

tae-miR399-UBC24 Module Enhances Freezing Tolerance in Winter Wheat via a CBF Signaling Pathway

Although the regulation of Pi homeostasis by miR399 has been studied in various plant species, its underlying molecular mechanism in response to freezing stress is still poorly understood. In this work, we found that the expression of tae-miR399 and its target gene TaUBC24 in the tillering nodes of the strong cold-resistant winter wheat cultivar Dongnongdongmai1 (Dn1) was not only significantly altered after severe winters but also responsive to short-term freezing stress. TaUBC24 physically interacted with TaICE1. Enhanced freezing tolerance was observed for tae-miR399-overexpressing Arabidopsis lines. Under freezing stress, overexpression of tae-miR399 ultimately decreased the expression of AtUBC24, inhibiting the degradation of AtICE1, which increased the expression of genes involved in the CBF signaling pathway and starch metabolism and promoted the activities of antioxidant enzymes. These results will improve our understanding of the molecular mechanism through which the miR399-UBC24 module plays a cardinal role in regulating plant freezing stress tolerance through mediation of downstream pathways.

Transcriptomes and DNA methylomes in apomictic cells delineate nucellar embryogenesis initiation in citrus

Citrus nucellar poly-embryony (NPE) is a mode of sporophytic apomixis that asexual embryos formed in the seed through adventitious embryogenesis from the somatic nucellar cells. NPE allows clonal propagation of rootstocks, but it impedes citrus cross breeding. To understand the cellular processes involved in NPE initiation, we profiled the transcriptomes and DNA methylomes in laser microdissection captured citrus apomictic cells. In apomictic cells, ribosome biogenesis and protein degradation were activated, whereas auxin polar transport was repressed. Reactive oxygen species (ROS) accumulated in the poly-embryonic ovules, and response to oxidative stress was provoked. The global DNA methylation level, especially that of CHH context, was decreased, whereas the methylation level of the NPE-controlling key gene CitRWP was increased. A C2H2 domain-containing transcription factor gene and CitRWP co-expressed specifically in apomictic cells may coordinate to initiate NPE. The activated embryogenic development and callose deposition processes indicated embryogenic fate of nucellar embryo initial (NEI) cells. In our working model for citrus NPE initiation, DNA hyper-methylation may activate transcription of CitRWP, which increases C2H2 expression and ROS accumulation, triggers epigenetic regulation and regulates cell fate transition and NEI cell identity in the apomictic cells.

ZmFAR1 and ZmABCG26 Regulated by microRNA Are Essential for Lipid Metabolism in Maize Anther

The function and regulation of lipid metabolic genes are essential for plant male reproduction. However, expression regulation of lipid metabolic genic male sterility (GMS) genes by noncoding RNAs is largely unclear. Here, we systematically predicted the microRNA regulators of 34 maize white brown complex members in ATP-binding cassette transporter G subfamily (WBC/ABCG) genes using transcriptome analysis. Results indicate that the ZmABCG26 transcript was predicted to be targeted by zma-miR164h-5p, and their expression levels were negatively correlated in maize B73 and Oh43 genetic backgrounds based on both transcriptome data and qRT-PCR experiments. CRISPR/Cas9-induced gene mutagenesis was performed on ZmABCG26 and another lipid metabolic gene, ZmFAR1. DNA sequencing, phenotypic, and cytological observations demonstrated that both ZmABCG26 and ZmFAR1 are GMS genes in maize. Notably, ZmABCG26 proteins are localized in the endoplasmic reticulum (ER), chloroplast/plastid, and plasma membrane. Furthermore, ZmFAR1 shows catalytic activities to three CoA substrates in vitro with the activity order of C12:0-CoA > C16:0-CoA > C18:0-CoA, and its four key amino acid sites were critical to its catalytic activities. Lipidomics analysis revealed decreased cutin amounts and increased wax contents in anthers of both zmabcg26 and zmfar1 GMS mutants. A more detailed analysis exhibited differential changes in 54 monomer contents between wild type and mutants, as well as between zmabcg26 and zmfar1. These findings will promote a deeper understanding of miRNA-regulated lipid metabolic genes and the functional diversity of lipid metabolic genes, contributing to lipid biosynthesis in maize anthers. Additionally, cosegregating molecular markers for ZmABCG26 and ZmFAR1 were developed to facilitate the breeding of male sterile lines.

Small RNAs, Degradome, and Transcriptome Sequencing Provide Insights into Papaya Fruit Ripening Regulated by 1-MCP

As an inhibitor of ethylene receptors, 1-methylcyclopropene (1-MCP) can delay the ripening of papaya. However, improper 1-MCP treatment will cause a rubbery texture in papaya. Understanding of the underlying mechanism is still lacking. In the present work, a comparative sRNA analysis was conducted after different 1-MCP treatments and identified a total of 213 miRNAs, of which 44 were known miRNAs and 169 were novel miRNAs in papaya. Comprehensive functional enrichment analysis indicated that plant hormone signal pathways play an important role in fruit ripening. Through the comparative analysis of sRNAs and transcriptome sequencing, a total of 11 miRNAs and 12 target genes were associated with the ethylene and auxin signaling pathways. A total of 1741 target genes of miRNAs were identified by degradome sequencing, and nine miRNAs and eight miRNAs were differentially expressed under the ethylene and auxin signaling pathways, respectively. The network regulation diagram of miRNAs and target genes during fruit ripening was drawn. The expression of 11 miRNAs and 12 target genes was verified by RT-qPCR. The target gene verification showed that cpa-miR390a and cpa-miR396 target CpARF19-like and CpERF RAP2-12-like, respectively, affecting the ethylene and auxin signaling pathways and, therefore, papaya ripening.

MicroRNAs Involved in Regulatory Cytoplasmic Male Sterility by Analysis RNA-seq and Small RNA-seq in Soybean

Cytoplasmic male sterility (CMS) is an important plant characteristic for exploiting heterosis to enhance crop traits during breeding. However, the CMS regulatory network remains unclear in plants, even though researchers have attempted to isolate genes associated with CMS. In this study, we performed high-throughput sequencing and degradome analyses to identify microRNAs (miRNAs) and their targets in a soybean CMS line (JLCMS9A) and its maintainer line (JLCMS9B). Additionally, the differentially expressed genes during reproductive development were identified using RNA-seq data. A total of 280 miRNAs matched soybean miRNA sequences in miRBase, including mature miRNAs and pre-miRNAs. Of the 280 miRNAs, 30, 23, and 21 belonged to the miR166, miR156, and miR171 families, respectively. Moreover, 410 novel low-abundant miRNAs were identified in the JLCMS9A and JLCMS9B flower buds. Furthermore, 303 and 462 target genes unique to JLCMS9A and JLCMS9B, respectively, as well as 782 common targets were predicted based on the degradome analysis. Target genes differentially expressed between the CMS line and the maintainer line were revealed by an RNA-seq analysis. Moreover, all target genes were annotated with diverse functions related to biological processes, cellular components, and molecular functions, including transcriptional regulation, the nucleus, meristem maintenance, meristem initiation, cell differentiation, auxin-activated signaling, plant ovule development, and anther development. Finally, a network was built based on the interactions. Analyses of the miRNA, degradome, and transcriptome datasets generated in this study provided a comprehensive overview of the reproductive development of a CMS soybean line. The data presented herein represent useful information for soybean hybrid breeding. Furthermore, the study results indicate that miRNAs might contribute to the soybean CMS regulatory network by modulating the expression of CMS-related genes. These findings lay the foundation for future studies on the molecular mechanisms underlying soybean CMS.

Heat-Responsive miRNAs Participate in the Regulation of Male Fertility Stability in Soybean CMS-Based F1 under High Temperature Stress

MicroRNAs (miRNAs), a class of noncoding small RNAs (sRNAs), are widely involved in the response to high temperature (HT) stress at both the seedling and flowering stages. To dissect the roles of miRNAs in regulating male fertility in soybean cytoplasmic male sterility (CMS)-based F1 under HT, sRNA sequencing was performed using flower buds from HT-tolerant and HT-sensitive CMS-based F1 combinations (NF1 and YF1, respectively). A total of 554 known miRNAs, 59 new members of known miRNAs, 712 novel miRNAs, and 1145 target genes of 580 differentially expressed miRNAs (DEMs) were identified under normal temperature and HT conditions. Further integrated analysis of sRNA and transcriptome sequencing found that 21 DEMs and 15 differentially expressed target genes, such as gma-miR397a/Laccase 2, gma-miR399a/Inorganic phosphate transporter 1-4, and gma-miR4413a/PPR proteins, mitochondrial-like, were negatively regulated under HT stress. Furthermore, all members of the gma-miR156 family were suppressed by HT stress in both NF1 and YF1, but were highly expressed in YF1 under HT condition. The negative correlation between gma-miR156b and its target gene squamosa promoter-binding protein-like 2b was confirmed by expression analysis, and overexpression of gma-miR156b in Arabidopsis led to male sterility under HT stress. With these results, we proposed that miRNAs play an important role in the regulation of male fertility stability in soybean CMS-based F1 under HT stress.

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

A photoperiod- and thermo-sensitive genic male sterile (PTGMS) line is the basic material for two-hybrid rice and is an important genetic breeding resource. Peiai64S (PA64S) is an important germplasm resource of PTGMS rice, and it has been applied to two-line hybrid rice systems in China. Pollen fertility in PA64S is regulated by the temperature and photoperiod, but the mechanism of the fertility transition is unclear. In this study, we obtained the male fertile plant PA64S(F) and the male sterile plant PA64S(S) by controlling different temperatures under long light conditions and used the male fertile and sterile plants to investigate the role of microRNAs (miRNAs) in regulating male fertility in rice. We performed the small RNA library sequencing of anthers from PA64S(S) and PA64S(F). A total of 196 miRNAs were identified-166 known miRNAs among 27 miRNA families and 30 novel miRNAs. In the transcriptome analysis, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of differentially expressed genes revealed significant enrichment in the synthesis and metabolism of fatty acids and some secondary metabolism pathways such as fatty acid metabolism and phenylalanine metabolism. With a comprehensive analysis of miRNA, transcriptome, and degradome sequencing, we identified that 13 pairs of miRNA/target genes regulated male fertility in rice by responding to temperature change, among which the miR156, miR5488, and miR399 affect the male fertility of PA64S by influencing SPLs, the lignin synthesis of anther walls, and the flavonoid metabolism pathway. The results provide a new understanding of PTGMS rice, which will help us better understand the potential regulatory mechanisms of male sterility in the future.