Integrated bioinformatics analysis of miRNA expression in Ewing sarcoma and potential regulatory effects of miR-21 via targeting ALCAM/CD166

mRNA expression profile reveals differentially expressed genes in splenocytes of experimental autoimmune encephalomyelitis model

Experimental autoimmune encephalomyelitis (EAE) is a mouse model that can be used to investigate aetiology, pathogenesis, and treatment approaches for multiple sclerosis (MS). A novel integrated bioinformatics approach was used to understand the involvement of differentially expressed genes (DEGs) in the spleen of EAE mice through data mining of existing microarray and RNA-seq datasets. We screened differentially expressed mRNAs using mRNA expression profile data of EAE spleens taken from Gene Expression Omnibus (GEO). Functional and pathway enrichment analyses of DEGs were performed by Database for Annotation, Visualization, and Integrated Discovery (DAVID). Subsequently, the DEGs-encoded protein-protein interaction (PPI) network was constructed. The 784 DEGs in GSE99300 A.SW PP-EAE mice spleen mRNA profiles, 859 DEGs in GSE151701 EAE mice spleen mRNA profiles, and 646 DEGs in GSE99300 SJL/J PP-EAE mice spleen mRNA profiles were explored. Functional enrichment of 55 common DEGs among 3 sub-datasets revealed several immune-related terms, such as neutrophil extravasation, leucocyte migration, antimicrobial humoral immune response mediated by an antimicrobial peptide, toll-like receptor 4 bindings, IL-17 signalling pathway, and TGF-beta signalling pathway. In the screening of 10 hub genes, including MPO, ELANE, CTSG, LTF, LCN2, SELP, CAMP, S100A9, ITGA2B, and PRTN3, and in choosing and validating the 5 DEGs, including ANK1, MBOAT2, SLC25A21, SLC43A1, and SOX6, the results showed that SLC43A1 and SOX6 were significantly decreased in EAE mice spleen. Thus this study offers a list of genes expressed in the spleen that might play a key role in the pathogenesis of EAE.

Oncomeric Profiles of microRNAs as New Therapeutic Targets for Treatment of Ewing's Sarcoma: A Composite Review

Ewing's sarcoma is a rare type of cancer that forms in bones and soft tissues in the body, affecting mostly children and young adults. Current treatments for ES are limited to chemotherapy and/or radiation, followed by surgery. Recently, microRNAs have shown favourable results as latent diagnostic and prognostic biomarkers in various cancers. Furthermore, microRNAs have shown to be a good therapeutic agent due to their involvement in the dysregulation of various molecular pathways linked to tumour progression, invasion, angiogenesis, and metastasis. In this review, comprehensive data mining was employed to explore various microRNAs that might have therapeutic potential as target molecules in the treatment of ES.

MicroRNA Content of Ewing Sarcoma Derived Extracellular Vesicles Leads to Biomarker Potential and Identification of a Previously Undocumented EWS-FLI1 Translocation

Objective: Ewing Sarcoma Family of Tumors (ESFT) are a highly aggressive pediatric bone and soft tissue malignancy with poor outcomes in the refractory and recurrent setting. Over 90% of Ewing Sarcoma (ES) tumors are driven by the pathognomonic EWS-ETS chimeric transcripts and their corresponding oncoproteins. It has been suggested that the EWS-ETS oncogenic action can mediate microRNA (miRNA) processing. Importantly, small extracellular vesicles (sEVs), including those frequently referred to as exosomes have been shown to be highly enriched with tumor-derived small RNAs such as miRNAs. We hypothesized that ESFT-specific sEVs are enriched with certain miRNAs which could be utilized toward an exo-miRNA biomarker signature specific to this disease. Methods: We performed miRNAseq to compare both the exo-derived and cell-derived miRNA content from 8 ESFT, 2 osteosarcoma, 2 non-cancerous cell lines, and pediatric plasma samples. Results: We found that sEVs derived from ESFT cells contained nearly 2-fold more number of unique individual miRNAs as compared to non-ESFT samples. Quantitative analysis of the differential enrichment of sEV miRNAs resulted in the identification of 62 sEV-miRNAs (exo-miRNAs) with significant (P < .05) enrichment variation between ESFT and non-ESFT sEV samples. To determine if we could utilize this miRNA signature to diagnose ESFT patients via a liquid biopsy, we analyzed the RNA content of total circulating sEVs isolated from 500 µL plasma from 5 pediatric ESFT patients, 2 pediatric osteosarcoma patients, 2 pediatric rhabdomyosarcoma patients, and 4 non-cancer pediatric controls. Pearson's clustering of 60 of the 62 candidate exo-miRNAs correctly identified 80% (4 of 5) of pathology confirmed ESFT patients. Importantly, RNAseq analysis of tumor tissue from the 1 outlier, revealed a previously uncharacterized EWS-FLI1 translocation.Conclusions: Taken together, these findings support the development and validation of an exo-miRNA-based liquid biopsy to aid in the diagnosis and monitoring of ESFT.

