高曉艷 董迎輝 姚韓韓 林志華（浙江萬里學(xué)院，浙江省水產(chǎn)種質(zhì)資源高效利用技術(shù)研究重點實驗室，寧波 315100）

文蛤HDAC1基因克隆、時空表達(dá)及生長相關(guān)SNP位點篩查

高曉艷 董迎輝 姚韓韓 林志華
（浙江萬里學(xué)院，浙江省水產(chǎn)種質(zhì)資源高效利用技術(shù)研究重點實驗室，寧波 315100）

摘要：為探索HDAC1基因在文蛤生長發(fā)育中的作用，研究通過已構(gòu)建的文蛤轉(zhuǎn)錄組文庫，利用SMART RACE技術(shù)擴(kuò)增得到文蛤HDAC1（Mm-HDAC1）基因的cDNA全長序列，分析了其生物信息學(xué)、組織及發(fā)育階段表達(dá)特征，并用直接測序法在外顯子中篩查生長相關(guān)的SNP位點。結(jié)果顯示，Mm-HDAC1基因的cDNA全長3065 bp，開放閱讀框1599 bp，編碼532個氨基酸； 氨基酸序列比對發(fā)現(xiàn)，文蛤與其他物種同源性為74.3%—78.7%。熒光定量PCR（qRT-PCR）結(jié)果表明，Mm-HDAC1基因在文蛤6個組織中均有表達(dá)，其中外套膜中的表達(dá)量相對最高，并與其他組織間有極顯著差異（P<0.01）； 不同發(fā)育時期的表達(dá)差異結(jié)果表明，Mm-HDAC1基因在殼頂幼蟲期表達(dá)量最高，顯著高于其他發(fā)育時期（P<0.05）。SNP位點篩查結(jié)果表明，在Mm-HDAC1基因的外顯子區(qū)域發(fā)現(xiàn)了19個SNP位點，其中有3個SNP位點（627A>T、924T>C和1266T>C）與文蛤生長相關(guān)。

關(guān)鍵詞：文蛤； 組蛋白去乙?；?； cDNA； 熒光定量； SNP

組蛋白去乙?；?（Histone deacetylase 1，HDAC1）是HDACs家族中的重要一員，在卵母細(xì)胞成熟、胚胎發(fā)育、神經(jīng)系統(tǒng)發(fā)生、骨骼分化、心臟發(fā)育、表皮分化、干細(xì)胞分化等過程中都發(fā)揮重要作用［1］。HDAC1是發(fā)現(xiàn)的第一個哺乳動物組蛋白脫乙酰化酶，推測其可能是基因轉(zhuǎn)錄的關(guān)鍵調(diào)節(jié)因子［2］，可介導(dǎo)位點特異DNA結(jié)合的轉(zhuǎn)錄抑制子［3］。在高等脊椎動物中，Sun等［4］純化并鑒定了雞HDAC1蛋白，劉峰等［5］、柳小春等［6］證明豬HDAC1基因與生長及胴體性狀相關(guān)，劉芳等［7］研究表明HDAC1與藏羚羊（Pantholops hodgsonii）的高原適應(yīng)性可能有關(guān)，Pillai等［8］證明HDAC1基因與斑馬魚（Danio rerio）顱面軟骨及魚鰭的形成有關(guān)，李慧等［9］證明HDAC1基因在牙鲆（Paralichthys olivaceus）變態(tài)前以及變態(tài)階段均有表達(dá)。而在海洋貝類中，僅在泥蚶（Tegillarca granosa）中見有基因克隆及時空表達(dá)研究的報道［10］。

