Home Life Sciences Identification of lysosome-related genes in connection with prognosis and immune cell infiltration for drug candidates in head and neck cancer
Article Open Access

Identification of lysosome-related genes in connection with prognosis and immune cell infiltration for drug candidates in head and neck cancer

  • Ye Shen , Haibin Chen , Xiaoyang Gong , Ziyi Wang , Mengjie Chen EMAIL logo and Donghui Chen EMAIL logo
Published/Copyright: August 9, 2023

Abstract

Lysosome dysfunction has been shown to play an important role in cancer progression. However, few research studies have reported the role of lysosomes in head and neck squamous cell carcinoma (HNSCC) progression. Lysosome-related genes (LRGs) were collected from the Molecular Signatures Database. Differentially expressed lysosome-related genes (DELRGs) were identified from the TCGA-HNSCC dataset. The least absolute shrinkage and selection operator and multivariate Cox regression analysis were used to identify the prognostic genes. The prognostic values and expression of hub DELRGs were further validated by GEO datasets. Estimation of STromal and Immune cells in MAlignant Tumors using Expression data and the single-sample gene set enrichment analysis were applied to evaluate the correlation between cathepsin G (CTSG) and immune infiltrates. Twenty-two DELRGs were identified. Among them, CTSG was an independent prognostic biomarker for HNSCC patients. Gene set enrichment analysis indicated that the potential mechanism of CTSG in regulating HNSCC was associated with the immune- and inflammation-related pathways. CTSG expression was highly correlated with immune cell infiltration. Finally, two potential compounds (CH and MAN) targeting CTSG protein were identified, and their reliability was validated through molecular docking analysis. CTSG was associated with immune infiltration and had prognostic value in HNSCC patients, which may be a potential biomarker for predicting the outcome of immunotherapy.

Graphical Abstract

1 Introduction

Head and neck cancer is the world’s sixth leading cause of cancer [1]. It accounts for 7% of all cancers, with approximately 344,000 deaths and 693,000 new cases reported worldwide each year [1]. The majority of head and neck cancers arise in the larynx, throat, and mucous membrane epithelium of the mouth and are known as head and neck squamous cell carcinoma (HNSC) [2]. Despite the considerable advances made in the treatment of malignancies, there has been no significant improvement in head and neck cancer mortality rates in recent decades [3]. In addition, the 5-year survival rate for those diagnosed with HNSC is only about 50%. Immunotherapy has gained much attention in recent years; however, it is only effective in between 20 and 30% of patients [4]. In view of the increase in incidence and the low survival rates of HNSC, and to better tailor diagnostic, preventive, and therapeutic approaches for HNSC, there is an urgent need to develop biomarkers that can accurately predict tumor prognosis [5,6].

Lysosomes play a vital role in cellular aging, development and homeostasis as signaling hubs, and degradation centers within cells [7]. Moreover, complex transcriptional programs regulated lysosomal composition and biogenesis, tailoring their activity to changing cellular needs [8]. It is not surprising that alterations in these fundamental activities are critical to the pathophysiology of diseases, including neurodegenerative disorders, pancreatitis, autoimmune disorders, atherosclerosis, and cancer [9,10,11,12,13]. Thus, unraveling lysosome-related genes (LRGs) will help us to understand the underlying pathological mechanisms and thereby develop new therapeutic approaches relevant to human diseases. Furthermore, lysosomal dysfunction is associated with cancer-associated vascularization, tumor spread, immune escape, cancer cell metastasis, invasion, energy metabolism, and proliferation [14]. Therefore, knowing what influences lysosomes and how they might be regulated will help in cancer diagnosis and prognosis prediction. However, limited studies have reported the role of genes involved in lysosome functions in tumorigenesis and progression.

In the present study, we systematically analyzed LRGs in HNSC. The least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analyses were used to identify the potential prognostic genes for HNSC patients. Furthermore, a series of bioinformatics approaches were used to systematically analyze the prognosis and potential functions of LRGs in HNSC patients.

2 Materials and methods

2.1 Collection of microarray data

The transcriptome and relevant clinical data of HNSC patients (44 normal samples and 504 tumor samples) were obtained from the TCGA dataset. The gene expression profile and survival data from the GSE6631 (22 normal samples and 22 tumor samples) and GSE41613 (97 tumor samples) datasets were downloaded from the GEO database. Lysosome-related data sets were collected from the Molecular Signatures Database. After removing duplicates, 163 LRGs were obtained.

2.2 Identification of the differentially expressed lysosome-related genes (DELRGs) in HNSC patients

DELRGs were filtered between the normal and tumor groups using the “limma” package of R. The conditions for filtering were set to false discovery rate (FDR) <0.05 and |log FC|> 1. The results were visualized using the “ggplot2” package.

2.3 Screening of potential prognostic biomarkers in HNSC patients

The “glmnet” and “survival” packages were used to perform the LASSO analysis based on the above DELRGs. Finally, multivariate Cox regression analysis was performed to further screen signature genes. The relationship between DELRGs and overall survival was assessed and survival curves were plotted using the “ggplot2” and “survminer” packages.

2.4 Construction of nomogram

Based on the results of multivariate analysis, a nomogram was established to predict overall survival in HNSC patients. Nomogram was generated and calibration curves were plotted using the “ggplot2” and “survival” R packages.

2.5 Identification of differentially expressed genes (DEGs) and functional enrichment analysis between the low- and high-cathepsin G (CTSG) HNSC subgroups

Based on the median CTSG expression value, we divided HNSC patients into low- and high-CTSG subgroups. DEGs were identified between the two subgroups using the “limma” package of R. The conditions for the filtering were set to FDR <0.05 and |log FC|> 1. The results were visualized using the “ggplot2” package. To assess the potential biological functions and signaling pathways of DEGs, gene set enrichment analysis (GSEA) was performed using the “clusterProfiler” package of R. The “ggplot2” package was used to visualize enrichment analysis results.

2.6 Immune infiltration analysis

