miR-33a-3p regulates METTL3-mediated AREG stability and alters EMT to inhibit pancreatic cancer invasion and metastasis

Recent studies have shown that amphoteric regulatory protein (AREG), a member of the epidermal growth factor (EGF) family, is expressed in many cancers and is an independent prognostic indicator for patients with pancreatic cancer, but whether AREG is regulated at the epigenetic level to promote the development of pancreatic cancer (PC) has not been elucidated. Our results support the notion that AREG is overexpressed in pancreatic cancer tissues and cell lines. Functionally, the deletion of AREG impedes pancreatic cancer (PC) cell proliferation, migration, and invasion. In addition, we identified and validated that methyltransferase-like 3 (METTL3) induced the m6A modification on AREG and facilitated the stability of AREG mRNA after sequencing. Additionally, we obtained experimental evidence that miR-33a-3p targets and inhibits METTL3 from taking action, as predicted by using the miRDB and RNAinter. Remediation experiments showed that miR-33a-3p inhibits PC progression through METTL3. In summary, this research reveals that miR-33a-3p inhibits m6A-induced stabilization of AREG by targeting METTL3, which plays a key role in the aggressive progression of PC. AREG could be a potential target for PC treatment.

Patients and clinical samples.Fresh tumor tissues and paired paracancerous tissue specimens were collected from 36 PC patients who underwent surgical resection from June 2021 to October 2022 at the Affiliated Hospital of Jiangnan University.All patients did not receive preoperative chemotherapy or radiotherapy.Tissue specimens were preserved in RNA Keeper Tissue Stabilizer (Vazyme, China).

Stable cell-line construction.
A lentivirus encoding a small hairpin RNA (shRNA) targeting METTL3 was created by Gene Chem (Shanghai, China) (targeting sequence: 5′-A TGC ACA TCC TAC TCT TGT AA-3′).Cells (2 × 10 4 ) were incubated in 24-well plates.After a 24-h culture, cells were infected with lentivirus at an infection multiple of 10.Stably transfected cells were then selected using puromycin (0.5 μg/mL) (MCE, USA).Analysis of the effect was carried out using qPCR.The qPCR primer sequences are listed in Supplementary Table S1.
RNA interference and plasmid transfection.Each well was seeded uniformly at a density of (5 × 10 5 ) cells in a six-well plate and incubated for 24 h to reach a cell density of approximately 70-80%.For RNA interference, a mixture of 50 nM siRNAs and 8 µL RFect Small Nucleic Acid Transfection Reagent (BAIDAI, China) was Protein isolation and western blotting.The protein was extracted by lysing each group of cells with RIPA (ProteinBio, China).Protein concentration was determined with a Bicinchoninic Acid Assay (BCA) protein analysis kit (Vazame, China).Protein samples were isolated on 10% SDS polyacrylamide gels and then electronically transferred to polyvinylidene fluoride (PVDF) membranes (Millipore, USA).The polyvinylidene fluoride membranes were confined and then incubated with anti-METTL3 (Abcam, ab195352, USA), amphiregulin (Proteintech, 16036-1-AP, USA), E-cadherin (Abcam, ab40772, USA), vimentin (Abcam, ab92547, USA), and GAPDH (ProteinBio, PA11002, China) at 4 °C for one night, followed by incubation with anti-IgG secondary antibodies (abclonal, AS014, USA) at room temperature for 1 h.GAPDH was used as an internal control for protein normalization.

Luciferase reporter assay.
To validate the association of miR-33a-3p with METTL3, the predicted miR-33a-3p binding sequence of the 3 ʹUTR region of METTL3 was wild type (WT), and the corresponding 3 ʹUTR domain of mutant METTL3 was MUT.These segments were inserted into the pmirGlo dual luciferase vector to generate separate WT and MUT METTL3 constructs.The METTL3 of pmirGlo WT or MUT was then cotransfected with miR-33a-3p mimics or miR-negative control (NC) in cells by using the lipofectamine 2000.Forty-eight h later, we measured the double luciferase reporter gene assay kit (Beyotime, China), used the sea kidney fluorophore enzyme as an internal reference, and calculated the ratio.

