# The Rice Serine/Arginine Splicing Factor RS33 Regulates Pre-mRNA Splicing during Abiotic Stress Responses

Haroon Butt, Jeremie Bazin, Kasavajhala Prasad, Nourelislam Awad, Martin Crespi, Anireddy Reddy, Magdy Mahfouz

► To cite this version:

Haroon Butt, Jeremie Bazin, Kasavajhala Prasad, Nourelislam Awad, Martin Crespi, et al.. The Rice Serine/Arginine Splicing Factor RS33 Regulates Pre-mRNA Splicing during Abiotic Stress Responses. Cells, 2022, 11 (11), pp.1796. 10.3390/cells11111796. hal-04313167

HAL Id: hal-04313167

https://hal.inrae.fr/hal-04313167v1

Submitted on 17 Sep 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# The Rice Serine/Arginine Splicing Factor RS33 Regulates Pre-mRNA Splicing during Abiotic Stress Responses

Haroon Butt <fnref n="1" />, Jeremie Bazin <fnref n="2" />, Kasavajhala V. S. K. Prasad <fnref n="3" />, Nourelislam Awad <fnref n="4" />, Martin Crespi <fnref n="2" />,
Anireddy S. N. Reddy <fnref n="3" /> and Magdy M. Mahfouz <fnref n="1" />,*

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Butt, H.; Bazin, J.; Prasad, K.V.S.K.; Awad, N.; Crespi, M.; Reddy, A.S.N.; Mahfouz, M.M. The Rice Serine/Arginine Splicing Factor RS33 Regulates Pre-mRNA Splicing during Abiotic Stress Responses.
*Cells* **2022**, *11*, 1796. https://doi.org/10.3390/cells11111796

Academic Editor: Suleyman
Allakhverdiev

Received: 18 March 2022

Accepted: 7 May 2022

Published: 30 May 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Abiotic stresses profoundly affect plant growth and development and limit crop productivity. Pre-mRNA splicing is a major form of gene regulation that helps plants cope with various stresses. Serine/arginine (SR)-rich splicing factors play a key role in pre-mRNA splicing to regulate different biological processes under stress conditions. Alternative splicing (AS) of SR transcripts and other transcripts of stress-responsive genes generates multiple splice isoforms that contribute to protein diversity, modulate gene expression, and affect plant stress tolerance. Here, we investigated the function of the plant-specific SR protein RS33 in regulating pre-mRNA splicing and abiotic stress responses in rice. The loss-of-function mutant *rs33* showed increased sensitivity to salt and low-temperature stresses. Genome-wide analyses of gene expression and splicing in wild-type and *rs33* seedlings subjected to these stresses identified multiple splice isoforms of stress-responsive genes whose AS are regulated by RS33. The number of RS33-regulated genes was much higher under low-temperature stress than under salt stress. Our results suggest that the plant-specific splicing factor RS33 plays a crucial role during plant responses to abiotic stresses.

**Keywords:** pre-mRNA splicing; alternative splicing; SR proteins; genome engineering; abiotic stress

## 1. Introduction

The splicing of pre-mRNAs, a key step in regulating gene expression, helps plants respond to stressful conditions such as salt, drought, and high or low temperatures [1–8]. The spliceosome is a megadalton protein complex that cleaves the introns, ligates the exons, and releases a mature mRNA. Chemical or genetic disruption of spliceosome function inhibits splicing in plants [9–12]. Spliceosome formation requires five small nuclear ribonucleoproteins (snRNPs: U1, U2, U4/U6, and U5) and ~200 additional proteins [13,14]. Splicing factors recognize the *cis*-regulatory elements of their target pre-mRNA and stabilize snRNP assembly across the introns and exons [14–17]. The splicing reaction starts with the binding of the U1 snRNP to the 5' splice site and the U2 auxiliary factor (U2AF) to the 3' splice site. This early complex (complex E) recruits the U2 snRNP to the branch-point (BP) sequence to form a pre-spliceosome complex (complex A). Complex A associates with the preassembled U5-U4/U6 tri-snRNP complex to form the first fully established spliceosome, known as the precatalytic spliceosome (complex B). Complex B is further converted into the activated spliceosome (complex B*) and catalytically activated spliceosome (complex B*), where the

<footnote n="1">¹ Laboratory for Genome Engineering and Synthetic Biology, Division of Biological and Environmental Sciences and Engineering, 4700 King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia; haroon.butt.1@kaust.edu.sa</footnote>

<footnote n="2">² CNRS, INRA, Institute of Plant Sciences Paris-Saclay IP52, University Paris-Saclay and University of Paris Bâtiment 630, 91192 Gif sur Yvette, France; jeremie.bazin@ips2.universite-paris-saclay.fr (J.B.); martin.crespi@ips2.universite-paris-saclay.fr (M.C.)</footnote>

<footnote n="3">³ Department of Biology, Program in Cell and Molecular Biology, Colorado State University, Fort Collins, CO 80523, USA; kasavajhala.prasad@colostate.edu (K.V.S.K.P.); anireddy.reddy@colostate.edu (A.S.N.R.)</footnote>

<footnote n="4">⁴ Helmy Institute of Biomedical Science, Zewail City of Science and Technology, Ahmed Zewail Road, Giza 12578, Egypt; nawad@zewailcity.edu.eg</footnote>

<footnote>* Correspondence: magdy.mahfouz@kaust.edu.sa</footnote>
first transesterification reaction takes place. There are additional distinct, catalytically active spliceosome complexes, C and C*, which catalyze the second transesterification reaction that leads to the excision of the intervening intron and ligation of the adjoining exons [14,17,18].

Serine/arginine (SR)-rich proteins are conserved RNA-binding proteins that function as splicing factors and regulate pre-mRNA splicing by determining the splice sites [18–24]. SR proteins have one or two characteristic N-terminal RNA recognition motifs (RRMs) that determine their RNA-binding properties. The C-terminus of SR proteins is rich in alternating arginine and serine residues, referred to as the RS domain, and this domain facilitates protein-protein interactions [23,25–27]. The C-terminus of SR proteins can be reversibly phosphorylated, which affects the functions of SR proteins, including altering their splicing function [28,29]. In animal systems, the functions of SR proteins are not limited to pre-mRNA splicing and include mRNA nuclear export, mRNA stability, genome maintenance, and translation [19,20,30,31]. Plants have three distinct SR protein subfamilies that are not observed in the animal kingdom, suggesting that SR proteins likely have unique roles during plant growth and development [25,26].

Pre-mRNAs transcribed from multi-exon genes are either constitutively spliced or undergo alternative splicing (AS). AS produces different isoforms with alternative 5' and 3' splice sites, skipped exons, or retained introns. AS thus increases the plasticity of the genome by producing multiple isoforms from a single locus and thereby enhancing the protein diversity [1,4,32]. AS takes place when splicing factors select different splice sites in a closely regulated process [33–35]. Environmental stresses and developmental cues greatly influence the AS patterns of pre-mRNAs in plants [1,7,36–38]. The loss-of-function mutants of different SR family members, such as *rs40*, *rs41*, *scl30a*, and *sr45*, showed hypersensitivity to salt and abscisic acid (ABA) stress in *Arabidopsis thaliana* [39–42]. The pre-mRNAs of SR proteins are also extensively alternatively spliced [43,44], and splicing factor SR45 regulates the salt stress response in an isoform-dependent manner in Arabidopsis [45].

Rice (*Oryza sativa*) SR proteins also play important roles in mineral nutrient homeostasis, and SR40, SCL25, and SCL57 regulate inorganic phosphate (Pi) uptake and mobilization in shoots [46]. Similarly, *SR34b* expression increased under cadmium (Cd) stress, and the *sr34b* mutant showed increased Cd sensitivity [47]. Another SR-like protein, REDUCED RED-LIGHT RESPONSES IN CRY1CRY2 BACKGROUND1 (RRC1), together with SPLICING FACTOR FOR PHYTOCHROME SIGNALING (SFPS), interact with photoreceptor phytochrome B and regulates light-dependent pre-mRNA splicing to promote photomorphogenesis [48]. These studies show that SR proteins regulate AS in response to environmental stresses; however, the precise role of individual SR proteins and the splicing events they regulate are not clear.

The clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 (Cas9) adaptive immunity system of bacteria and archaea has been adapted for genome editing applications [49–54]. In plants, this molecular tool is utilized for trait engineering, directed evolution, and crop improvement [9,55–59]. Recently, we employed the CRISPR-Cas9 system to target all members of the SR gene family and produced single- and multigene knockout mutants in rice [24].

In this study, we investigated the function of RS33, a plant-specific splicing factor, in abiotic stress responses and AS in rice using a loss-of-function mutant. We found that the *rs33* mutation caused hypersensitivity to salt and low-temperature stresses. The RNA deep-sequencing data identified a subset of genes involved in salt and low-temperature stress responses, which are regulated by RS33. Our results provide a link between RS33-regulated AS and plant stress responses.

**2. Materials and Methods**

**2.1. Plant Materials, Rice Transformation, and Genotyping of the Mutant Plants**

*Oryza sativa* L. ssp. *japonica* cv. Nipponbare was used for all the experiments. The *rs33-1* mutant progeny was used for the phenotypic analysis [24]. For CRISPR/Cas9 editing, the expression of *Cas9* was driven by the *OsUBIQUITIN* promoter, and the sgRNAs
were designed to target the 1st exon of the RS33 locus *LOC_Os02g03040*. The sgRNA was synthesized as a polycistronic tRNA-gRNA (PTG) fragment with BsaI overhangs for cloning under the OsU3 promoter in pRGEB32.

