Purple acid phosphatases: roles in phosphate utilization and new emerging functions

Extracellular ATP: an emerging multifaceted regulator of plant fitness

Adenosine 5'-triphosphate (ATP) is the energy currency of living organisms and the primary form of organic phosphate (Po) involved in cellular metabolism. In plants, some ATP is released into the extracellular matrix (ECM) in response to various stimuli, where it functions as extracellular ATP (eATP), a key signalling molecule. Recent advances have shed light on the mechanisms of eATP signalling in plants. This review consolidates these findings, beginning with the role of eATP in regulating plant growth, development and responses to biotic and abiotic stresses. It further summarizes the pathways of eATP accumulation and degradation in the ECM and introduces the cellular signalling pathways mediating eATP responses, as reported in key studies. Finally, perspectives on future research directions in this field are presented.

A lipid synthase maintains metabolic flux for jasmonate synthesis to regulate root growth and phosphate homeostasis

Plants require phosphate (Pi) for proper growth and development but often face scarcity of this vital nutrient in the soil. Pi starvation triggers membrane lipid remodeling to utilize the membrane phospholipid-bound Pi in plants. In this process, phospholipids are replaced by non-Pi-containing galactolipids (monogalactosyldiacylglycerol, MGDG; digalactosyldiacylglycerol, DGDG) and sulfolipids. The galactolipids ratio (MGDG:DGDG) is suggested to influence jasmonic acid (JA) biosynthesis. However, how the MGDG:DGDG ratio, JA levels, and root growth are coordinated under Pi deficiency in rice (Oryza sativa) remains unknown. Here, we characterized DGDG synthase 1 (OsDGD1) for its role in regulating root development by maintaining metabolic flux for JA biosynthesis. We showed that OsDGD1 is responsive under low Pi and is under the direct control of Phosphate Starvation Response 2, the master regulator of low Pi adaptations. Further, OsDGD1 knockout (KO) lines showed marked phenotypic differences compared to the wild type, including a significant reduction in root length and biomass, leading to reduced Pi uptake. Further, lipidome analyses revealed reduced DGDG levels in the KO line, leading to reduced membrane remodeling, thus affecting P utilization efficiency. We also observed an increase in the MGDG:DGDG ratio in KO lines, which enhanced the endogenous JA levels and signaling. This imbalance of JA in KO plants led to changes in auxin levels, causing drastic root growth inhibition. These findings indicate the critical role of OsDGD1 in maintaining optimum levels of JA during Pi deficiency for conducive root growth. Besides acting as signaling molecules and structural components, our study widens the role of lipids as metabolic flux controllers for phytohormone biosynthesis.

Identification and Expression Analysis of Acid Phosphatase Gene (PAP) in Brassica napus: Effects of cis-Acting Elements on Two BnaPAP10 Genes in Response to Phosphorus Stress

Purple acid phosphatases (PAPs) play a key role in phosphorus (P) assimilation and redistribution in plants, catalyzing the hydrolysis of phosphate esters to produce inorganic phosphate (Pi). In this study, a total of 77 PAP genes were identified in B. napus. The candidate genes were divided into three groups and ten subgroups based on the phylogenetic analyses and exon-intron organization. Among these 77 BnaPAP proteins, 35 exhibit typical metal-ligating residues characteristic of known PAPs, whereas certain unaltered amino acid residues were absent or displaced in other BnaPAPs. A computational prediction was conducted, revealing that the majority of PAPs contain signal peptide motifs and display a range of N-glycosylation levels, as well as transmembrane helix motifs. An analysis of previously obtained RNA-seq data revealed that 55.84% (43 of 77) of the BnaPAPs responded to Pi deficiency. Moreover, we conducted a preliminary examination of the expression profiles of BnaPAP genes in response to salt stress, and discovered that 42.86% (33 of 77) of these genes were induced under salt stress, either in the shoots or in the roots. Further qRT-PCR and GUS analyses revealed that BnaC9.PAP10 and BnaA7.PAP10, two paralogs of BnaPAP10s, were induced by Pi deficiency. Notably, BnaC9.PAP10 exhibits robust induction, compared to the relatively mild induction observed in BnaA7.PAP10. Our research shows that BnaA7.PAP10 uniquely responds to Pi stress via the W-box, while BnaA7.PAP10 predominantly responds via the P1BS element, and the differences in cis-regulatory elements (CREs) within their promoter regions specifically contribute to their distinct expression levels under Pi stress. Our findings provide valuable insights and establish a foundation for future functional studies of BnaPAPs.

Combined Analysis of the Leaf Metabolome, Lipidome, and Candidate Gene Function: Insights into Genotypic Variation in Phosphorus Utilization Efficiency in Stylosanthes guianensis

Stylo (Stylosanthes guianensis) exhibits excellent tolerance to low-phosphate (Pi) availability, but the underlying mechanisms responsible for improving the phosphorus (P) utilization efficiency (PUE) remain unclear. This study employed metabolomics, lipidomics, and gene expression analyses to investigate the differential responses to low-Pi stress between the high-PUE genotype CF047827 and the cultivar Reyan No. 2. Results showed that CF047827 had higher expression levels of membrane lipid remodeling-related genes in its leaves compared to Reyan No. 2 under low-Pi conditions. This was accompanied by greater phospholipid degradation and non-P-containing lipid biosynthesis in the leaves of CF047827. Furthermore, the purple acid phosphatase gene SgPAP27a, which is more highly expressed in the leaves of CF047827 than in Reyan No. 2 under low-Pi conditions, was identified and functionally characterized. Its role in promoting phospholipid degradation and enhancing PUE was confirmed through heterologous expression in Arabidopsis. These findings provide insights and identify potential candidate genes for breeding high-PUE crop cultivars.

Genome-Wide Analysis, Identification, and Transcriptional Profile of the Response to Abiotic Stress of the Purple Acid Phosphatases (PAP) Gene Family in Apple

