Dependence of Intracellular and Exosomal microRNAs on Viral E6/E7 Oncogene Expression in HPV-positive Tumor Cells

Specific types of human papillomaviruses (HPVs) cause cervical cancer. Cervical cancers exhibit aberrant cellular microRNA (miRNA) expression patterns. By genome-wide analyses, we investigate whether the intracellular and exosomal miRNA compositions of HPV-positive cancer cells are dependent on endogenous E6/E7 oncogene expression. Deep sequencing studies combined with qRT-PCR analyses show that E6/E7 silencing significantly affects ten of the 52 most abundant intracellular miRNAs in HPV18-positive HeLa cells, downregulating miR-17-5p, miR-186-5p, miR-378a-3p, miR-378f, miR-629-5p and miR-7-5p, and upregulating miR-143-3p, miR-23a-3p, miR-23b-3p and miR-27b-3p. The effects of E6/E7 silencing on miRNA levels are mainly not dependent on p53 and similarly observed in HPV16-positive SiHa cells. The E6/E7-regulated miRNAs are enriched for species involved in the control of cell proliferation, senescence and apoptosis, suggesting that they contribute to the growth of HPV-positive cancer cells. Consistently, we show that sustained E6/E7 expression is required to maintain the intracellular levels of members of the miR-17~92 cluster, which reduce expression of the anti-proliferative p21 gene in HPV-positive cancer cells. In exosomes secreted by HeLa cells, a distinct seven-miRNA-signature was identified among the most abundant miRNAs, with significant downregulation of let-7d-5p, miR-20a-5p, miR-378a-3p, miR-423-3p, miR-7-5p, miR-92a-3p and upregulation of miR-21-5p, upon E6/E7 silencing. Several of the E6/E7-dependent exosomal miRNAs have also been linked to the control of cell proliferation and apoptosis. This study represents the first global analysis of intracellular and exosomal miRNAs and shows that viral oncogene expression affects the abundance of multiple miRNAs likely contributing to the E6/E7-dependent growth of HPV-positive cancer cells.     

Profiles of serum microRNAs; miR-125b-5p and miR223-3p serve as novel biomarkers for HBV-positive hepatocellular carcinoma

Recently, circulating miRNAs have been reported as promising biomarkers for various pathologic conditions including cancer. Certain microRNAs (miRNAs) have been shown early diagnostic potential for many types of cancer. The objective of this study was to investigate the potential of certain serum/plasma miRNAs as novel non-invasive biomarkers for early diagnosis of hepatitis B virus (HBV) related hepatocellular carcinoma (HCC). For this reason, the expression levels of 24 miRNA (let-7c, miR-92a-3p, 423-5p, 150-5p, 223-3p, 125b-5p, 342-3p, miR-206, 122-5p, 375, 223-5p, 10a-5p, 23b-5p, 99a-5p, 23a-5p, 10a-3p, 122-3p, 125b-1-3p, 23b-3p, 125b-2-3p, 23a-3p, 92a-1-5p, 92a-2-5p, 99a-3p) were analyzed in plasma of patients with chronic hepatitis B, HBV-positive cirrhosis and HBV-positive HCC and compared with control group samples. Totally 94 plasma samples; 28 control and 66 patient plasma (24 CHB, 22 HBV-positive cirrhosis, 20 HBV-positive HCC) and were included in this study. The expression levels of 24 miRNAs were detected for all control and patient group plasma samples by qRT-PCR using BioMark? 96.96 Dynamic Array (Fluidigm Corporation) system. The expression levels of miR-125b-5p were detected 2.85 fold, 2.46 fold and 1.89 fold (p = 0.01513, p = 0.0009440, p = 0.0001446) up regulated in CHB, HBV-positive cirrhosis and HBV-positive HCC, respectively when compared versus control group individually by Mann–Whitney U test. The expression levels of miR-223-3p were detected 5.55 fold, 13.88 fold and 12.65 fold (p = 0.01513, p = 0.0009440, p = 0.0001446) down regulated in same comparisons. When all groups were compared versus control group by one-way ANOVA test, the expression levels of miR-223-3p were also found statistically significant (p < 0.05). Although not statistically significant, miR-125b-5p tended to be upregulated. (p = 0.07192). These results significantly imply that miR-125b-5p and miR223-3p could be used as novel non-invasive biomarkers of HBV-positive HCC in very early, even at CHB stage of liver disease.     

