Until recently, post-transcriptional gene regulation (PTGR), in contrast to transcriptional regulation, was not extensively explored in cancer, even though it seems to be highly important. PUM proteins are well described in the PTGR of several organisms and contain the PUF RNA-binding domain that recognizes the UGUANAUA motif, located mostly in the 3′ untranslated region (3′UTR) of target mRNAs. Depending on the protein cofactors recruited by PUM proteins, target mRNAs are directed towards translation, repression, activation, degradation, or specific localization. Abnormal profiles of PUM expression have been shown in several types of cancer, in some of them being different for PUM1 and PUM2.
genes are ubiquitously expressed in human tissues. Taking advantage of the Cancer Genome Atlas database, we analyzed RNA expression levels of
and
from human cancer samples and compared them to
expression in healthy tissues collected within the Genotype-Tissue Expression (GTEx) project. The accumulated data show that the expression levels of
and
are significantly altered in 17 types of cancer tissues (
1). Notably, in the majority of cancer tissues, PUM2 RNA expression was considerably higher than PUM1 expression. Interestingly, in almost all of the samples, the PUM1 level was increased compared to healthy tissues (except adrenal gland and bladder cancers). In the case of PUM2, it was overexpressed in almost all samples, except ovarian and uterus cancer tissues, where the RNA expression level was lower than in healthy tissues. These differences in expression patterns between PUM1 and PUM2 may indicate that they have divergent functions in adrenal gland, bladder, ovarian, and uterus cancers. Further studies are required for the proper interpretation of PUM1 and PUM2 expression patterns. In particular, alterations of RNA expression of PUMs require validation at the protein level.
Comparison of PUM1 and PUM2 RNA expression levels in 17 types of cancers and corresponding healthy tissues. Data were obtained for cancers from The Cancer Genome Atlas (October 2020) (PUM1/2 cancer) and for healthy tissues from the Genotype-Tissue Expression project (October 2020) (PUM1/2 tissue). Statistical significance was calculated using
-tests. Violin-box-whisker plots were generated using ggplot2 library in R. NS—non-significant, TPM—transcript per million (reads), *
< 0.05, **
< 0.01, ***
< 0.005.
To gain comprehensive insights into the importance of PUM-controlled PTGR in cancer as well as into potential functional differences between PUM proteins, the direct identification of PUM-bound RNAs by protein-RNA co-immunoprecipitation (co-IP) followed by RNA-Seq analysis was performed in TCam-2 cells. This cell line represents seminoma, exemplifying a type of testis germ cell tumor[1]. Considering that RNA-PUM binding does not imply regulation per se, transcriptome-wide identification of quantitative changes of mRNAs upon knock-down of PUM proteins was conducted in parallel. In this study, only mRNAs fulfilling two criteria were considered as regulated by PUM1 or PUM2 proteins, namely (1) binding to PUM, and (2) changing levels, as detected by differential gene expression analysis upon PUM knockdown. In total, 346 PUM1-regulated and 141 PUM2-regulated targets were identified in that study. About 90% of PUM-regulated targets were different for PUM1 and PUM2, and nearly 100% of all identified targets contained PBEs, thus validating the results[2]. Several other studies aimed at the global identification of human PUM1 and PUM2-regulated targets were performed in other human cell lines and using some methodological modifications compared to the study in TCam-2 cells. They all showed that the number of non-overlapping PUM1 and PUM2 targets was significant, ranging from 57% to 95%[3][4][5][6].
The high numbers of PUM1 and PUM2-specific targets that have been identified in those studies is puzzling, considering that despite the fact that PUF1 domain is less curved than PUF2 (
), these proteins recognize the same PBE motif. Notably, in contrast to the highly conserved C-terminal, the N-terminal regions of PUM1 and PUM2 exhibit a much lower degree of evolutionary conservation and contain low-complexity regions with insertions and deletions[7]. It was demonstrated that the N-terminal region of human PUM2 serves as an interface for interaction with other proteins, such as DDX20 (GEMIN3)[8]. Additionally, based on studies of the PUM cofactors—the NANOS family—it was demonstrated that only NANOS3 protein, and not NANOS1 or NANOS2, cooperates with PUM2 in repression of mRNA encoding the SIAH1 tumor suppressor[9]. This study showed that PUM may have different mRNA targets in cooperation with different combinations of cofactors. Therefore, it was proposed that PUM1 and PUM2 regulate separate mRNA pools by co-operating with different proteins. A global immunoprecipitation (IP) and mass spectrometry (MS)-based identification of PUM1- and PUM2-binding proteins showed that PUM1 and PUM2 indeed interact mainly with different groups of proteins, a majority of them representing RBPs. Combinatorial analysis of RNA immunoprecipitation (RIP) and RNA-Seq MS, as well as motif enrichment analysis for RNA-binding proteins in PUM1 and PUM2 mRNA targets, revealed that PUM1 and PUM2 form separate ribonucleoprotein networks[2], which resemble so-called “regulons” or “posttranscriptional operons”, as previously suggested for RBPs that recognize specific motifs in RNAs[10]. According to those networks, PUM1 and PUM2 may cooperate with varied protein cofactors to regulate mRNA target sub-pools responsible for unique pathways. For example, PUM1 with IGF2BP3 cofactor may co-regulate mRNAs implicated in regulation of cell division, while PUM1 with PABPC4 and MATR3 may co-regulate mRNAs involved in the epidermal growth factor receptor-signaling pathway. On the other hand, PUM2 with G3BP2, HNRNPA1, FXR2, and SFPQ co-factors may co-regulate mRNAs involved in the regulation of Rho protein signal transduction, while with MATR3, PUM2 may affect mRNAs that negatively regulate cell development[2]; however, such scenarios require validation at the protein level. One important aspect that should be taken into account when discussing the functional relationship between PUM1 and PUM2 interactomes is that PUM1 and PUM2 mRNAs contain functional PBEs in their 3′UTRs. Consequently, as has already been shown at the RNA and protein levels, they repress each other[11][12].
In addition to the regulation of stem cell fate, growth, and development, PUMs are involved in the control of a myriad of other biological functions such as hematopoiesis, neurogenesis, and gametogenesis, while their dysfunction contributes to several diseases, which is in line with their ubiquitous pattern of expression, as they have been identified in many human tissues.
Notably, there is an emerging field of research concerning the involvement of PUM proteins in specific types of cancer. It has also been demonstrated that PUM proteins control the levels of several mRNAs that encode proteins involved in processes that are often disrupted in cancerogenesis, such as apoptosis, proliferation, and the cell cycle.
Hematopoiesis is the process of differentiation of hematopoietic stem cells (HSCs) residing in the bone marrow, giving rise to various types of blood cell. Leukemia occurs when cancer-driving mutations cause conversion of some HSCs to leukemia stem cells (LSCs), which then abnormally proliferate and grow in size, resulting in bone marrow malfunction. New evidence suggests that PUM1 and PUM2 play vital regulatory roles in the maintenance and proliferation of normal human and mouse stem cells (SCs)[13][14]. It has also been demonstrated that both PUMs are overexpressed in a majority of acute myeloid leukemia (AML) samples, as well as in the cell-lines derived from those pathogenic samples[13]. Importantly, PUM1 and PUM2 influence the cell cycle, proliferation, and apoptosis of normal human and mouse HSCs, as well as AML cells[13]. PUMs induce these effects by activating the expression of Forkhead box protein P1 (FOXP1) transcription factor mediated by direct binding PBE motifs located in the 3′UTR of FOXP1 mRNA. This PUM-mediated FOXP1 activation suppresses the expression of cell cycle inhibitors, such as CDKN1B, thereby promoting proliferation. This is in line with the findings that FOXP1 itself regulates cell proliferation and differentiation and is crucial for the regulation of hematopoietic stem progenitor cells during leukemic cell growth[13]. Importantly, PUM1/2 can also directly suppress the expression of a tumor suppressor, CDKN1B[15][16][12], potentially enhancing the FOXP1 pathway in HSCs and AML cells. Additionally, the expression levels of multiple other proteins involved in HSCs and LSCs are likely to be influenced by PUM1/2 3′UTR-mediated regulation, which is a topic for future study[13].
Seminoma is the most common type of testis germ cell tumor among western populations and is increasing in frequency, with the highest incidence in men at reproductive age [1]. In the TCam-2 seminoma cell line, it has recently been demonstrated that PUM1 and PUM2 regulate several mRNAs functionally linked to cancer. Among such PUM1 and PUM2-repressed targets, mRNA encoding SPINDLIN1 (SPIN1) is known for its role in mammalian gametogenesis. This spindle-binding protein is necessary for the meiotic progression of germ cells, while abnormal overexpression of SPIN1 is associated with human ovarian cancers and has been observed in some other cancer cell lines[17][18][19][20]. SPIN1 protein is expressed in TCam-2 cells and its overexpression causes a significant increase in proliferation and a reduced level of apoptosis. Those two features indicate that SPIN1 plays a proto-oncogenic role in TCam-2 cells. PUM1 and PUM2 repress a SPIN1 homologue called SPIN3, whereas its overexpression elicits a decrease in proliferation and an increase in apoptosis of TCam-2 cells[11]. Taken together, SPIN1 demonstrated proto-oncogenic properties, while SPIN3 showed tumor suppressor features. PUM1 itself, but not PUM2, strongly stimulated apoptosis and moderately slowed down cell cycle progression, suggesting that PUM1, similarly to SPIN3, plays the role of a tumor suppressor in TCam-2 cells. Altogether, by acting as SPIN1 and SPIN3 repressors, PUM proteins may promote a normal human male germ cell apoptotic status and thus prevent cancer [11].
In TCam-2 cells, PUMs also cause repression of mRNA encoding kinesin KIF18A[21]. This kinesin exists on positively charged ends of microtubules in the vicinity of the kinetochore, and it regulates the dynamics of mitotic cell division. The knockout of Kif18a in a mouse causes mitotic arrest and apoptosis of the male germ cells, resulting in infertility. This kinesin is an important regulator of the cell cycle and apoptosis in human germ cells as well[22][23]. Notably, it also plays the role of a proto-oncogene, as it is overexpressed in many cancer types. It has been reported that KIF18A positively influences TCam-2 cell proliferation, downregulates apoptosis, and promotes cell cycle progression, these effects being opposite to the effects of PUMs. Therefore, repression by PUM proteins may represent one of the mechanisms affecting KIF18A levels in regulating proliferation, the cell cycle, and apoptosis in TCam-2 cells[21].
The findings concerning PUM1 and PUM2 mRNA targets and their functional relations in TCam-2 cells are of interest and should be validated in patients suffering from testis germ cell tumors and several types of cancer affecting other human tissues.
The functional relationship between 3′UTR-binding RBPs and non-coding RNAs is an emerging topic in PTGR research. Until recently, this issue had not been sufficiently explored due to methodological limitations. Several studies on the mechanisms of the interplay between PUM and non-coding RNAs have been published in recent years. These studies provide valuable data on the interplay between PUM and non-coding RNA in the context of cancer.
3.1. PUM Proteins Cooperate with microRNA Machinery in CancerFirst, in global screens of HeLa cells for PUM mRNA targets, it was documented that microRNA complementary sites were enriched in the vicinity of PBE motifs within PUM target mRNAs[3]. The significance of that vicinity has partly been clarified in studies describing the repression of mRNA encoding cyclin-dependent kinase inhibitor 1B (CDKN1B) by PUM1/2 proteins[12]. CDKN1B downregulates cell cycle progression by blocking CDK2 activity. It plays the role of a tumor suppressor, its level being lowered in several types of cancer, while becoming increased in non-dividing cells[24]. CDKN1B mRNA is a target for miR-221- and miR-222-mediated repression, and high levels of these two miRNAs correlate inversely with low CDKN1B levels in cancer, for example in glioblastoma[25]. It was shown in human fibroblasts that PUM binding to 3′UTR of CDKNB1 mRNA induces a local conformational rearrangement, making the miRNA complementary site accessible for miRNA hybridization. Through that mechanism, PUM1/2 enables the miRNA-induced repression of CDKN1B mRNA translation[12].Cooperation between the repression by PUM proteins and miRNAs-induced repression was also demonstrated in the regulation of E2F3 transcription activator, a well-described human oncogene. E2F3 belongs to a family of E2F transcription factors, which regulate cell proliferation as well as apoptosis, while their dysregulation frequently accompanies human malignancies [26]. The 3′UTR of mRNA encoding E2F3 contains two PBE motifs, which have been shown to be active in PUM-induced repression, as well as several complementary sites for miRNAs. It has been demonstrated that PUM1/2 proteins downregulate E2F3 expression cooperatively with mir-503 and mir-125pp. Interestingly, several cancer cell lines escape the PUM-mediated regulation by shortening the E2F3 3′UTR, eliminating PBEs, and thus disabling PUM-mediated repression[27]. In some cancer types, e.g., in bladder carcinoma, mir-503, and mir-125p, levels are significantly decreased, indicating their importance and their anti-oncogenic role[28].
Long non-coding RNAs (lncRNAs) are molecules longer than 200 nucleotides and lack an open reading frame. lncRNAs often contain microRNA complementary sites and act as miRNA sponges, thereby regulating their levels inside the cell. This feature of lncRNAs may result in upregulation of several oncogenic miRNA targets, thus contributing to carcinogenesis[29][30]. Interestingly, two simultaneous studies showed that a lncRNA, known as non-coding RNA activated by DNA damage (NORAD), contains a number of motifs that correspond to several RBPs, including 15 PUM-corresponding PBEs [31][5], suggesting a sponge-activity-based PUM sequestration. Given the diverse roles of PUM proteins in regulating proliferation, the cell cycle, and apoptosis-related mRNAs, as discussed in the previous sections, the regulation of PUM activity inside the cells by NORAD-mediated sequestration exerts a significant impact on carcinogenesis. NORAD is a ubiquitous, abundant, and conserved lncRNA in mammals. Both PUM1 and PUM2 bind to NORAD lncRNA. Compared to PUM1, PUM2 has a higher affinity for NORAD. The identification of PUM2 mRNA targets performed upon NORAD knockdown demonstrated the activity of the NORAD–PUM interaction and its importance for fine-tuning the levels of mRNAs encoding proteins related to DNA repair and the cell cycle [5]. Importantly, inactivation of NORAD in the HTC116 cell line leads to chromosomal instability by increasing chromosomal and mitotic abnormalities. These studies brought attention to the role of the NORAD–PUM interaction in regulating genomic stability. Mass spectrometry analysis identified over 500 RBPs that bind to NORAD fragments in the U2OS cell line. Interestingly, one of the identified RBPs is SAM68, a protein required for PUM activity regulation by NORAD. It was additionally demonstrated that SAM68 and PUM2 interact independently of NORAD, an interaction that is required for PUM2 binding to NORAD [32].Further research aimed at elucidating the exact molecular mechanism of the control of genome stability by the NORAD–PUM axis showed opposing findings. Previous findings showed that NORAD is predominantly located in the cytoplasm and acts as a PUM sponge, but single-molecule RNA fluorescent in situ hybridization (smRNA FISH) demonstrated that NORAD is also located in the nucleus. In addition, upon stress, NORAD localization in the nucleus substantially increases. Finally, NORAD was shown to regulate genomic stability by modulating the availability of the RNA Binding Motif Protein X-Linked (RBMX) protein inside the nucleus for topoisomerase complex formation, and PUM was shown to be dispensable for NORAD function[33]. A follow-up study addressed these findings and showed that NORAD indeed acts as a PUM sponge inside the cytoplasm, independently of stress, and regulates genomic stability and proper mitotic divisions, although PUM sequestration and RBMX are dispensable for this function[34]. The question of physiological significance of the NORAD–PUM interaction was addressed by using a mouse model. The research demonstrated that removal of mouse NORAD lncRNA, having 61% sequence identity to human NORAD, resulted in a phenotype that resembles premature aging due to genomic instability and mitochondrial dysfunction[34]. PUM2 overexpression resulted in a phenotype similar to NORAD-deficient mice.Although there are opposing results regarding NORAD localization, its role in regulating RBMX protein availability for topoisomerase complex formation, and whether this is required for NORAD function in genomic stability, the clear evidence presented by multiple studies shows that NORAD is a conserved lncRNA that regulates PUM proteins in mammals[35] [
60]. A follow-up study addressed these findings and showed that NORAD indeed acts as a PUM sponge inside the cytoplasm, independently of stress, and regulates genomic stability and proper mitotic divisions, although PUM sequestration and RBMX are dispensable for this function [61]. The question of physiological significance of the NORAD–PUM interaction was addressed by using a mouse model. The research demonstrated that removal of mouse NORAD lncRNA, having 61% sequence identity to human NORAD, resulted in a phenotype that resembles premature aging due to genomic instability and mitochondrial dysfunction [62]. PUM2 overexpression resulted in a phenotype similar to NORAD-deficient mice.]. This regulation is critical for genomic stability and proper mitotic division. These studies served as a proof that lncRNAs function as RBP sponges inside the cells and, most importantly, they demonstrated the critical role of the PUM sequester by NORAD in carcinogenesis.