CHERP Regulates the Alternative Splicing of pre-mRNAs in the Nucleus

Calcium homeostasis endoplasmic reticulum protein (CHERP) is colocalized with the inositol 1,4,5-trisphosphate receptor (IP3R) in the endoplasmic reticulum or perinuclear region, and has been involved in intracellular calcium signaling. Structurally, CHERP carries the nuclear localization signal and arginine/serine-dipeptide repeats, like domain, and interacts with the spliceosome. However, the exact function of CHERP in the nucleus remains unknown. Here, we showed that poly(A)+ RNAs accumulated in the nucleus of CHERP-depleted U2OS cells. Our global analysis revealed that CHERP regulated alternative mRNA splicing events by interaction with U2 small nuclear ribonucleoproteins (U2 snRNPs) and U2 snRNP-related proteins. Among the five alternative splicing patterns analyzed, intron retention was the most frequently observed event. This was in accordance with the accumulation of poly(A)+ RNAs in the nucleus. Furthermore, intron retention and cassette exon choices were influenced by the strength of the 5′ or 3′ splice site, the branch point site, GC content, and intron length. In addition, CHERP depletion induced anomalies in the cell cycle progression into the M phase, and abnormal cell division. These results suggested that CHERP is involved in the regulation of alternative splicing.


Introduction
Genetic information is stored in genomic DNA. In higher eukaryotes, immature pre-mRNA is transcribed from genomic DNA by RNA polymerase II and undergoes a series of processing events, including capping, splicing and polyadenylation. The cap is added at the 5 terminus while the RNA polymerase II pauses. Then, the transcription elongation continues. Introns are sometimes removed concomitantly to elongation in order to generate spliced mRNA. The 3 terminus becomes polyadenylated once the RNA polymerase II has reached it. Mature mRNAs translocate to the cytosol through the nuclear pore complex to be translated into proteins.
It is important to recognize exons, introns and their boundaries, which are regulated by cis-regulatory elements [1][2][3]. These elements are classified into splicing signals and splicing regulatory elements. Splicing signals, such as the 5 splice site, the 3 splice site, the branch point site and the polypyrimidine tract are associated with the spliceosome [4]. Splicing regulatory elements, such as the exonic splicing enhancer (ESE), the exonic splicing silencer (ESS), the intronic splicing enhancer (ISE) and the intronic splicing silencer (ISS) are associated with specific trans-acting factors, like U small nuclear ribonucleoproteins (snRNPs), serine and arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs) [5].
Splicing requires five U small nuclear RNAs (snRNAs), U1, U2, U4, U5 and U6, and more than 240 proteins [4,6,7]. These components form several sub-complex snRNPs that regulate splicing by binding to specific pre-mRNA sites within the spliceosome. There are two types of mRNA splicing, e.g., constitutive and alternative splicing [8]. The former generates one transcript by excising the introns and combining the exons, whereas the latter generates several transcripts obtained from multiple exon combinations. Thus, alternative splicing produces a variety of mRNAs from a pre-mRNA [9]. It is regulated in tissue, cell and developmental stage-specific manners and is involved in producing the transcript variants from about 90% of human genes with multiple exons. Therefore, this system is essential in humans. The pattern of alternative splicing can be primarily categorized into five groups: cassette exon, alternative 5 splice site, alternative 3 splice site, mutually exclusive exons and intron retention [9]. Beside U snRNPs, numerous cis-regulatory elements and trans-acting factors regulate precisely alternative splicing [1]. For example, SR proteins and hnRNPs are implicated in alternative splicing [10,11]. SR proteins are generally thought to activate splicing by binding to ESE and ISE, whereas hnRNPs inhibit splicing by binding to ESS and ISS [2]. Other RNA binding proteins have also been shown to affect alternative splicing [2,12].
Previously, it has been reported that calcium homeostasis endoplasmic reticulum protein (CHERP) is localized in the endoplasmic reticulum and the perinuclear region, and colocalizes with the inositol 1,4,5-trisphosphate receptor (IP3R) [13]. CHERP interacts with a calcium channel of the endoplasmic reticulum, the ryanodine receptor 1 (RyR1), in the skeletal muscle [14], suggesting that it has a function in cell cycle control, cell growth and intracellular calcium signaling by regulating intracellular calcium concentrations [14]. On the other hand, CHERP harbors a nuclear localization signal and multiple structural features related to RNA processing, including an arginine/serine-dipeptide repeat (RS) domain and a glycine-rich motif (G-patch). It also interacts with the spliceosome, especially the U2 snRNP-related factor, implying that CHERP participates in the alternative splicing as an SR-like protein [15]. In addition, CHERP is associated with the apoptosis-linked gene 2 (ALG-2) in nuclear speckles in a Ca 2+ concentration-dependent manner, and affects the alternative splicing of exon 41 and 42 of the inositol triphosphate receptor type 1 (IP 3 R1) [16]. Recently, it has been found that CHERP is involved in the regulation of splicing by interacting with RNA-binding motif protein 17 (RBM17) and the U2 snRNPassociated SURP domain-containing protein (U2SURP) [17]. However, the precise function of CHERP in the nucleus remains to be determined.
Here, we investigated the function of CHERP in the nucleus and found that CHERP depletion induced an accumulation of poly(A) + RNAs in the nucleus. To evaluate the causes of this nuclear accumulation, we performed RNA sequencing (RNA-seq) analyses, which indicated that the depletion of CHERP induced a dramatic change in alternative splicing, in particular, resulting in an increased retention of introns in transcripts. The properties of the intron retention and other alternative splicing events as well as the cell fate were also analyzed.

