Abstract

WD-repeat protein 72, a scaffolding protein lacking intrinsic enzymatic activity, produces numerous β-propeller blade formations, serves as a binding platform to assemble protein complexes and is critical for cell growth, differentiation, adhesion, and migration. Despite evidence supporting a basic role of WDR72 in the tumorigenesis of particular cancers, the value of WDR72 in non-small-cell lung cancer (NSCLC), the tumor with the highest mortality rate globally, is undocumented. We investigated the prognostic value of WDR72 in NSCLC and studied its potential immune function and its correlation with ferroptosis. According to The Cancer Genome Atlas, Cancer Cell Line Encyclopedia, Genotype-Tissue Expression, and Gene Set Cancer Analysis, we used multiple bioinformatic strategies to investigate the possible oncogenic role of WDR72, analyze WDR72 and prognosis, and immune cell infiltration in different tumors correlation. WDR72 exhibited a high expression in NSCLC and a positive association with prognosis. WDR72 expression was related to immune cell infiltration and tumor immune microenvironment in NSCLC. Finally, we validated WDR72 in human NSCLC; it has a predictive value in NSCLC related to its function in tumor progression and immunity. The significance of our study is that WDR72 can be used as a potential indicator of lung cancer prognosis. Helping physicians more accurately predict patient survival and risk of disease progression.

1. Introduction

Lung cancer is a deadly malignancy widespread worldwide, and non-small-cell lung cancer (NSCLC) represents approximately 85% of all lung cancers [1]. NSCLC onset is devious, with invisible clinical manifestations in the early stages, making early diagnosis hard. Furthermore, due to the highly aggressive nature of NSCLC, many patients have distant metastases at diagnosis [2]. Despite recent advances, early detection and treatment of NSCLC remain inadequate [3]. Gene mutation and aberrant expression lead to malignant alteration of airway epithelial cells, eventually leading to lung cancer [4]. Discovering new markers and driver genes is essential for the early detection, diagnosis, targeted treatment, and prognosis assessment of NSCLC.

WD-repeat protein 72 is a scaffold lacking intrinsic enzymatic activity. It produces several β-propeller blade formations, acting as binding platforms to assemble protein complexes and form tooth enamel, as well as causing autosomal recessive developmental defects [5]. Ibrahim-Verbaas et al. [6] identified a significant genome-wide association between a single nucleotide polymorphism (SNP) in the WDR72 gene (chromosome 15) and a cognitive test comparable to the Stroop interference test, demonstrating that at the level of executive function, WDR72 gene is also related to renal function. Howles et al. [7] demonstrated that WDR72 could involve in clathrin-mediated endocytosis, a critical process to persist intracellular CaSR signaling. WDR72 has a significant role in clear renal cell carcinoma [8], esophageal cancer [9], and colorectal cancer [10] but undetected in NSCLC. Accordingly, this paper intended to further analyze the specific mechanism of WDR72 and NSCLC.

Ferroptosis is a sort of controlled necrosis that performs a part in neurodegenerative cancers. Lipid peroxidation inhibitors, lipophilic antioxidants, and iron chelators can prevent ferroptosis [11]. Initially, ferroptosis was only recognized as erastin-induced death; however, additional studies revealed its function in various pathophysiological processes and diseases. Since certain oncogenically mutated tumor cells are more susceptible to ferroptosis, ferroptosis can be triggered, and ferroptosis-sensitive tumor cells may have significant therapeutic potential [1214]. Therefore, modulation of ferroptosis may offer therapeutic potential for certain disorders associated with ferroptosis. Many key targets are important mediators of ferroptosis induction; however, the relationship between WDR72 and ferroptosis is still unclear. Herein, we further study whether WDR72 is sensitive to ferroptosis-induced cell death. Studies have shown that excessive iron overload may suppress the function of the immune system, thereby increasing the risk of infection and tumor development [15]. Immune infiltration refers to the presence of immune cells (T cells, B cells, macrophages, natural killer cells, etc.) in tumor tissues [16]. The immune infiltration of lung cancer patients is closely related to the survival and prognosis of patients [17]. Immune cells can slow the growth and spread of tumors by attacking tumor cells [18]. However, studies have also shown that tumor cells can evade the attack of immune cells through different mechanisms, thereby suppressing the immune response [19]. Overall, immune infiltration and ferroptosis are associated with the prognosis and survival of lung cancer patients, but more studies are needed to confirm the direct relationship between the two.

2. Materials and Methods

2.1. Gene Expression and Clinical Data Collection

