A Novel Immune-Related Gene Signature Predicts Prognosis of Lung Adenocarcinoma

Unraveling the copper-death connection: Decoding COVID-19's immune landscape through advanced bioinformatics and machine learning approaches

This study aims to analyze Coronavirus Disease 2019 (COVID-19)-associated copper-death genes using the Gene Expression Omnibus (GEO) dataset and machine learning, exploring their immune microenvironment correlation and underlying mechanisms. Utilizing GEO, we analyzed the GSE217948 dataset with control samples. Differential expression analysis identified 16 differentially expressed copper-death genes, and Cell type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) quantified immune cell infiltration. Gene classification yielded two copper-death clusters, with Weighted Gene Co-expression Network Analysis (WGCNA) identifying key module genes. Machine learning models (random forest, Support Vector Machine (SVM), Generalized Linear Model (GLM), eXtreme Gradient Boosting (XGBoost)) selected 6 feature genes validated by the GSE213313 dataset. Ferredoxin 1 (FDX1) emerged as the top gene, corroborated by Area Under the Curve (AUC) analysis. Gene Set Enrichment Analysis (GSEA) and Gene Set Variation Analysis (GSVA) revealed enriched pathways in T cell receptor, natural killer cytotoxicity, and Peroxisome Proliferator-Activated Receptor (PPAR). We uncovered differentially expressed copper-death genes and immune infiltration differences, notably CD8 T cells and M0 macrophages. Clustering identified modules with potential implications for COVID-19. Machine learning models effectively predicted COVID-19 risk, with FDX1's pivotal role validated. FDX1's high expression was associated with immune pathways, suggesting its role in COVID-19 pathogenesis. This comprehensive approach elucidated COVID-19-related copper-death genes, their immune context, and risk prediction potential. FDX1's connection to immune pathways offers insights into COVID-19 mechanisms and therapy.

Construction of an immune-related gene prognostic model with experimental validation and analysis of immune cell infiltration in lung adenocarcinoma

There is a correlation between tumors and immunity with the degree of immune cell infiltration in tumors being closely related to tumor growth and progression. Therefore, the present study identified immune-related prognostic genes and evaluated the immune infiltration level in lung adenocarcinoma (LUAD). This study performed Kyoto Encyclopedia of Genes and Genomes, Gene Ontology, and Gene Set Enrichment Analysis (GSEA) enrichment analyses on differential immune-associated genes. A risk model was created and validated using six immune-related prognostic genes. Reverse transcription-quantitative PCR was used to assess the prognostic gene expression in non-small cell lung cancer cells. Immune cell infiltration in LUAD was analyzed using the CIBERSORT method. Single sample GSEA was used to compare Tumor Immune Dysfunction and Exclusion (TIDE) scores between high and low-risk groups and to assess the activation of thirteen immune-related pathways. Multifactor Cox proportional hazards model analysis identified six prognostic risk genes (S100A16, FURIN, FGF2, LGR4, TNFRSF11A and VIPR1) to construct a risk model. The survival and receiver operating characteristic curves indicated that patients with higher risk scores had lower overall survival rates. The expression levels of prognostic genes S100A16, FURIN, LGR4, TNFRSF11A and VIPR1 were significantly increased in LUAD. B cells naive, plasma cells, T cells CD4 memory activated, T cells follicular helper, T cells regulatory, NK cells activated, macrophages M1, macrophages M2, and Dendritic cells resting cells showed elevated expression in LUAD. The prognostic genes were differentially associated with individual immune cells. Immune-related function scores, such as those for antigen presenting cell (APC) co-stimulation, APC co-inhibition, check-point, Cytolytic-activity, chemokine receptor, parainflammation, major histocompatibility complex-class-I, type-I-IFN-reponse and T-cell-co-inhibition, were higher in the high-risk group compared with the low-risk group. Furthermore, the TIDE score of the high-risk group was significantly lower than the low-risk group. This immune-related gene prognostic model has the potential to predict the prognosis of LUAD patients, supporting the development of a personalized clinical diagnosis and treatment plan.

Unraveling Differences in Molecular Mechanisms and Immunological Contrasts between Squamous Cell Carcinoma and Adenocarcinoma of the Cervix

