首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 765 毫秒
1.
应用Illumina HiSeq测序平台,对过表达FST基因的细胞和正常细胞进行转录组测序,并对筛选的差异表达基因进行基因功能注释及信号通路分析。结果表明:过表达FST后共有286个基因的转录发生了变化,其中176个上调,110个下调。这些差异基因共得到44个GO功能注释,其中单一生物过程、代谢过程、细胞过程;细胞器、细胞区域、细胞和黏合terms下聚集的差异基因数量最多,均>100个。共发现9条显著富集的Pathway,包括与细胞增殖相关的细胞周期、DNA复制、碱基切除修复、p53信号通路等,与卵母细胞发育相关的卵母细胞减数分裂和孕激素介导的卵母细胞成熟通路。说明在成纤维细胞中过表达FST后,最显著的特征是与细胞增殖分化相关的基因和生物学通路被富集。  相似文献   

2.
氮素是小麦生长发育过程中重要的营养元素,利用转录组技术鉴定参与小麦低氮胁迫分子调控网络的基因,对揭示小麦耐低氮分子机理、开展耐低氮育种具有重要参考意义。采用Illumina HiseqTM 2500高通量测序技术,对小麦品种晋麦47正常生长和低氮胁迫下的叶片进行转录组测序,筛选差异表达基因并在GO、KEGG 数据库中进行比对注释,分析小麦低氮胁迫响应相关基因。结果显示,对照组与低氮组分别获得高质量序列52 383 726和52 192 061条,检测到差异表达基因1 267个,其中上调表达基因179个,下调表达基因1 088个。差异基因GO功能注释到3个大类的44个功能组。差异表达基因被注释到178个途径上,主要富集于氨基酸代谢、碳水化合物代谢、脂类代谢和信号传导等途径。转录因子分析发现,在低氮条件下变化明显的转录因子家族包括WRKY、MYB和NAC等。  相似文献   

3.
【目的】对盐胁迫下海马齿根系进行转录组测序分析,挖掘海马齿根系耐盐相关基因,为揭示海马齿耐盐的分子机制提供参考。【方法】利用Illumina测序技术对0 mmol/L NaCl (对照组)和400 mmol/L NaCl胁迫处理(盐胁迫处理组)下的海马齿根系进行转录组测序分析,从中筛选出差异表达基因,选取13个基因进行实时荧光定量PCR (qRTPCR)检测,以验证转录组数据的可靠性。【结果】在海马齿根系转录组中共鉴定出305145个转录本,平均长度为622 bp,其中,对照组有146177个长度>300 bp的转录本,盐胁迫处理组有72173个长度>300 bp的转录本;共有65535条Unigenes在Nr、GO、Swiss-Prot、COG和KEGG五大数据库注释成功,占Unigenes总数的52.36%。对照组和盐胁迫处理组共有65535个差异Unigenes,其中,有182个热休克蛋白基因。对照组和盐胁迫处理组间共有24042个差异表达基因,从中选取13个基因进行qRT-PCR检测,结果显示,9个基因表达上调,其余4个基因表达下调,与转录组测序结果一致。24042个差异表达基因中,共有10106个显著差异基因富集到129条代谢通路,其中富集程度排名前10的代谢途径为核糖体、次级代谢生物合成、RNA转运、内吞作用、剪接体、甘油磷脂代谢、内质网加工、吞噬、醚脂类代谢和植物-病原体相互作用,参与盐胁迫相关的硫代谢、脯氨酸积累、活性氧(ROS)代谢、与盐胁迫相关的钙信号通路和过氧化氢代谢等途径的差异基因上调。【结论】在盐胁迫下海马齿差异表达基因如小分子量热激蛋白基因、抗氧化酶相关基因及与离子交换相关基因发挥了重要调控作用。  相似文献   

