东湖杂交羊产羔性状的GWAS分析及候选基因GRID2的验证

doi:10.11843/j.issn.0366-6964.2023.11.018

摘要/Abstract

摘要： 旨在研究影响东湖杂交羊产羔性状的关键基因并挑选关键基因的变异位点进行验证,以期为东湖杂交奶绵羊多羔性状的选育提供分子标记。本研究从西北农林科技大学金昌奶绵羊试验示范基地采集身体健康无疾病的东湖杂交二代奶绵羊168只进行全基因组重测序,GEMMA和SnpEff软件进行全基因组关联分析,KOBAS 3.0数据库对注释的基因进行富集分析,LDBlockShow软件对显著性最强SNPs周边的连锁区域进行连锁不平衡分析。使用竞争性等位基因特异性PCR(kompetitive allele-specific PCR,KASP)技术对GRID2基因进行SNP的分型检测,统计该基因两个多态位点的遗传多样性参数,并将其与产羔性状进行关联分析。结果显示,全基因组关联分析确定了6号染色上1 163个达到显著性水平的SNPs位点,通过对候选位点的注释,确定与产羔性状相关的候选基因BMPR1B、PDLIM5、PDHA2、UNC5C、GRID2、NFKB1、TSPAN5、STPG2、PPP3CA。同时结合已有报道和最强关联区域的注释,选择GRID2上的两个SNPs(g.35685218A＞G和g.35685499C＞T)作为候选位点。通过统计学分析发现,本研究中涉及的两个位点均处于P>0.05的哈代-温伯格平衡状态,表明这两个位点均没有受到过强选择,并且这两个位点与东湖杂交羊的产羔性状有显著关联,两个位点的野生型均处于优势状态,且两个位点在胎次增加时均值也增加。综上所述,GRID2基因2个SNPs对东湖杂交奶绵羊产羔性状有显著影响,可以作为分子标记辅助选择的候选位点。本研究为提高东湖杂交奶绵羊产羔率,促进奶绵羊的遗传育种进展提供了坚实的理论依据。

关键词: 东湖杂交奶绵羊, 产羔数, 全基因组关联分析, GRID2基因

Abstract: The experiment aimed to study the key genes affecting lambing traits in East Friesian and Hu crossbred sheep and select the mutation loci of the key genes for validation, to provide molecular markers for the selection and breeding of multiple lambing traits in East Friesian and Hu crossbred sheep. In this study, 168 healthy and disease-free East Friesian and Hu crossbred sheep were collected from the Jinchang Milk Sheep Experimental Demonstration Base of Northwest Agriculture and Forestry University for whole-genome resequencing, genome-wide association study was performed by GEMMA and SnpEff softwares, enrichment analysis of annotated genes was conducted by KOBAS 3.0 database, and analysis of the chain region around the most significant SNPs with LDBlockShow software. Linkage disequilibrium analysis was performed on the linkage regions around the SNPs by LDBlockShow. SNPs were typed in the GRID2 gene using kompetitive allele-specific PCR (KASP), and the genetic diversity parameters of the two polymorphic loci of the gene were counted and analyzed in association with lambing traits. The results showed that genome-wide association analysis identified 1 163 SNP loci at the significance level on chromosome 6, and candidate genes BMPR1B, PDLIM5, PDHA2, UNC5C, GRID2, NFKB1, TSPAN5, STPG2, and PPP3CA were identified as associated with lambing traits through the annotation of the candidate loci. The two SNPs on GRID2 (g.35685218A>G and g.35685499C>T) were also selected as candidate loci by combining the existing reports and annotations of the most strongly associated regions. Statistical analysis revealed that both loci involved in this study were in Hardy-Weinberg equilibrium at P>0.05, indicating that neither locus had been subjected to over-selection and that both loci were significantly associated with lambing traits in East Friesian and Hu crossbred sheep, with wild-type being dominant at both loci and the means of the two loci increasing with the increasing parities. In summary, the two SNPs in the GRID2 gene have significant effects on lambing number traits in East Friesian and Hu crossbred sheep and can be used as candidate loci for molecular marker-assisted selection. This study provides a solid theoretical basis for improving the lambing rate of East Friesian and Hu crossbred sheep and promoting the progress of genetic breeding of milk sheep.

