金针菇单、双核菌丝差异表达基因分析

对金针菇Flammulina velutipes单核菌丝W23的菌丝体以及与L11质配后的双核菌丝H1123菌丝体进行了转录组测序,以本实验室已获得的W23基因组为参考基因组研究两样本间差异基因,并对这些差异基因进行了GO功能和Pathway显著性富集分析。差异基因分析显示,两个样本中共有显著性差异表达的基因3 504个,其中在双核菌丝中上调、下调的基因数分别为2 151和1 353个。研究发现差异表达基因含有很多的转录因子基因、蛋白激酶以及WD40 repeat-like蛋白。Gene Ontology（GO）功能分析结果表明,extracellular region和membrane-enclosed lumen条目下的差异基因全部为上调表达,而envelope下的差异基因全部为下调表达,以利于双核菌丝分裂时锁状联合的形成而便于核的迁移。Pathway 功能富集分析结果表明,脂肪酸、氨基酸以及大部分糖类合成相关基因具有比较活跃的上调表达。说明双核菌丝主要进行营养物质的富集,为下一步在合适条件下分化成原基,进入生殖生长阶段储备物质基础。     

Silencing SOCS3 could inhibit TNF-α induced apoptosis in 3T3-L1 and mouse preadipocytes

Tumor necrosis factor-alpha (TNF-α) is a pro-inflammatory cytokine involved in the apoptosis of many types of cells. In this study we demonstrated the effect of (suppressor of cytokine signalling-3) SOCS3 siRNA on TNF-α induced apoptosis in 3T3-L1 preadipocytes and mouse preadipocytes. 3T3-L1 preadipocytes and mouse preadipocytes were transfected with SOCS3 siRNA, and then the cells were treated with TNF-α at 100 ng/mL for 24 h. We used fluorescence microscope to observe morphological changes during apoptosis after Hoechst 33258 and PI staining. Quantitative PCR and Western blotting were used to measure the expression of apoptosis-associated gene c-myc, survivin, mcl-1, bcl-2, bax, NF-κB, and the key genes of the JAK/STAT3 pathway including SOCS1, SOCS2, JAK2, STAT3. Compared with control group, the number of cells apoptosis was decreased remarkably in SOCS3 siRNA group (P < 0.01). The expression of apoptotic suppressor genes c-myc, survivin, mcl-1, bcl-2 and NF-κB were up-regulated markedly (P < 0.01); in contrast, apoptotic gene bax was down-regulated (P < 0.05). Western blotting showed that the protein expressions of bcl-2 and NF-κB were increased remarkably (P < 0.01), while the protein expression of bax was decreased remarkably (P < 0.05). The expression of the JAK/STAT3 pathway key gene SOCS1 mRNA was down-regulated markedly (P < 0.05), but the key protein p-STAT3 was up-regulated (P < 0.05). Taken together, our data established that silenced SOCS3 can regulate the expression of apoptosis-associated genes via the JAK/STAT3 pathway, and effectively inhibit TNF-α induced apoptosis in 3T3-L1 preadipocytes and mouse preadipocytes.     

萱草叶片响应低温胁迫的转录组分析

低温胁迫是萱草（Hemerocallis fulva）生长过程中经常会遭遇的一种非生物胁迫。比较了萱草叶片在低温处理（10、5、0 ℃）下转录组与对照（15 ℃）数据的差异,共筛选出差异表达基因2 457个,其中上调基因1 253个,下调基因1 204个。差异表达基因主要富集在细胞过程、代谢过程和催化活性等49个GO过程,代谢途径、次生代谢产物的生物合成、植物激素信号转导等42条KEGG代谢途径中。其中参与植物激素信号转导通路的差异表达基因发生了不同程度的变化,GH3.10基因上调至对照组的13.624倍,IAA1基因下调0.120倍;参与可溶性糖合成通路的差异基因发生了0.076~28.114倍不同程度的变化。随后对3个低温处理组共有的29个差异表达基因进行热图和网络调控分析,基于基因在网络调控中的位置,对ABCF5、OFPs和SWEETs等基因在冷应答的作用进行了分析。本研究结果为进一步挖掘萱草低温响应的关键基因及耐寒萱草种质开发、分子育种提供了一定的理论支撑。     