文蛤（Meretrix meretrix）作為我國重要的經(jīng)濟(jì)養(yǎng)殖貝類之一，近些年來在分子遺傳和功能基因研究方面已經(jīng)取得了重要研究進(jìn)展。免疫相關(guān)基因的報道有血清補體C1q［11］、脂多糖誘導(dǎo)的腫瘤壞死因子（LITAF）［12］、銅鋅超氧化物歧化酶（SOB）［13］、金屬硫蛋白（MT）［14］、脂多糖和β-1，3葡聚糖結(jié)合蛋白（LGBP）［15］、富含半胱氨酸小腸蛋白（CRIP）［16］等；生長相關(guān)基因方面有鐵蛋白（Ferritin）［17］、組織蛋白酶B（Cathepsin B）［18］、磺基轉(zhuǎn)移酶（SULT）［19］和果糖-1，6-二磷酸醛縮酶（FBA）［20］等，但HDAC1基因卻未見報道。本研究首次克隆獲得文蛤HDAC1（Mm-HDAC1）的cDNA全長，研究了其在不同組織、不同發(fā)育階段中的表達(dá)特征，并在基因外顯子中篩選與生長相關(guān)的SNP位點，以期為文蛤生長相關(guān)基因資源的挖掘利用奠定基礎(chǔ)。

1 材料與方法

1.1 實驗材料

文蛤親貝取自寧波市鄞州區(qū)丹艷水產(chǎn)養(yǎng)殖基地，性腺發(fā)育成熟后進(jìn)行催產(chǎn)，通過隔離產(chǎn)卵和人工授精技術(shù)獲得同步發(fā)育的未受精卵、受精卵、4細(xì)胞胚胎、囊胚、原腸胚、擔(dān)輪幼蟲、D形幼蟲、殼頂幼蟲、眼點幼蟲和稚貝10個發(fā)育時期樣品，液氮速凍，存于-80℃冰箱，用于不同發(fā)育時期總RNA的提取。取3粒成貝，解剖取6個組織（外套膜、閉殼肌、內(nèi)臟團(tuán)、斧足、水管和鰓）樣品，用于不同組織總RNA提取。用于SNP位點篩選的文蛤自然群體共500粒，隨機(jī)取樣100粒，分別測量總重、殼寬、殼長和殼高等生長性狀，解剖取閉殼肌，液氮速凍，備用。

1.2 實驗方法

Mm-HDAC1基因cDNA全長克隆 用Trizol法提取文蛤閉殼肌總RNA，微量紫外分光光度儀檢測RNA濃度，1%瓊脂糖凝膠檢測RNA完整性。根據(jù)SMART RACE試劑盒（（Clontech）說明，合成cDNA第一鏈，用于全長cDNA 3′和 5′端序列快速擴(kuò)增（RACE）的模板。根據(jù)文蛤轉(zhuǎn)錄組文庫（登錄號:SRA021052）的注釋信息，檢索到Mm-HDAC1基因的EST片段，據(jù)此分別設(shè)計3′-RACE和5′-RACE特異性引物HDAC1-F和HDAC1-R（表 1），采用touch down PCR，分別進(jìn)行3′-和5′-RACE擴(kuò)增。PCR產(chǎn)物進(jìn)行割膠純化，純化產(chǎn)物與pMD18-T載體（TaKaRa）連接，并將重組質(zhì)粒轉(zhuǎn)化到大腸桿菌DH-5α感受態(tài)細(xì)胞（TIANGEN）中，涂平板，進(jìn)行藍(lán)白斑篩選，陽性克隆經(jīng)菌液PCR鑒定后送上海桑尼生物技術(shù)有限公司測序。

為保證Mm-HDAC1基因cDNA全長的準(zhǔn)確性，以拼接3′-和5′-RACE測序結(jié)果獲得的片段設(shè)計引物HDAC1-F1和HDAC1-R1，2×Taq MasterMix（北京康為世紀(jì)生物科技有限公司）進(jìn)行PCR擴(kuò)增，產(chǎn)物經(jīng)割膠純化，連接轉(zhuǎn)化后測序，具體步驟同上。

