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

    茶樹谷丙轉(zhuǎn)氨酶基因的克隆及其表達(dá)分析

    2016-02-06 03:25:10劉志薇吳致君王文麗
    西北植物學(xué)報(bào) 2016年12期

    崔 新,劉志薇,吳致君,李 輝,王文麗,莊 靜

    (南京農(nóng)業(yè)大學(xué) 園藝學(xué)院,茶葉科學(xué)研究所,南京 210095)

    ?

    茶樹谷丙轉(zhuǎn)氨酶基因的克隆及其表達(dá)分析

    崔 新,劉志薇,吳致君,李 輝,王文麗,莊 靜*

    (南京農(nóng)業(yè)大學(xué) 園藝學(xué)院,茶葉科學(xué)研究所,南京 210095)

    該研究基于茶樹的轉(zhuǎn)錄組數(shù)據(jù),采用RT-PCR方法從茶樹‘黃金芽’cDNA中克隆獲得茶樹谷丙轉(zhuǎn)氨酶基因(CsAlaAT),利用熒光定量PCR方法,對(duì)CsAlaAT在茶樹材料‘迎霜’和‘黃金芽’不同組織、溫度脅迫與激素處理的表達(dá)進(jìn)行分析。結(jié)果顯示:CsAlaAT基因開放閱讀框長(zhǎng)度為1 662 bp,編碼553個(gè)氨基酸,含有天冬氨酸轉(zhuǎn)氨酶家族(Aspartate aminotransferase family)典型的AAT-like保守結(jié)構(gòu)域。多序列比對(duì)顯示,該序列與多個(gè)相關(guān)物種的序列一致性達(dá)78.83%,與磷酸吡哆醛(PLP)結(jié)合的10個(gè)氨基酸殘基以及第358位賴氨酸催化位點(diǎn)在物種間高度保守。茶樹CsAlaAT蛋白屬親水性蛋白,相對(duì)分子質(zhì)量為60 877.5 D,等電點(diǎn)為6.11,堿性、酸性、脂肪族和芳香族氨基酸比例分別為12%、11%、22%和8%,無序化特征不明顯,與大麥HvAlaAT具有相似的三維結(jié)構(gòu)。實(shí)時(shí)定量PCR分析表明,CsAlaAT在茶樹‘迎霜’和‘黃金芽’中的表達(dá)均具有組織特異性,且均在根部的表達(dá)量最高;CsAlaAT響應(yīng)高溫(38 ℃)和低溫(4 ℃)脅迫的表達(dá)上調(diào);外源施用脫落酸(ABA)和赤霉素(GA)能夠抑制茶樹中CsAlaAT基因的表達(dá)。

    茶樹;谷丙轉(zhuǎn)氨酶;多序列比對(duì);溫度脅迫;激素;表達(dá)分析

    茶是世界三大無酒精飲料之一,茶樹〔Camelliasinensis(Linn.) O.Kuntze.〕為山茶科(Theaceae)山茶屬(CamelliaLinn.)多年生常綠木本植物[1-2]。作為一種重要的商業(yè)飲料作物,茶樹在中國(guó)、肯尼亞、斯里蘭卡、印度等地均有種植,栽培歷史悠久[3-4]。茶葉中含有許多對(duì)人體健康有益的次級(jí)代謝物質(zhì),包括茶氨酸、茶多酚、咖啡因和芳香油等[5]。其中茶氨酸(Theanine,N-乙基-γ-L-谷氨酰胺)是茶樹特征性非蛋白質(zhì)氨基酸,由Sakato首次從綠茶葉片中提取出并命名為茶氨酸[6]。

    茶氨酸具有焦糖的香味和類似于味精的鮮爽味,能中和茶湯中的苦澀感,是茶葉品質(zhì)的重要影響因素之一。茶樹氮素代謝途徑中,氮素大部分以茶氨酸的形式進(jìn)行儲(chǔ)存和代謝[7-8]。茶氨酸的合成與分解,與茶樹的呼吸代謝和某些物質(zhì)的代謝有關(guān),并與茶樹碳氮代謝的調(diào)節(jié)和控制有關(guān)[9]。通過同位素標(biāo)記研究證實(shí)茶樹茶氨酸的代謝前體為谷氨酸和乙胺[7],通過L-茶氨酸合成酶(TS)在ATP、mg2+、K+的催化下合成[10]。目前,與茶氨酸生物合成直接相關(guān)的酶已經(jīng)有很多研究報(bào)道,然而其中鮮有對(duì)參與茶氨酸代謝前體合成的谷丙轉(zhuǎn)氨酶的報(bào)道。

    谷丙轉(zhuǎn)氨酶(Alanine aminotransferase,AlaAT)又稱為丙氨酸轉(zhuǎn)移酶,是一種依賴磷酸吡哆醛(PLP)作為輔酶的轉(zhuǎn)氨酶,廣泛分布于動(dòng)物、植物、真菌和一些細(xì)菌中[11]。在不同的物種中,谷丙轉(zhuǎn)氨酶可分為2~6個(gè)亞型,通過亞細(xì)胞定位發(fā)現(xiàn)其亞型可分布于細(xì)胞質(zhì)、線粒體和過氧化酶體中[12]。谷丙轉(zhuǎn)氨酶可逆地催化丙酮酸和谷氨酸生成丙氨酸和戊酮二酸[13],在茶樹中這個(gè)可逆反應(yīng)是茶氨酸代謝前體丙氨酸合成的主要途徑,丙氨酸再通過丙氨酸脫羧酶(ADC)脫去羧基形成乙胺用于合成茶氨酸。

    谷丙轉(zhuǎn)氨酶在植物氮素同化、蛋白質(zhì)合成和碳代謝中發(fā)揮重要作用。有研究表明,過表達(dá)OsAlaAT基因的水稻可顯著提高產(chǎn)量[14];Yang等[15]證實(shí)在低氧環(huán)境下,OsAlaAT基因參與調(diào)控水稻種子胚乳中淀粉的合成。谷丙轉(zhuǎn)氨酶也與植物抵御非生物脅迫有關(guān),在擬南芥中缺氧環(huán)境會(huì)誘導(dǎo)AtAlaAT的表達(dá)[16];也有研究指出AlaAT的活性在植物經(jīng)受干旱脅迫后的復(fù)水階段輕微上升,協(xié)助植物恢復(fù)正常代謝[17]。

    本研究從茶樹品種‘黃金芽’(C.sinensisvar. Huangjinya)中克隆獲得茶樹谷丙轉(zhuǎn)氨酶基因(CsAlaAT),并采用生物信息學(xué)方法對(duì)其氨基酸序列進(jìn)行序列比對(duì)等分析。同時(shí)采用熒光定量PCR方法,對(duì)茶樹CsAlaAT基因在2個(gè)茶樹材料‘迎霜’(C.sinensisvar. Yingshuang)和‘黃金芽’不同組織、不同溫度脅迫及不同激素處理的表達(dá)分析,探討了谷丙轉(zhuǎn)氨酶的組織表達(dá)差異與溫度脅迫下的響應(yīng),為后續(xù)茶樹中谷丙轉(zhuǎn)氨酶功能以及茶氨酸代謝的深入研究奠定基礎(chǔ)。

    1 材料和方法

    1.1 試驗(yàn)材料及處理

    供試茶樹材料‘黃金芽’和‘迎霜’種植于南京農(nóng)業(yè)大學(xué)茶葉科學(xué)研究所,為2年生扦插幼苗,于2016年春季選取正常生長(zhǎng)無病害茶樹植株的根、莖、嫩葉和成熟葉進(jìn)行總RNA提取及cDNA合成。對(duì)正常植株分別進(jìn)行溫度脅迫處理與外源激素處理, 4 ℃、38 ℃、0.1 mmol/L ABA、1 mmol/L GA分別處理2、8、和24 h,提取葉片組織總RNA并反轉(zhuǎn)錄成cDNA,以未處理植株作為對(duì)照。另外,以茶樹品種‘黃金芽’葉片的cDNA為模板進(jìn)行茶樹CsAlaAT基因克隆。

    1.2 方 法

    1.2.1 茶樹CsAlaAT基因克隆 茶樹總RNA按照Quick RNA Isolation Kit(北京華越洋生物科技有限公司)試劑盒操作說明提取,提取RNA樣品濃度由微量紫外檢測(cè)儀NanoDrop測(cè)定,RNA質(zhì)量用12 g·L-1變性瓊脂糖凝膠電泳來檢測(cè)。cDNA合成按照 Prime Script RT reagent Kit(大連TaKaRa公司) 試劑盒操作說明完成?;诒緦?shí)驗(yàn)室獲得的茶樹轉(zhuǎn)錄組數(shù)據(jù)[18]設(shè)計(jì)1對(duì)引物(5′-ATGTGGAAATTCGTAGCCGAC-3′和3′- ATCACGAAATTCATCCATGAA-5′)。以‘黃金芽’cDNA第一鏈為模板進(jìn)行PCR擴(kuò)增。擴(kuò)增體系為20 μL體系:ddH2O 7 μL,ExTaqMix 10 μL,上下游引物各1 μL,模板1 μL。反應(yīng)條件為:95 ℃預(yù)變性5 min,95 ℃變性30 s,58 ℃退火30 s,72 ℃延伸1 min 20 s,共35個(gè)循環(huán);72 ℃延伸10 min。PCR產(chǎn)物用12 g·L-1瓊脂糖凝膠電泳分離,切取目標(biāo)膠塊,參照DNA回收試劑盒說明書將其回收純化,連接至pMD19-T載體,轉(zhuǎn)入感受態(tài)大腸桿菌DH5α菌株,挑取陽性克隆菌液送至南京金斯瑞公司測(cè)序。

    1.2.2 序列分析 利用NCBI網(wǎng)站相關(guān)程序(http://blast.ncbi.nlm.nih.gov/Blast.cgi)對(duì)獲得的核苷酸和蛋白序列進(jìn)行BLAST比對(duì)搜索和保守域預(yù)測(cè)。序列多重比對(duì)、組成成分和理化性質(zhì)用DNAMAN 6.0和EXPASY(http://web.expasy.org/protparam/)軟件完成;使用MEGA6[19]完成進(jìn)化樹測(cè)試和編輯并生成報(bào)告圖形;利用Foldindex[20]程序和SOPMA[21]程序完成蛋白質(zhì)無序化特性分析和二級(jí)結(jié)構(gòu)預(yù)測(cè)。采用Swiss-model[22]軟件進(jìn)行蛋白質(zhì)三維分子結(jié)構(gòu)分析。熒光定量PCR的數(shù)據(jù)分析采用IBM SPSS Statistic 20和 Microsoft Excel2010軟件制作完成。

    1.2.3 茶樹CsAlaAT基因的表達(dá)分析 熒光定量PCR按照SYBR PremixExTaq試劑盒(大連TaKaRa公司)操作說明進(jìn)行。根據(jù)本克隆的‘黃金芽’CsAlaAT基因序列設(shè)計(jì)表達(dá)檢測(cè)引物(5′-CGAGTCCTACGAGTCTTATTATGC-3′和3′-GGAGGCGTTGACAATAGAATG-5′)。每個(gè)處理設(shè)置3個(gè)生物學(xué)重復(fù),并分別用茶樹Actin基因和TBP基因作為內(nèi)參,與目標(biāo)基因一起擴(kuò)增[23]。在Bio-CFX96 Real-time PCR System 中完成熒光定量PCR。擴(kuò)增程序?yàn)椋?5 ℃預(yù)變性30 s;95 ℃變性5 s,60 ℃退火30 s,40個(gè)循環(huán),然后利用 Bio-cfx manager進(jìn)行溶解曲線分析,采用2-ΔΔCt法[24]進(jìn)行結(jié)果分析。

    2 結(jié)果與分析

    2.1 茶樹CsALaAT基因克隆

    以‘黃金芽’葉片的cDNA為模板,擴(kuò)增得到長(zhǎng)約1 700 bp片段。序列分析結(jié)果表明,該片段含有1個(gè)長(zhǎng)度為1 662 bp開放閱讀框(ORF),共編碼553個(gè)氨基酸(圖1)。

    2.2 茶樹CsAlaAT序列比對(duì)和進(jìn)化分析

    對(duì)茶樹CsAlaAT基因推導(dǎo)氨基酸的保守域預(yù)測(cè)結(jié)果(圖3,A)顯示,第154~528位序列含有1個(gè)典型AAT-like保守結(jié)構(gòu)域,該序列還包含磷酸吡哆醛(PLP)結(jié)合位點(diǎn)[25]和一個(gè)參與形成席夫堿(Schiffbase)或醛亞胺配體(aldimine intermediate)的催化殘基第358位賴氨酸(Lys)[26],以及形成二聚體多肽結(jié)合位點(diǎn)的8個(gè)氨基酸殘基,以上特征表明CsAlaAT屬于天冬氨酸轉(zhuǎn)氨酶家族。

    此外,將茶樹CsAlaAT同大麥(Hordeumvulgare)、擬南芥(Arabidopsisthaliana)、油菜(Brassicanapus)、煙草(Nicotianatabacum)、水稻(Oryzasativa)、馬鈴薯(Solanumtuberosum)、玉米(Zeamays)的AlaAT氨基酸序列進(jìn)行多重對(duì)比,結(jié)果顯示一致性為78.83%。其中上述物種的AlaAT均含有AAT-like結(jié)構(gòu)域,其中構(gòu)成磷酸吡哆醛(PLP)結(jié)合位點(diǎn)的10個(gè)氨基酸殘基,包括一個(gè)與PLP形成醛亞胺鍵的關(guān)鍵賴氨酸殘基,在上述物種間一致,具有高度的保守性(圖3,B)。

    為了進(jìn)一步分析茶樹CsAlaAT的進(jìn)化關(guān)系,對(duì)其進(jìn)行Blast同源檢索與比對(duì),得到不同物種間與其相似程度較高的AlaAT蛋白氨基酸序列,選取毛果楊(Populustrichocarpa)、可可(Theobromacacao)、甜椒(Capsicumannuum)、煙草(N.tabacum)、蕪菁(Brassicarapa)、咖啡(Coffeacanephora)、番茄(Solanumlycopersicum)、罌粟(Papaversomniferum)、棉花(Gossypiumhirsutum)、甜瓜(Cucumismelo)、馬鈴薯(S.tuberosum)、大豆(Glycinemax)、擬南芥(A.thaliana)、水稻(O.sativa)、大麥(H.vulgare)、油菜(B.napus)、玉米(Z.mays)的AlaAT氨基酸序列構(gòu)建同源進(jìn)化樹。結(jié)果顯示,茶樹CsAlaAT氨基酸序列與咖啡、甜瓜等的進(jìn)化關(guān)系較近,與水稻、玉米等的進(jìn)化關(guān)系較遠(yuǎn)(圖2)。

    2.3 茶樹CsAlaAT氨基酸組成及理化性質(zhì)分析

    根據(jù)茶樹CsAlaAT氨基酸序列在BLAST同源檢索與比對(duì),獲得不同植物中相似度較高的AlaAT氨基酸序列,利用ExPASy-ProtParam對(duì)其進(jìn)行氨基酸序列組成及理化性質(zhì)分析。結(jié)果顯示(表1),在這些植物中AlaAT蛋白殘基數(shù)為481~575。相對(duì)分子質(zhì)量為(5.2~6.3)×104;堿性氨基酸比重略多于酸性氨基酸,等電點(diǎn)范圍為5.3~7.5之間。脂肪族氨基酸占比22%左右,芳香族占比8%;總平均疏水性(GRAVY)均為負(fù)值,表明這些蛋白均屬于親水性蛋白。

    圖1 茶樹CsALaAT核苷酸和推導(dǎo)氨基酸序列Fig.1 Nucleotide acid and deduced amino acid sequences of CsAlaAT from tea plant

    2.4 茶樹CsAlaAT氨基酸序列的無序化分析

    利用FoldIndex程序?qū)sAlaAT氨基酸序列進(jìn)行了折疊的無序化分析。結(jié)果顯示(圖4):序列中無序化區(qū)域共有2個(gè),最大無序化區(qū)域含28個(gè)氨基酸,總的無序化氨基酸數(shù)目為42,無序化比例為7.59%;總體而言CsAlaAT氨基酸序列無序化不明顯。

    2.5 茶樹CsAlaAT的二級(jí)結(jié)構(gòu)和三級(jí)結(jié)構(gòu)預(yù)測(cè)

    從茶樹CsAlaAT全序列的二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果

    可知(圖5),CsAlaAT蛋白由41.41%的α螺旋,16.09%的β-折疊,12.84%的β-轉(zhuǎn)角和26.66%的隨機(jī)卷曲組成。結(jié)果表明,α-螺旋是CsAlaAT的主要組成部分,而β轉(zhuǎn)角和β-折疊則散布于CsAlaAT蛋白序列中。

    圖2 茶樹CsAlaAT與其他物種氨基酸序列的系統(tǒng)進(jìn)化樹Fig.2 The phylogenetic tree of amino acid sequences of the CsAlaAT

    以大麥HvAlaAT(PDB ID:3TCM)為模型,對(duì)CsAlaAT進(jìn)行蛋白質(zhì)三級(jí)結(jié)構(gòu)建模分析,結(jié)果顯示HvAlaAT與CsAlaAT的三級(jí)結(jié)構(gòu)非常相似(圖6),一致性為78.38%。茶樹CsAlaAT在細(xì)胞內(nèi)以二聚體的形式存在,其結(jié)構(gòu)域與HVAlaAT相比均包括蛋白質(zhì)N端(3~63)、一個(gè)較小的αβ-結(jié)構(gòu)域(氨基酸殘基64~107,396~517)和一個(gè)較大的αβα-結(jié)構(gòu)域(氨基酸殘基107~395),活化位點(diǎn)以及PLP加合物(PLP adduct)處于這兩個(gè)結(jié)構(gòu)域中間。

    上劃線部分為AAT-like結(jié)構(gòu)域,下標(biāo)三角形的為PLP結(jié)合位點(diǎn),下標(biāo)矩形的為催化位點(diǎn)圖3 茶樹CsAalAT保守域(A)及與其他物種氨基酸序列的多重對(duì)比(B)Upper score represents the AAT-like conserved domain; triangles represent the PLP binding sites; rectangle represents the catalytic residue respectivelyFig.3 The conserved domain (A) of CsAalAT and alignment of amino acid sequences from tea plant and other different plants (B)

    序列圖中下劃線代表無序狀態(tài)的氨基酸,無下劃線代表有序狀態(tài)的氨基酸,坐標(biāo)圖中正負(fù)數(shù)值分別表示蛋白的有序和無序狀態(tài)圖4 CsAlaAT折疊狀態(tài)的分析Amino acids represented ordered and non-ordered regions are shown in non-underline and underline, respectively; Positive and negative numbers represent ordered and disordered protein, respectivelyFig.4 Prediction of the folding state of CsAlaAT

    植物Plant名稱Name登錄號(hào)GenBankNo.總氨基酸數(shù)目Totalaminoacid/aa相對(duì)分子質(zhì)量Molecularweight/D理論等電點(diǎn)pI各氨基酸比例Compositionofaminoacid/%堿性Basic酸性Acid脂肪族Aliphatic芳香族Aromatics總平均疏水性Grandaverageofhydropathicity茶樹C.sinensisCsAlaAT55360877.56.111211228-0.209擬南芥A.thalianaAtAlaAT2NP_565040.254059510.65.951312218-0.269咖啡C.canephoraCaAlaATCDO97792.148753518.95.441212228-0.213甜瓜C.meloCmAlaAT2XP_008460348.151356307.16.251211228-0.176大豆G.maxGmAlaAT1ABW17196.148653320.65.311112218-0.213甜椒C.a(chǎn)nnuumCaAlaATAAR05449.148152793.25.291112238-0.147煙草N.tabacumNtAlaAT2XP_016473669.154760127.46.011211218-0.233番茄S.lycopersicumSlAlaAT2XP_004235527.157563179.27.541311228-0.184馬鈴薯S.tuberosumStAlaAT2XP_006342897.154559613.75.891211218-0.219水稻O.sativaOsAlaATAAK52114.148452692.36.221111238-0.110油菜B.napusBnAlaATXP_013744571.154059573.06.251312218-0.241玉米Z.maysZmAlaATAAC62456.148252997.76.291211229-0.165大麥H.vulgareHvAlaAT39575927350055135.76.061412228-0.214

    藍(lán)色. α-螺旋;紅色. β-折疊;綠色. β-轉(zhuǎn)角;粉色. 隨機(jī)卷曲圖5 茶樹CsAlaAT的二級(jí)結(jié)構(gòu)預(yù)測(cè)Blue. Alpha helix, Red. Extended strand, Green. Beta turn, Pink. Random coilFig.5 The secondary structure of the deduced CsAlaAT polypeptide

    圖6 茶樹CsAlaAT的三維結(jié)構(gòu)預(yù)測(cè)Fig.6 The three-dimension structures of CsAlaAT from tea plant

    2.6 茶樹CsAlaAT基因在茶樹不同部位的表達(dá)分析

    通過實(shí)時(shí)熒光定量PCR檢測(cè)CsAlaAT基因在茶樹‘迎霜’和‘黃金芽’不同部位的表達(dá)情況。結(jié)果顯示,茶樹CsAlaAT基因在‘迎霜’和‘黃金芽’中的表達(dá)模式較為相似,均在根部表達(dá)量最高(圖7)?!汀S金芽’CsAlaAT基因表達(dá)量由高到低均為根、莖、成熟葉、幼葉;根的表達(dá)量顯著高于其他3個(gè)部位,幼葉和成熟葉表達(dá)量較為接近。此外,CsAlaAT基因在‘黃金芽’根部的表達(dá)量顯著高于在‘迎霜’根部的表達(dá)量,為‘迎霜’根部的26.3倍。

    2.7 CsAlaAT基因不同溫度和激素處理下的表達(dá)情況

    茶樹‘迎霜’在4 ℃低溫、ABA處理和GA處理下,CsAlaAT表達(dá)量較對(duì)照均表現(xiàn)出一定程度下調(diào),基本趨勢(shì)為處理初的2 h左右表達(dá)量最低,隨處理時(shí)間延長(zhǎng)表達(dá)量有一定回升,但仍顯著低于對(duì)照(圖8,A)?!S金芽’與‘迎霜’的CsAlaAT基因在38 ℃高溫處理下,表達(dá)量相對(duì)于對(duì)照均有上調(diào),在處理2、8、24 h后的表達(dá)量分別為對(duì)照的1.7、1.1、

    2.4倍和4.2、2.5、2.1倍(圖8,B)。茶樹品種‘黃金芽’,在4 ℃低溫處理下,表達(dá)量相較對(duì)照有一定上調(diào),在處理2、8、24 h后分別為對(duì)照的1.2、1.7和1.6倍(圖8,C)。在ABA處理下,‘黃金芽’的CsAlaAT基因表達(dá)與‘迎霜’的表達(dá)模式相似,即處理最初下調(diào)最顯著,隨處理時(shí)間延長(zhǎng)開始有一定回升。在GA處理下,‘黃金芽’的CsAlaAT在處理2 與8 h情況下,下調(diào)量并未達(dá)到顯著水平,在持續(xù)處理24 h后,下調(diào)量達(dá)到顯著水平(圖8,D)。

    不同小寫字母表示同一品種不同組織的差異顯著(P<0.05)圖7 CsAlaAT基因在茶樹不同組織的表達(dá)Different letters represent significant difference among tissues of two cultivarsFig.7 Expression analysis of AlaAT in different tissues

    不同小寫字母表示同一品種不同處理的差異顯著(P<0.05)圖8 CsAlaAT基因在不同處理下的表達(dá)Different letters represent significant difference among different treatments of two cultivars (P<0.05)Fig.8 Expression analysis of CsAlaAT under different stress treatments

    3 討 論

    谷丙轉(zhuǎn)氨酶是一類依賴磷酸吡哆醛作為輔酶的轉(zhuǎn)氨酶類,其催化的可逆反應(yīng)中可利用的底物有碳代謝中關(guān)鍵的戊酮二酸、丙酮酸,以及參與氮素同化及信號(hào)轉(zhuǎn)導(dǎo)的谷氨酸。谷丙轉(zhuǎn)氨酶在植物的碳氮代謝有著不可或缺的作用,研究也證實(shí)其可協(xié)助植物抵御脅迫[27-28]。本實(shí)驗(yàn)從茶樹品種‘黃金芽’中克隆出長(zhǎng)約1 700 bp片段,該片段含有1個(gè)典型的AAT-like保守結(jié)構(gòu)域,屬于天冬氨酸轉(zhuǎn)氨酶超家族。與多種親緣關(guān)系相近的植物AlaAT氨基酸序列具有較高的相似性和相近的理化性質(zhì)。推導(dǎo)的三級(jí)結(jié)構(gòu)與大麥HvAlaAT三級(jí)結(jié)構(gòu)一致性較高,以二聚體的形式進(jìn)行催化,參與形成希夫堿的活化位點(diǎn)第358位賴氨酸(Lys)殘基在2個(gè)物種間保持一致。

    ‘黃金芽’和‘迎霜’不同組織和不同逆境處理下表達(dá)分析顯示,茶樹CsAlaAT基因表達(dá)具有組織特異性,均在根部表達(dá)量最高。其中‘黃金芽’根部表達(dá)量顯著地高于‘迎霜’根部表達(dá)量。Rocha等[29]測(cè)量了谷丙轉(zhuǎn)氨酶在大豆(Glycinemax)中的組織特異性表達(dá),GmAlaAT在接種了根瘤菌的大豆植株的根部表達(dá)量最高,高于沒接種根瘤菌的根部表達(dá)量,也高于在葉和豆莢中的表達(dá)量。這個(gè)結(jié)果同CSAlaAT在‘黃金芽’中的表達(dá)模式相似。接種根瘤菌可提高大豆的氮素同化率[30],‘黃金芽’屬白化茶樹品種,具有氨基酸含量高,茶多酚和咖啡因量適中的特點(diǎn)[31],推測(cè)其根部大量的茶氨酸合成和旺盛的氮素同化是CsAlaAT在根部表達(dá)量高的可能原因。

    目前,對(duì)茶氨酸代謝相關(guān)的研究較多,然而對(duì)茶樹谷丙轉(zhuǎn)氨酶研究報(bào)道較少,對(duì)其在茶樹碳氮代謝以及抵御脅迫中的作用和機(jī)理有待深入研究。對(duì)茶樹谷丙轉(zhuǎn)氨酶及其家族成員的研究,可作為一個(gè)新的角度來理解茶樹茶氨酸代謝與碳氮代謝之間的聯(lián)系以及其響應(yīng)茶樹非生物脅迫的分子機(jī)理。

    [1] 童啟慶. 茶樹栽培學(xué)(3版)[M].北京:中國(guó)農(nóng)業(yè)出版社,2007:4-33.

    [2] 王永鑫,劉志薇,吳致君,等. 茶樹中2個(gè)NAC轉(zhuǎn)錄因子基因的克隆及溫度脅迫的響應(yīng)[J]. 西北植物學(xué)報(bào), 2015, 35(11):2 148-2 156.

    WANG Y X, LIU Z W, WU Z J,etal. Isolation and expression profiles analysis of NAC transcription factor genes under temperature stresses inCamelliasinensis[J].ActaBotanicaBoreali-OccidentaliaSinica,2015, 35(11):2 148-2 156.

    [3] CHEN L, ZHOU Z X, YANG Y J. Genetic improvement and breeding of tea plant (Camelliasinensis) in China: from individual selection to hybridization and molecular breeding[J].Euphytica, 2007, 154(1-2): 239-248.

    [4] 張明澤,姚玉仙,陳世軍. 黔南60份茶樹種質(zhì)資源遺傳多樣性的SSR分析[J]. 西北植物學(xué)報(bào), 2016, 36(6):1 117-1 124.

    ZHANG M Z,YAO Y X,CHEN S J. Genetic diversity analysis of tea germplasm in Qiannan prefecture by SSR markers[J].ActaBotanicaBoreali-OccidentaliaSinica, 2016, 36(6):1 117-1 124.

    [5] ROGERS P J, SMITH J E, HEATHERLEY S V,etal. Time for tea: mood, blood pressure and cognitive performance effects of caffeine and theanine administered alone and together[J].Psychopharmacology, 2008, 195(4): 569-577.

    [6] VUONG Q V, BOWYER M C, ROACH P D. L-Theanine: properties, synthesis and isolation from tea[J].JournaloftheScienceofFoodandAgriculture, 2011, 91(11): 1 931-1 939.

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

    [8] 陳 琪. 茶樹體內(nèi)茶氨酸合成酶的克隆與異源表達(dá)及一氧化氮信號(hào)對(duì)其調(diào)控的研究[D]. 合肥:安徽農(nóng)業(yè)大學(xué), 2011.

    [9] KONISHI S, KASAI Z. Metabolism and regulation of theanine and related compounds in tea plants. II. Synthesis of theanine from carbon-14C dioxide in tea plants and sites of the synthesis[J].NipponDojo-HiryogakuZasshi, 1968, 39(9): 439-443.

    [10] SASAOKA K., KITO M., INAGAKI H. Studies on the biosynthesis of theanine in tea seedlings: synthesis of theanine by the homogenate of tea seedlings [J].AgriculturalandBiologicalChemistry, 1963, 27(6):467-468.

    [12] LIEPMAN A H, OLSEN L J. Genomic analysis of aminotransferases inArabidopsisthaliana[J].CriticalReviewsinPlantSciences, 2004, 23(1): 73-89.

    [13] MCALLISTER C H, FACETTE M, HOLT A,etal. Analysis of the enzymatic properties of a broad family of alanine aminotransferases[J].PloSone, 2013, 8(2): e55032.

    [14] BEATTY P H, SHRAWAT A K, CARROLL R T,etal. Transcriptome analysis of nitrogen-efficient rice over-expressing alanine aminotransferase[J].PlantBiotechnologyJournal, 2009,7(6): 562-576.

    [15] YANG J, KIM S R,etal. Alanine aminotransferase 1 (OsAlaAT1) plays an essential role in the regulation of starch storage in rice endosperm[J].PlantScience, 2015, 240: 79-89.

    [16] MIYASHITA Y, DOLFERUS R, ISMOND K P,etal. Alanine aminotransferase catalyses the breakdown of alanine after hypoxia inArabidopsisthaliana[J].ThePlantJournal, 2007, 49(6): 1 108-1 121.

    [17] MUENCH D G, GOOD A G. Hypoxically inducible barley alanine aminotransferase: cDNA cloning and expression analysis[J].PlantBolecularBiology, 1994, 24(3): 417-427.

    [18] WU Z J, LI X H,etal. De novo assembly and transcriptome characterization: novel insights into catechins biosynthesis inCamelliasinensis[J].BMCPlantBiology, 2014, 14(1): 1.

    [19] TAMURA K, STECHER G, PETERSON D,etal. MEGA6: molecular evolutionary genetics analysis version 6.0[J].MolecularBiologyandEvolution, 2013, 30(12): 2 725-2 729.

    [20] PRILUSKY J, FELDER C E, ZEEV-BEN-MORDEHAI T,etal. FoldIndex?: a simple tool to predict whether a given protein sequence is intrinsically unfolded[J].Bioinformatics, 2005, 21(16): 3 435-3 438.

    [21] COMBET C, BLANCHET C, GEOURJON C,etal. NPS@: network protein sequence analysis[J].TrendsinBiochemicalSciences, 2000, 25(3): 147-150.

    [22] SCHWEDE T, KOPP J, GUEX N,etal. SWISS-MODEL: an automated protein homology-modeling server[J].NucleicAcidsResearch, 2003, 31(13): 3 381-3 385.

    [23] WU Z J, LI X H, LIU Z W,etal. Transcriptome-wide identification ofCamelliasinensisWRKY transcription factors in response to temperature stress[J].MolecularGeneticsandGenomics, 2016, 291(1): 255-269.

    [24] PFAFFL M W. A new mathematical model for relative quantification in real-time RT-PCR[J].NucleicAcidsResearch, 2001, 29(9): e45-e45.

    [25] FAN T W M, HIGASHI R M, LANE A N. An in vivo 1 H and 31 P NMR investigation of the effect of nitrate on hypoxic metabolism in maize roots[J].ArchivesofBiochemistryandBiophysics, 1988, 266(2): 592-606.

    [26] CAPITANI G, HOHENESTER E, FENG L,etal. Structure of 1-aminocyclopropane-1-carboxylate synthase, a key enzyme in the biosynthesis of the plant hormone ethylene[J].JournalofMolecularBiology, 1999, 294(3): 745-756.

    [27] SIVARAMAN J, Li Y, LAROCQUE R,etal. Crystal structure of histidinol phosphate aminotransferase (HisC) from Escherichia coli, and its covalent complex with pyridoxal-5′-phosphate and l-histidinolphosphate[J].JournalofMolecularBiology, 2001, 311(4): 761-776.

    [28] DESOUSA C A F, SODEK L. Alanine metabolism and alanine aminotransferase activity in soybean (Glycinemax) during hypoxia of the root system and subsequent return to normoxia[J].EnvironmentalandExperimentalBotany, 2003, 50(1): 1-8.

    [29] ROCHA M, SODEK L, LICAUSI F,etal. Analysis of alanine aminotransferase in various organs of soybean (Glycinemax) and in dependence of different nitrogen fertilisers during hypoxic stress[J].AminoAcids, 2010, 39(4): 1 043-1 053.

    [30] KIM K J, PARK C J, AN J M,etal. CaAlaAT1 catalyzes the alanine: 2-oxoglutarate aminotransferase reaction during the resistance response against Tobacco mosaic virus in hot pepper[J].Planta, 2005, 221(6): 857-867.

    [30] 王 浩,劉 偉,姜 妍,等. 不同氮素水平下接種根瘤菌對(duì)大豆生長(zhǎng)的影響[J]. 大豆科技,2012,1:14-17.

    WANG H, LIU W, JIANG Y,etal. Effect of Rhizobium Inoculation and nitrugen level on growth of soybean[J].SoybeanScienceandTechnology, 2012, 1 :14-17.

    [31] 李娜娜. 新梢白化茶樹生理生化特征及白化分子機(jī)理研究[D]. 杭州: 浙江大學(xué),2015.

    [33] 潘瑞熾,王小菁,李娘輝. 植物生理學(xué)[M].北京:高等教育出版社,2001: 288-290.

    [34] 張桂蓮,陳立云,張順堂,等.高溫脅迫對(duì)水稻劍葉氮代謝的影響[J].雜交水稻,2007,22(4): 57-61.

    ZHANG G L, CHEN L Y, ZHANG S T,etal. Effect of high temperature stress on nitrogen metabolism of flag leaves in Rice[J].HybridRice,2007,22(4): 57-61.

    [35] ZHOLKEVICH V N, ZUBKOVA N K, MAEVSKAYA S N,etal. Interaction between heat shock and water stress in plants: 2. Osmoregulation in the leaves of cotton plants successively treated with short-term hyperthermia and soil drought[J].RussianJournalofPlantPhysiology, 1997, 44: 533-542.

    [36] VANLERBERGHE G C, JOY K W, TURPIN D H. Anaerobic metabolism in the N-limited green alga selenastrumminutum III. Alanine is the product of anaerobic ammonium assimilation[J].PlantPhysiology, 1991, 95(2): 655-658.

    (編輯:宋亞珍)

    Cloning and Expression Analysis of the Gene Encoding Alanine Aminotransferase inCamelliasinensis

    CUI Xin, LIU Zhiwei, WU Zhijun, LI Hui, WANG Wenli, ZHUANG Jing*

    (College of Horticulture, Tea Science Research Institute, Nanjing Agricultural University, Nanjing 210095, China)

    Alanine aminotransferase (AlaAT) is a pyridoxal-5′-phosphate-dependent (PLP) enzyme which plays a critical role in linking carbon and nitrogen metabolism. AlaAT is also involved in plant responses to abiotic stress. In this study, based on the transcriptome data of tea plant, we clonedCsAlaATgene encoding alanine aminotransferase from cDNA of ‘Huangjinya’ by RT-PCR method. In addition, the expression profiles of theCsAlaATgene in different tissues and under extreme temperatures and hormone treatments in two tea plant varieties (‘Huangjinya’ and ‘Yingshuang’) were analyzed by quantitative real time PCR. Results showed that the length of open reading frame (ORF) ofCsAlaATgene was 1 662 bp, encoding 553 amino acids. CsAlaAT contains typical conserved AAT-like domain which belongs to Aspartate aminotransferase family. Multiple alignments of CsAlaAT with related plant species showed that the identity of them was 78.83%. The PLP binding sites and the catalytic residue of 358 Lys were highly conserved. The CsAlaAT is hydrophilic protein. Its theoretical relative molecular weight is 60 877.5 D. Theoretical isoelectric point is 6.11. Percentages of basic, acidic, aliphatic and aromatic amino acids were 12%, 11%, 22% and 8%, respectively. A three-dimension crystal structure of CsAlaAT was established on the basis of the HvAlaAT which share identity of 78.38% with CsAlaAT. Quantitative real-time PCR analysis of the expression profiles showed that theCsAlaATgene was tissue-specific expressed in two tea plant varieties ‘Huangjinya’ and ‘Yingshuang’. The highest expression level was found in the root. TheCsAlaATgene could respond to high temperature (38 ℃) and low temperature (4 ℃) treatments. Exogenous application of abscisic acid (ABA) and gibberellins (GA) could down-regulated theCsAlaAT.

    Camelliasinensis; alanine aminotransferase; multiple alignment; temperature stress; hormone; expression profiles

    1000-4025(2016)12-2361-09

    10.7606/j.issn.1000-4025.2016.12.2361

    2016-10-18;修改稿收到日期:2016-12-09

    國(guó)家自然科學(xué)基金(31570691)

    崔 新(1992-),男,碩士研究生,主要從事茶樹分子生物學(xué)研究。E-mail: 2015104091@njau.edu.cn

    *通信作者:莊 靜,教授,博士生導(dǎo)師,主要從事茶樹遺傳育種與茶樹分子生物學(xué)研究。E-mail: zhuangjing@njau.edu.cn

    Q785;Q786

    A

    黑人欧美特级aaaaaa片| 欧美乱色亚洲激情| 级片在线观看| 欧美久久黑人一区二区| 亚洲国产中文字幕在线视频| 99久久综合精品五月天人人| 中文字幕色久视频| 久久精品国产99精品国产亚洲性色 | 久久久国产精品麻豆| 啦啦啦 在线观看视频| 99热只有精品国产| 黄片小视频在线播放| 亚洲成人国产一区在线观看| 国产精品二区激情视频| 天堂√8在线中文| 日韩三级视频一区二区三区| 国产精品九九99| 91av网站免费观看| 一夜夜www| 99国产精品免费福利视频| 久久久国产欧美日韩av| 久久久久久人人人人人| 亚洲全国av大片| 国产欧美日韩一区二区三区在线| 国产精品综合久久久久久久免费 | 老汉色av国产亚洲站长工具| 日本精品一区二区三区蜜桃| 天堂√8在线中文| 91成年电影在线观看| 淫秽高清视频在线观看| 久久中文字幕一级| 亚洲精品中文字幕一二三四区| 18禁美女被吸乳视频| 人成视频在线观看免费观看| 国产精品久久久久久人妻精品电影| 一级,二级,三级黄色视频| 国产成人影院久久av| 在线观看日韩欧美| 两人在一起打扑克的视频| 欧美中文日本在线观看视频| 成熟少妇高潮喷水视频| 免费久久久久久久精品成人欧美视频| svipshipincom国产片| 久久人妻av系列| 看免费av毛片| 亚洲av成人一区二区三| 久久精品亚洲精品国产色婷小说| 国产精品国产高清国产av| 国产精品亚洲一级av第二区| 少妇粗大呻吟视频| 99国产极品粉嫩在线观看| 亚洲自拍偷在线| 日本 av在线| 88av欧美| 亚洲人成77777在线视频| 国产精品av久久久久免费| 日韩欧美免费精品| 久久久久久久精品吃奶| 1024视频免费在线观看| a在线观看视频网站| 性少妇av在线| 亚洲欧美精品综合久久99| 国产精品一区二区在线不卡| 真人做人爱边吃奶动态| 一本综合久久免费| 在线观看免费视频日本深夜| 国产精品 欧美亚洲| 禁无遮挡网站| 国产成人一区二区三区免费视频网站| 91麻豆av在线| 黄片播放在线免费| 中文字幕久久专区| 女人爽到高潮嗷嗷叫在线视频| 免费无遮挡裸体视频| 国产99白浆流出| 嫩草影视91久久| 精品久久久久久,| 国产免费av片在线观看野外av| 国产欧美日韩精品亚洲av| 久久久国产精品麻豆| 人成视频在线观看免费观看| 国产欧美日韩一区二区精品| 欧美丝袜亚洲另类 | 母亲3免费完整高清在线观看| 久久性视频一级片| 欧美日韩一级在线毛片| 亚洲成国产人片在线观看| 国产片内射在线| 免费一级毛片在线播放高清视频 | 色综合站精品国产| 国产成人精品无人区| 啦啦啦免费观看视频1| 咕卡用的链子| 免费久久久久久久精品成人欧美视频| 欧美成人一区二区免费高清观看 | 国产亚洲精品一区二区www| 真人一进一出gif抽搐免费| 亚洲av熟女| 亚洲七黄色美女视频| a在线观看视频网站| 日韩一卡2卡3卡4卡2021年| 亚洲精品国产色婷婷电影| 精品高清国产在线一区| 国产亚洲精品av在线| 亚洲 欧美 日韩 在线 免费| 韩国精品一区二区三区| av欧美777| 99在线人妻在线中文字幕| 欧美激情高清一区二区三区| 亚洲国产欧美一区二区综合| 一级作爱视频免费观看| 成年女人毛片免费观看观看9| 国产片内射在线| 中文字幕久久专区| 久久影院123| 日日爽夜夜爽网站| 亚洲午夜精品一区,二区,三区| 亚洲中文字幕一区二区三区有码在线看 | 久久久精品欧美日韩精品| 日日干狠狠操夜夜爽| 国产色视频综合| 很黄的视频免费| 国内久久婷婷六月综合欲色啪| 国产精品电影一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 国产高清激情床上av| 国产黄a三级三级三级人| 欧美在线一区亚洲| 久9热在线精品视频| 老司机午夜福利在线观看视频| 级片在线观看| 天堂影院成人在线观看| 最新美女视频免费是黄的| 女人被狂操c到高潮| 校园春色视频在线观看| 国产激情欧美一区二区| 免费无遮挡裸体视频| 日韩精品青青久久久久久| 欧美黑人精品巨大| 欧美日韩中文字幕国产精品一区二区三区 | 变态另类丝袜制服| а√天堂www在线а√下载| 亚洲人成电影观看| 免费在线观看亚洲国产| 亚洲第一电影网av| 波多野结衣巨乳人妻| 一级a爱视频在线免费观看| 无遮挡黄片免费观看| 亚洲熟妇中文字幕五十中出| 大型av网站在线播放| 国产精品美女特级片免费视频播放器 | 精品久久久久久久毛片微露脸| 中国美女看黄片| 亚洲欧美激情综合另类| 久久影院123| 日韩国内少妇激情av| 午夜精品国产一区二区电影| 国产高清激情床上av| av在线天堂中文字幕| 欧美乱妇无乱码| 亚洲国产精品久久男人天堂| 在线观看舔阴道视频| 国产精品1区2区在线观看.| 国产精品影院久久| 在线观看一区二区三区| 久久久水蜜桃国产精品网| 国产成人精品久久二区二区91| 一区二区三区精品91| 很黄的视频免费| 国产日韩一区二区三区精品不卡| 色婷婷久久久亚洲欧美| 大型黄色视频在线免费观看| 男女做爰动态图高潮gif福利片 | 好男人在线观看高清免费视频 | 成人亚洲精品av一区二区| 午夜免费鲁丝| 精品国产超薄肉色丝袜足j| 99久久久亚洲精品蜜臀av| 波多野结衣av一区二区av| 亚洲精品粉嫩美女一区| 99国产精品一区二区蜜桃av| 亚洲人成伊人成综合网2020| 国产aⅴ精品一区二区三区波| 亚洲国产欧美网| 亚洲精品国产区一区二| 日本欧美视频一区| 亚洲精品国产精品久久久不卡| 人成视频在线观看免费观看| 黄色片一级片一级黄色片| 亚洲欧洲精品一区二区精品久久久| 老司机深夜福利视频在线观看| 午夜福利高清视频| 国产免费av片在线观看野外av| 日韩欧美免费精品| 国产麻豆成人av免费视频| 91大片在线观看| 超碰成人久久| 日韩免费av在线播放| 国产精品免费一区二区三区在线| 高清黄色对白视频在线免费看| 夜夜夜夜夜久久久久| 精品国产乱子伦一区二区三区| 欧美日韩精品网址| 久久久国产欧美日韩av| 亚洲成av片中文字幕在线观看| 国产成人精品无人区| 欧美日韩福利视频一区二区| av天堂久久9| 亚洲一区中文字幕在线| 免费少妇av软件| 亚洲熟妇中文字幕五十中出| 淫秽高清视频在线观看| 国产成人影院久久av| 国产单亲对白刺激| 免费一级毛片在线播放高清视频 | 国产精品一区二区三区四区久久 | 99精品欧美一区二区三区四区| 日韩av在线大香蕉| 99久久精品国产亚洲精品| 精品国产乱子伦一区二区三区| 又黄又爽又免费观看的视频| 国产野战对白在线观看| 啦啦啦 在线观看视频| 黄片大片在线免费观看| 天堂影院成人在线观看| 欧美精品啪啪一区二区三区| 国产野战对白在线观看| 欧美日本中文国产一区发布| 最近最新免费中文字幕在线| 精品午夜福利视频在线观看一区| 国产免费男女视频| 99国产极品粉嫩在线观看| 久久久水蜜桃国产精品网| 在线观看免费日韩欧美大片| 看黄色毛片网站| 一区二区三区激情视频| 在线观看免费视频网站a站| а√天堂www在线а√下载| 精品一区二区三区四区五区乱码| 啦啦啦观看免费观看视频高清 | cao死你这个sao货| 日韩一卡2卡3卡4卡2021年| 成人精品一区二区免费| ponron亚洲| 国产三级黄色录像| 午夜精品在线福利| 我的亚洲天堂| 国产成人欧美| 99精品欧美一区二区三区四区| 18禁裸乳无遮挡免费网站照片 | 又大又爽又粗| 欧洲精品卡2卡3卡4卡5卡区| 亚洲最大成人中文| 国产xxxxx性猛交| 久久人妻福利社区极品人妻图片| 99久久国产精品久久久| 亚洲欧美日韩高清在线视频| 女同久久另类99精品国产91| 日韩欧美一区二区三区在线观看| 日本a在线网址| 成年女人毛片免费观看观看9| 国产伦一二天堂av在线观看| 97碰自拍视频| 日韩三级视频一区二区三区| 大香蕉久久成人网| 国产亚洲欧美精品永久| 亚洲av成人一区二区三| 12—13女人毛片做爰片一| 日本精品一区二区三区蜜桃| 18禁裸乳无遮挡免费网站照片 | 人成视频在线观看免费观看| 免费不卡黄色视频| 美国免费a级毛片| 成人三级黄色视频| 搡老妇女老女人老熟妇| 久久久久国产一级毛片高清牌| 日韩精品中文字幕看吧| av在线播放免费不卡| av超薄肉色丝袜交足视频| 校园春色视频在线观看| 日韩国内少妇激情av| 亚洲国产看品久久| а√天堂www在线а√下载| 女人被狂操c到高潮| 亚洲国产日韩欧美精品在线观看 | 精品久久久久久成人av| 搞女人的毛片| 亚洲国产看品久久| 欧美成人性av电影在线观看| 性少妇av在线| 搡老熟女国产l中国老女人| 每晚都被弄得嗷嗷叫到高潮| 久久精品国产亚洲av香蕉五月| 一级片免费观看大全| 日本在线视频免费播放| 一边摸一边抽搐一进一出视频| 在线观看一区二区三区| 电影成人av| 天堂影院成人在线观看| 日本黄色视频三级网站网址| 十八禁网站免费在线| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲情色 制服丝袜| 老司机靠b影院| 成人亚洲精品一区在线观看| 黄色片一级片一级黄色片| 欧美乱码精品一区二区三区| 亚洲国产精品999在线| 国产精品自产拍在线观看55亚洲| 亚洲成a人片在线一区二区| 女人爽到高潮嗷嗷叫在线视频| 国产三级在线视频| av免费在线观看网站| 亚洲专区字幕在线| 欧美日韩黄片免| 黄色视频不卡| www国产在线视频色| 好看av亚洲va欧美ⅴa在| 女人高潮潮喷娇喘18禁视频| 国产精品电影一区二区三区| 日韩精品中文字幕看吧| 国产又色又爽无遮挡免费看| 夜夜爽天天搞| 久久国产精品人妻蜜桃| 亚洲国产欧美日韩在线播放| 日韩三级视频一区二区三区| 日本撒尿小便嘘嘘汇集6| av电影中文网址| 精品国产一区二区三区四区第35| 母亲3免费完整高清在线观看| 高清在线国产一区| 成人特级黄色片久久久久久久| 久久中文看片网| 69精品国产乱码久久久| 久久国产亚洲av麻豆专区| 久久久久亚洲av毛片大全| 90打野战视频偷拍视频| e午夜精品久久久久久久| 一本大道久久a久久精品| 国产在线观看jvid| 亚洲av熟女| 欧美中文综合在线视频| 亚洲午夜精品一区,二区,三区| 国产91精品成人一区二区三区| 给我免费播放毛片高清在线观看| 一级毛片女人18水好多| 午夜精品国产一区二区电影| 免费看美女性在线毛片视频| 成人三级做爰电影| 老司机在亚洲福利影院| 亚洲成人国产一区在线观看| 好男人在线观看高清免费视频 | 日韩欧美在线二视频| 亚洲熟女毛片儿| 长腿黑丝高跟| 美女 人体艺术 gogo| 91九色精品人成在线观看| 久久国产精品男人的天堂亚洲| 一级片免费观看大全| 最近最新中文字幕大全电影3 | 午夜久久久久精精品| 亚洲电影在线观看av| 久久精品影院6| 亚洲人成伊人成综合网2020| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美日本中文国产一区发布| 99精品久久久久人妻精品| 国产1区2区3区精品| 欧美av亚洲av综合av国产av| 最近最新中文字幕大全免费视频| 成人18禁高潮啪啪吃奶动态图| 精品电影一区二区在线| 亚洲精品粉嫩美女一区| 很黄的视频免费| 777久久人妻少妇嫩草av网站| 日本精品一区二区三区蜜桃| x7x7x7水蜜桃| 曰老女人黄片| 亚洲精品粉嫩美女一区| 老司机午夜福利在线观看视频| 国产精品一区二区精品视频观看| 亚洲成人国产一区在线观看| bbb黄色大片| 亚洲午夜精品一区,二区,三区| 亚洲av成人av| av片东京热男人的天堂| 黄色成人免费大全| 久久香蕉国产精品| 夜夜夜夜夜久久久久| 一边摸一边做爽爽视频免费| 99精品欧美一区二区三区四区| 极品人妻少妇av视频| 在线观看66精品国产| 亚洲精品美女久久av网站| av片东京热男人的天堂| 亚洲欧美日韩另类电影网站| 免费高清在线观看日韩| 欧洲精品卡2卡3卡4卡5卡区| 黄色 视频免费看| 欧美成人性av电影在线观看| 欧美久久黑人一区二区| 欧美+亚洲+日韩+国产| 免费观看人在逋| 两个人看的免费小视频| 女人高潮潮喷娇喘18禁视频| 午夜福利在线观看吧| 久久 成人 亚洲| 精品一品国产午夜福利视频| 别揉我奶头~嗯~啊~动态视频| 老司机深夜福利视频在线观看| 国产精品免费一区二区三区在线| 性少妇av在线| 成人特级黄色片久久久久久久| 免费高清在线观看日韩| svipshipincom国产片| 大型黄色视频在线免费观看| 狠狠狠狠99中文字幕| 国产精品电影一区二区三区| 精品国产一区二区三区四区第35| tocl精华| 午夜福利影视在线免费观看| 美女高潮喷水抽搐中文字幕| x7x7x7水蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 女生性感内裤真人,穿戴方法视频| 狂野欧美激情性xxxx| 亚洲中文字幕一区二区三区有码在线看 | 国产精品久久久久久亚洲av鲁大| 老司机午夜十八禁免费视频| 亚洲国产精品久久男人天堂| 12—13女人毛片做爰片一| 制服丝袜大香蕉在线| 人人澡人人妻人| 中文字幕人成人乱码亚洲影| 一区在线观看完整版| 757午夜福利合集在线观看| 19禁男女啪啪无遮挡网站| 免费久久久久久久精品成人欧美视频| 老熟妇仑乱视频hdxx| 日本 av在线| 国产亚洲av嫩草精品影院| 国产精品亚洲美女久久久| 88av欧美| 免费久久久久久久精品成人欧美视频| 身体一侧抽搐| 免费搜索国产男女视频| 窝窝影院91人妻| 欧美黑人精品巨大| 久久久精品国产亚洲av高清涩受| 变态另类成人亚洲欧美熟女 | 丁香六月欧美| 国产高清激情床上av| 非洲黑人性xxxx精品又粗又长| 在线观看66精品国产| 成在线人永久免费视频| 黄片播放在线免费| 在线观看一区二区三区| 久久这里只有精品19| 在线观看www视频免费| 51午夜福利影视在线观看| 可以在线观看毛片的网站| 在线观看免费视频网站a站| 精品无人区乱码1区二区| 日韩中文字幕欧美一区二区| 精品免费久久久久久久清纯| 精品国内亚洲2022精品成人| 男女做爰动态图高潮gif福利片 | 色综合婷婷激情| 亚洲成av片中文字幕在线观看| 可以免费在线观看a视频的电影网站| 久久久久久人人人人人| 又大又爽又粗| 老司机深夜福利视频在线观看| 久久久水蜜桃国产精品网| 亚洲天堂国产精品一区在线| 精品国产一区二区三区四区第35| 亚洲人成伊人成综合网2020| 色哟哟哟哟哟哟| 午夜福利在线观看吧| 免费在线观看视频国产中文字幕亚洲| 亚洲精品美女久久久久99蜜臀| 亚洲精品国产一区二区精华液| 日本五十路高清| 欧美人与性动交α欧美精品济南到| 麻豆成人av在线观看| 亚洲七黄色美女视频| 欧美乱码精品一区二区三区| 国产一区二区在线av高清观看| 久久国产精品人妻蜜桃| 在线观看免费视频网站a站| 国产亚洲欧美98| 日韩欧美国产一区二区入口| 99在线视频只有这里精品首页| 国产精品免费视频内射| 黄频高清免费视频| 91精品国产国语对白视频| 在线免费观看的www视频| 免费久久久久久久精品成人欧美视频| av中文乱码字幕在线| 久9热在线精品视频| 亚洲欧美激情在线| 国产成人精品久久二区二区91| 9191精品国产免费久久| 在线天堂中文资源库| 亚洲精品国产色婷婷电影| 两个人免费观看高清视频| 91麻豆av在线| 久久久久亚洲av毛片大全| 国产精品美女特级片免费视频播放器 | 最新美女视频免费是黄的| 欧美在线黄色| 操美女的视频在线观看| 亚洲全国av大片| 国产熟女午夜一区二区三区| 欧美成人午夜精品| 日韩av在线大香蕉| 一本综合久久免费| 国产一区二区三区在线臀色熟女| 亚洲欧美日韩高清在线视频| 欧美午夜高清在线| 欧美另类亚洲清纯唯美| 国产又爽黄色视频| 狂野欧美激情性xxxx| 最好的美女福利视频网| 国产熟女午夜一区二区三区| 亚洲电影在线观看av| 国产精品久久视频播放| 女性被躁到高潮视频| 搡老妇女老女人老熟妇| 女人爽到高潮嗷嗷叫在线视频| 女性生殖器流出的白浆| 久久青草综合色| 亚洲av美国av| 最近最新中文字幕大全电影3 | 一区二区三区激情视频| 亚洲精品久久国产高清桃花| 777久久人妻少妇嫩草av网站| 啪啪无遮挡十八禁网站| 女性被躁到高潮视频| 两个人视频免费观看高清| 欧美 亚洲 国产 日韩一| 亚洲性夜色夜夜综合| 在线国产一区二区在线| 女性生殖器流出的白浆| 国产精品野战在线观看| 99在线人妻在线中文字幕| 欧美最黄视频在线播放免费| 亚洲av成人av| 啦啦啦韩国在线观看视频| 岛国视频午夜一区免费看| 国产日韩一区二区三区精品不卡| 亚洲,欧美精品.| 成人18禁高潮啪啪吃奶动态图| 91麻豆精品激情在线观看国产| 在线播放国产精品三级| 色综合亚洲欧美另类图片| 国产精品 欧美亚洲| 一a级毛片在线观看| 亚洲一区高清亚洲精品| 国产在线观看jvid| 精品日产1卡2卡| 日韩欧美一区视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲视频免费观看视频| www.熟女人妻精品国产| 亚洲第一青青草原| 黄色a级毛片大全视频| 91精品三级在线观看| 欧美最黄视频在线播放免费| 欧美激情极品国产一区二区三区| 91精品国产国语对白视频| 久久久久久久午夜电影| 色播在线永久视频| 国产精品av久久久久免费| 国产成+人综合+亚洲专区| 中文字幕另类日韩欧美亚洲嫩草| 一级毛片高清免费大全| 国产欧美日韩精品亚洲av| 国产精品 欧美亚洲| 搡老岳熟女国产| 日韩有码中文字幕| 国产视频一区二区在线看| 免费女性裸体啪啪无遮挡网站| 亚洲色图综合在线观看| 午夜福利一区二区在线看| 热99re8久久精品国产| 麻豆成人av在线观看| 亚洲av电影在线进入| 免费观看精品视频网站| 国产精品一区二区三区四区久久 | 午夜精品在线福利| 色哟哟哟哟哟哟| 国产欧美日韩综合在线一区二区| 亚洲伊人色综图| www国产在线视频色| 久久精品aⅴ一区二区三区四区| 成人三级做爰电影| 久久精品国产亚洲av香蕉五月| 9色porny在线观看| 男女床上黄色一级片免费看| 久久青草综合色| 成人亚洲精品一区在线观看| 亚洲av成人一区二区三| 天天添夜夜摸| 又黄又粗又硬又大视频| 高清黄色对白视频在线免费看| av免费在线观看网站| 国产色视频综合| 国产高清视频在线播放一区| 国内精品久久久久久久电影| 日韩成人在线观看一区二区三区| 欧美另类亚洲清纯唯美| av欧美777| 成人18禁高潮啪啪吃奶动态图| 90打野战视频偷拍视频|