Research Paper Volume 12, Issue 20 pp 20366—20379

miR-140-3p inhibits bladder cancer cell proliferation and invasion by targeting FOXQ1

Yuan Wang1, , Junwen Chen2, , Xia Wang2, , Kefeng Wang2, ,

  • 1 Department of General Surgery, Shengjing Hospital of China Medical University, Shenyang 110004, China
  • 2 Department of Urology, Shengjing Hospital of China Medical University, Shenyang 110004, China

Received: March 12, 2020       Accepted: July 7, 2020       Published: October 24, 2020      

https://doi.org/10.18632/aging.103828
How to Cite

Copyright: © 2020 Wang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Upregulation of the forkhead box protein Q1 (FOXQ1) promotes bladder cancer (BCa) cell growth and metastasis. Factors affecting FOXQ1 expression at the post-transcriptional level have not yet been identified. We performed cell proliferation, cell invasion, and tumorigenesis experiments to characterize the relationship between FOXQ1 and miR-140-3p. We found that FOXQ1 was significantly upregulated and miR-140-3p was significantly downregulated in BCa tissues. We also identified an inverse correlation between miR-140-3p and FOXQ1 expression in BCa tissues. Overexpression of miR-140-3p reduced FOXQ1 expression, suppressing BCa cell proliferation and invasion. A luciferase assay confirmed that miR-140-3p bound to the 3’-UTR of FOXQ1 mRNA and decreased its expression. In addition, we used a mouse xenograft model to demonstrate that miR-140-3p suppressed tumor cell growth in vivo. Our findings suggest that miR-140-3p suppresses BCa cell proliferation and invasion by directly decreasing FOXQ1 expression.

Introduction

Bladder cancer (BCa) is one of the most prevalent malignant urological cancers worldwide. In 2019, approximately 80,470 patients were diagnosed with BCa in the United States, resulting in approximately 17,670 deaths [1]. BCa can be categorized into two groups, non-muscle-invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC). These cancers have 5-year survival rates of 60% and 90%, respectively [2]. The recurrence rate of BCa is high [3]. Approximately 80% of patients will die within 5 years after distant metastases are discovered [4]. Transurethral resection of bladder tumor is considered to be an effective treatment for NMIBC [5], while radical cystectomy and/or systemic chemotherapy provides limited benefits to late stage MIBC [6, 7].

The majority of transcripts in the human genome are derived from non-coding RNA (ncRNA) [8]. MicroRNA (miRNA) is a class of ncRNA that is 19-22 nucleotides long. These transcripts act as endogenous, post-transcriptional regulators of gene expression by binding to the 3’-untranslated region (3’-UTR), 5’-UTR or coding sequence region of the target mRNA. miRNA performs its functions through either inhibition of protein translation or through mRNA degradation. Microarray expression profiles and miRNA bioinformatics reveal aberrant expression in BCa tissues [9]. miRNA can act as either a tumor suppressor or an oncogene in different types of cancer. It has been reported that miR-140-3p functions as a tumor suppressor in breast cancer [10], lung cancer [11] and lymphoma [12].

The forkhead box protein Q1 (FOXQ1) gene is a member of the FOX gene family, which encode proteins characterized by a conserved 110-amino acid DNA-binding motif called the forkhead or winged helix domain. Recent studies revealed that upregulation of FOXQ1 was associated with progression of many human tumors, including those found in hepatocellular carcinoma [13], gastric cancer [14], colorectal carcinoma [15], prostate cancer [16], lung cancer [17], laryngeal carcinoma [18], esophageal cancer [19], pancreatic cancer [20], breast cancer [21], and BCa [22].

In our study, we detected the expression of miR-140-3p and FOXQ1 in BCa tissues. We investigated the regulation of FOXQ1 by miR-140-3p, as well as the biological functions of FOXQ1 targeted by miR-140-3p in BCa cells.

Results

FOXQ1 promotes BCa cell proliferation and invasion

Knockdown of FOXQ1 inhibits the migration and invasion of BCa cells [22]. We conducted a western blot assay and confirmed that FOXQ1 was effectively upregulated or downregulated after transfection of functional FOXQ1-cDNA or FOXQ1-shRNA in T24 and UMUC3 cells, respectively (Figure 1A). Furthermore, a cell proliferation assay demonstrated that overexpression of FOXQ1 increased cell proliferation, while knockdown of FOXQ1 inhibited cell proliferation in T24 and UMUC3 cells (Figure 1B). A transwell invasion assay in T24 and UMUC3 cells revealed similar results (Figure 1C).

FOXQ1 promotes BCa cell proliferation and invasion. (A) Verification of FOXQ1 overexpression and knockdown in T24 and UMUC3 cells by western blot assay. Approximately 50 ug of protein was loaded into each lane 2 to 4 days after transfection. (B) T24 and UMUC3 cells were transfected with pWPI-FOXQ1-cDNA, pWPI-vector, pLKO-FOXQ1-shRNA or pLKO-vector. Cell growth was measured by MTT assay. (C) Transwell invasion assays were performed using UMUC3 cells (pWPI and oeFOXQ1) and T24 cells (pLKO and shFOXQ1). The invaded cells were counted in 10 randomly chosen microscopic fields (100X) of each experiment and pooled. (D) Analysis of microarray sequencing from NCBI GEO Datasets (GSE40355, GPL13497) showed the FOXQ1 mRNA level in 16 BCa and 8 nonmalignant bladder tissue samples. (E) Curves of BCa patient OS were analyzed according to FOXQ1 expression (data were download from TCGA). (F) FOXQ1 expression in 30 paired human primary BCa tissues and adjacent normal bladder tissues. (B–D) Each sample was run in triplicate and used in multiple experiments to determine the mean ± SD. *P P

