TNFSF9 Is a Prognostic Biomarker and Correlated with Immune Infiltrates in Pancreatic Cancer

TNFSF9 Is Associated with Favorable Tumor Immune Microenvironment in Patients with Renal Cell Carcinoma Who Are Treated with the Combination Therapy of Nivolumab and Ipilimumab

Combination therapy of nivolumab and ipilimumab (NIVO + IPI) for metastatic renal cell carcinoma (mRCC) has shown efficacy, but approximately 20% of patients experience disease progression in the early stages of treatment. No useful biomarkers have been reported to date. Therefore, it is desirable to identify biomarkers to predict treatment responses in advance. We examined the tumor microenvironment (TME)-related gene expression in mRCC patients treated with NIVO + IPI, between the response and non-response groups, using tumor tissues, before administering NIVO + IPI. In TME-related genes, TNFSF9 expression was identified as a candidate for the predictive biomarker. Its expression discriminated between the response and non-response groups with 88.89% sensitivity and 87.50% specificity (AUC = 0.9444). We further analyzed the roles of TNFSF9 in TME using bioinformatics from The Cancer Genome Atlas (TCGA) cohort. An adaptive immune response was activated in the TNFSF9-high-expression tumors. Indeed, T follicular helper cells, plasma B cells, and tumor-infiltrating CD8+ T cells were increased in the tumors, which indicates the promotion of humoral immunity due to enhanced T-B interactions. However, as the number of regulatory T cells (Treg) increased in the tumors, the percentage of dysfunctional T cells also increased. This suggests that not only PD-1 but also CTLA-4 inhibition may have suppressed Treg activation and improved the therapeutic effect in the TNFSF9 high-expression tumors. Therefore, TNFSF9 may predict the therapeutic efficacy of NIVO + IPI for mRCC and allow more appropriate patient selection.

Ferroptosis mediated by TNFSF9 interferes in acute ischaemic stroke reperfusion injury with the progression of acute ischaemic stroke

In this study, we investigated the potential involvement of TNFSF9 in reperfusion injury associated with ferroptosis in acute ischaemic stroke patients, mouse models and BV2 microglia. We first examined TNFSF9 changes in peripheral blood from stroke patients with successful reperfusion, and constructed oxygen-glucose deprivation-reperfusion (OGD-R) on BV2 microglia, oxygen-glucose deprivation for 6 h followed by reoxygenation and re-glucose for 24 h, and appropriate over-expression or knockdown of TNFSF9 manipulation on BV2 cells and found that in the case of BV2 cells encountering OGD-R over-expression of TNFSF9 resulted in increased BV2 apoptosis. Still, the knockdown of TNFSF9 ameliorated apoptosis and ferroptosis. In an in vivo experiment, we constructed TNFSF9 over-expression or knockout mice by intracerebral injection of TNFSF9-OE or sh-TNFSF9 adenovirus. We performed the middle cerebral artery occlusion (MCAO) model on day four, 24 h after ligation of the proximal artery, for half an hour to recanalize. As luck would have it, over-expression of TNFSF9 resulted in increased brain infarct volumes, neurological function scores and abnormalities in TNFSF9-related TRAF1 and ferroptosis-related pathways, but knockdown of TNFSF9 improved brain infarcts in mice as well as reversing TNFSF9-related signalling pathways. In conclusion, our data provide the first evidence that TNFSF9 triggers microglia activation by activating the ferroptosis signalling pathway following ischaemic stroke, leading to brain injury and neurological deficits.

A 15-Inflammation-Related Gene Signature Predicts the Prognosis of Patients With Pancreatic Ductal Adenocarcinoma

