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

    梧州茶樹種質(zhì)資源的遺傳多樣性及親緣關(guān)系分析

    2022-10-31 02:16:10王留彬黃麗蘊滕翠琴吳立赟成浩于翠平王麗鴛
    茶葉科學(xué) 2022年5期
    關(guān)鍵詞:南渡六堡茶樹

    王留彬,黃麗蘊,滕翠琴,吳立赟,成浩,于翠平*,王麗鴛*

    梧州茶樹種質(zhì)資源的遺傳多樣性及親緣關(guān)系分析

    王留彬1,黃麗蘊2,滕翠琴2,吳立赟1,成浩1,于翠平2*,王麗鴛1*

    1. 中國農(nóng)業(yè)科學(xué)院茶葉研究所/國家茶樹改良中心,浙江 杭州 310008;2. 梧州市農(nóng)業(yè)科學(xué)研究所,廣西 梧州 543003

    基于SSR標(biāo)記對梧州六堡鎮(zhèn)群體種和南渡鎮(zhèn)野生大茶樹種質(zhì)資源的遺傳多樣性和親緣關(guān)系進行了分析,并篩選出用于該種質(zhì)資源鑒別的核心分子標(biāo)記。研究結(jié)果如下:(1)17對SSR引物在供試材料中共擴增得到98個等位基因,每對SSR引物擴增的等位位點為3~8個,平均每個位點5.764?7個等位基因;(2)從17個SSR分子標(biāo)記中篩選出8個核心標(biāo)記組合即可區(qū)分每份種質(zhì)資源;(3)六堡鎮(zhèn)茶樹種質(zhì)資源的平均等位基因數(shù)(A)、基因型數(shù)、基因多樣性(H)、多態(tài)性信息含量(PIC)分別為4.647?1、7.000?0、0.675?4、0.628?3,高于南渡鎮(zhèn)野生大茶樹種質(zhì)資源,與栽培種茶樹群體接近;(4)聚類分析表明,六堡鎮(zhèn)茶樹群體部分種質(zhì)資源單獨聚為一類,部分與云南的大葉種茶樹,少量與浙江、貴州地方栽培種聚為一類;而南渡鎮(zhèn)野生茶樹種質(zhì)資源單獨聚為一類,僅有2個六堡鎮(zhèn)的種質(zhì)材料散落其間。綜上所述,梧州茶樹種質(zhì)資源豐富,遺傳多樣性較高,研究結(jié)果為進一步開發(fā)和利用該資源奠定了基礎(chǔ)。

    SSR標(biāo)記;遺傳多樣性;梧州市;茶樹

    廣西地處華南地區(qū),與云南、貴州相鄰地區(qū)被認為是茶樹的次生起源中心,憑借優(yōu)越的生態(tài)和地理條件,蘊含十分豐富的茶樹資源[1]。近年來,廣西茶園面積和茶葉產(chǎn)量增長迅速,2019年,茶園面積為77.3?khm2,茶葉總產(chǎn)量為8.28萬t[2],2020年,茶園面積達到91.3?khm2,茶葉總產(chǎn)量9.13萬t[3]。課題組前期調(diào)查發(fā)現(xiàn)廣西梧州市具有較多種類的群體種茶樹資源,伴隨產(chǎn)業(yè)的崛起,六堡茶群體種的重要性越發(fā)凸顯。邱瑞瑾等[4]和覃秀菊等[5]分別研究了梧州市六堡茶群體種資源的分布和葉片顯微結(jié)構(gòu)。滕翠琴等[6]在梧州市南渡鎮(zhèn)發(fā)現(xiàn)了集中成片的野生大茶樹群落,并對其表型特征、生長環(huán)境等進行了調(diào)查。這些獨特的種質(zhì)資源在茶葉生產(chǎn)應(yīng)用、品種創(chuàng)新以及維持物種多樣性等方面具有重要的意義。目前,使用DNA分子標(biāo)記對上述種質(zhì)資源的研究尚少。

    SSR分子標(biāo)記作為群體遺傳多樣性和親緣關(guān)系分析中的有效工具,已在茶樹遺傳學(xué)分析中廣泛使用。如:Tan等[7]使用30對SSR標(biāo)記分析了128份中國主栽無性系茶樹品種間的遺傳關(guān)系,證實了40余對茶樹品種的親子關(guān)系。周斌等[8]利用15對SSR引物對50份雷波縣野生茶樹資源進行了遺傳多樣性和親緣關(guān)系分析,推測雷波地區(qū)附近可能是四川野生茶樹的發(fā)源中心。余書平等[9]使用14對SSR標(biāo)記對浙江開化縣不同種質(zhì)資源進行了分析,發(fā)現(xiàn)開化縣茶樹種質(zhì)資源遺傳多樣性豐富。本研究通過17對SSR引物,分析了梧州六堡鎮(zhèn)和南渡鎮(zhèn)兩個地區(qū)茶樹種質(zhì)資源遺傳多樣性,篩選出鑒定該種質(zhì)資源的核心SSR標(biāo)記。比較了兩地區(qū)茶樹種質(zhì)資源的遺傳信息和親緣關(guān)系,并與多個栽培種進行比較,旨在解析梧州茶樹種質(zhì)資源的分子遺傳背景,為將來梧州市優(yōu)異資源的收集、保護和科學(xué)利用提供研究基礎(chǔ)。

    1 材料與方法

    1.1 材料

    19個無性系栽培(ZP)品種(表1)種植于中國農(nóng)業(yè)科學(xué)院茶葉研究所嵊州基地,分別編號1~19,作為對照品種。以梧州市六堡鎮(zhèn)(LB)大中村(111°22′E,23°48′N)的13個六堡茶群體種資源和梧州市南渡鎮(zhèn)(ND)天龍頂(110°51′E,22°42′N,海拔600~700?m)采集的12個野生大茶樹種質(zhì)資源為試驗材料,梧州市六堡鎮(zhèn)大中村13個六堡茶群體種資源編號為20~32,葉片多樣性較高,大、小葉類均有,以中葉類為主。梧州市南渡鎮(zhèn)天龍頂?shù)?2個野生大茶樹種質(zhì)資源編號為33~44,樹型為喬木型。所用材料詳細信息及編號見表1。

    1.2 DNA提取及檢測

    DNA的提取使用植物基因組提取試劑盒[DP350,天根生化科技(北京)有限公司],并使用超微量分光光度計(Nano Drop ND-1000,賽默飛世爾科技公司)檢測DNA濃度和純度,使用1%瓊脂糖凝膠電泳檢測DNA的完整性,對DNA質(zhì)量不高或條帶不清楚的DNA樣品重新提取。將每份樣品DNA稀釋到30?ng·μL-1,放置于–20℃冰箱保存?zhèn)溆谩?/p>

    1.3 候選SSR引物的PCR擴增

    從課題組已設(shè)計的SSR引物中[10]挑選多態(tài)性高、條帶清晰的17對SSR引物(引物序列均已上傳至NCBI)用于本試驗,具體引物信息見表2。PCR反應(yīng)體系:Reaction Mix 5.0?μL(含DNA Polymerase、2×PCR Buffer、MgCl2和dNTP),模板DNA 2.0?μL,Primer-F、Primer-R(10?ng·μL-1)各0.2?μL,加ddH2O至10?μL。PCR循環(huán)參數(shù):94℃預(yù)變性4?min;94℃變性30?s,52~58℃退火30 s,72℃延伸30?s,共計35個循環(huán);72℃延伸10?min;4℃保存。

    表1 供試材料信息

    1.4 擴增產(chǎn)物的PAGE電泳和銀染顯色

    PCR擴增使用垂直電泳儀(CBS MGV-202-33,電壓120?V,時間90?min)進行10%的聚丙烯酰胺凝膠電泳(Polyacrylamide gel electrophoresis,PAGE),具體配制方法參照文獻[11]。電泳結(jié)束后進行銀染顯色,并用保鮮膜將膠片保存,拍照記錄。

    1.5 數(shù)據(jù)處理

    將電泳圖中位于目的片段范圍內(nèi)的清晰條帶,依據(jù)分子量大小依次記作A、B、C、D……,因條帶大小不同,基因型分別記錄為AA、AB、BB、AC、AD、CD等(單一條帶為純合基因型,兩個條帶為雜合基因型),條帶缺失記為“-”,計入Excel表,建立起原始基因型數(shù)據(jù)矩陣。利用GenAlEx 6.5.0.3軟件[12]進行基因頻率、雜合度、標(biāo)記鑒別力等一系列遺傳參數(shù)的運算。具體方法參照文獻[9]處理方法。根據(jù)樣本的基因型數(shù)據(jù),將其轉(zhuǎn)換成PowerMarker v3.25軟件[13]要求的數(shù)據(jù)格式?;凇癗ei.1983”計算樣本間的遺傳距離并進行UPGMA進化樹構(gòu)建。使用FigTreev1.2.4軟件繪制UPGMA進化樹。

    表2 17對SSR引物信息

    2 結(jié)果與分析

    2.1 44份茶樹種質(zhì)SSR標(biāo)記多態(tài)性分析

    本研究使用的17對引物呈現(xiàn)出較高水平的多態(tài)性(表3)。除3對引物(CsFM1051、CsFM1839、CsFM1562)擴增目的片段的成功率較低外,其余引物幾乎均能在所有樣品中成功擴增出目的片段??傮w來看,17對SSR引物共獲得98條等位基因,每對SSR引物擴增出的等位位點為3~8個,平均每個位點檢測到5.764?7個等位基因,平均有效等位位點數(shù)為3.772?2個。在這些SSR分子標(biāo)記中,CsFM1076擴增出的等位基因最少,僅3個。引物CsFM1131、CsFM1801、CsFM1585、CsFM1651、CsFM1715擴增的等位基因為7~8個,這5對引物的有效位點(Ne)均超過4個,說明它們在樣品中擴增多態(tài)性較好,具有較強的鑒別能力。17對引物擴增位點的多態(tài)性信息含量(PIC)為0.315?7~0.803?9,平均值為0.655?1,也說明本研究選用的SSR引物多態(tài)性較好。

    2.2 標(biāo)記的鑒別能力

    對基因分型數(shù)據(jù)進行復(fù)合位點分析,非同一概率值(PE)使用GenAlEx 6.5.0.3軟件分析得到。結(jié)果如表4所示,17個SSR位點的PE-1值范圍在0.172?8~0.654?7,平均值為0.478?2;PE-2值的范圍在0.062?7~0.481?5,平均值為0.312?2;PE-3的范圍為0.285?9~0.830?6,平均值為0.651?3。如圖1所示,當(dāng)17個SSR位點組合后,2個隨機樣本基因型組合不同的PE-1為1.000?0;若缺失1個位點和去除1個純合位點時,2個隨機樣本基因型組合不同的概率PE-2和PE-3分別為0.998?7和1.000?0,該結(jié)果表明本研究選用的標(biāo)記多態(tài)性較高,用17個SSR標(biāo)記完全可以鑒定出所選用的每一份種質(zhì)材料。

    表3 17對SSR標(biāo)記的遺傳參數(shù)

    研究結(jié)果顯示,選用的SSR標(biāo)記增多(圖1),其PE值顯著增加,材料的鑒定準確性亦顯著提高。數(shù)據(jù)顯示當(dāng)選用8個核心標(biāo)記時,PE-1、PE-2、PE-3分別超過了0.99、0.96、0.99,即可成功鑒別出每一份材料,因此將PE值排名最高的8個標(biāo)記作為鑒定該種質(zhì)資源的核心標(biāo)記,分別是CsFM1828、CsFM1715、CsFM1131、CsFM1068、CsFM1599、CsFM1801、CsFM1585、CsFM1651。

    2.3 2個不同地區(qū)茶樹群體的遺傳多樣性

    梧州的六堡鎮(zhèn)大中村和南渡鎮(zhèn)云龍頂相距約133?km。為探明兩地茶樹種質(zhì)資源遺傳多樣性與遺傳距離的差異,本研究將上述二地取樣采集的茶樹資源分別作為單獨群體,并以19個栽培品種作為對照,進行遺傳多樣性分析。結(jié)果顯示(表5),六堡鎮(zhèn)茶樹群體平均等位基因數(shù)為4.647?1個,基因型數(shù)為7.000?0,基因型多樣性為0.675?4,雜合度為0.664?7,多態(tài)性信息含量為0.628?3。六堡鎮(zhèn)茶樹群體的平均等位基因、基因型數(shù)、基因多樣性和多態(tài)性信息含量均接近于栽培種。南渡鎮(zhèn)野生茶樹群體的平均等位基因數(shù)(4.294?1)、基因型數(shù)(5.823?5)、基因多樣性(0.600?9)、多態(tài)性信息含量(0.558?0)均小于六堡鎮(zhèn)茶樹群體,僅雜合度(0.664?9)略高于六堡鎮(zhèn)茶樹群體。

    基于群體之間的遺傳距離(表6),利用UPGMA聚類分析發(fā)現(xiàn),六堡鎮(zhèn)茶樹群體與栽培種遺傳距離最近為0.129?3,被聚為一類;南渡鎮(zhèn)的野生茶樹群體與他們的遺傳距離較遠(圖2)。由此表明六堡鎮(zhèn)茶樹群體與栽培種親緣關(guān)系較近,南渡鎮(zhèn)茶樹群體與六堡鎮(zhèn)茶樹群體、栽培種均有一定的遺傳距離,這與南渡鎮(zhèn)的茶樹群體是野生大茶樹的研究結(jié)果相一致。

    表4 17個SSR數(shù)據(jù)的PE值

    注:PE-1表示位點未缺失時的PE值,PE-2表示缺失1個位點時的PE值,PE-3表示去除1個純合位點的PE值。下同

    Note: PE-1 represents PE value when one of locus is missing, PE-2 represents PE value when one locus is missing, PE-3 represents PE value when one of homozygous locus is removed. The same below

    圖1 不同標(biāo)記組合數(shù)的PE值

    表5 3個茶樹群體的遺傳多樣性分析

    注:LB,六堡鎮(zhèn)茶樹群體;ND,南渡鎮(zhèn)茶樹群體;ZP,栽培種。下同

    Note: LB, Tea plant populations in Liubao town. ND, Tea plant populations in Nandu town. ZP, Cultivar. The same below

    表6 3個茶樹群體的遺傳距離

    圖2 3個茶樹群體UPGMA聚類圖

    2.4 茶樹種質(zhì)資源個體間的聚類分析

    為進一步了解梧州市茶樹種質(zhì)資源的遺傳學(xué)地位,根據(jù)17對SSR引物擴增數(shù)據(jù)統(tǒng)計結(jié)果,基于Nei's遺傳距離,采用UPGMA法聚類得到茶樹聚類樹狀圖,將茶樹資源大致分為5個類群(圖3)。其中類群1主要是六堡鎮(zhèn)的LP2、LP3、LP4、LP5,說明它們具有較為相近的遺傳背景。類群2主要為茶樹栽培種,迎霜是福鼎大白茶和云南大葉種的雜交子代,它們具有親緣關(guān)系[14],在本研究中聚在一類。類群3由部分六堡鎮(zhèn)的種質(zhì)資源和幾個栽培種聚類而成,如LP9、LP1、LP8和大葉種的紫娟、英紅1號、英紅9號。類群4中,LP11、LP12和云抗10號,LP6、LP7、黔眉419和中黃1號分別聚在一類。類群5主要是來自南渡鎮(zhèn)的野生大茶樹,僅有2個六堡鎮(zhèn)的LP10、LP13和南渡鎮(zhèn)群體聚在一類。

    3 討論

    3.1 梧州茶樹種質(zhì)資源的遺傳多樣性

    茶樹種質(zhì)資源及遺傳育種研究以遺傳多樣性為基礎(chǔ),茶樹的遺傳多樣性研究已從最早的形態(tài)學(xué)標(biāo)記轉(zhuǎn)變?yōu)榉肿訕?biāo)記[15]。SSR分子標(biāo)記技術(shù)已廣泛使用在茶樹遺傳多樣性及種質(zhì)鑒定研究等方面[9]。徐禮羿等[16]使用8個核心SSR標(biāo)記即可成功識別200份以上的茶樹種質(zhì)資源。本研究使用的17對標(biāo)記即部分來源于上述SSR標(biāo)記,部分來源于Tan等[7]開發(fā)的SSR標(biāo)記,選取的17對引物均具有多態(tài)性。基于復(fù)合位點的非同一性分析發(fā)現(xiàn),選用8個核心SSR標(biāo)記組合時,便可實現(xiàn)對本研究的44份材料的準確識別。多態(tài)性信息量(PIC)為評價引物多態(tài)性的重要指標(biāo),當(dāng)PIC>0.5時,該位點為高度多態(tài)位點;當(dāng)0.25<PIC≤0.5時,該位點為中度多態(tài)位點;當(dāng)PIC≤0.25時,該位點為低度多態(tài)位點[17];本研究使用的引物平均PIC為0.655?1,PIC>0.5的有14對,占總引物的82.35%;其余3對引物為0.25<PIC≤0.5,說明篩選的SSR標(biāo)記具有較高的多態(tài)性,可較好地應(yīng)用于茶樹遺傳圖譜的構(gòu)建、遺傳多樣性和親緣關(guān)系分析等研究。

    本研究以梧州市六堡鎮(zhèn)群體種和南渡鎮(zhèn)野生大茶樹的種質(zhì)為研究材料,通過DNA分子標(biāo)記和聚類分析較為全面地研究了兩地茶樹資源的遺傳多樣性水平和親緣關(guān)系。目前針對上述種質(zhì)資源的研究較少,僅有周炎花等[18]選取了數(shù)個六堡鎮(zhèn)茶樹品種進行了研究,發(fā)現(xiàn)其遺傳多樣性較高。本研究發(fā)現(xiàn)六堡鎮(zhèn)茶樹群體種遺傳多樣性豐富(H:0.675?4,PIC:0.628?3),甚至超過栽培品種(H:0.665?0,PIC:0.617?8),這與本研究調(diào)查的六堡鎮(zhèn)茶樹群體表型特征多樣性相符合;在聚類分析中,4個六堡鎮(zhèn)群體種(LP2、LP3、LP4、LP5)單獨聚為類群1,推測它們具有獨特的遺傳背景。以上結(jié)果表明,六堡鎮(zhèn)群體種遺傳多樣性豐富,具有優(yōu)異資源選育的潛力。

    圖3 茶樹種質(zhì)資源聚類分析

    滕翠琴等[6]通過對野生大茶樹的資源調(diào)查,發(fā)現(xiàn)南渡鎮(zhèn)的野生大茶樹是梧州市內(nèi)面積最大的野生茶樹資源。我們將從該地收集到的12個種質(zhì)資源單獨分析,發(fā)現(xiàn)其遺傳多樣性水平(H:0.600?9,PIC:0.558?0)接近云南茶樹資源的遺傳多樣性水平(H:0.609,PIC:0.578),高于四川茶樹資源的遺傳多樣性水平(H:0.525,PIC:0.492)[19],且雜合度(He:0.664?9)高于云南西雙版納地區(qū)野生茶樹的雜合度(He:0.353)[20]。據(jù)此推斷南渡鎮(zhèn)野生大茶樹的遺傳多樣性較高,具有非常重要的利用價值,為深入研究該地區(qū)野生大茶樹提供了良好的種質(zhì)資源。該地區(qū)野生大茶樹的起源、分類是有待今后進一步探究的科學(xué)問題。

    3.2 梧州市茶樹種質(zhì)資源的親緣關(guān)系分析

    基于Nei′s遺傳距離,本研究中的44份供試茶樹種質(zhì)資源被分為5個類群,其中六堡鎮(zhèn)茶樹群體在5個類群中均有分布。如在類群3中,LP9、LP1、LP8、紫娟、英紅1號和英紅9號聚在了一起。在類群4中LP11、LP12和云抗10號聚在一起。紫娟和云抗10號均是大葉種茶樹品種[14],英紅1號和英紅9號是由廣東省農(nóng)業(yè)科學(xué)院從云南省大葉種種質(zhì)資源()中選育出來的茶樹品種,同屬于大葉種茶樹品種[21]。在形態(tài)調(diào)查中也發(fā)現(xiàn)LP9、LP1、LP8、LP11、LP12具有大葉種茶樹品種的表型特征,表明它們具有一定的親緣關(guān)系。在其他類群中,六堡鎮(zhèn)茶樹群體與浙江、貴州等地的中小葉栽培種茶樹具有較近的親緣關(guān)系。

    聚類分析結(jié)果表明,南渡鎮(zhèn)天龍頂野生古茶樹群體單獨聚在一起,并與六堡鎮(zhèn)的茶樹種質(zhì)資源和栽培種茶樹具有一定的遺傳距離,這與前人研究結(jié)果一致[22-23]。在個體聚類分析中,僅有來自六堡鎮(zhèn)的LP10、LP13和南渡鎮(zhèn)的天龍頂野生古茶樹聚在了一起。這些研究結(jié)果表明野生種與栽培茶樹品種在親緣關(guān)系上存在較大的差異性。周斌等[8]同樣發(fā)現(xiàn)雷波野生茶樹與栽培種在親緣關(guān)系上存在較大差異。前人的研究表明,基因交流是導(dǎo)致居群間茶樹種質(zhì)資源遺傳分化的主要原因,長期的基因交流,會導(dǎo)致茶樹資源分化程度較低,遺傳多樣性下降[24]。從本研究的結(jié)果看,南渡鎮(zhèn)野生大茶樹與六堡鎮(zhèn)茶樹群體、栽培種親緣關(guān)系都差異較大,說明群體之間的基因交流較少,仍存在一定程度的分化。南渡鎮(zhèn)野生大茶樹群落周圍樹林密布,環(huán)境復(fù)雜,而茶樹主要依靠昆蟲授粉[23],推測復(fù)雜的環(huán)境和地理距離使該野生大茶樹群落內(nèi)長期自交。這有可能是該群體遺傳多樣性低于六堡鎮(zhèn)茶樹群體和栽培種遺傳多樣性的原因。

    3.3 梧州市茶樹種質(zhì)資源的保護

    本研究表明,梧州市茶樹種質(zhì)資源豐富,遺傳多樣性高,具有非常重要的利用價值。如六堡鎮(zhèn)的部分茶樹種質(zhì)和南渡鎮(zhèn)野生茶樹資源與栽培種茶樹親緣關(guān)系較遠,使用該群體的材料與栽培種茶樹進行雜交育種可能會產(chǎn)生較好的雜種優(yōu)勢[25],這在茶樹品種多樣性方面具有積極的作用。此外由于市場需求,野生茶價格日益增高,導(dǎo)致對野生茶樹過度和破壞性采摘,從而影響野生茶樹資源的遺傳多樣性和種群規(guī)模。因此,需要進一步對梧州市內(nèi)的茶樹種質(zhì)資源進行保護,如對該地區(qū)的野生優(yōu)質(zhì)資源進行扦插擴繁種植,建立資源圃持久保存等。

    [1] 王新超, 陳亮, 楊亞軍. 廣西茶樹資源生化成分多樣性分析[J]. 植物遺傳資源學(xué)報, 2010, 11(3): 309-314.

    Wang X C, Chen L, Yang Y J. Biochenical diversity analysis of tea germplasms in Guangxi [J]. Journal of Plant Genetic Resources, 2010, 11(3): 309-314.

    [2] 韋靜峰, 劉曉東. 廣西茶產(chǎn)業(yè)發(fā)展歷程與展望[J]. 廣西農(nóng)學(xué)報, 2019, 34(4): 65-73.

    Wei J F, Liu X D. The development course and prospect of Guangxi tea industry [J]. Journal of Guangxi Agriculture, 2019, 34(4): 65-73.

    [3] 寧夏, 孔令政, 黃艷芳, 等. 鄉(xiāng)村振興背景下廣西茶葉產(chǎn)業(yè)發(fā)展研究[J]. 熱帶農(nóng)業(yè)科學(xué), 2021, 41(8): 125-132.

    Ning X, Kong L Z, Huang Y F, et al. Development of tea industry in Guangxi under the background of rural revitalization [J]. Chinese Journal of Tropical Agriculture, 2021, 41(8): 125-132.

    [4] 邱瑞瑾, 龍志榮, 于翠平, 等. 六堡茶群體種資源的調(diào)查、開發(fā)和保護[J]. 中國茶葉加工, 2020(3): 15-19, 26.

    Qiu R J, Long Z R, Yu C P, et al. Investigation, development and protection of population resources of Liupao tea [J]. China Tea Processing, 2020(3): 15-19, 26.

    [5] 覃秀菊, 李鳳英, 何建棟, 等. 廣西茶樹新品種品系葉片解剖結(jié)構(gòu)特征與特性關(guān)系的研究[J]. 中國農(nóng)學(xué)通報, 2009, 25(10): 36-39.

    Qin X J, Li F Y, He J D, et al. The relations between the blade anatomical structure and characters of the new tea varieties in Guangxi [J]. Chinese Agricultural Science Bulletin, 2009, 25(10): 36-39.

    [6] 滕翠琴, 邱瑞瑾, 于翠平, 等. 岑溪市野生古茶樹資源特征與開發(fā)利用現(xiàn)狀調(diào)查[J]. 南方農(nóng)業(yè), 2020, 14(4): 45-48.

    Teng C Q, Qiu R J, Yu C P, et al. Investigation on the characteristics and exploitation of wild ancient tea tree resources in Cenxi city [J]. South China Agriculture, 2020, 14(4): 45-48.

    [7] Tan L Q, Peng M, Xu L Y, et al. Fingerprinting 128 Chinese clonal tea cultivars using SSR markers provides new insights into their pedigree relationships [J]. Tree Genetics & Genomes, 2015, 11(5): 1-12.

    [8] 周斌, 王留彬, 徐新濤, 等. 雷波野生茶樹遺傳多樣性及親緣關(guān)系分析[J]. 云南農(nóng)業(yè)大學(xué)學(xué)報(自然科學(xué)), 2020, 35(1): 122-129.

    Zhou B, Wang L B, Xu X T, et al. Analysis of genetic diversity and genetic relationship of Leibo wild tea resources [J]. Journal of Yunnan Agricultural University (Natural Science), 2020, 35(1): 122-129.

    [9] 余書平, 徐禮羿, 吳榮梅, 等. 浙江開化縣茶樹種質(zhì)資源的遺傳多樣性及親緣關(guān)系分析[J]. 茶葉科學(xué), 2020, 40(3): 341-351.

    Yu S P, Xu L Y, Wu R M, et al. Genetic and phylogenetic analysis for resources offrom Kaihua County in Zhejiang Province [J]. Journal of Tea Science, 2020, 40(3): 341-351.

    [10] Tan L Q, Wang L Y, Xu L Y, et al. SSR-based genetic mapping and QTL analysis for timing of spring bud flush, young shoot color, and mature leaf size in tea plant () [J]. Tree Genetics & Genomes, 2016, 12(3): 1-13.

    [11] 黃丹娟. 我國茶樹優(yōu)良品種遺傳多樣性分析及指紋圖譜構(gòu)建[D]. 北京: 中國農(nóng)業(yè)科學(xué)院, 2016.

    Huang D J. Genetic diversity and fingerprinting of Chinese improved tea cultivars [D]. Beijing: Chinese Academy of Agricultural Sciences, 2016.

    [12] Smouse P E, Whitehead M R, Peakall R. An informational diversity framework, illustrated with sexually deceptive orchids in early stages of speciation [J]. Molecular Ecology Resources, 2015, 15(6): 1375-1384.

    [13] Liu K, Muse S V. PowerMarker: an integrated analysis environment for genetic marker analysis [J]. Bioinformatics, 2005, 21(9): 2128-2129.

    [14] 黃海濤, 余繼忠, 鄭旭霞, 等. 茶樹雜交新品系‘杭茶15號’的選育[J]. 茶葉, 2018, 44(3): 121-124.

    Huang H T, Yu J Z, Zheng X X, et al. New tea cultivar ‘Hangcha 15’ breeding by hybrid [J]. Journal of Tea, 2018, 44(3): 121-124.

    [15] 羅亦紓, 張霞, 毛君林, 等. 35份重慶大茶樹種質(zhì)資源遺傳多樣性的SSR分析[J]. 分子植物育種, 2019, 17(19): 6549-6557.

    Luo Y S, Zhang X, Mao J L, et al. Analysis of genetic diversity in 35 ancient tea accessions from chongqing with SSR markers [J]. Molecular Plant Breeding, 2019, 17(19): 6549-6557.

    [16] 徐禮羿, 王麗鴛, 蘇靜靜, 等. 云南雙江勐庫大葉種茶樹基因型和種群結(jié)構(gòu)分析[J]. 植物遺傳資源學(xué)報, 2019, 20(4): 1052-1064.

    Xu L Y, Wang L Y, Su J J, et al. Genotyping and genetic structure analysis for 235 tea plant (var. assamica cv. MengKu-dayecha) in Shuangjiang of Yunnan Province [J]. Journal of Plant Genetic Resources, 2019, 20(4): 1052-1064.

    [17] 李賽君, 雷雨, 段繼華, 等. 基于EST-SSR的祁門種群體遺傳多樣性和親緣關(guān)系分析[J]. 茶葉科學(xué), 2015, 35(4): 329-335.

    Li S J, Lei Y, Duan J H, et al. Analysis of genetic diversity and relationship of Qimen population of tea plants based on EST-SSR markers [J]. Journal of Tea Science, 2015, 35(4): 329-335.

    [18] 周炎花, 喬小燕, 馬春雷, 等. 廣西茶樹地方品種遺傳多樣性和遺傳結(jié)構(gòu)的EST-SSR分析[J]. 林業(yè)科學(xué), 2011, 47(3): 59-67.

    Zhou Y H, Qiao X Y, Ma C L, et al. Genetic diversity and structure of tea landraces from Guangxi based on EST-SSR analysis [J]. Scientia Silvae Sinicae, 2011, 47(3): 59-67.

    [19] Yao M Z, Ma C L, Qiao T T, et al. Diversity distribution and population structure of tea germplasms in China revealed by EST-SSR markers [J]. Tree Genetics & Genomes, 2012, 8(1): 205-220.

    [20] 楊金玲. 西雙版納地區(qū)古茶樹種質(zhì)資源的遺傳多樣性分析[D]. 昆明: 云南大學(xué), 2018.

    Yang J L. Genetic diversity analysis of ancient tea plants in Xishuangbanna area [D]. Kunming: Yunnan University, 2018.

    [21] 盧嘉麗, 王冬梅, 苗愛清, 等. 英紅1號、英紅9號和祁門茶樹芽葉中嘌呤生物堿和茶多酚的HPLC-DAD-MS/MS分析[J]. 中山大學(xué)學(xué)報(自然科學(xué)版), 2009, 48(1): 72-75.

    Lu J L, Wang D M, Miao A Q, et al. HPLC-DAD-MS/MS analysis of purine alkaloids and tea polyphenols in young leaves of Yinghong1, Yinghong 9 and Qimen cultivars [J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 2009, 48(1): 72-75.

    [22] 彭靖茹, 李朝昌, 檀業(yè)維, 等. 基于SSR分子標(biāo)記的廣西德??h和隆林縣野生古茶樹聚類分析[J]. 南方農(nóng)業(yè)學(xué)報, 2019, 50(1): 1-7.

    Peng J R, Li C C, Tan Y W, et al. Clustering analysis for wild ancient tea germplasm resources in Debao County and Longlin County, Guangxi based on SSR molecular markers [J]. Journal of Southern Agriculture, 2019, 50(1): 1-7.

    [23] 黃厚宸, 蔣維昕, 梁彩霞, 等. 廣西六堡茶古茶樹及其子代遺傳多樣性的EST-SSR標(biāo)記分析[J]. 分子植物育種, 2021, 19(7): 2410-2418.

    Huang H C, Jiang W X, Liang C X, et al. EST-SSR-based analysis on genetic diversity of ancient Liupao tea trees and their progeny in Guangxi [J]. Molecular Plant Breeding, 2021, 19(7): 2410-2418.

    [24] Yeh F C. Population genetic analysis of co-dominant and dominant marker and quantitative traits [J]. Belgian Journal of Botany, 1997, 130: 129-157.

    [25] 姚雪倩, 葉乃興. 雜種優(yōu)勢在茶樹育種中的應(yīng)用[J]. 茶葉學(xué)報, 2016, 57(3): 113-118.

    Yao X Q, Ye N X. Application of heterosis for breeding tea plants [J]. Acta Tea Sinica, 2016, 57(3): 113-118.

    Genetic and Phylogenetic Analysis for Germplasm Resources offrom Wuzhou City

    WANG Liubin1, HUANG Liyun2, TENG Cuiqin2, WU Liyun1, CHENG Hao1, YU Cuiping2*, WANG Liyuan1*

    1. National Centre for Tea Improvement, Tea Research Institute, Chinese Academy of Agricultural Sciences, Hangzhou 310008, China; 2. Wuzhou Institute of Agricultural Science, Wuzhou 543003, China

    Based on the SSR markers, the genetic diversity and genetic relationship of the germplasm resources of tea plants from Liubao town and Nandu town were fully analyzed in this study. The core molecular markers for the efficient identification of these germplasm resources were successfully screened. The main results show that: (1) 98 alleles were amplified from 17 pairs of SSR primers, and each pair of SSR primers amplified 3-8 alleles, with an average of 5.764?7 alleles per locus. (2) Totally 8 core markers were selected from 17 SSR markers to distinguish each germplasm resource. (3) The average number of alleles (4.647?1), genotypes per marker (7.000?0), genetic diversity (0.675?4), and the polymorphic information content (0.628?3) of native tea trees in Liubao town were higher than wild tea trees in Nandu, and close to the cultivated population. (4) Cluster analysis shows that the majority germplasm resources from Liubao town were clustered together except for several tea plants grouped with large-leaf tea cultivars from Yunnan province and a few resources were grouped into the same cluster with tea cultivars from Zhejiang and Guizhou provinces. The wild tea germplasm resources from Nandu town were grouped into the same cluster with two germplasm materials from Liubao town.In conclusion, it was showed that there are rich tea germplasmresources with high genetic diversity in Wuzhou city. This study might lay a solid foundation for the further studies to develop and utilize these tea resources.

    SSR marker, genetic diversity, Wuzhou city,

    S571.1

    A

    1000-369X(2022)05-601-09

    2022-03-23

    2022-06-13

    財政部和農(nóng)業(yè)農(nóng)村部:國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系(CARS-19)、浙江省農(nóng)業(yè)新品種選育重大專項(2021C02067-7)、創(chuàng)新驅(qū)動發(fā)展專項(桂科AA20302018-5)

    王留彬,男,博士研究生,主要從事茶樹種質(zhì)資源與育種研究。*通信作者:yucp2018@163.com;wangly@tricaas.com

    猜你喜歡
    南渡六堡茶樹
    溧陽市南渡鎮(zhèn)推進“聯(lián)防聯(lián)治聯(lián)建” 提升本質(zhì)安全水平
    參觀六堡開茶節(jié)游記
    西江月(2021年2期)2021-11-12 20:15:42
    他睡在九月的曇花里
    花火A(2019年8期)2019-11-05 09:44:20
    六堡茶機采茶園建設(shè)與管理關(guān)鍵技術(shù)
    山茶樹變身搖錢樹
    與茶相逢在六堡
    西江月(2017年1期)2017-01-17 03:31:34
    兩個推薦茶樹品種
    茶樹灣
    李清照詞在不同時期所表達的情感與心境
    淺說李清照后期詞的特點
    文教資料(2014年13期)2014-08-26 05:07:52
    国产高清视频在线播放一区| 国模一区二区三区四区视频| 亚洲乱码一区二区免费版| 人妻久久中文字幕网| 嫩草影院新地址| 日韩中字成人| 午夜福利在线观看免费完整高清在 | 一本一本综合久久| 国产精品一区二区三区四区久久| 国产精品av视频在线免费观看| 欧美日韩国产亚洲二区| 欧美丝袜亚洲另类 | 免费黄网站久久成人精品| 一边摸一边抽搐一进一小说| 国产成人福利小说| 村上凉子中文字幕在线| 成人永久免费在线观看视频| 成年免费大片在线观看| 久久亚洲真实| a级毛片a级免费在线| 黄片wwwwww| 身体一侧抽搐| 别揉我奶头 嗯啊视频| 国产一区二区在线av高清观看| 国产三级中文精品| 亚洲色图av天堂| 日韩大尺度精品在线看网址| 日本撒尿小便嘘嘘汇集6| 毛片一级片免费看久久久久 | 亚洲av第一区精品v没综合| av国产免费在线观看| 女同久久另类99精品国产91| 久久久久久久久中文| 直男gayav资源| 中亚洲国语对白在线视频| 国产成人福利小说| 久久精品国产自在天天线| 国产熟女欧美一区二区| 欧美绝顶高潮抽搐喷水| 国内久久婷婷六月综合欲色啪| 99久久精品国产国产毛片| 男女边吃奶边做爰视频| 男人舔女人下体高潮全视频| 日本成人三级电影网站| 麻豆一二三区av精品| 国产伦精品一区二区三区四那| 变态另类丝袜制服| 亚洲美女黄片视频| 国语自产精品视频在线第100页| 国产探花在线观看一区二区| 久久久久久久久久久丰满 | 亚洲男人的天堂狠狠| 精品欧美国产一区二区三| 久久精品综合一区二区三区| avwww免费| .国产精品久久| 天堂av国产一区二区熟女人妻| 女人被狂操c到高潮| 日韩亚洲欧美综合| 国产精品久久久久久av不卡| 深夜精品福利| 精华霜和精华液先用哪个| 国产免费男女视频| 天堂网av新在线| 国产在视频线在精品| 午夜久久久久精精品| 欧美中文日本在线观看视频| 一本久久中文字幕| 一个人观看的视频www高清免费观看| 露出奶头的视频| 欧美性猛交黑人性爽| 国产欧美日韩精品亚洲av| 欧美bdsm另类| 大型黄色视频在线免费观看| 天天一区二区日本电影三级| 国产精品一区二区免费欧美| 色av中文字幕| 18禁黄网站禁片午夜丰满| 中国美女看黄片| 欧美性猛交╳xxx乱大交人| 国产私拍福利视频在线观看| 亚洲精品成人久久久久久| 精品人妻偷拍中文字幕| h日本视频在线播放| netflix在线观看网站| 欧美最黄视频在线播放免费| 国产色爽女视频免费观看| 在线播放国产精品三级| 欧美另类亚洲清纯唯美| 亚洲精品456在线播放app | 国产熟女欧美一区二区| 国产白丝娇喘喷水9色精品| 国产免费男女视频| 91狼人影院| netflix在线观看网站| 免费观看精品视频网站| 国产av麻豆久久久久久久| 国产精品野战在线观看| 黄色欧美视频在线观看| 99久久成人亚洲精品观看| 亚洲最大成人手机在线| 亚洲内射少妇av| 国产极品精品免费视频能看的| 国产毛片a区久久久久| 国产免费男女视频| 中国美白少妇内射xxxbb| 尾随美女入室| 91在线观看av| 日韩在线高清观看一区二区三区 | 国产黄片美女视频| 久久久久久久久中文| 真人做人爱边吃奶动态| 久久久久久久亚洲中文字幕| av黄色大香蕉| 久久久久久久久久久丰满 | 人人妻人人看人人澡| 国产在线男女| 亚洲三级黄色毛片| 国产色爽女视频免费观看| 十八禁国产超污无遮挡网站| 日本在线视频免费播放| 又爽又黄a免费视频| 久久久久精品国产欧美久久久| 国产av一区在线观看免费| 级片在线观看| 中文字幕免费在线视频6| 亚洲av不卡在线观看| 国产v大片淫在线免费观看| 女人被狂操c到高潮| 免费看日本二区| 久久精品国产亚洲网站| 禁无遮挡网站| 亚洲精华国产精华精| 两个人视频免费观看高清| 精品午夜福利在线看| 国产伦精品一区二区三区视频9| 一本久久中文字幕| 搡老妇女老女人老熟妇| 一夜夜www| 欧美性感艳星| 我要看日韩黄色一级片| 一级毛片久久久久久久久女| 狂野欧美激情性xxxx在线观看| 听说在线观看完整版免费高清| 两个人视频免费观看高清| 国产精品自产拍在线观看55亚洲| 亚洲av成人av| 久久久久久国产a免费观看| 国产一区二区在线av高清观看| 亚洲精品粉嫩美女一区| 亚洲精品国产成人久久av| 久久精品综合一区二区三区| av在线天堂中文字幕| 舔av片在线| 狂野欧美激情性xxxx在线观看| 成年女人永久免费观看视频| 琪琪午夜伦伦电影理论片6080| 观看美女的网站| 亚洲精品影视一区二区三区av| 嫩草影院入口| 精品免费久久久久久久清纯| 综合色av麻豆| 制服丝袜大香蕉在线| 一个人免费在线观看电影| 嫩草影院入口| a在线观看视频网站| 18禁裸乳无遮挡免费网站照片| 亚洲av第一区精品v没综合| 精品久久久久久久久av| 国产一区二区激情短视频| 国产精品一区二区三区四区久久| 十八禁网站免费在线| 日本精品一区二区三区蜜桃| 久久99热这里只有精品18| 精品国产三级普通话版| 一级a爱片免费观看的视频| 国产蜜桃级精品一区二区三区| 亚洲综合色惰| 五月玫瑰六月丁香| 日韩强制内射视频| 成人午夜高清在线视频| 久久6这里有精品| 国产精品久久视频播放| 最近最新免费中文字幕在线| 春色校园在线视频观看| 午夜老司机福利剧场| 国产欧美日韩一区二区精品| 国产69精品久久久久777片| 国产av一区在线观看免费| 中文字幕久久专区| 联通29元200g的流量卡| 18+在线观看网站| 欧美日韩瑟瑟在线播放| 国产精品久久久久久久电影| 老熟妇仑乱视频hdxx| 国产午夜精品久久久久久一区二区三区 | 国模一区二区三区四区视频| 最近最新免费中文字幕在线| 男人的好看免费观看在线视频| 少妇人妻精品综合一区二区 | 欧美绝顶高潮抽搐喷水| 伦理电影大哥的女人| 精品一区二区三区视频在线观看免费| 亚洲 国产 在线| 国产主播在线观看一区二区| www.色视频.com| 亚洲熟妇中文字幕五十中出| 国产精品自产拍在线观看55亚洲| 一个人观看的视频www高清免费观看| 春色校园在线视频观看| 亚洲欧美激情综合另类| 久久草成人影院| 99久久成人亚洲精品观看| 午夜精品一区二区三区免费看| 国产精品电影一区二区三区| 成人高潮视频无遮挡免费网站| 亚洲专区国产一区二区| 欧美+亚洲+日韩+国产| 精华霜和精华液先用哪个| 成年女人永久免费观看视频| 免费搜索国产男女视频| 午夜激情欧美在线| 欧美日韩亚洲国产一区二区在线观看| 国产高清三级在线| 色5月婷婷丁香| 亚洲色图av天堂| 久久热精品热| 久久婷婷人人爽人人干人人爱| 色哟哟·www| 午夜福利在线观看吧| 亚洲人成伊人成综合网2020| 国产又黄又爽又无遮挡在线| netflix在线观看网站| 变态另类成人亚洲欧美熟女| 国产爱豆传媒在线观看| 久久这里只有精品中国| 亚洲av电影不卡..在线观看| 91精品国产九色| 亚洲av不卡在线观看| av在线老鸭窝| 麻豆一二三区av精品| 国产精品99久久久久久久久| 国产亚洲精品久久久久久毛片| 一区二区三区激情视频| 日本-黄色视频高清免费观看| av在线观看视频网站免费| 免费看a级黄色片| 日本与韩国留学比较| 亚洲av.av天堂| or卡值多少钱| 日韩高清综合在线| 日韩 亚洲 欧美在线| 精品午夜福利视频在线观看一区| 亚洲avbb在线观看| av在线亚洲专区| 动漫黄色视频在线观看| 亚洲av中文av极速乱 | 99riav亚洲国产免费| 女人被狂操c到高潮| 成人美女网站在线观看视频| 全区人妻精品视频| 嫩草影视91久久| 热99在线观看视频| 欧美一区二区国产精品久久精品| 日本爱情动作片www.在线观看 | 床上黄色一级片| 俄罗斯特黄特色一大片| 久久久久九九精品影院| 又爽又黄无遮挡网站| 淫妇啪啪啪对白视频| 日韩在线高清观看一区二区三区 | 亚洲精华国产精华精| 午夜福利高清视频| 99热这里只有精品一区| 毛片一级片免费看久久久久 | av黄色大香蕉| 久久精品国产99精品国产亚洲性色| 国产日本99.免费观看| 久久久久国内视频| 亚洲最大成人手机在线| 亚洲欧美日韩卡通动漫| 久久久久久国产a免费观看| 综合色av麻豆| 听说在线观看完整版免费高清| 天美传媒精品一区二区| 俄罗斯特黄特色一大片| 制服丝袜大香蕉在线| 少妇猛男粗大的猛烈进出视频 | 91午夜精品亚洲一区二区三区 | 91久久精品国产一区二区成人| 一夜夜www| 国产精品久久久久久av不卡| 免费观看的影片在线观看| 性色avwww在线观看| 99久久成人亚洲精品观看| 免费看日本二区| 日本三级黄在线观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲男人的天堂狠狠| 久久久精品大字幕| 在线播放无遮挡| 免费无遮挡裸体视频| 国产一区二区三区视频了| 日韩欧美 国产精品| 大型黄色视频在线免费观看| 五月伊人婷婷丁香| 久久精品国产清高在天天线| 久久中文看片网| 亚洲av成人av| 丰满乱子伦码专区| 日本免费a在线| 日韩av在线大香蕉| 精品久久久久久久末码| 禁无遮挡网站| 成人性生交大片免费视频hd| 两人在一起打扑克的视频| 国产欧美日韩一区二区精品| 日本三级黄在线观看| 免费观看在线日韩| 亚洲av.av天堂| 国产精品一区二区三区四区免费观看 | av中文乱码字幕在线| 亚洲人与动物交配视频| 午夜免费成人在线视频| 嫩草影视91久久| 午夜a级毛片| 精品免费久久久久久久清纯| 亚洲一区高清亚洲精品| 亚洲精华国产精华液的使用体验 | 亚洲欧美日韩卡通动漫| 18禁黄网站禁片午夜丰满| av在线老鸭窝| 在线观看66精品国产| 日韩精品有码人妻一区| 综合色av麻豆| 午夜爱爱视频在线播放| 国产激情偷乱视频一区二区| 欧美三级亚洲精品| 亚洲自偷自拍三级| 禁无遮挡网站| 最近最新中文字幕大全电影3| 日本熟妇午夜| 免费看a级黄色片| 嫩草影院精品99| 少妇人妻精品综合一区二区 | 熟女电影av网| 国产精品一及| 一区二区三区免费毛片| 免费黄网站久久成人精品| 天天一区二区日本电影三级| 日本与韩国留学比较| 成人特级av手机在线观看| 国产亚洲精品久久久com| 国产精品一区二区免费欧美| 成年女人看的毛片在线观看| 久久精品久久久久久噜噜老黄 | 国产精品爽爽va在线观看网站| 免费观看的影片在线观看| 大又大粗又爽又黄少妇毛片口| av.在线天堂| 亚洲熟妇中文字幕五十中出| 亚洲va在线va天堂va国产| 国产aⅴ精品一区二区三区波| 麻豆成人av在线观看| 天堂√8在线中文| 亚洲avbb在线观看| 欧美绝顶高潮抽搐喷水| 日本黄大片高清| 观看免费一级毛片| ponron亚洲| 久久久久九九精品影院| 人妻制服诱惑在线中文字幕| 少妇的逼好多水| 精品久久久久久久久久久久久| 精品久久久久久,| 三级毛片av免费| 波野结衣二区三区在线| 麻豆国产97在线/欧美| 国产aⅴ精品一区二区三区波| 深夜a级毛片| 国产人妻一区二区三区在| 日韩中文字幕欧美一区二区| 亚洲精品影视一区二区三区av| 日日摸夜夜添夜夜添av毛片 | 欧美日韩亚洲国产一区二区在线观看| 嫩草影视91久久| 久9热在线精品视频| 三级男女做爰猛烈吃奶摸视频| 超碰av人人做人人爽久久| 国产综合懂色| 成人三级黄色视频| 高清毛片免费观看视频网站| 成人永久免费在线观看视频| 亚洲精品粉嫩美女一区| 九色国产91popny在线| 久久精品夜夜夜夜夜久久蜜豆| 麻豆精品久久久久久蜜桃| 久久热精品热| 亚洲精华国产精华精| 久久久久久久久久成人| 99在线人妻在线中文字幕| 别揉我奶头 嗯啊视频| 精品人妻视频免费看| 亚洲成人久久爱视频| 国产精品1区2区在线观看.| 我要搜黄色片| 国产亚洲精品综合一区在线观看| 波多野结衣高清作品| 国产精品久久久久久久久免| 亚洲成人久久性| 亚洲国产精品成人综合色| 尾随美女入室| 变态另类丝袜制服| 一进一出抽搐gif免费好疼| 国产精品国产高清国产av| 变态另类成人亚洲欧美熟女| 亚洲熟妇熟女久久| 非洲黑人性xxxx精品又粗又长| 国产免费男女视频| 国产精品国产三级国产av玫瑰| 麻豆成人午夜福利视频| 免费一级毛片在线播放高清视频| 亚洲国产精品sss在线观看| 色精品久久人妻99蜜桃| 成人美女网站在线观看视频| 亚洲国产精品sss在线观看| 国产亚洲av嫩草精品影院| 中国美白少妇内射xxxbb| 草草在线视频免费看| 欧美xxxx性猛交bbbb| h日本视频在线播放| 最近视频中文字幕2019在线8| 99久久成人亚洲精品观看| 男人舔奶头视频| 国产伦人伦偷精品视频| 国产高清三级在线| 日本五十路高清| 欧美在线一区亚洲| 亚洲三级黄色毛片| 少妇丰满av| 国产亚洲av嫩草精品影院| 18+在线观看网站| 精品国产三级普通话版| 欧美色视频一区免费| 久久久久免费精品人妻一区二区| 久久久成人免费电影| 免费在线观看成人毛片| 国产老妇女一区| 欧美黑人欧美精品刺激| 亚洲av不卡在线观看| 亚洲av日韩精品久久久久久密| videossex国产| 赤兔流量卡办理| 99久久精品热视频| 欧美xxxx性猛交bbbb| 亚洲av二区三区四区| 欧美性猛交黑人性爽| 久久精品国产鲁丝片午夜精品 | 色哟哟·www| 日本在线视频免费播放| 老熟妇仑乱视频hdxx| 69av精品久久久久久| 可以在线观看毛片的网站| 亚洲美女视频黄频| 亚洲欧美激情综合另类| 99视频精品全部免费 在线| 欧美一区二区亚洲| 九色国产91popny在线| 欧美日韩精品成人综合77777| 日韩亚洲欧美综合| 国产毛片a区久久久久| 亚洲最大成人av| 国产又黄又爽又无遮挡在线| av天堂在线播放| 99久久精品国产国产毛片| 国产毛片a区久久久久| 国产一区二区激情短视频| 久久精品国产鲁丝片午夜精品 | 久久午夜亚洲精品久久| 看免费成人av毛片| 女生性感内裤真人,穿戴方法视频| 欧美高清性xxxxhd video| 神马国产精品三级电影在线观看| 国产亚洲精品久久久com| 少妇高潮的动态图| 亚洲狠狠婷婷综合久久图片| 国产色爽女视频免费观看| 亚洲成人中文字幕在线播放| 一级黄片播放器| 国产老妇女一区| 欧美成人性av电影在线观看| 成人精品一区二区免费| 日韩欧美三级三区| 欧美极品一区二区三区四区| 人人妻,人人澡人人爽秒播| 亚洲最大成人手机在线| h日本视频在线播放| 我要搜黄色片| 免费av观看视频| 十八禁网站免费在线| 九九爱精品视频在线观看| 午夜免费成人在线视频| 欧美色视频一区免费| 麻豆成人av在线观看| 一夜夜www| 亚洲精华国产精华液的使用体验 | 中文字幕人妻熟人妻熟丝袜美| 久久久午夜欧美精品| 久9热在线精品视频| 全区人妻精品视频| 欧美成人性av电影在线观看| 亚洲精华国产精华液的使用体验 | 精品久久久久久久久久免费视频| 日本色播在线视频| 国产成人av教育| 永久网站在线| videossex国产| 亚洲欧美日韩东京热| a级毛片a级免费在线| 亚洲欧美日韩高清在线视频| 国内揄拍国产精品人妻在线| 狂野欧美白嫩少妇大欣赏| 三级国产精品欧美在线观看| 午夜福利在线观看吧| 国产精品乱码一区二三区的特点| 色综合站精品国产| 欧美潮喷喷水| 天堂动漫精品| 观看免费一级毛片| 九九久久精品国产亚洲av麻豆| 美女高潮的动态| 91久久精品电影网| 91久久精品国产一区二区三区| 国产精品一区二区免费欧美| 人人妻,人人澡人人爽秒播| 亚洲国产色片| 亚洲成人久久性| 久久九九热精品免费| 日韩大尺度精品在线看网址| 亚洲国产精品合色在线| 舔av片在线| 欧美日本视频| 色综合站精品国产| 麻豆一二三区av精品| 99在线人妻在线中文字幕| 亚洲三级黄色毛片| 99热只有精品国产| 极品教师在线视频| 91av网一区二区| 亚洲欧美精品综合久久99| 91av网一区二区| 男人狂女人下面高潮的视频| 婷婷精品国产亚洲av在线| 简卡轻食公司| 国产女主播在线喷水免费视频网站 | 热99re8久久精品国产| 91久久精品电影网| 日韩 亚洲 欧美在线| 一级黄色大片毛片| 一本久久中文字幕| 亚洲成人免费电影在线观看| 午夜福利视频1000在线观看| 午夜福利在线观看免费完整高清在 | 亚洲图色成人| 亚洲精品成人久久久久久| 国内精品久久久久久久电影| 国产乱人视频| 午夜激情欧美在线| 91午夜精品亚洲一区二区三区 | 国产午夜福利久久久久久| 久9热在线精品视频| 国产高清视频在线播放一区| 91久久精品国产一区二区成人| 亚洲精华国产精华液的使用体验 | 久久人妻av系列| 我的女老师完整版在线观看| 国产精品一区二区免费欧美| 国产精品国产三级国产av玫瑰| 亚洲三级黄色毛片| 精品久久国产蜜桃| 真人做人爱边吃奶动态| 中文字幕熟女人妻在线| 尾随美女入室| 亚洲av第一区精品v没综合| 日韩在线高清观看一区二区三区 | 国产精品女同一区二区软件 | 国产黄a三级三级三级人| 久久久久久久久中文| 露出奶头的视频| 亚洲性久久影院| 最近视频中文字幕2019在线8| 欧美人与善性xxx| 亚洲av熟女| 变态另类成人亚洲欧美熟女| 俺也久久电影网| 日本-黄色视频高清免费观看| 精品一区二区三区视频在线| 91午夜精品亚洲一区二区三区 | 夜夜看夜夜爽夜夜摸| 精品欧美国产一区二区三| 中文字幕av成人在线电影| 桃色一区二区三区在线观看| 美女免费视频网站| 亚洲男人的天堂狠狠| 色吧在线观看| 窝窝影院91人妻| 在线免费观看的www视频| 在线看三级毛片| 五月伊人婷婷丁香| 日韩欧美国产在线观看| 国产av不卡久久| 日韩在线高清观看一区二区三区 | 啦啦啦啦在线视频资源| 无遮挡黄片免费观看| 久久草成人影院|