Figure 1. FOXQ1 promotes BCa cell proliferation and invasion. (A) Verification of FOXQ1 overexpression and knockdown in T24 and UMUC3 cells by western blot assay. Approximately 50 ug of protein was loaded into each lane 2 to 4 days after transfection. (B) T24 and UMUC3 cells were transfected with pWPI-FOXQ1-cDNA, pWPI-vector, pLKO-FOXQ1-shRNA or pLKO-vector. Cell growth was measured by MTT assay. (C) Transwell invasion assays were performed using UMUC3 cells (pWPI and oeFOXQ1) and T24 cells (pLKO and shFOXQ1). The invaded cells were counted in 10 randomly chosen microscopic fields (100X) of each experiment and pooled. (D) Analysis of microarray sequencing from NCBI GEO Datasets (GSE40355, GPL13497) showed the FOXQ1 mRNA level in 16 BCa and 8 nonmalignant bladder tissue samples. (E) Curves of BCa patient OS were analyzed according to FOXQ1 expression (data were download from TCGA). (F) FOXQ1 expression in 30 paired human primary BCa tissues and adjacent normal bladder tissues. (BD) Each sample was run in triplicate and used in multiple experiments to determine the mean ± SD. *P < 0.05; **P < 0.01 compared to controls.

Analysis of GEO DataSets (GSE40355, GPL13497) showed that FOXQ1 expression was significantly higher in 16 BCa tissues than in 8 nonmalignant bladder tissues (Figure 1D). Data from The Cancer Genome Atlas (TCGA) showed that BCa patients with higher FOXQ1 expression had a shorter overall survival (OS) than patients with lower FOXQ1 expression (Figure 1E). We then investigated FOXQ1 expression in BCa samples and adjacent normal bladder samples from 30 patients from our department. The results revealed that the expression of FOXQ1 was significantly higher in BCa samples compared with adjacent normal bladder samples (Figure 1F). Together, results from Figure 1A1F indicated that FOXQ1 expression was higher in BCa samples than in adjacent normal bladder samples.

miR-140-3p downregulates FOXQ1 expression and suppresses BCa cell proliferation and invasion

miRNA regulates the targeted mRNA at the post-transcriptional level to depress the mRNA degradation or protein translation. Therefore, miRNA might play a vital role in regulating FOXQ1 expression in BCa. We focused our attention on miRNAs that could decrease FOXQ1 expression and inhibit BCa cell proliferation and invasion.

Using miRNA target-prediction databases (miRDB, Targetscan and MicroCosm), we determined that miR-140-3p might target FOXQ1 mRNA. We explored the effects of miR-140-3p on FOXQ1 expression and found that miR-140-3p decreased FOXQ1 expression in T24 and UMUC3 cells (Figure 2A). We then quantified miR-140-3p expression in 30 pairs of BCa samples and adjacent normal bladder samples. The expression of miR-140-3p was significantly lower in BCa samples compared with adjacent normal bladder samples (Figure 2B). An inverse correlation was found between the expression of miR-140-3p and FOXQ1 in the 30 paired samples above (Figure 2C). We also observed the expression of FOXQ1 in patients tissue samples using immunohistochemical staining. The results showed that cells from normal bladder samples expressed lower level of FOXQ1 than cells from BCa samples. What’s more, we found more positive FOXQ1 expression in patients with lower miR-140-3p expression than in patients with higher miR-140-3p expression (Figure 2D). Analysis of GEO DataSets (GSE40355, GPL8227) also verified that the expression of miR-140-3p was significantly lower in 16 BCa tissues than in 8 nonmalignant bladder tissues (Figure 2E). We found from TCGA data analysis that patients with lower miR-140-3p expression had a shorter OS (HR=1.521, P=0.0136) than patients with higher miR-140-3p expression (Figure 2F).

miR-140-3p downregulates FOXQ1 expression and suppresses BCa cell proliferation and invasion. (A) Overexpression of miR-140-3p in T24 and UMUC3 cells to determine its effect on FOXQ1 expression by western blot assay. Approximately 50 ug of protein was loaded into each lane 2 to 4 days after transfection. (B) miR-140-3p expression in 30 paired human primary BCa and adjacent normal bladder tissues. (C) Correlation analysis of miR-140-3p and FOXQ1 mRNA level was performed using the Pearson correlation coefficient. (D) Representative immunohistochemical staining of FOXQ1 in BCa and nonmalignant bladder tissue samples. (E) Analysis of microarray sequencing from NCBI GEO Datasets (GSE40355, GPL8227) showed miR-140-3p expression in 16 BCa and 8 nonmalignant bladder tissue samples. (F) OS curves of BCa patients were analyzed according to miR-140-3p expression. Data were download from TCGA. (G) T24 cells were transfected with NC and miR-140-3p mimics. (H) UMUC3 cells were transfected with miR-140-3p inhibitor. Cell growth was measured by MTT assay. (I, J) Transwell invasion assays were performed by transfecting T24 cells with a negative control (NC) or miR-140-3p mimics (I) and UMUC3 cells with NC or miR-140-3p inhibitor (J). The invaded cells were counted in 10 randomly chosen microscopic fields (100X) of each experiment and pooled. (G–J) Each sample was run in triplicate and used in multiple experiments to determine the mean ± SD. *P P

Figure 2. miR-140-3p downregulates FOXQ1 expression and suppresses BCa cell proliferation and invasion. (A) Overexpression of miR-140-3p in T24 and UMUC3 cells to determine its effect on FOXQ1 expression by western blot assay. Approximately 50 ug of protein was loaded into each lane 2 to 4 days after transfection. (B) miR-140-3p expression in 30 paired human primary BCa and adjacent normal bladder tissues. (C) Correlation analysis of miR-140-3p and FOXQ1 mRNA level was performed using the Pearson correlation coefficient. (D) Representative immunohistochemical staining of FOXQ1 in BCa and nonmalignant bladder tissue samples. (E) Analysis of microarray sequencing from NCBI GEO Datasets (GSE40355, GPL8227) showed miR-140-3p expression in 16 BCa and 8 nonmalignant bladder tissue samples. (F) OS curves of BCa patients were analyzed according to miR-140-3p expression. Data were download from TCGA. (G) T24 cells were transfected with NC and miR-140-3p mimics. (H) UMUC3 cells were transfected with miR-140-3p inhibitor. Cell growth was measured by MTT assay. (I, J) Transwell invasion assays were performed by transfecting T24 cells with a negative control (NC) or miR-140-3p mimics (I) and UMUC3 cells with NC or miR-140-3p inhibitor (J). The invaded cells were counted in 10 randomly chosen microscopic fields (100X) of each experiment and pooled. (GJ) Each sample was run in triplicate and used in multiple experiments to determine the mean ± SD. *P < 0.05; **P < 0.01 compared to controls.

