An integrated approach for mining precise RNA-based cervical cancer staging biomarkers

The microRNA Let-7 and its exosomal form: Epigenetic regulators of gynecological cancers

Many types of gynecological cancer (GC) are often silent until they reach an advanced stage, and are therefore often diagnosed too late for effective treatment. Hence, there is a real need for more efficient diagnosis and treatment for patients with GC. During recent years, researchers have increasingly studied the impact of microRNAs cancer development, leading to a number of applications in detection and treatment. MicroRNAs are a particular group of tiny RNA molecules that regulate regular gene expression by affecting the translation process. The downregulation of numerous miRNAs has been observed in human malignancies. Let-7 is an example of a miRNA that controls cellular processes as well as signaling cascades to affect post-transcriptional gene expression. Recent research supports the hypothesis that enhancing let-7 expression in those cancers where it is downregulated may be a potential treatment option. Exosomes are tiny vesicles that move through body fluids and can include components like miRNAs (including let-7) that are important for communication between cells. Studies proved that exosomes are able to enhance tumor growth, angiogenesis, chemoresistance, metastasis, and immune evasion, thus suggesting their importance in GC management.

Multidimensional computational study to understand non-coding RNA interactions in breast cancer metastasis

Metastasis is a major breast cancer hallmark due to which tumor cells tend to relocate to regional or distant organs from their organ of origin. This study is aimed to decipher the interaction among 113 differentially expressed genes, interacting non-coding RNAs and drugs (614 miRNAs, 220 lncRNAs and 3241 interacting drugs) associated with metastasis in breast cancer. For an extensive understanding of genetic interactions in the diseased state, a backbone gene co-expression network was constructed. Further, the mRNA-miRNA-lncRNA-drug interaction network was constructed to identify the top hub RNAs, significant cliques and topological parameters associated with differentially expressed genes. Then, the mRNAs from the top two subnetworks constructed are considered for transcription factor (TF) analysis. 39 interacting miRNAs and 1641 corresponding TFs for the eight mRNAs from the subnetworks are also utilized to construct an mRNA-miRNA-TF interaction network. TF analysis revealed two TFs (EST1 and SP1) from the cliques to be significant. TCGA expression analysis of miRNAs and lncRNAs as well as subclass-based and promoter methylation-based expression, oncoprint and survival analysis of the mRNAs are also done. Finally, functional enrichment of mRNAs is also performed. Significant cliques identified in the study can be utilized for identification of newer therapeutic interventions for breast cancer. This work will also help to gain a deeper insight into the complicated molecular intricacies to reveal the potential biomarkers involved with breast cancer progression in future.

Exploring the methylation status of CFTR and PKIA genes as potential biomarkers for lung adenocarcinoma

BACKGROUND

One of the most prevalent cancers in the world is lung cancer, with adenocarcinoma (LUAD) making up a significant portion of cases. According to the National Cancer Institute (NCI), there are new cases and fatality rates per 100,000 individuals as follows: New instances of lung and bronchial cancer occur annually at a rate of 50.0 per 100,000 persons. The yearly death rate for men and women is 35.0 per 100,000. DNA methylation is one of the earliest discovered and widely studied epigenetic regulatory mechanisms, and its abnormality is closely related to the occurrence and development of cancer. However, the prognostic value of DNA methylation and LUAD needs to be further explored to improve the survival prediction of LUAD patients.

METHODS

The transcriptome data and clinical data of LUAD were downloaded from TCGA and GEO databases, and the Illumina Human Methylation450 array (450k array) data were downloaded from the TCGA database. Firstly, the intersection of the expressed genes of the two databases is corrected, the differential analysis is performed, and the methylation data is evaluated by the MethylMix package to obtain differentially methylated genes. Independent prognostic genes were screened out using univariate and multivariate Cox regression analysis, and a methylation prognostic model was developed using univariate Cox analysis and validated with the GSE30219 dataset in the GEO database. Survival analysis between methylation high-risk and low-risk groups was performed and a methylation-based gene prognostic model was constructed. Finally, the prediction of potential drugs associated with the LUAD gene signature using Drug Sensitivity Genomics in Cancer (GDSC).

RESULTS

