中国农业科学 ›› 2024, Vol. 57 ›› Issue (1): 204-215.doi: 10.3864/j.issn.0578-1752.2024.01.014
姬改革1(), 陈智武2(
), 单艳菊1, 刘一帆1, 屠云洁1, 邹剑敏1, 章明1, 巨晓军1, 束婧婷1(
), 张海涛3, 唐燕飞4, 蒋华莲4
收稿日期:
2023-05-09
接受日期:
2023-10-31
出版日期:
2024-01-01
发布日期:
2024-01-10
通信作者:
联系方式:
姬改革,E-mail:jigaige@126.com。陈智武,E-mail:ubtech@163.com。姬改革与陈智武为同等贡献者。
基金资助:
JI GaiGe1(), CHEN ZhiWu2(
), SHAN YanJu1, LIU YiFan1, TU YunJie1, ZOU JianMin1, ZHANG Ming1, JU XiaoJun1, SHU JingTing1(
), ZHANG HaiTao3, TANG YanFei4, JIANG HuaLian4
Received:
2023-05-09
Accepted:
2023-10-31
Published:
2024-01-01
Online:
2024-01-10
摘要:
【目的】通过比较干毛和未干毛羽毛毛囊的形态学和基因表达差异,挖掘调控黄羽肉鸡羽毛干毛性状的重要候选基因和信号通路。【方法】分别采集背部未干毛和干毛羽毛皮肤组织样本各3个,运用组织切片技术,比较干毛和未干毛羽毛毛囊形态学差异;运用RNA-seq技术,比较两组样本之间的差异表达基因(differentially expressed genes,DEGs),并对其进行功能富集分析和蛋白互作网络构建;采用荧光定量PCR技术(RT-qPCR)对测序结果准确性进行验证。【结果】证实未干毛羽毛的毛囊处于生长期,已干毛的羽毛毛囊处于静止期。以未干毛毛囊(生长期)皮肤样本为对照,在干毛毛囊(静止期)皮肤样本中发现了942个DEGs(|fold-change|>2和P<0.05),其中384个基因表达下调,558个基因表达上调。Go功能分析显示细胞分裂、周期调控等相关生物过程被显著富集(P<0.05)。KEGG分析发现,MAPK、TGF-β、p53及DNA复制等相关信号通路被显著富集(P<0.05)。构建差异蛋白相互作用(PPI)网络,通过CytoHubba分析获得6个hub基因,分别为CDK1、MAD2L1、BUB1、CCNB2、PLK1和BUB1B。GSEA富集分析筛选到紧密连接、胰岛素、MAPK、TGF-β和细胞周期等信号通路与鸡羽毛毛囊生长周期显著关联(|NES|>1, FDR<0.25)。RT-qPCR结果显示8个DEGs表达趋势与RNA-seq结果基本一致。【结论】鸡羽毛的干毛性状与毛囊周期发育相关,MAPK和TGF-β等信号通路可能通过调控细胞周期相关基因的表达在羽毛生长发育中发挥重要作用;结果为进一步深入了解黄羽肉鸡羽毛干毛性状分子调控机制提供了基础资料。
姬改革, 陈智武, 单艳菊, 刘一帆, 屠云洁, 邹剑敏, 章明, 巨晓军, 束婧婷, 张海涛, 唐燕飞, 蒋华莲. 基于转录组测序研究调控黄羽肉鸡干毛性状关键基因和信号通路[J]. 中国农业科学, 2024, 57(1): 204-215.
JI GaiGe, CHEN ZhiWu, SHAN YanJu, LIU YiFan, TU YunJie, ZOU JianMin, ZHANG Ming, JU XiaoJun, SHU JingTing, ZHANG HaiTao, TANG YanFei, JIANG HuaLian. Study of Key Genes and Signaling Pathways Regulating Dry Feather Traits in Yellow-Feathered Broiler Chickens Based on Transcriptome Analysis[J]. Scientia Agricultura Sinica, 2024, 57(1): 204-215.
表1
引物序列表"
基因 Gene symbol | 序列号 Accession number | 引物序列 Primer sequence (5′→3′) | 产物长度 Product length (bp) |
---|---|---|---|
β-actin | NM_205518.2 | F: GCTGAACTCCATCTTGGTTCCC R: TGCTTCTAGGGTCACTCGCA | 150 |
BUB1 | NM_001012870.1 | F: GCTGAACTCCATCTTGGTTCCC R: TGCTTCTAGGGTCACTCGCA | 112 |
BIRC5 | NM_001012319.1 | F:CTGGATAAAAAGCGGACCAA R: AACCTAAGGGCCCATGTTCT | 246 |
CDK1 | NM_205314.1 | F:GGAACGCATGTCCAAAACTT R: GCAGGCAGGCAAAGATAAAG | 236 |
NEK2 | NM_001031050.2 | F:ATGTCTTCCTGGATGGCAAG R: TCTGCCAACTCCTTCTGGTT | 241 |
CDC6 | XM_040691823.1 | F:GTGTCCTCTCGGAGGTGTTT R: AGGCACTCTGTCTGGTCCAC | 214 |
ANXA1 | NM_206906.1 | F:AAAAACTCCACCTGGCAATG R: CCACATAAAGCGACCAGGAT | 193 |
BCL2 | NM_205339.2 | F: TGACACCTCGATCTCACAGG R: CATGCAACACACTGGGATTC | 247 |
EDA2R | NM_001083360.1 | F: ACCCTCATCAACCGCATCCA R: GCCTTGCGGTAAAACCCTGG | 89 |
FZD4 | NM_204099.1 | F:AATGTCACCAAGATGCCCAACC R:AGACCGAACAAAGGAAGAACTGGA | 130 |
表2
测序信息统计及与参考基因组比对"
样品 Sample | 原始序列 Raw reads | 有效序列 Clean reads | Q30 (%) | GC (%) | 比对率 Mapped on reference | 未比对率 Unmapped | 唯一位置比对 Uniquely mapped |
---|---|---|---|---|---|---|---|
H1-1 | 40968048 | 38142462 | 92.33 | 48.605 | 36255024 (95.05%) | 1887438 (4.95%) | 35196169 (92.28%) |
H1-2 | 46296028 | 43781388 | 93.07 | 49.900 | 41547151 (94.9%) | 2234237 (5.1%) | 40274977 (91.99%) |
H1-3 | 48737026 | 45301826 | 92.01 | 48.995 | 42668140 (94.19%) | 2633686 (5.81%) | 41786822 (92.24%) |
H2-1 | 46250350 | 43867410 | 93.33 | 49.985 | 41559772 (94.74%) | 2307638 (5.26%) | 40597907 (92.55%) |
H2-2 | 43917200 | 40983494 | 92.28 | 49.945 | 38739383 (94.52%) | 2244111 (5.48%) | 37804671 (92.24%) |
H2-3 | 42372910 | 39985456 | 93.17 | 49.750 | 38040749 (95.14%) | 1944707 (4.86%) | 36991013 (92.51%) |
表3
通过CytoHubba插件4种算法在PPI网络中排名前10的基因"
基因 Gene | MCC | 基因 Gene | Degree | 基因 Gene | MNC | 基因 Gene | EPC | |||
---|---|---|---|---|---|---|---|---|---|---|
CDK1 | 5.11E+19 | CDK1 | 56 | CDK1 | 55 | APITD1 | 85.564 | |||
BUB1 | 5.11E+19 | PLK1 | 46 | PLK1 | 46 | AURKA | 85.564 | |||
PLK1 | 5.11E+19 | CCNB2 | 44 | CCNB2 | 43 | CDK1 | 85.564 | |||
CCNB2 | 5.11E+19 | CENPE | 43 | CENPE | 43 | TPX2 | 85.564 | |||
NDC80 | 5.11E+19 | BUB1 | 41 | BUB1 | 41 | BIRC5 | 85.564 | |||
CENPE | 5.11E+19 | NDC80 | 41 | NDC80 | 41 | MAD2L1 | 85.564 | |||
BUB1B | 5.11E+19 | BUB1B | 37 | BUB1B | 37 | BUB1 | 85.564 | |||
MAD2L1 | 5.11E+19 | MAD2L1 | 34 | KIF11 | 34 | CCNB2 | 85.564 | |||
CENPF | 5.11E+19 | KIF11 | 34 | MAD2L1 | 33 | PLK1 | 85.564 | |||
INCENP | 5.11E+19 | AURKA | 32 | AURKA | 31 | BUB1B | 85.564 |
[1] |
辛翔飞, 郑麦青, 文杰, 王济民. 2021年我国肉鸡产业形势分析、未来展望与对策建议. 中国畜牧杂志, 2022, 58(3): 222-226.
|
|
|
[2] |
张德祥, 黄军. 土鸡羽毛成熟度的选育效果分析//中国家禽科学研究进展——第十四次全国家禽科学学术讨论会论文集. 哈尔滨, 2009: 339-342.
|
|
|
[3] |
李培周. 鸡羽毛不同成熟度初级毛囊组织学观察及其与VEGF、VEGFR-2基因的相关性研究[D]. 湛江: 广东海洋大学, 2013.
|
|
|
[4] |
doi: 10.1007/s12565-018-0449-7 |
[5] |
doi: 10.1016/j.jmb.2015.07.009 |
[6] |
doi: 10.1016/j.ceb.2006.10.009 |
[7] |
doi: 10.1016/j.ydbio.2014.01.010 pmid: 24463139 |
[8] |
|
[9] |
doi: 10.1186/1471-2164-16-1 |
[10] |
doi: 10.1146/animal.2015.3.issue-1 |
[11] |
doi: 10.3389/fphys.2019.01609 |
[12] |
doi: 10.3390/genes13071249 |
[13] |
doi: S0012-1606(16)30600-5 pmid: 27847323 |
[14] |
doi: 10.1111/jbg.v138.1 |
[15] |
doi: 10.1093/bioinformatics/btp120 pmid: 19289445 |
[16] |
|
[17] |
doi: 10.1089/omi.2011.0118 pmid: 22455463 |
[18] |
|
[19] |
|
[20] |
doi: 10.1111/exd.14070 pmid: 31903650 |
[21] |
doi: 10.1080/00071668.2018.1564244 pmid: 30602288 |
[22] |
doi: 10.1002/jcp.v234.1 |
[23] |
doi: 10.1007/s10549-020-05962-2 |
[24] |
doi: 10.1016/j.bcp.2021.114747 |
[25] |
doi: 10.1016/j.celrep.2021.109740 |
[26] |
|
[27] |
doi: 10.18632/aging.v9i12 |
[28] |
|
[29] |
doi: 10.1038/cdd.2011.197 pmid: 22223105 |
[30] |
doi: 10.7150/jca.11404 |
[31] |
doi: 10.1016/j.tranon.2019.09.011 |
[32] |
doi: 10.1074/jbc.M110.106575 |
[33] |
陈兴勇, 谢珊珊, 周丽, 姜润深, 耿照玉. 皖西白鹅换羽前后基因表达谱差异分析. 畜牧兽医学报, 2013, 44(7): 1030-1036.
|
|
|
[34] |
doi: 10.3390/ijms15011647 pmid: 24451143 |
[35] |
doi: 10.1093/jmcb/mjy051 |
[36] |
doi: S0888-7543(18)30576-7 pmid: 30779940 |
[37] |
屠云洁, 姬改革, 章明, 刘一帆, 巨晓军, 单艳菊, 邹剑敏, 李华, 陈智武, 束婧婷. 鸡Wnt3a的SNPs筛选及其与皮肤毛囊密度性状关联分析. 中国农业科学, 2022, 55(23): 4769-4780.
doi: 10.3864/j.issn.0578-1752.2022.23.016 |
|
|
[38] |
|
[39] |
doi: 10.1186/s12864-022-09060-z pmid: 36510127 |
[40] |
doi: 10.1186/s12864-020-6725-7 pmid: 32393236 |
[41] |
doi: 10.1186/s12864-018-4773-z pmid: 29801437 |
[42] |
doi: 10.1186/s12861-017-0153-1 pmid: 28728543 |
[43] |
doi: S2211-1247(15)01026-8 pmid: 26456821 |
[44] |
doi: 10.5194/aab-63-155-2020 pmid: 32490151 |
[45] |
冯云奎, 王健, 马金亮, 张柳明, 李拥军. miR-31-5p对山羊毛囊干细胞增殖和凋亡的影响. 中国农业科学, 2021, 54(23): 5132-5143.
doi: 10.3864/j.issn.0578-1752.2021.23.017 |
doi: 10.3864/j.issn.0578-1752.2021.23.017 |
|
[46] |
doi: 10.3390/cells11142216 |
[1] | 张颖, 石婷瑞, 曹瑞, 潘文秋, 宋卫宁, 王利, 聂小军. ICARDA引进-小麦苗期抗旱性的全基因组关联分析[J]. 中国农业科学, 2024, 57(9): 1658-1673. |
[2] | 戚仁洁, 宁宇, 刘静, 刘之洋, 徐海, 罗志丹, 陈龙正. 基于转录组测序的苦瓜皂苷合成相关基因的鉴定和分析[J]. 中国农业科学, 2024, 57(9): 1779-1793. |
[3] | 黄立强, 江如, 朱波汁, 彭焕, 许翀, 宋家雄, 陈敏, 李永青, 黄文坤, 彭德良. 马铃薯主栽品种抗马铃薯金线虫鉴定及抗性分子标记检测[J]. 中国农业科学, 2024, 57(8): 1506-1516. |
[4] | 许娜, 唐颖, 徐正进, 孙健, 徐铨. 籼粳杂种不育的遗传分析和候选基因鉴定[J]. 中国农业科学, 2024, 57(8): 1417-1429. |
[5] | 陈兵先, 张琪, 戴彰言, 周旭, 刘军. 水杨酸引发提高低温下水稻种子萌发活力的生理与分子效应[J]. 中国农业科学, 2024, 57(7): 1220-1236. |
[6] | 王程泽, 张燕, 付伟, 贾京哲, 董金皋, 申珅, 郝志敏. 玉米ACO基因家族生物信息学及表达模式分析[J]. 中国农业科学, 2024, 57(7): 1308-1318. |
[7] | 袁苗, 周娟, 党仕卓, 汤学燊, 张亚红. ‘红地球’葡萄VvARF18功能分析[J]. 中国农业科学, 2024, 57(7): 1363-1376. |
[8] | 赵真坚, 王凯, 陈栋, 申琦, 余杨, 崔晟頔, 王俊戈, 陈子旸, 禹世欣, 陈佳苗, 王翔枫, 唐国庆. 基因组和DNA甲基化组联合分析筛选猪肉质性状关键基因[J]. 中国农业科学, 2024, 57(7): 1394-1406. |
[9] | 吴期滨, 谢婉婕, 钟惠, 冯春燕, 潘浩然, 齐浥颖, 张积森, 王恒波. 甘蔗抗褐锈病Bru1区域的鉴定及候选抗病基因的功能分析[J]. 中国农业科学, 2024, 57(5): 855-867. |
[10] | 胡衡亮, 顾天宇, 陈思颖, 王垚, 彭佳师. 莲中介导矿质元素胁迫耐受性基因的筛选与功能验证[J]. 中国农业科学, 2024, 57(5): 980-988. |
[11] | 张必东, 林泓, 朱思颖, 李忠成, 庄慧, 李云峰. 水稻颖壳异常突变体ah1的鉴定与候选基因分析[J]. 中国农业科学, 2024, 57(3): 429-441. |
[12] | 熊楚雯, 郭智滨, 周强华, 程艳波, 马启彬, 蔡占东, 年海. 大豆转录因子NAC1耐低磷胁迫的功能研究[J]. 中国农业科学, 2024, 57(3): 442-453. |
[13] | 李博天, 李春琪, 刘国平, 谢军, 曾攀, 赵润泽, 李桐, 裴洁, 郭利伟, 伍锐, 谭磊. 2020-2022年华中地区部分猪场猪圆环病毒3型的分子流行特征与遗传变异分析[J]. 中国农业科学, 2024, 57(3): 613-626. |
[14] | 亓晓蕾, 王君, 吕广德, 牟秋焕, 米勇, 孙盈盈, 尹逊栋, 钱兆国, 王瑞霞, 吴科. 优质高产小麦新品种泰科麦33的遗传构成分析[J]. 中国农业科学, 2024, 57(22): 4391-4401. |
[15] | 吴传磊, 胡晓渝, 王伟, 苗龙, 白鹏宇, 王郭伋, 李娜, 舒阔, 邱丽娟, 王晓波. 大豆油分相关功能基因分子标记开发与鉴定及优异等位变异聚合分析[J]. 中国农业科学, 2024, 57(22): 4402-4415. |
|