育成期高能饲喂下开产与未开产蛋鸡肝脏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是通过调控脂质代谢及其相关基因表达而影响蛋鸡开产。     

育成期高能饮食对蛋鸡肝脏转录组的影响

[目的]明确相关基因在育成期蛋鸡肝脏脂质代谢调控过程中的作用机理,为改善蛋鸡的生产性能打下基础.[方法]分别选取高能量(12.14 MJ/kg)和低能量(10.88 MJ/kg)饲料饲喂育成期蛋鸡,利用Illumina NextSeq 500平台对蛋鸡肝脏进行转录组测序,通过HTSeq和DESeq等生物信息学方法筛选出差异表达基因,并以超几何分布对差异表达基因进行GO功能富集分析和KEGG信号通路富集分析.[结果]从低能组和高能组蛋鸡肝脏样品RNA-seq测序获得的原始序列均超过8000万条,且能比对上基因组和基因的序列均在80.00%以上.其中,基因序列占87.32%~89.78%,基因间序列占10.22%~12.68%;外显子序列在基因序列中的所占比例很高,为78.10%~82.88%,内含子序列所占比例在17.12%~21.90%.相对于低能组,高能组蛋鸡肝脏中有82个基因上调表达、106个基因下调表达.差异表达基因GO功能富集分析发现有314个显著富集的GO功能(P<0.05,下同),其中,241个富集在生物过程(BP),18个富集在细胞组分(CC),55个富集在分子功能(MF).脂质代谢过程和细胞脂质代谢过程富集的基因分别有16和14个,占差异表达基因总数的8.51%和7.45%;脂质代谢过程较细胞脂质代谢过程多出的基因是PLIN1和FGF19.差异表达基因KEGG信号通路富集分析发现有4条显著富集的KEGG信号通路,且均与蛋鸡肝脏的脂质代谢相关.[结论]育成期对蛋鸡进行高能饮食饲喂,主要是通过上调或下调脂质代谢相关基因表达而影响肝脏脂质代谢,最终实现蛋鸡生产性能调控.     

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

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

高脂饲喂对肉鸡腹脂和肝脏组织RNA 可变剪接的影响

【目的】腹脂过度沉积会降低肉鸡抗病力和屠宰率。RNA可变剪接是影响基因功能的重要调控方式,对肉鸡腹脂沉积过程的RNA可变剪接事件进行系统分析,有助于进一步了解肉鸡腹脂沉积的分子调控机理。【方法】利用高脂饲喂肉鸡和普通饲喂肉鸡肝脏组织和腹脂组织的转录组测序结果,鉴定高脂饲喂过程中肉鸡腹脂和肝脏组织的差异表达和差异剪接基因。采用GO功能注释、KEGG通路富集分析和Metascape富集分析对显著差异剪接基因进行分析,并对显著差异剪接事件进行可视化处理,以及分析调控这些差异剪接事件的剪接因子。【结果】对高脂饲喂肉鸡与普通饲喂肉鸡的腹脂和肝脏组织进行转录组测序,在腹脂组织中发现233个显著差异表达基因和349个显著差异剪接基因,对显著差异剪接基因进行富集分析,发现显著差异剪接基因主要富集于细胞分裂、细胞生长等细胞增殖相关生物学进程,还富集于甘油磷脂代谢和胰岛素信号通路等脂肪生成相关通路。在肝脏组织中发现276个显著差异表达基因和224个显著差异剪接基因,对显著差异剪接基因进行富集分析,发现显著差异剪接基因主要富集于蛋白质乙酰化、蛋白质内部氨基酸乙酰化和m RNA加工等细胞代谢相关生物进程,还富...     

巴马小型猪高脂诱导模型建立及其肝脏转录组学研究

利用高脂饲料诱导建立巴马小型猪代谢疾病模型并分析肝脏转录组变化。选取1月龄巴马小型猪24头,随机分为对照组与高脂诱导组,对照组饲喂普通饲料,高脂诱导组饲喂含2.0%胆固醇饲料,定时监测体重、血常规及生化指标,6个月后试验结束,作心脏冠状动脉造影、器官组织学检测及肝脏转录组分析。结果表明,高脂诱导组个体饲喂4个月体重显著高于对照组,诱导期后血液生化指标中甘油三酯显著高于对照组,饲喂6个月时冠状动脉造影及器官组织切片未见显著病变,肝脏转录组分析显示高脂诱导组与对照组相比有2 468个差异表达基因,主要与代谢途径相关。为期6个月的高脂饲喂未致巴马小型猪明显代谢性疾病及心血管病变,但血液甘油三酯水平显著升高,显著影响肝脏代谢途径相关基因表达。文章系统探讨高脂诱导模型,结合血液生化、组织学、冠脉造影及肝脏转录组分析,进一步获得代谢通路相关潜在标记基因,为建立理想小型猪代谢疾病模型提供参考。     