We quantified the expression of miR-140-3p in T24 and UMUC3 cells. We found that miR-140-3p expression was higher in UMUC3 cells compared with T24 cells (Supplementary Figure 1A). miR-140-3p mimics was transfected into T24 cells, while an miR-140-3p inhibitor was transfected into UMUC3 cells. We explored whether miR-140-3p regulated FOXQ1 expression in BCa cells. The results of western blot analyses verified that miR-140-3p mimics reduced FOXQ1 expression in T24 cells and the miR-140-3p inhibitor increased FOXQ1 expression in UMUC3 cells (Figure 4A, 4B). qRT-PCR was conducted to verify that miR-140-3p was effectively overexpressed in T24 cells and knocked down in UMUC3 cells after transfection (Supplementary Figure 1B). A cell proliferation assay revealed that miR-140-3p overexpression significantly decreased the growth rate of T24 cells (Figure 2G). In contrast, miR-140-3p knockdown had the opposite effect on the growth rate in UMUC3 cells (Figure 2H). Moreover, we found that the cell invasion activity was significantly inhibited by miR-140-3p overexpression in T24 cells (Figure 2I). In contrast, miR-140-3p knockdown had the opposite effect on the cell invasion activity in UMUC3 cells (Figure 2J). These data suggested that miR-140-3p functioned as a tumor suppressor to inhibit BCa cell proliferation and invasion.

miR-140-3p directly regulates FOXQ1 expression by targeting the 3'-UTR

To verify that FOXQ1 is a potential downstream target of miR-140-3p, we used online bioinformatics databases (miRDB, Targetscan and MicroCosm) for analysis. The predicted interactions between miR-140-3p and the binding sites on the 3’-UTR of FOXQ1 mRNA are illustrated in Figure 3A. Results from luciferase reporter assays showed that in T24 and UMUC3 cells, transfection of miR-140-3p and wild-type sequences of the FOXQ1 3’-UTR decreased luminescence intensity. Altering miR-140-3p expression did not affect the luminescence intensity of the mutant FOXQ1 3’-UTR in either cell type (Figure 3B, 3C). Collectively, these data supported the hypothesis that miR-140-3p inhibited FOXQ1 expression by directly interacting with the 3'-UTR.

miR-140-3p directly regulates FOXQ1 expression by targeting the 3'-UTR. (A) Sequence alignment of FOXQ1 3’-UTR with wild-type (WT) versus mutant potential miR-140-3p targeting sites using bioinformatics online databases (miRDB, Targetscan, and MicroCosm). (B, C) Co-transfection of wild-type or mutant seed regions of FOXQ1 3’-UTR constructs with miR-140-3p in T24 (B) and UMUC3 cells (C). The luciferase assay was applied to detect the luciferase activity. *P

Figure 3. miR-140-3p directly regulates FOXQ1 expression by targeting the 3'-UTR. (A) Sequence alignment of FOXQ1 3’-UTR with wild-type (WT) versus mutant potential miR-140-3p targeting sites using bioinformatics online databases (miRDB, Targetscan, and MicroCosm). (B, C) Co-transfection of wild-type or mutant seed regions of FOXQ1 3’-UTR constructs with miR-140-3p in T24 (B) and UMUC3 cells (C). The luciferase assay was applied to detect the luciferase activity. *P < 0.05 compared to controls.

miR-140-3p suppresses proliferation and invasion of BCa cells by reducing FOXQ1 expression

We explored whether miR-140-3p regulated FOXQ1 expression in BCa cells. The results of western blot analyses verified that miR-140-3p mimics reduced FOXQ1 expression in T24 cells and miR-140-3p inhibitor increased FOXQ1 expression in UMUC3 cells (Figure 4A, 4B). To determine whether miR-140-3p is involved in the FOXQ1-induced promotion of BCa cell proliferation and invasion, we conducted rescue experiments in T24 and UMUC3 cells. The results showed that miR-140-3p mimics partially reversed the increase mediated by overexpression of FOXQ1 on T24 cell proliferation (Figure 4C). In UMUC3 cells, a miR-140-3p inhibitor partially reversed the proliferation suppression mediated by knockdown of FOXQ1 (Figure 4D). We obtained similar results in both T24 and UMUC3 cells using an invasion assay in place of the cell proliferation assay (Figure 4E, 4F). These data verified that miR-140-3p suppressed FOXQ1-induced proliferation and invasion of BCa cells.

miR-140-3p suppresses proliferation and invasion of BCa cells by reducing FOXQ1 expression. (A, B), Western blot assay was performed to detect FOXQ1 expression. Approximately 50 ug of protein was loaded into each lane 2 to 4 days after transfection. (A) T24 cells were transfected with pWPI+NC, oeFOXQ1+NC, oeFOXQ1+miR-140-3p mimics or pWPI+miR-140-3p mimics. (B) UMUC3 cells were transfected with pLKO+NC, shFOXQ1+NC, shFOXQ1+inh-miR-140-3p or pLKO+ inh-miR-140-3p (C) An MTT rescue assay revealed that FOXQ1-increased cell proliferation could be reversed after adding miR-140-3p mimics to T24 cells. (D) An MTT rescue assay revealed that shFOXQ1-decreased cell proliferation could be reversed after adding miR-140-3p inhibitor to UMUC3 cells. (E) A transwell invasion assay revealed that FOXQ1-increased cell invasion could be reversed after adding miR-140-3p mimics to T24 cells. (F) A transwell invasion assay revealed that shFOXQ1-decreased cell invasion could be reversed after adding miR-140-3p inhibitor to UMUC3 cells. (C–F) Each sample was run in triplicate and used in multiple experiments to determine the mean ± SD. *P P

Figure 4. miR-140-3p suppresses proliferation and invasion of BCa cells by reducing FOXQ1 expression. (A, B), Western blot assay was performed to detect FOXQ1 expression. Approximately 50 ug of protein was loaded into each lane 2 to 4 days after transfection. (A) T24 cells were transfected with pWPI+NC, oeFOXQ1+NC, oeFOXQ1+miR-140-3p mimics or pWPI+miR-140-3p mimics. (B) UMUC3 cells were transfected with pLKO+NC, shFOXQ1+NC, shFOXQ1+inh-miR-140-3p or pLKO+ inh-miR-140-3p (C) An MTT rescue assay revealed that FOXQ1-increased cell proliferation could be reversed after adding miR-140-3p mimics to T24 cells. (D) An MTT rescue assay revealed that shFOXQ1-decreased cell proliferation could be reversed after adding miR-140-3p inhibitor to UMUC3 cells. (E) A transwell invasion assay revealed that FOXQ1-increased cell invasion could be reversed after adding miR-140-3p mimics to T24 cells. (F) A transwell invasion assay revealed that shFOXQ1-decreased cell invasion could be reversed after adding miR-140-3p inhibitor to UMUC3 cells. (CF) Each sample was run in triplicate and used in multiple experiments to determine the mean ± SD. *P < 0.05; **P < 0.01 compared to controls.

