青蛤(Cyclina sinensis)IRAK-4基因的克隆及其组织间的表达分析

利用转录组测序获得了青蛤(Cyclina sinensis)白介素-1受体相关激酶-4(Interleukin-1 receptor-associated kinase 4,IRAK-4)基因的序列信息,利用巢式及荧光定量PCR克隆并分析了Cs IRAK-4基因在青蛤各组织中的表达情况,进一步在鳗弧菌的胁迫下检测了Cs IRAK-4基因在青蛤血淋巴中的时序表达情况。结果显示,Cs IRAK-4基因的序列全长共2 687 bp,其开放阅读框长1 926 bp,共编码641个氨基酸,Cs IRAK-4基因在青蛤的血淋巴、肝脏、外套膜、闭壳肌、性腺和鳃中均表达,其中在血淋巴中表达量最高,肝脏中表达量最低。青蛤在鳗弧菌胁迫下该基因的表达量在3 h即达到最大值,与对照组差异极显著(P0.01),证明该基因所表达的产物是青蛤重要的免疫信号通路蛋白。     

青蛤MAPK/p38基因的克隆及Poly∶Ic胁迫下的表达分析

[目的]克隆青蛤MAPK/p38基因,并分析其在Poly:Ic胁迫下的表达情况。[方法]利用转录组测序获得青蛤MAPK/p38基因的序列信息,采用生物信息学软件分析该基因在近缘生物间的进化关系;利用巢式及荧光定量PCR技术克隆青蛤MAPK/p38基因的结构域序列,并在Poly:Ic胁迫后立即分析该基因在青蛤各组织中的表达情况以及血淋巴中的时序性表达情况。[结果]克隆得到青蛤MAPK/p38基因的结构域序列,分析发现该基因在物种进化中很保守,其在青蛤的血淋巴、肝脏、外套膜、闭壳肌和鳃中均有表达,其中在血淋巴中表达量最高,腮中表达量最低。Poly:Ic胁迫3 h后该基因的表达量达到最大值,与对照组差异极显著(P0.01)。[结论]MAPK/p38基因所表达的产物是青蛤重要的免疫信号通路蛋白。     

3种白鲫杂交子代的转录组学分析

[目的]研究3种白鲫杂交子代转录组学特征,为揭示鲫鲤杂交优势的分子机理提供理论依据,同时为在生产上培育出生长速度快、肉质好、适应能力强的杂交品种提供技术参考.[方法]以白鲫(♀)×黑龙江野鲤(♂)杂交子代(简称HB)、白鲫(♀)×散鳞镜鲤(♂)杂交子代(简称SB)和白鲫(♀)×兴国红鲤(♂)杂交子代(简称XB)为研究对象,利用RNA-seq高通量测序技术构建3种白鲫杂交子代转录组文库,以HiSeq PE150进行测序分析,原始序列经Trinity组装后进行功能注释(E-value<1e-5);以DESeq2 R鉴定差异表达基因,利用GOseq R和KOBAS分别对差异表达的基因进行GO和KEGG富集分析;并采用MicroSAtellite对转录本中的SSR位点进行挖掘.[结果]共组装得225858条unigenes,平均长度为668 bp,N50为938 bp,有171461条unigenes可注释到蛋白质数据库(Nr)、非冗余核苷酸数据库(Nt)、蛋白质序列数据库(SwissPort)、基因本体论(GO)、直系同源基因簇(COG/KOG)和京都基因与基因组百科全书(KEGG)数据库中,注释比例为75.92%.其中,52630条unigenes注释到NR数据库,43659条unigenes注释到SwissPort数据库,35756条unigenes注释到COG/KOG数据库,包括生化代谢、信号转导机制、防御系统和细胞结构等.差异表达基因KEGG分析结果显示,较多的差异表达基因注释到内吞作用、Jak-STAT信号通路、溶酶体、吞噬体和Wnt信号通路等免疫相关及与生长发育相关的MAPK信号通路、Hippo信号通路和背腹轴形成等通路中.此外,从获得的转录组序列中共鉴定出20272个SSR位点,大多数为二核苷酸重复基元(占62.15%).[结论]不同白鲫杂交子代间存在较多的差异表达基因,从中获得参与抗氧化、免疫和生长发育相关的通路和基因序列,且挖掘出20272个SSR位点,有助于选择性育种、分子标记开发及开展遗传多样性、遗传图谱构建和QTL定位等研究.     

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

