- 1Shandong Key Laboratory of Radiation Oncology, Shandong Cancer Hospital and Institute, Shandong First Medical University and Shandong Academy of Medical Sciences, Jinan, China
- 2Medical Science and Technology Innovation Center, Shandong First Medical University, Jinan, China
- 3Key Laboratory of Cancer Proteomics of Chinese Ministry of Health, Xiangya Hospital, Central South University, Changsha, China
Abnormal ubiquitination is extensively associated with cancers. To investigate human lung cancer ubiquitination and its potential functions, quantitative ubiquitinomics was carried out between human lung squamous cell carcinoma (LSCC) and control tissues, which characterized a total of 627 ubiquitin-modified proteins (UPs) and 1209 ubiquitinated lysine sites. Those UPs were mainly involved in cell adhesion, signal transduction, and regulations of ribosome complex and proteasome complex. Thirty three UPs whose genes were also found in TCGA database were significantly related to overall survival of LSCC. Six significant networks and 234 hub molecules were obtained from the protein-protein interaction (PPI) analysis of those 627 UPs. KEGG pathway analysis of those UPs revealed 47 statistically significant pathways, and most of which were tumor-associated pathways such as mTOR, HIF-1, PI3K-Akt, and Ras signaling pathways, and intracellular protein turnover-related pathways such as ribosome complex, ubiquitin-mediated proteolysis, ER protein processing, and proteasome complex pathways. Further, the relationship analysis of ubiquitination and differentially expressed proteins shows that ubiquitination regulates two aspects of protein turnover - synthesis and degradation. This study provided the first profile of UPs and molecular networks in LSCC tissue, which is the important resource to insight into new mechanisms, and to identify new biomarkers and therapeutic targets/drugs to treat LSCC.
Introduction
Lung squamous cell carcinoma (LSCC) is closely related to smoking (1), which is 25%-30% of non-small cell lung cancer (NSCLC), and >70% LSCC patients were diagnosed in late stage (2, 3), with high mortality (~0.4 million/year) in the world (4). Although the incidence has decreased due to tobacco control and changes in lifestyle in recent decades, LSCC still brings huge burdens to society (5). In terms of treatment, radical thoracic surgery is preferred for early-stage cases of LSCC; however, unfortunately the proportion of this type of patients is very small. Besides, chemotherapy and radiotherapy are also indispensable treatments, especially in advanced (or metastatic) LSCC subpopulation (6). Recent years, although some encouraging targeted drugs such as anaplastic lymphoma kinase inhibitors and epithelial growth factor receptor (EGFR) tyrosine kinase inhibitors have been used to significantly improve the therapy of late-stage lung adenocarcinoma (4), there is still lack of FDA-approved targeted therapy for advanced LSCC patients. For most patients with advanced LSCC, first-line standard chemotherapy is a 4-6-cycle platinum-based third-generation cytotoxic drug with an effective rate of 20%-30%. The median survival time (MST) is 8-10 months, the time to progression (TTP) is 3-5 months (7). In addition, LSCC still lacks universally accepted biomarkers for early diagnosis and prognostic evaluation. Obviously, the current predictive, preventive, and personalized medicine (PPPM) of LSCC is still unmet; therefore, it is necessary to clarify molecular mechanism and discover new biomarkers and therapeutic targets/drugs to manage LSCC.
Protein ubiquitination is an important molecular event, which covalently link the C-terminus of ubiquitin (8.5 kDa, 76 amino acids, coded by four genes RPS27A, UBA52, UBC, and UBB) to the ϵ-amino group at lysine residue in a protein (8, 9). Ubiquitination is a multi-step reaction process, which is catalyzed by enzymes E1, E2, and E3 (8). Protein ubiquitinations are classified into monoubiquitination that means one ubiquitin in one protein, multiubiquitination that means several ubiquitination sites in a protein and each ubiquitination site links one ubiquitin, and polyubiquitination that means a polyubiquitin chain to be added to lysine residues or N-terminus of the previous ubiquitin (10); and different types of ubiquitinations have different biological functions (11). Also, there are seven lysine residues in an ubiquitin, which make polyubiquitin chain more complicated. Moreover, ubiquitination process can be reversed, which is catalyzed by more than 100 deubiquitination enzymes (12). The reversible ubiquitination reaction regulates multiple biological processes including protein degradation (13), and associated with a wide range of diseases including cancers, and inflammatory diseases (14). Ubiquitin-proteasome system (UPS) is the important protein degradation pathway, which has important roles in tumorigenesis. Targeting the UPS system is the promising anti-tumor strategies. Thereby, several effective anti-tumor drugs targeting UPS have been approved by FDA, including bortezomib (15), carfilzomib (16), thalidomide, pomalidomide, and lenalidomide (17), for different cancer treatment. At present, LSCC still lacks specific early diagnostic biomarkers, and in addition to surgery, radiotherapy, and chemotherapy, there is lack of effective molecular therapeutic drugs. Therefore, the study of ubiquitinome might have important scientific merits for discovery of effective biomarkers and novel therapeutic targets/drugs for LSCC.
Liquid chromatography in combination with tandem mass spectrometry (LC-MS/MS) is an effective method to identify ubiquitination sites and quantify the abundance of ubiquitination (18). However, endogenous ubiquitination is low abundant event, it is necessary to enrich the ubiquitinated tryptic peptides prior to MS/MS analysis. The commercial anti-K-ϵ-GG antibodies is able to effectively enrich ubiquitinated tryptic peptides (19). For antibody-based enrichment of ubiquitinated tryptic peptides, it is the best to avoid any salt, acid, and basic factors to negatively impact this enrichment. Compared to the isotope-labeling quantitative proteomics such as iTRAQ (isobaric tags for relative and absolute quantification) and TMT (tandem mass tags) that introduce above interfered factors to negatively impact antibody-antigen reaction, label-free quantitative proteomics is a preferred approach to identify and quantify endogenous ubiquitination, which does not label the analyzed protein component in a proteome (20) but quantify the protein with area under the curve (AUC) and signal intensity in the MS spectrum, and spectral counting based on MS/MS analysis (21). Thus, anti-ubiquitin antibody in combination with label-free LC-MS/MS is an effective quantitative ubiquitinomics, which has identified up to 10,000 ubiquitination sites in an experiment (22). For lung cancer analysis, currently researchers mainly focused on the study on ubiquitinome of lung cancer cells (23, 24), and these studies were not based strictly on LSCC cell lines.
This study used anti-ubiquitin antibody coupled with label-free LC-MS/MS to detect, identify, and quantify ubiquitinated proteins (UPs) and ubiquitination sites in fresh human LSCC tissues compare to controls. Gene ontology (GO) enrichment analysis was used to determine the functional characteristics of UPs, and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was used to mine the ubiquitination-involved pathway networks. This study provides the holistic profile of ubiquitination in LSCC tissue, and underscores ubiquitination may promote the occurrence and development of LSCC by affecting the intracellular protein turnover, and paving the way for further research on specific ubiquitination regulatory mechanisms.
Materials and methods
Tissue samples and protein preparation
Human LSCC tissues (n =5) and tumor-adjacent control (TAC) tissues (n=5) (Supplementary Table 1; these samples were at early stage of LSCC) were removed by thoracic surgery, and stored in –80°C, for protein extraction. A total of 750 mg LSCC tissues (n = 5; 150 mg per patient) and 750 mg TAC tissues (n=5; 150 mg per patient) derived from 5 patients were homogenized, after removal of blood contamination, in solution that included 2 M thiourea, 7 M urea, 1 mM phenylmethyl-sulfonyl fluoride (PMSF), and 100 mM dithiothreitol (DTT); followed by sonication (80 W, 10 s, and interval 15 s; 10 x), and centrifugation at 4°C (15,000xg, and 20 min). The supernatant was the extracted proteins, whose protein content was determined with Brandford assay.
Preparation of tryptic peptides and enrichment of ubiquitinated tryptic peptides
The protein sample (LSCC; TAC) was reduced with final concentration of 10 mM DTT (37 °C, 1.5 h), and alkylated with final concentration of 50 mM iodoacetamide (30 min, dark). The mixture was diluted to 2 M uranyl acetate with 50 mM Tris-HCl (pH 8.0), follwed by addition of trypsin (1:50 for the ratio of trypsin to protein), and incubation (15-18 h; 37 °C) to digest the proteins. Then, the digestion reaction was stopped with addition of trifluoroacetic acid (TFA) to adjust pH ≤ 3. The tryptic peptide mixture was desalted, and lyophilized. The lyophilized tryptic peptides were redisolved in 1.4 mL solution that contained 50 mM NaCl, 10 mM Na2HPO4, and 50 mM MOPS/NaOH pH 7.2, followed by addition of anti-K-ϵ-GG antibody beads [PTMScan ubiquitin remnant motif (K-ϵ-GG) kit, Cell Signal Technology], incubation (1.5 h, 4 °C), and centrifugation (30 s, 2,000×g). The beads with antibody-binding ubiquitinated peptides were washed for 3 times with 1 mL solution that contained 50 mM NaCl, 10 mM Na2HPO4, and 50 mM MOPS/NaOH pH 7.2, followed by washing with water for 3 times. The ubiquitinated tryptic peptides were eluted with 40 μL 0.15% TFA, followed by centrifugation (30 s, 2,000×g), and desaltation with C18 STAGE Tips.
LC-MS/MS analysis of enriched ubiquitinated peptides
The prepared ubiquitinated tryptic peptides were treated with a reverse-phase trap column (nanoViper C18, 2 cm x 100 μm), and then were separated in a reverse-phase analytical column (C18, 10 cm length x 75 μm i.d., and 3 μm resin) with buffers A (0.1% formic acid) and B (0.1% formic acid + 84% acetonitrile) for 2 h (flow-rate =300 nl/min) on the Q Exactive mass spectrometer (Thermo Scientific) to obtain MS/MS raw data. MaxQuant software (version 1.5.3.17) was used to search MS/MS data against human protein database uniprot_human_156639_20170105.fasta, with database search parameters (4 missed cleavages, trypsin, 6 ppm for MS, 20 ppm for MS/MS, GlyGly at Lys, acetylation at protein N-term, oxidation at Met, and carbamidomethyl at Cys). The peptide was determined by MS/MS data with false discovery rate (FDR) ≤0.01, ubiquitination site with FDR ≤0.01, and protein with FDR ≤0.01. The abundance of ubiquitination was determined with MaxQuant label-free calculation based on peptide intensity and peptide counting.
Bioinformatics and statistical analysis
The functional characteristics of the identified UPs, including biological processes (BPs), molecular functions (MFs), and cellular components (CCs) were analyzed with DAVID software (version 6.8, https://david.ncifcrf.gov/) (25). Ubiquitination-mediated signaling pathways were analyzed with the KEGG tool KOBAS (http://kobas.cbi.pku.cn) (26). The protein-protein interaction (PPI) networks that UPs are involved in were analyzed with STRING software (https://string-db.org/; version 10.0) (27), with statistical significance of a combined score > 0.4, in the Cytoscape program (version 3.6.0) (28). Also, molecular complex detection (MCODE) in the Cytoscape program was used to identify significant modules in the PPI networks with cutoff values (max depth=100, node score=0.2, degree=2, k-score=2, and MCODE score>5) (29). The hub nodes of PPI network were determined with cytoHubba (version 0.1) in Cytoscape package through 6 different topological algorithms (30), including Maximal Clique Centrality (MCC; cutoff score=100 with this node that had at least 5 molecules linked), Maximum Neighborhood Component (MNC; cutoff score=10 with this node that had at least 5 molecules linked), Degree (cutoff score=10 with this node that had at least 5 molecules linked), Edge Percolated Component (EPC; cutoff score=10 with this node that had at least 5 molecules linked), Closeness (cutoff=150 with this node that had at least 5 molecules linked), and Raliality (cutoff score=5 with this node that had at least 5 molecules linked). The final hub nodes were identified with overlapped analysis of six sets of hub nodes that were derived from 6 topological algorithms. The transcriptomic data (level 2 count data) and corresponding clinical survival data of LSCC patients were obtained from TCGA database with R-languages “RTCGA.mRNA” (version 1.12.0) and “RTCGA.clinical” (version 20151101.14.0), respectively. The survival analysis was performed with R-languages “survival” (version 2.44-1.1) and “survminer” (version 0.4.6).
The overall flow-chart to identify and analyze the ubiquitination profile and the corresponding functional characteristics in LSCC was summarized (Figure 1).
Results
Ubiquitination profile in human LSCC
In total, 627 UPs and 1209 ubiquitination sites were identified with MS/MS data, which referred to 1086 ubiquitinated tryptic peptides, in human LSCC tissues (Supplementary Table 2). For example, the MS/MS spectrum of tryptic peptide 55TLSDYNIQK*ESTLHLVLR72 ([M+2H]2+ m/z =1122.6031, RT = 140.86 min) was derived from 40S ribosomal protein S27a (RPS27A; P62987) (Figure 2A), with high S/N (signal to noise) ratio, excellent y-ion (y1-y17) and b-ion (b2-b10, b13-b16) series, localization of ubiquitination site at residue K63, and a decreased ubiquitination level at residue K63 (Ratio T/N = 0.425; P=0.0529) in LSCCs compared to TAC controls (Supplementary Table 2). In addition, another six ubiquitination sites were also identified in this RPS27A (P62987) in LSCC tissues compare to TAC tissues, including ubiquitinations at residues K6 (Ratio T/N = 0.271; P=5.07E-03), K11 (Ratio T/N=0.496; P=1.28E-02), K27 (Ratio T/N=3.812; P=4.53E-05), K29 (Ratio T/N=0.933; P=8.55E-01), K33 (Ratio T/N=2.903; P=3.75E-04), and K48 (Ratio T/N=0.839; P=5.47E-01) (Supplementary Table 2). It clearly demonstrated that 7 ubiquitination sites were identified in this RPS27A, with decreased ubiquitination levels at residues K6, K11, K29, K48, and K63 (of them, significantly decreased level at residues K6 and K11), and significantly increased ubiquitination levels at residues K27 and K33. This is an interesting result, because, in the same protein, the ubiquitination level at some Lys sites were decreased (K6, K11, K29, K48, and K63), and whereas the ubiquitination levels at some Lys sites were increased (K27 and K33). This phenomenon might be easily explained from the angle of proteoforms encoded by the same gene RPS27A. Another reason might be due to the kinetics change of ubiquitin linkages in LSCC because there are four ubiquitin-encoding genes, including UBB, UBC, UBA52 (which encodes RPL40), and UBA80 (also known as RPS27A; RPS27A encodes a fusion protein consisting of ubiquitin at the N terminus and ribosomal protein S27a at the C terminus), and all of which produce the identical ubiquitin protein; if S27a protein does not change, but the ubiquitin chain linking pattern increases or decreases, it is not due to the RPS27A gene alone; thereby thus results in the alteration of ubiquitin ligation dynamics mediated by ubiquitin ligases in LSCC.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970843/fendo-13-970843-HTML-r1/image_m/fendo-13-970843-g002.jpg)
Figure 2 Representative MS/MS spectra and holistic characteristics of ubiquitination sites in LSCC. (A) The identified ubiquitinated peptide 55TLSDYNIQK*ESTLHLVLR72 from 40S ribosomal protein S27a (RPS27A, P62987). K* = ubiquitinated lysine residue. (B) The length distribution of identified ubiquitinated peptides. (C) Pie-chart demonstrating the number of ubiquitination sites per protein. (D) Volcano plot demonstrated the changes in the degree of ubiquitination at a special lysine site among 642 ubiquitination sites (T+/N+) according to the fold change and p value between the LSCC and control groups. X axis: the log of the fold change between the two conditions (log base 2). Y axis: the negative log of the p value (log base 10). Red dot indicates that the degree of ubiquitination of the site is significantly increased (Ratio T/N>2.0; p<0.05) and blue indicates significantly decreased (Ratio T/N < 0.5; p<0.05).
For 1086 ubiquitinated tryptic peptides, their lengths were in the range of 5-40 amino acids (Figure 2B). Analysis of the number of ubiquitination sites in human LSCC UPs found that the number of putative ubiquitination sites per protein was arranged from 1 to 15. A total of 388 (388/627 = 62%) UPs had only one ubiquitination site, 18% (111/627) UPs had two sites, 10% (63/627) UPs had three sites, 4% (23/627) UPs had four sites, and 7% (42/627) UPs had five sites (Figure 2C). A total of 1133 (1133/1209 = 93.7%) ubiquitination sites had quantitative information, including 642 (642/1133 = 56.7%) sites with T+/N+, 395 (395/1133 = 34.9%) with T+/N-, and 96 (96/1133 = 8.4%) with T-/N+; and among these 1133 sites, 699 (61.7%) sites had increased ubiquitination levels and 434 (38.3%) sites had decreased ubiquitination levels in LSCC compared to TAC tissues. Another 76 (76/1209 = 6.3%) sites did not have quantitative information between LSCC and TAC tissues. Here showed the distribution of intensity changes of these 642 ubiquitination sites (T+/N+), including 146 sites had significantly increased ubiquitination levels (Ratio T/N >2; p<0.05), and 152 sites had significantly decreased ubiquitination levels (Ratio T/N <0.5; p<0.05) (Figure 2D), the corresponding detailed information was listed in Supplementary Table 2.
Functional characteristics of indentified UPs in LSCC
The functional characteristics of 627 UPs were enriched with GO analysis, including BPs, MFs, and CCs. For BP analysis, UPs were mainly enriched in cellular process (258 UPs), metabolic process (147 UPs), biological regulation (100 UPs), and localization (85 UPs) (Figure 3A). For MF analysis, UPs were mainly enriched in binding activity (222 UPs), catalytic activity (180 UPs), transporter activity (44 UPs), structural molecule activity (38 UPs), molecular function regulator (26 UPs), and transcription regulator activity (16 UPs) (Figure 3B). For CC analysis, UPs were mainly enriched in cell (258 UPs), organelle (134 UPs), protein-containing complex (50 UPs), membrane (38 UPs), extracellular region (17 UPs), and cell junction (13 UPs) (Figure 3C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970843/fendo-13-970843-HTML-r1/image_m/fendo-13-970843-g003.jpg)
Figure 3 Functional characteristics of 627 Ups identified with GO enrichment analysis. (A) The BP profile of UPs. (B) The MF profile of UPs. (C) The CC profile of UPs.
These BPs, MFs, and CCs were further grouped into 7 functional clusters (Supplementary Table 3), including cell adhesion, multiple signaling pathways that were closely related to tumorigenesis, such as Wnt signaling, tumor necrosis factor-mediated signaling, NF-kappa B signaling, assembly and functional regulation of proteasome, assembly and functional regulation of ribosome, haptoglobin hemoglobin complex, DNA damage recognition, and hemoglobin complex.
Survival-related UPs in LSCC
A total of 1030 overall survival-related genes were obtained from LSCC data in TCGA database (Supplementary Tables 4-6). A total of 33 overlapped molecules (UPs; survival-related genes) were obtained when overlapping analysis was performed between the identified 627 UPs and 1030 overall survival-related genes (Figure 4A), which showed that those 33 overlapped molecules were survival-related UPs. GO enrichment analysis of these 33 survival-related UPs showed that these UPs were significantly involved in different BPs, and were closely related to tumorigenesis, including cell growth, blood vessel development, and regulation of intracellular signal transduction (Figure 4B). Among them, the higher expressions of 11 overlapped molecules (UPs; survival-related genes), including SPG20, PFN2, SET, CBX1, HISTIH1B, DSC3, ENO1, NONO, HISTIH4H, ACTA2 and PGK showed a better prognosis, which suggested their protective effects; while the higher expressions of 22 overlapped molecules (UPs; survival-related genes), including PRKCKBP, DDX5, ATP11A, EPB2, RPS2, ITGB1, TANK, GPRC5A, TRIM47, MVP, RNF213, SAP18, HM13, ICAM1, SLC16A4, LRSAM1, TM9SF3, CALD1, SYNPO, CLTB, CA1 and GJA5 demonstrated a worse prognosis, which suggested their tumor-promoting effects. Here showed the overall survival curve of RPS2 (Figure 4C), and RPS2 was further identified as a hub molecule.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970843/fendo-13-970843-HTML-r1/image_m/fendo-13-970843-g004.jpg)
Figure 4 Overall survival-related UPs identified with overlapping analysis between 627 UPs and survival-related genes derived from TCGA database. (A) Venn diagram showed 33 molecules were associated with prognosis. (B) The biological processes involved in these prognostic related molecules. (C) Overall survival curve of RPS2.
Ubiquitination-mediated molecular networks and network-based hub molecules derived from 627 UPs in LSCCs
To deeper understand cellular processes regulated by ubiquitination in LSCC, 627 UPs were input into STRING online software, 457/627 UPs were legible for PPI network analysis, which identified 6 significantly molecular networks (Figures 5A–F and Supplementary Table 7). The functional analyses of genes within each network were carried out with DAVID, which found that network 1 (31 nodes and MCODE score=27.55) was mainly involved in multiple proteasome and ribosomal subunits to regulate UPS and ribosome large subunits assembly, translation, protein processing, and intracellular signaling pathway; network 2 (27 nodes and MCODE score=8.08) was mainly involved in intranuclear regulation of gene expression; network 3 (40 nodes and MCODE score=7.44) was mainly referred to cell movement; network 4 (43 nodes and MCODE score=5.24) mainly participated in cell adhesion and membrane organization; network 5 (16 nodes and MCODE score=5.20) did not have significant biological processes; and network 6 (13 nodes and MCODE score=5.17) played important roles in protein transport and small GTPase-mediated signal transduction (Supplementary Table 8).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970843/fendo-13-970843-HTML-r1/image_m/fendo-13-970843-g005.jpg)
Figure 5 Molecular interaction modules and hub nodes based on PPI molecular networks of 627 UPs. (A–F) Six molecular interaction modules that meet screening criteria derived from the entire UPs PPI network. Red nodes represented the hub molecules. Blue node means it was not recognized as hub molecule due to low score.
Moreover, 627 UPs were input into cytoHubba to determine hub molecules with 6 different topological algorithms, respectively; which identified 263 hub molecules with MCC, 232 with MNC, 268 with Degree, 278 with EPC, 347 with Closeness, and 334 with Raliality (Supplementary Table 9). Overlapped analysis of these six sets of hub molecules finally obtained 234 hub molecules (Supplementary Table 9). These finally determined 234 hub molecules were clustered into 5 functional categories, including cell-cell adhesion, ribosome-related biological process, proteasome-related biological process, cellular signaling pathways, and cellular energy metabolism (Supplementary Table 10). Interestingly, these results were also consistent with those results of GO enrichment analysis of 627 UPs (Supplementary Table 3).
Ubiquitination-mediated signaling pathways in LSCC
To reveal ubiquitination-mediated signaling pathway changes in LSCCs, 47 statistically significant pathways were revealed with KEGG pathway analysis of 627 UPs (Figure 6 and Supplementary Table 11). Multiple signaling pathways were cancer-related pathways, including (i) ubiquitin-mediated proteolysis (10 UPs, and FDR =0.0184), ribosome complex (14 UPs, and FDR<0.01), ER-mediated protein processing (13 UPs, and FDR<0.01), proteasome complex (12 UPs, and FDR<0.01), and biosynthesis of amino acids (10UPs, and FDR<0.01), which was intracellular protein turnover-related pathways, and might disturb the synthesis-degradation balance of protein in tumorigenesis; (ii) PI3K-Akt signaling pathway (25 UPs, and FDR < 0.01), Ras signaling pathway (15 UPs, and FDR < 0.01), mTOR signaling pathway (10 UPs, and FDR = 0.03), HIF-1 signaling pathway (10 UPs, and FDR < 0.01), cell cycle (9 UPs, and FDR = 0.03), and apoptosis pathway (10 UPs, and FDR = 0.02), which might function in tumorigenesis; (iii) glycolysis/gluconeogenesis (10 UPs, and FDR<0.01), carbon metabolism (12 UPs, and FDR<0.01), central carbon metabolism in cancer (9 UPs, and FDR<0.01), Fructose and mannose metabolism (5 UPs, and FDR=0.019), metabolomic pathways (46 UPs, and FDR=0.031), drug metabolism-other enzymes (5 UPs, and FDR=0.0457), and sphingolipid signaling pathway (11 UPs, and FDR<0.01), which were obviously metabolism-related pathways; (iv) Tight junction (16 UPs, and FDR<0.01), adherens junction (8 UPs, and FDR=0.0126), gap junction (8 UPs, and FDR<0.01), and non-homologous end-joining (3 UPs, and FDR<0.0356), which were involved in cell adherens and junction; (v) endocytosis (18 UPs, and FDR<0.01), and SNARE interactions in vesicular transport (5 UPs, and FDR=0.0118). which might be substance transport-related pathways; (vi) dopaminergic synapse (11 UPs, and FDR=0.0115), and GABAergic synapse (7 UPs, and FDR=0.0415), which function in synapse pathway; and (vii) others pathways such as chemokine signaling pathway (12 UPs, and FDR=0.0287), Rap1 signaling pathway (14 UPs, and FDR=0.0124), proteoglycans in cancer (12 UPs, and FDR=0.0466), pathways in cancer (20 UPs, and FDR=0.0287), microRNAs in cancer (11 UPs, and FDR=0.0208), leukocyte transendothelial migration (16 UPs, and FDR<0.01), regulation of actin cytoskeleton (18 UPs, and FDR<0.01), and platelet activation (10 UPs, and FDR=0.0129). These ubiquitination-mediated cancer-related signaling pathways clearly demonstrated that ubiquitination played important roles in LSCC pathophysiological processes, and that ubiquitination was not only involved in the synthesis-degradation process of protein, but also in other multiple cancer-related signaling pathways.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970843/fendo-13-970843-HTML-r1/image_m/fendo-13-970843-g006.jpg)
Figure 6 Statistically significant signaling pathways identified with KEGG pathway analysis of 627 UPs. A total of 47 KEGG pathways was significantly enriched from those UPs (p < 0.05; FDR < 0.05).
As described in Table 1, multiple proteins are decreased or increased in LSCC tissues relative to control tissues, which might be due to the alteration of ubiquitination-mediated protein synthesis-degradation system. (i) Ribosome, the central site for protein synthesis or translation, is a dense ribonucleoprotein particle that includes two subunits (large and small), and each subunit contains one rRNA and many different protein molecules. The small subunit is responsible for sequence-specific recognition of the template mRNA, such as the recognition of the initial part, the interaction of codon with anti-codon, and the binding site of mRNA is also on this subunit. The large subunit is responsible for carrying amino acids and tRNAs, including peptide bond formation, and peptidyl-RAN binding. In this study, 7 UPs [ribosomal proteins L3e (P39023), L11e (P62913), L9e (P32969), L7Ae (P62424), L40e (P62987), L18Ae (Q02543), and L12e (A0A024RBS2)] were identified in large subunit with significantly increased ubiquitination level in ribosomal proteins L11e, L12e, L18Ae, and L40e in LSCC tissues relative to controls, and 6 UPs [ribosomal proteins S20e (P60866), S3e (P23396), S2e (P15880), S16e (P62249), S27Ae (B2RDW1), and S9e (A0A024R4M0)] in small subunit with significantly increased ubiquitination level in ribosomal proteins S3e, S9e, S16e, and S27e in LSCC tissues relative to controls (Figure 7). (ii) Endoplasmic reticulum (ER) is the cellular organelle, and is responsible for the assembly of multi-subunit proteins, formation of proper conformation of protein, protein secretion, lipid biosynthesis, and calcium homeostasis. The correctly folded proteins are transited into the Golgi complex. ER stress occurs under normal or pathophysiological conditions, which can lead to accumulation of misfolded proteins in the ER lumen to restore the correct fold. Proteins that ultimately failed to recover to normal folding were transported to the proteasome for ubiquitination-mediated degradation, which is named as ER-associated degradation (ERAD) (31). Accumulation of misfolded proteins can also lead to ER stress and ultimately lead to the activation of a range of signaling pathways, which was called UPR (unfolded protein response). In the case of severely damaged cells, UPR is insufficient to restore the function of ER, when the cells undergo apoptosis. In this study, 13 UPs [EIF2AK2 (P19525), CAPN1 (P07384), BCAP31 (P51572), EL52 (K9JA46), YOD1 (Q5VVQ6), UBE4B (O95155), RAD23B (P54727), NGLY1 (Q96IV0), VCP (P55072), HSPA8 (P11142), DNAJB1 (P25685), HSP90AB1 (A0A024RD80), and RAD23A (A0A024R7G8)] with 39 ubiquitination sites were identified in protein ER processing pathway, wth significantly decreased ubiquitination levels in proteins YOD1, and VCP (K8, K18, K109, K389, K651, K658, and K668), and with significantly increased ubiquitination levels in proteins BCAP31, EL52, RAD23B, VCP(K217, and K614), HSPA8, and DNAJB1 (Figure 8). (iii)Ubiquitination was multiple enzymatic reaction process that requires enzymes E1, E2, and E3. In this study, 9 Ups, including six E3s [UBE4B (O95155), CUL5 (Q93034), ITCH (Q96J02), CUL4A (Q13619), PML (P29590), and HUWE1 (A0A024R9W5)], two E2s [UBE2O (Q9C0C9), and UBE2N (P61088)], and one E1 [UBE1 (A0A024RDB0)], were identified, with significantly decreased ubiquitination levels in E2 UBE2N, and with significantly increased ubiquitination levels at E3s ITCH, HUWE1, and CUL4A, in LSCC tissues relative to controls (Figure 9). For these six E3s, HUWE1 and ITCH were HECT type E3, UBE4B were U-box type E3, PML was single RING-finger type E3, and CUL5 and CUL4A were multi subunit RING-finger type E3. (iv) Proteasome is the crucial machine to response for ubiquitin-mediated proteolysis, which includes one 20S core, one 19S regulatory lid (PA700), and one 19S regulatory base (PA700) in the 26S proteasome complex. The 20S core particle provides an enclosed cavity where proteins are degraded. This present study identified 7 UPs (Rpn13, and Rpt 1-6) in PA700 (Base) with significantly increased ubiquitination levels in proteins Rpn13, Rpt1, Rpt3, Rpt4, Rpt5 (K276 and K372), and Rpt6, and significantly decreased ubiquitination level in protein Rpt 5 (K53); and 5 UPs (Rpn 3, 5, 6, 10, and 12) in PA700 (Lid) with significantly increased ubiquitination level, in human LSCC tissues relative to controls (Figure 10), but no any ubiquitination was found in the 20S core.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970843/fendo-13-970843-HTML-r1/image_m/fendo-13-970843-t001.jpg)
Table 1 Differentially expressed proteins with significantly ubiquitination alterations in LSCC tissues compared to controls.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970843/fendo-13-970843-HTML-r1/image_m/fendo-13-970843-g007.jpg)
Figure 7 Ubiquitination-mediated ribosome pathway. Pink = the identified ubiquitinated subunits. Green=non-ubiquitinated subunits. K* = ubiquitinated lysine residue. Red indicates the degree of ubiquitination of this site is upregulated, and blue indicates downregulated. Ratio(T/N) = Ratio of tumor to control. T+/N- = ubiquitination only in tumor. NA* means no quantitative information in both tumor and control. # means statistically significantly altered ubiquitination level in tumor compared to controls.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970843/fendo-13-970843-HTML-r1/image_m/fendo-13-970843-g008.jpg)
Figure 8 Ubiquitination-mediated endoplasmic reticulum protein processing pathway. Pink = the identified ubiquitinated subunits. Green=non-ubiquitinated subunits. K* = ubiquitinated lysine residue. Red indicates the degree of ubiquitination of this site is upregulated, and blue indicates downregulated. Ratio(T/N) = Ratio of tumor to control. T+/N- = ubiquitination only in tumor. T-/N+ = ubiquitination only in control. NA* means no quantitative information in both tumor and control. # means statistically significantly altered ubiquitination level in tumor compared to controls.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970843/fendo-13-970843-HTML-r1/image_m/fendo-13-970843-g009.jpg)
Figure 9 Ubiquitination-mediated UPS-related enzymes. Pink = the identified ubiquitinated subunits. Green = non-ubiquitinated subunits. K* = ubiquitinated lysine residue. Red indicates the degree of ubiquitination of this site is upregulated, and blue indicates downregulated. Ratio(T/N) = Ratio of tumor to control. T+/N- = ubiquitination only in tumor. NA* means no quantitative information in both tumor and control. # means statistically significantly altered ubiquitination level in tumor compared to controls.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970843/fendo-13-970843-HTML-r1/image_m/fendo-13-970843-g010.jpg)
Figure 10 Ubiquitination-mediated proteasome pathway. Pink=the identified ubiquitinated subunits. Green = non-ubiquitinated subunits. K* = ubiquitinated lysine residue. Red indicates the degree of ubiquitination of this site is upregulated, and blue indicates downregulated. Ratio(T/N) = Ratio of tumor to control. T+/N- = ubiquitination only in tumor. NA* means no quantitative information in both tumor and control. # means statistically significantly altered ubiquitination level in tumor compared to controls.
Overlapping analysis of differentially ubiquitinated proteins and differentially expressed proteins in LSCC
One of the consequences of ubiquitination is to degrade proteins. Those DEPs with differentially ubiquitinated modifications might have more relevance to tumorigenesis. A set of 265 DEP data in LSCC tissues compared to controls extracting from the published data (32–36) (Supplementary Table 12) were compared to this set of 464 DUP data with 789 differentially ubiquitinated sites in LSCC tissues compared to controls (Supplementary Table 13). A total of 45 DEPs in LSCC were found to have significantly differentially ubiquitinated modifications (Table 1). When one analyzed the relationship between DUP and DEP in LSCC, those 45 DEPs were grouped into 5 categories (Table 2): (i) When the ubiquitination level was decreased, the protein expression level was increased, including APOA1, HBD, CA1, VIM, IDE, and GNG5; (ii) When the ubiquitination level was increased, the protein expression level was decreased, including KRT18, HSPA1B, SNAP23, FADS2, ALDOA, UCHL1, LMNB1, MAP2K3, RPS16, EEF1A1, and MYO1E; (iii) When the ubiquitination level was decreased, the protein expression level was decreased, including MYO6, PRKAR1A, EZR, HBA2, DES, YWHAB, and erythrocytic spectrin beta; (iv) When the ubiquitination level was increased, the protein expression level was increased, including MYH9, CYBB, ALDOC, GAPDH, HSPB1, GNAS, PTMA, ATP5B, IGF1R, CD44, CBR1, CAD, SHC1, GCLM, TMP4, EFNB1, DRIP4, and NADSYN1; and (v) when a protein existed both increased and descresed ubiquitination sites, the protein expression level was increased, including ANXA1, S100A9, and HNRPK.
Discussion
Protein ubiquitination profile and ubiquitination-mediated signaling pathways in LSCC
To investigate ubiquitination profile and its potential functions in LSCC, and to promote PPPM of LSCC, label-free quantitative proteomics was used to investigate ubiquitination in LSCC tissues. This present study provided the first comprehensive, quantitative ubiquitination profile in LSCC tissues. A total of 627 UPs with 1209 ubiquitination sites were identified, and 93.7% of identified ubiquitination sites had quantitative intensity. These UPs were involved in multiple functional categories, cellular biological processes, and signaling pathways in LSCC. In total, 47 statistically significant pathways were identified to involve UPs in LSCC, and most of these signaling pathways were cancer-related pathways, including ubiquitin-mediated proteolysis, ribosome complex, ER-mediated protein processing, proteasome complex, biosynthesis of amino acids, PI3K-Akt signaling pathway, Ras signaling pathway, mTOR signaling pathway, HIF-1 signaling pathway, cell cycle, apoptosis pathway, glycolysis/gluconeogenesis, carbon metabolism, central carbon metabolism in cancer, fructose and mannose metabolism, metabolomic pathways, drug metabolism-other enzymes, sphingolipid signaling pathway, tight junction, adherens junction, gap junction, non-homologous end-joining, endocytosis, SNARE interactions in vesicular transport, dopaminergic synapse, GABAergic synapse, chemokine signaling pathway, Rap1 signaling pathway, proteoglycans in cancer, pathways in cancer, microRNAs in cancer, leukocyte transendothelial migration, regulation of actin cytoskeleton, and platelet activation. These ubiquitinaiton-mediated signaling pathways provide the functional profiles of ubiquitination in LSCC.
Ubiquitination-mediated tumor protein turnover in LSCC
It is well-known that one of main functions of ubiquitination is involved in protein degradation, and participates in protein turnover in tumor pathophysiological processes. Interestingly, this study found that ubiquitination was significantly involved in four protein turnover-associated pathways, including ribosome complex, ubiquitin-mediated proteolysis, ER-mediated protein processing, and proteasome complex pathways, which are discussed in detailed here.
Proteins are vital parts of living organisms, and are in dynamic balance between generation and degradation. Protein turnover is the protein synthesis-degradation balance, which gives cells the flexibility to adjust the abundance and function of the protein in response to various intracellular or extracellular stimuli (37). Protein synthesis is determined by both transcription and translation. The mRNA transcribed into the cytoplasm and binds to the ribosome to begin the process of protein synthesis. The overall activity of ribosome and the rate of translation initiation and elongation have an important effect on the rate of protein synthesis. The newly synthesized protein does not have a spatial structure and corresponding functions, which needs to be folded and assembled in ER. The processing time of different proteins in ER is also different. The UPS and autophagy are the main pathways for protein degradation (38, 39). The former selectively degrades the misfolded proteins and short half-life proteins that have been labeled with ubiquitin in the proteasome, while the latter transports denatured long-lived proteins and damaged or excess organelles to the lysosome for degradation. The ability to adjust protein abundances in a timely and precise manner based on external environmental stimuli is critical to maintain the normal function of cells. Protein turnover abnormalities are associated with the development of multiple diseases, including neurodegenerative diseases (40, 41) and tumors (42, 43). Targeting altered protein turnover provides new opportunities and challenges to develop anti-tumor drugs (44).
Studies demonstrate that multiple ribosomal subunits have been ubiquitinated in eukaryotes (45–47). This study also found 13 ubiquitinated ribosome subunits in LSCC, which demonstrates the important regulatory roles of ubiquitination in ribosomal complex. Ubiquitination affects protein turnover by affecting the function of ribosomes in several ways.
First, ubiquitination regulates multiple steps of ribosome biogenesis (48), which involves the production and correct assembly of rRNAs and ribosomal proteins. Abnormal ribosome biogenesis inevitably affects the turnover of intracellular proteins. However, this study did not find the ribosome biogenesis pathway in KEGG enrichment analysis, which suggests that ubiquitination might not affect protein turnover by participating in regulating ribosome biogenesis in LSCC. Moreover, ubiquitination of ribosomal proteins regulates ribosome-associated quality control (RQC) (49, 50). Ribosomes are not only responsible for protein synthesis, but also execute RQC to minimize production of aberrant proteins. Many factors cause abnormal protein synthesis in ribosomes. RQC is a biologically evolved surveillance mechanism that timely terminates the protein synthesis of ribosomes, and initiates corresponding pathways to degrade nascent polypeptides when there exist interrupted translation (51). Obviously, dysfunction of RQC results in accumulating abnormal proteins in cells, which eventually affects the protein turnover. Study showed that the poly (A) tails of mature mRNA stalled ribosomes, repressed downstream translation, and initiated RQC; this is a process known as ribosome stall resolution (49). One study showed that in the initiation of this process, it is essential for ubiquitin ligase ZNF598 to catalyze ubiquitination at residue K8 in RPS20 (P60866) and at residues K138 and K139 in RPS10, and the ubiquitination failure of RPS20 or RPS10 would lead to defective resolution of stalled ribosomes (50). Another study also found that ZNF598 primarily mono-ubiquitinated two lysine residues K4 and K8 in RPS20, which was required to stall ribosomes during poly (A) translation (49). In the yeast model, study found that ubiquitination at residue K8 in 40S ribosomal protein uS10 (corresponding to RPS20 in human) catalyzed by E3 ubiquitin ligase Hel2 switched on RQC (52). Those studies showed that ubiquitination at residue K8 in RPS20 was species conservation, and played a pivotal role in RQC pathways. Besides RPS20, recent study demonstrated that the balance between mono-ubiquitination at residue K214 in RPS3 (P23396) catalyzed by RNF123 E3 ligase and deubiquitination by USP10, also participated in the regulation of RQC (53). This study found ubiquitination level was decreased at residue K8 (T/N ratio=0.95, P value = 0.78) in RPS20 and increased at residue K214 (T/N ratio=5.03, P value = 0.32) in RPS3, which means that the impaired RQC in LSCC tissue may contribute to the abnormal protein turnover. Further, ribosome ubiquitination participated in reprogramming of ribosomal protein translation that was induced by unfolded protein reaction (UPR). UPR not only induced ER stress and disturbed protein homeostasis, but also reprogramed translation; and persistent UPR might lead to cell death. Reprogramming of ribosomal protein translation is an intracellular response to protein homeostasis, and it can also affect the turnover of intracellular proteins. One study demonstrated that ubiquitination at residues K58 or K275 in RPS2 (P15880) and at residue K8 in RPS20 participate in UPR-induced ribosomal translation reprogram, and these two proteins without ubiquitination at these residues promoted UPR-induced cell death (54). This study identified ubiquitination at residues K58 (T/N = 1.4; P value = 0.06) and K275 (T/N = 2.21; P value = 0.21) in RPS2 increased in LSCC. Therefore, ubiquitination in proteins RPS3, RPS20, and RPS2 in 40S subunits play an important role in ribosome-related protein turnover regulation.
Both ribosome and ER are indispensable organelles for intracellular protein homeostasis and protein turnover, and this study demonstrates that ubiquitination also has regulatory effects on ER. Among the various biological processes of ER, UPR and ERAD are essential for maintaining protein homeostasis and turnover. The former regulates the expression of multiple downstream genes through at least three pathways to recover UPR-induced ER stress (55). The latter relieves ER stress through the process of retro-translocation (terminally misfolded proteins were translocated from the ER lumen into the cytoplasm through the translocation pore in the ER membrane) and ultimately degraded by proteasome (56). Our results showed that UPs mainly concentrated in the above two biological processes. For example, study demonstrated that EIF2AK2 (P19525) was responsible for starting a branch of the UPR. When ER stress was activated, EIF2AK2 itself was oligomerized and phosphorylated, and also ubiquitously translational initiation factor eIF2a was indirectly inactivated to inhibit translation of mRNA. Thus, EIF2AK2 inhibits the flux of protein to enter ER to alleviate ER stress and restore protein homeostasis (55). This study qualitatively identified the ubiquitination at residue K408 in EIF2AK2. Currently, the exact effect of ubiquitination on the function of EIF2AK2 is not clear, but study reported that the residues K69 and K159 have been ISGylated, one type of ubiquitin-like modification, that down-regulated protein translation (57). During ERAD, studies showed that molecular chaperones not only assist in folding of the protein in ER, but also play a variety of important functions to ensure ERAD (58). These molecular chaperones assist degradation through substrate recognition and preventing substrate aggregation (59). This study identified the ubiquitination of three molecular chaperones (HSPA8, DNAJB1 and K9JA46) indicating that ubiquitination can affect the function of molecular chaperones, which in turn affects ERAD and protein turnover.
The UPS is another determinant of intracellular protein turnover, because the degradation of many proteins is carried out by UPS in mammalian cells (60). Under physiological conditions, the activity of UPS is dynamically regulated by different signaling pathways in response to various stimuli (60). As mentioned above, UPS contains a large number of components, so dysfunction of any one component can lead to abnormal protein degradation, which eventually affects protein turnover. In this study, it was found that 11 proteasome subunits (all belong to the 19S regulatory subunit) and 9 UPS-related enzymes were ubiquitinated, which revealed that ubiquitination can regulate protein turnover by regulating proteasome activity and UPS-related enzymes.
The 26S proteasome is a core component of UPS, and proteins labeled with K48-linked polyubiquitin chains are recognized by its 19S regulatory subunit (RS) and degraded by its 20S core subunit (CS) (61, 62). A variety of factors affect proteasome activity, such as proteasome biogenesis and PTMs (63). For example, the RS has two subunits (lid an base), among which the base assembly is associated with 4 RS assembly chaperones (RACs): Rpn14 (PAAF1), S5b (PSMD5), p28 (PSMD10), and p27 (PSMD9) (64). Loss of these chaperones leads to base subcomplex of RS assembly defect (65). This study identified ubiquitination at residues K23 (no quantitative information) and K30 (T+/N-) in PSMD10 (O75832) in LSCC for the first time. We hypothesize that ubiquitination may affect the intracellular abundance or function of PSMD10, which finally affects proteasome biogenesis and the activity of proteasome. The specific regulatory mechanism needs further research. PTMs of various proteasome subunits make the regulation of proteasome function more complicated. Ubiquitination is one of the 11 PTMs of proteasome subunits that have been identified so far (66). At present, phosphorylation has been deeply studied in the regulation mechanism of proteasome activity, and researches had revealed that a variety of kinases and phosphatases regulate proteasome (67). The functional consequences of other PTMs (including ubiquitination) in regulation of proteasome are presently lacking. In our study, all ubiquitinated proteasome subunits belong to the RS, which suggests that ubiquitination may affect the activity of the proteasome by regulating the function of the RS not the CS in LSCC. Under normal circumstances, Rpn13 and Rpn10/S5a are responsible for identification of the ubiquitin chains on the protein committed to degradation, after that protein is translocated through the ATPase ring into the CS (68). However, when the proteasome is impaired under various circumstances (such as proteotoxic stresses and treated with bortezomib), Ube3c ubiquitinates residues K21 and K34 in Rpn13 and inhibits the binding of 26S particles to ubiquitin conjugates (69). Our study found the increased ubiquitination level at residue K34 (T/N=2.27; p = 0.01) in Rpn13, which can be seen as a sign of impaired proteasome function in LSCC. More importantly, if we can confirm the above results in LSCC, the degree of ubiquitination of Rpn13 can be used as a biomarker to predict the severity of impaired proteasome function, and the biological effects produced by proteasome activity inhibitors in vivo. In addition to the ubiquitination of Rpn13, monoubiquitination of Rpn10 induced by Rsp5, a member of NEDD4 ubiquitin protein ligase family, significantly inhibits the interaction of Rpn10 with substrates to decrease the activities of proteasome complex (70). Different from the previous yeast-based research that discovered four ubiquitination sites (K71, K84, K99, and K268), this study identified significantly increased ubiquitination levels at residues K40 (N/T = 33.8; p < 0.01) and K74 (N/T = 5.2, p = 0.01) in Rpn10 in LSCC tissue. Both K40 and K74 in Rpn10 are located in the VWFA domain that has been implicated to be responsible for stabilizing the lid-base association of RS. In LSCC, monoubiquitination of K40 and K74 in Rpn10 may participate in the regulation of proteasome activity.
In addition to proteasome biogenesis and PTMs, the amount of some proteasome subunits can also affect the activity of proteasome. For example, the activity of the proteasome is positively correlated with the expression level of PSMD11, and overexpression of PSMD11 increases the activity of the proteasome in stem cells (71). Specially, the dysregulated proteasomal activities were showed in the muscles of PSMC4-knockout animals (72). The ubiquitination of both PSMD11 and PSMC4 have been identified in this study. If these proteins are degraded by UPS, which means ubiquitination of these proteins will affect the activity of proteasome. So it is especially important to clarify the specific ubiquitination regulatory mechanisms of these proteins. As mentioned above, few studies are involved in the impact of ubiquitination on proteasome subunits, ubiquitination proteomics can be used as a powerful screening tool to provide direction for further studies.
Above description discussed the effect of ubiquitination on proteasome activity. Further, we discuss the ubiquitination regulation on the UPS-related enzymes. Our study found 9 UPS-related enzymes including E1, E2 and E3, which indicate that enzymes in the process of protein ubiquitination are also regulated by ubiquitination. Currently, studies have focused on the roles of ubiquitin in these enzymes; however, just relatively few studies were focused on the roles of the ubiquitination in these 9 enzymes. For instance, PML (P29590) was ubiquitinated by several E3s to result in the degradation of proteasome complex (73, 74). ITCH (Q96J02) can be ubiquitinated by itself via lysine-63 linkages to control the cytoplasmic-nuclear shuffling of ITCH (75). Of course, ubiquitination and even other PTMs might extensively affect these enzyme activities, currently known ubiquitination are only the tip of the iceberg.
The relationship between ubiquitination and protein expression level in LSCC
One of main functions of ubiquitination is to degrade proteins. A total of 265 DEP data in LSCC tissues was obtained from the published data (32–36) (Supplementary Table 12). This study identified 464 DUP data with 789 differentially ubiquitinated sites in LSCC tissues compared to controls (Supplementary Table 13). Overlapping analysis of the DUP data and DEP data found 45 DEPs with significantly altered ubiquitination level (Table 1). Five types of relationships between DUP and DEP were found in LSCC (Table 2), including: (i) the protein expression level was increased with the decrease of ubiquitination level, which might be due to the constant synthesis velocity and the decreased degradation velocity; (ii) the protein expression level was decreased with the increase of ubiquitination level, which might be due to the constant synthesis velocity and the increased degradation velocity; (iii) the protein expression level was decreased with decreased of the ubiquitination level, which might be because the protein degradation velocity was faster than its synthesis velocity; (iv) the protein expression level was increased with the increase of the ubiquitination level, which might be because protein degradation velocity was slower its synthesis velocity; (v) the protein expression level was increased with both increase and decrease of ubiquitination sites in a protein, this type of reseason was complex. These results clearly demonstrate the complexity of protein synthesis-degradation system. Further experiment studies would be needed to confirm the real status. Moreover, one should note that only 265 DEPs are obtained due to a limited list of DEPs in the published literature, an expanded quantitative proteomics would be needed to identify more DEPs in LSCC tissues compared to controls, which will reveal more DEPs with significant ubiquitination alteration.
Conclusions
This study provides the first quantitative ubiquitination proteomics analysis of LSCC tissue, and shows that ubiquitination can affect the occurrence and development of LSCC in multiple biological functions and pathways. Here focuses on the regulation mechanism of ubiquitination on intracellular protein turnover-related pathways. The results showed that ubiquitination affected UPS functions by regulating proteasome activity and UPS-related enzymes. In addition, ubiquitination also affected protein turnover by regulating ribosome assembly, ribosome-associated quality control (RQC), ER stress response, and ER-associated degradation (ERAD). Although ubiquitination studies in tumors have made significant progress, many anti-tumor drugs have been clinically applied and have achieved excellent therapeutic effects, yet considering that ubiquitination has extensive regulatory networks and the complexity of ubiquitination itself, so there is still a lot of unknown things waiting for ones to explore. This study paved the way for further exploration of the specific regulatory mechanisms of ubiquitination in LSCC.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving human participants were reviewed and approved by Xiangya Hospital Medical Ethics Committee of Central South University, China. The patients/participants provided their written informed consent to participate in this study.
Author contributions
XHZ conceived the concept, designed experiments and manuscript, instructed experiments and data analysis, supervised results, coordinated, wrote and critically revised manuscript, and was responsible for its financial supports and the corresponding works. ML collected tissue samples, analyzed ubiquitinomic data, prepared figures and tables, and wrote the manuscript. XZ and LY analyzed the relationship between ubiquitination and protein differential expression. JY, XHZ, SZ, YG, BL, SW, JL, and NL participated in partial experiments. All authors approved the final manuscript.
Funding
This work was supported by the Shandong First Medical University Talent Introduction Funds (to XZ) , Shandong First Medical University High-level Scientific Research Achievement Cultivation Funding Program (to XZ), the Shandong Provincial Natural Science Foundation (ZR202103020356 or ZR2021MH156 to XZ), and the Academic Promotion Program of Shandong First Medical University (2019ZL002).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2022.970843/full#supplementary-material
Abbreviations
ALK, anaplastic lymphoma kinase; AUC, area under the curve; BP, biological processes; CC, cellular components; CS, core subunit; DTT, dithiothreitol; E1, ubiquitin-activating enzyme; E2, ubiquitin-conjugating enzyme; E3, ubiquitin ligase; EPC, edge percolated component; ER, endoplasmic reticulum; GO, gene ontology; IAP, immunoaffinity purification; ICAT, isotope-coded affinity tags; iTRAQ, isobaric tags for relative and absolute quantification; KEGG, Kyoto Encyclopedia of Genes and Genomes; LC, liquid chromatography; LSCC, lung squamous cell cancer; MCC, maximal clique centrality; MCODE, molecular complex detection; MF, molecular functions; MNC, maximum neighborhood component; MS/MS, tandem mass spectrometry; NSCLC, non-small-cell lung cancer; PMSF, phenylmethyl-sulfonyl fluoride; PPI, protein-protein interaction; PPPM, predictive, preventive and personalized medicine; PTM, post-translational modification; RACs, RS assembly chaperones; RQC, ribosome-associated quality control; RS, regulatory subunit; SILAC, isotope labels with amino acids in cell culture; TAC, tumor-adjacent control; TFA, Trifluoroacetic acid; TKI, tyrosine kinase inhibitors; TTP, time to progression; UPs, ubiquitinated proteins; 2D-DIGE, 2D fluorescence difference in-gel electrophoresis; 2D-PAGE, two-dimensional polyacrylamide gel electrophoresis.
References
1. Gandara DR, Hammerman PS, Sos ML, Lara PN, Hirsch FR. Squamous cell lung cancer: from tumor genomics to cancer therapeutics. Clin Cancer Res (2015) 21(10):2236–43. doi: 10.1158/1078-0432.CCR-14-3039
2. Socinski MA, Obasaju C, Gandara D, Hirsch FR, Bonomi P, Bunn PA Jr., et al. Current and emergent therapy options for advanced squamous cell lung cancer. J Thorac Oncol (2018) 13:165–83. doi: 10.1016/j.jtho.2017.11.111
3. Zhang X-C, Wang J, Shao G-G, Wang Q, Qu X, Wang B, et al. Comprehensive genomic and immunological characterization of Chinese non-small cell lung cancer patients. Nat Commun (2019) 10:1772. doi: 10.1038/s41467-019-09762-1
4. Network CGAR. Comprehensive genomic characterization of squamous cell lung cancers. Nature. (2012) 489:519. doi: 10.1038/nature11404
5. Janssen-Heijnen ML, Coebergh J-WW. The changing epidemiology of lung cancer in Europe. Lung cancer. (2003) 41:245–58. doi: 10.1016/S0169-5002(03)00230-7
6. Perez-Moreno P, Brambilla E, Thomas R, Soria J-C. Squamous cell carcinoma of the lung: molecular subtypes and therapeutic opportunities. Clin Cancer Res (2012) 18:2443–51. doi: 10.1158/1078-0432.CCR-11-2370
7. Zhang T, Shan L, Zhang D. Advances of molecular targeted drugs used in maintenance therapy of non-small cell lung cancer. Zhongguo fei ai za zhi= Chin J Lung cancer. (2010) 13:1070–3. doi: 10.3779/j.issn.1009-3419.2010.11.14
8. Ciechanover A. The ubiquitin-proteasome proteolytic pathway. Cell. (1994) 79:13–21. doi: 10.1016/0092-8674(94)90396-4
9. Kimura Y, Tanaka K. Regulatory mechanisms involved in the control of ubiquitin homeostasis. J Biochem (2010) 147:793–8. doi: 10.1093/jb/mvq044
11. Pickart CM, Fushman D. Polyubiquitin chains: polymeric protein signals. Curr Opin Chem Biol (2004) 8:610–6. doi: 10.1016/j.cbpa.2004.09.009
12. Wilkinson KD. Ubiquitination and deubiquitination: targeting of proteins for degradation by the proteasome. Semin Cell Dev Biol (2000) 11 (3):141–8. doi: 10.1006/scdb.2000.0164
13. Komander D, Rape M. The ubiquitin code. Annu Rev Biochem (2012) 81:203–29. doi: 10.1146/annurev-biochem-060310-170328
14. Kessler BM. Ubiquitin–omics reveals novel networks and associations with human disease. Curr Opin Chem Biol (2013) 17:59–65. doi: 10.1016/j.cbpa.2012.12.024
15. Adams J, Kauffman M. Development of the proteasome inhibitor Velcade™(Bortezomib). Cancer Invest (2004) 22:304–11. doi: 10.1081/CNV-120030218
16. Vij R, Siegel DS, Jagannath S, Jakubowiak AJ, Stewart AK, McDonagh K, et al. An open-label, single-arm, phase 2 study of single-agent carfilzomib in patients with relapsed and/or refractory multiple myeloma who have been previously treated with bortezomib. Br J haematology. (2012) 158:739–48. doi: 10.1111/j.1365-2141.2012.09232.x
17. Morrow JK, Lin H-K, Sun S-C, Zhang S. Targeting ubiquitination for cancer therapies. Future medicinal Chem (2015) 7:2333–50. doi: 10.4155/fmc.15.148
18. Danielsen JM, Sylvestersen KB, Bekker-Jensen S, Szklarczyk D, Poulsen JW, Horn H, et al. Mass spectrometric analysis of lysine ubiquitylation reveals promiscuity at site level. Mol Cell Proteomics. (2011) 10:M110.003590. doi: 10.1074/mcp.M110.003590
19. Udeshi ND, Mertins P, Svinkina T, Carr SA. Large-Scale identification of ubiquitination sites by mass spectrometry. Nat Protoc (2013) 8:1950. doi: 10.1038/nprot.2013.120
20. Zhu W, Smith JW, Huang C-M. Mass spectrometry-based label-free quantitative proteomics. J BioMed Biotechnol (2010) 2010:840518. doi: 10.1155/2010/840518
21. Neilson KA, Ali NA, Muralidharan S, Mirzaei M, Mariani M, Assadourian G, et al. Less label, more free: approaches in label-free quantitative mass spectrometry. Proteomics. (2011) 11:535–53. doi: 10.1002/pmic.201000553
22. Low TY, Magliozzi R, Guardavaccaro D, Heck AJ. Unraveling the ubiquitin-regulated signaling networks by mass spectrometry-based proteomics. Proteomics. (2013) 13:526–37. doi: 10.1002/pmic.201200244
23. Wu Q, Cheng Z, Zhu J, Xu W, Peng X, Chen C, et al. Suberoylanilide hydroxamic acid treatment reveals crosstalks among proteome, ubiquitylome and acetylome in non-small cell lung cancer A549 cell line. Sci Rep (2015) 5:9520. doi: 10.1038/srep09520
24. Qin X, Chen S, Qiu Z, Zhang Y, Qiu F. Proteomic analysis of ubiquitination-associated proteins in a cisplatin-resistant human lung adenocarcinoma cell line. Int J Mol Med (2012) 29:791–800. doi: 10.3892/ijmm.2012.912
25. Wei Huang D, Sherman BT, Tan Q, Collins JR, Alvord WG, Roayaei J, et al. The DAVID gene functional classification tool: a novel biological module-centric algorithm to functionally analyze large gene lists. Genome Biol (2007) 8:1–16. doi: 10.1186/gb-2007-8-9-r183
26. Wu J, Mao X, Cai T, Luo J, Wei L. KOBAS server: a web-based platform for automated annotation and pathway identification. Nucleic Acids Res (2006) 34:W720–W4. doi: 10.1093/nar/gkl167
27. Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, et al. STRING v9. 1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res (2012) 41:D808–D15. doi: 10.1093/nar/gks1094
28. Smoot ME, Ono K, Ruscheinski J, Wang P-L, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics (2010) 27:431–2. doi: 10.1093/bioinformatics/btq675
29. Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinf (2003) 4:2. doi: 10.1186/1471-2105-4-2
30. Chin C-H, Chen S-H, Wu H-H, Ho C-W, Ko M-T, Lin C-Y. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol (2014) 8:S11. doi: 10.1186/1752-0509-8-S4-S11
31. McCracken AA, Brodsky JL. Assembly of ER-associated protein degradation in vitro: dependence on cytosol, calnexin, and ATP. J Cell Biol (1996) 132:291–8. doi: 10.1083/jcb.132.3.291
32. Zeng GQ, Zhang PF, Li C, Peng F, Li MY, Xu Y, et al. Comparative proteome analysis of human lung squamous carcinoma using two different methods: two-dimensional gel electrophoresis and iTRAQ analysis. Technol Cancer Res Treat (2012) 11(4):395–408. doi: 10.7785/tcrt.2012.500287
33. Li B, Chang J, Chu Y, Kang H, Yang J, Jiang J, et al. Membrane proteomic analysis comparing squamous cell lung cancer tissue and tumour-adjacent normal tissue. Cancer Lett (2012) 319(1):118–24. doi: 10.1016/j.canlet.2011.12.037
34. Xu Y, Cao LQ, Jin LY, Chen Z, Zeng G, Tang CE, et al. Quantitative proteomic study of human lung squamous carcinoma and normal bronchial epithelial acquired by laser capture microdissection. J BioMed Biotechnol (2012) 2012:510418. doi: 10.1155/2012/510418
35. Poschmann G, Sitek B, Sipos B, Ulrich A, Wiese S, Stephan C, et al. Identification of proteomic differences between squamous cell carcinoma of the lung and bronchial epithelium. Mol Cell Proteomics. (2009) 8(5):1105–16. doi: 10.1074/mcp.M800422-MCP200
36. Dowling P, O’Driscoll L, Meleady P, Henry M, Roy S, Ballot J, et al. 2-d difference gel electrophoresis of the lung squamous cell carcinoma versus normal sera demonstrates consistent alterations in the levels of ten specific proteins. Electrophoresis. (2007) 28(23):4302–10. doi: 10.1002/elps.200700246
37. Waterlow JC, Garlick PJ MILL W. Protein turnover in mammalian tissues and in the whole body. 335 Jan van Galenstraat, PO Box 211: Elsevier/North-Holland Biomedical Press (1978).
38. Ciechanover A. Proteolysis: from the lysosome to ubiquitin and the proteasome. Nat Rev Mol Cell Biol (2005) 6:79. doi: 10.1038/nrm1552
39. Klionsky DJ. Autophagy: from phenomenology to molecular understanding in less than a decade. Nat Rev Mol Cell Biol (2007) 8:931. doi: 10.1038/nrm2245
40. Riederer IM, Schiffrin M, Kövari E, Bouras C, Riederer BM. Ubiquitination and cysteine nitrosylation during aging and alzheimer’s disease. Brain Res bulletin. (2009) 80:233–41. doi: 10.1016/j.brainresbull.2009.04.018
41. Lim K-L. Ubiquitin–proteasome system dysfunction in parkinson’s disease: current evidence and controversies. Expert Rev proteomics. (2007) 4:769–81. doi: 10.1586/14789450.4.6.769
42. Deng S, Zhou H, Xiong R, Lu Y, Yan D, Xing T, et al. Over-expression of genes and proteins of ubiquitin specific peptidases (USPs) and proteasome subunits (PSs) in breast cancer tissue observed by the methods of RFDD-PCR and proteomics. Breast Cancer Res Treat (2007) 104:21–30. doi: 10.1007/s10549-006-9393-7
43. Kammerl IE, Caniard A, Merl-Pham J, Ben-Nissan G, Mayr CH, Mossina A, et al. Dissecting the molecular effects of cigarette smoke on proteasome function. J proteomics. (2019) 193:1–9. doi: 10.1016/j.jprot.2018.12.015
44. Demarchi F, Brancolini C. Altering protein turnover in tumor cells: new opportunities for anti-cancer therapies. Drug resistance updates. (2005) 8:359–68. doi: 10.1016/j.drup.2005.12.001
45. Xu G, Paige JS, Jaffrey SR. Global analysis of lysine ubiquitination by ubiquitin remnant immunoaffinity profiling. Nat Biotechnol (2010) 28:868. doi: 10.1038/nbt.1654
46. Kim W, Bennett EJ, Huttlin EL, Guo A, Li J, Possemato A, et al. Systematic and quantitative assessment of the ubiquitin-modified proteome. Mol Cell (2011) 44:325–40. doi: 10.1016/j.molcel.2011.08.025
47. Matsumoto M, Hatakeyama S, Oyamada K, Oda Y, Nishimura T, Nakayama KI. Large-Scale analysis of the human ubiquitin-related proteome. Proteomics. (2005) 5:4145–51. doi: 10.1002/pmic.200401280
48. Stavreva DA, Kawasaki M, Dundr M, Koberna K, Müller WG, Tsujimura-Takahashi T, et al. Potential roles for ubiquitin and the proteasome during ribosome biogenesis. Mol Cell Biol (2006) 26:5131–45. doi: 10.1128/MCB.02227-05
49. Juszkiewicz S, Hegde RS. Initiation of quality control during poly (A) translation requires site-specific ribosome ubiquitination. Mol Cell (2017) 65:743–50.e4. doi: 10.1016/j.molcel.2016.11.039
50. Sundaramoorthy E, Leonard M, Mak R, Liao J, Fulzele A, Bennett EJ. ZNF598 and RACK1 regulate mammalian ribosome-associated quality control function by mediating regulatory 40S ribosomal ubiquitylation. Mol Cell (2017) 65:751–60.e4. doi: 10.1016/j.molcel.2016.12.026
51. Brandman O, Hegde RS. Ribosome-associated protein quality control. Nat Struct Mol Biol (2016) 23:7. doi: 10.1038/nsmb.3147
52. Matsuo Y, Ikeuchi K, Saeki Y, Iwasaki S, Schmidt C, Udagawa T, et al. Ubiquitination of stalled ribosome triggers ribosome-associated quality control. Nat Commun (2017) 8:159. doi: 10.1038/s41467-017-00188-1
53. Jung Y, Kim HD, Yang HW, Kim HJ, Jang C-Y, Kim J. Modulating cellular balance of Rps3 mono-ubiquitination by both Hel2 E3 ligase and Ubp3 deubiquitinase regulates protein quality control. Exp Mol Med (2017) 49:e390. doi: 10.1038/emm.2017.128
54. Higgins R, Gendron JM, Rising L, Mak R, Webb K, Kaiser SE, et al. The unfolded protein response triggers site-specific regulatory ubiquitylation of 40S ribosomal proteins. Mol Cell (2015) 59:35–49. doi: 10.1016/j.molcel.2015.04.026
55. Walter P, Ron D. The unfolded protein response: from stress pathway to homeostatic regulation. Science. (2011) 334:1081–6. doi: 10.1126/science.1209038
56. Christianson JC, Ye Y. Cleaning up in the endoplasmic reticulum: ubiquitin in charge. Nat Struct Mol Biol (2014) 21:325. doi: 10.1038/nsmb.2793
57. Okumura F, Okumura AJ, Uematsu K, Hatakeyama S, Zhang D-E, Kamura T. Activation of double-stranded RNA-activated protein kinase (PKR) by interferon-stimulated gene 15 (ISG15) modification down-regulates protein translation. J Biol Chem (2013) 288:2839–47. doi: 10.1074/jbc.M112.401851
58. Bukau B, Weissman J, Horwich A. Molecular chaperones and protein quality control. Cell. (2006) 125:443–51. doi: 10.1016/j.cell.2006.04.014
59. Nishikawa S-I, Brodsky JL, Nakatsukasa K. Roles of molecular chaperones in endoplasmic reticulum (ER) quality control and ER-associated degradation (ERAD). J Biochem (2005) 137:551–5. doi: 10.1093/jb/mvi068
60. Huang X, Dixit VM. Drugging the undruggables: exploring the ubiquitin system for drug development. Cell Res (2016) 26:484. doi: 10.1038/cr.2016.31
61. Yu H, Matouschek A. Recognition of client proteins by the proteasome. Annu Rev biophysics. (2017) 46:149–73. doi: 10.1146/annurev-biophys-070816-033719
62. Finley D. Recognition and processing of ubiquitin-protein conjugates by the proteasome. Annu Rev Biochem (2009) 78:477–513. doi: 10.1146/annurev.biochem.78.081507.101607
63. Livneh I, Cohen-Kaplan V, Cohen-Rosenzweig C, Avni N, Ciechanover A. The life cycle of the 26S proteasome: from birth, through regulation and function, and onto its death. Cell Res (2016) 26:869. doi: 10.1038/cr.2016.86
64. Rousseau A, Bertolotti A. An evolutionarily conserved pathway controls proteasome homeostasis. Nature. (2016) 536:184. doi: 10.1038/nature18943
65. Le Tallec B, Barrault M-B, Guérois R, Carré T, Peyroche A. Hsm3/S5b participates in the assembly pathway of the 19S regulatory particle of the proteasome. Mol Cell (2009) 33:389–99. doi: 10.1016/j.molcel.2009.01.010
66. Kikuchi J, Iwafune Y, Akiyama T, Okayama A, Nakamura H, Arakawa N, et al. Co-And post-translational modifications of the 26S proteasome in yeast. Proteomics. (2010) 10:2769–79. doi: 10.1002/pmic.200900283
67. Guo X, Huang X, Chen MJ. Reversible phosphorylation of the 26S proteasome. Protein Cell (2017) 8:255–72. doi: 10.1007/s13238-017-0382-x
68. Peth A, Kukushkin N, Bossé M, Goldberg AL. Ubiquitinated proteins activate the proteasomal ATPases by binding to Usp14 or Uch37 homologs. J Biol Chem (2013) 288:7781–90. doi: 10.1074/jbc.M112.441907
69. Besche HC, Sha Z, Kukushkin NV, Peth A, Hock EM, Kim W, et al. Autoubiquitination of the 26S proteasome on Rpn13 regulates breakdown of ubiquitin conjugates. EMBO J (2014) 33:1159–76. doi: 10.1002/embj.201386906
70. Isasa M, Katz EJ, Kim W, Yugo V, González S, Kirkpatrick DS, et al. Monoubiquitination of RPN10 regulates substrate recruitment to the proteasome. Mol Cell (2010) 38:733–45. doi: 10.1016/j.molcel.2010.05.001
71. Vilchez D, Boyer L, Morantte I, Lutz M, Merkwirth C, Joyce D, et al. Increased proteasome activity in human embryonic stem cells is regulated by PSMD11. Nature. (2012) 489:304. doi: 10.1038/nature11468
72. Kitajima Y, Tashiro Y, Suzuki N, Warita H, Kato M, Tateyama M, et al. Proteasome dysfunction induces muscle growth defects and protein aggregation. J Cell Sci (2014) 127:5204–17. doi: 10.1242/jcs.150961
73. Yuan W-C, Lee Y-R, Huang S-F, Lin Y-M, Chen T-Y, Chung H-C, et al. A Cullin3-KLHL20 ubiquitin ligase-dependent pathway targets PML to potentiate HIF-1 signaling and prostate cancer progression. Cancer Cell (2011) 20:214–28. doi: 10.1016/j.ccr.2011.07.008
74. Lim JH, Liu Y, Reineke E, Kao H-Y. Mitogen-activated protein kinase extracellular signal-regulated kinase 2 phosphorylates and promotes Pin1 protein-dependent promyelocytic leukemia protein turnover. J Biol Chem (2011) 286:44403–11. doi: 10.1074/jbc.M111.289512
Keywords: lung squamous cell carcinoma, ubiquitination, ubiquitinated protein, signaling pathway, ubiquitin-proteasome system, biomarker
Citation: Zhan X, Lu M, Yang L, Yang J, Zhan X, Zheng S, Guo Y, Li B, Wen S, Li J and Li N (2022) Ubiquitination-mediated molecular pathway alterations in human lung squamous cell carcinomas identified by quantitative ubiquitinomics. Front. Endocrinol. 13:970843. doi: 10.3389/fendo.2022.970843
Received: 16 June 2022; Accepted: 25 August 2022;
Published: 15 September 2022.
Edited by:
Tadashi Nakagawa, Tohoku University, JapanReviewed by:
Yanfen Liu, ShanghaiTech University, ChinaFuminori Tokunaga, Osaka Metropolitan University, Japan
Copyright © 2022 Zhan, Lu, Yang, Yang, Zhan, Zheng, Guo, Li, Wen, Li and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xianquan Zhan, yjzhan2011@gmail.com