ADAMTS9-AS1 Constrains Breast Cancer Cell Invasion and Proliferation via Sequestering miR-301b-3p

Introduction

As the most prevalent female malignant tumor globally, morbidity of breast cancer has been steadily increasing since 1987, and its incidence accounts for approximately 24.5% of all women’s malignancies (Sung et al., 2021). Not a single-gene disorder, breast cancer is a diversified disease classified by histology, immunopathology, mRNA expression profile and miRNA expression features rather than a monogenic disease (Laulin and Brudieux, 1990). As molecular biological research deepens, abnormal signaling pathways of tumor inhibitor genes and oncogenes manifest connection with breast cancer development (Boo et al., 2016; Gilam et al., 2016; Nassar et al., 2017; Sun et al., 2017). So far, the specific molecular mechanism for breast cancer progression has not been completely clarified.

Long non-coding RNAs (lncRNAs) exert a role in expression modulation at genetic levels and also exert functions on tumor-relevant cancer cell behaviors (Mercer et al., 2009; Fatica and Bozzoni, 2014). More and more evidence demonstrated association lncRNA-miRNA-mRNA axis maintains in tumor development. As a critical post-transcriptional regulator, the activity of miRNA can be sequestered by lncRNA (Wang et al., 2018; Peng et al., 2019; Shen et al., 2019). Such lncRNA is called competitive endogenous RNA (ceRNA). ADAMTS9-AS1 is an antisense lncRNA generated by reverse transcription of protein-coding gene ADAMTS9 and closely linked with various tumors. ADAMTS9 is a tumor inhibitor in assorted cancer cells: esophagus carcinoma (Lo et al., 2007), gastric cancer (Du et al., 2013), breast cancer (Ocak et al., 2013) and nasopharyngeal carcinoma (Lung et al., 2008; Sheu et al., 2009). Nevertheless, there is a shortage of deep research about the function of lncRNA ADAMTS9-AS1 in tumors. ADAMTS9-AS1 was found in esophageal squamous cell carcinoma at first (Li et al., 2017). It is greatly important for ectoderm and epithelial cells, as well as a novel biomarker for esophageal squamous cell carcinoma (Li et al., 2017). A study comprehensively analyzed lncRNA and found that ADAMTS9-AS1 is central in differential co-expression network in ovarian cancer (Wang et al., 2016). Few references have deeply researched ADAMTS9-AS1 in breast cancer yet. Here, such work was performed, which might conduce to biological understanding of this lethal disease.

Here, on the basis of the downloaded bioinformatics data, LncRNA ADAMTS9-AS1 was determined using differential analysis as well as survival analysis. Its role in breast cancer was explored by cell functional experiments. In an effort to scrutinize its downstream mechanism, bioinformatics databases were applied to predict differential miRNA and its target gene that had binding sites with ADAMTS9-AS1. The evidence proved that ADAMTS9-AS1 restrained invasion and proliferation of cancer cells through sequestering miR-301b-3p to target TGFBR2 and regulate JAK STAT signaling pathway. The results provide a foundation for further revealing of the mechanism of ADAMTS9-AS1 in breast cancer.

Materials and Methods

Bioinformatics Analysis

The Cancer Genome Atlas (TCGA)¹ provided expression data of lncRNAs, mRNAs and miRNAs related to breast cancer. Clinical information was offered in Supplementary Table 1. Differential expression analysis was performed with edgeR (| logFC| > 2, padj < 0.05). Differentially expressed lncRNAs (DElncRNAs) and miRNAs (DEmiRNAs) that had binding sites were screened through miRcode.² Downstream regulatory target genes of DEmiRNAs were predicted by miRDB,³ TargetScan,⁴ along with miRTarBase.⁵ The predicted mRNAs were intersected with DEmRNAs to obtain target DEmRNAs of DEmiRNAs. Pathway enrichment analysis of unigenes was conducted by using Gene Set Enrichment Analysis (GSEA).

Cell Culture and Transfection

Mammary epithelial cells MCF10A (CBP60419) and breast cancer cells MDA-MB-231 (CBP60382), BT-549 (CBP60357), MDA-MB-468 (CBP60387) and MCF-7 (CBP60380) were cultured in complete mediums. All mediums and cell lines were purchased from Nanjing Cobioer Biotech. Cells were inoculated in 6-well plates at 3 × 10⁵ each well until cell growth density reached 70–90%. Then, cells were subjected to transfection utilizing Lipofectamine 2000 (11668-019, Invitrogen, California, United States). Next, 250 μl Opti-MEM (11668-019, Invitrogen, California, United States) medium was used to release 4 μg target plasmids and 10 μl Lipofectamine 2000, respectively, and mixed well by flicking. After being maintained at room temperature for 5 min, two liquids were mixed evenly, and added into cell culture wells after 20 min. Culture plates were shaken gently and mixed evenly. Then plates were generally cultured in an incubator. After 24 h, medium was changed and cells were collected 36–48 h later.

The si-ADAMTS9-AS1 and si-TGFBR2 were chemically synthesized by Genepharma (Shanghai, China). The siRNA transfection was undertaken with Lipofectamine 2000 (Invitrogen) as the standard plan. The oe-negative control (NC) (NC plasmids of overexpressing ADAMTS9-AS1, the vector used was pcDNA™3.1(+) vector (V79020), Invitrogen), oe-ADAMTS9-AS1 (overexpressing ADAMTS9-AS1 plasmids, the vector used was pcDNA™3.1(+) vector), inhibitor NC (NC of miR-301b-3p inhibitor) and miR-301b-3p inhibitor were procured from Genepharma (Shanghai, China).

Quantitative Reverse Transcription-Polymerase Chain Reaction

Total RNA was detached from cells (1 × 10⁶ per well) with Trizol (No. 16096020, Thermo Fisher Scientific, New York, United States). Next, 5 μg samples were chosen for reverse transcription to synthesis into cDNA following the instruction of cDNA kit (K1622; Fermentas Inc., Ontario, CA, United States). MiR-301b-3p expression was detected by qPCR on TaqMan MicroRNA Assay with cDNA as the template. The reaction of qRT-PCR was as follows: 95^°C for 2 min, then 45 cycles of 95^°C for 15 s, 60^°C for 45 s, 72^°C for 45 s. U6 was taken as the internal reference to normalize results.