In this study, a total of 555 samples from the TCGA database and 307 samples from GSE30219 were included, and a total of 24 differential methylation driver genes were identified. Univariate and multivariate Cox regression analyzes were used to screen out independent prognostic genes, involving 2 genes: CFTR, PKIA. Survival analysis was different between the methylation high-risk group and the low-risk group, the CFTR high methylation group and the low methylation group were poor, and the opposite was true for PKIA.

CONCLUSIONS

Our study revealed that the methylation status of CFTR and PKIA can serve as potential prognostic biomarkers and therapeutic targets in lung cancer.

Antitumor activity of miR-188-3p in gastric cancer is achieved by targeting CBL expression and inactivating the AKT/mTOR signaling

BACKGROUND

Altered miR-188-3p expression has been observed in various human cancers.

AIM

To investigate the miR-188-3p expression, its roles, and underlying molecular events in gastric cancer.

METHODS

Fifty gastric cancer and paired normal tissues were collected to analyze miR-188-3p and CBL expression. Normal and gastric cancer cells were used to manipulate miR-188-3p and CBL expression through different assays. The relationship between miR-188-3p and CBL was predicted bioinformatically and confirmed using a luciferase gene reporter assay. A Kaplan-Meier analysis was used to associate miR-188-3p or CBL expression with patient survival. A nude mouse tumor cell xenograft assay was used to confirm the in vitro data.

RESULTS

MiR-188-3p was found to be lower in the plasma of gastric cancer patients, tissues, and cell lines compared to their healthy counterparts. It was associated with overall survival of gastric cancer patients (P < 0.001), tumor differentiation (P < 0.001), lymph node metastasis (P = 0.033), tumor node metastasis stage (I/II vs III/IV, P = 0.024), and American Joint Committee on Cancer stage (I/II vs III/IV, P = 0.03). Transfection with miR-188-3p mimics reduced tumor cell growth and invasion while inducing apoptosis and autophagy. CBL was identified as a direct target of miR-188-3p, with its expression antagonizing the effects of miR-188-3p on gastric cancer (GC) cell proliferation by inducing tumor cell apoptosis and autophagy through the inactivation of the Akt/mTOR signaling pathway. The in vivo data confirmed antitumor activity via CBL downregulation in gastric cancer.

CONCLUSION

The current data provides ex vivo, in vitro, and in vivo evidence that miR-188-3p acts as a tumor suppressor gene or possesses antitumor activity in GC.

Single-nucleus RNA sequencing and deep tissue proteomics reveal distinct tumour microenvironment in stage-I and II cervical cancer

BACKGROUND

Cervical cancer (CC) is the 3rd most common cancer in women and the 4th leading cause of deaths in gynaecological malignancies, yet the exact progression of CC is inconclusive, mainly due to the high complexity of the changing tumour microenvironment (TME) at different stages of tumorigenesis. Importantly, a detailed comparative single-nucleus transcriptomic analysis of tumour microenvironment (TME) of CC patients at different stages is lacking.

METHODS

In this study, a total of 42,928 and 29,200 nuclei isolated from the tumour tissues of stage-I and II CC patients and subjected to single-nucleus RNA sequencing (snRNA-seq) analysis. The cell heterogeneity and functions were comparatively investigated using bioinformatic tools. In addition, label-free quantitative mass spectrometry based proteomic analysis was carried out. The proteome profiles of stage-I and II CC patients were compared, and an integrative analysis with the snRNA-seq was performed.

RESULTS

Compared with the stage-I CC (CCI) patients, the immune response relevant signalling pathways were largely suppressed in various immune cells of the stage-II CC (CCII) patients, yet the signalling associated with cell and tissue development was enriched, as well as metabolism for energy production suggested by the upregulation of genes associated with mitochondria. This was consistent with the quantitative proteomic analysis that showed the dominance of proteins promoting cell growth and intercellular matrix development in the TME of CCII group. The interferon-α and γ responses appeared the most activated pathways in many cell populations of the CCI patients. Several collagens, such as COL12A1, COL5A1, COL4A1 and COL4A2, were found significantly upregulated in the CCII group, suggesting their roles in diagnosing CC progression. A novel transcript AC244205.1 was detected as the most upregulated gene in CCII patients, and its possible mechanistic role in CC may be investigated further.

