白来航鸡STING基因克隆、生物信息学及组织表达特性分析

【目的】克隆白来航鸡干扰素基因刺激因子基因(stimulator of interferon genes, STING)CDS区序列并进行生物信息学和组织表达分析,为阐明STING基因在抗病毒免疫应答中的作用奠定基础。【方法】采用PCR扩增并克隆白来航鸡STING基因CDS区,测序后对其编码氨基酸序列进行相似性比对及系统进化树构建,利用生物信息学预测STING蛋白的理化特性及结构功能,并利用实时荧光定量PCR技术检测STING基因在鸡心脏、肝脏等14个组织中的表达情况。【结果】白来航鸡STING基因CDS区序列全长1 140 bp,编码379个氨基酸。相似性比对和系统进化树分析结果表明,白来航鸡STING基因与原鸡的相似性最高(99.7%),亲缘关系最近,与冠小嘴乌鸦亲缘关系最远。STING蛋白为酸性、亲水性蛋白,分子质量为42.625 ku,等电点(pI)为6.67,不稳定系数为69.26,脂肪系数为105.01。该蛋白大部分在线粒体和内质网上合成,含有跨膜结构,不含信号肽。STING蛋白二级结构包括α-螺旋(54.62%)、延伸链(10.29%)、β-转角(3.43%)及无规则卷曲...     

无量山乌骨鸡HMOX1基因克隆、生物信息学分析及组织表达研究

【目的】扩增无量山乌骨鸡血红素加氧酶1(heme oxygenase 1,HMOX1)基因CDS区并进行生物信息学分析,检测其组织表达特征,为后续该基因的功能研究奠定基础。【方法】以无量山乌骨鸡皮下脂肪组织cDNA为模板,PCR扩增HMOX1基因,连接PESI-T载体后进行测序分析,使用DNAMAN软件对所获序列进行翻译并与原鸡进行序列比对;使用BLAST、Mega 7.0、SOPMA、ProtParam和ExPASy等在线软件进行相似性比对、系统进化树构建、编码蛋白的理化性质及蛋白结构预测;利用STRING软件预测HMOX1互作蛋白;通过实时荧光定量PCR检测HMOX1基因在无量山乌骨鸡不同组织中的表达水平。【结果】无量山乌骨鸡HMOX1基因CDS区长为891 bp,编码296个氨基酸,与原鸡HMOX1基因(登录号：NM＿205344.1)核苷酸序列相比存在3处突变,其中2处为同义突变,第217 bp处G>C为错义突变,导致第73位的天冬氨酸突变为甘氨酸。HMOX1基因与原鸡和日本鹌鹑的氨基酸序列相似性较高(99.6%和95.0%),与马的相似性最低(60.2%)。系统进化树分...     

家鸽及边鸡的MC1R基因片段的克隆及分析

为了研究黑素皮质素受体1(MC1R)基因与家鸡、家鸽羽色的相关关系,根据Gen Bank上原鸡MC1R基因序列,利用在线引物设计软件Primer3(V.0.4.0)在编码区设计1对引物PE,试验以家鸽和边鸡DNA为模板,PCR扩增,测序,分别获得长度为303 bp和335 bp的基因片段。结果表明:家鸽与原鸡MC1R基因核苷酸序列存在10处变异,边鸡与原鸡存在3处变异,其序列相似性分别为95.71%、99.10%。二级结构分析结果表明,家鸽和边鸡α螺旋、β折叠、无规则卷曲所占比例分别为47.92%、31.25%、20.83%,64.80%、21.90%、13.30%,均有2个跨膜区。针对2个物种氨基酸序列同源与变异预测其三级结构,结果表明,MC1R基因可能是表达禽类羽色的主效基因。     

静原鸡DERA基因生物信息学分析及组织表达研究

