芋不同组织淀粉生物合成和代谢相关基因的差异表达分析

【目的】开展芋（Colocasia esculenta）不同组织的转录组试验,为芋基因功能和代谢调控网络的分子机制研究提供基因资源。【方法】采集播种90 d的荔浦芋1号芋球茎、叶片、叶柄和根进行生理检测和转录组测序,对不同组织的差异表达基因（DEGs）进行GO和KEGG功能富集分析,并对淀粉生物合成和代谢相关基因进行发掘和表达模式分析。【结果】芋球茎的淀粉和可溶性糖含量分别为134.85和23.92 mg/g,在4种组织中含量最高。转录组测序共获得85.41 Gb Clean data,与参考基因组比对效率为90.30%～92.56%,共33828个基因获得功能注释。不同组织间DEGs数量为5625～10562个。球茎与根、叶片和叶柄的比较组中,在生物过程中分别富集到与淀粉生物合成和代谢相关的GO功能2、5和5个;KEGG富集分析显示,在淀粉和蔗糖代谢途径分别富集到146、161和126个DEGs。通过富集分析发掘涉及淀粉生物合成和代谢相关DEGs共82个,其在4种组织中呈不同表达模式,其中ADP-葡萄糖焦磷酸化酶基因（AGPase基因）（Taro＿004018和Taro＿026553...     

小麦幼苗根系应答重金属Pb胁迫的转录组分析

为了探究小麦对重金属铅(Pb)胁迫应答的分子机制,采用水培法对小麦品种矮抗58进行不同质量浓度[0(对照)、40、80、160 mg/L]Pb胁迫处理,并对幼苗根系进行转录组测序,筛选分析差异表达基因,进行GO分类及KEGG富集分析。结果表明,不同质量浓度Pb处理下,小麦幼苗的根长和根数均受到抑制,且随着Pb质量浓度升高抑制作用增强。在Pb胁迫处理与CK间共获得了38 904个差异表达基因。其中,40、80、160 mg/L Pb条件下分别有6 072、16 581、16 251个差异表达基因。选取80 mg/L Pb处理的差异表达基因作为研究重点分别进行GO分类和KEGG通路富集分析。GO分类发现,上调差异表达基因主要富集在免疫系统过程、运动过程、代谢过程、节律性过程及催化活性和电子载体活性等方面,下调差异表达基因主要富集在发育过程、生长过程、定位过程、复制过程、生殖过程及转运活性和抗氧化活性等方面;KEGG富集分析发现,上调差异表达基因主要涉及植物激素信号转导途径、植物-病原互作途径、药物代谢途径、MAPK信号通路等方面,下调差异表达基因主要涉及次生代谢物的生物合成途径、苯丙烷类生物合成途径、抗生素生物合成途径、碳代谢和蔗糖代谢等方面。选取6个响应Pb胁迫的差异表达基因进行RT-PCR验证,结果表明6个差异表达基因的表达模式与RNA-Seq分析结果一致,进一步验证了RNA-Seq结果的准确性。     

酸雨胁迫下琯溪蜜柚叶片转录组差异表达分析

