Prognostic value of SEC61G in lung adenocarcinoma: a comprehensive study based on bioinformatics and in vitro validation

Background Studies have shown that the Sec61 gamma subunit (SEC61G) is overexpressed in several tumors and could serve as a potential prognostic marker. However, the correlation between SEC61G and lung adenocarcinoma (LUAD) remains unclear. In the current study, we aimed to demonstrate the prognostic value and potential biological function of the SEC61G gene in LUAD. Methods Public datasets were used for SEC61G expression analyses. The prognostic value of SEC61G in LUAD was investigated using the Kaplan–Meier survival and Cox analyses. The correlation between the methylation level of SEC61G and its mRNA expression was evaluated via cBioPortal. Additionally, MethSurv was used to determine the prognostic value of the SEC61G methylation levels in LUAD. Functional enrichment analysis was conducted to explore the potential mechanism of SEC61G. Also, single sample GSEA (ssGSEA) and TIMER online tool were applied to identify the correlation between SEC61G and immune filtration. Furthermore, cell functional experiments were conducted to verify the biological behavior of SEC61G in lung adenocarcinoma cells (LAC). Results SEC61G was upregulated in pan-cancers, including LUAD. High SEC61G expression was significantly correlated with worse prognosis in LUAD patients. Multivariate analysis demonstrated that high SEC61G expression was an independent prognostic factor in the TCGA cohort. (HR = 1.760 95% CI: 1.297–2.388, p < 0.001). The methylation level of SEC61G negatively correlated with the SEC61G expression (R = − 0.290, p < 0.001), and patients with low SEC61G methylation had worse overall survival. (p = 0.0014). Proliferation-associated terms such as cell cycle and cell division were significantly enriched in GO and KEGG analysis. Vitro experiments demonstrated that knockdown of SEC61G resulted in decreased cell proliferation, invasion and facilitated apoptosis in LAC. GSEA analysis found that SEC61G expression was associated with the E2F targets. Moreover, SEC61G expression was negatively correlated with the immune cell infiltration including CD4+ T cell, CD8+ T cell, B cell, macrophage, neutrophil, and dendritic cell. Conclusion Our study indicated that overexpression of SEC61G was significantly associated with poor prognosis of LUAD patients and the malignant phenotypes of LUAD cells, suggesting that it could be a novel prognostic biomarker and potential therapeutic target of LUAD. Supplementary Information The online version contains supplementary material available at 10.1186/s12885-021-08957-4.

