趙斌 李挺 趙博昊 穆琳 胡帥帥 陳陽 吳信生

摘要[目的]對家兔Myh7基因進行生物信息學分析及表達規(guī)律研究。[方法]利用RT-PCR技術檢測了Myh7在肌肉發(fā)育過程中的表達規(guī)律，并結合生物信息學方法對家兔Myh7的序列進行深入分析。[結果]Myh7蛋白為不穩(wěn)定親水性蛋白，有Myosin_N和MYSc_class_II這2個結構區(qū)域，其中MYSc_class_II是肌球蛋白重鏈家族特有結構域。MEGA軟件分析家兔Myh7分子進化，發(fā)現(xiàn)其與人、猩猩等物種間同源性較高，與其親緣關系的遠近一致。定量分析顯示：在整個肌肉發(fā)育過程中，Mhy7表達量在腿肌中基本維持在較低水平，在背肌中的表達量相對波動較大，而在11周齡時出現(xiàn)表達量上升的趨勢。[結論]該研究為進一步研究Myh7基因在肌肉生長過程中的功能奠定了分子基礎，為研究家兔的肌肉發(fā)育分子機制提供了一定的參考依據(jù)。

關鍵詞兔；肌肉；Mhy7基因；生物信息學分析；表達規(guī)律

中圖分類號S829.1文獻標識碼

A文章編號0517-6611（2017）16-0130-03

Bioinformatics Analysis and Expression Characteristics of Myh7 Gene in Rabbit

ZHAO Bin，LI Ting，ZHAO Bohao，WU Xinsheng* et al（

Jiangsu Key Laboratory of Animal Genetics & Breeding and Mdecular Designing，College of Animal Science and Technology，Yangzhou University，Yangzhou，Jiangsu 225009）

Abstract[Objective]Bioinformatics analysis and expnession characteristics of Myh7 gene in rabbit were studied.[Method]The expression characteristics of Myh7 was detected by RTPCR technique in the process of muscle development，and we further used bioinformatics methods to analyze rabbit Myh7 sequence.[Result]Myh7 protein was unstable and hydrophilic，with the Myosin_N and MYSc_class_II structure domain，MYSc_class_II was the special structure domain of myosin heavy chain.We used MEGA software to analyze Myh7 molecular evolution of rabbit，the result showed that homology was high among rabbit，human and orangutan，which was conforming to the kinship.Quantitative analysis showed that Mhy7 in leg muscle expressed in the low level，and the expression was fluctuant in the back muscle，increased in 11th week old rabbit.[Conclusion] The study provided the molecular basis for exploring the Myh7 gene function in the process of muscle growth，and it could also provide certain reference basis for studying the molecular mechanism of rabbit muscle development.

Key wordsRabbit；Muscle；Mhy7 gene；Bioinformatics analysis；Expression characteristics

兔肉是一種具有高營養(yǎng)價值的肉類，其含有高氨基酸、高卵磷脂、高蛋白和低膽固醇等，在心腦血管病、肝病、糖尿病患者中較受歡迎。肌肉發(fā)育過程影響著兔肉的產(chǎn)量和質量，肌肉發(fā)育包含胚胎期肌纖維的形成、肌纖維的發(fā)育和成年期肌肉的再生等，此發(fā)育過程受到許多基因的相關調控。目前，在人和小鼠上已發(fā)現(xiàn)在成肌定向分化早期，Myf5基因和MyoD基因起著關鍵作用[1]；肌肉干細胞最終分化為肌管肌纖維的能力與Mrf4基因和MyoG基因調控密切相關[2]，這些基因都屬于成肌調節(jié)因子基因家族。除此之外，肌細胞增強因子基因家族、Wnt基因家族也是能直接調控肌肉發(fā)育的一類家族基因[3-4]。筆者所在課題組前期通過Agilent兔全基因組表達譜芯片技術篩選出肌肉中具有差異表達的基因，在這些差異表達的基因中，發(fā)現(xiàn)Myh7基因與新西蘭白兔肌肉的生長發(fā)育密切相關。該研究利用生物信息學分析方法，對Myh7基因的理化性質和進化進行了分析和預測，同時利用RT-PCR方法檢測不同部位不同時期肌肉中Myh7基因的表達規(guī)律，旨在探究Myh7基因在肌肉發(fā)育過程中作用，為研究家兔的肌肉發(fā)育分子機制提供一定的參考依據(jù)。

