PDGFRα signaling regulates Srsf3 transcript binding to affect PI3K signaling and endosomal trafficking

  1. Thomas E Forman
  2. Marcin P Sajek
  3. Eric D Larson
  4. Neelanjan Mukherjee
  5. Katherine A Fantauzzo  Is a corresponding author
  1. Department of Craniofacial Biology, School of Dental Medicine, University of Colorado Anschutz Medical Campus, United States
  2. Medical Scientist Training Program, University of Colorado Anschutz Medical Campus, United States
  3. Department of Biochemistry and Molecular Genetics, School of Medicine, University of Colorado Anschutz Medical Campus, United States
  4. RNA Bioscience Initiative, University of Colorado Anschutz Medical Campus, United States
  5. Institute of Human Genetics, Polish Academy of Sciences, Poland
  6. Department of Otolaryngology – Head and Neck Surgery, University of Colorado Anschutz Medical Campus, United States
  7. Basic and Translational Sciences, Penn Dental Medicine, University of Pennsylvania, United States
6 figures, 1 table and 1 additional file

Figures

Figure 1 with 2 supplements
PDGFRα signaling for 1 hr minimally affects gene expression.

(A) Schematic of RNA-seq experimental design. Immortalized mouse embryonic palatal mesenchyme (iMEPM) cells were transduced to stably express a scramble shRNA (scramble) or shRNA targeting the 3’ UTR of Srsf3 (shSrsf3). iMEPM cells expressing either scramble or shSrsf3 were left unstimulated or stimulated with 10 ng/mL PDGF-AA for 1 hr and RNA was isolated for RNA-seq analysis. (B) Western blot (WB) analysis of whole-cell lysates (WCL) from scramble and shSrsf3 cell lines with anti-Srsf3 and anti-Gapdh antibodies. The percentage of Srsf3 expression normalized to Gapdh expression is indicated below. (C) Volcano plots depicting differentially expressed genes in scramble versus shSrsf3 cell lines in the absence (left) or presence (right) of PDGF-AA stimulation. Log2 (fold change) (FC) values represent log2 (shSrsf3 normalized counts/scramble normalized counts). Significant changes in gene-level expression are reported for genes with adjusted p (padj)<0.05 and fold change |FC| ≥ 2. (D) Venn diagram of significant genes from (C). (E) Volcano plots depicting differentially expressed genes in the absence versus presence of PDGF-AA ligand in scramble (left) or shSrsf3 (right) cell lines. Log2 (FC) values represent log2 (+PDGF-AA normalized counts/-PDGF-AA normalized counts). (F) Venn diagram of significant genes from (E).

Figure 1—figure supplement 1
High correlation of Srsf3-dependent differentially expressed genes across ligand treatment conditions.

(A, B) Scatter dot plots depicting Srsf3-dependent (A) and PDGF-AA-dependent (B) differentially expressed genes. Log2 (fold change) (FC) values represent log2 (shSrsf3 normalized counts/scramble normalized counts) (A) or log2 (+PDGF-AA normalized counts/-PDGF-AA normalized counts) (B). Spearman correlation values and approximate p-values are listed. Immediate early genes are represented in red in (B).

Figure 1—figure supplement 2
Gene ontology (GO) analysis of differentially expressed genes across treatment comparisons.

(A, B) Bubble plots depicting up to 10 of the most significant GO terms for biological process (A) and molecular function (B) for Srsf3-dependent and PDGF-AA-dependent differentially expressed genes. Colors correspond to -log10 (adjusted p-value); sizes correspond to the number of genes.

Figure 2 with 2 supplements
PDGFRα signaling for 1 hr has a more pronounced effect on alternative RNA splicing.

