PRLR 基因InDel检测及其与陕北白绒山羊生长性状的关联分析

王聪亮 张政轩 付琪 魏宇新 白晶晶 宋晓越 李陇平 屈雷 朱海鲸

摘要 [目的]探究催乳素受體（prolactin receptor，PRLR）基因InDel突变对绒山羊种群生长特性的影响。[方法]以陕北白绒山羊为研究对象，采用PCR扩增和测序技术法检测分析陕北白绒山羊 PRLR 基因的多态性，利用独立样本 t 检验和单因素方差法分析 PRLR基因 不同基因型与陕北白绒山羊生长性状之间的相关性。[结果]在所检测的1 176只陕北白绒山羊群体中， PRLR 基因第2内含子存在1个16 bp InDel突变，形成了纯合插入型（II）、杂合型（ID）与纯合缺失型（DD）3种基因型;在育成羊群体中，该InDel突变与生长性状无显著关联（ P >0.05），在成年羊群体中，InDel突变与管围性状显著相关（ P< 0.05）;在育成羊和成年羊全部群体中，InDel突变仍与管围性状显著相关（ P< 0.05）。[结论] PRLR 基因16 bp InDel突变可能会对陕北白绒山羊生长发育产生显著影响，该位点可作为陕北白绒山羊辅助育种的候选分子标记。

关键词 绒山羊;PRLR基因;InDel;生长性状;关联分析

中图分类号 S827  文献标识码 A  文章编号 0517-6611（2022）11-0085-04

doi：10.3969/j.issn.0517-6611.2022.11.022

开放科学（资源服务）标识码（OSID）：

InDel Detection of  PRLR  Gene and Its Association Analysis with Growth Traits of Shanbei White Cashmere Goats

WANG Cong-liang，ZHANG Zheng-xuan，FU Qi et al

（Shaanxi Province Engineering & Technology Research Center of Cashmere Goats，Yulin University，Yulin，Shaanxi 719000）

Abstract [Objective]In order to investigate the effect of prolactin receptor（PRLR） gene InDel mutation on the growth characteristics of cashmere goat population.[Method]The Shanbei white cashmere goats were used as the research object，and used PCR amplification and direct sequencing techniques to detect the polymorphism of  PRLR  gene in Shanbei white cashmere goats，and analyzed the correlation between different genotypes of  PRLR  gene and growth traits of Shanbei white cashmere goats using independent sample  t -test and one-way ANOVA.[Result]The results showed that there was one 16 bp InDel mutation in intron 2 of the  PRLR  gene in the 1 176 Shanbei white cashmere goats population，and there were three genotypes：homozygous insertion type（II），heterozygous type（ID） and homozygous deletion type（DD），respectively.The InDel mutation was not significantly correlated with growth traits in the yearling goats population（ P >0.05），and was significantly correlated with the cannon circumference trait in the adult goats population（ P< 0.05）;In all of the goat population this InDel mutation was still significantly correlated with the cannon circumference trait（ P< 0.05）.[Conclusion]The 16 bp InDel mutation in  PRLR  gene may impose significant effect on the growth and development of Shanbei white cashmere goats，and this locus can be used as a candidate molecular marker for the auxiliary breeding of high-yielding Shanbei white cashmere goats.

Key words Cashmere goats; PRLR  gene;InDel;Growth traits;Association analysisB38838FF-7371-4061-8A06-23ACEAFF21C7

催乳素（Prolactin，PRL）主要是垂体前叶分泌的一种肽类激素[1]，仅包含1个跨膜区域的膜结合蛋白，最早发现于乳腺、卵巢、睾丸等组织，PRL通过与其受体（PRLR）结合作用于靶器官，参与调控乳腺发育、黄体生成、调节渗透压等及繁殖相关激素分泌等多种生理活动[2-3]。按照PRLR受体胞内域氨基酸数目和胞质区域组成可分为长型PRLR、中型PRLR和短型PRLR，其功能各不相同[4]。据报道，PRL在大鼠、绵羊和牛等多个物种内分布广泛，通过与不同PRLR受体亚型结合介导，作用于相应靶细胞质膜，进而发挥其生物学作用[5-6]。