【目的】研究琯溪蜜柚受酸雨胁迫的内在分子机制,为琯溪蜜柚科学种植提供基础资料,也为酸雨逆境生理提供理论基础。【方法】以模拟酸雨胁迫24 h的琯溪蜜柚叶片进行Illumina HiSeq TM 4000高通量转录组测序分析,将组装得到的基因在参考基因组、Nr和KEGG数据库进行比对;利用FDR与log2(FC)来筛选差异基因,筛选条件为FDR0.05且|log2(FC)|1,将筛选的差异基因做GO和KEGG富集分析。【结果】模拟酸雨喷淋24 h后,琯溪蜜柚嫩叶出现明显块状伤斑;共得到21 497个基因,全部得到注释,与对照相比,酸雨处理组中有879个基因显著上调,588个基因显著下调;筛选出样本中前50个DEGs(差异基因),均为上调表达基因,大部分涉及代谢途径、次生代谢、苯丙基类丙烷生物合成和萜类物质合成等相关基因。GO富集分析表明差异表达基因主要位于细胞外区域;执行分子功能中催化剂活性是最为显著富集的GO term,其次是氧化还原酶活性;生物过程中差异表达基因最显著的GO term是DNA代谢过程。KEGG富集分析表明DNA复制是差异表达基因中最显著富集的Pathway,其次是次生代谢的生物合成,再次是苯丙素类合成途径。对次生代谢合成途径中4个差异表达基因[POD同工酶cg3g018770和cg2g001440、肉桂酰辅酶A还原酶(CCR)同工酶cg1g021310、4-香豆酸-辅酶a连接酶(4CL)同工酶cg3g029290]进行PCR荧光定量分析,验证了转录组数据的可靠性。【结论】琯溪蜜柚对酸雨耐性较强,对酸雨胁迫的响应是多基因参与、多生物过程协同调控的过程,次生代谢的调节可能是应对酸雨胁迫的主要方式。     

泡桐丛枝病相关lncRNAs的表达和调控

为揭示与泡桐丛枝病相关的长链非编码RNA(lncRNA),探讨其调控模式,利用高通量测序平台Illumina HiseqTM2000对泡桐的健康苗、丛枝病苗和利福平处理的丛枝病苗进行测序,筛选出差异表达的lncRNAs和靶基因,并对靶基因进行基因本体论(GO)功能及京都基因与基因组百科全书(KEGG)通路富集性分析。结果表明:共鉴定出的3 237个候选lncRNAs,其中PT/PTI、PT/PTIL-100和PTI/PTIL-100中差异表达并且相同的lncRNAs有147个,靶基因富集于40个GO和119条代谢通路。通过比较分析,发现72个lncRNAs可能与泡桐丛枝病发生相关,其靶基因主要参与的KEGG代谢途径有磷脂酶D信号通路、Wnt信号通路、ABC转运、糖基磷脂酰肌醇(GPI)-锚定生物合成途径、次生代谢产物的生物合成、泛素介导的蛋白水解、玉米素生物合成、甘氨酸、丝氨酸和苏氨酸代谢等。预测lncRNAs(TCONS_00009507、TCONS_00019521、TCONS_00012917、TCONS_00030963和TCONS_0002705)调控靶基因参与泡桐丛枝植原体与泡桐互作和影响细胞周期改变泡桐形态建成,为泡桐丛枝病发生机理的研究奠定基础。     

小麦响应禾谷镰刀菌侵染的核心基因筛选与鉴定

【目的】筛选鉴定小麦品种连麦12响应禾谷镰刀菌侵染过程中的核心基因,为小麦赤霉病防御分子网络提供理论参考。【方法】对成功感染禾谷镰刀菌的连麦12小穗进行转录组测序,通过DESeq筛选差异表达基因（DEGs）,对DEGs进行GO功能注释和KEGG信号通路富集分析,构建DEGs的蛋白质互作网络,通过MCODE和cytoHubba筛选核心基因。【结果】对小麦品种连麦12接种禾谷镰刀菌（处理组）和接种无菌去离子水（对照,CK）72 h后的DEGs进行筛选,结果共筛选出14180个DEGs,其中有10966个DEGs上调表达,有3214个DEGs下调表达。经GO功能注释分成三大类,共51个条目,其中,细胞组分包含14个条目,主要富集在细胞、细胞部件和膜3个条目;分子功能包含11个条目,主要富集在催化活性和结合2个条目;生物学过程包含26个条目,主要富集在细胞过程、代谢过程和刺激反应3个条目。KEGG信号通路富集分析结果显示,富集到代谢、次级代谢物生物合成、苯丙烷类生物合成、植物信号激素转导、植物MAPK信号通路的DEGs数量最多;富集到光合作用—天线蛋白通路,苯丙氨酸、酪氨酸和色氨酸生物合成,谷...     

