Abstract
Background:
Elevated PPP4C expression has been associated with poor prognostic implications for patients suffering from lung adenocarcinoma (LUAD). The extent to which PPP4C affects immune cell infiltration in LUAD, as well as the importance of associated genes in clinical scenarios, still requires thorough investigation.
Methods:
In our investigation, we leveraged both single-cell and comprehensive RNA sequencing data, sourced from LUAD patients, in our analysis. This study also integrated datasets of immune-related genes from InnateDB into the framework. Our expansive evaluation employed various analytical techniques; these included pinpointing differentially expressed genes, constructing WGCNA, implementing Cox proportional hazards models. We utilized these methods to investigate the gene expression profiles of PPP4C within the context of LUAD and to clarify its potential prognostic value for patients. Subsequent steps involved validating the observed enhancement of PPP4C expression in LUAD samples through a series of experimental approaches. The array comprised immunohistochemistry staining, Western blotting, quantitative PCR, and a collection of cell-based assays aimed at evaluating the influence of PPP4C on the proliferative and migratory activities of LUAD cells.
Results:
In lung cancer, elevated expression levels of PPP4C were observed, correlating with poorer patient prognoses. Validation of increased PPP4C levels in LUAD specimens was achieved using immunohistochemical techniques. Experimental investigations have substantiated the role of PPP4C in facilitating cellular proliferation and migration in LUAD contexts. Furthermore, an association was identified between the expression of PPP4C and the infiltration of immune cells in these tumors. A prognostic framework, incorporating PPP4C and immune-related genes, was developed and recognized as an autonomous predictor of survival in individuals afflicted with LUAD. This prognostic tool has demonstrated considerable efficacy in forecasting patient survival and their response to immunotherapeutic interventions.
Conclusion:
The involvement of PPP4C in LUAD is deeply intertwined with the tumor’s immune microenvironment. PPP4C’s over-expression is associated with negative clinical outcomes, promoting both tumor proliferation and spread. A prognostic framework based on PPP4C levels may effectively predict patient prognoses in LUAD, as well as the efficacy of immunotherapy strategy. This research sheds light on the mechanisms of immune interaction in LUAD and proposes a new strategy for treatment.
1 Introduction
Lung adenocarcinoma (LUAD), a prevalent form of lung cancer, remains a leading cause of cancer-associated mortality globally (1), displaying a five-year survival rate of merely 15% (2–4). Among the primary therapeutic targets for LUAD are PD-1 and PD-L1 (5–8). The field has seen considerable interest in immune checkpoint blockade (ICB), noted for substantially enhancing survival rates in cancer patients. Despite these advances, the effectiveness of immunotherapy is not universal among LUAD patients (9). Factors such as the PD-1 expression, MSI, TMB, and tumor microenvironment play crucial roles in determining the response to immune treatments (10–13). Nevertheless, the prognostic tools currently available lack precision, underscoring the imperative need for research into vital genes and biomarkers that influence both the prognosis of LUAD patients and the response to immunological therapies (14–17).
Nowadays, high-throughput sequencing tools such as second-generation sequencing (18, 19), single-cell transcriptomics sequencing (20), spatial transcriptomics sequencing, proteomics sequencing, and metabolomics sequencing are rapidly evolving (21). These technological advances have brought new perspectives to the field of genetic research. Utilizing patterns of gene expression, researchers aim to uncover novel biomarkers within various cancers (22, 23). By rigorously analyzing high-throughput RNA sequencing data, it is possible to support initiatives in personalized treatment and precision medicine. This includes identifying emerging prognostic markers and therapeutic targets, elucidating principal genes that affect the immune infiltration status of patients, and delineating the molecular pathways that promote the progression of LUAD (24, 25). These objectives can be fulfilled through the application of differential expression analysis and functional enrichment studies (26).
Protein phosphatases from the PPP family play crucial roles in a variety of physiological and pathological contexts, notably in oncological disorders (27). In particular, the enzyme PPP4C has been linked to the advancement of numerous malignancies, such as breast and pancreatic cancers, and glioblastomas (28). The connection between the expression levels of PPP4C, immune cell infiltration, and the efficacy of chemoimmunotherapy treatments has been established (29–32). However, the specific role of PPP4C within the framework of LUAD has yet to be clarified.
To address the challenge of identifying potential therapeutic targets and prognostic biomarkers for lung adenocarcinoma, our study focused on analyzing the differential immunogenomic profiles in both tumor and non-tumor tissues from LUAD patients. Employing methodologies such as gene set enrichment analysis, survival analysis, WGCNA, single-cell transcriptomics analysis and assessment of differentially expressed immune-related genes, we pinpointed PPP4C as a pivotal subject for further investigation (33). Our research further explored how PPP4C expression within the tumor immune microenvironment influences LUAD prognosis. By constructing a risk assessment model incorporating PPP4C alongside other immunologically relevant genes, we successfully predicted outcomes and immune response efficacy in LUAD therapy. Furthermore, our investigation extended to analyzing clinical and immune profiles across various risk categories, culminating in the development of diagnostic charts that offer novel insights into managing LUAD and evaluating immunotherapy responses in treated patients.
2 Materials and methods
2.1 Participant recruitment and data collection
Between June and August 2023, individuals undergoing surgical procedures at Harbin Medical University’s Second Affiliated Hospital contributed samples from three distinct pairs of lung adenocarcinoma and adjacent non-malignant tissues. These specimens encompassed paired samples of lung adenocarcinoma tumors and the corresponding adjacent healthy tissues. Authorization for this research was granted by the Ethics Committee of the Second Affiliated Hospital at Harbin Medical University under the protocol number KY2023–042.
Clinical records and RNA sequencing data from 594 LUAD patients were retrieved from the TCGA repository (https://portal.gdc.cancer.gov/). Furthermore, RNA sequencing and survival statistics for an additional 246 LUAD patients were accessed through GEO dataset GSE31210 (34). The dataset GSE123902, which includes RNA sequencing data at the single-cell level, comprised thirteen tumor and four normal tissue samples. An immunotherapy validation cohort was established using data from the R program IMvigor210. Additionally, gene sets pertinent to immune responses were obtained from both InnateDB and ImmPort databases.
2.2 Identification of immune-related genes
To ascertain immune-associated genes, differential expression analysis was conducted on LUAD cohort data from the TCGA database, employing the limma software (|log2FC| > 0.585) (35). Subsequently, lists of immune genes sourced from the InnateDB and ImmPort databases were amalgamated to facilitate functional enrichment analysis. Genes within the yellow and brown modules were selected based on a significance threshold of P < 0.05. The identification of pivotal genes in immune regulation was achieved through the application of the WGCNA algorithm (33, 36). Functional enrichment analysis was subsequently conducted using the clusterProfiler software (37).
2.3 Single-cell RNA sequencing
In this investigation, data from single-cell RNA sequencing was analyzed using the Seurat package within R software (version 4.4.0) (38). Initially, rigorous quality control protocols were employed to eliminate substandard cells, which required the establishment of specific criteria such as nFeature between 300 and 7,500, nCount between 300 and 100,000, mitochondrial gene expression ratio of less than 20%, ribosomal gene expression ratio of greater than 3 and erythrocyte gene expression ratio of less than 0.1. Subsequently, the SCTransform method was applied for data normalization, followed by dimensionality reduction via the RunTSNE function, facilitating easier clustering and visualization. Knowledge from prior studies and the CellMarker database (http://xteam.xbio.top/CellMarker) aided in cell type annotation (39). Furthermore, normal epithelial cells in the sample were used as a control group, and the inferCNV (https://github.com/broadinstitute/inferCNV/) algorithm was utilized to differentiate between benign and malignant epithelial cells based on differences in chromosome copy number variants and to explore the differences in the expression of PPP4C in the two types of cells, Then, tumor cells were categorized according to the median PPP4C expression level, and differential pathway enrichment was explored using GSVA analysis, and CellPhoneDB (version 2.0) elucidated cell-cell interactions, with higher ligand-receptor interaction scores suggesting more robust intercellular communication.
2.4 Immune infiltration
LUAD patients were stratified into groups with high and low levels of PPP4C gene expression based on median values. Subsequent analysis of immune cell infiltration differences between these groups employed the CIBERSORT algorithm (40). Comparative immunological assessments utilized the ESTIMATE method (41, 42), alongside Spearman’s rank correlation analysis. To delineate key genetic markers for constructing a prognostic model, the LASSO and univariate regression methodologies isolated 14 critical genes (43, 44). Furthermore, multivariable Cox regression analysis facilitated the establishment of risk scores, categorizing patients into cohorts with elevated or diminished risk profiles.
2.5 The immunological landscapes, clinical manifestations, and immunotherapy
To elucidate the interplay between immunological profiles, clinical characteristics, and responses to immunotherapy, our study employed MCPcounter to quantify immune cell infiltration within distinct risk strata. This approach facilitated an exploration of the correlations between risk assessments and clinical manifestations. Furthermore, the ESTIMATE algorithm was applied to evaluate differences in immunological scores across the cohorts (45). We investigated disparities in the expression of genes associated with immune checkpoints across these risk groups. Subsequent analysis focused on contrasting the tumor mutational burden between individuals classified as high-risk and those deemed low-risk, alongside investigating variations in prognostic outcomes, responses to immune treatments, and staging utilizing the R IMvigor210 toolset. In pursuit of developing a predictive framework, we conducted both univariate and multivariate Cox regression analyses. The decision curve analysis (DCA) was subsequently utilized to determine the predictive precision and clinical utility of the developed model.
2.6 Real-time PCR
In this investigation, the methodology of real-time quantitative polymerase chain reaction (RT-qPCR) was utilized. For the synthesis of cDNA, cellular RNA was initially isolated using Trizol (Sigma), followed by reverse transcription using PrimeScriptTM RT Master Mix (TaKaRa). Subsequently, the ChamQ Universal SYBR qPCR Master Mix (Vazyme) facilitated the RT-qPCR analysis. The quantification of mRNA expression levels of the target genes was performed utilizing the 2-ΔΔCT approach, wherein β-actin served as the internal standard (46). Data analysis and the creation of graphical outputs were conducted with the aid of Prism software.
The primer sequences used were as follows: for PPP4C, 5’- GGTCTATGGCTTCTACGATG -3’; and for β-actin, 5’-GAAGAGCTACGAGCTGCCTGA-3’.
2.7 Western blot
Proteins were extracted from washed tissues using a RIPA buffer enhanced with PMSF and a cocktail of protease inhibitors. Subsequently, protein concentrations were quantified employing the BCA assay. Proteins were resolved by SDS-PAGE and transferred onto PVDF membranes. The membranes underwent a blocking process using 5% non-fat milk before antibody incubation. Primary antibodies were applied overnight at 4°C, followed by a one-hour incubation with secondary antibodies at ambient temperature (47). In this study, the primary antibodies utilized were anti-β-actin and anti-PPP4C, while Goat anti-Rabbit HRB served as the secondary antibody. Visualization of the protein bands on the Western blots was facilitated using an ECL detection system.
2.8 Immunohistochemistry
In the course of this study, the immunohistochemistry protocol was meticulously followed through several distinct methodological phases (48). Initially, tissue sections were subjected to a pre-treatment phase within a temperature-controlled oven. Ensuring cellular exposure to the antigen involved multiple steps including the removal of paraffin and the retrieval of antigens. Overnight incubation with primary antibodies, sourced from Abcam, was then conducted. Visualization of the staining process was facilitated using the Dako EnVisionTM FLEX+ kit. Subsequently, an Aperio digital pathology slide scanner was employed to capture the resultant images, and additional staining was performed using hematoxylin provided by Sigma Aldrich. The semi-quantitative analysis of immunohistochemical staining was performed using the IHC Profiler plugin in Image J software, and the IHC score was calculated based on the percentage of area of different staining intensities. (IHC sore = High positive areas percentage × 3 + Positive areas percentage × 2 + Low positive areas percentage × 1 + Negative areas percentage × 0)
2.9 Cultivation of lung cancer cells
Cultures of the A549, H1299, PC9 and H1975 cell lines (sourced from PLST Co., Ltd., China) were propagated using RPMI-1640 medium (manufactured by Gibco, USA), which was enriched with 10% fetal bovine serum. BEAS-2B cell line culture was performed in serum-free BEpiCM complete medium (manufactured by SclenCell, USA). These cultures were sustained at a constant temperature of 37°C within a humidified environment, supplemented by an atmosphere containing 5% CO2.
2.10 PPP4C influences on the migration and invasion of LUAD cells
In the investigation of PPP4C’s role in inhibiting LUAD cell migration, cells were evenly distributed across six-well plates for a wound healing assay, achieving the necessary confluency. Subsequently, serum-free DMEM was employed to foster cell cultures post-creation of uniform scratches using 200 μl pipette tips. After 24 hours of cell culture, the scratched areas were systematically photographed to monitor the healing progress, thus helping to analyze the migration dynamics of the cells.
To determine the invasive properties of the tumor cells, experiments utilized a 24-well, 8μm Transwell setup (NEST Biotechnology Co. LTD., Wuxi, China). Single-cell suspensions were introduced into the upper compartment at a concentration of 1×10^5 cells per well, with serum-free media above and 10% FBS-enriched media below. Following a 24-hour incubation, the cells underwent fixation with 4% paraformaldehyde post-staining with 1% crystal violet (49). Invasion assessments were conducted using an inverted microscope to capture images from at least three random fields.
2.11 Cell colonies formation
To investigate colony formation, 500 cells were uniformly distributed into each well of a six-well plate. These cells were incubated at 37°C for two weeks until distinct colonies, comprising no fewer than 50 cells each, could be discerned microscopically. Post-incubation, colonies were fixed using 4% paraformaldehyde and subsequently stained with 1% crystal violet. Photographs documenting colony growth across all wells were subsequently captured (50).
2.12 Statistical analysis
To conduct the statistical evaluations, we utilized versions 4.2.2 and 9.5.1 of Prism and R, respectively. Continuous variables were described using means and standard deviations of independent samples. Survival differences were assessed employing Kaplan-Meier curves and the Log-rank test, alongside Spearman’s correlation coefficient and the t-test for comparisons between two groups. A significance level was established at P < 0.05.
3 Results
3.1 Analysis of differentially expressed immune genes in LUAD
Utilizing the TCGA-LUAD data set for analysis, a survey of gene expression variances unveiled 13,618 genes with fluctuating expression levels. Among these, 4,237 were noted as down-regulated and 9,381 as up-regulated (Figure 1A). Subsequent scrutiny against immune-specific databases from InnateDB and ImmPort revealed 1,011 genes linked to immune functionalities showing disparate expression patterns in LUAD, divided into 499 down-regulated and 512 up-regulated genes (Figure 1B). Advanced probing into these genes’ roles identified enrichment in 132 KEGG pathways and 2,720 GO terms, with emphasis on the foremost 30 GO terms and KEGG pathways (Figures 1C, D).
Figure 1
3.2 Characterization of immune-related gene targets in LUAD
In an endeavor to identify pivotal hub genes associated with immunity in lung adenocarcinoma (LUAD), researchers conducted a comprehensive analysis using Weighted Gene Co-expression Network Analysis (WGCNA) on 1011 genes known for their variable expression in immune responses. This analysis constructed a scale-free network, which facilitated the identification of an optimal soft-thresholding power set at 2, as depicted in Figure 2A. This specific thresholding power, when applied alongside hierarchical clustering employing average linkage techniques, classified these genes into seven distinct clusters (Figures 2B, C). Subsequent analysis calculated Pearson correlation coefficients to link these gene clusters with particular traits of LUAD tumor samples, adopting a significance level of P < 0.05. Notably, genes within the brown and yellow clusters were selected for further detailed scrutiny. The intent of this deeper investigation was to elucidate the relationship between immune genes and patient prognosis in LUAD. Prior to a focused analysis of the selected gene clusters, a preliminary survey was conducted to assess their functional involvement and pathway engagement. This led to the revelation of considerable enrichment in the top 8 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, shedding light on their biological functions and interconnected roles in LUAD progression (Figures 2D, E).
Figure 2
3.3 Probing PPP4C expression patterns in lung cancer
To ascertain immune-associated genes of independent prognostic value, a meticulous univariate Cox regression coupled with Kaplan-Meier survival analysis was conducted for genes within the brown and yellow modules (Figures 3A-G). This analysis brought to light 58 immune genes of prognostic significance. Moreover, earlier investigations have illuminated distinctive expression profiles of the PPP gene family in breast cancer, underscoring its link to both prognosis and the infiltration states of immune cells (21). Notably, PPP4C has been distinguished as a valuable prognostic indicator and as a potential therapeutic target in breast cancer. However, the expression dynamics, biological functions, and impact on immune cell infiltration by PPP4C in lung adenocarcinoma (LUAD) necessitate further inquiry. Thus, our focus has shifted toward the PPP4C gene. Through the analysis of the GSE31210 dataset from the GEO database, we confirmed PPP4C’s unique expression in LUAD and its prognostic ramifications. The findings reveal a significant upregulation of PPP4C in LUAD, correlating with unfavorable outcomes (Figures 3H, I). Subsequently, an exhaustive pan-cancer assessment of PPP4C was undertaken to define its biological ramifications across various cancers (Figure 3J). This broad evaluation disclosed a uniform pattern of elevated PPP4C expression in numerous cancers, notably associated with adverse results, particularly in patients with liver hepatocellular carcinoma (LIHC) and kidney renal clear cell carcinoma (KIRC) (Figure 3K). In conclusion, the pervasive increase in PPP4C expression in LUAD and other cancers highlights its critical role as an indicator of poor prognosis, thereby solidifying PPP4C as a paramount prognostic marker and a prospective therapeutic target.
Figure 3
3.4 Delineation of PPP4C expression via single-cell profiling
A single-cell transcriptomic approach to investigating PPP4C expression has uncovered a considerable variation in expression levels and functional roles across cells. By examining the GSE123902 dataset, which encompasses primary, metastatic lung adenocarcinomas, and non-tumorous lung tissues, the study distinguished twenty distinct cellular clusters. These clusters were systematically categorized as epithelial, immune, or stromal cells through integration with the CellMarker database (Figures 4A-C). Deeper scrutiny revealed seven principal subtypes within the immune cells (Figures 4D-F). The classification strategy hinged on the use of defined cellular markers (Figures 4G, H). Subsequent analysis using inferCNV highlighted a conspicuous escalation in PPP4C expression among cancerous epithelial cells (Figure 5A), compared to their normal counterparts (Figure 5B). Further subdivision of the tumor epithelial cells into high and low PPP4C expression groups, followed by GSVA enrichment analysis, identified significant upregulation of pathways related to the cell cycle, DNA replication, oxidative phosphorylation, and metabolism of glycolysis/TCA cycle in cells with elevated PPP4C expression (Figure 5C). Moreover, a comparative analysis of intercellular communication between different factions of cancerous epithelial and immune cells was conducted. This investigation not only casts light on the crucial role of PPP4C in lung adenocarcinoma pathogenesis but also lays the groundwork for further explorations (Figures 5D, E).
Figure 4
Figure 5
3.5 Exploration of immunological subtypes based on PPP4C gene expression
Current research explores the dynamic interaction between PPP4C gene expression and immune cell infiltration, highlighting its importance in modern scientific investigations. This study stratifies patients into categories based on either elevated or reduced PPP4C gene expression levels. Utilizing CIBERSORT software to analyze their sequencing data, the study assesses how PPP4C gene expression impacts the distribution of twenty-two distinct immune cell types (Figure 6A). Results reveal that patients with increased PPP4C expression show higher levels of follicular helper T cells, regulatory T cells (Tregs), and both M0 and M1 macrophages. In contrast, individuals with lower expression levels display a higher prevalence of memory B cells, CD4 memory resting T cells, monocytes, and inactive mast cells (Figure 6B). Additionally, employing the ESTIMATE algorithm to determine stromal and immune scores within these gene expression subclasses demonstrates a negative correlation between PPP4C gene expression and the immune score (Figures 6C, D). Such findings provide crucial insights into the role of the PPP4C gene in immune regulation.
Figure 6
3.6 Constructing and validating a predictive model involving PPP4C
Given the correlation between elevated expression of the PPP4C gene and reduced survival rates in patients with lung adenocarcinoma (LUAD), as well as its ties to immune cell infiltration in the tumor environment, we developed a prognostic model based on PPP4C. This model aims to predict survival and immunotherapy outcomes for LUAD patients. We initially conducted an analysis to identify genes correlated with PPP4C expression and immune parameters, revealing a set of 267 genes (Figure 7A). Further enrichment analysis linked these genes to GTPases, pathways of non-small cell lung cancer, PD-L1 expression, and the PD-1 checkpoint, underscoring their potential impact on LUAD prognosis and the efficacy of PD-L1-targeted therapies (Supplementary Figure S1). A selected group of 100 genes with significant prognostic value was further narrowed down to 14 key genes through LASSO regression analysis (Figures 7B, C). Utilizing these key genes, a risk stratification model was constructed, and risk scores for individual samples were determined via multivariate COX regression analysis. The effectiveness of this model was confirmed through PCA analysis (Figure 7D), and the creation of Kaplan-Meier survival and ROC curves demonstrated significantly lower survival in the high-risk group compared to the low-risk group (Figures 7E, F). Further validation of the model’s robustness and precision was conducted using a GEO dataset, reinforcing the prognostic model’s accuracy and dependability in predicting survival and responses to immunotherapy in LUAD patients (Figures 7G, H).
Figure 7
3.7 Assessing the clinical impact of a PPP4C-centric prognostic framework
The prognostic model based on PPP4C demonstrates significant potential for predicting survival rates in patients with lung adenocarcinoma. This approach involves stratifying patients into distinct groups based on clinical features including TNM classification and age. The analysis aims to explore the relationship between these clinical variables and prognostic scores. Results indicate a positive correlation between higher prognostic scores and advanced T, N (Figures 8A, B), and pTNM stages (Figures 8C, D), particularly in patients under the age of 65 (Figure 8E). An extensive analysis was conducted to assess differences in immune cell composition among cohorts displaying varied prognostic scores. Notably, the high-risk group exhibited a reduced presence of T cells CD4 memory activated, Tregs, activated NK cells, and M0 macrophages compared to the low-risk group, which presented an inverse pattern (Figure 9A). The variation in immune cell and score metrics among these groups was evaluated using the ESTIMATE algorithm and MCP counting method (Figures 9B, C). Furthermore, an analysis comparing immune checkpoint gene expression across different prognostic groups was undertaken (Figure 9D). The model’s accuracy in predicting outcomes after immunotherapy was validated (Figure 9E), with the high-risk group showing a decreased likelihood of favorable response (Figure 9F). The study also included an examination of the disparity in tumor mutation burden between high and low-risk groups (Figures 9G, H). These findings underscore the critical clinical importance of utilizing a PPP4C-focused prognostic model for lung adenocarcinoma patients, providing key insights for optimizing therapeutic strategies.
Figure 8
Figure 9
3.8 Formulation of a nomogram for precise prognosis prediction in LUAD patients
To enhance the precision of predictive frameworks in evaluating patient risks, our research integrated clinical-pathological factors with risk scores utilizing both single-factor and multifactorial Cox regression analyses, considering elements like patient demographics. This examination demonstrated that in lung adenocarcinoma (LUAD) scenarios, the derived risk score is an essential independent prognostic factor (P < 0.001), whereas the clinical-pathological stage was of secondary importance (Figures 10A, B). Subsequently, a Nomogram was developed, merging both staging data and risk assessments (Figures 10C). The prognostic accuracy of this Nomogram was subsequently confirmed through methods of calibration and decision curve analysis (DCA). Calibration diagrams confirmed the Nomogram’s ability to accurately reflect the actual survival rates of patients across timeframes of 1, 3, and 5 years (Figures 10C). Moreover, the DCA underscored that the Nomogram markedly surpasses either the risk scores or the pathological staging alone in identifying high-risk patients, thus delivering an enhanced net clinical benefit.
Figure 10
3.9 Investigating the impact of PPP4C gene expression on cellular dynamics in LUAD
In investigating lung adenocarcinoma, a comprehensive study was conducted to assess the differential expression of PPP4C and its effects on cellular dynamics. The initial stage involved quantifying PPP4C expression in lung adenocarcinoma tissue as well as in adjacent non-cancerous tissue, utilizing methods such as Western blot analysis (Figures 11A, B), RT-qPCR (Figure 11C), and immunohistochemistry (Figures 11D, E). The data revealed a pronounced elevation of PPP4C levels in the tumor tissues compared with the normal tissues, confirmed at both protein and mRNA levels. Meanwhile, we also compared the expression of PPP4C in normal lung epithelial cell lines and lung adenocarcinoma cell lines by Western blot experiments. The results showed that PPP4C was significantly overexpressed in all four lung adenocarcinoma cell lines, which also laid the foundation for our subsequent cell function experiments (Figures 11F, G). Subsequent experiments across various lung adenocarcinoma cell lines were designed to clarify the role of PPP4C. By creating cellular models that either overexpressed (Figures 12A, C) or suppressed PPP4C expression (Figures 12B, D), it was observed that increased PPP4C levels markedly enhanced the proliferation (Figures 13A, B), and also the migration and invasion capabilities of A549 and H1299 cells (Figures 13C, E). Conversely, a reduction in PPP4C expression led to decreased these cellular processes (Figures 13D, F). These results underscore the pivotal influence of PPP4C in lung adenocarcinoma progression and suggest its potential association with patient prognosis.
Figure 11
Figure 12
Figure 13
4 Discussion
PPP4C, a member of the phosphatase enzyme family, demonstrates widespread expression across diverse human tissues, suggesting its implication in crucial biological processes. Its conservation throughout evolution underscores its putative involvement in fundamental physiological mechanisms (51). Prior investigations have noted an upregulation of PPP4C expression in numerous cancer types, including colorectal, breast, and pulmonary malignancies (52–55). Within the oncogenic landscape, heightened PPP4C expression correlates with the modulation of pivotal signaling cascades such as mTOR, JNK, and NF-κB, while its depletion may instigate cellular apoptosis (56–58). Notably, PPP4C-mediated augmentation of the ERK pathway fosters lung cancer cell proliferation and impedes apoptotic mechanisms, thereby exacerbating clinical prognosis (59). These insights underscore the potential of PPP4C as a promising therapeutic target across diverse cancer types. Nonetheless, further investigation is imperative to delineate the impact of PPP4C on the immune milieu in lung adenocarcinoma and elucidate its interplay with other oncogenic determinants.
This investigation employed diverse methodologies, encompassing RNA sequencing data analysis, weighted gene co-expression network analysis (WGCNA), and single-cell transcriptome sequencing, to identify pivotal genes within patients afflicted with LUAD and assess their impact on the tumor’s immune milieu. The study revealed a conspicuous escalation in the levels of PPP4C gene among individuals with LUAD, a phenomenon associated with an unfavorable prognosis, a conclusion substantiated by experimental verification. Moreover, an exploration into the nexus between PPP4C gene expression and the tumor’s immune milieu was conducted. The results elucidate a direct correlation between heightened PPP4C expression and diminished immune scores, underscored by a negative correlation between PPP4C expression and immune score. Furthermore, an augmented presence of regulatory T cells (Tregs) was observed in specimens exhibiting elevated PPP4C expression. The dampening of immune response by Tregs emerges as a pivotal tactic employed by tumors to evade immune surveillance, a factor potentially underpinning the adverse prognosis for LUAD patients with heightened PPP4C expression (60). These revelations furnish crucial insights into comprehending the role of PPP4C in LUAD and its ramifications for therapeutic interventions.
Given the significance of PPP4C in prognosis and its impact on the tumor immune microenvironment, we have devised a risk evaluation framework integrating PPP4C with genes associated with immune scoring. This framework serves as a robust prognostic tool, providing insights into patients’ probable response to immunotherapy. Integrating this risk assessment with conventional TNM classifications substantially enhances the predictive capacity of the model. Our investigation encompasses 14 genes pivotal in lung adenocarcinoma (LUAD). Notably, DOCK4 emerges as a suppressor of tumor growth by modulating tumor cell adhesion and invasiveness (61). Moreover, EFHD2 garners attention for its indispensable role in activating immune cells and promoting cancer dissemination (62). METTL7A and MT2A, among other genes, correlate with the prognostic landscape of LUAD patients, implicating their potential impact on disease outcomes (63–65). These revelations underscore a nuanced interplay between these genes and PPP4C’s involvement in LUAD, underscoring the necessity for further exploration to elucidate their functions.
In summary, this study elucidates key genetic markers that impact prognostic outcomes and the immunological microenvironment in patients diagnosed with Lung Adenocarcinoma (LUAD) by integrating single-cell analysis and multi-omics strategies. And the regulatory effect of PPP4C on lung adenocarcinoma cells was verified by differential expression assay and cell function assay. These findings pave the way for identifying new therapeutic targets and prognostic indicators for managing LUAD. Despite the reliance on publicly available datasets, these findings provide valuable insights for further investigations into the oncological significance of PPP4C. Future studies should consider integrating advanced techniques such as machine learning algorithms, spatial transcriptome analysis, and further experimental validation to unravel the biological mechanisms of PPP4C in cancer. These efforts aim to strengthen the basis of precision oncology.
5 Conclusion
In patients diagnosed with LUAD, a robust correlation has been established between the levels of PPP4C and the intricacies of the tumor microenvironment’s immune landscape. Elevated PPP4C concentrations serve as harbingers of unfavorable prognostic outcomes, thereby fostering the proliferation and metastasis of lung carcinoma cells. The development and implementation of a risk assessment paradigm centered on PPP4C afford precise prognostic capabilities for individuals afflicted with LUAD, concurrently facilitating the meticulous evaluation of the efficacy of immunotherapeutic interventions. This scholarly inquiry has not only elucidated the intricate pathways governing immune responses in LUAD but also delineated a strategic blueprint for the efficacious management of cancer.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving humans were approved by The Ethics Committee of the Second Affiliated Hospital of Harbin Medical University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
KW: Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft. BP: Data curation, Formal analysis, Software, Visualization, Writing – review & editing. RX: Data curation, Formal analysis, Software, Visualization, Writing – review & editing. TL: Data curation, Formal analysis, Software, Visualization, Writing – review & editing. XC: Data curation, Software, Visualization, Writing – review & editing. ZS: Data curation, Software, Visualization, Writing – review & editing. JS: Data curation, Software, Visualization, Writing – review & editing. ML: Investigation, Resources, Writing – review & editing. CW: Investigation, Resources, Writing – review & editing. XZ: Investigation, Resources, Writing – review & editing. CX: Investigation, Resources, Writing – review & editing. HC: Conceptualization, Investigation, Supervision, Validation, Writing – review & editing. LZ: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
All authors acknowledge the contributions from the TCGA portal.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2024.1416632/full#supplementary-material
Supplementary Figure 1Enrichment analysis of PPP4C and immune score intersection genes. (A-C) GO enrichment analysis of intersecting genes (P <0.05). (D) KEGG enrichment analysis of intersecting genes (P <0.05).
References
1
SiegelRLMillerKDJemalA. Cancer statistics, 2019. CA Cancer J Clin. (2019) 69:7–34. doi: 10.3322/caac.21551
2
SchmidtLEskiocakBKohnRDangCJoshiNSDuPageMet al. Enhanced adaptive immune responses in lung adenocarcinoma through natural killer cell stimulation. Proc Natl Acad Sci. (2019) 116:17460–9. doi: 10.1073/pnas.1904253116
3
LavinYKobayashiSLeaderAAmirE-aElefantNBigenwaldCet al. Innate immune landscape in early lung adenocarcinoma by paired single-cell analyses. Cell. (2017) 169:750–765.e717. doi: 10.1016/j.cell.2017.04.014
4
ChenHCarrot-ZhangJZhaoYHuHFreemanSSYuSet al. Genomic and immune profiling of pre-invasive lung adenocarcinoma. Nat Commun. (2019) 10:5472. doi: 10.1038/s41467-019-13460-3
5
MotzerRJEscudierBMcDermottDFGeorgeSHammersHJSrinivasSet al. Nivolumab versus everolimus in advanced renal-cell carcinoma. N Engl J Med. (2015) 373:1803–13. doi: 10.1056/NEJMoa1510665
6
LarkinJChiarion-SileniVGonzalezRGrobJJCoweyCLLaoCDet al. Combined nivolumab and ipilimumab or monotherapy in untreated melanoma. N Engl J Med. (2015) 373:23–34. doi: 10.1056/NEJMoa1504030
7
BellmuntJde WitRVaughnDJFradetYLeeJLFongLet al. Pembrolizumab as second-line therapy for advanced urothelial carcinoma. N Engl J Med. (2017) 376:1015–26. doi: 10.1056/NEJMoa1613683
8
KumarSChatterjeeMGhoshPGangulyKKBasuMGhoshMK. Targeting PD-1/PD-L1 in cancer immunotherapy: An effective strategy for treatment of triple-negative breast cancer (TNBC) patients. Genes Dis. (2023) 10:1318–50. doi: 10.1016/j.gendis.2022.07.024
9
BrahmerJReckampKLBaasPCrinòLEberhardtWEPoddubskayaEet al. Nivolumab versus docetaxel in advanced squamous-cell non-small-cell lung cancer. N Engl J Med. (2015) 373:123–35. doi: 10.1056/NEJMoa1504627
10
NishinoMRamaiyaNHHatabuHHodiFS. Monitoring immune-checkpoint blockade: response evaluation and biomarker development. Nat Rev Clin Oncol. (2017) 14:655–68. doi: 10.1038/nrclinonc.2017.88
11
SamsteinRMLeeCHShoushtariANHellmannMDShenRJanjigianYYet al. Tumor mutational load predicts survival after immunotherapy across multiple cancer types. Nat Genet. (2019) 51:202–6. doi: 10.1038/s41588-018-0312-8
12
PatelSPKurzrockR. PD-L1 expression as a predictive biomarker in cancer immunotherapy. Mol Cancer Ther. (2015) 14:847–56. doi: 10.1158/1535-7163.MCT-14-0983
13
XiaZChenSHeMLiBDengYYiLet al. Editorial: Targeting metabolism to activate T cells and enhance the efficacy of checkpoint blockade immunotherapy in solid tumors. Front Immunol. (2023) 14:1247178. doi: 10.3389/fimmu.2023.1247178
14
ChenJWDhahbiJ. Lung adenocarcinoma and lung squamous cell carcinoma cancer classification, biomarker identification, and gene expression analysis using overlapping feature selection methods. Sci Rep. (2021) 11:13323. doi: 10.1038/s41598-021-92725-8
15
LiuXJiaYShiCKongDWuYZhangTet al. CYP4B1 is a prognostic biomarker and potential therapeutic target in lung adenocarcinoma. PloS One. (2021) 16:e0247020. doi: 10.1371/journal.pone.0247020
16
LiuXSZhouLMYuanLLGaoYKuiXYLiuXYet al. NPM1 is a prognostic biomarker involved in immune infiltration of lung adenocarcinoma and associated with m6A modification and glycolysis. Front Immunol. (2021) 12:724741. doi: 10.3389/fimmu.2021.724741
17
YuYWangZZhengQLiJ. FAM72 serves as a biomarker of poor prognosis in human lung adenocarcinoma. Aging (Albany NY). (2021) 13:8155–76. doi: 10.18632/aging.v13i6
18
ZengJWuZLuoMXuXBaiWXieGet al. Development and validation of an endoplasmic reticulum stress long non-coding RNA signature for the prognosis and immune landscape prediction of patients with lung adenocarcinoma. Front Genet. (2023) 14:1024444. doi: 10.3389/fgene.2023.1024444
19
SunZWangJFanZYangYMengXMaZet al. Investigating the prognostic role of lncRNAs associated with disulfidptosis-related genes in clear cell renal cell carcinoma. J Gene Med. (2024) 26:e3608. doi: 10.1002/jgm.3608
20
LiuYWuZFengYGaoJWangBLianCet al. Integration analysis of single-cell and spatial transcriptomics reveal the cellular heterogeneity landscape in glioblastoma and establish a polygenic risk model. Front Oncol. (2023) 13:1109037. doi: 10.3389/fonc.2023.1109037
21
PanCDengDWeiTWuZZhangBYuanQet al. Metabolomics study identified bile acids as potential biomarkers for gastric cancer: A case control study. Front Endocrinol (Lausanne). (2022) 13:1039786. doi: 10.3389/fendo.2022.1039786
22
FuDZhangBZhangYFengJJiangH. Immunogenomic classification of lung squamous cell carcinoma characterizes tumor immune microenvironment and predicts cancer therapy. Genes Dis. (2023). doi: 10.1016/j.gendis.2023.01.022
23
ZhangYYangZTangYGuoCLinDChengLet al. Hallmark guided identification and characterization of a novel immune-relevant signature for prognostication of recurrence in stage I–III lung adenocarcinoma. Genes Dis. (2023) 10:1657–74. doi: 10.1016/j.gendis.2022.07.005
24
WangYZhouZChenLLiYZhouZChuX. Identification of key genes and biological pathways in lung adenocarcinoma via bioinformatics analysis. Mol Cell Biochem. (2021) 476:931–9. doi: 10.1007/s11010-020-03959-5
25
RuizEKandilEAlhassanSToraihEErramiYElmageedZYAet al. An integrative multi-omics analysis of the molecular links between aging and aggressiveness in thyroid cancers. Aging Dis. (2023) 14:992–1012. doi: 10.14336/AD.2022.1021
26
ChengFHongHYangSWeiY. Individualized network-based drug repositioning infrastructure for precision oncology in the panomics era. Brief Bioinform. (2017) 18:682–97.
27
NasaIKettenbachAN. Coordination of protein kinase and phosphoprotein phosphatase activities in mitosis. Front Cell Dev Biol. (2018) 6:30. doi: 10.3389/fcell.2018.00030
28
KhannaAPimandaJEWestermarckJ. Cancerous inhibitor of protein phosphatase 2A, an emerging human oncoprotein and a potential cancer therapy target. Cancer Res. (2013) 73:6548–53. doi: 10.1158/0008-5472.CAN-13-1994
29
XieWSunYZengYHuLZhiJLingHet al. Comprehensive analysis of PPPCs family reveals the clinical significance of PPP1CA and PPP4C in breast cancer. Bioengineered. (2022) 13:190–205. doi: 10.1080/21655979.2021.2012316
30
HangJLauSYYinRZhuLZhouSYuanXet al. The role of phosphoprotein phosphatases catalytic subunit genes in pancreatic cancer. Biosci Rep. (2021) 41. doi: 10.1042/BSR20203282
31
YeSYangBZhangTWeiWLiZChenJet al. Identification of an immune-related prognostic signature for glioblastoma by comprehensive bioinformatics and experimental analyses. Cells. (2022) 11. doi: 10.3390/cells11193000
32
RajaRWuCBassoyEYRubinoTEJr.UtagawaECMagtibayPMet al. PP4 inhibition sensitizes ovarian cancer to NK cell-mediated cytotoxicity via STAT1 activation and inflammatory signaling. J Immunother Cancer. (2022) 10. doi: 10.1136/jitc-2022-005026
33
LangfelderPHorvathS. WGCNA: an R package for weighted correlation network analysis. BMC Bioinf. (2008) 9:559. doi: 10.1186/1471-2105-9-559
34
OkayamaHKohnoTIshiiYShimadaYShiraishiKIwakawaRet al. Identification of genes upregulated in ALK-positive and EGFR/KRAS/ALK-negative lung adenocarcinomas. Cancer Res. (2012) 72:100–11. doi: 10.1158/0008-5472.CAN-11-1403
35
ZhangJPengGChiHYangJXieXSongGet al. CD8 + T-cell marker genes reveal different immune subtypes of oral lichen planus by integrating single-cell RNA-seq and bulk RNA-sequencing. BMC Oral Health. (2023) 23:464. doi: 10.1186/s12903-023-03138-0
36
ZhangTWongG. Gene expression data analysis using Hellinger correlation in weighted gene co-expression networks (WGCNA). Comput Struct Biotechnol J. (2022) 20:3851–63. doi: 10.1016/j.csbj.2022.07.018
37
ZhangSJiangCJiangLChenHHuangJGaoXet al. Construction of a diagnostic model for hepatitis B-related hepatocellular carcinoma using machine learning and artificial neural networks and revealing the correlation by immunoassay. Tumour Virus Res. (2023) 16:200271. doi: 10.1016/j.tvr.2023.200271
38
MaBQinLSunZWangJTranLJZhangJet al. The single-cell evolution trajectory presented different hypoxia heterogeneity to reveal the carcinogenesis of genes in clear cell renal cell carcinoma: Based on multiple omics and real experimental verification. Environ Toxicol. (2024) 39:869–81. doi: 10.1002/tox.24009
39
ZhangBSunJGuanHGuoHHuangBChenXet al. Integrated single-cell and bulk RNA sequencing revealed the molecular characteristics and prognostic roles of neutrophils in pancreatic cancer. Aging (Albany NY). (2023) 15:9718–42. doi: 10.18632/aging.v15i18
40
ChiHGaoXXiaZYuWYinXPanYet al. FAM family gene prediction model reveals heterogeneity, stemness and immune microenvironment of UCEC. Front Mol Biosci. (2023) 10:1200335. doi: 10.3389/fmolb.2023.1200335
41
GuanHChenXLiuJSunJGuoHJiangYet al. Molecular characteristics and therapeutic implications of Toll-like receptor signaling pathway in melanoma. Sci Rep. (2023) 13:13788. doi: 10.1038/s41598-023-38850-y
42
JiangSYangXLinYLiuYTranLJZhangJet al. Unveiling Anoikis-related genes: A breakthrough in the prognosis of bladder cancer. J Gene Med. (2024) 26:e3651. doi: 10.1002/jgm.3651
43
WuZZengJWuMLiangQLiBHouGet al. Identification and validation of the pyroptosis-related long noncoding rna signature to predict the prognosis of patients with bladder cancer. Med (Baltimore). (2023) 102:e33075. doi: 10.1097/MD.0000000000033075
44
ZhangBLiuJLiHHuangBZhangBSongBet al. Integrated multi-omics identified the novel intratumor microbiome-derived subtypes and signature to predict the outcome, tumor microenvironment heterogeneity, and immunotherapy response for pancreatic cancer patients. Front Pharmacol. (2023) 14:1244752. doi: 10.3389/fphar.2023.1244752
45
SongGPengGZhangJSongBYangJXieXet al. Uncovering the potential role of oxidative stress in the development of periodontitis and establishing a stable diagnostic model via combining single-cell and machine learning analysis. Front Immunol. (2023) 14:1181467. doi: 10.3389/fimmu.2023.1181467
46
ZhuCSunZWangJMengXMaZGuoRet al. Exploring oncogenes for renal clear cell carcinoma based on G protein-coupled receptor-associated genes. Discovery Oncol. (2023) 14:182. doi: 10.1007/s12672-023-00795-z
47
XiaoJLinHLiuBXiaZZhangJJinJ. Decreased S1P and SPHK2 are involved in pancreatic acinar cell injury. Biomark Med. (2019) 13:627–37. doi: 10.2217/bmm-2018-0404
48
ZhangHXiaTXiaZZhouHLiZWangWet al. KIF18A inactivates hepatic stellate cells and alleviates liver fibrosis through the TTC3/Akt/mTOR pathway. Cell Mol Life Sci. (2024) 81:96. doi: 10.1007/s00018-024-05114-5
49
ZhaiXXiaZDuGZhangXXiaTMaDet al. LRP1B suppresses HCC progression through the NCSTN/PI3K/AKT signaling axis and affects doxorubicin resistance. Genes Dis. (2023) 10:2082–96. doi: 10.1016/j.gendis.2022.10.021
50
LiZZhouHXiaZXiaTDuGFranziskaSDet al. HMGA1 augments palbociclib efficacy via PI3K/mTOR signaling in intrahepatic cholangiocarcinoma. Biomark Res. (2023) 11:33. doi: 10.1186/s40364-023-00473-w
51
BrewisNDCohenPT. Protein phosphatase X has been highly conserved during mammalian evolution. Biochim Biophys Acta. (1992) 1171:231–3. doi: 10.1016/0167-4781(92)90129-N
52
WangBZhaoASunLZhongXZhongJWangHet al. Protein phosphatase PP4 is overexpressed in human breast and lung tumors. Cell Res. (2008) 18:974–7. doi: 10.1038/cr.2008.274
53
LiXLiangLHuangLMaXLiDCaiS. High expression of protein phosphatase 4 is associated with the aggressive Malignant behavior of colorectal carcinoma. Mol Cancer. (2015) 14:95. doi: 10.1186/s12943-015-0356-7
54
LiMLiXXuSXuePLiQLuQet al. Protein phosphatase 4 catalytic subunit is overexpressed in glioma and promotes glioma cell proliferation and invasion. Tumour Biol. (2016) 37:11893–901. doi: 10.1007/s13277-016-5054-6
55
WengSWangHChenWKatzMHChatterjeeDLeeJEet al. Overexpression of protein phosphatase 4 correlates with poor prognosis in patients with stage II pancreatic ductal adenocarcinoma. Cancer Epidemiol Biomarkers Prev. (2012) 21:1336–43. doi: 10.1158/1055-9965.EPI-12-0223
56
HuMCTang-OxleyQQiuWRWangYPMihindukulasuriyaKAAfsharRet al. Protein phosphatase X interacts with c-Rel and stimulates c-Rel/nuclear factor kappaB activity. J Biol Chem. (1998) 273:33561–5. doi: 10.1074/jbc.273.50.33561
57
ZhouGBoomerJSTanTH. Protein phosphatase 4 is a positive regulator of hematopoietic progenitor kinase 1. J Biol Chem. (2004) 279:49551–61. doi: 10.1074/jbc.M410317200
58
BertramPGChoiJHCarvalhoJAiWZengCChanTFet al. Tripartite regulation of Gln3p by TOR, Ure2p, and phosphatases. J Biol Chem. (2000) 275:35727–33. doi: 10.1074/jbc.M004235200
59
WangBZhuXXPanLYChenHFShenXY. PP4C facilitates lung cancer proliferation and inhibits apoptosis via activating MAPK/ERK pathway. Pathol Res Pract. (2020) 216:152910. doi: 10.1016/j.prp.2020.152910
60
DunnGPOldLJSchreiberRD. The immunobiology of cancer immunosurveillance and immunoediting. Immunity. (2004) 21:137–48. doi: 10.1016/j.immuni.2004.07.017
61
YajnikVPauldingCSordellaRMcClatcheyAISaitoMWahrerDCet al. DOCK4, a GTPase activator, is disrupted during tumorigenesis. Cell. (2003) 112:673–84. doi: 10.1016/S0092-8674(03)00155-7
62
ThylurRPGowdaRMishraSJunCD. Swiprosin-1: its expression and diverse biological functions. J Cell Biochem. (2018) 119:150–6. doi: 10.1002/jcb.26199
63
YingXLiuBYuanZHuangYChenCJiangXet al. METTL1-m(7) G-EGFR/EFEMP1 axis promotes the bladder cancer development. Clin Transl Med. (2021) 11:e675. doi: 10.1002/ctm2.675
64
ChenZZhuWZhuSSunKLiaoJLiuHet al. METTL1 promotes hepatocarcinogenesis via m(7) G tRNA modification-dependent translation control. Clin Transl Med. (2021) 11:e661. doi: 10.1002/ctm2.661
65
HeJZhouMYinJWanJChuJJiaJet al. METTL3 restrains papillary thyroid cancer progression via m(6)A/c-Rel/IL-8-mediated neutrophil infiltration. Mol Ther. (2021) 29:1821–37. doi: 10.1016/j.ymthe.2021.01.019
Summary
Keywords
lung adenocarcinoma, PPP4C, immunotherapy, prognosis, single-cell
Citation
Wang K, Peng B, Xu R, Lu T, Chang X, Shen Z, Shi J, Li M, Wang C, Zhou X, Xu C, Chang H and Zhang L (2024) Comprehensive analysis of PPP4C’s impact on prognosis, immune microenvironment, and immunotherapy response in lung adenocarcinoma using single-cell sequencing and multi-omics. Front. Immunol. 15:1416632. doi: 10.3389/fimmu.2024.1416632
Received
12 April 2024
Accepted
17 June 2024
Published
04 July 2024
Volume
15 - 2024
Edited by
Yuan Li, Shandong University, China
Reviewed by
Ziwei Yin, Central South University, China
Jing Zhang, University of South Dakota, United States
Zhenyu Wu, First People’s Hospital of Foshan, China
Updates
Copyright
© 2024 Wang, Peng, Xu, Lu, Chang, Shen, Shi, Li, Wang, Zhou, Xu, Chang and Zhang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Linyou Zhang, lyzhang@hrbmu.edu.cn; Hao Chang, doctorhaochang@sohu.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.