急性髓系白血病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诊断和治疗的分子生物学研究提供了新的依据。     

Investigation of key microRNAs associated with hepatocellular carcinoma using small RNA-seq data

To identify key microRNAs (miRNAs) associated with hepatocellular carcinoma (HCC) using small RNA-seq data. Small RNA-seq data for two HCC samples and two normal samples were downloaded from NCBI Gene Expression Omnibus. MiRNAs were identified through database search. Differentially expressed miRNAs were screened out with t test and their target genes were retrieved. Functional enrichment analysis was performed to uncover their biological functions. Regulatory networks and core metabolic networks were also constructed to present the global patterns. In addition, new miRNAs and their target genes were predicted. A total of 59 differentially expressed miRNAs were obtained, 12 up-regulated and 47 down-regulated. A total of 3,306 target genes were retrieved for eight miRNAs. Pathway enrichment analysis for the target genes showed that “pathways in cancer” and “MAPK signaling pathway” were significantly over-represented. Functional enrichment analysis found that “biological regulation” and “macromolecule modification” were significantly related to the target genes. Two regulatory networks were constructed for up- and down-regulated differentially expressed miRNAs with information from Ingenuity Pathway Analysis database. Two metabolic networks were also established based upon “pathways in cancer” and “MAPK signaling pathway”, consisting of miRNAs, target genes, compounds and others genes. Moreover, a number of new miRNAs and relevant target genes were predicted. Our study discloses a number of miRNAs as well as genes which may be involved in the development of HCC and these findings are beneficial in guiding future researches.     

Analysis of connection networks among miRNAs differentially expressed in early gastric cancer for disclosing some biological features of disease development

This paper first identified differentially expressed miRNAs associated with early gastric cancer and then respectively constructed relevant connection networks among the identified differentially expressed miRNAs that corresponded to early gastric cancer and control tissues. Twenty-three differentially expressed miRNAs were identified, 18 of which were different with the related results on the same data, and they provide great discriminatory power between patients and controls. There are not only conserved unchangeable sub-networks but also different sub-networks between the two connection networks. From the consistency and differences between two connection networks, we disclosed several new biological features that promote early gastric cancer development. This study shows 23 miRNAs that are early gastric cancer-specific and are worthy to do further experimental studies. The revealed biological features for early gastric cancer will provide new insights into improved understanding of the molecular mechanisms of this disease.     

EA15, MIR22, LINC00472 as diagnostic markers for diabetic kidney disease

This study aimed to investigate the molecular mechanisms of diabetic kidney disease (DKD) and to explore new potential therapeutic strategies and biomarkers for DKD. First we analyzed the differentially expressed changes between patients with DKD and the control group using the chip data in Gene Expression Omnibus (GEO) database. Then the gene chip was subjected to be annotated again, so as to screen long noncoding RNAs (lncRNAs) and study expression differences of these lncRNAs in DKD and controlled samples. At last, the function of the differential lncRNAs was analyzed. A total of 252 lncRNAs were identified, and 14 were differentially expressed. In addition, there were 1,629 differentially expressed messenger RNAs (mRNAs) genes, and proliferation and apoptosis adapter protein 15 (PEA15), MIR22, and long intergenic nonprotein coding RNA 472 ( LINC00472) were significantly differentially expressed in DKD samples. Through functional analysis of the encoding genes coexpressed by the three lncRNAs, we found these genes were mainly enriched in type 1 diabetes and autoimmune thyroid disease pathways, whereas in Gene Ontology (GO) function classification, they were also mainly enriched in the immune response, type I interferon signaling pathways, interferon-γ mediated signaling pathways, and so forth. To summary, we identified EA15, MIR22, and LINC00472 may serve as the potential diagnostic markers of DKD.     

Pathways related to PMA-differentiated THP1 human monocytic leukemia cells revealed by RNA-Seq

Previous analyses have reported that the human monocytic cell line THP1 can be differentiated into cells with macrophage-like characteristics by phorbol 12-myristate 13-acetate(PMA). However, little is known about the mechanism responsible for regulating this differentiation process. Here, we performed high-throughput RNA-Seq analysis to investigate the genes differently expressed in THP1 cells treated with and without PMA and examined those that may be responsible for the PMA-induced differentiation of monocytes into macrophages. We found 3,000 genes to be differentially expressed after PMA treatment. Gene ontology analysis revealed that genes related to cellular processes and regulation of biological processes were significantly enriched. KEGG analysis also demonstrated that the differentially expressed genes(DEGs) were significantly enriched in the PI3K/AKT signaling pathway and phagosome pathway. Importantly, we reveal an important role of the PI3K/AKT pathway in PMA-induced THP1 cell differentiation. The identified DEGs and pathways may facilitate further study of the detailed molecular mechanisms of THP1 differentiation. Thus, our results provide numerous potential therapeutic targets for modulation of the differentiation of this disease.     