波纹龙虾对亚硝酸盐胁迫响应的转录组学分析

为了探究波纹龙虾Panulirus homarus对亚硝酸盐的胁迫响应机理,采用转录组学测序分析方法,进行了正常养殖条件下波纹龙虾(体质量为39.00 g±5.31 g)肝胰腺和质量浓度为80 mg/L亚硝酸盐胁迫7 d的肝胰腺转录组比较研究。结果表明：转录组共筛选得到264个差异表达基因,其中,86个基因上调表达,178个基因下调表达;GO功能注释显示,被注释的差异表达基因主要与结合、催化和代谢等功能有关;KEGG通路富集分析显示,差异表达基因在黏着斑、白细胞跨内皮层迁移、胰岛素抵抗、紧密连接、淀粉与蔗糖代谢等通路中显著富集;随机挑选9个差异表达基因进行qRT-PCR验证,发现荧光定量结果与测序结果基本一致。研究表明,在亚硝酸盐胁迫下,通过GO和KEGG分析发现差异表达基因主要富集在能量代谢、生长和免疫等通路中,挖掘出的数据可为后续研究提供科学参考。     

秋茄(Kandelia candel L.)响应重金属胁迫的数字基因表达谱

基于高通量测序的数字基因表达谱(DGE)技术,对重金属污水胁迫处理组和对照组秋茄材料进行转录组测序分析。结果表明,重金属污水胁迫处理样品与对照样品相比,共筛选出的3 097个差异表达基因,其中2 202个表达上调,895个表达下调。GO功能显著性富集分析出38个GO分类条目,大量与生长发育,代谢调控,环境刺激响应相关的基因表现为富集。进一步的KEGG代谢途径分析表明,差异表达基因主要与糖类、核酸、蛋白质和脂质等生物大分子代谢、能量代谢以及次生代谢过程有关。转录因子分析发现b HLH,MYB,NAC和WRKY在重金属胁迫中发挥重要作用。此外,重金属胁迫还促进萜类、黄酮类化合物生物合成的关键酶基因(牻牛儿基焦磷酸合酶(Unigene0029352)、黄酮醇合成酶(Unigene0010384))的表达,进而促进秋茄有效成分的积累;显著诱导细胞分裂素相关基因type-b响应调节子ARR2(Unigene0033740)的表达、抑制油菜素内酯的信号转导组分基因BSK(Unigene0008986)的表达,进而提高秋茄对重金属胁迫的适应能力。实验选取了5个与环境刺激响应密切相关的基因,通过qRT-PCR验证了它们在重金属污水胁迫处理下的表达差异,结果与数字基因表达谱分析的结果较为一致,证实了差异表达基因数据的有效性。说明在重金属胁迫处理下,秋茄通过调节基因表达提高自身对污染胁迫的耐受能力。     

基于转录组测序挖掘小麦叶片抗旱相关基因

【目的】通过转录组测序分析干旱胁迫的小麦,以期挖掘小麦响应干旱胁迫的关键基因。 【方法】采用高通量 Illumina Hiseq 2500 测序平台对郑麦 366 号正常生长与自然干旱处理进行测序,并对测序数 据进行分析,筛选响应干旱胁迫的关键基因。【结果】两组样品中分别获得 45 471 018、45 133 470 个 clean read 片段,包含 6.82 Gb、6.39 Gb 碱基信息,拼接组装后,得到 188 334 个转录本和 119 588 个 Unigenes。通过筛选 共获得 1 207 个 DEGs,其中 752 个上调,455 个下调。GO 功能富集分析显示,上调的 DEGs 主要富集在水解酶 活性、催化活性、代谢过程、甘油代谢过程、膜和细胞部分等;下调的 DEGs 富集在细胞器、细胞和氧化还原 酶活性的最多。KEGG 途径富集分析发现,DEGs 显著富集于淀粉和蔗糖代谢、卟啉和叶绿素代谢、光合作用 - 天线蛋白、光合作用、苯丙烷生物合成、乙醛酸和二羧酸代谢、甘油磷脂代谢、氨基酸代谢和光合生物中的碳 固定等代谢通路。5 个抗氧化酶基因的 qRT-PCR 结果与 RNA-Seq 测序结果变化趋势相一致。【结论】利用转录 组测序技术获得了大量小麦抗旱相关基因,为深入研究小麦响应干旱胁迫的分子机制提供更多的基因遗传资源。     

