Next Article in Journal
Deletion of Mitochondrial Translocator Protein (TSPO) Gene Decreases Oxidative Retinal Pigment Epithelial Cell Death via Modulation of TRPM2 Channel
Previous Article in Journal
Comparative Analysis of Periodontal Pain According to the Type of Precision Orthodontic Appliances: Vestibular, Lingual and Aligners. A Prospective Clinical Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biology
Volume 10
Issue 5
10.3390/biology10050378
biology-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ubiquitin Conjugating Enzyme E2 H (UBE2H) Is Linked to Poor Outcomes and Metastasis in Lung Adenocarcinoma

by
Meng-Chi Yen
1,2,
Kuan-Li Wu
2,3,4,
Yu-Wei Liu
5,
Yung-Yun Chang
6,
Chao-Yuan Chang
2,7,
Jen-Yu Hung
3,4,8,9,10,
Ying-Ming Tsai
2,3,4,9,10,* and
Ya-Ling Hsu
2,4
1
Department of Emergency Medicine, Kaohsiung Medical University Hospital, Kaohsiung Medical University, Kaohsiung 807, Taiwan
2
Graduate Institute of Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung 807, Taiwan
3
Division of Pulmonary and Critical Care Medicine, Kaohsiung Medical University Hospital, Kaohsiung Medical University, Kaohsiung 807, Taiwan
4
Drug Development and Value Creation Research Center, Kaohsiung Medical University, Kaohsiung 807, Taiwan
5
Division of Thoracic Surgery, Department of Surgery, Kaohsiung Medical University Hospital, Kaohsiung Medical University, Kaohsiung 807, Taiwan
6
Division of Internal Medicine, Kaohsiung Medical University Hospital, Kaohsiung Medical University, Kaohsiung 807, Taiwan
7
Department of Anatomy, Kaohsiung Medical University, Kaohsiung 807, Taiwan
8
Department of Internal Medicine, Kaohsiung Municipal Ta-Tung Hospital, Kaohsiung Medical University, Kaohsiung 807, Taiwan
9
Department of Respiratory Care, College of Medicine, Kaohsiung Medical University, Kaohsiung 807, Taiwan
10
Department of Internal Medicine, School of Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung 807, Taiwan
*
Author to whom correspondence should be addressed.
Biology 2021, 10(5), 378; https://doi.org/10.3390/biology10050378
Submission received: 7 April 2021 / Revised: 20 April 2021 / Accepted: 23 April 2021 / Published: 28 April 2021
(This article belongs to the Section Cancer Biology)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

The development of novel treatments for metastatic lung adenocarcinoma is an important issue because some patients do not respond to current standard therapies. Our study aimed to investigate the gene expression profiles in non-tumor tissue, primary tumor tissue, and the metastatic lung tumor tissue in the pleura. After RNA sequencing and bioinformatic analysis from a patient with lung adenocarcinoma (LUAD), ubiquitin conjugating enzyme E2 H (UBE2H) was identified. Compared with normal tissue, a higher expression of UBE2H was observed in the tumor tissue. The high UBE2H expression was significantly associated with poor survival. Suppressing UBE2H in cell lines of lung adenocarcinoma inhibited metastatic capacity and reversed epithelial–mesenchymal transition signaling pathway. Five microRNAs, including miR-101, miR-30a, miR-30b, miR-328, and miR-497, predicted to target UBE2H might be potential prognostic biomarkers for survival in lung adenocarcinoma. The copy number variation may be involved in the regulation of the UBE2H expression. Our observations show that UBE2H is a novel regulatory molecule of metastasis, and may be a prognostic biomarker and a therapeutic target in lung adenocarcinoma.

Abstract

The prognosis of patients with metastatic lung adenocarcinoma (LUAD) is poor. Although novel lung cancer treatments have been developed for metastatic LUAD, not all patients are fit to receive these treatments. The present study aimed to identify the novel regulatory genes in metastatic LUAD. Because the pleural cavity is a frequent metastasis site of LUAD, the adjacent non-tumor tissue, primary tumor tissue, and metastatic lung tumor tissue in the pleura of a single patient with LUAD were collected. The gene expression profiles of the collected samples were further analyzed via RNA sequencing and bioinformatic analysis. A high expression level of ubiquitin conjugating enzyme E2 H (UBE2H), a hypoxia-mediated gene, was identified in the metastatic malignant pleural tumor. After accessing the survival data in patients with lung adenocarcinoma through online databases, a high UBE2H expression was associated with poor survival for LUAD. UBE2H knockdown in two lung adenocarcinoma cell lines suppressed the cell migration capacity and reversed the epithelial–mesenchymal transition (EMT) signaling pathway. A high expression of UBE2H-targeting microRNAs, including miR-101, miR-30a, miR-30b, miR-328, and miR-497, were associated with a favorable prognosis. Moreover, the UBE2H expression revealed a significant correlation with the copy number variation. Taken together, the presence of UBE2H regulated the EMT program and metastasis in LUAD.

1. Introduction

