Identification of key genes associated with secondary metabolites biosynthesis by system network analysis in Valeriana officinalis

Valeriana officinalis is a medicinal plant, a source of bioactive chemical compounds and secondary metabolites which are applied in pharmaceutical industries. The advent of ethnomedicine has provided alternatives for disease treatment and has increased demands for natural products and bioactive compounds. A set of preliminary steps to answers for such demands can include integrative omics for systems metabolic engineering, as an approach that contributes to the understanding of cellular metabolic status. There is a growing trend of this approach for genetically engineering metabolic pathways in plant systems, by which natural and synthetic compounds can be produced. As in the case of most medicinal plants, there are no sufficient information about molecular mechanisms involved in the regulation of metabolic pathways in V. officinalis. In this research, systems biology was performed on the RNA-seq transcriptome and metabolome data to find key genes that contribute to the synthesis of major secondary metabolites in V. officinalis. The R Package Weighted Gene Co-Expression Network Analysis (WGCNA) was employed to analyze the data. Based on the results, some major modules and hub genes were identified to be associated with the valuable secondary metabolites. In addition, some TF-encoding genes, including AP2/ERF-ERF, WRKY and NAC TF families, as well as some regulatory factors including protein kinases and transporters were identified. The results showed that several novel hub genes, such as PCMP-H24, RPS24B, ANX1 and PXL1, may play crucial roles in metabolic pathways. The current findings provide an overall insight into the metabolic pathways of V. officinalis and can expand the potential for engineering genome-scale pathways and systems metabolic engineering to increase the production of bioactive compounds by plants.