【目的】研究半番鸭与番鸭精巢组织转录组差异表达基因,为进一步阐明半番鸭不育的遗传机制奠定理论基础。【方法】利用转录组测序方法对半番鸭和番鸭的精巢组织进行研究,筛选其差异表达基因,并对其功能进行注释和荧光定量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信号通路在鸭的生殖活动中发挥了重要作用,为进一步探索半番鸭生殖系统的分化机理提供可靠的参考依据。     

青蛤(Cyclina sinensis)ERK基因的克隆及其在Poly I:C刺激下的表达分析

[目的]研究青蛤(Cyclina sinensis) ERK基因的克隆及在Poly I:C刺激下的表达情况。[方法]在已建立的青蛤转录组文库中筛选得到青蛤细胞外调节蛋白激酶(extracellular regulated protein kinases)基因类似序列。运用生物信息学软件在线对该基因进行结构分析。采用PCR技术克隆基因,并使用实时荧光定量PCR技术克隆得到Cs ERK基因在青蛤5个不同组织中的表达情况及在干扰素诱导剂PolyI:C的刺激下ERK基因在青蛤血淋巴中的时序性表达情况。[结果]ERK基因序列全长为2 407 bp,开放阅读框长1 338 bp,共编码445个氨基酸。ERK基因在青蛤的血淋巴、外套膜、闭壳肌、肝脏和鳃5个组织中均表达,在血淋巴中表达量最高,闭壳肌中表达量最低。青蛤ERK基因在Poly I:C胁迫下表达量在6 h时达到最大值,与对照组相比差异极显著(P<0.01)。[结论]ERK基因所指导合成的蛋白是青蛤重要的免疫信号通路蛋白。     

基于转录组测序对紫花苜蓿细胞质雄性不育系相关代谢通路的鉴定

【目的】研究紫花苜蓿(Medicago sative L.)细胞质雄性不育(CMS)的机理,为紫花苜蓿利用不育系的育种工作提供理论基础。【方法】本研究以紫花苜蓿细胞质雄性不育系MSGN1A及其保持系MSGN1B为材料,进行转录组测序(RNA-Seq)以及基因功能注释,筛选与不育系机理相关的差异表达基因及代谢通路。【结果】Illumina测序共得到紫花苜蓿95 679条功能基因,经过差异表达分析筛选出187个显著差异表达的基因(DEGs),其中包括18个上调基因,169个下调基因。采用GO、KEGG、COG注释库分别对187个DEGs进行功能注释,其主要涉及催化活性、代谢活动、信号传导机制、合成、膜功能、碳水化合物运输和代谢、能量代谢等相关通路。【结论】筛选出与紫花苜蓿细胞质雄性不育性相关的代谢通路,初步探讨了紫花苜蓿细胞质雄性不育系的分子机理。     

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

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

罗氏沼虾转录组免疫相关SNP的挖掘与分析

为了挖掘罗氏沼虾免疫相关单核苷酸多态性(SNP),采用第二代高通量测序技术,对罗氏沼虾的肝胰腺进行转录组深度测序分析。结果发现SNP位点37 589个,其中转换25 188个,颠换12 401个。再对SNP所在的基因进行GO(gene ontology)分类和KOG注释,分别筛选到2 451个有GO注释和2 458个有KOG注释的SNP,并重点对4个重要免疫通路的145个免疫相关基因,筛选出129个SNP位点。相关研究结果将为进一步系统挖掘免疫相关SNP的分子遗传标记奠定基础,为罗氏沼虾疾病防控和选育抗病新品种提供科学参考。     

2种不同生长规格墨瑞鳕肌肉组织转录组测序分析

【目的】掌握不同生长规格墨瑞鳕个体间差异表达基因的表达特点,为其功能相关基因深度挖掘及分子遗传育种提供科学依据。【方法】挑选同一养殖条件下极大个体和极小个体的墨瑞鳕,构建肌肉组织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且|log₂ 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信号通路及细胞凋亡等能量代谢相关通路上,通过共同协调而对墨瑞鳕的生长发育起调控作用。     

FST基因的过表达对猪成纤维细胞内基因及生物学通路的影响

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

