內(nèi)蒙古絨山羊 PRDM6基因插入/缺失(InDel)檢測及其 與生長性狀的關(guān)聯(lián)分析

翟軍軍 王聰亮 劉曉宇 宋曉越 史雷 李新春 朱海鯨 藍賢勇 屈雷

摘要：? PRDM6屬于PRDM蛋白家族的一員，? PRDM6? 基因編碼1個轉(zhuǎn)錄抑制因子，具有與DNA、RNA、蛋白質(zhì)結(jié)合的能力，同時還具有組蛋白賴氨酸 N -甲基轉(zhuǎn)移酶活性。全基因組關(guān)聯(lián)分析結(jié)果顯示， ?PRDM6? 基因遺傳變異與骨密度和體質(zhì)量等性狀顯著相關(guān)，且該基因插入/缺失（Indel）遺傳變異還被發(fā)現(xiàn)與陜北白絨山羊的多個生長性狀顯著相關(guān)。因此，本研究將? PRDM6? 基因作為研究對象，以期挖掘到與內(nèi)蒙古絨山羊生長性狀相關(guān)的關(guān)鍵遺傳變異位點。結(jié)合Ensembl數(shù)據(jù)庫提供的山羊? PRDM6? 遺傳變異信息和已發(fā)表文章中的遺傳變異位點，本研究探究了2個InDel位點（? rs656578433? 、? rs651603667? ）多態(tài)性及其在內(nèi)蒙古絨山羊群體中的分布規(guī)律。發(fā)現(xiàn)在本研究檢測的群體中，只有第一內(nèi)含子區(qū)的12 bp InDel （? rs651603667? ）存在多態(tài)性。不同基因型與生長性狀的關(guān)聯(lián)分析結(jié)果顯示，? rs651603667? 存在3種基因型，即純合插入型（II）、雜合型（ID）和純合缺失型（DD），頻率分別為0.399、0.444、0.159。在育成羊群體中，此? rs651603667? 突變與體質(zhì)量、胸深顯著相關(guān) ?（ P < 0.05）;在成年羊群體中，? rs651603667? 突變與體長、髖寬顯著相關(guān) （ P < 0.05），和體高、胸深極極顯著相關(guān) ?（ P < 0.01）;此外，皮爾遜 （Pearson） 相關(guān)分析結(jié)果顯示，該突變位點顯著影響的生長性狀均與體質(zhì)量存在極顯著相關(guān) ?（ P < 0.01）。綜上，? PRDM6? 基因12 bp InDel突變位點與內(nèi)蒙古絨山羊體長、髖寬和體高、胸深具有顯著或極顯著相關(guān)性，可作為內(nèi)蒙古絨山羊生長性狀選育的有效分子標(biāo)記。

關(guān)鍵詞：? 山羊;?? PRDM6? 基因; 插入/缺失（InDel）; 生長性狀; 關(guān)聯(lián)分析

中圖分類號：? S827??? 文獻標(biāo)識碼： A??? 文章編號：? 1000-4440（2021）05-1244-07

Detection of insertion/deletion （InDel） of?? PRDM6?? gene in Inner Mongolia cashmere goats and its association analysis with growth traits

ZHAI Jun-jun? 1 ， WANG Cong-liang? 1 ， LIU Xiao-yu? 1 ， SONG Xiao-yue? 1 ， SHI Lei? 1 ， LI Xin-chun? 1 ， ?ZHU Hai-jing? 1，2，3 ，? LAN Xian-yong? 4 ， QU Lei? 1，2，3

（1.Shaanxi Province Engineering & Technology Research Center of Cashmere Goat， Yulin University， Yulin 719000， China; 2.Shaanxi Province ‘Four Subjects and One Joint Sheep and Goat Engineering & Technology University & Enterprise Alliance Research Center， Shenmu 719318， China; 3.Shaanxi Haoli Cashmere Goat Technology Development Co.， Ltd.， Zizhou 718499， China; 4.College of Animal Science and Technology， Northwest A&F University， Yangling 712000， China）