First, we downloaded gene expression and clinical transcriptome profiling information for NSCLC patients from The Cancer Genome Atlas (TCGA) website (http://portal.gdc.cancer. gov/). The R tool “limma” was utilized to normalize raw gene expression information. GSE19804 [20, 21], GSE118370 [22], GSE19188 [23], GSE27262 [24, 25], and GSE33532 [26] were obtained from the GEO database GPL570 platform ([HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array) through the “GEOquery” package for subsequent analysis (Table 1).

2.2. Identification DEGS

The selected datasets above were compared and examined utilizing the NCBI (https://www.ncbi.nlm.nih.gov/) toolGEO2R (http://www.ncbi.nlm.nih.gov/geo/geo2r). We applied adjusted (adj.) values and thresholds to the Benjamini and Hochberg false discovery rates (FDR) are utilized to develop a balance between the statistical thresholds of finding significant genes and false positives. Probe sets without matching gene symbols or genes with several probe sets were excluded or averaged accordingly. Log fold change and indicated statistical significance.

2.3. Analysis of Association between WDR72 and Prognosis of Patients with Small-Cell Lung Cancer (SCLC)

In the TCGA database, a Kaplan-Meier (KM) analysis was conducted to examine the overall survival (OS), disease-specific survival (DSS), disease-free survival (DFS), and progression-free survival (PFI) of patients with NSCLC. To show the connection between WDR72 expression and survival in NSCLC, a Cox regression analysis was done utilizing the R tools “survival” and “forest plot.”

2.4. Association between WDR72 Expression and Immunity

Employing the ESTIMATE algorithm, immune and stromal scores of NSCLC samples were measured to investigate the association between WDR72 expression and TME (tumor microenvironment). The relationship between WDR72 expression and score was assessed utilizing the R tools “estimate” and “limma” based on the degree of immune infiltration. We downloaded the TCGA immune cell infiltration value from the CIBERSORT database (https://cibersortx.stanford.edu/) and calculated the relative score of 24 immune cells in NSCLC utilizing the CIBERSORT algorithm. Then, the Spearman rank correlation coefficient was employed to assess the association between the level of WDR72 and the infiltration level of each immune cell in NSCLC. R tools “limma,” “reshape2,” and “RColorBrewer” were utilized to present the visualizations.

2.5. Gene Set Enrichment Analysis

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) gene sets were downloaded from the Gene Set Enrichment Analysis (GSEA) website. The biological role of WDR72 in NSCLC was investigated using the GSEA method. Both analyses were performed using the R “Cluster Profiler” tool.

2.6. Clinical Samples

We should be confirmed non-small-cell lung cancer at Nantong University from 2012 to 2020. A total of 158 matched NSCLC and paracancerous tissues were obtained through radical resection in the Affiliated Hospital of Nantong University. All diagnoses were confirmed by histopathological examination. For further study, the specimens were frozen in a -80°C freezer after surgery. All 158 patients had signed written informed consent and underwent a standardized ethnic review.

2.7. Reverse Transcription-Polymerase Chain Reaction (RT-qPCR)

Four matched lung adenocarcinoma and paracancerous tissues were randomly selected to extract total RNA. RNA was reverse transcript into cDNA utilizing a reverse transcription kit (transgene Biotech, China). NanoDrop 2000 (Thermo Fisher Scientific, Waltham, USA) was utilized to detect its fluorescent expression. GAPDH (glyceraldehyde-3-phosphate dehydrogenase) served as an internal control. WDR72 relative expression was computed employing the 2-ΔΔCT approach.WDR72 F:GCAACTCAAACTCGGCAAACTTCC, R:GGCTCACCTGGACTCTCAGACTC.

2.8. Immunohistochemistry Staining Analysis

Paraffin-embedded NSCLC tissue sections were deparaffinized, hydrated, antigen retrieved, and goat serum blocked. Therefore, samples were incubated overnight at 4°C with an anti-WDR72 antibody (thermofisher PA5-63780 1 : 200) and species-specific secondary antibody for 30 min at 37°C. The immunosignal of the samples was observed utilizing a DAB solution and double-stained with hematoxylin in turn. A light microscope was used to obtain representative immunohistochemical (IHC) staining images.

2.9. Western Blot Analysis

Protein specimens from NSCLC and paracancerous tissue were isolated and transmitted to polyvinylidene fluoride (PVDF) membranes for the next steps. Following blocking with 5% BSA, PVDF membranes were incubated overnight at 4°C with primary antibodies to WDR72 (Thermo Fisher PA5-637801 : 2000) and GAPDH (1 : 1000) and then were incubated at room temperature (15-25°C) with HRP-labeled secondary antibodies for 1 h (1 : 1000) 2000). ECL (enhanced chemiluminescence) is used to capture images. GAPDH was chosen as the internal control.

2.10. Statistical Analyses

All data were normalized on gene expression using log2 transformation. Analyses of the correlation between two variables were carried out utilizing Spearman or Pearson test. was judged statistically significant. Differences between adjacent tissues and cancerous tissues were conducted employing the Wilcox test. was judged statistically significant. Kaplan-Meier curves and Cox proportional hazards regressions were employed for all survival analyses. R program (Ver. 4.1.2) was used to process the statistical significance of bioinformatic results, and GraphPad Prism (Ver. 9) was used to analyze the statistical importance of the in vitro data.

3. Results

3.1. Differential Expression of WDR72 between NSCLC Tumor and Normal Tissue Samples

We selected five datasets of GPL570 platform in GEO database, GSE19804, GSE118370, GSE19188, GSE27262, and GSE33532; “limma” R package; and “umap” R package. By drawing the PC and UMAP graphs, the samples of the two groups were separated, showing a significant variation between the two groups. The subsequent difference analysis may have more meaningful results (Figures 1(a) and 1(b)). Using the ggplot2 package of classical difference analysis-volcano plot, the threshold was selected as , value < 0.05. The figure is highlighted in blue or red, and the location of WDR72 is highlighted (Figure 1(c)). Figure 1(d) depicts the visual expression profile of the ComplexHeatmap package and shows the expression of each top 20 genes with high and low expressions. Concurrently, the Venn diagram of Figure 1(e) plots the intersecting genes of the five datasets, showing 94 intersecting genes. Then, we evaluated the expression degrees of WDR72 in 33 cancers from TCGA data. The “ggplot2” package was utilized to examine the differential expression of WDR72 in pan-cancer tissues and corresponding adjacent tissues and the expression in paired and unpaired NSCLC tissues and corresponding adjacent tissues, ns, , ; ; (Figures 2(a)2(d)). Radar visualization using the “ggradar” and “ggplot2” packages was utilized to identify WDR72 expression in pan-cancer tissues (Figure 2(e)) and adjacent pan-cancer tissues (Figure 2(f)). WDR72 was differentially expressed in 27 tumors except ACC, CESC, CHOL, DLBC, KICH, and PCPG.

3.2. Correlation Analysis of WDR72 with Clinical Factors in TCGA Database

RNAseq and clinical data in level 3 HTSeq-FPKM form in the TCGA (https://portal.gdc.cancer.gov/), LUAD (lungadenocarcinoma) and LUSC (lung squamous cell carcinoma) (lung cancer) experiments were excluded from control/normal (unapplied in projects without control/normal)+retained after clinical data availability. R package “ggplot” (Ver. 3.6.3) was used to evaluate the clinical correlations. Clinical analysis showed that the expression of WDR72 gene was related to whether smoking, age, lesion location, and TNM stage. M grade correlation is presented in Figure 3(a). PFI event correlation is presented in Figure 3(b). Smoking years correlation is presented in Figure 3(c). Tumor distribution (central or peripheral) correlation is presented in Figure 3(d). Tumor distribution (left lung or right lung) correlation is presented in Figure 3(e). T grade correlation is presented in Figure 3(f). N grade correlation is presented in Figure 3(g). Smoking correlation is presented in Figure 3(h). Age correlation is presented in Figure 3(i). Gender correlation is presented in Figure 3(j). Significance indicator: ns, ; ; ; ; Figure 4(a) is the survival analysis of WDR72 in NSCLC on the Kaplan-Meier Plotter website (https://kmplot.com/analysis). The 5-year survival rate was significantly decreased for higher WDR72 expression patients compared to those with lower WDR72 expression. Control/normal (unapplied in projects without control/normal) was removed, and the clinical data was retained. Figure 1(d) shows genes with higher expression differences in the heat map employing the R (ver. 3.6.3) “glmnet” tool (ver. 4.1-2) and survival tool (ver. 3.2-10) to draw the LASSO analysis utilizing the RMS tool (Ver. 6.2-0) and the survival package (Ver. 3.2-10) in R (ver. 3.6.3, Figure 4(b)). The COX statistical method incorporated TMN and WDR72 parameters to visualize the prognosis calibration analysis (Figure 4(c)).

3.3. PPI Network Map, GO/KEGG, and GSEA Enrichment Analysis of WDR72-Related Differential Genes

TCGA (https://portal.gdc.cancer.gov/) was utilized to obtain RNAseq information in level 3 HTSeq-Counts formed in the LUAD-LUSC (lung cancer) experiment after eliminating controls/normals (unapplied in projects without control/normal) using R (3.6.3 ver.) and “DESeq2” (ver. 1.26.0, Love MI et al. 2014) packages to target the molecule WDR72 [ENSG00000166415] in NSCLC by the reduce expression group: 0–50% and the increased expression group: 50–100% with analysis of differential genes of high correlation with WDR72. R (ver. 3.6.3) “ggplot2” tool (Ver. 3.3.3) was employed for visualization, and “clusterProfiler” package (Ver. 3.14.3) was used for WDR72 differential genes data analysis. (Figures 5(a) and 5(b)) show GO/KEGG enrichment analysis and protein interaction network, respectively. Figure 5(c) depicts the GSEA mountain map, while Figures 5(d)5(f) depict the GSEA enrichment analysis map. KEGG/GO and PPI analyses were shown in adrenergic signaling in cardiomyocytes, cytokine−cytokine receptor interaction, organelle fission, extracellular structure organization, biomineral tissue development, odontogenesis spindle, apical part of the cell, collagen−containing extracellular matrix, glycosaminoglycan binding, enzyme inhibitor activity, amide binding, receptor-ligand activity, protein digestion, absorption, and mitotic nuclear. GSEA enrichment analysis showed enrichment in reactome mitotic metaphase, anaphase, reactome RHO GTPase effectors, reactome RHO GTPases activate formins, reactome cell cycle checkpoints, reactome mitotic prometaphase, reactome mitotic spindle checkpoint, reactome separation of sister chromatids, reactome M phase, reactome signaling by RHO GTPases, reactome cell cycle mitotic, reactome resolution of sister chromatid cohesion, reactome cell cycle, reactome innate immune system, reactome rna polymerase II transcription, and reactome extracellular matrix organization.

3.4. Coexpression of WDR72 and Ferroptosis-Related Genes in NSCLC

TCGA (https://portal.gdc.cancer.gov/) was utilized to acquire RNAseq data in level 3 HTSeq-FPKM form in the LUADLUSC (lung cancer) experiment after removing control/normal (unapplied in projects without control/normal) using R (ver. 3.6.3) in the “ggplot2” package to analyze the target molecule: WDR72 (ENSG00000166415) coexpression heat map of ferroptosis-related genes in NSCLC, significance indicator: ns, ; ; ; (Figures 6(a)6(m)).

3.5. Association of WDR72 with 24 Types of Immune Cells in NSCLC

It is well known that there are typical 24 types of immune infiltrating cells, namely, aDC cells, B cells, CD8 cells, T cells, cytotoxic cells, DC cells, eosinophils cells, iDC cells, macrophages cells, mast cells, neutrophils cells, NK CD56 bright cells, NK CD56dim cells, NK cells, pDC cells, T cells, T helper cells, Tcm cells, Tem cells, TFH cells, Tgd cells, Th1 cells, Th17 cells, Th2 cells, and TReg cells. TCGA (https://portal.gdc.cancer.gov/) was employed to recover RNAseq information in level 3 HTSeq-FPKM form in LUAD-LUSC (lung cancer) experiment after removing control/normal (unapplied in projects without control/normal) using R software (GSVA package ver. 3.6.3) to examine the immune infiltration method: ssGSEA (the built-in method of the GSVA) grouped WDR72 (ENSG00000166415) by the median and calculated the higher and lower expression of WDR72 and 24 forms of immunity cells, significant threshold: ns; ; ; ; . Except for Tcm and NK CD56, 24 immune cells had a nonsignificant correlation, while the values of the other 22 cells were significant (Figures 7(a)7(m)).

3.6. Clinical Correlation of WDR72 in NSCLC

WDR72 was upregulated in NSCLC. To further validate the results of the previous analysis, we found the RNA and protein expressions of WDR72 in NSCLC. (Figures 8(a) and 8(b)) show a significantly greater level of WDR72 expression in NSCLC tissues compared to the matching normal tissues. Immunohistochemistry verified that WDR72 was strongly produced in NSCLC tissue proteins (Figures 9(a) and 9(b)). The results demonstrated that WDR72 expression level was greater in NSCLC tissues than in surrounding normal tissues. We collected clinical reports of 158 patients with NSCLC and discovered that patients with higher WDR72 levels had a poor prognosis (Figure 9(c)). Similarly, we found that WDR72 expression and smoking, tumor size, TNM staging, tumor grade, and metastasis were correlated (, Table 2). Multivariate analysis revealed that tumor size, TNM staging, tumor grade, and metastasis were associated with the expression of WDR72 (, Table 3).

4. Discussion

Lung cancer has two major histological forms: small-cell and non-small-cell. Genomic alterations have been found in LUSC and SCC. Lung cancer causes the greatest cancer-related deaths among males and females and inflammation and environmental risk factors such as smoking significantly contribute to lung cancer progression [27]. We analyzed the clinical significance of WDR72 in NSCLC data from the TCGA database WDR72. We discovered a close correlation with smoking, and the duration of smoking in years is also positively related to WDR72 expression. Concurrently, we used our clinical data to conduct univariate and multivariate analysis and found consistent results with the bioinformatic analysis; accordingly, WDR72 can be used as an NSCLC-target gene.

Ferroptosis was first identified as erastin-induced death and found to be involved in various pathophysiological processes and diseases. Since certain oncogenically mutated tumor cells have increased susceptibility to ferroptosis, ferroptosis can be triggered, and ferroptosis-sensitive tumor cells may also have significant therapeutic potential. Therefore, modulation of ferroptosis may have clinical potential in disorders related to ferroptosis [28]. The high-iron concentration of cancer cells and their increased susceptibility to developing ferroptosis are promising for cancer therapy [29]. Herein, we analyzed the correlation between WDR72 and ferroptosis-related genes using the R language and found that WDR72 was related to most ferroptosis-related genes. It is speculated that WDR72 may affect non-small-cells through the ferroptosis pathway of lung cancer progression.

TME is critical in cancer development, and different microenvironment signals have tumor-promoting and tumor-suppressing effects [30]. In TME, injury-related signals can influence the phenotype and activation condition of tumor cells and infiltrate immune cells [31]. The connection between the inflammatory microenvironment and tumors is bidirectional and dynamic and has both tumor-promoting and tumor-suppressing properties that it is not only the basis for the onset and progression of tumors [32]. Correlation analysis was performed in the NSCLC data, except for NK CD56dim cells and Tcm. In contrast, the other 22 cells were related, so we confirmed that WDR72 changes the TME through immune infiltration and impacts the progression of NSCLC.

Finally, our enrichment analysis indicated that WDR72 might be enriched through the cell cycle, immune system, RNA polymerase II transcription, and extracellular matrix organization. This suggested that it may affect tumor growth and metastasis through these pathways. Moreover, we conducted a series of tests to examine WDR72 expression in NSCLC tissues. WDR72 was elevated in NSCLC tissues relative to surrounding normal tissue.

Our study has some limitations. First, in vitro cell experiments are required to study WDR72 expression and function further. Second, our data suggest that WDR72 can act as a prognostic factor in NSCLC, which needs additional validation. Third, the presence of WDR72 on the TME and immunotherapy requires in vitro and in vivo experimental and clinical validation. Fourth, although we confirmed WDR72 expression in human NSCLC tissues, its exact regulatory mechanism remains unclear. WDR72 can participate in some biological processes in cells, including cell death, cell division, protein synthesis, and membrane transport. Ferroptosis is also a kind of programmed death. According to our research and analysis, WDR72 has a high correlation with many molecules of the ferroptosis pathway. Is it possible that WDR72 affects NSCLC through the ferroptosis pathway? In the next study, we will focus on the mechanism research.

5. Conclusion

Herein, the analysis of WDR72 suggested that WDR72 may be a prognostic factor in NSCLC, which was expressed at the mRNA and protein levels among NSCLC tumors and normal tissues. We initially revealed the association between WDR72 and ferroptosis and immune infiltration. Furthermore, WDR72 expression was linked to MSI, TMB, and immune cell infiltration in NSCLC. These findings may elucidate the aim of WDR72 in the incidence and growth of NSCLC and provide a reference for future patients with NSCLC to receive more accurate and personalized immunotherapy.

Data Availability

The utilized and examined datasets (Table 1) are available upon reasonable request from the corresponding author.

Ethical Approval

The ethics committee of The Sixth People’s Hospital of Nantong and Affiliated Hospital of Nantong University, Jiangsu, China, approved the research. All human research was done following the most recent iteration of the Declaration of Helsinki.

Participants provided an informed permission before this research.

Conflicts of Interest

These authors declare no conflict of interest.

Authors’ Contributions

Guanglin Shi did the experiments and wrote the document. Hua Huang and Zhongbo Yin greatly contributed to pathological specimens and immunohistochemistry examination. Qinggan Ni helped significantly with the analysis and document preparation. Yuqing Miao helped with the analysis with productive discussions and contributed reagents and materials. Weirong Shi did the data analyses. Minhua Shi designed the experiments. All researchers read and approved the final document. Guanglin Shi and Qinggan Ni contributed equally to this work.

Acknowledgments

This work was supported by the Nantong Municipal Commission of Health and Family Planning (grant number MA2021030) and the Science and Technology Project of Nantong City (grant number jcz18109).