褪黑激素对PDGFA基因在绒山羊皮肤毛囊中表达模式的影响

【目的】探讨血小板衍生生长因子A（platelet-derived growth factor A,PDGFA）在内蒙古阿尔巴斯绒山羊皮肤毛囊上的表达规律和褪黑激素对其的影响。【方法】选择体重和产绒量相当的阿尔巴斯绒山羊成年母羊,分为4组,即从不埋植的对照组和3组不同时间埋植褪黑激素的试验组。试验组在埋植期间每隔一个月按照2 mg?kg-1 BW的剂量在绒山羊耳后皮下埋植褪黑激素。采集12个月（从2009年8月至2010年7月）的受试羊体侧部皮肤样品,用定量逆转录PCR （qRT-PCR）检测皮肤中PDGFA基因的表达量。【结果】各组间PDGFA基因的表达模式不同。对照组在11月时其表达量为全年最低,在12月,随着毛囊周期性生长缓慢进入退行期,其PDGFA基因的表达量也缓慢上升,并在休止期（1-4月）高度表达。而开始埋植褪黑激素后,仅第二年埋植组在翌年1月时表达量显著提高（P<0.05）;在2月与连续两年埋植组一样,随着绒毛开始脱落,PDGFA基因的表达量也有下降趋势;到5月,其表达量显著低于1月（P<0.05）,此时新一轮的绒毛已长出体表。停止埋植褪黑激素后,仅第一年埋植组翌年PDGFA基因的表达模式和毛囊生长周期与对照组基本一致。【结论】绒山羊皮肤组织当中PDGFA基因在毛囊生长期低表达,退行期和休止期高表达。埋植褪黑激素可缩短毛囊生长周期。PDGFA基因是一个控制毛囊从兴盛期向退行期和休止期转化的因子。     

光控和埋置褪黑激素促进绒山羊长绒的皮肤毛囊差异表达基因网络分析

为研究光控增绒和埋置褪黑激素促进绒山羊在非产绒季节(5—8月份)长绒的皮肤毛囊差异表达基因之间的互作关系,并筛选出皮肤毛囊发育相关的重要功能基因,以Agilent绵羊的全基因组表达谱芯片杂交获得的光控增绒羊皮肤差异表达的99个基因及埋置褪黑激素绒山羊皮肤差异表达的83个基因作为研究对象,结合NCBI数据库,利用RStudio程序包和Cytoscape 3.3.0构建差异基因网络并分析。光控增绒试验组差异基因构建了1个主网络和4个小网络;对照组差异基因构建了1个主网络和5个小网络。与表达量结合分析基因网络发现,K25、KAP16-2基因上调促进ISG17、RSRP1、GIMAP1和CYTIP基因的下调,间接促进TNFAIP6基因的下调;CYP17A1、V15和K2-12(K85)基因的上调促进ddit3基因的下调;KAP16-3基因促进CYP17A1基因的同时抑制CYP1A1基因。埋置褪黑激素试验组差异基因构建了2个主网络和7个小网络,对照组差异基因构建了3个主网络和1个小网络。与表达量结合分析基因网络发现,CTNNB1基因的上调促进COL3A1基因的下调;FLT-1基因的上调促进CGA基因的上调和COL1A2基因的下调;PAX6基因的上调促进DOK4基因的上调和LOC101110099基因的下调;P450基因的上调促进BLG基因的下调。光控增绒羊皮肤中筛选出来的基因中角蛋白、角蛋白关联蛋白、CYP17A1、SP-D等23个显著上调基因和CYP1A1、COL6A5、LOC100101238、ddit3等13个显著下调基因对光控增绒羊皮肤毛囊提前进入兴盛期有重要作用;埋置褪黑激素的绒山羊皮肤中筛选出来的基因中CTNNB1、PAX6、DOK4、LOC105602529等13个显著上调基因和酪蛋白家族、胶原蛋白类、LOC101110099等12个显著下调基因之间的信号传导增强,从而促进绒山羊绒毛进入兴盛期。     

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

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

褪黑激素对Bcl-2和Bax基因在绒山羊皮肤中表达的影响

【目的】研究褪黑激素对内蒙古绒山羊皮肤Bcl-2和Bax基因表达的影响,探讨Bcl-2和Bax基因以及褪黑激素与绒山羊绒毛生长和凋亡的关系。【方法】选用8只性别、体重、年龄一致的内蒙古绒山羊母羊,随机分为埋植组和对照组,提取皮肤总RNA,对Bcl-2和Bax基因mRNA表达进行相对定量分析。【结果】①Bcl-2和Bax基因在内蒙古绒山羊皮肤细胞中相对表达量比值（Bcl-2/Bax）的增减与绒山羊一年中绒毛生长、退行以及休止时间特点基本吻合;②褪黑激素对Bcl-2基因在内蒙古绒山羊皮肤中的表达无显著影响（P＞0.05）,但显著上调了Bax基因在各月份的表达（P＜0.01）。【结论】褪黑激素显著下调了Bcl-2和Bax基因在内蒙古绒山羊皮肤细胞中相对表达量的比值（Bcl-2/Bax）,可促进细胞凋亡,加快机体新陈代谢。     

