# Overlapping roles of spliceosomal components SF3B1 and PHF5A in rice splicing regulation

Haroon Butt, Jeremie Bazin, Sahar Alshareef, Ayman Eid, Moussa Benhamed, Anireddy S N Reddy, Martin Crespi, Magdy M Mahfouz

► To cite this version:

Haroon Butt, Jeremie Bazin, Sahar Alshareef, Ayman Eid, Moussa Benhamed, et al.. Overlapping roles of spliceosomal components SF3B1 and PHF5A in rice splicing regulation. Communications Biology, 2021, 4 (1), 10.1038/s42003-021-02051-y. hal-03331718

HAL Id: hal-03331718

https://hal.inrae.fr/hal-03331718

Submitted on 2 Sep 2021

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
<figure><img src="image_2.png" /></figure>

# Overlapping roles of spliceosomal components SF3B1 and PHF5A in rice splicing regulation

Haroon Butt ID <fnref n="1" />,⁴, Jeremie Bazin²,⁴, Sahar Alshareef ID <fnref n="1" />, Ayman Eid ID <fnref n="1" />, Moussa Benhamed², Anireddy S. N. Reddy<fnref n="3" />, Martin Crespi² & Magdy M. Mahfouz ID <fnref n="1" />✉

The SF3B complex, a multiprotein component of the U2 snRNP of the spliceosome, plays a crucial role in recognizing branch point sequence and facilitates spliceosome assembly and activation. Several chemicals that bind SF3B1 and PHF5A subunits of the SF3B complex inhibit splicing. We recently generated a splicing inhibitor-resistant SF3B1 mutant named *SF3B1 GEX1A RESISTANT 4* (SGR4) using CRISPR-mediated directed evolution, whereas splicing inhibitor-resistant mutant of *PHF5A* (*Overexpression-PHF5A GEX1A Resistance, OGR*) was generated by expressing an engineered version PHF5A-Y36C. Global analysis of splicing in wild type and these two mutants revealed the role of SF3B1 and PHF5A in splicing regulation. This analysis uncovered a set of genes whose intron retention is regulated by both proteins. Further analysis of these retained introns revealed that they are shorter, have a higher GC content, and contain shorter and weaker polypyrimidine tracts. Furthermore, splicing inhibition increased seedlings sensitivity to salt stress, consistent with emerging roles of splicing regulation in stress responses. In summary, we uncovered the functions of two members of the plant branch point recognition complex. The novel strategies described here should be broadly applicable in elucidating functions of splicing regulators, especially in studying the functions of redundant paralogs in plants.

<footnote n="1">¹Laboratory for Genome Engineering and Synthetic Biology, King Abdullah, University of Science and Technology (KAUST), Thuwal, Saudi Arabia. ²CNRS, INRA, Institute of Plant Sciences Paris-Saclay IPS2, Univ Paris Sud, Univ Evry, Univ Paris-Diderot, Sorbonne Paris-Cite, Universite Paris-Saclay, Orsay, France.</footnote>

<footnote n="3">³Department of Biology and Program in Cell and Molecular Biology, Colorado State University, Fort Collins, CO, USA. ⁴These authors contributed equally: Haroon Butt, Jeremie Bazin. ✉email: Magdy.mahfouz@kaust.edu.sa</footnote>
RNA splicing is an essential process in eukaryotes whereby the spliceosome excises the intron sequences and ligates the exon sequences together to produce mature mRNAs. The spliceosome is a macromolecular megaDalton machine comprising of small nuclear RNA molecules (snRNAs U1, U2, U4, U5, and U6) and associated proteins. The spliceosome undergoes a complex assembly process, performs two sequential translocation reactions, and is regulated in response to cellular, developmental, and environmental cues¹,². The pre-mRNA contains consensus sequences to identify exon/intron junctions and cis-sequences to guide the splicing machinery and other trans factors for successful splicing. Conserved cis-splicing consensus sequences include the 5' and 3' splice sites (SS), a polypyrimidine tract and the branch point sequence (BP), which are required for recognition and spliceosome assembly via numerous RNA-RNA, RNA-protein, and protein-protein interactions. Moreover, the spliceosome undergoes dynamic conformational rearrangements to mediate the splicing process³,⁴. Splicing must exhibit high fidelity as a single base pair shift in the SS can lead to a coding sequence frameshift that ultimately disrupts translation. Thus, aberrant splicing causes many genetic diseases in humans⁵ and affects a plethora of stress responses in plants⁶. For example, erroneous recognition and selection of the BP can lead to the accumulation of aberrant transcripts in the cell, potentially leading to cancerous phenotypes⁷,⁸. Alternative splicing (AS), which involves the inclusion or skipping of exonic or intronic sequences and differential selection of 5' or 3' Ss, enhances mRNA coding capacity⁹,¹⁰ and also regulates the levels of functional mRNAs. AS constitutes a co-/post-transcriptional regulatory layer that is activated by growth, developmental, and environmental cues to ensure appropriate molecular responses¹⁰,¹¹.

The SF3B complex of the U2 snRNP contains several spliceosome-associated proteins (SAPs), namely SF3B1/SAP155, SF3B2/SAP145, SF3B3/SAP130, SF3B4/SAP49, SF3B5/SAP10, SF3B6/SAP14a, and PHF5A/SF3B7/SAP14b, and this heptameric complex is essential for splicing⁸,¹². The SF3B complex is mainly involved in recognition of branch point adenosine (BPA) and promotion of stable interaction of U2 with pre-mRNA¹³. This complex is a target of natural compounds with antitumor activities, including Herboxidiene (GEX1A)¹⁴, Spliceostatin A (SSA)¹⁵, and Pladienolide B (PB)¹⁶,¹⁷. These compounds have anticancer activities due to their preferential cytotoxic effects on cancer cells. Structural and biochemical studies have shown that these compounds bind to the SF3B1-PHF5A complex, occupying the pocket of the SF3B complex and blocking SF3B complex binding to the BP¹⁸⁻²⁰. As a result, these splicing modulators perturb BP recognition and selection, preventing the stable formation of the U2 snRNP-BP duplex¹⁸,²¹. These splicing modulators share common binding sites in SF3B1 and PHF5A proximal to the binding pocket of the BP¹⁸,¹⁹,²². SF3B1 plays a role in the recruitment or stabilization of the U2 snRNP complex at the BP sequence during the formation of the spliceosome A complex²³. SF3B1 (SF3b155) is the largest subunit of the SF3B complex and contains an unstructured N-terminal domain (NTD) and C-terminal HEAT repeat domain (HD). The NTD and HD play essential roles at different stages of the splicing process²⁴,²⁵. The splicing factor PHF5A, characterized by its PHD-domain, is a ubiquitously expressed nuclear protein that is highly conserved among eukaryotes²⁶. PHF5A plays an essential role during embryo formation and tissue morphogenesis by modulating pluripotency and cellular differentiation²⁷. Dysregulation of PHF5A has different effects on the progression of different cancers²⁸⁻³⁰. PHF5A is a component of the SF3B complex, essential for spliceosome structural stability, and can link the spliceosome to histones¹⁸,²⁹. Point mutations in key residues of the SF3B complex subunits, including SF3B1-K1071E, SF3B1-

R1074H, SF3B1-V1078A/I, and PHF5A-Y36C, confer tolerance to chemical compounds that affect splicing¹⁸⁻²⁰,²².

In plants, chemical modulators of splicing such as GEX1A and PB affect the splicing machinery resulting in genome-wide inhibition of splicing³¹,³². Although they have considerable effects on all modes of splicing, these compounds show some level of specificity for stress-related genes. Furthermore, in plants, functions of SF3B1 and PHF5A subunits of the SF3B complex, the known targets of inhibitors of splicing are largely unknown. Therefore, the use of splicing modulators may assist in understanding the functions of SF3B1 and PHF5A, molecular mechanisms underlying stress responses and splicing regulation in plants³¹,³². Recently, we employed the CRISPR-Cas9 system for targeted engineering of the splicing protein SF3B1 in rice (Oryza sativa) under selective pressure³³. In this CRISPR-directed evolution approach, we recovered different SF3B1 variants capable of conferring tolerance to the splicing inhibitors GEX1A and PB³³. The mutant lines carrying these SF3B1 variants were termed SGR (SF3B1 GEX1A Resistant). However, the global impact of these SF3B1 mutant variants on gene expression and splicing, as well as the molecular responses of these mutant variants to splicing modulators, were not analyzed in our previous work. Similarly, our understanding of the molecular function of PHF5A, a SF3B1 interactor, in splicing is based primarily on studies in mammalian cell lines. Hence, the roles of PHF5A and SF3B1 in splicing regulation remains largely unknown in plants.