Lung cancer is one of the leading cancer types for the estimated new cancer cases and deaths according to the cancer statistics [1]. Non-small cell lung cancer (NSCLC; approximately 85% of patients) and small cell lung cancer (SCLC; approximately 15% of patients) are the two main forms of lung cancer, and NSCLC is classified into three types, namely adenocarcinoma, squamous cell carcinoma, and large cell lung cancer [2]. A World Health Organization report reveals adenocarcinoma is the most common type of NSCLC [3]. When patients are diagnosed with metastatic NSCLC, the five-year overall survival rate is less than 5% [4]. Immune checkpoint inhibitors (ICIs) and targeted therapies have emerged in the last decade. The specific genetic alterations, including epidermal growth factor receptor (EGFR) gene mutation (30–40% of Asians and 10–20% of Caucasians), echinoderm microtubule associated protein like 4(EML4)-anaplastic lymphoma kinase (ALK) gene fusions (up to 7%), ROS proto-oncogene receptor tyrosine kinase 1 (ROS1) gene rearrangements (1.7%), BRAF V600E gene mutation (2%), and Programmed death-ligand 1 (PD-L1) status, are examined in patients with metastatic diseases in order to choose an adequate ICI [5,6]. These treatments improve the clinical outcomes in some patients with metastatic NSCLC [6]. However, some patients are not able to receive these treatments of they do not have genetic alterations in above targeted genes [6]. Therefore, investigating the molecular mechanisms of metastatic lung cancer is still a crucial issue.
Besides the molecular alterations of NSCLC, the metastasis of NSCLC is driven by various factors. For example, interactions between cancer cells and fibroblasts, stromal cells, and endothelial cells remodels the extracellular matrix [7]. The tumor-infiltrating immune cells secret soluble factors and regulate the immune environment [8]. Hypoxic and metabolic stress also alter the metastatic cascade in malignant cells [9,10]. These factors are thought to trigger the signaling pathways of epithelial–mesenchymal transition (EMT) and to prime the premetastatic niche [11]. There are some limitations in the current targeted therapies and immunotherapies. To develop novel treatments for metastatic NSCLC, investigation of the detailed mechanisms in metastatic NSCLC is essential.
Cancer is one of the causes of pleural effusion [12]. Some soluble molecules in pleural effusion are used as diagnostic markers to identify whether the fluid is malignant pleural effusion [13,14]. The pleural cavity is a frequent metastatic site of lung cancer [15]. Malignant pleural effusion is observed in approximately 50% of patients with advanced NSCLC [16]. The molecular profiling of the mutation status in NSCLC and malignant pleural effusions with next generation sequencing (NGS) provides useful information for cancer management [17,18]. Furthermore, transcriptome analysis in tumor samples via RNA sequencing can identify novel biomarkers and predict signaling pathways. The present study sought to identify novel regulatory genes in metastatic lung adenocarcinoma by utilizing the tissues of the adjacent non-tumor tissue, primary tumor tissue, and metastatic lung tumor tissue in the pleura from a patient. The collected samples were further analyzed for their gene expression profiles via RNA sequencing, bioinformatic analysis, and experimental validation.

2. Materials and Methods

2.1. Sample Collection

The present study was approved by the Institutional Review Board of Kaohsiung Medical University Hospital (IRB no. KMUHIRB-G(II)-20180021). The samples were collected from non-tumor tissue and the tumor tissue of the primary tumor, as well as tumor cells in pleural effusion from an Asian male patient in his 50s with lung adenocarcinoma, staged as T4N0M1a.

2.2. RNA Sequencing

The RNA of each sample was extracted using TRIzol® reagent (Thermo Fisher Scientific, Inc., Madison, WI, USA). The RNA quality, RNA sequencing, small RNA sequencing, sequencing quality trimming, and mRNA alignment were accomplished by Welgene Biotech Co., Ltd. (Welgene Biotech, Taipei, Taiwan). The p-value of the differential genes was calculated using cuffdiff with a non-grouped sample using “blind mode” [19].

2.3. Bioinformatic Analysis

The genes with a differential expression (fold change > 2; p value < 0.05) were subjected to a gene set enrichment analysis (http://www.broad.mit.edu/gsea/; GSEA version 6.2; 1 October 2018) [20,21]. The UBE2H expression in the different stages of the lung adenocarcinoma (LUAD) samples of the Cancer Genome Atlas (TCGA) was analyzed from UALCAN website (http://ualcan.path.uab.edu; 1 November 2020) [22]. The correlation of UBE2H with various proteins was assessed using the Clinical Proteomic Tumor Analysis Consortium (CPTAC) database (https://proteomics.cancer.gov/programs/cptac; 1 November 2020). The association between the UBE2H expression and survival rate, as well as the microRNA expression and survival rate, was analyzed via the Kaplan–Meier plotter (KM plotter) (http://kmplot.com/; 1 October 2018) [23] and SurvExpress (http://bioinformatica.mty.itesm.mx:8080/Biomatec/SurvivaX.jsp; 1 October 2018) [24]. UBE2H targeted microRNAs were predicted through the miRWalk2.0 website (http://zmf.umm.uni-heidelberg.de/apps/zmf/mirwalk2/; 1 October 2018) [25]. The association between the UEB2H expression and each microRNA expression was evaluated via the dataset of TCGA-LUAD.

2.4. Cell Culture and Transfection

Human cell line A549 (ATCC® CCL-185™, Manassas, Virginia, USA) and EML4-ALK Fusion-A549 Isogenic Cell Line Human (ATCC® CCL-185IG™, Manassas, VA, USA) were obtained from the American Type Culture Collection (ATCC; Manassas, VA, USA). Both cells were maintained in an F-12K medium supplied with 10% fetal bovine serum (FBS) and 1% penicillin-streptomycin (Thermo Fisher Scientific, Madison, WI, USA) in a 37 °C, 5% CO2 incubator. The cells were transfected with of 1 μg UBE2H shRNA (clone ID, TRCN0000314525, targeting sequence: 5′-TACGATCCTGGGAGGACTTAA-3′) or control shRNA (non-specific shRNA) and were obtained from the National RNAi Core Facility (Taipei, Taiwan) via lipofectamine 2000 (Invitrogen), according to the manufacturer’s instructions.

2.5. qPCR

The RNA extraction was done using TRIzol reagent (Thermo Fisher Scientific, Waltham, MA, USA). The RNA reverse transcription was performed with a PrimeScript RT reagent kit (Clontech Laboratories, Inc., Mountainview, CA, USA). The PCR product was detected using the specific primers of UBE2H forward, 5′-ATGTACCTCCACCGACCAGA-3′ and reverse, 5′-ACTCCATATCCTGGGCCTCA-3′; GAPDH forward, 5′-GAGTCAACGGATTTGGTCGT-3′ and reverse; and 5′-TTGATTTTGGAGGGATCTCG-3′ on a StepOne Plus Real-Time PCR System (Applied Biosystems; Thermo Fisher Scientific, Inc., Austin, TX, USA) using the Fast SYBR Green Master Mix (Applied Biosystems; Thermo Fisher Scientific, Inc., Austin, TX, USA)). The relative mRNA expression levels were calculated using the 2 −ΔΔCt method [26] and GAPDH served as an internal control.