miR-140-3p suppresses the growth of BCa cells in vivo

To identify the effects of miR-140-3p on BCa cell growth, we transfected T24 cells with either miR-140-3p mimics or a negative control. The cells were injected into nude mice subcutaneously. Tumor volume was monitored weekly. The mice were sacrificed 4 weeks after injection. The weights and volumes of the xenografted tumors were measured. As expected, miR-140-3p mimics dramatically decreased the tumor weights and tumor volumes compared with those of the negative control group (Figure 5A5C). Our immunohistochemical staining demonstrated that both Ki-67 and FOXQ1 expression was significantly lower in the miR-140-3p mimics group than in negative control group (Figure 5D). We performed qRT-PCR to confirm that the expression of miR-140-3p increased in the miR-140-3p mimics group more than in negative control group (Figure 5E). The result of qRT-PCR also suggested that FOXQ1 expression was lower in the miR-140-3p overexpression group (Figure 5F). Together these data demonstrated that miR-140-3p suppressed the growth of BCa cells in vivo and might serve as a therapeutic marker in bladder cancer.

miR-140-3p suppresses the growth of BCa cells in vivo. (A) Macroscopic appearance of murine tumor xenografts. (B) Weights of tumors in 2 groups were measured using electronic scales. (C) Summary of tumor volume, which were measured weekly. (D) Representative immunohistochemical staining of Ki-67 and FOXQ1 in murine BCa cell xenografts. (E) The expression of miR-140-3p in xenografts was detected using qRT-PCR. (F) The expression of FOXQ1 in xenografts was detected using qRT-PCR. (B–F) Each sample was run in triplicate and used in multiple experiments to determine the mean ± SD. *P P

Figure 5. miR-140-3p suppresses the growth of BCa cells in vivo. (A) Macroscopic appearance of murine tumor xenografts. (B) Weights of tumors in 2 groups were measured using electronic scales. (C) Summary of tumor volume, which were measured weekly. (D) Representative immunohistochemical staining of Ki-67 and FOXQ1 in murine BCa cell xenografts. (E) The expression of miR-140-3p in xenografts was detected using qRT-PCR. (F) The expression of FOXQ1 in xenografts was detected using qRT-PCR. (BF) Each sample was run in triplicate and used in multiple experiments to determine the mean ± SD. *P < 0.05; **P < 0.01 compared to controls.

Discussion

The forkhead gene in Drosophila and the hepatocyte nuclear factor 3 alpha gene in rats were the first FOX family genes to be identified [23]. Members of the FOX gene family exist in a large range of organisms, functioning in a variety of tissues and playing vital roles in various biological processes. The FOX genes are involved in embryonic development [24], cell-cycle regulation [25], tissue-specific gene expression [26], cell signaling [27] and tumorigenesis [28].

The human FOXQ1 gene was first isolated and characterized in 2001 [29]. The FOXQ1 gene is expressed strongly in the trachea and stomach, and to a lesser degree in the salivary gland and bladder. In 2010, researchers found that in colorectal cancer, FOXQ1 promoted tumorigenicity and tumor growth by upregulating p21 expression, which in turn enhanced angiogenesis and anti-apoptosis. It has since been revealed that FOXQ1 functions as an oncogene in many human cancer types, including BCa [22, 30].

In our study, we found that the expression of FOXQ1 was higher in BCa tissues than adjacent normal bladder tissues. FOXQ1 promoted BCa cell proliferation and invasion. We also found that miRNA regulates FOXQ1 mRNA at a post-transcriptional level to repress mRNA degradation or protein translation. Investigators have discovered that miR-124, miR-506 and miR-342-3p suppressed nasopharyngeal carcinoma cell proliferation and metastasis by targeting FOXQ1 [3133]. Other researchers have revealed that miR-1271, miR-345 and miR-519 inhibited FOXQ1 in gastric cancer [14, 34, 35]. In colorectal cancer, miR-320, miR-342 and miR-106a restrained tumor cell growth by targeting FOXQ1 [3638]. Finally, lncRNA MALAT1 promoted BCa proliferation and metastasis by targeting miR-124/FOXQ1 expression [30].

A single gene can be targeted by multiple miRNAs to initiate various functions. We explored whether other miRNAs regulated FOXQ1 expression to influence BCa proliferation and invasion. We used in silico analysis to identify miR-140-3p as a candidate that regulated FOXQ1 expression in BCa cells. By literature review, we found that miR-140 had been shown to perform anti-cancer functions in many types of human cancer. Xie et al. reported that miR-140 suppressed PD-L1 and cyclin E expression to inhibit cell proliferation in non-small cell lung cancer (NSCLC) [39]. Fang et al. found that miR-140-5p downregulated YES1 expression to inhibit proliferation, invasion and migration of gastric cancer [40]. Lv et al. revealed that lncRNA-Unigene56159 silenced miR-140-5p to de-repress the expression of Slug, promoting cell invasion and EMT in hepatocellular carcinoma cells [41]. Wang et al. reported that miR-140-5p inhibited the proliferation, migration and invasion of BCa cells [42]. Here we found that miR-140-3p directly inhibited FOXQ1 expression in BCa cells. In addition, we determined that miR-140-3p inhibited FOXQ1 expression by targeting the 3'-UTR. Gain- and loss-of-function analysis revealed that miR-140-3p suppressed proliferation and invasion of BCa cells by reducing FOXQ1 expression. Furthermore, miR-140-3p inhibited the growth of BCa cells in a xenograft mouse model. These results revealed a network involving miR-140-3p and FOXQ1 that fine tunes the invasion and proliferation of BCa cells.

We found that the expression level of miR-140-3p was significantly lower in BCa tissues than in adjacent normal bladder tissues. We found similar results in GEO DataSets (GSE40355, GPL8227). Using TCGA database analysis, we found that patients with higher miR-140-3p expression had longer disease-free survival periods than patients with lower miR-140-3p expression. miR-140 was also downregulated in types of cancer [4345]. The above results suggested that miR-140 might work as a tumor suppressor miRNA in cancer.