Conclusion: Our study indicated that overexpression of SEC61G was significantly associated with poor prognosis of LUAD patients and the malignant phenotypes of LUAD cells, suggesting that it could be a novel prognostic biomarker and potential therapeutic target of LUAD.
Keywords: SEC61G, Lung adenocarcinoma, Novel biomarkers, Prognosis, Proliferation, Invasion Background Lung cancer, one of the most common malignant tumors worldwide, remains the first leading cause of cancer deaths [1]. Lung adenocarcinoma (LUAD) is the most frequently diagnosed histological subtype of nonsmall cell lung cancer (NSCLC), followed by squamous cell carcinoma [2]. Various factors such as cigarette smoking, second-hand or passive smoking, air pollution, genetic alteration, asbestos, and radon put individuals under the risk of LUAD [3]. Although alternative treatments such as targeted therapy and immune therapy for LUAD patients have been progressed rapidly over the past decades, the average 5-year survival rate of patients with LUAD remains less than 20% [4,5]. A better understanding of the biological processes and molecular mechanisms underlying lung cancer pathogenesis will be of great significance to clinicians and improve patients' outcomes. Thus it is of vital importance to keep searching for new tumor biomarkers and other potential genetic targets [6].
SEC61G, also known as Sec61 gamma subunit, is a central member of the SEC61 complex, a heterotrimeric protein channel formed by three subunits, SEC61 α, β, and γ [7]. Sec61 combined Sec62 and Sec63 serves as core component of the protein translocation machinery in the endoplasmic reticulum (ER) membrane and participates in protein folding, post-translational modification, translocation and unfolded protein response (UPR), especially under conditions of ER stress such as hypoxia and glucose deprivation in the tumor microenvironment [8,9]. High frequency of mutation and overexpression of SEC62 and SEC63 have been observed in kinds of cancers, suggesting the potential role of ER protein in tumor development [10,11]. Several studies have demonstrated that SEC61G was overexpressed in Glioblastoma [12], gastric cancer [13], hepatocellular carcinoma [14,15] and breast carcinomas [16]. SEC61G gene was also found to coamplify with epidermal growth factor receptor in patients with glioblastoma and served as a potential prognostic marker [12,17]. However, the potential correlation between SEC61G and lung adenocarcinoma has not been characterized.
In this present study, we comprehensively investigated the prognostic impact of SEC61G expression in LUAD patients through the gene expression profile and the matching clinical information of LUAD patients from the Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database. A nomogram based on several independent risk factors was constructed. Functional enrichment and Gene set enrichment analysis (GSEA) were performed to explore the underlying mechanism of SEC61G involved in LUAD pathogenesis. We also identified the association between SEC61G expression and genetic alteration and methylation. Additionally, single-sample Gene Set Enrichment Analysis (ssGSEA) and TIMER were used to explore the correlation between the infiltration of immune cells and the expression of SEC61G. Finally, cell functional experiments were performed to verify the biological behavior of SEC61G in lung adenocarcinoma cells (LAC). We hope our study will help others get further insight into understanding the potential role of SEC61G in tumor pathogenesis and contribute to the improvement of molecular targeted therapy and prognosis for LUAD patients.

