• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    半滑舌鰨(Cynoglossus semilaevis) 2種視黃酸受體RARα和RARγ克隆及組織表達(dá)特性*

    2018-12-19 08:28:58宋雪松徐永江柳學(xué)周劉永山張雅星
    漁業(yè)科學(xué)進(jìn)展 2018年6期
    關(guān)鍵詞:黑化舌鰨牙鲆

    宋雪松 徐永江 柳學(xué)周 史 寶 王 濱 劉永山 張雅星

    ?

    半滑舌鰨() 2種視黃酸受體RARα和RARγ克隆及組織表達(dá)特性*

    宋雪松1,2徐永江1柳學(xué)周1①史 寶1王 濱1劉永山1,2張雅星1,2

    (1. 青島海洋科學(xué)與技術(shù)試點(diǎn)國家實(shí)驗(yàn)室海洋漁業(yè)科學(xué)與食物產(chǎn)出過程功能實(shí)驗(yàn)室 中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所 青島 266071;2. 上海海洋大學(xué)水產(chǎn)與生命學(xué)院 上海 201306)

    利用RT-PCR和RACE方法獲得了半滑舌鰨() 2種視黃酸受體RARα和RARγ的cDNA全長序列,并采用定量PCR技術(shù)分析了其組織表達(dá)特性。半滑舌鰨RARα cDNA序列全長為1823 bp,編碼443個(gè)氨基酸;RARγ cDNA序列全長為1959 bp,編碼489個(gè)氨基酸。同源性分析顯示,半滑舌鰨RARα和RARγ同源性高達(dá)60.8%,與牙鲆()同源性高達(dá)97.0%,具有較強(qiáng)的進(jìn)化保守性。系統(tǒng)進(jìn)化分析顯示,半滑舌鰨RARα和RARγ分別歸屬于單獨(dú)的分支,且與其他魚類聚合成簇。組織表達(dá)分析顯示,RARα mRNA在腎中表達(dá)量最高,而RARγ mRNA在脾中表達(dá)量最高,RARα和RARγ mRNA在其他組織中均有表達(dá),表明半滑舌鰨2種RAR都可能參與多種生理過程調(diào)控。半滑舌鰨2種RAR mRNA在有眼側(cè)皮膚、無眼側(cè)黑化皮膚和無眼側(cè)正常皮膚中的表達(dá)量依次升高,且發(fā)現(xiàn)RARα在正常有眼側(cè)皮膚的表達(dá)高于RARγ,而RARγ在無眼側(cè)正常皮膚中的表達(dá)顯著高于RARα,無眼側(cè)黑化皮膚中RARα表達(dá)高于RARγ。RAR基因在有眼側(cè)和無眼側(cè)皮膚組織中的差異表達(dá)可能和RA/RAR系統(tǒng)調(diào)節(jié)體色有關(guān)。

    半滑舌鰨;視黃酸受體;基因克??;表達(dá)模式;無眼側(cè)黑化

    視黃酸(Retinoic acid, RA)是一種脂溶性小分子物質(zhì),又稱維甲酸,是維生素A(VA)在機(jī)體內(nèi)經(jīng)過一系列水解和酶催化作用生成的最終代謝產(chǎn)物。視黃酸合成進(jìn)入細(xì)胞后先與重組人胞內(nèi)維甲酸結(jié)合蛋白-1 (CRABP1)和CRΑBP2蛋白結(jié)合,轉(zhuǎn)運(yùn)進(jìn)入細(xì)胞核,形成核受體二聚體,調(diào)節(jié)相關(guān)基因的表達(dá)(Budhu, 2002)。RA的作用是通過定位在靶細(xì)胞核內(nèi)的特異性受體介導(dǎo)的。目前,已發(fā)現(xiàn)兩類視黃酸受體:視黃酸受體(RAR)和視黃酸X受體(RXR),作為配體激活的轉(zhuǎn)錄因子(Chambon, 1996)。RAR是一類細(xì)胞核受體,屬于類固醇和甲狀腺激素受體超家族,其通過與RXR形成異源二聚體(RAR/RXR),增加與視黃素應(yīng)答元件的親和力,由此加強(qiáng)本身的轉(zhuǎn)錄活性和對(duì)配體的敏感性(Kliewer, 1992; Hallenbeck, 1992)。哺乳動(dòng)物RAR有3種亞型:RARα、RARβ和RARγ,配體為順式視黃酸(9-cis-RA)和反式視黃酸(ATRA)(Brand, 1988),但魚類中目前沒有發(fā)現(xiàn)RARβ亞型的存在(Jones, 1995)。研究發(fā)現(xiàn),RA/RAR系統(tǒng)具有抑制腫瘤生長、保護(hù)上皮組織、促進(jìn)胚胎骨骼心臟功能和視覺發(fā)育功能(De Luca, 1991; Zhou, 2017; Isojima, 2014),并且在抑制細(xì)胞凋亡(Herget,1998)和免疫調(diào)節(jié)(Erkelens, 2017)也有重要生理作用,是機(jī)體維持穩(wěn)態(tài)的重要因子。

    RA及其受體RAR系統(tǒng)在魚類中的研究較少(Faehnrich, 2016; Zhang, 2013; 冷向軍, 2017)。在鲆鰈類中,RA/RAR系統(tǒng)對(duì)牙鲆()骨骼發(fā)育和眼睛位移等變態(tài)發(fā)育過程(Haga, 2002; Martinez, 2007),對(duì)塞內(nèi)加爾鰨()(Ignacio, 2009)和大西洋庸鰈()(Lewis-McCrea, 2010)骨骼發(fā)育畸形都具有重要的生理調(diào)控作用。對(duì)牙鲆的研究表明,利用9-cis-RA或者ATRA處理變態(tài)前仔魚,可以通過RA/RAR信號(hào)通路的介導(dǎo)作用抑制魚苗眼睛完成位移。同時(shí)發(fā)現(xiàn),牙鲆和半滑舌鰨()皮膚的RA濃度在有眼側(cè)皮膚高于無眼側(cè)皮膚,且光照可通過RA/RAR信號(hào)通路誘導(dǎo)牙鲆無眼側(cè)色素過度沉著(黑化),表明RA/RAR信號(hào)通路在半滑舌鰨、牙鲆的體色左右不對(duì)稱發(fā)育模式中起重要調(diào)控功能(Shao, 2017),但RA/RAR系統(tǒng)對(duì)半滑舌鰨體色調(diào)控的具體機(jī)制仍待進(jìn)一步研究。

    半滑舌鰨屬鰈形目、舌鰨科、舌鰨屬,為我國近海自然分布的重要經(jīng)濟(jì)魚類,自人工繁育技術(shù)取得突破以來,其養(yǎng)殖產(chǎn)業(yè)得到快速發(fā)展,已成為鲆鰈類三大主導(dǎo)養(yǎng)殖品種之一(鄧景耀等, 1988; 柳學(xué)周等, 2006、2014)。養(yǎng)殖生產(chǎn)中發(fā)現(xiàn),半滑舌鰨無眼側(cè)黑化(色素沉積過多)問題日益凸顯,而無眼側(cè)發(fā)生黑化的商品魚市場價(jià)格較無眼側(cè)正常魚價(jià)格低20%以上,成為制約半滑舌鰨養(yǎng)殖產(chǎn)業(yè)經(jīng)濟(jì)效益的瓶頸之一。目前,國內(nèi)外對(duì)半滑舌鰨無眼側(cè)黑化發(fā)生的相關(guān)機(jī)制鮮有報(bào)道,本實(shí)驗(yàn)室前期研究了半滑舌鰨體色相關(guān)功能基因POMC、MCHR與MCH等的克隆、表達(dá)調(diào)控及其與體色的關(guān)系(史學(xué)營等, 2015、2017; 朱學(xué)武等, 2016; 徐永江等, 2017),為認(rèn)識(shí)半滑舌鰨無眼側(cè)體色調(diào)控機(jī)制積累了資料。本研究擬開展半滑舌鰨RAR結(jié)構(gòu)及其表達(dá)特性研究,以期為探究RA/RAR系統(tǒng)在養(yǎng)殖半滑舌鰨無眼側(cè)黑化調(diào)控中的作用機(jī)制提供基礎(chǔ)資料。

    1 材料與方法

    1.1 實(shí)驗(yàn)用魚及樣品處理

    實(shí)驗(yàn)用半滑舌鰨于2017年6~8月取自山東省海陽市黃海水產(chǎn)有限公司。取樣實(shí)驗(yàn)魚3尾(都存在一定程度無眼側(cè)黑化),體長為(33±3) cm,體重為(237±30) g,用于RAR基因克隆與組織表達(dá)特性分析。實(shí)驗(yàn)魚以MS-222 (280 mg/L)麻醉后,快速取性腺、肝臟、心臟、胃、腸、脾、腎、垂體、腦、鰓、肌肉、有眼側(cè)正常皮膚、無眼側(cè)黑化皮膚和無眼側(cè)正常皮膚組織投入液氮速凍后,轉(zhuǎn)入–80℃保存,用于總RNA的提取。

    1.2 總RNA提取和cDNA第1鏈合成

    利用RNAiso Plus (TaKaRa,日本)試劑盒并按照操作說明提取各組織樣品總RNA,通過1%瓊脂糖凝膠電泳檢測RNA質(zhì)量,NanoDrop 2000 (Thermo,美國)測定RNA濃度。取適量鰓組織總RNA,以PrimeScriptTMⅡ1st strand cDNA Synthesis Kit (TaKaRa,日本)合成cDNA第1鏈。以SMARTerTMRACE cDNA Amplification Kit(Clontech,美國)合成5¢-RACE及3¢-RACE cDNA第1鏈,用于RAR基因RACE全長克隆。取等量各組織樣品的總RNA,用PrimeScript RT Reagent Kit with gDNA Eraser反轉(zhuǎn)錄試劑盒(TaKaRa,日本)合成cDNA第1鏈,用于RAR mRNA組織表達(dá)特性及分析。各操作步驟均嚴(yán)格按照使用說明書進(jìn)行。

    1.3 總RNA提取和中間片段擴(kuò)增

    根據(jù)GenBank登記的XM_017034299.1、XM_ 008318896.2預(yù)測半滑舌鰨RARs序列保守區(qū)設(shè)計(jì)特異引物(表1),以肝臟組織為模板,擴(kuò)增RARα基因的核心序列,PCR反應(yīng)體系(25 μl):0.2 μl酶、2.5 μl 10×PCR Buffer、2 μl dNTP Mixture、0.5 μl模板、1 μl RARα-F、1 μl RARα-R、17.8 μl ddH2O。反應(yīng)條件:94℃ 5 min;94℃ 30 s,55℃ 30 s,72℃ 2 min,38個(gè)循環(huán);72℃延伸10 min。擴(kuò)增RARγ基因的模板為脾臟,PCR反應(yīng)體系和條件同RARα;PCR產(chǎn)物經(jīng)1.2%瓊脂糖凝膠電泳分離后,切膠回收目的條帶并純化?;厥誔CR產(chǎn)物與pEASY-T1載體(北京全式金生物公司)連接,轉(zhuǎn)化至Trans1-T1感受態(tài)細(xì)胞(全式金),LB固體培養(yǎng)基37℃培養(yǎng)過夜,挑取陽性克隆送至生工生物工程(上海)股份有限公司測序;RAR的中間序列已上傳NCBI數(shù)據(jù)庫(RARα登錄號(hào):MG596268;RARγ登錄號(hào):MG596269)。

    表1 半滑舌鰨RAR基因克隆使用的PCR擴(kuò)增引物

    Tab.1 Primers used for PCR amplification of RAR of C. semilaevis

    1.4 RAR的RACE擴(kuò)增

    根據(jù)克隆驗(yàn)證的中心片段設(shè)計(jì)RACE引物。用Smart RACE Advantage 2 PCR試劑盒(Clontech,美國)進(jìn)行梯度PCR擴(kuò)增。第1次PCR,反應(yīng)體系:17 μl ddH2O、2.5 μl Buffer、2 μl 50×dNTP Mix、0.5 μl 50×Advantage 2 Polymerase Mix、1 μl cDNA、引物RARα1 1 μl和1 μl UPM,共計(jì)25 μl。設(shè)計(jì)Touchdown PCR,反應(yīng)條件為94℃ 30 s,66℃30 s,72℃ 2 min,15個(gè)循環(huán),T值每5個(gè)循環(huán)降低2℃;94℃ 30 s,55℃ 30 s,72℃ 2 min,20個(gè)循環(huán);最后72℃延伸5 min。以第1次PCR產(chǎn)物稀釋10倍為模板,進(jìn)行巢式PCR,反應(yīng)體系:17 μl ddH2O、2.5 μl Buffer、2 μl 50×dNTP Mix、0.5 μl 50×Advantage 2 Polymerase Mix、1 μl cDNA、1 μl NPM和1 μl RARα2引物,共計(jì)25 μl。PCR反應(yīng)條件同第1次PCR。PCR產(chǎn)物于1.2%瓊脂糖凝膠電泳檢測后,對(duì)目的條帶進(jìn)行膠回收、載體連接、轉(zhuǎn)化、篩選陽性克隆并測序。RARγ反應(yīng)體系與條件同RARα。

    1.5 RAR mRNA定量表達(dá)分析

    根據(jù)獲得的半滑舌鰨RARα和RARγ的cDNA序列設(shè)計(jì)定量PCR引物qRARα和qRARγ(表1),以18S為內(nèi)參。利用Mastercycler ep realplexreal-time PCR儀(Eppendorf,德國),使用SYBR Premix ExTMⅡ試劑盒(TaKaRa)進(jìn)行定量擴(kuò)增,PCR體系(20 μl):1 μl cDNA模板、上下游引物各0.8μl (10 μmol/L)、10μl SYBR Premix ExTMⅡ和7.4 μl dd H2O.。PCR反應(yīng)條件:95℃預(yù)變性30 s,95℃5 s,58℃20 s、共40個(gè)循環(huán)。每個(gè)樣品測試設(shè)置3個(gè)重復(fù)。RAR mRNA的表達(dá)量以18S mRNA表達(dá)量為基礎(chǔ),利用2-DD方法計(jì)算獲得(Livak, 2001)。

    1.6 序列分析及數(shù)據(jù)處理

    半滑舌鰨RAR基因的結(jié)構(gòu)、分子量預(yù)測、等電點(diǎn)預(yù)測使用ExPASy在線數(shù)據(jù)庫預(yù)測(www.expasy. org/tools/protparam.html);氨基酸序列推導(dǎo)、序列拼接和氨基酸同源性分析均使用軟件DNAMAN 6.0,信號(hào)肽預(yù)測使用SignalP 4.1 (http://www.cbs.dtu.dk/ services/SignalP/)。亞細(xì)胞定位使用PSORT Ⅱ軟件(https://psort.hgc.jp/form2.html)。結(jié)構(gòu)域預(yù)測使用NCBI數(shù)據(jù)庫(https://www.ncbi.nlm.nih.gov/Structure/ cdd/wrpsb.cgi),氨基酸序列比對(duì)和系統(tǒng)進(jìn)化分析使用ClustalW在線軟件(http://www.genome.jp/tools-bin/ clustalw)和MEGA 7軟件。通過SOPMA軟件(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html)分析蛋白質(zhì)二級(jí)結(jié)構(gòu),通過SWISS-MODEL在線軟件(http://www.swi-ssmodel.expasy.org/)分析預(yù)測蛋白質(zhì)三級(jí)結(jié)構(gòu)。

    實(shí)驗(yàn)數(shù)據(jù)均以平均值±標(biāo)準(zhǔn)差(Mean±SD)表示,多組數(shù)據(jù)間比較采用SPSS 19.0統(tǒng)計(jì)軟件進(jìn)行單因素方差分析(One-way ANOVA)、Duncan和SNK多重比較分析,當(dāng)<0.05時(shí)表示差異顯著。

    2 結(jié)果

    2.1 RAR cDNA序列結(jié)構(gòu)

    半滑舌鰨RARα cDNA序列全長為1823 bp (圖1),包括17 bp的5¢非編碼區(qū)(UTR)、1332 bp的開放閱讀框(ORF)和474 bp的3¢非編碼區(qū)(UTR),編碼443個(gè)氨基酸,預(yù)測編碼蛋白分子量為49 kb,等電點(diǎn)為8.47。半滑舌鰨RARγ cDNA序列全長為1959 bp (圖2),包括305 bp的5¢(UTR)、1497 bp的ORF和98 bp的3¢(UTR),編碼498個(gè)氨基酸,預(yù)測編碼蛋白的分子量為55.8 kb,等電點(diǎn)為4.99。預(yù)測2種RAR基因的亞細(xì)胞定位均位于細(xì)胞核。

    2.2 RAR編碼蛋白質(zhì)空間結(jié)構(gòu)預(yù)測

    通過SOPMA軟件分析RAR編碼蛋白的空間二級(jí)結(jié)構(gòu):在RARα成熟蛋白的二級(jí)結(jié)構(gòu)中,α-螺旋占34.31%,β-轉(zhuǎn)角占8.35%,無規(guī)則卷曲占43.12%,延伸鏈占14.22%。在RARγ成熟蛋白的二級(jí)結(jié)構(gòu)中,α-螺旋占36.75%,β-轉(zhuǎn)角占9.44%,無規(guī)則卷曲占40.36%,延伸鏈占13.45%。通過SWISS-MODEL網(wǎng)站同源建模方法構(gòu)建了半滑舌鰨RAR編碼的蛋白質(zhì)可能的三級(jí)結(jié)構(gòu),以同源建模的方法與各自的50個(gè)RA核受體家族模板構(gòu)建預(yù)測可能的蛋白質(zhì)三級(jí)結(jié)構(gòu),RARα獲得6種三級(jí)結(jié)構(gòu)模型,選取GMQE評(píng)價(jià)值0.55、QMEAN穩(wěn)定系數(shù)-2.79(負(fù)值越大越穩(wěn)定)構(gòu)建模型;RARγ有5種,選取GMQE評(píng)價(jià)值0.61,QMEAN穩(wěn)定系數(shù)-2.13構(gòu)建模型。三級(jí)結(jié)構(gòu)模型顯示(圖3),RAR蛋白質(zhì)分為3部分,藍(lán)色的DBD區(qū)和LBD區(qū)通過橘黃色的D區(qū)鉸鏈區(qū)相連。

    2.3 RAR的氨基酸序列同源性比較

    同源性分析顯示,半滑舌鰨RARα的氨基酸序列與同屬鰈形目的牙鲆的同源性最高,為97.0%,其次為銀大麻哈魚()92.2%、鱸魚() 87.9%、斑馬魚()86.8%、紅鰭東方鲀()64.6%,且與兩棲類、爬行類、嚙齒類、鳥類和人()的相似度分別為79.6%、82.2%、80.5%、81.5%和80.7%。半滑舌鰨RARγ的氨基酸序列同樣與牙鲆同源性最高,達(dá)97%,與其他魚類的同源性均達(dá)90.5%以上,與兩棲類、嚙齒類和人的相似度分別為71.4%、84.5%和80.7%。另外,半滑舌鰨RARα和RARγ的氨基酸序列相似度為60.8% (表2)。

    通過結(jié)構(gòu)域預(yù)測發(fā)現(xiàn),半滑舌鰨與鱸魚、人、斑馬魚具有相似的功能結(jié)構(gòu)域,都具有DNA結(jié)合區(qū)(DBD區(qū))和配體結(jié)合區(qū)(LBD)區(qū)。利用ClustalW對(duì)半滑舌鰨RAR的氨基酸序列與其他物種的RAR氨基酸序列進(jìn)行了比較(圖4)。結(jié)果發(fā)現(xiàn),半滑舌鰨與其他魚類RAR的氨基酸序列整體保守性較強(qiáng),除在N端的A/B區(qū)和C端的F區(qū)末尾保守型較差外,中間的DNA與配體結(jié)合區(qū)保守度較高。

    2.4 RARs系統(tǒng)進(jìn)化分析

    利用NJ法構(gòu)建了基于氨基酸序列的半滑舌鰨RARα、RARγ和其他脊椎動(dòng)物的系統(tǒng)進(jìn)化樹(圖5),半滑舌鰨RARα和RARγ都分別與鰈形目、鱸形目、鯉形目等其他魚類形成獨(dú)立的分支,而兩棲類、哺乳類和爬行類形成獨(dú)立的分支。

    2.5 RARmRNA的組織表達(dá)特性

    半滑舌鰨2種RAR基因在所有檢測組織中都有表達(dá)。RARα mRNA在腎臟和眼中表達(dá)量最高,與除胃的其他組織差異顯著,在胃中的表達(dá)也較高,在腦、脾臟、鰓、無眼側(cè)白皮膚、性腺、有眼側(cè)肌肉、腸、心臟、無眼側(cè)黑化皮膚、無眼側(cè)肌肉、有眼側(cè)皮膚等其他組織中也檢測到一定的表達(dá)量。半滑舌鰨RARγmRNA在脾臟和鰓中表達(dá)量最高,在心臟和腎臟表達(dá)量也較高,與其他組織差異顯著;而在無眼側(cè)白皮膚、腦、眼、胃、性腺、無眼側(cè)黑化皮膚、無眼側(cè)肌肉、垂體等其他組織表達(dá)量相對(duì)較低。在腎臟中,2種RAR mRNA表達(dá)水平都很高。脾臟、心臟、鰓、腎臟、無眼側(cè)白皮膚等組織中RARγmRNA表達(dá)量高于RARα,而在眼、胃、腸、性腺、肌肉、肝臟、有眼側(cè)皮膚、無眼側(cè)黑化皮膚等組織中,RARα mRNA表達(dá)量高于RARγ,脾臟、心臟、鰓和腎臟組織中RARγ與RARα表達(dá)差異顯著,而其他組織中RARγ與RARα表達(dá)差異不顯著。腦中2種RAR mRNA表達(dá)量基本一致,表明這2種RAR基因在不同組織中的生理功能可能存在差異,即使在同一組織中其生理功能也多存在一定的差異(圖6)。

    圖1 半滑舌鰨RARα基因cDNA全長序列及推導(dǎo)的氨基酸序列

    推導(dǎo)的氨基酸序列用單字母表示,從陰影顯示的起始甲硫氨酸開始計(jì)數(shù)。終止密碼子用*表示。下同

    The deduced amino acid residues were represented as single letter abbreviations and numbered from the initiating methionine which was shadowed. Termination codon was marked with *. The same as below

    圖2 半滑舌鰨RARγ基因cDNA全長序列及推導(dǎo)的氨基酸序列

    圖3 SWISS-MODEL預(yù)測的半滑舌鰨RARα(左)和RARγ(右)蛋白質(zhì)三級(jí)結(jié)構(gòu)

    分析有眼側(cè)皮膚、無眼側(cè)黑化皮膚和無眼側(cè)白皮膚中2種RAR基因mRNA的表達(dá)情況(圖7),發(fā)現(xiàn)無眼側(cè)白皮膚的2種RAR mRNA的表達(dá)量最高,其次為無眼側(cè)黑化皮膚,最低的為有眼側(cè)皮膚,其中,RARγ表達(dá)差異顯著。RARγ在無眼側(cè)未黑化皮膚中的表達(dá)顯著高于RARα (<0.05),而在正常有眼側(cè)皮膚和無眼側(cè)黑化皮膚中,RARα的表達(dá)高于RARγ,表明可能RARα和RARγ對(duì)皮膚組織中色素的調(diào)控作用具有不同的調(diào)控機(jī)制。

    3 討論

    本研究獲得了半滑舌鰨RAR的2個(gè)亞型RARα和RARγcDNA序列全長,并研究了其組織表達(dá)特性,為研究RA/RAR系統(tǒng)對(duì)體色異常調(diào)控提供了基礎(chǔ)材料。本研究獲得了半滑舌鰨2個(gè)RAR亞型的結(jié)構(gòu),預(yù)測未發(fā)現(xiàn)跨膜結(jié)構(gòu)和信號(hào)肽(Napoli, 1996)。同其他脊椎動(dòng)物一致,半滑舌鰨RAR具備A~F共6個(gè)功能域,在進(jìn)化過程中高度保守,其中,C區(qū)最為保守,具有鋅指結(jié)構(gòu)功能的為DBD區(qū)。E區(qū)在配體結(jié)合中起輔助作用,形成疏水氨基酸殘基。D區(qū)作為鉸鏈連接DBD和LBD (Leid, 1992)。同源性分析表明,半滑舌鰨RARα和RARγ的氨基酸序列與牙鲆的同源性最高,達(dá)97%,與其他魚類、兩棲類、爬行類、嚙齒類和人的氨基酸同源性也處于較高水平,表明其在進(jìn)化過程中保守性較強(qiáng)。同時(shí),本研究顯示,半滑舌鰨RARα和RARγ的氨基酸同源性達(dá)60.8%,但在進(jìn)化樹上屬于2個(gè)不同的進(jìn)化分支,表明在進(jìn)化過程中2種RAR受體基因出現(xiàn)了進(jìn)化差異。

    表2 半滑舌鰨RAR氨基酸序列與其他脊椎動(dòng)物的同源性比較

    Tab.2 Comparison of homology of the precursor peptide sequences of RAR gene between C. semilaevis and other vertebrates

    圖4 半滑舌鰨與其他物種的RAR氨基酸序列比較

    “*”表示一致的氨基酸;“:”表示高度保守度的氨基酸;“.”表示低保守度的氨基酸;陰影部分表示DBD和 LBD功能結(jié)構(gòu)域;RAR氨基酸序列號(hào)見表2;CS:半滑舌鰨;DR:斑馬魚LJ:鱸魚;HS:人

    Asterisks (*) indicated identical amino acid sequences; Dot (:) indicated highly conserved amino acid sequences; Dot (.) indicatedamino acid sequences of low degree conserved; GenBank accession numbers were shown in Tab.2. The shadow part represents the DBD and LBD functional domains CS:;DR:; LJ:;HS:

    圖5 基于RAR氨基酸序列的NJ系統(tǒng)進(jìn)化樹

    圖6 半滑舌鰨RAR mRNA在不同組織中的相對(duì)表達(dá)量

    BR:腦;EYE:眼;GI:鰓;H:心臟;L:肝臟;SP:脾臟;K:腎臟;ST:胃;I:腸;GO:性腺;EM:有眼側(cè)肌肉;BM:無眼側(cè)肌肉;ES:有眼側(cè)皮膚;BHS:無眼側(cè)黑化皮膚;BWS:無眼側(cè)白皮膚;P:垂體。不同組織表達(dá)差異分析中Y和y字母代表單獨(dú)分析RARγ的2個(gè)顯著差異集合,RARα單獨(dú)分析使用a、b、c,代表3個(gè)顯著差異集合,同一組織2種基因表達(dá)顯著差異(<0.05)用大括號(hào)表示,下同

    B: Brain; EYE: Eye; GI: Gill; H: Heart; L: Liver; SP: Spleen; K: Kindey; ST: Stomach; I: Intestine; GO: Gonad; EM: Eye-side muscle; BM: Blind-side muscle; ES: Eye-side skin; BHS: Blind-side hypermelanosis skin; BWS: Blind-side white skin; P: Pituitary. In different tissue expression analysis, Y and y letters represent two distinct different sets of RARγ analysis separately, RARα analysis alone, a, b and c represent three significant different sets. The significant difference in the expression of two genes in the same tissue is expressed in braces (<0.05), the same as below

    圖7 半滑舌鰨RAR mRNA在皮膚組織中的相對(duì)表達(dá)量

    ES:有眼側(cè)皮膚;BHS:無眼側(cè)黑化皮膚;BWS:無眼側(cè)白皮膚

    ES: Eye-side skin; BHS: Blind-side hypermelanosis skin; BWS: Blind-side white skin

    本研究發(fā)現(xiàn),半滑舌鰨2種RAR基因mRNA在所有檢測組織中均有表達(dá),表明其均具有廣泛的生理作用。RARα mRNA在腎臟中表達(dá)量最高,表明腎臟可能為其主要靶器官,而同時(shí)在眼、腦和胃中的表達(dá)量高于其他組織,表明RARα在這些組織器官的生理功能中都可能起重要的調(diào)控作用,這種組織表達(dá)分布特征與斑馬魚(Joore, 1994)、鱸魚(錢云霞等, 2012)以及哺乳動(dòng)物(Meng, 2011) RARα的組織表達(dá)特性相似。另外,半滑舌鰨RARγmRNA的主要靶器官為脾臟和鰓,同時(shí),在心臟和腎臟中均有高表達(dá)量,表明RARγ主要在這些器官中起重要的表達(dá)調(diào)控作用。2種RAR受體mRNA都在腎臟中具有高表達(dá),推測腎臟可能是RA/RAR信號(hào)通路的重要作用靶點(diǎn),具體作用方式有待于下一步深入研究。對(duì)虹鱒()的研究結(jié)果表明,鰓是RAR的主要表達(dá)器官,推測與其較高的新陳代謝速率有關(guān),也可能與RAR調(diào)節(jié)細(xì)胞生長與凋亡功能相關(guān)(Alsop, 2001)。本研究也發(fā)現(xiàn),鰓中2種RAR受體基因的表達(dá)量較高,但其具體的生理功能有待于進(jìn)一步研究。在對(duì)雞()的RARγ研究中,發(fā)現(xiàn)雞與哺乳動(dòng)物有相似的表達(dá)特異性,即RARγ表達(dá)高度限制在皮膚中(Michaille, 1994),然而,半滑舌鰨RAR表達(dá)在皮膚和肌肉相對(duì)其他組織較少,這種組織分布的差異主要可能是由種的差異引起的。

    本研究中,比較了有眼側(cè)皮膚、無眼側(cè)白皮膚和無眼側(cè)黑化皮膚中2種RAR基因mRNA的表達(dá)情況和兩者的相互關(guān)系,發(fā)現(xiàn)在正常有眼側(cè)皮膚中和在無眼側(cè)黑化皮膚中RARα的表達(dá)略高于RARγ,而在無眼側(cè)未發(fā)生黑化的皮膚中,RARγ的表達(dá)卻顯著高于RARα。先前對(duì)半滑舌鰨色素細(xì)胞的研究表明,黑色素與虹彩細(xì)胞等的數(shù)量分布以有眼側(cè)皮膚中最多,無眼側(cè)黑化皮膚中其次,而無眼側(cè)白皮膚中無黑色素細(xì)胞分布(史學(xué)營等, 2015)。這種不同狀態(tài)的皮膚組織中RAR轉(zhuǎn)錄產(chǎn)物表達(dá)與黑色素細(xì)胞分布的關(guān)系表明,RARγ與RARα均參與了半滑舌鰨皮膚組織中黑色素細(xì)胞的生長發(fā)育與分布調(diào)控過程,但其在皮膚組織中黑色素細(xì)胞的形成方面具有差異表達(dá)調(diào)控作用。與RARγ相比,RARα與皮膚中黑色素細(xì)胞的生長及數(shù)量分布調(diào)控可能具有更為密切的關(guān)系。今后應(yīng)結(jié)合鲆鰈類黑色素相關(guān)基因(如MCH、MCHR、POMC等)的表達(dá)和作用機(jī)理進(jìn)行深入研究,探討它們之間的相互作用關(guān)系。

    在牙鲆研究中,RARs與配體ATRA結(jié)合作用于苗種骨骼變態(tài)和體色沉著(Haga, 2003)。Shao等(2017)通過轉(zhuǎn)錄組研究發(fā)現(xiàn),牙鲆和半滑舌鰨有眼側(cè)皮膚ATRA和9-cis-RA濃度均高于無眼側(cè)皮膚,但兩側(cè)皮膚中RAR和RXR的基因表達(dá)沒有明顯差異。本研究未能測定不同類型皮膚中的RAR配體濃度,是否是配體濃度梯度的差異導(dǎo)致受體基因表達(dá)量的差異有待于今后深入研究。結(jié)合本研究RAR在皮膚中的分布,推測半滑舌鰨RAR可能與黑色素細(xì)胞的發(fā)育及分布具有不同的關(guān)系。已有研究證明,在斑馬魚中存在RARα-a和RARα-b亞型(Hale, 1993),而半滑舌鰨基因組預(yù)測其具有RARα、RARβ和RARγ三種亞型,同時(shí),每一種亞型又有數(shù)種不同的分子形式(Chen, 2014),開展RAR基因在不同組織器官中的表達(dá)特征和可能的生理功能研究有助于闡明RA/RAR系統(tǒng)在無眼側(cè)黑化發(fā)生的機(jī)理。本研究中,還發(fā)現(xiàn)半滑舌鰨RARα和RARγ mRNA在腎臟、脾臟、腦、眼、胃、無眼側(cè)白皮膚、性腺等組織中同時(shí)具有高表達(dá),且在不同的組織中2種受體基因表達(dá)量不同,表明半滑舌鰨2種RAR基因在不同的組織中生理功能不同。Shao等(2017)在成年牙鲆兩側(cè)皮膚中發(fā)現(xiàn),感光視蛋白對(duì)光照刺激有應(yīng)答反應(yīng),在光刺激下會(huì)影響配體RA在體兩側(cè)的濃度差,作用于RA/RAR系統(tǒng)導(dǎo)致無眼側(cè)發(fā)生黑化。工廠化養(yǎng)殖環(huán)境中發(fā)現(xiàn),幼魚通過補(bǔ)充RAR的配體RA控制了鲆鰈類有眼側(cè)無黑色素沉著現(xiàn)象(白化),而過量的RA補(bǔ)充卻會(huì)導(dǎo)致無眼側(cè)黑色素沉著過多(黑化)(Miwa, 1999)。研究發(fā)現(xiàn),改變光照(Shao, 2017)、鋪沙(Estevez, 2001)、棲息環(huán)境顏色(Takahashi, 2004)和養(yǎng)殖密度(Bolker, 2000)等也會(huì)影響無眼側(cè)黑色素沉著,且這些外源刺激都與攝食RA或通過眼和皮膚的感光視蛋白接收刺激產(chǎn)生RA,從而調(diào)控RAR表達(dá),并最終作用于鲆鰈類的體色沉著有一定關(guān)系。Isojima等(2014)研究表明,鲆鰈類體色發(fā)生與分布受到多種基因表達(dá)調(diào)控。深入開展RA/RAR系統(tǒng)對(duì)體色的調(diào)控作用及其信號(hào)通路的系統(tǒng)研究將有助于揭示半滑舌鰨無眼側(cè)黑化的分子調(diào)控機(jī)制。

    Alsop D, Brown S, Kraak GVD. Development of a retinoic acid receptor-binding assay with rainbow trout tissue: Characterization of retinoic acid binding, receptor tissue distribution, and developmental changes. General and Comparative Endocrinology, 2001, 123(3): 254–267

    Bolker JA, Hill CR. Pigmentation development in hatchery- reared ?at?shes. Journal of Fish Biology, 2000, 56(5): 1029–1052

    Brand N, Petkovich M, Krust A,. Identification of a 2nd human retinoic acid receptor. Nature, 1988, 332(6167): 850– 853

    Budhu AS, Noy N. Direct channeling of retinoic acid between cellular retinoic acid-binding protein II and retinoic acid receptor sensitizes mammary carcinoma cells to retinoic acid-induced growth arrest. Molecular & Cellular Biology, 2002, 22(8): 2632–2641

    Chambon P. A decade of molecular biology of retinoic acid receptors. FASEB Journal, 1996, 10(9): 940–954

    Chen SL, Zhang GJ, Shao CW,. Whole-genome sequence of a flatfish provides insights into ZW sex chromosome evolution and adaptation to a benthic lifestyle. Nature Genetics, 2014, 46(3): 253–260

    De Luca LM. Retinoids and their receptors in differentiation, embryogenesis, and neoplasia. FASEB Journal, 1991, 5(14): 2924–2933

    Deng JY, Meng TX, Ren SM,. Composition and quantity distribution of the Bohai Sea fish species. Marine Fisheries Research, 1988(9): 11–89 [鄧景耀, 孟田湘, 任勝民, 等. 渤海魚類種類組成及數(shù)量分布. 海洋水產(chǎn)研究, 1988(9): 11–89]

    Erkelens MN, Mebius RE. Retinoic acid and immune homeostasis: A balancing act. Trends in Immunology, 2017, 38(3): 168–180

    Estevez A, Kaneko T, Seikai T,. Ontogeny of ACTH and MSH cells in Japanese flounder () in relation to albinism. Aquaculture, 2001, 202(1): 131–143

    Faehnrich B, Lukas B, Humer E,. Phytogenic pigments in animal nutrition: Potentials and risks. Journal of the Science of Food and Agriculture, 2016, 96(5): 1420–1430

    Haga Y, Suzuki T, Takeuchi T,. A retinoic acid receptor- selective agonist causes jaw deformity in the Japanese flounder,. Aquaculture, 2003, 221(1): 381–392

    Haga Y, Suzuki T, Takeuchi T,. Retinoic acid isomers produce malformations in postembryonic development of the Japanese flounder,. Zoological Science, 2002, 19(10): 1105–1112

    Hale LA, Tallafuss A, Yan YL,. Characterization of the retinoic acid receptor genes RARaa, RARab and RARg during zebrafish development. Gene Expression Patterns, 2006, 6(5): 546–555

    Hallenbeck PL, Marks MS, Lippoldt RE,. Heterodimerization of thyroid hormone (TH) receptor with H-2RIIBP (RXR beta) enhances DNA binding and TH-dependent transcriptional activation. Proceedings of the National Academy of Science of the United States of America, 1992, 89(12): 5572–5576

    Herget T, Specht H, Esdar C,. Retinoic acid induces apoptosis-associated neural differentiation of a murine teratocarcinoma cell line. Journal of Neurochemistry, 1998, 70(1): 47–58

    Ignacio F, Matas P, Juanb OD,. Effect of dietary vitamin A on Senegalese sole () skeletogenesis and larval quality. Aquaculture, 2009, 295(3– 4): 250–265

    Isojima T, Tagawa M. Transplantation of pigmented and non-pigmented scales into the ocular and blind sides of the Japanese flounder, suggesting the presence of ocular-side characteristic inducer in pigmented scales. Fisheries Science, 2014, 80(5): 1027–1035

    Jones BB, Ohno CK, Allenby G,. New retinoid X receptor subtypes in zebra fish () differentially modulate transcription and do not bind 9-cis retinoic acid. Molecular and Cellular Biology, 1995, 15(10): 5226–5234

    Joore J, van der Lans GB, Lanser PH,. Effects of retinoic acid on the expression of retinoic acid receptors during zebrafish embryogenesis. Mechanisms of Development, 1994, 46(2): 137–150

    Kliewer SA, Umesono K, Noonnan DJ,. Convergence of 9-cis retinoic acid and peroxisome proliferator signalling pathways through heterodimer formation of their receptors. Nature, 1992, 358(6389): 771–774

    Leid M, Kastner P, Chambon P,. Multiplicity generates diversity in the retinoic acid signalling pathways. Trends in Biochemical Sciences, 1992, 17(10): 427–460

    Leng XJ. Research advance in vitamin nutrition of aquatic animals. Feed Industry, 2017, 38(16): 1–6 [冷向軍. 水產(chǎn)動(dòng)物維生素營養(yǎng)研究進(jìn)展. 飼料工業(yè), 2017, 38(16): 1–6]

    Lewis-McCrea LM, Lall SP. Effects of phosphorus and vitamin C deficiency, vitamin A toxicity, and lipid peroxidation on skeletal abnormalities in Atlantic halibut (). Journal of Applied Ichthyology, 2010, 26(2): 334–343

    Liu XZ, Zhuang ZM, Ma AJ,. Operative technologies for seedling rearing ofGünther. Marine Fisheries Research, 2006, 27(2): 17–24 [柳學(xué)周, 莊志猛, 馬愛軍, 等. 半滑舌鰨苗種生產(chǎn)技術(shù)的開發(fā)研究. 海洋水產(chǎn)研究, 2006, 27(2): 17–24]

    Liu XZ, Zhuang ZM. Reproduction biology and culture technology of half-smooth tongue sole. Beijing: China Agriculture Press, 2014, 1–10 [柳學(xué)周, 莊志猛. 半滑舌鰨繁育理論與養(yǎng)殖技術(shù). 北京: 中國農(nóng)業(yè)出版社, 2014, 1–10]

    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-DDCTMethod. Methods, 2001, 25(4): 402–408

    Martinez GM, Baron MP, Bolker JA. Skeletal and pigmentation defects following retinoic acid exposure in larval summer flounder,. Journal of the World Aquaculture Society, 2007, 38(3): 353–366

    Meng QY, Chen XN, Zhao J,. Distribution of retinoic acid receptor-α immunoreactivity in the human hypothalamus. Neuroscience, 2011, 174: 132–142

    Michaille JJ, Blanchet S, Kanzler B,. Characterization of cDNAs encoding the chick retinoic acid receptor gamma 2 and preferential distribution of retinoic acid receptor gamma transcripts during chick skin development. Developmental Dynamics, 1994, 201(4): 334–343

    Miwa S, Yamano K. Retinoic acid stimulates development of adult-type chormatophores in the flounder. Journal of Experimental Zoology Part A, 1999, 284(3): 317–324

    Napoli JL. Biochemical pathways of retinoid transport, metabolism, and signal transduction. Clinical Immunology and Immunopathology, 1996, 80(3): 52–62

    Qian YX, Han L, Tong CH,. Retinoic acid receptor α and retinoic acid receptor γ in: cDNA clonging and gene expression analysis. Acta Nutrimenta Sinica, 2012, 24(4): 712–721 [錢云霞, 韓柳, 童彩環(huán), 等. 鱸魚視黃酸受體α和視黃酸受體γ cDNA克隆和基因表達(dá)分析. 動(dòng)物營養(yǎng)學(xué)報(bào), 2012, 24(4): 712–721]

    Shao CW, Chen SL, Bao BL,. The genome and transcriptome of Japanese flounder provide insights into flatfish asymmetry. Nature Genetics, 2017, 49(1): 119–124

    Shi XY, Liu XZ, Shi Y,. Molecular characterization ofand its corelation with blind-side hypermelanosis inProgress in Fishery Sciences, 2017. 38(1): 91–102 [史學(xué)營, 柳學(xué)周, 石瑩, 等. 半滑舌鰨()黑色素聚集素受體()表達(dá)特性及其與無眼側(cè)黑化的關(guān)系. 漁業(yè)科學(xué)進(jìn)展, 2017, 38(1): 91–102]

    Shi XY, Liu XZ, Xu YJ,. Preliminary studies on blind-side hypermelanosis of: Chromatophores observation and expression of proopiomelanocortin. Progress in Fishery Sciences, 2015, 36(2): 45–54 [史學(xué)營, 柳學(xué)周, 徐永江, 等. 半滑舌鰨()體表色素細(xì)胞觀察及POMC表達(dá)特性分析. 漁業(yè)科學(xué)進(jìn)展, 2015, 36(2): 45–54]

    Takahashi A, Tsuchiya K, Yamanome T,. Possible involvement of melanin concentrating hormone in food intake in a teleost fish, barfin flounder. Peptides, 2004, 25(10): 1613–1622

    Xu YJ, Zhu XW, Liu XZ,. Cloning and expression of melanin-concentrating hormone in half-smooth tongue sole (). Progress in Fishery Sciences, 2017, 38(1): 81–90 [徐永江, 朱學(xué)武, 柳學(xué)周, 等.半滑舌鰨()黑色素富集激素基因的克隆和表達(dá). 漁業(yè)科學(xué)進(jìn)展, 2017, 38(1): 81–90]

    Zhang JL, Zuo ZH, Zhu WW,. Sex-different effects of tributyltin on brain aromatase, estrogen receptor and retinoid X receptor gene expression in rockfish (). Marine Environmental Research, 2013, 90(3): 113–118

    Zhou WD, Wang LL, Zhou LB,. All-trans retinoic acid upregulates the expression of ciliary neurotrophic factor in retinal pigment epithelial cells. Cell Biochemistry and Function, 2017, 35(4): 202–208

    Zhu XW, Xu YJ, Liu XY,. In vitro expression and bioactivity analysis of melanin concentration hormone from. Journal of Fisheries of China, 2016, 40(10): 1595–1605 [朱學(xué)武, 徐永江, 柳學(xué)周, 等. 半滑舌鰨黑色素聚集激素重組制備與生物活性分析. 水產(chǎn)學(xué)報(bào), 2016, 40(10): 1595–1605]

    Molecular Cloning and Spatial Expression of Two Retinoic Acid Receptors RARalpha and RARgammafrom

    SONG Xuesong1,2, XU Yongjiang1, LIU Xuezhou1①, SHI Bao1, WANG Bin1, LIU Yongshan1,2, ZHANG Yaxing1,2

    (1. Laboratory for Marine Fisheries Science and Food Production Processes, Pilot National Laboratory for Marine Science and Technology (Qingdao), Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071; 2. College of Fisheries and Life Science, Shanghai Ocean University, Shanghai 201306)

    Two retinoic acid receptors, RARalpha and RARgamma, were cloned fromusing RT-PCR and RACE methods, and their spatial expression patterns were investigated using a quantitative PCR assay. The results showed thatfull-length cDNA sequence encoding the RARalpha gene is 1823 bp in length, its open reading frame (ORF) length is 1332 bp, encoding 443 amino acids; the RARgamma cDNA sequence is 1959 bp in length, and the length of ORF is 1497 bp, encoding 489 amino acids. Homology analysis showed thatRARalpha and RARgamma have homology identity of 60.8%, and both have 97% homology identity with the Japanese flounder. Phylogenetic analysis showed thatRARalpha and RARgamma clustered into a separate branch with other fish counterparts. Spatial expression analysis showed that the highest expression level of RARalpha mRNA occurred in the kidney, whereas the highest expression level of RARgamma mRNA was in the spleen. RARgamma was also highly expressed in the gill, kidney, and heart. Furthermore, RARalpha and RARgamma mRNA expression were detected in all examined tissues, which indicated that these two retinoic acid receptors were both involved in multiple physiological regulation processes. In addition, the differential expression of these two RAR genes were found in the blind side and ocular skins, wherein they both had highest expression levels in the normal blind side skin, followed by the blacking blind side skin and had the lowest expression levels in the ocular side skin. In ocular skin and blind-side blacking skin, the RARalpha expression levels were higher than RARgamma but without significant difference, wherein in blind-side normal skin, the RARgamma expressed significantly higher than RARalpha. This differential expression pattern indicated that they might play important physiological roles in blind-side hypermelanosis regulation of.

    ; Retinoic acid receptor; Gene cloning; Expression pattern; Hyperpigmentation

    LIU Xuezhou, E-mail: liuxz@ysfri.ac.cn

    宋雪松, 徐永江, 柳學(xué)周, 史寶, 王濱, 劉永山, 張雅星. 半滑舌鰨() 2種視黃酸受體RARα和RARγ克隆及組織表達(dá)特性. 漁業(yè)科學(xué)進(jìn)展, 2018, 39(6): 52–64

    Song XS, Xu YJ, Liu XZ, Shi B, Wang B, Liu YS, Zhang YX. Molecular cloning and spatial expression of two retinoic acid receptors RARalpha and RARgammafrom. Progress in Fishery Sciences, 2018, 39(6): 52–64

    * 中國水產(chǎn)科學(xué)研究院中央級(jí)公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金(2017GH05; 2017GH17)、國家海水魚產(chǎn)業(yè)技術(shù)體系(CARS-47)、國家自然科學(xué)基金(31502145; 31602133)、中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所基本科研業(yè)務(wù)費(fèi)項(xiàng)目(20603022017016)共同資助 [This work was supported by Central Public-Interest Scientific Institution Basal Research Fund, Chinese Academy of Fishery Sciences (CAFS) (2017GH05; 2017GH17), China Agricultural Research System (CARS-47), National Natural Science Foundation of China (31502145; 31602133), and Central Public-interest Scientific Institution Basal Research Fund, YSFRI, CAFS (20603022017016)]. 宋雪松,E-mail: 746284973@qq.com

    柳學(xué)周,研究員,E-mail: liuxz@ysfri.ac.cn

    2017-12-20,

    2018-01-29

    10.19663/j.issn2095-9869.20171220003

    TS201.4;S917.4

    A

    2095-9869(2018)06-0052-13

    (編輯 馮小花)

    猜你喜歡
    黑化舌鰨牙鲆
    高黏高彈改性瀝青在城市橋面黑化中的應(yīng)用
    溶解氧對(duì)半滑舌鰨呼吸頻率、生理生化指標(biāo)和組織結(jié)構(gòu)的影響
    短吻三線舌鰨、褐斑三線舌鰨和紫斑舌鰨種類的有效性厘定
    北戴河站在牙鲆幼魚培育光色選擇研究方面取得新進(jìn)展
    “黑化”的白云
    我們不是被網(wǎng)絡(luò)黑化的小怪獸
    雌核發(fā)育、野生及養(yǎng)殖牙鲆形態(tài)差異分析
    黑化處理對(duì)2A14鋁合金激光-TIG復(fù)合焊焊縫成形的影響
    飼料維生素C含量對(duì)半滑舌鰨抵抗遲緩愛德華氏菌感染的影響
    雄性半滑舌鰨GH基因多態(tài)性與激素水平及生長性狀相關(guān)性分析
    久久久国产成人免费| 一本一本综合久久| 国产精品伦人一区二区| 搞女人的毛片| 欧美黑人巨大hd| 天堂√8在线中文| 亚洲午夜理论影院| 成年女人毛片免费观看观看9| 婷婷色综合大香蕉| 天堂影院成人在线观看| 亚洲精华国产精华精| 91狼人影院| 在线免费观看不下载黄p国产 | 狂野欧美白嫩少妇大欣赏| 我要看日韩黄色一级片| 3wmmmm亚洲av在线观看| 欧美另类亚洲清纯唯美| 午夜福利在线观看吧| 免费大片18禁| 国产伦精品一区二区三区视频9| 一区二区三区激情视频| 搡老熟女国产l中国老女人| 夜夜爽天天搞| 中文字幕熟女人妻在线| 久久国产乱子免费精品| 成人特级黄色片久久久久久久| 免费黄网站久久成人精品 | 免费人成在线观看视频色| 亚洲人成伊人成综合网2020| 亚洲aⅴ乱码一区二区在线播放| 国产私拍福利视频在线观看| 亚洲电影在线观看av| 免费看光身美女| 免费观看精品视频网站| 欧美色欧美亚洲另类二区| 十八禁人妻一区二区| 国产成人影院久久av| 床上黄色一级片| 国产一区二区在线观看日韩| 国产精品一区二区性色av| 成人永久免费在线观看视频| 成熟少妇高潮喷水视频| 成年人黄色毛片网站| 国产黄色小视频在线观看| 五月玫瑰六月丁香| 在线观看一区二区三区| 欧美+日韩+精品| 99久久成人亚洲精品观看| 91麻豆av在线| 亚洲国产色片| 亚洲aⅴ乱码一区二区在线播放| 国产欧美日韩精品亚洲av| 无遮挡黄片免费观看| 动漫黄色视频在线观看| 国产成+人综合+亚洲专区| av在线老鸭窝| 久久久国产成人精品二区| 国产人妻一区二区三区在| 亚洲国产色片| 女人十人毛片免费观看3o分钟| 国产主播在线观看一区二区| 国产精品综合久久久久久久免费| 黄色配什么色好看| 国产欧美日韩精品一区二区| 午夜福利成人在线免费观看| 麻豆国产av国片精品| 69av精品久久久久久| 亚洲熟妇熟女久久| 久久精品国产亚洲av香蕉五月| 日韩欧美国产一区二区入口| 精品不卡国产一区二区三区| 国产精品99久久久久久久久| 有码 亚洲区| 99久久久亚洲精品蜜臀av| 欧美不卡视频在线免费观看| 韩国av一区二区三区四区| 色av中文字幕| 人妻久久中文字幕网| 欧美色欧美亚洲另类二区| 亚洲成a人片在线一区二区| 免费观看的影片在线观看| 俄罗斯特黄特色一大片| 国产av不卡久久| 久久精品国产清高在天天线| 精品人妻视频免费看| 久久人人精品亚洲av| 亚洲 国产 在线| 成年版毛片免费区| 日韩精品中文字幕看吧| 成人精品一区二区免费| 在线播放国产精品三级| 特大巨黑吊av在线直播| 国内精品久久久久久久电影| www.www免费av| 一区二区三区四区激情视频 | 日韩欧美国产一区二区入口| 在线天堂最新版资源| 在线十欧美十亚洲十日本专区| 国产精品人妻久久久久久| 国产大屁股一区二区在线视频| 特大巨黑吊av在线直播| 中出人妻视频一区二区| 麻豆国产av国片精品| 男女视频在线观看网站免费| 舔av片在线| 久久人人爽人人爽人人片va | 午夜亚洲福利在线播放| 久久久久久国产a免费观看| 美女黄网站色视频| 国产高清激情床上av| 国产欧美日韩一区二区三| 国产成年人精品一区二区| 无遮挡黄片免费观看| 亚洲精品亚洲一区二区| 每晚都被弄得嗷嗷叫到高潮| 亚洲中文字幕日韩| 国产精品98久久久久久宅男小说| 国产精品99久久久久久久久| 午夜日韩欧美国产| 精品久久久久久久人妻蜜臀av| 亚洲人成网站高清观看| 日韩欧美国产在线观看| 午夜福利成人在线免费观看| 99久国产av精品| 亚洲成a人片在线一区二区| 麻豆av噜噜一区二区三区| 成人性生交大片免费视频hd| 99精品久久久久人妻精品| 欧美色视频一区免费| 欧美色欧美亚洲另类二区| 久久性视频一级片| а√天堂www在线а√下载| 亚洲黑人精品在线| 欧美成人一区二区免费高清观看| 琪琪午夜伦伦电影理论片6080| 丰满人妻一区二区三区视频av| 国产不卡一卡二| 国产三级在线视频| 日韩人妻高清精品专区| 在线观看av片永久免费下载| 精品一区二区三区视频在线观看免费| 免费无遮挡裸体视频| 黄色一级大片看看| 免费看光身美女| 少妇的逼水好多| 国产日本99.免费观看| 美女xxoo啪啪120秒动态图 | 成人鲁丝片一二三区免费| 亚洲一区高清亚洲精品| aaaaa片日本免费| 美女 人体艺术 gogo| 中文字幕精品亚洲无线码一区| 99精品久久久久人妻精品| 精品熟女少妇八av免费久了| 午夜精品久久久久久毛片777| 国产欧美日韩一区二区精品| 亚洲五月婷婷丁香| 桃红色精品国产亚洲av| 午夜免费激情av| 欧美乱妇无乱码| 久久性视频一级片| 最近在线观看免费完整版| 麻豆成人午夜福利视频| 欧美日韩福利视频一区二区| 亚洲狠狠婷婷综合久久图片| 国产精品免费一区二区三区在线| 少妇丰满av| 精品国产三级普通话版| 黄色女人牲交| 哪里可以看免费的av片| 亚洲精品色激情综合| 婷婷色综合大香蕉| 3wmmmm亚洲av在线观看| 午夜福利欧美成人| 国产91精品成人一区二区三区| 日韩精品青青久久久久久| 少妇人妻一区二区三区视频| 日本 av在线| 国产91精品成人一区二区三区| 宅男免费午夜| 九色成人免费人妻av| 亚洲国产欧洲综合997久久,| 国产成人福利小说| 88av欧美| 久久久久久国产a免费观看| 日本免费a在线| 欧美zozozo另类| 久久天躁狠狠躁夜夜2o2o| 成人无遮挡网站| 成年人黄色毛片网站| 一级作爱视频免费观看| 一本久久中文字幕| 禁无遮挡网站| 欧美xxxx性猛交bbbb| 国产视频内射| 亚洲第一欧美日韩一区二区三区| 成人亚洲精品av一区二区| av黄色大香蕉| 国产激情偷乱视频一区二区| 国产大屁股一区二区在线视频| 淫妇啪啪啪对白视频| 赤兔流量卡办理| 少妇高潮的动态图| 国产精品久久久久久久电影| 老熟妇仑乱视频hdxx| 自拍偷自拍亚洲精品老妇| 亚洲激情在线av| 欧美性猛交╳xxx乱大交人| 精品久久久久久久久av| 国产精品精品国产色婷婷| 国产亚洲精品久久久com| 麻豆av噜噜一区二区三区| 一级黄片播放器| 成年女人毛片免费观看观看9| 九色国产91popny在线| 亚洲av成人精品一区久久| 少妇的逼水好多| 老熟妇乱子伦视频在线观看| 亚洲精品在线美女| 自拍偷自拍亚洲精品老妇| 国产91精品成人一区二区三区| 久久久精品欧美日韩精品| 午夜亚洲福利在线播放| 欧美xxxx黑人xx丫x性爽| 国产精品嫩草影院av在线观看 | 男人舔女人下体高潮全视频| 男女床上黄色一级片免费看| 亚洲精品粉嫩美女一区| 偷拍熟女少妇极品色| 国产亚洲精品久久久com| 精品日产1卡2卡| 亚洲美女搞黄在线观看 | 久久这里只有精品中国| 99在线人妻在线中文字幕| 男女之事视频高清在线观看| 亚洲av免费高清在线观看| 两人在一起打扑克的视频| 一卡2卡三卡四卡精品乱码亚洲| 18美女黄网站色大片免费观看| 男人的好看免费观看在线视频| 最近最新免费中文字幕在线| 亚洲人成伊人成综合网2020| 欧美高清性xxxxhd video| 99精品在免费线老司机午夜| 国产午夜精品论理片| 亚洲av二区三区四区| 婷婷精品国产亚洲av在线| 婷婷六月久久综合丁香| 久久久久久国产a免费观看| 免费搜索国产男女视频| 午夜老司机福利剧场| 亚洲成人久久爱视频| 午夜福利成人在线免费观看| 窝窝影院91人妻| 国产毛片a区久久久久| 三级国产精品欧美在线观看| 2021天堂中文幕一二区在线观| 国产私拍福利视频在线观看| 男女之事视频高清在线观看| a在线观看视频网站| 中文字幕精品亚洲无线码一区| 日韩欧美在线二视频| 真人做人爱边吃奶动态| 老熟妇乱子伦视频在线观看| 91久久精品电影网| 男女下面进入的视频免费午夜| 久久精品久久久久久噜噜老黄 | 精品一区二区三区av网在线观看| 日韩欧美在线乱码| 国产视频一区二区在线看| 亚洲精品日韩av片在线观看| 麻豆国产97在线/欧美| 欧美黑人巨大hd| 亚洲精品久久国产高清桃花| 国产单亲对白刺激| 直男gayav资源| 日日摸夜夜添夜夜添av毛片 | 又黄又爽又刺激的免费视频.| 少妇的逼水好多| 一个人看的www免费观看视频| 亚洲 欧美 日韩 在线 免费| 亚洲精品成人久久久久久| 国产亚洲欧美在线一区二区| 岛国在线免费视频观看| 在现免费观看毛片| 老司机福利观看| 小蜜桃在线观看免费完整版高清| 麻豆成人av在线观看| 亚洲真实伦在线观看| 波多野结衣高清作品| 99久久无色码亚洲精品果冻| 国产中年淑女户外野战色| 久久国产精品人妻蜜桃| 99国产精品一区二区三区| 午夜福利高清视频| 国产精品一及| 51午夜福利影视在线观看| 很黄的视频免费| 女人被狂操c到高潮| 成人欧美大片| 亚洲人成网站在线播放欧美日韩| 日本 av在线| 免费人成视频x8x8入口观看| 别揉我奶头~嗯~啊~动态视频| 99热这里只有精品一区| 成人永久免费在线观看视频| 三级国产精品欧美在线观看| 久久精品国产亚洲av香蕉五月| 一级毛片久久久久久久久女| bbb黄色大片| 中文在线观看免费www的网站| 国产精品一区二区三区四区免费观看 | 欧美+日韩+精品| 成年版毛片免费区| 日韩欧美一区二区三区在线观看| 一区二区三区四区激情视频 | 亚洲人成电影免费在线| 乱码一卡2卡4卡精品| 黄色配什么色好看| 精品久久久久久久久亚洲 | 亚洲,欧美精品.| 啪啪无遮挡十八禁网站| 午夜久久久久精精品| 国产午夜福利久久久久久| 俄罗斯特黄特色一大片| 最近最新免费中文字幕在线| 女人被狂操c到高潮| 国产伦在线观看视频一区| 97人妻精品一区二区三区麻豆| 国产亚洲欧美在线一区二区| 久久性视频一级片| 高潮久久久久久久久久久不卡| 波野结衣二区三区在线| 国产一区二区三区视频了| 久久久久久久久大av| 一夜夜www| 国产69精品久久久久777片| 久久久久免费精品人妻一区二区| 97碰自拍视频| 国内毛片毛片毛片毛片毛片| 色视频www国产| 在线观看66精品国产| 久久精品综合一区二区三区| 亚洲精品久久国产高清桃花| 亚洲第一欧美日韩一区二区三区| 亚洲国产日韩欧美精品在线观看| 日本免费一区二区三区高清不卡| 国产高清视频在线观看网站| 十八禁人妻一区二区| 久久伊人香网站| 在线观看免费视频日本深夜| 他把我摸到了高潮在线观看| 亚洲黑人精品在线| 美女高潮的动态| 韩国av一区二区三区四区| 最好的美女福利视频网| 免费搜索国产男女视频| 直男gayav资源| 日韩中文字幕欧美一区二区| 宅男免费午夜| 青草久久国产| av专区在线播放| 亚洲av成人不卡在线观看播放网| 琪琪午夜伦伦电影理论片6080| 国产欧美日韩一区二区精品| 国内精品久久久久精免费| 国内精品一区二区在线观看| 在线观看一区二区三区| 亚洲五月天丁香| 乱码一卡2卡4卡精品| .国产精品久久| 国产伦精品一区二区三区视频9| 欧美成人a在线观看| 国产精品久久电影中文字幕| 婷婷精品国产亚洲av| 亚洲av.av天堂| 青草久久国产| 国产欧美日韩精品一区二区| 国产高清视频在线播放一区| 亚洲av免费在线观看| 十八禁人妻一区二区| av在线老鸭窝| 身体一侧抽搐| 男女视频在线观看网站免费| 日韩有码中文字幕| 国产在线男女| 午夜福利在线在线| 日本免费a在线| 亚洲男人的天堂狠狠| 三级男女做爰猛烈吃奶摸视频| 我的女老师完整版在线观看| 亚洲最大成人手机在线| 99精品在免费线老司机午夜| 国产av在哪里看| 嫩草影院新地址| 亚洲不卡免费看| 久久香蕉精品热| 国产欧美日韩精品亚洲av| 成人av在线播放网站| 日日摸夜夜添夜夜添小说| 国产三级中文精品| 欧美一区二区亚洲| 内地一区二区视频在线| 亚洲国产高清在线一区二区三| 一二三四社区在线视频社区8| 国产中年淑女户外野战色| 伊人久久精品亚洲午夜| 丰满的人妻完整版| 给我免费播放毛片高清在线观看| 看片在线看免费视频| 亚洲av成人不卡在线观看播放网| 免费大片18禁| 亚洲国产精品久久男人天堂| 不卡一级毛片| 亚洲av电影在线进入| 久久精品国产99精品国产亚洲性色| 欧美绝顶高潮抽搐喷水| 国产免费av片在线观看野外av| 男人的好看免费观看在线视频| 精品国内亚洲2022精品成人| 午夜免费成人在线视频| а√天堂www在线а√下载| 直男gayav资源| 国产免费一级a男人的天堂| 成年女人永久免费观看视频| 桃色一区二区三区在线观看| 成人国产综合亚洲| 国产成年人精品一区二区| 国产成人aa在线观看| 国产伦一二天堂av在线观看| 日韩欧美一区二区三区在线观看| 99riav亚洲国产免费| 一级黄片播放器| 欧美高清性xxxxhd video| 国产极品精品免费视频能看的| 天美传媒精品一区二区| 午夜视频国产福利| 一a级毛片在线观看| 精品无人区乱码1区二区| 久久热精品热| 久久午夜亚洲精品久久| 国产老妇女一区| 久久亚洲真实| 黄色一级大片看看| 欧美一区二区亚洲| 真人做人爱边吃奶动态| 黄片小视频在线播放| 一级黄片播放器| 精品一区二区三区视频在线观看免费| 免费观看人在逋| 国产在线精品亚洲第一网站| 亚洲国产精品sss在线观看| 一边摸一边抽搐一进一小说| 看片在线看免费视频| 久久亚洲精品不卡| a在线观看视频网站| 好看av亚洲va欧美ⅴa在| 宅男免费午夜| 国产精品嫩草影院av在线观看 | 69人妻影院| av在线老鸭窝| 日韩成人在线观看一区二区三区| 国模一区二区三区四区视频| 日日摸夜夜添夜夜添小说| 一a级毛片在线观看| 国产人妻一区二区三区在| 亚洲va日本ⅴa欧美va伊人久久| 五月玫瑰六月丁香| 一级av片app| 三级国产精品欧美在线观看| 亚洲七黄色美女视频| 欧美潮喷喷水| 亚洲乱码一区二区免费版| 十八禁人妻一区二区| 日日摸夜夜添夜夜添av毛片 | 97人妻精品一区二区三区麻豆| 亚洲五月婷婷丁香| 国产伦一二天堂av在线观看| 免费电影在线观看免费观看| 国产日本99.免费观看| 午夜两性在线视频| 国产精品女同一区二区软件 | 中文字幕免费在线视频6| 国产精品综合久久久久久久免费| 欧美成人性av电影在线观看| 国产探花在线观看一区二区| 精品久久久久久久久亚洲 | 给我免费播放毛片高清在线观看| 亚洲,欧美,日韩| 亚洲精品色激情综合| 精品99又大又爽又粗少妇毛片 | 女同久久另类99精品国产91| 免费大片18禁| 人妻丰满熟妇av一区二区三区| 亚洲精品成人久久久久久| 精品人妻熟女av久视频| 在线观看免费视频日本深夜| 日韩欧美免费精品| 99热精品在线国产| 九九久久精品国产亚洲av麻豆| 久久久久性生活片| 999久久久精品免费观看国产| 欧美在线一区亚洲| 亚洲国产日韩欧美精品在线观看| 美女高潮喷水抽搐中文字幕| 中文字幕免费在线视频6| 国产av在哪里看| 国产真实乱freesex| 亚洲人成网站在线播放欧美日韩| 免费观看人在逋| 国产av麻豆久久久久久久| 亚洲激情在线av| 男人狂女人下面高潮的视频| 国产爱豆传媒在线观看| 国产成人影院久久av| 少妇熟女aⅴ在线视频| 久久精品综合一区二区三区| 国产在视频线在精品| 18禁黄网站禁片免费观看直播| 国产三级中文精品| 欧美乱妇无乱码| 国产乱人伦免费视频| 免费看美女性在线毛片视频| 免费在线观看日本一区| 国产精品一及| 人妻久久中文字幕网| 此物有八面人人有两片| 亚洲最大成人手机在线| 黄色日韩在线| 欧美性猛交╳xxx乱大交人| 欧美一级a爱片免费观看看| 一区福利在线观看| 国产成人欧美在线观看| 亚洲 国产 在线| 亚洲欧美清纯卡通| 黄色一级大片看看| 免费大片18禁| 少妇的逼水好多| 亚洲熟妇中文字幕五十中出| 天堂影院成人在线观看| 3wmmmm亚洲av在线观看| 国产亚洲欧美98| 亚洲精品亚洲一区二区| 欧美色视频一区免费| 中文亚洲av片在线观看爽| 亚洲欧美日韩卡通动漫| 国产伦一二天堂av在线观看| 精品久久国产蜜桃| 亚洲在线观看片| 国产私拍福利视频在线观看| 亚洲人成网站在线播| 亚洲国产精品999在线| 无人区码免费观看不卡| 人妻夜夜爽99麻豆av| 午夜福利视频1000在线观看| av福利片在线观看| 免费搜索国产男女视频| 伊人久久精品亚洲午夜| 美女xxoo啪啪120秒动态图 | 国产高清激情床上av| 91麻豆av在线| 最好的美女福利视频网| 三级毛片av免费| 成人毛片a级毛片在线播放| 欧美成狂野欧美在线观看| 九九在线视频观看精品| 久久天躁狠狠躁夜夜2o2o| 在线十欧美十亚洲十日本专区| 黄色视频,在线免费观看| 波多野结衣巨乳人妻| 国产国拍精品亚洲av在线观看| 国产成人a区在线观看| 免费电影在线观看免费观看| 乱码一卡2卡4卡精品| 香蕉av资源在线| 级片在线观看| 国产欧美日韩一区二区三| 别揉我奶头 嗯啊视频| 97人妻精品一区二区三区麻豆| 麻豆av噜噜一区二区三区| 精品人妻1区二区| 在线看三级毛片| 国产野战对白在线观看| 亚洲成a人片在线一区二区| 亚洲 国产 在线| 国产成人aa在线观看| 午夜福利高清视频| 国产精品美女特级片免费视频播放器| 国产高潮美女av| 日韩有码中文字幕| 日本免费一区二区三区高清不卡| 简卡轻食公司| 久久久久久九九精品二区国产| 午夜精品久久久久久毛片777| 我的老师免费观看完整版| 真人一进一出gif抽搐免费| 国内精品美女久久久久久| 午夜免费激情av| 永久网站在线| 能在线免费观看的黄片| 午夜免费成人在线视频| 精品99又大又爽又粗少妇毛片 | 12—13女人毛片做爰片一| 成人av在线播放网站| 又粗又爽又猛毛片免费看| 国产成人av教育| 精品人妻1区二区| 日本黄色片子视频| 日韩成人在线观看一区二区三区| 成人鲁丝片一二三区免费| 亚洲乱码一区二区免费版| 国产成人欧美在线观看| 婷婷色综合大香蕉| 日韩欧美在线二视频| 他把我摸到了高潮在线观看|