It is well known that one miRNA can target multiple genes. For example, miR-140 suppresses tumor proliferation and metastasis by targeting insulin-like growth factor 1 receptor in NSCLC [46]. miR-140 also suppresses cell proliferation and invasion by targeting ATP8A1 in NSCLC [11]. To explore the role of this newly identified pathway in bladder tumorigenesis signaling, we upregulated miR-140-3p expression and found that the growth and invasion of BCa cells were inhibited, which imitated the function of FOXQ1 reduction by targeted siRNA. Furthermore, the overexpression of FOXQ1 partially attenuated the anti-proliferative and anti-invasive effects of miR-140-3p on BCa cells.

miRNA plays important roles in tumor onset and progression. Recent studies revealed that replacement of tumor suppressive miRNA or inhibition of oncogenic miRNA were possible strategies in cancer therapy [47]. In our study, we found that miR-140-3p downregulated FOXQ1 expression both in vitro and in vivo. Flamini et al. found that miR-140 replacement treatment combined with other drugs enhanced drug efficacy by reducing the invasion and migration ability of NSCLC [48]. Hence, miR-140-3p mimic replacement therapy is a candidate for BCa treatment.

Materials and Methods

Human tissue specimens

Human BCa tissues and normal bladder tissues were acquired from patients who had undergone surgery and were diagnosed with BCa by pathologists in Shengjing Hospital of China Medical University (Shenyang, China) between 2014 and 2017. Informed consents were signed by each patient before using the tissues. The ethics approval was authorized by the Ethics Committee of Shengjing Hospital of China Medical University.

Reagents

FOXQ1 antibody was purchased from Biorbyt Ltd (host: rabbit; catalog number: orb77456 for western blot and orb53843 for IHC). GAPDH antibody (0411) was purchased from Santa Cruz Biotechnology (host: mouse; catalog: sc-47724). Anti-mouse/rabbit second antibodies were from Invitrogen (Grand Island, NY). miR-140-3p mimic (sequence: UACCACAGGGUAGAACCACGG) and miR-140-3p inhibitor (sequence: AGGCGAAGGAUGACAAAGGGAA) were purchased from Biomics Biotechnologies. The antibodies were kept at -20°C.

Cell culture and transfection

The BCa cell lines (T24 and UMUC3) were purchased from the American Type Culture Collection (ATCC; Manassas, VA) and cultured with DMEM (Invitrogen, Grand Island, NY) containing 10% fetal bovine serum (FBS), penicillin (25 units/ml), 1% L-glutamine, and streptomycin (25 g/ml). The cells were maintained at 37°C in a 5% CO2 humidified atmosphere. Cells were detected and identified as mycoplasma and bacteria free for 3 months following the ATCC's instructions.

Lentivirus packaging

The pWPI; pWPI-FOXQ1-cDNA (F-primer: GTGAGGAATTTCGACATTTAAATTTAATGAAGTTGGAGGTGTTCGT, R-primer: TCCTGCAGCCCGTAGTTTTCAGGCTAGGAGCGTCTC); pLKO.1; pLKO.1-FOXQ1-shRNA (F-primer: CCGGCGAGTACCTCATGGGCAAGTTCTCGAGAACTTGCCCATGAGGTACTCGTTTTTG, R-primer: AATTCAAAAACGAGACCTCATGGGCAAGTTCTCGAGAACTGCCCATGAGGTACTCG); pMD2G envelope plasmid and psPAX2 packaging plasmid were transfected into HEK 293 cells using the standard calcium chloride transfection method. The lentivirus soups were collected after incubation 48 hours and 72 hours and frozen at -80°C for use.

RNA isolation and quantitative real-time polymerase chain reaction (qRT-PCR)

Total RNA was extracted from cells and tissues by Trizol reagent (Invitrogen, Grand Island, NY) according to the manufacturer’s protocols [49]. Complementary DNA was synthesized from RNA reverse transcription by Superscript III transcriptase (Invitrogen, Grand Island, NY). qRT-PCR was conducted using a Bio-Rad CFX96 system. GAPDH was used as an internal control. miRNA was isolated with a PureLink miRNA kit (Invitrogen, Grand Island, NY). RPL32 and U6 were used as endogenous controls. All primers were designed as follows: FOXQ1: F-primer: CTACTCGTACATCGCGCTCA; R-primer: ACCTTGACGAAGCAGTCGTT; miR-140-3p: TACCACAGGGTAGAACCACGG.

Cell proliferation assay

Transfected BCa cells were plated in 24-well plates (2000 cells/well). At day 1, 2, 3, and 4, 50 μL of 10 mg/mL MTT reagent was added to the wells. After 2 hours of incubation, the medium was sucked out and 500 μL DMSO was added. The absorbance was measured at 570 nm.

Cell invasion assay

The upper chamber was coated with 100 μL Matrigel (Corning), which was diluted at 1:20 for T24 cells and 1:30 for UMUC3 cells. The chamber was incubated for 2 hours before cells were plated. The transfected BCa cells were suspended in serum-free medium and seeded into the upper chambers of each transwell (8.0 μM pore size) at a 1x105/mL concentration. Then 750 μL medium containing 10% FBS was added to the bottom chambers. The cells were incubated at 37°C with 5% (v/v) CO2 for 18 hours (for T24 cells) or 24 hours (for UMUC3 cells). The invaded cells on the lower surface were permeabilized with methanol and stained with 0.1% crystal violet in the dark. The stained cells were photographed and counted under a microscope.

Western blot assay

Cells were lysed with RIPA lysis buffer on the ice. Proteins were collected and quantified with BCA analysis. Approximately 50 ug protein was loaded into each lane 2 to 6 days after transfection, separated by a 10% SDS/PAGE gel and transferred onto a PVDF membrane (Millipore, Billerica, MA). After incubation with FOXQ1 antibody (dilution, 1:1000) or GAPDH antibody (dilution, 1:1000) overnight at 4°C, the proteins were then incubated with secondary antibodies (dilution, 1:5000) at room temperature for 1 hour. After samples were washed with TBST, the signals were visualized using a chemiluminescent detection system (ThermoFisher Scientific, Rochester, NY) and analyzed with Image Lab software.

Luciferase reporter assay

