ALDH1A3–Linc00284 Axis Mediates the Invasion of Colorectal Cancer by Targeting TGFβ Signaling via Sponging miR-361-5p

ALDH1A3 and Linc00284 involve in colorectal cancer (CRC) development; however, the regulatory mechanism is still unclear. In this study, we collected clinicopathological characteristics and tissue samples from 73 CRC patients to analyze the expression of ALDH1A3, Linc00284, TGFβ signaling and miR-361-5p using qPCR, Western blotting, and ELISA. Multiple CRC cell lines were evaluated in this study, and the highest level of ALDH1A3 was observed in SW480 cells. To investigate the regulatory mechanism, RIP and luciferase assays were used to validate the interaction between Linc00284, miR-361-5p, and TGFβ. Proliferation, viability, migration, and invasion assays were performed to profile the effects of the ALDH1A3–Linc00284 axis in CRC cell functions, which was upregulated in CRC tissues. Knockdown ALDH1A3 or Linc00284 significantly reduced TGFβ expression and suppressed the EMT process, while overexpression had opposite effects. miR-361-5p targeted TGFβ directly, which negatively correlated with ALDH1A3–Linc00284 expression and CRC progression. Mechanistically, upregulation of ALDH1A3–Linc00284 promotes colorectal cancer invasion and migration by regulating miR-361-5p/TGFβ signaling pathway. Dysregulation of the ALDH1A3–Linc00284-miR-361-5p-TGFβ axis causes CRC invasion, which might provide a new insight into the treatment of CRC.


