Next Article in Journal
Understanding the Role of Radio-Sensitizing Nanoparticles in Enhancing Pathologic Response in Soft Tissue Sarcomas
Next Article in Special Issue
RNA-Independent Regulatory Functions of lncRNA in Complex Disease
Previous Article in Journal
Quality of Life and Independent Factors Associated with Poor Digestive Function after Ivor Lewis Esophagectomy
Previous Article in Special Issue
Functional In Vivo Screening Identifies microRNAs Regulating Metastatic Dissemination of Prostate Cancer Cells to Bone Marrow
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Tumor-Suppressive miR-139-3p-Regulated Genes: TRIP13 as a Therapeutic Target in Lung Adenocarcinoma

1
Department of Pulmonary Medicine, Graduate School of Medical and Dental Sciences, Kagoshima University, Kagoshima 890-8520, Japan
2
Department of Functional Genomics, Chiba University Graduate School of Medicine, Chuo-ku, Chiba 260-8670, Japan
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Cancers 2023, 15(23), 5571; https://doi.org/10.3390/cancers15235571
Submission received: 25 October 2023 / Revised: 17 November 2023 / Accepted: 22 November 2023 / Published: 24 November 2023

Abstract

:

Simple Summary

Based on the miRNA expression signature of LUAD, we focused on miR-139-3p, a passenger strand, and clarified its tumor-suppressive function in lung adenocarcinoma (LUAD) cells. A total of the 21 target genes (KRT80, CENPM, SPC24, ORC1, MYEOV, TRIP13, GPX8, ARHGEF39, MKI67, KIF18B, CHAF1B, CP, GPRIN1, UCK2, CHEK1, HELLS, CTSV, FAM111B, SLC16A3, MELK, and CENPF) were identified as miR-139-3p targets in LUAD, and the expression of these genes was as independent prognostic factor for patient survival. Moreover, inhibition of TRIP13 using a specific inhibitor (DCZ0415) enhanced the sensitivity of LUAD cells to anticancer drugs.

Abstract

Analyses of our microRNA (miRNA) expression signature combined with The Cancer Genome Atlas (TCGA) data revealed that both strands of pre-miR-139 (miR-139-5p, the guide strand, and miR-139-3p, the passenger strand) are significantly downregulated in lung adenocarcinoma (LUAD) clinical specimens. Functional analyses of LUAD cells ectopically expressing miR-139-3p showed significant suppression of their aggressiveness (e.g., cancer cell proliferation, migration, and invasion). The involvement of the passenger strand, miR-139-3p, in LUAD pathogenesis, is an interesting finding contributing to the elucidation of unknown molecular networks in LUAD. Of 1108 genes identified as miR-139-3p targets in LUAD cells, 21 were significantly upregulated in LUAD tissues according to TCGA analysis, and their high expression negatively affected the prognosis of LUAD patients. We focused on thyroid hormone receptor interactor 13 (TRIP13) and investigated its cancer-promoting functions in LUAD cells. Luciferase assays showed that miR-139-3p directly regulated TRIP13. siRNA-mediated TRIP13 knockdown and TRIP13 inhibition by a specific inhibitor (DCZ0415) attenuated the malignant transformation of LUAD cells. Interestingly, when used in combination with anticancer drugs (cisplatin and carboplatin), DCZ0415 exerted synergistic effects on cell proliferation suppression. Identifying the molecular pathways regulated by tumor-suppressive miRNAs (including passenger strands) may aid in the discovery of diagnostic markers and therapeutic targets for LUAD.

1. Introduction

Lung cancer is the leading cause of cancer-related death worldwide, with approximately 2.3 million new cases of lung cancer and 1.8 million lung cancer-related deaths each year [1]. Patient prognosis is extremely poor, with a 5-year overall survival rate of only 5% for patients with metastatic cancer. Even in patients whose tumors are localized to the lungs, the overall survival rate ranges from 33% to 60% [2]. According to the histological classification of lung cancer, approximately 85% of all lung cancers are non-small-cell lung cancer (NSCLC), the majority of which are lung adenocarcinoma (LUAD) [3]. The treatment methods for LUAD include surgery, radiotherapy, chemotherapy, molecular targeted therapy, and immunotherapy, and the treatment strategy is determined by the disease stage and the presence or absence of driver gene mutations [3]. Molecularly targeted drugs that address genetic mutations in driver genes are being developed rapidly, and many patients will benefit from these drugs [3]. However, in the past, cytotoxic anticancer drugs were the mainstay of drug therapy for advanced-stage LUAD in which molecularly targeted drugs were not applicable. Recently, the effectiveness of immunotherapy has been confirmed, and immunotherapy alone or in combination with anticancer drugs is being used to treat LUAD [4]. It is essential to develop new treatment regimens that increase the effectiveness of anticancer drugs for advanced-stage patients without driver gene mutations.
In the post-genomic era, it has been discovered that the human genome contains an extremely large number of functional RNAs that do not encode proteins. We realized the importance of clarifying their functions in normal and diseased cells [5]. microRNAs (miRNAs) are short (~22 nucleotides) non-coding RNAs that act mainly as gene expression regulators at the post-transcriptional level in a sequence-dependent manner [6]. A single miRNA controls the expression of many different types of genes, and miRNAs play important roles in physiological cellular processes [7].
Numerous studies have shown that the aberrant expression of miRNAs is closely involved in cancer cell malignant transformation, e.g., proliferation, invasion, metastasis, and drug resistance [8]. Traditionally, it was thought that only the guide strand of miRNAs derived from pre-miRNAs was functional in cells, and therefore, cancer research has been focused on the guide strands. Our recent studies revealed that the passenger strand of miRNAs derived from pre-miRNAs is intricately involved in the molecular pathogenesis of human cancers, including lung cancer [9,10]. We can uncover new molecular pathways in cancer cells by examining the molecules regulated by both strands of miRNAs derived from pre-miRNAs.
Recently, we created a miRNA expression signature based on LUAD tissues using RNA sequencing [10]. Analysis of this signature showed that both the guide and passenger strands of pre-miR-139 (miR-139-5p and miR-139-3p, respectively) were downregulated in LUAD tissues. Focusing on miR-139-3p, we have identified its tumor-suppressive roles and cancer-promoting target genes in various cancer types (e.g., bladder cancer, renal cell carcinoma, oral cancer, and colorectal cancer) [11,12,13,14]. A unique feature of miRNAs is that the target genes of a specific miRNA differ depending on the cell. Our previous studies confirmed that the same tumor-suppressive miRNA regulates different genes depending on the cancer type [11,12,13,14]. The aim of this study was to confirm the tumor-suppressive function of miR-139-3p and to clarify the molecular networks controlled by miR-139-3p in LUAD cells. Ectopic expression of miR-139-3p in cancer cells significantly blocked proliferation and induced cell cycle arrest and apoptosis. Regarding miR-139-3p target genes, thyroid hormone receptor interactor 13 (TRIP13) was found to be directly regulated by tumor-suppressive miR-139-3p in LUAD cells. Aberrant expression of TRIP13 facilitated the malignant transformation of LUAD cells. Importantly, a TRIP13-specific inhibitor (DCZ0415) exerted synergistic effects on the suppression of cell proliferation when used in combination with various anticancer drugs (e.g., cisplatin and carboplatin).
Exploring the oncogenic networks controlled by tumor-suppressive miRNAs, including the passenger strands, will facilitate the identification of therapeutic target molecules for LUAD.

2. Materials and Methods

2.1. miRNA Expression and Clinical Significance in Patients with LUAD Using an In Silico Database Analysis