新疆野苹果幼苗叶片响应NaCl胁迫的miRNAs及靶基因筛选

【目的】 筛选新疆野苹果响应NaCl胁迫相关差异miRNA及靶基因。【方法】 以新疆野苹果组培苗为材料,对150 mmol/L NaCl处理6和48 h时的新疆野苹果幼苗叶片进行small RNA测序。【结果】 NaCl处理6 h时3个miRNA差异表达,其中2个上调表达,1个下调表达,miRNA对应的靶基因数目为64个;NaCl处理48 h时13个miRNA差异表达,其中4个上调表达,9个下调表达,miRNA对应的靶基因数目为108个。对差异靶基因进行GO和KEGG分类注释,NaCl处理6 h时,GO主要富集包括:肌醇磷酸2激酶活性、寡糖基转移酶活性、蛋白丝氨酸/苏氨酸磷酸酶抑制剂活性、木聚糖生物合成过程等;NaCl处理48 h时,GO主要富集包括:氧氧化还原酶活性、DNA结合、质外体、次生代谢过程等。KEGG共同富集在N-glycan生物合成通路。【结论】 mdm-miR168b、mdm-miR159c、mdm-miR827、mdm-miR390f、mdm-miR171i、mdm-miR399j基因在NaCl处理6和48 h时表达量存在差异,其中mdm-miR390f及其靶基因HF10976、mdm-miR171i及其靶基因HF33844、mdm-miR399j及其靶基因HF11095表达量呈负相关,3个miRNA参与了新疆野苹果对NaCl胁迫的响应过程。     

镉胁迫下姬松茸菌株转录组测序及分析（英文）

为探明姬松茸菌株对重金属镉的富集机理,以前期筛选出的镉高、低富集姬松茸菌株（J11、J4）为材料进行2 mg/L的镉胁迫处理,以未添加外源镉作对照,采用Illumina高通量测序技术对其菌丝样品进行转录组测序分析。结果表明：差异表达基因主要富集在类固醇生物合成、抗生素的生物合成、碳代谢、谷胱甘肽代谢、内质网的蛋白质加工、甲烷代谢、淀粉和蔗糖的代谢等通路;其中,类固醇生物合成、抗生素的生物合成途径上调;甲烷代谢、淀粉和蔗糖代谢、碳代谢、谷胱甘肽、乙醛酸和二羧酸等代谢途径下调。分析认为碳代谢、谷胱甘肽代谢可能在姬松茸菌株响应镉胁迫方面起重要的作用。     

转录组分析二马牙和长脖类型林下参表型差异

为分析二马牙和长脖类型林下参不同组织中人参皂苷含量存在差异的分子机制，测量同一环境下14年生两种林下参表型及人参皂苷含量，利用高通量测序技术对两种类型林下参的叶片、芦头及主根进行转录组测序，筛选差异表达基因，并对芦头的7个差异基因进行了qRT-PCR验证。表型分析结果表明，与长脖相比，二马牙芦头短粗、须根发达、根部及地上部分均更粗壮、产量高。转录组分析结果显示，2种类型林下参叶片、芦头、主根中分别有124、604、89个差异基因，上调基因分别有48、283、31个，下调基因分别有76、321、58个；通过KEGG代谢通路富集分析，差异基因分别富集到16、70、17个通路中，主要注释到淀粉和蔗糖代谢、植物激素信号转导通路以及苯乙醇苷生物合成通路中。研究结果为林下参的品种选择以及揭示林下参表型差异的分子机制提供理论参考。     

