馮青輝 斯琴其木格 蘇和 朝格賽 薛寶玲 格日勒



摘 要:研究不同生長階段肉羊骨骼肌肌球蛋白重鏈(MyHC)基因相對表達量的變化規律,實驗選擇自然放牧條件下的1、3、6、9、12、18月齡的烏珠穆沁羊,屠宰后采取股二頭肌、背最長肌,采用十二烷基磺酸鈉-聚丙烯酰胺凝膠電泳(sodium dodecyl sulfate polyacrylamide gel electropheresis,SDS-PAGE)法分離骨骼肌中肌球蛋白重鏈異構體,利用實時熒光定量聚合酶鏈式反應(polymerase chain reaction,PCR)技術測定MyHC基因的相對表達量并探討其變化規律。結果表明:從SDS-PAGE圖中可以清晰地觀察到4條條帶分別為Ⅱa、Ⅱx、Ⅱb型與Ⅰ型。其中烏珠穆沁羊股二頭肌和背最長肌中MyHC的Ⅰ型和Ⅱa型基因的相對表達量均高于Ⅱb型基因的相對表達量。
關鍵詞:烏珠穆沁羊;骨骼肌;肌球蛋白重鏈基因;實時熒光定量聚合酶鏈式反應
Changes in Myosin Heavy Chain (MyHC) Gene from Skeletal Muscle during the Growth of Wuzhumuqin Sheep
FENG Qinghui1, Siqinqimuge1, SU He2, CHAO Gesai3, XUE Baoling1, BORJIGIN Gerelt1,*
(1.College of Food Science and Engineering, Inner Mongolia Agricultural University, Hohhot 010018, China;
2.Inner Mongolia Xilinhot Husbandry Bureau, Xilinhot 026000, China;
3.Inner Mongolia Xilinhgol Farming and Animal Husbandry Management, Xilinhot 026000, China)
Abstract: This study was aimed to investigate the changes in the relative expression level of myosin heavy chain (MyHC) gene from skeletal muscle during the growth of Wuzhumuqin sheep. Free-grazing Wuzhumuqin sheep were slaughtered at 1, 3, 6, 9, 12 and 18 months of age and the biceps femoris and longissimus dorsi muscles were excised. MyHC isoforms were separated from the muscles by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), and the relative expression level of MyHC gene was measured by real-time polymerase chain reaction (PCR). The results showed that four clear bands representing type Ⅱa, type Ⅱx, type Ⅱb and type Ⅰ, respectively, were observed in the electrophoretic pattern. The expression of MyHC type Ⅰ gene and type Ⅱa in biceps femoris and longissimus dorsi was higher than that of MyHC type Ⅱb gene. Understanding the MyHC gene expression in skeletal muscle during the growth of sheep is of great practical importance for improving meat quality and post-mortem physiological and biochemical characteristics.
Key words: Wuzhumuqin sheep; skeletal muscle; myosin heavy chain; real time-polymerase chain reaction (PCR)
中圖分類號:TS251.53 文獻標志碼:A 文章編號:1001-8123(2015)07-0020-05
doi: 10.7506/rlyj1001-8123-201507005
骨骼肌是由多核細胞的肌纖維(肌細胞)組成。肌纖維由很多肌原纖維組成,而肌原纖維是由很多不同蛋白質而構成。其中肌球蛋白是一個由高度保守的多基因編碼,含有多種同工蛋白的家族,它也是骨骼肌的主要收縮蛋白[1]。肌球蛋白分子由2個肌球蛋白重鏈和兩對起調節作用的肌球蛋白輕鏈組成,具有ATPase活性[2]。在多種肌球蛋白同工型基因共同表達模式的作用下,肌纖維的表型呈多樣性[3]。依據肌球蛋白重鏈(MyHC)基因的類型可將肌纖維分為4種,即MyHCⅠ、MyHCⅡa、MyHCⅡb和MyHCⅡx,它們分別特征性的對應4種不同類型的肌纖維Ⅰ、Ⅱa、Ⅱb和Ⅱx。在代謝上分別與慢速氧化型、快速氧化型、快速酵解型和中間類型肌纖維相對應[4-5]。肌纖維的類型及組成不僅決定了骨骼肌的生長發育,還會影響肉的生理生化特性及感官品質[6-7]。目前,國內外關于肌球蛋白重鏈基因表達規律的研究有很多。Changa等[8]指出肌球蛋白的分子類型依靠肌球蛋白重鏈相關基因來有效表達,通過MyHCⅠ、MyHCⅡa、MyHCⅡb、MyHCⅡx在骨骼肌上的表達和依據肌纖維與肌球蛋白重鏈的相關基因的特異性,動物骨骼肌肌纖維可對應的進行分型。楊曉靜[9]對豬骨骼肌生長及肌纖維類型分布的分子機理進行了研究。楊飛云[10]
對豬骨骼肌肌纖維類型分布及轉化的分子機理進行了研究。郭佳等[11]研究了金華豬背最長肌中不同肌球蛋白重鏈基因表達的發育性變化。張輝等[12]研究了FoxO1、肌纖維類型標志基因和MyHCⅠ、MyHCⅡb、MyHCⅡa和MyHCⅡx在特定骨骼肌中的表達規律,結果顯示,FoxO1在不同骨骼肌類型中mRNA表達差異不顯著(P>0.05),其蛋白表達與MyHC各亞型顯著相關。沈潔等[13]
研究了家兔骨骼肌肌球蛋白重鏈構成的電泳。范曉華等[14]
研究了大鼠脊髓橫斷早期比目魚肌質量與肌球蛋白重鏈亞型mRNA表達的變化。魯建清等[15]研究了骨骼肌肌球蛋白重鏈與不同類型訓練的影響。辛雪等[16]研究了巴美肉羊肌球蛋白重鏈基因的表達規律。此外有關魚、牛和豬等動物身上也有較多的肌球蛋白重鏈基因研究報道[17-18]。
而有關對肉羊身上的研究很少,特別是肉羊生長過程中的肌球蛋白重鏈基因的表達研究極少。烏珠穆沁羊是我國優良肉用品種,集中在烏珠穆沁區域,近幾年主要分布在錫林郭勒盟的東、西烏珠穆沁旗、阿巴哈納爾旗部分地區[19]。烏珠穆沁羊是經過自然放牧和人工選擇逐步繁育而形成的內蒙古地方優良品種之一,屬肉用粗毛羊;具有體大、肉質嫩、肉中脂肪含量適中、含有多種對人體有益的礦物質并且發育生長迅速等特點。本實驗對烏珠穆沁肉羊隨著月齡的增加MyHCⅠ、MyHCⅡa、MyHCⅡb和MyHCⅡx各基因的相對表達量進行分析,旨在掌握肉羊隨月齡的增加骨骼肌肌球蛋白重鏈基因的相對表達量的變化規律。
1 材料與方法
1.1 材料與試劑
內蒙古錫林郭勒盟自然放牧條件下的同一個羊群的烏珠穆沁羊各3 只,分別是1、3、6、9、12、18 月齡。將其屠宰后取股二頭肌和背最長肌,樣品分裝,然后液氮中速凍,帶回實驗室放入-80 ℃保存。
Tris試劑 南京博爾迪生物科技有限公司;丙烯酰胺、蔗糖 天津永大化學試劑廠;雙丙烯酰胺 上海紫一試劑廠;三氯甲烷 天津市科盟化工工貿有限公司;異丙醇 南京市新源化工;核酸染料 北京全式金生物技術有限公司;焦碳酸二乙酯、反轉錄試劑盒 大連寶生物有限公司。
1.2 儀器與設備
BG-power600型穩流恒壓電泳儀 北京白晶生物技術有限公司;LHS-50CH恒溫恒濕箱 上海一恒科學儀器有限公司;TU-1810紫外-可見分光光度計 北京普析通用儀器有限責任公司;ChampGel6000凝膠成像儀北京賽智創業科技有限公司;BG-Qspin600轉微型離心機 北京百晶生物技術有限公司;CFX96TM實時熒光定量
PCR儀 美國Bio-Rad公司;DSH-300A型回旋式恒溫水浴振蕩器 上海亞榮生化儀器廠。
1.3 方法
1.3.1 肌球蛋白的提取
骨骼肌肌球蛋白的提取方法參考Light等[20]的方法。取去掉脂肪以及大塊結締組織的肌肉組織切成碎塊,稱取10 mg。向10 mg的肌肉組織中加入500 μL蛋白質提取試劑Pro-Prep,3000V、0.3 min,勻漿3 次。將樣品冰上靜置30 min。4 ℃、13 000 r/min離心5min,分離上清液,并將其裝好放入-20 ℃環境下保存備用。
1.3.2 SDS-PAGE電泳分離肌球蛋白重鏈異構體
SDS-PAGE電泳分離肌球蛋白重鏈異構體參考Light等[20]的方法進行分離。進行SDS-PAGE電泳,待藍紫色線距離凝膠板底部2 cm處停止電泳。最后4 ℃染色約為2 h,4 ℃脫色,每1.5~2 h更換一次脫色液直至電泳條帶清晰。
1.3.3 總RNA的提取
總RNA的提取方法參考Trizol法[21]。取100~200 mg的烏珠穆沁羊的股二頭肌和背最長肌置于液氮中研磨,加液氮時應迅速、少量多次,并且肌肉組織應始終處于冷凍狀態下。將標號1.5 mL無酶無菌EP管稱量、去皮,然后拿小藥匙取100 mg肌肉粉末放入其中。在1.5 mL無酶無菌EP管中加1 mL的Trizol,冰浴5 min。4 ℃、12 000 r/min,離心5 min,吸取上清液(約1 000 μL)。加200 μL氯仿(1/5體積),蓋緊管蓋,劇烈振蕩15 s,無分相沉淀后室溫靜置5 min。4 ℃、12 000 r/min,離心15 min(離心后勿振蕩或混勻)這時樣品分為3 層,上層為RNA、中間層為DNA、下層為蛋白質。之后小心吸取上清液約500 μL轉移到新的無酶無菌EP管中。加入等體積異丙醇,上下混合均勻,室溫放置10 min。4 ℃、12 000 r/min,離心10 min后,得白色沉淀,丟棄上清。離心后可在管側壁和管底看見少量的RNA膠狀沉淀。管壁中加入焦碳酸二乙酯處理的75%乙醇1 mL,然后輕輕上下顛倒洗滌沉淀。4 ℃、12 000 r/min,離心5 min,棄上清。輕輕吹打沉淀,等RNA完全溶解后置于-80 ℃保存。分別取5 μLRNA溶液、2 μL緩沖液混勻,然后用1%瓊脂糖凝膠電泳來檢測總RNA質量的完整性,分析RNA的濃度和純度。
1.3.4 反轉錄
利用大連寶生物工程有限公司反轉錄試劑盒的指導說明進行反轉錄操作,進而合成cDNA。
1.3.5 引物設計及合成
在National Center For Biotechnology Information(NCFI)網中查找到綿羊肌球蛋白重鏈MyHCⅠ、MyHCⅡa、MyHCⅡb和MyHCⅡx各基因的核苷酸序列,然后應用Primer Premier 5.0軟件來對引物進行分析,之后由博覽特生物技術有限公司進行引物的設計和合成。
1.3.6 PCR條件優化
由于目的基因PCR引物長度的最適退火溫度與PCR反應體系的最適退火溫度不同,cDNA濃度也不同,導致在做實時熒光定量PCR實驗前,對PCR反應體系來進行溫度梯度實驗,以此來探究最佳退火溫度。
1.3.7 實時熒光定量PCR
利用CFX96TMReal Time System實時PCR系統來進行實時定量擴增。用反轉錄的cDNA為模板,看家基因β-actin為內參,每個樣品進行3個重復。Real-Time PCR體系中熒光物質的特異性與目的基因的目的片段來進行結合,CFX96TMReal Time System實時的PCR系統收集到熒光信號量與PCR產物量成正比。實時熒光定量PCR反應總體系為25μL,其中:滅菌ddH2O 8.5μL,cDNA的模板2μL,上游引物(10μmol/L)1μL,下游引物(10μmol/L)1μL,SYBR Premix Ex TaqⅡ(2×)12.5μL。
采用兩步法來進行PCR反應程序:95℃預變性30s,60℃退火30s,總共40個循環。基因相對表達采用2-ΔΔCt公式計算。
2 結果與分析
2.1 RNA的質量評定
用1%瓊脂糖凝膠電泳來檢測總RNA質量的完整性。圖1是股二頭肌總RNA瓊脂糖凝膠電泳,圖2是背最長肌總RNA瓊脂糖凝膠電泳。電泳檢測的結果顯示28S、18S、5S條帶較為清晰,說明RNA幾乎沒有被降解,可在cDNA的合成和RT-PCR中使用。
圖 1 股二頭肌總RNA質量鑒定結果
Fig.1 Identification of total RNA isolated from biceps femoris
圖 2 背最長肌總RNA質量鑒定結果
Fig.2 Identification of total RNA isolated from longissimus dorsi
2.2 PCR產物擴增結果
圖 3 普通RT-PCR產物瓊脂糖電泳圖
Fig.3 Agarose gel electrophoresis of RT-PCR products
由圖3可知,普通RT-PCR擴增出的MyHC基因的Ⅰ、Ⅱa、Ⅱb、Ⅱx型和β-actin基因長度分別為165、200、216、284、233bp。能夠清晰地觀察到普通RT-PCR擴增出的目的片段預期的大小,并且沒有其他條帶,這說明目的基因片段產生明顯的特異性擴增條帶,所設計引物的特異性較好,符合后期RT-PCR實驗要求。
2.3 MyHCⅠ型基因在股二頭肌和背最長肌中的相對表達量
大寫字母不同,表示不同月齡、相同部位之間差異顯著(P<0.05)。下同。
圖 4 不同月齡烏珠穆沁羊股二頭肌和背最長肌中MyHCⅠ基因
mRNA相對表達量
Fig.4 Relative expression of MyHCⅠ mRNA in biceps femoris and longissimus dorsi of Wuzhumuqin sheep at different ages
由圖4可知,在股二頭肌中MyHCⅠ型基因的mRNA相對表達量有不同變化現象。1月齡>3月齡>6月齡>9月齡≈12月齡>18月齡,總體上1~18月齡有明顯的下降趨勢,但9、12月齡之間的差異不顯著。1、3、6、9、18月齡之間的差異顯著(P<0.05)。在不同的生長階段中背最長肌的MyHCⅠ型基因mRNA在6月齡時的相對表達量顯著較高于其他月齡段(P<0.05)。1、3、9月齡之間差異均不顯著。
2.4 MyHCⅡb型基因在股二頭肌和背最長肌中的相對表達量
由圖5可知,在不同生長階段MyHCⅡb型基因mRNA相對表達量在不同骨骼肌中的表達有所不同。在股二頭肌中,12月齡的相對表達量與其他月齡相比有顯著性差異(P<0.05),在3、6、18月齡之間的差異不顯著。不同生長階段背最長肌中MyHCⅡb型基因mRNA相對表達量變化在3月齡時有較高的表達量,與其他月齡相比有顯著性差異(P<0.05)。1月齡和6月齡之間差異不顯著,9、12、18月齡之間差異均不顯著。
圖 5 不同月齡烏珠穆沁羊股二頭肌和背最長肌中MyHCⅡb基因mRNA相對表達量
Fig.5 Relative expression of MyHCⅡb mRNA in biceps femoris and longissimus dorsi of Wuzhumuqin sheep at different ages
2.5 MyHCⅡa型基因在股二頭肌和背最長肌中的相對
表達量
圖 6 不同月齡烏珠穆沁羊股二頭肌和背最長肌中MyHCⅡa基因mRNA相對表達量
Fig.6 Relative expression of MyHCⅡa mRNA in biceps femoris and longissimus dorsi of Wuzhumuqin sheep at different ages
由圖6可知,在不同生長階段MyHCⅡa型基因的mRNA相對表達量在股二頭肌中的變化為18月齡>9月齡≈6月齡>3月齡≈12月齡>1月齡。1~18月齡有逐步升高的趨勢。18月齡相對表達量與其他月齡的表達量之間差異顯著(P<0.05)。3、12月齡差異不顯著。6、9月齡之間差異不顯著。不同生長階段MyHCⅡa型的基因mRNA相對表達量在背最長肌中的變化為9月齡>1月齡>12月齡≈18月齡>3月齡>
6月齡,9月齡與其他月齡之間差異顯著(P<0.05),12、18月齡之間差異均不顯著。
2.6 MyHCⅡx型基因在股二頭肌和背最長肌中的相對表達量
由圖7可知,在不同生長階段MyHCⅡx型基因mRNA的相對表達量在股二頭肌中的變化為12月齡>3月齡≈
18月齡>6月齡>1月齡>9月齡。12月齡與其他月齡之間差異顯著(P<0.05),3、18月齡之間差異不顯著
(P<0.05)。不同生長階段背最長肌中MyHCⅡx型基因mRNA相對表達量變化為9月齡>1月齡>12月齡≈18月齡>3月齡≈6月齡,9月齡的表達與其他月齡相比有顯著性差異(P<0.05),3、6、12、18月齡之間差異均不顯著。
圖 7 不同月齡烏珠穆沁羊股二頭肌和背最長肌中MyHCⅡx基因mRNA相對表達量
Fig.7 Relative expression of MyHCⅡx mRNA in biceps femoris and longissimus dorsi of Wuzhumuqin sheep at different ages
3 討 論
總體上看,不同類型的肌球蛋白重鏈基因的特異性表達與不同類型的肌纖維(不同部位骨骼肌)之間存在著對應關系。本結果與Devoto等[22]的實驗結論一致。家畜在它的生長以及運動等過程中骨骼肌中不同肌球蛋白重鏈相關基因之間也存在著相互轉化的關系,這種轉化關系也是導致不同類型肌纖維之間相互轉化的誘導因素[23]。動物體肌纖維在色澤、直徑、供能方式、酶活性、收縮速度上存在差異。不同類型肌纖維的生長具有特異性,而且其受遺傳、營養和生活環境等因素的影響在生長過程中不斷發生轉化。事實上,各種肌球蛋白重鏈和肌球蛋白輕鏈異形體在肌纖維中的分布并不是嚴格的局限于某種肌纖維或肌肉發育的特定階段。肌球蛋白重鏈和肌球蛋白輕鏈異形體有多種多樣的結合形式,構成了大量的肌球蛋白分子。肌肉中肌纖維類型的構成、機能、形態學上的差異與酮體肉質性狀、產肉量以及脂肪交雜、風味等性能方面有著密切關系。例如,日本和牛是以優良肉質而聞名世界。研究表明,日本國內產的肉用牛因品種的不同,肌纖維類型在特征上有很大的區別,即黑色和牛中荷爾斯泰因品種的βR型(紅色肌纖維)比褐色和牛的多,這可能與肉質的改良有關系。另外,同一種黑色和牛中因飼養條件、產地等的不同,肌纖維類型也不同。肌纖維類型與產肉性質,特別是體質量增加和脂肪交雜的關系更值得重視;到目前為止還未明確的解析它們之間的關系。有研究表明動物出生后隨著生長發育肌纖維增大肌肉質量也增加,白色肌纖維的直徑比紅色肌纖維變化大,從αR型肌纖維轉移αW型。有關動物出生后肌纖維的分化方面Ashmore等[24]指出,牛出生時所有肌纖維是紅色肌纖維,純紅色肌纖維的βR型不引起出生后分化,而αR型的一部分分化為αW型。這些因素可能關系到動物產肉量的增加。研究者們認為通過育種、改良來增加肌肉質量主要關系到白色肌纖維的構成比率。肉質,特別是脂肪交雜的肌肉,其中紅色肌纖維的比例越高優質的脂肪交雜越高,例如黑色和牛。因此,不同種類骨骼肌骨骼肌之間,肌內蓄積脂肪含量比率有顯著差別。因此研究肌纖維的生成、結構、類型及其轉化規律有著重要的實際意義。肌肉中肌纖維類型的組成對宰后肉品質、生理生化等特性具有重要影響。
4 結 論
總體上烏珠穆沁羊股二頭肌中MyHCⅠ型基因的相對表達量高于MyHCⅡb型基因的相對表達量,并MyHCⅠ型和MyHCⅡa型基因的相對表達量高于MyHCⅡb型基因的相對表達量。烏珠穆沁羊不同生長階段背最長肌中除了3月齡外,其他月齡MyHCⅠ型和MyHCⅡa型基因的相對表達量都要高于MyHCⅡb型基因的相對表達量。
參考文獻:
[1] CHIKUNI K, TANABE R, MUROYA S, et al. Difference inmolecular structure among the procine myosin heavy chain-2a, -2x and -2b isoforms[J]. Meat Science, 2001, 57(3): 311-317.
[2] WEISS A, LEINWAND L A. The mammalian myosin heavychain gene family[J]. Annual Reviews of Cell & Development Biology, 1996, 12(1): 417-439.
[3] EGGERT J M, DEPREUX F F S, SCHINEKEL A P, et al. Myosin heavychain isoforms account for variation in pork quality[J]. Meat Science, 2002, 61(2): 117-126.
[4] SCHIAFFINO S, REGGIANI C. Molecular diversity of myofibrillar proteins: gene regulation and functional significance[J]. Physiological Reviews, 1996, 76(2): 371-423.
[5] CHANG K C, FEMANDES K. Developmental expression and 5cloning of the porcine 2x and 2b myosin heavy chaingenes[J]. DNA and Cell Biology, 1997, 16(12): 1429-1437.
[6] 尹靖東. 動物肌肉生物學與肉品科學[M]. 北京: 中國農業大學出版社, 2011: 164-183.
[7] LEFAUCHEUR L. A second look into fibre typing-relation tomeat quality[J]. Meat Science, 2010, 84(2): 257-270.
[8] CHANGA K C, COSTA N D, BLACKLEY R, et al. Relationships of myosin heavy chain fibre types to meat quality traits intraditional and modem pigs[J]. Meat Science, 2003, 64(1): 93-103.
[9] 楊曉靜. 豬骨骼肌生長及肌纖維類型分布的分子機理研究[D]. 南京: 南京農業大學, 2004: 28-35.
[10] 楊飛云. 豬骨骼肌肌纖維類型分布及轉化的分子機理研究[D]. 雅安: 四川農業大學, 2008: 33-38.
[11] 郭佳, 徐娥, 黃明, 等. 金華豬背最長肌中不同肌球蛋白重鏈基因表達的發育性變化[J]. 中國畜牧雜志, 2012, 48(5): 14-17.
[12] 張輝, 史新娥, 袁媛, 等. FoxO1抑制豬骨骼肌MyHCⅠ的表達[J]. 中國生物化學與分子生物學報, 2010, 26(3): 283-289.
[13] 沈潔, 徐大剛, 朱道立, 等. 家兔骨骼肌肌球蛋白重鏈構成的電泳[J]. 安徽農業科學, 2008, 36(1): 46-49.
[14] 范曉華, 紀樹榮, 周紅俊, 等. 大鼠脊髓橫斷早期比目魚肌重量與肌球蛋白重鏈亞型mRNA表達的變化[J]. 中國康復理論與實踐, 2006, 12(4): 298-300.
[15] 魯健清, 湯長發, 尹麗琴, 等. 骨骼肌肌球蛋白重鏈與不同類型訓練的影響[J]. 中國組織工程研究與臨床康復, 2010, 14(28): 5289-5292.
[16] 辛雪, 蘇琳, 趙雅娟, 等. 巴美肉羊肌球蛋白重鏈基因表達規律的研究[J]. 肉類研究, 2014, 26(10): 134-138.
[17] FRANCISCO C L, JORGE A M. Muscle fiber type characterization and myosin heavy chain (MyHC) isoform expression in Mediterranean buffaloes[J]. Meat Science, 2011, 88(3): 535-541.
[18] 佚名. 烏珠穆沁羊繁育基地建設情況[J]. 新疆畜牧業, 2009(增刊2): 43-44.
[19] 劉震乙, 文浴蘭, 沙里, 等. 烏珠穆沁羊品種志[J]. 內蒙古農牧學院學報, 1982(1): 1-6.
[20] LIGHT N, CHAMPION A E, VOYLE C, et al. The role of epimysial, perimysial and endomysial collagen in determining texture in six bovine muscles[J]. Meat Science, 1985, 13(3): 137-149.
[21] 于寒松, 彭帥, 謝遠紅, 等. 一種RNA提取試劑盒: Trizol的使用方法初探[J]. 食品科學, 2005, 26(11): 39-42.
[22] DEVOTO S H, MELANCON E, EISEN J S, et al. Identification of separate slow and fast muscle precursor cells in vivo prior to somite formation[J]. Development, 1996, 122(11): 3371-3380.
[23] 高強. 肌球蛋白重鏈基因MyHC與中國地方雞肌肉品質之間的相關研究[D]. 北京: 中國農業大學, 2007: 20-28.
[24] ASHMORE C R, DOERR L. Comparative aspects of muscle fiber types in different species[J]. Experimental Neurology, 1971, 31(3): 408-418.