In this study, we determined the molecular function and physiological roles of these two proteins of the branch point recognition complex in plants. We report the detailed phenotypic and molecular analyses of the SF3B1 mutant variant SGR4, insensitive to splicing modulators. Compared with WT plants, SGR4 did not exhibit disturbed pre-mRNA splicing under splicing inhibition by GEX1A. Moreover, we also engineered OsPHF5A to become resistant to the splicing inhibitory drug and showed that heterologous expression of PHF5A-Y36C in rice confers tolerance to splicing modulators. Global analysis of splicing in wild-type and these two mutants in the presence and absence of a splicing inhibitor revealed the role of SF3B1 and PHF5A in splicing regulation and its impact on rice stress responses. We discovered that the retained introns associated with the inhibition SF3B1 and PHF5A activity are shorter, have higher GC content, and have shorter and weaker polypyrimidine tracts. The GO terms enriched under splicing inhibition conditions are mainly 'response to chemical' and 'response to stress'. Furthermore, splicing inhibition increased seedlings sensitivity to salt stress. Collectively, our results uncovered the functions of two members of the branch point recognition complex. These novel approaches should be largely useful in revealing functions of splicing regulators and to study the role of redundant homologs in plants under normal and stress conditions.

## Results

**SGR4 displays insensitivity to the splicing-inhibitor GEX1A.** The SF3B1 protein has U2AF65 interaction and SF3B14 interaction domains in the N-terminal region and HEAT repeat domain (HD) and CTD domains in the C-terminal region (Fig. 1a). SGR4 was generated using CRISPR-mediated directed evolution platform and carries K1049R, K1050E, G1051H substitutions (Fig. 1a). We have previously shown that SGR4 is resistant to the GEX1A³³. To investigate the effect of GEX1A on the growth and development of SGR4, we conducted a detailed phenotypic analysis. We applied different concentrations of GEX1A to WT and the SGR4 and observed the effects on seed germination and seedling growth. Our analysis indicated that the germination of SGR4 is not affected even at 10 µM GEX1A while
<figure><img src="image_3.png" /><figcaption>**Fig. 1 SGR4 exhibits insensitivity to splicing inhibition under GEX1A treatment.** a SGR4, an SF3B1 mutant, was evolved through the CRISPR-based directed evolution (CDE) platform, has coding changes of c3144t, a3145c, a3146g, g3147t, a3148g, g3150a, g3151c, g3152a, t3153c and AA changes K1049R, K1050E, G1051H. b, c Seeds of SGR4 and WT were germinated for five days in water supplemented with 0, 5, and 10 µM GEX1A. The emergence of plumule was considered as germination. The germination is not only delayed but the shoot and root growth is severely inhibited in WT seeds. SGR4 germination was not affected by GEX1A treatment. d, e Seedlings of SGR4 and WT were germinated for three days on ½ MS media and then transferred to ½ MS media supplemented with 0.3 µM GEX1A for a further three days. Root tips were marked after transfer to new media to observe growth. The root growth of WT seedlings is seized while SGR4 is unaffected. (n = 5). f Seedlings were germinated for three days and later transferred to ½ MS media supplemented with 0.5 µM or 1 µM GEX1A. Lateral root growth was observed under a microscope. Lateral roots of WT are severely inhibited whereas the roots of SGR4 are insensitive. g cDNA was prepared from 1-week-old rice seedlings following treatment with 0.3 µM GEX1A for 6-h. RT-PCR was performed using primers that flank introns subject to AS in selected genes. Intron inhibition is observed in WT rice plants under GEX1A treatment. Arrowheads indicate splicing variants that changed following GEX1A treatment. The gene structures and retained introns are shown. Red boxes indicate the PCR fragments. M Mock treatment, G GEX1A treatment.</figcaption></figure>

WT germination is severely inhibited at 5 µM GEX1A (Fig. 1b, c). Consistent with the germination assays, SGR4 has a sustained primary root length in the presence of 0.3 µM GEX1A whereas WT was completely arrested (Fig. 1d, e). Next, we investigated the effect of the GEX1A splicing modulator on lateral root growth in WT and SGR4. We conducted a lateral root assay using 0.5 µM and 1 µM GEX1A treatments. The WT plants exhibited sensitivity to 0.5 µM and 1 µM GEX1A treatments, leading to inhibition of lateral root formation, confirming that splicing regulation is an important component of LR formation and development. However, SGR4 shows increased LR density, manifested as complete insensitivity to the GEX1A treatments. (Fig. 1f). These data indicate that the SGR4 may have different structural features³³ that affect its binding to the GEX1A splicing modulator, resulting in GEX1A-insensitive phenotypes in seed germination, primary root length, and lateral root growth.

GEX1A and other splicing modulators inhibit splicing leading to considerable levels of intron retention (IR) in plants³¹,³². We hypothesized that the GEX1A-tolerance phenotype was linked to splicing regulation. In our previous study, we have tested the SGR4 exhibited GEX1A-resistant molecular phenotypes on IR and splicing inhibition³³. To show that this effect is wider and applicable to a large set of genes, we selected another group of rice genes known to undergo AS and used RT-PCR to determine the level of IR under GEX1A treatment. Specifically, RT-PCR primer sets were designed

for genes encoding the AP2-like ethylene-responsive transcription factor AINTEGUMENTA (LOC_Os04g55970), oryzain beta chain precursor (LOC_Os04g57440), phosphoenolpyruvate carboxykinase (LOC_Os03g15050), and phospho-2-dehydro-3-deoxyheptonate aldolase, chloroplast precursor (LOC_Os07g42960). We treated 1-week-old seedlings of WT and the SGR4 mutant with 0.3 µM GEX1A for 6 h. Subsequently, using cDNA from each line, we performed RT-PCR targeting exons encompassing the intron involved in an AS event. Our RT-PCR data corroborate our phenotypic data showing that AS in SGR4 is insensitive to GEX1A treatment and remains similarly efficient as the AS observed under control conditions suggesting that the splicing-dependent gene expression pattern induced by the drug can be overcome by the SGR4 mutation (Fig. 1g).

These data link the GEX1A-tolerance phenotype to splicing and indicate that the structure of the SF3B1 variants and their binding to the GEX1A splicing modulator determine the resulting splicing efficiency and response to splicing inhibition phenotypes. Therefore, the SGR4 mutant provides a valuable resource to investigate the molecular mechanisms underlying splicing inhibition by splicing modulators and its impact on plant growth phenotypes.

**SGR4 inhibits GEX1A-induced change in global gene expression and splicing.** To determine the effect of the SF3B1
<figure><img src="image_4.png" /><figcaption>**Fig. 2 SGR4 protects against GEX1A-induced gene expression and splicing inhibition.** a Overlap of differentially expressed genes (DEGs) between WT and SGR4 under mock or GEX1A treatment. One-week-old rice seedlings were treated with 0.3 µM GEX1A for 6 h. A high number of DEGs were observed after GEX1A treatment in WT while SGR4 appears insensitive to GEX1A. SGR4 comparison with WT under mock treatment show very few number of DEGs. b GO enrichment analysis of the DEGs regulated by GEX1A. c Venn diagram showing the overlap of differentially spliced events (DAS) between WT and SGR4. d GO enrichment of GEX1A-responsive DAS genes. e Overlap of DEG with DAS events for SGR4 and WT under GEX1A treatment.</figcaption></figure>

<table><caption>**Table 1 Number of DEG per comparison.**</caption><thead><tr><th>Comparison</th><th>Direction</th><th>DEG number</th></tr></thead><tbody><tr><td rowspan="2">SGR4X vs SGR4C</td><td>Up</td><td>17</td></tr><tr><td>Down</td><td>4</td></tr><tr><td rowspan="2">WTC vs SGR4C</td><td>Up</td><td>21</td></tr><tr><td>Down</td><td>47</td></tr><tr><td rowspan="2">WTX vs SGR4X</td><td>Up</td><td>1115</td></tr><tr><td>Down</td><td>1127</td></tr><tr><td rowspan="2">WTX vs WTC</td><td>Up</td><td>1041</td></tr><tr><td>Down</td><td>859</td></tr><tr><td rowspan="2">OGR1X vs OGR1C</td><td>Up</td><td>207</td></tr><tr><td>Down</td><td>217</td></tr><tr><td rowspan="2">WTC vs OGR1C</td><td>Up</td><td>306</td></tr><tr><td>Down</td><td>438</td></tr><tr><td rowspan="2">WTX vs OGR1X</td><td>Up</td><td>1272</td></tr><tr><td>Down</td><td>663</td></tr><tr><td rowspan="2">WTX vs WTC</td><td>Up</td><td>1345</td></tr><tr><td>Down</td><td>736</td></tr></tbody></table>

mutations in SGR4 on genome-wide gene expression and pre-mRNA splicing, we performed RNA-seq analysis of 1-week-old rice seedlings of WT and SGR4 using polyA+ RNA. We identified 68 differentially expressed genes (DEGs) under control conditions (Fig. 2a, Table 1).

To further dissect the effect of SGR4 on AS, we analyzed differential splicing between these lines using rMATS. Similarly, to differential expression analysis, SGR4 had a minor impact on all types of splicing events in control conditions (Fig. 2c, Table 2). To investigate the genome-wide effects of SGR4 and the resulting molecular responses to splicing-inhibition, 1-week-old seedlings of WT and SGR4 were treated with 0.3 µM GEX1A for 6 h and RNA deep sequencing was performed to do a similar analysis on differential gene expression and AS. Indeed, our previous studies have shown that GEX1A treatment results in significant inhibition of constitutive splicing and AS³¹,³². Splicing inhibition is likely due to the binding of GEX1A to the SF3B1 proteins and the perturbation of the U2 snRNP machinery. We recently showed

