Mining and Analysis of Key Genes in the Biosynthetic Pathway of Kinsenoside in Anoectochilus roxburghii

Technical Summary: Mining and Analysis of Key Genes in the Biosynthetic Pathway of Kinsenoside in Anoectochilus roxburghii

Problem Statement
Anoectochilus roxburghii is a valuable medicinal herb where Kinsenoside (KD) serves as the primary bioactive component and a core indicator of medicinal quality. Despite its pharmacological importance, the biosynthetic pathway of KD remains unclear, limiting efforts to enhance production via genetic engineering. Furthermore, while Methyl jasmonate (MeJA) is known to induce secondary metabolite accumulation, the specific regulatory mechanisms of MeJA on KD synthesis in intact A. roxburghii plants, particularly regarding ploidy-specific responses (diploid vs. autotetraploid), have not been elucidated. Existing studies have largely focused on rhizome cultures or metabolic content changes without exploring the underlying molecular mechanisms or ploidy-dependent transcriptional networks.

Methodology
The study employed a multi-omics approach combining physiological induction, transcriptomics, and network analysis:

1. Plant Material and Induction: Diploid and autotetraploid A. roxburghii plants were treated with foliar sprays of MeJA at varying concentrations (50, 100, 500, and 1000 μmol/L). KD content was quantified using UPLC-TQ-MS at multiple time points (1, 3, and 5 days post-spraying) to determine optimal induction conditions for each ploidy.
2. Transcriptome Sequencing (RNA-seq): Based on optimal induction parameters, leaf samples were collected at specific time intervals (diploids: 0, 12, 24, 48 hours; tetraploids: 0, 36, 72, 96 hours). RNA was extracted, libraries constructed, and sequenced on the Illumina NovaSeq 6000 platform.
3. Bioinformatic Analysis:
   • Differential Expression: Differentially expressed genes (DEGs) were identified using DESeq2 with thresholds of |log₂FC| ≥ 2 and FDR < 0.01.
   • Functional Annotation: Genes were annotated against GO, KEGG, and KOG databases.
   • Weighted Gene Co-expression Network Analysis (WGCNA): A gene co-expression network was constructed to correlate gene modules with KD content phenotypes. The module most significantly correlated with KD accumulation was selected for further analysis.
   • Key Gene Screening: Core genes within the target module were identified using the Maximal Clique Centrality (MCC) algorithm in Cytoscape.
4. Validation: The reliability of RNA-seq data was confirmed via quantitative real-time PCR (qRT-PCR) on nine selected DEGs.

Key Results

• Ploidy-Specific Induction: The study established distinct optimal MeJA induction protocols for the two ploidies.
  • Diploids: Optimal at 50 μmol/L MeJA, peaking at 1 day post-spraying with a KD content of 178.92 mg/g (42% increase over control).
  • Tetraploids: Optimal at 500 μmol/L MeJA, peaking at 3 days post-spraying with a KD content of 242.35 mg/g (30% increase over control). High concentrations (1000 μmol/L) inhibited KD synthesis in tetraploids.
• Transcriptional Responses: RNA-seq revealed distinct transcriptional regulation patterns. Diploids exhibited a higher number of DEGs compared to tetraploids. GO and KEGG analyses indicated that diploids were enriched in RNA modification and transcriptional regulation, while tetraploids showed enrichment in abiotic stress responses and DNA-binding transcription factors.
• UGT Identification: Six UDP-glycosyltransferase (UGT) genes responsive to MeJA were identified. Notably, UGT88F3 was significantly upregulated in both diploid and tetraploid lines, suggesting it is a conserved core gene for KD glycosylation.
• WGCNA and Core Gene Mining: WGCNA identified 11 gene co-expression modules. The MEbrown module showed the strongest positive correlation with KD content in tetraploids. Within this module, 15 core genes were screened, including:
  • Transcription Factors: TIFY10B and MYC2 (JA signaling pathway).
  • Enzymes: 2ODD19 (2-oxoglutarate-dependent dioxygenase), OPR7, TLP1, SEOB, CIPK18, and various Cytochrome P450s (CYP74A2, CYP74B2, CYP75B137, CYP92C6, CYP71D7).
  • Other: At1g47380, At3g11180, At5g64030.
• Validation: qRT-PCR results for nine selected genes aligned with RNA-seq expression trends, confirming data reliability.

Significance and Claims
The authors claim this study provides a theoretical basis for dissecting the KD biosynthetic pathway in A. roxburghii. Specifically, the work:

1. Defines Optimal Induction: Establishes the first practical, ploidy-specific foliar spraying protocols for maximizing KD accumulation in intact plants, moving beyond previous rhizome culture limitations.
2. Elucidates Molecular Mechanisms: Reveals that diploid and tetraploid lines possess distinct transcriptional regulatory networks and secondary metabolic activation patterns in response to JA signaling.
3. Identifies Candidate Genes: Successfully screens a set of key candidate genes, including the conserved UGT88F3 and a network of 15 core regulatory genes involving JA signaling and secondary metabolism.
4. Future Utility: The identified genes and induction conditions lay a foundation for subsequent functional verification of the KD pathway and potential genetic engineering strategies to improve KD content for the efficient development of this medicinal plant.

The study maintains a modest tone, framing these findings as a "theoretical basis" and a "preliminary regulatory network" that requires further functional verification, rather than claiming the complete elucidation of the entire biosynthetic pathway.