An Integrated Bioinformatics Analysis of the Potential Regulatory Effects of miR-21 on T-cell Related Target Genes in Multiple Sclerosis

Background:

Overexpression of miR-21 is a characteristic feature of patients with Multiple Sclerosis (MS) and is involved in gene regulation and the expression enhancement of pro-inflammatory factors including IFNγ and TNF-α following stimulation of T-cells via the T Cell Receptor (TCR). In this study, a novel integrated bioinformatics analysis was used to obtain a better understanding of the involvement of miR-21 in the development of MS, its protein biomarker signatures, RNA levels, and drug interactions through existing microarray and RNA-seq datasets of MS.

Methods:

In order to obtain data on the Differentially Expressed Genes (DEGs) in patients with MS and normal controls, the GEO2R web tool was used to analyze the Gene Expression Omnibus (GEO) datasets, and then Protein-Protein Interaction (PPI) networks of co-expressed DEGs were designed using STRING. A molecular network of miRNA-genes and drugs based on differentially expressed genes was created for T-cells of MS patients to identify the targets of miR-21, that may act as important regulators and potential biomarkers for early diagnosis, prognosis and, potential therapeutic targets for MS.

Results:

It found that seven genes (NRIP1, ARNT, KDM7A, S100A10, AK2, TGFβR2, and IL-6R) are regulated by drugs used in MS and miR-21. Finally, three overlapping genes (S100A10, NRIP1, KDM7A) were identified between miRNA-gene-drug network and nineteen genes as hub genes which can reflect the pathophysiology of MS.

Conclusion:

Our findings suggest that miR-21 and MS-related drugs can act synergistically to regulate several genes in the existing datasets, and miR-21 inhibitors have the potential to be used in MS treatment.

Identification of CDK2-Related Immune Forecast Model and ceRNA in Lung Adenocarcinoma, a Pan-Cancer Analysis

Background

Tumor microenvironment (TME) plays important roles in different cancers. Our study aimed to identify molecules with significant prognostic values and construct a relevant Nomogram, immune model, competing endogenous RNA (ceRNA) in lung adenocarcinoma (LUAD).

Methods

“GEO2R,” “limma” R packages were used to identify all differentially expressed mRNAs from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. Genes with P-value <0.01, LogFC>2 or <-2 were included for further analyses. The function analysis of 250 overlapping mRNAs was shown by DAVID and Metascape software. By UALCAN, Oncomine and R packages, we explored the expression levels, survival analyses of CDK2 in 33 cancers. “Survival,” “survminer,” “rms” R packages were used to construct a Nomogram model of age, gender, stage, T, M, N. Univariate and multivariate Cox regression were used to establish prognosis-related immune forecast model in LUAD. CeRNA network was constructed by various online databases. The Genomics of Drug Sensitivity in Cancer (GDSC) database was used to explore correlations between CDK2 expression and IC50 of anti-tumor drugs.

Results

A total of 250 differentially expressed genes (DEGs) were identified to participate in many cancer-related pathways, such as activation of immune response, cell adhesion, migration, P13K-AKT signaling pathway. The target molecule CDK2 had prognostic value for the survival of patients in LUAD (P = 5.8e-15). Through Oncomine, TIMER, UALCAN, PrognoScan databases, the expression level of CDK2 in LUAD was higher than normal tissues. Pan-cancer analysis revealed that the expression, stage and survival of CDK2 in 33 cancers, which were statistically significant. Through TISIDB database, we selected 13 immunodepressants, 21 immunostimulants associated with CDK2 and explored 48 genes related to these 34 immunomodulators in cBioProtal database (P < 0.05). Gene Set Enrichment Analysis (GSEA) and Metascape indicated that 49 mRNAs were involved in PUJANA ATM PCC NETWORK (ES = 0.557, P = 0, FDR = 0), SIGNAL TRANSDUCTION (ES = –0.459, P = 0, FDR = 0), immune system process, cell proliferation. Forest map and Nomogram model showed the prognosis of patients with LUAD (Log-Rank = 1.399e-08, Concordance Index = 0.7). Cox regression showed that four mRNAs (SIT1, SNAI3, ASB2, and CDK2) were used to construct the forecast model to predict the prognosis of patients (P < 0.05). LUAD patients were divided into two different risk groups (low and high) had a statistical significance (P = 6.223e-04). By “survival ROC” R package, the total risk score of this prognostic model was AUC = 0.729 (SIT1 = 0.484, SNAI3 = 0.485, ASB2 = 0.267, CDK2 = 0.579). CytoHubba selected ceRNA mechanism medicated by potential biomarkers, 6 lncRNAs-7miRNAs-CDK2. The expression of CDK2 was associated with IC50 of 89 antitumor drugs, and we showed the top 20 drugs with P < 0.05.

