An Integrated Analysis of WRKY Genes in Autotetraploid Bupleurum chinense: Evolution, Stress Response, and Impact on Saikosaponin Biosynthesis

WRKY transcription factors play critical roles in plant growth, development, metabolism, and stress responses. In this study, we performed the first genome-wide characterization of the WRKY gene family in Bupleurum chinense, using a T2T-level assembly of the autotetraploid genome. A total of 303 BcWRKY genes were identified and found to be unevenly distributed across four subgenomes. Phylogenetic and structural analyses revealed that segmental duplications after polyploidization drove lineage-specific expansion of the family. Meta-transcriptome analysis demonstrated that BcWRKY genes exhibited tissue-specific expression patterns and dynamic responses to stress, suggesting functional diversification. Under drought, waterlogging, methyl jasmonate, and ABA treatments, the contents of saikosaponins A and D significantly increased. This increase was accompanied by transcriptional activation of multiple BcWRKY genes. Correlation analysis between ten BcWRKYs and ten saikosaponins biosynthetic associated genes (BcBASs, BcCYPs, and BcUGTs) identified BcWRKY22, BcWRKY33, and BcWRKY46 as potential regulators of saikosaponin metabolism under stress conditions. Our study provided a comprehensive framework for understanding BcWRKY gene evolution and secondary metabolic regulation in polyploid medicinal plants. It also offered candidate genes for breeding B. chinense cultivars with high saikosaponin content.