4.
氮素是小麦生长发育过程中重要的营养元素,利用转录组技术鉴定参与小麦低氮胁迫分子调控网络的基因,对揭示小麦耐低氮分子机理、开展耐低氮育种具有重要参考意义。采用Illumina HiseqTM 2500高通量测序技术,对小麦品种晋麦47正常生长和低氮胁迫下的叶片进行转录组测序,筛选差异表达基因并在GO、KEGG 数据库中进行比对注释,分析小麦低氮胁迫响应相关基因。结果显示,对照组与低氮组分别获得高质量序列52 383 726和52 192 061条,检测到差异表达基因1 267个,其中上调表达基因179个,下调表达基因1 088个。差异基因GO功能注释到3个大类的44个功能组。差异表达基因被注释到178个途径上,主要富集于氨基酸代谢、碳水化合物代谢、脂类代谢和信号传导等途径。转录因子分析发现,在低氮条件下变化明显的转录因子家族包括WRKY、MYB和NAC等。  相似文献   

5.
为探究clpV基因在禽致病性大肠杆菌(Avian pathogenic Escherichia coli, APEC)感染雏鸡气管黏膜过程中发挥的作用及机制。将APEC野生株(DE17)及其基因缺失株(DE17ΔclpV)感染7日龄雏鸡气管,12、24 h后收集雏鸡气管黏膜细胞进行转录组学测序,筛选野生株和缺失株的差异表达基因,进行GO和KEGG富集分析。结果显示,在感染12 h后,共筛选到108个差异表达基因(36个上调,72个下调);在感染24 h后,共筛选到59个差异表达基因(34个上调,25个下调)。GO分析表明,差异表达基因主要富集在细胞内信号转导、RNA聚合酶Ⅱ的转录正调控等生物学过程;蛋白质结合、ATP合成、DNA结合等分子功能类;生物膜及膜的组成成分、细胞质等细胞组分功能;KEGG分析表明,差异表达基因主要富集在紧密连接通路、内质网中的蛋白质加工通路、细胞因子与细胞因子受体相互作用通路、内吞通路、Wnt信号通路等。结论为clpV基因缺失后,会使部分差异表达基因的表达量下调,对细胞因子-细胞因子受体相互作用等通路产生影响。  相似文献   

6.
  目的  夏季高温可导致牡丹Paeonia suffruticosa生长受限,花期缩短,观赏品质降低。研究牡丹耐热基因对牡丹热胁迫的分子机制。  方法  以牡丹品种‘羽红’‘Yuhong’为材料,对高温(40 ℃)和对照(25 ℃)处理后的叶片进行转录组测序,分析响应高温表达的差异基因;同时利用实时荧光定量PCR(qRT-PCR)验证差异基因,并对热激蛋白基因HSPs进行时空表达分析。  结果  测序共获得45.97 Gb数据,与对照组(25 ℃)相比,高温处理后4 220个基因上调表达,3 453个基因下调表达。通过基因本体(GO)分析发现,这些差异基因主要富集于代谢等生物学过程、细胞和膜结构等细胞成分以及结合和催化活性等分子功能条目;另外,通过全基因组及代谢途径数据库(KEGG)分析发现,碳代谢途径的差异基因数量最多。对牡丹热激蛋白PsHSP基因定量分析发现,PsHSP基因随着高温处理时间呈上升表达趋势,24 h达到最高,之后便呈现下降趋势。  结论  高温显著影响牡丹的代谢和生物合成等过程,进而影响牡丹的生长和发育。PsHSP基因可在短期内迅速响应高温胁迫并参与牡丹的耐热调控。图 4表6参34  相似文献   

