佘平昌,吳曉云,梁春年,禇 敏,丁學智,閻 萍*
(1.中國農業科學院蘭州畜牧與獸藥研究所,蘭州 730050; 2.甘肅省牦牛繁育工程重點實驗室,蘭州 730050)
?
牦牛角性狀候選基因的篩選
佘平昌1,2,吳曉云1,2,梁春年1,2,禇敏1,2,丁學智1,2,閻萍1,2*
(1.中國農業科學院蘭州畜牧與獸藥研究所,蘭州 730050; 2.甘肅省牦牛繁育工程重點實驗室,蘭州 730050)
本研究旨在篩選牦牛角性狀候選基因,為角發育機制的研究和無角牦牛的培育提供科學依據。采集牦牛角發育初期胎齡90 d的有角和無角牛的角基及額部皮膚組織,利用實時熒光定量PCR技術(qRT-PCR),比較牛角性狀候選基因OLIG1、OLIG2、IFNAR1、IFNAR2、C1H21orf62、GCFC1、IFNGR2、SYNJ1、IL10RB、GART、RXFP2、FOXL2、TWIST1、TWIST2、ZEB2、E-cadherin和P-cadherinmRNA在角基間及皮膚間的相對表達量。結果顯示,候選基因OLIG1、IFNAR1、C1H21orf62、SYNJ1、IL10RB、GART、RXFP2、TWIST2和E-cadherin在有角和無角牦牛角基間和皮膚間轉錄量差異不顯著(P>0.05)。IFNGR2、FOXL2、TWIST1和P-cadherin在有角和無角牦牛角基間轉錄量差異顯著(P<0.05)。IFNAR2和GCFC1在有角和無角牦牛角基間和皮膚間轉錄量差異都顯著(P<0.05)。OLIG2和ZEB2在有角和無角牦牛角基間轉錄量差異極顯著(P<0.01)。推測OLIG2和FOXL2可能是牛角性狀的關鍵候選基因,同時ZEB2、TWIST1和IFNGR2基因可能對牛角的形成具有重要作用,而P-cadherin基因可能參與了牛角的發育。
牛角;無角性狀;牦牛;候選基因;qRT-PCR
牛角是??苿游锾赜械慕M織結構。在男耕女織的古代原始生產過程中,牛角是進行耕地過程中牽牛的基礎,所以牛角在原始耕作時代中是一種生產力。但是進入集約規?;曫B后,牛群常發生爭斗及傷害飼養人員的行為,造成了動物的傷害及養殖場的混亂,牛角變成了一種制約家畜集約化和設施養殖的性狀。所以,很多養牛場將犢牛進行去角[1],然而去角增加了養殖成本,同時容易引起犢牛的應激反應,是違背動物福利的行為。動物福利問題已經在發達國家引起高度重視,歐盟啟動了關于禁止動物閹割和去角的ALCASDE工程[2],而德國的動物福利組織也聯合政府發表申明將積極推進無角公牛精液的使用。因此,無角牛品種的培育是未來牛育種的主要方向之一。
牛角根據生物學分類屬于洞角,由角質鞘和骨心組成。牛角的發育較早,60胎齡的胎牛就細微可見角的形成,90胎齡時發育更加明顯,基本可以通過肉眼判定角的有無,在角的生長部位可見真皮乳頭狀的突起形成角竇[3],同時出現表皮增厚的現象[4],角竇和皮膚間通過真皮乳頭進行信號傳導[5]。但是只有在出生時才可見螺紋狀角的出現[6],出生后6月齡才開始出現角的空腔結構[3]。組織學觀察發現胚胎期牛角與皮膚間的差異主要有4點:腺導管真皮細胞的集群不同、囊泡化的角質細胞、沒有毛囊的發育[7]、真皮中出現大量的神經束[8],推測頭皮層的神經纖維對牛角形成具有觸動作用,所以神經發育相關基因可能會在角基部位高表達。在牛角性狀候選基因的研究中,根據目前的研究進展,很多研究將無角突變定位在1號染色體(BTA1,Bostaurusautosome1)1.265~2.198 M上1 M左右的一段序列內,包括C1H21orf62、OLIG1、OLIG2、SYNJ1、GCFC4、GART、IL10RB、IFNAR1、IFNAR2、IFNGR2等基因[9-14]。而不在該序列內的一些基因也可能對牛角形成和發育具有重要作用,FOXL2和RXFP2基因在牛角形成和發育過程中高度表達[7,15],而ZEB2、TWIST1和TWIST2是畸角牛的候選基因[16-18],E-cadherin和P-cadherin對胚胎分化及表皮、毛囊、神經等的形成具有重要作用[19-21]。牛角是一種質量性狀(Qualitative trait),其遺傳方式可能為多基因遺傳(Polygenic trait)。目前對牛角性狀相關基因還沒有進行精確定位,還需進一步研究。本研究通過熒光定量PCR (Real-time quantitative PCR,qRT-PCR)方法分析牛角性狀候選基因mRNA的表達水平,初步揭示牛角發育和形成的分子機制,以期為無角牦牛品種的培育提供科學依據和理論支撐。
在青海省大通牦牛種牛場采集90 d胎齡的有角和無角牦牛的角基部及額部皮膚樣品,液氮保存,并對其角基部及全身進行拍照,測量其頂臀長,根據胎兒頂臀長與胎齡的線性關系推算胎齡[22]。
1.1RNA提取及cDNA合成
采用Trizol (Invitrogen)法提取有角和無角牦牛的角基及額部皮膚組織的總RNA,按照TaKaRa公司試劑盒Prime Script RT reagent Kit with gDNA Eraser(Perfect Real Time)說明書合成cDNA。
1.2熒光定量PCR
使用Premier 5.0設計候選基因mRNA熒光定量PCR引物(表1),選擇GAPDH作為內參基因。引物由上海英駿生物公司合成,產物經2%瓊脂糖凝膠電泳檢測。
定量PCR(CFX96)反應使用25 μL體系:SYBR?Premix Ex TaqTMⅡ (Tli RNaseH Plus) 12.5 μL,上下游引物各1 μL,cDNA 2 μL,ddH2O 8.5 μL。熒光定量 PCR 的反應程序:95 ℃預變性30 s;95 ℃變性5 s,53.5 ℃退火30 s,40個循環;之后進行熔解曲線分析,反應程序:65 ℃ 逐漸升溫至95 ℃,每升溫0.5 ℃讀板1次。每個組織4個個體重復,每個樣本3次試驗重復。反應結束后以熔解曲線來判定qRT-PCR反應的特異性;獲得每個樣品的Ct值,運用2-ΔΔCt法計算相對轉錄水平。
1.3圖像數據處理
本試驗采用Excel 2013對數據進行統計,SPSS 17.0進行one-way ANOVA分析,Origin 7.5作柱狀圖。
2.1BTA1上候選區域內基因的相對表達量分析
為了分析候選基因的轉錄情況,本研究對比候選基因在有角牛和無角牛的角基間及額部皮膚間的相對轉錄水平,據此分析該基因是否參與了角的形成和發育(圖1)。OLIG1、OLIG2、IFNAR1、IFNAR2、C1H21orf62、GCFC1、IFNGR2、SYNJ1、IL10RB和GART是1號染色體上1 M左右片段上的基因。研究發現,OLIG1、IFNAR1、C1H21orf62、SYNJ1、IL10RB、GART、PXFP2、TWIST2和E-cadherin在有角和無角牦牛的角基間及皮膚間的轉錄量差異都不顯著(P>0.05);OLIG2在有角牛角基的轉錄量極顯著高于無角牛角基(P<0.01),而在皮膚間差異不顯著;IFNAR2在有角牛的角基和皮膚的轉錄量都要顯著高于無角牛(P<0.05);相反,GCFC1在有角牛的角基和皮膚的轉錄量都要顯著低于無角牛(P<0.05)。IFNGR2在無角牛角基中的轉錄量顯著高于有角牛(P<0.05),而在皮膚間的轉錄差異不顯著(P>0.05)。
表1熒光定量PCR引物序列
Table 1Primer sequences used for the real-time quantitative PCR

