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

    基于微衛(wèi)星標記的不同殼型湖北釘螺小尺度景觀群體遺傳結(jié)構(gòu)研究

    2014-04-02 11:43:20李石柱
    中國人獸共患病學報 2014年7期
    關(guān)鍵詞:螺殼釘螺遺傳變異

    崔 斌,官 威,尤 平,李石柱

    湖北釘螺(Oncomelaniahupensis) 是日本血吸蟲唯一的中間宿主,我國大陸僅有一種,廣泛分布于長江中下游地區(qū)[1-2]。由于其分布廣泛、孳生環(huán)境復雜以及地理隔離的影響,湖北釘螺不同地理群體發(fā)生了顯著的遺傳分化和變異。近年來,隨著湖北釘螺的群體遺傳學、景觀生態(tài)學等深入研究,我國大陸分布的釘螺地理種群主要劃分為長江中下游的湖沼型景觀種群、高山型景觀種群、喀斯特地貌景觀種群和濱海山丘型景觀種群[3],其中湖北釘螺湖沼型景觀種群主要分布在長江中下游的湖南、湖北、江西、安徽和江蘇等省的湖沼、水網(wǎng)和丘陵等地區(qū),這一地區(qū)景觀范圍廣,水網(wǎng)密布,植被豐富,且雨量充沛,是遺傳多樣性最為復雜的一個種群,甚至同一地區(qū)不同微生境種群已表現(xiàn)出明顯的遺傳分化[4]。

    湖北松滋地區(qū)位于湖北省西南部,地形西高東低,從西南海拔500~800 m較高區(qū)的山地向海拔50 m以下的江漢平原湖區(qū)呈四階梯遞降,境內(nèi)河渠縱橫,間有湖泊,其境內(nèi)有螺環(huán)境多樣,主要包括垸內(nèi)溝渠型和山丘型等,且不同類型環(huán)境下釘螺隨生境變化呈現(xiàn)出不同螺殼形態(tài),這種在小尺度生境湖北釘螺種群中呈現(xiàn)明顯表型差異的現(xiàn)象并不多見。微衛(wèi)星DNA標記是廣泛分布于真核生物基因組中的一類簡單的短核苷酸串聯(lián)重復序列,具有多態(tài)性高、重復性高、呈共顯性表達、易檢測和符合孟德爾分離等優(yōu)點,因此被廣泛應用于群體遺傳多樣性分析[5-6]。國內(nèi)應用微衛(wèi)星標記進行湖北釘螺遺傳多樣性的研究開展較晚,且僅有少量應用微衛(wèi)星錨定進行群體遺傳變異研究的報告[7-9]。本研究應用多態(tài)性微衛(wèi)星DNA標記位點,分析了湖北釘螺松滋地區(qū)小尺度生境種群遺傳結(jié)構(gòu),以探討其遺傳分化程度。

    1 材料和方法

    1.1釘螺樣本來源 湖北釘螺采自松滋市境內(nèi)的4個鄉(xiāng)鎮(zhèn)10個生境,采集地分布見表1。每個采集地釘螺在解剖鏡下觀察并記錄螺殼形態(tài),所有釘螺經(jīng)室內(nèi)飼養(yǎng)約1周后,用逸蚴法鑒定釘螺是否感染血吸蟲,取陰性釘螺作為實驗材料,75%乙醇浸泡,4 ℃保存待用。

    1.2方法

    1.2.1形態(tài)學鑒定 目測釘螺螺殼形態(tài)初步分類后,采用解剖鏡下觀察鑒定并拍照。

    1.2.2分子檢測

    1.2.2.1DNA抽提 取單只湖北釘螺的腹足肌肉組織30 mg,采用軟體動物DNA小量提取試劑盒(購自美國OMEGA生物技術(shù)公司)提取基因組DNA,-20 ℃保存待用,操作方法參考文獻[10]。

    1.2.2.2熒光標記引物 多態(tài)性微衛(wèi)星位點的篩選參照文獻[11-12],獲得的多態(tài)性微衛(wèi)星DNA位點的正向引物進行6-FAM、HEX和NED熒光標記(均購自上海生工生物技術(shù)有限公司)(表2)。

    1.2.2.3PCR擴增 擴增反應體系為25 μL,包括20 ng模板DNA、0.5 μmol/L引物、1U的Taq酶、0.15 mmol/L的dNTPs、1.5 mmol/L的MgCl2和1×PCR緩沖液。PCR擴增條件為:95 ℃ 5 min;95 ℃ 1 min,50~55 ℃ 45 s,72 ℃ 1.5 min,共35個循環(huán);72 ℃ 5 min。

    1.2.2.4擴增產(chǎn)物檢測 取1 μL PCR產(chǎn)物,0.6 μL ROX標準物和8.4 μL Hi-Di甲酰胺混勻,95 ℃變性5 min,隨后用全自動遺傳分析儀檢測PCR產(chǎn)物(3730XL,USA, ABI )。

    表1 湖北釘螺樣本采集地分布

    表2 湖北釘螺6個微衛(wèi)星位點的擴增引物

    1.2.3數(shù)據(jù)分析 基因掃描結(jié)果經(jīng)Genescan和Genotype分析軟件獲得擴增片段的精確長度,結(jié)果以Excel表格形式輸出,使用MS-Tools的Excel Microsatellite Toolkit轉(zhuǎn)化原始數(shù)據(jù)為FSTAT2.9.3.2可識別文件,計算群體內(nèi)的遺傳分化指標平均等位基因數(shù)(mean number alleles,MNA)、期望雜合度(expected heterozygosity,He)和觀察雜合度(observed heterozygosity,Ho)。應用POPGENE1.3.2軟件分別檢驗哈代-溫伯格平衡(Hardy-Weinberg equilibrium,HWE)、F-statistics固定指數(shù)(包括Fit、Fst和Fis)和遺傳距離。應用ARLEQUIN3.11[13]軟件計算成對群體間的基因流(Nm),利用AWOVA模塊分析群體遺傳變異的組成,再用POPGENE1.3.2軟件以非加權(quán)組平均法(Unweighted Pair Group Method with Arith-metic means,UPGMA)進行聚類分析[14]。在GENEPOP中對遺傳距離與地理距離的相關(guān)性進行回歸的Mantel檢驗,根據(jù)Botstein等[15]公式計算多態(tài)信息含量(Polymorphism Information Content,PIC)。

    2 結(jié) 果

    2.1形態(tài)學鑒定結(jié)果 本研究共調(diào)查了松滋市10個生態(tài)環(huán)境,獲得湖北釘螺樣本400余只,解剖鏡下鑒別不同生境下的釘螺螺殼形態(tài)可見,螺殼形態(tài)從山丘環(huán)境的光殼逐步向溝渠環(huán)境的淺肋、深肋演變,其中王家橋鎮(zhèn)、陳店鎮(zhèn)4個山丘環(huán)境采獲的釘螺為光殼,新江口鎮(zhèn)3個環(huán)境采獲的釘螺同時存在光殼和肋殼(淺肋和深肋),老城鎮(zhèn)3個溝渠環(huán)境為肋殼(圖1)。

    圖1本研究中不同采集地湖北釘螺形態(tài)

    Fig.1ShelltypesofO.hupensisinthisstudy

    From left in line 1: SA, SB, SC, SD, SE; line 2: SF, SG, SH, SI, SK.

    2.2基因掃描結(jié)果

    2.2.1微衛(wèi)星DNA位點多態(tài)信息含量 6個微衛(wèi)星位點在10個群體中一共檢測到141個等位基因,其中B14位點上最多,為34個,D11和T617位點最少,為20個,平均每個位點檢測到23.5個。6個微衛(wèi)星位點的平均有效等位基因數(shù)為1.575,各位點的有效等位基因數(shù)并不平衡,且數(shù)值差異較大,范圍從0.445至3.060。

    2.2.2群體內(nèi)遺傳變異 10個群體中,觀測雜合度較高,在0.59~0.89之間,其中溝渠環(huán)境肋殼釘螺群體觀察雜合度顯著高于山丘型環(huán)境的光殼釘螺群體,最高者為老城鎮(zhèn)鎮(zhèn)SH群體(0.892),最低者為陳店鎮(zhèn)SC群體(0.430)。各群體基因多態(tài)信息含量較高,PIC在0.528-0.857之間,平均PIC為0.732。其中溝渠環(huán)境肋殼釘螺群體顯著高于山丘型環(huán)境的光殼釘螺群體,最高者為老城鎮(zhèn)鎮(zhèn)SH群體(0.857),最低者為新江口鎮(zhèn)SF群體(0.528)。各群體內(nèi)近交系數(shù)Fis較低,范圍在0.056-0.638間,平均為0.301,山丘型光殼釘螺群體顯著高于溝渠環(huán)境肋殼釘螺群體。各群體的遺傳變異系數(shù)見表3。

    表3 松滋釘螺遺傳變異系數(shù)

    2.2.3群體間遺傳差異

    2.2.3.1F-統(tǒng)計量 F-統(tǒng)計量中3 個固定指數(shù)群體內(nèi)近交系數(shù)(Fis)、群體間基因分化系數(shù)(Fst)和總?cè)后w的近交系數(shù)(Fit)均是反映群體近交程度或群體間遺傳分化的指標。本研究F-統(tǒng)計量結(jié)果見表4,總?cè)后w和群體內(nèi)近交系數(shù)較高,分別為40.35%和30.32%,表明群體內(nèi)近交嚴重;各亞群體間基因分化系數(shù)(Fst)為14.40%,表明群體間存在一定遺產(chǎn)分化,但分化程度并不平衡,其中SC與SG群體間Fst值最大,為25.25%,SF與SI群體間遺傳分化最小,為-1.56%。群體SE與SK間的基因流最大,為19.78,群體SI和SF間基因流最小,為-16.23。

    表4 所有樣本F-statistics固定指數(shù)

    表5 遺傳分化程度(FST)(下)和基因流(Nm)(上)

    Note: *P<0.05; **P<0.01。群體代號同表1。Population codes see Table 1.

    表6 遺傳距離和地理距離

    注: 對角線右上:地理距離(KM);對角線左下;遺傳距離。Diagonal right: geographical distance (KM); diagonal left; genetic distance.

    2.2.3.2遺傳距離和地理距離 對群體間遺傳距離與地理距離的相關(guān)性進行回歸的Mantel檢驗,結(jié)果顯示,兩者間呈現(xiàn)正相關(guān)關(guān)系,(R2=0.248,P<0.05)(圖2),群體遺傳分布符合距離隔離模型,相關(guān)性為15.75%。

    圖2遺傳距離與地理距離間的相關(guān)性分析

    Fig.2Analysisontherelationshipbetweengeneticdistanceandgeographicdistance

    分子變異等級分析的結(jié)果顯示,松滋釘螺的變異主要存在于個體間,占總變異的86.48%,群體間變異和群體內(nèi)個體間變異僅占總變異的7.02%和6.51%(表7)。群體的Fst為0.144,Nm為1.4857,顯示群體間遺傳差異小,分化不顯著。

    2.3聚類分析 基于10個群體間的遺傳距離,用非加權(quán)平均法(UPGMA)聚類顯示,肋殼(深肋)釘螺SH、SI、SK與光殼螺SE聚為一支后,再與光殼螺SA、SB形成的分支聚為一支;淺肋殼螺SF、SG先聚一起后,再與光殼SD、SC聚為第二分支。

    圖310個湖北釘螺群體的UPGMA聚類分析

    Fig.3UPGMAclusteranalysisoftenO.hupensispopulations

    表7 釘螺種群分子變異分析結(jié)果

    3 討 論

    湖北松滋地區(qū)位于長江中下游地區(qū),該地處巫山山系荊門分支余脈和武陵山系石門分支余脈向江漢平原延伸的過渡地帶,地勢西高東低,是山丘型和平原溝渠交界地帶,生態(tài)景觀多樣,不同景觀下湖北釘螺表現(xiàn)為不同的螺殼形態(tài),即光殼和肋殼(淺肋和深肋)[16]。盡管如此,線粒體基因COI序列分析表明該地不同螺殼形態(tài)的湖北釘螺,并未表現(xiàn)明顯的遺傳差異[17]。由于單個線粒體基因序列標記的局限性和不能對全基因組遺傳多態(tài)性進行分析[18],因此,有必要對該地釘螺的遺傳多樣性進行進一步深入分析。本研究調(diào)查了山丘型和垸內(nèi)溝渠的10個地理群體,螺殼形態(tài)分別表現(xiàn)為光殼和肋殼,其中地處2種環(huán)境交界地帶的新江口鎮(zhèn)同時存在光殼和肋殼釘螺,肋殼則表現(xiàn)為過渡類型的淺肋型。為進一步探討這種螺殼形態(tài)演變趨勢的遺傳本質(zhì),本研究選擇6個微衛(wèi)星DNA位點對該地區(qū)10個湖北釘螺群體進行遺傳多樣性分析,探討松滋地區(qū)不同螺殼形態(tài)的湖北釘螺群體遺傳分化情況。

    本研究6個微衛(wèi)星位點在10個群體中共檢測到141個等位基因,平均每個位點檢測到23.5個,表明各位點在不同群體間多態(tài)性較為明顯,并據(jù)此計算各群體的多態(tài)信息含量(PIC)。PIC是描述群體在微衛(wèi)星座位上的遺傳變異程度的重要指標,PIC 值越大,說明群體內(nèi)基因一致性差,變異越大,選擇潛力也大;反之,PIC 越小,說明群體變異性越小,選擇潛力也小[19]。本研究的10個群體中,PIC值為肋殼螺群體SH、SI、Sk大于光殼螺群體SA、SB、SC、SD、SE,光殼螺群體大于淺殼螺群體SF、SG,提示地處垸內(nèi)溝渠環(huán)境下的肋殼釘螺群體的遺傳變異較大,光殼釘螺群體遺傳變異次之,而淺肋型釘螺群體遺傳變異最小。盡管如此,10個群體的期望雜合度均大于0.5,表明該地區(qū)釘螺群體遺傳分化明顯,具有豐富的遺傳多樣性。

    利用F統(tǒng)計量評估該地區(qū)遺傳多樣性表明,各位點的的FST平均值為0.1440,即14.4%的遺傳變異來自群體間,85.6%的遺傳變異來自群體內(nèi)個體間,這表明松滋地區(qū)湖北釘螺雖然受到地理環(huán)境影響而分割成一些亞群,但各群體間遺傳分化主要是由于釘螺個體間的遺傳變異形成,而Fis、Fit值均為正值,且群體近交系數(shù)為30.32%,也表明釘螺群體間基因流不高[20]。AMOVA分析結(jié)果顯示,10地釘螺個體間遺傳變異(86.48%)遠遠高于群體間遺傳變異(7.02%)和群體內(nèi)遺傳變異(6.51%),也說明遺傳變異主要來自群體個體間,群體間遺傳分化很小,其原因是由于采集點分布于同一水系,遺傳分化并不明顯[20]。因此,盡管在遺傳距離與地理距離表現(xiàn)出一定的相關(guān)性,但由于本研究為小尺度下景觀群體,因此,群體間遺傳分化指數(shù)相差并不大,這可能是由于10個景觀環(huán)境距離并不遠,雖然殼型不同,但是都屬于湖沼型群體的原因。利用UPGMA法聚類樹顯示,肋殼中的淺肋和深肋型釘螺群體各自形成明顯的聚類后,分別與不同的光殼釘螺群體形成聚類,表明光殼釘螺群體中遺傳多樣性更為豐富。

    綜上所述,盡管松滋地區(qū)湖北釘螺螺殼形態(tài)表現(xiàn)出多種形態(tài),但是微衛(wèi)星DNA掃描各釘螺群體顯示,雖然群體內(nèi)遺傳多樣性較高,但遺傳變異主要來自于群體個體間,各景觀環(huán)境的釘螺群體間遺傳分化不明顯,這與石朝輝等人對廟河地區(qū)釘螺不同殼型間基因差異[17]、邱持平等對安徽地區(qū)釘螺不同生態(tài)環(huán)境下CoI基因差異的研究結(jié)果相一致[21]。目前推測可能是由于不同的生態(tài)環(huán)境的影響,如溫度、濕度、海拔等,造成殼型有所不同。

    參考文獻:

    [1]Wang L, Utzinger J, Zhou XN. Schistosomiasis control: experiences and lessons from China[J]. Lancet, 2008, 372(9652): 1793-1795. DOI: 10.1016/S0140-6736(08)61358-6

    [2]Utzinger J, Zhou XN, Chen MG, et al. Conquering schistosomiasis in China: the long march[J]. Acta Trop, 2005, 96(2-3): 69-96. DOI: 10.1016/j.actatropica.2005.08.004

    [3]Li SZ, Wang YX, Yang K, et al. Landscape genetics: the correlation of spatial and genetic distances ofOncomelaniahupensis, the intermediate host snail ofSchistosomajaponicumin mainland China[J]. Geospatial Health, 2009, 3(2): 221-231.

    [4]Shi CH, Wilke T, Davis GM, et al. Population genetics, micro-phylogeography, ecology, and susceptibility to schistosome infection of Chinese Oncomelania hupensis hupensis (Gastropoda: Rissooidea: Pomatiopsidae) in the Miao River system[J]. Malacologia, 2002, 44(2): 333-347.

    [5]Sun XW, Zhang XF, Zhao YY, et al. Development and application of microsatellite markers in aquatic species[J]. J Fishery Sci China, 2008, 15(4): 689-703. (in Chinese)

    孫效文, 張曉鋒, 趙瑩瑩, 等.水產(chǎn)生物微衛(wèi)星標記技術(shù)研究進展及其應用[J]. 中國水產(chǎn)科學, 2008, 15(4): 689-703.

    [6]Guo W, Shen ZR. Microsatellite marker and its application in entomology[J]. Biotechnology, 2004,14(2): 60-61. (in Chinese)

    國偉,沈佐銳. 微衛(wèi)星DNA標記技術(shù)及其在昆蟲學上的應用[J]. 生物技術(shù),2004,14(2):60-61.

    [7]Niu AO, Xiong YW. Studies on the genetic variation ofOncomelaniahupensiswith SSR-PCR[J]. J Parasit Biol, 2002, 15(4): 230-233. (in Chinese)

    牛安歐,熊衍文. 微衛(wèi)星錨定PCR研究湖北釘螺的遺傳變異[J]. 中國寄生蟲病防治雜志,2002,15(4):230-233.

    [8]Zhou YB, Zhao GM, Wei JG, et al. Genetic diversity in 19 Chinese populations ofOncomelaniahupensis(Gastropoda: Rissooidea) detected by simple sequence repeat-anchored polymerase chain reaction amplification[J]. Chin J Epidemiol, 2007, 28(9): 859-862. (in Chinese)

    周藝彪,趙根明,韋建國,等. 微衛(wèi)星錨定PCR分析19個湖北釘螺種群之間的遺傳變異關(guān)系[J]. 中華流行病學雜志,2007,28(9):859-862.

    [9]Guo JT, Zhou YB, Zhang ZJ, et al. Simple sequence repeat variation and small-scale spatial autocorrelatian analysis on smooth-shell populations ofOncomelaniahupensisin Sichuan province[J]. Chin J Epidemiol, 2009, 30 (5): 497-501. (in Chinese)

    郭俊濤, 周藝彪, 張志杰, 等. 四川光殼釘螺微衛(wèi)星遺傳變異的小尺度空間自相關(guān)分析[J]. 中華流行病學雜志, 2009, 30 (5): 497-501.

    [10]Ma L, Li SZ, Yang P, et al. Comparison of five different methods of extracting genomic DNA fromOncomelaniahupensis[J]. J Pathog Biol, 2011, 6(2): 129-132. (in Chinese)

    馬琳,李石柱,楊頻,等. 湖北釘螺基因組5種抽提方法的比較[J]. 中國病原生物學雜志,2011,6(2):129-132.

    [11]Li SZ, Wang YX, Ma YJ, et al. Isolation and characterization of polymorphic microsatellite markers ofOncomelaniahupensis[J]. Chin J Schistosomiasis Ctrl, 2010, 22(2): 122-125. (in Chinese)

    李石柱,王藝秀,馬雅軍,等. 湖北釘螺多態(tài)微衛(wèi)星DNA位點篩選和特征初步分析[J].中國血吸蟲病防治雜志,2010,22(2):122-125.

    [12]Zhang L, Li SZ, Wang Q, et al. Isolation and characterization of 15 new microsatellite markers inOncomelaniahupensis, the snail intermediate host ofSchistosomajaponicumin Mainland China[J]. Int J Mol Sci, 2012, 13(5): 5844-5850. DOI: 10.3390/ijms13055844

    [13]Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): an integrated software package for population genetics data analysis[J]. Evol Bioinformat Online, 2005, 1: 47-50.

    [14]Stefan Van Dongen T, Winnepenninckx B. Multiple UPGMA and neighbor-joining trees and the performance of some computer packages[J]. Mol Biol Evol, 1996, 13(2): 309-313.

    [15]Khan HA, Arif IA, Bahkali AH, et al. Bayesian, maximum parsimony and UPGMA models for inferring the phylogenies of antelopes using mitochondrial markers[J]. Evol Bioinformat Online, 2008, 4: 263-270.

    [16]Yang JY. Study on technical characteristics and regional applicability of rice Green control in Songzi City[J]. Hubei Plant Protect, 2010,(5):7-9. (in Chinese)

    楊繼云.松滋市水稻綠色防控技術(shù)特點與區(qū)域適用性研究[J]. 湖北植保,2010, (5): 7-9.

    [17]Shi CH, Qiu CP, Xia MY, et al. Preliminary study on cytochrome C oxidase 1 gene ofOncomelaniahupensisfrom Miao River Area in Hubei Province[J]. Chin J Parasitol Parasit Dis, 2010, 19(1): 41-44. (in Chinese)

    石朝輝, 邱持平, 夏明儀, 等. 湖北省廟河地區(qū)釘螺細胞色素C氧化酶1基因差異的研究[J].中國寄生蟲學與寄生蟲病雜志, 2010, 19(1):41-44.

    [18]Hu Y, Li XM, Lin R, et al. Genetic variation of two mitochondrial DNA molecules from three isolates ofOncomelaniahupensis[J]. Chin J Parasitol Parasit Dis, 2007, 25(6): 474-477. (in Chinese)

    胡纓, 黎學銘, 林睿, 等. 三地釘螺線粒體 DNA 兩個分子的遺傳變異研究[J]. 中國寄生蟲學與寄生蟲病雜志, 2007, 25(6):474-477.

    [19]Botstein D, White RL, Skolnick M, et al. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. Am J Hum Genet, 1980, 32(3): 314-331.

    [20]Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure[J]. Evolution, 1984: 1358-1370.

    [21]Qiu CP, Chris S, Xia MY, et al. Study on genomic diversity among different isolates ofSchistosomajaponicumin China[J]. Chin J Parasitol Parasit Dis, 2002, 20(2): 79-82. (in Chinese)

    邱持平, Chris S, 夏明儀, 等. 中國日本血吸蟲地域株基因差異的研究[J]. 中國寄生蟲學與寄生蟲病雜志, 2002, 20(2):79-82.

    猜你喜歡
    螺殼釘螺遺傳變異
    2020年洪澇災害后湖北省釘螺擴散情況調(diào)查
    我不要做垃圾蟹!
    2016—2021年池州市貴池區(qū)新發(fā)和復現(xiàn)釘螺情況分析
    為什么寄居蟹居住在螺殼里
    中外文摘(2021年14期)2021-11-11 18:42:17
    一種快速鑒別湖北釘螺(Oncomelania hupensis)雌雄性的方法
    先導編輯技術(shù)可編輯近90%的人類遺傳變異
    基于改進遺傳變異算子的海島算法
    電子制作(2019年24期)2019-02-23 13:22:18
    火力楠子代遺傳變異分析及優(yōu)良家系選擇
    GABABR2基因遺傳變異與肥胖及代謝相關(guān)表型的關(guān)系
    螺殼漂流記
    亚洲人成网站在线观看播放| 国产黄色小视频在线观看| 国产熟女欧美一区二区| 亚洲人成网站在线观看播放| 日韩一区二区三区影片| 成人午夜高清在线视频| 在线免费观看的www视频| 91午夜精品亚洲一区二区三区| 人妻制服诱惑在线中文字幕| 成人亚洲精品av一区二区| 水蜜桃什么品种好| 嫩草影院精品99| 精品久久久久久久久亚洲| 日产精品乱码卡一卡2卡三| 免费大片黄手机在线观看| 免费少妇av软件| 晚上一个人看的免费电影| 91在线精品国自产拍蜜月| 男女边吃奶边做爰视频| xxx大片免费视频| 两个人的视频大全免费| 国产亚洲精品av在线| 中国美白少妇内射xxxbb| 久久人人爽人人片av| 国产精品久久久久久精品电影小说 | 天天躁夜夜躁狠狠久久av| 日日摸夜夜添夜夜添av毛片| 伊人久久国产一区二区| 97精品久久久久久久久久精品| 国产免费视频播放在线视频 | 国产单亲对白刺激| 18禁裸乳无遮挡免费网站照片| av在线天堂中文字幕| 婷婷六月久久综合丁香| 日韩精品青青久久久久久| 欧美区成人在线视频| 免费av观看视频| 亚洲人成网站高清观看| 国产精品熟女久久久久浪| 成人漫画全彩无遮挡| 欧美另类一区| 午夜激情久久久久久久| 丰满人妻一区二区三区视频av| 日韩av不卡免费在线播放| 欧美日韩一区二区视频在线观看视频在线 | 免费看a级黄色片| 伊人久久精品亚洲午夜| 真实男女啪啪啪动态图| 人妻系列 视频| 六月丁香七月| 91久久精品国产一区二区三区| 精品久久久久久电影网| 日韩中字成人| 国产精品国产三级专区第一集| 我的老师免费观看完整版| 日韩电影二区| 亚洲av电影不卡..在线观看| 韩国高清视频一区二区三区| 毛片女人毛片| 国产成人a∨麻豆精品| 在线观看美女被高潮喷水网站| 国产一区亚洲一区在线观看| 成人国产麻豆网| 亚洲精品日韩av片在线观看| 纵有疾风起免费观看全集完整版 | 亚洲国产精品sss在线观看| av在线观看视频网站免费| 你懂的网址亚洲精品在线观看| 精品酒店卫生间| 国产成人freesex在线| 国产成人一区二区在线| 在线播放无遮挡| a级一级毛片免费在线观看| 亚洲天堂国产精品一区在线| 亚洲av不卡在线观看| 麻豆精品久久久久久蜜桃| 少妇熟女欧美另类| 亚洲成人久久爱视频| 春色校园在线视频观看| 亚洲高清免费不卡视频| 久久综合国产亚洲精品| 日韩强制内射视频| 日韩成人av中文字幕在线观看| 日韩欧美三级三区| 搡女人真爽免费视频火全软件| 精品国产三级普通话版| 在线观看一区二区三区| 精品一区二区免费观看| 亚洲国产精品sss在线观看| 蜜臀久久99精品久久宅男| 日韩av在线免费看完整版不卡| 日本免费在线观看一区| 精品国产露脸久久av麻豆 | 肉色欧美久久久久久久蜜桃 | 精品熟女少妇av免费看| 成年版毛片免费区| 久久久久久久久久成人| 国产又色又爽无遮挡免| 精品久久久久久久人妻蜜臀av| 十八禁国产超污无遮挡网站| 一级爰片在线观看| 小蜜桃在线观看免费完整版高清| 夜夜爽夜夜爽视频| 精品久久久噜噜| 成人美女网站在线观看视频| 亚洲在线观看片| 小蜜桃在线观看免费完整版高清| 免费av观看视频| 久久久久免费精品人妻一区二区| 欧美日本视频| 国产精品人妻久久久久久| 国产在线男女| 国产高清三级在线| 精品亚洲乱码少妇综合久久| 久久人人爽人人片av| 禁无遮挡网站| 国内精品宾馆在线| 好男人视频免费观看在线| 午夜福利在线观看吧| 日韩欧美精品免费久久| 成人av在线播放网站| 能在线免费观看的黄片| 美女高潮的动态| 免费在线观看成人毛片| 亚洲成人久久爱视频| 久99久视频精品免费| 亚洲国产精品成人久久小说| 18禁裸乳无遮挡免费网站照片| 最后的刺客免费高清国语| 国产av国产精品国产| 最近中文字幕高清免费大全6| 亚洲精品久久午夜乱码| 少妇裸体淫交视频免费看高清| 亚洲av一区综合| 观看美女的网站| 在线观看美女被高潮喷水网站| 深爱激情五月婷婷| av在线天堂中文字幕| 亚洲av一区综合| 观看美女的网站| 97热精品久久久久久| 美女cb高潮喷水在线观看| 精品久久久精品久久久| 午夜精品在线福利| 国产片特级美女逼逼视频| 午夜免费激情av| 能在线免费观看的黄片| 午夜激情福利司机影院| 久99久视频精品免费| 亚洲成人久久爱视频| 女的被弄到高潮叫床怎么办| 欧美日本视频| 亚洲丝袜综合中文字幕| 联通29元200g的流量卡| 国产午夜精品一二区理论片| 亚洲三级黄色毛片| 80岁老熟妇乱子伦牲交| 男女那种视频在线观看| 欧美潮喷喷水| 青春草视频在线免费观看| 亚洲,欧美,日韩| 久久久久久久久久人人人人人人| 日韩欧美 国产精品| 在线观看av片永久免费下载| 自拍偷自拍亚洲精品老妇| 少妇人妻一区二区三区视频| 婷婷色综合大香蕉| 在线观看一区二区三区| 国产精品av视频在线免费观看| 久久久久精品久久久久真实原创| 日韩精品有码人妻一区| 色综合亚洲欧美另类图片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美 日韩 精品 国产| 免费黄网站久久成人精品| 午夜久久久久精精品| 99热这里只有是精品50| 国产精品99久久久久久久久| 免费看不卡的av| 国产伦理片在线播放av一区| 91精品一卡2卡3卡4卡| 精品久久久久久久人妻蜜臀av| 国内精品一区二区在线观看| 亚洲欧美一区二区三区黑人 | 成人性生交大片免费视频hd| 亚洲精品日本国产第一区| 亚洲人与动物交配视频| 午夜精品国产一区二区电影 | 色综合站精品国产| 婷婷色综合www| 国产美女午夜福利| 久久久久久久国产电影| 国产色婷婷99| 黄色日韩在线| 大话2 男鬼变身卡| 日韩人妻高清精品专区| 国产午夜精品久久久久久一区二区三区| 日韩伦理黄色片| av又黄又爽大尺度在线免费看| 久久久精品94久久精品| 超碰av人人做人人爽久久| 毛片一级片免费看久久久久| 毛片女人毛片| 国产精品国产三级国产专区5o| 国产精品蜜桃在线观看| 日韩一区二区视频免费看| 不卡视频在线观看欧美| 网址你懂的国产日韩在线| av.在线天堂| 黄片wwwwww| 日本免费a在线| 日本一二三区视频观看| 超碰97精品在线观看| av黄色大香蕉| 国产成人精品一,二区| 十八禁网站网址无遮挡 | 国产伦在线观看视频一区| 国产精品一区二区在线观看99 | av在线观看视频网站免费| 免费看a级黄色片| 精品一区二区免费观看| 97在线视频观看| 久久99热这里只频精品6学生| 欧美日韩国产mv在线观看视频 | 国产精品.久久久| 欧美xxⅹ黑人| 人人妻人人澡欧美一区二区| 精品午夜福利在线看| 亚洲av电影在线观看一区二区三区 | 日本免费在线观看一区| 边亲边吃奶的免费视频| 深夜a级毛片| 男女国产视频网站| 大香蕉97超碰在线| 日本黄色片子视频| 国产一区二区三区综合在线观看 | 白带黄色成豆腐渣| 日韩一本色道免费dvd| av在线观看视频网站免费| 又大又黄又爽视频免费| 日韩av免费高清视频| av在线播放精品| 国产成人福利小说| 久久久精品94久久精品| 国模一区二区三区四区视频| 99热6这里只有精品| 亚洲av在线观看美女高潮| 人人妻人人看人人澡| 亚洲色图av天堂| 国产乱人视频| 亚洲欧美成人综合另类久久久| 少妇的逼好多水| 日日摸夜夜添夜夜爱| 国产精品美女特级片免费视频播放器| 精品久久久久久成人av| 99久久精品一区二区三区| 欧美不卡视频在线免费观看| 国内少妇人妻偷人精品xxx网站| 男人舔女人下体高潮全视频| 亚洲内射少妇av| 天天躁日日操中文字幕| 少妇人妻精品综合一区二区| 日本三级黄在线观看| 波多野结衣巨乳人妻| 99视频精品全部免费 在线| 人人妻人人澡欧美一区二区| 黑人高潮一二区| 国产精品久久久久久av不卡| 国产乱来视频区| 国产淫语在线视频| 日本免费在线观看一区| 最新中文字幕久久久久| 69人妻影院| 三级毛片av免费| 免费人成在线观看视频色| 成人毛片60女人毛片免费| 亚洲av.av天堂| 久久综合国产亚洲精品| 最近2019中文字幕mv第一页| 日韩欧美国产在线观看| 欧美激情久久久久久爽电影| 亚洲精品久久久久久婷婷小说| 我的女老师完整版在线观看| 久久6这里有精品| 国产精品一二三区在线看| 日韩成人伦理影院| 91午夜精品亚洲一区二区三区| or卡值多少钱| 午夜激情福利司机影院| a级一级毛片免费在线观看| 三级国产精品欧美在线观看| 色综合亚洲欧美另类图片| 一区二区三区四区激情视频| 三级毛片av免费| 少妇的逼水好多| 亚洲欧美精品专区久久| 青青草视频在线视频观看| 亚洲成人一二三区av| 国产伦在线观看视频一区| 国产黄频视频在线观看| 日本猛色少妇xxxxx猛交久久| 国产精品熟女久久久久浪| 国产亚洲av片在线观看秒播厂 | 国产探花极品一区二区| 成人毛片60女人毛片免费| 精品久久国产蜜桃| 精品欧美国产一区二区三| 天堂网av新在线| 国产一级毛片七仙女欲春2| 97超视频在线观看视频| 亚洲高清免费不卡视频| 久久久久久久国产电影| 熟妇人妻久久中文字幕3abv| 免费高清在线观看视频在线观看| 久久久久久国产a免费观看| 好男人在线观看高清免费视频| 夜夜爽夜夜爽视频| ponron亚洲| 国产 亚洲一区二区三区 | 日本av手机在线免费观看| 国产精品国产三级国产专区5o| kizo精华| freevideosex欧美| 国产亚洲最大av| 天堂俺去俺来也www色官网 | 亚洲欧美一区二区三区黑人 | 少妇的逼水好多| 久久精品人妻少妇| 一级片'在线观看视频| 五月天丁香电影| 大又大粗又爽又黄少妇毛片口| 三级国产精品欧美在线观看| 国产 亚洲一区二区三区 | 一个人观看的视频www高清免费观看| 超碰97精品在线观看| 国产伦一二天堂av在线观看| 国产伦精品一区二区三区视频9| 国产精品爽爽va在线观看网站| 免费在线观看成人毛片| 尾随美女入室| 国产一区二区在线观看日韩| 插阴视频在线观看视频| 欧美一级a爱片免费观看看| 网址你懂的国产日韩在线| 久久精品夜夜夜夜夜久久蜜豆| 久久99热这里只频精品6学生| 免费看av在线观看网站| av在线亚洲专区| 国产av在哪里看| 亚洲av电影在线观看一区二区三区 | 日韩精品有码人妻一区| 久久99热这里只频精品6学生| 女人被狂操c到高潮| 中文资源天堂在线| 神马国产精品三级电影在线观看| 91久久精品国产一区二区三区| 国产美女午夜福利| 人人妻人人看人人澡| 亚洲精品亚洲一区二区| www.色视频.com| 日韩av在线大香蕉| 欧美成人精品欧美一级黄| 久久久久久久久久成人| 亚洲欧美日韩无卡精品| 免费观看性生交大片5| 又爽又黄无遮挡网站| 伊人久久精品亚洲午夜| 久久精品国产自在天天线| 欧美日韩国产mv在线观看视频 | 成人国产麻豆网| 欧美日韩国产mv在线观看视频 | 色综合亚洲欧美另类图片| 亚洲欧美清纯卡通| 日韩电影二区| 91午夜精品亚洲一区二区三区| 亚洲av成人av| av女优亚洲男人天堂| 啦啦啦韩国在线观看视频| 天天躁夜夜躁狠狠久久av| 乱人视频在线观看| 欧美丝袜亚洲另类| av在线亚洲专区| 美女国产视频在线观看| 国产精品久久久久久精品电影小说 | 高清午夜精品一区二区三区| 亚洲精品国产成人久久av| 韩国av在线不卡| 亚洲性久久影院| 久久精品国产自在天天线| 午夜福利在线观看吧| 国产高清三级在线| 免费看a级黄色片| 久久久国产一区二区| 国产精品蜜桃在线观看| 久久久久久久久久久免费av| 久久久成人免费电影| 午夜精品在线福利| 成人午夜精彩视频在线观看| 亚洲精品,欧美精品| 两个人视频免费观看高清| 久久精品国产亚洲av天美| 99视频精品全部免费 在线| 国产一区亚洲一区在线观看| 国产成人aa在线观看| 色5月婷婷丁香| 免费播放大片免费观看视频在线观看| 嫩草影院入口| 99热这里只有是精品50| 久久久久久九九精品二区国产| 淫秽高清视频在线观看| 亚洲精品成人av观看孕妇| 免费看光身美女| 99久久精品热视频| 最近2019中文字幕mv第一页| 麻豆成人av视频| 久久久久久伊人网av| 99热6这里只有精品| 色网站视频免费| 91久久精品国产一区二区三区| 一级av片app| 全区人妻精品视频| 99九九线精品视频在线观看视频| 久久精品国产亚洲av天美| 久久精品人妻少妇| 国产精品久久视频播放| 联通29元200g的流量卡| 91aial.com中文字幕在线观看| 色播亚洲综合网| 久久久久久久久久久丰满| 国产成人精品久久久久久| 国产亚洲午夜精品一区二区久久 | 少妇猛男粗大的猛烈进出视频 | 久久久久久久国产电影| 国产精品一区二区性色av| 91在线精品国自产拍蜜月| 又黄又爽又刺激的免费视频.| 人人妻人人澡欧美一区二区| 午夜福利视频1000在线观看| 中文字幕久久专区| 国产爱豆传媒在线观看| 一夜夜www| 又黄又爽又刺激的免费视频.| 国产精品国产三级专区第一集| 国产精品久久久久久av不卡| 2021天堂中文幕一二区在线观| 超碰av人人做人人爽久久| 99久久人妻综合| 国产爱豆传媒在线观看| 成人无遮挡网站| 精品国内亚洲2022精品成人| 亚洲欧洲国产日韩| 能在线免费观看的黄片| 国产白丝娇喘喷水9色精品| 床上黄色一级片| 久久久久久久国产电影| 大话2 男鬼变身卡| 69av精品久久久久久| 99久久精品一区二区三区| 国产一区二区在线观看日韩| 最近手机中文字幕大全| 国产精品人妻久久久久久| 日韩大片免费观看网站| 女的被弄到高潮叫床怎么办| 久久热精品热| 久久午夜福利片| 欧美日本视频| 中文精品一卡2卡3卡4更新| 久久久久久久大尺度免费视频| 国产熟女欧美一区二区| 狂野欧美激情性xxxx在线观看| 亚洲人成网站在线播| 男女啪啪激烈高潮av片| 国产在视频线精品| 亚洲在久久综合| av免费在线看不卡| 美女高潮的动态| 亚洲激情五月婷婷啪啪| 久久久色成人| www.av在线官网国产| 床上黄色一级片| 国产v大片淫在线免费观看| 永久免费av网站大全| 国模一区二区三区四区视频| 国产伦理片在线播放av一区| 亚洲欧美中文字幕日韩二区| 国内精品一区二区在线观看| 国产成人福利小说| 九九爱精品视频在线观看| 久久热精品热| 亚洲高清免费不卡视频| 在线播放无遮挡| 爱豆传媒免费全集在线观看| 久久精品久久久久久久性| 欧美日韩一区二区视频在线观看视频在线 | 非洲黑人性xxxx精品又粗又长| 国产黄片美女视频| 嫩草影院精品99| av在线老鸭窝| 欧美日韩精品成人综合77777| 亚洲国产精品成人久久小说| 亚洲av二区三区四区| 中国国产av一级| 久久精品国产亚洲av天美| 亚洲精品日韩在线中文字幕| 午夜精品在线福利| 不卡视频在线观看欧美| 国产成人福利小说| 国产久久久一区二区三区| 国产真实伦视频高清在线观看| 国产精品一区www在线观看| 精品不卡国产一区二区三区| 成人午夜高清在线视频| 日韩电影二区| 色吧在线观看| 如何舔出高潮| 午夜福利在线在线| 成年人午夜在线观看视频 | 亚洲国产精品sss在线观看| 人妻少妇偷人精品九色| 在线a可以看的网站| or卡值多少钱| 免费人成在线观看视频色| 亚州av有码| videos熟女内射| 亚洲欧美成人综合另类久久久| 色视频www国产| 两个人的视频大全免费| 两个人视频免费观看高清| 天堂av国产一区二区熟女人妻| 久久久久久久久久黄片| 国产在视频线在精品| 日本-黄色视频高清免费观看| 99久久精品国产国产毛片| 久久久久国产网址| 免费看美女性在线毛片视频| 日韩成人伦理影院| 国产高清国产精品国产三级 | 亚洲自偷自拍三级| 亚洲av不卡在线观看| 亚洲真实伦在线观看| 大香蕉97超碰在线| 国产午夜福利久久久久久| 午夜老司机福利剧场| www.av在线官网国产| 99九九线精品视频在线观看视频| 国产欧美日韩精品一区二区| 日韩三级伦理在线观看| 伊人久久国产一区二区| 日韩伦理黄色片| 嘟嘟电影网在线观看| 国产三级在线视频| 国产中年淑女户外野战色| 中文天堂在线官网| 亚洲最大成人手机在线| 欧美日韩综合久久久久久| 午夜福利视频精品| 亚洲欧美日韩东京热| 免费大片18禁| kizo精华| 激情 狠狠 欧美| 波多野结衣巨乳人妻| 日本猛色少妇xxxxx猛交久久| 一级爰片在线观看| 亚洲av.av天堂| 伊人久久国产一区二区| 26uuu在线亚洲综合色| 国产精品伦人一区二区| 汤姆久久久久久久影院中文字幕 | 别揉我奶头 嗯啊视频| 亚洲欧美成人综合另类久久久| 久热久热在线精品观看| 中文欧美无线码| 久热久热在线精品观看| 婷婷色综合大香蕉| 久久久久久久大尺度免费视频| 欧美日韩精品成人综合77777| 特大巨黑吊av在线直播| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av电影不卡..在线观看| 热99在线观看视频| 色尼玛亚洲综合影院| 婷婷六月久久综合丁香| 天堂影院成人在线观看| 日韩大片免费观看网站| 又黄又爽又刺激的免费视频.| 亚洲av免费在线观看| 99热这里只有是精品50| 亚洲av福利一区| 国产成人a区在线观看| 97超视频在线观看视频| 在线观看人妻少妇| 天堂√8在线中文| 日韩电影二区| 综合色av麻豆| 国产成人一区二区在线| 最近中文字幕高清免费大全6| 久久鲁丝午夜福利片| 卡戴珊不雅视频在线播放| 精品国产露脸久久av麻豆 | 亚洲av电影不卡..在线观看| 亚洲怡红院男人天堂| 最近最新中文字幕大全电影3| 亚洲av成人av| 99久久精品一区二区三区| 男女边吃奶边做爰视频| a级一级毛片免费在线观看| 亚洲精品国产成人久久av| 特级一级黄色大片| 啦啦啦中文免费视频观看日本| 伦理电影大哥的女人| 成人无遮挡网站| 精品国产露脸久久av麻豆 | 亚洲av在线观看美女高潮| 亚洲精品乱久久久久久| 青春草视频在线免费观看|