miR-184及miR-205在绵羊不同生长时期毛囊及其他组织中的表达研究

为了研究miR-184及miR-205在绵羊毛囊发育过程中的作用,本研究以绵羊为研究材料,检测了miR-184及miR-205在毛囊3个发育时期及不同组织中的表达情况,进行了基因组定位及前体预测,并预测了可能的靶基因及调控的基因通路.结果发现,miR-184的表达量从毛囊生长期到衰退期呈逐步上升的趋势,而miR-205则表现为先上升后下降,其在衰退期表达量最高.2个microRNA在绵羊的多个组织中均有表达.两者的候选靶基因涉及15个基因通路,可能通过Wnt通路调控毛囊的生长周期.     

乾华肉用美利奴羊次级毛囊小RNA文库构建及miR-1-3p靶基因预测

本试验旨在构建乾华肉用美利奴羊次级毛囊不同发育时期的小RNA文库和预测筛选miRNA的候选靶基因。通过高通量测序技术构建乾华肉用美利奴羊次级毛囊小RNA文库|采用荧光定量PCR对高通量测序结果进行验证|利用miRDeep2软件对miRNA（长度为20~24 nt的小RNA）进行鉴定|利用Targetscan v7.0、RNA22 v2.0和MiRanda三个在线靶基因预测软件对miR-1-3p的靶基因进行预测分析。成功获得乾华肉用美利奴羊次级毛囊生长期、退行期、休止期小RNA纯净reads分别为10279515、11712215和10926258。分别基于荧光定量与高通量测序的5个miRNA的相对表达量结果趋势基本一致。从次级毛囊生长期、退行期和休止期分别鉴定获得miRNA 1250、1304个和1266个。确定成纤维细胞生长因子14（FGF14）和类胰岛素生长因子Ⅰ 受体（IGF1R）为miR-1-3p的候选靶基因。研究结果为后续靶基因验证及确定调控乾华肉用美利奴羊毛囊周期性发育的候选基因奠定基础。 [关键词] 乾华肉用美利奴羊|小RNA文库|miR-1-3p|靶基因预测     

miRNA与毛囊相关信号通路关系研究进展

毛囊是哺乳动物皮肤中复杂的微小器官，具有明显的生长周期，包括生长期、退行期和休止期。许多细胞因子和信号通路相互作用共同调节毛囊的发育。微小RNA(micro RNA,miRNA)是由22个核苷酸组成的非编码RNA,通过降解靶基因mRNA或抑制靶基因mRNA翻译来负向调节基因的表达。miRNA通过调控毛囊关键的信号通路，参与细胞增殖、凋亡等多种生理过程，从而促进毛囊生长。笔者综述了miRNA调控毛囊发育相关的Wnt、Notch、骨形态发生蛋白(bone morphogenetic protein, BMP)、丝裂原活化蛋白激酶(mitogen-activated protein kinase, MAPK)信号通路的最新研究进展，旨在为后续探索动物毛囊发育相关的miRNA调控机制提供理论基础。     

绒山羊皮肤毛囊miR-1298-5P靶基因预测及表达载体构建

探讨辽宁绒山羊皮肤毛囊不同发育时期miR-1298-5P的调控作用及其与靶基因的潜在关系,为miR-1298-5P与其靶基因对皮肤毛囊发育作用提供理论依据。以辽宁绒山羊皮肤毛囊组织为材料,通过miRNA的分离、表达、靶基因预测筛选表达及靶位点3'UTR表达载体构建,采用RT-PCR进行表达检测,用生物信息学方法预测miR-1298-5P的靶基因,并对其靶基因Fg F2的靶位点进行克隆测序、缺失突变。miR-1298-5P及其靶基因在皮肤毛囊不同发育时期均呈现差异性表达,二者之间呈现一种负调控趋势,其可能对皮肤毛囊周期性发育起到一定的调控作用。靶基因预测发现Fg F2的3'UTR存在mi R-1298-5P种子区结合位点,并成功构建了Fg F2 3'UTR表达载体,该载体的成功构建为后期过表达转染体系的建立和功能基因的验证奠定了基础。     

WNT信号通路关键基因WNT5a和WNT10b在猪胚胎期毛囊的表达规律