CHERP Is Involved in Splicing in the Nucleus
To determine CHERP localization, cellular fractionation was performed and analyzed by Western blot. CHERP was found in the nuclear fraction of U2OS cells (Figure 1a), prompting us to investigate its function in mRNA metabolism. To this aim, we used siRNAs, which induced a marked decrease of CHERP expression (Figures 1b and S1a), and assessed poly(A) + RNA localization (Figures 1c,d and S1b). Poly(A) + RNAs clearly accumulated in the nucleus after the depletion of CHERP to a similar extent as the accumulation induced by silencing the U2 auxiliary factors 65-associated protein (UAP56) or the nuclear RNA export factor 1 (NXF1), which are both involved in mRNA export. These results suggested that CHERP was involved in mRNA processing or nuclear export. To elucidate the function of CHERP in mRNA nuclear metabolism, we identified proteins interacting with CHERP by immunoprecipitating CHERP from cells stably expressing FLAG-tagged CHERP with an anti-FLAG antibody, and analyzing the interaction partners by liquid chromatography-tandem mass spectrometry (LC-MS/MS). The experiments were performed twice ( Figure S1c,d). The proteins detected are listed in Table  S1. These were mostly categorized according to functions related to mRNA processing, especially to splicing (Figures 1e and S1e). Among them, components of U2 snRNPs were the most frequently detected. We confirmed that the U2 auxiliary factors U2AF1 (U2AF35) and U2AF2 (U2AF65) interacted with CHERP by Western blot ( Figure S1f). These results implicated that CHERP likely had a role in an early step of splicing. To evaluate the exact influence of CHERP on splicing events, transcriptome analyses of CHERP-depleted U2OS cells were performed (Figures 1f,g and S1g,h). We classified the alternative splicing events into five categories using rMATS. The frequency of splicing changes was measured, and an event was considered significant for false discovery rate (FDR) < 0.05. CHERP depletion elicited an alteration of all alternative splicing events, with retained introns being the most frequent.