Cells were co-transfected with miR-140-3p mimics and plasmid containing 3’-UTR sequences of FOXQ1 wild-type or mutant fragments using the Lipofectamine 3000 reagent (Invitrogen, Carlsbad, CA) according to the manufacturer’s instruction. Dual-luciferase activity was measured 36-48 hours after transfection by a dual-luciferase reporter assay system (Promega, Madison, WI) according to the manufacturer’s instruction.

Xenografts in nude mice

The animal experiments were performed according to the institutional ethics guidelines approved by the Animal Care Committee of China Medical University. Eight nude mice were purchased from Shanghai Laboratory Animal Center Co. Ltd (China). T24 cells were transfected with miR-140-3p mimics or a negative control. T24 cells (2×106) were injected subcutaneously into the posterior flank of the mice (5-6 weeks, 4 mice per group). Tumor sizes were monitored once a week by measuring the length and width with calipers. Tumor volume was calculated by the following formula: volume = (length x width2) x 0.52 [50]. Mice were sacrificed 4 weeks after injection, and tumor weights were measured. qRT-PCR was conducted to quantify miR-140-3p and FOXQ1 expression in xenograft tumors.

Immunohistochemical staining

Each tissue was cut to 5 μm section. The slide was treated for antigen retrieval and incubated with primary antibody. After rinsing with tris-buffered saline, the slide was incubated with secondary antibody, and then incubated with enzyme conjugate horseradish peroxidase. Finally, the slide was counter-stained with hematoxylin.

Statistical analysis

Statistical analyses were performed using SPSS 17.0 (SPSS, Chicago, IL). Data were expressed as mean ± standard deviation (SD) from at least 3 independent experiments. P < 0.05 was considered statistically significant. Data differences were identified by chi-square test or Student’s t-test. Pearson correlation analysis was used for comparing data sets. Kaplan–Meier curves were used for survival analysis.

Supplementary Materials

Supplementary Figure 1

Author Contributions

YKW and YW conceived the experiments; YW, JC and WK performed the experiments; KW and XW analyzed the data; KW wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgments

We thank Yun Cui from Department of Urology, Peking University Third Hospital and Institute of Urology for helping us preparing the manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Funding

This work was supported by the Key Research and Development Joint Program of Liaoning Province (Grant No. 2020JH 2/10300139), the Natural Science Foundation of Liaoning Province (Grant No. 2019-MS-360) and the 345 Talent Project of Shengjing Hospital of China Medical University.