*Agrobacterium tumefaciens*-mediated rice transformation was performed using strain EHA105, as described previously [24]. For genotyping, DNA was extracted from leaf samples, and PCR was conducted using gene-specific primers. Purified PCR products were cloned using the CloneJET PCR Cloning Kit (K1231). Sanger sequencing was used to analyze the mutations.

## 2.2. *Phenotypic Analysis*

Rice WT, *Cas9*, and homozygous mutant *rs33* seedlings were sown in a 2:1 mixture of ProMix BX and Green grade in 4-inch pots. The seedlings were allowed to grow until the seed-filling stage under controlled conditions in a growth chamber maintained at 60% relative humidity, temperatures of 27/25 °C, and a light intensity of 400 µmoles/m²/s, with a 16/8-h day/night photoperiod. Freshly harvested seeds of each genotype were used for carrying out stress tolerance assays.

Salt sensitivity assays: To study the effect of salt on seed germination, seeds of individual genotypes were dehusked and surface-sterilized with 70% ethanol for 1 min, followed by a 2.5% bleach solution (with a few drops of Tween 20) for an additional 40 min. Subsequently, the seeds were thoroughly washed five times for 5 min each with sterile deionized water. The sterilized seeds were transferred to a Petri dish with sterile filter paper presoaked in ½-strength Murashige and Skoog (½ MS) medium supplemented with and without 100 mM NaCl, along with 0.5 g/L of MES. The Petri dishes were covered and incubated in a plant growth chamber maintained at 27 °C with 200 µmoles/m²/s of light under a 16/8-h day/night photoperiod. Seed germination was recorded manually every day for the stipulated time. The seeds were considered germinated after the emergence of the radicle. To test the effect of salt stress on seedling growth, the surface-sterilized seeds were placed on filter paper bridges in glass test tubes (20 cm long and 2.5 cm in diameter), each containing 10 mL of ½ MS medium supplemented with and without 100 mM NaCl. The test tubes with seeds were transferred to a growth chamber (under the above-mentioned conditions). The seedlings were allowed to be grown vertically. After the specified time, the seedlings were removed from the test tubes and dry blotted on filter paper, and the fresh weight and shoot length were measured.

Low-temperature sensitivity assays: To test the low-temperature tolerance of these genotypes, surface-sterilized seeds were plated on ~70 mL of ½ MS medium supplemented with 0.5 g/L of MES and 0.8% Phytoblend in a square Petri dish. The Petri dishes were transferred to different growth chambers maintained at the indicated temperatures for a specified time. The Petri dishes were placed vertically to enable the equitable temperature treatment. Seed germination was recorded manually every day. The seeds were considered germinated after the emergence of the radicle.

No seed germination was noted at and below 14 °C. Hence, the seeds that were initially subjected to the low-temperature treatments of 14 and 4 °C were subsequently transferred to 27 °C after 7 and 4 days, respectively, to enable seed germination. To test the effect of low temperature on seed germination and seedling growth, the surface-sterilized seeds were placed on filter paper bridges in glass test tubes, each containing 10 mL of ½ MS medium. The test tubes with seeds were transferred to a growth chamber maintained at the specified temperatures. The seeds were germinated, and the seedlings were allowed to be grown vertically. At the indicated times, the seedlings were removed from the test tubes and dry blotted on filter paper, and the fresh weight and shoot length were measured.

## 2.3. *RNA Isolation and RNA-Seq*

Wild-type and mutant rice seeds were sterilized and transferred to ½ MS agar media (without sucrose) in glass jars for 7 days. For salt stress, 7-day-old seedlings were treated with 200 mM NaCl in ½ MS liquid media or mock-treated without NaCl for 3 and 6 h. For
cold stress, 7-day-old seedlings were treated at 4 °C (in precooled 1/2 MS liquid media) for 3 and 6 h.

After the indicated treatments, total RNA was extracted from the seedlings using the Direct-zol RNA Miniprep Plus kit (Zymo Research, 17062 Murphy Ave., Irvine, CA 92614, USA). RNA samples with an RNA integrity number (RIN) ≥7.0 were selected for library construction and examined through a 2100 Bioanalyzer (Agilent Technologies, 5301 Stevens Creek Blvd. Santa Clara, CA 95051, USA). The TruSeq mRNA stranded kit was used for RNA seq library construction. The HiSeq-4000 platform is used to generate high-quality paired-end reads.

## 2.4. Analysis of RNA-Seq Data and Gene Functional Classification

RNA-seq data were analyzed as previously described [10]. Briefly, Illumina adapters were removed from RNA-seq libraries, and read pairs were trimmed to the same length using Trimmomatic. TrimmomaticPE—LEADING:25, TRAILING:25, CROP:120, and MINLEN:120 arguments were used. Further alignment of the processed reads was done using STAR. The following processing pipeline was used: -outSAMtype BAM SortedByCoordinate -outReadsUnmapped Fastx -readFilesCommand zcat, -quantMode GeneCounts, -outFilterMultimapNmax 1. DEseq2 was used to estimate the significance of the differential gene expression, and FDR was used to correct the *p*-value during the pairwise comparison between genotypes. A gene was counted as differentially expressed if its adjusted *p*-value (FDR) was ≤0.01. Gene ontology enrichment analysis was done as described here [60] and simplified using GO slim terms.

rMATS v4.02 (Multivariate Analysis of Transcript Splicing (MATS): http://rnaseq-mats.sourceforge.net/; accessed on 24 June 2021) was used for differential splicing. The arguments used were: -nthread 4 -read-length 120-t paired -libType fr-firststrand. The gene annotation GTF file described here [46] was used as a reference annotation. Differentially spliced events were defined with FDR < 0.001 and | IncLevelDifference | > 0.

The UpSetR and eulerr R packages were used to compare the sets of DEGs or DAS. The PheatMap package was used for hierarchical clustering, and clusters were extracted using the cutree function. The log2 fold change (log2FC) values of the DEGs and IncLevelDifference of differentially retained introns were used to cluster DEGs and DAS. Lists of all comparisons of DEGs and DAS are listed in Supplementary Tables S2 and S3, respectively.

## 2.5. RT-PCR

DNA digestion of the total RNA was done using an RNase-Free DNase Set (Invitrogen Cat. No. 18068-015). SuperScript First-Strand Synthesis System (Invitrogen, 168 Third Avenue Waltham, MA 02451, USA) was used to generate cDNA from the DNA-digested total RNA. PCR conditions were: initial denaturation at 95 °C for 2 min, then 40 cycles of 95 °C for 30 s, annealing temperature according to the primers for 30 s, and 72 °C for a variable time (according to the amplicon length); then, there was a final elongation at 72 °C for 5 min. The primer sequences used for RT-PCR are listed in Supplementary Table S1.

# 3. Results

## 3.1. OsRS33 Regulates Gene Expression and Pre-mRNA Splicing of a Small Set of Genes

OsRS33, a plant-specific SR protein, contains two RRMs at the N-terminus and one RS domain at the C-terminus (Figure 1a). Using the CRISPR/Cas9 genome editing tool, we targeted the first exon of the OsRS33 locus and recovered a knockout mutant with an 11-nt deletion and a 1-nt substitution (Figure 1a) [24]. The rs33 mutant did not exhibit any obvious growth or developmental defects compared to the wild type (WT), except for a subtle reduction in seedling growth (Figure S1).

To identify RS33-regulated genes, we performed a deep RNA sequencing analysis of seven-day-old WT and rs33 seedlings. A differential gene expression analysis between the WT and rs33 plants revealed a reduced expression of ~150 genes, while ~100 genes were upregulated (Figure 1b). To determine the role of RS33 in regulating pre-mRNA splicing,
we analyzed the genome-wide RNA deep sequencing data to identify altered splicing events in the *rs33* mutant compared with the WT. We identified ~600 differential splicing events between the WT and *rs33* and found a high frequency of intron retention (IR) events, followed by alternative 3' splice sites (A3'S), alternative 5' splice sites (A5'S), and exon skipping (SE) events (Figure 1c).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Genome-wide effects of the loss of OsRS33 on alternative splicing in rice. (a) Targeted mutagenesis of the OsRS33 (LOC_Os02g03040) locus in rice [24]. RS33 is a plant-specific SR protein characterized by the RNA Recognition Motif (RRM) and arginine–serine (RS) repeat domains. The *rs33* knockout mutant has an 11-nt deletion and a 1-nt substitution. (b) The number of up- and downregulated differentially expressed genes (DEGs) between the WT and *rs33*. (c) Bar plot showing the number of each type of differential alternative splicing (DAS) event mediated by RS33; the majority of these AS events were intron retention (IR) events. SE, exon skipping; A5'S, 5' alternative splice site; A3'S, 3' alternative splice site. (d) Venn diagram displaying the overlap between DEGs and DAS genes. (e) cDNA was prepared from one-week-old WT and *rs33* seedlings. RT-PCR was performed using primers that flank introns subject to AS in selected genes. Intron retention was observed in *rs33* mutant rice plants. Arrowheads indicate splicing variants. The gene structures and retained introns are shown. Red boxes indicate the PCR fragments.</figcaption></figure>

To further validate these differential splicing events, we analyzed the IR events of a set of genes whose splicing is regulated by *RS33* using reverse transcription PCR (RT-PCR). We observed the IR of *OsRS29* transcripts specifically in the *rs33* mutant as compared to the WT (Figure 1e). Similarly, the IR events for transcripts of putative pentatricopeptide genes with 949 bp were enriched in the *rs33* background (Figure 1e). A comparison of the differentially expressed genes (DEGs) with the differentially alternatively spliced (DAS) genes revealed that only eight genes were both differentially expressed and differentially spliced (Figure 1d). Our analysis did not show enrichment of any specific biological processes for the DEGs and DAS genes. Overall, our results indicate that *OsRS33* regulates the expression and splicing of a small set of genes during early seedling growth under the control conditions.
3.2. *The Osrs33 Mutant Is Hypersensitive to Salt and Cold Stresses*

To characterize the function of OsRS33 in abiotic stress responses in rice, we subjected WT, transgenic *Cas9* (vector control), and *rs33* plants to salt or low-temperature stress. To examine the effect of salt stress on seed germination and coleoptile greening, WT, *Cas9*, and *rs33* seeds were exposed to 100 mM NaCl for ten days. The germination rates of all three lines were affected by the salt stress. However, *rs33* seeds exhibited significantly enhanced sensitivity to salt stress compared to WT and *Cas9* seeds (indicated by a lower germination rate) (Figure 2a,b). Furthermore, the inhibition of coleoptile greening was also higher in *rs33* seeds (Figure 2a), and coleoptile emergence was significantly delayed in *rs33* seeds compared to WT seeds. Salt stress also significantly reduced the shoot growth and fresh weight of *rs33* seedlings compared to WT and *Cas9* plants (Figure 2c,d).

<figure><img src="image_5.png" /><figcaption>**Figure 2.** The *rs33* mutant is hypersensitive to salt stress. (a,b) Seeds of *rs33* and WT were germinated for ten days in ½ MS media supplemented with 0 mM and 100 mM NaCl. Germination was severely inhibited in *rs33* compared to WT seeds. Other growth parameters—shoot length (c) and seedling fresh weight (d)—were also severely affected in *rs33* as compared to WT seedlings. In c, d; significance is represented by letters. All the variables with the same letters are not statistically significant. If two variables have different letters, they are significantly different.</figcaption></figure>

We also examined the WT, *Cas9*, and *rs33* phenotypes during seed germination under low-temperature stress. Based on our preliminary germination studies, we designed four temperature treatments (27 °C for 9 days (control condition), 20 °C for 9 days, 14 °C for seven days, followed by 27 °C for four days, and 4 °C for four days, followed by 27 °C for seven days).

At 27 °C, WT and *Cas9* seeds exhibited rapid germination within 3 days, while *rs33* mutants seeds exhibited slightly delayed germination (Figure 3a, top left panel), with 100% germination by day 4. At 20 °C, germination was significantly delayed in all three genotypes. For example, over 60% and 50% of WT and *Cas9* seeds, respectively, germinated
by the end of the third day and reached 100% germination on the fifth day, whereas the *rs33* mutant showed only 40% germination on the fifth day and took 9 days to obtain 100% germination (Figure 3a, bottom left panel). When seeds were exposed to 14 °C for seven days and transferred to 27 °C for four days, both WT and *Cas9* seeds exhibited rapid germination (within 1 day) and reached 100% on the 8th day, while *rs33* seeds exhibited delayed germination by 3 days and reached 100% on the 11th day (Figure 3a, top right panel). After exposure to 4 °C for four days and transferring to 27 °C, both WT and *Cas9* seeds exhibited >90% germination within one day at 27 °C. However, the germination of *rs33* seeds was significantly delayed with a germination rate of 10% on day 7, 60% on day 8, and 80% on day 9 (Figure 3a, bottom right panel).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** The *rs33* mutant is hypersensitive to cold stress. (a) The *rs33* mutant and WT seeds were germinated for ten days at the control temperature (27 °C) or different low temperatures (20 °C, 14 °C, and 4 °C). For 14 °C, the seeds were first kept at 14 °C for seven days and then transferred to 27 °C for 4 days. For 4 °C, the seeds were initially treated at 4 °C for four days and then germinated at 27 °C for seven days. At 20 °C, the germination of the *rs33* seeds was significantly delayed compared to that of the WT. The shoot length (b) and seedling fresh weight (c) were significantly affected in *rs33* as compared to WT seedlings. In b, c; all the variables with the same letters are not statistically significant. If two variables have different letters, they are significantly different.</figcaption></figure>

We further analyzed the effect of low temperature on seedling growth by estimating the shoot length and fresh weight. In general, irrespective of the genotype, exposure to the continuous low temperature of 20 °C caused a significant reduction in seedling growth (Figure 3b). However, the *rs33* seedlings exhibited a significant reduction in growth compared to WT and *Cas9* plants, indicating their enhanced sensitivity to low temperatures (Figure 3b). Taken together, our results indicate that the loss of RS33 significantly impacts the salt and cold tolerance of rice seedlings, highlighting the role of OsRS33 in abiotic stress responses.
### 3.3. Loss of OsRS33 Leads to Genome-Wide Splicing Defects and Major Gene Expression Changes under Salt Stress

To further determine the function of OsRS33 during salt stress at the gene expression level, 7-day-old rs33 and WT seedlings were subjected to salt stress for 3 and 6 h, followed by an RNA-seq analysis. After 3 h of salt stress treatment, we identified 148 and 164 DEGs in rs33 and WT plants, respectively (Figure 4a). However, after 6 h of salt stress treatment, there was a slight increase in the number of DEGs in the WT and a large increase in the number of DEGs (2344) in the rs33 mutant (Figure 4a). A majority of these DEGs (2058) in the rs33 mutant were unique and did not overlap with those in the WT (Figure S2a).

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Genome-wide effects of the loss of OsRS33 on gene expression and RNA splicing under salt stress. (a) Number of differentially expressed genes (DEGs) between the rs33 mutant and the WT. One-week-old rice seedlings were treated with 200 mM NaCl for 3 and 6 h and used for RNA extraction. A high number of DEGs was observed for rs33 seedlings after a 6-h treatment with NaCl. (b) Clustering analysis of the DEGs. Cluster 1 clearly shows the downregulation of translation/ribosome biogenesis genes in rs33 in the 6-h NaCl treatment. Cluster 3 shows the upregulation of stress-responsive genes in rs33 in the 6-h salt treatment as compared to the WT. (c) Bar plot showing the number of each type of differential alternative splicing (DAS) event induced by salt treatment in the rs33 mutant and the WT. (d) Heatmap and hierarchical clustering of DAS genes after salt treatment in rs33 and WT. (e) Venn diagram displaying the overlap between DEGs and DAS genes. (f) cDNA was prepared from one-week-old WT and rs33 seedlings treated with 200 mM NaCl for 3 and 6 h. Mock-treated samples were used as controls. RT-PCR was performed using primers that flank introns subject to AS in the selected genes. Arrowheads indicate splicing variants. The gene structures and retained introns are shown. Red boxes indicate the PCR fragments. C3; control 3 h, C6; control 6 h, S3; salt 3 h, and S6; salt 6 h.</figcaption></figure>

Hierarchical clustering revealed the expression profiles of the DEGs in response to salt stress, which were grouped into four co-expression modules (Figure 4b). Strikingly, the largest cluster (cluster 1) corresponded to a set of genes strongly associated with protein translation (padj < 1.6 × 10⁻¹⁴⁸) and associated processes such as ribosome biogenesis (padj < 3.7 × 10⁻¹²¹). This cluster corresponded to genes that were repressed by salt stress
in the *rs33* background as compared to the WT. By contrast, the genes of cluster 3 were specifically induced after 6 h of salt stress in the *rs33* mutant and were enriched for the genes associated with the response to water deprivation (p< 2.41 × 10⁻⁶) and response to ABA (p.adj < 1.9 × 10⁻⁴) (Figure S2b).

To further examine the global defects in pre-mRNA splicing, we used rMATS and defined AS events with FDR < 0.001 and |IncLevelDifference| > 0.2 as DAS. After comparing the control and salt stress conditions, we detected 1158 and 1116 nonredundant DAS events significantly regulated by salt stress in the WT and the *rs33* mutant, respectively. All types of AS events, such as IR, SE, A5'S, and A3'S events, were detected. Among these, differential IR represented the majority of events (Figure 4c). These AS events were mostly specific to the treatment time and the genotype and showed limited overlap (Figure 3a). Clustering of the DAS events in response to salt stress segregated them into five coregulated modules (Figures 4d and S3b). We then compared the DEGs with DAS events in *rs33* and WT plants after 3 h and 6 h of salt treatment (Figure 4e). Most of the DEGs were unique, and very few of them were alternatively spliced (Figure 4e), uncoupling the transcription from AS regulation.

Moreover, we tested the IR of a few selected genes that were alternatively spliced in *rs33* compared to the WT after salt stress treatment by semi-quantitative RT-PCR using primers flanking the retained introns in these target genes. The pre-mRNA splicing factor SF2 was differentially spliced, and isoforms of 500 bp and 592 bp were enriched in the *rs33* mutant after salt treatment (Figure 4f, top left panel). For the helix-loop-helix DNA-binding domain-containing protein, the isoforms of 1807 bp and 3178 bp were differentially enriched between the WT and *rs33* plants (Figure 4f). The presence of both isoforms of *D-MANNOSE-BINDING LECTIN FAMILY PROTEIN* was significantly enhanced in the *rs33* mutant; similarly, the *UBIQUITIN CARBOXYL-TERMINAL HYDROLASE, FAMILY 1* isoform of 582 bp was increased in the *rs33* mutant (Figure 4f) under salt stress. These results indicate that RS33 regulates the splicing patterns of specific genes in response to salt stress.