基因名稱Genename上游引物序列(5'-3')Forwardprimer下游引物序列(5'-3')ReverseprimerOLIG1GTGGCGATCTTGGAGAGCGAGAGGAAGCGGATGCACOLIG2CCGGCATCTACCAGTACATTATAGCGCACTCGTTGAGGCTGAGGTIFNAR1TCTTGACAGACTCACTGCTGCTGAAGCAAATCTGAAGCCTGAAIFNAR2CTCAAGGGCTCTCACACACAACTTCCGAGGAGGACAGTGAC1H21orf62ACCAGCTGTCAAGCCTGAGTCAGAAGGAGGTGCCATTCTCGCFC1GATGATGAGAAACGCCGCATAGGCATTGTCTGGTAGGTGTIFNGR2CTTGGACCAGCACTTTGGTTCTCGCCAAAGGATGACACTTSYNJ1TGGTACGGTGTTGGTCTCAAATTCAAGGCAGAGCTTCCCTIL10RBTCTGTGATGACGCTCAGTTTGTTTATGAGCCACCCTCATCACGARTGCCCGAGTACTTGTCATTGGCCTGGGGTAACCAACACTTGRXFP2AGAACCCCACAATCCAGATGTGACCCTGGAGAAGTTCCTGFOXL2CCGGCATCTACCAGTACATTATAGCGCACTCGTTGAGGCTGAGGTTWIST1TGTGAGTCAGTTTGATCCCAATGGCATCATTATGGACTTTCTCCTWIST2TTCAACTGGACCAAGGCTCTCCAGGCTAACTTGAGACAGTCAZEB2ACACGGGTTCTGAAACTGATGAGGCACGCTAGCTGGACTTCTE-cadherinGTACACCTTCATCGTCCAGAGCTAAGCTCTTCAATGGCTTGTCCATTTGAP-cadherinGCAGCCTGATACGGTAGAGCAGATCGAACCGGGTACTGGGAPDHATGCAACCTGGTGATAAAAGAGGAGGGCGCTAAGAGAAGAGTCAT
2.2其他重要候選基因的相對表達量分析
RXFP2、FOXL2、TWIST1、TWIST2和ZEB2也是重要的角性狀候選基因。研究發現,RXFP2和TWIST2在有角和無角牦牛的角基間和皮膚間轉錄量差異都不顯著(P>0.05);FOXL2在有角牛角基的轉錄量顯著高于無角牛(P<0.05),而在皮膚間差異不顯著;TWIST1卻在無角牛角基的轉錄量顯著高于有角牛(P<0.05),皮膚間的轉錄量差異不顯著。ZEB2在有角牛角基的轉錄量極顯著低于無角牛(P<0.01),而在皮膚間同樣差異不顯著(P>0.05)。
2.3角發育相關基因的相對表達量分析
關于E-cadherin和P-cadherin與角性狀的相關研究報道較少,但E-cadherin是表皮及毛囊的發育關鍵基因,P-cadherin的表達對表皮發育及神經分化具有重要作用,而據研究進展可知,角的發育是角基表皮分化的過程,中間也伴隨著毛囊的發育,引入分析為試驗提供參考。研究發現,E-cadherin在角基間和皮膚間的轉錄差異都不顯著(P>0.05),而P-cadherin在有角牛角基中的轉錄量顯著高于無角牛(P<0.05),皮膚間的轉錄量差異不顯著(P>0.05)。


