Integrated multi-omics analysis provides insights into genome evolution and phosphorus deficiency adaptation in pigeonpea (Cajanus cajan)

Abstract

Pigeonpea (Cajanus cajan) is an important legume food crop and plays a crucial role in a secure food supply in many developing countries. Several previous studies have suggested that pigeonpea has great potential for phosphorus (P) deficiency tolerance, but little is known about the underlying mechanism. In this study, the physiological and molecular responses of pigeonpea roots to phosphate (Pi) starvation were investigated through integrating phenotypic, genomic, transcriptomic, metabolomic, and lipidomic analyses. The results showed that low-Pi treatment increased total root length, root surface area, and root acid phosphatase activity, and promoted the secretion of organic acids (e.g. citric acids, piscidic acids, and protocatechuic acids) and the degradation of phospholipids and other P-containing metabolites in the roots of pigeonpea. Consistent with the morphological, physiological, and biochemical changes, a large number of genes involved in these Pi-starvation responses were significantly upregulated in Pi-deficient pigeonpea roots. Among these Pi-starvation response genes upregulated by low-Pi treatment, four gene families were expanded through recent tandem duplication in the pigeonpea genome, namely phosphate transporter 1 (PHT1), phosphoethanolamine/ phosphocholine phosphatase (PECP), fasciclin-like arabinogalactan protein (FLA), and glutamate decarboxylase (GAD). These gene families may be associated with Pi uptake from the soil, phospholipid recycling, root morphological remodeling, and regulation of organic acid exudation. Taken together, our results suggest that pigeonpea employs complex Pi-starvation responses to strengthen P acquisition and utilization during low-Pi stress. This study provides new insights into the genome evolution and P deficiency adaptation mechanism of pigeonpea.

References

[1] de Bang TC, Husted S, Laursen KH et al. The molecular-physiological functions of mineral macronutrients and their consequences for deficiency symptoms in plants. New Phytol 2021; 229: 2446-69.

[2] Zhu J, Li M, Whelan M . Phosphorus activators contribute to legacy phosphorus availability in agricultural soils: a review. Sci Total Environ 2018; 612: 522-37.

[3] López-Arredondo DL, Leyva-González MA, González-Morales SI et al. Phosphate nutrition: improving low-phosphate tolerance in crops. Annu Rev Plant Biol. 2014; 65: 95-123.

[4] Weeks JJ, Hettiarachchi GM . A review of the latest in phosphorus fertilizer technology: possibilities and pragmatism. J Environ Qual 2019; 48: 1300-13.

[5] Lambers H. Phosphorus acquisition and utilization in plants. Annu Rev Plant Biol 2022; 73: 17-42.

[6] Cong WF, Suriyagoda LDB, Lambers H . Tightening the phosphorus cycle through phosphorus-efficient crop genotypes. Trends Plant Sci. 2020; 25: 967-75.

[7] Liu D. Root developmental responses to phosphorus nutrition. J Integr Plant Biol. 2021; 63: 1065-90.

[8] Dissanayaka DMSB, Ghahremani M, Siebers M et al. Recent insights into the metabolic adaptations of phosphorus-deprived plants. J Exp Bot. 2021; 72: 199-223.

[9] Paz-Ares J, Puga MI, Rojas-Triana M et al. Plant adaptation to low phosphorus availability: core signaling, crosstalks, and applied implications. Mol Plant. 2022; 15: 104-24.

[10] Lan P, Li W, Schmidt W et al. Complementary proteome and transcriptome profiling in phosphate-deficient Arabidopsis roots reveals multiple levels of gene regulation. Mol Cell Proteomics 2012; 11: 1156-66.

[11] Secco D, Jabnoune M, Walker H et al. Spatio-temporal transcript profiling of rice roots and shoots in response to phosphate starvation and recovery. Plant Cell. 2013; 25: 4285-304.

[12] Mo X, Zhang M, Liang C et al. Integration of metabolome and transcriptome analyses highlights soybean roots responding to phosphorus deficiency by modulating phosphorylated metabolite processes. Plant Physiol Biochem 2019; 139: 697-706.

[13] Dokwal D, Cocuron JC, Alonso A et al. Metabolite shift in Medicago truncatula occurs in phosphorus deprivation. J Exp Bot. 2022; 73: 2093-111.

[14] Varshney RK, Chen W, Li Y et al. Draft genome sequence of pigeonpea (Cajanus cajan), an orphan legume crop of resource-poor farmers. Nat Biotechnol 2011; 30: 83-9.

[15] Singh N, Rai V, Singh NK . Multi-omics strategies and prospects to enhance seed quality and nutritional traits in pigeonpea. Nucleus. 2020; 63: 249-56.

[16] Ishikawa S, Adu-Gyamfi JJ, Nakamura T et al. Genotypic variability in phosphorus solubilizing activity of root exudates by pigeonpea grown in low-nutrient environments. Plant Soil 2002; 245: 71-81.

[17] Tian J, Liao H . The role of intracellular and secreted purple acid phosphatases in plant phosphorus scavenging and recycling. In: Plaxton WC, Lambers H (eds). Annual Plant Reviews, Volume 48: Phosphorus Metabolism in Plants.Hoboken: John Wiley & Sons, 2015, 265-87.

[18] Ae N, Arihara J, Okada K et al. Phosphorus uptake by pigeon pea and its role in cropping systems of the Indian subcontinent. Science. 1990; 248: 477-80.

[19] Chai YN, Schachtman DP . Root exudates impact plant performance under abiotic stress. Trends Plant Sci. 2022; 27: 80-91.