7.
为了解青田田鱼(Cyprinus carpio var. qingtianensis, Paddy Field Carp)对稻田急性溶氧变化的水环境适应性, 依据青田田鱼对低氧耐受性参数及当地稻田水环境溶氧日变化的实际情况模拟急性低氧和复氧环境,基于RNA-seq技术对青田田鱼幼鱼常氧对照组[温度:(25.37±0.76)℃,溶氧:(7.07±0.39)mg/L,CB组],低氧胁迫6 h组[温度:(25.65±0.76)℃,溶氧:(0.57±0.16)mg/L,HB组]和复氧恢复6 h组[温度:(25.78±0.58)℃,溶氧:(7.21±0.53)mg/L,RB组]脑组织进行转录组分析。结果表明:HB组与CB组文库比较发现存在1 094个差异基因,其中,730个上调,364个下调;RB组与HB组文库比较发现2 288个差异基因,其中,803个上调,1 485个下调;RB组与CB组文库比较发现201个差异基因,其中,80个上调,121个下调。进一步将差异表达基因进行 GO 功能注释和 KEGG 富集分析发现差异基因主要富集在糖酵解/糖异生、HIF-1信号、mTOR信号及VEGF信号等通路上。通过RT-qPCR对糖酵解/糖异生及相关通路中的6个差异基因进行验证,证实了转录组测序结果的可靠性。这些结果为研究青田田鱼在急性溶氧变化环境下的生理调节机制提供了分子生物学的基础研究。  相似文献   

8.
夏季奶牛热应激问题给奶业生产造成了巨大的经济损失。为揭示奶牛发生热应激后乳腺组织应答反应的分子机制,通过研究热应激奶牛乳腺组织差异表达基因,探究差异基因对乳腺组织转录调控的影响。以8头中国荷斯坦牛为研究对象,分别采集热应激期(8月,温湿指数THI=83.8)和非热应激期(3月,THI=65.8)各4头奶牛乳腺组织,利用HiSeq2000进行转录组测序,并进行生物信息学分析,共发现96个差异表达基因(差异倍数>2),其中上调表达46个,下调表达50个。GO分类结果表明,注释到细胞组分、分子功能和生物学过程的差异基因数量分别为41、38和36个(P<0.05)。KEGG通路分析结果表明,差异基因共富集到18条信号通路(P<0.05),其中6条与疾病有关,9条与代谢有关,此外,与免疫应答相关的细胞因子-细胞因子受体相互作用和NOD样受体信号通路明显富集。通过比较热应激和非热应激奶牛乳腺组织转录组,分析差异基因相关信号通路,为进一步了解奶牛发生热应激反应的分子机制奠定基础。  相似文献   

9.
为深入分析花生抗旱分子机制提供理论依据,采用转录组测序和小RNA测序技术,对桂花37花生在15% PEG6000模拟干旱胁迫下基因和miRNA的表达模式进行分析。结果表明:转录组测序发现有 13个显著持续上调的差异基因,其中7个为HSP分子伴侣家族基因,2个RING E3泛素连接酶,1个为谷胱甘肽-S-转移酶,1个柱头特异性STIG1样蛋白1,1个免疫球蛋白结合蛋白和1个 DnaJ同源B亚家族成员13;同时筛选出13个显著下调的差异基因;在小RNA测序数据中筛选到8个持续下调表达的miRNA,分别为ahy-mir156-3p、 ahy-miR159、ahy-miR167-5p、ahy-miR156a、ahy-miR156b-5p、novel_mir103、novel_mir83和novel_mir70, 持续下调表达的miRNA通过负调控方式参与花生干旱胁迫过程。  相似文献   

10.
【 目 的】 探 究 铁 皮 石 斛(Dendrobium catenatum Lindl.)DcCIPK24 的 下 游 响 应 基 因, 为 研 究DcCIPK24 提高耐盐性的分子调控途径提供指导。【方法】利用 Illumina NovaSeq 6000 测序平台对 200 mmol/L NaCl 处理和对照组的转 DcCIPK24 拟南芥、野生型拟南芥进行测序,分析差异表达基因的富集通路,并对通路中的差异基因进行荧光定量 PCR 验证。【结果】盐胁迫下,从 2 个转基因拟南芥株系与野生型拟南芥的比较组合中共鉴定出 96 个差异表达基因;GO 注释分析发现 96 个差异基因主要被富集在分子功能相关通路;KEGG 富集分析发现差异基因主要被注释在氨基糖和核苷酸糖代谢、植物 MAPK 信号通路和甘油酯代谢通路中,选取上调表达的差异基因 AtGPAT5、AtSIRK、AtMEE25、AtUGE5 和下调表达基因 AtAMT1-2、AtATTI3、AtCYP71A25、AtLHCB4.3 进行实时荧光定量 PCR 验证,结果表明盐胁迫下 8 个差异基因表达趋势与转录组数据基本一致。【结论】盐胁迫下,DcCIPK24 主要通过氨基糖和核苷酸糖代谢途径、植物 MAPK 信号途径和甘油酯代谢通路响应耐盐。  相似文献   