凡纳滨对虾转录组测序分析及肌肉生长发育相关基因的筛选

【目的】筛选出凡纳滨对虾生长发育的功能基因及代谢调控网络,揭示其生长发育的分子机制,为后续开展凡纳滨对虾分子生物学研究提供宝贵的基因数据来源。【方法】以快速生长群体和慢速生长群体的凡纳滨对虾肌肉组织为研究材料,通过Illumina HiSeqTM2500平台对构建的cDNA文库进行高通量测序分析,以StringTie进行拼接组装后利用DESeq2筛选差异表达基因,并基于KOG、GO、nr、COG、Swiss-Pro、KEGG和Pfam等数据库进行差异表达基因功能注释分析。【结果】经拼接组装共获得53458844条Clean reads,各样品Clean reads的Q30均在93.00%以上;Illumina测序获得的Clean reads与参考基因组的比对效率在87.75%~87.80%,说明转录组测序数据真实可靠。采用SnpEff进行SNP/InDel变异注释分析,结果显示,SNP位点中以A>G、G>A、C>T和T>C等4种类型的数量较多（14950~21562个）,InDel位点则以SYNONYMOUS_CODING的数量最多（33630个）。使用StringTie对Mapped reads（比对到参考基因组的Reads）进行拼接,共发掘到4607个新基因,分别输入COG、GO、KEGG、KOG、Pfam、Swiss-Prot、eggNOG和nr数据库中进行序列比对,最终发现共有1098个新基因被注释,以被nr数据库注释的新基因数量最多（1077个）,而被COG数据库注释的新基因数量最少（416个）。基于Q-value<0.05且Fold Change>2的筛选条件,共获得1408个差异表达基因（661个为显著上调表达基因,747个为显著下调表达基因）;1408个差异表达基因被注释到53个GO功能条目中,其中,22条被注释到生物学过程（Biological process）,16条被注释到细胞组分（Cellular component）,15条被注释到分子功能（Molecular function）;KEGG信号通路富集分析发现3条重要的信号通路,分别是溶酶体通路（Lysosome）、氨基糖和核苷酸糖新陈代谢（Amino sugar and nucleotide sugar metabolism）及鞘脂类代谢（Sphingolipid metabolism）。在溶酶体通路中,CTSL基因、Nramp基因、MyoG基因和Myf5基因在凡纳滨对虾快速生长群体肌肉组织中呈上调表达,而Trypsin基因呈下调表达。【结论】通过转录组测序分析从快速生长群体和慢速生长群体的凡纳滨对虾肌肉组织中筛选出1408个差异表达基因（661个为显著上调表达基因,747个为显著下调表达基因）,主要富集在溶酶体、氨基糖和核苷酸糖新陈代谢及鞘脂类代谢等通路上,在对虾肌肉生长发育过程中发挥重要作用。     

二倍体及其同源四倍体酸枣的生理特征和转录组分析