禽致病性大肠杆菌clpV基因缺失对雏鸡气管黏膜 细胞因子-细胞因子受体相互作用通路的影响

为探究clpV基因在禽致病性大肠杆菌（Avian pathogenic Escherichia coli, APEC）感染雏鸡气管黏膜过程中发挥的作用及机制。将APEC野生株（DE17）及其基因缺失株（DE17ΔclpV）感染7日龄雏鸡气管,12、24 h后收集雏鸡气管黏膜细胞进行转录组学测序,筛选野生株和缺失株的差异表达基因,进行GO和KEGG富集分析。结果显示,在感染12 h后,共筛选到108个差异表达基因（36个上调,72个下调）;在感染24 h后,共筛选到59个差异表达基因（34个上调,25个下调）。GO分析表明,差异表达基因主要富集在细胞内信号转导、RNA聚合酶Ⅱ的转录正调控等生物学过程;蛋白质结合、ATP合成、DNA结合等分子功能类;生物膜及膜的组成成分、细胞质等细胞组分功能;KEGG分析表明,差异表达基因主要富集在紧密连接通路、内质网中的蛋白质加工通路、细胞因子与细胞因子受体相互作用通路、内吞通路、Wnt信号通路等。结论为clpV基因缺失后,会使部分差异表达基因的表达量下调,对细胞因子-细胞因子受体相互作用等通路产生影响。     

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

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

亚慢性砷暴露对小鼠心肌细胞凋亡的影响

为研究砷暴露对小鼠心肌组织损伤及细胞凋亡的影响,试验选用40只5周龄C57BL/6小鼠,随机分为4组(对照组饮用去离子水,低、中、高剂量砷暴露组分别饮用含0.15、1.50、15.00 mg/L三氧化二砷的去离子水),90 d后检测各组心脏脏器系数、血清心肌酶谱活力、心肌细胞凋亡率及凋亡通路相关基因的表达。结果表明,不同浓度砷暴露组的心脏脏器系数均有所下降,其中,中、高剂量组与对照组相比差异显著;与对照组相比,中、高剂量组血清肌酸激酶(CK)、谷草转氨酶(AST)、乳酸脱氢酶(LDH)、α-羟基丁酸脱氢酶(HBDH)活力显著升高,低剂量组各项指标变化不明显。TUNEL法检测显示,对照组中未发现有凋亡心肌细胞,各剂量砷暴露组心内膜和心外膜下均可见不同数量的凋亡心肌细胞分布。与对照组相比,心肌组织Bax的表达在中、高剂量组显著上调;Bcl-2的表达在中、高剂量组显著下调;Caspase-3的表达在高剂量组显著上调。亚慢性砷暴露对小鼠心肌造成特异性损伤,且随着染毒剂量的增加损伤程度加重;诱导心肌细胞凋亡是砷暴露致心脏毒性的重要机制。     

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

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

持续感染口蹄疫病毒VP1蛋白对细胞转录组的影响

【目的】通过转录组测序（RNA-Seq）方法分析持续感染口蹄疫病毒结构蛋白VP1诱导的PK-15细胞基因表达谱的变化,探索VP1蛋白对宿主细胞信号通路的影响。为后续开展VP1影响FMDV持续感染相关机制研究提供基础数据和理论支持。【方法】PK-15细胞分别转染pCAGGS-HA（对照载体）、pCAGGS-HA-Pi-VP1质粒,转染36 h后,Western blot检测VP1蛋白表达,同时提取两组细胞总RNA进行转录组测序（每组3个重复,2组6个样本）。将测序所得的序列进行过滤比对,获得差异表达基因后,对差异表达基因进行了生物信息学分析。采用qRT-PCR验证关键差异表达基因的表达。【结果】RNA-Seq共鉴定出3 376个差异表达基因（DEGs）|log2(Fold Change)|>0和padj<0.05,其中上调基因1 744个,下调基因1 632个。许多参与免疫反应和炎症的效应基因差异表达。GO和KEGG富集分析显示,差异基因GO条目绝大部分与核糖体有关,而KEGG富集中涉及了与天然免疫相关的信号通路,如TNF信号通路、Toll受体信号通路。【结论】pCAGGS-...     