Key words: East Friesian and Hu crossbred sheep, litter size, genome-wide association studies, GRID2 gene

中图分类号:

S826.2

赵雪洋, 李丹妮, 王钰晨, 郭磊, 王丽, 黄杰, 焦仪强, 安小鹏, 张希云, 张磊, 宋宇轩. 东湖杂交羊产羔性状的GWAS分析及候选基因GRID2的验证[J]. 畜牧兽医学报, 2023, 54(11): 4625-4635.

ZHAO Xueyang, LI Danni, WANG Yuchen, GUO Lei, WANG Li, HUANG Jie, JIAO Yiqiang, AN Xiaopeng, ZHANG Xiyun, ZHANG Lei, SONG Yuxuan. GWAS Analysis of Lambing Traits in East Friesian and Hu Crossbred Sheep and Verification of Candidate Gene GRID2[J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(11): 4625-4635.

参考文献 37

[1]	种玉晴.绵羊产羔数性状的分子遗传机理研究[D].武汉:华中农业大学,2021.CHONG Y Q.Study on molecular genetic mechanism of litter size in sheep[D].Wuhan:Huazhong Agricultural University,2021.(in Chinese)
[2]	张仁森,贺建宁,魏彩虹,等.不同繁殖季节绵羊卵巢组织中piRNAs序列特征及表达差异分析[J].中国畜牧兽医,2020, 47(7): 2081-2092.ZHANG R S,HE J N,WEI C H,et al.Characterization and differentially expression analysis of ovine Ovary piRNAs during different reproductive seasons[J].China Animal Husbandry & Veterinary Medicine,2020,47(7):2081-2092.(in Chinese)
[3]	李亭亭,刘秋月,李向臣,等.绵羊经济性状相关基因研究进展及其应用[J].畜牧兽医学报,2022,53(8):2417-2434.LI T T,LIU Q Y,LI X C,et al.Research progress and applications of genes associated with economic traits in sheep[J].Acta Veterinaria et Zootechnica Sinica,2022,53(8):2417-2434.(in Chinese)
[4]	张亮.营养物质对绵羊繁殖的影响[J].黑龙江动物繁殖,2013,21(5):20-21.ZHANG L.Effects of nutrients on sheep reproduction[J].Heilongjiang Journal of Animal Reproduction,2013,21(5):20-21.(in Chinese)
[5]	狄冉,郭晓飞,刘秋月,等.生殖激素对绵羊繁殖性能影响的研究进展[J].家畜生态学报,2015,36(9):1-6.DI R,GUO X F,LIU Q Y,et al.Research progress in effect of reproductive hormones on ovine reproductive performance[J].Acta Ecologae Animalis Domastici,2015,36(9):1-6.(in Chinese)
[6]	贺小云,狄冉,胡文萍,等.绵羊季节性繁殖的神经内分泌研究进展[J].畜牧兽医学报,2018,49(1):18-25.HE X Y,DI R,HU W P,et al.Research progress on neuroendocrinology controlling sheep seasonal reproduction[J].Acta Veterinaria et Zootechnica Sinica,2018,49(1):18-25.(in Chinese)
[7]	常志,汪淼.提高母羊产羔率的措施分析[J].畜禽业,2022,33(6):61-63.CHANG Z,WANG M.Analysis of measures to increase lambing rate of ewe sheep[J].Livestock and Poultry Industry,2022,33(6): 61-63.(in Chinese)
[8]	吕桂萍,文亚洲.提高绵羊繁殖力的关键措施[J].养殖与饲料,2017(6):37-38.LV G P,WEN Y Z.Key measures to improve sheep fertility[J].Animals Breeding and Feed,2017(6):37-38.(in Chinese)
[9]	胡瑞雪,乐祥鹏,李万宏,等.绵羊产羔数相关基因的研究进展[J].基因组学与应用生物学,2017,36(10):4141-4148.HU R X,LE X P,LI W H,et al.A review on the genes related to litter size in sheep[J].Genomics and Applied Biology,2017,36(10): 4141-4148.(in Chinese)
[10]	李晓杰,韩建刚,梁奔梦,等.绵羊重要经济性状可用遗传标记研究进展[J].中国畜牧兽医,2023,50(8):3142-3156.LI X J,HAN J G,LIANG B M,et al.Research advance on available genetic markers for important economic traits in sheep[J].China Animal Husbandry & Veterinary Medicine,2023,50(8):3142-3156.(in Chinese)
[11]	ZHU C J,LI N,CHENG H P,et al.Genome wide association study for the identification of genes associated with tail fat deposition in Chinese sheep breeds[J].Biol Open,2021,10(5):bio054932.
[12]	NICIURA S C M,BENAVIDES M V,OKINO C H,et al.Genome-wide association study for haemonchus contortus resistance in morada nova sheep[J].Pathogens,2022,11(8):939.
[13]	汪文强,赵生国,马利青,等.动物基因组学重测序的应用研究进展[J].畜牧兽医学报,2016,47(10):1947-1953.WANG W Q,ZHAO S G,MA L Q,et al.The research progress and application of resequencing based on animal genomics[J].Acta Veterinaria et Zootechnica Sinica,2016,47(10):1947-1953.(in Chinese)
[14]	苏鹏,黄洋铭,徐天琪,等.全基因组关联分析在绵羊遗传育种研究中的应用[J].中国畜牧杂志,2023,doi:10.19556/j.0258-7033.20220805-07.SU P,HUANG Y M,XU T Q,et al.Application of genome-wide association study in sheep genetics and breeding[J].Chinese Journal of Animal Science,2023,doi:10.19556/j.0258-7033.20220805-07.(in Chinese)
[15]	GEBRESELASSIE G,BERIHULAY H,JIANG L,et al.Review on genomic regions and candidate genes associated with economically important production and reproduction traits in sheep (Ovies aries)[J].Animals,2020,10(1):33.
[16]	KIJAS J W,LENSTRA J A,HAYES B,et al.Genome-wide analysis of the World’s sheep breeds reveals high levels of historic mixture and strong recent selection[J].PLoS Biol,2012,10(2):e1001258.
[17]	欧阳清渊,胡深强,王继文.家禽重要性状的基因组学研究与应用现状[J].畜牧兽医学报,2022,53(3):663-679.OUYANG Q Y,HU S Q,WANG J W.Current genomics research and application of important traits in poultry[J].Acta Veterinaria et Zootechnica Sinica,2022,53(3):663-679.(in Chinese)
[18]	丁鸿祥,王净,刘彦,等.组学技术在绵、山羊繁殖性状研究中的应用[J].中国草食动物科学,2022,42(5):56-62.DING H X,WANG J,LIU Y,et al.The application of omics technologies in relevant studies of sheep and goat reproductive traits[J].China Herbivore Science,2022,42(5):56-62.(in Chinese)
[19]	SMOȽUCHA G,GURGUL A,JASIELCZUK I,et al.A genome-wide association study for prolificacy in three Polish sheep breeds[J]. J Appl Genet,2021,62(2):323-326.
[20]	寸静宇,刘秋月,王翔宇,等.小尾寒羊LHR基因组织表达、多态性及其与产羔性状的关联分析[J].畜牧兽医学报,2019,50(3): 495-506.CUN J Y,LIU Q Y,WANG X Y,et al.Analysis of tissue expression,polymorphism of LHR gene and its association with litter size in small tail han sheep[J].Acta Veterinaria et Zootechnica Sinica,2019,50(3):495-506.(in Chinese)
[21]	周梅,曹晓涵,贺小云,等.绵羊FGF7基因组织表达及其多态性与产羔数之间的关系[J].畜牧兽医学报,2018,49(3):525-533.ZHOU M,CAO X H,HE X Y,et al.Tissue expression and polymorphism of sheep FGF 7 gene and their association with litter size[J].Acta Veterinaria et Zootechnica Sinica,2018,49(3):525-533.(in Chinese)
[22]	王宏浩,靳光,李博,等.数量遗传学在肉牛育种中的应用[J].中国畜禽种业,2023,19(7):74-79.WANG H H,JIN G,LI B,et al.Application of quantitative genetics in beef cattle breeding[J].The Chinese Livestock and Poultry Breeding, 2023,19(7):74-79.(in Chinese)
[23]	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.
[24]	SOUZA C J,MACDOUGALL C,MACDOUGALL C,et al.The booroola (FecB) phenotype is associated with a mutation in the bone morphogenetic receptor type 1 B (BMPR1B) gene[J].J Endocrinol,2001,169(2):R1-R6.
[25]	CHONG Y Q,LIU G Q,JIANG X P.Effect of BMPRIB gene on litter size of sheep in China:A meta-analysis[J].Anim Reprod Sci,2019, 210: 106175.
[26]	GUO X F,WANG X Y,DI R,et al.Metabolic effects of FecB gene on follicular fluid and ovarian vein serum in sheep (Ovis aries)[J]. Int J Mol Sci,2018,19(2):539.
[27]	杨初蕾,张译元,唐红,等.FecB基因在绵羊繁殖性能上的研究进展[J].畜牧兽医学报,2023,54(4):1370-1380.YANG C L,ZHANG Y Y,TANG H,et al.Research progress of FecB gene on reproductive performance of sheep[J].Acta Veterinaria et Zootechnica Sinica,2023,54(4):1370-1380.(in Chinese)
[28]	LI D N,ZHANG L,WANG Y C,et al.FecB mutation and litter size are associated with a 90-base pair deletion in BMPR1B in East Friesian and Hu crossbred sheep[J].Anim Biotechnol,2023,34(4):1314-1323.
[29]	NIU W B,SPRADLING A C.Two distinct pathways of pregranulosa cell differentiation support follicle formation in the mouse ovary[J].Proc Natl Acad Sci USA,2020,117(33):20015-20026.
[30]	DOVERE L,FERA S,GRASSO M,et al.The niche-derived glial cell line-derived neurotrophic factor (GDNF) induces migration of mouse spermatogonial stem/progenitor cells[J].PLoS One,2013,8(4):e59431.
[31]	CHAOUAT G.The Th1/Th2 paradigm:still important in pregnancy?[J].Semin Immunopathol,2007,29(2):95-113.
[32]	王步军,郑飞云,杨孝军,等.HIF-1α和VEGF蛋白在子宫内膜异位症中的表达及临床意义[J].现代实用医学,2011, 23(10):1185-1188.WANG B J,ZHENG F Y,YANG X J,et al.Expression and clinical significance of HIF-1 and VEGF proteins in endometriosis[J]. Modern Practical Medicine,2011,23(10):1185-1188.(in Chinese)
[33]	PANDA P K,SHARAWAT I K,DAWMAN L.GRID2 mutation-related spinocerebellar ataxia type 18:a new report and literature review[J].J Pediatr Genet,2022,11(2):99-109.
[34]	CEYLAN A C,ACAR ARSLAN E,ERDEM H B,et al.Autosomal recessive spinocerebellar ataxia 18 caused by homozygous exon 14 duplication in GRID2 and review of the literature[J].Acta Neurol Belg,2021,121(6):1457-1462.
[35]	E G X,DUAN X H,ZHANG J H,et al.Genome-wide selection signatures analysis of litter size in Dazu black goats using single-nucleotide polymorphism[J].3 Biotech,2019,9(9):336.
[36]	CHEN S,GUO X F,HE X Y,et al.Insight into pituitary lncRNA and mRNA at two estrous stages in small tail han sheep with different FecB genotypes[J].Front Endocrinol (Lausanne),2021,12:789564.
[37]	CHERMUȽA B,BR[XCA.TIF;S*4,JZ]ZERT M,JESETA M,et al.The unique mechanisms of cellular proliferation,migration and apoptosis are regulated through oocyte maturational development-a complete transcriptomic and histochemical study[J].Int J Mol Sci,2018, 20(1):84.