肺腺癌预后相关miRNA生物信息学筛选及其临床意义

目的：探讨肺腺癌预后相关miRNA组学特征及其临床意义。方法：应用癌症基因组图谱（TCGA）数据库,检索人肺腺癌miRNA表达谱数据,进行差异分析,再利用Cox风险回归模型筛选预后相关miRNA;利用mirwalk分析平台,对筛选出的miRNA进行靶向调控基因预测、KEGG功能富集分析,最后,预测出预后相关miRNA的功能。结果：共筛选肺腺癌差异miRNA46个,其中,上调19个、下调27个;通过Cox生存分析筛选到预后相关miRNA有6个,即hsa-mir-21、hsa-mir-142、hsa-mir-200a高表达,hsa-mir-101、hsa-let-7c、hsa-mir-378e低表达,其中hsa-mir-21、hsa-mir-378e与肺腺癌患者不良预后有关,生存期显著缩短（P<0.05,AUC=0.618）。KEGG分析上述预后相关miRNA靶向调控基因与免疫反应通路、miRNA与癌症通路、代谢通路等有关。结论：hsa-mir-21、hsa-mir-378e与肺腺癌预后不良有关,未来经进一步临床验证有可能作为肺腺癌预后相关的分子标记物。     

Identification of a 4-miRNA signature as a potential prognostic biomarker for pancreatic adenocarcinoma

An microRNA (miRNA) signature to predict the clinical outcome of pancreatic adenocarcinoma (PAAD) is still lacking. In the current study, we aimed at identifying and evaluating a prognostic miRNA signature for patients with PAAD. The miRNA expression profile and the clinical information regarding patients with PAAD were recruited from The Cancer Genome Atlas database. Differentially expressed miRNAs were identified between normal and tumor samples. By means of survival analysis, a 4-miRNA signature for predicting patients' with PAAD overall survival (OS) was constructed. Receiver operating characteristic (ROC) analysis was applied to determine the efficiency of survival prediction. Furthermore, the biological function of the predicted miRNAs was evaluated using a bioinformatics approach. Four (hsa-mir-126, hsa-mir-3613, hsa-mir-424, and hsa-mir-4772) out of 17 differentially expressed miRNAs were associated to the OS of patients with PAAD. Moreover, the area under the curve (AUC) of the constructed 4-miRNA signature associated to patients' with PAAD 2-year survival was 0.789. The multivariate Cox's proportional hazards regression model suggested that this 4-miRNA signature was an independent prognostic factor of other clinical parameters in patients with PAAD. Further pathway enrichment analyses revealed that the miRNAs in the 4-miRNA signature might regulate genes that affect focal adhesion, Wnt signaling pathway, and PI3K-Akt signaling pathway. Thus, these findings indicated that the 4-miRNA signature might be an effective independent prognostic biomarker in the prediction of PAAD patients' survival.     

A Seven-microRNA Expression Signature Predicts Survival in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is the fifth common cancer. The differential expression of microRNAs (miRNAs) has been associated with the prognosis of various cancers. However, limited information is available regarding genome-wide miRNA expression profiles in HCC to generate a tumor-specific miRNA signature of prognostic values. In this study, the miRNA profiles in 327 HCC patients, including 327 tumor and 43 adjacent non-tumor tissues, from The Cancer Genome Atlas (TCGA) Liver hepatocellular carcinoma (LIHC) were analyzed. The associations of the differentially expressed miRNAs with patient survival and other clinical characteristics were examined with t-test and Cox proportional regression model. Finally, a tumor-specific miRNA signature was generated and examined with Kaplan–Meier survival, univariate\multivariate Cox regression analyses and KEGG pathway analysis. Results showed that a total of 207 miRNAs were found differentially expressed between tumor and adjacent non-tumor HCC tissues. 78 of them were also discriminatively expressed with gender, race, tumor grade and AJCC tumor stage. Seven miRNAs were significantly associated with survival (P value <0.001). Among the seven significant miRNAs, six (hsa-mir-326, hsa-mir-3677, hsa-mir-511-1, hsa-mir-511-2, hsa-mir-9-1, and hsa-mir-9-2) were negatively associated with overall survival (OS), while the remaining one (hsa-mir-30d) was positively correlated. A tumor-specific 7-miRNAs signature was generated and validated as an independent prognostic predictor. Collectively, we have identified and validated an independent prognostic model based on the expression of seven miRNAs, which can be used to assess patients’ survival. Additional work is needed to translate our model into clinical practice.     

Identifying a Polymorphic ‘Switch’ That Influences miRNAs' Regulation of a Myasthenia Gravis Risk Pathway