We recently created a miRNA expression signature based on LUAD clinical specimens using RNA sequencing [10]. Using this signature, we selected miRNAs with suppressed expression in LUAD tissues. The expression of miRNAs and genes in LUAD clinical tissues was analyzed using the following databases: The Cancer Genome Atlas (TCGA) (https://www.cancer.gov/tcga, accessed on 17 January 2023), FIREBROWSE (http://firebrowse.org/, accessed on 17 January 2023), and Genomic Data Commons Data Portal (https://portal.gdc.cancer.gov/, accessed on 17 January 2023). Overall survival data were obtained from cBioPortal (https://www.cbioportal.org/, accessed on 17 January 2023) and OncoLnc (http://www.oncolnc.org/) (data downloaded on 17 January 2023).

2.2. Functional Assay of miRNAs and miRNA Target Genes in LUAD Cells

Two LUAD cell lines, A549 and H1299, were used in the functional assays in this study (American Type Culture Collection, Manassas, VA, USA). Small RNAs (miRNAs and siRNAs) were transfected into LUAD cell lines, and cell proliferation, migration, and invasion were evaluated and compared with those of control cells. The transfection procedures of miRNAs and siRNAs were described in our previous studies [10,15,16,17]. All miRNA precursors were transfected at 10 nM, and all siRNAs were transfected at 5 nM into A549 and H1299 cell lines using RNAiMAX (Invitrogen, Carlsbad, CA, USA). Mock was a group without precursors or siRNAs. Cell cycle analysis was performed using a flow cytometer (BD FACSCelestaTM Flow Cytometer, BD Biosciences, Franklin Lakes, NJ, USA). Details of the cell functional assays are described in our previous paper [10,15,16,17]. The reagents used for these analyses are listed in Table S1.

2.3. Identification of Oncogenic Targets Regulated by miR-139-3p in LUAD Cells

To identify putative gene targets controlled by miR-139-3p in LUAD, we used TargetScanHuman v8.0 (https://www.targetscan.org/vert_80/, accessed on 24 May 2023) and a gene expression profile from the GEO database (GEO accession number: GSE242241). This expression profile comprises genes whose expression was altered after miR-139-3p transfection into A549 cells.
We used GeneCodis4 software to infer the molecular functions of the miR-139-3p target genes [18]. Gene set enrichment analysis software was used to infer the molecular pathways controlled by these genes [19,20].

2.4. Dual-Luciferase Reporter Assay

We conducted a dual-luciferase reporter assay to confirm that miR-139-3p binds directly to the 3′UTR of the TRIP13 gene. The miR-139-3p-binding sequence cloned into the psiCHEK2 vector (C8021; Promega, Madison, WI, USA) is shown in Figure S1. The dual-luciferase reporter assay procedure was described in our previous studies [10,15,16,17]. Transfection of the purified plasmid vectors into LUAD cells was performed using Lipofectamine 2000 (Invitrogen) at 50 ng/well. After 72 h of transfection, we conducted dual-luciferase reporter assays using the Dual Luciferase Reporter Assay System (Promega). The reagents used for these analyses are listed in Table S1.

2.5. Anticancer Effects of DCZ0415 in LUAD Cells

To determine the anticancer effects of cisplatin, carboplatin, and DCZ0415 (TRIP13 inhibitor) in LUAD cells, XTT assays were conducted. XTT assays using Cell Proliferation Kits (catalog no.: 20-300-1000, Biological Industries, Beit-Haemek, Israel) were performed to assess cell proliferation. LUAD cells were plated at 4.0 × 104 cells per well in 96-well plates. Then, the half-maximal inhibitory concentration (IC50) values of these drugs were calculated using GraphPad Prism8 software.

2.6. Western Blotting and Immunohistochemistry

The procedures for Western blotting and immunohistochemistry have been described in our previous studies [10,15,16,17]. LUAD cell lysates were prepared using a RIPA Lysis Buffer System (catalog no.: sc-24948, Santa Cruz Biotechnology Inc., Dallas, TX, USA). Protein concentrations were measured using a PierceTM BCA Protein Assay Kit (catalog no.: 23227, Thermo Fisher Scientific, Rockford, IL, USA). Protein (20 μg) was injected into each well of SuperSepTM Ace (7.5%, 13 well) (catalog no.: 198-14941, FUJIFILM Wako Pure Chemical Corporation, Osaka, Japan) and electrophoresis was performed. Precision Plus ProteinTM Dual Color Standards (catalog no.: #1610374, Bio-Rad Laboratories, Inc., Hercules, CA, USA) were used as the standard. Proteins were transferred to polyvinylidene fluoride membranes (catalog no.: PPVH00010, Merck KGaA, Darmstadt, Germany). The membranes were incubated with blocking buffer (5% skimmed milk) (catalog no.: 190-12865, FUJIFILM Wako Cells 2023, 12, 1885 4 of 20 Pure Chemical Corporation, Osaka, Japan) in TBST. The signal was detected using Amersham ECL Prime Western Blotting Detection Reagent (Cytiva, Marlborough, MA, USA). The antibodies used in the present study are listed in Table S1.

2.7. Statistical Analysis

In this study, statistical analyses were performed using JMP Pro 16 (SAS Institute Inc., Cary, NC, USA). Analysis of differences between the two groups used Welch’s t-test. Differences among multiple groups were analyzed using Dunnett’s test. Patient survival rates were analyzed using Kaplan–Meier survival curves and the log–rank test.

3. Results

3.1. Tumor-Suppressive Function of miR-139-3p in LUAD Cells

Analysis of our miRNA expression signature created by RNA sequencing showed that both strands of pre-miR-139 were downregulated in LUAD tissues (Figure 1A). According to the miRNA database miRBase, miR-139-5p is annotated as the guide strand derived from pre-miR-139 and miR-139-3p as the passenger strand (Figure 1B). Downregulation of these miRNAs in LUAD clinical specimens was confirmed by analysis of TCGA-LUAD datasets from TCGA (Figure 1C). A positive correlation was detected between miR-139-5p and miR-139-3p expression levels by Spearman’s rank analysis (r = 0.513, p < 0.001; Figure 1D).
Our research interests were to elucidate the function of the passenger strands of miRNAs and identify novel oncogenic pathways regulated by tumor-suppressive miRNAs in LUAD cells. The tumor-suppressive roles of miR-139-3p were evaluated by transient transfection assays using two LUAD cell lines, A549 and H1299. We also investigated the functional significance of miR-139-5p (the guide strand) in LUAD cells. Surprisingly, compared to miR-139-3p, miR-139-5p had poor tumor suppressive function in LUAD cells (Figure S2).
Cancer cell proliferation was significantly suppressed by ectopic expression of miR-139-3p in A549 and H1299 cells (Figure 2A). Cell cycle assays demonstrated an increased proportion of cells in the G0/G1 phase after the induction of miR-139-3p expression (Figure 2B). Furthermore, cell invasion and migration were significantly suppressed by ectopic expression of miR-139-3p in LUAD cells (Figure 2C,D). Typical images from the invasion and migration assays after miR-139-3p transfection are shown in Figure S3.
These results suggest that miR-139-3p acts as a tumor-suppressive miRNA in LUAD cells.

3.2. Identification of Cancer-Promoting Genes Regulated by miR-139-3p in LUAD Cells

Next, we aimed to clarify the molecular network regulated by tumor-suppressive miR-139-3p in LUAD cells. Our strategy for identifying miR-139-3p targets is shown in Figure 3.
In a search of the TargetScanHuman database (release 8.0), 3145 genes were identified as having miR-139-3p-binding sequences within their 3′UTR. We applied microarray analysis to search for genes with suppressed expression after miR-139-3p transfection in A549 cells. A total of 7608 genes were downregulated by miR-139-3p transfection. Gene expression data were deposited in the GEO database (accession number GSE242241). We integrated these two datasets to narrow down the gene list, and 1108 genes remained as candidate miR-139-3p targets in LUAD cells.
Among these target genes, 52 were upregulated in LUAD tissues according to the analysis of the TCGA-LUAD dataset (Table 1), of which 21 genes (Bold) were significantly upregulated (Figure 4). Notably, high expression of these 21 genes was negatively associated with the prognosis of LUAD patients (Figure 5).
Moreover, we investigated the extent to which expression of a negative correlation between miR-139-3p and their target genes in LUAD clinical specimens (Figure S4). Spearman’s rank test indicated negative correlations between the expression levels of miR-139-3p and 13 target genes (CENPM, SPC24, ORC1, TRIP13, ARHGEF39, MKI67, KIF18B, CHAF1B, CHEK1, HELLS, FAM111B, MELK, and CENPF).

3.3. Direct Regulation of TRIP13 by miR-139-3p in LUAD Cells

To identify new therapeutic targets for LUAD from among the 52 genes (Table 1), we functionally classified these genes using GeneCodis4 software (https://genecodis.genyo.es/ (accessed on 24 May 2023)) [17]. Among these 52 genes, 13 are involved in “cell cycle” and “DNA replication” (Table 2), and 2 (TRIP13 and CENPF) are involved in cell division checkpoints. Some genes (e.g., MCM2, CHEK1, GINS2, ADAM8, and THY1) were involved in multiple molecular pathways. Analysis of these genes is a critical issue for understanding the malignant transformation of LUAD. However, genes that did not affect prognosis were excluded from further study. Since an inhibitor is available for TRIP13, we conducted further analyses using the TRIP13 inhibitor DCZ0415.
First, we confirmed that TRIP13 mRNA and protein levels were suppressed by ectopic expression of miR-139-3p in LUAD cells (Figure 6A,B). Full-size images of the Western blots are shown in Figure S5.
Next, we confirmed by luciferase reporter assay that miR-139-3p binds directly to the 3’UTR of the TRIP13 gene. The putative miR-139-3p-binding sequence within TRIP13 is shown in Figure 6C. Luciferase activity was significantly suppressed when miR-139-3p and a vector containing a miR-139-3p-binding sequence were simultaneously transfected into LUAD cells (Figure 6D). In contrast, when a vector lacking the miR-139-3p-binding sequence was used, no reduction in luciferase activity was observed (Figure 6D).
These results revealed that miR-139-3p binds directly to the predicted binding sites in TRIP13 genes and regulates its expression in LUAD cells.

3.4. Expression and Clinical Significance of TRIP13 in LUAD

Immunostaining was performed to verify the localization of TRIP13 expression in LUAD clinical specimens. Stronger immunostaining of the TRIP13 protein was observed in cancer than in normal lung tissues (Figure 7A).
Multivariate analysis identified TRIP13 expression as an independent prognostic factor for LUAD after adjusting for the clinical prognostic factors of Stage, T-factor, N-factor, age, and gender (Figure 7B). Specifically, high TRIP13 expression was associated with a poorer 5-year overall survival rate.
Furthermore, we performed gene set enrichment analysis using TCGA-LUAD data to investigate which molecular pathways were dysregulated in LUAD patients with high TRIP13 expression. The “cell cycle”, “DNA repair”, “proteasome”, “P53 signaling pathway”, homologous recombination”, and “mismatch repair” pathways were found to be enriched in patients with high TRIP13 expression (Table 3, Figure 7C). It has been revealed that high expression of TRIP13 affected the prognosis of LUAD patients. An important question is which molecular pathway is affected by the high expression of the TRIP13 gene in lung cancer cells.

3.5. Effects of TRIP13 Knockdown by siRNAs and TRIP13 Inhibition by a Specific Inhibitor in LUAD Cells

To investigate the functional significance of TRIP13 in LUAD cells, we performed analysis using siRNA to knock down TRIP13. The two siRNAs used in this study (siTRIP13-1 and siTRIP13-2) were confirmed to significantly suppress TRIP13 expression (both at the mRNA and protein levels) in LUAD cells (Figure 8A,B). Full-size images of the Western blots are shown in Figure S6.
Cell proliferation was significantly suppressed after transfection of both siRNAs in the LUAD cells (Figure 8C). Cell cycle assays demonstrated an increased proportion of cells in the G0/G1 phase after transfection of both siRNAs in LUAD cells (Figure 8D).

3.6. Combined Effects of a TRIP13 Inhibitor (DCZ0415) and Anticancer Drugs on LUAD Cells

We also investigated the anticancer effects of a TRIP13 inhibitor (DCZ0415) combined with anticancer drugs (cisplatin and carboplatin) on LUAD cells. We found that DCZ0415 inhibited the proliferation of LUAD cells in a concentration-dependent manner (Figure 9).
Furthermore, the sensitivity of LUAD cells to anticancer drugs (cisplatin and carboplatin) was increased by DCZ0415. Co-treatment of DCZ0415 together with cisplatin or carboplatin resulted in increased sensitivity of A549 and H1299 cells to cisplatin or carboplatin (Figure 10). The IC50 of cisplatin significantly decreased from 3.55 µM to 2.17 µM in A549 cells and from 2.38 µM to 2.03 µM in H1299 cells. The IC50 of carboplatin significantly decreased from 40.98 µM to 21.22 µM in A549 cells and from 59.42 µM to 28.45 µM in H1299 cells.
Additionally, the Chou–Talalay method was used to determine the synergistic effect of two anticancer drugs (cisplatin and carboplatin) and DCZ0415. The results revealed that simultaneous administration of anticancer drugs and DCZ0415 had synergistic effects (Figure 11).

4. Discussion

In miRNA biogenesis, miRNAs, which function as single-stranded RNAs, are derived from miRNA precursors. Of the two single strands derived from a miRNA precursor, the functional strand is defined as the guide strand and the non-functional strand as the passenger strand [21,22]. Generally, the miRNAs evaluated in cancer research have been guide strands.
The precursor-miR-139 we focused on in this study is located within the second intron of the phosphodiesterase 2A (PDE2A) gene on chromosome 11q13.4 [23]. The guide strand derived from precursor-miR-139 is miR-139-5p, on which many studies have been conducted [23,24]. In cancer research, the downregulation of miR-139-5p has been reported in most cancer types (e.g., breast, gastric, colorectal, bladder, and head and neck cancers, among others) [25,26,27,28,29]. In contrast, a high expression of miR-139-5p has been reported in cancer tissues compared with normal tissues in gastrointestinal stromal tumors [23,30]. In lung cancer, several studies showed that miR-139-5p was downregulated in LUAD tissues and that it functions as a tumor suppressor by targeting several cancer-related genes [23,31,32,33,34].
From previous reports, it is clear that aberrant expression of miR-139-5p plays important roles in various types of cancer cells. In contrast, there have not been many studies on the role of miR-139-3p, the passenger strand, in cancer cells.
We have generated miRNA expression signatures for various cancer types using RNA sequencing and selected tumor-suppressive miRNA candidates based on those signatures [10,35,36,37]. Of these candidates, we identified downregulation and tumor-suppressive functions of miR-139-3p in various cancers, including bladder cancer, renal cell carcinoma, head and neck squamous cell carcinoma, and colorectal cancer [11,12,13,14]. Here, we showed that miR-139-3p functions as a tumor-suppressive miRNA in LUAD in addition to other cancer types.
Previous studies showed that miR-139-3p is downregulated in NSCLC tissues and that its expression suppresses cancer cell aggressiveness in vitro and in vivo [38,39]. Furthermore, ELAVL1 is directly regulated by miR-139-3p in NSCLC cells [38]. ELAVL1 is a member of the ELAVL family of RNA-binding proteins, which contain several RNA recognition motifs, and ELAVL1 overexpression has been observed in several cancers [40,41,42,43,44]. More recently, downregulation of miR-139-3p was detected in lung squamous cell carcinoma, and expression of miR-139-3p restrained cell growth and accelerated the cell cycle in this cancer [45]. Moreover, the expression of checkpoint kinase 1 (CHEK1) was regulated directly by miR-139-3p in lung squamous cell carcinoma cells [45]. CHEK1 is a serine/threonine protein kinase and a key modulator of DNA damage checkpoints [46,47]. Our present analysis and previous reports revealed that miR-139-3p has tumor-suppressive functions in lung cancer.
Next, we aimed to elucidate the oncogenic network that is controlled by tumor-suppressive miR-139-3p in LUAD cells. In this study, 52 genes were identified as putative targets of miR-139-3p in LUAD cells. Among these target genes, 21 were closely involved in the molecular pathogenesis of LUAD. Functional classification of these genes revealed that 11 (HELLS, MELK, CHEK1, CHAF1B, CENPF, KIF18B, MKI67, SPC24, FAM111B, ORC1, and TRIP13) are involved in the cell cycle, DNA replication, mitotic checkpoint, and cell division. A very interesting finding is that miR139-3p, a passenger strand, controls genes essential for cell division and maintenance.
The AAA (ATPases associated with various cellular activities) protein family is widely conserved from bacteria to humans, and these ATPases use the energy from ATP hydrolysis to unfold proteins and dissociate complexes [48,49]. TRIP13 is a member of the AAA protein family, and it functions as a spindle assembly checkpoint for accurate chromosome segregation during cell division [50]. A recent report showed that TRIP13 plays an important role as a gatekeeper for cell maintenance and management (e.g., DNA break repair, recombination, and chromosome synapsis) [51,52].
Aberrant expression of TRIP13 has been reported in a wide range of cancers, including lung cancer [50,53,54]. The knockdown of TRIP13 in cancer cells significantly suppressed the malignant phenotypes of the cancer cells (i.e., cell proliferation, invasion, and drug resistance) [55,56]. In bladder cancer, overexpression of TRIP13 enhanced the resistance of cancer cells to cisplatin and doxorubicin [56]. In head and neck cancer, epidermal growth factor receptor (EGFR)-mediated phosphorylation of TRIP13 (pY56) enhanced non-homologous end-joining repair and induced radiation resistance in cancer cells [57]. More recently, overexpression of TRIP13 promoted gefitinib (EGFR inhibitor) resistance in NSCLC cells by controlling autophagy and activating EGFR-mediated oncogenic signaling pathways [58]. Since TRIP13 overexpression is involved in drug resistance in various cancer cells, TRIP13 might be a potential target molecule for cancer therapy. TRIP13 overexpression is common in BRCA1-deficient cancers, confers PARP inhibitor resistance, and correlates with poor prognosis [59].
Recent studies have reported that various noncoding RNA molecules are involved in the aberrant expression of TRIP13 in cancer cells. miR-515-5p was found to be downregulated in prostate cancer tissues, and TRIP13 was regulated directly by miR-515-5p [60]. In hepatocellular carcinoma (HCC), upregulated TRIP13 induced malignant transformation of HCC in vitro and lung metastasis in vivo [61,62]. In addition, miR-192-5p regulated TRIP13 expression in HCC cells [63]. In colorectal cancer, aberrant expression of TRIP13 contributed markedly to the aggressive phenotype of cancer cells, and its overexpression induced the downregulation of miR-129-5p and miR-4693-5p [64,65]. More recently, miR-30c-1-3p and miR-30c-2-3p were found to be significantly downregulated in breast cancer (BrCa) tissues, and their target gene TRIP13 was overexpressed in BrCa tissues [55]. Moreover, overexpression of TRIP13 facilitated the aggressiveness of BrCa cells [55]. Based on this study and previous data, TRIP13 overexpression is partially due to the downregulation of several tumor-suppressive miRNAs in cancer cells.

5. Conclusions

Based on the miRNA expression signature of LUAD, we focused on miR-139-3p, a passenger strand, and clarified its tumor-suppressive function in LUAD cells. We revealed that miR-139-3p target molecules are intricately involved in the molecular pathogenesis of LUAD. Inhibition of TRIP13 using a specific inhibitor (DCZ0415) attenuated the malignant transformation of LUAD cells. Interestingly, DCZ0415 exhibited a synergistic effect on the suppression of cell proliferation when used in combination with anticancer drugs. TRIP13 is a potential therapeutic target for LUAD.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15235571/s1, Figure S1: The vector insertion sequence by luciferase-reporter assay; Figure S2: Effects of ectopic expression of miR-139-5p and miR-139-3p in LUAD cells (A549 and H1299). *, p < 0.05; ****, p < 0.0001; N.S., not significant; Figure S3: Atipical images of cell migration or cell invasion assays; Figure S4: Expression correlation between miR-139-3p and their target genes in LUAD clinical specimens; Figure S5: full-sized images of Western blot analysis (TRIP13 antibody) following ectopic expressions of miR-139-3p in A549 and H1299 cells; Figure S6: full-sized images of Western blotting (TRIP13 antibody) following the transfection of siRNAs (siTRIP13-1, siTRIP13-2) in A549 and H1299 cells; Table S1: reagents used in this study; Table S2: information of tissues by immunostaining.

Author Contributions

Conceptualization, T.S., K.M. and N.S.; methodology, N.S.; validation, H.I., K.M. and N.S.; formal analysis, Y.H., Y.T., T.S., M.S., Y.G., N.K. and M.K.; investigation, Y.H., Y.T., T.S., M.S., Y.G., N.K. and M.K.; resources, Y.G.; data curation, Y.H., Y.T. and T.S.; writing—original draft preparation, N.S.; writing—review and editing, Y.H., Y.T., T.S. and N.S.; visualization, Y.H., Y.T., T.S. and N.S.; supervision, N.S.; project administration, H.I., K.M. and N.S.; funding acquisition, S.M., N.K., M.K., K.M. and N.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by KAKENHI; grant numbers 22K16196, 22K09679, 21K09367, 22K08260, and 21K09577.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee on Epidemiological and its related Studies, Sakuragaoka Campus, Kagoshima University (approval no. 210101 eki-kai 2, 31 August 2023).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Publicly available datasets were analyzed in this study. These data can be accessed here: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE230229 (accessed on 19 September 2023) and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE242241 (accessed on 19 September 2023).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLO-BOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
  2. Siegel, R.L.; Miller, K.D.; Fuchs, H.E.; Jemal, A. Cancer statistics, 2022. CA Cancer J. Clin. 2022, 72, 7–33. [Google Scholar] [CrossRef] [PubMed]
  3. Thai, A.A.; Solomon, B.J.; Sequist, L.V.; Gainor, J.F.; Heist, R.S. Lung cancer. Lancet 2021, 398, 535–554. [Google Scholar] [PubMed]
  4. Alexander, M.; Kim, S.Y.; Cheng, H. Update 2020: Management of Non-Small Cell Lung Cancer. Lung 2020, 198, 897–907. [Google Scholar] [CrossRef] [PubMed]
  5. Fatica, A.; Bozzoni, I. Long non-coding RNAs: New players in cell differentiation and development. Nat. Rev. Genet. 2014, 15, 7–21. [Google Scholar] [CrossRef] [PubMed]
  6. Bartel, D.P. MicroRNAs: Target recognition and regulatory functions. Cell 2009, 136, 215–233. [Google Scholar] [CrossRef] [PubMed]
  7. Anastasiadou, E.; Jacob, L.S.; Slack, F.J. Non-coding RNA networks in cancer. Nat. Rev. Cancer 2018, 18, 5–18. [Google Scholar] [CrossRef] [PubMed]
  8. Wu, K.L.; Tsai, Y.M.; Lien, C.T.; Kuo, P.L.; Hung, A.J. The Roles of MicroRNA in Lung Cancer. Int. J. Mol. Sci. 2019, 20, 1611. [Google Scholar] [CrossRef]
  9. Matranga, C.; Tomari, Y.; Shin, C.; Bartel, D.P.; Zamore, P.D. Passenger-strand cleavage facilitates assembly of siRNA into Ago2-containing RNAi enzyme complexes. Cell 2005, 123, 607–620. [Google Scholar] [CrossRef]
  10. Tomioka, Y.; Suetsugu, T.; Seki, N.; Tanigawa, K.; Hagihara, Y.; Shinmura, M.; Asai, S.; Kikkawa, N.; Inoue, H.; Mizuno, K. The Molecular Pathogenesis of Tumor-Suppressive miR-486-5p and miR-486-3p Target Genes: GINS4 Facilitates Aggressiveness in Lung Adenocarcinoma. Cells 2023, 12, 1885. [Google Scholar] [CrossRef]
  11. Yonemori, M.; Seki, N.; Yoshino, H.; Matsushita, R.; Miyamoto, K.; Nakagawa, M.; Enokida, H. Dual tumor-suppressors miR-139-5p and miR-139-3p targeting matrix metalloprotease 11 in bladder cancer. Cancer Sci. 2016, 107, 1233–1242. [Google Scholar] [PubMed]
  12. Okada, R.; Goto, Y.; Yamada, Y.; Kato, M.; Asai, S.; Moriya, S.; Ichikawa, T.; Seki, N. Regulation of Oncogenic Targets by the Tumor-Suppressive miR-139 Duplex (miR-139-5p and miR-139-3p) in Renal Cell Carcinoma. Biomedicines 2020, 8, 599. [Google Scholar]
  13. Koma, A.; Asai, S.; Minemura, C.; Oshima, S.; Kinoshita, T.; Kikkawa, N.; Koshizuka, K.; Moriya, S.; Kasamatsu, A.; Hanazawa, T.; et al. Impact of Oncogenic Targets by Tumor-Suppressive miR-139-5p and miR-139-3p Regulation in Head and Neck Squamous Cell Carcinoma. Int. J. Mol. Sci. 2021, 22, 9947. [Google Scholar] [PubMed]
  14. Yasudome, R.; Seki, N.; Asai, S.; Goto, Y.; Kita, Y.; Hozaka, Y.; Wada, M.; Tanabe, K.; Idichi, T.; Mori, S.; et al. Molecular Pathogenesis of Colorectal Cancer: Impact of Oncogenic Targets Regulated by Tumor Suppressive miR-139-3p. Int. J. Mol. Sci. 2022, 23, 11616. [Google Scholar]
  15. Asai, S.; Koma, A.; Nohata, N.; Kinoshita, T.; Kikkawa, N.; Kato, M.; Minemura, C.; Uzawa, K.; Hanazawa, T.; Seki, N. Impact of miR-1/miR-133 Clustered miRNAs: PFN2 Facilitates Malignant Phenotypes in Head and Neck Squamous Cell Carcinoma. Biomedicines 2022, 10, 663. [Google Scholar]
  16. Suetsugu, T.; Koshizuka, K.; Seki, N.; Mizuno, K.; Okato, A.; Arai, T.; Misono, S.; Uchida, A.; Kumamoto, T.; Inoue, H. Down-regulation of matrix metalloproteinase 14 by the antitumor miRNA, miR-150-5p, inhibits the aggressiveness of lung squamous cellcarcinoma cells. Int. J. Oncol. 2018, 52, 913–924. [Google Scholar] [PubMed]
  17. Tanigawa, K.; Tomioka, Y.; Misono, S.; Asai, S.; Kikkawa, N.; Hagihara, Y.; Suetsugu, T.; Inoue, H.; Mizuno, K.; Seki, N. Mini-chromosome maintenance proteins in lung adenocarcinoma: Clinical significance and therapeutic targets. FEBS Open Bio. 2023, 13, 1737–1755. [Google Scholar] [PubMed]
  18. Garcia-Moreno, A.; López-Domínguez, R.; Villatoro-García, J.A.; Ramirez-Mena, A.; Aparicio-Puerta, E.; Hackenberg, M.; Pas-cual-Montano, A.; Carmona-Saez, P. Functional Enrichment Analysis of Regulatory Elements. Biomedicines 2022, 10, 590. [Google Scholar]
  19. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar]
  20. Mootha, V.K.; Lindgren, C.M.; Eriksson, K.F.; Subramanian, A.; Sihag, S.; Lehar, J.; Puigserver, P.; Carlsson, E.; Ridderstråle, M.; Laurila, E.; et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in humandiabetes. Nat. Genet. 2003, 34, 267–273. [Google Scholar]
  21. Mitra, R.; Sun, J.; Zhao, Z. microRNA regulation in cancer: One arm or two arms? Int. J. Cancer 2015, 137, 1516–1518. [Google Scholar] [PubMed]
  22. Mitra, R.; Adams, C.M.; Jiang, W.; Greenawalt, E.; Eischen, C.M. Pan-cancer analysis reveals cooperativity of both strands of microRNA that regulate tumorigenesis and patient survival. Nat. Commun. 2020, 11, 968. [Google Scholar] [CrossRef] [PubMed]
  23. Khalili, N.; Nouri-Vaskeh, M.; Hasanpour Segherlou, Z.; Baghbanzadeh, A.; Halimi, M.; Rezaee, H.; Baradaran, B. Diagnostic, prognostic, and therapeutic significance of miR-139-5p in cancers. Life Sci. 2020, 256, 117865. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, Y.H.; Ji, J.; Weng, H.; Wang, B.C.; Wang, F.B. MiR-139 in digestive system tumor diagnosis and detection: Bioinformatics and meta-analysis. Clin. Chim. Acta 2018, 485, 33–41. [Google Scholar] [CrossRef] [PubMed]
  25. Sun, H.; Dai, J.; Chen, M.; Chen, Q.; Xie, Q.; Zhang, W.; Li, G.; Yan, M. miR-139-5p Was Identified as Biomarker of Different Molecular Subtypes of Breast Carcinoma. Front. Oncol. 2022, 12, 857714. [Google Scholar] [CrossRef]
  26. Li, Y.; Sun, Y.; Li, Z.; Li, S.; Wu, C. MiR-139-5p Inhibits the Development of Gastric Cancer through Targeting TPD52. J. Healthc. Eng. 2022, 2022, 4033373. [Google Scholar] [CrossRef]
  27. Xu, K.; Shen, K.; Liang, X.; Li, Y.; Nagao, N.; Li, J.; Liu, J.; Yin, P. MiR-139-5p reverses CD44+/CD133+-associated multidrug resistance by downregulating NOTCH1 in colorectal carcinoma cells. Oncotarget 2016, 7, 75118–75129. [Google Scholar] [CrossRef]
  28. Luo, H.; Yang, R.; Li, C.; Tong, Y.; Fan, L.; Liu, X.; Xu, C. MicroRNA-139-5p inhibits bladder cancer proliferation and self-renewal by targeting the Bmi1 oncogene. Tumor Biol. 2017, 39, 1010428317718414. [Google Scholar] [CrossRef]
  29. Ferreira, L.A.M.; Bezerra, M.; Kawasaki-Oyama, R.S.; Fernandes, G.M.M.; Castanhole-Nunes, M.M.U.; Serafim Junior, V.; Cas-tilho, R.M.; Pavarino, É.C.; Maniglia, J.V.; Goloni-Bertollo, E.M. Effect of ZEB1 Associated with microRNAs on Tumor Stem Cells in Head and Neck Cancer. Int. J. Mol. Sci. 2023, 24, 5916. [Google Scholar] [CrossRef]
  30. Pantaleo, M.A.; Ravegnini, G.; Astolfi, A.; Simeon, V.; Nannini, M.; Saponara, M.; Urbini, M.; Gatto, L.; Indio, V.; Sammarini, G.; et al. Integrating miRNA and gene expression profiling analysis revealed regulatory networks in gastrointestinal stromal tumors. Epigenomics 2016, 8, 1347–1366. [Google Scholar] [CrossRef]
  31. Xu, W.; Hang, M.; Yuan, C.Y.; Wu, F.L.; Chen, S.B.; Xue, K. MicroRNA-139-5p inhibits cell proliferation and invasion by tar-geting insulin-like growth factor 1 receptor in human non-small cell lung cancer. Int. J. Clin. Exp. Pathol. 2015, 8, 3864–3870. [Google Scholar] [PubMed]
  32. Shao, Y.; Liang, B.; Long, F.; Jiang, S.J. Diagnostic MicroRNA Biomarker Discovery for Non-Small-Cell Lung Cancer Adeno-carcinoma by Integrative Bioinformatics Analysis. Biomed. Res. Int. 2017, 2017, 2563085. [Google Scholar] [CrossRef] [PubMed]
  33. Huang, N.; Guo, W.; Ren, K.; Li, W.; Jiang, Y.; Sun, J.; Dai, W.; Zhao, W. LncRNA AFAP1-AS1 Supresses miR-139-5p and Promotes Cell Proliferation and Chemotherapy Resistance of Non-small Cell Lung Cancer by Competitively Upregulating RRM2. Front. Oncol. 2019, 9, 1103. [Google Scholar] [CrossRef] [PubMed]
  34. Zhang, Z.; Li, W.; Jiang, D.; Liu, C.; Lai, Z. MicroRNA-139-5p inhibits cell viability, migration and invasion and suppresses tumor growth by targeting HDGF in non-small cell lung cancer. Oncol. Lett. 2020, 19, 1806–1814. [Google Scholar] [CrossRef] [PubMed]
  35. Kurozumi, A.; Goto, Y.; Okato, A.; Ichikawa, T.; Seki, N. Aberrantly expressed microRNAs in bladder cancer and renal cell carcinoma. J. Hum. Genet. 2017, 62, 49–56. [Google Scholar] [CrossRef] [PubMed]
  36. Koshizuka, K.; Hanazawa, T.; Fukumoto, I.; Kikkawa, N.; Okamoto, Y.; Seki, N. The microRNA signatures: Aberrantly ex-pressed microRNAs in head and neck squamous cell carcinoma. J. Hum. Genet. 2017, 62, 3–13. [Google Scholar] [CrossRef] [PubMed]
  37. Nepal, P.; Hozaka, Y.; Tanaka, T.; Wada, M.; Asai, S.; Minemura, C.; Idichi, T.; Arigami, T.; Kurahara, H.; Seki, N.; et al. Impact of Oncogenic Targets Controlled by Tumor-Suppressive miR-30a-5p in Pancreatic Ductal Adenocarcinoma. Anticancer Res. 2021, 41, 4821–4836. [Google Scholar] [CrossRef] [PubMed]
  38. Ni, Z.Z.; He, J.K.; Tang, X.; Tao, Z.; Zhang, Y.; Xie, B. Identification of ELAVL1 gene and miRNA-139-3p involved in the ag-gressiveness of NSCLC. Eur. Rev. Med. Pharmacol. Sci. 2020, 24, 9453–9464. [Google Scholar]
  39. Ma, Y.; Gong, Z.; Wang, H.; Liang, Y.; Huang, X.; Yu, G. Anti-cancer effect of miR-139-3p on laryngeal squamous cell carcinoma by targeting rab5a: In vitro and in vivo studies. Pathol. Res. Pract. 2020, 216, 153194. [Google Scholar] [CrossRef]
  40. Lebedeva, S.; Jens, M.; Theil, K.; Schwanhäusser, B.; Selbach, M.; Landthaler, M.; Rajewsky, N. Transcriptome-wide analysis of regulatory interactions of the RNA-binding protein HuR. Mol. Cell 2011, 43, 340–352. [Google Scholar] [CrossRef]
  41. Wu, X.; Xu, L. The RNA-binding protein HuR in human cancer: A friend or foe? Adv. Drug Deliv. Rev. 2022, 184, 114179. [Google Scholar] [CrossRef] [PubMed]
  42. Gu, C.; Zhang, M.; Sun, W.; Dong, C. Upregulation of miR-324-5p Inhibits Proliferation and Invasion of Colorectal Cancer Cells by Targeting ELAVL1. Oncol. Res. 2019, 27, 515–524. [Google Scholar] [CrossRef] [PubMed]
  43. Zhang, F.; Cai, Z.; Lv, H.; Li, W.; Liang, M.; Wei, X.; Zhou, F. Multiple functions of HuR in urinary tumors. J. Cancer Res. Clin. Oncol. 2019, 145, 11–18. [Google Scholar] [CrossRef] [PubMed]
  44. Chou, S.D.; Murshid, A.; Eguchi, T.; Gong, J.; Calderwood, S.K. HSF1 regulation of β-catenin in mammary cancer cells through control of HuR/elavL1 expression. Oncogene 2015, 34, 2178–2188. [Google Scholar] [CrossRef] [PubMed]
  45. Zheng, X.; Zhang, Y.; Wu, S.; Jiang, B.; Liu, Y. MiR-139-3p Targets CHEK1 Modulating DNA Repair and Cell Viability in Lung Squamous Carcinoma Cells. Mol. Biotechnol. 2022, 64, 832–840. [Google Scholar] [CrossRef] [PubMed]
  46. Sanchez, Y.; Wong, C.; Thoma, R.S.; Richman, R.; Wu, Z.; Piwnica-Worms, H.; Elledge, S.J. Conservation of the Chk1 check-point pathway in mammals: Linkage of DNA damage to Cdk regulation through Cdc25. Science 1997, 277, 1497–1501. [Google Scholar] [CrossRef]
  47. Walworth, N.; Davey, S.; Beach, D. Fission yeast chk1 protein kinase links the rad checkpoint pathway to cdc2. Nature 1993, 363, 368–371. [Google Scholar] [CrossRef] [PubMed]
  48. Dimroth, P. Primary sodium ion translocating enzymes. Biochim. Biophys. Acta 1997, 1318, 11–51. [Google Scholar] [CrossRef]
  49. Cerullo, F.; Filbeck, S.; Patil, P.R.; Hung, H.C.; Xu, H.; Vornberger, J.; Hofer, F.W.; Schmitt, J.; Kramer, G.; Bukau, B.; et al. Bacterial ribosome collision sensing by a MutS DNA repair ATPase paralogue. Nature 2022, 603, 509–514. [Google Scholar] [CrossRef]
  50. Lu, S.; Qian, J.; Guo, M.; Gu, C.; Yang, Y. Insights into a Crucial Role of TRIP13 in Human Cancer. Comput. Struct. Biotechnol. J. 2019, 17, 854–861. [Google Scholar] [CrossRef]
  51. Miniowitz-Shemtov, S.; Eytan, E.; Kaisari, S.; Sitry-Shevah, D.; Hershko, A. Mode of interaction of TRIP13 AAA-ATPase with the Mad2-binding protein p31comet and with mitotic checkpoint complexes. Proc. Natl. Acad. Sci. USA 2015, 112, 11536–11540. [Google Scholar] [PubMed]
  52. Vader, G. Pch2(TRIP13): Controlling cell division through regulation of HORMA domains. Chromosoma 2015, 124, 333–339. [Google Scholar] [PubMed]
  53. Li, W.; Zhang, G.; Li, X.; Wang, X.; Li, Q.; Hong, L.; Shen, Y.; Zhao, C.; Gong, X.; Chen, Y.; et al. Thyroid hormone receptor interactor 13 (TRIP13) overexpression associated with tumor progression and poor prognosis in lung adenocarcinoma. Biochem. Biophys Res. Commun. 2018, 499, 416–424. [Google Scholar] [PubMed]
  54. Kang, J.U.; Koo, S.H.; Kwon, K.C.; Park, J.W.; Kim, J.M. Gain at chromosomal region 5p15.33, containing TERT, is the most frequent genetic event in early stages of non-small cell lung cancer. Cancer Genet. Cytogenet. 2008, 182, 1–11. [Google Scholar] [PubMed]
  55. Mitsueda, R.; Toda, H.; Shinden, Y.; Fukuda, K.; Yasudome, R.; Kato, M.; Kikkawa, N.; Ohtsuka, T.; Nakajo, A.; Seki, N. Onco-genic Targets Regulated by Tumor-Suppressive miR-30c-1-3p and miR-30c-2-3p: TRIP13 Facilitates Cancer Cell Aggressiveness in Breast Cancer. Cancers 2023, 15, 4189. [Google Scholar] [PubMed]
  56. Lu, S.; Guo, M.; Fan, Z.; Chen, Y.; Shi, X.; Gu, C.; Yang, Y. Elevated TRIP13 drives cell proliferation and drug resistance inbladder cancer. Am. J. Transl. Res. 2019, 11, 4397–4410. [Google Scholar]
  57. Banerjee, R.; Liu, M.; Bellile, E.; Schmitd, L.B.; Goto, M.; Hutchinson, M.N.D.; Singh, P.; Zhang, S.; Damodaran, D.P.V.; Nyati, M.K.; et al. Phosphorylation of TRIP13 at Y56 induces radiation resistance but sensitizes head and neck cancer to cetuximab. Mol. Ther. 2022, 30, 468–484. [Google Scholar] [PubMed]
  58. Xiao, Z.; Li, M.; Zhang, X.; Rong, X.; Xu, H. TRIP13 overexpression promotes gefitinib resistance in non-small cell lung cancer via regulating autophagy and phosphorylation of the EGFR signaling pathway. Oncol. Rep. 2023, 49, 84. [Google Scholar]
  59. Clairmont, C.S.; Sarangi, P.; Ponnienselvan, K.; Galli, L.D.; Csete, I.; Moreau, L.; Adelmant, G.; Chowdhury, D.; Marto, J.A.; D’Andrea, A.D. TRIP13 regulates DNA repair pathway choice through REV7 conformational change. Nat. Cell Biol. 2020, 22, 87–96. [Google Scholar]
  60. Zhang, X.; Zhou, J.; Xue, D.; Li, Z.; Liu, Y.; Dong, L. MiR-515-5p acts as a tumor suppressor via targeting TRIP13 in prostate cancer. Int. J. Biol. Macromol. 2019, 129, 227–232. [Google Scholar]
  61. Yao, J.; Zhang, X.; Li, J.; Zhao, D.; Gao, B.; Zhou, H.; Gao, S.; Zhang, L. Silencing TRIP13 inhibits cell growth and metastasis of hepatocellular carcinoma by activating of TGF-β1/smad3. Cancer Cell Int. 2018, 18, 208. [Google Scholar] [CrossRef] [PubMed]
  62. Cai, W.; Ni, W.; Jin, Y.; Li, Y. TRIP13 promotes lung cancer cell growth and metastasis through AKT/mTORC1/c-Myc signaling. Cancer Biomark. 2021, 30, 237–248. [Google Scholar] [PubMed]
  63. Zhu, M.X.; Wei, C.Y.; Zhang, P.F.; Gao, D.M.; Chen, J.; Zhao, Y.; Dong, S.S.; Liu, B.B. Elevated TRIP13 drives the AKT/mTOR pathway to induce the progression of hepatocellular carcinoma via interacting with ACTN4. J. Exp. Clin. Cancer Res. 2019, 38, 409. [Google Scholar] [CrossRef] [PubMed]
  64. Cao, Y.; Huang, F.; Liu, J.; Qi, H.; Xiao, J. MiR-129-5p/TRIP13 affects malignant phenotypes of colorectal cancer cells. Histol. Histopathol. 2022, 37, 879–888. [Google Scholar]
  65. Chen, Y.; Chen, D.; Qin, Y.; Qiu, C.; Zhou, Y.; Dai, M.; Li, L.; Sun, Q.; Jiang, Y. TRIP13, identified as a hub gene of tumor progression, is the target of microRNA-4693-5p and a potential therapeutic target for colorectal cancer. Cell Death Discov. 2022, 8, 35. [Google Scholar]
Figure 1. Expression levels of miR-139-5p and miR-139-3p in LUAD clinical tissues. (A) Volcano plot of the miRNA expression signature based on miRNA sequencing (GEO accession number: GSE230229). The log2 fold change (FC) in expression is plotted on the x-axis and the log10 p-value on the y-axis. The blue and red dots represent the downregulated (log2FC < −2.0 and p < 0.05) and upregulated (log2FC > 2.0 and p < 0.05) miRNAs, respectively. (B) Chromosomal location of pre-miR-139 within the human genome. The mature sequences of miR-139-5p (guide strand) and miR-139-3p (passenger strand) are shown. (C) Expression levels of miR-139-5p and miR-139-3p validated in LUAD clinical specimens. The expression of both miRNAs was significantly downregulated in cancer tissues (p < 0.001). (D) Positive correlations (Spearman’s rank test) between miR-139-5p and miR-139-3p expression levels in clinical specimens (r = 0.513, p < 0.001).
Figure 1. Expression levels of miR-139-5p and miR-139-3p in LUAD clinical tissues. (A) Volcano plot of the miRNA expression signature based on miRNA sequencing (GEO accession number: GSE230229). The log2 fold change (FC) in expression is plotted on the x-axis and the log10 p-value on the y-axis. The blue and red dots represent the downregulated (log2FC < −2.0 and p < 0.05) and upregulated (log2FC > 2.0 and p < 0.05) miRNAs, respectively. (B) Chromosomal location of pre-miR-139 within the human genome. The mature sequences of miR-139-5p (guide strand) and miR-139-3p (passenger strand) are shown. (C) Expression levels of miR-139-5p and miR-139-3p validated in LUAD clinical specimens. The expression of both miRNAs was significantly downregulated in cancer tissues (p < 0.001). (D) Positive correlations (Spearman’s rank test) between miR-139-5p and miR-139-3p expression levels in clinical specimens (r = 0.513, p < 0.001).
Cancers 15 05571 g001
Figure 2. Effects of ectopic expression of miR-139-3p in LUAD cells (A549 and H1299). (A) Cell proliferation assessed by XTT assay. At 72 h after transient transfection of miRNAs, cancer cell viability was analyzed. (B) Cell cycle status at 72 h after transfection with miR-139-3p assessed using flow cytometry. (C) Cell invasion assessed using Matrigel invasion assays at 48 h after seeding miR-139-3p-transfected cells into the chambers. (D) Cell migration assessed using a membrane culture system at 48 h after seeding miR-139-3p transfected cells into the chambers. ****, p < 0.0001.
Figure 2. Effects of ectopic expression of miR-139-3p in LUAD cells (A549 and H1299). (A) Cell proliferation assessed by XTT assay. At 72 h after transient transfection of miRNAs, cancer cell viability was analyzed. (B) Cell cycle status at 72 h after transfection with miR-139-3p assessed using flow cytometry. (C) Cell invasion assessed using Matrigel invasion assays at 48 h after seeding miR-139-3p-transfected cells into the chambers. (D) Cell migration assessed using a membrane culture system at 48 h after seeding miR-139-3p transfected cells into the chambers. ****, p < 0.0001.
Cancers 15 05571 g002
Figure 3. Flowchart for identification of miR-139-3p targets in LUAD cell. To identify putative targets of miR-139-3p in LUAD cells, the following two datasets were merged: the TargetScanHuman database (release 8.0) and our original mRNA expression profile (miR-139-3p-transfected A549 cells; GEO accession number: GSE242241). A total of 1108 genes were identified as putative miR-139-3p targets. Furthermore, we searched for genes that were associated with the prognosis of LUAD patients using two databases: GEPIA (http://gepia2.cancer-pku.cn/#analysis (accessed on 17 January 2023)) and OncoLnc (http://www.oncolnc.org (accessed on 17 January 2023)). Of the miR-139-3p target genes, 21 were upregulated in LUAD tissues, and these 21 genes were analyzed further.
Figure 3. Flowchart for identification of miR-139-3p targets in LUAD cell. To identify putative targets of miR-139-3p in LUAD cells, the following two datasets were merged: the TargetScanHuman database (release 8.0) and our original mRNA expression profile (miR-139-3p-transfected A549 cells; GEO accession number: GSE242241). A total of 1108 genes were identified as putative miR-139-3p targets. Furthermore, we searched for genes that were associated with the prognosis of LUAD patients using two databases: GEPIA (http://gepia2.cancer-pku.cn/#analysis (accessed on 17 January 2023)) and OncoLnc (http://www.oncolnc.org (accessed on 17 January 2023)). Of the miR-139-3p target genes, 21 were upregulated in LUAD tissues, and these 21 genes were analyzed further.
Cancers 15 05571 g003
Figure 4. Expression levels of the 21 target genes regulated by miR-139-3p in LUAD. The expression levels of the 21 miR-139-3p target genes (KRT80, CENPM, SPC24, ORC1, MYEOV, TRIP13, GPX8, ARHGEF39, MKI67, KIF18B, CHAF1B, CP, GPRIN1, UCK2, CHEK1, HELLS, CTSV, FAM111B, SLC16A3, MELK, and CENPF) in LUAD clinical specimens were analyzed using the TCGA-LUAD dataset. All genes were upregulated in LUAD tissues (n = 499) compared with normal tissues (n = 58) (p < 0.001).
Figure 4. Expression levels of the 21 target genes regulated by miR-139-3p in LUAD. The expression levels of the 21 miR-139-3p target genes (KRT80, CENPM, SPC24, ORC1, MYEOV, TRIP13, GPX8, ARHGEF39, MKI67, KIF18B, CHAF1B, CP, GPRIN1, UCK2, CHEK1, HELLS, CTSV, FAM111B, SLC16A3, MELK, and CENPF) in LUAD clinical specimens were analyzed using the TCGA-LUAD dataset. All genes were upregulated in LUAD tissues (n = 499) compared with normal tissues (n = 58) (p < 0.001).
Cancers 15 05571 g004
Figure 5. Gene expression in, and 5-year overall survival rate of, patients with LUAD. Kaplan–Meier curves of the 5-year overall survival rates according to the expression of the 21 target genes (KRT80, CENPM, SPC24, ORC1, MYEOV, TRIP13, GPX8, ARHGEF39, MKI67, KIF18B, CHAF1B, CP, GPRIN1, UCK2, CHEK1, HELLS, CTSV, FAM111B, SLC16A3, MELK, and CENPF) are shown. Low expression of all 21 genes was significantly predictive of poorer overall survival in patients with LUAD. The patients (n = 487) were divided into high- and low-expression groups according to the median gene expression level. The red and blue lines represent the high and low expression groups, respectively.
Figure 5. Gene expression in, and 5-year overall survival rate of, patients with LUAD. Kaplan–Meier curves of the 5-year overall survival rates according to the expression of the 21 target genes (KRT80, CENPM, SPC24, ORC1, MYEOV, TRIP13, GPX8, ARHGEF39, MKI67, KIF18B, CHAF1B, CP, GPRIN1, UCK2, CHEK1, HELLS, CTSV, FAM111B, SLC16A3, MELK, and CENPF) are shown. Low expression of all 21 genes was significantly predictive of poorer overall survival in patients with LUAD. The patients (n = 487) were divided into high- and low-expression groups according to the median gene expression level. The red and blue lines represent the high and low expression groups, respectively.
Cancers 15 05571 g005
Figure 6. Direct regulation of TRIP13 by miR-139-3p expression in LUAD cells. (A) Significant reduction in the TRIP13 mRNA level by ectopic expression of miR-139-3p in LUAD cells (A549 and H1299). Total RNA was extracted 72 h after miR-139-3p transfection into LUAD cells, and expression levels were analyzed by real-time PCR. For miRNA expression, GAPDH was used as an internal control. (B) Significant reduction in the TRIP13 protein level by ectopic expression of miR-139-3p in LUAD cells (A549 and H1299). Protein level expression was determined by Western blotting. Proteins were collected 72 h after miR-139-3p transfection. GAPDH was used as an internal control. (C) Putative miR-139-3p binding sites in the 3′UTR of the TRIP13 gene detected in the TargetScanHuman database (release 8.0). (D) Direct binding of miR-139-3p to target sequences was analyzed by dual luciferase reporter assays. These data showed that miR-139-3p bound directly to the target sequence. ***, p < 0.001; ****, p < 0.0001; N.S., not significant.
Figure 6. Direct regulation of TRIP13 by miR-139-3p expression in LUAD cells. (A) Significant reduction in the TRIP13 mRNA level by ectopic expression of miR-139-3p in LUAD cells (A549 and H1299). Total RNA was extracted 72 h after miR-139-3p transfection into LUAD cells, and expression levels were analyzed by real-time PCR. For miRNA expression, GAPDH was used as an internal control. (B) Significant reduction in the TRIP13 protein level by ectopic expression of miR-139-3p in LUAD cells (A549 and H1299). Protein level expression was determined by Western blotting. Proteins were collected 72 h after miR-139-3p transfection. GAPDH was used as an internal control. (C) Putative miR-139-3p binding sites in the 3′UTR of the TRIP13 gene detected in the TargetScanHuman database (release 8.0). (D) Direct binding of miR-139-3p to target sequences was analyzed by dual luciferase reporter assays. These data showed that miR-139-3p bound directly to the target sequence. ***, p < 0.001; ****, p < 0.0001; N.S., not significant.
Cancers 15 05571 g006
Figure 7. Clinical significance of TRIP13 expression in LUAD. (A) Immunohistochemical staining of TRIP13. Immunostaining showed that the TRIP13 protein was strongly expressed in cancer lesions and less in non-cancerous tissues. Scale bar: 200 µm (low magnification); 50 µm (high magnification). (B) Forest plot showing the results of multivariate Cox proportional hazards regression analysis of the 5-year overall survival rate. Patients with high TRIP13 expression had a significantly lower overall survival rate. These data were obtained from TCGA-LUAD datasets. (C) Gene set enrichment analysis (GSEA) was applied to explore molecular pathways mediated by TRIP13 in LUAD cells. The top six pathways enriched in LUAD patients with high TRIP13 expression were cell cycle, DNA replication, proteasome, P53 signaling pathway, homologous recombination, and mismatch repair.
Figure 7. Clinical significance of TRIP13 expression in LUAD. (A) Immunohistochemical staining of TRIP13. Immunostaining showed that the TRIP13 protein was strongly expressed in cancer lesions and less in non-cancerous tissues. Scale bar: 200 µm (low magnification); 50 µm (high magnification). (B) Forest plot showing the results of multivariate Cox proportional hazards regression analysis of the 5-year overall survival rate. Patients with high TRIP13 expression had a significantly lower overall survival rate. These data were obtained from TCGA-LUAD datasets. (C) Gene set enrichment analysis (GSEA) was applied to explore molecular pathways mediated by TRIP13 in LUAD cells. The top six pathways enriched in LUAD patients with high TRIP13 expression were cell cycle, DNA replication, proteasome, P53 signaling pathway, homologous recombination, and mismatch repair.
Cancers 15 05571 g007aCancers 15 05571 g007b
Figure 8. Effects of knockdown of TRIP13 by siRNAs in LUAD cells (A549 and H1299). (A) TRIP13 mRNA levels were effectively blocked by each siRNA in LUAD cells (A549 and H1299). (B) TRIP13 protein levels were effectively inhibited by two siRNAs (siTRIP13-1 and siTRIP13-2) in LUAD cells (A549 and H1299). (C) Cell proliferation was assessed using XTT assays 72 h after siRNA transfection into LUAD cells. Cell proliferation was significantly blocked after transient transfection of siRNAs. (D) Flow cytometry analysis of cell cycle status 72 h after transfection with siTRIP13-1 and siTRIP13-2. ****, p < 0.0001.
Figure 8. Effects of knockdown of TRIP13 by siRNAs in LUAD cells (A549 and H1299). (A) TRIP13 mRNA levels were effectively blocked by each siRNA in LUAD cells (A549 and H1299). (B) TRIP13 protein levels were effectively inhibited by two siRNAs (siTRIP13-1 and siTRIP13-2) in LUAD cells (A549 and H1299). (C) Cell proliferation was assessed using XTT assays 72 h after siRNA transfection into LUAD cells. Cell proliferation was significantly blocked after transient transfection of siRNAs. (D) Flow cytometry analysis of cell cycle status 72 h after transfection with siTRIP13-1 and siTRIP13-2. ****, p < 0.0001.
Cancers 15 05571 g008
Figure 9. Effects of treatment with DCZ0415 (TRIP13 inhibitor) on LUAD cells (A549 and H1299). The proliferation of LUAD cells was significantly inhibited by DCZ0415 in a concentration-dependent manner. ****, p < 0.0001; N.S., not significant.
Figure 9. Effects of treatment with DCZ0415 (TRIP13 inhibitor) on LUAD cells (A549 and H1299). The proliferation of LUAD cells was significantly inhibited by DCZ0415 in a concentration-dependent manner. ****, p < 0.0001; N.S., not significant.
Cancers 15 05571 g009
Figure 10. Effects of co-treatment of DCZ0415 (TRIP13 inhibitor) with anticancer drugs (cisplatin and carboplatin) in LUAD cells (A549 and H1299). LUAD cells showed increased sensitivity to anticancer drugs (cisplatin and carboplatin) when co-treated with DCZ0415. CDDP, cisplatin; CBDCA, carboplatin.
Figure 10. Effects of co-treatment of DCZ0415 (TRIP13 inhibitor) with anticancer drugs (cisplatin and carboplatin) in LUAD cells (A549 and H1299). LUAD cells showed increased sensitivity to anticancer drugs (cisplatin and carboplatin) when co-treated with DCZ0415. CDDP, cisplatin; CBDCA, carboplatin.
Cancers 15 05571 g010
Figure 11. Synergistic effects between two anticancer drugs (cisplatin and carboplatin) and the TRIP13 inhibitor DCZ0415. The Chou–Talalay method was used to determine the synergistic effects between two anticancer drugs (cisplatin and carboplatin) and DCZ0415 in LUAD cells (A549 and H1299). CDDP, cisplatin; CBDCA, carboplatin.
Figure 11. Synergistic effects between two anticancer drugs (cisplatin and carboplatin) and the TRIP13 inhibitor DCZ0415. The Chou–Talalay method was used to determine the synergistic effects between two anticancer drugs (cisplatin and carboplatin) and DCZ0415 in LUAD cells (A549 and H1299). CDDP, cisplatin; CBDCA, carboplatin.
Cancers 15 05571 g011
Table 1. Putative target genes regulated by miR-139-3p in A549 cells.
Table 1. Putative target genes regulated by miR-139-3p in A549 cells.
Gene IDGene SymbolGene NamemiR-139-3p
Total Sites
miR-139-3p
Transfectant Log2 FC
8038ADAM12ADAM metallopeptidase domain 121−4.61
3898LAD1ladinin 11−3.28
144501KRT80keratin 801−2.86
4171MCM2minichromosome maintenance complex component 21−2.84
79019CENPMcentromere protein M1−2.78
147841SPC24SPC24, NDC80 kinetochore complex component1−2.70
4998ORC1origin recognition complex, subunit 11−2.70
26579MYEOVmyeloma overexpressed1−2.57
201266SLC39A11solute carrier family 39, member 111−2.56
9319TRIP13thyroid hormone receptor interactor 131−2.53
493869GPX8glutathione peroxidase 8 (putative)2−2.42
84904ARHGEF39Rho guanine nucleotide exchange factor (GEF) 392−2.35
1734DIO2deiodinase, iodothyronine, type II1−2.27
84733CBX2chromobox homolog 21−2.24
4288MKI67antigen identified by monoclonal antibody Ki-671−2.18
51659GINS2GINS complex subunit 2 (Psf2 homolog)1−2.14
116372LYPD1LY6/PLAUR domain containing 11−2.10
202915TMEM184Atransmembrane protein 184A1−2.10
146909KIF18Bkinesin family member 18B1−1.96
26256CABYRcalcium-binding tyrosine-(Y)-phosphorylation regulated1−1.96
8270LAGE3L antigen family, member 31−1.89
25837RAB26RAB26, member RAS oncogene family1−1.82
8645KCNK5potassium channel, subfamily K, member 53−1.79
8208CHAF1Bchromatin assembly factor 1, subunit B (p60)1−1.78
1763DNA2DNA replication helicase/nuclease 21−1.70
4234METTL1methyltransferase like 11−1.65
1356CPceruloplasmin (ferroxidase)1−1.61
51237MZB1marginal zone B and B1 cell-specific protein1−1.57
6723SRMspermidine synthase1−1.52
51629SLC25A39solute carrier family 25, member 391−1.46
114787GPRIN1G protein regulated inducer of neurite outgrowth 11−1.42
27286SRPX2sushi-repeat containing protein, X-linked 21−1.41
55612FERMT1fermitin family member 11−1.22
7371UCK2uridine-cytidine kinase 23−1.21
1111CHEK1nudix hydrolase 211−1.19
154467CCDC167coiled-coil domain containing 1671−1.18
3070HELLShelicase, lymphoid-specific1−1.15
1515CTSVcathepsin V1−1.14
101ADAM8ADAM metallopeptidase domain 81−1.13
374393FAM111Bfamily with sequence similarity 111, member B1−1.10
7070THY1Thy-1 cell surface antigen1−1.08
128272ARHGEF19Rho guanine nucleotide exchange factor (GEF) 191−1.03
9123SLC16A3solute carrier family 16 (monocarboxylate transporter),
member 3
1−1.02
65268WNK2WNK lysine deficient protein kinase 21−0.97
5163PDK1pyruvate dehydrogenase kinase, isozyme 14−0.96
9833MELKmaternal embryonic leucine zipper kinase1−0.75
1063CENPFcentromere protein F, 350/400 kDa1−0.72
6659SOX4SRY (sex determining region Y)-box 41−0.67
146857SLFN13schlafen family member 131−0.63
51087YBX2Y box binding protein 22−0.62
55502HES6hairy and enhancer of split 6 (Drosophila)1−0.60
51114ZDHHC9zinc finger, DHHC-type containing 91−0.57
Table 2. Significantly enriched annotations of target genes regulated by miR-139-3p.
Table 2. Significantly enriched annotations of target genes regulated by miR-139-3p.
Descriptionp-ValueFDRGenes
DNA replication<0.001<0.001MCM2, CHEK1, CHAF1B,
FAM111B, GINS2, ORC1, DNA2
cell cycle<0.0010.006HELLS, MELK, MCM2, CHEK1,
CHAF1B, CENPF, KIF18B, MKI67, SPC24
positive regulation of cellular extravasation<0.0010.013ADAM8, THY1
integrin-mediated signaling pathway<0.0010.013ADAM12, FERMT1, ADAM8,
THY1
double-strand break repair via break-induced rep-
lication
<0.0010.027MCM2, GINS2
intrinsic apoptotic signaling pathway in response
to oxidative stress
0.0010.041MELK, PDK1
negative regulation of the timing of anagen0.0030.050FERMT1
DNA replication, Okazaki fragment processing0.0030.050DNA2
mitotic spindle assembly checkpoint signaling0.0020.050CENPF, TRIP13
DNA unwinding involved in DNA replication0.0010.050MCM2, GINS2
negative regulation of neuron projection regener-
ation
0.0030.050THY1
positive regulation of fibronectin-dependent thy-
mocyte migration
0.0030.050ADAM8
regulation of transcription from RNA polymerase
II promoter in response to UV-induced DNA
damage
0.0030.050CHEK1
positive regulation of N-terminal peptidyl-lysine
acetylation
0.0030.050SOX4
DNA replication initiation0.0020.050MCM2, ORC1
development of primary sexual characteristics0.0030.050CBX2
exocrine system development0.0030.050RAB26
FDR, false discovery rate.
Table 3. TRIP13-mediated pathways by Gene Set Enrichment Analysis (GSEA).
Table 3. TRIP13-mediated pathways by Gene Set Enrichment Analysis (GSEA).
PathwayEnrichment ScoreNormalized Enrichment Scorep-ValueFDR
KEGG_CELL_CYCLE0.702.60<0.001<0.001
KEGG_DNA_REPLICATION0.752.26<0.001<0.001
KEGG_PROTEASOME0.682.08<0.001<0.001
KEGG_P53_SIGNALING_PATHWAY0.581.96<0.0010.002
KEGG_HOMOLOGOUS_RECOMBINATION0.671.93<0.0010.002
KEGG_MISMATCH_REPAIR0.711.92<0.0010.002
FDR, false discovery rate.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hagihara, Y.; Tomioka, Y.; Suetsugu, T.; Shinmura, M.; Misono, S.; Goto, Y.; Kikkawa, N.; Kato, M.; Inoue, H.; Mizuno, K.; et al. Identification of Tumor-Suppressive miR-139-3p-Regulated Genes: TRIP13 as a Therapeutic Target in Lung Adenocarcinoma. Cancers 2023, 15, 5571. https://doi.org/10.3390/cancers15235571

AMA Style

Hagihara Y, Tomioka Y, Suetsugu T, Shinmura M, Misono S, Goto Y, Kikkawa N, Kato M, Inoue H, Mizuno K, et al. Identification of Tumor-Suppressive miR-139-3p-Regulated Genes: TRIP13 as a Therapeutic Target in Lung Adenocarcinoma. Cancers. 2023; 15(23):5571. https://doi.org/10.3390/cancers15235571

Chicago/Turabian Style

Hagihara, Yoko, Yuya Tomioka, Takayuki Suetsugu, Masahiro Shinmura, Shunsuke Misono, Yusuke Goto, Naoko Kikkawa, Mayuko Kato, Hiromasa Inoue, Keiko Mizuno, and et al. 2023. "Identification of Tumor-Suppressive miR-139-3p-Regulated Genes: TRIP13 as a Therapeutic Target in Lung Adenocarcinoma" Cancers 15, no. 23: 5571. https://doi.org/10.3390/cancers15235571

APA Style

Hagihara, Y., Tomioka, Y., Suetsugu, T., Shinmura, M., Misono, S., Goto, Y., Kikkawa, N., Kato, M., Inoue, H., Mizuno, K., & Seki, N. (2023). Identification of Tumor-Suppressive miR-139-3p-Regulated Genes: TRIP13 as a Therapeutic Target in Lung Adenocarcinoma. Cancers, 15(23), 5571. https://doi.org/10.3390/cancers15235571

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop