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

    番茄SlMAPK6基因克隆及其表達(dá)特性分析

    2020-11-09 03:09:31岳寧波李云洲李玉龍潘鵬程潘寅濤王勇須文吳浪閆見敏
    南方農(nóng)業(yè)學(xué)報 2020年7期
    關(guān)鍵詞:非生物脅迫基因克隆番茄

    岳寧波 李云洲 李玉龍 潘鵬程 潘寅濤 王勇 須文 吳浪 閆見敏

    摘要:【目的】克隆番茄促分裂原活化蛋白激酶(MAPK)基因SlMAPK6,并分析其在不同非生物脅迫及信號物質(zhì)處理下的表達(dá)模式,為深入探究SlMAPK6基因在番茄逆境響應(yīng)中調(diào)控機(jī)制提供理論參考?!痉椒ā坎捎猛纯寺〖夹g(shù)從矮番茄中克隆SlMAPK6基因cDNA全長序列,并進(jìn)行生物信息學(xué)分析,采用實時熒光定量PCR檢測矮番茄幼苗在干旱(20% PEG-6000)、鹽(300 mmol/L NaCl)及信號物質(zhì)[0.50 mmol/L水楊酸(SA)和0.10 mmol/L 2,1,3-苯并噻二唑(BTH)]處理下SlMAPK6基因的表達(dá)模式?!窘Y(jié)果】克隆獲得SlMAPK6基因cDNA全長序列為1360 bp,開放閱讀框(ORF)為1131 bp,編碼376個氨基酸殘基,其編碼蛋白理論等電點(diǎn)(pI)6.76,不穩(wěn)定指數(shù)44.36,脂肪系數(shù)90.00,為不穩(wěn)定的強(qiáng)脂溶性蛋白,親水性平均值(GRAVY)-0.396,親水性氨基酸均勻分布在整個肽鏈中,故該蛋白為親水性蛋白;亞細(xì)胞定位于細(xì)胞核及胞漿質(zhì)中,無跨膜區(qū),含有36個磷酸化位點(diǎn),其中絲氨酸(Ser)磷酸位點(diǎn)19個、蘇氨酸(Thr)磷酸位點(diǎn)9個。SlMAPK6蛋白二級結(jié)構(gòu)中,α-螺旋占38.03%,延伸鏈占10.37%,無規(guī)則卷曲占51.60%。SlMAPK6蛋白與馬鈴薯StMAPK4、風(fēng)鈴辣椒CbMAPK4和中華辣椒CcMAPK4蛋白同源性均為99.47%。SlMAPK6和SlMAPK5蛋白的親緣關(guān)系最近,與擬南芥MAPK家族B亞族成員AtMAPK4、AtMPK5、AtMPK11、AtMPK12和AtMPK13蛋白的親緣關(guān)系也較近。經(jīng)SA、BTH、NaCl和PEG-6000處理后,SlMAPK6基因均出現(xiàn)抑制表達(dá),說明SlMAPK6基因負(fù)向調(diào)控番茄植株的逆境響應(yīng)。【結(jié)論】SlMAPK6基因可能參與負(fù)調(diào)控番茄植株的抗病和抗逆脅迫機(jī)制。SlMAPK6蛋白通過Ser或Thr磷酸化識別細(xì)胞信號的傳導(dǎo)功能,在矮番茄逆境脅迫響應(yīng)中發(fā)揮重要作用。

    關(guān)鍵詞: 番茄;促分裂原活化蛋白激酶6(MAPK6);基因克隆;信號物質(zhì);非生物脅迫;表達(dá)特性

    中圖分類號: S641.203.6? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼: A 文章編號:2095-1191(2020)07-1625-09

    Abstract:【Objective】 To clone the tomato mitogen-activated protein kinase(MAPK) gene SlMAPK6, to analyze its expression level under different abiotic stresses and signal substances, which? provided theoretical reference for exploration of the regulatory mechanism of the SlMAPK6 gene in tomato stress response. 【Method】Using homologous cloning technology to clone the full-length cDNA sequence of SlMAPK6 gene from dwarf tomato, conducting bioinformatics analysis, and real-time fluorescent quantitative PCR was used to detect the expression pattern of SlMAPK6 gene under the drought(20% PEG-6000), salt(300 mmol/L NaCl), signal substances[0.50 mmol/L salicylic acid (SA) and 0.10 mmol/L 2,1,3-benzothiadiazole(BTH)] treatments. 【Result】The full-length cDNA of SlMAPK6 was 1360 bp and contained an open reading frame(ORF) of 1131 bp that encoded 376 amino acids residues, the theoretical isoelectric point(pI) was 6.76, instability index? 44.36, and the fat coefficient 90.00,it was an unstable strong fat-soluble protein. The average hydrophilicity(GRAVY) was -0.396, hydrophilic amino acids were evenly distributed in the entire peptide chain, so the protein was a hydrophilic protein, subcellularly located in the nucleus and cytoplasm, without a transmembrane region.It contained 36 phosphorylation sites, including 19 serine(Ser) phosphate sites and 9 threonine(Thr) phosphate sites. In the secondary structure of SlMAPK6 protein, α-helix accounted for 38.03%, extended chain accounted for 10.37%, and random coil accounted for 51.60%. The homology of SlMAPK6 protein with potato StMAPK4, capsicum CbMAPK4 and Chinese capsicum CcMAPK4 were all 99.47%. SlMAPK6 and SlMAPK5 had the closest relationship, and it was also closely related to members of the arabidopsis MAPK family B subfamily: AtMAPK4, AtMPK5, AtMPK11, AtMPK12 and AtMPK13 proteins. After SA, BTH, NaCl and PEG-6000 treatments, the expression of SlMAPK6 gene was suppressed, indicating that SlMAPK6 gene negatively regulated the stress response of tomato plants. 【Conclusion】The SlMAPK6 gene may be involved in negatively regulating the disease resistance and stress resistance of tomato plants. The SlMAPK6 protein recognizes the transduction function of cell signals through Ser or Thr phosphorylation, and plays an important role in the response of dwarf tomato to stress.

    Key words: tomato; mitogen-activated protein kinase 6(MAPK6);gene cloning; signal substance response; abio-tic stress; expression characters

    Foundation item: National Natural Science Foundation of China(31960604); Guizhou Science and Technology Plan Project(QKHPTRC〔2017〕5788-28); Research Project on Talent Introduction of Guizhou University(KDRJHZ〔2017〕50); Laboratory Open Project of Guizhou University(SYSKF2019-55)

    0 引言

    【研究意義】番茄(Solanum lycopersicum L.)為茄科番茄屬一年生或多年生蔬菜作物,是我國重要的園藝蔬菜,其富含維生素、蘿卜素及具有抗氧化作用的番茄紅素等營養(yǎng)物質(zhì),可減緩衰老和預(yù)防癌癥(Perveen et al.,2015;Petyaev,2016;Crowe-White et al.,2019)。番茄在生長過程中常受到不利生長環(huán)境的影響,包括病原微生物侵染、干旱及鹽等生物和非生物脅迫,對番茄生產(chǎn)造成嚴(yán)重危害(李翠,2014)。經(jīng)研究發(fā)現(xiàn),促分裂原活化蛋白激酶(Mitogen-activated protein kinases,MAPK)在植物的生長、發(fā)育和脅迫應(yīng)答中發(fā)揮重要作用(Xu and Zhang,2015),因此,克隆番茄MAPK基因并分析其表達(dá)特性,深入探究其在番茄響應(yīng)生物和非生物脅迫中的調(diào)控機(jī)制,對提高番茄植株抗逆境脅迫能力具有重要意義?!厩叭搜芯窟M(jìn)展】MAPK在真核生物中普遍存在且高度保守,是一類絲氨酸(Ser)/蘇氨酸(Thr)蛋白激酶,能將外源或內(nèi)源信號分子導(dǎo)入細(xì)胞內(nèi)(Smekalova et al.,2013;Lee et al.,2016)。MAPK級聯(lián)信號傳導(dǎo)由MAPK、MAPK激酶(MAPKK)和MAPKK激酶(MAPKKK)3級信號傳導(dǎo)組成(de Zelicourt et al.,2016)。MAPK是末端級聯(lián)反應(yīng)信號激酶,通過上游MAPKK和MAPKKK對Ser和Thr殘基進(jìn)行磷酸化來傳遞信號(Chen et al.,2001;Hame et al.,2006;Rodriguez et al.,2010),從而參與植物的信號轉(zhuǎn)導(dǎo)應(yīng)答過程(Kalapos et al.,2019)。目前,已有大量的MAPK家族成員從不同植物物種中被鑒定出來,其中擬南芥、水稻、玉米和小麥分別鑒定出20、17、21和15個MAPK蛋白(Lian et al.,2012;Liu et al.,2013;Min and Zhao,2018)。番茄中已鑒定出16個MAPK蛋白(Wang et al.,2018a,2018b)。根據(jù)MAPK結(jié)構(gòu)差異,可將其分為A、B、C和D 4個亞族(Mohanta et al.,2015;王潔等,2016),其中,A、B和C亞族含TEY磷酸化位點(diǎn),D亞族含TDY磷酸化位點(diǎn)。在擬南芥中,A亞族以AtMPK3、AtMPK6和AtMPK10等為代表,主要與環(huán)境和激素反應(yīng)相關(guān),其中,AtMAPK3和AtMPK6在植株發(fā)育過程中發(fā)揮協(xié)同作用,MAPK3突變體和MAPK6突變體表現(xiàn)為胚胎異常,且AtMAPK3和AtMAPK6還參與脫落酸(ABA)、茉莉酸(JA)等激素響應(yīng)(Ortiz-Masia et al.,2007;Mine et al.,2017);B亞族以AtMPK4、AtMPK5、AtMPK11、AtMPK12和AtMPK13等為代表,主要參與細(xì)胞分裂和環(huán)境應(yīng)激反應(yīng),其中,AtMPK4還參與調(diào)節(jié)先天免疫(Ichimura et al.,2000;Petersen et al.,2000;Kosetsu et al.,2011);C亞族以AtMPK1、AtMPK2、AtMPK7和AtMPK14等為代表,主要參與應(yīng)激反應(yīng);D亞族包含AtMPK8、AtMPK9、AtMPK15、AtMPK16、AtMPK17、AtMPK18、AtMPK19和AtMPK20,主要參與逆境脅迫響應(yīng)和生長發(fā)育(孔福苓,2012)??梢?,各MAPK亞族均有成員參與非生物和生物應(yīng)激反應(yīng)(Wang et al.,2018a)?!颈狙芯壳腥朦c(diǎn)】在番茄MAPK家族中,以A亞族如SlMAPK1、SlMAPK2和SlMAPK3研究相對較多(Singh et al.,2018;Zhang et al.,2018),而B亞族如SlMAPK6(Kong et al.,2012)等研究甚少,其功能尚不清楚,至今鮮見有關(guān)SlMAPK6基因克隆及表達(dá)的研究報道?!緮M解決的關(guān)鍵問題】采用同源克隆技術(shù)從矮番茄中克隆SlMAPK6基因cDNA全長序列,并進(jìn)行生物信息學(xué)分析,采用實時熒光定量PCR檢測矮番茄幼苗在干旱、鹽脅迫及信號物質(zhì)處理下SlMAPK6基因的表達(dá)模式,為深入探究SlMAPK6基因在番茄逆境中的響應(yīng)調(diào)控機(jī)制提供理論參考。

    1 材料與方法

    1. 1 試驗材料

    以矮番茄為供試材料,引自西北農(nóng)林科學(xué)大學(xué)園藝學(xué)院,現(xiàn)保存于貴州大學(xué)農(nóng)學(xué)院番茄種質(zhì)資源庫。RNA isolater Total RNA Extraction Reagent細(xì)胞/組織中總RNA提取試劑盒、HiScript II 1st Strand cDNA Synthesis Kit試劑盒、Taq DNA聚合酶和SYBR Master Mix(2×)購自南京諾維贊生物科技有限公司;水楊酸(SA)、PEG-6000和NaCl購自生工生物工程(上海)股份有限公司;2,1,3-苯并噻二唑(BTH)購自上海源葉生物科技有限公司。

    主要儀器設(shè)備:人工智能光照培養(yǎng)箱(寧波市科技園區(qū)新江南儀器有限公司)、CFX96TM Real-time System熒光定量PCR儀(美國Bio-Rad公司)、T100TM PCR儀(美國Bio-Rad公司)、Tannon3500凝膠成像儀(上海天能科技有限公司)、低溫離心機(jī)(德國Eppendorf公司)和NanoDrop One微量核酸蛋白濃度測定儀[賽默飛世爾科技(中國)有限公司]。

    1. 2 樣品脅迫處理

    挑選飽滿的矮番茄種子,經(jīng)1.0% NaClO消毒15 min后,用無菌水沖洗殘留的NaClO,播種在裝有滅菌土壤的營養(yǎng)缽中,置于人工智能光照培養(yǎng)箱進(jìn)行育苗培養(yǎng)[光照強(qiáng)度240 mol/(m?s),光/暗條件16 h 28 ℃/8 h 18 ℃,相對濕度70%],待苗長至5片真葉時進(jìn)行非生物脅迫(干旱、鹽脅迫)及信號物質(zhì)處理。其中,干旱和鹽脅迫處理,分別采用20% PEG-6000溶液和300 mmol/L NaCl溶液澆灌幼苗根部,以無菌水澆灌的幼苗根部為對照;信號物質(zhì)處理:分別用0.10 mmol/L BTH和0.50 mmol/L SA噴灑幼苗葉面,以無菌水噴灑幼苗葉面為對照。每處理設(shè)3次重復(fù),分別于0、1、3、6、12和24 h采集0.50 g葉片組織,用液氮速凍后,置于-80 ℃保存?zhèn)溆谩?/p>

    1. 3 總RNA提取

    采用RNA Isolater Total RNA Extraction Reagent細(xì)胞/組織中總RNA提取試劑盒提取矮番茄葉片總RNA,分別用1.0%瓊脂糖凝膠電泳和核酸濃度檢測儀檢測RNA完整性和濃度,采用HiScript II 1st Strand cDNA Synthesis Kit試劑盒以DNA酶處理的1.0 μg總RNA為模板反轉(zhuǎn)錄合成cDNA,-80 ℃保存?zhèn)溆谩?/p>

    1. 4 基因克隆

    根據(jù)Sol Genomics Network網(wǎng)站公布的SlMAPK6基因序列(登錄號Solyc05G049970),利用Primer Premier 5.0設(shè)計其全長序列擴(kuò)增引物SlMAPK6-F和SlMAPK6-R(表1),委托生工生物工程(上海)股份有限公司合成。以cDNA為模板,用Taq DNA聚合酶克隆SlMAPK6基因cDNA全長序列。反應(yīng)體系20.0 μL:2×Taq PCR Master Mix 10.0 μL,10 μmol/L正、反向引物各1.0 μL,cDNA模板2.0 μL,ddH2O補(bǔ)足至20.0 μL。擴(kuò)增程序:94 ℃預(yù)變性2 min;94 ℃ 30 s,55 ℃ 30 s,72 ℃ 10 min,進(jìn)行35個循環(huán);72 ℃延伸10 min。PCR產(chǎn)物送往北京諾賽基因組研究中心有限公司測序。

    1. 5 生物信息學(xué)分析

    將測序結(jié)果提交至NCBI數(shù)據(jù)庫進(jìn)行BLAST比對,以獲得SlMAPK6蛋白的同源氨基酸序列。采用ProtParam對SlMAPK6蛋白進(jìn)行理化性質(zhì)分析;分別利用PSIPRED和SWISS-MOODEL預(yù)測SlMAPK6蛋白二、三級結(jié)構(gòu);采用Softberry-Protcomp進(jìn)行SlMAPK6蛋白亞細(xì)胞定位;采用TMHMM預(yù)測SlMAPK6蛋白跨膜區(qū);采用NetPhos 3.1預(yù)測SlMAPK6蛋白磷酸化位點(diǎn);通過ProtScale預(yù)測SlMAPK6蛋白親/疏水性;使用DNAMAN進(jìn)行氨基酸序列多重比對?;赟lMAPK6蛋白、番茄MAPK家族成員和擬南芥MAPK家族成員的氨基酸序列,用MEGA 7.0的鄰接法(Neighbor-joining,NJ)構(gòu)建系統(tǒng)發(fā)育進(jìn)化樹,Booststrap參數(shù)為1000。

    1. 6 實時熒光定量PCR檢測

    根據(jù)擴(kuò)增獲得的SlMAPK6基因序列設(shè)計定量引物qSlMAPK6-F和qSlMAPK6-R(表1)。以矮番茄葉片cDNA為模板,采用CFX96TM Real-time System熒光定量PCR儀檢測SlMAPK6基因在干旱、鹽脅迫及信號物質(zhì)處理下的相對表達(dá)量。反應(yīng)體系20.0 μL:SYBR Master Mix(2×) 10.0 μL,10 μmol/L正、反向引物各0.8 μL,cDNA模板4.0 μL,DEPC H2O補(bǔ)足至20.0 μL。每個樣品設(shè)3個重復(fù)。擴(kuò)增程序:95 ℃預(yù)變性10 min;95 ℃ 15 s,60 ℃ 15 s,進(jìn)行40個循環(huán);95 ℃ 10 s,65 ℃ 10 s,95 ℃ 5 s。以SlEF1α基因為內(nèi)參,其引物見表1。

    1. 7 統(tǒng)計分析

    使用Excel 2010整理處理數(shù)據(jù),采用2-△△Ct法計算目的基因的相對表達(dá)量(Livak and Schmittegen,2001),并以O(shè)rigin 2017制圖。

    2 結(jié)果與分析

    2. 1 SlMAPK6基因克隆結(jié)果

    以矮番茄葉片cDNA為模板進(jìn)行PCR擴(kuò)增,獲得1條長度約1300 bp的條帶(圖1),與預(yù)期結(jié)果相符。測序結(jié)果顯示,該片段全長1360 bp,開放閱讀框(ORF)為1131 bp,編碼376個氨基酸殘基。

    2. 2 SlMAPK6蛋白氨基酸序列同源性分析結(jié)果

    根據(jù)克隆獲得的SlMAPK6基因核苷酸序列推導(dǎo)其編碼蛋白(SlMAPK6)的氨基酸序列,通過NCBI數(shù)據(jù)庫的BLAST獲得其同源蛋白氨基酸序列,經(jīng)同源性分析發(fā)現(xiàn),SlMAPK6蛋白與馬鈴薯(Solanum tuberosum)StMAPK4、風(fēng)鈴辣椒(Capsicum baccatum)CbMAPK4和中華辣椒(C. chinense)CcMAPK4蛋白的同源性均為99.47%。對這4種同源蛋白進(jìn)行氨基酸序列多重比對分析,結(jié)果(圖2)顯示,SlMAPK6蛋白僅第149和311位氨基酸與StMAPK4蛋白存在差異,但與CbMAPK4和CcMAPK4蛋白無差異;CbMAPK4蛋白僅第43和226位氨基酸與CcMAPK4蛋白存在差異,但CcMAPK4蛋白第43位氨基酸與StMAPK4和SlMAPK6無差異;StMAPK4蛋白第138位氨基酸與其他3種蛋白均存在差異。綜上所述,SlMAPK6蛋白與馬鈴薯、風(fēng)鈴辣椒和中華辣椒MAPK4蛋白的同源性極高,序列高度保守。

    2. 3 SlMAPK6蛋白理化性質(zhì)預(yù)測結(jié)果

    SlMAPK6蛋白由376個氨基酸組成,分子量43.09 kD,總原子數(shù)6040,理論等電點(diǎn)(pI)6.76,帶負(fù)電荷氨基酸殘基總數(shù)[天冬氨酸(Asp)+谷氨酸(Glu)]47個,帶正電荷氨基酸殘基總數(shù)[精氨酸(Arg)+賴氨酸(Lys)]45個,不穩(wěn)定指數(shù)44.36,為不穩(wěn)定蛋白;脂肪系數(shù)為90.00,為強(qiáng)脂溶性蛋白。SlMAPK6蛋白親水性平均值(GRAVY)-0.396,多肽鏈第114位的Lys為最低分值(-2.578),親水性最強(qiáng);第234位亮氨酸(Leu)為最高分值(2.111),疏水性最強(qiáng),從整體來看,親水性氨基酸均勻分布在整個肽鏈中,故該蛋白為親水性蛋白(圖3)。Softberry-Protcomp亞細(xì)胞定位預(yù)測結(jié)果顯示,SlMAPK6蛋白可能存在于細(xì)胞核及胞漿質(zhì)中。TMHMM預(yù)測結(jié)果顯示,SlMAPK6蛋白無跨膜區(qū)(圖4)。

    MAPK作為末端級聯(lián)反應(yīng)信號激酶,磷酸化位點(diǎn)預(yù)測分析有助于了解其在細(xì)胞信號傳導(dǎo)過程中的作用。NetPhos 3.1預(yù)測結(jié)果顯示,SlMAPK6蛋白含有36個磷酸化位點(diǎn),其中,Ser磷酸位點(diǎn)19個,Thr磷酸位點(diǎn)9個(圖5)。

    2. 4 SlMAPK6蛋白二、三級結(jié)構(gòu)預(yù)測結(jié)果

    SlMAPK6蛋白二級結(jié)構(gòu)預(yù)測結(jié)果(圖6)顯示,α-螺旋占38.03%,延伸鏈占10.37%,無規(guī)則卷曲占51.60%,說明α-螺旋和無規(guī)則卷曲是SlMAPK6蛋白的主要結(jié)構(gòu)部件。SlMAPK6蛋白三級結(jié)構(gòu)預(yù)測結(jié)果(圖7)與二級結(jié)構(gòu)預(yù)測結(jié)果基本相符。圖6和圖7均顯示,SlMAPK6蛋白中無規(guī)則卷曲的占比最高。由于酶的功能部位常位于無規(guī)則卷曲中,可能是SlMAPK6在MAPK級聯(lián)反應(yīng)中發(fā)揮功能的主要原因。

    2. 5 系統(tǒng)發(fā)育進(jìn)化樹分析結(jié)果

    由圖8可知,SlMAPK6蛋白與SlMAPK5蛋白的親緣關(guān)系最近,與擬南芥MAPK家族成員AtMAPK4、AtMPK5、AtMPK11、AtMPK12和AtMPK13蛋白等親緣關(guān)系也較近。據(jù)前人報道,AtMAPK4蛋白在免疫調(diào)節(jié)中起負(fù)反饋調(diào)節(jié)(Ichimura et al.,2000),故推測SlMAPK6蛋白在番茄的免疫調(diào)節(jié)中發(fā)揮相似功能。

    2. 6 SlMAPK6基因在非生物脅迫及信號物質(zhì)處理下的表達(dá)模式

    對矮番茄幼苗進(jìn)行干旱(PEG-6000)、鹽(NaCl)及信號物質(zhì)(SA和BTH)處理,并采用實時熒光定量PCR檢測不同處理下SlMAPK6基因的表達(dá)模式。由圖9-A和圖9-B可知,SA處理后1和3 h,SlMAPK6基因表達(dá)量均較對照顯著升高(P<0.05,下同);BTH處理后3 h,SlMAPK6基因表達(dá)量較對照顯著升高,但SA和BTH處理后6~24 h,SlMAPK6基因的表達(dá)趨勢相似,均為抑制表達(dá)。由圖9-C可知,NaCl處理后,SlMAPK6基因表達(dá)較對照整體上無明顯變化,僅在處理后1和12 h與對照存在顯著差異,但處理后1 h為誘導(dǎo)表達(dá),處理后12 h為抑制表達(dá)。由圖9-D可知,PEG-6000處理后,SlMAPK6基因的整體表達(dá)趨勢為抑制表達(dá),僅在處理后6 h為誘導(dǎo)表達(dá)。綜上所述,經(jīng)SA、BTH、NaCl和PEG-6000處理后,SlMAPK6基因均出現(xiàn)抑制表達(dá),表明SlMAPK6基因負(fù)向調(diào)控番茄植株的逆境脅迫響應(yīng)。

    3 討論

    MAPK級聯(lián)信號參與調(diào)控植物的多種生命活動,如生長增殖、分化、運(yùn)動和凋亡及抗病、抗逆防御反應(yīng)(Abdukhi,2014)。MAPK處于MAPK級聯(lián)信號系統(tǒng)的下游,在整個信號傳導(dǎo)中發(fā)揮關(guān)鍵作用。番茄生長過程中常遭遇多種逆境脅迫,嚴(yán)重影響其生長發(fā)育而導(dǎo)致大幅減產(chǎn),因此研究MAPK級聯(lián)信號相關(guān)基因在番茄逆境中的調(diào)控功能尤為重要。目前,已從番茄中鑒定出16個MAPK家族成員,其中,A亞族的SlMAPK1、SlMAPK2和SlMAPK3基因研究較多,研究結(jié)果顯示,三者參與調(diào)控番茄植株的抗病、抗逆防御反應(yīng)(Li et al.,2017a,2017b;Wang et al.,2018b;楊楊等,2018)。經(jīng)研究證實,SlMAPK1和SlMAPK2蛋白與擬南芥AtMAPK6蛋白同源,而SlMAPK3蛋白與擬南芥AtMAPK3蛋白同源(孔福苓,2012)。AtMAPK3和AtMAPK6蛋白參與調(diào)控擬南芥多種生物與非生物脅迫(Beckers et al.,2009;Yu et al.,2010;Meng et al.,2013),即上述研究相互印證結(jié)論的準(zhǔn)確性。迄今,有關(guān)B亞族成員的研究報道極少,僅證實SlMAPK6屬于B亞族成員(Kong et al.,2012),而鮮見有關(guān)SlMAPK6基因克隆及表達(dá)的研究報道。本研究通過同源克隆技術(shù)獲得SlMAPK6基因序列,并對其進(jìn)行同源性分析,結(jié)果顯示,SlMAPK6蛋白與擬南芥AtMAPK4蛋白高度同源,推測SlMAPK6蛋白與擬南芥AtMAPK4蛋白功能相似。前人研究表明,AtMAPK4基因負(fù)調(diào)控擬南芥的自然免疫系統(tǒng)(Ichimura et al.,2000;Petersen et al.,2000;Kosetsu et al.,2011),故推測SlMAPK6基因負(fù)調(diào)控番茄植株的自然免疫系統(tǒng)。

    本研究通過生物信息學(xué)分析發(fā)現(xiàn),SlMAPK6是一個親水性蛋白,存在于細(xì)胞核及胞漿質(zhì)中,含有36個磷酸化位點(diǎn),其中Ser磷酸位點(diǎn)19個,Thr磷酸位點(diǎn)9個。經(jīng)前人研究證實,MAPK級聯(lián)反應(yīng)的磷酸化過程中絲氨酸和蘇氨酸發(fā)揮重要功能(Min and Zhao,2018),故推測SlMAPK6通過Ser磷酸位點(diǎn)和Thr磷酸位點(diǎn)激活下游蛋白發(fā)揮作用,后續(xù)可通過Western blotting技術(shù)進(jìn)行驗證。

    本研究對矮番茄幼苗分別進(jìn)行干旱(PEG-6000)、鹽(NaCl)及信號物質(zhì)(SA和BTH)處理,結(jié)果發(fā)現(xiàn)SlMAPK6蛋白參與番茄植株的多種非生物脅迫及信號物質(zhì)響應(yīng)。BTH是SA類似物,屬于人工合成的植物誘抗劑(于力等,2013),在SA和BTH的誘導(dǎo)下,SlMAPK6基因的表達(dá)趨勢整體相似,但又存在差異,暗示SlMAPK6基因參與SA和BTH誘導(dǎo)抗逆境響應(yīng)的作用機(jī)制存在差異。前人研究結(jié)果顯示,BTH能誘導(dǎo)植物抵抗多種逆境(Ramasamy et al.,2015;Frackowiak et al.,2019),而SA是植物體內(nèi)普遍存在的內(nèi)源信號分子,在植物的多種逆境信號中發(fā)揮重要作用(Rivas-San and Plasencia,2011;Oracz and Karpinski,2016;Li et al.,2017b;Wu et al.,2019),暗示SA和BTH誘導(dǎo)的SlMAPK6通路參與調(diào)控植物的抗逆響應(yīng)。經(jīng)NaCl處理后,SlMAPK6基因表達(dá)較對照整體上無明顯變化,僅在處理后1和12 h與對照存在顯著差異,但處理后1 h為誘導(dǎo)表達(dá),處理后12 h為抑制表達(dá),其原因可能是NaCl對植物的滲透脅迫需要時間,處理后12 h基因下調(diào)表達(dá)可能是植物對高濃度滲透脅迫的反饋調(diào)節(jié);同樣,PE-G6000處理6 h為誘導(dǎo)表達(dá),其他時間點(diǎn)均為抑制表達(dá),推測處理后6 h基因高效表達(dá)是植株對PEG-6000模擬干旱的逆境反饋調(diào)節(jié)。結(jié)合生物信息學(xué)分析結(jié)果,推測SlMAPK6基因可能受SA和BTH信號調(diào)節(jié)參與番茄的抗逆向應(yīng),且在番茄植株非生物逆境脅迫中發(fā)揮負(fù)調(diào)控作用。

    4 結(jié)論

    SlMAPK6基因可能參與負(fù)調(diào)控番茄植株的抗病和抗逆脅迫機(jī)制。SlMAPK6蛋白可能通過Ser或Thr磷酸化識別細(xì)胞信號的傳導(dǎo)功能,在矮番茄逆境脅迫響應(yīng)中發(fā)揮重要作用。

    參考文獻(xiàn):

    孔福苓. 2012. 番茄促分裂原活化蛋白激酶(SlMAPK)家族基因的分離及表達(dá)特征分析[D]. 杭州:浙江大學(xué). [Kong F L. 2012. Genome-wide identification and expression profiles analysis of mitogen-activated protein kinase family genes in Solanum lycopersicum[D]. Hangzhou:Zhejiang University.]

    李翠. 2014.番茄SpMPKs基因響應(yīng)非生物脅迫的功能分析[D]. 楊凌:西北農(nóng)林科技大學(xué). [Li C. 2014. The function of SpMPKs on tomato tolerance to abiotic stresses[D]. Yangling:Northwest A & F University.]

    王潔,王燕,潘長田,何艷軍,劉雪,盧鋼. 2016. 番茄SlMAPK9-2基因分離及表達(dá)分析[J]. 核農(nóng)學(xué)報,30(8):1480-1490. [Wang J,Wang Y,Pan C T,He Y J,Liu X,Lu G. 2016. Cloning and expression analysis of SlMAPK9-2 in Solanum lycopersicum L[J]. Journal of Nuclear Agricultural Sciences,30(8):1480-1490.]

    楊楊,劉燦,鄭鄢燕,生吉萍,申琳. 2018. SlMAPK1/2/3和SlACS2參與NO對抗病酶CHI和PAL的誘導(dǎo)[J]. 中國食品學(xué)報,18(7):18-25. [Yang Y,Liu C,Zheng Y R,Sheng J P,Shen L. 2018. SlMAPK1/2/3 and SlACS2 involved in NO-induced activity of phenylalnine ammonialyase and chitinase[J]. Journal of Chinese Institute of Food Science and Technology,18(7):18-25.]

    于力,孫錦,郭世榮,閻君,朱為民. 2013. 苯并噻二唑誘導(dǎo)番茄對番茄黃化曲葉病毒病的抗性[J]. 江蘇農(nóng)業(yè)學(xué)報, 29(1):71-75. [Yu L,Sun J,Guo S R,Yan J,Zhu W M. 2013. Resistance to tomato yellow leaf curl virus in tomato induced by benzothiadiazole[J]. Jiangsu Journal of Agricultural Sciences,29(1):71-75.]

    Abdukhi N. 2014. Functional analysis of MAP kinase pathways in plant defense responses in Solanum lycopersicum[D]. Ottawa:Carleton University.

    Beckers G J M,Jaskiewicz M,Liu Y D,Underwood W R,He S Y,Zhang S Q,Conrath U. 2009. Mitogen-activated protein kinases 3 and 6 are required for full priming of stress responses in Arabidopsis thaliana[J]. The Plant Cell,21(3):944-953.

    Chen Z,Gibson T B,Robinson F,Silvestro L,Pearson G, Xu B E,Wright A,Vanderbilt C,Cobb M H. 2001. MAP kinases[J]. Chemical Reviews,101(8):2449-2476.

    Crowe-White K M,Phillips T A,Ellis A C. 2019. Lycopene and cognitive function[J]. Journal of Nutritional Science.doi:10.1017/jns.2019.16.

    de Zelicourt A,Colcombet J,Hirt H. 2016. The role of MAPK modules and ABA during abiotic stress signaling[J]. Trends in Plant Science,21(8):677-685.

    Frackowiak P,Pospieszny H,Smiglak M,Obrepalska-Steplowska A. 2019. Assessment of the efficacy and mode of action of Benzo(1,2,3)-thiadiazole-7-carbothioic acid s-methyl ester(BTH) and its derivatives in plant protection against viral disease[J]. International Journal of Molecular Scien-ces,20(7):1-27.

    Hame L P,Nicole M C,Sritubtim S,Morency M J,Ellis M,Ehlting J,Beaudoin N,Barbazuk B,Klessig D,Lee J,Martin G,Mundy J,Ohashi Y,Scheel D,Sheen J,Xing T,Zhang S Q,Seguin A,Brian E,Ellis B E. 2006. Ancient signals:Comparative genomics of plant MAPK and MAPKK gene families[J]. Trends in Plant Science,11(4):1360-1385.

    Ichimura K,Mizoguchi T,Yoshida R,Yuasa T,Shinozaki K. 2000. Various abiotic stresses rapidly activate Arabidopsis MAP kinases ATMPK4 and ATMPK6[J]. The Plant Journal,24(5):655-665.

    Kalapos B,Hlavová M,Nádai T V,Galiba G, Bi?ová K,Dóczi R. 2019. Early evolution of the mitogen-activated protein kinase family in the plant kingdom[J]. Scientific Reports,9:1-14.

    Kong F L,Wang J,Cheng L,Liu S Y,Wu J,Peng Z,Lu G. 2012. Genome-wide analysis of the mitogen-activated protein kinase gene family in Solanum lycopersicum[J]. Gene,499(1):108-120.

    Kosetsu K,Matsunaga S,Nakagami H,Colcombet J,Sasabe M,Soyano T,Takahashi Y,HirtH,Machida Y. 2011. The MAP kinase MPK4 is required for cytokinesis in Arabidopsis thaliana[J]. The Plant Cell,22(11):3778-3790.

    Lee Y, Kim Y J,Kim M H,Kwak J M. 2016. MAPK cascades in guard cell signal transduction[J]. Frontiers in Plant Science. doi:10.3389/fpls.2016.00080.

    Li Y Z,Qin L,Zhao J J,Muhammad T,Cao H H,Li H L,Zhang Y,Liang Y. 2017a. SlMAPK3 enhances tolerance to tomato yellow leaf curl virus(TYLCV) by regulating salicylic acid and jasmonic acid signaling in tomato(Solanum lycopersicum)[J]. PLoS One,12(2):e0172466.

    Li Z,Yu J J,Peng Y,Huang B R. 2017b. Metabolic pathways regulated by abscisic acid,salicylic acid and γ-aminobutyric acid in association with improved drought tolerance in creeping bentgrass(Agrostis stolonifera)[J]. Physiologia Plantarum,159(1):42-58.

    Lian W W,Tang Y M,Gao S Q,Zhang Z,Zhao X,Zhao C P. 2012. Phylogenetic analysis and expression patterns of the MAPK gene family in wheat(Triticum aestivum L.)[J]. Agricultural Sciences in China,11(8):1227-1235.

    Liu Y K,Zhang D,Wang L,Li D Q. 2013. Genome-wide analysis of mitogen-activated protein kinase gene family in maize[J]. Plant Molecular Biology Reporter,31(6):1446-1460.

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

    Meng X Z,Xu J,He Y X,Yang K Y,Mordorski B,Liu Y D,Zhang S Q. 2013. Phosphorylation of an ERF transcription factor by Arabidopsis MPK3/MPK6 regulates plant defense gene induction and fungal resistance[J]. The Plant Cell,25(3):1126-1142.

    Min J, Zhao Q C. 2018. Comparative analysis of plant MKK gene family reveals novel expansion mechanism of the members and sheds new light on functional conservation[J]. BMC Genomics,19(1):407-424.

    Mine A,Berens M L,Nobori T,Anver S,F(xiàn)ukumoto K,Winkelmuller T M,Takeda A,Becker D,Tsuda K. 2017. Pathogen exploitation of an abscisic acid- and jasmonate-inducible MAPK phosphatase and its interception by Arabidopsis immunity[J]. Proceedings of the National Aca-demy of Sciences of the United States of America,114(28):7456-7461.

    Mohanta T K,Arora P K,Mohanta N,Parida P,Bae H. 2015. Identification of new members of the MAPK gene family in plants shows diverse conserved domains and novel activation loop variants[J]. BMC Genomics,16(1):58-77.

    Oracz K,Karpinski S. 2016. Phytohormones signaling pathways and ROS involvement in seed germination[J]. Frontiers in Plant Science,7:864. doi:10.3389/fpls.2016.00864.

    Ortiz-Masia D,Perez-Amador M A,Carbonell J,Marcote M J. 2007. Diverse stress signals activate the C1 subgroup MAP kinases of Arabidopsis[J]. FEBS Letters,581(9):1834-1840.

    Perveen R,Suleria H A R,Anjum F M,Butt M S,Pasha I,Ahmad S. 2015. Tomato(Solanum lycopersicum) carote-noids and lycopenes chemistry; metabolism,absorption,nutrition,and allied health claims—A comprehensive review[J]. Critical Reviews in Food Sicence and Nutrition,55(7):919-929.

    Petersen M,Brodersen P,Naested H,Andreasson E,Lindhart U,Johansen B,Nielsen H B,Lacy M,Austin M J,Parker J E,Sharma S B,Klessig D F,Martienssen R,Mattsson O,Jensen A B,Mundy J. 2000. Arabidopsis MAP kinase 4 negatively regulates systemic acquired resistance[J]. Cell,103(7):1111-1120.

    Petyaev I M. 2016. Lycopene deficiency in ageing and cardiovascular disease[J]. Oxidative Medicine and Cellular Longevity. doi:10.1155/2016/3218605.

    Ramasamy A D,Bokshi A I,Phan-Thien K,Mcconchie R M. 2015. Seed treatment with benzothiadiazole induces resistance against powdery mildew disease caused by sphaerotheca fuliginea and increases the activities of pathogenesis-related enzymes in cucumber plants[J]. Journal of Horticultural Science & Biotechnology,90(1):63-70.

    Rivas-San V M,Plasencia J. 2011. Salicylic acid beyond defence:Its role in plant growth and development[J]. Journal of Experimental Botany,62(10):3321-3338.

    Rodriguez M C S,Petersen M,Mundy J. 2010. Mitogen-activated protein kinase signaling in plants[J]. Annual Review of Plant Biology,61:621-649. doi:10.1146/annure-varplant-042809-112252.

    Singh A,Nath O,Singh S,Kumar S,Singh I K. 2018. Genome-wide identification of the MAPK gene family in chickpea and expression analysis during development and stress response[J]. Plant Gene,13:25-35. doi:10.1016/j.plgene.2017.12.001.

    Smekalova V,Doskocilova A,Komis G,Samaj J. 2013. Crosstalk between secondary messengers,hormones and MAPK modules during abiotic stress signalling in plants[J]. Biotechnology Advances,32(1):2-11.

    Wang H,Gong M,Guo J Y,Xin H,Gao Y,Chao Liu C,Dai D Q,Tang L Z. 2018. Genome-wide identification of Jatropha curcas MAPK,MAPKK,and MAPKKK gene fa-milies and their expression profile under cold stress[J]. Scientific Reports,8(1):16163. doi:10.1038/s41598-018-34614-1.

    Wang L,Zhao R R,Li R,Yu W Q,Yang M J,Sheng J P,Shen L. 2018. Enhanced drought tolerance in tomato plants by overexpression of SlMAPK1[J]. Plant Cell Ti-ssue and Organ Culture,133(1):27-38.

    Wu Z J,Han S M,Zhou H D,Tuang Z H,Wang Y Z, Jin Y, Shi H Z, Yang W N. 2019. Cold stress activates disease resistance in Arabidopsis thaliana through a salicylic acid dependent pathway[J]. Plant Cell and Environment,42(9):2645-2663.

    Xu J,Zhang S Q. 2015. Mitogen-activated protein kinase cascades in signaling plant growth and development[J]. Trends in Plant Science,20(1):56-64.

    Yu L J,Nie J N,Cao C Y,Jin Y K,Yan M,Wang F Z,Liu J,Xiao Y,Liang Y H,Zhang W H. 2010. Phosphatidic acid mediates salt stress response by regulation of MPK6 in Arabidopsis thaliana[J]. New Phytologist,188(3):762-773.

    Zhang S J,Wang L,Zhao R R,Yu W Q,Li R,Li Y J,Sheng J P,Shen L. 2018. Knockout of SlMAPK3 reduced di-sease resistance to botrytis cinerea in tomato plants[J]. Journal of Agricultural and Food Chemistry,66(34):8949-8956.

    (責(zé)任編輯 陳 燕)

    猜你喜歡
    非生物脅迫基因克隆番茄
    番茄炒蛋
    秋茬番茄“疑難雜癥”如何挽救
    番茄果實“起棱”怎么辦
    白木香轉(zhuǎn)錄因子AsMYB1和AsMYB2克隆及表達(dá)分析
    三個小麥防御素基因的克隆及序列分析
    非生物脅迫對擬南芥IQM4基因表達(dá)的影響
    科技視界(2016年16期)2016-06-29 11:55:38
    水稻富亮氨酸類受體蛋白激酶參與外界脅迫應(yīng)答的研究進(jìn)展
    山桃醇腈酶基因PdHNL1的克隆、序列分析與原核表達(dá)
    紅花低分子量油體蛋白基因的克隆與系統(tǒng)進(jìn)化分析
    葡萄查耳酮合酶基因克隆及其進(jìn)化分析
    亚洲五月天丁香| 国产片内射在线| 久久精品国产综合久久久| 成人特级黄色片久久久久久久| 国产av在哪里看| www.999成人在线观看| av超薄肉色丝袜交足视频| 首页视频小说图片口味搜索| 国产亚洲精品一区二区www| 国产高清视频在线播放一区| 国产午夜精品久久久久久| 亚洲五月婷婷丁香| 久久人妻av系列| 久久久久国产精品人妻aⅴ院| 欧美日韩亚洲高清精品| 日本欧美视频一区| 黑人巨大精品欧美一区二区蜜桃| 在线观看午夜福利视频| 日日干狠狠操夜夜爽| 精品久久久久久久毛片微露脸| 天堂√8在线中文| 神马国产精品三级电影在线观看 | 亚洲欧美一区二区三区久久| 国产三级在线视频| 嫩草影院精品99| 国产欧美日韩综合在线一区二区| 热re99久久国产66热| 999精品在线视频| 国产91精品成人一区二区三区| 啪啪无遮挡十八禁网站| 丰满人妻熟妇乱又伦精品不卡| 女性被躁到高潮视频| 亚洲美女黄片视频| 51午夜福利影视在线观看| 免费不卡黄色视频| 99国产精品一区二区三区| 女人高潮潮喷娇喘18禁视频| 国产亚洲精品久久久久5区| 久久久久久亚洲精品国产蜜桃av| 99国产精品一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 高清黄色对白视频在线免费看| 三级毛片av免费| 成人国语在线视频| 多毛熟女@视频| 最近最新中文字幕大全电影3 | 精品久久久精品久久久| 精品少妇一区二区三区视频日本电影| 亚洲久久久国产精品| 高清毛片免费观看视频网站 | 亚洲第一av免费看| 一级片免费观看大全| 熟女少妇亚洲综合色aaa.| 老司机亚洲免费影院| 亚洲国产精品999在线| 叶爱在线成人免费视频播放| 自拍欧美九色日韩亚洲蝌蚪91| 国产免费男女视频| 久热这里只有精品99| 交换朋友夫妻互换小说| 天堂影院成人在线观看| 久久婷婷成人综合色麻豆| 99精品欧美一区二区三区四区| 久久久久国内视频| 午夜免费鲁丝| 久久人人97超碰香蕉20202| 成人免费观看视频高清| 成年女人毛片免费观看观看9| 久久精品亚洲熟妇少妇任你| av视频免费观看在线观看| 免费观看人在逋| 麻豆av在线久日| 久久九九热精品免费| 超色免费av| 操出白浆在线播放| www国产在线视频色| 多毛熟女@视频| 狂野欧美激情性xxxx| 国产国语露脸激情在线看| 亚洲成国产人片在线观看| 首页视频小说图片口味搜索| 欧美日韩亚洲国产一区二区在线观看| 国产单亲对白刺激| 日本三级黄在线观看| 熟女少妇亚洲综合色aaa.| 自拍欧美九色日韩亚洲蝌蚪91| 成人黄色视频免费在线看| 两人在一起打扑克的视频| 免费在线观看日本一区| 九色亚洲精品在线播放| 亚洲午夜理论影院| 水蜜桃什么品种好| 成人特级黄色片久久久久久久| 无遮挡黄片免费观看| 亚洲一区二区三区色噜噜 | 久久精品91蜜桃| 女人被躁到高潮嗷嗷叫费观| 午夜两性在线视频| 亚洲va日本ⅴa欧美va伊人久久| 一本综合久久免费| 精品一区二区三卡| 日韩大码丰满熟妇| 免费观看精品视频网站| 无遮挡黄片免费观看| 久久久久久久精品吃奶| 亚洲在线自拍视频| av中文乱码字幕在线| 精品免费久久久久久久清纯| 9色porny在线观看| 大型黄色视频在线免费观看| 满18在线观看网站| 午夜两性在线视频| 欧美黑人欧美精品刺激| aaaaa片日本免费| 精品久久久精品久久久| 男女做爰动态图高潮gif福利片 | 欧美色视频一区免费| 国产99久久九九免费精品| 黑人巨大精品欧美一区二区蜜桃| 老鸭窝网址在线观看| 亚洲精品国产色婷婷电影| 亚洲精品国产区一区二| 国产一区二区在线av高清观看| 成年人黄色毛片网站| 在线观看免费高清a一片| 日本五十路高清| 性色av乱码一区二区三区2| 成人亚洲精品av一区二区 | 热99re8久久精品国产| 午夜91福利影院| 在线观看一区二区三区激情| 亚洲五月婷婷丁香| 国产成+人综合+亚洲专区| 国产一区在线观看成人免费| 成年人黄色毛片网站| 亚洲中文日韩欧美视频| 亚洲精品久久成人aⅴ小说| 在线观看午夜福利视频| 欧美久久黑人一区二区| a级毛片黄视频| 国产一区二区三区视频了| 丁香欧美五月| 一二三四社区在线视频社区8| 欧美日韩中文字幕国产精品一区二区三区 | 久久天堂一区二区三区四区| 免费搜索国产男女视频| 一区二区三区国产精品乱码| 国产高清激情床上av| 一区二区三区精品91| 国产99久久九九免费精品| 国产亚洲精品第一综合不卡| 1024视频免费在线观看| 在线国产一区二区在线| 久久天堂一区二区三区四区| 99热国产这里只有精品6| 国产伦一二天堂av在线观看| av欧美777| 久久精品亚洲av国产电影网| 亚洲一区二区三区不卡视频| 国产一卡二卡三卡精品| e午夜精品久久久久久久| a级毛片黄视频| 99精品欧美一区二区三区四区| 免费在线观看亚洲国产| 亚洲 欧美 日韩 在线 免费| 久久精品国产亚洲av高清一级| 久久影院123| 三上悠亚av全集在线观看| 国产精品av久久久久免费| 欧美另类亚洲清纯唯美| 亚洲成人久久性| 久久人妻熟女aⅴ| 亚洲国产中文字幕在线视频| 伊人久久大香线蕉亚洲五| 午夜亚洲福利在线播放| 可以在线观看毛片的网站| 亚洲成人免费av在线播放| 免费看十八禁软件| 无限看片的www在线观看| 久久精品国产清高在天天线| 国产97色在线日韩免费| 国产一区二区三区视频了| 国产有黄有色有爽视频| 成人影院久久| 99国产极品粉嫩在线观看| 母亲3免费完整高清在线观看| 欧美最黄视频在线播放免费 | 操出白浆在线播放| 国产视频一区二区在线看| 一个人免费在线观看的高清视频| 久久午夜综合久久蜜桃| 亚洲国产精品sss在线观看 | 一边摸一边抽搐一进一出视频| 精品乱码久久久久久99久播| 精品欧美一区二区三区在线| 午夜精品国产一区二区电影| 一个人观看的视频www高清免费观看 | av网站在线播放免费| 激情视频va一区二区三区| 免费高清在线观看日韩| 久久久久久久久免费视频了| 亚洲少妇的诱惑av| 如日韩欧美国产精品一区二区三区| 日韩免费高清中文字幕av| 黄片小视频在线播放| 国产区一区二久久| 国产精品综合久久久久久久免费 | 神马国产精品三级电影在线观看 | 国产av在哪里看| 国产91精品成人一区二区三区| 免费高清视频大片| 高清黄色对白视频在线免费看| 日韩一卡2卡3卡4卡2021年| 国产高清国产精品国产三级| 色播在线永久视频| 欧美激情久久久久久爽电影 | 亚洲专区国产一区二区| 午夜影院日韩av| 男女之事视频高清在线观看| av片东京热男人的天堂| 欧美人与性动交α欧美软件| 波多野结衣高清无吗| 在线天堂中文资源库| 交换朋友夫妻互换小说| bbb黄色大片| 黑人操中国人逼视频| 精品国产超薄肉色丝袜足j| 久久精品国产亚洲av香蕉五月| 成在线人永久免费视频| 日本一区二区免费在线视频| 成人三级做爰电影| 国产1区2区3区精品| 精品一区二区三区视频在线观看免费 | 在线观看免费高清a一片| 亚洲欧美激情在线| 久久久久久大精品| 欧美乱码精品一区二区三区| 亚洲精品在线观看二区| 在线观看日韩欧美| 亚洲欧美精品综合一区二区三区| 久久精品国产亚洲av高清一级| av国产精品久久久久影院| 日本 av在线| 久久久精品国产亚洲av高清涩受| 亚洲精品在线美女| www.www免费av| 看片在线看免费视频| 在线永久观看黄色视频| 欧美另类亚洲清纯唯美| 成人三级黄色视频| 日韩有码中文字幕| 亚洲熟女毛片儿| 中文字幕人妻丝袜制服| 欧美最黄视频在线播放免费 | x7x7x7水蜜桃| 老司机靠b影院| 人成视频在线观看免费观看| 日韩国内少妇激情av| 久久香蕉国产精品| 国产精品一区二区在线不卡| 国产av一区二区精品久久| 国产区一区二久久| 波多野结衣高清无吗| 国产色视频综合| 曰老女人黄片| 精品免费久久久久久久清纯| 美国免费a级毛片| 国产精品二区激情视频| 国产日韩一区二区三区精品不卡| 午夜福利一区二区在线看| 亚洲一区高清亚洲精品| 欧美+亚洲+日韩+国产| 大香蕉久久成人网| 一二三四在线观看免费中文在| 中文字幕人妻丝袜制服| 一级毛片高清免费大全| 男女床上黄色一级片免费看| 亚洲欧洲精品一区二区精品久久久| 国产野战对白在线观看| 99香蕉大伊视频| 久久国产精品影院| 啪啪无遮挡十八禁网站| 一边摸一边抽搐一进一出视频| 亚洲中文av在线| 老司机福利观看| 国产免费男女视频| 亚洲成人免费av在线播放| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品国产色婷婷电影| 亚洲av片天天在线观看| 日本黄色视频三级网站网址| 精品乱码久久久久久99久播| 久久天躁狠狠躁夜夜2o2o| 免费av中文字幕在线| 免费日韩欧美在线观看| 女人高潮潮喷娇喘18禁视频| 久久人妻福利社区极品人妻图片| 午夜免费激情av| av天堂久久9| 一进一出抽搐gif免费好疼 | 欧美乱码精品一区二区三区| 80岁老熟妇乱子伦牲交| 欧美人与性动交α欧美精品济南到| 久久人妻av系列| 亚洲一区二区三区色噜噜 | 国产成人啪精品午夜网站| 中文亚洲av片在线观看爽| 国产av精品麻豆| 老司机午夜十八禁免费视频| 国产乱人伦免费视频| 另类亚洲欧美激情| 精品无人区乱码1区二区| 啦啦啦在线免费观看视频4| 婷婷精品国产亚洲av在线| 亚洲一码二码三码区别大吗| 我的亚洲天堂| 每晚都被弄得嗷嗷叫到高潮| 深夜精品福利| e午夜精品久久久久久久| 亚洲三区欧美一区| 久久中文看片网| 波多野结衣一区麻豆| 日日夜夜操网爽| 婷婷精品国产亚洲av在线| 午夜福利在线免费观看网站| 黑人猛操日本美女一级片| 国产亚洲精品久久久久久毛片| 午夜福利免费观看在线| 日本wwww免费看| 一本综合久久免费| 人人妻人人爽人人添夜夜欢视频| 亚洲一区高清亚洲精品| 亚洲三区欧美一区| 国产成人精品久久二区二区91| 免费在线观看日本一区| 在线观看免费视频网站a站| 国产精品偷伦视频观看了| a级毛片黄视频| 老汉色av国产亚洲站长工具| www.精华液| 国产精品影院久久| 最近最新免费中文字幕在线| 国产熟女xx| 亚洲精品一二三| 美女 人体艺术 gogo| 国产精品1区2区在线观看.| 久久久精品欧美日韩精品| 国产日韩一区二区三区精品不卡| 国产精品香港三级国产av潘金莲| 久久香蕉激情| 久久国产精品影院| 天天添夜夜摸| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品sss在线观看 | 日本wwww免费看| www日本在线高清视频| 国产av一区二区精品久久| 免费在线观看黄色视频的| 日韩精品中文字幕看吧| 亚洲成av片中文字幕在线观看| 十分钟在线观看高清视频www| 999久久久精品免费观看国产| 香蕉丝袜av| 久久国产精品人妻蜜桃| 免费看a级黄色片| 十八禁人妻一区二区| 精品高清国产在线一区| 亚洲精品久久成人aⅴ小说| 99热国产这里只有精品6| 曰老女人黄片| 色综合站精品国产| 女人精品久久久久毛片| 身体一侧抽搐| 高清黄色对白视频在线免费看| 日日摸夜夜添夜夜添小说| 亚洲精品在线观看二区| 中亚洲国语对白在线视频| 国产精品一区二区免费欧美| 成人av一区二区三区在线看| 中文字幕人妻丝袜制服| 亚洲在线自拍视频| 欧美av亚洲av综合av国产av| 日韩高清综合在线| 国产熟女xx| 日本黄色日本黄色录像| 女人被狂操c到高潮| 久久久水蜜桃国产精品网| 精品久久久久久成人av| 国产亚洲av高清不卡| 黑丝袜美女国产一区| 正在播放国产对白刺激| 欧美日本亚洲视频在线播放| 男人的好看免费观看在线视频 | 欧美av亚洲av综合av国产av| 久9热在线精品视频| 欧美人与性动交α欧美软件| 亚洲精品粉嫩美女一区| 熟女少妇亚洲综合色aaa.| 神马国产精品三级电影在线观看 | 亚洲一卡2卡3卡4卡5卡精品中文| 国产成+人综合+亚洲专区| 在线天堂中文资源库| 成熟少妇高潮喷水视频| 18禁观看日本| 大陆偷拍与自拍| 最新美女视频免费是黄的| 亚洲国产精品999在线| 真人一进一出gif抽搐免费| 高清欧美精品videossex| 日韩欧美三级三区| 女人被狂操c到高潮| 久久婷婷成人综合色麻豆| 男女下面插进去视频免费观看| 精品国产乱子伦一区二区三区| 精品国产一区二区三区四区第35| 99精品久久久久人妻精品| av免费在线观看网站| 国产免费现黄频在线看| 成人影院久久| 亚洲自拍偷在线| 久久久国产成人免费| 老司机在亚洲福利影院| 国产精品爽爽va在线观看网站 | 国产成人系列免费观看| 一区在线观看完整版| 精品福利观看| 极品教师在线免费播放| 夜夜躁狠狠躁天天躁| 久久国产精品男人的天堂亚洲| 国产欧美日韩综合在线一区二区| 一边摸一边做爽爽视频免费| 黄频高清免费视频| 窝窝影院91人妻| 亚洲激情在线av| 亚洲精品国产色婷婷电影| 中文字幕人妻丝袜一区二区| 99久久99久久久精品蜜桃| 亚洲av第一区精品v没综合| 久久精品亚洲熟妇少妇任你| 黑丝袜美女国产一区| 亚洲午夜精品一区,二区,三区| 美女高潮喷水抽搐中文字幕| 老司机午夜福利在线观看视频| 亚洲欧美激情在线| 亚洲av电影在线进入| 午夜久久久在线观看| 最好的美女福利视频网| 一级作爱视频免费观看| 一二三四在线观看免费中文在| av在线天堂中文字幕 | 亚洲欧美一区二区三区久久| 日本黄色日本黄色录像| xxx96com| 亚洲激情在线av| 亚洲男人的天堂狠狠| 精品久久久久久,| aaaaa片日本免费| 在线观看66精品国产| 亚洲性夜色夜夜综合| 午夜免费鲁丝| 丁香欧美五月| 欧美成人免费av一区二区三区| 精品久久久久久,| 十分钟在线观看高清视频www| 女人被躁到高潮嗷嗷叫费观| 男人的好看免费观看在线视频 | 99国产极品粉嫩在线观看| 午夜免费激情av| 亚洲自偷自拍图片 自拍| 一进一出好大好爽视频| www国产在线视频色| 好男人电影高清在线观看| 欧美色视频一区免费| 一边摸一边做爽爽视频免费| 日日干狠狠操夜夜爽| 午夜福利影视在线免费观看| 久久久久久久午夜电影 | 他把我摸到了高潮在线观看| 校园春色视频在线观看| 国产av在哪里看| 精品一区二区三区四区五区乱码| 91字幕亚洲| 国产欧美日韩一区二区精品| 身体一侧抽搐| 制服诱惑二区| 色哟哟哟哟哟哟| 后天国语完整版免费观看| 超碰97精品在线观看| 日韩高清综合在线| 亚洲五月色婷婷综合| 色精品久久人妻99蜜桃| 777久久人妻少妇嫩草av网站| 亚洲aⅴ乱码一区二区在线播放 | 欧美日韩av久久| 亚洲七黄色美女视频| 亚洲av五月六月丁香网| 丝袜人妻中文字幕| 黄网站色视频无遮挡免费观看| 国产精品久久久久久人妻精品电影| 欧美日韩精品网址| 91av网站免费观看| 亚洲欧美一区二区三区久久| 嫩草影院精品99| 亚洲一区二区三区不卡视频| 国产精品日韩av在线免费观看 | 午夜日韩欧美国产| 视频区欧美日本亚洲| 好看av亚洲va欧美ⅴa在| 欧美日韩乱码在线| 日韩国内少妇激情av| 国产成人精品久久二区二区91| 老司机亚洲免费影院| 国产亚洲精品久久久久5区| 日韩免费高清中文字幕av| 日本免费a在线| 99国产精品一区二区三区| 国产成人精品在线电影| 久久国产亚洲av麻豆专区| 亚洲性夜色夜夜综合| 丝袜美足系列| 欧美日本中文国产一区发布| 色尼玛亚洲综合影院| 国产欧美日韩一区二区三| 亚洲欧美精品综合久久99| 精品久久久久久,| 99香蕉大伊视频| 国产成人欧美| 欧美日韩一级在线毛片| 欧美+亚洲+日韩+国产| 老司机午夜十八禁免费视频| 日本免费一区二区三区高清不卡 | 夜夜躁狠狠躁天天躁| 精品一品国产午夜福利视频| 成熟少妇高潮喷水视频| 日韩欧美一区视频在线观看| 久久人妻av系列| 19禁男女啪啪无遮挡网站| 久久 成人 亚洲| 90打野战视频偷拍视频| 国产精品久久久久成人av| 免费女性裸体啪啪无遮挡网站| 一区二区三区国产精品乱码| 日本精品一区二区三区蜜桃| 我的亚洲天堂| 亚洲av电影在线进入| 在线十欧美十亚洲十日本专区| 国产成人精品久久二区二区91| 欧美日韩精品网址| 成人手机av| 每晚都被弄得嗷嗷叫到高潮| 国产精品秋霞免费鲁丝片| 亚洲成人免费av在线播放| 成年女人毛片免费观看观看9| 天天影视国产精品| 多毛熟女@视频| 99热只有精品国产| 欧美最黄视频在线播放免费 | 亚洲成国产人片在线观看| 女人高潮潮喷娇喘18禁视频| 老司机福利观看| 久久精品国产清高在天天线| 精品国产乱子伦一区二区三区| 久久久精品欧美日韩精品| 99国产精品一区二区三区| 日韩大码丰满熟妇| 国产成人欧美| 女人精品久久久久毛片| 国产高清激情床上av| 一级毛片女人18水好多| 日韩视频一区二区在线观看| 亚洲av五月六月丁香网| 青草久久国产| 亚洲精品一区av在线观看| svipshipincom国产片| 国产精品九九99| www日本在线高清视频| 在线观看午夜福利视频| 一区二区三区激情视频| 中文字幕人妻熟女乱码| 日韩成人在线观看一区二区三区| 亚洲avbb在线观看| av视频免费观看在线观看| 欧美最黄视频在线播放免费 | 高清av免费在线| a级片在线免费高清观看视频| 看免费av毛片| 淫秽高清视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 男女床上黄色一级片免费看| 亚洲中文日韩欧美视频| 久久精品国产亚洲av高清一级| 一级a爱视频在线免费观看| 男人操女人黄网站| 亚洲精品在线观看二区| 国内久久婷婷六月综合欲色啪| 国产欧美日韩精品亚洲av| 美女福利国产在线| 夜夜爽天天搞| 国产精品综合久久久久久久免费 | 欧美日韩福利视频一区二区| 亚洲第一av免费看| 91大片在线观看| 亚洲人成77777在线视频| 女性被躁到高潮视频| 久久性视频一级片| 伊人久久大香线蕉亚洲五| 欧美丝袜亚洲另类 | 三级毛片av免费| 男女做爰动态图高潮gif福利片 | 法律面前人人平等表现在哪些方面| 免费av毛片视频| 精品国产乱子伦一区二区三区| 国产亚洲欧美在线一区二区| 亚洲中文日韩欧美视频| 少妇裸体淫交视频免费看高清 |