The significant roles of genetic variants in myasthenia gravis (MG) pathogenesis have been demonstrated in many studies, and recently it has been revealed that aberrant level/regulation of microRNAs (miRNAs) might contribute to the initiation and progression of MG. However, the dysfunction of miRNA associated with single nucleotide polymorphisms (miRSNPs) has not been well investigated in MG. In this study, we created a contemporary catalog of 89 MG risk genes via manual literature-mining. Based on this risk gene catalog, we obtained 18 MG risk pathways. Furthermore, we identified 93 miRNAs that target MG risk pathways and revealed the miRSNPs ‘switches’ in miRNA regulation in the MG risk pathways by integrating the database information of miRSNPs. We also constructed a miRNA-mediated SNP switching pathway network (MSSPN) to intuitively analyze miRNA regulation of MG risk pathways and the relationship of the polymorphism ‘switch’ with these changes in regulation. Moreover, we carried out in-depth dissection on the correlation between hsa05200 (pathway in cancer) and MG development, and elaborated the significance of 4 high-risk genes. By network analysis and literature mining, we proposed a potential mechanism of miRSNPs→gene→pathway effects on MG pathogenesis, especially for rs28457673 (miR-15/16/195/424/497 family)→IGF1R→hsa05200 (pathway in cancer). Therefore, our studies have revealed a functional role for genetic modulators in MG pathogenesis at a systemic level, which could be informative for further miRNA and miRSNPs studies in MG.     

Embryonic miRNA Profiles of Normal and Ectopic Pregnancies

Our objective was to investigate the miRNA profile of embryonic tissues in ectopic pregnancies (EPs) and controlled abortions (voluntary termination of pregnancy; VTOP). Twenty-three patients suffering from tubal EP and twenty-nine patients with a normal ongoing pregnancy scheduled for a VTOP were recruited. Embryonic tissue samples were analyzed by miRNA microarray and further validated by real time PCR. Microarray studies showed that four miRNAs were differentially downregulated (hsa-mir-196b, hsa-mir-30a, hsa-mir-873, and hsa-mir-337-3p) and three upregulated (hsa-mir-1288, hsa-mir-451, and hsa-mir-223) in EP compared to control tissue samples. Hsa-miR-196, hsa-miR-223, and hsa-miR-451 were further validated by real time PCR in a wider population of EP and control samples. We also performed a computational analysis to identify the gene targets and pathways which might be modulated by these three differentially expressed miRNAs. The most significant pathways found were the mucin type O-glycan biosynthesis and the ECM-receptor-interaction pathways. We also checked that the dysregulation of these three miRNAs was able to alter the expression of the gene targets in the embryonic tissues included in these pathways such as GALNT13 and ITGA2 genes. In conclusion, analysis of miRNAs in ectopic and eutopic embryonic tissues shows different expression patterns that could modify pathways which are critical for correct implantation, providing new insights into the understanding of ectopic implantation in humans.     

Downregulation of miR-31, miR-155, and miR-564 in chronic myeloid leukemia cells

Background/Aims

MicroRNAs (miRNAs) are short non-coding regulatory RNAs that control gene expression and play an important role in cancer development and progression. However, little is known about the role of miRNAs in chronic myeloid leukemia (CML). Our objective is to decipher a miRNA expression signature associated with CML and to determine potential target genes and signaling pathways affected by these signature miRNAs.

Results

Using miRNA microarrays and miRNA real-time PCR we characterized the miRNAs expression profile of CML cell lines and patients in reference to non-CML cell lines and healthy blood. Of all miRNAs tested, miR-31, miR-155, and miR-564 were down-regulated in CML cells. Down-regulation of these miRNAs was dependent on BCR-ABL activity. We next analyzed predicted targets and affected pathways of the deregulated miRNAs. As expected, in K562 cells, the expression of several of these targets was inverted to that of the miRNA putatively regulating them. Reassuringly, the analysis identified CML as the main disease associated with these miRNAs. MAPK, ErbB, mammalian target of rapamycin (mTOR) and vascular endothelial growth factor (VEGF) were the main molecular pathways related with these expression patterns. Utilizing Venn diagrams we found appreciable overlap between the CML-related miRNAs and the signaling pathways-related miRNAs.

Conclusions

The miRNAs identified in this study might offer a pivotal role in CML. Nevertheless, while these data point to a central disease, the precise molecular pathway/s targeted by these miRNAs is variable implying a high level of complexity of miRNA target selection and regulation. These deregulated miRNAs highlight new candidate gene targets allowing for a better understanding of the molecular mechanism underlying the development of CML, and propose possible new avenues for therapeutic treatment.     