Abstract：? PRDM6 is a member of PRDM protein family.?? PRDM6?? gene encodes a transcriptional repressor， which can bind to DNA， RNA and proteins， and also has histone lysine ?N -methyltransferase activity. The results of genome-wide association analysis showed that genetic variations of?? PRDM6?? gene were significantly associated with bone mineral density and body weight. In addition， the insertion/deletion （InDel） genetic variations of?? PRDM6?? were significantly associated with multiple growth traits in Shanbei white cashmere goats. Therefore， this study took?? PRDM6?? as the research object in order to dig out the key genetic variation sites related to the growth traits of Inner Mongolia cashmere goats. Based on the genetic variation information of goat?? PRDM6?? provided by Ensembl database and the genetic variation loci in the published articles， the polymorphism of two InDel loci （ rs656578433， rs651603667 ） and their distribution in Inner Mongolia cashmere goat population were investigated in this study. Only the 12 bp InDel （? rs651603667? ） in the first intron region was polymorphic in the population tested in this study. Different genotypes of 638 Inner Mongolia cashmere goats?? rs651603667?? were identified， and the correlation between different genotypes and growth traits was analyzed. The results showed that there were three genotypes of ?rs651603667， ?namely homozygous insertion type （II）， heterozygous type （ID） and homozygous deletion type （DD）， with frequencies of 0.399（II）， 0.444（ID） and 0.159（DD）. In yearling goat population， the?? rs651603667?? mutation was significantly correlated with body weight and chest depth ?（ P < 0.05）. In adult goat population，?? rs651603667?? mutation was significantly correlated with body length and hip width ?（ P < 0.05）， and extremely significantly correlated with body height and chest depth ?（ P < 0.01）. In addition， the results of Pearson correlation analysis indicated that the growth traits significantly affected by the mutation site were significantly correlated with body weight ?（ P < 0.01）. In conclusion， the 12 bp InDel mutation site of?? PRDM6?? gene is significantly or extremely significantly correlated with body length， hip width， body height and chest depth of Inner Mongolia cashmere goats， which can be used as an effective molecular marker for the breeding of growth traits of Inner Mongolia cashmere goats.

Key words：? goats;?? PRDM6?? gene; insertion/deletion （InDel）; growth traits; association analysis

內(nèi)蒙古絨山羊（ Capra hircus ）作為自然選擇和人工選擇培育而成的地方優(yōu)質(zhì)絨肉兼用山羊，因其柔軟纖細(xì)的羊絨與口感鮮嫩、風(fēng)味獨特的肉質(zhì)而在國內(nèi)外享有盛譽? [1-2] 。長期以來，內(nèi)蒙古絨山羊的選育主要集中在絨毛性狀，忽略了其優(yōu)秀的肉用價值，導(dǎo)致其體型偏小? [3] 。近年來山羊絨價格持續(xù)走低，而山羊肉需求量和價格穩(wěn)定攀升，使得山羊肉逐漸成為養(yǎng)殖戶的重要收入來源? [4] 。研究結(jié)果表明，生長性狀往往與胴體質(zhì)量、產(chǎn)肉量緊密相關(guān)? [5] ，但體型小、產(chǎn)肉性能較差等缺點是目前內(nèi)蒙古絨山羊養(yǎng)殖的主要問題? [6] 。因此，利用有效、實用的現(xiàn)代分子標(biāo)記輔助選擇技術(shù)（MAS）提高和改善此品種的產(chǎn)肉性能及生長性狀是十分必要的。MAS技術(shù)以其準(zhǔn)確、高效等優(yōu)點在畜禽育種工作中可發(fā)揮較大作用? [7] ，此技術(shù)可快速培育高產(chǎn)優(yōu)質(zhì)品種，提高經(jīng)濟性狀? [8-9] 。故利用MAS技術(shù)尋找重要的候選基因并探索其與山羊生長性狀的相關(guān)性具有重要意義。

