• <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
    99视频精品全部免费 在线| 久久6这里有精品| 狂野欧美白嫩少妇大欣赏| 真实男女啪啪啪动态图| .国产精品久久| 免费观看a级毛片全部| 国产成人福利小说| 亚洲图色成人| 国产精品99久久久久久久久| 亚洲精品亚洲一区二区| 久久久久久久午夜电影| 亚洲欧美日韩另类电影网站 | 精品人妻一区二区三区麻豆| 欧美最新免费一区二区三区| 少妇人妻久久综合中文| 欧美成人一区二区免费高清观看| 涩涩av久久男人的天堂| 全区人妻精品视频| 亚洲精品乱久久久久久| 国产在视频线精品| 久久精品国产a三级三级三级| 亚洲国产日韩一区二区| 久久精品国产自在天天线| 一级毛片久久久久久久久女| 一级av片app| 另类亚洲欧美激情| 久久99精品国语久久久| 18禁在线播放成人免费| 极品教师在线视频| 伦理电影大哥的女人| 精华霜和精华液先用哪个| 18禁在线播放成人免费| 欧美成人精品欧美一级黄| 日韩成人伦理影院| 久久国产乱子免费精品| 深爱激情五月婷婷| 91精品国产九色| av在线亚洲专区| 亚洲真实伦在线观看| a级毛片免费高清观看在线播放| 久久精品久久久久久噜噜老黄| 国产精品三级大全| 国产av国产精品国产| 亚洲人成网站高清观看| 三级国产精品片| 国产亚洲一区二区精品| 久久久久久久久久久丰满| 国产日韩欧美在线精品| 中文字幕制服av| videossex国产| 亚洲国产欧美人成| 久久久久久伊人网av| 久久精品综合一区二区三区| 亚洲av免费在线观看| 视频中文字幕在线观看| 亚洲精品色激情综合| 国产黄色免费在线视频| 亚洲av二区三区四区| 美女被艹到高潮喷水动态| 又黄又爽又刺激的免费视频.| 国产在线男女| 国产免费一区二区三区四区乱码| 亚洲精品视频女| 久久久久久久久久久丰满| 国国产精品蜜臀av免费| 伦理电影大哥的女人| 精品人妻熟女av久视频| 最新中文字幕久久久久| 国产淫语在线视频| 一级二级三级毛片免费看| 久久女婷五月综合色啪小说 | 美女内射精品一级片tv| 中文天堂在线官网| 特大巨黑吊av在线直播| 国产欧美另类精品又又久久亚洲欧美| 日日摸夜夜添夜夜添av毛片| 日韩欧美精品免费久久| 汤姆久久久久久久影院中文字幕| 久久精品国产a三级三级三级| 亚洲精品国产成人久久av| 欧美最新免费一区二区三区| 最近的中文字幕免费完整| 成人国产麻豆网| 免费看不卡的av| 男插女下体视频免费在线播放| 亚洲精品,欧美精品| av在线app专区| 激情 狠狠 欧美| 欧美区成人在线视频| 精品人妻熟女av久视频| 久久99热这里只有精品18| 80岁老熟妇乱子伦牲交| 欧美成人a在线观看| 亚洲精品亚洲一区二区| 乱码一卡2卡4卡精品| 久久人人爽av亚洲精品天堂 | 可以在线观看毛片的网站| 免费观看无遮挡的男女| 欧美精品人与动牲交sv欧美| 久久精品国产鲁丝片午夜精品| 国产精品精品国产色婷婷| 在线观看人妻少妇| 亚洲成人精品中文字幕电影| 99热这里只有是精品50| 最近中文字幕高清免费大全6| 国产乱人视频| 老师上课跳d突然被开到最大视频| 成人漫画全彩无遮挡| 国产免费一级a男人的天堂| 高清毛片免费看| 天美传媒精品一区二区| av网站免费在线观看视频| 又爽又黄a免费视频| av在线播放精品| 久久久午夜欧美精品| 99热这里只有是精品50| 大片免费播放器 马上看| 精品人妻视频免费看| 亚洲无线观看免费| 久久久久久久亚洲中文字幕| av卡一久久| 亚洲天堂av无毛| 小蜜桃在线观看免费完整版高清| 国产在视频线精品| 日韩不卡一区二区三区视频在线| 欧美人与善性xxx| 97精品久久久久久久久久精品| 最近最新中文字幕免费大全7| 人妻制服诱惑在线中文字幕| 99热这里只有是精品在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久成人免费电影| 黄色一级大片看看| 国产精品久久久久久精品古装| 国产黄a三级三级三级人| 亚洲国产高清在线一区二区三| av线在线观看网站| 亚洲国产精品成人综合色| 2021天堂中文幕一二区在线观| 免费av不卡在线播放| 日韩人妻高清精品专区| 亚洲欧美精品自产自拍| 99精国产麻豆久久婷婷| 国产白丝娇喘喷水9色精品| 十八禁网站网址无遮挡 | 精品人妻偷拍中文字幕| 免费人成在线观看视频色| 精品一区二区三区视频在线| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲内射少妇av| 欧美精品一区二区大全| 天堂网av新在线| 国产精品熟女久久久久浪| 久久久精品94久久精品| 免费观看的影片在线观看| 午夜激情福利司机影院| 99热全是精品| 天堂俺去俺来也www色官网| 99热6这里只有精品| 亚洲欧美一区二区三区国产| 天堂俺去俺来也www色官网| 国产高清有码在线观看视频| 成人国产麻豆网| 夜夜看夜夜爽夜夜摸| 精品久久久噜噜| 国产 精品1| 亚洲欧洲日产国产| 男女啪啪激烈高潮av片| 人人妻人人看人人澡| 建设人人有责人人尽责人人享有的 | 亚洲精品国产av成人精品| 日日撸夜夜添| 免费观看在线日韩| 成人综合一区亚洲| 丝袜脚勾引网站| av一本久久久久| 中文欧美无线码| 一个人看视频在线观看www免费| 久久久久久久大尺度免费视频| 六月丁香七月| 国产男人的电影天堂91| 大话2 男鬼变身卡| 免费av不卡在线播放| 国产 一区 欧美 日韩| 好男人在线观看高清免费视频| 精品久久国产蜜桃| 国产精品偷伦视频观看了| 搞女人的毛片| 直男gayav资源| 日本欧美国产在线视频| 亚洲内射少妇av| 欧美亚洲 丝袜 人妻 在线| 91狼人影院| 欧美日韩亚洲高清精品| 国产欧美日韩一区二区三区在线 | 最新中文字幕久久久久| 国产有黄有色有爽视频| 成人欧美大片| 婷婷色av中文字幕| 天天躁日日操中文字幕| 黑人高潮一二区| 国内少妇人妻偷人精品xxx网站| 午夜免费男女啪啪视频观看| 91精品国产九色| 亚洲精品国产av成人精品| 久久久久久久午夜电影| 亚洲成人中文字幕在线播放| 99精国产麻豆久久婷婷| 日韩一区二区视频免费看| 超碰97精品在线观看| 王馨瑶露胸无遮挡在线观看| 成人国产av品久久久| 麻豆国产97在线/欧美| 国产精品福利在线免费观看| 白带黄色成豆腐渣| 男人添女人高潮全过程视频| 日韩一区二区三区影片| 18+在线观看网站| 涩涩av久久男人的天堂| 最近2019中文字幕mv第一页| 午夜亚洲福利在线播放| 亚洲自拍偷在线| av福利片在线观看| 青春草视频在线免费观看| 亚洲成人一二三区av| 一区二区三区乱码不卡18| 亚洲最大成人av| www.色视频.com| 午夜日本视频在线| h日本视频在线播放| 欧美潮喷喷水| 久久久久久久国产电影| 久久久色成人| 一区二区三区四区激情视频| 亚洲av免费在线观看| 国产精品无大码| 午夜福利在线观看免费完整高清在| 91狼人影院| tube8黄色片| 美女高潮的动态| 2021少妇久久久久久久久久久| 欧美一区二区亚洲| 亚洲精品乱久久久久久| 欧美日韩在线观看h| 免费观看av网站的网址| 久久这里只有精品19| 嫩草影院入口| 老司机在亚洲福利影院| 欧美成人精品欧美一级黄| 午夜福利视频在线观看免费| 亚洲精华国产精华液的使用体验| 亚洲成人av在线免费| 人妻一区二区av| 国产精品香港三级国产av潘金莲 | 我要看黄色一级片免费的| 看免费av毛片| 欧美国产精品一级二级三级| e午夜精品久久久久久久| 午夜av观看不卡| 男的添女的下面高潮视频| 成人亚洲精品一区在线观看| 亚洲,欧美精品.| 久久精品熟女亚洲av麻豆精品| 色网站视频免费| 色吧在线观看| 亚洲,一卡二卡三卡| 五月天丁香电影| 在线亚洲精品国产二区图片欧美| 国产又爽黄色视频| 一级,二级,三级黄色视频| 亚洲情色 制服丝袜| 黑人巨大精品欧美一区二区蜜桃| 叶爱在线成人免费视频播放| 色综合欧美亚洲国产小说| 一本大道久久a久久精品| 精品国产超薄肉色丝袜足j| 国产又爽黄色视频| 精品一区二区三卡| tube8黄色片| 国产日韩欧美亚洲二区| 国产免费又黄又爽又色| 亚洲图色成人| 中文字幕亚洲精品专区| 亚洲成人一二三区av| 国产日韩欧美视频二区| 黄色视频在线播放观看不卡| 超碰成人久久| 午夜福利视频精品| 亚洲 欧美一区二区三区| 9色porny在线观看| 国产片内射在线| www.精华液| 亚洲精品国产一区二区精华液| 成人国产av品久久久| 一个人免费看片子| 亚洲国产精品国产精品| 丝袜在线中文字幕| 超碰97精品在线观看| 亚洲三区欧美一区| 欧美成人午夜精品| 亚洲欧美一区二区三区国产| 亚洲国产精品成人久久小说| 超碰成人久久| av不卡在线播放| av卡一久久| 在线观看免费高清a一片| 啦啦啦中文免费视频观看日本| 丰满迷人的少妇在线观看| 中文字幕最新亚洲高清| 男女国产视频网站| av一本久久久久| 男女高潮啪啪啪动态图| 黑人猛操日本美女一级片| 国产成人啪精品午夜网站| 高清av免费在线| 国产男人的电影天堂91| 中文精品一卡2卡3卡4更新| 国产精品人妻久久久影院| av又黄又爽大尺度在线免费看| 国产一区亚洲一区在线观看| 国产97色在线日韩免费| 中国三级夫妇交换| 大码成人一级视频| 大码成人一级视频| 欧美日韩亚洲国产一区二区在线观看 | 午夜影院在线不卡| 在线观看一区二区三区激情| 日韩视频在线欧美| 一区二区日韩欧美中文字幕| 少妇被粗大的猛进出69影院| kizo精华| 91精品三级在线观看| 精品少妇黑人巨大在线播放| 日韩,欧美,国产一区二区三区| 黄色视频不卡| 亚洲成人av在线免费| 狠狠婷婷综合久久久久久88av| 制服丝袜香蕉在线| 国产女主播在线喷水免费视频网站| av网站免费在线观看视频| 亚洲av国产av综合av卡| 如日韩欧美国产精品一区二区三区| 观看av在线不卡| 少妇人妻 视频| 成人国产av品久久久| 尾随美女入室| 捣出白浆h1v1| 国产高清不卡午夜福利| av.在线天堂| 老鸭窝网址在线观看| 美女大奶头黄色视频| av在线观看视频网站免费| 婷婷色综合大香蕉| 亚洲av男天堂| 亚洲欧洲日产国产| 久久人人爽av亚洲精品天堂| 久久 成人 亚洲| 最近手机中文字幕大全| 中文乱码字字幕精品一区二区三区| 丝袜脚勾引网站| 妹子高潮喷水视频| 国产深夜福利视频在线观看| 日本欧美视频一区| 1024香蕉在线观看| 国产在线免费精品| 成年美女黄网站色视频大全免费| 电影成人av| av福利片在线| 精品久久蜜臀av无| 成人18禁高潮啪啪吃奶动态图| 亚洲精品av麻豆狂野| 免费在线观看视频国产中文字幕亚洲 | 丝袜美足系列| 在现免费观看毛片| 国产亚洲一区二区精品| 天天躁夜夜躁狠狠躁躁| 免费少妇av软件| 亚洲一码二码三码区别大吗| 啦啦啦在线免费观看视频4| 亚洲国产精品国产精品| 国产免费现黄频在线看| 精品久久久久久电影网| 亚洲国产欧美网| 亚洲少妇的诱惑av| 欧美日韩亚洲综合一区二区三区_| 亚洲免费av在线视频| 成人18禁高潮啪啪吃奶动态图| 一区二区三区乱码不卡18| 建设人人有责人人尽责人人享有的| 人人妻,人人澡人人爽秒播 | 天堂8中文在线网| 一本久久精品| 免费看av在线观看网站| 国产精品麻豆人妻色哟哟久久| 亚洲成人国产一区在线观看 | 精品一品国产午夜福利视频| 精品酒店卫生间| 老鸭窝网址在线观看| 久久ye,这里只有精品| 熟妇人妻不卡中文字幕| 90打野战视频偷拍视频| 日日撸夜夜添| 亚洲精品日本国产第一区| 久久精品aⅴ一区二区三区四区| 麻豆av在线久日| 久久综合国产亚洲精品| 丝袜在线中文字幕| 精品免费久久久久久久清纯 | 欧美久久黑人一区二区| 婷婷色综合www| 国产欧美日韩一区二区三区在线| 欧美黑人精品巨大| 人体艺术视频欧美日本| 亚洲熟女精品中文字幕| 亚洲欧洲日产国产| av电影中文网址| 亚洲国产欧美日韩在线播放| 精品一品国产午夜福利视频| videosex国产| 街头女战士在线观看网站| 久久久国产一区二区| 亚洲婷婷狠狠爱综合网| 超碰97精品在线观看| 在线观看免费视频网站a站| 免费高清在线观看日韩| 波多野结衣一区麻豆| 黑人猛操日本美女一级片| 夜夜骑夜夜射夜夜干| 亚洲精品国产av蜜桃| 国产伦理片在线播放av一区| 天天操日日干夜夜撸| 七月丁香在线播放| 久久97久久精品| 亚洲精品成人av观看孕妇| 搡老岳熟女国产| 夫妻性生交免费视频一级片| 欧美日韩视频精品一区| 国产精品欧美亚洲77777| 高清欧美精品videossex| 伊人久久大香线蕉亚洲五| 国产精品二区激情视频| 日韩人妻精品一区2区三区| 日韩中文字幕视频在线看片| 精品少妇内射三级| 中文字幕亚洲精品专区| 深夜精品福利| 免费高清在线观看视频在线观看| 在线观看三级黄色| 天天躁日日躁夜夜躁夜夜| 久久人人97超碰香蕉20202| 精品少妇黑人巨大在线播放| 婷婷色麻豆天堂久久| 日韩av免费高清视频| 久久精品亚洲熟妇少妇任你| av在线app专区| 午夜福利免费观看在线| 国产欧美日韩综合在线一区二区| 成人18禁高潮啪啪吃奶动态图| 国产精品一区二区精品视频观看| 日韩视频在线欧美| 国产免费福利视频在线观看| 亚洲av男天堂| 一区二区三区乱码不卡18| 国产深夜福利视频在线观看| 伦理电影大哥的女人| 中文字幕人妻熟女乱码| 色视频在线一区二区三区| 在线观看www视频免费| 亚洲三区欧美一区| h视频一区二区三区| 国产精品久久久久久人妻精品电影 | 亚洲国产欧美日韩在线播放| 国产毛片在线视频| 这个男人来自地球电影免费观看 | 黄片播放在线免费| 精品国产超薄肉色丝袜足j| 久久av网站| 亚洲美女黄色视频免费看| 国产精品女同一区二区软件| 激情视频va一区二区三区| 熟女av电影| 街头女战士在线观看网站| 人人澡人人妻人| 美女视频免费永久观看网站| 啦啦啦 在线观看视频| 中文字幕人妻熟女乱码| 久久国产亚洲av麻豆专区| 日韩制服丝袜自拍偷拍| 久久国产精品大桥未久av| 一级毛片 在线播放| 久久精品亚洲熟妇少妇任你| 久久狼人影院| 亚洲精品视频女| 精品一区二区三区av网在线观看 | 国产精品国产三级专区第一集| 男女国产视频网站| 天堂俺去俺来也www色官网| 精品人妻熟女毛片av久久网站| 99精品久久久久人妻精品| 毛片一级片免费看久久久久| 亚洲av成人不卡在线观看播放网 | 国产成人欧美在线观看 | 乱人伦中国视频| 精品免费久久久久久久清纯 | 2021少妇久久久久久久久久久| 巨乳人妻的诱惑在线观看| 亚洲国产欧美一区二区综合| 青春草亚洲视频在线观看| 欧美成人精品欧美一级黄| 国产精品一区二区在线观看99| 色吧在线观看| 日韩 亚洲 欧美在线| 18禁国产床啪视频网站| 哪个播放器可以免费观看大片| 少妇 在线观看| 免费高清在线观看视频在线观看| 欧美日韩av久久| 叶爱在线成人免费视频播放| 国产99久久九九免费精品| 久久久精品免费免费高清| 七月丁香在线播放| 如何舔出高潮| 亚洲av国产av综合av卡| 亚洲四区av| 日本黄色日本黄色录像| 中国三级夫妇交换| 麻豆精品久久久久久蜜桃| 波野结衣二区三区在线| av网站在线播放免费| 啦啦啦啦在线视频资源| 国产熟女欧美一区二区| av一本久久久久| 免费黄网站久久成人精品| 成人毛片60女人毛片免费| 亚洲国产精品一区三区| 亚洲国产精品国产精品| 日韩人妻精品一区2区三区| 啦啦啦视频在线资源免费观看| 亚洲国产毛片av蜜桃av| 只有这里有精品99| netflix在线观看网站| 考比视频在线观看| 亚洲欧美成人综合另类久久久| av女优亚洲男人天堂| 久久久久精品人妻al黑| 国产人伦9x9x在线观看| 亚洲,欧美精品.| 成人影院久久| 国产av一区二区精品久久| 高清av免费在线| 午夜影院在线不卡| e午夜精品久久久久久久| 国产欧美日韩一区二区三区在线| 欧美日韩成人在线一区二区| 水蜜桃什么品种好| 满18在线观看网站| 欧美 亚洲 国产 日韩一| 天堂中文最新版在线下载| 丁香六月欧美| 午夜免费男女啪啪视频观看| av片东京热男人的天堂| 丰满饥渴人妻一区二区三| 精品免费久久久久久久清纯 | 亚洲欧美激情在线| 在线天堂最新版资源| 国产日韩一区二区三区精品不卡| 午夜福利视频在线观看免费| 国产一区二区三区av在线| 超碰成人久久| 国产在线免费精品| av一本久久久久| 午夜福利视频在线观看免费| 制服诱惑二区| 丝袜脚勾引网站| 亚洲成人av在线免费| 青青草视频在线视频观看| 免费少妇av软件| 国产精品秋霞免费鲁丝片| 黄网站色视频无遮挡免费观看| 天堂中文最新版在线下载| 久久 成人 亚洲| 欧美老熟妇乱子伦牲交| 中文精品一卡2卡3卡4更新| 自线自在国产av| 成人手机av| 最近2019中文字幕mv第一页| 99国产综合亚洲精品| 久久青草综合色| 免费看不卡的av| 亚洲av福利一区| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩成人在线一区二区| 亚洲精品日本国产第一区| 免费高清在线观看日韩| 黄色怎么调成土黄色| 久久天堂一区二区三区四区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av在线观看美女高潮| 中文欧美无线码| 久热爱精品视频在线9| 最近最新中文字幕大全免费视频 | 男女边吃奶边做爰视频| av在线老鸭窝| 国产成人91sexporn| 国产欧美日韩一区二区三区在线| 久久99精品国语久久久| 水蜜桃什么品种好| 亚洲七黄色美女视频| 一本久久精品| 国产男女超爽视频在线观看| 女的被弄到高潮叫床怎么办| 超色免费av| 韩国av在线不卡| 亚洲精品自拍成人| 亚洲欧美精品综合一区二区三区|