Common miRNAs,candidate genes and their interaction network across four subtypes of epithelial ovarian cancer

Epithelial ovarian cancer (EOC) is categorized into four major histological subtypes such as clear cell carcinoma (CCC), endometrioid carcinoma (EC), mucinous carcinoma (MC), and serous carcinoma (SC). Heterogeneity of the EOC leads to different clinical outcomes of the disease, although all the subtypes are originated from the same layer of tissue. Therefore, it is of interest to identify the common candidate genes, miRNA and their interaction network in four the subtypes of EOC. A comparative gene expression analysis identified 248 common differentially expressed genes (DEGs) in the four subtypes of EOC. Identified common DEGs were found to be enriched in cancer specific pathways. A protein-protein interaction (PPI) network of the common DEGs were constructed, and subsequent module and survival analyses identified seven key candidate genes (CCNB1, CENPM, CEP55, RACGAP1, TPX2, UBE2C, and ZWINT). We also documented 10 key candidate miRNAs (hsa-mir-16-5p, hsa-mir-23b-3p, hsa-mir-34a-5p, hsa-mir-103a-3p, hsa-mir-107, hsa-mir-124-3p, hsa-mir-129-2-3p, hsa-mir-147a, hsa-mir-205-5p, and hsa-mir-195-5p) linked to the candidate genes. These derived data find application in the understanding of EOC.     

MicroRNA regulation of molecular networks mapped by global microRNA, mRNA, and protein expression in activated T lymphocytes

MicroRNAs (miRNAs) regulate specific immune mechanisms, but their genome-wide regulation of T lymphocyte activation is largely unknown. We performed a multidimensional functional genomics analysis to integrate genome-wide differential mRNA, miRNA, and protein expression as a function of human T lymphocyte activation and time. We surveyed expression of 420 human miRNAs in parallel with genome-wide mRNA expression. We identified a unique signature of 71 differentially expressed miRNAs, 57 of which were previously not known as regulators of immune activation. The majority of miRNAs are upregulated, mRNA expression of these target genes is downregulated, and this is a function of binding multiple miRNAs (combinatorial targeting). Our data reveal that consideration of this complex signature, rather than single miRNAs, is necessary to construct a full picture of miRNA-mediated regulation. Molecular network mapping of miRNA targets revealed the regulation of activation-induced immune signaling. In contrast, pathways populated by genes that are not miRNA targets are enriched for metabolism and biosynthesis. Finally, we specifically validated miR-155 (known) and miR-221 (novel in T lymphocytes) using locked nucleic acid inhibitors. Inhibition of these two highly upregulated miRNAs in CD4(+) T cells was shown to increase proliferation by removing suppression of four target genes linked to proliferation and survival. Thus, multiple lines of evidence link top functional networks directly to T lymphocyte immunity, underlining the value of mapping global gene, protein, and miRNA expression.     

Folate deficiency disturbs hsa‐let‐7 g level through methylation regulation in neural tube defects

Folic acid deficiency during pregnancy is believed to be a high‐risk factor for neural tube defects (NTDs). Disturbed epigenetic modifications, including miRNA regulation, have been linked to the pathogenesis of NTDs in those with folate deficiency. However, the mechanism by which folic acid‐regulated miRNA influences this pathogenesis remains unclear. It is believed that DNA methylation is associated with dysregulated miRNA expression. To clarify this issue, here we measured the methylation changes of 22 miRNAs in 57 human NTD cases to explore whether such changes are involved in miRNA regulation in NTD cases through folate metabolism. In total, eight of the 22 miRNAs tested reduced their methylation modifications in NTD cases, which provide direct evidence of the roles of interactions between DNA methylation and miRNA level in these defects. Among the findings, there was a significant association between folic acid concentration and hsa‐let‐7 g methylation level in NTD cases. Hypomethylation of hsa‐let‐7 g increased its own expression level in both NTD cases and cell models, which indicated that hsa‐let‐7 g methylation directly regulates its own expression. Overexpression of hsa‐let‐7 g, along with its target genes, disturbed the migration and proliferation of SK‐N‐SH cells, implying that hsa‐let‐7 g plays important roles in the prevention of NTDs by folic acid. In summary, our data suggest a relationship between aberrant methylation of hsa‐let‐7 g and disturbed folate metabolism in NTDs, implying that improvements in nutrition during early pregnancy may prevent such defects, possibly via the donation of methyl groups for miRNAs.     

Differential expression of basal microRNAs’ patterns in human dental pulp stem cells

