Construction of Bone Metastasis-Specific Regulation Network Based on Prognostic Stemness-Related Signatures in Prostate Cancer.

Background: We planned to uncover the cancer stemness-related genes (SRGs) in prostate cancer (PCa) and its underlying mechanism in PCa metastasis. Methods: We acquired the RNA-seq data of 406 patients with PCa from the TCGA database. Based on the mRNA stemness index (mRNAsi) calculated by one-class logistic regression (OCLR) algorithm, SRGs in PCa were extracted by WGCNA. Univariate and multivariate regression analyses were applied to uncover OS-associated SRGs. Gene Set Variation Analysis (GSVA), Gene Set Enrichment Analysis (GSEA), and Pearson's correlation analysis were performed to discover the possible mechanism of PCa metastasis. The significantly correlated transcription factors of OS-associated SRGs were also identified by Pearson's correlation analysis. ChIP-seq was applied to validate the binding relationship of TFs and OS-associated SRGs and spatial transcriptome and single-cell sequencing were performed to uncover the location of key biomarkers expression. Lastly, we explored the specific inhibitors for SRGs using CMap algorithm. Results: We identified 538 differentially expressed genes (DEGs) between non-metastatic and metastatic PCa. Furthermore, OS-associated SRGs were identified. The Pearson correlation analysis revealed that FOXM1 was significantly correlated with NEIL3 (correlation efficient =0.89, p < 0.001) and identified hallmark_E2F_targets as the potential pathway mechanism of NEIL3 promoting PCa metastasis (correlation efficient =0.58, p < 0.001). Single-cell sequencing results indicated that FOXM1 regulating NEIL3 may get involved in the antiandrogen resistance of PCa. Rottlerin was discovered to be a potential target drug for PCa. Conclusion: We constructed a regulatory network based on SRGs associated with PCa metastasis and explored possible mechanism.

Identification of prognostic and bone metastasis​-related alternative splicing signatures in mesothelioma.

Mesothelioma (MESO) is an infrequent tumor derived from mesothelial cells of pleura, peritoneum, pericardium, and tunica vaginalis testis. Despite advancement in technologies and better understanding of tumor progression mechanism, the prognosis of MESO remains poor. The role of alternative splicing events (ASEs) in the oncogenesis, tumor metastasis and drug resistance has been widely discussed in multiple cancers. But the prognosis and potential therapeutic value of ASEs in MESO were not clearly studied by now. We constructed a prognostic model using RNA sequencing data and matched ASE data of MESO patients obtained from the TCGA and TCGASpliceSeq database. A total of 3,993 ASEs were identified associated with overall survival using Cox regression analysis. Eight of them were finally figured out to institute the model by lasso regression analysis. The risk score of the model can predict the prognosis independently. Among the identified 390 splicing factors (SF), HSPA1A and DDX3Y was significantly associated with 43 OS-SEs. Among these OS-SEs, SNX5-58744-AT (p = 0.048) and SNX5-58745-AT (p = 0.048) were significantly associated with bone metastasis. Co-expression analysis of signal pathways and SNX5-58744-AT, SNX5-58745-AT was also depicted using GSVA. Finally, we proposed that splicing factor (SF) HSPA1A could regulate SNX5-58744-AT (R = -0.414) and SNX5-58745-AT (R = 0.414) through the pathway "Class I MHC mediated antigen processing and presentation" (R = 0.400). In this way, tumorigenesis and bone metastasis of MESO were controlled.

Targeting Lymphotoxin Beta and Paired Box 5: a potential therapeutic strategy for soft tissue sarcoma metastasis.

BACKGROUND: Soft tissue sarcomas (STS) has a high rate of early metastasis. In this study, we aimed to uncover the potential metastasis mechanisms and related signaling pathways in STS with differentially expressed genes and tumor-infiltrating cells. METHODS: RNA-sequencing (RNA-seq) of 261 STS samples downloaded from the Cancer Genome Atlas (TCGA) database were used to identify metastasis-related differentially expressed immune genes and transcription factors (TFs), whose relationship was constructed by Pearson correlation analysis. Metastasis-related prediction model was established based on the most significant immune genes. CIBERSORT algorithm was performed to identify significant immune cells co-expressed with key immune genes. The GSVA and GSEA were performed to identify prognosis-related KEGG pathways. Ultimately, we used the Pearson correlation analysis to explore the relationship among immune genes, immune cells, and KEGG pathways. Additionally, key genes and regulatory mechanisms were validated by single-cell RNA sequencing and ChIP sequencing data. RESULTS: A total of 204 immune genes and 12 TFs, were identified. The prediction model achieved a satisfactory effectiveness in distant metastasis with the Area Under Curve (AUC) of 0.808. LTB was significantly correlated with PAX5 (P < 0.001, R = 0.829) and hematopoietic cell lineage pathway (P < 0.001, R = 0.375). The transcriptional regulatory pattern between PAX5 and LTB was validated by ChIP sequencing data. CONCLUSIONS: We hypothesized that down-regulated LTB (immune gene) modulated by PAX5 (TF) in STSs may have the capability of inducing cancer cell metastasis in patients with STS.

Identification of prognostic and metastasis-related alternative splicing signatures in hepatocellular carcinoma.

As the most common neoplasm in digestive system, hepatocellular carcinoma (HCC) is one of the most important leading cause of cancer deaths worldwide. Its high-frequency metastasis and relapse rate lead to the poor survival of HCC patients. However, the mechanism of HCC metastasis is still unclear. Alternative splicing events (ASEs) have a great effect in cancer development, progression and metastasis. We downloaded RNA sequencing and seven types of ASEs data of HCC samples, in order to explore the mechanism of ASEs underlying tumorigenesis and metastasis of HCC. The data were taken from the The Cancer Genome Atlas (TCGA) and TCGASpliceSeq databases. Univariate Cox regression analysis was used to determine a total of 3197 overall survival-related ASEs (OS-SEs). And based on five OS-SEs screened by Lasso regression, we constructed a prediction model with the Area Under Curve of 0.765. With a good reliability of the model, the risk score was also proved to be an independent predictor. Among identified 390 candidate SFs, Y-box protein 3 (YBX3) was significantly correlated with OS and metastasis. Among 177 ASEs, ATP-binding cassette subfamily A member 6 (ABCA6)-43162-AT and PLIN5-46808-AT were identified both associated with OS, bone metastasis and co-expressed with SFs. Then we identified primary bile acid biosynthesis as survival-related (KEGG) pathway by Gene Set Variation Analysis (GSVA) and univariate regression analysis, which was correlated with ABCA6-43162-AT and PLIN5-46808-AT. Finally, we proposed that ABCA6-43162-AT and PLIN5-46808-AT may contribute to HCC poor prognosis and metastasis under the regulation of aberrant YBX3 through the pathway of primary bile acid biosynthesis.