Inheritable changes in miRNAs expression in HeLa cells after X-ray and mitomycin C treatment

We identified 40 miRNAs with inherited aberrant expression by multiple parallel sequencing of human HeLa cells irradiated with X rays and mitomycin C. Twenty-two miRNAs were repressed and 15 miRNAs were induced after radiation and mytomycin C treatment. The expression of three miRNAs (miR-10b-5p, miR-148a-3p, and miR-340-5p) decreased after X-ray exposure and increased after mitomycin C treatment. The spectrum of aberrantly expressed miRNAs after X-ray and mitomycin C treatment is different, except for three miRNAs (mir-100-5p, miR-99b-5p, miR-501-3p), which showed the inherited decreased expression after both mutagens. It has been ascertained that for five miRNAs (miR-21-3p, miR-182-5p, miR-19b-3p, miR-30a-3p, and miR-30e-3p) with increased inherited expression, the targets are well-described tumor suppressor genes. For 9 miRNAs (miR-99b-5p, miR-148a-3p, miR-365a-3p, miR-193a-3p, miR-100-5p, miR-99a-5p, miR-29b-3p, miR-340-5p, and miR-23b-3p) with reduced inherited expression, the targets are oncogenes. The obtained results provide further support of the idea that induced epigenetic changes in the genome should be considered when assessing the long-term genetic effects of ionizing radiation and chemical compounds.     

Circulating miR-23b-3p,miR-145-5p and miR-200b-3p are potential biomarkers to monitor acute pain associated with laminitis in horses

Circulating microRNAs (miRNAs) are emerging as promising biomarkers for several disorders and related pain. In equine practice, acute laminitis is a common disease characterised by intense pain that severely compromises horse welfare. Recently, the Horse Grimace Scale (HGS), a facial expression-based pain coding system, was shown to be a valid welfare indicator to identify pain linked to acute laminitis. The present study aimed to: determine whether miRNAs can be used as biomarkers for acute pain in horses (Equus caballus) affected by laminitis; integrate miRNAs to their target genes and to categorise target genes for biological processes; gather additional evidence on concurrent validity of HGS by investigating how it correlates to miRNAs. Nine horses presenting acute laminitis with no prior treatment were recruited. As control group, nine healthy horses were further included in the experimental design. Samples were collected from horses with laminitis at admission before any treatment (‘pre-treatment’) and 7 days after routine laminitis treatment (‘post-treatment’). The expression levels of nine circulating miRNAs, namely hsa-miR-532-3p, hsa-miR-219-5p, mmu-miR-134-5p, mmu-miR-124a-3p, hsa-miR-200b-3p, hsa-miR-146a-5p, hsa-miR-23b-3p, hsa-miR-145-5p and hsa-miR-181a-5p, were detected and assessed as potential biomarkers of pain by quantitative PCR using TaqMan® probes. The area under the receiver operating curve (AUC) was then used to evaluate the diagnostic performance of miRNAs. Molecular data were integrated with HGS scores assessed by one trained treatment and time point blind veterinarian. The comparative analysis demonstrated that the levels of miR-23b-3p (P=0.029), miR-145-5p (P=0.015) and miR-200b-3p (P=0.023) were significantly higher in pre-treatment and the AUCs were 0.854, 0.859 and 0.841, respectively. MiR-200b-3p decreased after routine laminitis treatment (P=0.043). Combining two miRNAs in a panel, namely miR-145-5p and miR-200b-3p, increased efficiency in distinguishing animals with acute pain from controls. In addition, deregulated miRNAs were positively correlated to HGS scores. Computational target prediction and functional enrichment identified common biological pathways between different miRNAs. In particular, the glutamatergic pathway was affected by all three miRNAs, suggesting a crucial role in the pathogenesis of pain. In conclusion, the dynamic expression of circulating miR-23b-3p, miR-145-5p and miR-200b-3p was detected in horses with acute laminitis and miRNAs can be considered potentially promising pain biomarkers. Further studies are needed in order to assess their relevancy in other painful conditions severely compromising horse welfare. An important implication would be the possibility to use them for the concurrent validation of non-invasive indicators of pain in horses.     

