Next Article in Journal
Antiviral Effects of ABMA and DABMA against Influenza Virus In Vitro and In Vivo via Regulating the Endolysosomal Pathway and Autophagy
Previous Article in Journal
The Characterization of the Phloem Protein 2 Gene Family Associated with Resistance to Sclerotinia sclerotiorum in Brassica napus
Previous Article in Special Issue
An Influence of Modification with Phosphoryl Guanidine Combined with a 2′-O-Methyl or 2′-Fluoro Group on the Small-Interfering-RNA Effect
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 7
10.3390/ijms23073937
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Antisense Oligonucleotide Induction of the hnRNPA1b Isoform Affects Pre-mRNA Splicing of SMN2 in SMA Type I Fibroblasts

by
Jarichad Toosaranont
1,
Sukanya Ruschadaariyachat
1,
Warasinee Mujchariyakul
2,
Jantarika Kumar Arora
2,
Varodom Charoensawan
2,3,4,
Bhoom Suktitipat
1,3,
Thomas N. Palmer
5,
Sue Fletcher
5,6,
Steve D. Wilton
5,6 and
Chalermchai Mitrpant
1,5,6,*
1
Department of Biochemistry, Faculty of Medicine, Siriraj Hospital, Mahidol University, Bangkok 10700, Thailand
2
Department of Biochemistry, Faculty of Science, Mahidol University, Bangkok 10700, Thailand
3
Integrative Computational BioScience (ICBS) Center, Mahidol University, Nakhon Pathom 73170, Thailand
4
Systems Biology of Diseases Research Unit, Faculty of Science, Mahidol University, Bangkok 10400, Thailand
5
Perron Institute for Neurological and Translational Science, The University of Western Australia, Perth, WA 6009, Australia
6
Centre for Molecular Medicine and Innovative Therapeutics, Health Futures Institute, Murdoch University, Perth, WA 6150, Australia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(7), 3937; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms23073937
Submission received: 28 December 2021 / Revised: 22 March 2022 / Accepted: 28 March 2022 / Published: 1 April 2022
(This article belongs to the Special Issue Precision Nucleic Acid Therapeutics)
Browse Figures
Versions Notes

Abstract

:
Spinal muscular atrophy (SMA) is a severe, debilitating neuromuscular condition characterised by loss of motor neurons and progressive muscle wasting. SMA is caused by a loss of expression of SMN1 that encodes the survival motor neuron (SMN) protein necessary for the survival of motor neurons. Restoration of SMN expression through increased inclusion of SMN2 exon 7 is known to ameliorate symptoms in SMA patients. As a consequence, regulation of pre-mRNA splicing of SMN2 could provide a potential molecular therapy for SMA. In this study, we explored if splice switching antisense oligonucleotides could redirect the splicing repressor hnRNPA1 to the hnRNPA1b isoform and restore SMN expression in fibroblasts from a type I SMA patient. Antisense oligonucleotides (AOs) were designed to promote exon 7b retention in the mature mRNA and induce the hnRNPA1b isoform. RT-PCR and western blot analysis were used to assess and monitor the efficiency of different AO combinations. A combination of AOs targeting multiple silencing motifs in hnRNPA1 pre-mRNA led to robust hnRNPA1b induction, which, in turn, significantly increased expression of full-length SMN (FL-SMN) protein. A combination of PMOs targeting the same motifs also strongly induced hnRNPA1b isoform, but surprisingly SMN2 exon 5 skipping was detected, and the PMO cocktail did not lead to a significant increase in expression of FL-SMN protein. We further performed RNA sequencing to assess the genome-wide effects of hnRNPA1b induction. Some 3244 genes were differentially expressed between the hnRNPA1b-induced and untreated SMA fibroblasts, which are functionally enriched in cell cycle and chromosome segregation processes. RT-PCR analysis demonstrated that expression of the master regulator of these enrichment pathways, MYBL2 and FOXM1B, were reduced in response to PMO treatment. These findings suggested that induction of hnRNPA1b can promote SMN protein expression, but not at sufficient levels to be clinically relevant.

1. Introductions

Spinal muscular atrophy (SMA), a severe, debilitating neuromuscular condition, is the leading genetic cause of death in children. SMA is caused by insufficient expression of survival motor neuron (SMN) protein, most commonly stemming from genomic deletions of the entire SMN1 gene. The complete absence of the SMN protein is incompatible with life, as established by the unviability of Smn-knockout models. The SMA condition manifests in humans due to the presence of a near identical gene, SMN2, that could potentially encode a normal SMN but which is inappropriate for splicing of the majority of SMN2 transcripts, and hence, low expression of the full-length SMN (FL-SMN). The aberrant splicing of SMN2 pre-mRNA leads mainly to mature SMN2 mRNA transcripts missing exon 7 that are translated into a nonfunctional protein. hnRNPA1 was identified as a crucial splice repressor regulating inclusion of exon 7 in the SMN2 transcript [1,2,3,4]. Depletion of the hnRNPA1 expression was shown to promote inclusion of SMN2 exon 7 into mature transcript, and therefore, an increase in SMN protein [5,6]. Antisense oligonucleotides (AOs), synthetic short single-stranded nucleotide with chemical modifications on the sugar backbone of AOs, were designed to target the binding site of this repressor (ISS-N1) [2]. Treatment of this AO efficiently increased SMN protein in cultured fibroblast from a patient with SMA [7] and animal models [8,9]. This compound (Spinraza®) is now commercially available for SMA patients. The development of this treatment has highlighted the fact that AO mediated splice switching is a viable molecular treatment option for patients with SMA [8,10]. The class of the AO backbone to produce Spinraza® is the methoxyethyl modified AO chemistry; this class of backbone chemistry was previously associated with proinflammatory effect, thrombocythemia or kidney complications in some of the clinical trials [11,12,13]. Herein, we propose the use of AO to alter the splice repressor hnRNPA1 to another isoform as a means to induce inclusion of SMN2 exon 7 into mature transcript, in case Spinraza® causes major adverse effects later on.
hnRNPA1 encodes two main isoforms: hnRNPA1 (320 amino acid) and hnRNPA1b (372 amino acid). The hnRNPA1 isoform is predominantly expressed in most tissues (GTEx portal). The hnRNPA1b transcript is an alternatively spliced isoform with the inclusion of exon 7b encoding a protein with 52 additional amino acids in a glycine/arginine rich region, a domain that is reportedly important in protein–protein interactions [14]. Interestingly, the overexpression of hnRNPA1b is reported to yield an alternative 5′ splice site selection and overexpression of hnRNPA1 protein, as determined by in vitro splicing assay [15].
This study investigates whether AO-mediated shifting of the hnRNPA1 isoforms towards hnRNPA1b expression can alter the splicing pattern of SMN2 and augment SMN protein expression. We employed an RNA sequencing approach to determine the effects and consequences of this hnRNPA1 isoform shifting on the transcriptome in SMA patient fibroblasts.