Introduction
Colorectal cancer is the second highest mortality tumor worldwide [1,2]. By 2030, the number of global cases of CRC is expected to increase by 60%, with more than 2.25 million new cases and 1.15 million deaths from CRC cancer each year [3,4]. In the past few decades, CRC has become one of the most common cancers, the incidence has increased from 1.0% annually to 2.4% since 1974 in the USA [5]. In most Europe, the incidence of CRC is increasing every year, ranging from 0.4% to 3.6% [6]. The existing therapies for CRC include surgery, radiotherapy, chemotherapy, targeted drug therapy, and so on. Immunotherapy is an encouraging weapon to treat CRC, but so far it has only succeeded in a small proportion of CRC [7]. The prognosis of metastatic CRC is worse, with an overall survival time of only about 30 months [8]. Therefore, it is very urgent to develop new effective strategies to control CRC, especially metastatic CRC.
Cumulative evidence indicates that ALDH1A3 and Linc00284 play an important role in the progression of CRC [9][10][11][12]. ALDH1A3 is a known marker of cancer stem cells which has been shown to be important for the proliferation, migration, and maintenance of the mesenchymal cancer stem cell phenotype [13]. Durinikova et al. demonstrated that ALDH1A3 was increased in CRC tissues, which promotes spontaneous metastasis formation and associates with acquired chemoresistance of colorectal cancer [9]. Long noncoding RNAs (lncRNAs) Linc00284 has been reported to be involved in the initiation and progression of many cancers, including oral squamous cell carcinoma, ovarian carcinoma, papillary thyroid cancer, lung cancer, CRC, and so on [12,[14][15][16][17]. Our recent study demonstrated that Linc00284 exhibits oncogenic function and promotes the progression of CRC through upregulating the expression of c-Met [12]. However, the underlying molecular mechanisms of ALDH1A3 and Linc00284 on the growth and metastasis of CRC are still unclear. Here, we investigated the role of the ALDH1A3-Linc00284 signal in CRC progression and in the tumor microenvironment through in-depth analysis of clinical data and in vitro experiments. This study proves that ALDH1A3-Linc00284 upregulates TGFβ signaling through miR-361-5p, and then promotes the epithelial-mesenchymal transition (EMT) process and CRC tumor metastasis. This research will provide a key theoretical basis for the CRC treatment and drug development.

Materials and Methods
2.1. Human Study. CRC tissues and the paired adjacent nontumor samples were derived from 73 CRC patients, which are consistent with the patients and samples in our previous study [12], including the patient's information, inclusion and exclusion criteria, and metastatic and relapse records. The TNM stages of CRC in the above patients were classified based on the American Joint Committee on Cancer (AJCC) tumor, lymph node, metastasis (TNM) system. The study was obtained the written informed consent from all participants and approved by the ethics committee of the First Affiliated Hospital of Fujian Medical University.
2.4. qRT-PCR. Total RNA from homogenize tissue samples and lysed cells was extracted using Trizol reagent (Invitrogen, Waltham, MA) according to the protocol provided by the manufacturer. TaqMan® Reverse Transcription Reagents and random primers were used to synthesize the cDNA. Ste-pOne™ PCR System (Thermo Fisher Scientific Inc) was used to perform qPCR reaction. The expression of target genes was normalized by GAPDH or U6. All primers used in this study were listed in Table 1.

Forward primer
Reverse primer International Journal of Genomics 2.5. Western Blotting. Proteins were extracted from homogenized tissue samples by using RIPA buffer containing PMSF and protease inhibitor cocktails (KeyGEN BioTECH, Nanjing, China). Protein levels of ALDH1A3 and β-Actin were determined using Western blotting as described previously [12]. The primary antibodies and information were listed below: β-Actin (Proteintech, 60008- WT, and the corresponding mutated fragments were cloned into an overexpression vector (pcDNA3.1). The above constructs or control vectors were co-transfected with the luciferase reporter into HEK293T cells using Lipofectamine™ 2000 (Invitrogen). Two days after transfection, HEK293T cells were lysed and luciferase activity was determined by the Dual-Luciferase® Reporter (DLR™) Assay System (Promega) following the manufacturer's instruction.

Migration and Invasion Assays.
To analyze cell migration, 5 × 10 6 SW480 cells were seeded in 6-well culture plate. The cell monolayer was scratched by a sterile tip when they reached 90% confluence. The wound healing rate was imaged and calculated every 24 hours. CRC cell invasion was performed by trans-well assay as described previously [12].
2.11. Statistical Analysis. SPSS 22.0 statistical software (SPSS Inc.) was used to analyze the data in this study. A two-sided unpaired t test was used to compare differences between the two groups. One-way analysis of variance (ANOVA) and post-hoc least significant difference (LSD) tests were used to analyze the differences among multiple groups. The sig-nificant difference or correlation among groups was analyzed by two-way ANOVA and/or two-sided paired t-test, and Pearson's correlation analysis. Kaplan-Meier curves with the log rank test was used to calculate the overall survival (OS) of CRC patients. All data are presented as mean ± SD.

ALDH1A3 Upregulation Is Associated with Poor
Prognosis of CRC Patients. The expression level of ALDH1A3 in CRC tissues (n = 73) was measured by qRT-PCR, in comparison with paired adjacent normal tissues, ALDH1A3 mRNA was significantly increased in tumor samples ( Figure 1(a)). Western blot results showed markedly elevated protein level of ALDH1A3 in CRC tissues compared to that in controls ( Figure 1(b)). Moreover, ALDH1A3 expression was higher in tumors of patients with metastatic CRC (n = 43) than those without metastasis (n = 30) ( Figure 1(c)). In addition, patients with relapse CRC (n = 49) had higher ALDH1A3 expression compared with patients without relapse (n = 24) (Figure 1(d)). Besides, higher ALDH1A3 expression positively correlated with the TNM stages of CRC and poor overall survival. ALDH1A3 expression was higher in tumors with higher CRC TNM grade ( Figure 1(e)). Patients with higher ALDH1A3 expression in had a shorter survival time compared to patients with lower expression of ALDH1A3 (Figure 1(f)).

Knockdown of ALDH1A3 Represses CRC Proliferation and Metastasis In Vitro.
To investigate the role of ALDH1A3 in CRC, we examined the mRNA level of ALDH1A3 in six CRC cell lines (SW480, HCT116, LS174T, DLD-1, HCT15, and SW620), and found that ALDH1A3 was highest expressed in SW480 cells, which was used in further studies (Figure 2(a)). ALDH1A3 knockdown SW480 cell line was established through shRNA method. Among three shRNAs, shRNA TRCN0000027144 showed the highest efficacy and was used for silencing ALDH1A3 in following experiments ( Figure 2(b)). CCK-8 assay indicated that ALDH1A3 knockdown reduced cell viability of CRC cells significantly (Figure 2(c)). We performed flow cytometry analysis using Ki 67 staining to examine the cancer cell proliferation with and without ALDH1A3. Upon ALDH1A3 knockdown, Ki67 positive SW480 cell was decreased markedly relative to that in control group (Figure 2(d)). In addition, migration and trans-well assays showed that ALDH1A3 silencing decreased the ability of wound healing and invasion significantly in SW480 cells (Figures 2(e) and 2(f)).

ALDH1A3 Mediates the Expression of Linc00284
In Vitro and In Vivo. Our previous study indicated that Linc00284 plays an important role in CRC progression [12], qPCR analysis showed that knockdown of ALDH1A3 significantly reduced the expression level of Linc00284 (Figure 3(a)), while overexpression of ALDH1A3 upregulated Linc00284 expression markedly (Figure 3(b)). Interestingly, ALDH1A3 expression was correlated with Linc00284 expression positively (Figure 3(c)). The above data indicated that   International Journal of Genomics of EMT-related genes, as indicated by deceased mRNA and protein expression of E-cadherin, and increased mRNA and protein expression of Vimentin and N-cadherin (Figures 4(j)-4(p)), which suggests the inhibition of EMT process by overexpression of ALDH1A3 or Linc00284. TGFβ signaling is a key protein factor that mediates the process of EMT [19]. Both qPCR and ELISA data showed that knockdown ALDH1A3 or Linc00284 significantly reduced the expression level of TGFβ, while overexpression had the opposite effect (Figures 4(h)-4(i) and 4(q)-4(r)). Moreover, overexpression of Linc00284 could reverse the downregulation of ALDH1A3 knockdown on TGFβ mRNA and protein (Figures 4(s)-4(t)). These results suggested that ALDH1A3-Linc00284 might promote EMT process by regulating TGFβ signal, and then promote CRC tumor metastasis.     International Journal of Genomics database and identified the microRNA candidate that can bind to Linc00284 and TGFβ. Sequence analysis showed that miR-361-5p has the complementary sequences that match Linc00284 and TGFβ, respectively (Figures 5(a) and 5(g)).

Inhibition of miR-361-5p
Rescues the Effect ALDH1A3-Linc00284 on CRC Progression. In human CRC tissues, the expression of miR-361-5p was decreased significantly relative to that in adjacent normal tissues (Figure 6(a)). Interestingly, miR-361-5p was reduced in distant metastasis ( Figure 6(b)), relapse (Figure 6(c)), and high TNM staging  (Figure 6(d)). Furthermore, miR-361-5p expression negatively correlated with the expression of both ALDH1A3 and Linc00284 in CRC tissues (Figures 6(e) and 6(f)). Next, we silenced miR-361-5p in ALDH1A3-or Linc00284knockdown SW480 cells and examined the proliferation and metastasis of CRC cells. As expected, inhibition of miR-361-5p reversed the effects of ALDH1A3-or Linc00284knockdown on cell viability (Figure 7(a)). Identical effects of miR-361-5p inhibition on CRC cell migration and invasion were observed by wound healing assay and trans-well assay (Figures 7(b) and 7(c)). The above data further indicated that ALDH1A3-Linc00284 mediates CRC invasion through regulating the expression miR-361-5p.

Discussion
Long noncoding RNAs play an important role in the tumor initiation and progression [11,[20][21][22][23][24]. Multiple studies showed higher Linc00284 expression in breast cancer and liver cancer, which can promote the tumor progression and associated with poor overall survival [11,25,26]. Our recent study demonstrated that the upregulation of Linc00284 in tumor samples of CRC patients; in addition, Linc00284 expression positively correlated with metastasis, recurrence, and poor survival [12]. However, the underlying mechanism of Linc00284-mediated CRC progression is still unclear. In this study, we revealed that Linc00284 was upregulated by ALDH1A3 in CRC tissues and cells, and the ALDH1A3-Linc00284 axis promoted the invasion of CRC through activation of TGFβ signaling and downstream EMT process. Mechanistically, upregulation of ALDH1A3-Linc00284 promotes colorectal cancer invasion and migration by regulating the miR-361-5p/TGFβ signaling pathway, which might provide a new insight into the treatment of CRC.
EMT not only plays a critical role in tumor stemness, proliferation, migration, and invasion but also associates with the tumor microenvironment to induce immunosuppression and cause resistance to therapy [27,28]. TGFβ signaling can trigger EMT when cells are in the certain disease microenvironment, such as fibrosis and cancer [28,29]. In CRC, TGFβ and downstream factors are mobilized by integrin, which promotes the EMT processed in both cell and animal models [30]. On the contrary, decreasing TGFβ expression and activity in the tumor microenvironment leads to potent immune responses of CRC in rodent models [31,32]. In line with the above findings, knockdown ALDH1A3-Linc00284 axis significantly affects the expression level of EMT-related genes E-cadherin, N-cadherin, and Vimentin, and inhibits the EMT process, while overexpression of each of them has opposite effects. Moreover, changing the expression levels of the ALDH1A3-Linc00284 axis also affects the viability and invasion of CRC cells. Importantly, we found that higher expression of ALDH1A3-Linc00284 positively correlates with the TNM staging and poor survival of CRC patients. These findings suggested that ALDH1A3-Linc00284 promotes the EMT process by regulating the TGFβ signal, and then promotes CRC tumor metastasis. microRNAs are small noncoding RNAs that function as either tumor suppressors or oncogenes under certain conditions [33]. More and more studies have identified miRNAs as potential biomarkers for human cancer diagnosis, prognosis and treatment targets. The miR-361-5p has been reported widely expressed in many tissues. Higher expression of miR-361-5p indicates better prognosis of breast cancer patients [34]. In gastric cancer, miR-361-5p can suppress chemoresistance of SGC-7901 and MKN-28 cells through inhibition of the expression of FOXM1 [35]. Here, we identified that TGFβ is a direct target of miR-361-5p in CRC cells, silencing miR-361-5p can induce the expression of TGFβ and promote the EMT process. Interestingly, sequence analysis revealed that miR-361-5p can be bound by Linc00284 in SW480 cells. In tumor samples of CRC patients, miR-361- 10 International Journal of Genomics 5p expression is negatively correlated with the expression of both ALDH1A3 and Linc00284. More importantly, inhibition of miR-361-5p can rescue the effect ALDH1A3 or Linc00284 knockdown in CRC cells. Our findings indicate that the ALDH1A3-Linc00284 axis mediates the progression of CRC by targeting TGFβ signaling via sponging miR-361-5p. Animal models need to be established to further investigate and confirm the regulatory effect of ALDH1A3-Linc00284-miR-361-5p in the CRC tumor microenvironment in the future studies.

Conclusions
Our findings indicate that the ALDH1A3-Linc00284 axis mediates the progression of CRC by targeting TGFβ signaling via sponging miR-361-5p in CRC cells, providing new insight into the pathogenesis and treatment of colorectal cancer.

Data Availability
The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.