藏鸡与白羽肉鸡垂体转录组差异表达研究

试验旨在通过对藏鸡和白羽肉鸡垂体转录组测序和生物信息学分析,筛选出与鸡生长发育相关的差异表达基因（differentially expression genes,DEGs）及信号通路。本研究选用42日龄健康藏鸡和白羽肉鸡各3只,分别采集垂体组织利用Illumina HiSeq 2000平台进行转录组mRNA测序,对筛选到的DEGs筛选DEGs,GO和KEGG数据库功能注释和富集分析,实时荧光定量PCR验证随机挑选的DEGs表达水平。结果显示,藏鸡和白羽肉鸡分别获得126 094 302和125 666 442条clean reads。与白羽肉鸡相比,藏鸡垂体组织中共有DEGs 392个,其中196个为上调基因,196个为下调基因（P<0.05）,其中包括生长激素（GH）基因、生长激素释放激素受体（GHRHR）基因和胰岛素样生长因子2 mRNA结合蛋白1（IGF2BP1）基因。GO和KEGG富集分析发现,DEGs显著富集于细胞黏附分子信号通路和神经活性配体-受体相互作用信号通路等细胞通讯相关的信号通路,GH、GHRHR和IGF2BP1基因也显著富集于神经活性配体-受体相互作用信号通路。实时荧光定量PCR结果显示,转录组测序结果准确可靠。综上研究,本研究初步揭示了影响藏鸡和白羽肉鸡生长发育速度差异的关键候选基因和信号通路,为了解鸡垂体组织调节生长发育的分子机制提供了理论依据。     

基于RNA-Seq技术的塔里木马鹿毛色相关差异表达基因筛选及初步分析

旨在基于RNA-Seq技术对塔里木马鹿毛色相关基因进行筛选及分析。采用Illumina Hi Seq TM2000测序平台对塔里木马鹿和天山马鹿的皮肤组织进行转录组测序,所得序列经质控、组装后比对到NR、Swiss-Prot、COG、KOG、KEGG、GO和Pfam数据库中注释,并对差异表达基因进行筛选、功能注释和富集分析。结果表明,测序获得25 038个有注释信息的Unigenes,比对分析显示,塔里木马鹿与天山马鹿有922个差异表达基因,其中上调表达基因495个,下调表达基因427个;GO功能富集分析结果显示,568个差异表达基因富集到61个GO条目上,分别参与了生物学过程、细胞组分及分子功能;KEGG代谢通路富集分析发现,在差异表达基因中富集最显著的代谢通路是ECM-受体相互作用。利用实时荧光定量PCR（qRT-PCR）方法分析与塔里木马鹿毛色相关的7个候选基因的转录水平变化来验证转录组测序结果的准确性和可靠性,这些基因的表达趋势与转录组测序结果相一致。ECM-受体相互作用、蛋白质消化与吸收、PI3K-Akt信号通路及与黑色素合成相关的酪氨酸等通路可能与塔里木马鹿的毛色有关;候选基因MITF、Ggt1、VDR、PTPRF、CⅡTA、ARPC5L、POMC等可能在塔里木马鹿毛色形成过程中发挥重要作用。本研究结果为今后塔里木马鹿毛色相关基因的分子调控机制方面及挖掘潜在的新基因提供了丰富的试验数据。     

共轭亚油酸影响边鸡胸肌脂类代谢的相关差异基因的筛选