2. Results

2.1. Individual AO Induction of hnRNPA1 Exon 7b Inclusion

In this study, ten 2′ O-methyl modified AO candidates were evaluated for their effectiveness in modifying the hnRNPA1 splicing pattern. The panel of ten AOs was designed to target either the predicted ESE motif (AOs 1 and 2) or previously reported ISSs (AOs 3–10) (Figure 1) [16]. The splicing pattern of hnRNPA1 transcripts was determined by RT-PCR analysis, with a typical gel image displayed in the upper panel of Figure 2. A densitometric analysis of hnRNPA1 isoforms was performed, and is displayed as mean percentages of hnRNPA1 exon 7b inclusion (Figure 2B). AOs 1 and 2 were designed to mask predicted exonic splice enhancers (ESE) and, as expected, this led to substantial skipping of exon 7b with reduced overall amount of hnRNPA1 transcript. Three AOs (AOs 4, 7 and 8) in particular promoted an approximately 40% increase in hnRNPA1 exon 7b inclusion at all tested concentrations.

2.2. Systematic Assessment of Different AO Combinations Targeting Multiple Silencing Motifs to Induce hnRNPA1b Exon 7b Inclusion

We explored whether combinations of AOs could improve the efficiency of hnRNPA1b induction, as “looping out” had been proposed as a mechanism for hnRNPA1 exon 7b skipping [15]. This mechanism indicates that hnRNPA1 protein binds to splicing motifs located in adjacent introns of hnRNPA1 exon 7b to induce exon 7b skipping. This led us to focus on using AOs to sterically and simultaneously block these conserved elements (CE) by using combinations of AOs. One AO, AO9 targets a motif in intron 7, ISS(CE1A), whereas AOs 8, 7 and 4, targeting the other three ISS motifs, i.e., ISS(CE4), ISS(CE6) or ISS(CE9) respectively, were chosen based upon the efficiency of single AO treatment (Figure 1). We then combined AO9, that masks ISS(CE1A), with AOs targeting the other three motifs (viz. ISS(CE4), ISS(CE6) or ISS(CE9)), giving a total of seven possible combinations (see Table 1).
These combinations were then tested at three different concentrations, i.e., 500 nM, 250 nM or 125 nM, in three separate experiments, and a RT-PCR analysis of hnRNPA1 and SMN splice isoforms was undertaken to test the efficiency of the various AO cocktails (Figure 3). The AOs9 + 8 + 7, AOs9 + 8 + 4 and AOs9 + 8 + 7 + 4 combination yielded approximately 10–15% more of the hnRNPA1b transcript compared to the two AO cocktails targeting a similar motif (viz. AOs9 + 8 or AOs9 + 7). Moreover, increased levels of the hnRNPA1b transcript with higher concentrations of AOs9 + 8 + 7 correlated with changes in inclusion of SMN2 exon 7. AOs9 + 4 led to decreased overall amount of hnRNPA1 transcripts with substantially increased inclusion of SMN2 exon 7, despite only 5–10% of hnRNPA1 exon 7b being induced (Figure 3).

2.3. AO Combinations Targeting Multiple Silencing Motifs as a Means of Inducing hnRNPA1b Protein

Transfection with each AO cocktail at a concentration of 500 nM was examined to establish whether induced transcript levels are associated with increased protein expression. Most combinations significantly increased hnRNPA1b expression with slight increase in the overall amount of hnRNPA1 transcript. Levels of induction ranged from 20% to 35%, compared to 5% of hnRNPA1b expression in untreated and sham transfected fibroblasts. The AOs9 + 8 + 7 cocktail also significantly increased FL-SMN expression in SMA fibroblasts compared to sham treated control (Figure 4).
We synthesised AOs 9, 4, 7, and 8 as phosphorodiamidate morpholino oligomers to further enhance efficiency of exon skipping, and the different backbone chemistries showed greater biological stability. Most combinations, except for PMOs9 + 4, strongly induced 50 to 60% of hnRNPA1b protein in a cultured SMA fibroblast. PMO cocktails that contain the PMO 7 in the mixture, i.e., PMOs9 + 7, PMOs9 + 8 + 7, PMOs9 + 7 + 4, and PMOs9 + 8 + 7 + 4, caused skipping of SMN2 exon 5 (Figure 5C). None of the PMO cocktails evaluated was found to augment FL-SMN protein expression (Figure 5).

2.4. Transcriptomic Analysis of PMOs 9 + 7 + 8 Induced hnRNPA1b Isoform in SMA Fibroblasts

We further investigated the effect of PMOs9 + 7 + 8 induced hnRNPA1b expression at the transcriptomics level in cultured SMA fibroblast. A panel of differentially expressed genes between the combined PMOs9 + 7 + 8 treated and untreated SMA fibroblasts was detected and is displayed as a volcano plot in Figure 6A. A panel of differentially expressed 3244 genes were identified, as the absolute log2 fold change was greater than 1.5 with p-value lesser than 0.05. Expression of Survivin (BIRC5), Aurora kinase B (AURKB) and marker of proliferation ki-67 (MKI67) were approximately four-fold less in the PMOs9 + 7 + 8 treated SMA fibroblasts compared to untreated fibroblast with the adjusted p-value of 1.52 × 10−24, 3.58 × 10−6 or 2.91 × 10−6, respectively (Supplementary Table S3). A gene ontology analysis of these genes indicated that these differentially expressed genes involved the cycle control, mitosis, nuclear division and chromosome segregation (Figure 6B,C). We further performed a gene set enrichment analysis (GSEA) on this gene set and found that these genes were enriched in the cell cycle process, with a normalised enrichment score of −1.8477 and a FDR q-value of 0.00 (Figure 6D). FOXM1 and MYBL2 are strong candidates for regulating the cell cycle, nuclear division and chromosome segregation controlling protein, as either FOXM1 or MYBL2 was present in nine out of twenty enriched pathways identified by the GSEA analysis. FOXM1 has four isoforms, i.e., FOXM1A, FOXM1B, FOXM1C and FOXM1D, and the isoform specific RT-PCR indicated that only FOXM1B was reduced by transfection with PMOs9 + 7 + 8, but not with PMOs 9 + 4 or in untreated SMA fibroblasts, while MYBL2 was also reduced in SMA fibroblasts treated with this PMO cocktail (Figure 7).

3. Discussion