Purple acid phosphatases (PAPs) play a significant role in plant phosphorus nutrition and can not only release phosphorus from the soil but also regulate the distribution of phosphorus in plants throughout their entire growth and development process. Moreover, members of the PAP protein family exert a more extensive influence on plant mineral homeostasis, developmental processes, and stress responses. Three clusters of purple acid phosphatases, including 31 putative genes, were identified in apples (Malus domestica) by searching the Genome Database for Rosaceae. The structure, chromosomal distribution and location, phylogeny, motifs, and cis-acting elements in the gene promoter regions of the MdPAP gene family are reviewed. These genes exhibit different expression patterns in different tissues. For example, almost all MdPAP genes are strongly expressed in the roots, except for MdPAP10, MdPAP12, and MdPAP27. Similarily, all MdPAPs were expressed in the leaves while the transcript levels of MdPAP7, MdPAP10, MdPAP15, MdPAP21, MdPAP24, MdPAP26, MdPAP29, and MdPAP30 were highest in apple flowers. Overall, the expression of the 31 genes significantly changed in either the roots or leaves following the application of phosphorus and/or drought stress. These results indicate that MdPAP family members play a role in plant adaptation to adverse environments. This work explores the adaptative responses to phosphorus and/or drought conditions in apple and establishes a foundation for an enhanced comprehension of the evolution of PAP families and the exploration of the genes of interest.

A Phosphate-Starvation Enhanced Purple Acid Phosphatase, GmPAP23 Mediates Intracellular Phosphorus Recycling and Yield in Soybean

Plant internal phosphorus (P) recycling is a complex process, which is vital for improving plant P use efficiency. However, the mechanisms underlying phosphate (Pi) release from internal organic-P form remains to be deciphered in crops. Here, we functionally characterised a Pi-starvation responsive purple acid phosphatase (PAP), GmPAP23 in soybean (Glycine max). GmPAP23 could hydrolyse a series of Pi-containing compounds in vitro, such as trehalose-6-phosphate and glucose-l-phosphate. Moreover, GmPAP23 overexpression led to less P distribution in soybean source organs, including mature leaves and pod shells, but more P distribution in seeds under P sufficient conditions, although no effect was observed for transgenic soybean lines with its suppression. Metabolomic analysis found that a group of P-containing metabolites exhibited differential accumulations in mature leaves between wild type (WT) and GmPAP23 overexpression lines, such as glucose-l-phosphate and trehalose-6-phosphate. Moreover, a MYB transcription factor, GmPHR14 was subsequently found to activate the transcription of GmPAP23 via directly binding to its promoter. Collectively, these findings could highlight that the GmPHR14-GmPAP23 pathway, which controls internal P recycling in soybean, and thus affect yield.

Analysis and profiling of the purple acid phosphatase gene family in wheat (Triticum aestivum L.)

Purple acid phosphatases (PAPs) play a vital role in plant phosphorus nutrition, serving as a crucial family of metallo-phosphoesterase enzymes. This research aimed to identify the PAP genes from the A/B/D genomes of Triticum aestivum to elucidate evolutionary mechanisms of the gene family in plants and provide genomic information for subsequent research on phosphorous-use efficiency in wheat crops. In total, 105 PAP genes (TaPAPs) were identified from the A/B/D genomes by using the Arabidopsis thaliana and Oryza sativa PAP protein sequences as queries for BLASTP against the wheat protein database. The TaPAPs were grouped into six subfamilies, Ia (17), Ib (26), IIa (11), IIb (30), IIIa (12), and IIIb (9), based on their similarities in the structure of genes and the presence of conserved protein motifs. A majority of TaPAPs were derived from tandemly (20) or segmentally (87) duplicated, with the homoeologous chromosomes 5A/B/D harboring the most duplicated PAP genes. Further analysis indicated that TaPAPs were responsible for the modulation of seed, root, and leaf development and hormone synthesis and signaling, as well as plant responses to abiotic stresses, including low temperatures, drought, and anaerobic conditions. Nine TaPAPs (TaPAP9-4A/4B/4D, TaPAP24-6A/6B/6D, and TaPAP28-7A/7B/7D) were constitutively expressed in diverse tissues such as root, shoot, leaf, spike, and seed, while the remaining genes exhibited tissue-specific expression patterns. Concerning the response to phosphate (Pi) deprivation, 57 TaPAPs were highly expressed in roots under Pi stress, including TaPAP31-4A, 4B, and 4D homeologs from the subfamily IIIb. A TaPAP31-4A transgene in A. thaliana promoted plant growth and development while increasing plant resistance to Pi-deficiency stress by enhancing the secretion of phosphatase. These discoveries provide a scientific foundation for comprehending the role of TaPAPs, offering valuable insights for identifying additional candidate genes and fostering the development of new wheat varieties with enhanced tolerance to low phosphorus conditions.

The vacuolar phosphatases purple acid phosphatase 26 and haloacid dehalogenase IIA2.1 hydrolyze 5'-, 3'-, and 2'-nucleotides derived from RNA degradation

The vacuole is an important site for RNA degradation. Autophagy delivers RNA to the vacuole, where the vacuolar T2 RNase ribonuclease 2 (RNS2) plays a major role in RNA catabolism. The presumed products of RNS2 activity are 3'-nucleoside monophosphates (3'-NMPs). Vacuolar phosphatases that carry out 3'-NMP hydrolysis are required to metabolize 3'-NMPs, but the specific players remain unknown. Using a mutant of RNS2 and mutants of the autophagy-related genes 5 and 9 (atg5 and atg9), we confirmed that 3'-NMPs are products of vacuolar RNS2-mediated RNA degradation in Arabidopsis (Arabidopsis thaliana). Moreover, we identified purple acid phosphatase 26 (PAP26) and haloacid dehalogenase IIA2.1 (HIIA2.1) as vacuolar 3'-NMP phosphatases. Based on phylogenetic analysis, we propose systematic nomenclature for HADIIA enzymes, some of which were previously named vegetative storage proteins, which we critically discuss. PAP26 and HIIA2.1 differ in their NMP specificity and activity in vitro. However, hiia2.1 pap26 double mutant plants, but generally not the respective single mutants, accumulate 3'-NMPs in addition to 5'-NMPs and, surprisingly, also 2'-NMPs. These findings suggest that PAP26 and HIIA2.1 have overlapping NMP substrate spectra in vivo. Excess 3'- and 2'-NMPs accumulate in plants exposed to a prolonged night, presumably because carbon limitation enhances autophagy-mediated vacuolar RNA degradation. We conclude that vacuolar RNA catabolism releases 3'-NMPs and 2'-NMPs through RNS2 and other RNases that also generate 5'-NMPs. PAP26 and HIIA2.1 are required to dephosphorylate these NMPs, so that they can enter general nucleotide metabolism outside the vacuole.

Genome-wide association study identified BnaPAP17 genes involved in exogenous ATP utilization and regulating phosphorous content in Brassica napus