Conclusion

In conclusion, our study identified CDK2 related immune forecast model, Nomogram model, forest map, ceRNA network, IC50 of anti-tumor drugs, to predict the prognosis and guide targeted therapy for LUAD patients.

microRNAs and bone tumours: Role of tiny molecules in the development and progression of chondrosarcoma, of giant cell tumour of bone and of Ewing's sarcoma

The increasing interest on microRNAs (miRNAs), small non-coding RNA molecules containing about 22 nucleotides, about their biological functions led researchers to discover that they are actively involved in several biological processes. In the last decades, miRNAs become one of the most topic of cancer research. miRNAs, thanks to their function, are the perfect molecules to modulate multiple signaling pathways and gene expression in cancer, with the consequent capacity to modulate cancerous processes, such as cellular proliferation, invasion, metastasis and chemoresistance in various tumours. In the last years, several studies have demonstrated the role of miRNAs in their pathophysiology, but little we know about the underlying mechanism that lead to bone tumours like chondrosarcoma (COS), giant cell tumour of bone (GCTB) and Ewing sarcoma (EWS) to still be highly aggressive and resistant tumours. An exploration of the role of miRNAs in the biology of them will permit to researchers to find new molecular mechanisms that can be used to develop new and more effective therapies against these bone tumours. Here we present a comprehensive study of the latest discoveries which have been performed in relation to the role of miRNAs in the neoplastic processes which characterize COS, EWS and GCTB, demonstrating how these tiny molecules can act as tumour promoters or as tumour suppressors and how they can be used for improving therapeutic approaches.

The Landscape of Regulatory Noncoding RNAs in Ewing's Sarcoma

Ewing’s sarcoma (ES) is a pediatric sarcoma caused by a chromosomal translocation. Unlike in most cancers, the genomes of ES patients are very stable. The translocation product of the EWS-FLI1 fusion is most often the predominant genetic driver of oncogenesis, and it is pertinent to explore the role of epigenetic alterations in the onset and progression of ES. Several types of noncoding RNAs, primarily microRNAs and long noncoding RNAs, are key epigenetic regulators that have been shown to play critical roles in various cancers. The functions of these epigenetic regulators are just beginning to be appreciated in ES. Here, we performed a comprehensive literature review to identify these noncoding RNAs. We identified clinically relevant tumor suppressor microRNAs, tumor promoter microRNAs and long noncoding RNAs. We then explored the known interplay between different classes of noncoding RNAs and described the currently unmet need for expanding the noncoding RNA repertoire of ES. We concluded the review with a discussion of epigenetic regulation of ES via regulatory noncoding RNAs. These noncoding RNAs provide new avenues of exploration to develop better therapeutics and identify novel biomarkers.

The role of microRNAs in the osteogenic and chondrogenic differentiation of mesenchymal stem cells and bone pathologies

Mesenchymal stem cells (MSCs) have been identified in many adult tissues. MSCs can regenerate through cell division or differentiate into adipocytes, osteoblasts and chondrocytes. As a result, MSCs have become an important source of cells in tissue engineering and regenerative medicine for bone tissue and cartilage. Several epigenetic factors are believed to play a role in MSCs differentiation. Among these, microRNA (miRNA) regulation is involved in the fine modulation of gene expression during osteogenic/chondrogenic differentiation. It has been reported that miRNAs are involved in bone homeostasis by modulating osteoblast gene expression. In addition, countless evidence has demonstrated that miRNAs dysregulation is involved in the development of osteoporosis and bone fractures. The deregulation of miRNAs expression has also been associated with several malignancies including bone cancer. In this context, bone-associated circulating miRNAs may be useful biomarkers for determining the predisposition, onset and development of osteoporosis, as well as in clinical applications to improve the diagnosis, follow-up and treatment of cancer and metastases. Overall, this review will provide an overview of how miRNAs activities participate in osteogenic/chondrogenic differentiation, while addressing the role of miRNA regulatory effects on target genes. Finally, the role of miRNAs in pathologies and therapies will be presented.

Bioinformatics Analysis Reveals MicroRNA-193a-3p Regulates ACTG2 to Control Phenotype Switch in Human Vascular Smooth Muscle Cells

