Construction of circRNA-miRNA-mRNA network and identification of novel potential biomarkers for non-small cell lung cancer

Identification of Circular RNAs as Biomarker Candidates in Lung Cancer Treatment

OBJECTIVE

Lung cancer is the most common malignancy and among the leading cause of cancer death worldwide. Therefore, there is an important need for biomarkers that can be used in the early diagnosis of the disease and in the follow-up of treatment. Circular RNAs (circRNAs) have a covalently closed circular structure that lacks 3' and 5' polar ends and is resistant to RNAase enzymes. Due to these properties, they can be stably found in body fluids. Therefore, they can serve as potential biomarkers in the diagnosis, monitoring of therapeutic response and prognosis of cancer. In our study, we aimed to investigate the expression levels of circRNA molecules in the treatment of lung cancer and to determine those that have the potential to be biomarkers.

METHODS

In this in vitro study, expression levels of 163 circRNAs were investigated in A549 cells, a non-small cell lung cancer cell line, before and after treatment with carboplatin and pemetrexed. Total RNA isolation and cDNA synthesis were performed after treatments. Expression levels of circRNA genes were determined by RT-qPCR method with the designed divergent primer sequences.

RESULTS

The study revealed the characterisation of differentially expressed circRNAs by treatment in lung cancer cells. Of them, hsa_circ_0001320 is not expressed in cancer cells, is expressed only after treatment, and increased the level of its expression in response to combination therapy.

CONCLUSION

As a result, while carboplatin, pemetrexed, and combined drug applications changed the expression levels of some circRNAs in lung cancer cells, some circRNAs were expressed only after treatment. In treatment follow-up and management, hsa_circ_0001320 has been identified as potential biomarker candidate.

Hsa_circ_0102899 promotes epithelial-mesenchymal transition in non-small cell lung cancer

OBJECTIVE

The significance of circular RNAs (circRNAs) has been identified in the progression of non-small cell lung cancer (NSCLC). Consistently, our study probed the functional actions of hsa_circ_0102899 (circ_0102899) in NSCLC cells.

METHODS

circ_0102899 expression was checked in NSCLC tissues, as well as its correlation with clinical characteristics of patients, Using A459 cells, transfection to alter gene expression was performed, thus measuring the changes of proliferation, apoptosis, migration, and apoptosis, as well as epithelial-mesenchymal transition (EMT)-related proteins. circ_0102899's effects in vivo were validated by tumor xenograft assay. Finally, the regulatory mechanism of circ_0102899 was investigated.

RESULTS

circ_0102899 indicated a high-expression level in NSCLC tissues which was associated with NSCLC tumor characteristics. Functionally, circ_0102899 knockdown not only inhibited the growth and EMT process of NSCLC cells, but also inhibited tumor formation in vivo. In terms of the regulatory mechanism, circ_0102899 had a binding to miR-885-5p to target eukaryotic translation initiation factor 4γ2 (EIF4G2). circ_0102899 mediated miR-885-5/EIF4G2 axis to accelerate the process of cell malignant behavior in NSCLC.

CONCLUSION

circ_0102899 promotes EMT and metastasis in NSCLC by regulating the miR-885-5p/EIF4G2 axis.

Inhibitory role of bone marrow mesenchymal stem cells-derived exosome in non-small-cell lung cancer: microRNA-30b-5p, EZH2 and PI3K/AKT pathway

Exosomal microRNA (miRNA) exerts potential roles in non-small-cell lung cancer (NSCLC). The current study elucidated the role of miR-30b-5p shuttled by bone marrow mesenchymal stem cells (BMSCs)-derived exosomes in treating NSCLC. Bioinformatics analysis was performed with NSCLC-related miRNA microarray GSE169587 and mRNA data GSE74706 obtained for collection of the differentially expressed miRNAs and mRNAs. The relationship between miR-30b-5p and EZH2 was predicted and confirmed. Exosomes were isolated from BMSCs and identified. BMSCs-derived exosomes overexpressing miR-30b-5p were used to establish subcutaneous tumorigenesis models to study the effects of miR-30b-5p, EZH2 and PI3K/AKT signalling pathway on tumour growth. A total of 86 BMSC-exo-miRNAs were differentially expressed in NSCLC. Bioinfomatics analysis found that BMSC-exo-miR-30b-5p could regulate NSCLC progression by targeting EZH2, which was verified by in vitro cell experiments. Besides, the target genes of miR-30b-5p were enriched in PI3K/AKT signalling pathway. Animal experiments validated that BMSC-exo-miR-30b-5p promoted NSCLC cell apoptosis and prevented tumorigenesis in nude mice via EZH2/PI3K/AKT axis. Collectively, the inhibitory role of BMSC-derived exosomes-loaded miR-30b-5p in NSCLC was achieved through blocking the EZH2/PI3K/AKT axis.