11.
【目的】对高温胁迫下萝卜幼苗进行转录组测序分析,筛选不同耐热性萝卜品种间的差异表达基因,为阐明萝卜幼苗响应高温胁迫的分子机制提供理论参考。【方法】以耐热萝卜品种1116T和热敏感萝卜品种Wr129为对象,对其进行高温(40℃)胁迫处理0(常温,对照)和24 h,提取各时间点样品总RNA进行转录组测序,以|log2 FoldChange|≥1,P<0.05且TPM≥10为标准,筛选出差异表达基因(DEGs),并进行GO功能注释和KEGG代谢通路富集分析。【结果】随着高温胁迫时间的增加,Wr129和1116T幼苗茎叶逐渐发生萎蔫,胁迫24 h后Wr129叶片萎蔫下垂且干枯变黄,而1116T叶片轻度萎蔫下垂,仍保持绿色。12个文库共获得81.45 Gb Clean data,平均有62.30%的reads比对到萝卜参考基因组的一个位置。从2个不同耐热性品种中共鉴定到2701个差异表达基因。GO功能注释结果显示,2629个差异表达基因在GO数据库中得到注释,其中,有559个差异表达基因被注释为细胞组分类别中的叶绿体被膜、叶绿体基质和叶绿体类囊体3个条目,470个差异表达基因被注释为与胁迫响应相关的生物学过程。KEGG代谢通路富集结果显示,1166个差异表达基因富集到124条代谢通路,其中,有74个差异表达基因富集到光合作用和光合作用-天线蛋白代谢通路,23个基因富集到卟啉和叶绿素代谢通路,大部分基因在高温胁迫下被诱导; 35个差异表达基因富集到谷胱甘肽代谢途径,其中14个谷胱甘肽-S-转移酶基因(GSTs)在高温胁迫下上调表达。高温胁迫24 h,Wr129和1116T幼苗叶片PSII的最大光化学效率(Fv/Fm)和总叶绿素含量均下降,但1116T的下降幅度较小。【结论】 Wr129和1116T的耐热性存在明显差异,其原因是1116T植株可通过调控光合作用和抗氧化系统相关基因的表达,减缓Fv/Fm和叶绿素含量下降速率,从而提高植株对高温胁迫的耐受能力。  相似文献   

12.
【目的】探究脊尾白虾(Exopalaemon carinicauda)鳃组织在渐变式低氧—复氧胁迫下的分子调控机制,为今后开展脊尾白虾耐低氧品系(种)的选育提供理论指导。【方法】通过模拟自然低氧环境的形成过程,分别于低氧处理0(对照)、3和6 h及复氧后1和8 h采集脊尾白虾鳃组织,利用Illumina HiSeqTM 4000测序平台进行转录组测序分析,经过滤和Trinity组装获得Unigenes,选取Nr、Swiss-Prot、KEGG和COG/KOG等数据库进行注释分析,在Omicsmart平台上完成差异表达基因筛选及其表达趋势分析,然后进行GO功能注释分析和KEGG信号通路富集分析,并随机选取5个差异表达基因进行实时荧光定量PCR验证。【结果】脊尾白虾鳃组织样本转录组测序数据经过滤后的Cleanreads进行Trinity组装共获得93227条Unigenes,其长度范围在201~35402 bp,平均长度为834 bp,N50长度为1352 bp。通过组间两两比较分析,共鉴定出4750个差异表达基因,其中上调差异表达基因3557个、下调差异表达基因2829个;超过50%的差异表达基因被显著富集到6种基因表达趋势模式中(P<0.01),具体表现为:Profile 0模式富集到415个基因,Profile 5模式富集到201个基因,Profile 11模式富集到371个基因,Profile 13模式富集到841个基因,Profile 17模式富集到387个基因,Profile 18模式富集到411个基因。6种基因表达趋势模式中的差异表达基因被注释到代谢进程、细胞进程、单一有机体进程、细胞、细胞零件、大分子复合物及催化活性等GO功能条目上;而KEGG信号通路富集分析结果显示,以Profile 13模式中的差异表达基因富集到最多信号通路(86条),其中呈显著富集的有8条,分别为核糖体、碳代谢、氧化磷酸化、氨基酸生物合成、内质网蛋白质加工、糖酵解/糖异生、谷胱甘肽代谢和蛋白输出。【结论】脊尾白虾鳃组织在受低氧胁迫早期通过合成蛋白质及提高代谢能力来抵御低氧环境,随着低氧胁迫时间的延长,物质合成和能量代谢活动均显著下降;但在复氧后随着复氧时间的延长,其蛋白质合成和能量代谢水平又逐渐升高恢复至常氧水平。  相似文献   