Gene expression profiles and pathway enrichment analysis to identification of differentially expressed gene and signaling pathways in epithelial ovarian cancer based on high-throughput RNA-seq data

Epithelial ovarian cancer (EOC) can be considered as a stressful and challenging disease among all women in the world, which has been associated with a poor prognosis and its molecular pathogenesis has remained unclear. In recent years, RNA Sequencing (RNA-seq) has become a functional and amazing technology for profiling gene expression. In the present study, RNA-seq raw data from Sequence Read Archive (SRA) of six tumor and normal ovarian sample was extracted, and then analysis and statistical interpretation was done with Linux and R Packages from the open-source Bioconductor. Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were applied for the identification of key genes and pathways involved in EOC. We identified 1091 Differential Expression Genes (DEGs) which have been reported in various studies of ovarian cancer as well as other types of cancer. Among them, 333 genes were up-regulated and 273 genes were down-regulated. In addition, Differentially Expressed Genes (DEGs) including RPL41, ALDH3A2, ERBB2, MIEN1, RBM25, ATF4, UPF2, DDIT3, HOXB8 and IL17D as well as Ribosome and Glycolysis/Gluconeogenesis pathway have had the potentiality to be used as targets for EOC diagnosis and treatment. In this study, unlike that of any other studies on various cancers, ALDH3A2 was most down-regulated gene in most KEGG pathways, and ATF4 was most up-regulated gene in leucine zipper domain binding term. In the other hand, RPL41 as a regulatory of cellular ATF4 level was up-regulated in many term and pathways and augmentation of ATF4 could justify the increase of RPL41 in the EOC. Pivotal pathways and significant genes, which were identified in the present study, can be used for adaptation of different EOC study. However, further molecular biological experiments and computational processes are required to confirm the function of the identified genes associated with EOC.     

Gene Expression Profile of Bombyx mori Hemocyte under the Stress of Destruxin A

Destruxin A (DA) is a cyclo-peptidic mycotoxin from the entomopathogenic fungus Metarhizium anisopliae. To uncover potential genes associated with its molecular mechanisms, a digital gene expression (DGE) profiling analysis was used to compare differentially expressed genes in the hemocytes of silkworm larvae treated with DA. Ten DGE libraries were constructed, sequenced, and assembled, and the unigenes with least 2.0-fold difference were further analyzed. The numbers of up-regulated genes were 10, 20, 18, 74 and 8, as well as the numbers of down-regulated genes were 0, 1, 8, 13 and 3 at 1, 4, 8, 12 and 24 h post treatment, respectively. Totally, the expression of 132 genes were significantly changed, among them, 1, 3 and 12 genes were continually up-regulated at 4, 3 and 2 different time points, respectively, while 1 gene was either up or down-regulated continually at 2 different time points. Furthermore, 68 genes were assigned to one or multiple gene ontology (GO) terms and 89 genes were assigned to specific Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology. In-depth analysis identified that these genes putatively involved in insecticide resistance, cell apoptosis, and innate immune defense. Finally, twenty differentially expressed genes were randomly chosen and validated by quantitative real-time PCR (qRT-PCR). Our studies provide insights into the toxic effect of this microbial insecticide on silkworm''s hemocytes, and are helpful to better understanding of the molecular mechanisms of DA as a biological insecticide.     

Identification of potential biomarkers associated with immune infiltration in the esophageal carcinoma tumor microenvironment