Hsa_circ_0041150 serves as a novel biomarker for monitoring chemotherapy resistance in small cell lung cancer patients treated with a first-line chemotherapy regimen

PURPOSE

To explore the potential of circRNAs as biomarkers in non-invasive body fluids for monitoring chemotherapy resistance in SCLC patients.

METHODS

CircRNAs were screened and characterized using transcriptome sequencing, Sanger sequencing, actinomycin D treatment, and Ribonuclease R assay. Our study involved 174 participants, and serum samples were collected from all chemotherapy-resistant patients (n = 54) at two time points: stable disease and progressive disease. We isolated and identified serum extracellular vesicles (EVs) from the patients using ultracentrifugation, transmission electron microscopy, nanoflow cytometry, and western blotting analysis. The expression levels of serum and serum EVs circRNAs were determined by quantitative real-time polymerase chain reaction (qRT-PCR). The impact of circRNA on the function of SCLC cells was assessed through various assays, including proliferation assay, scratch assay, transwell assay, and cisplatin resistance assay.

RESULTS

Hsa_circ_0041150 was found to be upregulated in chemoresistant SCLC cells and played a role in promoting proliferation, invasion, migration, and cisplatin resistance. Furthermore, the expression levels of hsa_circ_0041150 in serum and serum EVs increased when SCLC patients developed resistance after a first-line chemotherapy regimen. When combined with NSE, the monitoring sensitivity (70.37%) and specificity (81.48%) for chemotherapy resistance significantly improved. Moreover, the expression level of hsa_circ_0041150 showed significant associations with time to progression from SD to PD, and high hsa_circ_0041150 levels after drug resistance were more likely to cause chemotherapy resistance. Additionally, hsa_circ_0041150 demonstrated valuable potential in monitoring the progression from initial diagnosis to chemotherapy resistance in SCLC patients.

CONCLUSION

Thus, EVs hsa_circ_0041150 holds promise as a biomarker for monitoring chemotherapy resistance in SCLC patients.

Identification and Validation of Signature Genes and Potential Therapy Targets of Inflammatory Bowel Disease and Periodontitis

Background

Inflammatory bowel disease (IBD) and periodontitis (PD) are correlated, although the pathogenic mechanism behind their correlation has not been clarified. This study aims to explore the common signature genes and potential therapeutic targets of IBD and PD using transcriptomic analysis.

Methods

The GEO database was used to download datasets of IBD and PD, and differential expression analysis was used to identify DEGs. We then conducted GO and KEGG enrichment analyses of the shared genes. Next, we applied 4 machine learning (ML) algorithms (GLM, RF, GBM, and SVM) to select the best prediction model for diagnosing the disease and obtained the hub genes of IBD and PD. The diagnostic value of the signature genes was verified by a validation set and qRT‒PCR experiments. Subsequently, immune cell infiltration in IBD samples and PD samples was analyzed by ssGSEA. Finally, we investigated and validated the response of hub genes to infliximab therapy.

Results

We identified 43 upregulated genes as shared genes by intersecting the DEGs of IBD and PD. Functional enrichment analysis suggested that the shared genes were closely associated with immunity and inflammation. The ML algorithm and qRT‒PCR results indicated that IGKC and COL4A1 were the hub genes with the most diagnostic value for IBD and PD. Subsequently, through immune infiltration analysis, CD4 T cells, NK cells and neutrophils were identified to play crucial roles in the pathogenesis of IBD and PD. Finally, through in vivo and in vitro experiments, we found that IGKC and COL4A1 were significantly downregulated during the treatment of patients with IBD using infliximab.

