Alternative splicing shapes the transcriptome complexity in blackgram [Vigna mungo (L.) Hepper]

Vigna mungo, a highly consumed crop in the pan-Asian countries, is vulnerable to several biotic and abiotic stresses. Understanding the post-transcriptional gene regulatory cascades, especially alternative splicing (AS), may underpin large-scale genetic improvements to develop stress-resilient varieties. Herein, a transcriptome based approach was undertaken to decipher the genome-wide AS landscape and splicing dynamics in order to establish the intricacies of their functional interactions in various tissues and stresses. RNA sequencing followed by high-throughput computational analyses identified 54,526 AS events involving 15,506 AS genes that generated 57,405 transcripts isoforms. Enrichment analysis revealed their involvement in diverse regulatory functions and demonstrated that transcription factors are splicing-intensive, splice variants of which are expressed differentially across tissues and environmental cues. Increased expression of a splicing regulator NHP2L1/SNU13 was found to co-occur with lower intron retention events. The host transcriptome is significantly impacted by differential isoform expression of 1172 and 765 AS genes that resulted in 1227 (46.8% up and 53.2% downregulated) and 831 (47.5% up and 52.5% downregulated) transcript isoforms under viral pathogenesis and Fe²⁺ stressed condition, respectively. However, genes experiencing AS operate differently from the differentially expressed genes, suggesting AS is a unique and independent mode of regulatory mechanism. Therefore, it can be inferred that AS mediates a crucial regulatory role across tissues and stressful situations and the results would provide an invaluable resource for future endeavours in V. mungo genomics.