Identification of potential miRNA biomarkers for traumatic osteonecrosis of femoral head

Traumatic osteonecrosis of femoral head (TONFH) is a common orthopedic disease caused by physical injury in hip. However, the unclear pathogenesis mechanism of TONFH and lacking of simple noninvasive early diagnosis method cause the necessity of hip replacement for most patients with TONFH. In this study, we aimed to identify circulating microRNAs (miRNAs) by integrated bioinformatics analyses as potential biomarker of TONFH. mRNA expression profiles were downloaded from the Gene Expression Omnibus database. Then we combined two miRNA screen methods: Weighted gene co-expression network analysis and fold change based differentially expressed miRNAs analysis. As a result, we identified 14 key miRNAs as potential biomarkers for TONFH. Besides, 302 target genes of these miRNAs were obtained and the miRNA–mRNA interaction network was constructed. Furthermore, the results of Kyoto Encyclopedia of Gene and Genome pathway analysis, Gene Ontology function analysis, protein–protein interaction (PPI) network analysis and PPI network module analysis showed close correlation between these 14 key miRNAs and TONFH. Then we established receiver operating characteristic curves and identified 6-miRNA signature with highly diagnosis value including miR-93-5p (area under the curve [AUC] = 0.93), miR-1324 (AUC = 0.92), miR-4666a-3p (AUC = 0.92), miR-5011-3p (AUC = 0.92), and miR-320a (AUC = 0.89), miR-185-5p (AUC = 0.89). Finally, the results of quantitative real-time polymerase chain reaction confirmed the significantly higher expression of miR-93-5p and miR-320a in the serum of patients with ONFH. These circulating miRNAs could serve as candidate early diagnosis markers and potential treatment targets of TONFH.     

Identification of key microRNAs and hub genes in non-small-cell lung cancer using integrative bioinformatics and functional analyses

Non-small-cell lung cancer (NSCLC) is an extremely debilitating respiratory malignancy. However, the pathogenesis of NSCLC has not been fully clarified. The main objective of our study was to identify potential microRNAs (miRNAs) and their regulatory mechanism in NSCLC. Using a systematic review, two NSCLC-associated miRNA data sets (GSE102286 and GSE56036) were obtained from Gene Expression Omnibus, and the differentially expressed miRNAs (DE-miRNAs) were accessed by GEO2R. Survival analysis of candidate DE-miRNAs was conducted using the Kaplan-Meier plotter database. To further illustrate the roles of DE-miRNAs in NSCLC, their potential target genes were predicted by miRNet and were annotated by the Database for Annotation, Visualization and Integrated Discovery (DAVID) program. Moreover, the protein-protein interaction (PPI) and miRNA-hub gene regulatory network were established using the STRING database and Cytoscape software. The function of DE-miRNAs in NSCLC cells was evaluated by transwell assay. Compared with normal tissues, a total of eight DE-miRNAs was commonly changed in two data sets. The survival analysis showed that six miRNAs (miR-21-5p, miR-31-5p, miR-708-5p, miR-30a-5p, miR-451a, and miR-126-3p) were significantly correlated with overall survival. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that target genes of upregulated miRNAs were enriched in pathways in cancer, microRNAs in cancer and proteoglycans in cancer, while the target genes of downregulated miRNAs were mainly associated with pathways in cancer, the PI3K-Akt signaling pathway and HTLV-I infection. Based on the miRNA-hub gene network and expression analysis, PTEN, EGFR, STAT3, RHOA, VEGFA, TP53, CTNNB1, and KRAS were identified as potential target genes. Furthermore, all six miRNAs exhibited significant effects on NSCLC cell invasion. These findings indicate that six DE-miRNAs and their target genes may play important roles in the pathogenesis of NSCLC, which will provide novel information for NSCLC treatments.     

A panel of two miRNAs correlated to systolic blood pressure is a good diagnostic indicator for stroke

