首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 828 毫秒
1.
【目的】对青花菜花蕾进行转录组测序分析,并挖掘与蜡粉合成相关基因,为探明青花菜花球表面蜡粉形成的分子机制提供理论参考。【方法】分别提取野生型和蜡粉缺失型青花菜花球总RNA,采用Illumina HiSeqTM 2500平台进行转录组测序,获得高质量Clean reads,采用Trinity进行序列组装后获得青花菜Unigene库,将获得的Unigene序列与Nr、Nt、KEGG、Pfam、KOG/COG、Swiss-Prot和GO数据库比对,获得基因功能注释信息;使用DESeq2进行差异表达分析。【结果】共获得44.68 Gb Clean data,De novo组装得到41244条Unigenes,N50长度为1847 bp。从所获得的Unigenes中筛选出8685个差异表达基因(DEGs)(上调基因5747个,下调基因2938个),共有8038个基因被注释到不同数据库,其中,5220个基因注释到Pfam数据库; 2066个基因注释到COG数据库,3866个基因注释到KOG数据库; 2580个差异表达基因被注释到75个转录因子家族中,注释最多的是MYB家族(235个);GO数据库中6095个差异表达基因注释到细胞组分、分子功能和生物学过程三大类的52个功能分类; KEGG数据库中,1671个差异表达基因富集到138条代谢通路,其中13个差异表达基因与脂肪酸合成有关,7个差异表达基因与蜡粉生物合成途径有关。【结论】转录因子MYB家族在调控青花菜蜡粉合成中发挥重要作用。蜡粉合成过程中相关酶基因的差异表达是调控青花菜蜡粉合成的关键,尤其是野生型和蜡粉缺失突变体中特异性表达的差异表达基因,可作为后续研究青花菜花球表面蜡粉形成分子机制的对象。  相似文献   

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 HiseqTM 2500高通量测序技术,对小麦品种晋麦47正常生长和低氮胁迫下的叶片进行转录组测序,筛选差异表达基因并在GO、KEGG 数据库中进行比对注释,分析小麦低氮胁迫响应相关基因。结果显示,对照组与低氮组分别获得高质量序列52 383 726和52 192 061条,检测到差异表达基因1 267个,其中上调表达基因179个,下调表达基因1 088个。差异基因GO功能注释到3个大类的44个功能组。差异表达基因被注释到178个途径上,主要富集于氨基酸代谢、碳水化合物代谢、脂类代谢和信号传导等途径。转录因子分析发现,在低氮条件下变化明显的转录因子家族包括WRKY、MYB和NAC等。  相似文献   

4.
应用SSH 技术构建了CTV诱导的枳的正向和反向cDNA文库,经菌液PCR检测,从两个文库中分别随机 挑选了100个阳性克隆测序,共获得169条表达序列标签(ESTs).用每条EST序列对应的克里曼丁基因组预测转 录序列进行GO 功能注释,去除在两个文库中均有出现的ESTs后,共得到与已知基因具有较高相似度,涉及39个 特异性表达基因的EST83条,并对4个特异性表达基因进行了实时定量PCR验证.结果表明:CTV侵染之后,枳 启动了相关抗逆防御机制来减少伤害,木质素代谢、RNA干扰、赤霉素、生长素以及光敏色素调控可能在抗CTV 中发挥了一定作用.  相似文献   

5.
6.
采用Illumina高通量测序平台对黑壳与金壳葡萄牙牡蛎(Crassostrea angulata)外套膜组织开展转录组测序研究,其中黑壳组总计产出5 857 718 092 nt数据,筛选后的干净阅读子58 004 560个;金壳组产出6 727 447 326 nt数据,干净阅读子66 614 854个。将发掘的新基因与COG、GO、KEGG、Swiss-Prot、NR数据库进行序列比对,最终得到各数据库注释的新基因数量分别为147、397、203、567、1 347个。通过对差异表达基因的序列进行生物学功能注释和聚类分析以及利用已知贝壳基质蛋白序列进行搜索,鉴定出14条可能与壳色表达有关的贝壳基质蛋白基因,其中5条下调基因可能与金壳基因表达相关。  相似文献   