RNA immunoprecipitation (RIP).
RIP was undertaken using the Magna RIP RNA-binding protein immunoprecipitation kit (Millipore, MA, USA) following the manufacturer's instructions.In short, cells were collected and lysed with RIP lysis buffer.Magnetic beads were incubated with the METTL3 antibody at room temperature for 1 h.Immunoglobulin IgG was used as a NC.After washing, the spheres were incubated with cell lysis at 4 °C for 3 h.The spheres were gathered and washed, and the RNA composite was isolated by phenolchloroform extraction.A quantitative polymerase chain reaction was carried out on the RNA-rich segments.
m6A RNA immunoprecipitation (MeRIP) assay.MeRIP was performed using the Magna MeRIP m6A kit (EPIGENTEK, USA) according to the manufacturer's instructions.Additionally, 2 ug of anti-m6A antibody was simply vortexed with protein A/G beads at room temperature for 90 min at 4 °C.The antibody-coupled pellets were washed three times and then incubated in a mixture of proteinase K and protein digestion buffer.The RNA was next eluted, and the RNA that interacts with it was extracted and measured by qPCR.

Statistical analysis.
All data were expressed as mean ± standard deviation (SD).Statistical analysis was carried out using GraphPad Prism 9 software.The t-test (two-tailed) and ANOVA were used to examine differences between two groups or between more than two groups.The correlation between AREG and clinicopathological features was analyzed using the Fisher's exact test.Spearman's correlation analysis was used to analyze the correlation between METTL3-AREG gene expression.Western blotting and transwell assay results were quantified using ImageJ software.A P-value < 0.05 was considered statistically significant (*P < 0.05).

