Identification of potential micro-messenger RNAs (miRNA-mRNA) interaction network of osteosarcoma

A comprehensive overview of liquid biopsy applications in pediatric solid tumors

Liquid biopsies are emerging as an alternative source for pediatric cancer biomarkers with potential applications during all stages of patient care, from diagnosis to long-term follow-up. While developments within this field are reported, these mainly focus on dedicated items such as a specific liquid biopsy matrix, analyte, and/or single tumor type. To the best of our knowledge, a comprehensive overview is lacking. Here, we review the current state of liquid biopsy research for the most common non-central nervous system pediatric solid tumors. These include neuroblastoma, renal tumors, germ cell tumors, osteosarcoma, Ewing sarcoma, rhabdomyosarcoma and other soft tissue sarcomas, and liver tumors. Within this selection, we discuss the most important or recent studies involving liquid biopsy-based biomarkers, anticipated clinical applications, and the current challenges for success. Furthermore, we provide an overview of liquid biopsy-based biomarker publication output for each tumor type based on a comprehensive literature search between 1989 and 2023. Per study identified, we list the relevant liquid biopsy-based biomarkers, matrices (e.g., peripheral blood, bone marrow, or cerebrospinal fluid), analytes (e.g., circulating cell-free and tumor DNA, microRNAs, and circulating tumor cells), methods (e.g., digital droplet PCR and next-generation sequencing), the involved pediatric patient cohort, and proposed applications. As such, we identified 344 unique publications. Taken together, while the liquid biopsy field in pediatric oncology is still behind adult oncology, potentially relevant publications have increased over the last decade. Importantly, steps towards clinical implementation are rapidly gaining ground, notably through validation of liquid biopsy-based biomarkers in pediatric clinical trials.

MicroRNA-98: the multifaceted regulator in human cancer progression and therapy

MicroRNA-98 (miR-98) stands as an important molecule in the intricate landscape of oncology. As a subset of microRNAs, these small non-coding RNAs have accompanied a new era in cancer research, underpinning their significant roles in tumorigenesis, metastasis, and therapeutic interventions. This review provides a comprehensive insight into the biogenesis, molecular properties, and physiological undertakings of miR-98, highlighting its double-edged role in cancer progression-acting both as a tumor promoter and suppressor. Intriguingly, miR-98 has profound implications for various aspects of cancer progression, modulating key cellular functions, including proliferation, apoptosis, and the cell cycle. Given its expression patterns, the potential of miR-98 as a diagnostic and prognostic biomarker, especially in liquid biopsies and tumor tissues, is explored, emphasizing the hurdles in translating these findings clinically. The review concludes by evaluating therapeutic avenues to modulate miR-98 expression, addressing the challenges in therapy resistance, and assessing the efficacy of miR-98 interventions. In conclusion, while miR-98's involvement in cancer showcases promising diagnostic and therapeutic avenues, future research should pivot towards understanding its role in tumor-stroma interactions, immune modulation, and metabolic regulation, thereby unlocking novel strategies for cancer management.

Hsa-miR-301a-3p inhibited the killing effect of natural killer cells on non-small cell lung cancer cells by regulating RUNX3

BACKGROUND

Non-small cell lung cancer (NSCLC) is the most commonly diagnosed solid tumor. Natural killer (NK) cell-based immunotherapy is a promising anti-tumor strategy in various cancers including NSCLC.

OBJECTIVE

We aimed to investigate the specific mechanisms that regulate the killing effect of NK cells to NSCLC cells.

METHODS

Reverse transcription-quantitative PCR (RT-qPCR) assay was applied to measure the levels of hsa-microRNA (miR)-301a-3p and Runt-related transcription factor 3 (RUNX3). Enzyme-linked immunosorbent assay (ELISA) was used to measure the levels of IFN-γ and TNF-α. Lactate dehydrogenase assay was applied to detect the killing effect of NK cells. Dualluciferase reporter assay and RNA immunoprecipitation (RIP) assay were carried out to confirm the regulatory relationship between hsa-miR-301a-3p and RUNX3.