Tumor immune cell infiltration was significantly correlated with the progression and the effect of immunotherapy in cancers including esophageal carcinoma (ESCA). However, no biomarkers were identified which were associated with immune infiltration in ESCA. In the present study, a total of 128 common differentially expressed genes (DEGs) were identified between esophageal squamous cell carcinomas (ESCC) and esophageal adenocarcinomas (EAC). The results of gene ontology (GO) enrichment and Reactome pathway analysis displayed that the up-regulated DEGs were mainly involved in the regulation of extracellular matrix (ECM), while the down-regulated DEGs were mainly involved in the regulation of cornification and keratinocyte differentiation. The most significant module of up-regulated DEGs was selected by Molecular Complex Detection (MCODE). Top ten similar genes of COL1A2 were explored, then validation and the prognostic analysis of these genes displayed that COL1A2, COL1A1, COL3A1, ZNF469 and Periostin (POSTN) had the prognostic value which were up-regulated in ESCA. The expressions of COL1A2 and its four similar genes were mainly correlated with infiltrating levels of macrophages and dendritic cells (DCs) and showed strong correlations with diverse immune marker sets in ESCA. To summarize, COL1A2 and its four similar genes were identified as the potential biomarkers associated with immune infiltration in ESCA. These genes might be applied to immunotherapy for ESCA.     

Microarray analysis of the gene expression profile in the midgut of silkworm infected with cytoplasmic polyhedrosis virus

In order to obtain an overall view on silkworm response to Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) infection, a microarray system comprising 22,987 oligonucluotide 70-mer probes was employed to compare differentially expressed genes in the midguts of BmCPV-infected and normal silkworm larvae. At 72 h post-inoculation, 258 genes exhibited at least 2.0-fold differences in expression level. Out of these, 135 genes were up-regulated, while 123 genes were down-regulated. According to gene ontology (GO), 140 genes were classified into GO categories. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicates that 35 genes were involved in 10 significant (P < 0.05) KEGG pathways. The expressions of genes related to valine, leucine, and isoleucine degradation, retinol metabolism, and vitamin B6 metabolism were all down-regulated. The expressions of genes involved in ribosome and proteasome pathway were all up-regulated. Quantitative real-time polymerase chain reaction was performed to validate the expression patterns of 13 selected genes of interest. The results suggest that BmCPV infection resulted in the disturbance of protein and amino acid metabolism and a series of major physiological and pathological changes in silkworm. Our results provide new insights into the molecular mechanism of BmCPV infection and host cell response.     

基于转录组的肝豆状核变性调控网络的构建和分析

本研究旨在利用生物信息学方法构建经铜诱导的ATP7B基因敲除HepG2细胞系的转录调控网络。探讨关键转录因子在肝豆状核变性发生、发展中的潜在作用机制。收集公共基因表达数据库(gene expression omnibus, GEO)中包含野生型、ATP7B基因敲除型、铜诱导的野生型和铜诱导的ATP7B基因敲除型HepG2细胞系数据。筛选由铜诱导产生的差异表达基因(differentially expressed genes,DEGs)后进行基因本体论(gene ontology,GO)、京都基因和基因组百科全书(Kyoto encyclopedia of genes and genomes, KEGG)富集分析。基于蛋白相互作用网络,识别疾病关键基因和功能模块,并对关键功能模块中的基因进行富集分析。最后,构建转录调控网络,筛选核心转录因子。共筛选出1 034个差异表达基因,其中上调525个,下调509个。上、下调关键功能模块分别包括了3785个和3931个基因。关键功能模块中的基因主要定位于细胞-基质连接、染色体、剪接复合体、核糖体等区域,共同参与了mRNA加工、组蛋白修饰、RNA剪切...     

Genome-wide gene expression analysis suggests an important role of hypoxia in the pathogenesis of endemic osteochondropathy Kashin-Beck disease