PRLR 基因作为泌乳相关基因，研究人员发现， PRLR基 因与母猪总产仔数和产活仔数显著相关，可作为猪高产仔数性状的候选基因[7]，国内外研究人员就 PRLR 基因对畜禽泌乳和繁殖性状的影响展开大量研究。如， PRLR 基因侧翼区一处SNP突变与疆岳驴平均日泌乳量和乳蛋白率显著相关，可作为乳用型疆岳驴选育的有效DNA标记[8]。此外， PRLR基因 也可作为高产仔数和海门山羊选育的重要候选基因[9-10]。但 PRLR 基因多态并不是都与动物繁殖相关，如 PRLR基 因多个SNP突变位点与湖羊和济宁青山羊等数个品种的产羔性状无显著关联[11-12]。陕北白绒山羊作为优秀绒肉兼用品种，目前关于 PRLR基 因与陕北白绒山羊生产相关性状的研究较少[13]。为此，笔者采用PCR扩增和测序技术探究 PRLR 基因InDel多态性及其与不同陕北白绒山羊群体生产性状的相关性，以期为绒山羊的高效选育提供参考。

1 材料与方法

1.1 试验动物与基因组DNA

试验羊由陕西应马安绒山羊养殖公司提供，包括618只育成羊（1.0～1.5岁）和558只成年羊（1.5～4.0岁），均为雌性。剪取试验羊耳尖部约0.4 cm2耳组织储存在含有70 %乙醇的1.5 mL Eppedorff管中，同时测量试验羊体重、体高、体长、胸围、髋宽、管围、十字部高、胸深、胸宽、绒细度10个性状数据[14]。同时提取试验羊耳组织DNA[15]，-80 ℃保存备用。

1.2 主要试剂 琼脂糖、Marker 1和PCR Master Mix，均购自天根生化科技（北京）有限公司;EB核酸染料，购自生工生物工程（上海）股份有限公司。

1.3  PRLR 基因InDel目的片段的PCR扩增 参照文献[13]，选取 PRLR 基因第2内含子16 bp InDel突变位点（NC_030827.1.g.61200-61206，61209-61214，61298-61300 del GAGGGC-GGGAGGG-AAG）作为候选变异位点，并利用文献中的引物对候选变异位点进行多态性检测。 PRLR 基因的PCR扩增体系：PCR Master Mix 6.5 μL，上下游引物各0.3 μL，模板DNA 0.7 μL，ddH2O补至13 μL。PCR扩增程序：95 ℃预变性5 min;94 ℃变性30 s，56 ℃退火30 s，72 ℃延伸20 s，共38个循环;72 ℃再延伸10 min，4 ℃结束反应保存，扩增目的片段长度分别为143 bp和127 bp。PCR扩增产物进行3%琼脂糖凝胶电泳并鉴定分型，引物合成与PCR产物测序均由生工生物工程（上海）股份有限公司完成。

1.4 数据统计

利用Microsoft Excel 2018软件计算 PRLR 基因InDel突变位点的基因型频率和等位基因频率，使用Nei方法计算基因纯合度（Ho）、杂合度（He）、有效等位基因數（Ne）和多态信息含量（PIC）等群体遗传学参数，使用在线平台SHEsis（https：//analysis.bio-x.cn）计算 PRLR 基因多态性是否处于哈代温伯格（HWE）平衡[16]。其生物统计模型为 Yijk=μ+Gi+Eij，其中μ为总体平均值，Gi为基因型固定效应，Eij 为随机误差。利用SPSS 23.0软件，分别运用独立样本 t 检验和单因素方差方法完成 PRLR基 因陕北白绒山羊育成羊和成年羊群体不同基因型与生长性状的关联分析;再利用Microsoft Excel 2018软件计算不同群体各个性状表型值的

平均值，用个体表型值减去平均值得到差值，运用单因素方

差方法分析不同基因型与陕北白绒山羊群体生长性状之间的差异显著性，其中育成羊和成年羊为一个整体，结果以“平均值±标准误”表示。

2 结果与分析

2.1  PRLR 基因InDel检测与分型 扩增产物经3.5%琼脂糖凝胶电泳，电泳结果见图1。陕北白绒山羊 PRLR 基因16 bp InDel位点（NC_030827.1.g.61200-61206，61209-61214，61298-61300 del GAGGGC-GGGAGGG-AAG）具有多态性，分别产生了纯合插入型（II，143 bp）、杂合型（ID，143 bp和127 bp）与纯合缺失型（DD，127 bp）3种基因型，产物测序结果见图2，测序结果与PCR扩增结果一致。

2.2  PRLR 基因遗传多样性  PRLR 基因16 bp InDel突变位点等位基因频率、基因型频率、多态性信息含量（PIC）等群体遗传学参数见表1。结果显示， PRLR 基因16 bp InDel突变位点在陕北白绒山羊育成羊、成年羊和全部群体中均属于低度多态（PIC<0.25），I型等位基因频率均高于D型;该突变位点在成年羊群体处于哈代温伯格平衡状态（ P >0.05），但在育成羊和全部群体中偏离了哈代温伯格平衡（ P< 0.05）。

