Advertisement
Cancer Genes and Networks

Histone Demethylase JMJD1A Promotes Tumor Progression via Activating Snail in Prostate Cancer

Dong-E Tang, Yong Dai, Ling-Ling Fan, Xin-Yan Geng, De-Xue Fu, Hao-Wu Jiang and Song-Hui Xu
DOI: 10.1158/1541-7786.MCR-19-0889 Published May 2020

Abstract

The histone demethylase JMJD1A plays a key functional role in spermatogenesis, sex determination, stem cell renewal, and cancer via removing mono- and di-methyl groups from H3K9 to epigenetically control gene expression. However, its role in prostate cancer progression remains unclear. Here, we found JMJD1A was significantly elevated in prostate cancer tissue compared with matched normal tissue. Ectopic JMJD1A expression in prostate cancer cells promoted proliferation, migration, and invasion in vitro, and tumorigenesis in vivo; JMJD1A knockdown exhibited the opposite effects. Mechanically, we revealed that JMJD1A directly interacted with the Snail gene promoter and regulated its transcriptional activity, promoting prostate cancer progression both in vitro and in vivo. Furthermore, we found that JMJD1A transcriptionally activated Snail expression via H3K9me1 and H3K9me2 demethylation at its special promoter region. In summary, our studies reveal JMJD1A plays an important role in regulating proliferation and progression of prostate cancer cells though Snail, and thus highlight JMJD1A as potential therapeutic target for advanced prostate cancer.

Implications: Our studies identify that JMJD1A promotes the proliferation and progression of prostate cancer cells through enabling Snail transcriptional activation, and thus highlight JMJD1A as potential therapeutic target for advanced prostate cancer.

This article is featured in Highlights of This Issue, p. 669

Introduction

Prostate cancer is one of the most commonly diagnosed malignancies and the second leading cause of cancer-related death in American men (1). Androgen deprivation therapy (ADT) remains the primary clinical treatment for patients in the early stage of prostate cancer (2). However, most prostate cancer tumors become resistant to ADT and progress to a lethal stage called castration-resistant prostate cancer (CRPC), which is characterized by aggressive growth and distal organ metastasis, and is incurable (3). Current therapies for CRPC, including treatment with new-generation androgen receptor (AR)-pathway inhibitors (abiraterone or enzalutamide) and chemotherapy extend patient life by only a few month and benefit a little for the patients with metastatic CRPC (4, 5). Thus, exploring the underlying mechanism in metastasis of prostate cancer is crucial to develop therapeutic targets to treat or prevent prostate cancer metastasis.

H3K9 methylation is the mark of heterochromatin and often functions as a repressive epigenetic modification to regulate gene expression (6). JMJD1A (also named KDM3A) demethylates H3K9me1 and H3K9Me2 to increase gene transcription, and with wide functional roles in spermatogenesis (7), sex determination (8), stem cell renewal (9), and cancer (6, 10). JMJD1A overexpression has been observed in a subset of primary human cancers including liver cancer and sarcoma (11). We previously found that JMJD1A promoted prostate cancer cells' growth as an AR coactivator or by upregulating levels of AR-V7, an AR splice variant 7 (6, 11). However, the tumor-promoting role and molecular mechanisms of JMJD1A in prostate cancer remains largely unknown.

Metastasis is a process by which cancer cells invade and migrate, and can be induced by regulators such as twist, Snail, and slug (12). Snail, a snail family zinc finger protein, has been found to promote cell survival and movement through induction of epithelial–mesenchymal transition (EMT) or independent of EMT (13–15). Snail has been found to overexpress in a set of cancers, one study reported that overexpression of Snail was related with poor prognosis in patient with breast cancer (16). Moreover, nuclear Snail expression was found to be related to tumour progression (17, 18). In prostate cancer, several studies have reported that Snail expression was elevated in metastatic prostate cancer cell lines, Snail transcriptional regulatory network was important for the progression and metastasis of cancer (19, 20) and Snail overexpression could induce EMT in prostate cancer cells (21). Importantly, the mechanism for upregulating Snail in prostate cancers is still unclear.

Here, we found that JMJD1A protein and mRNA levels are elevated in prostate cancer compared with normal tissues. Furthermore, JMJD1A overexpression significantly promotes prostate cancer cell proliferation, migration, and invasion, while JMJD1A knockdown exhibits the opposite effects. More importantly, these biologic effects mediated by JMJD1A might be dependent on the Snail transcription factor. Identifying these mechanisms highlights JMJD1A is closed with prostate cancer development through Snail, suggesting that JMJD1A–Snail axis could be targeted for prostate cancer therapies.

Materials and Methods

Cell culture, constructs, and antibodies