KEY MESSAGE

BnaPAP17s associated with root-secreted APases activity were identified by genome-wide association study, and those were induced by Pi-deficiency. BnaPAP17s were involved in improving exogenous organophosphorus utilization as secreted APases. Deficiency of available phosphorus (P) in soil has become an important limiting factor for yield and quality in oilseed rape (Brassica napus). In many soils, organic P (Po) is the main component of the soil P pool. Po must be hydrolyzed to inorganic P (Pi) through acid Phosphatase (APases), and then taken up by plants. However, root-secreted APases (SAP) activity, as a quantitative trait, plays an important role in soil Po utilization; those genetic loci are not clear in B. napus. In this study, we performed a genome-wide association study for SAP activity under Pi-deficiency using a panel of 350 accessions of B. napus and more than 4.5 million polymorphic single nucleotide polymorphisms (SNPs). Thirty-five significant SNPs associated with SAP activity were identified. BnaA01.PAP17 (BnaA01g27810D) was a candidate gene underlying lead SNP (ChrA01_19576615). We experimentally verified that both BnaA01.PAP17 and its three homologous genes had similar expression pattern in response to Pi-deficiency. The dynamic changes in BnaPAP17s expression level were opposite to those of Pi concentration in both roots and leaves, suggesting their potential utility as Pi marker genes in B. napus. Transient expression of BnaPAP17s in tobacco leaves proved that BnaPAP17s were located in the apoplast as secreted APases. The overexpression of BnaPAP17s enhanced SAP activity in response to Pi-deficiency and resulting in increased P content in plants when ATP was supplied as the sole P resource. Taken together, these results suggest that BnaPAP17s contributed to SAP activity, thus having a function in extracellular Po utilization in B. napus.

Enhanced biodegradation of glyphosate by Chlorella sorokiniana engineered with exogenous purple acid phosphatase

Organophosphate pesticides, particularly glyphosate, persist in aquatic environments due to widespread agricultural usage, posing substantial environmental and health risks. This study explores the bioremediation potential of genetically engineered Chlorella sorokiniana, expressing purple acid phosphatase (PAP) from Phaeodactylum tricornutum, for glyphosate biodegradation. The engineered strain (OE line) demonstrated complete glyphosate biodegradation at concentrations below 10 ppm within 4-6 days, surpassing the wild type (WT). Enhanced biodegradation in the OE line was attributed to increased growth and ATP levels due to the release of inorganic phosphate, indicating enhanced metabolic efficiency. Photosynthetic parameters, as well as chlorophyll, and carotenoid contents, were significantly improved, driving higher biomass accumulation. Metabolic shifts toward lipogenesis were observed, supported by the upregulation of triacylglycerol-related genes. Additionally, antioxidant enzyme activities (GPx, SOD, CAT) were elevated in the OE line, mitigating oxidative stress. Importantly, the overexpression of PAP activated and upregulated the level of endogenous CsPAP18, which displayed stable binding with glyphosate and its metabolite aminomethylphosphonic acid, highlighting the synergistic role of PAP and CsPAP18 in glyphosate biodegradation. The OE line effectively treated glyphosate-contaminated real wastewater, confirming the feasibility of engineered strain for environmental remediation. This study provides valuable insights into the potential of engineered microalgae for effective and sustainable wastewater treatment, specifically targeting the removal of organophosphate contaminants in freshwater environments.

Novel phosphatase PvPAP1 from the As-hyperaccumulator Pteris vittata promotes organic P utilization and plant growth: Extracellular exudation and phytate hydrolysis

Organic phosphorus (Po) is a large component of soil P, but it is often unavailable for plant uptake. Purple acid phosphatases (PAP) can hydrolyze a wide range of Po, playing an important role in Po utilization by plants. In this study, we investigated a novel secretary PvPAP1 from the As-hyperaccumulator Pteris vittata, which can effectively utilize exogenous Po, including adenosine triphosphate (ATP) and phytate. Unlike other PAP, PvPAP1 was abundantly-expressed in P. vittata roots, which was upregulated 3.5-folds under P-deprivation than P-sufficient conditions. When expressed in tobacco, its activity in the roots of PvPAP1-Ex lines was ∼8 folds greater than that in wild-type (WT) plants. Besides, PvPAP1 exhibited its secretory ability as evidenced by the sapphire-blue color on the root surface after treating with 5-bromo-4-chloro-3-indolyl phosphate. In a long-term experiment using sand media, PvPAP1-expressing tobacco plants showed 25-30 % greater root biomass than WT plants when using ATP as the sole P source. This is because PvPAP1-expression enhanced its phosphatase activity by 6.5-9.2 folds in transgenic tobacco, thereby increasing the P contents by 39-41 % in its roots under ATP treatment and 9.4-30 % under phytate treatment. The results highlight PvPAP1 as a novel secreted phosphatase crucial for external Po utilization in P. vittata, suggesting that PvPAP1 has the potential to serve as a valuable gene resource for enhancing Po utilization by crop plants.

Exploring the dynamic adaptive responses of Epimedium pubescens to phosphorus deficiency by Integrated transcriptome and miRNA analysis

Phosphorus, a crucial macronutrient essential for plant growth and development. Due to widespread phosphorus deficiency in soils, phosphorus deficiency stress has become one of the major abiotic stresses that plants encounter. Despite the evolution of adaptive mechanisms in plants to address phosphorus deficiency, the specific strategies employed by species such as Epimedium pubescens remain elusive. Therefore, this study observed the changes in the growth, physiological reponses, and active components accumulation in E. pubescensunder phosphorus deficiency treatment, and integrated transcriptome and miRNA analysis, so as to offer comprehensive insights into the adaptive mechanisms employed by E. pubescens in response to phosphorus deficiency across various stages of phosphorus treatment. Remarkably, our findings indicate that phosphorus deficiency induces root growth stimulation in E. pubescens, while concurrently inhibiting the growth of leaves, which are of medicinal value. Surprisingly, this stressful condition results in an augmented accumulation of active components in the leaves. During the early stages (30 days), leaves respond by upregulating genes associated with carbon metabolism, flavonoid biosynthesis, and hormone signaling. This adaptive response facilitates energy production, ROS scavenging, and morphological adjustments to cope with short-term phosphorus deficiency and sustain its growth. As time progresses (90 days), the expression of genes related to phosphorus cycling and recycling in leaves is upregulated, and transcriptional and post-transcriptional regulation (miRNA regulation and protein modification) is enhanced. Simultaneously, plant growth is further suppressed, and it gradually begins to discard and decompose leaves to resist the challenges of long-term phosphorus deficiency stress and sustain survival. In conclusion, our study deeply and comprehensively reveals adaptive strategies utilized by E. pubescens in response to phosphorus deficiency, demonstrating its resilience and thriving potential under stressful conditions. Furthermore, it provides valuable information on potential target genes for the cultivation of E. pubescens genotypes tolerant to low phosphorus.