Mm-HDAC1基因序列特征分析 將5′-RACE片段和3′-RACE片段序列進(jìn)行拼接，用DNAMAN軟件預(yù)測氨基酸序列，并推測氨基酸基本理化性質(zhì)； ORF Finder程序搜索開放閱讀框； 信號肽區(qū)域用SignalP4.1（http://www.cbs.dtu.dk/services/ SignalP/）預(yù)測； NetPhos2.0工具（http//www.cbs.dtu.dk/services/NetPhos/）預(yù)測氨基酸磷酸化位點； Prot Scale軟件（http://web.expasy.org/protscale/）分析蛋白質(zhì)疏水性； TMHMM軟件（http://www.cbs.dtu.dk/services/TMHMM）預(yù)測蛋白質(zhì)跨膜區(qū)； SMART 4.0（http://smart.emblheidel-berg.de/）預(yù)測蛋白結(jié)構(gòu)域；Swiss Model軟件（http://swissmodel.expasy.org/ workspace/）預(yù)測蛋白質(zhì)高級結(jié)構(gòu)； MEGA 6.0軟件NJ法構(gòu)建系統(tǒng)進(jìn)化樹； Clustal X進(jìn)行氨基酸多序列比對。

表 1 實驗所用引物及其序列Tab.1 Primers used in the study

Mm-HDAC1基因不同組織、不同發(fā)育時期的表達(dá)差異分析 Trizol法提取上述文蛤6個組織和10個發(fā)育時期樣品的總RNA，反轉(zhuǎn)錄合成cDNA。由已獲得的cDNA全長序列設(shè)計特異性引物Real-H-F和Real-H-R（表 1），18S rRNA 基因作為內(nèi)參，利用ABI 7500Fast進(jìn)行qRT-PCR，每個樣本做3個平行。數(shù)據(jù)處理采用相對值2-ΔΔCt進(jìn)行計算，SPSS 20.0進(jìn)行統(tǒng)計分析。

Mm-HDAC1基因外顯子SNP位點篩選 提取隨機(jī)取樣的100個文蛤樣品的總RNA，反轉(zhuǎn)錄合成cDNA。根據(jù)獲得的Mm-HDAC1基因的cDNA全長，設(shè)計3對特異性PCR引物HDAC1-F2和HDAC1-R2、HDAC1-F3 和HDAC1-R3、 HDAC1-F4和HDAC1-R4，以反轉(zhuǎn)錄合成的100個文蛤樣本的cDNA為模板，進(jìn)行PCR擴(kuò)增，PCR產(chǎn)物直接送公司測序。結(jié)果用MEGA6.0進(jìn)行序列比對，查找SNP位點，按照突變位點相對于起始密碼子的位置分別命名，并用SPSS 20.0統(tǒng)計分析SNP位點與生長的相關(guān)性。

2 結(jié)果

2.1 Mm-HDAC1基因cDNA全長克隆

Mm-HDAC1基因cDNA全長3065 bp（GenBank登錄號:KP250877），開放閱讀框1599 bp，編碼532個氨基酸，5′UTR長35 bp，3′UTR為1431 bp，包含一個終止密碼子TAG，加尾信號AATAAA和26 bp的polyA尾巴。軟件預(yù)測顯示，Mm-HDAC1基因編碼蛋白質(zhì)的分子量為60.16 kD，理論等電點為5.55，表現(xiàn)為親水性，沒有明顯的信號肽和跨膜區(qū)域。

SMART軟件預(yù)測顯示，Mm-HDAC1蛋白含有一段HDAC1高度保守的氨基酸序列（79 Ile-311Cys）。蛋白質(zhì)高級結(jié)構(gòu)預(yù)測顯示，Mm-HDAC1蛋白由17個α-螺旋和17個β-折疊片組成，β-折疊片基本平行排列，這種排列方式可能對維持蛋白的穩(wěn)定性具有重要作用（圖 1）。NetPhos2.0軟件預(yù)測顯示，Mm-HDAC1蛋白有21個絲氨酸磷酸化位點，11個蘇氨酸磷酸化位點和9個酪氨酸磷酸化位點。

