Identification of novel biomarkers affecting the metastasis of colorectal cancer through bioinformatics analysis and validation through qRT-PCR.

BACKGROUND: Tumor progression and distant metastasis are the main causes of deaths in colorectal cancer (CRC) patients, and the molecular mechanisms in CRC metastasis have not been completely discovered. METHODS: We identified differentially expressed genes (DEGs) and lncRNAs (DELs) of CRC from The Cancer Genome Atlas (TCGA) database. Then we conducted the weighted gene co-expression network analysis (WGCNA) to investigate co-expression modules related with CRC metastasis. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, DEG-DEL co-expression network and survival analyses of significant modules were also conducted. Finally, the expressions of selected biomarkers were validated in cell lines by quantitative real-time PCR (qRT-PCR). RESULTS: 2032 DEGs and 487 DELs were involved the construction of WGCNA network, and greenyellow, turquoise and brown module were identified to have more significant correlation with CRC metastasis. GO and KEGG pathway analysis of these three modules have proven that the functions of DEGs were closely involved in many important processes in cancer pathogenesis. Through the DEG-DEL co-expression network, 12 DEGs and 2 DELs were considered as hub nodes. Besides, survival analysis showed that 30 DEGs were associated with the overall survival of CRC. Then 10 candidate biomarkers were chosen for validation and the expression of CA2, CHP2, SULT1B1, MOGAT2 and C1orf115 were significantly decreased in CRC cell lines when compared to normal human colonic epithelial cells, which were consistent with the results of differential expression analysis. Especially, low expression of SULT1B1, MOGAT2 and C1orf115 were closely correlated with poorer survival of CRC. CONCLUSION: This study identified 5 genes as new biomarkers affecting the metastasis of CRC. Besides, SULT1B1, MOGAT2 and C1orf115 might be implicated in the prognosis of CRC patients.

Integrated analysis of dysregulated long non-coding RNAs/microRNAs/mRNAs in metastasis of lung adenocarcinoma.

BACKGROUND: Lung adenocarcinoma (LUAD), largely remains a primary cause of cancer-related death worldwide. The molecular mechanisms in LUAD metastasis have not been completely uncovered. METHODS: In this study, we identified differentially expressed genes (DEGs), miRNAs (DEMs) and lncRNAs (DELs) underlying metastasis of LUAD from The Cancer Genome Atlas database. Intersection mRNAs were used to perform gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and co-expression network analysis. In addition, survival analyses of intersection mRNAs were conducted. Finally, intersection mRNAs, miRNAs and lncRNAs were subjected to construct miRNA-mRNA-lncRNA network. RESULTS: A total of 1015 DEGs, 54 DEMs and 22 DELs were identified in LUAD metastasis and non-metastasis samples. GO and KEGG pathway analysis had proven that the functions of intersection mRNAs were closely related with many important processes in cancer pathogenesis. Among the co-expression interactions network, 22 genes in the co-expression network were over the degree 20. These genes imply that they have connections with many other gene nodes. In addition, 14 target genes (ARHGAP11A, ASPM, HELLS, PRC1, TMPO, ARHGAP30, CD52, IL16, IRF8, P2RY13, PRKCB, PTPRC, SASH3 and TRAF3IP3) were found to be associated with survival in patients with LUAD significantly (log-rank P < 0.05). Two lncRNAs (LOC96610 and ADAM6) acting as ceRNAs were identified based on the miRNA-mRNA-lncRNA network. CONCLUSIONS: Taken together, the results may provide a novel perspective to develop a multiple gene diagnostic tool for LUAD prognosis, which might also provide potential biomarkers or therapeutic targets for LUAD.

The prognostic value and response to immunotherapy of immunogenic cell death-associated genes in breast cancer.