1材料與方法

1.1材料與試劑采集1、3、5、7、9、11、13、15、17、19、21、23周齡新西蘭白兔背部和腿部肌肉，公母各1/2，分別取0.5 g 置于無酶管中，液氮凍存帶回實驗室-70 ℃保存?zhèn)溆谩?/p>

熒光定量試劑盒（CW0956A）、SuperRT cDNA第一鏈合成試劑盒（CW0741A）等購自北京康為世紀生物科技有限公司；上海生工生物工程有限公司提供DNA marker系列、無水乙醇等常規(guī)試劑；RNAstore樣本保存液購自TIANGEN公司。

1.2方法

1.2.1生物信息學分析。

蛋白的氨基酸序列組成等理化性質采用ProtParam軟件進行分析，基因編碼蛋白質的性質通過ExPASy工具對進行分析，運用NCBI的CDD軟件（http：//www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi）預測Myh7蛋白結構域；氨基酸序列的親水性分析通過ProtScale（http：//web.expasy.org/protscale/）進行，運用MEGA軟件分析Myh7蛋白分子進化；運用TMPRED（http：//www.ch.embnet.org/software/）分析預測蛋白質的二級結構，運用PSORT II Prediction（http：//www.psort.org/）對Myh7的亞細胞定位進行預測分析。

1.2.2總RNA提取及cDNA 合成。

用傳統(tǒng) Trizol 法提取樣品的總RNA，分光光度計測定總 RNA 的純度和濃度，立即對提取后的總RNA進行反轉錄。取0.5 mL Ependorf 管，以總RNA 提取液2 μg為模板， oligo d（T）為反轉錄引物，在反轉錄酶 M-M LVRT 作用下進行反轉錄，反應體系為25 μL。

1.2.3RT-PCR法對Mhy7基因表達差異進行檢測。根據(jù)已報道的Mhy7基因cDNA 部分序列設計引物，該引物由生工生物工程（上海）有限公司合成，引物信息見表1。北京康為世紀生物科技有限公司提供熒光定量PCR預混體系，表達水平以GAPDH基因作為內參進行均一化處理。

2結果與分析

2.1家兔Mhy7的生物信息學分析

2.1.1家兔Mhy7蛋白理化性質分析。

根據(jù) GenBank 上已發(fā)表的家兔的Myh7基因序列可知，家兔的Myh7基因全長20 869 bp。其中，編碼序列由5 808個堿基組成，編碼1 935個氨基酸殘基，是一個較為保守的基因。經(jīng)ProtParam在線軟件預測，分子式為C9 722H15 810N2 772O3 090S66，原子總數(shù)為31 460個，分子量為223 087.1 ku，理論等電點為5.61。在所有編碼氨基酸中，Glu （E）含量最高，為13.2%，Trp （W） 只占0.5%（表2）。帶負電荷殘基的總數(shù)量（ASP + Glu）：359，帶正電荷殘基的總數(shù)量（Arg+Lys）：315，不穩(wěn)定指數(shù)是47.38，說明這個蛋白質不穩(wěn)定。

由圖1可知，Myh7蛋白疏水性最小值為-3.111，最大值為2.189。其主要親水平均數(shù)（GRAVY）值是-0.421，因此Myh7蛋白為親水蛋白。Myh7主要分布在細胞核（69.6 %）中，少量定位在細胞質中（17.4 %），說明Myh7主要在細胞核中發(fā)揮生物學作用。

2.1.2家兔Myh7蛋白的分子進化分析。

將人、豬、牛、羊、馬、猩猩、貓、狗、小鼠、雞的Myh7基因編碼蛋白及兔Myh7蛋白序列調入MEGA中，進行系統(tǒng)進化分析。由圖2可知，進化樹總體分為兩支，家禽雞本身獨立成為一支，其余的10個物種成立為一個獨立大分支，而家兔Myh7基因與人、猩猩Myh7基因的分子進化距離最近，親緣關系最近。

