Large-scale analysis of the cassava transcriptome reveals the impact of cold stress on alternative splicing

Alternative splicing is an essential post-transcriptional regulatory mechanism that can impact mRNA stability and protein diversity of eukaryotic genomes. Although numerous forms of stress-responsive alternative splicing have been identified in model plants, a large-scale study of alternative splicing dynamics under abiotic stress conditions in cassava has not been conducted. Here, we report the parallel employment of isoform-Seq, ssRNA-Seq, and Degradome-Seq to investigate the diversity, abundance, and fate of alternatively spliced isoforms in response to cold and drought stress. We identified 38 164 alternative splicing events, among which 3292 and 1025 events were significantly regulated by cold and drought stress, respectively. Intron retention was the most abundant subtype of alternative splicing. Global analysis of splicing regulators revealed that the number of their alternatively spliced isoforms and the corresponding abundance were specifically modulated by cold stress. We found that 58.5% of cold-regulated alternative splicing events introduced a premature termination codon into the transcripts, and 77.6% of differential alternative splicing events were detected by Degradome-Seq. Our data reveal that cold intensely affects both quantitative and qualitative aspects of gene expression via alternative splicing pathways, and advances our understanding of the high complexity and specificity of gene regulation in response to abiotic stresses.