This study aims to refine our understanding of the inherent heterogeneity in cervical cancer by exploring differential gene expression profiles, immune cell infiltration dynamics, and implicated signaling pathways in the two predominant histological types of cervix carcinoma, Squamous Cell Carcinoma (SCC) and Adenocarcinoma (ADC). Targeted gene expression data that were previously generated from samples of primary cervical cancer were re-analyzed. The samples were grouped based on their histopathology, comparing SCC to ADC. Each tumor in the study was confirmed to be high risk human papilloma virus (hrHPV) positive. A total of 21 cervical cancer samples were included, with 11 cases of SCC and 10 of ADC. Data analysis revealed a total of 26 differentially expressed genes, with 19 genes being overexpressed in SCC compared to ADC (Benjamini-Hochberg (BH)-adjusted p-value < 0.05). Importantly, the immune checkpoint markers CD274 and CTLA4 demonstrated significantly higher expression in SCC compared to ADC. In addition, SCC showed a higher infiltration of immune cells, including B and T cells, and cytotoxic cells. Higher activation of a variety of pathways was found in SCC samples including cytotoxicity, interferon signaling, metabolic stress, lymphoid compartment, hypoxia, PI3k-AKT, hedgehog signaling and Notch signaling pathways. Our findings show distinctive gene expression patterns, signaling pathway activations, and trends in immune cell infiltration between SCC and ADC in cervical cancer. This study underscores the heterogeneity within primary cervical cancer, emphasizing the potential benefits of subdividing these tumours based on histological and molecular differences.

A robust six-gene prognostic signature based on two prognostic subtypes constructed by chromatin regulators is correlated with immunological features and therapeutic response in lung adenocarcinoma

Accumulating evidence has demonstrated that chromatin regulators (CRs) regulate immune cell infiltration and are correlated with prognoses of patients in some cancers. However, the immunological and prognostic roles of CRs in lung adenocarcinoma (LUAD) are still unclear. Here, we systematically revealed the correlations of CRs with immunological features and the survival in LUAD patients based on a cohort of gene expression datasets from the public TCGA and GEO databases and real RNA-seq data by an integrative analysis using a comprehensive bioinformatics method. Totals of 160 differentially expressed CRs (DECRs) were identified between LUAD and normal lung tissues, and two molecular prognostic subtypes (MPSs) were constructed and evaluated based on 27 prognostic DECRs using five independent datasets (p =0.016, <0.0001, =0.008, =0.00038 and =0.00055, respectively). Six differentially expressed genes (DEGs) (CENPK, ANGPTL4, CCL20, CPS1, GJB3, TPSB2) between two MPSs had the most important prognostic feature and a six-gene prognostic model was established. LUAD patients in the low-risk subgroup showed a higher overall survival (OS) rate than those in the high-risk subgroup in nine independent datasets (p <0.0001, =0.021, =0.016, =0.0099, <0.0001, =0.0045, <0.0001, =0.0038 and =0.00013, respectively). Six-gene prognostic signature had the highest concordance index of 0.673 compared with 19 reported prognostic signatures. The risk score was significantly correlated with immunological features and activities of oncogenic signaling pathways. LUAD patients in the low-risk subgroup benefited more from immunotherapy and were less sensitive to conventional chemotherapy agents. This study provides novel insights into the prognostic and immunological roles of CRs in LUAD.

Association of Immune Semaphorins with COVID-19 Severity and Outcomes

Semaphorins have recently been recognized as crucial modulators of immune responses. In the pathogenesis of COVID-19, the activation of immune responses is the key factor in the development of severe disease. This study aimed to determine the association of serum semaphorin concentrations with COVID-19 severity and outcomes. Serum semaphorin concentrations (SEMA3A, -3C, -3F, -4D, -7A) were measured in 80 hospitalized adult patients with COVID-19 (moderate (n = 24), severe (n = 32), critical, (n = 24)) and 40 healthy controls. While SEMA3C, SEMA3F and SEMA7A serum concentrations were significantly higher in patients with COVID-19, SEMA3A was significantly lower. Furthermore, SEMA3A and SEMA3C decreased with COVID-19 severity, while SEMA3F and SEMA7A increased. SEMA4D showed no correlation with disease severity. Serum semaphorin levels show better predictive values than CRP, IL-6 and LDH for differentiating critical from moderate/severe COVID-19. SEMA3F and SEMA7A serum concentrations were associated with the time to recovery, requirement of invasive mechanical ventilation, development of pulmonary thrombosis and nosocomial infections, as well as with in-hospital mortality. In conclusion, we provide the first evidence that SEMA3A, SEMA3C, SEMA3F and SEMA7A can be considered as new biomarkers of COVID-19 severity.

Identification of age- and immune-related gene signatures for clinical outcome prediction in lung adenocarcinoma

BACKGROUND