2.6. Western Blot

The protein lysates of A549 and EML4-ALK fusion-A549 Isogenic cell were collected in a radioimmunoprecipitation lysis (RIPA) buffer with a protease inhibitor (Millipore) and the protein concentration was determined with a bicinchoninic acid (BCA) protein assay kit (Thermo Fisher Scientific, Rockford, IL, USA). Each protein was incubated with a primary antibody, including anti-GAPDH (1:5000, cat. no. MAB374; Millipore), anti- N-cadherin (1:1000, cat. no. 610920), anti-E-cadherin (1:1000, cat. no. 610182), anti-Vimentin (1:1000, cat. no.550513) BD Biosciences, San Jose, CA, USA), anti-Snail (1:1000, cat. no. #3879, Cell Signaling Technology, Danvers, MA, USA), and anti-UBE2H (1:500, cat. no. E2-607, BostonBiochem) at 4 °C overnight. The blots were visualized with the Alpha Innotech FluorChem FC2 imaging system (ProteinSimple; Bio-Techne, Minneapolis, MN, USA).

2.7. Wound-Healing Assay

The cells were seeded into a 24 well plate and a scratch was made by a 200 μL pipette tip when the cells reached a 90–100% confluent monolayer. After removing the cell debris through phosphate-buffered saline washing, the cells were incubated in a culture medium with 2% FBS for 24 h.

2.8. Transwell Migration Assay

The cells were seeded into a QCM™ 24-well Cell Migration Assay and Invasion System (Merck KGaA, Darmstadt, Germany) with uncoated 8 μm pore size polycarbonate membranes (Millipore) in 300 μL of a serum free medium. Then, 500 μL of a medium with 10% FBS was placed in the lower chamber of a 24 well plate. After 24 h, the bottom of the membrane was fixed in 4% formaldehyde and then stained in 1% crystal violet. A cotton swab was used to remove the cells on the upper surface. The stained cells in the bottom of the membrane were then visualized using the Olympus inverted microscopes. Four random fields of view were counted and the relative fold of migration in each group was compared with the control group.

2.9. Statistics

The bar graphs and statistics were performed by GraphPad Prism 7 (GraphPad Software, San Diego, CA, USA). An unpaired, two-tailed Student’s t test was used for analyzing the difference between the two groups. A p-value < 0.05 was considered to indicate a statistically significant difference.

3. Results

3.1. Identification of Ubiquitin Conjugating Enzyme E2 H (UBE2H) in the Malignant Pleural of Lung Adenocarcinoma

The samples were collected from a patient with LUAD. The RNA of the adjacent non-tumor tissue, primary tumor tissue, and metastatic lung tumor tissue in the pleura were performed by RNA sequencing. The genes with a differential expression are shown in Figure 1A. The full gene list is shown in Supplementary Table S1. After the gene set enrichment analysis, the enriched pathways were observed (Figure 1B,C). Hypoxic stress is an important factor to induce the metastasis of cancer [9], and the list and fold change of the genes involved in the hypoxia-related gene set are shown in Table 1. Ubiquitin conjugating enzyme E2 H (UBE2H) is observed in Table 1. Currently, the role of UBE2H is almost unknown in hypoxic-related pathways. However, the correlation analysis showed a high positive association of UBE2H with genes of a hypoxic pathway (AP2B1, BNIP3L, ENOSF1, GPR31, HNRNPC, KCNMA1, KPNB1, LGALS1, LOX, NCL, NDRG1, and RPS3; Figure 1D). The present study focused on the role of UBE2H in LUAD. According to the expression data (fragments per kilo base per million mapped reads (FPKM)) from the RNA sequencing, the highest expression and the lowest expression of UBE2H were observed in the malignant pleural and the primary tumor tissue, respectively (Figure 1E). The UBE2H expression in other patients with LUAD in different tumor stages was evaluated in the LUAD samples of the Cancer Genome Atlas (TCGA). In Figure 1F, a significantly higher UBE2H expression in stages 1 to 4 was observed when compared with the normal tissues (normal vs. stage 1: p-value = 2.7 × 10−13; normal vs. stage 2: p-value = 2.4 × 10−11; normal vs. stage 3: p-value = 2.6 × 10−10; and normal vs. stage 4: p-value = 1.6 × 10−4). In addition, a statistical significance between stage 1 vs. stage 2 was observed (p-value = 3.9 × 10−2), but no statistical significance was observed between the other tumor stages. The results indicated that the UBE2H expression was higher in all LUAD samples when compared with the normal samples.

3.2. Investigation of the Association between UBE2H Expression and Survival in Patients with LUAD

As high expression levels of UBE2H were observed in the tumor samples of LUAD, we further investigated whether the expression levels of UBE2H were associated with the survival of patients through the online database of the KM plotter [23]. The results are shown in Figure 2A. UBE2H detected by three different probes was significantly associated with a poor survival. In another database (the data adapted from GSE31210) [24,27], a similar trend was observed (Figure 2B). This evidence supports that UBE2H may be associated with a poor prognosis in LUAD.

3.3. Investigating the Role of UBE2H in Migration Capacity of LUAD Cell Lines

To determine the role of UBE2H in LUAD, two LUAD cell lines, A549 and EML4-ALK fusion-A549, were used. The EML4-ALK fusion gene, which is identified in approximately 4–6% of LUAD patients, enhances neoplastic phenotypes because the mutation causes continuous ALK activation [28]. EML4-ALK-A549 is derived from the A549 parental cell line containing EML4-ALK fusion variant 1. The effect of UBE2H shRNA is shown in Figure 3A. The wound healing assay and transwell migration assay revealed that the migration capacity of A549 and EML4-ALK-A549 was suppressed after the transfection of UBE2H shRNA (Figure 3B,C). The epithelial–mesenchymal transition (EMT) is an important signaling pathway for regulating cell migration [29]. The key regulatory molecules were detected by Western blot. The mesenchymal phenotype-associated molecules including N-cadherin and Snail in UBE2H knockdown A549 and EML4-ALK-A549 cells was lower than that in the control cells (Figure 3D). By contrast, the expression of Vimentin was not significantly changed. The results suggest that UBE2H knockdown can suppress the metastatic capacity in LUAD through the downregulation of Snail, N-cadherin, and E-cadherin.

3.4. The Potential Regulatory Mechanisms for UBE2H Expression in LUAD

As the upregulation of UBE2H was associated with a poor prognosis and metastatic function in LUAD, the mechanism of UBE2H gene expression was further investigated. Because gene expression can be regulated by microRNA (miRNAs), the potential UBE2H-targeted miRNAs were predicted (Figure 4A). Here, 2390 UBE2H-targeted miRNAs were predicted by the miRWalk database [25]; besides, the results of RNA sequencing in the present study identified 37 downregulated miRNAs (Figure 4A). Of these, 11 shared miRNAs were observed, and the list and fold change are shown in Table 2. Among the 11 miRNAs, a relatively high expression of miR-101, miR-30a, miR-30b, miR-328, and miR-497 was significantly associated with good survival in patients with LUAD through the KM plotter database (Figure 4B). However, only a very weak correlation was observed, even though the p-value reached a statistically significant value after evaluating the expression profiles in the LUAD samples of TCGA (Figure 4C). Furthermore, the UBE2H expression was not correlated with the methylation status of the UBE2H promoter (Figure 4D), but was significantly correlated with the copy number variation (CNV; Figure 4E). These results might suggest that CNV, but not miRNAs or methylation status, is an important factor to up-regulate UBE2H expression.

4. Discussion

Lung cancer has had a notoriously high fatality rate for decades. The novel strategy, RNA sequencing, provides an advanced technique to discover the genes of interest. In this study, we utilized a lung cancer patient with pleural metastasis to discover a gene of interest, UBE2H. Through a public database, a higher expression of UBE2H conferred a poor survival outcome. The progression and invasion of lung cancer were attenuated via the knockdown of UBE2H. Furthermore, the regulatory mechanisms in the UBE2H expression were considered via epigenetic regulation. The findings suggest that UBE2H plays a role in tumor metastasis and can be of therapeutic application in lung cancer.
RNA sequencing was used to identify the novel regulatory genes for the metastatic phenotypes of LUAD. Hundreds of genes with a differential gene expression between primary tumor and malignant pleural were identified, and these genes are involved in multiple biological pathways. The hypoxia-related pathway was focused on it is an important factor affecting the behavior of tumor cells in the tumor microenvironment [9,10]. Among the 14 genes of the hypoxia-related gene set (Table 1), the UBE2H gene was chosen. UBE2H belongs to the ubiquitin conjugating enzyme (UBE2) family, which are involved in ubiquitination and proteasome-mediated protein degradation [30]. There are 40 members in the UBE2 family. Emerging evidence reveals that the aberrant expression of some UBE2 family members is observed in tumor samples. A high expression of UBE2S is associated with malignant breast cancer [31], and an overexpression of UBE2S promotes the proliferation and survival in LUAD cells [32]. The transfection of the UBE2C gene enhances cell proliferation and invasion in lung cancer cells [33], and the expression of UBE2C serves as a prognostic gene for lung adenocarcinoma [34]. The roles of UBE2H in human disease are not well-known. The association of UBE2H and sporadic amyotrophic lateral sclerosis has been reported [35]. Circulating UBE2H mRNA is a potential diagnostic marker or therapeutic target of treatment for Alzheimer’s disease [36]. UBE2H is one of the ubiquitination-related genes used to predict the prognosis of pancreatic cancer patients [37]. In addition, the overexpression of UBE2H enhances proliferation in murine hepatoma cells [38]. Currently, to the best of our knowledge, only one case report indicates that the MET-UBE2H gene fusion is associated with epidermal growth factor receptor (EGFR) resistance in LUAD [39]. Therefore, UBE2H was focused on in the present study.
The RNA sequencing data showed a relatively high UBE2H expression in the metastatic lung tumor tissue in the pleura (Figure 1D). In Figure 1E, the UBE2H expressions in stages 1 to 4 of the LUAD samples of TCGA are higher than in the normal samples. A high expression of UBE2H was associated with poor survival in patients with LUAD (Figure 2). Moreover, the UBE2H expression in stage 2 was significantly higher than that in stage 1 lung cancer. Although the expression of UBE2H did not reveal a statistical difference between stages 3 or 4 vs. stage 1, the mean UBE2H expression in stages 3 and 4 was higher compared with stage 1 (Figure 1E). This might imply that the upregulation of UBE2H is associated with tumor metastasis, even though it is in the early stage of metastasis. UBE2H has been observed in 75-gene signatures associated with lymph node metastasis in lung adenocarcinoma [40]. UBE2H has been identified in metastasis-associated gene sets in colorectal cancer [41]. Our results demonstrated that suppressing the UBE2H expression of A549 and EML4-ALK-A549 reduced the metastatic capacity. The regulatory molecules of the EMT pathways, including Snail, N-cadherin, and E-cadherin, were affected after silencing UBE2H (Figure 3). The ubiquitin–proteasome signaling pathway is an important factor for regulating the EMT phenotype [42], but how UBE2H affects EMT signaling is still unknown. This detailed mechanism should be investigated in the future.
miRNAs are involved in the regulation of the gene expression and various biological processes including metastasis [43]. According to the RNA sequencing in the present study and candidate UBE2H targeting miRNAs, we observed that the expression miR-101, miR-30a, miR-30b, miR-328, and miR-497 was associated with a good survival. mir-101, mir-30a, and miR-30b have been demonstrated to be suppressor miRNAs for lymph node metastasis in NSCLC [44,45]. However, no significant correlation between the UBE2H expression and miRNAs were observed in the LUAD samples of TCGA (Figure 4C,D). Our results implied that the copy number variation is a potential factor to increase UBE2H expression in TCGA samples (Figure 4E). Currently, no further experimental evidence or studies provide a regulatory mechanism between the copy number variation and the expression of UBE2H. This should be a critical issue to be investigated in the future.

5. Conclusions

Collectively, the present study identified a novel hypoxia-related gene UBE2H from the malignant pleural of a patient with LUAD, investigated the role of UBE2H in metastatic capacity in two LUAD cell lines, and predicted a potential mechanism for the regulation of the UBE2H expression (Figure 5). The current findings demonstrate that UBE2H serves as a potential prognostic biomarker and a therapeutic target for metastatic LUAD.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/biology10050378/s1, Table S1. The full list of genes with significant changes.

Author Contributions

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

Funding

This study was supported by grants from the Ministry of Science and Technology (grant nos. MOST 108-2320-B-037-024 -MY3, MOST 108-2314-B-037-094-, MOST 108-2314-B-037-095-, and 107-2314-B-037-107-MY3) and Kaohsiung Medical University (grant nos. KMU-DK 109002, KMU-DK108003, KMU-DK108008, and KMU-TC108A03-5).

Institutional Review Board Statement

The present study was approved by the Institutional Review Board of Kaohsiung Medical University Hospital (IRB no. KMUHIRB-G (II)-20180021).

Informed Consent Statement

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

Acknowledgments

We thank the National RNAi Core Facility at Academia Sinica in Taiwan for providing shRNA reagents and related services. The authors thank the Center for Research Resources and Development of Kaohsiung Medical University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics, 2020. CA Cancer J. Clin. 2020, 70, 7–30. [Google Scholar] [CrossRef]
  2. Duma, N.; Santana-Davila, R.; Molina, J.R. Non-Small Cell Lung Cancer: Epidemiology, Screening, Diagnosis, and Treatment. Mayo Clin. Proc. 2019, 94, 1623–1640. [Google Scholar] [CrossRef]
  3. Travis, W.D.; Brambilla, E.; Burke, A.P.; Marx, A.; Nicholson, A.G. Introduction to The 2015 World Health Organization Classification of Tumors of the Lung, Pleura, Thymus, and Heart. J. Thorac. Oncol. 2015, 10, 1240–1242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Arbour, K.C.; Riely, G.J. Systemic Therapy for Locally Advanced and Metastatic Non-Small Cell Lung Cancer: A Review. JAMA 2019, 322, 764–774. [Google Scholar] [CrossRef] [PubMed]
  5. Elkrief, A.; Joubert, P.; Florescu, M.; Tehfe, M.; Blais, N.; Routy, B. Therapeutic landscape of metastatic non-small-cell lung cancer in Canada in 2020. Curr. Oncol. 2020, 27, 52–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Reckamp, K.L. Targeted Therapy for Patients With Metastatic Non-Small Cell Lung Cancer. J. Natl. Compr. Cancer Netw. 2018, 16, 601–604. [Google Scholar] [CrossRef] [PubMed]
  7. Tarbe, N.; Losch, S.; Burtscher, H.; Jarsch, M.; Weidle, U.H. Identification of rat pancreatic carcinoma genes associated with lymphogenous metastasis. Anticancer Res. 2002, 22, 2015–2027. [Google Scholar] [PubMed]
  8. Kitamura, T.; Qian, B.Z.; Pollard, J.W. Immune cell promotion of metastasis. Nat. Rev. Immunol. 2015, 15, 73–86. [Google Scholar] [CrossRef]
  9. Rankin, E.B.; Giaccia, A.J. Hypoxic control of metastasis. Science 2016, 352, 175–180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Rankin, E.B.; Nam, J.M.; Giaccia, A.J. Hypoxia: Signaling the Metastatic Cascade. Trends Cancer 2016, 2, 295–304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Fares, J.; Fares, M.Y.; Khachfe, H.H.; Salhab, H.A.; Fares, Y. Molecular principles of metastasis: A hallmark of cancer revisited. Signal. Transduct. Target. Ther. 2020, 5, 28. [Google Scholar] [CrossRef]
  12. Jany, B.; Welte, T. Pleural Effusion in Adults-Etiology, Diagnosis, and Treatment. Dtsch. Arztebl. Int. 2019, 116, 377–386. [Google Scholar] [CrossRef] [PubMed]
  13. Wu, D.W.; Chang, W.A.; Liu, K.T.; Yen, M.C.; Kuo, P.L. Vascular endothelial growth factor and protein level in pleural effusion for differentiating malignant from benign pleural effusion. Oncol. Lett. 2017, 14, 3657–3662. [Google Scholar] [CrossRef] [Green Version]
  14. Epelbaum, O.; Rahman, N.M. Contemporary approach to the patient with malignant pleural effusion complicating lung cancer. Ann. Transl Med. 2019, 7, 352. [Google Scholar] [CrossRef] [PubMed]
  15. Agalioti, T.; Giannou, A.D.; Stathopoulos, G.T. Pleural involvement in lung cancer. J. Thorac. Dis. 2015, 7, 1021–1030. [Google Scholar] [CrossRef]
  16. Roberts, M.E.; Neville, E.; Berrisford, R.G.; Antunes, G.; Ali, N.J.; Group BTSPDG. Management of a malignant pleural effusion: British Thoracic Society Pleural Disease Guideline 2010. Thorax 2010, 65 (Suppl. 2), ii32–ii40. [Google Scholar] [CrossRef] [Green Version]
  17. Liu, L.; Shao, D.; Deng, Q.; Tang, H.; Wang, J.; Liu, J.; Guo, F.; Lin, Y.; Peng, Z.; Mao, M.; et al. Next generation sequencing-based molecular profiling of lung adenocarcinoma using pleural effusion specimens. J. Thorac. Dis. 2018, 10, 2631–2637. [Google Scholar] [CrossRef]
  18. Grigoriadou, G.; Esagian, S.M.; Ryu, H.S.; Nikas, I.P. Molecular Profiling of Malignant Pleural Effusions with Next Generation Sequencing (NGS): Evidence that Supports Its Role in Cancer Management. J. Pers. Med. 2020, 10, 206. [Google Scholar] [CrossRef] [PubMed]
  19. Trapnell, C.; Roberts, A.; Goff, L.; Pertea, G.; Kim, D.; Kelley, D.R.; Pimentel, H.; Salzberg, S.L.; Rinn, J.L.; Pachter, L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 2012, 7, 562–578. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. 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] [CrossRef] [Green Version]
  21. Mootha, V.K.; Lindgren, C.M.; Eriksson, K.F.; Subramanian, A.; Sihag, S.; Lehar, J.; Puigserver, P.; Carlsson, E.; Ridderstrale, M.; Laurila, E.; et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat. Genet. 2003, 34, 267–273. [Google Scholar] [CrossRef]
  22. Chandrashekar, D.S.; Bashel, B.; Balasubramanya, S.A.H.; Creighton, C.J.; Ponce-Rodriguez, I.; Chakravarthi, B.; Varambally, S. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia 2017, 19, 649–658. [Google Scholar] [CrossRef] [PubMed]
  23. Gyorffy, B.; Surowiak, P.; Budczies, J.; Lanczky, A. Online survival analysis software to assess the prognostic value of biomarkers using transcriptomic data in non-small-cell lung cancer. PLoS ONE 2013, 8, e82241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Aguirre-Gamboa, R.; Gomez-Rueda, H.; Martinez-Ledesma, E.; Martinez-Torteya, A.; Chacolla-Huaringa, R.; Rodriguez-Barrientos, A.; Tamez-Pena, J.G.; Trevino, V. SurvExpress: An online biomarker validation tool and database for cancer gene expression data using survival analysis. PLoS ONE 2013, 8, e74250. [Google Scholar] [CrossRef] [Green Version]
  25. Dweep, H.; Gretz, N. miRWalk2.0: A comprehensive atlas of microRNA-target interactions. Nat. Methods 2015, 12, 697. [Google Scholar] [CrossRef]
  26. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  27. Okayama, H.; Kohno, T.; Ishii, Y.; Shimada, Y.; Shiraishi, K.; Iwakawa, R.; Furuta, K.; Tsuta, K.; Shibata, T.; Yamamoto, S.; et al. Identification of genes upregulated in ALK-positive and EGFR/KRAS/ALK-negative lung adenocarcinomas. Cancer Res. 2012, 72, 100–111. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Bayliss, R.; Choi, J.; Fennell, D.A.; Fry, A.M.; Richards, M.W. Molecular mechanisms that underpin EML4-ALK driven cancers and their response to targeted drugs. Cell. Mol. Life Sci. 2016, 73, 1209–1224. [Google Scholar] [CrossRef] [Green Version]
  29. Xiao, D.; He, J. Epithelial mesenchymal transition and lung cancer. J. Thorac. Dis. 2010, 2, 154–159. [Google Scholar] [CrossRef]
  30. Clague, M.J.; Heride, C.; Urbe, S. The demographics of the ubiquitin system. Trends Cell Biol. 2015, 25, 417–426. [Google Scholar] [CrossRef]
  31. Ayesha, A.K.; Hyodo, T.; Asano, E.; Sato, N.; Mansour, M.A.; Ito, S.; Hamaguchi, M.; Senga, T. UBE2S is associated with malignant characteristics of breast cancer cells. Tumour Biol. 2016, 37, 763–772. [Google Scholar] [CrossRef]
  32. Liu, Z.; Xu, L. UBE2S promotes the proliferation and survival of human lung adenocarcinoma cells. BMB Rep. 2018, 51, 642–647. [Google Scholar] [CrossRef]
  33. Tang, X.K.; Wang, K.J.; Tang, Y.K.; Chen, L. Effects of ubiquitin-conjugating enzyme 2C on invasion, proliferation and cell cycling of lung cancer cells. Asian Pac. J. Cancer Prev. 2014, 15, 3005–3009. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Jiang, Y.; Huang, Y.; Du, Y.; Zhao, Y.; Ren, J.; Ma, S.; Wu, C. Identification of Prognostic Genes and Pathways in Lung Adenocarcinoma Using a Bayesian Approach. Cancer Inform. 2017, 16, 1176935116684825. [Google Scholar] [CrossRef] [PubMed]
  35. Martin, I.; Vourc’h, P.; Mahe, M.; Thepault, R.A.; Antar, C.; Vedrine, S.; Praline, J.; Camu, W.; Andres, C.R.; Corcia, P.; et al. Association study of the ubiquitin conjugating enzyme gene UBE2H in sporadic ALS. Amyotroph. Lateral Scler. 2009, 10, 432–435. [Google Scholar] [CrossRef]
  36. Lim, K.H.; Joo, J.Y. Predictive Potential of Circulating Ube2h mRNA as an E2 Ubiquitin-Conjugating Enzyme for Diagnosis or Treatment of Alzheimer’s Disease. Int. J. Mol. Sci 2020, 21, 3398. [Google Scholar] [CrossRef] [PubMed]
  37. Zuo, H.; Chen, L.; Li, N.; Song, Q. Identification of a Ubiquitination-Related Gene Risk Model for Predicting Survival in Patients With Pancreatic Cancer. Front. Genet. 2020, 11, 612196. [Google Scholar] [CrossRef] [PubMed]
  38. Keng, V.W.; Villanueva, A.; Chiang, D.Y.; Dupuy, A.J.; Ryan, B.J.; Matise, I.; Silverstein, K.A.; Sarver, A.; Starr, T.K.; Akagi, K.; et al. A conditional transposon-based insertional mutagenesis screen for genes associated with mouse hepatocellular carcinoma. Nat. Biotechnol. 2009, 27, 264–274. [Google Scholar] [CrossRef]
  39. Zhu, Y.C.; Wang, W.X.; Song, Z.B.; Zhang, Q.X.; Xu, C.W.; Chen, G.; Zhuang, W.; Lv, T.; Song, Y. MET-UBE2H Fusion as a Novel Mechanism of Acquired EGFR Resistance in Lung Adenocarcinoma. J. Thorac. Oncol. 2018, 13, e202–e204. [Google Scholar] [CrossRef] [Green Version]
  40. Inamura, K.; Shimoji, T.; Ninomiya, H.; Hiramatsu, M.; Okui, M.; Satoh, Y.; Okumura, S.; Nakagawa, K.; Noda, T.; Fukayama, M.; et al. A metastatic signature in entire lung adenocarcinomas irrespective of morphological heterogeneity. Hum. Pathol. 2007, 38, 702–709. [Google Scholar] [CrossRef]
  41. Li, X.; Peng, S. Identification of metastasis-associated genes in colorectal cancer through an integrated genomic and transcriptomic analysis. Chin. J. Cancer Res. 2013, 25, 623–636. [Google Scholar] [CrossRef]
  42. Voutsadakis, I.A. The ubiquitin-proteasome system and signal transduction pathways regulating Epithelial Mesenchymal transition of cancer. J. Biomed. Sci 2012, 19, 67. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Kim, J.; Yao, F.; Xiao, Z.; Sun, Y.; Ma, L. MicroRNAs and metastasis: Small RNAs play big roles. Cancer Metastasis Rev. 2018, 37, 5–15. [Google Scholar] [CrossRef] [PubMed]
  44. Wu, S.G.; Chang, T.H.; Liu, Y.N.; Shih, J.Y. MicroRNA in Lung Cancer Metastasis. Cancers (Basel) 2019, 11, 265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Han, L.; Chen, W.; Xia, Y.; Song, Y.; Zhao, Z.; Cheng, H.; Jiang, T. MiR-101 inhibits the proliferation and metastasis of lung cancer by targeting zinc finger E-box binding homeobox 1. Am. J. Transl. Res. 2018, 10, 1172–1183. [Google Scholar]