13.
【目的】基于转录组学对西方蜜蜂工蜂不同虫态间的差异表达基因(DEGs)进行筛选和功能注释分析,揭示与工蜂生长发育相关的信号通路,为深入解析工蜂生长发育的分子调控机理提供基础数据。【方法】以西方蜜蜂工蜂的3日龄幼虫、1日龄蛹和1日龄羽化工蜂3个虫态为研究对象,利用llumina NovaSeq 6000平台进行转录组测序,采用DESeq2筛选不同虫态样品间的表达差异基因,然后分别进行GO功能注释分析及KEGG信号通路富集分析,并通过实时荧光定量PCR进行验证。【结果】经转录组测序,在西方蜜蜂工蜂3日龄幼虫与1日龄蛹间筛选出4823个差异表达基因(51.86%上调,48.14%下调),在1日龄蛹与1日龄羽化工蜂间筛选出3295个差异表达基因(57.51%上调,42.49%下调),在3日龄幼虫与1日龄羽化工蜂间筛选出5267个差异表达基因(52.95%上调,47.05%下调)。GO功能注释分析结果显示,3日龄幼虫与1日龄蛹间的差异表达基因注释到43个GO功能条目,1日龄蛹与1日龄羽化工蜂间的差异表达基因注释到45个GO功能条目,3日龄幼虫与1日龄羽化工蜂间的差异表达基因注释到44个GO功能条目,主要涉及细胞过程、细胞部分及结合等。KEGG信号通路富集分析发现,3日龄幼虫与1日龄蛹间有2905个差异表达基因富集到332条KEGG信号通路上,其中17条KEGG信号通路呈显著富集,涉及核糖体、氧化磷酸化和昆虫激素生物合成等;1日龄蛹与1日龄羽化工蜂间有1644个差异表达基因富集到331条KEGG信号通路上,其中45条KEGG信号通路呈显著富集,涉及氧化磷酸化、生热作用和胰岛素分泌等;3日龄幼虫与1日龄羽化工蜂间有2958个差异表达基因富集到337条KEGG信号通路上,其中14条KEGG信号通路呈显著富集,涉及核糖体、蛋白酶体和胰岛素分泌等。6个随机挑选差异表达基因的实时荧光定量PCR检测结果与转录组测序结果相符,即转录组测序结果可靠。【结论】昆虫激素生物合成通路相关差异表达基因调控与西方蜜蜂工蜂各虫态JH滴度变化规律一致,氧化磷酸化信号通路则与各虫态的营养摄入和活动行为相关,而胰岛素分泌通路涉及各虫态的营养调控、脂肪体合成及细胞凋亡。可见,昆虫激素生物合成、胰岛素分泌和氧化磷酸化3种信号通路在西方蜜蜂工蜂幼虫、蛹和成虫的发育调控中发挥着重要作用。  相似文献   