7.
【目的】掌握不同生长规格墨瑞鳕个体间差异表达基因的表达特点,为其功能相关基因深度挖掘及分子遗传育种提供科学依据。【方法】挑选同一养殖条件下极大个体和极小个体的墨瑞鳕,构建肌肉组织cDNA文库后,采用Illumina HiSeqTM 4000测序平台对存在生长差异的墨瑞鳕肌肉组织进行转录组测序分析,获得的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等数据库中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行GO功能注释及KEGG信号通路富集分析,并采用MISA进行SSR鉴定分析。【结果】从墨瑞鳕肌肉组织中共测序获得39749条Unigenes,其长度范围在301~55230 bp,平均长度为1705 bp。注释到Nt、Nr、Swiss-Prot、Pfam数据库的Unigenes分别有27046、20824、18268和17772条,在7个数据库中均得到注释的Unigenes共计6742条,占Unigenes总数的16.96%。根据差异表达基因筛选条件P<0.05且|log2 Fold Change|>1,共筛选出722个差异表达基因,其中上调基因308个、下调基因414个。差异表达基因GO功能注释分析结果表明,注释基因数目较多的GO功能条目包括细胞过程、代谢过程、膜、细胞器及结合等;KEGG信号通路富集分析发现,差异表达基因被成功富集到234条信号通路上,主要涉及磷脂酰肌醇3激酶/蛋白激酶信号通路、MAPK信号通路、胰岛素信号通路及FoxO信号通路等。在39749条Unigenes中鉴定筛选出22120个SSRs,占Unigene总数的55.65%,SSR的平均间距为3063 bp。【结论】基于转录组测序分析获得的墨瑞鳕肌肉组织差异表达基因以发挥结合、细胞过程及代谢过程等功能为主,且主要富集在PI3K-Akt信号通路、核糖体信号通路、FoxO信号通路及细胞凋亡等能量代谢相关通路上,通过共同协调而对墨瑞鳕的生长发育起调控作用。  相似文献   

8.
旨在比较分析捻转血矛线虫敏感虫株和耐药虫株miRNA表达谱,探讨miRNA与捻转血矛线虫丙硫咪唑耐药性相关基因之间的调控机制。运用Illumina Hiseq2000平台进行测序并进行cDNA文库的构建,利用RNAhybrid、Miranda和TargetScan对测序得到的数据进行靶基因预测,并取交集作为miRNA的靶基因预测结果,使用Bowtie、RepeatMasker、MIREAP、Rfam、miRBase和DAVID等生物信息学分析软件和权威数据库系统,筛选出差异的miRNA以及与捻转血矛线虫耐药相关的miRNA并进行靶基因预测,同时对KEGG pathway和GO功能富集到的靶基因和可能参与调控的通路进行预测分析。结果表明,敏感虫株和耐药虫株共筛选显著差异表达的miRNA共294个,其中113个上调,181个下调。对miRNA和其靶向调节的mRNA进行关联分析,共关联到1 770个miRNA-mRNA对为负调控关系,其中涉及到274个差异的miRNA和603个差异的mRNA。显著差异表达的miRNA靶基因被注释到了1 430条GO terms和327条KEGG pathway,富集到FoxO信号通路(FoxO signaling pathway),mTOR信号通路(mTOR signaling pathway),PI3K-Akt信号通路(PI3K-Akt signaling pathway)等与耐药相关的一些抗性通路。综上,通过对miRNA表达谱分析以及miRNA-mRNA靶向分析,发现捻转血矛线虫敏感虫株和耐药虫株显著差异表达的miRNA和其对应的靶基因,并对部分耐药相关通路中存在的靶基因及显著差异的靶基因所在的通路进行筛选和分析,为后续进一步通过转录组学深层次挖掘捻转血矛线虫产生耐药性的关键基因及相关调控分子提供科学依据。  相似文献   

