Identification of core genes associated with different phosphorus levels in quinoa seedlings by weighted gene co-expression network analysis

Abstract

Background: Quinoa is a highly nutritious and novel crop that is resistant to various abiotic stresses. However, its growth and development is restricted due to its limited utilization of soil phosphorus. Studies on the levels of phosphorus in quinoa seedlings are limited; therefore, we analyzed transcriptome data from quinoa seedlings treated with different concentrations of phosphorus. Results: To identify core genes involved in responding to various phosphorus levels, the weighted gene co-expression network analysis method was applied. From the 12,085 expressed genes, an analysis of the gene co-expression network was done. dividing the expressed genes into a total of twenty-five different modules out of which two modules were strongly correlated with phosphorus levels. Subsequently we identified five core genes that correlated strongly either positively or negatively with the phosphorus levels. Gene ontology and assessments of the Kyoto Encyclopedia of Genes and Genomes have uncovered important biological processes and metabolic pathways that are involved in the phosphorus level response. Conclusions: We discovered crucial new core genes that encode proteins from various transcription factor families, such as MYB, WRKY, and ERF, which are crucial for abiotic stress resistance. This new library of candidate genes associated with the phosphorus level responses in quinoa seedlings will help in breeding varieties that are tolerant to phosphorus levels.