• <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在线天堂中文字幕| 色综合站精品国产| 91在线精品国自产拍蜜月| 亚洲精品影视一区二区三区av| 午夜老司机福利剧场| 国产精品美女特级片免费视频播放器| 能在线免费观看的黄片| 嘟嘟电影网在线观看| 国产久久久一区二区三区| 国产综合懂色| 国产日韩欧美在线精品| 真实男女啪啪啪动态图| 午夜视频国产福利| 伊人久久精品亚洲午夜| 在线观看66精品国产| 一进一出抽搐gif免费好疼| 久久99蜜桃精品久久| 国产精品久久久久久精品电影| 麻豆国产97在线/欧美| 春色校园在线视频观看| 久久久久久久亚洲中文字幕| 麻豆久久精品国产亚洲av| 91麻豆精品激情在线观看国产| 亚洲aⅴ乱码一区二区在线播放| 夜夜夜夜夜久久久久| 精品久久久久久久人妻蜜臀av| 永久网站在线| 久久久久久久久久黄片| 久久久色成人| 99久久精品一区二区三区| 不卡一级毛片| 一边亲一边摸免费视频| 国产精品美女特级片免费视频播放器| 亚洲国产日韩欧美精品在线观看| 自拍偷自拍亚洲精品老妇| 国产精品日韩av在线免费观看| 日韩亚洲欧美综合| 亚洲丝袜综合中文字幕| 1000部很黄的大片| 亚州av有码| 18禁裸乳无遮挡免费网站照片| 免费黄网站久久成人精品| 日本熟妇午夜| 五月伊人婷婷丁香| 亚洲五月天丁香| 欧美极品一区二区三区四区| 国产精品一区二区在线观看99 | 亚洲精华国产精华液的使用体验 | 欧美另类亚洲清纯唯美| av福利片在线观看| 一区福利在线观看| 成人无遮挡网站| 亚洲图色成人| 变态另类成人亚洲欧美熟女| www.色视频.com| 卡戴珊不雅视频在线播放| 桃色一区二区三区在线观看| 看免费成人av毛片| 边亲边吃奶的免费视频| 国产精品,欧美在线| 午夜亚洲福利在线播放| 免费大片18禁| 18禁在线无遮挡免费观看视频| 精品一区二区三区视频在线| 日韩高清综合在线| 禁无遮挡网站| 性色avwww在线观看| 一区二区三区免费毛片| 激情 狠狠 欧美| 国产黄a三级三级三级人| 久久久久久久午夜电影| 爱豆传媒免费全集在线观看| 国产高清有码在线观看视频| 春色校园在线视频观看| 男人的好看免费观看在线视频| 中文字幕精品亚洲无线码一区| 99热全是精品| 蜜桃久久精品国产亚洲av| 精品久久久久久久人妻蜜臀av| 欧美人与善性xxx| 精品不卡国产一区二区三区| 亚洲精品乱码久久久久久按摩| 日韩欧美精品v在线| 日日撸夜夜添| 国产一区二区亚洲精品在线观看| 欧美激情在线99| 国产精品一区www在线观看| 国产单亲对白刺激| 又黄又爽又刺激的免费视频.| 久久久欧美国产精品| 级片在线观看| 国产精品野战在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产老妇女一区| av天堂在线播放| avwww免费| 国产精品一及| 中国国产av一级| 亚洲精品久久国产高清桃花| 国产人妻一区二区三区在| 中文字幕制服av| 亚洲丝袜综合中文字幕| 亚洲性久久影院| 91av网一区二区| 亚洲欧美日韩卡通动漫| 日本三级黄在线观看| 黄色一级大片看看| 国产高清不卡午夜福利| 精品日产1卡2卡| 免费黄网站久久成人精品| 中文字幕久久专区| 国语自产精品视频在线第100页| 天堂中文最新版在线下载 | 大香蕉久久网| 亚洲一区二区三区色噜噜| 久久精品国产自在天天线| 国产一区亚洲一区在线观看| 欧美zozozo另类| 老司机影院成人| 国产在线精品亚洲第一网站| 中文欧美无线码| 夜夜爽天天搞| 少妇的逼水好多| av天堂在线播放| 岛国在线免费视频观看| 尤物成人国产欧美一区二区三区| 欧美又色又爽又黄视频| 能在线免费看毛片的网站| 亚洲三级黄色毛片| 日本五十路高清| 一夜夜www| 国产久久久一区二区三区| 成年女人看的毛片在线观看| 美女xxoo啪啪120秒动态图| 国产精品蜜桃在线观看 | 国产精品精品国产色婷婷| 国产一级毛片在线| 午夜福利在线观看免费完整高清在 | 日韩成人av中文字幕在线观看| 亚洲人成网站在线播| 小说图片视频综合网站| 国产视频首页在线观看| 99久久人妻综合| 国产免费一级a男人的天堂| 亚洲精品国产成人久久av| 男人狂女人下面高潮的视频| 在线天堂最新版资源| 国产精品日韩av在线免费观看| 免费看av在线观看网站| 婷婷精品国产亚洲av| 国产成年人精品一区二区| 国产中年淑女户外野战色| 欧美一级a爱片免费观看看| 岛国毛片在线播放| 又爽又黄无遮挡网站| 亚洲欧美精品综合久久99| 国产探花极品一区二区| 亚洲欧洲国产日韩| 亚洲精品影视一区二区三区av| 蜜桃久久精品国产亚洲av| 精品不卡国产一区二区三区| 国产伦精品一区二区三区视频9| 综合色丁香网| 中文字幕精品亚洲无线码一区| 亚洲欧美成人精品一区二区| av天堂在线播放| 亚洲七黄色美女视频| 高清毛片免费看| 亚洲中文字幕日韩| 中文精品一卡2卡3卡4更新| 天天躁日日操中文字幕| 晚上一个人看的免费电影| 久久鲁丝午夜福利片| 黄片无遮挡物在线观看| 我要看日韩黄色一级片| 日本撒尿小便嘘嘘汇集6| 赤兔流量卡办理| 一个人看视频在线观看www免费| 亚洲第一区二区三区不卡| 国产午夜精品久久久久久一区二区三区| 亚州av有码| 成人无遮挡网站| 国产精品久久久久久久久免| 啦啦啦韩国在线观看视频| 欧美日韩综合久久久久久| 亚洲国产精品sss在线观看| 国国产精品蜜臀av免费| 国产淫片久久久久久久久| 精品久久久噜噜| 亚洲第一区二区三区不卡| 狂野欧美激情性xxxx在线观看| 国产一级毛片在线| 中文字幕av在线有码专区| 国产v大片淫在线免费观看| 麻豆乱淫一区二区| or卡值多少钱| 色播亚洲综合网| 男人和女人高潮做爰伦理| 国产精品麻豆人妻色哟哟久久 | 一级毛片我不卡| 久久久久久久久久久丰满| 寂寞人妻少妇视频99o| 在线免费观看的www视频| 一本一本综合久久| 日韩一区二区视频免费看| 成年女人看的毛片在线观看| 日日啪夜夜撸| 舔av片在线| 在线观看午夜福利视频| 久久草成人影院| 国产成人午夜福利电影在线观看| 欧美精品一区二区大全| 亚洲精品成人久久久久久| 久久久精品94久久精品| 在线观看午夜福利视频| 蜜桃久久精品国产亚洲av| 91狼人影院| 99国产极品粉嫩在线观看| 最近中文字幕高清免费大全6| 国产 一区 欧美 日韩| 黄片无遮挡物在线观看| 亚洲四区av| 久久久久久久久久成人| 色综合站精品国产| 国产综合懂色| 精品国产三级普通话版| a级毛色黄片| 99视频精品全部免费 在线| 国产一区二区三区在线臀色熟女| 国产 一区 欧美 日韩| 欧美日韩在线观看h| 日本欧美国产在线视频| 免费无遮挡裸体视频| a级毛片a级免费在线| 久久精品国产鲁丝片午夜精品| 亚洲精品久久久久久婷婷小说 | 亚洲国产色片| 特大巨黑吊av在线直播| 18+在线观看网站| 国产乱人视频| 亚洲av成人av| 秋霞在线观看毛片| 免费av观看视频| 亚洲精品日韩av片在线观看| 日本爱情动作片www.在线观看| 国产91av在线免费观看| 夫妻性生交免费视频一级片| 91麻豆精品激情在线观看国产| 日本欧美国产在线视频| 免费一级毛片在线播放高清视频| 中文字幕av成人在线电影| 综合色av麻豆| 色视频www国产| 麻豆成人av视频| 亚洲精品日韩av片在线观看| 女人十人毛片免费观看3o分钟| 久久草成人影院| 国内精品一区二区在线观看| 国内揄拍国产精品人妻在线| 91麻豆精品激情在线观看国产| 欧美最新免费一区二区三区| 亚洲欧美日韩无卡精品| 卡戴珊不雅视频在线播放| 久久久精品大字幕| 亚洲精品亚洲一区二区| 美女xxoo啪啪120秒动态图| 欧美色欧美亚洲另类二区| 搞女人的毛片| 哪个播放器可以免费观看大片| 日日干狠狠操夜夜爽| 欧美高清性xxxxhd video| 亚洲国产精品合色在线| 99久久久亚洲精品蜜臀av| 国产成年人精品一区二区| 麻豆国产av国片精品| 看黄色毛片网站| 亚洲精品亚洲一区二区| 欧美性感艳星| 欧美成人免费av一区二区三区| 日本免费一区二区三区高清不卡| 久久久久网色| 午夜福利在线观看免费完整高清在 | 亚洲人与动物交配视频| 12—13女人毛片做爰片一| 18禁裸乳无遮挡免费网站照片| 免费看光身美女| 狂野欧美激情性xxxx在线观看| 天堂网av新在线| 国产精品久久久久久精品电影| 中出人妻视频一区二区| 亚洲国产精品sss在线观看| 午夜激情福利司机影院| 成人毛片a级毛片在线播放| 久久久久九九精品影院| 国产亚洲精品久久久久久毛片| 最近中文字幕高清免费大全6| 久久久欧美国产精品| 国产国拍精品亚洲av在线观看| 久久久久久久久久久丰满| 亚洲精品成人久久久久久| 国产日本99.免费观看| 91久久精品电影网| av女优亚洲男人天堂| 小说图片视频综合网站| 中国美女看黄片| 免费看光身美女| 欧美性猛交╳xxx乱大交人| 亚洲三级黄色毛片| 亚洲18禁久久av| 少妇熟女欧美另类| 亚洲色图av天堂| 狂野欧美激情性xxxx在线观看| 亚洲欧美精品综合久久99| 成年女人看的毛片在线观看| 给我免费播放毛片高清在线观看| 精品熟女少妇av免费看| 一本久久中文字幕| 人体艺术视频欧美日本| 中国美女看黄片| 免费无遮挡裸体视频| 狠狠狠狠99中文字幕| 国产精品久久久久久精品电影小说 | 亚洲精品色激情综合| 女人被狂操c到高潮| 国产成人精品一,二区 | 免费看光身美女| 床上黄色一级片| 亚洲精品国产av成人精品| 国产91av在线免费观看| 欧美三级亚洲精品| 12—13女人毛片做爰片一| 激情 狠狠 欧美| 国产黄色小视频在线观看| 精品久久久久久久久久久久久| 一级黄色大片毛片| 国产成人福利小说| 如何舔出高潮| 2022亚洲国产成人精品| 一本久久中文字幕| 国产黄片美女视频| 国产一区二区三区在线臀色熟女| 99久久精品一区二区三区| 97超视频在线观看视频| 亚洲国产精品久久男人天堂| 18禁黄网站禁片免费观看直播| 1024手机看黄色片| 色哟哟哟哟哟哟| 成人无遮挡网站| 亚洲最大成人中文| 久久综合国产亚洲精品| 国产 一区 欧美 日韩| 男女下面进入的视频免费午夜| 国产精品一二三区在线看| 在线天堂最新版资源| 国产精品久久久久久av不卡| 人妻制服诱惑在线中文字幕| 69人妻影院| 中国美白少妇内射xxxbb| 亚洲精品色激情综合| 国产乱人偷精品视频| 少妇丰满av| 国产精品久久久久久久久免| 日韩av不卡免费在线播放| 好男人在线观看高清免费视频| 成人亚洲欧美一区二区av| 欧美在线一区亚洲| 国产成人freesex在线| 国产一区二区激情短视频| 国产成人精品一,二区 | 欧美三级亚洲精品| 久久久色成人| 国产女主播在线喷水免费视频网站 | 久久久久久久久久成人| 美女xxoo啪啪120秒动态图| 少妇的逼水好多| 麻豆国产97在线/欧美| 国内精品久久久久精免费| 97超视频在线观看视频| 啦啦啦观看免费观看视频高清| 亚洲av免费在线观看| 日韩欧美国产在线观看| 亚州av有码| 久久久a久久爽久久v久久| 国产精品久久久久久精品电影| 欧美高清性xxxxhd video| www日本黄色视频网| 又爽又黄无遮挡网站| 久久人人精品亚洲av| 麻豆精品久久久久久蜜桃| 日韩精品青青久久久久久| 精品无人区乱码1区二区| 日本免费一区二区三区高清不卡| 国产麻豆成人av免费视频| 午夜福利高清视频| 美女 人体艺术 gogo| 国产精品久久久久久亚洲av鲁大| 麻豆精品久久久久久蜜桃| 自拍偷自拍亚洲精品老妇| 人妻系列 视频| 国产极品天堂在线| 干丝袜人妻中文字幕| 日本与韩国留学比较| 99热这里只有是精品在线观看| 日本五十路高清| 欧美色欧美亚洲另类二区| 亚洲真实伦在线观看| 久久精品国产亚洲网站| 国产午夜精品久久久久久一区二区三区| 99热这里只有是精品在线观看| 国产午夜精品久久久久久一区二区三区| 干丝袜人妻中文字幕| av国产免费在线观看| 亚洲欧美精品综合久久99| 婷婷色综合大香蕉| 国产乱人偷精品视频| 丰满人妻一区二区三区视频av| 国产又黄又爽又无遮挡在线| 亚洲国产欧美在线一区| av专区在线播放| 国产精品三级大全| 日日摸夜夜添夜夜添av毛片| 校园人妻丝袜中文字幕| 看十八女毛片水多多多| 97热精品久久久久久| 九九在线视频观看精品| 亚洲国产精品成人久久小说 | 女的被弄到高潮叫床怎么办| 热99在线观看视频| 国产亚洲av片在线观看秒播厂 | 寂寞人妻少妇视频99o| 一本一本综合久久| 中文字幕av成人在线电影| 在线免费观看不下载黄p国产| 免费看a级黄色片| 啦啦啦啦在线视频资源| 国产一区二区在线观看日韩| 色尼玛亚洲综合影院| 高清日韩中文字幕在线| 精品国内亚洲2022精品成人| a级毛片a级免费在线| 在线播放无遮挡| 免费大片18禁| 成年版毛片免费区| 午夜福利视频1000在线观看| 日韩欧美 国产精品| 99热全是精品| 免费观看a级毛片全部| 国产v大片淫在线免费观看| 少妇被粗大猛烈的视频| 又爽又黄a免费视频| 国产成人精品婷婷| 免费av观看视频| 综合色丁香网| 亚洲经典国产精华液单| 久久国内精品自在自线图片| 此物有八面人人有两片| 蜜桃亚洲精品一区二区三区| 最新中文字幕久久久久| 久久久a久久爽久久v久久| 国产黄a三级三级三级人| 午夜免费男女啪啪视频观看| 一个人免费在线观看电影| 少妇高潮的动态图| 高清毛片免费观看视频网站| 一进一出抽搐gif免费好疼| 特级一级黄色大片| 成人午夜高清在线视频| 人体艺术视频欧美日本| 在线免费十八禁| 日日啪夜夜撸| 人妻久久中文字幕网| 美女国产视频在线观看| 国产精品电影一区二区三区| 一级二级三级毛片免费看| 少妇熟女aⅴ在线视频| 亚洲五月天丁香| 国产午夜福利久久久久久| 午夜爱爱视频在线播放| 亚洲成人中文字幕在线播放| 高清在线视频一区二区三区 | 狠狠狠狠99中文字幕| 波野结衣二区三区在线| 精品欧美国产一区二区三| 91久久精品电影网| 久久精品久久久久久久性| 在线观看一区二区三区| 国产成人精品一,二区 | 成人亚洲欧美一区二区av| 麻豆成人午夜福利视频| 日本色播在线视频| 级片在线观看| 国产一区二区激情短视频| 国产探花极品一区二区| 能在线免费看毛片的网站| 国产爱豆传媒在线观看| 一级av片app| 高清日韩中文字幕在线| 在线国产一区二区在线| 性插视频无遮挡在线免费观看| 精品熟女少妇av免费看| 亚洲无线观看免费| 亚洲精品久久国产高清桃花| ponron亚洲| 中出人妻视频一区二区| 美女黄网站色视频| 亚洲精华国产精华液的使用体验 | 国产极品精品免费视频能看的| 一进一出抽搐动态| 午夜a级毛片| 联通29元200g的流量卡| 一级毛片久久久久久久久女| 一级黄色大片毛片| 非洲黑人性xxxx精品又粗又长| 久久精品久久久久久噜噜老黄 | 国产高潮美女av| 美女高潮的动态| 久久99热这里只有精品18| 大又大粗又爽又黄少妇毛片口| 99九九线精品视频在线观看视频| .国产精品久久| 可以在线观看的亚洲视频| 亚洲国产精品sss在线观看| 国产色爽女视频免费观看| 亚洲精品国产成人久久av| 久久久久国产网址| 午夜精品在线福利| 欧美日本亚洲视频在线播放| 日本在线视频免费播放| 床上黄色一级片| 日韩国内少妇激情av| 亚洲最大成人手机在线| 国产 一区精品| 国产成人a区在线观看| 国产精品精品国产色婷婷| 国产精品1区2区在线观看.| 国产真实乱freesex| 国产真实伦视频高清在线观看| 日韩欧美国产在线观看| 久久精品综合一区二区三区| 黄色日韩在线| 少妇熟女欧美另类| 99久久精品热视频| 麻豆一二三区av精品| 女人十人毛片免费观看3o分钟| 国产极品精品免费视频能看的| 性欧美人与动物交配| 久久中文看片网| 在线播放无遮挡| 久久久久九九精品影院| av在线观看视频网站免费| 直男gayav资源| 成年女人永久免费观看视频| 亚洲av二区三区四区| 亚洲欧美中文字幕日韩二区| 日韩中字成人| 国产精品野战在线观看| 69av精品久久久久久| 精品人妻偷拍中文字幕| 精品99又大又爽又粗少妇毛片| 国产亚洲91精品色在线| 大香蕉久久网| 国产美女午夜福利| 搡女人真爽免费视频火全软件| 成人午夜精彩视频在线观看| 男女视频在线观看网站免费| 日本免费一区二区三区高清不卡| 亚洲不卡免费看| 色综合站精品国产| 亚洲av二区三区四区| 99热全是精品| 在线免费观看的www视频| av在线亚洲专区| 久久久久九九精品影院| 成人三级黄色视频| 黄色欧美视频在线观看| 欧美在线一区亚洲| 日韩精品青青久久久久久| 亚洲天堂国产精品一区在线| 一级二级三级毛片免费看| 天堂中文最新版在线下载 | 男女啪啪激烈高潮av片| 别揉我奶头 嗯啊视频| 欧美性猛交黑人性爽| 夫妻性生交免费视频一级片| 久久精品国产亚洲网站| 亚洲精品色激情综合| 九草在线视频观看| 精品人妻熟女av久视频| 亚洲av男天堂| av在线亚洲专区| 好男人在线观看高清免费视频| 午夜a级毛片| 国产精品不卡视频一区二区| 亚洲欧美清纯卡通| 中文欧美无线码| 99精品在免费线老司机午夜| 亚洲精品影视一区二区三区av| 亚洲av中文av极速乱| av在线老鸭窝| 69人妻影院| 日韩精品青青久久久久久| av在线老鸭窝| 欧美在线一区亚洲| 2022亚洲国产成人精品| 亚洲人成网站高清观看| 欧美在线一区亚洲| 亚洲精品久久久久久婷婷小说 | 日日摸夜夜添夜夜添av毛片| 在现免费观看毛片| 3wmmmm亚洲av在线观看| kizo精华| 久久九九热精品免费|