Next Article in Journal
Phytochemical Profiling, Biological Activities, and In Silico Molecular Docking Studies of Causonis trifolia (L.) Mabb. & J.Wen Shoot
Previous Article in Journal
Meta-Analysis of the Effect of Nitric Oxide Application on Heavy Metal Stress Tolerance in Plants
Previous Article in Special Issue
Ecophysiology of Endangered Plant Species Saussurea esthonica: Effect of Mineral Nutrient Availability and Soil Moisture
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 12
Issue 7
10.3390/plants12071496
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Genome-Wide Identification of SPX Family Genes and Functional Characterization of PeSPX6 and PeSPX-MFS2 in Response to Low Phosphorus in Phyllostachys edulis

by
Jiali Luo
,
Zhihui Liu
,
Jiawen Yan
,
Wenhui Shi
* and
Yeqing Ying
*
State Key Laboratory of Subtropical Silviculture, Zhejiang A&F University, Hangzhou 311300, China
*
Authors to whom correspondence should be addressed.
Plants 2023, 12(7), 1496; https://doi.org/10.3390/plants12071496
Submission received: 19 February 2023 / Revised: 10 March 2023 / Accepted: 14 March 2023 / Published: 29 March 2023
(This article belongs to the Special Issue Impact of Abiotic Stresses on Plant-Soil Nutrient Status, Metabolism, Growth, and Ecosystem Functions)
Browse Figures
Review Reports Versions Notes

Abstract

:
Moso bamboo (Phyllostachys edulis) is the most widely distributed bamboo species in the subtropical regions of China. Due to the fast-growing characteristics of P. edulis, its growth requires high nutrients, including phosphorus. Previous studies have shown that SPX proteins play key roles in phosphorus signaling and homeostasis. However, the systematic identification, molecular characterization, and functional characterization of the SPX gene family have rarely been reported in P. edulis. In this study, 23 SPXs were identified and phylogenetic analysis showed that they were classified into three groups and distributed on 13 chromosomes. The analysis of conserved domains indicated that there was a high similarity between PeSPXs among SPX proteins in other species. RNA sequencing and qRT-PCR analysis indicated that PeSPX6 and PeSPX-MFS2, which were highly expressed in roots, were clearly upregulated under low phosphorus. Co-expression network analysis and a dual luciferase experiment in tobacco showed that PeWRKY6 positively regulated the PeSPX6 expression, while PeCIGR1-2, PeMYB20, PeWRKY6, and PeWRKY53 positively regulated the PeSPX-MFS2 expression. Overall, these results provide a basis for the identification of SPX genes in P. edulis and further exploration of their functions in mediating low phosphorus responses.

1. Introduction

Moso bamboo (Phyllostachys edulis) is the most widely distributed bamboo species in the subtropical regions of China. The planting area is 4,677,800 hectares, accounting for 70% of the total bamboo forest area in China [1,2,3]. P. edulis not only has the highest ecological and economic value, but also has the highest cultural value of all bamboos. Due to the fast-growing characteristics of P. edulis, its growth requires a high amount of nutrients [4]. Thus, an adequate supply of soil nutrients is necessary to ensure the normal growth of P. edulis.
As an essential mineral element during plant growth and development, phosphorus (P) plays key roles in energy transmission, signal transduction, photosynthesis, and respiration [5,6,7]. Plants prefer to obtain inorganic soluble phosphate in the soil through the root system, and cannot directly assimilate the other forms of P, such as organic phosphate (Pi) [8]. However, the content of inorganic phosphorus in soil is very low, and it does not diffuse easily in the soil. As a result, plants have evolved mechanisms to adapt to Pi-deficient conditions [9]. The phosphorus starvation signaling pathway in plants has been well documented [10].
Among them, SPX proteins play very important roles in P signaling and homeostasis [11,12]. The SPX domain is a conservative domain named after SYG1 (suppressor of yeast gpa1), Pho81 (CDK inhibitor in yeast PHO pathway), and XPR1 (xenotropic and polytropic retrovirus receptor). Based on the presence of additional domains, SPX proteins can be further classified into four subfamilies: SPX proteins, which contain only an SPX domain; SPX-EXS proteins, which have an SPX domain and an EXS (ERD1, XPR1, and SYG1) domain; SPX-MFS proteins, which have an SPX domain and a major facility superfamily (MFS) domain; and SPX-RING proteins, which contain an SPX domain and a RING-type zinc finger domain [13,14,15]. SPX proteins have been identified in many plant species, including Arabidopsis, rice, and maize [16]. For example, there are 20 SPX domain-containing proteins in Arabidopsis. AtSPX1 participates in the transcription regulation of Pi-responsive genes, while AtSPX3 functions in the potential negative feedback regulation of the P signal network [17]. 33 SPX gene family members have been found in maize. Xiao et al. [18] found that ZmSPX4.1 and ZmSPX4.2 showed strong responses to low Pi stress and exhibited remarkably different expression patterns in low-Pi sensitive and insensitive cultivars of maize [18]. Currently, the role of SPX genes in response to low phosphorus stress in P. edulis is still poorly understood. Thus, the systematic identification, molecular characterization, and functional characterization of the SPX gene family in response to low Pi in P. edulis are of great importance to improve the absorption and utilization of Pi in P. edulis.
In this study, the whole genome of SPX members in P. edulis was first identified and characterized. The chromosomal mapping, phylogenetic relationships, and conserved motifs of SPX genes were analyzed. Then, the expression levels of SPX members exposed to low Pi were analyzed using RNA sequencing (RNA-seq) and quantitative real-time PCR (qRT-PCR). Furthermore, several transcriptional factors were identified that modulated the expression of PeSPX6 and PeSPX-MFS2 in P. edulis.

2. Result

2.1. Identification of SPX-Domain-Containing Proteins in P. edulis