旨在通过转录组测序技术筛选出共轭亚油酸(conjugated linoleic acids,CLA)影响边鸡胸肌脂类代谢的相关基因,为日粮中CLA调控边鸡脂类代谢的作用机制提供理论依据。本试验选取90日龄健康边鸡180只,随机分为5组,每组3个重复,每组日粮中分别添加不同比例的CLA 0%(对照组)、0.5%、1.0%、1.5%和2.0%,预饲期1周,正饲期6周。饲养试验结束后采集胸肌组织进行转录组测序,使用DESeq2软件对测序数据进行差异表达分析,对差异表达基因进行GO功能和KEGG通路富集分析,筛选出与胸肌脂类代谢相关的差异表达基因,利用qPCR对差异表达基因进行验证。通过转录组测序数据分析,共筛选出1 229个差异表达基因,其中594个基因上调,635个基因下调。随机选取9个差异表达基因进行qPCR验证,其相对表达量变化趋势与测序结果一致。GO功能分析发现,差异表达基因主要集中在生物学过程中的细胞过程、单一生物过程、生物调节和代谢过程。KEGG通路富集分析发现,各试验组中差异表达基因显著富集在不同的信号通路中,主要富集在ECM-受体相互作用、黏着斑、吞噬体、细胞凋亡、肌动蛋白细胞骨架调节和细胞因子-细胞因子受体相互作用等通路中。通过GO和KEGG功能注释分析,共筛选出MCAT、APOA1、PTGDS、ALDH3A2、PLTP、FABP3、FOXO1、UCP3和FABP4等18个主要的参与脂类代谢相关的差异表达基因,可能在调控边鸡胸肌脂类代谢中发挥重要的作用。本研究通过转录组测序筛选出CLA影响边鸡胸肌脂类代谢过程中的主要调控基因,发现了日粮中CLA影响差异表达基因作用的主要生物学过程和信号通路,为今后研究CLA影响边鸡胸肌脂类代谢的分子作用机制奠定了基础。     

云南独龙鸡和白来航鸡卵巢组织转录组水平的比较分析

旨在利用RNA-seq技术筛选和分析云南独龙鸡和白来航鸡卵巢组织的差异表达基因,为独龙鸡后续基因功能、品种改良和种质特性等研究奠定基础。本研究采集了3只43周龄健康状况良好的独龙鸡卵巢组织进行转录组测序,并与相似饲养条件和周龄的白来航鸡卵巢组织进行比较。结果表明,与白来航鸡相比,独龙鸡中共有1 273个差异表达基因,其中上调基因522个,下调基因751个。GO和KEGG富集分析显示,差异表达基因主要参与蛋白水解、翻译、细胞外泌体、膜组分、结合和催化等过程,在蛋白合成、ECM-受体相互作用、氧化磷酸化、粘着斑等通路中起到重要作用。对显著性通路分析发现独龙鸡在蛋品质和免疫抗性上具有一定优势,但在产蛋性能上与白来航鸡还存在一定差距。本研究筛选到的RPL22、ATP5I、COX5A、JAKMIP1、CATH2基因主要富集在抗性、产蛋性能和蛋品质相关通路中,可作为重要候选基因进一步深入研究。     

基于转录组测序挖掘广西麻鸡生长发育相关基因

【目的】通过对1和60日龄广西麻鸡腿肌进行转录组测序(RNA-Seq),筛选出广西麻鸡生长发育关键基因。【方法】选择1和60日龄健康母鸡各3只,分别采集腿肌组织,利用Illumina Novaseq^TM 6000平台进行mRNA转录组测序,利用DESeq2软件进行差异表达基因的分析,并对差异表达基因进行GO功能和KEGG通路富集分析,进一步筛选出生长发育相关基因,并用实时荧光定量PCR验证转录组测序结果。【结果】转录组测序结果显示,6个样本分别获得34 296 198~41 647 642条clean reads。与1日龄腿肌相比,60日龄腿肌共获得2 304个差异表达基因,其中998个基因上调,1 306个基因下调。GO功能富集分析发现,共获得富集条目572条,其中生物过程381条,细胞组分70条,分子功能121条。KEGG通路富集分析发现,共获得34个显著富集的信号通路,其中心肌收缩、MAPK信号通路、紧密连接、肌动蛋白细胞骨架信号通路、心肌细胞中的肾上腺素能信号传导等与生长发育相关。通过GO功能和KEGG通路富集分析获得5个生长发育相关基因,分别为肌球蛋白重链10(MYH10)、肌球蛋白链15(MYH15)、成纤维细胞生长因子10(FGF10)、成纤维细胞生长因子16(FGF16)、肌肉生长抑制素(GDF8)基因,其中MYH10、MYH15、FGF10为上调差异表达基因,FGF16和GDF8为下调差异表达基因。所选差异表达基因实时荧光定量PCR结果与转录组测序的基因表达水平一致,说明转录组测序结果可靠。【结论】本试验通过对广西麻鸡1和60日龄2个不同生长阶段进行转录组测序分析,获得MYH10、MYH15、FGF10、FGF16、GDF8共5个与生长发育相关的关键基因,为后续进一步探讨广西麻鸡生长发育的分子调控机制提供参考。     

