Comprehensive characterization of stemness-related lncRNAs in triple-negative breast cancer identified a novel prognostic signature related to treatment outcomes, immune landscape analysis and therapeutic guidance: a silico analysis with in vivo experiments.

BACKGROUND: Cancer stem cells (CSCs) and long non-coding RNAs (lncRNAs) are known to play a crucial role in the growth, migration, recurrence, and drug resistance of tumor cells, particularly in triple-negative breast cancer (TNBC). This study aims to investigate stemness-related lncRNAs (SRlncRNAs) as potential prognostic indicators for TNBC patients. METHODS: Utilizing RNA sequencing data and corresponding clinical information from the TCGA database, and employing Weighted Gene Co-expression Network Analysis (WGCNA) on TNBC mRNAsi sourced from an online database, stemness-related genes (SRGs) and SRlncRNAs were identified. A prognostic model was developed using univariate Cox and LASSO-Cox analysis based on SRlncRNAs. The performance of the model was evaluated using Kaplan-Meier analysis, ROC curves, and ROC-AUC. Additionally, the study delved into the underlying signaling pathways and immune status associated with the divergent prognoses of TNBC patients. RESULTS: The research identified a signature of six SRlncRNAs (AC245100.6, LINC02511, AC092431.1, FRGCA, EMSLR, and MIR193BHG) for TNBC. Risk scores derived from this signature were found to correlate with the abundance of plasma cells. Furthermore, the nominated chemotherapy drugs for TNBC exhibited considerable variability between different risk score groups. RT-qPCR validation confirmed abnormal expression patterns of these SRlncRNAs in TNBC stem cells, affirming the potential of the SRlncRNAs signature as a prognostic biomarker. CONCLUSION: The identified signature not only demonstrates predictive power in terms of patient outcomes but also provides insights into the underlying biology, signaling pathways, and immune status associated with TNBC prognosis. The findings suggest the possibility of guiding personalized treatments, including immune checkpoint gene therapy and chemotherapy strategies, based on the risk scores derived from the SRlncRNA signature. Overall, this research contributes valuable knowledge towards advancing precision medicine in the context of TNBC.

Zinc metabolism and its role in immunity status in subjects with trisomy 21: chromosomal dosage effect.

Introduction: Trisomy 21 (T21), which causes Down syndrome (DS), is the most common chromosomal aneuploidy in humankind and includes different clinical comorbidities, among which the alteration of the immune system has a heavy impact on patient's lives. A molecule with an important role in immune response is zinc and it is known that its concentration is significantly lower in children with T21. Different hypotheses were made about this metabolic alteration and one of the reasons might be the overexpression of superoxide dismutase 1 (SOD1) gene, as zinc is part of the SOD1 active enzymatic center. Methods: The aim of our work is to explore if there is a linear correlation between zinc level and immune cell levels measured in a total of 217 blood samples from subjects with T21. Furthermore, transcriptome map analyses were performed using Transcriptome Mapper (TRAM) software to investigate whether a difference in gene expression is detectable between subjects with T21 and euploid control group in tissues and cells involved in the immune response such as lymphoblastoid cells, thymus and white blood cells. Results: Our results have confirmed the literature data stating that the blood zinc level in subjects with T21 is lower compared to the general population; in addition, we report that the T21/control zinc concentration ratio is 2:3, consistent with a chromosomal dosage effect due to the presence of three copies of chromosome 21. The transcriptome map analyses showed an alteration of some gene's expression which might explain low levels of zinc in the blood. Discussion: Our data suggest that zinc level is not associated with the levels of immunity cells or proteins analyzed themselves and rather the main role of this ion might be played in altering immune cell function.

AIOLOS-Associated Inborn Errors of Immunity.