RESULTS

A low expression of hsa-miR-301a-3p was observed in NK cells stimulated by IL-2. The levels of IFN-γ and TNF-α were increased in NK cells of the IL-2 group. Overexpression of hsa-miR-301a-3p reduced the levels of IFN-γ and TNF-α as well as the killing effect of NK cells. Furthermore, RUNX3 was identified to be a target of hsamiR-301a-3p. hsa-miR-301a-3p suppressed the cytotoxicity of NK cells to NSCLC cells by inhibiting the expression of RUNX3. We found hsa-miR-301a-3p promoted tumor growth by suppressing the killing effect of NK cells against NSCLC cells in vivo.

CONCLUSIONS

Hsa-miR-301a-3p suppressed the killing effect of NK cells on NSCLC cells by targeting RUNX3, which may provide promising strategies for NK cell-based antitumor therapies.

Integration of Human and Viral miRNAs in Epstein-Barr Virus-Associated Tumors and Implications for Drug Repurposing

Epstein-Barr virus (EBV) is associated with several tumors, and has substantial relevance for public health. Therapeutics innovation for EBV-related disorders is much needed. In this context, miRNAs are noncoding RNA molecules that play vital roles in EBV infection. miRNA-Seq and RNA-Seq data for EBV-associated clinical samples and cell lines have been generated, but their detailed integrative analyses, and exploitation for drug repurposing against EBV are lacking. Hence, we identified and analyzed the differentially expressed miRNAs (DEmiRs) in EBV-infected cell lines (28) and infected (28) and uninfected human tissue (20) samples using an in-house pipeline. We found significantly enriched host miRNAs like hsa-mir-3651, hsa-mir-1248, and hsa-mir-29c-3p in EBV-infected samples from EBV-associated nasopharyngeal carcinoma and Hodgkin's lymphoma, among others. Furthermore, we also identified significantly enriched novel miRNAs such as hsa-mir-29c-3p, hsa-mir-3651, and hsa-mir-98-3p, which were not previously reported in EBV-related tumors. Differentially expressed mRNAs (DEMs) were identified in EBV-infected cell lines (21) and uninfected human tissue (14) samples. We predicted and selected 1572 DEMs (upregulated) that are targeted by 547 DEmiRs (downregulated). These were further classified into essential (870) and nonessential (702) genes. Moreover, a miRNA-mRNA network was developed for the hub miRNAs. Importantly, we used the DEMs during EBV latent infection types I, II, and III to identify the candidate drugs for repurposing: Glyburide, Levodopa, Nateglinide, and Stiripentol, among others. To the best of our knowledge, this is the first integrative analyses that identified DEmiRs and DEMs as potential therapeutic targets and predicted drugs as potential candidates for repurposing against EBV-related tumors.

Construction of Osteosarcoma Diagnosis Model by Random Forest and Artificial Neural Network

Osteosarcoma accounts for 28% of primary bone malignancies in adults and up to 56% in children and adolescents (<20 years). However, early diagnosis and treatment are still inadequate, and new improvements are still needed. Missed diagnoses exist due to fewer traditional diagnostic methods, and clinical symptoms are often already present before diagnosis. This study aimed to develop novel and efficient predictive models for the diagnosis of osteosarcoma and to identify potential targets for exploring osteosarcoma markers. First, osteosarcoma and normal tissue expression microarray datasets were downloaded from the Gene Expression Omnibus (GEO). Then we screened the differentially expressed genes (DEGs) in the osteosarcoma and normal groups in the training group. Next, in order to explore the biologically relevant role of DEGs, Metascape and enrichment analyses were also performed on DEGs. The “randomForest” and “neuralnet” packages in R software were used to select representative genes and construct diagnostic models for osteosarcoma. The next step is to validate the model of the artificial neural network. Then, we performed an immune infiltration analysis by using the training set data. Finally, we constructed a prognostic model using representative genes for prognostic analysis. The copy number of osteosarcoma was also analyzed. A random forest classifier identified nine representative genes (ANK1, TGFBR3, RSF21, HSPB8, ITGA7, RHD, AASS, GREM2, NFASC). HSPB8, RHD, AASS, and NFASC were genes we identified that have not been previously reported to be associated with osteosarcoma. The osteosarcoma diagnostic model we constructed has good performance with areas under the curves (AUCs) of 1 and 0.987 in the training and validation groups, respectively. This study opens new horizons for the early diagnosis of osteosarcoma and provides representative markers for the future treatment of osteosarcoma. This is the first study to pioneer the establishment of a genetic diagnosis model for osteosarcoma and advance the development of osteosarcoma diagnosis and treatment.