The Myb73-GDPD2-GA2ox1 transcriptional regulatory module confers phosphate deficiency tolerance in soybean

Phosphorus is indispensable in agricultural production. An increasing food supply requires more efficient use of phosphate due to limited phosphate resources. However, how crops regulate phosphate efficiency remains largely unknown. Here, we identified a major quantitative trait locus, qPE19, that controls 7 low-phosphate (LP)-related traits in soybean (Glycine max) through linkage mapping and genome-wide association studies. We identified the gene responsible for qPE19 as GLYCEROPHOSPHORYL DIESTER PHOSPHODIESTERASE2 (GmGDPD2), and haplotype 5 represents the optimal allele favoring LP tolerance. Overexpression of GmGDPD2 significantly affects hormone signaling and improves root architecture, phosphate efficiency and yield-related traits; conversely, CRISPR/Cas9-edited plants show decreases in these traits. GmMyb73 negatively regulates GmGDPD2 by directly binding to its promoter; thus, GmMyb73 negatively regulates LP tolerance. GmGDPD2 physically interacts with GA 2-oxidase 1 (GmGA2ox1) in the plasma membrane, and overexpressing GmGA2ox1 enhances LP-associated traits, similar to GmGDPD2 overexpression. Analysis of double mutants for GmGDPD2 and GmGA2ox1 demonstrated that GmGDPD2 regulates LP tolerance likely by influencing auxin and gibberellin dose-associated cell division in the root. These results reveal a regulatory module that plays a major role in regulating LP tolerance in soybeans and is expected to be utilized to develop phosphate-efficient varieties to enhance soybean production, particularly in phosphate-deficient soils.

Trade-offs between root-secreted acid phosphatase and root morphology traits, and their contribution to phosphorus acquisition in Brassica napus

Oilseed rape (Brassica napus) is one of the most important oil crops in the world and shows sensitivity to low phosphorus (P) availability. In many soils, organic P (Po) is the main component of the soil P pool. Po must be mineralised to Pi through phosphatases, and then taken up by plants. However, the relationship between root-secreted acid phosphatases (APase) and root morphology traits, two important P-acquisition strategies in response to P deficiency, is unclear among B. napus genotypes. This study aimed to understand their relationship and how they affect P acquisition, which is crucial for the sustainable utilisation of agricultural P resources. This study showed significant genotypic variations in root-secreted APase activity per unit root fresh weight (SAP) and total root-secreted APase activity per plant (total SAP) among 350 B. napus genotypes. Seed yield was positively correlated with total SAP but not significantly correlated with SAP. Six root traits of 18 B. napus genotypes with contrasting root biomass were compared under normal Pi, low Pi and Po. Genotypes with longer total root length (TRL) reduced SAP, but those with shorter TRL increased SAP under P deficiency. Additionally, TRL was important in P-acquisition under three P treatments, and total SAP was also important in P-acquisition under Po treatment. In conclusion, trade-offs existed between the two P-acquisition strategies among B. napus genotypes under P-deficient conditions. Total SAP was an important root trait under Po conditions. These results might help to breed B. napus with greater P-acquisition ability under low P availability conditions.

Multi-omics analysis reveals the roles of purple acid phosphatases in organic phosphorus utilization by the tropical legume Stylosanthes guianensis

Stylo (Stylosanthes guianensis) is a tropical legume known for its exceptional tolerance to low phosphate (Pi), a trait believed to be linked to its high acid phosphatase (APase) activity. Previous studies have observed genotypic variations in APase activity in stylo; however, the gene encoding the crucial APase responsible for this variation remains unidentified. In this study, transcriptomic and proteomic analyses were employed to identify eight Pi starvation-inducible (PSI) APases belonging to the purple APase (PAP) family in the roots of stylo and seven in the leaves. Among these PSI-PAPs, SgPAP7 exhibited a significantly positive correlation in its expression levels with the activities of both internal APase and root-associated APase across 20 stylo genotypes under low-Pi conditions. Furthermore, the recombinant SgPAP7 displayed high catalytic activity toward adenosine 5'-diphosphate (ADP) and phosphoenolpyruvate (PEP) in vitro. Overexpression (OE) of SgPAP7 in Arabidopsis facilitated exogenous organic phosphorus utilization. Moreover, SgPAP7 OE lines showed lower shoot ADP and PEP levels than the wild type, implying that SgPAP7 is involved in the catabolism and recycling of endogenous ADP and PEP, which could be beneficial for plant growth in low-Pi soils. In conclusion, SgPAP7 is a key gene with a major role in stylo adaptation to low-Pi conditions by facilitating the utilization of both exogenous and endogenous organic phosphorus sources. It may also function as a PEP phosphatase involved in a glycolytic bypass pathway that minimizes the need for adenylates and Pi. Thus, SgPAP7 could be a promising target for improving tolerance of crops to low-Pi availability.

SlPHL1 positively modulates acid phosphatase in response to phosphate starvation by directly activating the genes SlPAP10b and SlPAP15 in tomato