Background: We aimed to develop a diagnostic indicator of stroke based on serum miRNAs correlated to systolic blood pressure.Methods: Using miRNA expression profiles in {"type":"entrez-geo","attrs":{"text":"GSE117604","term_id":"117604"}}GSE117604 from the Gene Expression Omnibus (GEO), we utilized the WGCNA to identify hub miRNAs correlated to systolic blood pressure (SBP). Differential analysis was applied to highlight hub differentially expressed miRNAs (DE-miRNAs), whereby we built a miRNA-based diagnostic indicator for stroke using bootstrap ranking Least Absolute Shrinkage and Selection Operator (LASSO) regression with 10-fold cross-validation. The classification value of the indicator was validated with receiver operating characteristic (ROC) analysis in both the training set and test set, as well as quantitative real-time PCR (qRT-PCR) for the feature miRNAs. Further, target genes of hub miRNAs and hub DE-miRNAs were retrieved for functional enrichment.Results: A total of 447 hub miRNAs in the blue modules were significantly correlated with systolic blood pressure (r = 0.32, false discovery rate = 10⁻⁶). Target genes predicted with the hub miRNAs were mostly implicated in the Kyoto Encyclopedia of Genes and Genomes (KEGG) terms including mitogen-activated protein kinase (MAPK) pathway, senescence, and TGF-β signaling pathway. The diagnostic indicator with miR-4420 and miR-6793-5p showed remarkable performance in the training set (area under curve [AUC]= 0.953), as well as in the test set (AUC = 0.894). Results of qRT-PCR validated the diagnostic value of the two miRNAs embedded in the proposed indicator.Conclusions: We developed a panel of two miRNAs, which is a good diagnostic indicator for stroke. These results require further investigation.     

Construction of miRNA-mRNA network for the identification of key biological markers and their associated pathways in IgA nephropathy by employing the integrated bioinformatics analysis

BackgroundAbout half-century ago, Immunoglobulin A nephropathy (IgAN) was discovered as a complicated disease with frequent clinical symptoms. Until now, exact mechanism underlying the pathogenesis of IgAN is poorly known. Therefore, current study was aimed to understand the molecular mechanism of IgAN by identifying the key miRNAs and their targeted hub genes. The key miRNAs might contribute to the diagnosis and therapy of IgAN, and could turn out to be a new star in the field of IgAN.MethodsThe microarray datasets were downloaded from Gene Expresssion Omnibus (GEO) database and analyzed using R package (LIMMA) in order to obtain differential expressed genes (DEGs). Then, the hub genes were identified using cytoHubba plugin of cytoscpae tool and other bioinformatics approaches including protein-protein interaction (PPI) network analysis, module analysis, and miRNA-hub gene network construction was also performed.ResultsA total of 348 DEGs were identified, of which 107 were upregulated genes and 241 were downregulated genes. Subsequently, the 12 overlapped genes were predicted from cytoHubba, and considered as hub genes. Moreover, a network among miRNA-hub genes was created to explore the correlation between the hub genes and their targeted miRNAs. Network construction ultimately lead to the identification of nine gene named FN1, EGR1, FOS, JUN, SERPINE1, MMP2, ATF3, MYC, and IL1B and one novel key miRNA namely, has-miR-144-3p as biomarker for diagnosis and therapy of IgAN.ConclusionThis study updates the information and yield a new perspective in context of understanding the pathogenesis and development of IgAN. In future, key miRNAs might be capable of improving the personalized detection and therapies for IgAN. In vivo and in vitro investigation of miRNAs and pathway interaction is essential to delineate the specific roles of the novel miRNAs, which may help to further reveal the mechanisms underlying IgAN.     

The mRNA–miRNA–lncRNA Regulatory Network and Factors Associated with Prognosis Prediction of Hepatocellular Carcinoma

The aim of this study was to identify novel prognostic mRNA and microRNA (miRNA) biomarkers for hepatocellular carcinoma (HCC) using methods in systems biology. Differentially expressed mRNAs, miRNAs, and long non-coding RNAs (lncRNAs) were compared between HCC tumor tissues and normal liver tissues in The Cancer Genome Atlas (TCGA) database. Subsequently, a prognosis-associated mRNA co-expression network, an mRNA–miRNA regulatory network, and an mRNA–miRNA–lncRNA regulatory network were constructed to identify prognostic biomarkers for HCC through Cox survival analysis. Seven prognosis-associated mRNA co-expression modules were obtained by analyzing these differentially expressed mRNAs. An expression module including 120 mRNAs was significantly correlated with HCC patient survival. Combined with patient survival data, several mRNAs and miRNAs, including CHST4, SLC22A8, STC2, hsa-miR-326, and hsa-miR-21 were identified from the network to predict HCC patient prognosis. Clinical significance was investigated using tissue microarray analysis of samples from 258 patients with HCC. Functional annotation of hsa-miR-326 and hsa-miR-21-5p indicated specific associations with several cancer-related pathways. The present study provides a bioinformatics method for biomarker screening, leading to the identification of an integrated mRNA–miRNA–lncRNA regulatory network and their co-expression patterns in relation to predicting HCC patient survival.     

Exosomal small RNA sequencing uncovers the microRNA dose markers for power frequency electromagnetic field exposure