Regulation of SMN2 exon 7 inclusion has become a most promising focus as a molecular therapy for SMA patients. Despite potentially encoding an identical SMN protein, the predominant SMN2 gene transcripts in SMA are missing exon 7 and encode a defective unstable protein. Shifting the SMN2 splicing pattern towards the normal SMN full-length product has been proven to be effective in increasing SMN protein levels with concomitant improvements in phenotype [2]. Different strategies (viz. small molecules and AOs) have also been utilised to modify the process of SMN2 exon 7 inclusion [9,17,18,19].
SF2/ASF (SRSF1) and hnRNPA1 are two crucial splicing factors that are responsible for control of differential splicing of the SMN1 and SMN2 gene transcripts. SF2/ASF binds to an ESE located in SMN1 exon 7, [5,20] while a C>T polymorphism at the sixth nucleotide of SMN2 exon 7 creates an ESS that is recognised by hnRNPA1, thus leading to omission of exon 7 from the mature SMN2 mRNA [3,4]. In addition, one silencer element located in SMN2 intron 7 (ISS-N1) that has been used as a target for FL-SMN induction exerts its SMN2 exon 7 skipping effect through its binding to hnRNPA1 [3,4].
As a splice repressor, hnRNPA1 binds to pre-mRNA through two RNA recognition motifs (RRMs) located in the amino terminal region of hnRNPA1 protein, [16] and the glycine/arginine rich carboxyl terminal region is reported to be critical in protein–protein interactions and pre-mRNA splicing [15]. hnRNPA1b contains an extra 52 amino acids inserted into the glycine/arginine rich carboxyl terminal region, and overexpression of hnRNPA1b is known to lead to alternative 5′ splice site selection rather than overexpression of the full-length isoform [21].
Eight AOs targeting four silencing motifs [16] were designed and tested individually for their effectiveness in inducing hnRNPA1b exon 7b inclusion. Transfection of SMA fibroblasts with individual AOs raised hnRNPA1b exon 7b inclusion by approximately 20–40% based upon semiquantitative RT-PCR analysis (Figure 2). As the proposed model for hnRNPA1 exon 7b exclusion involves simultaneous binding of splice repressor to multiple splice silencers, we explored the relative effectiveness of combinations of AOs targeting four different silencing motifs. We combined AO9, which targets the proximal CE1A silencer, with other AOs masking either the CE4 (AO8), CE6 (AO7) and CE9 (AO 4) silencers. Interaction of CE1A and CE4 motifs was demonstrated to facilitate the exclusion of hnRNPA1 exon 7b, [16] and our AO cocktails (9 + 8) may destabilise the hnRNPA1 pre-mRNA loop, promoting inclusion of hnRNPA1 exon 7b and thereby increasing production of the hnRNPa1b isoform.
Seven possible combinations were tested for hnRNPA1b induction (Table 1). These induced 35–40% hnRNPA1b exon 7b inclusion, as assessed at the protein level, and the AOs 9 + 8 + 7 combination was the most effective AO cocktail at inducing hnRNPA1b protein expression (Figure 4). Increased hnRNPA1b expression was associated with inclusion of exon 7 in up to 80% of the SMN2 transcripts (Figure 3), and this effect translated into significantly increased FL-SMN protein expression (Figure 4).
Antisense oligomers with the uncharged morpholino backbone have proved to be promising for clinical use as steric blockers to redirect splicing, as these compounds possesses excellence safety profiles and are not degraded by known nucleolytic enzymes [22,23,24]. Four PMOs, i.e., PMOs 9, 8, 7, 4 were assessed in seven combinations and tested for their splice switching activity in SMA fibroblasts. Treatment of selected combination of PMO strongly increased expression of the hnRNPA1b isoform, but also induced SMN2 exon 5 skipping without substantial induction in SMN protein expression. We postulated that this effect was caused by disproportionate ratio of hnRNPA1/hnRNPA1b (50–60% hnRNPA1b) and SMN2 exon 5 skipping. This result prompted us to investigate whether induction of hnRNPA1b may have a global effect on gene expression in SMA fibroblasts. The transcriptome of cells transfected with PMOs9 + 8 + 7 combination was chosen for comparison to that of untreated SMA fibroblasts. Differential Expression analysis indicated decreased expression of survivin (BIRC5), Aurora kinase B (AURKB) and marker of proliferation ki-67 (MKI67) in PMOs9 + 8 + 7 cocktail treated SMA fibroblast. In keeping with this result, cardiomyocytes from the severe mouse model of SMA underwent cell cycle arrest and these three genes were also upregulated in response to increased SMN expression [25]. PMOs9 + 8 + 7 may decrease BIRC5, AURKB and MKI67 through modulating SMN expression through SMN exon 5 skipping. Gene ontology (GO) and gene set enrichment analysis (GSEA) on differentially expressed gene list were enriched for pathways involving cell cycle control, mitosis, nuclear division and chromosome segregation. The top twenty pathways identified by GSEA were statistically significant and nine out of those twenty pathways include the key master regulators, i.e., FOXM1 or MYBL2. RT-PCR demonstrated that treatment with the PMO cocktail reduced MYBL2 and FOXM1b expression in a type I SMA patient cells. MYBL2 and FOXM1, form complex with MuvB, and facilitate cell cycle progression through transcriptional activation of the cell cycle gene homology region (CHR) motif [26]. Cell cycle arrest and senescence is reported to be induced in MYBL2 or FOXM1 deficient cells, [27,28] while overexpression of FOXM1b contributes to the pathogenesis of many cancers i.e., colorectal, [29] nasopharyngeal [30] and renal cancer [31]. hnRNPA1 was reported to antagonise senescence in human lung fibroblasts, [32] and our PMO cocktail mediated switching between two hnRNPA1 isoforms, thereby reducing the expression of key regulators involved in the cell senescence process.

4. Materials and Methods

4.1. Culture Conditions of SMA Cell Line

A fibroblast cell line derived from a type I SMA patient (GM03813) was purchased from the Coriell Cell Repository (Camden, NJ, USA). The fibroblasts were maintained in Dulbecco’s modified Eagle medium (DMEM) with high glucose (Invitrogen, Waltham, MA, USA) containing penicillin G (100 units/mL), streptomycin (100 µg/mL) and 10% (v/v) foetal bovine serum (FBS). Fibroblasts were trypsinised and passaged two to three times per week using 0.25% (w/v) trypsin-EDTA, and were seeded into either a 24-well plate at a density of 1.5 × 104 cells or a 25 cm2 cultured flask at a density of 0.2 × 106 cells for AO transfection.

4.2. Antisense Oligonucleotides Design and Transfection into SMA Fibroblasts In Vitro