氨基酸多序列比對顯示，Mm-HDAC1基因與其他物種的同源性較高，為74.3%—78.7%，N端去乙?；富钚越Y(jié)構(gòu)域（54Try-325Ala）與其他物種的同源性則高達(dá)92%以上（表 2）。從構(gòu)建的系統(tǒng)進(jìn)化樹來看，文蛤與同為無脊椎動物的紫色球海膽（Strongylocentrotus purpuratus）親緣關(guān)系較近，脊椎動物中的哺乳綱、鳥綱、兩棲綱和魚綱聚成一大支，與文蛤親緣關(guān)系較遠(yuǎn)（圖 2）。

圖 1 文蛤Mm-HDAC1蛋白的二級結(jié)構(gòu)預(yù)測

2.2 Mm-HDAC1基因在文蛤不同組織和不同發(fā)育時期的表達(dá)差異分析

利用qRT-PCR技術(shù)，檢測了Mm-HDAC1基因在6個組織和10個發(fā)育時期的時空表達(dá)特征。結(jié)果表明，Mm-HDAC1基因在文蛤6個組織中均有表達(dá)，其中外套膜中的表達(dá)量相對最高，并與其他組織之間有極顯著性差異（P<0.01），足中的表達(dá)量僅次于外套膜，與其他組織之間也有極顯著差異（P<0.01）（圖3）。

Mm-HDAC1基因在文蛤幼蟲的10個發(fā)育時期均有表達(dá)，從未受精卵開始表量呈逐漸上升的趨勢，在殼頂幼蟲中達(dá)到最高，并顯著高于其他發(fā)育時期（P<0.05），到眼點幼蟲和稚貝時期表達(dá)量降低，并與其他時期之間存在顯著差異（P<0.05），但兩組之間并無顯著差異（P>0.05）（圖 4）。

2.3 Mm-HDAC1基因外顯子SNP位點篩查與分析

表 2 文蛤Mm-HDAC1氨基酸序列與其他物種HDAC1的同一性比對Tab.2 Comparison of amino acid sequences of Mm-HDAC1 of M.meretrix with others

經(jīng)序列比對，共篩查到19個SNP位點，將這些位點分別與生長性狀（總重、殼長、殼高、殼寬）進(jìn)行相關(guān)性分析。結(jié)果發(fā)現(xiàn)，3個位點627A>T、924T>C和1266T>C與生長性狀顯著相關(guān)，其中位點627A>T 的AT型個體總重和殼寬均顯著高于AA型個體（P<0.05），位點924T>C的CC型個體殼寬顯著高于TT型個體（P<0.05），位點1266T>C的TC型個體殼寬顯著高于TT型個體（P<0.05），位點231T>G、255A>T/G、588A>T/C、591A>G、1114A>C、1158A>G、1194C>T、1221T>C、1266T>C、1281A>G和1356A>G中，雜合子的各生長性狀均高于純合子的平均值，但兩組數(shù)據(jù)間無顯著差異（P>0.05）（表 3）。

圖 2 用MEGA6.0軟件鄰接法構(gòu)建的文蛤與其他物種HDAC1的系統(tǒng)進(jìn)化樹

圖 3 Mm-HDAC1基因在文蛤不同組織中的表達(dá)差異分析

圖 4 Mm-HDAC1基因在文蛤不同發(fā)育時期的表達(dá)差異分析

3 討論

