邵 芳，費曉燕，盧祥云，徐建榮，顧志良

（常熟理工學院 生物與食品工程學院，江蘇 常熟 215500）

松江鱸（Trachidermus fasciatus）屬鲉形目，杜父魚科，松江鱸屬.近年來由于資源枯竭，種群瀕危，已被列為國家二級保護動物.目前對松江鱸的研究主要集中在生態(tài)、生理、胚胎發(fā)育與繁殖習性等方面[1].為有效保護和合理利用松江鱸資源，對其種群遺傳結(jié)構(gòu)和遺傳多樣性的研究就顯得尤為重要．

肌肉生長抑制素(Myostatin，MSTN)作為骨骼肌生長發(fā)育的負調(diào)節(jié)因子[2]，引起生物學家的廣泛關(guān)注.弄清MSTN基因的結(jié)構(gòu)和功能對闡明骨骼肌生長發(fā)育的調(diào)控機理具有重要意義.研究表明魚MSTN基因由三個外顯子和兩個內(nèi)含子組成[3-5]，我們對松江鱸MSTN基因的全長cDNA和組織表達特異性作了研究[6]，但目前對松江鱸MSTN基因的內(nèi)含子序列尚未有報道.mtDNA是獨立于細胞核基因組外的環(huán)狀雙鏈DNA分子，在mtDNA上存在一個獨特的D-loop結(jié)構(gòu)，是整個線粒體基因組序列和長度變異最大的區(qū)域，其進化速度最快.線粒體DNA的D-loop區(qū)可用于研究物種起源和進化以及種內(nèi)種群的系統(tǒng)進化分析，探討物種在進化過程中可能發(fā)生的變異情況.魚類mtDNA是一個相對獨立的復(fù)制單位，具有分子小、結(jié)構(gòu)簡單、重組率低、進化速度快、不同區(qū)域進化速度有差異等特點，是魚類進化遺傳學、分子生態(tài)學、遺傳多樣性及其保護生物學等研究的重要標記[7]．

本研究對MSTN基因內(nèi)含子1、2以及長江、丹東和葫蘆島的松江鱸線粒體D-loop區(qū)擴增和測序，檢測不同地理群體松江鱸線粒體D-loop區(qū)的多態(tài)性，探索種群內(nèi)和種群間的遺傳多樣性狀況，為今后松江鱸的資源保護、開發(fā)和養(yǎng)殖提供相應(yīng)的理論依據(jù)．

1 材料與方法

1.1 實驗材料：

本實驗所用的松江鱸來自長江（C）、丹東（D）、葫蘆島（H）三個地理群體.實驗所用材料由蘇州市長江特色水產(chǎn)繁殖工程中心提供.將魚宰殺后取肌肉組織在-80℃冰箱中保存?zhèn)溆?，用?氯仿法提取基因組DNA．

1.2 引物設(shè)計與合成

根據(jù)我們克隆的松江鱸Myostatin基因cDNA序列（Genbank登錄號：GU198192)，以及斑馬魚和鯉魚Myostatin基因DNA序列（Genbank登錄號：AY693972，GQ214770），分別設(shè)計引物Sjl intr1 F1/R1，Sjl intr2 F1/R1以擴增松江鱸Myostatin基因的內(nèi)含子1和內(nèi)含子2序列.根據(jù)GenBank數(shù)據(jù)庫中已公布松江鱸線粒體D-loop部分序列（Genbank登錄號：AB188172），設(shè)計引物SjldpF1/R1，擴增松江鱸D-Loop部分片段，以分析不同地理群體的關(guān)系.引物由上海生工生物工程有限公司合成．

表1 實驗所用的引物序列

1.3 松江鱸Myostatin基因內(nèi)含子的PCR擴增和克隆

PCR擴增在25μL反應(yīng)體系中進行：0.2μL ExTaq酶（5U·μL-1，TaKaRa Tokyo，Japan），2.5μL10×PCR ExTaq Buffer（with Mg2+），2μL dNTP（2.5mmol·L-1each），Sjl intr1 F1和Sjl intr1 R1各1μL（10μM），基因組DNA（100ng/μL）1μL，加滅菌水至25μL.PCR擴增條件為95℃變性3min；然后在95℃ 30 Sec，60℃30Sec，72℃1.5min，進行30次循環(huán)；72℃延伸10min；4℃保存.將PCR擴增產(chǎn)物于1%瓊脂糖凝膠電泳，切膠回收目的片段.將回收片段克隆到PCR2.1 Vector（Invitrogen，美國）中，挑取陽性菌落，于LB液體培養(yǎng)基（Amp+）中搖菌培養(yǎng)，提取質(zhì)粒，經(jīng)鑒定后將重組質(zhì)粒送往上海桑尼生物科技公司進行測序．