9.
【目的】基于转录组学对西方蜜蜂工蜂不同虫态间的差异表达基因(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种信号通路在西方蜜蜂工蜂幼虫、蛹和成虫的发育调控中发挥着重要作用。  相似文献   

10.
C2H2型锌指转录因子在植物、动物和部分真菌的生长发育过程中发挥着重要的作用,但其在偏肿革裥菌(Lenzites gibbosa)中的作用尚不清楚.通过对木屑处理后不同时期的偏肿革裥菌菌丝进行提取和转录组分析,为偏肿革裥菌在营养基质条件下菌丝的生长研究提供支持.运用COG、GO和KEGG数据库进行注释分析,筛选出与偏肿革裥菌生长发育相关的C2H2型锌指基因,并对差异表达基因进行GO Pathway富集分析以及实时荧光定量PCR分析.COG注释分析表明,C2H2型锌指蛋白差异基因大部分与复制、重组和修复功能有关.GO富集分析表明,C2H2型锌指蛋白差异基因参与了细胞过程、细胞部分和结合功能.将差异基因比对到KEGG数据库中,共得到35条注释的序列,其中4条初生代谢通路得到注释,对其进行Pathway富集分析,发现差异基因显著富集于内质网蛋白加工通路中,表明偏肿革裥菌的生长发育与内质网蛋白加工通路密切相关.在木屑诱导下,随着木屑诱导时间的延长,9406基因相对表达量呈逐渐增长的趋势,11332、7835和5947基因相对表达量呈"先升后降"的趋势,11819和4346基因相对表达量呈"升-降-升"的趋势,3913和5738基因相对表达量呈"降-升-降"的趋势.此研究结果对偏肿革裥菌菌丝体生长研究具有重要的参考价值,为深入研究木质纤维降解关键基因的功能与作用机制奠定基础.  相似文献   

11.
为探明骆驼斯氏副柔线虫雌虫和雄虫的转录组差异,挖掘出与性别决定相关的功能基因,采用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种雌性性别特异表达基因,这些基因在虫体性别决定和性别特异性表达中发挥重要作用。本研究利用生物信息学技术预测出雌虫和雄虫差异表达基因在功能分类和代谢通路等方面的生物学特征,挖掘出性别决定基因和性别特异性表达相关基因,为进一步开展斯氏副柔线虫的功能基因组学研究、药物靶点预测和疫苗候选基因筛选提供有价值的数据资源。  相似文献   

12.
为解析日本医蛭(Hirudo nippnica Whitman)繁殖的分子机理,研究不同温度下饲养的日本医蛭各组织基因表达差异情况。利用转录组学测序检测不同温度饲养条件下日本医蛭肌肉、环带和头部组织的基因表达情况,筛选出11个差异表达基因,并对其进行qPCR验证。qPCR结果显示一种具有去饱和酶功能的膜蛋白编码基因G6,在日本医蛭环带中的表达随温度的升高而降低,与转录测序结果一致,暗示G6基因可能为日本医蛭繁殖相关的关键基因。进一步对G6蛋白进行功能域预测,结果显示其属于一种跨膜蛋白。其他10个差异表达基因的转录组测序结果和qPCR结果基本保持一致,通过蛋白功能域预测结果显示预测出假定蛋白均含有的功能结构域;通过miRNA亲和预测,非编码RNA(G13)具有高亲和miRNAs特性。研究发现了10个潜在的功能基因和一个非编码RNA基因。  相似文献   

13.
14.
为了探究生物炭添加前后枯草芽孢杆菌应对镉胁迫的分子机制,本研究基于转录组测序,对比分析了枯草芽孢杆菌在4种体系(正常体系、镉胁迫体系、生物炭体系、生物炭-镉胁迫体系)中的基因表达差异。结果表明:在镉胁迫压力下,枯草芽孢杆菌共产生1 932个差异表达基因,其中上调基因967个,下调基因965个。枯草芽孢杆菌的代谢活性受到明显抑制,但是其鞭毛组装和细菌趋化性通路明显增强,离子转运和阳离子跨膜转运蛋白功能条目的基因表达也呈现上调,说明枯草芽孢杆菌通过上调细胞运动及离子转运相关基因的表达来规避镉毒性压力。相比于镉胁迫体系,在生物炭-镉胁迫体系中,细菌共产生1 541个差异表达基因,其中上调基因691个,下调基因850个。生物炭的存在增强了枯草芽孢杆菌生物过程、细胞组分和分子功能基因的表达,枯草芽孢杆菌物质和能源代谢相关通路的基因表达均呈现上调,并且细菌体内镉抗性基因czcD、cadA和膜功能相关差异基因的表达也以上调为主,说明生物炭可以缓解镉胁迫对枯草芽孢杆菌的毒性效应,并可刺激枯草芽孢杆菌的金属抗性能力。  相似文献   

15.
为探讨花魔芋抵御软腐病害的分子机制,分别以染病初期和健康的花魔芋球茎为材料进行转录组测序。结果表明,健康和染病初期花魔芋两组出现3 222个差异表达基因,其中2 660个基因上调表达,562个基因下调表达,差异表达基因主要涉及细胞组分、生物学过程及分子功能等生理生化过程。表达量上调和下调最大的前10个基因包含编码MiAMP1抗菌蛋白、热激蛋白、胰蛋白酶与蛋白酶抑制剂、胰凝乳蛋白酶抑制剂等与抗病相关的基因。GO功能分类和KEGG通路分析表明,类黄酮和苯丙烷类物质在花魔芋抵御软腐病菌侵染初期发挥重要作用。值得关注的是,硒化合物代谢、维生素B6代谢、类黄酮生物合成、N-聚糖生物合成及链霉素生物合成通路等抗病相关的差异表达基因在花魔芋染病初期全部或大部分受诱导表达。利用高通量测序获得大量花魔芋转录组信息,有助于挖掘与花魔芋抗软腐病相关的关键基因,也为魔芋抗病分子育种提供理论参考。  相似文献   

16.
旨在通过miRNA测序技术对捻转血矛线虫阿苯达唑给药前后的耐药虫株差异表达miRNA及网络调控分析,获取其转录本中全miRNA表达谱及耐药相关差异表达基因和miRNA,探究miRNA在捻转血矛线虫阿苯达唑耐药虫株给药前后的差异。采用Illumina Hiseq2 000平台对捻转血矛线虫阿苯达唑耐药虫株给药前和给药后的虫体进行Small RNA-seq测序,将Raw Data与参考基因组和miRBase数据库进行比对分析,使用RNAhybrid、Miranda和TargetScan对miRNA靶基因进行预测,筛选捻转血矛线虫耐药虫株给药前后表达量显著差异的miRNA并通过Stem-loop qRT-PCR验证,同时对差异显著的miRNA-mRNA调控关系采用Cytoscape软件进行可视化分析。将靶基因向KEGG和GO两大权威数据库映射,预测和分析靶基因功能注释情况和参与调控的通路情况。结果显示:1)捻转血矛线虫阿苯达唑耐药虫株给药前和给药后共筛选出显著差异表达的miRNA 188个,其中125个上调,63个下调,且qRT-PCR结果与转录组测序结果表达情况一致。2)显著差异表达的144个miRNA共关联到353个miRNA-mRNA负调控关系,注释到46条GO功能条目以及代谢,运输和降解等327条通路。3)筛选出的候选靶基因主要富集在ABC转运蛋白,HIF-1信号通路,cGMP - PKG 信号通路以及P450药物代谢通路,证实在药物胁迫的影响下捻转血矛线虫主要通过代谢和转运的方式参与耐药性调控。本研究为捻转血矛线虫耐药性诊断辅助分子标记的筛选和虫体耐药性调控因素及作用机制的研究提供参考。  相似文献   

