急性髓系白血病mRNA与非编码RNA差异表达谱及CeRNA调控网络分析

利用GEO数据库(gene expression omnibus database)通过生物信息学分析方法探讨急性髓系白血病(acute myelogenous leukemia,AML)的发病机制。检索GEO数据库中AML相关芯片数据集GSE142698、GSE142699和GSE96535。利用GEO2R分析得到差异mRNAs、miRNAs以及差异lncRNAs。利用在线生物信息学分析工具DAVID对差异mRNAs进行GO富集分析和KEGG通路分析。利用miRWalk数据库预测AML相关miRNAs的靶向mRNAs,利用Spongescan数据库预测AML相关miRNAs的靶向lncRNAs,构建lncRNA-miRNA-mRNA竞争性内源RNA （competing endogenous RNA,ceRNA）调控网络。共筛选出29个显著差异mRNAs、70个显著差异miRNAs和20 005个显著差异lncRNAs。GO富集分析和KEGG通路分析显示,差异表达基因主要涉及蛋白磷酸化、细胞分裂、细胞增殖的负调控、基因表达的正向调节、周期蛋白依赖的丝氨酸/苏氨酸激酶活性的调节等生物过程以及细胞周期、细胞衰老、癌症通路、PI3K-Akt通路等信号通路。将miRWalk数据库预测的靶向mRNAs与差异mRNAs取交集,Spongescan数据库预测的靶向lncRNAs与差异lncRNAs取交集,分别确定了25个mRNAs、6个lncRNAs参与AML相关ceRNA调控网络的构建。结果表明,lncRNAs可能作为关键的ceRNA,通过调控miRNA和相关靶基因参与AML的发生与发展,研究结果为AML诊断和治疗的分子生物学研究提供了新的依据。     

桑树桑黄菌丝体和子实体基因差异表达分析

通过对桑树桑黄Sanghuangporus sanghuang菌丝体和子实体2个不同生长阶段的转录组进行分析,为研究桑黄子实体生长发育相关机制奠定基础。采用Illumina测序技术,对桑树桑黄菌株S23菌丝体和子实体2个不同生长发育阶段进行了全转录组测序。将转录组测序reads比对到参考序列上,菌丝体测序样本的reads比对率为82.89%;子实体测序样本的reads比对率为83%。基因差异表达分析显示,与菌丝体相比,子实体中显著上调表达基因为2 898个,显著下调表达基因为1 965个。经过Blast nr比对发现,桑黄菌在子实体阶段表达量上升的基因主要与各种氧化酶活性、疏水蛋白等相关;表达量下降的基因主要与糖类、氨基酸结合、运输等相关。基因本体（gene ontology,GO）富集分析表明,菌丝体及子实体两个阶段与跨膜转运相关的差异表达基因富集明显。代谢通路（pathway）富集分析表明,类固醇生物合成、精氨酸生物合成、丝裂原活化蛋白激酶（mitogen-activated protein kinase,MAPK）信号通路等差异基因富集明显。     

肾透明细胞癌预后相关miRNAs的生物信息学分析

目的寻找可作为肾透明细胞癌（ccRCC）生物标志物的miRNA,以及ccRCC与正常组织间miRNA差异表达情况。 方法利用TCGA数据库下载ccRCC中miRNA表达数据,分析肿瘤与正常组织间差异表达miRNA。使用Kaplan-Meier曲线对患者进行生存分析,筛选出表达情况与临床预后相关的miRNA。通过生物信息学对miRNA的靶基因进行预测,然后运用FunRich软件和ClueGO对靶基因进行GO和KEGG富集分析。 结果通过TCGA数据库分析发现,ccRCC较正常组织差异表达miRNA共54个,其中上调33个,下调21个。通过生存分析发现hsa-miR-21和hsa-miR-155与患者预后相关,P≤0.05。进一步通过Perl软件在Targetscan、miRDB、miRTarBase、miRPath这四个数据库中预测miRNA靶基因并将结果取交集,共发现129个靶基因。GO和KEGG分析结果表明,这些靶基因主要与转录因子活性、信号转导以及FoxO、TNF等信号通路密切相关。 结论通过生物信息学分析发现了ccRCC与正常组织的差异表达miRNA;其中hsa-miR-21和hsa-miR-155与患者总体生存率相关,并通过调控靶基因参与相关的信号通路进而影响ccRCC的发生发展进程,提示hsa-miR-21和hsa-miR-155可能是ccRCC潜在的生物标志物。     