The Estimation of STromal and Immune cells in MAlignant Tumors using Expression data (ESTIMATE) Score, Immune Score, and Stromal Score were compared between the low- and high-CTSG subgroups using the ESTIMATE software. In addition, we used the single-sample gene set enrichment analysis (ssGSEA) to compare the number of 24 immune cell types between the two subgroups. The “ggplot2” package was used to perform the correlation analysis and visualize the results. Tumor Immune Estimation Resource (TIMER) is an online analysis tool for comprehensive analysis of tumor-infiltrating immune cells (https://cistrome.shinyapps.io/timer/) [15]. The TISCH (http://tisch.comp-genomics.org) database is a comprehensive collection of single-cell RNA sequencing data that provide insights into the heterogeneity of the tumor microenvironment across different datasets and cell types [16]. We used this database to explore the heterogeneity of the tumor microenvironment in various cells.

2.7 Molecular docking

First, we screened the candidate drugs targeting CTSG genes by using the Dgibd database (https://dgidb.org/) [17]. The 3D protein structure of CTSG was downloaded from the PDB database (www.rcsb.org). The selection of suitable target proteins was based on the following criteria: (1) the crystal resolution was below 3.5 Å, (2) they must be obtained from Homo sapiens, and (3) the structure with the best resolution. Solvent and organics were removed from the protein conformations using PyMOL software. Subsequently, we performed molecular docking analysis using AutoDock software (1.5.6), and molecular docking results were visualized using PyMOL software.

3 Results

3.1 Identification of DELRGs in HNSC

DELRGs were presented in the volcano plot (Figure 1a). There were 11 downregulated LRGs and 11 upregulated LRGs between the normal and tumor groups. As shown in Figure 1b and c, the mRNA expression of ACP5, ATP6V0D2, CLTCL1, CTSC, CTSV, GALNS, LAMP3, LAPTM4B, SLC11A1, IFI30, and MYO7A was significantly upregulated in HNSC patients, whereas the mRNA expression of ASAH1, ATP6V0A4, CTSG, DNASE2B, SLC11A2, SORT1, MPO, NEU4, TSPAN8, and TYP was significantly downregulated in HNSC patients (p < 0.05).

Figure 1 
                  Identification of DELRGs in HNSC. (a) The volcano plot presents the LRGs. Blue dots indicate downregulated genes and red dots indicate upregulated genes. (b and c) The bar charts present the DELRG expression levels between the normal and tumor groups. *p < 0.05, **p < 0.01, ***p < 0.001, and ns indicates no significant difference.
Figure 1

Identification of DELRGs in HNSC. (a) The volcano plot presents the LRGs. Blue dots indicate downregulated genes and red dots indicate upregulated genes. (b and c) The bar charts present the DELRG expression levels between the normal and tumor groups. *p < 0.05, **p < 0.01, ***p < 0.001, and ns indicates no significant difference.

3.2 Screening of potential prognostic biomarkers in HNSC patients

In the present study, the above 22 DELRGs were used to perform the LASSO regression analysis. Three DELRGs (LAPTM4B, CTSG, and CLTCL1) with non-zero coefficients were screened in HNSC patients (Figure 2a and b). Subsequently, a multivariate regression analysis was performed on these three DELRGs. The results showed that M stage, CTSG, and CLTCL1 were the independent prognostic factors in HNSC patients (Table 1).

Figure 2 
                  Screening of potential prognostic biomarkers for HNSC patients. (a and b) LASSO analysis identifies the prognostic biomarkers in the TCGA dataset. Effect of CLTCL1 genes on overall survival (c), disease-specific survival (d), and progress-free interval (e). Effect of CTSG genes on overall survival (f), disease-specific survival (g), and progress-free interval (h).
Figure 2

Screening of potential prognostic biomarkers for HNSC patients. (a and b) LASSO analysis identifies the prognostic biomarkers in the TCGA dataset. Effect of CLTCL1 genes on overall survival (c), disease-specific survival (d), and progress-free interval (e). Effect of CTSG genes on overall survival (f), disease-specific survival (g), and progress-free interval (h).

Table 1

Results of univariate and multivariate regression analyses

Characteristics Total (N) Univariate analysis Multivariate analysis
Hazard ratio (95% CI) p value Hazard ratio (95% CI) p value
T stage 486
 T1 33 Reference
 T2 143 1.086 (0.568–2.074) 0.803
 T3 131 1.461 (0.769–2.773) 0.247
 T4 179 1.249 (0.665–2.344) 0.490
N stage 479
 N0 238 Reference
 N1 80 1.058 (0.728–1.539) 0.768 1.020 (0.696–1.494) 0.919
 N2&N3 161 1.404 (1.038–1.900) 0.028 1.369 (1.005–1.865) 0.047
M stage 476
 M0 471 Reference
 M1 5 4.745 (1.748–12.883) 0.002 3.945 (1.434–10.852) 0.008
LAPTM4B 501
 Low 251 Reference
 High 250 1.161 (0.888–1.518) 0.275
CTSG 501
 Low 251 Reference
 High 250 0.656 (0.499–0.861) 0.002 0.651 (0.493–0.860) 0.002
CLTCL1 501
 Low 251 Reference
 High 250 1.316 (1.007–1.720) 0.044 1.389 (1.055–1.830) 0.019

Bold values represents a significant difference.

In addition, we performed a comprehensive survival analysis of prognostic genes. As shown in Figure 2c–h, only CLTCL1 genes could predict overall survival (p = 0.04), while CTSG genes could predict overall survival (p = 0.002), disease-specific survival (p = 0.023), and progress-free interval (p = 0.01).

3.3 Validation of expression and prognostic values of CLTCL1 and CTSG

As shown in Figure 3a and b, the gene expression of CLTCL1 and CTSG was significantly downregulated in HNSC patients in the TCGA-HNSC dataset. However, in the GSE6631 dataset, CLTCL1 expression was not significantly different between normal and tumor groups (Figure 3c). CTSG expression was significantly downregulated in the cancer group (p < 0.01, Figure 3d). We also validated the prognostic values of CLTCL1 and CTSG in the GSE41613 dataset. As shown in Figure 3e and f, CLTCL1 genes could not predict overall survival (p = 0.84), while CTSG genes could predict overall survival in HNSC patients (p = 0.0378). These findings further demonstrated that low expression of CTSG was associated with poor prognosis in HNSC patients.

Figure 3 
                  Validation of expression and prognostic values of CLTCL1 and CTSG. The scatter plots present the CLTCL1 (a) and CTSG (b) expression levels between the normal and tumor groups in the TCGA database. The scatter plots present the CLTCL1 (c) and CTSG (d) expression levels between the normal and tumor groups in the GSE6631 database.*p < 0.05, **p < 0.01, ***p < 0.001, and ns indicates no significant difference. Effect of CLTCL1 (e) and CTSG (f) genes on overall survival in the GSE41613 dataset.
Figure 3

Validation of expression and prognostic values of CLTCL1 and CTSG. The scatter plots present the CLTCL1 (a) and CTSG (b) expression levels between the normal and tumor groups in the TCGA database. The scatter plots present the CLTCL1 (c) and CTSG (d) expression levels between the normal and tumor groups in the GSE6631 database.*p < 0.05, **p < 0.01, ***p < 0.001, and ns indicates no significant difference. Effect of CLTCL1 (e) and CTSG (f) genes on overall survival in the GSE41613 dataset.

As shown in Table 2, downregulation of CTSG in HNSC was significantly associated with the primary therapy outcome (p = 0.003), race (p = 0.044), anatomic neoplasm subdivision (p = 0.005), and overall survival (p = 0.015). In addition, as shown in Figure S1, the downregulation of CTSG was associated with a poor prognosis in the subgroup of CR therapy outcome (p = 0.005), radiation therapy (p = 0.008), males (p = 0.021), M0 stage (p = 0.02), N0 stage (p = 0.007), T2 stage (p = 0.041), lymph-node neck dissection (p = 0.002), smokers (p = 0.001), and age >60 (p = 0.002).

Table 2

Relationship of CTSG expression with clinicopathological features in HNSC

Characteristics Low expression of CTSG High expression of CTSG p
n 251 251
T stage, n (%) 0.110
 T1 12 (2.5%) 21 (4.3%)
 T2 69 (14.2%) 75 (15.4%)
 T3 60 (12.3%) 71 (14.6%)
 T4 100 (20.5%) 79 (16.2%)
N stage, n (%) 0.396
 N0 116 (24.2%) 123 (25.6%)
 N1 34 (7.1%) 46 (9.6%)
 N2 83 (17.3%) 71 (14.8%)
 N3 3 (0.6%) 4 (0.8%)
M stage, n (%) 0.682
 M0 232 (48.6%) 240 (50.3%)
 M1 3 (0.6%) 2 (0.4%)
Clinical stage, n (%) 0.268
 Stage I 6 (1.2%) 13 (2.7%)
 Stage II 44 (9%) 51 (10.5%)
 Stage III 49 (10%) 53 (10.9%)
 Stage IV 143 (29.3%) 129 (26.4%)
Radiation therapy, n (%) 0.061
 No 68 (15.4%) 86 (19.5%)
 Yes 155 (35.1%) 132 (29.9%)
Primary therapy outcome, n (%) 0.003
 PD 29 (6.9%) 12 (2.9%)
 SD 1 (0.2%) 5 (1.2%)
 PR 5 (1.2%) 1 (0.2%)
 CR 176 (42.1%) 189 (45.2%)
Gender, n (%) 0.190
 Female 60 (12%) 74 (14.7%)
 Male 191 (38%) 177 (35.3%)
Race, n (%) 0.044
 Asian 3 (0.6%) 7 (1.4%)
 Black or African-American 31 (6.4%) 16 (3.3%)
 White 210 (43.3%) 218 (44.9%)
Age, n (%) 1.000
 ≤60 123 (24.6%) 122 (24.4%)
 >60 128 (25.5%) 128 (25.5%)
Histologic grade, n (%) 0.300
 G1 33 (6.8%) 29 (6%)
 G2 139 (28.8%) 161 (33.3%)
 G3 66 (13.7%) 53 (11%)
 G4 1 (0.2%) 1 (0.2%)
Anatomic neoplasm subdivision, n (%) 0.005
 Alveolar ridge 12 (2.4%) 6 (1.2%)
 Base of tongue 10 (2%) 13 (2.6%)
 Buccal mucosa 12 (2.4%) 10 (2%)
 Floor of mouth 39 (7.8%) 22 (4.4%)
 Hard palate 2 (0.4%) 5 (1%)
 Hypopharynx 8 (1.6%) 2 (0.4%)
 Larynx 53 (10.6%) 58 (11.6%)
 Lip 0 (0%) 3 (0.6%)
 Oral cavity 32 (6.4%) 40 (8%)
 Oral tongue 52 (10.4%) 74 (14.7%)
 Oropharynx 8 (1.6%) 1 (0.2%)
 Tonsil 23 (4.6%) 17 (3.4%)
Smoker, n (%) 1.000
 No 56 (11.4%) 55 (11.2%)
 Yes 190 (38.6%) 191 (38.8%)
Alcohol history, n (%) 0.406
 No 83 (16.9%) 75 (15.3%)
 Yes 160 (32.6%) 173 (35.2%)
Lymphovascular invasion, n (%) 0.379
 No 108 (31.7%) 111 (32.6%)
 Yes 67 (19.6%) 55 (16.1%)
Lymphnode neck dissection, n (%) 0.924
 No 46 (9.2%) 44 (8.8%)
 Yes 204 (40.9%) 205 (41.1%)
OS event, n (%) 0.015
 Alive 128 (25.5%) 156 (31.1%)
 Dead 123 (24.5%) 95 (18.9%)
 Age, median (IQR) 61 (53, 69) 61 (54, 68.75) 0.789

Bold values represents a significant difference.

3.4 Pan-cancer analysis of CTSG expression

As shown in Figure S2, we found that CTSG was significantly downregulated in a variety of cancers, such as bladder urothelial carcinoma, breast invasive carcinoma, cervical squamous cell carcinoma and endocervical adenocarcinoma, cholangio carcinoma, colon adenocarcinoma, esophageal carcinoma, HNSC, kidney chromophobe, liver hepatocellular carcinoma, lung adenocarcinoma, lung squamous cell carcinoma, prostate adenocarcinoma, rectum adenocarcinoma, stomach adenocarcinoma, and endometrial carcinoma. CTSG was significantly upregulated in thyroid carcinoma.

3.5 Construction of a nomogram for HNSC patients

As shown in Figure 4a, we constructed a nomogram and it could predict 1-, 3-, and 5-year overall survival of HNSC patients. In addition, the prediction result of the nomogram was highly consistent with the observation of HNSC patients, based on the nomogram calibration curve (Figure 4b). The diagnostic value of CTSG was assessed by the receiver operating characteristic (ROC) curve, and the area under the curve of CTSG was 0.694 (Figure 4c).

Figure 4 
                  Establishment of a nomogram for HNSC patients. (a) A nomogram for evaluating the survival probability of 1-, 3-, and 5 years for HNSC. (b) Calibration curve of the nomogram. (c) ROC curve of CTSG.
Figure 4

Establishment of a nomogram for HNSC patients. (a) A nomogram for evaluating the survival probability of 1-, 3-, and 5 years for HNSC. (b) Calibration curve of the nomogram. (c) ROC curve of CTSG.

3.6 Screening of CTSG-associated DEGs in HNSC patients

As shown in Figure 5a, a total of 641 CTSG-related DEGs were identified, which included 62 downregulated genes and 579 upregulated genes. The heatmap showed the expression of the top 50 DEGs (Figure 5b). Furthermore, the GSEA was performed to further explore the function of DEGs associated with CTSG expression in HNSC patients. CTSG-related DEGs were involved in immune-related pathways, including neutrophil degranulation (NES = 2.08, p-adj < 0.05), signaling by interleukins (NES = 1.74, p -adj < 0.05), cytokine–cytokine receptor interaction (NES = 2.16, p-adj < 0.05), signaling by B cell receptor BCR (NES = 2.7, p-adj < 0.05), natural killer cell-mediated cytotoxicity (NES = 1.72, p.adj < 0.05), leukocyte transendothelial migration (NES = 2.16, p-adj < 0.05), T-cell receptor signaling pathway (NES = 1.87, p-adj < 0.05), B-cell receptor signaling pathway (NES = 1.92, p-adj < 0.05), acute myeloid leukemia (NES = 1.65, p-adj < 0.05), inflammatory response pathway (NES = 1.92, p-adj < 0.05), cancer immunotherapy by PD1 blockade (NES = 1.86, p-adj < 0.05), etc. (Figure 5c).

Figure 5 
                  Screening of CTSG-associated DEGs in HNSC. (a) The volcano plot presents the DEGs. Blue dots indicate down-regulated genes and red dots indicate up-regulated genes. (b) Co-expression heatmap of top 50 DEGs. (c) Enrichment analysis of DEGs between the low- and high-CTSG subgroups.
Figure 5

Screening of CTSG-associated DEGs in HNSC. (a) The volcano plot presents the DEGs. Blue dots indicate down-regulated genes and red dots indicate up-regulated genes. (b) Co-expression heatmap of top 50 DEGs. (c) Enrichment analysis of DEGs between the low- and high-CTSG subgroups.

3.7 Correlation between the CTSG expression and the infiltration of immune cells

The ESTIMATE algorithm results showed that the low-CTSG subgroup had a lower ESTIMATE score (Figure 6a, p < 0.001), immune score (Figure 6b, p < 0.001), and Sromal score (Figure 6c, p < 0.001) compared to those in the high-CTSG subgroup. In addition, the ssGESA algorithm results indicated that the proportion of dendritic cells (DCs) (p < 0.001), aDC (p < 0.01), B cells (p < 0.001), CD8 T cells (p < 0.05), cytotoxic cells (p < 0.001), eosinophils (p < 0.001), iDC (p < 0.001), macrophages (p < 0.001), mast cells (p < 0.001), neutrophils (p < 0.001), NK CD56dim cells (p < 0.01), NK cells (p < 0.001), pDC (p < 0.001), and T cells (p < 0.001) was significantly decreased in the low-CTSG subgroup compared with those in the high-CTSG subgroup (Figure 6d). The correlation analysis revealed positive correlations between CTSG expression and most immune cells, including mast cells, DC, eosinophils, pDC, B cells, TFH, macrophages, Th17 cells, Th1 cells, T cells, TReg, Tem, cytotoxic cells, NK cells, T helper cells, aDC, NK CD56dim cells, Th2 cells, and CD8 T cells (Figure 6e).

Figure 6 
                  Correlation between the expression of CTSG and the infiltration of immune cells. ESTIMATE algorithm shows the relationship between CTSG expression and ESTIMATE score (a), immune score (b), and Stromal score (c). (d) ssGSEA algorithm showed the relationship between CTSG expression and immune cell infiltration. (e) The correlation results between CTSG genes and immune cell infiltration were visualized by a lollipop plot. *p < 0.05, **p < 0.01, and ***p < 0.001.
Figure 6

Correlation between the expression of CTSG and the infiltration of immune cells. ESTIMATE algorithm shows the relationship between CTSG expression and ESTIMATE score (a), immune score (b), and Stromal score (c). (d) ssGSEA algorithm showed the relationship between CTSG expression and immune cell infiltration. (e) The correlation results between CTSG genes and immune cell infiltration were visualized by a lollipop plot. *p < 0.05, **p < 0.01, and ***p < 0.001.

In addition, we used the TIMER analysis tool to further validate the correlation between CTSG expression and immune cell infiltration. Consistent with the above results, the expression of CTSG was positively correlated with DC (p = 3.02 × 10−8), neutrophil (p = 4.41 × 10−5), macrophage (p = 1.44 × 10−8), CD4+ T cell (p = 2.81 × 10−8), CD8+ T cell (1.91 × 10−6), and B cell (p = 9.5 × 10−5) (Figure 7a). Furthermore, as shown in Figure 7b, the results of copy number alterations in CTSG expression could inhibit immune cell infiltration. We analyzed the level of CTSG genes in tumor immune microenvironment-associated cells, utilizing the TISCH database. Our findings revealed that CTSG showed a higher degree of infiltration in mast cells (Figure S3).

Figure 7 
                  TIMER analysis tool was used to assess the relationship between CTSG and immune cell infiltration. (a) The correlation results between CTSG genes and the proportion of DC, neutrophils, macrophage, CD4+ T cells, CD8+ T cells, and tumor purity. (b) The relationship between somatic copy number alterations of CTSG and immune cell infiltration. *p < 0.05, **p < 0.01, and ***p < 0.001.
Figure 7

TIMER analysis tool was used to assess the relationship between CTSG and immune cell infiltration. (a) The correlation results between CTSG genes and the proportion of DC, neutrophils, macrophage, CD4+ T cells, CD8+ T cells, and tumor purity. (b) The relationship between somatic copy number alterations of CTSG and immune cell infiltration. *p < 0.05, **p < 0.01, and ***p < 0.001.

3.8 Screening of potential drugs and molecular docking

In the present study, the Dgibd database was applied to screen the potential small molecular compounds for CTSG genes. Our findings showed that CHEMBL374027 (CH) and MANNITOL (MAN) were the targeted drugs of CTSG. Subsequently, the molecular docking analysis was performed to further confirm the above results. Our findings showed that CTSG binds to CH with a binding energy of −4.83 kcal/mol and MAN at −7.05 kcal/mol. As presented in Figure 8, there were interactions and binding modes between the CTSG protein and small molecules (CH and MAN). All results showed that CH and MAN had a good potential to target CTSG genes.

Figure 8 
                  Molecular docking analysis. CTSG protein binding to CH (a and b) and MAN (c and d) was visualized with PyMOL software. A 3D surface structure of the CTSG receptor and the small molecular ligand is presented on the left side of the figure (a and c). The binding pattern of the small molecular ligand to the CTSG receptor is presented on the right side of the figure (b and d).
Figure 8

Molecular docking analysis. CTSG protein binding to CH (a and b) and MAN (c and d) was visualized with PyMOL software. A 3D surface structure of the CTSG receptor and the small molecular ligand is presented on the left side of the figure (a and c). The binding pattern of the small molecular ligand to the CTSG receptor is presented on the right side of the figure (b and d).

4 Discussion

Despite advancements in cancer research and treatment, HNSC still has a bleak outlook with little improvement in prognosis over the last few decades [18]. Thus, the identification of new biomarkers and a deeper understanding of the molecular mechanisms involved in HNSC are essential for the development of more effective diagnostic and treatment approaches. In the present study, three LRGs (CLTCL1, LAPTM4B, and CTSG) were identified for predicting the prognosis of HNSC patients. Among these prognostic genes, only the CTSG gene proved to be an independent prognostic factor for HNSC patients. In addition, our findings indicated that low CTSG expression was associated with poor overall survival, disease-specific survival, and progress-free interval in HNSC patients. CTSG is a serine protease mainly found in azurophilic granules of neutrophils [19]. Few studies have shown that CTSG genes were associated with tumor progression and prognosis. However, a recent study has demonstrated that CTSG was a prognostic biomarker for oral squamous cell carcinoma patients [20]. The aim of this study was to examine the biological function and molecular mechanism of CTSG and to determine its potential as a cancer biomarker.

The tumor microenvironment plays a crucial role in the development and progression of cancer as it contains tumor cells that interact with surrounding cells through the lymphatic systems. This interaction influences tumorigenesis and highlights the importance of studying the tumor microenvironment in cancer research [21]. The progression of cancer is heavily influenced by changes in the tumor microenvironment, which includes non-immune and immune cells as well as the extracellular matrix. These changes can have a significant impact on tumor progression [22,23]. In the present study, we also performed DEG analysis and functional enrichment analysis to assess the potential function of CTSG expression in subgroups of HNSC patients. The GSEA results showed CTSG-associated DEGs were mainly implicated in inflammation- and immune-related pathways, such as neutrophil degranulation, signaling by interleukins, natural killer cell-mediated cytotoxicity, B-cell receptor signaling pathway, inflammatory response pathway, etc. Previous studies have revealed that CTSG could activate matrix metalloproteinases and regulate inflammation, and play an important role in a variety of inflammatory diseases such as acute respiratory distress syndrome, rheumatoid arthritis, chronic obstructive pulmonary disease, and periodontitis [19,24,25]. In addition, CTSG genes were shown to be important for maintaining neutrophil-medicated acute histopathology and fibrosis after renal ischemia/reperfusion injury [26]. CTSG could regulate the level of soluble forms of TNF-α and promote the pro-inflammatory effect of neutrophils [27]. Moreover, immune inflammation played a primary role in drug resistance, clinical outcomes, and tumor microenvironment in HNSC [28]. Specifically, interleukin 2 immunotherapy showed promising results in the treatment of HNSC [29,30].

Another important finding was that CTSG expression was positively correlated with the level of infiltration of immune cells (mast cells, DC, eosinophils, neutrophils, B cells, macrophages, CD4 T cells, and CD8 T cells) in HNSC patients. Our findings imply that CTSG plays a significant role in promoting tumor progression by facilitating communication between tumor cells and their inflammatory microenvironment. Increasing evidence suggests that immune cells and their biomarkers can act as therapeutic targets or prognostic markers or participate in tumor progression. An immune-related gene prognostic index differentiated prognosis and served as a potential biomarker in HNSC patients [31]. A recent study revealed the importance of B cells in prognosis and responses to PD-1 blockade in HPV-associated HNSC patients [32]. CD8 T cells were independent prognostic biomarkers in HNSC patients, and patients with low levels of B cells, CD8 T cells, and endothelial cells had a worse prognosis [33]. The poor prognosis of HNSC patients was associated with the presence of tumor-associated macrophages [34]. Neutrophils increased the invasiveness of oral squamous cell carcinoma by activating invadopodia [35]. All these pieces of evidence further showed that CTSG genes impacted the prognosis of HNSC patients by affecting tumor immune cell infiltration. However, the relationship between immunoregulatory cells and CTSG expression requires further investigation.

In our study, two potential drugs (CH and MAN) with high affinity to CTSG protein were identified. Among them, MAN exerted the nephroprotective effect in HNSC patients with cisplatin therapy [36]. Although these potential compounds should be investigated in depth for specific mechanisms, our findings implied their potential in tumor immunotherapy.

5 Conclusion

Our findings suggested that CTSG was an independent prognostic biomarker in HNSC patients and was associated with tumor immune cell infiltration in patients. These results laid a theoretical foundation for further exploration of the role of CTSG in HNSC patients.


# Ye Shen and Haibin Chen contributed equally to this work.


  1. Funding information: Authors state no funding involved.

  2. Author contributions: Ye Shen and Haibin Chen wrote the main manuscript. Xiaoyang Gong and Ziyi Wang prepared figures and performed bioinformatics analysis. Mengjie Chen and Donghui Chen edited the article. All authors reviewed the manuscript.

  3. Conflict of interest: Authors state no conflict of interest.

  4. Data availability statement: The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

References

[1] Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.10.3322/caac.21660Search in Google Scholar PubMed

[2] Johnson DE, Burtness B, Leemans CR, Lui VWY, Bauman JE, Grandis JR. Head and neck squamous cell carcinoma. Nat Rev Dis Primers. 2020;6(1):92.10.1038/s41572-020-00224-3Search in Google Scholar PubMed PubMed Central

[3] Svider PF, Blasco MA, Raza SN, Shkoukani M, Sukari A, Yoo GH, et al. Head and Neck Cancer. Otolaryngol–Head Neck Surg Off J Am Acad Otolaryngol–Head Neck Surg. 2017;156(1):10–3.10.1177/0194599816674672Search in Google Scholar PubMed

[4] Brahmer JR, Drake CG, Wollner I, Powderly JD, Picus J, Sharfman WH, et al. Phase I study of single-agent anti–programmed death-1 (MDX-1106) in refractory solid tumors: Safety, clinical activity, pharmacodynamics, and immunologic correlates. J Clin Oncol. 2023;41(4):715–23.10.1200/JCO.22.02270Search in Google Scholar PubMed

[5] Russell S, Angell T, Lechner M, Liebertz D, Correa A, Sinha U, et al. Immune cell infiltration patterns and survival in head and neck squamous cell carcinoma. Head Neck Oncol. 2013;5(3):24.Search in Google Scholar

[6] Leemans CR, Snijders PJF, Brakenhoff RH. The molecular landscape of head and neck cancer. Nat Rev Cancer. 2018;18(5):269–82.10.1038/nrc.2018.11Search in Google Scholar PubMed

[7] Yang C, Wang X. Lysosome biogenesis: Regulation and functions. The. J Cell Biol. 2021;220(6):e202102001.10.1083/jcb.202102001Search in Google Scholar PubMed PubMed Central

[8] Perera RM, Zoncu R. The lysosome as a regulatory hub. Annu Rev Cell Dev Biol. 2016;32:223–53.10.1146/annurev-cellbio-111315-125125Search in Google Scholar PubMed PubMed Central

[9] Cisneros J, Belton TB, Shum GC, Molakal CG, Wong YC. Mitochondria-lysosome contact site dynamics and misregulation in neurodegenerative diseases. Trends Neurosci. 2022;45(4):312–22.10.1016/j.tins.2022.01.005Search in Google Scholar PubMed PubMed Central

[10] Root J, Merino P, Nuckols A, Johnson M, Kukar T. Lysosome dysfunction as a cause of neurodegenerative diseases: Lessons from frontotemporal dementia and amyotrophic lateral sclerosis. Neurobiol Dis. 2021;154:105360.10.1016/j.nbd.2021.105360Search in Google Scholar PubMed PubMed Central

[11] Radisavljevic Z. Lysosome activates AKT inducing cancer and metastasis. J Cell Biochem. 2019;120(8):12123–7.10.1002/jcb.28752Search in Google Scholar PubMed

[12] Wang S, Ni HM, Chao X, Wang H, Bridges B, Kumer S, et al. Impaired TFEB-mediated lysosomal biogenesis promotes the development of pancreatitis in mice and is associated with human pancreatitis. Autophagy. 2019;15(11):1954–69.10.1080/15548627.2019.1596486Search in Google Scholar PubMed PubMed Central

[13] Ge W, Li D, Gao Y, Cao X. The roles of lysosomes in inflammation and autoimmune diseases. Int Rev Immunol. 2015;34(5):415–31.10.3109/08830185.2014.936587Search in Google Scholar PubMed

[14] Tang T, Yang ZY, Wang D, Yang XY, Wang J, Li L, et al. The role of lysosomes in cancer development and progression. Cell Biosci. 2020;10(1):131.10.1186/s13578-020-00489-xSearch in Google Scholar PubMed PubMed Central

[15] Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: A web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108.10.1158/0008-5472.CAN-17-0307Search in Google Scholar PubMed PubMed Central

[16] Sun D, Wang J, Han Y, Dong X, Ge J, Zheng R, et al. TISCH: A comprehensive web resource enabling interactive single-cell transcriptome visualization of tumor microenvironment. Nucleic Acids Res. 2021;49(D1):D1420–30.10.1093/nar/gkaa1020Search in Google Scholar PubMed PubMed Central

[17] Freshour SL, Kiwala S, Cotto KC, Coffman AC, McMichael JF, Song JJ, et al. Integration of the drug-gene interaction database (DGIdb 4.0) with open crowdsource efforts. Nucleic Acids Res. 2021;49(D1):D1144–51.10.1093/nar/gkaa1084Search in Google Scholar PubMed PubMed Central

[18] Bai S, Zhang P, Zhang JC, Shen J, Xiang X, Yan YB, et al. A gene signature associated with prognosis and immune processes in head and neck squamous cell carcinoma. Head Neck. 2019;41(8):2581–90.10.1002/hed.25731Search in Google Scholar PubMed

[19] Pérez-Is L, Ocaña MG, Montes AH, Carton JA, Álvarez V, Meana Á, et al. The N125S polymorphism in the cathepsin G gene (rs45567233) is associated with susceptibility to osteomyelitis in a Spanish population. PLoS One. 2019;14(10):e0220022.10.1371/journal.pone.0220022Search in Google Scholar PubMed PubMed Central

[20] Huang GZ, Wu QQ, Zheng ZN, Shao TR, Li F, Lu XY, et al. Bioinformatics analyses indicate that Cathepsin G (CTSG) is a potential immune-related biomarker in oral squamous cell carcinoma (OSCC). OncoTargets Ther. 2021;14:1275–89.10.2147/OTT.S293148Search in Google Scholar PubMed PubMed Central

[21] Arneth B. Tumor microenvironment. Medicina (Kaunas, Lithuania). 2019;56(1):15.10.3390/medicina56010015Search in Google Scholar PubMed PubMed Central

[22] Swartz MA, Iida N, Roberts EW, Sangaletti S, Wong MH, Yull FE, et al. Tumor microenvironment complexity: Emerging roles in cancer therapy. Cancer Res. 2012;72(10):2473–80.10.1158/0008-5472.CAN-12-0122Search in Google Scholar PubMed PubMed Central

[23] Junttila MR, de Sauvage FJ. Influence of tumour micro-environment heterogeneity on therapeutic response. Nature. 2013;501(7467):346–54.10.1038/nature12626Search in Google Scholar PubMed

[24] Korkmaz B, Moreau T, Gauthier F. Neutrophil elastase, proteinase 3 and cathepsin G: physicochemical properties, activity and physiopathological functions. Biochimie. 2008;90(2):227–42.10.1016/j.biochi.2007.10.009Search in Google Scholar PubMed

[25] Korkmaz B, Horwitz MS, Jenne DE, Gauthier F. Neutrophil elastase, proteinase 3, and cathepsin G as therapeutic targets in human diseases. Pharmacol Rev. 2010;62(4):726–59.10.1124/pr.110.002733Search in Google Scholar PubMed PubMed Central

[26] Shimoda N, Fukazawa N, Nonomura K, Fairchild RL. Cathepsin g is required for sustained inflammation and tissue injury after reperfusion of ischemic kidneys. Am J Pathol. 2007;170(3):930–40.10.2353/ajpath.2007.060486Search in Google Scholar PubMed PubMed Central

[27] Mezyk-Kopeć R, Bzowska M, Bzowska M, Mickowska B, Mak P, Potempa J, et al. Effects of elastase and cathepsin G on the levels of membrane and soluble TNFalpha. Biol Chem. 2005;386(8):801–11.10.1515/BC.2005.094Search in Google Scholar PubMed

[28] Zhu L, Wang Y, Yuan X, Ma Y, Zhang T, Zhou F, et al. Effects of immune inflammation in head and neck squamous cell carcinoma: Tumor microenvironment, drug resistance, and clinical outcomes. Front Genet. 2022;13:1085700.10.3389/fgene.2022.1085700Search in Google Scholar PubMed PubMed Central

[29] O’Malley BW Jr, Sewell DA, Li D, Kosai K, Chen SH, Woo SL, et al. The role of interleukin-2 in combination adenovirus gene therapy for head and neck cancer. Mol Endocrinol (Baltimore, MD). 1997;11(6):667–73.10.1210/mend.11.6.0012Search in Google Scholar PubMed

[30] O’Malley BW Jr, Li D, McQuone SJ, Ralston R. Combination nonviral interleukin-2 gene immunotherapy for head and neck cancer: From bench top to bedside. Laryngoscope. 2005;115(3):391–404.10.1097/00005537-200503000-00002Search in Google Scholar PubMed

[31] Chen Y, Li ZY, Zhou GQ, Sun Y. An immune-related gene prognostic index for head and neck squamous cell carcinoma. Clin Cancer Res Off J Am Assoc Cancer Res. 2021;27(1):330–41.10.1158/1078-0432.CCR-20-2166Search in Google Scholar PubMed

[32] Kim SS, Shen S, Miyauchi S, Sanders PD, Franiak-Pietryga I, Mell L, et al. B cells improve overall survival in HPV-associated squamous cell carcinomas and are activated by radiation and PD-1 blockade. Clin Cancer Res Off J Am Assoc Cancer Res. 2020;26(13):3345–59.10.1158/1078-0432.CCR-19-3211Search in Google Scholar PubMed PubMed Central

[33] Zhang Y, Li L, Zheng W, Zhang L, Yao N. CD8(+) T-cell exhaustion in the tumor microenvironment of head and neck squamous cell carcinoma determines poor prognosis. Ann Transl Med. 2022;10(6):273.10.21037/atm-22-867Search in Google Scholar PubMed PubMed Central

[34] Li B, Ren M, Zhou X, Han Q, Cheng L. Targeting tumor-associated macrophages in head and neck squamous cell carcinoma. Oral Oncol. 2020;106:104723.10.1016/j.oraloncology.2020.104723Search in Google Scholar PubMed

[35] Glogauer JE, Sun CX, Bradley G, Magalhaes MA. Neutrophils increase oral squamous cell carcinoma invasion through an invadopodia-dependent pathway. Cancer Immunol Res. 2015;3(11):1218–26.10.1158/2326-6066.CIR-15-0017Search in Google Scholar PubMed

[36] Hägerström E, Lindberg L, Bentzen J, Brødbæk K, Zerahn B, Kristensen B. The nephroprotective effect of mannitol in head and neck cancer patients receiving cisplatin therapy. Clin Med Insights Oncol. 2019;13:1179554918821320.10.1177/1179554918821320Search in Google Scholar PubMed PubMed Central

Received: 2023-02-19
Revised: 2023-05-28
Accepted: 2023-06-23
Published Online: 2023-08-09

© 2023 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Biomedical Sciences
  2. Systemic investigation of inetetamab in combination with small molecules to treat HER2-overexpressing breast and gastric cancers
  3. Immunosuppressive treatment for idiopathic membranous nephropathy: An updated network meta-analysis
  4. Identifying two pathogenic variants in a patient with pigmented paravenous retinochoroidal atrophy
  5. Effects of phytoestrogens combined with cold stress on sperm parameters and testicular proteomics in rats
  6. A case of pulmonary embolism with bad warfarin anticoagulant effects caused by E. coli infection
  7. Neutrophilia with subclinical Cushing’s disease: A case report and literature review
  8. Isoimperatorin alleviates lipopolysaccharide-induced periodontitis by downregulating ERK1/2 and NF-κB pathways
  9. Immunoregulation of synovial macrophages for the treatment of osteoarthritis
  10. Novel CPLANE1 c.8948dupT (p.P2984Tfs*7) variant in a child patient with Joubert syndrome
  11. Antiphospholipid antibodies and the risk of thrombosis in myeloproliferative neoplasms
  12. Immunological responses of septic rats to combination therapy with thymosin α1 and vitamin C
  13. High glucose and high lipid induced mitochondrial dysfunction in JEG-3 cells through oxidative stress
  14. Pharmacological inhibition of the ubiquitin-specific protease 8 effectively suppresses glioblastoma cell growth
  15. Levocarnitine regulates the growth of angiotensin II-induced myocardial fibrosis cells via TIMP-1
  16. Age-related changes in peripheral T-cell subpopulations in elderly individuals: An observational study
  17. Single-cell transcription analysis reveals the tumor origin and heterogeneity of human bilateral renal clear cell carcinoma
  18. Identification of iron metabolism-related genes as diagnostic signatures in sepsis by blood transcriptomic analysis
  19. Long noncoding RNA ACART knockdown decreases 3T3-L1 preadipocyte proliferation and differentiation
  20. Surgery, adjuvant immunotherapy plus chemotherapy and radiotherapy for primary malignant melanoma of the parotid gland (PGMM): A case report
  21. Dosimetry comparison with helical tomotherapy, volumetric modulated arc therapy, and intensity-modulated radiotherapy for grade II gliomas: A single‑institution case series
  22. Soy isoflavone reduces LPS-induced acute lung injury via increasing aquaporin 1 and aquaporin 5 in rats
  23. Refractory hypokalemia with sexual dysplasia and infertility caused by 17α-hydroxylase deficiency and triple X syndrome: A case report
  24. Meta-analysis of cancer risk among end stage renal disease undergoing maintenance dialysis
  25. 6-Phosphogluconate dehydrogenase inhibition arrests growth and induces apoptosis in gastric cancer via AMPK activation and oxidative stress
  26. Experimental study on the optimization of ANM33 release in foam cells
  27. Primary retroperitoneal angiosarcoma: A case report
  28. Metabolomic analysis-identified 2-hydroxybutyric acid might be a key metabolite of severe preeclampsia
  29. Malignant pleural effusion diagnosis and therapy
  30. Effect of spaceflight on the phenotype and proteome of Escherichia coli
  31. Comparison of immunotherapy combined with stereotactic radiotherapy and targeted therapy for patients with brain metastases: A systemic review and meta-analysis
  32. Activation of hypermethylated P2RY1 mitigates gastric cancer by promoting apoptosis and inhibiting proliferation
  33. Association between the VEGFR-2 -604T/C polymorphism (rs2071559) and type 2 diabetic retinopathy
  34. The role of IL-31 and IL-34 in the diagnosis and treatment of chronic periodontitis
  35. Triple-negative mouse breast cancer initiating cells show high expression of beta1 integrin and increased malignant features
  36. mNGS facilitates the accurate diagnosis and antibiotic treatment of suspicious critical CNS infection in real practice: A retrospective study
  37. The apatinib and pemetrexed combination has antitumor and antiangiogenic effects against NSCLC
  38. Radiotherapy for primary thyroid adenoid cystic carcinoma
  39. Design and functional preliminary investigation of recombinant antigen EgG1Y162–EgG1Y162 against Echinococcus granulosus
  40. Effects of losartan in patients with NAFLD: A meta-analysis of randomized controlled trial
  41. Bibliometric analysis of METTL3: Current perspectives, highlights, and trending topics
  42. Performance comparison of three scaling algorithms in NMR-based metabolomics analysis
  43. PI3K/AKT/mTOR pathway and its related molecules participate in PROK1 silence-induced anti-tumor effects on pancreatic cancer
  44. The altered expression of cytoskeletal and synaptic remodeling proteins during epilepsy
  45. Effects of pegylated recombinant human granulocyte colony-stimulating factor on lymphocytes and white blood cells of patients with malignant tumor
  46. Prostatitis as initial manifestation of Chlamydia psittaci pneumonia diagnosed by metagenome next-generation sequencing: A case report
  47. NUDT21 relieves sevoflurane-induced neurological damage in rats by down-regulating LIMK2
  48. Association of interleukin-10 rs1800896, rs1800872, and interleukin-6 rs1800795 polymorphisms with squamous cell carcinoma risk: A meta-analysis
  49. Exosomal HBV-DNA for diagnosis and treatment monitoring of chronic hepatitis B
  50. Shear stress leads to the dysfunction of endothelial cells through the Cav-1-mediated KLF2/eNOS/ERK signaling pathway under physiological conditions
  51. Interaction between the PI3K/AKT pathway and mitochondrial autophagy in macrophages and the leukocyte count in rats with LPS-induced pulmonary infection
  52. Meta-analysis of the rs231775 locus polymorphism in the CTLA-4 gene and the susceptibility to Graves’ disease in children
  53. Cloning, subcellular localization and expression of phosphate transporter gene HvPT6 of hulless barley
  54. Coptisine mitigates diabetic nephropathy via repressing the NRLP3 inflammasome
  55. Significant elevated CXCL14 and decreased IL-39 levels in patients with tuberculosis
  56. Whole-exome sequencing applications in prenatal diagnosis of fetal bowel dilatation
  57. Gemella morbillorum infective endocarditis: A case report and literature review
  58. An unusual ectopic thymoma clonal evolution analysis: A case report
  59. Severe cumulative skin toxicity during toripalimab combined with vemurafenib following toripalimab alone
  60. Detection of V. vulnificus septic shock with ARDS using mNGS
  61. Novel rare genetic variants of familial and sporadic pulmonary atresia identified by whole-exome sequencing
  62. The influence and mechanistic action of sperm DNA fragmentation index on the outcomes of assisted reproduction technology
  63. Novel compound heterozygous mutations in TELO2 in an infant with You-Hoover-Fong syndrome: A case report and literature review
  64. ctDNA as a prognostic biomarker in resectable CLM: Systematic review and meta-analysis
  65. Diagnosis of primary amoebic meningoencephalitis by metagenomic next-generation sequencing: A case report
  66. Phylogenetic analysis of promoter regions of human Dolichol kinase (DOLK) and orthologous genes using bioinformatics tools
  67. Collagen changes in rabbit conjunctiva after conjunctival crosslinking
  68. Effects of NM23 transfection of human gastric carcinoma cells in mice
  69. Oral nifedipine and phytosterol, intravenous nicardipine, and oral nifedipine only: Three-arm, retrospective, cohort study for management of severe preeclampsia
  70. Case report of hepatic retiform hemangioendothelioma: A rare tumor treated with ultrasound-guided microwave ablation
  71. Curcumin induces apoptosis in human hepatocellular carcinoma cells by decreasing the expression of STAT3/VEGF/HIF-1α signaling
  72. Rare presentation of double-clonal Waldenström macroglobulinemia with pulmonary embolism: A case report
  73. Giant duplication of the transverse colon in an adult: A case report and literature review
  74. Ectopic thyroid tissue in the breast: A case report
  75. SDR16C5 promotes proliferation and migration and inhibits apoptosis in pancreatic cancer
  76. Vaginal metastasis from breast cancer: A case report
  77. Screening of the best time window for MSC transplantation to treat acute myocardial infarction with SDF-1α antibody-loaded targeted ultrasonic microbubbles: An in vivo study in miniswine
  78. Inhibition of TAZ impairs the migration ability of melanoma cells
  79. Molecular complexity analysis of the diagnosis of Gitelman syndrome in China
  80. Effects of maternal calcium and protein intake on the development and bone metabolism of offspring mice
  81. Identification of winter wheat pests and diseases based on improved convolutional neural network
  82. Ultra-multiplex PCR technique to guide treatment of Aspergillus-infected aortic valve prostheses
  83. Virtual high-throughput screening: Potential inhibitors targeting aminopeptidase N (CD13) and PIKfyve for SARS-CoV-2
  84. Immune checkpoint inhibitors in cancer patients with COVID-19
  85. Utility of methylene blue mixed with autologous blood in preoperative localization of pulmonary nodules and masses
  86. Integrated analysis of the microbiome and transcriptome in stomach adenocarcinoma
  87. Berberine suppressed sarcopenia insulin resistance through SIRT1-mediated mitophagy
  88. DUSP2 inhibits the progression of lupus nephritis in mice by regulating the STAT3 pathway
  89. Lung abscess by Fusobacterium nucleatum and Streptococcus spp. co-infection by mNGS: A case series
  90. Genetic alterations of KRAS and TP53 in intrahepatic cholangiocarcinoma associated with poor prognosis
  91. Granulomatous polyangiitis involving the fourth ventricle: Report of a rare case and a literature review
  92. Studying infant mortality: A demographic analysis based on data mining models
  93. Metaplastic breast carcinoma with osseous differentiation: A report of a rare case and literature review
  94. Protein Z modulates the metastasis of lung adenocarcinoma cells
  95. Inhibition of pyroptosis and apoptosis by capsaicin protects against LPS-induced acute kidney injury through TRPV1/UCP2 axis in vitro
  96. TAK-242, a toll-like receptor 4 antagonist, against brain injury by alleviates autophagy and inflammation in rats
  97. Primary mediastinum Ewing’s sarcoma with pleural effusion: A case report and literature review
  98. Association of ADRB2 gene polymorphisms and intestinal microbiota in Chinese Han adolescents
  99. Tanshinone IIA alleviates chondrocyte apoptosis and extracellular matrix degeneration by inhibiting ferroptosis
  100. Study on the cytokines related to SARS-Cov-2 in testicular cells and the interaction network between cells based on scRNA-seq data
  101. Effect of periostin on bone metabolic and autophagy factors during tooth eruption in mice
  102. HP1 induces ferroptosis of renal tubular epithelial cells through NRF2 pathway in diabetic nephropathy
  103. Intravaginal estrogen management in postmenopausal patients with vaginal squamous intraepithelial lesions along with CO2 laser ablation: A retrospective study
  104. Hepatocellular carcinoma cell differentiation trajectory predicts immunotherapy, potential therapeutic drugs, and prognosis of patients
  105. Effects of physical exercise on biomarkers of oxidative stress in healthy subjects: A meta-analysis of randomized controlled trials
  106. Identification of lysosome-related genes in connection with prognosis and immune cell infiltration for drug candidates in head and neck cancer
  107. Development of an instrument-free and low-cost ELISA dot-blot test to detect antibodies against SARS-CoV-2
  108. Research progress on gas signal molecular therapy for Parkinson’s disease
  109. Adiponectin inhibits TGF-β1-induced skin fibroblast proliferation and phenotype transformation via the p38 MAPK signaling pathway
  110. The G protein-coupled receptor-related gene signatures for predicting prognosis and immunotherapy response in bladder urothelial carcinoma
  111. α-Fetoprotein contributes to the malignant biological properties of AFP-producing gastric cancer
  112. CXCL12/CXCR4/CXCR7 axis in placenta tissues of patients with placenta previa
  113. Association between thyroid stimulating hormone levels and papillary thyroid cancer risk: A meta-analysis
  114. Significance of sTREM-1 and sST2 combined diagnosis for sepsis detection and prognosis prediction
  115. Diagnostic value of serum neuroactive substances in the acute exacerbation of chronic obstructive pulmonary disease complicated with depression
  116. Research progress of AMP-activated protein kinase and cardiac aging
  117. TRIM29 knockdown prevented the colon cancer progression through decreasing the ubiquitination levels of KRT5
  118. Cross-talk between gut microbiota and liver steatosis: Complications and therapeutic target
  119. Metastasis from small cell lung cancer to ovary: A case report
  120. The early diagnosis and pathogenic mechanisms of sepsis-related acute kidney injury
  121. The effect of NK cell therapy on sepsis secondary to lung cancer: A case report
  122. Erianin alleviates collagen-induced arthritis in mice by inhibiting Th17 cell differentiation
  123. Loss of ACOX1 in clear cell renal cell carcinoma and its correlation with clinical features
  124. Signalling pathways in the osteogenic differentiation of periodontal ligament stem cells
  125. Crosstalk between lactic acid and immune regulation and its value in the diagnosis and treatment of liver failure
  126. Clinicopathological features and differential diagnosis of gastric pleomorphic giant cell carcinoma
  127. Traumatic brain injury and rTMS-ERPs: Case report and literature review
  128. Extracellular fibrin promotes non-small cell lung cancer progression through integrin β1/PTEN/AKT signaling
  129. Knockdown of DLK4 inhibits non-small cell lung cancer tumor growth by downregulating CKS2
  130. The co-expression pattern of VEGFR-2 with indicators related to proliferation, apoptosis, and differentiation of anagen hair follicles
  131. Inflammation-related signaling pathways in tendinopathy
  132. CD4+ T cell count in HIV/TB co-infection and co-occurrence with HL: Case report and literature review
  133. Clinical analysis of severe Chlamydia psittaci pneumonia: Case series study
  134. Bioinformatics analysis to identify potential biomarkers for the pulmonary artery hypertension associated with the basement membrane
  135. Influence of MTHFR polymorphism, alone or in combination with smoking and alcohol consumption, on cancer susceptibility
  136. Catharanthus roseus (L.) G. Don counteracts the ampicillin resistance in multiple antibiotic-resistant Staphylococcus aureus by downregulation of PBP2a synthesis
  137. Combination of a bronchogenic cyst in the thoracic spinal canal with chronic myelocytic leukemia
  138. Bacterial lipoprotein plays an important role in the macrophage autophagy and apoptosis induced by Salmonella typhimurium and Staphylococcus aureus
  139. TCL1A+ B cells predict prognosis in triple-negative breast cancer through integrative analysis of single-cell and bulk transcriptomic data
  140. Ezrin promotes esophageal squamous cell carcinoma progression via the Hippo signaling pathway
  141. Ferroptosis: A potential target of macrophages in plaque vulnerability
  142. Predicting pediatric Crohn's disease based on six mRNA-constructed risk signature using comprehensive bioinformatic approaches
  143. Applications of genetic code expansion and photosensitive UAAs in studying membrane proteins
  144. HK2 contributes to the proliferation, migration, and invasion of diffuse large B-cell lymphoma cells by enhancing the ERK1/2 signaling pathway
  145. IL-17 in osteoarthritis: A narrative review
  146. Circadian cycle and neuroinflammation
  147. Probiotic management and inflammatory factors as a novel treatment in cirrhosis: A systematic review and meta-analysis
  148. Hemorrhagic meningioma with pulmonary metastasis: Case report and literature review
  149. SPOP regulates the expression profiles and alternative splicing events in human hepatocytes
  150. Knockdown of SETD5 inhibited glycolysis and tumor growth in gastric cancer cells by down-regulating Akt signaling pathway
  151. PTX3 promotes IVIG resistance-induced endothelial injury in Kawasaki disease by regulating the NF-κB pathway
  152. Pancreatic ectopic thyroid tissue: A case report and analysis of literature
  153. The prognostic impact of body mass index on female breast cancer patients in underdeveloped regions of northern China differs by menopause status and tumor molecular subtype
  154. Report on a case of liver-originating malignant melanoma of unknown primary
  155. Case report: Herbal treatment of neutropenic enterocolitis after chemotherapy for breast cancer
  156. The fibroblast growth factor–Klotho axis at molecular level
  157. Characterization of amiodarone action on currents in hERG-T618 gain-of-function mutations
  158. A case report of diagnosis and dynamic monitoring of Listeria monocytogenes meningitis with NGS
  159. Effect of autologous platelet-rich plasma on new bone formation and viability of a Marburg bone graft
  160. Small breast epithelial mucin as a useful prognostic marker for breast cancer patients
  161. Continuous non-adherent culture promotes transdifferentiation of human adipose-derived stem cells into retinal lineage
  162. Nrf3 alleviates oxidative stress and promotes the survival of colon cancer cells by activating AKT/BCL-2 signal pathway
  163. Favorable response to surufatinib in a patient with necrolytic migratory erythema: A case report
  164. Case report of atypical undernutrition of hypoproteinemia type
  165. Down-regulation of COL1A1 inhibits tumor-associated fibroblast activation and mediates matrix remodeling in the tumor microenvironment of breast cancer
  166. Sarcoma protein kinase inhibition alleviates liver fibrosis by promoting hepatic stellate cells ferroptosis
  167. Research progress of serum eosinophil in chronic obstructive pulmonary disease and asthma
  168. Clinicopathological characteristics of co-existing or mixed colorectal cancer and neuroendocrine tumor: Report of five cases
  169. Role of menopausal hormone therapy in the prevention of postmenopausal osteoporosis
  170. Precisional detection of lymph node metastasis using tFCM in colorectal cancer
  171. Advances in diagnosis and treatment of perimenopausal syndrome
  172. A study of forensic genetics: ITO index distribution and kinship judgment between two individuals
  173. Acute lupus pneumonitis resembling miliary tuberculosis: A case-based review
  174. Plasma levels of CD36 and glutathione as biomarkers for ruptured intracranial aneurysm
  175. Fractalkine modulates pulmonary angiogenesis and tube formation by modulating CX3CR1 and growth factors in PVECs
  176. Novel risk prediction models for deep vein thrombosis after thoracotomy and thoracoscopic lung cancer resections, involving coagulation and immune function
  177. Exploring the diagnostic markers of essential tremor: A study based on machine learning algorithms
  178. Evaluation of effects of small-incision approach treatment on proximal tibia fracture by deep learning algorithm-based magnetic resonance imaging
  179. An online diagnosis method for cancer lesions based on intelligent imaging analysis
  180. Medical imaging in rheumatoid arthritis: A review on deep learning approach
  181. Predictive analytics in smart healthcare for child mortality prediction using a machine learning approach
  182. Utility of neutrophil–lymphocyte ratio and platelet–lymphocyte ratio in predicting acute-on-chronic liver failure survival
  183. A biomedical decision support system for meta-analysis of bilateral upper-limb training in stroke patients with hemiplegia
  184. TNF-α and IL-8 levels are positively correlated with hypobaric hypoxic pulmonary hypertension and pulmonary vascular remodeling in rats
  185. Stochastic gradient descent optimisation for convolutional neural network for medical image segmentation
  186. Comparison of the prognostic value of four different critical illness scores in patients with sepsis-induced coagulopathy
  187. Application and teaching of computer molecular simulation embedded technology and artificial intelligence in drug research and development
  188. Hepatobiliary surgery based on intelligent image segmentation technology
  189. Value of brain injury-related indicators based on neural network in the diagnosis of neonatal hypoxic-ischemic encephalopathy
  190. Analysis of early diagnosis methods for asymmetric dementia in brain MR images based on genetic medical technology
  191. Early diagnosis for the onset of peri-implantitis based on artificial neural network
  192. Clinical significance of the detection of serum IgG4 and IgG4/IgG ratio in patients with thyroid-associated ophthalmopathy
  193. Forecast of pain degree of lumbar disc herniation based on back propagation neural network
  194. SPA-UNet: A liver tumor segmentation network based on fused multi-scale features
  195. Systematic evaluation of clinical efficacy of CYP1B1 gene polymorphism in EGFR mutant non-small cell lung cancer observed by medical image
  196. Rehabilitation effect of intelligent rehabilitation training system on hemiplegic limb spasms after stroke
  197. A novel approach for minimising anti-aliasing effects in EEG data acquisition
  198. ErbB4 promotes M2 activation of macrophages in idiopathic pulmonary fibrosis
  199. Clinical role of CYP1B1 gene polymorphism in prediction of postoperative chemotherapy efficacy in NSCLC based on individualized health model
  200. Lung nodule segmentation via semi-residual multi-resolution neural networks
  201. Evaluation of brain nerve function in ICU patients with Delirium by deep learning algorithm-based resting state MRI
  202. A data mining technique for detecting malignant mesothelioma cancer using multiple regression analysis
  203. Markov model combined with MR diffusion tensor imaging for predicting the onset of Alzheimer’s disease
  204. Effectiveness of the treatment of depression associated with cancer and neuroimaging changes in depression-related brain regions in patients treated with the mediator-deuterium acupuncture method
  205. Molecular mechanism of colorectal cancer and screening of molecular markers based on bioinformatics analysis
  206. Monitoring and evaluation of anesthesia depth status data based on neuroscience
  207. Exploring the conformational dynamics and thermodynamics of EGFR S768I and G719X + S768I mutations in non-small cell lung cancer: An in silico approaches
  208. Optimised feature selection-driven convolutional neural network using gray level co-occurrence matrix for detection of cervical cancer
  209. Incidence of different pressure patterns of spinal cerebellar ataxia and analysis of imaging and genetic diagnosis
  210. Pathogenic bacteria and treatment resistance in older cardiovascular disease patients with lung infection and risk prediction model
  211. Adoption value of support vector machine algorithm-based computed tomography imaging in the diagnosis of secondary pulmonary fungal infections in patients with malignant hematological disorders
  212. From slides to insights: Harnessing deep learning for prognostic survival prediction in human colorectal cancer histology
  213. Ecology and Environmental Science
  214. Monitoring of hourly carbon dioxide concentration under different land use types in arid ecosystem
  215. Comparing the differences of prokaryotic microbial community between pit walls and bottom from Chinese liquor revealed by 16S rRNA gene sequencing
  216. Effects of cadmium stress on fruits germination and growth of two herbage species
  217. Bamboo charcoal affects soil properties and bacterial community in tea plantations
  218. Optimization of biogas potential using kinetic models, response surface methodology, and instrumental evidence for biodegradation of tannery fleshings during anaerobic digestion
  219. Understory vegetation diversity patterns of Platycladus orientalis and Pinus elliottii communities in Central and Southern China
  220. Studies on macrofungi diversity and discovery of new species of Abortiporus from Baotianman World Biosphere Reserve
  221. Food Science
  222. Effect of berrycactus fruit (Myrtillocactus geometrizans) on glutamate, glutamine, and GABA levels in the frontal cortex of rats fed with a high-fat diet
  223. Guesstimate of thymoquinone diversity in Nigella sativa L. genotypes and elite varieties collected from Indian states using HPTLC technique
  224. Analysis of bacterial community structure of Fuzhuan tea with different processing techniques
  225. Untargeted metabolomics reveals sour jujube kernel benefiting the nutritional value and flavor of Morchella esculenta
  226. Mycobiota in Slovak wine grapes: A case study from the small Carpathians wine region
  227. Elemental analysis of Fadogia ancylantha leaves used as a nutraceutical in Mashonaland West Province, Zimbabwe
  228. Microbiological transglutaminase: Biotechnological application in the food industry
  229. Influence of solvent-free extraction of fish oil from catfish (Clarias magur) heads using a Taguchi orthogonal array design: A qualitative and quantitative approach
  230. Chromatographic analysis of the chemical composition and anticancer activities of Curcuma longa extract cultivated in Palestine
  231. The potential for the use of leghemoglobin and plant ferritin as sources of iron
  232. Investigating the association between dietary patterns and glycemic control among children and adolescents with T1DM
  233. Bioengineering and Biotechnology
  234. Biocompatibility and osteointegration capability of β-TCP manufactured by stereolithography 3D printing: In vitro study
  235. Clinical characteristics and the prognosis of diabetic foot in Tibet: A single center, retrospective study
  236. Agriculture
  237. Biofertilizer and NPSB fertilizer application effects on nodulation and productivity of common bean (Phaseolus vulgaris L.) at Sodo Zuria, Southern Ethiopia
  238. On correlation between canopy vegetation and growth indexes of maize varieties with different nitrogen efficiencies
  239. Exopolysaccharides from Pseudomonas tolaasii inhibit the growth of Pleurotus ostreatus mycelia
  240. A transcriptomic evaluation of the mechanism of programmed cell death of the replaceable bud in Chinese chestnut
  241. Melatonin enhances salt tolerance in sorghum by modulating photosynthetic performance, osmoregulation, antioxidant defense, and ion homeostasis
  242. Effects of plant density on alfalfa (Medicago sativa L.) seed yield in western Heilongjiang areas
  243. Identification of rice leaf diseases and deficiency disorders using a novel DeepBatch technique
  244. Artificial intelligence and internet of things oriented sustainable precision farming: Towards modern agriculture
  245. Animal Sciences
  246. Effect of ketogenic diet on exercise tolerance and transcriptome of gastrocnemius in mice
  247. Combined analysis of mRNA–miRNA from testis tissue in Tibetan sheep with different FecB genotypes
  248. Isolation, identification, and drug resistance of a partially isolated bacterium from the gill of Siniperca chuatsi
  249. Tracking behavioral changes of confined sows from the first mating to the third parity
  250. The sequencing of the key genes and end products in the TLR4 signaling pathway from the kidney of Rana dybowskii exposed to Aeromonas hydrophila
  251. Development of a new candidate vaccine against piglet diarrhea caused by Escherichia coli
  252. Plant Sciences
  253. Crown and diameter structure of pure Pinus massoniana Lamb. forest in Hunan province, China
  254. Genetic evaluation and germplasm identification analysis on ITS2, trnL-F, and psbA-trnH of alfalfa varieties germplasm resources
  255. Tissue culture and rapid propagation technology for Gentiana rhodantha
  256. Effects of cadmium on the synthesis of active ingredients in Salvia miltiorrhiza
  257. Cloning and expression analysis of VrNAC13 gene in mung bean
  258. Chlorate-induced molecular floral transition revealed by transcriptomes
  259. Effects of warming and drought on growth and development of soybean in Hailun region
  260. Effects of different light conditions on transient expression and biomass in Nicotiana benthamiana leaves
  261. Comparative analysis of the rhizosphere microbiome and medicinally active ingredients of Atractylodes lancea from different geographical origins
  262. Distinguish Dianthus species or varieties based on chloroplast genomes
  263. Comparative transcriptomes reveal molecular mechanisms of apple blossoms of different tolerance genotypes to chilling injury
  264. Study on fresh processing key technology and quality influence of Cut Ophiopogonis Radix based on multi-index evaluation
  265. An advanced approach for fig leaf disease detection and classification: Leveraging image processing and enhanced support vector machine methodology
  266. Erratum
  267. Erratum to “Protein Z modulates the metastasis of lung adenocarcinoma cells”
  268. Erratum to “BRCA1 subcellular localization regulated by PI3K signaling pathway in triple-negative breast cancer MDA-MB-231 cells and hormone-sensitive T47D cells”
  269. Retraction
  270. Retraction to “Protocatechuic acid attenuates cerebral aneurysm formation and progression by inhibiting TNF-alpha/Nrf-2/NF-kB-mediated inflammatory mechanisms in experimental rats”
Downloaded on 27.1.2026 from https://www.degruyterbrill.com/document/doi/10.1515/biol-2022-0660/html
Scroll to top button