Conclusion

We investigated the potential association between IBD and PD using transcriptomic analysis. The IGKC and COL4A1 genes were identified as characteristic genes and novel intervention targets for these two diseases. Infliximab may be used to treat or prevent IBD and PD.

Circular RNA PGPEP1 induces colorectal cancer malignancy and immune escape

OBJECTIVE

Colorectal cancer (CRC) is a prevalent gastrointestinal tumor globally. Circular RNAs (circRNAs) have been identified as regulatory players in the pathogenesis of CRC. However, it is unclear whether hsa_circ_0050102 (circPGPEP1) affects the malignant progression and immune escape in CRC.

METHODS

Bioinformatics analysis and circRNA in vivo precipitation experiments were performed to analyze and identify circRNAs that mediate immune escape in CRC. Using luciferase reporter assay, RIP, RNA pull-down assay, and FISH, the interaction between circPGPEP1, miR-515-5p, and nuclear factor of activated T-cell 5 (NFAT5) was identified. The functional role of circPGPEP1/miR-515-5p/NFAT5 axis in CRC anti-tumor immunity was investigated by co-culture assay, CFSE assay, and flow cytometry of CRC cells and T cells.

RESULTS

circPGPEP1 was a stable circRNA that was highly expressed in CRC. Functionally, circPGPEP1 silencing not only effectively inhibited CRC cell proliferation, migration, EMT, and immune escape and promoted apoptosis in vitro, but also inhibited CRC tumor growth and immune escape in vivo. In terms of the regulatory mechanism, circIGF2BP3 competitively upregulated NFAT5 expression by sponging miR-515-5p. Furthermore, functional rescue experiments showed that circPGPEP1 acted in CRC by regulating the miR-515-5p/NFAT5 axis.

CONCLUSION

Collectively, circPGPEP1 exerts an oncogene role in CRC by regulating the miR-515-5p/NFAT5 axis.

In silico analysis revealed the potential circRNA-miRNA-mRNA regulative network of non-small cell lung cancer (NSCLC)

BACKGROUND

The primary source of death in the world is non-small cell lung cancer (NSCLC). However, NSCLCs pathophysiology is still not completely understood. The current work sought to study the differential expression of mRNAs involved in NSCLC and their interactions with miRNAs and circRNAs.

METHODS

We utilized three microarray datasets (GSE21933, GSE27262, and GSE33532) from the GEO NCBI database to identify the differentially expressed genes (DEGs) in NSCLC. We employed DAVID Functional annotation tool to investigate the underlying GO biological process, molecular functions, and KEGG pathways involved in NSCLC. We performed the Protein-protein interaction (PPI) network, MCODE, and CytoHubba analysis from Cytoscape software to identify the significant DEGs in NSCLC. We utilized miRnet to anticipate and build interaction between miRNAs and mRNAs in NSCLC and ENCORI to predict the miRNA-circRNA relationships and build the ceRNA regulatory network. Finally, we executed the gene expression and Kaplan-Meier survival analysis to validate the significant DEGs in the ceRNA network utilizing TCGA NSCLC and GEPIA data.

RESULTS

We revealed a total of 156 overlapped DEGs (47 upregulated and 109 downregulated genes) in NSCLC. The PPI network, MCODE, and CytoHubba analysis revealed 12 hub genes (cdkn3, rrm2, ccnb1, aurka, nuf2, tyms, kif11, hmmr, ccnb2, nek2, anln, and birc5) that are associated with NSCLC. We identified that these 12 genes encode 12 mRNAs that are strongly linked with 8 miRNAs, and further, we revealed that 1 circRNA was associated with this 5 miRNA. We constructed the ceRNAs network that contained 1circRNA-5miRNAs-7mRNAs. The expression of these seven significant genes in LUAD & LUSC (NSCLC) was considerably higher in the TCGA database than in normal tissues. Kaplan-Meier survival plot reveals that increased expression of these hub genes was related to a poor survival rate in LUAD.

CONCLUSION

Overall, we developed a circRNA-miRNA-mRNA regulation network to study the probable mechanism of NSCLC.

Overview of MicroRNAs as Diagnostic and Prognostic Biomarkers for High-Incidence Cancers in 2021