### 3.4. Low Temperature Differentially Affects Global Gene Expression and AS in *rs33* and WT Plants

In plants, cold stress impacts pre-mRNA splicing and gene expression. To investigate the role of *OsRS33* in cold-regulated gene expression we exposed 7-day-old *rs33* and WT seedlings to low-temperature (4 °C) stress for 3 h and 6 h, then performed RNA-seq and analyzed the data for DEGs and DAS events. Our analysis identified a large number of genes whose expression was regulated by low temperature. The number of DEGs was 1867 for WT-3h, 2153 for *rs33*-3h, 4972 for WT-6h, and 3777 for *rs33*-6h (Figure 5a). Interestingly, at the 3-h time point, *rs33* had a higher number of DEGs, while, at 6 h, the WT had a higher number of DEGs, indicating that the *rs33* mutant was less responsive to prolonged low-temperature stress. The unique DEGs numbered 1984 for WT-6h, 1230 for *rs33*-6h, 203 for WT-3h, and 347 for *rs33*-3h and showed substantial overlap (Figure 4a).

Hierarchical clustering revealed five co-expressed modules, and their functional annotations were associated with photosynthesis and reactive oxygen species metabolic processes (cluster 1), amino acid metabolic process (cluster 2), xylan catabolic process and cell communication (cluster 3), glycoprotein metabolic process, cell wall organization and biogenesis, protein folding (cluster 4), phosphorus metabolic process, xylem development, and phosphorylation (cluster 5) (Figure 4b). Strikingly, cluster 5 showed the upregulation of phosphorylation and phosphorus metabolism genes at both time points in the *rs33* mutant, indicating a strong role of *RS33* in phosphorus metabolism during cold stress (Figures 4b and 5b). Cluster 1 showed a strong upregulation of the oxidative stress response in *rs33* at 6 h of cold treatment as compared to the WT, while cluster 4 showed the downregulation of cell wall biogenesis in *rs33* at 6 h (Figures 4b and 5b).

To further dissect the role of *RS33* during pre-mRNA splicing regulation under cold stress, we performed a splicing analysis and identified DAS events for each genotype and time point. All types of AS events were detected, with IR representing the majority of AS
events (Figure 5c). Our analysis identified IR events numbering 1041 for WT-3h, 1091 for *rs33*-3h, 2148 for WT-6h, and 1868 for *rs33*-6h (Figure 5c). A large set of unique DAS events was identified for each genotype: 874 for WT-6h, 625 for *rs33*-6h, 370 for WT-3h, and 406 for *rs33*-3h (Figure 5a).

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Genome-wide effects of the loss of OsRS33 on gene expression and RNA splicing under low-temperature stress. (a) A number of differentially expressed genes (DEGs) between the *rs33* mutant and the WT. One-week-old rice seedlings were treated with low-temperature stress (4 °C) for 3 and 6 h and used for RNA extraction. (b) Clustering analysis of the DEGs. Cluster 1 shows the upregulation of oxidative stress response genes in *rs33* after 6 h of cold treatment as compared to the WT. Cluster 4 shows that the cell wall biogenesis genes were downregulated after 6 h of cold treatment in *rs33*. (c) Bar plot showing the number of each type of differential alternative splicing (DAS) event induced by low-temperature stress in the *rs33* mutant and the WT. (d) Heatmap and hierarchical clustering of DAS genes after the cold treatment in *rs33* and WT. Cluster 1 contains stress-related genes with quick AS regulation upon stress. Cluster 5 indicates that some genes of the lipid metabolism pathways do not respond to cold in *rs33*. (e) Venn diagram displaying the overlap between DEGs and DAS genes. (f) cDNA was prepared from one-week-old WT and *rs33* seedlings treated at 4 °C for 3 and 6 h. Mock-treated samples were used as the controls. RT-PCR was performed using primers that flank introns subject to AS in the selected genes. Arrowheads indicate splicing variants. The gene structures and retained introns are shown. Red boxes indicate the PCR fragments. C3; control 3 h, C6; control 6 h, F3; 4 °C 3 h, and F6; 4 °C 6 h.</figcaption></figure>

Functional annotation of the individual DAS clusters revealed associations with regulation of the response to stress and lipid oxidation (cluster 1), signal transduction by protein phosphorylation, intracellular signal transduction, the disaccharide biosynthetic process (cluster 2), positive regulation of molecular function and activation of GTPase activity (cluster 3), histone ubiquitination and cellular ketone body metabolic process (cluster 4), and phospholipid metabolic process and dephosphorylation (cluster 5) (Figure 5b). Cluster 4 was upregulated in *rs33* after 3 h and 6 h of cold treatment compared to the WT, while
cluster 5 was downregulated in *rs33* at 3 h and 6 h of cold treatment compared to the WT (Figure 5d).

We compared the DEGs with DAS events, and we found a big overlap between the DEGs and DAS in cold stress, indicating that many of the DEGs were also differentially spliced (Figure 5e). We validated the IR events for a set of genes via semi-quantitative RT-PCR and found an enrichment of introns in the *rs33* mutant seedlings after cold treatment (Figure 5f). We tested several genes with alternatively spliced pre-mRNAs, including those encoding WRKY1 (LOC_Os01g14440), bZIP transcription factor family protein (LOC_Os02g14910), RNA recognition motif-containing protein (LOC_Os03g17760), andcaleosin-related protein/ABA-induced protein (LOC_Os02g50140), to test whether they exhibit intron retention. Our data showed that the *rs33* mutant seedlings had increased IR events during the pre-mRNA splicing of these genes in comparison to WT seedlings (Figure 5f). Overall, these results demonstrate that RS33 regulates cold stress responses at the posttranscriptional level and highlights a link between AS and cold tolerance.

4. Discussion

Exploring adaptive strategies to withstand extreme environmental conditions is needed to develop new plant cultivars that can withstand climate change. Plants have evolved flexible and complex regulatory systems that allow them to adapt to abiotic stresses. AS is a complex, coordinated transcriptional regulatory mechanism, whereby multiple mRNA variants are produced from a single gene via differential intron removal and the retention of different exon combinations from pre-mRNAs [26,33,35]. AS is finely regulated by RNA-binding proteins that recognize sequence signals in RNA and regulate their splicing. SR proteins are RNA-binding proteins that contain one or two characteristic RRMs in the N-terminal portion and an RS domain at the C-terminal end.

Here, we identified the role of a plant-specific SR protein, RS33, with two tandem RRMs, in which one (RRM2) is different from that of human ASF/SF2 [25,61]. Previous studies have shown that RS33 is involved in zinc homeostasis [46]. Additionally, Isshiki et al. were unable to recover transgenic rice plants overexpressing RS33 from repeated transformations, which further emphasized the essential role of this splicing factor during rice growth and development [43]. In our analysis of the *rs33* mutant, we identified only a small set of genes that were differentially expressed or differentially spliced under control conditions (Figure 1). Interestingly, RS29 was alternatively spliced in the *rs33* mutant, suggesting that SR proteins regulate the splicing of other SR members [21,43,62,63]. We also observed IR in transcripts of *LOC_Os08g05750*, which encodes a putative pentatricopeptide (PPR) (Figure 1e). PPRs are targeted at chloroplasts and mitochondria and bind one or several organellar RNA transcripts and, thus, influence their expression [64–66]. A chloroplast retrograde signal was found to regulate the AS of the RS33 homolog in Arabidopsis and is necessary for proper responses to varying light conditions [67]. RS33-mediated regulation of splicing of the chloroplast-localized PPR transcripts may be another level of chloroplast-nucleus communication.

The regulation of gene expression is important for plant survival and acclimation under abiotic stress conditions, and genome-wide studies have shown that pre-mRNA splicing is affected by stress treatments [7,8,68,69]. Splicing inhibition initiates abiotic stress signaling in plants [11,12], and the repression of splicing inhibition can overcome abiotic stress, resulting in stress tolerance [9,10]. The activities of spliceosome proteins and splicing factors are important for plant stress tolerance [1,44,70]. We found that the germination of *rs33* seeds and seedlings was much more sensitive to the 100 mM NaCl treatment than the WT plants (Figure 2). Gene expression analysis indicated that, after 6 h of salt treatment, a large number of genes were differentially expressed between the *rs33* mutant and the WT. A major set of these DEGs was related to translation and RNA biogenesis (cluster 1), which were downregulated in *rs33* plants after 6 h of treatment. Another major portion of DEGs was related to water deprivation and water transport (cluster 3), which were upregulated in *rs33* plants after 6 h of treatment, indicating that RS33 regulates water stress-related
genes under salt stress conditions. However, the DAS genes that were highly expressed after 6 h of salt treatment were mainly related to the regulation of protein modification, histone ubiquitination (cluster 1), and DNA methylation (cluster 2).