Breast cancer (BRCA) remains the most prevalent cancer worldwide and the tumor microenvironment (TME) has been discovered to exert a wide influence on the overall survival and therapeutic response. Numerous lines of evidence reported that the effects of immunotherapy of BRCA were manipulated by TME. Immunogenic cell death (ICD) is a form of regulated cell death (RCD) that is capable of fueling adaptive immune responses and aberrant expression of ICD-related genes (ICDRGs) can govern the TME system by emitting danger signals or damage-associated molecular patterns (DAMPs). In the current study, we obtained 34 key ICDRGs in BRCA. Subsequently, using the transcriptome data of BRCA from the TCGA database, we constructed a risk signature based on 6 vital ICDRGs, which had a good performance in predicting the overall survival of BRCA patients. We also examined the efficacy of our risk signature in the validation dataset (GSE20711) in the GEO database and it performed excellently. According to the risk model, patients with BRCA were divided into high-risk and low-risk groups. Also, the unique immune characteristics and TME between the two subgroups and 10 promising small molecule drugs targeting BRCA patients with different ICDRGs risk have been investigated. The low-risk group had good immunity indicated by T cell infiltration and high immune checkpoint expression. Moreover, the BRCA samples could be divided into three immune subtypes according to immune response severity (ISA, ISB, and ISC). ISA and ISB predominated in the low-risk group and patients in the low-risk group exhibited a more vigorous immune response. In conclusion, we developed an ICDRGs-based risk signature that can predict the prognosis of BRCA patients and offer a novel therapeutic strategy for immunotherapy, which would be of great significance in the BRCA clinical setting.

Immune Infiltration-Related ceRNA Network Revealing Potential Biomarkers for Prognosis of Head and Neck Squamous Cell Carcinoma.

Background: Head and neck squamous cell carcinoma (HNSCC) is a frequently lethal malignancy, and the mortality is considerably high. The tumor microenvironment (TME) has been identified as a critical participation in cancer development, treatment, and prognosis. However, competing endogenous RNA (ceRNA) networks grouping with immune/stromal scores of HNSCC patients need to be further illustrated. Therefore, our study aimed to provide clues for searching promising prognostic markers of TME in HNSCC. Materials and Methods: ESTIMATE algorithm was used to calculate immune scores and stromal scores of the enrolled HNSCC patients. Differentially expressed genes (DEGs), lncRNAs (DELs), and miRNAs (DEMs) were identified by comparing the expression difference between high and low immune/stromal scores. Then, a ceRNA network and protein-protein interaction (PPI) network were constructed for selecting hub regulators. In addition, survival analysis was performed to access the association between immune scores, stromal scores, and differentially expressed RNAs in the ceRNA network and the overall survival (OS) of HNSCC patients. Then, the GSE65858 datasets from Gene Expression Omnibus (GEO) database was used for verification. At last, the difference between the clinical characteristics and immune cell infiltration in different expression groups of IL10RA, PRF1, and IL2RA was analyzed. Results: Survival analysis showed a better OS in the high immune score group, and then we constructed a ceRNA network composed of 97 DEGs, 79 DELs and 22 DEMs. Within the ceRNA network, FOXP3, IL10RA, STAT5A, PRF1, IL2RA, miR-148a-3p, miR-3065-3p, and lncRNAs, including CXCR2P1, HNRNPA1P21, CTA-384D8.36, and IGHV1OR15-2, were closely correlated with the OS of HNSCC patients. Especially, using the data from GSE65858, we successfully verified that IL10RA, PRF1, and IL2RA were not only significantly upregulated in patients high immune scores, but also their high expressions were associated with longer survival time. In addition, stratified analysis showed that PRF1 and IL2RA might be involved in the mechanism of tumor progress. Conclusion: In conclusion, we constructed a ceRNA network related to the TME of HNSCC, which provides candidates for therapeutic intervention and prognosis evaluation.

Comprehensive Characterization of the Function of Metabolic Genes and Establishment of a Prediction Model in Breast Cancer.