qRT-PCR was performed on ADAMTS9-AS1 and TGFBR2 according to the instruction of TaqMan Gene Expression Assays protocol (Applied Biosystems, Foster City, CA, United States). GAPDH was taken as the internal reference. PCR procedures were designed as follows: 95^°C for 10 min, 35 cycle of 95^°C for 15 s, 60^°C for 30 s, 72^°C for 45 s. All qRT-PCR was set with 3 parallel wells. Primers were designed as shown in Table 1 and synthesized by TAKARA Company (Beijing, China). Relative transcription level of target gene was formulated by relative quantitative method (2^–△△CT). △△Ct = △Ct (experimental group)−△Ct (control group), △Ct = Ct (target gene)−t (internal reference), relative transcription level of target mRNA = 2^{–△△ Ct}.

TABLE 1

Table 1. qRT-RCR primer sequences.

Western Blot

Total proteins were acquired from cells (1 × 10⁶ per well), and protein concentration was determined by bicinchoninic acid (BCA) kit (Thermo, United States). Next, 30 μg total proteins were loaded on sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) at constant 80 V for 35 min and then 120 V for 45 min. The proteins were transferred onto a polyvinylidene fluoride (PVDF) membrane (Amersham, United States). The membrane was blocked with 5% skim milk at normal temperature for 1 h. It was followed by an overnight incubation at 4°C overnight with rabbit anti-TGFBR2, Ki67, PCNA, MMP-2, MMP-9, JAK1, p-JAK1, JAK2, p-JAK2, STAT1, p-STAT1, STAT3, p-STAT3, β-actin. All antibodies were offered by Abcam (Cambridge, United Kingdom). The membrane was washed with Tris-buffered saline with 0.1% Tween^®20 detergent (TBST) three times (10 min each time). Afterward, secondary antibody goat anti-rabbit IgG H&L labeled by horseradish peroxidase (HRP) was added for 1 h of incubation at room temperature. Phosphate buffered saline with Tween 20 (PBST) solution was used for 3 washes of the membrane (each time 10 min). Chemiluminescence apparatus (GE, United States) was applied for scanning and development. Protein bands were gray-scale scanned with Image Pro Plus 6.0 (Media Cybernetics, United States) software to analyze relative protein expression.

In vitro Colony Formation Assay

The transfected breast cancer cells (1 × 10⁴ per well) were inoculated into 12-well plates and cultured at 37^°C with 5% CO₂ for 2 weeks. Then cells were fixed with 4% paraformaldehyde and stained with 0.1% crystal violet. The number of colonies, whose cells are more than 50, was manually counted.

EdU Assay

Cells (3 × 10⁵ per well) to be tested were seeded in 24-well plates with 3 parallel wells. EdU (C10310-1, RiboBio, Guangzhou, China) was added into the culture solution to achieve a concentration of 10 μmol/l for 2 h of incubation in the incubator. After removing medium, cells undergone 15 min of fixation with PBS + 4% paraformaldehyde at common temperature for and twice washes with PBS + 3% bovine serum albumin (BSA). PBS containing 0.5% Triton-100 was added for 20 min of incubation at room temperature and then PBS with 3% BSA was used to wash cells twice. Each well was added with 100 μl Apollo^® 567 (C10310-1, RiboBio, Guangzhou, China) for 30 min of incubation in the dark at room temperature. Thereafter, cells were washed with PBS containing 3% BSA. DIPA (C0065, Solarbio, Beijing, China) was added to stain nucleus for 5 min. Next, cells were washed with PBS in triplicate and sealed. Cells were observed under a fluorescent microscope (Type: FM-600, Shanghai PuDa Company) and positive cell number in each field was recorded. Under the microscope, total cells presented as blue and positive cells presented as red. Percentage of EdU stained positive cells was calculated in 3 random fields of each well. EdU positive rate = (EdU positive nucleus number/total nucleus) × 100%.

Transwell Assay

Extracellular matrix (ECM) gel was maintained at 4°C overnight and then attenuated with serum-free medium at 1:9 ratio to 1 mg/ml. Each polycarbonate membrane of the 24-well Transwell chamber was added with 40 μl ECM gel and maintained generally in the incubator for 5 h. Polymerized gel was prepared by ECM gel and redundant solution was removed. Then, 70 μl/chamber pure Dulbecco’s Modified Eagle Medium (DMEM) was added. Cells (1 × 10⁴ per well) were incubated in an incubator at 37^°C for 0.5 h to rehydrate Matrigel and redundant culture solution was removed. Cells in each group were starved for 24 h without serum, digested, centrifuged and resuspended with DMEM without FBS to a final concentration of 2.5 × 10⁵ cells/ml. The upper chamber where the basement membrane was hydrated was added with added 0.2 ml suspension, while the lower chamber was added with 700 μl pre-cooled DMEM containing 10% FBS. Then, plates were generally cultured in a moist incubator for 24 h. Thereafter, the chamber was discarded, and cells in the chamber and on the membrane were wiped by a wet cotton swab. The remaining cells were fixed with methyl alcohol for 30 min, stained with 0.1% crystal violet for 20 min and dried in inverted position. Cells were observed and photographed under an inverted microscope. Average number of cells that passed through the membrane was counted in 5 random fields.

RNA Binding Protein Immunoprecipitation

