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

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

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

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

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

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

1 材料與方法

1.1 試驗(yàn)材料

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

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

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

1.3?? PRDM6? 基因的PCR擴(kuò)增

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

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

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

2 結(jié)果與分析

2.1 PCR擴(kuò)增與測(cè)序結(jié)果

利用合成引物對(duì)? rs656578433? 和? rs651603667? 位點(diǎn)進(jìn)行PCR擴(kuò)增，結(jié)果顯示，在內(nèi)蒙古絨山羊群體中，位于基因上游區(qū)域的? rs656578433?? InDel位點(diǎn)不存在多態(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）。測(cè)序結(jié)果（圖2）顯示，與參考序列相比，? rs651603667? 突變型缺失了12 bp序列。綜上，? PRDM6? 基因的12 bp InDel位點(diǎn)在內(nèi)蒙古絨山羊群體中存在多態(tài)性，可進(jìn)行進(jìn)一步試驗(yàn)及關(guān)聯(lián)分析。

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

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

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

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

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

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

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

3 討 論

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

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

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

參考文獻(xiàn)：

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

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

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

[4] 吳鐵梅. 不同飼養(yǎng)模式對(duì)絨山羊羔羊育肥性能、屠宰性能及肉品質(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 基因啟動(dòng)子甲基化在非小細(xì)胞肺癌中的作用研究[D]. 長(zhǎng)沙： 中南大學(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]. 中國(guó)生物化學(xué)與分子生物學(xué)報(bào)， 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ìn)展[J]. 中國(guó)畜牧獸醫(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] 楊 韓，張陽(yáng)海，潘傳英，等. 反芻動(dòng)物? POU1F1? 基因多態(tài)性與經(jīng)濟(jì)性狀關(guān)聯(lián)研究進(jìn)展[J]. 中國(guó)牛業(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] 史旭升，孫永峰，楊童奧，等. 胰島素樣生長(zhǎng)因子Ⅰ基因在鵝不同發(fā)育時(shí)期肌肉組織中的表達(dá)[J]. 吉林農(nóng)業(yè)大學(xué)學(xué)報(bào)， 2011， 33（2）： 214-217.

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

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

[32] 馬曉萌，軒俊麗，王慧華，等.?? TXNRD1? 基因多態(tài)性與烏珠穆沁綿羊生長(zhǎng)性狀的關(guān)聯(lián)分析[J]. 畜牧獸醫(yī)學(xué)報(bào)， 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] 張 俊，李隱俠，錢(qián) 勇，等. 山羊? 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é)報(bào)，2019，50（4）：860-866.

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

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

[44] 王 真，王 敏，李 銘，等. 內(nèi)蒙古白絨山羊 GHR 和 ?GDF9? 基因多態(tài)性及其與生產(chǎn)性狀的相關(guān)性分析[J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào)（自然科學(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)性及其與生長(zhǎng)性狀的關(guān)聯(lián)分析[J]. 中國(guó)畜牧獸醫(yī)， 2019， 46（3）： 792-799.

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