Salinity stress significantly affects the contents of bioactive constituents in licorice . To elucidate the molecular mechanism underlying the difference in the accumulation of these constituents under sodium chloride (NaCl, salt) stress, licorice seedlings were treated with NaCl and then subjected to an integrated transcriptomic and metabolite profiling analysis. The transcriptomic analysis results identified 3,664 differentially expressed genes (DEGs) including transcription factor family MYB and basic helix-loop-helix (). Most DEGs were involved in flavonoid and terpenoid biosynthesis pathways. In addition, 121 compounds including a triterpenoid and five classes of flavonoids (isoflavone, flavone, flavanone, isoflavan, and chalcone) were identified, and their relative levels were compared between the stressed and control groups using data from the ultrafast liquid chromatography (UFLC)-triple quadrupole-time of flight-tandem mass spectrometry (TOF-MS/MS) analysis. Putative biosynthesis networks of the flavonoids and triterpenoids were created and combined with structural DEGs such as phenylalanine ammonia-lyase (), 4-coumarate-CoA ligase [], cinnamate 4-hydroxylase [], chalcone synthase [], chalcone-flavanone isomerase [], and flavonoid-3',5' hydroxylase ('') for flavonoids, and and for glycyrrhizin biosynthesis. Notably, significant upregulation of UDP-glycosyltransferase genes () in salt-stressed licorice indicated that postmodification of glycosyltransferase may participate in downstream biosynthesis of flavonoid glycosides and triterpenoid saponins. Accordingly, the expression trend of the DEGs is positively correlated with the accumulation of glycosides. Our study findings indicate that key DEGs and crucial genes co-regulate flavonoid and saponin biosynthesis in licorice under salt stress.