Results
AREG is highly expressed in PC and negatively correlates with prognosis.RNA microarray technology was applied to detect differentially expressed RNAs in three pairs of PC tissues and paraneoplastic noncancerous tissues.A total of 8345 different genes were identified by strict screening with fold change (FC) ≥ 2 and false discovery rate (FDR) < 0.05.Of these genes, there were 5131 upregulated RNAs and 3214 downregulated RNAs, with AREG expression increasing approximately 30-fold in PC tissues (FC = 39.63,P = 0.003, FDR = 0.033) (Fig. 1A).To identify the expression levels of AREG in PC tissues, we first assessed AREG expression in PC tissues using the GEPIA database (http:// gepia.cancer-pku.cn/) and observed that AREG mRNA levels were markedly upregulated in PC tissues (T = 179) when compared to adjacent normal tissues (N = 171) (Fig. 1B), and high AREG expression was related to poor prognosis (Fig. 1C).We next analyzed AREG mRNA levels in 36 pairs of PC tissues and revealed that AREG mRNA was upregulated in cancer tissues compared to normal tissues adjacent to the cancer (Fig. 1D).Then, PC tissues were divided into two groups based on the expression level of AREG, and it was found that the level of AREG expression correlated with clinical stage and pathological grade (Table 1).Furthermore, immunoblot analysis confirmed the expression level of AREG protein in PC tissues and paracancerous tissues, which was consistent with expression at the mRNA level (Fig. 1E).In PC cells, we showed that AREG was upregulated at the mRNA level in PC cells compared to pancreatic epithelial cell lines (HPNE) (Fig. 1F).MIA PaCa-2 and PANC-1 (which had a relatively high expression of AREG) were then chosen for the next stage of the research.

AREG knockdown suppresses proliferation, migration, invasion, and EMT of PC cells in vitro.
Secondly, we explored the biological functions of AREG.qPCR was used to verify that AREG was knocked down in pancreatic cancer cell lines (Fig. 2A) and that the expression of E-cadherin was upregulated and vimentin was decreased in AREG siRNA-transfected cells (Fig. 2B).The silencing of AREG decreased the growth, migration, and invasion of both PC cell lines, as shown by the results of CCK-8 assays and transwell assays (Fig. 2C-E).In general, the downregulation of AREG inhibited PC proliferation, migration, invasion, and EMT in vitro.

METTL3 regulates m 6 A upregulation on AREG mRNA.
To explore the mechanism of AREG upregulation in PC, we conducted methylation RNA immunoprecipitation (MeRIP) and RNA sequencing on METTL3 knockdown PANC-1 cells and NC cells.Heat map analysis from RNA-seq data listed genes that were differentially expressed in the NC or METTL3 knockdown-treated PANC-1 cells (Fig. 3A).A four-quadrant plot shows that METTL3 knockdown resulted in 357 upregulated hypermethylated RNAs, 295 downregulated hypermethylated RNAs, 244 upregulated demethylated RNAs, and 541 downregulated demethylated RNAs, with a fold change (FC) > 1 (Fig. 3B).And m6A signaling was abundant around the stop codon and 3'UTR of the mRNAs (Fig. 3C).Further screening was carried out to analyze genes with changes in the m6A peak detected by MeRIPseq (35 genes, fold change > 2) and genes whose expression was upregulated or downregulated in METTL3deficient PANC-1 cells using RNA-seq (44 genes, fold change > 2), where nine genes, including AREG (Fig. 3D), were overlapping.Then, further analysis showed that the m6A peak in AREG was located in chromosome 4 (74,444,099-74,455,243) (Fig. 3E).

METTL3-dependent upregulation of AREG mRNA by m 6 A suppressed AREG expression.
In addition, RIP assay demonstrated that METTL3 bound to AREG mRNA in PANC-1 and MIA PaCa-2 cells (Fig. 4A).MeRIP-qPCR assay indicated that in PANC-1 and MIA PaCa-2 cells, m 6 A levels of AREG mRNA were markedly reduced when METTL3 was knocked down (Fig. 4B).Furthermore, we assessed the rate of RNA decay in METTL3-knocked down PC cells and corresponding control cells.qPCR results showed that the knockdown of METTL3 significantly shortened the half-life of AREG mRNA (Fig. 4C).To confirm the clinical relevance of METTL3 and AREG expression, we performed data analysis on 36 pairs of tissue specimens, and the results showed that AREG mRNA expression was positively correlated with METTL3 mRNA expression in PC (Fig. 4D).Meanwhile, the mRNA and protein expressions of AREG in cells decreased after METTL3 were knocked down, and METTL3 overexpression increased the expression of AREG in cells (Fig. 4E, F).Overall, our data demonstrated that METTL3-mediated upregulation of AREG mRNA m 6 A contributed to the expression and stability of AREG mRNA in PC cells.
METTL3 was highly expressed in human PC tissues and cell lines.The above results indicated that METTL3 medicated the upregulation of AREG mRNA in an m 6 A-dependent manner.We then examined www.nature.com/scientificreports/ the expression in PC tissues using RNA microarray technology in three pairs of pancreatic cancer tissues and adjacent non-cancerous tissues.As shown in Fig. 5A, several differentially expressed RNAs were identified, among which the expression of METTL3 in PC tissues was increased by approximately threefold (FC = 3.00, P = 0.001, FDR = 0.03) (Fig. 5A).Meanwhile, we identified the expression of METTL3 in PC tissues and cells by qPCR and found that METTL3 was expressed at higher levels in 36 pairs of PC tissues than in paired adjacent normal tissues (Fig. 5B).Besides, the level of METTL3 protein was markedly raised in the tissues of representative PC patients compared to normal tissues (Fig. 5C).Consistent with this finding, METTL3 mRNA (Fig. 5D) and protein (Fig. 5E) were expressed at higher levels in three different PC cell lines than in normal pancreatic ductal epithelial cells.These results suggest that METTL3 expression is upregulated in PC.To determine the role of METTL3 in cell proliferation, we depleted METTL3 in PANC-1 and MIA PaCa-2 PC cells by expressing METTL3 short hairpin RNA (ShRNA) (Fig. 5F, G).We showed that depletion of METTL3 reduced the proliferation (Fig. 5H), migration (Fig. 5I) and invasion (Fig. 5J) of these cells.
METTL3 is a target gene for miR-3a-3p.We further explored the upstream mechanisms of METTL3/ AREG in PC.As miRNAs are widely recognized as players in cancer progression and regulators of gene expression, potential miRNAs targeting METTL3 were searched through databases such as RNAinter (http:// www.rnain ter.org/ search/) and miRDB (https:// mirdb.org/).We found that miR-33a-3p interacted with METTL3 at a high level (Fig. 6A, B).And reviewing the relevant literature, we found few studies on the function and role of miR-33a-3p in PC.Therefore, we selected miR-33a-3p for further study.Given that miR-33a-3p expression levels were significantly lower in PC tissues compared to normal tissues (Fig. 6C).Additionally, we selected PANC-1 and MIA PaCa-2 cells and transfected their mimics for subsequent experiments to explore the potential biological function of miR-33a-3p in vitro (Fig. 6D).As shown in the Fig. 6E, miR-33a-3p expression in PANC-1 and MIA PaCa-2 cells was significantly higher in the miR-33a-3p mock-transfected group than in the miR-NCtransfected group.www.nature.com/scientificreports/miR-33a-3p downregulates METTL3 expression by directly targeting its 3′-UTR and inhibits proliferation, migration, and invasion.We confirmed that transfection of miR-33a-3p mimics in PC cell lines attenuated the mRNA (Fig. 7A) and protein (Fig. 7B) level expression of METTL3 and AREG.To explore the relationship between METTL3 and miR-33a-3p, miR-33a-3p and METTL3 binding sequences and mutant sites are shown (Fig. 7C).Additionally, miR-33a-3p overexpression significantly reduced the luciferase activity of METTL3 WT but not the METTL3 MUT luciferase activity (Fig. 7D).The role of miR-33a-3p in cell proliferation, migration, and invasion in vitro was analyzed using CCK-8 and Transwell assays.The results showed that the miR-33a-3p mock-transfected group significantly inhibited the proliferative capacity of both PANC-1 and MIA PaCa-2 cells compared to miR-NC (Fig. 7E).In addition, miR-33a-3p overexpression significantly reduced the migratory and invasive abilities of cells (Fig. 7F, G).In summary, these data indicate that miR-33a-3p has an inhibitory effect on PC cell growth and metastasis in vitro.
miR-33a-3p suppresses PC progression via the miR-33a-3p/METTL3 axis.To exploit whether miR-33a-3p exerts its inhibitory effect by regulating METTL3 expression, a rescue assay between miR-33a-3p and METTL3 was performed in PC cells.We found that the expressions of METTL3 and AREG were largely boosted in PANC-1 and MIA PaCa-2 transfected with miR-33a-3p mimics by the introduction of pCMV3-METTL3 at both RNA and protein levels (Fig. 8A, B).CCK-8 assay showed that overexpression of METTL3 significantly reversed of the role of miR-33a-3p mimics on cell proliferation capacity (Fig. 8C).Similarly, transwell analysis showed that overexpression of METTL3 significantly reversed the role of miR-33a-3p mimics on cell migration and invasion capacity (Fig. 8D, E).These results suggest that miR-33a-3p may attenuate the migratory and invasive behavior of PC cells, at least in part, by down-modulating the expression of METTL3.

Discussion
In short, our findings revealed that low expression of miR-33a-3p in PC was correlated with high expression of METTL3.High expression of METTL3-mediated m6A methylation of AREG enhanced the stability of AREG RNA, led to the increased expression of AREG, and further induced EMT, which promoted the aggressive and metastatic processes of PC (Fig. 8F).Oncogenic or tumor suppressor genes that are dysregulated in PC have been suggested as potential biomarkers of PC progression.Previous studies have shown that AREG is highly regulated and acts as an oncogene in a range of cancers 23 .AREG has been described as being highly expressed in breast cancer and facilitating tumor progression.In the present study, we explored the role of AREG in PC in the same way and identified that AREG is highly expressed in PC and is strongly related to the prognosis of the disease.Functional assays indicated that silencing of AREG impeded PC cell proliferation, invasion, and migration, but the reasons for its high expression and the molecular mechanism that facilitates tumor progression remain unclear.
Recently, m 6 A modifications in cancer progression have become a focus of research 24 .The most frequent internal mRNA modifications involve N6-methyladenosine (m 6 A), N1-methyladenosine (m 1 A), and 5-methylcytosine (m 5 C) 18 .As the most common catalyst of m 6 A modification, METTL3 has been shown to promote the progression of various cancers, such as breast, lung and gastric cancers [25][26][27][28] .It was previously indicated that METTL3 could be induced by m 6 A modification of target mRNAs 29 .Here, we identified m 6 A modifications in AREG by sequencing and hypothesized that AREG overexpression in PC is mediated by its m 6 A modifications.We then confirmed that METTL3 was upregulated in PC cells and positively modulated AREG expression.In addition, we demonstrated that METTL3 promotes m 6 A on AREG mRNA.It is recognized that m 6 A modification can affect mRNA splicing, stability, and translation.Therefore, we identified that METTL3 stabilized AREG mRNA.However, the exact reading protein that recognizes the m 6 A site to better enhance the stability of AREG mRNA needs to be further determined.
The modulation of gene expression in PC by miRNAs has been well documented 30 .Axiomatically, miRNA interacts with mRNA at the 3′UTR to suppress gene expression 31 .In this study, we identified miR-33a-3p as a potential interactor with METTL3 through the miRDB and RNAinter databases.Previously, miR-33a-3p overexpression in hepatocellular carcinoma was found to repress the proliferation and migration of hepatocellular carcinoma cells, indicating that miR-3a-3p functions as a tumor suppressor gene 32 .Therefore, we found that miR-33a-3p could regulate METTL3 by bioinformatics prediction.We further investigated the regulatory relationship between miR-33a-3p and METTL3.Dual fluorescence assays indicate that miR-33a-3p can target METTL3 and inhibit its expression.Finally, the rescue assays showed that miR-33a-3p inhibited PC progression via METTL3.Considering that miRNAs usually have multiple downstream targets, we cannot exclude that miR-33a-3p may also mediate PC progression via other potential targets.
There are some limitations in this study.Firstly, the 36 patients with PC should be followed up for prognosis analysis.Secondly, animal study should be conducted to further confirm the results in vitro.
To conclude, our study reveals that miR-33a-3p inhibits m 6 A-induced stabilization of AREG by targeting METTL3, which plays a key role in the aggressive progression of PC.AREG could be a potential target for PC treatment.

Conclusion
Our results indicate that miR-33a-3p modulates METTL3-mediated m 6 A-AREG stability and thereby suppresses pancreatic cancer invasion and metastasis.The results of this study help provide new insights for the development of new therapeutic strategies for PC.

Figure 1 .
Figure 1.AREG is highly expressed in PC and negatively correlates with prognosis.(A) The expression level of AREG mRNA was analyzed using RNA microarray containing 3 pairs of PC tissues and paracancerous tissues.(B) AREG mRNA expression was analyzed in PC tissues using an online tool GEPIA.(C) The association between AREG mRNA expression and overall survival in PC was analyzed using an online tool Kaplan-Meier plotter.(D) AREG mRNA was analyzed by qPCR in 36 pairs of PC tissues collected from our hospital.(E) AREG protein expression was examined in four pairs of PC tissues by western blotting.(F) AREG expression was examined in PC cells and HPNE cells by qPCR.Data are expressed as the mean ± SD. *P < 0.05, **P < 0.01, ***P ≤ 0.001, and ****P ≤ 0.001 compared with the control.Abbreviations: PC pancreatic cancer, T tumor, N normal, qPCR quantitative polymerase chain reaction, HPNE pancreatic epithelial cell line.

Figure 2 .
Figure 2. AREG knockdown suppresses proliferation, migration, invasion, and EMT of PC cells in vitro.(A) AREG was knocked down using siRNA and the inhibitory efficacy was examined by qPCR.(B) The effects of AREG knockdown on the expression of E-cadherin and Vimentin in PC cells were examined by western blotting.(C) The effect of AREG knockdown on the proliferative capacity of PC cells was examined by CCK8 assay.(D, E) The effects of AREG knockdown on the migratory (D) and invasive (E) capacity of PC cells were examined by transwell assay.Scale bar: 100 μm.Data are expressed as the means ± SD. *P < 0.05, **P < 0.01, ***P ≤ 0.001, and ****P ≤ 0.001 compared with the control.Abbreviations: PC pancreatic cancer, CCK8 cell counting kit-8, qPCR quantitative polymerase chain reaction.

Figure 3 .
Figure 3. METTL3 mediates m 6 A upregulation on AREG mRNA.(A) Differentially expressed genes in PANC-1 cells transfected by NC or METTL3 shRNA were analyzed by RNA-sequencing.(B) An interaction analysis of MeRIP-sequencing and RNA-sequencing to show the distribution of m 6 A peaks (fold change > 1) in METTL3 knockdown PANC-1 cells compared to control cells.(C) The proportion of m 6 A peaks distributed in the 5′ untranslated regions (5′-UTR), the start codon, the coding region (CDS), the stop codon, and the 3′ untranslated regions(3′-UTR) was calculated.(D) Venn diagram was employed toidentify the mRNAs with both changes on m6A peaks and gene expression.(E) IGV was used to analyzed the location of m 6 A peak in AREG.Abbreviations: MeRIP Methylated RNA immunoprecipitation, NC negative control, m 6 A N6-methyladenosine, 5′-UTR 5′untranslated regions, 3′-UTR 3′untranslated regions, CDS the coding region.

Figure 4 .
Figure 4. METTL3-dependent upregulation of AREG mRNA by m 6 A suppressed AREG expression.(A) RIP assay was used to confirm the binding capacity between AREG mRNA and METTL3 in PC cells.(B) The effect of METTL3 knockdown on the m 6 A levels in AREG in PC cells was analyzed by MeRIP-qPCR.(C) The effect of METTL3 knockdown on the half-life of AREG in PC cells was examined by qPCR.(D) Correlation between AREG mRNA and METTL3 mRNA in PC tissues was analyzed by Spearman.(E) The effects of METTL3 overexpression or knockdown on the expression of AREG mRNA were analyzed by qPCR.(F) The effects of METTL3 overexpression or knockdown on the expression of AREG protein were analyzed by western blotting.Data are expressed as the means ± SD. *P < 0.05, **P < 0.01, and ****P ≤ 0.001 compared with the control.Abbreviations: PC pancreatic cancer, NC negative control, RIP RNA immunoprecipitation, MeRIP Methylated RNA immunoprecipitation, qPCR quantitative polymerase chain reaction, m 6 A N6-methyladenosine.

Figure 5 .
Figure 5. METTL3 was highly expressed in human PC tissues and cell lines.(A) The level of METTL3 mRNA was analyzed using RNA microarray in 3 pairs of PC tissues and paracancerous tissues.(B) METTL3 mRNA was analyzed in 36 pairs of PC tissues collected from our hospital by qPCR.(C) METTL3 protein expression was examined in four pairs of PC tissues by western blotting.(D, E) The expression of METTL3 in PC cells and HPNE cells was detected by qPCR (D) and western blotting (E).(F, G) METTL3 silencing cell lines were established by infecting PC cells with METTL3 shRNA lentivirus.The inhibitory efficiacy was analyzed by qPCR (F) and western blotiing (G).(H) The effect of METTL3 knockdown on the proliferative capacity of PC cells was examined by CCK8 assay.(I, J) The effects of METTL3 knockdown onthe migratory (I) and invasive (J) capacity of PC cells were examined by transwell assay.Scale bar: 100 μm.Data are expressed as the means ± SD. *P < 0.05, **P < 0.01, ***P ≤ 0.001, and ****P ≤ 0.001 compared with the control.Abbreviations: PC pancreatic cancer, NC negative control, T tumor, N normal, qPCR quantitative polymerase chain reaction, CCK8 cell counting kit-8, FC fold change.

Figure 6 .
Figure 6.METTL3 is a target gene for miR-3a-3p.(A, B) The target scores of selected miRNAs that bind to the METTL3 gene.miR-33a-3P was selected from the RNAinter and miRDB databases and predicted to have the highest target scores.(C) The expression of miR-33a-3p in 36 pairs of PC tissues collected from our hospital was examined by qPCR.(D) The expression of miR-33a-3p in PC cells and HPNE cells was examined by qPCR.(E) The overexpression of miR-33a-3p was conducted by transfecting PC cells with miR-33a-3p mimics.And the efficacy was analyzed by qPCR.Data are expressed as the means ± SD. *P < 0.05, **P < 0.01, and ***P ≤ 0.001 compared with the control.Abbreviations: PC pancreatic cancer, NC negative control, qPCR quantitative polymerase chain reaction.

Figure 8 .
Figure 8. miR-33a-3p suppresses PC progression via the miR-33a-3p/METTL3 axis.(A) The reversal effect of METTL3 overexpression on AREG mRNA level in miR-33a-3p mimics PC cells was examined by qPCR.(B) The reversal effect of METTL3 overexpression on AREG protein level in miR-33a-3p mimics PC cells was examined by western blotting.(C) The reversal effect of METTL3 overexpression on the proliferative capacity in miR-33a-3p mimics PC cells was examined by CCK8 assay.(D, E) T The reversal effects of METTL3 overexpression on the migratory (C) and invasive (D) capacity in miR-33a-3p PC cells were examined by transwell assay.Scale bar: 100 μm.(F) Schematic model on the role of miR-33a-3p in regulating METTL3mediated AREG stability in an m 6 A-dependent manner, thus inhibiting pancreatic cancer invasion and metastasis.Data are expressed as the means ± SD. *P < 0.05, **P < 0.01, ***P ≤ 0.001, and ****P ≤ 0.001 compared with the control.Abbreviations: PC pancreatic cancer, NC negative control, qPCR quantitative polymerase chain reaction.

Methylated RNA immunoprecipitation sequencing and RNA sequencing. METTL3
To measure the RNA stability of PC cells and METTL3 knockdown PC cells, Actinomycin D (MCE, HY-17559, USA) was added to the cells at 5 μg/ml.After incubation for a certain period of time, the cells were collected at different time points, and subsequently, RNA was extracted from each cell sample for real-time quantitative polymerase chain reaction and normalized by the GAPDH method.Vol:.(1234567890)