PRDM6? 基因編碼的PRDM6蛋白是一種轉(zhuǎn)錄抑制因子，屬于組蛋白甲基轉(zhuǎn)移酶超家族成員? [10] 。該蛋白質(zhì)家族成員（PRDM11除外）均包含1個保守的PR域和多個鋅指結(jié)構(gòu)域，其中PR域由于最早發(fā)現(xiàn)于PRDM1和PRDM2而得名? [11-12] 。包含鋅指結(jié)構(gòu)域的蛋白質(zhì)基序可結(jié)合DNA、RNA和蛋白質(zhì)，為調(diào)控基因表達的相關(guān)功能蛋白? [13] 。PR結(jié)構(gòu)域與組蛋白甲基轉(zhuǎn)移酶和組蛋白脫乙酰酶等蛋白質(zhì)相互作用? [14] ，或直接發(fā)揮組蛋白甲基轉(zhuǎn)移酶活性并使目標(biāo)基因的啟動子區(qū)域甲基化，從而影響目標(biāo)基因的轉(zhuǎn)錄和表達? [15-16] 。該蛋白質(zhì)家族獨特的結(jié)構(gòu)特點，使其能夠參與細(xì)胞調(diào)控與機體生長發(fā)育等生命活動的調(diào)節(jié)? [17-18] 。有研究發(fā)現(xiàn)PRDM蛋白家族的成員在胚胎發(fā)育過程中起重要作用，其中? PRDM6? 基因作為轉(zhuǎn)錄阻遏物，是心血管發(fā)育生理控制和維持血管平滑肌細(xì)胞增殖潛能必不可少的因素? [19-20] 。也有報道稱在胚胎和成人內(nèi)皮細(xì)胞中檢測到 ?PRDM6? ??[21] ，其可能是骨質(zhì)疏松和肥胖相關(guān)的潛在多活性基因? [22] 。近期有研究發(fā)現(xiàn) ?PRDM6? 基因的InDel變異與陜北白絨山羊的多個生長性狀顯著相關(guān)，可以用于陜北白絨山羊的選育? [23] 。然而，關(guān)于內(nèi)蒙古絨山羊 ?PRDM6? 基因的遺傳變異及其對生長性狀的影響尚未見報道。

本研究擬檢測? PRDM6? 基因InDel位點多態(tài)性，挖掘并探究與內(nèi)蒙古絨山羊生長性狀相關(guān)的遺傳變異，為內(nèi)蒙古絨山羊生長性狀的遺傳改良提供基礎(chǔ)材料。

1 材料與方法

1.1 試驗材料

試驗羊均來自陜西省榆林市子洲縣某養(yǎng)殖基地，均畜舍飼養(yǎng)，其中包括444只 （10～ 18月齡左右）育成羊和194只 （19～ 36月齡左右）成年羊，均為雌性且無親緣關(guān)系。分別收集耳組織樣品和生長性狀數(shù)據(jù)。育成羊和成年羊在采樣測量期間飼喂環(huán)境和飼喂日糧相同，1日飼喂2次，分別在 8∶00 飼喂、 10∶00 飲水和 16∶00 點飼喂、 18∶00 飲水。

1.2 樣品采集及體尺指標(biāo)測定

剪取試驗羊耳尖約0.4 m? 2 耳組織，放入加有1 ml 70%乙醇的離心管中， -80 ℃ 冰箱保存。采用高鹽法? [24] 從耳組織中提取DNA，采用Nanodrop 2000分光光度計測定DNA樣品的質(zhì)量，DNA樣本稀釋至20 ng/μl并在 -80 ℃ 冰箱保存? [25] 。采樣同時測量內(nèi)蒙古絨山羊的生長性狀數(shù)據(jù)，包括體高（cm）、體長（cm）、胸圍（cm）、髖寬（cm）、管圍（cm）、薦高（cm）、胸深（cm）、胸寬（cm）、體質(zhì)量（kg）。

1.3?? PRDM6? 基因的PCR擴增

參考NCBI數(shù)據(jù)庫和Ensembl數(shù)據(jù)庫（https：//asia.ensembl.org/index.html）中公布的山羊? PRDM6? 基因堿基序列（序列號：NC_030814.1）和已發(fā)表的關(guān)于? PRDM6? 山羊遺傳變異的文獻? [23] ，本研究選取 ?rs656578433? 和? rs651603667? 作為候選遺傳變異位點進行研究。參考文獻 [23] 中的引物，對候選InDel突變位點進行多態(tài)性檢測。以提取的內(nèi)蒙古絨山羊DNA為模板，利用合成引物對目的片段進行PCR擴增。PCR反應(yīng)體系：PCR Master Mix 6.5 μl，上下游引物（10 ?μmol/L ）各0.3 μl，模板DNA （20 ?ng/μl ） 0.7 μl，以ddH? 2 O補至13.0 μl。PCR擴增程序為95 ℃預(yù)變性5 min;94 ℃變性30 s， 68～ 50 ℃（每個循環(huán)遞減1 ℃）退火30 s，72 ℃延伸18 s，共18個循環(huán);然后94 ℃變性30 s，50 ℃退火30 s，72 ℃延伸18 s，共34個循環(huán)，最后72 ℃延伸10 min。PCR產(chǎn)物用30 ?g/L 瓊脂糖凝膠電泳和測序方法鑒定分型? [26] 。引物由西安擎科澤西生物科技有限責(zé)任公司合成。