To identify P. edulis SPXs, the conserved PF03105 domain was used as the probe to execute a genome-wide search of candidate genes with the HMMER tool, and a total of 30 putative protein hits were obtained. To verify the HMMER search results, domain analysis was further performed for the 30 putative proteins with CDD and SMART. Among the 30 putative proteins, 23 harboring the conserved SPX domains were identified using both CDD and SMART. These proteins were designated PeSPX1 through PeSPX1-EXS6 and were considered authentic SPX candidates in the P. edulis genome. The length of these protein-coding regions ranged from 738 bp to 2571 bp. These proteins consisted of 245–856 amino acids, and the corresponding molecular weights ranged from 27.97 to 97.55 kDa. The isoelectric point (pI) value of these PeSPX proteins ranged from 5.04 to 9.4 (Table 1).
The chromosomal distribution of the P. edulis SPX-domain-containing proteins genes is illustrated in Figure 1. Furthermore, 23 PeSPX genes were distributed on 13 P. edulis chromosomes. Chromosomes 4, 5, 10, 14, 20, 23, and 24 each contained one gene; chr 3, 15, and 21 each had two genes; chr 8 and 17 each had three genes; and chr 6 had four genes (Figure 1).

2.2. Phylogenetic Analysis of the SPX-Domain-Containing Protein Genes in P. edulis and Other Species

To evaluate the evolutionary relationships of SPX-domain-containing protein genes in P. edulis, this work analyzed the sequence features in three different species, including A. thaliana, O. sativa, and P. edulis, and a total of 56 SPX-domain-containing protein genes were used to construct a phylogenetic tree with the neighbor-joining (NJ) method using MEGA (version 7) (Figure 2). The phylogenetic tree indicated that PeSPXs could be divided into three subfamilies (the SPX, SPX-MFS, and SPX-EXS subfamilies). In addition, the numbers of SPX subfamily proteins in the three species were highly asymmetrical. For example, 11 PeSPXs, five OsSPXs, and four AtSPXs were classified in the SPX subfamilies, and six PeSPXs, 11 AtSPXs, and three OsSPXs were included in the SPX-EXS subfamilies. Only one OsSPX and two AtSPXs were classified in the SPX-RING subfamilies. These results indicate that the SPX gene family is highly conserved and diverse in different plants.

2.3. Structure Analysis of the SPX-Domain-Containing Genes and Proteins in P. edulis

All SPX-domain-containing proteins contain the SPX domain in the C-terminal portion (Figure 3). The three subfamilies of SPX-domain-containing proteins (SPX, SPX-EXS, and SPX-MFS) were found in P. edulis. The SPX subfamily possessed 11 members, the SPX-EXS subfamily contained six members, and the SPX-MFS subfamily had six members.
Structural features were then characterized for PeSPX genes, and a large divergence in exon number was observed: three exons were detected for PeSPX subfamily, ten for PeSPX-MFS subfamily, 11 for PeSPX-EXS1, 13 for PeSPX-EXS4, 14 for PeSPX-EXS2, PeSPX-EXS3, PeSPX-EXS5, 15 for PeSPX-EXS6 (Figure 3A). PeSPX-MFS3 had the greatest length (approximately 31 kb), with only 10 exons and nine introns.

2.4. Expression Analysis of P. edulis SPX-Domain-Containing Protein Genes

To explore the effects of low phosphorus on the growth and development of P. edulis, P. edulis plants were treated with 1/2 Kimura nutrient solution in different Pi levels. Under low phosphorus conditions, the root length was significantly lower than that under the normal Pi supply, but the number of lateral roots was significantly higher than that under the normal Pi supply (Figure 4A,B). Furthermore, the phosphorus content in the shoots and roots of P. edulis was also measured. As shown in Figure 4C, the phosphorus content in both shoots and roots of P. edulis treated with low phosphorus was much lower compared to those in shoots and roots of P. edulis with a normal phosphorus supply. To further explore the molecular mechanism of P. edulis in response to low phosphorus, RNA-seq was performed (Table S1). The expression data of 23 PeSPXs were clustered and displayed in a heat map (Figure 4). The expression levels of five genes were upregulated, while 18 genes were downregulated under low phosphorus stress.
Next, two PeSPXs (PeSPX6 and PeSPX-MFS2) with higher expression levels were selected for qRT-PCR analysis. Consistent with the RNA-seq result, the expression levels of PeSPX6 and PeSPX-MFS2 were clearly upregulated under low phosphorus (Figure 5A). Furthermore, the expression patterns of PeSPX6 and PeSPX-MFS2 were also analyzed. As shown in Figure 5B, PeSPX6 and PeSPX-MFS2 were more highly expressed in the roots, while the expression was lower in the leaves and stems of P. edulis.

2.5. Validation of SPX Gene Network Regulation

Because the expression levels of PeSPX6 and PeSPX-MFS2 were increased by low phosphorus, this study then investigated what transcriptional regulation of PeSPX6 and PeSPX-MFS2 might exist in response to low phosphorus. Pearson’s correlations were performed between the expression levels of PeSPX6 and PeSPX-MFS2 and transcription factors induced by low phosphorus (Figure 6, Table S3). Eight transcription factors with high expression levels under low phosphorus stress were selected to explore their function in regulating the expression levels of PeSPX6 and PeSPX-MFS2 using a dual-luciferase analysis (Table S4). The promoters of PeSPX6 (2000 bp) and PeSPX-MFS2 (2000 bp) were separately cloned and fused to the firefly luciferase protein (Fluc) at the N-terminus, which also had a Renilla luciferase (Rluc). The Fluc/Rluc ratio represents the ability of transcription factors to transcriptionally activate the PeSPX6 or PeSPX-MFS2 promoter. As shown in Figure 6, PeWRKY6 enhanced the activity of luciferase driven by the PeSPX6 promoter, while PeCIGR1-2, PeMYB20, PeWRKY6, and PeWRKY53 enhanced the activity of luciferase driven by the PeSPX-MFS2 promoter.