Increased acid phosphatase (APase) activity is a prominent feature of tomato (Solanum lycopersicum) responses to inorganic phosphate (Pi) restriction. SlPHL1, a phosphate starvation response (PHR) transcription factor, has been identified as a positive regulator of low Pi (LP)-induced APase activity in tomato. However, the molecular mechanism underlying this regulation remains to be elucidated. Here, SlPHL1 was found to positively regulate the LP-induced expression of five potential purple acid phosphatase (PAP) genes, namely SlPAP7, SlPAP10b, SlPAP12, SlPAP15, and SlPAP17b. Furthermore, we provide evidence that SlPHL1 can stimulate transcription of these five genes by binding directly to the PHR1 binding sequence (P1BS) located on their promoters. The P1BS mutation notably weakened SlPHL1 binding to the promoters of SlPAP7, SlPAP12, and SlPAP17b but almost completely abolished SlPHL1 binding to the promoters of SlPAP10b and SlPAP15. As a result, the transcriptional activation of SlPHL1 on SlPAP10b and SlPAP15 was substantially diminished. In addition, not only did transient overexpression of either SlPAP10b or SlPAP15 in tobacco leaves increase APase activity, but overexpression of SlPAP15 in Arabidopsis and tomato also increased APase activity and promoted plant growth. Subsequently, two SPX proteins, SlSPX1 and SlSPX4, were shown to physically interact with SlPHL1. Moreover, SlSPX1 inhibited the transcriptional activation of SlPHL1 on SlPAP10b and SlPAP15 and negatively regulated the activity of APase. Taken together, these results demonstrate that SlPHL1-mediated LP signaling promotes APase activity by activating the transcription of SlPAP10b and SlPAP15, which may provide valuable insights into the mechanisms of tomato response to Pi-limited stress.

Impaired glycosylation of GmPAP15a, a root-associated purple acid phosphatase, inhibits extracellular phytate-P utilization in soybean

Phosphorus (P) is an essential nutrient, but easily fixed in soils. Therefore, most of soil P exists in the form of inaccessible organic phosphorus (Po), particularly phytate-P. Root-associated purple acid phosphatases (PAPs) are considered to play a crucial role in phosphate (Pi) scavenging in soils. However, evidence for regulating root-associated PAPs in utilization of extracellular phytate-P remain largely unknown in plants at both transcriptional and posttranslational levels. In this study, a Pi-starvation responsive GmPAP15a was identified in soybean (Glycine max). Overexpressing GmPAP15a led to significant increases in root-associated phytase activities, as well as total P content when phytate-P was supplied as the sole P resource in soybean hairy roots. Meanwhile, mass spectrometry (MS) analysis showed GmPAP15a was glycosylated at Asn144 and Asn502 , and its glycan structures of N-linked oligosaccharide chains exhibited microheterogeneity. Moreover, two homologues of AtPHR1, GmPHR9 and GmPHR32 were found to activate GmPAP15a transcription through luciferase activity analysis. Taken together, it is strongly suggested that GmPAP15a plays a vital role in phytate-P utilization in soybean, which might be regulated at both transcriptional and glycosylation modification levels. Our results highlight the GmPHR9/GmPHR32-GmPAP15a signalling pathway might present, and control phytate-P utilization in soybean.

Can nutrients act as signals under abiotic stress?

Plant cells are in constant communication to coordinate development processes and environmental reactions. Under stressful conditions, such communication allows the plant cells to adjust their activities and development. This is due to intercellular signaling events which involve several components. In plant development, cell-to-cell signaling is ensured by mobile signals hormones, hydrogen peroxide (H2O2), nitric oxide (NO), or hydrogen sulfide (H2S), as well as several transcription factors and small RNAs. Mineral nutrients, including macro and microelements, are determinant factors for plant growth and development and are, currently, recognized as potential signal molecules. This review aims to highlight the role of nutrients, particularly calcium, potassium, magnesium, nitrogen, phosphorus, and iron as signaling components with special attention to the mechanism of response against stress conditions.

Phosphorus absorption kinetics and exudation strategies of roots developed by three lupin species to tackle P deficiency

MAIN CONCLUSION

Different lupin species exhibited varied biomass, P allocation, and physiological responses to P-deprivation. White and yellow lupins had higher carboxylate exudation rates, while blue lupin showed the highest phosphatase activity. White lupin (Lupinus albus) can produce specialized root structures, called cluster roots, which are adapted to low-phosphorus (P) soil. Blue lupin (L. angustifolius) and yellow lupin (L. luteus), which are two close relatives of white lupin, do not produce cluster roots. This study characterized plant responses to nutrient limitation by analyzing biomass accumulation and P distribution, absorption kinetics and root exudation in white, blue, and yellow lupins. Plants were grown in hydroponic culture with (64 µM NaH2PO4) or without P for 31 days. Under P limitation, more biomass was allocated to roots to improve P absorption. Furthermore, the relative growth rate of blue lupin showed the strongest inhibition. Under + P conditions, the plant total-P contents of blue lupin and yellow lupin were higher than that of white lupin. To elucidate the responses of lupins via the perspective of absorption kinetics and secretion analysis, blue and yellow lupins were confirmed to have stronger affinity and absorption capacity for orthophosphate after P-deprivation cultivation, whereas white lupin and yellow lupin had greater ability to secrete organic acids. The exudation of blue lupin had higher acid phosphatase activity. This study elucidated that blue lupin was more sensitive to P-scarcity stress and yellow had the greater tolerance of P-deficient condition than either of the other two lupin species. The three lupin species have evolved different adaptation strategies to cope with P deficiency.

Predicting the allergenic risk of Phosphite-NAD+-Oxidoreductase and purple acid phosphatase 17 proteins in genetically modified canola using bioinformatic approaches

Recent advancements in the generation of high-throughput multi-omics data have provided a vast array of candidate genes for the genetic engineering of plants. However, as part of their safety assessment, newly expressed proteins in genetically modified crops must be evaluated for potential cross-reactivity with known allergens. In this study, we developed transgenic canola plants expressing the Arabidopsis thaliana PAP17 gene and a novel selectable marker composed of the ptxD gene from Pseudomonas stutzeri. To evaluate the potential allergenic cross-reactivity of the AtPAP17 and PTXD proteins expressed in transgenic canola, we applied a comprehensive approach utilizing sequence-based, motif-based, and 3D structure-based analyses. Our results demonstrate that the risk of conferring cross-reactivity with known allergens is negligible, indicating that the expression of these proteins in transgenic canola poses a low allergenic risk.

Imazethapyr disrupts plant phosphorus homeostasis and acquisition strategies