14.
为探明骆驼斯氏副柔线虫雌虫和雄虫的转录组差异,挖掘出与性别决定相关的功能基因,采用Illumina HiSeqTM2000分别对其雌虫和雄虫样本进行转录组测序,构建cDNA文库,通过建库质量评估和组装效率评估后对其进行功能注释及相关生物信息学分析。结果显示,比较雌虫和雄虫的转录本,获得6 430个差异表达基因,其中2 131个基因上调;4 299个基因下调。对差异表达基因进行GO功能分类,有341 328和714个Unigenes分别注释到生物学过程、细胞组成和分子功能三大类41个分支中。KEGG数据库分析,有1 116个差异表达基因参与到198条KEGG通路中,显著富集在机体免疫调节,卵母细胞成熟等通路。功能聚类分析中,雌虫样本在生殖发育、产卵、生殖行为中高表达;雄虫在主要精子蛋白、丝氨酸/苏氨酸特异性蛋白磷酸酶中高表达。此外,在雄虫和雌虫中筛选出Tab-3,Tab-5、Fem-1和Mog-3等10种线虫性别决定相关基因,以及8种雄性性别特异表达基因和4种雌性性别特异表达基因,这些基因在虫体性别决定和性别特异性表达中发挥重要作用。本研究利用生物信息学技术预测出雌虫和雄虫差异表达基因在功能分类和代谢通路等方面的生物学特征,挖掘出性别决定基因和性别特异性表达相关基因,为进一步开展斯氏副柔线虫的功能基因组学研究、药物靶点预测和疫苗候选基因筛选提供有价值的数据资源。  相似文献   

15.
16.
【目的】探明贵紫麦1号小麦灌浆期变紫后和变紫前2个时期籽粒的转录组差异,发掘影响贵紫麦1号花青素合成的关键基因和关键酶,丰富小麦籽粒色素转录组数据信息,为转录因子的克隆及表达提供参考。【方法】利用Illumina Hiseq 2000TM高通量测序技术对贵紫麦1号籽粒变紫前和变紫后2个时期进行转录组测序、文库构建及建库质量评估,对测序结果进行信息学分析。采用TTM对read count数据进行标准化处理,随后用DEGseq进行差异分析,设定q-value<0.005且|log2 (fold change)|>1为阈值。通过筛选分析,获得两者间差异表达基因,按照无参转录组分析方法,对差异表达基因进行BLAST搜索,Nr数据库比对,GO功能富集及KEGG pathway分析,找出与花青素相关的关键基因和关键酶,并结合qRT-PCR验证所找到的关键基因及关键酶在不同时期的表达水平,掌握这些关键基因的信息。【结果】测序结果表明,贵紫麦1号变紫后和变紫前分别获得13.36 G和12.69 G的clean bases,clean reads为106 906 108条和101 547 534条,占原始序列的93.73%和94.90%。通过Trinity软件对所得clean reads进行拼接,共获得170 396条转录本,长度为119 020 625。拼接clean reads后获得119 572条Unigenes。在BLAST搜索中,119 572个高质量独特序列中有86 004条(71.92%)Unigenes与现有基因模型具有至少1个显著匹配。在Nr数据库比对结果鉴定了至少5种具有与来自节节麦、乌拉尔图小麦、二穗短柄草、大麦、小麦等已知基因同一性且序列相似性高的Unigenes。KOG数据库比对结果显示,注释成功的基因按KOG的26个group进行分类,注释在一般功能基因,蛋白质翻译后修饰与转运、分子伴侣及翻译、核糖体结构与生物合成等类别基因所占比重较大,分别为15.79%、14.51%和10.54%。643个差异基因中,236个呈上调趋势,407个呈下调趋势。GO注释表明,按照基因参与的生物过程、所处的细胞组分、具有的分子功能下一层级分类,共44个分类,差异基因显著富集在碳水化合物代谢过程(GO:0005975,16.03%)、应激反应(GO:0006950,10.83%)和水解酶活性分子功能(GO:0016787,34.84%)等类别中。KEGG pathway富集分析可知,353个差异基因富集到153条相关通路上,其中淀粉与蔗糖代谢、苯丙素生物合成、类黄酮生物合成等通路富集显著。类黄酮生物合成途径相关基因共66个,2条相关上调表达Unigenes,涉及查尔酮酶、隐色花色素双加氧酶2个关键酶基因,log2(fold change)分别为3.4164和2.1258。对所得关键基因进行qRT-PCR验证,证实查尔酮酶、隐色花色素双加氧酶在贵紫麦中1号中表达量呈明显上调趋势,与转录组测序分析结果一致,测序结果可靠度高。【结论】比较分析贵紫麦1号籽粒变紫后和变紫前2个时期转录组测序结果,获得大量Unigenes数据及差异表达基因相关信息,明确类黄酮代谢途径中2个关键酶基因(CHS和ANS)在调控贵紫麦1号籽粒花青素合成过程中作用显著。  相似文献   