鹅产蛋性状全基因组选择信号分析

【目的】利用选择信号分析方法在伊犁鹅和霍尔多巴吉鹅中筛选与产蛋性状相关的候选基因。【方法】选取健康状况良好、饲养管理水平一致的2岁伊犁鹅母鹅和霍尔多巴吉鹅母鹅各24只,采集血液样本,提取基因组DNA,利用全基因组重测序技术进行选择信号检测分析,筛选出受到选择的候选区域,针对注释基因进行GO功能和KEGG通路富集分析,进一步筛选出与鹅产蛋性状相关的候选基因。【结果】全基因组重测序的平均测序深度为15.28×,与参考基因组比对率在97.31%以上。通过群体分化指数(Fst)对伊犁鹅和霍尔多巴吉鹅2个群体进行分析,共筛选到1 231个候选区域。与核苷酸多样性(Pi)进行联合分析后,伊犁鹅群体共筛选出10个候选区域,得到5个注释基因;霍尔多巴吉鹅群体中共筛选出353个候选区域,得到263个注释基因。GO功能和KEGG通路富集分析发现,初步筛选出6个可能与鹅产蛋性状相关的候选基因(BMP2、BMP6、MIS、ENO1、LIF、EP300),还发现了IL-18基因可能与禽类的免疫有关。【结论】本研究筛选出6个可能与鹅产蛋性状相关的候选基因,为揭示鹅产蛋性状的分子调控机制提供了参考。     

MSTN基因编辑牛肌肉转录组测序及生物信息学分析

为探究肌生长抑制素（myostatin,MSTN）基因对牛肌肉发育的具体调控机制,本研究选取同一牛场健康鲁西黄牛10头,其中通过转基因技术得到的基因编辑牛MSTN^-/-和同种非转基因野生型牛各5头。分别采集两组牛腿臀肌肉样品,利用IIlumina HiSeq高通量测序技术进行转录组测序分析,通过生物信息学方法比较两组样本间的差异表达基因,并进行GO和KEGG富集分析,最后利用实时荧光定量PCR验证转录组测序数据。结果显示,基因编辑型牛和野生型牛之间共检测到18 071个基因。在log₂|FoldChange|≥ 1.48条件下,筛选出406个差异表达基因,其中347个显著上调,59个显著下调。GO功能富集分析显示,MSTN基因编辑后显著影响915个功能类别（P<0.05）,差异基因主要参与结合、生物系统调节、免疫系统等相关功能。KEGG通路富集分析结果共涉及211个通路,差异基因主要富集在细胞黏附分子、趋化因子信号通路、细胞因子互作等信号通路上,进一步从中筛选出可能参与细胞生长、肌肉发育的差异基因（CD14、KIT、CSF1R、FBP1、DUSP4、ULBP21、PRKCB、SPN、CHAD、SRC）。实时荧光定量PCR检测结果显示,所选差异基因表达水平与转录组表达水平一致,证明测序结果的可靠性。本研究结果表明,MSTN基因发挥作用后可以介导多个下游基因表达,从而影响相关信号通路及生物学过程;同时,所筛选出的差异表达基因可作为进一步研究骨骼肌调控机制的候选靶标。     

基于单细胞转录组测序技术筛选牛体内囊胚性别特异性基因