1.4 數(shù)據(jù)統(tǒng)計與分析

利用在線網(wǎng)站SHEsis（http：//analysis.bio-x.cn）分析? PRDM6? 基因InDel多態(tài)位點在群體對象中的Hardy-Weinberg平衡常數(shù)（ HWE ），分析計算基因型頻率、基因純合度（ H o）、雜合度（ H e）、有效等位基因數(shù)（ N e）和多態(tài)信息含量（ PIC ）等。應(yīng)用一般線性模型分析不同參數(shù)對性狀的影響： Y? ijk ?=μ+G? i?? ??+E? ij? ，這里 μ 為總體平均值， G? i? 代表基因型固定效應(yīng)， E? ij? 代表隨機誤差。利用SPSS 23.0軟件的單因素方差分析對此InDel突變位點不同基因型與內(nèi)蒙古絨山羊不同群體體質(zhì)量及體尺性狀進行關(guān)聯(lián)分析，同時分析內(nèi)蒙古絨山羊的體質(zhì)量和體尺性狀的皮爾遜（Pearson）相關(guān)系數(shù)。

2 結(jié)果與分析

2.1 PCR擴增與測序結(jié)果

利用合成引物對? rs656578433? 和? rs651603667? 位點進行PCR擴增，結(jié)果顯示，在內(nèi)蒙古絨山羊群體中，位于基因上游區(qū)域的? rs656578433?? InDel位點不存在多態(tài)性，位于第一內(nèi)含子區(qū)域的? rs651603667? ，即12 bp InDel突變（NC_030814.1：?? rs651603667? ， g：79985625-79985636delTTGACTGATCCA）存在多態(tài)性。電泳圖結(jié)果（圖1）顯示，? rs651603667? 存在3種基因型，純合插入型II（287 bp），雜合型ID（287 bp和275 bp）和純合缺失型DD（275 bp）。測序結(jié)果（圖2）顯示，與參考序列相比，? rs651603667? 突變型缺失了12 bp序列。綜上，? PRDM6? 基因的12 bp InDel位點在內(nèi)蒙古絨山羊群體中存在多態(tài)性，可進行進一步試驗及關(guān)聯(lián)分析。

2.2?? PRDM6? 基因突變位點的遺傳參數(shù)多樣性

依照電泳分型結(jié)果對內(nèi)蒙古絨山羊育成羊和成年羊試驗群體進行分析，結(jié)果（表1）顯示，在育成羊和成年羊群體中，? PRDM6? 基因12 bp InDel的基因型分布均符合Hardy-Weinberg平衡 （ P > 0.05），多態(tài)性信息含量（ PIC ）表明此突變位點屬中度多態(tài) （0.25

2.3 內(nèi)蒙古絨山羊? PRDM6? 基因多態(tài)性與生長性狀的關(guān)聯(lián)性

試驗探究了? PRDM6? 基因12 bp InDel突變與內(nèi)蒙古絨山羊育成羊和成年羊體質(zhì)量及生長性狀之間的關(guān)系。結(jié)果（表2）發(fā)現(xiàn)，在育成羊群體中，12 bp InDel突變與體質(zhì)量 （ P = 0.027）、胸深 （ P = 0.038）性狀顯著相關(guān)，II基因型對應(yīng)較好的性狀。

對成年羊? PRDM6? 基因12 bp InDel突變與體質(zhì)量及生長性狀進行關(guān)聯(lián)分析，發(fā)現(xiàn)12 bp InDel突變與體高 （ P = 0.003）、胸深 （ P = 0.004）極顯著相關(guān)，與體長 （ P = 0.049）、髖寬 （ P = 0.024）顯著相關(guān)。在體高性狀上，II基因型對應(yīng)較好的性狀;在體長、髖寬和胸深性狀上，ID基因型對應(yīng)較好的生長性狀（表3）。

2.4 內(nèi)蒙古絨山羊體質(zhì)量與生長性狀的皮爾遜相關(guān)性

將內(nèi)蒙古絨山羊的體質(zhì)量和生長性狀進行皮爾遜相關(guān)性分析。結(jié)果（表4）顯示，12 bp InDel突變位點顯著影響的生長性狀與體質(zhì)量均極顯著相關(guān) （ P < 0.01）。此外，胸圍、管圍、薦高、胸寬性狀也均與體質(zhì)量存在極顯著的正相關(guān) （ P < 0.01）。

3 討 論

