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

    斑海豹線粒體基因組序列分析及比較研究

    2017-12-18 03:32:54高祥剛鹿志創(chuàng)田甲申韓家波赫崇波宋新然
    水產科學 2017年4期
    關鍵詞:斑海豹遼東灣阿拉斯加

    高祥剛,鹿志創(chuàng),田甲申,韓家波,赫崇波,宋新然

    ( 1.遼寧省海洋水產科學研究院,遼寧省海洋水產分子生物學重點實驗室,遼寧 大連 116023; 2. 大連圣亞旅游控股股份有限公司,遼寧 大連 116023 )

    斑海豹線粒體基因組序列分析及比較研究

    高祥剛1,鹿志創(chuàng)1,田甲申1,韓家波1,赫崇波1,宋新然2

    ( 1.遼寧省海洋水產科學研究院,遼寧省海洋水產分子生物學重點實驗室,遼寧 大連 116023; 2. 大連圣亞旅游控股股份有限公司,遼寧 大連 116023 )

    本研究采用LA-PCR擴增和引物步移相結合的方法對我國遼東灣斑海豹的線粒體基因組全序列進行了測定,并完成了基因組成和系統(tǒng)發(fā)生學分析。遼東灣斑海豹線粒體基因組長為16 754 bp,其基因構成與其他海洋哺乳動物基本一致,包括37個基因(2個rRNA基因、22個tRNA基因和13個蛋白質編碼基因)和1個控制區(qū)。在其37個基因中,ND6、tRNAAla、tRNAAsn、tRNACys、tRNATyr、tRNASer(UCN)、tRNAGlu、tRNAPro位于L鏈上,其余均位于H鏈上。對我國遼東灣、韓國和美國阿拉斯加斑海豹線粒體基因組進行了比較分析,結果顯示,三者線粒體基因組長度略有不同,差異主要體現(xiàn)在以重復片段形式存在的控制區(qū)序列中。值得注意的是,在編碼區(qū)中,3個目標群體的線粒體ND5基因核苷酸差異最為顯著,共檢測到了18個突變位點,造成4處氨基酸序列變異。斑海豹線粒體基因組蛋白編碼基因的變異分析及海豹科系統(tǒng)進化分析結果顯示,我國遼東灣和韓國斑海豹的序列同源性顯著高于美國阿拉斯加斑海豹,說明遼東灣和韓國斑海豹很可能來自同一個種群。

    斑海豹;線粒體基因組;結構比較;系統(tǒng)進化

    線粒體DNA(mtDNA)是動物細胞核外唯一的遺傳物質,與核DNA相比,具有母系遺傳、分子結構簡單、進化速度快、無組內特異性等特點。目前mtDNA已成為研究物種起源與系統(tǒng)發(fā)生、近緣物種和種內群體間遺傳分化、物種鑒定、遺傳多樣性分析等的有力工具[1]。

    海豹是對鰭足亞目種海豹科動物的統(tǒng)稱,現(xiàn)存的海豹種類較多,目前共有10屬19種。斑海豹(Phocalargha)屬食肉目、犬形亞目、海豹科、斑海豹屬,主要分布于北太平洋的北部和西部海域及其沿岸和島嶼。斑海豹在我國主要分布于渤海和黃海海域,偶見于南海。是唯一能在我國海域繁殖的鰭腳類動物,屬于國家二類重點保護野生動物,也是黃、渤海生態(tài)系統(tǒng)健康狀況的重要指示生物,因而受到國內外學者的普遍重視[2-3]。

    筆者通過LA-PCR擴增和引物步移相結合的方法對我國遼東灣斑海豹的線粒體基因組全序列進行了測定,初步分析了斑海豹線粒體基因組序列的基本特征,探討了我國、韓國、美國阿拉斯加斑海豹線粒體基因組基因排列、蛋白質編碼基因和多態(tài)位點等。同時利用線粒體基因組蛋白編碼序列,分析了海豹科的系統(tǒng)分類關系。本研究對于海豹的遺傳分化與系統(tǒng)發(fā)生研究、物種鑒定、遺傳多樣性的保護等諸多方面提供寶貴的遺傳信息。

    1 材料與方法

    1.1 試驗材料和基因組DNA提取

    解剖盤錦雙臺河口救助的當年死亡斑海豹幼崽的肝臟,清洗后立即保存于-20 ℃冰柜中。采用常規(guī)酚氯仿法提取全基因組DNA[4],溶于TE緩沖液中,4 ℃保存?zhèn)溆谩?/p>

    1.2 PCR擴增及序列測定

    參考韓國、美國阿拉斯加斑海豹(GenBank號分別為:FJ895151、NC_008430)線粒體基因組序列,應用FastPCR v3.6軟件設計引物COⅠ_F 和COⅠ_R。利用課題組已設計引物[5-6]COⅠ_F、COⅠ_R、L16274_F 、H34_R,形成引物組合(表1)L16274_F + COⅠ_R 、COⅠ_F + ND4L_R、 ND4L_F + H34_R進行遼東灣斑海豹的線粒體基因組Long-PCR擴增。反應條件和反應體系見表2。測序采用引物步移法進行雙向測通。因為PCR產物直接測序法不能準確讀出測序引物后面約50 bp長度的序列,所以在設計擴增片段的引物時,盡可能保證相鄰擴增片段之間包含80 bp以上的序列重疊,以保證測序后序列的準確性及正確拼接。由于在控制區(qū)測序時出現(xiàn)套峰現(xiàn)象,此段區(qū)域采用克隆后測序的方法。

    表1 用于擴增遼東灣斑海豹線粒體基因組的引物

    表2 遼東灣斑海豹線粒體基因組擴增的反應條件和反應體系

    1.3 序列分析、基因預測

    首先使用Chromas 2.4.1軟件打開和核實基因測序得到的峰圖,并用CONSED手工檢查以避免錯配[7]。利用DOGMA軟件對13個蛋白質編碼基因和2個核糖體RNA 的位置進行預測[8],與其他海豹屬物種的線粒體基因組比對后確定精確的起始位置。利用tRNAscan-SE 1.21軟件對轉運RNA及它們的二級結構特征進行預測[9]。

    1.4 斑海豹線粒體基因組的比較分析

    將所測得的遼東灣斑海豹mtDNA序列,應用Clustal 1.83和Mega 5.0生物信息學軟件,分別與韓國、美國阿拉斯加斑海豹線粒體基因組全序列進行比較分析,探討三者在蛋白編碼基因的氨基酸和核苷酸序列、tRNAs、rRNAs上的差異,以及在非編碼區(qū)的變異情況。

    1.5 海豹科線粒體基因組的系統(tǒng)進化

    對NCBI上已公布的海豹科的15個種、17條線粒體基因組DNA序列的13個蛋白編碼基因的氨基酸序列和核苷酸序列采用鄰接法分別構建系統(tǒng)關系樹。

    2 結 果

    2.1 遼東灣斑海豹線粒體基因組的組成

    將所測得的遼東灣斑海豹線粒體基因組全序列拼接完整后提交至GenBank(NCBI登錄號:KT818831),其長度為16 754 bp,呈現(xiàn)典型的海洋哺乳動物線粒體基因組的特點。即斑海豹線粒體基因組為環(huán)狀,由37個基因組成,其中的 13個蛋白質編碼基因包括:1個細胞色素b基因(CytB)、2個ATP酶亞基的基因(ATPase6, ATPase8)、3個細胞色素氧化酶亞基的基因(COⅠ,COⅡ,CO Ⅲ)和7個呼吸鏈NADH脫氫酶亞基的基因(ND1-6,ND4L);2個rRNA基因分別為1個12S rRNA和 1個16S rRNA;另外還有22個tRNA基因。

    在13個蛋白質編碼基因中,12個是由H 鏈編碼,僅ND6由L鏈編碼;在22個tRNA中,tRNAAla、tRNAAsn、tRNACys、tRNATyr、tRNASer(UCN)、tRNAGlu、tRNAPro的編碼基因位于L鏈上,其余則位于H 鏈上。

    2.2 與其他斑海豹的比較

    將所測得的中國遼東灣斑海豹線粒體基因組與韓國、美國阿拉斯加斑海豹線粒體基因組進行比較分析發(fā)現(xiàn),三者線粒體基因組長度略有不同,分別為16 754、16 728、16 701 bp,主要的變化出現(xiàn)在控制區(qū)(OH),OH長度分別為1316、1290、1264 bp。

    對上述3個線粒體基因組的37個基因分別進行了比較研究,在rRNA中,發(fā)現(xiàn)5個單核苷酸變異;在tRNA中,發(fā)現(xiàn)阿拉斯加斑海豹的tRNALeu出現(xiàn)1個堿基缺失,其他的tRNAPhe、tRNAAla、tRNATyr、tRNAThr等各有1個單核苷酸變異,除此之外的17個tRNA完全一致;13個蛋白質編碼基因中,均發(fā)現(xiàn)有核苷酸變異,其中8個基因(ATPase6、COⅠ、COⅡ、CytB、ND1、ND3、ND4、 ND4L)相對比較保守,其氨基酸序列在3個基因組中完全一致,5個基因出現(xiàn)氨基酸序列變化(表3),其中,ND5的核苷酸突變程度最大,總長1821 bp的核苷酸中有18個突變位點,變異率為0.99%,出現(xiàn)4個氨基酸序列突變。

    表3 斑海豹3個線粒體基因組蛋白編碼基因的突變比較分析

    2.3 海豹科線粒體基因組的系統(tǒng)進化

    以海洋哺乳動物的鯨目白鰭豚(Lipotesvexillifer,登錄號:NC_007629)、和海牛目儒艮(Dugongdugon,登錄號: NC_003314)作外群,根據(jù)鰭足亞目海豹科的2個亞科17條mtDNA序列,用鄰接法構建13個編碼蛋白的氨基酸序列系統(tǒng)關系樹(圖1),與最大似然法和最小進化法構建的系統(tǒng)樹基本一致。結果顯示遼東灣斑海豹首先和韓國斑海豹聚為一支,然后與阿拉斯加斑海豹聚為一支,3個斑海豹再與海豹亞科的其他海豹聚為一大支。僧海豹亞科的5個物種序列單獨聚為一大支。根據(jù)13個編碼蛋白的核苷酸與氨基酸序列構建的系統(tǒng)關系樹結果基本一致,并且均顯示遼東灣斑海豹與韓國斑海豹的關系較阿拉斯加斑海豹更近。

    圖1 基于13種線粒體蛋白質的氨基酸序列的鄰接法進化樹中國遼東灣斑海豹(KT818831)、韓國斑海豹(FJ895151)、美國阿拉斯加斑海豹(NC_008430)、港海豹(NC_001325)、環(huán)斑海豹(NC_008433)、貝加爾環(huán)海豹(NC_008432)、里海環(huán)斑海豹(NC_008431)、灰海豹(NC_001602)、鞍紋海豹(NC_008429)、帶紋海豹(NC_008428)、冠海豹(NC_008427)、髯海豹(NC_008426)、夏威夷僧海豹(NC_008421)、食蟹海豹(NC_008423)、韋德爾海豹(NC_008424)、豹形海豹(NC_008422)、南象海豹(NC_008422)、白鰭豚(NC_007629)和儒艮(NC_003314).

    3 討 論

    3.1 斑海豹線粒體控制區(qū)異質性現(xiàn)象

    異質性是指在有些脊椎動物,如中華鱘(Acipensersinensis)、弓鰭魚(Amiacalva),體內存在多種重復序列數(shù)目不同的線粒體基因組。至于控制區(qū)異質型和多態(tài)性水平高的原因,有學者指出線粒體基因異質性中的重復序列可能形成發(fā)卡結構,發(fā)卡結構引起高頻率回復突變,這也許是導致異質性形成的原因,當然也不排除父本mtDNA的滲漏造成異質性的可能性。同時異質性也是物種進化的原因,線粒體中的串聯(lián)重復、存在的大量插入與缺失等是物種進化的特征[10-12]。

    而本文通過比較分析我國遼東灣、韓國、美國阿拉斯加斑海豹線粒體基因組控制區(qū)發(fā)現(xiàn),三者主要變化出現(xiàn)在控制區(qū)CSB-Ⅲ上游約500 bp序列中,其上有大量重復片段,且重復片段長度有差異,這些重復片段主要為:ACGTAC、GCACAC、ACGT、ACGC、AC、GC。Bahn等[13-14]在比較韓國與日本斑海豹的控制區(qū)57個單倍型時,未發(fā)現(xiàn)完全一致的序列。孫凡越[15]也對7頭遼東灣斑海豹控制區(qū)序列進行了克隆測序,得到了146個單倍型,其中異質型數(shù)目占所測序列總數(shù)的92%,并發(fā)現(xiàn)序列差異的原因主要是存在6種重復單元(GTACACAC、GTACAC、GTAC、GCACACAC、GCACAC和GCAC),而各單倍型的重復單元及重復次數(shù)均不相等,由此造成了遼東灣斑海豹不同個體之間出現(xiàn)長度異質性。這一結論與本試驗的研究結果較類似,說明我國遼東灣斑海豹的線粒體DNA異質性程度比較高,因此,OH重復單元的特征以及線粒體ND5基因的序列多態(tài)性,都有可能成為區(qū)分斑海豹個體和不同種群的重要分子標記。

    3.2 我國遼東灣斑海豹與韓國、阿拉斯加斑海豹的親緣關系

    本文在比較分析斑海豹線粒體基因組蛋白編碼基因的突變時,發(fā)現(xiàn)在5個突變蛋白質編碼基因中,共有10個位點出現(xiàn)氨基酸變化,總突變核苷酸位點數(shù)為61個,由表3可見,遼東灣和韓國斑海豹的蛋白編碼基因中,只出現(xiàn)1個氨基酸突變,證明這兩個個體的序列比較保守;相對應的阿拉斯加與遼東灣、阿拉斯加與韓國的突變氨基酸位點數(shù)分別為6個、10個,突變核苷酸位點數(shù)分別為56個、58個。通過對斑海豹3個線粒體基因組蛋白編碼基因的突變比較分析(表3)及系統(tǒng)進化分析(圖1)可以看出,遼東灣和韓國斑海豹的序列同源性顯著高于阿拉斯加與遼東灣、阿拉斯加與韓國,說明遼東灣和韓國斑海豹很可能來自同一個斑海豹種群。

    北太平洋的北部和西部海域及其沿岸和島嶼都是斑海豹在世界上的主要分布區(qū)[16],美國阿拉斯加海域與中國、韓國沿海相距較遠,而從韓家波等[17]的衛(wèi)星信標跟蹤斑海豹研究可以看出,我國遼東灣的斑海豹可以到達韓國西海岸的白翎島,說明我國海域與韓國西海岸的斑海豹有交流,這進一步表明了我國遼東灣與韓國斑海豹較美國阿拉斯加斑海豹的親緣關系更近。

    [1] 邢晶晶. 分子遺傳標記及其技術在水產生物中的研究與應用[J]. 水產學雜志, 2002, 14(2):61-70.

    [2] Shaughnessy P D, Fay R H. A review of the taxonomy and nomenclature of North Pacific harbour seals[J]. J Zool, 1997,182(3):385-419.

    [3] Won C, Yoo B H. Abundance, seasonal haul-out patterns and conservation of spotted sealPhocalarghaalong the coast of Bak-ryoung Island[J]. Oryx,2004,38(1):109-112.

    [4] Sambrook J,F(xiàn)ritsch E F,Maniatis T. Molecular Cloning[M]. New York:Cold Spring Harbor Laboratory Press,1989.

    [5] 韓家波,赫崇波,王強,等. 遼東灣斑海豹線粒體ND4、tRNA Arg、ND4L和ND3基因序列分析[J]. 水產科學, 2006,25(10):500-504.

    [6] 韓家波,赫崇波,王效敏,等. 遼東灣斑海豹線粒體蘇氨酸和脯氨酸及控制區(qū)部分序列分析[J]. 水產科學,2007,26(2):74-78.

    [7] Gordon D, Abajian C, Green P. Consed: a graphical tool for sequence finishing [J]. Genome Res, 1998, 8(3):195-202.

    [8] Wyman S K, Jansen R K, Boore J L. Automatic annotation of organellar genomes with DOGMA[J]. Bioinformatics, 2004, 20(17):3252-3255.

    [9] Lowe T M, Eddy S R. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence [J]. Nucleic Acids Res, 1997, 25(5): 955-964.

    [10] 蘇瑛. 動物mtDNA控制區(qū)及保守與異質[J]. 四川動物, 2005, 24(4):669-672.

    [11] 王巍, 朱華, 胡紅霞,等. 五種鱘魚線粒體控制區(qū)異質性和系統(tǒng)發(fā)育分析[J]. 動物學研究, 2009, 30(5): 487-496.

    [12] 肖海霞, 帕熱哈提江·吾甫爾, 玉山江,等. 線粒體DNA在動物起源進化和遺傳多樣性中的應用[J]. 獸醫(yī)導刊, 2015(12):39-40.

    [13] Bahn S Y, Choi E H, Jang K H, et al. Complete mitochondrial genome of the Korean spotted sealPhocalargha(Mammalia, Pinnipedia, Phocidae): genetic differences betweenP.larghaandPhocavitulina[J]. Mitochondrial DNA, 2009, 20(5/6):105-106.

    [14] Mizuno A W, Onuma M, Takahashi M, et al. Population genetic structure of the spotted sealPhocalarghaalong the Coast of Hokkaido, based on mitochondrial DNA sequences[J]. Zool Sci, 2003, 20(6): 783-788.

    [15] 孫凡越. 遼東灣斑海豹種群遺傳特征的研究 [D]. 大連:遼寧師范大學,2009.

    [16] Rugh D J, Shelden K E W, Withrow D E. Spotted seals,Phocalargha, in Alaska[J]. Marine Fisheries Review, 1997, 59(1):1-18.

    [17] 韓家波,鹿志創(chuàng),田甲申,等. 基于衛(wèi)星信標跟蹤的斑海豹放流效果研究[J]. 獸類學報, 2013, 33(4): 300-307.

    AnalysisandComparisonofCompleteMitochondrialGenomeinSpottedSealPhocalargha

    GAO Xianggang1, LU Zhichuang1, TIAN Jiashen1, HAN Jiabo1, HE Chongbo1, SONG Xinran2

    ( 1. Key Lab of Marine Fishery Molecular Biology of Liaoning Province, Liaoning Ocean and Fisheries Science Research Institute, Dalian 116023, China; 2. Dalian Sun Asia Tourism Holding Co, Ltd, Dalian 116023, China )

    The complete sequence of the mitochondrial genome of the spotted seals (Phocalargha) in the Liaodong Gulf has been determined by long and accurate polymerase chain reaction (LA-PCR) and primer walking technology. The genome is 16 754 bp in size, containing 2 ribosomal RNA genes, 13 protein-coding genes, 22 transfer RNA (tRNA) genes and one control region, demonstrating a structure basically in accordance with other marine mammals. The ND6, tRNAAla, tRNAAsn, tRNACys, tRNATyr, tRNASer(UCN), tRNAGluand tRNAProgenes were located on the light strand (L strand), and the remainders lied on the heavy strand (H strand). The comparative analysis of the seal mitochondrial genome among Liaodong, Korean and Alaskan populations revealed the existence of length differences and that the major variations exhibited in the form of repeat motifs were located in the control region. Notably, in the coding region of the three targeted populations, ND5 showed the maximal mutation degree in total of 18 mutant bases was detected, which resulted in four amino acid changes. The variant and phylogenetic analysis of protein-coding genes indicated that the Liaodong seal shared higher homology with Korean seal than with Alaskan seal. It is inferred that Liaodong seal and Korean seal probably belong to one population.

    Phocalargha; mitochondrial genome; comparative structure;phylogeny

    10.16378/j.cnki.1003-1111.2017.04.018

    Q786

    A

    1003-1111(2017)04-0509-05

    2016-07-01;

    2016-08-16.

    遼寧省海洋與漁業(yè)科研項目(201417);香港海洋公園保育基金資助項目(MM04-1516).

    高祥剛(1980-),男,副研究員;研究方向:動物分子遺傳學. E-mail:xiangganggao@163.com. 通訊作者: 赫崇波(1961-),男,研究員,博士;研究方向:動物分子遺傳學. E-mail:chongbohe@163.com.

    猜你喜歡
    斑海豹遼東灣阿拉斯加
    2023年4月3日。在阿拉斯加聯(lián)合太平洋多國戰(zhàn)備中心舉行的阿拉斯加23-02演習期間,美軍士兵正扮演敵方調查周圍地區(qū)
    軍事文摘(2023年9期)2023-04-30 04:11:52
    2021年6月18日,一架F-16戰(zhàn)隼在聯(lián)合太平洋阿拉斯加靶場上空
    軍事文摘(2021年15期)2021-08-11 13:01:12
    大連百頭斑海豹幼崽被盜獵,沒有買賣才沒有傷害
    渤海精靈斑海豹
    天天愛科學(2017年2期)2017-04-29 11:25:40
    天然氣分輸站規(guī)劃選址可行性論證初探——以盤錦遼東灣新區(qū)榮興分輸站為例
    遼東灣西部海域表層沉積物重金屬的含量分布與污染評價
    遼東灣坳陷南部旅大10-2構造區(qū)源蓋動態(tài)耦合關系定量評價
    遼東灣地區(qū)濁積扇體識別及儲層地震描述方法研究
    斷塊油氣田(2012年6期)2012-03-25 09:53:57
    冰川的藍,極光的燦——阿拉斯加夢飛行
    50石油職工志愿者5年守護斑海豹
    亚洲美女搞黄在线观看| АⅤ资源中文在线天堂| 免费观看的影片在线观看| 成人二区视频| av在线亚洲专区| 久久精品人妻少妇| or卡值多少钱| 国产老妇伦熟女老妇高清| 国产av一区在线观看免费| 欧美+日韩+精品| 日韩一本色道免费dvd| 亚洲18禁久久av| 国模一区二区三区四区视频| 久久久久久久国产电影| www.av在线官网国产| 人人妻人人澡欧美一区二区| 免费观看人在逋| 国产中年淑女户外野战色| 十八禁国产超污无遮挡网站| 亚洲第一区二区三区不卡| 久久久久久伊人网av| 色播亚洲综合网| 国产精品伦人一区二区| 少妇熟女aⅴ在线视频| 村上凉子中文字幕在线| 我要搜黄色片| 久久草成人影院| 国产伦理片在线播放av一区| 日本色播在线视频| 观看美女的网站| 国产乱来视频区| 亚洲人成网站在线观看播放| 人人妻人人澡欧美一区二区| 伦精品一区二区三区| 国产免费一级a男人的天堂| 国产精品av视频在线免费观看| 日本免费a在线| 国产精品国产三级国产专区5o | 久久99热这里只有精品18| 嫩草影院新地址| 麻豆乱淫一区二区| 国产伦在线观看视频一区| 男人舔奶头视频| 欧美性猛交黑人性爽| 日韩 亚洲 欧美在线| 免费观看精品视频网站| av.在线天堂| 亚洲av.av天堂| 亚洲精品影视一区二区三区av| 精品久久久久久成人av| 熟女电影av网| av免费在线看不卡| 成人性生交大片免费视频hd| 国产成人午夜福利电影在线观看| 久久精品人妻少妇| 嫩草影院新地址| 超碰97精品在线观看| 国产午夜精品论理片| 国产精品国产三级专区第一集| 国国产精品蜜臀av免费| 国产单亲对白刺激| 免费无遮挡裸体视频| 久久精品夜夜夜夜夜久久蜜豆| 免费观看性生交大片5| 国产精品久久久久久久久免| 日韩,欧美,国产一区二区三区 | 久久久a久久爽久久v久久| 一级黄色大片毛片| 久久人人爽人人爽人人片va| av女优亚洲男人天堂| 久久精品熟女亚洲av麻豆精品 | 91aial.com中文字幕在线观看| 国产精品精品国产色婷婷| 精品无人区乱码1区二区| 国产又色又爽无遮挡免| 97人妻精品一区二区三区麻豆| 综合色丁香网| 亚洲欧美精品综合久久99| 日韩,欧美,国产一区二区三区 | 啦啦啦韩国在线观看视频| 我的老师免费观看完整版| 亚洲国产欧洲综合997久久,| 亚洲一区高清亚洲精品| av福利片在线观看| 高清av免费在线| 久99久视频精品免费| 中文字幕av在线有码专区| 国产三级中文精品| 日日摸夜夜添夜夜爱| 亚洲人成网站在线播| 国语自产精品视频在线第100页| 色综合亚洲欧美另类图片| 国模一区二区三区四区视频| 亚洲av中文字字幕乱码综合| 色噜噜av男人的天堂激情| 国产成人a∨麻豆精品| 99久久人妻综合| 美女脱内裤让男人舔精品视频| av免费在线看不卡| 免费在线观看成人毛片| 欧美激情国产日韩精品一区| 国产伦精品一区二区三区视频9| 亚洲va在线va天堂va国产| 丝袜喷水一区| 国产极品天堂在线| 六月丁香七月| 成人综合一区亚洲| 国产三级中文精品| 18禁在线无遮挡免费观看视频| 又爽又黄a免费视频| 在线天堂最新版资源| 欧美精品国产亚洲| 亚洲人与动物交配视频| 欧美不卡视频在线免费观看| 天堂√8在线中文| 亚洲成人久久爱视频| 日韩欧美在线乱码| 精品人妻视频免费看| 久久6这里有精品| 亚洲图色成人| 天天躁夜夜躁狠狠久久av| 噜噜噜噜噜久久久久久91| 久久99热6这里只有精品| www.av在线官网国产| 亚洲av男天堂| 好男人在线观看高清免费视频| 美女被艹到高潮喷水动态| 久久久午夜欧美精品| 日韩成人av中文字幕在线观看| 中文字幕制服av| 三级男女做爰猛烈吃奶摸视频| 亚洲精品乱久久久久久| 免费看av在线观看网站| 免费观看在线日韩| 搡老妇女老女人老熟妇| 成人综合一区亚洲| 色噜噜av男人的天堂激情| 午夜免费激情av| 国产色婷婷99| www日本黄色视频网| 欧美成人一区二区免费高清观看| 国产又色又爽无遮挡免| 国产成人免费观看mmmm| 国产精品久久久久久精品电影小说 | 男人舔女人下体高潮全视频| 一区二区三区乱码不卡18| 国产淫语在线视频| 亚洲美女视频黄频| 久久精品夜色国产| 国产精品1区2区在线观看.| 黑人高潮一二区| 欧美xxxx性猛交bbbb| 国产精品99久久久久久久久| 国内揄拍国产精品人妻在线| 2021少妇久久久久久久久久久| 免费一级毛片在线播放高清视频| 国产午夜精品久久久久久一区二区三区| 亚洲欧美日韩卡通动漫| 赤兔流量卡办理| 精华霜和精华液先用哪个| 午夜福利在线观看吧| 欧美97在线视频| 人人妻人人澡欧美一区二区| 国产亚洲精品久久久com| 午夜视频国产福利| 久久婷婷人人爽人人干人人爱| 亚洲最大成人手机在线| 亚洲成av人片在线播放无| 最新中文字幕久久久久| 熟女人妻精品中文字幕| 精品酒店卫生间| 亚洲18禁久久av| 一本一本综合久久| 在线播放无遮挡| 搡女人真爽免费视频火全软件| 久久国内精品自在自线图片| 伊人久久精品亚洲午夜| 亚洲精品乱久久久久久| 免费看av在线观看网站| 91久久精品电影网| 日韩av在线大香蕉| 中文字幕av在线有码专区| 久久久久精品久久久久真实原创| 国产成人aa在线观看| 免费看光身美女| 亚洲av免费在线观看| 在线免费观看不下载黄p国产| 九草在线视频观看| 精品久久国产蜜桃| 高清毛片免费看| 国产女主播在线喷水免费视频网站 | 亚洲国产精品成人综合色| 男女那种视频在线观看| 精华霜和精华液先用哪个| 亚洲最大成人手机在线| 国产黄片美女视频| 偷拍熟女少妇极品色| 高清在线视频一区二区三区 | 午夜免费激情av| 中国美白少妇内射xxxbb| 菩萨蛮人人尽说江南好唐韦庄 | 精品不卡国产一区二区三区| 天美传媒精品一区二区| 看片在线看免费视频| 国产一区二区在线观看日韩| 一级黄片播放器| 免费无遮挡裸体视频| 亚洲国产欧美人成| 直男gayav资源| 午夜日本视频在线| 99久国产av精品国产电影| 搞女人的毛片| 亚洲最大成人av| 内地一区二区视频在线| 欧美又色又爽又黄视频| 久久久国产成人免费| 内射极品少妇av片p| 国产亚洲91精品色在线| 99久久人妻综合| 亚洲精品久久久久久婷婷小说 | 伦理电影大哥的女人| 国产成人福利小说| 真实男女啪啪啪动态图| 男人舔女人下体高潮全视频| 亚洲一区高清亚洲精品| 日韩av在线免费看完整版不卡| 亚洲婷婷狠狠爱综合网| 国产成人一区二区在线| 有码 亚洲区| 国产精品永久免费网站| 91av网一区二区| 男人和女人高潮做爰伦理| 两性午夜刺激爽爽歪歪视频在线观看| 国产一区亚洲一区在线观看| 高清毛片免费看| 久久欧美精品欧美久久欧美| 亚洲四区av| 日韩国内少妇激情av| 国产成人福利小说| 亚洲在线观看片| 免费观看a级毛片全部| 男人和女人高潮做爰伦理| www.av在线官网国产| 色综合站精品国产| 丝袜喷水一区| 久久精品夜夜夜夜夜久久蜜豆| 在线免费观看不下载黄p国产| 国产av码专区亚洲av| 久久综合国产亚洲精品| 精品久久久久久久人妻蜜臀av| 婷婷六月久久综合丁香| 精品无人区乱码1区二区| 日本免费在线观看一区| 一区二区三区免费毛片| 免费无遮挡裸体视频| 乱人视频在线观看| 秋霞伦理黄片| 国产成年人精品一区二区| 国产精品伦人一区二区| 人人妻人人澡欧美一区二区| 在线免费观看不下载黄p国产| 日本免费在线观看一区| 欧美成人一区二区免费高清观看| 99久久九九国产精品国产免费| 九九热线精品视视频播放| 国产又黄又爽又无遮挡在线| 男女下面进入的视频免费午夜| 久久人人爽人人片av| 在线天堂最新版资源| 国产黄片视频在线免费观看| 亚洲天堂国产精品一区在线| 边亲边吃奶的免费视频| 午夜精品在线福利| 日韩在线高清观看一区二区三区| 亚洲成人中文字幕在线播放| 国产中年淑女户外野战色| 欧美三级亚洲精品| 18禁裸乳无遮挡免费网站照片| 免费看光身美女| 日韩强制内射视频| 欧美潮喷喷水| 免费电影在线观看免费观看| 国产精品久久电影中文字幕| 国产午夜精品久久久久久一区二区三区| 成人av在线播放网站| 精品不卡国产一区二区三区| 亚洲高清免费不卡视频| 熟妇人妻久久中文字幕3abv| 有码 亚洲区| 男女那种视频在线观看| 日韩大片免费观看网站 | 欧美日本亚洲视频在线播放| 三级毛片av免费| 亚洲电影在线观看av| 欧美成人a在线观看| 国产伦一二天堂av在线观看| 久久久久久久午夜电影| 麻豆一二三区av精品| av免费在线看不卡| 老师上课跳d突然被开到最大视频| 久久精品国产自在天天线| 国产黄色视频一区二区在线观看 | 麻豆av噜噜一区二区三区| 日本免费一区二区三区高清不卡| 99在线视频只有这里精品首页| 夫妻性生交免费视频一级片| 91久久精品国产一区二区成人| 中文资源天堂在线| 听说在线观看完整版免费高清| 国产探花在线观看一区二区| 久久人人爽人人爽人人片va| 免费观看精品视频网站| 国产精品国产三级国产av玫瑰| 亚洲美女搞黄在线观看| 最新中文字幕久久久久| 国产精品久久电影中文字幕| 亚洲av成人精品一二三区| 干丝袜人妻中文字幕| 亚洲精品自拍成人| 久久精品国产亚洲网站| 久久久午夜欧美精品| 精品人妻熟女av久视频| 91午夜精品亚洲一区二区三区| 毛片女人毛片| 精品一区二区三区人妻视频| 99久国产av精品国产电影| 九九在线视频观看精品| 男插女下体视频免费在线播放| 欧美不卡视频在线免费观看| 中国国产av一级| 又爽又黄a免费视频| 久久久欧美国产精品| 直男gayav资源| 国语自产精品视频在线第100页| or卡值多少钱| 国产色婷婷99| 久久这里只有精品中国| 麻豆成人午夜福利视频| videossex国产| 国产一区有黄有色的免费视频 | 九九在线视频观看精品| 国产av不卡久久| 99热全是精品| 国产黄片美女视频| 午夜福利在线观看吧| 一个人看的www免费观看视频| 欧美zozozo另类| 欧美激情国产日韩精品一区| 嫩草影院精品99| 成人午夜精彩视频在线观看| 久久久久久久午夜电影| 91精品一卡2卡3卡4卡| 九九热线精品视视频播放| 成人午夜精彩视频在线观看| 舔av片在线| 午夜免费男女啪啪视频观看| 噜噜噜噜噜久久久久久91| 久久久欧美国产精品| 欧美性猛交黑人性爽| 精品久久久久久久人妻蜜臀av| 亚洲18禁久久av| 日本一本二区三区精品| 国产精品一区二区三区四区久久| 中文字幕久久专区| 丰满乱子伦码专区| 中国美白少妇内射xxxbb| 听说在线观看完整版免费高清| 久久亚洲国产成人精品v| 青春草亚洲视频在线观看| 人人妻人人看人人澡| 麻豆国产97在线/欧美| 国产探花在线观看一区二区| 变态另类丝袜制服| 国产在线男女| 乱人视频在线观看| 日本熟妇午夜| 精品久久久久久久末码| 成年女人看的毛片在线观看| 麻豆成人午夜福利视频| 神马国产精品三级电影在线观看| 久久精品久久久久久久性| 日本与韩国留学比较| 成人鲁丝片一二三区免费| 亚洲成人精品中文字幕电影| 欧美成人精品欧美一级黄| 亚洲国产精品sss在线观看| 超碰97精品在线观看| 日本午夜av视频| 亚洲人成网站在线播| 99久久精品热视频| 亚洲在线观看片| 国产黄a三级三级三级人| 亚洲丝袜综合中文字幕| 91久久精品国产一区二区三区| kizo精华| 免费观看性生交大片5| 精品酒店卫生间| 好男人视频免费观看在线| 国产av码专区亚洲av| 2021少妇久久久久久久久久久| 三级国产精品欧美在线观看| 免费黄色在线免费观看| 91av网一区二区| 中文字幕免费在线视频6| 国内揄拍国产精品人妻在线| 久久久a久久爽久久v久久| 亚洲av成人av| 在线播放无遮挡| 少妇裸体淫交视频免费看高清| 熟妇人妻久久中文字幕3abv| 国产精品久久久久久av不卡| 午夜福利在线观看免费完整高清在| 国产精品熟女久久久久浪| av在线蜜桃| 最近2019中文字幕mv第一页| 亚洲精品,欧美精品| 亚洲欧美精品综合久久99| 亚洲av成人精品一二三区| 国产亚洲av片在线观看秒播厂 | 中国国产av一级| 亚洲在久久综合| 成人午夜高清在线视频| 一级黄色大片毛片| 国产精品久久久久久精品电影小说 | 国语对白做爰xxxⅹ性视频网站| 天堂中文最新版在线下载 | 噜噜噜噜噜久久久久久91| 亚洲性久久影院| 一个人免费在线观看电影| 国产毛片a区久久久久| 毛片一级片免费看久久久久| 日本黄大片高清| 国产美女午夜福利| 亚洲欧美精品专区久久| av线在线观看网站| 高清av免费在线| 日韩在线高清观看一区二区三区| 欧美zozozo另类| 国产av码专区亚洲av| 亚洲国产精品成人久久小说| 国产三级在线视频| 亚洲不卡免费看| 黑人高潮一二区| 欧美变态另类bdsm刘玥| 国产高清有码在线观看视频| 一级黄色大片毛片| 亚洲成人久久爱视频| 亚洲乱码一区二区免费版| 晚上一个人看的免费电影| 久久久久久久久大av| 国产在视频线在精品| 国产高清视频在线观看网站| 国产成人午夜福利电影在线观看| 日日摸夜夜添夜夜爱| 夜夜看夜夜爽夜夜摸| 国产淫片久久久久久久久| 国产亚洲5aaaaa淫片| 嘟嘟电影网在线观看| 热99在线观看视频| 久久精品熟女亚洲av麻豆精品 | 久久久久九九精品影院| 欧美一区二区精品小视频在线| 韩国高清视频一区二区三区| 亚洲自偷自拍三级| 啦啦啦韩国在线观看视频| 18禁在线播放成人免费| 老司机影院毛片| 麻豆成人午夜福利视频| 国产 一区精品| 国产精品永久免费网站| 亚洲欧洲日产国产| 九草在线视频观看| 亚洲av免费高清在线观看| 午夜福利在线在线| 九九爱精品视频在线观看| 美女国产视频在线观看| 成人欧美大片| 亚洲av熟女| 熟女电影av网| 日韩一区二区三区影片| 久久久成人免费电影| 欧美97在线视频| 免费看美女性在线毛片视频| 免费av不卡在线播放| 午夜爱爱视频在线播放| 亚洲成人中文字幕在线播放| 亚洲av成人精品一区久久| 97在线视频观看| 亚洲欧美日韩东京热| 春色校园在线视频观看| 人人妻人人澡欧美一区二区| av黄色大香蕉| 久久久久久久午夜电影| 免费无遮挡裸体视频| 国产老妇伦熟女老妇高清| 黄片无遮挡物在线观看| 久久人妻av系列| 中文字幕熟女人妻在线| 亚洲成人av在线免费| 简卡轻食公司| 黄片wwwwww| 夜夜爽夜夜爽视频| 综合色丁香网| 国产极品精品免费视频能看的| 国产精品国产三级国产专区5o | 小蜜桃在线观看免费完整版高清| 卡戴珊不雅视频在线播放| 午夜福利网站1000一区二区三区| 久热久热在线精品观看| 日韩大片免费观看网站 | 麻豆国产97在线/欧美| 少妇熟女aⅴ在线视频| 午夜精品在线福利| 美女大奶头视频| 最近最新中文字幕大全电影3| 国产高清国产精品国产三级 | 国产精品一区www在线观看| 97人妻精品一区二区三区麻豆| 国产高清三级在线| 少妇裸体淫交视频免费看高清| 在线天堂最新版资源| 建设人人有责人人尽责人人享有的 | 一级毛片我不卡| 日韩精品青青久久久久久| 九九爱精品视频在线观看| 色哟哟·www| 国产精品.久久久| 水蜜桃什么品种好| 婷婷六月久久综合丁香| 日韩欧美在线乱码| 18禁裸乳无遮挡免费网站照片| 国产精品国产三级专区第一集| .国产精品久久| 国产成人aa在线观看| 精品国产三级普通话版| 国产欧美日韩精品一区二区| 欧美激情在线99| 国产极品精品免费视频能看的| 男人狂女人下面高潮的视频| 亚洲av免费在线观看| 久久精品国产亚洲av涩爱| 亚洲va在线va天堂va国产| 久久99蜜桃精品久久| 国产一区有黄有色的免费视频 | 毛片女人毛片| 日本午夜av视频| 国产淫语在线视频| 毛片一级片免费看久久久久| 永久免费av网站大全| 国产v大片淫在线免费观看| 久久精品影院6| 日韩,欧美,国产一区二区三区 | 精品久久久久久久久亚洲| 欧美潮喷喷水| 美女内射精品一级片tv| 国产极品天堂在线| 久久人妻av系列| 校园人妻丝袜中文字幕| 日韩欧美 国产精品| 女的被弄到高潮叫床怎么办| 国产精品三级大全| 三级男女做爰猛烈吃奶摸视频| 人妻系列 视频| 精品一区二区免费观看| 色5月婷婷丁香| 国产精品熟女久久久久浪| 久久精品国产自在天天线| 看片在线看免费视频| 一边摸一边抽搐一进一小说| 村上凉子中文字幕在线| 大又大粗又爽又黄少妇毛片口| 免费搜索国产男女视频| 午夜福利在线在线| 99热这里只有是精品50| 最近中文字幕高清免费大全6| 男的添女的下面高潮视频| 99热这里只有精品一区| 欧美性猛交黑人性爽| 在线天堂最新版资源| 国产精品久久久久久久久免| 久久久久久国产a免费观看| 欧美日韩综合久久久久久| 青春草国产在线视频| 亚洲精品乱久久久久久| 国内精品一区二区在线观看| 三级毛片av免费| 免费不卡的大黄色大毛片视频在线观看 | 成人亚洲欧美一区二区av| 欧美丝袜亚洲另类| 爱豆传媒免费全集在线观看| 黄色配什么色好看| 最近2019中文字幕mv第一页| 国产精品,欧美在线| 免费看光身美女| 三级毛片av免费| 亚州av有码| 国内精品一区二区在线观看| 非洲黑人性xxxx精品又粗又长| 日韩三级伦理在线观看| 99久久九九国产精品国产免费| 免费在线观看成人毛片| 亚洲av福利一区| 亚洲av二区三区四区| 女人被狂操c到高潮| 麻豆国产97在线/欧美| 插逼视频在线观看| 蜜桃亚洲精品一区二区三区| a级毛片免费高清观看在线播放| 精品国产露脸久久av麻豆 | 免费观看精品视频网站| 午夜福利在线观看吧| 哪个播放器可以免费观看大片| 美女被艹到高潮喷水动态| 国产在视频线精品| 久久精品久久久久久噜噜老黄 |