AIOLOS, encoded by the IKZF3 gene, belongs to the Ikaros zinc finger transcription factor family and plays a pivotal role in regulating lymphocyte development. Recently, heterozygous missense loss-of-function variants within the DNA-binding domain of the IKZF3 gene (G159R, N160S, and G191R) have been identified in patients with inborn errors of immunity (IEI). Additionally, a missense and a truncating variant (E82K and Q402X) leading to the AIOLOS haploinsufficiency have been documented. The majority of individuals with AIOLOS-associated IEI manifest recurrent sinopulmonary infections, as well as various bacterial and viral infections. The patients carrying the AIOLOSN160S variant exhibit severe immunodeficient phenotypes. In contrast, patients harboring AIOLOS haploinsufficient variants predominantly present with clinical phenotypes associated with immune dysregulation. A varying degree of B-lymphopenia and hypoimmunoglobulinemia was noted in approximately half of the patients. Mouse models of AIOLOSG159R and AIOLOSN160S variants (AiolosG158R and AiolosN159S in mice, respectively) recapitulated most of the immune abnormalities observed in the patients. Among these models, AiolosG158R mice prominently exhibited defects in early B cell differentiation resulting from mutant Aiolos interfering with Ikaros function through heterodimer formation. In contrast, AiolosN159S mice did not manifest early B cell differentiation defects. However, they displayed a distinct immune abnormality characterized by impaired induction of CD62L expression in lymphocytes, which is likely attributable to dysfunction of Ikaros, leading to defective lymphocyte homing to lymph nodes. Considering the diverse clinical phenotypes observed in the reported cases and the distinct molecular pathogenesis associated with each variant, further studies with more patients with AIOLOS-associated IEI would contribute to a better understanding of the clinical spectrum and underlying molecular mechanisms associated with this disorder.

Shaping immunity: The influence of natural selection on population immune diversity.

Humans exhibit considerable variability in their immune responses to the same immune challenges. Such variation is widespread and affects individual and population-level susceptibility to infectious diseases and immune disorders. Although the factors influencing immune response diversity are partially understood, what mechanisms lead to the wide range of immune traits in healthy individuals remain largely unexplained. Here, we discuss the role that natural selection has played in driving phenotypic differences in immune responses across populations and present-day susceptibility to immune-related disorders. Further, we touch on future directions in the field of immunogenomics, highlighting the value of expanding this work to human populations globally, the utility of modeling the immune response as a dynamic process, and the importance of considering the potential polygenic nature of natural selection. Identifying loci acted upon by evolution may further pinpoint variants critically involved in disease etiology, and designing studies to capture these effects will enrich our understanding of the genetic contributions to immunity and immune dysregulation.

Epigenetic control of microglial immune responses.

Microglia, the major population of brain-resident macrophages, are now recognized as a heterogeneous population comprising several cell subtypes with different (so far mostly supposed) functions in health and disease. A number of studies have performed molecular characterization of these different microglial activation states over the last years making use of "omics" technologies, that is transcriptomics, proteomics and, less frequently, epigenomics profiling. These approaches offer the possibility to identify disease mechanisms, discover novel diagnostic biomarkers, and develop new therapeutic strategies. Here, we focus on epigenetic profiling as a means to understand microglial immune responses beyond what other omics methods can offer, that is, revealing past and present molecular responses, gene regulatory networks and potential future response trajectories, and defining cell subtype-specific disease relevance through mapping non-coding genetic variants. We review the current knowledge in the field regarding epigenetic regulation of microglial identity and function, provide an exemplary analysis that demonstrates the advantages of performing joint transcriptomic and epigenomic profiling of single microglial cells and discuss how comprehensive epigenetic analyses may enhance our understanding of microglial pathophysiology.

RNA Metabolism Governs Immune Function and Response.

Inflammation is a complex process that protects our body from various insults such as infection, injury, and stress. Proper inflammation is beneficial to eliminate the insults and maintain organ homeostasis, however, it can become detrimental if uncontrolled. To tightly regulate inflammation, post-transcriptional mechanisms governing RNA metabolism play a crucial role in monitoring the expression of immune-related genes, such as tumor necrosis factor (TNF) and interleukin-6 (IL-6). These mechanisms involve the coordinated action of various RNA-binding proteins (RBPs), including the Regnase family, Roquin, and RNA methyltransferases, which are responsible for mRNA decay and/or translation regulation. The collaborative efforts of these RBPs are essential in preventing aberrant immune response activation and consequently safeguarding against inflammatory and autoimmune diseases. This review provides an overview of recent advancements in our understanding of post-transcriptional regulation within the immune system and explores the specific roles of individual RBPs in RNA metabolism and regulation.

Aberrant H3K4me3 modification of immune response genes in CD4+ T cells of patients with systemic lupus erythematosus.