Chronic inflammation promotes the development of pancreatic ductal adenocarcinoma (PDAC) and PDAC-related inflammatory tumor microenvironment facilitates tumor growth and metastasis. Thus, we aimed to study the association between inflammatory response and prognosis in patients with PDAC. We conducted the whole transcriptomic sequencing using tissue samples collected from patients diagnosed with PDAC (n = 106) recruited from Shandong Cancer Hospital. We first constructed a prognostic signature using 15 inflammation-related genes in The Cancer Genome Atlas (TCGA) cohort (n = 177) and further validated it in an independent International Cancer Genome Consortium (ICGC) cohort (n = 90) and our in-house cohort. PDAC patients with a higher risk score had poorer overall survival (OS) (P < 0.001; HR, 3.02; 95% CI, 1.94-4.70). The association between the prognostic signature and OS remained significant in the multivariable Cox regression adjusting for age, sex, alcohol exposure, diabetes, and stage (P < 0.001; HR, 2.91; 95% CI, 1.73-4.89). This gene signature also robustly predicted prognosis in the ICGC cohort (P = 0.01; HR, 1.94; 95% CI, 1.14-3.30) and our cohort (P < 0.001; HR, 2.40; 95% CI, 1.45-3.97). Immune subtype C3 (inflammatory) was enriched and CD8+ T cells were higher in patients with a lower risk score (P < 0.05). Furthermore, PDAC patients with higher risk scores were more sensitive to chemotherapy, immunotherapy, and PARP inhibitors (P < 0.05). In sum, we identified a novel gene signature that was associated with inflammatory response for risk stratification, prognosis prediction, and therapy guidance in PDAC patients. Future studies are warranted to validate the clinical utility of the signature.

The potential use and experimental validation of genomic instability-related lncRNA in pancreatic carcinoma

This study explored the potential role of long noncoding RNA (lncRNAs) associated with genomic instability in the diagnosis and treatment of pancreatic adenocarcinoma (PAAD). Transcriptome and single-nucleotide variation data of PAAD samples were downloaded from the cancer genome atlas database to explore genomic instability-associated lncRNAs. We constructed a genomic instability-associated lncRNA prognostic signature. Then gene ontology and Kyoto encyclopedia of genes and genomes enrichment analyses were used to explore the physiological role of lncRNAs involved in genomic instability. Tumor microenvironments, immunotherapy response, immune cell infiltration, immune checkpoint, and drug sensitivity were compared between high-risk and low-risk groups. In vitro experiments were performed for external validation. Six lncRNAs associated with genomic instability were identified, capable of predicting the prognosis of PAAD. Patients were assigned to low-risk or high-risk groups using these biomarkers, with better or worse prognosis, respectively. The tumor immune score, immune cell infiltration, and efficacy of immunotherapy were worse in the high-risk group. A drug sensitivity analysis revealed the high- and low-risk groups had different half-maximal inhibitory concentrations. The expression of cancer susceptibility candidate 8 was significantly higher in tumor tissues than in normal tissues, while the expression of LYPLAL1-AS1 exhibited an opposite pattern. They may be potential diagnostic or prognostic biomarkers for patients with pancreatic cancer. Genomic instability-associated lncRNAs were explored in this study and predicted the prognosis of PAAD and stratified patients risk in PAAD. These lncRNAs also predicted the efficacy of immunotherapy and potential therapeutic targets in PAAD.

Spatiotemporally deciphering the mysterious mechanism of persistent HPV-induced malignant transition and immune remodelling from HPV-infected normal cervix, precancer to cervical cancer: Integrating single-cell RNA-sequencing and spatial transcriptome

BACKGROUND

The mechanism underlying cervical carcinogenesis that is mediated by persistent human papillomavirus (HPV) infection remains elusive.

AIMS

Here, for the first time, we deciphered both the temporal transition and spatial distribution of cellular subsets during disease progression from normal cervix tissues to precursor lesions to cervical cancer.

MATERIALS & METHODS

We generated scRNA-seq profiles and spatial transcriptomics data from nine patient samples, including two HPV-negative normal, two HPV-positive normal, two HPV-positive HSIL and three HPV-positive cancer samples.

RESULTS

