Integrated analysis of miRNA-mRNA regulatory networks of potato (Solanum tuberosum L.) in response to cadmium stress

Nano-Fe3O4: Enhancing the tolerance of Elymus nutans to Cd stress through regulating programmed cell death

Cadmium (Cd) poses a significant threat to plant growth and the environment. Nano-Fe3O4 is effective in alleviating Cd stress in plants. Elymus nutans Griseb. is an important fodder crop on the Qinghai-Tibetan Plateau (QTP). However, the potential mechanism by which nano-Fe3O4 alleviates Cd stress in E. nutans is not well understood. E. nutans were subjected to single Cd, single nano-Fe3O4, and co-treatment with nano-Fe3O4 and Cd, and the effects on morphology, Cd uptake, antioxidant enzyme activity, reactive oxygen species (ROS) levels and programmed cell death (PCD) were studied to clarify the regulatory mechanism of nano-Fe3O4. The results showed that Cd stress significantly decreased the germination percentage and biomass of E. nutans. The photosynthetic pigment content decreased significantly under Cd stress. Cd stress also caused oxidative stress and lipid peroxidation, accumulation of excessive ROS, resulting in PCD, but the effect of nano-Fe3O4 was different. Seed germination, seedling growth, and physiological processes were analyzed to elucidate the regulatory role of nano-Fe3O4 nanoparticles in promoting photosynthesis, reducing Cd accumulation, scavenging ROS, and regulating PCD, to promote seed germination and seedling growth in E. nutans. This report provides a scientific basis for improving the tolerance of Elymus to Cd stress by using nano-Fe3O4.

Integration analysis of miRNA-mRNA pairs between two contrasting genotypes reveals the molecular mechanism of jujube (Ziziphus jujuba Mill.) response to high-temperature stress

With global warming, high temperature (HT) has become one of the most common abiotic stresses resulting in significant crop yield losses, especially for jujube (Ziziphus jujuba Mill.), an important temperate economic crop cultivated worldwide. This study aims to explore the coping mechanism of jujube to HT stress at the transcriptional and post-transcriptional levels, including identifying differentially expressed miRNAs and mRNAs as well as elucidating the critical pathways involved. High-throughput sequencing analyses of miRNA and mRNA were performed on jujube leaves, which were collected from "Fucumi" (heat-tolerant) and "Junzao" (heat-sensitive) cultivars subjected to HT stress (42 °C) for 0, 1, 3, 5, and 7 days, respectively. The results showed that 45 known miRNAs, 482 novel miRNAs, and 13,884 differentially expressed mRNAs (DEMs) were identified. Among them, integrated analysis of miRNA target genes prediction and mRNA-seq obtained 1306 differentially expressed miRNAs-mRNAs pairs, including 484, 769, and 865 DEMIs-DEMs pairs discovered in "Fucuimi", "Junzao" and two genotypes comparative groups, respectively. Furthermore, functional enrichment analysis of 1306 DEMs revealed that plant-pathogen interaction, starch and sucrose metabolism, spliceosome, and plant hormone signal transduction were crucial pathways in jujube leaves response to HT stress. The constructed miRNA-mRNA network, composed of 20 DEMIs and 33 DEMs, displayed significant differently expressions between these two genotypes. This study further proved the regulatory role of miRNAs in the response to HT stress in plants and will provide a theoretical foundation for the innovation and cultivation of heat-tolerant varieties.

Effect of exogenous γ-aminobutyric acid on physiological property, antioxidant activity, and cadmium uptake of quinoa seedlings under cadmium stress