We observed a high number of DEGs (~2344) after only 6 h of salt treatment in the *rs33* mutant, but the number of DAS genes in the same treatment was similar to that in the WT. This large difference might be due to the aberrant splicing of some transcription factors in the *rs33* mutant after the salt treatment. We validated the aberrant splicing of some important targets, like pre-mRNA-splicing factor SF2, which has RNA-binding activity, and helix-loop-helix DNA-binding protein, which functions as a transcription factor and is characterized as *IRON-RELATED TRANSCRIPTION FACTOR 2 (IRO2)*. Iron homeostasis is essential for plant growth and development, and *IRO2* expression was elevated in response to glutamine in the regulation of the stress responses [71–73]. The phenotypic analysis indicated that the *IRO2* overexpression lines germinated and grew better than the WT, while the *IRO2-RNAi* lines grew slower than the WT [74]. The D-mannose-binding lectin family protein/S-Domain receptor-like protein-11 (SDRLP-11) is characterized by Ser/Thr kinase domains and induced by salicylic acid [75,76]. Ubiquitin carboxyl-terminal hydrolase family 1/ubiquitin c-terminal hydrolase 5 (UCH5) is an uncharacterized protein. The UCH protein family functions in protein deubiquitination [77]. The Arabidopsis homologs *UCH1* and *UCH2* influence the auxin-dependent developmental pathways, while *UCH3* acts to maintain the period of the circadian clock at high temperatures, redundantly with *UCH1* and *UCH2* [78].

Additionally, we found a number of other genes involved in salt stress responses and are aberrantly spliced in the *rs33* mutant compared with WT under salt stress conditions. For example, *OsSTL1* (*SALT TOLERANCE LEVEL 1*; *LOC_Os04g02000*), a zinc finger family member, was identified as a candidate gene for salt tolerance in a genome-wide association study (GWAS) [79]. Another important protein, *NAC2* (*NAM-AFAT-CUC*; *LOC_Os04g38720*), a member of the NAC transcription factor that regulates many aspects of plant growth and development, was mis-spliced in the *rs33* mutant. The *OsNAC2* expression was enhanced in rice seedlings exposed to high-salt concentrations, and *OsNAC2* accelerates salt-induced programmed cell death [80]. The *NAC2* gene is regulated by the microRNA miR164b. Interestingly, the transgenic plants overexpressing the miR164b-resistant form of *OsNAC2* had higher levels of drought and salt tolerance than the wild type [81]. Similarly, *OsNAC61*, which was aberrantly spliced in the *rs33* mutant, was upregulated after the salt stress treatment in elite rice cultivars [82]. The splicing of *OsGLYII-2* is also affected after the salt treatment in *rs33*. *OsGLYII-2* (*glutathione-responsive rice glyoxalase II-2*; *LOC_Os03g21460*) is the second enzyme of the glyoxalase pathway that detoxifies cytotoxic metabolite methylglyoxal (MG) and belongs to the superfamily of metallo-β-lactamases. *OsGLYII-2* is highly stress inducible and functions in salinity adaptation by maintaining a better photosynthesis efficiency [83,84].

We tested different low-temperature treatments (4 °C for four days, 14 °C for seven days, and 20 °C for nine days) and found that the *rs33* mutant was hypersensitive to these temperatures compared to the WT. We used seedlings treated at 4 °C for 3 h and 6 h for RNA sequencing. Our results revealed a high abundance of DEGs and DAS events in the WT as compared to the *rs33* mutant after cold stress treatment. However, a large number of these genes were unique and did not show overlap (Figures S4 and S5). Photosynthesis and reactive oxygen species metabolic genes were upregulated in *rs33*, while cell wall organization- and glycoprotein metabolic process-related genes were downregulated in *rs33* after 6 h of cold treatment. For splicing, a different set of genes was affected, as the genes responding to stress and involved in lipid oxidation were upregulated, whereas the genes related to protein phosphorylation and signal transduction were downregulated.

We validated certain target genes showing significant changes during splicing. For example, *WRKY1* was alternatively spliced after the low-temperature treatment in the *rs33* mutant. The WRKYS are a family of transcription factors involved in various physiological processes, including biotic and abiotic stress responses and developmental processes. Os-
WRKY62 and OsWRKY76, two genes of the WRKY IIa subfamily, undergo AS in response to pathogen infection [85,86]. The AS event we identified here would generate a truncated iso-form of WRKY1. In Arabidopsis, it was found that the expression of WRKY1 was induced by salt stress, and the *wrky1* mutant conferred increased salt sensitivity [87]. Another transcription factor, *bZIP19*, was alternatively spliced, which added 135 amino acids in-frame to produce a functional protein isoform. The bZIP19 and bZIP23 transcription factors are central regulators of the zinc deficiency response in Arabidopsis [88]. In Arabidopsis, *AtbZIP60*, an ortholog of *OsbZIP50*, was identified to be a target of *INOSITOL-REQUIRING ENZYME-1* (IRE1)-mediated splicing [89,90]. Overall, the loss of *rs33* generates more isoforms of the target genes as compared to WT. This might happen due to the lack of efficient splicing activity or aberrant splicing. In the absence of splicing factors, the spliceosome fails to distinguish the exon-intron junctions of the pre-mRNAs and produce multiple transcripts.

Recent developments in the genome-wide data analysis uncovered the roles of many spliceosome proteins that are important for plant stress tolerance. For example, *SUPERSENSITIVE TO ABSCISIC ACID AND DROUGHT 1 GENE (SAD1)*, which encodes *SM-LIKE PROTEIN 5 (LASM5)*, the *sad1/lsm5* mutation, leads to increased sensitivity to drought stress and ABA [91,92]. Another mutant of the same family, *lsm4*, showed oversensitivity to drought and ABA in Arabidopsis [93]. The Arabidopsis protein U1A controls pre-mRNA splicing under salt stress, and the *atu1a* mutant is hypersensitive to salt stress [94]. An RNA-directed DNA methylation factor, RDM16, is important for splicing, and loss-of-function mutants of RDM16 are hypersensitive to salt and ABA [95]. Similar examples were found for low-temperature stress responses for mutants of spliceosomal proteins. For example, PRP31 regulates the formation of the U4/U6.U5 snRNP complex during splicing. The Arabidopsis *prp31* mutant is hypersensitive to cold stress [96]. Arabidopsis STABILIZED 1, similar to the human U5 snRNP, is necessary for low-temperature tolerance [97,98]. Recently, it was found that a DEAD-Box RNA Helicase protein, OsRH42, plays a critical role under cold stress to regulate pre-mRNA splicing in rice [6]. Our work on RS33 reinforces the link between spliceosomal proteins and stress tolerance and how these SR proteins function to help plants rapidly cope with changing environmental conditions.

**5. Conclusions**

Our work revealed that *rs33* mutant seedlings are hypersensitive to abiotic stresses. The gene expression analysis showed that several genes are differentially spliced in *rs33* plants. Under normal conditions, the loss of *RS33* is irrelevant, possibly due to the existence of other SR proteins. By contrast, under stress conditions, like salt and low temperatures, *RS33* is needed for the splicing of a large number of pre-mRNAs to acclimate to stress. Therefore, in the *rs33* mutant, the absence of *RS33* resulted in aberrant pre-mRNA splicing of a large set of stress response genes and, subsequently, hypersensitivity to abiotic stress. The precise control of splicing factor homeostasis and manipulation of specific splicing machinery genes may be an effective approach to breed plants for tolerance to environmental stresses.

**Supplementary Materials:** The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/cells1111796/s1: Figure S1: The *rs33* mutant phenotype. Figure S2: Genome-wide transcriptional effects of salt stress on the *rs33* mutant. Figure S3: Impact of *OsRS33* on the splicing efficiency under salt stress. Figure S4: Effect of the *rs33* mutation on global gene expression under cold stress. Figure S5: Comparison of global AS events between the WT and the *rs33* mutant under cold stress. Table S1: List of primer sequences used in this study. Table S2: List of DEGs of all comparisons. Table S3: List of DAS of all comparisons.

