山麻鸭开产期和产蛋高峰期卵巢组织转录组分析

【目的】探明山麻鸭开产期和产蛋高峰期卵巢组织的转录组差异,丰富蛋鸭转录组数据信息。【方法】选取6份山麻鸭卵巢组织样品（开产期和产蛋高峰期各3份）,利用Illumina HiSeq^TM 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生物合成和核糖体生物合成,这两类代谢通路在卵巢生理生殖活动的作用尚不明确。【结论】利用高通量测序技术对山麻鸭开产期和产蛋高峰期卵巢组织的转录组进行测序和分析,揭示了山麻鸭不同生理状态下卵巢组织差异表达基因的数量,获得了差异表达基因的功能、分类和代谢通路。为丰富蛋鸭卵巢组织转录组信息,为开展蛋鸭产蛋性状相关基因的研究及分子调控机制研究奠定基础。     

基于转录组测序的金钱蒲类黄酮生物合成基因的表达分析

【目的】高通量测序获取金钱蒲（Acorus gramineus）7个不同组织转录组信息，为不同组织中类黄酮化学成分合成差异提供分子信息，进而在分子水平上研究金钱蒲不同组织内类黄酮化学成分合成差异。【方法】利用高通量测序技术平台完成金钱蒲7个不同部位的转录组测序，对unigenes进行功能注释，借助注释结果挖掘类黄酮生物合成通路，筛选通路中的差异表达基因（Differentially expressed genes,DEGs）并进行表达量分析。【结果】共获得高质量数据39.91～42.97 M，总碱基量为5.98～6.45 Gb,Q30碱基百分比大于94.05%,GC含量为47.81%～50.32%。18 616条unigenes在GO数据库中获得62 506个注释，根据功能划分为细胞组分、分子功能及生物过程三大类，分别对应6、14、21个亚类，大量基因分布在细胞解剖实体、连接、催化活性和细胞过程等亚类中。10 124条unigenes富集在KEGG数据库的5大类19个亚类中，从其中筛选出4条类黄酮生物合成途径，14个关键酶，63个差异表达酶基因。这些基因在金钱蒲7个组织中具有表达差异性...     

基于银耳转录组测序的多糖代谢途径分析

【目的】对银耳进行转录组测序分析,分析多糖生物合成途径和挖掘关键功能基因。【方法】采用Illumina Hiseq 2500测序平台对银耳菌丝体进行测序,组装得到unigene序列,并与公共数据库进行对比分析和注释,挖掘果糖和甘露糖代谢通路相关基因。【结果】本研究共得到4.72 Gb raw data,拼接后获得17 008条unigene序列,N50为2073 bp;注释结果表明,76.47%序列能够注释到Uniprot公共数据库,注释到GO分类和KEGG通路的unigene分别为10 152和5671条。进一步分析银耳果糖和甘露糖代谢途径相关基因,挖掘得到74条unigene;其中编码L-iditol 2-dehydrogenase的unigene最多,其次是butanol dehydrogenase。【结论】基于转录组测序数据解析及预测,为研究银耳多糖和其它产物的生物合成途径及分子机制提供了数据支持,也为银耳品质的形成提供理论依据。     

贵紫麦1号籽粒色素形成相关基因的差异表达

【目的】探明贵紫麦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号籽粒花青素合成过程中作用显著。     

转录组测序法研究草珊瑚叶和根的基因差异表达

【目的】从基因表达的水平,初步分析草珊瑚叶和根之间次生代谢差异的分子机制,为二者之间临床疗效差异形成的分子机制分析提供信息。【方法】以福建省福州市的草珊瑚作为样品,采用Illumina HiSeq TM高通量测序技术测定草珊瑚叶和根的转录组,然后经过滤和Trinity组装,得到的unigenes再通过blast与Nr、Nt、Pfam、KOG、Swiss-Prot、Kegg和GO进行比对注释,并对叶和根的基因差异表达进行分析,尤其是对KEGG代谢通路富集的差异基因进行分析。【结果】转录组测序结果共获得0.4亿多个clean reads,经Trinity组装后共得到508 271个unigenes,其平均长度为740 bp,最大长度为17.3 kb。基于blast分析,共有148 561个unigenes在七大功能注释数据库中得到成功注释,占总基因数的58.80%。在分析基因表达水平差异时,发现草珊瑚叶和根的共同基因有93 127个,叶和根的差异基因分别为36 327个和52 268个;同时还发现在29 732种不同表达的unigenes中,有12 511个上调基因和17 221个下调基因;代谢相关KEGG具有显著差异的通路有淀粉和蔗糖代谢、苯丙烷类生物合成、乙醛酸和二羧酸代谢、光合生物的固碳作用、吞噬体、谷胱甘肽代谢、光合作用、丙氨酸、天冬氨酸和谷氨酸代谢、倍半萜类和三萜类生物合成、卟啉和叶绿素代谢、氮素代谢、昼夜节律-植物、光合作用-天线蛋白、芪类、二芳基庚酸和姜酚生物合成、不饱和脂肪酸生物合成、柠檬烯和蒎烯降解、类胡萝卜素生物合成、二萜类生物合成、类黄酮生物合成、脂肪酸延伸等。其中与药效密切相关的次生代谢通路苯丙烷类、倍半萜类和三萜类、二萜类、类黄酮类生物合成等途径分别有193个、82个、40个、35个差异表达基因,而上调倍半萜合酶、ent-kaur-16-烯合酶、黄酮醇合酶/黄烷酮3-羟化酶等基因和下调8-羟基香叶醇脱氢酶、vinorine合酶、角鲨烯合酶等关键酶基因差异显著。【结论】草珊瑚叶和根中苯丙烷类、倍半萜类和三萜类、二萜类、类黄酮次生代谢途径的相关基因差异最为显著,其中差异显著的关键酶基因可为分析其叶和根之间次生代谢差异的分子机制提供重要信息。     