17.
为明确截形叶螨Tetranychus truncatus雌成螨应对热胁迫的差异表达基因,采用Illumina HiSeq 4000进行转录组测序。结果表明,截形叶螨雌成螨热胁迫后共检测到17 577个差异基因,其中12 831个上调,4 746个下调。GO功能富集发现,雌成螨热胁迫后的差异基因主要富集于细胞过程、催化活性和细胞,KEGG通路富集分析发现,差异基因主要富集在代谢途径和溶酶体等。从转录组测序数据中筛选出与高温存在关联的抗氧化酶和热激蛋白等基因。采用RT-qPCR技术检测这些基因在截形叶螨雌成螨上的表达特征,发现17条基因的上下调趋势与转录组测序结果一致。这为后续深入研究截形叶螨与热胁迫相关基因的功能奠定了基础。  相似文献   

18.
为解析宁夏枸杞响应NaCl胁迫的ABA代谢分子调控机制,检测NaCl胁迫下宁夏枸杞叶片脱落酸(ABA)含量变化,并利用RNA-seq和qRT-PCR技术研究NaCl胁迫下ABA代谢相关基因的差异表达规律。结果表明,不同浓度NaCl胁迫下,宁夏枸杞叶片中ABA含量随NaCl浓度的增加呈现增加趋势;通过转录组测序筛选得到21个ABA代谢相关的差异表达基因,从100 mmol/L、200 mmol/L、300 mmol/L NaCl处理下的宁夏枸杞叶片中检测到ABA代谢相关的差异表达基因分别有17、11、9个,qRT-PCR检测结果与转录组测序结果基本一致。宁夏枸杞通过相关基因的差异表达调控ABA代谢和信号通路,从而参与调控其响应NaCl胁迫。  相似文献   