【目的】探究牛体内雌性和雄性囊胚在mRNA水平上的特异性差异,挖掘雌性和雄性囊胚发育差异相关候选基因。【方法】参考GenBank中牛牙釉质基因（AMEL）序列（AMELX基因,登录号：NM_001014984.1;AMELY基因,登录号：NM_174240.2）设计引物,以胚胎内细胞DNA为模板,采用巢式PCR扩增技术鉴定单个牛体内囊胚的性别。选用确定性别的单个雌性和雄性囊胚为试验材料,采用Smart-Seq2扩增技术构建单个胚胎测序文库,应用Illumina HiSeq Xten高通量测序平台对单个胚胎进行单细胞转录组测序（scRNA-Seq）,并进行了基因差异表达、GO功能和KEGG通路分析。【结果】通过巢式PCR扩增胚胎内细胞中的AMEL基因,确定了胚胎性别;基因表达分析筛选出两组之间差异表达基因（DEGs）6 160个,其中雌性特异性基因675个,雄性特异性基因305个;GO功能注释发现雌性特异性基因显著注释在核苷酸结合、信号转导调控、多细胞生物发育、细胞骨架等条目,雄性特异性基因显著注释在氧化磷酸化、线粒体、线粒体内膜、核糖体等条目;KEGG通路富集分析发现7条与雌性和雄性囊胚发育差异相关的通路,分别为代谢途径、糖酵解、磷酸戊糖途径、细胞衰老、氧化磷酸化、调节干细胞多能性的信号通路和Wnt信号通路;并利用GO和KEGG结果筛选出5个在牛体内囊胚期胚胎发育过程中具有直接或间接作用的基因：FBP1、GADD45G、FHL2、FOSB和WNT2B。【结论】牛体内雌性和雄性囊胚之间存在广泛的转录差异,性别特异性基因FBP1、GADD45G、FHL2、FOSB和WNT2B可能是直接或间接影响胚胎发育的关键基因。     

鸭产蛋量相关遗传变异和功能基因筛选及鉴定

【目的】 通过高通量混池重测序和选择清除分析比较鸭不同产蛋量组间基因组显著差异区域内的SNP和基因差异,以筛选和鉴定出鸭产蛋量相关的遗传变异位点和功能基因,为通过分子遗传育种手段提高鸭产蛋性能提供依据。【方法】 根据金定鸭群体开产后150 d内个体产蛋量情况,选择2种极端表型,分为高产蛋组（CH）和低产蛋组（CL）。基于混池全基因组重测序和选择清除分析技术筛选不同鸭产蛋量组间基因组显著差异区域内的SNP及相关功能基因,通过单个样本PCR扩增子测序对筛选的产蛋量相关SNP进行验证,对筛选的候选基因进行GO功能和KEGG通路富集分析,确定候选基因参与的最主要生化代谢途径和信号转导途径,并分析不同基因型间产蛋量高低差异。【结果】 在低产蛋量组和高产蛋量组分别获得192 071 438和229 836 820条的clean reads,共定位到的SNP差异极显著区间为1 368个,受选择候选基因为214个,而且这些区间和基因主要位于Z号染色体,主要包括KDM4C、LURAP1L、PTCH1、PRUNE2、TRPM3和VPS13AD等基因。验证结果表明重测序结果准确。GO功能分析表明,受选择基因在分子功能、细胞组分和生物过程3个本体中均有富集。KEGG通路富集分析表明,受选择基因主要富集到代谢途径、肌动蛋白骨架调节、真核生物核糖体发生等信号通路。鉴定出候选基因KDM4C上Z-28286537、Z-28286879、Z-28288421、Z-28434122和Z-28436368位点,LURAP1L基因上Z-30802227位点、TRPM3基因上Z-36500134、Z-36503668、Z-36534782、Z-36684262、Z-36710928和Z-36732487位点,VPS13A基因上Z-37498270和Z-37513004位点,PTCH1基因上Z-41510597位点显著影响鸭的产蛋量（P<0.05）。【结论】 鸭Z号染色体上存在多个与鸭产蛋量显著相关的SNPs位点及相关基因,本研究结果为通过分子遗传育种手段提高蛋鸭产蛋性能提供了依据。     