Four Intronic splice silencers (ISS) identified by Blanchette et al. [16] were chosen for AO targeting. Ten 2′ O-methyl modified AOs on a phosphorothioate backbone were designed to anneal to these four ISSs, as well as an exonic sequence enhancer (ESE) that was predicted by ESE finder 3.0, [33,34], as shown in Figure 1. The nomenclature of the AOs contains the gene name followed by the number of the targeted exon and the coordinates of the AO annealing site relative to either acceptor (A) or donor splice site (D). The AO HBS1LE6A(+114+139) designed to target HBS1L was used as a negative/sham control, and had been shown to be inactive at modulating SMN splicing. The AO hSMNE7D(-10-29), designed to target a strong intronic silencer element ISS-N1, is a potent FL-SMN transcript inducer, and was used as a positive control to promote SMN expression. The AO coordinates used in this study are listed in the supplementary Table S1. Each AO was individually transfected into SMA fibroblasts to determine its effects on the splicing pattern of hnRNPA1. Lipofectin reagent (LifeTechnology, Waltham, MA, USA) was used for AO transfection at a ratio of 2:1 Lipofectin to AO, as previously reported [35]. Briefly, lipofectin was mixed with DMEM to a final volume of 100 µL. After ten minutes, the lipofectin mixture was combined with 100 µL of the AO diluted in DMEM. The mixture was then made up to a final volume of 1 mL using antibiotic-free DMEM. At the time of transfection, 500 µL aliquots of this mixture were transferred to each well of a 24-well plate containing SMA fibroblasts and incubated for 48 h.

4.3. Transfection of Phosphorodiamidate Morpholino Oligomers in SMA Fibroblasts

Six PMOs (PMOs 4, 7, 8, 9, anti-ISS-N1 and a negative control from Gene Tools) were purchased from Gene Tools LLC (Philomath, OR, USA). PMOs were transfected into the cell by nucleofection using the P2 primary kit (Lonza, Basel, Switzerland). Two microliters of stock PMOs were added into a cuvette, and 250,000 of SMA fibroblasts were resuspended into 18 µL of prewarmed transfection solution, and then added to the cuvette. The CA-137 protocol on the Amaxa 4D machine (Lonza, Basel, Switzerland) was used for nucleofection. Nucleofected fibroblasts were incubated at 37 °C for 3 days prior to harvesting for further analysis.

4.4. RNA Extraction and RT-PCR

Total RNA was harvested from transfected Type I SMA fibroblasts 48 h after treatment with AOs using Trizol reagent, according to the manufacturer’s protocol (Invitrogen, Waltham, MA, USA). The RNA solution was then stored at −20 °C until use. Primers used to assess expression of hnRNPA1, SMN and other transcripts are listed in the supplementary Table S2 [36,37,38,39,40,41,42,43]. Approximately 200 ng of total RNA from each extraction was used as templates for single-step RT-PCR (Invitrogen, Waltham, MA, USA). Amplicons from the hnRNPA1 and SMN RT-PCRs were fractionated on 2% agarose gels and images were captured after staining with ethidium bromide. Product intensity was measured on an ImageQuant LAS 4010 (GE Healthcare Life Sciences, Chicago, IL, USA) and densitometric analysis was determined using image J software [44]. Percentage of SMN2 exon 7 inclusion is calculated by FL-SMN signal divided by the summation of density values of FL-SMN and SMNΔ7. The percentage of hnRNPA1 inclusion was determined as the density of exon 7b inclusion divided by the density values of the hnRNPA1 band and the hnRNPA1b band.

4.5. RNA Extraction and Transcriptomic Analysis

Total RNA was extracted from type I SMA fibroblasts transfected with PMOs 9 + 8 + 7 or untreated SMA fibroblasts using the RNeasy Mini Kit (Qiagen, Hilden, Germany) according to manufacturer’s protocol. RNA Sequencing (RNA-Seq) was performed by Illumina using 150 bp pair-end with 6 GB of output on the Hiseq. FastQC was used to filter out poor quality reads from FASTQ files [45]. Filtered reads were mapped against human reference genome (GRCh38) using Hisat2 v2.1.0 [46]. Picard Tools was used to remove duplicates [47]. HTseq was employed to generate input for differential expression analysis [48]. DESeq2 v-1.24.0 was used for differential gene expression analysis [49]. gProfiler2 was used for gene ontology analysis, [49,50], while GSEA-P was used for gene set enrichment analysis [51].

4.6. Immunoblotting

Three days after transfection, treated SMA fibroblasts were lysed in TRIS glycine buffer at pH 6.8, 10% glycerol with 4% SDS and sonicated with six pulses of 30% amplitude for 30 s on ice using Vibra-cellTM (Sonics and Materials INC., Newtown, CT, USA). Protein concentrations in extracts were measured using a BCA kit (Piercenet). Extracts were centrifuged for 5 min at 10,000× g at 4 °C. Sodium dodecyl sulphate (SDS)-PAGE electrophoresis was performed using 10% Tris-glycine gels and run at 100 volts for one hour at 4 °C. After electrophoresis, the gel was removed from the glass plates and the proteins were transferred to FluoroTrans® (Pall, Victoria, Australia) PVDF transfer membranes at 100 volts for 2 h at 4 °C using Omnipage sub blot system transfer system (Cleaver Scientific, Rugby, UK). Two primary antibodies were used: MANSMA1 (Santa Cruz Biotechnology, Dallas, TX, USA), a mouse anti-SMN monoclonal antibody, was applied at a 1:100 dilution for SMN protein detection, whereas the 9H10 antibody (Abcam, Cambridge, UK), a mouse anti-hnRNPA1 monoclonal antibody, was applied at 1:1000 dilution for hnRNPA1 protein detection. A rabbit β-tubulin monoclonal antibody (Abcam, Cambridge, UK) was used as a reference protein loading control. Secondary horseradish peroxidase (HRP)-conjugated antimouse or rabbit IgG antibodies were used at the dilution of 1 in 10,000. Images of the membrane were obtained using ImageQuant LAS 4010. Band intensity was analysed using the ImageJ software [43].

4.7. Statistical Analysis

Data are expressed as means ± standard deviation based on three separate experiments and the student T-test was used for statistical analysis.

5. Conclusions

In summary, our data confirm the role of hnRNPA1 and its isoforms in FL-SMN transcript expression, and that induction of the hnRNPA1b isoform is readily achievable using an optimised combination of AO sequences. Different levels of hnRNPA1b protein induction regulated exon selection during SMN2 pre-mRNA splicing. Induction of hnRNPA1 using PMO cocktail in SMA fibroblasts led to reduced expression of the master regulator of cell cycle progression. However, the effective PMO concentration range that can rescue expression of SMN protein appears to be narrow, and optimisation of the PMO cocktail, i.e., different dosing and the chosen AO component in the cocktail, should be evaluated together to refine the appropriate dosage. Further studies that use PMO treatment in animal models are required to determine effective dosage windows. The induction of the hnRNPA1b isoform in SMA mice would be a crucial avenue of study to provide a good indicator of the effectiveness of AO cocktails in the treatment of SMA.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijms23073937/s1.