Identification and Characterization of Intron Retention after CHERP Depletion
To investigate the characteristics of introns retained in response to CHERP depletion, we considered separately increased intron inclusions (III) and decreased intron inclusions (DII) and found 290 III and 389 DII cases for FDR < 0.05 (Figures 1e and 2a), respectively. We measured the levels of mRNAs with III and DII events in the cytoplasm. As expected, CHERP depletion induced a tendency toward a decrease of the levels of transcripts with III, whereas the amounts of mRNAs with DII were clearly increased in the cytoplasm of CHERP-depleted cells (Figures 2b and S2a). Next, we validated the results regarding retained introns. To this aim, reverse transcription-polymerase chain reactions (RT-PCRs) were performed using a series of polymerase chain reaction (PCR) primer pairs targeting two adjacent constitutive exons of a selected intron. III and DII were observed in the CHERP-depleted condition (Figure 2c,d). The introns in calcium-activated nucleotidase 1 (CANT1) and arginine and serine rich protein 1 (RSRP1), which typically are III, were increased, and introns of disheveled segment polarity protein 3 (DVL3), like podocalyxin (PODXL) and STAGA complex 65 subunit gamma (SUPT7L), which are usually DII, were reduced. Gene ontology (GO) analysis of transcripts with III and DII in CHERP-depleted cells are shown in Figure 2e,f. The III-containing transcripts in CHERP-depleted cells were related to RNA metabolism, whereas a variety of gene functions were affected by DII. GOsenriched transcripts with retained introns were similar to those enriched in DII-containing mRNAs, rather than to those with III. Thus, poly(A) + RNAs accumulated in the nucleus of CHERP-depleted cells seemed to be linked to III events. Moreover, dividing retained intron events into III and DII allowed identification of the poly(A) + RNA that accumulated. cells are shown in Figure 2e,f. The III-containing transcripts in CHERP-depleted cells were related to RNA metabolism, whereas a variety of gene functions were affected by DII. GOs-enriched transcripts with retained introns were similar to those enriched in DII-containing mRNAs, rather than to those with III. Thus, poly(A) + RNAs accumulated in the nucleus of CHERP-depleted cells seemed to be linked to III events. Moreover, dividing retained intron events into III and DII allowed identification of the poly(A) + RNA that accumulated.  (e,f) Gene ontology (GO) groups enriched in transcripts with III (e) or DII (f) after CHERP depletion are shown. FDR < 0.05, inclusion or skipping count ≥ 10, and IncLevelDifference ≤ -0.1 for RNA sequencing (RNA-Seq) data was considered statistically significant.

CHERP Regulates Intron Retention in the Target Transcripts
To determine which mechanisms were triggered by CHERP to alter splicing events, we investigated the properties of introns. In both III and DII, the 5 and 3 splice site scores of introns affected by CHERP depletion were weaker than those obtained for the reference genes (Figure 3a,b). The branch point site scores were not significantly changed, and the polypyrimidine tract site scores in III were weaker than those of the reference genes ( Figure 3c,d). Furthermore, the length of introns affected by CHERP depletion was significantly shorter than that of the reference genes ( Figure 3e). The DII regulated by CHERP had a higher GC content in introns than that in reference genes ( Figure 3f). Short introns with high GC content have been associated with intron inclusion events [18]. As splicing efficiency is based on the splice site, branch point site and polypyrimidine tract site scores, we analyzed the nucleotides in CHERP-targeted 5 and 3 splice sites ( Figure 3g). Both sides of the intronic sequences in the reference genes were rich in thymine residues. In contrast, III and DII induced by CHERP depletion were slightly different from those of the reference genes. To examine whether the weak splice sites were targeted by CHERP, we constructed a retinitis pigmentosa GTPase regulator (RPGR) minigene as a target for CHERP. Additionally, mutant minigenes, containing 5 and 3 splice site mutations that rendered the splice sites stronger, were generated. Plasmids with wild-type or mutant minigenes were transfected to cells, and transgene products were analyzed by RT-PCR. The splicing in the wild-type minigenes, but not in the mutant minigenes, was dependent on CHERP ( Figure 3h). These data suggested that CHERP functioned as a regulator involved in the exact splice site selection.

Identification and Characterization of Cassette Exons Regulated by CHERP
The number of altered exon inclusions (EI) and exon skipping (ES) in cassette exons regulated by CHERP were counted. We found 2945 cases of EI and 3102 instances of ES with FDR < 0.05 (Figures 1e and 4a), implicating that CHERP did not selectively affect EI or ES. EI tended to be increased in cytoplasmic target transcripts, whereas ES was not dramatically affected (Figures 4b and S4a). To validate the cytoplasmic expression of alternatively spliced transcripts, we designed PCR primer pairs located in two exons adjacent to the target exon.   (a-f) FDR < 0.05, inclusion or skipping count ≥ 10, and IncLevelDifference ≤ -0.1 (III) or IncLevelDifference ≥ 0.1 (DII) were considered significantly different. All introns were used as reference introns. For statistical analysis, a Wilcoxon rank sum test was used to calculate p-values. * p < 0.05, *** p < 0.001. (g) Nucleotide sequences of exon-intron boundaries in III, DII and references, estimated with BedTools. (h) Splicing assay using retinitis pigmentosa GTPase regulator (RPGR) minigenes containing CHERP target intron. Red letters indicated bases and scores that introduce mutation strengthening 5 or 3 splice site.