[20] Showalter AM, Basu D . Extensin and arabinogalactan-protein biosynthesis: glycosyltransferases, research challenges, and biosensors. Front Plant Sci 2016; 7: 814.

[21] Ashagre HA, Zaltzman D, Idan-Molakandov A et al. FASCICLIN-LIKE 18 is a new player regulating root elongation in Arabidopsis thaliana. Front Plant Sci. 2021; 12: 645286.

[22] Kirchner TW, Niehaus M, Rössig KL et al. Molecular background of Pi deficiency-induced root hair growth in Brassica carinata - a fasciclin-like arabinogalactan protein is involved. Front Plant Sci. 2018; 9: 1372.

[23] Yang B, Li M, Phillips A et al. Nonspecific phospholipase C4 hydrolyzes phosphosphingolipids and sustains plant root growth during phosphate deficiency. Plant Cell 2021; 33: 766-80.

[24] Han Y, White PJ, Cheng L . Mechanisms for improving phosphorus utilization efficiency in plants. Ann Bot. 2022; 129: 247-58.

[25] Hanchi M, Thibaud MC, Légeret B et al. The phosphate fast-responsive genes PECP1 and PPsPase1 affect phosphocholine and phosphoethanolamine content. Plant Physiol. 2018; 176: 2943-62.

[26] Upadhyay N, Kar D, Deepak Mahajan B et al. The multitasking abilities of MATE transporters in plants. J Exp Bot 2019; 70: 4643-56.

[27] Du Z, Su Q, Wu Z et al. Genome-wide characterization of MATE gene family and expression profiles in response to abiotic stresses in rice (Oryza sativa). BMC Ecol Evol 2021; 21: 141.

[28] Khan MIR, Jalil SU, Chopra P et al. Role of GABA in plant growth, development and senescence. Plant Gene 2021; 26: 100283.

[29] Wang P, Dong Y, Zhu L et al. The role of γ-aminobutyric acid in aluminum stress tolerance in a woody plant, Liriodendron chinense × tulipifera. Hortic Res 2021; 8: 80.

[30] Sulieman S, Mühling KH . Utilization of soil organic phosphorus as a strategic approach for sustainable agriculture. J Plant Nutr Soil Sci 2021; 184: 311-9.

[31] Zhu S, Chen M, Liang C et al. Characterization of purple acid phosphatase family and functional analysis of GmPAP7a/7b involved in extracellular ATP utilization in soybean. Front Plant Sci 2020; 11: 661.

[32] Roch GV, Maharajan T, Ceasar SA et al. The role of PHT1 family transporters in the acquisition and redistribution of phosphorus in plants. Crit Rev Plant Sci 2019; 38: 171-98.

[33] Zhou S, Sauvé R, Thannhauser TW . Proteome changes induced by aluminium stress in tomato roots. J Exp Bot 2009; 60: 1849-57.

[34] Murphy J, Riley JP . A modified single solution method for the determination of phosphate in natural waters. Anal Chim Acta 1962; 27: 31-6.

[35] Wu W, Lin Y, Liu P et al. Association of extracellular dNTP utilization with a GmPAP1-like protein identified in cell wall proteomic analysis of soybean roots. J Exp Bot. 2018; 69: 603-17.

[36] Bradford MM . A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem 1976; 72: 248-54.

[37] Liu P, Xue Y, Chen Z et al. Characterization of purple acid phosphatases involved in extracellular dNTP utilization in Stylosanthes. J Exp Bot. 2016; 67: 4141-54.

[38] Tannert M, May A, Ditfe D et al. Pi starvation-dependent regulation of ethanolamine metabolism by phosphoethanolamine phosphatase PECP1 in Arabidopsis roots. J Exp Bot 2018; 69: 467-81.

[39] Kong L, Zhang Y, Ye ZQ et al. CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res 2007; 35: W345-9.

[40] Shen Y, Liu J, Geng H et al. De novo assembly of a Chinese soybean genome. Sci China Life Sci. 2018; 61: 871-84.

[41] Young ND, Debellé F, Oldroyd GE et al. The Medicago genome provides insight into the evolution of rhizobial symbioses. Nature. 2011; 480: 520-4.

[42] Vlasova A, Capella-Gutiérrez S, Rendón-Anaya M et al. Genome and transcriptome analysis of the Mesoamerican common bean and the role of gene duplications in establishing tissue and temporal specialization of genes. Genome Biol. 2016; 17: 32.

[43] Price MN, Dehal PS, Arkin AP . FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol 2009; 26: 1641-50.

[44] Liu C, Wu Y, Liu Y et al. Genome-wide analysis of tandem duplicated genes and their contribution to stress resistance in pigeonpea (Cajanus cajan). Genomics 2021; 113: 728-35.

[45] Wang Y, Tang H, Debarry JD et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res 2012; 40: e49.

[46] Nei M, Gojobori T . Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol 1986; 3: 418-26.

[47] Van K, Kim DH, Cai CM et al. Sequence level analysis of recently duplicated regions in soybean [Glycine max (L.) Merr.] genome. DNA Res. 2008; 15: 93-102.

[48] Chen X, Wang P, Wei M et al. Lipidomics analysis unravels changes from flavor precursors in different processing treatments of purple-leaf tea. J Sci Food Agric 2022; 102: 3730-41.

[49] Yang J, Fu M, Ji C et al. Maize oxalyl-CoA decarboxylase1 degrades oxalate and affects the seed metabolome and nutritional quality. Plant Cell. 2018; 30: 2447-62.

[50] Sun L, Liang C, Chen Z et al. Superior aluminium (Al) tolerance of Stylosanthes is achieved mainly by malate synthesis through an Al-enhanced malic enzyme, SgME1. New Phytol. 2014; 202: 209-19.