We not only identified three 'HPV-related epithelial clusters' that are unique to normal, high-grade squamous intraepithelial lesions (HSIL) and cervical cancer tissues but also discovered node genes that potentially regulate disease progression. Moreover, we observed the gradual transition of multiple immune cells that exhibited positive immune responses, followed by dysregulation and exhaustion, and ultimately established an immune-suppressive microenvironment during the malignant program. In addition, analysis of cellular interactions further verified that a 'homeostasis-balance-malignancy' change occurred within the cervical microenvironment during disease progression.

DISCUSSION

We for the first time presented a spatiotemporal atlas that systematically described the cellular heterogeneity and spatial map along the four developmental steps of HPV-related cervical oncogenesis, including normal, HPV-positive normal, HSIL and cancer. We identified three unique HPV-related clusters, discovered critical node genes that determined the cell fate and uncovered the immune remodeling during disease escalation.

CONCLUSION

Together, these findings provided novel possibilities for accurate diagnosis, precise treatment and prognosis evaluation of patients with precancer and cervical cancer.

Prognostic factor identification by screening changes in differentially expressed genes in oral squamous cell carcinoma

OBJECTIVE

This study was designed to identify changes in the expression of proteins occurring during the progression of oral squamous cell carcinoma (OSCC) and to validate their impact on patient prognosis.

MATERIALS AND METHODS

The human OSCC cell line UPCI-SCC-040 was treated in vitro with TGF-β1, and transcriptome analysis of differentially expressed genes (DEGs) revealed putative candidates relative to untreated cells. The respective protein expression levels of the most important genes were immunohistochemically validated on a tissue microarray (TMA) containing tissue samples from 39 patients with OSCC and were correlated with disease-free survival (DFS) as the primary clinical endpoint.

RESULTS

Our univariate Cox proportional hazard regression (CR) analysis revealed significant correlations among positive N stage (local lymph node metastasis, p = .04), stearoyl-CoA desaturase-1 (p < .01), sclerostin (p = .01), and CD137L expression (p = .04) and DFS. Stearoyl-CoA desaturase-1 and sclerostin remained the main prognostic factors (p < .01) in the multiple CR model.

CONCLUSION

We identified changes in differentially expressed genes during OSCC progression in vitro and translated the impact of the most deregulated genes on patient prognosis. Stearoyl-CoA desaturase-1 and sclerostin acted as independent prognostic factors in OSCC and could also be interesting candidates for new cancer targeted therapeutic approaches.

A Novel L-Phenylalanine Dipeptide Inhibits the Growth and Metastasis of Prostate Cancer Cells via Targeting DUSP1 and TNFSF9

Prostate cancer (PCa) is a common malignant cancer of the urinary system. Drug therapy, chemotherapy, and radical prostatectomy are the primary treatment methods, but drug resistance and postoperative recurrence often occur. Therefore, seeking novel anti-tumor compounds with high efficiency and low toxicity from natural products can produce a new tumor treatment method. Matijin-Su [N-(N-benzoyl-L-phenylalanyl)-O-acetyl-L-phenylalanol, MTS] is a phenylalanine dipeptide monomer compound that is isolated from the Chinese ethnic medicine Matijin (Dichondra repens Forst.). Its derivatives exhibit various pharmacological activities, especially anti-tumor. Among them, the novel MTS derivative HXL131 has a significant inhibitory effect against prostate tumor growth and metastasis. This study is designed to investigate the effects of HXL131 on the growth and metastasis of human PCa cell lines PC3 and its molecular mechanism through in vitro experiments combined with proteomics, molecular docking, and gene silencing. The in vitro results showed that HXL131 concentration dependently inhibited PC3 cell proliferation, induced apoptosis, arrested cell cycle at the G2/M phase, and inhibited cell migration capacity. A proteomic analysis and a Western blot showed that HXL131 up-regulated the expression of proliferation, apoptosis, cell cycle, and migration-related proteins CYR61, TIMP1, SOD2, IL6, SERPINE2, DUSP1, TNFSF9, OSMR, TNFRSF10D, and TNFRSF12A. Molecular docking, a cellular thermal shift assay (CETSA), and gene silencing showed that HXL131 had a strong binding affinity with DUSP1 and TNFSF9, which are important target genes for inhibiting the growth and metastasis of PC3 cells. This study demonstrates that HXL131 exhibited excellent anti-prostate cancer activity and inhibited the growth and metastasis of prostate cancer cells by regulating the expression of DUSP1 and TNFSF9.