MicroRNAs (miRNAs) are small non-coding RNAs (ncRNAs) about 22 nucleotides in size, which play an important role in gene regulation and are involved in almost all major cellular physiological processes. In recent years, the abnormal expression of miRNAs has been shown to be associated with human diseases including cancer. In the past ten years, the link between miRNAs and various cancers has been extensively studied, and the abnormal expression of miRNAs has been reported in various malignant tumors, such as lung cancer, gastric cancer, colorectal cancer, liver cancer, breast cancer, and prostate cancer. Due to the high malignancy grade of these cancers, it is more necessary to develop the related diagnostic and prognostic methods. According to the study of miRNAs, many potential cancer biomarkers have been proposed for the diagnosis and prognosis of diseases, especially cancer, thus providing a new theoretical basis and perspective for cancer screening. The use of miRNAs as biomarkers for diagnosis or prognosis of cancer has the advantages of being less invasive to patients, with better accuracy and lower price. In view of the important clinical significance of miRNAs in human cancer research, this article reviewed the research status of miRNAs in the above-mentioned cancers in 2021, especially in terms of diagnosis and prognosis, and provided some new perspectives and theoretical basis for the diagnosis and treatment of cancers.

Pan-cancer analysis identifies CDCA3 as a novel prognostic marker associated with immune infiltration in lung adenocarcinoma through bioinformatics analysis

Background

Lung adenocarcinoma (LUAD) is the most common subtype of lung malignancy. However, the expression of cell division cycle-associated protein-3 (CDCA3) and its significance in LUAD remain unclear. In this study, we investigated the functional role of CDCA3 in LUAD through bioinformatics analysis and expected to provide a new direction for clinical treatment.

Methods

The expression of CDCA3 was analyzed by online database. The association between the expression of CDCA3 and clinical parameters with LUAD was explored in TCGA. Survival and independent prognostic analysis were performed by TCGA database and the GSE30219 and GSE31210 datasets. Furthermore, Enrichment analyses were conducted to analyze the functions of CDCA3. Afterward, the relationship between CDCA3 and immune infiltration was investigated. Additionally, a competing endogenous RNA (ceRNA) regulatory network related to CDCA3 was constructed. Finally, CDCA3 expression was validated in clinical tissues by immunohistochemistry (IHC), real-time quantitative reverse transcriptase polymerase chain reaction (qRT-PCR), and western blotting (WB).

Results

CDCA3 expression was upregulated in 20 tumors and was significantly higher in LUAD compared with normal tissues in3 datasets. In addition, CDCA3 was significantly correlated with age, gender, stage, N, and smoking status. Kaplan-Meier survival curves showed that LUAD samples with higher CDCA3 expression were associated with poorer overall survival (OS) and disease-free survival (DFS). Univariate Cox regression analysis showed that the p value of CDCA3 expression was less than 0.05 (P<0.05) and it appeared in the results of multivariate Cox regression analysis (HR ≥1), indicating that CDCA3 can be used as an independent prognostic factor for LUAD. Intriguingly, Gene Set Enrichment Analysis (GSEA) suggested that CDCA3 was correlated with DNA-related terms and metabolic-related pathways in LUAD. CDCA3 expression was correlated with four immune scores and 14 immune cells in different groups. Next, a ceRNA network was constructed with CDCA3, and the experimental results of IHC, qRT-PCR, and WB were consistent with the bioinformatic analysis.

Conclusions

CDCA3 could serve as a prognostic biomarker for LUAD.

Identification of the Key miRNAs and Genes Associated with the Regulation of Non-Small Cell Lung Cancer: A Network-Based Approach