Expression profiles of microRNAs in oxidized low-density lipoprotein-stimulated RAW 264.7 cells

Macrophage-derived foam cells were one of the hallmarks of atherosclerosis, and microRNAs played an important role in the formation of foam cells. In order to explore the roles of miRNA in the formation of foam cells, we investigated miRNA expression profiles in foam cells through high-throughput sequencing technology. A total of 84 miRNAs were differentially expressed between RAW 264.7 macrophages and foam cells induced by ox-LDL. Thirty miRNAs were upregulated and 54 miRNAs were downregulated. GO terms and KEGG pathways analysis revealed that the target genes of most of DE miRNAs were mainly enriched in “cell differentiation,” “endocytosis,” “MAPK signaling pathway,” and “FoxO signaling pathway.” The target genes of some DE miRNAs were enriched in “Insulin signaling pathway,” “Hippo signaling pathway,” “TNF signaling pathway,” “NF-kappa B signaling pathway,” and “cell death.” Using bioinformatics analyses and dual-luciferase reporter assays, we found that miR-28a-5p and miR-30c-1-3p directly inhibited LRAD3 and LOX-1 mRNA expression through targeting the 3’UTR of LRAD3 and LOX-1 mRNA, respectively. Our study indicates that miRNAs are extensively involved in the formation of foam cells, and provides a valuable resource for further study the role of miRNAs in atherosclerosis.     

胃癌多药耐药相关microRNA的筛选和鉴定

目的:研究胃癌多药耐药相关microRNA并对其进行鉴定、靶基因预测和预测靶基因的生物信息学分析。方法:运用microRNA芯片对胃癌多药耐药细胞SGC7901/ADR和其亲本细胞SGC7901进行microRNA表达谱分析;采用实时定量PCR的方法对差异表达的miRNA进行验证;再运用生物信息学方法对差异表达的miRNA进行靶基因预测;再对预测的靶基因进行GO和KEGG通路分析。结果:与SGC7901相比SGC7901/ADR表达上调超过2倍的miRNA有6个,表达下调超过2倍的有11个。实时定量PCR对共同差异表达的microRNA进行验证显示与芯片结果的一致性。对这17个差异表达的miRNA进行靶基因预测,再对预测得到的靶基因进行GO和KEGG通路分析显示预测的靶基因参与了肿瘤相关通路、MAPK通路、Focal Adhesion通路等。结论:我们初步筛选得到了胃癌多药耐药相关miRNA并对其进行了生物信息学分析,为进一步地探索miRNA在胃癌多药耐药中的作用及其分子机制奠定了基础。     

Comparison of microRNAs in adipose and muscle tissue from seven indigenous Chinese breeds and Yorkshire pigs

Elucidation of the pig microRNAome is essential for interpreting functional elements of the genome and understanding the genetic architecture of complex traits. Here, we extracted small RNAs from skeletal muscle and adipose tissue, and we compared their expression levels between one Western breed (Yorkshire) and seven indigenous Chinese breeds. We detected the expression of 172 known porcine microRNAs (miRNAs) and 181 novel miRNAs. Differential expression analysis found 92 and 12 differentially expressed miRNAs in adipose and muscle tissue respectively. We found that different Chinese breeds shared common directional miRNA expression changes compared to Yorkshire pigs. Some miRNAs differentially expressed across multiple Chinese breeds, including ssc‐miR‐129‐5p, ssc‐miR‐30 and ssc‐miR‐150, are involved in adipose tissue function. Functional enrichment analysis revealed that the target genes of the differentially expressed miRNAs are associated mainly with signaling pathways rather than metabolic and biosynthetic processes. The miRNA–target gene and miRNA–phenotypic traits networks identified many hub miRNAs that regulate a large number of target genes or phenotypic traits. Specifically, we found that intramuscular fat content is regulated by the greatest number of miRNAs in muscle tissue. This study provides valuable new candidate miRNAs that will aid in the improvement of meat quality and production.     