Co-transfection was, respectively, undertaken on MCD-7 cells (3 × 10⁵ per well) with mimic NC/miR-301b-3p mimic along with oe-NC/oe-ADAMTS9-AS1. The binding between miR-301b-3p with ADAMTS9-AS1 and TGFBR2 was examined using RIP kit. MCF-7 cells were washed with pre-cooled PBS and then the supernatant was removed. Cells were lysed with equivalent volume of lysis solution (P0013B, Beyotime) in ice-bath for 5 min and centrifuged at 4^°C and 14,000 rpm for 10 min. The cell extract was incubated with antibodies for co-precipitation. Specific steps were: 50 μl magnetic beads were taken from each co-precipitation reaction system, washed and resuspended in 100 μl RIP Wash Buffer. According to experimental groups, 5 μg antibodies were added for binding. Magnetic beads-antibody complex was resuspended in 900 μl RIP Wash Buffer after being washed. Thereafter, 100 μl cell extract was added for an overnight incubation at 4^°C. Samples were placed on the magnetic drill to collect magnetic beads-protein complex. The antibodies used in RIP were: Ago2 (ab32381, 1:50, Abcam, United Kingdom) which was evenly shaken at room temperature for 30 min, and IgG (1:100, ab172730, Abcam, United Kingdom) as the NC. Samples were digested by protease K. Then, RNA was extracted from the samples for succeeding analyses.

Dual-Luciferase Reporter Gene Assay

The whole length of 3′UTR of TGFBR2 gene was amplified. PCR products were cloned into multiple cloning sites on the downstream of pmirGLO luciferase gene (Promega, WI, United States) by using restriction enzyme sites, to construct TGFBR2-wild type (WT) (CCAGCUAUGA CCACAUUGCACUU) group. Binding sites of miR-301b-3p and its target gene were predicted by target gene databases. TGFBR2-mutant (MUT) vector (CCAGCUAUGACCACAUAC GUGAU) was constructed by PCR site-directed mutation method. Renilla luciferase expression vector pRL-TK (TaKaRa, Dalian, China) was taken as the internal reference. HEK-293T cells (3 × 10⁵ per well) were co-transfected with mimic NC/miR-301b-3p mimic and luciferase reporter vectors, respectively. Dual Luciferase activity was tested by using Dual-Luciferase Reporter Assay System (Promega, Madison, WI, United States).

Statistical Analysis

Results were treated with SPSS 21.0 statistical software (SPSS, Inc., Chicago, IL, United States). Three times of biological repetitions and technical repetitions were required for each assay. Measurement data were shown as mean ± standard deviation. Comparison of two or more groups were, respectively, Student’s t-test and analysis of variance (ANOVA). Repeated-measures ANOVA was adopted to compare results at each time point. Patient’s overall survival curve was calculated by Kaplan-Meier. Patient’s survival differences were analyzed by log-rank. Analysis method of enumeration data was chi-square test. P < 0.05 stands for statistical significance.

Results

Inactivated ADAMTS9-AS1 in Breast Cancer

There were 1,017 DElncRNAs, 74 DEmiRNAs and 2,161 DEmRNAs, as offered by Figure 1A and detailed in Supplementary Tables 2–4. miRcode database was used to compare the results, finding that 60 lncRNAs and 18 miRNAs had binding sites (Supplementary Table 5). Low and high ADAMTS9-AS1 groups were classified by median ADAMTS9-AS1 level. Survival analysis of the 60 lncRNAs revealed notably lower survival of low ADAMTS9-AS1 expression group. Furthermore, ADAMTS9-AS1 levels noticeably pertained to clinical M stage, and patients in high ADAMTS9-AS1 group were mostly in M0 stage (Figure 1B). Hence the downregulation of ADAMTS9-AS1 remarkably impacted patient’s prognostic conditions. Meanwhile, qRT-PCR uncovered that level of ADAMTS9-AS1 was prominently low in cancer cell lines. MDA-MB-231 cells (p < 0.05) and MCF-7 cells (p < 0.05) showed the lowest and highest ADAMTS9-AS1 levels, respectively (Figure 1C).

FIGURE 1

Figure 1. ADAMTS9-AS1 is lowly expressed in breast cancer. (A) Volcano map of DElncRNAs, DEmiRNAs and DEmRNAs in breast cancer tumor group and normal group. (B) The survival curve of ADAMTS9-AS1 expression level on the prognosis of patients, with X-axis representing time (in years) and Y-axis representing survival. Red: high expression, blue: low expression. (C) The level of ADAMTS9-AS1 in breast cancer cells; *p < 0.05.

ADAMTS9-AS1 Restrains Breast Cancer Cell Proliferation and Invasion

After transfection of plasmids oe-NC and oe-ADAMTS9-AS1 was done, the latter exhibited a markedly increase in ADAMTS9-AS1 level, with respect to the former group (Figure 2A). Transfected cells were therefore employed in the succeeding analyses. Cell proliferation was examined through colony formation assay and EdU assay (Figures 2B,C). The results revealed that the cell proliferation was remarkably reduced by overexpressed ADAMTS9-AS1. Meanwhile, the result of Transwell assay uncovered that the overexpressed ADAMTS9-AS1 notably constrained the invasive ability of cells (Figure 2D). Western blot result revealed that oe-ADAMTS9-AS1 markedly decreased the expression of tumor markers for proliferation and invasion, such as Ki67, PCNA, MMP-2, MMP-9 (Figure 2E), suggesting that oe-ADAMTS9-AS1 could reduce tumor proliferative, and invasive abilities.

FIGURE 2

Figure 2. ADAMTS9-AS1 restrains cancer cell invasion and proliferation. ADAMTS9-AS1 was overexpressed in MDA-MB-231 and BT-549, while being silenced in MCF7. (A) Level of ADAMTS9-AS1 in breast cancer cells. (B) Results of in vitro colony formation method. (C) Results of EdU method (400×). (D) Cell invasion results (100×). (E) The expression of Ki67, PCNA, MMP-2 and MMP-9 after interfering with ADAMTS9-AS1 expression. (F) ADAMTS9-AS1 level in cancer cells. (G) Proliferation of cancer cells. (H) Results of EdU method (400×). (I) Transwell method is used to test cell invasion (100×). (J) The expression of Ki67, PCNA, MMP-2 and MMP-9. *p < 0.05.