Purpose: The potential health risks caused by power frequency electromagnetic field (PFEMF) have led to increase public health concerns. However, the diagnosis and prognosis remain challenging in determination of exact dose of PFEMF exposure.

Materials and methods: Mice were exposed to different magnetic doses of PFEMF for the following isolation of serum exosomes, microRNAs (miRNAs) extraction and small RNA sequencing. After small RNA sequencing, bioinformatic analysis, quantitative real-time PCR (qRT-PCR) validation and serum exosomal miRNA biomarkers were determined.

Results: Significantly changed serum exosomal miRNA as biomarkers of 0.1, 0.5, 2.5?mT and common PFEMF exposure were confirmed. Gene ontology (GO) and Kyoto encyclopaedia of genes and genomes (KEGG) pathway analysis of the downstream target genes of the above-identified exosomal miRNA markers indicated that, exosomal miRNA markers were predicted to be involved in critical pathophysiological processes of neural system and cancer- or other disease-related signalling pathways.

Conclusions: Aberrantly-expressed serum exosomal miRNAs, including miR-128-3p for 0.1?mT, miR-133a-3p for 0.5?mT, miR-142a-5p for 2.5?mT, miR-218-5p and miR-199a-3p for common PFEMF exposure, suggested a series of informative markers for not only identifying the exact dose of PFEMF exposure, also consolidating the base for future clinical intervention.     

Predicting functional circular RNA-based competitive endogenous RNA network in gastric carcinoma using novel bioinformatics analysis

Gastric cancer (GC) remains one of the most prevalent types of malignancies worldwide, and also one of the most reported lethal tumor-related diseases. Circular RNAs (circRNAs) have been certified to be trapped in multiple aspects of GC pathogenesis. Yet, the mechanism of this regulation is mostly undefined. This research is designed to discover the vital circRNA-microRNA (miRNA)-messenger RNA (mRNA) regulatory network in GC. Expression profiles with diverse levels including circRNAs, miRNAs, and mRNAs were all determined using microarray public datasets from Gene Expression Ominous (GEO). The differential circRNAs expressions were recognized against the published robust rank aggregation algorithm. Besides, a circRNA-based competitive endogenous RNA (ceRNA) interaction network was visualized via Cytoscape software (version 3.8.0). Functional and pathway enrichment analysis associated with differentially expressed targeted mRNAs were conducted using Cytoscape and an online bioinformatics database. Furthermore, an interconnected protein–protein interaction association network which consisted of 51 mRNAs was predicted, and hub genes were screened using STRING and CytoHubba. Then, several hub genes were chosen to explore their expression associated with survival rate and clinical stage in GEPIA and Kaplan-Meier Plotter databases. Finally, a carefully designed circRNA-related ceRNA regulatory subnetwork including four circRNAs, six miRNAs, and eight key hub genes was structured using the online bioinformatics tool.     

MicroRNA-384-5p regulates ischemia-induced cardioprotection by targeting phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit delta (PI3K p110δ)

MicroRNAs (miRNAs) are a novel class of powerful, endogenous regulators of gene expression. This study identified 16 differentially expressed miRNAs in ischemic myocardium of rats using TaqMan Low Density Array. In addition, bioinformatics analyses, such as Gene ontology and Pathway assays, were applied to determine the apoptosis pathway, only regulated by miR-384-5p, and all the associated target genes (PIK3CD, PPP3CA, PPP3CB, PPP3R1, CASP3 and IL1A). These target genes, besides PIK3CB, were shown to be significantly up-regulated by qRT-PCR assay, which further suggested that PIK3CD, PPP3CA, PPP3R1, CASP3, IL1A could be regulated by miR-384-5p. MTT, Western blot, qRT-PCR and luciferase assays were used to investigate the role of miR-384-5p in myocardial ischemia. We found that cleaved caspase3 expression was up-regulated by miR-384-5p and down-regulated by miR-384-5p inhibitor suggesting that apoptosis pathway was regulated by miR-384-5p. We also found that miR-384-5p suppressed cell viability while miR-384-5p inhibitor improved it, confirming H9c2 cell survival was affected by miR-384-5p. In addition, the PIK3CD protein level in H9c2 cells was up-regulated by miR-384-5p inhibitor. We found that miR-384-5p expression level decreased and PIK3CD protein level increased in both ischemic myocardium of rats and hypoxic H9c2 cells, and that miR-384-5p suppress PIK3CD expression through a miR-384-5p binding site within the 3′ untranslational region of PIK3CD. These results show that miR-384-5p, an important protecting factor, plays a significant role in cardioprotection by regulating PIK3CD in myocardial ischemia.     

