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

    鰳魚全基因組survey分析及微衛(wèi)星位點分布

    2024-01-01 00:00:00高魯修陳詩逸馮桃波劉炳艦柳意樊沈豪迪黃文化梁旭東
    水產(chǎn)科技情報 2024年5期
    關(guān)鍵詞:高通量測序

    摘" 要:為研究鰳魚(Ilisha elongata)的遺傳演化概況,對鰳魚基因組進(jìn)行了高通量測序并開展了survey分析和基因組范圍內(nèi)的遺傳標(biāo)記開發(fā)。K-mer分析結(jié)果顯示,鰳魚基因組大小約為692.8 Mb,雜合率為0.18%,重復(fù)序列比例為39.6%?;蚪M初步組裝Scaffold N50為27 694 bp,Contig N50為6 306 bp。利用MISA軟件對鰳魚基因組的微衛(wèi)星標(biāo)記(simple sequence repeat,SSR)進(jìn)行檢索和分析,總共檢測到786 123個SSR位點,相對豐度為1 204個/Mb。在所有類型中,重復(fù)最多的是二堿基,占SSR總量的75.94%,其SSR重復(fù)頻率主要集中在6~28次;其次為單堿基和四堿基,分別占SSR總量的15.31%和4.25%。對二堿基類型而言,AC型具有最多的重復(fù)數(shù)量,為131 732個,單堿基重復(fù)數(shù)量最多的則是A型,有48 775個。研究結(jié)果表明,鰳魚基因組經(jīng)組裝后可得到高質(zhì)量的全基因組序列,經(jīng)過篩選的SSR位點可為后續(xù)的遺傳分子標(biāo)記開發(fā)提供有力支持,研究結(jié)果可為鰳魚種質(zhì)資源管理和保護(hù)、生物進(jìn)化和群體遺傳等研究工作提供基礎(chǔ)資料。

    關(guān)鍵詞:鰳魚;全基因組survey;微衛(wèi)星標(biāo)記;高通量測序;生物遺傳

    doi:10.16446/j.fsti.20230700110

    收稿日期:2023-07-10

    作者簡介:高魯修(1998—),男,碩士研究生,研究方向為分子生態(tài)學(xué)。E-mail:736857968@qq.com

    通信作者:劉炳艦(1984—),男,副教授,研究方向為海洋生物學(xué)分子生態(tài)學(xué)。E-mail:liubingjian@zjou.edu.cn

    項目資助:國家自然科學(xué)基金項目(41806156);國家級大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計劃項目(202110340031);浙江省自然科學(xué)基金項目" (LY22D060001amp;LY20C190008);舟山市科技項目(2020C21016)。

    目前,人類已進(jìn)入了組學(xué)時代并見證了空前規(guī)模的遺傳信息積累。為深入了解和研究遺傳信息,研究人員采用了多種先進(jìn)的科學(xué)技術(shù),其中高通量測序技術(shù)[1](high-throughput sequencing)因其準(zhǔn)確度高、運(yùn)行成本低等特點而成為目前主流的測序技術(shù)之一,為大規(guī)模的組學(xué)研究提供了技術(shù)支持,并被稱作“下一代”測序技術(shù)(“next-generation” sequencing technology,NGS)。根據(jù)NCBI(national center for biotechnology information)數(shù)據(jù)庫(https://www.ncbi.nlm.nih.gov/genome)和IGSR(the international genome sample resource)數(shù)據(jù)庫(https://www.internationalgenome.org)中的數(shù)據(jù)顯示,截至目前,全球約有100 000種生物進(jìn)行了基因組測序,其中魚類有1 200多種?;贜GS的基因組survey分析,是指利用小片段文庫的低深度測序數(shù)據(jù)(50~100X),運(yùn)用K-mer分析方法對基因組展開研究,包括評估基因組大小、雜合度、GC堿基含量等?;蚪Msurvey分析是了解某一生物遺傳特性的有效方法,相關(guān)結(jié)果可為后續(xù)的建庫及高質(zhì)量基因組組裝策略的制定提供資料,并對全基因組微衛(wèi)星標(biāo)記(simple sequence repeats,SSR)的開發(fā)發(fā)揮重要作用。研究表明,SSR均勻分布在真核生物基因組中,目前已被應(yīng)用于種群遺傳多樣性、親緣關(guān)系識別、構(gòu)建遺傳圖譜等研究領(lǐng)域[2-4]。相對于常規(guī)的SSR開發(fā)流程,高通量測序數(shù)據(jù)下的SSR開發(fā)工作具備低成本、高效率、高生產(chǎn)率的優(yōu)勢。趙蕊蕊等[5]對絨杜父魚(Hemitripterus villosus)全基因組進(jìn)行了分析,獲得約710 Mb的基因組序列,其雜合度和重復(fù)序列占比分別為38.61%和0.26%,鑒定獲得SSR位點583 498個,為絨杜父魚進(jìn)化生物學(xué)和遺傳學(xué)等研究提供了基礎(chǔ)資料。Kim等[6]對南極冰魚(Chionobathysscus dewitti)基因組測序數(shù)據(jù)進(jìn)行了評估,得出南極冰魚基因組大小為880 Mb,其誤差率、雜合度和重復(fù)率依次為0.317%、0.421%和0.738%,同時從組裝數(shù)據(jù)中鑒定出2 252 265個SSR位點,其中二堿基重復(fù)序列比例最高。Xu等[7]對褐菖鲉(Sebastiscus marmoratus)的全基因組進(jìn)行了分析,獲得800 Mb左右的基因組序列,其中重復(fù)序列為39.65%,雜合度為0.17%,并通過全基因組數(shù)據(jù)開發(fā)出191 592個SSR位點,為其他相關(guān)研究提供了參考。Li等[8]的研究顯示,羅布麻(Apocynum venetum)基因組大小為254.40 Mb,重復(fù)序列占比為40.87%,雜合度約為0.63%,同時挖掘出101 918個SSR位點,為今后羅布麻遺傳信息識別以及分析棉麻纖維合成的關(guān)鍵基因等工作提供了基礎(chǔ)資料。

    鰳魚(Ilisha elongata)又稱白鰳魚,隸屬于鯡形目(Clupeiformes)、鋸腹鰳科(Pristigasteridae),主要分布于太平洋西部、爪哇海海域[9],在中國南海、東海、黃海以及渤海均有分布[10],為亞熱帶及暖溫帶近海洄游性中上層魚類,主要以蝦類、魚類和毛顎類等為食,是我國沿海重要的經(jīng)濟(jì)魚類[11]。鰳魚肉質(zhì)鮮美,富含蛋白質(zhì)、脂肪及鈣等營養(yǎng)物質(zhì),是我國漁業(yè)史上最早的捕撈對象之一,距今已有5 000多年的歷史[12]。近年來,由于過度捕撈以及海洋生態(tài)環(huán)境惡化[13]等因素,鰳魚資源量急劇下降,因此亟需開展鰳魚種質(zhì)資源調(diào)查和保護(hù)工作。目前,有關(guān)鰳魚的研究工作主要集中在食品加工[14-19]和形態(tài)特征[20-23]等方面,有關(guān)遺傳分析方面的研究[24-26]相對較少,限制了鰳魚生物進(jìn)化、遺傳多樣性研究和種質(zhì)資源保護(hù)等工作的開展。本研究對鰳魚進(jìn)行高通量測序,展開全基因組的survey分析和識別SSR,對鰳魚基因組大小、雜合度、GC含量等進(jìn)行有效評估,以期為后續(xù)鰳魚資源管理和保護(hù)、生物進(jìn)化和遺傳多樣性研究方法的選用提供基本資料。

    1" 材料和方法

    1.1" 樣本采集與測序

    本研究鰳魚樣品于2023年2月采自浙江舟山近海水域。經(jīng)形態(tài)學(xué)特征[27]進(jìn)行物種鑒定后,剪取適量樣品肌肉組織裝入含有95%(體積分?jǐn)?shù))乙醇溶液的凍存管中固定,然后于-80 ℃冰箱中保存待用。

    采用常規(guī)方法(苯酚-氯仿法)提取鰳魚基因組DNA,并委托北京百邁客生物科技有限公司進(jìn)行高通量測序文庫構(gòu)建和測序,原始數(shù)據(jù)(raw data)經(jīng)質(zhì)控和過濾后,所獲得的有效數(shù)據(jù)(clean data)可用于后續(xù)的相關(guān)分析。將原始數(shù)據(jù)上傳至NCBI數(shù)據(jù)庫,注冊登錄號為PRJNA967191。

    1.2" 全基因組數(shù)據(jù)處理

    測序后得到的鰳魚有效數(shù)據(jù)通過GCE 1.0.2軟件的K-mer分析工具來評估其基因組的大小、雜合度、GC含量等基本特征信息[28]。利用SOAPdenovo 2.0軟件進(jìn)行基因組初步組裝[29],并計算Contig和Scaffold(K=47)大小?;蚪M組裝完畢后,使用MISA腳本的默認(rèn)參數(shù)值(即單堿基至六堿基的最小重復(fù)序列次數(shù)分別為10、6、5、5、5和5)對序列中的SSR位點進(jìn)行篩選。利用Primer 3軟件進(jìn)行SSR引物批量設(shè)計[30],引物委托生工生物工程(上海)股份有限公司制備。經(jīng)聚合酶鏈?zhǔn)椒磻?yīng)(PCR)擴(kuò)增后,采用瓊脂糖凝膠電泳分析法對引物和SSR的有效性進(jìn)行檢驗。

    2" 結(jié)果

    2.1" 數(shù)據(jù)量統(tǒng)計與K-mer分析

    利用Illumina NovaSeq 6000平臺測序,經(jīng)質(zhì)控和過濾后得到有效數(shù)據(jù)為155.80 Gb,其中GC含量為42.92%,錯誤率為0.02%,Q20為99.57%,Q30為97.39%,測序質(zhì)量符合后續(xù)分析要求(見表1)。通過與NCBI數(shù)據(jù)庫中相關(guān)物種[太平洋鰳(Ilisha furthii)、縱帶鰳(Ilisha striatula)、西非鰳(Ilisha africana)等]序列進(jìn)行BLASTn分析,比對結(jié)果無異常,表明測序數(shù)據(jù)不存在污染情況。使用K-mer分析評估有效數(shù)據(jù),設(shè)定K-mer值為18,基因組的預(yù)估結(jié)果為:鰳魚基因組大小約為692.8 Mb,重復(fù)序列占比39.6%,雜合率0.18%,可對其展開進(jìn)一步研究和分析,包括基因組裝和SSR篩選等工作。

    2.2" 基因組初步組裝

    鰳魚全基因組初步組裝結(jié)果見表2。利用KmerGenie軟件進(jìn)行有效數(shù)據(jù)分析[31],得出最佳基因組組裝值K=47。通過SOAPdenovo軟件初步組裝鰳魚全基因組,最終獲得基因組大小約為556 914 471 bp。進(jìn)一步分析結(jié)果顯示,Scaffold N50長度為27 694 bp,Contig N50長度為6 306 bp。Scaffolds最長序列為47 402 bp,有110 255個序列不小于2 000 bp,有536 946個序列不小于100 bp,表明鰳魚全基因組初步組裝結(jié)果較好。

    2.3" 微衛(wèi)星位點篩選

    基因組初步組裝完成后,通過MISA軟件挖掘SSR位點,在全部Scaffolds中篩選獲得SSR位點786 123個,相對豐度為1 204個/Mb。進(jìn)一步分析結(jié)果顯示,在不同的堿基類型中,二堿基重復(fù)類型最多(596 970個),約占SSR總量的75.94%,六堿基重復(fù)最少(880個),約占SSR總量的0.11%。這與趙蕊蕊等[5]、Tóth等[32]、王耀嶸等[33]的研究結(jié)果類似,即在堿基類型所占比例中,二堿基占比最高,六堿基占比最低。其余的堿基類型占比依次是單堿基(120 368個,15.31%),四堿基(33 347個,4.25%),三堿基(28 249個,3.59%)和五堿基(6 309個,0.80%)。對SSR位點進(jìn)行分析,其重復(fù)范圍在5~103次,主要集中次數(shù)為6~28次。對單堿基而言,重復(fù)次數(shù)大多集中于10~28次,約占單堿基SSR總量的99.30%。二、三堿基重復(fù)主要集中在6~28次和5~15次,約占其SSR總量的99.58%和99.63%。四堿基、五堿基和六堿基的重復(fù)次數(shù)則分別集中于5~11次、5~9次和5~7次,約各占SSR總量的98.72%、99.06%和97.84%(見圖1)。

    鰳魚全基因組SSR重復(fù)次數(shù)分布情況見圖2。在單堿基重復(fù)類型中,A/T類型與C/G類型相比重復(fù)單元較高,為93 801個,約占單堿基SSR總量的77.93%,C/G類型為26 567個,約占單堿基SSR總量的22.07%。二堿基重復(fù)類型中,以AC/GT型重復(fù)單元最多,為449 474個,約占二堿基SSR總量的75.29%,占比第二的是AG/CT型,為136 861個,約占22.93%,其余依次是AT/AT型(10 123個,1.70%)、CG/CG型(512個,0.09%)。

    三堿基結(jié)果顯示(見圖3),AAT/ATT型重復(fù)單元最多,為10 347個,約占三堿基SSR總量的36.63%,重復(fù)單元最少的是ATC/ATG型,為2 629個,約占其總量的9.31%,其次是AGG/CCT型(4 261個,15.08%)、AAC/GTT型(2 685個,9.51%)。對四堿基重復(fù)而言,以AGAT/ATCT型重復(fù)單元最多,為7 756個,約占四堿基SSR總量的23.26%,其余依次是ACAG/CTGT型(4 863個,14.58%)、AAAG/CTTT型(4 276個,12.82%)、ACTC/AGTG型(2 853個,8.56%)。五堿基和六堿基中,重復(fù)單元最多的是AGAGG/CCTCT型(1 491個,23.63%)和ACACGC/CGTGTG型(86個,9.77%)(見圖4)。

    通過PCR擴(kuò)增和瓊脂糖凝膠電泳檢測法對隨機(jī)選取的15個SSR位點進(jìn)行有效性驗證,結(jié)果顯示,PCR擴(kuò)增成功率在95%以上,滿足后續(xù)的研究工作要求。

    3" 討論

    3.1" 鰳魚全基因組的基本特征

    NGS的發(fā)展為基因組和轉(zhuǎn)錄組研究提供了一種相對經(jīng)濟(jì)的方法。本研究通過survey分析對鰳魚全基因組數(shù)據(jù)進(jìn)行評估,開發(fā)出了豐富的SSR標(biāo)記,研究結(jié)果可為后續(xù)鰳魚漁業(yè)資源管理保護(hù)和生物遺傳等研究工作提供基礎(chǔ)資料。

    K-mer分析顯示,鰳魚的基因組大小約為692.8 Mb,小于南極冰魚[6](880 Mb)、斑魚祭(Setipinna tenuifilis)[34](797 Mb)、黃鯽(Setipinna tenuifilis)[35](815 Mb)和大西洋鯡(Clupea harengus)[36](850 Mb),與歐洲沙丁魚(Sardina pilchardus)[37]

    (625~637 Mb)、卵形鯧鲹(Trachinotus ovatus)[38](642.68 Mb)相似?;蚪M大小的不同是由物種間的差異性造成的,并且可能與物種的基因組重復(fù)比例相關(guān)。鰳魚基因組的重復(fù)序列比例為39.6%,處于中低水平,高于絨杜父魚[5](38.61%)和卵形鯧鲹[38](30.19%),與褐菖鲉[7](39.65%)、斑魚祭[34](39.22%)和黃鯽[35](39.69%)相似,但低于條紋斑竹鯊(Chiloscyllium plagiosum)[39](63.53%)。此外,鰳魚測序數(shù)據(jù)的雜合率為0.18%,GC含量為42.92%,低于絨杜父魚[5](0.26%,43.13%)和南極冰魚[6](0.412%,49.9%),但與褐菖鲉[7](0.17%,41.3%)和斑尾刺蝦虎魚(Acanthogobius ommaturus)[40](0.17%,40.88%)相當(dāng)。根據(jù)基因組大小、重復(fù)序列比、雜合度等survey分析指標(biāo)判斷,鰳魚的基因組應(yīng)屬于簡單型。這一結(jié)果可為后續(xù)鰳魚的高質(zhì)量基因組組裝工作提供參考。

    本研究初步組裝結(jié)果顯示,Scaffold N50大小為27 694 bp,Contig N50大小為6 306 bp。組裝后的Scaffold大小為556 914 471 bp,數(shù)目為536 946條,平均長度為1 038 bp。其長度與絨杜父魚[5](577 386 707 bp)相當(dāng),但低于南極冰魚[6](897 784 561 bp)、褐菖鲉[7](609 456 819 bp)、斑魚祭[34](800 444 663 bp)和黃鯽[35](798 382 266 bp)等??紤]到組裝的N50和N90長度較短,可能會導(dǎo)致鰳魚基因組初步組裝質(zhì)量、連續(xù)性和完整性較低,不能滿足其生物進(jìn)化相關(guān)的研究工作,建議通過多種技術(shù)相結(jié)合的方法(例如Illumina測序與PacBio和Hi-C相結(jié)合)來進(jìn)行基因組序列的構(gòu)建。盡管如此,本次組裝得到的序列仍然可以用于后續(xù)SSR位點、SNP位點開發(fā)以及鰳魚群體基因組學(xué)等研究。

    3.2" 鰳魚全基因組SSR分布特征

    重復(fù)序列比例(鰳魚為39.6%)通常被認(rèn)為與基因組所含的SSR數(shù)量成正比,本試驗從組裝好的鰳魚Scaffolds中共篩選出786 123個SSR位點,與其他物種[5-8,33-35,38-41]相比處于中等水平。其中,重復(fù)最少的是六堿基,僅為SSR總量的0.11%(880個),二堿基則重復(fù)最高,為596 970個,約占SSR總量的75.94%。這與其他海洋魚類,如絨杜父魚[5]、南極冰魚[6]、褐菖鲉[7]、金錢魚(Scatophagus argus)[33]、斑魚祭[34]、黃鯽[35]、卵形鯧鲹[38]、條紋斑竹鯊[39]、斑尾刺蝦虎魚[40]和巨魚丕(Bagarius yarrelli)[41]等的研究結(jié)果相一致,均體現(xiàn)出二堿基重復(fù)數(shù)量較多、六堿基最少的特征,可以作為后續(xù)相關(guān)研究的參考資料。在單堿基重復(fù)中以A/T類型為主,與前人[5-7,33-35,38-41]的研究結(jié)果相同。相關(guān)研究表明,A/T堿基含量與物種基因組中SSR數(shù)量呈正相關(guān),即A/T堿基所占比例越高,SSR數(shù)量越多。此外,還可能與DNA聚合酶在DNA復(fù)制過程中的錯誤率和甲基化的CpG位點相關(guān)。二堿基重復(fù)中,AC/GT類型占主要優(yōu)勢,與絨杜父魚[5]、南極冰魚[6]、條紋斑竹鯊[39]一致,但與褐菖鲉[7]、卵形鯧鲹[38](AC/TG為主)和巨魚丕[41]、斑點叉尾魚回(Ictalurus punctatus)[42](AC/AG為主)略有不同。這可能與物種差異或堿基重復(fù)頻率有關(guān),但與大多數(shù)魚類二堿基類型中都是以AC型為主[43]這一觀點相吻合。特定堿基重復(fù)的相對豐度在重復(fù)序列中具有較高的變異性,鰳魚全基因組中,二堿基、三堿基、四堿基、五堿基和六堿基的重復(fù)頻率分布范圍分別為6~28次、5~15次、5~11次、5~9次和5~7次,表明堿基重復(fù)次數(shù)隨重復(fù)長度的增加而逐漸減少,這可能是由重復(fù)長度越長、突變率越高所導(dǎo)致的。Wierdl等[44]和Chen等[45]也報道了重復(fù)次數(shù)與重復(fù)長度呈負(fù)相關(guān)的研究結(jié)果,本研究結(jié)果也證實了這一模式。

    4" 結(jié)論

    利用高通量測序技術(shù)結(jié)合survey分析得出,鰳魚基因組大小約為692.8 Mb,重復(fù)序列比例為39.6%,雜合率為0.18%,GC含量為42.92%,錯誤率為0.02%,組裝后的Scaffold大小為556 914 471 bp。識別鰳魚全基因組SSR序列共786 123個,蘊(yùn)含豐富的遺傳信息,不僅有助于開展鰳魚種質(zhì)資源評估保護(hù)和遺傳分子標(biāo)記優(yōu)選等工作,而且可以為其數(shù)量性狀基因座(QTL)定位、種群基因分析、遺傳基因多樣化水平鑒別等工作提供參考。

    參考文獻(xiàn)

    [1]SCHUSTER S C. Next-generation sequencing transforms today’s biology[J].Nature Methods,2008,5(1):16-18.

    [2]SONG H,GUO X M,SUN L N,et al. The hard clam genome reveals massive expansion and diversification of inhibitors of apoptosis in Bivalvia[J].BMC Biology,2021,19(1):15.

    [3]VIEIRA M L C,SANTINI L,DINIZ A L,et al. Microsatellite markers:what they mean and why they are so useful[J].Genetics and Molecular Biology,2016,39(3):312-328.

    [4]OLUBUNMI O O. Application of microsatellite in fish biotechnology:Prospects and drawback-review[J].Bioengineering and Biotechnology,2019,4(3):37-43.

    [5]趙蕊蕊,徐勝勇.絨杜父魚全基因組survey分析及微衛(wèi)星分布特征[J].中國水產(chǎn)科學(xué),2022,29(7):994-1001.

    [6]KIM J,LEE S J,JO E,et al.Whole-genome survey and microsatellite marker detection of Antarctic crocodile icefish,Chionobathyscus dewitti[J].Animals,2022,12(19):2598.

    [7]XU S Y,SONG N,XIAO S J,et al.Whole genome survey analysis and microsatellite motif identification of Sebastiscus marmoratus[J].Bioscience Reports,2020,40(2):20192252.

    [8]LI G Q,SONG L X,JIN C Q,et al.Genome survey and SSR analysis of Apocynum venetum[J].Bioscience Reports,2019,39(6):20190146.

    [9]BLABER S J M,STAUNTON-SMITH J,MILTON D A,et al.The biology and life-history strategies of Ilisha (Teleostei:Pristigasteridae) in the coastal waters and estuaries of Sarawak[J].Estuarine,Coastal and Shelf Science,1998,47(4):499-511.

    [10]謝漢陽.東海區(qū)兩種經(jīng)濟(jì)魚類的資源動態(tài)分析[D].寧波:寧波大學(xué),2011.

    [11]蔣日進(jìn),薛利建,張洪亮,等.東海區(qū)鰳魚的攝食習(xí)性[J].海洋漁業(yè),2013,35(2):168-175.

    [12]王世光,王瑞瑞,陳瑜,等.淺談舟山海鮮之舟山鰳魚[J].食品安全導(dǎo)刊,2019(15):173-174.

    [13]HALPERN B S,LONGO C,HARDY D,et al.An index to assess the health and benefits of the global ocean[J].Nature,2012,488(7413):615-620.

    [14]陳慧,黃建聯(lián),周文果,等.鰳魚魚糜與常見魚糜凝膠品質(zhì)對比[J].食品與機(jī)械,2022,38(5):19-23.

    [15]吳佳佳,王思齊,戴志遠(yuǎn).咸鰳魚中產(chǎn)組胺菌的分離與鑒定[J].食品與發(fā)酵工業(yè),2019,45(3):71-76.

    [16]劉檾菘,湯海青,歐昌榮.固態(tài)發(fā)酵過程中鰳魚品質(zhì)及抗氧化能力的研究[J].安徽農(nóng)業(yè)大學(xué)學(xué)報,2018,45(2):214-218.

    [17]湯海青,歐昌榮,宋玉玲,等.清酒乳桿菌固態(tài)發(fā)酵鰳魚過程中的品質(zhì)變化[J].核農(nóng)學(xué)報,2017,31(10):1938-1945.

    [18]劉敏,夏松養(yǎng).腌制條件對鰳魚腌制品品質(zhì)的影響研究[J].安徽農(nóng)業(yè)科學(xué),2014,42(35):12655-12657.

    [19]賀林娟.發(fā)酵鰳魚的制作工藝及揮發(fā)性風(fēng)味成分的研究[D].寧波:寧波大學(xué),2013.

    [20]盧振彬,周金福,顏尤明,等.閩江口鰳魚生殖群體的生態(tài)特性[J].水產(chǎn)科技情報,1998,25(4):158-162.

    [21]倪海兒.東海鰳魚雌雄性狀比較及其鑒別[J].浙江水產(chǎn)學(xué)院學(xué)報,1994,13(2):93-98.

    [22]陳必哲,張澄茂.廈門近海鰳魚生殖群體生長與資源狀況[J].福建水產(chǎn),1993,15(4):35-38.

    [23]陳大剛,王啟華,劉長安.山東南部沿海鰳魚產(chǎn)卵群體的生物學(xué)特征[J].海洋通報,1988,7(1):49-55.

    [24]王倩.基于核基因數(shù)據(jù)重建鯡形目魚類系統(tǒng)發(fā)育關(guān)系和太平洋西部鰳魚群體遺傳多樣性分析[D].上海:上海海洋大學(xué),2019.

    [25]李鵬飛,賀舟挺,徐開達(dá),等.基于COⅠ基因和D-loop序列的東海鰳魚種質(zhì)資源遺傳變異研究[J].浙江海洋學(xué)院學(xué)報(自然科學(xué)版),2013,32(2):93-98.

    [26]呂振明,許逸天,吳常文,等.中國沿海鰳不同地理群體16S rRNA基因的遺傳變異分析[J].中國水產(chǎn)科學(xué),2010,17(3):463-470.

    [27]陳大剛,張美昭.中國海洋魚類[M].青島:中國海洋大學(xué)出版社,2015.

    [28]LIU B H,SHI Y J,YUAN J Y,et al.Estimation of genomic characteristics by analyzing k-mer frequency in de novo genome projects[J].arXiv preprint:1308.2012,2013.

    [29]LUO R B,LIU B H,XIE Y L,et al.SOAPdenovo2:an empirically improved memory-efficient short-read de novo assembler[J].GigaScience,2012,1(1):18.

    [30]ROZEN S,SKALETSKY H.Primer3 on the WWW for general users and for biologist programmers[M].Bioinformatics Methods and Protocols.New Jersey:Humana Press,2003.

    [31]CHIKHI R,MEDVEDEV P.Informed and automated k-mer size selection for genome assembly[J].Bioinformatics,2014,30(1):31-37.

    [32]TTH G,GSPRI Z,JURKA J.Microsatellites in different eukaryotic genomes:survey and analysis[J].Genome Research,2000,10(7):967-981.

    [33]王耀嶸,楊尉,任席林,等.金錢魚基因組微衛(wèi)星分布特征分析及多態(tài)性標(biāo)記開發(fā)[J].廣東海洋大學(xué)學(xué)報,2020,40(4):7-14.

    [34]張坤.斑魚祭全基因組組裝及基于微衛(wèi)星的群體遺傳學(xué)研究[D].舟山:浙江海洋大學(xué),2022.

    [35]LIU B J,LI J S,PENG Y,et al.Chromosome-level genome assembly and population genomic analysis reveal evolution and local adaptation in common hairfin anchovy(Setipinna tenuifilis)[J].Molecular Ecology,2024,33(10):17067.

    [36]BARRIO A M,LAMICHHANEY S,F(xiàn)AN G Y,et al.The genetic basis for ecological adaptation of the Atlantic herring revealed by genome sequencing[J].eLife,2016,5:12081.

    [37]MACHADO A M,TRRESEN O K,KABEYA N,et al.“Out of the can”:a draft genome assembly,liver transcriptome,and nutrigenomics of the European sardine,Sardina pilchardus[J].Genes,2018,9(10):485.

    [38]張永德,文露婷,羅洪林,等.卵形鯧鲹基因組調(diào)研及其SSR分子標(biāo)記的開發(fā)應(yīng)用[J].南方農(nóng)業(yè)學(xué)報,2020,51(5):983-994.

    [39]ZHAO R R,CAI S S,LU D F,et al.Genomic comparison and genetic marker identification of the white-spotted bamboo shark Chiloscyllium plagiosum[J].Frontiers in Marine Science,2022,9:936681.

    [40]CHEN B J,SUN Z C,LOU F R,et al.Genomic characteristics and profile of microsatellite primers for Acanthogobius ommaturus by genome survey sequencing[J].Bioscience Reports,2020,40(11):20201295.

    [41]楊汶珊,唐榮葉,蘇孟園,等.巨魚丕(Bagarius yarrelli)全基因組微衛(wèi)星分布特征分析[J].南京師范大學(xué)學(xué)報(工程技術(shù)版),2021,21(3):62-68.

    [42]唐榮葉,蘇孟園,楊汶珊,等.斑點叉尾魚回全基因組微衛(wèi)星分布特征分析[J].漁業(yè)科學(xué)進(jìn)展,2022,43(2):89-97.

    [43]劉士力,劉一諾,李飛,等.3種鲌亞科魚類全基因組微衛(wèi)星分布特征分析[J].上海海洋大學(xué)學(xué)報,2023,32(6):1176-1183

    [44]WIERDL M,DOMINSKA M,PETES T D.Microsatellite instability in yeast:dependence on the length of the microsatellite[J].Genetics,1997,146(3):769-779.

    [45]CHEN M,TAN Z Y,ZENG G M,et al.Comprehensive analysis of simple sequence repeats in pre-miRNAs[J].Molecular Biology and Evolution,2010,27(10):2227-2232.

    Whole-genome analysis and microsatellite distribution of Ilisha elongata

    GAO Luxiu1, CHEN Shiyi2, FENG Taobo2, LIU Bingjian2, LIU Yifan2, SHEN Haodi2, HUANG Wenhua2, LIANG Xudong2

    (1. National Engineering Research Center for Marine Aquaculture,Zhejiang Ocean University,Zhoushan" 316022,China;

    2. Marine Science and Technology College,Zhejiang Ocean University,Zhoushan" 316022,China)

    Abstract: To study the genetic evolutionary profile of Ilisha elongata, high-throughput sequencing was performed on the I. elongata genome,and survey analysis and genome-wide genetic marker development were conducted.Results showed that the genome size of I. elongata,as determined by K-mer analysis,was approximately 692.8 Mb,with a heterozygosity rate of 0.18% and the repeated sequence rate of 39.6%.The preliminary genome assembly results showed that the Scaffold N50 was 27 694 bp and the Contig N50 was 6 306 bp.The simple sequence repeats(SSRs) of I. elongata genome were searched and analyzed using MISA software,and a total of 786 123 SSR loci were detected, with a relative abundance of 1 204 per Mb.Among all types,dinucleotides were the most abundant repeats,which accounting for 75.94% of the total SSRs and the SSR repeat frequencies mainly concentrated between 6 and 28 times.This was followed by mononucleotide and tetranucleotide repeats,accounting for 15.31% and 4.25% of the total SSRs,respectively.For dinucleotide repeats,type AC had the highest number of repeats at 131 732,while type A had the highest number of mononucleotide repeats at 48 775.The results showed that the genome of I. elongata could be assembled to obtain high-quality whole genome sequences,and the SSR loci could provide strong support for the subsequent genetic molecular marker development,as well as basic information for research work on germplasm resource management and conservation,biological evolution and population genetics of I. elongata.

    Key words: Ilisha elongata; whole-genome survey; microsatellite marker; high-throughput sequencing; biogenetics

    猜你喜歡
    高通量測序
    基于高通量測序的野生毛葡萄轉(zhuǎn)錄組SSR信息分析
    木質(zhì)纖維素分解復(fù)合菌系的分解特性與細(xì)菌組成多樣性分析
    基于高通量測序技術(shù)對三種太歲樣品細(xì)菌組成的分析
    基于高通量測序的玄參根部轉(zhuǎn)錄組學(xué)研究及萜類化合物合成相關(guān)基因的挖掘
    污水處理中壓力變化對污泥中微生物群落組成的影響研究
    石柱黃連根腐病根際土壤細(xì)菌微生態(tài)研究
    環(huán)狀RNA在疾病發(fā)生中的作用
    川明參輪作對煙地土壤微生物群落結(jié)構(gòu)的影響
    多穗柯轉(zhuǎn)錄組分析及黃酮類化合物合成相關(guān)基因的挖掘
    人參根際真菌群落多樣性及組成的變化
    日本-黄色视频高清免费观看| 亚洲伊人久久精品综合| 爱豆传媒免费全集在线观看| 亚洲精品国产一区二区精华液| 久久久精品94久久精品| 校园人妻丝袜中文字幕| 一本久久精品| 欧美成人精品欧美一级黄| 老熟女久久久| 亚洲国产精品999| 国产男女超爽视频在线观看| av在线观看视频网站免费| 日本色播在线视频| 一边亲一边摸免费视频| 好男人视频免费观看在线| 国产综合精华液| 日本-黄色视频高清免费观看| 中文字幕av电影在线播放| 大片电影免费在线观看免费| 免费在线观看黄色视频的| 超色免费av| 国产女主播在线喷水免费视频网站| 国产av码专区亚洲av| 亚洲精品久久成人aⅴ小说| 日本91视频免费播放| 一级a爱视频在线免费观看| 日韩三级伦理在线观看| 国产精品亚洲av一区麻豆 | 亚洲美女搞黄在线观看| 午夜91福利影院| 搡老乐熟女国产| 青春草亚洲视频在线观看| 七月丁香在线播放| 亚洲国产av影院在线观看| 五月伊人婷婷丁香| 香蕉丝袜av| 久久精品国产亚洲av涩爱| 美女午夜性视频免费| 成年av动漫网址| 免费黄频网站在线观看国产| 激情视频va一区二区三区| 九色亚洲精品在线播放| av不卡在线播放| 亚洲国产最新在线播放| 美女福利国产在线| av国产精品久久久久影院| 丰满乱子伦码专区| 久久久久久伊人网av| 一级毛片 在线播放| 欧美日韩亚洲国产一区二区在线观看 | 一区二区日韩欧美中文字幕| 天天影视国产精品| 国产精品.久久久| av福利片在线| 91成人精品电影| 久久综合国产亚洲精品| 午夜免费男女啪啪视频观看| 丰满乱子伦码专区| 亚洲第一青青草原| 久久av网站| 一级毛片我不卡| av国产久精品久网站免费入址| 国产成人精品婷婷| 欧美精品国产亚洲| 亚洲内射少妇av| 日韩中文字幕视频在线看片| 亚洲色图 男人天堂 中文字幕| 亚洲久久久国产精品| 丝袜在线中文字幕| 亚洲男人天堂网一区| 伦理电影大哥的女人| 亚洲av国产av综合av卡| 成人18禁高潮啪啪吃奶动态图| 成年人午夜在线观看视频| 18禁动态无遮挡网站| 亚洲第一av免费看| 男女边摸边吃奶| 国产成人免费观看mmmm| 中文字幕av电影在线播放| 国产成人aa在线观看| 精品亚洲乱码少妇综合久久| 欧美xxⅹ黑人| 亚洲一级一片aⅴ在线观看| 午夜福利影视在线免费观看| 男人舔女人的私密视频| 国产成人精品福利久久| 成人国产av品久久久| 免费黄网站久久成人精品| 不卡视频在线观看欧美| 久久精品亚洲av国产电影网| av天堂久久9| 春色校园在线视频观看| 中文字幕亚洲精品专区| 欧美黄色片欧美黄色片| 黄片无遮挡物在线观看| 两个人看的免费小视频| 建设人人有责人人尽责人人享有的| 国产爽快片一区二区三区| 国产成人午夜福利电影在线观看| freevideosex欧美| 黑人欧美特级aaaaaa片| 麻豆精品久久久久久蜜桃| 欧美精品高潮呻吟av久久| 国产欧美亚洲国产| 哪个播放器可以免费观看大片| 欧美少妇被猛烈插入视频| 99久国产av精品国产电影| 欧美成人午夜免费资源| 色94色欧美一区二区| 大香蕉久久成人网| 色哟哟·www| 亚洲精品国产av成人精品| 久久精品国产a三级三级三级| 老司机影院毛片| 欧美日韩视频精品一区| 亚洲成人手机| 三级国产精品片| 少妇精品久久久久久久| 日韩人妻精品一区2区三区| 一本—道久久a久久精品蜜桃钙片| 亚洲欧洲日产国产| 国产亚洲一区二区精品| 亚洲精品久久成人aⅴ小说| 日韩av在线免费看完整版不卡| 国产精品亚洲av一区麻豆 | 下体分泌物呈黄色| 少妇被粗大的猛进出69影院| 精品国产乱码久久久久久男人| 日韩大片免费观看网站| 高清欧美精品videossex| 自线自在国产av| 18禁国产床啪视频网站| 日产精品乱码卡一卡2卡三| 亚洲综合色网址| 一级毛片我不卡| 亚洲伊人色综图| 国产在线视频一区二区| 精品福利永久在线观看| 极品人妻少妇av视频| 卡戴珊不雅视频在线播放| 国产乱来视频区| 久久这里有精品视频免费| 亚洲欧美一区二区三区黑人 | 国产爽快片一区二区三区| 亚洲在久久综合| 2022亚洲国产成人精品| 成人漫画全彩无遮挡| 97人妻天天添夜夜摸| 亚洲第一av免费看| 亚洲av日韩在线播放| 极品人妻少妇av视频| 毛片一级片免费看久久久久| 亚洲精品一区蜜桃| 久久精品人人爽人人爽视色| 婷婷色麻豆天堂久久| 高清av免费在线| 日韩免费高清中文字幕av| 久久精品国产鲁丝片午夜精品| 亚洲中文av在线| 日本av免费视频播放| 国产视频首页在线观看| 国产亚洲一区二区精品| 久久久久精品人妻al黑| 搡女人真爽免费视频火全软件| 观看美女的网站| 久久精品夜色国产| 99国产精品免费福利视频| 亚洲精品久久久久久婷婷小说| 亚洲欧洲日产国产| 国产日韩欧美视频二区| 亚洲成av片中文字幕在线观看 | 99久久中文字幕三级久久日本| 黄色毛片三级朝国网站| 免费看av在线观看网站| 日韩,欧美,国产一区二区三区| 搡女人真爽免费视频火全软件| 亚洲精品av麻豆狂野| 久久精品国产亚洲av天美| 欧美激情极品国产一区二区三区| 美女国产视频在线观看| a级毛片黄视频| 美女中出高潮动态图| 欧美日韩一区二区视频在线观看视频在线| 国产综合精华液| 在线观看www视频免费| 亚洲成av片中文字幕在线观看 | 曰老女人黄片| 熟女av电影| 免费播放大片免费观看视频在线观看| 午夜91福利影院| 亚洲国产精品国产精品| 香蕉丝袜av| 国产97色在线日韩免费| 国产精品久久久久成人av| 男男h啪啪无遮挡| 巨乳人妻的诱惑在线观看| 尾随美女入室| 大香蕉久久成人网| 大片免费播放器 马上看| 黑人巨大精品欧美一区二区蜜桃| 亚洲av中文av极速乱| 成人18禁高潮啪啪吃奶动态图| 天堂俺去俺来也www色官网| 国产免费又黄又爽又色| 丝袜人妻中文字幕| 国产成人av激情在线播放| 中文字幕制服av| 色94色欧美一区二区| 美女国产视频在线观看| 五月天丁香电影| 亚洲三级黄色毛片| 这个男人来自地球电影免费观看 | 国产97色在线日韩免费| 少妇 在线观看| 亚洲精华国产精华液的使用体验| 免费观看a级毛片全部| 校园人妻丝袜中文字幕| 亚洲av在线观看美女高潮| 久久女婷五月综合色啪小说| 中国三级夫妇交换| 国产成人aa在线观看| 免费看不卡的av| 精品人妻偷拍中文字幕| 久久久久久久精品精品| videos熟女内射| 亚洲欧美一区二区三区黑人 | 精品一区二区三卡| 久久99精品国语久久久| 精品少妇久久久久久888优播| 超碰97精品在线观看| 日韩,欧美,国产一区二区三区| 免费播放大片免费观看视频在线观看| 9191精品国产免费久久| 美女午夜性视频免费| 国产精品 国内视频| 成人18禁高潮啪啪吃奶动态图| 久热久热在线精品观看| 叶爱在线成人免费视频播放| 国产亚洲欧美精品永久| 青春草国产在线视频| 欧美97在线视频| 一级片免费观看大全| 亚洲av福利一区| 免费在线观看视频国产中文字幕亚洲 | 亚洲,欧美,日韩| 亚洲av日韩在线播放| 国产一区二区在线观看av| 校园人妻丝袜中文字幕| 亚洲情色 制服丝袜| 精品国产国语对白av| 亚洲国产精品国产精品| 女性生殖器流出的白浆| 超色免费av| 在线看a的网站| 久久免费观看电影| 精品国产一区二区久久| 亚洲国产精品一区三区| 中文乱码字字幕精品一区二区三区| 国产成人aa在线观看| 国产精品蜜桃在线观看| 人妻 亚洲 视频| 人体艺术视频欧美日本| 视频区图区小说| 丰满迷人的少妇在线观看| 校园人妻丝袜中文字幕| 国产亚洲一区二区精品| 高清视频免费观看一区二区| 一区在线观看完整版| 日本爱情动作片www.在线观看| 国产欧美日韩一区二区三区在线| 久久鲁丝午夜福利片| 99热国产这里只有精品6| 国产成人免费无遮挡视频| 91成人精品电影| av.在线天堂| 日日摸夜夜添夜夜爱| 国产精品国产av在线观看| 各种免费的搞黄视频| 欧美 日韩 精品 国产| 日韩一本色道免费dvd| 高清欧美精品videossex| 国产精品.久久久| 丰满迷人的少妇在线观看| 日韩精品免费视频一区二区三区| 国产精品人妻久久久影院| 午夜精品国产一区二区电影| 如日韩欧美国产精品一区二区三区| 中文乱码字字幕精品一区二区三区| 免费高清在线观看日韩| 99久久综合免费| 亚洲男人天堂网一区| 夫妻午夜视频| 美女中出高潮动态图| 桃花免费在线播放| 男女午夜视频在线观看| 久久久久久伊人网av| 纯流量卡能插随身wifi吗| freevideosex欧美| 少妇人妻精品综合一区二区| 精品福利永久在线观看| 亚洲欧美中文字幕日韩二区| 久久精品国产亚洲av高清一级| 国精品久久久久久国模美| 久久亚洲国产成人精品v| 大片电影免费在线观看免费| 精品国产一区二区三区四区第35| 亚洲av中文av极速乱| 色婷婷久久久亚洲欧美| 国产成人精品久久久久久| 在线观看免费高清a一片| 欧美97在线视频| 五月伊人婷婷丁香| 少妇人妻精品综合一区二区| av.在线天堂| 色视频在线一区二区三区| 亚洲欧美色中文字幕在线| 如日韩欧美国产精品一区二区三区| 亚洲欧美清纯卡通| 亚洲久久久国产精品| 成年女人毛片免费观看观看9 | 亚洲精品在线美女| 久久久久久久久久久免费av| 欧美精品av麻豆av| 久久精品国产亚洲av天美| 免费女性裸体啪啪无遮挡网站| 国产激情久久老熟女| 熟妇人妻不卡中文字幕| 欧美日韩视频精品一区| 欧美成人精品欧美一级黄| 丝袜人妻中文字幕| 国产精品一二三区在线看| 老鸭窝网址在线观看| 性色avwww在线观看| 90打野战视频偷拍视频| 亚洲欧美一区二区三区久久| 久久久久久久久久人人人人人人| 精品卡一卡二卡四卡免费| 1024香蕉在线观看| 久久久精品国产亚洲av高清涩受| 一边摸一边做爽爽视频免费| 久久久国产精品麻豆| 多毛熟女@视频| 久久久精品区二区三区| 男女下面插进去视频免费观看| 国产1区2区3区精品| 人妻 亚洲 视频| 亚洲av男天堂| 精品99又大又爽又粗少妇毛片| 精品人妻偷拍中文字幕| 一本大道久久a久久精品| 亚洲在久久综合| 久久久国产一区二区| 啦啦啦啦在线视频资源| 永久免费av网站大全| 日日摸夜夜添夜夜爱| 另类亚洲欧美激情| 一本—道久久a久久精品蜜桃钙片| 丝袜在线中文字幕| 18在线观看网站| 日韩免费高清中文字幕av| 亚洲精品美女久久av网站| 国产激情久久老熟女| 伊人久久国产一区二区| 久久99蜜桃精品久久| 99久国产av精品国产电影| 在线天堂中文资源库| 国产精品 欧美亚洲| 老女人水多毛片| 国产 一区精品| 国产精品久久久久久精品电影小说| 国产一区二区激情短视频 | 精品视频人人做人人爽| 久久人妻熟女aⅴ| 亚洲精品国产色婷婷电影| 在线观看www视频免费| 永久网站在线| 成人影院久久| 免费看av在线观看网站| 色网站视频免费| 大码成人一级视频| 中文字幕精品免费在线观看视频| 欧美人与性动交α欧美软件| 久久久久精品久久久久真实原创| 国产欧美日韩综合在线一区二区| 免费观看性生交大片5| 韩国精品一区二区三区| 男女啪啪激烈高潮av片| 亚洲精品国产av成人精品| 乱人伦中国视频| 国产精品久久久久久av不卡| 2021少妇久久久久久久久久久| 蜜桃国产av成人99| 1024香蕉在线观看| 免费高清在线观看日韩| 啦啦啦在线观看免费高清www| 久久国内精品自在自线图片| 久久久久久久国产电影| 国产乱人偷精品视频| 亚洲成人av在线免费| 一本—道久久a久久精品蜜桃钙片| 下体分泌物呈黄色| 蜜桃在线观看..| 日韩一本色道免费dvd| 久久久久久伊人网av| 大陆偷拍与自拍| 一本色道久久久久久精品综合| 一级,二级,三级黄色视频| 日韩三级伦理在线观看| 亚洲欧美精品综合一区二区三区 | 丝袜脚勾引网站| 国产淫语在线视频| 日韩av在线免费看完整版不卡| 黑人欧美特级aaaaaa片| 在线亚洲精品国产二区图片欧美| 久久精品国产亚洲av涩爱| 一本—道久久a久久精品蜜桃钙片| 麻豆精品久久久久久蜜桃| 国产精品久久久久久久久免| 精品酒店卫生间| 大码成人一级视频| 欧美成人精品欧美一级黄| 国产福利在线免费观看视频| 一级毛片电影观看| 欧美成人精品欧美一级黄| 99国产精品免费福利视频| 桃花免费在线播放| 韩国高清视频一区二区三区| 免费女性裸体啪啪无遮挡网站| 人人澡人人妻人| 丝袜喷水一区| 欧美日韩成人在线一区二区| 国产1区2区3区精品| 国产精品久久久av美女十八| 国产免费现黄频在线看| 亚洲综合色网址| 国产一级毛片在线| 亚洲一码二码三码区别大吗| 久久99热这里只频精品6学生| 97精品久久久久久久久久精品| 久久鲁丝午夜福利片| 人人妻人人澡人人看| 最近的中文字幕免费完整| 欧美国产精品一级二级三级| 久久精品夜色国产| 秋霞伦理黄片| 午夜av观看不卡| 国产成人a∨麻豆精品| 久久久久久人人人人人| 久久精品熟女亚洲av麻豆精品| av在线观看视频网站免费| 国产综合精华液| 久久精品国产鲁丝片午夜精品| 满18在线观看网站| 久久久久久久久免费视频了| 天天躁夜夜躁狠狠躁躁| 国产一区二区激情短视频 | 一本大道久久a久久精品| 午夜福利一区二区在线看| 国产精品一国产av| 黄色视频在线播放观看不卡| 中文字幕另类日韩欧美亚洲嫩草| 99九九在线精品视频| 亚洲国产日韩一区二区| 母亲3免费完整高清在线观看 | 纯流量卡能插随身wifi吗| 亚洲欧美一区二区三区黑人 | 午夜福利网站1000一区二区三区| 老司机影院成人| 亚洲精品美女久久久久99蜜臀 | 多毛熟女@视频| 国产又爽黄色视频| 精品一区二区三区四区五区乱码 | 一区二区三区乱码不卡18| 欧美日本中文国产一区发布| 黄色怎么调成土黄色| 伦理电影大哥的女人| av国产久精品久网站免费入址| 黑丝袜美女国产一区| 久久韩国三级中文字幕| 寂寞人妻少妇视频99o| 亚洲精品日本国产第一区| 纵有疾风起免费观看全集完整版| 亚洲欧美精品综合一区二区三区 | 超碰成人久久| 边亲边吃奶的免费视频| 三级国产精品片| 美女主播在线视频| 女人高潮潮喷娇喘18禁视频| 十分钟在线观看高清视频www| 国产色婷婷99| 久久综合国产亚洲精品| 十分钟在线观看高清视频www| 欧美97在线视频| 国产片内射在线| 国产精品二区激情视频| 99久久综合免费| 高清视频免费观看一区二区| 男女国产视频网站| 涩涩av久久男人的天堂| 国产免费现黄频在线看| 青草久久国产| 丝袜脚勾引网站| 精品福利永久在线观看| 咕卡用的链子| 色94色欧美一区二区| 亚洲,一卡二卡三卡| 亚洲欧美成人精品一区二区| 大片电影免费在线观看免费| 日日啪夜夜爽| 搡女人真爽免费视频火全软件| 国产一级毛片在线| 18在线观看网站| 一本大道久久a久久精品| 人人妻人人爽人人添夜夜欢视频| 久久精品国产自在天天线| 人妻人人澡人人爽人人| 99re6热这里在线精品视频| 熟女少妇亚洲综合色aaa.| 久久精品国产鲁丝片午夜精品| 亚洲av国产av综合av卡| 亚洲欧美一区二区三区久久| 老汉色av国产亚洲站长工具| 2018国产大陆天天弄谢| 97在线人人人人妻| 日本vs欧美在线观看视频| 免费看av在线观看网站| 男女边摸边吃奶| 男人爽女人下面视频在线观看| 女性生殖器流出的白浆| 久久人人爽av亚洲精品天堂| 80岁老熟妇乱子伦牲交| 精品久久蜜臀av无| 日本91视频免费播放| 免费在线观看黄色视频的| 亚洲三区欧美一区| 只有这里有精品99| 97在线人人人人妻| 国产人伦9x9x在线观看 | 一区二区三区四区激情视频| 在线天堂最新版资源| 一级a爱视频在线免费观看| 26uuu在线亚洲综合色| 天天影视国产精品| 久久人人爽av亚洲精品天堂| 熟女电影av网| 亚洲精品成人av观看孕妇| 国产麻豆69| h视频一区二区三区| 久久97久久精品| 精品国产乱码久久久久久男人| 在线免费观看不下载黄p国产| 少妇人妻 视频| 成年动漫av网址| 免费观看无遮挡的男女| 亚洲av欧美aⅴ国产| 亚洲精品国产av蜜桃| 亚洲精品国产一区二区精华液| 91精品三级在线观看| 欧美 亚洲 国产 日韩一| 老鸭窝网址在线观看| 亚洲精品中文字幕在线视频| 国产亚洲欧美精品永久| 免费在线观看完整版高清| 亚洲精品国产av成人精品| 丝袜人妻中文字幕| 国产深夜福利视频在线观看| www.精华液| 日本黄色日本黄色录像| 美女国产高潮福利片在线看| 街头女战士在线观看网站| 欧美av亚洲av综合av国产av | 久久久精品区二区三区| 国产成人免费观看mmmm| 午夜免费男女啪啪视频观看| 大片免费播放器 马上看| 成人毛片a级毛片在线播放| 亚洲av综合色区一区| 久久久久久伊人网av| 中文字幕人妻丝袜一区二区 | 日韩一区二区视频免费看| 午夜免费鲁丝| 成人毛片60女人毛片免费| 国产精品久久久久久精品古装| 亚洲一级一片aⅴ在线观看| 久久久欧美国产精品| 亚洲婷婷狠狠爱综合网| 免费高清在线观看视频在线观看| 在线免费观看不下载黄p国产| 99热网站在线观看| 国产午夜精品一二区理论片| 日韩欧美精品免费久久| 亚洲精品自拍成人| 伊人久久大香线蕉亚洲五| 夫妻午夜视频| 大香蕉久久网| 一本大道久久a久久精品| 亚洲av日韩在线播放| 亚洲欧美日韩另类电影网站| 日本色播在线视频| 国产色婷婷99| 大片免费播放器 马上看| 电影成人av| 三上悠亚av全集在线观看| 97在线视频观看| 国产精品一区二区在线观看99| 人妻系列 视频| 又大又黄又爽视频免费| 国产成人午夜福利电影在线观看| 亚洲综合精品二区| 欧美bdsm另类| 国产又色又爽无遮挡免| 日韩欧美精品免费久久| 美女主播在线视频| 18禁国产床啪视频网站| 国语对白做爰xxxⅹ性视频网站|