2. 河南省家禽育种国际联合实验室，郑州 450046;
3. 河南省鸡种质资源创新与利用重点实验室，郑州 450046
, XING Yuxin1, MA Chenglin1, GUAN Hongbo1, JIA Qihui1, KANG Xiangtao1,2,3, TIAN Yadong1,2,3, LI Zhuanjian1,2,3, LIU Xiaojun1,2,3
, LI Hong1,2,3
2. International Joint Research Laboratory for Poultry Breeding of Henan, Zhengzhou 450046, China;
3. Henan Key Laboratory for Innovation and Utilization of Chicken Germplasm Resources, Zhengzhou 450046, China
成纤维细胞生长因子6(fibroblast growth factor 6, FGF6)是成纤维细胞生长因子家族(FGFs)的成员之一，主要通过与成纤维细胞生长因子受体(fibroblast growth factor receptor, FGFR)1和4结合后发挥其生物学活性[1-2]。有研究表明，FGF6是血管新生诱导因子之一，它可以促进血管内皮细胞的生成[3]。沉默FGF6会抑制膀胱癌(BC)细胞系有氧糖酵解，从而调节PI3K/Akt和MAPK信号通路抑制其血管生成[4]。FGF6基因可调节人舌头发育过程中的肌源性分化和成肌细胞融合[5]。Cai等[6]通过构建肌母细胞和周围神经损伤的体内模型，发现FGF6在去神经支配的大鼠肌肉中增强了MyHC-IIb肌纤维表型，减少了肌肉萎缩，进而促进神经支配后再生骨骼肌的功能恢复。Xu等[7]通过对男性骨骼肌样本进行甲基化DNA免疫沉淀测序(MeDIP-Seq)和RNA测序(RNA-Seq)发现FGF6基因表达下调，骨骼肌中过表达FGF6能够改善胰岛素抵抗，增加骨骼肌的AMPK活性，进而增加全身葡萄糖量和脂质代谢。Shamsi等[8]发现, FGF6作为一种脂肪因子可独立于脂肪生成机制调节解偶联蛋白1(uncoupling protein 1，UCP1)的表达，并调节全身能量代谢。Liu等[9]研究表明，FGF6是脂肪细胞祖细胞(APC)维持脂肪稳态和胰岛素敏感性的增殖因子。Hu等[10]研究发现, 重组人FGF6蛋白能够通过抑制Hippo途径改善心肌梗死(MI)小鼠的心脏功能并促进其心脏的修复能力。Creighton等[11]通过对非裔美国人前列腺癌(Pca)转录组学和表观基因组学全DNA甲基化的综合分析，发现了与Pca唯一相关的包括FGF6在内的12个候选基因。
上述关于FGF6基因的报道主要集中在哺乳动物上，表明该基因在血管形成、骨骼肌再生、葡萄糖代谢、脂肪生成及稳态、心肌发育和肿瘤细胞的发育等生理和病理过程中具有重要作用。在鸡上，Kumar和Chapman[12]对鸡FGF6基因进行了克隆，并通过原位杂交发现其在雏鸡的头部、咽部和内脏中胚层广泛表达。先前研究发现FGF6基因在固始鸡的十二指肠高表达[13]。除此之外，有关鸡FGF6基因的相关研究鲜有报道。
综上，FGF6基因在哺乳动物的细胞增殖、肌肉发育过程中发挥重要作用，在鸡上对其研究较少，尚未见对FGF6基因多态性及与鸡经济性状的关联分析。本研究旨在通过对鸡FGF6基因进行生物信息学分析，并对鸡FGF6基因的SNPs与鸡的生长、屠体、肉品质和血清生化指标性状进行关联分析，筛选与鸡经济性状关联的重要SNPs及单倍型组合，为进一步研究其生物学功能及调控机制提供理论基础，为鸡分子育种提供数据支撑。
1 材料与方法 1.1 试验材料本研究以河南省家禽种质资源创新工程研究中心构建的固始鸡×安卡鸡F2资源群768个个体为研究对象(该群体按F2远缘半同胞设计方案，有正反交7个家系)[14]。所有鸡饲养在同一环境，自由采食饮水，饲喂至第84天时，颈静脉采血，血液置于含/不含EDTA的抗凝管中，用于制备DNA/血清。768只鸡简化基因组测序(genotyping-by-sequencing, GBS)数据[15]及个体所对应的表型数据库由实验室保存。
1.2 测定表型指标固始鸡×安卡鸡F2资源群所有个体从孵化起，每隔2周称体重，每隔4周测定体尺指标。体重指标包括出生重(birth weight, BW0)、2周体重(BW2)、BW4、BW6、BW8、BW10和BW12。体尺指标包括0、4、8和12周龄的胫长(shank length, SL)、胫围(shank girth, SG)、胸深(chest depth, CD)、胸宽(chest width, CW)、胸骨长(breast bone length, BBL)、胸角(breast angle, BA)、体斜长(body slanting length, BSL)和骨盆宽(pelvis breadth, PB)等32个指标[14]。
84日龄屠宰后测定屠体、肉质和血清生化相关指标。屠体性状包括肝脏重(liver weight, LiW)、胸肌重(breast muscle weight, BMW)、腿肌重(leg muscle weight, LMW)、腿重(leg weight, LW)、腹脂重(abdominal fat weight, AFW)、屠体重(carcass weight, CW)、半净膛重(semi-evisceration weight, SEW)、全净膛重(evisceration weight, EW)及其对应的比率等16个指标[15]。肉质性状包含腿肌剪切力(leg shear force, LSF)、胸肌剪切力(breast shear force, BSF)、胸肌失水率(water loss rate of breast muscle, WLRB)和腿肌失水率(water loss rate of leg muscle, WLRL)4个指标[16]。血清生化指标包括谷丙转氨酶(alanine transaminase, ALT)、胆碱酯酶(choline esterase, ChE)、高密度脂蛋白(high density lipoprotein, HDL)、低密度脂蛋白(low density lipoprotein, LDL)、总胆固醇(total cholesterol, CHO)、甘油三酯(triglyceride, TG)等15个指标[17, 18]。
1.3 生物信息学分析从NCBI数据库中下载人(NP_066 276.2)、鼠(NP_571 983.1)、鸡(XP_001 232 071.1)、火鸡(XP_003 202 616.1)、绿海龟(XP_007 055 701.1)、热带爪蟾(NP_001 136 295.1)和斑马鱼(NP_001 001 398.2)FGF6的氨基酸序列进行相似性比对，利用MEGA10.0软件构建系统进化树。利用生物信息学在线软件对FGF6蛋白理化性质及功能进行分析和预测。相关分析项目及对应的工具见表 1。
|
表 1 生物信息学分析在线软件 Table 1 The bioinformatics analysis online softwares |
使用POPGENE Version32构建单倍型并计算SNPs位点的基因型频率、等位基因频率、有效等位基因数(Ne)、杂合度(He)、多态信息含量(PIC)，并检验Hardy-Weinberg平衡。
利用SPSS24.0软件的混合线性模型(最小二乘分析)对FGF6不同基因型和单倍型组合与F2资源群经济性状进行关联分析，数据的多重比较使用邦弗伦尼(Bonferroni)校正。关联分析模型如下：
Yijklm=μ+Gi+Sj+ Hk+fl+eijklm。
式中：Yijklm为个体性状的表型值；μ为性状总体均值；Gi为基因型的固定效应(1~3)；Sj为性别的固定效应(1, 2)；Hk为批次的固定效应(1, 2)；fl为家系的随机效应(1~7)；eijklm为随机误差。P < 0.05表示差异达到显著水平。
2 结果 2.1 FGF6系统进化树和保守性分析从NCBI数据库获得FGF6不同物种的氨基酸序列并进行分析(图 1)，利用CLUSTALW将获得的鸡FGF6基因序列与其他6个物种对比发现，鸡FGF6氨基酸序列与人、鼠、火鸡、非洲爪蟾、斑马鱼和绿海龟的序列相似性分别为65.05%、66.02%、91.26%、63.59%、68.29%、85.44%。不同物种FGF6氨基酸序列利用MEGA10.0软件进行多序列比对后构建系统进化树，结果显示鸡FGF6与火鸡亲缘关系最近，其次为绿海龟，与斑马鱼亲缘关系最远(图 1A)。FGF6保守性分析表明FGF6在不同物种中保守性较高，拥有相同的motif1~4、6基序(图 1B)。
|
A. 不同物种间FGF6系统进化树；B. 不同物种FGF6保守性分析 A. Phylogenetic tree of FGF6 among different species; B. Conservativeness analysis of FGF6 among different species 图 1 不同物种间FGF6氨基酸序列分析 Fig. 1 Amino acid sequence analysis of FGF6 among different species |
利用Expasy-Protparam在线软件对鸡FGF6基因编码的氨基酸序列进行理化性质分析，结果显示FGF6蛋白分子式为C1023H1617N287O298S9，分子量为22 993.38 u，等电点(pI)为9.28，其蛋白不稳定系数(Ⅱ)为35.68，说明该蛋白为稳定蛋白。脂肪系数为89.47，总平均亲水性系数为-0.163，说明FGF6蛋白是亲水性蛋白。信号肽预测结果显示FGF6在第36~37位存在信号肽序列，概率为62.74%，推测其为分泌蛋白(图 2A)。Expasy-ProtScale在线软件计算FGF6蛋白的疏水性，疏水性最大值为3.011(第27位氨基酸)，最小值为-2.400(第70位和71为氨基酸)，该蛋白为亲水性蛋白(图 2B)。FGF6蛋白二级结构分析显示无规则卷曲占比44.17%，α螺旋占26.70%，β转角为8.74%，延伸链为20.39%(图 2C)。FGF6蛋白三级结构预测结果与二级结构预测一致(图 2D)。
|
A.鸡FGF6蛋白质信号肽分析；B.鸡FGF6蛋白质疏水结构预测；C.鸡FGF6蛋白质二级结构分析(h.α-螺旋；t.β-转角；c.无规则卷曲；e.延伸链)；D.鸡FGF6蛋白质三级结构预测 A. Chicken FGF6 protein signal peptide analysis; B. Chicken FGF6 protein hydrophobic structure prediction; C. Chicken FGF6 protein secondary structure analysis (h.α-helix; t.β-turns; c.irregular coiling; e.elongated strand); D. Chicken FGF6 protein tertiary structure prediction 图 2 鸡FGF6蛋白质信号肽、亲疏水性及二、三级结构预测 Fig. 2 The prediction of chicken FGF6 protein signal peptide, hydrophilicity, potential secondary and tertiary structure |
利用TMHMM-2.0进行蛋白质跨膜结构预测，结果显示，FGF6第21~40位氨基酸处存在一个跨膜区域(图 3A)。利用UniProt进行亚细胞定位，GO注释显示FGF6位于细胞质和细胞外空间(图 3B)。利用NetPhos 3.1进行磷酸化预测，结果显示，FGF6蛋白有12个丝氨酸、5个苏氨酸和5个酪氨酸位点发生磷酸化(图 3C)。利用NetNGlyc 1.0进行糖基化预测，结果显示FGF6蛋白第44位氨基酸处有一个N糖基化位点，可能性为74.42%(图 3D)。
|
A.鸡FGF6蛋白质跨膜分析；B.鸡FGF6蛋白质亚细胞定位预测；C.鸡FGF6蛋白磷酸化位点预测；D.鸡FGF6蛋白糖基化位点预测 A. Chicken FGF6 protein transmembrane analysis; B. Prediction of subcellular localization of chicken FGF6 protein; C. Prediction of chicken FGF6 protein phosphorylation sites; D. Prediction of chicken FGF6 protein glycosylation sites 图 3 鸡FGF6蛋白质跨膜结构、亚细胞定位、磷酸化、糖基化位点预测 Fig. 3 Predicted transmembrane structure, subcellular localization, phosphorylation and glycosylation sites of chicken FGF6 protein |
利用STRING在线软件对FGF6的互作蛋白进行预测分析，结果显示FGF6与FGFR1、FGFR2、FGFR3、FGFR4和表皮生长因子(Epidermal growth factor，EGF)等5个蛋白可能存在相互作用(图 4)。
|
图 4 鸡FGF6蛋白与其他蛋白相互作用网络分析 Fig. 4 Network analysis of chicken FGF6 protein interaction with other proteins |
根据GenBank公布的鸡FGF6基因序列(XM_001 232 070.4)，下载获取FGF6基因组启动子区域2 kb(73397575-73410291(-))序列，基于实验室保存的固始鸡-安卡鸡F2资源群的GBS数据(GRCg6a版本)[17]，筛选FGF6基因组包含启动子区域2 kb序列在固始鸡-安卡鸡F2资源群GBS数据中存在的所有SNPs位点，共得到6个遗传变异位点，均位于第1内含子上，见表 2。
|
|
表 2 鸡FGF6基因上SNPs信息 Table 2 SNPs information on chicken FGF6 gene |
对鸡FGF6基因的多态性进行群体遗传学分析表明，在F2资源群中，6个SNPs位点均存在3种基因型，其中rs73399071位点的优势基因型为GG(n=467)，基因型频率为0.650；rs73399075位点的优势基因型为TT(n=590)，基因型频率为0.821；rs73399124位点的优势基因型为CC(n=541)，基因型频率为0.752；rs73399207位点的优势基因型为TT(n=367)，基因型频率为0.510；rs73399282位点的优势基因型为CC(n=438)，基因型频率为0.609；rs73399353位点的优势基因型为TA(n=313)，基因型频率为0.435。除rs73399353外，其余5个位点均处于Hardy-Weinberg平衡状态(P＜0.001)。6个多态位点的有效等位基因数(Ne)介于1.419 3~2.853 2，杂合度(He)介于0.295 4~0.649 5，多态信息含量(PIC)介于0.484 9~1.073 4。各位点的等位基因频率，Ne、He、PIC等群体遗传学指标统计结果见表 3。
|
|
表 3 F2资源群FGF6基因SNPs位点的遗传特性 Table 3 Genetic characterization of the SNPs loci of the FGF6 gene in the F2 resource population |
利用Haploview 4.2软件对rs73399071、rs73399075、rs73399124、rs73399207和rs73399282这5个处于哈代-温伯格平衡位点进行连锁不平衡(LD)分析(图 5)，结果显示，这5个位点之间完全连锁(D′=1，r2=1)(图 5A)，后续使用rs73399071位点进行关联分析；存在6种单倍型，其中GTCCC是主要单倍型，比例为0.257，其次是GTCTT和GTCTC，比例分别是0.202和0.199(图 5B)。
|
A.SNPs连锁不平衡分析；B.单倍型模块分析。颜色越深表示D′越高，数字表示R2 A. SNPs linkage disequilibrium analysis; B. Haplotype module analysis. Darker colors indicate higher D′, numbers indicate R2 图 5 鸡FGF6基因5个位点连锁不平衡分析 Fig. 5 Linkage disequilibrium analysis of 5 loci of chicken FGF6 gene |
由表 4可知，鸡FGF6基因5个处于哈代-温伯格平衡的SNPs位点共存在28种单倍型组合，其中以GCTTCCTCCC组合的频率最高(0.134 9)，GGTTCCTCCT组合的频率(0.126 6)次之，GCTGTCTTCC、GCTTCCCCCC、GCTTCCTCCT、GG-GGTTTCCC、GGGGTTTTCC、GGTGTCTCCT和GGTTTTTCCC组合的频率最低(0.001 4)。为避免因个体数太少引起的误差，单倍型关联分析时剔除个体数少于30的单倍型组合。
|
|
表 4 单倍型组合及其频率 Table 4 Haplotype combinations and their frequencies |
LD分析结果显示5个位点之间强连锁，且属于Block1中，因此仅对rs73399071位点的不同基因型与F2资源群体相关经济性状进行了关联分析。如表 5~8所示，rs73399071位点与鸡0周龄胫长、2周龄体重、4周龄体重、胫围、胸骨长、体斜长、骨盆宽、6周龄体重、8周龄体重、胸宽、胸骨长、胸角、10周龄体重和12周龄体重在内的17个性状显著相关(P＜0.05)；且大部分CC突变基因型个体的表型均值普遍高于野生基因型GG(表 5)。与其他生长性状关联不显著(数据未显示)。
|
|
表 5 鸡FGF6基因遗传变异位点rs73399071与F2资源群生长性状的关联分析 Table 5 Association analysis of genetic variant locus (rs73399071) of chicken FGF6 gene with growth traits in F2 resource population |
|
|
表 6 鸡FGF6基因遗传变异位点(rs73399071)与F2资源群屠体性状的关联分析 Table 6 Association analysis of genetic variant locus (rs73399071) of chicken FGF6 gene with carcass traits of F2 resource population |
|
|
表 7 鸡FGF6基因遗传变异位点(rs73399071)与F2资源群肉质性状的关联分析 Table 7 Association analysis of genetic variant locus (rs73399071) of chicken FGF6 gene with meat quality traits in F2 resource population |
|
|
表 8 鸡FGF6基因遗传变异位点(rs73399071)与F2资源群血清生化指标的关联分析 Table 8 Association analysis of genetic variant locus (rs73399071) of chicken FGF6 gene with serum biochemical indexs of F2 resource population |
rs73399071位点与鸡胸肌重、腿重、腿肌重、屠体重、半净膛重、全净膛重、肝脏率、腿肌重率、半净膛率和全净膛率10个屠体性状显著相关(P＜0.05)。除腿肌重率外，该位点与其他9个表型性状极显著相关(P＜0.01)。且突变基因型CC型个体表型均值普遍显著高于GG基因型(表 6)。与其他屠体性状关联不显著(数据未显示)。
2.11 rs73399071位点与F2资源群肉质性状的关联分析rs73399071位点与肉质性状的关联分析结果显示，如表 7所示，该位点与腿肌剪切力、胸肌剪切力、胸肌失水率和腿肌失水率的关联均未达到显著水平(P＞0.05)。
2.12 rs73399071位点与F2资源群血清生化指标的关联分析rs73399071位点与F2资源群血清生化指标的关联结果分析如表 8所示，该位点与胆碱酯酶和高密度脂蛋白显著相关(P＜0.05)，其中GG基因型和CC基因型普遍高于GC基因型。与其他血清生化指标关联不显著(数据未显示)。
2.13 不同单倍型组合与鸡生长性状的关联分析不同单倍型组合与生长性状关联分析结果如表 9和表 10所示，CCTTCCTTCC单倍型组合个体的4周体重、胸骨长，6周体重，8周体重、胫围，10周体重和12周体重、胸宽均最高，4周胫围、8周胸宽和胸骨长以GGTTCCCCCC组合最高，4周体斜长以GCTTCCTTCT组合最高；与其他体重性状关联不显著(数据未显示)。
|
|
表 9 SNPs位点单倍型组合与鸡不同阶段体重的关联分析 Table 9 Association analysis of haplotype combinations at SNPs loci and chicken body weight at different stages |
|
|
表 10 SNPs位点单倍型组合与鸡不同阶段体尺的关联分析 Table 10 Association analysis of haplotype combinations at SNPs loci and chicken body size at different stages |
本研究对不同物种FGF6基因编码的氨基酸序列进行同源性比对，结果显示，鸡与火鸡的同源性达91%以上，与非洲爪蟾同源性较低，表明FGF6基因在鸟类中的保守性较高。同源性在一定程度上反映亲缘关系的远近，通过构建系统进化树显示鸡FGF6与火鸡亲缘关系较近，与爬行类(绿海龟)和两栖类(热带爪蟾)次之，与鱼类(斑马鱼)最远。由于鸡FGF6基因能够与其他脊椎动物聚集在一起，说明这些物种的FGF6基因可能由同一祖先进化而来。同时经过序列比对发现, 鸡FGF6与其他6个物种存在相同的5个motif基序，推测其可能是在鸡FGF6执行主要生理功能的蛋白片段。通过分析其疏水性和跨膜结构得知鸡FGF6蛋白表现为亲水性蛋白，存在跨膜区。鸡FGF6蛋白的高级结构(二级、三级结构)主要由无规则卷曲和α-螺旋组成。
蛋白质磷酸化是一种广泛的蛋白质翻译后修饰(post-translational modifications，PTMs)方式，也是一种最有效的调控途径，是由蛋白质激酶催化的磷酸基转移反应，在调控细胞周期、信号传导和影响肉品质变化等生物过程中发挥重要作用[19-21]。鸡FGF6蛋白存在12个丝氨酸、5个苏氨酸和5个酪氨酸磷酸化位点，这些磷酸化位点对于FGF6的转运、接收和表达可能具有重要意义。鸡FGF6蛋白有1个N-糖基化位点。糖基化也是一种重要的蛋白质翻译后修饰方式，能够通过调控蛋白质的电荷态、结构及分子间相互作用影响其功能[22]，参与细胞免疫、信号传导、细胞壁的合成、蛋白质的翻译调控和蛋白降解等生物过程[23-24]。通过蛋白互作分析发现，鸡FGF6蛋白与FGFR1、FGFR2、FGFR3、FGFR4等成纤维细胞生长因子受体家族蛋白均有互作关系，有研究表明FGF成员的生物学活性是通过与FGFR1~4受体结合，启动MAPK(有丝分裂原激活的蛋白激酶)信号通路，从而刺激纤维细胞的生长、分裂增殖[25]，这也说明鸡FGF6蛋白与其受体存在互作，具有一定的研究价值。
相关研究表明，FGF6基因会对骨骼肌发育、修复再生等方面产生一定影响[26-28]。FGF6基因多次被鉴定为影响经济性状的候选基因[29-31]，但目前国内外在鸡FGF6基因多态性上未有相关研究报道。杂合度(He)、有效等位基因数(Ne)和多态信息量(PIC)常用来衡量一个群体的遗传变异程度, 其数值越大，表明试验群体遗传变异度越高，遗传丰富度越高[32]。本研究对鸡FGF6基因遗传特性分析发现，6个SNPs位点均处于高等多态性、5个SNPs处于Hardy-Weinberg平衡。可能是本试验中所用群体处于闭锁的生长繁育环境，导致了这种动态平衡。单位点分子标记通常适用于大规模群体的遗传分析，如种质资源的遗传多样性评估、品种鉴定和遗传改良中的辅助选择[33-34]。本研究通过对鸡FGF6基因多态性与经济性状进行关联分析，发现rs73399071位点(G＞C)与F2资源群不同阶段的体重、胫长、胫围、体斜长等17个生长性状表型指标显著相关(P＜0.05)，与全净膛重、半净膛重、屠体重等10个屠体性状表型指标显著相关(P＜0.05)，但与肉质性状的关联未达到显著水平(P＞0.05)，总体来看，CC基因型(突变型)普遍优于GG基因型(野生型)，故在生产实践上rs73399071位点可以作为鸡生长性状和屠体性状的候选分子标记，选留CC基因型个体用于辅助育种，提高经济效益。
在之前的研究中，一般认为基因突变发生在内含子区域不会直接改变蛋白质的氨基酸序列，故不会对其生物学特性产生影响[35]。但随着研究的深入，人们发现内含子可能潜在地影响剪切、mRNA加工、转录和表达调控等各种细胞过程，内含子的变异可能导致剪切准确性或效率的变化，从而导致氨基酸编码的变化，最终改变真核生物的表型[36-38]，如紧密蛋白2(TJP2)基因内含子突变(9q21.11位点)通过改变TJP2基因的表达和促进经表皮水分流失而赋予特应性皮炎(AD)易感性[39]。本研究在FGF6基因上共筛选到6个位于内含子区域的SNPs(rs73399071、rs73399075、rs73399124、rs73399207、rs73399282、rs73399353)，其中rs73399071突变位点与经济性状相关表型关联显著，该位点是否对FGF6基因的表达产生影响，以及如何调控基因来影响动物表型后续还需进一步验证。
动物的表型受单个突变位点和多个突变位点组合效应的影响，其后者影响较为重大[40]。因为在研究遗传变异对表型的影响时，单个SNP分析作用较小，单倍型分析可以提供更丰富的信息[41-43]。单倍型组合分析通常适用于复杂性状的遗传研究，例如复杂疾病的关联分析和数量性状的QTL定位[44]。本研究发现, FGF6基因的CCTTCCTTCC单倍型组合个体的BW4、BW6和BW10均高于其他9种单倍型组合；GGTTCCCCCC单倍型组合个体的SG4、CW8和BBL8高于其他9种单倍型组合，而其SG8和CW12低于CCTTCCTTCC单倍型组合；GCTTCCTTCT单倍型组合的BBL4和BSL4高于其他组合；GGTGTCTTCC单倍型组合的BA12高于其他单倍型。不同单倍型组合对不同阶段体重和体尺指标的关联效果不同。考虑到该研究的样本量，有些单倍型个体数较少，产生的试验结果具有一定的误差，后续还需要进一步扩大样本量来验证和探讨。
4 结论本研究发现，鸡FGF6氨基酸序列与火鸡亲缘关系最近，与斑马鱼亲缘关系最远；不同物种间保守性较高且其与FGFR1、FGFR2、FGFR3、FGFR4、EGF共5个蛋白可能存在相互作用。鸡FGF6基因内含子6个遗传变异位点均表现出高度多态性(PIC≥0.5)，且存在强连锁，GTCCC是主要单倍型；FGF6基因遗传变异位点rs73399071与固始鸡-安卡鸡F2资源群的生长性状、屠体性状和血清生化指标显著相关，且CCTTCCTTCC单倍型组合个体具有较高的体重。相关结果显示，rs73399071位点可以作为改良鸡生长性状的分子标记用于辅助选择，表明FGF6基因可能在鸡生长发育中发挥着重要的作用。本研究为进一步深入研究FGF6基因在鸡上的生物学功能提供重要的参考。
| [1] |
ARMAND A S, LAZIZ I, CHANOINE C. FGF6 in myogenesis[J]. Biochim Biophys Acta, 2006, 1763(8): 773-778. DOI:10.1016/j.bbamcr.2006.06.005 |
| [2] |
郑婕, 姜潮, 李校堃, 等. 成纤维细胞生长因子6(FGF6)的研究进展[J]. 中国生物工程杂志, 2017, 37(4): 110-114. ZHENG J, JIANG C, LI X K, et al. The progression of fibroblast growth factor 6[J]. China Biotechnology, 2017, 37(4): 110-114. (in Chinese) |
| [3] |
李乐乐. 成纤维细胞生长因子与脉络膜新生血管关系的研究进展[J]. 眼科学报, 2020, 35(2): 99-105. LI L L. Research progress on the relationship between fibroblast growth factor and choroidal neovascularization[J]. Eye Science, 2020, 35(2): 99-105. (in Chinese) |
| [4] |
ZHI Y L, CAI C K, XU T X, et al. Silencing of FGF6 hampers aerobic glycolysis and angiogenesis in bladder cancer by regulating PI3K/Akt and MAPK signaling pathways[J]. J Biochem Mol Toxicol, 2023, 37(8): e23399. DOI:10.1002/jbt.23399 |
| [5] |
HAN D, ZHAO H, PARADA C, et al. A TGFβ-Smad4-Fgf6 signaling cascade controls myogenic differentiation and myoblast fusion during tongue development[J]. Development, 2012, 139(9): 1640-1650. DOI:10.1242/dev.076653 |
| [6] |
CAI Q C, WU G B, ZHU M, et al. FGF6 enhances muscle regeneration after nerve injury by relying on ERK1/2 mechanism[J]. Life Sci, 2020, 248: 117465. DOI:10.1016/j.lfs.2020.117465 |
| [7] |
XU B, LIU C Z, ZHANG H, et al. Skeletal muscle-targeted delivery of Fgf6 protects mice from diet-induced obesity and insulin resistance[J]. JCI Insight, 2021, 6(19): e149969. DOI:10.1172/jci.insight.149969 |
| [8] |
SHAMSI F, XUE R D, HUANG T L, et al. FGF6 and FGF9 regulate UCP1 expression independent of brown adipogenesis[J]. Nat Commun, 2020, 11(1): 1421. DOI:10.1038/s41467-020-15055-9 |
| [9] |
LIU C Z, MENG M Y, XU B, et al. Fibroblast growth factor 6 promotes adipocyte progenitor cell proliferation for adipose tissue homeostasis[J]. Diabetes, 2023, 72(4): 467-482. DOI:10.2337/db22-0585 |
| [10] |
HU Z C, CHEN P, WANG L L, et al. FGF6 promotes cardiac repair after myocardial infarction by inhibiting the Hippo pathway[J]. Cell Prolif, 2022, 55(5): e13221. DOI:10.1111/cpr.13221 |
| [11] |
CREIGHTON C J, ZHANG F, ZHANG Y Q, et al. Comparative and integrative analysis of transcriptomic and epigenomic-wide DNA methylation changes in African American prostate cancer[J]. Epigenetics, 2023, 18(1): 2180585. DOI:10.1080/15592294.2023.2180585 |
| [12] |
KUMAR M, CHAPMAN S C. Cloning and expression analysis of Fgf5, 6 and 7 during early chick development[J]. Gene Expr Patterns, 2012, 12(7/8): 245-253. |
| [13] |
曹玉珠, 贾其辉, 邢雨欣, 等. 鸡成纤维细胞生长因子家族生物信息学特性及表达调控[J]. 河南农业大学学报, 2023, 57(3): 422-435. CAO Y Z, JIA Q H, XING Y X, et al. Bioinformatics characteristics and expression regulation of chicken fibroblast growth factor family[J]. Journal of Henan Agricultural University, 2023, 57(3): 422-435. (in Chinese) |
| [14] |
付金秀, 郭红威, 李超, 等. 鸡DGAT2基因第7外显子多态性与生产性能的关联[J]. 农业生物技术学报, 2016, 24(5): 689-696. FU J X, GUO H W, LI C, et al. Polymorphism of exon 7 of DGAT2 gene and its association with growth traits in chicken(Gallus gallus)[J]. Journal of Agricultural Biotechnology, 2016, 24(5): 689-696. (in Chinese) |
| [15] |
ZHANG Y H, WANG Y Z, LI Y Y, et al. Genome-wide association study reveals the genetic determinism of growth traits in a Gushi-Anka F2 chicken population[J]. Heredity (Edinb), 2021, 126(2): 293-307. DOI:10.1038/s41437-020-00365-x |
| [16] |
苏玲. 鸡PNPLA3基因多态性及组织表达规律的研究[D]. 郑州: 河南农业大学, 2012. SU L. Study of polymorphisms and gene expression profile of chicken PNPLA3 gene[D]. Zhengzhou: Henan Agricultural University, 2012. (in Chinese) |
| [17] |
SONG H J, LI W T, LI Y F, et al. Genome-wide association study of 17 serum biochemical indicators in a chicken F2 resource population[J]. BMC Genomics, 2023, 24(1): 98. DOI:10.1186/s12864-023-09206-7 |
| [18] |
梅承君. 鸡血清生化指标与部分脂肪性状相关性及谷草转氨酶同工酶的研究[D]. 郑州: 河南农业大学, 2007. MEI C J. Studies on the corralation between serum biochemical indexes and some fat traits and got isozymes in chickens[D]. Zhengzhou: Henan Agricultural University, 2007. (in Chinese) |
| [19] |
梁前进, 王鹏程, 白燕荣. 蛋白质磷酸化修饰研究进展[J]. 科技导报, 2012, 30(31): 73-79. LIANG Q J, WANG P C, BAI Y R. Summarization on the progress in protein phosphorylation[J]. Science & Technology Review, 2012, 30(31): 73-79. (in Chinese) |
| [20] |
谭锬, 柯柏怡, 梁前进. 蛋白质磷酸化修饰及其在细胞周期调控中的作用研究进展[J/OL]. 北京师范大学学报: 自然科学版, 1-8. http://kns.cnki.net/kcms/detail/11.1991.n.20230720.1204.002.html. TAN T, KE B Y, LIANG Q J. Research progress on protein phosphorylation modification and its role in regulating cell cycle progression[J]. Journal of Beijing Normal University: Natural Science, 1-8. http://kns.cnki.net/kcms/detail/11.1991.n.20230720.1204.002.html. (in Chinese) |
| [21] |
王文琪, 张雅玮, 李加慧, 等. 蛋白质磷酸化对宰后肉品质影响研究进展[J]. 食品科学, 2023, 44(9): 221-230. WANG W Q, ZHANG Y W, LI J H, et al. Research progress on the effect of protein phosphorylation on postmortem meat quality[J]. Food Science, 2023, 44(9): 221-230. (in Chinese) |
| [22] |
贾翼菲, 王雅梅, 李功玉. 蛋白质糖基化修饰的非变性构象分辨质谱研究进展[J/OL]. 中国药科大学学报, 1-13. http://kns.cnki.net/kcms/detail/32.1157.r.20230921.1134.004.html. JIA Y F, WANG Y M, LI G Y. Recent progress of protein glycosylation characterization utilizing native conformer-resolved mass spectrometry[J/OL]. Journal of China Pharmaceutical University, 1-13. http://kns.cnki.net/kcms/detail/32.1157.r.20230921.1134.004.html. (in Chinese) |
| [23] |
李军, 杜鑫, HOSSEINI M S H, 等. 蛋白质糖基化修饰研究进展[J]. 科技通报, 2009, 25(6): 773-778, 783. LI J, DU X, HOSSEINI M S H, et al. The research progress in protein glycosylation modification[J]. Bulletin of Science and Technology, 2009, 25(6): 773-778, 783. (in Chinese) |
| [24] |
王家红, 童玥, 朱玥, 等. 蛋白质糖基化的研究进展[J]. 药物生物技术, 2011, 18(1): 77-80. WANG J H, TONG Y, ZHU Y, et al. The research progress in protein glycosylation[J]. Pharmaceutical Biotechnology, 2011, 18(1): 77-80. (in Chinese) |
| [25] |
李伟. 酸性成纤维细胞因子的研究进展[J]. 科技信息, 2010(22): 493, 496. LI W. Advances in the study of acidic fibroblast factors[J]. Science & Technology Information, 2010(22): 493, 496. (in Chinese) |
| [26] |
张力青, 付平, 易四凤, 等. 成纤维细胞生长因子(FGF)在骨骼肌发育中的作用[J]. 养殖与饲料, 2016(1): 6-10. ZHANG L Q, FU P, YI S F, et al. Role of fibroblast growth factor (FGF) in skeletal muscle development[J]. Animals Breeding and Feed, 2016(1): 6-10. (in Chinese) |
| [27] |
ARMAND A S, LAUNAY T, PARISET C, et al. Injection of FGF6 accelerates regeneration of the soleus muscle in adult mice[J]. Biochim Biophys Acta, 2003, 1642(1/2): 97-105. |
| [28] |
HAN J K, MARTIN G R. Embryonic expression of Fgf-6 is restricted to the skeletal muscle lineage[J]. Dev Biol, 1993, 158(2): 549-554. DOI:10.1006/dbio.1993.1212 |
| [29] |
YIN T, KÖNIG S. Genome-wide associations and detection of potential candidate genes for direct genetic and maternal genetic effects influencing dairy cattle body weight at different ages[J]. Genet Sel Evol, 2019, 51(1): 4. DOI:10.1186/s12711-018-0444-4 |
| [30] |
GHOREISHIFAR S M, ERIKSSON S, JOHANSSON A M, et al. Signatures of selection reveal candidate genes involved in economic traits and cold acclimation in five Swedish cattle breeds[J]. Genet Sel Evol, 2020, 52(1): 52. DOI:10.1186/s12711-020-00571-5 |
| [31] |
WEI Y L, GUO D S, BAI Y B, et al. Transcriptome analysis of mRNA and lncRNA related to muscle growth and development in Gannan yak and jeryak[J]. Int J Mol Sci, 2023, 24(23): 16991. DOI:10.3390/ijms242316991 |
| [32] |
李青, 王慧华, 卢曾奎, 等. 绵羊PDGF-D基因多态性与尾型的关联[J]. 畜牧兽医学报, 2019, 50(4): 688-700. LI Q, WANG H H, LU Z K, et al. Association of the PDGF-D gene genetic variation with sheep tail traits[J]. Acta Veterinaria et Zootechnica Sinica, 2019, 50(4): 688-700. (in Chinese) |
| [33] |
武师良, 王昱彤, 党婉怡, 等. 基于SSR分子标记的辽西地区驴群体遗传多样性分析[J]. 中国畜禽种业, 2023, 19(5): 27-32. WU S L, WANG Y T, DANG W Y, et al. Genetic diversity analysis of donkey populations in western Liaoning based on SSR molecular markers[J]. The Chinese Livestock and Poultry Breeding, 2023, 19(5): 27-32. (in Chinese) |
| [34] |
王培国, 张君权, 贾鸿昌, 等. 大豆PRR家族基因KASP标记的开发[J]. 大豆科学, 2023, 42(4): 385-395. WANG P G, ZHANG J Q, JIA H C, et al. Development of KASP markers for soybean PRR family genes[J]. Soybean Science, 2023, 42(4): 385-395. (in Chinese) |
| [35] |
ZHANG Z C, LIU C, HAO W J, et al. Novel single nucleotide polymorphisms and haplotype of MYF5 gene are associated with body measurements and ultrasound traits in grassland short-tailed sheep[J]. Genes (Basel), 2022, 13(3): 483. DOI:10.3390/genes13030483 |
| [36] |
PENG X W, CHENG J L, LI H, et al. Whole-genome sequencing reveals adaptations of hairy-footed jerboas (Dipus, Dipodidae) to diverse desert environments[J]. BMC Biol, 2023, 21(1): 182. DOI:10.1186/s12915-023-01680-5 |
| [37] |
彭雅鑫, 刘军, 赵诗瑜, 等. 猪RXRB基因SNPs检测及其与生长育肥和繁殖性状的关联分析[J]. 畜牧兽医学报, 2021, 52(3): 596-609. PENG Y X, LIU J, ZHAO S Y, et al. Detection of SNPs in porcine RXRB gene and their association analysis with growth, fattening and reproduction traits[J]. Acta Veterinaria et Zootechnica Sinica, 2021, 52(3): 596-609. (in Chinese) |
| [38] |
JO B S, CHOI S S. Introns: the functional benefits of introns in genomes[J]. Genomics Inform, 2015, 13(4): 112-118. DOI:10.5808/GI.2015.13.4.112 |
| [39] |
ELIZA L Y Y, YIE S Y, SAY Y H, et al. Deep intronic 9q21.11 polymorphism contributes to atopic dermatitis risk through methylation regulated expression of tight junction protein 2[J]. J Investig Allergol Clin Immunol, 2023. DOI:10.18176/jiaci.0978 |
| [40] |
OROZCO G, HINKS A, EYRE S, et al. Combined effects of three independent SNPs greatly increase the risk estimate for RA at 6q23[J]. Hum Mol Genet, 2009, 18(14): 2693-2699. DOI:10.1093/hmg/ddp193 |
| [41] |
ZHAO P F, HE Z H, XI Q M, et al. Variations in HIF-1α contributed to high altitude hypoxia adaptation via affected oxygen metabolism in tibetan sheep[J]. Animals (Basel), 2021, 12(1): 58. |
| [42] |
DA Y. Multi-allelic haplotype model based on genetic partition for genomic prediction and variance component estimation using SNP markers[J]. BMC Genet, 2015, 16: 144. |
| [43] |
YUAN Y M, JIANG H O, KUANG J Y, et al. Genetic variations in ADIPOQ gene are associated with chronic obstructive pulmonary disease[J]. PLoS One, 2012, 7(11): e50848. DOI:10.1371/journal.pone.0050848 |
| [44] |
胡艳玲. 复杂性状与基因组多位点的关联分析方法研究[D]. 上海: 上海交通大学, 2009. HU Y L. Study on the methods of genome multi-loci association analysis for complex trait[D]. Shanghai: Shanghai Jiao Tong University, 2009. (in Chinese) |
(编辑 郭云雁)