红螯螯虾转录组高通量测序及分析

【目的】为发掘和研究红螯螯虾功能基因SSR标记奠定基础。【方法】对红螯螯虾肝脏、精巢以及卵巢转录组测序选用新一代高通量测序技术Illumina Hi Seq 2000,同时借助生物信息学法展开基因表达谱探究及功能基因预测。【结果】肝脏组织获得47 629 414条clean reads,精巢组织获得47 018 968条clean reads,卵巢组织获得53 267 362条clean reads。所有reads组装后获得了67 369个Unigene。通过GO分类,16 989个Unigene被分为3个种类,即生物学过程、细胞组分、分子功能。通过COG分类,4 697个Unigene被分为25个种类。通过KEGG分类,9 842个Unigene分属331个代谢途径。【结论】通过对红螯鳌虾肝脏等3个组织转录组的测序,获得了其相关基因的表达图谱及其功能分类。     

套袋对桃果皮生理生化指标的影响

摘 要：【研究目的】晚熟桃套袋是生产上常用技术，本研究通过对套袋后桃果实果皮生理生化指标的测定，为桃桃袋研究提供参考。【方法】以晚熟桃品种‘晚西菲’为试验材料，在幼果期对桃果实进行套袋处理，成熟期采收时测定桃果实果皮的可溶性糖、电导率、SOD、POD等生理生化指标。【结果】套袋处理后果皮的电导率、脯氨酸及可溶性糖含量有所降低；可溶性蛋白、丙二醛含量、超氧化物歧化酶（SOD）、过氧化物酶(POD)活性均有所增高。【结论】套袋环境对于晚西非桃果实发育是有利环境。     

套袋对桃果实品质影响的研究进展

本文简要综述了套袋对桃果实大小、着色、果面光洁度、果实硬度、可溶性固形物含量、可溶性糖、可滴定酸以及维生素C含量等果实品质影响的研究进展,以及套袋对桃果实花青苷、维生素C合成的影响机理,分析了桃套袋栽培存在的主要问题,并探讨了今后的研究与发展方向.     

荸荠优良品种‘桂蹄3号’球茎发育过程转录组测序分析

【目的】开展荸荠球茎发育过程转录组测序研究,为研究荸荠球茎发育过程相关基因表达信息提供参考。【方法】运用高通量测序技术,对荸荠球茎不同发育时期进行转录组测序研究。【结果】组装共得到223 182条转录本和90 542条Unigene,平均长度为809 bp,N50为1119。组装完整性较高,效果较好。对所得Unigene进行不同数据库注释,共有50 583条Unigene成功注释到7个数据库(NR、GO、NT、Pfam、KEGG、KOG及Swiss-prot)。对差异表达基因分析,发现球茎膨大初期的差异表达基因数目最多,为最活跃阶段。GO功能富集分析结果表明,33 205个基因获得功能注释,分为分子功能、细胞组分和生物学过程等3大类和54个亚类。COG功能分类结果表明,17 743个基因分布于25个功能区域,其中碳水化合物代谢占重要地位。KEGG代谢通路注释结果表明有20 667个基因获得功能注释,共有116条代谢途径,其中淀粉-蔗糖代谢占主要作用。【结论】利用高通量转录组测序技术首次建立了荸荠优良品种‘桂蹄3号’球茎的转录组数据库,为进一步研究荸荠球茎淀粉生物合成相关基因的功能及形成的分子机制提供了数据基础。     