【目的】分析2-脱氧-D-核糖-5-磷酸醛缩酶(2-deoxyriboaldolase, DERA)CDS区序列特征及其在静原鸡不同组织和不同生长发育时期的表达水平,为进一步研究DERA基因在静原鸡肌肉中的调控机制提供参考。【方法】以42日龄静原鸡胸肌cDNA为模板,克隆DERA基因完整CDS区,测序并进行生物信息学分析。利用实时荧光定量PCR技术检测DERA基因在静原鸡不同生长发育阶段(7、42、126、180日龄)胸肌、腿肌和42日龄不同组织(心脏、脾脏、肝脏、肾脏、肺脏、胸肌和腿肌)中的表达水平。【结果】静原鸡DERA基因CDS区全长699 bp,编码232个氨基酸,与环颈雉的相似性最高,为99.1%,表明DERA蛋白氨基酸序列在进化过程中保守性较高。DERA属于稳定蛋白,平均亲水系数为0.175,属于疏水性蛋白;存在跨膜结构,没有信号肽,属于跨膜蛋白。蛋白结构域是一种deoC/LacD家族醛缩酶域,二级结构包含α-螺旋、延伸链和无规则卷曲。蛋白互作网络分析表明,DERA蛋白与PGM2、TPI1、TKT等蛋白存在互作关系。DERA基因在静原鸡肾脏、胸肌、腿肌、心脏、肝脏、肺脏、脾...     

太行鸡FST基因克隆、序列特征及表达分析

旨在基于生物信息学方法分析太行鸡卵泡抑素基因（follistatin，FST），并研究其在不同组织及卵泡不同发育时期表达差异，解析FST对太行鸡产蛋性能的影响。本研究使用RT-PCR对43周龄健康太行鸡FST基因编码区序列进行克隆和测序，运用生物信息学软件对其功能结构进行预测，荧光定量PCR技术对43周龄健康麻羽太行鸡FST在不同组织及不同发育时期卵泡中的表达差异进行分析，每组设3个重复，进行3次平行试验。结果显示，太行鸡FST基因编码区为1 032 bp，编码343个氨基酸，太行鸡FST蛋白与日原鸡同源性最高（100.0%），存在Sec信号肽，成熟蛋白具有FOLN和KAZAL_FS 2个保守结构域。该蛋白为亲水蛋白，蛋白分子式为C₁₆₁₉H₂₅₆₇N₄₅₉O₅₂₀S₄₄，分子量为38.192 6 ku，理论等电点为5.59，其中含量最高的是Cys （C），为10.50%；蛋白质二级结构含有α-螺旋、延伸链、β转角、无规卷曲，占比分别为16.33%、16.62%、5.25%、61.81%。FST基因在不同组织中均有表达，在肝中表达量最高。FST基因在各级卵泡表达量存在差异，其中，在大白卵泡中表达量最高（P<0.01），在卵泡选择前后表达量存在1.6倍表达差异。上述结果为研究太行鸡FST基因结构和功能提供了重要的基础数据，也为进一步挖掘太行鸡产蛋性能候选高产基因提供了参考。     

水牛TRAIL基因克隆及生物信息学分析

【目的】对水牛肿瘤坏死因子(tumor necrosis factor,TNF)相关凋亡诱导配体(TNF-related apotosis-inducing ligand,TRAIL)基因CDS序列进行克隆及序列分析,并对其编码的蛋白进行生物信息学分析,为后期TRAIL蛋白调控水牛卵巢卵泡发育、颗粒细胞增殖及凋亡的研究奠定基础。【方法】利用RT-PCR方法克隆水牛TRAIL基因CDS序列,对所获序列进行核苷酸序列、氨基酸序列相似性比对,构建系统进化树,并通过生物信息学软件分析TRAIL基因编码蛋白的结构和功能。【结果】试验成功克隆水牛TRAIL基因CDS序列,长864 bp,编码287个氨基酸;水牛TRAIL基因与牦牛、普通牛、山羊、绵羊、野猪、马、人、黑猩猩和家鼠的核苷酸序列相似性分别为99.2%、99.3%、95.9%、96.3%、84.7%、84.8%、81.3%、81.3%和70.0%。系统进化树结果表明,水牛与牦牛、普通牛的亲缘关系最近,与家鼠亲缘关系最远。氨基酸序列比对结果表明,在不同物种间,其跨膜结构域和TNF结构域序列保守性较高。TRAIL蛋白属于亲水性蛋白,存在1个跨膜结构域,140―285位氨基酸处为TNF区,具有29个磷酸化位点,无信号肽和糖基化位点,主要定位于细胞质中。TRAIL蛋白二级结构主要以无规则卷曲为主,约占51.57%,其次为延伸链(24.39%)和α-螺旋(24.04%)。TRAIL蛋白三级结构与二级结构一致,且与模型蛋白人TRAIL蛋白的相似性为75.53%。【结论】本试验克隆得到水牛TRAIL基因CDS区序列,大小为864 bp,编码287个氨基酸,水牛与牦牛、普通牛亲缘关系最近,TRAIL蛋白跨膜结构域和TNF结构域在不同物种间序列保守性较高,这可能与其功能有关。     

