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

    茶樹地方群體種資源葉片表型及生化組分多樣性分析

    2023-08-30 14:19:02唐璐李長(zhǎng)樂葛悅王璞趙華王明樂王郁郭飛倪德江
    茶葉科學(xué) 2023年4期
    關(guān)鍵詞:咖啡堿表型茶樹

    唐璐,李長(zhǎng)樂,葛悅,王璞,趙華,王明樂,王郁,郭飛,倪德江

    茶樹地方群體種資源葉片表型及生化組分多樣性分析

    唐璐,李長(zhǎng)樂,葛悅,王璞,趙華,王明樂,王郁,郭飛*,倪德江

    華中農(nóng)業(yè)大學(xué)園藝林學(xué)學(xué)院果蔬園藝作物種質(zhì)創(chuàng)新與利用全國(guó)重點(diǎn)實(shí)驗(yàn)室,湖北 武漢 430070

    為保護(hù)茶樹種質(zhì)資源和促進(jìn)茶樹種質(zhì)創(chuàng)新,對(duì)來自全國(guó)12個(gè)省份的32份茶樹地方群體種資源表型性狀及生化組分多樣性進(jìn)行了研究。結(jié)果表明,16個(gè)表型性狀的遺傳多樣性指數(shù)變化范圍為0.81~1.94,供試群體表型遺傳多樣性較高;葉面積、葉形、葉面、芽葉茸毛、葉尖、葉身6個(gè)表型性狀變異的累計(jì)貢獻(xiàn)率為71.78%,是造成供試茶樹群體表型差異的主要因素;綜合兩年生化成分測(cè)定結(jié)果,供試樣品咖啡堿含量變異系數(shù)(17.95%和14.55%)最大,其次是茶多酚含量變異系數(shù)(13.61%和8.11%),游離氨基酸變異系數(shù)(5.62%和7.52%)最??;基于生化組分含量測(cè)定結(jié)果進(jìn)行聚類分析將供試茶樹群體劃分為3個(gè)類群,分別包含不同類型的種質(zhì)。從供試茶樹地方群體種資源的地理分布范圍來看,供試種質(zhì)生化成分含量呈現(xiàn)出地域差異。

    茶樹地方群體種;表型性狀;生化組分;多樣性

    茶樹[(L.)O. Kuntze],作為一種經(jīng)濟(jì)價(jià)值極高的葉用作物,被廣泛栽培于國(guó)內(nèi)多個(gè)省份。近年來,隨著茶產(chǎn)業(yè)的蓬勃發(fā)展,對(duì)優(yōu)質(zhì)、特異、多元茶樹種質(zhì)資源的需求日益增長(zhǎng),發(fā)掘和保護(hù)優(yōu)異茶樹種質(zhì)資源的重要性日益凸顯。茶樹地方群體種對(duì)其所在茶區(qū)生態(tài)條件和生產(chǎn)模式具有良好的適應(yīng)性,經(jīng)過長(zhǎng)期的自然或人工選擇育種,為茶樹品種改良提供了豐富的資源。當(dāng)前育成的無性系茶樹品種多是從茶樹地方群體種中采用單株系統(tǒng)選育而來,如龍井43[(L.)O. Kuntze cv.]是從浙江龍井群體種中單株系統(tǒng)選育而來。因此,茶樹地方群體種既可作為育種材料,從中篩選出具有特異性狀(表型性狀特異、經(jīng)濟(jì)性狀特異)的新品種應(yīng)用于生產(chǎn),也是研究不同地域茶樹親緣關(guān)系、傳播方向、演化規(guī)律的重要材料,對(duì)其開展遺傳多樣性研究具有重要意義。

    開展茶樹遺傳多樣性研究是保護(hù)、開發(fā)、利用茶樹種質(zhì)資源的前提,是茶樹創(chuàng)新育種和基因遺傳改良的基礎(chǔ)[1]。對(duì)茶樹表型性狀進(jìn)行測(cè)量、統(tǒng)計(jì)是最早的遺傳多樣性分析方法,該方法具有直觀、操作簡(jiǎn)單、成本低的優(yōu)勢(shì),從表型上便可直接選擇到與目標(biāo)性狀相關(guān)聯(lián)的個(gè)體或群體[2]。茶樹作為葉用經(jīng)濟(jì)作物,其葉片性狀觀測(cè)在茶樹種質(zhì)資源鑒定中具有重要地位[3]。植物在逆境脅迫下,出于自我保護(hù)會(huì)啟動(dòng)自身的抗逆機(jī)制來適應(yīng)環(huán)境變化,長(zhǎng)期的逆境脅迫會(huì)產(chǎn)生不可逆的植物表型變異,基于表型的遺傳多樣性研究對(duì)探究植物對(duì)環(huán)境的適應(yīng)性及其進(jìn)化方向都具有重要意義[4]。以植物主要代謝產(chǎn)物差異為基礎(chǔ)進(jìn)行遺傳多樣性研究的方法[5]有助于指導(dǎo)茶樹育種,培育特異種質(zhì),運(yùn)用的分析方法主要有光譜分析法和色譜分析法。茶多酚、咖啡堿、游離氨基酸是茶葉中主要的次生代謝產(chǎn)物,這些生化組分的含量與茶葉品質(zhì)息息相關(guān),各物質(zhì)之間的比例差異使不同茶樹品種適制性不同。劉彤[6]、莫嵐等[7]分別對(duì)廣西省內(nèi)柳州及來賓野生茶樹群體的主要生化特征進(jìn)行了鑒定和遺傳多樣性分析,通過聚類分析,揭示了不同茶樹群體的適制性和品質(zhì)特征。對(duì)群體種資源的研究有利于發(fā)掘生化成分特異的植株。有學(xué)者從廣東省茶樹群體種中發(fā)現(xiàn)了具有高茶多酚、高氨基酸及高兒茶素含量的群體[8-9];對(duì)四川地區(qū)茶樹資源的生化組分表型多樣性研究發(fā)現(xiàn),該地區(qū)茶樹資源的水浸出物含量較高[10];在對(duì)廣西省不同區(qū)域茶樹群體種的研究中也發(fā)現(xiàn)桂南區(qū)域種質(zhì)的茶多酚、水浸出物及咖啡堿含量較高,桂北區(qū)域種質(zhì)的氨基酸含量較高[11]。因此,通過研究不同茶樹品種間葉片表型與生化組分的差異,能從不同角度闡釋茶樹的遺傳多樣性水平[12]。

    目前茶樹種質(zhì)資源方面的研究主要集中于某一區(qū)域(省份),對(duì)不同區(qū)域(省份)間茶樹群體種資源的綜合比較研究較少,不利于直觀探究不同地域茶樹群體種之間的相關(guān)性與差異性[13]。本研究搜集了來自全國(guó)12個(gè)省份的32份茶樹群體品種資源,對(duì)其新梢芽葉、成熟葉片表型數(shù)據(jù)進(jìn)行測(cè)量統(tǒng)計(jì),采用主成分分析和聚類分析法對(duì)葉片表型性狀的多樣性展開研究,揭示了不同地區(qū)茶樹群體種可能的遺傳背景;測(cè)定樣品的主要生化成分游離氨基酸、咖啡堿、茶多酚含量,并進(jìn)行遺傳多樣性分析及聚類分析,揭示了不同地區(qū)茶樹種質(zhì)資源生化成分含量規(guī)律?;诒硇托誀罘治龊蜕M分含量的分析開展茶樹群體種資源遺傳多樣性研究,以期深入了解不同類型的茶樹地方群體種,明確國(guó)內(nèi)群體種質(zhì)的區(qū)域差異與遺傳結(jié)構(gòu),為其演化過程中茶樹表型性狀、生化成分含量的變化研究提供理論依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)材料

    供試材料為在全國(guó)12個(gè)省份收集的32份茶樹地方群體種資源,均為當(dāng)?shù)鼐哂幸欢ㄓ绊懥Φ娜后w種資源。供試地方群體種茶樹經(jīng)實(shí)地考察、多人鑒定后于2015年統(tǒng)一收集種子并播種在營(yíng)養(yǎng)缽內(nèi),隨后于2016年3月定植于華中農(nóng)業(yè)大學(xué)茶學(xué)試驗(yàn)基地茶園種質(zhì)資源圃。供試茶樹種質(zhì)的名稱、來源、編號(hào)等信息見表1。試驗(yàn)基地位于東經(jīng)114°35′,北緯30°47′,海拔29.5?m,年平均溫度17.9?℃,年均降雨量2?012?mm,相對(duì)濕度70.3%,屬亞熱帶季風(fēng)氣候。茶園種質(zhì)資源圃行距1?m,地勢(shì)平緩,立地條件和栽培管理措施基本一致。

    1.2 茶樹葉片表型性狀觀測(cè)方法

    參照陳亮等[14]編著的《茶樹種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)》,對(duì)32份供試群體種茶樹的新梢芽葉及成熟葉片表型性狀進(jìn)行觀測(cè)、統(tǒng)計(jì)。從每份群體種資源中隨機(jī)選取10株作為觀測(cè)材料,共計(jì)觀測(cè)320株茶樹種質(zhì)。性狀觀測(cè)由2人同時(shí)進(jìn)行,對(duì)茶樹的芽葉色澤、芽葉茸毛、葉片著生狀態(tài)、葉色、葉身、葉質(zhì)等16個(gè)表型性狀進(jìn)行觀測(cè)、記錄,12個(gè)質(zhì)量性狀描述代碼及4個(gè)數(shù)量性狀統(tǒng)計(jì)方法如表2所示。

    1.3 茶樹種質(zhì)生化成分檢測(cè)方法

    生化成分含量測(cè)定樣品取樣采用混樣法,分別于2019年4月和2020年7月在供試茶樹三分之二植株新梢發(fā)育為一芽二葉(平展)期采摘其新梢一芽二葉,為避免采摘樣品受種植區(qū)域影響,每份種質(zhì)分3個(gè)區(qū)域,每個(gè)區(qū)域內(nèi)隨機(jī)混合采摘1份樣品,即進(jìn)行3次生物學(xué)重復(fù),采摘下的樣品用液氮處理后于–80?℃超低溫冰箱預(yù)冷,隨后真空冷凍干燥固樣。

    生化成分檢測(cè)方法采用紫外吸光光度法,干物質(zhì)含量測(cè)定參照GB/T 8303—2013;游離氨基酸總量測(cè)定參照GB/T 8314—2013;咖啡堿含量測(cè)定參照GB/T 8312—2013;茶多酚含量測(cè)定參照GB/T 8313—2018。每份樣品均進(jìn)行3次技術(shù)重復(fù)。

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

    根據(jù)田間統(tǒng)計(jì)結(jié)果及質(zhì)量性狀的賦值標(biāo)準(zhǔn)將質(zhì)量性狀數(shù)值化,使用Excel 2019整理并計(jì)算表型性狀的分布頻率、Shannon-weaver多樣性指數(shù)、標(biāo)準(zhǔn)差、變異系數(shù)、平均值等;使用數(shù)據(jù)處理軟件IBM SPSS Statistics 26進(jìn)行主成分分析(PCA)及歐式遺傳距離(Euclidean distance)的計(jì)算;基于IBM SPSS Statistics 26計(jì)算得出的歐式遺傳距離矩陣,利用MEGA 7采用非加權(quán)組平均法(UPGMA)對(duì)供試材料進(jìn)行聚類分析。其中變異系數(shù)CV=SD/MN×100%,SD為標(biāo)準(zhǔn)差,MN為平均值;Shannon-weaver指數(shù)H′=?Σln,其中為該性狀第個(gè)描述值出現(xiàn)的頻率;數(shù)值型性狀經(jīng)質(zhì)量化處理后計(jì)算Shannon-weaver多樣性指數(shù),參照郭榮華等[15]的數(shù)值型性狀質(zhì)量化處理方法將每一性狀的觀測(cè)值賦予不同的等級(jí)。使用SPSS Statistics 26進(jìn)行數(shù)據(jù)處理,使用TB tools繪制熱圖(Heat map)。

    表1 供試材料的名稱及原產(chǎn)地

    表2 新梢和葉片的質(zhì)量性狀描述代碼及數(shù)量性狀統(tǒng)計(jì)方法

    2 結(jié)果與分析

    2.1 茶樹葉片表型性狀的遺傳多樣性分析

    2.1.1 茶樹葉片形態(tài)學(xué)性狀的遺傳多樣性分析

    基于供試茶樹表型性狀觀測(cè)結(jié)果,對(duì)描述型性狀分布頻率及遺傳多樣性指數(shù)進(jìn)行統(tǒng)計(jì),結(jié)果如表3所示。12個(gè)描述型性狀中遺傳多樣性最小的是葉緣(0.81),最大的是芽葉色澤(1.27)。根據(jù)多樣性指數(shù)從高到低排序,依次是芽葉色澤>葉片著生狀態(tài)>葉形>葉色>葉質(zhì)>葉身>葉齒>芽葉茸毛>葉面>葉尖>葉片大小>葉緣。供試群體的芽葉色澤以淡綠色為主,占比42.44%;大多數(shù)植株芽葉茸毛數(shù)量表現(xiàn)為中等,占比55.95%;成熟葉片著生狀態(tài)以水平著生為主,占44.69%,其次是稍上斜著生,占27.97%,下垂著生占比最少,為8.68%;葉色大多表現(xiàn)為綠色及深綠色,分別占比58.84%和23.15%;葉面主要表現(xiàn)為平和微隆起,分別占比47.59%和44.70%;葉身以平為主,占比54.02%;葉緣微波占比最高,為63.67%,平其次,占30.87%;葉尖主要表現(xiàn)為漸尖,占比67.85%;葉齒銳度多表現(xiàn)為中等銳度,占54.98%;葉質(zhì)中等硬度的占比最高,為51.77%;葉形以橢圓形居多,占52.09%,其次是長(zhǎng)橢圓形,占28.62%;葉片大小以中、小葉為主,其中小葉占比62.06%,中葉占比31.51%。

    對(duì)供試茶樹數(shù)值型性狀的遺傳多樣性進(jìn)行統(tǒng)計(jì),結(jié)果見表4。從整體來看,供試茶樹的數(shù)量性狀表現(xiàn)出豐富的多樣性:供試材料葉長(zhǎng)范圍為4.00~18.00?cm,平均為8.12?cm;葉寬的變化范圍為1.70~8.60?cm,平均為3.42?cm;長(zhǎng)寬比最小值為0.93,最大值為4.00,平均值為2.41;葉面積指數(shù)從5.60?cm2到73.71?cm2不等,平均值為20.48?cm2。供試植株各數(shù)量性狀遺傳多樣性指數(shù)范圍為1.77~1.94,皆高于質(zhì)量性狀。變異系數(shù)以葉面積(53.33%)最大,其次是葉寬(26.71%),再次為葉長(zhǎng)(26.12%),葉片長(zhǎng)寬比的變異系數(shù)最小,為17.09%。

    2.1.2 茶樹葉片形態(tài)學(xué)性狀的主成分分析

    對(duì)供試茶樹的表型性狀進(jìn)行主成分分析,得到各主成分特征值變化折線圖(圖1)。其中橫坐標(biāo)顯示的是16個(gè)主成分,縱坐標(biāo)顯示的是各主成分對(duì)應(yīng)的特征值,由圖可知前6個(gè)主成分都包含較大的特征值(特征值>1)、折線較為陡峭,第7個(gè)主成分之后折線趨于平緩,說明前6個(gè)主成分包含的信息量大。對(duì)各主成分的特征值、貢獻(xiàn)率、累計(jì)貢獻(xiàn)率進(jìn)行統(tǒng)計(jì)(表5),結(jié)果表明前6個(gè)主成分的累計(jì)貢獻(xiàn)率達(dá)71.78%,能較好地概括16個(gè)表型性狀的大部分信息。其中第一主成分貢獻(xiàn)率為25.49%,主要受葉長(zhǎng)、葉寬、葉面積和葉片大小的影響,載荷值分別為0.87、0.93、0.96和0.91,表現(xiàn)為高度正相關(guān),其中占主導(dǎo)地位的是葉面積;第二主成分主要受到長(zhǎng)寬比和葉形的影響,載荷值為0.83和0.88,貢獻(xiàn)率為13.16%,其中占主導(dǎo)地位的是葉形;第三主成分貢獻(xiàn)率為10.21%,與葉色、葉面、葉身、葉緣、葉質(zhì)呈正相關(guān),載荷值分別為0.47、0.55、0.52、0.44、0.48,與葉齒呈負(fù)相關(guān),載荷值為–0.43,其中占主導(dǎo)地位的是葉面及葉身;第四主成分與芽葉茸毛高度正相關(guān),載荷值為0.79,與葉緣高度負(fù)相關(guān),載荷值為–0.60,貢獻(xiàn)率為9.04%;第五主成分主要受葉面的影響,載荷值為0.60,貢獻(xiàn)率為7.05%;第六主成分主要受葉尖性狀的影響,載荷值為0.51,貢獻(xiàn)率為6.83%。綜合來看,葉面積、葉形、葉面、芽葉茸毛、葉尖、葉身這6個(gè)表型性狀可以反映出16個(gè)表型性狀71.78%的信息,是造成供試茶樹群體表型差異的主要因素。

    表3 描述型表型性狀的分布頻率和多樣性指數(shù)

    2.1.3 茶樹葉片形態(tài)學(xué)性狀的種質(zhì)聚類分析

    根據(jù)供試材料的16個(gè)表型性狀信息,采用非加權(quán)平均法(UPMGA)基于歐式遺傳距離對(duì)32份茶樹群體種質(zhì)進(jìn)行聚類分析(圖2)。聚類分析結(jié)果表明,供試茶樹種質(zhì)被劃分為4個(gè)類群。類群Ⅰ包括2個(gè)亞群,亞群a包含龍井種等9個(gè)群體種,亞群b包括石阡苔茶等6個(gè)群體種。其主要特征是葉片以中小葉為主、葉質(zhì)中等、葉緣微波、葉身內(nèi)折、芽葉茸毛少、芽葉色澤綠。類群Ⅱ包括黃山種等8個(gè)群體種,其主要特征是葉片以中小葉為主、芽葉色澤黃綠及淡綠、芽葉茸毛多、葉緣平、葉形橢圓形。類群Ⅲ包括黎家坪群體種等7個(gè)群體種,其主要特征是葉片以小葉為主、芽葉茸毛中等、葉片著生狀態(tài)上斜及稍上斜、葉色綠色、葉質(zhì)中等偏硬。類群Ⅳ包括勐庫(kù)大葉茶和鳳慶大葉茶,主要特征是大葉種、葉質(zhì)柔軟、葉齒較鈍、葉面微隆起、葉色淡綠色、芽葉茸毛中等。

    2.2 茶樹種質(zhì)主要生化成分的遺傳多樣性分析

    2.2.1 主要生化成分含量遺傳多樣性分布

    春季樣品測(cè)定結(jié)果顯示(圖3),供試茶樹種質(zhì)的游離氨基酸含量范圍為34~49?mg·g-1(干重),其中氨基酸含量在37~40?mg·g-1的茶樹材料有15份(占比46.88%),游離氨基酸含量在>40~43?mg·g-1的茶樹材料有8份(占比25.00%);供試茶樹材料的咖啡堿含量在20~45?mg·g-1,其中12份茶樹材料咖啡堿含量在30~35?mg·g-1(占比37.50%),9份茶樹材料咖啡堿含量在>35~40?mg·g-1(占比28.13%);供試茶樹材料的茶多酚含量位于212~352?mg·g-1,其中19份茶樹材料的茶多酚含量在268~296?mg·g-1(占比59.38%)。

    表4 數(shù)值型表型性狀的變異情況和分布特征

    圖1 主成分特征值變化折線圖

    表5 基于16個(gè)表型性狀的茶樹群體種質(zhì)主成分分析

    2.2.2 主要生化成分含量遺傳多樣性分析

    對(duì)32份茶樹地方群體種資源的兩年生化成分含量進(jìn)行測(cè)定、統(tǒng)計(jì),供試材料表現(xiàn)出豐富的遺傳多樣性,結(jié)果見表6。從整體上看,春季,供試材料的游離氨基酸含量平均值為32.28?mg·g-1(干重),變異系數(shù)為5.78%,其中鳥王茶含量最低(27.12?mg·g-1),宜春群體種含量最高(36.48?mg·g-1);咖啡堿含量平均值為21.54?mg·g-1,變異系數(shù)為18.63%,變化范圍12.29~31.19?mg·g-1,含量最低的是鳥王茶,含量最高的是遂川群體種;茶多酚含量平均值為200.40?mg·g-1,變異系數(shù)為13.61%,含量最低的是宜春群體種(153.71?mg·g-1),含量最高的是江華苦茶(277.07?mg·g-1)。2020年夏季供試材料的游離氨基酸含量平均值為40.52?mg·g-1,變異系數(shù)為7.52%,竹山大黑葉含量最低(36.49?mg·g-1),十里香含量最高,為46.90?mg·g-1;咖啡堿含量平均值為33.81?mg·g-1,變異系數(shù)為14.55%,鳳慶大葉茶含量最低(22.70?mg·g-1),遂川群體種含量最高,為42.11?mg·g-1);茶多酚含量平均值為282.00?mg·g-1,變異系數(shù)為8.11%,變化范圍為214.18~349.93?mg·g-1,龍井種含量最低,江華苦茶含量最高。

    注:1. 石阡苔茶,2. 宜春群體種,3. 鳩坑種,4. 都勻毛尖,5. 竹山大黑葉,6. 江華苦茶,7. 城步峒茶,8. 坦洋菜茶,9. 上猶群體種,10. 龍井種,11. 婺源種,12. 南山白毛茶,13. 修水群體種,14. 永修群體種,15. 遂川群體種,16. 黎家坪群體種,17. 浮梁種,18. 馮家臺(tái)群體種,19. 宜昌大葉茶,20. 汝城白毛茶,21. 云臺(tái)山種,22. 牙己茶,23. 信陽(yáng)種,24. 黃山種,25. 十里香,26. 鳥王茶,27. 祁門種,28. 紫陽(yáng)種,29. 勐庫(kù)大葉茶,30. 羅坑群體種,31. 六堡茶,32. 鳳慶大葉茶

    綜合兩季生化成分含量的測(cè)定結(jié)果對(duì)含量特異的茶樹種質(zhì)材料進(jìn)行挖掘,結(jié)果表明,春季樣品中宜春群體種的茶多酚含量最低,夏季樣品中其茶多酚含量亦低于大部分群體種;江華苦茶在春、夏兩季中咖啡堿和茶多酚含量都較高;龍井種的茶多酚含量在春、夏兩季中皆相對(duì)較低;遂川群體種春、夏兩季咖啡堿含量皆為最高;黎家坪群體種多酚含量相對(duì)較高,且游離氨基酸含量在春季樣品中最高,在夏季樣品中亦相對(duì)較高;十里香游離氨基酸含量在夏季樣品中最高,春季亦較高;鳥王茶在兩季節(jié)的氨基酸和咖啡堿含量均極低;羅坑群體種在兩季節(jié)中的氨基酸含量均極高;鳳慶大葉茶在兩季中的咖啡堿含量均極低。

    2.2.3 茶樹種質(zhì)主要生化成分的聚類分析

    使用TB tools對(duì)夏季茶樹地方群體種資源的生化成分含量數(shù)據(jù)進(jìn)行歸一化處理后,基于歐式遺傳距離采用最長(zhǎng)距離法(Complete)聚類并繪制熱圖,結(jié)果如圖4所示,供試茶樹群體種被劃分為3個(gè)類群。類群Ⅰ包含遂川群體種等9個(gè)群體種,其主要特征是氨基酸含量中等偏低,咖啡堿含量較高。類群Ⅱ包含3個(gè)亞群,亞群a包括江華苦茶、城步峒茶、黎家坪群體種和祁門種,其主要特征是咖啡堿和茶多酚含量較高,氨基酸含量中等,其中江華苦茶具有極高的茶多酚含量;亞群b包括石阡苔茶等6個(gè)群體種,其主要特征是氨基酸含量較低,咖啡堿和茶多酚含量居中;亞群c包括鳥王茶、宜昌大葉茶、云臺(tái)山種和紫陽(yáng)種,其主要特征是氨基酸含量偏低,咖啡堿和茶多酚含量低。類群Ⅲ包含鳳慶大葉茶等9個(gè)群體種,除龍井種外氨基酸含量皆高于其他群體種,鳳慶大葉種的咖啡堿含量在供試群體中最低。

    從供試種質(zhì)資源來源省份的地理分布范圍來看(表7),來自同一轄區(qū)內(nèi)的茶樹種質(zhì)資源多聚為一類,江西省有5份供試茶樹群體種聚在一起,來自湖南、湖北、安徽的全部種質(zhì)聚在類群Ⅱ中,來自云南、廣西的全部種質(zhì)聚在類群Ⅲ中。貴州省和江西省茶樹種質(zhì)資源表現(xiàn)出豐富的多樣性,分別屬于兩種及以上類群:貴州?。á?、Ⅱ)、江西?。á?、Ⅱ、Ⅲ)。

    3 討論

    3.1 茶樹葉片表型性狀的遺傳多樣性分析

    表型性狀是區(qū)分不同群體種茶樹的重要依據(jù)[16],本研究對(duì)32份供試茶樹種質(zhì)資源的表型性狀進(jìn)行遺傳多樣性分析,16個(gè)表型性狀的Shannon-weaver遺傳多樣性指數(shù)變化范圍為0.81~1.94。其中12個(gè)描述型性狀的遺傳多樣性指數(shù)最小值為0.81,最大值為1.27,遠(yuǎn)高于前人關(guān)于鳩坑茶樹資源及貴州古茶樹芽葉質(zhì)量性狀調(diào)查所得的遺傳多樣性指數(shù)[17-18],說明供試材料表型性狀變異較為豐富。整體看,供試茶樹種質(zhì)的芽葉色澤以淡綠色為主,芽葉茸毛數(shù)量表現(xiàn)為中等,成熟葉片著生狀態(tài)以水平為主,葉色大多表現(xiàn)為綠色,葉面主要表現(xiàn)為平,葉身以平為主,葉緣大多表現(xiàn)為微波,葉尖主要表現(xiàn)為漸尖,葉齒銳度多表現(xiàn)為中等銳度,葉質(zhì)主要表現(xiàn)為中等硬度,葉形以橢圓形和長(zhǎng)橢圓形為主,葉片大小以中小葉為主,其中小葉占比較多,與周炎花[19]關(guān)于國(guó)內(nèi)茶樹品種芽葉表型的調(diào)查結(jié)果較為一致。

    圖3 夏季樣品主要生化成分含量分布圖

    表6 32份種質(zhì)資源主要生化成分含量

    不同種質(zhì)表型性狀主成分分析也可作為種質(zhì)評(píng)價(jià)的一種手段,目前已被廣泛用于楊梅、獼猴桃、柑橘、蘋果等種質(zhì)資源的評(píng)價(jià)工作中[20-25],通過計(jì)算不同種質(zhì)的綜合得分,可以實(shí)現(xiàn)對(duì)優(yōu)異種質(zhì)資源的快速便捷篩選。茶樹16個(gè)表型性狀的主成分分析及聚類分析結(jié)果表明,葉面積、葉形、葉面、芽葉茸毛、葉尖、葉身這6個(gè)表型性狀是造成供試茶樹群體表型差異的主要因素;供試茶樹除類群Ⅳ表現(xiàn)為大葉種外,其余3個(gè)類群皆表現(xiàn)為中、小葉種,其中類群Ⅲ葉片偏小;類群Ⅱ葉形以橢圓形為主,其余類群葉形多在橢圓、長(zhǎng)橢圓間變異;類群Ⅳ葉面多表現(xiàn)為隆起,類群Ⅰ中亞群b多表現(xiàn)為平;芽葉茸毛以類群Ⅰ中亞群b最多,類群Ⅰ中亞群a最少;類群Ⅳ尖多表現(xiàn)為急尖和漸尖,其他類群多表現(xiàn)為漸尖;類群Ⅲ葉身包含平和稍背卷兩大類型,類群Ⅰ和類群Ⅳ葉身多表現(xiàn)為內(nèi)折、平。主成分分析篩選得到的6個(gè)性狀基本上可以將聚類分析劃分的四大類型區(qū)分開,且親緣關(guān)系近、地域位置近的材料大多聚在一起。林子翔[26]對(duì)觀賞辣椒種質(zhì)進(jìn)行了主成分分析,挖掘了影響觀賞辣椒價(jià)值的主要農(nóng)藝性狀,并在此基礎(chǔ)上進(jìn)行聚類分析,39份供試材料被劃分為3個(gè)類型,各類型之間特征顯著、且差異較大,與本研究結(jié)果相似。本研究主成分分析中第四主成分與芽葉茸毛高度正相關(guān)(載荷值為0.79),與葉緣高度負(fù)相關(guān)(載荷值為–0.6),推測(cè)葉緣與芽葉茸毛存在相關(guān)性。

    注:1. 石阡苔茶,2. 宜春群體種,3. 鳩坑種,4. 都勻毛尖,5. 竹山大黑葉,6. 江華苦茶,7. 城步峒茶,8. 坦洋菜茶,9. 上猶群體種,10. 龍井種,11. 婺源種,12. 南山白毛茶,13. 修水群體種,14. 永修群體種,15. 遂川群體種,16. 黎家坪群體種,17. 浮梁種,18. 馮家臺(tái)群體種,19. 宜昌大葉茶,20. 汝城白毛茶,21. 云臺(tái)山種,22. 牙己茶,23:信陽(yáng)種,24. 黃山種,25. 十里香,26. 鳥王茶,27. 祁門種,28. 紫陽(yáng)種,29. 勐庫(kù)大葉茶,30. 羅坑群體種,31. 六堡茶,32. 鳳慶大葉茶。

    表7 基于主要生化成分含量的32份茶樹群體種資源聚類族群

    3.2 茶樹種質(zhì)主要生化成分的遺傳多樣性分析

    兩年供試樣品生化成分測(cè)定結(jié)果表明,供試茶樹種質(zhì)資源的游離氨基酸、茶多酚、咖啡堿含量總體上符合正態(tài)分布,且表現(xiàn)出豐富的遺傳多樣性,供試樣品咖啡堿含量的變異系數(shù)(17.95%和14.55%)最大,其次為茶多酚含量的變異系數(shù)(13.61%和8.11%),再次為游離氨基酸的變異系數(shù)(5.62%和7.52%),說明供試樣品中資源類型豐富,可為選育生化成分含量特異(如低咖啡堿含量、低茶多酚含量等)的優(yōu)良茶樹單株提供參考。在立地條件和栽培管理措施基本一致的情況下,供試材料咖啡堿含量的差異較大,這與寧功偉等[27]、董方等[28]、丁帥濤等[29]的研究結(jié)果相似,推測(cè)原因可能是不同茶樹品種間遺傳背景差異所導(dǎo)致??Х葔A合成通路中-甲基轉(zhuǎn)移酶基因()[30]、咖啡堿合成酶基因()[31]等關(guān)鍵酶基因及[32]、[33]等轉(zhuǎn)錄因子的表達(dá)水平均會(huì)影響茶類植物中咖啡堿的含量。Jin等[34]研究表明,低咖啡堿茶樹品種中存在兩種控制咖啡堿生物合成的分子機(jī)制,即低轉(zhuǎn)錄水平的等位基因或其編碼的僅具有可可堿合成酶活性的蛋白質(zhì)。除此之外,植物體內(nèi)的生物堿與其抗逆性有關(guān),具有抵御病害、應(yīng)對(duì)環(huán)境脅迫等生理作用,咖啡堿作為茶樹體內(nèi)最主要的生物堿,其含量的差異與不同品種間抗逆性的聯(lián)系有待進(jìn)一步研究。綜合測(cè)定兩年的生化含量,結(jié)果發(fā)現(xiàn),生化成分含量特異的地方群體種資源有9份:宜春群體種(低茶多酚含量)、江華苦茶(高咖啡堿和茶多酚含量)、龍井種(低茶多酚含量)、遂川群體種(高咖啡堿含量)、黎家坪群體種(高茶多酚含量、高氨基酸含量)、十里香(高氨基酸含量)、鳥王茶(低咖啡堿含量)、羅坑群體種(高氨基酸含量)、鳳慶大葉茶(低咖啡堿含量),其中鳳慶大葉茶咖啡堿含量極低,且兩年測(cè)定結(jié)果一致,可以作為篩選低咖啡堿含量種質(zhì)或研究茶樹品種苦味改良的育種材料[35]。目前針對(duì)特定人群的低咖啡堿產(chǎn)品開發(fā)取得了快速進(jìn)展,但大多數(shù)脫咖啡堿工藝會(huì)影響產(chǎn)品的風(fēng)味且增加了額外的加工成本,因此選育低咖啡堿茶樹品種具有廣闊發(fā)展前景[36]。蔣會(huì)兵等[37]通過對(duì)云南省內(nèi)51份茶樹地方品種進(jìn)行生化成分測(cè)定,篩選出包括鳳慶大葉茶在內(nèi)的4份具有較低咖啡堿含量的茶樹地方群體種資源,以及昆明十里香在內(nèi)的適制綠茶的群體種,與本研究結(jié)果一致。對(duì)32份供試茶樹群體種的生化成分含量測(cè)定結(jié)果進(jìn)行聚類分析,供試茶樹群體被劃分為3類,分別包含不同類型的種質(zhì)資源:類群Ⅰ的主要特征是咖啡堿含量較高;類群Ⅱ中亞群c主要表現(xiàn)為較低的咖啡堿和茶多酚含量;類群Ⅲ包含具有較高氨基酸含量、較低咖啡堿含量的種質(zhì)資源。茶樹群體種變異類型多但自然變異具有不確定性,對(duì)各群體種基于生化成分含量劃分類群,了解群體種的共性,可為育種材料及相關(guān)代謝研究的材料選擇、核心種質(zhì)的篩選提供科學(xué)依據(jù)。

    從各類型茶樹種質(zhì)資源地理分布范圍來看,不同類型的茶樹種質(zhì)資源分布具有明顯的地域差異:類群Ⅲ除龍井種外主要分布在廣東、廣西、云南等茶區(qū),這些茶區(qū)位于中國(guó)西南及南方;而其他省份種質(zhì)資源除上猶群體種外皆屬于類群Ⅰ、類群Ⅱ。其中來自貴州省、江西省的供試茶樹群體種分布于多個(gè)類群,說明這些地區(qū)包含的茶樹群體種質(zhì)資源豐富。姚明哲[38]、陳勛等[39]基于SSR分子標(biāo)記技術(shù)對(duì)茶樹種質(zhì)進(jìn)行遺傳多樣性進(jìn)行了研究,結(jié)果亦表明,江西省、貴州省的茶樹種質(zhì)在中國(guó)各省份茶樹種質(zhì)資源中表現(xiàn)出中度偏高的遺傳多樣性水平。另外,隨著無性系品種的推廣、引進(jìn),大多栽植地方群體種的茶園被改植換代,群體種質(zhì)遭到破壞,若不及時(shí)加以保護(hù),會(huì)導(dǎo)致物種基因多樣性降低、大量種質(zhì)類型丟失[40]。中國(guó)西南地區(qū)作為茶樹起源的中心地帶,江西省作為茶樹傳播、演化的中間地帶,含有豐富的茶樹群體種質(zhì)資源[41],應(yīng)重點(diǎn)開展這些地區(qū)茶樹地方群體種表型性狀、生化性狀、遺傳背景、群體結(jié)構(gòu)研究工作,推進(jìn)茶樹群體種質(zhì)的發(fā)掘與保護(hù),加快核心種質(zhì)庫(kù)的建立,為茶樹遺傳育種提供更多寶貴的種質(zhì)材料[42]。

    3.3 基于表型與生化成分兩種聚類結(jié)果的比較

    利用聚類分析對(duì)茶樹品種資源進(jìn)行分類能夠快捷、準(zhǔn)確地得出資源間的相似性并展現(xiàn)出品種的綜合性狀,被靈活運(yùn)用于育種工作的親本選配環(huán)節(jié)。本研究基于表型性狀與主要生化成分分別對(duì)32份供試材料進(jìn)行聚類分析,不同分析標(biāo)準(zhǔn)聚類結(jié)果不完全一致,基于表型性狀的分析將供試材料分為4類,而基于主要生化成分的分析將供試材料分為3類?;诒硇偷木垲惤Y(jié)果中聚在類群Ⅰ的石阡苔茶、竹山大黑葉、上猶群體種、龍井種、南山白毛茶、羅坑群體種、六堡茶在基于生化成分的聚類結(jié)果中分離至類群Ⅱ與類群Ⅲ;基于表型的聚類結(jié)果中聚在類群Ⅱ的信陽(yáng)種、牙己茶、十里香在基于生化成分的聚類結(jié)果中分離至類群Ⅰ與類群Ⅲ;基于表型的聚類結(jié)果中聚在類群Ⅲ的遂川群體種、浮粱種在基于生化成分的聚類結(jié)果中分離至類群Ⅰ。兩種聚類結(jié)果也有相同的地方,例如基于表型性狀的聚類分析中聚為類群Ⅳ的勐庫(kù)大葉茶與鳳慶大葉茶在基于生化成分的分析結(jié)果中依舊聚在同一類中,且距離相近;黎家坪群體種、馮家臺(tái)群體種、汝城白毛茶、永修群體種、宜昌大葉茶始終在同一類群中;祁門種、黃山種、鳥王茶、云臺(tái)山種、紫陽(yáng)種總是聚在同一類群。

    綜合分析兩種聚類結(jié)果表明,表型性狀與生化成分間具有一定相關(guān)性,說明供試材料遺傳多樣性豐富、遺傳背景復(fù)雜。表型多樣性是遺傳多樣性與環(huán)境多樣性共同作用的結(jié)果[43],有研究表明,茶樹表型性狀易受環(huán)境條件影響,尤其是數(shù)量性狀變異較大,表現(xiàn)出豐富的遺傳多樣性[44]。生化成分含量水平的變異在長(zhǎng)期的自然選擇與人為馴化下發(fā)生了明顯的變異[45]。兩種聚類結(jié)果的差異表明,依據(jù)單一表型或生化成分的聚類分析只能作為參考,在實(shí)際應(yīng)用中還需結(jié)合更深層且靈敏度更高的分子標(biāo)記方法,綜合考慮表型、生化成分與DNA分子間的遺傳多樣性,才能高效精準(zhǔn)地指導(dǎo)茶樹育種工作。

    [1] 郭燕, 張樹航, 李穎, 等. 燕山板栗種質(zhì)資源葉片表型性狀多樣性研究[J]. 園藝學(xué)報(bào), 2022, 49(8): 1673-1688. Guo Y, Zhang S H, Li Y, et al. Diversity analysis of leaves phenotypic traits of Yanshan Chestnut [J]. Acta Horticulturae Sinica, 2022, 49(8): 1673-1688.

    [2] 姬廣磊. 黃鮟鱇()形態(tài)特征及其遺傳學(xué)研究[D]. 青島: 中國(guó)海洋大學(xué), 2007. Ji G L. Study on the morphological characters and genetics of[D]. Qingdao: Ocean University of China, 2007.

    [3] 陳琪予, 馬建強(qiáng), 陳杰丹, 等. 利用圖像特征分析茶樹成熟葉表型的遺傳多樣性[J]. 茶葉科學(xué), 2022, 42(5): 649-660. Chen Q Y, Ma J Q, Chen J D, et al. Genetic diversity of mature leaves of tea germplasms based on image features [J]. Journal of Tea Science, 2022, 42(5): 649-660.

    [4] Naegele R P, Mitchell J, Hausbeck M K. Genetic diversity, population structure, and heritability of fruit traits in[J]. PLos One, 2016, 11(7): e0156969. doi: 10.1371/journal.pone.0156969

    [5] 吳國(guó)林. 基于ISSR和SSR標(biāo)記的鄱陽(yáng)湖流域野生菰種質(zhì)資源遺傳多樣性研究[D]. 南昌: 江西農(nóng)業(yè)大學(xué), 2014. Wu G L. Studies on genetic diversity of germplasm resources of(Griseb.) Trucz. based on ISSR and SSR markers [D]. Nanchang: Jiangxi Agricultural University, 2014.

    [6] 劉彤. 柳州九萬(wàn)山野生茶樹種質(zhì)資源研究[D]. 長(zhǎng)沙: 湖南農(nóng)業(yè)大學(xué), 2016. Liu T. The research of Liuzhou Jiuwan Mountain wild tea germplasm resources [D]. Changsha: Hunan Agricultural University, 2016.

    [7] 莫嵐, 黃亞輝, 趙文芳, 等. 金秀野生茶樹群體的生化特征[J]. 福建農(nóng)林大學(xué)學(xué)報(bào)(自然科學(xué)版), 2018, 47(4): 409-414. Mo L, Huang Y H, Zhao W F, et al. Analysis of biochemical components of wild tea colonies in Jinxiu [J]. Journal of Fujian Agriculture and Forestry University (Natural Science Edition), 2018, 47(4): 409-414.

    [8] 晏嫦妤, 李家賢, 黃華林, 等. 廣東茶樹群體品種品質(zhì)成分差異初步研究[J]. 廣東農(nóng)業(yè)科學(xué), 2009, 59(12): 39-41. Yan C Y, Li J X, Huang H L, et al. Study on differences in quality composition of Guangdong tea germplasms [J]. Guangdong Agricultural Sciences, 2009, 59(12): 39-41.

    [9] 黃亞輝, 羅榮康, 晏嫦妤, 等. 連南茶樹種質(zhì)資源葉片表型性狀及生化成分的研究[J]. 廣東茶業(yè), 2017, 45(6): 11-13. Huang Y H, Luo R K, Yan C Y, et al. Study on leaf phenotypic traits and biochemical components of tea germplasm resources in Liannan [J]. Guangdong Tea Industry, 2017, 45(6): 11-13.

    [10] 王小萍, 唐曉波, 王迎春, 等. 52份茶樹資源生化組分的表型多樣性分析[J]. 茶葉科學(xué), 2012, 32(2):129-134. Wang X P, Tang X B, Wang Y C, et al. Phenotypic diversity analysis of biochemical composition in 52 tea germplasms [J]. Journal of Tea Science, 2012, 32(2): 129-134.

    [11] 甘春萍, 鄧慧群, 陳佳, 等. 廣西地方群體品種主要生化成分分析研究[J]. 農(nóng)業(yè)與技術(shù), 2012, 32(4): 95-99. Gan C P, Deng H Q, Chen J, et al. Analysis of main biochemical components of Guangxi local population varieties [J]. Agriculture and Technology, 2012, 32(4): 95-99.

    [12] 萬(wàn)明軒. 菠菜表型三維點(diǎn)云數(shù)據(jù)分析與形態(tài)品質(zhì)指標(biāo)遺傳多樣性研究[D]. 上海: 上海師范大學(xué), 2020. Wan M X. Study on analysis of three-dimensional point cloud data of spinach phenotype and genetic diversity of morphological quality index [D]. Shanghai: Shanghai Normal University, 2020.

    [13] 張悅, 索玉靜, 孫鵬, 等. 柿種質(zhì)資源果實(shí)形態(tài)多樣性分析[J]. 園藝學(xué)報(bào), 2022, 49(7): 1473-1490. Zhang Y, Suo Y J, Sun P, et al. Analysis on fruit morphological diversity of persimmon germplasm resources [J]. Acta Horticulturae Sinica, 2022, 49(7): 1473-1490.

    [14] 陳亮, 楊亞軍. 茶樹種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)[M]. 北京: 中國(guó)農(nóng)業(yè)出版社, 2005. Chen L, Yang Y J. Description specification and data standard of tea germplasm resources [M]. Beijing: China Agriculture Press, 2005.

    [15] 郭榮華, 劉志齋, 蔡一林, 等. 兩個(gè)歷史上重要的玉米地方品種群體的表型多樣性評(píng)價(jià)[J]. 玉米科學(xué), 2009, 17(2): 24-28. Guo R H, Liu Z Z, Cai Y L, et al. Evaluation of phenotypic diversity for two maize landrace population shaving historical importance in China [J]. Journal of Maize Sciences, 2009, 17(2): 24-28.

    [16] 丁洲, 李燁昕, 袁藝, 等. 安徽茶區(qū)優(yōu)良群體種的表型性狀和遺傳多樣性分析[J]. 茶葉科學(xué), 2018, 38(2):155-161. Ding Z, Li Y X, Yuan Y, et al. Phenotypic traits and genetic diversity of elite tea population in Anhui (China) [J]. Journal of Tea Science, 2018, 38(2): 155-161.

    [17] 黃海濤, 王風(fēng)雷, 王賢波, 等. 鳩坑茶樹資源芽葉表型性狀的多樣性分析[J]. 浙江農(nóng)業(yè)科學(xué), 2013, 64(3): 251-254. Huang H T, Wang F L, Wang X B, et al. Diversity analysis of bud and leaf phenotypic traits of Jiukeng tea resources [J]. Journal of Zhejiang Agricultural Sciences, 2013, 64(3): 251-254.

    [18] 劉青. 貴州三都縣古茶樹遺傳多樣性研究[D]. 貴陽(yáng): 貴州大學(xué), 2020. Liu Q. Genetic diversity of ancient tea trees in Sandu County, Guizhou Province [D]. Guiyang: Guizhou University, 2020.

    [19] 周炎花. 基于葉片形態(tài)和EST-SSR茶樹遺傳多樣性與遺傳演化研究[D]. 福州: 福建農(nóng)林大學(xué), 2010. Zhou Y H. Study of genetic diversity and genetic evolution of tea germplasms based on leaf morphology and EST-SSR markers [D]. Fuzhou: Fujian Agriculture and Forestry University, 2010.

    [20] 胡波, 李國(guó)元, 謝志兵, 等. 主成分分析法在冬棗果實(shí)貯藏特性分析中的應(yīng)用[J]. 食品科技, 2014, 39(2): 44-47. Hu B, Li G Y, Xie Z B, et al. Application of principle component analysis on winter jujube storage characteristic [J]. Food Science and Technology, 2014, 39(2): 44-47.

    [21] 楊玲, 張彩霞, 康國(guó)棟, 等. ‘華紅’蘋果果肉的流變特性及其主成分分析[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2015, 48(12): 2417-2427. Yang L, Zhang C X, Kang G D, et al. Rheologic properties of ‘Huahong’ apple pulp and their principal component analysis [J]. Scientia Agricultura Sinica, 2015, 48(12): 2417-2427.

    [22] 郎彬彬, 黃春輝, 朱博, 等. 基于果實(shí)相關(guān)性狀的江西野生毛花獼猴桃初級(jí)核心種質(zhì)的構(gòu)建方法研究[J]. 果樹學(xué)報(bào), 2016, 33(7): 794-803. Lang B B, Huang C H, Zhu B, et al. Study on the method of constructing a primary core collection of Jiangxi wildbased on fruit traits [J]. Journal of Fruit Science, 2016, 33(7): 794-803.

    [23] 李偉, 郜海燕, 陳杭君, 等. 基于主成分分析的不同品種楊梅果實(shí)綜合品質(zhì)評(píng)價(jià)[J]. 中國(guó)食品學(xué)報(bào), 2017, 17(6): 161-171. Li W, Gao H Y, Chen H J, et al. Evaluation of comprehensive quality of different varieties of bayberry based on principal components analysis [J]. Journal of Chinese Institute of Food Science and Technology, 2017, 17(6): 161-171.

    [24] 李勛蘭, 洪林, 王武, 等. 晚熟雜柑新品種果實(shí)品質(zhì)綜合評(píng)價(jià)[J]. 果樹學(xué)報(bào), 2018, 35(2): 195-203. Li X L, Hong L, Wang W, et al. Comprehensive evaluation of fruit quality of new late-maturing mandarin cultivars [J]. Journal of Fruit Science, 2018, 35(2): 195-203.

    [25] 呂正鑫, 賀艷群, 賈東峰, 等. 獼猴桃種質(zhì)資源表型性狀遺傳多樣性分析[J]. 園藝學(xué)報(bào), 2022, 49(7): 1571-1581. Lü Z X, He Y Q, Jia D F, et al. Genetic diversity analysis of phenotypic traits for kiwifruit germplasm resource [J]. Acta Horticulturae Sinica, 2022, 49(7): 1571-1581.

    [26] 林子翔. 觀賞辣椒種質(zhì)資源園藝性狀的多元統(tǒng)計(jì)分析與種質(zhì)創(chuàng)新研究[D]. 杭州: 浙江農(nóng)林大學(xué), 2020. Lin Z X. Multivariate statistical analysis and germplasm innovation of ornamental pepper germplasm resources [D]. Hangzhou: Zhejiang A&F University, 2020.

    [27] 寧功偉, 楊盛美, 段志芬, 等. 云南野生茶樹厚軸茶種質(zhì)資源化學(xué)成分多樣性分析[J/OL]. 分子植物育種, 2023: 1-20[2023-05-06]. https://kns.cnki.net/kcms/detail//46. 1068. S.20230118.1701.005.html. Ning G W, Yang S M, Duan Z F, et al. Diversity analysis of chemical components of wild tea germplasm resources () in Yunnan [J/OL]. Molecular Plant Breeding, 2023: 1-20[2023-05-06]. https://kns.cnki. net/ kcms/detail//46.1068.S.20230118.1701.005.html.

    [28] 董方, 李小飛, 沈思言, 等. 江西茶樹資源的遺傳多樣性分析及優(yōu)異種質(zhì)篩選[J]. 江西農(nóng)業(yè)大學(xué)學(xué)報(bào), 2022, 44(6): 1466-1477. Dong F, Li X F, Shen S Y, et al. Genetic diversity analysis and screening of excellent germplasm of tea plant resources in Jiangxi [J]. Acta Agriculturae Universitatis Jiangxiensis, 2022, 44(6): 1466-1477.

    [29] 丁帥濤, 程曉梅, 張亞, 等. 基于表型性狀和生化成分的陜西茶樹種質(zhì)資源遺傳多樣性研究[J]. 西北農(nóng)業(yè)學(xué)報(bào), 2019, 28(4): 607-619. Ding S T, Cheng X M, Zhang Y, et al. Genetic diversity in phenotypic traits and biochemical components of tea resources in Shanxi [J]. Acta Agriculturae Boreali-occidentalis Sinica, 2019, 28(4): 607-619.

    [30] Zhang Y R, Fu J M, Zhou Q Y, et al. Metabolite profiling and transcriptome analysis revealed the conserved transcriptional regulation mechanism of caffeine biosynthesis in tea and coffee germplasm plants [J]. Journal of Agricultural and Food Chemistry, 2022, 70(10): 3239-3251.

    [31] Li G D, Li Y, Yao X Z, et al. Establishment of a virus-induced gene-silencing (VIGS) system in tea plant and its use in the functional analysis of[J]. International Journal of Molecular Sciences, 2023, 24(1): 392. doi: 10.3390/ijms24010392.

    [32] Li P H, Ye Z L, Fu J M, et al. CsMYB184 regulates caffeine biosynthesis in tea plants [J]. Plant Biotechnology Journal, 2022, 20(6): 1012-1014.

    [33] Ma W H, Kang X, Liu P, et al. The NAC-like transcription factor CsNAC7 positively regulates the caffeine biosynthesis-related genein[J]. Horticulture Research, 2022, 9: uhab046. doi: 10.1093/hr/uhab046.

    [34] Jin J Q, Yao M Z, Ma C L, et al. Natural allelic variations of TCS1 play a crucial role in caffeine biosynthesis of tea plant and its related species [J]. Plant Physiology and Biochemistry, 2016, 100(12): 18-26.

    [35] 蔡志翔, 嚴(yán)娟, 宿子文, 等. 不同類型桃種質(zhì)資源主要酚類物質(zhì)含量評(píng)價(jià)[J]. 園藝學(xué)報(bào), 2022, 49(5): 1008-1022. Cai Z X, Yan J, Su Z W, et al. Genetic diversity analysis of phenotypic traits for kiwifruit germplasm resources [J]. Acta Horticulturae Sinica, 2022, 49(5): 1008-1022.

    [36] Zhu Q F, Liu L J, Lu X F, et al. The biosynthesis of EGCG, theanine and caffeine in response to temperature is mediated by hormone signal transduction factors in tea plant (L.) [J]. Frontiers in Plant Science, 2023, 14. doi: 10.3389/fpls.2023.1149182.

    [37] 蔣會(huì)兵, 田易萍, 陳林波, 等. 云南茶樹地方品種農(nóng)藝性狀與品質(zhì)性狀遺傳多樣性分析[J]. 植物遺傳資源學(xué)報(bào), 2013, 14(4): 634-640. Jiang H B, Tian Y P, Chen L B, et al. Diversity of tea landraces based on agronomic and qualitytraits in Yunnan Province [J]. Journal of Plant Genetic Resources, 2013, 14(4): 634-640.

    [38] 姚明哲. 利用ISSR和EST-SSR標(biāo)記研究中國(guó)茶樹資源的遺傳多樣性和遺傳結(jié)構(gòu)[D]. 杭州: 浙江大學(xué), 2009. Yao M Z. Studies on genetic diversity and structure of tea germplasm in China based on ISSR and EST-SSR markers [D]. Hangzhou: Zhejiang University, 2009.

    [39] 陳勛, 龔自明. 基于EST-SSR標(biāo)記的湖北茶樹種質(zhì)資源遺傳多樣性分析[J]. 分子植物育種, 2017, 15(5): 1831-1838. Chen X, Gong Z M. Analysis of genetic diversity with EST-SSR markers for tea germplasm in Hubei Province [J]. Molecular Plant Breeding, 2017, 15(5): 1831-1838.

    [40] 張曉麗, 郭輝, 王海崗, 等. 中國(guó)普通野生稻與栽培稻種SSR多樣性的比較分析[J]. 作物學(xué)報(bào), 2008, 105(4): 591-597. Zhang X L, Guo H, Wang H G, et al. Comparative assessment of SSR allelic diversity in wild and cultivated rice in China [J]. Acta Agronomica Sinica, 2008, 105(4): 591-597.

    [41] 劉振, 趙洋, 楊培迪, 等. 湖南省茶樹種質(zhì)資源現(xiàn)狀及研究進(jìn)展[J]. 茶葉通訊, 2011, 38(3): 7-10, 18. Liu Z, Zhao Y, Yang P D, et al. The status and research progress of Hunan tea germplasm [J]. Tea Communication, 2011, 38(3): 7-10, 18.

    [42] 馬建強(qiáng), 姚明哲, 陳亮. 茶樹種質(zhì)資源研究進(jìn)展[J]. 茶葉科學(xué), 2015, 35(1):11-16. Ma J Q, Yao M Z, Chen L. Research progress on germplasms of tea plant () [J]. Journal of Tea Science, 2015, 35(1):11-16.

    [43] 張見明. 26份武夷名叢茶樹種質(zhì)資源秋茶生化成分多樣性分析[J]. 武夷學(xué)院學(xué)報(bào), 2016, 35(9): 1-7. Zhang J M. Diversity analysis of biochemical components in autumn shoots of 26 Wuyi Mingcong tea plant germplasm resources [J]. Journal of Wuyi University, 2016, 35(9): 1-7.

    [44] 喬婷婷. 茶樹資源遺傳多樣性及其表型性狀關(guān)聯(lián)EST-SSR位點(diǎn)的初步鑒定[D]. 北京: 中國(guó)農(nóng)業(yè)科學(xué)院, 2010. Qiao T T. Genetic diversity of tea ((L.)O. Kuntze) and association analysis on phenotypic traits with EST-SSR markers [D]. Beijing: Chinese Academy of Agricultural Sciences, 2010.

    [45] 齊冬晴. 不同茶樹品種(品系)鮮葉香氣及其表型差異的研究[D]. 長(zhǎng)沙: 湖南農(nóng)業(yè)大學(xué), 2015. Qi D Q. Studies on the tea leaves aroma components and phenotypic differences of different tea varieties [D]. Changsha: Hunan Agricultural University, 2015.

    Diversity Analysis of Leaf Phenotype and Biochemical Components in Tea Local Population Resources

    TANG Lu, LI Changle, GE Yue, WANG Pu, ZHAO Hua, WANG Mingle, WANG Yu, GUO Fei*, NI Dejiang

    National Key Laboratory for Germplasm Innovation & Utilization of Horticultural Crops, College of Horticulture and Forestry Sciences, Huazhong Agricultural University, Wuhan 430070, China

    In order to protect tea germplasm resources and promote tea germplasm innovation, the phenotypic character and biochemical component diversity of 32 tea germplasm resources from 12 provinces were studied. The results showed that the genetic diversity index of 16 phenotypic character ranged from 0.81 to 1.94, and the phenotypic genetic diversity of the test population was high. The cumulative contribution rate of the six phenotypic charactersincluding leaf area, leaf shape, leaf upper surface, young shoot pubescence, leaf tip shape and leaf cross section were 71.78%, which was the main factor causing the phenotypic difference of the tested tea population. Based on the determination results of biochemical components in two years, the variation coefficient of caffeine content (17.95% and 14.55%) was the largest, followed by the variation coefficient of tea polyphenol content (13.61% and 8.11%).The variation coefficient of free amino acid (5.62% and 7.52%) was the lowest. Cluster analysis based on the results of biochemical component content determination was conducted to divide the tested tea tree populations into three groups, each including different types of germplasm. From the geographical distribution of the tested tea populations, the content of biochemical components showed regional differences.

    tea local population species, phenotypic character, biochemical components, diversity

    S571.1

    A

    1000-369X(2023)04-473-16

    2023-05-06

    2023-07-17

    農(nóng)業(yè)生物種質(zhì)資源挖掘與利用(2022YFD1200505)、臨滄市建設(shè)國(guó)家可持續(xù)發(fā)展議程創(chuàng)新示范區(qū)科技專項(xiàng)(202204AC100001-A01)

    唐璐,女,碩士研究生,主要從事茶樹種質(zhì)資源方面研究,Tanglu@webmail.hzau.edu.cn。*通信作者:guofei@mail.hzau.edu.cn

    猜你喜歡
    咖啡堿表型茶樹
    安徽農(nóng)業(yè)大學(xué)揭示茶樹咖啡堿合成調(diào)控機(jī)制
    山茶樹變身?yè)u錢樹
    19份茶樹資源鮮葉中咖啡堿合成酶的活性
    建蘭、寒蘭花表型分析
    大孔吸附樹脂富集速溶茶中咖啡堿的研究
    兩個(gè)推薦茶樹品種
    晚上喝茶頭遍要倒掉
    茶樹灣
    GABABR2基因遺傳變異與肥胖及代謝相關(guān)表型的關(guān)系
    慢性乙型肝炎患者HBV基因表型與血清學(xué)測(cè)定的臨床意義
    国产成人精品福利久久| 久久久久精品人妻al黑| 亚洲伊人色综图| 国产无遮挡羞羞视频在线观看| 国产亚洲一区二区精品| 热re99久久精品国产66热6| 99九九在线精品视频| 国产精品女同一区二区软件| 免费久久久久久久精品成人欧美视频 | 少妇人妻 视频| 久久久久人妻精品一区果冻| 韩国av在线不卡| 男人添女人高潮全过程视频| 日本-黄色视频高清免费观看| 欧美人与性动交α欧美精品济南到 | 精品国产乱码久久久久久小说| av一本久久久久| 亚洲精品国产av成人精品| 婷婷成人精品国产| 精品少妇黑人巨大在线播放| 欧美亚洲 丝袜 人妻 在线| 久久人人爽人人爽人人片va| 国产亚洲精品久久久com| 18禁观看日本| 免费看光身美女| 最新的欧美精品一区二区| 精品久久久久久电影网| 又大又黄又爽视频免费| 国产色婷婷99| 久久久精品免费免费高清| 国产成人午夜福利电影在线观看| 国产精品不卡视频一区二区| 国产一区二区激情短视频 | 精品人妻一区二区三区麻豆| 久久亚洲国产成人精品v| 国产一区二区三区综合在线观看 | 爱豆传媒免费全集在线观看| 日韩大片免费观看网站| 26uuu在线亚洲综合色| 久久99蜜桃精品久久| 婷婷色麻豆天堂久久| 22中文网久久字幕| 赤兔流量卡办理| 在线精品无人区一区二区三| 免费人妻精品一区二区三区视频| xxxhd国产人妻xxx| 啦啦啦啦在线视频资源| 性高湖久久久久久久久免费观看| 一区二区av电影网| 性色avwww在线观看| 麻豆乱淫一区二区| 啦啦啦啦在线视频资源| 丰满饥渴人妻一区二区三| 精品亚洲乱码少妇综合久久| 午夜视频国产福利| 一级毛片黄色毛片免费观看视频| 亚洲一级一片aⅴ在线观看| 久久午夜福利片| 久久精品国产亚洲av天美| 久久精品人人爽人人爽视色| 日韩大片免费观看网站| 97人妻天天添夜夜摸| 人妻 亚洲 视频| 亚洲av欧美aⅴ国产| 午夜久久久在线观看| 人妻少妇偷人精品九色| 一区二区日韩欧美中文字幕 | 最近中文字幕高清免费大全6| 美女视频免费永久观看网站| 精品酒店卫生间| 久久精品久久久久久噜噜老黄| 日本爱情动作片www.在线观看| 欧美亚洲日本最大视频资源| 亚洲av电影在线进入| 国产麻豆69| 毛片一级片免费看久久久久| 精品一区二区三卡| 日本与韩国留学比较| 多毛熟女@视频| 久久影院123| 亚洲第一av免费看| 国产色婷婷99| 欧美激情极品国产一区二区三区 | 日韩精品有码人妻一区| 久久99蜜桃精品久久| 久久久久久伊人网av| 亚洲综合精品二区| 午夜视频国产福利| 26uuu在线亚洲综合色| 国产精品久久久久久av不卡| 成人综合一区亚洲| 汤姆久久久久久久影院中文字幕| 亚洲国产精品一区三区| 亚洲图色成人| 香蕉丝袜av| 中文精品一卡2卡3卡4更新| 久久久久精品性色| av视频免费观看在线观看| 国产黄色免费在线视频| 青青草视频在线视频观看| 欧美日韩亚洲高清精品| 高清黄色对白视频在线免费看| 啦啦啦在线观看免费高清www| 老司机亚洲免费影院| 亚洲精品久久成人aⅴ小说| 看免费成人av毛片| 91国产中文字幕| 中文精品一卡2卡3卡4更新| 美女内射精品一级片tv| 亚洲 欧美一区二区三区| 人妻系列 视频| 亚洲在久久综合| 国产精品一国产av| 成人影院久久| 欧美亚洲日本最大视频资源| 成年人免费黄色播放视频| 国产精品秋霞免费鲁丝片| 色视频在线一区二区三区| 狠狠精品人妻久久久久久综合| 午夜免费男女啪啪视频观看| 日韩不卡一区二区三区视频在线| 国产av国产精品国产| 熟妇人妻不卡中文字幕| 黄色配什么色好看| 色吧在线观看| 水蜜桃什么品种好| 午夜免费观看性视频| 最新的欧美精品一区二区| 少妇的逼好多水| 少妇 在线观看| 国产色婷婷99| av在线观看视频网站免费| 国产国拍精品亚洲av在线观看| 亚洲欧美中文字幕日韩二区| 午夜视频国产福利| 26uuu在线亚洲综合色| 免费大片18禁| 亚洲美女搞黄在线观看| 黄色配什么色好看| 国产精品三级大全| 我要看黄色一级片免费的| av黄色大香蕉| 欧美成人精品欧美一级黄| 建设人人有责人人尽责人人享有的| 国产欧美日韩一区二区三区在线| 精品99又大又爽又粗少妇毛片| 中文字幕亚洲精品专区| 精品亚洲乱码少妇综合久久| 国产欧美日韩一区二区三区在线| 久久久久网色| a级毛色黄片| h视频一区二区三区| 汤姆久久久久久久影院中文字幕| 视频中文字幕在线观看| 日日啪夜夜爽| 亚洲经典国产精华液单| 丰满乱子伦码专区| 中文字幕精品免费在线观看视频 | 国产成人精品无人区| 搡老乐熟女国产| 婷婷色综合www| 国产精品国产av在线观看| 最近最新中文字幕大全免费视频 | 成年av动漫网址| 久久av网站| 日本wwww免费看| 色5月婷婷丁香| 精品国产乱码久久久久久小说| 亚洲综合精品二区| 不卡视频在线观看欧美| 久久久久久伊人网av| 最近手机中文字幕大全| 亚洲,欧美,日韩| 国产黄色免费在线视频| 婷婷色综合www| 丁香六月天网| 丝袜人妻中文字幕| 黄色怎么调成土黄色| 高清在线视频一区二区三区| 性高湖久久久久久久久免费观看| 伦理电影免费视频| 欧美 亚洲 国产 日韩一| 九九爱精品视频在线观看| 寂寞人妻少妇视频99o| 国产男人的电影天堂91| 亚洲美女搞黄在线观看| 国产有黄有色有爽视频| 黄色怎么调成土黄色| 国产精品久久久久久久久免| 黑丝袜美女国产一区| 制服诱惑二区| 老司机亚洲免费影院| 欧美激情国产日韩精品一区| 免费在线观看黄色视频的| 热99国产精品久久久久久7| 亚洲第一av免费看| 又大又黄又爽视频免费| 99九九在线精品视频| 午夜免费男女啪啪视频观看| 汤姆久久久久久久影院中文字幕| 欧美 亚洲 国产 日韩一| 午夜福利乱码中文字幕| 久久久久精品人妻al黑| 两个人免费观看高清视频| 看免费成人av毛片| 日韩熟女老妇一区二区性免费视频| 久久久久久久久久人人人人人人| freevideosex欧美| 欧美bdsm另类| 国产精品欧美亚洲77777| 美女内射精品一级片tv| 黄色视频在线播放观看不卡| av黄色大香蕉| 亚洲在久久综合| 国产 一区精品| 最近手机中文字幕大全| 日韩在线高清观看一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 色婷婷久久久亚洲欧美| 久久国产精品大桥未久av| 久久久精品94久久精品| 婷婷色av中文字幕| 母亲3免费完整高清在线观看 | 亚洲国产av影院在线观看| 熟女人妻精品中文字幕| 伦精品一区二区三区| 国产欧美另类精品又又久久亚洲欧美| a级毛片黄视频| 女人久久www免费人成看片| 国产精品秋霞免费鲁丝片| 一级黄片播放器| 色5月婷婷丁香| 日韩中文字幕视频在线看片| 国产老妇伦熟女老妇高清| 免费看光身美女| 久久久久精品久久久久真实原创| 王馨瑶露胸无遮挡在线观看| 国产色婷婷99| 男女啪啪激烈高潮av片| 久久韩国三级中文字幕| 青春草视频在线免费观看| 最新中文字幕久久久久| 秋霞在线观看毛片| 日韩,欧美,国产一区二区三区| 欧美精品国产亚洲| 麻豆精品久久久久久蜜桃| 熟女电影av网| 色哟哟·www| 十八禁高潮呻吟视频| 久久久国产精品麻豆| 亚洲国产精品专区欧美| 欧美xxxx性猛交bbbb| 色网站视频免费| 97在线视频观看| 欧美亚洲日本最大视频资源| 一级片免费观看大全| 国产高清国产精品国产三级| 精品亚洲成国产av| 中文天堂在线官网| 春色校园在线视频观看| 狠狠精品人妻久久久久久综合| 亚洲激情五月婷婷啪啪| av不卡在线播放| 久热这里只有精品99| 亚洲国产av影院在线观看| 2021少妇久久久久久久久久久| 少妇高潮的动态图| 久久人人爽人人爽人人片va| 日日啪夜夜爽| 男女下面插进去视频免费观看 | 久久久精品免费免费高清| www.色视频.com| 午夜福利影视在线免费观看| 性高湖久久久久久久久免费观看| 一级毛片 在线播放| 另类精品久久| 久久精品熟女亚洲av麻豆精品| 亚洲精品久久久久久婷婷小说| av不卡在线播放| 欧美激情 高清一区二区三区| 国产一区亚洲一区在线观看| 国产精品一二三区在线看| 国产免费福利视频在线观看| 精品国产一区二区久久| 一二三四中文在线观看免费高清| 国内精品宾馆在线| 国产1区2区3区精品| 成年av动漫网址| 热99久久久久精品小说推荐| freevideosex欧美| 少妇人妻 视频| 青青草视频在线视频观看| 激情视频va一区二区三区| 又黄又爽又刺激的免费视频.| 色94色欧美一区二区| 久久影院123| 亚洲国产精品国产精品| 国产欧美日韩综合在线一区二区| 又黄又粗又硬又大视频| 国产精品一区二区在线观看99| av在线观看视频网站免费| 五月天丁香电影| 女人久久www免费人成看片| 午夜激情久久久久久久| 久久久久久久亚洲中文字幕| 日本与韩国留学比较| 最近的中文字幕免费完整| 亚洲精品一二三| 啦啦啦在线观看免费高清www| 久久久久久久久久成人| 午夜久久久在线观看| 又黄又爽又刺激的免费视频.| 久久久久精品久久久久真实原创| 欧美xxxx性猛交bbbb| 最近最新中文字幕大全免费视频 | 中文字幕制服av| 男女边摸边吃奶| 国内精品宾馆在线| 制服丝袜香蕉在线| 青春草亚洲视频在线观看| 亚洲国产看品久久| 国产成人免费观看mmmm| 激情视频va一区二区三区| 亚洲精品自拍成人| 国语对白做爰xxxⅹ性视频网站| 久久鲁丝午夜福利片| 91aial.com中文字幕在线观看| 亚洲成人手机| 欧美成人精品欧美一级黄| 色94色欧美一区二区| 国产 一区精品| av天堂久久9| 如何舔出高潮| 久久精品国产亚洲av涩爱| 精品一区二区三区四区五区乱码 | 香蕉丝袜av| 一本色道久久久久久精品综合| 秋霞在线观看毛片| kizo精华| 亚洲欧洲国产日韩| 啦啦啦视频在线资源免费观看| 亚洲国产欧美日韩在线播放| 熟女电影av网| 在线观看国产h片| 另类亚洲欧美激情| 哪个播放器可以免费观看大片| 国产69精品久久久久777片| 亚洲中文av在线| 欧美日韩视频高清一区二区三区二| 飞空精品影院首页| 久久免费观看电影| 韩国av在线不卡| 99视频精品全部免费 在线| av卡一久久| 99国产精品免费福利视频| 精品久久蜜臀av无| 极品少妇高潮喷水抽搐| 亚洲色图 男人天堂 中文字幕 | 天天躁夜夜躁狠狠躁躁| 观看av在线不卡| 九色亚洲精品在线播放| 一级毛片黄色毛片免费观看视频| 亚洲四区av| 大香蕉久久网| 久久97久久精品| 精品一品国产午夜福利视频| 日本vs欧美在线观看视频| 亚洲国产精品国产精品| 91aial.com中文字幕在线观看| 在线观看www视频免费| 性高湖久久久久久久久免费观看| 9191精品国产免费久久| 久久女婷五月综合色啪小说| 9色porny在线观看| 久久久久久久久久成人| 97在线视频观看| 赤兔流量卡办理| 国产黄频视频在线观看| 亚洲国产色片| 国产精品 国内视频| 好男人视频免费观看在线| videosex国产| 亚洲国产av影院在线观看| 韩国av在线不卡| 国产精品久久久久久久电影| 校园人妻丝袜中文字幕| 在线观看人妻少妇| 一级片'在线观看视频| 久久99热6这里只有精品| 国产极品天堂在线| 国产乱来视频区| 国产极品天堂在线| 久热久热在线精品观看| 日本wwww免费看| av黄色大香蕉| 日韩中字成人| 韩国精品一区二区三区 | 欧美 亚洲 国产 日韩一| 肉色欧美久久久久久久蜜桃| 一区二区三区四区激情视频| 大香蕉久久网| 精品国产一区二区久久| 青春草国产在线视频| 亚洲人成77777在线视频| 成人影院久久| 飞空精品影院首页| 久久人妻熟女aⅴ| 国产精品欧美亚洲77777| 又黄又粗又硬又大视频| 在线精品无人区一区二区三| 免费播放大片免费观看视频在线观看| 97人妻天天添夜夜摸| 国产伦理片在线播放av一区| 我的女老师完整版在线观看| 美女内射精品一级片tv| 热99久久久久精品小说推荐| 国产一区二区三区av在线| 我要看黄色一级片免费的| 久久这里有精品视频免费| 国产麻豆69| 欧美日韩av久久| videosex国产| 国产精品人妻久久久影院| 超碰97精品在线观看| 99视频精品全部免费 在线| 人人澡人人妻人| 久久久久久人妻| 日产精品乱码卡一卡2卡三| 啦啦啦中文免费视频观看日本| 久久狼人影院| 欧美日韩成人在线一区二区| 一本久久精品| 国产淫语在线视频| 精品人妻偷拍中文字幕| 黄片播放在线免费| 久久精品aⅴ一区二区三区四区 | 最近中文字幕高清免费大全6| 丝袜在线中文字幕| 久久久亚洲精品成人影院| 91久久精品国产一区二区三区| 精品亚洲成a人片在线观看| 成人亚洲精品一区在线观看| 视频中文字幕在线观看| 亚洲丝袜综合中文字幕| 两性夫妻黄色片 | 成年女人在线观看亚洲视频| 成人手机av| 天天躁夜夜躁狠狠躁躁| xxx大片免费视频| 欧美日韩亚洲高清精品| 一本色道久久久久久精品综合| 在线观看www视频免费| 久久久久久久久久久久大奶| a级毛色黄片| 国产精品 国内视频| 亚洲精品自拍成人| 亚洲精品国产av成人精品| 成人免费观看视频高清| 一本大道久久a久久精品| 久久久久久久久久久久大奶| 日韩欧美精品免费久久| 永久免费av网站大全| 建设人人有责人人尽责人人享有的| 成人毛片a级毛片在线播放| 少妇熟女欧美另类| 91国产中文字幕| 麻豆乱淫一区二区| 日韩不卡一区二区三区视频在线| 亚洲欧美色中文字幕在线| 国产国语露脸激情在线看| 国产高清三级在线| 国产亚洲最大av| 亚洲国产看品久久| 免费大片18禁| 男女无遮挡免费网站观看| 大码成人一级视频| 国产成人一区二区在线| av天堂久久9| 国产白丝娇喘喷水9色精品| 日韩免费高清中文字幕av| 久久久久久久久久久免费av| 男男h啪啪无遮挡| 免费大片黄手机在线观看| av一本久久久久| av视频免费观看在线观看| 啦啦啦中文免费视频观看日本| 国产精品久久久久久久电影| 国产精品国产三级国产专区5o| 日韩不卡一区二区三区视频在线| 免费人成在线观看视频色| 天堂8中文在线网| 精品国产一区二区三区久久久樱花| 亚洲精品日本国产第一区| 90打野战视频偷拍视频| videossex国产| 日韩人妻精品一区2区三区| 最近手机中文字幕大全| 伊人亚洲综合成人网| 飞空精品影院首页| 国产精品成人在线| 亚洲伊人色综图| 国产精品欧美亚洲77777| 久久精品熟女亚洲av麻豆精品| 日韩成人av中文字幕在线观看| 成年人免费黄色播放视频| 日韩av免费高清视频| 国产色婷婷99| 欧美bdsm另类| 激情视频va一区二区三区| 亚洲国产精品一区二区三区在线| 毛片一级片免费看久久久久| h视频一区二区三区| av片东京热男人的天堂| 少妇熟女欧美另类| 男女边摸边吃奶| 精品久久久精品久久久| 成人综合一区亚洲| 看非洲黑人一级黄片| 青青草视频在线视频观看| 九九爱精品视频在线观看| 亚洲精品自拍成人| 黄色视频在线播放观看不卡| 中国美白少妇内射xxxbb| 免费大片18禁| 欧美激情国产日韩精品一区| 国产亚洲欧美精品永久| 欧美97在线视频| 观看av在线不卡| 久久久久久人人人人人| 欧美精品人与动牲交sv欧美| 亚洲欧美中文字幕日韩二区| 日韩av在线免费看完整版不卡| 欧美精品亚洲一区二区| 国产精品国产三级专区第一集| 久久久久国产精品人妻一区二区| a级毛片在线看网站| 校园人妻丝袜中文字幕| 国产日韩欧美在线精品| xxx大片免费视频| 国产爽快片一区二区三区| 黄片无遮挡物在线观看| av不卡在线播放| 亚洲欧美一区二区三区黑人 | 国产精品 国内视频| 欧美精品亚洲一区二区| 欧美最新免费一区二区三区| 91精品伊人久久大香线蕉| 欧美老熟妇乱子伦牲交| 狂野欧美激情性xxxx在线观看| 亚洲欧美色中文字幕在线| 亚洲综合色网址| 亚洲美女搞黄在线观看| 极品人妻少妇av视频| 丰满饥渴人妻一区二区三| 如日韩欧美国产精品一区二区三区| 国产精品久久久久成人av| 亚洲欧美精品自产自拍| av在线播放精品| 成人亚洲欧美一区二区av| 又黄又粗又硬又大视频| 国产白丝娇喘喷水9色精品| 最近2019中文字幕mv第一页| 欧美日韩成人在线一区二区| 欧美精品高潮呻吟av久久| 欧美日韩视频高清一区二区三区二| 久久鲁丝午夜福利片| 欧美精品亚洲一区二区| 日本vs欧美在线观看视频| av在线老鸭窝| 国产av国产精品国产| 一级片免费观看大全| 下体分泌物呈黄色| 又大又黄又爽视频免费| kizo精华| 啦啦啦中文免费视频观看日本| 母亲3免费完整高清在线观看 | 免费黄色在线免费观看| 日韩av不卡免费在线播放| 女人被躁到高潮嗷嗷叫费观| 免费久久久久久久精品成人欧美视频 | 精品人妻在线不人妻| 九色成人免费人妻av| 免费av不卡在线播放| 激情视频va一区二区三区| av福利片在线| 国产又色又爽无遮挡免| 国产精品久久久久成人av| 亚洲av日韩在线播放| 自线自在国产av| 亚洲av电影在线观看一区二区三区| 亚洲性久久影院| 亚洲国产精品成人久久小说| 少妇人妻 视频| 伊人久久国产一区二区| 日韩成人av中文字幕在线观看| av网站免费在线观看视频| 久久久久国产精品人妻一区二区| 最近最新中文字幕免费大全7| 欧美日本中文国产一区发布| 在线观看免费高清a一片| 亚洲精品视频女| 精品亚洲成国产av| 亚洲国产欧美在线一区| 久久久国产精品麻豆| 性色av一级| 1024视频免费在线观看| 人妻系列 视频| 欧美日韩视频精品一区| av在线app专区| av在线播放精品| 最新的欧美精品一区二区| 久久国内精品自在自线图片| 汤姆久久久久久久影院中文字幕| 久久久久久人人人人人| 国产 一区精品|