CHERP Regulates Cassette Exon Inclusion and Skipping in Target Genes
To elucidate the specificity of CHERP-regulated cassette exons, we performed a correlation analysis between the splice site properties and CHERP-dependent or -independent genes. In EI and ES, the 5′ splice sites in introns upstream and downstream of targeted

CHERP Regulates Cassette Exon Inclusion and Skipping in Target Genes
To elucidate the specificity of CHERP-regulated cassette exons, we performed a correlation analysis between the splice site properties and CHERP-dependent or -independent genes. In EI and ES, the 5 splice sites in introns upstream and downstream of targeted exons were not significantly different from those of reference genes, except for 5 splice sites in introns upstream of ES (Figure 5a). By contrast, 3 splice sites in introns upstream and downstream of targeted exons were weaker than those of reference genes (Figure 5b), implicating that U2 snRNP recruitment to 3 splice sites of upstream and downstream introns might be reduced. Branch point site scores were significantly weaker for upstream and downstream introns targeted by CHERP. Polypyrimidine tract site scores were not altered compared with those of reference genes (Figure 5c,d). These results suggest that the environment around CHERP-targeted exons was less likely to be spliced than reference exons. CHERP-regulated introns upstream and downstream from targeted exons were significantly longer than those of reference genes (Figure 5e). Furthermore, in introns upstream and downstream from targeted exons, the GC content was lower in EI and ES than that of reference genes (Figure 5f), suggesting that exons surrounded by longer introns with lower GC content tended to be targets of CHERP in skipped exon (SE). These intron profiles were different from those of III and DII (compare with Figure 3e,f). Then, we calculated the nucleotide usage of 5 and 3 splice sites for introns upstream and downstream from CHERP-targeted exons. CHERP-targeted nucleotide sequences were not much different from that in reference genes, although 5 and 3 splice site scores seemed lower (Figure 5g). Therefore, a slight difference in nucleotides compared with reference genes might affect CHERP-dependent EI and ES. To examine this hypothesis, we constructed autophagy related 16 like 1 (ATG16L1) and serine/threonine-protein kinase 3 (SIK3) minigenes for EI and ES analysis, respectively. Mutations leading to stronger or weaker splice sites were introduced to 5 or 3 splice sites. When strengthening splice site mutations were introduced in the 3 splice sites, the target exon was fully recognized independently of CHERP depletion (Figure 5h). By contrast, weakening mutations in the 3 splice sites of the target exons led to the exons being fully skipped, regardless of CHERP depletion. Similar results were obtained with mutations introduced in the 5 splice sites. Among these mutations, four mutants consisted of a single nucleotide change, supporting the hypothesis that nucleotide sequences around the 5 and 3 sites of target exons were sensitive to CHERP.

CHERP Regulates Cell Survival and Cell Death
The GO analysis revealed an enrichment of cell cycle genes in transcripts with EI induced by CHERP depletion, and those with ES were enriched in cell cycle or division, implying that CHERP played a role in cell cycle progression and cell viability. To investigate the role of CHERP in cell function, the cytoplasmic mRNA transcriptome was analyzed in cells with or without CHERP depletion. CHERP depletion induced the increase of the expression of 3934 transcripts and a decreased expression of 2822 mRNAs (Figure 6a). To determine the cellular function of CHERP depletion, the GO analysis of transcripts with an expression increased more than 1.5 fold and less than 0.67 fold is shown (Figure S6a,b). Several transport functions were enriched, whereas the representation of cell cycle-related genes was reduced when CHERP was depleted. We then investigated whether the function of CHERP was associated with cell cycle or cell division. U2OS cells stably expressing histone 2B-green fluorescent protein (H2B-GFP) were used to analyze cell division by live cell imaging (Figure 6b). In CHERP-depleted cells, centromeres were not well aligned, and the time of mitosis was often delayed during the prometaphase. aurora kinase B (AURKB) has a key role in the alignment of chromosomes at prometaphase and AURKB mRNA expression was markedly reduced in CHERP-depleted cells (Figure 6c). CHERP depletion also reduced the cytoplasmic expression of the Survivin/baculoviral IAP repeat containing 5 (BIRC5), AURKB, centromere protein A (CENPA) and breast cancer type 1 (BRCA1) mRNAs, which function in the progression of the M phase, prompting us to investigate the corresponding protein expression. BIRC5, AURKB, CENPA and BRCA1 protein levels were reduced in CHERP-depleted cells (Figure 6d). Moreover, CHERP depletion not only resulted in a weak expression of AURKB, but also in its aberrant distribution, as shown by immunostaining (Figure 6e). CHERP depletion also led to the aberrant distribution of α-Tubulin. These mislocalizations might be the cause of the delay of the M phase progression. The viability of CHERP-depleted cells was significantly reduced (Figure 6f), suggesting that CHERP was indispensable for cellular function, likely by regulating alternative splicing.