that GEX1A may not be capable of binding to the SF3B1 variants in the SGR mutants or that this binding may be compromised, which would thus compromise GEX1A-mediated splicing inhibition³³. We wanted to determine whether the SF3B1 mutant variant maintains its splicing functions or whether they result in global differential splicing due to its potentially compromised binding to the splicing modulator. As expected, when treated with 0.3 µM GEX1A, WT seedlings exhibit major changes in DEG number (Table 1) and in the retention of a large number (2832) of introns but also on all types of AS events (3' SS, 5' SS, and ES) (Table 2). Nevertheless, in line with the above-described phenotypic data, our analysis showed that gene expression was almost unaffected by GEX1A treatment in SGR4. After GEX1A treatment, we detected only 21 DEG (Table 1) and 320 AS events (Table 2) in SGR4, further corroborating our data that SGR4 is largely insensitive to GEX1A. To identify genes whose GEX1A response is affected at DEG and splicing level in WT in comparison to SGR4, we compared at genome-wide level SGR4 and WT under GEX1A treatment. We identified a large number of DEG and DAS events, which show strong overlap with GEX1A-induced DEG and DAS between WT and SGR4 respectively (Fig. 2a, c). Consistent with our previous work showing that GEX1A preferentially affected splicing of stress-related genes³¹,³², GO analysis of GEX1A-regulated DEG also indicated enrichment in transport, secondary metabolic process, response to stress, response to endogenous stimulus, response to chemical, response to biotic stimulus, pollination, and catabolic process (Fig. 2b). At the splicing level, DAS genes were enriched in multiple cellular pathways like cellular protein modification process, response to endogenous stimulus, lipid metabolic process, and signal transduction (Fig. 2d). Although the global GO is in the same pathways, we have observed very little correlation between DEG and DAS events (Fig. 2e), which indicate most of the DEG are not alternatively spliced and vice versa.

Overall, our data confirm that global GEX1A-induced DEG and splicing inhibition is apparent in WT whereas negligible effects on DEGs and DAS are observed in SGR4 following GEX1A treatment, which supports its highly GEX1A-insensitive phenotype.
<table><caption>Table 2 Number and type of AS events per comparison.</caption><thead><tr><th>Comparison</th><th>Type of event</th><th>Number of AS events</th></tr></thead><tbody><tr><td rowspan="4">SGR4C vs SGR4X</td><td>A3S</td><td>43</td></tr><tr><td>A5S</td><td>62</td></tr><tr><td>IR</td><td>175</td></tr><tr><td>ES</td><td>40</td></tr><tr><td rowspan="4">WTC vs SGR4C</td><td>A3S</td><td>28</td></tr><tr><td>A5S</td><td>12</td></tr><tr><td>IR</td><td>96</td></tr><tr><td>ES</td><td>15</td></tr><tr><td rowspan="4">WTX vs SGR4X</td><td>A3S</td><td>141</td></tr><tr><td>A5S</td><td>199</td></tr><tr><td>IR</td><td>2603</td></tr><tr><td>ES</td><td>143</td></tr><tr><td rowspan="4">WTX vs WTC</td><td>A3S</td><td>181</td></tr><tr><td>A5S</td><td>228</td></tr><tr><td>IR</td><td>2855</td></tr><tr><td>ES</td><td>173</td></tr><tr><td rowspan="4">OGRX vs OGRC</td><td>A3S</td><td>54</td></tr><tr><td>A5S</td><td>69</td></tr><tr><td>IR</td><td>469</td></tr><tr><td>ES</td><td>45</td></tr><tr><td rowspan="4">WTC vs OGRC</td><td>A3S</td><td>152</td></tr><tr><td>A5S</td><td>94</td></tr><tr><td>IR</td><td>1452</td></tr><tr><td>ES</td><td>82</td></tr><tr><td rowspan="4">WTX vs OGRX</td><td>A3S</td><td>170</td></tr><tr><td>A5S</td><td>159</td></tr><tr><td>IR</td><td>2828</td></tr><tr><td>ES</td><td>112</td></tr><tr><td rowspan="4">WTX vs WTC</td><td>A3S</td><td>272</td></tr><tr><td>A5S</td><td>300</td></tr><tr><td>IR</td><td>3841</td></tr><tr><td>ES</td><td>227</td></tr></tbody></table>

Expression of engineered PHF5A-Y36C, another member of the branch point recognition complex, confers tolerance to GEX1A. The SF3B complex consists of seven proteins with distinct domain structures and molecular sizes (Fig. 3a)³⁴, and none of these proteins is characterized in plants yet. The SF3B1 together with PHF5A recognize and bind the BP during splicing reaction. The SF3B1 HEAT repeats grab the U2/branch point sequence (BPS) duplex and BP Adenosine is bulged out and covered in a pocket between HEAT repeats and PHF5A³⁵. This structure permits the SF3b complex to protect the BP from premature nucleophilic attack of its 2'-OH and adopt a proper confirmation before the first transesterification reaction¹³. Splicing inhibitors like GEX1A disrupt the interactions of the SF3B1 and PHF5A with BP and impairs the pre-mRNA splicing³⁶. Recently, in an effort to identify clones tolerant to splicing modulators such as E707, whole-exome sequencing data identified a PHF5A mutant, specifically PHF5A-Y36C, as a molecular target in mammalian cells conferring tolerance to the E707 splicing modulator¹⁸,²².

To identify the homolog of human PHF5A in rice, we interrogated the rice genome database and identified two homologous OsPHF5A genes (*LOC_Os04g56760* and *LOC_Os05g30410*), which show little variation in DNA sequence and encode identical proteins (Fig. 3b). We named these loci OsPHF5A1 (*LOC_Os04g56760*) and OsPHF5A2 (*LOC_Os05g30410*). A protein alignment showed that the majority of protein regions share high sequence similarity between human and rice PHF5A homologs (Fig. 3b). In particular, the Y36 residue is conserved among these protein sequences. We used the OsPHF5A1 (*LOC_Os04g56760*) DNA

sequence and introduced a single nucleotide substitution mutation (TAT to TGT, Y36C) to create PHF5A-Y36C for expression in rice (Fig. 3c). The synthetic PHF5A-Y36C cassette was FLAG-tagged and expressed under the control of the OsUBIQUITIN promoter. In the T₀ plants, PHF5A-Y36C expression was analyzed using FLAG-tag, and confirmed lines were used for further experiments (Fig. 3d).

In the progeny, we examined the effects of the GEX1A splicing modulator on the growth and development of PHF5A-Y36C-OX plants. Our phenotypic analysis indicated that these PHF5A-Y36C-OX plants were partially tolerant to GEX1A- and PB-treatment (Fig. 3e, Supplementary Fig. 1), which is different from SGR4. We termed this line OGR (Overexpression-PHF5A-Y36C GEX1A Resistance). The overexpression of PHF5A-Y36C had no apparent effect on plant growth and development and OGR plants were phenotypically indistinguishable from WT plants.

To further dissect the role of OsPHF5A in rice growth and splicing, we employed CRISPR/Cas9-targeted mutagenesis to engineer functional knockout of the OsPHF5A genes. As we did not find a PAM site to cut exactly at Y36, we designed two sgRNAs closest to the available PAM site to target each of the OsPHF5A loci (Supplementary Fig. 2A, B). In the T₀ generation, we recovered several single knockout mutants. For OsPHF5A1 (*LOC_Os04g56760*), targeting via sgRNA-1 produced a homozygous monoallelic knockout mutant with two nucleotide deletions that caused a stop codon after K25 (Supplementary Fig. 2A). Similarly, targeting via sgRNA-2 also produced a homozygous monoallelic knockout mutant with one nucleotide insertion that produced a frameshift after C40 and a stop codon (Supplementary Fig. 2A). To target OsPHF5A2 (*LOC_Os05g30410*), we also used two sgRNAs. Targeting via sgRNA-1 produced a homozygous monoallelic knockout mutant with one nucleotide insertion that caused a frameshift after Y36 and a stop codon (Supplementary Fig. 2B). Using sgRNA-2, we produced a homozygous biallelic knockout mutant. One allele has a single nucleotide deletion and caused a frameshift after Y36 with 38 additional residues. The second allele has one nucleotide insertion and caused a frameshift after V37 and a stop codon (Supplementary Fig. 2B).

We tested the effects of GEX1A on the growth and development of OsPHF5A homozygous mutant lines. We treated 3-day-old seedlings with 0.3 µM GEX1A and observed subsequent primary root growth. Our data show that these knockout mutants are partially tolerant to GEX1A treatment (Supplementary Fig. 2C). We termed these OsPHF5A mutants PGR (PHF5A GEX1A Resistance). The degree of GEX1A tolerance was similar in all PGR mutants but slightly weaker than that in OGR seedlings (Supplementary Fig. 2C).

Overall, these results indicate that the expression of mutant variant PHF5A-Y36C at a single locus is sufficient to perturb GEX1A binding. Similarly, the two rice PHF5A genes may have additive effects and knocking out either of them perturbed the binding with splicing inhibitors. However, it seems that OsPHF5A is essential for plant survival as we did not recover a functional double knockout mutant.

**Expression of PHF5A-Y36C represses GEX1A-mediated splicing inhibition.** To determine the role of PHF5A in genome-wide pre-mRNA splicing under splicing inhibition, we treated 7-day-old seedlings of WT and OGR with 0.3 µM GEX1A and performed RNA-seq analysis using polyA+ RNA. First, we analyzed gene expression in OGR as compared to WT in control condition. We identified 744 DEGs in OGR (Table 1) and 1780 AS events in OGR as compared to WT in control condition, among which were 152 3' SS, 94 5' SS, 82 ES, and 1452 IR genes (Table 2). The comparison of
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Generation of PHF5A-Y36C heterologous expression in rice. **a** The structure of SF3B complex binding to U2 and the pre-mRNA intron. The branch point is in contact with SF3B1 and PHF5A. The branch point recognition sequence (BPRS) of U2 snRNA binds to branch point and flanking nucleotides. **b** GEX1A interrupts these interactions and inhibit this early spliceosome assembly. The human PHF5A protein sequence is aligned with sequences of rice PHF5A homologs. PHF5A-Y36C (highlighted in bold) confers tolerance to splicing modulators (Teng et al., 2017) and is highly conserved among eukaryotes. **c** The coding sequence of rice OsPHF5A (LOC_Os4g56760) was mutated (Y36C; TAT to TGT) and cloned into a construct under the control of the Ubiquitin promoter with a 3XFLAG-tag. **d** Confirmation of PHF5A-Y36C heterologous expression in transgenic lines. Total protein was extracted from independent lines of *O. sativa* cv *japonica* transformed with UBIp::PHF5A-Y36C. Anti-FLAG antibody was used to detect FLAG-tagged PHF5A-Y36C. The arrow indicates the presence of PHF5A. WT was used as a negative control. **e, f** Heterologous expression of PHF5A-Y36C confers tolerance to GEX1A. Rice seeds were germinated on 1/2 MS media for 3 days, then transferred onto 1/2 MS media supplemented with 0.3 µM GEX1A for 3 days. Root tips were marked to observe post-transfer growth. The transgenic lines were termed OGR (Overexpression-PHF5A GEX1A Resistance).</figcaption></figure>

WT and OGR under GEX1A treatment show a high number of DEG (1935) and most of these overlap with the ones identified as GEX1A responsive in WT (Fig. 4a). But the treatment with GEX1A of OGR led to only 424 DEG as compared to mock conditions showing that OGR is less sensitive to GEX1A, in agreement with the phenotyping data (Table 1).

AS analysis showed that OGR exhibited a much lower number of differential splicing events compared to WT following GEX1A treatment, confirming that the Y36C mutation provides tolerance to splicing inhibition by GEX1A (Fig. 4c). The numbers of DAS events were 4640 in WT and only 637 in OGR (Table 2). A large number of DAS events were identified when comparing WT and OGR after GEX1A splicing inhibition. Most of these DAS events overlap with GEX1A-responsive events in WT (Fig. 4c). IR events were the most abundant type of AS event affected by OGR under GEX1A but all types of AS events were also represented. As previously shown for SGR4, the GO analysis shows that most of the DEG identified during comparison of WT and OGR under GEX1A treatment are enriched in response to chemical, response to stress, secondary metabolic process, response to endogenous stimulus and response to biotic stress (Fig. 4b). Whereas, DAS were enriched in genes involved in various functional categories such as cellular protein modification process, transport, carbohydrate and lipid metabolic processes and response to an endogenous stimulus (Fig. 4d). We did not observe a very high correlation between DEG and DAS events and most of these DEG and DAS are unique compared to those under GEX1A treatment conditions (Fig. 4e).

Taken together our data show that splicing-inhibition treatment affects genome-wide gene expression in WT. However, global GEX1A-induced splicing inhibition is clearly evident in WT whereas minimal effects on DEGs and DAS are observed in OGR following GEX1A treatment, which supports its highly GEX1A-tolerant phenotype.

**SF3B1 and PHF5A-regulated genes and AS events showed significant overlap.** SF3B1 and PHF5A are partner proteins and function in the SF3B complex as part of the U2 snRNP. The structural and functional affinity emphasizes a common downstream effect on differentially expressed and differentially spliced genes. To identify the similar and differential effect of these partner proteins, we compared the data sets of SGR4, carrying an SF3B1 mutant variant insensitive to GEX1A, with OGR, a PHF5A-Y36C overexpression line highly tolerant to GEX1A. We compared the DEGs and DAS genes significantly affected in GEX1A-treated SGR4 or OGR with those of their respective WT (Fig. 5). This identified the genes that do not respond to GEX1A inhibition in these mutants but that are targeted by GEX1A in WT plants. The data show that OGR and SGR4 are impaired in the GEX1A response involving 822 DEGs and 1134 DAS events (Fig. 5a, b). Heatmap and hierarchical clustering of DEG and splicing of GEX1A-sensitive genes in WT further revealed that SGR4 and OGR similarly affect GEX1A-dependent splicing and DEG in a way that is opposite to what is observed in WT (Fig. 5c, d). However, there is a large number of DAS events that
<figure><img src="image_6.png" /><figcaption>**Fig. 4 OGR impedes the effects of GEX1A on global gene expression and splicing patterns.** a Overlap of differentially expressed genes (DEGs) between WT and OGR under mock or GEX1A treatment. The number of DEGs much higher in WT compared to OGR. Heterologous expression of PHF5A-Y36C affects the global gene expression under mock conditions. b GO enrichment analysis of the DEGs regulated by GEX1A. c Overlap of differentially spliced events (DAS) between WT and OGR. d GO enrichment of GEX1A-sensitive DAS genes. e Overlap of DEG with DAS events for OGR and WT under GEX1A treatment.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 5 SF3B1- and PHF5A-regulated genes and AS events displayed significant overlap.** Overlap of DEGs a and DAS events b of GEX1A-responsive gene between SGR4 and OGR. Heatmap and hierarchical clustering of log2FC of DEG. Heatmap and hierarchical clustering of log2FC of DEG c and delta Ψ (change in splicing) of each intron-retention event d significantly affected by GEX1A in WT_sfb1, WT_phf5, SGR4, and PHF5 overexpressing line (OGR). WT of the SF3B1 (WT_sfb1) and PHF5A (WT_phf5).</figcaption></figure>
<figure><img src="image_8.png" /></figure>

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

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