CONCLUSIONS

Our study provides important resources for decoding the progression of CC and set the foundation for developing novel approaches for diagnosing CC and tackling the immunosuppressive TME.

An integrated approach to understand fluid shear stress-driven and reactive oxygen species-mediated metastasis of colon adenocarcinoma through mRNA-miRNA-lncRNA-circRNA networks

Development of colon adenocarcinoma (COAD) metastasis involves several mediators including fluid shear stress (FSS), intracellular ROS levels, and non-coding RNAs. In our present study, we identified and investigated the role of regulatory non-coding RNA molecules specifically involved in COAD metastasis and their association with FSS and ROS. Interactions between the mRNAs associated with FSS and ROS, the corresponding microRNAs (miRNAs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) in COAD metastasis were used to generate the mRNA-miRNA-lncRNA-circRNA network. Experimental validation of the identified RNA hubs using quantitative real-time PCR demonstrated a direct effect of the FSS on their expression levels in cancer cells. FSS resulted in the downregulation of HMGA1 and RAN, as well as the upregulation of HSP90AA1, PMAIP1 and BIRC5. Application of shear stress also led to downregulation of hsa-miR-26b-5p and hsa-miR-34a-5p levels in HCT116 cells. Further, functional enrichment and survival analysis of the significant miRNAs, as well as the OncoPrint and the survival analyses of the selected mRNAs were performed. Subsequently, their functional role was also corroborated with existing literature. Ten significant miRNA hubs were identified, out of which hsa-miR-17-5p and hsa-miR-20a-5p were found to interact with lncRNA (CCAT2) while hsa-miR-335 was found to interact with four circRNAs. Fifteen significant miRNAs were identified in 10 different modules suggesting their importance in FSS and ROS-mediated COAD metastasis. Finally, 10 miRNAs and 3 mRNAs associated with FSS and/or ROS were identified as significant overall survival markers; 33 mRNAs were also identified as metastasis-free survival markers whereas 15 mRNAs showed > 10% gene alterations in TCGA-COAD data and may serve as promising therapeutic biomarkers in the COAD metastasis.

Construction of competing endogenous RNA interaction network as prognostic markers in metastatic melanoma

Malignant melanoma (MM) is a malignant tumor originating from melanocytes, with high aggressiveness, high metastasis and extremely poor prognosis. MM accounts for 4% of skin cancers and 80% of mortality, and the median survival of patients with metastatic melanoma is only about 6 months, with a five-year survival rate of less than 10%. In recent years, the incidence of melanoma has gradually increased and has become one of the serious diseases that endanger human health. Competitive endogenous RNA (ceRNA) is the main model of the mechanism by which long chain non-coding RNAs (lncRNAs) play a regulatory role in the disease. LncRNAs can act as a "sponge", competitively attracting small RNAs (micoRNAs; miRNAs), thus interfering with miRNA function, and affect the expression of target gene messenger RNAs (mRNAs), ultimately promoting tumorigenesis and progression. Bioinformatics analysis can identify potentially prognostic and therapeutically relevant differentially expressed genes in MM, finding lncRNAs, miRNAs and mRNAs that are interconnected through the ceRNA network, providing further insight into gene regulation and prognosis of metastatic melanoma. Weighted co-expression networks were used to identify lncRNA and mRNA modules associated with the metastatic phenotype, as well as the co-expression genes contained in the modules. A total of 17 lncRNAs, six miRNAs, and 11 mRNAs were used to construct a ceRNA interaction network that plays a regulatory role in metastatic melanoma patients. The prognostic risk model was used as a sorter to classify the survival prognosis of melanoma patients. Four groups of ceRNA interaction triplets were finally obtained, which miR-3662 might has potential implication for the treatment of metaststic melanoma patients, and futher experiments confirmed the regulating relationship and phenotype of this assumption. This study provides new targets to regulate metastatic process, predict metastatic potential and indicates that the miR-3662 can be used in the treatment of melanoma.

Construction of circRNA-based ceRNA network and its prognosis-associated subnet of clear cell renal cell carcinoma