The deficiency of essential mineral nutrients caused by xenobiotics often results in plant mortality or an inability to complete its life cycle. Imazethapyr, a widely utilized imidazolinone herbicide, has a long-lasting presence in the soil-plant system and can induce toxicity in non-target plants. However, the effects of imazethapyr on mineral nutrient homeostasis remain poorly comprehended. In this study, Arabidopsis seedlings exposed to concentrations of 4 and 10 μg/L imazethapyr showed noticeable reductions in shoot development and displayed a distinct dark purple color, which is commonly associated with phosphorus (P) deficiency in crops. Additionally, the total P content in both the shoots and roots of Arabidopsis significantly decreased following imazethapyr treatment when compared to the control groups. Through the complementary use of physiological and molecular analyses, we discovered that imazethapyr hinders the abundance and functionality of inorganic phosphorus (Pi) transporters and acid phosphatase. Furthermore, imazethapyr impairs the plant's Pi-deficiency adaptation strategies, such as inhibiting Pi transporter activities and impeding root hair development, which ultimately exacerbate P starvation. These results provide compelling evidence that residues of imazethapyr have the potential to disrupt plant P homeostasis and acquisition strategies. These findings offer valuable insights for risk assessment and highlight the need to reconsider the indiscriminate use of imazethapyr, particularly under specific scenarios such as nutrient deficiency.

Phosphorus Plays Key Roles in Regulating Plants' Physiological Responses to Abiotic Stresses

Phosphorus (P), an essential macronutrient, plays a pivotal role in the growth and development of plants. However, the limited availability of phosphorus in soil presents significant challenges for crop productivity, especially when plants are subjected to abiotic stresses such as drought, salinity and extreme temperatures. Unraveling the intricate mechanisms through which phosphorus participates in the physiological responses of plants to abiotic stresses is essential to ensure the sustainability of agricultural production systems. This review aims to analyze the influence of phosphorus supply on various aspects of plant growth and plant development under hostile environmental conditions, with a special emphasis on stomatal development and operation. Furthermore, we discuss recently discovered genes associated with P-dependent stress regulation and evaluate the feasibility of implementing P-based agricultural practices to mitigate the adverse effects of abiotic stress. Our objective is to provide molecular and physiological insights into the role of P in regulating plants' tolerance to abiotic stresses, underscoring the significance of efficient P use strategies for agricultural sustainability. The potential benefits and limitations of P-based strategies and future research directions are also discussed.

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.

Multi-omics-based identification of purple acid phosphatases and metabolites involved in phosphorus recycling in stylo root exudates

Stylo (Stylosanthes guianensis) is a tropical forage and cover crop that possesses low phosphate (Pi) tolerance traits. However, the mechanisms underlying its tolerance to low-Pi stress, particularly the role of root exudates, remain unclear. This study employed an integrated approach using physiological, biochemical, multi-omics, and gene function analyses to investigate the role of stylo root exudates in response to low-Pi stress. Widely targeted metabolomic analysis revealed that eight organic acids and one amino acid (L-cysteine) were significantly increased in the root exudates of Pi-deficient seedlings, among which tartaric acid and L-cysteine had strong abilities to dissolve insoluble-P. Furthermore, flavonoid-targeted metabolomic analysis identified 18 flavonoids that were significantly increased in root exudates under low-Pi conditions, mainly belonging to the isoflavonoid and flavanone subclasses. Additionally, transcriptomic analysis revealed that 15 genes encoding purple acid phosphatases (PAPs) had upregulated expression in roots under low-Pi conditions. Among them, SgPAP10 was characterized as a root-secreted phosphatase, and overexpression of SgPAP10 enhanced organic-P utilization by transgenic Arabidopsis. Overall, these findings provide detailed information regarding the importance of stylo root exudates in adaptation to low-Pi stress, highlighting the plant's ability to release Pi from organic-P and insoluble-P sources through root-secreted organic acids, amino acids, flavonoids, and PAPs.

Allelic Variation in GmPAP14 Alters Gene Expression to Affect Acid Phosphatase Activity in Soybean

Improvement in acid phosphatase (APase) activity is considered as an important approach to enhance phosphorus (P) utilization in crops. Here, GmPAP14 was significantly induced by low P (LP), and its transcription level in ZH15 (P efficient soybean) was higher than in NMH (P inefficient soybean) under LP conditions. Further analyses demonstrated that there were several variations in gDNA (G-GmPAP14Z and G-GmPAP14N) and the promoters (P-GmPAP14Z and P-GmPAP14N) of GmPAP14, which might bring about differential transcriptional levels of GmPAP14 in ZH15 and NMH. Histochemical staining measurements revealed that a stronger GUS signal was present in transgenic Arabidopsis with P-GmPAP14Z under LP and normal P (NP) conditions compared with the P-GmPAP14N plant. Functional research demonstrated that transgenic Arabidopsis with G-GmPAP14Z had a higher level of GmPAP14 expression than the G-GmPAP14N plant. Meanwhile, higher APase activity was also observed in the G-GmPAP14Z plant, which led to increases in shoot weight and P content. Additionally, validation of variation in 68 soybean accessions showed that varieties with Del36 displayed higher APase activities than the del36 plant. Thus, these results uncovered that allelic variation in GmPAP14 predominantly altered gene expression to influence APase activity, which provided a possible direction for research of this gene in plants.

Biochemical and Molecular Responses Underlying the Contrasting Phosphorus Use Efficiency in Ryegrass Cultivars

Improving plant ability to acquire and efficiently utilize phosphorus (P) is a promising approach for developing sustainable pasture production. This study aimed to identify ryegrass cultivars with contrasting P use efficiency, and to assess their associated biochemical and molecular responses. Nine ryegrass cultivars were hydroponically grown under optimal (0.1 mM) or P-deficient (0.01 mM) conditions, and P uptake, dry biomass, phosphorus acquisition efficiency (PAE) and phosphorus utilization efficiency (PUE) were evaluated. Accordingly, two cultivars with high PAE but low PUE (Ansa and Stellar), and two cultivars with low PAE and high PUE (24Seven and Extreme) were selected to analyze the activity and gene expression of acid phosphatases (APases), as well as the transcript levels of P transporters. Our results showed that ryegrass cultivars with high PAE were mainly influenced by root-related responses, including the expression of genes codifying for the P transporter LpPHT1;4, purple acid phosphatase LpPAP1 and APase activity. Moreover, the traits that contributed greatly to enhanced PUE were the expression of LpPHT1;1/4 and LpPHO1;2, and the APase activity in shoots. These outcomes could be useful to evaluate and develop cultivars with high P-use efficiency, thus contributing to improve the management of P in grassland systems.

Integrated multi-omics reveals the molecular mechanisms underlying efficient phosphorus use under phosphate deficiency in elephant grass (Pennisetum purpureum)