Aberrant expression of microRNA induced by high-fructose diet: implications in the pathogenesis of hyperlipidemia and hepatic insulin resistance

Fructose is a highly lipogenic sugar that can alter energy metabolism and trigger metabolic disorders. In the current study, microRNAs (miRNAs) altered by a high-fructose diet were comprehensively explored to elucidate their significance in the pathogenesis of chronic metabolic disorders. miRNA expression profiling using small noncoding RNA sequencing revealed that 19 miRNAs were significantly upregulated and 26 were downregulated in the livers of high-fructose-fed mice compared to chow-fed mice. Computational prediction and functional analysis identified 10 miRNAs, miR-19b-3p, miR-101a-3p, miR-30a-5p, miR-223-3p, miR-378a-3p, miR-33-5p, miR-145a-3p, miR-128-3p, miR-125b-5p and miR-582-3p, assembled as a regulatory network to potentially target key genes in lipid and lipoprotein metabolism and insulin signaling at multiple levels. qRT-PCR analysis of their potential target genes [IRS-1, FOXO1, SREBP-1c/2, ChREBP, insulin-induced gene-2 (Insig-2), microsomal triglyceride transfer protein (MTTP) and apolipoprotein B (apoB)] demonstrated that fructose-induced alterations of miRNAs were also reflected in mRNA expression profiles of their target genes. Moreover, the miRNA profile induced by high-fructose diet differed from that induced by high-fat diet, indicating that miRNAs mediate distinct pathogenic mechanisms in dietary-induced metabolic disorders. This study presents a comprehensive analysis of a new set of hepatic miRNAs, which were altered by high-fructose diet and provides novel insights into the interaction between miRNAs and their target genes in the development of metabolic syndrome.     

miRNA expression profiling uncovers a role of miR-302b-3p in regulating skin fibroblasts senescence

Numbers of emerging evidence suggest that variable microRNA (miRNA) expression facilitates the aging process. In this study, we distinguished aberrant miRNA expression in aged skin and explored the biological functions and potential mechanism of upregulated miR-302b-3p. At first, miRNA microarray analysis was examined to explore miRNA expression profiling in the skin of aging mice model by D -galactose (d -gal) injection. We identified 29 aberrant miRNAs in aged mice skin. Next, KEGG enrichment analysis was conducted with DIANA-miPath v3.0, which was revealed that enrichment pathways involved in such processes as extracellular matrix-receptor interaction, MAPK signaling pathway, and mammalian target of rapamycin (mTOR) signaling pathway. The target genes of deregulated miRNAs were predicted from four bioinformatic algorithms (miRDB, Targetscan, miRwalk, and Tarbase). The interaction network of miRNAs and their targets were visualized using Cytoscape software. As a result, we found that some hub genes (including JNK2, AKT1/2/3, PAK7, TRPS1, BCL2L11, and IKZF2) were targeted by 12 potential miRNAs (including miR-302b-3p, miR-291a-5p, miR-139-3p, miR-467c-3p, miR-186-3p, etc.). Subsequently, we identified five upregulated miRNA via quantitative polymerase chain reaction and all of them were confirmed increased significantly in aged skin tissues compared with young control tissues. Among them, high expression of miR-302b-3p was verified in both aged skin tissues and senescence fibroblasts. Furthermore, miR-302b-3p mimic accelerated skin fibroblast senescence and suppressed the longevity-associated gene Sirtuin 1(Sirt1) expression, whereas miR-302b-3p inhibitor could delay skin fibroblast senescence and contribute Sirt1 expression. In addition, we demonstrated that c-Jun N-terminal kinase 2(JNK2) is a direct target of miR-302b-3p by a luciferase reporter assay. An inverse correlation was verified in fibroblasts between miR-302b-3p and JNK2. Most importantly, siRNA JNK2 confirmed that low expression of JNK2 could accelerate fibroblasts senescence. In conclusion, our results indicated that overexpressed miR-302b-3p plays an important biological role in accelerating skin aging process via directly targeting JNK2 gene.     