HDAC1基因可通過調(diào)節(jié)染色質(zhì)結(jié)構(gòu)和抑制特異性轉(zhuǎn)錄因子活性來調(diào)控細(xì)胞生長和分化，故在動物各個組織中均有表達(dá)。本研究通過qRT-PCR技術(shù)檢測了Mm-HDAC1基因在文蛤各組織中的表達(dá)情況，結(jié)果表明，Mm-HDAC1基因在各組織中均有表達(dá)，其中外套膜中表達(dá)量最高，文蛤外套膜是由內(nèi)外兩層上皮細(xì)胞及其間的結(jié)締組織和肌肉纖維組成［21］，有研究表明HDAC1基因?qū)Σ溉閯游锷掀ぜ?xì)胞分化起到很重要的作用［2 2］，因此推斷Mm-HDAC1基因在文蛤外套膜中上皮細(xì)胞的分化中發(fā)揮作用。另外Mm-HDAC1在足中表達(dá)量也很高，僅次于外套膜，并與其他組織間有顯著差異。有研究顯示，組蛋白乙?；福℉ATs）和組蛋白去乙?；福℉DACs）以重塑染色質(zhì)的方式直接或間接影響成肌分化抗原（MyoD）或肌細(xì)胞增強(qiáng)因子2（MEF2）及其下游基因的表達(dá)，進(jìn)而調(diào)節(jié)肌肉的生肌過程［23，24］，而文蛤足中肌肉含量很高，導(dǎo)致HDAC1在該組織中大量表達(dá)。

HDAC1基因在多種動物胚胎發(fā)育中的表達(dá)情況都有研究，研究表明幾乎在所有物種的早期胚胎中都有表達(dá)。敲除HDAC1基因或HDAC1碳末端發(fā)生突變，都會導(dǎo)致小鼠早期胚胎發(fā)育異常，表現(xiàn)為胚胎生長阻滯、頭部和尿囊發(fā)育異常，胚胎在妊娠10.5d死亡［25，26］，說明HDAC1基因在生物的發(fā)育過程中發(fā)揮著重要的作用。本研究結(jié)果表明，Mm-HDAC1基因在各發(fā)育時期中均有表達(dá)，在未受精的卵子中就有表達(dá)，說明HDAC1基因是母源性的，Dufourcq等［27］在模式動物線蟲的早期胚胎發(fā)育過程中發(fā)現(xiàn)HDAC1基因是母源的，胚胎基因組開始轉(zhuǎn)錄時就能檢測到。而在研究動物后期胚胎發(fā)育的過程中，發(fā)現(xiàn)HDAC1基因在果蠅、斑馬魚、兩棲動物的頭部區(qū)域，包括中樞神經(jīng)系統(tǒng)和大腦發(fā)育過程中都有所表達(dá)［28—30］。Mm-HDAC1基因在文蛤殼頂幼蟲時期表達(dá)量達(dá)到最高，可能是殼頂幼蟲時期是文蛤快速生長的時期，細(xì)胞分裂和組織器官分化速度加快，生命活動旺盛［31］，需要更多HDAC1基因來調(diào)節(jié)胚胎發(fā)育過程中的基因轉(zhuǎn)錄。到眼點幼蟲和稚貝時期，Mm-HDAC1的表達(dá)量又逐漸下降，可能是由于此時文蛤各組織、器官已基本形成，Mm-HDAC1基因的表達(dá)量逐漸達(dá)到穩(wěn)定。

表 3 不同基因型文蛤的生長性狀Tab.3 Growth traits of M.meretrix with different genotypes

續(xù)表

本研究利用直接測序的方法在Mm-HDAC1基因的外顯子區(qū)域共篩查到19個SNP位點，3個位點（627A>T、924T>C 和1266T>C）與生長性狀有相關(guān)性，其中位點627A>T 的AT型個體總重、殼寬和殼高均顯著高于AA型個體（P<0.05），位點924T>C的CC型個體殼寬顯著高于TT型個體（P<0.05），位點1266T>C的TC型個體殼寬顯著高于TT型個體（P<0.05），提示HDAC1基因可能是貝類生長相關(guān)的潛在基因。有研究表明，在杜洛克豬種的HDAC1基因中，GG基因型的生長速度和生長性能指數(shù)明顯優(yōu)于AA型，等位基因G對生長具有正效應(yīng)［5］。進(jìn)一步分析發(fā)現(xiàn)，位點627A>T和924T>C位于Mm-HDAC1的乙?；富钚越Y(jié)構(gòu)區(qū)域，該位點的突變可能與HDAC1發(fā)揮其活性有關(guān)，但其具體機(jī)制還需進(jìn)一步深入研究。

參 考 文 獻(xiàn)：