基于RNA-Seq技术分析微囊藻毒素-LR对草鱼肝脏脂质代谢的影响

[目的]为探讨微囊藻毒素-LR(MC-LR)对草鱼肝脏脂质代谢的影响机制,对已获得的草鱼肝脏转录组数据进行了挖掘分析.[方法]草鱼经腹腔注射不同剂量(0,25,75,100μg/kg)MC-LR 96 h后,分离肝脏提取RNA,依托IlluminaHiseq-2500平台进行转录组测序分析,将筛选到的差异表达基因(DE...     

荷花转录组测序及花青素苷合成相关基因表达分析

以荷花不同花色品种‘青玉'和‘白洋淀红莲'为研究对象,运用高通量测序技术,对松蕾期的花瓣进行转录组测序,并对2个测序文库进行差异表达及实时荧光定量PCR分析。结果表明:共获得1 142条差异表达基因。筛选出与花青素苷合成相关的关键基因ANS、CHS、DFR、UF3GT。4个基因在两种荷花花瓣不同时期中的表达量存在显著差异,在松蕾期或初花期的表达量最高。     

芒果2个不同花芽分化时期转录组分析

[目的]分析芒果不同花芽分化时期的转录组,了解成花过程相关基因的表达情况,为芒果开花时间的分子调控机制研究提供理论参考.[方法]以芒果品种南逗迈4号为材料,采用Illumina高通量测序技术,分别对其新梢停长期(I期)和基部膨大期(II期)2个不同花芽分化时期顶芽进行转录组测序,并利用基因功能注释、差异表达基因筛选等生物信息学方法对测序获得的高质量序列进行分析.通过实时荧光定量PCR(qPCR)检测差异表达基因的表达情况以验证转录组测序结果的可信度.[结果]共获得85691条Unigenes,平均长度为870 bp,N50为1077 bp,与UniProt数据库比对,发现48589条Unigenes有同源信息,注释比例为56.7%,广泛涉及细胞组分、生物过程和分子功能三大类,共55个小类,其中参与生物过程的Unigene数量最多(有21个小类),以参与分子功能的Unigene数量最少(有14个小类).以Fold-Change≥2为条件,筛选出花芽分化Ⅰ和Ⅱ期转录组间的2031个差异表达基因,其中1073个基因表达上调,958个基因表达下调,涉及247条代谢通路,富集的KEGG通路为内质网蛋白加工、次生物质的生物合成和积累、糖代谢和光合作用等,其中注释为内质网蛋白加工的基因数量最多.从2个不同时期顶芽的转录组数据中获得了春化、光周期、赤霉素(GA)、成花抑制、自主和年龄等途径的开花相关基因.将MiSOC1、MiVIN3、MiDof和MiMADS1的qPCR检测结果与其转录组测序结果进行比对,发现这4个基因虽然在II期的表达量比Ⅰ期上调差异倍数上存在一定差异,但表达量变化趋势一致,说明转录组测序结果的可信度较高.[结论]芒果花芽分化过程与内质网蛋白加工、次生代谢生物合成、植物激素信号转导、淀粉和蔗糖代谢等途径密切相关,可为芒果花期人工调控和产期调节提供参考.     

利用高通量测序技术分析IHHNV感染凡纳滨对虾的基因差异表达

[目的]研究传染性皮下及造血组织坏死病毒（Infectious hypodermal andhematopoietic necrosis virus,IHH-NV）感染凡纳滨对虾的基因差异表达,为进一步了解IHHNV与凡纳滨对虾的相互作用分子机制提供参考.[方法]采用454高通量测序技术对IHHNV感染和对照凡纳滨对虾进行转录组测序,测序获得的序列除去接头后用iAssembler软件进行拼接得到unigene（单一基因序列）,通过计算各unigene的转录本数目（RPKM）确定差异表达基因,并以实时荧光定量PCR验证基因的差异表达.[结果]采用454高通量测序技术对IHHNV感染和对照凡纳滨对虾进行转录组测序,共获得208690条短序列,序列拼接共得到21912条tmigenes;通过对IHHNV感染和对照凡纳滨对虾的测序数据比较分析,结果发现221个差异表达基因,包括81个上调表达基因和140个下调表达基因.经实时荧光定量PCR验证,454高通量测序数据统计的基因表达差异结果可靠.[结论]IHHNV感染显著影响凡纳滨对虾免疫相关基因的表达,而利用高通量测序技术可有效检测出这些差异表达基因.     

火龙果果肉颜色相关MYB转录因子的表达分析

对红肉与白肉火龙果的果实进行转录组测序,通过数据的分析与筛选,得到30个差异表达的MYB类转录因子。运用生物信息学进行了MYB转录因子序列的比对、基因注释分析、进化树分析、基因表达情况分析等,通过qRT-PCR验证,探究红肉与白肉火龙果果实颜色的代谢差异与MYB转录因子之间的调控关系。结果表明:火龙果果肉颜色差异可能受MYB转录因子家族部分成员的调控。     

中华蜜蜂幼虫肠道响应球囊菌早期胁迫的转录组学

【目的】白垩病是困扰养蜂生产的顽疾。目前,尚无利用二代测序技术研究中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫白垩病的报道。本研究利用RNA-seq技术对健康(Ac CK)及球囊菌(Ascosphaera apis)胁迫的中蜂4日龄幼虫肠道(Ac T)进行深度测序,在转录组水平研究中蜂幼虫在球囊菌胁迫早期的胁迫应答。【方法】通过Illumina Hi Seq 2500平台对Ac CK和Ac T进行双端(PE125)测序,首先对测序数据进行质控和评估,利用edge R软件进行差异表达基因(DEG)分析,进而对DEGs进行GO富集分析及KEGG代谢通路富集分析,最后,利用实时荧光定量PCR(q RT-PCR)验证测序数据的可靠性。【结果】Ac CK和Ac T的转录组测序共得到188 457 338条原始读段(raw reads),经过滤得到182 088 448条有效读段(clean reads),两端Q20与Q30均在97.96%和94.97%以上,说明测序数据质量良好;主成分分析(PCA)结果显示第一与第二主成分可分别解释基因表达总体差异的75.8%和10.7%;DEG分析结果显示上调基因与下调基因的数量分别为344和239个;GO富集分析结果显示DEGs共富集在36个GO分类(term)上,其中基因富集数最多的细胞(106 unigenes)、细胞组件(106 unigenes)和代谢进程(104 unigenes);KEGG代谢通路富集分析结果显示上调与下调基因分别富集在72个和45个代谢通路上,其中上调基因富集数最多的是核糖体(72 unigenes)、碳代谢(16 unigenes)和糖酵解(14 unigenes),而下调基因富集数最多的是碳代谢(9 unigenes)、二羧酸代谢(8 unigenes)和氨基酸生物合成(7 unigenes),进一步分析表明中蜂幼虫肠道的部分细胞免疫对球囊菌的胁迫产生应答,而体液免疫不产生应答,宿主的代谢相关基因受到球囊菌的显著抑制。【结论】揭示了中蜂幼虫在球囊菌入侵早期的胁迫应答,为深入解析中蜂幼虫的胁迫应答机制提供了重要信息,也为在分子水平研究关键应答基因打下了基础。     

氮胁迫下小麦叶片转录组分析

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

氮胁迫下小麦叶片转录组分析

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

基于比较转录组学分析揭示中华蜜蜂及意大利蜜蜂幼虫的球囊菌抗性差异机制

为探究中华蜜蜂(Apis cerana cerana,以下简称中蜂)与意大利蜜蜂(Apis mellifera ligustica,以下简称意蜂)幼虫的球囊菌(Ascospahaera apis)抗性差异产生的原因,利用Venn分析、GO分类分析、KEGG代谢通路分析以及Ka/Ks分析对二者的转录组进行比较研究。Venn分析结果显示,中蜂与意蜂的同源基因有17 656个,归属于8 111个基因家族,二者的特有基因分别有1 158和241个,分别归属于468和86个基因家族。GO分类结果显示,中蜂和意蜂的特有基因分别富集在38和28个GO条目。KEGG代谢通路富集分析结果显示,中蜂和意蜂的特有基因分别富集于96和21个代谢通路;进一步分析发现中蜂的特有基因富集在内吞作用、溶酶体、泛素介导的蛋白水解和黑化作用等4个细胞免疫通路,以及MAPK信号通路和Toll-like受体信号通路等2个体液免疫通路,而意蜂仅有1个特有基因富集在MAPK信号通路。Ka/Ks分析结果显示受到强烈正向选择、弱正向选择、中性选择和负向选择的单拷贝同源基因分别有37、281、2 630和3 269个。对受到强烈正向选择的基因进行GO分类和KEGG代谢通路富集分析,GO结果显示中蜂和意蜂的单拷贝同源基因富集的GO条目类型相同;KEGG结果显示中蜂的单拷贝基因仅富集在氧化磷酸化。上述结果表明免疫相关基因数量的差异是二者的球囊菌抗性差异产生的重要原因之一,中蜂在与球囊菌的协同进化过程中可能通过提高能量利用从而限制球囊菌的增殖。本研究结果可为阐明中蜂及意蜂幼虫的球囊菌抗性差异产生的分子机制提供参考。     

截形叶螨雌成螨应对高温胁迫的转录组分析

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

低氧—复氧胁迫下脊尾白虾鳃组织差异表达基因趋势分析

【目的】探究脊尾白虾（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条,分别为核糖体、碳代谢、氧化磷酸化、氨基酸生物合成、内质网蛋白质加工、糖酵解/糖异生、谷胱甘肽代谢和蛋白输出。【结论】脊尾白虾鳃组织在受低氧胁迫早期通过合成蛋白质及提高代谢能力来抵御低氧环境,随着低氧胁迫时间的延长,物质合成和能量代谢活动均显著下降;但在复氧后随着复氧时间的延长,其蛋白质合成和能量代谢水平又逐渐升高恢复至常氧水平。     

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

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

卵形鲳鲹卵巢不同发育时期的转录组测序分析

【目的】鉴定筛选出与卵形鲳鲹卵巢发育相关的候选基因及信号通路,为揭示其卵巢性成熟过程的分子机制打下基础。【方法】挑选卵巢发育处于?期和Ш期的雌性卵形鲳鲹,分别构建卵形鲳鲹卵巢?期和Ш期的cDNA文库,采用Illumina HiSeqTM 2500进行转录组测序,经过滤、质量控制及拼接组装后获得的Unigenes在七大数据库（Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO）中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行功能注释及信号通路富集分析,并采用MISA和GATK3进行SSR鉴定及SNP分析。【结果】卵形鲳鲹卵巢组织转录组测序获得的325156432条Raw reads,经过滤筛选得到317206752条Clean reads,拼接组装后得到59554条Unigenes;69.65%的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等七大数据库中注释成功,其中有24599条Unigenes被注释到GO数据库,15997条Unigenes被注释到KEGG数据库。在卵形鲳鲹卵巢组织的2个发育时期共鉴定获得56115个基因,经差异表达分析后获得17737个差异基因,其中8169个基因在卵巢Ш期上调表达、9568个基因在卵巢Ш期下调表达。GO功能注释分析发现,卵形鲳鲹卵巢差异表达基因主要注释在细胞过程、氮化合物代谢过程、初级代谢过程、核、核部分、离子结合及水解酶活性等条目上;而KEGG信号通路富集分析结果显示,17737个差异表达基因显著富集在318条代谢途径上,其中前20条KEGG信号通路包括2-氧代羧酸代谢、PI3K-Akt信号通路、甲状腺激素信号通路、磷脂酶D信号通路、Fc εRI信号通路和细胞周期等。卵形鲳鲹卵巢转录组（59554条Unigenes）中共存在30133个SSRs和82490个SNPs。【结论】GnRHR、FSHR、FSHβ、CYP11A、SIRT3和PEG3等差异表达基因及PI3K-Akt信号通路和VEGF信号通路等与卵形鲳鲹卵巢的发育密切相关,共同调节卵巢的发育与成熟,在卵巢性成熟过程中发挥重要作用。     

基于RNA-Seq技术的黄连根茎与须根转录组分析

采用RNA-Seq技术，对黄连根茎与须根进行转录组分析，经过trinity软件组装得到130 592条Unigenes，平均长度为700 bp，其中68 976条Unigenes在NR、GO、KEGG、eggNOG、Swiss-Prot和Pfam数据库获得基因注释信息，仅占全部Unigenes的52.82%。39 572个Unigene被注释到GO数据库的生物学进程、细胞组分和分子功能三大类的52个小类中；参与KEGG的5大类代谢通路，34个代谢路径，共搜索到24 999个SSR位点，单核苷酸、两核苷酸和三核苷酸为黄连根部SSR的主要重复基序。黄连根茎与须根的差异表达基因有13 496个，在黄连须根中上调的差异基因有8 375个，下调的有5 121个。两者差异基因GO富集分析表明，跟次生代谢相关进程的前5项为次生代谢进程、次生代谢物生物合成过程、生物碱代谢过程、异喹啉生物碱生物合成过程和异喹啉生物碱代谢过程。在细胞组分分类中，膜结构、细胞壁和外封装结构为差异基因富集的前三；在分子功能分类中，跟氧化还原酶活性相关条目占差异基因富集显著的前四。在KEGG数据库中，共有3 749个差异基因定位到125个代谢途径中，其中，24个差异基因注释到异喹啉生物碱生物合成，在关联的8条KEGG通路中，筛选得到11个关键酶基因；同时，有39个差异基因注释到络氨酸代谢途径中，在关联的21条KEGG通路中，筛选得到10种关键酶基因。黄连根部转录组数据的获得为研究黄连异喹啉类生物碱的合成及黄连根茎与须根代谢差异奠定了基础。