19.
山麻鸭开产期和产蛋高峰期卵巢组织转录组分析   总被引:3,自引:3,他引:0  
【目的】探明山麻鸭开产期和产蛋高峰期卵巢组织的转录组差异,丰富蛋鸭转录组数据信息。【方法】选取6份山麻鸭卵巢组织样品(开产期和产蛋高峰期各3份),利用Illumina HiSeqTM 2000进行高通量测序,构建山麻鸭开产期和产蛋高峰期卵巢组织转录组文库,并用测序评估、基因功能注释等生物信息学方法进行分析。【结果】经过测序获得质量不低于20的碱基比例(Q20)均高于93%;开产期获得了43 489 798条reads,4 380 847 061 bp数据量;高峰期获得了42 782 676条reads,4 307 499 083 bp数据量;分别有70.92%和72.70%的reads能比对到鸭参考基因组序列上。对开产期与产蛋高峰期转录组进行比较,发现开产期有26 808个表达基因,产蛋高峰期有27 013个表达基因,与开产期为参考,在高峰期卵巢组织中共获得1 929个差异表达基因,其中上调基因989个,下调基因940个;进一步将这些差异基因在Nr数据库进行注释,发现获得注释的基因有1 423个,其中上调基因695个,下调基因728个;COG功能注释分析发现这些差异表达基因共获得663个功能注释,涉及25个功能分类;GO 功能注释分析表明,有1 122个基因获得GO 功能注释,可以分为61个功能分类,这些分类主要涉及到分子绑定、催化活性、细胞过程、生物调节等诸多生理生化过程,其中涉及到发育繁殖生物学过程的相关注释基因有406个;KEGG分析发现共有425个基因注释到160个代谢通路中,其中GnRH信号通路、GPI-anchor生物合成、TGF-β信号通路、核糖体生物合成等通路显著富集。GnRH信号通路和TGF-β信号通路在卵巢的生理生殖活动发挥了重要作用,而GPI-anchor生物合成和核糖体生物合成,这两类代谢通路在卵巢生理生殖活动的作用尚不明确。【结论】利用高通量测序技术对山麻鸭开产期和产蛋高峰期卵巢组织的转录组进行测序和分析,揭示了山麻鸭不同生理状态下卵巢组织差异表达基因的数量,获得了差异表达基因的功能、分类和代谢通路。为丰富蛋鸭卵巢组织转录组信息,为开展蛋鸭产蛋性状相关基因的研究及分子调控机制研究奠定基础。  相似文献   

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

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