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

    中國虎紋蛙多堿基重復(fù)微衛(wèi)星位點的多態(tài)性

    2024-04-30 00:37:50喬芬邵偉偉馬力林植華韋力
    安徽農(nóng)業(yè)科學(xué) 2024年7期
    關(guān)鍵詞:分子標記遺傳多樣性多態(tài)性

    喬芬 邵偉偉 馬力 林植華 韋力

    摘要 ?[目的]篩選出中國虎紋蛙高質(zhì)量多態(tài)性的微衛(wèi)星標記,用于該物種的鑒定以及遺傳多樣性、種群結(jié)構(gòu)、遺傳連鎖圖譜分析。[方法]利用新一代測序技術(shù),從中國虎紋蛙部分基因組中篩選出33個三核苷酸重復(fù)、42個四核苷酸重復(fù)和1個五核苷酸重復(fù)微衛(wèi)星位點,共76個多堿基重復(fù)微衛(wèi)星位點。所有位點對浙江麗水學(xué)院兩棲爬行動物實驗室人工飼養(yǎng)的種群(n=30)進行了擴增和基因型分析。[結(jié)果]每個位點觀察到的等位基因數(shù)和多態(tài)信息含量(PIC)的平均值分別為4個(2~7個)和0.514(0.164~0.822),觀測雜合度(HO)和期望雜合度(HE)分別為0.00~0.97和0.18~0.82。有2個位點(THC141和THC377)顯著偏離哈迪溫伯格平衡,在11個位點中發(fā)現(xiàn)了顯著的連鎖不平衡。[結(jié)論]這些新篩選的多態(tài)性微衛(wèi)星標記有助于進一步研究中國虎紋蛙的遺傳多樣性和開展遺傳育種,為制定有效的保護策略和開展分子輔助選擇提供依據(jù)。

    關(guān)鍵詞? 中國虎紋蛙;多堿基重復(fù)微衛(wèi)星;遺傳多樣性;分子標記;多態(tài)性

    中圖分類號? S917.4?? 文獻標識碼? A?? 文章編號? 0517-6611(2024)07-0082-05

    doi: 10.3969/j.issn.0517-6611.2024.07.020

    Polymorphic Multiple-base Repeat Microsatellite Loci of Chinese Tiger Frog Hoplobatrachus chinensis

    QIAO Fen, SHAO Wei-wei, MA Li et al

    (College of Ecology, Lishui University, Lishui, Zhejiang 323000)

    Abstract? [Objective] Screening of highquality polymorphism microsatellite markers in the Chinese tiger frog Hoplobatrachus chinensis can be used for identification, genetic diversity, population structure, and genetic linkage map analysis of this species. [Method] In this study, we identified 76 multiplebase repeat microsatellite loci, including 33 trinucleotide, 42 tetranucleotide and 1 pentanucleotide repeat microsatellite loci from partial genome of Chinese tiger frog Hoplobatrachus chinensis using the nextgeneration sequencing technology. We genotyped 30 individuals from a population in an artificial pond at the herpetological laboratory of Lishui University, Zhejiang Province, China. [Result] The results showed that the mean number of observed alleles and polymorphic information content (PIC) per locus was 4 (range from 2 to 7) and 0.514 (range from 0.164 to 0.822), respectively. The observed and expected heterozygosity values ranged from 0.00 to 0.97 and from 0.18 to 0.82, respectively. Two loci (THC141 and THC377) showed significant departure from HardyWeinberg equilibrium and significant linkage disequilibrium was found among 11 pairwise loci. [Conclusion] These novel polymorphic microsatellite markers are expected to be valuable for future studies on genetic diversity and conducting genetic breeding, enabling the formulation of effective conservation strategies and the use of molecular assisted selection for H. chinensis.

    Key words? Hoplobatrachus chinensis;Multiplebase repeat microsatellite;Genetic diversity;Molecular marker;Polymorphism

    中國虎紋蛙(Hoplobatrachus chinensis)之前被稱為H.rugulosus,屬于無尾目(Anura)蛙科(Ranidae),主要棲息于農(nóng)田及其附近區(qū)域[1]。與其他兩棲動物一樣,中國虎紋蛙在害蟲防治、水生和陸生食物網(wǎng)等方面發(fā)揮著重要作用。關(guān)于中國虎紋蛙的研究主要包括生長發(fā)育[2-3]、聲通訊[4]、繁殖特征[5-6]、人工養(yǎng)殖[7]等。作為可食用蛙類,中國虎紋蛙在過去幾十年因過度捕捉而造成野生種群資源急劇減少,目前已被列為國家二級保護動物,因此,許多農(nóng)民在我國南方地區(qū)開展了中國虎紋蛙的人工養(yǎng)殖[8]。由于許多養(yǎng)殖場近交現(xiàn)象嚴重,中國虎紋蛙種質(zhì)資源多樣性嚴重下降。因此,開展中國虎紋蛙優(yōu)質(zhì)人工育種具有重要意義。其中,微衛(wèi)星標記輔助育種是提高種質(zhì)資源質(zhì)量的可行方法之一[9]。迄今為止,雖然已經(jīng)報道了幾十個中國虎紋蛙的微衛(wèi)星位點[10-11],但其中大多數(shù)位點是二核苷酸重復(fù)類型,且不同等位基因之間通常沒什么差異,影響種群研究中基因分型的準確性[12]。相比之下,應(yīng)用多堿基重復(fù)微衛(wèi)星可以提高基因分型的準確性,從而在群體遺傳結(jié)構(gòu)和保護遺傳學(xué)方面具有更高的效 率[12]。為了進一步研究中國虎紋蛙的保護遺傳學(xué),筆者篩選并鑒定了76個新的多堿基重復(fù)微衛(wèi)星標記,為進一步研究該物種的遺傳多樣性和人工育種提供了有力的工具。

    1 材料與方法

    采用組織提取試劑盒從中國虎紋蛙的肌肉組織中分離得到基因組DNA。利用新一代測序技術(shù),在北京諾禾致源科技股份有限公司進行部分基因組DNA的測序。測序后,利用SSR Hunter軟件(Version 1.3)對已組裝的序列進行搜索和微衛(wèi)星位點識別[13]。選定的微衛(wèi)星序列使用 primer 5.0軟件進行引物設(shè)計。對于設(shè)計的每一對微衛(wèi)星引物,通過生工生物工程(上海)股份有限公司進行引物合成,然后以6個中國虎紋蛙個體的基因組DNA為模板,優(yōu)化擴增條件。

    引物擴增條件優(yōu)化后,對來自麗水學(xué)院兩棲爬行動物實驗室人工養(yǎng)殖的30個中國虎紋蛙DNA樣品進行擴增,擴增產(chǎn)物送往武漢天一輝遠生物科技有限公司進行測序和分型。采用 M13尾端(5′-TGTAAAACGACGGCCAGT-3′)標記引物進行巢式PCR擴增。所有微衛(wèi)星的正向引物的5′端都加入了M13的尾端標記,PCR混合物中加入了熒光標記的M13引物[11]。每個PCR擴增體系的總體積為15.0 μL,包含100 ng 模板DNA 1.0 μL,0.5 μmol/L的正向引物1.0 μL,0.5 μmol/L的反向引物1.0 μL,0.5 mmol/L的dNTP、0.1 U 熱啟動酶(QIAGEN)和1.2 mmol/L 的MgCl2 7.5 μL,不足用滅菌水補足。PCR擴增的條件是:先在95 ℃的初始變性溫度下進行5 min,然后進行34次95 ℃的30 s循環(huán),Ta(最佳退火溫度)30 s,72 ℃延伸30 s,最后進一步擴增72 ℃ 20 min。PCR產(chǎn)物用ABI 3730測序儀(applied biosystems)進行基因分型,并用Genemarker v1.8軟件進行分析。用Genepop 4.0軟件[14]對觀測雜合度(observed heterozygosity,HO)和期望雜合度(expected heterozygosity,HE)進行計算,用Fstat 2.9.3.2軟件[15]對哈迪溫伯格(hardy-weinberg equilibrium,HWE)平衡和連鎖不平衡(linkage disequilibrium,LD)偏差進行檢驗。用Cervus 3.0軟件計算多態(tài)信息含量(polymorphic information content,PIC)[16]。

    2 結(jié)果與分析

    利用新一代測序技術(shù),從部分基因組中獲得了450條潛在多堿基重復(fù)微衛(wèi)星DNA序列,共設(shè)計并檢測了80對微衛(wèi)星引物。其中4對引物為單帶或擴增質(zhì)量差,其余76對引物(包括33個三核苷酸、42個四核苷酸和1個五核苷酸重復(fù)微衛(wèi)星引物)在30個中國虎紋蛙檢測樣品中均呈現(xiàn)多態(tài)性(表1)。共觀察到290個不同的等位基因位點,每個位點平均為4個(范圍為2~7個等位基因)。等位基因大小在99(THC353)到300(THC346)之間。

    微衛(wèi)星的多態(tài)信息含量(PIC)平均值為0.514,范圍在0.164(THC098和THC197)和0.822(THC226)之間。

    觀測雜合度(HO)范圍為0.000~0.967,期望雜合度(HE)范圍為0.184~0.822。經(jīng)過Bonferroni校對后,發(fā)現(xiàn)有2個微衛(wèi)星位點顯著偏離哈迪溫伯格平衡(HWE)。經(jīng)過連鎖不平衡(linkage disequilibrium)檢測試驗結(jié)果,發(fā)現(xiàn)有11對位點間存在連鎖不平衡,即位點THC100與THC105之間、THC100與THC250之間、THC100與THC320之間、THC031與THC105之間、THC031與THC320之間、THC063與THC204之間、THC075與THC281之間、THC105與THC320之間、THC106與THC161之間、THC130與THC352之間以及THC312與THC346之間。

    3 討論

    微衛(wèi)星在原核生物和真核生物中普遍存在[17-18]。由于其高重復(fù)性、多等位基因性、共顯性遺傳模式以及在基因組中的廣泛覆蓋率,微衛(wèi)星已成為最受歡迎的遺傳標記類別之一[19]。傳統(tǒng)的微衛(wèi)星篩選方法包括磁珠富集法、選擇性雜交和表達序列標簽等,這些篩選方法雖有效但數(shù)量有限、耗時長且價格昂貴[20]。近年來,利用計算機軟件從DNA序列數(shù)據(jù)庫中對微衛(wèi)星序列進行篩選,已經(jīng)迅速取代了從基因組文庫中篩選微衛(wèi)星標記的傳統(tǒng)方法[21]。該研究通過新一代測序技術(shù)對中國虎紋蛙進行簡化測序,從其部分基因組中篩選微衛(wèi)星序列。結(jié)果表明,通過簡化測序能夠獲得足夠多的潛在微衛(wèi)星序列。對設(shè)計的80對微衛(wèi)星引物進行檢測,發(fā)現(xiàn)有76對呈現(xiàn)出多態(tài)性,進一步說明基因組篩選微衛(wèi)星的有效性。用76對引物對30個檢測樣品進行擴增和分型,共觀察到290個不同的等位基因位點,每個位點平均為4個,等位基因大小在99到300之間。Shao等[10]利用磁珠富集法獲得9對中國虎紋蛙微衛(wèi)星位點,共觀察到77個不同的等位基因,平均為8.6個等位基因。Chen等[11]通過轉(zhuǎn)錄組測序法獲得33對中國虎紋蛙微衛(wèi)星位點,共觀察到174個不同的等位基因,平均為5.3個等位基因。造成這些差異的原因可能與檢測樣品個體的親緣關(guān)系有關(guān)。

    該研究中,中國虎紋蛙微衛(wèi)星的多態(tài)信息含量(PIC)平均值為0.514,根據(jù)Botstein等[22]描述的PIC值,發(fā)現(xiàn)該研究有43個微衛(wèi)星位點為高信息量位點(PIC>0.5),31個微衛(wèi)星位點為中等信息量位點(0.5>PIC>0.25),2個微衛(wèi)星位點為低信息量位點(PIC<0.25),該結(jié)果與Chen等[11]對該物種的轉(zhuǎn)錄組測序篩選的微衛(wèi)星結(jié)果相似,即PIC為0.552。類似的結(jié)果見海陸蛙(Fejervarya cancrivora)[23]的研究報道。但是,與科基蛙(Eleutherodactylus coqui)[24]的PIC(均值為0743)和飾紋姬蛙(Microhyla fissipes)[25]的PIC(均值為0.679)相比,該研究的PIC偏低。

    該研究的中國虎紋蛙的平均觀測雜合度(HO)和期望雜合度(HE)分別為0.460和0.582,該兩項指標與Chen等[11]對該物種的研究結(jié)果相似(HO=0.502,HE=0.609),但均低于Shao等[10]對該物種的研究結(jié)果(HO=0.796,HE=0.837),也低于其他無尾兩棲類(如沼水蛙[26]和海陸蛙[23]的HO和HE分別為0.488、0.635和0.789、0.669),進一步表明微衛(wèi)星位點的多態(tài)性與檢測樣品個體的親緣關(guān)系有關(guān)。另外,該研究發(fā)現(xiàn),有2個微衛(wèi)星位點顯著偏離哈迪溫伯格平衡,相似的研究結(jié)果也在海陸蛙[23]和飾紋姬蛙[25]等上發(fā)現(xiàn),這可能是由于2個位點之間的等位基因大小相似引起的[14,27]。

    4 結(jié)論

    該研究篩選的76個多堿基重復(fù)微衛(wèi)星位點豐富了中國虎紋蛙微衛(wèi)星文庫資源,研究結(jié)果將有助于開展中國虎紋蛙的遺傳多樣性、遺傳圖譜、種群遺傳結(jié)構(gòu)和人工育種研究。

    參考文獻

    [1]? 費梁,葉昌媛,江建平.中國兩棲動物及其分布彩色圖鑒[M].成都:四川科學(xué)技術(shù)出版社,2012.

    [2] WEI L,SHAO W W,DING G H,et al.Density but not kinship regulates the growth and developmental traits of Chinese tiger frog(Hoplobatrachus chinensis)Tadpoles[J].Asian herpetological research,2014,5(2):113-118.

    [3] TANG Y,CHEN Z Q,LIN Y F,et al.The combined effects of temperature and aromatase inhibitor on metamorphosis,growth,locomotion,and sex ratio of tiger frog(Hoplobatrachus rugulosus)tadpoles[J].PeerJ,2020,8:1-17.

    [4] WEI L,LIN Z H,MA X M,et al.Acoustic characteristics of the tiger frog,Hoplobatrachus rugulosus,during the breeding season[J].Zoological research,2011,32(4):456-460.

    [5] 李艷艷.親緣關(guān)系對虎紋蛙繁殖及其蝌蚪同類相食的影響[D].杭州:杭州師范大學(xué),2016.

    [6] 馬小浩.虎紋蛙的繁殖特征及其蝌蚪攻擊行為[D].杭州:杭州師范大學(xué),2012.

    [7] 符史杭,汪繼超,傅麗蓉,等.海南省人工養(yǎng)殖虎紋蛙的現(xiàn)狀及其對策[J].湖北農(nóng)業(yè)科學(xué),2013,52(5):1118-1120.

    [8] YU D N,ZHANG J Y,LI P,et al.Do cryptic species exist in Hoplobatrachus rugulosus? An examination using four nuclear genes,the Cyt b gene and the complete MT genome[J].PLoS One,2015,10(4):1-22.

    [9] MUNDKUR R,MUNIRAJU E.Molecular markerassisted selection breeding in silkworm,Bombyx mori[M]//KUMAR D,GONG C L.Trends in insect molecular biology and biotechnology.Cham:Springer International Publishing,2018.

    [10]? SHAO C,WANG Y,QIAO N N.Isolation and characterization of microsatellite loci in tiger frog(Hoplobatrachus rugulosus)[J].Conservation genetics,2009,10(5):1601-1603.

    [11] CHEN Z Q,DING G H,WANG Y M,et al.Development and characterization of 33 microsatellite loci for the tiger frog Hoplobatrachus rugulosus(Wiegmann 1834)through transcriptome sequencing[J].Journal of genetics,2018,97(5):e147-e151.

    [12] DIWAN N,CREGAN P B.Automated sizing of fluorescentlabeled simple sequence repeat(SSR)markers to assay genetic variation in soybean[J].Theoretical and applied genetics,1997,95:723-733.

    [13] 李強,萬建民.SSRHunter:一個本地化的SSR位點搜索軟件的開發(fā)[J].遺傳,2005,27(5):808-810.

    [14] ROUSSET F.Genepop007:A complete reimplementation of the genepop software for Windows and Linux[J].Molecular ecology resources,2008,8(1):103-106.

    [15] GOUDET J.FSTAT(version 1.2):A computer program to calculate Fstatistics[J].Journal of heredity,1995,86(6):485-486.

    [16] KALINOWSKI S T,TAPER M L,MARSHALL T C.Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment[J].Molecular ecology,2007,16(5):1099-1106.

    [17] SREENU V B,ALEVOOR V,NAGARAJU J,et al.MICdb:Database of prokaryotic microsatellites[J].Nucleic acids research,2003,31(1):106-108.

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

    [19] SHARMA P C,GROVER A,KAHL G.Mining microsatellites in eukaryotic genomes[J].Trends in biotechnology,2007,25(11):490-498.

    [20] ZANE L,BARGELLONI L,PATARNELLO T.Strategies for microsatellite isolation:A review[J].Molecular ecology,2002,11(1):1-16.

    [21] BEIER S,THIEL T,MNCH T,et al.MISAweb:A web server for microsatellite prediction[J].Bioinformatics,2017,33(16):2583-2585.

    [22] BOTSTEIN D,WHITE R L,SKOLNICK M,et al.Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J].American journal of human genetics,1980,32(3):314-331.

    [23] ZHENG Y T,WANG L J,HONG M L,et al.Isolation and characterization of 20 microsatellite markers for the crabeating frog Fejervarya cancrivora[J].Conservation genetics resources,2016,8(3):251-253.

    [24] PETERS M B,BEARD K H,HAGEN C,et al.Isolation of microsatellite loci from the coqui frog,Eleutherodactylus coqui[J].Molecular ecology resources,2008,8(1):139-141.

    [25] ZHANG J,WANG J,JIANG J P.Development and characterization of 15 polymorphic microsatellite markers for a highly adaptive and widerange frog(Microhyla fissipes)[J].Asian herpetological research,2013,4(3):202-206.

    [26] SHAO W W,LIN Z H,F(xiàn)AN X L,et al.Characterization of 25 microsatellite loci for the Guenthers frog Hylarana guentheri[J].Conservation genetics resources,2018,10(1):31-33.

    [27] WATTIER R,ENGEL C R,SAUMITOULAPRADE P,et al.Short allele dominance as a source of heterozygote deficiency at microsatellite loci:Experimental evidence at the dinucleotide locus Gv1CT in Gracilaria gracilis(Rhodophyta)[J].Molecular ecology,1998,7(11):1569-1573.

    基金項目?? 浙江省自然科學(xué)基金項目(LY19C040001);麗水市重點研究項目(2021ZDYF09,2020ZDYF07);百山祖國家公園科學(xué)研究項目(2021KFLY01)。

    作者簡介?? 喬芬(1983—),女,河北保定人,講師,博士,從事生態(tài)學(xué)研究。

    通信作者,教授,博士,從事動物學(xué)研究。

    鳴? 謝?? 該研究工作得到楊寧同學(xué)的幫助,在此表示感謝。

    猜你喜歡
    分子標記遺傳多樣性多態(tài)性
    單核苷酸多態(tài)性與中醫(yī)證候相關(guān)性研究進展
    軟棗獼猴桃性別相關(guān)的SRAP分子標記
    大白菜種質(zhì)資源抗根腫病基因CRa和CRb的分子標記鑒定與分析
    茄子種質(zhì)資源農(nóng)藝性狀遺傳多樣性分析
    玉米大斑病的研究進展
    淺析田間水稻紋枯病抗性鑒定體系的確立與完善
    西藏野核桃的表型特征及其保育措施
    水稻紋枯病抗性鑒定體系的確立與遺傳多樣性研究
    馬鈴薯cpDNA/mtDNA多態(tài)性的多重PCR檢測
    GlobalFiler~? PCR擴增試劑盒驗證及其STR遺傳多態(tài)性
    色婷婷久久久亚洲欧美| 色吧在线观看| 亚洲三区欧美一区| 精品国产国语对白av| 欧美久久黑人一区二区| 免费不卡黄色视频| 高清av免费在线| 国产男人的电影天堂91| 青青草视频在线视频观看| 我要看黄色一级片免费的| 久久久精品区二区三区| 伊人久久大香线蕉亚洲五| 国产亚洲精品第一综合不卡| 欧美日韩亚洲高清精品| 美女国产高潮福利片在线看| 国产精品嫩草影院av在线观看| 777久久人妻少妇嫩草av网站| av不卡在线播放| 国产日韩欧美视频二区| 涩涩av久久男人的天堂| 欧美黄色片欧美黄色片| 丝袜美腿诱惑在线| 一级,二级,三级黄色视频| 精品国产一区二区久久| 狂野欧美激情性xxxx| 欧美中文综合在线视频| 黄片播放在线免费| www.熟女人妻精品国产| 男女高潮啪啪啪动态图| 中国国产av一级| 久久久久久久精品精品| 久久青草综合色| 色播在线永久视频| 夜夜骑夜夜射夜夜干| 久久精品国产综合久久久| 咕卡用的链子| 韩国av在线不卡| 亚洲三区欧美一区| 亚洲国产成人一精品久久久| 国产野战对白在线观看| 婷婷色av中文字幕| 男女高潮啪啪啪动态图| 亚洲欧美成人精品一区二区| 啦啦啦 在线观看视频| 国产女主播在线喷水免费视频网站| 久久精品亚洲av国产电影网| 日韩大码丰满熟妇| 久久久久视频综合| 亚洲,一卡二卡三卡| 欧美国产精品va在线观看不卡| 亚洲欧美中文字幕日韩二区| 国产精品一区二区精品视频观看| 午夜精品国产一区二区电影| 国产一卡二卡三卡精品 | 午夜影院在线不卡| 国语对白做爰xxxⅹ性视频网站| 看十八女毛片水多多多| 久久精品人人爽人人爽视色| 亚洲激情五月婷婷啪啪| 久久久久精品久久久久真实原创| 亚洲美女视频黄频| 国产精品欧美亚洲77777| 久久青草综合色| 日日摸夜夜添夜夜爱| 大片免费播放器 马上看| 热99国产精品久久久久久7| 悠悠久久av| 巨乳人妻的诱惑在线观看| 精品一区二区三卡| 黑人欧美特级aaaaaa片| 精品视频人人做人人爽| 亚洲国产毛片av蜜桃av| 人妻一区二区av| 大香蕉久久网| 久久精品国产亚洲av高清一级| 欧美日韩一级在线毛片| 久久久久精品国产欧美久久久 | 热99国产精品久久久久久7| 亚洲美女搞黄在线观看| 成年动漫av网址| 五月开心婷婷网| 久久亚洲国产成人精品v| 极品少妇高潮喷水抽搐| 久久久久久久国产电影| 一级a爱视频在线免费观看| 男人爽女人下面视频在线观看| 韩国高清视频一区二区三区| 精品亚洲乱码少妇综合久久| 亚洲,欧美精品.| 狠狠婷婷综合久久久久久88av| 国产精品无大码| 亚洲国产中文字幕在线视频| 丁香六月欧美| 少妇猛男粗大的猛烈进出视频| 精品酒店卫生间| 精品一品国产午夜福利视频| 久热爱精品视频在线9| 青春草国产在线视频| 国产亚洲av高清不卡| 丁香六月天网| av在线老鸭窝| 国产精品国产三级国产专区5o| 中文乱码字字幕精品一区二区三区| 亚洲国产欧美日韩在线播放| 日韩视频在线欧美| 亚洲成av片中文字幕在线观看| 两性夫妻黄色片| 国产亚洲欧美精品永久| 丰满迷人的少妇在线观看| 亚洲精品国产av成人精品| 久久精品国产a三级三级三级| 国产精品久久久人人做人人爽| 在线观看免费视频网站a站| 国产免费又黄又爽又色| 宅男免费午夜| 午夜福利网站1000一区二区三区| 黄片无遮挡物在线观看| 中国三级夫妇交换| 一区二区三区四区激情视频| 狂野欧美激情性bbbbbb| 97在线人人人人妻| 一区福利在线观看| 老汉色av国产亚洲站长工具| 波野结衣二区三区在线| 一本大道久久a久久精品| 如日韩欧美国产精品一区二区三区| 18在线观看网站| 中文欧美无线码| 免费不卡黄色视频| 一级毛片 在线播放| 国产激情久久老熟女| 国产男人的电影天堂91| 久久久国产一区二区| 国产av国产精品国产| av电影中文网址| 中文字幕最新亚洲高清| 日韩,欧美,国产一区二区三区| 免费高清在线观看日韩| 大话2 男鬼变身卡| 久久天躁狠狠躁夜夜2o2o | 国产淫语在线视频| 亚洲av成人不卡在线观看播放网 | 天天躁夜夜躁狠狠久久av| 国产精品 国内视频| 中国三级夫妇交换| 久久精品国产综合久久久| 国产亚洲av高清不卡| av在线app专区| 不卡视频在线观看欧美| 亚洲av欧美aⅴ国产| 成人国产av品久久久| 天美传媒精品一区二区| 99热全是精品| 午夜福利视频在线观看免费| 日韩熟女老妇一区二区性免费视频| 久久99精品国语久久久| 国产激情久久老熟女| 色婷婷久久久亚洲欧美| 亚洲成人免费av在线播放| 日韩熟女老妇一区二区性免费视频| 亚洲熟女精品中文字幕| 无限看片的www在线观看| 久久婷婷青草| 日韩精品有码人妻一区| 最近2019中文字幕mv第一页| 国产探花极品一区二区| 在线观看三级黄色| 久久久久久人妻| 亚洲国产精品一区二区三区在线| 日本av手机在线免费观看| 国产乱人偷精品视频| 交换朋友夫妻互换小说| 精品卡一卡二卡四卡免费| 精品一区二区三卡| 精品国产一区二区久久| 亚洲欧美日韩另类电影网站| 色网站视频免费| 激情五月婷婷亚洲| 天天躁夜夜躁狠狠躁躁| 亚洲精品自拍成人| av卡一久久| 精品一品国产午夜福利视频| 性少妇av在线| 亚洲欧美一区二区三区国产| 久久国产精品大桥未久av| 欧美另类一区| 亚洲精品一区蜜桃| 国产一区二区三区av在线| 夜夜骑夜夜射夜夜干| 亚洲国产精品999| av片东京热男人的天堂| 另类精品久久| 男女免费视频国产| 国产1区2区3区精品| 国产精品久久久av美女十八| 国产一区二区三区av在线| 又黄又粗又硬又大视频| 捣出白浆h1v1| 午夜久久久在线观看| 女的被弄到高潮叫床怎么办| 日韩人妻精品一区2区三区| 免费观看性生交大片5| 国产xxxxx性猛交| 成人黄色视频免费在线看| 亚洲欧美清纯卡通| 性高湖久久久久久久久免费观看| www.熟女人妻精品国产| 日本欧美视频一区| 日本wwww免费看| 美女国产高潮福利片在线看| 看十八女毛片水多多多| 91老司机精品| 久久精品aⅴ一区二区三区四区| 最黄视频免费看| a 毛片基地| 19禁男女啪啪无遮挡网站| 在线观看三级黄色| 女人高潮潮喷娇喘18禁视频| 99香蕉大伊视频| 一边摸一边抽搐一进一出视频| 大香蕉久久网| 两性夫妻黄色片| 国产亚洲欧美精品永久| 老司机深夜福利视频在线观看 | 在线免费观看不下载黄p国产| 人人妻,人人澡人人爽秒播 | 精品亚洲成国产av| 在线天堂中文资源库| 少妇被粗大猛烈的视频| 国精品久久久久久国模美| 丝袜脚勾引网站| 巨乳人妻的诱惑在线观看| 波野结衣二区三区在线| 亚洲一区二区三区欧美精品| 国产精品亚洲av一区麻豆 | 精品亚洲成a人片在线观看| 午夜免费男女啪啪视频观看| 汤姆久久久久久久影院中文字幕| 国产视频首页在线观看| 国产黄色视频一区二区在线观看| 日韩制服丝袜自拍偷拍| 校园人妻丝袜中文字幕| 97精品久久久久久久久久精品| 美女中出高潮动态图| 男女之事视频高清在线观看 | 亚洲人成电影观看| 黄色毛片三级朝国网站| tube8黄色片| 高清黄色对白视频在线免费看| 久久久精品国产亚洲av高清涩受| 在线观看国产h片| 97人妻天天添夜夜摸| 亚洲视频免费观看视频| 久久99精品国语久久久| 亚洲 欧美一区二区三区| 欧美成人精品欧美一级黄| 美女高潮到喷水免费观看| 亚洲av男天堂| 日韩 亚洲 欧美在线| 欧美老熟妇乱子伦牲交| 欧美亚洲 丝袜 人妻 在线| 亚洲欧洲日产国产| 精品国产超薄肉色丝袜足j| 男女床上黄色一级片免费看| 亚洲精品,欧美精品| 多毛熟女@视频| 亚洲av综合色区一区| 伊人亚洲综合成人网| av国产久精品久网站免费入址| 尾随美女入室| 欧美中文综合在线视频| 日本91视频免费播放| 少妇被粗大猛烈的视频| 老司机靠b影院| 一级片免费观看大全| 欧美日韩国产mv在线观看视频| 尾随美女入室| 久久午夜综合久久蜜桃| 国产成人午夜福利电影在线观看| 熟女少妇亚洲综合色aaa.| 嫩草影视91久久| 日本av免费视频播放| 亚洲欧美一区二区三区国产| 国产精品久久久久久人妻精品电影 | 在现免费观看毛片| 狠狠婷婷综合久久久久久88av| 国产成人免费无遮挡视频| 久久av网站| 欧美中文综合在线视频| 成年动漫av网址| 最新在线观看一区二区三区 | 国产色婷婷99| 日韩中文字幕视频在线看片| 国产一区有黄有色的免费视频| 亚洲av成人精品一二三区| 建设人人有责人人尽责人人享有的| 亚洲一码二码三码区别大吗| 无限看片的www在线观看| 性高湖久久久久久久久免费观看| 高清不卡的av网站| 不卡视频在线观看欧美| 精品少妇久久久久久888优播| 午夜福利一区二区在线看| 中文字幕色久视频| 免费黄频网站在线观看国产| 七月丁香在线播放| 男女午夜视频在线观看| 最近最新中文字幕免费大全7| 欧美日韩精品网址| 国产不卡av网站在线观看| 国产黄色免费在线视频| 欧美久久黑人一区二区| 精品少妇一区二区三区视频日本电影 | 最近的中文字幕免费完整| 午夜免费观看性视频| 久久精品亚洲av国产电影网| 欧美乱码精品一区二区三区| 久久亚洲国产成人精品v| 一本—道久久a久久精品蜜桃钙片| a 毛片基地| 五月开心婷婷网| 777米奇影视久久| 亚洲欧美激情在线| 国产一区二区激情短视频 | 美女中出高潮动态图| 街头女战士在线观看网站| 大片免费播放器 马上看| h视频一区二区三区| 日本色播在线视频| 午夜福利影视在线免费观看| kizo精华| 日韩伦理黄色片| 啦啦啦视频在线资源免费观看| 啦啦啦在线观看免费高清www| 欧美久久黑人一区二区| 久久人人97超碰香蕉20202| 啦啦啦 在线观看视频| 亚洲美女黄色视频免费看| 精品少妇内射三级| 不卡视频在线观看欧美| 久久久国产一区二区| 精品人妻一区二区三区麻豆| 日韩大码丰满熟妇| 日韩欧美精品免费久久| 中文字幕人妻丝袜一区二区 | 99热全是精品| 免费在线观看完整版高清| 街头女战士在线观看网站| 美国免费a级毛片| 久久久精品免费免费高清| 一区二区av电影网| 一区二区三区激情视频| 叶爱在线成人免费视频播放| 欧美另类一区| 乱人伦中国视频| 亚洲成人国产一区在线观看 | 国产熟女欧美一区二区| 国产一区二区三区av在线| 国产男女内射视频| 日韩一区二区三区影片| 国产精品一区二区在线观看99| 十分钟在线观看高清视频www| 色视频在线一区二区三区| 亚洲第一av免费看| 亚洲av福利一区| 亚洲成人一二三区av| 国产福利在线免费观看视频| 九色亚洲精品在线播放| 日韩人妻精品一区2区三区| 国产精品久久久久久精品电影小说| 男女下面插进去视频免费观看| 久久久久网色| 视频在线观看一区二区三区| 亚洲精品在线美女| 美女主播在线视频| 亚洲伊人色综图| 国产成人免费观看mmmm| 美女午夜性视频免费| 高清av免费在线| 欧美黑人精品巨大| 精品一区在线观看国产| 中文字幕人妻丝袜一区二区 | 久久人妻熟女aⅴ| 日本欧美国产在线视频| 电影成人av| 亚洲五月色婷婷综合| 超色免费av| 国产精品国产三级国产专区5o| 久久久久久久久久久久大奶| 国产精品99久久99久久久不卡 | 国产欧美日韩一区二区三区在线| 亚洲国产av新网站| 天天影视国产精品| 欧美少妇被猛烈插入视频| 国产一卡二卡三卡精品 | 日韩欧美一区视频在线观看| 伦理电影免费视频| 在线精品无人区一区二区三| 日日爽夜夜爽网站| 国产精品av久久久久免费| av电影中文网址| 制服丝袜香蕉在线| 水蜜桃什么品种好| 只有这里有精品99| 18禁观看日本| 桃花免费在线播放| 亚洲精品久久成人aⅴ小说| av一本久久久久| 热99久久久久精品小说推荐| 国产精品女同一区二区软件| 九草在线视频观看| 欧美日韩国产mv在线观看视频| 国产熟女欧美一区二区| 欧美日韩精品网址| 香蕉国产在线看| 女人被躁到高潮嗷嗷叫费观| 精品一区二区三区四区五区乱码 | 国产伦理片在线播放av一区| 国产日韩一区二区三区精品不卡| 国产成人精品福利久久| 又黄又粗又硬又大视频| 一级,二级,三级黄色视频| 韩国高清视频一区二区三区| 视频区图区小说| 丁香六月天网| 国产精品欧美亚洲77777| 亚洲第一青青草原| 美女福利国产在线| 男女之事视频高清在线观看 | 蜜桃国产av成人99| 少妇人妻久久综合中文| 欧美在线黄色| 久久午夜综合久久蜜桃| 女性生殖器流出的白浆| 日本欧美国产在线视频| 18在线观看网站| 老司机靠b影院| 国产一区二区在线观看av| 搡老岳熟女国产| 欧美另类一区| 国产欧美亚洲国产| 一区二区三区乱码不卡18| 午夜福利网站1000一区二区三区| 狂野欧美激情性bbbbbb| 麻豆av在线久日| 99热全是精品| 久久99精品国语久久久| 韩国高清视频一区二区三区| 91国产中文字幕| 免费黄色在线免费观看| 国产黄频视频在线观看| 亚洲熟女毛片儿| 国产成人精品无人区| 国产精品熟女久久久久浪| 黄频高清免费视频| 欧美日韩精品网址| 91精品伊人久久大香线蕉| 亚洲欧美成人精品一区二区| 国产日韩一区二区三区精品不卡| 尾随美女入室| 午夜免费鲁丝| 国产免费现黄频在线看| 精品久久久精品久久久| 丝袜美腿诱惑在线| 亚洲国产精品成人久久小说| 母亲3免费完整高清在线观看| 下体分泌物呈黄色| 亚洲av国产av综合av卡| 青春草国产在线视频| 久久99精品国语久久久| 女的被弄到高潮叫床怎么办| 久久精品久久久久久噜噜老黄| 国产在线一区二区三区精| 国产精品亚洲av一区麻豆 | 可以免费在线观看a视频的电影网站 | 一本一本久久a久久精品综合妖精| 国产成人av激情在线播放| 黄色一级大片看看| 一区二区三区乱码不卡18| 美女脱内裤让男人舔精品视频| 女人爽到高潮嗷嗷叫在线视频| 51午夜福利影视在线观看| 在线免费观看不下载黄p国产| 最黄视频免费看| 一级爰片在线观看| 亚洲国产精品一区三区| 亚洲成人av在线免费| 日韩 亚洲 欧美在线| 婷婷色综合大香蕉| 精品一区二区三卡| 久久久久久久精品精品| 天堂中文最新版在线下载| 国产精品国产三级国产专区5o| 少妇人妻 视频| 日韩精品免费视频一区二区三区| 久久精品久久精品一区二区三区| 亚洲欧美一区二区三区黑人| 一区二区三区激情视频| e午夜精品久久久久久久| 另类亚洲欧美激情| 亚洲婷婷狠狠爱综合网| 久久久久久久久久久久大奶| 久久久久久久国产电影| 午夜久久久在线观看| 啦啦啦 在线观看视频| 99久国产av精品国产电影| 国产精品国产三级专区第一集| 亚洲图色成人| √禁漫天堂资源中文www| 久久鲁丝午夜福利片| 免费日韩欧美在线观看| 国产精品无大码| 无遮挡黄片免费观看| 欧美少妇被猛烈插入视频| 性高湖久久久久久久久免费观看| 午夜激情av网站| 国产极品粉嫩免费观看在线| 秋霞伦理黄片| 中文字幕av电影在线播放| 亚洲精品aⅴ在线观看| 婷婷色av中文字幕| xxx大片免费视频| 国产伦理片在线播放av一区| 无限看片的www在线观看| 国产精品成人在线| 欧美av亚洲av综合av国产av | 叶爱在线成人免费视频播放| 下体分泌物呈黄色| 成人黄色视频免费在线看| 搡老岳熟女国产| 久久久久久久精品精品| 国产福利在线免费观看视频| 亚洲国产成人一精品久久久| 亚洲国产中文字幕在线视频| videosex国产| 亚洲精品成人av观看孕妇| 国产成人91sexporn| 亚洲视频免费观看视频| 国语对白做爰xxxⅹ性视频网站| av在线观看视频网站免费| 日本91视频免费播放| 欧美久久黑人一区二区| 男女之事视频高清在线观看 | 熟妇人妻不卡中文字幕| 国产亚洲欧美精品永久| 伦理电影大哥的女人| 99久国产av精品国产电影| 亚洲天堂av无毛| 日韩电影二区| 日日爽夜夜爽网站| 色精品久久人妻99蜜桃| 高清在线视频一区二区三区| 国产高清国产精品国产三级| 香蕉丝袜av| 蜜桃国产av成人99| 熟女少妇亚洲综合色aaa.| 啦啦啦 在线观看视频| 97人妻天天添夜夜摸| 亚洲国产av影院在线观看| 免费看不卡的av| 亚洲精品国产av成人精品| 午夜av观看不卡| 天天影视国产精品| 国产成人免费无遮挡视频| 亚洲一区二区三区欧美精品| 在线观看人妻少妇| 91老司机精品| 欧美在线一区亚洲| av线在线观看网站| 丰满迷人的少妇在线观看| 亚洲欧美精品自产自拍| 欧美 日韩 精品 国产| 国产精品三级大全| 色网站视频免费| 国产无遮挡羞羞视频在线观看| 中文精品一卡2卡3卡4更新| 成人免费观看视频高清| 国产日韩欧美视频二区| 精品久久蜜臀av无| 伦理电影大哥的女人| 国产亚洲av高清不卡| 成年女人毛片免费观看观看9 | 亚洲精品aⅴ在线观看| 伊人久久大香线蕉亚洲五| 国产精品一区二区在线不卡| 成年美女黄网站色视频大全免费| 亚洲国产欧美日韩在线播放| 日韩av免费高清视频| 黄色一级大片看看| 韩国av在线不卡| 国产精品人妻久久久影院| 伦理电影大哥的女人| 国产乱来视频区| 国产成人欧美在线观看 | 午夜福利网站1000一区二区三区| 在线 av 中文字幕| 男女下面插进去视频免费观看| 亚洲五月色婷婷综合| 国产熟女午夜一区二区三区| 免费黄网站久久成人精品| 国产成人免费观看mmmm| 午夜免费观看性视频| 9色porny在线观看| videosex国产| 无限看片的www在线观看| 日韩av在线免费看完整版不卡| 日韩中文字幕欧美一区二区 | 精品一区二区免费观看| 一区福利在线观看| 精品一区二区三区av网在线观看 | 欧美在线一区亚洲| a级毛片黄视频| av在线观看视频网站免费| 三上悠亚av全集在线观看| 各种免费的搞黄视频|