Genome Instability and Long Noncoding RNA Reveal Biomarkers for Immunotherapy and Prognosis and Novel Competing Endogenous RNA Mechanism in Colon Adenocarcinoma

Combining single-cell and bulk RNA sequencing, NK cell marker genes reveal a prognostic and immune status in pancreatic ductal adenocarcinoma

The NK cell is an important component of the tumor microenvironment of pancreatic ductal adenocarcinoma (PDAC), also plays a significant role in PDAC development. This study aimed to explore the relationship between NK cell marker genes and prognosis, immune response of PDAC patients. By scRNA-seq data, we found the proportion of NK cells were significantly downregulated in PDAC and 373 NK cell marker genes were screened out. By TCGA database, we enrolled 7 NK cell marker genes to construct the signature for predicting prognosis in PDAC patients. Cox analysis identified the signature as an independent factor for pancreatic cancer. Subsequently, the predictive power of signature was validated by 6 GEO datasets and had an excellent evaluation. Our analysis of relationship between the signature and patients' immune status revealed that the signature has a strong correlation with immunocyte infiltration, inflammatory reaction, immune checkpoint inhibitors (ICIs) response. The NK cell marker genes are closely related to the prognosis and immune capacity of PDAC patients, and they have potential value as a therapeutic target.

Construction and Validation of a Novel Prognosis Model in Colon Cancer Based on Cuproptosis-Related Long Non-Coding RNAs

Colon cancer (CC) is one of the most common (6%) malignancies and leading cause of cancer-associated death (more than 0.5 million) worldwide, which demands reliable prognostic biomarkers. Cuproptosis is a novel modality of regulated cell death triggered by the accumulation of intracellular copper. LncRNAs have been reported as prognostic signatures in different types of tumors. However, the correlation between cuproptosis-related lncRNAs (CRLs) and CC remains unclear. Data of CC patients were downloaded from public databases. The prognosis-associated CRLs were identified by co-expression analysis and univariate Cox. Least absolute shrinkage and selection operator were utilized to construct the CRLs-based prognostic signature in silico for CC patients. CRLs level was validated in human CC cell lines and patient tissues. ROC curve and Kaplan-Meier curve results revealed that high CRLs-risk score was associated with poor prognosis in CC patients. Moreover, the nomogram revealed that this model possessed a steady prognostic prediction capability with C-index as 0.68. More importantly, CC patients with high CRLs-risk score were more sensitive to eight targeted therapy drugs. The prognostic prediction power of the CRLs-risk score was further confirmed by cell lines, tissues and two independent CC cohorts. This study constructed a novel ten-CRLs-based prognosis model for CC patients. The CRLs-risk score is expected to serve as a promising prognostic biomarker and predict targeted therapy response in CC patients.

Circulating miRNA Expression Profiles and Machine Learning Models in Association with Response to Irinotecan-Based Treatment in Metastatic Colorectal Cancer

Colorectal cancer represents a leading cause of cancer-related morbidity and mortality. Despite improvements, chemotherapy remains the backbone of colorectal cancer treatment. The aim of this study is to investigate the variation of circulating microRNA expression profiles and the response to irinotecan-based treatment in metastatic colorectal cancer and to identify relevant target genes and molecular functions. Serum samples from 95 metastatic colorectal cancer patients were analyzed. The microRNA expression was tested with a NucleoSpin miRNA kit (Machnery-Nagel, Germany), and a machine learning approach was subsequently applied for microRNA profiling. The top 10 upregulated microRNAs in the non-responders group were hsa-miR-181b-5p, hsa-miR-10b-5p, hsa-let-7f-5p, hsa-miR-181a-5p, hsa-miR-181d-5p, hsa-miR-301a-3p, hsa-miR-92a-3p, hsa-miR-155-5p, hsa-miR-30c-5p, and hsa-let-7i-5p. Similarly, the top 10 downregulated microRNAs were hsa-let-7d-5p, hsa-let-7c-5p, hsa-miR-215-5p, hsa-miR-143-3p, hsa-let-7a-5p, hsa-miR-10a-5p, hsa-miR-142-5p, hsa-miR-148a-3p, hsa-miR-122-5p, and hsa-miR-17-5p. The upregulation of microRNAs in the miR-181 family and the downregulation of those in the let-7 family appear to be mostly involved with non-responsiveness to irinotecan-based treatment.

Identification of an m6A-Related Long Noncoding RNA Risk Model for Predicting Prognosis and Directing Treatments in Patients With Colon Adenocarcinoma

N6-methyladenosine (m6A) and lncRNAs have been implicated in the development of colon cancer, including tumorigenesis, migration, and invasion. However, the specific effect of m6A regulators on lncRNAs is not clear, and m6A-related lncRNAs may be new prognostic biomarkers and may help direct treatment and medication. We identified 29 prognostic m6A-related lncRNAs and constructed a risk model using 12 lncRNAs. The model was an independent prognostic factor and could accurately predict the prognosis. A stable and robust nomogram that combined the model and pathologic stage was constructed. A total of 2,424 differentially expressed genes (DEGs) were identified based on the model. Functional analysis of the DEGs showed that they were associated with tumor progression, helping investigate the underlying biological functions and signaling pathways of the risk model. In addition, the low-risk group based on the risk model had more sensitivity to afatinib, metformin, and GW.441756, and patients with low risk would more likely respond to immunotherapy. Moreover, patients with higher risk were more sensitive to olaparib, bexarotene, and doxorubicin.

Genomic instability genes in lung and colon adenocarcinoma indicate organ specificity of transcriptomic impact on Copy Number Alterations