horn.角基組織;skin.皮膚組織;horned.有角牛群體;polled.無角牛群體;*.P<0.05;**.P<0.01horn.Horn bud;skin.Forehead skin;horned.Horned yaks;polled.Polled yaks;*.P<0.05;**.P<0.01圖1 候選基因在有角和無角牦牛的角基間和皮膚間的相對轉錄水平Fig.1 Relative expression analyses of candidate genes in horn bud and skin tissues between horned and polled yaks by qRT-PCR
本研究對BAT1上定位的無角性狀基因座區域內的候選基因的表達量進行了分析。研究結果顯示,OLIG1、IFNAR1、C1H21orf62、SYNJ1、IL10RB和GART在有角和無角牦牛的角基間及皮膚間的轉錄量差異不顯著(P>0.05);OLIG2、IFNAR2、GCFC1和IFNGR2的mRNA轉錄量在角基間具有顯著(P<0.05)或極顯著(P<0.01)差異。候選基因中,OLIG2為少突膠質細胞轉錄因子,其主要參與少突膠質細胞的形成分化及髓鞘的形成,同時對哺乳動物的組織分化起重要調節作用[23]。本研究發現,OLIG2在角基間轉錄差異極顯著,研究結果與在荷斯坦牛[7]和西門塔爾牛[15]中的研究結果一致。近期研究發現,OLIG2參與神經分化,推測該基因對牛角形成的作用機制可能與神經調控有關,說明OLIG2對角的發育有重要作用,因此將該基因作為荷斯坦奶牛、西門塔爾牛及牦牛無角性狀的候選基因[7,13,15]。IFNAR2和GCFC1基因在角基間和皮膚間轉錄差異均顯著,說明IFNAR2和GCFC1 的表達模式不具有角組織特異性,IFNAR2和GCFC1可能對表皮的發育過程具有關鍵作用。IFNGR2是參與T細胞活動的一個細胞因子,在無角荷斯坦牛中該基因存在突變位點AC000158:g.1390292G>A[24],但是S.Rothammer等[12]的研究又否定了該突變與無角性狀的相關性,本研究發現,有角牛中IFNGR2基因的表達量顯著低于無角牛,所以還不能排除IFNGR2可能為角性狀候選基因。
本研究發現,其他候選基因RXFP2、TWIST2和E-cadherin在有角和無角牦牛的角基間和皮膚間轉錄量差異不顯著(P>0.05);FOXL2、TWIST1、ZEB2和P-cadherin的mRNA轉錄量在角基間具有顯著(P<0.05)或極顯著(P<0.01)差異。候選基因中,FOXL2又名叉頭框因子2,是叉頭框轉錄因子家族成員之一,對哺乳動物早期胚胎的發育和生殖系統功能維持具有重要作用[25],本研究發現,FOXL2在角基間的表達量存在顯著差異,在西門塔爾牛有角牛角基中該基因的表達量也顯著高于無角牛[15],同時FOXL2基因是羊角性狀的主效基因之一[26],所以FOXL2可能參與了牛角的形成和發育。TWIST1[16]和ZEB2[18]都是目前研究畸角牛的關鍵基因,TWIST1是在胚胎發育的一個轉錄因子,對誘導細胞移動及組織朔形具有重要作用,參與骨骼、肌肉和神經等的形成。ZEB2屬于E盒結合鋅指蛋白家族,在胚胎發育過程中對神經系統的發育具有重要作用。而牦牛中沒有畸角牛,但是本研究發現其在角基中表達量差異顯著,推測TWIST1和ZEB2可能與牛角的形成或發育具有調控作用,但是其功能仍需進一步明確。P-cadherin基因為鈣粘蛋白家族成員,在胚胎發育過程中對表皮及神經的發育發揮著重要的作用,而角的發育為表皮增厚的過程,本研究中P-cadherin在角基間轉錄差異顯著,推測該基因可能參與了牛角的發育。
本研究中,候選基因OLIG1、IFNAR1、C1H21orf62、SYNJ1、GART和RXFP2雖然在有角和無角牦牛的角基間和皮膚間的轉錄均不存在顯著差異(P>0.05),但是根據本研究的內容不能完全排除這些基因可能參與角的形成或發育的可能。而IL10RB、TWIST2和E-cadherin基因的相關研究較少,且本試驗結果顯示這些基因在有角和無角牛角基及皮膚間的表達量都不存在差異,故推測IL10RB、TWIST2和E-cadherin基因可能不是牛角早期形成和發育的關鍵基因。
牦牛角形狀候選基因的表達研究發現,候選基因OLIG1、IFNAR1、C1H21orf62、SYNJ1、IL10RB、GART、RXFP2、TWIST2和E-cadherin在有角和無角牦牛的角基和皮膚間相對轉錄水平差異不顯著(P>0.05)。IFNGR2、FOXL2、TWIST1和P-cadherin基因在角基間轉錄量差異顯著(P<0.05)。IFNAR2和GCFC1基因在角基和皮膚間轉錄量差異都顯著(P<0.05)。OLIG2和ZEB2基因在角基間轉錄量差異極顯著(P<0.01)。結合前期研究推測OLIG2和FOXL2基因可能是牛角性狀的關鍵候選基因,同時ZEB2、TWIST1和IFNGR2基因可能對牛角的形成具有重要作用,而P-cadherin基因可能參與了牛角的發育。
[1]GRAF B,SENN M.Behavioural and physiological responses of calves to dehorning by heat cauterization with and without local anesthesia[J].ApplAnimBehavSci,1999,62(2):153-171.
[2]OLIVER M A,MIRABITO L,VEISSIER I,et al.Study on the improved methods for animal-friendly production in particular on alternatives to the castration of pigs and on alternatives to the dehorning of cattle[D].England:University of the West of England,2009:1-65.
[3]DOVE W F.The physiology of horn growth:A study of the morphogenesis,the interaction of tissues,and the evolutionary processes of a mendelian recessive character by means of transplantation of tissues[J].JExpZool,2005,69(3):347-450.
[4]任強,郁楓,王旭強,等.牦牛角的形態發生及Shh和β-catenin的時空表達[J].西北農業學報,2009,18(1):24-29.
REN Q,YU F,WANG X Q,et al.Temporal-spatial expression of Shh and β-catenin during horn morphogenesis in yak[J].ActaAgriculturaeBoreali-occidentalisSinica,2009,18(1):24-29.(in Chinese)
[5]HILZHEIMER M,LENGERKEN H.Natürliche rassengeschichte der hauss?ugetiere[M].Berlin:De Gruyter,1926:25-183.
[6]RUSSE I,SINOWATZ F.Lehrbuch der Embryologie der Haustiere[M].Berlin:Paul Parey,1991:304-344.
[7]ALLAIS-BONNET A,GROHS C,MEDUGORAC I,et al.Novel insights into the bovine polled phenotype and horn ontogenesis in Bovidae[J].PLoSOne,2013,8(5):e63512.
[8]WIENER D J,WIEDEMAR N,WELLE M M,et al.Novel features of the prenatal horn bud development in cattle(Bostaurus)[J].PLoSOne,2015,10(5):e0127691.
[9]GEORGES M,DRINKWATER R,KING T,et al.Microsatellite mapping of a gene affecting horn development inBostaurus[J].NatGenet,1993,4(2):206-210.
[10]SEICHTER D,RUSS I,ROTHAMMER S,et al.SNP-based association mapping of the polled gene in divergent cattle breeds[J].AnimGenet,2012,43(5):595-598.
[11]MEDUGORAC I,SEICHTER D,GRAF A,et al.Bovine polledness--an autosomal dominant trait with allelic heterogeneity[J].PLoSOne, 2012,7(6):e39477.[12]ROTHAMMER S,CAPITAN A,MULLAART E,et al.The 80-kb DNA duplication on BTA1 is the only remaining candidate mutation for the polled phenotype of Friesian origin[J].GenetSelEvol,2014,46:44.
[13]LIU W B,LIU J,LIANG C N,et al.Associations of single nucleotide polymorphisms in candidate genes with the polled trait in Datong domestic yaks[J].AnimGenet,2014,45(1):138-141.
[14]趙娟花,梁春年,裴杰,等.牦牛C1H21orf62 基因的克隆、生物信息學及表達分析[J].廣東農業科學,2015(5):98-102.
ZHAO J H,LIANG C N,PEI J,et al.Cloning,bioinformatics and expression analysis ofC1H21orf62 gene in yak[J].GuangdongAgriculturalSciences,2015(5):98-102.(in Chinese)
[15]WIEDEMAR N,TETENS J,JAGANNATHAN V,et al.Independent polled mutations leading to complex gene expression differences in cattle[J].PLoSOne,2014,9(3):e93435.
[16]CAPITAN A,GROHS C,GAUTIER M,et al.The scurs inheritance:new insights from the French Charolais breed[J].BMCGenet,2009,10:33.
[17]ASAI M,BERRYERE T G,SCHMUTZ S M.The scurs locus in cattle maps to bovine chromosome 19[J].AnimGenet,2004,35(1):34-39.
[18]CAPITAN A,ALLAIS-BONNET A,PINTON A,et al.A 3.7 Mb deletion encompassing ZEB2 causes a novel polled and multisystemic syndrome in the progeny of a somatic mosaic bull[J].PLoSOne,2012,7(11):e49084.
[19]LARUE L,OHSUGI M,HIRCHENHAIN J,et al.E-cadherin null mutant embryos fail to form a trophectoderm epithelium[J].ProcNatlAcadSciUSA,1994,91(17):8263-8267.
[20]MCMILLEN P,HOLLEY S A.Integration of cell-cell and cell-ECM adhesion in vertebrate morphogenesis[J].CurrOpinCellBiol,2015,36:48-53.
[21]吳興龍,王棉娟,胡春超,等.E-鈣粘蛋白抗體對小鼠早期胚胎發育的影響[J].畜牧獸醫學報,2012,43(8):1222-1229.
WU X L,WANG M J,HU C C,et al.Effects of E-cadherin antibody on mouse early embryo development [J].ActaVeterinariaetZootechnicaSinica,2012,43(8):1222-1229.(in Chinese)
[22]LIU B,CUI Y,YANG B,et al.Morphometric analysis of yak placentomes during gestation[J].AnatRec(Hoboken),2010,293(11):1873-1879.
[23]ZHOU Q,ANDERSON D J.The bHLH transcription factors OLIG2 and OLIG1 couple neuronal and glial subtype specification[J].Cell,2002,109(1):61-73.
[24]GLATZER S,MERTEN N J,DIERKS C,et al.A single nucleotide polymorphism within the interferon gamma receptor 2 gene perfectly coincides with polledness in Holstein cattle[J].PLoSOne,2013,8(6):e67992.
[25]BATISTA F,VAIMAN D,DAUSSET J,et al.Potential targets of FOXL2,a transcription factor involved in craniofacial and follicular development,identified by transcriptomics[J].ProcNatlAcadSciUSA,2007,104(9):3330-3335.
[26]PAILHOUX E,VIGIER B,CHAFFAUX S,et al.A 11.7-kb deletion triggers intersexuality and polledness in goats[J].NatGenet,2001,29(4):453-458.
(編輯郭云雁)
Screening of Candidate Genes for Horn Trait of Yak
SHE Ping-chang1,2,WU Xiao-yun1,2,LIANG Chun-nian1,2,CHU Min1,2,DING Xue-zhi1,2,YAN Ping1,2*
(1.LanzhouInstituteofHusbandryandPharmaceuticalSciencesofChineseAcademyofAgriculturalSciences,Lanzhou730050,China;2.KeyLaboratoryofYakBreedingEngineeringofGansuProvince,Lanzhou730050,China)
This study screened the candidate genes of horn trait which provides scientific basis for the horn developmental mechanism and the polled yak breeding.The horn bud and forehead skin tissue from horn and polled yaks with 90 d gestational age were collected,and the mRNA different expression ofOLIG1,OLIG2,IFNAR1,IFNAR2,C1H21orf62,GCFC1,IFNGR2,SYNJ1,IL10RB,GART,RXFP2,FOXL2,TWIST1,TWIST2,ZEB2,E-cadherin,P-cadheringenes were analyzed using real-time quantitative PCR(qRT-PCR).The results showed that expressions ofOLIG1,IFNAR1,C1H21orf62,SYNJ1,IL10RB,GART,RXFP2,TWIST2 andE-cadheringenes showed no significant differences in horn and skin between horned and polled yaks (P>0.05).There were significant difference at transcript levels ofIFNGR2,FOXL2,TWIST1 andP-cadheringenes in horn buds between horned and polled yaks(P<0.05).IFNAR2 andGCFC1 had significant transcription difference in the expression of horn bud and skin between horned and polled yaks (P<0.05).OLIG2 andZEB2 had highly significant transcription difference in horn bud expression between horned and polled yaks (P<0.01).We concluded thatOLIG2 andFOXL2 might be the key candidate genes.ZEB2,TWIST1 andIFNGR2 might also involve in the formation of horn.P-cadherinmight participate in the development of horn.
horn;polled trait;yak;candidate genes;qRT-PCR
10.11843/j.issn.0366-6964.2016.06.009
2015-11-02
牦牛遺傳資源與育種(CAAS-ASTIP-2014-LIHPS-01);科技部科技支撐計劃(2012BAD13B05)
佘平昌(1989-),男,云南普洱人,碩士生,主要從事家畜分子生物學研究,E-mail:sheemail@163.com
閻萍,研究員,博士生導師,主要從事動物遺傳育種與繁殖研究,E-mail:pingyanlz@163.com
S823.85;S813.3
A
0366-6964(2016)06-1147-07