**Author Contributions:** M.M.M. conceived the project. M.M.M., H.B. and K.V.S.K.P. designed and conducted the experiments. H.B., J.B. and N.A. analyzed the data. H.B., J.B., M.C., A.S.N.R. and M.M.M. designed and analyzed the splicing experiments. M.M.M., H.B., J.B., K.V.S.K.P., M.C. and A.S.N.R. wrote the paper. All authors have read and agreed to the published version of the manuscript.
**Funding:** This work was funded by KAUST baseline funding to Magdy Mahfouz and by the Agri-culture and Food Research Initiative competitive grant (2019-67013-29239) of the USDA National Institute of Food and Agriculture to Anireddy S. N. Reddy.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The RAW data and processed data files were deposited in GEO (GEO; https://www.ncbi.nlm.nih.gov/geo/; accessed on 24 January 2022) under accession number GSE194283. The following secure token was created to allow a review of the record GSE194283, while it remains in a private status: etopyiycdhcllgb. Please note the following points: This token allows anonymous, read-only access to GSE194283 and its associated accessions while they are private. Treat the token as you would a password and realize that the token provides access to GSE194283 to anyone who uses it.

**Acknowledgments:** We would like to thank the members of the genome engineering and synthetic biology laboratory at KAUST for their critical discussion and technical help in this work.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

**References**

1. Laloum, T.; Martin, G.; Duque, P. Alternative Splicing Control of Abiotic Stress Responses. *Trends Plant Sci.* **2018**, *23*, 140–150.
[CrossRef] [PubMed]

2. Bashir, K.; Matsui, A.; Rasheed, S.; Seki, M. Recent advances in the characterization of plant transcriptomes in response to drought, salinity, heat, and cold stress. *F1000Res* **2019**, *8*, F1000 Faculty Rev-658. [CrossRef]

3. Chechanovsky, N.; Hovav, R.; Frenkel, R.; Faigenboim, A.; Eselson, Y.; Petreikov, M.; Moy, M.; Shen, S.; Schaffer, A.A. Low temperature upregulates cwp expression and modifies alternative splicing patterns, increasing the severity of cwp-induced tomato fruit cuticular microfissures. *Hortic. Res.* **2019**, *6*, 122. [CrossRef] [PubMed]

4. Matsui, A.; Nakaminami, K.; Seki, M. Biological Function of Changes in RNA Metabolism in Plant Adaptation to Abiotic Stress. *Plant Cell Physiol.* **2019**, *60*, 1897–1905. [CrossRef] [PubMed]

5. Ule, J.; Blencowe, B.J. Alternative Splicing Regulatory Networks: Functions, Mechanisms, and Evolution. *Mol. Cell* **2019**, *76*, 329–345. [CrossRef] [PubMed]

6. Lu, C.A.; Huang, C.K.; Huang, W.S.; Huang, T.S.; Liu, H.Y.; Chen, Y.F. DEAD-Box RNA Helicase 42 Plays a Critical Role in Pre-mRNA Splicing under Cold Stress. *Plant Physiol.* **2020**, *182*, 255–271. [CrossRef]

7. Ganie, S.A.; Reddy, A.S.N. Stress-Induced Changes in Alternative Splicing Landscape in Rice: Functional Significance of Splice Isoforms in Stress Tolerance. *Biology* **2021**, *10*, 309. [CrossRef]

8. Yu, H.; Du, Q.; Campbell, M.; Yu, B.; Walia, H.; Zhang, C. Genome-wide discovery of natural variation in pre-mRNA splicing and prioritising causal alternative splicing to salt stress response in rice. *New Phytol.* **2021**, *230*, 1273–1287. [CrossRef]

9. Butt, H.; Eid, A.; Momin, A.A.; Bazin, J.; Crespi, M.; Arold, S.T.; Mahfouz, M.M. CRISPR directed evolution of the spliceosome for resistance to splicing inhibitors. *Genome Biol.* **2019**, *20*, 73. [CrossRef]

10. Butt, H.; Bazin, J.; Alshareef, S.; Eid, A.; Benhamed, M.; Reddy, A.S.N.; Crespi, M.; Mahfouz, M.M. Overlapping roles of spliceosomal components SF3B1 and PHF5A in rice splicing regulation. *Commun. Biol.* **2021**, *4*, 529. [CrossRef]

11. AlShareef, S.; Ling, Y.; Butt, H.; Mariappan, K.G.; Benhamed, M.; Mahfouz, M.M. Herboxidiene triggers splicing repression and abiotic stress responses in plants. *BMC Genom.* **2017**, *18*, 260. [CrossRef] [PubMed]

12. Ling, Y.; Alshareef, S.; Butt, H.; Lozano-Juste, J.; Li, L.; Galal, A.A.; Moustafa, A.; Momin, A.A.; Tashkandi, M.; Richardson, D.N.; et al. Pre-mRNA splicing repression triggers abiotic stress signaling in plants. *Plant J.* **2017**, *89*, 291–309. [CrossRef]

13. Wahl, M.C.; Will, C.L.; Luhrmann, R. The Spliceosome: Design Principles of a Dynamic RNP Machine. *Cell* **2009**, *136*, 701–718.
[CrossRef] [PubMed]

14. Kastner, B.; Will, C.L.; Stark, H.; Luhrmann, R. Structural Insights into Nuclear p e-mRNA Splicing in Higher Eukaryotes. *Cold Spring Harb. Perspect. Biol* **2019**, *11*, ARTN a032417. [CrossRef] [PubMed]

15. Kornbliytt, A.R.; Schor, I.E.; Allo, M.; Dujardin, G.; Petrillo, E.; Munoz, M.J. Alternative splicing: A pivotal step between eukaryotic transcription and translation. *Nat. Rev. Mol. Cell Biol.* **2013**, *14*, 153. [CrossRef] [PubMed]

16. Lee, Y.; Rio, D.C. Mechanisms and Regulation of Alternative Pre-mRNA Splicing. *Annu. Rev. Biochem.* **2015**, *84*, 291–323.
[CrossRef] [PubMed]

17. Shi, Y.G. Mechanistic insights into precursor messenger RNA splicing by the spliceosome. *Nat. Rev. Mol. Cell Biol.* **2017**, *18*, 655–670. [CrossRef]

18. Long, J.C.; Caceres, J.F. The SR protein family of splicing factors: Master regulators of gene expression. *Biochem. J.* **2009**, *417*, 15–27. [CrossRef]
19. Shepard, P.J.; Hertel, K.J. The SR protein family. *Genome Biol.* **2009**, *10*, 242. [CrossRef]

20. Twyffels, L.; Gueydan, C.; Kruys, V. Shuttling SR proteins: More than splicing factors. *FEBS J.* **2011**, *278*, 3246–3255. [CrossRef]

21. Bradley, T.; Cook, M.E.; Blanchette, M. SR proteins control a complex network of RNA-processing events. *Rna* **2015**, *21*, 75–92. [CrossRef] [PubMed]

22. Howard, J.M.; Sanford, J.R. The RNAissance family: SR proteins as multifaceted regulators of gene expression. *Wiley Interdiscip. Rev. RNA* **2015**, *6*, 93–110. [CrossRef] [PubMed]

23. Jeong, S. SR Proteins: Binders, Regulators, and Connectors of RNA. *Mol. Cells* **2017**, *40*, 1–9. [CrossRef] [PubMed]

24. Butt, H.; Piatek, A.; Li, L.; Reddy, A.S.N.; Mahfouz, M.M. Multiplex CRISPR Mutagenesis of the Serine/Arginine-Rich (SR) Gene Family in Rice. *Genes* **2019**, *10*, 596. [CrossRef]

25. Barta, A.; Kalyna, M.; Reddy, A.S.N. Implementing a rational and consistent nomenclature for serine/arginine-rich protein splicing factors (SR proteins) in plants. *Plant Cell* **2010**, *22*, 2926–2929. [CrossRef]

26. Morton, M.; AlTamimi, N.; Butt, H.; Reddy, A.S.N.; Mahfouz, M. Serine/Arginine-rich protein family of splicing regulators: New approaches to study splice isoform functions. *Plant Sci.* **2019**, *283*, 127–134. [CrossRef]

27. Melo, J.P.; Kalyna, M.; Duque, P. Current Challenges in Studying Alternative Splicing in Plants: The Case of Physcomitrella patens SR Proteins. *Front. Plant Sci.* **2020**, *11*, 286. [CrossRef]

28. de la Fuente van Bentem, S.; Anrather, D.; Roitinger, E.; Djamei, A.; Hufnagl, T.; Barta, A.; Csaszar, E.; Dohnal, I.; Lecourieux, D.; Hirt, H. Phosphoproteomics reveals extensive in vivo phosphorylation of Arabidopsis proteins involved in RNA metabolism. *Nucleic Acids Res.* **2006**, *34*, 3267–3278. [CrossRef]

29. Carvalho, R.F.; Szakonyi, D.; Simpson, C.G.; Barbosa, I.C.; Brown, J.W.; Baena-Gonzalez, E.; Duque, P. The Arabidopsis SR45 Splicing Factor, a Negative Regulator of Sugar Signaling, Modulates SNF1-Related Protein Kinase 1 Stability. *Plant Cell* **2016**, *28*, 1910–1925. [CrossRef]

30. Huang, Y.; Steitz, J.A. SRprises along a messenger's journey. *Mol. Cell* **2005**, *17*, 613–615. [CrossRef]

31. Wu, H.; Sun, S.Y.; Tu, K.; Gao, Y.; Xie, B.; Krainer, A.R.; Zhu, J. A Splicing-Independent Function of SF2/ASF in MicroRNA Processing. *Mol. Cell* **2010**, *38*, 67–77. [CrossRef] [PubMed]

32. Chaudhary, S.; Jabre, I.; Reddy, A.S.N.; Staiger, D.; Syed, N.H. Perspective on Alternative Splicing and Proteome Complexity in Plants. *Trends Plant Sci.* **2019**, *24*, 496–506. [CrossRef] [PubMed]

33. Staiger, D.; Brown, J.W. Alternative splicing at the intersection of biological timing, development, and stress responses. *Plant Cell* **2013**, *25*, 3640–3656. [CrossRef] [PubMed]

34. Chen, H.C.; Cheng, S.C. Functional roles of protein splicing factors. *Biosci. Rep.* **2012**, *32*, 345–359. [CrossRef] [PubMed]

35. Reddy, A.S.; Marquez, Y.; Kalyna, M.; Barta, A. Complexity of the alternative splicing landscape in plants. *Plant Cell* **2013**, *25*, 3657–3683. [CrossRef]

36. Shang, X.D.; Cao, Y.; Ma, L.G. Alternative Splicing in Plant Genes: A Means of Regulating the Environmental Fitness of Plants. *Int. J. Mol. Sci.* **2017**, *18*, 432. [CrossRef]

37. Szakonyi, D.; Duque, P. Alternative Splicing as a Regulator of Early Plant Development. *Front. Plant Sci.* **2018**, *9*, 1174. [CrossRef]

38. Jabre, I.; Reddy, A.S.N.; Kalyna, M.; Chaudhary, S.; Khokhar, W.; Byrne, L.J.; Wilson, C.M.; Syed, N.H. Does co-transcriptional regulation of alternative splicing mediate plant stress responses? *Nucleic Acids Res.* **2019**, *47*, 2716–2726. [CrossRef]

39. Carvalho, R.F.; Carvalho, S.D.; Duque, P. The plant-specific SR45 protein negatively regulates glucose and ABA signaling during early seedling development in Arabidopsis. *Plant Physiol.* **2010**, *154*, 772–783. [CrossRef]

40. Chen, T.; Cui, P.; Chen, H.; Ali, S.; Zhang, S.; Xiong, L. A KH-domain RNA-binding protein interacts with FIERY2/CTD phosphatase-like 1 and splicing factors and is important for pre-mRNA splicing in Arabidopsis. *PLoS Genet.* **2013**, *9*, e1003875. [CrossRef]

41. Cruz, T.M.; Carvalho, R.F.; Richardson, D.N.; Duque, P. Abscisic acid (ABA) regulation of Arabidopsis SR protein gene expression. *Int. J. Mol. Sci.* **2014**, *15*, 17541–17564. [CrossRef] [PubMed]

42. Xing, D.; Wang, Y.; Hamilton, M.; Ben-Hur, A.; Reddy, A.S. Transcriptome-Wide Identification of RNA Targets of Arabidopsis SERINE/ARGININE-RICH45 Uncovers the Unexpected Roles of This RNA Binding Protein in RNA Processing. *Plant Cell* **2015**, *27*, 3294–3308. [CrossRef] [PubMed]

43. Isshiki, M.; Tsumoto, A.; Shimamoto, K. The serine/arginine-rich protein family in rice plays important roles in constitutive and alternative splicing of pre-mRNA. *Plant Cell* **2006**, *18*, 146–158. [CrossRef] [PubMed]

44. Palusa, S.G.; Ali, G.S.; Reddy, A.S.N. Alternative splicing of pre-mRNAs of Arabidopsis serine/arginine-rich proteins: Regulation by hormones and stresses. *Plant J.* **2007**, *49*, 1091–1107. [CrossRef] [PubMed]

45. Albaqami, M.; Laluk, K.; Reddy, A.S.N. The Arabidopsis splicing regulator SR45 confers salt tolerance in a splice isoform-dependent manner. *Plant Mol. Biol.* **2019**, *100*, 379–390. [CrossRef]

46. Dong, C.; He, F.; Berkowitz, O.; Liu, J.; Cao, P.; Tang, M.; Shi, H.; Wang, W.; Li, Q.; Shen, Z.; et al. Alternative Splicing Plays a Critical Role in Maintaining Mineral Nutrient Homeostasis in Rice (Oryza sativa). *Plant Cell* **2018**, *30*, 2267–2285. [CrossRef] [PubMed]

47. Zhang, W.; Du, B.; Liu, D.; Qi, X. Splicing factor SR34b mutation reduces cadmium tolerance in Arabidopsis by regulating iron-regulated transporter 1 gene. *Biochem. Biophys. Res. Commun.* **2014**, *455*, 312–317. [CrossRef]

48. Xin, R.; Kathare, P.K.; Huq, E. Coordinated Regulation of Pre-mRNA Splicing by the SFPS-RRC1 Complex to Promote Photomorphogenesis. *Plant Cell* **2019**, *31*, 2052–2069. [CrossRef]
49. Haurwitz, R.E.; Jinek, M.; Wiedenheft, B.; Zhou, K.; Doudna, J.A. Sequence- and structure-specific RNA processing by a CRISPR endonuclease. *Science* **2010**, *329*, 1355–1358. [CrossRef]

50. Jinek, M.; Chylinski, K.; Fonfara, I.; Hauer, M.; Doudna, J.A.; Charpentier, E. A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity. *Science* **2012**, *337*, 816–821. [CrossRef]

51. Li, J.F.; Norville, J.E.; Aach, J.; McCormack, M.; Zhang, D.; Bush, J.; Church, G.M.; Sheen, J. Multiplex and homologous recombination-mediated genome editing in Arabidopsis and Nicotiana benthamiana using guide RNA and Cas9. *Nat. Biotechnol.* **2013**, *31*, 688–691. [CrossRef] [PubMed]

52. Mali, P.; Yang, L.; Esvelt, K.M.; Aach, J.; Guell, M.; DiCarlo, J.E.; Norville, J.E.; Church, G.M. RNA-guided human genome engineering via Cas9. *Science* **2013**, *339*, 823–826. [CrossRef] [PubMed]

53. Butt, H.; Eid, A.; Ali, Z.; Atia, M.A.M.; Mokhtar, M.M.; Hassan, N.; Lee, C.M.; Bao, G.; Mahfouz, M.M. Efficient CRISPR/Cas9-Mediated Genome Editing Using a Chimeric Single-Guide RNA Molecule. *Front. Plant Sci.* **2017**, *8*, 1441. [CrossRef]

54. Butt, H.; Rao, G.S.; Sedeek, K.; Aman, R.; Kamel, R.; Mahfouz, M. Engineering herbicide resistance via prime editing in rice. *Plant Biotechnol. J.* **2020**, *18*, 2370. [CrossRef]

55. Butt, H.; Jamil, M.; Wang, J.Y.; Al-Babili, S.; Mahfouz, M. Engineering plant architecture via CRISPR/Cas9-mediated alteration of strigolactone biosynthesis. *BMC Plant Biol.* **2018**, *18*, 174. [CrossRef] [PubMed]

56. Butt, H.; Zaidi, S.S.; Hassan, N.; Mahfouz, M. CRISPR-Based Directed Evolution for Crop Improvement. *Trends Biotechnol.* **2020**, *38*, 236–240. [CrossRef] [PubMed]

57. Shan, Q.; Wang, Y.; Li, J.; Zhang, Y.; Chen, K.; Liang, Z.; Zhang, K.; Liu, J.; Xi, J.J.; Qiu, J.L.; et al. Targeted genome modification of crop plants using a CRISPR-Cas system. *Nat. Biotechnol.* **2013**, *31*, 686–688. [CrossRef]

58. Zong, Y.; Wang, Y.; Li, C.; Zhang, R.; Chen, K.; Ran, Y.; Qiu, J.L.; Wang, D.; Gao, C. Precise base editing in rice, wheat and maize with a Cas9-cytidine deaminase fusion. *Nat. Biotechnol.* **2017**, *35*, 438–440. [CrossRef]

59. Rao, G.S.; Jiang, W.; Mahfouz, M. Synthetic directed evolution in plants: Unlocking trait engineering and improvement. *Synth. Biol.* **2021**, *6*, ysab025. [CrossRef]

60. Horan, K.; Jang, C.; Bailey-Serres, J.; Mittler, R.; Shelton, C.; Harper, J.F.; Zhu, J.K.; Cushman, J.C.; Gollery, M.; Girke, T. Annotating genes of known and unknown function by large-scale coexpression analysis. *Plant Physiol.* **2008**, *147*, 41–57. [CrossRef]

61. Manley, J.L.; Krainer, A.R. A rational nomenclature for serine/arginine-rich protein splicing factors (SR proteins). *Genes Dev.* **2010**, *24*, 1073–1074. [CrossRef] [PubMed]

62. Ali, G.S.; Palusa, S.G.; Golovkin, M.; Prasad, J.; Manley, J.L.; Reddy, A.S.N. Regulation of Plant Developmental Processes by a Novel Splicing Factor. *PLoS ONE* **2007**, *2*, e471. [CrossRef] [PubMed]

63. Zhang, X.N.; Mount, S.M. Two alternatively spliced isoforms of the Arabidopsis SR45 protein have distinct roles during normal plant development. *Plant Physiol.* **2009**, *150*, 1450–1458. [CrossRef] [PubMed]

64. Saha, D.; Prasad, A.M.; Srinivasan, R. Pentatricopeptide repeat proteins and their emerging roles in plants. *Plant Physiol. Biochem.* **2007**, *45*, 521–534. [CrossRef] [PubMed]

65. Barkan, A.; Small, I. Pentatricopeptide repeat proteins in plants. *Annu. Rev. Plant Biol.* **2014**, *65*, 415–442. [CrossRef]

66. Manna, S. An overview of pentatricopeptide repeat proteins and their applications. *Biochimie* **2015**, *113*, 93–99. [CrossRef] [PubMed]

67. Petrillo, E.; Godoy Herz, M.A.; Fuchs, A.; Reifer, D.; Fuller, J.; Yanovsky, M.J.; Simpson, C.; Brown, J.W.; Barta, A.; Kalyna, M.; et al. A chloroplast retrograde signal regulates nuclear alternative splicing. *Science* **2014**, *344*, 427–430. [CrossRef]

68. Filichkin, S.A.; Priest, H.D.; Givan, S.A.; Shen, R.; Bryant, D.W.; Fox, S.E.; Wong, W.K.; Mockler, T.C. Genome-wide mapping of alternative splicing in Arabidopsis thaliana. *Genome Res.* **2010**, *20*, 45–58. [CrossRef]

69. Ding, F.; Cui, P.; Wang, Z.Y.; Zhang, S.D.; Ali, S.; Xiong, L.M. Genome-wide analysis of alternative splicing of pre-mRNA under salt stress in Arabidopsis. *BMC Genom.* **2014**, *15*, 431. [CrossRef]

70. Tanabe, N.; Yoshimura, K.; Kimura, A.; Yabuta, Y.; Shigeoka, S. Differential expression of alternatively spliced mRNAs of Arabidopsis SR protein homologs, atSR30 and atSR45a, in response to environmental stress. *Plant Cell Physiol.* **2007**, *48*, 1036–1049. [CrossRef]

71. Kan, C.C.; Chung, T.Y.; Juo, Y.A.; Hsieh, M.H. Glutamine rapidly induces the expression of key transcription factor genes involved in nitrogen and stress responses in rice roots. *BMC Genom.* **2015**, *16*, 731. [CrossRef] [PubMed]

72. Wang, S.; Li, L.; Ying, Y.; Wang, J.; Shao, J.F.; Yamaji, N.; Whelan, J.; Ma, J.F.; Shou, H. A transcription factor OsbHLH156 regulates Strategy II iron acquisition through localising IRO2 to the nucleus in rice. *New Phytol.* **2020**, *225*, 1247–1260. [CrossRef] [PubMed]

73. Wang, W.; Ye, J.; Ma, Y.; Wang, T.; Shou, H.; Zheng, L. OsIRO3 Plays an Essential Role in Iron Deficiency Responses and Regulates Iron Homeostasis in Rice. *Plants* **2020**, *9*, 1095. [CrossRef] [PubMed]

74. Ogo, Y.; Itai, R.N.; Kobayashi, T.; Aung, M.S.; Nakanishi, H.; Nishizawa, N.K. OsIRO2 is responsible for iron utilization in rice and improves growth and yield in calcareous soil. *Plant Mol. Biol.* **2011**, *75*, 593–605. [CrossRef]

75. Blanco, F.; Salinas, P.; Cecchini, N.M.; Jordana, X.; Van Hummelen, P.; Alvarez, M.E.; Holuigue, L. Early genomic responses to salicylic acid in Arabidopsis. *Plant Mol. Biol.* **2009**, *70*, 79–102. [CrossRef]

76. Naithani, S.; Dikeman, D.; Garg, P.; Al-Bader, N.; Jaiswal, P. Beyond gene ontology (GO): Using biocuration approach to improve the gene nomenclature and functional annotation of rice S-domain kinase subfamily. *PeerJ* **2021**, *9*, e11052. [CrossRef]
77. Wang, D.H.; Song, W.; Wei, S.W.; Zheng, Y.F.; Chen, Z.S.; Han, J.D.; Zhang, H.T.; Luo, J.C.; Qin, Y.M.; Xu, Z.H.; et al. Characterization of the Ubiquitin C-Terminal Hydrolase and Ubiquitin-Specific Protease Families in Rice (Oryza sativa). *Front. Plant Sci.* **2018**, *9*, 1636. [CrossRef]

78. Hayama, R.; Yang, P.; Valverde, F.; Mizoguchi, T.; Furutani-Hayama, I.; Vierstra, R.D.; Coupland, G. Ubiquitin carboxyl-terminal hydrolases are required for period maintenance of the circadian clock at high temperature in Arabidopsis. *Sci. Rep.* **2019**, *9*, 17030. [CrossRef]

79. Yuan, J.; Wang, X.; Zhao, Y.; Khan, N.U.; Zhao, Z.; Zhang, Y.; Wen, X.; Tang, F.; Wang, F.; Li, Z. Genetic basis and identification of candidate genes for salt tolerance in rice by GWAS. *Sci. Rep.* **2020**, *10*, 9958. [CrossRef]

80. Mao, C.; Ding, J.; Zhang, B.; Xi, D.; Ming, F. OsNAC2 positively affects salt-induced cell death and binds to the OsAP37 and OsCOX11 promoters. *Plant J.* **2018**, *94*, 454–468. [CrossRef]

81. Jiang, D.; Zhou, L.; Chen, W.; Ye, N.; Xia, J.; Zhuang, C. Overexpression of a microRNA-targeted NAC transcription factor improves drought and salt tolerance in Rice via ABA-mediated pathways. *Rice* **2019**, *12*, 76. [CrossRef] [PubMed]

82. Garcia-Morales, S.; Gomez-Merino, F.C.; Trejo-Tellez, L.I. NAC transcription factor expression, amino acid concentration and growth of elite rice cultivars upon salt stress. *Acta Physiol. Plant* **2014**, *36*, 1927–1936. [CrossRef]

83. Mustafiz, A.; Singh, A.K.; Pareek, A.; Sopory, S.K.; Singla-Pareek, S.L. Genome-wide analysis of rice and Arabidopsis identifies two glyoxalase genes that are highly expressed in abiotic stresses. *Funct. Integr. Genom.* **2011**, *11*, 293–305. [CrossRef]

84. Ghosh, A.; Pareek, A.; Sopory, S.K.; Singla-Pareek, S.L. A glutathione responsive rice glyoxalase II, OsGLYII-2, functions in salinity adaptation by maintaining better photosynthesis efficiency and anti-oxidant pool. *Plant J.* **2014**, *80*, 93–105. [CrossRef] [PubMed]

85. Wu, K.L.; Guo, Z.J.; Wang, H.H.; Li, J. The WRKY family of transcription factors in rice and Arabidopsis and their origins. *DNA Res.* **2005**, *12*, 9–26. [CrossRef]

86. Liu, J.; Chen, X.; Liang, X.; Zhou, X.; Yang, F.; Liu, J.; He, S.Y.; Guo, Z. Alternative Splicing of Rice WRKY62 and WRKY76 Transcription Factor Genes in Pathogen Defense. *Plant Physiol.* **2016**, *171*, 1427–1442. [CrossRef]

87. Wu, X.; Xu, J.; Meng, X.; Fang, X.; Xia, M.; Zhang, J.; Cao, S.; Fan, T. Linker histone variant HIS1-3 and WRKY1 Oppositely Regulate Salt Stress Tolerance in Arabidopsis. *Plant Physiol.* **2022**. [CrossRef] [PubMed]

88. Assuncao, A.G.; Herrero, E.; Lin, Y.F.; Huettel, B.; Talukdar, S.; Smaczniak, C.; Immink, R.G.; van Eldik, M.; Fiers, M.; Schat, H.; et al. Arabidopsis thaliana transcription factors bZIP19 and bZIP23 regulate the adaptation to zinc deficiency. *Proc. Natl. Acad. Sci. USA* **2010**, *107*, 10296–10301. [CrossRef]

89. Deng, Y.; Humbert, S.; Liu, J.X.; Srivastava, R.; Rothstein, S.J.; Howell, S.H. Heat induces the splicing by IRE1 of a mRNA encoding a transcription factor involved in the unfolded protein response in Arabidopsis. *Proc. Natl. Acad. Sci. USA* **2011**, *108*, 7247–7252. [CrossRef]

90. Nagashima, Y.; Mishiba, K.; Suzuki, E.; Shimada, Y.; Iwata, Y.; Koizumi, N. Arabidopsis IRE1 catalyses unconventional splicing of bZIP60 mRNA to produce the active transcription factor. *Sci. Rep.* **2011**, *1*, 29. [CrossRef]

91. Xiong, L.; Gong, Z.; Rock, C.D.; Subramanian, S.; Guo, Y.; Xu, W.; Galbraith, D.; Zhu, J.K. Modulation of abscisic acid signal transduction and biosynthesis by an Sm-like protein in Arabidopsis. *Dev. Cell* **2001**, *1*, 771–781. [CrossRef]

92. Cui, P.; Zhang, S.; Ding, F.; Ali, S.; Xiong, L. Dynamic regulation of genome-wide pre-mRNA splicing and stress tolerance by the Sm-like protein LSm5 in Arabidopsis. *Genome Biol.* **2014**, *15*, R1. [CrossRef] [PubMed]

93. Zhang, Z.; Zhang, S.; Zhang, Y.; Wang, X.; Li, D.; Li, Q.; Yue, M.; Li, Q.; Zhang, Y.E.; Xu, Y.; et al. Arabidopsis floral initiator SKB1 confers high salt tolerance by regulating transcription and pre-mRNA splicing through altering histone H4R3 and small nuclear ribonucleoprotein LSM4 methylation. *Plant Cell* **2011**, *23*, 396–411. [CrossRef]

94. Gu, J.; Xia, Z.; Luo, Y.; Jiang, X.; Qian, B.; Xie, H.; Zhu, J.K.; Xiong, L.; Zhu, J.; Wang, Z.Y. Spliceosomal protein U1A is involved in alternative splicing and salt stress tolerance in Arabidopsis thaliana. *Nucleic Acids Res.* **2018**, *46*, 1777–1792. [CrossRef] [PubMed]

95. Huang, C.F.; Miki, D.; Tang, K.; Zhou, H.R.; Zheng, Z.; Chen, W.; Ma, Z.Y.; Yang, L.; Zhang, H.; Liu, R.; et al. A Pre-mRNA-splicing factor is required for RNA-directed DNA methylation in Arabidopsis. *PLoS Genet.* **2013**, *9*, e1003779. [CrossRef] [PubMed]

96. Du, J.L.; Zhang, S.W.; Huang, H.W.; Cai, T.; Li, L.; Chen, S.; He, X.J. The Splicing Factor PRP31 Is Involved in Transcriptional Gene Silencing and Stress Response in Arabidopsis. *Mol. Plant* **2015**, *8*, 1053–1068. [CrossRef] [PubMed]

97. Lee, B.H.; Kapoor, A.; Zhu, J.; Zhu, J.K. STABILIZED1, a stress-upregulated nuclear protein, is required for pre-mRNA splicing, mRNA turnover, and stress tolerance in Arabidopsis. *Plant Cell* **2006**, *18*, 1736–1749. [CrossRef]

98. Kim, G.D.; Cho, Y.H.; Lee, B.H.; Yoo, S.D. STABILIZED1 Modulates Pre-mRNA Splicing for Thermotolerance. *Plant Physiol.* **2017**, *173*, 2370–2382. [CrossRef]