2.1.3家兔Mhy7保守結構域預測分析。

1 935個氨基酸由家兔Mhy7基因的CDS區(qū)編碼，利用NCBI在線工具CDD（http：//www.ncbi.nlm.nih.gov/Structure/cdd/）進行保守結構域分析，結果顯示：Myh7蛋白有2個特殊的區(qū)域，N端34～75 aa區(qū)段為Myosin_N結構域，79～799 aa區(qū)段為MYSc_class_II結構域。其中N端的Myosin_N區(qū)是肌凝蛋白SH3樣結構域，而MYSc_class_II區(qū)由2個重鏈組成，頭和尾域的卷曲螺旋把2個重鏈結合在一起，它有ATP的結合位點和P環(huán)結構（圖3）。

2.1.4家兔Mhy7蛋白結構預測分析。

對蛋白質二級結構進行預測分析有助于認識蛋白質的空間結構，因為二級結構主要是多肽鏈主鏈骨架中局部的構象[5]。使用TMPRED和Swiss-Model軟件分別對家兔Myh7的二級和三級結構進行預測。預測發(fā)現(xiàn)Myh7蛋白二級結構以α-螺旋為主，其中α螺旋（h）74.94%，無規(guī)卷曲（c）19.17%，延伸鏈（e）5.89%；而Myh7蛋白三級結構由2條互相交織程度不明顯、較為簡單且長度不等的螺旋鏈組成，呈現(xiàn)出鏈狀分布，該蛋白與二級結構預測結果基本一致（圖4）。

2.2不同生長階段Mhy7基因表達規(guī)律

分析新西蘭白兔在生長發(fā)育不同階段背、腿肌Myh7基因的表達規(guī)律，從圖5可以看出，背肌中Myh7的表達量在1～23周齡這一生長發(fā)育時期，所呈現(xiàn)出的波動較大，而在腿肌中Myh7基因的表達量基本處在較低水平；1～9周齡Myh7基因在背肌中的表達量呈現(xiàn)出下降趨勢，11周齡Myh7基因的表達量在腿肌中開始出現(xiàn)上升趨勢；而腿肌中的表達量在13周齡以后出現(xiàn)下降趨勢，最后維持在較低水平，Myh7基因在背肌中表達量相對于腿肌而言總體偏高，且呈現(xiàn)出不規(guī)則變化。

3討論

肌肉的生長發(fā)育包含胚胎期肌纖維的形成、出生后肌纖維的發(fā)育和成年期肌肉的再生等步驟，是一個復雜的過程。從整體上看，整個過程伴隨大量的調控基因在轉錄組水平和通路水平等多層次的調控，為成肌始祖細胞發(fā)育成為適應各種生命過程的終末肌纖維提供保障。這些調控包含成肌調節(jié)因子、miRNA、表觀遺傳、肌肉生長抑制素、GH-IGFs通路、Wnt/β-catenin 通路等[6]。

Myh7基因是肌球蛋白重鏈家族基因中的一員，肌球蛋白重鏈家族基因在肌肉的生長發(fā)育中扮演重要角色，通過前期差異表達基因篩選，發(fā)現(xiàn)Myh7基因與新西蘭白兔肌肉發(fā)育密切相關。已有研究表明，Myh7基因的突變能夠引發(fā)人類肥厚型心肌病，30%～50%的家族性肥厚型心肌病由Myh7基因的突變所致[7-9]。王虎等[10]對一個中國人肥厚型心肌病家族進行研究，發(fā)現(xiàn)該病與Myh7基因的第14號外顯子雜合突變c.1273G>A相關，且攜帶者的臨床表型存在明顯差異。Myh7基因nt2025C插入和nt2013c缺失被楊春等[11]在一家族中發(fā)現(xiàn)，屬國內罕有的移碼突變，與該家族肥厚型心肌病患者中存在的致病基因的突變相關。楊忠偉等[12]對3個HCM家族成員進行研究，發(fā)現(xiàn)在1個家族中，