Screening for MicroRNA combination with engineered exosomes as a new tool against osteosarcoma in elderly patients

The most common primary malignant bone sarcoma is Osteogenic sarcoma (OS) which has a bimodal age distribution. Unfortunately, the treatment of OS was less effective for elderly patients than for younger ones. The study aimed to explore a new microRNA (miRNA) which can bind to combining engineered exosomes for treatment of older OS patients. Based on GSE65071 and miRNet 2.0, two up-regulated miRNAs (miR-328, miR-107) and seven down-regulated miRNAs (miR-133b, miR-206, miR-1-3p, miR-133a, miR-449a, miR-181daysay, miR-134) were selected. Next, we used FunRich software to predict the up-stream transcription factors (TFs) of differentially expressed miRNAs (DE-miRNAs). By comparing target genes predicted from DE-miRNAs with differentially expressed genes, we identified 12 down-regulated and 310 up-regulated mRNAs. For KEGG analysis, the most enriched KEGG pathway was Cell cycle, Spliceosome, and Protein digestion and absorption. By using protein-protein interactions network, topological analysis algorithm and GEPIA database, miR-449a /CCNB1 axis was identified. Experiments in vitro were conducted to confirm the results too. MiRNA-449a is down-regulated in osteosarcoma and suppresses cell proliferation by targeting CCNB1. Our findings not only reveal a novel mechanism of miR-449a /CCNB1 in OS but also had laid the groundwork for further investigation and analysis in the field of exosome engineering.

Bioinformatics analysis of lncRNAs in the occurrence and development of osteosarcoma

BACKGROUND

Osteosarcoma (OS) is a disease with high mortality in children and adolescents, and metastasis is one of its important clinical features. However, the molecular mechanism of OS occurrence is not completely clear. Thus, we screened potential biomarkers of OS and analyze their prognostic value.

METHODS

The Cancer Genome Atlas (TCGA) datasets were used to analyze the differential lncRNAs in patients with OS of different immune score and the lncRNAs expressed by immune cells. Cox regression was used to develop the prognosis prediction model and specify the prognosis outcomes. Risk-proportional regression model was constructed, and the samples were divided into high and low groups based on the risk scores for the survival analysis. The areas under the receiver operating characteristic (ROC) curve were calculated and the risk-score model was verified. Finally, using 4 gene sets (comprising chemokines, immune checkpoint blockades, immune activity-related genes, and immune cells), and 4 analysis tools (CIBERSORT, TIMER, XCELL and MCP) to evaluated tumor immune infiltration.

RESULTS