毛发是人和动物的重要生理特征，研究表明，猪可以作为人类毛发研究的良好模式动物。本实验通过采集猪胚胎期毛囊发育重要时期的皮肤样本，利用实时荧光定量PCR技术，检测WNT(Wingless-Type MMTVIntegrationSiteFamily)信号通路中关键基因WNT5a和WNT10b在猪胚胎期毛囊发育关键时间节点的mRNA表达量，探究WNT信号通路在猪胚胎期毛囊发生发育中的作用。结果显示，WNT5a基因在猪毛囊发育诱导期高表达，其中在胚胎期39 d表达量极显著升高，在胚胎期41 d(基板形成时期)表达量极显著降低；在器官形成期和细胞分化期表达量均低于诱导期。WNT10b基因在诱导期表达量较低，在胚胎期41 d表达量极显著升高，在器官形成期和细胞分化期的临界点胚胎期85 d表达量极显著升高，之后又极显著降低。本研究结论表明，WNT5a,WNT10b参与了猪毛囊早期发育的调控过程，其中WNT5a、WNT10b均促进毛囊基板形成，WNT10b促进毛囊细胞分化。     

细毛羊皮肤毛囊不同发育时期miR-1298-5p靶基因预测及验证

为了探讨miR-1298-5p及其靶基因与细毛羊毛囊生长发育的关系,以细毛羊不同发育时期皮肤毛囊组织为材料,通过生物信息学方法预测miR-1298-5p的靶基因,采用RT-PCR对miR-1298-5p与其潜在靶基因FGF_2进行核酸水平相对定量检测,利用Western Blot对潜在靶基因FGF_2进行蛋白相对定量检测。生物信息学预测结果表明,FGF_2的3'UTR存在miR-1298-5p种子区结合位点;miR-1298-5p在皮肤毛囊不同发育时期呈现差异性表达;miR-1298-5p与其潜在靶基因之间呈现一种负调控趋势,说明miR-1298-5p可能通过负调控FGF_2从而调控皮肤毛囊生长发育。表明结合靶基因预测结果、miR-1298-5p在核酸水平表达规律及FGF_2在核酸水平表达规律和蛋白表达水平表达规律初步确定FGF_2为miR-1298-5p的靶基因。     

绵羊miR-150靶基因预测及生物信息学分析

为探究miR-150在绵羊卵巢中的调控机制,本研究利用生物信息学方法对miR-150进行靶基因预测与分析。通过miRBase数据库获取miR-150成熟序列,进行序列比对和保守性分析;利用Ensemble数据库搜索miR-150前体序列,并利用PROMO在线软件对其上游2 000 bp的启动子区进行转录因子结合位点预测;通过TargetScan和miRwalk在线软件预测并取靶基因的交集,采用DAVID和KOBAS在线软件对预测靶基因进行GO功能和KEGG通路富集分析;利用Networkanalyst在线软件对靶基因进行功能性蛋白的互作分析,绘制miR-150靶基因参与的调控网络;利用YM500V2在线软件分析miR-150靶基因在不同组织和疾病中的差异表达水平。结果显示,miR-150成熟序列在不同物种间具有高度保守性;miR-150启动子区存在p53、ID1、FOXO4等转录因子结合位点。GO功能分析结果显示,主要富集在RNA聚合酶Ⅱ启动子转录的正调控等生物过程,锌离子结合、DNA结合、转录DNA模板等分子功能,以及核、核质、胞质、胞质核周区域等细胞组分。KEGG通路富集结果显示,主要富集在癌症通路、癌症中蛋白聚糖、人乳头瘤病毒感染、乳腺癌、调节干细胞多能性等信号通路。miR-150靶基因互作分析发现,靶基因PIK3R1、PIK3CB、CTNNB1与其他蛋白存在较多靶向关系,参与了雌激素信号通路、细胞衰老、磷脂酰肌醇信号系统、TGF-β等通路。根据YM500V2在线软件分析miR-150在不同组织和疾病中的表达发现,在血液中表达水平最高,在卵巢、胰腺、淋巴和皮肤中表达水平相对较高,在脑和肾上腺皮质中表达水平很低或不表达。通过对绵羊miR-150靶基因的生物信息学分析,为其功能和调控机制研究提供参考依据,为绵羊miR-150在卵巢发育中的调控机制奠定基础。     

辽宁绒山羊皮肤毛囊miR-1298-5p靶基因预测及验证

