No abstract

Objective. The aim of the study is to explore the relationship between lymphatic metastasis genes, prognosis, and immune cell infiltration in patients with colon cancer. Methods. Based on the Cancer Genome Atlas Program (TCGA) database, differentially expressed genes and prognostic genes related to colon adenocarcinoma (COAD) lymphatic metastasis were screened and intersected. We used lasso and univariate Cox regression analysis to screen core genes and establish a preliminary prediction model. GO and KEGG enrichment analysis was used for lymphatic metastasis-related genes, and single GSEA was used for the final screening results. Finally, we evaluated the relationship between identified genes and immune cell infiltration. Results. A total of 1727 genes were differentially expressed between COAD patients with TNM stages of N0 and N1. After further screening, six core genes (RNU4-2, ZNF556, RNVU1-15, NSA2P6, RN7SL767P, and RN7SL473P) were obtained, and a preliminary prediction model was established, in which ZNF556 was a risk factor, and the rest were protective factors. Single GSEA showed that pathways such as systemic lupus erythematosus might play an important role in the initial lymphatic metastasis of COAD. GO and KEGG enrichment analysis of 1727 genes supported this result. Immune infiltration analysis showed that six genes were significantly correlated with T cell and NK cell families. Conclusion. Six core genes may affect COAD initial lymphatic metastasis through the systemic lupus erythematosus pathway and immune cell infiltration.

N6-methyladenosine (m6A) modification is a common epigenetic modification. It is reported that lncRNA can be regulated by m6A modification. Previous studies have shown that lncRNAs associated with m6A regulation (m6A-lncRNAs) serve as ideal prognostic biomarkers. However, whether lncRNAs are involved in m6A modification in colon adenocarcinoma (COAD) needs further exploration. The objective of this study was to construct an m6A-lncRNAs-based signature for patients with COAD. We obtained the RNA sequencing data and clinical information from The Cancer Genome Atlas (TCGA). Pearson correlation analysis was employed to recognize lncRNAs associated with m6A regulation (m6A-lncRNAs). 24 prognostic m6A-lncRNAs was identified by univariate Cox regression analysis. Gene set enrichment analysis (GSAE) was used to investigate the potential cellular pathways and biological processes. We have also explored the relationship between immune infiltrate levels and m6A-lncRNAs. Then, a predictive signature based on the expression of 13 m6A-lncRNAs was constructed by the Lasso regression algorithm, including UBA6-AS1, AC139149.1, U91328.1, AC138207.5, AC025171.4, AC008760.1, ITGB1-DT, AP001619.1, AL391422.4, AC104532.2, ZEB1-AS1, AC156455.1, and AC104819.3. ROC curves and K M survival curves have shown that the risk score has a well-predictive ability. We also set up a quantitative nomogram on the basis of risk score and prognosis-related clinical characteristics. In summary, we have identified some m6A-lncRNAs that correlated with prognosis and tumor immune microenvironment in COAD. In addition, a potential alternative signature based on the expression of m6A-lncRNAs was provided for the management of COAD patients.