(A) Volcano plots depicting alternatively spliced transcripts in scramble versus shSrsf3 cell lines in the absence (left) or presence (right) of PDGF-AA stimulation. Difference in percent spliced in (ΔPSI) values represent scramble PSI – shSrsf3 PSI. Significant changes in alternative RNA splicing are reported for events with a false discovery rate (FDR) ≤ 0.05 and a difference in percent spliced in (|ΔPSI|) ≥ 0.05. (B) Venn diagram of significant transcripts from (A), filtered to include events detected in at least 10 reads in either condition. (C) Volcano plots depicting alternatively spliced transcripts in the absence versus presence of PDGF-AA ligand in scramble (left) or shSrsf3 (right) cell lines. Difference in percent spliced in (ΔPSI) values represent -PDGF-AA PSI –+PDGF-AA PSI. (D) Venn diagram of significant transcripts from (C), filtered to include events detected in at least 10 reads in either condition. (E) Bar graphs depicting alternative RNA splicing events in scramble versus shSrsf3 cell lines in the absence or presence of PDGF-AA stimulation (left) or in the absence versus presence of PDGF-AA ligand in scramble or shSrsf3 cell lines (right).

Figure 2—source data 1

rMATS output for scramble (-PDGF-AA) versus shSrsf3 (-PDGF-AA) RNA-seq analysis.

https://cdn.elifesciences.org/articles/98531/elife-98531-fig2-data1-v1.xlsx
Figure 2—source data 2

rMATS output for scramble (+PDGF-AA) versus shSrsf3 (+PDGF-AA) RNA-seq analysis.

https://cdn.elifesciences.org/articles/98531/elife-98531-fig2-data2-v1.xlsx
Figure 2—source data 3

rMATS output for -PDGF-AA (scramble) versus +PDGF-AA (scramble) RNA-seq analysis.

https://cdn.elifesciences.org/articles/98531/elife-98531-fig2-data3-v1.xlsx
Figure 2—source data 4

rMATS output for -PDGF-AA (shSrsf3) versus +PDGF-AA (shSrsf3) RNA-seq analysis.

https://cdn.elifesciences.org/articles/98531/elife-98531-fig2-data4-v1.xlsx
Figure 2—figure supplement 1
qPCR validation of differential AS between scramble and shSrsf3 samples.

(A, B) Representative qPCR gels (left) with depictions of differentially alternatively spliced exon (gray), and upstream and downstream sequences (white) that were assessed by qPCR in scramble (sc) versus shSrsf3 (sh) samples for Arhgap12 (A) and Cep55 (B). Scatter dot plots (right) depicting the percent spliced in from n = 3 biological replicates as at left. Data are mean ± s.e.m. *p<0.05 (two-tailed, unpaired t-test with Welch’s correction). Shaded circles correspond to independent experiments.

Figure 2—figure supplement 2
Gene ontology (GO) analysis of alternatively spliced transcripts across treatment comparisons.

(A, B) Bubble plots depicting up to 10 of the most significant GO terms for biological process (A) and molecular function (B) for Srsf3-dependent and PDGF-AA-dependent alternatively spliced transcripts. Colors correspond to -log10 (adjusted p-value); sizes correspond to number of genes.

Figure 3 with 1 supplement
Srsf3 exhibits differential transcript binding upon PDGFRα signaling.

(A) Schematic of enhanced UV-crosslinking and immunoprecipitation (eCLIP) experimental design. Immortalized mouse embryonic palatal mesenchyme (iMEPM) cells were left unstimulated or stimulated with 10 ng/mL PDGF-AA for 1 hr and processed for eCLIP analysis. (B) Immunoprecipitation (IP) of Srsf3 from cells that were UV-crosslinked or not UV-crosslinked with IgG or an anti-Srsf3 antibody followed by western blotting (WB) of input, supernatant (Sup), and IP samples with an anti-Srsf3 antibody. (C) Mapping of eCLIP peaks to various transcript locations in the absence or presence of PDGF-AA stimulation. 5’ UTR, 5’ untranslated region; CDS, coding sequence; 3’ UTR, 3’ untranslated region. (D, E) Mean coverage of eCLIP peaks across various transcript locations (D) and surrounding the 5’ and 3’ splice sites (E) in the absence or presence of PDGF-AA stimulation. (F, G) Top three motifs enriched in eCLIP peaks in the absence (F) or presence (G) of PDGF-AA stimulation with associated p-values (t-test).