1.4 不同地理群體線粒體D-loop序列擴增和序列測定

PCR擴增反應(yīng)體系和反應(yīng)條件同上述1.3.將PCR產(chǎn)物進行瓊脂糖凝膠電泳檢測，回收純化目的片段，直接測序．

1.5 數(shù)據(jù)分析

將獲得的序列在DNAMAN5.0軟件上進行分析，并進行物種間同源性比較.利用MEGA4.1軟件中的Kimura 2-parameter法計算凈遺傳距離矩陣，以距離矩陣鄰接法（Neighbor-Joining，NJ）建系統(tǒng)發(fā)生樹，通過自引導(dǎo)檢驗（bootstrap）獲得系統(tǒng)分支的置信度（重復(fù)次數(shù)為1000）．

2 結(jié)果與分析

2.1 松江鱸MSTN基因內(nèi)含子的擴增

以松江鱸基因組DNA為模板，用引物Sjl intr1F1/R1和Sjl intr2F1/R1進行PCR擴增，PCR產(chǎn)物在1%瓊脂糖凝膠中電泳檢測（圖1）.電泳檢測結(jié)果顯示，分別得到了450bp和1000bpPCR擴增產(chǎn)物，與設(shè)計的目的片段長度相同.將這兩個片段分別回收后克隆到PCR2.1載體中，經(jīng)測序后獲得了兩個片段的核苷酸序列，與松江鱸的MSTN基因cDNA序列比對，確認該序列為松江鱸MSTN基因的內(nèi)含子1和內(nèi)含子2序列．

圖1 松江鱸MSTN基因內(nèi)含子PCR擴增產(chǎn)物的1%瓊脂糖凝膠電泳圖

2.2 松江鱸MSTN基因內(nèi)含子的序列分析

通過PCR擴增并克隆測序后獲得了松江鱸MSTN基因內(nèi)含子1和內(nèi)含子2的序列，內(nèi)含子1長度為325bp，內(nèi)含子2長度為835bp（圖2和3），兩個內(nèi)含子的兩端都符合GT-AG規(guī)則.在內(nèi)含子2中發(fā)現(xiàn)了一個（AC）14的微衛(wèi)星序列．

圖2 松江鱸MSTN基因的內(nèi)含子1序列

圖3 松江鱸MSTN基因的內(nèi)含子2序列

根據(jù)對其他物種MSTN基因序列的分析，發(fā)現(xiàn)魚類中存在幾種MSTN基因，包括MSTN1a，MSTN1b，MSTN2a，MSTN2b4種基因，且每個物種的MSTN基因都是含有兩個內(nèi)含子，三個外顯子，不同物種的內(nèi)含子長度不盡相同[6，8-10].表2為不同物種MSTN基因的內(nèi)含子長度，可以發(fā)現(xiàn)有的物種中內(nèi)含子1比內(nèi)含子2長，有的物種內(nèi)含子2比內(nèi)含子1長.斑點叉尾鮰的MSTN內(nèi)含子1最長為1748bp，瘤棘鲆MSTN內(nèi)含子1最短為244bp，斑馬魚MSTN-2的內(nèi)含子2最長為2201bp，尖吻鱸的MSTN內(nèi)含子2最短為188bp.而我們得到的松江鱸MSTN基因的內(nèi)含子2比內(nèi)含子1長．

2.3 不同地理群體松江鱸的線粒體D-loop序列比較

對分別來自長江、丹東、葫蘆島三個不同地理群體的松江鱸18個個體線粒體D-loop片段進行了PCR擴增和測序.用DNAMAN軟件對這18個個體進行序列比對后，共發(fā)現(xiàn)19個多態(tài)位點，但是未發(fā)現(xiàn)具有地理群體特征的單倍型.將這18個個體用MEGA軟件繪制N-J進化樹，也未發(fā)現(xiàn)三個地理群體松江鱸的特征性差異．

表2 不同物種MSTN基因內(nèi)含子的長度比較

3 討論