BACKGROUND: Increasing evidence has highlighted the significant role of histone modifications in pathogenesis of systemic lupus erythematosus (SLE). However, few studies have comprehensively analyzed trimethylation of histone H3 lysine 4 (H3K4me3) features at specific immune gene loci in SLE patients. METHODS: We conducted H3K4me3 chromatin immunoprecipitation sequencing (ChIP-seq) on CD4+ T cells from SLE patients and healthy controls (HC). Differential H3K4me3 peaks were identified, followed by enrichment analysis. We integrated online RNA-seq and DNA methylation datasets to explore the relationship between H3K4me3 modification, DNA methylation and gene expression. We validated several upregulated peak regions by ChIP-qPCR and confirmed their impact on gene expression using RT-qPCR. Finally, we investigated the impact of H3K4 methyltransferases KMT2A on the expression of immune response genes. RESULTS: we identified 147 downregulated and 2701 upregulated H3K4me3 peaks in CD4+ T cells of SLE. The upregulated peaks primarily classified as gained peaks and enriched in immune response genes such as FCGR2A, C5AR1, SERPING1 and OASL. Genes with upregulated H3K4me3 and downregulated DNA methylations in the promoter were highly expressed in SLE patients. These genes, including OAS1, IFI27 and IFI44L, were enriched in immune response pathways. The IFI44L locus also showed increased H3K27ac modification, chromatin accessibility and chromatin interactions in SLE. Moreover, knockdown of KMT2A can downregulate the expression of immune response genes in T cells. CONCLUSION: Our study uncovers dysregulated H3K4me3 modification patterns in immune response genes loci, which also exhibit downregulated DNA methylation and higher mRNA expression in CD4+ T cells of SLE patients.

Downregulation of CPSF6 leads to global mRNA 3' UTR shortening and enhanced antiviral immune responses.

Alternative polyadenylation (APA) is a widespread mechanism of gene regulation that generates mRNA isoforms with alternative 3' untranslated regions (3' UTRs). Our previous study has revealed the global 3' UTR shortening of host mRNAs through APA upon viral infection. However, how the dynamic changes in the APA landscape occur upon viral infection remains largely unknown. Here we further found that, the reduced protein abundance of CPSF6, one of the core 3' processing factors, promotes the usage of proximal poly(A) sites (pPASs) of many immune related genes in macrophages and fibroblasts upon viral infection. Shortening of the 3' UTR of these transcripts may improve their mRNA stability and translation efficiency, leading to the promotion of type I IFN (IFN-I) signalling-based antiviral immune responses. In addition, dysregulated expression of CPSF6 is also observed in many immune related physiological and pathological conditions, especially in various infections and cancers. Thus, the global APA dynamics of immune genes regulated by CPSF6, can fine-tune the antiviral response as well as the responses to other cellular stresses to maintain the tissue homeostasis, which may represent a novel regulatory mechanism for antiviral immunity.

Role of sex in immune response and epigenetic mechanisms.

The functioning of the human immune system is highly dependent on the sex of the individual, which comes by virtue of sex chromosomes and hormonal differences. Epigenetic mechanisms such as X chromosome inactivation, mosaicism, skewing, and dimorphism in X chromosome genes and Y chromosome regulatory genes create a sex-based variance in the immune response between males and females. This leads to differential susceptibility in immune-related disorders like infections, autoimmunity, and malignancies. Various naturally available immunomodulators are also available which target immune pathways containing X chromosome genes.

A2BAR Antagonism Decreases the Glomerular Expression and Secretion of Chemoattractants for Monocytes and the Pro-Fibrotic M2 Macrophages Polarization during Diabetic Nephropathy.

Some chemoattractants and leukocytes such as M1 and M2 macrophages are known to be involved in the development of glomerulosclerosis during diabetic nephropathy (DN). In the course of diabetes, an altered and defective cellular metabolism leads to the increase in adenosine levels, and thus to changes in the polarity (M1/M2) of macrophages. MRS1754, a selective antagonist of the A2B adenosine receptor (A2BAR), attenuated glomerulosclerosis and decreased macrophage-myofibroblast transition in DN rats. Therefore, we aimed to investigate the effect of MRS1754 on the glomerular expression/secretion of chemoattractants, the intraglomerular infiltration of leukocytes, and macrophage polarity in DN rats. Kidneys/glomeruli of non-diabetic, DN, and MRS1754-treated DN rats were processed for transcriptomic analysis, immunohistopathology, ELISA, and in vitro macrophage migration assays. The transcriptomic analysis identified an upregulation of transcripts and pathways related to the immune system in the glomeruli of DN rats, which was attenuated using MRS1754. The antagonism of the A2BAR decreased glomerular expression/secretion of chemoattractants (CCL2, CCL3, CCL6, and CCL21), the infiltration of macrophages, and their polarization to M2 in DN rats. The in vitro macrophages migration induced by conditioned-medium of DN glomeruli was significantly decreased using neutralizing antibodies against CCL2, CCL3, and CCL21. We concluded that the pharmacological blockade of the A2BAR decreases the transcriptional expression of genes/pathways related to the immune response, protein expression/secretion of chemoattractants, as well as the infiltration of macrophages and their polarization toward the M2 phenotype in the glomeruli of DN rats, suggesting a new mechanism implicated in the antifibrotic effect of MRS1754.