Myh7基因14外顯子存在Thr441Met突變，該突變首次在我國被發(fā)現(xiàn)。故可以通過Myh7基因突變的篩查使提前知曉家族性肥厚型心肌病病變的存在成為可能。

該試驗通過生物信息學方法對家兔Myh7基因進行分析研究， 結果表明Myh7蛋白為不穩(wěn)定親水性蛋白，有MYSc_class_II和Myosin_N這2個結構區(qū)域，其中MYSc_class_II是肌球蛋白重鏈家族特有結構域。通過MEGA軟件對家兔Myh7蛋白進行分子進化分析，發(fā)現(xiàn)家禽雞單獨成一個分支，而家兔與人、猩猩等物種間同源性較高，與其親緣關系遠近相一致。利用RT-PCR方法研究了1～23周齡不同生長階段背、腿肌Mhy7基因的表達規(guī)律，結果顯示，Mhy7基因在家兔整個生長發(fā)育時期因不同部位表現(xiàn)出不完全一致的表達規(guī)律，其在腿肌中的表達量基本處在較低的水平；而在背肌中，Mhy7基因的表達量波動較大，提示Myh7基因在背肌中的調控作用更大。在11周齡時出現(xiàn)明顯上升，此時背肌和腿肌中Myh7基因的表達量均有上升。該研究初步探索了Mhy7基因在肌肉發(fā)育過程中的表達規(guī)律，并結合生物信息學方法進行探究，能為研究家兔的肌肉發(fā)育分子機制提供一定的參考依據(jù)。

安徽農業(yè)科學2017年

參考文獻

[1]

ARNOLD H H，WINTER B.Muscle differentiation：More complexity to the network myogenic regulators[J].EMBOJ，1989，8（12）：987-1000.

[2] ARNOLD H H，WINTER B.Muscle differentiation：More complexity to the network of myogenic regulators[J].Current opinion in genetics & development，1998，8（5）：539-544.

[3] GEETHALOGANATHAN P，NIMMAGADDA S，HUANG R，et al.Role of wnt6 in limb myogenesis[J].Anatomy and embryology，2006，211（3）：183-188.

[4] ANAKWE K，ROBSON L，HADLEY J，et al.Wnt signaling regulates myogenic differentiation in the developing avian wing[J].Development，2003，130（15）：3503-3514.

[5] KLONT R E，BROCKS L，EIKELENBOOM G.Musele fibre type and meat quality[J].Meat science，1998，49（S1）：219-229.

[6] 余章斌，韓樹萍，朱春，等.基因表達譜芯片篩選室間隔缺損胎兒心肌組織差異表達基因[J].實用兒科臨床雜志，2011，26（18）：1425-1428.

[7] PERROT A，SCHMIDTTRAUB H，HORRMANN B，et al.Prevalence of cardiac betamyosin heavy chain gene mutations in patients with hypertrophic cardiomyopathy[J].Journal of molecular medicine，2005，83（6）：468-477.

[8] ALCALA R，SEIDMAN J G，SEIDMAN C E.Genetic basis of hypertrophic cardiomyopathy：From bench to the clinics[J].Journal of cardiovascular electrophysiology，2008，19（1）：104-110.

[9] 楊忠偉，馮秀麗，王忠.利用變性高效液相色譜方法進行HCM家系MYH7基因突變篩查研究[J].農墾醫(yī)學，2012，34（2）：127-129.

[10] 王虎，鄒玉寶，宋雷，等.心臟肌球蛋白重鏈基因c.1273G>A 突變與肥厚型心肌病的關聯(lián)分析[J].遺傳，2009，31（5）：485-488.

[11] 楊春，王愛玲，余元勛，等.漢族家族性肥厚型心肌病MYH7基因突變及相應臨床特征[J].安徽醫(yī)科大學學報，2008，43（3）：336-339.

[12] 楊忠偉，馮秀麗，王偉，等.新疆3個HCM 家系MYH7 基因突變DHPLC 檢測和DNA 測序結果分析[J].基礎醫(yī)學與臨床，2012，32（12）：1437-1441.