The SF3B1R625H Mutation Promotes Prolactinoma Tumor Progression Through Aberrant Splicing Of DLG1

The SF3B1R625H Mutation Promotes Prolactinoma Tumor Progression Through Aberrant Splicing Of DLG1 | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article The SF3B1 R625H Mutation Promotes Prolactinoma Tumor Progression Through Aberrant Splicing Of DLG1 Jing Guo, Qiuyue Fang, Yulou Liu, Dawei Wang, Chuzhong Li, Weiyan Xie, and 1 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-934219/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 17 Jan, 2022 Read the published version in Journal of Experimental & Clinical Cancer Research → Version 1 posted 10 You are reading this latest preprint version Abstract Background Recently, a hotspot mutation in prolactinoma was observed in splicing factor 3b subunit 1 ( SF3B1 R625H ), but its functional effects and mechanisms are poorly understood. Methods Using the CRISPR/Cas9 genome editing system and rat pituitary GH3 cells, we generated heterozygous Sf3b1 R625H mutant cells. Sanger and whole-genome sequencing were conducted to verify the introduction of this mutation. Transcriptome analysis was performed in SF3B1 -wild-type versus mutant human prolactinoma samples and GH3 cells. Quantitative PCR and minigene reporter assays were conducted to verify aberrant splicing. The functional consequences of SF3B1 R625H were evaluated in vitro and in vivo . Critical makers of epithelial-mesenchymal transition and key components of relevant signaling pathways were detected by western blot, immunohistochemistry, and immunofluorescence, and were knocked down by siRNA-mediated silencing. Results Transcriptomic analysis of prolactinomas and heterozygous mutant cells revealed that the SF3B1 R625H allele led to different alterations in splicing properties, affecting different genes in different species. Consistently between rat cells and human tumor samples, mutant SF3B1 promoted aberrant splicing and the suppression of DLG1. Additionally, mutant SF3B1 with knockdown of DLG1 expression promoted cell migration, invasion, and epithelial-mesenchymal transition by activating the PI3K/Akt pathway. Conclusions Our findings elucidate a mechanism through which mutant SF3B1 promotes tumor progression and may provide a potent molecular therapeutic target for prolactinomas with the SF3B1 R625H mutation. Cancer Biology Oncology SF3B1 mutation Prolactinomas Alternative splicing DLG1 Invasion Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Background Pituitary adenomas account for between 10% and 20% of intracranial neoplasms, with prolactinomas being the most common subtype[ 1 , 2 ]. Except for hyperprolactinemia, a subset of prolactinomas are characterized by tumor invasion, resistance to conventional therapy, and high recurrence[ 3 ]. Therefore, deciphering the underlying pathogenesis of prolactinoma and determining efficient treatment targets are of great significance to the clinical management of prolactinoma patients. Alternative splicing is an essential step in the posttranscriptional regulation of gene expression that is a highly regulated and complex process that diversifies the proteome by creating multiple proteins from the same gene[ 4 , 5 ]. Mutations in splicing factor cause aberrant alternative splicing, which is a key molecular characteristic of tumorigenesis[ 6 ]. Splicing factor 3b subunit 1 (SF3B1) is a core component of the U2 small nuclear ribonucleoprotein complex (U2 snRNP), which is essential for pre-mRNA splicing[ 7 ]. Recurrent somatic mutations in SF3B1 have been reported in myelodysplastic syndrome (MDS), chronic lymphocytic leukemia (CLL), and some solid tumors such as uveal melanoma, breast carcinoma, pancreas adenocarcinoma, and others[ 8 – 13 ]. Mutant SF3B1 leads to splicing defects that result in aberrantly spliced transcripts[ 14 ]. Approximately half of the aberrantly spliced mRNAs are subject to nonsense-mediated decay (NMD), which causes gene and protein production to be downregulated[ 15 ]. Our previous study found a recurrent SF3B1 R625H mutation in prolactinoma, which was associated with poor progression-free survival (PFS) and higher levels of prolactin (PRL)[ 16 ]. However, the underlying molecular mechanisms of this SF3B1 mutation and its downstream cellular processes in prolactinoma remain unclear. Herein, we generated a heterozygous Sf3b1 R625H mutant rat cell line using CRISPR/Cas9 genome editing technology. RNA sequencing analysis revealed that while SF3B1 was conserved between rat and human, the aberrantly spliced mRNAs and affected genes differed significantly. Our results suggest that SF3B1 mutations in prolactinoma stimulate PI3K/AKT signaling by downregulating Discs large 1 (DLG1), which is induced by aberrant splicing and enhances the invasion and migration of tumor cells. Methods Cell culture and generation of Sf3b1 mutant R625H cells The GH3 rat pituitary cell line was purchased from the American Type Culture Collection (CLL-82.1; Manassas, VA, USA) and was cultured in Ham’s F12K medium with 2.5% fetal bovine serum (FBS) and 15% horse bovine serum (Gibco, Waltham, MA, USA). MCF7 and HEK293T cells were obtained from the National Infrastructure of Cell Line Resource (Beijing, China) and were cultured in DMEM (Gibco) supplemented with 10% FBS. PRL levels in cell culture supernatant were detected using an ELISA kit (BioVision, Milpitas, CA, USA) according to the manufacturer’s instructions. The CRISPR/Cas9 gene editing system was used to generate the Sf3b1 p.R625H (c.G1874A, c.A1875T) mutant GH3 cell line. A Sanger Centre CRISPR webtool ( http://www.sanger.ac.uk/htgt/wge/ ) was used to identify two small guide RNAs (sgRNAs). Cas9/sgRNA-mediated DNA double-strand break and homologous recombination contributed to generated the specific mutation. The 5′ sgRNA sequence was 5′-GGCAAGATTCCTTCCTCA-3′ and the 3′ sgRNA sequence was 5′-GGTTAAGAGTACTGTTGTC-3′. The sgRNA pairs were cloned into a wild-type spCas9 and sgRNA expression plasmid. The donor plasmid (CL-GJ-016-puro-ΔTK) consisted of a puromycin resistance cassette flanked by two loxP sites and ~ 1 kb homology arms at both ends of the c. 1874 G > A, c. 1875 A > T mutant exon 14 (Wuhan Genecreate Biological Engineering Co. Ltd., Wuhan, China) (Fig. 1 ). The plasmid was constructed according to previously published protocols and confirmed by sequencing. The Neon Transfection System (Thermo Fisher Scientific, Waltham, MA, USA) with pulse voltage of 1000 V for 40 ms was used to co-transfect GH3 cells on a 10-cm plate with 1 µg of CL-GJ-016-puro-TK and 3 µg of Cas9/sgRNA. The cells were screened with puromycin (0.5 µg/mL) starting 48-h post-transfection. Thereafter, single-cell cloning was performed in a 96-well plate. Sanger sequencing was performed using genomic DNA to confirm the specific mutations (DIA-UP Biotech, Beijing, China). Primers were designed to amplify exon 14 of the Sf3b1 gene (Table S1). Adenoviral constructs and primary cell culture of human prolactinomas The adenoviral constructs for mutant SF3B1 R625H and wild-type SFB31 were generated by BAC Biological Technology (Beijing, China), and primary cultures of human prolactinoma cells were prepared as described in our previous publication[ 16 ]. Tumor cells were infected with adenovirus at a multiplicity of infection of 100, and then harvested for quantitative (q)PCR 48-h later. RNA sequencing For GH3 cells, total RNA was extracted using the AllPrep® DNA/RNA Mini kit (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions. Sequencing libraries were generated using NEBNext® Ultra™ Directional RNA Library Prep Kit for Illumina® (NEB, Ipswich, MA, USA). The libraries were sequenced on an Illumina platform, and 150-bp paired-end reads were generated. Clean reads were mapped to the Rat Rnor_6.0 genome using Hisat2 (v2.0.5) to obtain read counts, fragments per kilobase per million mapped fragments (FPKM), and transcripts per million (TPM) for each sequenced gene. Differential alternative splicing and differential gene expression analyses Differentially expressed genes (DEGs) were identified using DESeq2 using P < 0.05 and absolute value of fold change ≥ 1.5. For alterative splicing analysis, rMATS (version 4.0.2)[ 17 ] was used to identify alternative splicing events by quantifying exon-exon junction spanning reads on annotated splice junctions present in the rat GENCODE Rnor_6.0 assembly. Differentially spliced mRNAs were defined as FDR 10% or < − 10%. Three mutant GH3 replicates and three wild-type replicates were compared. Data of two prolactinoma patients with the SF3B1 R625H mutation and two wild-type cases were selected for alternative splicing analysis, and two mutant and 13 wild-type cases from our previously reports were selected for differential expression analysis[ 16 ]. Basic clinic data of these patients are listed in Table S2. Functional enrichment analysis was performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) to predict the biological functions of DEGs and differentially spliced transcripts. The GO terms analyzed included biological process, cellular components, and molecular functions. A KEGG pathway with a P value < 0.05 was considered to be significant. Whole-genome sequencing and variant identification Genomic DNA was extracted from wild-type and mutant GH3 cells using the Blood & Cell Culture DNA Mini Kit (QIAGEN). DNA was sequenced using an Illumina HiSeq 2000, generating 150-bp paired reads. Reads were aligned to the Rat Rnor_6.0 genome using Burrows–Wheeler Aligner (BWA)[ 18 ]. SAMtools was used to generate BAM files[ 19 ]. Visual inspection of the Sf3b1 R625H mutation was performed using the Integrative Genomics Viewer[ 20 ]. Minigene assay A DNA fragment containing the DLG1 exon 22 genomic sequence with 151 bp of flanking intron 21 and 271 bp of flanking intron 22 was inserted between the KpnI and EcoRI restriction sites of the pcMINI vector to produce the DLG1 minigene construct. Sanger sequencing was performed to confirm the sequence of the inserted fragment. Briefly, the DLG1 minigene and an adenoviral vector were co-transfected into 293T cells. After 48 h, RNA was harvested for PCR analysis. PCR products were separated by 2% agarose gels and confirmed by Sanger sequencing. The primers used for the spliced products were: 5′-CTAGAGAACCCACTGCTTAC-3′ (forward) and 5′-TAGAAGGCACAGTCGAGG-3′ (reverse). Transfection RIBOBIO (Guangdong, China) synthesized small interfering (si)RNA duplexes; siRNA sequences of rat and human DLG1 are listed in Table S2. The Dlg1 overexpressed plasmid, pLV-hef1a-Puro-WPRE-CMV-Dlg1-3×FLAG, was constructed by Beijing Syngentech Co., Ltd. (Beijing, China). All transfections (siRNAs and overexpression plasmids) were performed using Lipofectamine® 3000 Transfection Reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s protocols. Cells were transfected with siRNA or plasmids for 48–72 h. and then collected for qPCR and western blot assays. Reverse transcription PCR and qPCR The RNeasy Mini Kit (QIAGEN) was used to extract total RNA and the High Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific) was used to generate cDNA, both according to the manufacturers’ instructions. For reverse transcription PCR, amplification was performed with I-5 High-Fidelity Master Mix (MCLAB, San Francisco, CA, USA) and specific primers. PCR products were separated on 1–3% agarose gels. All qPCR assays were performed using Power SYBR™ Green PCR Master Mix (Thermo Fisher Scientific) and were analyzed using QuantStudio 3 and 5 systems (Applied Biosystems, Waltham, MA, USA). The comparative Ct method was used to evaluate relative gene expression. The primers used in this study are listed in Table S3. Scanning electron microscopy (SEM) SEM was used to examine the external morphological alterations of Sf3b1 wild-type and mutant GH3 cells. Cells were collected and fixed with 2.5% glutaraldehyde (Solarbio, Beijing, China) at 4°C overnight for SEM preparation. After rinsing with PBS and sterile water, the cells were dehydrated using an ethanol gradient. The samples were coated with gold after critical point drying, and electron micrographs were taken on a Hitachi SU8020 SEM (Tokyo, Japan). Cell migration and wound-healing assays Cell migration culture dish inserts from Ibidi (Martinsried, Germany) were used to conduct wound-healing assays. After 24 h of transfection, cells were seeded into the chambers of the culture dish inserts. The inserts were removed on the next day, and fresh culture medium was added to each well. Scratches were photographed at different points in time using a Zeiss microscope (Oberkochen, Germany). Sf3b1 mutant and wild-type GH3 cells were seeded onto Imagelock 96-well plates (Essen Bioscience, Ann Arbor, MI, USA). The IncuCyte® Wound Maker (Essen Bioscience) was used to make uniform wounds in a monolayer of confluent cells. Phase contrast imaging was performed every 12 h for a total of 96 h. Images were analyzed using IncuCyte® S3 2018B-2019A software (Essen Bioscience), and data were analyzed using GraphPad Prism7 (GraphPad Software, Inc., La Jolla, CA, USA). Transwell assays Transwell plates and Matrigel-coated transwell plates (Corning-Costar, Corning, NY, USA) were used to determine cell migration and invasion capabilities, respectively. Briefly, HEK293T cells resuspended in serum-free medium were inoculate into the upper chamber, and DMEM with 10% FBS was placed in the bottom chamber. Following 24 h of culture, cells on the bottom surface of the chamber were stained with crystal violet and counted under a microscope (Zeiss). Western blot analysis Protein samples were separated on 8–10% Bis-Tris SDS-PAGE gels and transferred to polyvinylidene fluoride membranes (Merk, Kenilworth, NJ, USA). All primary antibodies (Table S4) were diluted in TBST containing 1% bovine serum albumin (BSA) and were incubated with the membranes overnight at 4°C. Immunoreactive bands were visualized by chemiluminescence. Immunohistochemistry Human prolactinoma tumor specimens were used to examine DLG1, E-cadherin, and Snail protein levels in SF3B1 wild-type and mutant tumor tissues. The primary antibodies used are summarized in Table S4. Immunohistochemistry was performed by the Leica Bond Polymer Refine Detection system (Leica Biosystems, Wetzlar, Germany). All slides were scanned into digital images, and expression was examined using Aperio AT2 (Leica Biosystems). Staining intensity was scored as 0 (negative), 1 (weak), 2 (moderate), or 3 (strong). The percentage of immunostaining was recorded, and H-scores were calculated using the formula: H-score = 1 × (% weakly stained cells) + 2 × (% moderately stained cells) + 3 × (% strongly stained cells)], ranging from 0 to 300. Rat prolactinoma model Rat pituitary tumors were induced by subcutaneously implanting 1-cm silastic capsules containing 10 mg of 17-β estradiol in 4-week-old female F344 rats. Prolactinomas were induced by 17β-estradiol for 5 weeks, as described in our previous report[ 21 ]. All experimental protocols were approved by the Animal Use and Care Committee of Beijing Tiantan Hospital. Five weeks later, the prolactinomas were validated via 7.0-T magnetic resonance imaging (MRI) before intra-pituitary injection. The rats were anaesthetized and then 1 µl of adenovirus vector control, wild-type SFB31 , or SF3B1 R625H was stereotactically injected into each bilateral tumor. A second MRI examination was performed 2 weeks later, and all rats were sacrificed to collect tumor tissues for further analysis. Immunofluorescence staining of rat tumor tissues First, rats were anaesthetized and heart-perfused with 4% formalin. Tumor tissues were then collected and immersed in 10% sucrose for 2 h, and then incubated in 30% sucrose overnight before freezing in OCT compound. Prior to immunofluorescence staining, frozen 5-µm-thick sections were fixed in ice-cold acetone. Tissue slides were blocked with goat serum and incubated with primary antibodies overnight at 4°C (summarized in Table S4). After washing with PBS, the slides were incubated in Alexa Fluor 488 and 594 secondary antibodies (Invitrogen) for 1 h at room temperature; DAPI was used to visualize nuclei. Phalloidin staining and confocal microscopy Cells were plated on confocal dishes coated with poly-L-lysine. After incubation at 37°C for 24 hours, the cells were fixed with 4% paraformaldehyde for 30 min, and then stained with 5 µg/mL Alexa Fluor 488-phalloidin (Invitrogen) to reveal filamentous actin (F-actin) in PBS for 40 min at 37°C. DAPI was used to visualize nuclei, and images were captured by confocal laser-scanning microscopy (Zeiss). Statistical analysis All statistical analyses were performed using R v3.4.1 ( https://www.r-project.org/ ) and Prism 7 (GraphPad Software, Inc.). All experiments were performed with at least three biological replicates, and all quantitative data represent mean ± standard deviation (SD). Statistical significance was determined by unpaired Student’s t test (two groups) or one-way ANOVA (multiple groups). P < 0.05 was considered statistically significant. Results Generation of Sf3b1 R625H-mutant cells Our previous study identified the somatic hotspot mutation SF3B1 R625H in 19.8% of prolactinomas [ 16 ]. To investigate the biological role of the SF3B1 R625H mutation on prolactinoma development and progression as well as similarities and differences between human and rats, we introduced the Sf3b1 R625H missense mutation into the rat pituitary cell line GH3 using the CRISPR/Cas9 gene editing system (Figure S1A). After isolating and expanding single cell clones, a heterozygous Sf3b1 -R625H mutation in GH3 cells was confirmed by Sanger sequencing (Fig. 1 A). Furthermore, whole-exome sequencing analyses of the mutant cell line also validated the heterozygous Sf3b1 mutation, p.R625H (NC_005108.4:g.61608130T > A, g.61608131C > T, NM_053426.1:c.1874G > A, c.1875A > T), and the frequency of the mutant allele obtained from read counts of whole-exome sequencing was 28.57% (18 A&T and 45 G&A, Figure S1B). Additionally, to determine whether heterozygous Sf3b1 -R625H mutation status affected PRL secretions, we performed ELISA to test PRL levels in culture media from wild-type and mutant cells. Consist with our previously reported results of mutant human prolactinoma[ 16 ], PRL levels in supernatant from Sf3b1 mutant cells was higher compared with wild-type cells (Fig. 1 B). Together, these data demonstrate successful construction of a rat cell line model harboring the Sf3b1 -R625H mutation to model human SF3B1 -mutant prolactinoma. Heterozygous Sf3b1 -R625H status altered the transcriptome of GH3 cells To explore transcriptome-level alterations in Sf3b1 -R625H mutant cells, we conducted RNA-seq analysis from wild-type and mutant cultured cells. Using P 1.5, we identified 3058 DEGs, with 1271 upregulated and 1787 downregulated genes (Figure S2A; Table S5). Clustering analysis showed distinctly different gene expression patterns between wild-type and Sf3b1 -mutant cells (Figure S2B). Moreover, compared with SF3B1 -wild-type patient samples ( n = 13), we identified 1062 downregulated and 599 upregulated genes in SF3B1 -mutant human prolactinomas ( n = 2) with P 1.5 (Figure S2C, 2D; Table S6). Comparing the two datasets, we found that only approximately 9% (285/3055) of the DEGs overlapped between rat GH3 cells and human prolactinoma samples (Fig. 1 C). Previous studies have described that the alternative splicing defects caused by cancer-associated SF3B1 mutations involve incorrect recognition of the 3′ splice site, producing aberrant transcripts [ 15 , 22 , 23 ]. To determine whether Sf3b1 -R625H is associated with mutation-specific alterations in pre-mRNA splicing, we performed percent spliced-in (PSI) analyses of alternative splicing using the rMATS tool[ 17 ]. At a False Discovery Rate (FDR) of 0.1, we obtained 143 alternative splicing events in Sf3b1 -R625H cells (Fig. 1 D, Table S7), with exon skipping being the most frequent event (Fig. 1 E). GO functional analysis of the alternatively spliced genes revealed enrichment of ontologies including lamellipodium organization, regulation of cell morphogenesis, and motor activity (Fig. 1 F, Table S8). We also performed qPCR to confirm several alternative splicing-associated changes (Figure S3A). Additionally, re-analysis of RNA-seq data generated from two SF3B1 -mutant and two SF3B1 -wild-type prolactinoma cases[ 16 ] uncovered 2287 alternative splicing events in SF3B1 -mutant cases, with a higher proportion of inclusion alternative 3′ splice site (ΔPSI > 10%) and inclusion retained intron (ΔPSI > 10%) events (Fig. 1 G; Figure S3B-C; Table S9-10). Similar to human data, Sf3b1 mutant GH3 cells had more inclusion ( n = 83, ΔPSI > 10%) than exclusion ( n = 60, ΔPSI < − 10%) alternative splicing events compared with wild-type cells (Fig. 1 E and G). Even so, there was minimal overlap between aberrantly spliced genes, and only approximately 15% (18/120) of the aberrantly spliced genes in rat Sf3b1 -mutant GH3 cells were common to human samples (Fig. 1 H). These data indicate that the limited conservation between rat and human intronic sequences caused the mutant SF3B1 to have different effects on alternative mRNA splicing in the two species. The R625H mutation induced epithelial-mesenchymal transition (EMT) phenotypes in tumor cells Our previous study found that most SF3B1 -mutant tumors exhibit invasive behavior, which is associated with poor PFS in prolactinoma[ 24 ]. Furthermore, GO enrichment analysis of the DEGs and alternative splicing events derived from GH3 cells revealed that cell adhesion molecules and lamellipodium organization were significantly altered in Sf3b1 -mutant GH3 cells (Figs. 2 A and 1 F, Tables S8 and S11). These results indicated that SFB31 -mutant tumor cells have invasive characteristics. To examine changes in the adhesive properties of mutant cells, SEM was performed. Compared with Sf3b1 wild-type cells, which showed a relatively rounded morphology, Sf3b1 mutant cells showed prominent lamellipodial extensions (Fig. 2 B). Furthermore, F-actin polymerization and filopodium formation were observed in Sf3b1 mutant cells (Fig. 2 B). Lamellipodia formation and dramatic reorganization of the actin cytoskeleton are involved in EMT, which is associated with tumor cell invasion and metastasis[ 25 , 26 ]. To explore whether Sf3b1 mutation promotes cell migration, wound-healing assays was conducted. The migration ability of Sf3b1 mutant cells was significantly enhanced compared with wild-type cells (Fig. 2 C). Moreover, the influence of Sf3b1 mutation on cell migration was also verified by the IncuCyte ZOOM® 96-well Scratch Wound cell migration assay (Fig. 2 D). To determine whether the Sf3b1 mutation modulated cell migration and invasion through the EMT pathway, we analyzed EMT maker expression. We observed a significant decrease in mesenchymal maker expression and an increase in epithelial maker expression in mutant cells. Immunoblotting of Sf3b1 mutant cells revealed downregulated expression of E-cadherin along with elevated expression of N-cadherin, Vimentin, and Snail (Fig. 2 E). Consistently, immunohistochemistry of consecutive sections from wild-type and mutant prolactinoma tissues supported these findings (Fig. 2 F). E-cadherin expression levels yielded an H-score of 146.7 ± 17.64 in wild-type samples and 53.3 ± 14.53 in mutant samples (P = 0.015). The H-scores of Snail were 20 ± 5.78 and 196.7 ± 20.28, respectively (P = 0.0011). To explore if SF3B1 -R625H promoted EMT of prolactinoma in vivo , F344 rat prolactinomas were induced by 6-weel 17β-estradiol treatment. Then the rat prolactinomas were stereotactically injected with wild-type and R625H SF3B1 encoding adenoviruses and negative controls. After 2 weeks, tumors were collected and immunofluorescence was used to examine EMT maker expression. Decreased levels of E-cadherin and increased levels of Vimentin were observed in the R625H groups (Fig. 2 G). Overall, these data demonstrated that mutant SF3B1 promoted EMT. To explore the molecular mechanisms through which the SF3B1 mutation induced EMT, we analyzed the intersection of aberrantly spliced genes and DEGs in rat GH3 cells and human prolactinoma samples (Fig. 2 H). From this intersection, DLG1 was the only identified gene (Fig. 2 H), and was significantly downregulated and mis-spliced in SF3B1 mutants (Tables S5–S7 and S9). Mutant SF3B1 caused aberrant splicing of DLG1 DLG1 is a member of the molecular scaffold protein family knows as membrane associated guanylate kinases (MAGUKs)[ 27 ]. DLG1 is a component of the conserved Scribble polarity complex, which is associated with the establishment and maintenance of cell polarity and is defined as a potential tumor suppressor[ 28 , 29 ]. To confirm that aberrant splicing of DLG1 was caused by SF3B1 mutation, RNA-seq data from human prolactinoma and GH3 cells was examined. Across human prolactinoma, aberrant splicing of DLG1 was induced by mutant SF3B1 through usage of a cryptic 3′ splice site (Fig. 3 A; Table S9), which was confirmed by qPCR in prolactinoma patient samples (Fig. 3 B). To verify the effect of mutant SF3B1 on DLG1 in vitro, we infected primary cultured pituitary tumor cells with adenovirus carrying the SF3B1 -R625H mutation, and the results showed that the cryptic DLG1 transcript was observed in the Ad-SF3B1-R625H group (Fig. 3 C). The same results were also observed in HEK293T cells (Fig. 3 D and S4A). To verify differences between the wild-type and mutant alleles, we performed TA cloning. Through reverse transcription PCR and sequencing clones, we confirmed that the Ad-SF3B1-R625H group included more abnormal PCR products than the Ad-SF3B1-WT group (Figure S4B). The expected fragment with a 16-bp extension of exon 22 was identified by Sanger sequencing of PCR products (Fig. 3 E). We then generated a minigene construct that contained sequences from the misspliced intron and flanking exon to explore alternative splicing (Fig. 3 F). After co-transfection into HEK293T cells, a complete spliced RNA containing exon 22 was observed in the SF3B1 -WT and control groups (Fig. 3 G, left and middle lanes); however, a larger transcript was observed in the SF3B1 -R625H group that retained an extra 16 bp (Fig. 3 G, right lane). The extra 16-bp fragment was confirmed by Sanger sequencing (Fig. 3 H). Additionally, mutually exclusive splicing of Dlg1 in Sf3b1 -mutant GH3 cells was confirmed by reverse transcription PCR (Fig. 3 I; Table S7). These data identified that mutant SF3B1 induced aberrant splicing of DLG1 . To verify whether the SF3B1 mutation altered DLG1 expression, RNA-seq data was examined and qPCR was performed in SF3B1 mutant samples. These data revealed that DLG1 was downregulated in SF3B1 -mutant human prolactinoma samples and GH3 cells (Fig. 4 A, B). Furthermore, immunohistochemistry and western blotting revealed that DLG1 protein expression was decreased in SF3B1 -mutant human samples and GH3 cells (Fig. 4 C and 2 E). The H-scores of DLG1 were 153 ± 17.64 in human SF3B1 -wild-type and 40 ± 20.82 in human SF3B1 -mutant samples (P = 0.0142). To further confirm the effect of the SF3B1 mutation on DLG1 expression, immunofluorescence was performed on tumors stereotactically injected with adenoviruses. Decreased DLG1 expression was observed in the Ad-SF3B1 R625H group (Fig. 4 D). Taken together, these data demonstrate that mutant SF3B1 induced aberrant splicing and reduced expression of DLG1 . Knocking down DLG1 promoted cell migration, invasion, and EMT To evaluate the role of DLG1 in tumor cell invasion and migration, we knocked down DLG1 by transfecting specific siRNAs. Both the mRNA and protein levels of DLG1 were efficiently depleted in GH3, 293T, and MCF7 cells (Fig. 5 A-B). Western blot results revealed that knocking down DLG1 promoted EMT by decreasing the expression of E-cadherin and increasing the expression of N-cadherin, Vimentin, and Snail (Fig. 5 B). Furthermore, restoring DLG1 expression in wild-type and mutant GH3 cells inhibited the EMT program (Fig. 2 E). Wound healing assays revealed that si-DLG1-transfected cells traversed wounds significantly faster than the si-NC-transfected cells (Fig. 5 C). Comparable migration differences were also observed in the transwell assay (Fig. 5 D). In transwell Matrigel invasion assays, significantly more si-DLG1-transfected cells migrated through the Matrigel layer (Fig. 5 D). Furthermore, phalloidin staining showed differences in F-actin organization at the wound margin. Compared with si-NC-transfected cells, si-DLG1-transfected cells showed polarized lamellipodia formation via F-actin (phalloidin) staining (Fig. 5 E). These results suggested that depleting DLG1 enhanced tumor cell migration and invasion, and induced EMT. The SFB31/DLG1 axis promoted tumor invasion via the PI3K/AKT pathway To determine the key biological mechanism through which mutant SF3B1 promoted prolactinoma progression, we performed a function enrichment analysis of the DEGs from human tumors and rat GH3 cells. This revealed that PI3K/AKT signaling was significantly altered in mutant SF3B1 samples compared with wild type samples (Fig. 6 A, Tables S12-13). We therefore detected the relevant effectors within this pathway through western blotting and found increased phosphorylation of AKT (T308) and GSK-3β in Sf3b1 -mutant GH3 cells (Fig. 6 B). To clarify the specific role of DLG1 in activation of the PI3K/AKT pathway in tumor cells, we manipulated DLG1 activity in GH3, 293T, and MCF7 cells. Knocking down DLG1 in cells by siRNA increased phosphorylation of AKT and GSK-3β (Fig. 6 C). Furthermore, restoration of DLG1 expression in wild-type and mutant GH3 cells showed reduced phosphorylation of AKT and GSK-3β (Fig. 6 B). Thus, these data strongly suggest that altering DLG1 levels due to alternative splicing induced by mutant SF3B1 promotes PI3K/AKT signaling in prolactinoma. Discussion SF3B1 is the most frequently mutated splicing gene in cancer, with three major hotspots: the codon positions K700, R625, and R622[ 30 ]. Among these mutated codons, K700E is the most common mutation in CLL, MDS, acute myeloid leukemia, and invasive breast carcinoma, while R625 is the predominant mutation in uveal melanoma and skin cutaneous melanoma[ 31 ]. Previous studies have explored how the SF3B1 K700E mutation contributes to aberrant splicing and downstream pro-tumorigenesis mechanisms in CLL and breast cancer using cell lines and mouse disease models[ 31 , 32 ]. To understand the molecular and phenotypic consequences of the SF3B1 R625H mutation on prolactinoma tumorigenesis, we generated a heterozygous Sf3b1 -R625H mutant rat pituitary cell line. We then demonstrated that heterozygosity for Sf3b1 -R625H in GH3 cells was sufficient to cause characteristic features similar to human prolactinoma, including higher PRL levels and invasive behaviors. Although phylogenetic analysis revealed that the SF3B1 R625H locus is highly conserved across species[ 16 ], only a few events were shared between mutant human prolactinoma and rat cells. These results may be due to the frequency of the mutant allele in heterozygous mutant cells and relatively poor interspecies conservation of intronic sequences between rats and humans. Moreover, our study found that the majority of aberrant splicing events in SF3B1 mutant human prolactinoma were retained intron followed by alternative 3′ splice site. This is different from previous observations in SF3B1 -mutant cancers[ 15 , 33 ], which suggested a higher proportion of alternative 3′ splice site events, but were similar to a recent report in SF3B1 -mutant MDS[ 14 ]. Considering the limited number SF3B1 -mutant human prolactinoma tumor samples, further studies with larger cohorts are needed to confirm these results. In particular, we consistently found that the DLG1 transcript was shared between rat cells and human prolactinoma samples, with mRNA downregulation and differential missplicing. Interestingly, mutant SF3B1 induced the inclusion of 16 intronic nucleotides in exon 22 by using an upstream cryptic 3′ splice site in human sequences (Fig. 3 A), which may lead to NMD and reduced DLG1 expression, because the shifted reading frame produces a premature termination codon in this alternative spliceform[ 32 ]. Despite mutant SF3B1 inducing different splicing events in rats and humans, DLG1 mRNA and protein were downregulated in both models. Our data demonstrated that downregulating DLG1 expression promote tumor cell migration and invasion. As previously reported, DLG1 was involved in mammalian tumorigenesis and was localized to adherens junctions in epithelial cells[ 34 ]. In mammals, DLG1 appears to act as a tumor suppressor by affecting other polarity complexes, Myosin II activity, the actin cytoskeleton, and/or other signaling pathways[ 35 – 37 ]. Our studies further indicated that knocking down DLG1 induced EMT in tumor cells. EMT is well known to be related to the loss of cell polarity, specifically the loss of adherens junctions and enables the development of an invasive phenotype[ 38 ]. Consist with our results, the study by Sugihara, T. et al. found that loss of DLG1 was associated with poor prognosis in endometrial cancer and that knocking down DLG1 accelerated tumor migration and invasion in vitro [ 29 ]. Previous studies have revealed that the AKT/GSK-3β/Snail axis is critical for the induction and maintenance of EMT [ 39 , 40 ]. In this study, we show that mutant SF3B1 and decreased DLG1 expression activate the PI3K/AKT pathway, which results in tumor cells undergoing EMT and gaining migration and invasive properties (Fig. 7 ). It is known that p-AKT can suppress GSK-3β activity (through phosphorylation of Ser 9) and stabilize Snail[ 41 ]. Furthermore, DLG1 has been demonstrated to modulate the PI3K/Akt pathway to precisely regulate regulatory T cell activity[ 42 ]. Consistent with this, we found that low DLG1 expression increased p-AKT, p-G3K-3β, and Snail expression levels, which could be rescued by overexpressing DLG1. Future work comprehensively defining the how DLG1 activates the PI3K/Akt signaling pathway in tumor cell may therefore be very important. Conclusions In summary, we confirmed the molecular mechanism through which the SF3B1 -R625H mutation actives the PI3K/Akt pathway in prolactinoma by aberrant splicing of DLG1 , which promotes tumor invasion progression. Our findings provide the rationale for investigating new therapeutic strategies in patients with SF3B1 -mutant prolactinoma. Abbreviations MDS Myelodysplastic syndrome CLL Chronic lymphocytic leukemia NMD Nonsense-mediated decay PFS Progression-free survival FBS Fetal bovine serum FPKM Fragments per kilobase per million mapped fragments TPM Transcripts per million DEGs Differentially expressed genes GO Gene Ontology KEGG Kyoto Encyclopedia of Genes and Genomes BSA Bovine serum albumin MRI Magnetic resonance imaging MAGUKs Membrane associated guanylate kinases Declarations Ethics approval and consent to participate The present study was approved by The Ethics Committees of Beijing Tiantan Hospital. All subjects provided written informed consent. The animal experiments were approved by the Animal Use and Care Committee of Beijing Tiantan Hospital. Consent for publication Not applicable. Availability of data and materials The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. Competing interests The authors declare that they have no competing interests. Funding This study was supported by the National Natural Science Foundation of China (Grant codes: 81672495, 81771489, 82072804, 82071559, and 82071558). Authors' contributions WX and YZ worked on the conception and designed the research. CL and DW were involved in the collection and analysis of patients’ clinical data. JG, QF and YL performed the experiments. JG was dedicated to data analysis, interpretation, and drafting. All authors read and approved the final manuscript. Acknowledgements Not applicable. References Asa SL, Ezzat S: The cytogenesis and pathogenesis of pituitary adenomas . Endocr Rev 1998, 19 (6):798-827. Daly AF, Rixhon M, Adam C, Dempegioti A, Tichomirowa MA, Beckers A: High prevalence of pituitary adenomas: a cross-sectional study in the province of Liege, Belgium . J Clin Endocrinol Metab 2006, 91 (12):4769-4775. Vroonen L, Daly AF, Beckers A: Epidemiology and Management Challenges in Prolactinomas . Neuroendocrinology 2019, 109 (1):20-27. Nilsen TW, Graveley BR: Expansion of the eukaryotic proteome by alternative splicing . Nature 2010, 463 (7280):457-463. Lee Y, Rio DC: Mechanisms and Regulation of Alternative Pre-mRNA Splicing . Annu Rev Biochem 2015, 84 :291-323. Sveen A, Kilpinen S, Ruusulehto A, Lothe RA, Skotheim RI: Aberrant RNA splicing in cancer; expression changes and driver mutations of splicing factor genes . Oncogene 2016, 35 (19):2413-2427. Kramer A: The structure and function of proteins involved in mammalian pre-mRNA splicing . Annu Rev Biochem 1996, 65 :367-409. Papaemmanuil E, Cazzola M, Boultwood J, Malcovati L, Vyas P, Bowen D, Pellagatti A, Wainscoat JS, Hellstrom-Lindberg E, Gambacorti-Passerini C et al : Somatic SF3B1 mutation in myelodysplasia with ring sideroblasts . N Engl J Med 2011, 365 (15):1384-1395. Wang L, Lawrence MS, Wan Y, Stojanov P, Sougnez C, Stevenson K, Werner L, Sivachenko A, DeLuca DS, Zhang L et al : SF3B1 and other novel cancer genes in chronic lymphocytic leukemia . N Engl J Med 2011, 365 (26):2497-2506. Banerji S, Cibulskis K, Rangel-Escareno C, Brown KK, Carter SL, Frederick AM, Lawrence MS, Sivachenko AY, Sougnez C, Zou L et al : Sequence analysis of mutations and translocations across breast cancer subtypes . Nature 2012, 486 (7403):405-409. Biankin AV, Waddell N, Kassahn KS, Gingras MC, Muthuswamy LB, Johns AL, Miller DK, Wilson PJ, Patch AM, Wu J et al : Pancreatic cancer genomes reveal aberrations in axon guidance pathway genes . Nature 2012, 491 (7424):399-405. Furney SJ, Pedersen M, Gentien D, Dumont AG, Rapinat A, Desjardins L, Turajlic S, Piperno-Neumann S, de la Grange P, Roman-Roman S et al : SF3B1 mutations are associated with alternative splicing in uveal melanoma . Cancer Discov 2013, 3 (10):1122-1129. Armenia J, Wankowicz SAM, Liu D, Gao J, Kundra R, Reznik E, Chatila WK, Chakravarty D, Han GC, Coleman I et al : The long tail of oncogenic drivers in prostate cancer . Nat Genet 2018, 50 (5):645-651. Pellagatti A, Armstrong RN, Steeples V, Sharma E, Repapi E, Singh S, Sanchi A, Radujkovic A, Horn P, Dolatshad H et al : Impact of spliceosome mutations on RNA splicing in myelodysplasia: dysregulated genes/pathways and clinical associations . Blood 2018, 132 (12):1225-1240. Darman RB, Seiler M, Agrawal AA, Lim KH, Peng S, Aird D, Bailey SL, Bhavsar EB, Chan B, Colla S et al : Cancer-Associated SF3B1 Hotspot Mutations Induce Cryptic 3' Splice Site Selection through Use of a Different Branch Point . Cell reports 2015, 13 (5):1033-1045. Li C, Xie W, Rosenblum JS, Zhou J, Guo J, Miao Y, Shen Y, Wang H, Gong L, Li M et al : Somatic SF3B1 hotspot mutation in prolactinomas . Nat Commun 2020, 11 (1):2506. Shen S, Park JW, Lu ZX, Lin L, Henry MD, Wu YN, Zhou Q, Xing Y: rMATS: robust and flexible detection of differential alternative splicing from replicate RNA-Seq data . Proc Natl Acad Sci U S A 2014, 111 (51):E5593-5601. Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform . Bioinformatics 2009, 25 (14):1754-1760. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Genome Project Data Processing S: The Sequence Alignment/Map format and SAMtools . Bioinformatics 2009, 25 (16):2078-2079. Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP: Integrative genomics viewer . Nat Biotechnol 2011, 29 (1):24-26. Lei C, Jing G, Jichao W, Xiaohui L, Fang Q, Hua G, Yazhou M, Zhang Y: MiR-137's Tumor Suppression on Prolactinomas by Targeting MITF and Modulating Wnt Signaling Pathway . J Clin Endocrinol Metab 2019, 104 (12):6391-6402. Wang L, Brooks AN, Fan J, Wan Y, Gambe R, Li S, Hergert S, Yin S, Freeman SS, Levin JZ et al : Transcriptomic Characterization of SF3B1 Mutation Reveals Its Pleiotropic Effects in Chronic Lymphocytic Leukemia . Cancer cell 2016, 30 (5):750-763. Shiozawa Y, Malcovati L, Galli A, Sato-Otsubo A, Kataoka K, Sato Y, Watatani Y, Suzuki H, Yoshizato T, Yoshida K et al : Aberrant splicing and defective mRNA production induced by somatic spliceosome mutations in myelodysplasia . Nat Commun 2018, 9 (1):3649. Kim K, Park YW, Kim D, Ahn SS, Moon JH, Kim EH, Lee EJ, Ku CR: Biochemical Remission after Cabergoline Withdrawal in Hyperprolactinemic Patients with Visible Remnant Pituitary Adenoma . J Clin Endocrinol Metab 2021, 106 (2):e615-e624. Yilmaz M, Christofori G: EMT, the cytoskeleton, and cancer cell invasion . Cancer Metastasis Rev 2009, 28 (1-2):15-33. Aiello NM, Maddipati R, Norgard RJ, Balli D, Li J, Yuan S, Yamazoe T, Black T, Sahmoud A, Furth EE et al : EMT Subtype Influences Epithelial Plasticity and Mode of Cell Migration . Dev Cell 2018, 45 (6):681-695 e684. Ivanova S, Gregorc U, Vidergar N, Javier R, Bredt DS, Vandenabeele P, Pardo J, Simon MM, Turk V, Banks L et al : MAGUKs, scaffolding proteins at cell junctions, are substrates of different proteases during apoptosis . Cell Death Dis 2011, 2 :e116. Facciuto F, Cavatorta AL, Valdano MB, Marziali F, Gardiol D: Differential expression of PDZ domain-containing proteins in human diseases - challenging topics and novel issues . FEBS J 2012, 279 (19):3538-3548. Sugihara T, Nakagawa S, Sasajima Y, Ichinose T, Hiraike H, Kondo F, Uozaki H, Fukusato T, Ayabe T: Loss of the cell polarity determinant human Discs-large is a novel molecular marker of nodal involvement and poor prognosis in endometrial cancer . Br J Cancer 2016, 114 (9):1012-1018. Alsafadi S, Houy A, Battistella A, Popova T, Wassef M, Henry E, Tirode F, Constantinou A, Piperno-Neumann S, Roman-Roman S et al : Cancer-associated SF3B1 mutations affect alternative splicing by promoting alternative branchpoint usage . Nat Commun 2016, 7 :10615. Liu Z, Yoshimi A, Wang J, Cho H, Chun-Wei Lee S, Ki M, Bitner L, Chu T, Shah H, Liu B et al : Mutations in the RNA Splicing Factor SF3B1 Promote Tumorigenesis through MYC Stabilization . Cancer Discov 2020, 10 (6):806-821. Liu B, Liu Z, Chen S, Ki M, Erickson C, Reis-Filho JS, Durham BH, Chang Q, de Stanchina E, Sun Y et al : Mutant SF3B1 promotes AKT- and NF-kappaB-driven mammary tumorigenesis . J Clin Invest 2021, 131 (1). DeBoever C, Ghia EM, Shepard PJ, Rassenti L, Barrett CL, Jepsen K, Jamieson CH, Carson D, Kipps TJ, Frazer KA: Transcriptome sequencing reveals potential mechanism of cryptic 3' splice site selection in SF3B1-mutated cancers . PLoS Comput Biol 2015, 11 (3):e1004105. Humbert PO, Grzeschik NA, Brumby AM, Galea R, Elsum I, Richardson HE: Control of tumourigenesis by the Scribble/Dlg/Lgl polarity module . Oncogene 2008, 27 (55):6888-6907. Peng CY, Manning L, Albertson R, Doe CQ: The tumour-suppressor genes lgl and dlg regulate basal protein targeting in Drosophila neuroblasts . Nature 2000, 408 (6812):596-600. Bilder D, Schober M, Perrimon N: Integrated activity of PDZ protein complexes regulates epithelial polarity . Nat Cell Biol 2003, 5 (1):53-58. Hawkins ED, Russell SM: Upsides and downsides to polarity and asymmetric cell division in leukemia . Oncogene 2008, 27 (55):7003-7017. Lamouille S, Xu J, Derynck R: Molecular mechanisms of epithelial-mesenchymal transition . Nat Rev Mol Cell Biol 2014, 15 (3):178-196. Lan Y, Han J, Wang Y, Wang J, Yang G, Li K, Song R, Zheng T, Liang Y, Pan S et al : STK17B promotes carcinogenesis and metastasis via AKT/GSK-3beta/Snail signaling in hepatocellular carcinoma . Cell Death Dis 2018, 9 (2):236. Stegeman H, Span PN, Kaanders JH, Bussink J: Improving chemoradiation efficacy by PI3-K/AKT inhibition . Cancer Treat Rev 2014, 40 (10):1182-1191. Zhou BP, Deng J, Xia W, Xu J, Li YM, Gunduz M, Hung MC: Dual regulation of Snail by GSK-3beta-mediated phosphorylation in control of epithelial-mesenchymal transition . Nat Cell Biol 2004, 6 (10):931-940. Zanin-Zhorov A, Lin J, Scher J, Kumari S, Blair D, Hippen KL, Blazar BR, Abramson SB, Lafaille JJ, Dustin ML: Scaffold protein Disc large homolog 1 is required for T-cell receptor-induced activation of regulatory T-cell function . Proc Natl Acad Sci U S A 2012, 109 (5):1625-1630. Supplementary Files FigureS1.pdf Figure S1. Generation of heterozygous Sf3b1 R625H mutant GH3 cells. (A) Schematic of the Sf3b1 R625H CRISPR/Cas9-mediated strategy. (B) Visualization of the engineered mutation in the Integrative Genomic Viewer (IGV) browser: whole-exome sequencing reads overlapping the heterozygous mutation in the Sf3b1 gene (c.1874G>A, c.1875A>T, p.R625H, frequency: 28.57%). FigureS2.pdf Figure S2. SF3B1-R625H alters the transcriptome in GH3 cells and human prolactinoma. (A) Volcano plot comparing mRNA expression between Sf3b1 mutant vs. wild-type GH3 cells. (B) Heatmap of unsupervised clustering of differentially expressed genes in Sf3b1 mutant vs. wild-type GH3 cells.. (C) Volcano plot comparing mRNA expression between SF3B1 mutant vs. wild-type human prolactinoma.. (B) Heatmap of unsupervised clustering of differentially expressed genes in SF3B1 mutant vs. wild-type human prolactinoma. FigureS3.pdf Figure S3. Alternative splicing in Sf3b1 mutant GH3 cells. (A-N) PCR products of the 13 genes (not including Dlg1) – Pkd1l2, Myo15a, Tgif1, Rest, Mroh2a, Pipox, Emc1, Mapk15, Ccdc136, Tfb2m, Parp16, Parp4, and Cd44 were amplified from Sf3b1-mutant and WT GH3 cells. The cryptic and canonical transcripts were indicated by the arrow. Gapdh was used as a loading control. FigureS4.pdf Figure S4. AS events in SF3B1 mutant human prolactinoma. (A) Volcano plot showing the differential alternative splicing events in SF3B1 mutant vs. wild-type human prolactinoma. Top 20 significant mis-spliced genes are indicated. (B) GO analysis of mis-spliced genes in SF3B1 mutant vs. wild-type human prolactinoma; the top five ranked terms are shown. FigureS5.pdf Figure S5. Mutant SF3B1 caused aberrant splicing of DLG1 in 293T and GH3 cells. (A) Schematic representation of the validated aberrant splicing of DLG1 gene fragments by PCR. Schematic a represents the canonical DLG1 transcript and schematic b represents the aberrant DLG1 transcript. (B) The number of abnormal variant clones in the Ad-SF3B1-WT (2/48) and Ad-SF3B1-MT (7/43) groups, out of a total of 50 TA clones. (C) Sashimi plots of aberrant Dlg1 splicing event in GH3 cells with or without Sf3b1 mutations. SupplementalTables113.xlsx Cite Share Download PDF Status: Published Journal Publication published 17 Jan, 2022 Read the published version in Journal of Experimental & Clinical Cancer Research → Version 1 posted Editorial decision: Major revision 01 Nov, 2021 Review # 1 received at journal 31 Oct, 2021 Review # 2 received at journal 07 Oct, 2021 Reviewer # 2 agreed at journal 02 Oct, 2021 Reviewer # 1 agreed at journal 01 Oct, 2021 Reviewers invited by journal 24 Sep, 2021 Editor assigned by journal 23 Sep, 2021 Submission checks completed at journal 23 Sep, 2021 Editor invited by journal 23 Sep, 2021 First submitted to journal 23 Sep, 2021 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-934219","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":54094275,"identity":"8994b1d5-989f-4bdb-bfc1-8473f733194c","order_by":0,"name":"Jing Guo","email":"","orcid":"","institution":"Beijing Neurosurgical Institute","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Jing","middleName":"","lastName":"Guo","suffix":""},{"id":54094276,"identity":"1609c187-c7c4-46b2-85da-1e6390a21259","order_by":1,"name":"Qiuyue Fang","email":"","orcid":"","institution":"Beijing Neurosurgical Institute","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Qiuyue","middleName":"","lastName":"Fang","suffix":""},{"id":54094277,"identity":"66a0958a-962c-456e-a567-84c07d975424","order_by":2,"name":"Yulou Liu","email":"","orcid":"","institution":"Beijing Neurosurgical Institute","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Yulou","middleName":"","lastName":"Liu","suffix":""},{"id":54094278,"identity":"323a5936-23d8-4945-af21-8f92072182bd","order_by":3,"name":"Dawei Wang","email":"","orcid":"","institution":"Beijing Neurosurgical Institute","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Dawei","middleName":"","lastName":"Wang","suffix":""},{"id":54094279,"identity":"036c783a-f80f-40cb-9624-ae27b2834ee2","order_by":4,"name":"Chuzhong Li","email":"","orcid":"","institution":"Beijing Tiantan Hospital","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Chuzhong","middleName":"","lastName":"Li","suffix":""},{"id":54094280,"identity":"00dec2a3-d442-40ae-b86a-904632b136e4","order_by":5,"name":"Weiyan Xie","email":"","orcid":"","institution":"Beijing Neurosurgical Institute","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Weiyan","middleName":"","lastName":"Xie","suffix":""},{"id":54094281,"identity":"6796c254-1cc8-4a52-8ad2-914e8e855cc6","order_by":6,"name":"Yazhuo Zhang","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA8klEQVRIiWNgGAWjYDACCShtwMDA+ADBJlILswGCTaQWNiQ2HiA/u/nZw69tdvbmEsnHKn7UbEtsYG/eJsFQcwenFsY5x8yNZduSmS1npKXd7Dl2O7GB51iZBMOxZzi1MEskmElLtjGzGdzIMbvN2ADUIpFjJsHYcBinFjaJ9G9ALfU8IC3FYC3yb/Br4QGaKfmx7bAESAszxBYe/FokJHLKpBnOHTcwOPMsWRLoF+M2nrRii4RjuLXIz0jfJvmjrNre4HjywQ8/am7L9rMf3njjQw1uLeAg4GVD9h2ISMCrARjQP/4QUDEKRsEoGAUjGwAAg2pRxgTFOA4AAAAASUVORK5CYII=","orcid":"","institution":"Beijing Neurosurgical Institute","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Yazhuo","middleName":"","lastName":"Zhang","suffix":""}],"badges":[],"createdAt":"2021-09-24 10:34:03","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-934219/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-934219/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1186/s13046-022-02245-0","type":"published","date":"2022-01-17T09:21:26+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":14083230,"identity":"64fbcf72-034c-42e0-bf8d-256994c8e00e","added_by":"auto","created_at":"2021-09-28 19:24:04","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":461699,"visible":true,"origin":"","legend":"Generation of heterozygous Sf3b1-R625H mutant GH3 cells. \n(A) Sequencing chromatogram demonstrating the heterozygous Sf3b1-R625H mutation. The homozygous nucleotide in Sf3b1 wild-type GH3 cells and heterozygous nucleotides in Sf3b1 mutant GH3 cells are shown. (B) PRL secretion in Sf3b1 wild-type and mutant GH3 cells are shown. Data are represented as mean ± SD. (P-values by two-way ANOVA, ****P\u003c0.0001). (C) Venn diagram showing the overlapping DEGs identified in SF3B1 mutant human prolactinoma and GH3 cells. (D) Volcano plot showing the differential alternative splicing events in Sf3b1 mutant vs. wild-type GH3 cells. Significant alternative splicing events identified at FDR \u003c0.05 and |∆PSI| \u003e0.1. PSI: Percent Splice In. Top 30 significant mis-spliced genes are indicated. (E) Number and type of mis-splicing events in Sf3b1 mutant vs. wild-type GH3 cells are shown. Inclusion: events with ΔPSI\u003e10%; exclusion: events with ΔPSI\u003c−10% (F) GO analysis of mis-spliced genes in Sf3b1 mutant vs. wild-type GH3 cells; the top five ranked terms are shown. (G) Number and type of mis-splicing events in SF3B1 mutant vs. wild-type human prolactinoma. (H) Venn diagram showing the overlap of mis-spliced genes identified in SF3B1 mutant human prolactinoma and GH3 cells.\n","description":"","filename":"Figure1.png","url":"https://assets-eu.researchsquare.com/files/rs-934219/v1/b77287bf8b66908331dab9f3.png"},{"id":14083232,"identity":"cb92197c-b832-4b6c-962f-e53905f05716","added_by":"auto","created_at":"2021-09-28 19:24:04","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":1619202,"visible":true,"origin":"","legend":"Mutant SF3B1 promoted EMT phenotypes in prolactinoma. \n(A) GO functional enrichment analysis of DEGs in Sf3b1 mutant vs. wild-type GH3 cells. (B) Scanning electron microscopy (SEM; upper) and confocal images of the actin cytoskeleton stained with rhodamine–phalloidin (lower) of Sf3b1 wild-type and mutant GH3 cells; scale bar: 10 µm. (C) Scratch wound healing assays were performed using Sf3b1 wild-type and mutant GH3 cells; scale bar: 200 µm. Relative scratch widths are shown over time (n=3; mean ± SD shown below; P values by Student’s test). (D) Kinetics of wound confluence over 96 h were analyzed using IncuCyte 2018B software (Essen Bioscience). The results are expressed as mean ± SD, n=3. P values by two-way ANOVA; ****P\u003c0.0001. (E) DLG1, E-cadherin, N-cadherin, Vimentin, and Snail levels in Sf3b1 wild-type and mutant GH3 cells transduced with or without Dlg1-FLAG cDNA. (G) Immunofluorescence staining of E-cadherin and Vimentin in rat prolactinoma tumors stereotactically injected with Ad-null, Ad-SF3B1WT, and Ad-SF3B1R625H; scale bar: 50 µm. (F) Representative immunohistochemical staining for E-cadherin and Snail in SF3B1 mutant and wild-type human prolactinoma tumors; scale bar: 50 µm. (H) Venn diagram of the numbers of differentially spliced (q\u003c0.05; t test) and differentially expressed genes (q\u003c0.05; DESeq2) in SF3B1 mutant human prolactinoma and GH3 samples.\n","description":"","filename":"Figure2.png","url":"https://assets-eu.researchsquare.com/files/rs-934219/v1/a276f67282f33f4bfb0ead0b.png"},{"id":14083105,"identity":"593126b3-c548-46ba-93cf-660a25a08ce7","added_by":"auto","created_at":"2021-09-28 19:21:04","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":662920,"visible":true,"origin":"","legend":"Mutant SF3B1 caused aberrant splicing of DLG1. \n(A) Sashimi plots of alternative 3′ splice sites of DLG1 in representative human prolactinoma with or without SF3B1 mutations. (B) Reverse transcription PCR of DLG1 splicing events in human prolactinoma samples with or without SF3B1 mutation. The 3′ splice sites used are shown in the schematic. (C) Reverse transcription PCR results of aberrant DLG1 splicing in primary human prolactinoma cells infected with adenovirus carrying SF3B1-WT and SF3B1-R625H vectors. (D) Reverse transcription PCR of DLG1 splicing in 293T cells infected with adenovirus carrying SF3B1-WT and SF3B1-R625H vectors. Band a: canonical DLG1 transcript; Band b: aberrant DLG1 transcript. (E) The sequences of gel-purified fragments showing the canonical DLG1 transcript (a) and aberrant transcript (b). (F) Schematic diagram of the minigene construct and schematic diagram of Sanger sequencing of reverse transcription PCR products. (G) Reverse transcription PCR results of DLG1 splicing in 293T cells co-transduced with the minigene reporter and adenovirus carrying SF3B1-WT and SF3B1-R625H vectors. (H) The sequence of gel-purified fragments. The canonical DLG1 transcript (a) and aberrant DLG1 transcript (b) were observed in the SF3B1-WT and SF3B1-R625 groups, respectively. (I) Reverse transcription PCR results of Dlg1 splicing in Sf3b1 mutant GH3 cells. The mutually exclusive splicing type is shown in the schematic.\n","description":"","filename":"Figure3.png","url":"https://assets-eu.researchsquare.com/files/rs-934219/v1/5cda1d55f2292e9bf1d09e60.png"},{"id":14083107,"identity":"e2d5cd5d-91c7-48fa-829d-11025d38414e","added_by":"auto","created_at":"2021-09-28 19:21:04","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":501471,"visible":true,"origin":"","legend":"Mutant SF3B1 altered DLG1 expression. \n(A) DLG1 mRNA expression (RNA-seq data and quantitative PCR) in SF3B1 mutant (n=2) vs. wild-type (n=13) human prolactinoma tumors. (B) DLG1 mRNA expression in Sf3b1 mutant (n=3) versus wild-type (n=3) GH3 cells. TPM: Transcripts per million kilobase. The results are expressed as mean ± SD (Student’s t test, *P\u003c0.05, ****P\u003c0.0001). (C) Representative immunohistochemical staining for DLG1 in SF3B1 mutant and wild-type human prolactinoma tumors; scale bar: 50 µm. (D) Immunofluorescence staining for DLG1 in rat prolactinoma tumors stereotactically injected with Ad-null, Ad-SF3B1WT, and Ad-SF3B1R625H, respectively; scale bar: 50 µm. \n","description":"","filename":"Figure4.png","url":"https://assets-eu.researchsquare.com/files/rs-934219/v1/3cabe0401c8d976162fe8f07.png"},{"id":14083231,"identity":"a6a63214-fd8e-4be1-871d-841f8dd770ff","added_by":"auto","created_at":"2021-09-28 19:24:04","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":1156001,"visible":true,"origin":"","legend":"Knocking down DLG1 promoted cell migration, invasion, and EMT. \n(A) Quantitative PCR analysis of DLG1 knockdown efficiency in GH3, MCF-7, and 293T cells. Results are expressed as mean ± SD; P values by one-way ANOVA compared with control (si-NC). (B) DLG1 knockdown significantly increased N-cadherin, Vimentin, and Snail expression but decreased E-cadherin expression in GH3, MCF-7, and 293T cells. (C) Scratch assay of 293T and MCF7 cells with or without DLG1 siRNA; scale bar: 200 µm. Relative scratch width over time shown below (n=3; mean ± SD shown below; P values by two-way ANOVA). (D) Transwell assay of 293T cells transfected with negative control (si-NC) and si-DLG1; scale bar: 200 µm. Data are presented as the mean ± SD from three independent experiments. P values by one-way ANOVA; **P\u003c0.01; ***P\u003c0.001; ****P\u003c0.0001, relative to control (si-NC). (E) Confocal immunofluorescence images in wound healing migration assays of MCF7 cells. After transfection with siRNAs for 48 h, MCF7 cells were scratched to allow cell migration. Cells were fixed after migrating for 12 h, and then stained with phalloidin (green) and DAPI (blue) to reveal F actin and nuclei, respectively; scale bar: 50 µm.\n","description":"","filename":"Figure5.png","url":"https://assets-eu.researchsquare.com/files/rs-934219/v1/3c6651269c051f9757658d07.png"},{"id":14083228,"identity":"a97c4f00-8444-4da1-b5ae-5b63bcf48fbd","added_by":"auto","created_at":"2021-09-28 19:24:04","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":529820,"visible":true,"origin":"","legend":"The SFB31/DLG1 axis promoted tumor invasion via the PI3K/AKT signaling pathway\n(A) KEGG enrichment analysis of DEGs in SF3B1 mutant vs. wild-type human prolactinoma and GH3 cells. (B) Phosphorylated and total AKT and GSK-3β levels in Sf3b1 WT and mutant GH3 cells transduced with or without Dlg1-FLAG cDNA. (C) Phosphorylated and total AKT and GSK-3β levels in GH3, 293T, and MCF7 cells transduced with or without DLG1 siRNA.\n","description":"","filename":"Figure6.png","url":"https://assets-eu.researchsquare.com/files/rs-934219/v1/245cc84d9deab512f56b66e2.png"},{"id":14083110,"identity":"a728efb6-5fc8-44bf-be4b-2c6cf88c9eb9","added_by":"auto","created_at":"2021-09-28 19:21:04","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":353618,"visible":true,"origin":"","legend":"Schematic representation showing the proposed mechanisms through which the SF3B1 mutation promotes the prolactinoma progression.","description":"","filename":"Figure7.png","url":"https://assets-eu.researchsquare.com/files/rs-934219/v1/03a32115644ecc63dbfade91.png"},{"id":17371911,"identity":"d87d4d28-8db7-4fa2-822b-8c76a8b4ae85","added_by":"auto","created_at":"2022-01-17 09:21:30","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":5441070,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-934219/v1/1ee44ee7-8e8b-4fea-8726-c19351451408.pdf"},{"id":14083233,"identity":"347dc5de-b254-46c5-9f18-c88634bd79e2","added_by":"auto","created_at":"2021-09-28 19:24:04","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":616446,"visible":true,"origin":"","legend":"Figure S1. Generation of heterozygous Sf3b1 R625H mutant GH3 cells.\n(A) Schematic of the Sf3b1 R625H CRISPR/Cas9-mediated strategy. (B) Visualization of the engineered mutation in the Integrative Genomic Viewer (IGV) browser: whole-exome sequencing reads overlapping the heterozygous mutation in the Sf3b1 gene (c.1874G\u003eA, c.1875A\u003eT, p.R625H, frequency: 28.57%).\n","description":"","filename":"FigureS1.pdf","url":"https://assets-eu.researchsquare.com/files/rs-934219/v1/67f325c0898ee3cdeaa860d7.pdf"},{"id":14083116,"identity":"ff6e5d0e-c54b-4e0f-9b44-bfcbd27f894f","added_by":"auto","created_at":"2021-09-28 19:21:04","extension":"pdf","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":10967268,"visible":true,"origin":"","legend":"Figure S2. SF3B1-R625H alters the transcriptome in GH3 cells and human prolactinoma.\n(A) Volcano plot comparing mRNA expression between Sf3b1 mutant vs. wild-type GH3 cells. (B) Heatmap of unsupervised clustering of differentially expressed genes in Sf3b1 mutant vs. wild-type GH3 cells.. (C) Volcano plot comparing mRNA expression between SF3B1 mutant vs. wild-type human prolactinoma.. (B) Heatmap of unsupervised clustering of differentially expressed genes in SF3B1 mutant vs. wild-type human prolactinoma.\n","description":"","filename":"FigureS2.pdf","url":"https://assets-eu.researchsquare.com/files/rs-934219/v1/7f1b959a4d367ed621f0e0dd.pdf"},{"id":14083338,"identity":"6ab07a0b-dcec-46aa-a36e-4e69bcbde21b","added_by":"auto","created_at":"2021-09-28 19:27:04","extension":"pdf","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":1258565,"visible":true,"origin":"","legend":"Figure S3. Alternative splicing in Sf3b1 mutant GH3 cells.\n(A-N) PCR products of the 13 genes (not including Dlg1) – Pkd1l2, Myo15a, Tgif1, Rest, Mroh2a, Pipox, Emc1, Mapk15, Ccdc136, Tfb2m, Parp16, Parp4, and Cd44 were amplified from Sf3b1-mutant and WT GH3 cells. The cryptic and canonical transcripts were indicated by the arrow. Gapdh was used as a loading control. \n","description":"","filename":"FigureS3.pdf","url":"https://assets-eu.researchsquare.com/files/rs-934219/v1/6f248238f7db3b0af2ed4bd0.pdf"},{"id":14083113,"identity":"60435b36-542d-4cb9-84d4-bf22644befc0","added_by":"auto","created_at":"2021-09-28 19:21:04","extension":"pdf","order_by":4,"title":"","display":"","copyAsset":false,"role":"supplement","size":1676186,"visible":true,"origin":"","legend":"Figure S4. AS events in SF3B1 mutant human prolactinoma.\n(A) Volcano plot showing the differential alternative splicing events in SF3B1 mutant vs. wild-type human prolactinoma. Top 20 significant mis-spliced genes are indicated.\n\n(B) GO analysis of mis-spliced genes in SF3B1 mutant vs. wild-type human prolactinoma; the top five ranked terms are shown.\n","description":"","filename":"FigureS4.pdf","url":"https://assets-eu.researchsquare.com/files/rs-934219/v1/a7619b7f104c3b3657066033.pdf"},{"id":14083111,"identity":"22a0a7fd-c514-437b-9639-3b09892267ca","added_by":"auto","created_at":"2021-09-28 19:21:04","extension":"pdf","order_by":5,"title":"","display":"","copyAsset":false,"role":"supplement","size":1513473,"visible":true,"origin":"","legend":"Figure S5. Mutant SF3B1 caused aberrant splicing of DLG1 in 293T and GH3 cells.\n(A) Schematic representation of the validated aberrant splicing of DLG1 gene fragments by PCR. Schematic a represents the canonical DLG1 transcript and schematic b represents the aberrant DLG1 transcript. (B) The number of abnormal variant clones in the Ad-SF3B1-WT (2/48) and Ad-SF3B1-MT (7/43) groups, out of a total of 50 TA clones. (C) Sashimi plots of aberrant Dlg1 splicing event in GH3 cells with or without Sf3b1 mutations.\n","description":"","filename":"FigureS5.pdf","url":"https://assets-eu.researchsquare.com/files/rs-934219/v1/2fb89c6446d24e300f5f2e46.pdf"},{"id":14083115,"identity":"2d35810d-0057-499b-8ef5-c90b7db2f91f","added_by":"auto","created_at":"2021-09-28 19:21:04","extension":"xlsx","order_by":6,"title":"","display":"","copyAsset":false,"role":"supplement","size":788818,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementalTables113.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-934219/v1/905ad12d11280e7c0f26b5fe.xlsx"}],"financialInterests":"","formattedTitle":"\u003cp\u003eThe \u003cem\u003eSF3B1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e Mutation Promotes Prolactinoma Tumor Progression Through Aberrant Splicing Of\u003cem\u003e DLG1\u003c/em\u003e\u003c/p\u003e","fulltext":[{"header":"Background","content":"\u003cp\u003ePituitary adenomas account for between 10% and 20% of intracranial neoplasms, with prolactinomas being the most common subtype[\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e, \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e]. Except for hyperprolactinemia, a subset of prolactinomas are characterized by tumor invasion, resistance to conventional therapy, and high recurrence[\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e]. Therefore, deciphering the underlying pathogenesis of prolactinoma and determining efficient treatment targets are of great significance to the clinical management of prolactinoma patients.\u003c/p\u003e \u003cp\u003eAlternative splicing is an essential step in the posttranscriptional regulation of gene expression that is a highly regulated and complex process that diversifies the proteome by creating multiple proteins from the same gene[\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e, \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e]. Mutations in splicing factor cause aberrant alternative splicing, which is a key molecular characteristic of tumorigenesis[\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. Splicing factor 3b subunit 1 (SF3B1) is a core component of the U2 small nuclear ribonucleoprotein complex (U2 snRNP), which is essential for pre-mRNA splicing[\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e]. Recurrent somatic mutations in \u003cem\u003eSF3B1\u003c/em\u003e have been reported in myelodysplastic syndrome (MDS), chronic lymphocytic leukemia (CLL), and some solid tumors such as uveal melanoma, breast carcinoma, pancreas adenocarcinoma, and others[\u003cspan additionalcitationids=\"CR9 CR10 CR11 CR12\" citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]. Mutant SF3B1 leads to splicing defects that result in aberrantly spliced transcripts[\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e]. Approximately half of the aberrantly spliced mRNAs are subject to nonsense-mediated decay (NMD), which causes gene and protein production to be downregulated[\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e]. Our previous study found a recurrent \u003cem\u003eSF3B1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e mutation in prolactinoma, which was associated with poor progression-free survival (PFS) and higher levels of prolactin (PRL)[\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]. However, the underlying molecular mechanisms of this SF3B1 mutation and its downstream cellular processes in prolactinoma remain unclear.\u003c/p\u003e \u003cp\u003eHerein, we generated a heterozygous \u003cem\u003eSf3b1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e mutant rat cell line using CRISPR/Cas9 genome editing technology. RNA sequencing analysis revealed that while \u003cem\u003eSF3B1\u003c/em\u003e was conserved between rat and human, the aberrantly spliced mRNAs and affected genes differed significantly. Our results suggest that \u003cem\u003eSF3B1\u003c/em\u003e mutations in prolactinoma stimulate PI3K/AKT signaling by downregulating Discs large 1 (DLG1), which is induced by aberrant splicing and enhances the invasion and migration of tumor cells.\u003c/p\u003e"},{"header":"Methods","content":"\u003ch2\u003eCell culture and generation of \u003cem\u003eSf3b1\u003c/em\u003e mutant R625H cells\u003c/h2\u003e \u003cp\u003eThe GH3 rat pituitary cell line was purchased from the American Type Culture Collection (CLL-82.1; Manassas, VA, USA) and was cultured in Ham\u0026rsquo;s F12K medium with 2.5% fetal bovine serum (FBS) and 15% horse bovine serum (Gibco, Waltham, MA, USA). MCF7 and HEK293T cells were obtained from the National Infrastructure of Cell Line Resource (Beijing, China) and were cultured in DMEM (Gibco) supplemented with 10% FBS. PRL levels in cell culture supernatant were detected using an ELISA kit (BioVision, Milpitas, CA, USA) according to the manufacturer\u0026rsquo;s instructions.\u003c/p\u003e \u003cp\u003eThe CRISPR/Cas9 gene editing system was used to generate the \u003cem\u003eSf3b1\u003c/em\u003e p.R625H (c.G1874A, c.A1875T) mutant GH3 cell line. A Sanger Centre CRISPR webtool (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.sanger.ac.uk/htgt/wge/\u003c/span\u003e\u003c/span\u003e) was used to identify two small guide RNAs (sgRNAs). Cas9/sgRNA-mediated DNA double-strand break and homologous recombination contributed to generated the specific mutation. The 5\u0026prime; sgRNA sequence was 5\u0026prime;-GGCAAGATTCCTTCCTCA-3\u0026prime; and the 3\u0026prime; sgRNA sequence was 5\u0026prime;-GGTTAAGAGTACTGTTGTC-3\u0026prime;. The sgRNA pairs were cloned into a wild-type spCas9 and sgRNA expression plasmid. The donor plasmid (CL-GJ-016-puro-ΔTK) consisted of a puromycin resistance cassette flanked by two \u003cem\u003eloxP\u003c/em\u003e sites and ~\u0026thinsp;1 kb homology arms at both ends of the c. 1874 G\u0026thinsp;\u0026gt;\u0026thinsp;A, c. 1875 A\u0026thinsp;\u0026gt;\u0026thinsp;T mutant exon 14 (Wuhan Genecreate Biological Engineering Co. Ltd., Wuhan, China) (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). The plasmid was constructed according to previously published protocols and confirmed by sequencing. The Neon Transfection System (Thermo Fisher Scientific, Waltham, MA, USA) with pulse voltage of 1000 V for 40 ms was used to co-transfect GH3 cells on a 10-cm plate with 1 \u0026micro;g of CL-GJ-016-puro-TK and 3 \u0026micro;g of Cas9/sgRNA. The cells were screened with puromycin (0.5 \u0026micro;g/mL) starting 48-h post-transfection. Thereafter, single-cell cloning was performed in a 96-well plate. Sanger sequencing was performed using genomic DNA to confirm the specific mutations (DIA-UP Biotech, Beijing, China). Primers were designed to amplify exon 14 of the \u003cem\u003eSf3b1\u003c/em\u003e gene (Table S1).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eAdenoviral constructs and primary cell culture of human prolactinomas\u003c/h2\u003e \u003cp\u003eThe adenoviral constructs for mutant \u003cem\u003eSF3B1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e and wild-type \u003cem\u003eSFB31\u003c/em\u003e were generated by BAC Biological Technology (Beijing, China), and primary cultures of human prolactinoma cells were prepared as described in our previous publication[\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]. Tumor cells were infected with adenovirus at a multiplicity of infection of 100, and then harvested for quantitative (q)PCR 48-h later.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003eRNA sequencing\u003c/h2\u003e \u003cp\u003eFor GH3 cells, total RNA was extracted using the AllPrep\u0026reg; DNA/RNA Mini kit (QIAGEN, Hilden, Germany) according to the manufacturer\u0026rsquo;s instructions. Sequencing libraries were generated using NEBNext\u0026reg; Ultra\u0026trade; Directional RNA Library Prep Kit for Illumina\u0026reg; (NEB, Ipswich, MA, USA). The libraries were sequenced on an Illumina platform, and 150-bp paired-end reads were generated. Clean reads were mapped to the Rat Rnor_6.0 genome using Hisat2 (v2.0.5) to obtain read counts, fragments per kilobase per million mapped fragments (FPKM), and transcripts per million (TPM) for each sequenced gene.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003eDifferential alternative splicing and differential gene expression analyses\u003c/h2\u003e \u003cp\u003eDifferentially expressed genes (DEGs) were identified using DESeq2 using P\u0026thinsp;\u0026lt;\u0026thinsp;0.05 and absolute value of fold change\u0026thinsp;\u0026ge;\u0026thinsp;1.5. For alterative splicing analysis, rMATS (version 4.0.2)[\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e] was used to identify alternative splicing events by quantifying exon-exon junction spanning reads on annotated splice junctions present in the rat GENCODE Rnor_6.0 assembly. Differentially spliced mRNAs were defined as FDR\u0026thinsp;\u0026lt;\u0026thinsp;0.05 and a minimum inclusion level difference\u0026thinsp;\u0026gt;\u0026thinsp;10% or \u0026lt;\u0026thinsp;\u0026minus;\u0026thinsp;10%. Three mutant GH3 replicates and three wild-type replicates were compared. Data of two prolactinoma patients with the \u003cem\u003eSF3B1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e mutation and two wild-type cases were selected for alternative splicing analysis, and two mutant and 13 wild-type cases from our previously reports were selected for differential expression analysis[\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]. Basic clinic data of these patients are listed in Table S2. Functional enrichment analysis was performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) to predict the biological functions of DEGs and differentially spliced transcripts. The GO terms analyzed included biological process, cellular components, and molecular functions. A KEGG pathway with a P value\u0026thinsp;\u0026lt;\u0026thinsp;0.05 was considered to be significant.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003eWhole-genome sequencing and variant identification\u003c/h2\u003e \u003cp\u003eGenomic DNA was extracted from wild-type and mutant GH3 cells using the Blood \u0026amp; Cell Culture DNA Mini Kit (QIAGEN). DNA was sequenced using an Illumina HiSeq 2000, generating 150-bp paired reads. Reads were aligned to the Rat Rnor_6.0 genome using Burrows\u0026ndash;Wheeler Aligner (BWA)[\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e]. SAMtools was used to generate BAM files[\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e]. Visual inspection of the \u003cem\u003eSf3b1\u003c/em\u003e R625H mutation was performed using the Integrative Genomics Viewer[\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e].\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003eMinigene assay\u003c/h2\u003e \u003cp\u003eA DNA fragment containing the \u003cem\u003eDLG1\u003c/em\u003e exon 22 genomic sequence with 151 bp of flanking intron 21 and 271 bp of flanking intron 22 was inserted between the \u003cem\u003eKpnI\u003c/em\u003e and \u003cem\u003eEcoRI\u003c/em\u003e restriction sites of the pcMINI vector to produce the DLG1 minigene construct. Sanger sequencing was performed to confirm the sequence of the inserted fragment. Briefly, the \u003cem\u003eDLG1\u003c/em\u003e minigene and an adenoviral vector were co-transfected into 293T cells. After 48 h, RNA was harvested for PCR analysis. PCR products were separated by 2% agarose gels and confirmed by Sanger sequencing. The primers used for the spliced products were: 5\u0026prime;-CTAGAGAACCCACTGCTTAC-3\u0026prime; (forward) and 5\u0026prime;-TAGAAGGCACAGTCGAGG-3\u0026prime; (reverse).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003eTransfection\u003c/h2\u003e \u003cp\u003eRIBOBIO (Guangdong, China) synthesized small interfering (si)RNA duplexes; siRNA sequences of rat and human \u003cem\u003eDLG1\u003c/em\u003e are listed in Table S2. The \u003cem\u003eDlg1\u003c/em\u003e overexpressed plasmid, pLV-hef1a-Puro-WPRE-CMV-Dlg1-3\u0026times;FLAG, was constructed by Beijing Syngentech Co., Ltd. (Beijing, China). All transfections (siRNAs and overexpression plasmids) were performed using Lipofectamine\u0026reg; 3000 Transfection Reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer\u0026rsquo;s protocols. Cells were transfected with siRNA or plasmids for 48\u0026ndash;72 h. and then collected for qPCR and western blot assays.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec9\" class=\"Section2\"\u003e \u003ch2\u003eReverse transcription PCR and qPCR\u003c/h2\u003e \u003cp\u003eThe RNeasy Mini Kit (QIAGEN) was used to extract total RNA and the High Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific) was used to generate cDNA, both according to the manufacturers\u0026rsquo; instructions. For reverse transcription PCR, amplification was performed with I-5 High-Fidelity Master Mix (MCLAB, San Francisco, CA, USA) and specific primers. PCR products were separated on 1\u0026ndash;3% agarose gels. All qPCR assays were performed using Power SYBR\u0026trade; Green PCR Master Mix (Thermo Fisher Scientific) and were analyzed using QuantStudio 3 and 5 systems (Applied Biosystems, Waltham, MA, USA). The comparative Ct method was used to evaluate relative gene expression. The primers used in this study are listed in Table S3.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003eScanning electron microscopy (SEM)\u003c/h2\u003e \u003cp\u003eSEM was used to examine the external morphological alterations of \u003cem\u003eSf3b1\u003c/em\u003e wild-type and mutant GH3 cells. Cells were collected and fixed with 2.5% glutaraldehyde (Solarbio, Beijing, China) at 4\u0026deg;C overnight for SEM preparation. After rinsing with PBS and sterile water, the cells were dehydrated using an ethanol gradient. The samples were coated with gold after critical point drying, and electron micrographs were taken on a Hitachi SU8020 SEM (Tokyo, Japan).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003eCell migration and wound-healing assays\u003c/h2\u003e \u003cp\u003eCell migration culture dish inserts from Ibidi (Martinsried, Germany) were used to conduct wound-healing assays. After 24 h of transfection, cells were seeded into the chambers of the culture dish inserts. The inserts were removed on the next day, and fresh culture medium was added to each well. Scratches were photographed at different points in time using a Zeiss microscope (Oberkochen, Germany).\u003c/p\u003e \u003cp\u003e \u003cem\u003eSf3b1\u003c/em\u003e mutant and wild-type GH3 cells were seeded onto Imagelock 96-well plates (Essen Bioscience, Ann Arbor, MI, USA). The IncuCyte\u0026reg; Wound Maker (Essen Bioscience) was used to make uniform wounds in a monolayer of confluent cells. Phase contrast imaging was performed every 12 h for a total of 96 h. Images were analyzed using IncuCyte\u0026reg; S3 2018B-2019A software (Essen Bioscience), and data were analyzed using GraphPad Prism7 (GraphPad Software, Inc., La Jolla, CA, USA).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003ch2\u003eTranswell assays\u003c/h2\u003e \u003cp\u003eTranswell plates and Matrigel-coated transwell plates (Corning-Costar, Corning, NY, USA) were used to determine cell migration and invasion capabilities, respectively. Briefly, HEK293T cells resuspended in serum-free medium were inoculate into the upper chamber, and DMEM with 10% FBS was placed in the bottom chamber. Following 24 h of culture, cells on the bottom surface of the chamber were stained with crystal violet and counted under a microscope (Zeiss).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003eWestern blot analysis\u003c/h2\u003e \u003cp\u003eProtein samples were separated on 8\u0026ndash;10% Bis-Tris SDS-PAGE gels and transferred to polyvinylidene fluoride membranes (Merk, Kenilworth, NJ, USA). All primary antibodies (Table S4) were diluted in TBST containing 1% bovine serum albumin (BSA) and were incubated with the membranes overnight at 4\u0026deg;C. Immunoreactive bands were visualized by chemiluminescence.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec14\" class=\"Section2\"\u003e \u003ch2\u003eImmunohistochemistry\u003c/h2\u003e \u003cp\u003eHuman prolactinoma tumor specimens were used to examine DLG1, E-cadherin, and Snail protein levels in \u003cem\u003eSF3B1\u003c/em\u003e wild-type and mutant tumor tissues. The primary antibodies used are summarized in Table S4. Immunohistochemistry was performed by the Leica Bond Polymer Refine Detection system (Leica Biosystems, Wetzlar, Germany). All slides were scanned into digital images, and expression was examined using Aperio AT2 (Leica Biosystems). Staining intensity was scored as 0 (negative), 1 (weak), 2 (moderate), or 3 (strong). The percentage of immunostaining was recorded, and H-scores were calculated using the formula: H-score\u0026thinsp;=\u0026thinsp;1 \u0026times; (% weakly stained cells)\u0026thinsp;+\u0026thinsp;2 \u0026times; (% moderately stained cells)\u0026thinsp;+\u0026thinsp;3 \u0026times; (% strongly stained cells)], ranging from 0 to 300.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec15\" class=\"Section2\"\u003e \u003ch2\u003eRat prolactinoma model\u003c/h2\u003e \u003cp\u003eRat pituitary tumors were induced by subcutaneously implanting 1-cm silastic capsules containing 10 mg of 17-β estradiol in 4-week-old female F344 rats. Prolactinomas were induced by 17β-estradiol for 5 weeks, as described in our previous report[\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]. All experimental protocols were approved by the Animal Use and Care Committee of Beijing Tiantan Hospital. Five weeks later, the prolactinomas were validated via 7.0-T magnetic resonance imaging (MRI) before intra-pituitary injection. The rats were anaesthetized and then 1 \u0026micro;l of adenovirus vector control, wild-type \u003cem\u003eSFB31\u003c/em\u003e, or \u003cem\u003eSF3B1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e was stereotactically injected into each bilateral tumor. A second MRI examination was performed 2 weeks later, and all rats were sacrificed to collect tumor tissues for further analysis.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec16\" class=\"Section2\"\u003e \u003ch2\u003eImmunofluorescence staining of rat tumor tissues\u003c/h2\u003e \u003cp\u003eFirst, rats were anaesthetized and heart-perfused with 4% formalin. Tumor tissues were then collected and immersed in 10% sucrose for 2 h, and then incubated in 30% sucrose overnight before freezing in OCT compound. Prior to immunofluorescence staining, frozen 5-\u0026micro;m-thick sections were fixed in ice-cold acetone. Tissue slides were blocked with goat serum and incubated with primary antibodies overnight at 4\u0026deg;C (summarized in Table S4). After washing with PBS, the slides were incubated in Alexa Fluor 488 and 594 secondary antibodies (Invitrogen) for 1 h at room temperature; DAPI was used to visualize nuclei.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec17\" class=\"Section2\"\u003e \u003ch2\u003ePhalloidin staining and confocal microscopy\u003c/h2\u003e \u003cp\u003eCells were plated on confocal dishes coated with poly-L-lysine. After incubation at 37\u0026deg;C for 24 hours, the cells were fixed with 4% paraformaldehyde for 30 min, and then stained with 5 \u0026micro;g/mL Alexa Fluor 488-phalloidin (Invitrogen) to reveal filamentous actin (F-actin) in PBS for 40 min at 37\u0026deg;C. DAPI was used to visualize nuclei, and images were captured by confocal laser-scanning microscopy (Zeiss).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec18\" class=\"Section2\"\u003e \u003ch2\u003eStatistical analysis\u003c/h2\u003e \u003cp\u003eAll statistical analyses were performed using R v3.4.1 (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.r-project.org/\u003c/span\u003e\u003c/span\u003e) and Prism 7 (GraphPad Software, Inc.). All experiments were performed with at least three biological replicates, and all quantitative data represent mean\u0026thinsp;\u0026plusmn;\u0026thinsp;standard deviation (SD). Statistical significance was determined by unpaired Student\u0026rsquo;s \u003cem\u003et\u003c/em\u003e test (two groups) or one-way ANOVA (multiple groups). P\u0026thinsp;\u0026lt;\u0026thinsp;0.05 was considered statistically significant.\u003c/p\u003e \u003c/div\u003e"},{"header":"Results","content":"\u003ch2\u003eGeneration of \u003cem\u003eSf3b1\u003c/em\u003e R625H-mutant cells\u003c/h2\u003e\n\u003cp\u003eOur previous study identified the somatic hotspot mutation \u003cem\u003eSF3B1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e in 19.8% of prolactinomas [\u003cspan class=\"CitationRef\"\u003e16\u003c/span\u003e]. To investigate the biological role of the \u003cem\u003eSF3B1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e mutation on prolactinoma development and progression as well as similarities and differences between human and rats, we introduced the \u003cem\u003eSf3b1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e missense mutation into the rat pituitary cell line GH3 using the CRISPR/Cas9 gene editing system (Figure S1A). After isolating and expanding single cell clones, a heterozygous \u003cem\u003eSf3b1\u003c/em\u003e-R625H mutation in GH3 cells was confirmed by Sanger sequencing (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eA). Furthermore, whole-exome sequencing analyses of the mutant cell line also validated the heterozygous \u003cem\u003eSf3b1\u003c/em\u003e mutation, p.R625H (NC_005108.4:g.61608130T\u0026thinsp;\u0026gt;\u0026thinsp;A, g.61608131C\u0026thinsp;\u0026gt;\u0026thinsp;T, NM_053426.1:c.1874G\u0026thinsp;\u0026gt;\u0026thinsp;A, c.1875A\u0026thinsp;\u0026gt;\u0026thinsp;T), and the frequency of the mutant allele obtained from read counts of whole-exome sequencing was 28.57% (18 A\u0026amp;T and 45 G\u0026amp;A, Figure S1B). Additionally, to determine whether heterozygous \u003cem\u003eSf3b1\u003c/em\u003e-R625H mutation status affected PRL secretions, we performed ELISA to test PRL levels in culture media from wild-type and mutant cells. Consist with our previously reported results of mutant human prolactinoma[\u003cspan class=\"CitationRef\"\u003e16\u003c/span\u003e], PRL levels in supernatant from \u003cem\u003eSf3b1\u003c/em\u003e mutant cells was higher compared with wild-type cells (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eB). Together, these data demonstrate successful construction of a rat cell line model harboring the \u003cem\u003eSf3b1\u003c/em\u003e-R625H mutation to model human \u003cem\u003eSF3B1\u003c/em\u003e-mutant prolactinoma.\u003c/p\u003e\n\u003ch2\u003eHeterozygous \u003cem\u003eSf3b1\u003c/em\u003e-R625H status altered the transcriptome of GH3 cells\u003c/h2\u003e\n\u003cp\u003eTo explore transcriptome-level alterations in \u003cem\u003eSf3b1\u003c/em\u003e-R625H mutant cells, we conducted RNA-seq analysis from wild-type and mutant cultured cells. Using P\u0026thinsp;\u0026lt;\u0026thinsp;0.05 and |fold change|\u0026gt;1.5, we identified 3058 DEGs, with 1271 upregulated and 1787 downregulated genes (Figure S2A; Table S5). Clustering analysis showed distinctly different gene expression patterns between wild-type and \u003cem\u003eSf3b1\u003c/em\u003e-mutant cells (Figure S2B). Moreover, compared with \u003cem\u003eSF3B1\u003c/em\u003e-wild-type patient samples (\u003cem\u003en\u003c/em\u003e\u0026thinsp;=\u0026thinsp;13), we identified 1062 downregulated and 599 upregulated genes in \u003cem\u003eSF3B1\u003c/em\u003e-mutant human prolactinomas (\u003cem\u003en\u003c/em\u003e\u0026thinsp;=\u0026thinsp;2) with P\u0026thinsp;\u0026lt;\u0026thinsp;0.05 and |fold change|\u0026gt;1.5 (Figure S2C, 2D; Table S6). Comparing the two datasets, we found that only approximately 9% (285/3055) of the DEGs overlapped between rat GH3 cells and human prolactinoma samples (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eC).\u003c/p\u003e\n\u003cp\u003ePrevious studies have described that the alternative splicing defects caused by cancer-associated \u003cem\u003eSF3B1\u003c/em\u003e mutations involve incorrect recognition of the 3\u0026prime; splice site, producing aberrant transcripts [\u003cspan class=\"CitationRef\"\u003e15\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003e22\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003e23\u003c/span\u003e]. To determine whether \u003cem\u003eSf3b1\u003c/em\u003e-R625H is associated with mutation-specific alterations in pre-mRNA splicing, we performed percent spliced-in (PSI) analyses of alternative splicing using the rMATS tool[\u003cspan class=\"CitationRef\"\u003e17\u003c/span\u003e]. At a False Discovery Rate (FDR) of \u0026lt;\u0026thinsp;0.05 and |\u0026Delta;PSI|\u0026gt;0.1, we obtained 143 alternative splicing events in \u003cem\u003eSf3b1\u003c/em\u003e-R625H cells (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eD, Table S7), with exon skipping being the most frequent event (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eE). GO functional analysis of the alternatively spliced genes revealed enrichment of ontologies including lamellipodium organization, regulation of cell morphogenesis, and motor activity (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eF, Table S8). We also performed qPCR to confirm several alternative splicing-associated changes (Figure S3A). Additionally, re-analysis of RNA-seq data generated from two \u003cem\u003eSF3B1\u003c/em\u003e-mutant and two \u003cem\u003eSF3B1\u003c/em\u003e-wild-type prolactinoma cases[\u003cspan class=\"CitationRef\"\u003e16\u003c/span\u003e] uncovered 2287 alternative splicing events in \u003cem\u003eSF3B1\u003c/em\u003e-mutant cases, with a higher proportion of inclusion alternative 3\u0026prime; splice site (\u0026Delta;PSI\u0026thinsp;\u0026gt;\u0026thinsp;10%) and inclusion retained intron (\u0026Delta;PSI\u0026thinsp;\u0026gt;\u0026thinsp;10%) events (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eG; Figure S3B-C; Table S9-10). Similar to human data, \u003cem\u003eSf3b1\u003c/em\u003e mutant GH3 cells had more inclusion (\u003cem\u003en\u003c/em\u003e\u0026thinsp;=\u0026thinsp;83, \u0026Delta;PSI\u0026thinsp;\u0026gt;\u0026thinsp;10%) than exclusion (\u003cem\u003en\u003c/em\u003e\u0026thinsp;=\u0026thinsp;60, \u0026Delta;PSI\u0026thinsp;\u0026lt;\u0026thinsp;\u0026minus;\u0026thinsp;10%) alternative splicing events compared with wild-type cells (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eE and G). Even so, there was minimal overlap between aberrantly spliced genes, and only approximately 15% (18/120) of the aberrantly spliced genes in rat \u003cem\u003eSf3b1\u003c/em\u003e-mutant GH3 cells were common to human samples (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eH). These data indicate that the limited conservation between rat and human intronic sequences caused the mutant SF3B1 to have different effects on alternative mRNA splicing in the two species.\u003c/p\u003e\n\u003cdiv class=\"Section2\" id=\"Sec20\"\u003e\n \u003ch2\u003eThe R625H mutation induced epithelial-mesenchymal transition (EMT) phenotypes in tumor cells\u003c/h2\u003e\n \u003cp\u003eOur previous study found that most \u003cem\u003eSF3B1\u003c/em\u003e-mutant tumors exhibit invasive behavior, which is associated with poor PFS in prolactinoma[\u003cspan class=\"CitationRef\"\u003e24\u003c/span\u003e]. Furthermore, GO enrichment analysis of the DEGs and alternative splicing events derived from GH3 cells revealed that cell adhesion molecules and lamellipodium organization were significantly altered in \u003cem\u003eSf3b1\u003c/em\u003e-mutant GH3 cells (Figs.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eA and \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eF, Tables S8 and S11). These results indicated that \u003cem\u003eSFB31\u003c/em\u003e-mutant tumor cells have invasive characteristics. To examine changes in the adhesive properties of mutant cells, SEM was performed. Compared with \u003cem\u003eSf3b1\u003c/em\u003e wild-type cells, which showed a relatively rounded morphology, \u003cem\u003eSf3b1\u003c/em\u003e mutant cells showed prominent lamellipodial extensions (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eB). Furthermore, F-actin polymerization and filopodium formation were observed in \u003cem\u003eSf3b1\u003c/em\u003e mutant cells (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eB). Lamellipodia formation and dramatic reorganization of the actin cytoskeleton are involved in EMT, which is associated with tumor cell invasion and metastasis[\u003cspan class=\"CitationRef\"\u003e25\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003e26\u003c/span\u003e]. To explore whether \u003cem\u003eSf3b1\u003c/em\u003e mutation promotes cell migration, wound-healing assays was conducted. The migration ability of \u003cem\u003eSf3b1\u003c/em\u003e mutant cells was significantly enhanced compared with wild-type cells (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eC). Moreover, the influence of \u003cem\u003eSf3b1\u003c/em\u003e mutation on cell migration was also verified by the IncuCyte ZOOM\u0026reg; 96-well Scratch Wound cell migration assay (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eD).\u003c/p\u003e\n \u003cp\u003eTo determine whether the \u003cem\u003eSf3b1\u003c/em\u003e mutation modulated cell migration and invasion through the EMT pathway, we analyzed EMT maker expression. We observed a significant decrease in mesenchymal maker expression and an increase in epithelial maker expression in mutant cells. Immunoblotting of \u003cem\u003eSf3b1\u003c/em\u003e mutant cells revealed downregulated expression of E-cadherin along with elevated expression of N-cadherin, Vimentin, and Snail (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eE). Consistently, immunohistochemistry of consecutive sections from wild-type and mutant prolactinoma tissues supported these findings (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eF). E-cadherin expression levels yielded an H-score of 146.7\u0026thinsp;\u0026plusmn;\u0026thinsp;17.64 in wild-type samples and 53.3\u0026thinsp;\u0026plusmn;\u0026thinsp;14.53 in mutant samples (P\u0026thinsp;=\u0026thinsp;0.015). The H-scores of Snail were 20\u0026thinsp;\u0026plusmn;\u0026thinsp;5.78 and 196.7\u0026thinsp;\u0026plusmn;\u0026thinsp;20.28, respectively (P\u0026thinsp;=\u0026thinsp;0.0011). To explore if \u003cem\u003eSF3B1\u003c/em\u003e-R625H promoted EMT of prolactinoma \u003cem\u003ein vivo\u003c/em\u003e, F344 rat prolactinomas were induced by 6-weel 17\u0026beta;-estradiol treatment. Then the rat prolactinomas were stereotactically injected with wild-type and R625H \u003cem\u003eSF3B1\u003c/em\u003e encoding adenoviruses and negative controls. After 2 weeks, tumors were collected and immunofluorescence was used to examine EMT maker expression. Decreased levels of E-cadherin and increased levels of Vimentin were observed in the R625H groups (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eG). Overall, these data demonstrated that mutant SF3B1 promoted EMT.\u003c/p\u003e\n \u003cp\u003eTo explore the molecular mechanisms through which the \u003cem\u003eSF3B1\u003c/em\u003e mutation induced EMT, we analyzed the intersection of aberrantly spliced genes and DEGs in rat GH3 cells and human prolactinoma samples (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eH). From this intersection, \u003cem\u003eDLG1\u003c/em\u003e was the only identified gene (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eH), and was significantly downregulated and mis-spliced in \u003cem\u003eSF3B1\u003c/em\u003e mutants (Tables S5\u0026ndash;S7 and S9).\u003c/p\u003e\n \u003ch2\u003eMutant SF3B1 caused aberrant splicing of \u003cem\u003eDLG1\u003c/em\u003e\u003c/h2\u003e\n \u003cp\u003eDLG1 is a member of the molecular scaffold protein family knows as membrane associated guanylate kinases (MAGUKs)[\u003cspan class=\"CitationRef\"\u003e27\u003c/span\u003e]. DLG1 is a component of the conserved Scribble polarity complex, which is associated with the establishment and maintenance of cell polarity and is defined as a potential tumor suppressor[\u003cspan class=\"CitationRef\"\u003e28\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003e29\u003c/span\u003e]. To confirm that aberrant splicing of \u003cem\u003eDLG1\u003c/em\u003e was caused by \u003cem\u003eSF3B1\u003c/em\u003e mutation, RNA-seq data from human prolactinoma and GH3 cells was examined. Across human prolactinoma, aberrant splicing of \u003cem\u003eDLG1\u003c/em\u003e was induced by mutant SF3B1 through usage of a cryptic 3\u0026prime; splice site (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eA; Table S9), which was confirmed by qPCR in prolactinoma patient samples (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eB). To verify the effect of mutant SF3B1 on \u003cem\u003eDLG1\u003c/em\u003e in vitro, we infected primary cultured pituitary tumor cells with adenovirus carrying the \u003cem\u003eSF3B1\u003c/em\u003e-R625H mutation, and the results showed that the cryptic \u003cem\u003eDLG1\u003c/em\u003e transcript was observed in the Ad-SF3B1-R625H group (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eC). The same results were also observed in HEK293T cells (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eD and S4A). To verify differences between the wild-type and mutant alleles, we performed TA cloning. Through reverse transcription PCR and sequencing clones, we confirmed that the Ad-SF3B1-R625H group included more abnormal PCR products than the Ad-SF3B1-WT group (Figure S4B). The expected fragment with a 16-bp extension of exon 22 was identified by Sanger sequencing of PCR products (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eE). We then generated a minigene construct that contained sequences from the misspliced intron and flanking exon to explore alternative splicing (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eF). After co-transfection into HEK293T cells, a complete spliced RNA containing exon 22 was observed in the \u003cem\u003eSF3B1\u003c/em\u003e-WT and control groups (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eG, left and middle lanes); however, a larger transcript was observed in the \u003cem\u003eSF3B1\u003c/em\u003e-R625H group that retained an extra 16 bp (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eG, right lane). The extra 16-bp fragment was confirmed by Sanger sequencing (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eH). Additionally, mutually exclusive splicing of \u003cem\u003eDlg1\u003c/em\u003e in \u003cem\u003eSf3b1\u003c/em\u003e-mutant GH3 cells was confirmed by reverse transcription PCR (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eI; Table S7). These data identified that mutant \u003cem\u003eSF3B1\u003c/em\u003e induced aberrant splicing of \u003cem\u003eDLG1\u003c/em\u003e.\u003c/p\u003e\n \u003cp\u003eTo verify whether the \u003cem\u003eSF3B1\u003c/em\u003e mutation altered \u003cem\u003eDLG1\u003c/em\u003e expression, RNA-seq data was examined and qPCR was performed in \u003cem\u003eSF3B1\u003c/em\u003e mutant samples. These data revealed that \u003cem\u003eDLG1\u003c/em\u003e was downregulated in \u003cem\u003eSF3B1\u003c/em\u003e-mutant human prolactinoma samples and GH3 cells (Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eA, B). Furthermore, immunohistochemistry and western blotting revealed that DLG1 protein expression was decreased in \u003cem\u003eSF3B1\u003c/em\u003e-mutant human samples and GH3 cells (Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eC and \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eE). The H-scores of DLG1 were 153\u0026thinsp;\u0026plusmn;\u0026thinsp;17.64 in human \u003cem\u003eSF3B1\u003c/em\u003e-wild-type and 40\u0026thinsp;\u0026plusmn;\u0026thinsp;20.82 in human \u003cem\u003eSF3B1\u003c/em\u003e-mutant samples (P\u0026thinsp;=\u0026thinsp;0.0142). To further confirm the effect of the \u003cem\u003eSF3B1\u003c/em\u003e mutation on DLG1 expression, immunofluorescence was performed on tumors stereotactically injected with adenoviruses. Decreased DLG1 expression was observed in the Ad-SF3B1 R625H group (Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eD). Taken together, these data demonstrate that mutant SF3B1 induced aberrant splicing and reduced expression of \u003cem\u003eDLG1\u003c/em\u003e.\u003c/p\u003e\n \u003ch2\u003eKnocking down \u003cem\u003eDLG1\u003c/em\u003e promoted cell migration, invasion, and EMT\u003c/h2\u003e\n \u003cp\u003eTo evaluate the role of DLG1 in tumor cell invasion and migration, we knocked down \u003cem\u003eDLG1\u003c/em\u003e by transfecting specific siRNAs. Both the mRNA and protein levels of DLG1 were efficiently depleted in GH3, 293T, and MCF7 cells (Fig. \u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003eA-B). Western blot results revealed that knocking down \u003cem\u003eDLG1\u003c/em\u003e promoted EMT by decreasing the expression of E-cadherin and increasing the expression of N-cadherin, Vimentin, and Snail (Fig. \u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003eB). Furthermore, restoring DLG1 expression in wild-type and mutant GH3 cells inhibited the EMT program (Fig. \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eE). Wound healing assays revealed that si-DLG1-transfected cells traversed wounds significantly faster than the si-NC-transfected cells (Fig. \u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003eC). Comparable migration differences were also observed in the transwell assay (Fig. \u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003eD). In transwell Matrigel invasion assays, significantly more si-DLG1-transfected cells migrated through the Matrigel layer (Fig. \u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003eD). Furthermore, phalloidin staining showed differences in F-actin organization at the wound margin. Compared with si-NC-transfected cells, si-DLG1-transfected cells showed polarized lamellipodia formation \u003cem\u003evia\u003c/em\u003e F-actin (phalloidin) staining (Fig. \u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003eE). These results suggested that depleting DLG1 enhanced tumor cell migration and invasion, and induced EMT.\u003c/p\u003e\n \u003ch2\u003eThe SFB31/DLG1 axis promoted tumor invasion \u003cem\u003evia\u003c/em\u003e the PI3K/AKT pathway\u003c/h2\u003e\n \u003cp\u003eTo determine the key biological mechanism through which mutant SF3B1 promoted prolactinoma progression, we performed a function enrichment analysis of the DEGs from human tumors and rat GH3 cells. This revealed that PI3K/AKT signaling was significantly altered in mutant SF3B1 samples compared with wild type samples (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e6\u003c/span\u003eA, Tables S12-13). We therefore detected the relevant effectors within this pathway through western blotting and found increased phosphorylation of AKT (T308) and GSK-3\u0026beta; in \u003cem\u003eSf3b1\u003c/em\u003e-mutant GH3 cells (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e6\u003c/span\u003eB). To clarify the specific role of DLG1 in activation of the PI3K/AKT pathway in tumor cells, we manipulated DLG1 activity in GH3, 293T, and MCF7 cells. Knocking down \u003cem\u003eDLG1\u003c/em\u003e in cells by siRNA increased phosphorylation of AKT and GSK-3\u0026beta; (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e6\u003c/span\u003eC). Furthermore, restoration of DLG1 expression in wild-type and mutant GH3 cells showed reduced phosphorylation of AKT and GSK-3\u0026beta; (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e6\u003c/span\u003eB). Thus, these data strongly suggest that altering DLG1 levels due to alternative splicing induced by mutant SF3B1 promotes PI3K/AKT signaling in prolactinoma.\u003c/p\u003e\n \u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003e \u003cem\u003eSF3B1\u003c/em\u003e is the most frequently mutated splicing gene in cancer, with three major hotspots: the codon positions K700, R625, and R622[\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e]. Among these mutated codons, K700E is the most common mutation in CLL, MDS, acute myeloid leukemia, and invasive breast carcinoma, while R625 is the predominant mutation in uveal melanoma and skin cutaneous melanoma[\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e]. Previous studies have explored how the \u003cem\u003eSF3B1\u003c/em\u003e\u003csup\u003e\u003cem\u003eK700E\u003c/em\u003e\u003c/sup\u003e mutation contributes to aberrant splicing and downstream pro-tumorigenesis mechanisms in CLL and breast cancer using cell lines and mouse disease models[\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e, \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e]. To understand the molecular and phenotypic consequences of the \u003cem\u003eSF3B1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e mutation on prolactinoma tumorigenesis, we generated a heterozygous \u003cem\u003eSf3b1\u003c/em\u003e-R625H mutant rat pituitary cell line. We then demonstrated that heterozygosity for \u003cem\u003eSf3b1\u003c/em\u003e-R625H in GH3 cells was sufficient to cause characteristic features similar to human prolactinoma, including higher PRL levels and invasive behaviors. Although phylogenetic analysis revealed that the \u003cem\u003eSF3B1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e locus is highly conserved across species[\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e], only a few events were shared between mutant human prolactinoma and rat cells. These results may be due to the frequency of the mutant allele in heterozygous mutant cells and relatively poor interspecies conservation of intronic sequences between rats and humans. Moreover, our study found that the majority of aberrant splicing events in \u003cem\u003eSF3B1\u003c/em\u003e mutant human prolactinoma were retained intron followed by alternative 3\u0026prime; splice site. This is different from previous observations in \u003cem\u003eSF3B1\u003c/em\u003e-mutant cancers[\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e, \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e], which suggested a higher proportion of alternative 3\u0026prime; splice site events, but were similar to a recent report in \u003cem\u003eSF3B1\u003c/em\u003e-mutant MDS[\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e]. Considering the limited number \u003cem\u003eSF3B1\u003c/em\u003e-mutant human prolactinoma tumor samples, further studies with larger cohorts are needed to confirm these results. In particular, we consistently found that the \u003cem\u003eDLG1\u003c/em\u003e transcript was shared between rat cells and human prolactinoma samples, with mRNA downregulation and differential missplicing. Interestingly, mutant \u003cem\u003eSF3B1\u003c/em\u003e induced the inclusion of 16 intronic nucleotides in exon 22 by using an upstream cryptic 3\u0026prime; splice site in human sequences (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eA), which may lead to NMD and reduced DLG1 expression, because the shifted reading frame produces a premature termination codon in this alternative spliceform[\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e]. Despite mutant SF3B1 inducing different splicing events in rats and humans, \u003cem\u003eDLG1\u003c/em\u003e mRNA and protein were downregulated in both models.\u003c/p\u003e \u003cp\u003eOur data demonstrated that downregulating \u003cem\u003eDLG1\u003c/em\u003e expression promote tumor cell migration and invasion. As previously reported, DLG1 was involved in mammalian tumorigenesis and was localized to adherens junctions in epithelial cells[\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e]. In mammals, \u003cem\u003eDLG1\u003c/em\u003e appears to act as a tumor suppressor by affecting other polarity complexes, Myosin II activity, the actin cytoskeleton, and/or other signaling pathways[\u003cspan additionalcitationids=\"CR36\" citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e]. Our studies further indicated that knocking down \u003cem\u003eDLG1\u003c/em\u003e induced EMT in tumor cells. EMT is well known to be related to the loss of cell polarity, specifically the loss of adherens junctions and enables the development of an invasive phenotype[\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e]. Consist with our results, the study by Sugihara, T. \u003cem\u003eet al.\u003c/em\u003e found that loss of DLG1 was associated with poor prognosis in endometrial cancer and that knocking down \u003cem\u003eDLG1\u003c/em\u003e accelerated tumor migration and invasion \u003cem\u003ein vitro\u003c/em\u003e[\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e].\u003c/p\u003e \u003cp\u003ePrevious studies have revealed that the AKT/GSK-3β/Snail axis is critical for the induction and maintenance of EMT [\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e, \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e]. In this study, we show that mutant SF3B1 and decreased DLG1 expression activate the PI3K/AKT pathway, which results in tumor cells undergoing EMT and gaining migration and invasive properties (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e). It is known that p-AKT can suppress GSK-3β activity (through phosphorylation of Ser 9) and stabilize Snail[\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e]. Furthermore, DLG1 has been demonstrated to modulate the PI3K/Akt pathway to precisely regulate regulatory T cell activity[\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e]. Consistent with this, we found that low DLG1 expression increased p-AKT, p-G3K-3β, and Snail expression levels, which could be rescued by overexpressing DLG1. Future work comprehensively defining the how DLG1 activates the PI3K/Akt signaling pathway in tumor cell may therefore be very important.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e"},{"header":"Conclusions","content":"\u003cp\u003eIn summary, we confirmed the molecular mechanism through which the \u003cem\u003eSF3B1\u003c/em\u003e-R625H mutation actives the PI3K/Akt pathway in prolactinoma by aberrant splicing of \u003cem\u003eDLG1\u003c/em\u003e, which promotes tumor invasion progression. Our findings provide the rationale for investigating new therapeutic strategies in patients with \u003cem\u003eSF3B1\u003c/em\u003e-mutant prolactinoma.\u003c/p\u003e "},{"header":"Abbreviations","content":"\u003cp\u003eMDS \u0026nbsp; \u0026nbsp; \u0026nbsp;Myelodysplastic syndrome\u003c/p\u003e\n\u003cp\u003eCLL \u0026nbsp; \u0026nbsp; \u0026nbsp; Chronic lymphocytic leukemia\u003c/p\u003e\n\u003cp\u003eNMD \u0026nbsp; \u0026nbsp; \u0026nbsp;Nonsense-mediated decay\u003c/p\u003e\n\u003cp\u003ePFS \u0026nbsp; \u0026nbsp; \u0026nbsp; Progression-free survival\u003c/p\u003e\n\u003cp\u003eFBS \u0026nbsp; \u0026nbsp; \u0026nbsp; Fetal bovine serum\u003c/p\u003e\n\u003cp\u003eFPKM \u0026nbsp; \u0026nbsp; Fragments per kilobase per million mapped fragments\u003c/p\u003e\n\u003cp\u003eTPM \u0026nbsp; \u0026nbsp; \u0026nbsp;Transcripts per million\u003c/p\u003e\n\u003cp\u003eDEGs \u0026nbsp; \u0026nbsp; Differentially expressed genes\u003c/p\u003e\n\u003cp\u003eGO \u0026nbsp; \u0026nbsp; \u0026nbsp; Gene Ontology\u003c/p\u003e\n\u003cp\u003eKEGG \u0026nbsp; \u0026nbsp; Kyoto Encyclopedia of Genes and Genomes\u003c/p\u003e\n\u003cp\u003eBSA \u0026nbsp; \u0026nbsp; \u0026nbsp; Bovine serum albumin\u003c/p\u003e\n\u003cp\u003eMRI \u0026nbsp; \u0026nbsp; \u0026nbsp; Magnetic resonance imaging\u003c/p\u003e\n\u003cp\u003eMAGUKs \u0026nbsp;Membrane associated guanylate kinases\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe present study was approved by The Ethics Committees of Beijing Tiantan Hospital. All subjects provided written informed consent. The animal\u0026nbsp;experiments were approved by the Animal Use and Care Committee of Beijing Tiantan Hospital.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no competing interests.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis study was supported by the National Natural Science Foundation of China (Grant codes: 81672495, 81771489, 82072804, 82071559, and 82071558).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors\u0026apos; contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWX and YZ worked on the conception and designed the research. CL and DW were involved in the collection and analysis of patients\u0026rsquo; clinical data. JG, QF and YL performed the experiments.\u0026nbsp;JG was dedicated to data analysis, interpretation, and drafting. All authors read and approved the final manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eAsa SL, Ezzat S: \u003cstrong\u003eThe cytogenesis and pathogenesis of pituitary adenomas\u003c/strong\u003e. \u003cem\u003eEndocr Rev \u003c/em\u003e1998, \u003cstrong\u003e19\u003c/strong\u003e(6):798-827.\u003c/li\u003e\n\u003cli\u003eDaly AF, Rixhon M, Adam C, Dempegioti A, Tichomirowa MA, Beckers A: \u003cstrong\u003eHigh prevalence of pituitary adenomas: a cross-sectional study in the province of Liege, Belgium\u003c/strong\u003e. \u003cem\u003eJ Clin Endocrinol Metab \u003c/em\u003e2006, \u003cstrong\u003e91\u003c/strong\u003e(12):4769-4775.\u003c/li\u003e\n\u003cli\u003eVroonen L, Daly AF, Beckers A: \u003cstrong\u003eEpidemiology and Management Challenges in Prolactinomas\u003c/strong\u003e. \u003cem\u003eNeuroendocrinology \u003c/em\u003e2019, \u003cstrong\u003e109\u003c/strong\u003e(1):20-27.\u003c/li\u003e\n\u003cli\u003eNilsen TW, Graveley BR: \u003cstrong\u003eExpansion of the eukaryotic proteome by alternative splicing\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e2010, \u003cstrong\u003e463\u003c/strong\u003e(7280):457-463.\u003c/li\u003e\n\u003cli\u003eLee Y, Rio DC: \u003cstrong\u003eMechanisms and Regulation of Alternative Pre-mRNA Splicing\u003c/strong\u003e. \u003cem\u003eAnnu Rev Biochem \u003c/em\u003e2015, \u003cstrong\u003e84\u003c/strong\u003e:291-323.\u003c/li\u003e\n\u003cli\u003eSveen A, Kilpinen S, Ruusulehto A, Lothe RA, Skotheim RI: \u003cstrong\u003eAberrant RNA splicing in cancer; expression changes and driver mutations of splicing factor genes\u003c/strong\u003e. \u003cem\u003eOncogene \u003c/em\u003e2016, \u003cstrong\u003e35\u003c/strong\u003e(19):2413-2427.\u003c/li\u003e\n\u003cli\u003eKramer A: \u003cstrong\u003eThe structure and function of proteins involved in mammalian pre-mRNA splicing\u003c/strong\u003e. \u003cem\u003eAnnu Rev Biochem \u003c/em\u003e1996, \u003cstrong\u003e65\u003c/strong\u003e:367-409.\u003c/li\u003e\n\u003cli\u003ePapaemmanuil E, Cazzola M, Boultwood J, Malcovati L, Vyas P, Bowen D, Pellagatti A, Wainscoat JS, Hellstrom-Lindberg E, Gambacorti-Passerini C\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eSomatic SF3B1 mutation in myelodysplasia with ring sideroblasts\u003c/strong\u003e. \u003cem\u003eN Engl J Med \u003c/em\u003e2011, \u003cstrong\u003e365\u003c/strong\u003e(15):1384-1395.\u003c/li\u003e\n\u003cli\u003eWang L, Lawrence MS, Wan Y, Stojanov P, Sougnez C, Stevenson K, Werner L, Sivachenko A, DeLuca DS, Zhang L\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eSF3B1 and other novel cancer genes in chronic lymphocytic leukemia\u003c/strong\u003e. \u003cem\u003eN Engl J Med \u003c/em\u003e2011, \u003cstrong\u003e365\u003c/strong\u003e(26):2497-2506.\u003c/li\u003e\n\u003cli\u003eBanerji S, Cibulskis K, Rangel-Escareno C, Brown KK, Carter SL, Frederick AM, Lawrence MS, Sivachenko AY, Sougnez C, Zou L\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eSequence analysis of mutations and translocations across breast cancer subtypes\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e2012, \u003cstrong\u003e486\u003c/strong\u003e(7403):405-409.\u003c/li\u003e\n\u003cli\u003eBiankin AV, Waddell N, Kassahn KS, Gingras MC, Muthuswamy LB, Johns AL, Miller DK, Wilson PJ, Patch AM, Wu J\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003ePancreatic cancer genomes reveal aberrations in axon guidance pathway genes\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e2012, \u003cstrong\u003e491\u003c/strong\u003e(7424):399-405.\u003c/li\u003e\n\u003cli\u003eFurney SJ, Pedersen M, Gentien D, Dumont AG, Rapinat A, Desjardins L, Turajlic S, Piperno-Neumann S, de la Grange P, Roman-Roman S\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eSF3B1 mutations are associated with alternative splicing in uveal melanoma\u003c/strong\u003e. \u003cem\u003eCancer Discov \u003c/em\u003e2013, \u003cstrong\u003e3\u003c/strong\u003e(10):1122-1129.\u003c/li\u003e\n\u003cli\u003eArmenia J, Wankowicz SAM, Liu D, Gao J, Kundra R, Reznik E, Chatila WK, Chakravarty D, Han GC, Coleman I\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eThe long tail of oncogenic drivers in prostate cancer\u003c/strong\u003e. \u003cem\u003eNat Genet \u003c/em\u003e2018, \u003cstrong\u003e50\u003c/strong\u003e(5):645-651.\u003c/li\u003e\n\u003cli\u003ePellagatti A, Armstrong RN, Steeples V, Sharma E, Repapi E, Singh S, Sanchi A, Radujkovic A, Horn P, Dolatshad H\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eImpact of spliceosome mutations on RNA splicing in myelodysplasia: dysregulated genes/pathways and clinical associations\u003c/strong\u003e. \u003cem\u003eBlood \u003c/em\u003e2018, \u003cstrong\u003e132\u003c/strong\u003e(12):1225-1240.\u003c/li\u003e\n\u003cli\u003eDarman RB, Seiler M, Agrawal AA, Lim KH, Peng S, Aird D, Bailey SL, Bhavsar EB, Chan B, Colla S\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eCancer-Associated SF3B1 Hotspot Mutations Induce Cryptic 3\u0026apos; Splice Site Selection through Use of a Different Branch Point\u003c/strong\u003e. \u003cem\u003eCell reports \u003c/em\u003e2015, \u003cstrong\u003e13\u003c/strong\u003e(5):1033-1045.\u003c/li\u003e\n\u003cli\u003eLi C, Xie W, Rosenblum JS, Zhou J, Guo J, Miao Y, Shen Y, Wang H, Gong L, Li M\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eSomatic SF3B1 hotspot mutation in prolactinomas\u003c/strong\u003e. \u003cem\u003eNat Commun \u003c/em\u003e2020, \u003cstrong\u003e11\u003c/strong\u003e(1):2506.\u003c/li\u003e\n\u003cli\u003eShen S, Park JW, Lu ZX, Lin L, Henry MD, Wu YN, Zhou Q, Xing Y: \u003cstrong\u003erMATS: robust and flexible detection of differential alternative splicing from replicate RNA-Seq data\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2014, \u003cstrong\u003e111\u003c/strong\u003e(51):E5593-5601.\u003c/li\u003e\n\u003cli\u003eLi H, Durbin R: \u003cstrong\u003eFast and accurate short read alignment with Burrows-Wheeler transform\u003c/strong\u003e. \u003cem\u003eBioinformatics \u003c/em\u003e2009, \u003cstrong\u003e25\u003c/strong\u003e(14):1754-1760.\u003c/li\u003e\n\u003cli\u003eLi H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Genome Project Data Processing S: \u003cstrong\u003eThe Sequence Alignment/Map format and SAMtools\u003c/strong\u003e. \u003cem\u003eBioinformatics \u003c/em\u003e2009, \u003cstrong\u003e25\u003c/strong\u003e(16):2078-2079.\u003c/li\u003e\n\u003cli\u003eRobinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP: \u003cstrong\u003eIntegrative genomics viewer\u003c/strong\u003e. \u003cem\u003eNat Biotechnol \u003c/em\u003e2011, \u003cstrong\u003e29\u003c/strong\u003e(1):24-26.\u003c/li\u003e\n\u003cli\u003eLei C, Jing G, Jichao W, Xiaohui L, Fang Q, Hua G, Yazhou M, Zhang Y: \u003cstrong\u003eMiR-137\u0026apos;s Tumor Suppression on Prolactinomas by Targeting MITF and Modulating Wnt Signaling Pathway\u003c/strong\u003e. \u003cem\u003eJ Clin Endocrinol Metab \u003c/em\u003e2019, \u003cstrong\u003e104\u003c/strong\u003e(12):6391-6402.\u003c/li\u003e\n\u003cli\u003eWang L, Brooks AN, Fan J, Wan Y, Gambe R, Li S, Hergert S, Yin S, Freeman SS, Levin JZ\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eTranscriptomic Characterization of SF3B1 Mutation Reveals Its Pleiotropic Effects in Chronic Lymphocytic Leukemia\u003c/strong\u003e. \u003cem\u003eCancer cell \u003c/em\u003e2016, \u003cstrong\u003e30\u003c/strong\u003e(5):750-763.\u003c/li\u003e\n\u003cli\u003eShiozawa Y, Malcovati L, Galli A, Sato-Otsubo A, Kataoka K, Sato Y, Watatani Y, Suzuki H, Yoshizato T, Yoshida K\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eAberrant splicing and defective mRNA production induced by somatic spliceosome mutations in myelodysplasia\u003c/strong\u003e. \u003cem\u003eNat Commun \u003c/em\u003e2018, \u003cstrong\u003e9\u003c/strong\u003e(1):3649.\u003c/li\u003e\n\u003cli\u003eKim K, Park YW, Kim D, Ahn SS, Moon JH, Kim EH, Lee EJ, Ku CR: \u003cstrong\u003eBiochemical Remission after Cabergoline Withdrawal in Hyperprolactinemic Patients with Visible Remnant Pituitary Adenoma\u003c/strong\u003e. \u003cem\u003eJ Clin Endocrinol Metab \u003c/em\u003e2021, \u003cstrong\u003e106\u003c/strong\u003e(2):e615-e624.\u003c/li\u003e\n\u003cli\u003eYilmaz M, Christofori G: \u003cstrong\u003eEMT, the cytoskeleton, and cancer cell invasion\u003c/strong\u003e. \u003cem\u003eCancer Metastasis Rev \u003c/em\u003e2009, \u003cstrong\u003e28\u003c/strong\u003e(1-2):15-33.\u003c/li\u003e\n\u003cli\u003eAiello NM, Maddipati R, Norgard RJ, Balli D, Li J, Yuan S, Yamazoe T, Black T, Sahmoud A, Furth EE\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eEMT Subtype Influences Epithelial Plasticity and Mode of Cell Migration\u003c/strong\u003e. \u003cem\u003eDev Cell \u003c/em\u003e2018, \u003cstrong\u003e45\u003c/strong\u003e(6):681-695 e684.\u003c/li\u003e\n\u003cli\u003eIvanova S, Gregorc U, Vidergar N, Javier R, Bredt DS, Vandenabeele P, Pardo J, Simon MM, Turk V, Banks L\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eMAGUKs, scaffolding proteins at cell junctions, are substrates of different proteases during apoptosis\u003c/strong\u003e. \u003cem\u003eCell Death Dis \u003c/em\u003e2011, \u003cstrong\u003e2\u003c/strong\u003e:e116.\u003c/li\u003e\n\u003cli\u003eFacciuto F, Cavatorta AL, Valdano MB, Marziali F, Gardiol D: \u003cstrong\u003eDifferential expression of PDZ domain-containing proteins in human diseases - challenging topics and novel issues\u003c/strong\u003e. \u003cem\u003eFEBS J \u003c/em\u003e2012, \u003cstrong\u003e279\u003c/strong\u003e(19):3538-3548.\u003c/li\u003e\n\u003cli\u003eSugihara T, Nakagawa S, Sasajima Y, Ichinose T, Hiraike H, Kondo F, Uozaki H, Fukusato T, Ayabe T: \u003cstrong\u003eLoss of the cell polarity determinant human Discs-large is a novel molecular marker of nodal involvement and poor prognosis in endometrial cancer\u003c/strong\u003e. \u003cem\u003eBr J Cancer \u003c/em\u003e2016, \u003cstrong\u003e114\u003c/strong\u003e(9):1012-1018.\u003c/li\u003e\n\u003cli\u003eAlsafadi S, Houy A, Battistella A, Popova T, Wassef M, Henry E, Tirode F, Constantinou A, Piperno-Neumann S, Roman-Roman S\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eCancer-associated SF3B1 mutations affect alternative splicing by promoting alternative branchpoint usage\u003c/strong\u003e. \u003cem\u003eNat Commun \u003c/em\u003e2016, \u003cstrong\u003e7\u003c/strong\u003e:10615.\u003c/li\u003e\n\u003cli\u003eLiu Z, Yoshimi A, Wang J, Cho H, Chun-Wei Lee S, Ki M, Bitner L, Chu T, Shah H, Liu B\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eMutations in the RNA Splicing Factor SF3B1 Promote Tumorigenesis through MYC Stabilization\u003c/strong\u003e. \u003cem\u003eCancer Discov \u003c/em\u003e2020, \u003cstrong\u003e10\u003c/strong\u003e(6):806-821.\u003c/li\u003e\n\u003cli\u003eLiu B, Liu Z, Chen S, Ki M, Erickson C, Reis-Filho JS, Durham BH, Chang Q, de Stanchina E, Sun Y\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eMutant SF3B1 promotes AKT- and NF-kappaB-driven mammary tumorigenesis\u003c/strong\u003e. \u003cem\u003eJ Clin Invest \u003c/em\u003e2021, \u003cstrong\u003e131\u003c/strong\u003e(1).\u003c/li\u003e\n\u003cli\u003eDeBoever C, Ghia EM, Shepard PJ, Rassenti L, Barrett CL, Jepsen K, Jamieson CH, Carson D, Kipps TJ, Frazer KA: \u003cstrong\u003eTranscriptome sequencing reveals potential mechanism of cryptic 3\u0026apos; splice site selection in SF3B1-mutated cancers\u003c/strong\u003e. \u003cem\u003ePLoS Comput Biol \u003c/em\u003e2015, \u003cstrong\u003e11\u003c/strong\u003e(3):e1004105.\u003c/li\u003e\n\u003cli\u003eHumbert PO, Grzeschik NA, Brumby AM, Galea R, Elsum I, Richardson HE: \u003cstrong\u003eControl of tumourigenesis by the Scribble/Dlg/Lgl polarity module\u003c/strong\u003e. \u003cem\u003eOncogene \u003c/em\u003e2008, \u003cstrong\u003e27\u003c/strong\u003e(55):6888-6907.\u003c/li\u003e\n\u003cli\u003ePeng CY, Manning L, Albertson R, Doe CQ: \u003cstrong\u003eThe tumour-suppressor genes lgl and dlg regulate basal protein targeting in Drosophila neuroblasts\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e2000, \u003cstrong\u003e408\u003c/strong\u003e(6812):596-600.\u003c/li\u003e\n\u003cli\u003eBilder D, Schober M, Perrimon N: \u003cstrong\u003eIntegrated activity of PDZ protein complexes regulates epithelial polarity\u003c/strong\u003e. \u003cem\u003eNat Cell Biol \u003c/em\u003e2003, \u003cstrong\u003e5\u003c/strong\u003e(1):53-58.\u003c/li\u003e\n\u003cli\u003eHawkins ED, Russell SM: \u003cstrong\u003eUpsides and downsides to polarity and asymmetric cell division in leukemia\u003c/strong\u003e. \u003cem\u003eOncogene \u003c/em\u003e2008, \u003cstrong\u003e27\u003c/strong\u003e(55):7003-7017.\u003c/li\u003e\n\u003cli\u003eLamouille S, Xu J, Derynck R: \u003cstrong\u003eMolecular mechanisms of epithelial-mesenchymal transition\u003c/strong\u003e. \u003cem\u003eNat Rev Mol Cell Biol \u003c/em\u003e2014, \u003cstrong\u003e15\u003c/strong\u003e(3):178-196.\u003c/li\u003e\n\u003cli\u003eLan Y, Han J, Wang Y, Wang J, Yang G, Li K, Song R, Zheng T, Liang Y, Pan S\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eSTK17B promotes carcinogenesis and metastasis via AKT/GSK-3beta/Snail signaling in hepatocellular carcinoma\u003c/strong\u003e. \u003cem\u003eCell Death Dis \u003c/em\u003e2018, \u003cstrong\u003e9\u003c/strong\u003e(2):236.\u003c/li\u003e\n\u003cli\u003eStegeman H, Span PN, Kaanders JH, Bussink J: \u003cstrong\u003eImproving chemoradiation efficacy by PI3-K/AKT inhibition\u003c/strong\u003e. \u003cem\u003eCancer Treat Rev \u003c/em\u003e2014, \u003cstrong\u003e40\u003c/strong\u003e(10):1182-1191.\u003c/li\u003e\n\u003cli\u003eZhou BP, Deng J, Xia W, Xu J, Li YM, Gunduz M, Hung MC: \u003cstrong\u003eDual regulation of Snail by GSK-3beta-mediated phosphorylation in control of epithelial-mesenchymal transition\u003c/strong\u003e. \u003cem\u003eNat Cell Biol \u003c/em\u003e2004, \u003cstrong\u003e6\u003c/strong\u003e(10):931-940.\u003c/li\u003e\n\u003cli\u003eZanin-Zhorov A, Lin J, Scher J, Kumari S, Blair D, Hippen KL, Blazar BR, Abramson SB, Lafaille JJ, Dustin ML: \u003cstrong\u003eScaffold protein Disc large homolog 1 is required for T-cell receptor-induced activation of regulatory T-cell function\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2012, \u003cstrong\u003e109\u003c/strong\u003e(5):1625-1630.\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":true,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"journal-of-experimental-and-clinical-cancer-research","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"jecc","sideBox":"Learn more about [Journal of Experimental \u0026 Clinical Cancer Research](http://jeccr.biomedcentral.com)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/jecc/default.aspx","title":"Journal of Experimental \u0026 Clinical Cancer Research","twitterHandle":"@OncoBioMed","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"SF3B1 mutation, Prolactinomas, Alternative splicing, DLG1, Invasion","lastPublishedDoi":"10.21203/rs.3.rs-934219/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-934219/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003ch2\u003eBackground\u003c/h2\u003e \u003cp\u003eRecently, a hotspot mutation in prolactinoma was observed in splicing factor 3b subunit 1 (\u003cem\u003eSF3B1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e), but its functional effects and mechanisms are poorly understood.\u003c/p\u003e\u003ch2\u003eMethods\u003c/h2\u003e \u003cp\u003eUsing the CRISPR/Cas9 genome editing system and rat pituitary GH3 cells, we generated heterozygous \u003cem\u003eSf3b1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e mutant cells. Sanger and whole-genome sequencing were conducted to verify the introduction of this mutation. Transcriptome analysis was performed in \u003cem\u003eSF3B1\u003c/em\u003e-wild-type versus mutant human prolactinoma samples and GH3 cells. Quantitative PCR and minigene reporter assays were conducted to verify aberrant splicing. The functional consequences of \u003cem\u003eSF3B1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e were evaluated \u003cem\u003ein vitro\u003c/em\u003e and \u003cem\u003ein vivo\u003c/em\u003e. Critical makers of epithelial-mesenchymal transition and key components of relevant signaling pathways were detected by western blot, immunohistochemistry, and immunofluorescence, and were knocked down by siRNA-mediated silencing.\u003c/p\u003e\u003ch2\u003eResults\u003c/h2\u003e \u003cp\u003eTranscriptomic analysis of prolactinomas and heterozygous mutant cells revealed that the \u003cem\u003eSF3B1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e allele led to different alterations in splicing properties, affecting different genes in different species. Consistently between rat cells and human tumor samples, mutant SF3B1 promoted aberrant splicing and the suppression of DLG1. Additionally, mutant SF3B1 with knockdown of DLG1 expression promoted cell migration, invasion, and epithelial-mesenchymal transition by activating the PI3K/Akt pathway.\u003c/p\u003e\u003ch2\u003eConclusions\u003c/h2\u003e \u003cp\u003eOur findings elucidate a mechanism through which mutant \u003cem\u003eSF3B1\u003c/em\u003e promotes tumor progression and may provide a potent molecular therapeutic target for prolactinomas with the \u003cem\u003eSF3B1\u003c/em\u003e\u003csup\u003e\u003cem\u003eR625H\u003c/em\u003e\u003c/sup\u003e mutation.\u003c/p\u003e","manuscriptTitle":"The SF3B1R625H Mutation Promotes Prolactinoma Tumor Progression Through Aberrant Splicing Of DLG1","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2021-09-28 19:21:02","doi":"10.21203/rs.3.rs-934219/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Major revision","date":"2021-11-02T03:06:30+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2021-11-01T00:00:00+00:00","index":1,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"editorInvitedReview","content":"","date":"2021-10-08T00:00:00+00:00","index":2,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"reviewerAgreed","content":"","date":"2021-10-03T00:00:00+00:00","index":2,"fulltext":""},{"type":"reviewerAgreed","content":"","date":"2021-10-02T00:00:00+00:00","index":1,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2021-09-24T10:18:05+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2021-09-24T00:00:00+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2021-09-23T23:00:00+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2021-09-23T23:00:00+00:00","index":"","fulltext":""},{"type":"submitted","content":"Journal of Experimental \u0026 Clinical Cancer Research","date":"2021-09-23T12:49:59+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"journal-of-experimental-and-clinical-cancer-research","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"jecc","sideBox":"Learn more about [Journal of Experimental \u0026 Clinical Cancer Research](http://jeccr.biomedcentral.com)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/jecc/default.aspx","title":"Journal of Experimental \u0026 Clinical Cancer Research","twitterHandle":"@OncoBioMed","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"0b4ec4da-c893-4abc-9ec6-66692ef756e8","owner":[],"postedDate":"September 28th, 2021","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":7504107,"name":"Cancer Biology"},{"id":7504108,"name":"Oncology"}],"tags":[],"updatedAt":"2022-01-17T09:21:26+00:00","versionOfRecord":{"articleIdentity":"rs-934219","link":"https://doi.org/10.1186/s13046-022-02245-0","journal":{"identity":"journal-of-experimental-and-clinical-cancer-research","isVorOnly":false,"title":"Journal of Experimental \u0026 Clinical Cancer Research"},"publishedOn":"2022-01-17 09:21:26","publishedOnDateReadable":"January 17th, 2022"},"versionCreatedAt":"2021-09-28 19:21:02","video":"","vorDoi":"10.1186/s13046-022-02245-0","vorDoiUrl":"https://doi.org/10.1186/s13046-022-02245-0","workflowStages":[]},"version":"v1","identity":"rs-934219","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-934219","identity":"rs-934219","version":["v1"]},"buildId":"7rjqhiLT3MXkJMwkYKINL","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}