Aortic dissection (AD) is among the most fatal cardiovascular diseases. However, the pathogenesis of AD remains poorly understood. This study aims to integrate the microRNAs (miRNA) and mRNA profiles and use bioinformatics analyses with techniques in molecular biology to delineate the potential mechanisms involved in the development of AD. We used the human miRNA and mRNA microarray datasets GSE98770, GSE52093, and GEO2R, Venn diagram analysis, gene ontology, and protein–protein interaction networks to identify target miRNAs and mRNAs involved in AD. RNA interference, western blotting, and luciferase reporter assays were performed to validate the candidate miRNAs and mRNAs in AD tissues and human vascular smooth muscle cells (VSMCs). Furthermore, we studied vascular smooth muscle contraction in AD. In silico analyses revealed that miR-193a-3p and ACTG2 were key players in the pathogenesis of AD. miR-193a-3p was upregulated in the AD tissues. We also found that biomarkers for the contractile phenotype in VSMCs were downregulated in AD tissues. Overexpression and depletion of miR-193a-3p enhanced and suppressed VSMC proliferation and migration, respectively. Dual luciferase reporter assays confirmed that ACTG2 was a target of miR-193a-3p. ACTG2 was also downregulated in human AD tissues and VMSCs overexpressing miR-193a-3p. Taken together, miR-193a-3p may be a novel regulator of phenotypic switching in VSMCs and the miR-193a-3p/ACTG2 axis may serve as a promising diagnostic biomarker and therapeutic candidate for AD.

Systems Biology Analysis for Ewing Sarcoma

Ewing sarcoma (EwS) is a highly aggressive pediatric bone cancer that is defined by a somatic fusion between the EWSR1 gene and an ETS family member, most frequently the FLI1 gene, leading to expression of a chimeric transcription factor EWSR1-FLI1. Otherwise, EwS is one of the most genetically stable cancers. The situation when the major cancer driver is well known looks like a unique opportunity for applying the systems biology approach in order to understand the EwS mechanisms as well as to uncover some general mechanistic principles of carcinogenesis. A number of studies have been performed revealing the direct and indirect effects of EWSR1-FLI1 on multiple aspects of cellular life. Nevertheless, the emerging picture of the oncogene action appears to be highly complex and systemic, with multiple reciprocal influences between the immediate consequences of the driver mutation and intracellular and intercellular molecular mechanisms, including regulation of transcription, epigenome, and tumoral microenvironment. In this chapter, we present an overview of existing molecular profiling resources available for EwS tumors and cell lines and provide an online comprehensive catalogue of publicly available omics and other datasets. We further highlight the systems biology studies of EwS, involving mathematical modeling of networks and integration of molecular data. We conclude that despite the seeming simplicity, a lot has yet to be understood on the systems-wide mechanisms connecting the driver mutation and the major cellular phenotypes of this pediatric cancer. Overall, this chapter can serve as a guide for a systems biology researcher to start working on EwS.

MicroRNA expression profiling and the role of ALCAM modulating tumor growth and metastasis in benzo[a]pyrene-transformed 16HBE cells

Benzo[a]pyrene (BaP) is a potent carcinogen and microRNAs (miRNAs) may play an important role in carcinogenesis. Activated leukocyte cell adhesion molecule (ALCAM) was up-regulated in BaP-transformed 16HBE cell line (THBEc1), and may be a key molecule for THBEc1 cells to gain and maintain the malignant phenotype. Here we screened the differentially expressed miRNAs which resulted in up-regulation of ALCAM in THBEc1 cells by comparing miRNA expression profiles between THBEc1 and 16HBE (HBE) cells. Results showed that a total of 555 miRNAs differentially expressed between THBEc1 and HBE cells, of which 351 miRNAs were down-regulated and 204 miRNAs were up-regulated in THBEc1 cells. MiR-152-3p, miR-142-5p and miR-211-5p down-regulated in THBEc1 cells were demonstrated to participate in the regulation of ALCAM. With THBEc1 as a tumor cell model, we determined the role of ALCAM in tumor growth and metastasis employing two ALCAM knockout THBEc1 cell lines via CRISPR/Cas9 technology. Results showed that ALCAM knockout inhibited colony formation and tumor growth, but enhanced cell migration and lung metastasis of THBEc1 cells. In conclusion, miR-152-3p/ALCAM, miR-142-5p/ALCAM and miR-211-5p/ALCAM axes may be involved in BaP-induced carcinogenesis. BaP might induce up-regulation of ALCAM via inhibiting miR-152-3p, miR-142-5p and miR-211-5p, which in turn allows ALCAM to exert its role promoting cell proliferation and tumor growth, and suppressing cell migration and metastasis.