基于RAD-seq静原鸡保种群体的遗传变异分析

doi:10.11843/j.issn.0366-6964.2022.07.008

摘要/Abstract

摘要： 旨在分析静原鸡白羽、麻羽、黑羽保种群之间的遗传变异，挖掘调控特征性状的关键候选基因。本研究利用RAD-seq技术对180日龄特征明显、发育正常且健康的白羽、麻羽、黑羽静原鸡(每个羽色选取60个个体，40只母鸡，20只公鸡)进行测序，基于SNP标记计算观察杂合度(Ho)、群体内核酸多态性(Pi)、群体的平均近交系数(Fis)等指标，分析静原鸡3个类群的遗传多样性和群体结构，并通过选择性清除分析和全基因组关联分析(GWAS)筛选出候选基因，利用KOBAS对候选基因进行KEGG (Kyoto Encyclopedia of Genes and Genomes)通路富集,最终筛选出调控静原鸡羽色的候选基因。测序结果表明，静原鸡180个样品共产生198.83 Gb Clean Data，Q30达到93%以上。遗传多样性分析表明，白羽、麻羽、黑羽静原鸡分别鉴定出238 533、233 562和240 820个SNPs标记，Ho、Pi和Fis分别在0.273 2~0.278 2、0.304 9~0.309 6和0.096 1~0.109 8之间。群体结构分析表明，静原鸡根据不同羽色分为不同类群。通过选择性清除分析和全基因组关联分析共筛选出11个(FZD4、WNT16、EDNRB、TYR、KRAS、CTNNB1、DDC、MC1R、CAMK2A、PRKCB、PRKCA)与静原鸡羽色相关的候选基因，富集结果显示这些基因主要与黑色素生成、酪氨酸代谢和Wnt信号传导等通路相关。综上所述，本研究利用SNPs标记信息可以全面地评价静原鸡的保种现状，为静原鸡的遗传资源保护奠定理论基础。同时，筛选出了与静原鸡羽色性状相关的候选基因，为静原鸡不同羽色的品系化培育提供新的遗传标记和基因靶点。

关键词: 静原鸡, RAD-seq, 遗传多样性分析, 选择性清除分析, 全基因组关联分析

Abstract: The aim of this study was to analyze the genetic variation among the white-, hemp- and black-feathered breeding populations of Jingyuan chickens, and to explore the key candidate genes regulating the characteristic traits. In this study, the RAD-seq technology was used to sequence 180-day-old white-, hemp- and black-feathered Jingyuan chickens (60 individuals of each plumage color, 40 hens and 20 roosters) with distinctive features, normal development and health, and the observed heterozygosity (Ho), nucleic acid polymorphism (Pi) within the population and the average inbreeding coefficient (Fis) of the population based on SNPs markers were calculated to analyze the genetic diversity and population structure of Jingyuan chickens. The candidate genes were screened by selective sweep analysis and genome-wide association analysis (GWAS), and enriched by KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways using KOBAS. The candidate genes regulating plumage color of Jingyuan chickens were finally screened. The sequencing results showed that a total of 198.83 Gb of clean data was obtained from 180 samples of Jingyuan chickens, with Q30 reaching more than 93%. Genetic diversity analysis showed that 238 533, 233 562 and 240 820 SNPs markers were identified for the white-, hemp-, and black- feathers of Jingyuan chickens, with Ho, Pi and Fis range of 0.273 2-0.278 2, 0.304 9-0.309 6 and 0.096 1-0.109 8, respectively. The analysis of population structure showed that Jingyuan chickens were divided into different groups according to different feather colors. A total of 11 candidate genes (FZD4, WNT16, EDNRB, TYR, KRAS, CTNNB1, DDC, MC1R, CAMK2A, PRKCB, PRKCA) associated with plumage color in Jingyuan chickens were screened by selective sweep analysis and genome-wide association analysis, and the enrichment results showed that these genes were related to the pathways of melanogenesis, tyrosine metabolism, and Wnt signaling. In summary, the breeding status of the Jingyuan chicken can be comprehensively evaluated using SNPs marker information and lay a theoretical foundation for the conservation of the genetic resources of the Jingyuan chicken in this study. At the same time, the candidate genes related to the plumage color trait of Jingyuan chicken were screened to provide new genetic markers and gene targets for the breeding of different plumage colors of Jingyuan chicken.

Key words: Jingyuan chicken, RAD-seq, genetic diversity analysis, selective sweep analysis, genome-wide association analysis

中图分类号:

S831.2

马丽霞, 曹国伟, 朱红芳, 邓占钊, 蔡正云, 周成浩, 韩威, 顾亚玲, 张娟. 基于RAD-seq静原鸡保种群体的遗传变异分析[J]. 畜牧兽医学报, 2022, 53(7): 2104-2117.

MA Lixia, CAO Guowei, ZHU Hongfang, DENG Zhanzhao, CAI Zhengyun, ZHOU Chenghao, HAN Wei, GU Yaling, ZHANG Juan. Analysis of Genetic Variation in a Conserved Population of Jingyuan Chickens Based on RAD-seq[J]. Acta Veterinaria et Zootechnica Sinica, 2022, 53(7): 2104-2117.

参考文献 34

[1]	母童.静原鸡ELOVL2和ELOVL5基因组织表达特性及生物信息学分析[D].银川:宁夏大学,2019.MU T.Tissue-specific expression and bioinformatics analysis of ELOVL2 and ELOVL5 gene strains in Jingyuan chicken[D]. Yinchuan:Ningxia University,2019.(in Chinese)
[2]	额尔和花,丁伟,杨文清,等.宁夏地方品种鸡的开发与利用[J].畜牧与饲料科学,2013,34(7-8):106-107.E E H H,DING W,YANG W Q,et al.Development and utilization of local breed chicken in Ningxia[J].Animal Husbandry and Feed Science,2013,34(7-8):106-107.(in Chinese)
[3]	OLIEHOEK P A,BIJMA P.Effects of pedigree errors on the efficiency of conservation decisions[J].Genet Sel Evol,2009,41(1):9.
[4]	李国辉,魏清宇,张会永,等.基于RAD-seq技术分析边鸡保种群体的遗传进化[J].农业生物技术学报,2019,27(12):2198-2206.LI G H,WEI Q Y,ZHANG H Y,et al.The Analysis on genetic evolution of Bian chicken (Gallus gallus) population based on RAD-seq technology[J].Journal of Agricultural Biotechnology,2019,27(12):2198-2206.(in Chinese)
[5]	韩威,朱云芬,殷建玫,等.基于RAD-seq简化基因组测序的19个地方鸡种遗传进化研究[J].畜牧兽医学报,2020,51(4):670-678.HAN W,ZHU Y F,YIN J M,et al.Study on genetic evolution of 19 indigenous chicken breeds based on RAD-seq[J].Acta Veterinaria et Zootechnica Sinica,2020,51(4):670-678.(in Chinese)
[6]	LIU C M,WANG S H,DONG X G,et al.Exploring the genomic resources and analysing the genetic diversity and population structure of Chinese indigenous rabbit breeds by RAD-seq[J].BMC Genomics,2021,22(1):573.
[7]	LIU C M,CHEN H L,REN Z J,et al.Development of Genomic resources and identification of genetic diversity and genetic structure of the domestic Bactrian Camel in China by RAD sequencing[J].Front Genet,2020,11:797.
[8]	WANG M S,THAKUR M,PENG M S,et al.863 genomes reveal the origin and domestication of chicken[J].Cell Res,2020,30(8):693-701.
[9]	ZHANG M M,HAN W,TANG H,et al.Genomic diversity dynamics in conserved chicken populations are revealed by genome-wide SNPs[J].BMC Genomics,2018,19(1):598.
[10]	史良玉,王立刚,张鹏飞,等.不同来源大白猪总产仔数近交衰退评估[J].畜牧兽医学报,2021,52(10):2772-2782.SHI L Y,WANG L G,ZHANG P F,et al.Evaluation of inbreeding depression on the total numbers of piglets born in different groups of large white pigs[J].Acta Veterinaria et Zootechnica Sinica,2021,52(10):2772-2782.(in Chinese)
[11]	STEPHAN W.Selective sweeps[J].Genetics,2019,211(1):5-13.
[12]	吕世杰,陈付英,金磊,等.利用简化基因组测序筛选安格斯牛生长相关的受选择基因[J].畜牧兽医学报,2020,51(4):713-721.LV S J,CHEN F Y,JIN L,et al.Identification of growth-related genes under selection in angus cattle using SLAF-seq[J].Acta Veterinaria et Zootechnica Sinica,2020,51(4):713-721.(in Chinese)
[13]	YU S,LIAO J,TANG M,et al.A functional single nucleotide polymorphism in the tyrosinase gene promoter affects skin color and transcription activity in the black-boned chicken[J].Poult Sci,2017,96(11):4061-4067.
[14]	李洪林.贵州地方鸡种隐性白羽与慢羽分子标记研究与应用[D].贵州:贵州大学,2019.LI H L.Study and Application of recessive white feather and slow feather molecular markers in Guizhou local chicken breeds[D]. Guizhou:Guizhou University,2019.(in Chinese)
[15]	WU N,QIN H,WANG M,et al.Variations in endothelin receptor B subtype 2(EDNRB2) coding sequences and mRNA expression levels in 4 Muscovy duck plumage colour phenotypes[J].Br Poult Sci,2017,58(2):116-121.
[16]	SALDANA-CABOVERDE A,KOS L.Roles of endothelin signaling in melanocyte development and melanoma[J].Pigment Cell Melanoma Res,2010,23(2):160-170.
[17]	DAVOODI P,EHSANI A,TORSHIZI R V,et al.New insights into genetics underlying of plumage color[J].Anim Genet,2022, 53(1):80-93.
[18]	BELLEI B,PITISCI A,CATRICALÀ C,et al.Wnt/β-catenin signaling is stimulated by α-melanocyte-stimulating hormone in melanoma and melanocyte cells:implication in cell differentiation[J].Pigment Cell Melanoma Res,2011,24(2):309-325.
[19]	RABBANI P,TAKEO M,CHOU W,et al.Coordinated activation of Wnt in epithelial and melanocyte stem cells initiates pigmented hair regeneration[J].Cell,2011,145(6):941-955.
[20]	NAGARATHINAM S,BAALANN K P.Nevus pigmentosus et pilosus[J].Pan Afr Med J,2021,40:107.
[21]	LEE J E,KIM S Y,JEONG Y M,et al.The regulatory mechanism of melanogenesis by FTY720,a sphingolipid analogue[J].Exp Dermatol,2011,20(3):237-241.
[22]	XIAO D,CHEN X,TIAN R B,et al.Molecular and potential regulatory mechanisms of melanin synthesis in Harmonia axyridis[J]. Int J Mol Sci, 2020,21(6):2088.
[23]	MUNDY N I.A window on the genetics of evolution:MC1R and plumage colouration in birds[J].Proc Biol Sci,2005, 272(1573):1633-1640.
[24]	NAM I S,OH M G,NAM M S,et al.Specific mutations in the genes of MC1R and TYR have an important influence on the determination of pheomelanin pigmentation in Korean native chickens[J].J Adv Vet Anim Res,2021,8(2):266-273.
[25]	XU Q,LIU X M,CHAO Z,et al.Transcriptomic analysis of coding genes and non-coding RNAs reveals complex regulatory networks underlying the black back and white belly coat Phenotype in Chinese Wuzhishan pigs[J].Genes (Basel),2019,10(3):201.
[26]	LI L F,CHANG B D,RONINSON I B,et al.Alteration of skin protein kinase C α protein and mRNA levels during induced mouse hair growth[J].J Dermatol,1999,26(4):203-209.
[27]	周庆红,周灿,郑伟,等.甘蓝型油菜角果长度全基因组关联分析[J].中国农业科学,2017,50(2):228-239.ZHOU Q H,ZHOU C,ZHENG W,et al.Genome wide association analysis of silique length in Brassica napus L[J].Scientia Agricultura Sinica,2017,50(2):228-239.(in Chinese)
[28]	LI Z Z,ZHENG M,ALI ABDALLA B,et al.Genome-wide association study of aggressive behaviour in chicken[J].Sci Rep,2016, 6(1):30981.
[29]	段星海,安炳星,杜丽丽,等.中国肉用西门塔尔牛生长曲线参数的全基因组关联分析[J].畜牧兽医学报,2021,52(5):1267-1277.DUAN X H,AN B X,DU L L,et al.Genome-wide association study of growth curve parameters in Chinese Simmental beef cattle[J].Acta Veterinaria et Zootechnica Sinica,2021,52(5):1267-1277.(in Chinese)
[30]	周李生,杨杰,刘先先,等.苏太猪宰后72 h pH和肉色性状的全基因组关联分析[J].中国农业科学,2014,47(3):564-573.ZHOU L S,YANG J,LIU X X,et al.Genome-wide association analyses for musle pH 72 h value and meat color traits in Sutai Pigs[J].Scientia Agricultura Sinica,2014,47(3):564-573.(in Chinese)
[31]	王晓薇,马青.滩羊毛色的全基因组关联分析[J].浙江农业学报,2020,32(1):28-34.WANG X W,MA Q.Genome-wide association studies for coat color in Tan sheep[J].Acta Agriculturae Zhejiangensis,2020,32(1):28-34.(in Chinese)
[32]	ZHANG H H,XIE T,SHUI Y J,et al.Knockdown of PLCB2 expression reduces melanoma cell viability and promotes melanoma cell apoptosis by altering Ras/Raf/MAPK signals[J].Mol Med Rep,2019,21(1):420-428.
[33]	KORSTANJE R,DESAI J,LAZAR G,et al.Quantitative trait loci affecting phenotypic variation in the vacuolated lens mouse mutant,a multigenic mouse model of neural tube defects[J].Physiol Genomics,2008,35(3):296-304.
[34]	FERNSTROM J D,FERNSTROM M H.Tyrosine,phenylalanine,and catecholamine synthesis and function in the brain[J].J Nutr,2007,137(6 Suppl 1):1539S-1547S.