藏羊MSTN基因克隆、生物信息学及组织表达分析

【目的】对藏羊肌肉生长抑制素(myostatin, MSTN)基因进行克隆和生物信息学分析,检测其在藏羊不同组织中的表达,为探究MSTN基因在藏羊中的生物学功能提供参考。【方法】以藏羊背最长肌组织cDNA为模板,克隆藏羊MSTN基因完整CDS区序列并测序,用SeqMan程序对测序结果进行拼接,并用BLAST在线程序对组装后的序列进行分析鉴定。用生物信息学软件进行相似性比对、系统进化树构建及生物信息学分析,用实时荧光定量PCR检测MSTN基因在藏羊不同组织中的表达量。【结果】藏羊MSTN基因CDS全长为1 128 bp,编码375个氨基酸。藏羊MSTN基因氨基酸序列与绵羊、牦牛、牛、猪、恒河猴、人、黑猩猩、犬、鸡及斑马的相似性依次为100.0%、93.4%、93.4%、95.2%、94.4%、94.2%、94.4%、93.1%、87.8%和87.5%;系统进化树分析结果表明,藏羊与绵羊的亲缘关系最近,与斑马和鸡的亲缘关系最远。藏羊MSTN蛋白属于亲水性分泌蛋白,且具有不稳定性,不含跨膜结构,含1个信号肽,存在31个潜在的磷酸化位点、2个N-糖基化修饰位点,主要分布在线粒体和细胞质中;MS...     

大围山微型鸡POU1F1基因编码区的克隆及序列分析

大围山微型鸡是我国目前报道的体形最小、体重最轻的地方优良品种,具有较高的屠宰率和饲料报酬。研究通过设计特定引物对大围山微型鸡POU1F1基因的编码区进行PCR分段扩增并测序。利用分子生物学软件对其基因序列及编码产物的结构、功能进行了生物信息学分析,并构建其系统发育关系。结果表明,大围山微型鸡POU1F1基因编码区全长984bp,与原鸡的POU1F1基因（NM〈sub〉2〈/sub〉04319）比对,在836位处发生一次碱基转换（C/T）,但未造成氨基酸改变;其编码产物由327个氨基酸残基组成,分子质量为38801.19Da,理论等电点为8.67;系统发育树结果表明,不同物种间在该基因编码区核苷酸序列和氨基酸序列上有较高的相似性,大围山微型鸡与原鸡POU1F1基因编码区核苷酸序列和氨基酸的相似性分别达到了99%和100%。大围山微型鸡POU1F1基因编码区的成功克隆,为进一步研究POU1F1基因的遗传特性和大围山微型鸡的矮小机理奠定了基础。     

白来航鸡半乳糖凝集素-1基因克隆、生物信息学及组织表达特性分析

【目的】克隆白来航鸡半乳糖凝集素-1(galectin-1,Gal-1)基因,对其编码蛋白进行生物信息学分析,并检测其在不同组织中的表达情况,为进一步阐明其抗病毒功能提供科学依据。【方法】以鸡脾脏cDNA为模板,通过PCR扩增鸡Gal-1基因完整CDS区序列,并进行相似性比对及系统进化树构建;运用生物信息学软件对其编码蛋白的理化性质、亲/疏水性、跨膜区、信号肽、修饰结构、保守结构域及高级结构进行预测。利用实时荧光定量PCR检测Gal-1基因在白来航鸡心脏、肝脏、脾脏、肺脏、肾脏、脑、腺胃、肌胃、十二指肠、空肠、盲肠、直肠、胸肌和腿肌组织中的表达情况。【结果】白来航鸡Gal-1基因CDS区序列长度为408 bp,编码135个氨基酸。相似性比对结果表明,白来航鸡Gal-1基因核苷酸序列与火鸡、绿头鸭和珍珠鸟的相似性分别为97.1%、88.4%和82.2%;系统进化树结果表明,白来航鸡与火鸡亲缘关系最近。Gal-1蛋白分子质量为15.06 ku,理论等电点为6.57,不稳定系数为36.05,脂肪系数为74.30,平均亲水指数为-0.259。Gal-1蛋白无信号肽,不存在跨膜区;存在2个明显的...     