ADAMTS9-AS1 was silenced in MCF-7 and the expression was examined afterward, finding a prominent reduction in si-ADAMTS9-AS1-1/si-ADAMTS9-AS1-2 groups (Figure 2F). The si-ADAMTS9-AS1-1 with the highest interference efficiency was employed hereinafter and was uniformly named as si-ADAMTS9-AS1. Further, the changes of breast cancer proliferative ability (Figures 2G,H) and invasive ability (Figure 2I) were examined by multiple methods. The results uncovered that cancer cell proliferation and invasion were prominently boosted after the downregulation of ADAMTS9-AS1. Western blot result uncovered that levels of Ki67, MMP-9, MMP-2 and PCNA were notably boosted after the downregulation of ADAMTS9-AS1 (Figure 2J). All the above, ADAMTS9-AS1 repressed the proliferation and invasion of breast cancer cells. Morphology of cells after ADAMTS9-AS1 silence or overexpression manifested no noticeable differences between groups (Supplementary Figure 2).

ADAMTS9-AS1 Traps miR-301b-3p to Regulate TGFBR2

ADAMTS9-AS1 targeted 6 DEmiRNAs (miR-301b-3p, miR-96-3p, miR-144-3p, miR-145-3p, miR-182-3p, and miR-21-3p), among which miR-301b-3p markedly pertained to patient’s prognosis (Figure 3A). Median miR-301b-3p was utilized to classify groups. As displayed in Supplementary Figure 1, miR-301b-3p levels were irrelevant to TNM stage. 8 DEmRNAs binding with miR-301b-3p were discovered in the overlap of target genes and DEmRNAs (Figure 3B). These 8 genes were subjected to Pearson correlation. Pearson correlation coefficient between TGFBR2 and ADAMTS9-AS1 was the highest (Figure 3C). Median TGFBR2 level was utilized for grouping, followed by survival analysis and expression differences. Level of TGFBR2 was of notably difference in varying pathological stages (Figure 3D). Moreover, TGFBR2 participates in the proliferation and invasion of tumor cells (Liu et al., 2018; Zhou B. et al., 2018). Therefore, we speculated that ADAMTS9-AS1 mediated TGFBR2 levels via sequestering miR-301b-3p, thus regulating cancer proliferation and invasion.

FIGURE 3

Figure 3. ADAMTS9-AS1 competitively binds to miR-301b-3p to promote TGFBR2. (A) Survival curve of miR-301b-3p level and prognosis of patients. (B) Venn plot of predicted mRNAs and DEmRNAs; diffSig denotes differentially expressed genes screened from TCGA database. (C) Association between ADAMTS9-AS1 and 8 target mRNAs. (D) The survival curve of TGFBR2 expression level on the prognosis of patients. (E,F) Results of RIP method. (G) Result of dual-luciferase method. (H) The expression of ADAMTS9-AS1, miR-301b-3p and TGFBR2. *p < 0.05.

Further mechanism verification was conducted. RIP assay exhibited stimulated enrichment of ADAMTS9-AS1/TGFBR2 upon miR-301b-3p mimic treatment, implying that miR-301b-3p bound to ADAMTS9-AS1 and TGFBR2, respectively (p < 0.05) (Figure 3E). Moreover, similar trend was observed concerning enrichment of miR-301b-3p/TGFBR2 in MCF-7 cell with oe-ADAMTS9-AS1, implying that ADAMTS9-AS1/TGFBR2 competitively bound miR-301b-3p (Figure 3F). Dual-luciferase result uncovered that miR-301b-3p mimic markedly repressed luciferase intensity of TGFBR2-WT (p < 0.05) but did not impact TGFBR2-MUT (Figure 3G). TGFBR2 expression was prominently decreased after ADAMTS9-AS1 was interfered (p < 0.05), while miR-301b-3p inhibitor could reverse TGFBR2 expression (p < 0.05) (Figure 3H). Altogether, ADAMTS9-AS1 competitively bound miR-301b-3p, so as to promote the expression of TGFBR2.

ADAMTS9-AS1 Traps miR-301b-3p to Regulate TGFBR2, Thereby Repressing Breast Cancer Cell Proliferation and Invasion

We first divided MCF-7 cells into varying transfecting groups. Further, changes in cell proliferation (Figures 4A,B) and invasion (Figure 4C) of breast cancer cells were detected by in vitro colony formation, EdU and Transwell assays, respectively. The results revealed that the invasion (p < 0.05) as well as proliferation (p < 0.05) of breast cancer were markedly boosted after the downregulation of ADAMTS9-AS1 expression, while these abilities reversed upon co-transfection with si-ADAMTS9-AS1 and miR-301b-3p inhibitor. Western blot result uncovered that Ki67, MMP-9, MMP-2, PCNA and other tumor markers were prominently activated after the downregulation of ADAMTS9-AS1, indicating increased tumor proliferative activity. However, in si-ADAMTS9-AS1 + miR-301b-3p inhibitor group, the expression of those tumor makers was reversed (p < 0.05) (Figure 4D), indicating a reduced tumor proliferative activity. To sum up, downregulated ADAMTS9-AS1 facilitated cancer cell invasion and proliferation through miR-301b-3p.

FIGURE 4

Figure 4. ADAMTS9-AS1 promotes TGFBR2 through miR-301b-3p to inhibit breast cancer cell invasion and proliferation. (A) The proliferation of MCF-7 cells. (B) EdU assay used to detect MCF-7 cell proliferation (400×). (C) Invasion of MCF-7 cells (100×). (D) Levels of Ki67, PCNA, MMP-2 and MMP-9. (E) Level of TGFBR2 in si-NC, si-TGFBR2-1 and si-TGFBR2-2 groups. The silencing sequence with the highest silencing efficiency was named as si-TGFBR2. (F) ADAMTS9-AS1 and TGFBR2 expression in MDA-MB-231 cells. (G) Proliferation of MDA-MB-231. (H) Proliferation of MDA-MB-231 (400×). (I) Cancer cell invasion (100×). (J) Levels of Ki67, PCNA, MMP-2 and MMP-9. *p < 0.05.