西方蜜蜂工蜂不同虫态发育的转录组学分析

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

萝卜蚜有翅型和无翅型转录组差异分析

采用Illumina测序平台获得了萝卜蚜转录组数据,有翅蚜和无翅蚜混合样本共得到107 992 806条Cleaning Reads,组装得到69 588条Unigene序列,平均长度为716 bp。有翅蚜测序结果共得到44 620 184条Cleaning Reads;无翅萝卜蚜测序结果共得到47 113 716条Cleaning Reads。2组样品差异表达基因的分析表明共有显著性差异表达基因104个,其中表达上调基因有74个,表达下调基因有30个。对差异表达基因进行GO注释及KEGG富集分析显示,Hippo信号通路、甲状腺激素信号通路和催产素信号通路与翅型分化及与之紧密联系的生殖率差异相关性较强。     

松乳菇子实体两个发育时期的转录组分析

以松乳菇子实体为研究材料,采用de novo测序对松乳菇子实体进行测序与生物信息学分析,探究松乳菇子实体生长发育过程中幼菇期(LDG-1)和成熟期(LDG-2)基因表达的差异,并探寻影响松乳菇生长发育相关的主要基因和代谢通路。测序结果显示,共获得7.44G(LDG-1)和7.21G(LDG-2)的分析数据(clean reads)。KOG功能注释结果显示,松乳菇在幼菇期(LDG-1)和成熟期(LDG-2)具有全面而复杂的基因功能类别。GO富集结果提示,松乳菇在不同生长阶段(LDG-1和LDG-2)的生物代谢功能存在一定的差异。基因表达水平与聚类分析显示,ATP酶、多功能伴侣蛋白等家族基因为LDG-1和LDG-2主要高表达基因,在能量代谢与参与调控细胞周期中起重要作用。分析差异表达基因发现,相较于LDG-1期,LDG-2期有16 789个差异表达基因,其中7 635个基因上调,9 154个基因下调。差异表达基因KEGG富集通路分析结果显示,氧化磷酸化途径是松乳菇幼菇期生长发育过程中能量代谢的关键途径。进一步分析显示,MAPK信号通路是影响松乳菇子实体发育的关键信号通路,在促进松乳菇子实体进行有性生殖和调控菌丝体的极性生长中起重要作用。     

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

【目的】鉴定筛选出与卵形鲳鲹卵巢发育相关的候选基因及信号通路,为揭示其卵巢性成熟过程的分子机制打下基础。【方法】挑选卵巢发育处于?期和Ш期的雌性卵形鲳鲹,分别构建卵形鲳鲹卵巢?期和Ш期的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信号通路等与卵形鲳鲹卵巢的发育密切相关,共同调节卵巢的发育与成熟,在卵巢性成熟过程中发挥重要作用。     

低温胁迫下牛皮杜鹃MAPK级联参与ABA信号转导的基因表达分析

为了探究牛皮杜鹃在低温胁迫下MAPK级联途径参与ABA信号转导的分子机制,通过转录组测序的方法对4℃低温处理组和25℃正常对照组牛皮杜鹃进行研究。结果表明,转录组测序共得到6.40 Gb clean data,低温组筛选出12 261个差异表达基因,其中上调基因和下调基因的数量分别为6 811个和5 450个。对关键基因进行KEGG注释发现,共有228个差异表达基因富集到MAPK信号通路-植物通路,其中ABA信号转导通路中大部分差异表达基因发生上调,说明该信号通路在低温胁迫下被激活。推测牛皮杜鹃在低温胁迫下可激活MAPK级联途径中相关基因的表达来参与ABA信号转导过程,进而应对低温环境所带来的不利影响。     

基于RNA-seq技术分析意大利蜜蜂幼虫肠道响应球囊菌胁迫的高表达基因