Immune responses of the Asian onion moth, Acrolepiopsis sapporensis, and their genetic factors from RNA-Seq analysis.

A nonmodel insect, Acrolepiopsis sapporensis, has been analyzed in immune responses. The total hemocytes in the fifth instar larvae were 2.33 × 106 cells/mL. These hemocytes comprised at least five different types and different relative ratios: 47% granulocytes, 26% plasmatocytes, 11% oenocytoid, 8% prohemocytes, and 5% spherulocytes. Upon bacterial challenge, some of the hemocytes exhibited typical hemocyte-spreading behaviors, such as focal adhesion, and filopodial and lamellipodial cytoplasmic extensions. The hemocyte behaviors induced cellular immune responses demonstrated by nodule formation. In addition, the plasma collected from the immune-challenged larvae exhibited humoral immune responses by bacterial growth inhibition along with enhanced phenoloxidase enzyme activity. These cellular and humoral immune responses were further analyzed by determining the immune-associated genes from a transcriptome generated by RNA-Seq. A total of about 12 Gb sequences led to about 218,116 contigs, which were predicted to encode about 46,808 genes. Comparative expression analysis showed 8392 uniquely expressed genes in the immune-challenged larvae. Differentially expressed gene (DEG) analysis among the commonly expressed genes indicated that 782 genes were upregulated and 548 genes were downregulated in the expressions after bacterial challenge. These immune-associated genes included pattern recognition receptors, immune mediation/signaling genes, and various immune effectors. Specifically, the genetic components of the Toll, IMD, and JAK/STAT immune signaling pathways were included in the DEG database. These results demonstrate the immune responses of A. sapporensis larvae and suggest the genes associated with the immune responses in this nonmodel insect.

A pangenome reference of 36 Chinese populations.

Human genomics is witnessing an ongoing paradigm shift from a single reference sequence to a pangenome form, but populations of Asian ancestry are underrepresented. Here we present data from the first phase of the Chinese Pangenome Consortium, including a collection of 116 high-quality and haplotype-phased de novo assemblies based on 58 core samples representing 36 minority Chinese ethnic groups. With an average 30.65× high-fidelity long-read sequence coverage, an average contiguity N50 of more than 35.63 megabases and an average total size of 3.01 gigabases, the CPC core assemblies add 189 million base pairs of euchromatic polymorphic sequences and 1,367 protein-coding gene duplications to GRCh38. We identified 15.9 million small variants and 78,072 structural variants, of which 5.9 million small variants and 34,223 structural variants were not reported in a recently released pangenome reference1. The Chinese Pangenome Consortium data demonstrate a remarkable increase in the discovery of novel and missing sequences when individuals are included from underrepresented minority ethnic groups. The missing reference sequences were enriched with archaic-derived alleles and genes that confer essential functions related to keratinization, response to ultraviolet radiation, DNA repair, immunological responses and lifespan, implying great potential for shedding new light on human evolution and recovering missing heritability in complex disease mapping.

Spleen Transcriptome Profiling Reveals Divergent Immune Responses to LPS and Poly (I:C) Challenge in the Yellow Drum (Nibea albiflora).

The yellow drum (Nibea albiflora) is a marine teleost fish with strong disease resistance, yet the understanding of its immune response and key functional genes is fragmented. Here, RNA-Seq was used to investigate the regulation pathways and genes involved in the immune response to infection with lipopolysaccharide (LPS) and polyinosinic-polycytidylic acid (poly (I:C)) on the spleen of the yellow drum. There were fewer differentially expressed genes (DEGs) in the LPS-infected treatment group at either 6 or 48 h. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that these DEGs were mainly significantly enriched in c5-branching dibasic acid metabolic and complement and coagulation cascades pathways. The yellow drum responded more strongly to poly (I:C) infection, with 185 and 521 DEGs obtained under 6 and 48 h treatments, respectively. These DEGs were significantly enriched in the Toll-like receptor signaling pathway, RIG-I-like receptor signaling pathway, Jak-STAT signaling pathway, NOD-like signaling pathway, and cytokine-cytokine receptor interaction. The key functional genes in these pathways played important roles in the immune response and maintenance of immune system homeostasis in the yellow drum. Weighted gene co-expression network analysis (WGCNA) revealed several important hub genes. Although the functions of some genes have not been confirmed, our study still provides significant information for further investigation of the immune system of the yellow drum.

