Alternative splicing (AS) is a post-transcriptional level of gene expression regulation that increases transcriptome and proteome diversity. How the AS landscape of rapeseed (Brassica napus L.) changes in response to the fungal pathogen Sclerotinia sclerotiorum is unknown. Here, we analyzed 18 RNA-seq libraries of mock-inoculated and S. sclerotiorum-inoculated susceptible and tolerant B. napus plants. We found that infection increased AS, with intron retention being the main AS event. To determine the key genes functioning in the AS response, we performed a differential AS (DAS) analysis. We identified 79 DAS genes, including those encoding splicing factors, defense response proteins, crucial transcription factors and enzymes. We generated coexpression networks based on the splicing isoforms, rather than the genes, to explore the genes’ diverse functions. Using this weighted gene coexpression network analysis alongside a gene ontology enrichment analysis, we identified 11 modules putatively involved in the pathogen defense response. Within these regulatory modules, six DAS genes (ascorbate peroxidase 1, ser/arg-rich protein 34a, unknown function 1138, nitrilase 2, v-atpase f, and amino acid transporter 1) were considered to encode key isoforms involved in the defense response. This study provides insight into the post-transcriptional response of B. napus to S. sclerotiorum infection.
Alternative splicing (AS) is a post-transcriptional regulatory process that enhances transcriptome diversity, thereby affecting plant growth, development, and stress responses. To identify the new transcripts and changes in the isoform-level AS landscape of rapeseed (Brassica napus) infected with the fungal pathogen Leptosphaeria maculans, we compared eight RNA-seq libraries prepared from mock-inoculated and inoculated B. napus cotyledons and stems. The AS events that occurred in stems were almost the same as those in cotyledons, with intron retention representing the most common AS pattern. We identified 1892 differentially spliced genes between inoculated and uninoculated plants. We performed a weighted gene co-expression network analysis (WGCNA) to identify eight co-expression modules and their Hub genes, which are the genes most connected with other genes within each module. There are nine Hub genes, encoding nine transcription factors, which represent key regulators of each module, including members of the NAC, WRKY, TRAF, AP2/ERF-ERF, C2H2, C2C2-GATA, HMG, bHLH, and C2C2-CO-like families. Finally, 52 and 117 alternatively spliced genes in cotyledons and stems were also differentially expressed between mock-infected and infected materials, such as HMG and C2C2-Dof; which have dual regulatory mechanisms in response to L. maculans. The splicing of the candidate genes identified in this study could be exploited to improve resistance to L. maculans.