［1］Wang Z W.Function of HDAC1gene during the development of bovine oocytes and cloned embryos［D］.Thesis for Doctor of Science，Jilin University，Changchun.2011［王中偉.HDAC1基因在牛卵母細(xì)胞和早期克隆胚胎發(fā)育中的作用.博士學(xué)位論文，吉林大學(xué)，長春.2011］

［2］Taunton J，Hassig C A，Schreiber S L.A mammalian histone deacetylase related to the yeast transcriptional regulator Rpd3p［J］.Science，1996，272（8）:408—411

［3］Hassiga C A，F(xiàn)leischer T C，Billin AN，et al.Histone deacety-lase activity is required for full transcriptional repression by mSin3A［J］.Cell，1997，89（3）:341

［4］Sun J M，Chen H Y，Moniwa M，et al.Purification and characterization of chicken erythrocyte histone deacetylase 1［J］.Biochemistry，1999，38（18）:5939—5947

［5］Liu F.Study of on effects of HDAC1 gene and PIT-1 gene on partial economic traits in pigs［D］.Thesis for Master of Science，Hunan Agricultural University，Changsha.2005［劉峰.HDAC1基因、PIT-1基因與豬部分經(jīng)濟(jì)性狀相關(guān)性研究.碩士學(xué)位論文，湖南農(nóng)業(yè)大學(xué)，長沙.2005］

［6］Liu X C，Xiao T Y，He W G，et al.Genetic effect of three loci related with growth and carcass traits in swine［J］.Scientia Agricultura Sinica，2009，42（2）:742—747［柳小春，肖調(diào)義，何偉光，等.豬生長及胴體性狀3個相關(guān)基因座遺傳效應(yīng).中國農(nóng)業(yè)科學(xué)，2009，42（2）:742—747］

［7］Liu F，Ge R L.Cloning and sequencing of HDAC1 coding cDNA sequence from Tibetan antelope（Pantholops hodgsonii）［J］.Acta Theriologica Sinica，2011，31（4）:396—403［劉芳，格日力.藏羚羊組蛋白去乙?；?基因編碼區(qū)的克隆與序列分析.獸類學(xué)報，2011，31（4）:396—403］

［8］Pillai R，Coverdale L E，Dubey G，et al.Histone deacetylase 1（HDAC-1）required for the normal formation of craniofacial cartilage and pectoral fins of the zebrafish［J］.Developmental Dynamics，2004，231（3）:647—654

［9］Li H，Xu Y P，Bao Y P.Different expression pattern of HDAC1 before and during the process of metamorphosis in Paralichthys olivaceus［J］.Journal of Shanghai Ocean University，2014，23（5）:641—649［李慧，徐義平，鮑寶龍.牙鲆胚后發(fā)育階段HDAC1基因的空間表達(dá).上海海洋大學(xué)學(xué)報，2014，23（5）:641—649］

［10］Zhou X L.Development and application of SSR and cloning of growth-related gene HDAC1 in Tegillarca granosa［D］.Thesis for Master of Science，Shanghai Ocean University，Shanghai.2013［周小龍.泥蚶微衛(wèi)星分子標(biāo)記的開發(fā)與應(yīng)用及生長相關(guān)基因HDAC1的克隆.碩士學(xué)位論文，上海海洋大學(xué)，上海.2013］

［11］Sui L J，Liu W D，Li Y F，et al.Molecular cloning，mRNA expression and bioassay of recombinant protein of C1q gene in hard clam（Meretrix meretrix）［J］.Fisheries Science，2012，31（6）:311—319

［12］Li H J，Yang Q，Gao X G，et al.Identification and expression of a putative LPS-induced TNF-a factor from Asiatic hard clam Meretrix meretrix［J］.Molecular BiologyReports，2012，39（2）:865—871