RV1 and LNCAP cells were purchased from ATCC and cultured in DMEM (Gibco, 110 Life Technologies) supplemented with 10% FBS (Biological Industries), 100 U/mL penicillin/streptomycin (Gibco, Life Technologies). Cells were grown at 37°C with 5% CO2. All cell lines were authenticated by short tandem repeat analysis at China Center for Type Culture Collection and all the cell lines were tested and authenticated by karyotyping analysis on January 3, 2018, and confirmed by National Infrastructure of Cell Line Resources of China. Cell lines were tested Mycoplasma free. Cells used for experiments were between 20–30 passages from thawing. Flag-JMJD1A and Myc-JMJD1A were described previously (6). Flag-Snail (Plasmid #16218) was purchased from Addgene.

JMJD1A (12835-1-AP) was purchased from Proteintech. EMT marker antibodies (#9782) were purchased from Cell Signaling Technology. H3K9me2 (07-441) was purchased from EMD Millipore. H3K9me1 (ab9045) antibodies, anti-Ki67 antibody (ab15580), and anti-Cleaved Caspase-3 antibody (ab2302) were purchased from Abcam. Anti-β-actin antibodies (A1978) were purchased from Sigma.

Soft-agar colony formation assays

The indicated cells (RV1, 1 × 104 cells and LNCAP, 1.5 × 104 cells) were mixed with agar to a final concentration of 0.4% and layered on top of 0.8% agar in 6-well plates. Cells were maintained in the normal growth media. Triplicate plates were incubated at 37°C for 3 weeks. Photographs of colonies were obtained using a microscope. The independent experiments were repeated at least three times. The results are shown as means ± SD (n = 3), and ANOVA was used for statistical analysis.

Chromatin immunoprecipitation assay

Cells were cross-linked using 1% formaldehyde for 10 minutes at room temperature, and then quenched with 5 mol/L glycine. Cell nuclei were extracted with buffer A (10 mmol/L HEPES, pH 7.4, 10 mmol/L KCl, 1.5 mmol/L MgCl2, 10% glycerol, 1 mmol/L DTT, and 1× protease cocktail). Nuclear extracts were solubilized in the lysis buffer (50 mmol/L Tris-HCl, pH 8.0, 1% SDS, 10 mmol/L EDTA, and 1× protease cocktail), and sonicated to obtain 500-bp chromatin fragments. Chromatin (100 μg) was 10-fold diluted and incubated with 5 μg of antibodies overnight followed by incubation with 30 μL of protein A/G beads for 4 hours at 4°C. After four washes, cross-linking was reversed, and DNA was purified using spin columns and subjected to qPCR analysis for the indicated genes. The primers are listed in Supplementary Materials and Methods. The chromatin immunoprecipitation (ChIP) assays were performed in triplicate and independent experiments were repeated at least three times. Data were calculated as percentage of input and presented as means ± SD (n = 3), and ANOVA was used for statistical analysis.

IHC

Prostate tumor tissue samples and matched adjacent nontumoral prostate tissue samples were collected at The Second Clinical Medical College of Jinan University (Guangdong, P.R. China). Tissue sample collection was approved by the medical ethics committee of Shenzhen People's Hospital (Shenzhen, P.R. China). Tissue microarray chips containing adjacent nontumoral prostate tissue samples and prostate tumor tissue samples were obtained from Shanghai OUTDO Biotech Co., Ltd.. IHC staining and quantification were conducted as described previously (6). The score of JMJD1A expression at the protein level was classified semiquantitatively as follows: no staining, 0 points; weak staining, 1 point; moderate staining, 2 points; and strong staining, 3 points; while 2 or 3 points were defined as high expression. The proportion of tumor cells was graded from 0 to 2: 0 points (<5%), 1 point (5%–50%), and 2 points (>50%). The sum of scores represented the expression levels: 0 to 2 points as low expression and 3 to 4 points as high expression. The mean score from two pathologists was used as the final immunostaining score. The IHC score was calculated as described previously (6). The χ2 test was used for statistical analysis of JMJD1A expression.

RV1 cell xenograft model

Seven-week-old female BALB/c nude mice were purchased from Charles River Laboratories in China and housed under pathogen-free conditions in the animal experiment center of The Second Clinical Medical College of Jinan University (Shenzhen, China). RV1 cells (1 × 106) stably expressing the control vector, stably expressing JMJD1A, control short hairpin RNA (shRNA), or shJMJD1A#1 were injected subcutaneously into female BALB/c nude mice. Every 4 days, the tumor volume and weight were measured following the standard protocol. ANOVA test was used to analysis the data. For the blinding procedures, two researchers in each group performed all the mice experiments. All protocols involving live mice were approved by the Animal Care and Use Committee of The Second Clinical Medical College of Jinan University. Mice were sacrificed after subcutaneous injection on day 22 (control vector and stably expressing JMJD1A) or day 26 (control shRNA and shJMJD1A#1).

Statistical analysis

Data are represented as the mean ± SEM of three independent experiments for cell proliferation and cell clone assays. Data are represented as the mean ± SD of 5 mice for the tumor growth study in xenograft. The two-tailed, unpaired Student t test, ANOVA, and χ2 test were utilized for the statistical analyses (*, P < 0.05; **, P < 0.01). Statistical analysis was performed using Prism 5.0 (GraphPad Software). Data were considered significant when P < 0.05.

Results

JMJD1A expression is elevated in prostate cancer tissue

We previously found that JMJD1A played a key role in AR-FL coactivation and AR-V7 alternative splicing (6, 11). Here, we analyzed the data from Oncomine and found JMJD1A mRNA levels were significantly upregulated in prostate cancer tissues compared with normal tissues (Fig. 1A). Surprisingly, we found elevated JMJD1A mRNA was positively associated with more advanced tumor (T) stage in prostate cancer specimens (Fig. 1B). To confirm this finding, JMJD1A levels in 14 pairs of prostate tumor tissue samples and matched adjacent nontumoral prostate tissue samples were analyzed. As shown in Fig. 1C–E, JMJD1A protein levels and mRNA levels were elevated in most prostate cancer tissues. Supplementary Table S1 shows the clinical pathologic variables and JMJD1A expression in 14 prostate cancer cases for detection by Western blot and RT-PCR analyses. To confirm the relevance of our findings to prostate cancer tissues, we evaluated the expression of JMJD1A in a human prostate cancer tissue microarray (TMA) containing 77 tumor samples and 25 adjacent nontumoral prostate tissue samples. As shown in Fig. 1F, JMJD1A was primarily nuclear localized. IHC results were consistent with prior observations that JMJD1A is significantly increased in prostate cancer samples (Fig. 1G). On the basis of patient data from TMA, we found that a high JMJD1A expression was related with a higher preoperative PSA level (P = 0.00589), Gleason score (P = 0.028), and or lymph node metastasis (P = 0.00585). The results are shown in Supplementary Table S2.

Figure 1.
Figure 1.

JMJD1A expression is elevated in prostate cancer tissue. A, JMJD1A expression was analyzed in normal and prostate cancer tissues from the Oncomine public dataset. *, P = 0.02. B, Correlation between the expression of JMJD1A and the clinical pathologic feature T stage in 102 prostate cancer cases from the Oncomine public dataset. C and D, JMJD1A protein levels in adjacent nontumoral prostate tissue (NT) and prostate cancer tissue samples (T) were detected by Western blot analysis (C). D, JMJD1A protein expression was quantified using Image J software. **, P = 0.0081. E, JMJD1A mRNA expression was detected in adjacent nontumoral prostate tissues and prostate cancer tissues. **, P = 0.0004. F, Representative IHC images of JMJD1A expression in adjacent nontumoral prostate tissue and prostate cancer tissue samples. G, The association between JMJD1A expression and prostate tissue samples was assessed using Pearson χ2 test (**, P = 8.43684E-05).

JMJD1A overexpression promotes prostate cancer cell proliferation, migration, and invasion

To investigate the function of JMJD1A in prostate cancer progression, JMJD1A was overexpressed in two prostate cancer cell lines. As shown in Fig. 2A, Western blot results revealed JMJD1A was overexpressed in RV1 and LNCAP cells. Next, we evaluated the role of JMJD1A in prostate cancer proliferation. We found, using MTT assays, that the cells of JMJD1A overexpression promoted cell growth (Fig. 2B). To further determine the role of JMJD1A in prostate cancer proliferation, we overexpressed JMJD1A in prostate cancer cells and measured the cell proliferation by colony formation assay and soft-agar assay. Concordantly, JMJD1A overexpression also increased the number and size of cell colonies, as determined by colony formation assay or soft-agar assay (Fig. 2C and D). These data demonstrate that JMJD1A overexpression promotes prostate cancer cell proliferation.

Figure 2.
Figure 2.

JMJD1A overexpression promotes prostate cancer cell proliferation, migration, and invasion. A, JMJD1A protein levels were analyzed after JMJD1A overexpression in RV1 cells and LNCAP cells. B, The viability of RV1 and LNCAP cells stably expressing vector or JMJD1A was detected by MTT assay. C, The cells generated as above were seeded on 6-well plates for 2D colony formation assays. Colony numbers of cellular clones with more than 100 cells was measured (mean ± SEM of three independent experiments). The statistical analyses were performed with the ANOVA (right). **, P = 0.0004 (RV1 cells); **, P = 0.00041 (LNCAP cells). D, The cells described above were maintained in soft agar for 3 weeks, and colony number per field was determined (mean ± SEM of three independent experiments). Quantitative analysis of colony formation is presented in the bottom panel. **, P = 0.00043 (RV1 cells); **, P = 0.0022 (LNCAP cells). E and F, Both RV1 cells and LNCAP cells stably expressing vector or JMJD1A were plated in the top chamber with/without Matrigel for 24 hours. Then, cell migration (E) and invasion (F) were assessed by counting the cells that migrated to the underside of the Transwell insert. Mean ± SEM of three independent experiments[**, P = 0.001 (E: RV1 cells); **, P = 0.0005 (E: LNCAP cells); **, P = 0.0002 (F: RV1 cells); **, P = 0.004 (F: LNCAP cells)].

Given that overexpression of JMJD1A was associated with T stage in prostate cancer (Fig. 1B) and metastasis seriously limits the treatments available for prostate cancer. Thus, we assessed whether JMJD1A plays a role in prostate cancer progression. Overexpression of JMJD1A markedly stimulated prostate cancer cell migration by transwell assays (Fig. 2E) and promoted cell invasion through Matrigel (Fig. 2F). These results suggest that JMJD1A participates in the regulation of prostate cancer cell motility and invasion, which is consistent with the functions of JMJD1A that we observed.

JMJD1A knockdown antagonizes prostate cancer cells proliferation, migration, and invasion

To study JMJD1A's role in prostate cancer, we silenced JMJD1A in RV1 cells and LNCAP cells by two different JMJD1A shRNAs. Our results revealed that the shRNAs effectively knocked down JMJD1A in sh-JMJD1A cells compared with control shRNA-pLKO.1 (Fig. 3A). To determine whether JMJD1A knockdown inhibits cell proliferation, we performed an MTT assay to assess cell proliferation. Our results found that silencing of JMJD1A significantly decreased RV1 and LNCAP cell growth (Fig. 3B). We further found JMJD1A knockdown decreased the number and size of cell colonies (Fig. 3C) and reduced cell colony formation using a soft-agar assay (Fig. 3D). To further investigate the role of JMJD1A in prostate cancer metastasis, JMJD1A was silenced in two prostate cancer cell line models. As shown in Fig. 3E and F, we demonstrated that silencing of JMJD1A markedly abolished RV1 and LNCAP cell migration and invasion. Taken together, our studies reveal that JMJD1A downregulation decreases survival of prostate cancer cells and inhibits cell migration and invasion in vitro.

Figure 3.
Figure 3.

JMJD1A knockdown antagonizes prostate cancer cell proliferation, migration, and invasion. A, RV1 cells and LNCAP cells stably expressing control vector pLKO.1 or shJMJD1A #1 or #2 were examined by Western blot analysis. B, The viability of the above cells was detected by MTT assay. **, P < 0.01. C, The cells generated as above were seeded on 6-well plates for 2D colony formation assays. Colony numbers of cellular clones with more than 100 cells was measured (mean ± SEM of three independent experiments). The statistical analyses were performed with the ANOVA (right). **, P = 0.0013 (RV1 cells: pLK0.1 vs. shJMJD1A#1); **, P = 0.0019 (RV1 cells: pLK0.1 vs. shJMJD1A#2); **, P = 0.0048 (LNCAP cells: pLK0.1 vs. shJMJD1A#1); **, P = 0.0027 (LNCAP cells: pLK0.1 vs. shJMJD1A#2). D, The cells described above were maintained in soft agar for 3 weeks, and colony number per field was determined (mean ± SEM of three independent experiments). Quantitative analysis of colony formation is presented in the bottom panel. **, P = 0.002 (RV1 cells: pLK0.1 vs. shJMJD1A#1); **, P = 0.001 (RV1 cells: pLK0.1 vs. shJMJD1A#2); **, P = 0.0004 (LNCAP cells: pLK0.1 vs. shJMJD1A#1); **, P = 0.0041 (LNCAP cells: pLK0.1 vs. shJMJD1A#2). E and F, Both RV1 cells and LNCAP cells stably expressing vector or JMJD1A were plated in the top chamber with/without Matrigel for 24 hours. Then, cell migration (E) and invasion (F) were assessed by counting the cells that migrated to the underside of the Transwell insert. Mean ± SEM of three independent experiments.**, P = 0.0053 (E: RV1 cells pLK0.1 vs. shJMJD1A#1); **, P = 0.0013 (E: RV1 cells: pLK0.1 vs. shJMJD1A#2); **, P = 0.0019 (E: LNCAP cells PLK0.1 vs. shJMJD1A#1); **, P = 0.0007 (E: LNCAP cells pLK0.1 vs. shJMJD1A#2); **, P = 0.0001 (F: RV1 cells pLK0.1 vs. shJMJD1A#1); **, P = 0.0048 (F: RV1 cells: pLK0.1 vs. shJMJD1A#2); **, P = 0.0003 (F: LNCAP cells: pLK0.1 vs. shJMJD1A#1); **, P = 0.0004 (F: LNCAP cells: pLK0.1 vs. shJMJD1A#2).

JMJD1A impedes prostate cancer progression through a key transcription factor Snail

Many reports showed that JMJD1A affects cancer metastasis (22–26). However, the mechanism of the histone demethylase JMJD1A in prostate cancer metastasis is poorly understood. EMT is a crucial step during tumor metastasis. E-cadherin and vimentin were widely used as classic EMT markers (27–29), while Snail, ZEB1, and Slug are the key transcription factors involved in the regulation of metastasis (30–32). Here, we revealed that JMJD1A knockdown dramatically reduced Snail, but not others (Fig. 4A and B). Then, we asked whether Snail was able to reverse the effects caused by JMJD1A silencing in two prostate cancer cells. As shown in Fig. 4C–F and Supplementary Fig. S2, depletion of JMJD1A dramatically inhibited cell proliferation, migration, and invasion, while the reexpression of Snail was able to reverse the effects. Collectively, our findings show that Snail is a key downstream effector of JMJD1A in the regulation of prostate cancer progression.

Figure 4.
Figure 4.

JMJD1A promotes prostate cancer migration and invasion mainly by upregulating Snail. A and B, RV1 cells and LNCAP cells stably expressing the indicated JMJD1A shRNA, and EMT markers were examined by Western blot analysis. C and D, RV1 cells and LNCAP cells stably expressing pLKO.1 or JMJD1A shRNA together with or without Flag-Snail were subjected to Western blot analysis to examine the indicated protein levels. E and F, The cells descripted as above were plated in the top chamber with/without Matrigel for 24 hours. Then, cell migration (E) and invasion (F) were assessed by counting the cells that migrated to the underside of the Transwell insert. Mean ± SEM of three independent experiments. **, P = 0.0001 (E: column 1 vs. column 2); **, P = 0.0003 (E: column 1 vs. column 3); **, P = 0.0013 (E: column 2 vs. column 4); **, P = 0.0012 (E: column 3 vs. column 4); **, P = 0.0001 (E: column 5 vs. column 6); **, P = 0.0002 (E: column 5 vs. column 7); **, P = 0.00021 (E: column 6 vs. column 8); **, P = 0.0001 (E: column 7 vs. column 8); **, P = 0.0036 (F: column 1 vs. column 2); **, P = 0.0008 (F: column 1 vs. column 3); **, P = 0.008 (F: column 2 vs. column 4); **P, = 0.0038 (F: column 3 vs. column 4); **, P = 0.0014 (F: column 5 vs. column 6); **, P = 0.0002 (F: column 5 vs. Column 7); **, P = 0.002 (F: column 6 vs. column 8); **, P = 0.00021 (F: column 7 vs. column 8).

Snail transcription activation by JMJD1A is closely correlated with H3K9me1 and H3K36me2 demethylation

To further investigate how JMJD1A reduces Snail expression, the influence of JMJD1A on the Snail mRNA level was determined. First, we found that JMJD1A overexpression increased Snail mRNA levels (Supplementary Fig. S1A), while silencing of JMJD1A downregulated Snail mRNA levels (Supplementary Fig. S1B), suggesting that JMJD1A downregulates Snail expression at the transcriptional level, because JMJD1A demethylates H3K9me1 and H3K9Me2 to enable transcriptional activation. Thus, to reveal the mechanism of JMJD1A-mediated regulation of Snail expression, we first detected methylation levels of histone3 protein. As shown in Fig. 5A, JMJD1A overexpression significantly reduced the methylation levels of both H3K9me1 and H3K9me2. In contrast, JMJD1A dramatically increased H3K9me1 and H3K9me2 methylation (Fig. 5B). This promoted us to test whether JMJD1A protein removes mono- and di-methyl groups from H3K9 on the promoter of Snail gene and then directly interacts with Snail gene promoter. Here, using ChIP assay, we revealed that JMJD1A obviously binds to two regions of the Snail gene promoter (Fig. 5C) and the Snail gene promoter was directly immunoprecipitated with anti-JMJD1A antibodies (Fig. 5D; Supplementary Fig. S3). Moreover, ChIP assay results demonstrated that upregulated JMJD1A significantly reduced the interaction between the Snail gene promoter and both H3K9me1 and H3K9me2 compared with control (Fig. 5E and F), while JMJD1A knockdown enhanced the interaction (Fig. 5H and I). Furthermore, the luciferase reporter assays, told us that Snail gene promoter activity was activated in JMJD1A-overexpressing cells compared with control cells (Fig. 5G). In contrast, JMJD1A significantly attenuated Snail gene promoter activity compared with control cells (Fig. 5J). Taken together, these results reveal that JMJD1A transcriptionally activates Snail expression via demethylation of H3K9me1 and H3K9me2 at its special promoter region.

Figure 5.
Figure 5.

Snail transcription activation by JMJD1A is closely correlated with H3K9me1 and H3K36me2 demethylation. A, Western blot analysis detected the indicated proteins between vector and JMJD1A in RV1 cells. B, Western blot analysis measured the indicated proteins among control pLKO.1, shJMJD1A#1, and shJMJD1A#2 in RV1 cells. C, A schematic illustration of the Snail promoter regions (1–8) with (+) or without (−) binding affinity for JMJD1A (top). The arrow indicates the transcriptional start site. The result of ChIP-qPCR analysis of JMJD1A binding to the distinct regions in the Snail promoter, showing enrichment with JMJD1A antibody compared with IgG control (bottom). The p16 and GAPDH promoters as the positive and negative controls, respectively. **, P = 0.005 (column 9 vs. column 10); **, P = 0.0043 (column 11 vs. column 12). D, CHIP assay was used to examine the interaction between JMJD1A and the Snail promoter in RV1 cells. Anti-RNA Polymerase II antibody as the positive control and IgG as negative control, respectively. E and H, RV1 cells stably expressing the indicated shRNA were analyzed by ChIP assays using antibodies for control, JMJD1A, H3K9me1, or H3K9me2. The precipitated chromatin fragments were analyzed by the qPCR using the primers for the promoter region of Snail gene. Data were calculated as the percentage of input. **, P = 0.0006 (E: column 3 vs. column 4); **, P = 0.0032 (E: column 5 vs. column 6); **, P = 0.0024 (E: column 7 vs. column 8); **, P = 0.0055 (H: column 4 vs. column 5); **, P = 0.0008 (H: column 4 vs. column 6); **, P = 0.0006 (H: column 7 vs. column 8); **, P = 0.0005 (H: column 7 vs. column 9); **, P = 0.0034 (H: column 10 vs. column 11); **, P = 0.0019 (H: column 10 vs. column 12). F and I, CHIP assays were used to examine the interaction between the Snail gene promoter 4 and H3K9me1 and/or H3K9me2 among indicated cells. G and J, Dual luciferase assays were performed to detect gene luciferase activity of the Snail promoter 4 among indicated cells. **, P = 0.0003 (G: vector vs. JMJD1A); **, P = 0.00031 (J: pLK0.1 vs. shJMJD1A#1); **, P = 0.0002 (J: pLK0.1 vs. shJMJD1A#2). The results in C, E, G, H, and J are expressed as the mean ± SD of individual samples from three independent experiments. *, P < 0.05 and **, P < 0.01, Student t test.

JMJD1A overexpression promotes prostate cancer cell proliferation and tumor growth in vivo

To further test the role of JMJD1A in prostate tumor growth in vivo, control vector or stably expressed JMJD1A RV1 cells were injected into subcutaneous tissues of Balb/c nude mice. As shown in Fig. 6A, compared with control cells, JMJD1A -overexpressing cells showed an approximately 4-fold increase in tumor volume (Fig. 6A). Moreover, the weight and size of tumors formed from control cells were decreased compared with tumors formed from JMJD1A-overexpressing cells (Fig. 6B and C). We next performed IHC staining for the proliferation marker Ki67 on xenograft tumor sections. As shown in Fig. 6D, JMJD1A overexpression significantly increased the percentage of Ki67-positive cells. Consistent with our previous finding, Western blot and qRT-PCR results also demonstrated Snail expression and mRNA levels were positively associated with JMJD1A (Fig. 6E–G).

Figure 6.
Figure 6.

JMJD1A overexpression promotes prostate cancer cell survival in vivo. A–C, Tumor volumes at the indicated dates and tumor weight as well as images at day 22 (vector vs. JMJD1A) were measured. The average values are present in the bar graphs (means ± SD; n = 5 for each pair). **, P = 0.002 (A: vector vs. JMJD1A); **, P = 0.00012 (B: vector vs. JMJD1A). D, Representative IHC staining of Ki67 on tumor sections derived from the above RV1 cells. The staining was developed by DAB (brown) and counterstained by hematoxylin (blue). The image width is 0.5 mm. **, P = 0.0002. E–G, JMJD1A mRNA levels, Snail mRNA levels, and their protein levels in the above tumor sections were analyzed as the indicated methods (E: **, P = 0.00013; F: **, P = 0.0003).

Silencing of JMJD1A antagonizes tumorigenesis of prostate cancer cells

To further determine the functional roles of JMJD1A in prostate tumorigenesis, control cells (pLKO.1) or stably expressing shJMJD1A#1 RV1 cells were injected into the subcutaneous tissues of Balb/c nude mice, separately. As shown in Fig. 7A–C, JMJD1A knockdown significantly retarded tumor volume and suppressed weight and size of tumors formed from stably expressing shJMJD1A#1 cells compared with tumors formed from control cells. Next, we performed IHC staining for the proliferation marker Ki67 and the apoptosis marker cleaved caspase-3 on xenograft tumor sections. JMJD1A knockdown reduced the percentage of Ki67-positive cells (Fig. 7D), concomitant with an increased number of cells positive for active caspase-3 (Fig. 7E). Furthermore, we also examined the relationship among JMJD1A and Snail in xenograft mouse models. Western blot analysis results and qRT-PCR results revealed that JMJD1A and Snail were decreased in the shJMJD1A#1 group compared with the control group (Fig. 7F–H). In addition, we investigated whether JMJD1A knockdown decreased bone tumorigenesis of prostate cancer cells. As Supplementary Fig. S4 shows, JMJD1A knockdown inhibit bone tumorigenesis and decreased the expression of Snail. Taken together, these results suggest that JMJD1A knockdown inhibits proliferation and promotes apoptosis in xenograft prostate tumors, indicating that JMJD1A participates in prostate cancer tumorigenesis.

Figure 7.
Figure 7.

JMJD1A knockdown suppresses prostate cancer cell survival in vivo. A–C, Tumor volumes at the indicated dates and tumor weight as well as images at day 22 (pLKO.1 vs. shJMJD1A#1) were measured. The average values are present in the bar graphs (means ± SD; n = 5 for each pair). A: **, P = 0.0014; B: **, P = 0.0039. D, Representative IHC staining of Ki67 on tumor sections derived from the above RV1 cells. The staining was developed by DAB (brown) and counterstained by hematoxylin (blue). The image width is 0.5 mm. **, P = 0.0002. E, Representative IHC staining of cleaved caspase-3 and number of positive cells for active caspase-3 per field on staining described in D. **, P = 0.0012. F–H, JMJD1A mRNA levels, Snail mRNA levels, and their protein levels in the above tumor sections were analyzed as the indicated methods (F: **, P = 0.000014; G: **, P = 0.0032).

Discussion

We previously demonstrated a tumor-promoting role of JMJD1A by regulating the activities of AR and c-Myc, and thus promoting the prostate cancer cell survival (6, 11). Herein, we found that JMJD1A protein and mRNA levels were upregulated and positively associated with more advanced stage in prostate cancer specimens. In addition, a novel molecular mechanism for JMJD1A in migration and invasion was fully elucidated in prostate cancer, an AR-driven cancer, whereby JMJD1A promotes prostate cancer cell migration and invasion through enabling Snail transcriptional activation.

Many studies have reported that histone demethylases were upregulated in human cancers (11, 33–35). However, it is still unclear whether this is a bystander event or a driver of tumorigenesis. Therefore, we first evaluated the roles of JMJD1A in prostate cancer cell proliferation, migration, and invasion. Surprisingly, we found that upregulation of JMJD1A significantly promotes prostate cancer cell proliferation, migration, and invasion, while JMJD1A knockdown exhibits the opposite effects. Our mice data also support JMJD1A exerts a tumor-driver role in prostate cancer. Targeting JMJD1A might offer an effective method to control the development of prostate cancer.

The primary cause of prostate cancer–related death is metastasis such as bone metastasis and heart metastasis, which is regulated by signaling pathways such as EMT, a dynamic process that promotes cell motility with decreased adhesive ability (36–39). Lots of studies demonstrate Snail was associated with prostate cancer cell migration (40–43) and high levels of Snail expression promoted prostate cancer cell metastasis (21). However, which histone demethylase regulates Snail in prostate cancer is still unknown. Herein, we surprisingly reveal JMJD1A promotes the prostate cancer cell metastasis partially through Snail, but does not affect E-cadherin expression, consistent with a report that Snail could promote cell survival through induction of EMT or independent of EMT (13–15).

JMJD1A regulates gene expression by modulating the methylation levels of both H3K9me1 and H3K9me2 (44). Consistent with previous reports, upregulation of JMJD1A dramatically enhanced H3K9me1 and H3K9me2 methylation levels. Interestingly, our ChIP assays showed JMJD1A interacts with the Snail special promoter region. Moreover, JMJD1A overexpression or knockdown dramatically affect the interaction between the Snail promoter and both H3K9me1 and H3K9me2, suggesting that JMJD1A induces tumor progression through regulating Snail transcription. Considering we do not perform JMJD1A ChIP-seq to show whether there are binding peak(s) of JMJD1A or at least H3K9me1 and H3K9me2 at the snail gene locus, this suggests that JMJD1A may probably regulate snail expression through indirect mechanism. Importantly, compared with the control Balb/c nude mice, the growth of JMJD1A-knockdown cells was decreased along with reducing Snail expression, while JMJD1A overexpression exhibits the opposite effects. Taken together, these results suggest a tumor-promoting role of JMJD1A in tumor metastasis. Future works are needed to determine the mechanism of JMJD1A upregulation in the setting of prostate cancer progression.

To sum up, our current studies reveal a new role of JMJD1A in promoting prostate cancer cell proliferation and migration. Even the mechanism for JMJD1A upregulation in prostate cancer cells need further study, we demonstrated that JMJD1A promoted EMT, migration, and invasion through regulating Snail expression in prostate cancer cells, and thus highlight JMJD1A as potential therapeutic target for advanced prostate cancer.

Disclosure of Potential Conflicts of Interest

No potential conflicts of interest were disclosed.

Authors' Contributions

Conception and design: S.-H. Xu

Development of methodology: X.-Y. Geng, H.-W. Jiang

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): D.-E. Tang, Y. Dai, D.-X. Fu

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): D.-E. Tang, L.-L. Fan

Writing, review, and/or revision of the manuscript: S.-H. Xu

Acknowledgments

This work was supported by the National Natural Science Foundation of China (31700795), the Science & Technology Planning Project of Guangdong Province of China (2017B020209001), the Natural Science Foundation of Guangdong Province of China (2017A030310629), and the Science & Technology Plan of Shenzhen (JCYJ20170307095606266).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Footnotes

  • Note: Supplementary data for this article are available at Molecular Cancer Research Online (http://mcr.aacrjournals.org/).

  • Mol Cancer Res 2020;18:698–708

  • Received August 29, 2019.
  • Revision received November 27, 2019.
  • Accepted January 31, 2020.
  • Published first February 4, 2020.

Reference

  1. 1.↵
    1. Siegel R,
    2. Naishadham D,
    3. Jemal A
    . Cancer statistics, 2012. CA Cancer J Clin 2012;62:10–29.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Chien MH,
    2. Lin YW,
    3. Wen YC,
    4. Yang YC,
    5. Hsiao M,
    6. Chang JL,
    7. et al.
    Targeting the SPOCK1-snail/slug axis-mediated epithelial-to-mesenchymal transition by apigenin contributes to repression of prostate cancer metastasis. J Exp Clin Cancer Res 2019;38:246.
    OpenUrl
  3. 3.↵
    1. Qin Y,
    2. Dey A,
    3. Purayil HT,
    4. Daaka Y
    . Maintenance of androgen receptor inactivation by S-nitrosylation. Cancer Res 2013;73:6690–9.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Giacinti S,
    2. Bassanelli M,
    3. Aschelter AM,
    4. Milano A,
    5. Roberto M,
    6. Marchetti P
    . Resistance to abiraterone in castration-resistant prostate cancer: a review of the literature. Anticancer Res 2014;34:6265–9.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Vander Ark A,
    2. Cao J,
    3. Li X
    . Mechanisms and approaches for overcoming enzalutamide resistance in prostate cancer. Front Oncol 2018;8:180.
    OpenUrl
  6. 6.↵
    1. Fan L,
    2. Zhang F,
    3. Xu S,
    4. Cui X,
    5. Hussain A,
    6. Fazli L,
    7. et al.
    Histone demethylase JMJD1A promotes alternative splicing of AR variant 7 (AR-V7) in prostate cancer cells. Proc Natl Acad Sci U S A 2018;115:E4584–93.
    OpenUrlAbstract/FREE Full Text
  7. 7.↵
    1. Okada Y,
    2. Scott G,
    3. Ray MK,
    4. Mishina Y,
    5. Zhang Y
    . Histone demethylase JHDM2A is critical for Tnp1 and Prm1 transcription and spermatogenesis. Nature 2007;450:119–23.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Kuroki S,
    2. Matoba S,
    3. Akiyoshi M,
    4. Matsumura Y,
    5. Miyachi H,
    6. Mise N,
    7. et al.
    Epigenetic regulation of mouse sex determination by the histone demethylase Jmjd1a. Science 2013;341:1106–9.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Loh YH,
    2. Zhang W,
    3. Chen X,
    4. George J,
    5. Ng HH
    . Jmjd1a and Jmjd2c histone H3 Lys 9 demethylases regulate self-renewal in embryonic stem cells. Genes Dev 2007;21:2545–57.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    1. Kaukonen R,
    2. Mai A,
    3. Georgiadou M,
    4. Saari M,
    5. De Franceschi N,
    6. Betz T,
    7. et al.
    Normal stroma suppresses cancer cell proliferation via mechanosensitive regulation of JMJD1a-mediated transcription. Nat Commun 2016;7:12237.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Fan L,
    2. Peng G,
    3. Sahgal N,
    4. Fazli L,
    5. Gleave M,
    6. Zhang Y,
    7. et al.
    Regulation of c-Myc expression by the histone demethylase JMJD1A is essential for prostate cancer cell growth and survival. Oncogene 2016;35:2441–52.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Hay ED
    . An overview of epithelio-mesenchymal transformation. Acta Anat 1995;154:8–20.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Barrallo-Gimeno A,
    2. Nieto MA
    . The Snail genes as inducers of cell movement and survival: implications in development and cancer. Development 2005;132:3151–61.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Martinez-Alvarez C,
    2. Blanco MJ,
    3. Perez R,
    4. Rabadan MA,
    5. Aparicio M,
    6. Resel E,
    7. et al.
    Snail family members and cell survival in physiological and pathological cleft palates. Dev Biol 2004;265:207–18.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Emadi Baygi M,
    2. Soheili ZS,
    3. Schmitz I,
    4. Sameie S,
    5. Schulz WA
    . Snail regulates cell survival and inhibits cellular senescence in human metastatic prostate cancer cell lines. Cell Biol Toxicol 2010;26:553–67.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Moody SE,
    2. Perez D,
    3. Pan TC,
    4. Sarkisian CJ,
    5. Portocarrero CP,
    6. Sterner CJ,
    7. et al.
    The transcriptional repressor Snail promotes mammary tumor recurrence. Cancer Cell 2005;8:197–209.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Tuhkanen H,
    2. Soini Y,
    3. Kosma VM,
    4. Anttila M,
    5. Sironen R,
    6. Hamalainen K,
    7. et al.
    Nuclear expression of Snail1 in borderline and malignant epithelial ovarian tumours is associated with tumour progression. BMC Cancer 2009;9:289.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Franci C,
    2. Gallen M,
    3. Alameda F,
    4. Baro T,
    5. Iglesias M,
    6. Virtanen I,
    7. et al.
    Snail1 protein in the stroma as a new putative prognosis marker for colon tumours. PLoS One 2009;4:e5595.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Dhanasekaran SM,
    2. Barrette TR,
    3. Ghosh D,
    4. Shah R,
    5. Varambally S,
    6. Kurachi K,
    7. et al.
    Delineation of prognostic biomarkers in prostate cancer. Nature 2001;412:822–6.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Beach S,
    2. Tang H,
    3. Park S,
    4. Dhillon AS,
    5. Keller ET,
    6. Kolch W,
    7. et al.
    Snail is a repressor of RKIP transcription in metastatic prostate cancer cells. Oncogene 2008;27:2243–8.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Odero-Marah VA,
    2. Wang R,
    3. Chu G,
    4. Zayzafoon M,
    5. Xu J,
    6. Shi C,
    7. et al.
    Receptor activator of NF-kappaB Ligand (RANKL) expression is associated with epithelial to mesenchymal transition in human prostate cancer cells. Cell Res 2008;18:858–70.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Peng K,
    2. Su G,
    3. Ji J,
    4. Yang X,
    5. Miao M,
    6. Mo P,
    7. et al.
    Histone demethylase JMJD1A promotes colorectal cancer growth and metastasis by enhancing Wnt/beta-catenin signaling. J Biol Chem 2018;293:10606–19.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    1. Sechler M,
    2. Parrish JK,
    3. Birks DK,
    4. Jedlicka P
    . The histone demethylase KDM3A, and its downstream target MCAM, promote Ewing sarcoma cell migration and metastasis. Oncogene 2017;36:4150–60.
    OpenUrl
  24. 24.↵
    1. Osawa T,
    2. Tsuchida R,
    3. Muramatsu M,
    4. Shimamura T,
    5. Wang F,
    6. Suehiro J,
    7. et al.
    Inhibition of histone demethylase JMJD1A improves anti-angiogenic therapy and reduces tumor-associated macrophages. Cancer Res 2013;73:3019–28.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    1. Zhan M,
    2. Wen F,
    3. Liu L,
    4. Chen Z,
    5. Wei H,
    6. Zhou H
    . JMJD1A promotes tumorigenesis and forms a feedback loop with EZH2/let-7c in NSCLC cells. Tumour Biol 2016;37:11237–47.
    OpenUrl
  26. 26.↵
    1. Tee AE,
    2. Ling D,
    3. Nelson C,
    4. Atmadibrata B,
    5. Dinger ME,
    6. Xu N,
    7. et al.
    The histone demethylase JMJD1A induces cell migration and invasion by up-regulating the expression of the long noncoding RNA MALAT1. Oncotarget 2014;5:1793–804.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Dong C,
    2. Yuan T,
    3. Wu Y,
    4. Wang Y,
    5. Fan TW,
    6. Miriyala S,
    7. et al.
    Loss of FBP1 by Snail-mediated repression provides metabolic advantages in basal-like breast cancer. Cancer Cell 2013;23:316–31.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Pastushenko I,
    2. Brisebarre A,
    3. Sifrim A,
    4. Fioramonti M,
    5. Revenco T,
    6. Boumahdi S,
    7. et al.
    Identification of the tumour transition states occurring during EMT. Nature 2018;556:463–8.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Ye X,
    2. Brabletz T,
    3. Kang Y,
    4. Longmore GD,
    5. Nieto MA,
    6. Stanger BZ,
    7. et al.
    Upholding a role for EMT in breast cancer metastasis. Nature 2017;547:E1–3.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Persa OD,
    2. Niessen CM
    . Epithelial polarity limits EMT. Nat Cell Biol 2019;21:299–300.
    OpenUrl
  31. 31.↵
    1. Zheng X,
    2. Carstens JL,
    3. Kim J,
    4. Scheible M,
    5. Kaye J,
    6. Sugimoto H,
    7. et al.
    Epithelial-to-mesenchymal transition is dispensable for metastasis but induces chemoresistance in pancreatic cancer. Nature 2015;527:525–30.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Inoue A,
    2. Seidel MG,
    3. Wu W,
    4. Kamizono S,
    5. Ferrando AA,
    6. Bronson RT,
    7. et al.
    Slug, a highly conserved zinc finger transcriptional repressor, protects hematopoietic progenitor cells from radiation-induced apoptosis in vivo. Cancer Cell 2002;2:279–88.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Wang Z,
    2. Yang X,
    3. Liu C,
    4. Li X,
    5. Zhang B,
    6. Wang B,
    7. et al.
    Acetylation of PHF5A modulates stress responses and colorectal carcinogenesis through alternative splicing-mediated upregulation of KDM3A. Mol Cell 2019;74:1250–63.
    OpenUrl
  34. 34.↵
    1. Kawamura A,
    2. Munzel M,
    3. Kojima T,
    4. Yapp C,
    5. Bhushan B,
    6. Goto Y,
    7. et al.
    Highly selective inhibition of histone demethylases by de novo macrocyclic peptides. Nat Commun 2017;8:14773.
    OpenUrl
  35. 35.↵
    1. Kim TD,
    2. Jin F,
    3. Shin S,
    4. Oh S,
    5. Lightfoot SA,
    6. Grande JP,
    7. et al.
    Histone demethylase JMJD2A drives prostate tumorigenesis through transcription factor ETV1. J Clin Invest 2016;126:706–20.
    OpenUrl
  36. 36.↵
    1. Liu X,
    2. Chen L,
    3. Fan Y,
    4. Hong Y,
    5. Yang X,
    6. Li Y,
    7. et al.
    IFITM3 promotes bone metastasis of prostate cancer cells by mediating activation of the TGF-beta signaling pathway. Cell Death Dis 2019;10:517.
    OpenUrl
  37. 37.↵
    1. Jiang S,
    2. Wang X,
    3. Song D,
    4. Liu X,
    5. Gu Y,
    6. Xu Z,
    7. et al.
    Cholesterol induces epithelial-to-mesenchymal transition of prostate cancer cells by suppressing degradation of EGFR through APMAP. Cancer Res 2019;79:3063–75.
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    1. Stylianou N,
    2. Lehman ML,
    3. Wang C,
    4. Fard AT,
    5. Rockstroh A,
    6. Fazli L,
    7. et al.
    A molecular portrait of epithelial-mesenchymal plasticity in prostate cancer associated with clinical outcome. Oncogene 2019;38:913–34.
    OpenUrlCrossRef
  39. 39.↵
    1. Grant CM,
    2. Kyprianou N
    . Epithelial mesenchymal transition (EMT) in prostate growth and tumor progression. Transl Androl Urol 2013;2:202–11.
    OpenUrl
  40. 40.↵
    1. Zhau HE,
    2. Odero-Marah V,
    3. Lue HW,
    4. Nomura T,
    5. Wang R,
    6. Chu G,
    7. et al.
    Epithelial to mesenchymal transition (EMT) in human prostate cancer: lessons learned from ARCaP model. Clin Exp Metastasis 2008;25:601–10.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Gan Y,
    2. Shi C,
    3. Inge L,
    4. Hibner M,
    5. Balducci J,
    6. Huang Y
    . Differential roles of ERK and Akt pathways in regulation of EGFR-mediated signaling and motility in prostate cancer cells. Oncogene 2010;29:4947–58.
    OpenUrlCrossRefPubMed
  42. 42.↵
    1. Mak P,
    2. Leav I,
    3. Pursell B,
    4. Bae D,
    5. Yang X,
    6. Taglienti CA,
    7. et al.
    ERbeta impedes prostate cancer EMT by destabilizing HIF-1alpha and inhibiting VEGF-mediated snail nuclear localization: implications for Gleason grading. Cancer Cell 2010;17:319–32.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Bokobza SM,
    2. Ye L,
    3. Kynaston H,
    4. Jiang WG
    . Growth and differentiation factor 9 (GDF-9) induces epithelial-mesenchymal transition in prostate cancer cells. Mol Cell Biochem 2011;349:33–40.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Yamane K,
    2. Toumazou C,
    3. Tsukada Y,
    4. Erdjument-Bromage H,
    5. Tempst P,
    6. Wong J,
    7. et al.
    JHDM2A, a JmjC-containing H3K9 demethylase, facilitates transcription activation by androgen receptor. Cell 2006;125:483–95.
    OpenUrlCrossRefPubMed
Back to top
View this article with LENS

Open full page PDF
Article Alerts
Email Article
Citation Tools
Histone Demethylase JMJD1A Promotes Tumor Progression via Activating Snail in Prostate Cancer
Dong-E Tang, Yong Dai, Ling-Ling Fan, Xin-Yan Geng, De-Xue Fu, Hao-Wu Jiang and Song-Hui Xu
Mol Cancer Res May 1 2020 (18) (5) 698-708; DOI: 10.1158/1541-7786.MCR-19-0889
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo

Jump to section

Advertisement