Comprehensive analysis of co-expressed genes with TDP-43: prognostic and therapeutic potential in lung adenocarcinoma

Abstract

Abstract Background Transactivating DNA-binding protein 43 (TDP-43) is intimately associated with tumorigenesis and progression by regulating mRNA splicing, transport, stability, and non-coding RNA molecules. The exact role of TDP-43 in lung adenocarcinoma (LUAD) has not yet been fully elucidated, despite extensive research on its function in various cancer types. An imperative aspect of comprehending the underlying biological characteristics associated with TDP-43 involves investigating the genes that are co-expressed with this protein. This study assesses the prognostic significance of these co-expressed genes in LUAD and subsequently explores potential therapeutic strategies based on these findings. Methods Transcriptomic and clinical data pertaining to LUAD were retrieved from open-access databases to establish an association between mRNA expression profiles and the presence of TDP-43. A risk-prognosis model was developed to compare patient survival rates across various groups, and its accuracy was also assessed. Additionally, differences in tumor stemness, mutational profiles, tumor microenvironment (TME) characteristics, immune checkpoints, and immune cell infiltration were analyzed in the different groups. Moreover, the study entailed predicting the potential response to immunotherapy as well as the sensitivity to commonly employed chemotherapeutic agents and targeted drugs for each distinct group. Results The TDP-43 Co-expressed Gene Risk Score (TCGRS) model was constructed utilizing four genes: Kinesin Family Member 20A (KIF20A), WD Repeat Domain 4 (WDR4), Proline Rich 11 (PRR11), and Glia Maturation Factor Gamma (GMFG). The value of this model in predicting LUAD patient survival is effectively illustrated by both the Kaplan–Meier (K–M) survival curve and the area under the receiver operating characteristic curve (AUC-ROC). The Gene Set Enrichment Analysis (GSEA) revealed that the high TCGRS group was primarily enriched in biological pathways and functions linked to DNA replication and cell cycle; the low TCGRS group showed primary enrichment in immune-related pathways and functions. The high and low TCGRS groups showed differences in tumor stemness, mutational burden, TME, immune infiltration level, and immune checkpoints. The predictions analysis of immunotherapy indicates that the Tumor Immune Dysfunction and Exclusion (TIDE) score (p < 0.001) and non-response rate (74% vs. 51%, p < 0.001) in the high TCGRS group are higher than those in the low TCGRS group. The Immune Phenotype Score (IPS) in the high TCGRS group is lower than in the low TCGRS group (p < 0.001). The drug sensitivity analysis revealed that the half-maximal inhibitory concentration (IC50) values for cisplatin, docetaxel, doxorubicin, etoposide, gemcitabine, paclitaxel, vincristine, erlotinib, and gefitinib (all p < 0.01) in the high TCGRS group are lower than those in the low TCGRS group. Conclusions The TCGRS derived from the model exhibits a reliable biomarker for evaluating both prognosis and treatment effectiveness among patients with LUAD. This study is anticipated to offer valuable insights into developing effective treatment strategies for this patient population. It is believed that this study is anticipated to contribute significantly to clinical diagnostics, the development of therapeutic drugs, and the enhancement of patient care.