高通量测序技术分析急性髓系白血病血浆外泌体微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的方法具有十分重要和深远的意义。     

Comparative analysis of neurological disorders focuses genome-wide search for autism genes

The behaviors of autism overlap with a diverse array of other neurological disorders, suggesting common molecular mechanisms. We conducted a large comparative analysis of the network of genes linked to autism with those of 432 other neurological diseases to circumscribe a multi-disorder subcomponent of autism. We leveraged the biological process and interaction properties of these multi-disorder autism genes to overcome the across-the-board multiple hypothesis corrections that a purely data-driven approach requires. Using prior knowledge of biological process, we identified 154 genes not previously linked to autism of which 42% were significantly differentially expressed in autistic individuals. Then, using prior knowledge from interaction networks of disorders related to autism, we uncovered 334 new genes that interact with published autism genes, of which 87% were significantly differentially regulated in autistic individuals. Our analysis provided a novel picture of autism from the perspective of related neurological disorders and suggested a model by which prior knowledge of interaction networks can inform and focus genome-scale studies of complex neurological disorders.     

RNA-seq用于获得共轭亚油酸对贵妃鸡肌内脂肪代谢调控的关键基因

通过转录物组测序获得在贵妃鸡基础日粮中添加共轭亚油酸(CLA)对肌内脂肪代谢的差异表达基因,经生物信息学分析获得相关的信号通路及可能发挥重要作用的候选基因,为CLA对肌内脂肪沉积的分子机制奠定基础。本研究选用55日龄健康的贵妃鸡为试验动物,在基础日粮中添加CLA 0%、1%和2%,预饲期1周,正饲期6周。屠宰采集胸肌组织进行转录物组测序,对测序数据进行差异表达分析,差异表达基因GO功能和差异表达基因KEGG通路富集分析,筛选出与胸肌脂类代谢相关的差异表达基因,利用qRT-PCR对差异表达基因进行验证。结果显示,共获得1 065个差异表达基因,其中上调基因703个,下调基因362个。GO富集结果显示,差异表达基因主要富集在生物过程的细胞过程、单一生物过程、生物调节和代谢过程。KEGG信号通路富集显示,差异表达基因显著富集在黏着斑、不饱和脂肪酸生物合成、脂肪酸生物合成和类固醇生物合成等信号通路中,发现11个主要与肌内脂肪代谢相关的候选基因,分别是FADS1、FADS2、ELOVL5、ACOX2、SLC27A1、FABP5、LPL、LOC107050163、ENSGALG00000030996、ENSGALG00000005043和ENSGALG00000048882。并随机选取6个基因进行qRT-PCR验证,其相对表达量变化趋势与测序结果一致。本研究筛选到CLA影响贵妃鸡胸肌脂类代谢相关的差异表达基因,并对11个主要参与脂肪代谢相关的基因进行分析,为揭示CLA调控肌内脂肪沉积的分子机制奠定基础。     

茶足柄瘤蚜茧蜂蛹滞育过程中胰岛素信号通路及其相关途径的初探

本实验通过探索胰岛素信号通路及其相关途径对茶足柄瘤蚜茧蜂蛹滞育的影响,从而方便寻找胰岛素替代物,为害虫防治提供新思路。利用RNA-Seq,对滞育组与非滞育组的茶足柄瘤蚜茧蜂进行转录组测序,结合生物信息学方法对转录组中胰岛素信号通路及其相关途径的差异表达基因进行了分析。与胰岛素信号通路相关差异表达基因共31个,重点分析的PI3K-Akt, FoxO, MAPK三条途径,差异表达基因分别为55, 21和28个。这些滞育关联基因呈现不同程度的上调或下调表达,发现Sos, FASN, TSC1, PRKAB等基因与茶足柄瘤蚜茧蜂滞育密切相关,共同影响茶足柄瘤蚜茧蜂的滞育。胰岛素信号通路及其相关途径对茶足柄瘤蚜茧蜂的滞育起着非常重要的作用,主要体现在影响虫体能量代谢、脂质积累、细胞增殖等方面。     

Differential gene expression profiling of human epidermal growth factor receptor 2-overexpressing mammary tumor