Background: Breast cancer (BC) is a highly heterogeneous disease with high morbidity and mortality. Its subtypes may have distinctly different biological behaviors, clinical outcomes, and therapeutic responses. The metabolic status of BC tissue is closely related to its progress. Therefore, we comprehensively characterized the function of metabolic genes in BC and identified new biomarkers to predict BC patients' prognoses. Methods: Metabolic genes were identified by intersecting genes obtained from two published pieces of literature. The function of metabolic genes in BC was determined by extracting differentially expressed genes (DEGs), performing functional enrichment analyses, analyzing the infiltrating proportion of immune cells, and conducting metabolic subgroup analyses. A risk score model was constructed to assess the prognoses of BC patients by performing the univariate Cox regression, LASSO algorithm, multivariate Cox regression, Kaplan-Meier survival analyses, and ROC curve analyses in the training set. The prognostic model was then validated on the testing dataset, external dataset, the whole TCGA-BC database, and our clinical specimens. Finally, a nomogram was constructed for clinical prognostic prediction based on the risk score model and other clinicopathological parameters. Results: 955 metabolic genes were obtained. Among these, 157 metabolic DEGs were identified between BC and normal tissues for subsequent GO and KEGG pathway enrichment analyses. 5 metabolic genes were negatively correlated with CD8+ T cells, while 49 genes were positively correlated with CD8+ T cells. Furthermore, 5 metabolic subgroups with varying proportions of PAM50 subtypes, TNM classification, and immune cell infiltration were obtained. Finally, a risk score model was constructed to predict the prognoses of BC patients, and a nomogram incorporating the risk score model was established for clinical application. Conclusion: In this study, we elucidated tumor heterogeneity from metabolite profiling of BC. The roles of metabolic genes in the occurrence of BC were comprehensively characterized, clarifying the relationship between the tumor microenvironment (TME) and metabolic genes. Meanwhile, a concise prediction model was also constructed based on metabolic genes, providing a convenient and precise method for the individualized diagnosis and treatment of BC patients.

Identification and validation of immunotherapy for four novel clusters of colorectal cancer based on the tumor microenvironment.

The incidence and mortality of colorectal cancer (CRC) are increasing year by year. The accurate classification of CRC can realize the purpose of personalized and precise treatment for patients. The tumor microenvironment (TME) plays an important role in the malignant progression and immunotherapy of CRC. An in-depth understanding of the clusters based on the TME is of great significance for the discovery of new therapeutic targets for CRC. We extracted data on CRC, including gene expression profile, DNA methylation array, somatic mutations, clinicopathological information, and copy number variation (CNV), from The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO) (four datasets-GSE14333, GSE17538, GSE38832, and GSE39582), cBioPortal, and FireBrowse. The MCPcounter was utilized to quantify the abundance of 10 TME cells for CRC samples. Cluster repetitive analysis was based on the Hcluster function of the Pheatmap package in R. The ESTIMATE package was applied to compute immune and stromal scores for CRC patients. PCA analysis was used to remove batch effects among different datasets and transform genome-wide DNA methylation profiling into methylation of tumor-infiltrating lymphocyte (MeTIL). We evaluated the mutation differences of the clusters using MOVICS, DeconstructSigs, and GISTIC packages. As for therapy, TIDE and SubMap analyses were carried out to forecast the immunotherapy response of the clusters, and chemotherapeutic sensibility was estimated based on the pRRophetic package. All results were verified in the TCGA and GEO data. Four immune clusters (ImmClust-CS1, ImmClust-CS2, ImmClust-CS3, and ImmClust-CS4) were identified for CRC. The four ImmClusts exhibited distinct TME compositions, cancer-associated fibroblasts (CAFs), functional orientation, and immune checkpoints. The highest immune, stromal, and MeTIL scores were observed in CS2, in contrast to the lowest scores in CS4. CS1 may respond to immunotherapy, while CS2 may respond to immunotherapy after anti-CAFs. Among the four ImmClusts, the top 15 markers with the highest mutation frequency were acquired, and CS1 had significantly lower CNA on the focal level than other subtypes. In addition, CS1 and CS2 patients had more stable chromosomes than CS3 and CS4. The most sensitive chemotherapeutic agents in these four ImmClusts were also found. IHC results revealed that CD29 stained significantly darker in the cancer samples, indicating that their CD29 was highly expressed in colon cancer. This work revealed the novel clusters based on TME for CRC, which would guide in predicting the prognosis, biological features, and appropriate treatment for patients with CRC.