Increasing cadmium (Cd) pollution has negative effects on quinoa growth and production. Gamma-aminobutyric acid (GABA) confers plants with stress resistance to heavy metals; however, the mechanism remains unclear. We explored the effects of exogenous GABA on the physiological characteristics, antioxidant capacity, and Cd accumulation of quinoa seedlings under Cd stress using hydroponic experiments. Partial least-squares regression was used to identify key physical and chemical indices of seedlings affecting Cd accumulation. Compared with those of the CK group, exposure to 10 and 25 µmol·L-1 Cd significantly reduced the photosynthetic pigment contents, photosynthesis, and biomass accumulation of quinoa seedlings; resulted in shorter and thicker roots; decreased the length of the lateral roots; decreased the activities of superoxide dismutase (SOD) and peroxide (POD); and increased H2O2 and malondialdehyde (MDA) contents. Exogenous GABA reduced the Cd content in the stem/leaves and roots of quinoa seedlings under Cd stress by 13.22-21.63% and 7.92-28.32%, decreased Cd accumulation by 5.37-6.71% and 1.91-4.09%, decreased the H2O2 content by 38.21-47.46% and 45.81-55.73%, and decreased the MDA content by 37.65-48.12% and 29.87-32.51%, respectively. GABA addition increased the SOD and POD activities in the roots by 2.78-5.61% and 13.81-18.33%, respectively, under Cd stress. Thus, exogenous GABA can reduce the content and accumulation of Cd in quinoa seedlings by improving the photosynthetic characteristics and antioxidant enzyme activity and reducing the degree of lipid peroxidation in the cell membrane to alleviate the toxic effect of Cd stress on seedling growth.

Multiomics and biotechnologies for understanding and influencing cadmium accumulation and stress response in plants

Cadmium (Cd) is one of the most toxic heavy metals faced by plants and, additionally, via the food chain, threatens human health. It is principally dispersed through agro-ecosystems via anthropogenic activities and geogenic sources. Given its high mobility and persistence, Cd, although not required, can be readily assimilated by plants thereby posing a threat to plant growth and productivity as well as animal and human health. Thus, breeding crop plants in which the edible parts contain low to zero Cd as safe food stuffs and harvesting shoots of high Cd-containing plants as a route for decontaminating soils are vital strategies to cope with this problem. Recently, multiomics approaches have been employed to considerably enhance our understanding of the mechanisms underlying (i) Cd toxicity, (ii) Cd accumulation, (iii) Cd detoxification and (iv) Cd acquisition tolerance in plants. This information can be deployed in the development of the biotechnological tools for developing plants with modulated Cd tolerance and detoxification to safeguard cellular and genetic integrity as well as to minimize food chain contamination. The aim of this review is to provide a current update about the mechanisms involved in Cd uptake by plants and the recent developments in the area of multiomics approach in terms of Cd stress responses, as well as in the development of Cd tolerant and low Cd accumulating crops.

Gibberellin Positively Regulates Tomato Resistance to Tomato Yellow Leaf Curl Virus (TYLCV)

Tomato yellow leaf curl virus (TYLCV) is a prominent viral pathogen that adversely affects tomato plants. Effective strategies for mitigating the impact of TYLCV include isolating tomato plants from the whitefly, which is the vector of the virus, and utilizing transgenic lines that are resistant to the virus. In our preliminary investigations, we observed that the use of growth retardants increased the rate of TYLCV infection and intensified the damage to the tomato plants, suggesting a potential involvement of gibberellic acid (GA) in the conferring of resistance to TYLCV. In this study, we employed an infectious clone of TYLCV to inoculate tomato plants, which resulted in leaf curling and growth inhibition. Remarkably, this inoculation also led to the accumulation of GA3 and several other phytohormones. Subsequent treatment with GA3 effectively alleviated the TYLCV-induced leaf curling and growth inhibition, reduced TYLCV abundance in the leaves, enhanced the activity of antioxidant enzymes, and lowered the reactive oxygen species (ROS) levels in the leaves. Conversely, the treatment with PP333 exacerbated TYLCV-induced leaf curling and growth suppression, increased TYLCV abundance, decreased antioxidant enzyme activity, and elevated ROS levels in the leaves. The analysis of the gene expression profiles revealed that GA3 up-regulated the genes associated with disease resistance, such as WRKYs, NACs, MYBs, Cyt P450s, and ERFs, while it down-regulated the DELLA protein, a key agent in GA signaling. In contrast, PP333 induced gene expression changes that were the opposite of those caused by the GA3 treatment. These findings suggest that GA plays an essential role in the tomato's defense response against TYLCV and acts as a positive regulator of ROS scavenging and the expression of resistance-related genes.