Figure 3—figure supplement 1
PDGFRα signaling influences Srsf3 binding specificity.

Top 10 motifs enriched in enhanced UV-crosslinking and immunoprecipitation (eCLIP) peaks in the absence (left) or presence (right) of PDGF-AA stimulation with associated p-values (t-test).

Srsf3 and PDGFRα signaling are associated with differential GC content and length of alternatively spliced exons.

(A) Box and whisker plot depicting the percentage of exon GC content in exons that are not differentially alternatively spliced, and exons that are included or skipped when Srsf3 is present from the rMATS analysis. (B) Box and whisker plot depicting the ratio of downstream intron to exon GC content in exons that are not differentially alternatively spliced, and exons that are included or skipped when Srsf3 is present from the rMATS analysis. (C, D) Box and whisker plots depicting the ratio of upstream intron to exon length (C) and downstream intron to exon length (D) in exons that are not differentially alternatively spliced, and exons that are included or skipped when Srsf3 is present from the rMATS analysis. (E) Violin and box and whisker (inset) plots depicting the percentage of exon GC content in exons that are not bound by Srsf3, and exons that are bound in the absence and/or presence of PDGF-AA stimulation from the enhanced UV-crosslinking and immunoprecipitation (eCLIP) analysis. *p<0.05; **p<0.01; ***p<0.001 (Mann–Whitney U-test).

Figure 5 with 1 supplement
Transcripts bound by Srsf3 that undergo alternative splicing upon PDGFRα signaling encode regulators of PI3K signaling.

(A) Venn diagram of genes with differential expression (DE) or transcripts subject to alternative RNA splicing (AS) across the four treatment comparisons that overlap with transcripts with Srsf3 enhanced UV-crosslinking and immunoprecipitation (eCLIP) peaks in the absence or presence of PDGF-AA stimulation. (B, C) Top 10 (B) and PI3K-related (C) biological process gene ontology (GO) terms for transcripts from the high-confidence, overlapping datasets. p.val, p. (D) Difference in percent spliced in (ΔPSI) values for PI3K/endosome-related transcripts of interest. ΔPSI values represent -PDGF-AA PSI –+PDGF-AA PSI. FDR, false detection rate. (E) Peak visualization for input and eCLIP samples in the absence or presence of PDGF-AA stimulation from Integrative Genomics Viewer (left) with location of motifs from Figure 5—figure supplement 1 indicated below for PI3K/endosome-related transcripts of interest. Predicted alternative splicing outcomes for PI3K/endosome-related transcripts of interest (right). (F) Scatter dot plot depicting the percent spliced in as assessed by qPCR analysis of Wdr81 exon 9 in the scramble cell line following 1 hr of PDGF-AA stimulation quantified from n = 3 biological replicates. Data are mean ± s.e.m. Shaded circles correspond to independent experiments. (G) Western blot (WB) analysis of whole-cell lysates (WCL) from the scramble cell line following 24 hr of PDGF-AA stimulation with anti-Wdr81 and anti-Gapdh antibodies. Bar graph depicting relative ratios of larger/smaller Wdr81 isoforms quantified from n = 3 biological replicates as above. Data are mean ± s.e.m. p=0.052 (two-tailed, ratio paired t-test).

Figure 5—source data 1

List of transcripts and genes from Venn diagram in Figure 5A.

https://cdn.elifesciences.org/articles/98531/elife-98531-fig5-data1-v1.xlsx
Figure 5—source data 2

High-confidence, overlapping dataset output correlating eCLIP with scramble (-PDGF-AA) versus shSrsf3 (-PDGF-AA) rMATS RNA-seq analysis.

https://cdn.elifesciences.org/articles/98531/elife-98531-fig5-data2-v1.xlsx
Figure 5—source data 3

