• <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è)定的臨床意義
    一区在线观看完整版| 嫁个100分男人电影在线观看| av福利片在线| 两个人看的免费小视频| 色婷婷久久久亚洲欧美| 老熟女久久久| 国产一区二区三区在线臀色熟女 | 他把我摸到了高潮在线观看 | 精品少妇一区二区三区视频日本电影| 免费少妇av软件| 黑丝袜美女国产一区| 亚洲一卡2卡3卡4卡5卡精品中文| 搡老乐熟女国产| 美女国产高潮福利片在线看| 国产真人三级小视频在线观看| 国产精品自产拍在线观看55亚洲 | 国产精品电影一区二区三区 | 亚洲精品久久成人aⅴ小说| 大陆偷拍与自拍| 免费在线观看黄色视频的| 黑人猛操日本美女一级片| 免费日韩欧美在线观看| av天堂在线播放| 啦啦啦 在线观看视频| 久久精品亚洲av国产电影网| 黄色a级毛片大全视频| 亚洲五月色婷婷综合| 久久精品熟女亚洲av麻豆精品| 99久久国产精品久久久| 美女扒开内裤让男人捅视频| 九色亚洲精品在线播放| 99re在线观看精品视频| 精品国产乱码久久久久久男人| bbb黄色大片| av国产精品久久久久影院| 久久久精品区二区三区| 他把我摸到了高潮在线观看 | 少妇粗大呻吟视频| 18禁国产床啪视频网站| 电影成人av| 亚洲综合色网址| 我要看黄色一级片免费的| av天堂在线播放| 日本黄色日本黄色录像| 久久免费观看电影| 91麻豆av在线| 一边摸一边抽搐一进一出视频| 国产一区有黄有色的免费视频| 日韩视频一区二区在线观看| 狂野欧美激情性xxxx| 搡老熟女国产l中国老女人| 日韩制服丝袜自拍偷拍| av有码第一页| e午夜精品久久久久久久| 色播在线永久视频| 亚洲全国av大片| 亚洲精品美女久久久久99蜜臀| 大香蕉久久成人网| 国产高清国产精品国产三级| 日韩一区二区三区影片| 国产不卡av网站在线观看| 最近最新免费中文字幕在线| 精品国产乱子伦一区二区三区| 久久婷婷成人综合色麻豆| 亚洲av美国av| 热99re8久久精品国产| 一级毛片精品| 丁香六月天网| 国内毛片毛片毛片毛片毛片| 亚洲avbb在线观看| 国产精品国产高清国产av | 99国产精品一区二区蜜桃av | 下体分泌物呈黄色| 男女无遮挡免费网站观看| 高清黄色对白视频在线免费看| 婷婷丁香在线五月| 热99国产精品久久久久久7| 欧美精品高潮呻吟av久久| 国产伦人伦偷精品视频| 亚洲国产av新网站| 王馨瑶露胸无遮挡在线观看| 新久久久久国产一级毛片| 黄色视频,在线免费观看| 性少妇av在线| 在线十欧美十亚洲十日本专区| 欧美精品高潮呻吟av久久| 男女免费视频国产| 国产精品免费大片| 精品福利观看| www.999成人在线观看| 亚洲天堂av无毛| 亚洲 欧美一区二区三区| 12—13女人毛片做爰片一| 日韩三级视频一区二区三区| 久久精品国产99精品国产亚洲性色 | 真人做人爱边吃奶动态| 中文字幕制服av| 亚洲精品国产精品久久久不卡| 精品熟女少妇八av免费久了| 精品少妇内射三级| 亚洲精品在线观看二区| 成人免费观看视频高清| 欧美精品一区二区大全| 天堂中文最新版在线下载| 大香蕉久久网| videosex国产| 视频区欧美日本亚洲| 黄色视频不卡| 最近最新免费中文字幕在线| 午夜老司机福利片| 免费av中文字幕在线| 亚洲午夜理论影院| 成人特级黄色片久久久久久久 | 51午夜福利影视在线观看| 久久九九热精品免费| 婷婷丁香在线五月| 亚洲成人手机| 国产视频一区二区在线看| 国产野战对白在线观看| 一二三四社区在线视频社区8| 久久亚洲真实| 桃花免费在线播放| 欧美 亚洲 国产 日韩一| 少妇精品久久久久久久| 一级片免费观看大全| 国产日韩一区二区三区精品不卡| 久久精品国产a三级三级三级| 91成年电影在线观看| 国产精品亚洲av一区麻豆| 757午夜福利合集在线观看| 热99久久久久精品小说推荐| 国产精品美女特级片免费视频播放器 | 757午夜福利合集在线观看| 久久久精品区二区三区| 国产欧美日韩一区二区三区在线| 老汉色av国产亚洲站长工具| 亚洲人成电影观看| 美女扒开内裤让男人捅视频| 国产成人av教育| 久久久久精品国产欧美久久久| 欧美日韩国产mv在线观看视频| 12—13女人毛片做爰片一| avwww免费| 黄色视频不卡| av福利片在线| 王馨瑶露胸无遮挡在线观看| 黄色成人免费大全| 亚洲,欧美精品.| 成人永久免费在线观看视频 | 亚洲精品美女久久久久99蜜臀| 国产精品成人在线| 亚洲欧美日韩高清在线视频 | 久久亚洲真实| 日本wwww免费看| 我的亚洲天堂| 欧美日韩精品网址| 99久久人妻综合| 免费女性裸体啪啪无遮挡网站| 大香蕉久久成人网| 自拍欧美九色日韩亚洲蝌蚪91| 国产人伦9x9x在线观看| 正在播放国产对白刺激| 高清欧美精品videossex| 欧美激情久久久久久爽电影 | 国产在线一区二区三区精| 日韩视频在线欧美| 91精品三级在线观看| 国产av精品麻豆| 一级片免费观看大全| 久久久久久久久久久久大奶| 热99re8久久精品国产| 夜夜夜夜夜久久久久| 法律面前人人平等表现在哪些方面| 久久人人爽av亚洲精品天堂| 成人18禁在线播放| 久久精品亚洲精品国产色婷小说| 看免费av毛片| 十八禁人妻一区二区| 三级毛片av免费| 人人澡人人妻人| 岛国在线观看网站| 十八禁高潮呻吟视频| 成年动漫av网址| 国产成人精品在线电影| 欧美亚洲 丝袜 人妻 在线| 高清av免费在线| 中文字幕精品免费在线观看视频| 丝瓜视频免费看黄片| 日韩三级视频一区二区三区| 亚洲一区二区三区欧美精品| 色婷婷av一区二区三区视频| 一本综合久久免费| 三上悠亚av全集在线观看| 国产三级黄色录像| 精品国产乱码久久久久久小说| 9191精品国产免费久久| av线在线观看网站| 俄罗斯特黄特色一大片| 色综合欧美亚洲国产小说| 热99国产精品久久久久久7| 国产精品久久久av美女十八| 日本精品一区二区三区蜜桃| 亚洲人成77777在线视频| 老司机福利观看| 欧美久久黑人一区二区| 色综合婷婷激情| 色老头精品视频在线观看| 国产在线一区二区三区精| 国产三级黄色录像| 黄色片一级片一级黄色片| 久久久精品94久久精品| 人人妻人人澡人人爽人人夜夜| 亚洲av国产av综合av卡| 搡老乐熟女国产| 建设人人有责人人尽责人人享有的| 色在线成人网| 天天操日日干夜夜撸| 国产日韩一区二区三区精品不卡| 亚洲精品中文字幕在线视频| 亚洲精品美女久久av网站| 久久久精品国产亚洲av高清涩受| 99精品在免费线老司机午夜| 国产视频一区二区在线看| 精品视频人人做人人爽| 精品福利永久在线观看| 成人三级做爰电影| 免费一级毛片在线播放高清视频 | 丰满饥渴人妻一区二区三| 国产免费av片在线观看野外av| 国产精品一区二区在线观看99| 69精品国产乱码久久久| 精品国内亚洲2022精品成人 | 国产精品欧美亚洲77777| 午夜91福利影院| 欧美日韩黄片免| 亚洲视频免费观看视频| a级毛片在线看网站| 国产av国产精品国产| 欧美成人午夜精品| 天天躁夜夜躁狠狠躁躁| 精品人妻熟女毛片av久久网站| 国产午夜精品久久久久久| 99香蕉大伊视频| 亚洲精品久久午夜乱码| 久热爱精品视频在线9| 巨乳人妻的诱惑在线观看| av有码第一页| 女人精品久久久久毛片| 激情在线观看视频在线高清 | 丁香欧美五月| 亚洲欧美一区二区三区黑人| 法律面前人人平等表现在哪些方面| www日本在线高清视频| 国产精品国产av在线观看| 欧美乱妇无乱码| 国产高清videossex| 99热网站在线观看| 91麻豆av在线| 国产av国产精品国产| 欧美亚洲日本最大视频资源| 亚洲视频免费观看视频| 国产黄频视频在线观看| 丰满迷人的少妇在线观看| 老司机在亚洲福利影院| 国产淫语在线视频| 亚洲欧美一区二区三区黑人| 高清av免费在线| 精品国产超薄肉色丝袜足j| 深夜精品福利| 国产成人影院久久av| 国产精品一区二区免费欧美| 色老头精品视频在线观看| 丁香六月天网| 精品福利永久在线观看| 大片电影免费在线观看免费| 热99国产精品久久久久久7| 少妇精品久久久久久久| 啪啪无遮挡十八禁网站| 黑人猛操日本美女一级片| 每晚都被弄得嗷嗷叫到高潮| 日本黄色日本黄色录像| 国产极品粉嫩免费观看在线| 久久青草综合色| 999久久久国产精品视频| 久久中文字幕一级| 99久久精品国产亚洲精品| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美激情在线| 久久亚洲精品不卡| 亚洲精品自拍成人| 亚洲精品乱久久久久久| 午夜福利影视在线免费观看| 亚洲av成人不卡在线观看播放网| 国产精品自产拍在线观看55亚洲 | 一区在线观看完整版| 少妇猛男粗大的猛烈进出视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品美女久久av网站| 国产精品偷伦视频观看了| 热re99久久国产66热| 亚洲欧美色中文字幕在线| 51午夜福利影视在线观看| svipshipincom国产片| 少妇 在线观看| 男女床上黄色一级片免费看| 欧美日韩成人在线一区二区| 男女高潮啪啪啪动态图| 国产深夜福利视频在线观看| 最近最新中文字幕大全免费视频| 正在播放国产对白刺激| 午夜精品国产一区二区电影| 人妻一区二区av| 日韩大片免费观看网站| 激情视频va一区二区三区| 久久精品人人爽人人爽视色| 精品欧美一区二区三区在线| 在线 av 中文字幕| 99在线人妻在线中文字幕 | 精品久久蜜臀av无| 国产成+人综合+亚洲专区| 国产亚洲精品第一综合不卡| a在线观看视频网站| 午夜免费鲁丝| 国产深夜福利视频在线观看| 69av精品久久久久久 | 日韩成人在线观看一区二区三区| 国产精品一区二区在线观看99| 精品国产一区二区久久| 亚洲一区中文字幕在线| 精品国产一区二区三区四区第35| 国产免费福利视频在线观看| 日韩中文字幕视频在线看片| 久久婷婷成人综合色麻豆| 国产精品久久久久久人妻精品电影 | 丰满少妇做爰视频| 精品一品国产午夜福利视频| 亚洲一区二区三区欧美精品| 91精品国产国语对白视频| 大型黄色视频在线免费观看| 日韩熟女老妇一区二区性免费视频| 精品熟女少妇八av免费久了| 咕卡用的链子| 人妻 亚洲 视频| a在线观看视频网站| 99久久99久久久精品蜜桃| tube8黄色片| 亚洲国产av影院在线观看| 18禁观看日本| 亚洲伊人久久精品综合| 国产精品影院久久| 日韩一卡2卡3卡4卡2021年| 日韩免费av在线播放| av福利片在线| 乱人伦中国视频| 日本wwww免费看| 欧美+亚洲+日韩+国产| 久久久国产成人免费| 热re99久久精品国产66热6| 亚洲五月婷婷丁香| 女同久久另类99精品国产91| 国产黄色免费在线视频| 国产视频一区二区在线看| 国产黄色免费在线视频| 欧美日韩精品网址| 国产免费现黄频在线看| 精品国产一区二区久久| 亚洲av日韩在线播放| 精品一区二区三区av网在线观看 | 人妻 亚洲 视频| 久久久久国内视频| 精品国内亚洲2022精品成人 | 日韩免费高清中文字幕av| 精品亚洲成国产av| 国产精品av久久久久免费| 欧美另类亚洲清纯唯美| 日韩免费av在线播放| 夜夜爽天天搞| av国产精品久久久久影院| av欧美777| 在线观看免费午夜福利视频| 午夜精品国产一区二区电影| aaaaa片日本免费| 精品亚洲成a人片在线观看| 99热网站在线观看| 日韩成人在线观看一区二区三区| 欧美另类亚洲清纯唯美| 80岁老熟妇乱子伦牲交| 欧美大码av| 一本一本久久a久久精品综合妖精| 淫妇啪啪啪对白视频| 亚洲,欧美精品.| 黄色a级毛片大全视频| 亚洲精品美女久久av网站| 脱女人内裤的视频| 人妻一区二区av| 日本av手机在线免费观看| 十八禁网站免费在线| 亚洲精品av麻豆狂野| 欧美日韩中文字幕国产精品一区二区三区 | 男女无遮挡免费网站观看| 国产成人啪精品午夜网站| 欧美黑人精品巨大| 91大片在线观看| 十八禁网站免费在线| 国产av国产精品国产| 国产免费现黄频在线看| 黄色毛片三级朝国网站| 久久精品人人爽人人爽视色| 老司机深夜福利视频在线观看| 国产成人精品无人区| 下体分泌物呈黄色| 十八禁人妻一区二区| 国产精品久久久久久精品古装| 成人黄色视频免费在线看| 天天躁夜夜躁狠狠躁躁| 亚洲免费av在线视频| 天天影视国产精品| 欧美另类亚洲清纯唯美| 亚洲视频免费观看视频| 日本av手机在线免费观看| 国产欧美亚洲国产| 欧美日韩国产mv在线观看视频| 国内毛片毛片毛片毛片毛片| 欧美日韩精品网址| 国产色视频综合| 999久久久国产精品视频| 国产一区二区在线观看av| 国产无遮挡羞羞视频在线观看| 中亚洲国语对白在线视频| 国产亚洲av高清不卡| 国产黄频视频在线观看| 久久久久久久国产电影| 久久久久久亚洲精品国产蜜桃av| 久久亚洲精品不卡| 人人澡人人妻人| 女人久久www免费人成看片| 国产成人精品久久二区二区免费| 亚洲av日韩在线播放| 在线观看66精品国产| 国产av又大| 亚洲三区欧美一区| 1024视频免费在线观看| 妹子高潮喷水视频| 日本一区二区免费在线视频| 国产成人系列免费观看| 欧美中文综合在线视频| 男女床上黄色一级片免费看| av一本久久久久| 亚洲伊人色综图| 精品国产乱码久久久久久小说| 青青草视频在线视频观看| 我要看黄色一级片免费的| 捣出白浆h1v1| 日韩中文字幕视频在线看片| 99久久人妻综合| 成年人黄色毛片网站| 天堂俺去俺来也www色官网| 最近最新免费中文字幕在线| 国产精品国产高清国产av | 欧美激情高清一区二区三区| 热re99久久国产66热| 亚洲,欧美精品.| 欧美另类亚洲清纯唯美| 丁香六月天网| 老司机午夜十八禁免费视频| 超碰97精品在线观看| 男女免费视频国产| 国产成+人综合+亚洲专区| 国产亚洲欧美在线一区二区| 一级毛片女人18水好多| 欧美老熟妇乱子伦牲交| 国产激情久久老熟女| 亚洲国产毛片av蜜桃av| 国产欧美日韩综合在线一区二区| 国产亚洲av高清不卡| 中文字幕人妻丝袜制服| 91大片在线观看| h视频一区二区三区| 一本—道久久a久久精品蜜桃钙片| 亚洲av成人不卡在线观看播放网| 久9热在线精品视频| 久久久久久免费高清国产稀缺| 又大又爽又粗| 一级a爱视频在线免费观看| 亚洲av欧美aⅴ国产| 亚洲国产欧美网| 91成年电影在线观看| 99国产精品一区二区三区| 高清毛片免费观看视频网站 | 免费一级毛片在线播放高清视频 | 美女视频免费永久观看网站| 亚洲成国产人片在线观看| 日韩视频在线欧美| 人人澡人人妻人| 国产av精品麻豆| 久热爱精品视频在线9| 国产极品粉嫩免费观看在线| 久久精品熟女亚洲av麻豆精品| 人人妻人人添人人爽欧美一区卜| 天天添夜夜摸| 亚洲欧美精品综合一区二区三区| 别揉我奶头~嗯~啊~动态视频| 国产成人系列免费观看| 淫妇啪啪啪对白视频| 亚洲九九香蕉| 一区二区三区激情视频| 99久久精品国产亚洲精品| 免费黄频网站在线观看国产| 丰满迷人的少妇在线观看| 国产男女超爽视频在线观看| 亚洲欧美日韩另类电影网站| 女警被强在线播放| 欧美黑人精品巨大| 精品高清国产在线一区| 男女边摸边吃奶| 欧美+亚洲+日韩+国产| 亚洲国产中文字幕在线视频| av有码第一页| 国产精品98久久久久久宅男小说| 三级毛片av免费| 两个人免费观看高清视频| 少妇猛男粗大的猛烈进出视频| 国产精品久久久久久人妻精品电影 | 18禁美女被吸乳视频| 久久久久久免费高清国产稀缺| 十分钟在线观看高清视频www| 18禁国产床啪视频网站| 精品卡一卡二卡四卡免费| 久久天躁狠狠躁夜夜2o2o| 狠狠精品人妻久久久久久综合| 国产男女内射视频| 美女高潮到喷水免费观看| 岛国在线观看网站| 操出白浆在线播放| av在线播放免费不卡| 日日摸夜夜添夜夜添小说| 夜夜爽天天搞| 亚洲av成人不卡在线观看播放网| 国产片内射在线| 久久久国产精品麻豆| 亚洲专区国产一区二区| 日本撒尿小便嘘嘘汇集6| 亚洲精品一二三| 美国免费a级毛片| av超薄肉色丝袜交足视频| 麻豆成人av在线观看| 日本a在线网址| 十八禁人妻一区二区| 国产精品久久久久成人av| 亚洲一码二码三码区别大吗| 免费在线观看视频国产中文字幕亚洲| 男女之事视频高清在线观看| 国产老妇伦熟女老妇高清| 交换朋友夫妻互换小说| 在线观看舔阴道视频| 国产成人精品在线电影| 亚洲色图综合在线观看| 99国产精品免费福利视频| 2018国产大陆天天弄谢| 午夜福利欧美成人| 法律面前人人平等表现在哪些方面| 精品人妻在线不人妻| 国产精品一区二区在线不卡| 国产精品二区激情视频| 中文字幕av电影在线播放| videosex国产| 国产亚洲精品一区二区www | 久久久水蜜桃国产精品网| 黄色毛片三级朝国网站| 大陆偷拍与自拍| 欧美精品啪啪一区二区三区| 1024视频免费在线观看| 国产无遮挡羞羞视频在线观看| 国产欧美日韩综合在线一区二区| 男女午夜视频在线观看| 欧美另类亚洲清纯唯美| 亚洲成人国产一区在线观看| 国产精品电影一区二区三区 | 国产成人精品久久二区二区91| 亚洲色图av天堂| 狠狠狠狠99中文字幕| 一边摸一边抽搐一进一出视频| 色精品久久人妻99蜜桃| 国产亚洲一区二区精品| 日本撒尿小便嘘嘘汇集6| 人人妻人人爽人人添夜夜欢视频| 在线 av 中文字幕| 亚洲国产av新网站| 高清视频免费观看一区二区| 亚洲午夜精品一区,二区,三区| 久久这里只有精品19| 亚洲av美国av| 下体分泌物呈黄色| 亚洲成a人片在线一区二区| 亚洲精品美女久久久久99蜜臀| 欧美日韩精品网址| e午夜精品久久久久久久| 国产欧美日韩精品亚洲av| 九色亚洲精品在线播放| 成人国产av品久久久| 91大片在线观看| 51午夜福利影视在线观看| 伊人久久大香线蕉亚洲五| 国产人伦9x9x在线观看| 久久精品国产亚洲av高清一级| 丝瓜视频免费看黄片| 久久久久久久久免费视频了| 人成视频在线观看免费观看| 国产深夜福利视频在线观看| 精品一区二区三区av网在线观看 | 久久精品亚洲精品国产色婷小说| 久久香蕉激情| kizo精华| 国产精品国产高清国产av | 99国产精品一区二区三区|