Two TGFBR2 interference sequences were constructed, and TGFBR2 expression was examined. The silencing sequence with the highest silencing efficiency, si-TGFBR2-2 (p < 0.05), was selected and named as si-TGFBR2 (Figure 4E). Three transfecting groups were generated. qRT-PCR indicated that TGFBR2 level was notably boosted after the overexpression of ADAMTS9-AS1 (p < 0.05). Relative to the oe-ADAMTS9-AS1 + si-NC group, the level of ADAMTS9-AS1 in the oe-ADAMTS9-AS1 + si-TGFBR2 group had no notable change (p < 0.05) while TGFBR2 was prominently constrained (p < 0.05) (Figure 4F). Further, cell proliferation (Figures 4G,H) and invasion (Figure 4I) were detected to be markedly reduced upon ADAMTS9-AS1 overexpression(p < 0.05), while overexpressed ADAMTS9-AS1 and silencing TGFBR2 together could reverse that effect (p < 0.05). Western blot result uncovered that Ki67, PCNA, MMP-2, and MMP-9 were notably repressed upon ADAMTS9-AS1 overexpression (p < 0.05), indicating that the tumor proliferation was decreased. However, overexpressed ADAMTS9-AS1 and silencing TGFBR2 together could reverse this regulatory effect (p < 0.05), indicating increased tumor proliferation (Figure 4J). Altogether, ADAMTS9-AS1 restrained the proliferation and invasion of breast cancer via stimulating the level of TGFBR2 through miR-301b-3p.

ADAMTS9-AS1 Traps miR-301b-3p to Target TGFBR2, Thereby Regulating the JAK STAT Signaling Pathway

GSEA analysis result revealed that TGFBR2 was notably enriched in the JAK STAT signaling (Figure 5A), suggesting that ADAMTS9-AS1/miR-301b-3p/TGFBR2 signaling axis may regulate breast cancer cells through the JAK STAT signaling pathway. For further verification, ADAMTS9-AS1 was firstly overexpressed in cell line MDA-MB-231. Corresponding detection uncovered that the protein levels of p-STAT3, p-JAK2, p-STAT1 and p-JAK1 in the oe-ADAMTS9-AS1 group were notably repressed (p < 0.05) (Figure 5B), implying that ADAMTS9-AS1 constrained the JAK STAT signaling pathway. Further analyses depicted that the protein levels of p-STAT3, p-STAT1, p-JAK2, and p-JAK1 were notably boosted after silencing ADAMTS9-AS1, while the regulatory effects were reversed in si-ADAMTS9-AS1 + miR-301b-3p inhibitor (p < 0.05) (Figure 5C). Generally, silencing ADAMTS9-AS1 activated the JAK STAT signaling pathway, while silencing miR-301b-3p could reverse this impact.

FIGURE 5

Figure 5. ADAMTS9-AS1 regulates the JAK STAT signaling pathway by targeting TGFBR2 through miR-301b-3p. (A) GSEA enrichment pathway of TGFBR2 (JAK STAT signaling pathway). (B) The expression of JAK STAT signaling-relevant proteins. (C) The expression of JAK STAT signaling-relevant proteins in si-NC + inhibitor NC, si-ADAMTS9-AS1 + inhibitor NC and si-ADAMTS9-AS1 + miR-301b-3p inhibitor. (D) Levels of JAK STAT signaling-relevant proteins in oe-NC + si-NC, oe-ADAMTS9-AS1 + si-NC and oe-ADAMTS9-AS1 + si-TGFBR2.

MDA-MB-231 cells divided into the following transfecting groups: oe-NC + si-NC, oe-ADAMTS9-AS1 + si-NC and oe-ADAMTS9-AS1 + si-TGFBR2. Western blot examined that p-JAK2, p-STAT3, p-JAK1 and p-STAT1 protein levels were notably reduced in oe-ADAMTS9-AS1 + si-NC group, while co-transfection of oe-ADAMTS9-AS1 and si-TGFBR2 reversed the regulatory effect caused by overexpressed ADAMTS9-AS1 (p < 0.05) (Figure 5D), suggesting that ADAMTS9-AS1 inhibited the JAK STAT signaling pathway, while silencing TGFBR2 reversed that inhibitory effect. Overall, ADAMTS9-AS1 regulated the JAK STAT signaling pathway by targeting TGFBR2 through miR-301b-3p.

ADAMTS9-AS1 Constrains Breast Cancer Cell Proliferation and Invasion via Regulating JAK STAT Signaling Pathway

MCF-7 cells were simultaneously transfected with si-NC and si-ADAMTS9-AS1, and treated with JAK STAT signaling pathway inhibitor INCB018424 (HY-50856, MedchemExpress; Punwani et al., 2012). qRT-PCR result indicated that the expression of ADAMTS9-AS1 in si-ADAMTS9-AS1 + DMSO group was markedly lower than that in si-NC + DMSO group (p < 0.05). While there was no significant difference in ADAMTS9-AS1 expression between the si-ADAMTS9-AS1 + DMSO group and the si-ADAMTS9-AS1 + INCB018424 group (Figure 6A). The expression of JAK STAT signaling pathway-related proteins was detected by western blot. The result suggested that compared with si-NC + DMSO group, the protein expression levels of p-JAK1, p-JAK2, p-STAT1 and p-STAT3 in si-ADAMTS9-AS1 + DMSO group were remarkably increased. Compared with si-ADAMTS9-AS1 + DMSO group, the protein expression levels in si-ADAMTS9-AS1 + INCB018424 group were notably decreased (p < 0.05) (Figure 6B). In vitro colony formation (Figure 6C), EdU (Figure 6D), and Transwell (Figure 6E) assays indicated that compared with the si-NC + DMSO group, the proliferative ability and invasive ability of breast cancer cells in the si-ADAMTS9-AS1 + DMSO group were prominently increased. Compared with the si-ADAMTS9-AS1 + DMSO group, the proliferative and invasive abilities of breast cancer cells in si-ADAMTS9-AS1 + INCB018424 group were markedly decreased (p < 0.05). Western blot result exhibited that compared with si-NC + DMSO group, the expression levels of Ki67, PCNA, MMP-2 and MMP-9 in si-ADAMTS9-AS1 + DMSO group were markedly boosted, indicating increased tumor proliferative activity. Compared with si-ADAMTS9-AS1 + DMSO group, the expression levels of Ki67, PCNA, MMP-2 and MMP-9 in si-ADAMTS9-AS1 + INCB018424 group were prominently constrained, indicating decreased tumor proliferative activity (p < 0.05) (Figure 6F). In general, ADAMTS9-AS1 could suppress the proliferation and invasion of breast cancer cells by regulating the JAK STAT signaling pathway. Besides, JAK inhibitor did not impact TGFBR2 levels according to detection results (Supplementary Figure 3).