Phosphorus (P) is an essential macronutrient element for plant growth, and deficiency of inorganic phosphate (Pi) limits plant growth and yield. Elephant grass (Pennisetum purpureum) is an important fodder crop cultivated widely in tropical and subtropical areas throughout the world. However, the mechanisms underlying efficient P use in elephant grass under Pi deficiency remain poorly understood. In this study, the physiological and molecular responses of elephant grass leaves and roots to Pi deficiency were investigated. The results showed that dry weight, total P concentration, and P content decreased in Pi-deprived plants, but that acid phosphatase activity and P utilization efficiency (PUE) were higher than in Pi-sufficient plants. Regarding Pi starvation-responsive (PSR) genes, transcriptomics showed that 59 unigenes involved in Pi acquisition and transport (especially 18 purple acid phosphatase and 27 phosphate transporter 1 unigenes) and 51 phospholipase unigenes involved in phospholipids degradation or Pi-free lipids biosynthesis, as well as 47 core unigenes involved in the synthesis of phenylpropanoids and flavonoids, were significantly up-regulated by Pi deprivation in leaves or roots. Furthermore, 43 unigenes related to Pi-independent- or inorganic pyrophosphate (PPi)-dependent bypass reactions were markedly up-regulated in Pi-deficient leaves, especially five UDP-glucose pyrophosphorylase and 15 phosphoenolpyruvate carboxylase unigenes. Consistent with PSR unigene expression changes, metabolomics revealed that Pi deficiency significantly increased metabolites of Pi-free lipids, phenylpropanoids, and flavonoids in leaves and roots, but decreased phospholipid metabolites. This study reveals the mechanisms underlying the responses to Pi starvation in elephant grass leaves and roots, which provides candidate unigenes involved in efficient P use and theoretical references for the development of P-efficient elephant grass varieties.

Genome wide association studies for acid phosphatase activity at varying phosphorous levels in Brassica juncea L

Acid phosphatases (Apases) are an important group of enzymes that hydrolyze soil and plant phosphoesters and anhydrides to release Pi (inorganic phosphate) for plant acquisition. Their activity is strongly correlated to the phosphorus use efficiency (PUE) of plants. Indian mustard (Brassica juncea L. Czern & Coss) is a major oilseed crop that also provides protein for the animal feed industry. It exhibits low PUE. Understanding the genetics of PUE and its component traits, especially Apase activity, will help to reduce Pi fertilizer application in the crop. In the present study, we evaluated 280 genotypes of the diversity fixed foundation set of Indian mustard for Apase activity in the root (RApase) and leaf (LApase) tissues at three- low (5µM), normal (250µM) and high (1mM) Pi levels in a hydroponic system. Substantial effects of genotype and Pi level were observed for Apase activity in both tissues of the evaluated lines. Low Pi stress induced higher mean RApase and LApase activities. However, mean LApase activity was relatively more than mean RApase at all three Pi levels. JM06016, IM70 and Kranti were identified as promising genotypes with higher LApase activity and increased R/S at low Pi. Genome-wide association study revealed 10 and 4 genomic regions associated with RApase and LApase, respectively. Annotation of genomic regions in the vicinity of peak associated SNPs allowed prediction of 15 candidates, including genes encoding different family members of the acid phosphatase such as PAP10 (purple acid phosphatase 10), PAP16, PNP (polynucleotide phosphorylase) and AT5G51260 (HAD superfamily gene, subfamily IIIB acid phosphatase) genes. Our studies provide an understanding of molecular mechanism of the Apase response of B. juncea at varying Pi levels. The identified SNPs and candidate genes will support marker-assisted breeding program for improving PUE in Indian mustard. This will redeem the crop with enhanced productivity under restricted Pi reserves and degrading agro-environments.

Characterization of the PHOSPHATE RESPONSE 2-dependent and -independent Pi-starvation response secretome in rice

Many proteins secreted from plant cells into the surrounding extracellular space help maintain cell structure and regulate stress responses in the external environment. In this study, under Pi-replete and depleted conditions, 652 high-confidence secreted proteins were quantified from wild-type (WT) and PHOSPHATE RESPONSE 2 (OsPHR2)-overexpressing suspension-cultured cells (SCCs). These proteins were functionally grouped as phosphatases, signal transduction proteins, pathogen-related (PR) proteins, cell wall-remodeling proteins, and reactive oxygen species (ROS) metabolism proteins. Although PHOSPHATE RESPONSE (PHR) transcription factors regulate two-thirds of Pi-responsive genes at the transcriptional level, only 30.6% of the Pi-starvation-regulated secreted proteins showed significant changes in OsPHR2-overexpressing SCCs. The OsPHR2-dependent systemic Pi signaling pathway mainly regulates phosphatases and PR proteins, which are involved in the utilization of organophosphate, pathogen resistance, and colonization by rhizosphere microorganisms. The OsPHR2-independent local Pi signaling pathway, on the other hand, largely regulated ROS metabolism proteins, cell wall-remodeling proteins, and signal transduction proteins, which are involved in modifying cell wall structure and root architecture. The functions of differentially expressed secreted proteins between WT and OsPHR2-overexpressing plants under Pi-sufficient and Pi-deficient conditions were further confirmed by analysis of the acid phosphatase activity, ROS content, and cell wall composition.

Root hair-specific transcriptome reveals response to low phosphorus in Cicer arietinum

Root hairs (RH) are a single-cell extension of root epidermal cells. In low phosphorus (LP) availability, RH length and density increase thus expanding the total root surface area for phosphate (Pi) acquisition. However, details on genes involved in RH development and response to LP are missing in an agronomically important leguminous crop, chickpea. To elucidate this response in chickpea, we performed tissue-specific RNA-sequencing and analyzed the transcriptome modulation for RH and root without RH (Root-RH) under LP. Root hair initiation and cellular differentiation genes like RSL TFs and ROPGEFs are upregulated in Root-RH, explaining denser, and ectopic RH in LP. In RH, genes involved in tip growth processes and phytohormonal biosynthesis like cell wall synthesis and loosening (cellulose synthase A catalytic subunit, CaEXPA2, CaGRP2, and CaXTH2), cytoskeleton/vesicle transport, and ethylene biosynthesis are upregulated. Besides RH development, genes involved in LP responses like lipid and/or pectin P remobilization and acid phosphatases are induced in these tissues summarizing a complete molecular response to LP. Further, RH displayed preferential enrichment of processes involved in symbiotic interactions, which provide an additional benefit during LP. In conclusion, RH shows a multi-faceted response that starts with molecular changes for epidermal cell differentiation and RH initiation in Root-RH and later induction of tip growth and various LP responses in elongated RH.