2.3 不同陕北白绒山羊群体 PRLR 基因多态性与生长性状的相关性分析B38838FF-7371-4061-8A06-23ACEAFF21C7

通过SPSS 23.0软件对 PRLR 基因InDel位点多态性与不同陕北白绒山羊群体生长性状进行相关性分析，结果见表2～4。在陕北白绒山羊育成羊群体中 （n =618），由于体重和绒细度性状DD基因型只有2个样本，故运用独立样本 t 检验法对其进行相关性分析，结果显示不同基因型与生长性状均无显著关联（表2）;在成年羊群体中 （n =558），ID基因型山羊的管围显著高于II基因型（ P< 0.05）（表3）;在育成羊和成年羊全部群体中（ n =1 176），ID基因型山羊的管围仍显著高于II基因型（ P< 0.05），优势基因型与成年羊群体相同，均为ID型（表4）。

3 讨论

陕北白绒山羊作为国内知名绒肉兼用山羊品种，其中陕北地区是重要养殖区域，产肉性能及遗传稳定性差等是养殖过程中的主要问题[17]。相对传统育种方法而言，分子标记辅助选择技术（MAS）具有选择准确度高、育种周期短等优势[18]，目前已广泛应用于家畜育种，用于改良家畜经济性状，提高生产性能[19-20]。

多态信息含量等与基因遗传多样性往往具有直接关系[21]。该突变多态信息含量显示，在陕北白绒山羊不同群体中均属于低度多态水平（PIC<0.25），说明群体遗传变异较小，可适当加强选择强度，但在育成羊和全部群体中偏离了哈代温伯格平衡（ P< 0.05），说明可能受到遗传漂变及人工选育的影响，而陕北白绒山羊正是人工选育与杂交形成的新品种，等位基因流失可能是其主要影响因素[22]。将 PRLR 基因16 bp InDel突变位点多态性与陕北白绒山羊的生长性状进行相关性分析发现，该突变位点与陕北白绒山羊育成羊生长性状无显著关联，但与成年羊管围性状显著相关;当把育成羊和成年羊作为一个群体分析时，该InDel突变位点仍与管围性状显著相关，优势基因型与成年羊群体相同，均为ID型，具有很好的一致性。在成年羊群体中，该突变位点能显著影响山羊的管围性状，但对育成羊生长性状无显著影响，推测可能存在以下原因：①生长发育阶段不同，育成羊尚未发育完善;② PRLR 基因表达存在时序性，如 PRLR 基因在牦牛生长发育后期睾丸组织中的表达量极显著高于前期[23]，同时5日龄与8周龄雌性大鼠肝脏组织中 PRLR 基因表达量有较大差异[24]，故 PRLR 基因可能在陕北白绒山羊不同生长发育阶段存在差异表达。同时研究显示[25]，羊体质量与管围性状存在显著关联，表明 PRLR 基因该突变位点可能会间接影响陕北白绒山羊体重性状，其中具体机制仍需要进一步探究。

有报道称， PRLR 基因分别定位于绵羊、山羊的16、20号染色体[26-27]，其中在绵羊腹股沟窦、脂肪和山羊肌肉、肝脏及皮下脂肪中存在不同程度的 PRLR 基因表達[28-29]。目前研究人员普遍认为 PRLR 主要通过JAK2-STAT信号转导途径参与机体调控，通过JAK2酪氨酸激酶的活化对STAT5转录因子进行磷酸化，从而形成同型二聚体，与靶基因启动子上的γ-干扰素激活位点序列特异性结合，激活基因转录，诱导并调控蛋白表达[20-31]，除此之外，有丝分裂原激活蛋白（MAP）激酶途径也被认为是参与介导PRL调控作用的重要途径[32]。研究显示，催乳素合成分泌与下丘脑多巴胺的调控有直接关系，多巴胺与催乳素细胞表面的多巴胺D2特异性结合，在小鼠中，缺乏多巴胺D2R受体亚型与其体重增加有直接关系[33]，小鼠敲除 PRLR 基因后骨形成速率和骨密度较敲除前有所下降[34]。研究证实， PRLR 在调控棕色脂肪细胞分化过程中发挥重要作用[35]，在刺激NIH3T3细胞快速增殖的同时还能有效促进细胞分裂蛋白激酶的激活[36]。同时研究人员发现， PRLR在 牦牛等动物肝脏、骨骼肌和肾脏均有较高表达水平[23，37]，并 且PRLR 基因p.H406R和p.M480T 2个错义突变位点多态与大白猪背部脂肪有显著影响[38]，说 明PRLR 基因可能参与早期动物胚胎及机体的生长发育。