以辽宁绒山羊不同发育时期皮肤毛囊组织为材料,用生物信息学方法预测miR-1298-5p的靶基因,采用RTPCR对miR-1298-5p与其潜在靶基因TGF-βR1进行核酸水平表达检测,利用Weston blot对潜在靶基因TGF-βR1进行蛋白水平表达检测。生物信息学预测结果表明TGF-βR1的3′UTR存在miR-1298-5p种子区结合位点;miR-1298-5p在皮肤毛囊不同发育时期呈现差异性表达;TGF-βR1在皮肤毛囊不同发育时期无论在核酸水平还是蛋白水平均呈现差异性表达;miR-1298-5p与其潜在靶基因之间呈现一种负调控趋势,说明miR-1298-5p可能通过负调控TGF-βR1,从而对毛囊周期性发育起到调节作用。结合靶基因预测结果、miR-1298-5p在核酸水平表达规律及TGFβR1在核酸水平表达规律和蛋白水平表达规律初步确定TGF-βR1为miR-1298-5p的靶基因。本实验旨在探讨辽宁绒山羊皮肤毛囊不同发育时期miR-1298-5p与其靶基因的调控作用及其与靶基因的潜在关系,为miR-1298-5p与其靶基因对皮肤毛囊发育作用提供理论研究数据。     

成年辽宁绒山羊皮肤毛囊周期性发育不同时期FGF_5 mRNA的表达分析

成纤维细胞生长因子(FGF_5)基因在次级毛囊的周期性发育过程中具有重要的调控作用。利用RT-PCR技术检测了FGF_5基因mRNA在成年辽宁绒山羊皮肤毛囊不同时期(兴盛前期、兴盛期、退行期和休止期)的表达规律。结果表明:FGF_5基因mRNA在不同时期的辽宁绒山羊皮肤组织中均有表达,表达量由高到低顺序依次为休止期、兴盛期、退行期和兴盛前期;在休止期表达量最高,兴盛前期表达量最低;其中休止期表达量与兴盛前期和退行期间差异均显著(P0.05),兴盛前期与其余3个时期间也差异显著(P0.05),但休止期与兴盛期、兴盛期与退行期间差异均不显著(P0.05)。因此,FGF_5基因表达可能控制着绒山羊的毛囊从兴盛期向退行期转换。     

香猪睾丸miR-449b和miR-222分子的发育性变化

为了阐明miR-449b和miR-222对贵州从江香猪睾丸发育的调节作用。运用qPCR方法检测香猪睾丸组织中miR-449b和miR-222表达量的发育性变化,通过miRSystem数据库、Kobas数据库对预测的靶基因进行信号转导通路富集分析。经克隆测序,得到香猪miR-222序列,与猪的相似性为100%,而且物种之间miR-449b和miR-222两个分子的保守性都很高。香猪睾丸中miR-449b于1~4月龄期间的表达量逐渐上升,4月龄时表达量达到峰值,之后下降;miR-222的表达量呈现相反的变化趋势,于1~4月龄期间表达量逐渐下降,1月龄时最高,4月龄时最低,之后的表达量有所回升;已知物种的miR-449b和miR-222序列高度保守;经预测miR-449b靶基因有420个,富集于卵母细胞减数分裂、Notch信号通路、Ras信号通路等100个途径(P0.05),miR-222的靶基因有429个,富集于细胞周期、Fox O信号通路、TGF-β信号通路等69个通路(P0.05);miR-449b和miR-222的靶基因共享的通路为Wnt信号通路、MAPK信号通路、雌激素信号通路等42个(P0.05)。     

苏博美利奴羊毛囊发育不同时期ncRNA调控网络的构建

为探索苏博美利奴羊毛囊发育相关ncRNA分子调控网络，本研究通过挖掘苏博美利奴羊皮肤组织样本不同时期转录组数据，利用生物信息学方法构建lncRNA-miRNA-mRNA分子调控网络。结果表明：不同时期mRNA与lncRNA及miRNA具有复杂的网络连接。通过GO富集分析发现靶基因主要参与细胞增殖调控、凋亡过程及毛囊形态发生等生物学过程；参与细胞外空间、细胞外泌体、转录因子复合体、质膜的整体成分等细胞组分；参与转录因子活性、结构分子活性、细胞外基质结构组成等分子功能。通过KEGG通路分析发现差异的靶基因主要参与TGF-beta信号通路、癌症途径、MAPK信号通路、PI3K-Akt信号通路、细胞因子与细胞因子受体的相互作用等多种途径。筛选寻找调控网络中的关键基因，发现G2 vs.G1组中miR-133、miR-433-5p的靶基因TGFβ2和NOTCH1参与毛囊形态发生（GO:0031069）的生物学过程，因此TGFβ2和NOTCH1为毛囊发育的关键基因。通过构建苏博美利奴羊毛囊发育相关ncRNA分子调控网络，解析毛囊发生发育的分子调控机制，从转录组学层面对毛囊生长发育进行探讨，能够对苏博...     

EDA基因及其信号通路在动物皮肤毛囊发育中作用研究进展

EDA基因的突变会引起多种外胚层附属物的减少或是缺乏,像头发、牙齿、汗腺、皮脂腺和乳腺等外胚层附属物的生长发育,都需要EDA的参与。EDA信号通路在多种类型毛囊的发生、毛干的形成和皮脂腺的形态发生中发挥作用,它可控制胚胎上皮细胞的命运,调控毛囊细胞的分化过程。作者简述了EDA基因及其信号通路与动物毛囊发育的关系,介绍了EDA基因的功能和蛋白分子结构,以及它主要的两种转录体EDA1和EDA2的功能与作用机制,EDA1-EDAR系统在皮肤衍生结构的发育中发挥主要作用,EDA1和EDAR的突变均会导致HED。阐述了EDA基因信号通路调控毛囊组织生长发育机理以及与多个信号通路的关联,EDA信号可控制毛囊生长期到退行期的转换,对毛囊循环周期中退行期毛囊角质细胞的细胞凋亡起调控作用。EDA-EDAR系统是对皮肤附属器官发育有重要作用的Wnt信号通路的一个重要的效应器,其紧接着出现在Wnt诱导信号的下游。它依靠下游的NF-κB通路的活动来调控靶基因,通过刺激NF-κB调整Wnt、SHH、FGF、和TGF-β信号通路的效应器和抑制因子的转录来发挥作用,调控上皮细胞和间充质细胞还包括组织内的相互作用。结论认为EDA信号通路参与毛囊细胞的发育与分化调控,对动物毛囊生长发育起重要影响作用。     

牦牛miR-381的靶基因预测及生物信息学分析

旨在对牦牛miR-381的靶基因进行预测和生物信息学分析,探索其影响牦牛肌肉发育的作用机制。将高通量测序获得的牦牛miR-381序列与miRbase数据库中已有物种miR-381成熟序列的保守性进行比对分析,利用TargetScan、miRDB和miRanda对miR-381靶基因进行预测,交集的基因与miRTarbase数据库中已证实的靶基因合并后进行基因本体论富集和信号通路富集分析。结果显示,miR-381序列在各物种间高度保守,其靶基因主要参与转录调控过程、横纹肌发育的负调节、细胞周期、心脏发育、肌肉收缩和细胞增殖等生物学过程。信号通路分析结果显示靶基因显著富集于Wnt、TGF-β、MAPK和Notch信号通路等与肌肉发育相关的通路中。由此推测,miR-381可能通过抑制Wnt、MAPK和TGF-β等信号通路的靶基因,进而调控牦牛骨骼肌的发育。研究结果将为进一步探讨miR-381在牦牛肌肉发育中的作用奠定理论基础,同时也为解析牦牛肌肉发育的分子机制提供新的思路。     

绒山羊毛囊发育中非编码RNA的研究进展

非编码RNA(non-coding RNA,ncRNA)是指从基因组上转录而来,不编码蛋白质,但在RNA水平上能行使一定生物学功能的RNA。非编码RNA的种类较多,种类的不同造成功能的差异。毛囊是绒山羊皮肤特殊的附属物,位于皮肤的真皮层,发生于绒山羊胚胎期。由于真皮细胞和上皮细胞间的信号分子传递,使其发育呈周期性变化,历经生长期、退行期和休止期。近年来,有关ncRNA在绒山羊毛囊发育中作用报道较多。文章就其中的微小RNA(miRNA)、长链非编码RNA(lncRNA)及环状RNA(circRNA)的生物发生及其在绒山羊毛囊发育中的作用机制进行了归纳、总结。miRNA是真核生物中存在的一种长18～25 nt的ncRNA分子,可通过与靶信使核糖核酸mRNA特异结合,抑制转录后基因表达。在绒山羊毛囊发育中,miRNA通过抑制相关基因及相关信号通路中关键分子的表达而调控毛囊的发育周期以及毛囊的再生能力。lncRNA是长度超过200 nt的ncRNA,经过剪接,具有polyA尾巴与启动子结构,在绒山羊毛囊发育中能促进毛囊细胞增殖与分化,通过调控基因的靶向表达与多种信号通路互作调节毛囊发育及绒毛生...     