Author Contributions

Conceptualisation and Funding acquisition, C.M. and S.D.W.; Methodology, J.T. and C.M.; Data Analysis, S.R., C.M., W.M., J.K.A. and V.C.; Original draft preparation, C.M.; Writing, T.N.P., S.F., S.D.W. and C.M.; Review and editing, W.M., J.K.A., V.C., B.S., T.N.P., S.F., S.D.W. and C.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by Thailand Research Fund (Grant number: TRG5880047) and Research Development fund, Faculty of Medicine, Siriraj Hospital, Mahidol University (Grant number: R0158333001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

RNAseq Data presented in this study are submitted to Sequence Read Archive database (SRA) at Bioproject: PRJNA821645.

Acknowledgments

C.M. is supported by Chalermprakiat fund, Faculty of Medicine, Siriraj Hospital. We thanks to Wipada Promkhuntod for the supporting role as laboratory manager.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lefebvre, S.; Burlet, P.; Liu, Q.; Bertrandy, S.; Clermont, O.; Munnich, A.; Dreyfuss, G.; Melki, J. Correlation between severity and SMN protein level in spinal muscular atrophy. Nat. Genet. 1997, 3, 265–269. [Google Scholar] [CrossRef] [PubMed]
  2. Prior, T.W.; Krainer, A.R.; Hua, Y.; Swoboda, K.J.; Snyder, P.C.; Bridgeman, S.J.; Burghes, A.H.; Kissel, J.T. A positive modifier of spinal muscular atrophy in the SMN2 gene. Am. J. Hum. Genet. 2009, 85, 408–413. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Singh, N.K.; Singh, N.N.; Androphy, E.J.; Singh, R.N. Splicing of a critical exon of human Survival Motor Neuron is regulated by a unique silencer element located in the last intron. Mol. Cell. Biol. 2006, 26, 1333–1346. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Singh, N.N.; Androphy, E.J.; Singh, R.N. In vivo selection reveals combinatorial controls that define a critical exon in the spinal muscular atrophy genes. RNA 2004, 10, 1291–1305. [Google Scholar] [CrossRef] [Green Version]
  5. Kashima, T.; Manley, J.L. A negative element in SMN2 exon 7 inhibits splicing in spinal muscular atrophy. Nat. Genet. 2003, 4, 460–463. [Google Scholar] [CrossRef]
  6. Kashima, T.; Rao, N.; David, C.J.; Manley, J.L. hnRNP A1 functions with specificity in repression of SMN2 exon 7 splicing. Hum. Mol. Genet. 2007, 16, 3149–3159. [Google Scholar] [CrossRef] [Green Version]
  7. Singh, N.N.; Hollinger, K.; Bhattacharya, D.; Singh, R.N. An antisense microwalk reveals critical role of an intronic position linked to a unique long-distance interaction in pre-mRNA splicing. RNA 2010, 16, 1167–1181. [Google Scholar] [CrossRef] [Green Version]
  8. Porensky, P.N.; Mitrpant, C.; McGovern, V.L.; Bevan, A.K.; Foust, K.D.; Kaspar, B.K.; Wilton, S.D.; Burghes, A.H. A single administration of morpholino antisense oligomer rescues spinal muscular atrophy in mouse. Hum. Mol. Genet. 2012, 21, 1625–1638. [Google Scholar] [CrossRef]
  9. Zhou, H.; Janghra, N.; Mitrpant, C.; Dickinson, R.L.; Anthony, K.; Price, L.; Eperon, I.C.; Wilton, S.D.; Morgan, J. A novel morpholino oligomer targeting ISS-N1 improves rescue of severe spinal muscular atrophy transgenic mice. Hum. Gene Ther. 2013, 24, 331–342. [Google Scholar] [CrossRef] [Green Version]
  10. Mitrpant, C.; Porensky, P.; Zhou, H.; Price, L.; Muntoni, F.; Fletcher, S.; Wilton, S.D.; Burghes, A.H. Improved antisense oligonucleotide design to suppress aberrant SMN2 gene transcript processing: Towards a treatment for spinal muscular atrophy. PLoS ONE 2013, 8, e62114. [Google Scholar] [CrossRef] [Green Version]
  11. Crooke, S.T.; Baker, B.F.; Pham, N.C.; Hughes, S.G.; Kwoh, T.J.; Cai, D.; Tsimikas, S.; Geary, R.S.; Bhanot, S. The effects of 2′-O-Methoxyethyl Oligonucleotides on renal function in humans. Nucleic Acid Ther. 2018, 28, 10–22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Crooke, S.T.; Baker, B.F.; Witztum, J.L.; Kwoh, T.J.; Pham, N.C.; Salgado, N.; McEvoy, B.W.; Cheng, W.; Hughes, S.G.; Bhanot, S.; et al. The effects of 2’-O-Methoxyethyl containing antisense oligonucleotides on platelets in human clinical trials. Nucleic Acid Ther. 2017, 27, 121–129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Frazier, K.S. Antisense oligonucleotide therapies: The promise and the challenges from a toxicologic pathologist’s perspective. Toxicol. Pathol. 2015, 43, 78–89. [Google Scholar] [CrossRef] [PubMed]
  14. Mayeda, A.; Munroe, S.H.; Caceres, J.F.; Krainer, A.R. Function of conserved domains of hnRNP A1 and other hnRNP A/B proteins. EMBO J. 1994, 13, 5483–5495. [Google Scholar] [CrossRef]
  15. Chabot, B.; Blanchette, M.; Lapierre, I.; La Branche, H. An intron element modulating 5′ splice site selection in the hnRNP A1 pre-mRNA interacts with hnRNP A1. Mol. Cell. Biol. 1997, 17, 1776–1786. [Google Scholar] [CrossRef] [Green Version]
  16. Blanchette, M.; Chabot, B. Modulation of exon skipping by high-affinity hnRNP A1-binding sites and by intron elements that repress splice site utilization. EMBO J. 1999, 18, 1939–1952. [Google Scholar] [CrossRef] [Green Version]
  17. Bernal, S.; Also-Rallo, E.; Martinez-Hernandez, R.; Alias, L.; Rodriguez-Alvarez, F.J.; Millan, J.M.; Hernandez-Chico, C.; Baiget, M.; Tizzano, E.F. Plastin 3 expression in discordant spinal muscular atrophy (SMA) siblings. Neuromuscul. Disord. 2011, 21, 413–419. [Google Scholar] [CrossRef]
  18. Foust, K.D.; Nurre, E.; Montgomery, C.L.; Hernandez, A.; Chan, C.M.; Kaspar, B.K. Intravascular AAV9 preferentially targets neonatal neurons and adult astrocytes. Nat. Biotechnol. 2009, 27, 59–65. [Google Scholar] [CrossRef] [Green Version]
  19. Foust, K.D.; Wang, X.; McGovern, V.L.; Braun, L.; Bevan, A.K.; Haidet, A.M.; Le, T.T.; Morales, P.R.; Rich, M.M.; Burghes, A.H.; et al. Rescue of the spinal muscular atrophy phenotype in a mouse model by early postnatal delivery of SMN. Nat. Biotechnol. 2010, 28, 271–274. [Google Scholar] [CrossRef] [Green Version]
  20. Cartegni, L.; Krainer, A.R. Disruption of an SF2/ASF-dependent exonic splicing enhancer in SMN2 causes spinal muscular atrophy in the absence of SMN1. Nat. Genet. 2002, 30, 377–384. [Google Scholar] [CrossRef]
  21. Blanchette, M.; Chabot, B. A highly stable duplex structure sequesters the 5’ splice site region of hnRNP A1 alternative exon 7B. RNA 1997, 3, 405–419. [Google Scholar] [PubMed]
  22. Cirak, S.; Arechavala-Gomeza, V.; Guglieri, M.; Feng, L.; Torelli, S.; Anthony, K.; Abbs, S.; Garralda, M.E.; Bourke, J.; Wells, D.J.; et al. Exon skipping and dystrophin restoration in patients with Duchenne muscular dystrophy after systemic phosphorodiamidate morpholino oligomer treatment: An open-label, phase 2, dose-escalation study. Lancet 2011, 378, 595–605. [Google Scholar] [CrossRef] [Green Version]
  23. Kinali, M.; Arechavala-Gomeza, V.; Feng, L.; Cirak, S.; Hunt, D.; Adkin, C.; Guglieri, M.; Ashton, E.; Abbs, S.; Nihoyannopoulos, P.; et al. Local restoration of dystrophin expression with the morpholino oligomer AVI-4658 in Duchenne muscular dystrophy: A single-blind, placebo-controlled, dose-escalation, proof-of-concept study. Lancet Neurol. 2009, 8, 918–928. [Google Scholar] [CrossRef] [Green Version]
  24. Moulton, J.D. Using Morpholinos to Control Gene Expression. Curr. Protoc. Nucleic Acid Chem. 2017, 68, 4–30. [Google Scholar] [CrossRef]
  25. Sheng, L.; Wan, B.; Feng, P.; Sun, J.; Rigo, F.; Bennett, C.F.; Akerman, M.; Krainer, A.R.; Hua, Y. Downregulation of Survivin contributes to cell-cycle arrest during postnatal cardiac development in a severe spinal muscular atrophy mouse model. Hum. Mol. Genet. 2018, 27, 486–498. [Google Scholar] [CrossRef]
  26. Engeland, K. Cell cycle arrest through indirect transcriptional repression by p53: I have a DREAM. Cell Death Differ. 2018, 25, 114–132. [Google Scholar] [CrossRef] [Green Version]
  27. Huang, Y.; Wu, J.; Li, R.; Wang, P.; Han, L.; Zhang, Z.; Tong, T. B-MYB delays cell aging by repressing p16 (INK4alpha) transcription. Cell Mol. Life Sci. 2011, 68, 893–901. [Google Scholar] [CrossRef]
  28. Smirnov, A.; Panatta, E.; Lena, A.; Castiglia, D.; Di Daniele, N.; Melino, G.; Candi, E. FOXM1 regulates proliferation, senescence and oxidative stress in keratinocytes and cancer cells. Aging 2016, 8, 1384–1397. [Google Scholar] [CrossRef] [Green Version]
  29. Yang, K.; Jiang, B.; Lu, Y.; Shu, Q.; Zhai, P.; Zhi, Q.; Li, Q. FOXM1 promotes the growth and metastasis of colorectal cancer via activation of beta-catenin signaling pathway. Cancer Manag. Res. 2019, 11, 3779–3790. [Google Scholar] [CrossRef] [Green Version]
  30. Luo, W.; Gao, F.; Li, S.; Liu, L. FoxM1 promotes cell proliferation, invasion, and stem cell properties in Nasopharyngeal Carcinoma. Front. Oncol. 2018, 8, 483. [Google Scholar] [CrossRef]
  31. Xue, Y.J.; Xiao, R.H.; Long, D.Z.; Zou, X.F.; Wang, X.N.; Zhang, G.X.; Yuan, Y.H.; Wu, G.Q.; Yang, J.; Wu, Y.T.; et al. Overexpression of FoxM1 is associated with tumor progression in patients with clear cell renal cell carcinoma. J. Transl. Med. 2012, 10, 200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Wang, H.; Han, L.; Zhao, G.; Shen, H.; Wang, P.; Sun, Z.; Xu, C.; Su, Y.; Li, G.; Tong, T.; et al. hnRNP A1 antagonizes cellular senescence and senescence-associated secretory phenotype via regulation of SIRT1 mRNA stability. Aging Cell 2016, 15, 1063–1073. [Google Scholar] [CrossRef] [PubMed]
  33. Cartegni, L.; Wang, J.; Zhu, Z.; Zhang, M.Q.; Krainer, A.R. ESEfinder: A web resource to identify exonic splicing enhancers. Nucleic Acids Res. 2003, 31, 3568–3571. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Smith, P.J.; Zhang, C.; Wang, J.; Chew, S.L.; Zhang, M.Q.; Krainer, A.R. An increased specificity score matrix for the prediction of SF2/ASF-specific exonic splicing enhancers. Hum. Mol. Genet. 2006, 15, 2490–2508. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Mitrpant, C.; Adams, A.M.; Meloni, P.L.; Muntoni, F.; Fletcher, S.; Wilton, S.D. Rational design of antisense oligomers to induce dystrophin exon skipping. Molecular therapy. J. Am. Soc. Gene Ther. 2009, 17, 1418–1426. [Google Scholar] [CrossRef] [PubMed]
  36. Bonomi, S.; di Matteo, A.; Buratti, E.; Cabianca, D.S.; Baralle, F.E.; Ghigna, C.; Biamonti, G. HnRNP A1 controls a splicing regulatory circuit promoting mesenchymal-to-epithelial transition. Nucleic Acids Res. 2013, 41, 8665–8679. [Google Scholar] [CrossRef] [Green Version]
  37. Flynn, L.L.; Mitrpant, C.; Pitout, I.L.; Fletcher, S.; Wilton, S.D. Antisense oligonucleotide-mediated terminal intron retention of the SMN2 transcript. Mol. Ther. Nucleic Acids 2018, 11, 91–102. [Google Scholar] [CrossRef] [Green Version]
  38. Ghigna, C.; Giordano, S.; Shen, H.; Benvenuto, F.; Castiglioni, F.; Comoglio, P.M.; Green, M.R.; Riva, S.; Biamonti, G. Cell motility is controlled by SF2/ASF through alternative splicing of the Ron protooncogene. Mol. Cell 2005, 20, 881–890. [Google Scholar] [CrossRef]
  39. Hatfield, J.T.; Rothnagel, J.A.; Smith, R. Characterization of the mouse hnRNP A2/B1/B0 gene and identification of processed pseudogenes. Gene 2002, 295, 33–42. [Google Scholar] [CrossRef]
  40. Karni, R.; de Stanchina, E.; Lowe, S.W.; Sinha, R.; Mu, D.; Krainer, A.R. The gene encoding the splicing factor SF2/ASF is a proto-oncogene. Nat. Struct. Mol. Biol. 2007, 14, 185–193. [Google Scholar] [CrossRef] [Green Version]
  41. Päth, G.; Bornstein, S.R.; Ehrhart-Bornstein, M.; Scherbaum, W.A. Interleukin-6 and the Interleukin-6 receptor in the human adrenal gland: Expression and effects on steroidogenesis1. J. Clin. Endocrinol. Metab. 1997, 82, 2343–2349. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Ren, F.; Wang, L.; Shen, X.; Xiao, X.; Liu, Z.; Wei, P.; Wang, Y.; Qi, P.; Shen, C.; Sheng, W.; et al. MYBL2 is an independent prognostic marker that has tumor-promoting functions in colorectal cancer. Am. J. Cancer Res. 2015, 5, 1542–1552. [Google Scholar] [PubMed]
  43. Zhang, X.; Zhang, L.; Du, Y.; Zheng, H.; Zhang, P.; Sun, Y.; Wang, Y.; Chen, J.; Ding, P.; Wang, N.; et al. A novel FOXM1 isoform, FOXM1D, promotes epithelial-mesenchymal transition and metastasis through ROCKs activation in colorectal cancer. Oncogene 2017, 36, 807–819. [Google Scholar] [CrossRef] [Green Version]
  44. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 2012, 9, 671–675. [Google Scholar] [CrossRef] [PubMed]
  45. Andrew, S. FastQC: A Quality Control Tool for High Throughput Sequence Data. 2010. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc (accessed on 1 February 2022).
  46. Kim, D.; Paggi, J.M.; Park, C.; Bennett, C.; Salzberg, S.L. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 2019, 37, 907–915. [Google Scholar] [CrossRef]
  47. Broad Institute. Picard Tools. Available online: https://github.com/broadinstitute/picard (accessed on 15 January 2021).
  48. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef]
  49. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  50. Reimand, J.; Kull, M.; Peterson, H.; Hansen, J.; Vilo, J. g:Profiler—A web-based toolset for functional profiling of gene lists from large-scale experiments. Nucleic Acids Res. 2007, 3, W193–W200. [Google Scholar] [CrossRef]
  51. Subramanian, A.; Kuehn, H.; Gould, J.; Tamayo, P.; Mesirov, J.P. GSEA-P: A desktop application for Gene Set Enrichment Analysis. Bioinformatics 2007, 23, 3251–3253. [Google Scholar] [CrossRef]