综上所述， PRLR 基因16 bp InDel突变对陕北白绒山羊管围性状有显著影响，可作为陕北白绒山羊优良生长性状选育的重要候选基因。

4 结论

该研究验证了 PRLR 基因在陕北白绒山羊群体中存在1个16 bp的插入缺失，并发现该突变位点对陕北白绒山羊群体部分生长性状存在显著关联，提示该位点可作为陕北白绒山羊分子标记辅助育种的有效DNA标记。

参考文献

[1]

邵顺成，李新海，闫背背，等.5个绵羊群体PRLR基因多态性与产羔数的关联分析[J].中国草食动物科学，2021，41（4）：14-18.

[2] BACHELOT A，BINART N.Reproductive role of prolactin[J].Reproduction，2007，133（2）：361-369.

[3] TERMAN A.Effect of the polymorphism of prolactin receptor（PRLR） and leptin（LEP） genes on litter size in Polish pigs[J].J Anim Breed Genet，2005，122（6）：400-404.

[4] 张芬芬，王长法，黄金明，等.催乳素受体的研究进展[J].家畜生态学报，2009，30（3）：6-9.

[5] DALRYMPLE A，EDERY M，JABBOUR H N.Sequence and functional characterisation of the marmoset monkey（ Callithrix jacchus）  prolactin receptor：Comparative homology with the human long-form prolactin receptor[J].Mol Cell Endocrinol，2000，167（1/2）：89-97.B38838FF-7371-4061-8A06-23ACEAFF21C7

[6] COSENZA G，IANNACCONE M，AUZINO B，et al.Remarkable genetic diversity detected at river buffalo prolactin receptor（ PRLR）  gene and association studies with milk fatty acid composition[J].Anim Genet，2018，49（3）：159-168.

[7] VINCENT A L，TUGGLE C K，ROTHSCHILD M F，et al.The prolactin receptor gene is associated with increased litter size in pigs[J].Swine Res Rep，1998，27：15-18.

[8] 毕兰舒，肖海霞，臧长江，等.疆岳驴催乳素受体基因多态性及其与泌乳性状间的相关性分析[J].中国畜牧兽医，2017，44（1）：180-185.

[9] 范一萍，王晓梅，贾启涛，等.猪PRLR和FSHβ基因多态性检测及其与繁殖性状的关联分析[J].中国畜牧兽医，2019，46（3）：774-781.

[10] LI Y J，ZHANG L，SHANG L Q，et al.Genetic polymorphisms at three loci of  PRLR  and  FSHR  gene correlate with litter size in Chinese haimen goat[J].J Anim Vet Adv，2010，9（22）：2835-2838.

[11] 郭海燕，金银，李拥军，等. PRLR 基因多态性与湖羊产羔性能的关系[J].扬州大学学报（农业与生命科学版），2017，38（1）：57-61.

[12] 张圆，袁婧，高一珊，等.济宁青山羊 PRLR 基因的多态性及其与主要经济性状的关联分析[J].安徽农业科学，2020，48（16）：94-97.

[13] LIU X F，MA L，WANG M，et al.Two indel variants of prolactin receptor（ PRLR）  gene are associated with growth traits in goat[J].Anim Biotechnol，2020，31（4）：314-323.

[14] WANG Z，ZHANG X，JIANG E，et al.InDels within caprine  IGF2BP1  intron 2 and the 3′-untranslated regions are associated with goat growth traits[J].Anim Genet，2020，51（1）：117-121.

[15] LAN X Y，PAN C Y，CHEN H，et al.An AluI PCR-RFLP detecting a silent allele at the goat POU1F1 locus and its association with production traits[J].Small Ruminant Res，2007，73（1/2/3）：8-12.

[16] CUI Y，YAN H L，WANG K，et al.Insertion/deletion within the  KDM6A  gene is significantly associated with litter size in goat[J].Front Genet，2018，9：1-11.

[17] 周寒，陈继明，刘瑞霞，等.陕北白绒山羊产业发展现状调查研究[J].畜牧与饲料科学，2012，33（4）：64-65.

[18] 馬彬云，吴建平.遗传标记技术在动物育种中的研究进展[J].甘肃农业大学学报，2004，39（1）：92-96，100.