苏博美利奴羊皮肤组织中WNT2基因的DNA甲基化与表达量分析

毛囊是皮肤的衍生物,WNT信号通路是已知与毛囊发育密切相关的通路之一,为研究WNT信号通路中关键基因WNT2的DNA甲基化和基因表达对羊毛性状的调控作用,本研究以周岁苏博美利奴羊为研究对象,采用亚硫酸氢盐测序法(BSP),并结合前期MeDIP-seq结果检测了WNT2基因在第5外显子及前后100bp处的DNA甲基化模式,通过实时荧光定量PCR法检测WNT2基因在苏博美利奴羊皮肤组织中的表达量。结果显示:WNT2基因在2组中的甲基化水平均较高,且2组的表达水平与DNA甲基化水平变化趋势一致,同时,CpG11位点的甲基化水平与表达量之间呈现显著正相关,该位点可能与羊毛纤维直径有关,表明WNT2基因的DNA甲基化水平对羊毛纤维直径有一定作用。本实验为后期深入研究DNA甲基化对毛囊的影响提供了重要参考,也为该基因的进一步研究提供理论基础。     

细毛羊皮肤毛囊miR-1298-5p的靶基因预测及验证

为探讨miR-1298-5p及其靶基因与细毛羊毛囊生长发育的关系,本实验以乾华肉用美利奴羊不同发育时期皮肤毛囊组织为材料,通过生物信息学方法预测miR-1298-5p的靶基因,采用RT-PCR和Western Blot对miR-1298-5p及其靶基因进行定量检测;扩增靶基因3'UTR区,构建野生型和突变型靶基因表达载体后与miR-1298-5p mimics共同转染到HEK293T/17细胞中,利用双荧光素酶报告基因实验验证miR-1298-5p与靶基因的靶向关系。结果表明:FGF_2基因3'UTR区含有miR-1298-5p结合位点,FGF_2为miR-1298-5p潜在靶基因,通过核酸水平和蛋白水平定量表达规律发现miR-1298-5p与靶基因FGF_2呈负相关,初步确定FGF_2为miR-1298-5p的靶基因。双荧光素酶报告基因实验显示,miR-1298-5p mimics能够抑制野生型FGF_2的靶位点报告载体的荧光素酶活性;但当靶基因结合位点突变后,FGF_2被抑制的荧光素酶活性能够得到恢复。综上,miR-1298-5p与FGF_2有靶向关系,miR-1298-5p能通过负调控FGF_2的表达调控羊毛囊发育。     

基于表达谱芯片筛选鸡不同部位皮肤组织差异表达基因

旨在筛选鸡不同部位(背部和腿部)皮肤组织中的差异表达基因,探讨二者皮肤毛囊出现表型差异的分子机制。以处于上市日龄的(63日龄)商品代肉鸡为素材,利用Agilent鸡全基因组表达谱芯片对皮肤厚度、毛囊密度和直径出现显著差异的背部和腿部皮肤组织mRNA表达水平进行检测,对与皮肤毛囊发生发育相关的候选基因及信号通路进行系统筛查,并运用实时荧光定量PCR技术对部分差异表达基因进行验证。结果显示,芯片分析共筛选到差异表达基因676个(|FC|≥2,P0.05),以腿部皮肤为参照,背部皮肤组织中上调基因223个,下调基因453个。荧光定量PCR验证部分差异表达基因的表达趋势与芯片结果基本一致。GO分析表明,差异表达基因参与细胞增殖和凋亡调控、Wnt信号通路、神经元分化、细胞间黏附调节、细胞命运的确定等生物学过程。KEGG信号通路分析发现,除了部分已知的与皮肤毛囊生长发育相关的信号通路(Wnt和Hedgehog信号通路),ECM受体相互作用、黏着、细胞黏附分子、黏合连接等信号通路也被富集。蛋白互作网络分析表明,这些差异表达基因相互作用形成一个调控网络,可能通过影响皮肤毛囊的生长发育和周期变化来影响皮肤毛囊的性状,推测ECM受体相互作用、黏着、运输等调控通路及DKK1、WNT3A、SHH、FGF1、IGF1等基因可能是参与鸡不同部位皮肤毛囊性状出现差异的重要调控通路和候选基因。     