Twenty-nine long non-coding ribonucleic acids (lncRNAs) were obtained from the intersection of the screened lncRNAs. Caspase recruitment domain-containing protein 8-antisense RNA 1 (CARD8-AS1), lncRNA five prime to Xist (FTX), KAT8 regulatory NSL complex unit 1-antisense RNA 1 (KANSL1-AS1), Neuroplastin Intronic Transcript 1 (NPTN-IT1), oligodendrocyte maturation-associated long intervening non-coding RNA (OLMALINC) and RPARP Antisense RNA 1 (RPARP-AS1) were found to be correlated with survival. Univariate and multivariate regression analysis showed risk score [HR (hazard ratio) 3.5, P value 0.0043; HR 3.7, P value 0.0033] and metastasis (HR 4.7, P value 6.60E-05; HR 4.8, P value 8.36E-05) were the key factors of patients with OS. The areas under curves (AUCs) of the 1-, 3-, and 5-year ROC curves of the prognostic model were 0.715, 0.729, and 0.771. The low-risk patients tended to have a high abundance of immune cells.

CONCLUSIONS

This study showed that a risk score based on 6 lncRNAs has potential value in the prognosis of OS, and patients with low-risk scores have high immune cell infiltration and good prognosis. This study may enrich understandings of underlying mechanisms related to the occurrence and development of OS.

Necdin, one of the important pathway proteins in the regulation of osteosarcoma progression by microRNA-200c

MicroRNA-200c (miR-200c) generally acts as a tumor suppressor in multiple cancer types and a promising therapeutic target in tumorigenesis. However, only a few studies have explained the role of miR-200c in the development of osteosarcoma (OS). In this study, we investigated the role of miR-200c in OS progression and identified the regulatory pathway protein NDN involved in inhibiting the occurrence and development of OS. Firstly, we found that miR-200c is downregulated in OS cells and tissues. As well, in vitro and in vivo experiments showed that upregulating miR-200c inhibits the proliferation, invasion, metastasis of Saos-2 cells, promotes the apoptosis of Saos-2 cells and suppresses tumor growth in mice, indicating miR-200c plays a major role in regulating the OS progression. Furthermore, bioinformatics analysis showed that an anti-tumor protein, necdin (NDN), might be a potential target by miR-200c. To verify this hypothesis, we measured the expression level of NDN in OS cells and tissues and found NDN is downregulated, suggesting NDN is functional in OS progression. Moreover, we found that the expression levels of NDN and miR-200c in in vivo and in vitro experiments were positively correlated. However, the results of dual-luciferase reporter gene experiment showed miR-200c does not directly act on the 3ʹ untranslated region (UTR) of NDN gene, indicating that NDN might be an important pathway protein which regulates OS progression in the presence of miR-200c. Therefore, miR-200c/NDN could be potential targets for developing effective treatment against OS.

The effect and mechanism of miR-30e-5p targeting SNAI1 to regulate epithelial-mesenchymal transition on pancreatic cancer

Accumulating evidence indicates that abnormally expressed microRNAs (miRNAs, miRs) contribute to cancer progression. Nonetheless, the role of miR-30e-5p in pancreatic cancer (PCa) remains unclear. In this study, using quantitative real-time polymerase chain reaction analysis, we found that miR-30e-5p expression was downregulated in human PCa tissues compared with that in normal para-cancerous tissues. After transfecting with miR-30e-5p inhibitors, miR-30e-5p mimics, or empty vectors in the BxPC-3 and PANC-1 cells, respectively, the experiments revealed that the upregulation of miR-30e-5p expression inhibited cell growth, invasion, migration and epithelial-mesenchymal transition (EMT), and promoted apoptosis, while miR-30e-5p downregulation had the opposite effects. RNA sequencing of miR-30e-5p inhibitor-, miR-30e-5p mimic-, and the negative control (NC)-treated groups revealed that miR-30e-5p may affect epithelial cell differentiation, cell growth and death. Next, the snail family transcriptional repressor 1 (SNAI1) was predicted and verified as the target gene of miR-30e-5p using bioinformatics analysis and luciferase assays. SNAI1 expression levels were decreased in the PCa cells transfected with miR-30e-5p mimics, whereas the opposite was observed in the cells transfected with miR-30e-5p inhibitors. Subsequently, PCa cells were transfected with a vector overexpressing SNAI1 (OE-SNAI1) and miR-30e-5p mimics, miR-30e-5p inhibitors, or empty vectors. Compared with that in the OE-SNAI1 + miR-30e-5p NC group, transfection with OE-SNAI1 + miR-30e-5p mimics inhibited the PCa cell growth, migration, and increased apoptosis, whereas transfection with OE-SNAI1 + miR-30e-5p inhibitors had the opposite effects. In conclusion, miR-30e-5p potentially inhibits PCa cell proliferation, migration, and invasion via the SNAI1/EMT axis.