[19] YANG Q，ZHANG S H，LIU L L，et al.Application of mathematical expectation（ME） strategy for detecting low frequency mutations：An example for evaluating 14-bp insertion/deletion（indel） within the bovine  PRNP  gene[J].Prion，2016，10（5）：409-419.

[20] 田佺，王思寒，刘磊，等.延边黄牛 MYH3和MYH8基 因多态性及其与生长性状的关联分析[J].中国畜牧兽医，2021，48（6）：2122-2128.

[21] 唐晓琴，陈平博，袁婷婷，等.同羊 ETAA1 基因InDel检测及其与生长性状的关联分析[J].安徽农业大学学报，2021，48（2）：245-249.

[22] 刘燕，李学武，王瑞军，等.周岁内蒙古绒山羊不同部位羊绒长度、细度变化规律及其相关性分析[J].黑龙江畜牧兽医，2019（17）：58-62.

[23] 马富龙，褚敏，常永芳，等. LHR、PRLR基 因在不同发育阶段牦牛睾丸中的表达与定位[J].农业生物技术学报，2020，28（4）：672-680.

[24] TANAKA M，SUZUKI M，KAWANA T，et al.Differential effects of sex steroid hormones on the expression of multiple first exons including a novel first exon of prolactin receptor gene in the rat liver[J].J Mol Endocrinol，2005，34（3）：667-673.B38838FF-7371-4061-8A06-23ACEAFF21C7

[25] 赵森，王真，白洋洋，等.绒山羊 MARCH1 基因InDel位点与体尺性状的关联研究[J].家畜生态学报，2020，41（3）：27-32.

[26] JENKINS Z A，HENRY H M，SISE J A，et al.Follistatin（FST），growth hormone receptor（GHR） and prolactin receptor（PRLR） genes map to the same region of sheep chromosome 16[J].Anim Genet，2000，31（4）：280.

[27] MOODY D E，POMP D，BARENDSE W，et al.Assignment of the growth hormone receptor gene to bovine chromosome 20 using linkage analysis and somatic cell mapping[J].Anim Genet，1995，26（5）：341-343.

[28] ALEXANDRE-PIRES G，MARTINS C，GALVO A M，et al.Understanding the inguinal sinus in sheep（ Ovis aries） -morphology，secretion，and expression of progesterone，estrogens，and prolactin receptors[J].Int J Mol Sci，2017，18（7）：1-20.

[29] SHI H P，ZHANG T Y，YI Y Q，et al.Long form PRLR （lPRLR） regulates genes involved in the triacylglycerol synthesis in goat mammary gland epithelial cells[J].Small Rumin Res，2016，139：7-14.

[30] RADHAKRISHNAN A，RAJU R，TULADHAR N，et al.A pathway map of prolactin signaling[J].J Cell Commun Signal，2012，6（3）：169-173.

[31] GALLEGO M I，BINART N，ROBINSON G W，et al.Prolactin，growth hormone，and epidermal growth factor activate Stat5 in different compartments of mammary tissue and exert different and overlapping developmental effects[J].Dev Biol，2001，229（1）：163-175.

[32] DOBOLYI A，OLH S，KELLER D，et al.Secretion and function of pituitary prolactin in evolutionary perspective[J].Front Neurosci，2020，14：1-16.

[33] PEREZ-MILLAN M I，LUQUE G M，RAMIREZ M C，et al.Selective disruption of dopamine D2 receptors in pituitary lactotropes increases body weight and adiposity in female mice[J].Endocrinology，2014，155（3）：829-839.

[34] KELLY P A，BINART N，FREEMARK M，et al.Prolactin receptor signal transduction pathways and actions determined in prolactin receptor knockout mice[J].Biochem Soc Trans，2001，29（Pt2）：48-52.

[35] KIM H，TOYOFUKU Y，LYNN F C，et al.Serotonin regulates pancreatic beta cell mass during pregnancy[J].Nat Med，2010，16（7）：804-808.

[36] DAS R，VONDERHAAR B K.Transduction of prolactin's （PRL） growth signal through both long and short forms of the PRL receptor[J].Mol Endocrinol，1995，9（12）：1750-1759.

[37] 邢光東.鹅催乳素受体基因克隆和多克隆抗体制备及组织表达分析[D].南京：南京农业大学，2011.

[38] FONTANESI L，GALIMBERTI G，CAL D G，et al.Identification and association analysis of several hundred single nucleotide polymorphisms within candidate genes for back fat thickness in Italian Large White pigs using a selective genotyping approach[J].J Anim Sci，2012，90（8）：2450-2464.B38838FF-7371-4061-8A06-23ACEAFF21C7