Skip to main content
Log in

Integrated mRNA and miRNA Expression Analyses of Pinus massoniana Roots and Shoots in Long-Term Response to Phosphate Deficiency

  • Published:
Journal of Plant Growth Regulation Aims and scope Submit manuscript

Abstract

Masson pine (Pinus massoniana) is primarily present in southwest of China, which are severely deficient in inorganic phosphate (Pi). Although some studies identified transcriptomic and proteomic responses to Pi deficiency in Masson pine seedlings, miRNAs and molecular responses in different tissues have not been well studied. To shed further light on the complex responses of Masson pine to Pi deficiency, a spatiotemporal experiment was performed to identify differentially expressed mRNAs and miRNAs under Pi deficiency. Spatiotemporal analyses of 72 RNA sequencing libraries provided a comprehensive overview of the dynamic responses of Masson pine to low-Pi stress. Differentially expressed gene analysis revealed several high-affinity phosphate transporter genes (PHT1-1, PHT1-4, PHT1-5 and PHT1-12) and a nitrate transporter gene (NRT), reflecting the crosstalk between nitrate and Pi homeostasis in plants. miRNA differential expression analysis identified several families that were associated with Pi deficiency, such as miR399. In addition, some other families were dramatically changed in response to Pi starvation, such as miR156, miR169, and some novel miRNAs. Integrated analysis of DE miRNAs and mRNAs indicated that “amino acid metabolism”, “energy metabolism” and “lipid metabolism” were most enriched under Pi deficiency. This study provided essential regulation information between miRNAs and target genes on the response of Masson pine seedlings to Pi deficiency, which will aid in further elucidation of the biological regulatory mechanisms of pines in response to low-Pi stress.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8

Similar content being viewed by others

Abbreviations

P:

Phosphorus

Pi:

Inorganic phosphate

PHR1:

Phosphate starvation response1

SPX:

SYG1/PHO81/XPR1

PHO2:

Phosphate2

PHTs:

Phosphate transporters

PAPs:

Purple acid phosphatases

DEGs:

Differentially expressed genes

ERF:

Ethylene responsive factor

PEPC:

Phosphoenolpyruvate carboxylase

GDPDs:

Glycerophosphodiester phosphodiesterases

SQD1:

UDP-sulfoquinovose synthase

MGDs:

Monogalactosyldiacylglycerol synthases

AP2 :

APETALA2

RSA:

Root system architecture

RPKM:

Reads per kilobase per million reads

PCC:

Pearson correlation coefficient

RT-qPCR:

Reverse transcription quantitative real-time polymerase chain reaction

References

Download references

Funding

This study was funded by the National Key Research and Development Project, P. R. China (2017YFD0600302); National Natural Science Foundation of China (31660185); Science and Technology of Guizhou Province, P. R. China (20175788); Innovation Driven Development Special Fund of Guangxi, P. R. China (AA17204087-4); and the Post-National Key Research and Development Project, P. R. China (20185261).

Author information

Authors and Affiliations

Authors

Contributions

FF conceived, designed the experiments and wrote the paper. XS performed experiments. GD contributed reagents and materials, revised manuscript. ZZ performed experiments. JT collected data. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Fuhua Fan or Guijie Ding.

Ethics declarations

Conflict of interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Additional information

Handling Editor: Mikihisa Umehara.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Figure S1

Length distribution of transcripts using next-generation sequencing. (A)Length distribution of Masson pine unigenes using next-generation sequencing. (B) miRNA length distribution of Masson pine.

Figure S2

Venn diagram of the unigenes annotated by four databases (KEGG, KOG, Nr and Swissprot).

Supplementary file3 (XLSX 19 kb)

Supplementary file4 (XLSX 21 kb)

Supplementary file5 (XLSX 20 kb)

Supplementary file6 (XLSX 156 kb)

Supplementary file7 (XLSX 167 kb)

Supplementary file8 (XLSX 19 kb)

Supplementary file9 (XLSX 17 kb)

Supplementary file10 (XLSX 32 kb)

Supplementary file11 (XLSX 56 kb)

Supplementary file12 (XLSX 32 kb)

Supplementary file13 (XLSX 100 kb)

Supplementary file14 (XLSX 23 kb)

Appendix A: Supplementary Data

Appendix A: Supplementary Data

Illumina reads of all samples have been submitted to the Sequence Read Archive at the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/sra) under Accession Number PRJNA641031.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fan, F., Shang, X., Ding, G. et al. Integrated mRNA and miRNA Expression Analyses of Pinus massoniana Roots and Shoots in Long-Term Response to Phosphate Deficiency. J Plant Growth Regul 41, 2949–2966 (2022). https://doi.org/10.1007/s00344-021-10486-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s00344-021-10486-0

Keywords

Navigation