Analysis of the expression,function, and evolution of miR-27 isoforms and their responses in metabolic processes

This study aimed to discuss the potential roles of isomiRs of miR-27 family in metabolisms associated with disease via analyses of their evolution, expression, and function. miR-27b-3p was relatively highly expressed in liver cancer samples compared to miR-27a-3p and miR-27-5p loci. The diversity of isomiRs in miR-27-3p locus is similar to that of miRNAs among homologous genes. IsomiRs exhibited variable expression across different cancer tissue types, and some of them were abnormally expressed in ob/ob mice. Further experimental validation indicated that the protein expression of metabolism-related proteins, including PEPCK, G6Pase, FAS, and CPT1A, were significantly suppressed when canonical miR-27b was transfected into AML-12 cells. In contrast, the expression of these proteins was only slightly inhibited by isomiR-27b-1 or isomiR-27b-2 after transfection into AML-12 cells. These observations support that isomiRs exhibiting sequence divergence are functional regulatory molecules, and that they may contribute to biological processes via coordinated interactions in regulatory networks.     

MicroRNA Profiling Implies New Markers of Chemoresistance of Triple-Negative Breast Cancer

Objective

Triple-negative breast cancer (TNBC) patients with truly chemosensitive disease still represent a minority among all TNBC patients. The aim of the present study is to identify microRNAs (miRNAs) that correlate with TNBC chemoresistance.

Methods

In this study, we conducted miRNAs profile comparison between triple-negative breast cancer (TNBCs) and normal breast tissues by microRNA array. Quantitative real-time PCR (qRT-PCR) was utilized to confirm the specific deregulated miRNAs change trend. We used starBase 2.1 and GOrilla to predict the potential targets of the specific miRNAs. Cells viability and apoptosis assays were employed to determine the effect of alteration of the specific miRNAs in TNBC cells on the chemosensitivity.

Results

We identified 11 specific deregulated miRNAs, including 5 up-regulated miRNAs (miR-155-5p, miR-21-3p, miR-181a-5p, miR-181b-5p, and miR-183-5p) and 6 down-regulated miRNAs (miR-10b-5p, miR-451a, miR-125b-5p, miR-31-5p, miR-195-5p and miR-130a-3p). Thereafter, this result was confirmed by qRT-PCR. We predicted the potential targets of the candidate miRNAs and found that they are involved in cancer-associated pathways. For the first time, we found that miR-130a-3p and miR-451a were down-regulated in TNBC. 9 of the 11 specific deregulated miRNAs were found to be associated with chemoresistance. In vitro assays, we found that up-regulation of either miR-130a-3p or miR-451a in MDA-MB-231 cells significantly changed the cells sensitivity to doxorubicin. The results suggest that TNBC chemotherapy might be affected by a cluster of miRNAs.

Conclusion

The abnormal expression miRNAs in TNBC are mainly chemoresistance related. This might be part of reason that TNBC likely to evade from chemotherapy resulting in early relapse and high risk of death. To alter their expression status might be a potential therapeutic strategy to improve the outcome of chemotherapy for TNBC patients.     

Licochalcone A reverses NNK-induced ectopic miRNA expression to elicit in vitro and in vivo chemopreventive effects