Identification of necroptosis-related subtypes and prognosis model in triple negative breast cancer

Background

Necroptosis is considered to be a new form of programmed necrotic cell death, which is associated with metastasis, progression and prognosis of various types of tumors. However, the potential role of necroptosis-related genes (NRGs) in the triple negative breast cancer (TNBC) is unclear.

Methods

We extracted the gene expression and relevant clinicopathological data of TNBC from The Cancer Genome Atlas (TCGA) databases and the Gene Expression Omnibus (GEO) databases. We analyzed the expression, somatic mutation, and copy number variation (CNV) of 67 NRGs in TNBC, and then observed their interaction, biological functions, and prognosis value. By performing Lasso and COX regression analysis, a NRGs-related risk model for predicting overall survival (OS) was constructed and its predictive capabilities were verified. Finally, the relationship between risk_score and immune cell infiltration, tumor microenvironment (TME), immune checkpoint, and tumor mutation burden (TMB), cancer stem cell (CSC) index, and drug sensitivity were analyzed.

Results

A total 67 NRGs were identified in our analysis. A small number of genes (23.81%) detected somatic mutation, most genes appeared to have a high frequency of CNV, and there was a close interaction between them. These genes were remarkably enriched in immune-related process. A seven-gene risk_score was generated, containing TPSG1, KRT6A, GPR19, EIF4EBP1, TLE1, SLC4A7, ESPN. The low-risk group has a better OS, higher immune score, TMB and CSC index, and lower IC50 value of common therapeutic agents in TNBC. To improve clinical practicability, we added age, stage_T and stage_N to the risk_score and construct a more comprehensive nomogram for predicting OS. It was verified that nomogram had good predictive capability, the AUC values for 1-, 3-, and 5-year OS were 0.847, 0.908, and 0.942.

Conclusion

Our research identified the significant impact of NRGs on immunity and prognosis in TNBC. These findings were expected to provide a new strategy for personalize the treatment of TNBC and improve its clinical benefit.

Long noncoding RNAs (lncRNAs) in pancreatic cancer progression

Long noncoding RNAs (lncRNAs) are RNA molecules involved in gene regulation at transcriptional, post-transcriptional, and epigenetic levels. LncRNAs participate in regulating apoptosis and autophagy in pancreatic cancer (PCa) and can promote and/or decrease the proliferation rate of tumor cells. The metastasis of PCa cells is tightly regulated by lncRNAs and they can affect the mechanism of epithelial-mesenchymal transition (EMT) to modulate metastasis. The drug resistance of PCa cells, especially to gemcitabine, can be affected by lncRNAs. In addition, lncRNAs enriched in exosomes can be transferred among tumor cells to regulate their proliferation and metastasis. Antitumor compounds, such as curcumin and ginsenosides, can regulate lncRNA expression in PCa therapy. As we discuss here, the expression level of lncRNAs can be considered as both a diagnostic and prognostic tool in patients with PCa.

A Novel Immunogenomic Signature to Predict Prognosis and Reveal Immune Infiltration Characteristics in Pancreatic Ductal Adenocarcinoma

Background

The immune response in the tumor microenvironment (TME) plays a crucial role in cancer progression and recurrence. We aimed to develop an immune-related gene (IRG) signature to improve prognostic predictive power and reveal the immune infiltration characteristics of pancreatic ductal adenocarcinoma (PDAC).

Methods