17.
意大利蜜蜂响应东方蜜蜂微孢子虫胁迫的免疫应答   总被引:2,自引:2,他引:0  
【目的】通过对意大利蜜蜂(Apis mellifera ligustica,简称意蜂)工蜂中肠响应东方蜜蜂微孢子虫(Nosema ceranae)胁迫的差异表达基因(differentially expressed gene,DEG)及免疫通路进行深入细致的分析,揭示宿主的细胞和体液免疫应答,为关键免疫应答基因的筛选及功能研究打下基础。【方法】基于前期获得的正常及东方蜜蜂微孢子虫胁迫的意蜂7日龄和10日龄工蜂中肠(Am7CK、Am7T、Am10CK、Am10T)转录组数据,按照FDR≤1,P≤0.05和|log2 fold change|≥1的标准筛选出各组的DEG,利用相关生物信息学软件对DEG进行Pearson相关性分析、Venn分析、GO分类和KEGG代谢通路富集分析,进而对免疫通路富集基因进行统计和分析,利用实时荧光定量PCR(real-time quantitative PCR,RT-qPCR)验证转录组数据的可靠性。【结果】 差异分析结果显示,Am7CK vs Am7T比较组包含472个上调基因和385个下调基因;Am10CK vs Am10T比较组包含611个上调基因和360个下调基因。Venn分析结果显示,两个比较组特有的DEG分别为739和853个,共有的DEG为118个。GO分类结果显示,Am7CK vs Am7T比较组中上调和下调基因分别涉及23和29个功能条目,其中富集上调基因数最多的前5位分别为结合、催化活性、代谢进程、细胞进程和单组织进程;富集下调基因数最多的前5位分别为代谢进程、单组织进程、催化活性、细胞进程和结合。Am10CK vs Am10T比较组中上调和下调基因分别涉及36和26个功能条目,其中富集上调基因数最多的前5位分别为单组织进程、结合、细胞进程、催化活性和代谢活性;富集下调基因数最多的前5位分别为结合、细胞进程、催化活性、代谢进程和单组织进程。KEGG代谢通路富集分析结果显示,Am7CK vs Am7T比较组中上调和下调基因分别富集在38和33条代谢通路,富集上调基因数最多的前5条分别为胆汁分泌、内质网蛋白加工、泛素介导的蛋白水解、PI3K-Akt信号通路和神经营养因子信号通路;富集下调基因数最多的前5条分别为胞质DNA传感通路、嘌呤代谢、嘧啶代谢、RNA聚合酶和核糖体;涉及泛素介导的蛋白水解等3条细胞免疫通路,以及PI3K-Akt信号通路等7条体液免疫通路。Am10CK vs Am10T比较组中上调和下调基因分别富集在54和43条代谢通路,富集上调基因数最多的前5条分别为Hippo信号通路、药物代谢-细胞色素P450、P450对外源物质的代谢、泛素介导的蛋白水解和鞘脂类代谢;富集下调基因数最多的前5条分别为mRNA监测、鞘脂类信号通路、果糖和甘露糖代谢、半乳糖代谢和鞘脂类代谢;涉及泛素介导的蛋白水解等7条细胞免疫通路,以及NF-κB信号通路等2条体液免疫通路。RT-qPCR验证结果显示6个随机挑选的DEG的表达量变化趋势与测序结果一致,证实了本研究中测序数据的可靠性。进一步分析发现意蜂7日龄和10日龄工蜂中肠的NF-κB信号通路均被东方蜜蜂微孢子虫激活,随即启动了3种抗菌肽apidaecin、defensin-1和hymenoptaecin的合成,体现了它们在宿主抵御东方蜜蜂微孢子虫入侵中的重要性。【结论】在转录组水平解析了意蜂工蜂中肠对东方蜜蜂微孢子虫的免疫应答,揭示宿主在胁迫前期同时作出细胞和体液免疫应答,前者可能在抵御病原入侵方面发挥主要作用;宿主的细胞免疫在胁迫后期持续增强,但体液免疫较大程度地减弱;泛素介导的蛋白水解通路及富集DEG,NF-κB信号通路及富集DEG,以及抗菌肽编码基因apidaecindefensin-1hymenoptaecin可能在意蜂工蜂对东方蜜蜂微孢子虫的免疫应答中起到关键作用。  相似文献   