References

  • 1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019; 69:7–34. https://doi.org/10.3322/caac.21551 [PubMed]
  • 2. Zuiverloon TC, Nieuweboer AJ, Vékony H, Kirkels WJ, Bangma CH, Zwarthoff EC. Markers predicting response to bacillus calmette-guérin immunotherapy in high-risk bladder cancer patients: a systematic review. Eur Urol. 2012; 61:128–45. https://doi.org/10.1016/j.eururo.2011.09.026 [PubMed]
  • 3. Witjes JA, Compérat E, Cowan NC, De Santis M, Gakis G, Lebret T, Ribal MJ, Van der Heijden AG, Sherif A, and European Association of Urology. EAU guidelines on muscle-invasive and metastatic bladder cancer: summary of the 2013 guidelines. Eur Urol. 2014; 65:778–92. https://doi.org/10.1016/j.eururo.2013.11.046 [PubMed]
  • 4. Meeks JJ, Bellmunt J, Bochner BH, Clarke NW, Daneshmand S, Galsky MD, Hahn NM, Lerner SP, Mason M, Powles T, Sternberg CN, Sonpavde G. A systematic review of neoadjuvant and adjuvant chemotherapy for muscle-invasive bladder cancer. Eur Urol. 2012; 62:523–33. https://doi.org/10.1016/j.eururo.2012.05.048 [PubMed]
  • 5. Cambier S, Sylvester RJ, Collette L, Gontero P, Brausi MA, van Andel G, Kirkels WJ, Silva FC, Oosterlinck W, Prescott S, Kirkali Z, Powell PH, de Reijke TM, et al. EORTC nomograms and risk groups for predicting recurrence, progression, and disease-specific and overall survival in non-muscle-invasive stage ta-T1 urothelial bladder cancer patients treated with 1-3 years of maintenance bacillus calmette-guérin. Eur Urol. 2016; 69:60–69. https://doi.org/10.1016/j.eururo.2015.06.045 [PubMed]
  • 6. Suer E, Hamidi N, Gokce MI, Gulpinar O, Turkolmez K, Beduk Y, Baltaci S. Significance of second transurethral resection on patient outcomes in muscle-invasive bladder cancer patients treated with bladder-preserving multimodal therapy. World J Urol. 2016; 34:847–51. https://doi.org/10.1007/s00345-015-1710-5 [PubMed]
  • 7. Gupta S, Mahipal A. Role of systemic chemotherapy in urothelial urinary bladder cancer. Cancer Control. 2013; 20:200–10. https://doi.org/10.1177/107327481302000308 [PubMed]
  • 8. Djebali S, Davis CA, Merkel A, Dobin A, Lassmann T, Mortazavi A, Tanzer A, Lagarde J, Lin W, Schlesinger F, Xue C, Marinov GK, Khatun J, et al. Landscape of transcription in human cells. Nature. 2012; 489:101–08. https://doi.org/10.1038/nature11233 [PubMed]
  • 9. Hecker N, Stephan C, Mollenkopf HJ, Jung K, Preissner R, Meyer HA. A new algorithm for integrated analysis of miRNA-mRNA interactions based on individual classification reveals insights into bladder cancer. PLoS One. 2013; 8:e64543. https://doi.org/10.1371/journal.pone.0064543 [PubMed]
  • 10. Salem O, Erdem N, Jung J, Münstermann E, Wörner A, Wilhelm H, Wiemann S, Körner C. The highly expressed 5’isomiR of hsa-miR-140-3p contributes to the tumor-suppressive effects of miR-140 by reducing breast cancer proliferation and migration. BMC Genomics. 2016; 17:566. https://doi.org/10.1186/s12864-016-2869-x [PubMed]
  • 11. Dong W, Yao C, Teng X, Chai J, Yang X, Li B. MiR-140-3p suppressed cell growth and invasion by downregulating the expression of ATP8A1 in non-small cell lung cancer. Tumour Biol. 2016; 37:2973–85. https://doi.org/10.1007/s13277-015-3452-9 [PubMed]
  • 12. Lian L, Zhang D, Wang Q, Yang N, Qu L. The inhibitory effects of gga-miR-199-3p, gga-miR-140-3p, and gga-miR-221-5p in marek’s disease tumorigenesis. Poult Sci. 2015; 94:2131–35. https://doi.org/10.3382/ps/pev175 [PubMed]
  • 13. Luo Q, Wang CQ, Yang LY, Gao XM, Sun HT, Zhang Y, Zhang KL, Zhu Y, Zheng Y, Sheng YY, Lu L, Jia HL, Yu WQ, et al. FOXQ1/NDRG1 axis exacerbates hepatocellular carcinoma initiation via enhancing crosstalk between fibroblasts and tumor cells. Cancer Lett. 2018; 417:21–34. https://doi.org/10.1016/j.canlet.2017.12.021 [PubMed]
  • 14. Feng A, Yuan X, Li X. MicroRNA-345 inhibits metastasis and epithelial-mesenchymal transition of gastric cancer by targeting FOXQ1. Oncol Rep. 2017; 38:2752–60. https://doi.org/10.3892/or.2017.6001 [PubMed]
  • 15. Liu JY, Wu XY, Wu GN, Liu FK, Yao XQ. FOXQ1 promotes cancer metastasis by PI3K/AKT signaling regulation in colorectal carcinoma. Am J Transl Res. 2017; 9:2207–18. [PubMed]
  • 16. Zhang X, Wang L, Wang Y, Shi S, Zhu H, Xiao F, Yang J, Yang A, Hao X. Inhibition of FOXQ1 induces apoptosis and suppresses proliferation in prostate cancer cells by controlling BCL11A/MDM2 expression. Oncol Rep. 2016; 36:2349–56. https://doi.org/10.3892/or.2016.5018 [PubMed]
  • 17. Xiao B, Liu H, Gu Z, Ji C. Expression of microRNA-133 inhibits epithelial-mesenchymal transition in lung cancer cells by directly targeting FOXQ1. Arch Bronconeumol. 2016; 52:505–11. https://doi.org/10.1016/j.arbres.2015.10.016 [PubMed]
  • 18. Zhang J, Li W, Dai S, Tai X, Jia J, Guo X. FOXQ1 is overexpressed in laryngeal carcinoma and affects cell growth, cell cycle progression and cell invasion. Oncol Lett. 2015; 10:2499–504. https://doi.org/10.3892/ol.2015.3530 [PubMed]
  • 19. Pei Y, Wang P, Liu H, He F, Ming L. FOXQ1 promotes esophageal cancer proliferation and metastasis by negatively modulating CDH1. Biomed Pharmacother. 2015; 74:89–94. https://doi.org/10.1016/j.biopha.2015.07.010 [PubMed]
  • 20. Zhan HX, Xu JW, Wang L, Wu D, Zhang GY, Hu SY. FoxQ1 is a novel molecular target for pancreatic cancer and is associated with poor prognosis. Curr Mol Med. 2015; 15:469–77. https://doi.org/10.2174/1566524015666150630125247 [PubMed]
  • 21. Meng F, Speyer CL, Zhang B, Zhao Y, Chen W, Gorski DH, Miller FR, Wu G. PDGFRα and β play critical roles in mediating Foxq1-driven breast cancer stemness and chemoresistance. Cancer Res. 2015; 75:584–93. https://doi.org/10.1158/0008-5472.CAN-13-3029 [PubMed]
  • 22. Zhu Z, Zhu Z, Pang Z, Xing Y, Wan F, Lan D, Wang H. Short hairpin RNA targeting FOXQ1 inhibits invasion and metastasis via the reversal of epithelial-mesenchymal transition in bladder cancer. Int J Oncol. 2013; 42:1271–78. https://doi.org/10.3892/ijo.2013.1807 [PubMed]
  • 23. Kaestner KH. The hepatocyte nuclear factor 3 (HNF3 or FOXA) family in metabolism. Trends Endocrinol Metab. 2000; 11:281–85. https://doi.org/10.1016/s1043-2760(00)00271-x [PubMed]
  • 24. Ang SL, Rossant J. HNF-3 beta is essential for node and notochord formation in mouse development. Cell. 1994; 78:561–74. https://doi.org/10.1016/0092-8674(94)90522-3 [PubMed]
  • 25. Medema RH, Kops GJ, Bos JL, Burgering BM. AFX-like forkhead transcription factors mediate cell-cycle regulation by ras and PKB through p27kip1. Nature. 2000; 404:782–87. https://doi.org/10.1038/35008115 [PubMed]
  • 26. Costa RH, Grayson DR, Darnell JE Jr. Multiple hepatocyte-enriched nuclear factors function in the regulation of transthyretin and alpha 1-antitrypsin genes. Mol Cell Biol. 1989; 9:1415–25. https://doi.org/10.1128/mcb.9.4.1415 [PubMed]
  • 27. Ortiz L, Zannini M, Di Lauro R, Santisteban P. Transcriptional control of the forkhead thyroid transcription factor TTF-2 by thyrotropin, insulin, and insulin-like growth factor I. J Biol Chem. 1997; 272:23334–39. https://doi.org/10.1074/jbc.272.37.23334 [PubMed]
  • 28. Li J, Vogt PK. The retroviral oncogene qin belongs to the transcription factor family that includes the homeotic gene fork head. Proc Natl Acad Sci USA. 1993; 90:4490–94. https://doi.org/10.1073/pnas.90.10.4490 [PubMed]
  • 29. Bieller A, Pasche B, Frank S, Gläser B, Kunz J, Witt K, Zoll B. Isolation and characterization of the human forkhead gene FOXQ1. DNA Cell Biol. 2001; 20:555–61. https://doi.org/10.1089/104454901317094963 [PubMed]
  • 30. Jiao D, Li Z, Zhu M, Wang Y, Wu G, Han X. LncRNA MALAT1 promotes tumor growth and metastasis by targeting miR-124/foxq1 in bladder transitional cell carcinoma (BTCC). Am J Cancer Res. 2018; 8:748–60. [PubMed]
  • 31. Peng XH, Huang HR, Lu J, Liu X, Zhao FP, Zhang B, Lin SX, Wang L, Chen HH, Xu X, Wang F, Li XP. MiR-124 suppresses tumor growth and metastasis by targeting Foxq1 in nasopharyngeal carcinoma. Mol Cancer. 2014; 13:186. https://doi.org/10.1186/1476-4598-13-186 [PubMed]
  • 32. Zhang Z, Ma J, Luan G, Kang L, Su Y, He Y, Luan F. MiR-506 suppresses tumor proliferation and invasion by targeting FOXQ1 in nasopharyngeal carcinoma. PLoS One. 2015; 10:e0122851. https://doi.org/10.1371/journal.pone.0122851 [PubMed]
  • 33. Cui Z, Zhao Y. microRNA-342-3p targets FOXQ1 to suppress the aggressive phenotype of nasopharyngeal carcinoma cells. BMC Cancer. 2019; 19:104. https://doi.org/10.1186/s12885-018-5225-5 [PubMed]
  • 34. Xiang XJ, Deng J, Liu YW, Wan LY, Feng M, Chen J, Xiong JP. MiR-1271 inhibits cell proliferation, invasion and EMT in gastric cancer by targeting FOXQ1. Cell Physiol Biochem. 2015; 36:1382–94. https://doi.org/10.1159/000430304 [PubMed]
  • 35. Xu J, You Q, Wei Z, Fu H, Zhang Y, Hu Z, Cai Q. miR-519 inhibits epithelial-mesenchymal transition and biologic behavior of gastric cancer cells by down-regulating FOXQ1. Int J Clin Exp Pathol. 2020; 13:425–36. [PubMed]
  • 36. Vishnubalaji R, Hamam R, Yue S, Al-Obeed O, Kassem M, Liu FF, Aldahmash A, Alajez NM. MicroRNA-320 suppresses colorectal cancer by targeting SOX4, FOXM1, and FOXQ1. Oncotarget. 2016; 7:35789–802. https://doi.org/10.18632/oncotarget.8937 [PubMed]
  • 37. Weng W, Okugawa Y, Toden S, Toiyama Y, Kusunoki M, Goel A. FOXM1 and FOXQ1 are promising prognostic biomarkers and novel targets of tumor-suppressive miR-342 in human colorectal cancer. Clin Cancer Res. 2016; 22:4947–57. https://doi.org/10.1158/1078-0432.CCR-16-0360 [PubMed]
  • 38. Liu Z, Qin Y, Dong S, Chen X, Huo Z, Zhen Z. Overexpression of miR-106a enhances oxaliplatin sensitivity of colorectal cancer through regulation of FOXQ1. Oncol Lett. 2020; 19:663–70. https://doi.org/10.3892/ol.2019.11151 [PubMed]
  • 39. Xie WB, Liang LH, Wu KG, Wang LX, He X, Song C, Wang YQ, Li YH. MiR-140 expression regulates cell proliferation and targets PD-L1 in NSCLC. Cell Physiol Biochem. 2018; 46:654–63. https://doi.org/10.1159/000488634 [PubMed]
  • 40. Fang Z, Yin S, Sun R, Zhang S, Fu M, Wu Y, Zhang T, Khaliq J, Li Y. miR-140-5p suppresses the proliferation, migration and invasion of gastric cancer by regulating YES1. Mol Cancer. 2017; 16:139. https://doi.org/10.1186/s12943-017-0708-6 [PubMed]
  • 41. Lv J, Fan HX, Zhao XP, Lv P, Fan JY, Zhang Y, Liu M, Tang H. Long non-coding RNA Unigene56159 promotes epithelial-mesenchymal transition by acting as a ceRNA of miR-140-5p in hepatocellular carcinoma cells. Cancer Lett. 2016; 382:166–75. https://doi.org/10.1016/j.canlet.2016.08.029 [PubMed]
  • 42. Wang M, Du M, Ma L, Chu H, Lv Q, Ye D, Guo J, Gu C, Xia G, Zhu Y, Ding Q, Yuan L, Fu G, et al. A functional variant in TP63 at 3q28 associated with bladder cancer risk by creating an miR-140-5p binding site. Int J Cancer. 2016; 139:65–74. https://doi.org/10.1002/ijc.29978 [PubMed]
  • 43. Li Q, Yao Y, Eades G, Liu Z, Zhang Y, Zhou Q. Downregulation of miR-140 promotes cancer stem cell formation in basal-like early stage breast cancer. Oncogene. 2014; 33:2589–600. https://doi.org/10.1038/onc.2013.226 [PubMed]
  • 44. Li W, Jiang G, Zhou J, Wang H, Gong Z, Zhang Z, Min K, Zhu H, Tan Y. Down-regulation of miR-140 induces EMT and promotes invasion by targeting slug in esophageal cancer. Cell Physiol Biochem. 2014; 34:1466–76. https://doi.org/10.1159/000366351 [PubMed]
  • 45. Zhang Y, Eades G, Yao Y, Li Q, Zhou Q. Estrogen receptor α signaling regulates breast tumor-initiating cells by down-regulating miR-140 which targets the transcription factor SOX2. J Biol Chem. 2012; 287:41514–22. https://doi.org/10.1074/jbc.M112.404871 [PubMed]
  • 46. Yuan Y, Shen Y, Xue L, Fan H. miR-140 suppresses tumor growth and metastasis of non-small cell lung cancer by targeting insulin-like growth factor 1 receptor. PLoS One. 2013; 8:e73604. https://doi.org/10.1371/journal.pone.0073604 [PubMed]
  • 47. Rothschild SI. microRNA therapies in cancer. Mol Cell Ther. 2014; 2:7. https://doi.org/10.1186/2052-8426-2-7 [PubMed]
  • 48. Flamini V, Jiang WG, Cui Y. Therapeutic role of MiR-140-5p for the treatment of non-small cell lung cancer. Anticancer Res. 2017; 37:4319–27. https://doi.org/10.21873/anticanres.11825 [PubMed]
  • 49. Jung U, Jiang X, Kaufmann SH, Patzel V. A universal TaqMan-based RT-PCR protocol for cost-efficient detection of small noncoding RNA. RNA. 2013; 19:1864–73. https://doi.org/10.1261/rna.040501.113 [PubMed]
  • 50. Tomayko MM, Reynolds CP. Determination of subcutaneous tumor size in athymic (nude) mice. Cancer Chemother Pharmacol. 1989; 24:148–54. https://doi.org/10.1007/BF00300234 [PubMed]