Circular RNAs (circRNAs) are novel biomarkers of various cancers. CircRNAs can sponge miRNAs and regulate target mRNAs, which was called competing endogenous RNAs (ceRNA). This study was designed to identify circRNAs related to patients with clear cell renal cell carcinoma (ccRCC) and the first to select three independent Gene Expression Omnibus microarrays covering circRNAs, miRNAs, and mRNAs for multiple analyses. The data of clinical cases applied in our study were obtained from The Cancer Genome Atlas. We successfully conducted a circRNA/miRNA/mRNA ceRNA network related to ccRCC patients via R software and Cytoscape including 8 circRNAs, 6 miRNAs, and 49 mRNAs. The prognosis‐associated subnet covered 8 circRNAs, 6 miRNAs, and 22 mRNAs. Quantitative real‐time PCR was applied to measure our prediction in three renal cell lines and 23 pairs of tissues. Small interfering RNA targeting the back‐splice region of hsa_circ_0001167 was further implied to confirm the regulation. Ultimately, hsa_circ_0001167/hsa‐miR‐595/CCDC8 regulatory axis was identified in this study, which may serve as prognostic indicators. Lower levels of hsa_circ_0001167 and CCDC8 were potentially correlated with worse patient survival.

Profiling of Serum Exosome MiRNA Reveals the Potential of a MiRNA Panel as Diagnostic Biomarker for Alzheimer's Disease

Alzheimer's disease (AD) is the most common neurodegenerative disease in the older adults. Although much effort has been made in the analyses of diagnostic biomarkers, such as amyloid-β, tau, and neurofilament light chain, identifying peripheral blood-based biomarkers is in extremely urgent need for their minimal invasiveness and more convenience. Here we characterized the miRNA profile by RNA sequencing in human serum exosomes from AD patients and healthy controls (HC) to investigate its potential for AD diagnosis. Subsequently, Gene Ontology analysis and pathway analysis were performed for the targeted genes from the differentially expressed miRNAs. These basic functions were differentially enriched, including cell adhesion, regulation of transcription, and the ubiquitin system. Functional network analysis highlighted the pathways of proteoglycans in cancer, viral carcinogenesis, signaling pathways regulating pluripotency of stem cells, and cellular senescence in AD. A total of 24 miRNAs showed significantly differential expression between AD and HC with more than ± 2.0-fold change at p value < 0.05 and at least 50 reads for each sample. Logistic regression analysis established a model for AD prediction by serum exosomal miR-30b-5p, miR-22-3p, and miR-378a-3p. Sequencing results were validated using quantitative reverse transcription PCR. The data showed that miR-30b-5p, miR-22-3p, and miR-378a-3p were significantly deregulated in AD, with area under the curve (AUC) of 0.668, 0.637, and 0.718, respectively. The combination of the three miRs gained a better diagnostic capability with AUC of 0.880. This finding revealed a miR panel as potential biomarker in the peripheral blood to distinguish AD from HC.

Noble-metal nanoparticle labelling multiplex miRNAs by ICP-MS readout with internal standard isotopes of 115In and 209Bi

Inductively coupled plasma-mass spectrometry (ICP-MS) is one of the most powerful techniques for multiplex nucleotide assay owing to the virtue of the high resolution of multiple-elements' mass to charge ratio, in a mass spectrum. Here, a small sized (less than 20 nm) noble-metal nanoparticle labelled ICP-MS (NP-ICP-MS) is proposed for high-throughput microRNA (miRNA) determination. Three miRNA targets - miR-486-5p, miR-221, and miR-21 - in serum, were distinguished by single-stranded DNA (ssDNA) probes labelled with a small sized noble-metal nanoparticle - silver nanoparticles (AgNPs), platinum nanoparticles (PtNPs), and gold nanoparticles (AuNPs). The counting isotopes ion intensity per second (CPS) of the noble-metal label versus internal standard isotope intensity of ¹¹⁵In and ²⁰⁹Bi, exhibited good linearity in the range 0.25 pM to 100 pM with correlation coefficients (R²) of 0.9680, 0.9305, and 0.9418. The specific sandwich-type miRNA assay using the sensitive NP-ICP-MS readout pushed the detection limits down to 0.18 pM for miR-221, 0.23 pM for miR-486-5p, and 0.22 pM for miR-21. And the relative standard deviations (RSDs) for 10 pM target miRNA were less than 3.7%. This work promises a potential ultrasensitive ICP-MS bioassay of multiplex miRNA biomarkers for clinical serum diagnosis.