‘仓方早生’桃及其早熟芽变不同发育时期果实的转录组分析

【目的】通过对桃品种‘仓方早生’及其早熟芽变不同发育时期的果实进行转录组分析，挖掘参与调控桃果实成熟的关键因子，为深入研究桃果实成熟调控机理提供理论依据。【方法】以桃品种‘仓方早生’及其早熟芽变为试材，每个品种分别选择长势一致的样品树5株，分别于花后30 d（对应‘仓方早生’c1、早熟芽变y1）、45 d（对应c2、y2）、59 d（对应c3、y3）、71 d（对应c4、y4）及89 d（对应c5）对不同发育时期的桃去皮果肉进行取样和转录组测序，并利用实时荧光定量PCR(qRT-PCR)对筛选的差异表达基因进行定量验证；利用GO和KEGG对‘仓方早生’及其早熟芽变的差异表达基因进行分析；基于差异表达基因构建加权基因共表达网络分析（weighted gene co-expression network analysis,WGCNA），从中鉴定出与果实成熟密切相关的枢纽模块和枢纽基因。【结果】将处于果实相同发育时期的转录组数据进行比较，得到y1与c1、y2与c2、y3与c4和y4与c5四组对比数据，共筛选出差异表达基因4 395个，其中上调表达基因2 212个，下调表达基因2 183个。其...     

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

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

半番鸭与番鸭精巢组织差异表达转录组测序分析

【目的】研究半番鸭与番鸭精巢组织转录组差异表达基因,为进一步阐明半番鸭不育的遗传机制奠定理论基础。【方法】利用转录组测序方法对半番鸭和番鸭的精巢组织进行研究,筛选其差异表达基因,并对其功能进行注释和荧光定量PCR(quantitative real-time PCR,QRT-PCR)验证。【结果】测序共获得43.84Gb Clean Data,组装后共获得193 535条Unigene。DESeq分析发现3 597个基因在两个鸭品种间差异表达,其中上调基因1 194个和下调基因2 403个,包括与生殖功能相关的基因,如成纤维细胞生长因子(fibroblast growth factor,FGF)、蛋白激酶A(protein kinase A,PKA)、丝裂原活化蛋白激酶7(mitogen-activated protein kinase 7-like,partial,BMK)、生长因子受体结合蛋白2(growth factor receptor-bound protein 2,GRB2)、c-Jun氨基末端激酶(c-Jun N-terminal kinase,JNK)和肿瘤坏死因子受体超家族成员6(tumor necrosis factor receptor superfamily member 6,FAS)等。GO(gene Ontology)分析发现382个差异基因获得功能注释,其中97个基因涉及发育繁殖生物学过程。KEGG(kyoto Encyclopedia of Genes and Genomes)通路分析表明差异表达基因共富集到50条信号通路中,其中17个通路显著富集,包括丝裂原活化蛋白激酶信号转导通路(mitogen-activated protein kinase signaling pathway,MAPK)、甘油酯代谢(glycerolipid metabolism)以及钙信号途径(calcium signaling pathway)信号通路等,与生殖功能密切相关的有促性腺激素释放激素信号通路(gonadotropin releasing hormone(Gn RH)signaling pathway)和MAPK信号转导通路。经QRT-PCR验证,差异基因表达变化模式与转录组测序结果一致,测序结果可靠。【结论】在转录组水平上筛选出半番鸭和番鸭精巢组织差异表达基因,揭示了Gn RH和MAPK信号通路在鸭的生殖活动中发挥了重要作用,为进一步探索半番鸭生殖系统的分化机理提供可靠的参考依据。     

基于转录组分析油菜素内酯对高温胁迫下酿酒葡萄花色苷合成及果实品质的调控机制

【目的】分析高温胁迫下参与油菜素内酯调控葡萄花色苷及果实品质合成的相关基因,探讨油菜素内酯调控果实花色苷及品质合成的机制。【方法】以酿酒葡萄‘赤霞珠’为试材,转色前一周利用红外辐射器模拟高温环境,并全树喷施0.6 mg·L^-1的2,4-表油菜素内酯(2,4-Epibrassinolide,EBR),测定花色苷、总糖及相关品质指标,选择转色中期(花后70 d)的果实进行转录组测序,从分子水平阐述EBR对高温胁迫下花色苷合成的影响。【结果】从转色开始,各处理花色苷含量逐渐升高;成熟时,高温组(HT)花色苷总量显著低于对照组(CK),高温油菜素内酯组(HTE)花色苷含量高于HT组。总糖、还原糖、蔗糖变化规律与花色苷相似,HT组含量均在成熟期时低于CK组,成熟期各种糖含量为CK组>HTE组>HT组。分析3种处理下‘赤霞珠’果实基因水平的差异,通过GO和KEGG富集发现了14个与蔗糖和淀粉代谢途经相关的差异基因,其中HT和HTE处理显著上调了10个基因,显著下调了4个基因;苯丙氨酸代谢途径有11个差异基因,其中有7个参与花色苷合成的基因在HT处理中上调,有4个参...     

