RBM15B is upregulated in HCC and predicts poor prognosis.
Firstly, we evaluated the mRNA levels of RBM15B in human HCC using The Cancer Genome Atlas (TCGA) database. RBM15B mRNA levels were remarkably elevated in tumor tissues compared with those in normal tissues (Fig. 1A). Furthermore, RBM15B mRNA levels increased with increasing tumor grade and lymph node metastasis stage (Fig. 1B and 1C). We examined RBM15B mRNA levels in cohort 1 by reverse transcription-quantitative polymerase chain reaction (RT-qPCR). The paired t-test demonstrated that RBM15B mRNA levels were upregulated in tumors compared with those in para-cancerous tissues (Fig. 1D). Data from TCGA database showed that RBM15B protein levels were also prominently expressed in HCC tumor tissues (Fig. 1E). Western blot analysis of cohort 1 showed that RBM15B was highly expressed in tumor tissues (Fig. 1F). We investigated the relationship between RBM15B levels and clinical features in paired HCC MTA samples of cohort 2 and found that higher RBM15B levels were correlated with decreased tumor differentiation, increased microvascular infiltration, and increased AFP levels (Table 1). Immunohistochemistry (IHC) analysis of RBM15B in cohort 2 showed that tumor tissues had higher RBM15B protein expression than para-cancerous tissues (Fig. 1G). Kaplan–Meier analysis revealed that patients with high levels of RBM15B had shorter overall survival (OS) and disease-free survival (DFS) rates (Fig. 1H and 1I). These findings indicate that RBM15B is highly expressed in HCC and that patients with high RBM15B levels have worse outcomes than those with low RBM15B levels.
Table 1
Clinical characteristics of the 56 patients with HCC in Cohort 2
Features | | RBM15B-low | RBM15B-high | p-value |
Gender | Male | 19 | 2 | 0.696 |
| Female | 29 | 6 | |
Age | < 65 | 14 | 7 | 0.1077 |
| ≥ 65 | 30 | 5 | |
AFP | < 20 µg/L | 10 | 9 | 0.0032** |
| ≥ 20 µg/L | 4 | 28 | |
HbsAg | Negative | 5 | 30 | 0.2977 |
| Positive | 6 | 15 | |
Differentiation | Low | 14 | 13 | 0.0023** |
| High | 2 | 19 | |
Microvascular infiltration | Absent | 17 | 4 | 0.0021** |
| Present | 13 | 22 | |
Recurrence | No | 10 | 11 | 0.1411 |
| Yes | 8 | 23 | |
IHC scores were independently evaluated by two proficient pathologists. A semi-quantitative method was used to assess the staining intensity and percentage of positive cells. The staining intensity score was defined as follows:0 (no staining, negative), 1 (light yellow, weak), 2 (dark yellow, medium), and 3 (brown, yellow, strong). The percentage of positive cells was grouped as follows:0 (< 5%), 1 (5–25%), 2 (26–50%), 3 (51–75%), and 4 (> 75%). The intensity score was multiplied by the positive cell score to yield the total score:0–7 was defined as low expression and ≥ 8 was regarded as high expression. Some cases were excluded from statistical analysis due to missing clinical data. A X2 test was used to test the differences between the two variables. Statistical significance: **p < 0.01.
RBM15B strengthens tumor migration and invasion ability in vitro.
RBM15B was depleted and confirmed by RT-qPCR and western blot in MHCC97H and HCCLM3 cells, respectively. Both siRBM15B#1 and siRBM15B#2 showed excellent silencing efficiencies (Fig. 2A and 2B). We performed a Transwell assay and found that the deficiency of RBM15B restrained both the migration and invasion ability of MHCC97H cells. The histograms in Fig. 2C and 2D showed the number of penetrated cells on the right. Similar results were observed in the RBM15B-knockdown HCCLM3 cells (Fig. 2E and 2F). We performed a gain-of-function experiment and confirmed the overexpression of RBM15B using RT-qPCR and western blotting (Fig. 2G and 2H). Overexpression of RBM15B significantly enhanced the mobility and invasion capacity of MHCC97H and HCCLM3 cells (Fig. 2I and 2J). In summary, RBM15B promotes the invading abilities of HCC cells, and lack of RBM15B impairs the invasive abilities of HCC cells.
RBM15B facilitates HCC growth.
To verify the proliferative phenotype of RBM15B, cell counting kit (CCK)-8 and colony formation assays were performed. The CCK-8 assay demonstrated that depletion of RBM15B inhibited the proliferation of MHCC97H and HCCLM3 cells (Fig. 3A and 3B). A colony formation assay showed that the depletion of RBM15B inhibited the colony-forming ability of HCC cells (Fig. 3C and 3D). To further clarify the role of RBM15B in cell proliferation, gain-of-function experiments were performed. CCK-8 assay showed that overexpression of RBM15B strongly accelerated cell propagation (Fig. 3E and 3F). Moreover, the colony-forming assay showed that upregulated RBM15B levels enhanced the colony-forming ability of both MHCC97H and HCCLM3 cells (Fig. 3G and 3H). To verify these effects in vivo, lentivirus (LV)-shRBM15B or LV-shNC was transfected into MHCC97H and HCCLM3 cells. Stable shRBM15B knockdown efficiency in MHCC97H and HCCLM3 cells was confirmed by western blot analysis (Fig.S2A), which were subcutaneously implanted into 4-week-old male nude mice. The LV-shRBM15B group yielded smaller tumors than the LV-shNC group (Fig. 3I-L). A remarkable decrease in both tumor volume (1 / 2×[length × width2]) and weight was observed in the LV-shRBM15B group compared with that in the shNC group (Fig. 3M-P). IHC of RBM15B and Ki67 was performed for MHCC97H subcutaneous tumors. The RBM15B expression in the shRBM15B group was lower than that in the shNC group (Fig. S2B). Ki67 expression was lower in the shRBM15B group than that in the shNC group (Fig. 3Q). In conclusion, RBM15B facilitates cancer propagation and depletion of RBM15B greatly inhibits tumor growth in vitro and in vivo.
ITSN2 is a potential downstream target of RBM15B in HCC.
To determine the exact molecular mechanism underlying the tumor-promoting phenotype of RBM15B, methylated (m6A) RNA immunoprecipitation sequencing (meRIP-seq) and RNA-seq were performed to evaluate the m6A peaks and transcriptional alterations between shNC and shRBM15B in MHCC97H cells. Differential m6A peaks were determined at p < 0.05, and log fold change (FC) ≤ − 1. Differential gene expression were determined at p < 0.05, and log|FC| ≥ 1. To screen transcripts with significantly downregulated m6A peaks, we determined the intersection of meRIP-seq and RNA-seq data in a Venn plot, as indicated by the hypo-down and hypo-up sections in the density map (Fig. 4A). Consequently, 38 transcripts were identified in the present study. Transcripts with FPKM < 0.05 were excluded due to extremely low expression. Among these, there were two novel transcripts, one lincRNA, one antisense RNAs, and one pseudogene. We focused on the analysis of the mRNA transcripts (Fig. 4B). Theoretically, RBM15B can be combined with target transcripts to regulate m6A peaks and alter transcript expression. First, the possible candidates were verified by RT-qPCR. Among the 24 potential candidates, ITSN2 was the only unique alternative whose levels declined consistently upon depletion of RBM15B with two separate ITSN2 siRNAs in MHCC97H and HCCLM3 cells (Fig. 4C-F). Furthermore, RBM15B deficiency in MHCC97H and HCCLM3 cells induced a uniform decrease in the ITSN2 protein levels (Fig. 4G and 4H). Accordingly, we hypothesize that ITSN2 is a promising candidate target for RBM15B.
RBM15B promotes ITSN2 mRNA stabilization via IGF2BP1.
As RBM15B is a member of the MTC family, we hypothesized that RBM15B acts on the downstream target ITSN2 in an m6A-dependant manner. We determined whether depletion of RBM15B would downregulate global m6A levels in HCC by m6A dot blots. The results showed that global m6A peaks were reduced upon inhibition of RBM15B in both MHCC97H and HCCLM3 cells (Fig. 5A and 5B). Potential m6A sites of ITSN2 were estimated using RBM v2.0. As expected, the predicted m6A sites were predominantly distributed at the 3' UTR and CDS near the termination codon, as shown in Supplementary File 2. Primers targeting these predicted m6A sites were subjected to methylated RNA immunoprecipitation (meRIP) analysis, followed by RT-qPCR. Fragments of ITSN2 were strongly reinforced in the m6A-IP groups compared to those in the IgG-IP groups. More importantly, compared to those in the shNC group, fewer fragments of m6A-modified ITSN2 transcripts were observed in the shRBM15B group, as determined by an assay with primers assigned to the four predicted m6A sites located at the 3′ UTR and end of the CDS (Fig. 5C and 5D). To further clarify the essential mechanism of m6A, we produced ITSN2 luciferase plasmids, including 288 bps of 3' UTR and 126 bps of the CDS near the termination codon, or mutant (Mut) one, in which the four m6A sites were mutated. For the mutated one, the four m6A sites verified by meRIP were converted into cytosine (C) to eliminate the impact of m6A, whereas these four m6A sites were preserved in the wild type (WT) plasmid (Supplementary File 3). As expected, the luciferase activity of cells transfected with the ITSN2-WT plasmid was attenuated upon the elimination of RBM15B, whereas that of the ITSN2-Mut plasmid appeared to be insusceptible (Fig. 5E and 5F). The effect of m6A writers depends on their readers; therefore, we determined specific reader(s) for RBM15B. As common m6A readers, the IGF2BP family IGF2BP1, IGF2BP2, and IGF2BP3 are involved in the regulation of mRNA stabilization. Hence, we depleted IGF2BP1, IGF2BP2, and IGF2BP3 levels with the respective siRNAs and found that ITSN2 levels were distinctly decreased upon silencing of IGF2BP1 both in MHCC97H and HCCLM3 cells (Fig. 5G and 5H), but not with silencing of IGF2BP2 or IGF2BP3 (Fig. 5I-L). Western blotting analysis also indicated that ITSN2 levels decreased following IGF2BP1 depletion (Fig. 5M). RNA-binding protein immunoprecipitation (RIP) and RT-qPCR were performed to confirm the binding affinity of IGF2BP1 to ITSN2 mRNA. ITSN2 mRNA levels were evidently enriched by IGF2BP1 antibody compared to IgG antibody treatment (Fig. 5N and 5O). An RNA decay assay showed that the degradation of ITSN2 mRNA was strengthened upon reduction of RBM15B levels, thus reducing the levels of ITSN2 mRNA (Fig. 5P and 5Q). Therefore, we conclude that RBM15B facilitates the overall modification of m6A levels in HCC, especially at the 3′ UTR and CDS near the stop codon of ITSN2 mRNA. These results confirm the importance of the RBM15B-IGF2BP1-ITSN2 regulation axis in HCC; specifically, RBM15B promotes the installation of m6A at specific sites of ITSN2 mRNA and, upon transfer to the cytoplasm, IGF2BP1 recognizes the modified m6A and subsequently binds to it to protect ITSN2 mRNA from degradation.
ITSN2 plays a vital role in promoting HCC.
ITSN2 mRNA was prominently fortified in HCC tissues compared to that in normal liver tissues, according to TCGA (Fig. 6A). Moreover, ITSN2 mRNA levels increased with increasing tumor grade (Fig. 6B) and lymph node stage (Fig. 6C). We evaluated ITSN2 mRNA and protein levels in cohort 1 and found that tumor tissues had higher ITSN2 mRNA levels than paired para-cancerous tissues (Fig. 6D). Western blotting analysis also indicated that tumor tissues had enhanced protein expression of ITSN2 compared with para-cancerous tissues (Fig. 6E). Kaplan–Meier analysis revealed that individuals with higher expression of ITSN2 had worse OS (p = 0.046) than those with lower expression of ITSN2 (Fig. 6F). We performed TMA in cohort 2 to determine ITSN2 expression via IHC and found that tumor tissues had higher expression of ITSN2 in the cytoplasm (Fig. 6G). We explored the ITSN2 phenotype in HCC in vitro. CCK-8 and colony formation assays showed that exhaustion of ITSN2 levels inhibited the proliferation of MHCC97H and HCCLM3 cells (Fig. 6H-L). Transwell assays revealed that inhibition of ITSN2 expression decreased the migration and invasion abilities of MHCC97H (Fig. 6M and 6N) and HCCLM3 cells (Fig. 6O and 6P). Therefore, ITSN2 expedites HCC proliferation and invasion, and patients with higher levels of ITSN2 are predicted to have a worse clinical prognosis.
Depletion of ITSN2 abolishes tumor-promoting phenotype of RBM15B in HCC.
Rescue experiments were performed by overexpressing RBM15B, followed by ITSN2 depletion. Western blot analysis verified the overexpression of RBM15B and the silencing of ITSN2 (Fig. S1A). CCK-8 and colony formation assays showed that ITSN2 silencing impaired the propagation-promoting effects of RBM15B overexpression (Fig. S1B and S1C), and ITSN2 exhaustion reversed the enhanced colony-forming ability associated with RBM15B overexpression (Fig. S1D, and S1E, respectively). Transwell assays indicated that ITSN2 silencing reversed the invasion-promoting effect of RBM15B upregulation (Fig. S1F-I). These findings suggest that enhanced ITSN2 levels may contribute to HCC development through the RBM15B-ITSN2 regulatory axis. In summary, ITSN2 depletion reverses the tumor-promoting effects of RBM15B upregulation.
Combined high expression of RBM15B and ITSN2 predicts an inferior prognosis in HCC.
To assess the relationship between RBM15B and ITSN2, IHC for RBM15B and ITSN2 was performed using the TMA of cohort 2. IHC showed that RBM15B expression was positively correlated with ITSN2 expression (Fig. S2C) as show in the histogram in Fig. S2D, showing high and low expression ratios of RBM15B and ITSN2, respectively. We observed a positive correlation between RBM15B and ITSN2 expression in the TCGA database (Fig. S2E). In summary, RBM15B expression is positively correlated with ITSN2 expression, and high co-expression of these two molecules predicts an inferior outcome in patients with HCC.