Immune-check blocking combination multiple cytokines shown curative potential in mice tumor model.

OBJECTIVE: In order to ensure the stable transcription of target genes, we constructed a eukaryotic high expression vector carrying an immune-check inhibitor PD-1v and a variety of cytokines, and studied their effects on activating immune response to inhibit tumor growth. METHODS: A novel eukaryotic expression plasmid vector named pT7AMPCE containing T7RNA polymerase, T7 promoter, internal ribosome entry site (IRES), and poly A tailing signal was constructed by T4 DNA ligase, on which homologous recombination was used to clone and construct the vector carrying PD-1v, IL-2/15, IL-12, GM-CSF, and GFP. In vitro transfection of CT26 cells was performed, and the protein expression of PD-1v, IL-12 and GM-CSF was detected by Western blot and ELISA after 48 h. Mice were subcutaneously inoculated with CT26-IRFP tumor cells in the rib abdomen, and the tumor tissues were injected with PD-1v, IL-2/15, IL-12, and GM-CSF recombinant plasmids for treatment during the experimental period. The efficacy of the treatment was evaluated by assay tumor size and survival time of tumor-bearing mice during the experiment. Expression levels of IFN-γ, TNF, IL-4, IL-2, and IL-5 in mouse blood were measured using the CBA method. Tumor tissues were extracted and immune cell infiltration in tumor tissues was detected by HE staining and the IHC method. RESULTS: The recombinant plasmids carrying PD-1v, IL-2/15, IL-12, and GM-CSF were successfully constructed, and the Western blot and ELISA results showed that PD-1v, IL-12, and GM-CSF were expressed in the supernatant of CT26 cells 48 h after in vitro cell transfection. The combined application of PD-1v, IL-2/15, IL-12, and GM-CSF recombinant plasmids significantly inhibited tumor growth in mice, and the tumor growth rate was significantly lower than that in the blank control group and GFP plasmid control group (p < 0.05). Cytometric bead array data suggested that the combination of PD-1v and various cytokines can effectively activate immune cells. HE and IHC analysis revealed plenty of immune cell infiltrates in the tumor tissue, and a large proportion of tumor cells showed the necrotic phenotype in the combination treatment group. CONCLUSION: The combination of immune check blockade and multiple cytokine therapy can significantly activate the body's immune response and inhibit tumor growth.

Chicken adaptive response to nutrient density: immune function change revealed by transcriptomic analysis of spleen.

Feed accounts for the largest portion (65-70%) of poultry production costs. The feed formulation is generally improved to efficiently meet the nutritional needs of chickens by reducing the proportion of crude protein (CP) and metabolizable energy (ME) levels in the diet. Although many studies have investigated the production performance during dietary restriction, there is a lack of research on the mechanisms by which immune cell function is altered. This study examined the effects of ME and CP restriction in the chicken diet on serum immunoglobulins and expression of immune function genes in spleen. Changes in serum immunoglobulins and immune-related gene expression were analyzed in 216 YS-909 broilers fed with 9 different dietary treatments, including experimental treatment diets containing low, standard, and high levels of ME or CP in the diet. At 42 days of age, serum immunoglobulins and expression of spleen immune genes in 6 female chickens selected randomly from each dietary treatment (3×3 factorial arrangement) group were measured by enzyme-linked immunosorbent assay (ELISA) and transcriptomic analysis using RNA sequencing, respectively. The results showed that the IgM level in the low ME group chickens was significantly (p < 0.05) lower than that in other groups. In addition, immune-related genes, such as MX1, USP18, TLR4, IFNG and IL18 were significantly upregulated when the dietary nutrient density was reduced, which may put the body in an inflammatory state. This study provided general information on the molecular mechanism of the spleen immune response to variable nutrient density.

Transcriptomic and targeted immune transcript analyses confirm localized skin immune responses in Atlantic salmon towards the salmon louse.