The Cancer Genome Atlas (TCGA) PDAC was used to construct a prognostic model as a training cohort. The International Cancer Genome Consortium (ICGC) and the Gene Expression Omnibus (GEO) databases were set as validation datasets. Prognostic genes were screened by using univariate Cox regression. Then, a novel optimal prognostic model was developed by using least absolute shrinkage and selection operator (LASSO) Cox regression. Cell type identification by estimating the relative subsets of RNA transcripts (CIBERSORT) and estimation of stromal and immune cells in malignant tumors using expression data (ESTIMATE) algorithms were used to characterize tumor immune infiltrating patterns. The tumor immune dysfunction and exclusion (TIDE) algorithm was used to predict immunotherapy responsiveness.

Results

A prognostic signature based on five IRGs (MET, ERAP2, IL20RB, EREG, and SHC2) was constructed in TCGA-PDAC and comprehensively validated in ICGC and GEO cohorts. Multivariate Cox regression analysis demonstrated that this signature had an independent prognostic value. The area under the curve (AUC) values of the receiver operating characteristic (ROC) curve at 1, 3, and 5 years of survival were 0.724, 0.702, and 0.776, respectively. We further demonstrated that our signature has better prognostic performance than recently published ones and is superior to traditional clinical factors such as grade and tumor node metastasis classification (TNM) stage in predicting survival. Moreover, we found higher abundance of CD8+ T cells and lower M2-like macrophages in the low-risk group of TCGA-PDAC, and predicted a higher proportion of immunotherapeutic responders in the low-risk group.

Conclusions

We constructed an optimal prognostic model which had independent prognostic value and was comprehensively validated in external PDAC databases. Additionally, this five-genes signature could predict immune infiltration characteristics. Moreover, the signature helped stratify PDAC patients who might be more responsive to immunotherapy.

Time Series Ovarian Transcriptome Analyses of the Porcine Estrous Cycle Reveals Gene Expression Changes during Steroid Metabolism and Corpus Luteum Development

Simple Summary

The estrous cycle, which is divided into follicular and luteal phases based on ovulation, is influenced by reproductive hormones which affect reproduction and cause changes in the reproductive system of the pig. As the main reproductive organ, the ovary is involved in ovulation and changes in the corpus luteum. We aimed to identify dynamic changes in gene expression through differentially expressed gene profiling and to provide a comprehensive understanding of the molecular mechanisms that occur in the pig ovary during the estrous cycle. The transcriptome analysis revealed that the dynamic change in gene expression was more activated in the luteal phase than in the follicular phase. Functional analysis revealed that the metestrus and diestrus periods are important in preparation for pregnancy or the next estrous cycle after ovulation.

Abstract

The porcine estrous cycle is influenced by reproductive hormones, which affect porcine reproduction and result in physiological changes in the reproductive organs. The ovary is involved in ovulation, luteinization, corpus luteum development, and luteolysis. Here, we aimed to provide a comprehensive understanding of the gene expression patterns in porcine ovarian transcriptomes during the estrous cycle through differentially expressed genes profiling and description of molecular mechanisms. The transcriptomes of porcine ovary were obtained during the estrous cycle at three-day intervals from day 0 to day 18 using RNA-seq. At seven time points of the estrous cycle, 4414 DEG were identified; these were classified into three clusters according to their expression patterns. During the late metestrus and diestrus periods, the expression in cluster 1 increased rapidly, and steroid biosynthesis was significant in the pathway. Cluster 2 gene expression patterns represented the cytokine–cytokine receptor interaction in significant pathways. In cluster 3, the hedgehog signaling pathway was selected as the significant pathway. Our study exhibited dynamic gene expression changes with these three different patterns of cluster 1, 2, and 3. The results helped identify the functions and related significant genes especially during the late metestrus and diestrus periods in the estrous cycle.

TNFSF9 promotes metastasis of pancreatic cancer through Wnt/Snail signaling and M2 polarization of macrophages