Combined transcriptomics and metabolomics to analyse the response of Cuminum cyminum L. under Pb stress

Lead (Pb) can disrupt plant gene expression, modify metabolite contents, and influence the growth of plants. Cuminum cyminum L. is highly adaptable to adversity, but molecular mechanism by which it responds to Pb stress is unknown. For this study, transcriptomic and metabolomic sequencing was performed on root tissues of C. cyminum under Pb stress. Our results showed that high Pb stress increased the activity of peroxidase (POD), the contents of malondialdehyde (MDA) and proline by 80.03 %, 174.46 % and 71.24 %, respectively. Meanwhile, Pb stress decreased the activities of superoxide dismutase (SOD) and catalase (CAT) as well as contents of soluble sugars and GSH, which thus affected the growth of C. cyminum. In addition, Pb stress influenced the accumulation and transport of Pb in C. cyminum. Metabolomic results showed that Pb stress affected eight metabolic pathways involving 108 differentially expressed metabolites, primarily amino acids, organic acids, and carbohydrates. The differentially expressed genes identified through transcriptome analysis were mainly involved the oxidation reductase activity, transmembrane transport, phytohormone signaling, and MAPK signaling pathway. The results of this study will help to understand the molecular mechanisms of C. cyminum response to Pb stress, and provide a basis for screening seeds with strong resistance to heavy metals.

L-type lectin receptor-like kinase OsCORK1 as an important negative regulator confers copper stress tolerance in rice

Copper (Cu) is vital for plant growth but becomes toxic in excess, posing potential threats to human health. Although receptor-like kinases (RLKs) have been studied in plant response to abiotic stresses, their roles in Cu stress response remain poorly understood. Therefore, we aimed to evaluate Cu toxicity effects on rice and elucidate its potential molecular mechanisms. Specifically, rice lectin-type RLK OsCORK1 (Copper-response receptor-like kinase 1) function in Cu stress response was investigated. RNA sequencing and expression assays revealed that OsCORK1 is mainly expressed in roots and leaves, and its expression was significantly induced by Cu stress time- and dose-dependently. Kinase activity assays demonstrated OsCORK1 as a Mn2+-preferred functional kinase. Genetically, OsCORK1 gene-edited mutants exhibited increased tolerance to Cu stress and reduced Cu accumulation compared to the wild type (WT). Conversely, OsCORK1 overexpression compromised the Cu stress tolerance observed in OsCORK1 gene-edited mutants. OsCORK1 gene-edited mutants slightly damaged the root tips compared to the WT under Cu stress. Furthermore, OsCORK1 was demonstrated to modulate Cu stress tolerance by mainly altering cell wall components, particularly lignin, in rice. Overall, OsCORK1 is an important negative regulator of Cu stress tolerance, providing a potential gene target to reduce Cu pollution in rice production.

Unveiling the biosynthesis, mechanisms, and impacts of miRNAs in drought stress resilience in plants

Drought stress is one of the most serious threats to sustainable agriculture and is predicted to be further intensified in the coming decades. Therefore, understanding the mechanism of drought stress tolerance and the development of drought-resilient crops are the major goals at present. In recent years, noncoding microRNAs (miRNAs) have emerged as key regulators of gene expressions under drought stress conditions and are turning out to be the potential candidates that can be targeted to develop drought-resilient crops in the future. miRNAs are known to target and decrease the expression of various genes to govern the drought stress response in plants. In addition, emerging evidence also suggests a regulatory role of long non-coding RNAs (lncRNAs) in the regulation of miRNAs and the expression of their target genes by a process referred as miRNA sponging. In this review, we present the regulatory roles of miRNAs in the modulation of drought-responsive genes along with discussing their biosynthesis and action mechanisms. Additionally, the interactive roles of miRNAs with phytohormone signaling components have also been highlighted to present the global view of miRNA functioning under drought-stress conditions.