Human epidermal growth factor receptor 2 (HER2) is highly expressed in approximately 30% of breast cancer patients, and substantial evidence supports the relationship between HER2 overexpression and poor overall survival. However, the biological function of HER2 signal transduction pathways is not entirely clear. To investigate gene activation within the pathways, we screened differentially expressed genes in HER2-positive mouse mammary tumor using two-directional suppression subtractive hybridization combined with reverse dot-blotting analysis. Forty genes and expressed sequence tags related to transduction, cell proliferation/growth/apoptosis and secreted/extracellular matrix proteins were differentially expressed in HER2-positive mammary tumor tissue. Among these, 19 were already reported to be differentially expressed in mammary tumor, 11 were first identified to be differentially expressed in mammary tumor in this study but were already reported in other tumors, and 10 correlated with other cancers. These genes can facilitate the understanding of the role of HER2 signaling in breast cancer.     

Modulation of innate immune-related pathways in nicotine-treated SH-SY5Y cells

Although nicotine has a broad impact on both the central and peripheral nervous systems, the molecular mechanisms remain largely unknown, especially at the signaling pathway level. To investigate that aspect, we employed both conventional molecular techniques, such as quantitative real-time PCR and Western blotting analysis, and high-throughput microarray approach to identify the genes and signaling pathways that are modulated by nicotine. We found 14 pathways significantly altered in SH-SY5Y neuroblastoma cells. Of these, the Toll-like receptor pathway (TLR; p?=?2.57?×?10(-4)) is one of the most important innate immune pathways. The death receptor pathway (DR; p?=?8.71?×?10(-4)), whose transducers coordinate TLR signals and help conduct the host immune response to infection, was also significantly changed by nicotine. Furthermore, we found that several downstream pathways of TLR and DR signaling, such as PI3K/AKT signaling (p?=?9.55?×?10(-6)), p38 signaling (p?=?2.40?×?10(-6)), and ERK signaling (p?=?1.70?×?10(-4)), were also significantly modulated by nicotine. Interestingly, most of the differentially expressed genes in these pathways leading to nuclear factor κB (NF-κB) activation and those important inhibitors of pathways leading to apoptosis, including FLIP and Bcl-2, were up-regulated by nicotine. Taken together, our findings demonstrate that nicotine can regulate multiple innate immune-related pathways, and our data thus provide new clues to the molecular mechanisms underlying nicotine's regulatory effects on neurons.     

Genome-Wide Functional Analysis of Cotton (Gossypium hirsutum) in Response to Drought

Cotton is one of the most important crops for its natural textile fibers in the world. However, it often suffered from drought stress during its growth and development, resulting in a drastic reduction in cotton productivity. Therefore, study on molecular mechanism of cotton drought-tolerance is very important for increasing cotton production. To investigate molecular mechanism of cotton drought-resistance, we employed RNA-Seq technology to identify differentially expressed genes in the leaves of two different cultivars (drought-resistant cultivar J-13 and drought-sensitive cultivar Lu-6) of cotton. The results indicated that there are about 13.38% to 18.75% of all the unigenes differentially expressed in drought-resistant sample and drought-sensitive control, and the number of differentially expressed genes was increased along with prolonged drought treatment. DEG (differentially expression gene) analysis showed that the normal biophysical profiles of cotton (cultivar J-13) were affected by drought stress, and some cellular metabolic processes (including photosynthesis) were inhibited in cotton under drought conditions. Furthermore, the experimental data revealed that there were significant differences in expression levels of the genes related to abscisic acid signaling, ethylene signaling and jasmonic acid signaling pathways between drought-resistant cultivar J-13 and drought-sensitive cultivar Lu-6, implying that these signaling pathways may participate in cotton response and tolerance to drought stress.     

基因芯片技术筛选转SlNAC1基因拟南芥差异表达基因

研究已表明植物特有的一些NAC(NAM,ATAF1/2,CUC2)转录因子可提高植物抗逆性,利用基因芯片技术筛选转SlNAC1基因拟南芥与野生型拟南芥间差异表达基因,能够为研究转基因拟南芥非生物胁迫抗性相关基因提供依据。结果显示,在转SlNAC1基因拟南芥43 604个基因中有3 046个差异表达2倍以上的基因。对差异表达5倍以上基因经过GO富集度统计学分析表明,细胞组分相关基因占33.05%;分子功能相关基因占33.95%;生物学过程相关基因占33.00%。对差异表达2倍以上基因进行KEGG信号通路分析,结果表明有2 431个基因涉及到88个不同的信号通路。通过筛选获得转基因拟南芥非生物胁迫抗性相关候选基因,为后续研究NAC转录因子的下游基因及其调控网络的构建提供方向和理论支撑。