Uncovering the molecular mechanisms behind steroidal saponin accumulation in Liriope muscari (Decne.) Baily through transcriptome sequencing and bioinformatics analysis

Abstract

The leaves and roots of Liriope muscari (Decne.) Baily were subjected to high-throughput Illumina transcriptome sequencing. Bioinformatics analysis was used to investigate the enzyme genes and key transcription factors involved in regulating the accumulation of steroidal saponins, which are the main active ingredient in L. muscari. These analyses aimed to reveal the molecular mechanism behind steroidal saponin accumulation. The sequencing results of L. muscari revealed 31 enzymes, including AACT, CAS, DXS and DXR, that are involved in the synthesis of steroidal saponins. Among these enzymes, 16 were in the synthesis of terpenoid skeleton, 3 were involved in the synthesis of sesquiterpene and triterpene, and 12 were involved in the synthesis of steroidal compound. Differential gene expression identified 15 metabolic enzymes coded by 34 differentially expressed genes (DEGs) in the leaves and roots, which were associated with steroidal saponin synthesis. Further analysis using gene co-expression patterns showed that 14 metabolic enzymes coded by 31 DEGs were co-expressed. In addition, analysis using gene co-expression analysis and PlantTFDB's transcription factor analysis tool predicted the involvement of 8 transcription factors, including GAI, PIF4, PIL6, ERF8, SVP, LHCA4, NF-YB3 and DOF2.4, in regulating 6 metabolic enzymes such as DXS, DXR, HMGR, DHCR7, DHCR24, and CAS. These eight transcription factors were predicted to play important roles in regulating steroidal saponin accumulation in L. muscari. Promoter analysis of these transcription factors indicated that their main regulatory mechanisms involve processes such as abscisic acid response, drought-induction stress response and light response, especially abscisic acid responsive elements (ABRE) response and MYB binding site involved in drought-inducibility (MBS) response pathway. Furthermore, qRT-PCR analysis of these eight key transcription factors demonstrated their specific differences in the leaves and roots.