FIGURE 6

Figure 6. ADAMTS9-AS1 represses cancer invasion and proliferation by regulating JAK STAT signaling pathway. (A) ADAMTS9-AS1 expression. (B) Levels of JAK STAT signaling-relevant proteins. (C) The proliferation of breast cancer cells in vitro. (D) The proliferation of breast cancer cells (400×). (E) Cancer cell invasion (100×). (F) Levels of Ki67, PCNA, MMP-2 and MMP-9. *p < 0.05.

Discussion

Exposed to metastasis and recurrence, patients with breast cancer bear unsatisfied therapy efficacy and survival. Mechanism for invasion and proliferation is therefore urgently warranted for depiction. ADAMTS9-AS1 was noticeably inactivated in breast cancer cells here. ADAMTS9-AS1 is a dysregulated antisense lncRNA in many tumors and is involved in the pathogenesis of a variety of human malignancies (Wang et al., 2016; Xing et al., 2018). For example, in colorectal cancer, lncRNA ADAMTS9-AS1 constrains invasion and migration of colorectal cancer (Li et al., 2020). ADAMTS9-AS1 also represses tumor-relevant prostate cancer behaviors (Zhou et al., 2021). The present paper found that ADAMTS9-AS1 could repress cancer proliferation and invasion through in vitro functional experiments. Nevertheless, there are reports that are contrary to our research results. For instance, overexpression of ADAMTS9-AS1 aggravates hepatocellular carcinoma (Zhang et al., 2020) and colorectal cancer development (Chen et al., 2020). We speculated that the different roles of ADAMTS9-A1 may be attributed to different cancer types.

To date, lncRNA/miRNA axis has been a limelight in the field of cancer. LncRNA can trap miRNA, regulate mRNA, and modulate progression of various cancers (Ao et al., 2019; Tu et al., 2019). This paper verified miR-301b-3p/ADAMTS9-AS1/TGFBR2 interaction by bioinformatics method, RIP and dual luciferase assay. This paper discovered relevance of miR-301b-3p and poor prognosis of breast cancer patients, and the cancer-promoting effect of si-ADAMTS9-AS1 could be reversed when cells with miR-301b-3p inhibitor at the same time. TGFBR2 was reported as a functional mRNA. For instance, TGFBR2 down-regulation will promote the malignant progression of prostate cancer (Zhou H. et al., 2018). MiR-3139 targets TGFBR2 expression to promote cell migration and invasion in colorectal cancer (He et al., 2019). Here, inhibiting the expression of TGFBR2 could weaken the tumor suppressor effect of ADAMTS9-AS1, suggesting that TGFBR2 was an inhibitor in breast cancer, which is similar to the results of previous studies.

Subsequently, we conducted enrichment analysis on TGFBR2 and found that it was mainly enriched in the JAK STAT signaling. The STAT3 protein is widely present in a variety of tumor tissue, affecting many biological behaviors of tumor cells, and acting as an oncogene, while JAK is responsible for the activation of STAT (Bournazou and Bromberg, 2013; Panni et al., 2014; Wen et al., 2015). Dysregulation of JAK STAT signaling pertains to breast cancer metastasis along with high risk of recurrence (Khanna et al., 2018). In addition, activated STAT3 triggers characteristics of cancer, including angiogenesis (Chung et al., 2014; Banerjee and Resat, 2016). Herein, in combination with rescue experiment, we found that ADAMTS9-AS1 constrained the activation of the JAK STAT signaling through miR-301b-3p/TGFBR2 axis. Thus, cancer invasion and proliferation were repressed.

The paper confirmed that ADAMTS9-AS1 was inactivated in breast cancer cells. ADAMTS9-AS1 could sequester miR-301b-3p to mediate TGFBR2 protein, thereby affecting the JAK STAT signaling, and repressing tumor progression. ADAMTS9-AS1/miR-301b-3p/TGFBR2 regulatory network may provide a more effective clinical treatment strategy for breast cancer patients. ADAMTS9-AS1 might serve as a promising factor for breast cancer. However, this study also has some limitations. First, we lack tissue samples and in vivo experiments to fully confirm the accuracy of the results. Second, detailed investigations on the genes that make up the lncRNA ADAMTS9-AS1/miR-301b-3p/TGFBR2 axis can further understand the mechanism of lncRNA ADAMTS9-AS1 inhibiting breast cancer progression. In the future, we will add animal experiments, and also pay more attention to the relationship between ADAMTS9-AS1 and other potential target miRNAs. We hope to offer a more perfect basis for therapeutic strategy of breast cancer.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author Contributions

JC contributed to data analysis and manuscript writing. LC contributed to experiment concrete implementation. WZ acquired the data. RW looked up the literature. XW revised the manuscript. ZC designed the study and wrote the manuscript. All authors read and approved the final manuscript.

Funding

This study was supported by the funds from the Medical Science and Technology Program of Zhejiang Province (2013KYA026 and 2015DTA004). This study was supported by the funds from the National Nature Science Foundation of China (81672597).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.719993/full#supplementary-material

Supplementary Figure 1 | Association between miR-301b-3p and patient’s clinicalstage and TNM stage. (A) Box plot of miR-301b-3p expression at different clinical stages; (B) Box plot of miR-301b-3p expression at different T stages; (C) Box plot of miR-301b-3p expression at different N stages; (D) Box plot of miR-301b-3p expression at different M stages.

Supplementary Figure 2 | Cell morphology of MCF-7 upon overexpression or silence of ADAMTS9-AS1.

Supplementary Figure 3 | TGFBR2 protein and mRNA in MCF-10A and MCF-7 cells.

Supplementary Table 1 | Clinical information of patients.

Supplementary Table 2 | Detailed information of 1,017 DElncRNAs.

Supplementary Table 3 | Detailed information of 74 DEmiRNAs.

Supplementary Table 4 | Detailed information of 2,161 DEmRNAs.

