PRMT5 mRNA expression levels are upregulated in lung cancer cells. Lung cancer is a leading cause of cancer-related deaths, with more than 85% of lung cancer cases classified as non-small cell lung cancer (NSCLC)25. The most common genetic driver alteration in NSCLC is mutations of epidermal growth factor receptor (EGFR) that result in constitutive activation. STAT3 is one of the key signalling mediators downstream of EGFR, and STAT3 is indeed constitutively activated in about 50% of NSCLC tissues and cell lines26, 27, 28. In addition to this pathway, the HEDGEHOG (HH) pathway is known to be activated in NSCLC29. Experimental evidence suggests that STAT3 activation is enhanced downstream of the HH signal30. The HH signalling pathway plays a central role in the growth, patterning, and morphogenesis of organisms, and is also involved in adult organ homeostasis, stem cell maintenance, and cancer development31, 32, 33. We previously found that HH signalling in cancer cells can induce expression of PRMT5, with arginine methylation by PRMT5 being important for activation of the effector transcription factor GLI1 protein34. PRMTs are enzymes that can methylate arginine residues in certain proteins and have been reported to regulate gene expression, signal transduction, and mRNA splicing fidelity35. Of the nine members of the PRMT protein family, PRMT5 is believed to be involved in cancer through various functions, including tumour growth36, stem cell reprogramming and maintenance37, and immune evasion38. Furthermore, overexpression of PRMT5 has been observed in a variety of cancer types39. Therefore, we first analysed PRMT5 expression levels in NSCLC samples. As shown in Fig. 1a/b, PRMT5 mRNA expression levels were enhanced in the two most common NSCLC subtypes, lung adenocarcinoma and lung squamous cell carcinoma, as determined by RNA sequencing (RNA-seq) transcriptomic data from The Cancer Genome Atlas (TCGA) database40. PRMT5 expression levels and copy number were correlated, with higher PRMT5 expression being associated with worse overall survival in lung adenocarcinoma patients (Supplementary Fig. 1a, b). Moreover, mRNA expression of methylosome protein 50 (MEP50), an adaptor molecule between PRMT5 and its substrates41, 42, was also enhanced (Supplementary Fig. 1c, d). This suggests that the functions of methylosomes containing PRMT5 are increased in NSCLC. PRMT5 expression levels were higher in NSCLC compared with other well-known HH signalling-enhanced cancers, such as esophagogastric, colorectal, prostate, and pancreatic cancers43, 44, 45, 46 (Fig. 1c). Ranking of PRMT5 expression in cancer cell lines with the Cancer Cell Line Encyclopedia (https://sites.broadinstitute.org/ccle/) also showed high PRMT5 expression levels in NSCLC cells (Fig. 1d). These results suggest that methylation by PRMT5 plays an important role in NSCLC, including cancer progression.
STAT3 activates PRMT5 gene transcription. Because PRMT5 mRNA expression levels were upregulated in NSCLC samples, we next analysed the relationship between the transcription factor STAT3, which is constitutively activated in NSCLC, and PRMT5 expression. Using RNA interference (RNAi) to knock down STAT3 expression (siSTAT3), we observed reduced PRMT5 protein expression levels in normally growing murine diploid NIH3T3 cells. STAT3 expression was induced by v-SRC, a typical oncogenic signalling molecule known to activate STAT347 (Fig. 2a), suggesting that STAT3 is involved in the regulation of PRMT5 expression in cancer cells. Moreover, PRMT5 mRNA expression was upregulated by STAT3 activator IL-6, while this induction was suppressed by STAT3 knockdown (Fig. 2b). EGF-induced PRMT5 mRNA expression was also inhibited by STAT3 inhibitor LLL-12 (Fig. 2c), suggesting that PRMT5 gene transcription is activated by STAT3. Indeed, the PRMT5 promoter, which includes putative STAT3 binding sites48, was activated by a constitutively active mutant of STAT313 (Fig. 2d–f) and STAT3 protein directly bound with the PRMT5 promoter region (Fig. 2g, h). Binding of STAT3 to the PRMT5 promoter was also confirmed in the UCSC Genome Browser database (Fig. 2i). Finally, STAT3-dependent constitutive expression of PRMT5 was also observed in the NSCLC cell lines HCC827 (EGFR mutant, adenocarcinoma; Fig. 2j) and A549 (EGFR wild-type, adenocarcinoma; Fig. 2k).
PRMT5 induces tumour growth and CSC maintenance in NSCLC cells. We next analysed the role of PRMT5 in HCC827 and A549 NSCLC cells. The proliferation rates of cells with knockdown of PRMT5 (Fig. 3a) and its adaptor molecule MEP50 (Fig. 3b) were slightly attenuated compared with that of the parental cells (Fig. 3c, d). Cancer arises from CSCs, also known as cancer-initiating cells, which have tumourigenic potential, self-renewal properties, and long-term tumour repopulating activity14, 49, 50. It has been shown that CSCs can form spheres under low attachment culture conditions in media containing growth factors51, 52. Indeed, these cells formed primary spheres and secondary spheres, characteristic of stem cells53, but PRMT5 knockdown reduced their number (Fig. 3e–h). Certain reprogramming factors, such as SOX2 and OCT4, have induced pluripotent stem cell (iPSC)-inducing activity and are a hallmark of CSCs15. In sphere cells, the expression levels of these reprogramming factors were elevated (Supplementary Fig. 2). In addition, the population of highly activated ALDH cells, which is a hallmark of lung CSCs, was reduced in PRMT5 knockdown cells (Fig. 3i, j). These results suggest that PRMT5-mediated methylation is important for CSC maintenance. Like the decreased number of sphere-forming cells, the growth of A549 cell tumours transplanted in nude mice were also reduced in PRMT5 knockdown cells (Fig. 3k, l). Furthermore, this tumour growth inhibitory effect not only seen by constitutive knockdown, but also by inducible knockdown of PRMT5 and MEP50 (Fig. 3m, n).
STAT3 activity is upregulated by PRMT5-mediated methylation. Next, we analysed the functional effects of the PRMT5 and STAT3 interaction in cancer cells. The number of sphere-forming cells in HCC827 and A549 cells was increased by IL-6 stimulation, but this increase did not occur in PRMT5 knockdown cells (Fig. 4a, b). Related to this, the induction of MCP-1 mRNA, which is known to be regulated by the IL-6-STAT3 pathway54, was suppressed by PRMT5 knockdown (Supplementary Fig. 3a, b). Furthermore, PRMT5 knockdown also suppressed the increase of sphere-forming cell number following hepatocyte growth factor (HGF) treatment, which is known to induce oncogenic transformation via the MET-STAT3 pathway55 (Supplementary Fig. 3c, d). These results suggest that STAT3 is possibly involved in CSC maintenance downstream of PRMT5. Indeed, IL-6-induced Y705 phosphorylation was suppressed in PRMT5 knockdown cells (Fig. 4c), as was STAT3 transcriptional activity (Fig. 4d, e). However, PRMT5 knockdown did not affect the transcriptional activity of transcription factors NF-kB and HIF156, which are known to be activated in other cancer types (Supplementary Fig. 3e, f). In A549 cells, we found that PRMT5 was bound with STAT3 (Fig. 4f). Furthermore, STAT3 itself was arginine methylated in this complex, and its methylation was inhibited by the PRMT5 inhibitor (Fig. 4f). Formation of this complex with STAT3 and PRMT5 was suppressed by MEP50 knockdown (Fig. 4g), suggesting that STAT3 is methylated in the methylosome complex formed by PRMT5/MEP5041, 42. In addition, PRMT5 was mainly located in the cytoplasm, implying that STAT3 methylation occurs there (Supplementary Fig. 3g, h). Previously, it has been shown that STAT3 plays an important role in v-SRC transformation, as dominant negative STAT3 mutants can inhibit v-SRC transformation57. As shown in Supplementary Fig. 4a–f, tumour suppressor p53-deficient embryonic fibroblasts expressing v-SRC showed foci formation58, a hallmark of cancer cells. However, this ability was reduced by PRMT5, MEP50, or STAT3 knockdown, suggesting that the PRMT5/MEP50/STAT3 pathway is important for CSC generation.
PRMT5-mediated STAT3 methylation is required for tumourigenesis and CSC maintenance in NSCLC cells. Candidate arginine methylation sites in STAT3 were predicted by PRmePRed (Fig. 5a, b). Then, wild-type STAT3 and STAT3 mutants were transiently expressed in STAT3 knockout BEAS2B cells (Supplementary Fig. 5). Amino acid substitution (arginine to lysine) at candidate methylation sites revealed that arginine 609 (R609K) was clearly methylated and arginine 518 (R518K) was weakly methylated by PRMT5. Correlated with the decrease in methylation by the arginine mutant, phosphorylation of STAT3 Y705 was also decreased (Fig. 5c). Furthermore, R609K mutant STAT3 lost most of its transcriptional activity in A549 cells (Fig. 5d), and stably expressing STAT3-R609K could also suppress Y705 phosphorylation in STAT3 knockout HCC827 cells (Fig. 5e).
Finally, we examined the impact of STAT3-R609K on tumour growth and CSC maintenance. STAT3 knockdown HCC827 cells showed slightly decreased cell proliferation rates, but these were increased to levels similar to the parental cells after wild-type STAT3 was expressed. However, STAT3-R609K mutant-expressing cells did not show recovered proliferation rates (Fig. 5f). Moreover, STAT3 knockdown cells had a markedly decreased number of sphere-forming cells, and the R609K mutant showed a reduced recovery of sphere-forming cells compared with wild-type STAT3 (Fig. 5g). Furthermore, tumour growth in nude mice was markedly reduced in STAT3 knockdown cells and was restored by wild-type STAT3 expression, similar to the parental cells, while this recovery was not observed in mutant STAT3-R609K-expressing cells (Fig. 5h–j). These results suggest that methylation of STAT3 R609 by PRMT5 is important for CSC maintenance and tumour growth in NSCLC.