Biology 10 00378 g001 550
Figure 1. Identification of ubiquitin conjugating enzyme E2 H (UBE2H) from lung adenocarcinoma (LUAD) tumor samples. (A) The volcano plot (log2 fold-change vs. log10 p-value) of upregulated and downregulated genes showed the differential expression distribution for the sample pair. The green and red dots indicate significant down- and up-regulated genes, respectively. (B) The enriched biological processes by gene set enrichment analysis. (C) The histogram plot of hypoxia-related gene set. (D) The correlation of UBE2H with the genes of the hypoxic pathway. (E) The expression of UBE2H in the samples collected in the present study. (F) The expression of UBE2H in the lung adenocarcinoma (LUAD) samples of the Cancer Genome Atlas (TCGA). NES-normalized enrichment score; FDR-false discovery rate; FPKM-fragments per kilobase of transcript per million.
Figure 1. Identification of ubiquitin conjugating enzyme E2 H (UBE2H) from lung adenocarcinoma (LUAD) tumor samples. (A) The volcano plot (log2 fold-change vs. log10 p-value) of upregulated and downregulated genes showed the differential expression distribution for the sample pair. The green and red dots indicate significant down- and up-regulated genes, respectively. (B) The enriched biological processes by gene set enrichment analysis. (C) The histogram plot of hypoxia-related gene set. (D) The correlation of UBE2H with the genes of the hypoxic pathway. (E) The expression of UBE2H in the samples collected in the present study. (F) The expression of UBE2H in the lung adenocarcinoma (LUAD) samples of the Cancer Genome Atlas (TCGA). NES-normalized enrichment score; FDR-false discovery rate; FPKM-fragments per kilobase of transcript per million.
Biology 10 00378 g001
Biology 10 00378 g002 550
Figure 2. The correlation of UBE2H expression and survival in patients with LUAD. (A) The survival analysis via the KM plotter. The samples are divided into high and low expressions according to median expression value of UBE2H. The red and black lines indicate high and low UBE2H expressions, respectively. Three plots are analyzed by the different probes labeled on each plot. (B) The survival analysis via GSE31210 dataset on SurvExpress. The red and green lines indicate high and low UBE2H expressions, respectively. HR—hazard ratios.
Figure 2. The correlation of UBE2H expression and survival in patients with LUAD. (A) The survival analysis via the KM plotter. The samples are divided into high and low expressions according to median expression value of UBE2H. The red and black lines indicate high and low UBE2H expressions, respectively. Three plots are analyzed by the different probes labeled on each plot. (B) The survival analysis via GSE31210 dataset on SurvExpress. The red and green lines indicate high and low UBE2H expressions, respectively. HR—hazard ratios.
Biology 10 00378 g002
Biology 10 00378 g003 550
Figure 3. The effect of UBE2H silencing on the cell migratory capacity in LUAD cell lines. (A) The effect of UBE2H shRNA is examined in A549. (B) Wound-healing assay (n = 3). (C) Transwell migration assay and quantitative results (n = 3). (D) Western blot assay for determining the regulatory molecules of the epithelial–mesenchymal transition signaling pathway (n = 3). * p-value < 0.05, as compared with the vector control.
Figure 3. The effect of UBE2H silencing on the cell migratory capacity in LUAD cell lines. (A) The effect of UBE2H shRNA is examined in A549. (B) Wound-healing assay (n = 3). (C) Transwell migration assay and quantitative results (n = 3). (D) Western blot assay for determining the regulatory molecules of the epithelial–mesenchymal transition signaling pathway (n = 3). * p-value < 0.05, as compared with the vector control.
Biology 10 00378 g003
Biology 10 00378 g004 550
Figure 4. The potential regulatory mechanism of UBE2H. (A) 11 shared genes are identified from the 2390 candidate UBE2H-targeted miRNAs (blue circle) predicted by miRWalk and 37 down-regulated miRNAs from RNA sequencing (green circle). (B) The survival analysis via the KM plotter. The samples are divided by the median expression value of each miRNA. (C) The correlation between the expression of UBE2H and the expression of each miRNA. (D) The correlation between the expression of UBE2H and the methylation status of the UBE2H promoter. (E) The correlation between the expression of UBE2H and the copy number variation (CNV) among the LUAD samples of TCGA.
Figure 4. The potential regulatory mechanism of UBE2H. (A) 11 shared genes are identified from the 2390 candidate UBE2H-targeted miRNAs (blue circle) predicted by miRWalk and 37 down-regulated miRNAs from RNA sequencing (green circle). (B) The survival analysis via the KM plotter. The samples are divided by the median expression value of each miRNA. (C) The correlation between the expression of UBE2H and the expression of each miRNA. (D) The correlation between the expression of UBE2H and the methylation status of the UBE2H promoter. (E) The correlation between the expression of UBE2H and the copy number variation (CNV) among the LUAD samples of TCGA.
Biology 10 00378 g004
Biology 10 00378 g005 550
Figure 5. The proposed mechanisms for how UBE2H regulates lung cancer EMT and metastasis and UBE2H is regulated.
Figure 5. The proposed mechanisms for how UBE2H regulates lung cancer EMT and metastasis and UBE2H is regulated.
Biology 10 00378 g005
Table
Table 1. The expression of the specific genes in the hypoxia-related gene set.
Table 1. The expression of the specific genes in the hypoxia-related gene set.
Gene SymbolLog2 (Fold) (Pleural Malignancy/Primary Malignancy)Gene SymbolLog2 (Fold) (Pleural Malignancy/Primary Malignancy)
LGALS118.337 UBE2H4.889
AP2B14.034ENOSF13.074
NCL2.996KPNB12.816
RPS32.768KCNMA12.708
LOX2.304HNRNPC2.286
HIF1A2.114BNIP3L2.087
NDRG11.963GPR391.770
Table
Table 2. The candidate miRNAs regulating UBE2H.
Table 2. The candidate miRNAs regulating UBE2H.
miRNAFold (Pleural Malignant/Primary Malignant)miRNAFold (Pleural Malignant/Primary Malignant)
hsa-miR-101-3p0.442752hsa-miR-122-5p0.000752
hsa-miR-12750.486842hsa-miR-135a-5p0.035668
hsa-miR-143-3p0.323242hsa-miR-30a-3p0.167051
hsa-miR-30b-3p0.410777hsa-miR-3150b-5p0.000195
hsa-miR-497-3p0.142195hsa-miR-328-3p0.477685
hsa-miR-6724-5p0.413534
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yen, M.-C.; Wu, K.-L.; Liu, Y.-W.; Chang, Y.-Y.; Chang, C.-Y.; Hung, J.-Y.; Tsai, Y.-M.; Hsu, Y.-L. Ubiquitin Conjugating Enzyme E2 H (UBE2H) Is Linked to Poor Outcomes and Metastasis in Lung Adenocarcinoma. Biology 2021, 10, 378. https://doi.org/10.3390/biology10050378

AMA Style

Yen M-C, Wu K-L, Liu Y-W, Chang Y-Y, Chang C-Y, Hung J-Y, Tsai Y-M, Hsu Y-L. Ubiquitin Conjugating Enzyme E2 H (UBE2H) Is Linked to Poor Outcomes and Metastasis in Lung Adenocarcinoma. Biology. 2021; 10(5):378. https://doi.org/10.3390/biology10050378

Chicago/Turabian Style

Yen, Meng-Chi, Kuan-Li Wu, Yu-Wei Liu, Yung-Yun Chang, Chao-Yuan Chang, Jen-Yu Hung, Ying-Ming Tsai, and Ya-Ling Hsu. 2021. "Ubiquitin Conjugating Enzyme E2 H (UBE2H) Is Linked to Poor Outcomes and Metastasis in Lung Adenocarcinoma" Biology 10, no. 5: 378. https://doi.org/10.3390/biology10050378

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