目的研究酸枣基因组复制后分子水平的变异机制,为深入探讨枣树多倍体性状、关键调控基因的克隆及基因工程育种提供参考。方法本实验以二倍体及其同源四倍体酸枣为研究材料,比较了二者的叶片相对含水量、叶绿素含量、可溶性糖含量以及可溶性蛋白质含量,对两种材料进行了转录组测序,并参考GO Ontology、KEGG等数据库对差异基因和转录因子进行功能分类与富集分析。结果四倍体的叶绿素含量、可溶性糖含量和可溶性蛋白质含量均显著高于二倍体植株,叶片相对含水量在两种材料中没有显著差异。二倍体与四倍体酸枣共1 329个基因具有显著差异,GO功能分析表明这些差异基因主要参与生长发育和胁迫耐受相关功能。KEGG通路分析显示,大部分差异基因富集在碳水化合物代谢、氨基酸代谢和信号传递过程,其中16个关键基因参与糖和氨基酸的代谢与转运过程,如SPS2、GAE6和PGDH3,这些基因在四倍体中具有较高的表达量;23个基因参与激素信号传导通路,其中与生长素传导、应答相关的基因如ARG7,GH3.6和IAA26在四倍体植株中表达量较高,而与油菜素内酯合成酶（CYP）、乙烯不敏感蛋白（EIN）相关的基因在四倍体植株中表达量较低。在对二倍体与四倍体差异转录因子的分析中发现,四倍体MYB转录因子家族基因表达量高于二倍体。结论四倍体叶片叶绿素含量高是植株叶色加深的原因之一,可溶性糖、可溶性蛋白质含量高为四倍体叶片变大、茎加粗提供了更多能源物质。参与糖、氨基酸代谢和激素合成、信号传导的关键基因在二倍体与四倍体中差异表达,可能与四倍体植株体内能源物质含量高、生长势强的性状相关;具有渗透调节功能的基因在四倍体中表达量较高暗示着四倍体可能具有较强的抗性。进一步对差异转录因子的分析表明,MYB转录因子在二倍体与四倍体植株中的差异表达可能会导致植物在生长发育、形态建成和抗逆过程中的差异,但二者的具体差异性状还需进一步的实验研究。      

初情期山羊胸腺microRNAs的筛选

为分析初情期前和初情期雌性山羊胸腺中MicroRNA的差异表达情况,并探讨其在山羊初情期启动过程中所起的作用.对初情期前(n = 3)和初情期(n = 3)山羊的胸腺组织进行了 Solexa测序,采用GO富集和KEGG分析评估差异miRNA靶基因.在初情期前和初情期的山羊胸腺组织中检测到538个miRNA,其中64个miRNA显著差异表达.与初情期前的山羊胸腺样本相比,初情期胸腺样本中存在19个上调基因和45个下调基因.KEGG通路富集分析显示,差异miRNA靶基因主要富集在苯丙氨酸,酪氨酸和色氨酸的生物合成,丁螺菌素和新霉素的生物合成等通路以及胆碱能突触和过氧化物酶体增殖物激活受体(PPAR)等与初情期启动相关信号途径.结果提示胸腺中miRNA参与山羊初情期启动.     

育成期高能饲喂下开产与未开产蛋鸡肝脏miRNA表达谱差异分析

【目的】对育成期高能饲喂下开产与未开产蛋鸡的肝脏进行miRNA高通量测序分析,探明高能饮食状态下蛋鸡肝脏中影响其开产的miRNA,为提高产蛋性能打下基础。【方法】以代谢能水平为12.14 MJ/kg的高能饲粮饲喂育成期蛋鸡,通过Illumina NextSeq500高通量测序平台对开产与未开产蛋鸡肝脏进行small RNA测序,使用DESeq分析miRNA表达量,并以实时荧光定量PCR进行验证;采用miRanda对差异表达miRNA进行靶基因和靶位点预测,同时以超几何分布对差异表达的靶基因进行GO功能注释分析和KEGG信号通路富集分析。【结果】 6个样本（3个开产蛋鸡样本,3个未开产蛋鸡样本）能注释到的miRNA均超过300个,约占miRBase中已鉴定鸡miRNA的30.00%,且前体鉴定结果显示各样本注释到的miRNA有部分定位在同一前体上。与未开产组样本相比,开产组样本有9个miRNA表达上调、3个miRNA表达下调。其中,gga-miR-132c-5p、gga-miR-132b-5p、gga-miR-2184a-5p、gga-miR-132c-3p、gga-miR-132b-3p及gga-miR-2184a-3p属于同一miRNA基因簇,且gga-miR-1682、gga-miR-132b-3p和gga-miR-2184a-3p等3个上调miRNA在未开产蛋鸡样本中不表达。通过miRanda共预测得到129个差异表达潜在靶基因,以miR-203a预测得到的靶基因最多,为32个;其次是gga-miR-2184a-5p、gga-miR-148a-5p和gga-miR-12211-5p,分别预测到30、25和23个靶基因;而miR-132c-5p预测得到的靶基因最少,仅为9个。129个靶基因显著注释到8条GO功能条目上,生物学过程（Biological process）主要注释到脂质代谢过程、细胞脂质代谢过程、脂质生物合成过程、磷脂生物合成过程、磷脂代谢过程、甘油磷脂生物合成过程和解剖学结构发育,分子功能（Molecular funcion）仅注释到1条GO功能条目,即利钠肽受体活性,未发现涉及细胞组分（Cellular component）的GO功能条目;在注释到的8条GO功能条目中有6条与脂质代谢相关,涉及的靶基因有22个,占总潜在靶基因的17.10%。KEGG信号通路富集分析结果显示共显著富集到7条KEGG信号通路,其中脂肪酸降解、脂肪酸生物合成、酮体合成与降解及PPAR信号通路等4条信号通路与脂质代谢相关。【结论】育成期高能饲喂下开产与未开产蛋鸡中共存在12个差异表达miRNA,涉及129个差异表达潜在靶基因,且主要富集在肝脏脂质代谢相关过程和信号通路,说明miRNA是通过调控脂质代谢及其相关基因表达而影响蛋鸡开产。     