MicroRNAs (miRNAs) are small non‐coding RNAs that regulate translation of mRNA into protein and play a crucial role for almost all biological activities. However, the identification of miRNAs from mesenchymal stem cells (MSCs), especially from dental pulp, is poorly understood. In this study, dental pulp stem cells (DPSCs) were characterized in terms of their proliferation and differentiation capacity. Furthermore, 104 known mature miRNAs were profiled by using real‐time PCR. Notably, we observed 19 up‐regulated miRNAs and 29 significantly down‐regulated miRNAs in DPSCs in comparison with bone marrow MSCs (BM‐MSCs). The 19 up‐regulated miRNAs were subjected to ingenuity analysis, which were composed into 25 functional networks. We have chosen top 2 functional networks, which comprised 10 miRNA (hsa‐miR‐516a‐3p, hsa‐miR‐125b‐1‐3p, hsa‐miR‐221‐5p, hsa‐miR‐7, hsa‐miR‐584‐5p, hsa‐miR‐190a, hsa‐miR‐106a‐5p, hsa‐mir‐376a‐5p, hsa‐mir‐377‐5p and hsa‐let‐7f‐2‐3p). Prediction of target mRNAs and associated biological pathways regulated by each of this miRNA was carried out. We paid special attention to hsa‐miR‐516a‐3p and hsa‐miR‐7‐5p as these miRNAs were highly expressed upon validation with qRT‐PCR analysis. We further proceeded with loss‐of‐function analysis with these miRNAs and we observed that hsa‐miR‐516a‐3p knockdown induced a significant increase in the expression of WNT5A. Likewise, the knockdown of hsa‐miR‐7‐5p increased the expression of EGFR. Nevertheless, further validation revealed the role of WNT5A as an indirect target of hsa‐miR‐516a‐3p. These results provide new insights into the dynamic role of miRNA expression in DPSCs. In conclusion, using miRNA signatures in human as a prediction tool will enable us to elucidate the biological processes occurring in DPSCs.     

A family of human microRNA genes from miniature inverted-repeat transposable elements

While hundreds of novel microRNA (miRNA) genes have been discovered in the last few years alone, the origin and evolution of these non-coding regulatory sequences remain largely obscure. In this report, we demonstrate that members of a recently discovered family of human miRNA genes, hsa-mir-548, are derived from Made1 transposable elements. Made1 elements are short miniature inverted-repeat transposable elements (MITEs), which consist of two 37 base pair (bp) terminal inverted repeats that flank 6 bp of internal sequence. Thus, Made1 elements are nearly perfect palindromes, and when expressed as RNA they form highly stable hairpin loops. Apparently, these Made1-related structures are recognized by the RNA interference enzymatic machinery and processed to form 22 bp mature miRNA sequences. Consistent with their origin from MITEs, hsa-mir-548 genes are primate-specific and have many potential paralogs in the human genome. There are more than 3,500 putative hsa-mir-548 target genes; analysis of their expression profiles and functional affinities suggests cancer-related regulatory roles for hsa-mir-548. Taken together, the characteristics of Made1 elements, and MITEs in general, point to a specific mechanism for the generation of numerous small regulatory RNAs and target sites throughout the genome. The evolutionary lineage-specific nature of MITEs could also provide for the generation of novel regulatory phenotypes related to species diversification. Finally, we propose that MITEs may represent an evolutionary link between siRNAs and miRNAs.     

Interplay between SARS-CoV-2-derived miRNAs,immune system,vitamin D pathway and respiratory system

The new coronavirus pandemic started in China in 2019. The intensity of the disease can range from mild to severe, leading to death in many cases. Despite extensive research in this area, the exact molecular nature of virus is not fully recognized; however, according to pieces of evidence, one of the mechanisms of virus pathogenesis is through the function of viral miRNAs. So, we hypothesized that SARS-CoV-2 pathogenesis may be due to targeting important genes in the host with its miRNAs, which involved in the respiratory system, immune pathways and vitamin D pathways, thus possibly contributing to disease progression and virus survival. Potential miRNA precursors and mature miRNA were predicted and confirmed based on the virus genome. The next step was to predict and identify their target genes and perform functional enrichment analysis to recognize the biological processes connected with these genes in the three pathways mentioned above through several comprehensive databases. Finally, cis-acting regulatory elements in 5′ regulatory regions were analysed, and the analysis of available RNAseq data determined the expression level of genes. We revealed that thirty-nine mature miRNAs could theoretically derive from the SARS-CoV-2 genome. Functional enrichment analysis elucidated three highlighted pathways involved in SARS-CoV-2 pathogenesis: vitamin D, immune system and respiratory system. Our finding highlighted genes' involvement in three crucial molecular pathways and may help develop new therapeutic targets related to SARS-CoV-2.