藏羊BMPRⅡ基因多态性及其单倍型与产羔性状相关分析

doi:10.16431/j.cnki.1671-7236.2022.02.020

摘要/Abstract

摘要： 【目的】研究骨形态发生蛋白受体Ⅱ(bone morphogenetic protein receptor Ⅱ, BMPRⅡ)基因多态性及其单倍型与藏羊产羔性状的相关性。【方法】以433只藏羊母羊为研究对象, 采用改良多重高温连接酶检测反应技术(improved multiple ligase detection reaction, iMLDR)对BMPRⅡ基因9个单核苷酸多态性(single nucleotide polymorphism, SNP)位点在藏羊群体中的多态性进行检测, 使用Haploview 4.2软件分析其连锁不平衡性并构建单倍型, 应用基因关联分析探究BMPRⅡ基因多态性及单倍型与藏羊产羔性状的关联。【结果】藏羊BMPRⅡ基因外显子12中存在6个SNPs, 分别为: g.90192 T>C、g.142532 C>T、g.142614 A>G、g.142751 T>C、g.143138 A>G和g.143189 G>A, 外显子3、9和13中各存在1个SNP, 分别为: g.143570 C>G、g.124843 A>C和g.145233 A>G, 这9个SNPs均为同义突变。9个SNPs中, 除g.90192 T>C外, 均存在3种基因型。多态信息含量分析显示, 群体在g.90192 T>C、g.142614 A>G和g.145233 A>G位点处于低度多态(PIC<0.25), 在g.124843 A>C、g.142532 C>T、g.142751 T>C、g.143138 A>G、g.143189 G>A和g.143570 C>G位点均处于中度多态(0.25<PIC<0.5)。χ²适合性检验结果显示, 所有突变位点均未偏离哈代-温伯格平衡状态。相关性分析显示, 不同位点不同基因型与藏羊产羔性状均无显著相关(P>0.05), 但g.142532 C>T位点CT基因型平均产羔数高于CC和TT基因型, g.142751 T>C位点CC基因型平均产羔数高于TT和TC基因型, g.143189 G>A位点GA基因型平均产羔数高于GG和AA基因型, g.143570 C>G位点CG基因型平均产羔数高于CC和GG基因型, g.145233 A>G位点GG基因型平均产羔数高于AA和AG基因型, 差异均不显著(P>0.05)。BMPRⅡ基因中除g.90192 T>C位点外的8个SNPs位点在藏羊群体中共形成14种单倍型(H1~H14), 单倍型与藏羊产羔数间均无显著相关(P>0.05)。【结论】 BMPRⅡ基因g.142532 C>T、g.142751 T>C、g.143189 G>A、g.143570 C>G和g.145233 A>G位点对藏羊产羔数有一定潜在的影响。

关键词: 藏羊; BMPRⅡ基因; 单核苷酸多态性; 单倍型; 产羔数

Abstract: 【Objective】 The aim of this study was to investigate the relationship between the polymorphism of bone morphogenetic protein receptor Ⅱ (BMPRⅡ) gene and its haplotype lambing traits of Tibetan sheep.【Method】 433 Tibetan ewes were used to detect the polymorphisms of 9 single nucleotide polymorphisms (SNPs) of BMPRⅡ gene in Tibetan sheep population by improved multiple ligase detection reaction (iMLDR).Haploview 4.2 software was used to analyze the linkage disequilibrium and construct haplotypes.Gene association analysis was used to explore the association between BMPRⅡ gene polymorphism and haplotype and lambing traits of Tibetan sheep.【Result】 The results showed that there were g.90192 T>C, g.142532 C>T, g.142614 A>G, g.142751 T>C, g.143138 A>G and g.143189 G>A in exon 12 of BMPRⅡ gene in Tibetan sheep, and there were one mutation site each in exons 3, 9 and 13, g.143570 C>G, g.124843 A>C and g.145233 A>G, respectively.All the 9 SNPs were synonymous mutations.Among the 9 SNPs, there were 3 genotypes except g.90192 T>C.Polymorphism information content analysis showed that the population was low polymorphic at g.90192 T>C, g.142614 A>G and g.145233 A>G (PIC < 0.25).g.124843 A>C, g.142532 C>T, g.142751 T>C, g.143138 A>G, g.143189 G>A and g.143570 C>G were all moderately polymorphic (0.25 < PIC < 0.5).The results of Chi-square fitness test showed that all mutation sites did not deviate from Hardy-Weinberg equilibrium state.Correlation analysis showed that there was no significant correlation between different genotypes at different SNPs and lambing traits of Tibetan sheep (P > 0.05), but the average lambing number of CT genotype at g.142532 C>T was higher than CC and TT genotypes, the average lambing number of CC genotype at g.142751 T>C was higher than TT and TC genotypes, and the average lambing number of GA genotype at g.143189 G>A was higher than GG and AA genotypes, the average lambing number of CG genotype at 143570 C>G was higher than that of CC and GG genotypes, and the average lambing number of GG genotype at g.145233 A>G was higher than that of AA and AG genotypes, and these differences were not significant (P>0.05).8 SNPs in BMPRⅡ gene except for g.90192 T>C formed 14 haplotypes (H1-H14) in Tibetan sheep, and there was no significant correlation between haplotypes and lamb number (P > 0.05).【Conclusion】 BMPRⅡ gene g.142532 C>T, g.142751 T>C, g.143189 G>A, g.143570 C>G and g.145233 A>G might have some potential effects on the number of lambing in Tibetan sheep.

Key words: Tibetan sheep; BMPRⅡ gene; single nucleotide polymorphism; haplotype; litter size

中图分类号:

S826.3

李明明, 吴庆伟, 邓玉婷, 靳生伟, 贺娜, 张军霞. 藏羊BMPRⅡ基因多态性及其单倍型与产羔性状相关分析[J]. 中国畜牧兽医, 2022, 49(2): 587-597.

LI Mingming, WU Qingwei, DENG Yuting, JIN Shengwei, HE Na, ZHANG Junxia. Correlation Analysis of BMPRⅡ Gene Polymorphism and Its Haplotype with Lambing Traits in Tibetan Sheep[J]. China Animal Husbandry and Veterinary Medicine, 2022, 49(2): 587-597.

参考文献

[1] UPADHYAY D, KAUR T, KAPILA R, et al. Repertoire of structure-activity-based novel modified peptides elicits enhanced osteogenic potential[J]. Journal of Agricultural and Food Chemistry, 2020, 68(31): 8308-8320.
[2] KISHIGAMI S, MISHINA Y. BMP signaling and early embryonic patterning[J]. Cytokine & Growth Factor Reviews, 2005, 16(3): 265-278.
[3] BEYER T A, NARIMATSU M, WEISS A, et al. The TGFβ superfamily in stem cell biology and early mammalian embryonic development[J]. Biochimica et Biophysica Acta (BBA)-General Subjects, 2013, 1830(2): 2268-2279.
[4] CARREIRA A C, ALVES G G, ZAMBUZZI W F, et al. Bone morphogenetic proteins: Structure, biological function and therapeutic applications[J]. Archives of Biochemistry and Biophysics, 2014, 561: 64-73.
[5] BEPPU H, KAWABATA M, HAMAMOTO T, et al. BMP type Ⅱ receptor is required for gastrulation and early development of mouse embryos[J]. Developmental Biology, 2000, 221(1): 249-258.
[6] NAGASHIMA T, LI Q, CLEMENTI C, et al. BMPR2 is required for postimplantation uterine function and pregnancy maintenance[J]. Journal of Clinical Investigation, 2013, 123(6): 2539-2550.
[7] KUMAR S, LAKSHMI D H, SINGH J N, et al. Expression and functional role of bone morphogenetic proteins (BMPs) in placenta during different stages of pregnancy in water buffalo (Bubalus bubalis)[J]. General and Comparative Endocrinology, 2020, 285: 113249.
[8] SELVARAJU S, FOLGER J K, GUPTA P S P, et al. Stage-specific expression and effect of bone morphogenetic protein 2 on bovine granulosa cell estradiol production: Regulation by cocaine and amphetamine regulated transcript[J]. Domestic Animal Endocrinology, 2013, 44(3): 115-120.
[9] RAJESH G, PAUL A, MISHRA S R, et al. Expression and functional role of bone morphogenetic proteins (BMPs) in cyclical corpus luteum in buffalo (Bubalus bubalis)[J]. General and Comparative Endocrinology, 2017, 240: 198-213.
[10] FATEHI A N, VAN DEN HURK R, COLENBRANDER B, et al. Expression of bone morphogenetic protein 2 (BMP2), BMP4 and BMP receptors in the bovine ovary but absence of effects of BMP2 and BMP4 during IVM on bovine oocyte nuclear maturation and subsequent embryo development[J]. Theriogenology, 2005, 63(3): 872-889.
[11] COSTA J J N, PASSOS M J, LEITÃO C C F, et al. Levels of mRNA for bone morphogenetic proteins, their receptors and SMADs in goat ovarian follicles grown in vivo and in vitro[J]. Reproduction, Fertility and Development, 2012, 24(5): 723.
[12] LIMA I M T, BRITO I R, ROSSETTO R, et al. BMPRIB and BMPRⅡ mRNA expression levels in goat ovarian follicles and the in vitro effects of BMP-15 on preantral follicle development[J]. Cell and Tissue Research, 2012, 348(1): 225-238.
[13] VIANA A C N P, FREITAS J L S, MAGALHÃES L C, et al. Ovarian gene expression in collared peccary (Pecari tajacu Linnaeus, 1758) subjected to gonadotropin stimulation protocols[J]. Reproduction in Domestic Animals, 2021, 56(2): 351-359.
[14] MOROS N C, IZQUIERDO R M J, LI Y, et al. Relative transcript abundance in porcine cumulus cells collected from different sized follicles[J]. Reproduction in Domestic Animals, 2021, 56(2): 374-380.
[15] PARADIS F, NOVAK S, MURDOCH G K, et al. Temporal regulation of BMP2. BMP6. BMP15. GDF9. BMPR1A, BMPR1B, BMPR2 and TGFBR1 mRNA expression in the oocyte, granulosa and theca cells of developing preovulatory follicles in the pig[J]. Reproduction, 2009, 138(1): 115-129.
[16] HU X, POKHAREL K, PEIPPO J, et al. Identification and characterization of miRNAs in the ovaries of a highly prolific sheep breed[J]. Animal Genetics, 2016, 47(2): 234-239.
[17] BERTOLDO M J, DUFFARD N, BERNARD J, et al. Effects of bone morphogenetic protein 4 (BMP4) supplementation during culture of the sheep ovarian cortex[J]. Animal Reproduction Science, 2014, 149(3-4): 124-134.
[18] FEARY E S, JUENGEL J L, SMITH P, et al. Patterns of expression of messenger RNAs encoding GDF9. BMP15. TGFBR1. BMPR1B, and BMPR2 during follicular development and characterization of ovarian follicular populations in ewes carrying the woodlands FecX2W mutation1[J]. Biology of Reproduction, 2007, 77(6): 990-998.
[19] KUMAR S, RAJPUT P K, BAHIRE S V, et al. Differential expression of BMP/SMAD signaling and ovarian-associated genes in the granulosa cells of FecB introgressed GMM sheep[J]. Systems Biology in Reproductive Medicine, 2020, 66(3): 185-201.
[20] KYASARI O R, VALOJERDI M R, FARROKHI A, et al. Expression of maturation genes and their receptors during in vitro maturation of sheep COCs in the presence and absence of somatic cells of cumulus origin[J]. Theriogenology, 2012, 77(1): 12-20.
[21] ONGARO L, ZHOU X, CUI Y, et al. Gonadotrope-specific deletion of the BMP type 2 receptor does not affect reproductive physiology in mice[J]. Biology of Reproduction, 2020, 102(3): 639-646.
[22] WANG W, ZHANG H Y, HE Y, et al. Expression of BMP receptors in porcine granulosa cells (GCs) and their regulation by luteinizing hormone (LH)[J]. Agricultural Sciences in China, 2011, 10(2): 289-295.
[23] KAYANI A R, GLISTER C, KNIGHT P G. Evidence for an inhibitory role of bone morphogenetic protein(s) in the follicular-luteal transition in cattle[J]. Reproduction, 2009, 137(1): 67-78.
[24] SANTOS-BIASE W K F, BIASE F H, BURATINI J, et al. Single nucleotide polymorphisms in the bovine genome are associated with the number of oocytes collected during ovum pick up[J]. Animal Reproduction Science, 2012, 134(3-4): 141-149.
[25] 文禹粱. 绵羊繁殖相关基因的组织表达及其多态性与产羔数的关联分析[D]. 兰州: 甘肃农业大学, 2019.
WEN Y L. Correlation analysis of tissue expression and polymorphism of sheep reproductive related genes with lamb number[D]. Lanzhou: Gansu Agricultural University, 2019. (in Chinese)
[26] GUO X F, HU W P, LANG X Z, et al. Two single nucleotide polymorphisms sites in α1-AT gene and their association with somatic cell score in Chinese Holstein cows[J]. Journal of Biological Research, 2017, 24(1): 8.
[27] DA CUNHA E V, DE SOUZA G B, PASSOS J R S, et al. Effects of bone morphogenetic protein 4 (BMP4) on in vitro development and survival of bovine preantral follicles enclosed in fragments ovarian tissue[J]. Zygote, 2017, 25(3): 256-264.
[28] IKEDA Y, HASEGAWA A, TSUBAMOTO H, et al. Effects of gremlin-2 on the transition of primordial follicles during early folliculogenesis in the human ovary[J]. European Journal of Obstetrics & Gynecology and Reproductive Biolology 2016, 203: 72-77.
[29] YANG Z, YANG X, LIU G, et al. Polymorphisms in BMPR-1B gene and their association with litter trait in Chinese Hu sheep[J]. Animnal Biotechnology, 2020, 12: 1-10.
[30] KIMCHI-SARFATY C, OH J M, KIM I W, et al. A "Silent" polymorphism in the MDR1 gene changes substrate specificity[J]. Science, 2007, 315(5811): 525-528.
[31] 杜红丽, 崔建勋, 张细权. 利用DNA池和测序技术快速筛查候选效应单倍型[J]. 农业生物技术学报, 2008, 2: 281-285.
DU H L, CUI J X, ZHANG X Q. Rapid screening of candidate efficent haplotypes by DNA pooling and sequencing[J]. Journal of Agricultural Biotechnology, 2008, 2: 281-285. (in Chinese)
[32] AKEY J, JIN L, XIONG M. Haplotypes vs single marker linkage disequilibrium tests: What do we gain?[J]. European Journal of Human Genetics, 2001, 9(4): 291-300.