內(nèi)蒙古絨山羊作為優(yōu)秀的肉絨兼用型山羊，其體質(zhì)量和生長性狀顯得尤為重要，一般來說，通常體質(zhì)量增大，胴體質(zhì)量和凈肉率隨之增大? [27] 。 ?PRDM6? 作為一種與體質(zhì)量及生長相關(guān)的候選基因，? PRDM6-CEP120? 基因位點已被報道與人類體質(zhì)量和肥胖癥相關(guān)? [22，28] ，但其與內(nèi)蒙古絨山羊生長性狀的關(guān)聯(lián)性尚未有報道。本研究首先驗證了位于 ?PRDM6? 基因第一內(nèi)含子12 bp InDel變異并探究其與內(nèi)蒙古絨山羊體質(zhì)量和生長性狀的關(guān)系，發(fā)現(xiàn)此位點與內(nèi)蒙古絨山羊的育成羊體質(zhì)量和胸深顯著相關(guān)，與成年羊體高、體長、髖寬和胸深顯著或極顯著相關(guān)。在育成羊群體中，此突變位點能顯著影響山羊的體質(zhì)量，但在成年羊群體中并未達到顯著水平，可能存在以下原因：（1）成年羊樣本量較少導(dǎo)致兩群體結(jié)果差異，后期可通過增加樣本量進一步分析;（2）伴隨動物機體的生長發(fā)育，基因可能存在時空表達，從而導(dǎo)致不同結(jié)果。例如，胰島素樣生長因子I（IGF-I）基因在鵝肌肉組織中早期表達量高于晚期表達量? [29] ;不同生長階段豬骨骼肌中 ?PRKAG3? 基因的表達存在差異? [30] ;肌肉生長抑制素（MSTN）基因在3月齡的新疆也木勒白羊肌肉中的表達水平明顯低于其他月齡? [31] ，因此，內(nèi)蒙古絨山羊不同生長階段可能伴隨著? PRDM6? 基因的差異表達。一般來說，較高的群體雜合度代表著更為豐富的遺傳多樣性? [32] ，在本研究中，成年羊群體雜合度高于育成羊群體，說明成年羊群體的遺傳變異更大，這可能是造成該基因突變位點與成年羊群體顯著相關(guān)的性狀數(shù)量高于育成羊群體的原因之一。Hardy-Weinberg平衡常數(shù)表明，在2個不同的內(nèi)蒙古絨山羊種群中，突變位點均處于Hardy-Weinberg平衡 （ P > 0.05）狀態(tài)，表明該位點受選擇壓力和基因突變的影響較小，突變位點的選擇強度可以適當(dāng)加強。通過皮爾遜相關(guān)系數(shù)分析結(jié)果顯示，2個不同群體的突變位點顯著影響的生長性狀均與體質(zhì)量存在極顯著關(guān)聯(lián)，說明? PRDM6? 基因可能通過影響多種生長性狀進而間接影響山羊的體質(zhì)量。

PRDM家族成員具有一個與甲基轉(zhuǎn)移酶相關(guān)的胞間PR結(jié)構(gòu)域和多個鋅指結(jié)構(gòu)，其中鋅指結(jié)構(gòu)與DNA、RNA和蛋白質(zhì)特異性結(jié)合，通過改變目的基因啟動子的狀態(tài)來控制基因表達，一些家族成員可以通過內(nèi)源性甲基轉(zhuǎn)移酶活性控制與細(xì)胞完整性相關(guān)的蛋白質(zhì)，影響原始生殖細(xì)胞發(fā)育、神經(jīng)系統(tǒng)發(fā)育和造血功能等? [33-36] 。研究結(jié)果表明， ?PRDM6? 在小鼠胚胎內(nèi)皮細(xì)胞中過量表達，誘導(dǎo)DNA合成在前期停滯與細(xì)胞凋亡，可減緩內(nèi)皮細(xì)胞的形成速率? [21] 。同時 ?PRDM6? 富含于血管前體細(xì)胞，降低內(nèi)皮細(xì)胞增殖速率，突變會導(dǎo)致? PRDM6? 基因的異常表達，從而影響血管平滑肌細(xì)胞早期分化生成及功能，而基因正常表達和細(xì)胞內(nèi)環(huán)境穩(wěn)定在胚胎發(fā)育過程中至關(guān)重要? [37-38] ，這表明 ?PRDM6? 基因可能以各種方式參與機體生長發(fā)育。有報道證實，? PRDM6? 作為參與組蛋白甲基化的一種甲基轉(zhuǎn)移酶，可能通過表觀遺傳修飾參與并調(diào)節(jié)機體生長發(fā)育? [39] 。綜上所述， ?PRDM6? 基因可能通過甲基轉(zhuǎn)移酶參與山羊組蛋白甲基化修飾或維持調(diào)控細(xì)胞增殖等途徑影響山羊的體質(zhì)量和生長性狀，具體作用機制仍需進一步探究。