Beneficial Bacterium Azospirillum brasilense Induces Morphological, Physiological and Molecular Adaptation to Phosphorus Deficiency in Arabidopsis

Although most cultivated soils have high levels of total phosphorus (P), the levels of bioavailable inorganic P (Pi) are insufficient. The application of plant-growth-promoting rhizobacteria (PGPR) is an eco-friendly strategy for P utilization; however, PGPR-mediated plant responses that enhance Pi acquisition remain unexplored. Here, we investigated the effect of Azospirillum brasilense on Arabidopsis adaptation to Pi deficiency. Results showed that A. brasilense inoculation alleviated Pi-deficiency-induced growth inhibition and anthocyanin accumulation and increased the total P content in Arabidopsis plants. A comprehensive analysis of root morphology revealed that A. brasilense increased root hair density and length under Pi-limited conditions. We further demonstrated that A. brasilense enhanced the acid phosphatase activity and upregulated the expression of several Pi transporter genes, such as PHOSPHATE1 (PHO1), PHOSPHATE TRANSPORTER 1:(PHT1:1) and PHT1;4. However, A. brasilense did not enhance the growth o total P content in pht1;1, pht1;4 and pht1;1pht1;4 mutants. Moreover, A. brasilense could not increase the P content and PHT1;1 expression in the root hairless mutant rsl4rsl2, because of the occurrence of low-Pi-induced PHT1;1 and PHT1;4 in root hairs. These results indicate that A. brasilense can promote root hair development and enhance acid phosphatase activity and Pi transporter expression levels, consequently improving the Pi absorption capacity and conferring plant tolerance to Pi deficiency.

Monogalactosyl diacylglycerol synthase 3 affects phosphate utilization and acquisition in rice

Galactolipids are essential to compensate for the loss of phospholipids by 'membrane lipid remodelling' in plants under phosphorus (P) deficiency conditions. Monogalactosyl diacylglycerol (MGDG) synthases catalyse the synthesis of MGDG which is further converted into digalactosyl diacylglycerol (DGDG), later replacing phospholipids in the extraplastidial membranes. However, the roles of these enzymes are not well explored in rice. In this study, the rice MGDG synthase 3 gene (OsMGD3) was identified and functionally characterized. We showed that the plant phosphate (Pi) status and the transcription factor PHOSPHATE STARVATION RESPONSE 2 (OsPHR2) are involved in the transcriptional regulation of OsMGD3. CRISPR/Cas9 knockout and overexpression lines of OsMGD3 were generated to explore its potential role in rice adaptation to Pi deficiency. Compared with the wild type, OsMGD3 knockout lines displayed a reduced Pi acquisition and utilization while overexpression lines showed an enhancement of the same. Further, OsMGD3 showed a predominant role in roots, altering lateral root growth. Our comprehensive lipidomic analysis revealed a role of OsMGD3 in membrane lipid remodelling, in addition to a role in regulating diacylglycerol and phosphatidic acid contents that affected the expression of Pi transporters. Our study highlights the role of OsMGD3 in affecting both internal P utilization and P acquisition in rice.

Effect of Fluoride on Cytotoxicity Involved in Mitochondrial Dysfunction: A Review of Mechanism

Fluoride is commonly found in the soil and water environment and may act as chronic poison. A large amount of fluoride deposition causes serious harm to the ecological environment and human health. Mitochondrial dysfunction is a shared feature of fluorosis, and numerous studies reported this phenomenon in different model systems. More and more evidence shows that the functions of mitochondria play an extremely influential role in the organs and tissues after fluorosis. Fluoride invades into cells and mainly damages mitochondria, resulting in decreased activity of mitochondrial related enzymes, weakening of protein expression, damage of respiratory chain, excessive fission, disturbance of fusion, disorder of calcium regulation, resulting in the decrease of intracellular ATP and the accumulation of Reactive oxygen species. At the same time, the decrease of mitochondrial membrane potential leads to the release of Cyt c, causing a series of caspase cascade reactions and resulting in apoptosis. This article mainly reviews the mechanism of cytotoxicity related to mitochondrial dysfunction after fluorosis. A series of mitochondrial dysfunction caused by fluorosis, such as mitochondrial dynamics, mitochondrial Reactive oxygen species, mitochondrial fission, mitochondrial respiratory chain, mitochondrial autophagy apoptosis, mitochondrial fusion disturbance, mitochondrial calcium regulation are emphasized, and the mechanism of the effect of fluoride on cytotoxicity related to mitochondrial dysfunction are further explored.

Transcript Profiling Analysis and ncRNAs' Identification of Male-Sterile Systems of Brassica campestris Reveal New Insights Into the Mechanism Underlying Anther and Pollen Development

Male-sterile mutants are useful materials to study the anther and pollen development. Here, whole transcriptome sequencing was performed for inflorescences in three sterile lines of Chinese cabbage (Brassica campestris L. ssp. chinensis Makino, syn. B. rapa ssp. chinensis), the genic male-sterile line (A line), the Polima cytoplasmic male-sterile (CMS) line (P line), and the Ogura CMS line (O line) along with their maintainer line (B line). In total, 7,136 differentially expressed genes (DEGs), 361 differentially expressed long non-coding RNAs (lncRNAs) (DELs), 56 differentially expressed microRNAs (miRNAs) (DEMs) were selected out. Specific regulatory networks related to anther cell differentiation, meiosis cytokinesis, pollen wall formation, and tapetum development were constructed based on the abortion characteristics of male-sterile lines. Candidate genes and lncRNAs related to cell differentiation were identified in sporocyteless P line, sixteen of which were common to the DEGs in Arabidopsis spl/nzz mutant. Genes and lncRNAs concerning cell plate formation were selected in A line that is defected in meiosis cytokinesis. Also, the orthologs of pollen wall formation and tapetum development genes in Arabidopsis showed distinct expression patterns in the three different sterile lines. Among 361 DELs, 35 were predicted to interact with miRNAs, including 28 targets, 47 endogenous target mimics, and five precursors for miRNAs. Two lncRNAs were further proved to be functional precursors for bra-miR156 and bra-miR5718, respectively. Overexpression of bra-miR5718HG in B. campestris slowed down the growth of pollen tubes, caused shorter pollen tubes, and ultimately affected the seed set. Our study provides new insights into molecular regulation especially the ncRNA interaction during pollen development in Brassica crops.