Lung cancer is the major cause of cancer-associated deaths across the world in both men and women. Lung cancer consists of two major clinicopathological categories, i.e., small cell lung cancer (SCLC) and non-small cell lung cancer (NSCLC). Lack of diagnosis of NSCLC at an early stage in addition to poor prognosis results in ineffective treatment, thus, biomarkers for appropriate diagnosis and exact prognosis of NSCLC need urgent attention. The proposed study aimed to reveal essential microRNAs (miRNAs) involved in the carcinogenesis of NSCLC that probably could act as potential biomarkers. The NSCLC-associated expression datasets revealed 12 differentially expressed miRNAs (DEMs). MiRNA-mRNA network identified key miRNAs and their associated genes, for which functional enrichment analysis was applied. Further, survival and validation analysis for key genes was performed and consequently transcription factors (TFs) were predicted. We obtained twelve miRNAs as common DEMs after assessment of all datasets. Further, four key miRNAs and nine key genes were extracted from significant modules based on the centrality approach. The key genes and miRNAs reported in our study might provide some information for potential biomarkers profitable to increased prognosis and diagnosis of lung cancer.

Clinical significance of cyclin-dependent kinase inhibitor 2C expression in cancers: from small cell lung carcinoma to pan-cancers

Background

Cyclin-dependent kinase inhibitor 2C (CDKN2C) was identified to participate in the occurrence and development of multiple cancers; however, its roles in small cell lung carcinoma (SCLC) remain unclear.

Methods

Differential expression analysis of CDKN2C between SCLC and non-SCLC were performed based on 937 samples from multiple centers. The prognosis effects of CDKN2C in patients with SCLC were detected using both Kaplan–Meier curves and log-rank tests. Using receiver-operating characteristic curves, whether CDKN2C expression made it feasible to distinguish SCLC was determined. The potential mechanisms of CDKN2C in SCLC were investigated by gene ontology terms and signaling pathways (Kyoto Encyclopedia of Genes and Genomes). Based on 10,080 samples, a pan-cancer analysis was also performed to determine the roles of CDKN2C in multiple cancers.

Results

For the first time, upregulated CDKN2C expression was detected in SCLC samples at both the mRNA and protein levels (p of Wilcoxon rank-sum test < 0.05; standardized mean difference = 2.86 [95% CI 2.20–3.52]). Transcription factor FOXA1 expression may positively regulate CDKN2C expression levels in SCLC. High CDKN2C expression levels were related to the poor prognosis of patients with SCLC (hazard ratio > 1, p < 0.05) and showed pronounced effects for distinguishing SCLC from non-SCLC (sensitivity, specificity, and area under the curve ≥ 0.95). CDKN2C expression may play a role in the development of SCLC by affecting the cell cycle. Furthermore, the first pan-cancer analysis revealed the differential expression of CDKN2C in 16 cancers (breast invasive carcinoma, etc.) and its independent prognostic significance in nine cancers (e.g., adrenocortical carcinoma). CDKN2C expression was related to the immune microenvironment, suggesting its potential usefulness as a prognostic marker in immunotherapy.

Conclusions

This study identified upregulated CDKN2C expression and its clinical significance in SCLC and other multiple cancers, suggesting its potential usefulness as a biomarker in treating and differentiating cancers.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12890-022-02036-5.

Different lipid metabolic profiles and their associated genes in sessile serrated adenoma or polyps compared to hyperplastic polyps

The serrated pathway is important in the development of colorectal cancer; currently, knowledge about the lipid metabolism profiles of serrated lesions is limited. Clinical characteristics were compared via Pearson's chi-squared test, nonparametric Kruskal-Wallis test and ANOVA. For some missing values, the MCAR test and multiple imputations were performed. Compared to patients with HP, the rates of younger patients (<50) and male patients with SSA or SSP were increased (P<0.05). Additionally, the BMI index and triglyceride levels were increased in patients with SSA or SSP. Inversely, patients with SSA/P had lower levels of HDL (P<0.05). Interestingly, the value of uric acid and tumor size in SSA/P patients tended to be greater than those in HP patients, and the ratio of patients who smoked was also increased. Other characteristics, such as LDL, ALB, γ-GT, and the N/L ratio, were similar among the subtypes of serrated lesions. Analysis of GEO data (GSE43841) showed that 9 genes were associated with lipid metabolism, including ADRB3, DEGS2, PRKACB, SLC44A1, and CA4. PRKACB was downregulated in SSA/P tissue compared to HP tissue samples from the GSE76987 dataset and our hospital. In conclusion, compared to benign HP, lower HDL levels and higher triglyceride levels tended to occur in CRC precursor SSA/P lesions, and these factors may be associated with metabolic genomic markers, such as PRKACB.