A Pan-Cancer Analysis of Clinical Prognosis and Immune Infiltration of CKS1B in Human Tumors

Integrated bulk and single-cell RNA sequencing identifies an aneuploidy-based gene signature to predict sensitivity of lung adenocarcinoma to traditional chemotherapy drugs and patients' prognosis

Background

Patients with lung adenocarcinoma (LUAD) often develop a poor prognosis. Currently, researches on prognostic and immunotherapeutic capacity of aneuploidy-related genes in LUAD are limited.

Methods

Genes related to aneuploidy were screened based on bulk RNA sequencing data from public databases using Spearman method. Next, univariate Cox and Lasso regression analyses were performed to establish an aneuploidy-related riskscore (ARS) model. Results derived from bioinformatics analysis were further validated using cellular experiments. In addition, typical LUAD cells were identified by subtype clustering, followed by SCENIC and intercellular communication analyses. Finally, ESTIMATE, ssGSEA and CIBERSORT algorithms were employed to analyze the potential relationship between ARS and tumor immune environment.

Results

A five-gene ARS signature was developed. These genes were abnormally high-expressed in LUAD cell lines, and in particular the high expression of CKS1B promoted the proliferative, migratory and invasive phenotypes of LUAD cell lines. Low ARS group had longer overall survival time, higher degrees of inflammatory infiltration, and could benefit more from receiving immunotherapy. Patients in low ASR group responded more actively to traditional chemotherapy drugs (Erlotinib and Roscovitine). The scRNA-seq analysis annotated 17 cell subpopulations into seven cell clusters. Core transcription factors (TFs) such as CREB3L1 and CEBPD were enriched in high ARS cell group, while TFs such as BCLAF1 and UQCRB were enriched in low ARS cell group. CellChat analysis revealed that high ARS cell groups communicated with immune cells via SPP1 (ITGA4-ITGB1) and MK (MDK-NCl) signaling pathways.

Conclusion

In this research, integrative analysis based on the ARS model provided a potential direction for improving the diagnosis and treatment of LUAD.

Multi-omic profiling of simultaneous ductal carcinoma in situ and invasive breast cancer

PURPOSE

The progression of ductal carcinoma in situ (DCIS) to invasive breast carcinoma (IBC) in humans is highly variable. To better understand the relationship between them, we performed a multi-omic characterization of co-occurring DCIS and IBC lesions in a cohort of individuals.

METHODS

Formalin-fixed paraffin-embedded tissue samples from 50 patients with co-occurring DCIS and IBC lesions were subjected to DNA-seq and whole transcriptome RNA-seq. Paired DCIS and IBC multi-omics profiles were then interrogated for DNA mutations, gene expression profiles and pathway analysis.

RESULTS

Most small variants and copy number variations were shared between co-occurring DCIS and IBC lesions, with IBC exhibiting on average a higher degree of additional mutations. However, 36% of co-occurring lesions shared no common mutations and 49% shared no common copy number variations. The most frequent genomic variants in both DCIS and IBC were PIK3CA, TP53, KMT2C, MAP3K1, GATA3 and SF3B1, with KMT2C being more frequent in DCIS and TP53 and MAP3K1 more frequent in IBC, though the numbers are too small for definitive conclusions. The most frequent copy number variations were seen in MCL1, CKSB1 and ERBB2. ERBB2 changes were not seen in IBC unless present in the corresponding DCIS. Transcriptional profiles were highly distinct between DCIS and IBC, with DCIS exhibiting upregulation of immune-related signatures, while IBC showed significant overexpression in genes and pathways associated with cell division and proliferation. Interestingly, DCIS and IBC exhibited significant differential expression of different components of extracellular matrix (ECM) formation and regulation, with DCIS showing overexpression of ECM-membrane interaction components while IBC showed upregulation of genes associated with fibronectin and invadopodia.

CONCLUSION

While most co-occurring DCIS and IBC were mutationally similar and suggestive of a common clonal progenitor, transcriptionally the lesions are highly distinct, with IBC expressing key pathways that facilitate invasion and proliferation. These results are suggestive of additional levels of regulation, epigenetic or other, that facilitate the acquisition of invasive properties during tumor evolution.

Identification of novel potential homologous repair deficiency-associated genes in pancreatic adenocarcinoma via WGCNA coexpression network analysis and machine learning

Homologous repair deficiency (HRD) impedes double-strand break repair, which is a common driver of carcinogenesis. Positive HRD status can be used as theranostic markers of response to platinum- and PARP inhibitor-based chemotherapies. Here, we aimed to fully investigate the therapeutic and prognostic potential of HRD in pancreatic adenocarcinoma (PAAD) and identify effective biomarkers related to HRD using comprehensive bioinformatics analysis. The HRD score was defined as the unweighted sum of the LOH, TAI, and LST scores, and it was obtained based on the previous literature. To characterize PAAD immune infiltration subtypes, the "ConsensusClusterPlus" package in R was used to conduct unsupervised clustering. A WGCNA was conducted to elucidate the gene coexpression modules and hub genes in the HRD-related gene module of PAAD. The functional enrichment study was performed using Metascape. LASSO analysis was performed using the "glmnet" package in R, while the random forest algorithm was realized using the "randomForest" package in R. The prognostic variables were evaluated using univariate Cox analysis. The prognostic risk model was built using the LASSO approach. ROC curve and KM survival analyses were performed to assess the prognostic potential of the risk model. The half-maximal inhibitory concentration (IC50) of the PARP inhibitors was estimated using the "pRRophetic" package in R and the Genomics of Drug Sensitivity in Cancer database. The "rms" package in R was used to create the nomogram. A high HRD score indicated a poor prognosis and an advanced clinical process in PAAD patients. PAAD tumors with high HRD levels revealed significant T helper lymphocyte depletion, upregulated levels of cancer stem cells, and increased sensitivity to rucaparib, Olaparib, and veliparib. Using WGCNA, 11 coexpression modules were obtained. The red module and 122 hub genes were identified as the most correlated with HRD in PAAD. Functional enrichment analysis revealed that the 122 hub genes were mainly concentrated in cell cycle pathways. One novel HRD-related gene signature consisting of CKS1B, HJURP, and TPX2 were screened via LASSO analysis and a random forest algorithm, and they were validated using independent validation sets. No direct association between HRD and CKS1B, HJURP, or TPX2 has not been reported in the literature so far. Thus, these findings indicated that CKS1B, HJURP, and TPX2 have potential as diagnostic and prognostic biomarkers for PAAD. We constructed a novel HRD-related prognostic model that provides new insights into PAAD prognosis and immunotherapy. Based on bioinformatics analysis, we comprehensively explored the therapeutic and prognostic potential of HRD in PAAD. One novel HRD-related gene signature consisting of CKS1B, HJURP, and TPX2 were identified through the combination of WGCNA, LASSO analysis and a random forest algorithm. A novel HRD-related risk model that can predict clinical prognosis and immunotherapeutic response in PAAD patients was constructed.