发情周期不同阶段绵羊卵巢lncRNA的鉴定和功能分析

旨在分析不同繁殖周期绵羊卵巢组织中长链非编码RNA (long non-coding RNA,lncRNA)的表达谱,了解lncRNA表达及其调控机理,为绵羊繁育研究提供理论依据。本研究以湖羊为研究对象,选取年龄在1.5～2.5岁的黄体期和卵泡期母羊各3只,通过高通量测序筛选出卵巢组织中的lncRNA,运用生物信息学分析对差异表达的lncRNA进行靶基因预测,通过GO和KEGG富集分析找出与绵羊繁殖相关的通路。结果显示,本研究共获得1 379个差异表达的lncRNAs,其中1 158个表达上调,221个表达下调。GO和KEGG富集分析表明,差异表达的lncRNAs及其靶基因主要参与卵泡发育、排卵周期过程、钙离子信号通路、卵母细胞减数分裂、催产素信号通路、MAPK信号通路、甲状腺激素合成通路、雌激素信号通路。关键的lncRNAs可能通过调控参与这些信号通路和生物学过程的相关基因来调控生殖。其中,LNC_011239、LNC_012847、LNC_003902、LNC_003906、LNC_003907等靶向的MAPK1、ADCY1、ADCY5、PPP3CA和CDC23可能发挥关键的调控作用。qRT-PCR验证证明,随机选取的5个差异lncRNAs定量结果与测序结果基本一致。本研究利用RNA-Seq技术筛选出黄体期和卵泡期卵巢组织中的lncRNAs,并进行差异分析,为揭示绵羊繁殖能力的分子机制提供依据。     

鸡肌内脂肪沉积相关候选基因筛选

本研究对课题组前期在京星黄鸡转录组研究中筛选到的与肌内脂肪（IMF）差异沉积相关的14个候选基因进行验证，检测其在中外两个鸡种群体中与胸肌IMF沉积的关联性。以98日龄慢速型地方鸡京星黄鸡和42日龄快速型白羽科宝肉鸡胸肌组织为素材，通过胸肌甘油三酯（TG）含量区分高低表型组，并检测候选基因在组间的基因表达差异。结果表明，在京星黄鸡胸肌TG高、低组间ATP结合盒亚家族B成员8（ABCB8）、脂联素（ADIPOQ）、第6号染色体开放阅读框65（BEND6）、CD74分子（CD74）、核糖基5-磷酸转移酶（FKTN）、组蛋白乙酰转移酶1（HAT1）、硫酸乙酰肝素-氨基葡萄糖3-磺基转移酶5（HS3ST5）、介体复合物亚基4（MED4）、肿瘤坏死因子超家族成员8（TNFSF8）和TNFAIP3相互作用蛋白1（TNIP1）共10个基因表达差异显著（P<0.05）；在科宝肉鸡TG高、低组间ADIPOQ、BEND6、FKTN、HAT1、HS3ST5、MED4和TNIP1共7个基因表达差异显著（P<0.05）；ADIPOQ、FKTN、HAT1、HS3ST5、MED4和TNIP1共6个基因在两个品种中差异表达趋势一致（P<0.05）。本研究提供了鸡IMF沉积相关新候选基因，为IMF分子调控机理研究和相关分子标记筛选研究奠定了良好的基础。     

简州大耳羊肌内脂肪细胞成脂分化差异表达基因的筛选与鉴定