拟轮枝镰孢与玉米籽粒互作的差异基因筛选及代谢通路分析

【目的】拟轮枝镰孢(Fusarium verticillioides)引起的玉米穗腐病是我国玉米产区发生严重的病害之一，论文旨在了解病原菌与玉米籽粒互作过程中的基因表达差异，为揭示病原菌的致病机制和玉米抗病机制提供依据。【方法】对拟轮枝镰孢侵染玉米籽粒0、4、12和72 h的样品进行转录组测序，之后采用生物信息学分析，分别以玉米和拟轮枝镰孢基因组为参考，以|log₂FC|≥1,P-adjust<0.05为阈值筛选互作过程中玉米和拟轮枝镰孢的差异表达基因，利用GO和KEGG对其进行功能注释及富集分析。Goatools软件分析植物-病原互作、MAPK途径和植物激素信号转导通路相关差异基因的表达变化，采用实时荧光定量PCR(qRT-PCR)方法对测序差异基因进行验证。【结果】在互作4、12和72 h后拟轮枝镰孢分别有140、400和1 945个基因上调表达，有9、302和1 784个基因下调表达；玉米分别有293、692和1 426个基因上调表达，320、482和153个基因下调表达。GO和KEGG富集分析显示，侵染早期拟轮枝镰孢在细胞间隙生长，差异基因主要富集在...     

不同葡萄糖/果糖类型桃在果实发育期间果实和叶片中可溶性糖含量变化及其相关关系

 【目的】研究不同葡萄糖/果糖（glucose/fructose，G/F）类型桃果实内G/F差异的部位和时期。【方法】以不同G/F类型的6个桃品种（G/F≈1品种：‘冈山白’、‘山一白桃’和‘燕红’；高G/F品种：‘张黄7号’、‘龙246’和‘临白7号’）为试材，采用高效液相色谱法测定果实发育期果实和叶片中可溶性糖含量，并在盛花后74 d或101 d测定了‘冈山白’、‘山一白桃’、‘张黄 7号’和‘龙 246’新梢韧皮部中可溶性糖的含量。【结果】两类不同G/F桃果实中均以蔗糖作为主要碳水化合物积累形式，花后43～85 d蔗糖含量很低，随后持续快速积累直至果实成熟；花后43～85 d山梨醇有升高趋势，在果实成熟前40 d左右迅速降低；葡萄糖和果糖含量在果实发育早期较高，之后逐渐降低；但两类不同G/F桃在整个果实发育过程中G/F值与果实成熟时相似。叶片中贮藏的可溶性糖主要是蔗糖和山梨醇，在果实整个发育期间，G/F≈1品种叶片中G/F约为1～3，而高G/F品种叶片中G/F约为2～7。G/F≈1品种‘冈山白’和‘山一白桃’与高G/F品种‘张黄7号’和‘龙246’韧皮部中山梨醇占总可溶性糖47%～58％，显著高于蔗糖、葡萄糖和果糖的含量，G/F为0.8～0.91，且两类不同G/F桃品种间G/F不存在显著差异。【结论】光合产物在韧皮部的运输对桃果实的G/F没有显著影响，果实中G/F的差异主要由于果实内糖代谢差异所导致。     

纽荷尔脐橙果实套袋品质分析

【目的】探讨果实套袋技术对纽荷尔脐橙果实外观质量和内在品质的影响,为贵州南亚热带低热河谷地区纽荷尔脐橙无公害栽培提供参考依据。【方法】以纽荷尔脐橙果实为材料,分别采用黄色单层纸袋和内层黑色双层纸袋进行套袋,对套袋果实的主要品质指标进行分析。【结果】纽荷尔脐橙果实套袋后,能促进果实生长发育,明显改善果实外观;果形指数与未套袋果实无显著差异;套袋后果实可溶性固形物、总糖、维生素C和酸含量均有所下降,除酸含量外,内层黑色双层纸袋与未套袋处理之间均存在极显著差异。【结论】纽荷尔脐橙进行套袋后能促进果实生长发育,改善果实外观品质,减少病虫为害,增加果品商品价值,减少投入成本,且黄色单层纸袋套袋效果优于内层黑色双层纸袋。因此,在贵州南亚热带低热河谷地区生产纽荷尔脐橙时,宜采用黄色单层纸袋进行果实套袋。