［13］Zhu D.Molecular cloning and mRNA expression and recombinant protein of an Cu/Zn-superoxide dismutase（Cu/Zn-SOD）gene from hard clam（Meretrix meretrix）［D］.Thesis for Master of Science，Liaoning Normal University，Dalian.2011［朱丹.文蛤胞內(nèi)銅鋅超氧化物歧化酶的基因克隆與表達(dá)及重組蛋白活性的測定.碩士學(xué)位論文，遼寧師范大學(xué)，大連.2011］

［14］Gao X G，He C B，Li Y F，et al.Analysis on cloning and sequence Characteristics of metallothionein gene in hard clam Meretrix meretrix［J］.Chinese Journal Fisheries，2009，22（4）:8—11［高祥剛，赫崇波，李云峰，等.文蛤金屬硫蛋白基因的克隆與分析.水產(chǎn)學(xué)雜志，2009，22（4）:8—11］

［15］Liu S X，Qi Z H，Zhang J J，et al.Lipopolysaccharide and β-1，3-glucan binding protein in the hard clam（Meretrix meretrix）:molecular characterization and expression analysis［J］.Genetics and Molecular Research，2014，13（3）:4956—4966

［16］Chen H，Yang X，Tang T，et al.The involvement of cysteine-rich intestinal protein in early development and innate immunity of Asiatic hard clam，Meretrix meretrix［J］.Fish and Shellfish Immunology，2014，40（2）:435—440

［17］Wang X M.Molecular cloning and functional analysis of genes involved in larval development in clam Meretrix meretrix［D］.Chinese Academy of Sciences，Institute of Oceanology，Qingdao.2008［王曉梅.文蛤（Meretrix meretrix）幼蟲生長發(fā)育相關(guān)基因的克隆和功能分析.博士學(xué)位論文，中國科學(xué)院海洋研究所，青島.2008］

［18］Yao X L，Zhang J Q，Sun J S，et al.Recombinant expression，characterization and expressional analysis of clam Meretrix meretrix cathepsin B，an enzyme involved in nutrient digestion［J］.Molecular Biology Reports，2011，38（3）:1861—1868

［19］Wang C，You Y N，Wang H X，et al.Genetic diversity of the sulfotransferase-like gene and one nonsynonymous SNP associated with growth traits of clam，Meretrix meretrix［J］.Molecular Biology Reports，2012，39（2）:1323—1331

［20］Wang C，Wang H X，Li Y，et al.Identification of a fructose-1，6-bisphosphate aldolase gene and association of the single nucleotide polymorphisms with growth traits in the clam Meretrix meretrix［J］.Molecular Biology Reports，2012，39（4）:5017—5024

［21］Cheng L，Wang K L.Studies on the morphology and ultrastructure of Meretrix meretrix Linnaeus［J］.Zoological Research，1997，18（4）:351—357［程量，王康樂.文蛤形態(tài)學(xué)及其亞顯微結(jié)構(gòu)的研究.動物學(xué)研究，1997，18（4）:351—357］

［22］Tou L，Liu Q，Shivdasani R A.Regulation of mammalian epithelial differentiation and intestine development by class I histone deacetylases［J］.Molecular and Cellular Biology，2004，24（8）:3132—3139

［23］Puri P L，Iezzi S，Stiegler P，et al.Class I histone deacetylases sequentially interact with MyoD and pRb during skeletal myogenesis［J］.Molecular Cell，2001，8（4）:885—897

［24］Mal A，Harter M L.MyoD is functionally linked to the silencing of a muscle-specific regulatory gene prior to skeletal myogenesis［J］.Proceedings of the National Academy of Sciences of the USA，2003，100（4）:1735—1739

［25］Lagger G，O'Carroll D，Rembold M，et al.Essential function of histone deacetylase 1 in proliferation control and CDK inhibitor repression［J］.The EMBO Journal，2002，21（11）:2672—2681

［26］Montgomery R L，Davis C A，Potthoff M J，et al.Histone deacetylases 1 and 2 redundantly regulate cardiac morphogenesis，growth，and contractility［J］.Genes and Development，2007，21（14）:1790—1802.

