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

    茶樹R2R3-MYB轉(zhuǎn)錄因子CsTT2表達(dá)分析及功能初步鑒定

    2022-08-27 03:24:16王玉源劉任堅(jiān)劉少群舒燦偉孫彬妹鄭鵬
    茶葉科學(xué) 2022年4期
    關(guān)鍵詞:生物

    王玉源,劉任堅(jiān),劉少群,舒燦偉,孫彬妹,鄭鵬*

    茶樹R2R3-MYB轉(zhuǎn)錄因子CsTT2表達(dá)分析及功能初步鑒定

    王玉源1,劉任堅(jiān)1,劉少群1,舒燦偉2,孫彬妹1,鄭鵬1*

    1. 華南農(nóng)業(yè)大學(xué)園藝學(xué)院,廣東 廣州 510642;2. 華南農(nóng)業(yè)大學(xué)植物保護(hù)學(xué)院/廣東省微生物信號與作物病害防控重點(diǎn)實(shí)驗(yàn)室/群體微生物研究中心,廣東 廣州 510642

    兒茶素是茶樹中特色的次生代謝產(chǎn)物之一,是影響茶葉的品質(zhì)與風(fēng)味的主要組分,具有抗氧化、抗病毒、降脂減肥等藥理功效。通過系統(tǒng)發(fā)育進(jìn)化樹分析、基因表達(dá)模式分析和分子生物學(xué)試驗(yàn)對茶樹兒茶素生物合成相關(guān)調(diào)控因子CsTT2的功能進(jìn)行初步鑒定。結(jié)果顯示,CsTT2是R2R3-MYB轉(zhuǎn)錄因子,與擬南芥中調(diào)控次生代謝產(chǎn)物的MYB轉(zhuǎn)錄因子同在一個(gè)分支。在茶樹品種頂芽組織中總兒茶素含量較高,和兒茶素生物合成相關(guān)基因的表達(dá)水平也較高。亞細(xì)胞定位、酵母試驗(yàn)和雙熒光素酶報(bào)告系統(tǒng)試驗(yàn)結(jié)果表明,定位在細(xì)胞核中,其編碼的蛋白是具有轉(zhuǎn)錄激活能力的調(diào)控因子,可以結(jié)合兒茶素生物合成關(guān)鍵基因的啟動(dòng)子激活其表達(dá)。

    茶樹;R2R3-MYB轉(zhuǎn)錄因子;;兒茶素生物合成

    茶樹[(L.) O. Kuntze]屬于山茶屬(Theaceae)灌木或小喬木,主要包括中國種(var)和阿薩姆種(var)兩個(gè)栽培品種[1],是中國、日本、印度和肯尼亞等地區(qū)重要的經(jīng)濟(jì)作物[2]。兒茶素屬于類黃酮物質(zhì),是茶樹中重要的次生代謝產(chǎn)物,在茶樹的各個(gè)器官中均有分布。兒茶素含量占茶葉干重的12%~24%,占多酚類總量的70%~80%[3]。兒茶素作為茶樹中重要的次生代謝產(chǎn)物之一,不僅能夠影響茶葉的品質(zhì)風(fēng)味,還具有抗癌[4]、抗衰老[5-6]、抗氧化[7]、抗病毒[8-9]、保護(hù)神經(jīng)[10]、調(diào)節(jié)免疫力[11]及降脂減肥[12]等藥理功能。

    在植物中,許多次生代謝產(chǎn)物的生物合成途徑是保守的。研究表明,茶樹中類黃酮生物合成途徑可以參照大多數(shù)模式植物,特別是擬南芥中花青素的生物合成途徑[13]。與植物類黃酮生物合成途徑比較,兒茶素類化合物種類較多,其生物合成途徑涉及多條代謝途徑,是相對復(fù)雜、需要受多種因素協(xié)調(diào)調(diào)控的過程。前人研究發(fā)現(xiàn),茶樹兒茶素生物合成通路主要包括莽草酸途徑、苯丙烷途徑、非酯型兒茶素合成和酯型兒茶素合成等4個(gè)環(huán)節(jié)[13-17]。目前,茶樹兒茶素生物合成過程中主要的催化酶基因已經(jīng)被克隆,包括苯丙氨酸解氨酶、肉桂酸羥化酶、4-香豆酸CoA連接酶、查爾酮合成酶、查爾酮異構(gòu)酶、黃烷酮3-羥化酶、3',5'-類黃酮羥化酶、二氫黃酮醇還原酶、無色花色素還原酶、花青素合酶、花青素還原酶等[18]。茶樹中的兒茶素又分為非酯型兒茶素和酯型兒茶素。非酯型兒茶素生物合成的主要途徑已經(jīng)清晰[19-21],而酯型兒茶素的生物合成通路尚未明確。有研究者推測,酯型兒茶素是由非酯型兒茶素沒食子酯化得到,但這一猜測并未被試驗(yàn)證實(shí)[13]。Niemetz和Gross認(rèn)為1--沒食子酰--葡萄糖是單寧合成的有效?;w和受體[22-23]。Liu等[24]在此基礎(chǔ)上進(jìn)一步證實(shí),表兒茶素(EC)和表沒食子兒茶素(EGC)在沒食子?;?1----葡萄糖轉(zhuǎn)移酶(UGGT)和1--沒食子酰基---沒食子?;D(zhuǎn)移酶(ECGT)的催化下合成酯型兒茶素。與此同時(shí),酯型兒茶素也會在水解酶的作用下水解得到非酯型兒茶素EC和EGC[25]。也有研究發(fā)現(xiàn),基因家族能夠參與非酯型兒茶素轉(zhuǎn)化為酯型兒茶素的過程[26]。近年來,有關(guān)茶樹兒茶素生物合成的調(diào)控已有相關(guān)報(bào)道。Luo等[27]研究發(fā)現(xiàn),CsWRKY31和CsWRKY48能夠結(jié)合甲基化兒茶素生物合成相關(guān)基因、和的啟動(dòng)子,抑制它們的表達(dá)。Wang等[28]研究發(fā)現(xiàn),在煙草中過表達(dá),轉(zhuǎn)基因植株中兒茶素單體的含量上升,且煙草中相關(guān)合成基因和的表達(dá)上調(diào)。綜上所述,茶樹兒茶素的生物合成途徑已經(jīng)取得了許多研究進(jìn)展,然而其生物合成的分子調(diào)控機(jī)制還需進(jìn)一步探究。

    MYB是植物中最大的轉(zhuǎn)錄因子家族,它在植物的生長發(fā)育、植物響應(yīng)生物和非生物脅迫以及調(diào)控植物次生代謝產(chǎn)物生物合成方面起著重要作用[29]。它通常具有兩個(gè)不同的區(qū)域,一個(gè)是位于N端的高度保守的DNA結(jié)合域,另一個(gè)是位于C端的負(fù)責(zé)調(diào)節(jié)蛋白質(zhì)活性的調(diào)節(jié)區(qū)域。根據(jù)結(jié)構(gòu)域所含R結(jié)構(gòu)的數(shù)量,植物中的MYB轉(zhuǎn)錄因子可分為4類,分別是1R-MYB/MYB-related(1個(gè)R結(jié)構(gòu))、R2R3-MYB(2個(gè)R結(jié)構(gòu))、R1R2R3-MYB蛋白(3個(gè)R結(jié)構(gòu))和4R-MYB(4個(gè)R結(jié)構(gòu))[29]。在植物中,大多數(shù)MYB屬于R2R3-MYB類轉(zhuǎn)錄因子[30],通常作為轉(zhuǎn)錄激活子或者抑制子直接或間接調(diào)控類黃酮生物合成中相關(guān)結(jié)構(gòu)基因的表達(dá),從而影響物質(zhì)的積累[31-33]。TT2類轉(zhuǎn)錄因子是R2R3-MYB蛋白,通過激活相關(guān)合成酶基因的表達(dá)調(diào)控次生代謝產(chǎn)物的生物合成[34-37]。目前,在許多植物中,TT2類轉(zhuǎn)錄因子已被克隆和鑒定,包括擬南芥[38]、葡萄[39]、草莓和[33]、蘋果和[40]、楊樹[41]等,這些基因都參與調(diào)控植物中次生代謝產(chǎn)物的生物合成。

    本研究通過系統(tǒng)發(fā)育進(jìn)化樹分析、基因表達(dá)模式分析和分子生物學(xué)試驗(yàn)對茶樹兒茶素生物合成相關(guān)調(diào)控因子CsTT2的功能進(jìn)行初步鑒定,為研究茶樹兒茶素生物合成的分子調(diào)控機(jī)理奠定理論基礎(chǔ)。

    1 材料與方法

    1.1 試驗(yàn)材料

    試驗(yàn)材料為嶺頭單叢(LTDC)、紅柄單叢(HBDC)和龍井長葉(LJCY)茶樹品種的多年生灌木茶苗,種植于華南農(nóng)業(yè)大學(xué)茶樹種質(zhì)資源圃。于2020年,分別采摘3個(gè)茶樹品種一芽二葉頂芽組織,用液氮速凍后置于–80℃冰箱貯存,用于提取RNA和測定兒茶素含量;采摘LTDC品種的一芽二葉頂芽、嫩葉、成熟葉、老葉、根、莖、花和果實(shí)組織,用液氮速凍后置于–80℃冰箱貯存,用于提取RNA;以上每個(gè)樣本設(shè)置3個(gè)生物學(xué)重復(fù)。此外,試驗(yàn)用到的本氏煙草(),種子來自本實(shí)驗(yàn)室。

    1.2 同源比對分析及系統(tǒng)發(fā)育進(jìn)化樹構(gòu)建

    將擬南芥AtTT2(AT5G35550)蛋白序列比對到茶樹基因組數(shù)據(jù)庫(http://tpdb.shengxin.ren/index.html)中,找到茶樹同源基因(),并下載其蛋白序列。將CsTT2的蛋白序列放入到NCBI數(shù)據(jù)庫中進(jìn)行比對,下載其他物種中同源蛋白序列,并在DNAMAN軟件中進(jìn)行蛋白序列比對分析。從擬南芥基因組數(shù)據(jù)庫(https://www.arabidopsis.org)下載擬南芥所有MYB蛋白序列,用clustalW軟件對CsTT2蛋白序列和擬南芥所有MYB蛋白序列進(jìn)行完全比對分析。根據(jù)已有報(bào)道,使用Jalview軟件截取R2R3-MYB蛋白的保守結(jié)構(gòu)域區(qū)域[42],再使用MEGA-X軟件中的鄰接法進(jìn)行系統(tǒng)發(fā)育進(jìn)化樹的構(gòu)建。

    1.3 兒茶素含量測定

    兒茶素(C)、EC、沒食子兒茶素(GC)、EGC、表兒茶素沒食子酸酯(ECG)和表沒食子兒茶素沒食子酸酯(EGCG)標(biāo)準(zhǔn)樣品購自上海源葉生物科技有限公司。用70%甲醇溶液將標(biāo)準(zhǔn)樣品分別配制成質(zhì)量濃度為0.01、0.1、0.2、0.5?mg·mL-1的標(biāo)準(zhǔn)液,用于兒茶素單體標(biāo)準(zhǔn)曲線的構(gòu)建。

    使用冷凍干燥機(jī)將樣品凍干,用研缽將樣品磨成粉末狀,準(zhǔn)確稱取0.1?g粉末加入到含有5?mL 70%甲醇的10?mL離心管中,浸提5?min后(期間上下振蕩),70℃水浴超聲振蕩30?min,在12?000?r·min-1下離心5?min,用1?mL注射器吸取上清液,通過0.45?μm微孔過濾器過濾到棕色樣品瓶中,最后用Waters Alliance 2489高效液相色譜儀(Waters Technologies,美國)測定樣品中兒茶素的含量。上機(jī)條件:柱溫,(30±5)℃;流速,1?mg·mL-1;波長,280?nm;進(jìn)樣量,10?μL;上機(jī)程序如表1所示。

    1.4 實(shí)時(shí)熒光定量PCR(qRT-PCR)

    表1 梯度程序

    1.5 基因表達(dá)模式分析

    從茶樹基因組數(shù)據(jù)庫(http://tpdb.shengxin.ren/index.html)中下載茶樹頂芽、嫩葉、成熟葉、老葉、根、莖、花和果實(shí)等8個(gè)不同組織部位的轉(zhuǎn)錄組數(shù)據(jù),用TBtools軟件進(jìn)行相關(guān)基因表達(dá)量熱圖的構(gòu)建[45],分析基因表達(dá)模式。

    1.6 載體構(gòu)建

    根據(jù)的CDS序列以及載體酶切位點(diǎn),使用CE Design V1軟件設(shè)計(jì)同源重組引物,引物見表3,下劃線部分為載體的酶切位點(diǎn)。以白葉單叢的cDNA為模板,使用PrimeSTAR酶進(jìn)行PCR擴(kuò)增。PCR產(chǎn)物經(jīng)過瓊脂糖凝膠電泳后回收目的條帶,連接到對應(yīng)的線性化載體上,轉(zhuǎn)化至DH5后進(jìn)行鑒定并測序,將陽性克隆提取質(zhì)粒,用于轉(zhuǎn)化農(nóng)桿菌和酵母菌。

    1.7 亞細(xì)胞定位試驗(yàn)

    將的CDS序列(去除終止密碼子)構(gòu)建到pEAQ-EGFP載體上,以空白載體為對照,使用熱激法轉(zhuǎn)化農(nóng)桿菌感受態(tài)GV3101。用侵染液(200?μmol·L-1AS,10?mmol·L-1MgCl2,10?mmol·L-1MES,PH為5.6,現(xiàn)配現(xiàn)用)將活化后的菌體重懸至OD值為0.4~0.6,活化2?h后,將空載體(35S::CsTT2::EGFP)與細(xì)胞核marker(35S::NLS::Dsred)按體積比1∶1混合,使用無針頭注射器將菌液從葉片下表皮注射到本氏煙草中,黑暗保濕培養(yǎng)48?h后,用生物顯微鏡BX53(OLYMPUS,日本)進(jìn)行觀察和拍照。

    1.8 酵母自激活試驗(yàn)

    將全長CDS序列以及根據(jù)保守結(jié)構(gòu)域短截后的序列構(gòu)建到pGBKT7(BD)載體上,以空白載體為對照,將重組產(chǎn)物轉(zhuǎn)化到酵母感受態(tài)Y2HGold。用100?μL ddH2O重懸陽性單克隆菌,分別在SD/-Trp和SD/-Trp-His-Ade培養(yǎng)基上點(diǎn)5?μL菌液,30℃倒置培養(yǎng)3?d后,用X--gal對菌斑進(jìn)行染色,然后觀察和拍照。

    1.9 轉(zhuǎn)錄激活試驗(yàn)

    將全長CDS序列構(gòu)建到pEAQ-PBD載體上,以空白載體為陰性對照,pEAQ-PBD-VP16為陽性對照,將重組產(chǎn)物轉(zhuǎn)化到農(nóng)桿菌感受態(tài)GV3101。用侵染液(200?μmol·L-1AS,10?mmol·L-1MgCl2,10?mmol·L-1MES,PH為5.6,現(xiàn)配現(xiàn)用)將活化后的菌體重懸至OD值為0.2,活化2?h后,將陰性對照/陽性對照/pEAQ-PBD-CsTT2與報(bào)告載體(GAL-LUC)按體積比1∶1混合,使用無針頭注射器將菌液從葉片下表皮注射到本氏煙草中,黑暗保濕培養(yǎng)24?h,再正常培養(yǎng)48?h后,使用Promega公司的雙熒光素酶報(bào)告系統(tǒng)試劑盒檢測螢火蟲熒光素酶(Firefly luciferase,LUC)和海腎熒光素酶(Renilla luciferase,REN)的活性。

    表2 定量引物序列

    表3 載體引物序列

    1.10 酵母單雜交試驗(yàn)

    使用在線數(shù)據(jù)庫JASPAR(https://jaspar.genereg.net)對兒茶素合成基因啟動(dòng)子中可能存在的MYB結(jié)合元件進(jìn)行預(yù)測和統(tǒng)計(jì)。

    根據(jù)預(yù)測結(jié)果,將兒茶素合成基因啟動(dòng)子短截為4個(gè)片段,分別構(gòu)建到pAbAi載體上,并整合到Y(jié)1HGold酵母菌中,形成誘餌特異性報(bào)告菌株,篩選每個(gè)啟動(dòng)子的自激活抑制濃度(其中,ProLAR和ProSCPL的自激活活性在AbA質(zhì)量濃度為1?000?ng·mL-1下仍無法抑制,因此沒進(jìn)行后續(xù)試驗(yàn))。將基因的CDS全長序列構(gòu)建到pGADT7(AD)載體上,再轉(zhuǎn)化到誘餌特異性報(bào)告菌株中。用ddH2O重懸陽性單克隆菌,將濃度OD600調(diào)為0.2、0.02、0.002,分別在SD/-Leu和SD/-Leu+AbA培養(yǎng)基上點(diǎn)5?μL菌液,30℃倒置培養(yǎng)3?d后,進(jìn)行觀察和拍照。

    1.11 雙熒光素酶試驗(yàn)

    將全長CDS序列構(gòu)建到pEAQ載體上,將兒茶素合成基因啟動(dòng)子(、、)構(gòu)建到pGreenII-0800-LUC載體上,以空白載體pEAQ作為對照,將重組產(chǎn)物轉(zhuǎn)化到農(nóng)桿菌感受態(tài)GV3101。用侵染液(200?μmol·L-1AS,10?mmol·L-1MgCl2,10?mmol·L-1MES,pH為5.6,現(xiàn)配現(xiàn)用)將活化后的菌體重懸至OD值為0.2,活化2?h后,將pEAQ/pEAQ-CsTT2+pGreenII-0800-ProANR-LUC/pGreenII-0800-ProLAR-LUC/pGreenII-0800-ProSCPL-LUC分別按體積比20∶1混合,使用無針頭注射器將菌液從葉片下表皮注射到本氏煙草中,黑暗保濕培養(yǎng)24?h,再正常培養(yǎng)48?h后,使用Promega公司的雙熒光素酶報(bào)告系統(tǒng)試劑盒檢測Firefly和REN值。

    1.12 數(shù)據(jù)統(tǒng)計(jì)與分析

    數(shù)據(jù)方差分析使用IBM SPSS Statistics 22軟件完成,圖片處理和組合使用Adobe Illustrator 2020軟件完成。

    2 結(jié)果與分析

    2.1 CsTT2保守結(jié)構(gòu)域及系統(tǒng)發(fā)育進(jìn)化樹分析

    多序列比對結(jié)果顯示,CsTT2具有兩個(gè)串聯(lián)不完整重復(fù)序列R2和R3,表明CsTT2是R2R3-MYB轉(zhuǎn)錄因子(圖1)。在擬南芥中,根據(jù)DNA結(jié)合域和C端氨基酸基序的保守性,對MYB轉(zhuǎn)錄因子家族進(jìn)行了分組;其中,R2R3-MYB蛋白被分成23個(gè)亞組[46]。本研究對CsTT2和擬南芥中所有MYB蛋白進(jìn)行了系統(tǒng)發(fā)育進(jìn)化樹分析,將蛋白分為6個(gè)主要亞組,Group Ⅵ又進(jìn)一步分為Group Ⅵa、Group Ⅵb、Group Ⅵc、Group Ⅵd和Group Ⅵe。結(jié)果顯示,CsTT2和擬南芥S2、S3、S5、S6、S7和S15中的R2R3-MYB蛋白聚集在同一亞組(Group Ⅳ)中(圖2)。S7中的成員可以調(diào)控黃酮醇的合成[47],S6中的成員可以調(diào)控營養(yǎng)生殖過程中花青素的合成[48],S5中的成員可以調(diào)控?cái)M南芥種皮原花青素的合成[49],而S3中的成員可以激活導(dǎo)管中的木質(zhì)素生物合成[50-51]。因此,CsTT2可能參與調(diào)控茶樹類黃酮次生代謝物生物合成的過程。

    2.2 茶樹頂芽兒茶素含量與兒茶素生物合成相關(guān)基因表達(dá)水平分析

    兒茶素主要分布在茶樹的頂芽部位。本研究使用高效液相色譜法測定LTDC、HBDC和LJCY 3個(gè)茶樹品種頂芽部位中6種常見兒茶素的含量(圖3-A)。如圖3-B所示,在3個(gè)茶樹品種中,GC的含量沒有顯著差異,而HBDC中EGC和C的含量顯著高于其他兩個(gè)品種,EC的含量顯著高于LJCY。此外,LTDC和HBDC中EGCG的含量顯著高于LJCY,而ECG的含量則相反。

    注:葡萄:VvTT2(RVW44110.1),木豆:CcMYB6(XP_020227123.1),楊梅:MrTT2(KAB1209459.1),雞血藤:SsTT2(TKY54429.1),大豆:GmMYB23(XP_003541301.1),擬南芥AtTT2(AT5G35550)

    注:擬南芥中的R2R3-MYB蛋白被分為23個(gè)亞組[46],每個(gè)亞組對應(yīng)分支的顏色如圖右上方所示,而黑色分支為除R2R3-MYB外的其他MYB成員。本研究構(gòu)建的系統(tǒng)發(fā)育進(jìn)化樹主要分為6個(gè)亞組,對應(yīng)分組情況如圖紫色弧線所示

    注:A. 液相色譜峰圖。B. 3個(gè)茶樹品種中6種兒茶素單體的含量。數(shù)據(jù)是用平均值±標(biāo)準(zhǔn)誤來表示,圖中不同的字母表示顯著不同的值(P<0.05, Tukey’s test)

    Fig.3 Contents of catechincomponents in apical bud of different tea cultivars

    進(jìn)一步對3個(gè)茶樹品種頂芽的總兒茶素含量以及兒茶素生物合成相關(guān)基因的表達(dá)量進(jìn)行分析。結(jié)果顯示,LTDC和HBDC中總兒茶素含量顯著高于LJCY(圖4-A)。此外,在LTDC和HBDC中,除了()外,目的基因和兒茶素生物合成相關(guān)基因()、()、()、()和()的表達(dá)量均顯著高于LJCY(圖4-B和圖4-C)。

    2.3 CsTT2轉(zhuǎn)錄活性分析

    為了進(jìn)一步分析CsTT2的轉(zhuǎn)錄活性,本研究還進(jìn)行了亞細(xì)胞定位、雙熒光素酶報(bào)告系統(tǒng)試驗(yàn)和酵母自激活試驗(yàn)。結(jié)果顯示,轉(zhuǎn)入pEAQ-EGFP空載體的本氏煙草葉片中的細(xì)胞膜和細(xì)胞核均激發(fā)出綠色熒光,轉(zhuǎn)入pEAQ-CsTT2-EGFP的本氏煙草葉片中僅細(xì)胞核激發(fā)出綠色熒光(圖5-A),這表明定位在細(xì)胞核中。雙熒光素酶報(bào)告系統(tǒng)結(jié)果顯示,相對于pEAQ-PBD陰性對照,CsTT2可以顯著提高Firefly LUC的活性(圖5-B),表明CsTT2在植物中具有轉(zhuǎn)錄激活功能。為了進(jìn)一步研究CsTT2的激活區(qū)域,根據(jù)保守結(jié)構(gòu)域的位置將CsTT2序列短截后進(jìn)行酵母自激活試驗(yàn)。結(jié)果顯示,CsTT2蛋白的激活區(qū)域主要位于C端的113—308氨基酸(圖5-C)。上述結(jié)果表明,定位在細(xì)胞核中,是具有轉(zhuǎn)錄激活能力的調(diào)控基因,且轉(zhuǎn)錄激活區(qū)域主要位于C端的113—308氨基酸處。

    2.4 CsTT2結(jié)合調(diào)控兒茶素合成基因ANR啟動(dòng)子活性的分析

    對和兒茶素生物合成相關(guān)基因(、、、、和)在茶樹8個(gè)不同組織部位中的轉(zhuǎn)錄組數(shù)據(jù)進(jìn)行分析。結(jié)果顯示,()、()和()在茶樹8個(gè)不同組織部位中的表達(dá)模式和一致,主要在茶樹的幼嫩芽葉部位中表達(dá),且表達(dá)量在嫩葉組織中最高(圖6-A)。qRT-PCR結(jié)果進(jìn)一步證明,、()、()和()在LTDC 8個(gè)不同組織部位中的表達(dá)模式一致,均在頂芽和嫩葉中高表達(dá)(圖6-B)。

    注:A. 3個(gè)茶樹品種頂芽中總兒茶素的含量。B. CsTT2在3個(gè)茶樹品種頂芽中的表達(dá)量。C. 兒茶素生物合成相關(guān)基因(LAR、ANR、DFR、CHS、ANS和SCPL)在3個(gè)茶樹品種頂芽中的表達(dá)量。數(shù)據(jù)是用平均值±標(biāo)準(zhǔn)誤來表示,采用Student’s test與對照“LJCY”進(jìn)行顯著性分析(**,P<0.01)

    注:A. CsTT2的亞細(xì)胞定位分析。比例尺=50?μm。Bright field:明場通道;DsRed channel:紅色熒光通道;GFP channel:綠色熒光通道;Merged channel:融合通道。B. CsTT2在煙草中的轉(zhuǎn)錄活性分析。pEAQ-PBD和pEAQ-PBD-VP16為陰性和陽性對照。數(shù)據(jù)是用平均值±標(biāo)準(zhǔn)誤來表示,采用Student’s test與陰性對照進(jìn)行顯著性分析(**,P<0.01)。C. CsTT2在酵母中的轉(zhuǎn)錄活性分析。BD載體為陰性對照。圖上方為CsTT2蛋白結(jié)構(gòu)圖。SD/-Trp代表SD培養(yǎng)基缺少色氨酸;SD/-Trp-His-Ade/+X-α-Gal代表SD培養(yǎng)基缺少色氨酸、組氨酸、腺嘌呤且添加X-α-Gal

    Fig.5 Subcellular localization and transcription activity analysis of CsTT2

    進(jìn)一步通過酵母單雜交試驗(yàn)驗(yàn)證CsTT2與啟動(dòng)子的結(jié)合能力。如圖7-A所示,克隆出來的啟動(dòng)子區(qū)域(–1?bp至–1?423?bp)可能存在16個(gè)MYB-bingding-element。根據(jù)預(yù)測結(jié)果,將其啟動(dòng)子短截成4段約為350?bp的短片段,進(jìn)行后續(xù)試驗(yàn)。與AD+pAbAi-ProANR2酵母菌相比,轉(zhuǎn)化AD-CsTT2+pAbAi-ProANR2的菌株能在SD/Leu和SD/Leu+AbA的平板上生長,這說明CsTT2能夠與的啟動(dòng)子相結(jié)合(圖7-B)。此外,雙熒光素酶報(bào)告系統(tǒng)試驗(yàn)結(jié)果顯示,在植物體內(nèi),CsTT2能夠激活啟動(dòng)子,進(jìn)而驅(qū)動(dòng)下游螢火蟲熒光素酶報(bào)告基因的表達(dá),增強(qiáng)Firefly LUC的活性;而它不能激活和的啟動(dòng)子(圖7-C)。上述結(jié)果表明,CsTT2蛋白能夠結(jié)合兒茶素生物合成相關(guān)基因的啟動(dòng)子,進(jìn)而激活的表達(dá)。

    3 討論

    茶樹兒茶素不僅能決定茶葉的品質(zhì)和風(fēng)味,還具有抗氧化、抗癌和降脂減肥等藥理功效。目前,茶樹兒茶素生物合成途徑的研究已取得了重要的進(jìn)展,然而其分子調(diào)控機(jī)制尚未清晰。許多研究表明,R2R3-MYB轉(zhuǎn)錄因子在植物次生代謝產(chǎn)物的生物合成過程中起著重要作用[52-54]。本研究通過系統(tǒng)發(fā)育進(jìn)化樹分析、基因表達(dá)模式和分子生物學(xué)試驗(yàn)對茶樹兒茶素生物合成相關(guān)調(diào)控因子CsTT2的功能進(jìn)行初步鑒定。

    在系統(tǒng)發(fā)育進(jìn)化樹中,同一分支的基因往往具有相似的功能。本研究通過系統(tǒng)發(fā)育進(jìn)化樹分析發(fā)現(xiàn),CsTT2與擬南芥中S2、S3、S5、S6、S7和S15中的R2R3-MYB蛋白聚集在同一亞組中(圖2)。S3、S5、S6和S7中的成員均參與調(diào)控?cái)M南芥類黃酮次生代謝物的生物合成[47-51]。此外,S2中的成員AtMYB15能夠控制以及其他基因的表達(dá)去響應(yīng)低溫脅迫[55],S15中的成員參與調(diào)控?cái)M南芥毛狀體以及根毛的形成[56]。這表明CsTT2不僅可能參與調(diào)控茶樹類黃酮次生代謝生物合成,還可能響應(yīng)低溫脅迫以及參與調(diào)控茶樹毛狀體和根毛的形成?;虮磉_(dá)模式分析有利于挖掘該基因的潛在功能。本研究通過高效液相色譜法和qRT-PCR技術(shù)測定茶樹品種頂芽的兒茶素含量以及兒茶素生物合成相關(guān)基因的表達(dá)量。結(jié)果顯示,在兒茶素含量高的品種中,目的基因和兒茶素生物合成相關(guān)基因(、、、和)的表達(dá)量均上調(diào)(圖4)。

    注:A. CsTT2和兒茶素生物合成相關(guān)基因在茶樹8個(gè)不同組織部位的轉(zhuǎn)錄組數(shù)據(jù)分析。熱圖右側(cè)顯示每個(gè)基因的ID號。B. CsTT2和兒茶素生物合成相關(guān)基因在LTDC的8個(gè)不同組織部位的qRT-PCR分析。數(shù)據(jù)是用平均值±標(biāo)準(zhǔn)誤來表示,圖中不同的字母表示顯著不同的值(P<0.05,Tukey’s test)。RT:根,ST:莖,F(xiàn)L:花,F(xiàn)R:果實(shí),AB:頂芽,YL:嫩葉,ML:成熟葉,OL:老葉

    注:A. 兒茶素生物合成相關(guān)基因ANR啟動(dòng)子MYB轉(zhuǎn)錄因子結(jié)合元件的預(yù)測。紅色的三角符號表示MYB轉(zhuǎn)錄因子特異結(jié)合元件。B. CsTT2結(jié)合兒茶素生物合成相關(guān)基因ANR啟動(dòng)子的鑒定。ProANR1:–1?074?bp至–1?423?bp,ProANR2:–724?bp至–1?073?bp,ProANR3:–374?bp至–723?bp,ProANR4:–1?bp至–373?bp。SD/-Leu表示SD培養(yǎng)基缺乏亮氨酸,SD/-Leu+AbA表示SD/-Leu添加AbA,AbA使用濃度:ProANR1和ProANR3: 1?000?ng·mL-1,ProANR2:850?ng·mL-1,ProANR4:900?ng·mL-1;藍(lán)色三角符號代表菌液稀釋的倍數(shù)。C. CsTT2調(diào)控兒茶素生物合成相關(guān)基因(ANR、LAR和SCPL)啟動(dòng)子活性的分析。數(shù)據(jù)是用平均值±標(biāo)準(zhǔn)誤來表示,采用Student’s test與對照(pEAQ+ProANR/ProLAR/ProSCPL)進(jìn)行顯著性分析(**,P<0.01)

    在茶樹中,基因具有多個(gè)拷貝。其中,在煙草中過表達(dá),轉(zhuǎn)基因植株的花青素含量得到積累,而MYB5b能夠調(diào)控兒茶素單體C和EC的合成。在煙草中過表達(dá)(),煙草中類黃酮物質(zhì)生物合成相關(guān)基因(、、和)的表達(dá)量上調(diào),但是轉(zhuǎn)基因株系中的花青素含量并沒有明顯變化[28,57-58],它如何影響相關(guān)基因的表達(dá)以及是否參與調(diào)控茶樹中兒茶素生物合成還有待進(jìn)一步探究。因此,為了進(jìn)一步研究CsTT2的功能,本研究進(jìn)行了亞細(xì)胞定位、酵母自激活、酵母單雜交以及雙熒光素酶報(bào)告系統(tǒng)試驗(yàn)。結(jié)果表明,定位在細(xì)胞核中,是具有轉(zhuǎn)錄激活能力的調(diào)控因子(圖5)。進(jìn)一步試驗(yàn)分析發(fā)現(xiàn),CsTT2可以結(jié)合兒茶素生物合成重要基因的啟動(dòng)子,從而激活它的表達(dá)(圖7)。CsTT2有可能是通過直接激活基因的表達(dá),進(jìn)而調(diào)控茶樹兒茶素的生物合成。若后續(xù)能在茶樹本體上直接進(jìn)行基因功能驗(yàn)證,則更能進(jìn)一步鑒定CsTT2的調(diào)控功能。

    本研究通過基本的生物信息學(xué)分析和一系列分子生物學(xué)試驗(yàn)對CsTT2的功能進(jìn)行初步鑒定。結(jié)果表明,CsTT2可以結(jié)合茶樹兒茶素生物合成重要合成基因的啟動(dòng)子,并激活其表達(dá)。本研究為研究茶樹兒茶素生物合成的分子調(diào)控機(jī)制奠定理論基礎(chǔ)。

    [1] Ming T, Bartholomew B. Theaceae [J]. Flora China, 2007, 12: 366-478.

    [2] Zhao M, Zhang N, Gao T, et al. Sesquiterpene glucosylation mediated by glucosyltransferase UGT91Q2 is involved in the modulation of cold stress tolerance in tea plants [J]. New Phytol, 2020, 226(2): 362-372.

    [3] 宛曉春. 茶葉生物化學(xué)[M]. 3版. 北京: 中國農(nóng)業(yè)出版社, 2003.

    Wan X C. Biochemistry of tea: the third edition [M]. 3rd ed. Beijing: China Agriculture press, 2003.

    [4] Zhao M, Yu Y, Sun L, et al. GCG inhibits SARS-CoV-2 replication by disrupting the liquid phase condensation of its nucleocapsid protein [J]. Nature Communications, 2021, 12(1): 2114.doi: 10.1038/s41467-021-22297-8.

    [5] Xiong L G, Chen Y J, Tong J W, et al. Epigallocatechin-3-gallate promotes healthy lifespan through mitohormesis during early-to-mid adulthood in[J]. Redox Biology, 2018, 14: 305-315.

    [6] Yuan H, Li Y, Ling F, et al. The phytochemical epigallocatechin gallate prolongs the lifespan by improving lipid metabolism, reducing inflammation and oxidative stress in high-fat diet-fed obese rats [J]. Aging cell, 2020, 19(9): e13199. doi: 10.1111/acel.13199.

    [7] Lwxa B, Shang C, Tsza B, et al. Green tea derivative epigallocatechin-3-gallate (EGCG) confers protection against ionizing radiation-induced intestinal epithelial cell death bothand[J]. Free Radical Biology and Medicine, 2020, 161: 175-186.

    [8] Zhang Z, Zhang X, Bi K, et al. Potential protective mechanisms of green tea polyphenol EGCG against COVID-19 [J]. Trends in Food Science & Technology, 2021, 114: 11-24.

    [9] Zhao J, Blayney A, Liu X, et al. EGCG binds intrinsically disordered N-terminal domain of p53 and disrupts p53-MDM2 interaction [J]. Nature Communications, 2021, 12(1): 986. doi: 10.1038/s41467-021-21258-5.

    [10] Bernier L P, York E M, Kamyabi A, et al. Microglial metabolic flexibility supports immune surveillance of the brain parenchyma [J]. Nat Commun, 2020, 11(1): 1559. doi: 10.1038/s41467-020-15267-z.

    [11] Liu Z S, Cai H, Xue W, et al. G3BP1 promotes DNA binding and activation of cGAS [J]. Nat Immunol, 2019, 20(1): 18-28.

    [12] Yang C S, Hong J. Prevention of chronic diseases by tea: possible mechanisms and human relevance [J]. Annual Review of Nutrition, 2013, 33: 161-181.

    [13] 夏濤, 高麗萍. 類黃酮及茶兒茶素生物合成途徑及其調(diào)控研究進(jìn)展[J]. 中國農(nóng)業(yè)科學(xué), 2009, 42(8): 2899-2908.

    Xia T, Gao L P. Research progress on biosynthesis pathway and regulation of flavonoids and catechins [J]. Scientia Agricultura Sinica, 2009, 42(8): 2899-2908.

    [14] Weisshaar B, Jenkins G I. Phenylpropanoid biosynthesis and its regulation [J]. Current Opinion in Plant Biology, 1998, 1(3): 251-257.

    [15] Furukawa T, Eshima A, Koiya M, et al. Coordinate expression of genes involved in catechin biosynthesis incells [J]. Plant Cell Reports, 2002, 21(4): 385-389.

    [16] 夏濤, 高麗萍, 劉亞軍, 等. 茶樹酯型兒茶素生物合成及水解途徑研究進(jìn)展[J]. 中國農(nóng)業(yè)科學(xué), 2013, 46(11): 2307-2320.

    Xia T, Gao L P, Liu Y J, et al. Advances in biosynthesis and hydrolysis of catechins from tea tree[J]. Scientia Agricultura Sinica, 2013, 46(11): 2307-2320.

    [17] 宛曉春. 茶樹次生代謝[M]. 北京: 科學(xué)出版社, 2015.

    Wan X C. Secondary metabolism of tea plant [M]. Beijing: Science Press, 2015.

    [18] 陸建良, 林晨, 駱穎穎, 等. 茶樹重要功能基因克隆研究進(jìn)展[J]. 茶葉科學(xué), 2007, 27(2): 95-103.

    Lu J L, Lin C, Luo Y Y, et al. Advances in cloning of important functional genes from[J]. Journal of Tea Science, 2007, 27(2): 95-103.

    [19] Stafford H A. Flavonoid metabolism pathway to proanthocyanindins (condensed tannins), flavan-3-ols, and unsubstituted flavans [M]. New York: CRC Press, 1990.

    [20] Stafford H A, Lester H H. The conversion of ()-dihydromyricetin to its flavan-3, 4-diol (leucodelphinidin) and to ()-gallocatechin by reductase extracted from tissue cultures ofand[J]. Plant Physiology, 1985, 78(4): 791-794.

    [21] Punyasiri P, Abeysinghe I, Kumar V, et al. Flavonoid biosynthesis in the tea plant: properties of enzymes of the prominent epicatechin and catechin pathways [J]. Archives of Biochemistry & Biophysics, 2004, 431(1): 22-30.

    [22] Niemetz R, Gross G G. Enzymology of gallotannin and ellagitannin biosynthesis [J]. Phytochemistry, 2005, 66(17): 2001-2011.

    [23] Gross G G. From lignins to tannins: forty years of enzyme studies on the biosynthesis of phenolic compounds [J]. Phytochemistry, 2008, 69(18): 3018-3031.

    [24] Liu Y, Gao L, Liu L, et al. Purification and characterization of a novel galloyltransferase involved in catechin galloylation in the tea plant () [J]. Journal of Biological Chemistry, 2012, 287(53): 44406-44417.

    [25] Zhong K, Zhao S Y, J?nsson L J, et al. Enzymatic conversion of epigallocatechin gallate to epigallocatechin with an inducible hydrolase from[J]. Biocatalysis, 2009, 26(4): 306-312.

    [26] Wei C, Hua Y, Wang S, et al. Draft genome sequence ofvar.provides insights into the evolution of the tea genome and tea quality [J]. Proc Natl Acad Sci U S A, 2018, 115(18): E4151-E4158.

    [27] Luo Y, Yu S, Li J, et al. Molecular characterization of WRKY transcription factors that act as negative regulators of O-Methylated catechin biosynthesis in tea plants (L.) [J]. J Agric Food Chem, 2018, 66(43): 11234-11243.

    [28] Wang P, Zhang L, Jiang X, et al. Evolutionary and functional characterization of leucoanthocyanidin reductases from[J]. Planta, 2018, 247(1): 139-154.

    [29] 牛義嶺, 姜秀明. 植物轉(zhuǎn)錄因子MYB基因家族的研究進(jìn)展[J]. 分子植物育種, 2016, 14(8): 2050-2059.

    Niu Y L, Jiang X M. Research progress of plant transcription factor MYB gene family [J]. Molecular Plant Breeding, 2016, 14(8): 2050-2059.

    [30] Martin C, Paz-Ares J. MYB transcription factors in plants [J]. Trends in Genetics, 1997, 13(2): 67-73.

    [31] Verdonk J C, Haring M A, Tunen A J, et al. ODORANT1 regulates fragrance biosynthesis in[J]. Plant Cell, 2005, 17(5): 1612-1624.

    [32] Bomal C, Bedon F, Caron S, et al. Involvement ofMYB1 and MYB8 in phenylpropanoid metabolism and secondary cell wall biogenesis: a comparative in planta analysis [J]. Journal of Experimental Botany, 2008, 59(14): 3925-3939.

    [33] Schaart J G, Dubos C, Romero De La Fuente I, et al. Identification and characterization of MYB-bHLH-WD40 regulatory complexes controlling proanthocyanidin biosynthesis in strawberry [J]. The New Phytologist, 2013, 197(2): 454-467.

    [34] An X H, Tian Y, Chen K Q, et al. MdMYB9 and MdMYB11 are involved in the regulation of the JA-induced biosynthesis of anthocyanin and proanthocyanidin in apples [J]. Plant Cell Physiol, 2015, 56(4): 650-662.

    [35] Tian J, Zhang J, Han Z Y, et al. McMYB12 transcription factors co-regulate proanthocyanidin and anthocyanin biosynthesis in[J]. Scientific Reports, 2017, 7(1): 43715. doi: 10.1038/srep43715.

    [36] James A M, Ma D, Mellway R, et al. Poplar MYB115 and MYB134 transcription factors regulate proanthocyanidin synthesis and structure [J]. Plant Physiology, 2017, 174(1): 154-171.

    [37] Wang N, Qu C, Jiang S, et al. The proanthocyanidin-specific transcription factor MdMYBPA1 initiates anthocyanin synthesis under low temperature conditions in red-fleshed apple [J]. The Plant J, 2018, 96(1): 39-55.

    [38] Xu W, Dubos C, Lepiniec L. Transcriptional control of flavonoid biosynthesis by MYB-bHLH-WDR complexes [J]. Trends in Plant Science, 2015, 20(3): 176-185.

    [39] Terrier N, Torregrosa L, Ageorges A, et al. Ectopic expression of VvMybPA2 promotes proanthocyanidin biosynthesis in grapevine and suggests additional targets in the pathway [J]. Plant Physiol, 2009, 149(2): 1028-1041.

    [40] Gesell A, Yoshida K, Tran L T, et al. Characterization of an apple TT2-type R2R3 MYB transcription factor functionally similar to the poplar proanthocyanidin regulator PtMYB134 [J]. Planta, 2014, 240(3): 497-511.

    [41] Mellway R D, Tran L T, Prouse M B et al. The wound-, pathogen-, and ultraviolet B-responsivegene encodes an R2R3 MYB transcription factor that regulates proanthocyanidin synthesis in poplar [J]. Plant Physiology, 2009, 150(2): 924-941.

    [42] Stracke R, Werber M, Weisshaar B. The R2R3-MYB gene family in[J]. Curr Opin Plant Biol, 2001, 4(5): 447-456.

    [43] Liu R, Wang Y, Tang S, et al. Genome-wide identification of the tea plant bHLH transcription factor family and discovery of candidate regulators of trichome formation [J]. Sci Rep, 2021, 11(1): 10764. doi: 10.21203/rs.3.rs-148784/v1.

    [44] Livak K J, Schmittgen T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta DeltaC(T)) method [J]. Methods, 2001, 25(4): 402-408.

    [45] Chen C, Chen H, Zhang Y, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data [J]. Mol Plant, 2020, 13(8): 1194-1202.

    [46] Dubos C, Stracke R, Grotewold E, et al. MYB transcription factors in[J]. Trends Plant Sci, 2010, 15(10): 573-581.

    [47] Stracke R, Ishihara H, Huep G, et al. Differential regulation of closely related R2R3-MYB transcription factors controls flavonol accumulation in different parts of theseedling [J]. Plant J, 2007, 50(4): 660-677.

    [48] Gonzalez A, Zhao M, Leavitt J M, et al. Regulation of the anthocyanin biosynthetic pathway by the TTG1/bHLH/Myb transcriptional complex inseedlings[J]. Plant J, 2008, 53(5): 814-827.

    [49] Lepiniec L, Debeaujon I, Routaboul J M, et al. Genetics and biochemistry of seed flavonoids [J]. Annu Rev Plant Biol, 2006, 57: 405-430.

    [50] Zhong R, Lee C, Zhou J, et al. A battery of transcription factors involved in the regulation of secondary cell wall biosynthesis in[J]. Plant Cell, 2008, 20(10): 2763-2782.

    [51] Zhou J, Lee C, Zhong R, et al. MYB58 and MYB63 are transcriptional activators of the lignin biosynthetic pathway during secondary cell wall formation in[J]. Plant Cell, 2009, 21(1): 248-266.

    [52] Sun B, Zhu Z, Cao P, et al. Purple foliage coloration in tea (L.) arises from activation of the R2R3-MYB transcription factor CsAN1 [J]. Sci Rep, 2016, 6: 32534. doi: 10.1038/srep32534.

    [53] Wang X C, Wu J, Guan M L, et al.MYB4 plays dual roles in flavonoid biosynthesis [J]. Plant J, 2020, 101(3): 637-652.

    [54] Ma D, Constabel C P. MYB repressors as regulators of phenylpropanoid metabolism in plants [J]. Trends Plant Sci, 2019, 24(3): 275-289.

    [55] Agarwal M, Hao Y, Kapoor A, et al. A R2R3 type MYB transcription factor is involved in the cold regulation of CBF genes and in acquired freezing tolerance [J]. J Biol Chem, 2006, 281(49): 37636-37645.

    [56] Kang Y H, Kirik V, Hulskamp M, et al. Thegene provides a positive feedback loop for cell fate specification in theroot epidermis [J]. Plant Cell, 2009, 21(4): 1080-1094.

    [57] Jiang X, Huang K, Zheng G, et al. CsMYB5a and CsMYB5e fromdifferentially regulate anthocyanin and proanthocyanidin biosynthesis [J]. Plant Sci, 2018, 270: 209-220.

    [58] Wang P, Ma G, Zhang L, et al. A sucrose-induced MYB (SIMYB) transcription factor promoting proanthocyanidin accumulation in the tea plant () [J]. J Agric Food Chem, 2019, 67(5): 1418-1428.

    Expression Analysis and Functional Identification ofTranscription Factor in Tea Plants

    WANG Yuyuan1, LIU Renjian1, LIU Shaoqun1, SHU Canwei2, SUN Binmei1, ZHENG Peng1*

    1. College of Horticulture, South China Agricultural University, Guangzhou 510642, China; 2. College of Agriculture, South China Agricultural University/Guangdong Province Key Laboratory of Microbial Signals and Disease Control/Integrative Microbiology Research Center, Guangzhou 510642, China

    Catechin is one of the most characteristic secondary metabolites in tea plants and the main component affecting tea quality and flavor. It possesses rich pharmacological effects, such as anti-oxidation, anti-virus, lipid lowering and weight loss, etc. In this study, the function of a catechin biosynthetic regulator CsTT2 was preliminarily identified using phylogenetic analysis, gene expression pattern analysis and molecular biology experiments. The results show that CsTT2 was a R2R3-MYB transcription factor that shares a branch withthe MYB transcription factor inthat regulates secondary metabolites. The expression levels ofand catechin biosynthesis genes were relatively high in the apical bud tissue of tea plants with higher total catechin content. The results of subcellular localization, yeast assay and dual luciferase reporting system further reveal thatwas located into the nucleus and the protein it encodes possessed transcriptional activation capacity. CsTT2 could bind to the promoter of a key catechin biosynthetic geneto activate its expression.

    tea plant, R2R3-MYB transcription factor,, catechin biosynthesis

    S571.1

    A

    1000-369X(2022)04-463-14

    2021-12-27

    2022-01-19

    廣州市科技計(jì)劃項(xiàng)目(202102020290)、廣東省自然科學(xué)基金項(xiàng)目(2018A030313089、2021A1515012091)

    王玉源,女,碩士研究生,主要從事茶樹遺傳育種與分子生物學(xué),wangyy@stu.scau.edu.cn。*通信作者:zhengp@scau.edu.cn

    猜你喜歡
    生物
    生物多樣性
    生物多樣性
    上上生物
    發(fā)現(xiàn)不明生物
    史上“最黑暗”的生物
    軍事文摘(2020年20期)2020-11-28 11:42:50
    第12話 完美生物
    航空世界(2020年10期)2020-01-19 14:36:20
    最初的生物
    自然生物被直銷
    清晨生物初歷直銷
    生物的多樣性
    精品亚洲成国产av| 久久久久国产一级毛片高清牌| 日韩成人av中文字幕在线观看| 国产片特级美女逼逼视频| 99久久人妻综合| 深夜精品福利| 成年av动漫网址| 日本欧美视频一区| 亚洲欧美精品自产自拍| 日日撸夜夜添| 国产99久久九九免费精品| 韩国av在线不卡| 巨乳人妻的诱惑在线观看| 国产精品一区二区在线观看99| 亚洲天堂av无毛| 搡老乐熟女国产| 五月天丁香电影| 十八禁高潮呻吟视频| 日韩精品免费视频一区二区三区| 午夜福利网站1000一区二区三区| av网站在线播放免费| 日韩免费高清中文字幕av| 久久人人97超碰香蕉20202| 久久av网站| av.在线天堂| 精品一区二区免费观看| 毛片一级片免费看久久久久| 丝瓜视频免费看黄片| 久久鲁丝午夜福利片| 街头女战士在线观看网站| 精品一区在线观看国产| 免费黄色在线免费观看| 久久天堂一区二区三区四区| 汤姆久久久久久久影院中文字幕| www.熟女人妻精品国产| 美女国产高潮福利片在线看| 午夜福利乱码中文字幕| 久久天躁狠狠躁夜夜2o2o | 男男h啪啪无遮挡| 美女视频免费永久观看网站| 在线精品无人区一区二区三| 日本午夜av视频| 国产成人av激情在线播放| 久久精品久久久久久久性| 国产一卡二卡三卡精品 | 亚洲综合精品二区| 国产精品麻豆人妻色哟哟久久| 叶爱在线成人免费视频播放| 国产成人欧美在线观看 | 亚洲欧美色中文字幕在线| 99精国产麻豆久久婷婷| 国产精品亚洲av一区麻豆 | 亚洲av中文av极速乱| 亚洲欧洲日产国产| 秋霞在线观看毛片| 人人妻人人爽人人添夜夜欢视频| 日本色播在线视频| 国产成人免费无遮挡视频| 亚洲精品国产av蜜桃| 19禁男女啪啪无遮挡网站| 69精品国产乱码久久久| 午夜福利视频精品| 涩涩av久久男人的天堂| 亚洲av日韩在线播放| av又黄又爽大尺度在线免费看| 色播在线永久视频| 久热爱精品视频在线9| 在现免费观看毛片| 五月天丁香电影| 青春草国产在线视频| 自线自在国产av| av视频免费观看在线观看| videosex国产| av网站免费在线观看视频| 看非洲黑人一级黄片| 一级片'在线观看视频| 亚洲精品久久午夜乱码| 不卡av一区二区三区| 国产一级毛片在线| 国产日韩欧美在线精品| 99国产精品免费福利视频| 男女下面插进去视频免费观看| 欧美日韩福利视频一区二区| 高清av免费在线| 国产精品一区二区精品视频观看| 宅男免费午夜| 一区二区三区精品91| 男女边吃奶边做爰视频| 97人妻天天添夜夜摸| 国产xxxxx性猛交| 欧美日韩亚洲综合一区二区三区_| 亚洲国产av新网站| 亚洲欧美一区二区三区久久| tube8黄色片| 香蕉丝袜av| 精品久久久久久电影网| 亚洲欧美一区二区三区久久| 欧美av亚洲av综合av国产av | 丰满少妇做爰视频| 韩国精品一区二区三区| 看免费成人av毛片| 亚洲伊人久久精品综合| 国产在视频线精品| 男人添女人高潮全过程视频| 欧美精品一区二区大全| 国产人伦9x9x在线观看| 建设人人有责人人尽责人人享有的| 人妻 亚洲 视频| 最近最新中文字幕大全免费视频 | 中文天堂在线官网| 亚洲国产av影院在线观看| 激情五月婷婷亚洲| 一级黄片播放器| 国产av一区二区精品久久| 亚洲激情五月婷婷啪啪| 成人18禁高潮啪啪吃奶动态图| 免费黄色在线免费观看| 99国产精品免费福利视频| 成人黄色视频免费在线看| 欧美精品av麻豆av| 国产精品一二三区在线看| 亚洲成人国产一区在线观看 | 欧美人与性动交α欧美精品济南到| 人人澡人人妻人| 建设人人有责人人尽责人人享有的| 尾随美女入室| 波多野结衣一区麻豆| 色视频在线一区二区三区| 18禁观看日本| 你懂的网址亚洲精品在线观看| 一级毛片 在线播放| 无限看片的www在线观看| 91成人精品电影| 人体艺术视频欧美日本| 女人精品久久久久毛片| 老司机在亚洲福利影院| 考比视频在线观看| 波多野结衣一区麻豆| 午夜福利视频精品| 一级爰片在线观看| 欧美精品人与动牲交sv欧美| 亚洲欧洲精品一区二区精品久久久 | 亚洲,一卡二卡三卡| 999久久久国产精品视频| 九草在线视频观看| 精品国产乱码久久久久久小说| 午夜福利视频精品| 国产野战对白在线观看| 纵有疾风起免费观看全集完整版| √禁漫天堂资源中文www| 一级黄片播放器| 久久久欧美国产精品| 电影成人av| 久久久久精品国产欧美久久久 | 高清黄色对白视频在线免费看| 国产激情久久老熟女| 丝袜脚勾引网站| 日韩欧美一区视频在线观看| 日本色播在线视频| 不卡av一区二区三区| 亚洲精品一二三| 中文字幕人妻熟女乱码| 在线 av 中文字幕| 欧美日本中文国产一区发布| 一级片'在线观看视频| 国产亚洲一区二区精品| 亚洲第一区二区三区不卡| av一本久久久久| 女人被躁到高潮嗷嗷叫费观| 满18在线观看网站| 日本爱情动作片www.在线观看| 高清在线视频一区二区三区| 亚洲精品久久久久久婷婷小说| 久久97久久精品| 亚洲第一av免费看| 波多野结衣一区麻豆| 欧美人与性动交α欧美软件| 捣出白浆h1v1| 亚洲欧洲日产国产| 亚洲国产精品国产精品| 男女边摸边吃奶| 国产av国产精品国产| 免费黄频网站在线观看国产| 亚洲av中文av极速乱| 亚洲成人手机| 国产一卡二卡三卡精品 | 在线天堂最新版资源| 久久久久久久国产电影| 丝袜喷水一区| 免费观看av网站的网址| 操出白浆在线播放| 久久久久久久久免费视频了| 国产亚洲av高清不卡| 国产淫语在线视频| 欧美日韩一区二区视频在线观看视频在线| 精品久久久久久电影网| 丁香六月欧美| 亚洲精品国产av蜜桃| 啦啦啦 在线观看视频| 日本av免费视频播放| 亚洲情色 制服丝袜| 亚洲欧美激情在线| a级片在线免费高清观看视频| 一级,二级,三级黄色视频| 国产精品欧美亚洲77777| 黄色视频在线播放观看不卡| 久久精品国产a三级三级三级| 黄色怎么调成土黄色| 9191精品国产免费久久| 欧美成人精品欧美一级黄| av.在线天堂| 青青草视频在线视频观看| 国产一区有黄有色的免费视频| 久久99精品国语久久久| av不卡在线播放| 男女床上黄色一级片免费看| 成人三级做爰电影| 丝瓜视频免费看黄片| 亚洲精品日韩在线中文字幕| 一本大道久久a久久精品| 熟女av电影| 一级毛片 在线播放| 久久久久久免费高清国产稀缺| 人体艺术视频欧美日本| 日日爽夜夜爽网站| 国产成人午夜福利电影在线观看| 一本色道久久久久久精品综合| 久久久久国产精品人妻一区二区| 午夜免费观看性视频| 欧美激情高清一区二区三区 | av不卡在线播放| 国产淫语在线视频| 天天操日日干夜夜撸| 国产男女内射视频| 曰老女人黄片| 成人亚洲精品一区在线观看| 久久影院123| 丁香六月天网| 观看av在线不卡| 欧美精品亚洲一区二区| 人体艺术视频欧美日本| 91精品国产国语对白视频| 成人午夜精彩视频在线观看| 一本—道久久a久久精品蜜桃钙片| 黄网站色视频无遮挡免费观看| 国产伦人伦偷精品视频| 精品卡一卡二卡四卡免费| 久久天堂一区二区三区四区| 高清av免费在线| 老熟女久久久| 一区在线观看完整版| 国产片特级美女逼逼视频| 色94色欧美一区二区| 青春草亚洲视频在线观看| 免费黄频网站在线观看国产| 国产麻豆69| 麻豆乱淫一区二区| 欧美精品一区二区大全| 精品国产国语对白av| 亚洲,欧美,日韩| 自拍欧美九色日韩亚洲蝌蚪91| 天天躁夜夜躁狠狠躁躁| 国产精品无大码| 九色亚洲精品在线播放| 纵有疾风起免费观看全集完整版| 我要看黄色一级片免费的| 日本黄色日本黄色录像| 亚洲国产精品一区三区| 久久久欧美国产精品| 中文字幕色久视频| 国产精品国产av在线观看| 欧美人与性动交α欧美精品济南到| 国产精品一区二区精品视频观看| 永久免费av网站大全| av卡一久久| 久久亚洲国产成人精品v| 国产成人欧美在线观看 | 女人久久www免费人成看片| 伦理电影免费视频| 赤兔流量卡办理| 成人国语在线视频| 亚洲精品一区蜜桃| 99热国产这里只有精品6| 亚洲欧美激情在线| 这个男人来自地球电影免费观看 | 又大又爽又粗| 亚洲欧美清纯卡通| 国产高清不卡午夜福利| 最近中文字幕高清免费大全6| 成年动漫av网址| 国产成人91sexporn| 老汉色∧v一级毛片| a级毛片在线看网站| 精品久久久精品久久久| 精品国产一区二区三区久久久樱花| 最近中文字幕高清免费大全6| 亚洲精品一二三| 熟女少妇亚洲综合色aaa.| 国产亚洲精品第一综合不卡| 七月丁香在线播放| 狠狠精品人妻久久久久久综合| 亚洲人成电影观看| 少妇人妻精品综合一区二区| 纵有疾风起免费观看全集完整版| 美女国产高潮福利片在线看| 欧美在线一区亚洲| 久久久久精品性色| 满18在线观看网站| 国产片特级美女逼逼视频| 精品国产露脸久久av麻豆| 啦啦啦在线免费观看视频4| 亚洲精品一二三| 国产精品三级大全| 一区福利在线观看| 男女边吃奶边做爰视频| 欧美变态另类bdsm刘玥| 国产亚洲最大av| 国产色婷婷99| 大码成人一级视频| 久久国产精品男人的天堂亚洲| 久久这里只有精品19| 国产 精品1| 国产精品二区激情视频| 日本wwww免费看| 丝袜脚勾引网站| 国产成人精品久久久久久| 一本久久精品| 叶爱在线成人免费视频播放| 午夜福利影视在线免费观看| 色网站视频免费| 99久久综合免费| 女性生殖器流出的白浆| 亚洲 欧美一区二区三区| 日本欧美视频一区| 亚洲精品国产一区二区精华液| 国产野战对白在线观看| 爱豆传媒免费全集在线观看| 波多野结衣av一区二区av| e午夜精品久久久久久久| 久久久久久人妻| av福利片在线| 嫩草影视91久久| 日本一区二区免费在线视频| 亚洲av成人不卡在线观看播放网 | 肉色欧美久久久久久久蜜桃| 午夜福利影视在线免费观看| 制服丝袜香蕉在线| 五月天丁香电影| 亚洲av日韩在线播放| av不卡在线播放| 国产伦人伦偷精品视频| 老司机影院成人| videosex国产| 女人爽到高潮嗷嗷叫在线视频| kizo精华| 一级毛片黄色毛片免费观看视频| 少妇人妻久久综合中文| 一本大道久久a久久精品| 久久精品熟女亚洲av麻豆精品| 青草久久国产| 免费在线观看视频国产中文字幕亚洲 | 两性夫妻黄色片| 日本欧美国产在线视频| 国产探花极品一区二区| 成人国产av品久久久| 极品少妇高潮喷水抽搐| 亚洲人成电影观看| 亚洲欧美清纯卡通| 亚洲精品美女久久av网站| 啦啦啦在线观看免费高清www| 亚洲精品国产区一区二| 亚洲精品aⅴ在线观看| 国产av国产精品国产| 在线天堂中文资源库| 亚洲色图综合在线观看| 国产成人精品在线电影| 精品国产国语对白av| 最近中文字幕2019免费版| 卡戴珊不雅视频在线播放| 日韩熟女老妇一区二区性免费视频| 国产av国产精品国产| av有码第一页| 美国免费a级毛片| 国产av精品麻豆| 午夜日韩欧美国产| 男女之事视频高清在线观看 | 国产亚洲一区二区精品| 99香蕉大伊视频| 国产爽快片一区二区三区| 少妇人妻 视频| 啦啦啦啦在线视频资源| 日本av免费视频播放| 999精品在线视频| 成人黄色视频免费在线看| 晚上一个人看的免费电影| 久久精品久久久久久噜噜老黄| 嫩草影视91久久| 日本vs欧美在线观看视频| 在线观看人妻少妇| 日本猛色少妇xxxxx猛交久久| 国产av一区二区精品久久| 中文字幕最新亚洲高清| 深夜精品福利| 精品少妇内射三级| 制服丝袜香蕉在线| 亚洲精品第二区| av在线观看视频网站免费| 蜜桃国产av成人99| 国产一区二区在线观看av| 久久久国产欧美日韩av| 国产成人啪精品午夜网站| 80岁老熟妇乱子伦牲交| 色视频在线一区二区三区| 国产一区二区激情短视频 | 看免费av毛片| 亚洲久久久国产精品| svipshipincom国产片| 国产野战对白在线观看| 久久99精品国语久久久| 欧美日韩国产mv在线观看视频| 亚洲精品一二三| 各种免费的搞黄视频| 欧美日韩亚洲高清精品| 男女边摸边吃奶| av女优亚洲男人天堂| 国产高清不卡午夜福利| 久久99热这里只频精品6学生| kizo精华| 国产亚洲午夜精品一区二区久久| 日本91视频免费播放| 视频区图区小说| 一区二区三区激情视频| 一级片免费观看大全| 夜夜骑夜夜射夜夜干| 午夜福利免费观看在线| 日本爱情动作片www.在线观看| 亚洲精华国产精华液的使用体验| xxx大片免费视频| 亚洲伊人色综图| 日本色播在线视频| 一级片免费观看大全| 久久精品久久精品一区二区三区| 亚洲精品一二三| √禁漫天堂资源中文www| 国产日韩欧美在线精品| 99久久精品国产亚洲精品| 中文字幕亚洲精品专区| 中文欧美无线码| 最近2019中文字幕mv第一页| 久久精品国产亚洲av高清一级| 欧美日韩国产mv在线观看视频| 国产一区二区激情短视频 | 母亲3免费完整高清在线观看| 婷婷色综合www| 亚洲人成77777在线视频| 国产日韩一区二区三区精品不卡| 在线 av 中文字幕| 亚洲五月色婷婷综合| 亚洲av电影在线进入| 亚洲欧洲精品一区二区精品久久久 | 男男h啪啪无遮挡| 宅男免费午夜| 丰满乱子伦码专区| 日韩精品免费视频一区二区三区| 日韩大片免费观看网站| 女人高潮潮喷娇喘18禁视频| 成人国产麻豆网| 亚洲av综合色区一区| 午夜日本视频在线| 9热在线视频观看99| 狠狠婷婷综合久久久久久88av| netflix在线观看网站| www.av在线官网国产| 亚洲情色 制服丝袜| 亚洲欧美成人精品一区二区| 精品少妇一区二区三区视频日本电影 | 国产亚洲欧美精品永久| kizo精华| 好男人视频免费观看在线| 欧美精品一区二区大全| 欧美 日韩 精品 国产| 国产精品亚洲av一区麻豆 | 中文字幕高清在线视频| 国产成人精品久久二区二区91 | 国产精品 国内视频| 老司机靠b影院| 91aial.com中文字幕在线观看| 国产乱来视频区| 亚洲一区二区三区欧美精品| 亚洲一区中文字幕在线| 国产免费福利视频在线观看| 一区二区三区乱码不卡18| 亚洲成人av在线免费| 人妻人人澡人人爽人人| 99精国产麻豆久久婷婷| 人妻人人澡人人爽人人| 18在线观看网站| 国产深夜福利视频在线观看| √禁漫天堂资源中文www| 欧美日本中文国产一区发布| 国产精品熟女久久久久浪| 国产精品偷伦视频观看了| 国产乱人偷精品视频| 亚洲精品国产一区二区精华液| 中文字幕人妻丝袜一区二区 | 视频在线观看一区二区三区| 国产爽快片一区二区三区| 久久久国产一区二区| 精品少妇黑人巨大在线播放| 亚洲四区av| 两性夫妻黄色片| 日本91视频免费播放| 中国国产av一级| 成人亚洲欧美一区二区av| 97精品久久久久久久久久精品| 国产国语露脸激情在线看| 国产又爽黄色视频| 最近中文字幕高清免费大全6| 成人免费观看视频高清| 丁香六月欧美| 两个人免费观看高清视频| 日韩熟女老妇一区二区性免费视频| 日韩成人av中文字幕在线观看| 亚洲男人天堂网一区| 黄色怎么调成土黄色| 在线精品无人区一区二区三| 女人爽到高潮嗷嗷叫在线视频| 亚洲av中文av极速乱| 亚洲图色成人| 国产人伦9x9x在线观看| 99精国产麻豆久久婷婷| 黄网站色视频无遮挡免费观看| e午夜精品久久久久久久| 亚洲欧美成人精品一区二区| av电影中文网址| 18在线观看网站| 极品少妇高潮喷水抽搐| 欧美久久黑人一区二区| 成年美女黄网站色视频大全免费| 蜜桃国产av成人99| 大片电影免费在线观看免费| 日本黄色日本黄色录像| 丝袜脚勾引网站| 视频区图区小说| 考比视频在线观看| 欧美日韩国产mv在线观看视频| 国产极品天堂在线| 97在线人人人人妻| 中文字幕最新亚洲高清| 你懂的网址亚洲精品在线观看| 热re99久久精品国产66热6| 成年女人毛片免费观看观看9 | 欧美人与善性xxx| 欧美国产精品va在线观看不卡| 免费观看a级毛片全部| av在线老鸭窝| av国产精品久久久久影院| 在线看a的网站| 18禁观看日本| 国产人伦9x9x在线观看| 中文字幕另类日韩欧美亚洲嫩草| 久久精品久久久久久久性| 亚洲av福利一区| 五月开心婷婷网| 亚洲精品第二区| 日韩大片免费观看网站| 国产一区二区三区综合在线观看| a 毛片基地| 性少妇av在线| 男女床上黄色一级片免费看| 亚洲精品日本国产第一区| 日韩精品有码人妻一区| 少妇被粗大的猛进出69影院| 亚洲综合精品二区| 亚洲一卡2卡3卡4卡5卡精品中文| 精品人妻在线不人妻| 大片免费播放器 马上看| 大香蕉久久成人网| 国产又爽黄色视频| 美女福利国产在线| 亚洲精品,欧美精品| 国产色婷婷99| 大片免费播放器 马上看| 91国产中文字幕| 男女无遮挡免费网站观看| 午夜日韩欧美国产| 街头女战士在线观看网站| 国产精品香港三级国产av潘金莲 | 极品少妇高潮喷水抽搐| 久久久久久久精品精品| av网站免费在线观看视频| 十八禁高潮呻吟视频| 亚洲精品久久午夜乱码| 一级,二级,三级黄色视频| 亚洲成人免费av在线播放| av在线播放精品| 这个男人来自地球电影免费观看 | 国产精品女同一区二区软件| 亚洲精品aⅴ在线观看| 伊人久久国产一区二区| 精品少妇久久久久久888优播| 久久免费观看电影| av在线老鸭窝| 午夜福利视频精品| 国产精品偷伦视频观看了| 国产精品香港三级国产av潘金莲 | tube8黄色片| 国产免费视频播放在线视频| 免费人妻精品一区二区三区视频| 精品国产一区二区三区四区第35| 少妇猛男粗大的猛烈进出视频| 一区二区三区精品91| avwww免费| 久久精品国产综合久久久| 日本午夜av视频|