本研究根據(jù)松江鱸MSTN的cDNA序列設(shè)計引物，獲得了松江鱸MSTN基因的內(nèi)含子1和內(nèi)含子2序列，該序列符合真核基因內(nèi)含子的特征.在內(nèi)含子2中發(fā)現(xiàn)了1個（AC）14的重復(fù)序列.研究表明，已在多種哺乳動物的MSTN基因中發(fā)現(xiàn)微衛(wèi)星序列，De la Rosa-Reyna（2006）在牛MSTN內(nèi)含子1中的T-motif中發(fā)現(xiàn)了（TG）微衛(wèi)星[11]；姜運良（2004）在豬MSTN側(cè)翼區(qū)發(fā)現(xiàn)（TG）13的微衛(wèi)星標記[12].在魚類的MSTN基因中也發(fā)現(xiàn)了微衛(wèi)星序列：金鯛MSTN內(nèi)含子Ⅱ和3′-UTR區(qū)內(nèi)都存在(AC)n重復(fù)序列[4]，石鼓魚MSTN序列的3′-UTR區(qū)存在(AC)29微衛(wèi)星序列[13]；大黃魚肌肉生長抑制素基因的3′-UTR區(qū)存在（CA）n微衛(wèi)星序列，且該微衛(wèi)星座位與大黃魚體長和體重存在輕度的正相關(guān)[14]．

對物種進化的分子生物學研究目前主要集中在線粒體DNA方面，多數(shù)主要依賴于線粒體D-loop的信息，我們檢測到了松江鱸不同個體間的D-Loop區(qū)序列存在差別，但是三個地理群體未發(fā)現(xiàn)明顯的地理序列特征.可能是由于這些地理群體分隔時間短，沒有明顯大的差別.這一結(jié)果在對松江鱸MSTN基因序列多態(tài)性的檢測時一定程度上得到了驗證，在松江鱸MSTN基因的兩個內(nèi)含子和3′-UTR區(qū)未發(fā)現(xiàn)多態(tài)性位點（未發(fā)表資料），三個地理群體之間沒有明顯差異，一方面說明了MSTN基因比較保守，另一方面，可能也是由于這幾個地理群體的松江鱸分隔時間短.我們將進一步檢測其他線粒體基因中的變異情況，加入其他基因的信息，分析這些地理群體的松江鱸之間進化的關(guān)系．

[1]王金秋，潘連德，梁天紅，等.松江鱸魚(Trachidermus f asciatus)胚胎發(fā)育的初步觀察[J].復(fù)旦大學學報（自然科學版)，2004，43(2)：250-254．

[2]Mcpherron AC,Lawler AM,Lee SJ.Regulation of skeletal muscle mass in mice by a new TGF-beta superfamily membe[J].Nature,1997,387:83-90.

[3]趙浩斌，彭扣，王玉鳳，等.魚類肌肉生長抑制素研究進展[J].水生生物學報，2006，30（2）：227-231．

[4]Maccatrozzo L,Bargelloni L,Radaelli G,et al.Characterization of the myostatin gene in the gilthead seabream(Sparus aurata):sequence,genomic structure,and expression pattern[J].Mar Biotechnol(NY),2001,3(3):224-230.

[5]Garikipati DK,Gahr SA,Roalson EH,et al.Characterization of Rainbow Trout Myostatin-2 Genes(rtMSTN-2a and-2b):Genomic Organization,Differential Expression,and Pseudogenization[J].Endocrinology,2007,148(5):2106-2115.

[6]盧祥云，張營，王星果，等.松江鱸肌肉生長抑制素基因克隆和序列特征分析[J].動物學研究，2010，31(4)：387-394．

[7]郭新紅，劉少軍，劉巧，等.魚類線粒體DNA研究新進展[J].遺傳學報，2004，31（9）：983-1000．

[8]Maccatrozzo L,Bargelloni L,Cardazzo B,et al.A novel second myostatin gene is present in teleost fish[J].FEBS Lett,2001,509(1):36-40.

[9]Rodgers BD,Weber GM.Sequence conservation among fish myostatin orthologues and the characterization of two additional cDNA clones from Morone saxatilis and Morone americana[J].Comp Biochem Physiol B Biochem Mol Biol,2001,129(2-3):597-603.

[10]Rodgers BD,Garikipati DK.Clinical,agricultural,and evolutionary biology of myostatin:a comparative review[J].Endocr Rev,2008,29(5):513-534.

[11]De la Rosa-Reyna XF,Rodríguez Pérez MA,Sifuentes-Rincón AM.Microsatellite polymorphism in intron 1 of the bovine myostatin gene[J].J Appl Genet,2006,47(1):55-57.

[12]姜運良，李寧，趙興波，等.豬肌肉生長抑制素基因側(cè)翼區(qū)新微衛(wèi)星標記的鑒定及分析[J].遺傳學報，2004，31（5）:480-484.

[13]Maccatrozzo L,Bargelloni L,Patarnello P,et al.Characterization of the myostatin gene and a linked microsatellite marker in shi drum(Umbrina cirrosa,Sciaenidae)[J].Aquaculture,2002,205(122):49-60.

[14]薛良義，孫升，肖章奎，等.大黃魚肌肉生長抑制素基因微衛(wèi)星序列多態(tài)性分析[J].中國生物化學與分子生物學報，2008，24（10）：980-985．