旨在通过转录组测序技术获得简州大耳羊肌内前体脂肪细胞成脂分化前后的差异表达基因,并经生物信息学分析获得相关信号通路及可能发挥作用的关键功能候选基因。本研究以7日龄的健康简州大耳羊公羊为试验动物（n=3）,采用胶原酶消化法分离获得其肌内前体脂肪细胞;利用Illumina平台对肌内前体脂肪细胞和诱导分化5 d的肌内脂肪细胞cDNA样品（n=3）进行高通量测序;以|log₂fold change|>0和P<0.05为阈值筛选获得差异表达基因,并利用clusterProfiler R包对其进行GO功能和KEGG通路富集;最后利用实时荧光定量PCR（quantitative real-time PCR,qRT-PCR）技术检测功能候选基因在细胞分化前后的表达量变化。结果显示,共获得差异表达基因7 916个,其中4 143个为表达上调基因,主要富集到氧化磷酸化、核糖体合成和三羧酸循环通路等305条通路;3 773个为表达下调基因,且主要富集到甲状腺激素信号通路等303条通路;差异基因GO功能注释中52.8%为生物过程、13.4%是细胞组成和33.8%是分子功能。qRT-PCR结果表明,UCP3、ACACB、ACOT11、ACOX3、APOA1和WISP2在分化前后的山羊肌内脂肪细胞中的表达趋势与RNA-seq结果一致,提示这些基因适合作为下一步研究的功能候选基因。本研究筛选得到简州大耳羊肌内脂肪细胞成脂分化的差异基因,并确认了6个基因可作为功能候选基因。研究结果为阐明肉用山羊肌内脂肪细胞成脂分化的分子调控网络提供基础数据和系统资料。     

猪肌内脂肪沉积相关基因的筛选及其表达特性分析

【目的】 通过分析马身猪和大白猪背最长肌转录组测序数据,筛选肌内脂肪沉积相关基因,并对其表达特性进行分析。【方法】 采集180日龄马身猪和大白猪背最长肌,采用RNA-Seq技术对其进行转录组测序,筛选差异表达基因并进行GO功能和KEGG通路富集分析;采用GeneCards在线查询基因功能,进一步筛选与肌内脂肪沉积相关的差异基因;采用实时荧光定量PCR技术检测差异基因在2个品种猪不同组织以及肌内脂肪细胞成脂分化过程中的表达变化。【结果】 获得2个品种猪背最长肌差异表达基因共280个,其中128个表达显著上调,152个表达显著下调。GO功能富集分析注释到46个条目,其中生物过程24条,分子功能8条,细胞组分14条。KEGG通路分析显著富集到PPAR信号通路、MAPK信号通路和脂肪细胞因子信号等脂肪沉积相关通路。GeneCards功能查询获得TRAF2、DUSP1、ACOT4、NR4A1、SLC27A6、PLIN5共6个与脂肪沉积相关的基因。实时荧光定量PCR分析表明,马身猪和大白猪背最长肌中NR4A1、DUSP1、PLIN5基因表达量差异显著或极显著(P<0.05;P<0.01),TRAF2、ACOT4和SLC27A6基因表达量差异不显著(P>0.05),但表达变化趋势与测序结果一致;马身猪腹部皮下脂肪组织中NR4A1和PLIN5基因表达量极显著高于背部皮下脂肪组织(P<0.01),而大白猪背部皮下脂肪组织中NR4A1和DUSP1基因表达量极显著高于腹部皮下脂肪组织(P<0.01),PLIN5基因表达量显著高于腹部皮下脂肪组织(P<0.05)。在背部皮下脂肪组织中,NR4A1、DUSP1和PLIN5基因在大白猪中的表达量极显著高于马身猪(P<0.01);在腹部皮下脂肪组织中,NR4A1和DUSP1基因在大白猪中的表达量显著高于马身猪(P<0.05),PLIN5基因表达量在2个品种间无显著差异(P>0.05)。体外培养的猪肌内脂肪细胞分析表明,随着成脂天数的增加,NR4A1基因表达量呈下降趋势,DUSP1和PLIN5基因表达量均呈上升趋势。【结论】 本研究获得了NR4A1、DUSP1、PLIN5 3个与猪肌内脂肪沉积相关的候选基因,可为后续进一步探讨猪肌内脂肪沉积调控机制提供一定的参考。     

大型迪庆藏猪不同生长阶段肌内脂肪沉积差异表达基因及其调控通路分析

