2.1 PiR-26681 is decreased in ovarian cancer and correlated with primary stage and positive patient outcomes
We analyzed a pair of human ovarian cancer samples, comparing primary-stage samples with matched advanced-stage samples using Pandora sequencing to identify differentially expressed piRNAs. PiR-26681 (piRBase Id: piR-has-26681, NCBI GenBank: DQ596465.1) was selected for further analysis due to its low expression in the advanced stages (fold change value= -3.436), as shown in Fig. 1A. Subsequent verification on a larger dataset, which included 43 patients with primary-stage ovarian cancer, 42 patients with advanced-stage ovarian cancer, and 19 normal ovaries, revealed that piR-26681 expression is significantly lower in ovarian cancer tissues than in normal ovarian tissues (Fig. 1B). Additionally, piR-26681 levels were higher in primary-stage ovarian cancer compared to the advanced stage (Fig. 1C). Kaplan-Meier analysis showed that low levels of piR-26681 in patients with ovarian cancer are associated with poorer progression-free survival (PFS) (Fig. 1D).
Further, qRT-PCR studies confirmed that piR-26681 expression was lower in ovarian cancer cell lines compared to normal ovarian cell lines (IOSE-80), with the lowest expression observed in CAOV3 cells and the highest in OVCAR3 cells (Fig. 1E). Consequently, our forthcoming experiments will be conducted in CAOV3 and OVCAR3 cells.
2.2 Overexpression of piR-26681 inhibited the malignant biological behavior of ovarian cancer cells
To examine the physiological impacts of piR-26681 on ovarian cancer, we upregulated its expression in CAOV3 and OVCAR3 cells using the mimic technique and subsequently assessed its biological properties. qRT-PCR was used to validate the transfection efficiency of cells transfected with mimic (Fig. 1F, 1G). A series of assays, including CCK8, plate clone, EdU incorporation assay, apoptosis, wound healing, and transwell assays, were conducted to evaluate the impact of piR-26681 on the proliferation, apoptosis, migration, and invasion of ovarian cancer cells. Overexpressing piR-26681 resulted in the suppression of malignant biological activity in ovarian cancer cells, including proliferative ability as observed in CCK8 assay (Fig. 1H, 1I), clonal formation (Fig. 1J), and the EDU test (Fig. 1K, 1L). This was also accompanied by a reduction in migration (Fig. 1M, 1N) and invasion (Fig. 1O) and an increase in cell apoptosis (Fig. 1P).
2.3 Downregulation of piR-26681 promoted the malignant biological behavior of ovarian cancer cells
To further confirm the function of piR-26681 in ovarian cancer, we down-regulated piR-26681 by transfecting sh-piR-26681 in the CAOVR3 and OVCAR3 cells. qRT-PCR was used to validate the transfection efficiency of cells transfected with sh-piR-26681 (Fig. 2A). The malignant nature of cancer cells was enhanced by interference with piR-26681 expression, as evidenced by their proliferative ability in the CCK8 assay (Figs. 2B), clonal formation (Fig. 2C), and the EDU test (Figs. 2D). Alongside this, there was an increase in invasion (Fig. 2E) and migration (Fig. 2F) and a decrease in cell apoptosis (Fig. 2G). These findings suggest the crucial role of piR-26681 in suppressing the malignant behavior of ovarian cancer cells, suggesting that its manipulation could influence the progression of the diseases.
2.4 PiR-26681 forms a bond with METTL3 and METTL14
To explore the precise mechanisms through which piR-26681 exerts its anti-cancer effects in ovarian cancer, we utilized the RBPsuite website (http://www.csbio.sjtu.edu.cn/bioinf/RBPsuite/) 19 to predict potential binding protein. PiR-26681 is likely to interact with METTL3 (score: 0.997) and METTL14 (score: 0.904). Prior research has shown that many piRNAs are involved in regulating m6A-associated genes in cancer, playing either a pro-cancer or anti-cancer role. Thus, it is crucial to investigate the relationship between pir-26681 and METTL3 and METTL14, enzymes responsible for m6A modification. The interactions of piR-26681 with METTL3 and METTL14 were further validated by RNA immunoprecipitation (RIP) and pulldown assays (Fig. 3A, 3D). Additionally, immunofluorescence tests confirmed the co-localization of piR-26681 with METTL3 and METTL14, providing visual evidence of their interaction (Fig. 3B,3E). These findings highlight the potential regulatory pathways through which piR-26681 may influence the progression of ovarian cancer.
2.5 The upregulation of piR-26681 increases the expression of METTL3 and METTL14
Previous experiments have confirmed that piR-26681 can bind to METTL3 and METTL14. To investigate further the impact of piR-26681 on the biological function of these enzymes, we manipulated piR-26681 levels in CAOV3 and OVCAR3 cells. Specifically, we increased the level of piR-26681 in CAOV3 cells and decreased it in OVCAR3 cells. The result showed that overexpression of piR-26681 in CAOV3 cells led to an up-regulation of METTL3 and METTL14 expression. Conversely, suppressing piR-26681 expression in OVCAR3 cells resulted in a down-regulation of these enzymes. These findings suggest that piR-26681 directly influences the expression levels of METTL3 and METTL14, which may contribute to its effects on ovarian cancer cell behavior (Fig. 3C, 3F).
2.6 Simultaneous knockdown of METTL3 and METTL14 inhibits ovarian cancer development
Previous research indicates that METTL3 acts as a stimulator in ovarian cancer 20, whereas METTL14 has an anti-cancer function in various cancers, including ovarian cancer 21,22. The finding that piR-26681 suppresses ovarian cancer while increasing the expression of both METTL3 and METTL14 presents an apparent contradiction. Yang et al. reported that the simultaneous knockdown of METTL3 and METTL14 increased cell proliferation and invasion 23. We conducted separate and combined knockdowns of METTL3 and METTL14 in CAOV3 and OVCAR3 cells to further investigate these effects. The results revealed that knocking down METTL3 alone inhibited cell growth, whereas suppression of METTL14 alone and concurrent suppression of both METTL3 and METTL14 enhanced cell growth, as evidenced by the CCK-8 test and plate cloning experiment (Fig. 3G,3H). Likewise, as METTL3 was knocked down, cell migration decreased, METTL14 was knocked down, and both genes were knocked down together, increasing cell migration. (Fig. 3I,3K). Additionally, Transwell experiment results demonstrated that inhibiting METTL3 reduced the invasion of ovarian cancer cells. In contrast, knocking down METTL14, as well as both METTL3 and METTL14 together, enhanced invasion (Fig. 3J, 3L). Furthermore, cell apoptosis was observed to increase following METTL3 knockdown, while a decrease was noted following METTL14 knockdown and when both genes were knocked down simultaneously (Fig. 3M, 3N).
These findings highlight the complex role that METTL3 and METTL14 play in ovarian cancer progression and underscore the need for nuanced approaches to targeting these enzymes in cancer therapy.
2.7 PiR-26681 improves the protein stability of METTL3 and METTL14 and facilitates their binding
The catalytic functions of methyltransferases, such as their ability to bind to target RNA and transfer methyl groups, are largely contingent upon the formation of a stable METTL3-METTL14 heterodimer, as demonstrated in previous studies 24,25. METTL3 and METTL14 typically function as heterodimers 24. The co-localization of piR-26681 with METTL3 and METTL14 was observed in an immunofluorescence experiment (Fig. 4A).
To further analyze these interactions, we utilized ZDOCK 3.0.2 software to forecast the binding sites of piR-26681 on the METTL3-METTL14 protein complexes 26( Fig. 4B). The interaction score between piR-26681 and the METTL3-METTL14 heterodimer was significant, at scored as 1071.365. ZDOCK 3.0.2 prediction helped us to determine the binding pattern of piR-26681 with the METTL3-METTL14 complex, showing METTL14 in green and METTL3 in magenta. PiR-26681 is illustrated in yellow, and the yellow dashed line signifies hydrogen bonding. This visualization shows piR-26681 binding between two proteins and interacting with both simultaneously. Hydrogen bonding plays a crucial role in these interactions. For example, residues A-23, A-6, and T-22 on piR-26681 from hydrogen bonds with S299, R298, and P402 on the METTL14 proteins, respectively. Additionally, T-22, A-5, C-8, and G-20 on piR-26681 interact through hydrogen bonding with T-9, C-10, T-19, R471, H400, E403, E481, K459, Q462, T469, and H474 on the METTL3 protein. These interactions highlight the complex molecular interplay between piR-26681 and the METTL3-METTL14 complex, suggesting a potential mechanism by which piR-26881 may influence the methylation landscape in ovarian cancer cells.
The anticipated binding location of amino acid R298 on METTL14 is crucial for the activity of the heterodimer complex it forms with METTL3. Additionally, amino acids Q462, T469, and H474 on METTL13 are involved in the formation of an interfacial loop that interacts with METTL14 27. These structural details provide strong evidence that piR-26681 influences the interaction between METTL3 and METTL14. Our co-immunoprecipitation (Co-IP) test revealed that overexpression of piR-26681 led to an increase in the binding affinity between METTL3 and METTL14. Conversely, knocking down piR-26681 resulted in a decreased interaction between these proteins (Fig. 4C). These findings suggest that piR-26681 plays a significant role in regulating the interaction between METTL3 and METTL14, potentially affecting the function of their heterodimer complex.
Furthermore, the interplay between METTL13 and METTL14 contributes to the stability of each other’s proteins 28. Wang et al. demonstrated that the binding of METTL3 to METTL14 not only enhances the stability of METTL14 protein but also protects it from degradation29. Based on these observations, we speculate that piR-26681 not only promotes the binding between METTL3 and METTL14 but also increases their protein stability. To test this hypothesis, we conducted cycloheximide (CHX) chase experiments and MG132 treatments to assess the effect on protein stability. Our results demonstrated that overexpression of piR-26681 enhances the stability of the METTL3 protein (Fig. 4D, 4E). Similarly, there was a notable increase in the stability of METTL14 protein following piR-26681 overexpressing (Fig. 4G, 4H). Conversely, interference with piR-26681 expression resulted in decreased stability of both METTL3 and METTL14 protein (Fig. 4F, 4I).
2.8 Analysis of the m6A methylome and downstream targets of piR-26681 in ovarian cancer
The immunostaining data showed that the global m6A levels are increased in CAOV3 cells overexpressing piR-26681, a finding that is corroborated by the m6A dot blot assay (Fig. 4J). Conversely, in OVCAR3 cells, m6A levels decreased when piR-26681 expression was inhibited (Fig. 4K). The observation that cells with overexpressed piR-26681 exhibit higher levels of m6A suggests that piR-26681 may improve m6A methylation. This enhancement could play a crucial role in regulating gene expression and various cellular processes, highlighting the potential mechanistic impact of piR-26681 on cellular processes.
To clarify the molecular mechanisms by which piR-26681 influences m6A modifications, we conducted m6A methylation RNA immunoprecipitation sequencing (MeRIP-seq) on both control and piR-26681-overexpressing CAOV3 cells. The analysis of m6A peaks showed that the consensus motif GGAC was notably enriched in the immune-purified RNA (Fig. 3L). Additionally, MeRIP-seq confirmed an increase in mRNA m6A peaks in the piR-26681 overexpression group (Fig. 4M).
Analysis of all differentially methylated genes showed significant enrichment in Gene Ontology (GO) pathways, with the top two pathways being cadherin binding and cell adhesion molecule binding (Fig. 4N). Notably, METTL14 emerged as one of the top 10 genes with increased levels of m6A modification. Among the genes exhibiting high m6A modifications, FBXO16 showed the most significant fold change in differential expression (Fig. 4O). Consequently, METTL14 and FBXO16 were selected as the focal points of our investigation. Specifically, the m6A modification site of METTL14 was identified in the 3’UTR area (Fig. 4P), while the m6A modification sites of FBXO16 were identified in the CDS area (Fig. 4Q).
2.9 PiR-26681 promotes METLL3-METTL14 mediated m6A methylation of METTL14 and FBXO16 mRNA and increases their stability through the IGF2BP2 reader
Upregulation of piR-26681 was found to elevate mRNA levels and improve RNA stability in METTL14 and FBXO16(Fig. 5A,5C). In contrast, blocking piR-26681 reduced both the quantity of RNA and the stability of these genes (Fig. 5B,5D). Moreover, overexpression of piR-26681 led to an increase in FBXO16 protein levels, whereas suppression of piR-26681 resulted in a decrease in FBXO16 protein levels (Fig. 5E,5F).
IGF2BP2, a known m6A reader protein, binds to and stabilizes m6A-modified RNAs 30. Prediction tools from the RBP binding prediction site indicated that IGBF2BP2 would bind to the m6A modification sites identified on METTL14 and FBXO16 1931. Subsequent RIP experiments confirmed that IGF2BP2 significantly enriches the RNA of both METTL14 and FBXO16 (Fig. 5G). Knockdown of IGF2BP2 in CAOV3 cells overexpressing piR-26681 resulted in reduced stability of METTL14 and FBXO16 mRNAs (Fig. 5H), as well as decreased protein level of genes (Fig. 5I). Meanwhile, knocking down METTL3 and METTL14 reduced FBXO16 protein levels in cells overexpressing piR-26681 (Fig. 5J). Additionally, suppressing IGF2BP2 in these cells neglected the cancer-suppressing impact of piR-26681(Fig. 5K,5L,5M, 5N).
2.10 PiR-26681 downregulated Beta-Catenin and suppressed WNT and EMT pathway by increasing FBXO16 expression
Previous experiments have demonstrated that piR-26681 increases the m6A modification of FBXO16 mRNA through the METTL3-METTL14 complex, leading to increased FBXO16 expression. As an E3 ubiquitin ligase, FBXO16 acts as a tumor suppressor in ovarian cancer, with high levels of FBXO16 expression correlating with a better prognosis 32. Additionally, FBXO16 promotes nuclear β-catenin degradation, which suppresses Wnt signaling and inhibits the epithelial-to-mesenchymal transition (EMT) by reducing the amount of β-catenin 3334.
To investigate the effects of piR-26681 on β-catenin, we analyzed protein levels following the overexpression or knockdown of piR-26681. Overexpression of piR-26681 resulted in a decrease in β-catenin levels, mainly in the nucleus, whereas knockdown of piR-26681 resulted in increased β-catenin levels (Fig. 5O). We also assessed the expression of EMT markers in the presence of piR-26681. Cells overexpressing piR-26681 showed increased E-cadherin expression and decreased N-cadherin and vimentin levels. Conversely, the knockdown of piR-26681 led to the opposite effects (Fig. 5P). Further studies into the downstream targets of β-catenin in the Wnt pathway, such as cyclin D1 and C-Myc, revealed that these proteins were decreased in cells overexpressing piR-26681 and increased in which piR-26681 was knocked down (Fig. 5Q).
2.11 Knockdown of FBXO16 or concurrent knockdown of METTL3 and METTL14 rescues the cancer suppression function of piR-26681
In addition, we explored the impact of inhibiting FBXO16 or concurrently knocking down METTL3 and METTL14 in cells that overexpressed piR-26681. We observed that these manipulations suppressed the cancer-inhibitory impacts of piR-26681 (Fig. 6A-H). Specifically, the depletion of FBXO16 or IGF2BP2, or the simultaneous depletion of METTL3 and METTL14, in CAOV3 cells that overexpressed piR-26681 resulted in elevated levels of N-cadherin and vimentin proteins and decreased levels of E-cadherin protein. This indicates an up-regulation of the EMT pathway, which is blocked by piR-26681(Figure S1A). In addition, the knockdown of FBXO16 or IGF2BP2 or depletion of METTL3 and METTL14 also reduced inhibition of the WNT pathways in cells overexpressing piR-26681 (Figure S1B). These results suggest that FBXO16, METTL3, METTL14, and IGF2BP2 are crucial in how piR-26681 stops cancer from growing. In summary, piR-26681 increases the m6a modification of FBXO16 by stabilizing the METTL3-METTL14 methyltransferase complex, which in turn blocks the WNT and EMT pathways.
2.12 PiR-26681 as a promising agent in ovarian cancer
Previous experiments showed the potent cancer-inhibiting effect of piR-26681. For in vivo validation, we upregulated piR-26681 using agopiR-26681 in xenograft models and OC-PDOs. We created xenograft ovarian cancer models using OVCA3 cells to evaluate the impact of piR-26681 on ovarian cancer progression in vivo (Fig. 7A). Treatment with agopiR-26681 resulted in a significant reduction in tumor volume compared to the control group (Fig. 7B, 7C). Western blot and immunohistochemistry analysis of tumor tissue collected from nude mice further showed an increase in METTL3, METTL14, and FBXO16 proteins (Fig. 7D-F).
Consistent with the in vivo tumor model results, treatment with agopiR-26681 in OC-PDOs also reduced organoid diameter compared to the control group (Fig. 7H). These findings indicated that piR-26681 is a promising therapeutic target for ovarian cancer.