外源褪黑激素（MT）对内蒙古绒山羊皮肤中BMP2基因表达的影响

研究旨在检测骨形态发生蛋白2（bone morphogenetic protein2,BMP2）基因在内蒙古绒山羊皮肤毛囊发育一个完整周期内的表达规律及体内埋植褪黑激素（Melatonin,MT）后对其表达的影响。选取背景相同的内蒙古绒山羊成年母羊8只,随机分为4组,其中1～3组为埋植MT的试验组,第4组为不埋植MT的对照组,采集12个月的体侧部皮肤样品,通过实时定量PCR法检测BMP2基因在发育毛囊周期内的表达量。结果表明,在毛囊发育兴盛期和休止期,BMP2基因的表达量呈缓慢增加。且休止期的表达量高于生长期的表达量;埋植MT后降低了该基因在休止期的表达量,且在整个生长周期内表达量差异不显著。内蒙古绒山羊皮肤中BMP2基因的表达量与毛囊周期性生长密切相关,埋植MT后改变了BMP2基因在毛囊生长周期过程中的表达模式,推测该基因参与调节内蒙古绒山羊皮肤毛囊生长发育。     

成年滩羊和小尾寒羊皮肤毛囊差异表达miRNA的筛选

为比较miRNA在小尾寒羊和滩羊皮肤毛囊中的表达差异,本研究利用高通量测序技术分析了成年滩羊(TY_1)和小尾寒羊(XWHY_1)皮肤毛囊组织中miRNAs的表达谱,在2个品种绵羊皮肤毛囊组织中共鉴定出561个miRNAs,其中包括138个已知和423个新发现的miRNAs。鉴定出的miRNAs进行表达量差异分析发现,在TY_1与XWHY_1共筛选到63个上调和16个下调的miRNAs。对差异表达miRNA靶基因预测后与基因本体数据库(GO)和京都基因与基因组百科全书(KEGG)比对,分别获得靶基因的注释信息为3 886个和4 449个。GO统计发现,差异表达miRNA的靶基因主要参与代谢过程、催化活性、细胞进程和细胞组分等;而KEGG通路分析表明,4 449个靶基因富集到113个信号通路上,其中在嘌呤代谢、内吞作用和糖酵解/糖异生等信号通路上富集显著。综上,在小尾寒羊和滩羊皮肤毛囊中筛选到的差异表达miRNA可能通过调控其靶基因最终参与了绵羊皮肤毛囊的发育。     

黄芪多糖对种公鸡不同组织miR-16表达的影响及其功能预测分析

本试验研究黄芪多糖(APS)对种公鸡不同组织中miR-16表达的影响,结合miR-16预测靶基因的GO和KEGG功能富集性分析,旨在探讨APS对不同组织功能的潜在调控作用。选取1日龄科宝500父母代肉种鸡64只,随机分为对照组和APS组,每组4个重复,每个重复8只鸡。对照组饲喂玉米-豆粕型基础饲粮,APS组在基础饲粮中添加10 g/kg APS,试验期40周。结果表明:1)相比于对照组,APS组可以提高种公鸡的采精量和有效精子数,降低精子畸形率,但差异不显著(P0.05)。2)miR-16具有组织差异性,APS可显著或极显著上调肝脏、脾脏、十二指肠黏膜、回肠黏膜miR-16的表达(P0.05或P0.01),显著降低睾丸miR-16的表达(P0.05)。3)对miR-16预测靶基因进行GO富集分析,表明预测靶基因在跨膜运输、泛素依赖的蛋白质降解、胞内蛋白运输、糖酵解等物质代谢相关生物过程显著或极显著富集(P0.05或P0.01)。KEGG富集分析表明,miR-16靶基因在黏着斑、胰岛素信号通路、糖酵解/糖异生等与机体物质代谢、细胞增殖分化相关通路显著富集(P0.05),在Toll样受体信号通路、自然杀伤细胞介导的细胞毒性等免疫相关通路也有富集。由此可见,APS上调种公鸡脾脏、肝脏、十二指肠黏膜、回肠黏膜组织的miR-16表达,下调睾丸miR-16表达;miR-16预测靶基因与物质代谢、免疫调控密切相关。APS通过调控种公鸡不同组织miR-16表达,影响机体物质代谢和免疫调控功能。