A dictionary on microRNAs and their putative target pathways

While in the last decade mRNA expression profiling was among the most popular research areas, over the past years the study of non-coding RNAs, especially microRNAs (miRNAs), has gained increasing interest. For almost 900 known human miRNAs hundreds of pretended targets are known. However, there is only limited knowledge about putative systemic effects of changes in the expression of miRNAs and their regulatory influence. We determined for each known miRNA the biochemical pathways in the KEGG and TRANSPATH database and the Gene Ontology categories that are enriched with respect to its target genes. We refer to these pathways and categories as target pathways of the corresponding miRNA. Investigating target pathways of miRNAs we found a strong relation to disease-related regulatory pathways, including mitogen-activated protein kinase (MAPK) signaling cascade, Transforming growth factor (TGF)-beta signaling pathway or the p53 network. Performing a sophisticated analysis of differentially expressed genes of 13 cancer data sets extracted from gene expression omnibus (GEO) showed that targets of specific miRNAs were significantly deregulated in these sets. The respective miRNA target analysis is also a novel part of our gene set analysis pipeline GeneTrail. Our study represents a comprehensive theoretical analysis of the relationship between miRNAs and their predicted target pathways. Our target pathways analysis provides a ‘miRNA-target pathway’ dictionary, which enables researchers to identify target pathways of differentially regulated miRNAs.     

The RNA expression signature of the HepG2 cell line as determined by the integrated analysis of miRNA and mRNA expression profiles

Understanding miRNAs' regulatory networks and target genes could facilitate the development of therapies for human diseases such as cancer. Although much useful gene expression profiling data for tumor cell lines is available, microarray data for miRNAs and mRNAs in the human HepG2 cell line have only been compared with that of other cell lines separately. The relationship between miRNAs and mRNAs in integrated expression profiles for HepG2 cells is still unknown. To explore the miRNA–mRNA correlations in hepatocellular carcinoma (HCC) cells, we performed miRNA and mRNA expression profiling in HepG2 cells and normal liver HL-7702 cells at the genome scale using next-generation sequencing technology. We identified 193 miRNAs that are differentially expressed in these two cell lines. Of these, 89 miRNAs were down-regulated in HepG2 cells compared with HL-7702 cells, while 104 miRNAs were up-regulated. We also observed 3035 mRNAs that are significantly dys-regulated in HepG2 cells. We then performed an integrated analysis of the expression data for differentially expressed miRNAs and mRNAs and found several miRNA–mRNA pairs that are significantly correlated in HepG2 cells. Further analysis suggested that these differentially expressed genes were enriched in four tumorigenesis-related signaling pathways, namely, ErbB, JAK–STAT, mTOR, and WNT, which until now had not been fully reported. Our results could be helpful in understanding the mechanisms of HCC occurrence and development.     

Infection with street strain rabies virus induces modulation of the microRNA profile of the mouse brain

ABSTRACT: BACKGROUND: Rabies virus (RABV) causes a fatal infection of the central nervous systems (CNS) of warm-blooded animals. Once the clinical symptoms develop, rabies is almost invariably fatal. The mechanism of RABV pathogenesis remains poorly understood. Recent studies have shown that microRNA (miRNA) plays an important role in the pathogenesis of viral infections. Our recent findings have revealed that infection with laboratory-fixed rabies virus strain can induce modulation of the microRNA profile of mouse brains. However, no previous report has evaluated the miRNA expression profile of mouse brains infected with RABV street strain. RESULTS: The results of microarray analysis show that miRNA expression becomes modulated in the brains of mice infected with street RABV. Quantitative real-time PCR assay of the differentially expressed miRNAs confirmed the results of microarray assay. Functional analysis showed the differentially expressed miRNAs to be involved in many immune-related signaling pathways, such as the Jak-STAT signaling pathway, the MAPK signaling pathway, cytokine-cytokine receptor interactions, and Fc gamma R-mediated phagocytosis. The predicted expression levels of the target genes of these modulated miRNAs were found to be correlated with gene expression as measured by DNA microarray and qRT-PCR. CONCLUSION: RABV causes significant changes in the miRNA expression profiles of infected mouse brains. Predicted target genes of the differentially expression miRNAs are associated with host immune response, which may provide important information for investigation of RABV pathogenesis and therapeutic method.     