目前，大量研究結(jié)果表明，基因多態(tài)性與家畜生長性狀密切相關(guān)? [40-42] 。例如， GHR 基因9 bp InDel突變與陜北白絨山羊體質(zhì)量和髖寬、體高等生長性狀顯著相關(guān)? [43] ，與內(nèi)蒙古絨山羊體質(zhì)量和體長、胸圍等生長性狀存在顯著關(guān)聯(lián)? [44] ; ?PLAG1? 基因19 bp InDel突變對多個國內(nèi)本土品種牛的生長性狀影響顯著? [45] ; ?MC4R? 基因的SNP突變對約克夏豬的體質(zhì)量和日增質(zhì)量影響顯著，? CDC16? 基因第18外顯子上SNP突變位點與約克夏豬的體質(zhì)量、體高和體長顯著相關(guān)? [46] ，上述研究結(jié)果表明，基因多態(tài)性可能影響家畜的生長性狀。在本研究中，我們發(fā)現(xiàn) ?PRDM6? 基因在內(nèi)蒙古絨山羊種群中存在一個12 bp的InDel變異，并發(fā)現(xiàn)它與內(nèi)蒙古絨山羊的體質(zhì)量和多個生長性狀有顯著或極顯著關(guān)聯(lián)，可作為內(nèi)蒙古絨山羊育種中的有效候選分子標(biāo)記。

參考文獻：

[1]? 劉志紅，謝遇春，奈日樂，等. 內(nèi)蒙古白絨山羊不同部位肌肉脂肪酸的研究[J]. 中國畜牧雜志， 2020， 56（5）： 88-91.

[2] 烏力吉，成海榮，哈斯牧仁，等. 內(nèi)蒙古阿爾巴斯白絨山羊羔羊短期育肥效果與屠宰性能分析[J]. 畜牧與飼料科學(xué)， 2021， 42（3）： 21-27.

[3] 孟偉星，李玉榮，劉 斌，等. 內(nèi)蒙古白絨山羊經(jīng)濟性狀與候選基因分析[J]. 畜牧與飼料科學(xué)， 2008（4）： 35-38.

[4] 吳鐵梅. 不同飼養(yǎng)模式對絨山羊羔羊育肥性能、屠宰性能及肉品質(zhì)的影響[D]. 呼和浩特： 內(nèi)蒙古農(nóng)業(yè)大學(xué)， 2013.

[5] MA L， QIN Q， YANG Q， et al. Associations of six SNPs of?? POU1F1-PROP1-PITX1-SIX3?? pathway genes with growth traits in two Chinese indigenous goat breeds[J]. Annals of Animal Science， 2017， 17（2）： 399-411.

[6] 王振宇. 內(nèi)蒙古阿爾巴斯絨山羊與大青山山羊體尺與體重的差異比較分析[D]. 呼和浩特： 內(nèi)蒙古農(nóng)業(yè)大學(xué)， 2018.

[7] WANG X Y， YANG Q， WANG K， et al. Two strongly linked single nucleotide polymorphisms （Q320P and V397I） in ?GDF9 ?gene are associated with litter size in cashmere goats[J]. Theriogenology， 2019， 125： 115-121.

[8] YANG Q， YAN H L， LI J， et al. A novel 14-bp duplicated deletion within goat ?GHR ?gene is significantly associated with growth traits and litter size[J]. Animal Genetics， 2017， 48（4）： 499-500.

[9] KNORST V， BYRNE S， YATES S， et al. Pooled DNA sequencing to identify SNPs associated with a major QTL for bacterial wilt resistance in Italian ryegrass （ Lolium multiflorum ?Lam.）[J]. Theoretical and Applied Genetics， 2019， 132（4）： 947-958.

[10] 譚雙香. ?PRDM 基因啟動子甲基化在非小細(xì)胞肺癌中的作用研究[D]. 長沙： 中南大學(xué)， 2014.