Identification of an immune-based mRNA-lncRNA signature for overall survival in cervical squamous cell carcinoma

Aim: To better predict the survival of cervical squamous cell carcinoma (CESC) patients, we aimed to construct a signature according to different immune infiltration. Methods: We downloaded the RNA sequences of CESC patients from the Cancer Genome Atlas database. By using single-sample gene set enrichment analysis, we separated the samples into high- and low-immunity groups. Then we separated the samples into training and testing datasets and performed the following analyses: univariate, least absolute shrinkage and selection operator analysis, multivariate Cox regression analyses and weighted gene coexpression network analysis using R software. Gene ontology and Kyoto Encyclopedia of Genes and Genomes studies were performed using the Database for Annotation, Visualization and Integrated Discovery website. Results & conclusion: We finally identified a signature with three mRNAs and two lncRNAs: ADGRG5, HSH2D, ZMAT4, RBAKDN and LINC00200. In short, our study constructed an mRNA-lncRNA signature related to immune infiltration to better predict the survival of CESC patients.

Identification of mRNA and non-coding RNA hubs using network analysis in organ tropism regulated triple negative breast cancer metastasis

Triple negative breast cancer (TNBC) is aggressive in nature, resistant to conventional therapy and often ends in organ specific metastasis. In this study, publicly available datasets were used to identify miRNA, mRNA and lncRNA hubs. Using validated mRNA-miRNA, mRNA-mRNA and lncRNA-miRNA interaction information obtained from various databases, RNA interaction networks for TNBC and its subtype specific as well as organ tropism regulated metastasis were generated. Further, miRNA-mRNA-lncRNA triad classification was performed using social network analysis from subnetworks and visualized using Cytoscape. Survival analysis of the RNA hubs, oncoprint analysis for mRNAs and pathway analysis of the lncRNAs were also performed. Results indicated that two lncRNAs (NEAT1 and CASC7) and four miRNAs (hsa-miR-106b-5p, hsa-miR-148a-3p, hsa-miR-25-3p and hsa-let-7i-5p) were common between hubs identified in TNBC and TNBC associated metastasis. The exclusive hubs for TNBC associated metastasis were hsa-miR-200b-3p, SP1, HSPA4 and RAB1B. HMGA1 was the top ranked hub in mesenchymal subtype associated lung metastasis, while hsa-miR-27a-3p was identified as the top ranked hub mRNA in luminal androgen receptor subtype associated bone metastasis. When lncRNA associated pathway analysis was performed, Hs Cytoplasmic Ribosomal Protein pathway was found to be the most significant and among the selected hubs, CTNND1, SON and hsa-miR-29c emerged as TNBC survival markers. TP53, FOXA1, MTDH and HDGF were found as the top ranked mRNAs in oncoprint analysis. The pipeline proposed for the first time in this study with validated RNA interaction data integration and graph-based learning for miRNA-mRNA-lncRNA triad classification from RNA hubs may aid experimental cost reduction and its successful execution will allow it to be extended to other diseases too.

Integrated analysis of key microRNAs /TFs /mRNAs/ in HPV-positive cervical cancer based on microRNA sequencing and bioinformatics analysis

BACKGROUND

Cervical squamous cell carcinoma (CESC) is one of the most common malignancies associated with mortality in females. Its onset and prognosis are primarily concerned with persistent infection with high-risk types of human papillomavirus (HPV). However, the molecular mechanisms of HPV-positive CESC remain unclear.

METHODS

In this study, we conducted a high-throughput sequencing to identify differentially expressed miRNAs (DEMs). Besides, three series were selected from the Gene Expression Omnibus (GEO) database to identify differentially expressed genes (DEGs). Then the miRNA-TF-gene regulatory network was constructed using bioinformatic methods. Genes in the network were performed functional enrichment analysis and protein-protein interaction (PPI) network analysis. Ultimately, the expression levels of six key miRNAs, TFs, and mRNAs were validated by 20 HPV-positive CESC tissues and 15 normal cervical samples.

RESULTS