Early metastasis of pancreatic cancer (PC) leads to high mortality, and the underlying mechanism of metastasis remains unclear. Tumor necrosis factor superfamily member 9 (TNFSF9) is associated with poor prognosis in PC. Here, we investigated the effect of TNFSF9 on PC proliferation and apoptosis, and focused on the effect of TNFSF9 on PC metastasis and its potential mechanism. We found that TNFSF9 promotes PC metastasis in vivo and in vitro, and may be partially dependent on the Wnt/Snail signaling pathway. In addition, TNFSF9 also regulates the release of cytokines IL-10 and transforming growth factor-β (TGF-β) in pancreatic cancer cells through Wnt signaling to induce the M2 polarization of macrophages and promote the migration of PC cells. Overall, our study found that TNFSF9 may directly promote PC metastasis or indirectly promote PC metastasis through macrophage M2 polarization. Our study provides a new costimulatory target for the treatment of PC.

An Integrated Data Analysis of mRNA, miRNA and Signaling Pathways in Pancreatic Cancer

Although many genes and miRNAs have been reported for various cancers, pancreatic cancer's specific genes or miRNAs have not been studied precisely yet. Therefore, we have analyzed the gene and miRNA expression profile of pancreatic cancer data in the gene expression omnibus (GEO) database. The microarray-derived miRNAs and mRNAs were annotated by gene ontology (GO) and signaling pathway analysis. We also recognized mRNAs that were targeted by miRNA through the mirDIP database. An integrated analysis of the microarray revealed that only 6 out of 43 common miRNAs had significant differences in their expression profiles between the tumor and normal groups (P value < 0.05 and |log Fold Changes (logFC)|> 1). The hsa-miR-210 had upregulation, whereas hsa-miR-375, hsa-miR-216a, hsa-miR-217, hsa-miR-216b and hsa-miR-634 had downregulation in pancreatic cancer (PC). The analysis results also revealed 109 common mRNAs by microarray and mirDIP 4.1 databases. Pathway analysis showed that amoebiasis, axon guidance, PI3K-Akt signaling pathway, absorption and focal adhesion, adherens junction, platelet activation, protein digestion, human papillomavirus infection, extracellular matrix (ECM) receptor interaction, and riboflavin metabolism played important roles in pancreatic cancer. GO analysis revealed the significant enrichment in the three terms of biological process, cellular component, and molecular function, which were identified as the most important processes associated strongly with pancreatic cancer. In conclusion, DTL, CDH11, COL5A1, ITGA2, KIF14, SMC4, VCAN, hsa-mir-210, hsa-mir-217, hsa-mir-216a, hsa-mir-216b, hsa-mir-375 and hsa-mir-634 can be reported as the novel diagnostic or even therapeutic markers for the future studies. Also, the hsa-mir-107 and hsa-mir-125a-5p with COL5A1, CDH11 and TGFBR1 genes can be introduced as major miRNA and genes on the miRNA-drug-mRNA network. The new regulatory network created in our study could give a deeper knowledge of the pancreatic cancer.

TANK-Binding Kinase 1 (TBK1) Serves as a Potential Target for Hepatocellular Carcinoma by Enhancing Tumor Immune Infiltration

Background

Numerous cancer types present the aberrant TANK-binding kinase 1 (TBK1) expression, which plays an important role in driving inflammation and innate immunity. However, the prognostic role of TBK1 and its relationship with immune cell infiltration in hepatocellular carcinoma (HCC) remain unclear.

Methods

The expression and prognostic value of TBK1 was analyzed by Tumor Immune Estimation Resource (TIMER), Kaplan-Meier plotter and Gene Expression Profiling Interactive Analysis (GEPIA), Clinical Proteomic Tumor Analysis Consortium (CPTAC) and further confirmed in the present cohort of patients with HCC. The association between TBK1 and HCC immune infiltrates, and its potential mechanism were investigated via analyses of the Tumor Immune Estimation Resource, tumor-immune system interactions database (TISIDB), CIBERSORT, STRING, and Metascape. The effect of TBK1 on immune infiltrates and the therapeutic value of targeting TBK1 were further investigated in a HCC mouse model by treatment with a TBK1 antagonist.

Results