High-confidence, overlapping dataset output correlating eCLIP with scramble (+PDGF-AA) versus shSrsf3 (+PDGF-AA) rMATS RNA-seq analysis.

https://cdn.elifesciences.org/articles/98531/elife-98531-fig5-data3-v1.xlsx
Figure 5—source data 4

High-confidence, overlapping dataset output correlating eCLIP with -PDGF-AA (scramble) versus +PDGF-AA (scramble) rMATS RNA-seq analysis.

https://cdn.elifesciences.org/articles/98531/elife-98531-fig5-data4-v1.xlsx
Figure 5—source data 5

High-confidence, overlapping dataset output correlating eCLIP with -PDGF-AA (shSrsf3) versus +PDGF-AA (shSrsf3) rMATS RNA-seq analysis.

https://cdn.elifesciences.org/articles/98531/elife-98531-fig5-data5-v1.xlsx
Figure 5—source data 6

Wdr81 western blots.

https://cdn.elifesciences.org/articles/98531/elife-98531-fig5-data6-v1.zip
Figure 5—source data 7

Gapdh western blots.

https://cdn.elifesciences.org/articles/98531/elife-98531-fig5-data7-v1.zip
Figure 5—source data 8

Relative ratios of larger/smaller Wdr81 isoforms.

https://cdn.elifesciences.org/articles/98531/elife-98531-fig5-data8-v1.xlsx
Figure 5—figure supplement 1
Srsf3 exhibits differential transcript binding upon PDGFRα signaling in the subset of transcripts from the high-confidence, overlapping datasets.

(A, B) Mean coverage of enhanced UV-crosslinking and immunoprecipitation (eCLIP) peaks within the high-confidence, overlapping datasets across various transcript locations (A) and surrounding the 5’ and 3’ splice sites (B) in the absence or presence of PDGF-AA stimulation. (C, D) Top three motifs enriched in eCLIP peaks within the high-confidence, overlapping datasets in the absence (C) or presence (D) of PDGF-AA stimulation.

Srsf3 regulates early endosome size and phosphorylation of Akt downstream of PDGFRα signaling.

(A, B) Scatter dot plots depicting average size of Rab5 puncta per cell (A) and Pearson’s correlation coefficient of signals from anti-Rab5 and anti-PDGFRα antibodies (B) in scramble and shSrsf3 cell lines in the absence or presence (15–60 min) of PDGF-AA stimulation. Data are mean ± s.e.m. *p<0.05 (two-way ANOVA followed by uncorrected Fisher’s LSD test). Shaded shapes correspond to independent experiments. Summary statistics from biological replicates consisting of independent experiments (large shapes) are superimposed on top of data from all cells; n = 20 technical replicates across each of three biological replicates. (C–H”’) PDGFRα antibody signal (white or magenta) and Rab5 antibody signal (white or green) as assessed by immunofluorescence analysis of scramble and shSrsf3 cells in the absence or presence (15–60 min) of PDGF-AA stimulation. Nuclei were stained with DAPI (blue). White arrows denote regions of colocalization, which are expanded in (C’’’–H’’’). Scale bars: 20 μm (C–H’’), 3 μm (C’’’–H’’’). (I) Western blot (WB) analysis of whole-cell lysates (WCL) from scramble (left) and shSrsf3 (right) cell lines following a time course of PDGF-AA stimulation from 15 min to 4 hr, with anti-phospho-(p)-Akt and anti-Akt antibodies. Line graph depicting quantification of band intensities from n = 4 biological replicates as above. Data are mean ± s.e.m. *p<0.05; **p<0.01 (two-tailed, ratio paired t-test within each cell line and a two-tailed, unpaired t-test with Welch’s correction between each cell line). (J) Model of experimental results in which PI3K/Akt-mediated PDGFRα signaling results in the nuclear translocation of Srsf3 and the subsequent AS of transcripts to decrease levels of proteins that promote PDGFRα trafficking out of early endosomes.

Tables

