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

    ?基于轉(zhuǎn)錄組測序的紫薇SSR特征分析?

    2023-04-29 00:00:00秦波韋廣綏覃杰孫開道秦賜梅劉昊黃欣
    湖南農(nóng)業(yè)科學(xué) 2023年4期

    摘 要:以紫薇葉片為材料,通過Illumina Hiseq X Ten測序平臺進行高通量測序,利用MISA軟件分析轉(zhuǎn)錄組的SSR相關(guān)信息。結(jié)果表明:測序共獲得31 461條Unigene序列,總長度34 408 463 bp,其中有5 617條Unigene序列含有SSR位點,共發(fā)現(xiàn)SSR位點6 761個,出現(xiàn)頻率為21.49%,平均分布距離為5.09 kb;主要的優(yōu)勢重復(fù)類型為二核苷酸重復(fù)和三核苷酸重復(fù),分別占總SSR的44.95%和36.41%;AG/CT和A/T為SSR優(yōu)勢重復(fù)基元,出現(xiàn)頻率分別為40.23%(2 720個)和16.57%(1 120個)。分析表明,紫薇轉(zhuǎn)錄組中SSR類型豐富、位點出現(xiàn)頻率高、多態(tài)性潛力高,在紫薇種質(zhì)遺傳多樣性研究、品種鑒定、雜種鑒定等方面具有應(yīng)用價值。

    關(guān)鍵詞:紫薇;轉(zhuǎn)錄組;SSR

    中圖分類號:Q754 文獻標識碼:A 文章編號:1006-060X(2023)04-0001-04

    Abstract:Illumina Hiseq X Ten sequencing platform was used to conduct high-throughput sequencing on Lagerstroemia indica leaves, and SSR related information of transcriptome was analyzed by MISA software. The results showed that a total of 31 461 Unigene sequences with a total length of 34 408 463 bp were obtained, among which 5 617 Unigene sequences contained SSR loci. A total of 6 761 SSR loci were found, with a frequency of 21.49% and an average distribution distance of 5.09 kb. The dominant repeats were dinucleotide and trinucleotide, accounting for 44.95% and 36.41% of the total SSR repeats, respectively. AG/CT and A/T were SSR dominant repeat primitives, with occurrence frequencies of 40.23% (2 720) and 16.57% (1 120), respectively. The results showed that SSR types were abundant, loci appeared frequently, and polymorphism potential was high in the transcriptome of Lagerstroemia indica, which has application value in genetic diversity analysis, variety identification, hybrid identification and so on.

    Key words:Lagerstroemia indica; transcriptome; SSR

    紫薇(Lagerstroemia indica)屬千屈菜科紫薇屬落葉灌木或小喬木,具有花期長、花量大、抗逆性強等特點,是重要的夏季觀賞木本花卉,被廣泛應(yīng)用于世界各地的園林景觀中[1-2]。紫薇栽培歷史悠久,最早的栽培記錄可追溯至我國東晉王嘉所著的《拾遺記》。經(jīng)過長期選育,紫薇形成了豐富多樣的品種,目前全世界的紫薇品種已超過500個[3]。

    SSR(Simple Sequence Repeats)又稱為微衛(wèi)星DNA(Microsatellite DNA),是由1~6個核苷酸組成的短串聯(lián)重復(fù)序列,具有可變的重復(fù)次數(shù)和保守的側(cè)翼序列,因而可通過設(shè)計特異引物進行PCR以擴增串聯(lián)重復(fù)序列的多態(tài)性[4]。SSR具有共顯性遺傳、多態(tài)性豐富、重復(fù)性高等特點,被廣泛應(yīng)用于物種親緣關(guān)系分析、遺傳多樣性分析、品種鑒定、雜種鑒定等研究中。董俊美等[5]基于山藥轉(zhuǎn)錄組數(shù)據(jù)挖掘了大量的SSR位點,并通過開發(fā)分子標記對來自河北省、河南省、山東省及福建省等山藥主產(chǎn)區(qū)的43份山藥種質(zhì)資源進行了遺傳多樣性和親緣關(guān)系分析;胡光明等[6]基于獼猴桃全基因組數(shù)據(jù)篩選得到一批多態(tài)性高、通用性強的SSR引物,并對種質(zhì)資源進行了基因分型及親緣聚類分析;張成才等[7]通過篩選薄殼山核桃核心SSR引物,構(gòu)建了36個常見品種的特異分子指紋圖譜和分子身份證;谷艷鵬等[8]以花楸樹為母本、少葉花楸為父本進行種間雜交,利用篩選出的9對擴增結(jié)果穩(wěn)定且具多態(tài)性的EST-SSR引物,對雜交F1代群體進行雜種鑒定,并對親本和雜交F1代群體的基因分離規(guī)律以及遺傳多樣性和遺傳關(guān)系進行分析。

    隨著轉(zhuǎn)錄組測序技術(shù)的不斷進步,其成本越來越低,利用轉(zhuǎn)錄組開發(fā)SSR分子標記的技術(shù)也隨之得到普遍應(yīng)用。該研究基于紫薇轉(zhuǎn)錄組高通量測序結(jié)果,對分布于轉(zhuǎn)錄組上的SSR位點信息及其序列特征進行分析,以期為紫薇種質(zhì)資源親緣關(guān)系分析、遺傳多樣性分析、品種鑒定、雜種鑒定等提供基礎(chǔ)資料。

    1 材料與方法

    1.1 試驗材料

    所用材料為紫薇的新鮮幼嫩葉片,采集于廣西壯族自治區(qū)林業(yè)科學(xué)研究院紫薇種質(zhì)資源圃。采集時用干凈紙巾擦掉葉片表面的雜質(zhì),接著用離心管裝好迅速放入液氮中,帶回實驗室置于-80℃冰箱中保存。

    1.2 試驗方法

    1.2.1 RNA提取 使用Trizol法從紫薇樣品中提取總RNA,使用DNase I(TaKara)去除基因組DNA。文庫質(zhì)量檢測合格后進行測序。

    1.2.2 測序與拼接 通過Illumina Hiseq X Ten測序平臺進行高通量測序,獲得Raw reads后,去除基因列中包含接頭、含N和低質(zhì)量的參考基因合成序列,從而得到Clean reads。接著采用Trinity軟件進行拼接,將每條基因中獲得的最長轉(zhuǎn)錄本作為Unigene,用于后續(xù)分析。

    1.2.3 SSR位點搜索和數(shù)據(jù)分析 使用MISA軟件對獲得的Unigene進行SSR位點搜索,搜索標準為:含有1、2、3、4、5、6核苷酸重復(fù)的最小重復(fù)次數(shù)分別為10、6、5、5、5次。復(fù)合型核苷酸2個位點之間距離小于150 bp。將生成的數(shù)據(jù)用Microsoft Excel 2016軟件進行整理,對序列特征進行分析并繪制圖表。SSR的出現(xiàn)頻率為SSR位點總數(shù)與Unigene序列總數(shù)的比值;SSR位點的平均距離為序列總長度與SSR位點總數(shù)之比。

    2 結(jié)果與分析

    2.1 紫薇轉(zhuǎn)錄組SSR的分布特征

    通過高通量測序技術(shù)獲得紫薇葉片轉(zhuǎn)錄組數(shù)據(jù),接著利用MISA軟件進行SSR位點搜索,結(jié)果表明,在31 461條Unigenes序列中發(fā)現(xiàn)總計有5 617條序列含有SSR位點,其中含1個以上 SSR位點的序列有955條,搜索到的SSR位點總數(shù)為6 761個,SSR位點的出現(xiàn)頻率為21.49%。同時,搜索到包含復(fù)合型SSR位點的序列數(shù)為517個,占序列總數(shù)的9.2%。序列總長度為34 408 463 bp,SSR位點的分布距離為5.09 kb,即平均每5.09 kb出現(xiàn)一個SSR位點。

    2.2 紫薇轉(zhuǎn)錄組SSR重復(fù)類型分析

    紫薇轉(zhuǎn)錄組SSR重復(fù)類型情況如表1所示。各重復(fù)類型SSR按照豐度排序結(jié)果為二核苷酸重復(fù)>三核苷酸重復(fù)>單核苷酸重復(fù)>四核苷酸重復(fù)>五核苷酸重復(fù)>六核苷酸重復(fù)。紫薇轉(zhuǎn)錄組中的SSR主要重復(fù)類型為二核苷酸重復(fù),數(shù)量為3 039個,占比為44.95%,出現(xiàn)頻率為7.83%,平均距離13.98 kb;三核苷酸重復(fù)類型豐度次之,為2 462個,占比36.41%,出現(xiàn)頻率為9.66%,平均距離11.32 kb;第三是單核苷酸重復(fù),數(shù)量為1 132個,占比16.74%,出現(xiàn)頻率為3.60%,平均距離30.40 kb。其中,二核苷酸重復(fù)、三核苷酸重復(fù)和單核苷酸重復(fù)數(shù)量之和達到所有SSR數(shù)量的98.1%。

    2.3 紫薇轉(zhuǎn)錄組SSR重復(fù)基元分析

    由表2可知,由單核苷酸重復(fù)到六核苷酸重復(fù),重復(fù)基元的種類分別為2、4、10、19、18、16種,共69種。單核苷酸重復(fù)的優(yōu)勢重復(fù)單元為A/T,數(shù)量為1 120個,占單核苷酸重復(fù)的98.94%,出現(xiàn)頻率為16.57%;二核苷酸重復(fù)的優(yōu)勢重復(fù)單元為AG/CT,數(shù)量為2 720個,占二核苷酸重復(fù)的89.50%,出現(xiàn)頻率為40.23%;三核苷酸重復(fù)的優(yōu)勢重復(fù)單元為AGG/CCT,數(shù)量為751個,占三核苷酸重復(fù)的30.50%,出現(xiàn)頻率為11.11%。四、五、六核苷酸重復(fù)的優(yōu)勢重復(fù)單元分別為ACAT/ATGT、AGAGG/CCTCT、AAGGAG/CCTTCT,出現(xiàn)頻率均低于1%,分別為0.22%、0.09%和0.04%。轉(zhuǎn)錄組中主要SSR重復(fù)基元類型的和比例詳見圖1。

    2.4 紫薇轉(zhuǎn)錄組SSR重復(fù)次數(shù)分析

    SSR重復(fù)基元的重復(fù)次數(shù)決定了SSR序列長度的變化范圍,與其多態(tài)性密切相關(guān)。從表3可知,紫薇SSR位點重復(fù)基元的重復(fù)次數(shù)介于5~12次的比例達到90.56%,重復(fù)次數(shù)高于12次的比例為9.44%(638個)。頻率最高的是重復(fù)6次,共1 510個SSR位點,占比22.33%;其次是重復(fù)5次(1 288個,占比19.05%);按占比的高低將剩余的重復(fù)次數(shù)排列依次是重復(fù)7次(827個,占比12.23%)、重復(fù)10次(790個,占比11.68%)、重復(fù)8次(640個,占比9.47%)、重復(fù)11次(427個,占比6.32%)、重復(fù)9次(336個,占比4.97%)和重復(fù)12次(305個,占比4.51%)。

    3 討論與結(jié)論

    由于成本低、效率高等優(yōu)點,利用轉(zhuǎn)錄組開發(fā)SSR分子標記已成為常規(guī)手段。該研究基于紫薇轉(zhuǎn)錄組數(shù)據(jù),共獲得Unigene序列31 461條,發(fā)現(xiàn)有5 617條Unigene序列含有SSR位點,SSR位點數(shù)量總計6 761個,出現(xiàn)頻率為21.49%。該研究中紫薇SSR的出現(xiàn)頻率高于黃稈烏哺雞竹的14.07%[9]、八角的14.24%[10]、柴胡的12.24%[11],低于小果甜柿的31.52%[12]、睡蓮的33.98%[13],與多頭菊的19.77%[14]、藍莓的20.81%[15]接近。在已有的利用轉(zhuǎn)錄組開發(fā)紫薇SSR的研究中得到的發(fā)生頻率為21.53%和20.03%[16-17],與該研究的21.49%接近。通過搜索獲得序列的總長度為344 08 463 bp,SSR的分布密度為1/5.09 kb,與已有紫薇研究中的1/5.60、1/5.46 kb接近[16-17],遠大于辣椒的1/56 kb[18],低于野柿的1/2.94 kb[19]。

    SSR的出現(xiàn)頻率和分布密度在不同物種間存在差異,而紫薇中SSR分布較緊密,數(shù)量較豐富。紫薇SSR主要集中在二核苷酸重復(fù)類型(3 039個,占比44.95%)、三核苷酸重復(fù)類型(2 462個,占比36.41%)和單核苷酸重復(fù)(1 132個,占比16.74%),三者總和為98.1%。已有研究表明,植物轉(zhuǎn)錄組測序獲得的SSR以單、二、三核苷酸重復(fù)基因為主,該研究結(jié)果與前人的研究結(jié)果一致。同時前人研究發(fā)現(xiàn),SSR主導(dǎo)重復(fù)基元無明顯規(guī)律,分析認為這與物種差異、搜索參數(shù)設(shè)置差異等相關(guān)[20-21]。進一步分析發(fā)現(xiàn),AG/CT的重復(fù)頻率最高(40.23%),隨后依次是A/T(16.57%)、AGG/CCT(11.11%)、AAG/CTT(10.81%)、AGC/CTG(5.72%),而其他重復(fù)基元占比均小于2.5%。雙子葉植物常見的重復(fù)基元有AG/ CT、AT/AT、AAG/ CTT、AG/ CT、AGC/ TCG等,分析認為這可能與密碼子的偏好性有關(guān)[22]。該研究對紫薇轉(zhuǎn)錄組SSR不同重復(fù)類型在不同重復(fù)下的數(shù)量進行統(tǒng)計,結(jié)果表明,在各重復(fù)類型中,重復(fù)次數(shù)最高的為6次(1 510個,占比22.33%),隨后依次是5次(1 288個,占比19.05%)、7次(827個,占比12.23%)、10次(790個,占比11.68%),其余重復(fù)次數(shù)(不含12次以上)的占比介于4.51%~9.47%之間。SSR的多態(tài)性與重復(fù)基元類型、頻率、重復(fù)次數(shù)等密切相關(guān),該研究表明利用轉(zhuǎn)錄組開發(fā)具有多態(tài)性的紫薇SSR在理論上是可行的。

    參考文獻:

    [1] 張啟翔. 紫薇品種分類及其在園林中的應(yīng)用(英文)[J]. 北京林業(yè)大學(xué)學(xué)報,1991,13(4):57-66.

    [2] QIN B,SUN K D,HUANG X. The complete chloroplast genome of Lagerstroemia balansae,an endangered species of genus Lagerstroemia native to China[J]. Mitochondrial DNA Part B,2021,6(2):684-685.

    [3] CAI M,PAN H T,WANG X F,et al. Development of novel microsatellites in Lagerstroemia indica and DNA fingerprinting in Chinese Lagerstroemia cultivars[J]. Scientia Horticulturae,2011,131:88-94.

    [4] DAVIERWALA A P,RAMAKRISHNA W,RANJEKAR P K,et al. Sequence variations at a complex microsatellite locus in rice and its conservation in cereals[J]. Theoretical and Applied Genetics,2000,101(8):1291-1298.

    [5] 董俊美,李錦超,賈凱旋,等. 基于農(nóng)藝性狀的山藥種質(zhì)資源遺傳多樣性分析[J]. 江蘇農(nóng)業(yè)科學(xué),2022,50(18):136-143.

    [6] 胡光明,張 瓊,韓 飛,等. 獼猴桃屬植物通用型SSR分子標記引物的篩選及應(yīng)用[J]. 中國農(nóng)業(yè)科學(xué),2022,55(17):3411-3425.

    [7] 張成才,王亞萍,王開良,等. 利用SSR構(gòu)建薄殼山核桃主要品種的分子身份證[J]. 植物遺傳資源學(xué)報,2022,23(6):1835-1846.

    [8] 谷艷鵬,張澤人,孫 濤,等. 基于EST-SSR標記的花楸樹和少葉花楸雜交F1代群體的雜種鑒定及遺傳關(guān)系分析[J]. 植物資源與環(huán)境學(xué)報,2022,31(4):65-73.

    [9] 婁永峰,甘 然,朱成磊,等. 黃稈烏哺雞竹轉(zhuǎn)錄組EST-SSR分子標記開發(fā)與應(yīng)用[J]. 植物科學(xué)學(xué)報,2022,40(3):355-364.

    [10] 余興華,楊艷梅,李玉祥,等. 基于八角轉(zhuǎn)錄組序列的SSR分子標記引物開發(fā)[J]. 西南農(nóng)業(yè)學(xué)報,2022,35(6):1261-1268.

    [11] 劉 迪,歐陽艷飛,徐" 鵬,等. 柴胡轉(zhuǎn)錄組SSR的分布及序列特征分析[J]. 南方農(nóng)業(yè)學(xué)報,2022,53(3):676-683.

    [12] 王藝儒,索玉靜,傅建敏. 小果甜柿果實轉(zhuǎn)錄組的SSR、SNP和InDel特征分析[J]. 西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版),2022,50(7):147-154.

    [13] 蘇 群,田 敏,劉 俊,等. 基于生物信息學(xué)的睡蓮SSR位點特征分析[J]. 西南農(nóng)業(yè)學(xué)報,2021,34(10):2076-2083.

    [14] 吳 超,郭方其,徐丹彬,等. 多頭菊頭狀花序轉(zhuǎn)錄組SSR特點分析[J]. 分子植物育種,2022,20(18):6067-6075.

    [15] 王小敏,張春紅,胡麗超,等. 基于高通量測序的藍莓果實轉(zhuǎn)錄組SSR信息分析[J]. 北方園藝,2022(13):49-54.

    [16] 李素珍,胡 玲,蔡 明,等. 基于轉(zhuǎn)錄組的紫薇SSR分子標記開發(fā)[A]. 中國觀賞園藝研究進展2018[C]. 北京:中國林業(yè)出版社,2018. 526-532.

    [17] 王 鵬,李 亞,王淑安,等. 紫薇EST-SSR標記的開發(fā)和利用[A]. 中國觀賞園藝研究進展2016[C]. 北京:中國林業(yè)出版社,2016. 220-225.

    [18] 劉 峰,王運生,田雪亮,等. 辣椒轉(zhuǎn)錄組SSR挖掘及其多態(tài)性分析[J]. 園藝學(xué)報,2012,39(1):168-174.

    [19] 李樹戰(zhàn),王藝儒,孫 鵬,等. ‘龍巖野柿1號’雄花和兩性花花芽的轉(zhuǎn)錄組SSR特征分析[J]. 中南林業(yè)科技大學(xué)學(xué)報,2022,42(9):171-177,186.

    [20] 烏日娜,徐 舶,石鳳翎. 基于轉(zhuǎn)錄組數(shù)據(jù)的直立型扁蓿豆SSR序列特征分析[J]. 西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版),2022,50(5):1-8.

    [21] 公培民,寧書菊,葉 齊,等. 馬藍轉(zhuǎn)錄組中SSR位點信息特征分析[J]. 分子植物育種,2020,18(12):4000-4008.

    [22] CHAGNé D,CHAUMEIL P,RAMBOER A,et al. Cross-species transferability and mapping of genomic and cDNA SSRs in pines[J]. Theoretical and Applied Genetics,2004,109(6):1204-1214.

    (責(zé)任編輯:成 平)

    久久国产精品影院| 观看免费一级毛片| 国产黄a三级三级三级人| 精品国产亚洲在线| 日本撒尿小便嘘嘘汇集6| 97超视频在线观看视频| 一级毛片女人18水好多| 免费大片18禁| 少妇的逼水好多| 床上黄色一级片| 午夜日韩欧美国产| 午夜福利视频1000在线观看| 精品久久久久久,| 国产精品99久久久久久久久| 久久久国产成人精品二区| 欧美黑人巨大hd| 九九久久精品国产亚洲av麻豆 | 看片在线看免费视频| 性色avwww在线观看| 国产高清视频在线观看网站| 亚洲熟女毛片儿| 免费一级毛片在线播放高清视频| 美女 人体艺术 gogo| 琪琪午夜伦伦电影理论片6080| 午夜免费成人在线视频| 美女被艹到高潮喷水动态| 欧美色视频一区免费| 久久久久久国产a免费观看| 久久国产乱子伦精品免费另类| 日本黄色视频三级网站网址| 成人永久免费在线观看视频| 桃红色精品国产亚洲av| 国产蜜桃级精品一区二区三区| 一二三四在线观看免费中文在| 精品久久久久久久久久久久久| 免费观看的影片在线观看| 亚洲性夜色夜夜综合| 国产男靠女视频免费网站| 中文资源天堂在线| 成人18禁在线播放| 午夜精品久久久久久毛片777| a级毛片a级免费在线| 色在线成人网| 国产黄片美女视频| 9191精品国产免费久久| 国产综合懂色| 两性夫妻黄色片| 岛国在线免费视频观看| 香蕉av资源在线| 亚洲成人久久性| 欧美高清成人免费视频www| 亚洲精品456在线播放app | 久久精品影院6| 村上凉子中文字幕在线| 小说图片视频综合网站| 国产欧美日韩精品一区二区| 在线播放国产精品三级| 久久精品91无色码中文字幕| 美女高潮喷水抽搐中文字幕| ponron亚洲| 亚洲av五月六月丁香网| 日韩欧美国产在线观看| 在线观看免费午夜福利视频| 亚洲av电影在线进入| 久久99热这里只有精品18| 亚洲中文字幕一区二区三区有码在线看 | 国产亚洲欧美98| 一本久久中文字幕| or卡值多少钱| 久久婷婷人人爽人人干人人爱| 两个人的视频大全免费| 欧美日本视频| 色av中文字幕| 亚洲色图 男人天堂 中文字幕| 一进一出好大好爽视频| 久久久久久人人人人人| 美女被艹到高潮喷水动态| 欧美国产日韩亚洲一区| 成人鲁丝片一二三区免费| 国产黄a三级三级三级人| 日韩免费av在线播放| 岛国在线观看网站| 看免费av毛片| 欧美黄色淫秽网站| 久久午夜亚洲精品久久| 精品国产乱码久久久久久男人| 色噜噜av男人的天堂激情| 亚洲自偷自拍图片 自拍| 亚洲av成人一区二区三| 亚洲欧美日韩卡通动漫| 成年人黄色毛片网站| 99热6这里只有精品| 中文字幕人成人乱码亚洲影| 亚洲电影在线观看av| 夜夜夜夜夜久久久久| 日韩欧美一区二区三区在线观看| 国产精品久久久久久久电影 | 亚洲欧美日韩高清专用| 成人特级av手机在线观看| 国产一区二区三区视频了| 午夜a级毛片| 两个人的视频大全免费| 日韩三级视频一区二区三区| 久久精品国产99精品国产亚洲性色| 国产野战对白在线观看| 成人欧美大片| 亚洲欧美精品综合久久99| 欧美色视频一区免费| 亚洲国产欧美一区二区综合| 国产一区二区在线av高清观看| 2021天堂中文幕一二区在线观| 999久久久精品免费观看国产| 又黄又粗又硬又大视频| 精品一区二区三区av网在线观看| 欧美成人免费av一区二区三区| 别揉我奶头~嗯~啊~动态视频| 在线观看午夜福利视频| 搡老妇女老女人老熟妇| 久久久水蜜桃国产精品网| 在线永久观看黄色视频| 国产精品久久久av美女十八| 99久久久亚洲精品蜜臀av| 成人高潮视频无遮挡免费网站| 亚洲国产日韩欧美精品在线观看 | 又爽又黄无遮挡网站| 午夜激情欧美在线| 日韩有码中文字幕| 香蕉国产在线看| 成年人黄色毛片网站| 久久午夜亚洲精品久久| 中国美女看黄片| 俄罗斯特黄特色一大片| 国产一区二区在线av高清观看| 长腿黑丝高跟| 亚洲熟妇中文字幕五十中出| 黄片小视频在线播放| 久久精品夜夜夜夜夜久久蜜豆| 美女午夜性视频免费| 97碰自拍视频| 99国产极品粉嫩在线观看| 人人妻,人人澡人人爽秒播| 嫩草影视91久久| 999精品在线视频| 夜夜看夜夜爽夜夜摸| 亚洲人成网站在线播放欧美日韩| 精品福利观看| 国产精品国产高清国产av| 成年版毛片免费区| 在线a可以看的网站| 香蕉久久夜色| 少妇的逼水好多| 两个人看的免费小视频| 色视频www国产| 欧美成狂野欧美在线观看| 午夜免费激情av| 99国产精品一区二区蜜桃av| 亚洲中文av在线| 久久久水蜜桃国产精品网| 丰满人妻一区二区三区视频av | 可以在线观看毛片的网站| 欧美国产日韩亚洲一区| 精品久久蜜臀av无| 亚洲成人久久爱视频| 欧美性猛交╳xxx乱大交人| 黄频高清免费视频| av欧美777| 少妇熟女aⅴ在线视频| 国产精品野战在线观看| 国内久久婷婷六月综合欲色啪| 国产69精品久久久久777片 | 亚洲国产精品合色在线| 亚洲av片天天在线观看| 国产亚洲精品久久久久久毛片| 亚洲av中文字字幕乱码综合| 可以在线观看毛片的网站| 亚洲人成网站高清观看| 午夜福利成人在线免费观看| e午夜精品久久久久久久| 亚洲五月婷婷丁香| a级毛片在线看网站| 亚洲精品中文字幕一二三四区| 日韩欧美免费精品| 国产三级黄色录像| 免费高清视频大片| 1024香蕉在线观看| 亚洲激情在线av| 国产高清激情床上av| 亚洲精品一区av在线观看| 欧美日本亚洲视频在线播放| ponron亚洲| 不卡一级毛片| 精品一区二区三区视频在线观看免费| 日韩精品中文字幕看吧| 午夜日韩欧美国产| 亚洲人成网站高清观看| xxx96com| 日本在线视频免费播放| 人妻久久中文字幕网| www国产在线视频色| 97超级碰碰碰精品色视频在线观看| 国产69精品久久久久777片 | 女同久久另类99精品国产91| 亚洲精品一卡2卡三卡4卡5卡| 亚洲色图av天堂| 中文资源天堂在线| 国产熟女xx| 国产精品乱码一区二三区的特点| 在线十欧美十亚洲十日本专区| 欧美日韩精品网址| 网址你懂的国产日韩在线| 一个人看的www免费观看视频| 精品国产超薄肉色丝袜足j| 国产三级中文精品| 亚洲精品一区av在线观看| 男人舔女人下体高潮全视频| 亚洲国产精品sss在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 免费av毛片视频| 欧美性猛交黑人性爽| 亚洲欧美日韩卡通动漫| 麻豆av在线久日| tocl精华| 国产麻豆成人av免费视频| 欧美黄色淫秽网站| 久久人妻av系列| 亚洲午夜理论影院| 亚洲av五月六月丁香网| 日本成人三级电影网站| 国产精品99久久久久久久久| 国产亚洲精品综合一区在线观看| 久久伊人香网站| 国产精品亚洲av一区麻豆| 97超级碰碰碰精品色视频在线观看| 操出白浆在线播放| 嫁个100分男人电影在线观看| 亚洲精品美女久久久久99蜜臀| 国产欧美日韩一区二区三| 床上黄色一级片| 天天添夜夜摸| 制服丝袜大香蕉在线| 香蕉丝袜av| 色播亚洲综合网| 18美女黄网站色大片免费观看| 香蕉国产在线看| 老司机午夜十八禁免费视频| 亚洲一区高清亚洲精品| 久久久国产成人精品二区| 亚洲成人久久爱视频| 国产视频内射| 久久精品国产99精品国产亚洲性色| 18禁黄网站禁片午夜丰满| 欧美色视频一区免费| 国产97色在线日韩免费| 两个人看的免费小视频| 男女那种视频在线观看| 成人欧美大片| 91九色精品人成在线观看| 国产精品1区2区在线观看.| 全区人妻精品视频| 亚洲欧美日韩无卡精品| 国产主播在线观看一区二区| 欧美一区二区国产精品久久精品| 69av精品久久久久久| 小蜜桃在线观看免费完整版高清| 亚洲av五月六月丁香网| 白带黄色成豆腐渣| 九九在线视频观看精品| 91字幕亚洲| 97人妻精品一区二区三区麻豆| 欧美午夜高清在线| 亚洲狠狠婷婷综合久久图片| 男女视频在线观看网站免费| 真人一进一出gif抽搐免费| 黄色片一级片一级黄色片| 亚洲国产精品久久男人天堂| 黄色丝袜av网址大全| 美女 人体艺术 gogo| 91麻豆av在线| 蜜桃久久精品国产亚洲av| 亚洲在线观看片| 久久精品国产综合久久久| 99国产精品一区二区三区| 网址你懂的国产日韩在线| h日本视频在线播放| 亚洲国产高清在线一区二区三| www.自偷自拍.com| 999精品在线视频| 久久国产乱子伦精品免费另类| 香蕉av资源在线| 成年女人永久免费观看视频| 最近最新中文字幕大全免费视频| 午夜激情福利司机影院| 蜜桃久久精品国产亚洲av| 亚洲av成人一区二区三| 九九在线视频观看精品| 亚洲精品一区av在线观看| 天堂动漫精品| cao死你这个sao货| 亚洲国产色片| 国产精品亚洲av一区麻豆| 国产av一区在线观看免费| 久久人妻av系列| 久久精品亚洲精品国产色婷小说| 成人特级黄色片久久久久久久| 精华霜和精华液先用哪个| 国产一区二区激情短视频| 嫁个100分男人电影在线观看| 婷婷丁香在线五月| 欧美日韩国产亚洲二区| 国产一区二区三区视频了| 亚洲国产看品久久| 精品免费久久久久久久清纯| 亚洲精品一区av在线观看| 丰满人妻一区二区三区视频av | 老熟妇乱子伦视频在线观看| 午夜亚洲福利在线播放| 女警被强在线播放| 欧美性猛交黑人性爽| 欧美日韩福利视频一区二区| 老鸭窝网址在线观看| 日韩有码中文字幕| 99久久成人亚洲精品观看| 国产高清三级在线| 手机成人av网站| 男人舔女人的私密视频| 日韩欧美三级三区| 色播亚洲综合网| 性欧美人与动物交配| 亚洲人成网站高清观看| 婷婷亚洲欧美| 国产野战对白在线观看| 国产亚洲欧美在线一区二区| 日本 欧美在线| 午夜激情欧美在线| 国产在线精品亚洲第一网站| 国产免费男女视频| 午夜福利高清视频| 久久久成人免费电影| 国产成人精品无人区| 国产成人欧美在线观看| 免费在线观看日本一区| 麻豆成人av在线观看| 搡老岳熟女国产| 免费大片18禁| 午夜激情福利司机影院| 亚洲国产中文字幕在线视频| 中文资源天堂在线| 91字幕亚洲| 日韩欧美一区二区三区在线观看| 91在线精品国自产拍蜜月 | 午夜亚洲福利在线播放| 欧美色欧美亚洲另类二区| 91老司机精品| 亚洲午夜理论影院| 一区二区三区激情视频| 久久精品aⅴ一区二区三区四区| 日韩欧美免费精品| 99久久精品国产亚洲精品| 真人做人爱边吃奶动态| 亚洲成人久久性| 免费av毛片视频| 久久久久久久精品吃奶| 国产精品亚洲av一区麻豆| 窝窝影院91人妻| 亚洲成人免费电影在线观看| 久久婷婷人人爽人人干人人爱| 国产精品精品国产色婷婷| 一进一出抽搐动态| 亚洲精品在线观看二区| 午夜福利免费观看在线| www.www免费av| 精品日产1卡2卡| 麻豆一二三区av精品| 老汉色∧v一级毛片| 宅男免费午夜| 午夜免费激情av| 99热这里只有精品一区 | 欧美日韩黄片免| 久久国产乱子伦精品免费另类| 国产午夜精品久久久久久| 女人被狂操c到高潮| 成熟少妇高潮喷水视频| 国产aⅴ精品一区二区三区波| 成人国产综合亚洲| 亚洲国产精品999在线| 五月玫瑰六月丁香| 夜夜爽天天搞| 亚洲avbb在线观看| 国产精品乱码一区二三区的特点| 少妇人妻一区二区三区视频| a级毛片a级免费在线| 精品无人区乱码1区二区| 午夜日韩欧美国产| 91字幕亚洲| 国产精品九九99| 熟女人妻精品中文字幕| 久久久久久久久中文| 午夜影院日韩av| 国产精品精品国产色婷婷| 美女 人体艺术 gogo| 搞女人的毛片| av天堂在线播放| 国产成人福利小说| 中文字幕精品亚洲无线码一区| 亚洲国产欧美网| 婷婷精品国产亚洲av在线| 国产精品自产拍在线观看55亚洲| 国产单亲对白刺激| 白带黄色成豆腐渣| 国产91精品成人一区二区三区| 国产亚洲av嫩草精品影院| 精品久久久久久成人av| 婷婷亚洲欧美| 日韩欧美一区二区三区在线观看| 丝袜人妻中文字幕| 看免费av毛片| 国产男靠女视频免费网站| 亚洲欧美日韩卡通动漫| 99久久成人亚洲精品观看| 老汉色∧v一级毛片| 国产黄a三级三级三级人| 美女 人体艺术 gogo| 午夜日韩欧美国产| 亚洲乱码一区二区免费版| 免费人成视频x8x8入口观看| 精品一区二区三区四区五区乱码| 亚洲国产日韩欧美精品在线观看 | 成人一区二区视频在线观看| 女同久久另类99精品国产91| 午夜免费观看网址| 亚洲欧美日韩高清专用| 国产精品99久久99久久久不卡| 久久久水蜜桃国产精品网| tocl精华| 欧美一级a爱片免费观看看| 男人舔女人的私密视频| 极品教师在线免费播放| 久久精品国产99精品国产亚洲性色| 成年版毛片免费区| 美女黄网站色视频| 国产美女午夜福利| 国产亚洲精品综合一区在线观看| 欧美在线黄色| 桃红色精品国产亚洲av| 脱女人内裤的视频| 午夜福利在线在线| 国产高清有码在线观看视频| 99久久成人亚洲精品观看| 日日摸夜夜添夜夜添小说| 12—13女人毛片做爰片一| 日本精品一区二区三区蜜桃| 美女 人体艺术 gogo| 无限看片的www在线观看| 一a级毛片在线观看| 免费观看的影片在线观看| 在线国产一区二区在线| 亚洲国产精品999在线| 国产爱豆传媒在线观看| 亚洲国产看品久久| 舔av片在线| 人妻久久中文字幕网| 精品福利观看| 哪里可以看免费的av片| 18禁裸乳无遮挡免费网站照片| 国产精品久久电影中文字幕| 亚洲人成网站在线播放欧美日韩| 色在线成人网| 香蕉丝袜av| 久久热在线av| 日韩人妻高清精品专区| 国产高清有码在线观看视频| 黑人欧美特级aaaaaa片| 神马国产精品三级电影在线观看| 亚洲国产精品成人综合色| 淫妇啪啪啪对白视频| 此物有八面人人有两片| 免费人成视频x8x8入口观看| www.熟女人妻精品国产| 久久精品人妻少妇| 啪啪无遮挡十八禁网站| 真人做人爱边吃奶动态| 久久精品夜夜夜夜夜久久蜜豆| av女优亚洲男人天堂 | 中文字幕人妻丝袜一区二区| 中文字幕精品亚洲无线码一区| 老司机在亚洲福利影院| 男女床上黄色一级片免费看| 国内精品久久久久久久电影| 欧美日韩瑟瑟在线播放| 男人和女人高潮做爰伦理| 亚洲一区二区三区不卡视频| 久久天堂一区二区三区四区| 久久久国产成人免费| 一个人看视频在线观看www免费 | 欧美日本亚洲视频在线播放| 欧美成人一区二区免费高清观看 | 黄色女人牲交| 99热6这里只有精品| 老司机在亚洲福利影院| 亚洲,欧美精品.| 热99re8久久精品国产| 白带黄色成豆腐渣| а√天堂www在线а√下载| 中文字幕高清在线视频| 中文字幕人成人乱码亚洲影| 精品一区二区三区视频在线 | 窝窝影院91人妻| 精华霜和精华液先用哪个| 真人一进一出gif抽搐免费| 亚洲自偷自拍图片 自拍| 亚洲精品美女久久久久99蜜臀| 亚洲成a人片在线一区二区| 国产黄片美女视频| 天堂√8在线中文| 精品99又大又爽又粗少妇毛片 | 国产伦在线观看视频一区| 丰满的人妻完整版| 久久这里只有精品19| 国产毛片a区久久久久| 国产成人福利小说| 久久伊人香网站| 欧美又色又爽又黄视频| 色播亚洲综合网| 好看av亚洲va欧美ⅴa在| 国产成人福利小说| 久久午夜综合久久蜜桃| 欧美日本亚洲视频在线播放| 熟女人妻精品中文字幕| 成人特级黄色片久久久久久久| 久久久久久久精品吃奶| 亚洲一区二区三区不卡视频| 久久香蕉国产精品| 久久国产精品人妻蜜桃| 热99在线观看视频| 亚洲精品456在线播放app | 精品人妻1区二区| 午夜精品久久久久久毛片777| 女同久久另类99精品国产91| 日本a在线网址| 欧美日本视频| 国产精品日韩av在线免费观看| 真实男女啪啪啪动态图| 国产伦精品一区二区三区视频9 | 精品免费久久久久久久清纯| netflix在线观看网站| 久久香蕉国产精品| 国产在线精品亚洲第一网站| 国产精品久久久久久精品电影| 成人三级做爰电影| 欧美av亚洲av综合av国产av| 女警被强在线播放| 午夜视频精品福利| 欧美xxxx黑人xx丫x性爽| 成人国产综合亚洲| 制服丝袜大香蕉在线| 亚洲av片天天在线观看| 女生性感内裤真人,穿戴方法视频| 淫妇啪啪啪对白视频| 国产精品免费一区二区三区在线| 蜜桃久久精品国产亚洲av| 国产一区二区在线观看日韩 | 最近视频中文字幕2019在线8| 欧美一级a爱片免费观看看| 美女黄网站色视频| 国产aⅴ精品一区二区三区波| 99久久精品热视频| 亚洲男人的天堂狠狠| 国产真人三级小视频在线观看| 老司机福利观看| www日本黄色视频网| 国产v大片淫在线免费观看| 欧美av亚洲av综合av国产av| e午夜精品久久久久久久| 神马国产精品三级电影在线观看| 国产真实乱freesex| 黄色 视频免费看| 国产亚洲精品久久久com| 久久精品91蜜桃| 欧美乱妇无乱码| 国产极品精品免费视频能看的| 午夜福利在线观看吧| h日本视频在线播放| 99久久无色码亚洲精品果冻| 老熟妇乱子伦视频在线观看| 欧美黑人巨大hd| 亚洲专区国产一区二区| 午夜福利免费观看在线| 一级毛片高清免费大全| 97超视频在线观看视频| 亚洲人成伊人成综合网2020| 亚洲熟女毛片儿| 国产免费av片在线观看野外av| 男人的好看免费观看在线视频| 少妇熟女aⅴ在线视频| 国产黄色小视频在线观看| 日本在线视频免费播放| 少妇丰满av| 国产精品久久电影中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 99在线视频只有这里精品首页| 成人特级黄色片久久久久久久| 偷拍熟女少妇极品色| 欧美日韩国产亚洲二区| 亚洲精品国产精品久久久不卡| 免费在线观看日本一区| 免费高清视频大片| 国产综合懂色| 美女高潮喷水抽搐中文字幕| 免费在线观看日本一区| 久久天堂一区二区三区四区| 国产人伦9x9x在线观看| www日本黄色视频网| 色av中文字幕| 日本与韩国留学比较| 美女高潮的动态|