Atlantic salmon (Salmo salar) are highly susceptible to infestations with the ectoparasite Lepeophtheirus salmonis, the salmon louse. Infestations elicit an immune response in the fish, but the response does not lead to parasite clearance, nor does it protect against subsequent infestations. It is, however, not known why the immune response is not adequate, possibly because the local response directly underneath the louse has been poorly evaluated. The present study describes the transcriptomic response by RNA sequencing of skin at the site of copepodid attachment. Analysing differentially expressed genes, 2864 were higher and 1357 were lower expressed at the louse attachment site compared to uninfested sites in the louse infested fish, while gene expression at uninfested sites were similar to uninfested control fish. The transcriptional patterns of selected immune genes were further detailed in three skin compartments/types: Whole skin, scales only and fin tissue. The elevation of pro-inflammatory cytokines and immune cell marker transcripts observed in whole skin and scale samples were not induced in fin, and a higher cytokine transcript level in scale samples suggest it can be used as a nonlethal sampling method to enhance selective breeding trials. Furthermore, the immune response was followed in both skin and anterior kidney as the infestation developed. Here, newly moulted preadult 1 stage lice induced a higher immune response than chalimi and adult lice. Overall, infestation with salmon louse induce a modest but early immune response with an elevation of mainly innate immune transcripts, with the response primarily localized to the site of attachment.

Gene Signatures and Associated Transcription Factors of Allergic Rhinitis: KLF4 Expression Is Associated with Immune Response.

This study is aimed at investigating the potential molecular features of allergic rhinitis (AR) and identifying gene signatures and related transcription factors using transcriptome analysis and in silico datasets. Transcriptome profiles were obtained using three independent cohorts (GSE101720, GSE19190, and GSE46171) comprising healthy controls (HC) and patients with AR. The pooled dataset (n = 82) was used to identify the critical signatures of AR compared with HC. Subsequently, key transcription factors were identified by a combined analysis using transcriptome and in silico datasets. Gene ontology: bioprocess (GO: BP) analysis using differentially expressed genes (DEGs) revealed that immune response-related genes were significantly enriched in AR compared with HC. Among them, IL1RL1, CD274, and CD44 were significantly higher in AR patients. We also identified key transcription factors between HC and AR using the in silico dataset and found that AR samples frequently express KLF transcription factor 4 (KLF4), which regulates immune response-related genes including IL1RL1, CD274, and CD44 in human nasal epithelial cells. Our integrative analysis of transcriptomic regulation provides new insights into AR, which may help in developing precision management for patients with AR.

Genes associated with cellular senescence favor melanoma prognosis by stimulating immune responses in tumor microenvironment.

PURPOSE: Skin cutaneous melanoma (SKCM), a malignant tumor from melanocytes, is the fifth most prevalent tumor. Immune checkpoint inhibitor (ICI) immunotherapy improves prognosis of SKCM, but immune response varies for different populations. Cellular senescence in the tumor microenvironment (TME) promotes antitumor immunity, mediated by dendritic cells (DC) and CD8+ T cells. Therefore, we sought to explore the role of cellular senescence in the TME of SKCM through bioinformatics and machine learning. METHODS: First, we obtained 93 cellular senescence-prognosis genes (CSPGs) by univariate survival analysis. Thereafter, 23 optimal CSPGs were obtained by least absolute shrinkage and selection operator (lasso) analysis. Based on the riskscore obtained by lasso analysis and clinical information from multivariate cox, we obtained the nomogram of SKCM, which was validated in the validation cohort. Based on the riskscore, the patients were split into low- and high-risk groups. Functional differences between the two groups were analyzed using Metascape and GSEA, and immune infiltration differences were achieved by multiple algorithms. We obtained a risk prediction nomogram for the validated SKCM based on the lasso model by univariate and multivariate cox regression analysis. RESULTS: In the low-risk group, immune responses were in an active state. NK, CD8+ T, DC, macrophages, and neutrophils were significantly upregulated, and ICI-relevant genes were notably upregulated. With the differentially expressed genes (DEGs) and optimal CSPGs, we obtained the hub genes: NOX4, NTN4, PROX1, and TRPM8. The hub genes were mainly expressed by cancer-associated fibroblasts (CAFs) and endothelial cells by single cell analysis, which were mainly associated with angiogenesis. CONCLUSION: Genes associated with cellular senescence favor SKCM prognosis by stimulating immune responses in TME. Patients with high expression of cellular senescence associated genes in the TME might have better benefit from ICI immunotherapy. Cellular senescence functions as a pro-tumor agent in mesenchymal cells and needs further study.