Kashin-Beck Disease (KBD) is an endemic osteochondropathy, the pathogenesis of which remains unclear now. In this study, we compared gene expression profiles of articular cartilage derived respectively from KBD patients and normal controls. Total RNA were isolated, amplified, labeled and hybridized to Agilent human 1A 22 k whole genome microarray chip. qRT-PCR was conducted to validate our microarray data. We detected 57 up-regulated genes (ratios ≥2.0) and 24 down-regulated genes (ratios ≤0.5) in KBD cartilage. To further identify the key genes involved in the pathogenesis of KBD, Bayesian analysis of variance for microarrays (BAM) software was applied and identified 12 potential key genes with an average ratio 6.64, involved in apoptosis, metabolism, cytokine & growth factor and cytoskeleton & cell movement. Gene Set Enrichment Analysis (GSEA) software was used to identify differently expressed gene ontology categories and pathways. GSEA found that a set of apoptosis, hypoxia and mitochondrial function related gene ontology categories and pathways were significantly up-regulated in KBD compared to normal controls. Based on the results of this study, we suggest that chronic hypoxia-induced mitochondrial damage and apoptosis might play an important role in the pathogenesis of KBD. Our efforts may help to understand the pathogenesis of KBD as well as other osteoarthrosis with similar articular cartilage lesions.     

脂多糖调控对大鼠心脏微血管内皮细胞的转录组分析

目的探索脂多糖(LPS)对大鼠心脏微血管内皮细胞(rCMECs)转录组的调节作用。方法对照组(正常培养rCMECs),LPS组(100 ng/mL LPS处理6 h的rCMECs),每组进行3个生物学重复转录组测序。得到差异基因后,使用实时定量PCR对部分差异基因mRNA的表达进行验证。分别对上调和下调基因进行GO和KEGG富集,并对差异基因进行共表达网络分析。采用独立t检验进行统计学分析。结果LPS处理后,265个基因表达上调,118个基因的表达下调。前10个最显著上调基因为:Mt2a、Cyp7b1、Sod2、Icam1、Ccl2、AC128848.1、Mt1、Cebpd、Serpinb2和Tnfrsf11b。前10个最显著下调基因为:Cavin2、Ankrd1、Edn1、Prss35、Lmod1、Dhrs3、Ttc22、Sema6a、Map2k3和Sema7a。定量PCR的结果表明Mt2a、Sod、Ccl2、Cxcl1、Icam1和Vcaml基因的表达得到了上调(P均<0.01);而Cavin2、Ankrd1、Edn1和Prss35基因表达下调(P均<0.05)。GO和KEGG富集的结果表明,上调基因与内皮细胞对炎性免疫细胞的趋化作用和黏附作用密切相关;而下调基因则是与钙离子信号和G蛋白相关通路以及内皮通透性增加有关。此外,差异基因进行共表达网络分析发现Sod2处于核心位置,提示其可能与LPS诱导的rCMECs的各种变化密切相关。结论LPS调控了rCMECs中大量与炎性免疫细胞进入心肌组织相关基因的表达。     

The influence of adiponectin on the transcriptomic profile of porcine luteal cells

Reproductive functions are closely related to nutritional status. Recent studies suggest that adiponectin may be a hormonal link between them. Adiponectin is an adipocytokine, abundantly expressed in adipose tissues. It plays a dominant role in lipid and carbohydrate metabolism by stimulating fatty acid oxidation, decreasing plasma triglycerides, and increasing cells’ sensitivity to insulin and has direct antiatherosclerotic effects. The hormone is also postulated to play a modulatory role in the regulation of the reproductive system. The aim of this study was to identify differentially expressed genes (DE-genes) in response to adiponectin treatment of porcine luteal ovarian cells. The global expression of genes in the porcine ovary was investigated using the Porcine (V2) Two-color gene expression microarray, 4?×?44 (Agilent, USA). Analysis of the microarray data showed that 701 genes were differentially expressed and 389 genes showed a fold change greater than 1.2 (p?<?0.05). Among this number, 186 genes were up-regulated and 203 were down-regulated. The list of DE-genes was used for gene ontology analyses. The biological process list was generated from up-regulated and down-regulated DE-genes. We found that up-regulated products of DE-genes take part in 30 biological processes and down-regulated products in 9. Analysis of the interaction network among DE-genes showed that adiponectin interacts with genes involved in important processes in luteal cells. These results provide a basis for future work describing the detailed interactions and relationships explaining local regulation of adiponectin actions in the ovary of pigs.     