家蚕5龄幼虫精巢和卵巢microRNA芯片及转录组比较分析

【目的】本研究旨在通过对家蚕Bombyx mori 5龄幼虫精巢和卵巢组织微小RNA (microRNA, miRNA)基因芯片及转录组进行分析,找到参与家蚕性腺发育相关的miRNA分子及可能的靶基因。【方法】采用新一代高通量测序平台对家蚕5龄幼虫精巢和卵巢(分别定义为Test和Control)进行miRNA基因芯片检测及转录组测序分析,根据P＜0.05且log₂(fold change, FC)≥2的标准,通过比较筛选出Test vs Control的差异表达miRNA;根据q≤0.05且|log₂(fold change)|≥1的标准,通过比较筛选出Test vs Control的差异表达基因 (differentially expressed genes, DEGs);随机选取8个上调和12个下调差异表达miRNA,对其表达及其预测的5个靶基因进行qRT-PCR验证;对DEGs以及差异表达miRNA的靶基因进行KEGG通路富集分析。【结果】从精巢和卵巢样本中(Test vs Control)分别鉴定出68个差异表达miRNA和3 991个DEGs,其中上调和下调miRNA分别为36和32个,上调和下调DEGs分别为2 033和1 958个。差异表达miRNA的qRT PCR验证结果均与芯片数据一致。KEGG通路富集分析结果显示DEGs在新陈代谢及核糖体的信号通路显著富集。对差异表达miRNA在DEGs中的可能靶基因进行预测,结果找到了4组表达趋势相反的miRNA与靶基因：分别是bmo-miR-2774a与LOC101745556;bmo-miR-92b与LOC101735954以及bmo-miR-3266与LOC733130和LOC778467;1组表达趋势一致的miRNA与靶基因：bmo-miR-3321与LOC101744895。5个靶基因的qRT-PCR验证结果与转录组测序结果一致。【结论】本研究获得了家蚕5龄幼虫精巢和卵巢转录组及miRNA芯片数据,筛选并验证了4组差异表达和1组一致表达miRNA及潜在靶基因,为探究家蚕精巢和卵巢发育差异奠定了基础。     

miRNA Profiling Reveals Dysregulation of RET and RET-Regulating Pathways in Hirschsprung's Disease

Hirschsprung’s disease (HSCR), the most common congenital malformation of the gut, is regulated by multiple signal transduction pathways. Several components of these pathways are important targets for microRNAs (miRNAs). Multiple miRNAs have been associated with the pathophysiology of HSCR, and serum miRNAs profiles of HSCR patients have been reported, but miRNA expression in HSCR colon tissue is almost completely unexplored. Using microarray technology, we screened colon tissue to detect miRNAs whose expression profiles were altered in HSCR and identify targets of differentially expressed miRNAs. Following filtering of low-intensity signals, data normalization, and volcano plot filtering, we identified 168 differentially expressed miRNAs (104 up-regulated and 64 down-regulated). Fifty of these mRNAs represent major targets of dysegulated miRNAs and may thus important roles in the pathophysiology of HSCR. Pathway analysis revealed that 7 of the miRNA targets encode proteins involved in regulation of cell proliferation and migration via RET and related signaling pathways (MAPK and PI3K/AKT). Our results identify miRNAs that play key roles in the pathophysiology of the complex multi-factorial disease HSCR.     

n-3多不饱和脂肪酸调节饮食诱导肥胖大鼠miRNA表达变化(英文)

目的:研究n-3多不饱和脂肪酸(polyunsaturated fatty acids,PUFA)饮食对饮食诱导肥胖大鼠的miR NA表达影响。方法:将10只饮食诱导肥胖(diet induced obese,DIO)大鼠随机分成两组:n-3PUFA添加组和安慰剂添加组(对照组);每周记录两组老鼠的体重、体长和进食量。对外周血miR NA的表达并进行分析和预测。结果:两组老鼠Lee指数有统计学差异(P0.05);与对照组相比,在n-3组的外周血单核细胞中,29个miR NA上调,31个下调;其中rno-miR-200和rno-miR-211的表达量上调,rno-miR-29b和rno-miR-92b的表达量下调,其靶基因预测结果与神经营养因子,脂肪细胞因子,趋化因子和胰岛素信号通路有关。结论:n-3PUFA能够调节DIO大鼠的miR NA水平,其中有些与脂肪代谢相关。     

