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

    甘菊ClHSP70與ClHSP90基因的克隆及表達(dá)分析

    2016-09-02 02:49:14陳東亮黃叢林
    西北植物學(xué)報(bào) 2016年7期
    關(guān)鍵詞:甘菊擬南芥結(jié)構(gòu)域

    羅 昌,陳東亮,程 曦,黃叢林

    (北京市農(nóng)林科學(xué)院 北京農(nóng)業(yè)生物技術(shù)研究中心,北京100097)

    ?

    甘菊ClHSP70與ClHSP90基因的克隆及表達(dá)分析

    羅昌,陳東亮,程曦,黃叢林*

    (北京市農(nóng)林科學(xué)院 北京農(nóng)業(yè)生物技術(shù)研究中心,北京100097)

    該研究以甘菊(Chrysanthemumlavandulifolium)為實(shí)驗(yàn)材料,通過RT-PCR方法從甘菊轉(zhuǎn)錄組數(shù)據(jù)中分離出熱激蛋白合成相關(guān)基因,命名為ClHSP70和ClHSP90。序列分析表明,ClHSP70基因ORF全長(zhǎng)為2 559 bp,編碼852個(gè)氨基酸,蛋白功能區(qū)預(yù)測(cè)表明含有典型的HSP70蛋白NBD和SBD保守結(jié)構(gòu)域;ClHSP90基因ORF全長(zhǎng)為2 094 bp,編碼697個(gè)氨基酸,含有HATPase結(jié)構(gòu)域和HSP90保守結(jié)構(gòu)域。生物信息學(xué)分析表明,甘菊ClHSP70與大豆(Glycinemax)和煙草(Nicotianatomentosiformis)HSP70蛋白有較高的一致性,ClHSP90基因編碼的氨基酸序列與紫莖澤蘭(Ageratinaadenophora)HSP90高度相似;實(shí)時(shí)熒光定量表達(dá)分析表明,在42 ℃處理不同時(shí)間,甘菊葉片中ClHSP70和ClHSP90基因表達(dá)均在0.5 h時(shí)顯著增加,1 h達(dá)到最大值,2 h后緩慢下降;不同組織表達(dá)分析表明,甘菊在42 ℃處理1 h后,ClHSP70在成熟葉中的表達(dá)量顯著高于嫩葉和根等其他組織;ClHSP90在成熟莖中的表達(dá)量最高。研究說明,ClHSP70和ClHSP90基因具有熱激蛋白特征,參與了甘菊熱脅迫應(yīng)答過程,該研究結(jié)果為以后深入研究其基因功能奠定了基礎(chǔ)。

    甘菊;ClHSP70;ClHSP90; 熱脅迫

    熱激蛋白(heat shock protein, HSPs)是在熱激條件下被誘導(dǎo)且表達(dá)增加的一部分特異蛋白質(zhì),在植物應(yīng)對(duì)外界環(huán)境壓力,如干旱、鹽、低溫及高溫等時(shí)起著重要的作用[1]。HSPs的積累與生物耐熱性呈正相關(guān),也就是說,在熱脅迫下HSPs的積累可以提高植株的耐熱性。

    HSPs種類存在范圍很廣,根據(jù)其分子量大小,可分為5類:HSP110、HSP90、HSP70、HSP60和‘小HSPs’[2]。在植物熱激蛋白家族中,HSP70最為保守和普遍,HSP70在氨基酸序列上高度相似,含有核酸結(jié)合區(qū)和底物結(jié)合區(qū)2個(gè)功能區(qū),N-末端具有典型的約45 kD的ATPase區(qū)域。研究表明,HSP70在熱脅迫條件下能夠被顯著誘導(dǎo),小球藻HSP70B含量在熱脅迫2 h后增加最明顯[3],花生(ArachishypogaeaLinn.)AhHSP70在高溫脅迫3 h后表達(dá)明顯升高[4]。近年來在擬南芥[5]、煙草[1]、玉米[6]等研究中也證明了HSP70與植物耐熱性相關(guān)。HSP90是一個(gè)受ATP調(diào)節(jié)的二聚體分子伴侶,在分子進(jìn)化上高度保守,每個(gè)單體由ATPase結(jié)構(gòu)域、中間結(jié)構(gòu)域和介導(dǎo)HSP90二聚化的C端結(jié)構(gòu)域組成[7-8]。目前對(duì)于植物中HSP90蛋白的研究相對(duì)較少,擬南芥[9]和水稻[10]中的研究表明受到熱脅迫時(shí)HSP90表達(dá)量顯著增加,HSP90蛋白的積累量明顯提高,但其結(jié)構(gòu)和功能仍待進(jìn)一步探究。另外,有研究證實(shí)HSP90在植物生長(zhǎng)發(fā)育進(jìn)程中起著至關(guān)重要的作用,其功能的缺失會(huì)引起植物形態(tài)和生理特性發(fā)生改變[11-12]。此外,由于熱誘導(dǎo)和組成型表達(dá)的熱激蛋白可定位于細(xì)胞的細(xì)胞質(zhì)和多種細(xì)胞器中,包括葉綠體、內(nèi)質(zhì)網(wǎng)[13-14]等,且執(zhí)行著不同的功能,也可在不同的器官中被誘導(dǎo),可在根、莖、葉、種子以及幼苗中產(chǎn)生,因此高溫脅迫下植物HSPs基因的表達(dá)也呈現(xiàn)出明顯的組織特異性。

    菊花(Chrysanthemum×morifoliumRamat.)是多年生草本植物,起源于中國(guó),世界上均有廣泛栽培,在中國(guó)花卉生產(chǎn)中占有重要地位,周年生產(chǎn)鮮花是菊花生產(chǎn)的迫切需求。但是,周年生產(chǎn)過程中菊花在現(xiàn)蕾階段容易受到夏季高溫脅迫,導(dǎo)致花芽早熟或夭折。因此,了解菊花在受到熱脅迫時(shí)的耐熱機(jī)制對(duì)于培育適合周年生產(chǎn)的菊花品種很有必要。甘菊(Chrysanthemumlavandulifolium) 為菊科菊屬二倍體野生種(2n=2x=18),是菊花原始種之一,生長(zhǎng)于海拔200~2 800 m平原及山坡,有較強(qiáng)抗逆性,其生長(zhǎng)季節(jié)和開花習(xí)性與栽培菊花相似[15-16]。本研究參照甘菊均一化轉(zhuǎn)錄組數(shù)據(jù)庫(kù),利用RT-PCR技術(shù),分離得到ClHSP70和ClHSP90的cDNA全長(zhǎng)序列,利用實(shí)時(shí)熒光定量PCR技術(shù)分析了2個(gè)基因在高溫下不同處理時(shí)間和不同組織部位中的表達(dá)特性,為ClHSP70和ClHSP90基因的功能鑒定以及菊花耐熱育種奠定基礎(chǔ)。

    1 材料與方法

    1.1植物材料及處理

    從北京農(nóng)業(yè)生物技術(shù)研究中心組培室中選取15株生長(zhǎng)水平一致的甘菊組培苗為實(shí)驗(yàn)材料。將光照培養(yǎng)箱溫度設(shè)置為42 ℃,每3株作為一個(gè)重復(fù),分別處理0、0.5、1、2和4 h時(shí)選取上部新鮮葉片置于液氮中冷凍,并放入-80 ℃冰箱中備用;42 ℃處理1 h, 選取植株嫩葉、成熟葉、嫩莖、成熟莖、根置于液氮中冷凍,并放入-80 ℃冰箱中備用。

    1.2方法

    1.2.1總RNA提取及甘菊ClHSP70和ClHSP90基因克隆用植物RNA小量提取試劑盒(MiniBEST Plant RNA Extraction Kit,TaKaRa公司)提取甘菊葉片總RNA,再用PrimeScriptTMⅡ 1st Strand cDNA Synthesis Kit(TaKaRa公司)將所提取RNA反轉(zhuǎn)錄為cDNA并檢測(cè)其濃度。

    從甘菊均一化轉(zhuǎn)錄組數(shù)據(jù)中篩選ClHSP70和ClHSP90基因序列信息,根據(jù)在NCBI網(wǎng)站上預(yù)測(cè)2個(gè)基因的ORF設(shè)計(jì)引物(表1)。PCR反應(yīng)體系(50 μL):2 μL cDNA模板(100 mmol/L),0.5 μL PrimeSTAR HS DNA Polymerase (2.5 U/μL),10 μL 5×PrimeSTAR Buffer (Mg2+plus),4 μL dNTP Mixture (各2.5 mmol/L),正反向引物各2 μL 以及29.5 μL ddH2O。PCR擴(kuò)增條件為98 ℃預(yù)變性5 min;98 ℃變性10 s,58 ℃復(fù)性5 s,72 ℃延伸15 s,30個(gè)循環(huán);72 ℃延伸6 min。PCR產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳,回收純化目的片段后連接到pGEM-T Easy載體上,42 ℃熱激轉(zhuǎn)化大腸桿菌DH5α感受態(tài)細(xì)胞,進(jìn)行Amp抗性篩選和X-gal/IPTG藍(lán)白斑篩選。挑取白色菌落進(jìn)行PCR檢測(cè),陽(yáng)性克隆送往上海生物工程有限公司進(jìn)行測(cè)序。用Segman序列拼接軟件對(duì)測(cè)序結(jié)果進(jìn)行拼接,分別獲得ClHSP70和ClHSP90的ORF全長(zhǎng)序列。

    表1 本實(shí)驗(yàn)相關(guān)引物

    1.2.2ClHSP70和ClHSP90基因的生物信息學(xué)分析利用NCBI數(shù)據(jù)庫(kù)中Blast工具對(duì)ClHSP70和ClHSP90基因序列進(jìn)行比對(duì)和同源性分析(http://blast.ncbi.nlm.nih.gov/Blast.cgi);運(yùn)用Finder軟件(http://www.ncbi.nlm.nih.gov/gorf/g orf.html)進(jìn)行開放閱讀框分析;在NCBI Conserved Domain Database (http://www.ncbi.nlm. nih.gov/Structure/cdd/wrpsb.cgi)中進(jìn)行蛋白質(zhì)保守結(jié)構(gòu)域預(yù)測(cè);通過ExPASy-ProtParam軟件(http://web.expasy.org/protparam)預(yù)測(cè)蛋白的分子量及等電點(diǎn);SMART在線軟件對(duì)蛋白質(zhì)結(jié)構(gòu)域進(jìn)行分析預(yù)測(cè)。利用SignalP (http://www.cbs.dtu.dk/services/SignalP/)在線分析蛋白質(zhì)的信號(hào)肽。使用Clustal X軟件進(jìn)行氨基酸序列多重比對(duì)[17],利用MEGA5.0軟件構(gòu)建進(jìn)化樹(Neighbor- Joining)[18],并進(jìn)行Boot-strap檢測(cè)。

    1.2.3ClHSP70和ClHSP90基因表達(dá)分析用植物RNA小量提取試劑盒(MiniBEST Plant RNA Extraction Kit,TaKaRa公司)提取42 ℃不同時(shí)間(0、0.5、1、2、4 h)處理后的甘菊葉片RNA,已經(jīng)42 ℃處理1 h后不同組織部位(嫩葉、成熟葉、嫩莖、成熟莖、根)的RNA。用PrimeScriptTMⅡ 1st Strand cDNA Synthesis Kit(TaKaRa公司)將所提取RNA反轉(zhuǎn)錄為cDNA。根據(jù)ClHSP70和ClHSP90基因的全長(zhǎng)序列,用Premier 5.0設(shè)計(jì)Real-time PCR引物(表1),采用Actin基因作為內(nèi)標(biāo)。參照熒光定量試劑盒SYBR?Premix EX TaqTMⅡ(TilRNaseH Plus,TaKaRa)說明書進(jìn)行操作:SYBR Premix EX Taq (2×) 10 μL,上、下游引物(10 μmol/L)各0.4 μL,ROX Reference Dye (50×) 0.4 μL, cDNA模板 (100 mmol/L) 1 μL, ddH2O 7.8 μL,總體積20 μL。反應(yīng)程序?yàn)?5 ℃預(yù)變性10 min,95 ℃變性15 s,60 ℃退火1 min,40個(gè)循環(huán)后作熔解曲線。每個(gè)樣品均重復(fù)3次,在ABIStep-OnePlusTMReal-time PCR儀(美國(guó)ABI公司)上進(jìn)行,采用2-ΔΔCt法計(jì)算相對(duì)表達(dá)量。

    2 結(jié)果與分析

    2.1ClHSP70和ClHSP90基因克隆

    根據(jù)甘菊轉(zhuǎn)錄組數(shù)據(jù)Unigene blastx注釋結(jié)果篩選HSP70和HSP90序列,在NCBI中進(jìn)行開放閱讀框分析并設(shè)計(jì)特異性引物,以甘菊葉片cDNA為模板擴(kuò)增到ClHSP70和ClHSP90編碼框序列,擴(kuò)增產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳分析,2個(gè)基因擴(kuò)增條帶長(zhǎng)度分別為2 500 bp和2000 bp左右(圖 1),片段經(jīng)測(cè)序后正反向拼接得到基因的全長(zhǎng)。分析結(jié)果表明,ClHSP70含有1個(gè)2 559 bp開放閱讀框,編碼852個(gè)氨基酸(圖 2);ClHSP90含有一個(gè)2 094 bp開放閱讀框,編碼697個(gè)氨基酸(圖 3)。

    2.2ClHSP70和ClHSP90基因生物信息學(xué)分析

    經(jīng)Protparatam在線工具(http://web.expasy.org/protparam/)分析,結(jié)果顯示(表 2)ClHSP70蛋白的總平均疏水性為-0.444, ClHSP90蛋白的總平均疏水性為-0.604,均為親水性蛋白;總的帶負(fù)電荷的氨基酸比重均略高于帶正電荷的氨基酸,表現(xiàn)為酸性蛋白;從穩(wěn)定系數(shù)推測(cè)ClHSP70為不穩(wěn)定蛋白,ClHSP90為穩(wěn)定蛋白。通過二級(jí)結(jié)構(gòu)預(yù)測(cè)發(fā)現(xiàn)ClHSP70含有40.73%的α-螺旋,5.87%的β-折疊和36.85%的隨機(jī)卷曲;ClHSP90則含有51.51%的α-螺旋、6.31%的β-折疊和24.82%的隨機(jī)卷曲。SignalP在線分析ClHSP70和ClHSP90均不含信號(hào)肽序列,也無跨膜結(jié)構(gòu)。

    M. DL5000; 1. ClHSP70;2. ClHSP90圖1 甘菊ClHSP70和ClHSP90基因擴(kuò)增結(jié)果Fig. 1 PCR amplification results of ClHSP70 and ClHSP90 gene from C. lavandulifolium

    下劃線部分為N端NBD結(jié)構(gòu)域(A)和N端保守的ATPase結(jié)構(gòu)域(B);灰色陰影部分為高度保守的ClHSP90氨基酸序列(B);黑色陰影部分為C-末端胞質(zhì)亞型特殊信號(hào)的基序(B);*.表示終止密碼子圖2 甘菊ClHSP70 (A)和ClHSP90 (B)基因的核苷酸序列及其編碼的氨基酸序列Underline represents NBD domain of N-terminus domain(A)and the ATP binding domain of N-terminus domain (B); shaded gray represents highly conserved amino acid ClHSP90 signature sequence motifs characteristic (B); shaded black represents the C-terminus cytosolic isoform specific signature motif (MEEVD) (B); *. Stop codonFig. 2 The nucleotide and deduced amino acid sequences of ClHSP70 (A) and ClHSP90 (B) in C. lavandulifolium

    NtHSP70-15.煙草 (XP009624163.1);GmHSP70-14.大豆(XP 003546366.1);AtHSP70. 擬南芥(NP 850984.4);BoHSP70-14.甘藍(lán)(XP 013591718.1)圖3 甘菊ClHSP70與其他物種HSPs多序列比對(duì)NtHSP70-15.Nicotiana tomentosiformis (XP 009624163.1);GmHSP70-14. Glycine max(XP 003546366.1);AtHSP70. Arabidopsis thaliana (NP 850984.4);BoHSP70-14. Brassica oleracea var.oleracea (XP 01591718.1).Fig. 3 Multiple alignment of ClHSP70 with HSPs from other plant species

    蛋白性質(zhì)ProteincharacterClHSP70ClHSP90總氨基酸數(shù)目Numberofaminoacids/aa852697蛋白分子質(zhì)量Molecularweight/kD94.179.9理論等電點(diǎn)pI5.104.98蛋白分子式FormulaC4133H6584N1136O1302S34C3544H5649N917O1124S26總帶負(fù)電荷的殘基數(shù)Totalnumberofnegativelychargedresidues(Asp+Glu)135138總帶正電荷的殘基數(shù)Totalnumberofpositivelychargedresidues(Arg+Lys)103106脂肪族指數(shù)Aliphaticindex78.9182.50不穩(wěn)定系數(shù)Instabilityindex43.08,不穩(wěn)定Unstable37.55,穩(wěn)定Stable總平均疏水性Grandaverageofhydropathicity-0.444-0.604

    氨基酸序列比對(duì)結(jié)果顯示,甘菊ClHSP70與大豆(Glycinemax)和煙草(Nicotianatomentosiformis) HSP70蛋白有較高的一致性,分別為78%和76%, 其次是與擬南芥(Arabidopsisthaliana)、甘藍(lán)(Brassicaoleraceavar.oleracea)同為74%(圖 3);甘菊ClHSP90與紫莖澤蘭(Ageratinaadenophora)HSP90有顯著的一致性,高達(dá)93%,其次是與煙草為92%,與番茄(Solanumlycopersicum)、擬南芥一致性分別為91%和90%(圖 4)。說明ClHSP70和ClHSP90具有熱激蛋白家族高度保守的特征。保守域預(yù)測(cè)結(jié)果(圖 5)表明ClHSP70蛋白C端3~691aa為HSP70蛋白家族保守區(qū)域,782~818 aa為一段低復(fù)雜性區(qū)段(low-complexity region,LCR)。N端9~381aa是核酸結(jié)合區(qū)(nucleotide binding domain, NBD)(圖 2,A),該結(jié)合區(qū)氨基酸序列高度保守,具有結(jié)合并水解ATP的活性中心,它與ATP的結(jié)合-水解過程調(diào)控著它對(duì)多肽的親和力,34~219aa為底物結(jié)合區(qū)(substrate binding domain, SBD)。

    ClHSP90蛋白質(zhì)結(jié)構(gòu)中含有N端高度保守的ATPase結(jié)構(gòu)域29~178aa(圖 5),具有ATPase活性,ClHSP90蛋白質(zhì)184~697aa。另外,ClHSP90氨

    NtHSP80-like.煙草 (XP 009594054.1);AaHSP90.紫莖澤蘭 (ABX76302.1); SlHSP90-1.番茄 (NP 001234436.1);AtHSP81-2.擬南芥 (NP 200414.1)圖4 甘菊ClHSP90與其他物種HSPs多序列比對(duì)NtHSP80-like. Nicotiana tomentosiformis (XP 009594054.1);AaHSP90. Ageratina adenophora (ABX76302.1); SlHSP90-1. Solanum lycopersicum (NP 001234436.1);AtHSP81-2. Arabidopsis thaliana (NP 200414.1);GmHSP90-2. Glycine max (NP 001236599.1)Fig. 4 Multiple alignment of ClHSP90 with HSPs from other plant species

    基酸序列中還包括HSP90家族5個(gè)特征標(biāo)記序列:NKEIFL RELISNSSDALDKIR、LGTIARSGT、MIGQFGVGF YSAYLVA、IKLYVRRVFI、GIVDSEDLPLNIS RE,C端末尾基序?yàn)镸EEVD,推測(cè)該蛋白為細(xì)胞質(zhì)熱激蛋白(圖 2,B)。

    系統(tǒng)進(jìn)化樹分析表明,ClHSP70基因與擬南芥AtHSP70 (NP-850984.4)、甘藍(lán)BoHSP70-14 (XP-013591718.1)的關(guān)系較近,處于同一分支(圖6,A)。ClHSP90與同為菊科的紫莖澤蘭的AaHSP90(ABX76302.1)關(guān)系最近,紫莖澤蘭HSP90已驗(yàn)證參與脅迫應(yīng)答,其次是煙草、番茄,與小麥、擬南芥、甘藍(lán)等的HSP90類基因關(guān)系較遠(yuǎn) (圖6,B)。

    2.3ClHSP70和ClHSP90基因表達(dá)分析

    實(shí)時(shí)熒光定量PCR分析結(jié)果表明,甘菊在42 ℃不同時(shí)間處理?xiàng)l件下,葉片中ClHSP70和ClHSP90基因表達(dá)量均有顯著增加,其中ClHSP70和ClHSP90在處理0.5 h時(shí)表達(dá)量都開始明顯升高,分別為對(duì)照的5.8和5.1倍,1 h時(shí)表達(dá)量達(dá)到最高,分別為對(duì)照的12.2和9.7倍,2和4 h后ClHSP70和ClHSP90基因表達(dá)量有所下降(圖7)。說明ClHSP70和ClHSP90屬于熱誘導(dǎo)型表達(dá)的熱激蛋白,當(dāng)甘菊受到高溫脅迫時(shí),在很短時(shí)間內(nèi)就能夠參與熱脅迫響應(yīng)進(jìn)行大量表達(dá),增加植株體內(nèi)HSPs的積累,增強(qiáng)植株的耐熱性。

    根據(jù)ClHSP70和ClHSP90高溫處理不同時(shí)間的表達(dá)情況,將甘菊植株42 ℃處理1 h, 分別取植株的嫩葉、成熟葉、嫩莖、成熟莖、根提取RNA,實(shí)時(shí)熒光定量PCR分析結(jié)果表明,ClHSP70基因在成熟葉中的表達(dá)量顯著高于其他部位,其次是嫩葉、根、成熟莖、嫩莖。ClHSP90基因在成熟莖中的表達(dá)量顯著高于其他部位,嫩莖中表達(dá)最低,嫩葉、成熟葉、根之間的表達(dá)差異不明顯(圖7)。說明在高溫誘導(dǎo)下ClHSP70和ClHSP90的表達(dá)存在明顯的組織器官特異性,同時(shí)ClHSP70與ClHSP90之間表達(dá)水平最高的部位也不相同,推測(cè)上述2類熱激蛋白在植株中的定位和執(zhí)行功能有所不同。

    圖5 甘菊ClHSP70和ClHSP90保守域分析結(jié)果Fig. 5 Prediction of conserver domain of ClHSP70 and ClHSP90 from C. lavandulifolium

    節(jié)點(diǎn)數(shù)字表示Bootstrap驗(yàn)證中基于1 000次重復(fù)的可信度;括號(hào)內(nèi)編號(hào)為氨基酸登錄號(hào)圖6 熱激蛋白HSP70(A)和HSP90(B)基因編碼的氨基酸系統(tǒng)進(jìn)化樹Bootstrap values are shown at nodes based on 1 000 replications; The accession No. is given in bracketsFig. 6 Phylogenetic tree for HSP70 (A) and HSP90 (B) based on alignment of amino acid sequences from other species

    不同小寫字母表示同一基因不同處理時(shí)間和在不同組織中的差異顯著,P<0.05圖7 甘菊葉片中ClHSP70和ClHSP90 42 ℃高溫處理不同時(shí)間和在不同組織中的表達(dá)Different letters represent significant difference among different treatment times and different tissues under 42 ℃ (P<0.05)Fig. 7 Relative expression level of ClHSP70 and ClHSP90 gene in different treatment times and different tissues under 42 ℃ high temperature stress.

    3 討 論

    高溫是植物生命活動(dòng)中常會(huì)遭遇的非生物逆境,熱激蛋白(HSPs)是細(xì)胞受到高溫環(huán)境或其他脅迫環(huán)境誘導(dǎo)時(shí)產(chǎn)生的一類功能型蛋白質(zhì),可以作為分子伴侶參與細(xì)胞內(nèi)新合成蛋白質(zhì)的折疊、加工、轉(zhuǎn)運(yùn)以及蛋白質(zhì)變性后的復(fù)性[19],增強(qiáng)細(xì)胞對(duì)脅迫環(huán)境的忍耐力從而提高細(xì)胞的適應(yīng)能力[20]。菊花是中國(guó)十大名花,世界四大切花(菊花、月季、康乃馨、唐菖蒲)之一,產(chǎn)量居首,隨著國(guó)內(nèi)菊花市場(chǎng)規(guī)模的發(fā)展,其栽培面積在不斷擴(kuò)大,但夏季高溫影響了其產(chǎn)量和品質(zhì)。因此,研究熱激蛋白HSPs基因的結(jié)構(gòu)和功能,對(duì)于提高菊花耐熱性,培育耐高溫品種有長(zhǎng)遠(yuǎn)意義。

    本研究利用甘菊均一化轉(zhuǎn)錄組數(shù)據(jù),從甘菊葉片中克隆得到ClHSP70和ClHSP90基因,本實(shí)驗(yàn)中保守區(qū)預(yù)測(cè)分析表明,ClHSP70含有典型的NBD結(jié)構(gòu)域和SBD結(jié)構(gòu)域,其在HSP70與ATP結(jié)合-水解和發(fā)揮分子伴侶功能時(shí)起著重要作用。據(jù)報(bào)道,HSP90多肽N端結(jié)構(gòu)域包含1個(gè)ATP結(jié)合區(qū)域,在HSP90執(zhí)行生物功能時(shí)起著重要作用,真核生物細(xì)胞質(zhì)HSP90的C端含有保守的MEEVD基序[21],能夠與TPR(tetratricopeptide repeat)結(jié)構(gòu)域的蛋白識(shí)別并結(jié)合,從而使HSP90發(fā)揮不同的作用。本研究中ClHSP90 N端結(jié)構(gòu)域包含ATPase結(jié)合位點(diǎn),C端含有MEEVD基序,與萵苣[22]、番茄[23]和擬南芥[24]中的研究一致,是細(xì)胞質(zhì)HSP90。氨基酸序列比對(duì)分析表明:ClHSP70與大豆和煙草HSP70蛋白有較高的同源性,ClHSP90與紫莖澤蘭、煙草、番茄等植物HSP90有著很高的同源性。從系統(tǒng)進(jìn)化樹可以得出,ClHSP70與擬南芥、甘藍(lán)HSP70聚成一類,ClHSP90與同屬于菊科植物的紫莖澤蘭聚為一組,并與煙草、番茄HSP90聚成一大分支。因此,ClHSP70和ClHSP90分別屬于植物HSP70基因家族和HSP90基因家族成員。

    高溫脅迫下,熱激蛋白受到誘導(dǎo)且表達(dá)量明顯增加,同時(shí)在正常的環(huán)境下為維持細(xì)胞的正常運(yùn)轉(zhuǎn)HSP70和HSP90可以組成型表達(dá)。本研究中,ClHSP70和ClHSP90在25 ℃時(shí)均有表達(dá),且表達(dá)相對(duì)穩(wěn)定。研究表明,多數(shù)HSPs在受高溫脅迫后短時(shí)間內(nèi)就能誘導(dǎo)合成,熱激處理3~5 min內(nèi),HSPs mRNA的量增加,20 min時(shí)可檢測(cè)到新合成的HSPs, 隨著熱激時(shí)間的延長(zhǎng),HSP的合成能夠持續(xù)數(shù)小時(shí)。本研究中,甘菊ClHSP70和ClHSP90在42 ℃處理時(shí)其合成主要在1 h內(nèi)完成并達(dá)到最高,2 h后開始下降,4 h急劇下降。這與擬南芥中HSP70家族[25]和萵苣LsHsp90[22]在受到熱激處理時(shí)的表達(dá)模式類似。隨著熱激時(shí)間增加相對(duì)表達(dá)量下降的原因可能是因?yàn)闊峒せ虻拈L(zhǎng)時(shí)間表達(dá)對(duì)生物體不利,植物適應(yīng)高溫環(huán)境后將其表達(dá)下調(diào)[26]。根據(jù)亞細(xì)胞定位,植物HSP70主要分為四大類,分別為位于細(xì)胞質(zhì)HSP70、內(nèi)質(zhì)網(wǎng)HSP70、線粒體HSP70和葉綠體HSP70。HSP90主要定位于細(xì)胞質(zhì),也有少部分存在于內(nèi)質(zhì)網(wǎng)、線粒體和葉綠體等亞細(xì)胞中。研究表明,HSPs在熱脅迫下的表達(dá)與其定位以及執(zhí)行的功能有一定關(guān)系,即高溫脅迫下植物HSPs基因的表達(dá)有明顯的組織器官特異性。本實(shí)驗(yàn)中,42 ℃熱脅迫1 h,ClHSP70和ClHSP90在植株各個(gè)部位均有表達(dá),但表達(dá)量具有顯著差異性。其中ClHSP70在成熟葉中的表達(dá)量最高,與臘梅E-CpHSP70-1基因[27]、玉米熱激蛋白70[28]在各個(gè)組織器官中的表達(dá)分析一致,推測(cè)可能與其亞細(xì)胞定位主要分布在葉片中有關(guān),相對(duì)于嫩葉和其他組織部位,成熟葉細(xì)胞中葉綠體、內(nèi)質(zhì)網(wǎng)等細(xì)胞器含量更多。另外,當(dāng)植株受到高溫脅迫時(shí),首先會(huì)引起葉片萎蔫、葉色失綠等表型特征,ClHSP70在成熟葉片中的高表達(dá)表明其在植物響應(yīng)高溫脅迫和緩解植物體受到的熱損傷過程中起重要作用。ClHSP90在成熟莖中的表達(dá)量最高,在葉片中表達(dá)量較低,與劉云飛等在番茄中的研究結(jié)果類似[23],說明HSP90在植株受到熱脅迫后能在多個(gè)組織部位參與植物高溫脅迫的應(yīng)答。此外,HSP90除了參與高溫脅迫的應(yīng)答之外,其在植物生長(zhǎng)發(fā)育進(jìn)程中起著至關(guān)重要的作用,模式植物擬南芥[29-31]中的研究表明HSP90廣泛參與了植物的生長(zhǎng)發(fā)育過程。本實(shí)驗(yàn)中ClHSP90在成熟莖和嫩莖中的表達(dá)差異最為明顯,推測(cè)該基因可能參與莖尖分生組織相關(guān)調(diào)控活動(dòng)。甘菊熱激蛋白基因?qū)Ω邷孛{迫的表達(dá)模式研究為菊花耐熱性研究提供了分子基礎(chǔ),并為利用基因工程提高菊花耐熱性和耐熱種質(zhì)資源創(chuàng)新與應(yīng)用提供理論依據(jù)。關(guān)于ClHSP70和ClHSP90在菊花中如何行使其功能我們今后將做進(jìn)一步研究。

    [1]CHO E K, CHOI Y J. A nuclear-localized HSP70 confers thermoprotective activity and drought-stress tolerance on plants[J].BiotechnolLett, 2009, 31(4):597-606.

    [2]GUPTA S C, SHARMA A, MISHRA M. Heat shock proteins in toxicology: how close and how far[J].LifeSci, 2010, 86(11-12):377-384.

    [3]CHANKOVA S, MITROVSKA Z, MITEVA D. Heat shock protein HSP70B as a marker for genotype resistance to environmental stress in Chlorella species from contrasting habitats[J].Gene, 2013, 516(1):184-189.

    [4]SANGSTER T A, QUEITSCH C. The HSP90 chaperone complex, an emerging force in plant development and phenotypic plasticity[J].CurrOpinPlantBiol, 2005, 8(1):86-92.

    [5]DONG Y S, GUY C L. Physiological and molecular assessment of altered expression of HSP70-1 inArabidopsis. Evidence for pleiotropic,consequences[J].PlantPhysiology, 2003, 132(2):979-987.

    [6]HU X, LIU R, LI Y. Heat shock protein 70 regulates the abscisic acid-induced antioxidant response of maize to combined drought and heat stress[J].PlantGrowthRegulation, 2010, 60(3):225-235.

    [7]YOUNG J C, MOAREFI I, HARTL F U. HSP90: a specialized but essential protein-folding tool[J].JournalofCellBiology, 2001, 154(2):267-273.

    [8]TERASAWA K, MINAMI M, MINAMI Y. Constantly updated knowledge of Hsp90[J].JournalofBiochemistry, 2005, 137(4):443-447.

    [9]YABE N, TAKAHASHI T, KOMEDA Y. Analysis of tissue-specific expression ofArabidopsisthalianaHSP90-family gene HSP81[J].PlantCellPhysiol., 1994, 35(8):1 207-1 219.

    [10]LIU D, ZHANG X, CHENG Y. rHSP90 gene expression in response to several environmental stresses in rice (OryzasativaL.)[J].PlantPhysiol.Biochem., 2006, 44(5-6):380-386.

    [11]PICARD D. Heat-shock protein 90, a chaperone for folding and regulation[J].Cellular&MolecularLifeSciencesCmls, 2002, 59(10):1 640-1 648.

    [12]JACKSON S E, QUEITSCH C, TOFT D. HSP90: from structure to phenotype[J].NatureStructural&MolecularBiology, 2004, 11(12):1 152-1 155.

    [13]HAQ N U, RAZA S, LUTHE D S. A dual role for the chloroplast small heat shock protein ofChenopodiumalbum including protection from both heat and metal stress[J].PlantMolecularBiologyReporter, 2012, 31(2):398-408.

    [14]SIDDIQUE M, GERNHARD S, KOSKULL-D?RING P V. The plant sHSP superfamily: five new members inArabidopsisthalianawith unexpected properties[J].CellStress&Chaperones, 2007, 13(2):183-197.

    [15]DAI S L, WANG W K, LI M X. Phylogenetic relationship of dendranthema (DC.) Des Moul. revealed by fluorescentinsituhybridization[J].JournalofIntegrativePlantBiology, 2005, 47(7):783-791.

    [16]WANG W, DAI S, LI M. Physical mapping of rDNA inDendranthemanankingense and its close related species by fluorescent in situ hybridization[J].Cellular&MolecularBiologyLetters, 2002, 7(3):911-917.

    [17]CHENNA R, SUGAWARA H, KOIKE T. Multiple sequence alignment with the clustal series of programs[J].NucleicAcidsRes, 2003, 31(13):3 497-3 500.

    [18]TAMURA K, PETERSON D, PETERSON N. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods[J].Mol.Biol.Evol., 2011, 28(10):2 731-2 739.

    [19]HORWITZ J. Alpha-crystallin can function as a molecular chaperone[J].ProceedingsoftheNationalAcademyofSciencesoftheUnitedStatesofAmerica, 1992, 89(21):10 449-10 453.

    [20]BALLINGER D G, PARDUE M L. The control of protein synthesis during heat shock in drosophila cells involves altered polypeptide elongation rates[J].Cell, 1983, 33(33):103-113.

    [21]PRODROMOU C, ROE S M, O'BRIEN R. Identification and structural characterization of the ATP/ADP-binding site in the HSP90 molecular chaperone[J].Cell, 1997, 90(90):65-75.

    [22]任月,韓瑩琰,李婷,等。葉用萵苣熱激蛋白90(LsHSP90)基因的克隆及其在熱激下的表達(dá)[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2013, 46(16):3 514-3 522.

    REN Y, HAN Y Y, LI T,etal. Molecular cloning and expression analysis of heat-shock-protein90 (LsHSP90) gene from leaf lettuce (LactucasativaL.) under heat shock[J].ScientiaAgriculturaSinica, 2013, 46(16):3 514-3 522.

    [23]劉云飛, 萬紅建, 楊悅儉,等。番茄熱激蛋白90的全基因組鑒定及分析[J]. 遺傳,2014,36(10):1 043-1 052.

    LIU Y F, WAN H J, YANG Y J,etal. Genome-wide identification and analysis of heat shock protein 90 in tomato[J].Hereditas(Beijing), 2014,36(10):1 043-1 052.

    [24]MINAMI Y, KIMURA Y, KAWASAKI H. The carboxy-terminal region of mammalian HSP90 is required for its dimerization and functioninvivo[J].Molecular&CellularBiology, 1994, 14(2):1 459-1 464.

    [25]SUNG D Y, VIERLING E, GUY C L. Comprehensive expression profile analysis of theArabidopsisHSP70 gene family[J].PlantPhysiol., 2001, 126(2):789-800.

    [26]KREBS R A, HOLBROOK S H. Reduced enzyme activity following HSP70 overexpression in drosophila melanogaster[J].BiochemicalGenetics, 2001, 39(1-2):73-82.

    [27]阮文進(jìn), 門維婷, 馬婧,等。蠟梅熱激蛋白基因CpHSP70-1的克隆、亞細(xì)胞定位與表達(dá)分析[J]. 西南大學(xué)學(xué)報(bào)(自然科學(xué)版),2016, 38(1):43-52.

    TIEN V N, MEN W T, MA J,etal. Cloning, subcellular localization and expression analysis of heat shock protein geneCpHSP70-1 fromChimonanthuspraecox(L.)[J].JournalofSouthwestUniversity(Natural Science Edition), 2016,38(1):43-52.

    [28]李慧聰,郭秀林,王冬梅,等。玉米熱激蛋白70基因?qū)囟让{迫的響應(yīng)[J]. 河北農(nóng)業(yè)大學(xué)學(xué)報(bào), 2010, 33(6):12-15.

    LI H C, GUO X L, WANG D M,etal. Responses ofHSP70 gene expression to temperature stresses in maize (ZeamaysL.)[J].JournalofAgriculturalUniversityofHebei, 2010, 33(6):12-15.

    [29]QUEITSCH C, SANGSTER T A, LINDQUIST S. HSP90 as a capacitor phenotypic variation[J].Nature, 2002, 417(6 889): 618-624.

    [30]SANGSTER T A, BAHRAMI A, WILCZEK A. Phenotypic diversity and altered environmental plasticity inArabidopsisthalianawith reduced HSP90 levels[J].PlosOne, 2007, 2(7):e648.

    [31]PRASINOS C, KRAMPIS K, SAMAKOVLI D. Tight regulation of expression of twoArabidopsiscytosolicHSP90 genes during embryo development[J].JournalofExperimentalBotany, 2005, 56(412):633-644.

    (編輯:宋亞珍)

    Cloning and Expression ofClHSP70 andClHSP90 Gene fromChrysanthemumlavandulifolium

    LUO Chang, CHEN Dongliang, CHENG Xi, HUANG Conglin*

    (Beijing Academy of Agriculture and Forestry Sciences, Beijing Agricultural Biotechnology Research Center, Beijing 100097, China)

    Based onChrysanthemumlavandulifoliumRNA-seq data, we obtained two code DNA sequences namedClHSP70 andClHSP90 using RT-PCR. Sequence analysis showed the predicted open reading frame (ORF) ofClHSP70 is 2 559 bp, encoding a protein of 852 amino acids. The conserver domain search exhibited that N-teminar contains typical NBD domain which belongs to HSP70 family. The predicted ORF ofClHSP90 is 2 094 bp, encoding a protein of 697 amino acids. The conserver domain search showed that N-teminar contains ATPase-superfamily and HSP90 superfamily conservative structure. Bioinformatics showed that the amino acid sequences of ClHSP70 had high consistency withGlycinemaxandNicotianatomentosiformisHSP70, while ClHSP90 had high consistency withAgeratinaadenophoraHSP90. Real-time PCR demonstrated the expression patterns ofClHSP70 andClHSP90 in leaves at 42 ℃. The expression of these two genes were both significantly up-regulated after 0.5 h, reached the highest level after 1 h. Then the expression of both two genes decreased after 2 h, and remained at a relatively lower level after 4 h. In different plant tissues (root, stem, leaf) after 1 h at 42 ℃,ClHSP70 was highly expressed in mature leaves, lower in young, leaves, root, mature stems, young stems, andClHSP90 was the highest in mature stems. The results showed thatClHSP70 andClHSP90 participated the process of heat tolerance inChrysanthemumlavandulifoliumand laid a foundation for analyzingClHSP70 andClHSP90 gene function.

    Chrysanthemumlavandulifolium;ClHSP70;ClHSP90; heat stress

    1000-4025(2016)07-1321-10

    10.7606/j.issn.1000-4025.2016.07.1321

    2016-03-28;修改稿收到日期:2016-05-17

    北京市農(nóng)林科學(xué)院科技創(chuàng)新能力建設(shè)專項(xiàng)(KJCX20140109)

    羅昌(1982-),男,碩士,助理研究員,主要從事菊花遺傳育種研究。E-mail: changluo1983@126.com.

    黃叢林,研究員,主要從事菊花遺傳育種研究。E-mail: conglinh@126.com

    Q785;Q786

    A

    猜你喜歡
    甘菊擬南芥結(jié)構(gòu)域
    擬南芥:活得粗糙,才讓我有了上太空的資格
    薇甘菊生物防治的研究進(jìn)展
    薇甘菊光能利用及葉綠素合成在不同光照強(qiáng)度下的響應(yīng)
    薇甘菊的發(fā)生現(xiàn)狀及防治措施
    蛋白質(zhì)結(jié)構(gòu)域劃分方法及在線服務(wù)綜述
    可愛的德國(guó)甘菊
    尿黑酸對(duì)擬南芥酪氨酸降解缺陷突變體sscd1的影響
    兩種LED光源作為擬南芥生長(zhǎng)光源的應(yīng)用探究
    擬南芥干旱敏感突變體篩選及其干旱脅迫響應(yīng)機(jī)制探究
    重組綠豆BBI(6-33)結(jié)構(gòu)域的抗腫瘤作用分析
    中文字幕精品免费在线观看视频 | 国产极品粉嫩免费观看在线 | 精品国产露脸久久av麻豆| 亚洲精品国产av蜜桃| 久久久久精品性色| 日韩中字成人| 精品少妇黑人巨大在线播放| 亚洲美女搞黄在线观看| 18禁在线播放成人免费| 亚洲精品色激情综合| 亚洲色图 男人天堂 中文字幕 | 亚洲精品乱码久久久v下载方式| 亚洲成人手机| 国产片内射在线| 国产伦精品一区二区三区视频9| 成人毛片60女人毛片免费| 久久亚洲国产成人精品v| 香蕉精品网在线| av不卡在线播放| 男女高潮啪啪啪动态图| 精品国产乱码久久久久久小说| av.在线天堂| 欧美人与善性xxx| 久久狼人影院| 亚洲第一av免费看| 老熟女久久久| 一级黄片播放器| 91在线精品国自产拍蜜月| 黄色怎么调成土黄色| 亚州av有码| 哪个播放器可以免费观看大片| 制服丝袜香蕉在线| 男女啪啪激烈高潮av片| 如日韩欧美国产精品一区二区三区 | 一本久久精品| 免费播放大片免费观看视频在线观看| 精品熟女少妇av免费看| 亚洲国产日韩一区二区| 热99久久久久精品小说推荐| 亚洲国产成人一精品久久久| 狂野欧美激情性bbbbbb| 高清欧美精品videossex| av一本久久久久| 亚洲精品乱久久久久久| 99热国产这里只有精品6| 一区在线观看完整版| 国产精品一国产av| 男女啪啪激烈高潮av片| 亚洲国产精品一区二区三区在线| 岛国毛片在线播放| 人人妻人人爽人人添夜夜欢视频| 考比视频在线观看| 国产亚洲精品第一综合不卡 | 国产成人aa在线观看| 亚洲婷婷狠狠爱综合网| 久久综合国产亚洲精品| av又黄又爽大尺度在线免费看| 国产一区二区三区av在线| 欧美精品一区二区免费开放| 欧美人与善性xxx| 国产色爽女视频免费观看| 女性被躁到高潮视频| videossex国产| 亚洲精品久久午夜乱码| 久久ye,这里只有精品| 日本猛色少妇xxxxx猛交久久| 又粗又硬又长又爽又黄的视频| 少妇丰满av| 街头女战士在线观看网站| 亚洲国产精品一区二区三区在线| 大香蕉久久网| 考比视频在线观看| 少妇精品久久久久久久| 亚洲国产成人一精品久久久| 狂野欧美白嫩少妇大欣赏| 少妇的逼好多水| freevideosex欧美| 美女脱内裤让男人舔精品视频| 国产毛片在线视频| 极品少妇高潮喷水抽搐| 久热久热在线精品观看| 99热这里只有精品一区| 亚洲人成网站在线播| 一区二区日韩欧美中文字幕 | 黑人高潮一二区| 少妇人妻 视频| 久久精品国产自在天天线| 亚洲一级一片aⅴ在线观看| 91国产中文字幕| 最新的欧美精品一区二区| 22中文网久久字幕| 久久99蜜桃精品久久| av在线观看视频网站免费| 国产黄色免费在线视频| 美女大奶头黄色视频| 色94色欧美一区二区| 成年美女黄网站色视频大全免费 | 日韩大片免费观看网站| 又大又黄又爽视频免费| 一级毛片 在线播放| 亚洲欧洲国产日韩| 日韩制服骚丝袜av| 搡女人真爽免费视频火全软件| 99热国产这里只有精品6| 亚洲怡红院男人天堂| 视频区图区小说| av福利片在线| 男人操女人黄网站| 成人国产麻豆网| 人人妻人人澡人人看| 久久久久久久久大av| 成人影院久久| 黑人猛操日本美女一级片| 又粗又硬又长又爽又黄的视频| 成人18禁高潮啪啪吃奶动态图 | 日韩成人av中文字幕在线观看| 免费大片18禁| 亚洲,一卡二卡三卡| 午夜老司机福利剧场| 免费大片18禁| 有码 亚洲区| 蜜臀久久99精品久久宅男| 国产欧美亚洲国产| 亚洲三级黄色毛片| 精品一品国产午夜福利视频| 少妇的逼水好多| 午夜福利,免费看| 成人亚洲欧美一区二区av| 成人亚洲精品一区在线观看| 亚洲精品aⅴ在线观看| av播播在线观看一区| 少妇人妻 视频| 亚洲内射少妇av| 一二三四中文在线观看免费高清| 午夜激情福利司机影院| av黄色大香蕉| 青春草视频在线免费观看| 天天影视国产精品| 国产又色又爽无遮挡免| 亚洲,欧美,日韩| 51国产日韩欧美| xxx大片免费视频| 成人18禁高潮啪啪吃奶动态图 | 国产免费一区二区三区四区乱码| 久久久久精品性色| 亚洲av综合色区一区| 久久精品国产a三级三级三级| 亚洲图色成人| 18+在线观看网站| 久久女婷五月综合色啪小说| 人人妻人人澡人人看| 99九九在线精品视频| 特大巨黑吊av在线直播| 免费看光身美女| av播播在线观看一区| 国产高清国产精品国产三级| av国产久精品久网站免费入址| 春色校园在线视频观看| 少妇熟女欧美另类| 欧美 日韩 精品 国产| 日韩视频在线欧美| 色网站视频免费| 亚洲精品自拍成人| 欧美日韩精品成人综合77777| 视频区图区小说| 自拍欧美九色日韩亚洲蝌蚪91| 一级a做视频免费观看| 最近中文字幕高清免费大全6| 91精品国产九色| 99re6热这里在线精品视频| 成年人免费黄色播放视频| 看免费成人av毛片| 大陆偷拍与自拍| 亚洲av国产av综合av卡| 国产一区二区在线观看av| 日本av手机在线免费观看| 在线免费观看不下载黄p国产| 亚洲不卡免费看| 国产成人freesex在线| 久久热精品热| 自拍欧美九色日韩亚洲蝌蚪91| 国产片内射在线| 色视频在线一区二区三区| 丝袜脚勾引网站| 国产无遮挡羞羞视频在线观看| 97在线人人人人妻| 赤兔流量卡办理| 黄色配什么色好看| 两个人的视频大全免费| 亚洲天堂av无毛| 18禁在线播放成人免费| 下体分泌物呈黄色| 日日啪夜夜爽| 高清午夜精品一区二区三区| 亚洲av男天堂| 成人18禁高潮啪啪吃奶动态图 | 涩涩av久久男人的天堂| 亚洲国产精品999| 交换朋友夫妻互换小说| 七月丁香在线播放| 观看美女的网站| 99热全是精品| 国产不卡av网站在线观看| 2022亚洲国产成人精品| 中国三级夫妇交换| 午夜91福利影院| 国产精品久久久久久久久免| 国产精品99久久久久久久久| 国产女主播在线喷水免费视频网站| av免费观看日本| 国产男人的电影天堂91| 国产免费一区二区三区四区乱码| 青春草国产在线视频| 精品一品国产午夜福利视频| 满18在线观看网站| 大香蕉久久成人网| 亚洲精品色激情综合| 2018国产大陆天天弄谢| 91精品国产国语对白视频| 欧美日韩av久久| 丝瓜视频免费看黄片| 人妻制服诱惑在线中文字幕| 只有这里有精品99| 国产精品一区二区在线不卡| 国产精品熟女久久久久浪| 国产视频内射| 欧美精品高潮呻吟av久久| 五月开心婷婷网| 日韩一本色道免费dvd| 十分钟在线观看高清视频www| 在线观看国产h片| 狂野欧美白嫩少妇大欣赏| 新久久久久国产一级毛片| 国产男女超爽视频在线观看| 少妇丰满av| 黄色怎么调成土黄色| 亚洲av免费高清在线观看| 超碰97精品在线观看| 久久久欧美国产精品| 国产av精品麻豆| 99久久中文字幕三级久久日本| 在现免费观看毛片| 欧美精品高潮呻吟av久久| 欧美激情 高清一区二区三区| 亚洲av二区三区四区| 久久热精品热| 日本wwww免费看| 午夜福利视频精品| 亚洲精品成人av观看孕妇| 丁香六月天网| 在线精品无人区一区二区三| 亚洲欧美一区二区三区黑人 | 性色avwww在线观看| 欧美 日韩 精品 国产| 成人午夜精彩视频在线观看| 国产高清不卡午夜福利| 这个男人来自地球电影免费观看 | 七月丁香在线播放| 视频中文字幕在线观看| 欧美 亚洲 国产 日韩一| 欧美日韩亚洲高清精品| 久热这里只有精品99| 国产高清不卡午夜福利| 蜜臀久久99精品久久宅男| 男女高潮啪啪啪动态图| 99九九在线精品视频| 性色av一级| 99久久精品一区二区三区| 女性生殖器流出的白浆| 亚洲精品自拍成人| 久久久久视频综合| 日日爽夜夜爽网站| 亚洲人与动物交配视频| 新久久久久国产一级毛片| 欧美激情国产日韩精品一区| 永久免费av网站大全| 日韩制服骚丝袜av| 黄片播放在线免费| 黄色视频在线播放观看不卡| 亚洲欧美色中文字幕在线| 99精国产麻豆久久婷婷| 日本色播在线视频| 亚洲内射少妇av| 一区二区三区免费毛片| 搡女人真爽免费视频火全软件| 中文欧美无线码| 十分钟在线观看高清视频www| 中文字幕最新亚洲高清| 成人无遮挡网站| 精品久久久噜噜| 我的老师免费观看完整版| 国产精品秋霞免费鲁丝片| 男女边吃奶边做爰视频| 韩国av在线不卡| 欧美+日韩+精品| 午夜av观看不卡| 免费观看无遮挡的男女| 特大巨黑吊av在线直播| 日韩在线高清观看一区二区三区| 亚洲成色77777| 美女大奶头黄色视频| 欧美日韩精品成人综合77777| av卡一久久| 老女人水多毛片| 永久网站在线| 免费av中文字幕在线| 亚洲欧美中文字幕日韩二区| 欧美97在线视频| 成人亚洲精品一区在线观看| 黑人高潮一二区| 在线观看www视频免费| videos熟女内射| 男的添女的下面高潮视频| 免费不卡的大黄色大毛片视频在线观看| 免费日韩欧美在线观看| 国产色爽女视频免费观看| 美女脱内裤让男人舔精品视频| 免费播放大片免费观看视频在线观看| 水蜜桃什么品种好| a级毛色黄片| 久久久欧美国产精品| 精品人妻偷拍中文字幕| 男女高潮啪啪啪动态图| 午夜福利影视在线免费观看| 天天影视国产精品| 制服丝袜香蕉在线| av又黄又爽大尺度在线免费看| 人妻人人澡人人爽人人| 久久久久久久大尺度免费视频| 在线精品无人区一区二区三| 午夜日本视频在线| 一边亲一边摸免费视频| 午夜91福利影院| 亚洲av在线观看美女高潮| 黑丝袜美女国产一区| 亚洲怡红院男人天堂| 毛片一级片免费看久久久久| 91精品国产九色| 熟女av电影| 久久久久久久久久久久大奶| 在线观看三级黄色| 久久久久国产精品人妻一区二区| 久久国产精品男人的天堂亚洲 | 黑人巨大精品欧美一区二区蜜桃 | 国产精品99久久久久久久久| 日韩人妻高清精品专区| √禁漫天堂资源中文www| 青春草国产在线视频| 精品少妇黑人巨大在线播放| 日本色播在线视频| 成人无遮挡网站| 黄色欧美视频在线观看| 亚洲第一区二区三区不卡| 久久久久视频综合| 亚洲精品国产av成人精品| 欧美97在线视频| 纵有疾风起免费观看全集完整版| 又粗又硬又长又爽又黄的视频| 国产精品一区www在线观看| 97超碰精品成人国产| 免费人妻精品一区二区三区视频| 亚洲国产精品一区三区| 久久99热6这里只有精品| 人妻一区二区av| 国产精品一国产av| 亚洲精品av麻豆狂野| 黑人巨大精品欧美一区二区蜜桃 | 亚洲性久久影院| 国产伦精品一区二区三区视频9| 亚洲婷婷狠狠爱综合网| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品亚洲一区二区| 丰满乱子伦码专区| 三级国产精品片| videosex国产| 国产片内射在线| 国产高清三级在线| 人妻夜夜爽99麻豆av| 国产片内射在线| 日本黄色片子视频| 一个人看视频在线观看www免费| 亚洲精品亚洲一区二区| 国产精品女同一区二区软件| 亚洲av综合色区一区| videosex国产| 18禁裸乳无遮挡动漫免费视频| 午夜激情av网站| 熟女av电影| 最新中文字幕久久久久| 欧美日韩国产mv在线观看视频| 久久婷婷青草| 欧美日韩视频精品一区| av在线播放精品| 涩涩av久久男人的天堂| 亚洲av日韩在线播放| 狂野欧美激情性xxxx在线观看| 婷婷色av中文字幕| 伊人久久精品亚洲午夜| 国产 一区精品| 国产高清有码在线观看视频| 极品少妇高潮喷水抽搐| 永久网站在线| 在线观看美女被高潮喷水网站| 亚洲图色成人| 亚洲欧美一区二区三区黑人 | 日韩av免费高清视频| 日韩亚洲欧美综合| 精品人妻熟女毛片av久久网站| 观看av在线不卡| 午夜91福利影院| 精品视频人人做人人爽| 免费观看a级毛片全部| 最后的刺客免费高清国语| h视频一区二区三区| 亚洲av成人精品一二三区| 狂野欧美激情性bbbbbb| www.av在线官网国产| 亚洲一级一片aⅴ在线观看| 色婷婷久久久亚洲欧美| 精品人妻在线不人妻| 日韩av在线免费看完整版不卡| 国产成人a∨麻豆精品| 80岁老熟妇乱子伦牲交| 制服诱惑二区| 中文字幕精品免费在线观看视频 | 亚洲国产欧美日韩在线播放| 狠狠婷婷综合久久久久久88av| 久久影院123| 久久久久久久亚洲中文字幕| 精品人妻偷拍中文字幕| 久久99热这里只频精品6学生| 国产在线一区二区三区精| 夜夜看夜夜爽夜夜摸| 男女边吃奶边做爰视频| 亚洲国产欧美日韩在线播放| 九九在线视频观看精品| 黑人猛操日本美女一级片| 飞空精品影院首页| 色网站视频免费| 亚洲伊人久久精品综合| 婷婷色综合www| 久久久久久久久久久免费av| 九色成人免费人妻av| 久久久久久久大尺度免费视频| 日本黄色日本黄色录像| xxx大片免费视频| 麻豆精品久久久久久蜜桃| 人妻夜夜爽99麻豆av| 午夜福利视频精品| 热99久久久久精品小说推荐| 一本—道久久a久久精品蜜桃钙片| 性色av一级| 两个人的视频大全免费| 亚洲婷婷狠狠爱综合网| 国产成人午夜福利电影在线观看| 亚洲精品,欧美精品| 国产老妇伦熟女老妇高清| 大香蕉久久网| 久久影院123| 黄色欧美视频在线观看| 亚洲国产av影院在线观看| 亚洲欧美一区二区三区国产| 黄色一级大片看看| av一本久久久久| 亚洲图色成人| 国产午夜精品久久久久久一区二区三区| 国语对白做爰xxxⅹ性视频网站| 精品少妇黑人巨大在线播放| xxx大片免费视频| 涩涩av久久男人的天堂| 国产精品国产三级专区第一集| 久久久久久久精品精品| 99久国产av精品国产电影| 一区二区三区精品91| 王馨瑶露胸无遮挡在线观看| 国产午夜精品一二区理论片| 国产精品秋霞免费鲁丝片| 国产黄色视频一区二区在线观看| www.av在线官网国产| 亚洲欧洲日产国产| 18禁观看日本| 天天躁夜夜躁狠狠久久av| 日日摸夜夜添夜夜爱| 久久这里有精品视频免费| www.色视频.com| 热re99久久精品国产66热6| 蜜臀久久99精品久久宅男| 亚洲欧美清纯卡通| 久久精品国产亚洲av涩爱| 亚洲欧洲日产国产| 极品少妇高潮喷水抽搐| 91aial.com中文字幕在线观看| 热99国产精品久久久久久7| 精品酒店卫生间| 欧美国产精品一级二级三级| 97在线人人人人妻| 哪个播放器可以免费观看大片| 国产成人aa在线观看| 天堂中文最新版在线下载| 国产精品久久久久久av不卡| 亚洲在久久综合| 成人影院久久| 日本欧美视频一区| 极品人妻少妇av视频| 另类精品久久| 婷婷色麻豆天堂久久| 免费观看的影片在线观看| 久久精品人人爽人人爽视色| 黑人欧美特级aaaaaa片| 午夜免费观看性视频| 久久久国产一区二区| 丝袜喷水一区| 国产伦精品一区二区三区视频9| 美女视频免费永久观看网站| 男的添女的下面高潮视频| 精品国产国语对白av| 九草在线视频观看| 伦理电影免费视频| 精品亚洲乱码少妇综合久久| 久久精品国产亚洲网站| 午夜影院在线不卡| 免费不卡的大黄色大毛片视频在线观看| 免费黄频网站在线观看国产| 日日撸夜夜添| 国产日韩欧美亚洲二区| 亚洲怡红院男人天堂| 久久久久久久大尺度免费视频| 中文字幕人妻丝袜制服| 成人影院久久| av专区在线播放| 亚洲久久久国产精品| 午夜久久久在线观看| 亚洲国产毛片av蜜桃av| 少妇人妻 视频| 香蕉精品网在线| 久久99精品国语久久久| 久久久久精品性色| 寂寞人妻少妇视频99o| 欧美日韩国产mv在线观看视频| 女的被弄到高潮叫床怎么办| 亚洲精品日韩av片在线观看| 亚洲国产色片| 99国产综合亚洲精品| 亚洲高清免费不卡视频| av有码第一页| 国产高清有码在线观看视频| 久久av网站| 搡老乐熟女国产| 人妻 亚洲 视频| 国产男人的电影天堂91| 国产色爽女视频免费观看| 日本vs欧美在线观看视频| 久久精品国产a三级三级三级| 欧美日韩精品成人综合77777| 成人二区视频| 久久免费观看电影| 2021少妇久久久久久久久久久| 成年美女黄网站色视频大全免费 | 嘟嘟电影网在线观看| 亚洲精品视频女| 欧美xxxx性猛交bbbb| 人妻夜夜爽99麻豆av| 大片免费播放器 马上看| 久久精品人人爽人人爽视色| 久久国产亚洲av麻豆专区| 日韩,欧美,国产一区二区三区| 美女xxoo啪啪120秒动态图| 少妇被粗大的猛进出69影院 | 亚洲国产日韩一区二区| 欧美日韩视频高清一区二区三区二| 亚洲精华国产精华液的使用体验| 亚洲内射少妇av| 久久婷婷青草| 九色成人免费人妻av| 亚洲国产毛片av蜜桃av| 中国三级夫妇交换| 97在线人人人人妻| 久久久久久久久大av| 97在线人人人人妻| 亚洲国产精品专区欧美| 99久久人妻综合| 精品国产一区二区久久| 在线亚洲精品国产二区图片欧美 | 亚洲精品国产av成人精品| 色网站视频免费| 亚洲国产最新在线播放| 久久久久久久久久人人人人人人| 91精品伊人久久大香线蕉| 久久鲁丝午夜福利片| 18禁在线无遮挡免费观看视频| 亚洲欧美精品自产自拍| 最近2019中文字幕mv第一页| 男的添女的下面高潮视频| 国产片特级美女逼逼视频| 人妻夜夜爽99麻豆av| 在现免费观看毛片| 午夜福利网站1000一区二区三区| 成年人午夜在线观看视频| 国产免费一级a男人的天堂| 欧美日韩av久久| 22中文网久久字幕| 狂野欧美激情性bbbbbb| 男女啪啪激烈高潮av片| 97精品久久久久久久久久精品| 午夜日本视频在线| 麻豆乱淫一区二区| 91国产中文字幕| 欧美 日韩 精品 国产| 99热国产这里只有精品6| 久久精品国产亚洲av天美| 十分钟在线观看高清视频www| 午夜福利视频在线观看免费| 亚洲国产欧美在线一区| 一区二区日韩欧美中文字幕 | a级毛片在线看网站|