Comprehensive transcriptomic analysis of BjGL1-knockout Brassica juncea: novel insights into leaf trichome formation and defense-associated biological processes

Technical Summary: Comprehensive Transcriptomic Analysis of BjGL1-Knockout Brassica juncea

Problem Statement
Brassica juncea (mustard) is a vital crop for oil, condiments, and medicine, yet it faces significant challenges from biotic stresses, particularly aphid infestations. Leaf trichomes serve as a primary physical and chemical barrier against such pests. While previous research established that the simultaneous knockout of two homologous genes, BjA06.GL1 and BjB02.GL1, results in a complete loss of leaf trichomes and reduced aphid resistance, the downstream molecular regulatory networks and the broader biological processes influenced by these genes remain poorly understood. Specifically, there is a lack of transcriptomic data elucidating how the absence of these genes alters gene expression related to trichome formation, secondary metabolite biosynthesis, and defense responses in gene-edited mustard.

Methodology
The study employed a multi-omics approach combining CRISPR/Cas9 genome editing with RNA sequencing (RNA-seq) and quantitative real-time PCR (qRT-PCR) validation:

1. Plant Material and Gene Editing: The researchers utilized CRISPR/Cas9 to simultaneously target the first and third exons of BjA06.GL1 and BjB02.GL1 in B. juncea. T₁ generation plants with double-knockout mutations (designated BG) exhibiting a trichome-deficient phenotype were compared against a wild-type control (BT) with normal trichome density.
2. Genotype and Off-Target Verification: Genomic DNA was extracted for PCR amplification and Sanger sequencing to confirm mutations. Off-target effects were rigorously screened using BLAST alignment against the B. juncea reference genome, followed by PCR and sequencing of potential off-target loci, confirming high editing specificity with no detected off-target mutations.
3. Transcriptome Profiling: RNA was extracted from 21-day-old leaves of both BT and BG lines (three biological replicates each). Libraries were sequenced on an Illumina NovaSeq 6000 platform. Clean reads were aligned to the B. juncea (Xuecai) reference genome using HISAT2, and transcript assembly was performed with StringTie.
4. Differential Expression Analysis: Differentially expressed genes (DEGs) were identified using DESeq2 with thresholds of |log₂ fold change| ≥ 1 and adjusted P-value < 0.05. Gene expression was quantified as FPKM.
5. Functional Annotation and Validation: DEGs were subjected to Gene Ontology (GO) and KEGG pathway enrichment analyses. Additionally, K-means clustering was used to group genes by expression patterns. Key DEGs, including transcription factors and metabolic pathway genes, were validated via qRT-PCR.

Key Contributions and Results

• Phenotypic Confirmation: The double-knockout BG lines exhibited a complete absence of leaf trichomes on both adaxial and abaxial surfaces, confirmed by scanning electron microscopy (SEM). The mutations involved frameshifts caused by deletions and insertions at the target sites, stably inherited in the T₁ generation.
• Transcriptomic Landscape: The study identified 4,604 differentially expressed genes (DEGs) between the wild-type (BT) and knockout (BG) lines. Specifically, 1,831 genes were upregulated and 2,773 were downregulated in the trichome-bearing BT leaves compared to the BG mutant.
• Enriched Pathways: Functional enrichment analysis revealed that DEGs were significantly associated with:
  • Cytochrome P450 (CYP) metabolism: Crucial for the oxidative transformation of endogenous precursors and biosynthesis of defensive secondary metabolites.
  • Plant–pathogen interactions: Suggesting a direct link between trichome presence and immune response mechanisms.
  • Secondary metabolite biosynthesis: Including flavonoid and phenylpropanoid pathways, which are vital for plant defense and stress adaptation.
  • Environmental adaptation and transmembrane transport.
• Transcription Factor Networks: The analysis identified 357 differentially expressed transcription factors (TFs), with ERF, MYB, and NAC families being the most abundant.
  • Positive Regulators: Core genes such as BjGL1, BjGL2, BjZFP8, BjSPL9, and BjMYB106 were significantly upregulated in BT.
  • Negative Regulators: Genes known to suppress trichome development, such as BjTCP10, BjSPL5, and R3-type MYB genes (BjTRY, BjCPC, BjETC1), showed distinct expression patterns, with repressors generally downregulated or less active in the context of the active regulatory network in BT.
  • Notably, while four GL1 homologs exist in the B. juncea genome, only BjA06.GL1 and BjB02.GL1 showed significant differential expression; the other two homologs (BjuLA09G063230 and BjuLB07G024530) were not differentially expressed.
• Validation: qRT-PCR results confirmed the expression trends of key regulatory genes identified in the RNA-seq data, validating the reliability of the transcriptomic findings.

Significance and Claims
The authors claim that this study provides a critical theoretical basis for understanding the molecular mechanisms governing trichome-mediated resistance in Brassica juncea. The primary significance lies in demonstrating that BjA06.GL1 and BjB02.GL1 are not merely structural determinants of trichome density but are central regulators that coordinate a broad transcriptional network. This network extends beyond simple morphogenesis to include:

1. Defense Mechanisms: The upregulation of plant-pathogen interaction genes and secondary metabolite pathways in trichome-rich plants suggests that trichome development is intrinsically linked to the plant's chemical defense arsenal.
2. Metabolic Regulation: The findings indicate that these genes influence the phenylpropanoid pathway and cytochrome P450 metabolism, potentially altering the plant's chemical profile to enhance resistance against biotic and abiotic stresses.
3. Breeding Applications: By elucidating the specific regulatory genes and pathways involved, the study offers valuable targets for molecular breeding programs aimed at developing B. juncea cultivars with enhanced insect resistance and reduced reliance on chemical pesticides.

The paper concludes that while the study successfully maps the transcriptomic shifts associated with BjGL1 knockout, further research integrating metabolomics is necessary to fully clarify the downstream metabolic changes and the precise ecological roles of these regulatory networks.