高通量测序技术分析急性髓系白血病血浆外泌体微RNA表达谱差异

急性髓系白血病（acute myeloid leukemia,AML）是一类造血干细胞的恶性克隆性疾病,目前的诊断方法不利于疾病的早期发现,且诊断结果重复性较差。已有大量研究显示,细胞外microRNA（miRNA或miR）富集在外泌体（exosome）中,且受其表面膜的保护而具有很好的稳定性,是理想的分子标志物。目前,多种实体肿瘤均已检测到肿瘤特异性外泌体miRNA(exosomal miRNA)。然而,在AML患者中未见此外泌体miRNA报道。本研究探讨急性髓系白血病血浆外泌体miRNA表达谱差异及新miRNA序列。采用solexa高通量测序技术对7例AML患者（AML组）及7例健康对照者（对照组）血浆外泌体miRNAs进行测序,利用Mireap预测软件进行新miRNAs分析,通过edger差异分析软件筛选组间差异miRNA,获得211个已知的差异表达miRNAs以及2个新miRNAs,选择4个差异表达的miRNAs：miR-155-5p、miR-335-5p、miR-451a及xxx-m0038 5p（新miRNA）,在两组（各23例）的血浆外泌体样本中,进行实时荧光定量PCR（qRT-PCR）验证,验证结果与测序结果一致。对差异表达的外泌体miRNA进行靶基因预测及其GO（Gene Ontology）和信号通路富集分析,发现靶基因聚集的生物学功能多数参与生物进展过程的调控。靶基因主要富集在FoxO、MAPK、Hippo信号通路以及HTLV-I感染等。结果显示,AML患者与健康对照者的血浆外泌体miRNA存在着差异性表达。差异性表达的miRNA特异性很高,对进一步阐明AML白血病发生与发展的分子机制、研发新的无创诊断方法、新的诊断标记物和有效治疗AML的方法具有十分重要和深远的意义。     

Computational methods for analysis of cellular functions and pathways collectively targeted by differentially expressed microRNA.

This report presents computational methods of analysis of cellular processes, functions, and pathways affected by differentially expressed microRNA, a statistical basis of the gene enrichment analysis method, a modification of enrichment analysis method accounting for combinatorial targeting of Gene Ontology categories by multiple miRNAs and examples of the global functional profiling of predicted targets of differentially expressed miRNAs in cancer. We have also summarized an application of Ingenuity Pathway Analysis tools for in depth analysis of microRNA target sets that may be useful for the biological interpretation of microRNA profiling data. To illustrate the utility of these methods, we report the main results of our recent computational analysis of five published datasets of aberrantly expressed microRNAs in five human cancers (pancreatic cancer, breast cancer, colon cancer, lung cancer, and lymphoma). Using a combinatorial target prediction algorithm and statistical enrichment analysis, we have determined Gene Ontology categories as well as biological functions, disease categories, toxicological categories, and signaling pathways that are: targeted by multiple microRNAs; statistically significantly enriched with target genes; and known to be affected in specific cancers. Our recent computational analysis of predicted targets of co-expressed miRNAs in five human cancers suggests that co-expressed miRNAs provide systemic compensatory response to the abnormal phenotypic changes in cancer cells by targeting a broad range of functional categories and signaling pathways reportedly affected in a particular cancer.     

参与调控意大利蜜蜂工蜂中肠基因表达的东方蜜蜂微孢子虫miRNA的组学解析及其调控网络

[目的]本研究结合前期已获得的miRNA和mRNA组学数据对东方蜜蜂微孢子虫Nosema ceranae 的差异表达 miRNA(differentially expressed miRNA,DEmiRNA)靶向意大利蜜蜂 Apis mellifera ligustica 工蜂中肠的 mRNA 和差异表达 mRNA(d...