BackgroundChemoprevention is the best cost-effective way regarding cancers. MicroRNAs (miRNAs) have been reported to be differentially expressed during the development of lung cancer. However, if lung cancer prevention can be achieved through modulating miRNAs expression so far remains unknown.PurposeTo discover ectopically expressed miRNAs in NNK-induced lung cancer and clarify whether Licochalcone A (lico A) can prevent NNK-induced lung cancer by modulating miRNA expression.Study design and methodsA/J mice were used to construct a lung cancer model by intraperitoneal injection with physiological saline NNK (100 mg/kg). Chemopreventive effects of lico A against lung cancer at 2 mg/kg and 20 mg/kg doses were evaluated in vivo. MicroRNA array and RT-qPCR were used to assess the expression levels of miRNAs. MLE-12 cells were treated with 0.1 mg/ml NNK, stimulating the ectopic expression pattern of miR-144-3p, miR-20a-5p, miR-29c-3p, let-7d-3p, and miR-328-3p. miR-144-3p mimics and inhibitors were used to manipulate miR-144-3p levels. The effects of lico A (10 μM) on cell cycle distribution, apoptosis, and the expression of CK19, RASA1, miR-144-3p, miR-20a-5p, miR-29c-3p, let-7d-3p, and miR-328-3p in NNK-treated MLE-12 cells were studied.ResultsThe expression levels of miR-144-3p, miR-20a-5p, and miR-29c-3p increased, while those of let-7d-3p and miR-328-3p decreased in both NNK-induced A/J mice and MLE-12 cells. Lico A could reverse the NNK-induced ectopic miRNA (miR-144-3p, miR-20a-5p, miR-29c-3p, let-7d-3p, and miR-328-3p) expression both in vivo and in vitro and elicit in vivo lung cancer chemopreventive effect against NNK. In MLE-12 cells, the overexpression of miR-144-3p elicited the same effect as NNK regarding the expression of lung cancer biomarker CK19; the silencing of miR-144-3p reversed the effect of NNK on cell cycle distribution and apoptosis. Lico A could reverse the effect of NNK on the expression of miR-144-3p, CK19, and RASA1 (predicted target of miR-144-3p).ConclusionThe present study suggests that miR-144-3p, miR-20a-5p, miR-29c-3p, let-7d-3p, and miR-328-3p were involved in the in vivo pathogenesis of NNK-induced lung cancer, and lico A could reverse the effect of NNK both in vivo and in vitro to elicit lung cancer chemopreventive effects through, at least partially, these five ectopically expressed miRNAs, especially miR-144-3p.     

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.     

Changes in microRNA abundance may regulate diapause in the flesh fly,Sarcophaga bullata

Diapause, an alternative developmental pathway characterized by changes in developmental timing and metabolism, is coordinated by molecular mechanisms that are not completely understood. MicroRNA (miRNA) mediated gene silencing is emerging as a key component of animal development and may have a significant role in initiating, maintaining, and terminating insect diapause. In the present study, we test this possibility by using high-throughput sequencing and qRT-PCR to discover diapause-related shifts in miRNA abundance in the flesh fly, Sarcophaga bullata. We identified ten evolutionarily conserved miRNAs that were differentially expressed in diapausing pupae compared to their nondiapausing counterparts. miR-289-5p and miR-1-3p were overexpressed in diapausing pupae and may be responsible for silencing expression of candidate genes during diapause. miR-9c-5p, miR-13b-3p, miR-31a-5p, miR-92b-3p, miR-275-3p, miR-276a-3p, miR-277-3p, and miR-305-5p were underexpressed in diapausing pupae and may contribute to increased expression of heat shock proteins and other factors necessary for the enhanced environmental stress-response that is a feature of diapause. In S. bullata, a maternal effect blocks the programming of diapause in progeny of females that have experienced pupal diapause, and in this study we report that several miRNAs, including miR-263a-5p, miR-100-5p, miR-125-5p, and let-7-5p were significantly overexpressed in such nondiapausing flies and may prevent entry into diapause. Together these miRNAs appear to be integral to the molecular processes that mediate entry into diapause.     

外周血miR-17-92基因簇对于早期胃癌的诊断及预测价值

目的:探讨外周血mi R-17-92簇对早期胃癌的诊断价值,为胃癌的早期诊断及治疗提供参考依据。方法:收集胃癌125例(Ⅰ期35例,Ⅱ期28例,Ⅲ期39例,Ⅳ期23例)和癌前病变24例(包括肠化生及上皮内瘤变),同时选择65例慢性胃炎作为对照组。采用实时荧光定量PCR技术(Real-time quantitative PCR,RT-qPCR)检测患者血清中的mi R-17-92基因簇的表达水平。通过受试者工作曲线(Receiver Operating Curve, ROC)及曲线下的面积(Area Under the Curve,AUC)评估mi R-17-92基因簇表达水平诊断早期胃癌的敏感性和特异性。结果:(1)慢性胃炎与癌前病变mi R-17-92基因簇表达比较无显著差异(P0.05);(2)早期胃癌及进展期胃癌mi R-17-5p表达明显高于慢性胃炎(P0.05),mi R-19a-3p、mi R-19b-3p、mi R-20a-5p和mi R-92a-3p表达则显著低于慢性胃炎及进展期胃癌(P0.05);(3)miR-17-5p诊断早期胃癌的曲线下面积较mi R-19a-3p、mi R-19b-3p、mi R-20a-5p、mi R-92a-3p及CEA更高;(4)miR-19a-3p、mi R-19b-3p、mi R-20a-5p、mi R-92a-3p高低表达组与在胃癌的浸润深度间有显著性差异(P0.05),mi R-19b-3p高低表达组在胃癌的临床分期间有显著性差异(P0.05);(5)miR-17-5p、mi R-19a-3p、mi R-19b-p、mi R-20a-5p、mi R-92a-3p诊断早期胃癌的阳性率较CEA、CA199高。结论:外周血mi R-17-92基因簇对于早期胃癌的诊断价值明显优于CEA和CA199,这可能为胃癌的早诊早治提供新的策略。     