Appendix 1—key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Cell line (Mus musculus, male)iMEPMFantauzzo and Soriano, 2017
Cell line (M. musculus)Scramble iMEPMThis paperCell line maintained in K. Fantauzzo lab
Cell line (M. musculus)Srsf3 shRNA iMEPMThis paperCell line maintained in K. Fantauzzo lab
AntibodySrsf3 (rabbit polyclonal)AbcamAbcam:ab73891IP (10 μg/sample), WB (1:1000)
AntibodyGapdh (mouse monoclonal)Proteintech Group, IncProteintech:60004-1-IgWB (1:50,000)
AntibodyWdr81 (rabbit polyclonal)Proteintech Group, IncProteintech:24874-1-APWB (1:1000)
AntibodyPhospho-Akt (rabbit polyclonal)Cell Signaling TechnologyCell Signaling Technology:9271SWB (1:1000)
AntibodyAkt (rabbit polyclonal)Cell Signaling TechnologyCell Signaling Technology:9272SWB (1:1000)
AntibodyPeroxidase AffiniPure Goat Anti-Mouse IgG (H+L) (goat polyclonal)Jackson ImmunoResearch IncJackson ImmunoResearch:115-035-003WB (1:20,000)
AntibodyPeroxidase AffiniPure Goat Anti-Rabbit IgG (H+L) (goat polyclonal)Jackson ImmunoResearch IncJackson ImmunoResearch:111-035-003WB (1:20,000)
AntibodyRab5 (rabbit monoclonal)Cell Signaling TechnologyCell Signaling Technology:3547SIF (1:200)
AntibodyPDGFRα (mouse polyclonal)R&D SystemsR&D Systems:AF1062IF (1:20)
AntibodyDonkey anti-Rabbit IgG (H+L) Highly Cross-Absorbed Secondary Antibody, Alexa Fluor 488 (donkey polyclonal)InvitrogenInvitrogen:A21206IF (1:1000)
AntibodyDonkey anti-Mouse IgG (H+L) Highly Cross-Absorbed Secondary Antibody, Alexa Fluor 546 (donkey polyclonal)InvitrogenInvitrogen:A10036IF (1:1000)
Recombinant DNA reagentpLV[shRNA]-EGFP:T2A:Puro-U6>Scramble_shRNAVectorBuilderVector Builder:VB010000-0009mxcCCTAAGGTTAAGTCGCCCTCGCTCGAGCGAGGGCGACTTAACCTTAGG
Recombinant DNA reagentpLV[shRNA]-EGFP:T2A:Puro-U6>mSrsf3[shRNA#1]VectorBuilderVector Builder:VB90060-7699yyhGAATGATAAAGCGGTGTTTACTCGAGTAAACACCGCTTTATCATTCC
Sequence-based reagentRandom primersInvitrogenInvitrogen:48190011
Sequence-based reagentoligo (dT)12-18 primerInvitrogenInvitrogen:18418012
Sequence-based reagentArhgap12_FThis paperPCR primerGGAGACATAGCACCATTGTG
Sequence-based reagentArhgap12_RThis paperPCR primerGCACTGCCCAAGAAGACAAC
Sequence-based reagentCep55_FThis paperPCR primerCCTTTCGGCTCCTTTGAACT
Sequence-based reagentCep55_RThis paperPCR primerGCAGTGTCTGACTTGGAGCT
Sequence-based reagentWdr81_FThis paperPCR primerGCTTTGTGGACTGCAGGAAG
Sequence-based reagentWdr81_RThis paperPCR primerGCAGGGAACAGACACCAAT
Sequence-based reagentRiL19This paperBarcoded RNA adapters/5phos/rArGrArUrCrGrGrArArGrArGrCrGrUrCrGrUrG/3SpC3/
Sequence-based reagentRand103Tr3This paperBarcoded DNA adapters/5Phos/NNNNNNNNNNAGATCGGAAGAGCACACGTCTG/3SpC3/
Peptide, recombinant proteinPDGF-AAR&D SystemsR&D:1055AA05010 ng/mL
Peptide, recombinant proteinProtein A DynabeadsInvitrogenInvitrogen:10002DIP (125 μL/sample)
Peptide, recombinant proteinDynabeads M-280 sheep anti-rabbit IgGInvitrogenInvitrogen:11203DIP (125 μL/sample)
Peptide, recombinant proteinSuperScript II RTInvitrogenInvitrogen:18064014
Peptide, recombinant proteinRNase IInvitrogenInvitrogen:AM2294
Peptide, recombinant proteinFastAPThermo Fisher ScientificThermo Fisher Scientific:EF0654
Peptide, recombinant proteinT4 PNKNew England BiolabsNew England Biolabs:M0201S
Peptide, recombinant proteinT4 RNA ligaseNew England BiolabsNew England Biolabs:M0437M
Peptide, recombinant proteinProteinase KNew England BiolabsNew England Biolabs:P8107S
Peptide, recombinant proteinSuperScript IVInvitrogenInvitrogen:18090010
Chemical compound, drugDAPISigma-AldrichSigma-Aldrich:D9542-10MG2 μg/mL
Commercial assay or kitMycoAlert Mycoplasma Detection KitLonza Group LtdLonza:LT07-218
Commercial assay or kitRNeasy Mini KitQIAGEN, IncQIAGEN:74104
Commercial assay or kitUniversal Plus mRNA-Seq kit with NuQuantTecan Group LtdTecan Group:0361-A01
Commercial assay or kitAMPure XP ReagentBeckman CoulterBeckman Coulter:A63880
Commercial assay or kitExoSAP-ITApplied BiosystemsApplied Biosystems:78200.200.UL
Commercial assay or kitDynabeads MyOne SilaneInvitrogenInvitrogen:37002D
Commercial assay or kitQ5 PCR mixNew England BiolabsNew England Biolabs:M0492S
Software, algorithmImageJNIHRRID:SCR_003070
Software, algorithmPrism 10GraphPad Software Inc.RRID:SCR_000306
Software, algorithmbcl2fastqIlluminaRRID:SCR_015058
Software, algorithmBBDukhttps://sourceforge.net/projects/bbmap/
Software, algorithmSalmonPatro et al., 2017
Software, algorithmGENCODEFrankish et al., 2019
Software, algorithmtximportSoneson et al., 2015
Software, algorithmDESeq2Love et al., 2014
Software, algorithmSTARDobin et al., 2013
Software, algorithmrMATSShen et al., 2014
Software, algorithmEnrichrChen et al., 2013; Kuleshov et al., 2016
Software, algorithmcutadaptMartin, 2011
Software, algorithmFASTX-Toolkithttp://hannonlab.cshl.edu/fastx_toolkit
Software, algorithmseqtkhttps://github.com/lh3/seqtk; Li, 2018
Software, algorithmmsRepDBLiao et al., 2022
Software, algorithmomniCLIPDrewe-Boss et al., 2018
Software, algorithmRCASUyar et al., 2017
Software, algorithmdeepToolsRamírez et al., 2016
Software, algorithmIntegrative Genomics ViewerRobinson et al., 2011
Software, algorithmMattGohr and Irimia, 2019
Software, algorithmrMAPSPark et al., 2016
Software, algorithmvalrRiemondy et al., 2017
Software, algorithmZEN BlueCarl Zeiss Microscopy LLCRRID:SCR_013672
OtherNormal donkey serumJackson ImmunoResearch IncJackson ImmunoResearch:017-000-121
OtherVECTASHIELD HardSet Antifade Mounting MediumVector Laboratories, IncVector Laboratories:H-1400–10

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Thomas E Forman
  2. Marcin P Sajek
  3. Eric D Larson
  4. Neelanjan Mukherjee
  5. Katherine A Fantauzzo
(2024)
PDGFRα signaling regulates Srsf3 transcript binding to affect PI3K signaling and endosomal trafficking
eLife 13:RP98531.
https://doi.org/10.7554/eLife.98531.3