A global survey of bicarbonate stress-induced pre-mRNA alternative splicing in soybean via integrative analysis of Iso-seq and RNA-seq

Abstract

Alternative splicing (AS) plays important roles in modulating environmental stress responses in plants. However, little is known about the functions of bicarbonate-induced AS in cultivated soybean (Glycine max L. Merr.). In this study, we combined PacBio isoform sequencing (Iso-seq) and Illumina RNA sequencing (RNA-seq) to elucidate the bicarbonate-induced AS events in soybean root and leaf tissues. Compared to RNA-seq, Iso-seq identified more novel genes and transcripts, as well as more AS events, indicating that Iso-seq is more efficient in AS detection. Combining these two technologies, we found that intron retention (IR) is the most frequent AS event type. We identified a total of 913 and 1974 bicarbonate stress-responsive differentially alternative spliced genes (DAGs) in soybean leaves and roots respectively, from our RNA-seq results. Additionally, we determined a transcription factor (GmNTL9) and a splicing factor (GmRSZ22), and validated their roles in bicarbonate stress response by AS. Overall, our study opens an avenue for evaluating plant AS regulatory networks, and the obtained global landscape of alternative splicing provides valuable insights into the AS-mediated bicarbonate-responsive mechanisms in plant species.