Skip to main content

Comparative expression patterns and diagnostic efficacies of SR splicing factors and HNRNPA1 in gastric and colorectal cancer

Abstract

Background

Serine/arginine-rich splicing factors (SRSFs) and HNRNPA1 have oncogenic properties. However, their proteomic expressions and practical priority in gastric cancer (GC) and colorectal cancer (CRC) are mostly unknown. To apply SFs in clinics, effective marker selection and characterization of properties in the target organ are essential.

Methods

We concurrently analyzed SRSF1, 3, and 5–7, and HNRNPA1, together with the conventional tumor marker carcinoembryonic antigen (CEA), in stomach and colorectal tissue samples (n = 420) using semiquantitative immunoblot, subcellular fractionation, and quantitative real-time polymerase chain reaction methods.

Results

In the semiquantitative immunoblot analysis, HNRNPA1 and SRSF7 levels were significantly higher in GC than in gastric normal mucosa, and SRSF7 levels were higher in intestinal-type compared with diffuse-type of gastric adenocarcinoma. Of the SFs, only HNRNPA1 presented greater than 50 % upregulation (cancer/normal mucosa > 2-fold) incidences and CEA-comparable, acceptable (>70 %) detection accuracy (74 %) for GC. All SF protein levels were significantly higher in CRC than in colorectal normal mucosa, and HNRNPA1 levels were higher in low-stage CRC compared with high-stage CRC. Among the SFs, HNRNPA1 and SRSF3 presented the two highest upregulation incidences (88 % and 74 %, respectively) and detection accuracy (90 % and 84 %, respectively) for CRC. The detection accuracy of HNRNPA1 was comparable to that of CEA in low (≤ II)-stage CRC but was inferior to that of CEA in high (>II)-stage CRC. Extranuclear distributions of HNRNPA1 and SRSF6 (cytosol/microsome) differed from those of other SRSFs (membrane/organelle) in both cancers. In an analysis of the six SF mRNAs, all mRNAs presented unacceptable detection accuracies (≤70 %) in both cancers, and all mRNAs except SRSF6 were disproportionate to the corresponding protein levels in GC.

Conclusion

Our results provide a comprehensive insight into the six SF expression profiles in GC and indicate that, among the SFs, HNRNPA1, but not HNRNPA1 mRNA, is the most effective, novel GC marker. Regardless of the good to excellent detection accuracy of SRSF3 and HNRNPA1 in CRC, the SFs have lower practical priority than CEA, especially for high-stage CRC detection.

Peer Review reports

Background

Alternative splicing is a ubiquitous post-transcriptional process that leads to proteomic diversity and the disruption of splicing regulatory networks is a critical component of carcinogenesis.

The serine/arginine (SR) protein family is an important class of splicing regulators and its members, including SR splicing factor (SF) 1, SRSF3, and SRSF6, have shown multiple proto-oncogenic properties and aberrant expressions in various cancer cells [14]. SRSF1, 3, and 7 shuttle between nucleus and cytoplasm, and their subcellular distribution is associated with various cellular functions and reactive responses [57]. For instance, in the nuclear compartment of cells, SFs promotes the splicing process, while in the extranuclear compartment, they regulate protein translation [5]. In this sense, SFs are found in ribosomes, the translation machinery [5], and are expected to colocalize with SF-binding translation regulatory proteins such as the target of rapamycin complex 1 (TORC1), which locates in lysosomes or in cytosol [6]. Additionally, specific cellular stresses or conditions induce post-translational modifications of SFs which associated with cytoplasmic localization and the stability of SFs and/or inhibit general splicing process in nucleus [58].

Heterogeneous nuclear ribonucleoprotein A1 (HNRNPA1) is also a well-known splicing regulator with effects antagonistic to SR proteins [9]. Upregulated expression and aberrant cytoplasmic localization of HNRNPA1, as determined by immunohistochemical staining, were noted in colorectal cancer (CRC) [10]. Recently, HNRNPA1 has emerged as a plausible biomarker of CRC [11].

However, in gastric cancer (GC), proteomic expressions of these SFs are unclear and their diagnostic values have not been determined. The comparative practical priority among the SR proteins and conventional tumor markers in GC and CRC is unknown. Moreover, it is unclear as to whether the expression patterns and practical priority of these SFs differ between the cancers.

We aimed to: (1) identify the translational and transcriptional profiles of five SRSFs (SRSF1, SRSF3, SRSF5, SRSF6, and SRSF7) and HNRNPA1 in GC; and (2) compare the detection accuracy (DA) among the SFs and the currently used tumor marker carcinoembryonic antigen (CEA) [12] and determine specific features in these parameters in GC and CRC.

Methods

Subjects and sample preparation

A total of 420 fresh stomach and colon biopsy samples were obtained from 224 patients (Table 1) who had undergone surgical resection.

Table 1 Demographic characteristics of subjects and tumors

All tumor samples and patients were defined and diagnosed, respectively, at the Department of Pathology and Surgery, Wonkwang University Hospital. Each fresh biopsy sample was aliquoted into two or four tubes. The first sub-sample was stored frozen in liquid nitrogen until immunoblot analysis; the second was immediately fixed with formalin, paraffin embedded, then stored for hematoxylin and eosin staining and immunohistochemical staining; the third was immediately homogenized for subcellular fractionation; and for the fourth, RNA was immediately extracted and the RNA was reverse transcripted. The synthesized cDNA sample was stored frozen at −75 °C until real-time polymerase chain reaction (PCR) analysis. All paraffin-embedded tissue samples were sectioned, stained with hematoxylin and eosin and evaluated twice by two pathologists. All histological findings were consistent with the diagnosis, and all were concordant between the two pathologists. All sections of cancer tissue included at least 30 % cancer cells, and none showed light microscopically-detectable degeneration or necrosis. Aneuploidy status was determined using routine clinical laboratory diagnostic methods consisting of propidium iodide staining followed by flow cytometric analysis. The SW1116 cell line (a CEA-producing colon cancer cell line) was obtained from the American Type Culture Collection and used between passage 3 and 8 after they were obtained.