3. Discussion

SPX family genes are widely found in eukaryotes, including plants, fungi, and multicellular animals [19,20]. At present, the SPX family genes have been extensively studied in Arabidopsis [17] and rice [21,22]. In rice, there are six SPX subfamily genes. Studies on OsSPX1 and OsSPX genes in rice have found that these two genes can participate in the negative regulation of PHR2 under the condition of rich phosphorus [21]. The SPX-EXS family, involved in phosphorus transport and signaling from root to terminal bud in rice and Arabidopsis thaliana [23,24,25]. Members of the SPX family, including 20 AtSPXs, 69 BnaSPXs, 9 GmSPXs, 14 OsSPXs, and 46 TaSPXs, have been reported and characterized by bioinformatics analysis [17,19,26,27]. However, the functional analysis of the SPX family in P. edulis has not been reported. In the present study, systematic identification, molecular characterization, and functional characterization of the SPX gene family were performed.
Previous studies have shown that SPX family genes play an important role in the sensing, signal transduction, and transport of inorganic phosphate (Pi) in eukaryotes, such as PHO1 (SPX-EXS) [28], PHT5 (SPX-MFS) [29] and AtNLA (SPX-RING) [30]. AtSPX1 was expressed most rapidly in Pi starvation induction, which indicated that AtSPX1 has a potential transcriptional regulation effect on Pi starvation [19]. In rice, eight SPX genes (OsSPX1, OsPHO2) were significantly increased under exposure to low phosphorus [19,31]. The expression of TaSPX2 was significantly induced in wheat under low P stress [32]. It was found that SPX2 was induced to express in tea plants under low phosphorus stress [33]. Bn SPX2 was continuously induced to express in various tissue parts of Brassica napus under low P stress [27]. In the present study, only five SPX genes were clearly upregulated under low phosphorus. Among them, two genes (PeSPX6 and PeSPX-MFS2) showed higher expression levels compared to the others in the roots. This strongly implies that PeSPX6 and PeSPX-MFS2 might regulate the low phosphorus response in P. edulis. In addition, the tissue-specific expression indicated that PeSPX6 and PeSPX-MFS2 were more highly expressed in roots, which further supported the hypothesis that PeSPX6 and PeSPX-MFS2 were involved in phosphorus nutrient uptake in the roots of P. edulis. In future research, the biological function of these genes in P. edulis needs to be explored.
Increasing evidence shows that proteins containing the SPX domain are key players in a series of processes that control the dynamic balance of phosphorus in plants. In Arabidopsis, the physical interaction between AtSPXs and AtPHRs under sufficient Pi conditions prevents AtPHRs from binding to the promoter of AtPSI, thus inhibiting the PHR transcriptional activity [34,35,36]. AtSPX3 is involved in responses to low Pi stress [23]. In rice, phosphorus starvation induces the accumulation of OsSPX3 to restore the phosphorus balance [37]. PeSPX6 has the highest similarity to OsSPX3 and AtSPX3, suggesting that PeSPX6 might be involved in the response to low phosphorus and play a role in the Pi starvation signal transduction pathway of P. edulis. As a volatile phosphate efflux transporter, OsSPX-MFS3 is involved in maintaining phosphate homeostasis in rice. Phylogenetic analysis showed that PeSPX-MFS2 and OsSPX-MFS3 had high homology. Thus, this study speculates that PeSPX-MFS2 might be a low-affinity Pi transporter that mediates Pi efflux from the vacuole into the cytosol. In general, the plant SPX family can be divided into four subclasses based on the presence of additional protein domains, namely SPX, SPX⁃EXS, SPX⁃MFS, and SPX⁃RING [13]. Among them, SPX-RING is the smallest subfamily in the SPX family. For example, there is only one member of the SPX-RING subfamily in rice and tow in Arabidopsis. Surprisingly, this subfamily was not found in the P. edulis genome. This demonstrates the complexity of the SPX family in different plant species.
Wang et al. [35] showed that OsSPX1 and OsSPX2 inhibited phosphate starvation in rice by interacting with PHR2 in a phosphate-dependent manner. AtSPX4 regulates the PHR1-dependent and independent regulation of stem phosphorus status in Arabidopsis [38]. In this study, several transcription factors, namely PeMYB20, PeMYB30, PeWRKY6, PeWYKY53, PeERF110, PeNAC030, PeCIGR1-1, and PeCIGR1-2, were upregulated up by low phosphorus. Furthermore, the dual luciferase experiment in tobacco showed that PeWRKY6 positively regulated the expression of PeSPX6, while PeCIGR1-2, PeMYB20, PeWRKY6, and PeWRKY53 positively regulated the expression of PeSPX-MFS2. These results imply that these transcription factors might also play a vital role in the adaptation of P. edulis to low-phosphorus environments.
In this study, phylogenetic analysis showed that 23 SPXs were classified into three groups and distributed on 13 chromosomes. The analysis of conserved domains indicated that there was a high similarity between PeSPXs among SPX proteins in other species. PeSPX6 and PeSPX-MFS2, which were highly expressed in roots, were obviously upregulated under low phosphorus. PeWRKY6 positively regulated the expression of PeSPX6, while PeCIGR1-2, PeMYB20, PeWRKY6, and PeWRKY53 positively regulated the expression of PeSPX-MFS2. Based on this, it is speculated that these genes play different roles in various biological processes, laying a theoretical foundation for elucidating the functions of the SPX family genes in P. edulis.

4. Materials and Methods

4.1. Identification and Classification of SPX Genes

