全基因组关联分析揭示白羽肉鸡孵化性状的遗传基础

doi:10.11843/j.issn.0366-6964.2023.02.011

摘要/Abstract

摘要： 旨在揭示白羽肉鸡孵化性状的遗传基础。本试验以白羽肉鸡A、B两个公鸡群体为素材，测定A群体5个世代（556只）、B群体3个世代（398只）的40周龄种蛋受精率和孵化率，并利用55K SNP芯片对两个群体共954个健康个体进行基因分型。基于系谱和基因型信息计算的A群体的受精率遗传力分别为0.21和0.14。利用混合线性模型对受精率和孵化率进行全基因组关联分析（genome-wide association study，GWAS）。全基因组关联分析共关联到12个显著SNPs，分布在1、9、11、12、20号染色体上，其中有9个SNPs与受精率显著关联，有3个SNPs与孵化率显著关联，对A群体受精率高、低组各4只公鸡睾丸组织进行转录组测序，实时荧光定量PCR验证两个差异表达基因。全基因组关联分析筛选出9个与受精率相关的候选基因COPG1、ADAMTS18、FABP2、CDH8、SLCO4A1、ATP13A3、NELL2、VEZT和IGHMBP2；3个与孵化率相关的候选基因TMEM、SCO1和MYH1A。转录组测序与荧光定量PCR验证了两个与受精率相关的候选基因HMGCLL1和COA6。本研究利用GWAS和转录组测序筛选到孵化性状的相关候选基因，为白羽肉鸡孵化性状的遗传改良提供了重要的分子标记。

关键词: 白羽肉鸡, 受精率, 孵化率, 全基因组关联分析, 转录组测序

Abstract: The study aimed to reveal the genetic basis of hatching traits of white-feathered broilers. In this experiment, two groups of White-feather broilers A and B were used as materials to determine the fertilization rate and hatching rate of 40-week-old eggs in 5 generations (556) of group A and 3 generations (398) of group B. The 55K SNP chip was used to genotype 954 healthy individuals in the two groups. Genetic parameter estimation based on pedigree and single nucleotide polymorphisms (SNPs) showed that the heritability of fertilization rate of group A was 0.21 and 0.14, respectively. Genome-wide association study (GWAS) of fertilization rate and hatching rate was performed using mixed linear model. A total of 12 significant SNPs were obtained by genome-wide association study, distributed on chromosome 1, 9, 11, 12, 20. Among them, 9 SNPs were significantly associated with fertilization rate, 3 SNPs were significantly associated with hatching rate, RNA-seq was performed on testis tissues of 4 roosters in each group with high or low fertilization rate, 2 differentially expressed genes were identified by qPCR. Nine genes related to fertilization rate were identified by genome-wide association study, including COPG1, ADAMTS18, FABP2, CDH8, SLCO4A1, ATP13A3, NELL2, VEZT and IGHMBP2,3 genes related to hatching rate, including TMEM, SCO1 and MYH1A, RNA-seq and qPCR identified two candidate genes, HMGCLL1 and COA6, which were related to fertilization rate. In this study, GWAS and RNA-seq were used to identified candidate genes related to hatching traits, which provided important molecular markers for the genetic improvement of hatching traits of White-feathered broilers.

Key words: White-feather broiler, fertilization rate, hatching rate, GWAS, RNA-seq

中图分类号:

S831.2

张高猛, 丁纪强, 刘昱宏, 郑麦青, 文杰, 赵桂苹, 李庆贺. 全基因组关联分析揭示白羽肉鸡孵化性状的遗传基础[J]. 畜牧兽医学报, 2023, 54(2): 534-544.