A total of 52 DEMs and 300 DEGs differed between the HPV-positive CESC and normal cervical samples. Then the miRNA-TF-gene regulatory network was constructed consisting of 22 miRNAs, 6 TFs, and 76 corresponding genes, among which miR-149-5p, miRNA-1248 and E2F4 acted as key regulators. PPI network analysis showed that ten genes including TOP2A, AURKA, CHEK1, KIF11, MCM4, MKI67, DTL, FOXM1, SMC4, and FBXO5 were recognized as hub genes with the highest connectivity degrees. Besides, five key molecules miRNA-149-5p, E2F4, KIF11, DTL, and SMC4 were suggested to play crucial roles in the development of HPV-positive CESC.

CONCLUSION

These results present a unique insight into the pathological mechanisms of HPV-positive CESC and possibly provides potential therapeutic targets.

Identification of a Recurrence Signature and Validation of Cell Infiltration Level of Thyroid Cancer Microenvironment

Though many patients with thyroid cancer may be indolent, there are still about 50% lymph node metastases and 20% the recurrence rates. There is still no ideal method to predict its relapse. In this study, we analyzed the gene transcriptome profiles of eight Gene Expression Omnibus (GEO), and next screened 77 commonly differential expressed genes. Next, Least Absolute Shrinkage and Selection Operator (LASSO) regression model was performed and seven genes (i.e., FN1, PKIA, TMEM47, FXYD6, SDC2, CD44, and GGCT) were then identified, which is highly associated with recurrence data from the Cancer Genome Atlas (TCGA) database. These patients were then divided into low and high-risk groups with specific risk-score formula. Univariate and multivariate Cox regression further revealed that the 7-mRNA signature plays a functional causative role independent of clinicopathological characteristics. The 7-mRNA-signature integrated nomogram showed better discrimination, and decision curve analysis demonstrated that it is clinically useful. Besides, patient with lower risk score shows a relatively lower level of activated dendritic cells (DCs), resting DCs, regulatory T cells and γδT cells, and process of DCs apoptotic. In conclusion, our present immune-related classifier could produce a potential tool for predicting early-relapse.

HAND2-AS1 inhibits invasion and metastasis of cervical cancer cells via microRNA-330-5p-mediated LDOC1

Background

Cervical cancer is a serious disease with complicated pathogenesis and thus there is an urgent need to find novel targets for the treatment. Recently, long non-coding RNAs (lncRNAs) have emerged as critical factors in tumorigenesis. In this study, we aimed to investigate the mechanism of HAND2 antisense RNA 1 (HAND2-AS1) on the invasion and metastasis of cervical cancer cells.

Methods

The expression patterns of HAND2-AS1, microRNA-330-5p (miR-330-5p) and leucine zipper down-regulated in cancer 1 (LDOC1) in cervical cancer were characterized by RT-qPCR and western blot analysis. Dual luciferase reporter assay and RIP were applied to verify relationship between HAND2-AS1, miR-330-5p and LDOC1. Fluorescence in situ hybridization (FISH) was used to detect the subcellular localization of HAND2-AS1. Besides, viability, invasion and migration ability of HeLa cells were investigated by cell counting kit-8 (CCK-8) and Transwell assays respectively. Hematoxylin-eosin staining was performed for lymph node metastasis detection. In addition, the tumor growth in nude mice was evaluated.

Results

Low expression of HAND2-AS1 and LDOC1, and high expression of miR-330-5p were detected in cervical cancer tissues and cells. It was found that binding of HAND2-AS1 to miR-330-5p results in upregulation of LDOC1 expression. Also, overexpressed HAND2-AS1 and LDOC1 or down-regulated miR-330-5p inhibited expression of proliferation-associated proteins Ki-67, PCNA, migration-associated proteins N-cad and invasion-related proteins MMP-2, MMP-9 as well as lymph node metastasis. Moreover, HAND2-AS1 inhibited tumor formation and lymph node metastasis by binding to miR-330-5p in vivo.

Conclusion

HAND2-AS1 promotes LDOC1 expression by competitively binding to miR-330-5p and consequently inhibiting cervical cancer cell invasion and metastasis. This could facilitate development of therapeutic strategies against cervical cancer.