琅琊鸡胚MSX2基因克隆及发育性表达变化

【目的】 克隆鸡Msh同源框2(Msh-homeobox 2，MSX2)基因，并对其进行生物信息学和胚胎期表达模式分析，为进一步研究MSX2基因的结构和功能提供支持。【方法】 对80枚琅琊鸡种蛋进行孵化，分别于孵化期第1、2、3、4、5、6、9、12、15、18天采集样品(1~6胚龄采集整胚，9、12、15、18胚龄分别采集心脏和肝脏)，提取总RNA，利用RT-PCR方法扩增并克隆MSX2基因，对其进行生物信息学分析;利用实时荧光定量PCR技术分析MSX2基因在琅琊鸡鸡胚和组织中的表达水平。【结果】 成功克隆了琅琊鸡胚MSX2基因，序列长度为818 bp，开放阅读框为780 bp，编码259个氨基酸。琅琊鸡MSX2氨基酸序列与日本鹌鹑、秃鹰和仓鸮的相似性最高，均为99.23%，与山羊的相似性最低，为74.54%;系统进化树分析结果显示，琅琊鸡与日本鹌鹑聚为一类。生物信息学分析表明，MSX2蛋白分子质量为28 235.18 u，等电点(pI)为9.71，没有信号肽和跨膜域，为不稳定的亲水性蛋白，主要分布在细胞核(60.9%);存在37个磷酸化位点、8个O-糖基化位点和1个N-糖基化位点;二级结构主要由无规则卷曲(63.32%)和α-螺旋(28.19%)组成。实时荧光定量PCR分析表明，MSX2基因在整胚(1~6胚龄)中的表达水平呈现先上升后下降趋势，且在4胚龄时表达水平最高，极显著高于1、2、3胚龄(P<0.01);MSX2基因在12胚龄鸡心脏中的表达水平极显著低于9胚龄(P<0.01);在9~18胚龄鸡肝脏中的表达量呈逐渐下降趋势，且在9胚龄鸡肝脏中表达水平最高(P<0.05)。【结论】 试验成功克隆了琅琊鸡MSX2基因序列，其表达模式在不同胚龄和同一胚龄不同组织间存在差异，为进一步探索琅琊鸡MSX2基因的结构和功能奠定了基础。     

鸡Wnt基因家族生物信息学分析及表达

【目的】试验旨在对鸡Wnt基因家族进行全面系统地调查分析,并进一步明确鸡感染新城疫病毒(Newcastle disease virus, NDV)后Wnt基因家族成员的表达模式。【方法】基于鸡基因组测序结果,通过查询Pfam和UniPort数据库获取鸡Wnt基因家族全基因组序列,构建系统进化树,分析Wnt蛋白理化性质、亚细胞蛋白定位、染色体分布、保守结构域等,并进一步研究鸡感染NDV后Wnt基因家族成员的表达模式。【结果】试验共鉴定到39个鸡Wnt基因家族成员。系统进化树分析发现,鸡与人和鼠的Wnt成员相似性较高,均有12个亚家族,即Wnt1、Wnt2(Wnt2a和Wnt2b)、Wnt3a、Wnt4、Wnt5(Wnt5a和Wnt5b)、Wnt6、Wnt7(Wnt7a和Wnt7b)、Wnt8(Wnt8a和Wnt8b)、Wnt9(Wnt9a和Wnt9b)、Wnt10a、Wnt11及Wnt16成员,但鸡缺少Wnt10b。鸡Wnt成员分布于11条染色体上,分别为Chr1、Chr2、Chr4、Chr6、Chr7、Chr12、Chr13、Chr21、Chr26、Chr27和Chr34,成员间出现1...     