Data acquisition
The gene expression data, phenotype data, and the corresponding clinicopathological information of the TCGA-LUAD project and other tumor were acquired from the UCSC Xena browser (version: 2019-07-20,http://xenabrowser.net/datapages/). Clinical parameters, including age, gender, TNM stage and pathological stage, were evaluated. Transcriptome profiling data of patients with LUAD in the GSE11969 dataset from the Gene Expression Omnibus (https://www/ncbi.nlm.nih. gov/geo/) database were used for external validation of survival analyses. The two databases' exclusion criteria were as follows: cases without complete gene expression data and survival information. Finally, 497 cases were extracted for further analysis. Patients with LUAD were classified into low-and high-expression groups according to the optimal cut-off value of SEC61G.
The expression levels of SEC61G in various cancers including LUAD were acquired from the TIMER database [18]. (https://cistrome.shinyapps.io/timer/). The Broad Institute Cancer Cell Line Encyclopedia Database (CCLE; http://www.broadinstitute.org/ccle) was used to validate the expression levels of SEC61G in different types of cancer cell lines [19]. The datasets of IMvigor 210 Clinical Trial was also download via the IMvigor210CoreBiologies R package and used for the survival analyses [20].

Co-expressed and functional enrichment analysis
We used R (version 3.6.1) to identify the genes coexpressed with SEC61G in the TCGA-LUAD cohort. Spearman's correlation coefficient (r) was calculated to evaluate the correlation between SEC61G and co-expressed genes, of which |r| > 0.35 and P < 0.001 were selected. Meanwhile, co-expressed analysis was also performed to explore the correlation between SEC61G and ER-stress related genes which was obtained from the Molecular Signatures Database (https://www.gsea-msigdb.org/gsea/ msigdb/genesets.jsp). To explore the underlying biological mechanism of SEC61G in LUAD pathogenesis, Metascape [21](https://metascape.org) was used conducted gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis for genes co-expressed with SEC61G as previously selected.

Gene set enrichment analysis
Gene Set Enrichment Analysis (GSEA) is a computational method that determines whether a set of a priori defined genes show statistically significant and consistent differences between biological states [22]. In this study, We conducted GSEA using the "clusterProfiler" R package (3.8.0, [23]) to elucidate the statistically significant function and pathway difference between high and low SEC61G expression groups of LUAD. H.all.v7.0.symbols.gmt in the MSigDB Collections was used as the reference gene collection. The expression level of SEC61G was regarded as a phenotype label. Adjusted P-value < 0.001, FDR q-value < 0.001 and |NES| > 1.5 was considered as statistically significant.

Analysis between SEC61G expression and its correlation with immune infiltration
We conducted the ssGSEA (single-sample Gene Set Enrichment Analysis) method from the GSVA package [24] in R software to investigate the correlation between SEC61G expression and the immune cell infiltration level based on the published signature gene lists [25]. Spearman correlation was carried out to evaluate the correlation between SEC61G and immune cell infiltration. TIMER online tool was applied to validate the association between the SEC61G expression and the level of immune cell infiltration.

SEC61G methylation level and its prognostic analysis
The copy number variation (CNV) and methylation level data of SEC61G were acquired from the cBioPortal web platform (https://www.cbioportal.org/) and a comparison of the varying SEC61G gene expressions in SEC61G copy number variation groups (Kruskal-Wallis test) and the correlation between SEC61G methylation level and SEC61G gene expression (Spearman correlation) was conducted. The SMART online platform (http://www. bioinfo-zs.com/smartapp/) was used to visualize the methylation levels of SEC61G in pan-cancer and normal samples from the TCGA database. The UALCAN online tool (http://ualcan.path.uab.edu/) was used to identify the differences in the promoter methylation level of SEC61G between LUAD and normal tissues from TCGA data. MethSurv online tool (https://biit.cs.ut.ee/ methsurv/) was used to explore the prognostic value of the SEC61G methylation level in the TCGA-LUAD cohort.

Prognostic model generation and prediction
Univariate analysis and multivariate cox regression analysis were used to determine the optimal prognostic model. Then, a nomogram was constructed to predict the prognosis by R packages rms. The patients were stratified into high and low risk groups based on the optimal cut-off value. The difference in OS between the high-risk group and low-risk group was analyzed by the Kaplan-Meier method with a two-sided log-rank test. The concordance index(C-index) and calibration curves were used to evaluate nomogram models' quality. The C-index is between 0.5 and 1.0, where 1.0 indicates the model has a perfect capacity to distinguish outcomes with the model correctly, and 0.5 indicates random probability. The calibration curve is evaluated graphically by plotting the nomogram's predicted probability against the observed rates. Overlap with the reference line indicates that the consistency of the model is perfect.

Cell lines and cell culture
Human lung carcinoma cell lines A549 and H1299 were obtained from ATCC. A549 cells were cultured in F12K medium supplemented with 10% foetal bovine serum (FBS, AusGeneX, Australia), 100 U/ml penicillin, and 100 mg/ml streptomycin. H1299 cells were cultured in RPMI 1640 medium supplemented with 10% FBS, 100 U/ml penicillin, and 100 mg/ml streptomycin. Cell cultures were maintained in a humidified incubator consisting 5% CO2 at 37°C.

Cell counting kit-8 (CCK8) assay
After transfection, cells were seeded into a 96-well plate according to the standard of 5000 cells per well. The proliferation of A549 or H1299 cells in the si-SEC61G and NC group was detected respectively. Cell viability was measured every 12 h. Briefly, 10 μl CCK8 reagent (MCE, Shanghai, China) was added daily to each hole in the 96-well plates and incubated at 37°C for 1.5 h. The microplate reader (Infinite® M1000 PRO, TECAN, Switzerland) was used to measure the OD value at 450 nm.

Colony formation assay
Briefly, 1000 cells were seeded into a 6-well plate and incubated in culture medium with 10% foetal bovine serum at 37°C. The plates were incubated at 37°C, 5% CO2 incubator for 2 weeks with frequent observation. Then the supernatant was removed and cells were washed carefully twice with PBS. Cells were fixed with 4% paraformaldehyde for 0.5 h before dyeing with 0.1% crystal violet. The number of clones was counted directly by the naked eye. The size and number of clones were compared finally.

Cell apoptosis
The culture medium was replaced with a serum-free medium after 48 h of transfection., Cells were obtained by centrifugation after another 24 h of starvation and then resuspended by cold PBS (4°C) and then centrifuged again. The supernatant was carefully removed. Cells were spread around in the binding buffer and complied with the Annexin V-FITC staining manufacturer's instructions (US Everbright® Inc., Suzhou, China). FACS Calibur flow cytometer (BD Biosciences, Mountain View, CA) was used to determine the effects of SEC61G on lung cancer cell apoptosis.

Cell invasion assay
Cell invasion capability was conducted using the 24-well Transwell chambers (8 μm pore size; Corning, NY, USA). The transwell chambers were coated with 100 μl matrigel (5× dilution; 100 μL/well; BD Biosciences, Bedford, MA) in a 24-well plate at 37°C for 4 h. Then the upper compartment containing 100 μL of serum-free medium was added with 2 × 104 cells. Meanwhile, the lower compartment was added with 600 μL of medium with 10% FBS. After being cultured for 48 h in a 37°C incubator, cotton sticks were used to wipe off the cells remaining in the upper chamber, and the cells penetrating the filter were fixed with 4% paraformaldehyde and mixed with 0.4% crystal violet solution. Subsequently, the stained cells were observed and photographed under a microscope at 100× magnification (Olympus, Tokyo, Japan). At least four fields in each chamber were observed under the microscope. The average number of invasion cells in each microscopic fields was counted.

Statistical analysis
All statistical analyses were performed with the statistical package IBM SPSS Statistics software (SPSS26.0), Graph-Pad Prism 8.0 and R (version 3.6.1). The Wilcoxon signed-rank test and Wilcoxon rank-sum test were used to analyze the expression of SEC61G in paired and nonpaired samples. X-tile software (Version 3.6.1) was conducted to identify the optimal cut-off value. The relationships between clinicopathological features and SEC61G expression were evaluated using Wilcoxon signed-rank test and logistic regression. Univariate analysis and multivariate analysis were done by applying Cox logistic regression model to identify independent variables, including age, gender, T stage, N stage, M stage, pathological stage, and SEC61G expression. The 95% confidence interval (Cl) of HR was calculated to assess individual factors hazard risk. The "survival" R package (version:0.1.3) and the "survminer" R package (version 0.4.8) were used to draw survival curves. The DESeq2 package (version: 3.26.5) was used to identify differential expression analysis between low and high SEC61G expression groups. The student's t-test was carried out for comparisons between each group in vitro experiments. All tests were two-sided, and a P-value < 0.05 was considered as statistically significant.

SEC61G is overexpressed in lung adenocarcinoma
The results showed that SEC61G was highly expressed in lung adenocarcinoma tissue compared with normal tissue (p < 0.001) (Fig. 1A). In the paired specimens, SEC61G expression level in the LUAD group was significantly higher than that of the adjacent normal tissues (p < 0.001) (Fig. 1B). The ROC curve presented that the expression of SEC61G in LUAD was 0.887 (95% CI: 0.859-0.916) (Fig. 1C).
To further evaluate SEC61G expression in human cancers, we used TIMER to identify the SEC61G expression in multiple malignancies. The differential expression of SEC61G between the tumor and adjacent normal tissues is shown in Fig. 1D. SEC61G mRNA expression was significantly higher in bladder, breast, colorectal, esophageal, head and neck, kidney, liver, gastric, lung, prostate, gastric cancers and cholangiocarcinoma compared with the corresponding normal tissues (Fig. 1D). All tumor tissues presented in the TIMER database except the thyroid carcinoma showed higher SEC61G expression compared with the corresponding normal tissues. Then we investigated the expression levels of SEC61G in different cancer cell lines through CCLE database, as shown in the additional file 1.

High SEC61G mRNA expression correlates with adverse outcome in LUAD patients
We used the TCGA-LUAD cohort (n = 497) to investigate the correlation between SEC61G expression and LUAD patients' prognosis. Survival analyses was also performed in the other tumor types whose SEC61G expression was differential expressed with normal tissue, as shown in additional file 2. Baseline characteristics of lung adenocarcinoma patients in the TCGA dataset were shown in Table 1. Kaplan-Meier analysis demonstrated that high SEC61G expression group is significantly related to shorter overall survival (OS) (HR=1.74 (1.30-2.34), p < 0.001) (Fig. 2A). We also used the GSE11969 datasets from the GEO database for validation, which was consistent with the result from the TCGA cohort (HR = 1.64 (1.04-2.59), p = 0.033) (Fig. 2B). Also, Univariate and multivariate Cox regression analysis was conducted to examine whether the SEC61G expression was an independent factor in LUAD patients. (Table 2) Finally, age, pathological stage, and SEC61G expression were identified as the independent prognostic factors. This indicated a role of SEC61G in the prognosis of LUAD. The above results demonstrated that SEC61G is a prognostic factor, and increased SEC61G expression level correlates with poor OS.

Development of predictive nomogram based on SEC61G and clinicopathological factors
A nomogram integrating those independent clinical risk factors (age, pathological stage and SEC61G expression) was constructed (Fig. 2C). A high total score predicted the low 1-, 3-and 5-year survival. And a low total score showed the opposite. The C-index for OS prediction was 0.696, with 1000 bootstrap resamples for the nomogram. And the calibration plots (Fig. 2D) showed good agreement compared with the ideal curves, indicating that our assembled nomogram has the stability for predicting LUAD patient prognosis in clinical practice. The correlation between the mRNA expression of SEC61G and its copy number variation (CNV) was analyzed using cBioPortal. Patients with the amplification of CNV of SEC61G had a higher level of SEC61G expression in LUAD, but only 5.3% of patients (27/512) exhibited this (Fig. 3A). This indicated that CNV might not be the primary cause behind the high expressed SEC61G gene. We further investigated the association between SEC61G methylation and gene expression, and these results showed that SEC61G methylation negatively associated with SEC61G gene expression (R = − 0.290, p < 0.001) (Fig. 3B). Multiple malignancies showed a lower methylation level compared with normal tissues in the TCGA database, including lung adenocarcinoma (Fig. 3C). The UALCAN showed that the promoter methylation of SEC61G in lung adenocarcinoma tissues was significantly lower than that of adjacent normal tissues (p < 0.001) (Fig. 3D). In addition, the MethSurv online tool revealed that patients with lower SEC61G methylation had a worse overall survival. (p = 0.0014) (Fig. 3E). Furthermore, we also found that the expression of SEC61G was positively correlated with gene copy number and demethylation level in CCLE database, which was consistent with the results of our bioinformatics analysis, as shown in additional file 1.

Functional enrichment analysis of SEC61G co-expressed genes
A total of 973 co-expression genes were identified, of which 560 genes were positively correlated and 413 were negatively correlated. Following this, the functions of co-expression in patients with LUAD were predicted using GO enrichment and KEGG pathway analysis. The top 10 GO enrichment items (Fig. 4A) were including mitochondrial protein complex, purine nucleoside triphosphate metabolic process, mitotic cell cycle phase transition, cell division, chromosome, centromeric region, translational termination, DNA repair, transferase complex, mitochondrial intermembrane space, DNA replication. KEGG pathway analysis (Fig. 4B) showed enrichment in the proteasome, cell cycle, oxidative phosphorylation, glycolysis/gluconeogenesis, DNA replication, Homologous recombination, Pyrimidine metabolism, Human T-cell leukemia virus one infection, spliceosome, adherens junction, ribosome, N-Glycan biosynthesis, ECM-receptor interaction, ubiquitin mediated proteolysis, and RNA transport pathways. We also performed a GSEA analysis to identify the potential pathways related to SEC61G. The most significantly enriched pathways were E2F targets, G2M checkpoint, glycolysis, MTORC1 signaling, MYC targets, oxidative phosphorylation, DNA repair, hypoxia, unfolded protein response. (Fig. 5) The ER-stress related genes with statistical differences in high and low SEC61G expression groups were shown in additional file 3.

The correlation between SEC61G expression and immune cell infiltration
We further conducted ssGSEA to identify the potential relationship between the expression of SEC61G and immune cell infiltration. The correlation between immune cell infiltration and SEC61G expression is shown in Table 3 and Fig. 6A. The results showed that SEC61G expression positively correlated with the infiltration of Type 2    (Fig. 6B) What's more, The IMvigor 210 study showed that metastatic urothelial cancer patients with higher SEC61G expression level had a better survival benefit from the atezolizumab immunotherapy than those patients with low SEC61G expression (p = 0.006, breslow test) (Fig. 2C).

Knockdown of SEC61G in LAC cell lines
The SEC61G expression was silenced by three small interfering RNAs; as a result, the qRT-PCR assay showed (See figure on previous page.) Fig. 3 The copy number variation (CNV) and methylation of SEC61G in LUAD. A the expression level in different CNV of SEC61G; B the correlation between SEC61G methylation and its expression level; C The methylation levels of SEC61G in pan-cancer and normal tissues from TCGA data; D The promoter methylation of SEC61G in tumor tissues from TCGA-LUAD data; E the Kaplan-Meier survival of the promoter methylation of SEC61G in LUAD Fig. 4 Functional enrichment of SEC61G co-expressed genes in LUAD. A Gene Ontology (GO) enrichment analysis; B Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis that si-1 had the best silencing effect among the three siRNAs against SEC61G in both A549 and H1299 cell lines (Fig. 7A). Then in western blotting assay, We found that the protein expression of SEC61G decreased significantly in the two cell lines transfected with si-1 compared to that in the same cell lines transfected with si-2, si-3 and NC (Fig. 7B). Taken together, these results demonstrate that si-1 was effective in inhibiting the expression of SEC61G, thereby si-1 was selected in subsequent experiments.

Knockdown of SEC61G inhibited the malignant biological behaviors of LAC cells
To explore the regulatory effects of SEC61G on LAC cells' proliferation, apoptosis and invasion, a series of In vitro experiments were conducted. CCK-8 assays evinced that compared to the NC group, knockdown of SEC61G could significantly inhibit the multiplication of A549 and H1299 cells (Fig. 7C-D). Colony formation assay also showed that the number of clones of A549 cells and H1299 cells transfected with si-SEC61G was significantly less than that of the NC group ( Fig. 7E-F). Transwell assay suggested that in comparison to those in the NC group, cell invasion capacities were significantly reduced in the si-SEC61G groups ( Fig. 8A: t = 16.91, P < 0.0001; Fig. 8B: t = 8.469, P = 0.0011). PI-FITC-annexin assay demonstrated that knockdown of SEC61G significantly induced cell apoptosis, compared with the NC group. (Fig. 9A: t = 6.291, P = 0.0033; Fig.  9B: t = 4.901, P = 0.0080). These findings indicated that knockdown of SEC61G promoted apoptosis in lung adenocarcinoma cells, which could partly responsible for SEC61G depletion-induced cell proliferation suppression. The above evidence confirmed that knockdown of SEC61G could remarkably inhibit the malignant phenotypes of LAC cells.

Discussion
Lung adenocarcinoma, which constitutes almost 50% of NSCLCs, is a severe global public health problem with high mortality and morbidity in cancer patients [1,2]. Although newly molecular targeted therapy and immunotherapy have shed lights on NSCLC treatment, the therapeutic efficacy is still limited due to the high heterogeneity of lung cancer and lung cancer patients' outcome is still far from satisfactory [26]. This study aimed to identify a novel prognostic biomarker for LUAD patients and help to guide individualized treatments and optimize the therapeutic strategy. Several studies have demonstrated that SEC61G was overexpressed in various malignancies. However, the potential correlation between SEC61G and lung adenocarcinoma remains unclear. Our study showed that SEC61G was upregulated in LUAD tissues compared with adjacent normal tissues in the TCGA-LUAD cohort, which was consistent with previous studies founded in glioblastoma and hepatocellular carcinoma [14,17]. Our analysis also found that SEC61G is significantly overexpressed in multiple malignancies in the TCGA database, such as gastric, liver cancer and breast carcinomas. These findings indicated that SEC61G could be a potential diagnostic marker in some type of cancers. In our study, the SEC61G gene was demonstrated as a potential diagnostic marker in LUAD, and its AUC exceeded 0.85. In the meanwhile, we also found that SEC61G correlated with the T stage, N stage and pathological grade of LUAD, further corroborating that the expression of SEC61G might be associated with the degree of malignancy of LUAD.
Then we demonstrated that patients with high SEC61G expression were significantly related to shorter overall survival (OS). Multivariate Cox regression analyses revealed that SEC61G expression served as an independent prognostic factor of OS. This result was also validated in LUAD patients from GEO database. Considering that SEC61G is a strong prognostic factor, we constructed a nomogram based on age, pathological stage and the SEC61G expression to predict the 1-, 3-and 5year survival probability in LUAD patients. Calibration curves suggested that the actual prognosis was closely corresponded to the 1-and 3-year predicted prognosis, indicating good prediction performance of the nomogram. The nomogram could help to identify the highrisk patients and choose the more aggressive therapeutic strategy, which had better predictability than individual prognostic factors.
Although various mechanisms can contribute to elevated gene expression levels, DNA methylation and CNVs are two of the most common situations. DNA methylation can not only regulate gene expression but also plays a key role in tumorigenesis [27]. Previous studies demonstrated that several highly-expressed genes due to hypomethylation were associated with poor prognosis in NSCLC [28,29]. Our results showed that SEC61G expression was associated with SEC61G hypomethylation (R = − 0.290, p < 0.001) and the promoter methylation level of SEC61G in LUAD is lower than that in normal tissue. More importantly, SEC61G methylation was related to the prognosis of LUAD, and hypomethylated patients have worse OS, which is in line with the prognostic value of the mRNA expression of SEC61G. Hence, SEC61G hypomethylation might partly contribute to SEC61G overexpression and was associated with poor prognosis in LUAD.
SEC61G encodes the core subunit of SEC61 complex in the endoplasmic reticulum membrane and plays a critical role in protein translocation and cellular calcium homeostasis [8]. Currently, the biological mechanism of SEC61G in tumors is still under exploration. In this study, GO and KEGG analysis indicated that SEC61G was significantly associated with the proliferationassociated biological process such as DNA replication, cell cycle and cell division, which is consistent with previous research in hepatocellular carcinoma [14]. To further verify the biological mechanism of SEC61G, we conducted a series of In vitro experiments in A427 and H1299 cells. Our experiments here demonstrated that knockdown of SEC61G inhibits LAC cell proliferation, invasion, and favorited apoptosis, which is in accordance with our bioinformatic prediction. It is worth noting that GSEA analysis found SEC61G was significantly related to the E2F targets pathway. E2Fs are a complex family of transcriptional regulators and play a key role in protecting cells from cell cycle-generated genomic errors and abnormal proliferation [30]. Li et al. found that overexpressed E2F genes were associated with poor prognosis in lung cancer patients [31]. Therefore, SEC61G might participate in the E2F-related pathway to regulate the cell cycle of lung cancer cells. Based on the C-D CCK8 assay demonstrated the proliferation capacity of A549 and H1299 cells in si-1 group was significantly weaker than that of NC group. * P < 0.05 vs control. E-F Colony formation assay indicated that the clone number of A549 and H1299 cells transfected with si-1 was significantly less than that of NC group above findings, overexpression of SEC61G might take an active part in cell cycle, cell division and E2F-related pathway in LUAD and leading to the occurrence and progression of lung cancer. Another important finding was that SEC61G was significantly related to the unfolded protein response (UPR), which is an adaptive mechanism to reinstate ER proteostasis under the ER stress. Potent ER stress responses have been reported in the majority of human cancer, including brain, lung, breast, colon, gastric, pancreatic, prostate skin. Productive, non-lethal ER stress or activation of UPR-mediated cytoprotective functions, could promote tumorigenesis, mediate resistance to treatment and orchestrate various immune-evasive mechanisms [32]. UPR modulators such as IRE1α kinase inhibitors, PERK inhibitors and eIF2α inhibitors have shown notable anti-tumor efficacy in preclinical cancer models and hence provide a new insight of targets therapies [33][34][35]. It's worth noting that SEC61G might confers a selective growth advantage by facilitating a cytoprotective response to ER stress [17]. Given that Our study showed significant association between SEC61G expression and ER stress response, we suspected that high expression of SEC61G play a critical role in malignant tumor behavior of LUAD via adaptive ER stress responses and SEC61G-related inhibitors might encouraging for the potential treatment of LUAD. Further studies are needed to explore the detailed molecular mechanisms through which ER stress response pathways stimulates the expression of SEC61G and induce its prosurvival effects. With the advent of immunotherapy, the tumor microenvironment has received more and more attention. Studies have shown that immune cells account for a large proportion of the tumor microenvironment and play a critical role in tumor development [25]. Our study demonstrated that the SEC61G expression was negatively correlated with CD4+ T cells, CD8+ T cells, NK cells, DCs and B cells, indicating that SEC61G might play an inhibitory role in both innate immunity and adaptive immunity [36]. There is accumulating evidence indicating that high levels of CD4+ T cells, CD8+ T cells, NK cells, DCs, and B cells infiltration are correlated with better prognosis in LUAD patients [37][38][39]. Additionally, patients with higher levels of CD8+ T cells and CD4+ cells infiltration were more likely to benefit from immunotherapy [40,41]. Interestingly, We found that metastatic urothelial cancer patients with higher SEC61G expression level had a worse prognosis than those patients with low SEC61G expression who both received the atezolizumab immunotherapy in The Invasion assay demonstrated that in comparison to those in NC group, cell invasion capacities were significantly reduced in the si-SEC61G groups. Data are expressed as mean ± standard deviation (SD). ** p < 0.01 and *** p < 0.001 compared with NC group IMvigor 210 study. Hence, it is reasonable to speculate that the SEC61G-mediatedimmunosuppression was one of the underlying causes of poor outcomes in LUAD patients and could provide a reference for the efficacy of immunotherapy in LUAD patients.
Although this current study enhanced a better understanding of the relationship between SEC61G and LUAD, some limitations of our study needed to be noted. First, our investigations into the role of SEC61G in tumors were based on the TIMER, TCGA and GEO databases, which lacks verification from our own clinical samples. Second, given that our study design is limited, additional key signaling pathways associated with SEC61G might be missed, and specific details on these relevant pathways and SEC61G-mediated immunosuppression are still unclear. Traditional in-house experimental studies and prospective studies are needed for further validation.

Conclusion
In summary, using bioinformatic analysis, we systematically analyzed the expression pattern and prognostic value of SEC61G in LUAD patients from various databases. Our results indicated that high SEC61G expression correlated with worse prognosis and SEC61G was an independent prognostic factor for overall survival for LUAD patients. Additionally, vitro experiments verified the biological behaviors of SEC61G in lung cancer. Large-scale and comprehensive researches are needed to strengthen the findings before promoting the clinical efficacy of SEC61G as a prognostic biomarker and therapeutic target. Fig. 9 Knockdown of SEC61G enhanced cell apoptosis in lung adenocarcinoma cells. A-B PI-FITC-annexin assay demonstrated that knockdown of SEC61G significantly induced cell apoptosis. The statistic of apoptotic cells percentage in each group were quantified. Data are expressed as mean ± SD (n = 3). **p < 0.01; ****p < 0.0001