The genomic information of P. edulis was downloaded from the Gigascience database (accession numbers: PRJEB2955 and PRJEB2956). Then, using “hmmsearch” with an expected value (e-value) threshold of 0.5 × 10−3, the genomic protein sequences of P. edulis were searched for the Hidden Markov Model (HMM) profiles of the SPX domain (PF03105). Using the NCBI Conserved Domain Database (NCBI-CDD) (http://www.ncbi.nlm.nih.gov/cdd (accessed on 12 November 2022)) and the SMART database (http://smart.embl-heidelberg.de (accessed on 12 November 2022)) the outcomes were further validated (e-value = 1 × 10−2). The SPX protein sequences of Arabidopsis thaliana and Oryza sativa were searched using Phytozome 13 (https://phytozome-next.jgi.doe.gov/ (accessed on 12 November 2022)).

4.2. Phylogenetic Analysis of SPX Proteins

Multiple amino acid sequence alignments were performed using MAFFT with the default parameters utilizing A. thaliana, O. sativa, and P. edulis sequences to examine the evolutionary relationships and the full-length protein sequences of P. edulis SPXs. Then, using the auto mode in MEGA (version 7), a phylogenetic tree of amino acid sequences based on the neighbor-joining (NJ) was created.

4.3. Gene Structure and Chromosomal Location

Conserved motif structures were examined using the MEME online software (meme.nbcr.net/meme (accessed on 12 November 2022)) and the motif function was explored using NCBI-CDD. To display the locations of PeSPXs on chromosomes, TBtools was employed. The genomic sequence information applied to the analysis of both gene structure and chromosomal location was derived from NCBI [39].

4.4. Transcriptome Sequencing

High-quality RNA was used to create the cDNA library and the 150-bp paired-end reads that resulted were sequenced using the Illumina NovaSeq 6000 platform. After filtering, the cleaning data were compared to the reference genome of P. edulis through HISAT2, and the gene expression levels were counted as Fragments Per Kilobase per Million (FPKM) using StringTie. Based on p-values < 0.05 and fold changes ≥ 2, the R package Deseq2 identified differentially expressed genes (DEGs) in each sample group.

4.5. Determination of Phosphorus Content

The content of phosphorus in each organ of P. edulis was determined using the Mo-Sb colorimetric method [40]. The extracted solution was mixed with 2,4-dinitrophenol, NaOH, and 0.5 mol/L H2SO4 until the yellow was arised. Then, the mixed solution was reacted with Mo-Sb-Vc chromogenic agent at 25 °C for 30 min. The absorbance at 700 nm was measured using an ultraviolet-visible spectrophotometer (V 5600, Shang Hai, China), then the phosphorus content was calculated.

4.6. Real-Time RT-PCR Analysis

Utilizing the cDNA reverse transcription kit, whole pure RNA was transformed into cDNA (PrimeScriptTM RT Master Mix, Takara, Kusatsu, Japan). Then, qRT-PCR was carried out using the CFX96 TouchTM Thermal Cycler and the ChamQ SYBR qPCR Master Mix kit from Vazyme (Vazyme, Nanjing, China) (Bio-Rad, Hercules, CA, USA). The 2−ΔΔCT method was utilized to normalize for transcript levels using the housekeeping gene PeTIP41 (Table S2) [41].

4.7. Luciferase Assays

The full length of transcription factors were separately cloned into pGreenII 62-SK vector. The promoter of PeSPX6 and PeSPX-MFS2 (2000 bp) was cloned into pGreenII 0800 LUC vector (Table S2). Agrobacterium tumefaciens GV3101 was used to convert the aforementioned vectors into four-week-old Nicotiana benthamiana leaves [42]. After infiltration for 4 days, the activity of LUC and REN were assessed in accordance with the instructions provided in the Dual Luciferase Reporter Gene Assay Kit (Beyotime, Nantong, China).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12071496/s1, Table S1. Summary of the RNA-Seq data. Table S2. Primers used in this study. Table S3. Transcription factor correlation analysis. Table S4. Transcription factors gene ID.

Author Contributions

Conceptualization, J.Y.; Data curation, Z.L.; Funding acquisition, W.S. and Y.Y.; Methodology, J.L., Z.L. and J.Y.; Project administration, W.S. and Y.Y.; Supervision, W.S. and Y.Y.; Writing—original draft, J.L.; Writing—review & editing, W.S. and Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

Supported by the National Natural Science Foundation of China (NSFC): “Mechanism of Root Activation and Uptake of Mineral Phosphorus mediated by Phosphorus solubilizing Microorganisms” (32171879), and the National Natural Science Foundation of China (31770645): Effects of Nitrogen Deposition on Response of Phyllostachys pubesculus to Low Phosphorus Stress.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or Supplementary Material.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

SPX:SYG1: PHO81, and XPR1
SPX-EXS:SPX-ERD1/XPR1/SYG1
SPX-MFS:SPX-Major Facilitator Superfamily
SPX-RING:SPX-Really Interesting New Gene
qRT-PCR:quantitative real-time PCR
NJ:Neighbor-joining
aa:Amino acids; in length
MWs:Molecular weights
Pi:Isoelectric points

References

  1. Song, X.; Zhou, G.; Jiang, H.; Yu, S.; Fu, J.; Li, W.; Wang, W.; Ma, Z.; Peng, C. Carbon sequestration by Chinese bamboo forests, and their ecological benefits: Assessment of potential, problems, and future challenges. Environ. Rev. 2011, 19, 418–428. [Google Scholar] [CrossRef] [Green Version]
  2. SFAPRC (State Forestry Administration, P.R. China). Statistics of Forest Resources in China (2009–2013); China Forestry Publishing House: Beijing, China, 2014. [Google Scholar]
  3. Peng, Z.H.; Zhang, C.; Zhang, Y.; Hu, T.; Mu, S.H.; Li, X.; Gao, J. Transcriptome sequencing and analysis of the fast growing shoots of moso bamboo (Phyllostachys edulis). PLoS ONE 2013, 8, e78944. [Google Scholar] [CrossRef]
  4. Tang, Y.L. Study on Phosphorus Utilization Mechanism of Phyllostachys Pubesculi in Response to Different Phosphorus Environments under Simulated Nitrogen Deposition; Zhejiang Agriculture and Forestry University: Hangzhou, China, 2020. [Google Scholar]
  5. Feder, D.; McGeary, R.P.; Mitić, N.; Lonhienne, T.; Furtado, A.; Schulz, B.L.; Henry, R.J.; Schmidt, S.; Guddat, L.W.; Schenk, G. Structural elements that modulate the substrate specificity of plant purple acid phosphatases: Avenues for improved phosphorus acquisition in crops. Plant Sci. 2020, 294, 110445. [Google Scholar] [CrossRef]
  6. Rubio, V. A conserved MYB transcription factor involved in phosphate starvation signaling both in vascular plants and in unicellular algae. Genes Dev. 2001, 15, 2122–2133. [Google Scholar] [CrossRef] [Green Version]
  7. Abel, S.; Ticconi, C.A.; Delatorre, C.A. Phosphate sensing in higher plants. Physiol. Plant. 2002, 115, 1–8. [Google Scholar] [CrossRef]
  8. Robinson, W.D.; Carson, I.; Ying, S.; Ellis, K.; Plaxton, W.C. Eliminating the purple acid phosphatase At PAP 26 in A. thaliana delays leaf senescence and impairs phosphorus remobilization. New Phytol. 2012, 196, 1024–1029. [Google Scholar] [CrossRef]
  9. Vance, C.P.; Uhde-Stone, C.; Allan, D.L. Phosphorus acquisition and use: Critical adaptations by plants for securing a nonrenewable resource. New Phytol. 2003, 157, 423–447. [Google Scholar] [CrossRef] [Green Version]
  10. Ham, B.K.; Chen, J.; Yan, Y.; Lucas, W.J. Insights into plant phosphate sensing and signaling. Curr. Opin. Biotechnol. 2018, 49, 1–9. [Google Scholar] [CrossRef]
  11. Chang, M.X.; Gu, M.; Xia, Y.W.; Dai, X.L.; Dai, C.R.; Zhang, J.; Wang, S.C.; Qu, H.Y.; Yamaji, N.; Feng, M.J.; et al. OsPHT1;3 Mediates Uptake, Translocation, and Remobilization of Phosphate under Extremely Low Phosphate Regimes. Plant Physiol. 2019, 179, 656–670. [Google Scholar] [CrossRef] [Green Version]
  12. Yang, J.; Zhao, X.; Chen, Y.; Li, G.; Li, X.; Xia, M.; Sun, Z.; Chen, Y.; Li, Y.; Yao, L.; et al. Identification, Structural, and Expression Analyses of SPX Genes in Giant Duckweed (Spirodela polyrhiza) Reveals Its Role in Response to Low Phosphorus and Nitrogen Stresses. Cells 2022, 11, 1167. [Google Scholar] [CrossRef]
  13. Chiou, T.J.; Lin, S.I. Signaling network in sensing phosphate availability in plants. Annu. Rev. Plant Biol. 2011, 62, 185–206. [Google Scholar] [CrossRef] [Green Version]
  14. Yang, J.; Wang, L.; Mao, C.; Lin, H. Characterization of the rice NLA family reveals a key role for OsNLA1 in phosphate homeostasis. Rice 2017, 10, 52. [Google Scholar] [CrossRef] [Green Version]
  15. Yue, W.H.; Ying, Y.H.; Wang, C.; Zhao, Y.; Dong, C.G.; Whelan, J.; Shou, H.X. OsNLA1, a RING-type ubiquitin ligase, maintains phosphate homeostasis in Oryza sativa via degradation of phosphate transporters. Plant J. 2017, 90, 1040–1051. [Google Scholar] [CrossRef] [Green Version]
  16. Liu, N.; Shang, W.Y.; Li, C.; Jia, L.H.; Wang, X.; Xing, G.Z.; Zheng, W.M. Evolution of the SPX gene family in plants and its role in the response mechanism to phosphorus stress. Open Biol. 2018, 8, 170231. [Google Scholar] [CrossRef] [Green Version]
  17. Duan, K.; Yi, K.K.; Dang, L.; Huang, H.J.; Wu, W.; Wu, P. Characterization of a sub-family of Arabidopsis genes with the SPX domain reveals their diverse functions in plant tolerance to phosphorus starvation. Plant J. 2008, 54, 965–975. [Google Scholar] [CrossRef]
  18. Xiao, J.B.; Xie, X.M.; Li, C.; Xing, G.Z.; Cheng, K.; Li, H.; Liu, N.; Tan, J.F.; Zheng, W.M. Identification of SPX family genes in the maize genome and their expression under different phosphate regimes. Plant Physiol. Biochem. 2021, 168, 211–220. [Google Scholar] [CrossRef]
  19. Secco, D.; Wang, C.; Arpat, B.A.; Wang, Z.Y.; Poirier, Y.; Tyerman, S.D.; Wu, P.; Shou, H.X.; Whelan, J. The emerging importance of the SPX domain-containing proteins in phosphate homeostasis. New Phytol. 2011, 193, 842–851. [Google Scholar] [CrossRef] [Green Version]
  20. Wild, R.; Gerasimaite, R.; Jung, J.Y.; Truffault, V.; Pavlovic, I.; Schmidt, A.; Saiardi, A.; Jessen, H.J.; Poirier, Y.; Hothorn, M.; et al. Control of eukaryotic phosphate homeostasis by inositol polyphosphate sensor domains. Science 2016, 352, 986–990. [Google Scholar] [CrossRef] [Green Version]
  21. Wang, C.; Ying, S.; Huang, H.; Li, K.; Wu, P.; Shou, H. Involvement of OsSPX1 in phosphate homeostasis in rice. Plant J. 2009, 57, 895–904. [Google Scholar] [CrossRef]
  22. Wang, Z.; Hu, H.; Huang, H.; Duan, K.; Wu, Z.; Wu, P. Regulation of OsSPX1 and OsSPX3 on expression of OsSPX domain genes and Pi-starvation signaling in rice. J. Integr. Plant Biol. 2009, 51, 663–674. [Google Scholar] [CrossRef]
  23. Hamburger, D.; Rezzonico, E.; MacDonald-Comber Petétot, J.; Somerville, C.; Poirier, Y. Identification and characterization of the Arabidopsis PHO1 gene involved in phosphate loading to the xylem. Plant Cell. 2002, 14, 889–902. [Google Scholar] [CrossRef] [Green Version]
  24. Stefanovic, A.; Ribot, C.; Rouached, H.; Wang, Y.; Chong, J.; Belbahri, L.; Delessert, S.; Poirier, Y. Members of the PHO1 gene family show limited functional redundancy in phosphate transfer to the shoot, and are regulated by phosphate deficiency via distinct pathways. Plant J. 2007, 50, 982–994. [Google Scholar] [CrossRef]
  25. Secco, D.; Baumann, A.; Poirier, Y. Characterization of the rice PHO1 gene family reveals a key role for OsPHO1;2 in phosphate homeostasis and the evolution of a distinct clade in dicotyledons. Plant Physiol. 2010, 152, 1693–1704. [Google Scholar] [CrossRef] [Green Version]
  26. Yao, Z.; Tian, J.; Liao, H. Comparative characterization of GmSPX members reveals that GmSPX3 is involved in phosphate homeostasis in soybean. Ann. Bot. 2014, 114, 477–488. [Google Scholar] [CrossRef] [Green Version]
  27. Du, H.; Yang, C.; Ding, G.; Shi, L.; Xu, F. Genome-Wide Identification and Characterization of SPX Domain-Containing Members and Their Responses to Phosphate Deficiency in Brassica napus. Front. Plant Sci. 2017, 8, 35. [Google Scholar] [CrossRef] [Green Version]
  28. Wege, S.; Khan, G.A.; Jung, J.Y.; Vogiatzaki, E.; Pradervand, S.; Aller, I.; Meyer, A.J.; Poirier, Y. The EXS Domain of PHO1 Participates in the Response of Shoots to Phosphate Deficiency via a Root-to-Shoot Signal. Plant Physiol. 2016, 170, 385–400. [Google Scholar] [CrossRef] [Green Version]
  29. Liu, J.; Yang, L.; Luan, M.; Wang, Y.; Zhang, C.; Zhang, B.; Shi, J.; Zhao, F.G.; Lan, W.; Luan, S. A vacuolar phosphate transporter essential for phosphate homeostasis in Arabidopsis. Proc. Natl. Acad. Sci. USA 2015, 112, E6571–E6578. [Google Scholar] [CrossRef] [Green Version]
  30. Lin, W.Y.; Huang, T.K.; Chiou, T.J. Nitrogen limitation adaptation, a target of microRNA827, mediates degradation of plasma membrane-localized phosphate transporters to maintain phosphate homeostasis in Arabidopsis. Plant Cell. 2013, 25, 4061–4074. [Google Scholar] [CrossRef] [Green Version]
  31. Lv, Q.; Zhong, Y.; Wang, Y.; Wang, Z.; Zhang, L.; Shi, J.; Wu, Z.; Liu, Y.; Mao, C.; Yi, K.; et al. SPX4 Negatively Regulates Phosphate Signaling and Homeostasis through Its Interaction with PHR2 in Rice. Plant Cell. 2014, 26, 1586–1597. [Google Scholar] [CrossRef] [Green Version]
  32. Kumar, A.; Sharma, M.; Gahlaut, V.; Nagaraju, M.; Chaudhary, S.; Kumar, A.; Tyagi, P.; Gajula, M.N.V.P.; Singh, K.P. Genome-wide identification, characterization, and expression profiling of SPX gene family in wheat. Int. J. Biol. Macromol. 2019, 140, 17–32. [Google Scholar] [CrossRef]
  33. Kc, S.; Liu, M.; Zhang, Q.; Fan, K.; Shi, Y.; Ruan, J. Metabolic Changes of Amino Acids and Flavonoids in Tea Plants in Response to Inorganic Phosphate Limitation. Int. J. Mol. Sci. 2018, 19, 3683. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Puga, M.I.; Mateos, I.; Charukesi, R.; Wang, Z.; Franco-Zorrilla, J.M.; de Lorenzo, L.; Irigoyen, M.L.; Masiero, S.; Bustos, R.; Rodríguez, J.; et al. SPX1 is a phosphate-dependent inhibitor of Phosphate Starvation Response 1 in Arabidopsis. Proc. Natl. Acad. Sci. USA 2014, 111, 14947–14952. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Wang, Z.; Ruan, W.; Shi, J.; Zhang, L.; Xiang, D.; Yang, C.; Li, C.; Wu, Z.; Liu, Y.; Yu, Y.; et al. Rice SPX1 and SPX2 inhibit phosphate starvation responses through interacting with PHR2 in a phosphate-dependent manner. Proc. Natl. Acad. Sci. USA 2014, 111, 14953–14958. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Zhong, Y.; Wang, Y.; Guo, J.; Zhu, X.; Shi, J.; He, Q.; Liu, Y.; Wu, Y.; Zhang, L.; Lv, Q.; et al. Rice SPX6 negatively regulates the phosphate starvation response through suppression of the transcription factor PHR2. New Phytol. 2018, 219, 135–148. [Google Scholar] [CrossRef] [Green Version]
  37. Shi, J.; Hu, H.; Zhang, K.; Zhang, W.; Yu, Y.; Wu, Z.; Wu, P. The paralogous SPX3 and SPX5 genes redundantly modulate Pi homeostasis in rice. J. Exp. Bot. 2014, 65, 859–870. [Google Scholar] [CrossRef] [Green Version]
  38. Osorio, M.B.; Ng, S.; Berkowitz, O.; De Clercq, I.; Mao, C.; Shou, H.; Whelan, J.; Jost, R. SPX4 Acts on PHR1-Dependent and -Independent Regulation of Shoot Phosphorus Status in Arabidopsis. Plant Physiol. 2019, 181, 332–352. [Google Scholar] [CrossRef] [Green Version]
  39. Zhao, H.; Gao, Z.; Wang, L.; Wang, J.; Wang, S.; Fei, B.; Chen, C.; Shi, C.; Liu, X.; Zhang, H.; et al. Chromosome-level reference genome and alternative splicing atlas of moso bamboo (Phyllostachys edulis). GigaScience 2018, 7, giy115. [Google Scholar] [CrossRef]
  40. Close, D.C.; Beadle, C.L. Total, and chemical fractions, of nitrogen and phosphorus in Eucalyptus seedling leaves: Effects of species, nursery fertiliser management and transplanting. Plant Soil 2004, 259, 85–95. [Google Scholar] [CrossRef]
  41. Schmittgen, T.D.; Livak, K.J. Schmittgen TD, Livak KJAnalyzing real-time PCR data by the comparative C(T) method. Nat. Protocols 2008, 3, 1101–1108. [Google Scholar] [CrossRef]
  42. Lou, H.Q.; Yang, Y.; Zheng, S.; Ma, Z.M.; Chen, W.J.; Yu, C.L.; Song, L.L.; Wu, J.S. Identification of key genes contributing to amino acid biosynthesis in Torreya grandis using transcriptome and metabolome analysis. Food Chem. 2022, 379, 132078. [Google Scholar] [CrossRef]
Plants 12 01496 g001 550
Figure 1. Chromosomal locations of the SPX-domain-containing proteins family genes in Phyllostachys edulis.
Figure 1. Chromosomal locations of the SPX-domain-containing proteins family genes in Phyllostachys edulis.
Plants 12 01496 g001
Plants 12 01496 g002 550
Figure 2. Phylogenetic analysis of the SPX-domain-containing protein genes in Phyllostachys edulis, Arabidopsis thaliana, and Oryza sativa. The phylogenetic tree was constructed with the neighbor-joining (NJ) algorithm and 1000 bootstrap replicates. Different colors indicate the different subfamilies of SPX-domain-containing proteins.
Figure 2. Phylogenetic analysis of the SPX-domain-containing protein genes in Phyllostachys edulis, Arabidopsis thaliana, and Oryza sativa. The phylogenetic tree was constructed with the neighbor-joining (NJ) algorithm and 1000 bootstrap replicates. Different colors indicate the different subfamilies of SPX-domain-containing proteins.
Plants 12 01496 g002
Plants 12 01496 g003 550
Figure 3. Gene structure, motifs, and conserved structural domains analysis of the PeSPXs in Phyllostachys edulis. (A) Intron-exon structure of the PeSPX genes. Yellow boxes, black lines, and green boxes represent exons (CDS), introns, and the 5′ and 3′ untranslated regions, respectively. (B) Predictions of the conserved domain in 23 PeSPX proteins. The length of each protein sequence is represented by the gray bars, and colored boxes represent conserved domains. (C) Distribution of the conserved motifs in PeSPX proteins. The scale bar at the bottom indicates the protein lengths, and sequence logos for each conserved motif are shown on the right. Visualized by TBtools (v1.108) (https://github.com/CJ-Chen/TBtools/releases (accessed on 12 November 2022)).
Figure 3. Gene structure, motifs, and conserved structural domains analysis of the PeSPXs in Phyllostachys edulis. (A) Intron-exon structure of the PeSPX genes. Yellow boxes, black lines, and green boxes represent exons (CDS), introns, and the 5′ and 3′ untranslated regions, respectively. (B) Predictions of the conserved domain in 23 PeSPX proteins. The length of each protein sequence is represented by the gray bars, and colored boxes represent conserved domains. (C) Distribution of the conserved motifs in PeSPX proteins. The scale bar at the bottom indicates the protein lengths, and sequence logos for each conserved motif are shown on the right. Visualized by TBtools (v1.108) (https://github.com/CJ-Chen/TBtools/releases (accessed on 12 November 2022)).
Plants 12 01496 g003
Plants 12 01496 g004 550
Figure 4. Expression profile of PeSPXs exposed to low phosphorus. (A) Primary root length of Phyllostachys edulis seedlings exposed to low phosphorus. (B) Number of lateral roots of P. edulis seedlings exposed to low phosphorus. (C) Phosphorus content in P. edulis roots and shoots exposed to low phosphorus. (D) Heatmap of PeSPXs exposed to low phosphorus via transcriptome analysis. Error bars in (A–C) indicate SD (n = 3). The asterisk shows a significant difference compared to the control using the unpaired Student’s t test (* p < 0.05).
Figure 4. Expression profile of PeSPXs exposed to low phosphorus. (A) Primary root length of Phyllostachys edulis seedlings exposed to low phosphorus. (B) Number of lateral roots of P. edulis seedlings exposed to low phosphorus. (C) Phosphorus content in P. edulis roots and shoots exposed to low phosphorus. (D) Heatmap of PeSPXs exposed to low phosphorus via transcriptome analysis. Error bars in (A–C) indicate SD (n = 3). The asterisk shows a significant difference compared to the control using the unpaired Student’s t test (* p < 0.05).
Plants 12 01496 g004
Plants 12 01496 g005 550
Figure 5. Expression patterns of PeSPX6 and PeSPX-MFS2 in Phyllostachys edulis. (A) Quantitative RT-PCR analysis of PeSPX6 and PeSPX-MFS2 exposed to low phosphorus in roots. (B) Quantitative RT-PCR analysis of PeSPX6 and PeSPX-MFS2 expression in various tissues (leaf, stem, and root tissues). Error bars in (A,B) indicate SD (n = 3). The asterisks in (A,B) show a significant difference compared to the control using the unpaired Student’s t test (* p < 0.05).
Figure 5. Expression patterns of PeSPX6 and PeSPX-MFS2 in Phyllostachys edulis. (A) Quantitative RT-PCR analysis of PeSPX6 and PeSPX-MFS2 exposed to low phosphorus in roots. (B) Quantitative RT-PCR analysis of PeSPX6 and PeSPX-MFS2 expression in various tissues (leaf, stem, and root tissues). Error bars in (A,B) indicate SD (n = 3). The asterisks in (A,B) show a significant difference compared to the control using the unpaired Student’s t test (* p < 0.05).
Plants 12 01496 g005
Plants 12 01496 g006 550
Figure 6. Validation of the regulatory networks of PeSPX6 and PeSPX-MFS2. (A) Pearson’s correlation analysis between the expression levels of PeSPX6 and PeSPX-MFS2 and transcription factors induced by low phosphorus. (B) Schematic diagrams of the effector and reporter plasmids used for transient expression analysis. (C) Dual-luciferase analysis of the effects of potential transcription factors on LUC activity driven by the PeSPX6 and PeSPX-MFS2 promoter. The asterisk shows a significant difference compared to the control using the unpaired Student’s t test (* p < 0.05).
Figure 6. Validation of the regulatory networks of PeSPX6 and PeSPX-MFS2. (A) Pearson’s correlation analysis between the expression levels of PeSPX6 and PeSPX-MFS2 and transcription factors induced by low phosphorus. (B) Schematic diagrams of the effector and reporter plasmids used for transient expression analysis. (C) Dual-luciferase analysis of the effects of potential transcription factors on LUC activity driven by the PeSPX6 and PeSPX-MFS2 promoter. The asterisk shows a significant difference compared to the control using the unpaired Student’s t test (* p < 0.05).
Plants 12 01496 g006
Table
Table 1. Description of Phyllostachys edulis SPX-domain-containing proteins family genes.
Table 1. Description of Phyllostachys edulis SPX-domain-containing proteins family genes.
Gene IDGene NameLocationORF
(aa)		CDS
(bp)		MW
(KDa)		PiStrand
PH02Gene06029SPX1chr21:3680047-368278431594835.45.44-
PH02Gene33177SPX2chr15:98773223-9877665631695135.765.62+
PH02Gene04435SPX3chr15:32188589-3218952025275928.595.61-
PH02Gene50182SPX4chr21:57199806-5720079825476528.875.67-
PH02Gene21885SPX5chr4:8401540-840256825777429.535.25+
PH02Gene10585SPX6chr10:6035976-603710124573827.975.13+
PH02Gene17759SPX7chr5:30261388-3026370226078329.049.3-
PH02Gene49111SPX8chr20:37686660-3769021728385231.585.37-
PH02Gene12037SPX9chr17:86005776-8600997827783431.455.04-
PH02Gene39430SPX10chr6:41492472-4149733029588833.065.14+
PH02Gene13927SPX11chr8:59290518-5929528729388232.675.26-
PH02Gene14463SPX-MFS1chr8:3039865-3044869701210677.837.5+
PH02Gene13036SPX-MFS2chr6:7519110-7524176699210077.498.43+
PH02Gene00584SPX-MFS3chr24:57835339-57863679692207977.425.86+
PH02Gene01924SPX-MFS4chr23:16516769-16531212695208877.676.14-
PH02Gene01487SPX-MFS5chr3:82291560-82301443686206176.688.43+
PH02Gene03341SPX-MFS6chr17:22492959-22496359696209177.886.65-
PH02Gene41866SPX-EXS1chr6:67258648-67262438741222684.078.89+
PH02Gene41868SPX-EXS2chr6:67340589-67344713757227486.428.98+
PH02Gene10542SPX-EXS3chr8:35123167-35127993856257197.558.74-
PH02Gene19658SPX-EXS4chr14:106783817-106788584804241592.169.25-
PH02Gene33784SPX-EXS5chr3:100920611-100925663824247593.539.28-
PH02Gene05859SPX-EXS6chr17:2650733-2655194789237089.519.4+
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Luo, J.; Liu, Z.; Yan, J.; Shi, W.; Ying, Y. Genome-Wide Identification of SPX Family Genes and Functional Characterization of PeSPX6 and PeSPX-MFS2 in Response to Low Phosphorus in Phyllostachys edulis. Plants 2023, 12, 1496. https://doi.org/10.3390/plants12071496

AMA Style

Luo J, Liu Z, Yan J, Shi W, Ying Y. Genome-Wide Identification of SPX Family Genes and Functional Characterization of PeSPX6 and PeSPX-MFS2 in Response to Low Phosphorus in Phyllostachys edulis. Plants. 2023; 12(7):1496. https://doi.org/10.3390/plants12071496

Chicago/Turabian Style

Luo, Jiali, Zhihui Liu, Jiawen Yan, Wenhui Shi, and Yeqing Ying. 2023. "Genome-Wide Identification of SPX Family Genes and Functional Characterization of PeSPX6 and PeSPX-MFS2 in Response to Low Phosphorus in Phyllostachys edulis" Plants 12, no. 7: 1496. https://doi.org/10.3390/plants12071496

APA Style

Luo, J., Liu, Z., Yan, J., Shi, W., & Ying, Y. (2023). Genome-Wide Identification of SPX Family Genes and Functional Characterization of PeSPX6 and PeSPX-MFS2 in Response to Low Phosphorus in Phyllostachys edulis. Plants, 12(7), 1496. https://doi.org/10.3390/plants12071496

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop