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

    烏鱧和斑鱧EST序列微衛(wèi)星信息分析

    2016-06-14 09:42:42馮曉宇姚桂桂潘彬斌
    淡水漁業(yè) 2016年3期
    關(guān)鍵詞:烏鱧微衛(wèi)星堿基

    謝 楠,劉 凱,馮曉宇,姚桂桂,潘彬斌

    (杭州市農(nóng)業(yè)科學(xué)研究院,杭州 310024)

    ?

    烏鱧和斑鱧EST序列微衛(wèi)星信息分析

    謝楠,劉凱,馮曉宇,姚桂桂,潘彬斌

    (杭州市農(nóng)業(yè)科學(xué)研究院,杭州310024)

    摘要:烏鱧(Channa argus,♂)、斑鱧(Channa maculatus,♀)是雜交鱧“杭鱧1號(hào)”的親本。利用高通量測(cè)序技術(shù)對(duì)烏鱧和斑鱧的肝臟進(jìn)行轉(zhuǎn)錄組測(cè)序獲得大量EST序列后,利用MISA軟件進(jìn)行微衛(wèi)星信息分析,結(jié)果表明,通過(guò)轉(zhuǎn)錄組測(cè)序獲得烏鱧EST序列59 959條,長(zhǎng)度45 mb,發(fā)現(xiàn)15 428個(gè)SSR,出現(xiàn)頻率為25.73%;獲得斑鱧EST序列44 337條,長(zhǎng)度27 mb,發(fā)現(xiàn)8 439個(gè)SSR,出現(xiàn)頻率為19.03%。在烏鱧和斑鱧EST-SSR中,重復(fù)單元以1~2堿基重復(fù)為最多,并以長(zhǎng)度小于16 bp的短重復(fù)序列為主,間隔SSR和復(fù)合SSR的EST序列RPKM均值要低于單純型SSR的EST序列RPKM均值,且在單純型SSR中SSR長(zhǎng)度越長(zhǎng),其RPKM均值則越低。

    關(guān)鍵詞:烏鱧(Channa argus); 斑鱧(Channa maculatus);EST序列;微衛(wèi)星

    烏鱧(Channaargus)和斑鱧(Channamaculatus)是雜交鱧“杭鱧1號(hào)”(以下簡(jiǎn)稱“雜交鱧”)的親本,同屬于鱸形目鱧科鱧屬。通過(guò)烏鱧(♂)、斑鱧(♀)雜交獲得的F1代新品種雜交鱧“杭鱧1號(hào)”可采用人工飼料投喂模式養(yǎng)殖,兼有生長(zhǎng)速度快、養(yǎng)殖周期短、產(chǎn)量高、病害少、效益高等優(yōu)點(diǎn),有效地解決了傳統(tǒng)烏鱧投喂冰鮮魚(yú)模式造成的面源污染及產(chǎn)品食用安全問(wèn)題[1],并成為浙江省養(yǎng)殖漁業(yè)主推品種之一,替代傳統(tǒng)烏鱧的比例不斷增加,成為杭州地區(qū)的黑魚(yú)養(yǎng)殖主導(dǎo)品種。

    目前,針對(duì)雜交鱧的研究在養(yǎng)殖密度、營(yíng)養(yǎng)成分、染色體核型以及雜交子代與親本的遺傳標(biāo)記等方面已有相關(guān)報(bào)道[2-4],但相關(guān)分子標(biāo)記的開(kāi)發(fā)量仍然較少,還不能完全滿足研究的需求。微衛(wèi)星( SSR) 作為第二代分子標(biāo)記,被越來(lái)越多地用于水生動(dòng)物群體遺傳學(xué)分析工作。本研究基于高通量測(cè)序技術(shù)對(duì)雜交鱧親本烏鱧和斑鱧進(jìn)行轉(zhuǎn)錄組測(cè)序,獲得相關(guān)EST序列并進(jìn)行SSR信息分析,以期能為雜交鱧的標(biāo)記開(kāi)發(fā)、分子標(biāo)記輔助育種等方面提供基礎(chǔ)。

    1材料與方法

    1.1實(shí)驗(yàn)材料

    實(shí)驗(yàn)用烏鱧為采自錢(qián)塘江群體的后代,經(jīng)篩選獲得的體質(zhì)健康的個(gè)體(體重1 000 g左右);實(shí)驗(yàn)用斑鱧為采自珠江流域的親本經(jīng)過(guò)人工繁殖培育的后代,經(jīng)篩選獲得的體質(zhì)健康的個(gè)體(體重750 g左右)。實(shí)驗(yàn)用樣品組織為肝臟,烏鱧和斑鱧各采集6尾,將采集到的6個(gè)肝臟組織混合后分別作為轉(zhuǎn)錄組建庫(kù)及測(cè)序樣本,先經(jīng)液氮瞬時(shí)冷凍,之后于-80 ℃低溫保存?zhèn)溆谩?/p>

    1.2EST序列的獲得

    從烏鱧和斑鱧的肝臟提取mRNA構(gòu)建cDNA文庫(kù),建好的文庫(kù)用Illumina Solexa進(jìn)行測(cè)序。烏鱧和斑鱧測(cè)序數(shù)據(jù)量分別達(dá)到8.4 G和5.7 G,Q30數(shù)據(jù)比率分別為88.13%和86.18%,Q20數(shù)據(jù)比率分別為94.55%和92.90%。經(jīng)過(guò)質(zhì)量處理后,烏鱧和斑鱧轉(zhuǎn)錄本留存率分別為97.89%和96.53%,可以滿足后續(xù)數(shù)據(jù)分析要求。序列拼接及去重復(fù)后,得到轉(zhuǎn)錄本。烏鱧和斑鱧樣本中分別得到85 802個(gè)和52 769個(gè)長(zhǎng)度大于200 bp的轉(zhuǎn)錄本,大小分別為113.75 mb和44.29 mb,截取每個(gè)基因座下最長(zhǎng)的轉(zhuǎn)錄本后分別得到59 959個(gè)和44 337個(gè)EST序列,大小分別為45 mb和27 mb,序列平均長(zhǎng)度為751.99 bp和608.88 bp。

    1.3基因表達(dá)豐度

    對(duì)拼接得到的轉(zhuǎn)錄本(≥200 bp)建庫(kù)后,使用軟件bwa(http://sourceforge.net/projects/bio-bwa)中single-end的mapping方法計(jì)算轉(zhuǎn)錄本在各樣本中的表達(dá)豐度,允許一個(gè)reads比對(duì)到多個(gè)轉(zhuǎn)錄本上,此時(shí)每個(gè)reads在每個(gè)轉(zhuǎn)錄本上計(jì)1/n,每個(gè)EST序列下面的所有轉(zhuǎn)錄本比對(duì)上的reads均計(jì)于該EST序列下。豐度值采用RPKM(Reads Per Kilobase of exon model per Million mapped reads)法計(jì)算。

    1.4SSR分析

    以Illumina Solexa測(cè)序獲得的EST序列作為參考序列,使用SSR篩選軟件MISA(http://pgrc.ipk-gatersleben.de/misa/misa.html)進(jìn)行 SSR篩選。篩選標(biāo)準(zhǔn)為1個(gè)堿基重復(fù)10次及10次以上,2個(gè)堿基重復(fù)6次及6 次以上,3~6個(gè)堿基重復(fù)5次及5次以上,兩個(gè)微衛(wèi)星之間的距離小于100 bp的時(shí)候,兩個(gè)微衛(wèi)星組成一個(gè)復(fù)合微衛(wèi)星。用SSR出現(xiàn)頻率和SSR平均分布距離來(lái)描述SSR,出現(xiàn)頻率前2位的重復(fù)單元定義為優(yōu)勢(shì)重復(fù)單元。計(jì)算公式分別為:SSR出現(xiàn)頻率=搜索到的SSR數(shù)量/EST序列數(shù)量;SSR平均分布頻率=EST序列總堿基數(shù)/搜索到的SSR數(shù)量。利用Primer3 interface modules ( http://pgrc.ipk-gatersleben.de/misa/primer3.html )對(duì)SSR篩選結(jié)果進(jìn)行預(yù)處理后,利用primer3 ( http://sourceforge.net/projects/primer3 )進(jìn)行SSR引物的批量設(shè)計(jì),引物設(shè)計(jì)的參數(shù)是Tm為60 ℃,引物長(zhǎng)度為20 bp。

    2結(jié)果

    2.1烏鱧和斑鱧EST序列中SSR重復(fù)單元分布

    利用MISA軟件對(duì)烏鱧和斑鱧的EST序列進(jìn)行SSR特征分析,其中檢測(cè)的烏鱧EST序列共59959條,長(zhǎng)度為45 mb,共發(fā)現(xiàn)15 428個(gè)SSR,出現(xiàn)頻率為25.73%,平均每2.92 kb出現(xiàn)1個(gè)SSR;檢測(cè)的斑鱧EST序列共44 337條,長(zhǎng)度為27 mb,共發(fā)現(xiàn)8 439個(gè)SSR,出現(xiàn)頻率為19.03%,平均每3.20 kb出現(xiàn)1個(gè)SSR。檢出的SSR 包含了1、2、3,4和5堿基重復(fù)單元,不同重復(fù)單元出現(xiàn)頻率隨堿基數(shù)增加而明顯降低,其中最常見(jiàn)的是1堿基和2堿基重復(fù)單元,各重復(fù)單元SSR出現(xiàn)頻率及分布頻率信息見(jiàn)表1。

    表1 EST序列中不同重復(fù)單元SSR分布情況

    在搜索到的烏鱧EST序列的SSR 中,共觀察到45種不同的重復(fù)單元類(lèi)型,其中1堿基重復(fù)2 種,2 堿基重復(fù)4種,3堿基重復(fù)10種,4堿基重復(fù)23種,5 堿基重復(fù)6種。1堿基重復(fù)中,A/T類(lèi)型出現(xiàn)的頻率最高,占94.47%;2堿基重復(fù)中,AC/GT出現(xiàn)的頻率最高,占78.61%;3堿基重復(fù)中,AGC/CTG出現(xiàn)的頻率最高,占24.94%;4堿基重復(fù)中,AAAC/GTTT出現(xiàn)的頻率最高,占25.29%;5堿基重復(fù)中,AAAAT/ATTTT出現(xiàn)的頻率最高,占28.57%,結(jié)果見(jiàn)表2。

    表2 烏鱧EST序列中不同重復(fù)單元SSR的出現(xiàn)頻率(考慮到堿基互補(bǔ)作用)

    在搜索到的斑鱧EST序列SSR 中,共觀察到40種不同的重復(fù)單元類(lèi)型,其中1堿基重復(fù)2 種,2 堿基重復(fù)4種,3堿基重復(fù)10種,4堿基重復(fù)17種,5 堿基重復(fù)7種。1堿基重復(fù)中,A/T類(lèi)型出現(xiàn)的頻率最高,占95.14%;2堿基重復(fù)中,AC/GT出現(xiàn)的頻率最高,占75.56%;3堿基重復(fù)中,AGC/CTG出現(xiàn)的頻率最高,占27.45%;4堿基重復(fù)中,AAAC/GTTT出現(xiàn)的頻率最高,占33.71%,結(jié)果見(jiàn)表3。

    表3 斑鱧EST序列中不同重復(fù)單元SSR的出現(xiàn)頻率(考慮到堿基互補(bǔ)作用)

    2.2烏鱧和斑鱧EST序列中SSR的長(zhǎng)度分布

    對(duì)烏鱧和斑鱧的EST序列中所發(fā)現(xiàn)的SSR長(zhǎng)度的變化進(jìn)行了分析,結(jié)果發(fā)現(xiàn)序列所含SSR在長(zhǎng)度上存在顯著變異,SSR長(zhǎng)度從10到237個(gè)堿基不等,烏鱧SSR的加權(quán)平均長(zhǎng)度為17.11 bp,斑鱧SSR的加權(quán)平均長(zhǎng)度為17.20 bp。烏鱧和斑鱧SSR以長(zhǎng)度小于16 bp的短重復(fù)序列最多,其次是長(zhǎng)度為16~25 bp的序列重復(fù),大于30 bp的較長(zhǎng)序列重復(fù)僅為4.0%左右,結(jié)果如圖1,利用SPSS進(jìn)行Pearson相關(guān)性分析表明,烏鱧和斑鱧微衛(wèi)星數(shù)量和長(zhǎng)度呈顯著負(fù)相關(guān)性,相關(guān)系數(shù)分別為-0.391和-0.447,P=0.000<0.01,具有顯著相關(guān)性。

    以SSR數(shù)量對(duì)堿基重復(fù)數(shù)量作圖,結(jié)果如圖2,可見(jiàn)SSR數(shù)量隨著堿基重復(fù)次數(shù)增加呈明顯下降趨勢(shì),而且堿基重復(fù)次數(shù)越少,SSR數(shù)量下降的速率就越快。僅當(dāng)2堿基重復(fù)從8次增加到10次時(shí),烏鱧和斑鱧SSR數(shù)量出現(xiàn)了增加的情況。1堿基重復(fù)次數(shù)達(dá)到22次、2堿基重復(fù)次數(shù)達(dá)到12次、其他堿基重復(fù)次數(shù)達(dá)到8次的時(shí)候,SSR數(shù)量的下降速率降低,進(jìn)入平臺(tái)期。

    圖2 EST序列中SSR數(shù)量隨重復(fù)次數(shù)變化曲線

    烏鱧59 959條和斑鱧44 337條EST序列中,含SSR的EST序列分別占總數(shù)的19.02%和15.08%。將含SSR的EST序列與不含SSR的EST序列的RPKM均值比較,烏鱧中含SSR的EST序列的RPKM均值為16.85,是不含SSR的EST序列的RPKM均值5.89的2.86倍,斑鱧中含SSR的EST序列的RPKM均值為42.48,是不含SSR的EST序列的RPKM均值6.18的6.87倍。Weber[5]將微衛(wèi)星分為3類(lèi),即單純(pure) SSR、復(fù)合(compound) SSR和間隔(interrupted) SSR。本實(shí)驗(yàn)發(fā)現(xiàn),烏鱧中含有間隔SSR和復(fù)合SSR的EST序列的RPKM均值分別為6.53、3.44,低于僅含有單純型SSR的EST序列的RPKM均值17.43;斑鱧中含有間隔SSR和復(fù)合SSR的EST序列的RPKM均值分別為20.56、7.07,低于僅含有單純型SSR的EST序列的RPKM均值43.60,結(jié)果見(jiàn)表4。對(duì)SSR長(zhǎng)度與相關(guān)EST序列的RPKM值進(jìn)行相關(guān)性分析,SSR長(zhǎng)度與EST序列的RPKM值相關(guān)系數(shù)為-0.124,無(wú)顯著相關(guān)性。

    表4 EST序列中SSR類(lèi)型與長(zhǎng)度對(duì)RPKM值的影響

    以烏鱧和斑鱧的EST序列為基礎(chǔ),利用primer3軟件進(jìn)行批量引物設(shè)計(jì)。烏鱧中14 561個(gè)含SSR序列,7 297條序列引物設(shè)計(jì)失敗,成功設(shè)計(jì)7 264對(duì)SSR引物,成功率為49.89%。斑鱧中8 017個(gè)含SSR序列,4 607條序列引物設(shè)計(jì)失敗,成功設(shè)計(jì)3 410對(duì)SSR引物,成功率為42.53%。

    3討論

    由重復(fù)單元的重復(fù)次數(shù)不同所造成的多態(tài)性常常表現(xiàn)為復(fù)等位性,在不同的基因型間存在廣泛的多態(tài)性。研究者認(rèn)為這種多態(tài)性起因于復(fù)制過(guò)程中的滑動(dòng)[16]。因此,SSR長(zhǎng)度的變化反映SSR位點(diǎn)獲得(或失去)重復(fù)單元的活躍程度。越來(lái)越多的研究表明,SSR參與基因的表達(dá)調(diào)控及基因重排、變異等活動(dòng)并與腫瘤發(fā)生及神經(jīng)系統(tǒng)疾病有關(guān)[17-20]。在羅非魚(yú)(Oreochromismossambicus)研究中發(fā)現(xiàn),SSR長(zhǎng)度與基因表達(dá)強(qiáng)度相關(guān),在無(wú)環(huán)境脅迫的條件下SSR長(zhǎng)度越長(zhǎng),基因表達(dá)水平越低,在有環(huán)境脅迫的條件下則反之[21]。本實(shí)驗(yàn)發(fā)現(xiàn),烏鱧和斑鱧測(cè)序所獲得的EST序列中,間隔SSR和復(fù)合SSR的EST序列的RPKM均值要低于單純型SSR的EST序列的RPKM均值,且在單純型SSR中SSR長(zhǎng)度越長(zhǎng),其RPKM均值則越低,這與羅非魚(yú)上的研究結(jié)論類(lèi)似。而在楊樹(shù)和茶樹(shù)EST-SSR的研究中認(rèn)為,含有SSR的基因總體表達(dá)水平偏低,這樣有利于相應(yīng)基因在基因組中的保存,規(guī)避部分選擇壓力[22-23]。但本實(shí)驗(yàn)中發(fā)現(xiàn),不論烏鱧和斑鱧,含SSR的EST序列的RPKM均值要明顯高于不含SSR的EST序列的RPKM均值,與楊樹(shù)和茶樹(shù)的研究結(jié)果不同,這也許由于物種差異造成或?qū)嶒?yàn)方法上的差別造成,有待進(jìn)一步研究。此外,本實(shí)驗(yàn)中使用primer3對(duì)轉(zhuǎn)錄本進(jìn)行批量SSR引物設(shè)計(jì),大大提高了工作效率,為下一步工作的開(kāi)展奠定了基礎(chǔ)。

    參考文獻(xiàn):

    [1]馮曉宇,王宇希,李行先,等.池塘養(yǎng)殖雜交鱧高效試驗(yàn)[J].科學(xué)養(yǎng)魚(yú),2008,(3):31-32.

    [2]劉新軼,王宇希,謝楠,等.雜交鱧及其親本肌肉營(yíng)養(yǎng)成分分析[J].齊魯漁業(yè),2010,27(5):24-26.

    [3]謝楠,馮曉宇,劉新軼,等.4個(gè)雜交鱧繁殖性能、F1生長(zhǎng)和形態(tài)的比較[J].浙江農(nóng)業(yè)科學(xué),2013,(10):1366-1369.

    [4]朱樹(shù)人,李忠全,馮曉宇,等.烏鱧和斑鱧微衛(wèi)星制備及其對(duì)自交與雜交F1代的初步鑒定[J].動(dòng)物學(xué)雜志,2011,46 (4):1-7.

    [5]Weber J L Informativeness of human (dC-dA)n(dG-dT)npoly-morphisms[J] Genomics,1990,7(4):524-530

    [6]Zhan A B,Bao Z M,Hu X L.Characterization of 95 novel microsatellite markers for Zhikong scallopChlamysfareriusing FIASCO colony hybridization and EST database mining[J].Fish Sci,2008,74(3):516-526.

    [7]Wang H X,Li F H,Xiang J H.Polymorphic EST-SSR markers and their mode of inheritance inFenneropenaeuschinensis[J].Aquaculture,2005,249(s1-4):107-114.

    [8]Wang D,Liao X L,Cheng L, et al.Development of novel EST-SSR markers in common carp by data mining from public EST sequences[J].Aquaculture,2007,271:558-574.

    [9]許曉軍,張海琪,張超,等.中華鱉表達(dá)序列標(biāo)簽資源中的微衛(wèi)星信息分析[J].經(jīng)濟(jì)動(dòng)物學(xué)報(bào),2013,17(1):15-18.

    [10]Serapion J,Kucuktas H,Feng J N,etal.Bioinformatic mining of type I microsatellites from expressed sequence tags of channel catfish(Ictaluruspunctatus)[J].Mar Biotechnol,2004,6(4):364-377.

    [11]陳松波,龔麗,劉海金.牙鲆EST資源的SSR信息分析[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào),2010,41(10):82-86.

    [12]王艷紅,胡超群,張呂平,等.凡納濱對(duì)蝦EST微衛(wèi)星標(biāo)記初步篩選[J].大連水產(chǎn)學(xué)院學(xué)報(bào),2008,23(4):304-308.

    [13]劉博,邵艷卿,滕爽爽,等.縊蟶(Sinonovaculaconstricta)EST-SSR分布特征及引物開(kāi)發(fā)利用[J].海洋與湖沼,2012,43(1):133-137.

    [14]Zhang A,Wang Y,Brown B,et al.Isolation and characterization of novel microsatellite markers for yellow perch(Percaflavescens)[J].Int Mol Sci,2009,10(1):18-27.

    [15]張琳琳,李莉,張國(guó)范.長(zhǎng)牡蠣(Crassostreagigas)EST串聯(lián)重復(fù)序列的組成和分布[J].海洋科學(xué),2011,35(4):9-14.

    [16]Sharopova N.Plant simple sequence repeats:distribution,variation,and effects on gene expression[J] Genome,2008,51:79-90

    [17]Thibodeau S N,Bren G,Schaid D.Microsatellite instability in cancer of the proximal colon[J].Science,1993,260(5109):816-819.

    [18]Aaltonen L A,Lauri A,Leach F S,et al.Clues to the pathogenesis of familial colorctal cancer[J].Science,1993,260(5109):812-816.

    [19]Ionov Y,Miguel A,Peinado,et al.Ubiquitous somatic mutations in simple repeated sequences reveal a new mechanism for colon carcinogenesis[J].Nature,1993,363:558-561.

    [20]Toth G,Gaspari Z,Jurka J.Microsatellites in different eukaryotic genomes:survey and analysis[J].Genome Research,2000,10(7):967-981.

    [21]Streelman J,Kocher T.Microsatellite variation associated with prolactin expression and growth of salt-challenged Tilapia[J].Physiol Genom,2002,9(1):1-4.

    [22]劉菁菁,戴曉港,王潔,等.楊樹(shù)微衛(wèi)星序列對(duì)基因表達(dá)頻率的影響及表達(dá)序列中微衛(wèi)星特征的分析[J].南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,35(1):11-14.

    [23]王麗鴛,韋康,張成才,等.茶樹(shù)花轉(zhuǎn)錄組微衛(wèi)星分布特征[J].作物學(xué)報(bào),2014,40(1):80-85.

    (責(zé)任編輯:陳細(xì)華)

    Analysis of simple sequence repeats information in expressed sequence tags of Channa argus and Channa maculatus

    XIE Nan,LIU Kai,FENG Xiao-yu,YAO Gui-gui,PAN Bin-bin

    (HangzhouAcademyofAgriculturalSciences,Hangzhou310024,China)

    Abstract:In this experiment,a large number of Expressed Sequence Tags(EST) were obtained from transcriptome sequences of livers in Channa argus and Channa maculatus by high-throughput sequencing technology,and Simple Sequence Repeats(SSR) Information was analysed.The result showed that,59 959 EST with 45 mb length were obtained from Channa argus,and 15428 EST-SSR were mined out by MISA,with the frequency of 25.73%;while 44 337 EST with 27mb length were obtained from Channa maculatus,and 8 439 EST-SSR were mined out by MISA,with the frequency of 19.03%.In the EST-SSR of C.argus and C.maculatus,nucleotide and dinucleotide repeats were the main types,while the repeats less than 16bp in length were the main types.Mean values of reads per kilobase per million (RPKM) of EST in both interrupted and compound SSR were lower than mean values of RPKM of EST in pure SSR,and in the pure EST-SSR,the longer SSR in length,the lower mean values of RPKM.

    Key words:Channa argus;Channa maculatus;expressed sequence tags;simple sequence repeats

    收稿日期:2015-08-11;

    修訂日期:2015-10-12

    第一作者簡(jiǎn)介:謝楠(1981-),男,碩士,主要從事魚(yú)類(lèi)遺傳育種研究工作。E-mail:n.xie@163.com通訊作者:馮曉宇。E-mail:fengxy1962@163.com

    中圖分類(lèi)號(hào):S917.4

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1000-6907-(2016)03-0016-07

    資助項(xiàng)目:浙江省水產(chǎn)農(nóng)業(yè)新品種選育重大科技專(zhuān)項(xiàng)(2012C12907-7)

    猜你喜歡
    烏鱧微衛(wèi)星堿基
    烏鱧
    垂釣(2024年2期)2024-04-29 00:00:00
    梅花鹿基因組微衛(wèi)星分布特征研究
    基于微衛(wèi)星標(biāo)記的不同地區(qū)烏鱧和白化烏鱧群體遺傳多樣性研究
    應(yīng)用思維進(jìn)階構(gòu)建模型 例談培養(yǎng)學(xué)生創(chuàng)造性思維
    中國(guó)科學(xué)家創(chuàng)建出新型糖基化酶堿基編輯器
    生命“字母表”迎來(lái)4名新成員
    烏鱧全人工飼料健康養(yǎng)殖技術(shù)
    生命“字母表”迎來(lái)4名新成員
    烏鱧與白烏鱧肌肉營(yíng)養(yǎng)成分分析與評(píng)價(jià)
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    夜夜爽夜夜爽视频| 亚洲成色77777| 精品不卡国产一区二区三区| 午夜福利网站1000一区二区三区| 哪个播放器可以免费观看大片| 国产一区二区亚洲精品在线观看| 久99久视频精品免费| 国产久久久一区二区三区| 欧美+日韩+精品| 欧美xxⅹ黑人| 99久久中文字幕三级久久日本| 午夜亚洲福利在线播放| 一级av片app| 人人妻人人看人人澡| 亚洲成人av在线免费| 亚洲精品久久久久久婷婷小说| 国产精品一区二区三区四区久久| 日本欧美国产在线视频| 99热这里只有是精品50| 熟女电影av网| 亚洲国产av新网站| 国产精品久久久久久av不卡| 别揉我奶头 嗯啊视频| 少妇熟女aⅴ在线视频| 免费看a级黄色片| 国产精品日韩av在线免费观看| 亚洲最大成人手机在线| 日韩视频在线欧美| 国产精品一区二区性色av| 99热网站在线观看| 精品久久久精品久久久| 真实男女啪啪啪动态图| 身体一侧抽搐| 五月玫瑰六月丁香| 国产精品久久久久久精品电影小说 | 高清日韩中文字幕在线| 久久久久九九精品影院| 亚洲人成网站高清观看| 永久免费av网站大全| 日韩大片免费观看网站| 最后的刺客免费高清国语| 成年av动漫网址| 亚洲国产日韩欧美精品在线观看| 夫妻性生交免费视频一级片| 免费看a级黄色片| 亚洲av男天堂| 国产男人的电影天堂91| 国产精品日韩av在线免费观看| 搡老乐熟女国产| 偷拍熟女少妇极品色| 国产精品女同一区二区软件| 久久久精品欧美日韩精品| 免费观看性生交大片5| 一区二区三区免费毛片| 欧美成人a在线观看| 欧美日韩国产mv在线观看视频 | 一级黄片播放器| 国产精品精品国产色婷婷| 国产高清国产精品国产三级 | 国模一区二区三区四区视频| 18禁在线无遮挡免费观看视频| 久久韩国三级中文字幕| 禁无遮挡网站| 亚洲精品成人av观看孕妇| 欧美丝袜亚洲另类| 午夜久久久久精精品| 老司机影院毛片| 亚洲怡红院男人天堂| 一个人观看的视频www高清免费观看| 久久久久久久国产电影| 插逼视频在线观看| 成人特级av手机在线观看| 成人av在线播放网站| 亚洲av在线观看美女高潮| 国产亚洲精品av在线| 2022亚洲国产成人精品| 综合色av麻豆| 如何舔出高潮| 国产精品一区二区在线观看99 | 成年免费大片在线观看| 免费无遮挡裸体视频| 欧美日韩精品成人综合77777| 亚洲aⅴ乱码一区二区在线播放| 成人亚洲精品av一区二区| 人妻制服诱惑在线中文字幕| 高清日韩中文字幕在线| 五月天丁香电影| 黄色欧美视频在线观看| 亚洲欧美精品专区久久| 国产精品久久久久久av不卡| 国产探花在线观看一区二区| 亚洲真实伦在线观看| 噜噜噜噜噜久久久久久91| 一级黄片播放器| 伦精品一区二区三区| 精品人妻视频免费看| 一个人免费在线观看电影| 一级av片app| 黄色一级大片看看| 亚洲自偷自拍三级| 亚洲精品成人av观看孕妇| 人妻一区二区av| 久久国产乱子免费精品| av卡一久久| 久久久久久久久久久丰满| 色网站视频免费| 美女黄网站色视频| 99久久精品国产国产毛片| 在线 av 中文字幕| 国内揄拍国产精品人妻在线| 少妇丰满av| 亚洲美女视频黄频| 丰满乱子伦码专区| 婷婷色综合www| 蜜桃亚洲精品一区二区三区| 99热这里只有是精品50| 免费大片黄手机在线观看| 中文天堂在线官网| 亚洲av成人精品一二三区| 色综合站精品国产| 一区二区三区免费毛片| 免费看日本二区| 国产日韩欧美在线精品| 十八禁网站网址无遮挡 | 亚洲成人中文字幕在线播放| 亚洲av成人av| 好男人视频免费观看在线| 国产黄色视频一区二区在线观看| 亚洲美女搞黄在线观看| 国产亚洲精品av在线| 身体一侧抽搐| 精品一区二区免费观看| 欧美一区二区亚洲| 成人av在线播放网站| 国产 一区 欧美 日韩| 少妇熟女aⅴ在线视频| 免费高清在线观看视频在线观看| av播播在线观看一区| 亚洲av中文av极速乱| 国产成年人精品一区二区| 69人妻影院| 国产美女午夜福利| 免费少妇av软件| 亚洲经典国产精华液单| 亚洲乱码一区二区免费版| 国产成人精品婷婷| 能在线免费看毛片的网站| 免费看av在线观看网站| 亚洲av成人精品一区久久| 欧美最新免费一区二区三区| 亚洲欧美日韩无卡精品| 欧美另类一区| 国产精品熟女久久久久浪| 国产乱人偷精品视频| 岛国毛片在线播放| eeuss影院久久| 国产国拍精品亚洲av在线观看| 国产成人a∨麻豆精品| 亚洲成人久久爱视频| 亚洲av成人精品一区久久| 国产精品1区2区在线观看.| 黄片无遮挡物在线观看| 夫妻午夜视频| 男插女下体视频免费在线播放| 又爽又黄无遮挡网站| 啦啦啦韩国在线观看视频| 欧美日韩一区二区视频在线观看视频在线 | 精品一区二区三区视频在线| 最近2019中文字幕mv第一页| 欧美丝袜亚洲另类| 欧美高清性xxxxhd video| 看黄色毛片网站| 九九在线视频观看精品| 日韩av不卡免费在线播放| 日韩大片免费观看网站| 欧美xxxx黑人xx丫x性爽| 日韩一区二区三区影片| 国产 一区 欧美 日韩| 一个人看的www免费观看视频| 97人妻精品一区二区三区麻豆| 久久久精品欧美日韩精品| 精品人妻熟女av久视频| 免费黄频网站在线观看国产| 国产成人精品福利久久| 国产黄频视频在线观看| 黑人高潮一二区| 好男人视频免费观看在线| 亚洲精品国产av蜜桃| 在线播放无遮挡| 免费无遮挡裸体视频| 亚洲精品456在线播放app| 99久久精品一区二区三区| 少妇熟女aⅴ在线视频| 亚洲精品,欧美精品| 亚洲精品国产成人久久av| 午夜爱爱视频在线播放| 日韩视频在线欧美| 亚洲av日韩在线播放| 国产一区亚洲一区在线观看| 精品久久久精品久久久| 色尼玛亚洲综合影院| 国产精品99久久久久久久久| 久久精品国产鲁丝片午夜精品| 亚洲欧美日韩东京热| 亚洲熟女精品中文字幕| 有码 亚洲区| 成人av在线播放网站| 欧美97在线视频| 伦理电影大哥的女人| 少妇人妻一区二区三区视频| 青春草亚洲视频在线观看| 51国产日韩欧美| 日韩欧美一区视频在线观看 | 久久99精品国语久久久| 精品一区在线观看国产| 少妇人妻精品综合一区二区| 成人一区二区视频在线观看| 亚洲电影在线观看av| 国产一区二区在线观看日韩| 视频中文字幕在线观看| 男插女下体视频免费在线播放| 伊人久久国产一区二区| 夫妻性生交免费视频一级片| 亚洲伊人久久精品综合| 丝袜美腿在线中文| 亚洲国产成人一精品久久久| 特级一级黄色大片| 亚洲乱码一区二区免费版| 国产精品女同一区二区软件| 99热这里只有是精品在线观看| 少妇的逼水好多| 日韩欧美 国产精品| 国产精品人妻久久久影院| 97热精品久久久久久| 99热这里只有精品一区| 国产伦理片在线播放av一区| 一区二区三区高清视频在线| 色综合亚洲欧美另类图片| av卡一久久| 直男gayav资源| 免费av不卡在线播放| 日日撸夜夜添| 91久久精品国产一区二区成人| 国产视频内射| 成人综合一区亚洲| 精品不卡国产一区二区三区| 欧美日韩精品成人综合77777| 亚洲av成人av| 久久久久精品性色| 亚洲丝袜综合中文字幕| 插逼视频在线观看| 婷婷六月久久综合丁香| 一夜夜www| 麻豆精品久久久久久蜜桃| 欧美潮喷喷水| 日韩,欧美,国产一区二区三区| 亚洲自拍偷在线| 我要看日韩黄色一级片| 国产精品国产三级国产专区5o| 国产成人a区在线观看| 欧美人与善性xxx| 高清午夜精品一区二区三区| 精品久久久久久久久亚洲| 舔av片在线| 99热6这里只有精品| 18+在线观看网站| 国产黄色免费在线视频| 又粗又硬又长又爽又黄的视频| 午夜老司机福利剧场| 日日啪夜夜爽| 你懂的网址亚洲精品在线观看| 国产男人的电影天堂91| 男人爽女人下面视频在线观看| 成人一区二区视频在线观看| 日日啪夜夜撸| 激情 狠狠 欧美| 国产亚洲一区二区精品| 国产一区二区三区综合在线观看 | 国产成人免费观看mmmm| 国产免费视频播放在线视频 | 国产精品一区二区三区四区免费观看| 五月伊人婷婷丁香| 亚洲国产精品成人综合色| 亚洲图色成人| 国产午夜福利久久久久久| 久久久久久国产a免费观看| 少妇人妻精品综合一区二区| 97超碰精品成人国产| 美女主播在线视频| 嫩草影院入口| 亚洲天堂国产精品一区在线| 欧美精品一区二区大全| 国产亚洲最大av| 麻豆精品久久久久久蜜桃| 大香蕉久久网| www.色视频.com| 午夜福利视频1000在线观看| 秋霞伦理黄片| 国产极品天堂在线| 搡老妇女老女人老熟妇| 亚洲图色成人| 人妻夜夜爽99麻豆av| 床上黄色一级片| 久久人人爽人人片av| 国产伦在线观看视频一区| 国产不卡一卡二| 一本久久精品| 亚洲在线观看片| 久久热精品热| 一个人免费在线观看电影| 精品99又大又爽又粗少妇毛片| 亚洲成人av在线免费| 大香蕉久久网| 97人妻精品一区二区三区麻豆| 少妇猛男粗大的猛烈进出视频 | 日韩电影二区| 免费高清在线观看视频在线观看| 久久6这里有精品| 午夜久久久久精精品| 男的添女的下面高潮视频| a级毛片免费高清观看在线播放| 成人鲁丝片一二三区免费| 国产精品.久久久| 国产精品一区二区三区四区久久| 韩国高清视频一区二区三区| 国产成人freesex在线| 亚洲伊人久久精品综合| 男人爽女人下面视频在线观看| 小蜜桃在线观看免费完整版高清| 久久99精品国语久久久| 丝袜美腿在线中文| 国产免费又黄又爽又色| 亚洲精品,欧美精品| 男女视频在线观看网站免费| 久久久a久久爽久久v久久| ponron亚洲| 观看美女的网站| 熟妇人妻不卡中文字幕| 身体一侧抽搐| 午夜亚洲福利在线播放| 久久久a久久爽久久v久久| 婷婷色麻豆天堂久久| 毛片一级片免费看久久久久| 免费人成在线观看视频色| av播播在线观看一区| 18禁裸乳无遮挡免费网站照片| 欧美潮喷喷水| 日韩伦理黄色片| 国产片特级美女逼逼视频| 熟妇人妻久久中文字幕3abv| 99久久中文字幕三级久久日本| 国产中年淑女户外野战色| 欧美不卡视频在线免费观看| av又黄又爽大尺度在线免费看| 少妇熟女欧美另类| 免费无遮挡裸体视频| 精品久久久久久久久av| 中文字幕制服av| 我的老师免费观看完整版| 国国产精品蜜臀av免费| 熟妇人妻不卡中文字幕| 亚洲国产av新网站| 久久久久九九精品影院| 午夜免费男女啪啪视频观看| 菩萨蛮人人尽说江南好唐韦庄| 最近的中文字幕免费完整| 国产成人精品福利久久| 精品久久国产蜜桃| 午夜精品国产一区二区电影 | 日韩视频在线欧美| 少妇人妻一区二区三区视频| 国产精品嫩草影院av在线观看| 天天一区二区日本电影三级| 黑人高潮一二区| 成人二区视频| 人妻夜夜爽99麻豆av| 国产av国产精品国产| 日韩欧美一区视频在线观看 | 亚洲婷婷狠狠爱综合网| kizo精华| 亚洲自偷自拍三级| 国产精品蜜桃在线观看| 两个人的视频大全免费| 91在线精品国自产拍蜜月| 联通29元200g的流量卡| 亚洲欧洲国产日韩| 国产精品国产三级国产专区5o| 国产有黄有色有爽视频| 亚洲精品久久午夜乱码| 国产成人精品一,二区| 免费高清在线观看视频在线观看| 一级毛片我不卡| 内射极品少妇av片p| 日本免费在线观看一区| 久久久久久久久久黄片| 少妇熟女欧美另类| 精品不卡国产一区二区三区| 视频中文字幕在线观看| 一级毛片aaaaaa免费看小| 亚洲成人中文字幕在线播放| 中文乱码字字幕精品一区二区三区 | 国产午夜精品论理片| 欧美zozozo另类| 成人欧美大片| 日本wwww免费看| 水蜜桃什么品种好| 久久久久国产网址| 男女边吃奶边做爰视频| 别揉我奶头 嗯啊视频| 国产精品国产三级专区第一集| 一区二区三区免费毛片| 国产午夜精品一二区理论片| 永久免费av网站大全| 全区人妻精品视频| 国产精品一区二区三区四区久久| 大香蕉97超碰在线| 51国产日韩欧美| 天天一区二区日本电影三级| 高清av免费在线| 日韩一区二区三区影片| 亚洲第一区二区三区不卡| 深爱激情五月婷婷| 免费黄频网站在线观看国产| av在线观看视频网站免费| 一本一本综合久久| 成人高潮视频无遮挡免费网站| 五月玫瑰六月丁香| 一级av片app| av.在线天堂| 国精品久久久久久国模美| 欧美一级a爱片免费观看看| 国产午夜精品久久久久久一区二区三区| 狂野欧美激情性xxxx在线观看| 又爽又黄a免费视频| 国产黄频视频在线观看| 欧美日韩国产mv在线观看视频 | 国产黄a三级三级三级人| 午夜精品国产一区二区电影 | 熟女人妻精品中文字幕| 三级经典国产精品| 久久99热这里只频精品6学生| www.av在线官网国产| 国产精品久久视频播放| 日韩欧美 国产精品| 日韩一区二区三区影片| av免费观看日本| 亚洲av在线观看美女高潮| 色吧在线观看| 亚洲av成人精品一区久久| 久久久久久国产a免费观看| 91狼人影院| 日韩av在线大香蕉| 亚洲经典国产精华液单| 在现免费观看毛片| 久久久亚洲精品成人影院| 99热6这里只有精品| 国产真实伦视频高清在线观看| 内地一区二区视频在线| 精品午夜福利在线看| 久久精品熟女亚洲av麻豆精品 | 日韩欧美国产在线观看| 精品一区二区三区视频在线| 日韩成人av中文字幕在线观看| 欧美97在线视频| 免费电影在线观看免费观看| 99热这里只有精品一区| 哪个播放器可以免费观看大片| 亚洲精品成人av观看孕妇| 在线观看av片永久免费下载| 大片免费播放器 马上看| 国产探花极品一区二区| 97在线视频观看| 午夜精品一区二区三区免费看| 亚洲综合色惰| 成人欧美大片| 免费观看性生交大片5| 美女脱内裤让男人舔精品视频| 国产精品熟女久久久久浪| 蜜臀久久99精品久久宅男| 女的被弄到高潮叫床怎么办| 免费电影在线观看免费观看| 人人妻人人看人人澡| 国产高清有码在线观看视频| 天堂中文最新版在线下载 | 亚洲国产色片| 自拍偷自拍亚洲精品老妇| 亚洲精品乱码久久久v下载方式| 免费电影在线观看免费观看| 久久久欧美国产精品| 国精品久久久久久国模美| 亚洲国产高清在线一区二区三| 午夜激情福利司机影院| 美女主播在线视频| 国产免费又黄又爽又色| 午夜亚洲福利在线播放| 高清午夜精品一区二区三区| 人妻系列 视频| 成人国产麻豆网| 联通29元200g的流量卡| 国产黄片视频在线免费观看| 成年女人看的毛片在线观看| 51国产日韩欧美| 少妇被粗大猛烈的视频| 日韩视频在线欧美| 久久99热这里只频精品6学生| 免费看光身美女| 一级二级三级毛片免费看| 欧美日韩亚洲高清精品| 日韩在线高清观看一区二区三区| 激情 狠狠 欧美| 国产成人午夜福利电影在线观看| 欧美另类一区| 亚洲性久久影院| 一区二区三区四区激情视频| 免费无遮挡裸体视频| 精品久久久久久久末码| 亚洲欧美精品自产自拍| 乱系列少妇在线播放| 人妻夜夜爽99麻豆av| 久久久久精品性色| 1000部很黄的大片| 亚洲精品日本国产第一区| 久久久亚洲精品成人影院| 麻豆av噜噜一区二区三区| 国产真实伦视频高清在线观看| 亚洲成人一二三区av| 美女大奶头视频| 精品人妻一区二区三区麻豆| 亚洲18禁久久av| 亚洲在线观看片| 国产亚洲午夜精品一区二区久久 | 日韩视频在线欧美| 高清av免费在线| 伊人久久精品亚洲午夜| 最近中文字幕高清免费大全6| 国产午夜精品论理片| 美女xxoo啪啪120秒动态图| 日日摸夜夜添夜夜添av毛片| 久久久精品欧美日韩精品| 国产伦理片在线播放av一区| 最近视频中文字幕2019在线8| 97超碰精品成人国产| av在线天堂中文字幕| 色5月婷婷丁香| 在线观看美女被高潮喷水网站| 搞女人的毛片| 视频中文字幕在线观看| 少妇熟女aⅴ在线视频| 少妇猛男粗大的猛烈进出视频 | 午夜福利在线观看免费完整高清在| 能在线免费看毛片的网站| av.在线天堂| 麻豆成人午夜福利视频| 亚洲激情五月婷婷啪啪| 建设人人有责人人尽责人人享有的 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲美女搞黄在线观看| 国产精品日韩av在线免费观看| 亚洲av二区三区四区| 最近2019中文字幕mv第一页| 久久99热这里只有精品18| 少妇的逼水好多| 一区二区三区乱码不卡18| 婷婷色麻豆天堂久久| 国产综合懂色| 一个人免费在线观看电影| 成人毛片60女人毛片免费| 亚洲人成网站在线播| 少妇的逼好多水| 成人亚洲精品一区在线观看 | 免费大片18禁| 日本免费在线观看一区| 国产精品蜜桃在线观看| 一级a做视频免费观看| 夫妻性生交免费视频一级片| 国产亚洲一区二区精品| 国产精品人妻久久久影院| 九草在线视频观看| 三级毛片av免费| 一级毛片 在线播放| 久久久久免费精品人妻一区二区| 欧美xxxx黑人xx丫x性爽| 免费看a级黄色片| 国产综合懂色| 亚洲精品久久午夜乱码| 精品人妻偷拍中文字幕| videos熟女内射| 久久精品久久精品一区二区三区| 中文天堂在线官网| 精品人妻一区二区三区麻豆| 欧美日韩视频高清一区二区三区二| 两个人的视频大全免费| 直男gayav资源| 一级二级三级毛片免费看| 国产免费一级a男人的天堂| 亚洲精品乱码久久久v下载方式| 狂野欧美激情性xxxx在线观看| 精品熟女少妇av免费看| 如何舔出高潮| 色综合站精品国产| 亚洲国产高清在线一区二区三| 亚洲国产成人一精品久久久| 久久精品综合一区二区三区| 久久久久久久久久人人人人人人| 婷婷色av中文字幕| 久久人人爽人人爽人人片va| av福利片在线观看| 一级av片app| 久久久久久国产a免费观看| 一级毛片黄色毛片免费观看视频| 男人狂女人下面高潮的视频| 少妇的逼好多水| 美女脱内裤让男人舔精品视频| 国产精品蜜桃在线观看| 夫妻午夜视频|