FABP6 Expression Correlates with Immune Infiltration and Immunogenicity in Colorectal Cancer Cells.

Background: Immune checkpoint inhibitors (ICIs) have rapidly revolutionized colorectal cancer (CRC) treatment, but resistance caused by the heterogeneous tumor microenvironment (TME) still presents a challenge. Therefore, it is necessary to characterize TME immune infiltration and explore new targets to improve immunotherapy. Methods: The compositions of 64 types of infiltrating immune cells and their relationships with CRC patient clinical characteristics were assessed. Differentially expressed genes (DEGs) between "hot" and "cold" tumors were used for functional analysis. A prediction model was constructed to explore the survival of CRC patients treated with and without immunotherapy. Finally, fatty acid-binding protein (FABP6) was selected for in vitro experiments, which revealed its roles in the proliferation, apoptosis, migration, and immunogenicity of CRC tissues and cell lines. Results: The infiltration levels of several immune cells were associated with CRC tumor stage and prognosis. Different cell types showed the synergistic or antagonism infiltration patterns. Enrichment analysis of DEGs revealed that immune-related signaling was significantly activated in hot tumors, while metabolic process pathways were altered in cold tumors. In addition, the constructed model effectively predicted the survival of CRC patients treated with and without immunotherapy. FABP6 knockdown did not significantly alter tumor cell proliferation, apoptosis, and migration. FABP6 was negatively correlated with immune infiltration, and knockdown of FABP6 increased major histocompatibility complex (MHC) class 1 expression and promoted immune-related chemokine secretion, indicating the immunogenicity enhancement of tumor cells. Finally, knockdown of FABP6 could promote the recruitment of CD8+ T cells. Conclusion: Collectively, we described the landscape of immune infiltration in CRC and identified FABP6 as a potential immunotherapeutic target for treatment.

The prognostic significance of high/positive expression of cyclin B1 in patients with three common digestive cancers: a systematic review and meta-analysis.

Recently, several studies have reported that the expression of cyclin B1 may be associated with the prognosis of cancer. Nevertheless, their conclusions were still controversial. The present was designed to analyze and evaluate the prognostic role of cyclin B1 expression in patients with digestive cancer. PubMed, Embase, Cochrane Library and Web of Science were searched to January, 2017. Pooled odds ratio (OR) with 95% confidence intervals (CIs) were estimated. For the pooled OR estimates of OS, we performed subgroup analysis. Besides, sensitivity analysis was performed to examine the stability of the combined results. All statistical analyses were performed using standard statistical procedures provided in RevMan 5.2. A total of 12 studies (N = 2080 participants) were included for this meta-analysis. The positive/high expression of cyclin B1 had an obvious association with both 3-year overall survival (OR 0.21, 95% CI 0.12-0.37; P < 0.00001) and 5-year overall survival (OR 0.20, 95% CI 0.12-0.34; P < 0.00001) in esophageal cancer, and 5-year overall survival of colorectal cancer (OR 2.01, 95% CI 1.32-3.08; P = 0.001). This meta-analysis indicated that positive/high expression of cyclin B1 may have a close association with worse survival in patients with esophageal cancer, but better prognosis in patients with colorectal cancer.

Self-assembled monolayers of DNA on cysteamine modified Au(111) surface: Atomic force microscopy study.

The lambda-DNA molecules self-assemble on cysteamine-modified gold (111) surface to form flat-lying self-assembled monolayers (SAMs). The formation kinetics of such DNA SAMs is studied by atomic force microscopy (AFM). AFM results show that DNA molecules do not arrange themselves on cysteamine-modified gold (111) surface into a well-ordered monolayer. It is also found that the surface density of DNA monolayer does not increase as the DNA concentration increases. The high temperature of DNA solution and the immersing in ultrapure water produce some obvious DNA bundles. Whereas divalent cations in DNA solution result in the formation of more compact DNA films. The obtained information may be useful for practical application of the DNA films and further theoretical studies.