小儿脓毒症休克相关基因的筛选及网络构建

为探究脓毒症休克与SIRS的差异表达基因及网络的构建,筛选潜在的核心基因,从GEO数据库下载相关基因表达谱GSE26378,数据分为脓毒症休克与SIRS各29个样本,通过在线软件GCBI对其进行标准化及差异基因筛选;对差异基因进行GO分析;基于KEGG进行功能通路分析以及基因信号网络分析;差异基因共表达网络分析。结果表明:两组中总共有1 456个基因被识别为差异基因(P0.05),与SIRS组相比,脓毒症休克组中有条859条下调基因,597条上调基因。GO功能富集分析显示差异基因主要参与了细胞周期、细胞免疫、细胞代谢。KEGG功能通路分析显示差异基因主要参与了MAPK信号通路、P53信号通路、wnt信号通路、细胞凋亡信号通路,细胞周期受体信号通路等。共表达分析发现基因CCNB1、NUSAP1、OIP5、SHCBP1、ZWINT、TOP2A、DLGAP5等位于网络中央部位,而基因信号网络分析发现基因PLCB1、PIK3CA、STAT3、CAMK2D、PRKCB、CREB1位于网络核心。基因芯片分析有助于发现脓毒症休克与SIRS患儿外周血单核细胞在转录组学上的改变,而生物信息学网络分析有助于发现潜在的靶点。     

Monitoring expression profiles of Arabidopsis genes during cold acclimation and deacclimation using DNA microarrays

A comparative analysis of gene expression profiles during cold acclimation and deacclimation is necessary to elucidate the molecular mechanisms of cold stress responses in higher plants. We analyzed gene expression profiles in the process of cold acclimation and deacclimation (recovery from cold stress) using two microarray systems, the 7K RAFL cDNA microarray and the Agilent 22K oligonucleotide array. By both microarray analyses, we identified 292 genes up-regulated and 320 genes down-regulated during deacclimation, and 445 cold up-regulated genes and 341 cold down-regulated genes during cold acclimation. Many genes up-regulated during deacclimation were found to be down-regulated during cold acclimation, and vice versa. The genes up-regulated during deacclimation were classified into (1) regulatory proteins involved in further regulation of signal transduction and gene expression and (2) functional proteins involved in the recovery process from cold-stress-induced damages and plant growth. We also applied expression profiling studies to identify the key genes involved in the biosynthesis of carbohydrates and amino acids that are known to play important roles in cold acclimation. We compared genes that are regulated during deacclimation with those regulated during rehydration after dehydration to discuss the similarity and difference of each recovery process.Electronic Supplementary Material Supplementary materials are available for this article at     

Functional dissection of drought-responsive gene expression patterns in Cynodon dactylon L.

Water deficit is one of the main abiotic factors that affect plant productivity in subtropical regions. To identify genes induced during the water stress response in Bermudagrass (Cynodon dactylon), cDNA macroarrays were used. The macroarray analysis identified 189 drought-responsive candidate genes from C. dactylon, of which 120 were up-regulated and 69 were down-regulated. The candidate genes were classified into seven groups by cluster analysis of expression levels across two intensities and three durations of imposed stress. Annotation using BLASTX suggested that up-regulated genes may be involved in proline biosynthesis, signal transduction pathways, protein repair systems, and removal of toxins, while down-regulated genes were mostly related to basic plant metabolism such as photosynthesis and glycolysis. The functional classification of gene ontology (GO) was consistent with the BLASTX results, also suggesting some crosstalk between abiotic and biotic stress. Comparative analysis of cis-regulatory elements from the candidate genes implicated specific elements in drought response in Bermudagrass. Although only a subset of genes was studied, Bermudagrass shared many drought-responsive genes and cis-regulatory elements with other botanical models, supporting a strategy of cross-taxon application of drought-responsive genes, regulatory cues, and physiological-genetic information.     