<figure><img src="image_11.png" /><figcaption>**Fig. 6 SGR4 and OGR show specificity for splicing-related features.** a Intron GC content, b Intron length, c Polypyrimidine tract length, and d Polypyrimidine tract score (as defined by SVM-BP, http://regulatorygenomics.upf.edu/Software/SVM_BP/) of introns not affected by GEX1A in SGR4 (SGR4 nondiff) or OGR (PHF5A-Y36C heterologous expressor, OGR_nondiff), introns differentially retained in GEX1A-treated SGR4 or OGR as compared to their corresponding WT (WTX-vs.SGR4X and WTX vs OGRX, respectively) in the same condition, and introns affected by GEX1A in the WT of the SGR4 or OGR experiment (WTC-vs.WTX (sfb1) and WTC-vs.WTX (phf5), respectively). The polypyrimidine tract score indicates C/T content of the polypyrimidine tract between Branch Point and 3' splice site.</figcaption></figure>

do not overlap (Fig. 5b) which might point toward independent roles during splicing of specific mRNA targets.

IR results from a pervasive splicing event induced by splicing inhibitors in plants. The AS events are known to be linked with specific splicing elements and nucleotide content, as well as with specific chromatin marks. In particular, intron GC content and the polypyrimidine tract downstream from the BP are important for splicing and may have evolved as recognition signals for the splicing machinery³⁷,³⁸. Therefore, we asked whether these sequence elements might also affect splice site recognition in SGR4 or OGR under splicing-inhibition conditions. We compared the sequence elements of retained introns (induced by drug treatment) in WT and SGR4/ OGR with those of nonaffected introns in SGR4/OGR. The retained introns have significantly enriched GC content (Fig. 6a). Similarly, the retained introns are shorter in length and displayed significantly weaker Py-tract length and score (Fig. 6b-d). Collectively, these data suggest an important role of these intronic features during GEX1A-mediated splicing inhibition.

**SGR4 and OGR show differing levels of sensitivity to combined salt stress and splicing inhibition.** Several reports have linked splicing mutations to sensitivity to abiotic stress factors including salt and drought stress³⁹⁻⁴⁵. To address the role of SF3B1- and PHF5A-dependent splicing events in relation to abiotic stress, we tested whether SGR4 and OGR had tolerant or sensitive responses to other abiotic stresses, as shown previously in Arabidopsis³¹,³². As mentioned before, GO analysis of GEX1A-induced DAS genes and DEGs revealed that GEX1A targets many abiotic stress-related genes. More specifically, 822 DEG in both SGR4 and OGR lines (Fig. 5a) were significantly enriched for transcription factors GO category including key regulatory players of the salt-stress

response in rice such as DREB2A, NAC5, and ERF922⁴⁵⁻⁴⁸. The 1134 common AS events (Fig. 5b) were also enriched for important abiotic stress responses such as members of the trehalose biosynthesis pathways which has been linked to salt-stress tolerance in rice⁴⁹.

Together, this suggested that, the splicing-inhibition treatment mimics an abiotic stress molecular response. We speculated that coupling splicing inhibition with salt stress would induce synergistic effects in stress responses. To test this hypothesis, we subjected 3-day-old seedlings to a 3-day combined treatment of 50 mM NaCl and 0.15 µM GEX1A. Treatment with 50 mM NaCl or 0.15 µM GEX1A alone slightly inhibits root growth in WT. We also tested the primary root growth of WT, OGR, and SGR4 at 100 mM and 200 mM NaCl and did not find any significant differences among these genotypes (Supplementary Fig. 3). However, following combined splicing inhibition at 0.15 µM GEX1A and salt-stress induction at 50 mM NaCl, the primary root growth of WT seedlings was significantly inhibited (Fig. 7a, b). The SGR4 and OGR showed variable levels of tolerance to combined splicing and salt-stress treatment, with OGR exhibiting slight root growth inhibition compared to SGR4, supporting a salt-stress dependent phenotype depending on splicing levels (Fig. 7a, b).

To further validate these responses, we treated the 7-day-old WT seedlings with 200 mM NaCl for 6 h and conducted genome-wide RNA sequencing to identify salt-responsive DEGs and DAS events. We compared the GEX1A-insensitive DEGs and DAS events in SGR4 and OGR to DEG and DAS events in response to salt (Fig. 7c, d). Our analysis showed that 318 out of 600 (53%) of salt-responsive DEGs and 261 out of 715 (36.5%) of salt-responsive AS events overlapped with GEX1A-responsive genes in SGR4 and OGR (Fig. 7c, d).
<figure><img src="image_12.png" /><figcaption>**Fig. 7** Splicing inhibition sensitizes the WT seedlings against salt stress. a, b Seedlings of WT, SGR4, and OGR were germinated for 3 days on ½ MS media and then transferred to ½ MS media supplemented with different concentrations of NaCl and/or GEX1A for a further 3 days. Root tips were marked to observe post-transfer growth. On 50 mM NaCl or 0.15 µM GEX1A alone, seedling root growth was slightly affected, but severe root growth inhibition was observed in WT seedlings on 50 mM NaCl and 0.15 µM GEX1A combined (n = 3). Overlap of DEG c and DAS d in response to salt treatment (7-day-old WT seedlings treated with 200 mM NaCl for 6-h) in WT (WTC vs WTS) with GEX1A-resistant DEG and DAS genes in SGR4 (WTX vs SGR4X) and OGR (WTX vs OGRX).</figcaption></figure>

We have previously shown that the splicing inhibition via PB and GEX1A triggers the expression of stress-responsive reporters *RD29A::LUC* and *MAPKKK18::uidA* and activates the ABA signaling pathway in *Arabidopsis*<sup>31,32</sup>. To better understand how the activation of GEX1A regulates stress signaling in rice, we investigated the transcriptional regulation of stress-responsive genes under GEX1A treatment. One-week-old seedlings of OGR, SGR4, and WT were treated with 0.3 µM GEX1A, 200 mM NaCl, and Mock for 6 h. The expression levels of rice stress-responsive transcription factors (*OsASR4*, *OsbZIP*, *OsERF*, *OsNAC*, *OsWRKY*) and several other genes were tested (Fig. 8a, Supplementary Fig. 4). These analyses showed that the GEX1A activates the expression of the stress-responsive genes similar to salt-stress treatment. The transcript levels of *OsASR4*, *OsERF124*, *OsNAC4*, and *OsWRKY27* were significantly elevated, while the transcript level of *OsbZIP* was significantly reduced in WT after GEX1A treatment similar to a salt stress (Fig. 8a). The expression levels of these transcripts were not affected in the OGR and SGR4 after GEX1A treatment (Fig. 8a, Supplementary Fig. 4). Interestingly, the expression of all stress-responsive genes was altered in OGR, and SGR4 like WT under salt-stress treatment. These results show that GEX1A triggers the expression of stress response in rice, a response lost in OGR and SGR4 mutants after GEX1A treatment only because of the inability to bind the GEX1A. This is further verified by the splicing patterns analysis of stress-responsive genes in OGR and SGR4 after GEX1A treatment (Fig. 8b, Supplementary Fig. 5). We selected a set of genes whose splicing is inhibited by the GEX1A in WT rice plants. In our analysis, we did not observe splicing inhibition in SGR4 for any target gene corroborating the GEX1A-insensitive phenotype of SGR4 (Fig. 8b, Supplementary Fig. 5). In OGR the intron retention is reduced for *OsSRDW*, *OsWAK5* and *OsERF109* as compared to wild-type. We did not observe the splicing inhibition after salt treatment for these sets of genes.

Together, these data indicate that GEX1A inhibits the splicing of stress-responsive genes and the GEX1A splicing-inhibition overlay salt-stress responses, highlighting a link between splicing and abiotic stress responses in plants that is dependent on the SF3B-PHF5A complex. The generation of gain of function mutants allowed us to detect this overlay stress response linked to splicing regulation in plants.

## Discussion

**SGR4 showed insensitivity to splicing-inhibitor GEX1A.** Structural studies have revealed that the SF3B complex plays a role during BP recognition<sup>12,18,23,38,46</sup>. The interaction of SF3B1 HEAT domains and BP adenosine activates the first catalytic step of the splicing reaction<sup>47,48</sup>. These HEAT domain repeats are the main target of small molecules that grasp the BP and have anti-tumor activities via inhibition of splicing. There are many classes of these small molecules, termed splicing inhibitors/modulators, of which pladienolides, herboxidiene, and their analogues are well-studied examples<sup>19,49–51</sup>. A number of mutations in the HEAT repeats of SF3B1 have been reported that result in tolerance to splicing inhibitors<sup>18,21,22,52,53</sup>. In a previous study, we evolved a series of splicing-inhibitor tolerant SF3B1 versions in rice; these versions carry mutations in the HEAT repeats or the unstructured part of SF3B1<sup>33,54,55</sup>.

The SGR4 carries the substitutions in the HEAT repeat (K1049R, K1050E, and G1051H). The detailed phenotypic analysis of SGR4 carrying the *SF3B1* mutations at different developmental stages showed that these plants are insensitive to splicing inhibitors. These inhibitors affect germination and primary root growth. Lateral root growth is not only controlled by a different mechanism than primary root<sup>56–58</sup> but also requires a functional splicing machinery. We showed here that lateral roots as well as the other phenotypes induced by splicing inhibitors were not observed in SGR4 lines. These observations
<figure><img src="image_13.png" /><figcaption>**Fig. 8** Expression and splicing pattern analysis of stress-responsive genes after GEX1A and salt treatment in OGR and SGR4. One-week-old rice seedlings of OGR, SGR4, and WT are treated with mock, 0.3 µM GEX1A, and 200 mM NaCl for 6-h. Total RNA extracted from the whole seedling was used for mRNA expression and splicing pattern analysis. **a** Expression of stress-responsive genes in OGR, SGR4, and WT. The GEX1A triggers the expression patterns similar to salt-stress treatment for all genes in WT but not in SGR4 and OGR. Bars represent the mean ± SEM of three replicates. *OsActin* was used as an internal control. (Student's t-test; \*P < 0.05, \*\*P < 0.01). **b** Semiquantitative RT-PCR analysis of alternative splicing patterns of stress-responsive genes in OGR, SGR4, and WT. No intron retention is observed in SGR4 under GEX1A treatment. Arrowheads indicate splicing variants that changed following GEX1A treatment. The gene structures and retained introns are shown. Red boxes indicate the PCR fragments.</figcaption></figure>

are further supported by global transcriptomic and splicing analyses. Our analyses indicate that these SGR mutants are barely affected in transcription and splicing compared to the large impact observed by the treatment with splicing inhibitors to the WT.

**Expression of engineered PHF5A confers tolerance to splicing inhibitors.** PHF5A is another protein of the SF3B complex and is involved in BP interaction together with SF3B1. Studies have shown that splicing inhibitors bind and fit in the pocket of the PHF5A-SF3B1 complex and consequently disrupt the interaction with BP. A single amino acid substitution in PHF5A (Y36C) interrupts the binding with splicing inhibitors and results in tolerance to these splicing inhibitors<sup>18,22</sup>. Our analyses also indicate that the overexpression of PHF5A-Y36C confers tolerance to splicing inhibition in rice. Rice has two identical homologs of human PHF5A, which were probably evolutionarily conserved to help the plant survive environmental challenges. We used CRISPR/Cas9 to target both loci simultaneously to produce a complete *OsPHF5A* knockout mutant in rice. Of note, we also used two sgRNAs to target the two *OsPHF5A* loci simultaneously in an attempt to generate a complete knockout of PHF5A function in plants, but we did not recover a homozygous double knockout mutant for both *OsPHF5A* loci. Recent

structural analysis data indicate that PHF5A-Y36C forms a direct contact with BPA, thus implicating PHF5A function in BP recognition<sup>18,30</sup>. Such an essential role of PHF5A in splicing may explain our inability to recover a double mutant using CRISPR/Cas9-targeted mutagenesis; it appears that an *OsPHF5A* double mutant is embryonic lethal. Interestingly, we produced single knockout mutants for both of the *PHF5A* loci in rice. These single knockout mutants also showed a certain degree of tolerance to splicing inhibitors, albeit weaker than the one conferred by PHF5A-Y36C. We performed global gene expression and splicing analyses and found that under control conditions, the number of DEGs and DAS genes are significantly higher in PHF5A-Y36C-expressing plants compared to WT plants. This observation is in contrast to what is observed in human cell lines where expression of PHF5A-Y36C does not significantly affect splicing<sup>18</sup>.

In yeast, the spliceosome undergoes multiple ATP-dependent conformational changes. After activation of the spliceosome, Prp2 removes U2 component SF3a/b to allow the interaction of the BP and the 5' SS so that the catalytic reaction can take place. The Prp2 activity is ATP-dependent and it moves in a 3' to 5' direction toward the branch site to displace SF3a/b from the spliceosome. The dissociation of SF3a/b is necessary for the binding of Cwc25 to facilitate the branching reaction<sup>13,59–61</sup>.
Based on the yeast splicing mechanism, we speculate that in PHF5A-Y36C-expressing plants under normal conditions, a high amount of PHF5A is produced and this increases the likelihood that PHF5A remain bound to the U2 snRNP during the catalytic phase. This delays the release of SF3a/b from the spliceosome and causes partial disruption of pre-mRNA splicing. However, under GEX1A treatment, the WT PHF5A is not functional and only PHF5A-Y36C is able to pass through the activation and catalytic phase and subsequent splicing steps. Taken together, accurate control of PHF5A homeostasis is needed for precise pre-mRNA splicing reaction under normal or stress conditions.

**Regulation of gene expression and AS to fine-tune plant responses to abiotic stresses.** Plants as sessile organisms have evolved complex structures and systems to cope with environmental stresses and adapt their growth and development. Rapid perception and consequent reprogramming of gene expression play a significant role in the adaptation of plants to abiotic stress and environmental constraints. The splicing process is accomplished by hundreds of proteins including the splicing factors (SFs) serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs). These SFs bind to the regulatory *cis*-elements in the pre-mRNA thus act as activators and repressors of splice site selection⁶²,⁶³. However, the pre-mRNA of these SFs can also be alternatively spliced under different environmental conditions emphasizing the occurrence of feedbacks in splicing regulation. Transcription factors (TFs) are the principal regulators of plant growth and development and stress-responsive TFs like AREB, DREB, NAC, and 14-3-3 also undergo alternative splicing under abiotic stress conditions⁶⁴⁻⁶⁷. Understanding the molecular basis of this regulation will unlock the potential to harness the splicing machinery to answer basic biological questions and for biotechnological applications. Small-molecule splicing inhibitors capable of targeting different splicing factors provide invaluable tools to underpin the molecular regulation of the splicing machinery³¹,³²,⁶⁸.

Our work also established a link between splicing inhibition and responses to abiotic stress. Based on our observations, we propose a model linking modulation of splicing factors action with adaptive stress responses (Fig. 9). When a plant submits an environmental stress there are several responses at transcriptional level but also at post-transcriptional levels such as AS. In mutants where splicing sensitivity to a particular drug is affected, the effect of stress is different than the one observed in wild-type pointing to an overlay regulatory level of AS acting during stress adaptation. Specific isoforms of transcription factors may have different roles on gene regulation and may also control the expression of splicing factors necessary for the induction of those isoform changes revealing a complex interaction between splicing regulation and TF action. Hence, TFs affect the expression of splicing factors and in turn, these splicing factors regulate their splicing patterns and indirectly their potential function. These regulatory networks involving TFs, SFs and their targets may help plants to cope with abiotic stresses (Fig. 9). From stress perception to adaptive stress response, the pathway is intertwined at many steps. This is further validated by previous work that indicates AS has regulatory functions during plant development⁶,⁵⁷,⁶⁹⁻⁷¹.

Our study sets a milestone for the use of splicing modulators to fine-tune the splicing mechanism to cope under abiotic stress conditions. Finally, there is a pressing need to develop cheaper analogues of splicing inhibitors to use as herbicides and our findings can have potential applications for generating crops with herbicide tolerance.

**Methods**

**Plant materials, chemicals, and vector construction.** *Oryza sativa* L. ssp. *japonica* cv Nipponbare was used for all experiments. The SGR4 mutant is described elsewhere³³.

GEX1A/ Herboxidiene (CAS: 142861-00-5) was purchased from BOC Sciences (45-16 Ramsey Road, Shirley, NY 11967, USA).

To target *OsPHF5A* (*LOC_Os04g56760*) and *OsPHF5A-Like* (*LOC_Os05g30410*), sgRNAs were cloned in the *pRGEB32* plasmid⁷². The expression of Cas9 was driven by *OsUbiquitin* and the sgRNA was expressed under the *OsU3* promoter. The *pRGEB32* plasmid was digested with *Bsal* and sgRNAs were synthesized as oligonucleotides with *Bsal* overhangs, **GGCA** in the forward and **AAC** in the reverse oligonucleotides. The oligonucleotides were annealed and ligated in the *Bsal*-digested vector. The oligonucleotide sequences are listed in Supplementary Table 1.

To generate the PHF5A-Y36C heterologous expression line, the synthetic fragment of rice *OsPHF5A* (*LOC_Os04g56760*) with Y36C modification and 3xFLAG-tag at N-terminus was cloned under the control of the *OsUbiquitin* promoter via a Gateway cloning strategy. The sequence of the synthetic fragment is given in Supplementary Table 1.

**Rice transformation and genotyping of transgenic plants.** Rice transformation was performed using the agrobacterium strain EHA105 and genotyped as described previously⁵⁶. Briefly, the sgRNA-targeted region was amplified by PCR using locus-specific primers. Purified PCR products were cloned using CloneJET PCR Cloning Kit (K1231). Sanger sequencing was conducted to analyze the mutation.

**Phenotypic analysis.** The germination inhibition assays were performed in six-well culture plates. Rice seeds were surface-sterilized and germinated in 4 mL of H₂O for controls and with 5 µM and 10 µM of GEX1A. The germination rate was calculated after 5 d.

For the primary root elongation inhibition assay, sterilized rice seeds were germinated on ½ MS agar media (without sucrose) square plates for three days in a vertical position. On the third day, seedlings with similar root lengths were transferred to ½ MS media (without sucrose) supplemented with different concentrations of GEX1A and/or NaCl and grown in a vertical position. The tip of the root was marked and subsequent root growth was calculated after 3 days.

The elongated root length after transfer to control plates was set at 1 (100%). The root elongation rates on the chemical-containing plates were calculated as (elongated root length on chemical)/ (elongated root length on control) × 100%. Values are means ± SE (n = 3). Significance (*P* < 0.05) was assessed by the Student's *t*-test.

For lateral root inhibition analysis, rice seeds were germinated on ½ MS agar media (without sucrose) square plates for three days. Seedlings with similar root size were transferred to ½ MS media (without sucrose) supplemented with 0.5 µM and 1 µM GEX1A plates. Lateral root growth was observed 3 days after transfer using a Nikon SMZ25 stereomicroscope.

**Immunoblot analysis.** Total proteins were extracted from 100 mg of sample using extraction buffer (100 mM Tris-Cl pH8, 150 mM NaCl, 0.6% IGEPAL, 1 mM EDTA, 3 mM DTT with protease inhibitors, PMSF, leupeptin, aprotinin, pepstatin, antipain, chymostatin, Na₂VO₃, NaF, MG132, and MG115. Proteins were separated on a 10% polyacrylamide gel. Immunoblot analysis was carried out using mouse α-FLAG M2 (dilution, 1:1000) antibody for FLAG-PHF5A from OGR transgenic rice plants. The antigens were detected by chemiluminescence using an ECL-detecting reagent (Thermo Scientific).

**RNA isolation and RNA-seq.** Total RNA was extracted from 7-day-old rice seedlings after 6 h of treatment with DMSO, 0.3 µM GEX1A or 200 mM NaCl using the Direct-zol RNA MiniPrep Plus (Zymo Research) according to the manufacturer's instructions. RNA was quantified using a Nanodrop and RNA quality was examined using a 2100 Bioanalyzer (Agilent Technologies). High-quality RNA samples with RNA integrity number (RIN) ≥ 7.0 were selected for library construction. The RNA-seq libraries were constructed using TruSeq mRNA stranded kit following the standard protocol and sequenced on the HiSeq-4000 platform to generate high-quality paired-end reads.

**RT-PCR and RT-qPCR.** For reverse-transcription PCR (RT-PCR), DNA digestion of total RNA samples was performed using an RNase-Free DNase Set (Invitrogen cat. No. 18068-015) following the manufacturer's protocol. The total RNA was reverse transcribed using a SuperScript First-Strand Synthesis System (Invitrogen) to generate cDNA. PCR conditions were: initial denaturation at 95 °C for 2 min, then 40 cycles of 95 °C for 30 sec, 55 °C for 30 sec, and 72 °C for 60 sec, then final elongation at 72 °C for 5 min. Primers used for RT-PCR are listed in Supplementary Table 1.

The qPCR was performed using the 50 ng RNA in 10 µl final reaction volume. The iTaq™ Universal SYBR® Green One-Step Kit (Bio-Rad, cat. no. 172-5150) was used with the following manufacturer's protocol conditions: reverse transcription 50 °C for 10 min, followed by 40 cycles of 95 °C for 15 s, 60 °C for 1 min, melt-curve analysis 65-95 °C, 0.5 °C increment. Primers used for qPCR are listed in Supplementary Table 2.
<figure><img src="image_14.png" /></figure>

**Fig. 9 Model for splicing-mediated adaptive stress responses.** Abiotic stress affects global gene expression by reprogramming it at the transcriptional and post-transcriptional levels. Plants perceive stress signals through multiple integrated signaling pathways, which regulate transcription and splicing processes. Transcription factors affect the expression of splicing factors, and in turn splicing factors regulate the splicing patterns of transcription factors. Splicing inhibitors (e.g., GEX1A) inhibit the splicing and activate the stress signaling in plants. The splicing factor protein variants insensitive to GEX1A (e.g., SGR4) impede the inhibitors' effects of splicing inhibitors. Stress-modulated transcripts and splice isoforms help plants to cope with abiotic stresses. Arrows indicate activation or repression, and T-bars indicate repression. CS constitutive splicing, ES exon skipping, 5' SS 5' splice site, 3' SS 3' splice site, IR intron retention.

**Analysis of RNA-seq data and gene functional classification.** RNA-seq pre-processing included trimming library adapters and quality controls with Trimmomatic using the following arguments; TrimmomaticPE - LEADING:25 TRAILING:25 CROP:120 MINLEN:120. Illumina adapters were removed and read pairs were trimmed to the same length. Processed reads were aligned using STAR with the following arguments: -outSAMtype BAM SortedByCoordinate-outReadsUnmapped Fastx-readFilesCommand zcat,-quantMode GeneCounts,-outFilterMultimapNmax 1. Reads overlappings exons per genes were counted using the FeatureCounts function of the Rssubreads package using the previously published GTF annotation files⁵⁷.

Significance of differential gene expression was estimated using DEseq2, and the FDR correction of the *p*-value was used during pairwise comparison between genotypes. A gene was declared differentially expressed if its adjusted *p*-value (FDR) was ≤ 0.01. Gene ontology enrichment analysis was done as in ref. ⁵⁸ and simplified using GO slim terms.

**Analysis of differential splicing.** Differential splicing was determined using rMATS v4.02 (http://rnaseq-mats.sourceforge.net/) using the following arguments -nthread 4 -readLength 120 -t paired -libType fr-firststrand. The gene annotation GTF file from⁵⁷ was as a reference annotation. Events with FDR < 0.001 and | IncLevelDifference| >0.2 were defined as differentially alternatively spliced²¹.

**Clustering analysis.** Comparison between sets of DEGs or DAS genes and events were done using the UpSetR and the eulerr R packages. Clustering was done using the hierarchical clustering function of the PheatMap package using euclidean distance. Clusters were extracted from the heatmap using the cutree function. From DEG, the log₂ fold change (log₂FC) values were used to cluster genes by their changes of expression. For splicing the IncLevelDifference of differentially retained introns was used to cluster splicing change profiles in the different genotypes.

**Intron and splice site analysis.** The length and the GC content of sequences of differentially retained introns in SGR4 and PHF5 mutants were analyzed using custom scripts. SVM-BPfinder (regulatorygenomics.upf.edu/Software/SVM_BP/) was used for scoring BP and Py-tracts.

**Reporting summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

**Data availability**

(1) The Raw sequencing data are available at SRA under the BioProject accession PRJNA636200. (2) The processed data (https://www.bam.fwd.bigWig) files are deposited in GEO (https://www.ncbi.nlm.nih.gov/geo/) under accession number GSE153234. The following secure token has been created to allow review of record GSE153234 while it remains in private status: egzpgmkkdfyhzuv. Please note the following points: (i) This token allows anonymous, read-only access to GSE153234, and associated accessions while they are private. (ii) Treat the token as you would a password and realize that the token provides access to GSE153234 to anyone who uses it. (3) The lists of DEGs and DAS events were attached with the manuscript as supplementary data files.

Received: 28 October 2020; Accepted: 26 March 2021; Published online: 05 May 2021

**References**

1. Will C. L. & Lührmann R. Spliceosome structure and function. *Cold Spring Harb Perspect. Biol.* **3**, a003707 (2011).

2. Plaschka, C., Pei-Chun, L. & Nagai, K. Structure of a pre-catalytic spliceosome. *Nature* **546**, 617 (2017).

3. MacMillan, A. M. et al. Dynamic association of proteins with the pre-mRNA branch region. *Genes Dev.* **8**, 3008–3020 (1994).

4. Wahl, M. C., Will, C. L. & Luhrmann, R. The spliceosome: design principles of a dynamic RNP machine. *Cell* **136**, 701–718 (2009).

5. Shi, Y. Mechanistic insights into precursor messenger RNA splicing by the spliceosome. *Nat. Rev. Mol. Cell Biol.* **18**, 655 (2017).

6. Staiger, D. & Brown, J. W. Alternative splicing at the intersection of biological timing, development, and stress responses. *Plant Cell* **25**, 3640–3656 (2013).

7. DeBoever, C. et al. Transcriptome sequencing reveals potential mechanism of cryptic 3' splice site selection in SF3B1-mutated cancers. *PLoS Comput. Biol.* **11**, e1004105 (2015).
8. Cretu, C. et al. Molecular architecture of SF3b and structural consequences of its cancer-related mutations. *Mol. Cell* **64**, 307–319 (2016).

9. Lee, B. M. & Christopher, L. A genomic view of alternative splicing. *Nat. Genet.* **30**, 13 (2002).

10. Graveley, T. W. N. & Brenton, R. Expansion of the eukaryotic proteome by alternative splicing. *Nature* **463**, 457 (2010).

11. Filichkin, S., Priest, H. D., Megraw, M. & Mockler, T. C. Alternative splicing in plants: directing traffic at the crossroads of adaptation and environmental stress. *Curr. Opin. Plant Biol.* **24**, 125–135 (2015).

12. Golas, M. M., Sander, B., Will, C. L., Luhrmann, R. & Stark, H. Molecular architecture of the multiprotein splicing factor SF3b. *Science* **300**, 980–984 (2003).

13. Lardelli, R. M., Thompson, J. X., Yates, J. R. 3rd & Stevens, S. W. Release of SF3 from the intron branchpoint activates the first step of pre-mRNA splicing. *RNA* **16**, 516–528 (2010).

14. Hasegawa, M. et al. Identification of SAP155 as the target of GEX1A (Herboxidiene), an antitumor natural product. *ACS Chem. Biol.* **6**, 229–233 (2011).

15. Kaida, D. et al. Spliceostatin A targets SF3b and inhibits both splicing and nuclear retention of pre-mRNA. *Nat. Chem. Biol.* **3**, 576–583 (2007).

16. Malcovati, L. et al. Clinical significance of SF3B1 mutations in myelodysplastic syndromes and myelodysplastic/myeloproliferative neoplasms. *Blood* **118**, 6239–6246 (2011).

17. Kotake, Y. et al. Splicing factor SF3b as a target of the antitumor natural product pladienolide. *Nat. Chem. Biol.* **3**, 570–575 (2007).

18. Teng, T. et al. Splicing modulators act at the branch point adenosine binding pocket defined by the PHF5A-SF3b complex. *Nat. Commun.* **8**, 15522 (2017).

19. Cretu, C. et al. Structural basis of splicing modulation by antitumor macrolide compounds. *Mol. Cell* **70**, 265–273 e268 (2018).

20. Hansen, S. R., Nikolai, B. J., Spreaker, P. J., Carrocci, T. J. & Hoskins, A. A. Chemical inhibition of pre-mRNA splicing in living Saccharomyces cerevisiae. *Cell Chem. Biol.* **26**, 443–448 e443 (2019).

21. Yokoi, A. et al. Biological validation that SF3b is a target of the antitumor macrolide pladienolide. *FEBS J.* **278**, 4870–4880 (2011).

22. Finci, L. I. et al. The cryo-EM structure of the SF3b spliceosome complex bound to a splicing modulator reveals a pre-mRNA substrate competitive mechanism of action. *Genes Dev.* **32**, 309–320 (2018).

23. Gozani, O., Potashkin, J. & Reed, R. A potential role for U2AF-SAP 155 interactions in recruiting U2 snRNP to the branch site. *Mol. Cell Biol.* **18**, 4752–4760 (1998).

24. Cass, D. M. & Berglund, J. A. The SF3b155 N-terminal domain is a scaffold important for splicing. *Biochemistry* **45**, 10092–10101 (2006).

25. Thickman, K. R., Swenson, M. C., Kabogo, J. M., Gryczynski, Z. & Kielkopf, C. L. Multiple U2AF65 binding sites within SF3b155: thermodynamic and spectroscopic characterization of protein-protein interactions among pre-mRNA splicing factors. *J. Mol. Biol.* **356**, 664–683 (2006).

26. Trappe, R. et al. Identification and characterization of a novel murine multigene family containing a PHD-finger-like motif. *Biochem. Biophys. Res. Commun.* **293**, 816–826 (2002).

27. Strikoudis, A. et al. Regulation of transcriptional elongation in pluripotency and cell differentiation by the PHD-finger protein Phf5a. *Nat. Cell Biol.* **18**, 1127–1138 (2016).

28. Yang, Y. et al. PHD-finger domain protein 5A functions as a novel oncoprotein in lung adenocarcinoma. *J. Exp. Clin. Cancer Res.* **37**, 65 (2018).

29. Zheng, Y. Z. et al. PHF5A epigenetically inhibits apoptosis to promote breast cancer progression. *Cancer Res.* **78**, 3190–3206 (2018).

30. Yang, Q. et al. Knockdown of PHF5A inhibits migration and invasion of HCC cells via downregulating NF-kappaB signaling. *Biomed. Res. Int.* **2019**, 1621854 (2019).

31. AlShareef, S. et al. Herboxidiene triggers splicing repression and abiotic stress responses in plants. *BMC Genomics* **18**, 260 (2017).

32. Ling, Y. et al. Pre-mRNA splicing repression triggers abiotic stress signaling in plants. *Plant J.* **89**, 291–309 (2017).

33. Butt, H. et al. CRISPR directed evolution of the spliceosome for resistance to splicing inhibitors. *Genome Biol.* **20**, 73 (2019).

34. Sun, C. The SF3b complex: splicing and beyond. *Cell Mol. Life Sci.* **77**, 3583–3595 (2020).

35. Schellenberg, M. J., Dul, E. L. & MacMillan, A. M. Structural model of the p14/SF3b155. branch duplex complex. *RNA* **17**, 155–165 (2011).

36. Bonnal, S., Vigevani, L. & Valcarcel, J. The spliceosome as a target of novel antitumour drugs. *Nat Rev. Drug Discov.* **11**, 847–859 (2012).

37. Amit, M. et al. Differential GC content between exons and introns establishes distinct strategies of splice-site recognition. *Cell Rep.* **1**, 543–556 (2012).

38. Vigevani, L., Gohr, A., Webb, T., Irimia, M. & Valcarcel, J. Molecular basis of differential 3' splice site sensitivity to anti-tumor drugs targeting U2 snRNP. *Nat. Commun.* **8**, 2100 (2017).

39. Feng, J. et al. SKIP confers osmotic tolerance during salt stress by controlling alternative gene splicing in Arabidopsis. *Mol. Plant* **8**, 1038–1052 (2015).

40. Zhan, X. et al. An Arabidopsis PWI and RRM motif-containing protein is critical for pre-mRNA splicing and ABA responses. *Nat. Commun.* **6**, 8139 (2015).

41. Thatcher, S. R. et al. Genome-wide analysis of alternative splicing during development and drought stress in maize. *Plant Physiol.* **170**, 586–599 (2016).

42. Jiang, J. et al. Integrating omics and alternative splicing reveals insights into grape response to high temperature. *Plant Physiol.* **173**, 1502–1518 (2017).

43. Carrasco-Lopez, C. et al. Environment-dependent regulation of spliceosome activity by the LSM2-8 complex in Arabidopsis. *Nucleic Acids Res.* **45**, 7416–7431 (2017).

44. Esteve-Bruna, D. et al. Prefoldins contribute to maintaining the levels of the spliceosome LSM2-8 complex through Hsp90 in Arabidopsis. *Nucleic Acids Res.* **48**, 6280–6293 (2020).

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.* **100**, 379–390 (2019).

46. Corrionero, A., Minana, B. & Valcarcel, J. Reduced fidelity of branch point recognition and alternative splicing induced by the anti-tumor drug spliceostatin A. *Genes Dev.* **25**, 445–459 (2011).

47. Bessonov, S. et al. Characterization of purified human Bact spliceosomal complexes reveals compositional and morphological changes during spliceosome activation and first step catalysis. *RNA* **16**, 2384–2403 (2010).

48. Rauhut, R. et al. Molecular architecture of the Saccharomyces cerevisiae activated spliceosome. *Science* **353**, 1399–1405 (2016).

49. Pham, D. & Koide, K. Discoveries, target identifications, and biological applications of natural products that inhibit splicing factor 3B subunit 1. *Nat. Prod. Rep.* **33**, 637–647 (2016).

50. Effenberger, K. A., Urabe, V. K. & Jurica, M. S. Modulating splicing with small molecular inhibitors of the spliceosome. *Wiley Interdiscip. Rev. RNA* **8**, e1381 (2017).

51. DeNicola, A. B. & Tang, Y. Therapeutic approaches to treat human spliceosomal diseases. *Curr. Opin. Biotechnol.* **6**, 72–81 (2019).

52. Jin, S. et al. Splicing factor SF3B1K700E mutant dysregulates erythroid differentiation via aberrant alternative splicing of transcription factor TAL1. *PLoS ONE* **12**, e0175523 (2017).

53. Serrat, X. et al. CRISPR editing of sftb-1/SF3B1 in Caenorhabditis elegans allows the identification of synthetic interactions with cancer-related mutations and the chemical inhibition of splicing. *PLoS Genet.* **15**, e1008464 (2019).

54. Butt, H., Zaidi, S. S., Hassan, N. & Mahfouz, M. CRISPR-based directed evolution for crop improvement. *Trends Biotechnol.* **38**, 236–240 (2019).

55. Zhang, Y. & Qi, Y. CRISPR enables directed evolution in plants. *Genome Biol.* **20**, 83 (2019).

56. 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 (Basel)* **10**, 596 (2019).

57. Dong, C. et al. Alternative splicing plays a critical role in maintaining mineral nutrient homeostasis in rice (Oryza sativa). *Plant Cell* **30**, 2267–2285 (2018).

58. Horan, K. et al. Annotating genes of known and unknown function by large-scale coexpression analysis. *Plant Physiol.* **147**, 41–57 (2008).

59. Warkocki, Z. et al. Reconstitution of both steps of Saccharomyces cerevisiae splicing with purified spliceosomal components. *Nat. Struct. Mol. Biol.* **16**, 1237–1243 (2009).

60. Liu, H. L. & Cheng, S. C. The interaction of Prp2 with a defined region of the intron is required for the first splicing reaction. *Mol. Cell Biol.* **32**, 5056–5066 (2012).

61. Liu, Y. C. & Cheng, S. C. Functional roles of DExD/H-box RNA helicases in pre-mRNA splicing. *J. Biomed. Sci.* **22**, 54 (2015).

62. Yeap, W. C., Namasivayam, P. & Ho, C. L. HnRNP-like proteins as post-transcriptional regulators. *Plant Sci.* **227**, 90–100 (2014).

63. Howard, J. M. & Sanford, J. R. The RNAissance family: SR proteins as multifaceted regulators of gene expression. *Wiley Interdiscip. Rev. RNA* **6**, 93–110 (2015).

64. Lata, C. & Prasad, M. Role of DREBs in regulation of abiotic stress responses in plants. *J. Exp. Bot.* **62**, 4731–4748 (2011).

65. Todaka, D., Nakashima, K., Shinozaki, K. & Yamaguchi-Shinozaki, K. Toward understanding transcriptional regulatory networks in abiotic stress responses and tolerance in rice. *Rice (NY)* **5**, 6 (2012).

66. Nakashima, K., Ito, Y. & Yamaguchi-Shinozaki, K. Transcriptional regulatory networks in response to abiotic stresses in Arabidopsis and grasses. *Plant Physiol.* **149**, 88–95 (2009).

67. Filichkin, S. A. et al. Abiotic Stresses Modulate Landscape of Poplar Transcriptome via Alternative Splicing, Differential Intron Retention, and Isoform Ratio Switching. *Front. Plant Sci.* **9**, 5 (2018).

68. Lee, S. C. W. & Abdel-Wahab, O. Therapeutic Targeting of Splicing in Cancer. *Nat. Med.* **22**, 976–986 (2016).

69. Shen, Y. et al. Global dissection of alternative splicing in paleopolyploid soybean. *Plant Cell* **26**, 996–1008 (2014).
70. 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.* **24**, 496–506 (2019).

71. Xin, R., Kathare, P. K. & Huq, E. Coordinated Regulation of Pre-mRNA Splicing by the SFPS-RRC1 Complex to Promote Photomorphogenesis. *Plant Cell* **31**, 2052–2069 (2019).

72. Xie, K., Minkenberg, B. & Yang, Y. Boosting CRISPR/Cas9 multiplex editing capability with the endogenous tRNA-processing system. *Proc. Natl Acad. Sci. USA* **112**, 3570–3575 (2015).

## Acknowledgements

We would like to thank members of the genome engineering and synthetic biology laboratory at KAUST for their critical discussion and technical help in this work. This work is funded by KAUST-baseline funding to M.M.

## Author contributions

M.M. conceived the project. M.M., H.B., S.A., and A.E. designed and conducted the experiments. J.B. and H.B. analyzed the data. H.B., J.B., M.C., A.R., and M.M. designed and analyzed the splicing experiments. M.M., H.B., J.B., M.C., and A.R. wrote the paper. All authors read and approved the final manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s42003-021-02051-y.

**Correspondence** and requests for materials should be addressed to M.M.M.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021