Integrative physiological, transcriptome and metabolome analysis reveals the involvement of carbon and flavonoid biosynthesis in low phosphorus tolerance in cotton

Phosphorus (P) is an essential nutrient controlling plant growth and development through the regulation of basic metabolic processes; however, the molecular details of these pathways remain largely unknown. In this study, physiological, transcriptome, and metabolome analysis were compared for two cotton genotypes with different low P tolerance under P starvation and resupply. The results showed that the glucose, fructose, sucrose, and starch contents increased by 18.2%, 20.4%, 20.2%, and 14.3% in the roots and 18.3%, 23.3%, 11.0%, and 13.6% in the shoot of Jimian169 than DES926, respectively. Moreover, the activities of enzymes related to carbon and phosphorus metabolism were higher in the roots and shoots of Jimian169 than DES926. In addition, transcriptome analysis revealed that the number of differentially expressed genes (DEGs) was higher in both roots (830) and shoots (730) under P starvation and the DEGs drastically reduced upon P resupply. The KEGG analysis indicated that DEGs were mainly enriched in phenylpropanoid biosynthesis, carbon metabolism, and photosynthesis. The metabolome analysis showed the enrichment of phenylpropanoid, organic acids and derivatives, and lipids in all the pairs at a given time point. The combined transcriptome and metabolome analysis revealed that carbon metabolism and flavonoid biosynthesis are involved in the P starvation response in cotton. Moreover, co-expression network analysis identified 3 hub genes in the roots and shoots that regulate the pathways involved in the P starvation response. This study provides the foundation for understanding the mechanisms of low P tolerance and the hub genes as a potential target for the development of low P tolerant genotypes.