［27］Dufourcq P，Victor M，Gay F，et al.Function requirement for histone deacetylase 1 in Caenorhabditis elegans gonadogenesis［J］.Molecular Cellular Biology，2002，22（9）:3024—3034

［28］Mannerrvik M，Levine M.The Rpd3 histone deacetylaseis required for segmentation of the Drosophila embryo［J］.Proceedings of the National Academy of Sciences of the USA，1999，96（12）:6797—6801

［29］Nambiar R M，Henion P D.Sequential antagonism of early and late Wnt-signaling by zebrafish colgate promotes dorsal and anterior fates［J］.Developmental Biology，2004，267（1）:165—180

［30］Damjanovski S，Sachs L M，Shi Y B.Multiple stage-dependent roles for histone deacetylases during amphibian embryogenesis:implications for the Involvement of extracellular matrix remodeling［J］.International Journal of Developmental Biology，2000，44（7）:769—776

［31］Wang R C，Wang Z P.Science of Marine Shellfish Culture［M］.Qingdao:China Ocean University Press.2008，28—29［王如才，王昭萍.海水貝類養(yǎng)殖學(xué).青島:中國海洋大學(xué)出版社.2008，28—29］

CLONING，SPATIOTEMPORAL EXPRESSION AND SNPS IDENTIFICATION OF HDAC1 GENE IN HARD CLAM MERETRIX MERETRIX

GAO Xiao-Yan，DONG Ying-Hui，YAO Han-Han and LIN Zhi-Hua
（Key Laboratory of Aquatic Germplasm Resource of Zhejiang，Zhejiang Wanli University，Ningbo 315100，China）

Abstract:Mm-HDAC1 cDNA was cloned by SMART RACE techniques using 454 cDNA library of Meretrix meretrix.The bioinformatics and expression profiles in different tissues and developmental stages were analyzed，and SNPs were screened.The results indicated that the full length cDNA of Mm-HDAC1 gene was 3065 bp，containing a complete 1599 bp ORF that encodes 532 amino acids.Comparison of animal acid sequence，M.meretrix shared 74.3%—78.5% identity with others，suggesting that Mm-HDAC1 was highly conservative.The result of qRT-PCR showed that Mm-HDAC1 expressed in all six tissues，and its expression was significantly higher in mantle than other tissues（P<0.01）.Moreover，the expression of Mm-HDAC1 gradually increased with the process of the development with the highest level in umbo larvae stage（P<0.01）.19 SNPs in exons of Mm-HDAC1 were identified and 3 SNPs were potentially associated with M.meretrix growth（627A>T，924T>C and 1266T>C）.

Key words:Meretrix meretrix； Histone deacetylase1（HDAC1）； cDNA； qRT-PCR； SNP

中圖分類號：Q785； S968.3

文獻(xiàn)標(biāo)識碼：A

文章編號：1000-3207（2016）03-0601-08

doi:10.7541/2016.81

收稿日期：2015-07-14；

修訂日期：2015-12-20

基金項目：國家自然科學(xué)基金項目（31372527）； 國家現(xiàn)代貝類產(chǎn)業(yè)技術(shù)體系項目（CARS-48）； 浙江省重大科技專項（2012C12907-4）； 國家水產(chǎn)種質(zhì)資源平臺項目（2015DKA30470）資助［Suppprted by National Natural Science Foundation of China（31372527）； Modern Agro-industry Technology Research Systerm（CARS-48）； Zhejiang Major Program of Science and Technology（2012C12907-4）and National Infrastructure of Fishery Germplasm Resource Programme（2015DKA30470）］

作者簡介：高曉艷（1989—），女，河南安陽人； 碩士研究生； 主要從事海洋貝類遺傳育種研究。E-mail:yanyuninan6264@163.com

通信作者：董迎輝（1980—），副研究員； 碩導(dǎo)； 主要從事貝類遺傳育種研究。E-mail:dongyinghui118@126.com； 林志華（1965—），研究員；博導(dǎo)； E-mail:zhihua9988@126.com