谷子萌发吸水期关键代谢途径的筛选与分析

【目的】谷子适应性强,抗旱耐瘠,是起源于中国的重要作物。通过转录组测序技术分析谷子萌发不同吸水期的转录组差异,以期获得谷子萌发过程中的差异表达基因,寻找调控谷子萌发的重要代谢途径和代谢物。【方法】以晋谷20为材料,构建谷子萌发过程中开始快速吸水期、滞缓吸水期和重新大量吸水期的cDNA文库,进行转录组分析;采用K-Means开展基因表达聚类分析;利用DESeq筛选差异表达基因;通过COG、GO、KEGG等对差异表达基因进行功能注释;利用KEGG富集挖掘不同吸水期调控种子萌发的关键代谢途径和关键基因;并采用qRT-PCR验证其可靠性;用HPLC分析关键代谢物含量。【结果】转录物测序分析获得谷子萌发开始快速吸水期、滞缓吸水期和重新大量吸水期覆盖整个基因组的基因表达谱,共获得33 643个基因,识别9个具有不同表达模式的共表达基因簇。比较种子萌发的开始快速吸水期与滞缓吸水期、滞缓吸水期与重新大量吸水期、开始快速吸水期与重新大量吸水期,分别筛选出3 893、4 612和8 472个差异表达基因。KEGG富集分析表明,3个比较的差异表达基因都显著富集到phenylpropanoid biosynthesis、phenylalanine metabolism、starch and sucrose metabolism代谢途径;开始快速吸水期与滞缓吸水期、开始快速吸水期与重新大量吸水期的差异表达基因还显著富集到plant hormone signal transduction途径。并且3个比较中富集到phenylpropanoid biosynthesis和phenylalanine metabolism代谢途径的差异表达基因数都最多,其中过氧化物酶基因（peroxidase）比例最高。通过qRT-PCR对4个苯丙烷生物合成途径相关基因的分析表明,其表达趋势与转录组分析结果基本一致,其中,4-香豆酸-CoA连接酶3（4-coumarate-CoA ligase 3）在谷子种子中存在已形成mRNA,萌发吸水过程中呈先下调后上调再下调的表达趋势。苯丙烷类相关代谢物含量分析显示,芥子酸在种子中大量储备,在萌发过程中呈下调趋势;阿魏酸、对香豆酸和咖啡酸呈先上调后下调趋势。【结论】谷子萌发过程中,不同吸水期的差异表达基因显著与苯丙烷生物合成途径和苯丙氨酸代谢途径相关;其上游基因4-香豆酰-辅酶A连接酶和下游基因过氧化物酶家族成员在谷子萌发响应水分过程中发挥调控作用;中间产物芥子酸可能参与种子的休眠与萌发。     