逆转座子LINE1在Hela细胞中的表达及转录组分析

为研究逆转座子LINE1(简称为L1)在Hela细胞内的作用机制,对转染L1的Hela细胞和正常细胞进行转录组测序,通过对转录组数据进行拼接组装,对差异表达基因进行了生物信息学分析。结果表明:试验共筛选出391个差异表达基因,其中45.6%的差异表达基因显著性上调,这些基因通过GO分类注释,可分为14类;通过KEGG对差异基因进行通路富集,表明差异基因主要富集在破骨细胞分化、色氨酸代谢、乙型肝炎、病毒致癌作用和可卡因成瘾等信号通路,发现了一些与细胞黏着连接、表皮相关的基因有差异性表达。研究表明,通过转录组测序,获得了丰富的关于L1在Hela细胞中表达的信息,为进一步研究L1在细胞中的作用机制供了思路和理论依据。     

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

[目的]明确相关基因在育成期蛋鸡肝脏脂质代谢调控过程中的作用机理,为改善蛋鸡的生产性能打下基础.[方法]分别选取高能量(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信号通路,且均与蛋鸡肝脏的脂质代谢相关.[结论]育成期对蛋鸡进行高能饮食饲喂,主要是通过上调或下调脂质代谢相关基因表达而影响肝脏脂质代谢,最终实现蛋鸡生产性能调控.     

高温胁迫下不同玉米材料间的基因差异表达分析

【目的】认识高温胁迫条件下不同玉米材料间的基因差异表达,为阐明玉米耐高温胁迫的分子遗传机制提供理论基础.【方法】以耐高温胁迫的中地88和高温敏感的先玉335两个玉米品种为试验材料,在42℃高温条件下取其花粉提取总RNA,再经mRNA纯化、反转录、文库构建及转录组测序等一系列过程,对其基因差异表达作生物信息学分析.【结果】共发现基因24 553个,以先玉335为对照,上调表达基因857个,下调表达基因965个,无显著差异表达基因22 731个,这些基因较均一地分布于玉米的10条染色体.对差异表达基因进行GO富集分析,结果显示:与细胞组分、分子功能及生物过程相关的GO条目数分别为491、1 062和2 320个;在GO注释的差异表达基因中,进一步分析发现直接应答于热激效应的基因有7个,其中2个上调表达,5个下调表达;与热激蛋白结合有关的基因有3个,其中1个上调表达,2个下调表达.对差异表达基因作KEGG富集分析发现有102个小通路,这些小通路可进一步归为6个大通路,包括细胞过程3个、环境信息处理4个、遗传信息处理20个、人类疾病1个、新陈代谢72个和有机系统2个.【结论】转录组测序是在基因组水平认识基因差异表达的有效方法,共发现12个与玉米耐高温胁迫相关的差异表达基因.     

小鼠2细胞胚胎基因表达谱分析

2细胞胚胎是哺乳动物植入前胚胎发育的关键时期。采用Oligo小鼠全基因组芯片对昆白小鼠2细胞胚的前期与后期胚胎总RNA进行芯片杂交,并分析基因表达谱。结果显示,2细胞胚胎表达的基因有5 136个;晚期和早期胚胎差异表达的基因1 248个,其中上调基因847个,下调基因401个;差异表达基因通路分析,多数基因参与细胞周期调控相关通路;基因本体论(GO)分类结果显示,参与细胞组成、细胞加工、催化活性和分子结合的基因占多数,部分基因参与细胞器组成、代谢加工、生物学调控和细胞信号转导等。     

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

采用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信号通路、甲状腺激素信号通路和催产素信号通路与翅型分化及与之紧密联系的生殖率差异相关性较强。     

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

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

基于RNA-Seq数据分析油酸对牛前体脂肪细胞分化及PPAR信号通路的影响

利用100μmol/L油酸对细胞进行诱导分化,通过RNA-Seq技术对分化后48,96 h的细胞进行转录组测序,对差异表达基因进行GO功能注释和KEGG富集分析。结果表明：与分化48 h时相比,分化96 h时差异表达基因数增加3 519个,差异表达基因的分子功能、参与的生物过程和细胞组成随分化时间不同均发生变化,生化代谢途径和信号转导途径也差异较大;由过氧化物酶体增殖物激活受体（Peroxisome proliferatorsactivated receptors,PPARs）信号通路的分析结果可以得出,分化48 h时,PPARγ的表达上调,信号通路成员PLIN2、CPT1β,SREBP1、FABP4上调,SCD表达下调;分化96 h时,PLIN2、SREBP1、CPT1β显著上调,PPARγ、SCD的表达显著下调,FABP4表达无显著变化。说明油酸对延边牛前体脂肪细胞分化过程中的转录组学及信号通路存在一定影响。     

桃嫩枝扦插不定根发育的分子机制

利用高通量测序技术对桃砧木‘GF677’生根发育的4个时期进行转录组测序和数据分析，为阐明桃不定根发育的分子机制奠定理论基础。以扦插后0(对照)、7(激活期)、14(愈伤形成期)、21(不定根形成期)、28 d(不定根伸长期)插穗基部的韧皮部为材料进行RNA-sep测序并研究其生根机理。结果表明，生根过程中不同时期共鉴定出25 656个差异表达基因，其中上调13 166个，下调12 490个。GO功能注释表明，差异表达基因主要富集在代谢过程、细胞过程、细胞、细胞组分、细胞器黏合物和催化活性等代谢途径中。KEGG富集结果显示，差异表达基因主要响应糖酵解、半胱氨酸和蛋氨酸、精氨酸和脯氨酸、核糖体和氨基酸生物合成。应用实时荧光定量qRT-PCR对主要通路的15个DEGs表达模式进行验证，其中13个基因表达水平的变化与转录组基因丰度的变化基本一致，表明转录组数据具有较高的可靠性。总之，糖酵解通路关键基因显著上调，参与激活期的能量补给；蛋氨酸代谢通路诱导基因的表达，参与韧皮部维管组织的分化与形成；精氨酸代谢激活NOS基因的表达，有利于NO的传递，实现根源信号向地上部的运输，进而调控桃的不定根发育...     

基于高通量测序分析肺炎支原体感染姜曲海猪肺组织miRNA表达谱

为探明姜曲海猪感染肺炎支原体后肺组织miRNA(microRNA)表达谱和分子机制,选取50日龄姜曲海猪为实验猪,随机分成感染组和对照组,人工感染肺炎支原体28 d后,采集肺部组织,进行高通量miRNA测序,采用生物信息学软件进行miRNA鉴定和靶基因预测。结果显示:感染组和对照组的肺组织分别筛选到14 265 786条和14 000 588条小RNA纯净序列(clean reads)。与对照组相比,感染组中筛选到73个显著差异表达的miRNAs,其中39个表达上调,34个表达下调,从中随机选取4个miRNAs进行定量PCR(quantitative real time PCR,qRT-PCR)验证,感染组和对照组的表达水平与测序结果基本一致。73个差异表达miRNAs预测到1 685个靶基因和4 220个靶位点,靶基因预测筛选到8个与免疫调控相关的miRNAs。靶基因GO(gene ontology)分析显示,miRNA广泛参与生物过程、细胞组成和分子功能的调控。靶基因KEGG(Kyoto encyclopedia of genes and genomes)分析显示,miRNA参与调控细胞凋亡、ECM受体相互作用、粘附斑通路等信号通路。本研究通过高通量测序获得姜曲海猪感染肺炎支原体后肺组织miRNA表达谱,通过生物信息学分析筛选到与免疫调控相关的miRNAs和信号通路,为进一步阐明姜曲海猪的肺炎支原体感染机制奠定了基础。     

绒山羊胚胎期次级毛囊发生发育相关miRNA的筛选及靶基因验证

为筛选绒山羊次级毛囊发生发育相关microRNA(miRNA)及靶基因验证,本研究基于课题组前期高通量测序技术得到绒山羊胎儿期(45、55、65、75、95、115和135 d)皮肤组织miRNA表达谱的基础上,鉴定各比对组中差异表达的miRNA并根据毛囊形态发生规律筛选次级毛囊发生发育相关的miRNA,运用RNAhybrid、miRanda和TargetScan软件预测筛选到的miRNA靶基因并与次级毛囊发育相关基因取交集,所得基因进行GO和KEGG富集分析,筛选毛囊发生发育重要通路中的关键基因构建miRNA-mRNA调控网络,进而运用双荧光素酶报告验证关键miRNA-mRNA的靶向关系。结果表明:1)筛选到次级毛囊发生发育相关miRNA 153个,共预测到32 978个靶基因;2)取交集后得到1 352个靶向的次级毛囊发生发育相关的基因,且显著富集在29条KEGG通路中,其中毛囊发生发育重要通路TGF-β同样显著富集;3)构建了chi-miR-26b-5p-SMAD1和 chi-miR-495-3p-SMAD1等232对次级毛囊发生发育可能相关的miRNA-mRNA关系对,双荧光素酶报告系统进一步验证了关键靶向关系。本研究为解析miRNA在绒山羊次级毛囊发育过程中的分子机制奠定了前期基础。     

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信号通路及细胞凋亡等能量代谢相关通路上,通过共同协调而对墨瑞鳕的生长发育起调控作用。