Andrographolide, which is enriched in the leaves of Andrographis paniculata, has been known as “natural antibiotic” due to its pharmacological activities such as anti-inflammatory, antimicrobial and antioxidant effects. Several key enzymes in andrographolide biosynthetic pathway have been studied since the genome sequences were released, but its regulatory mechanism remains unknown. WRKY transcription factors proteins have been reported to regulate plant secondary metabolism, development as well as biotic and abiotic stresses. Here, WRKY transcription factors related to andrographolide biosynthesis were systematically identified, including sequences alignment, phylogenetic analysis, chromosomal distribution, gene structure, conserved motifs, synteny, alternative splicing event and Gene ontology (GO) annotation. A total of 58 WRKYs were identified in Chuanxinlian genome and phylogenetically classified into three groups. Moreover, nine WRKY genes underwent alternative splicing events. Furthermore, the combination of binding site prediction, gene-specific expression patterns, and phylogenetic analysis suggested that 7 WRKYs (ApWRKY01, ApWRKY08, ApWRKY12, ApWRKY14, ApWRKY19, ApWRKY20, and ApWRKY50) might regulate andrographolide biosynthesis. This study laid a foundation for understanding the regulatory mechanism of andrographolide biosynthesis and the improvement and breeding of Andrographis paniculata varieties.

Timosaponins, as steroidal saponins, are the primary active constituents and quality biomarkers in Anemarrhena asphodeloides Bunge. Despite their significance, the biosynthetic pathways of timosaponins have not been thoroughly investigated. This study aims to delineate the biosynthetic pathway of timosaponins in A. asphodeloides, elucidate the catalytic mechanisms of the key cycloartenol synthase (CAS), and investigate the antifungal properties of timosaponins. Genes were cloned from A. asphodeloides and heterologous expressed in yeast, tobacco or bacillus coli. Site-directed mutagenesis and molecular docking were used to elucidate the catalytic mechanism of CAS. Antifungal assays were conducted to evaluate the antifungal activities of timosaponins. In this study, we elucidated the biochemical functions of seven genes involved in timosaponins biosynthesis in A. asphodeloides. Among three candidate OSC genes, AaOSCR12 was identified as the gene encoding cycloartenol synthase, which is responsible for the skeleton cyclization in timosaponin biosynthesis. Six residues (257H, 369N, 448T, 507V, 558P, 616Y) were identified as the critical catalytic active sites of CAS (AaOSCR12). Sterol methyltransferase (AaSMT1) and sterol side-chain reductase (AaSSR2) were found to be the subsequent enzymes and the branching points leading to phytosterol and cholesterol biosynthesis, respectively. Two oxide reductase genes, AaCYP90B27 and AaCYP90B2, were responsible for post-modification of cholesterol, serving as a precursor of timosaponins. A key 26-O-β-glucosidase (AaF26G1) was identified as facilitating the conversion of furostanol-type timosaponins into spirostanol-type timosaponins. Antifungal assays revealed that spirostanol-type timosaponin AⅢ exhibits superior antifungal activity compared to furostanol-type timosaponin BⅡ, potentially linked to plant defense mechanisms involving AaF26G1. This study utilized a multi-chassis cross-identification strategy, revealing key enzymes in the timosaponin biosynthetic pathway and offering novel insights into plant defense mechanisms against microbial pathogens.

Bufo bufo gargarizans Cantor are precious medicinal animals in traditional Chinese medicine (TCM). Bufadienolides as the major pharmacological components are generated from the venomous glands of B. bufo gargarizans. Bufadienolides are one type of cardiac aglycone with a six-member lactone ring and have properties of antitumor, cardiotonic, tonsillitis, and anti-inflammatory. The biosynthesis of bufadienolides is complex and unclear. This study explored the transcriptome of three different tissues (skin glands, venom glands, and muscles) of B. bufo gargarizans by high-throughput sequencing. According to the gene tissue–specific expression profile, 389 candidate genes were predicted possibly participating in the bufadienolides biosynthesis pathway. Then, BbgCYP11A1 was identified as a cholesterol side chain cleaving the enzyme in engineering yeast producing cholesterol. Furthermore, the catalytic activity of BbgCYP11A1 was studied with various redox partners. Interestingly, a plant NADPH-cytochrome P450 reductase (CPR) from Anemarrhena asphodeloides showed notably higher production than BbgAdx-2A-BbgAdR from B. bufo gargarizans. These results will provide certainly molecular research to reveal the bufadienolides biosynthesis pathway in B. bufo gargarizans.