槐猪和杜洛克猪背最长肌差异表达基因的筛选与注释

为了筛选槐猪和杜洛克猪背最长肌中的差异表达基因,对槐猪和杜洛克猪背最长肌进行转录组测序,并进行生物信息学分析.结果表明,槐猪和杜洛克猪分别获得36 177 844和34 496 450条可用读段,且与猪参考基因组的比对率分别为78.71%和80.09%.槐猪和杜洛克猪间共有953个差异表达基因,以杜洛克猪为对照组,槐猪与之相比,有589个上调基因和364个下调基因.槐猪和杜洛克猪的差异表达基因共分类到48个基因本体(GO)条目中,其中,细胞组分13个、分子功能12个和生物学过程23个.显著富集的KEGG通路共计6个,其中,PPAR信号通路是与调节脂类代谢显著相关的通路.     

桑葚幼果的落果与正常果的果柄转录组分析

  目的  果实的脱落是一种常见的生理现象，会在一定程度上限制果实的产量。探究桑葚幼果在脱落过程中的代谢及生化反应，可为桑葚的脱落机制研究提供参考。  方法  通过高通量测序对白桑Morus alba幼果在落果和正常果的果柄离区进行转录组测序，筛选差异表达基因，并对其进行功能注释和代谢通路分析，采用实时荧光定量PCR技术对转录组数据进行验证。  结果  经华大基因转录组测序分析共计获得262.92 Mb的原始序列；桑葚幼果果柄在2种状态下共筛选到10 481个差异表达基因，其中，有5 239个差异表达基因上调表达，5 242个差异表达基因下调表达。经基因本体数据库(GO)功能富集分析，共发现37个显著性条目，大多数差异基因具有催化活性、膜的组成成分、氧化还原酶活性、膜的内在成分等功能；京都基因与基金组百科全书(KEGG)富集分析发现:大多数差异表达基因集中在柠檬酸循环、植物激素信号转导途径、黄酮类生物合成等代谢通路上。  结论  桑葚在脱落过程中受植物激素的调控，氨基酸、黄酮类次生代谢物和碳水化合物等物质也能调控果实脱落。图9表2参35     

somatostatin1基因突变斑马鱼仔鱼转录组分析

生长抑素（Somatostatin,SST）是一类神经激素多肽,对脊椎动物的生长和代谢起着关键的调控作用。在斑马鱼（Danio rerio）和其它鱼类中,生长抑素由6个基因编码（sst1-sst6）,sst1存在于所有脊椎动物中,并且在进化上最为保守。本研究利用CRISPR/Cas9基因编辑技术在斑马鱼中构建了可稳定遗传的生长抑素基因sst1突变体鱼系。通过比较6 dpf（days post fertilization,dpf）的sst1-/-突变体和野生型仔鱼转录组,发现sst1-/-突变体相关生化过程和信号通路发生了显著变化。相比野生型,sst1-/-突变体仔鱼中有354个基因显著上调,504个基因显著下调。GO（Gene ontology）富集分析表明,大部分差异基因与氨基酸和蛋白质生物合成相关。KEGG（Kyoto encyclopedia of genes and genomes）通路富集显示代谢通路、氨基酸合成、氨基酰基- tRNA合成、核糖体合成和rRNA加工相关通路表达增强。结果表明,sst1-/-突变体蛋白质生物合成和代谢过程更加活跃。RT-qPCR验证结果显示,随机选取的10个基因mRNA表达变化趋势与转录组测序结果一致,说明转录组测序结果真实可靠。6月龄突变体体重体长和野生型相比并无显著差异。本研究结果为进一步挖掘生长抑素基因家族的潜在功能奠定基础。     

花魔芋响应软腐病菌侵染初期的转录组分析

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