Supplementary Table 5 | Binding between 60 lncRNA and 18 miRNAs.

Footnotes

1. ^ https://portal.gdc.cancer.gov/
2. ^ http://www.mircode.org/
3. ^ http://mirdb.org/
4. ^ http://www.targetscan.org/
5. ^ http://mirtarbase.mbc.nctu.EdU.tw/php/search.php

References

Ao, X., Jiang, M., Zhou, J., Liang, H., Xia, H., and Chen, G. (2019). lincRNA-p21 inhibits the progression of non-small cell lung cancer via targeting miR-17-5p. Oncol. Rep. 41, 789–800. doi: 10.3892/or.2018.6900

PubMed Abstract | CrossRef Full Text | Google Scholar

Banerjee, K., and Resat, H. (2016). Constitutive activation of STAT3 in breast cancer cells: a review. Int. J. Cancer 138, 2570–2578. doi: 10.1002/ijc.29923

PubMed Abstract | CrossRef Full Text | Google Scholar

Boo, L., Ho, W. Y., Ali, N. M., Yeap, S. K., Ky, H., Chan, K. G., et al. (2016). MiRNA transcriptome profiling of spheroid-enriched cells with cancer stem cell properties in human breast MCF-7 cell line. Int. J. Biol. Sci. 12, 427–445. doi: 10.7150/ijbs.12777

PubMed Abstract | CrossRef Full Text | Google Scholar

Bournazou, E., and Bromberg, J. (2013). Targeting the tumor microenvironment: JAK-STAT3 signaling. JAKSTAT 2:e23828. doi: 10.4161/jkst.23828

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W., Tu, Q., Yu, L., Xu, Y., Yu, G., Jia, B., et al. (2020). LncRNA ADAMTS9-AS1, as prognostic marker, promotes cell proliferation and EMT in colorectal cancer. Hum. Cell 33, 1133–1141. doi: 10.1007/s13577-020-00388-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Chung, S. S., Giehl, N., Wu, Y., and Vadgama, J. V. (2014). STAT3 activation in HER2-overexpressing breast cancer promotes epithelial-mesenchymal transition and cancer stem cell traits. Int. J. Oncol. 44, 403–411. doi: 10.3892/ijo.2013.2195

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, W., Wang, S., Zhou, Q., Li, X., Chu, J., Chang, Z., et al. (2013). ADAMTS9 is a functional tumor suppressor through inhibiting AKT/mTOR pathway and associated with poor survival in gastric cancer. Oncogene 32, 3319–3328. doi: 10.1038/onc.2012.359

PubMed Abstract | CrossRef Full Text | Google Scholar

Fatica, A., and Bozzoni, I. (2014). Long non-coding RNAs: new players in cell differentiation and development. Nat. Rev. Genet. 15, 7–21. doi: 10.1038/nrg3606

PubMed Abstract | CrossRef Full Text | Google Scholar

Gilam, A., Conde, J., Weissglas-Volkov, D., Oliva, N., Friedman, E., Artzi, N., et al. (2016). Local microRNA delivery targets Palladin and prevents metastatic breast cancer. Nat. Commun. 7:12868. doi: 10.1038/ncomms12868

PubMed Abstract | CrossRef Full Text | Google Scholar

He, H., Zhao, X., Zhu, Z., Du, L., Chen, E., Liu, S., et al. (2019). MicroRNA-3191 promotes migration and invasion by downregulating TGFBR2 in colorectal cancer. J. Biochem. Mol. Toxicol. 33:e22308. doi: 10.1002/jbt.22308

PubMed Abstract | CrossRef Full Text | Google Scholar

Khanna, P., Lee, J. S., Sereemaspun, A., Lee, H., and Baeg, G. H. (2018). GRAMD1B regulates cell migration in breast cancer cells through JAK/STAT and Akt signaling. Sci. Rep. 8:9511. doi: 10.1038/s41598-018-27864-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Laulin, J. P., and Brudieux, R. (1990). Atrial natriuretic factor is unlikely to be involved in the reduced aldosterone production in the Brattleboro rat. Regul. Pept. 31, 157–165. doi: 10.1016/0167-0115(90)90002-e

CrossRef Full Text | Google Scholar

Li, N., Li, J., Mi, Q., Xie, Y., Li, P., Wang, L., et al. (2020). Long non-coding RNA ADAMTS9-AS1 suppresses colorectal cancer by inhibiting the Wnt/β-catenin signalling pathway and is a potential diagnostic biomarker. J. Cell. Mol. Med. 24, 11318–11329. doi: 10.1111/jcmm.15713

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Yao, Q., Zhao, S., Wang, Y., Li, Y., and Wang, Z. (2017). Comprehensive analysis of differential co-expression patterns reveal transcriptional dysregulation mechanism and identify novel prognostic lncRNAs in esophageal squamous cell carcinoma. Oncotargets Ther. 10, 3095–3105. doi: 10.2147/OTT.S135312

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J. J., Zhang, X., and Wu, X. H. (2018). miR-93 promotes the growth and invasion of prostate cancer by upregulating its target genes TGFBR2, ITGB8, and LATS2. Mol. Ther. Oncol. 11, 14–19. doi: 10.1016/j.omto.2018.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Lo, P. H., Leung, A. C., Kwok, C. Y., Cheung, W. S., Ko, J. M., Yang, L. C., et al. (2007). Identification of a tumor suppressive critical region mapping to 3p14.2 in esophageal squamous cell carcinoma and studies of a candidate tumor suppressor gene, ADAMTS9. Oncogene 26, 148–157. doi: 10.1038/sj.onc.1209767

PubMed Abstract | CrossRef Full Text | Google Scholar

Lung, H. L., Lo, P. H. Y., Xie, D., Apte, S. S., Cheung, A. K. L., Cheng, Y., et al. (2008). Characterization of a novel epigenetically-silenced, growth-suppressive gene, ADAMTS9, and its association with lymph node metastases in nasopharyngeal carcinoma. Int. J. Cancer 123, 401–408. doi: 10.1002/ijc.23528

PubMed Abstract | CrossRef Full Text | Google Scholar