Discussion
Here, the function of CHERP in mRNA metabolism was investigated. RNA fluorescent in situ hybridization (RNA-FISH) analyses revealed that poly(A) + RNAs accumulated in the nucleus of CHERP-depleted cells, suggesting that CHERP was involved in mRNA metabolism in the nucleus [19,20]. We then performed LC-MS/MS analyses to assess CHERP functions by identifying its interaction partners. CHERP interacted strongly with U2 snRNP components, such as spliceosome factor 3b subunits SF3B1, SF3B2 and SF3B3, and U2-related proteins, such as U2AF2 and U2SURP, indicating that CHERP participated in mRNA splicing [21,22]. Our results are consistent with the results of Lin-Moshier et al., who used HEK293 cells stably expressing GFP-CHERP [15]. De Maio et al. also showed, by immunoprecipitation, that CHERP forms a complex with RBM17 and U2SURP to induce splicing changes [17]. RBM17 and U2SURP were detected as a CHERP interaction partner in our LC-MS/MS experiments, indicating that our data were solid.
Poly(A) + RNAs accumulated in the nucleus of CHERP-depleted cells and CHERP was associated with U2 snRNP components. This prompted us to investigate CHERP's role in mRNA splicing using next generation sequencing (NGS) and rMATS. Among alternative splicing events, retained introns were the most frequently observed, which was consistent with our observation that CHERP depletion induced poly(A) + RNA accumulation in the nucleus [18,23]. Although the transcriptome analysis revealed that all types of alternative splicing were affected by CHERP depletion, which was in agreement with the De Maio et al. paper, a detailed analysis showed that CHERP depletion caused retained introns more frequently, although this was not observed by De Maio et al. [17]. This discrepancy might be due to the different cell types used in both analyses. This possibility is supported by the recent findings that CHERP depletion in Hela cells frequently caused retained intron events [24] and implies that CHERP regulated the alternative splicing in a cell-specific manner. In our analysis, weak or suboptimal splice sites, high GC content in introns, and short intron lengths were important and common features of intron retention. Interestingly, the depletion of another RNA binding protein, SON, also resulted in the nuclear retention of introns with the same features [23,25]. It has been reported that CHERP forms a complex with RBM17, an RNA-binding protein, and U2SURP, an SR protein [17], suggesting that CHERP does not function as a core splicing factor but rather as a modulating factor [26,27].
Alternative splicing is, importantly, regulated by the sequences of 5 and 3 splicing sites and their combinations [6,7,23,[28][29][30]. In addition, cis-acting RNA elements are involved in splicing regulation [31][32][33][34][35]. Cis-acting RNA elements often become the target sites of trans-acting factors [36][37][38][39][40][41][42][43]. The status and positional relations of these signals have a significant effect on splicing. Thus, alternative splicing is intricately regulated. We have shown that the GC content is low in the upstream and downstream introns of EI induced by CHERP depletion. This seemed to be different from the data obtained for CHERP-regulated retained introns, which had high GC contents. Interestingly, it has been reported that CHERP binds to specific adenine-rich sequences on exons and causes splicing changes [44]. Thus, CHERP might participate in the regulation of EI and ES by binding to adenine-rich sequences.
The depletion of CHERP induced alternative splicing changes, in particular of factors promoting the entrance into the M phase, including BIRC5, AURKB, CENPA, and BRCA1, which results in the reduced cytoplasmic expression of these factors [45][46][47]. Our data showed that the cellular expression of these factors was reduced, and AURKB and α-Tubulin were totally mislocalized, implying that CHERP targets genes involved in the progression into the M phase. The implication of the CHERP association with RBM17 and U2SURP in the regulation of splicing events, such as RI and SE, will be investigated next by depleting these factors to uncover a more specific role of CHERP in splicing.

Reagents and Instruments
The following commercial chemicals were used: Dulbecco's Modified Eagle's Medium

Plasmid DNAs, siRNAs and Oligo DNA Primers
The plasmid DNAs, siRNAs and oligo DNA primers used in this study are listed in Tables S2-S4. We used pcDNA5/FRT/TO and siRNAs from Invitrogen, and siRNAs and oligo DNA primers from Integrated DNA Technologies (Coralville, IA, USA).

Construction of Plasmid DNAs
To generate inducible stable cell lines, wild-type human CHERP cDNA was inserted into pcDNA5 FLAG FRT/TO. Cloning sites of individual plasmid DNA are described in Table S2. The template DNA containing full-length CHERP cDNA was obtained from Dr. Masatoshi Maki. Wild-type human CHERP was generated by PCR. Sequences were confirmed by sequencing.
For minigene splicing assays, the target region was amplified from human genomic DNA by PCR using KOD Fx Neo according to the manufacturer's instructions, and cloned into pcDNA5 3 × FLAG FRT/TO. Mutant minigenes were obtained by mutation PCR. Sequences were confirmed by sequencing.

Establishment of Stable Cell Lines
To establish stable cell lines expressing FLAG-CHERP wild-type, CHERP expression vector, pcDNA5 FLAG FRT/TO, was transfected with Lipofectamine 2000 according to the manufacturer's instructions. Successfully transfected cells were selected using 100 ng/mL hygromycin in DMEM containing 10% FBS. A colony was picked up and expanded in a 6-well plate. To induce the expression of FLAG-CHERP, 1 µg/mL of doxycycline was added to the culture medium for at least 24 h. Cell lysates were separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and analyzed by Western blot. U2OS cells expressing H2B-GFP were described previously [48,49].

Immunostaining
Cells (1 × 10 5 cells/mL) were seeded and cultured on coverslips. After 24 h, cells were fixed with 4% formaldehyde in cultured medium. Then, cells were washed once with phosphate-buffered saline (PBS) and blocked with PBS containing 6% bovine serum albumin for 1 h. The coverslips were incubated with anti-CHERP antiserum in PBS containing 2% bovine serum albumin, washed with PBS three times, and incubated with goat anti-rabbit antibody labeled with Alexa Fluor 594. DNA was counterstained with DAPI.

RNA-FISH
RNA-FISH was performed as described previously [49,50]. Briefly, cells (1 × 10 5 cells/mL) were seeded and cultured on coverslips. After 24 h, cells were transfected with 20 nM siRNAs using Lipofectamine 2000, according to the manufacturer's instructions. After 48 h, cells were fixed with 10% formaldehyde and permeabilized with PBS containing 0.1% toctylphenoxypolyethoxyethanol/Triton TM X-100. The buffer was replaced with 2 × salinesodium citrate buffer (SSC). Then, cells were prehybridized with ULTRAhyb TM -Oligo Hybridization Buffer and incubated with 20 pmol Alexa Fluor 594-labeled oligo-d(T) 45 probe overnight. Cells were washed successively with 2 × SSC, 0.5 × SSC and 0.1 × SSC. Images were taken at random using a Zeiss Axioplan 2 (Carl Zeiss, Jena, Germany) with an OLYMPUS DP70 camera (OLYMPUS, Tokyo, Japan). Segmentation of the nuclei and the cytoplasm and quantification of poly(A) + RNA signals were performed using the CellProfiler cell image analysis software (version 3.1.5) [51].

Preparation of RNA from Cytoplasmic, Nuclear and Whole-Cell Fractions
For the preparation of cytoplasmic and nuclear fractions, cells were scraped, washed with PBS, and treated with a cell lysis buffer containing 20 mM tris(hydroxymethyl) aminomethane (Tris)-HCl (pH 8.0), 200 mM NaCl, 1 mM MgCl 2, and 1% octylphenoxy poly(ethyleneoxy)ethanol/Nonidet TM P-40 (NP-40) on ice for 5 min. After centrifugation at 10,000× g for 1 min, the supernatant represented the cytoplasmic fractions. The pellets represented the nuclear fractions. Whole-cell fractions were prepared by lysing the cells as described previously [20,49].

Protein Preparation from Cytoplasmic and Nuclear Extracts
Proteins were extracted from the cytoplasmic fraction and nuclear extraction as follows. Cells were centrifuged and the pellet was suspended very carefully in three times its volume of lysis solution containing 10 mM 4-(2-hydroxyethyl) piperazine-1-ethanesulfonic acid (HEPES)-KOH (pH 7.9), 10 mM KCl, 1.5 mM MgCl 2 , 0.2 mM phenylmethylsulfonyl fluoride (PMSF), and 0.5 mM dithiothreitol (DTT). The cell solutions were then incubated on ice for 10 min, vortexed for 5 s, and centrifuged at 10,000× g for 10 s. The supernatants were collected as cytoplasmic fractions. The pellets were resuspended in an equal volume of solution containing 20 mM HEPES-KOH (pH 7.9), 25% glycerol, 420 mM NaCl, 1.5 mM MgCl 2 , 0.2 mM EDTA, 0.2 mM PMSF and 0.5 mM DTT. The nuclei were lysed for 20 min on ice and centrifuged at 10,000× g for 10 min. The resulting supernatants were the nuclear fractions. The protein content in each fraction was determined by the Bradford assay (Nacalai tesque). Protein samples were mixed with a 4 × SDS buffer (190 mM Tris-HCl [pH 6.8], 40% glycerol, 0.8% SDS, 0.2% bromophenol blue and 40 mM DTT) and boiled for 2 min.

Immunoprecipitations
Protein nuclear extracts were prepared from HEK293 Flp-In TM T-REx TM cells stably expressing FLAG-CHERP and immunoprecipitation was performed as described previously [30]. Briefly, samples containing nuclear extract and RNase A, for the removal of RNA, were incubated for 20 min at 30 • C and then centrifuged to remove precipitants. The supernatants were mixed with 10 µL anti-FLAG M2 antibody beads. Beads were rotated overnight at 4 • C and washed with PBS containing 0.1% Triton TM X-100, 0.2 mM PMSF and 0.5 mM DTT. The proteins on the beads were recovered in SDS buffer containing 250 mM Tris-HCl, 1% SDS, 0.002% bromophenol blue and 40% glycerol. After boiling for 2 min, the eluates were transferred to a new tube and separated by SDS-PAGE.

Silver Staining
The protein samples were mixed with 4× SDS buffer and boiled for 2 min. Protein separation was performed using SuperSep™ Ace, 5-20%, 17-well. Silver staining of SDS-PAGE gels was carried out with Silver Staining Kit Protein, in accordance with the manufacturer's protocol. Pictures were captured using the image analyzer LAS-1000 (FUJIFILM, Tokyo, Japan).

LC-MS/MS Analyses
LC-MS/MS analyses were described previously [50]. Results were analyzed using the number of peptides identified prot_matches and prot_score, which was calculated with the Mascot software (Matrix Science, London, UK). Protein scores were calculated by subtracting the prot_score of the controls from that of the CHERP-expressing cells. The identified CHERP-interacting proteins are listed in Table S1.

Western Blots
Samples were separated by SDS-PAGE and transferred to a PVDF membrane. Membranes were blocked with 5% skim milk and then incubated with a primary antibody overnight at 4 • C. After washing, membranes were incubated with a secondary antibody conjugated with horseradish peroxidase for 2 h at room temperature. After careful washing, the signals were detected using ECL substrate with LAS 4000 mini (GE Healthcare).

RNA Isolation, cDNA Synthesis and PCR
Total or cytoplasmic RNAs were isolated from whole cells or the cytoplasmic fraction, respectively, with Sepasol ® I Super G, and then treated with RQ1 DNase-Free RNase to remove genomic DNA, according to the manufacturer's instructions. Then, cDNAs were synthesized by reverse transcription using Rever Tra Ace ® , according to the manufacturer's instructions. PCRs were performed using KOD FX NEO and analyzed by electrophoresis using either 1% agarose gel with a Tris-acetate-EDTA (TAE) buffer, or 8% polyacrylamide gel with a Tris-borate-EDTA (TBE) buffer. PCR products were visualized by a FAS-IV gel imaging system (Nippon Genetics, Tokyo, Japan). Amplified PCR products were confirmed by sequencing.

Preparation of RNA-seq Library
U2OS cells (30% confluent) were transfected with 20 nM EGFP or CHERP#2 siRNA for 48 h. Whole and cytoplasmic RNAs were prepared as described above. The RNA-seq library was prepared using Truseq Stranded mRNA Library Prep Kit (Illumina, San Diego, CA, USA).

Bioinformatic Analyses
RNA-seq analyses of two biological replicates for each group were performed using the NextSeq High Output sequencer (Illumina). Sequence reads were mapped to the hg19 human genome using STAR (version 2.5.2b) [52]. Mapped reads were transformed to expression counts for each transcript with RSEM (version 1.2.25) [53] and a gtf file of GENCODE (version 19). Normalized counts of gene expression for each replicate were obtained by processing EBSeq (version 1.2.25) software [54]. Alternative splicing events were analyzed by rMATS [55]. Significant splicing events were captured for FDR < 0.05. Mapping results were indexed using SAMtools (version 1.9) [56] and visualized using the Integrative Genomics Viewer (IGV: version 2.3) (http://software.broadinstitute.org/software/igv/home, accessed on 25 January 2022). Gene ontology (GO) analysis was performed using Database for Annotation, Visualization and Integrated Discovery (DAVID: version 6.7) [57,58]. The splice site score was calculated using MaxEntScan [59], and branch point and polypyrimidine tract score were determined using the SVM-BPfinder [60] program. To analyze the sequence usage, 15-and 20-nucleotide sequences of exon and intron boundaries at 5 and 3 splice sites, respectively, were obtained by BEDtools (version 2.24.0) [61]. Motifs were visualized using Weblogo (version 2.8.2) [62]. Sequence reads were deposited in NCBI under accession number GSE193399.

Quantitative PCR (qPCR)
qPCRs were performed with the Thermal Cycler Dice Real Time System (Takara BIO) using Thunderbird SYBR qPCR Mix. Each sample was prepared in triplicate. The mRNA amounts were calculated by threshold cycle (Ct) values.

Minigene Splicing Reporter Assay
Cells were cotransfected with siRNAs and minigene reporter plasmid DNAs. Briefly, cells were seeded at 15-20% confluency per well in 6-well plates. The next day, cells were transfected with siRNAs using Lipofectamine 2000. After 48 h of incubation, cells were transfected with minigene reporter plasmid DNAs (750 ng plasmid DNAs per well) using Lipofectamine 2000. After one day of incubation, cells were collected for the preparation of total RNA, as described above. Expression of the minigene reporter was analyzed using RT-PCR, followed by electrophoresis.

Live Cell Imaging
Cells expressing H2B-GFP were seeded on 35 mm-diameter glass bottom dishes (AGC techno glass, Shizuoka, Japan) and transfected with siRNAs using Lipofectamine 2000, as reported previously [20,48]. After siRNA transfection, cells were cultured for 36 h. Images were taken every 5 min using Biostation IMq (Nikon, Tokyo, Japan) at 37 • C in a humidified chamber with 5% CO 2 for 16 h.

MTT Assays
Cell viability was assayed by colorimetric MTT assay [50,63]. U2OS cells were seeded at 7.5 × 10 3 cells/mL in a 96-well plate and cultured for 24 h. siRNAs were introduced to each well and incubated for 48 h. Afterward, 5 µL of 5 mg/mL MTT was added to each well and incubated for 4 h. The cell culture plate was then centrifuged at 400× g for 5 min. The supernatant was removed and cells were solubilized with 10 mM NH 4 Cl containing 10% SDS (pH 7.0). Cell viability was estimated by measuring the optical density at 600 nm.