Upregulated microRNAs in membranous glomerulonephropathy are associated with significant downregulation of IL6 and MYC mRNAs

Membranous glomerulonephropathy (MGN) is a glomerulopathy characterized by subepithelial deposits of immune complexes on the extracapillary side of the glomerular basement membrane. Insertion of C5b-9 (complement membrane-attack complex) into the membrane leads to functional impairment of the glomerular capillary wall. Knowledge of the molecular pathogenesis of MGN is actually scanty. MicroRNA (miRNA) profiling in MGN and unaffected tissues was performed by TaqMan Low-Density Arrays. Expression of miRNAs and miRNA targets was evaluated in Real-Time polymerase chain reaction (PCR). In vitro transient silencing of miRNAs was achieved through transfection with miRNA inhibitors. Ten miRNAs (let-7a-5p, let-7b-5p, let-7c-5p, let-7d-5p, miR-107, miR-129-3p, miR-423-5p, miR-516-3p, miR-532-3p, and miR-1275) were differentially expressed (DE) in MGN biopsies compared to unaffected controls. Interleukin 6 (IL6) and MYC messenger RNAs (mRNAs; targets of DE miRNAs) were significantly downregulated in biopsies from MGN patients, and upregulated in A498 cells following let-7a-5p or let-7c-5p transient silencing. Gene ontology analysis showed that DE miRNAs regulate pathways associated with MGN pathogenesis, including cell cycle, proliferation, and apoptosis. A significant correlation between DE miRNAs and mRNAs and clinical parameters (i.e., antiphospholipid antibodies, serum creatinine, estimated glomerular filtration, proteinuria, and serum cholesterol) has been detected. Based on our data, we propose that DE miRNAs and their downstream network may be involved in MGN pathogenesis and could be considered as potential diagnostic biomarkers of MGN.     

miR-425-5p as an exosomal biomarker for metastatic prostate cancer

Prostate cancer-related deaths are mostly caused by metastasis, which indicates the importance of identifying clinical prognostic biomarkers. In this study, we evaluated the expression profile of exosomal microRNAs (miRNAs) derived from metastatic prostate cancer (mPCa) cell lines (LNCaP and PC-3). miRNA signatures in exosomes and cells were evaluated by miRNA microarray analysis. Fourteen miRNAs were identified as candidates for specific noninvasive biomarkers. The expression of five miRNAs was validated using RT-qPCR, which confirmed that miR-205-5p, miR-148a-3p, miR-125b-5p, miR-183-5p, and miR-425-5p were differentially expressed in mPCa exosomes. Bioinformatic analyses showed that miR-425-5p was associated with residual tumor, pathologic T and N stages, and TP53 status in PCa samples. Gene ontology analysis of negatively correlated and predicted targeted genes showed enrichment of genes related to bone development pathways. The LinkedOmics database indicated that the potential target HSPB8 has a significant negative correlation with miR-425-5p. In conclusion, this study identified a panel of exosomal miRNAs with potential value as prognostic biomarkers for prostate cancer.     

Screening candidate microRNA-mRNA network for predicting the response to chemoresistance in osteosarcoma by bioinformatics analysis

The search for biomarkers is important for providing more targeted treatments for osteosarcoma patients with chemoresistance. In this study, differentially expressed microRNAs (miRNAs) were identified from miRNA expression profiles. And the target messenger RNAs (mRNAs) of miRNA were obtained from two websites in public domains. Analysis of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway by these miRNA targets suggests that they may have potential links to osteosarcoma chemoresistance. In the protein-protein interaction (PPI) network, we screened three subnetworks and 10 hub RNAs, and analyzed through KEGG pathway and searched the PubMed database, indicating that they were significantly associated with drug resistance. Then we found 12 key mRNAs by analyzing the mRNA expression profile. Survival analyses showed that most of the 10 hub mRNAs and 12 key mRNAs had a significant influence on the prognosis of patients with chemoresistance osteosarcoma. A miRNA-mRNA network is constructed by integrating mRNAs and miRNAs information. The network biomarkers in this study have an advantage over traditional single-molecule biomarkers in terms of predictive power. And the mRNAs in this network biomarkers are supported by survival analysis or by existing theories. These results will contribute to the choice of chemotherapy before treatment and the prediction of patient prognosis.     

MicroRNA Profiles in Spontaneous Decidualized Menstrual Endometrium and Early Pregnancy Decidua with Successfully Implanted Embryos