β-Amyrin ameliorates diabetic nephropathy in mice and regulates the miR-181b-5p/HMGB2 axis in high glucose-stimulated HK-2 cells

Diabetic nephropathy (DN) is a diabetic complication that can cause renal failure. β-amyrin has been identified to possess anti-diabetic property. This study was designed to evaluate the potential role of β-amyrin in DN and its underlying mechanism. Streptozotocin-induced diabetic mice were used as the in vivo model, and high glucose (HG)-stimulated human proximal tubular HK-2 cells were utilized as the in vitro model. Renal histological changes in mice were assessed by hematoxylin-eosin and periodic acid-Schiff staining. HK-2 cell viability and apoptosis were detected by Cell Counting Kit-8 assay and flow cytometry analysis, respectively. β-amyrin was found to ameliorate kidney injury in DN mice and suppressed inflammatory response as well as apoptosis of HG-stimulated HK-2 cells. miR-181-5p expression in murine renal tissues and HK-2 cells was detected by in situ hybridization (ISH) and fluorescence in situ hybridization (FISH). MiR-181b-5p, a previously identified target for diabetic kidney disease, was downregulated in renal tissues and HG stimulated HK-2 cells, and β-amyrin induced the upregulation of miR-181b-5p. Binding relationship between miR-181b-5p and high mobility group box 2 (HMGB2) was confirmed by luciferase reporter assay. MiR-181b-5p bound to 3' untranslated region of HMGB2 to suppress its expression. As shown by immunohistochemical staining and immunofluorescence staining, HMGB2 was upregulated in the in vivo and in vitro models of DN, and β-amyrin induced the downregulation of HMGB2. Moreover, HMGB2 overexpression neutralized the suppressive effects of miR-181b-5p elevation on the inflammatory response and apoptosis of HG-treated HK-2 cells. Overall, β-amyrin ameliorates DN in mice and suppresses inflammatory response and apoptosis of HG-stimulated HK-2 cells via the miR-181b-5p/HMGB2 axis.

Identification and preliminary validation of a four-gene signature to predict metastasis and survival in osteosarcoma

Osteosarcoma is a primary malignant bone tumor that occurs frequently in children and adolescents and has a propensity for drug resistance, recurrence, and metastasis. The purpose of this study was to identify potential target genes to predict metastasis and survival in patients with osteosarcoma. We analyzed gene expression profiles and corresponding clinical data of patients with osteosarcoma in the Gene Expression Omnibus database and identified 202 genes that were differentially expressed between osteosarcoma cells and normal osteoblasts. Univariate and multivariable Cox regression analyses identified four risk genes that affected osteosarcoma prognosis: MCAM, ENPEP, LRRC1, and CPE. Independent prognostic analyses and clinical correlation studies showed that the four risk genes constituted an independent prognostic signature that correlated with survival and clinical parameters including age and distant metastasis. In a single-sample Gene Set Enrichment Analysis, risk scores based on the prognostic signature correlated with tumor infiltration by immune cells and immune functions in osteosarcoma. A subsequent analysis showed that the expression levels of the four genes in the prognostic signature were predictive of overall survival and metastasis-free survival of patients with osteosarcoma. Furthermore, Human Cancer Metastasis Database and qRT-PCR analyses demonstrated that the four risk genes are overexpressed in osteosarcoma tissues and cell lines. In summary, we developed and validated a four-gene prognostic signature that may be useful in osteosarcoma diagnosis and metastasis prediction.