Abstract Background Lung adenocarcinoma (LUAD) is the primary subtype of human lung cancer. The effectiveness of treatment and long-term survival of patients with LUAD are current suboptimal. Tripartite motif containing 56 (TRIM56) is a member of the TRIM protein family that have functions predominantly in immunity and cancer. Purpose To investigate the expression of TRIM56 in LUAD, and explore the potential regulatory role of TRIM56 in the invasion and migration of LUAD cells. Methods The Gene Expression Omnibus datasets and The Cancer Genome Atlas-LUAD cohort were used to analyze the mRNA expression of TRIM56 in LUAD. The differential expression profiles of miRNAs associated with TRIM56 were obtained from The Cancer Genome Atlas-LUAD cohort. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed to determine the principal functions of miRNAs and interacting proteins. Transwell and wound healing were used to detect the effect of overexpression of TRIM56 on the invasion and migration of LUAD cells. Results The expression of TRIM56 was decreased in LUAD and associated with poor prognosis. We determined the genome copy number, negatively correlated miRNA and potential transcription factors of TRIM56, and conducted enrichment analysis. Among them, hsa-mir-542 and hsa-mir-627 were the most likely to inhibit the expression of TRIM56. We also predicted the interacting proteins and potential ubiquitination substrate of TRIM56. Finally, we demonstrated that overexpression of TRIM56 inhibits the invasion and migration of LUAD cells. Conclusion This study is the first to analyze the expression of TRIM56 and its inhibitory effect on the invasion and migration of LUAD. This evidence provides a new direction for further study of the reasons for the low expression of TRIM56 in LUAD and its regulatory mechanism. Keywords: TRIM56, lung adenocarcinoma, bioinformatics analysis, invasion, migration Introduction Lung cancer is one of the most common tumors and associated with the highest rates of morbidity and mortality worldwide.[30]^1 In 2017, the death toll due to lung cancer exceeded the total caused by prostate, colorectal, brain, and breast cancer. Almost one-quarter of cancer-related deaths are due to lung cancer.[31]^2 Epidemiological statistics have shown that lung adenocarcinoma (LUAD) has become the primary subtype of lung cancer.[32]^2 Based on its characteristics, LUAD is prone to recurrence and metastasis; in addition, the therapeutic effect in patients with LUAD is currently suboptimal, and their long-term survival rate remains low.[33]^3 In general, the occurrence and development of tumors is a complex, multi-stage process, involving numerous signal transduction pathways, and a variety of factors in vitro and in vivo. Therefore, investigation of the abnormally altered genes in LUAD is crucial to further elucidate the relevant mechanisms involved in its occurrence and development. The human genome encodes more than 70 members of the tripartite motif containing (TRIM) family. In previous reports, almost all TRIM proteins exhibited E3 ubiquitin ligase activity depending on their RING finger domain, which mediates the ubiquitination of various substrates and themselves.[34]^4 TRIM proteins are involved in various cell processes, including differentiation, proliferation, apoptosis, and morphogenesis. Moreover, evidence on its role in immune signal transduction and antiviral functions is emerging.[35]^5^,[36]^6 TRIM56 is an important member of the TRIM family. At present, research on TRIM56 mainly focuses on the antiviral activity TRIM56 can exert by mediating the K63 ubiquitin link of stimulator of interferon genes (STING) (a type I interferon stimulating gene), as well as the induction process of type I interferon by double-stranded DNA of viruses and other pathogens.[37]^7 TRIM56 can also activate the toll like receptor 3 (TLR3) signaling pathway and promote the production of interferon. Of note, this process is independent of its E3 ubiquitin ligase activity.[38]^8^,[39]^9 In addition, TRIM56 can also inhibit the RNA replication of influenza viruses A and B through its specific C-terminal functional domain.[40]^10 However, there are few studies on TRIM56 in tumors. TRIM56 inhibits the proliferation of multiple myeloma cells and induces apoptosis.[41]^11 In ovarian cancer, TRIM56 can mediate the K48 ubiquitin connection of vimentin to promote its degradation, thus inhibiting the development of ovarian cancer.[42]^12 Interestingly, in breast cancer, TRIM56 mediates K63 ubiquitination with estrogen receptor-alpha (ER-α), upregulates the protein stability of ER-α, and promotes the proliferation of breast cancer cells.[43]^13 Therefore, the role of TRIM56 in tumors remains unclear, and alterations in its expression and function in lung cancer have not been reported.In our study, we will explore the relationship between abnormal expression of TRIM56 and lung adenocarcinoma through bioinformatics analysis, and verify that TRIM56 overexpression suppress the invasion and migration of lung adenocarcinoma cells. Materials and Methods Cells Culture and Plasmids H1299 cells were obtained from the American Type Culture Collection (Manassas, VA, USA), while 293T cells for cell transfection purchased from Gene Pharma (China). The cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM; high glucose) supplemented with 10% fetal bovine serum (Bioind, USA) and 1% penicillin/streptomycin (Solarbio, Beijing, China) in a humidified atmosphere with 5% CO[2] at 37°C. The transfection reagent Lipofectamine 2000 was purchased from Vigorous Biotechnology (China). Plasmids for cell transfection were provided and constructed by a laboratory. Western Blotting H1299 cells were collected and lysed with radioimmunoprecipitation assay buffer (Beyotime, Shanghai, China). Total cellular protein was collected using 5 × sodium dodecyl sulfate sample buffer and denatured for 10 min at 100°C. Proteins were separated by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis and electrotransferred to polyvinylidene difluoride membranes. After blocking with 5% skimmed milk (weight/volume), the polyvinylidene difluoride membranes were incubated with the following antibodies: mouse monoclonal anti-Myc antibody (Cell Signaling Technology, Danvers, MA, USA) and monoclonal rabbit anti-glyceraldehyde-3-phosphate dehydrogenase (Proteintech, Rosemont, IL, USA). Subsequently, the membranes were incubated with mouse and rabbit secondary antibodies for 1.5 h. A highly sensitive substrate (Millipore, Billerica, MA, USA) was used to detect specific bands on autoradiographic film. Lentivirus Packaging and Infection pCDH-Myc-TRIM56 or pCDH empty vector and two packaging plasmids (delta 8.9, pLP-VSVG) were stably transfected into 293T cells using the Lipofectamine 2000 transfection reagent (Thermo Scientific, Waltham, MA, USA). The concentration of lentivirus was determined after 48 h and 72 h using filtered virus supernatant. H1299 cells were infected with virus concentrate (5 μL) for 48 h, and the stability of transfection was assessed using puromycin. Polybrene (8 μg/mL) was used to infect cells with viruses. UALCAN Analysis UALCAN ([44]http://ualcan.path.uab.edu) is an open interactive website. The website includes RNA sequencing, clinical research, and DNA methylation data regarding 31 types of cancer, which are available in The Cancer Genome Atlas (TCGA) database.[45]^14 UALCAN can analyze the relative expression of genes in tumors and normal samples, as well as various tumor subgroups based on clinicopathological characteristics. We used UALCAN to analyze the relationship between the expression of TRIM56 and the clinical characteristics of LUAD, including the patient’s age, race, sex, lymph node metastasis, tumor stage, and tumor protein p53 (TP53) mutation. Furthermore, we analyzed the relationship between promoter methylation and TRIM56 expression in patients with LUAD. Gene Expression Profiling Interactive Analysis (GEPIA) Analysis GEPIA ([46]http://gepia.cancer-pku.cn/) is a web-based tool,[47]^15 providing customizable analysis functions. These functions include gene expression analysis in tumors and normal tissues, multi-gene expression correlation analysis, and patient survival analysis. The main data sources of GEPIA are TCGA and Genotype-Tissue Expression databases. We used GEPIA to analyze the expression of TRIM56 in three subtypes of LUAD, namely terminal respiratory unit, proximal hyperplasia, and proximal inflammatory. LinkedOmics Analysis LinkedOmics is a publicly available portal ([48]http://www.LinkedOmics.org/login.php) that includes multi-omics data regarding 32 types of cancer, which are available in TCGA.[49]^16 We used the LinkFinder module of LinkedOmics to explore the differentially expressed miRNA associated with TRIM56 in the TCGA-LUAD cohort. Pearson correlation coefficients were determined to explore the negative correlations with TRIM56. All results are graphically displayed in the volcano and heat map. A false discover ratio <0.05 was set as the threshold. cBioPortal Analysis The cBioPortal database ([50]http://cbioportal.org) is an open-source tool for the interactive exploration of multi-dimensional cancer genomics datasets and currently contains evidence obtained from 308 cancer studies.[51]^17^,[52]^18 We used the cBioPortal to analyze TRIM56 copy number alteration (CNA), mutations, and changes in DNA methylation using data from TCGA-LUAD. Functional Enrichment (FunRich) Analysis Tool We performed enrichment analysis of negatively correlated miRNAs using the FunRich tool ([53]http://www.funrich.org/). FunRich is a stand-alone software tool used mainly for functional enrichment and interaction network analysis of genes and proteins.[54]^19^,[55]^20 We used FunRich to analyze the biological process, cellular component, molecular function, and biological pathways. The Kaplan–Meier Plotter The prognostic value of the expression of TRIM56 in [56]GSE31210 and negatively correlated miRNAs was evaluated using the online database Kaplan–Meier Plotter ([57]http://www.kmplot.com).[58]^21^,[59]^22 To analyze the survival of patients with LUAD, patient samples were classified into two groups according to their median expression (high versus low expression). A Kaplan–Meier survival plot was produced, and the results included hazard ratios with 95% confidence intervals and log-rank p-values. Search Tool for the Retrieval of Interacting Genes (STRING) STRING is a protein–protein interaction database, involving known and predicted data ([60]https://string-db.org/). The protein interaction data available in STRING include direct (physical) and indirect (functional) associations; they stem from computational prediction, knowledge transfer between organisms, and interactions aggregated from other (primary) databases.[61]^23^,[62]^24 We used STRING to analyze proteins interacting with TRIM56 and the enrichment analysis data. Transwell Assays Corning Matrigel was purchased from Becton, Dickinson and Company (USA). H1299 cells were stably transfected with pCDH-Myc-TRIM56 or pCDH empty vector, and seeded at a density of 1×10^4 cells per chamber in 24-well plates. DMEM medium with 0.1% bovine serum albumin was placed in the upper chamber. DMEM medium with 25% fetal bovine serum was added in the lower chamber. Wound Healing Assays H1299 cells were stably transfected with pCDH-Myc-TRIM56 or pCDH empty vector, and seeded at a density of 1×10^5 cells per well in six-well plates. After drawing lines for 0 h, the DMEM high glucose medium containing 100× penicillin-streptomycin solution was replaced. qRT-PCR Total RNA was isolated using TRIzol reagent (Invitrogen, Carlsbad, CA). First strand cDNA was synthesized subsequently using the Fast Quant RT Kit (with gDNase) (TOYOBO,Japan) according to the product manual. Quantitative polymerase chain reaction (qPCR) was performed in three replicate wells on an ABI 7500 Real-Time PCR System (ThermoFisher Scientific, Waltham, MA) using SuperReal PreMix Plus (SYBR Green) (Tiangen Biotech, Beijing, China). GAPDH, forward:CATGAGAAGTATGACAACAGCCT; reverse: AGTCCTTCCACGATACAAAG. TRIM56,forward:GATCAGCTCTGGCTAGTTCTCACA;reverse:CTAGTCTTCTGAGGCCTTGGGC. The relative expression levels of TRIM56 were calculated using the 2^−ΔΔCT method and GAPDH was used as an internal control. Each experiment was repeated three times. Statistical Analysis The cut-off values for TRIM56 expression were determined by their median values. All analyses, including the t-test and correlation analysis were performed by GraphPad Prism. In our analysis, P values less than 0.05 were considered significant. Results TRIM56 Was Downregulated in LUAD and Associated with Poor Prognosis To determine the expression levels of TRIM56 in LUAD, we first analyzed its mRNA expression by searching the open databases Gene Expression Omnibus (GEO) and TCGA. The results showed that the expression of TRIM56 in LUAD samples was lower than that observed in normal tissues in GEO datasets ([63]GSE117049 and [64]GSE43458) ([65]Figure 1A and [66]B). Consistent with the GEO results, TRIM56 showed low expression in LUAD tissue from the TCGA-LUAD cohort ([67]Figure 1C). To verify the mRNA expression of TRIM56 in lung adenocarcinoma cells and normal lung cells, we performed qRT-PCR to detect the TRIM56 mRNA levels. We found that the mRNA levels of TRIM56 was down-regulated in A549 cells compared to Beas-2B cells ([68]Figure 1D). In addition, representative images from The Human Protein Atlas ([69]https://www.proteinatlas.org/)[70]^25 demonstrated that TRIM56 protein expression was lower in LUAD samples compared with normal lung tissues ([71]Figure 1E). Figure 1. [72]Figure 1 [73]Open in a new tab Low expression of TRIM56 in lung adenocarcinoma is associated with poor prognosis. (A and B) Analysis of TRIM56 expression in GEO datasets. (A) [74]GSE117049 (N=11, T=57), (B) [75]GSE43458 (N=30, T=80). *P<0.05, **P < 0.01. (C) Box plot showing the RNAseq dataset, LUAD tissue (N=515, T=59 in TCGA cohort via UALCAN. **P < 0.01. (D) mRNA of TRIM56 were detected by qRT-PCR in Beas-2B cells and A549 cells. **P < 0.01. (E) The protein levels of TRIM56 in lung adenocarcinoma tissues compared with normal tissues in The Human Protein Atlas database. Normal lung tissue (patient ID: 2268; staining: high; intensity: strong; quantity: >75%), lung Adenocarcinoma sample (patient ID: 3048; staining: Medium; intensity: Moderate; quantity: 75‑25%). (F) TRIM56 expression associated Lung adenocarcinoma patient Overall survival (OS), pre-tumor survival (FPS) and late-stage tumor survival (PPS), data from [76]GSE31210. Next, we evaluated the prognostic value of TRIM56 in LUAD. We used the [77]GSE31210 dataset to analyze the relationship between the mRNA expression of TRIM56 and overall, pre-tumor, and late-stage tumor survival. The analysis results showed that the downregulation of TRIM56 was associated with poor prognosis of LUAD in terms of overall, pre-tumor, and late-stage tumor survivals ([78]Figure 1F). Therefore, TRIM56 is poorly expressed in LUAD and related to prognosis. In addition, we analyzed the relationship between the expression of TRIM56 and the clinicopathological features of LUAD. Notably, the expression levels of TRIM56 did not show a difference between groups classified by sex, smoking habit, nodal metastasis status, age, and tumor stage ([79]Supplementary Figure 1A–[80]E). The expression of TRIM56 was significantly higher in Caucasian than African-American patients ([81]Supplementary Figure 1F). In addition, its expression was higher in the TP53 mutation group versus the TP53 non-mutation group ([82]Supplementary Figure 1G). TRIM56 was downregulated in the three subtypes of LUAD, namely proximal inflammatory, proximal proliferative, and terminal respiratory unit ([83]Supplementary Figure 1H–[84]J).[85]^26^,[86]^27 The reduced expression of TRIM56 was more pronounced in the proximal proliferative subtype ([87]Supplementary Figure 1I). This subtype is associated with LUAD metastases to the brain,[88]^28 suggesting that TRIM56 may affect tumor migration and invasion. The Expression of TRIM56 Was Regulated by DNA CNA in LUAD In this study, we used RNA sequencing, copy number, mutation, and DNA methylation data available in TCGA-LUAD to investigate the potential mechanism of TRIM56 dysregulation. Firstly, we detected the genetic changes of TRIM56 in the TCGA-LUAD cohort using the cancer genomics data available in cBioportal. Among 816 LUAD samples for which sequencing and DNA CNA data were available, changes in the DNA sequence of TRIM56 were observed in 27 (4%) samples ([89]Figure 2A); the main type of gene mutation was gene amplification. Subsequently, we used the UCSC genome browser database ([90]https://genome.ucsc.edu) to simultaneously detect mRNA expression, CNA, and DNA methylation of TRIM56 in a total of 611 patients ([91]Figure 2B). Next, we used the TCGA-LUAD dataset in cBioportal to analyze the effect of CNA on the expression of TRIM56 ([92]Figure 2C). The results showed a significant correlation between DNA deletion and the downregulation of TRIM56 expression. Furthermore, the probability of a missense mutation (Missense) and a truncating mutation (Truncating) in TRIM56 was 1.86% and 0.16%, respectively. However, the mutation of TRIM56 did not affect its mRNA levels ([93]Figure 2D). Although regression analysis showed that the expression of TRIM56 was negatively correlated with its methylation status ([94]Figure 2E), the level of TRIM56 promoter methylation in LUAD samples was lower than that recorded in normal tissues ([95]Figure 2F). Therefore, we hypothesized that DNA methylation does not play a major role in the regulatory mechanism for the expression of TRIM56. Figure 2. [96]Figure 2 [97]Open in a new tab TRIM56 DNA copy number alteration (CNA), mutation and methylation in LUAD. (A) Genetic alteration of TRIM56 in 816 cases of LUAD via cBioportal. (B) Heatmap of TRIM56 mRNA expression, copy number alteration and DNA methylation in patients with primary LUAD. −2: Homozygous deletion; −1: Heterozygous loss, 0: Copy neutral; +1: Low-level copy gain; +2: High-level amplification. Data from GDC TCGA Lung Adenocarcinoma (LUAD) of UCSC (C) Plot chart showing TRIM56 expression in LUAD tissues with copy-number alteration. *P<0.05. (D) Scatter plot showing mutation alteration of TRIM56 in LUAD patients via cBioPortal. ns, no significance. (E) Regression analysis of TRIM56 DNA methylation and TRIM56 expression. (F) Boxplot showing Promoter methylation of TRIM56 from LUAD patients via UALCAN. *P<0.05. miRNAs and Transcription Factors May Regulate the Expression of TRIM56 in LUAD In addition to changes in its DNA sequence, we hypothesized that TRIM56 may be regulated by miRNAs and transcription factors. Firstly, we used LinkedOmics to identify the miRNAs that are negatively related to TRIM56 expression. The volcano map is shown in [98]Figure 3A. The red and green dots represent miRNAs that are positively (n=407) and negatively (n=403) correlated with TRIM56 expression, respectively. Typically, miRNAs suppress the expression of their target genes. In this analysis, we explored the top 50 miRNAs most negatively correlated with TRIM56 expression ([99]Figure 3B). Figure 3. [100]Figure 3 [101]Open in a new tab miRNA-TRIM56 network in LUAD. (A) Pearson test was used to analyze correlations between TRIM56 and miRNA differentially expressed in LUAD. (B) Heat maps showing miRNAs negatively correlated with TRIM56 in LUAD (TOP 50). (C–F) Enrichment analysis of the miRNAs negatively correlated with TRIM56 in LUAD. The groups display the enrichment results of the top 30 miRNA that was negatively correlated with TRIM56 in LUAD. (C) Pathway. (D) Biological processes. (E) Molecular functions. (F) Cellular components. (G and H) Kaplan-Meier Plotter was used to analyze the effect of miRNAs expression on the overall survival of LUAD patients: (G) has-mir-542 (H) has-mir-627. Next, we used the FunRich tool to perform enrichment analysis on the target genes of the top 30 miRNAs that were screened. The results showed that the pathways for the enrichment of target genes of these miRNAs include the glypican pathway, proteoglycan syndecan-mediated signaling events, interferon-gamma (IFNG) pathway, alpha9 beta1 integrin signaling events, ErbB receptor signaling network, and vascular endothelial growth factor (VEGF) and VEFG receptor (VEGFR) signaling network ([102]Figure 3C). The main enrichment biological processes of these miRNA target genes are signal transduction, cell communication, transport, regulation of nucleobase, nucleoside, nucleotide and nucleic acid metabolism, regulation of cell growth, regulation of gene expression, and epigenetic regulation ([103]Figure 3D). The main enriched molecular functions of these miRNA target genes are transcription factor activity, ubiquitin-specific protease activity, protein serine/threonine kinase activity, GTPase activity, receptor binding, and receptor signaling complex scaffold activity ([104]Figure 3E). The target genes of these miRNAs are mainly located in the nucleus, cytoplasm, lysosome, endosome, Golgi apparatus, and receptor complex ([105]Figure 3F). In addition, we analyzed the expression and prognostic value of the top 20 miRNAs, and found that hsa-mir-542 and hsa-mir-627 were abnormally overexpressed and related to poor prognosis of LUAD ([106]Figure 3G and [107]H; [108]Table 1). These results imply that these miRNAs may affect the prognosis of LUAD by regulating the expression of TRIM56. Table 1. Overall Survival of miRNA That Target TRIM56 in LUAD miRNA Low Expression High Expression Prognosis P value HR hsa-mir-450b 253 251 NA 0.310 1.16(0.87–1.56) hsa-mir-556 145 359 Favorable 0.030 0.71(0.52–0.97) hsa-mir-340 317 187 NA 0.270 0.84(0.63–1.14) hsa-mir-454 215 289 NA 0.250 0.84(0.63–1.13) hsa-mir-96 181 323 NA 0.095 0.78(0.58–1.05) hsa-mir-592 245 259 NA 0.240 0.84(0.63–1.12) hsa-mir-362 179 325 NA 0.092 0.77(0.57–1.04) hsa-mir-542 360 144 Unfavorable 0.012 1.47(1.08–1.99) hsa-mir-3613 356 148 Favorable 0.046 0.72(0.52–1) hsa-mir-130b 247 257 NA 0.260 1.18(0.88–1.59) hsa-mir-424 244 260 NA 0.220 1.2(0.89–1.62) hsa-mir-3174 175 329 NA 0.180 0.81(0.6–1.1) hsa-mir-98 200 304 NA 0.190 0.82(0.61–1.11) hsa-mir-450a 364 140 NA 0.073 1.33(0.97–1.82) hsa-mir-301b 304 200 NA 0.068 0.75(0.56–1.02) hsa-mir-627 128 376 Unfavorable 0.027 1.53(1.05–2.25) hsa-mir-548u 154 350 NA 0.340 0.86(0.64–1.17) hsa-mir-938 316 188 NA 0.300 1.17(0.87–1.58) hsa-mir-18a 266 238 NA 0.052 1.33(1–1.79) hsa-mir-374a 242 262 Favorable 0.0039 0.62(0.49–0.87) hsa-mir-708 285 219 NA 0.086 0.77(0.57–1.04) hsa-mir-425 233 271 NA 0.160 0.81(0.61–1.09) hsa-mir-33a 358 146 NA 0.280 0.84(0.6–1.16) hsa-mir-195 203 301 Favorable 0.014 0.69(0.52–0.93) hsa-mir-1256 308 196 NA 0.260 1.19(0.88–1.6) hsa-mir-106a 281 223 NA 0.053 1.33(1–1.78) hsa-mir-19a 171 333 NA 0.190 0.82(0.61–1.1) hsa-mir-30b 195 309 Favorable 0.0003 0.58(0.44–0.78) hsa-mir-216b 371 133 NA 0.23 1.21(0.89–1.64) hsa-mir-548y 280 224 NA 0.11 1.26(0.95–1.69) [109]Open in a new tab The regulation of TRIM56 expression by transcription factors has not been previously reported. Thus, we analyzed the potential transcription factors of TRIM56. Firstly, we retrieved the TRIM56 promoter sequence through the UCSC Genome Browser, and subsequently used PROMO and ConTra V3 to identify possible transcription factors ([110]Table 2). There were 22 and four potential transcription factors identified from PROMO and ConTra v3, respectively. We also examined the correlation between the expression of these transcription factors and that of TRIM56 in LUAD. TRIM56 was significantly correlated with the following transcription factors: transcription factor II–I (TFII-I), RXR-alpha, AP-2 alpha, Pax-5, c-Ets-1, zinc finger and BTB domain containing 7B (ZBTB7B), and SP2 ([111]Table 2). Therefore, these transcription factors may regulate the expression of TRIM56 in LUAD. Table 2. Transcription Factor of TRIM56 in PROMO and ConTra V3 Database Gene Correlation P value PROMO YY1 −0.09 0.043 C/EBPbeta −0.19 1.106e-05 IRF-2 0.19 1.111e-05 ER-alpha −0.02 0.716 STAT4 0.01 0.762 TFII-I^# 0.60 3.763e-52 HNF-3alpha −4.87e-03 0.912 GR 0.18 2.633e-05 RXR-alpha^# 0.34 3.352e-15 HOXD9 −0.12 8.670e-03 HOXD10 −0.09 0.050 AP-2alphaA^# 0.31 4.134e-13 Pax-5^# 0.21 1.578e-06 GATA-1 0.08 0.088 TFIID −0.11 0.014 IRF-1 0.18 3.755e-05 Elk-1 0.06 0.188 c-Ets-1^# 0.34 5.727e-15 ConTra V3 ZBTB7B^# 0.24 7.973e-08 SP4 0.17 1.298e-04 SP2^# 0.49 3.318e-32 KLF15 −0.04 0.330 [112]Open in a new tab Note: ^#Indicate high correlation with TRIM56 expression, correlation> 0.5. Gene Ontology (GO) Functional Annotation and Pathway Enrichment of Proteins Interacting with TRIM56 As an E3 ubiquitin ligase, TRIM56 functions by binding to other proteins. Proteome analysis revealed 20 proteins interacting with TRIM56; these proteins were displayed in a physical interaction network using STRING and Cytoscape ([113]Figure 4A). The most reliable proteins were LSM2, TRIM38, Mab-21 domain containing 1 (MB21D1), T cell receptor associated transmembrane adaptor 1 (TRAT1), transmembrane protein 137 (TMEM137), and TRIM32 ([114]Figure 4A). Figure 4. [115]Figure 4 [116]Open in a new tab GO Functional Annotation and Pathway Enrichment of TRIM56 interaction proteins (A) Protein-protein interaction network of TRIM56 via STRING. (B–E) Enrichment analysis of the protein interacted with TRIM56. (B) Biological processes (C) Molecular functions. (D) Cellular components. (E) KEGG Pathway. Next, we performed GO and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis of the 20 aforementioned proteins interacting with TRIM56. These proteins are involved in numerous biological processes, such as negative regulation of gene expression, nuclear-transcribed mRNA catabolic process, nonsense-mediated decay, protein targeting, signal recognition particle-dependent cotranslational protein targeting to membrane, mRNA metabolic process, GO: 0044265cellular macromolecule catabolic process, and nuclear-transcribed mRNA catabolic process ([117]Figure 4B). Moreover, we noticed significant enrichment of functions corrected with U6 snRNA binding, structural constituent of ribosome, RNA binding, translation initiation factor binding, and rRNA binding ([118]Figure 4C). These proteins are mainly localized to the U6 snRNP, ribonucleoprotein complex, spliceosomal tri-snRNP complex, U4/U6 x U5 tri-snRNP complex, Lsm1-7-Pat1 complex, and ribosome ([119]Figure 4D). The Kyoto Encyclopedia of Genes and Genomes pathways in which these proteins are mainly involved are (hsa03010) ribosome, (hsa03018) RNA degradation, (hsa03040) spliceosome, and (hsa04623) cytosolic DNA-sensing pathway ([120]Figure 4E). The predicted ubiquitination substrates of TRIM56 are shown in [121]Table 3. The expression of TP53, H2A histone family member X (H2AFX), lymphoid enhancer binding factor 1 (LEF1), and hepatocyte growth factor-regulated tyrosine kinase substrate (HGS) is abnormal in LUAD. In summary, we analyzed the proteins interacting with TRIM56 and predicted the potential ubiquitination substrates of TRIM56. Table 3. Ubiquitination Substrates of TRIM56 in LUAD from UbiBrowser Genes Score Gene Expression Genes Score Gene Expression Log2(TPM + 1) for Log-Scale Log2(TPM + 1) for Log-Scale Normal Tumor Normal Tumor TP53^# 0.746 4.5 5.0 SOS1 0.645 3.4 3.3 DLG2 0.645 1.2 0.9 IL6R 0.645 4.6 3.5 H2AFX^# 0.645 4.0 5.0 TRIM3 0.645 3.3 3.0 TJP2 0.645 5.6 4.6 RYR2 0.645 1.3 0.5 MPP5 0.645 4.0 3.3 PML 0.645 5.4 5.3 TJP1 0.645 5.9 4.8 SHANK2 0.645 3.2 2.5 SHANK3 0.645 3.5 2.1 CRK 0.645 5.2 4.6 LEF1^# 0.645 2.5 3.0 TRIM8 0.645 6.2 6.2 SH3BP4 0.645 4.2 4.2 HGS^# 0.645 5.8 6.4 TRIM69 0.645 4.4 3.8 MPP7 0.645 3.5 3.6 [122]Open in a new tab Overexpression of TRIM56 Suppressed the Invasion and Migration of LUAD Cells To determine the effect of TRIM56 on the invasive and migratory ability of LUAD cells, we stably expressed Myc-TRIM56 in H1299 cells via a lentiviral system. As shown in [123]Figure 5A, the expression of TRIM56 was detected by Western blotting. Next, we examined the invasive and migratory ability of pCDH-TRIM56 H1299 cells using Transwell assays ([124]Figure 5B). We found that the invasion and migration of H1299 cells overexpressing TRIM56 were significantly reduced than those observed in the control group. Consistent with these results, wound healing experiments showed that the migratory ability of pCDH-TRIM56 H1299 cells was significantly downregulated ([125]Figure 5C). Therefore, overexpression of TRIM56 may suppress the invasion and migration of LUAD cells. Figure 5. [126]Figure 5 [127]Open in a new tab Overexpression TRIM56 could inhibits the invasion and migration of lung adenocarcinoma cells H1299 (A) Western blot confirmed TRIM56 overexpression effect. GAPDH was used as an internal loading control. (B) Cell invasion and migration were examined by transwell assay (scale bar, 200 μm). H1299 cell were cultured in the transwell culture chambers, after incubated 48 h, used crystal violet to stain cells attached to the lower surface of the filter and counted. The quantitative analysis are shown in right panels. The data are presented as the mean ± s.e.m. from three independent experiments. (n = 3). **P < 0.01, ***P < 0.001. (C) Wound healing assay (scale bar, 200 μm) for investigating migration potential of TRIM56 overexpression in H1299 cells. Right histogram represents quantification analysis. The data are presented as the mean ± s.e.m. from three independent experiments. (n = 3). *P < 0.05. Discussion TRIM56 catalyzes the transfer of ubiquitin molecules to specific substrates.[128]^4 TRIM56 was localized in the cytoplasm and failed to interact with double-stranded DNA. However, it is responsible for the induction of type I IFN by promoting the function of STING through ubiquitination.[129]^7 In addition, TRIM56 can activate TLR3 antiviral signaling pathway-induced interferon.[130]^9 These findings confirm that TRIM56 plays an important role in antiviral and immune regulation. Previous studies have shown that the expression of TRIM56 is reduced in a variety of human cancers, including multiple myeloma[131]^11 and ovarian cancer.[132]^29 In these tumors, the low expression of TRIM56 is associated with tumor development, metastasis, and poor prognosis. However, TRIM56 can promote breast cancer through ubiquitination and stabilize ER-α.[133]^13 In summary, TRIM56 plays various roles in different tumors. This is the first study to demonstrate that the expression of TRIM56 in LUAD was decreased and related to poor prognosis. The migration and invasion experiments in vitro were verified that TRIM56 overexpression suppress the invasion and migration of lung adenocarcinoma cells. It provides a basis for further in vivo experiments. In general, genetic and epigenetic alterations can regulate the tumor phenotype and patient survival.[134]^30 In this study, the potential mechanism underlying the dysregulation of TRIM56 was further investigated using deep sequencing data from the TCGA-LUAD cohort. The present findings showed that upregulated TRIM56 expression was significantly associated with DNA copy deletion rather than DNA mutations and methylation. Moreover, TRIM56 expression was negatively correlated with DNA methylation. Nevertheless, the methylation of the promoter of TRIM56 in patients with LUAD was reduced, implying the presence of an unknown mechanism for the regulation of TRIM56 methylation in LUAD. Previous research found that hsa-mir-9 promoted the development and progression of multiple myeloma by regulating the TRIM56/NF-κB pathway.[135]^31 In this study, we found two miRNAs negatively correlated with TRIM56 expression and associated with poor prognosis of LUAD, namely hsa-mir-542 and hsa-mir-627. These miRNAs exert important regulatory effects on tumors. Hsa-mir-542 can be used as a non-invasive blood biomarker for tumor monitoring and prognosis prediction in patients with osteosarcoma.[136]^32 Moreover, it can target CDK14 to inhibit the development of ovarian cancer,[137]^33 and inhibit the growth and invasion of colon cancer cells through PI3K/AKT/survivin signaling.[138]^34 Hence, these two miRNAs may affect the prognosis of LUAD by inhibiting the expression of TRIM56. The enrichment analysis of proteins interacting with TRIM56 showed that TRIM56 may be involved in numerous biological processes of mRNA splicing and ribosomes. During mRNA splicing, pre-mRNA is assembled into a large, dynamic spliceosome. This process requires the participation of five small nuclear RNAs (U1, U2, U4, U5, and U6) and more than 50 protein components.[139]^35 The splicing and maturation of mRNA are essential for cell growth and development. Therefore, it is possible that TRIM56 regulates tumors by regulating gene splicing and maturation. Conclusion In this study, we analyzed the expression of TRIM56 in LUAD using bioinformatics. The results demonstrated that the low expression of TRIM56 was related to poor prognosis of LUAD. In addition, we discussed the change in gene copy number and the regulatory effects of miRNA and transcription factors on the mRNA expression levels of TRIM56. Furthermore, we predicted the proteins interacting with TRIM56 and analyzed their enrichment. Finally, our experimental results showed that TRIM56 can inhibit the invasion and migration of LUAD. In summary, the present research showed, for the first time, that TRIM56 is a potential tumor suppressor. This evidence provides a new direction for analyzing the reasons for the low expression of TRIM56 in LUAD and its regulatory mechanism involved in tumor migration and invasion. Acknowledgments