转录组分析揭示东方蜜蜂微孢子虫侵染意大利蜜蜂的分子机制

【目的】本研究旨在通过差异表达基因(differentially expressed gene, DEG)分析以及毒力因子和其他侵染相关因子分析,在转录组水平揭示东方蜜蜂微孢子虫Nosema ceranae侵染意大利蜜蜂Apis mellifera ligustica的分子机制。【方法】基于前期已获得高质量的东方蜜蜂微孢子虫纯化孢子(NcCK)及侵染意大利蜜蜂工蜂7和10 d的东方蜜蜂微孢子虫(分别为NcT1和NcT2)转录组数据,根据P≤0.05且|log_2(Fold change)|≥1的标准,通过比较分析筛选出NcCK vs NcT1, NcCK vs NcT2和NcT1 vs NcT2比较组的DEG。通过相关生物信息学软件对上述DEG进行Venn分析、GO分类和KEGG代谢通路富集分析。根据Nr和KEGG数据库注释信息和相关文献进行对东方蜜蜂微孢子虫的毒力因子和侵染相关因子的统计和分析。通过RT-qPCR验证转录组数据及DEG表达趋势。【结果】从NcCK vs NcT1, NcCK vs NcT2和NcT1 vs NcT2比较组分别鉴定出1 397, 1 497和52个DEG。Venn分析结果显示各比较组共有的上调和下调基因分别为10和1个。GO分类结果显示,NcCK vs NcT1和NcCK vs NcT2中DEG富集数最多的功能条目为代谢进程、细胞进程、单组织进程、细胞、细胞组件、细胞器、催化活性和结合,而NcT1 vs NcT2中DEG富集数最多的是代谢进程、细胞进程、单组织进程、催化活性和结合。KEGG代谢通路富集分析结果显示,NcCK vs NcT1和NcCK vs NcT2中DEG分别富集到80和79条通路;富集在糖酵解/糖异生和MAPK信号通路的上调基因数量多于下调基因。毒力因子分析结果显示,孢壁蛋白9基因和孢壁蛋白12基因在NcCK vs NcT1和NcCK vs NcT2中均下调表达,孢壁蛋白8基因仅在NcCK vs NcT1中表达量下调;此外孢壁蛋白前体基因、孢壁和锚定盘复合蛋白基因、几丁质合酶基因、极管蛋白基因、蓖麻毒素B凝集素基因的表达水平在NcCK vs NcT1和NcCK vs NcT2中表现为上调。侵染相关因子分析结果表明,糖酵解途径的3个关键酶基因在NcCK vs NcT1和NcCK vs NcT2中上调表达;3个涉及ATP/ADP移位酶的基因在NcCK vs NcT1和NcCK vs NcT2中上调表达,但有1个表达量下调;2个涉及ABC转运蛋白的基因在NcCK vs NcT1和NcCK vs NcT2中上调表达,另有4个下调表达。RT-qPCR结果证实了本研究中转录组数据及DEG表达趋势的真实可靠性。【结论】本研究通过比较分析解析东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程的转录组动态,揭示了孢壁蛋白、孢壁和锚定盘复合蛋白、几丁质酶、极管蛋白和蓖麻毒素B凝集素等毒力因子编码基因,以及己糖激酶、丙酮酸激酶、6-磷酸果糖激酶、ATP/ADP移位酶和ABC转运蛋白等侵染相关因子编码基因在病原增殖中扮演重要角色,为阐明东方蜜蜂微孢子虫的侵染机制提供了基础。