The understanding of the factors causing decreased overall survival (OS) in older patients compared to younger patients in lung adenocarcinoma (LUAD) remains.

METHODS

Gene expression profiles of LUAD were obtained from publicly available databases by Kaplan-Meier analysis was performed to determine whether age was associated with patient OS. The immune cell composition in the tumor microenvironment (TME) was evaluated using CIBERSORT. The fraction of stromal and immune cells in tumor samples were also using assessed using multiple tools including ESTIMATE, EPIC, and TIMER. Differentially expressed genes (DEGs) from the RNA-Seq data that were associated with age and immune cell composition were identified using the R package DEGseq. A 22-gene signature composed of DEGs associated with age and immune cell composition that predicted OS were constructed using Least Absolute Shrinkage and Selection Operator (LASSO).

RESULTS

In The Cancer Genome Atlas (TCGA)-LUAD dataset, we found that younger patients (≤70) had a significant better OS compared to older patients (>70). In addition, older patients had significantly higher expression of immune checkpoint proteins including inhibitory T cell receptors and their ligands. Moreover, analyses using multiple bioinformatics tools showed increased immune infiltration, including CD4+ T cells, in older patients compared to younger patients. We identified a panel of genes differentially expressed between patients >70 years compared to those ≤70 years, as well as between patients with high or low immune scores and selected 84 common genes to construct a prognostic gene signature. A risk score calculated based on 22 genes selected by LASSO predicted 1, 3, and 5-year OS, with an area under the curve (AUC) of 0.72, 0.72, 0.69, receptively, in TCGA-LUAD dataset and an independent validation dataset available from the European Genome-phenome Archive (EGA).

CONCLUSION

Our results demonstrate that age contributes to OS of LUAD patients atleast in part through its association with immune infiltration in the TME.

An epithelial-mesenchymal transition-related mRNA signature associated with the prognosis, immune infiltration and therapeutic response of colon adenocarcinoma

Background: Epithelial-mesenchymal transition (EMT) is closely associated with cancer cell metastasis. Colon adenocarcinoma (COAD) is one of the most common malignancies in the world, and its metastasis leading to poor prognosis remains a challenge for clinicians. The purpose of this study was to explore the prognostic value of EMT-related genes (EMTRGs) by bioinformatics analysis and to develop a new EMTRGs prognostic signature for COAD. Methods: The TCGA-COAD dataset was downloaded from the TCGA portal as the training cohort, and the GSE17538 and GSE29621 datasets were obtained from the GEO database as the validation cohort. The best EMTRGs prognostic signature was constructed by differential expression analysis, Cox, and LASSO regression analysis. Gene set enrichment analysis (GSEA) is used to reveal pathways that are enriched in high-risk and low-risk groups. Differences in tumor immune cell levels were analyzed using microenvironmental cell population counter and single sample gene set enrichment analysis. Subclass mapping analysis and Genomics of Drug Sensitivity in Cancer were applied for prediction of immunotherapy response and chemotherapy response, respectively. Results: A total of 77 differentially expressed EMTRGs were identified in the TCGA-COAD cohort, and they were significantly associated with functions and pathways related to cancer cell metastasis, proliferation, and apoptosis. We constructed EMTRGs prognostic signature with COMP, MYL9, PCOLCE2, SCG2, and TIMP1 as new COAD prognostic biomarkers. The high-risk group had a poorer prognosis with enhanced immune cell infiltration. The GSEA demonstrated that the high-risk group was involved in "ECM Receptor Interaction," "WNT Signaling Pathway" and "Colorectal Cancer." Furthermore, patients with high risk scores may respond to anti-CTLA4 therapy and may be more resistant to targeted therapy agents BI 2536 and ABT-888. Conclusion: Together, we developed a new EMTRGs prognostic signature that can be an independent prognostic factor for COAD. This study has guiding implications for individualized counseling and treatment of COAD patients.

Prognosis and personalized treatment prediction in lung adenocarcinoma: An in silico and in vitro strategy adopting cuproptosis related lncRNA towards precision oncology