Mercer, T. R., Dinger, M. E., and Mattick, J. S. (2009). Long non-coding RNAs: insights into functions. Nat. rev. Genet. 10, 155–159. doi: 10.1038/nrg2521

PubMed Abstract | CrossRef Full Text | Google Scholar

Nassar, F. J., Nasr, R., and Talhouk, R. (2017). MicroRNAs as biomarkers for early breast cancer diagnosis, prognosis and therapy prediction. Pharmacol. Ther. 172, 34–49. doi: 10.1016/j.pharmthera.2016.11.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Ocak, Z., Acar, M., Gunduz, E., Gunduz, M., Demircan, K., Uyeturk, U., et al. (2013). Effect of hypericin on the ADAMTS-9 and ADAMTS-8 gene expression in MCF7 breast cancer cells. Eur. Rev. Med. Pharmacol. Sci. 17, 1185–1190.

Google Scholar

Panni, R. Z., Sanford, D. E., Belt, B. A., Mitchem, J. B., Worley, L. A., Goetz, B. D., et al. (2014). Tumor-induced STAT3 activation in monocytic myeloid-derived suppressor cells enhances stemness and mesenchymal properties in human pancreatic cancer. Cancer Immunol. Immunother. 63, 513–528. doi: 10.1007/s00262-014-1527-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, W., He, D., Shan, B., Wang, J., Shi, W., Zhao, W., et al. (2019). LINC81507 act as a competing endogenous RNA of miR-199b-5p to facilitate NSCLC proliferation and metastasis via regulating the CAV1/STAT3 pathway. Cell Death Dis. 10:533. doi: 10.1038/s41419-019-1740-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Punwani, N., Scherle, P., Flores, R., Shi, J., Liang, J., Yeleswaram, S., et al. (2012). Preliminary clinical activity of a topical JAK1/2 inhibitor in the treatment of psoriasis. J. Am. Acad. Dermatol. 67, 658–664. doi: 10.1016/j.jaad.2011.12.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, S. N., Li, K., Liu, Y., Yang, C. L., He, C. Y., and Wang, H. R. (2019). Down-regulation of long noncoding RNA PVT1 inhibits esophageal carcinoma cell migration and invasion and promotes cell apoptosis via microRNA-145-mediated inhibition of FSCN1. Mol. Oncol. 13, 2554–2573. doi: 10.1002/1878-0261.12555

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheu, J. J., Lee, C. H., Ko, J. Y., Tsao, G. S., Wu, C. C., Fang, C. Y., et al. (2009). Chromosome 3p12.3-p14.2 and 3q26.2-q26.32 are genomic markers for prognosis of advanced nasopharyngeal carcinoma. Cancer Epidemiol. Biomarkers Prev. 18, 2709–2716. doi: 10.1158/1055-9965.EPI-09-0349

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Y. S., Zhao, Z., Yang, Z. N., Xu, F., Lu, H. J., Zhu, Z. Y., et al. (2017). Risk factors and preventions of breast cancer. Int. J. Biol. Sci. 13, 1387–1397. doi: 10.7150/ijbs.21635

PubMed Abstract | CrossRef Full Text | Google Scholar

Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 71, 209–249. doi: 10.3322/caac.21660

PubMed Abstract | CrossRef Full Text | Google Scholar

Tu, J., Zhao, Z., Xu, M., Chen, M., Weng, Q., Wang, J., et al. (2019). LINC00707 contributes to hepatocellular carcinoma progression via sponging miR-206 to increase CDK14. J. Cell. Physiol. 234, 10615–10624.

Google Scholar

Wang, H., Fu, Z., Dai, C., Cao, J., Liu, X., Xu, J., et al. (2016). LncRNAs expression profiling in normal ovary, benign ovarian cyst and malignant epithelial ovarian cancer. Sci. Rep. 6:38983. doi: 10.1038/srep38983

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Sun, L., Wang, L., Liu, Z., Li, Q., Yao, B., et al. (2018). Long non-coding RNA DSCR8 acts as a molecular sponge for miR-485-5p to activate Wnt/beta-catenin signal pathway in hepatocellular carcinoma. Cell Death Dis. 9:851. doi: 10.1038/s41419-018-0937-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen, W., Wu, J., Liu, L., Tian, Y., Buettner, R., Hsieh, M. Y., et al. (2015). Synergistic anti-tumor effect of combined inhibition of EGFR and JAK/STAT3 pathways in human ovarian cancer. Mol. Cancer 14:100. doi: 10.1186/s12943-015-0366-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Xing, Y., Zhao, Z., Zhu, Y., Zhao, L., Zhu, A., Piao, D., et al. (2018). Comprehensive analysis of differential expression profiles of mRNAs and lncRNAs and identification of a 14-lncRNA prognostic signature for patients with colon adenocarcinoma. Oncol. Rep. 39, 2365–2375. doi: 10.3892/or.2018.6324

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Li, H., Hu, Y., and Wang, F. (2020). Long non-coding RNA ADAMTS9-AS1 exacerbates cell proliferation, migration, and invasion via triggering of the PI3K/AKT/mTOR pathway in hepatocellular carcinoma cells. Am. J. Transl. Res. 12, 5696–5707.

Google Scholar

Zhou, B., Guo, W., Sun, C., Zhang, B., and Zheng, F. (2018). Linc00462 promotes pancreatic cancer invasiveness through the miR-665/TGFBR1-TGFBR2/SMAD2/3 pathway. Cell Death Dis. 9:706. doi: 10.1038/s41419-018-0724-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, H., Wu, G., Ma, X., Xiao, J., Yu, G., Yang, C., et al. (2018). Attenuation of TGFBR2 expression and tumour progression in prostate cancer involve diverse hypoxia-regulated pathways. J. Exp. Clin. Cancer Res. 37:89. doi: 10.1186/s13046-018-0764-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Z., Wu, X., Zhou, Y., and Yan, W. (2021). LncRNA ADAMTS9-AS1 inhibits the progression of prostate cancer by modulating miR-142-5p/CCND1 axis. J. Gene Med. 23:e3331. doi: 10.1002/jgm.3331

PubMed Abstract | CrossRef Full Text | Google Scholar