[11] FUMASONI I， MEANI N， RAMBALDI D， et al. Family expansion and gene rearrangements contributed to the functional specialization of ?PRDM ?genes in vertebrates[J]. BMC Evolutionary Biology， 2007， 7（1）： 187-198.

[12] LI N， SUBRAHMANYAN L， SMITH E， et al. Mutations in the histone modifier?? PRDM6?? are associated with isolated nonsyndromic patent ductus arteriosus[J]. The American Journal of Human Genetics， 2016， 98（6）： 1082-1091.

[13] 李曉波，張俊武. 真核生物中鋅指蛋白的結(jié)構(gòu)與功能[J]. 中國生物化學(xué)與分子生物學(xué)報， 2009， 25（3）： 206-211.

[14] DAVIS C， HABERLAND M， AMOLD M， et al.?? PRISM/PRDM6? ， a transcriptional repressor that promotes the proliferative gene program in smooth muscle cells[J]. Molecular and Cellular Biology， 2006，26（7）： 2626-2636.

[15] HUANG S， SHAO G， LIU L. The PR domain of the Rb-binding zinc finger protein RIZ1 is a protein binding interface and is related to the SET domain functioning in chromatin-mediated gene expression.[J]. The Journal of Biological Chemistry， 1998， 273（26）： 15933-15939.

[16] HOHENAUER T， MOORE A W. The Prdm family： expanding roles in stem cells and development[J]. Development （Cambridge）， 2012， 139（13）： 2267-2282.

[17] 趙生軍，張 勇，閻 萍，等. PRDM家族蛋白結(jié)構(gòu)與功能研究進展[J]. 中國畜牧獸醫(yī)， 2016， 43（5）： 1188-1193.

[18] SLIM M， YING X T， DIANA L. The role of PRDMs in cancer： one family， two sides[J]. Current Opinion in Genetics & Development， 2016， 36： 83-91.

[19] GEWIES A， CASTINEIRAS-VILARION M， FERCH U， et al.?? PRDM6?? is essential for cardiovascular development ?in vivo [J]. PLoS One， 2013， 8（11）： e81833.

[20] HADZISELIMOVIC F， CATHOMAS G， VERKAUSKAS G， et al. PRDM histone methyltransferase mRNA levels increase in response to curative hormone treatment for cryptorchidism-dependent male infertility[J]. Genes， 2018， 9（8）： 391-403.

[21] WU Y X， FERGUSON J E， WANG H， et al.?? PRDM6?? is enriched in vascular precursors during development and inhibits endothelial cell proliferation， survival， and differentiation[J]. Journal of Molecular and Cellular Cardiology， 2008， 44（1）： 47-58.

[22] HU Y， TAN L J， CHEN X D， et al. Identification of novel potentially pleiotropic variants associated with osteoporosis and obesity using the cFDR method[J]. The Journal of Clinical Endocrinology & Metabolism， 2017， 103（1）： 125-138.

[23] WANG Z， WANG C L， GUO Y N， et al. Screening of deletion variants within the goat?? PRDM6?? gene and its effects on growth traits[J]. Animals， 2020， 10（2）： 208-219.

[24] ALJANABI S， MARTINEZ I. Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques[J]. Nucleic Acids Research， 1997， 25（22）： 4692-4693.

[25] ZHU H J， ZHANG Y H， BAI Y Y， et al. Relationship between SNPs of?? POU1F1?? gene and litter size and growth traits in Shaanbei white cashmere goats[J]. Animals， 2019， 9（3）： 114-125.

[26] 楊 韓，張陽海，潘傳英，等. 反芻動物? POU1F1? 基因多態(tài)性與經(jīng)濟性狀關(guān)聯(lián)研究進展[J]. 中國牛業(yè)科學(xué)， 2019， 45（2）： 42-48.

[27] MA L， QIN Q M， YANG Q， et al. Associations of six SNPs of?? POU1F1-PROP1-PITX1-SIX3?? pathway genes with growth traits in two Chinese indigenous goat breeds[J]. Annals of Animal Science， 2017， 17（2）： 399-411.

[28] RIVEROS-MCKAY F， MISTRY V， BOUNDS R， et al. Genetic architecture of human thinness compared to severe obesity[J]. PLoS Genetics， 2019， 15（1）： e1007603.

[29] 史旭升，孫永峰，楊童奧，等. 胰島素樣生長因子Ⅰ基因在鵝不同發(fā)育時期肌肉組織中的表達[J]. 吉林農(nóng)業(yè)大學(xué)學(xué)報， 2011， 33（2）： 214-217.