18.
为探究重金属污染地生态修复先锋植物——盐肤木对重金属胁迫的分子响应机制,设置3个铅胁迫水平(0、250、1000mg·kg~(-1)),应用Illumina HiSeq~(TM)2000进行盐肤木根系的转录组测序,并通过Gene Ontology(GO)与Pathway功能显著富集分析方法,分析说明盐肤木响应不同铅胁迫水平下的差异表达基因。结果表明,盐肤木在遭受轻度铅胁迫时,主要有4类水解酶与2类磷酸酶相关基因显著富集,说明盐肤木主要通过调节细胞内水解酶与磷酸酶相关基因来抵御胁迫,而在重度铅胁迫下,主要有6类氧化还原酶相关基因出现显著富集,此时氧化还原酶相关基因起到主要调节作用;盐肤木在铅胁迫下细胞受损,在轻、重度铅胁迫下,分别有24、16类细胞相关基因显著富集,但其自身可通过调节细胞内细胞器、细胞质等相关的细胞运动以应对逆境;核糖体代谢通路是盐肤木适应铅胁迫的主要代谢通路。研究表明,核糖体相关基因是盐肤木响应铅胁迫的主要应答与调节基因。  相似文献   

19.
  目的  木棉是我国南方地区重要的观赏植物,研究不同时期花蕾和花朵的花瓣基因表达差异,揭示花色变异的遗传调控机制,为建立花色定向育种技术提供科学依据。  方法  以木棉不同发育时期深红色花和黄色花的花瓣为研究对象,利用Illumina HiSeqTM 4000开展转录组测序;分别采用DESeq2和Kyoto Encyclopedia of Genes and Genomes(KEGG)数据库进行差异表达基因鉴定和通路富集分析。  结果  测序共获得75 190个单基因,4 772个单基因能被公共数据库注释;基因功能富集分析显示,基因主要富集在基本功能预测、信号转导机制和转录后修饰、蛋白折叠、分子伴侣等通路。共获得不同时期差异表达基因10 397个,显著富集在29个生物学通路,主要包括光合作用、代谢和植物激素信号转导通路等。其中,参与苯丙烷生物合成通路的差异表达基因有72个,参与类胡萝卜素、黄酮类生物合成通路和苯丙氨酸代谢通路的差异表达基因分别为25个,这些通路和基因均与花青素的生物合成相关;另有4个差异表达基因显著富集在甜菜碱的生物合成通路中。qRT-PCR数据验证了转录组数据的可靠性。  结论  (1)光合作用、代谢、植物激素信号转导、黄酮类生物合成、苯丙烷生物合成、苯丙氨酸代谢及能量代谢等通路相关基因可能参与花瓣的发育过程。(2)黄酮类生物合成通路相关基因在深红色花发育的花蕾中期与花朵期显著高表达,可能是花色呈现红色的主要原因。(3)类胡萝卜素生物合成通路的关键合成酶基因的部分家族成员在黄色花发育的花蕾期和花蕾中期显著高表达,可能是导致花色呈现黄色的主要原因。   相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号