An insight into the roles of regulatory ncRNAs in plants: An abiotic stress and developmental perspective

Different environmental cues lead to changes in physiology, biochemistry and molecular status of plant's growth. Till date, various genes have been accounted for their role in regulating plant development and response to abiotic stress. Excluding genes that code for a functional protein in a cell, a large chunk of the eukaryotic transcriptome consists of non-coding RNAs (ncRNAs) which lack protein coding capacity but are still functional. Recent advancements in Next Generation Sequencing (NGS) technology have led to the unearthing of different types of small and large non-coding RNAs in plants. Non-coding RNAs are broadly categorised into housekeeping ncRNAs and regulatory ncRNAs which work at transcriptional, post-transcriptional and epigenetic levels. Diverse ncRNAs play different regulatory roles in nearly all biological processes including growth, development and response to changing environments. This response can be perceived and counteracted by plants using diverse evolutionarily conserved ncRNAs like miRNAs, siRNAs and lncRNAs to participate in complex molecular regimes by activating gene-ncRNA-mRNA regulatory modules to perform the downstream function. Here, we review the current understanding with a focus on recent advancements in the functional studies of the regulatory ncRNAs at the nexus of abiotic stresses and development. Also, the potential roles of ncRNAs in imparting abiotic stress tolerance and yield improvement in crop plants are also discussed with their future prospects.

Multiomics reveals Claroideoglomus etunicatum regulates plant hormone signal transduction, photosynthesis and La compartmentalization in maize to promote growth under La stress

Rare earth elements (REEs) have been widely used in traditional and high-tech fields, and high doses of REEs are considered a risk to the ecosystem. Although the influence of arbuscular mycorrhizal fungi (AMF) in promoting host resistance to heavy metal (HM) stress has been well documented, the molecular mechanism by which AMF symbiosis enhances plant tolerance to REEs is still unclear. A pot experiment was conducted to investigate the molecular mechanism by which the AMF Claroideoglomus etunicatum promotes maize (Zea mays) seedling tolerance to lanthanum (La) stress (100 mg·kg^-1 La). C. etunicatum symbiosis significantly improved maize seedling growth, P and La uptake and photosynthesis. Transcriptome, proteome, and metabolome analyses performed alone and together revealed that differentially expressed genes (DEGs) related to auxin /indole-3-acetic acid (AUX/IAA) and the DEGs and differentially expressed proteins (DEPs) related to ATP-binding cassette (ABC) transporters, natural resistance-associated macrophage proteins (Nramp6), vacuoles and vesicles were upregulated. In contrast, photosynthesis-related DEGs and DEPs were downregulated, and 1-phosphatidyl-1D-myo-inositol 3-phosphate (PI(3)P) was more abundant under C. etunicatum symbiosis. C. etunicatum symbiosis can promote plant growth by increasing P uptake, regulating plant hormone signal transduction, photosynthesis and glycerophospholipid metabolism pathways and enhancing La transport and compartmentalization in vacuoles and vesicles. The results provide new insights into the promotion of plant REE tolerance by AMF symbiosis and the possibility of utilizing AMF-maize interactions in REE phytoremediation and recycling.

microRNAomic profiling of maize root reveals multifaceted mechanisms to cope with Cr (VI) stress