To comparatively analyze the human microRNA (miRNA) profiles between spontaneous decidualized menstrual endometrium and early pregnancy decidua by an in-depth sequencing of miRNAs. The specific miRNAs expressed at conception might be involved in pregnancy establishment and expression of let-7f-5p and let-7g-5p was experimentally up-regulated or inhibited to assess the effect on the expression of IGF2BP-1 and IGF2R in vitro, respectively. Samples of endometria and deciduas were obtained from 25 women who suffered from tubal or male factor subfertility and from 35 early pregnant women who underwent pregnancy termination at 6–8 weeks gestation were irrespectively collected and comparatively analyzed by miRNA sequencing and differential expression of known and novel miRNAs was analyzed using bioinformatics. The 2042 miRNA expression was analyzed in the study and the differential expression of six miRNAs was validated by qRT-PCR. The expression of four miRNAs in decidua samples was down-regulated (miR-34c, miR-92a, miR-181a-5p, and miR-191), whereas the expression of miR-10a-5p and let-7f-5p was significantly up-regulated. The expression of IGF2BP-1 and IGF2R declined and increased with overexpression and inhibition of let-7f-5p and let-7g-5p, respectively. Changes in the expression of particular miRNAs might play a role in the physiology of decidualization following successful embryo implantation, ultimately resulting in continuous decidualization.     

Comprehensive gene and microRNA expression profiling reveals the crucial role of hsa-let-7i and its target genes in colorectal cancer metastasis

Accumulating evidence has demonstrated that miRNAs play important roles in the occurrence and development of colorectal cancer (CRC). However, whether miRNAs are associated with the metastasis of CRC remains largely unexplored. The aim of the current study is to profile miRNAs in different CRC metastatic cell lines to identify the biomarkers in CRC metastasis. Gene and miRNA expression profiling was performed to analyze the global expression of mRNAs and miRNAs in the four human CRC cell lines (LoVo, SW480, HT29 and Caco-2) with different potential of metastasis. Expression patterns of mRNAs and miRNAs were altered in different CRC cell lines. By developing an integrated bioinformatics analysis of gene and miRNA expression patterns, hsa-let-7i was identified to show the highest degree in the microRNA-GO-network and microRNA-Gene-network. The expression level of hsa-let-7i was further validated by qRT-PCR in CRC cells. In addition, the targets of hsa-let-7i were predicted by two programs TargetScan and PicTar, and target genes were validated by expression profiling in the most epresentative LoVo and Caco-2 cell lines. Eight genes including TRIM41, SOX13, SLC25A4, SEMA4F, RPUSD2, PLEKHG6, CCND2, and BTBD3 were identified as hsa-let-7i targets. Our data showed the power of comprehensive gene and miRNA expression profiling and the application of bioinformatics tools in the identification of novel biomarkers in CRC metastasis.     

Identification of Differential MicroRNAs in Cerebrospinal Fluid and Serum of Patients with Major Depressive Disorder

Major depression is a debilitating disease. To date, the development of biomarkers of major depressive disorder (MDD) remains a challenge. Recently, alterations in the expression of microRNAs (miRNAs) from post-mortem brain tissue and peripheral blood have been linked to MDD. The goals of this study were to detect the differential miRNAs in cerebrospinal fluid (CSF) and serum of MDD patients. First, the relative expression levels of 179 miRNAs (relative high levels in serum) were analyzed by miRNA PCR Panel in the CSF of MDD patients. Then, the differentially altered miRNAs from CSF were further assessed by qRT-PCR in the serum of the same patients. Finally, the serum differentially altered miRNAs were further validated by qRT-PCR in the serum of another MDD patients. The CSF-results indicated that 11 miRNAs in MDD patients were significantly higher than these in control subjects, and 5 miRNAs were significantly lower than these in control subjects. The serum-results from the same patients showed that 3 miRNAs (miR-221-3p, miR-34a-5p, and let-7d-3p) of the 11 miRNAs were significantly higher than these in control subjects, and 1 miRNA (miR-451a) of 5 miRNAs was significantly lower than these in control subjects. The up-regulation of miR-221-3p, miR-34a-5p, let-7d-3p and down-regulation of miR-451a was further validated in another 32 MDD patients. ROC analysis showed that the area under curve of let-7d-3p, miR-34a-5p, miR-221-3p and miR-451a was 0.94, 0.98, 0.97 and 0.94, with specificity of 90.48%, 95.24%, 90.48% and 90.48%, and sensitivity of 93.75%, 96.88%, 90.63% and 84.85%, respectively. In addition, target gene prediction found that the altered miRNAs are involved in affecting some important genes and pathway related to MDD. Our results suggested that differentially altered miRNAs in CSF might be involved in MDD, and serum miR-221-3p, miR-34a-5p, let-7d-3p, and miR-451a might be able to serve as biomarkers for MDD.