对意蜂幼虫肠道组织基因进行深度RNA测序,并对高表达基因进行GO注释和KEGG分析,旨在揭示意蜂幼虫免疫抗病机制。分别取3、4、5、6日龄意蜂幼虫中肠组织抽提RNA,合成cDNA。得到的总读段数都在26 715 638条以上,有效数据达到97%以上,且测序饱和度良好,说明测序所得数据真实可信。将测序有效数据和意蜂公布DNA数据库进行比对,发现意蜂幼虫肠道中有13 789个基因有不同程度表达。选取各样本FPKM排序最前且表达水平最高的100个基因进行GO注释,这里主要涉及细胞组分、分子功能和生物学进程功能分析。利用KEGG分析显示这100个高表达基因分布在能量代谢、翻译、组装和降解以及信号转导等多条通路上。初步分析相应共有基因和特有基因在相关通路上的作用,并对意蜂幼虫肠道受球囊菌胁迫后的相关免疫机制进行一定分析。     

热应激奶牛血清miRNA表达谱及miR-181a靶基因分析

通过分析热应激奶牛和正常奶牛血清差异表达miRNA,及对重要的候选miRNA进行靶基因预测及相关生物信息学分析,探索miRNA在荷斯坦奶牛发生热应激过程中的作用及其调控机制。利用miRNA高通量测序技术对奶牛血清中miRNA表达谱进行检测和分析以筛选差异表达miRNA;用RT-qPCR方法验证4个在热应激奶牛血清与正常奶牛血清中显著差异表达的候选miRNAs;采用Targetscan生物信息学计算方法预测目标miR-181a(对应激以及免疫反应可能起重要调节作用)的靶基因,并对靶基因进行GO注释描述、富集分析及KEGG富集。结果表明:miRNA高通量测序分析发现共有52个差异表达极显著miRNA(P0.01)。候选的4个miRNAs(miR-19a、miR-19b、miR-181a、miR-1246)的RT-qPCR检测结果与测序分析结果具有很好的一致性。GO和KEGG分析结果显示,miR-181a的靶基因与应激反应以及免疫功能密切相关,并且靶基因在B细胞和T细胞受体信号通路中显著富集。结论:在热应激奶牛与正常奶牛血清中存在差异表达的miRNAs,某些重要的miRNA(如miR-181a)可能以独特的通路(如B细胞和T细胞受体信号通路)来调节奶牛热应激的发生过程。     

乌鳢感染鰤诺卡氏菌的转录组反应（英文）

利用转录组技术对鰤诺卡氏菌（Nocardia seriolae）感染的乌鳢（Channa argus）进行分析与比较,以揭示可能的免疫机制,尤其在Toll样受体和Nod样受体的免疫反应．从乌鳢的头肾提取总RNA,采用lllumina HiSeqTM2500进行测序和de novo转录组分析．乌鳢感染鰤诺卡氏菌后,对照组和试验组的clean reads分别是33556284（93.79%）和34202766（93.52%）．利用Trinity软件对unigenes进行注释,从133999条unigenes里,一共106319条unigenes得到注释,占79.34%;NR(NCBI non-redundant protein database)注释54886(40.97%)条unigenes;39795(29.69%)的Swiss-Prot注释;5885(4.39%)条unigenes被KEGG(Kyoto encyclopedia of genes and genomes)注释到335条通路;5753条(4.29%)被注释到GO（Gene ontology）里.3912条unigenes表达差异,其中1552条表达上调(39.67%)和2360条表达下调(60.33％).Toll样受体和Nod样受体的免疫信号通路分别有28和13条差别表达基因发现有调制,12和3个上调,11和6个下调,5和4个unigenes没有显著的变化.经过转录组测序分析显示,乌鳢感染鰤诺卡氏菌后的转录复杂程度高,导致感染个体重诱导表达的基因产品的验证发现,尤其在非特异性免疫系统方面．     

冷冻胁迫转金柑MLP2–1基因拟南芥的转录组测序和代谢通路

利用农杆菌介导的花序浸润法,将金柑(Fortunella margarita(Lour.)Swingle)MLP2–1基因转入野生型拟南芥(Col–1),对转MLP2–1基因拟南芥在冷冻胁迫下目的基因的转录组数据进行分析。结果显示:转MLP2–1基因拟南芥的抗冻能力强于野生型拟南芥;转录组测序分析表明,转MLP2–1拟南芥与野生型拟南芥植株差异表达基因为1 422个,其中551个表达上调,871个表达下调;差异表达基因数目在GO数据库注释到1 208个,在KEGG数据库中注释到474个,在COG数据库中注释到603个;KEGG代谢通路富集结果表明,差异表达基因主要集中在植物激素与信号转导代谢通路。