Chromium (Cr) contamination of soil and water poses serious threats to agricultural crop production. MicroRNAs (miRNAs) are conserved, non-coding small RNAs that play pivotal roles in plant growth, development and stress responses through fine-tuning of post-transcriptional gene expression. To better understand the molecular circuit of Cr-responsive miRNAs, two sRNA libraries were prepared from control and Cr (VI) [100 ppm] exposed maize roots. Using deep sequencing, we identified 80 known (1 up and 79 down) and 18 downregulated novel miRNAs from Cr (VI) challenged roots. Gene ontology (GO) analysis reveals that predicted target genes of Cr (VI) responsive miRNAs are potentially involved in diverse cellular and biological processes including plant growth and development (miR159c, miR164d, miR319b-3p and zma_25.145), redox homeostasis (miR528-5p, miR396a-5p and zma_9.132), heavy metal uptake and detoxification (miR159f-5p, 164e-5p, miR408a, miR444f and zma_2.127), signal transduction (miR159f, miR160a-5p, miR393a-5p, miR408-5p and zma_43.158), cell signalling (miR156j, 159c-5p, miR166c-5p and miR398b). Higher accumulation of Cr in maize roots might be due to upregulation of ABC transporter G family member 29 targeted by miR444f. Instead of isolated increase in SOD expression, significant decline in GSH:GSSH ratio and histochemical staining strongly suggest Cr (VI) stress mediated disruption of ROS scavenging machinery thus unbalancing normal cellular homeostasis. Moreover, miR159c-mediated enhanced expression of GAMYB might be a reason for impaired root growth under Cr (VI) stress. In a nutshell, the present microRNAomic study sheds light on the miRNA-target gene regulatory network involved in adaptive responses of maize seedlings to Cr (VI) stress.

Effects of cadmium on transcription, physiology, and ultrastructure of two tobacco cultivars

Cadmium (Cd) is one of the most toxic heavy metal pollutants worldwide. Tobacco is an important cash crop; however, the accumulation of Cd in its biomass is very high. Cadmium may enter the body of smokers with contaminated tobacco and the surrounding environment via smoke. Therefore, it is important to understand the mechanisms of Cd accumulation and tolerance in tobacco plants, especially in the leaves. In this study, the effects of Cd on the growth, accumulation, and biochemical indices of two tobacco varieties, K326 (Cd resistant) and NC55 (Cd sensitive), were studied through transcriptomic and physiological experiments. Transcriptome and physiological analyses showed differences in the expression of Cd transport and Cd resistance related genes between NC55 and K326 under Cd stress. The root meristem cells of NC55 were more severely damaged. The antioxidant enzyme activity, ABA and ZT content, chlorophyll content, photosynthetic rate, and nitrogen metabolism enzyme activity in K326 leaves were higher than in NC55. These data elucidate the mechanisms of low Cd accumulation and high Cd tolerance in K326 leaves and provide a theoretical basis for cultivating tobacco varieties with low Cd accumulation and high Cd resistance.

Integration of miRNA and mRNA analysis reveals the role of ribosome in to anti-artificial aging in sweetcorn

Sweetcorn is a kind of maize with high sugar content and has poor seed aging tolerance, which seriously limits its production. However, few studies have explored the artificial aging (AA) tolerance by miRNA-mRNA integration analysis in sweetcorn. Here, we characterized the physiological, biochemical and transcriptomic changes of two contrasting lines K62 and K107 treated with AA during time series. Both the germination indexes and antioxidant enzymes showed significant difference between two lines. The MDA content of AA-tolerant genotype K62 was significantly lower than that of K107 on the fourth and sixth day. Subsequently, 157 differentially expressed miRNAs (DEMIs) and 8878 differentially expressed mRNAs (DEMs) were identified by RNA-seq analysis under aging stress. The "ribosome" and "peroxisome" pathways were enriched to respond to aging stress, genes for both large units and small ribosomal subunits were significantly upregulated expressed and higher translation efficiency might exist in K62. Thirteen pairs of miRNA-target genes were obtained, and 8 miRNA-mRNA pairs might involve in ribosome protein and translation process. Our results elucidate the mechanism of sweetcorn response to AA at miRNA-mRNA level, and provide a new insight into sweetcorn AA response to stress.

Molecular mechanisms of atrazine toxicity on H19-7 hippocampal neurons revealed by integrated miRNA and mRNA "omics"