Background: There is a rapid increase in lung adenocarcinomas (LUAD), and studies suggest associations between cuproptosis and the occurrence of various types of tumors. However, it remains unclear whether cuproptosis plays a role in LUAD prognosis. Methods: Dataset of the TCGA-LUAD was treated as training cohort, while validation cohort consisted of the merged datasets of the GSE29013, GSE30219, GSE31210, GSE37745, and GSE50081. Ten studied cuproptosis-related genes (CRG) were used to generated CRG clusters and CRG cluster-related differential expressed gene (CRG-DEG) clusters. The differently expressed lncRNA that with prognosis ability between the CRG-DEG clusters were put into a LASSO regression for cuproptosis-related lncRNA signature (CRLncSig). Kaplan-Meier estimator, Cox model, receiver operating characteristic (ROC), time-dependent AUC (tAUC), principal component analysis (PCA), and nomogram predictor were further deployed to confirm the model's accuracy. We examined the model's connections with other forms of regulated cell death, including apoptosis, necroptosis, pyroptosis, and ferroptosis. The immunotherapy ability of the signature was demonstrated by applying eight mainstream immunoinformatic algorithms, TMB, TIDE, and immune checkpoints. We evaluated the potential drugs for high risk CRLncSig LUADs. Real-time PCR in human LUAD tissues were performed to verify the CRLncSig expression pattern, and the signature's pan-cancer's ability was also assessed. Results: A nine-lncRNA signature, CRLncSig, was built and demonstrated owning prognostic power by applied to the validation cohort. Each of the signature genes was confirmed differentially expressed in the real world by real-time PCR. The CRLncSig correlated with 2,469/3,681 (67.07%) apoptosis-related genes, 13/20 (65.00%) necroptosis-related genes, 35/50 (70.00%) pyroptosis-related genes, and 238/380 (62.63%) ferroptosis-related genes. Immunotherapy analysis suggested that CRLncSig correlated with immune status, and checkpoints, KIR2DL3, IL10, IL2, CD40LG, SELP, BTLA, and CD28, were linked closely to our signature and were potentially suitable for LUAD immunotherapy targets. For those high-risk patients, we found three agents, gemcitabine, daunorubicin, and nobiletin. Finally, we found some of the CRLncSig lncRNAs potentially play a vital role in some types of cancer and need more attention in further studies. Conclusion: The results of this study suggest our cuproptosis-related CRLncSig can help to determine the outcome of LUAD and the effectiveness of immunotherapy, as well as help to better select targets and therapeutic agents.

Single-cell RNA sequencing analysis to explore immune cell heterogeneity and novel biomarkers for the prognosis of lung adenocarcinoma

Background: Single-cell RNA sequencing is necessary to understand tumor heterogeneity, and the cell type heterogeneity of lung adenocarcinoma (LUAD) has not been fully studied.

Method: We first reduced the dimensionality of the GSE149655 single-cell data. Then, we statistically analysed the subpopulations obtained by cell annotation to find the subpopulations highly enriched in tumor tissues. Monocle was used to predict the development trajectory of five subpopulations; beam was used to find the regulatory genes of five branches; qval was used to screen the key genes; and cellchart was used to analyse cell communication. Next, we used the differentially expressed genes of TCGA-LUAD to screen for overlapping genes and established a prognostic risk model through univariate and multivariate analyses. To identify the independence of the model in clinical application, univariate and multivariate Cox regression were used to analyse the relevant HR, 95% CI of HR and p value. Finally, the novel biomarker genes were verified by qPCR and immunohistochemistry.

Results: The single-cell dataset GSE149655 was subjected to quality control, filtration and dimensionality reduction. Finally, 23 subpopulations were screened, and 11-cell subgroups were annotated in 23 subpopulations. Through the statistical analysis of 11 subgroups, five important subgroups were selected, including lung epithelial cells, macrophages, neuroendocrine cells, secret cells and T cells. From the analysis of cell trajectory and cell communication, it is found that the interaction of five subpopulations is very complex and that the communication between them is dense. We believe that these five subpopulations play a very important role in the occurrence and development of LUAD. Downloading the TCGA data, we screened the marker genes of these five subpopulations, which are also the differentially expressed genes in tumorigenesis, with a total of 462 genes, and constructed 10 gene prognostic risk models based on related genes. The 10-gene signature has strong robustness and can achieve stable prediction efficiency in datasets from different platforms. Two new molecular markers related to LUAD, HLA-DRB5 and CCDC50, were verified by qPCR and immunohistochemistry. The results showed that HLA-DRB5 expression was negatively correlated with the risk of LUAD, and CCDC50 expression was positively correlated with the risk of LUAD.

Conclusion: Therefore, we identified a prognostic risk model including CCL20, CP, HLA-DRB5, RHOV, CYP4B1, BASP1, ACSL4, GNG7, CCDC50 and SPATS2 as risk biomarkers and verified their predictive value for the prognosis of LUAD, which could serve as a new therapeutic target.