Ijms 23 03937 g001 550
Figure 1. Schematic diagram displaying the targeting splice motifs and AO coordinates of ten AOs induced hnRNPA1 exon 7b skipping (AOs1 and 2) and inclusion (AOs3–10).
Figure 1. Schematic diagram displaying the targeting splice motifs and AO coordinates of ten AOs induced hnRNPA1 exon 7b skipping (AOs1 and 2) and inclusion (AOs3–10).
Ijms 23 03937 g001
Ijms 23 03937 g002 550
Figure 2. RT-PCR analysis demonstrated splicing patterns of hnRNPA1, (A) and percentage of exon 7b inclusion in SMA fibroblasts transfected with different 2′O-methyl modified AOs (B).
Figure 2. RT-PCR analysis demonstrated splicing patterns of hnRNPA1, (A) and percentage of exon 7b inclusion in SMA fibroblasts transfected with different 2′O-methyl modified AOs (B).
Ijms 23 03937 g002
Ijms 23 03937 g003 550
Figure 3. Gel images and densitometric analysis of RT-PCR of hnRNPA1 (A,B) and SMN (C,D) splicing isoforms in SMA fibroblasts transfected by seven combinations of AO cocktails, as well as densitometric analysis of the percentage of FL-SMN. The indicated concentrations were a dose of combined AOs in the cocktail. The * symbol indicates the p-value of the average percentage of exon inclusion compared to untreated samples (viz ** = p < 0.01, *** = p < 0.001 and **** = p < 0.0001).
Figure 3. Gel images and densitometric analysis of RT-PCR of hnRNPA1 (A,B) and SMN (C,D) splicing isoforms in SMA fibroblasts transfected by seven combinations of AO cocktails, as well as densitometric analysis of the percentage of FL-SMN. The indicated concentrations were a dose of combined AOs in the cocktail. The * symbol indicates the p-value of the average percentage of exon inclusion compared to untreated samples (viz ** = p < 0.01, *** = p < 0.001 and **** = p < 0.0001).
Ijms 23 03937 g003
Ijms 23 03937 g004 550
Figure 4. Western blot analysis indicates the hnRNPA1 (34 kDa) and the induced isoform of hnRNPA1 proteins at approximately 38 kDa (A), SMN protein at approximately 38 kDa (C) and β-tubulin protein at approximately 50 kDa (C) expressed in SMA fibroblasts that were transfected with selected seven different AO cocktails, as well as densitometric analysis of percentage of hnRNPA1b (B), as compared to SMN protein expression normalised against β-tubulin (D). The * symbol indicates the p-value of the average relative protein expression compared to untreated samples (viz * = p < 0.05, ** = p < 0.01 and *** = p < 0.001).
Figure 4. Western blot analysis indicates the hnRNPA1 (34 kDa) and the induced isoform of hnRNPA1 proteins at approximately 38 kDa (A), SMN protein at approximately 38 kDa (C) and β-tubulin protein at approximately 50 kDa (C) expressed in SMA fibroblasts that were transfected with selected seven different AO cocktails, as well as densitometric analysis of percentage of hnRNPA1b (B), as compared to SMN protein expression normalised against β-tubulin (D). The * symbol indicates the p-value of the average relative protein expression compared to untreated samples (viz * = p < 0.05, ** = p < 0.01 and *** = p < 0.001).
Ijms 23 03937 g004
Ijms 23 03937 g005 550
Figure 5. Gel images of RT-PCR shows hnRNPA1 (827 bp) and hnRNPA1b (983 bp) transcripts (A). Panel (B),(C) shows gel images of different isoforms of SMN2 transcripts missing exon 5 (Δ5), missing exon 7 (Δ7) or missing both exons 5 and 7 (Δ5 + 7) with densitometric analysis of SMN2 transcript. Western blot analysis of hnRNPA1, FL-SMN and β-tubulin expression in SMA fibroblasts transfected with selected seven different PMO cocktails as well as densitometric analysis of percentage of hnRNPA1b compared to SMN protein expression normalised against β-tubulin (D–G). Numbers above each bar graph indicates the average values of percentage of hnRNPA1b isoform (E) and the average fold change of increased SMN protein (G). ** = p < 0.01, *** = p < 0.001.
Figure 5. Gel images of RT-PCR shows hnRNPA1 (827 bp) and hnRNPA1b (983 bp) transcripts (A). Panel (B),(C) shows gel images of different isoforms of SMN2 transcripts missing exon 5 (Δ5), missing exon 7 (Δ7) or missing both exons 5 and 7 (Δ5 + 7) with densitometric analysis of SMN2 transcript. Western blot analysis of hnRNPA1, FL-SMN and β-tubulin expression in SMA fibroblasts transfected with selected seven different PMO cocktails as well as densitometric analysis of percentage of hnRNPA1b compared to SMN protein expression normalised against β-tubulin (D–G). Numbers above each bar graph indicates the average values of percentage of hnRNPA1b isoform (E) and the average fold change of increased SMN protein (G). ** = p < 0.01, *** = p < 0.001.
Ijms 23 03937 g005
Ijms 23 03937 g006 550
Figure 6. Differential expression comparing between PMOs9 + 7 + 8 treated SMA fibroblast and untreated SMA fibroblasts. Heatmap of the differentially expressed genes in three replicates of both conditions. (A) Volcano plot of differentially expressed gene list comparing between cells transfected with PMOs9 + 8 + 7 and untreated SMA fibroblasts. (B) Selected genes with absolute log2 fold change greater than 1.5 folds and p-value lesser than 0.05 was analysed for its gene ontology. Top twenty of biological processes that were upregulated (C) or downregulated (D) in PMO treatment group compared to untreated control. Gene set enrichment analysis indicated that gene set regulating cell cycles were enriched at the normalised enrichment score (NES) of −1.8477 and false discovery rate (FDR) q-value of 0.00 (E).
Figure 6. Differential expression comparing between PMOs9 + 7 + 8 treated SMA fibroblast and untreated SMA fibroblasts. Heatmap of the differentially expressed genes in three replicates of both conditions. (A) Volcano plot of differentially expressed gene list comparing between cells transfected with PMOs9 + 8 + 7 and untreated SMA fibroblasts. (B) Selected genes with absolute log2 fold change greater than 1.5 folds and p-value lesser than 0.05 was analysed for its gene ontology. Top twenty of biological processes that were upregulated (C) or downregulated (D) in PMO treatment group compared to untreated control. Gene set enrichment analysis indicated that gene set regulating cell cycles were enriched at the normalised enrichment score (NES) of −1.8477 and false discovery rate (FDR) q-value of 0.00 (E).
Ijms 23 03937 g006
Ijms 23 03937 g007 550
Figure 7. Isoform specific RT-PCR to detect four different isoforms of FOXM1 (A) with the graph displaying band densitometry in the panels (B–E) (viz FOXM1a, FOXM1b, FOXM1c and FOXM1d) showed reduced expression of FOXM1b in cultured SMA fibroblasts treated with the combination of PMOs 9 + 8 + 7, while RT-PCR analysis of MYBL2 (F) was also reduced, as indicated in its band densitometry (G).
Figure 7. Isoform specific RT-PCR to detect four different isoforms of FOXM1 (A) with the graph displaying band densitometry in the panels (B–E) (viz FOXM1a, FOXM1b, FOXM1c and FOXM1d) showed reduced expression of FOXM1b in cultured SMA fibroblasts treated with the combination of PMOs 9 + 8 + 7, while RT-PCR analysis of MYBL2 (F) was also reduced, as indicated in its band densitometry (G).
Ijms 23 03937 g007
Table
Table 1. List of AO combinations targeting the four splice silencing motifs located in introns adjacent to hnRNPA1 exon 7b.
Table 1. List of AO combinations targeting the four splice silencing motifs located in introns adjacent to hnRNPA1 exon 7b.
AO CombinationSilencing Motifs Targeted
AOs9 + 8ISS(CE1A) + ISS(CE4)
AOs9 + 7ISS(CE1A) + ISS(CE6)
AOs9 + 4ISS(CE1A) + ISS(CE9)
AOs9 + 8 + 7ISS(CE1A) + ISS(CE4) + ISS(CE6)
AOs9 + 8 + 4 ISS(CE1A) + ISS(CE4) + ISS(CE9)
AOs9 + 7 + 4ISS(CE1A) + ISS(CE6) + ISS(CE9)
AOs9 + 8 + 7 + 4ISS(CE1A) + ISS(CE4) + ISS(CE6) + ISS(CE9)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Toosaranont, J.; Ruschadaariyachat, S.; Mujchariyakul, W.; Arora, J.K.; Charoensawan, V.; Suktitipat, B.; Palmer, T.N.; Fletcher, S.; Wilton, S.D.; Mitrpant, C. Antisense Oligonucleotide Induction of the hnRNPA1b Isoform Affects Pre-mRNA Splicing of SMN2 in SMA Type I Fibroblasts. Int. J. Mol. Sci. 2022, 23, 3937. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms23073937

AMA Style

Toosaranont J, Ruschadaariyachat S, Mujchariyakul W, Arora JK, Charoensawan V, Suktitipat B, Palmer TN, Fletcher S, Wilton SD, Mitrpant C. Antisense Oligonucleotide Induction of the hnRNPA1b Isoform Affects Pre-mRNA Splicing of SMN2 in SMA Type I Fibroblasts. International Journal of Molecular Sciences. 2022; 23(7):3937. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms23073937

Chicago/Turabian Style

Toosaranont, Jarichad, Sukanya Ruschadaariyachat, Warasinee Mujchariyakul, Jantarika Kumar Arora, Varodom Charoensawan, Bhoom Suktitipat, Thomas N. Palmer, Sue Fletcher, Steve D. Wilton, and Chalermchai Mitrpant. 2022. "Antisense Oligonucleotide Induction of the hnRNPA1b Isoform Affects Pre-mRNA Splicing of SMN2 in SMA Type I Fibroblasts" International Journal of Molecular Sciences 23, no. 7: 3937. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms23073937

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop