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

    福建省茶樹品種資源初步關聯(lián)分析

    2017-11-23 08:03:03孔祥瑞王讓劍鄭國華
    茶葉學報 2017年3期
    關鍵詞:烏龍茶農(nóng)藝綠茶

    楊 軍,孔祥瑞,王讓劍,鄭國華

    (福建省農(nóng)業(yè)科學院茶葉研究所/國家茶樹改良中心福建分中心,福建 福安 355015)

    福建省茶樹品種資源初步關聯(lián)分析

    楊 軍,孔祥瑞,王讓劍,鄭國華

    (福建省農(nóng)業(yè)科學院茶葉研究所/國家茶樹改良中心福建分中心,福建 福安 355015)

    為挖掘茶樹農(nóng)藝性狀相關的分子標記位點,選用38對SSR引物對43份福建茶樹品種進行多態(tài)性掃描與群體結構分析。在群體遺傳結構分析的基礎上采用Tassel GLM(general linear model)和 MLM(mixed linear model)方法進行標記與農(nóng)藝性狀的關聯(lián)分析。共檢測出238個等位變異,每個位點6.26個等位變異。供試群體的 Shannon 指數(shù)平均值為1.261。通過群體遺傳結構分析將供試材料劃分成3個亞群。以 GLM 分析,發(fā)現(xiàn)8個與種徑、茶類、百粒重和花瓣數(shù)相關聯(lián)的標記,對表型變異的解釋率分別54.3%、13.3%~21.8%、57.2%和50.9%;以MLM分析,發(fā)現(xiàn) 6 個與種徑、茶類和百粒重相關的標記,各標記對表型變異的解釋率分別為54.2%、19.7%~24.2%、31.1%。

    茶;SSR;群體遺傳結構;關聯(lián)分析

    茶樹是多年生常異花授粉植物,構建群體進行QTL定位難度較大,且時間較長。關聯(lián)分析僅需以自然群體為材料,通過分子標記進行全基因組多態(tài)性掃描,結合目標性狀調查,就可以發(fā)掘相關聯(lián)的優(yōu)異等位基因。因此,關聯(lián)分析方法比較適合茶樹上目標性狀等位基因的挖掘與輔助育種。關聯(lián)分析的材料一般選擇代表性強,來源廣泛,其中有選擇芝麻、可可核心種質[1-2]、秈稻微核心種質[3]、全國范圍內具代表性的豇豆種質[4]、國內外大麥品種[5]等進行關聯(lián)分析。同時,近年關聯(lián)分析在植物的抗性、農(nóng)藝性狀、產(chǎn)量、品質等性狀上有較多研究,檢測到煙草赤星病抗性[6]、水稻紋枯病[7]、青稞的穗粒數(shù)、株高等[8]、向日葵分枝性狀[9]、玉米穗行數(shù)[10]、陸地棉葉枝數(shù)[11]、大麥的株高、穗長等[12]、半野生棉仁含油量[13]等顯著關聯(lián)的分子標記。目前使用SSR標記對茶樹農(nóng)藝性狀進行關聯(lián)分析研究的報道相對較少,本研究對43份福建審(認)定茶樹品種材料進行關聯(lián)分析,篩選描述茶樹品種特征特性密切相關的性狀,花、果、常規(guī)生化成分等性狀相關的分子標記,為后期茶樹等位基因發(fā)掘與分子標記輔助育種提供參考。

    1 材料與方法

    1.1 試驗材料

    試驗材料來源于福建省農(nóng)業(yè)科學院茶葉研究所,福建省審(認)定茶樹品種43個,采集各個材料一芽二葉鮮葉保存在-70℃冰箱中備用。

    1.2 試驗方法

    1.2.1 基因組DNA的提取 采用改進的CTAB法[14]提取茶樹基因組DNA。用1.0%的瓊脂糖凝膠電泳進行茶樹基因組分子量大小檢測,用756-MC型紫外分光光度計測定茶樹基因組DNA的純度(OD260/280比值)。

    1.2.2 農(nóng)藝性狀調查數(shù)據(jù) 農(nóng)藝性狀來源于《福建省茶樹品種圖志》[15]調查數(shù)據(jù),詳細見表1。

    1.2.3 引物合成 參照文獻[16-17]引物序列,序列由上海Sangon公司合成。

    1.2.4 PCR擴增和產(chǎn)物鑒定 PCR反應體系為:ddH2O 18.8 μL,10×Buffer 2.5 μL(Mg2+),dNTP(10 mmol·L-1)0.5 μL,上、下游primer(10 μmol·L-1)各0.5 μL,Taq polymerase 0.2 μL(0.5 U),模板DNA 1 μL。PCR反應于美國ABI-9600型擴增儀上進行,熱循環(huán)程序為:94℃預變性4 min,使模板DNA 充分變性,然后進入下列溫度循環(huán):94℃變性45 s,不同溫度條件退火60 s,72℃延伸75 s,重復35個熱循環(huán);72℃延伸10 min,最后4℃保存。

    引物統(tǒng)一在反向引物(R)的5′段標記熒光(FAM/TAMRA),由上海百力格生物科技有限公司合成。擴增產(chǎn)物0.5 μL,GeneScanTM500 ROXTM0.5 μL,HiDi 9 μL,混勻后使用ABI公司3730XL進行毛細管電泳。

    1.2.5 數(shù)據(jù)處理 采用ABI公司的GeneMapper 4.0軟件,選擇GeneScanTM500 ROXTM作為分子量標準,對每一個擴增出的條帶記錄大小。運用PopGen3.2軟件計算Shannon信息指數(shù)、觀測雜合度、期望雜合度。應用PIC_CALC引物的PIC值。方差分析采用SAS。

    1.2.6 關聯(lián)分析 使用Struture 2.3.4[18]進行群體結構分析,先計算各材料歸屬第k亞群的概率Q值。使用混合模型K值從2設置到9模擬群體結構,每一個K值重復6次。將MCMC迭代設為100000次,初始burn-in次數(shù)為100000。依據(jù)似然值最大原則選取合適的K值為群體數(shù)目,計算Q參數(shù),將其作為協(xié)變量,并CLUMPP對多次運行結果的進行整合。以SPAGeDi-1.3d處理基因型數(shù)據(jù)獲得個體間親緣關系Kinship矩陣[19]。

    采用Tassel 3軟件一般線性模型(GLM)和混合線性模型(MLM)兩種程序進行關聯(lián)分析。GLM分析中,以各親本材料的對應Q值作為協(xié)變量,將SSR標記與花冠直徑、百粒重、氨基酸、種徑等表型性狀進行回歸分析,尋找與之相關聯(lián)的標記,并確定其解釋率;MLM分析中,采用Q+K方法分析,性狀與GLM分析相同,確定關聯(lián)位點,并計算標記對表型變異的解釋率。

    表1 茶樹品種農(nóng)藝性狀數(shù)據(jù)

    (續(xù)表1)

    品種名稱花冠直徑(cm)花瓣數(shù)(片)果實直徑(cm)果皮厚度(cm)一芽三葉重(g)種徑(cm)百粒重(g)茶多酚(%)氨基酸(%)咖啡堿(%)水浸出物(%)茶類大葉烏龍4.662.200.06075.01.6276.2017.54.23.448.3烏龍茶福云20號4.472.150.06796.51.36135.8018.83.93.450.6綠茶春蘭4.271.920.06758.01.41102.4015.63.93.751.4烏龍茶鳳圓春3.372.080.096148.81.4963.4014.55.23.742.7烏龍茶鐵觀音3.261.680.07060.51.25123.9417.44.73.751.0烏龍茶黃觀音3.962.140.08858.01.28112.2019.44.83.448.4烏龍茶杏仁茶3.262.030.050139.61.39139.5015.64.73.745.4烏龍茶紫玫瑰3.262.020.04062.01.31110.7016.34.53.148.4烏龍茶九龍袍4.161.660.07083.01.25126.5818.84.13.249.9烏龍茶福安大白茶3.87NANA98.0NANA15.56.13.451.3綠茶黃奇4.462.120.07165.01.2287.9019.64.24.050.2烏龍茶歌樂茶3.872.080.13180.01.6170.0013.55.73.646.2綠茶黃棪2.9562.100.07059.01.4895.8016.23.53.648.0烏龍茶紫牡丹4.062.020.04054.01.31110.7018.43.94.348.6烏龍茶丹桂4.261.940.06866.01.2390.0017.73.33.249.9烏龍茶大紅袍3.462.350.070NA1.18101.1617.15.03.551.0烏龍茶福云6號3.562.180.08369.01.32108.5814.94.72.945.1綠茶梅占4.161.800.080103.01.35105.3016.54.13.951.7烏龍茶福鼎大毫茶4.7571.800.090104.01.5190.5017.35.33.247.2綠茶福鼎大白茶3.8581.630.06063.01.33118.5014.84.03.349.8綠茶金萱4.072.080.079NA1.4990.5019.24.43.349.9烏龍茶福建水仙4.0571.530.150112.01.47109.4417.63.34.050.5烏龍茶毛蟹4.462.110.10068.51.44119.8614.74.23.248.2烏龍茶金觀音3.672.000.06050.01.51163.6219.04.43.845.6烏龍茶福云5953.472.320.060111.01.6995.4012.54.74.747.3綠茶綠芽佛手4.082.900.096147.01.78100.5016.23.13.149.0烏龍茶肉桂371.860.07053.01.1475.8817.73.83.152.3烏龍茶悅茗香4.161.900.07060.01.26125.8221.43.63.949.4烏龍茶政和大白茶4.87NANA123.0NANA13.55.93.346.8綠茶白芽奇蘭4.071.810.118139.01.2852.0016.43.63.948.2烏龍茶榕春早3.571.850.063NA1.0690.0017.74.63.245.9綠茶春閨3.072.510.090NA1.3575.8013.34.23.841.4烏龍茶

    注:NA代表數(shù)據(jù)缺失。

    2 結果與分析

    2.1 引物多態(tài)性與等位位點分析

    從表1中看出,38對SSR標記在參試材料中共檢測到等位位點238個,平均每對引物檢測到6.26個,其中最多為SSR2、SSR16為13個,多態(tài)性十分豐富;38對引物的PIC值在0.255~0.841之間,平均為0.580,篩選出來的引物具有較強鑒定能力;SSR引物的Shannon指數(shù)變幅為0.506~2.158,平均值為1.261,參試茶樹種質材料群體間的遺傳變異較大,具有較為豐富的遺傳多樣性。計算觀測雜合度與期望雜合度平均值比較接近,說明篩選的引物檢測到的等位基因在參試材料中分布相對均勻。

    2.2 群體遺傳結構分析

    利用Structure 2.3.1軟件,基于數(shù)學模型分析由43份材料構成的參試材料的遺傳結構,發(fā)現(xiàn)樣本的等位變異頻率特征類型數(shù)K呈持續(xù)增大趨勢,當K=4時出現(xiàn)峰值,并呈現(xiàn)明顯的拐點(圖1)。由此將43份茶樹親本分為3個類群(圖2),分別包含9、11和23份。該自然群體結構比較簡單,能夠有效地降低群體結構對關聯(lián)分析產(chǎn)生的影響。

    表2 38對引物的擴增結果

    圖1 Structure分析中△K隨K值的變化Fig.1 △K value from structure analysis

    2.3 SSR標記與部分農(nóng)藝性狀的關聯(lián)分析

    GLM分析結果顯示,在所檢測的38個標記中有8個與種徑、茶類、花瓣數(shù)、百粒重顯著相關(P<0.01)。同時,將P值進行蒙特卡羅模擬顯著性檢測(防止偽關聯(lián)),SSR14、SSR15、SSR16、SSR19、SSR26、SSR37達到顯著性標準(P<0.05),其中5個與茶類相關,1個與種徑相關。各標記對表型變異的解釋率為13.3%~54.3%,解釋率最大的標記是SSR37為54.3%,與種徑相關(見表3)。

    MLM分析結果顯示,在所檢測的38個標記中有6個與種徑、茶類、百粒重數(shù)顯著相關(P<0.05),其中4個與茶類相關,1個與種徑相關,1個與百粒重相關。各標記對表型變異的解釋率為19.7%~54.3%。MLM分析檢測到與農(nóng)藝性狀相關的SSR比GLM分析結果少2個,其中有5個標記與GLM 分析得到的相同,但解釋率均有所上升,但P值標準變大,準確性不如GLM。

    3 討論與結論

    3.1 GLM與MLM分析結果相比較

    本試驗GLM模型關聯(lián)顯著性選擇的P<0.01,而MLM模型關聯(lián)顯著性選擇的P<0.05,與前人研究相一致[20-21]。從P值選擇標準來看,是GLM模型分析結果相對更準確,但李曉娜等[13]認為,GLM模型以各親本材料的對應Q值作為協(xié)變量,MLM模型采用Q+K方法分析,不僅考慮群體結構Q的影響,還考慮兩兩材料間的親緣關系(Kinship)系數(shù),減少“偽關聯(lián)”的數(shù)目,分析獲得的結果更加準確可靠。5對標記在MLM模型與GLM模型中都檢測達到顯著差異,因此,這5對標記的關聯(lián)分析結果可靠性相對較高。單一模型中達到顯著差異花瓣數(shù)、百粒重性狀,需要進一步檢測。

    圖2 K=3時,Structure的運算結果Fig. 2 Calculated results of structure(K=3)

    表3 與農(nóng)藝性狀顯著(P<0.05)相關的位點及其對表型變異的解釋率

    3.2 關聯(lián)分析結果

    本試驗選擇的目標性狀主要參考茶樹品種志中對茶樹品種特征特性的描述,其中,氨基酸、咖啡堿等生化成分未檢測到顯著相關的分子標記,推測是參試材料較少、引物數(shù)量不足與常規(guī)生化成分受年份、取樣等因素影響較大。茶樹品質相關基因的挖掘一直是茶樹輔助育種與資源鑒定的重點,尤其是烏龍茶香氣、滋味相關品質基因。烏龍茶香氣、滋味品質是在加工過程中形成,目前還沒有明確受主要影響因子控制,推測是多基因共同起作用的。烏龍茶與綠茶間最大的品質特征差異是香氣、滋味品質,茶類性狀有可能與烏龍茶滋味、香氣品質相關,參考前人研究對福建省茶樹品種進行烏龍茶與綠茶分類,利用分子標記進行掃描與關聯(lián)分析。從中篩選到5對與茶類顯著相關的標記,連鎖群參考Ma Jian-qiang[16]等的研究結果,5對標記位于不同的連鎖群,初步說明茶類性狀是多位點,多基因共同控制,與實驗預期一致。以這些標記作為組合,可以初步鑒定福建野生資源、地方資源的遺傳組成偏向綠茶或烏龍茶,進行雜交親本輔助選擇。同時,進行綠茶與烏龍茶雜交時,在幼齡期具有綠茶表型性狀群體中,選擇具有烏龍茶遺傳組成的后代,可以初步輔助選擇高香綠茶新品系。

    致謝:感謝試驗過程中陳常頌研究員等人的支持。

    [1]王蕾,黎冬華,齊小瓊,等.芝麻核心種質芝麻素和芝麻酚林的關聯(lián)分析[J].中國油料作物學報,2014,36(1):32-37.

    [2]李付鵬,秦曉威,郝朝運,等.可可核心種質遺傳多樣性及果實性狀與 SSR 標記關聯(lián)分析[J].熱帶作物學報,2016,37(2):226-233.

    [3]邱先進,袁志華,陳凱,等.用全基因組關聯(lián)分析解析秈稻堊白的遺傳基礎[J].作物學報,2015,41(7):1007-1016.

    [4]潘磊,李依,余曉露,等.豇豆產(chǎn)量性狀與 SSR 分子標記的關聯(lián)分析[J].湖北農(nóng)業(yè)科學,2015,54(16):3952-3962.

    [5]司二靜,張宇,汪軍成,等.大麥農(nóng)藝性狀與 SSR 標記的關聯(lián)分析[J].作物學報,2015,41(7):1064-1072.

    [6]朱承廣,任民,蔣彩虹,等.以關聯(lián)分析發(fā)掘煙草抗赤星病基因分子標記[J].中國煙草科學,2017,38(1):68-72.

    [7]孫曉棠,盧冬冬,歐陽林娟,等.水稻紋枯病抗性關聯(lián)分析及抗性等位變異發(fā)掘[J].作物學報,2014,40(5):779-787.

    [8]孟亞雄,孟林祎,汪軍成,等.青稞遺傳多樣性及其農(nóng)藝性狀與 SSR 標記的關聯(lián)分析[J].作物學報,2016,42(2):180-189.

    [9]劉勝利,段維,王鵬,等.向日葵分枝和株高性狀的關聯(lián)分析定位研究[J].海南熱帶海洋學院學報,2017,24(2):80-84.

    [10]張煥欣,翁建峰,張曉聰,等.玉米穗行數(shù)全基因組關聯(lián)分析[J].作物學報,2014,40(1):1-6.

    [11]梁冰,范術麗,宋美珍,等.陸地棉農(nóng)藝性狀與SSR標記的關聯(lián)分析[J].棉花學報,2014,26(5):387-395.

    [12]賴勇,王鵬喜,范貴強,等.大麥SSR標記遺傳多樣性及其與農(nóng)藝性狀關聯(lián)分析[J].中國農(nóng)業(yè)科學2013,46(2):233-242.

    [13]李曉娜,蔡小彥,Kiflom Weldu OkubazghP,等.半野生棉棉仁含油量與SSR標記的關聯(lián)分析[J].棉花學報,2015,27(2):104-110.

    [14]金基強,崔海瑞,龔曉春,等.用EST-SSR標記對茶樹種質資源的研究[J].遺傳,2007,29(1):103-108.

    [15]陳常頌、余文權、王慶森等等.福建省茶樹品種圖志[M].北京:中國農(nóng)業(yè)科學技術出版社,2016:16-104.

    [16]Ma Jian-qiang,Yao Ming-zhe,Ma Chun-lei,etal. Construction of a SSR-based genetic map and identification of QTLs for catechins content in tea plant(Camelliasinensis)[J].PLOS ONE,2014,9(3):1-11.

    [17]段云裳,姜燕華,王麗鴛,等.中國紅、綠茶適制品種(系)遺傳多樣性與親緣關系的SSR分析[J].中國農(nóng)業(yè)科學,2011,44(1):99-109.

    [18]Prichard J K, Stepphen M, Donnelly P. Inference of population structure using multilocus genotype data[J]. Genetics, 2000(155):945-959.

    [19]Earl D A. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method[J]. Conserv Genet Resour, 2012(4):359-361.

    [20]郭志軍,趙云雷,陳偉,等.陸地棉SSR標記遺傳多樣性及其與農(nóng)藝性狀的關聯(lián)分析[J].棉花學報,2014,26(5):420-430.

    [21]馮英娜,柳李旺,劉衛(wèi)東,等.茄子SSR遺傳多樣性及其農(nóng)藝性狀的關聯(lián)分析[J].江蘇農(nóng)業(yè)學報,2014,30(4):839-847.

    APreliminaryReportonAssociationAnalysisofTeaVarietiesinFujian

    YANG Jun, WANG Rang-jian, KONG Xiang-rui, ZHENG Guo-hua

    (TeaResearchInstitute,FujianAcademyofAgriculturalSciences/FujianBranch,NationalCenterforTeaImprovement,Fu’an,Fujian355015,China)

    To identify the molecular marker loci related to specific tea agronomic traits, 38 simple sequence repeat (SSR) markers were selected for the polymorphism detection and population structure analysis on 43 tea varieties found in Fujian. The population structure was analyzed based on the polymorphic SSR markers, while the association between the markers and agronomic traits was determined by fitting the data to TASSEL general linear model (GLM) and mixed linear model (MLM). A total of 238 alleles were detected in the 43 tea varieties with an average count of 6.26. The averaged Shannon’s index was 1.261. The structure analysis showed 3 genetic subpopulations for the varieties tested. The association of markers and traits based on GLM showed 8 SSR markers relating to the seed stem, tea type, hundred-grain weight and flower petal count of a plant, with the phenotypic contributions of 54.3%, 13.3%~21.8%, 57.2%, and 50.9%, respectively. Whereas, MLM indicated that 6 SSR markers closely related to the seed stem, tea type, and hundred-grain weight, with the phenotypic contributions of 54.2%, 19.7%~24.2% and 31.1%, respectively.

    tea plant; SSR; population structure; association analysis

    S571.1

    A

    2096-0220(2017)03-0115-06

    2017-06-27初稿;2017-08-11修改稿

    福建省現(xiàn)代農(nóng)業(yè)(茶葉)產(chǎn)業(yè)技術體系建設項目(2014~2017);福建省屬公益類科研院所基本科研專項(2015R1012-7、2017R1012-6);福建省自然科學基金項目(2015J01098);福建省農(nóng)業(yè)科學院茶樹種質資源創(chuàng)新團隊項目(STIT 2017-3-12)。

    楊軍(1981-),男,助理研究員,從事茶樹種質資源與遺傳育種研究。E-mail:yangjun007@qq.com

    猜你喜歡
    烏龍茶農(nóng)藝綠茶
    THE PRICE OF BEAUTY
    漢語世界(2022年3期)2022-06-19 05:49:22
    The Price of Beauty
    綠茶或能預防和延緩阿爾茨海默病
    農(nóng)機需要農(nóng)藝“標準”,農(nóng)藝需要農(nóng)機“靈活”——2021國際農(nóng)機展不容錯過的果蔬茶機械化論壇
    落葉果樹(2021年6期)2021-02-12 01:29:20
    水稻主要農(nóng)藝性狀與產(chǎn)量的相關性分析
    烏龍茶秋季采摘時期調節(jié)技術研究
    14份蔞蒿種質資源主要農(nóng)藝性狀及營養(yǎng)成分評價
    中國蔬菜(2016年8期)2017-01-15 14:23:37
    近年烏龍茶做青工藝的發(fā)展研究
    烏龍茶連續(xù)自動化初制生產(chǎn)線的現(xiàn)狀與展望
    90團舉辦初級農(nóng)藝工培訓班
    一个人看视频在线观看www免费| 亚洲伊人久久精品综合| 日韩av免费高清视频| 最后的刺客免费高清国语| 欧美成人精品欧美一级黄| 日日啪夜夜爽| 99re6热这里在线精品视频| 熟女电影av网| 国产成人免费观看mmmm| 99九九在线精品视频| 人妻系列 视频| 一级毛片aaaaaa免费看小| 精品国产露脸久久av麻豆| 少妇丰满av| 国产成人av激情在线播放 | 久久av网站| 日本免费在线观看一区| 亚州av有码| 18禁在线无遮挡免费观看视频| 国产精品久久久久久久久免| 亚洲久久久国产精品| 亚洲一区二区三区欧美精品| 9色porny在线观看| 日韩一区二区三区影片| 久久久久久久久久成人| 少妇被粗大猛烈的视频| 亚洲,一卡二卡三卡| 少妇高潮的动态图| av天堂久久9| 99久久综合免费| 日韩欧美精品免费久久| 高清av免费在线| 亚洲精品自拍成人| 高清av免费在线| 日韩免费高清中文字幕av| 欧美xxxx性猛交bbbb| 国产乱人偷精品视频| 80岁老熟妇乱子伦牲交| 肉色欧美久久久久久久蜜桃| 搡女人真爽免费视频火全软件| 国产精品一区二区在线不卡| 久久久久精品性色| 久久久精品免费免费高清| 欧美日韩综合久久久久久| 欧美日韩精品成人综合77777| 伦理电影免费视频| 三级国产精品片| 亚洲精品久久成人aⅴ小说 | 午夜福利视频精品| 久久久久久久国产电影| 黄色一级大片看看| 一区二区av电影网| 日韩欧美精品免费久久| 国产免费又黄又爽又色| 在线看a的网站| 七月丁香在线播放| 99久久精品国产国产毛片| 亚洲精品一区蜜桃| 久久国产精品大桥未久av| 久久精品国产自在天天线| 人人澡人人妻人| 国产成人aa在线观看| 精品久久久噜噜| 全区人妻精品视频| 亚洲国产av影院在线观看| 天美传媒精品一区二区| 欧美 日韩 精品 国产| 欧美人与善性xxx| 一区二区三区精品91| 国产亚洲最大av| 欧美精品亚洲一区二区| 亚洲精华国产精华液的使用体验| 久久99蜜桃精品久久| 成年美女黄网站色视频大全免费 | 午夜91福利影院| 777米奇影视久久| 只有这里有精品99| 亚洲欧洲精品一区二区精品久久久 | 欧美国产精品一级二级三级| 久久人人爽人人片av| 亚洲精品久久成人aⅴ小说 | 国语对白做爰xxxⅹ性视频网站| 内地一区二区视频在线| 99精国产麻豆久久婷婷| 亚州av有码| 精品久久久精品久久久| 又黄又爽又刺激的免费视频.| 精品国产一区二区三区久久久樱花| www.色视频.com| 美女内射精品一级片tv| 熟妇人妻不卡中文字幕| 成人亚洲欧美一区二区av| 草草在线视频免费看| 色哟哟·www| 国产男人的电影天堂91| 国产老妇伦熟女老妇高清| 国产日韩欧美在线精品| 永久免费av网站大全| 精品午夜福利在线看| 美女cb高潮喷水在线观看| 久久久a久久爽久久v久久| 夫妻性生交免费视频一级片| 国产欧美另类精品又又久久亚洲欧美| 男人爽女人下面视频在线观看| 日本欧美视频一区| 日韩av在线免费看完整版不卡| 多毛熟女@视频| 亚洲人与动物交配视频| 亚洲久久久国产精品| 欧美日韩成人在线一区二区| 午夜福利视频精品| 狠狠精品人妻久久久久久综合| 欧美激情极品国产一区二区三区 | 啦啦啦视频在线资源免费观看| 中文字幕亚洲精品专区| 九九在线视频观看精品| 午夜久久久在线观看| 丝袜脚勾引网站| 考比视频在线观看| 精品国产乱码久久久久久小说| 免费观看无遮挡的男女| 日本91视频免费播放| 久久人妻熟女aⅴ| 少妇的逼好多水| 国产精品一区www在线观看| 久久毛片免费看一区二区三区| 亚洲av国产av综合av卡| 91久久精品国产一区二区成人| 亚洲精品一区蜜桃| 久久免费观看电影| 国产成人av激情在线播放 | 国产av码专区亚洲av| 一级毛片 在线播放| 综合色丁香网| 久久久亚洲精品成人影院| 国产av码专区亚洲av| 欧美精品一区二区免费开放| 黑人高潮一二区| 狂野欧美激情性xxxx在线观看| 最近中文字幕高清免费大全6| 高清欧美精品videossex| 免费黄色在线免费观看| 飞空精品影院首页| 一级毛片aaaaaa免费看小| 中国美白少妇内射xxxbb| 美女脱内裤让男人舔精品视频| 欧美日韩在线观看h| 亚洲精品乱久久久久久| xxxhd国产人妻xxx| 午夜老司机福利剧场| 亚洲成人手机| 夜夜骑夜夜射夜夜干| 国产乱来视频区| 亚洲精品一区蜜桃| av有码第一页| 亚洲国产精品一区二区三区在线| 交换朋友夫妻互换小说| 日本与韩国留学比较| 精品久久久久久电影网| 国产伦精品一区二区三区视频9| 国产成人免费观看mmmm| 日本猛色少妇xxxxx猛交久久| 国产免费一区二区三区四区乱码| a 毛片基地| 三上悠亚av全集在线观看| 高清在线视频一区二区三区| 亚洲av成人精品一区久久| 精品亚洲乱码少妇综合久久| 亚洲av在线观看美女高潮| 黑人巨大精品欧美一区二区蜜桃 | 欧美日韩视频高清一区二区三区二| 国精品久久久久久国模美| 视频中文字幕在线观看| 黄色毛片三级朝国网站| 一本一本综合久久| 国产精品久久久久久久久免| 亚洲国产精品国产精品| 国产精品蜜桃在线观看| 99热这里只有是精品在线观看| 国语对白做爰xxxⅹ性视频网站| 91精品伊人久久大香线蕉| 2021少妇久久久久久久久久久| 久久99热6这里只有精品| 蜜桃国产av成人99| 日韩电影二区| 精品国产露脸久久av麻豆| 亚洲天堂av无毛| 91久久精品国产一区二区三区| 美女cb高潮喷水在线观看| 18禁观看日本| 午夜福利,免费看| 亚洲精品乱码久久久久久按摩| 99久久中文字幕三级久久日本| 免费观看无遮挡的男女| 亚洲成人一二三区av| 久久久久久久亚洲中文字幕| 成人免费观看视频高清| 日本与韩国留学比较| 91久久精品国产一区二区三区| 亚洲欧美中文字幕日韩二区| 国产亚洲精品第一综合不卡 | 免费大片黄手机在线观看| 亚洲精品视频女| 国产男女内射视频| 亚洲精品国产av成人精品| 亚洲美女搞黄在线观看| www.av在线官网国产| av播播在线观看一区| 久久久精品94久久精品| 日韩亚洲欧美综合| 乱人伦中国视频| 国产精品.久久久| 日本av免费视频播放| av天堂久久9| 麻豆乱淫一区二区| 久热这里只有精品99| 狠狠婷婷综合久久久久久88av| 国产极品粉嫩免费观看在线 | 国产在线一区二区三区精| 人妻制服诱惑在线中文字幕| 制服诱惑二区| 日本av手机在线免费观看| 国产乱来视频区| 日韩在线高清观看一区二区三区| 亚洲精品自拍成人| 国产亚洲最大av| 亚洲精品成人av观看孕妇| h视频一区二区三区| 三级国产精品欧美在线观看| 99久久中文字幕三级久久日本| 欧美日韩国产mv在线观看视频| 老熟女久久久| 黑丝袜美女国产一区| 观看av在线不卡| 久久人人爽av亚洲精品天堂| 99久久精品一区二区三区| 蜜桃国产av成人99| 亚洲人成网站在线观看播放| 22中文网久久字幕| 亚洲精品色激情综合| 性色avwww在线观看| 国产在线一区二区三区精| 免费高清在线观看日韩| 69精品国产乱码久久久| 国产黄片视频在线免费观看| 婷婷成人精品国产| 日韩,欧美,国产一区二区三区| 亚洲熟女精品中文字幕| 校园人妻丝袜中文字幕| 18禁动态无遮挡网站| 国产成人精品福利久久| 精品一区二区免费观看| 精品国产国语对白av| 精品亚洲成国产av| 国产熟女欧美一区二区| 国产一区亚洲一区在线观看| 免费观看的影片在线观看| 男女无遮挡免费网站观看| 免费观看在线日韩| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品一区www在线观看| 国产午夜精品一二区理论片| 久久99精品国语久久久| 在线观看美女被高潮喷水网站| 妹子高潮喷水视频| 日韩,欧美,国产一区二区三区| 视频区图区小说| 午夜91福利影院| 国产日韩欧美在线精品| 日韩人妻高清精品专区| 91久久精品国产一区二区三区| 日本黄色片子视频| 日本黄色日本黄色录像| 最新中文字幕久久久久| 丰满少妇做爰视频| 亚洲av成人精品一区久久| 久久久亚洲精品成人影院| 国产视频首页在线观看| 天堂俺去俺来也www色官网| www.av在线官网国产| 女人精品久久久久毛片| 亚洲国产色片| 熟女av电影| 欧美人与性动交α欧美精品济南到 | 一本一本综合久久| 赤兔流量卡办理| 国产一区二区三区av在线| 青青草视频在线视频观看| 中文字幕制服av| 美女脱内裤让男人舔精品视频| 2018国产大陆天天弄谢| 天堂俺去俺来也www色官网| 国产精品一区二区三区四区免费观看| 成人国产麻豆网| 2021少妇久久久久久久久久久| 久久久精品免费免费高清| 日本av免费视频播放| 亚洲婷婷狠狠爱综合网| tube8黄色片| 婷婷色麻豆天堂久久| 欧美国产精品一级二级三级| 国产精品久久久久久精品电影小说| 国产精品人妻久久久久久| 黄片无遮挡物在线观看| 国精品久久久久久国模美| 国产一级毛片在线| 亚洲精品乱码久久久v下载方式| 亚洲欧美精品自产自拍| 日本猛色少妇xxxxx猛交久久| 99热这里只有是精品在线观看| 99re6热这里在线精品视频| 国产午夜精品久久久久久一区二区三区| 99久久综合免费| 日韩,欧美,国产一区二区三区| 91久久精品国产一区二区三区| 欧美97在线视频| 女人精品久久久久毛片| 涩涩av久久男人的天堂| 色吧在线观看| tube8黄色片| 99热网站在线观看| 久久精品国产鲁丝片午夜精品| 看非洲黑人一级黄片| 国产视频内射| 欧美日韩成人在线一区二区| 国产精品一二三区在线看| 母亲3免费完整高清在线观看 | 少妇人妻久久综合中文| 十八禁高潮呻吟视频| 少妇被粗大猛烈的视频| 99热这里只有精品一区| 国产欧美日韩综合在线一区二区| 特大巨黑吊av在线直播| 美女主播在线视频| 18禁动态无遮挡网站| 亚洲天堂av无毛| 妹子高潮喷水视频| 日本色播在线视频| 国产一级毛片在线| av福利片在线| 久久久久久人妻| 一区在线观看完整版| 麻豆精品久久久久久蜜桃| 久久久久久久久久久久大奶| 最近的中文字幕免费完整| 伊人亚洲综合成人网| 在线 av 中文字幕| 22中文网久久字幕| 在线免费观看不下载黄p国产| 久久免费观看电影| 韩国av在线不卡| 亚洲av中文av极速乱| 黄色视频在线播放观看不卡| 精品人妻偷拍中文字幕| 欧美xxxx性猛交bbbb| 精品人妻偷拍中文字幕| 一边摸一边做爽爽视频免费| 街头女战士在线观看网站| 亚洲精品成人av观看孕妇| 午夜日本视频在线| 国产69精品久久久久777片| 日产精品乱码卡一卡2卡三| 国产精品99久久99久久久不卡 | 狠狠精品人妻久久久久久综合| 少妇被粗大猛烈的视频| 老熟女久久久| 中国三级夫妇交换| 亚洲第一av免费看| 国产国拍精品亚洲av在线观看| 日本欧美国产在线视频| 最近中文字幕高清免费大全6| 亚洲精品av麻豆狂野| 亚洲少妇的诱惑av| 久久久午夜欧美精品| 最后的刺客免费高清国语| 国产乱来视频区| 国产免费一级a男人的天堂| 人人妻人人添人人爽欧美一区卜| 国产一区有黄有色的免费视频| 亚洲激情五月婷婷啪啪| 亚洲av电影在线观看一区二区三区| 国产一区二区在线观看日韩| 91精品伊人久久大香线蕉| 亚洲av成人精品一区久久| 免费观看a级毛片全部| 亚洲色图综合在线观看| 一级爰片在线观看| 欧美日韩视频精品一区| 老司机影院成人| 亚洲中文av在线| 中文精品一卡2卡3卡4更新| 亚洲精品456在线播放app| 欧美精品人与动牲交sv欧美| 成人无遮挡网站| 熟女av电影| 2018国产大陆天天弄谢| 亚洲精品乱码久久久v下载方式| 国产极品天堂在线| 美女内射精品一级片tv| 国产又色又爽无遮挡免| 一级毛片 在线播放| 精品99又大又爽又粗少妇毛片| 国产精品偷伦视频观看了| a 毛片基地| av.在线天堂| av国产久精品久网站免费入址| 精品人妻在线不人妻| 亚洲av成人精品一二三区| 女性被躁到高潮视频| 大陆偷拍与自拍| 日日摸夜夜添夜夜爱| 久热这里只有精品99| 精品久久蜜臀av无| 一边摸一边做爽爽视频免费| a级毛片在线看网站| 成人无遮挡网站| 精品酒店卫生间| 久久青草综合色| 在线观看人妻少妇| 日本色播在线视频| 男人爽女人下面视频在线观看| 日本黄色日本黄色录像| 日韩强制内射视频| 日日啪夜夜爽| 99热国产这里只有精品6| 99九九在线精品视频| 久久久欧美国产精品| 国产国拍精品亚洲av在线观看| 中文字幕人妻熟人妻熟丝袜美| av播播在线观看一区| 久久久久国产精品人妻一区二区| 国产高清不卡午夜福利| 日韩制服骚丝袜av| 亚洲国产成人一精品久久久| 国产精品国产三级专区第一集| 99久国产av精品国产电影| 十八禁网站网址无遮挡| av一本久久久久| 老熟女久久久| 欧美日韩精品成人综合77777| 妹子高潮喷水视频| 日本与韩国留学比较| 国产精品一二三区在线看| 久久av网站| 熟女电影av网| 久久午夜综合久久蜜桃| 国产精品不卡视频一区二区| xxxhd国产人妻xxx| 国产av码专区亚洲av| 简卡轻食公司| 制服诱惑二区| 国产亚洲av片在线观看秒播厂| av又黄又爽大尺度在线免费看| 亚洲丝袜综合中文字幕| 一级二级三级毛片免费看| 少妇高潮的动态图| 老司机亚洲免费影院| 一级,二级,三级黄色视频| 一二三四中文在线观看免费高清| 国产免费现黄频在线看| 亚洲av免费高清在线观看| 亚洲激情五月婷婷啪啪| 国产黄色视频一区二区在线观看| 高清黄色对白视频在线免费看| 国产精品99久久99久久久不卡 | 国产男女内射视频| 久久久久久久久久成人| 亚洲高清免费不卡视频| 欧美精品一区二区大全| 热re99久久精品国产66热6| 搡女人真爽免费视频火全软件| 免费少妇av软件| 制服人妻中文乱码| 国产精品一区二区在线观看99| 亚洲四区av| 亚洲精品久久久久久婷婷小说| 午夜精品国产一区二区电影| 啦啦啦在线观看免费高清www| 欧美激情 高清一区二区三区| 91精品伊人久久大香线蕉| 大话2 男鬼变身卡| 欧美3d第一页| 中国三级夫妇交换| 亚洲精品一二三| 国产成人a∨麻豆精品| 久久99一区二区三区| 韩国av在线不卡| 免费播放大片免费观看视频在线观看| 欧美97在线视频| 青青草视频在线视频观看| 黑丝袜美女国产一区| 国产成人aa在线观看| av福利片在线| 男女无遮挡免费网站观看| 视频中文字幕在线观看| 免费观看av网站的网址| 麻豆乱淫一区二区| 国产成人91sexporn| 日韩欧美一区视频在线观看| 啦啦啦视频在线资源免费观看| 国产免费福利视频在线观看| 韩国高清视频一区二区三区| 女人精品久久久久毛片| 日韩在线高清观看一区二区三区| videossex国产| 91午夜精品亚洲一区二区三区| 中文字幕免费在线视频6| 亚洲综合色网址| 国产毛片在线视频| 成人无遮挡网站| 亚洲精品乱码久久久久久按摩| 中文字幕免费在线视频6| 久久久国产一区二区| 久久精品国产亚洲av涩爱| 一区二区三区乱码不卡18| 18在线观看网站| 蜜桃久久精品国产亚洲av| 亚洲内射少妇av| 免费大片黄手机在线观看| 老司机亚洲免费影院| 欧美日韩视频精品一区| 欧美一级a爱片免费观看看| 午夜免费鲁丝| 天堂中文最新版在线下载| 天堂俺去俺来也www色官网| 国产成人免费观看mmmm| 国产伦理片在线播放av一区| 国模一区二区三区四区视频| 国产极品天堂在线| 新久久久久国产一级毛片| av天堂久久9| 最近中文字幕2019免费版| 97精品久久久久久久久久精品| 精品国产一区二区久久| 另类亚洲欧美激情| 国产免费一区二区三区四区乱码| 最新的欧美精品一区二区| 精品人妻一区二区三区麻豆| 亚洲av电影在线观看一区二区三区| 99国产综合亚洲精品| 亚洲精品第二区| 国产成人91sexporn| 97超碰精品成人国产| 午夜老司机福利剧场| 久久99精品国语久久久| av在线老鸭窝| 日本爱情动作片www.在线观看| 久久99一区二区三区| 久久久久久久大尺度免费视频| 九色成人免费人妻av| 国产白丝娇喘喷水9色精品| 狂野欧美激情性bbbbbb| 精品久久久久久电影网| 中国美白少妇内射xxxbb| 成人国产av品久久久| 高清av免费在线| 大香蕉久久网| 久久久久国产精品人妻一区二区| 99热网站在线观看| 麻豆成人av视频| 最后的刺客免费高清国语| av播播在线观看一区| 一级爰片在线观看| 亚洲,一卡二卡三卡| 美女大奶头黄色视频| 乱人伦中国视频| 成人漫画全彩无遮挡| 欧美xxⅹ黑人| 亚洲精品国产av蜜桃| a级毛片免费高清观看在线播放| 成人亚洲欧美一区二区av| a级毛片黄视频| 只有这里有精品99| 成年av动漫网址| 91久久精品国产一区二区三区| 九九在线视频观看精品| 国产 精品1| 国产精品无大码| 国产在线一区二区三区精| 成人18禁高潮啪啪吃奶动态图 | 国产精品久久久久久av不卡| 久久久午夜欧美精品| 国产精品麻豆人妻色哟哟久久| 亚洲色图综合在线观看| 91午夜精品亚洲一区二区三区| 精品亚洲成a人片在线观看| 午夜免费鲁丝| 久久热精品热| 乱码一卡2卡4卡精品| 久久青草综合色| a级毛片在线看网站| 在线观看人妻少妇| 国产高清有码在线观看视频| 亚洲精品自拍成人| 国产高清三级在线| 黑人欧美特级aaaaaa片| 两个人的视频大全免费| 乱人伦中国视频| 青青草视频在线视频观看| 中文天堂在线官网| 亚洲人成网站在线播| 国产一区亚洲一区在线观看| 少妇被粗大的猛进出69影院 | 精品人妻熟女av久视频| 亚洲精品视频女| 午夜精品国产一区二区电影| 各种免费的搞黄视频| 在线播放无遮挡| 精品卡一卡二卡四卡免费| 91久久精品电影网| 三上悠亚av全集在线观看| 又粗又硬又长又爽又黄的视频| 午夜免费男女啪啪视频观看| 校园人妻丝袜中文字幕| 亚洲成色77777|