绵羊KAP11.1基因克隆、原核表达及皮肤毛囊表达研究

【目的】 对毛囊角蛋白关联蛋白11.1（keratin associated protein 11.1,KAP11.1）基因进行克隆及原核表达,并对KAP11.1基因在不同品种绵羊皮肤毛囊中表达量进行比较,探究KAP11.1基因在南疆地方绵羊品种间表达差异及其对羊毛品质的影响。【方法】 以平原型和田羊、山区型和田羊和卡拉库尔羊体侧皮肤毛囊为研究材料,以GenBank中绵羊KAP11.1基因序列（登录号：HQ595347.1）为参照设计引物,对KAP11.1基因进行PCR扩增,构建pMD19-T-KAP11.1克隆质粒,双酶切鉴定后构建pET-28a （+）-KAP11.1原核重组表达质粒,经PCR和双酶切鉴定后测序并进行序列分析,转化大肠杆菌BL21（DE3）感受态细胞中表达,采用SDS-PAGE和Western blotting检测;利用实时荧光定量PCR技术检测KAP11.1基因在不同绵羊皮肤毛囊中的表达情况。【结果】 3种绵羊KAP11.1基因CDS区序列为480 bp,编码159个氨基酸,为不稳定的疏水蛋白。相似性比对结果发现,与参照基因相比,2种类型和田羊基因序列相似性均为99.79％,均在423 bp处发生突变,由C变为T,卡拉库尔羊基因序列相似性为99.38％,其69 bp处G变为T、93 bp处C变为T、423 bp处C变为T。系统进化树分析发现,3种绵羊和山羊亲缘关系最近,和瘤牛亲缘关系最远。KAP11.1蛋白二级结构主要由无规则卷曲组成。试验成功构建了pET-28a （+）-KAP11.1原核重组表达质粒,并纯化得到19 ku的KAP11.1蛋白。KAP11.1基因在山区型和田羊和卡拉库尔羊皮肤毛囊中的表达量均显著高于平原型和田羊（P<0.05）,在山区型和田羊和卡拉库尔羊中差异不显著（P>0.05）。【结论】 克隆获得480 bp的绵羊KAP11.1基因CDS区序列,成功构建了pET-28a （+）-KAP11.1原核重组表达质粒,并获得19 ku的KAP11.1蛋白,且KAP11.1基因在3个品种绵羊皮肤毛囊中均有表达。     

黔东南小香鸡MEF2C基因外显子多态性及其与体重、体尺性状的关联分析

【目的】探究肌细胞增强因子2C(myocyte enhancer factor 2C,MEF2C)基因多态性对黔东南小香鸡体重和体尺性状的影响。【方法】试验以251羽300日龄黔东南小香鸡为研究对象,利用PCR扩增和直接测序技术对MEF2C基因全部(15个)外显子进行单核苷酸多态性(single nucleotide polymorphism, SNP)位点筛选,并运用SPSS 18.0统计分析软件对不同SNPs位点所对应的基因型与黔东南小香鸡的体重、体尺性状进行相关性分析。【结果】黔东南小香鸡群体中仅在MEF2C基因外显子12上发现3个SNPs位点：g.1819 T>C、g.2426 T>C和g.2543 C>T,在其他外显子中未发现多态位点。χ²检验结果表明,3个SNPs在黔东南小香鸡群体中均偏离Hardy-Weinberg平衡状态(P>0.05)。关联分析结果显示,g.1819 T>C位点TC基因型个体体斜长、胸深和胫围均极显著高于TT和CC基因型,体重极显著高于TT基因型(P<0.01);TC基因型个体胫长显著高于TT和...     

鸡繁殖性状近交衰退相关lncRNA的筛选

