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

    石榴MAPK家族基因鑒定及其響應(yīng)冷脅迫的表達(dá)分析

    2023-10-27 09:00:25陳利娜曹尚銀唐麗穎李好先嚴(yán)瓊李松開楊慶華魯振華
    果樹學(xué)報 2023年10期
    關(guān)鍵詞:表達(dá)分析抗寒性石榴

    陳利娜 曹尚銀 唐麗穎 李好先 嚴(yán)瓊 李松開 楊慶華 魯振華

    DOI:10.13925/j.cnki.gsxb.20230136

    摘? ? 要:【目的】評價不同石榴種質(zhì)資源籽粒硬度及抗寒性,篩選可能參與調(diào)控石榴抗寒性的MAPK家族基因?!痉椒ā恳?1份石榴種質(zhì)資源為試材,進行抗寒性及籽粒硬度評價;基于全基因組篩選石榴MAPK家族基因,對其進行進化、基因結(jié)構(gòu)和蛋白理化性質(zhì)分析,同時利用實時熒光定量PCR(real-time quantitative PCR,qRT-PCR)分析冷脅迫對石榴MAPK家族基因表達(dá)模式的影響。【結(jié)果】31個石榴品種籽粒硬度及半致死溫度測定結(jié)果表明,嶧城粉紅牡丹、淮北六棱甜和魯白榴2號等硬籽石榴抗寒性較強,華紫、以3和瑪麗斯等軟籽石榴抗寒性較弱。基于石榴全基因組鑒定出17個MAPK家族基因,廣泛分布于8條染色體上;MAPK家族所有成員主要分為3個亞類,其中,A和B亞類成員主要包含PKc_MAPKK_plant_like和PTZ00024結(jié)構(gòu)域,C亞類主要包含PLN00034結(jié)構(gòu)域,所有成員均具有S_TKc結(jié)構(gòu)域;各成員氨基酸殘基數(shù)量分布在314~860 aa,外顯子數(shù)目1~18個,蛋白分子質(zhì)量為34 910.05~97 965.26 u,等電點4.94~9.35;PgMKK2、PgMPK6、PgMPK9、PgMPK16和PgMPK13在嶧城粉紅牡丹響應(yīng)冷脅迫過程中表現(xiàn)為顯著上調(diào),PgMKK8、PgMPK1-1和PgMKK4在瑪麗斯響應(yīng)冷脅迫過程中表現(xiàn)為顯著上調(diào);PgMKK2、PgMPK6、PgMPK9、PgMPK16和PgMPK13在嶧城粉紅牡丹響應(yīng)冷脅迫過程中的表達(dá)量顯著高于瑪麗斯,PgMKK8和PgMPK1-1在瑪麗斯響應(yīng)冷脅迫過程中的表達(dá)量顯著性高于嶧城粉紅牡丹;PgMKK3在嶧城粉紅牡丹不同時間均未檢測到表達(dá),在瑪麗斯中表現(xiàn)為先升高后降低的趨勢;PgMPK12-2在瑪麗斯不同時間均未檢測到表達(dá),在嶧城粉紅牡丹中表現(xiàn)為逐漸升高的趨勢?!窘Y(jié)論】石榴MAPK家族基因響應(yīng)冷脅迫信號,其中,PgMKK2、PgMPK6、PgMPK12-2和PgMPK9可能參與正調(diào)控石榴的抗寒性。

    關(guān)鍵詞:石榴;籽粒硬度;抗寒性;MAPK;表達(dá)分析

    中圖分類號:S665.4 文獻標(biāo)志碼:A 文章編號:1009-9980(2023)10-2076-13

    Identification of MAPK family genes and analysis of their expression patterns in response to cold stress in pomegranate

    CHEN Lina1, CAO Shangyin1#, TANG Liying2, LI Haoxian1, YAN Qiong3, LI Songkai3, YANG Qinghua3, LU Zhenhua1*

    (1Zhengzhou Fruit Research Institute, Chinese Academy of Agricultural Sciences, Zhengzhou 450009, Henan, China; 2Hebei Chengde Academy of Agriculture and Forestry Sciences, Chengde 067000, Hebei, China; 3Mengzi Fruit and Vegetable Technology Promotion Station, Mengzi 661100, Yunnan, China)

    Abstract: 【Objective】 Cold stress is one of the most important factors limiting the progress of pomegranate production. The objectives of this study were to compare the cold hardiness among 31 pomegranate cultivars with wide distribution of seed hardness, and analyze the expression pattern of Mitogen-activated protein kinase (MAPK) family genes in response to cold stress. The results can pave the way for shedding light to the function of Mitogen-activated protein kinase cascades in pomegranate cod tolerance. 【Methods】Plant materials were cultivated in the National Horticulture Germplasm Resources Center of China (NHGRC), with conventional cultivation management. Juvenile branches of 31 pomegranate cultivars were sampled in mid-January, the relative electrolyte conductivity (REC) was measured, and the semi-lethal temperature (LT50) was calculated for the evolution of cold hardiness. Seed hardness was detected by the TA-XT texture apparatus. MAPK family genes of Arabidopsis were used as queries to search in the whole pomegranate genome database, and reference genome sequence of ‘Tunisia was obtained from NCBI database. Potential members of MAPK family were identified. Phylogeny relationship, gene structure and protein physicochemical properties were analyzed. Moreover, we performed real-time quantitative PCR (qRT-PCR) to analyze the expression pattern of 17 MAPK family genes in response to cold stress. 【Results】 The seed hardness and cold tolerance were evaluated. A total of 17 MAPK family genes were identified, which widely distributed on different chromosomes. All the members of the MAPK family could be mainly divided into three sub-classes, among which the members of subclass A and B mainly included PKc_MAPKK_plant_like and PTZ00024 Domain, subclass C mainly contained the PLN00034 domain, and all members contained a S_TKc domain. The results of analysis of physicochemical properties of proteins showed that the number of amino acid residues of each member distributed from 314 to 860 aa, the number of exons was from 1 to 18, the molecular weight of the protein was from 34 910.05 to 97 965.26 u, and the isoelectric point was from 4.94 to 9.35. Specific primer was designed for each PgMAPKs and PgMAPKKs, and their expression patterns were detected. The results showed that 12 out of 17 members were activated after low-temperature treatment. Yichengfenhongmudan that showed strong cold tolerance, and Malisi with weak cold tolerance, were both selected to compare the expression pattern of MAPKs. PgMKK2, PgMPK6, PgMPK9, PgMPK16 and PgMPK13 were all rapidly activated after low-temperature treatment in Yichengfenhongmudan. PgMKK8, PgMPK1-1 and PgMKK4 were rapidly activated in Malisi. PgMKK2, PgMPK6, PgMPK9, PgMPK16 and PgMPK13 showed significantly higher expression level in Yichengfenhongmudan than those in Malisi after low-temperature treatment. The expression of PgMPK3, PgMPK12-1, PgMPK20, PgMPK18 and PgMKK6 was not affected by low-temperature treatment. PgMKK8, PgMPK1-1 and PgMKK4 were up-regulated in Malisi than Yichengfenhongmudan. The expression of PgMPK12-2 was not activated after low-temperature treatment in Malisi, but showed gradual increase in Yichengfenhongmudan. The expression level of PgMKK3 increased at first and then decreased during low-temperature treatment in Malisi, whereas it was not detected in Yichengfenhongmudan. 【Conclusion】 Pomegranate MAPK family genes responded to cold stress signals. PgMKK2, PgMPK6, PgMPK12-2 and PgMPK9 might positively regulate cold tolerance.

    Key words: Pomegranate (Punica granatum L.); Seed hardness; Cold resistance; MAPK; Expression analysis

    石榴(Punica granatum L.)屬千屈菜科(Lythraceae)石榴屬(Punica L.)落葉果樹,是中國重要的經(jīng)濟作物。石榴抗寒性差,硬籽石榴一般遇-17 ℃以下低溫會出現(xiàn)凍害,以突尼斯軟籽為代表的軟籽石榴遇-10 ℃以下低溫即出現(xiàn)凍害[1-2]。軟籽石榴商品價值高,約占石榴總市場的80%,已成為多個石榴主產(chǎn)區(qū)的主栽品種。但是,目前主栽的軟籽石榴突尼斯軟籽和中農(nóng)紅等品種在河南、山東、陜西等主產(chǎn)區(qū)頻發(fā)凍害,給果農(nóng)造成了極大的損失,這已成為制約石榴產(chǎn)業(yè)健康和可持續(xù)發(fā)展的主要問題。

    為緩解石榴凍害問題,國內(nèi)外學(xué)者通過冬季扣棚或埋土、抗寒砧木嫁接和冬季噴施防凍劑等栽培方式防寒[3],并在部分地區(qū)取得了緩解凍害的效果。但與此同時,培育優(yōu)良抗寒軟籽石榴新品種仍是解決凍害問題的根本途徑。目前,國內(nèi)軟籽石榴資源較少,收集或創(chuàng)制并篩選抗寒性強的軟籽石榴資源可加快優(yōu)良抗寒軟籽石榴培育進程。本研究團隊對經(jīng)過多年收集或創(chuàng)制獲得的軟籽石榴種質(zhì)資源進行抗寒性評價,其結(jié)果可為親本選擇提供依據(jù)。石榴抗寒性評價方法主要包括組織褐變法、生理生化指標(biāo)測定法和電導(dǎo)法等,其中電導(dǎo)法是石榴抗寒性最可行的評價方法[4-6]。

    截至目前,已有大量報道植物抗寒調(diào)控關(guān)鍵基因的研究,其中CBF(C-repeat binding factors)基因是調(diào)控植物抗寒性的重要開關(guān)基因[7],而促分裂原活化蛋白激酶(Mitogen-activated protein kinases, MAPK)可通過級聯(lián)反應(yīng)調(diào)控CBF基因的表達(dá)從而調(diào)控植物抗寒性。MAPK是一類保守的絲氨酸-蘇氨酸蛋白激酶,主要通過逐級磷酸化放大和傳遞細(xì)胞外來刺激[8],在植物生長發(fā)育和脅迫響應(yīng)過程中發(fā)揮重要作用[9-12]。MEKK1-MKK1/2-MPK4級聯(lián)反應(yīng)可通過與MPK3和MPK6拮抗,激活CBF基因表達(dá)從而提高植物抗寒性[13-15];MPK3和MPK6等主要通過調(diào)控質(zhì)膜H+-ATP酶活性,磷酸化ICE1蛋白,從而抑制CBF基因表達(dá),負(fù)調(diào)控植物抗寒性[16-18]。

    關(guān)于石榴抗寒性機制研究的相關(guān)報道較少。劉貝貝[19]研究表明,CBF1基因是參與調(diào)控石榴抗寒性的關(guān)鍵基因。同時,MAPK級聯(lián)途徑相關(guān)基因在硬籽群體內(nèi)受強烈選擇,這從進化角度初步解釋了大部分硬籽石榴抗逆能力強于軟籽石榴的原因,暗示了MAPK級聯(lián)途徑對石榴抗逆過程中的作用[20]。Ren等[21]的研究發(fā)現(xiàn),MAPK家族基因可能參與石榴頂端分生組織、花和果實發(fā)育過程。但MAPK家族基因?qū)κ窨购缘挠绊懮胁幻鞔_。

    因此,為研究MAPK家族基因在石榴冷脅迫過程中表達(dá)變化,筆者在本研究中對31份種質(zhì)資源進行籽粒硬度及抗寒性評價。同時,基于全基因組挖掘MAPK家族基因,對MAPK家族基因成員進行系統(tǒng)進化、基因結(jié)構(gòu)和蛋白理化性質(zhì)分析,明確MAPK家族基因在石榴響應(yīng)冷脅迫過程中的表達(dá)模式,進而為石榴抗寒機制的研究提供理論支持。

    1 材料和方法

    1.1 植物材料及試驗地概況

    試驗以31個石榴品種的果實和枝條為試材,進行籽粒硬度和抗寒性檢測。試驗材料均取自國家園藝種質(zhì)資源庫(鄭州)內(nèi)7年生石榴樹,株行距2 m × 4 m,樹體主干開心形,管理方式采用常規(guī)栽培管理,土壤為壤砂土。

    1.2 試驗方法

    1.2.1 不同品種石榴籽粒硬度的測定 成熟期取石榴果實籽粒,每個果實取上中下部各10粒種子混合測定(每個果實測定30個籽粒),每個品種取3個果實分別測定。測定時使用紗布去除籽粒外種皮,取光潔未被破壞的籽粒進行硬度測定,使用TA-XT質(zhì)構(gòu)儀(英國SMS)測定,選擇P2探頭(位移1.3~2.0 mm),測前測后速率5 mm·s-1,測中速率1 mm·s-1。

    1.2.2 不同品種石榴枝條電導(dǎo)率的測定 2021年1月中旬,選取每個品種1年生健壯枝條20~30根,用去離子水洗凈,吸干水分后剪成8~10 cm小段,每組5段,分成3組,捆綁并用保鮮膜包裹。將分組的枝條于4 ℃冰箱預(yù)冷24 h,之后用低溫培養(yǎng)箱進行梯度冷處理,處理溫度分別為-4、-8、-12、-16、-20 ℃,降溫速率4 ℃·h-1,達(dá)到目的溫度后保持24 h,以4 ℃為對照。

    在50 mL離心管中加入25 mL去離子水,將低溫處理后的枝條剪成0.5 cm的小段,各稱取1 g,加入去離子水中,置于25 ℃搖床中,90 min,用電導(dǎo)率儀(雷磁 DDS-307型號)測量初始電導(dǎo)值,然后沸水浴20 min,自然冷卻至室溫,測量其終電導(dǎo)值。相對電導(dǎo)率=初始電導(dǎo)值/終電導(dǎo)值[6]。Logistic回歸方程y=k/(1+ ae-bx),其中y為相對電導(dǎo)率,x為處理溫度,k為當(dāng)x趨于無窮大時的值,a和b為方程參數(shù)。計算半致死溫度(LT50)。

    1.2.3 MAPK家族基因鑒定及生物信息分析 擬南芥MAPK家族的30個成員的氨基酸序列(包括10個MKK和20個MPK成員)來源于擬南芥基因數(shù)據(jù)庫TAIR(https://www.arabidopsis.org/)[22]。以突尼斯軟籽石榴基因組為參考基因組通過NCBI blast+-2.9.0進行本地blast獲得。氨基酸序列通過在線軟件emboss[23](https://www.ebi.ac.uk/Tools/st/emboss_transeq/)翻譯獲得。各編碼蛋白的分子質(zhì)量和理論等電點通過在線軟件ProtParam(https://web.expasy.org/protparam/)進行預(yù)測[24]。MAPK家族基因系統(tǒng)發(fā)育樹的構(gòu)建和多序列比對使用本地軟件MEGA 7.0[25],采用鄰接法,bootstrap值為1000,進化樹可視化通過Evolview在線軟件[26]。氨基酸序列的保守結(jié)構(gòu)域通過在線軟件Batch CD-Search(https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi)進行預(yù)測[27],通過DNAMAN進行MAPK家族基因多序列比對。通過在線網(wǎng)站String(https://cn.string-db.org/)在線網(wǎng)站進行蛋白互作網(wǎng)絡(luò)預(yù)測。

    1.2.4 MAPK響應(yīng)冷脅迫過程中表達(dá)模式分析 取嶧城粉紅牡丹、瑪麗斯1年生枝條,8~10 cm小段,每組5段,分成3組。4 ℃開始降溫,溫度達(dá)到-12 ℃時,分別處理0(對照)、5、10、20、30、60、120 min后取出枝條剪成0.5 cm小塊,液氮速凍,-80 ℃保存。每個處理設(shè)3次生物學(xué)重復(fù)。

    以石榴PgActin為內(nèi)參基因[28],實時熒光定量PCR(real-time quantitative PCR,RT-qPCR)測定各基因在不同脅迫條件下的相對表達(dá)量。利用NCBI在線軟件Primer-BLAST設(shè)計引物(表1)。選用Universal SYBR Green Master熒光定量試劑盒(Roche)在LightCycler? 480 Ⅱ (Roche)進行檢測。Real Time-qPCR反應(yīng)程序:95 ℃ 5 min; 95 ℃ 10 s、56 ℃ 10 s、72 ℃ 10 s,共 45 個循環(huán);95 ℃ 5 s、65 ℃ 1 min,10個循環(huán);40 ℃冷卻。

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

    采用Excel 2013對枝條電導(dǎo)率進行計算和Logistic方程對31個石榴品種的電導(dǎo)率進行回歸分析,熒光定量數(shù)據(jù)計算采用2-ΔΔCT法[29]。采用Excel 2013、SPSS 25.0對籽粒硬度與相對電導(dǎo)率相關(guān)性和熒光定量結(jié)果進行數(shù)據(jù)整理與分析。利用R語言gplots繪制不同基因表達(dá)熱圖。使用T檢驗和多重比較法進行差異顯著性分析。

    2 結(jié)果與分析

    2.1 石榴籽粒硬度與抗寒性分析

    籽粒硬度測定結(jié)果表明,31個石榴品種籽粒硬度范圍主要分布在(1.21±0.48)~(1.96±0.44) kg·cm-2和(4.10±0.80)~(7.29±2.14) kg·cm-2;依據(jù)Zarei等[30]分的籽粒硬度等級標(biāo)準(zhǔn),將籽粒硬度分布在(1.21±0.48)~(1.96±0.44) kg·cm-2的品種歸類為軟籽,將籽粒硬度分布在(4.10±0.80)~(7.29±2.14) kg·cm-2的品種歸類為硬籽;31份石榴種質(zhì)資源共包括17份軟籽和14份硬籽資源(表2)。

    對31個石榴品種半致死溫度測定結(jié)果表明,半致死溫度分布范圍在-13.78~-32.79 ℃;其中,嶧城粉紅牡丹、淮北六棱甜和魯白榴2號等硬籽石榴抗寒性較強,華紫、以3和土庫曼斯坦等軟籽石榴抗寒性較弱;除蒙自甜光顏和蒙自厚皮沙籽外,硬籽石榴半致死溫度低于大部分軟籽石榴;同時,除中石榴4號外,軟籽石榴半致死溫度高于大部分硬籽石榴(表2)。

    2.2 石榴MAPK家族基因的進化分析

    基于擬南芥MAPK家族成員氨基酸序列從突尼斯軟籽石榴基因組數(shù)據(jù)庫中篩選出17個候選的MAPK同源基因(包括5個MAKK和12個MAPK成員)(表3),根據(jù)與擬南芥MAPK家族基因同源性進行命名,廣泛分布于8條染色體上(圖1)。通過擬南芥與石榴MAPK家族基因系統(tǒng)發(fā)育分析,石榴MAPK家族基因可分為3個亞類(A、B、C),其中,A亞類包含PgMPK12-1、PgMPK12-2、PgMPK6、PgMPK3、PgMPK13、PgMPK1-1和PgMPK1-2共7個成員,B亞類包含PgMPK18、PgMPK16、PgMPK20、PgMPK8和PgMPK9共5個成員,C亞類包含PgMKK2、PgMKK3、PgMKK4、PgMKK6和PgMKK8共5個成員(圖2)。

    2.3 石榴MAPK家族的結(jié)構(gòu)分析

    對石榴和擬南芥MAPK家族的功能域進行分析,發(fā)現(xiàn)這些基因均具有促分裂原活化蛋白激酶(Pkinase)的功能結(jié)構(gòu)域,其中,A和B亞類成員主要包含PKc_MAPKK_plant_like、PTZ00024結(jié)構(gòu)域,C亞類主要包含PLN00034結(jié)構(gòu)域(圖2)?;谑馦APK家族基因氨基酸多序列比對結(jié)果,發(fā)現(xiàn)其具有MAPK家族特有的保守結(jié)構(gòu)域S_TKc(圖3)。

    石榴MAPK家族結(jié)構(gòu)分析結(jié)果顯示,石榴MAPK家族各成員氨基酸殘基數(shù)量分布在314~860 aa(表3),外顯子數(shù)目1~18個,同一亞類成員間外顯子數(shù)目、位置及大小類型相近;A亞類成員外顯子數(shù)目主要為6個,PgMAK1-1和PgMAK1-2外顯子數(shù)目為1個,PgMAK12-1外顯子數(shù)目為18個;B亞類成員外顯子數(shù)目主要為10個;C亞類成員外顯子數(shù)目主要為6~9個,PgMKK8外顯子數(shù)目為1個(表3)。利用ProtParam對石榴MAPK家族基因蛋白理化性質(zhì)進行分析,結(jié)果顯示,石榴MAPK家族蛋白分子質(zhì)量為34 910.05~97 965.26 u,等電點4.94~9.35(表3)。

    2.4 石榴MAPK家族基因冷脅迫過程中表達(dá)模式分析

    在抗寒性強的嶧城粉紅牡丹石榴響應(yīng)冷脅迫的過程中,隨著冷脅迫時間推移,MAPK家族基因表達(dá)模式主要分為3類(Ⅰ、Ⅱ和Ⅲ類)。Ⅰ類基因隨著冷脅迫時間推移表達(dá)量表現(xiàn)為先升高后降低的趨勢,主要包括PgMKK2、PgMPK13、PgMPK1-2、PgMPK8、PgMKK4、PgMPK20、PgMPK18共7個基因;Ⅱ類基因隨著冷脅迫時間推移表達(dá)量表現(xiàn)為逐漸升高的趨勢,主要包括PgMPK3、PgMKK6、PgMPK1-1、PgMPK16共4個基因;Ⅲ類基因表達(dá)量表現(xiàn)為逐漸降低的趨勢,主要包括PgMPK12-1、PgMPK12-2、PgMPK6、PgMPK9和PgMKK8共5個基因;PgMKK3在嶧城粉紅牡丹石榴冷脅迫不同時間均未檢測到表達(dá)(圖4-A)。

    在抗寒性差的瑪麗斯石榴響應(yīng)冷脅迫過程中,隨著時間推移,MAPK家族基因表達(dá)模式同樣主要分為3類(Ⅰ、Ⅱ和Ⅲ類)。Ⅰ類主要包括PgMKK3、PgMKK2、PgMPK20、PgMPK13和PgMPK1-1共5個基因;Ⅱ類主要包括PgMPK6、PgMKK6、PgMPK9、PgMPK16、PgMPK18和PgMPK3共6個基因;Ⅲ類主要包括PgMPK8、PgMKK4、PgMKK8、PgMPK12-1和PgMPK1-2共5個基因;PgMPK12-1在瑪麗斯響應(yīng)冷脅迫過程中未檢測到表達(dá)(圖4-B)。

    其中,PgMKK2、PgMPK20和PgMPK13在嶧城粉紅牡丹石榴和瑪麗斯石榴冷脅迫過程中表達(dá)量均表現(xiàn)為先升高后降低的趨勢;PgMPK3、PgMKK6和PgMPK16表達(dá)量均表現(xiàn)為逐漸降低的趨勢;PgMPK12-1和PgMKK8表達(dá)量均表現(xiàn)為逐漸升高的趨勢(圖4)。PgMPK9和PgMPK6在嶧城粉紅牡丹石榴響應(yīng)冷脅迫過程中表現(xiàn)為逐漸升高的趨勢,在瑪麗斯石榴中表現(xiàn)為逐漸降低的趨勢;PgMPK18在嶧城粉紅牡丹石榴不同時間表現(xiàn)為先升高后降低的趨勢,在瑪麗斯石榴中表現(xiàn)為逐漸降低的趨勢;PgMPK1-2、PgMPK8和PgMKK4在嶧城粉紅牡丹石榴不同時間表現(xiàn)為先升高后降低的趨勢,在瑪麗斯石榴中表現(xiàn)為逐漸升高的趨勢;PgMPK1-1在嶧城粉紅牡丹石榴不同時間表現(xiàn)為逐漸降低,在瑪麗斯石榴中表現(xiàn)為先升高后降低;PgMKK3在嶧城粉紅牡丹石榴不同時間均未檢測到表達(dá),在瑪麗斯石榴中表現(xiàn)為先升高后降低的趨勢;PgMPK12-2在瑪麗斯石榴不同時間均未檢測到表達(dá),在嶧城粉紅牡丹石榴中表現(xiàn)為逐漸升高的趨勢(圖4)。

    MAPK家族基因在嶧城粉紅牡丹石榴和瑪麗斯石榴間的相對表達(dá)量差異結(jié)果表明,PgMKK2、PgMPK6、PgMPK9、PgMPK16在響應(yīng)冷脅迫過程中,在嶧城粉紅牡丹石榴中的相對表達(dá)量顯著高于瑪麗斯石榴;PgMPK13相對表達(dá)量在響應(yīng)脅迫后5和10 min時在嶧城粉紅牡丹石榴和瑪麗斯石榴間不存在顯著性差異表達(dá),20、30和60 min時在嶧城粉紅牡丹石榴中的相對表達(dá)量顯著高于瑪麗斯石榴;PgMKK8和PgMPK1-1在瑪麗斯石榴中的相對表達(dá)量顯著性高于嶧城粉紅牡丹石榴;PgMKK6、PgMPK12-1和PgMPK3的相對表達(dá)量在瑪麗斯和嶧城粉紅牡丹石榴間無顯著差異;PgMKK4、PgMPK8、PgMPK1-2響應(yīng)冷脅迫前期在嶧城粉紅牡丹石榴中的相對表達(dá)量顯著高于瑪麗斯石榴、后期在瑪麗斯石榴中的相對表達(dá)量顯著高于嶧城粉紅牡丹石榴;PgMPK20和PgMPK18的相對表達(dá)量與石榴對冷脅迫的響應(yīng)無顯著相關(guān)性(圖5)。

    3 討 論

    3.1 枝條取樣時間不同影響檢測的半致死溫度

    筆者在本研究中測定的不同石榴品種半致死溫度范圍在-13.77~-32.79 ℃,與羅華等[31]發(fā)表的石榴半致死溫度存在范圍差異(-9.04~-11.77 ℃)。Ghasemi等[4]對不同取樣時間的石榴枝條半致死溫度檢測的差異進行分析,發(fā)現(xiàn)同一品種1月份取石榴枝條進行測定獲得的半致死溫度顯著低于3月份取枝條所獲得的半致死溫度。筆者取樣時間為石榴枝條休眠期(1月中旬),羅華等[31]取樣時間為萌芽前(3月中旬),這一結(jié)論解釋了本研究測定的石榴半致死溫度范圍顯著低于其發(fā)表范圍的原因。但是,半致死溫度的范圍不影響品種間抗寒性強弱的差異比較。

    3.2 石榴抗寒性可能與籽粒硬度呈正相關(guān)

    大部分硬籽石榴果實內(nèi)種皮次生細(xì)胞壁厚度顯著高于軟籽石榴,枝條較軟籽石榴硬。石榴抗寒性和籽粒硬度性狀可能存在連鎖遺傳現(xiàn)象。例如,NAC家族基因通過調(diào)控木質(zhì)素、纖維素和半纖維素合成從而參與調(diào)控石榴籽粒硬度形成[28],在硬籽石榴中的表達(dá)量顯著高于軟籽石榴。同時,NAC家族基因參與植物對抗寒性等非生物脅迫的響應(yīng)[32]。同時,已有研究報道生長環(huán)境溫度對植物抗寒性具有顯著影響,來源于云南產(chǎn)區(qū)的蒙自甜光顏石榴和蒙自滑皮沙籽石榴抗寒性、顯著低于新疆、山東、河南和安徽產(chǎn)區(qū)的石榴,長期自然環(huán)境馴化對石榴抗寒性具有顯著影響[33]。

    3.3 石榴MAPK家族基因響應(yīng)冷脅迫

    目前,水稻[34]、玉米[35]、葡萄[36]、棗[37]、麻風(fēng)樹[38]等多個物種上MAPK家族基因被挖掘,在植物生長發(fā)育、抗逆等過程中發(fā)揮重要功能。筆者在本研究中分析石榴MAPK家族基因功能性結(jié)構(gòu)域,發(fā)現(xiàn)在不同物種間MAPK家族基因功能性結(jié)構(gòu)域具有保守性。在擬南芥和玉米中MAPK家族基因MPK3、MPK4、MPK6和MKK2等在冷處理30 min內(nèi)表達(dá)量顯著上調(diào)[34-35]。MKK2基因在馬鈴薯響應(yīng)冷脅迫過程中顯著上調(diào)表達(dá)[39]。MPK5[40]和MPK3[41]響應(yīng)香蕉冷脅迫處理過程,并通過調(diào)控NAC042和ICE1基因表達(dá)從而參與調(diào)控香蕉抗寒性。石榴17個MAPK家族基因中12個均能不同程度地響應(yīng)石榴冷脅迫信號。PgMKK2、PgMPK6、PgMPK9、PgMPK16、PgMPK13在嶧城粉紅牡丹石榴響應(yīng)冷脅迫過程中表現(xiàn)為顯著上調(diào),PgMKK8、PgMPK1-1和PgMKK4在瑪麗斯石榴響應(yīng)冷脅迫過程中表現(xiàn)為顯著上調(diào),MAPK家族基因參與石榴響應(yīng)冷脅迫過程。

    3.4 PgMKK2-MPK6級聯(lián)反應(yīng)途徑可能參與調(diào)控石榴抗寒性

    已有研究表明,低溫誘導(dǎo)MKK2蛋白磷酸化[15],MKK2激發(fā)MPK4/MPK6蛋白磷酸化[14],MPK4/MPK6通過調(diào)控ICE1和CBF基因表達(dá)從而調(diào)控植物抗寒性[14];同時MKK2-MAPK4/7通過調(diào)控ICE1基因表達(dá)從而參與調(diào)控植物抗寒性;MEKK1-MKK2-MPK4通過調(diào)控CBFs基因表達(dá)從而正調(diào)控植物抗寒性。而低溫處理可顯著激活石榴PgMKK2和PgMPK6基因表達(dá),冷脅迫處理后PgMKK2和PgMPK6在嶧城粉紅牡丹石榴中的表達(dá)量顯著高于瑪麗斯石榴,同時,PgMPK6表達(dá)趨勢與PgMKK2一致,推測PgMPK6基因表達(dá)上調(diào)可能與PgMKK2級聯(lián)反應(yīng)相關(guān),PgMKK2-MPK6級聯(lián)反應(yīng)途徑可能參與正調(diào)控石榴抗寒性。

    3.5 PgMPK12-2可能參與正調(diào)控石榴抗寒性

    油菜、擬南芥等物種中MPK12與MPK9功能冗余,正向調(diào)控脫落酸、茉莉酸、水楊酸[42-43]、低溫和鹽[44]等逆境脅迫過程。而石榴PgMPK12-2在瑪麗斯石榴不同時間均未檢測到表達(dá),在嶧城粉紅牡丹石榴中表現(xiàn)為逐漸升高的趨勢。PgMPK9在響應(yīng)冷脅迫過程中在嶧城粉紅牡丹石榴中的表達(dá)量顯著高于瑪麗斯石榴。PgMPK12-2和PgMPK9可能參與石榴響應(yīng)低溫脅迫的過程。

    4 結(jié) 論

    筆者在本研究中基于石榴全基因組共挖掘到17個MAPK家族基因成員,分布于不同染色體,根據(jù)系統(tǒng)發(fā)育樹將其分為3個亞類,所有成員均含有S-TKc結(jié)構(gòu)域,同一亞類成員間保守結(jié)構(gòu)域、外顯子數(shù)目具有保守性。在響應(yīng)冷脅迫過程中,PgMKK2、PgMPK6、PgMPK9、PgMPK16和PgMAPK13在嶧城粉紅牡丹石榴中的表達(dá)量顯著高于瑪麗斯石榴,PgMKK8和PgMPK1-1在瑪麗斯石榴中的表達(dá)量顯著高于嶧城粉紅牡丹石榴,PgMKK2、PgMPK6、PgMPK12-2和PgMPK9可能正調(diào)控石榴抗寒性。

    參考文獻 References:

    [1] 楊雪梅,苑兆和,尹燕雷,馮立娟,招雪晴. 不同石榴品種抗寒性綜合評價[J]. 山東農(nóng)業(yè)科學(xué),2014,46(2):46-51.

    YANG Xuemei,YUAN Zhaohe,YIN Yanlei,F(xiàn)ENG Lijuan,ZHAO Xueqing. Comprehensive evaluation on cold resistance of different pomegranate cultivars[J]. Shandong Agricultural Sciences,2014,46(2):46-51.

    [2] 姚方,王寧,曹尚銀,馬貫羊,司守霞,姚海雷.不同軟籽品種石榴抗寒性綜合評價[J]. 森林與環(huán)境學(xué)報,2016,36(3):373-379.

    YAO Fang,WANG Ning,CAO Shangyin,MA Guanyang,SI Shouxia,YAO Hailei.Comprehensive evaluation on cold resistance of different pomegranate cultivars[J]. Journal of Forest and Environment,2016,36(3):373-379.

    [3] 焦其慶,馮立娟,尹燕雷,崔洪濤. 石榴凍害及抗寒評價研究進展[J]. 植物生理學(xué)報,2019,55(4):425-432.

    JIAO Qiqing,F(xiàn)ENG Lijuan,YIN Yanlei,CUI Hongtao. Research progress on evaluation of freezing injury and cold resistance of pomegranate[J]. Plant Physiology Journal,2019,55(4):425-432.

    [4] SOLOKLUI A A G,ERSHADI A,F(xiàn)ALLAHI E. Evaluation of cold hardiness in seven Iranian commercial pomegranate (Punica granatum L.) cultivars[J]. HortScience,2012,47(12):1821-1825.

    [5] 畢潤霞,郝兆祥,侯樂峰,王慶軍,李昭慧,王艷琴. 電導(dǎo)法評價石榴抗寒性方法的探討[J]. 山東農(nóng)業(yè)科學(xué),2015,47(2):38-41.

    BI Runxia,HAO Zhaoxiang,HOU Lefeng,WANG Qingjun,LI Zhaohui,WANG Yanqin. Evaluation on cold resistance of pomegranate by conductivity method[J]. Shandong Agricultural Sciences,2015,47(2):38-41.

    [6] 劉貝貝,陳利娜,牛娟,李好先,張杰,曹尚銀. 6個石榴品種抗寒性評價及方法篩選[J]. 果樹學(xué)報,2018,35(1):66-73.

    LIU Beibei,CHEN Lina,NIU Juan,LI Haoxian,ZHANG Jie,CAO Shangyin. Selection of methods for evaluation on cold tolerance of six pomegranate varieties[J]. Journal of Fruit Science,2018,35(1):66-73.

    [7] ZHANG F,JIANG Y,BAI L P,ZHANG L,CHEN L J,LI H G,YIN Y,YAN W W,YI Y,GUO Z F. The ICE-CBF-COR pathway in cold acclimation and AFPs in plants[J]. Middle-East Journal of Scientific Research,2011,8(2):493-498.

    [8] TIAN K,LI Q,ZHANG X M,GUO H Y,WANG Y H,CAO P L,XU S Y,LI W Y. Analysis of the expression and function of the CBL-CIPK network and MAPK cascade genes in Kandelia obovata seedlings under cold stress[J]. Frontiers in Marine Science,2023,10:1113278.

    [9] PITZSCHKE A,SCHIKORA A,HIRT H. MAPK cascade signalling networks in plant defence[J]. Current Opinion in Plant Biology,2009,12(4):421-426.

    [10] BERGMANN D C,LUKOWITZ W,SOMERVILLE C R. Stomatal development and pattern controlled by a MAPKK kinase[J]. Science,2004,304(5676):1494-1497.

    [11] MOUSTAFA K,ABUQAMAR S,JARRAR M,AL-RAJAB A J,TR?MOUILLAUX-GUILLER J. MAPK cascades and major abiotic stresses[J]. Plant Cell Reports,2014,33(8):1217-1225.

    [12] ZHANG M M,SU J B,ZHANG Y,XU J,ZHANG S Q. Conveying endogenous and exogenous signals:MAPK cascades in plant growth and defense[J]. Current Opinion in Plant Biology,2018,45:1-10.

    [13] BANERJEE G,SINGH D,SINHA A K. Plant cell cycle regulators:Mitogen-activated protein kinase,a new regulating switch?[J]. Plant Science,2020,301:110660.

    [14] TEIGE M,SCHEIKL E,EULGEM T,D?CZI R,ICHIMURA K,SHINOZAKI K,DANGL J L,HIRT H. The MKK2 pathway mediates cold and salt stress signaling in Arabidopsis[J]. Molecular Cell,2004,15(1):141-152.

    [15] FURUYA T,MATSUOKA D,NANMORI T. Phosphorylation of Arabidopsis thaliana MEKK1 via Ca2+ signaling as a part of the cold stress response[J]. Journal of Plant Research,2013,126(6):833-840.

    [16] LIU Y K,ZHOU J. MAPping kinase regulation of ICE1 in freezing tolerance[J]. Trends in Plant Science,2018,23(2):91-93.

    [17] ZHAO C Z,WANG P C,SI T,HSU C C,WANG L,ZAYED O,YU Z P,ZHU Y F,DONG J,TAO W A,ZHU J K. MAP kinase cascades regulate the cold response by modulating ICE1 protein stability[J]. Developmental Cell,2017,43(5):618-629.

    [18] PONCE-PINEDA I G,CARMONA-SALAZAR L,SAUCEDO-GARC?A M,CANO-RAM?REZ D,MORALES-CEDILLO F,PE?A-MORAL A,GUEVARA-GARC?A ? A,S?NCHEZ-NIETO S,GAVILANES-RU?Z M. MPK6 kinase regulates plasma membrane H+-ATPase activity in cold acclimation[J]. International Journal of Molecular Sciences,2021,22(12):6338.

    [19] 劉貝貝. 石榴抗寒品種篩選及轉(zhuǎn)錄因子PgCBF1功能分析[D]. 北京:中國農(nóng)業(yè)科學(xué)院,2018.

    LIU Beibei. Screening cold-tolerante cultivars of pomegranate and analyzing the function of PgCBF1 transcription factor[D]. Beijing:Chinese Academy of Agricultural Sciences,2018.

    [20] LUO X,LI H X,WU Z K,YAO W,ZHAO P,CAO D,YU H Y,LI K D,POUDEL K,ZHAO D G,ZHANG F H,XIA X C,CHEN L N,WANG Q,JING D,CAO S Y. The pomegranate (Punica granatum L.) draft genome dissects genetic divergence between soft-and hard-seeded cultivars[J]. Plant Biotechnology Journal,2020,18(4):955-968.

    [21] REN Y A,GE D P,DONG J M,GUO L H,YUAN Z H. Genome-wide identification and expression analysis of MAPK and MAPKK gene family in pomegranate (Punica granatum L.)[J]. Agronomy,2020,10(7):1015.

    [22] LAMESCH P,BERARDINI T Z,LI D H,SWARBRECK D,WILKS C,SASIDHARAN R,MULLER R,DREHER K,ALEXANDER D L,GARCIA-HERNANDEZ M,KARTHIKEYAN A S,LEE C H,NELSON W D,PLOETZ L,SINGH S,WENSEL A,HUALA E. The Arabidopsis Information Resource (TAIR):Improved gene annotation and new tools[J]. Nucleic Acids Research,2012,40(Database issue):1202-1210.

    [23] MADEIRA F,PEARCE M,TIVEY A R N,BASUTKAR P,LEE J,EDBALI O,MADHUSOODANAN N,KOLESNIKOV A,LOPEZ R. Search and sequence analysis tools services from EMBL-EBI in 2022[J]. Nucleic Acids Research,2022,50(W1):276-279.

    [24] ARTIMO P,JONNALAGEDDA M,ARNOLD K,…,STOCKINGER H. ExPASy:SIB bioinformatics resource portal[J]. Nucleic Acids Research,2012,40(Web Server issue):597-603.

    [25] KUMAR S,STECHER G,TAMURA K. MEGA7:Molecular evolutionary genetics analysis version 7.0 for bigger datasets[J]. Molecular Biology and Evolution,2016,33(7):1870-1874.

    [26] HE Z L,ZHANG H K,GAO S H,LERCHER M J,CHEN W H,HU S N. Evolview v2:An online visualization and management tool for customized and annotated phylogenetic trees[J]. Nucleic Acids Research,2016,44(W1):236-241.

    [27] LU S N,WANG J Y,CHITSAZ F,DERBYSHIRE M K,GEER R C,GONZALES N R,GWADZ M,HURWITZ D I,MARCHLER G H,SONG J S,THANKI N,YAMASHITA R A,YANG M Z,ZHANG D C,ZHENG C J,LANCZYCKI C J,MARCHLER-BAUER A. CDD/SPARCLE:The conserved domain database in 2020[J]. Nucleic Acids Research,2020,48(D1):265-268.

    [28] QIN G H,LIU C Y,LI J Y,QI Y J,GAO Z H,ZHANG X L,YI X K,PAN H F,MING R,XU Y L. Diversity of metabolite accumulation patterns in inner and outer seed Coats of pomegranate:Exploring their relationship with genetic mechanisms of seed coat development[J]. Horticulture Research,2020,7:10.

    [29] LIVAK K J,SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method[J]. Methods,2001,25(4):402-408.

    [30] ZAREI A, ZAMANI Z, FATAHI R, MOUSAVI A, SALAMI S A, AVILA C , C?NOVAS F M. Differential expression of cell wall related genes in the seeds of soft- and hard-seeded pomegranate genotypes[J]. Scientia Horticulturae, 2016, 205: 7-16.

    [31] 羅華,王慶軍,郝兆祥,趙麗娜,陳穎,趙登超,畢潤霞,孟健,馬敏. 石榴抗寒種質(zhì)篩選研究[J]. 中國果樹,2018(4):46-52.

    LUO Hua,WANG Qingjun,HAO Zhaoxiang,ZHAO Lina,CHEN Ying,ZHAO Dengchao,BI Runxia,MENG Jian,MA Min. Study on cold resistance of pomegranate germplasms[J]. China Fruits,2018(4):46-52.

    [32] SINGH S,KOYAMA H,BHATI K K,ALOK A. The biotechnological importance of the plant-specific NAC transcription factor family in crop improvement[J]. Journal of Plant Research,2021,134(3):475-495.

    [33] DUNN S. Factors affecting cold resistance in plants[J]. Plant Physiology,1937,12(2):519-526.

    [34] RAO K P,RICHA T,KUMAR K,RAGHURAM B,SINHA A K. In silico analysis reveals 75 members of mitogen-activated protein kinase kinase kinase gene family in rice[J]. DNA Research,2010,17(3):139-153.

    [35] KONG X P,L? W,ZHANG D,JIANG S S,ZHANG S Z,LI D Q. Genome-wide identification and analysis of expression profiles of maize mitogen-activated protein kinase kinase kinase[J]. PLoS One,2013,8(2):e57714.

    [36] ?AKIR B,KILI?KAYA O. Mitogen-activated protein kinase cascades in Vitis vinifera[J]. Frontiers in Plant Science,2015,6:556.

    [37] LIU Z G,ZHANG L M,XUE C L,F(xiàn)ANG H,ZHAO J,LIU M J. Genome-wide identification and analysis of MAPK and MAPKK gene family in Chinese jujube (Ziziphus jujuba Mill.)[J]. BMC Genomics,2017,18(1):855.

    [38] WANG H B,GONG M,GUO J Y,XIN H,GAO Y,LIU C,DAI D Q,TANG L Z. Genome-wide identification of Jatropha curcas MAPK,MAPKK,and MAPKKK gene families and their expression profile under cold stress[J]. Scientific Reports,2018,8:16163.

    [39] CHEN Y,CHEN L,SUN X M,KOU S,LIU T T,DONG J K,TU W,ZHANG Y L,SONG B T. The mitogen-activated protein kinase kinase MKK2 positively regulates constitutive cold resistance in the potato[J]. Environmental and Experimental Botany,2022,194:104702.

    [40] TAK H,NEGI S,RAJPUROHIT Y S,RAJPUROHIT Y S,MISRA H S,GANAPATHI T R. MusaMPK5,a mitogen activated protein kinase is involved in regulation of cold tolerance in banana[J]. Plant Physiology and Biochemistry,2020,146:112-123.

    [41] GAO J,DOU T X,HE W D,SHENG O,BI F C,DENG G M,GAO H J,DONG T,LI C Y,ZHANG S,YI G J,HU C H,YANG Q S. MaMAPK3-MaICE1-MaPOD P7 pathway,a positive regulator of cold tolerance in banana[J]. BMC Plant Biology,2021,21(1):97.

    [42] JAMMES F,SONG C,SHIN D,MUNEMASA S,TAKEDA K,GU D,CHO D S,LEE S M,GIORDO R,SRITUBTIM S,LEONHARDT N,ELLIS B E,MURATA Y,KWAK J M. MAP kinases MPK9 and MPK12 are preferentially expressed in guard cells and positively regulate ROS-mediated ABA signaling[J]. Proceedings of the National Academy of Sciences of the United States of America,2009,106(48):20520-20525.

    [43] KHOKON M A R,SALAM M A,JAMMES F,YE W X,HOSSAIN M A,OKUMA E,NAKAMURA Y,MORI I C,KWAK J M,MURATA Y. MPK9 and MPK12 function in SA-induced stomatal closure in Arabidopsis thaliana[J]. Bioscience,Biotechnology,and Biochemistry,2017,81(7):1394-1400.

    [44] 李萍,聶亭亭,張騰國,鄭晟,王娟,毛玉珊. 白菜型油菜MPK12基因克隆及表達(dá)分析[J]. 西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版),2017,45(5):26-33.

    LI Ping,NIE Tingting,ZHANG Tengguo,ZHENG Sheng,WANG Juan,MAO Yushan. Clone and characterization of MPK12 gene from Brassica campestris[J]. Journal of Northwest A & F University (Natural Science Edition),2017,45(5):26-33.

    收稿日期:2023-04-12 接受日期:2023-06-27

    基金項目:國家重點研發(fā)計劃項目(2021YFD1600802);中國農(nóng)業(yè)科學(xué)院科技創(chuàng)新工程(CAAS-ASTIP-2023-ZFRI);河南省科技攻關(guān)項目(212102110421)

    作者簡介:陳利娜,女,助理研究員,研究方向為果樹遺傳育種。Tel:13283811852, E-mail:1571863765@qq.com。#為共同第一作者。

    通信作者 Author for correspondence. Tel:0371-65330990,E-mail:luzhenhua@caas.cn

    猜你喜歡
    表達(dá)分析抗寒性石榴
    石榴紅了
    石榴籽
    草原歌聲(2021年1期)2021-07-16 06:58:06
    緊緊抱在一起的“石榴籽”
    少先隊活動(2021年1期)2021-03-29 05:27:06
    石榴
    紅花生育酚環(huán)化酶基因的克隆及表達(dá)分析
    棗樹抗寒性檢測方法的篩選
    巴梨的抗寒性試驗情況初報
    中國果菜(2016年9期)2016-03-01 01:28:40
    馬鈴薯普通栽培種雜交后代抗寒性分析
    中國馬鈴薯(2015年5期)2016-01-09 06:11:11
    玉米紋枯病病菌y—谷氨酰轉(zhuǎn)肽酶基因克隆與表達(dá)分析
    信號分子與葉銹菌誘導(dǎo)下小麥病程相關(guān)蛋白1基因的表達(dá)分析
    国产精品女同一区二区软件 | 日本黄大片高清| 99热6这里只有精品| 国产精品伦人一区二区| 精品一区二区免费观看| 久久久久九九精品影院| 99精品久久久久人妻精品| 精品久久久久久久久久免费视频| 久久久久久国产a免费观看| 婷婷精品国产亚洲av在线| 成人国产综合亚洲| 成人欧美大片| 又爽又黄a免费视频| 久久久久性生活片| 禁无遮挡网站| 亚洲av一区综合| 久久99热这里只有精品18| 免费在线观看影片大全网站| 一区二区三区高清视频在线| 99久久精品一区二区三区| 中国美女看黄片| 十八禁网站免费在线| 中文字幕人成人乱码亚洲影| 成人av一区二区三区在线看| av福利片在线观看| 精品福利观看| 国内精品一区二区在线观看| 在线播放无遮挡| 国产伦精品一区二区三区四那| 国产aⅴ精品一区二区三区波| 日韩欧美国产在线观看| 亚洲不卡免费看| 久久久久性生活片| 精品久久久久久久人妻蜜臀av| 熟女电影av网| 亚洲av免费在线观看| 欧美激情在线99| 又爽又黄无遮挡网站| 人人妻人人看人人澡| 午夜两性在线视频| 我的女老师完整版在线观看| 午夜老司机福利剧场| 精品午夜福利视频在线观看一区| 久久久久国产精品人妻aⅴ院| 欧美又色又爽又黄视频| 国产男靠女视频免费网站| 九九在线视频观看精品| 国产激情偷乱视频一区二区| 窝窝影院91人妻| 波野结衣二区三区在线| 婷婷精品国产亚洲av在线| 欧美黄色片欧美黄色片| 91在线精品国自产拍蜜月| 欧美日韩福利视频一区二区| 别揉我奶头 嗯啊视频| 中文字幕高清在线视频| 小说图片视频综合网站| 亚洲自拍偷在线| 欧美精品国产亚洲| 色播亚洲综合网| 亚洲人成伊人成综合网2020| 久久人妻av系列| 亚洲国产色片| 久久天躁狠狠躁夜夜2o2o| 90打野战视频偷拍视频| 久久精品综合一区二区三区| 淫秽高清视频在线观看| 国产av不卡久久| 国产一区二区在线av高清观看| 国产一区二区在线av高清观看| 舔av片在线| 深爱激情五月婷婷| 在线观看美女被高潮喷水网站 | 久久精品国产亚洲av涩爱 | 一个人看视频在线观看www免费| 欧美+日韩+精品| 成人高潮视频无遮挡免费网站| 婷婷亚洲欧美| 一级黄色大片毛片| 熟女电影av网| 麻豆成人av在线观看| 国产精品99久久久久久久久| av专区在线播放| 欧美日韩瑟瑟在线播放| 国产毛片a区久久久久| 757午夜福利合集在线观看| 国产精品人妻久久久久久| 又紧又爽又黄一区二区| 久久久久亚洲av毛片大全| 午夜福利欧美成人| 变态另类丝袜制服| 天美传媒精品一区二区| 天美传媒精品一区二区| 青草久久国产| 中出人妻视频一区二区| 成人美女网站在线观看视频| 亚洲,欧美,日韩| 999久久久精品免费观看国产| 久久国产乱子免费精品| 哪里可以看免费的av片| 性色av乱码一区二区三区2| 日韩亚洲欧美综合| x7x7x7水蜜桃| 中文字幕av成人在线电影| 白带黄色成豆腐渣| 国产一级毛片七仙女欲春2| 日韩欧美国产一区二区入口| 又紧又爽又黄一区二区| 每晚都被弄得嗷嗷叫到高潮| 90打野战视频偷拍视频| 深夜精品福利| a在线观看视频网站| 国产在线精品亚洲第一网站| 久久热精品热| 国产精品精品国产色婷婷| xxxwww97欧美| 婷婷亚洲欧美| 欧美精品国产亚洲| 99在线视频只有这里精品首页| 99久国产av精品| 人妻久久中文字幕网| 级片在线观看| 免费一级毛片在线播放高清视频| 欧美+日韩+精品| 两个人的视频大全免费| 无遮挡黄片免费观看| 热99在线观看视频| 我要看日韩黄色一级片| 日韩高清综合在线| 99热这里只有是精品在线观看 | 免费看a级黄色片| 国产私拍福利视频在线观看| 免费人成视频x8x8入口观看| 国产色爽女视频免费观看| 不卡一级毛片| 听说在线观看完整版免费高清| 夜夜夜夜夜久久久久| 成年版毛片免费区| 麻豆国产97在线/欧美| 少妇人妻精品综合一区二区 | 亚洲精品影视一区二区三区av| 精品人妻熟女av久视频| 精品久久久久久久久久免费视频| 亚洲色图av天堂| 老熟妇仑乱视频hdxx| 久久久成人免费电影| 亚洲七黄色美女视频| 极品教师在线免费播放| or卡值多少钱| 亚洲av免费高清在线观看| 欧美黄色淫秽网站| 国内精品久久久久精免费| 国产白丝娇喘喷水9色精品| 91字幕亚洲| 国产精品久久久久久精品电影| 黄色视频,在线免费观看| 俺也久久电影网| 欧美极品一区二区三区四区| 18+在线观看网站| 成人特级黄色片久久久久久久| 久久精品久久久久久噜噜老黄 | 欧美黄色片欧美黄色片| 麻豆国产97在线/欧美| 十八禁人妻一区二区| 精品久久久久久久久久久久久| 一个人看的www免费观看视频| 自拍偷自拍亚洲精品老妇| 美女高潮喷水抽搐中文字幕| 国产精品一区二区三区四区久久| 国产成人影院久久av| 国产v大片淫在线免费观看| 人人妻人人看人人澡| 日本黄大片高清| 亚洲精品日韩av片在线观看| 国产高清视频在线观看网站| 男女视频在线观看网站免费| 午夜福利成人在线免费观看| 少妇高潮的动态图| 午夜福利18| 久久久精品欧美日韩精品| 亚洲人成电影免费在线| 窝窝影院91人妻| 大型黄色视频在线免费观看| 亚洲熟妇中文字幕五十中出| 毛片女人毛片| 亚洲在线观看片| 国产精品一区二区三区四区免费观看 | 国产亚洲精品久久久久久毛片| 五月伊人婷婷丁香| 国产精品日韩av在线免费观看| 免费人成视频x8x8入口观看| 男人狂女人下面高潮的视频| 欧美+日韩+精品| 久久国产乱子伦精品免费另类| 免费观看人在逋| 成年人黄色毛片网站| 欧美一级a爱片免费观看看| 国产一级毛片七仙女欲春2| 日韩欧美一区二区三区在线观看| 日本成人三级电影网站| 三级国产精品欧美在线观看| 久久久成人免费电影| 一个人观看的视频www高清免费观看| 国产精品久久久久久久久免 | 国产视频一区二区在线看| 国内揄拍国产精品人妻在线| 日韩欧美免费精品| 亚洲欧美激情综合另类| 两个人视频免费观看高清| 中国美女看黄片| 精品午夜福利视频在线观看一区| 搡老熟女国产l中国老女人| 不卡一级毛片| 国产精品日韩av在线免费观看| 99久久无色码亚洲精品果冻| 午夜精品在线福利| 成人性生交大片免费视频hd| 又黄又爽又刺激的免费视频.| 国产午夜福利久久久久久| 亚洲精品色激情综合| 国产精品嫩草影院av在线观看 | 校园春色视频在线观看| 久久久久久久久中文| 99久久久亚洲精品蜜臀av| 欧美性猛交黑人性爽| 直男gayav资源| 久久久久久九九精品二区国产| 亚洲av熟女| 亚洲三级黄色毛片| 神马国产精品三级电影在线观看| 淫秽高清视频在线观看| 午夜老司机福利剧场| 人妻久久中文字幕网| 亚洲专区国产一区二区| 久久久久久久久久黄片| 久久午夜亚洲精品久久| 成人精品一区二区免费| 精品一区二区三区人妻视频| 国产精品一区二区三区四区久久| 51国产日韩欧美| 亚洲一区二区三区不卡视频| 亚洲片人在线观看| 最近最新免费中文字幕在线| 美女高潮的动态| 久久久久免费精品人妻一区二区| av在线观看视频网站免费| 极品教师在线免费播放| 少妇人妻精品综合一区二区 | 成人av一区二区三区在线看| 一级a爱片免费观看的视频| 最近中文字幕高清免费大全6 | 国产麻豆成人av免费视频| h日本视频在线播放| 亚洲最大成人av| 老女人水多毛片| 欧美在线黄色| 国产69精品久久久久777片| 午夜两性在线视频| 亚洲七黄色美女视频| 欧美一级a爱片免费观看看| 日韩欧美国产一区二区入口| 亚洲av.av天堂| 欧美性猛交╳xxx乱大交人| 可以在线观看的亚洲视频| 成人亚洲精品av一区二区| 免费看光身美女| 国产高潮美女av| 琪琪午夜伦伦电影理论片6080| 欧美乱妇无乱码| 亚洲久久久久久中文字幕| 欧美日韩瑟瑟在线播放| 2021天堂中文幕一二区在线观| 最近中文字幕高清免费大全6 | 国产亚洲欧美在线一区二区| 老女人水多毛片| 一个人看的www免费观看视频| 波多野结衣高清作品| 非洲黑人性xxxx精品又粗又长| 日韩精品中文字幕看吧| 亚洲av免费在线观看| 狂野欧美白嫩少妇大欣赏| 91在线观看av| 少妇人妻一区二区三区视频| 久久6这里有精品| 亚洲天堂国产精品一区在线| 成人特级黄色片久久久久久久| 女生性感内裤真人,穿戴方法视频| 桃色一区二区三区在线观看| netflix在线观看网站| 亚洲精品一区av在线观看| 国内毛片毛片毛片毛片毛片| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品一卡2卡三卡4卡5卡| 97超视频在线观看视频| www.www免费av| 最近最新中文字幕大全电影3| 久久精品91蜜桃| 久久人人爽人人爽人人片va | 在线播放无遮挡| 国产男靠女视频免费网站| av视频在线观看入口| 亚洲国产精品久久男人天堂| 精品国产亚洲在线| 国产一区二区亚洲精品在线观看| 亚洲av免费高清在线观看| 搡老熟女国产l中国老女人| 国产欧美日韩一区二区三| 性欧美人与动物交配| 午夜久久久久精精品| 免费无遮挡裸体视频| 国产欧美日韩精品亚洲av| 一本精品99久久精品77| 久久精品国产自在天天线| 两人在一起打扑克的视频| 日韩人妻高清精品专区| 久久性视频一级片| 午夜两性在线视频| 午夜福利视频1000在线观看| 精品国产亚洲在线| 成年人黄色毛片网站| 麻豆成人av在线观看| 国产一区二区亚洲精品在线观看| 搡老岳熟女国产| 12—13女人毛片做爰片一| 久久亚洲精品不卡| 非洲黑人性xxxx精品又粗又长| 老熟妇乱子伦视频在线观看| 国产精品一区二区三区四区久久| 日韩 亚洲 欧美在线| 乱码一卡2卡4卡精品| 亚洲专区中文字幕在线| 制服丝袜大香蕉在线| 乱码一卡2卡4卡精品| 欧美色欧美亚洲另类二区| 亚洲成人久久性| 欧美成狂野欧美在线观看| 在线天堂最新版资源| 97热精品久久久久久| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av电影在线进入| 亚洲av中文字字幕乱码综合| 国产一区二区激情短视频| 国产探花在线观看一区二区| 久久精品国产清高在天天线| 欧美日韩亚洲国产一区二区在线观看| 久久国产乱子免费精品| 成熟少妇高潮喷水视频| avwww免费| 一卡2卡三卡四卡精品乱码亚洲| 人人妻人人澡欧美一区二区| 亚洲成a人片在线一区二区| 久久精品人妻少妇| 国产美女午夜福利| 国产精品久久久久久亚洲av鲁大| 欧美日韩福利视频一区二区| 波多野结衣高清无吗| 能在线免费观看的黄片| 国产精品久久久久久亚洲av鲁大| 亚洲成a人片在线一区二区| 国产精品影院久久| 99久久精品国产亚洲精品| 午夜免费成人在线视频| 国产成人欧美在线观看| av在线观看视频网站免费| 日本撒尿小便嘘嘘汇集6| 亚洲 国产 在线| 亚洲成人久久爱视频| 91字幕亚洲| 少妇人妻精品综合一区二区 | or卡值多少钱| 熟妇人妻久久中文字幕3abv| 男女之事视频高清在线观看| 国产熟女xx| 偷拍熟女少妇极品色| 欧美激情国产日韩精品一区| 69av精品久久久久久| 757午夜福利合集在线观看| 亚洲欧美日韩卡通动漫| 三级男女做爰猛烈吃奶摸视频| 在线观看午夜福利视频| 日本黄色视频三级网站网址| 成人av在线播放网站| 99久久久亚洲精品蜜臀av| 能在线免费观看的黄片| 亚洲真实伦在线观看| 国产麻豆成人av免费视频| 久久欧美精品欧美久久欧美| 亚洲成a人片在线一区二区| 精品久久久久久,| 极品教师在线视频| 国产高清激情床上av| 国内久久婷婷六月综合欲色啪| 亚洲综合色惰| 日韩大尺度精品在线看网址| 国产一区二区激情短视频| 午夜免费男女啪啪视频观看 | 99热只有精品国产| 亚洲三级黄色毛片| 久久久精品欧美日韩精品| 久久久久免费精品人妻一区二区| 亚洲av熟女| 国产精品98久久久久久宅男小说| 久久性视频一级片| ponron亚洲| 国产美女午夜福利| 精品人妻视频免费看| 如何舔出高潮| 国产精品久久久久久精品电影| 精品人妻一区二区三区麻豆 | 国内精品美女久久久久久| 男人和女人高潮做爰伦理| 精品人妻偷拍中文字幕| 99热这里只有是精品在线观看 | 天堂影院成人在线观看| 国产精品一区二区三区四区久久| 人人妻人人看人人澡| 偷拍熟女少妇极品色| 国产精品乱码一区二三区的特点| 国模一区二区三区四区视频| 在线天堂最新版资源| 最近在线观看免费完整版| 三级国产精品欧美在线观看| 亚洲欧美激情综合另类| 精华霜和精华液先用哪个| 91麻豆av在线| 久久婷婷人人爽人人干人人爱| 老司机福利观看| 免费av毛片视频| 久久国产精品人妻蜜桃| 床上黄色一级片| 18+在线观看网站| 99国产综合亚洲精品| 日韩精品中文字幕看吧| 久久性视频一级片| 此物有八面人人有两片| 老司机午夜福利在线观看视频| 搡老岳熟女国产| 高清在线国产一区| 国产黄a三级三级三级人| 国产av麻豆久久久久久久| 噜噜噜噜噜久久久久久91| 一本综合久久免费| 久久午夜福利片| 亚洲人成网站在线播放欧美日韩| 国产亚洲精品久久久com| 日本免费一区二区三区高清不卡| 日本黄色片子视频| 一区福利在线观看| 高潮久久久久久久久久久不卡| 丰满乱子伦码专区| 怎么达到女性高潮| 舔av片在线| 日韩精品中文字幕看吧| 国产成人影院久久av| 老熟妇仑乱视频hdxx| 日本撒尿小便嘘嘘汇集6| 五月伊人婷婷丁香| 又爽又黄a免费视频| 少妇的逼好多水| 中文字幕高清在线视频| 欧美激情在线99| 在线观看免费视频日本深夜| 午夜精品久久久久久毛片777| 日韩精品青青久久久久久| 久久久国产成人精品二区| 一区福利在线观看| 脱女人内裤的视频| 亚洲最大成人av| 少妇高潮的动态图| 成年人黄色毛片网站| 亚洲精品在线观看二区| 亚洲第一欧美日韩一区二区三区| 国产精品,欧美在线| 亚洲三级黄色毛片| 国产黄色小视频在线观看| 白带黄色成豆腐渣| 成人av在线播放网站| 亚洲av成人不卡在线观看播放网| 特大巨黑吊av在线直播| 亚洲av.av天堂| 天堂av国产一区二区熟女人妻| 亚洲av日韩精品久久久久久密| 亚洲欧美日韩东京热| 性色avwww在线观看| 91久久精品电影网| 两人在一起打扑克的视频| 少妇丰满av| 最近最新中文字幕大全电影3| 国产欧美日韩精品亚洲av| 国产伦在线观看视频一区| 日日干狠狠操夜夜爽| 成年女人看的毛片在线观看| 欧美日韩福利视频一区二区| av国产免费在线观看| 99久久成人亚洲精品观看| 亚洲国产欧美人成| 一本综合久久免费| 桃红色精品国产亚洲av| 久久久久亚洲av毛片大全| 欧美高清性xxxxhd video| 在线免费观看的www视频| 9191精品国产免费久久| 69av精品久久久久久| 国产免费一级a男人的天堂| 桃色一区二区三区在线观看| 91麻豆精品激情在线观看国产| 99久久精品一区二区三区| 性色av乱码一区二区三区2| 在线a可以看的网站| 成人性生交大片免费视频hd| 午夜福利18| 久久久久精品国产欧美久久久| 久久人人爽人人爽人人片va | 十八禁网站免费在线| 免费av观看视频| 亚洲精华国产精华精| 三级国产精品欧美在线观看| 久久久国产成人精品二区| 久久99热这里只有精品18| 波多野结衣巨乳人妻| 日本a在线网址| 身体一侧抽搐| 如何舔出高潮| 变态另类丝袜制服| 中文字幕高清在线视频| 日韩欧美在线二视频| 国产主播在线观看一区二区| 亚洲三级黄色毛片| 99国产精品一区二区蜜桃av| 我的老师免费观看完整版| 99久国产av精品| 久久精品影院6| 欧美日韩瑟瑟在线播放| 亚洲人成电影免费在线| 久久中文看片网| 亚洲午夜理论影院| 一区二区三区免费毛片| 国产av在哪里看| 亚洲天堂国产精品一区在线| 99国产综合亚洲精品| 日韩中文字幕欧美一区二区| 美女高潮的动态| 亚洲精品色激情综合| 两性午夜刺激爽爽歪歪视频在线观看| 好男人在线观看高清免费视频| www日本黄色视频网| 国产精品电影一区二区三区| 免费搜索国产男女视频| 日韩欧美国产一区二区入口| 午夜久久久久精精品| 深爱激情五月婷婷| 日本黄大片高清| 人妻久久中文字幕网| 国产三级黄色录像| 亚洲专区国产一区二区| 97超级碰碰碰精品色视频在线观看| 国产私拍福利视频在线观看| 欧美激情国产日韩精品一区| 麻豆国产av国片精品| av欧美777| 国产伦精品一区二区三区四那| 日本 欧美在线| 婷婷精品国产亚洲av在线| 老女人水多毛片| 国产伦人伦偷精品视频| 亚洲一区二区三区色噜噜| 亚洲第一欧美日韩一区二区三区| 床上黄色一级片| .国产精品久久| 亚洲成a人片在线一区二区| 欧美黑人欧美精品刺激| 亚洲va日本ⅴa欧美va伊人久久| 别揉我奶头~嗯~啊~动态视频| 日本 av在线| netflix在线观看网站| 毛片女人毛片| 男女下面进入的视频免费午夜| 久久久久久久久中文| 桃红色精品国产亚洲av| 3wmmmm亚洲av在线观看| 少妇的逼水好多| 男插女下体视频免费在线播放| 一级a爱片免费观看的视频| 99久久精品一区二区三区| 18美女黄网站色大片免费观看| 久久人人精品亚洲av| 中文字幕人妻熟人妻熟丝袜美| 国产伦人伦偷精品视频| 在线观看av片永久免费下载| 亚洲第一电影网av| 欧美日韩乱码在线| 欧美一区二区精品小视频在线| 国产av一区在线观看免费| 俄罗斯特黄特色一大片| 12—13女人毛片做爰片一| 色尼玛亚洲综合影院| 国产精品精品国产色婷婷| 精品一区二区三区人妻视频| 中文字幕免费在线视频6| 成人特级av手机在线观看| 亚洲 欧美 日韩 在线 免费| 蜜桃亚洲精品一区二区三区| a级毛片免费高清观看在线播放| 婷婷亚洲欧美| 久久99热6这里只有精品| 少妇的逼好多水| 变态另类成人亚洲欧美熟女| 久久精品国产99精品国产亚洲性色| 国产成年人精品一区二区| 国产白丝娇喘喷水9色精品| 亚洲欧美激情综合另类| 黄片小视频在线播放| 99久久九九国产精品国产免费| 亚洲人成伊人成综合网2020| 亚州av有码|