Genomic instability (GI) in cancer facilitates cancer evolution and is an exploitable target for therapy purposes. However, specific genes involved in cancer GI remain elusive. Causal genes for GI via expressions have not been comprehensively identified in colorectal cancers (CRCs). To fill the gap in knowledge, we developed a data mining strategy (Gene Expression to Copy Number Alterations; "GE-CNA"). Here we applied the GE-CNA approach to 592 TCGA CRC datasets, and identified 500 genes whose expression levels associate with CNA. Among these, 18 were survival-critical (i.e., expression levels correlate with significant differences in patients' survival). Comparison with previous results indicated striking differences between lung adenocarcinoma and CRC: (a) less involvement of overexpression of mitotic genes in generating genomic instability in the colon and (b) the presence of CNA-suppressing pathways, including immune-surveillance, was only partly similar to those in the lung. Following 13 genes (TIGD6, TMED6, APOBEC3D, EP400NL, B3GNT4, ZNF683, FOXD4, FOXD4L1, PKIB, DDB2, MT1G, CLCN3, CAPS) were evaluated as potential drug development targets (hazard ratio [> 1.3 or < 0.5]). Identification of specific CRC genomic instability genes enables researchers to develop GI targeting approach. The new results suggest that the "targeting genomic instability and/or aneuploidy" approach must be tailored for specific organs.

Amino Acid Metabolism-Related lncRNA Signature Predicts the Prognosis of Breast Cancer

Background and Purpose: Breast cancer (BRCA) is the most frequent female malignancy and is potentially life threatening. The amino acid metabolism (AAM) has been shown to be strongly associated with the development and progression of human malignancies. In turn, long noncoding RNAs (lncRNAs) exert an important influence on the regulation of metabolism. Therefore, we attempted to build an AAM-related lncRNA prognostic model for BRCA and illustrate its immune characteristics and molecular mechanism.

Experimental Design: The RNA-seq data for BRCA from the TCGA-BRCA datasets were stochastically split into training and validation cohorts at a 3:1 ratio, to construct and validate the model, respectively. The amino acid metabolism-related genes were obtained from the Molecular Signature Database. A univariate Cox analysis, least absolute shrinkage and selection operator (LASSO) regression, and a multivariate Cox analysis were applied to create a predictive risk signature. Subsequently, the immune and molecular characteristics and the benefits of chemotherapeutic drugs in the high-risk and low-risk subgroups were examined.

Results: The prognostic model was developed based on the lncRNA group including LIPE-AS1, AC124067.4, LINC01655, AP005131.3, AC015802.3, USP30-AS1, SNHG26, and AL589765.4. Low-risk patients had a more favorable overall survival than did high-risk patients, in accordance with the results obtained for the validation cohort and the complete TCGA cohort. The elaborate results illustrated that a low-risk index was correlated with DNA-repair–associated pathways; a low TP53 and PIK3CA mutation rate; high infiltration of CD4⁺ T cells, CD8⁺ T cells, and M1 macrophages; active immunity; and less-aggressive phenotypes. In contrast, a high-risk index was correlated with cancer and metastasis-related pathways; a high PIK3CA and TP53 mutation rate; high infiltration of M0 macrophages, fibroblasts, and M2 macrophages; inhibition of the immune response; and more invasive phenotypes.

Conclusion: In conclusion, we attempted to shed light on the importance of AAM-associated lncRNAs in BRCA. The prognostic model built here might be acknowledged as an indispensable reference for predicting the outcome of patients with BRCA and help identify immune and molecular characteristics.

Construction of immune-related lncRNA signature to predict aggressiveness, immune landscape, and drug resistance of colon cancer

BACKGROUND

Colon cancer remains one of the most common malignancies across the world. Thus far, a biomarker, which can comprehensively predict the survival outcomes, clinical characteristics, and therapeutic sensitivity, is still lacking.

METHODS

We leveraged transcriptomic data of colon cancer from the existing datasets and constructed immune-related lncRNA (irlncRNA) pairs. After integrating with clinical survival data, we performed differential analysis and identified 11 irlncRNAs signature using Lasso regression analysis. We next plotted the 1-, 5-, and 10-year curve lines of receiver operating characteristics, calculated the areas under the curve, and recognized the optimal cutoff point. Then, we validated the pair-risk model in terms of the survival outcomes of the patients involved. Moreover, we tested the reliability of the model for predicting tumor aggressiveness and therapeutic susceptibility of colon cancer. Additionally, we reemployed the 11 of irlncRNAs involved in the pair-risk model to construct an expression-risk model to predict the prognostic outcomes of the patients involved.

RESULTS

We recognized a total of 377 differentially expressed irlncRNAs (DEirlcRNAs), including 28 low-expressed and 349 high-expressed irlncRNAs in colon cancer patients. After performing a univariant Cox analysis, we identified 115 risk irlncRNAs that were significantly correlated with survival outcomes of patients involved. By taking the overlap of the DEirlcRNAs and the risk irlncRNAs, we ultimately recognized 55 irlncRNAs as core irlncRNAs. Then, we established a Cox HR model (pair-risk model) as well as an expression HR model (exp-risk model) based on 11 of the 55 core irlncRNAs. We found that both of the two models significantly outperformed the commonly used clinical characteristics, including age, T, N, and M stages when predicting survival outcomes. Moreover, we validated the pair-risk model as a potential tool for studying the tumor microenvironment of colon cancer and drug susceptibility. Additionally, we noticed that combinational use of the pair-risk model and the exp-risk model yielded a more robust approach for predicting the survival outcomes of patients with colon cancer.

CONCLUSIONS

We recognized 11 irlncRNAs and created a pair-risk model and an exp-risk model, which have the potential to predict clinical characteristics of colon cancer, either solely or conjointly.