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

    河流對黃胸鼠遺傳格局的影響

    2016-12-12 10:22:29趙芳林恭華吳愛國
    四川動物 2016年2期
    關(guān)鍵詞:云南省河流種群

    趙芳, 林恭華, 吳愛國

    (1.中國科學(xué)院西北高原生物研究所,高原生物適應(yīng)與進(jìn)化重點(diǎn)實(shí)驗(yàn)室,西寧810008;2. 云南省地方病防治所,云南大理671000)

    ?

    DOI:10.11984/j.issn.1000-7083.20150286

    河流對黃胸鼠遺傳格局的影響

    趙芳1, 林恭華1, 吳愛國2*

    (1.中國科學(xué)院西北高原生物研究所,高原生物適應(yīng)與進(jìn)化重點(diǎn)實(shí)驗(yàn)室,西寧810008;2. 云南省地方病防治所,云南大理671000)

    河流可能會阻礙陸生動物的遷移并導(dǎo)致其遺傳分化,很多研究證明了河流屏障假說,但是也有少數(shù)研究不支持這個假說。本研究采集了云南省384只黃胸鼠Rattustanezumi樣本,分別來自于被5條主要河流隔開的6個地理區(qū)域。通過對mtDNA和IRBP兩個基因的測序分析,檢驗(yàn)云南省主要河流對黃胸鼠種群遺傳格局的影響。共鑒定得到81個mtDNA單倍型和135個IRBP單倍型,其中很多單倍型都被2個或2個以上區(qū)域共享。AMOVA分析顯示,大多數(shù)的變異源于種群內(nèi)部或種群之間,只有少數(shù)變異(<15%)源于6個區(qū)域之間。SAMOVA分析將mtDNA種群分成了3組,但每一組都被2個或2個以上的區(qū)域共享;對IRBP的分析則沒有檢測到明顯的譜系地理結(jié)構(gòu)。結(jié)果表明,河流屏障對黃胸鼠的遺傳格局幾乎無影響。

    河流;黃胸鼠;mtDNA;IRBP;遺傳分化

    河流屏障假說(riverine barrier hypothesis)認(rèn)為,河流作為一種地理隔離,可能會阻礙陸生動物的遷移以及基因流動,從而引起種群分化,促進(jìn)種群格局的形成(Gasconetal.,2000;Hayes & Sewlal,2004;Lietal.,2009)。針對這一假說,前人的研究得到了2種完全對立的結(jié)論。通過對哺乳動物、鳥類、爬行動物、兩棲動物和昆蟲的反復(fù)驗(yàn)證,大多數(shù)的研究支持這一假說(Knoppetal.,2011;Díaz-Muoz,2012;Fernandesetal.,2012)。然而,少數(shù)研究并不認(rèn)為河流能構(gòu)成影響陸生動物種群格局的屏障。例如,對法國阿爾卑斯山脈中山溪兩岸的一些鼩鼱Sorexaraneus遺傳結(jié)構(gòu)的研究證明,這些山溪并沒有明顯影響鼩鼱的基因流(Lugon-moulinetal.,1999)。

    大家鼠屬Rattus現(xiàn)存66個物種,廣泛分布于全球各地,通常被認(rèn)為是遷移擴(kuò)散能力最強(qiáng)的哺乳動物之一(Wilson & Reeder,2005)。作為各種鼠傳病菌(如鼠疫桿菌Yersiniapestis)和體外寄生蟲(如跳蚤)最重要的宿主,它們是各種重大動物疾病的元兇(Pagèsetal.,2010),而且還會嚴(yán)重破壞生活環(huán)境內(nèi)的建筑物,并且污染和毀壞人類的糧食儲存(Hindsetal.,2003)。因此,了解大家鼠屬生物多樣性形成的機(jī)制對疾病的預(yù)防和鼠害的防控至關(guān)重要。

    黃胸鼠Rattustanezumi是主要分布在東南亞的大家鼠,其生境廣泛,包括城區(qū)、農(nóng)村、農(nóng)田以及森林邊緣(Wilson & Reeder,2005),因攜帶很多鼠傳病菌和對人類設(shè)施巨大的破壞力而臭名昭著(鄭智民等,2008),其中最為人們所熟悉的是,它們是東南亞鼠疫的主要宿主(賀雄,王虎,2010)。黃胸鼠主要分布在熱帶和亞熱帶地區(qū),這些地區(qū)的大江大河常年不結(jié)冰。然而,這些河流是否能影響黃胸鼠的種群遺傳分化,目前還沒有系統(tǒng)的調(diào)查和分析資料。本研究中,我們從云南省捕獲了一定數(shù)量的黃胸鼠,以線粒體基因和核基因共同作為分子標(biāo)記來驗(yàn)證云南省內(nèi)5條大江(怒江、瀾滄江、李仙江、元江和南盤江)對黃胸鼠種群遺傳格局的影響。

    1 材料和方法

    1.1 樣品采集、DNA提取、基因擴(kuò)增和測序

    用活捕籠(20 cm×12 cm×9 cm)捕獲黃胸鼠,乙醚麻醉,斷頸法處死并收集大腿內(nèi)側(cè)肌肉樣品,保存在95%乙醇中,帶回實(shí)驗(yàn)室用于分子實(shí)驗(yàn)。

    采用動物組織基因組DNA提取試劑盒(Qiagen,Germany)提取基因組DNA。根據(jù)黃胸鼠線粒體基因組序列(GenBank登錄號:KF011916)設(shè)計(jì)擴(kuò)增線粒體部分序列的引物:mtDNA-F:5’-ACCCTACTTACTGGCTTCA-3’;mtDNA-R:5’-TAATGGCGAATACTGCTC-3’。mtDNA包含完整的tRNA-Tyr基因和細(xì)胞色素C氧化酶Ⅰ基因(cytochrome c oxidase subunit Ⅰ,COⅠ)的部分序列。類視黃醇結(jié)合蛋白(interphotoreceptor retinoid binding protein,IRBP)部分序列的引物根據(jù)已發(fā)表的黃胸鼠IRBP序列(GenBank登錄號:HM217693)設(shè)計(jì):IRBP-F:5’-CTATGAGCCCAGTACCCTCG-3’;IRBP-R:5’-AGCACGGACACCTGAAACA-3’。

    PCR擴(kuò)增使用40 μL反應(yīng)體系,包括:40~60 ng的基因組DNA,0.6 mmol·L-1dNTPs,0.2 μmol·L-1的正/反向引物,1 U Taq酶,以及DNA buffer和去離子水的混合反應(yīng)緩沖液。PCR反應(yīng)程序是:95 ℃預(yù)變性5 min;94 ℃變性45 s,55 ℃退火45 s,72 ℃延伸75 s,35個循環(huán);最后72 ℃延伸7 min。PCR產(chǎn)物用CASpure PCR Purification Kit(Casarray,Shanghai,China)純化,用PCR擴(kuò)增引物進(jìn)行Sanger法測序[由生工生物工程(上海)有限公司完成]。

    1.2 譜系地理學(xué)分析

    用MEGA 5(Tamuraetal.,2011)中內(nèi)嵌的Clustal W(Tompsonetal.,1997)程序比對測序得到的序列,輔以人工校對。用DnaSP 5(Librado & Rozas,2009)分析mtDNA的變異位點(diǎn)和單倍型的分布。由于IRBP基因?yàn)槎扼w序列,其中有很多單核苷酸多態(tài)性(single nucleotide polymorphisms,SNPs)位點(diǎn),我們首先使用DnaSP 5中的PHASE 2.1(Stephensetal.,2001)程序構(gòu)建單倍型,然后基于這些單倍型進(jìn)行變異位點(diǎn)和遺傳結(jié)構(gòu)分析。

    為了推斷河流對黃胸鼠譜系地理學(xué)特征的影響,按照5條河流的分隔,將采樣地點(diǎn)劃分為6個地理區(qū)域分組(G1~G6,圖1)。用Arlequin 3.11(Excoffier & Lischer,2005)進(jìn)行分子方差分析(analyses of molecular variance,AMOVA),分別計(jì)算區(qū)域之間、區(qū)域內(nèi)種群之間和種群內(nèi)部的方差百分比。

    同時用SAMOVA v1.0(Dupanloupetal.,2002)對黃胸鼠種群進(jìn)行譜系地理結(jié)構(gòu)聚類。組數(shù)K值的設(shè)定從2到6,置換500次。判定標(biāo)準(zhǔn)為:在避免某一組只有單個種群的前提下,F(xiàn)CT值最大且P值有統(tǒng)計(jì)學(xué)意義的分組為最好的種群分組(Magrietal.,2006;Shietal.,2012)。

    2 結(jié)果

    從云南省38個縣50個地點(diǎn)共捕獲384只黃胸鼠,地點(diǎn)信息詳見表1和圖1。測序結(jié)果顯示,比對好的mtDNA序列總長1 182 bp(File S1),有93個多態(tài)位點(diǎn),81個單倍型(mH1~mH81;File S2);IRBP基因總長932 bp(File S3),包括45個變異位點(diǎn)和135個通過PHASE軟件分離的SNP單倍型(iH1~iH135;File S4)。mtDNA和IRBP基因的單倍型分布情況詳見Table S1。以上所有附屬文件可在Dryad數(shù)據(jù)庫中下載(http://dx.doi.org/10.5061/dryad.584kb),或聯(lián)系作者郵件索取。

    AMOVA分析顯示,對于mtDNA而言,超過50%的變異源于種群內(nèi)部,只有13.33%源于區(qū)域之間(即G1~G6之間)。對于IRBP,大多數(shù)變異(88.12%)源于種群內(nèi)部,僅有4.62%源于區(qū)域之間(表2)。

    表1 黃胸鼠樣本的大小和地理信息

    圖1 黃胸鼠樣本及地理分布

    Fig. 1 Geographic distribution and grouping ofRattustanezumi

    縮寫見表1中地點(diǎn)編號;下圖同。

    Abbreviations are the same as the locality ID in table 1; the same below.

    表2 黃胸鼠mtDNA和IRBP的AMOVA分析結(jié)果

    SAMOVA分析顯示,對mtDNA而言,K值越大,F(xiàn)CT值越大。但是,當(dāng)K>3時,會存在某些分組內(nèi)只有一個種群的情況。因此,K=3時的模型屬于最佳模型,此模型的3個分組(S1~S3)分別包含32個、16個和2個種群(圖2)。對IRBP而言,F(xiàn)CT值隨著K值的增大而減小,同時也存在某些分組只有一個種群的情況,因此,基于IRBP數(shù)據(jù)無法得到明顯的種群結(jié)構(gòu)(表3)。

    3 討論

    本研究檢驗(yàn)了河流屏障對黃胸鼠種群遺傳的影響。單倍型分布分析顯示,河流對地理區(qū)域(G1~G6)組間的單倍型分布沒有明顯的影響??偟膩碚f,81個mtDNA單倍型中有15個是2個或2個以上區(qū)域共享的。例如,G2~G6共享mH3單倍型。同時,在135個IRBP單倍型中,有49個被2個或2個以上區(qū)域共享,其中,iH47和iH65為所有區(qū)域共享(詳見Table S1)。

    表3 黃胸鼠mtDNA和IRBP的SAMOVA分析結(jié)果

    圖2 基于mtDNA 的50個黃胸鼠種群的SAMOVA分組

    Fig. 2 SAMOVA grouping of the 50Rattustanezumipopulations based on mtDNA sequences

    AMOVA分析顯示,無論是mtDNA還是IRBP,都只有少數(shù)變異源于6個區(qū)域之間(分別為13.33%和4.67%),表明5條河流對黃胸鼠種群遺傳變異的屏障作用不明顯。同時,基于mtDNA的SAMOVA分析也顯示所有種群被分成了3組(S1~S3),且每一個SAMOVA分組都被2條或以上河流分隔開,換句話說,每個SAMOVA分組都不是單獨(dú)位于某一區(qū)域內(nèi)。此外,基于IRBP的SAMOVA分析證明沒有明顯的種群結(jié)構(gòu)存在。綜合來看,所有結(jié)果都表明云南省5條主要河流對黃胸鼠種群遺傳格局的形成影響很小。

    由于對各種生境的適應(yīng)能力很強(qiáng),以及高效的遷移和繁殖能力,大家鼠一直被認(rèn)為是最具有侵略性擴(kuò)張能力的哺乳動物之一。本屬的很多動物在有需要時都會游泳,褐家鼠R.norvegicus的游泳能力(能在水中呆3 d或游0.75 km)在很多文獻(xiàn)中被反復(fù)證明(Russelletal.,2008)。也有文獻(xiàn)(Innes,2005;Shielsetal.,2014)證明屋頂鼠R.rattus能游300~750 m到鄰近的島嶼上。作為與屋頂鼠親緣關(guān)系最近的黃胸鼠,我們有理由相信它們也是游泳能手。事實(shí)上,在福貢縣采樣時,當(dāng)?shù)鼐用窀嬖V我們黃胸鼠有時候甚至能橫游過怒江(個人交流)。云南省5條主要河流的寬度一般都不超過200 m,我們假設(shè)黃胸鼠能夠輕易游過這幾條河流的話,也就證明了對黃胸鼠而言,河流不存在屏障效應(yīng)。

    除了游泳外,黃胸鼠還可能通過船只往來進(jìn)行遷移(Russelletal.,2008)。由于交通業(yè)的發(fā)展,目前世界上大部分地區(qū)都有褐家鼠和屋頂鼠的分布(Krinke,2000;Aplinetal.,2011)。在中國,火車運(yùn)輸對黃胸鼠在各個省份之間的遷移扮演了重要的角色(張美文等,2000)。同時,橋梁也為黃胸鼠橫跨河流兩岸提供了路徑,然而這些跨河交通及設(shè)施究竟如何影響黃胸鼠的種群遺傳結(jié)構(gòu),則需要進(jìn)一步研究。

    賀雄, 王虎. 2010. 現(xiàn)代鼠疫概論[M]. 北京: 科學(xué)出版社.

    張美文, 郭聰, 王勇, 等. 2000. 我國黃胸鼠的研究現(xiàn)狀[J]. 動物學(xué)研究, 21(6): 487-497.

    鄭智民, 姜志寬, 陳安國. 2008. 嚙齒動物學(xué)[M]. 上海: 上海交通大學(xué)出版社.

    Aplin KP, Suzuki H, Chinen AA,etal. 2011. Multiple geographic origins of commensalism and complex dispersal history of black rats[J]. PLoS ONE, 6(11): e26357.

    Dupanloup I, Schneider S, Excoffier L. 2002. A simulated annealing approach to define the genetic structure of populations[J]. Molecular Ecology, 11: 2571-2581.

    Excoffier L, Lischer HL. 2010. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows[J]. Molecular Ecology Resources, 10: 564-567.

    Fernandes AM, Wink M, Aleixo A. 2012. Phylogeography of the chestnut-tailed antbird (Myrmecizahemimelaena) clarifies the role of rivers in Amazonian biogeography[J]. Journal of Biogeography, 39(8):1524-1535.

    Gascon C, Malcolm JR, Patton JL,etal. 2000. Riverine barriers and the geographic distribution of Amazonian species[J]. Proceedings of the National Academy of Sciences of the United States of America, 97(25): 13672-13677.

    Hayes FE, Sewlal JAN. 2004. The Amazon River as a dispersal barrier to passerine birds: effects of river width, habitat and taxonomy[J]. Journal of Biogeography, 31(11): 1809-1818.

    Hinds LA, Krebs CJ, Spratt DM. 2003. Rats, mice and people: rodent biology and management[M]. Canberra, Australia: Australian Centre for International Agricultural Research.

    Innes JG. 2005. Ship rat[M]// King CA. The handbook of New Zealand mammals 2nd edition. Oxford, UK: Oxford University Press: 187-203.

    Knopp T, Rahagalala P, Miinala M,etal. 2011. Current geographical ranges of Malagasy dung beetles are not delimited by large rivers[J]. Journal of Biogeography, 38(6): 1098-1108.

    Krinke GJ. 2000. The laboratory rat[M]. San Diego, USA: Academic Press.

    Li R, Chen W, Tu L,etal. 2009. Rivers as barriers for high elevation amphibians: a phylogeographic analysis of the alpine stream frog of the Hengduan Mountains[J]. Journal of Zoology, 277(4): 309-316.

    Librado P, Rozas J. 2009. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data[J]. Bioinformatics, 25: 1451-1452.

    Lugon-moulin N, Bruenner H, Balloux F,etal. 1999. Do riverine barriers, history or introgression shape the genetic structuring of a common shrew (Sorexaraneus) population?[J]. Heredity, 83: 155-161.

    Magri D, Vendramin GG, Comps B,etal. 2006. A new scenario for the quaternary history of European beech populations: palaeobotanical evidence and genetic consequences[J]. New Phytologist, 171: 199-221.

    Pagès M, Chaval Y, Herbreteau V,etal. 2010. Revisiting the taxonomy of the Rattini tribe: a phylogeny-based delimitation of species boundaries[J]. BMC Evolutionary Biology, 10(1): 184.

    Russell JC, Towns DR, Clout MN. 2008. Review of rat invasion biology: implications for island biosecurity (Science for Conservation 286) [M]. Wellington, NZ: Science & Technical Publishing.

    Shi W, Kerdelhué C, Ye H. 2012. Genetic structure and inferences on potential source areas forBactroceradorsalis(Hendel) based on mitochondrial and microsatellite markers[J]. PLoS ONE, 7(5): e37083.

    Shiels AB, William CP, Robert TS,etal. 2014. Biology and impacts of pacific island invasive species. 11.Rattusrattus, the black rat (Rodentia: Muridae)[J]. Pacific Science, 68(2): 145-184.

    Stephens M, Smith NJ, Donnelly P. 2001. A new statistical method for haplotype reconstruction from population data[J]. The American Journal of Human Genetics, 68: 978-989.

    Tamura K, Peterson D, Peterson N,etal. 2011. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods[J]. Molecular Biology and Evolution, 28(10): 2731-2739.

    Tompson JD, Gibson TJ, Plewniak F,etal. 1997. The Clustal X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools[J]. Nucleic Acids Research, 25(24): 4876-4882.

    Wilson DE, Reeder DA. 2005. Mammal species of the world: a taxonomic and geographic reference[M]. Baltimore, USA: John Hopkins University Press.

    Influence of Riverine Barriers on Genetic Structure ofRattustanezumi

    ZHAO Fang1, LIN Gonghua1, WU Aiguo2*

    (1. Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences,Xining 810008, China; 2. Yunnan Institute of Endemic Disease Control and Prevention, Dali, Yunnan Province 671000, China)

    Rivers may hinder population migration and promote genetic differentiation of terrestrial animals. Most studies testing the riverine barrier hypothesis have yielded positive results, while very few studies do not support this hypothesis. In this study, 384 Asian house rats (Rattustanezumi) were collected from Yunnan, China. By combined use of the sequence markers of mitochondrial (mtDNA) and nuclear (IRBP), the influence of five local rivers on the population genetics ofR.tanezumiwas tested. A total of 81 mtDNA haplotypes and 135 IRBP haplotypes were determined, and many of which were shared by two or more of the six geographic regions. AMOVA analyses showed that most variations happened either within populations or among populations, while very few variations (<15%) happened among the six geographic regions formed by the five rivers. SAMOVA analysis classified the mtDNA sequences into three groups, while no genetic differentiation was detected on IRBP. Our results concluded that riverine barriers in the study area had little influence on the genetic structure ofR.tanezumi.

    river;Rattustanezumi; mtDNA; IRBP; genetic differentiation

    2015-09-18 接受日期:2015-11-04

    國家自然科學(xué)基金面上項(xiàng)目(81260416, 31101628); 中國科學(xué)院青年創(chuàng)新促進(jìn)會項(xiàng)目(2015352)

    趙芳(1987—), 女, 博士研究生, 研究方向: 動物生態(tài)學(xué), E-mail:zf_lgh@163.com

    *通信作者Corresponding author, E-mail:wuaiguo698@163.com

    Q959.8; Q3

    A

    1000-7083(2016)02-0167-06

    猜你喜歡
    云南省河流種群
    邢氏水蕨成功繁衍并建立種群 等
    云南省安寧市老年大學(xué)之歌
    山西省發(fā)現(xiàn)刺五加種群分布
    民國初期云南省議會述論
    河流
    小太陽畫報(2019年4期)2019-06-11 10:29:48
    流放自己的河流
    散文詩(2018年20期)2018-05-06 08:03:44
    云南省瑞麗市老年大學(xué)之歌
    近31年來云南省潛在蒸散量的時空演變
    當(dāng)河流遇見海
    崗更湖鯉魚的種群特征
    精品久久久久久久久av| 国产亚洲精品久久久com| 久久综合国产亚洲精品| 天堂8中文在线网| xxxhd国产人妻xxx| 夜夜骑夜夜射夜夜干| 欧美 日韩 精品 国产| 久久毛片免费看一区二区三区| 国产成人精品在线电影| 街头女战士在线观看网站| 国产精品成人在线| 精品人妻一区二区三区麻豆| 午夜福利视频精品| 黑人猛操日本美女一级片| 国产精品三级大全| 3wmmmm亚洲av在线观看| 制服诱惑二区| 成人黄色视频免费在线看| 视频区图区小说| 国产一区二区在线观看av| 精品久久蜜臀av无| 夫妻午夜视频| 免费观看性生交大片5| 丰满少妇做爰视频| 2018国产大陆天天弄谢| 男男h啪啪无遮挡| 嫩草影院入口| 特大巨黑吊av在线直播| 热99久久久久精品小说推荐| av播播在线观看一区| 丰满乱子伦码专区| 国产精品一二三区在线看| 国产av一区二区精品久久| 国产女主播在线喷水免费视频网站| 一本—道久久a久久精品蜜桃钙片| 精品卡一卡二卡四卡免费| 久久久久久久大尺度免费视频| 久久亚洲国产成人精品v| 亚洲精品自拍成人| 免费大片18禁| 中文字幕免费在线视频6| 精品少妇黑人巨大在线播放| 欧美bdsm另类| 午夜老司机福利剧场| 乱人伦中国视频| 亚洲经典国产精华液单| 亚洲国产av新网站| 国产免费一区二区三区四区乱码| 久久久a久久爽久久v久久| 美女脱内裤让男人舔精品视频| 亚洲伊人久久精品综合| 王馨瑶露胸无遮挡在线观看| av黄色大香蕉| 国产成人精品无人区| 又粗又硬又长又爽又黄的视频| 国产精品久久久久成人av| 涩涩av久久男人的天堂| 高清黄色对白视频在线免费看| xxx大片免费视频| 国内精品宾馆在线| 80岁老熟妇乱子伦牲交| 只有这里有精品99| av在线播放精品| 国产一级毛片在线| 国产男女超爽视频在线观看| 又粗又硬又长又爽又黄的视频| 中文字幕精品免费在线观看视频 | 伦精品一区二区三区| 三上悠亚av全集在线观看| 日本欧美视频一区| av网站免费在线观看视频| 亚洲国产av影院在线观看| 亚洲婷婷狠狠爱综合网| 国产亚洲精品第一综合不卡 | 国产成人91sexporn| 男男h啪啪无遮挡| 国产乱人偷精品视频| 久久久久久久久久人人人人人人| 亚洲国产av影院在线观看| 一级二级三级毛片免费看| 亚洲精品成人av观看孕妇| 熟女av电影| 91午夜精品亚洲一区二区三区| 一级毛片电影观看| 国产精品蜜桃在线观看| 免费av不卡在线播放| 亚洲综合色惰| 国产精品99久久99久久久不卡 | 精品少妇内射三级| 久久久久久人妻| 女人久久www免费人成看片| 免费观看a级毛片全部| av国产精品久久久久影院| 国产爽快片一区二区三区| 国产成人精品在线电影| 色哟哟·www| 亚洲人与动物交配视频| 日韩,欧美,国产一区二区三区| 美女主播在线视频| 欧美 亚洲 国产 日韩一| 国产高清三级在线| 一区二区三区精品91| 日韩成人av中文字幕在线观看| 特大巨黑吊av在线直播| 亚洲国产精品一区二区三区在线| 亚洲欧洲日产国产| 欧美 日韩 精品 国产| 高清视频免费观看一区二区| 久久久久久伊人网av| 99久久精品国产国产毛片| 国产淫语在线视频| av免费观看日本| 少妇猛男粗大的猛烈进出视频| 亚洲不卡免费看| av.在线天堂| 久久人人爽人人爽人人片va| 国产视频内射| 女性生殖器流出的白浆| 亚洲伊人久久精品综合| 国产精品一区www在线观看| 国产成人午夜福利电影在线观看| 一区二区av电影网| 最近的中文字幕免费完整| 高清视频免费观看一区二区| 18禁裸乳无遮挡动漫免费视频| 美女国产高潮福利片在线看| 国产在线视频一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 51国产日韩欧美| 欧美+日韩+精品| 成人国产麻豆网| 一本大道久久a久久精品| 秋霞伦理黄片| 日韩不卡一区二区三区视频在线| 最近2019中文字幕mv第一页| 伦精品一区二区三区| 一区二区av电影网| 免费日韩欧美在线观看| 伊人久久精品亚洲午夜| 日本wwww免费看| 永久网站在线| 午夜影院在线不卡| 黄色一级大片看看| 黑人高潮一二区| 男人爽女人下面视频在线观看| 亚洲av成人精品一区久久| 中文天堂在线官网| 亚洲成人手机| 欧美变态另类bdsm刘玥| 在线观看免费日韩欧美大片 | 三级国产精品欧美在线观看| 久久久久久久精品精品| 十八禁高潮呻吟视频| 三级国产精品欧美在线观看| 日日爽夜夜爽网站| 欧美变态另类bdsm刘玥| 亚洲av男天堂| 人妻系列 视频| 欧美日本中文国产一区发布| 国产亚洲午夜精品一区二区久久| av线在线观看网站| 国产精品偷伦视频观看了| av卡一久久| 国产免费又黄又爽又色| 又大又黄又爽视频免费| 国产成人aa在线观看| 伊人久久精品亚洲午夜| 高清av免费在线| 亚洲成人一二三区av| 国产男女超爽视频在线观看| 99九九线精品视频在线观看视频| 在线观看人妻少妇| 只有这里有精品99| 2021少妇久久久久久久久久久| 日韩欧美一区视频在线观看| 亚洲天堂av无毛| 午夜福利影视在线免费观看| 久久鲁丝午夜福利片| 桃花免费在线播放| 蜜桃国产av成人99| 久久免费观看电影| av电影中文网址| 少妇被粗大的猛进出69影院 | 亚洲综合色网址| 一边摸一边做爽爽视频免费| 你懂的网址亚洲精品在线观看| 亚洲精品国产色婷婷电影| 午夜精品国产一区二区电影| 人妻夜夜爽99麻豆av| 欧美日韩视频高清一区二区三区二| 涩涩av久久男人的天堂| 国产熟女午夜一区二区三区 | 中文欧美无线码| 久久午夜福利片| 一区二区日韩欧美中文字幕 | 欧美精品国产亚洲| 老熟女久久久| 亚洲美女黄色视频免费看| 能在线免费看毛片的网站| 国产午夜精品一二区理论片| 国产欧美日韩综合在线一区二区| 麻豆精品久久久久久蜜桃| 免费av不卡在线播放| 丰满乱子伦码专区| 国产成人精品一,二区| 只有这里有精品99| a 毛片基地| 少妇高潮的动态图| 国产精品三级大全| 国产精品麻豆人妻色哟哟久久| 欧美日韩成人在线一区二区| 在线亚洲精品国产二区图片欧美 | 亚洲人与动物交配视频| 人体艺术视频欧美日本| 91在线精品国自产拍蜜月| 亚洲欧美色中文字幕在线| 制服诱惑二区| 国产片特级美女逼逼视频| 岛国毛片在线播放| 男女高潮啪啪啪动态图| 午夜影院在线不卡| 视频区图区小说| 51国产日韩欧美| 久久精品国产亚洲av涩爱| 自拍欧美九色日韩亚洲蝌蚪91| 久热久热在线精品观看| 国产免费一级a男人的天堂| 久久97久久精品| 精品人妻在线不人妻| 视频区图区小说| 免费黄色在线免费观看| 国产女主播在线喷水免费视频网站| av在线播放精品| 日韩电影二区| 一级毛片 在线播放| 韩国高清视频一区二区三区| 91久久精品电影网| 插阴视频在线观看视频| 麻豆成人av视频| 午夜激情久久久久久久| 91成人精品电影| 国产精品一国产av| 51国产日韩欧美| 99视频精品全部免费 在线| 成人18禁高潮啪啪吃奶动态图 | 久久精品国产自在天天线| 黄色视频在线播放观看不卡| 日本-黄色视频高清免费观看| 日本爱情动作片www.在线观看| 日韩精品免费视频一区二区三区 | 国产精品免费大片| 老女人水多毛片| 精品久久久久久久久亚洲| 欧美日韩av久久| 大陆偷拍与自拍| 天美传媒精品一区二区| 中国美白少妇内射xxxbb| 亚洲av福利一区| 高清毛片免费看| 久久99精品国语久久久| 五月天丁香电影| 久久久久久久久久久免费av| 99九九线精品视频在线观看视频| 午夜视频国产福利| 午夜久久久在线观看| 高清在线视频一区二区三区| 99热国产这里只有精品6| 国产乱来视频区| 肉色欧美久久久久久久蜜桃| 大片免费播放器 马上看| 午夜福利,免费看| 亚洲综合色网址| 久久久久久久久久久丰满| 久久99精品国语久久久| 99热6这里只有精品| 老司机影院毛片| 国产综合精华液| 中文字幕久久专区| 精品午夜福利在线看| 九九爱精品视频在线观看| 欧美日韩在线观看h| 亚洲国产色片| 在线看a的网站| 18禁在线无遮挡免费观看视频| 人人妻人人爽人人添夜夜欢视频| 狂野欧美激情性xxxx在线观看| 熟妇人妻不卡中文字幕| 精品久久久久久久久亚洲| 国产乱来视频区| 男的添女的下面高潮视频| 亚洲精品国产av成人精品| 国产在线免费精品| 亚洲精品久久午夜乱码| 久久久久久久久久成人| 国产综合精华液| 另类亚洲欧美激情| 精品久久国产蜜桃| 热re99久久精品国产66热6| 97在线人人人人妻| 考比视频在线观看| 成人毛片60女人毛片免费| 大片电影免费在线观看免费| 国产黄色视频一区二区在线观看| 欧美激情 高清一区二区三区| 观看av在线不卡| 久久国内精品自在自线图片| 久久久欧美国产精品| 国产 一区精品| 中文字幕av电影在线播放| 国产一区有黄有色的免费视频| 国产女主播在线喷水免费视频网站| 18禁在线播放成人免费| 自线自在国产av| 国产免费一级a男人的天堂| 久久人人爽av亚洲精品天堂| 91精品国产九色| 亚洲欧美精品自产自拍| 我要看黄色一级片免费的| .国产精品久久| 中文字幕av电影在线播放| 欧美人与善性xxx| 狂野欧美激情性xxxx在线观看| 精品一品国产午夜福利视频| 国产亚洲午夜精品一区二区久久| 成人午夜精彩视频在线观看| 国产片内射在线| 美女中出高潮动态图| 国产一区亚洲一区在线观看| 黄色欧美视频在线观看| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品国产精品| 美女国产视频在线观看| 精品人妻一区二区三区麻豆| 十分钟在线观看高清视频www| 成人二区视频| 亚洲天堂av无毛| 热re99久久精品国产66热6| 国产精品国产三级专区第一集| 亚洲欧美一区二区三区国产| 亚洲精品日韩在线中文字幕| 99热这里只有是精品在线观看| 婷婷色麻豆天堂久久| 黄色视频在线播放观看不卡| 亚洲精品自拍成人| 日韩一区二区三区影片| 欧美亚洲 丝袜 人妻 在线| 亚洲五月色婷婷综合| 男女啪啪激烈高潮av片| 在线免费观看不下载黄p国产| 两个人的视频大全免费| 日韩 亚洲 欧美在线| 午夜福利在线观看免费完整高清在| 丰满少妇做爰视频| 人体艺术视频欧美日本| 国内精品宾馆在线| 中文天堂在线官网| 黄色视频在线播放观看不卡| 亚洲av国产av综合av卡| 日本黄色片子视频| 久久99蜜桃精品久久| 韩国高清视频一区二区三区| 日本色播在线视频| 熟妇人妻不卡中文字幕| 亚洲av.av天堂| 伊人久久国产一区二区| 爱豆传媒免费全集在线观看| 欧美97在线视频| 大码成人一级视频| av专区在线播放| 国产精品欧美亚洲77777| 亚洲精品久久成人aⅴ小说 | 日本91视频免费播放| 色婷婷av一区二区三区视频| 亚洲av成人精品一区久久| 五月玫瑰六月丁香| 人妻 亚洲 视频| 国产精品免费大片| 国产精品不卡视频一区二区| 观看美女的网站| 国产探花极品一区二区| 男男h啪啪无遮挡| 国产成人精品在线电影| 高清欧美精品videossex| av.在线天堂| 热99国产精品久久久久久7| 午夜福利在线观看免费完整高清在| 国产精品熟女久久久久浪| 狠狠婷婷综合久久久久久88av| 精品国产国语对白av| 大陆偷拍与自拍| 亚洲欧美清纯卡通| 亚洲国产精品一区二区三区在线| 亚洲欧美一区二区三区国产| 91国产中文字幕| 亚洲av综合色区一区| 亚洲熟女精品中文字幕| 女的被弄到高潮叫床怎么办| 日韩精品有码人妻一区| 丁香六月天网| 色婷婷av一区二区三区视频| 18禁动态无遮挡网站| 黄片无遮挡物在线观看| 久久久久久人妻| 免费人妻精品一区二区三区视频| 99久久综合免费| 国产成人精品在线电影| 日本欧美国产在线视频| 91久久精品电影网| 国产视频内射| 国产欧美日韩综合在线一区二区| 丁香六月天网| 色哟哟·www| 日韩在线高清观看一区二区三区| 最近中文字幕高清免费大全6| 国产精品一区二区三区四区免费观看| 一级毛片我不卡| 日本-黄色视频高清免费观看| 久久久久网色| freevideosex欧美| 少妇 在线观看| 在线观看免费日韩欧美大片 | 五月玫瑰六月丁香| 国产又色又爽无遮挡免| av国产精品久久久久影院| 国产免费现黄频在线看| av不卡在线播放| 永久免费av网站大全| 女性被躁到高潮视频| 日韩av在线免费看完整版不卡| 秋霞在线观看毛片| 久久精品国产自在天天线| 精品人妻一区二区三区麻豆| 国产国拍精品亚洲av在线观看| 一级毛片aaaaaa免费看小| av天堂久久9| 成年人免费黄色播放视频| 美女福利国产在线| 午夜福利在线观看免费完整高清在| 黄色视频在线播放观看不卡| 午夜av观看不卡| 亚洲精品国产av蜜桃| 极品人妻少妇av视频| 免费少妇av软件| 国产精品欧美亚洲77777| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | av在线观看视频网站免费| 晚上一个人看的免费电影| 日本vs欧美在线观看视频| 精品一区二区三卡| 久久综合国产亚洲精品| 国产男女超爽视频在线观看| 丝袜脚勾引网站| 只有这里有精品99| 久久精品国产a三级三级三级| 各种免费的搞黄视频| 午夜视频国产福利| 亚洲欧美日韩卡通动漫| 午夜福利在线观看免费完整高清在| 少妇被粗大的猛进出69影院 | 亚洲久久久国产精品| 制服诱惑二区| 一本一本综合久久| 日韩精品免费视频一区二区三区 | 免费看光身美女| 亚洲精品国产色婷婷电影| 精品久久久噜噜| 一级毛片 在线播放| 亚洲av不卡在线观看| 视频在线观看一区二区三区| 精品人妻熟女av久视频| 精品少妇黑人巨大在线播放| 大片电影免费在线观看免费| 午夜福利视频在线观看免费| 久久精品久久久久久噜噜老黄| 国产无遮挡羞羞视频在线观看| 亚洲av中文av极速乱| 国产 精品1| 边亲边吃奶的免费视频| 色婷婷久久久亚洲欧美| 成年美女黄网站色视频大全免费 | 久久国产精品大桥未久av| 国产免费视频播放在线视频| 看免费成人av毛片| 观看美女的网站| 亚洲久久久国产精品| 青青草视频在线视频观看| 国产成人一区二区在线| 久久久久人妻精品一区果冻| 午夜免费男女啪啪视频观看| 黑人猛操日本美女一级片| 亚洲丝袜综合中文字幕| 日韩三级伦理在线观看| 新久久久久国产一级毛片| 国产欧美日韩综合在线一区二区| 亚洲欧美日韩另类电影网站| 久久久亚洲精品成人影院| 插阴视频在线观看视频| kizo精华| 亚洲精品国产色婷婷电影| 纯流量卡能插随身wifi吗| 伊人亚洲综合成人网| 嫩草影院入口| 国产又色又爽无遮挡免| 一级二级三级毛片免费看| 一个人看视频在线观看www免费| 91在线精品国自产拍蜜月| 日日摸夜夜添夜夜添av毛片| av一本久久久久| 久久这里有精品视频免费| 午夜激情久久久久久久| 夜夜爽夜夜爽视频| 国产精品人妻久久久影院| 在线观看一区二区三区激情| 亚洲综合色网址| 精品视频人人做人人爽| 久久久久久久久大av| 亚洲图色成人| 精品久久久久久电影网| 精品99又大又爽又粗少妇毛片| 超色免费av| videossex国产| 亚洲中文av在线| 最近中文字幕高清免费大全6| 97超碰精品成人国产| 一级毛片电影观看| 精品一区二区三卡| 七月丁香在线播放| 妹子高潮喷水视频| 日韩av在线免费看完整版不卡| 美女cb高潮喷水在线观看| 亚洲中文av在线| 美女内射精品一级片tv| 成人免费观看视频高清| 天美传媒精品一区二区| 伊人久久精品亚洲午夜| 亚洲av电影在线观看一区二区三区| 久久精品熟女亚洲av麻豆精品| 亚洲国产欧美在线一区| 精品亚洲乱码少妇综合久久| 最新的欧美精品一区二区| 我要看黄色一级片免费的| 亚洲成人av在线免费| 免费观看无遮挡的男女| 国产精品.久久久| 最后的刺客免费高清国语| 夜夜看夜夜爽夜夜摸| 亚洲四区av| 高清黄色对白视频在线免费看| 插阴视频在线观看视频| 精品熟女少妇av免费看| 久久国内精品自在自线图片| 夫妻性生交免费视频一级片| 久久国产亚洲av麻豆专区| 黄色一级大片看看| 毛片一级片免费看久久久久| 在线天堂最新版资源| 热re99久久精品国产66热6| 一本色道久久久久久精品综合| 欧美+日韩+精品| 色婷婷久久久亚洲欧美| av.在线天堂| 亚洲精品中文字幕在线视频| 一级片'在线观看视频| 亚洲中文av在线| 少妇人妻精品综合一区二区| av免费在线看不卡| 99re6热这里在线精品视频| 中文乱码字字幕精品一区二区三区| 日韩伦理黄色片| 欧美变态另类bdsm刘玥| 人妻人人澡人人爽人人| 国产 一区精品| 在线观看免费视频网站a站| 亚洲婷婷狠狠爱综合网| 大香蕉久久成人网| 精品久久久噜噜| 久久韩国三级中文字幕| 日产精品乱码卡一卡2卡三| 国产一区有黄有色的免费视频| xxxhd国产人妻xxx| 精品人妻偷拍中文字幕| videos熟女内射| 天美传媒精品一区二区| 亚洲av欧美aⅴ国产| av线在线观看网站| 黄片播放在线免费| 欧美精品一区二区大全| 国产乱人偷精品视频| 观看美女的网站| 22中文网久久字幕| 国产亚洲精品第一综合不卡 | 高清欧美精品videossex| 欧美日韩精品成人综合77777| 欧美变态另类bdsm刘玥| 高清欧美精品videossex| 少妇被粗大猛烈的视频| 18+在线观看网站| 99久国产av精品国产电影| 少妇被粗大猛烈的视频| kizo精华| 在线 av 中文字幕| 搡女人真爽免费视频火全软件| 亚洲国产精品999| av女优亚洲男人天堂| 极品人妻少妇av视频| 精品少妇黑人巨大在线播放| 亚洲欧美一区二区三区国产| 成人国产av品久久久| 搡老乐熟女国产| 午夜福利在线观看免费完整高清在| 国产爽快片一区二区三区| 18禁在线播放成人免费| 91久久精品电影网| 国产亚洲av片在线观看秒播厂| 亚洲精品av麻豆狂野|