Isoflurane and low-level carbon monoxide exposures increase expression of pro-survival miRNA in neonatal mouse heart

Anesthetics such as isoflurane are known to cause apoptosis in the developing mammalian brain. However, isoflurane may have protective effects on the heart via relieving ischemia and downregulating genes related to apoptosis. Ischemic preconditioning, e.g. through the use of low levels of carbon monoxide (CO), has promise in preventing ischemia-reperfusion injury and cell death. However, it is still unclear how it either triggers the stress response in neonatal hearts. For this reason, thirty-three microRNAs (miRNAs) known to be differentially expressed following anesthesia and/or ischemic or hypoxic heart damage were investigated in the hearts from neonatal mice exposed to isoflurane or low level of CO, using an air-exposed control group. Only miR-93-5p increased with isoflurane exposure, which may be associated with the suppression of cell death, autophagy, and inflammation. By contrast, twelve miRNAs were differentially expressed in the heart following CO treatment. Many miRNAs previously shown to be responsible for suppressing cell death, autophagy, and myocardial hypertrophy were upregulated (e.g., 125b-3p, 19-3p, and 21a-5p). Finally, some miRNAs (miR-103-3p, miR-1a-3p, miR-199a-1-5p) which have been implicated in regulating energy balance and cardiac contraction were also differentially expressed. Overall, this study demonstrated that CO-mediated miRNA regulation may promote ischemic preconditioning and cardioprotection based on the putative protective roles of the differentially expressed miRNAs explored herein and the consistency of these results with those that have shown positive effects of CO on heart viability following anesthesia and ischemia-reperfusion stress.Supplementary InformationThe online version contains supplementary material available at 10.1007/s12192-021-01199-0.     

miRNA调控恶性黑色素瘤细胞上皮-间充质转化的分子网络及机制

黑色素瘤是一种极易发生转移的恶性皮肤肿瘤,具有高度的致死性。上皮-间充质细胞转化(Epithelial-mesenchymal transition, EMT)在胚胎发育过程中起到非常重要的作用,同时在肿瘤的发生和恶化过程中也扮演着重要的角色。miRNA具有广谱的调节能力,对于肿瘤发生和EMT形成都能产生不同程度的影响。本文整合黑色素瘤细胞系转录组和miRNA组测序数据,在转录组数据中筛选得到参与肿瘤EMT过程的基因,通过Mirsystem软件预测并从miRNA组数据中筛选出与之负相关的11个miRNA,包括miR-130a-3p、miR-130b-3p、miR-125a-5p、miR-30a-3p、miR-195-5p、miR-345-5p、miR-509-3-5p、miR-374a-5p、miR-509-5p、miR-148a-3p和miR-330-3p。经过生物信息学分析miRNA靶基因富集的分子网络和信号途径,发现了两个与细胞发育和细胞间相互作用密切相关的网络,以及多个参与调控EMT过程的信号通路。对11个miRNA进行分子生物学验证,发现miR-195-5p、miR-130a-3p、miR-509-5p和miR-509-3-5p共4个可以调节重要肿瘤基因的miRNA。本研究运用mRNA和miRNA两种转录组的测序数据筛选EMT相关miRNA的方法,为肿瘤多组学数据整合分析提供了新的研究思路,并以期能为肿瘤精准基因组学的发展发挥重要的推进作用。