【目的】 筛选与鸡繁殖性状近交衰退相关的长链非编码RNA （long non-coding RNA，lncRNA），为深入探究lncRNA对鸡繁殖性能近交衰退调控机制提供参考。【方法】 以狼山鸡高、低近交群体为试验素材，通过转录组测序技术分析狼山鸡下丘脑和卵巢中lncRNA表达情况，筛选高、低近交组间差异表达lncRNA，并对其进行顺式调控靶基因预测及功能分析。【结果】 狼山鸡下丘脑和卵巢中共鉴定出4 222个lncRNAs，高、低近交组间比较发现，下丘脑中差异表达lncRNAs 35个，卵巢中差异表达lncRNAs 215个。下丘脑中差异lncRNAs中预测到顺式调控靶基因98个，这些靶基因显著富集于出生或孵化时胚胎发育终止、胚胎心导管发育、视黄酸应答等繁殖相关生物过程（P<0.05），涉及lncRNA MSTRG.9196.4、MSTRG.9195.2、MSTRG.6254.2以及相应靶基因DNAAF2和FKBP1B等胚胎发育相关基因；卵巢中差异lncRNAs预测到顺式调控靶基因414个，这些靶基因富集到卵母细胞减数分裂、MAPK、叶酸合成等信号通路，包括MSTRG.7683.1、MSTRG.13604.4、MSTRG.16570.1、MSTRG.8330.5、MSTRG.8330.4等神经内分泌调节及配子生成相关的lncRNAs。【结论】 本研究在下丘脑和卵巢中筛选到了一系列鸡胚胎发育及配子生成调控相关差异lncRNAs，这些lncRNAs可作为鸡繁殖性能近交衰退候选lncRNAs，为进一步揭示鸡繁殖性能近交衰退调控机制提供线索。     

犬ANP32蛋白多态性及其对A型流感病毒RNA聚合酶活性的影响

【目的】探索犬酸性核磷蛋白32(canis acidic (leucine-rich) nuclear phosphoprotein 32 ku, caANP32)对不同物种A型流感病毒(Influenza A virus, IAV) RNA聚合酶活性的影响。【方法】利用实时荧光定量PCR方法分析caANP32家族基因(caANP32A、caANP32B及caANP32E)的组织分布并对其进行扩增和克隆,评估caANP32家族蛋白对不同物种A型流感病毒RNA聚合酶活性的支持作用。【结果】caANP32家族基因在不同组织中分布相似,其中caANP32B基因组织丰度显著或极显著高于caANP32A和caANP32E基因(P<0.05;P<0.01),在心脏、盲肠和大脑中caANP32E基因组织丰度极显著或显著高于caANP32A基因(P<0.01;P<0.05),在肝脏、肺脏等组织中caANP32E与caANP32A基因丰度均无显著差异(P>0.05)。在犬心脏组织和MDCK细胞中发现,caANP32A基因仅存在1个转录本,caANP32B基因存在3个不同剪接...     

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

【目的】通过对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个生长发育相关基因,分别为肌球蛋白重...     

儋州鸡体尺性状全基因组关联分析

【目的】 挖掘影响地方鸡体尺性状的有效SNP位点及功能基因, 给儋州鸡育种工作提供有效的数据基础和理论支撑。【方法】 共采集200只儋州鸡血样并提取基因组DNA, 利用10×全基因组重测序技术获得全基因组SNP标记并对试验个体基因型进行分型。使用EMMAX软件基于混合线性模型对70日龄的儋州鸡体尺性状(胫长、胫围、体斜长、胸宽、髋骨宽、胸深、龙骨长)进行全基因组关联分析。【结果】 共发现与胫长性状和胫围性状基因组水平显著相关的SNPs位点有12和8个, 与胫长性状相关SNPs分别定位于1、2、4和8号染色体上; 与胫围性状相关的SNPs定位于2、4、8和13号染色体上。预测与胫长相关的候选基因为KCNA1、TPK1、EZH2、FSTL5和AMY2A基因, 与胫围相关的候选基因为TPK1、FSTL5、AMY2A、TGFBI、LECT2和IL-9。通过KEGG通路分析和GO注释发现, 8个基因参与钾离子跨膜转运、硫胺素新陈代谢、细胞增殖、钙离子结合、骨骼肌卫星细胞维持与骨骼肌再生、细胞受体相互作用、生长因子活性等生物学进程。【结论】 本研究发现了20个与儋州鸡体尺性状关联的SNPs位点, 并筛选到8个目标性状候选基因, 为儋州鸡育种提供候选的分子标记, 为地方鸡标记辅助选择提供新的思路。