The level of TBK1 expression in HCC was higher than that measured in normal tissues, and associated with poorer overall survival (GEPIA: hazard ratio [HR]=1.80, P=0.038; Kaplan-Meier plotter: HR=1.87, P<0.001; CPTAC: HR=2.23, P=0.007; Our cohort: HR=2.92, P=0.002). In addition, high TBK1 expression was found in HCC with advanced TNM stage and identified as an independent poor prognostic factor for overall survival among patients with HCC. In terms of immune infiltration, tumor tissues from HCC patients with high TBK1 expression had a low proportion of CD8+ T cells, and TBK1 expression did not show prognostic value in HCC patients with enriched CD8+ T cells. Furthermore, TBK1 expression was positively correlated with the markers of T cell exhaustion and immunosuppressive cells in the HCC microenvironment. Mechanistically, the promotion of HCC immunosuppression by TBK1 was involved in the regulation of inflammatory cytokines. In vivo experiments revealed that treatment with a TBK1 antagonist delayed HCC growth by increasing the number of tumor-infiltrating CD8+ T cells.

Conclusions

The up-regulated expression of TBK1 may be useful in predicting poor prognosis of patients with HCC. In addition, TBK1, which promotes the HCC immunosuppressive microenvironment, may be a potential immunotherapeutic target for patients with HCC.

Identification of ribosomal protein family in triple-negative breast cancer by bioinformatics analysis

Triple-negative breast cancer (TNBC) accounts for ∼20% of all breast cancer (BC) cases. The management of TNBC represents a challenge due to its worse prognosis, heterogeneity and lack of targeted therapy. Moreover, its mechanisms are not fully clear. The aim of the study is to identify crucial genes between TNBC and non-TNBC for underlying targets for diagnostic and therapeutic methods of TNBC. The differentially expressed genes (DEGs) between TNBC and non-TNBC were selected from the Gene Expression Omnibus (GEO) database after the integrated analysis of two datasets (GSE65194 and GSE76124). Then Gene ontology (GO) and KEGG analysis were performed by DAVID database, protein-protein interaction (PPI) of DEGs was constructed by Search Tool for the Retrieval of Reciprocity Genes (STRING) database. Furthermore, centrality analysis and module analysis were carried out by Cytoscape to analyze the TNBC-related PPI. Subsequently, overall survival (OS) analysis was performed by GEPIA. Finally, the expressions of these key genes in TNBC and non-TNBC tissues were tested by qRT-PCR. The results showed that 955 DEGs were obtained, which were mainly enriched in ribosome, ribosomal subunit, and so on. Moreover, 19 candidate genes were focused on by centrality analysis and module analysis. Furthermore, we found the low expressions of ribosomal protein S9 (RPS9), ribosomal protein S14 (RPS14), ribosomal protein S27 (RPS27), ribosomal protein L11 (RPL11) and ribosomal protein L14 (RPL14) were related to a poor OS in BC patients. Additionally, qRT-PCR results suggested that these five genes were notably down-regulated in TNBC tissues. In summary, the present study suggests that ribosomal proteins are related to TNBC, and they may play an important role in the diagnosis, treatment and prognosis of TNBC.

FunRes: resolving tissue-specific functional cell states based on a cell-cell communication network model

The functional specialization of cell types arises during development and is shaped by cell-cell communication networks determining a distribution of functional cell states that are collectively important for tissue functioning. However, the identification of these tissue-specific functional cell states remains challenging. Although a plethora of computational approaches have been successful in detecting cell types and subtypes, they fail in resolving tissue-specific functional cell states. To address this issue, we present FunRes, a computational method designed for the identification of functional cell states. FunRes relies on scRNA-seq data of a tissue to initially reconstruct the functional cell-cell communication network, which is leveraged for partitioning each cell type into functional cell states. We applied FunRes to 177 cell types in 10 different tissues and demonstrated that the detected states correspond to known functional cell states of various cell types, which cannot be recapitulated by existing computational tools. Finally, we characterize emerging and vanishing functional cell states in aging and disease, and demonstrate their involvement in key tissue functions. Thus, we believe that FunRes will be of great utility in the characterization of the functional landscape of cell types and the identification of dysfunctional cell states in aging and disease.