[30] 李夢云，鄭 萍，李婉濤，等.?? PRKAG3? 基因在不同品種豬不同生長階段骨骼肌中的表達差異及其表達量與肉質(zhì)的關(guān)系[J]. 動物營養(yǎng)學(xué)報， 2017， 29（5）： 1661-1669.

[31] 祖玲玲，吐來力江·哈木太，安外爾·熱合曼，等. 新疆也木勒白羊肌肉生長抑制素基因在不同生長時期的表達研究[J]. 中國畜牧獸醫(yī)， 2014， 41（12）： 98-101.

[32] 馬曉萌，軒俊麗，王慧華，等.?? TXNRD1? 基因多態(tài)性與烏珠穆沁綿羊生長性狀的關(guān)聯(lián)分析[J]. 畜牧獸醫(yī)學(xué)報， 2016， 47（5）： 909-921.

[33] ED Z， ROSA C D， ABBONDANZA C， et al. PRDM proteins： molecular mechanism in signal transduction and transcriptional regulation[J]. Biology， 2013， 2（1）： 107-141.

[34] HAYASHI K， YOSHIDA K， MATSUI Y. A histone H3 methyltransferase controls epigenetic events required for meiotic prophase[J]. Nature， 2005， 438（7066）： 374-378.

[35] DENG Q D， HUANG S.?? PRDM5?? is silenced in human cancers and has growth suppressive activities[J]. Oncogene， 2004， 23（28）： 4903-4910.

[36] FOG C K， GALLI G G， LUND A H. PRDM proteins： important players in differentiation and disease[J]. BioEssays， 2015， 34（1）： 50-60.

[37] MEIJER I A， SIMOES-LOPES A A， LAURENT S， et al. A novel duplication confirms the involvement of 5q23.2 in autosomal dominant leukodystrophy[J]. Archives of Neurology， 2008， 65（11）： 1496-1501.

[38] PHILIP R， TOWBIN J A， SATHANANDAM S， et al. Effect of patent ductus arteriosus on the heart in preterm infants[J]. Congenital Heart Disease， 2019， 14（1）： 33-36.

[39] FELLOUS A， EARLEY R L， SILVESTRE F. The ?Kdm/Kmt ?gene families in the self-fertilizing mangrove rivulus fish， Kryptolebias marmoratus， suggest involvement of histone methylation machinery in development and reproduction[J]. Gene， 2019， 687： 173-187.

[40] 張 俊，李隱俠，錢 勇，等. 山羊? MC4R? 基因多態(tài)性與體質(zhì)量性狀的關(guān)聯(lián)性分析[J]. 江蘇農(nóng)業(yè)科學(xué)，2019，47（19）：53-56.

[41] 鄒 輝，瞿秋紅，夏 琴，等. 努比亞山羊? BMPR-IB? 基因多態(tài)性與其產(chǎn)羔性狀的關(guān)聯(lián)分析[J].南方農(nóng)業(yè)學(xué)報，2019，50（4）：860-866.

[42] 梅秀麗，李叢艷，鄺良德，等. 胰島素樣生長因子1受體多態(tài)性與生產(chǎn)性能關(guān)系的研究進展[J]. 江蘇農(nóng)業(yè)科學(xué)，2019，47（2）：29-32.

[43] 閆海龍，白洋洋，王 真，等. ?GHR 基因9-bp InDel與陜北白絨山羊體重和生長性狀的關(guān)聯(lián)研究[J]. 家畜生態(tài)學(xué)報， 2019， 40（2）： 27-31.

[44] 王 真，王 敏，李 銘，等. 內(nèi)蒙古白絨山羊 GHR 和 ?GDF9? 基因多態(tài)性及其與生產(chǎn)性狀的相關(guān)性分析[J]. 西北農(nóng)林科技大學(xué)學(xué)報（自然科學(xué)版）， 2020， 48（3）： 1-8.

[45] XU W， HE H， ZHENG L， et al. Detection of 19-bp deletion within?? PLAG1 ??gene and its effect on growth traits in cattle[J]. Gene， 2018， 675： 144-149.

[46] 趙聰哲，羅曉彤，李兆華，等. 大白豬? MC4R、CDC16? 基因多態(tài)性及其與生長性狀的關(guān)聯(lián)分析[J]. 中國畜牧獸醫(yī)， 2019， 46（3）： 792-799.

（責(zé)任編輯：陳海霞）