Immunoblot and semiquantitative analysis of SF proteins and CEA

Tissue samples (40 mg) were homogenized in RIPA buffer containing protease inhibitors using a Bullet Blender homogenizer (Next Advance; Averill Park, NY, USA), and whole-cell lysate was obtained by sequential centrifugations. Proteins (~30 μg) in the whole-cell lysate were separated on 10 % sodium dodecyl sulfate–polyacrylamide gels with SW1116 cell protein (~30 μg). The proteins were transferred onto polyvinylidene difluoride membranes. The membranes were blocked with 5 % skim milk in Tris-buffered saline containing 0.1 % Tween 20 (TBS-T), rinsed, and incubated with the appropriate antibodies in TBS-T containing 3 % skim milk. Excess primary antibody was then removed by washing the membrane four times in TBS-T. The membranes were then incubated with horseradish peroxidase-conjugated secondary antibody (anti-rabbit or anti-mouse). After three washes in TBS-T, bands were visualized using Clarity Western ECL Substrate (Biorad; Hercules, CA, USA) on the FluorChem E System (Protein Simple; Santa Clara, CA, USA). The following primary antibodies were used in the immunoblot analysis of whole- or fractionated cell lysates: anti-HNRNPA1 (dilution 1:2000; cat. sc-32301; Santa Cruz Biotechnology; Santa Cruz, CA, USA); anti-SRSF1 (dilution 1:1000; cat. 324600; Invitrogen; Carlsbad, CA, USA); anti-SRSF3 (dilution 1:1000; cat. RN080PW; MBL; Nagoya, Japan); anti-SRSF5, anti-SRSF6, anti-SRSF7 (dilution 1:1000; cat. HPA 043484, HPA029005, HPA043850; Sigma-Aldrich; St. Louis, MO, USA); anti-CEA (dilution 1:3000; cat. MS-613-P; Thermo Fisher Scientific; Fremont, CA, USA); anti-ACTB (dilution 1:5000; cat. MA5-15739; Invitrogen); anti-poly (ADP-ribose) polymerase (PARP; dilution 1:1000; cat. sc-7150; Santa Cruz Biotechnology); anti-histone H3 (dilution 1:1000; cat. 9715 s; Cell Signaling; Beverly, MA, USA); and anti-prohibitin (dilution 1:1000; cat. AB28172; Abcam; Cambridge, UK).