【目的】筛选大型迪庆藏猪不同生长阶段肌内脂肪含量差异的关键基因并分析其调控途径。【方法】选择胎次相同、出生日期及体重相近的大型迪庆藏猪36头,随机分为3组,在相同条件下进行育肥试验,分别在体重达40、80及120 kg左右时屠宰,每组采集3头猪的背最长肌,采用RNA-Seq技术进行转录组测序,测序数据进行拼接、比对和注释,筛选与脂肪沉积相关的差异显著基因进行功能富集、STEM分析,构建基因互作网络,并进行实时荧光定量PCR验证。【结果】40 kg vs 80 kg、80 kg vs 120 kg与40 kg vs 120 kg阶段分别检测到730、981及735个基因差异显著表达;STEM分析共有4个差异显著模块,模块11、14的基因集先显著上调后轻微下调;模块9、10的基因集先显著上调后显著下调。差异基因互作网络显示,40 kg vs 80 kg阶段,EGR1、EGR2、PRKAG2、NOR-1和ATF3基因位于网络核心,EGR1和EGR2基因表达下调,PRKAG2、NOR-1和ATF3基因表达上调;80 kg vs 120 kg阶段,FOXO1、PDK4、PPARD、PPARGC-1、LIPE、ATF3和STAT1基因位于网络核心,FOXO1、PPARD、PPARGC-1基因表达下调,PPARG基因通过级联调控引起STAT1基因表达下调,LIPE和ATF3基因表达上调;40 kg vs 120 kg阶段,ATF3、NOR-1、EGR1、EGR2和STAT1基因位于网络核心,EGR1、EGR2和STAT1基因表达下调,ATF3和NOR-1基因表达上调。ATF3、FOXO1等10个基因的实时荧光定量PCR验证结果与转录组测序结果一致。【结论】EGR1、FOXO1等基因作为核心基因参与大型迪庆藏猪肌内脂肪调控,不同生长阶段参与调控的核心基因并不完全相同,结果可丰富中国地方猪肌内脂肪调控基础数据,为大型迪庆藏猪肌内脂肪含量的遗传改良提供参考。     

多浪羊与小尾寒羊皮下脂肪组织转录组分析

旨在探索多浪羊（D组）与小尾寒羊（X组）皮下脂肪组织中的特异性表达基因,并研究其潜在的作用,为理解绵羊脂肪组织发育规律以及对脂代谢相关疾病的预防和治疗研究提供依据。本研究选取脂肪沉积能力存在差异、健康无病、体况良好、种内个体体重相近（约50 kg）的雌性成年多浪羊和小尾寒羊为试验材料,分为D组（试验组）和X组（对照组）,每组3个重复,采集位于背最长肌的皮下脂肪组织,应用RNA-Seq技术和生物信息学方法进行转录组测序并对结果进行分析。以|Fold change|≥2,P adjust≤0.05为标准筛选差异表达基因,通过对差异表达基因进行GO功能注释和KEGG通路富集分析,得到与脂肪沉积和脂代谢有关的差异基因。为了验证测序数据的可靠性,本研究随机选取了6个差异表达基因进行qRT-PCR验证。结果显示,在6个样本中共检测到38 672个已知的mRNAs,新的mRNAs为1 606个,在两组中共有839个差异表达基因,其中有320个差异基因在多浪羊组中上调表达,有519个差异基因在多浪羊组中下调表达。通过GO功能注释分析发现,差异表达基因主要参与脂质分解代谢过程、脂质生物合成过程、脂质分解代谢负调控过程、MAPK级联反应调控、对甘油三酯的反应等生物学过程。KEGG通路富集结果显示,差异表达基因显著富集到了PI3K-Akt、MAPK、胰岛素以及PPAR等信号通路中。qRT-PCR结果与测序结果一致,表明测序结果可靠。通过对多浪羊和小尾寒羊皮下脂肪组织进行转录组测序以及生物信息学分析,筛选到与脂肪沉积和脂代谢相关的差异表达基因,这些基因主要参与脂质生物合成、脂质代谢等过程,其中COL1A1、AKT2、SCD、LPL、PCK1与PPP2R5A可能在多浪羊与小尾寒羊的皮下脂肪组织的沉积与代谢中发挥重要作用。