Atrazine (ATR) is a widely applied herbicide in Asia and South America with slow natural degradation and documented deleterious effects on human and animal health, including hippocampal toxicity. However, relatively little is known about the molecular mechanisms responsible for ATR-induced hippocampal damage. Screening for differentially expressed mRNAs and microRNAs (miRNAs), and construction of potential miRNA-mRNA regulatory networks can reveal such mechanisms, so we analyzed the mRNA and miRNA expression profiles of rat hippocampus-derived H19-7 cells in response to ATR (500 μM) and conducted Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes enrichment (KEGG) analyses. Integration of miRNA sequencing (miRNA-seq) and mRNA sequencing (mRNA-seq) results identified 114 differentially expressed miRNAs (DEMIs, 40 upregulated and 74 downregulated), and 510 differentially expressed mRNAs (DEMs, 177 upregulated and 333 downregulated) targeted by these DEMIs. The top 10 hub mRNAs (Fos, Prkcb, Ncf1, Vcam1, Atf3, Pak3, Pak1, Cacna1s, Junb, and Ccl2) and 19 related miRNAs (rno-miR-194-5p, rno-miR-24-3p, rno-miR-3074, rno-miR-1949, rno-miR-218a-1-3p, rno-miR-1843a-5p, rno-miR-1843b-5p, rno-miR-296-3p, rno-miR-320-3p, rno-miR-219a-1-3p, rno-miR-122-5p, rno-miR-1839-5p, rno-miR-1843a-3p, rno-miR-215, rno-miR-3583-3p, rno-miR-194-3p, rno-miR-128-1-5p, rno-miR-1956-5p, and rno-miR-466b-2-3p) were validated by quantitative real-time PCR. GO analysis indicated that these DEMs were enriched in genes associated with synaptic plasticity and antioxidant capacity, while KEGG analysis suggested that enriched DEMs were involved in calcium signaling, axon guidance, MAPK signaling, and glial carcinogenesis. The miRNA-mRNA regulatory network identified here may provide potential biomarkers and novel strategies for the treatment of hippocampal neurotoxicity induced by ATR.

Integration of mRNA and microRNA analysis reveals the molecular mechanisms underlying drought stress tolerance in maize (Zea mays L.)

Drought is among the most serious environmental issue globally, and seriously affects the development, growth, and yield of crops. Maize (Zea mays L.), an important crop and industrial raw material, is planted on a large scale worldwide and drought can lead to large-scale reductions in maize corn production; however, few studies have focused on the maize root system mechanisms underlying drought resistance. In this study, miRNA-mRNA analysis was performed to deeply analyze the molecular mechanisms involved in drought response in the maize root system under drought stress. Furthermore, preliminary investigation of the biological function of miR408a in the maize root system was also conducted. The morphological, physiological, and transcriptomic changes in the maize variety "M8186" at the seedling stage under 12% PEG 6000 drought treatment (0, 7, and 24 h) were analyzed. With prolonged drought stress, seedlings gradually withered, the root system grew significantly, and abscisic acid, brassinolide, lignin, glutathione, and trehalose content in the root system gradually increased. Furthermore, peroxidase activity increased, while gibberellic acid and jasmonic acid gradually decreased. Moreover, 32 differentially expressed miRNAs (DEMIRs), namely, 25 known miRNAs and 7 new miRNAs, and 3,765 differentially expressed mRNAs (DEMRs), were identified in maize root under drought stress by miRNA-seq and mRNA-seq analysis, respectively. Through combined miRNA-mRNA analysis, 16 miRNA-target gene pairs, comprising 9 DEMIRs and 15 DEMRs, were obtained. In addition, four metabolic pathways, namely, "plant hormone signal transduction", "phenylpropane biosynthesis", "glutathione metabolism", and "starch and sucrose metabolism", were predicted to have important roles in the response of the maize root system to drought. MiRNA and mRNA expression results were verified by real-time quantitative PCR. Finally, miR408a was selected for functional analysis and demonstrated to be a negative regulator of drought response, mainly through regulation of reactive oxygen species accumulation in the maize root system. This study helps to elaborate the regulatory response mechanisms of the maize root system under drought stress and predicts the biological functions of candidate miRNAs and mRNAs, providing strategies for subsequent mining for, and biological breeding to select for, drought-responsive genes in the maize root system.

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.