Mol. Cells 2014; 37(1): 81-87
Published online January 27, 2013
https://doi.org/10.14348/molcells.2014.2337
© The Korean Society for Molecular and Cellular Biology
Correspondence to : *Correspondence:
Chronic pressure-overload cardiac hypertrophy is associated with an increased risk of morbidity/mortality, largely due to maladaptive remodeling and dilatation that progresses to dilated cardiomyopathy. Alternative splicing is an important biological mechanism that generates proteomic complexity and diversity. The recent development of next-generation RNA sequencing has improved our understanding of the qualitative signatures associated with alternative splicing in various biological conditions. However, the role of alternative splicing in cardiac hypertrophy is yet unknown. The present study employed RNA-Seq and a bioinformatic approach to detect the RNA splicing regulatory elements involved in alternative splicing during pressure-overload cardiac hypertrophy. We found GC-rich exonic motifs that regulate intron retention in 5′ UTRs and AT-rich exonic motifs that are involved in exclusion of the AT-rich elements that cause mRNA instability in 3′ UTRs. We also identified motifs in the intronic regions involved in exon exclusion and inclusion, which predicted splicing factors that bind to these motifs. We found, through Western blotting, that the expression levels of three splicing factors, ESRP1, PTB and SF2/ASF, were significantly altered during cardiac hypertrophy. Collectively, the present results suggest that chronic pressure-overload hypertrophy is closely associated with distinct alternative splicing due to altered expression of splicing factors.
Keywords cardiac hypertrophy, ESRP1, PTB, RNA-Seq, SF2/ASF, splicing factor
Cardiac hypertrophy is an adaptive response to various stress conditions such as pressure-overload. Untreated pathological cardiac hypertrophy is associated with an increased risk of morbidity and mortality, largely due to maladaptive remodeling and dilatation that progresses to dilated cardiomyopathy and heart failure.
Alternative splicing is a key biological event that contributes to proteome diversity and complexity by generating distinct mRNAs in different circumstances, which can influence gene expression. Aberrant splicing can have serious pathological consequences (Kim et al., 2008; Zhou et al., 2012). Our understanding of transcriptomic signatures, including qualitative changes due to alternative splicing, have been vastly increased by the recent development of next-generation RNA sequencing (RNA-Seq) methods. RNA-Seq on various cancers has helped improve the strategies for diagnosis and treatment (Eswaran et al., 2013; Kalari et al., 2012; Ren et al., 2012). Although high-throughput analyses of alternative splicing in cardiac diseases have recently been performed (Ames et al., 2013; Park et al., 2011; Song et al., 2012), the detailed mechanisms of altered alternative splicing during cardiac hypertrophy have not yet been reported.
In the present study, we performed an in silico analysis to detect the binding motifs of splicing factors that regulate alternative splicing during pressure-overload cardiac hypertrophy using RNA-Seq datasets. Through this analysis, we found GC-rich exonic motif-mediated intron retention in 5′ untranslated regions (UTRs) and AT-rich exonic motif-mediated exclusion of AREs that conferred mRNA instability in 3′ UTRs. We also identified 20 intronic motifs with a strong positional bias in the introns surrounding variant exons, and experimentally confirmed the differential expression of three splicing factors corresponding to particular motifs. Collectively, our study provides the first report of the unique mechanisms of alternative splicing controlled by various regulatory elements, and identifies the splicing factors that are associated with cardiac hypertrophy.
Our previous RNA-Seq datasets for cardiac hypertrophy were utilized (Song et al., 2012). Exon variants were analyzed based on RefSeq structure and less expressed exons < 10 reads per kilobase per million reads (RPKM) were discarded. Fisher’s exact test (p < 0.1) and Bayesian error rate (Vencio et al., 2004) (e ≤ 0.1) were used to determine the exon variants. Non-regulated exons were decided by p > 0.5.
Given exon variants, 500 bp intronic sequences (R1-R4), defined as Fig. 2A, were collected. If the intron size was smaller than 500 bp, whole intron was accepted. Considering the previously reported approaches for examining the regulatory motifs by enriched hexanucleotides in a position-dependent manner (Erkelenz et al., 2013; Fairbrother et al., 2002; Hirose et al., 2006; Huh and Hynes, 1994; Kuroyanagi et al., 2013; Lim et al., 2011; Yamashita et al., 2012), we counted the occurrence of each hexanucleotide at the position [ith ∼ i+5th] of the introns. The enrichment of hexanucleotides at position i was analyzed by Fisher’s exact test (p < 0.05). In case of exon inclusion, two independent tests were performed for exon exclusion and non-regulated exons as follows. The enrichment of the hexanucleo-tides at position i of exon inclusion was compared with the number of the same motif at the same position i observed in exon exclusion, and vice versa for exon exclusion.
The enrichment of genes for KEGG pathways was analyzed through Fisher’s exact test (p < 0.05) as follows.
The 8 weeks old male (C57BL/6J) mice (body weight 28?33 g) were used. Cardiac hypertrophy was induced by transverse aortic constriction (TAC) operation under anesthesia with intraperitoneal injection of avertin, 2-2-2 tribromoethanol (Sigma, USA) dissolved in tert-amyl alcohol (Sigma, USA). The procedure of operation was followed as previously described (Cha et al., 2008). As a control group, sham operation (same procedure except for tying) was done. One week after the operation, mice were sacrificed, and hearts were removed and stored in deep freezer at ?80°C. Validation of cardiac hypertrophy model was performed by the measurement of heart weight and body weight ratio and quantification for mRNA expression of hypertrophy markers (ANP and BNP; Fig. 3A).
First-strand cDNA was synthesized from 2 μg of total RNA with Random hexamer using Omniscript? reverse transcription (Qiagen, USA) according to the manufacturer’s instruction. qRT-PCR assays were followed as previously described (Hong et al., 2008). Briefly, qRT-PCR assays were performed using TOPreal™ qPCR premix (Enzynomics, Korea) under the following two-step conditions: denaturation at 95°C for 15 s followed by annealing and extension at 60°C for 40 s, for a total of 40 cycles. The 18S transcript was used as an endogenous reference to assess the relative level of mRNA transcript. The following primers were used:
ANP, 5′-TCGTCTTGGCCTTTTGGC-3′ (sense)
and 5′-TCCAGGTGGTCTAGCAGGTTCT-3′ (antisense)
BNP, 5′-AGGGAGAACACGGCATCATT-3′ (sense)
and 5′-GACAGCACCTTCAGGAGAT-3′ (antisense)
18S, 5′-TTCTGGCGAACGGTCTAGACAAC-3′ (sense)
and 5′-CCACTGGTCTTGGTGTGCTGA-3′ (antisense)
Frozen mouse hearts were homogenized in liquid nitrogen using a mortar and pestle. Homogenized heart powder was lysed in 1 % SDS lysis buffer containing protease inhibitor cocktail, Na3VO4 and NaF. Proteins were estimated by using BCA assay kit (Thermo scientific) and 30 μg of proteins were resolved electrophoretically on 10% SDS polyacrylamide gels, transferred to PVDF membrane. The blots were probed with following antibodies and detected by chemiluminescent substrate (Thermo scientific) and Biomolecular imager, Image-Quant LAS 4000 mini (GE). The band intensities were quantified using NIH ImageJ software.
The following primary antibodies were used for Western blot analysis: ESRP1 (210-301-B89) was purchased from Rockland. PTB, SF2/ASF, SRp30c, and SRp40 antibodies were kind gifts from Dr. Hai-Hong Shen (Gwangju Institute of Science and Technology, Korea).
Our previous study showed the transcriptomic signatures of cardiac hypertrophy and identified hundreds of exon variants that were included or excluded (Song et al., 2012). In the present study, the putative motifs that recruit splicing factors and regulate alternative splicing during cardiac hypertrophy were investigated. For this study, we first collected 3,345 and 3,101 included exons and excluded exons, respectively, and then identified 4,125 that were not regulated (p > 0.5) during pressure-overload cardiac hypertrophy (See “Materials and Methods”). The distributions of alternative splicing types for the exons of interest (i.e., included exons, excluded exons and non-regulated exons) are shown in Supplementary Fig. 1. The results indicate that alternative promoter (p = 7.8 × 10?11), alternative termination (p = 3.7 × 10?22) and cassette exons (p = 0.049) are differentially distributed between exon exclusion and exon inclusion.
The exon variants were characterized in terms of their nucleotide composition, length, and the distribution of 5′- or 3′ UTRs. As shown in Supplementary Fig. 2A, the nucleotide composition was somewhat distinct across the exon variants; an even distribution of A, T, G, and C (ranging from 24% to 25%) was observed in the non-regulated exons, whereas a significantly biased distribution of GC (∼28%) and AT (∼26%) was observed in included exons and excluded exons, respectively (p = 2.2 × 10?16; X2-test). Interestingly, the distribution of 5′ UTRs and 3′ UTRs in the included exons and excluded exons exhibited an obvious dichotomous pattern between; exons containing 5′ UTRs were more likely to be included, whereas exons containing 3′ UTRs tended to be excluded during pressure-overload cardiac hypertrophy (Supplementary Fig. 2B). Longer exons were more likely to be excluded; however, no significant difference was observed in median length across all exons (Supplementary Fig. 2C).
The biased distribution of GC/AT between included and excluded exons and the dichotomous patterns of 5′ UTRs and 3′ UTRs suggest specific architectures (e.g., motifs) that recruit specific splicing factors for either exon inclusion or exon exclusion. To predict the putative motifs in the exon variants, we searched for motifs composed of hexanucleotides that were significantly enriched in an inclusion- or exclusion-specific manner. We found that, among all possible hexanucleotides (i.e., 46 = 4,096), 29 and 24 hexanucleotides were abundantly distributed in excluded exons and included exons, respectively (Supplementary Table 1). Consistent with the high AT-content pattern observed in the excluded exons (Supplementary Fig. 2A), motifs with a high proportion of AT were enriched in excluded exons compared to included exons or non-regulated exons (Fig. 1A). The enrichment of AT-rich motifs implies that AT-rich elements (AREs) contribute to the instability of mRNA in the 3′ UTRs (Barreau et al., 2005). We also found that ‘alternative terminal exons’, defined by the ‘knownAlt’ track of the UCSC genome browser, and 3′ UTRs were significantly enriched in the exons containing AT-rich motifs (p < 0.05), implying the removal of AREs in 3′ UTRs (Fig. 1C and Supplementary Table 1). In contrast, motifs with high GC content were observed in included exons (Fig. 1B and Supplementary Table 1). Consistent with the high proportion of 5′ UTRs observed in the included exons, 5′ UTRs participated in the inclusion of exons with anchoring GC-rich motifs (Fig. 1D). The exon inclusion capability of 5′ UTR content does not appear to be related to an alternative promoter, as no known alternative promoters were significantly enriched (Supplementary Table 1). Rather, it seems to be related to intron retention (e.g., p-value = 8.19 × 10?5?6.16 × 10?8 for GGCCGC, GGCGGC, and CGGCGG; Fisher’s exact test). In light of previous reports that exon borders with higher GC content tended to retain introns (Galante et al., 2004; Wong et al., 2013), this result suggests that the GC-rich motifs located in the 5′ UTRs might be strongly related to intron retention during cardiac hypertrophy.
In summary, we found that exons that contain a GC-rich motif in the 5′ UTRs were significantly included and that AT-rich exons with 3′ UTRs were highly excluded during pressure-overload cardiac hypertrophy. These findings also suggest possible roles for GC-rich and AT-rich motifs in intron retention and exclusion of ARE, respectively.
Similar to what was observed in the exonic motif analysis, the enriched intronic motifs associated with exon variants were also investigated. Due to the importance of proximity to exons, the distribution of all possible hexanucleotides within the 500 base pairs (bp) flanking the intronic sequences (i.e., R1?R4) were analyzed (Fig. 2A). We found 13 and 7 motifs with strong positional bias within the 70 bp adjacent to the introns of excluded exons and included exons, respectively (Figs. 2B and 2C, and Supplementary Table 2). Interestingly, R1 was the most important region, anchoring eight motifs (52%) of the excluded exons. For the included exons, the ∼20 bp upstream of the exons (i.e., R2) were likely to be critical for the binding of splicing factors.
For intronic motifs, we found a correlation between G-rich motifs and intron length. Introns containing G-rich motifs were much shorter than introns containing other motifs. The average intron size in the mouse ranges from 3.8 to 5 kb (Taft et al., 2007; Waterston et al., 2002), whereas the length of the introns upstream and downstream of the exon variants containing G-rich motifs, such as GGAGGG, were less than 1 kb, suggesting a role for these G-rich motifs in the exclusion of exons flanking short introns (Supplementary Fig. 3). Since ESRP1, SF2/ASF, SRp30c, and SRp40 were reported to bind to G-rich motifs (such as GGTGGG, GGAGAG, and GGAGGG) (Supplementary Table 2), we experimentally examined the expression level of splicing factors in Sham/TAC-operated heart sample and foundthe differential expression of SF2/ASF and ESRP1 (Fig. 3). G-rich intronic element at 5′ splice site was known to regulate alternative splicing of thyroid hormone receptor by SF2/ASF (Hastings et al., 2001) and we also found the significantly enriched G-rich motif at 5’ splice site along with significantly increased level of SF2/ASF during cardiac hypertrophy. G-rich motifs are also involved in the regulation of alternative splicing by ESRP (Dittmar et al., 2012; Warzecha et al., 2010). Of the G-rich motifs, GTGGG-containing motifs (i.e., GGTGGG and GTGGGG) were significantly enriched in multiple sites for R1 and R3 of 57 exons, which were highly involved in cellular movement (p = 7.62 × 10?5; Fisher’s exacttest). ESRP is known to bind to TGG motifs, which are located within these GTGGG motifs, and play an important role in cellular movements such as cell adhesion, cell polarity, and cell migration (Warzecha et al., 2010). Interestingly, the GTGGG motifs were found along with TGCCTG motifs in the adjacent 50 bp of R2 (p = 1.49 × 10?4), suggesting close interplay between these two motifs in exon exclusion. The co-occurrence of TGCCTG at the 3′ splicing acceptor of the exons (i.e., R2) and GTGGG at the 5′ splicing donor of excluded exons (i.e., R3) were observed in Shank3, Vps16, Rhbdl1, Arhgef2, Nbeal2, and Brd9 (data not shown). TGC/CTG repeats are known binding motifs for CELF proteins, which are critical splicing factors during cardiac development (Ladd et al., 2001). Besides the G-rich motifs, we also found the significantly increased PTB known to bind to TCTT motif, during cardiac hypertrophy (Das et al., 2007).
Taken together, we predicted that intronic motifs were significantly enriched at specific positions in introns, and experimentally confirmed the expression of the splicing factors that bind to these motifs. The role of ESRP1 is notable, inducing the exclusion of exons flanking a short intron in genes involved in cellular movement.
During pressure overload cardiac hypertrophy, the heart undergoes rapid enlargement, remodeling and irreversible maladaptation accompanied by dramatically increased energy consumption due to rapid cardiac growth and the synthesis of cellular building blocks. Therefore, systematic cellular and metabolic transitions are necessary to meet the substantial demands during cardiac remodeling. Even though the quantitative signatures of pathological cardiac hypertrophy have been well studied (Park et al., 2011; Song et al., 2012), the qualitative changes during cardiac remodeling have remained to be investigated. Here, we show the evidence for the association of alternative splicing with pressure-overload hypertrophy as follows: 1) GC-rich exonic motifs regulate intron retention in 5′ UTRs and AT-rich exonic motifs are involved in exclusion of the AT-rich elements in 3′ UTRs during cardiac hypertrophy (Figs. 1 and 4, and Supplementary Table 1); 2) 20 putative splicing factor binding motifs enriched in intronic regions were identified (Fig. 2 and Supplementary Table 2); 3) Splicing factors such as PTB and ESRP1 that are known to bind to the putative motifs were differentially expressed during cardiac hypertrophy (Fig. 3); 4) Collectively, the splicing regulatory elements are involved in post-transcriptional regulation via the 3′ or 5′ UTRs or by recruiting splicing factors during cardiac hypertrophy (Fig. 4).
As summarized in Fig. 4A, alternative splicing for both 5′ and 3′ UTRs may play an important role in regulating translation during cardiac remodeling. The observed retention of introns in 5′ UTRs is important, because the mechanism could be involved in post-transcriptional regulation by inserting regulatory elements such as internal ribosome entry site (IRES) and G-quadruplex into mRNA. According to the previous reports, the intron retention led to the acquisition of IRES in Zeb2 and to the insertion of G-quadruplex in p53 (Beltran et al., 2008; Marcel et al., 2011). In the case of alternative splicing for 3′ UTR, we observed the removal of AREs that make the mRNA unstable. These findings suggest that post-transcriptional regulation by alternative UTRs might increase mRNAstability, thereby maximizing protein synthesis without significantly increasing transcription. Based on the rapid increase in heart volume during cardiac hypertrophy, this mechanism is an excellent way to adapt against urgent cardiac overload.
As shown in Fig. 4B, we found the altered expression of splicing factors that bind to the intronic motifs. Among them, PTB and ESRP1 are already known to be involved in alternative splicing of smooth muscle and metastasis in cancer, respectively (Gooding et al., 2013; Reinke et al., 2012; Yae et al., 2012), but their roles in cardiac hypertrophy are not known. They likely bind to the G-rich regions and are widely distributed nearby exons excluded (Supplementary Table 2). Warzecha et al. (2010) reported that UGG is a core motif of ESRP1 and the motif located in upstream intron possibly promotes exon exclusion, whereas the motifs in downstream of the regulated exon enhance exon inclusion. Their following study on the ESRP1 binding motifs using SELEX-Seq predicted the enriched GU-rich motifs in the flanking upstream (R2) of the excluded exons. For the included exons, they found the motifs significantly enriched in the distal upstream (R1) and flanking downstream (R3). During cardiac hypertrophy, the decreased level of ESRP1 might result in the weak binding at R1 and R3 as reported recently (Dittmar et al., 2012).
Intronic GU-rich repeats at 5′ splice site are also known to be the binding region for SF2/ASF and hnRNP H (Arrisi-Mercado et al., 2004; Hastings et al., 2001). Similarly, we also found the significantly increased level of SF2/ASF during cardiac hypertrophy. However, the detailed mechanistic relationship between the differential expression of the splicing factors (e.g., ESRP1 and SF2/ASF) and the G-rich motifs remains to be elucidated, as the splicing factors usually play roles as both enhancer and repressor depending on the position, concentration and the location of the branchpoint (Lim et al., 2011; Taggart et al., 2012; Zuo and Manley, 1993).
Consistent with the known role of ESRP1 in the alternative splicing of genes involved in cellular movement such as cell adhesion, polarity and migration (Warzecha et al., 2010), we also identified 57 genes with the ESRP1 binding motif that are involved in cellular movement. The significantly decreased expression of ESRP1 (Fig. 3) also suggests the importance of ESRP1 in alternative splicing during cardiac hypertrophy. We also identified a significant co-occurrence of the binding motifs for ESRP1 and CELF. However, the expression of the gene encoding CELF did not change during cardiac hypertrophy (data not shown). Although there was no change in the expression of CELF genes (e.g., Cugbp1, Cugbp2, Brunol4, and Brunol5), CELF expression can possibly be regulated post-transcriptionally through microRNAs (Kalsotra et al., 2010).
Here, we report the splicing regulatory elements and the putative splicing factors that may play important roles for pathological cardiac hypertrophy. The molecular mechanisms for alternative splicing regulated by the splicing regulatory elements and the functional roles of the individual splicing factors identified from this study will await the future molecular studies.
SUPPLEMENTARY FILE
Mol. Cells 2014; 37(1): 81-87
Published online January 31, 2014 https://doi.org/10.14348/molcells.2014.2337
Copyright © The Korean Society for Molecular and Cellular Biology.
Taeyong Kim, Jin Ock Kim, Jae Gyun Oh1, Seong-Eui Hong*, and Do Han Kim*
School of Life Sciences and Systems Biology Research Center, Gwangju Institute of Science and Technology, Gwangju 500-712,Korea, 1Cardiovascular Research Center, Icahn School of Medicine at Mount Sinai, New York, 10029,USA
Correspondence to:*Correspondence:
Chronic pressure-overload cardiac hypertrophy is associated with an increased risk of morbidity/mortality, largely due to maladaptive remodeling and dilatation that progresses to dilated cardiomyopathy. Alternative splicing is an important biological mechanism that generates proteomic complexity and diversity. The recent development of next-generation RNA sequencing has improved our understanding of the qualitative signatures associated with alternative splicing in various biological conditions. However, the role of alternative splicing in cardiac hypertrophy is yet unknown. The present study employed RNA-Seq and a bioinformatic approach to detect the RNA splicing regulatory elements involved in alternative splicing during pressure-overload cardiac hypertrophy. We found GC-rich exonic motifs that regulate intron retention in 5′ UTRs and AT-rich exonic motifs that are involved in exclusion of the AT-rich elements that cause mRNA instability in 3′ UTRs. We also identified motifs in the intronic regions involved in exon exclusion and inclusion, which predicted splicing factors that bind to these motifs. We found, through Western blotting, that the expression levels of three splicing factors, ESRP1, PTB and SF2/ASF, were significantly altered during cardiac hypertrophy. Collectively, the present results suggest that chronic pressure-overload hypertrophy is closely associated with distinct alternative splicing due to altered expression of splicing factors.
Keywords: cardiac hypertrophy, ESRP1, PTB, RNA-Seq, SF2/ASF, splicing factor
Cardiac hypertrophy is an adaptive response to various stress conditions such as pressure-overload. Untreated pathological cardiac hypertrophy is associated with an increased risk of morbidity and mortality, largely due to maladaptive remodeling and dilatation that progresses to dilated cardiomyopathy and heart failure.
Alternative splicing is a key biological event that contributes to proteome diversity and complexity by generating distinct mRNAs in different circumstances, which can influence gene expression. Aberrant splicing can have serious pathological consequences (Kim et al., 2008; Zhou et al., 2012). Our understanding of transcriptomic signatures, including qualitative changes due to alternative splicing, have been vastly increased by the recent development of next-generation RNA sequencing (RNA-Seq) methods. RNA-Seq on various cancers has helped improve the strategies for diagnosis and treatment (Eswaran et al., 2013; Kalari et al., 2012; Ren et al., 2012). Although high-throughput analyses of alternative splicing in cardiac diseases have recently been performed (Ames et al., 2013; Park et al., 2011; Song et al., 2012), the detailed mechanisms of altered alternative splicing during cardiac hypertrophy have not yet been reported.
In the present study, we performed an in silico analysis to detect the binding motifs of splicing factors that regulate alternative splicing during pressure-overload cardiac hypertrophy using RNA-Seq datasets. Through this analysis, we found GC-rich exonic motif-mediated intron retention in 5′ untranslated regions (UTRs) and AT-rich exonic motif-mediated exclusion of AREs that conferred mRNA instability in 3′ UTRs. We also identified 20 intronic motifs with a strong positional bias in the introns surrounding variant exons, and experimentally confirmed the differential expression of three splicing factors corresponding to particular motifs. Collectively, our study provides the first report of the unique mechanisms of alternative splicing controlled by various regulatory elements, and identifies the splicing factors that are associated with cardiac hypertrophy.
Our previous RNA-Seq datasets for cardiac hypertrophy were utilized (Song et al., 2012). Exon variants were analyzed based on RefSeq structure and less expressed exons < 10 reads per kilobase per million reads (RPKM) were discarded. Fisher’s exact test (p < 0.1) and Bayesian error rate (Vencio et al., 2004) (e ≤ 0.1) were used to determine the exon variants. Non-regulated exons were decided by p > 0.5.
Given exon variants, 500 bp intronic sequences (R1-R4), defined as Fig. 2A, were collected. If the intron size was smaller than 500 bp, whole intron was accepted. Considering the previously reported approaches for examining the regulatory motifs by enriched hexanucleotides in a position-dependent manner (Erkelenz et al., 2013; Fairbrother et al., 2002; Hirose et al., 2006; Huh and Hynes, 1994; Kuroyanagi et al., 2013; Lim et al., 2011; Yamashita et al., 2012), we counted the occurrence of each hexanucleotide at the position [ith ∼ i+5th] of the introns. The enrichment of hexanucleotides at position i was analyzed by Fisher’s exact test (p < 0.05). In case of exon inclusion, two independent tests were performed for exon exclusion and non-regulated exons as follows. The enrichment of the hexanucleo-tides at position i of exon inclusion was compared with the number of the same motif at the same position i observed in exon exclusion, and vice versa for exon exclusion.
The enrichment of genes for KEGG pathways was analyzed through Fisher’s exact test (p < 0.05) as follows.
The 8 weeks old male (C57BL/6J) mice (body weight 28?33 g) were used. Cardiac hypertrophy was induced by transverse aortic constriction (TAC) operation under anesthesia with intraperitoneal injection of avertin, 2-2-2 tribromoethanol (Sigma, USA) dissolved in tert-amyl alcohol (Sigma, USA). The procedure of operation was followed as previously described (Cha et al., 2008). As a control group, sham operation (same procedure except for tying) was done. One week after the operation, mice were sacrificed, and hearts were removed and stored in deep freezer at ?80°C. Validation of cardiac hypertrophy model was performed by the measurement of heart weight and body weight ratio and quantification for mRNA expression of hypertrophy markers (ANP and BNP; Fig. 3A).
First-strand cDNA was synthesized from 2 μg of total RNA with Random hexamer using Omniscript? reverse transcription (Qiagen, USA) according to the manufacturer’s instruction. qRT-PCR assays were followed as previously described (Hong et al., 2008). Briefly, qRT-PCR assays were performed using TOPreal™ qPCR premix (Enzynomics, Korea) under the following two-step conditions: denaturation at 95°C for 15 s followed by annealing and extension at 60°C for 40 s, for a total of 40 cycles. The 18S transcript was used as an endogenous reference to assess the relative level of mRNA transcript. The following primers were used:
ANP, 5′-TCGTCTTGGCCTTTTGGC-3′ (sense)
and 5′-TCCAGGTGGTCTAGCAGGTTCT-3′ (antisense)
BNP, 5′-AGGGAGAACACGGCATCATT-3′ (sense)
and 5′-GACAGCACCTTCAGGAGAT-3′ (antisense)
18S, 5′-TTCTGGCGAACGGTCTAGACAAC-3′ (sense)
and 5′-CCACTGGTCTTGGTGTGCTGA-3′ (antisense)
Frozen mouse hearts were homogenized in liquid nitrogen using a mortar and pestle. Homogenized heart powder was lysed in 1 % SDS lysis buffer containing protease inhibitor cocktail, Na3VO4 and NaF. Proteins were estimated by using BCA assay kit (Thermo scientific) and 30 μg of proteins were resolved electrophoretically on 10% SDS polyacrylamide gels, transferred to PVDF membrane. The blots were probed with following antibodies and detected by chemiluminescent substrate (Thermo scientific) and Biomolecular imager, Image-Quant LAS 4000 mini (GE). The band intensities were quantified using NIH ImageJ software.
The following primary antibodies were used for Western blot analysis: ESRP1 (210-301-B89) was purchased from Rockland. PTB, SF2/ASF, SRp30c, and SRp40 antibodies were kind gifts from Dr. Hai-Hong Shen (Gwangju Institute of Science and Technology, Korea).
Our previous study showed the transcriptomic signatures of cardiac hypertrophy and identified hundreds of exon variants that were included or excluded (Song et al., 2012). In the present study, the putative motifs that recruit splicing factors and regulate alternative splicing during cardiac hypertrophy were investigated. For this study, we first collected 3,345 and 3,101 included exons and excluded exons, respectively, and then identified 4,125 that were not regulated (p > 0.5) during pressure-overload cardiac hypertrophy (See “Materials and Methods”). The distributions of alternative splicing types for the exons of interest (i.e., included exons, excluded exons and non-regulated exons) are shown in Supplementary Fig. 1. The results indicate that alternative promoter (p = 7.8 × 10?11), alternative termination (p = 3.7 × 10?22) and cassette exons (p = 0.049) are differentially distributed between exon exclusion and exon inclusion.
The exon variants were characterized in terms of their nucleotide composition, length, and the distribution of 5′- or 3′ UTRs. As shown in Supplementary Fig. 2A, the nucleotide composition was somewhat distinct across the exon variants; an even distribution of A, T, G, and C (ranging from 24% to 25%) was observed in the non-regulated exons, whereas a significantly biased distribution of GC (∼28%) and AT (∼26%) was observed in included exons and excluded exons, respectively (p = 2.2 × 10?16; X2-test). Interestingly, the distribution of 5′ UTRs and 3′ UTRs in the included exons and excluded exons exhibited an obvious dichotomous pattern between; exons containing 5′ UTRs were more likely to be included, whereas exons containing 3′ UTRs tended to be excluded during pressure-overload cardiac hypertrophy (Supplementary Fig. 2B). Longer exons were more likely to be excluded; however, no significant difference was observed in median length across all exons (Supplementary Fig. 2C).
The biased distribution of GC/AT between included and excluded exons and the dichotomous patterns of 5′ UTRs and 3′ UTRs suggest specific architectures (e.g., motifs) that recruit specific splicing factors for either exon inclusion or exon exclusion. To predict the putative motifs in the exon variants, we searched for motifs composed of hexanucleotides that were significantly enriched in an inclusion- or exclusion-specific manner. We found that, among all possible hexanucleotides (i.e., 46 = 4,096), 29 and 24 hexanucleotides were abundantly distributed in excluded exons and included exons, respectively (Supplementary Table 1). Consistent with the high AT-content pattern observed in the excluded exons (Supplementary Fig. 2A), motifs with a high proportion of AT were enriched in excluded exons compared to included exons or non-regulated exons (Fig. 1A). The enrichment of AT-rich motifs implies that AT-rich elements (AREs) contribute to the instability of mRNA in the 3′ UTRs (Barreau et al., 2005). We also found that ‘alternative terminal exons’, defined by the ‘knownAlt’ track of the UCSC genome browser, and 3′ UTRs were significantly enriched in the exons containing AT-rich motifs (p < 0.05), implying the removal of AREs in 3′ UTRs (Fig. 1C and Supplementary Table 1). In contrast, motifs with high GC content were observed in included exons (Fig. 1B and Supplementary Table 1). Consistent with the high proportion of 5′ UTRs observed in the included exons, 5′ UTRs participated in the inclusion of exons with anchoring GC-rich motifs (Fig. 1D). The exon inclusion capability of 5′ UTR content does not appear to be related to an alternative promoter, as no known alternative promoters were significantly enriched (Supplementary Table 1). Rather, it seems to be related to intron retention (e.g., p-value = 8.19 × 10?5?6.16 × 10?8 for GGCCGC, GGCGGC, and CGGCGG; Fisher’s exact test). In light of previous reports that exon borders with higher GC content tended to retain introns (Galante et al., 2004; Wong et al., 2013), this result suggests that the GC-rich motifs located in the 5′ UTRs might be strongly related to intron retention during cardiac hypertrophy.
In summary, we found that exons that contain a GC-rich motif in the 5′ UTRs were significantly included and that AT-rich exons with 3′ UTRs were highly excluded during pressure-overload cardiac hypertrophy. These findings also suggest possible roles for GC-rich and AT-rich motifs in intron retention and exclusion of ARE, respectively.
Similar to what was observed in the exonic motif analysis, the enriched intronic motifs associated with exon variants were also investigated. Due to the importance of proximity to exons, the distribution of all possible hexanucleotides within the 500 base pairs (bp) flanking the intronic sequences (i.e., R1?R4) were analyzed (Fig. 2A). We found 13 and 7 motifs with strong positional bias within the 70 bp adjacent to the introns of excluded exons and included exons, respectively (Figs. 2B and 2C, and Supplementary Table 2). Interestingly, R1 was the most important region, anchoring eight motifs (52%) of the excluded exons. For the included exons, the ∼20 bp upstream of the exons (i.e., R2) were likely to be critical for the binding of splicing factors.
For intronic motifs, we found a correlation between G-rich motifs and intron length. Introns containing G-rich motifs were much shorter than introns containing other motifs. The average intron size in the mouse ranges from 3.8 to 5 kb (Taft et al., 2007; Waterston et al., 2002), whereas the length of the introns upstream and downstream of the exon variants containing G-rich motifs, such as GGAGGG, were less than 1 kb, suggesting a role for these G-rich motifs in the exclusion of exons flanking short introns (Supplementary Fig. 3). Since ESRP1, SF2/ASF, SRp30c, and SRp40 were reported to bind to G-rich motifs (such as GGTGGG, GGAGAG, and GGAGGG) (Supplementary Table 2), we experimentally examined the expression level of splicing factors in Sham/TAC-operated heart sample and foundthe differential expression of SF2/ASF and ESRP1 (Fig. 3). G-rich intronic element at 5′ splice site was known to regulate alternative splicing of thyroid hormone receptor by SF2/ASF (Hastings et al., 2001) and we also found the significantly enriched G-rich motif at 5’ splice site along with significantly increased level of SF2/ASF during cardiac hypertrophy. G-rich motifs are also involved in the regulation of alternative splicing by ESRP (Dittmar et al., 2012; Warzecha et al., 2010). Of the G-rich motifs, GTGGG-containing motifs (i.e., GGTGGG and GTGGGG) were significantly enriched in multiple sites for R1 and R3 of 57 exons, which were highly involved in cellular movement (p = 7.62 × 10?5; Fisher’s exacttest). ESRP is known to bind to TGG motifs, which are located within these GTGGG motifs, and play an important role in cellular movements such as cell adhesion, cell polarity, and cell migration (Warzecha et al., 2010). Interestingly, the GTGGG motifs were found along with TGCCTG motifs in the adjacent 50 bp of R2 (p = 1.49 × 10?4), suggesting close interplay between these two motifs in exon exclusion. The co-occurrence of TGCCTG at the 3′ splicing acceptor of the exons (i.e., R2) and GTGGG at the 5′ splicing donor of excluded exons (i.e., R3) were observed in Shank3, Vps16, Rhbdl1, Arhgef2, Nbeal2, and Brd9 (data not shown). TGC/CTG repeats are known binding motifs for CELF proteins, which are critical splicing factors during cardiac development (Ladd et al., 2001). Besides the G-rich motifs, we also found the significantly increased PTB known to bind to TCTT motif, during cardiac hypertrophy (Das et al., 2007).
Taken together, we predicted that intronic motifs were significantly enriched at specific positions in introns, and experimentally confirmed the expression of the splicing factors that bind to these motifs. The role of ESRP1 is notable, inducing the exclusion of exons flanking a short intron in genes involved in cellular movement.
During pressure overload cardiac hypertrophy, the heart undergoes rapid enlargement, remodeling and irreversible maladaptation accompanied by dramatically increased energy consumption due to rapid cardiac growth and the synthesis of cellular building blocks. Therefore, systematic cellular and metabolic transitions are necessary to meet the substantial demands during cardiac remodeling. Even though the quantitative signatures of pathological cardiac hypertrophy have been well studied (Park et al., 2011; Song et al., 2012), the qualitative changes during cardiac remodeling have remained to be investigated. Here, we show the evidence for the association of alternative splicing with pressure-overload hypertrophy as follows: 1) GC-rich exonic motifs regulate intron retention in 5′ UTRs and AT-rich exonic motifs are involved in exclusion of the AT-rich elements in 3′ UTRs during cardiac hypertrophy (Figs. 1 and 4, and Supplementary Table 1); 2) 20 putative splicing factor binding motifs enriched in intronic regions were identified (Fig. 2 and Supplementary Table 2); 3) Splicing factors such as PTB and ESRP1 that are known to bind to the putative motifs were differentially expressed during cardiac hypertrophy (Fig. 3); 4) Collectively, the splicing regulatory elements are involved in post-transcriptional regulation via the 3′ or 5′ UTRs or by recruiting splicing factors during cardiac hypertrophy (Fig. 4).
As summarized in Fig. 4A, alternative splicing for both 5′ and 3′ UTRs may play an important role in regulating translation during cardiac remodeling. The observed retention of introns in 5′ UTRs is important, because the mechanism could be involved in post-transcriptional regulation by inserting regulatory elements such as internal ribosome entry site (IRES) and G-quadruplex into mRNA. According to the previous reports, the intron retention led to the acquisition of IRES in Zeb2 and to the insertion of G-quadruplex in p53 (Beltran et al., 2008; Marcel et al., 2011). In the case of alternative splicing for 3′ UTR, we observed the removal of AREs that make the mRNA unstable. These findings suggest that post-transcriptional regulation by alternative UTRs might increase mRNAstability, thereby maximizing protein synthesis without significantly increasing transcription. Based on the rapid increase in heart volume during cardiac hypertrophy, this mechanism is an excellent way to adapt against urgent cardiac overload.
As shown in Fig. 4B, we found the altered expression of splicing factors that bind to the intronic motifs. Among them, PTB and ESRP1 are already known to be involved in alternative splicing of smooth muscle and metastasis in cancer, respectively (Gooding et al., 2013; Reinke et al., 2012; Yae et al., 2012), but their roles in cardiac hypertrophy are not known. They likely bind to the G-rich regions and are widely distributed nearby exons excluded (Supplementary Table 2). Warzecha et al. (2010) reported that UGG is a core motif of ESRP1 and the motif located in upstream intron possibly promotes exon exclusion, whereas the motifs in downstream of the regulated exon enhance exon inclusion. Their following study on the ESRP1 binding motifs using SELEX-Seq predicted the enriched GU-rich motifs in the flanking upstream (R2) of the excluded exons. For the included exons, they found the motifs significantly enriched in the distal upstream (R1) and flanking downstream (R3). During cardiac hypertrophy, the decreased level of ESRP1 might result in the weak binding at R1 and R3 as reported recently (Dittmar et al., 2012).
Intronic GU-rich repeats at 5′ splice site are also known to be the binding region for SF2/ASF and hnRNP H (Arrisi-Mercado et al., 2004; Hastings et al., 2001). Similarly, we also found the significantly increased level of SF2/ASF during cardiac hypertrophy. However, the detailed mechanistic relationship between the differential expression of the splicing factors (e.g., ESRP1 and SF2/ASF) and the G-rich motifs remains to be elucidated, as the splicing factors usually play roles as both enhancer and repressor depending on the position, concentration and the location of the branchpoint (Lim et al., 2011; Taggart et al., 2012; Zuo and Manley, 1993).
Consistent with the known role of ESRP1 in the alternative splicing of genes involved in cellular movement such as cell adhesion, polarity and migration (Warzecha et al., 2010), we also identified 57 genes with the ESRP1 binding motif that are involved in cellular movement. The significantly decreased expression of ESRP1 (Fig. 3) also suggests the importance of ESRP1 in alternative splicing during cardiac hypertrophy. We also identified a significant co-occurrence of the binding motifs for ESRP1 and CELF. However, the expression of the gene encoding CELF did not change during cardiac hypertrophy (data not shown). Although there was no change in the expression of CELF genes (e.g., Cugbp1, Cugbp2, Brunol4, and Brunol5), CELF expression can possibly be regulated post-transcriptionally through microRNAs (Kalsotra et al., 2010).
Here, we report the splicing regulatory elements and the putative splicing factors that may play important roles for pathological cardiac hypertrophy. The molecular mechanisms for alternative splicing regulated by the splicing regulatory elements and the functional roles of the individual splicing factors identified from this study will await the future molecular studies.
SUPPLEMENTARY FILE
Hye Kyeong Kwon, Hyunwoo Choi, Sung-Gyoo Park, Woo Jin Park, Do Han Kim, and Zee-Yong Park*
Mol. Cells 2021; 44(7): 500-516 https://doi.org/10.14348/molcells.2021.4002Jeesun Chun, Kum-Kang So, Yo-Han Ko, Jung-Mi Kim, and Dae-Hyuk Kim
Mol. Cells 2019; 42(4): 363-375 https://doi.org/10.14348/molcells.2019.0019Ki-Beom Moon, Dong-Joo Ahn, Ji-Sun Park, Won Yong Jung, Hye Sun Cho, Hye-Ran Kim, Jae-Heung Jeon, Youn-il Park, and Hyun-Soon Kim
Mol. Cells 2018; 41(11): 979-992 https://doi.org/10.14348/molcells.2018.0312