Then, the optical density of the region of the target molecular weight (±20 %), as described by the respective antibody manufacturers, was analyzed using ImageJ (http://imagej.nih.gov/ij/). Relative protein levels of the samples were determined after normalization to the β-actin band and calibrated using bands from SW1116 cells (a value of 10 was assigned to the SW1116 cell) on each membrane. As presented in previous reports [1315] and/or manufacturers’ instructions, SRSF1, SRSF6, SRSF7, HNRNPA1, and CEA showed multiple bands in each target molecular weight area. For the five proteins, multiple bands in each target molecular weight area were calculated and the sum of the distinct band/bands was then analyzed for the proteins. Values of 0.2 (HNRNPA1 and CEA), 0.7 (SRSF1), 0.8 (SRSF3, SRSF5, and SRSF7) or 1.0 (SRSF6) were assigned to the undetected bands of target regions based on band density of SW1116 cell lysate.

Subcellular fractionation and immunoblot

Ten paired [normal mucosa (NM) and cancer] samples (n = 20) were fractionated using a modification of the method described previously [1618]. Tissue samples (70 mg) were homogenized in 800 μL of homogenization buffer [0.25 M sucrose, 10 mM EDTA, 10 mM EGTA, 2 mM MgCl2, 20 mM Tris–HCl (pH 7.4) and protease inhibitors]. The homogenized lysate was centrifuged at 1500 × g for 5 min at 4 °C, and the resulting pellet, containing nuclei, and the supernatant (post-nuclear supernatant) were separated. Subsequently, the pellet, containing nuclei, was resuspended in nuclear extraction buffer [2.5 % glycerol, 1 mM EDTA, 1 mM EGTA, 1.5 mM MgCl2, 0.42 M NaCl, and 20 mM HEPES (pH 7.6)] for 1 h at 4 °C, centrifuged again at 20000 × g for 30 min at 4 °C. The supernatant was then collected and referred to as the nuclear extraction fraction. In parallel, the post-nuclear supernatant was centrifuged at 20000 × g for 30 min at 4 °C and the resulting supernatant was saved as the cytosol/microsome fraction. The pellet was resuspended in RIPA buffer for 10 min at 4 °C followed by centrifugation at 1500 × g for 5 min at 4 °C, and the resulting supernatant was collected and referred to as the membrane/organelle fraction, including intracellular membranes and some plasma membrane. Immunoblot analysis was performed on these three fractions (nuclear extract, membrane/organelle, and cytosol/microsome) as described above.

The resultant three fractions differed in terms of the total protein amounts (membrane/organelle < nuclear extract < cytosol/microsome) for all tissue samples. In addition, the subcellular fraction indicator proteins (PARP, histone H3, and prohibitin) showed inter-individual and intra-individual (cancer vs. NM) variations. Accordingly, to compare the relative SF expression levels for the respective fractions, an equal amount of total protein per lane was loaded.

Enzymatic immunohistochemical staining

Formalin-fixed, paraffin-embedded tissue was sectioned and placed on slides. Sections were stained using the Discovery XT automated immunohistochemistry (IHC) stainer (Ventana Medical Systems; Tucson, AZ, USA) and Ventana Chromo Map Kit (Ventana Medical Systems) according to the manufacturer’s instructions. The sections were deparaffinized using EZ prep solution; the antigen was retrieved in cell a conditioning solution (CC1; Ventana Medical Systems) under standard conditions (100 °C, 60 min), and endogenous peroxidase was inhibited by treatment with 3 % H2O2 for 4 min. Then, the sections on slides were incubated with the primary antibody that had been used in the immunoblot analysis (dilution 1:300 for anti-HNRNPA1 and anti-CEA and dilution 1:100 for other antibodies) for 60 min at 37 °C, and then with a secondary antibody (UltraMap anti-RB HRP or UltraMap anti-MS HRP; Ventana Medical Systems) for 28 min at 37 °C. The sections were incubated in diamidobenzidine and H2O2 for 8 min at 37 °C followed by counterstaining with hematoxylin and treatment with bluing solution. On completion of staining, sections were dehydrated in alcohol, cleared in xylene, and mounted in synthetic resin.

mRNA quantification using quantitative real-time PCR

Total RNA from 40 mg of tissue sample was isolated using 1 mL of TRIzol (Life technologies; Carlsbad, CA, USA) in accordance with the manufacturer’s instructions. RNA (500 ng) was reverse transcribed using ReverTra Ace qPCR RT Kits (Toyobo; Osaka, Japan). Quantitative real-time PCR was performed in a StepOnePlus Real-Time PCR System (Applied Biosystems; Foster City, CA, USA) using a SYBR Green Realtime PCR Master Mix (Toyobo) in accordance with the manufacturer’s instructions. Each assay was performed in triplicate, and results were normalized to 18S rRNA levels. Relative mRNA level was calculated using StepOne software v.2.2.2 and calibrated to that of SW1116 cell sample (A value of 10 was assigned to the mRNA level of SW1116 cell) in each batch test.

Primer sequences were as follows: 5'-TGGATTTGGTAATGATGGAA-3' and 5'-TCTCTGGCTCTCCTCTCCTG-3' (HNRNPA1); 5'-TGCCTACATCCGGGTTAAAG-3' and 5'-CTGCTGTTGCTTCTGCTACG-3' (SRSF1); 5'-TCTTGGAAACAATGGCAACA-3' and 5'-CTCGGGGATCTTCAAATTCA-3' (SRSF3); 5'-GAGGCTTTGGTTTTGTGGAA-3' and 5'-CGAGCCCTAGCATGTTCAAT-3' (SRSF5); 5'-AAATACGGACCACCTGTTCG-3' and 5'-CTTCACCTGCTTGTCGCATA-3' (SRSF6); 5'-CGCTGGCAAAGGAGAGTTAG-3' and 5'-CGAATTCCACAAAGGCAAAT-3' (SRSF7); 5'-GTAACCCGTTGAACCCCATT-3' and 5'-CCATCCAATCGGTAGTAGCG-3' (18S rRNA).

Statistical analysis

All group values were non-normally distributed in the Kolmogorov-Smirnov test (p < 0.05). Thus, the values were compared using the Mann–Whitney U-test (between two groups) or Kruskal-Wallis test (among more than two groups). If values were significantly different in the Kruskal-Wallis test, Conover’s post-hoc tests were performed. Ratios were compared using the chi-square test for comparison of incidences between GC and CRC. Spearman correlations were used to assess relationships among the levels of SF proteins or CEA and levels of SF mRNAs. The DAs of each SF for discriminating between cancer and NM were obtained by constructing receiver operating characteristic curves. The cut-offs for defining cancer were determined by the respective highest Youden indexes. The statistical differences between the DAs were determined using the DeLong method (for same-sample-derived comparison) or a comparison of the areas under independent receiver operating characteristic curves (for two independent sample-derived comparisons).

Data were analyzed with MedCalc version 12.7 (MedCalc Software, Mariakerke, Belgium) and StatsDirect version 2.7.8 (StatsDirect Ltd, Cheshire, UK). A two-tailed p value of less than 0.05 was considered statistically significant.

Results

Relative protein levels of the SFs and CEA in GC and CRC

In a comparison of relative immunoblot band intensity (Fig. 1a and b) of stomach samples, the median levels of HNRNPA1 (2.2 vs. 0.5; 4.4-fold difference; p < 0.001), SRSF7 (4.8 vs. 3.0; 1.6-fold difference; p = 0.006), and CEA (2.6 vs. 0.6; 4.3-fold difference; p < 0.001) were significantly higher in GC samples than in gastric NM (GNM) samples (Fig. 1b). Other SF protein median levels were not significantly (p > 0.05) different between GC and GNM samples.

Fig. 1
figure 1

Proteomic expression of six SFs and CEA. a Representative semiquantitative immunoblot analysis of gastric and colorectal (CR) cancer (C) and their adjacent normal mucosa (NM) tissue samples. β-actin (ACTB) was used as a loading control. SW1116 cell lysate was loaded as a semiquantitative calibrator. For each SF and CEA, distinct bands in the following target molecular weights were analyzed: 29–39 kDa (HNRNPA1), 27–33 kDa (SRSF1), 20 kDa (SRSF3), 36–40 kDa (SRSF6), 33–15 kDa (SRSF7), and 180–220 kDa (CEA). b Relative quantitation of SF proteins and CEA in immunoblot analysis for gastric normal mucosa (GNM), gastric cancer (GC), colorectal NM (CRNM), CR cancer (CRC) and CR-adenoma. Each bar indicates the median and the green horizontal lines indicate the interquartile ranges. The color (compared group) matched small blocks above each median bar indicate the p value < 0.05 (by Kruskal-Wallis test followed by Conover’s post-hoc tests). c Representative subcellular distribution analysis using biochemical fractionation and immunoblot method. Whole cell (wc); nuclear extract (nu); cytosol/microsome (cyt); and membrane/organelle (me/og) fractions were loaded. The nuclear extract fraction was identified with poly (ADP-ribose) polymerase (PARP) and histone H3 (H3) and membrane/organelle fraction was identified with prohibitin. Respective subcellular fractions were normalized using the total protein quantity/ACTB. The actual experiment involved 10 paired samples (n = 20; 5 GC/GNM pairs and 5 CRC/CRNM pairs). d Representative immunohistochemical staining analysis. The actual experiment involved 10 paired samples (n = 20; 5 GC/GNM pairs and 5 CRC/CRNM pairs)

In immunoblot analyses of CR samples (Fig. 1a and b), the median levels of all six SFs and CEA were significantly higher in CRC than in CRNM. Also, the median levels of all SFs except SRSF1 were significantly higher in CR-adenoma than in CRNM. HNRNPA1 and SRSF3 median levels in CRC were most markedly higher than in CRNM samples among the SFs. Median levels of the six SFs were not lower or even higher in CR-adenoma samples compared with CRC samples, but the median level of CEA was lower in CR-adenoma than in CRC.

When GC and CRC were compared (Fig. 1b), median levels of all SFs except SRSF1 were significantly higher in CRC than in GC samples.

Association between the clinicopathologic factors and SF protein levels

In an association analysis between the clinicopathologic factors and SF proteins (Table 1), none of the SF protein levels were significantly different based on patient age, gender, or cancer location in both cancers. Most GCs and CRCs were adenocarcinomas; in gastric adenocarcinoma, the median level of SRSF7 (7.0 vs. 2.9; p = 0.049) was significantly higher in intestinal-type than in diffuse-type. In CR-adenocarcinoma, the median level of HNRNPA1 was significantly higher in well-differentiated type (7.0 vs. 5.0; p = 0.020), no lymph node metastasis (7.0 vs. 4.3; p = 0.003), or low (≤II)-stage (7.0 vs. 4.3; p = 0.003) groups than in the other groups. Median levels of other SF proteins, except HNRNPA1 and SRSF7, were not significantly different based on tumor status, lymph node metastasis, TNM stage, or aneuploidy status in both gastric and CR adenocarcinoma.

Comparison of upregulation incidences of the SF and CEA proteins in GC and CRC

In paired sample (cancer/NM from each same patient) comparison (Table 2) of stomach samples, only HNRNPA1 showed >50 % upregulation (cancer/NM >2) incidence (UI), followed by the other UIs in the order of HNRNPA1 (52 %), SRSF7 (42 %), and other SF proteins (22 %–33 %). In CRC, all SFs except SRSF1 and SRSF5 showed UIs of >50 % in the order of HNRNPA1 (88 %), SRSF3 (74 %), and other SFs (31 %–57 %). When UIs in GC and CRC were compared, the UIs of all SFs were lower (SRSF3, HNRNPA1, and SRSF6) or tended to be lower in GC than in CRC; the differences between GC and CRC were greatest for SRSF3 (40 %; p < 0.001), followed by HNRNPA1 (36 %; p < 0.001), SRSF6 (24 %; p = 0.004), and the other SFs. The UIs of CEA were >50 % in both cancers; it was also lower in GC (53 %) compared with CRC (92 %), and the results were similar to HNRNPA1 values.

Table 2 Comparison of upregulation incidences and detection accuracies of SFs and CEA in GC and CRC

Comparison of DAs of the SF and CEA proteins in GC and CRC

In the DA analysis (Table 2) for GC, only DA-HNRNPA1 (74 %) was acceptable (>70 %), while the other SFs presented poor (≤70 %) DAs. The DAs of all SFs were not significantly different between low- and high-stage gastric adenocarcinoma. In CRC, all SFs except SRSF1 and SRSF5 presented acceptable DAs in the order of HNRNPA1 (90 %), SRSF3 (84 %), and the other SFs (62–76 %). When the DAs for GC and CRC were compared, the DAs of all SFs were lower (SRSF3, HNRNPA1, and SRSF6) or tended to be lower in GC than in CRC; the differences between GC and CRC were greatest for SRSF3 (32 %; p < 0.001), followed by HNRNPA1 (17 %; p < 0.001), SRSF6 (14 %; p = 0.006), and the other SFs. The DAs of all SFs were not significantly different based on the stage of CR-adenocarcinoma, and the greatest difference between GC and CRC being in SRSF3 was consistently observed regardless of the stage. In the comparison with CEA (Table 2), DA-HNRNPA1 was not different from DA-CEA in both low- and high-stage gastric adenocarcinoma, whereas in CRC, DA-HNRNPA1 was not different from DA-CEA in low-stage CR-adenocarcinoma, but it was significantly lower than DA-CEA in high-stage CR-adenocarcinoma.

Subcellular distribution of SF proteins in GC and CRC

Previously, several studies showed SF distribution using enzymatic IHC [7].

Generally, this method provides information about both cellular morphology and tissue structural changes, but it is limited when used for a clear discrimination of organelles (the major location of translation regulatory molecules) from cytosol. So far, no study has demonstrated the SF subcellular distribution in primary cancer cells using biochemical fractionation and immunoblotting. Using the biochemical fractionation approach (Fig. 1c), we were able to more precisely distinguish the extranuclear fractions. HNRNPA1 and SRSF6 in both cancers frequently (≥50 %) showed different distribution patterns from those of the other SFs. When HNRNPA1 and SRSF6 were upregulated in whole cancer-cell lysates, they were predominantly distributed in nuclear and/or cytosol/microsome fractions; the upregulated SFs in the cytosol/microsome fraction were detectable in >50 % of upregulated cases in both cancers. Alternatively, when SRSF1, 3, 5, and 7 were upregulated in whole cancer-cell lysates, they were predominantly distributed in nuclear and/or membrane/organelle fractions; the upregulated SFs in the membrane/organelle fraction were detectable in >50 % of upregulated cases in both cancers. There were no remarkable differences in subcellular distributions between GC and CRC. With regard to methodological principle and analytical points, biochemical fractionation followed by immunoblotting differed from IHC. The former was normalized using equal amounts of total protein per lane, while IHC images were analyzed directly. During the former method process, the nuclear and extranuclear fractions were transiently treated in a different buffer, which may affect the affinity of the antibody, and a certain amount of the protein in the border zone of the different fractions was discarded to prevent contamination of a adjacent fraction. Regarding these concerns, we only were able crudely to match the nuclear and extranuclear fraction quantities, and were not able to precisely match the distribution of four SF proteins (SRSF 1, 3, 5, and 7) in the membrane/organelle fraction to the corresponding fraction of the enzymatic IHC analysis (Fig. 1d). Nonetheless, the two methods showed an approximate match with regard to whole-cell SF protein intensities.

Relative mRNA levels of the SFs in GC and CRC

In the SF mRNA level analysis (Fig. 2a), none of the six SF median levels significantly differed between GC and GNM. In contrast, all SFs except SRSF5 were significantly higher in CRC than in CRNM.

Fig. 2
figure 2

Relative mRNA levels of six SFs (a) and correlation analysis among the SF proteins and mRNAs (b). a Each bar indicates the median, and the green horizontal lines indicate the interquartile ranges for gastric normal mucosa (GNM), gastric cancer (GC), colorectal NM (CRNM), CR cancer (CRC), and CR-adenoma. The color (compared group) matched small blocks above each median bar indicate the p value < 0.05 (by Kruskal-Wallis test followed by Conover’s post-hoc tests). b R1, S1, S3, S5, S6, and S7 indicate HNRNPA1, SRSF1, SRSF3, SRSF5, SRSF6, and SRSF7, respectively. The numbers in colored cells (p < 0.05) or uncolored cells (p > 0.05) indicate the correlation coefficient (r)

None of the SF mRNA median levels were significantly different based on histological features, lymph node metastasis, or TNM stage in gastric or CR-adenocarcinoma (Table 1).

In paired sample comparisons of mRNA (Table 2), all SF mRNAs showed <50 % UI in GC, and only HNRNPA1 showed >50 % UI in CRC.

In DA analysis of mRNA (Table 2), including DA-HNRNPA1 (69 %) for CRC, which was the highest SF mRNA DA, all of the SF mRNAs showed poor DAs for both cancers. The UIs and DAs of all of the SF mRNAs were not significantly different between GC and CRC (Table 2).

Correlation between the SF proteins and their respective mRNA levels

In correlation analysis (Fig. 2b), all SF proteins were significantly correlated with each other and all SF mRNAs were correlated with each other in both stomach and CR samples.

None of the SF protein levels, except SRSF6, were correlated with their respective mRNA levels in stomach samples. Whereas in CR samples, levels of three SF proteins (HNRNPA1, SRSF3, and SRSF6), which had relatively high DAs, were correlated with their respective mRNA levels.

Discussion

Our study demonstrated, for the first time, the following: (1) HNRNPA1 is the most useful, CEA-comparable marker for GC among the six SF proteins and their mRNAs, and HNRNPA1 and SRSF7 have significantly elevated levels in GC tissue compared with GNM; (2) comparative diagnostic priority of HNRNPA1 for CRC is dependent on stage (i.e., DA-HNRNPA1 for CR-adenocarcinoma is comparable to that of CEA in low-stage, but inferior to high-stage, CR-adenocarcinoma) (3) unlike in GC, SRSF3 presents relatively high UI and very good DA in CRC; and (4) extranuclear distribution patterns of HNRNPA1 and SRSF6 (cytosol/microsome) differ from those of the other SFs (membrane/organelle) in both cancers.

Although strong evidence supports the role of several SFs in tumorigenesis [14, 14] and suggests their potential as diagnostic markers of cancers [11, 19], to apply the SFs in clinics, identification of their DAs and specific characteristics in the target organ is essential. In this regard, we concurrently compared the six SFs and found that only the levels of SRSF7 and HNRNPA1 are statistically significantly higher in GC than in GNM. Currently, no reports have shown, with statistical significance, elevated levels of SRSF7 in any type of primary cancer or elevated levels of HNRNPA1 in primary GC samples. Our results present robust evidence for the relationship between two SF proteins and GC using a statistically valid sample number. In addition, we found significantly higher levels of six SF proteins in CRC than in CRNM. Previously, a study reported amplified SRSF6 expression in primary CRC [3], but the data was limited to the gene level. It is difficult to find any report that shows significantly elevated levels of the five SRSF proteins in primary CRC. Here, we provide new evidence of a positive relationship between the five SRSFs and CRC in terms of the protein levels using primary CRC samples.

Even though most GC and CRC originate from gastrointestinal tract mucosa and consist of adenocarcinomas, when the SF protein levels were compared between GC and CRC, CRC generally tended to show higher levels than GC. In addition, UI and DA analyses demonstrated significant differences between GC and CRC for HNRNPA1, SRSF3, and SRSF6, and the difference in SRSF3 was most marked. Although it is difficult to pinpoint the exact reason for this, we postulated that it is related to the previously known different etiologies and pathogenetic factors of GC and CRC [2022]. These different findings for GC and CRC also hold true for our CEA levels in GC and CRC.

None of the SRSF 5–7 and HNRNPA1 proteins showed significant relationships with unfavorable histopathologic or aneuploidy status [23] in gastric adenocarcinoma or CR-adenocarcinoma. Instead, SRSF7 in gastric adenocarcinoma was higher in intestinal-type, which has a higher survival rate than diffuse-type [24], and HNRNPA1 levels in CRC were higher in well-differentiated, without nodal involvement, or low-stage groups. Furthermore, all SF protein levels in CR-adenoma were either not different from or even higher than those of CRC, which has a poorer prognosis than adenoma. Collectively, the elevated levels of all SFs, especially the two SF proteins, do not seem to correlate with poor prognostic factors in either cancer. For the association between HNRNPA1 and prognostic factors of CRC, previous reports showed contradictory findings. One presented a higher upregulation incidence of HNRNPA1 in low-stage group than in high-stage group [25], while another [11] presented a lower survival rate in patients who had higher HNRNPA1 levels in CRC tissue. Our HNRNPA1 result was in line with the former finding. As a concurrent comparison study of the diagnostic performance levels of multiple SF proteins, our study was designed in a cross-sectional view for selection of a diagnostic marker rather than a longitudinal view for identification of a prognostic/predictive marker, which generally requires a 3–5-year follow-up period. In this respect, although we do not formally report inconclusive results at this point, we believe that our results are also in line with the former one, since in our preliminary survival analysis for 0.6–2.5 years, patients with higher HNRNPA1 levels tended to have longer relapse-free intervals and overall survival times.

Considering the UIs and DAs among the six SFs, HNRNPA1 was the most effective marker for both cancers. Analytical results in many cancer databases, including The Cancer Genome Atlas (TCGA) data [26, 27], partly supported our findings. Specifically, the SF mRNA median levels in the TCGA data (<1.9-fold) were only marginally upregulated in both cancers as in our mRNA results. In contrast, in this study, HNRNPA1 protein levels presented remarkably elevated levels (>4.4-fold) in both cancers. Together, not only the direct comparison within our data (protein vs. mRNA) but also the indirect comparison with the results of the cancer database (marked fold change of HNRNPA1 in this study vs. marginal fold change of HNRNPA1 in the database), to a degree, validated the relative excellence of HNRNPA1 compared with HNRNPA1. In addition, because the DA-HNRNPA1 for GC and low-stage CRC were comparable with the DA of CEA, one of the most commonly used tumor markers for gastrointestinal cancers in clinical practice, our results indicate that HNRNPA1 is applicable as a novel GC marker as well as an early CRC marker. However, in high-stage CR-adenocarcinoma, DA-HNRNPA1 was significantly inferior to DA-CEA. For this, we assumed that HNRNPA1 production reaches a plateau during the premalignant period in the course of CRC progression. Alternatively, it may be due to microenvironmental factors or the instability of HNRNPA1 in a large tumor mass. Among the SF proteins, SRSF7 in GC and SRSF3 in CRC presented the second highest UI and DA, respectively. However, these were less reliable detection markers of both cancers than CEA.

Although the aberrant cytoplasmic localization of HNRNPA1 has already been observed in CRC tissue [10], subcellular distribution analyses of SRSF3 and SRSF 5–7 in primary cancer cells are rare in the literature, and no studies further describe the extranuclear fractions (cytosol/microsome vs. membrane/organelle). We employed a pilot approach using biochemical fractionation and immunoblot assay for the SF distribution analysis in GC and CRC tissue. We detected aberrantly upregulated proteins in the extranuclear fractions in both types of cancer cells and noted that the major distribution patterns of proteins in the extranuclear fractions differ based on the SF type [HNRNPA1 and SRSF6 (cytosol/microsome) vs. SRSF1, 3, 5, and 7 (membrane/organelle)] rather than by cancer type (GC vs. CRC). The cytoplasmic expression of SFs in our IHC analysis, especially prominent in cancer samples, could be explained by two plausible mechanisms: (1) tumor-specific translation control by SFs and (2) tumorigenesis-prone or tumor-associated conditions. Regarding the first mechanism, SFs interact with translation regulators (tumor-promoting or -suppressing factors), which are mainly found in cellular organelles (membrane/organelle fraction) and cytosol (cytosol/microsome fraction) [5, 6]. The second mechanism is associated with the post-translational modification of SFs; rapidly growing tumor cell-related or tumorigenesis-prone stress/status (e.g., metaphase, hypoxia or ultraviolet light) promotes the phosphorylation-related extranuclear localization and/or induction of SF proteins [5, 7, 8].

Because HNRNPA1 and SRSFs have a well-known antagonistic relationship [9], we expected an inverse correlation between HNRNPA1 and SRSF proteins, but all of the SF proteins showed positive correlations with each other. All of the SF mRNAs also showed positive correlations with each other in GC or in CRC. This may be due to collapse of the normal homeostatic relationship among the SFs in cancer cells and/or serial escalation within the SF proteins or within the SF mRNAs. Meanwhile, five SF proteins in GC and three SF proteins in CRC showed expression levels disproportional to those of their respective mRNAs. This could be explained by the specific kinetics of each SF or various post-transcriptional modifications [2830] that affect the SF protein production rate or stability.

The major purpose of this study was to select an effective marker for GC. Through the comparison of proteins and mRNAs, aside from selection of a specific SF, our results indicate that the HNRNPA1 protein itself, rather than its mRNA, is the most reliable tool for GC detection among the SFs tested here.

Conclusions

Apart from indicating that among the six SFs, HNRNPA1, but not HNRNPA1 mRNA, is the most clinically applicable marker for GC, our results offer new insight into the expression profiles and diagnostic efficacies of the six SFs in GC that is required for their future practical implementations. Despite the frequent upregulation and very good to excellent detection accuracy of SRSF3 and HNRNPA1 in CRC, we still recommend the use of CEA, rather than the SFs, especially for high-stage CRC detection.

Abbreviations

CEA, carcinoembryonic antigen; CR, colorectal; CRC, colorectal cancer; CRNM, colorectal normal mucosa; DA, detection accuracy; GC, gastric cancer; GNM, gastric normal mucosa; HNRNPA1, heterogeneous nuclear ribonucleoprotein A1; IHC, immunohistochemistry; NM, normal mucosa; PCR, polymerase chain reaction; SF, splicing factor; SRSF, serine/arginine-rich splicing factor; TBS-T, Tris-buffered saline containing 0.1 % Tween 20, TCGA, The Cancer Genome Atlas; UI, upregulation incidence.

References

  1. Das S, Anczuków O, Akerman M, Krainer AR. Oncogenic splicing factor SRSF1 is a critical transcriptional target of MYC. Cell Rep. 2012;1:110–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. He X, Arslan AD, Pool MD, Ho TT, Darcy KM, Coon JS, et al. Knockdown of splicing factor SRp20 causes apoptosis in ovarian cancer cells and its expression is associated with malignancy of epithelial ovarian cancer. Oncogene. 2011;30:356–65.

    Article  CAS  PubMed  Google Scholar 

  3. Cohen-Eliav M, Golan-Gerstl R, Siegfried Z, Andersen CL, Thorsen K, Ørntoft TF, et al. The splicing factor SRSF6 is amplified and is an oncoprotein in lung and colon cancers. J Pathol. 2013;229:630–9.

    Article  CAS  PubMed  Google Scholar 

  4. Zhu S, Chen Z, Katsha A, Hong J, Belkhiri A, El-Rifai W. Regulation of CD44E by DARPP-32-dependent activation of SRp20 splicing factor in gastric tumorigenesis. Oncogene. 2015. doi:10.1038/onc.2015.250.

    Google Scholar 

  5. Twyffels L, Gueydan C, Kruys V. Shuttling SR proteins: more than splicing factors. FEBS J. 2011;278:3246–55.

    Article  CAS  PubMed  Google Scholar 

  6. Karni R, Hippo Y, Lowe SW, Krainer AR. The splicing-factor oncoprotein SF2/ASF activates mTORC1. Proc Natl Acad Sci U S A. 2008;105:15323–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Sanford JR, Ellis JD, Cazalla D, Cáceres JF. Reversible phosphorylation differentially affects nuclear and cytoplasmic functions of splicing factor 2/alternative splicing factor. Proc Natl Acad Sci U S A. 2005;102:15042–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Yang R, Weber DJ, Carrier F. Post-transcriptional regulation of thioredoxin by the stress inducible heterogenous ribonucleoprotein A18. Nucleic Acids Res. 2006;34:1224–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Zhu J, Mayeda A, Krainer AR. Exon identity established through differential antagonism between exonic splicing silencer-bound hnRNP A1 and enhancer-bound SR proteins. Mol Cell. 2001;8:1351–61.

    Article  CAS  PubMed  Google Scholar 

  10. Hope NR, Murray GI. The expression profile of RNA-binding proteins in primary and metastatic colorectal cancer: relationship of heterogeneous nuclear ribonucleoproteins with prognosis. Hum Pathol. 2011;42:393–402.

    Article  CAS  PubMed  Google Scholar 

  11. Ma YL, Peng JY, Zhang P, Huang L, Liu WJ, Shen TY, et al. Heterogeneous nuclear ribonucleoprotein A1 is identified as a potential biomarker for colorectal cancer based on differential proteomics technology. J Proteome Res. 2009;8:4525–35.

    Article  CAS  PubMed  Google Scholar 

  12. Duffy MJ. Carcinoembryonic antigen as a marker for colorectal cancer: is it clinically useful? Clin Chem. 2001;47:624–30.

    CAS  PubMed  Google Scholar 

  13. Zou L, Zhang H, Du C, Liu X, Zhu S, Zhang W, et al. Correlation of SRSF1 and PRMT1 expression with clinical status of pediatric acute lymphoblastic leukemia. J Hematol Oncol. 2012;5:42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. David CJ, Chen M, Assanah M, Canoll P, Manley JL. HnRNP proteins controlled by c-Myc deregulate pyruvate kinase mRNA splicing in cancer. Nature. 2010;463:364–8.

    Article  CAS  PubMed  Google Scholar 

  15. Garcia M, Seigner C, Bastid C, Choux R, Payan MJ, Reggio H. Carcinoembryonic antigen has a different molecular weight in normal colon and in cancer cells due to N-glycosylation differences. Cancer Res. 1991;51:5679–86.

    CAS  PubMed  Google Scholar 

  16. Brown AJ, Sun L, Feramisco JD, Brown MS, Goldstein JL. Cholesterol addition to ER membranes alters conformation of SCAP, the SREBP escort protein that regulates cholesterol metabolism. Mol Cell. 2002;10:237–45.

    Article  CAS  PubMed  Google Scholar 

  17. Wu CC, MacCoss MJ, Howell KE, Yates 3rd JR. A method for the comprehensive proteomic analysis of membrane proteins. Nat Biotechnol. 2003;21:532–8.

    Article  CAS  PubMed  Google Scholar 

  18. Wei Z, Mousseau DD, Richardson JS, Dyck LE, Li XM. Atypical antipsychotics attenuate neurotoxicity of beta-amyloid(25–35) by modulating Bax and Bcl-X(l/s) expression and localization. J Neurosci Res. 2003;74:942.

    Article  CAS  PubMed  Google Scholar 

  19. Wang J, Su L, Chen X, Li P, Cai Q, Yu B, et al. MALAT1 promotes cell proliferation in gastric cancer by recruiting SF2/ASF. Biomed Pharmacother. 2014;68:557–64.

    Article  CAS  PubMed  Google Scholar 

  20. Oliveira C, Pinto M, Duval A, Brennetot C, Domingo E, Espín E, et al. BRAF mutations characterize colon but not gastric cancer with mismatch repair deficiency. Oncogene. 2003;22:9192–6.

    Article  CAS  PubMed  Google Scholar 

  21. Tajima K, Tominaga S. Dietary habits and gastro-intestinal cancers: a comparative case–control study of stomach and large intestinal cancers in Nagoya. Japan Jpn J Cancer Res. 1985;76:705–16.

    CAS  PubMed  Google Scholar 

  22. Nagini S. Carcinoma of the stomach: A review of epidemiology, pathogenesis, molecular genetics and chemoprevention. World J Gastrointest Oncol. 2012;4:156–69.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Russo A, Bazan V, Migliavacca M, Zanna I, Tubiolo C, Tumminello FM, et al. Prognostic significance of DNA ploidy, S-phase fraction, and tissue levels of aspartic, cysteine, and serine proteases in operable gastric carcinoma. Clin Cancer Res. 2000;6:178–84.

    CAS  PubMed  Google Scholar 

  24. Setälä LP, Kosma VM, Marin S, Lipponen PK, Eskelinen MJ, Syrjänen KJ, et al. Prognostic factors in gastric cancer: the value of vascular invasion, mitotic rate and lymphoplasmacytic infiltration. Br J Cancer. 1996;74:766–72.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Ushigome M, Ubagai T, Fukuda H, Tsuchiya N, Sugimura T, Takatsuka J, et al. Up-regulation of hnRNP A1 gene in sporadic human colorectal cancers. Int J Oncol. 2005;26:635–40.

    CAS  PubMed  Google Scholar 

  26. The Cancer Genome Atlas. National Cancer Institute/National Human Genome Research Institute, Bethesda. 2016. https://tcga-data.nci.nih.gov/tcga/ Accessed 15 Apr 2016.

  27. Oncomine database. Compendia Bioscience, Ann Arbor. 2016. https://www.oncomine.org Accessed 15 Apr 2016.

  28. Naro C, Sette C. Phosphorylation-mediated regulation of alternative splicing in cancer. Int J Cell Biol. 2013;2013:151839.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Wu H, Sun S, Tu K, Gao Y, Xie B, Krainer AR, et al. A splicing-independent function of SF2/ASF in microRNA processing. Mol Cell. 2010;38:67–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Sun S, Zhang Z, Sinha R, Karni R, Krainer AR. SF2/ASF autoregulation involves multiple layers of post-transcriptional and translational control. Nat Struct Mol Biol. 2010;17:306–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The biospecimens and clinical data used in this study were provided by the Biobank of Wonkwang University Hospital, Korea Biobank Network.

This research was supported by the Basic Science Research Program through the National Research Foundation of Korea funded by the Ministry of Education (2012R1A1A2041914 and 2014R1A1A2054908) and ICT & Future Planning (2011–0030130 and NRF-2015M3A9E3051054) and Korea Health Industry Development Institute, Ministry of Health and Welfare (HI12C0110).

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article.

Authors’ contributions

WCP and DSP conceived and designed the study. Experiments were performed by WCP, HRK, DBK, JSR, KHC, GOL, KJY, KYK, RP, KHY, JHC, SCC, and DSP. Data were analyzed by WCP, YJL, MCP and DSP. The first two authors (WCP and HRK) contributed equally to this work. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

“Not applicable” in this section.

Ethics approval and consent to participate

The study was approved by the Wonkwang university hospital Institutional Review Board (WKUHIRB2014-HRE-021). All participants gave written informed consent in accordance with the Declaration of Helsinki.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Do-Sim Park.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Park, W.C., Kim, HR., Kang, D.B. et al. Comparative expression patterns and diagnostic efficacies of SR splicing factors and HNRNPA1 in gastric and colorectal cancer. BMC Cancer 16, 358 (2016). https://0-doi-org.brum.beds.ac.uk/10.1186/s12885-016-2387-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12885-016-2387-x

Keywords