ZHANG Gaomeng, DING Jiqiang, LIU Yuhong, ZHENG Maiqing, WEN Jie, ZHAO Guiping, LI Qinghe. Genome-wide Association Study Reveals the Genetic Basis of Hatching Traits in White-feathered Broilers[J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(2): 534-544.

参考文献 34

[1]	RIAZ M F, MAHMUD A, HUSSAIN J, et al.Impact of light stimulation during incubation on hatching traits and post-hatch performance of commercial broilers[J].Trop Anim Health Prod, 2021, 53(1):107.
[2]	LIU K, CAO H Y, DONG X Y, et al.Polymorphisms of pro-opiomelanocortin gene and the association with reproduction traits in chickens[J].Anim Reprod Sci, 2019, 210:106196.
[3]	赵秀华, 王金玉, 李满雨, 等.STATTAT5b基因多态性及其与京海黄鸡繁殖性状的相关分析[J].中国家禽, 2015, 37(22):4-9.ZHAO X H, WANG J Y, LI M Y, et al.Polymorphisms of STAT5b gene and its relationship with reproduction traits in Jinghai Yellow chicken[J].China Poultry, 2015, 37(22):4-9.(in Chinese)
[4]	张涛, 樊庆灿, 张向前, 等.京海黄鸡体组成性状的全基因组关联分析[J].畜牧兽医学报, 2015, 46(9):1502-1514.ZHANG T, FAN Q C, ZHANG X Q, et al.A genome-wide association study on body composition of Jinghai Yellow chicken[J].Acta Veterinaria et Zootechnica Sinica, 2015, 46(9):1502-1514.(in Chinese)
[5]	HU Z L, PARK C A, WU X L, et al.Animal QTLdb:an improved database tool for livestock animal QTL/association data dissemination in the post-genome era[J].Nucleic Acids Res, 2013, 41(D1):D871-D879.
[6]	LIU R R, XING S Y, WANG J, et al.A new chicken 55K SNP genotyping array[J].BMC Genomics, 2019, 20(1):410.
[7]	LIU Y F, TU Y J, ZHANG M, et al.Identification of molecular pathways and candidate genes associated with cocks' comb size trait by genome-wide transcriptome analysis[J].Sci Rep, 2018, 8(1):2015.
[8]	WU N, GAUR U, ZHU Q, et al.Expressed microRNA associated with high rate of egg production in chicken ovarian follicles[J].Anim Genet, 2017, 48(2):205-216.
[9]	SUN Y Y, XUE F G, LI Y L, et al.Differences in semen quality, testicular histomorphology, fertility, reproductive hormone levels, and expression of candidate genes according to sperm motility in Beijing-You chickens[J].Poult Sci, 2019, 98(9):4182-4189.
[10]	LI W, ZHENG M Q, ZHAO G P, et al.Identification of QTL regions and candidate genes for growth and feed efficiency in broilers[J].Genet Sel Evol, 2021, 53(1):13.
[11]	PURCELL S, NEALE B, TODD-BROWN K, et al.PLINK:a tool set for whole-genome association and population-based linkage analyses[J].Am J Hum Genet, 2007, 81(3):559-575.
[12]	ZHOU X, STEPHENS M.Genome-wide efficient mixed-model analysis for association studies[J].Nat Genet, 2012, 44(7):821-824.
[13]	ZHAO Y D, LI M C, KONATÉ M M, et al.TPM, FPKM, or normalized counts? A comparative study of quantification measures for the analysis of RNA-seq data from the NCI patient-derived models repository[J].J Transl Med, 2021, 19(1):269.
[14]	PRICE A L, PATTERSON N J, PLENGE R M, et al.Principal components analysis corrects for stratification in genome-wide association studies[J].Nat Genet, 2006, 38(8):904-909.
[15]	孙艳发, 刘冉冉, 郑麦青, 等.鸡胫长和胫围的全基因组关联分析[J].畜牧兽医学报, 2013, 44(3):358-365.SUN Y F, LIU R R, ZHENG M Q, et al.Genome-wide association study on shank length and shank girth in chicken[J].Acta Veterinaria et Zootechnica Sinica, 2013, 44(3):358-365.(in Chinese)
[16]	PEARSON T A, MANOLIO T A.How to interpret a genome-wide association study[J].JAMA, 2008, 299(11):1335-1344.
[17]	VISSCHER P M, BROWN M A, MCCARTHY M I, et al.Five years of GWAS discovery[J].Am J Hum Genet, 2012, 90(1):7-24.
[18]	DING J Q, YING F, LI Q H, et al.A significant quantitative trait locus on chromosome Z and its impact on egg production traits in seven maternal lines of meat-type chicken[J].J Anim Sci Biotechnol, 2022, 13(1):96.
[19]	ZHAO X R, NIE C S, ZHANG J X, et al.Identification of candidate genomic regions for chicken egg number traits based on genome-wide association study[J].BMC Genomics, 2021, 22(1):610.
[20]	ZHANG H, SHEN L Y, XU Z C, et al.Haplotype-based genome-wide association studies for carcass and growth traits in chicken[J].Poult Sci, 2020, 99(5):2349-2361.
[21]	RUAN D L, ZHUANG Z W, DING R R, et al.Weighted single-step GWAS identified candidate genes associated with growth traits in a duroc pig population[J].Genes (Basel), 2021, 12(1):117.
[22]	ZHANG Y F, ZHANG J J, GONG H F, et al.Genetic correlation of fatty acid composition with growth, carcass, fat deposition and meat quality traits based on GWAS data in six pig populations[J].Meat Sci, 2019, 150:47-55.
[23]	VAN DEN BERG I, HO P N, NGUYEN T V, et al.GWAS and genomic prediction of milk urea nitrogen in Australian and New Zealand dairy cattle[J].Genet Sel Evol, 2022, 54(1):15.
[24]	KEOGH K, CARTHY T R, MCCLURE M C, et al.Genome-wide association study of economically important traits in Charolais and Limousin beef cows[J].Animal, 2021, 15(1):100011.
[25]	ZEPEDA-BATISTA J L, NÚÑEZ-DOMINGUEZ R, RAMÍREZ-VALVERDE R, et al.Discovering of genomic variations associated to growth traits by GWAS in braunvieh cattle[J].Genes (Basel), 2021, 12(11):1666.
[26]	MEITEI H T, JADHAV N, LAL G.CCR6-CCL20 axis as a therapeutic target for autoimmune diseases[J].Autoimmun Rev, 2021, 20(7):102846.
[27]	BELLO S F, ADEOLA A C, NIE Q H.The study of candidate genes in the improvement of egg production in ducks-a review[J].Poult Sci, 2022, 101(7):101850.
[28]	TSAI I T, WU C C, HUNG W C, et al.FABP1 and FABP2 as markers of diabetic nephropathy[J].Int J Med Sci, 2020, 17(15):2338-2345.
[29]	LIN X T, WU T J, WANG L Y, et al.ADAMTS18 deficiency leads to preputial gland hypoplasia and fibrosis in male mice[J]. Reprod Biol, 2021, 21(3):100542.
[30]	KIYOZUMI D, NODA T, YAMAGUCHI R, et al.NELL2-mediated lumicrine signaling through OVCH2 is required for male fertility[J].Science, 2020, 368(6495):1132-1135.
[31]	XU W W, WANG Z Z, QU Y Q, et al.Genome-wide association studies and haplotype-sharing analysis targeting the egg production traits in Shaoxing duck[J].Front Genet, 2022, 13:828884.
[32]	ZENG M R, LIU J L, YANG W X, et al.Multiple-microarray analysis for identification of hub genes involved in tubulointerstial injury in diabetic nephropathy[J].J Cell Physiol, 2019, 234(9):16447-16462.
[33]	KIM J, LEE T, KIM T H, et al.An integrated approach of comparative genomics and heritability analysis of pig and human on obesity trait:evidence for candidate genes on human chromosome 2[J].BMC Genomics, 2012, 13:711.
[34]	ZHONG O, JI L, WANG J Y, et al.Association of diabetes and obesity with sperm parameters and testosterone levels:a meta-analysis[J].Diabetol Metab Syndr, 2021, 13(1):109.