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

    掌葉大黃4個(gè)MYB轉(zhuǎn)錄因子的分子克隆及脅迫表達(dá)

    2024-07-12 13:39:22趙霞李元敏李依民胡曉晨高靜顏永剛張崗
    關(guān)鍵詞:脅迫實(shí)時(shí)熒光定量PCR表達(dá)分析

    趙霞 李元敏 李依民 胡曉晨 高靜 顏永剛 張崗

    doi:10.7606/j.issn.1004-1389.2024.07.017

    https://doi.org/10.7606/j.issn.1004-1389.2024.07.017

    收稿日期:2023-01-05? 修回日期:2023-03-07

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(82104334,81973430);咸陽市中青年科技創(chuàng)新領(lǐng)軍人才(L2022CXNLRC009);陜西中醫(yī)藥大學(xué)學(xué)科創(chuàng)新團(tuán)隊(duì)項(xiàng)目(2019-QN01)。

    第一作者:趙? 霞,女,碩士研究生,研究方向?yàn)榉肿由帉W(xué)。E-mail:2443284316@qq.com

    李元敏,并列第一作者,女,碩士研究生,研究方向?yàn)榉肿由帉W(xué)。E-mail:1272151264@qq.com

    通信作者:張? 崗,男,教授,研究方向?yàn)榉肿由帉W(xué)。E-mail:jay_gumling2003@aliyun.com

    李依民,女,副教授,研究方向?yàn)榉肿由帉W(xué)。E-mail:2051058@sntcm.edu.cn

    摘? 要? 旨在克隆獲得? RpMYB2、? RpMYB3、? RpMYB5、? RpMYB6基因ORF,并對(duì)其理化特性、系統(tǒng)進(jìn)化關(guān)系、不同組織表達(dá)、激素與非生物脅迫下的表達(dá)等展開分析。研究結(jié)果表明,? RpMYB2、? RpMYB3、? RpMYB5、? RpMYB6依次編碼333、275、293、291個(gè)氨基酸,分子質(zhì)量分別為34.781、30.525、31.243、33.313 ku,4個(gè)RpMYBs蛋白均為親水性蛋白并具有較高的熱穩(wěn)定性,除? RpMYB2外其余3個(gè)基因結(jié)構(gòu)主要為? α-螺旋和無規(guī)則卷曲,且均定位于細(xì)胞核。由蛋白結(jié)構(gòu)域和保守基序分析發(fā)現(xiàn),? RpMYB2、? RpMYB5含有1個(gè)HTH-MYB結(jié)構(gòu)域,? RpMYB3、? RpMYB6含有2個(gè)HTH-MYB結(jié)構(gòu)域,且不同MYB轉(zhuǎn)錄因子包含的保守元件數(shù)量及種類存在差異,系統(tǒng)進(jìn)化分析顯示? RpMYB2、? RpMYB5屬于R1-MYB亞家族,? RpMYB3、? RpMYB6屬于R2R3-MYB亞家族,與其結(jié)構(gòu)域分類一致。實(shí)時(shí)熒光定量結(jié)果顯示,根中? RpMYB2、? RpMYB6轉(zhuǎn)錄本豐度最高,?? RpMYB3、? RpMYB5基因分別在葉片、葉柄中高表達(dá)。經(jīng)200 μmol·L-1脫落酸 (abscisic acid, ABA)處理,? RpMYB5于24 h達(dá)峰值;200 μmol·L-1茉莉酸甲酯 (methyl jasmonate, MeJA)誘導(dǎo)?? RpMYB6 表達(dá),表達(dá)量隨時(shí)間推移呈逐漸上升趨勢(shì)且于24 h最高(為0 h的10.30倍),抑制? RpMYB3表達(dá);200 μmol·L-1水楊酸 (salicylic acid, SA)誘導(dǎo)? RpMYB6表達(dá)最為明顯,3 h表達(dá)量上調(diào)至? 0 h的21.84 倍。 RpMYBs 在非生物脅迫處理下表達(dá)也各有差異,? RpMYB2受高溫脅迫誘導(dǎo)表達(dá),? RpMYB3在損傷、高溫脅迫下表達(dá)下調(diào),? RpMYB5在損傷脅迫下表達(dá)受抑制,? RpMYB6響應(yīng)高溫、鹽脅迫。

    關(guān)鍵詞? 掌葉大黃;MYB;表達(dá)分析;脅迫;實(shí)時(shí)熒光定量PCR

    轉(zhuǎn)錄因子在植物生長(zhǎng)發(fā)育過程中發(fā)揮著不可或缺的作用,通常通過與靶基因啟動(dòng)子區(qū)域中的相關(guān)順式作用元件相互作用,參與基因表達(dá)調(diào)控[1]。在眾多的轉(zhuǎn)錄因子中,MYB存在于所有真核生物中,功能多樣且廣泛參與到植物生長(zhǎng)、發(fā)育、防御及次生代謝調(diào)控等生理過程。 ZmMYBC1是植物中第一個(gè)被發(fā)現(xiàn)的MYB轉(zhuǎn)錄因子[2],與玉米花青素的生物合成密切相關(guān)。近年來,越來越多的植物MYB家族成員相繼被報(bào)道,如擬南芥[3]、過路黃[4]、辣椒[5]中分別鑒定出197、51、172個(gè)MYB轉(zhuǎn)錄因子。

    通常,MYB蛋白在N端有一段高度保守的DNA結(jié)構(gòu)域,該結(jié)構(gòu)域由50~53個(gè)氨基酸殘基的重復(fù)片段(R)組成[6-7],依據(jù)包含的R結(jié)構(gòu)數(shù)目又分為4個(gè)亞類:4R-MYB,最小的MYB亞家族,相關(guān)報(bào)道較少;3R-MYB(R1R2R3-MYB),植物中家族成員較少,主要參與調(diào)控細(xì)胞周期[8-9];R2R3-MYB,植物MYB基因編碼蛋白中最豐富的一類,在初生與次生代謝、激素應(yīng)答以及植物對(duì)逆境脅迫的響應(yīng)等方面具有重要功能[10];R1-MYB是植物中第二大類MYB亞家族,在細(xì)胞形態(tài)建成、染色體結(jié)構(gòu)維持、根毛形成、生物鐘調(diào)節(jié)等生命過程具有重要調(diào)控作用[11]。諸多研究表明,MYB轉(zhuǎn)錄因子參與調(diào)控苯丙烷類次生代謝,如? McMYB4通過調(diào)節(jié)苯丙烷代謝和激素信號(hào)傳導(dǎo),參與蘋果響應(yīng)溫度變化的非生物抗性和生長(zhǎng)[12];白楊中? PtoMYB156參與調(diào)控苯丙烷途徑和次生細(xì)胞壁的初步形成[13];? AtMYB4通過苯丙烷代謝途徑參與調(diào)節(jié)擬南芥中花粉壁的形成[14]。

    掌葉大黃(Rheum palmatum? L.)是目前中國(guó)最為正宗、名貴的中藥材大黃屬三基源植物之一,以其干燥根、根莖為藥用,味苦、性寒,主治瀉下,清熱瀉火,涼血解毒,利濕祛黃[15],主要分布于甘肅、青海等高海拔地區(qū)。蒽醌類是大黃的主要活性物質(zhì),具有保肝、抗腫瘤、抗氧化等作用,其成分積累是影響大黃品質(zhì)的重要因素[16]。研究發(fā)現(xiàn)部分MYB轉(zhuǎn)錄因子負(fù)調(diào)控苯丙烷代謝[17],Perassolo等[18]闡明有效抑制莽草酸途徑中苯丙氨酸解氨酶活性可改變代謝流,活化蒽醌類成分合成途徑進(jìn)而導(dǎo)致蒽醌類成分累積。但是,MYB基因家族在掌葉大黃中是否參與調(diào)控蒽醌類化合物的合成還有待研究。李元敏[19]對(duì)1 a生掌葉大黃根、根莖和葉片進(jìn)行比較轉(zhuǎn)錄組測(cè)序分析(SRR10855670),篩選到6條組織差異性表達(dá)且包含完整開放閱讀框(ORF)的MYB轉(zhuǎn)錄因子序列。本研究通過RT-PCR克隆得到4個(gè)MYB家族成員基因,并對(duì)其進(jìn)行生物信息學(xué)、激素和非生物脅迫響應(yīng)分析,有助于深入探討MYB基因在植物防御反應(yīng)中的重要作用,為揭示大黃中蒽醌類活性成分合成與積累的調(diào)控機(jī)制及通過遺傳工程改良藥材品質(zhì)奠定基礎(chǔ)。

    1? 材料與方法

    1.1? 材? 料

    植物材料為1 a生掌葉大黃 (R. palmatum L.)及其成熟種子,2019年10月采集于甘肅省隴南市宕昌縣阿塢鄉(xiāng)麻界村,陜西中醫(yī)藥大學(xué)胡本祥教授對(duì)其進(jìn)行鑒定。取大黃根、根莖、葉、葉柄及種子等樣品并置于液氮中速凍后,于-80 ℃冰箱保存。

    基于黑小斌等[20]建立的大黃無菌苗體系培育出掌葉大黃無菌苗,參考李元敏等[21]的方法進(jìn)行植物激素與非生物脅迫處理。選取生長(zhǎng)旺盛、生長(zhǎng)均勻的1月齡幼苗,激素處理組:噴施200 μmol·L-1茉莉酸甲酯 (MeJA)、200 μmol·L-1水楊酸 (SA)、200 μmol·L-1脫落酸 (ABA);非生物脅迫處理組:包括干旱(10% PEG6000)、機(jī)械損傷(針刺葉片)、鹽(100?? mmol·L-1 NaCl)、低溫(4? ℃)、高溫(40? ℃);此外,噴施純?nèi)軇┠M激素對(duì)照組,噴灑無菌水作為脅迫對(duì)照組;所有樣品均于處理后1、3、6、12和24? h進(jìn)行取樣,0 h為空白對(duì)照,每個(gè)處理組的每個(gè)處理時(shí)間作3個(gè)樣本重復(fù),用液氮速凍后于-80? ℃冰箱保存。

    1.2? 總RNA提取及cDNA合成

    采用RN38-EASYspin Plus RNA試劑盒(艾德萊,中國(guó)北京)提取大黃各樣品總RNA,K5800自動(dòng)檢測(cè)超微量分光光度計(jì)(凱奧,中國(guó)北京)測(cè)定質(zhì)量,并用1.0%瓊脂糖凝膠電泳檢測(cè)RNA完整性。按說明書操作,利用PrimeScriptTM?? RT Master Mix反轉(zhuǎn)錄試劑盒(TaKaRa,日本)合成cDNA,-20 ℃保存。

    1.3? 基因克隆

    利用Primer Premier 5.0對(duì)篩選到的4條MYB轉(zhuǎn)錄因子序列設(shè)計(jì)跨ORF引物 (表1)。25 μL標(biāo)準(zhǔn)體系為:12.5 μL 2×Rapid Taq Master Mix、1 μL Forward Primer、1 μL Reverse Primer和1 μL cDNA,并用9.5 μL ddH2O補(bǔ)齊。PCR條件為:95? ℃ 5 min;95? ℃ 30 s,60? ℃?? 30 s,72? ℃ 2 min,40個(gè)循環(huán);72 ℃ 5 min,15? ℃保溫。用1.0%瓊脂糖凝膠電泳檢測(cè)擴(kuò)增產(chǎn)物。隨后用50 μL體系進(jìn)行擴(kuò)增、DNA膠回收試劑盒(TianGen,中國(guó)北京)純化回收目的片段,與pUC19連接后并轉(zhuǎn)化大腸桿菌感受態(tài)細(xì)胞DH5α,挑取陽性克隆送上海生工完成測(cè)序。

    1.4? 生物信息學(xué)分析

    利用軟件及在線工具對(duì)掌葉大黃 RpMYBs 進(jìn)行生物信息分析預(yù)測(cè)(表2)。

    1.5? 基因表達(dá)模式

    利用實(shí)時(shí)熒光定量PCR檢測(cè) RpMYBs 在大黃不同樣品中基因表達(dá)。以β-actin為內(nèi)參基因[22],qPCR引物序列見表1。反應(yīng)體系20 μL:包括10 μL TB Green○R? Premix Ex TaqTM? Ⅱ (TliRNaseHPlus) (2×)、Forward/Reverse Primer各0.8 μL、2 μL cDNA、0.4 μL ROX Reference Dye (50×)、6 μL ddH2O。反應(yīng)程序?yàn)椋?5? ℃ 30 s;95? ℃ 5 s,60? ℃ 30 s,60? ℃ 34 s,40個(gè)循環(huán);之后繪制熔解曲線,條件為95? ℃ 15 s,? 60? ℃ 1 min,95? ℃ 15 s。設(shè)置3次技術(shù)和試驗(yàn)重復(fù),包括不加模板的對(duì)照。用Step One Plus○R Real-Time (Applied Biosystems,美國(guó))儀器進(jìn)行qPCR檢測(cè),數(shù)據(jù)分析采用2-△△Ct計(jì)算基因相對(duì)表達(dá)量[23],以SPSS? 24軟件進(jìn)行統(tǒng)計(jì)分析,P<0.05代表有顯著性差異。

    2? 結(jié)果與分析

    2.1? 基因克隆

    基于陜西中醫(yī)藥大學(xué)藥學(xué)院中藥資源與開發(fā)課題組前期篩選到的MYB基因家族成員,本研究對(duì)其中4條包含完整開放閱讀框(ORF)的MYB轉(zhuǎn)錄因子序列進(jìn)行基因克隆研究。提取掌葉大黃全株的總RNA,以反轉(zhuǎn)錄的cDNA為模板,采用高保真DNA聚合酶進(jìn)行RT-PCR擴(kuò)增,產(chǎn)物各取5 μL經(jīng)1.0%瓊脂糖凝膠電泳檢測(cè),分別獲得4條清晰且特異性較好的條帶 (圖1),轉(zhuǎn)錄組數(shù)據(jù)與PCR產(chǎn)物大小相符且與克隆測(cè)序ORF序列吻合,基因命名為? RpMYB2、? RpMYB3、? RpMYB5、? RpMYB6 (注冊(cè)號(hào) MW269444、MW269445、MW269447、MW269448) (表3)。

    2.2? 蛋白理化特性及結(jié)構(gòu)域分析

    利用ProtParam分析RpMYB2、RpMYB3、RpMYB5、RpMYB6蛋白序列的一級(jí)結(jié)構(gòu)發(fā)現(xiàn)(表4),蛋白的分子質(zhì)量在30.525 (RpMYB3)~34.781 (RpMYB2) ku,等電點(diǎn)為6.18~9.14,

    其中,RpMYB2、RpMYB3屬堿性蛋白,RpMYB5、RpMYB6屬酸性蛋白。4個(gè)MYB蛋白的平均親水性是負(fù)值,說明其均為親水性蛋白。二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果顯示,RpMYB2、RpMYB3、RpMYB5、RpMYB6均由α-螺旋(alphahelix, Hh)、延伸鏈 (extendedstrand,Ee)、β-轉(zhuǎn)角 (betaturn,Tt) 和無規(guī)卷曲 (random coil,Cc) 組成,其中,RpMYB3、RpMYB5、RpMYB6主要由α-螺旋和無規(guī)卷曲構(gòu)成,而RpMYB2主要由延伸鏈和無規(guī)卷曲構(gòu)成。亞細(xì)胞定位結(jié)果顯示其均定位于細(xì)胞核,SignalP-5.0和TMHMM在線分析表明RpMYB2屬于跨膜蛋白,其余3個(gè)均無跨膜結(jié)? 構(gòu)域。

    借助InterProScan與PROSITE對(duì)RpMYBs蛋白結(jié)構(gòu)域分析預(yù)測(cè)結(jié)果一致(圖2),顯示RpMYB2、RpMYB5含有1個(gè)HTH-MYB結(jié)構(gòu)域,屬于R1-MYB亞族,RpMYB3、RpMYB6是典型的R2R3-MYB,含有2個(gè)HTH-MYB結(jié)構(gòu)。

    2.3? 系統(tǒng)進(jìn)化分析

    從NCBI和PlantTFDB數(shù)據(jù)庫下載擬南芥、大豆、丹參、水稻、玉米、苦蕎麥、羊草7種植物MYB蛋白序列進(jìn)行多序列比對(duì),篩選出相似度較高的序列,采用MEGA7軟件中的鄰接法進(jìn)行MYB蛋白分子演化分析。圖3結(jié)果表明,RpMYB2與RpMYB5被劃入R1-MYB家族,且RpMYB2與GmMYB176、RpMYB5與OsMYBS1蛋白序列相似性最高;RpMYB3與RpMYB6屬于R2R3-MYB類型的家族成員,RpMYB3與FtMYB15、RpMYB6與擬南芥S20亞家族中AtMYB62/116同源性較高,可能具有相似的生物學(xué)功能。

    2.4? 保守基序分析

    利用在線軟件MEME對(duì)掌葉大黃4個(gè)MYB蛋白序列進(jìn)行保守基序分析(圖4)。結(jié)果表明,不同亞類的MYB基因包含的保守元件數(shù)量及種類存在差異。4個(gè)MYB序列中均包括Motif1,且存在3個(gè)均勻分布的氨基酸 (W) 殘基,可能是掌葉大黃MYB轉(zhuǎn)錄因子的特征基序,但每條序列中Motif所處位置不同,RpMYB3和RpMYB6中Motif1均位于近N端且與Motif2相鄰,而RpMYB2和RpMYB5中Motif1單獨(dú)位于序列中間區(qū)段內(nèi)。Motif3是R2R3亞族的特征基序,Motif4是R1亞族的特征基序,二者均含有6個(gè)氨基酸殘基,是最短基序。Motif2是最長(zhǎng)基序,包括48個(gè)氨基酸殘基。

    2.5? 組織表達(dá)特異性

    以根為參考樣本,qPCR檢測(cè)4個(gè) RpMYBs 基因在大黃根、根莖、葉片、葉柄和種子中的表達(dá)情況。從圖5可以看出,? RpMYB2在根中高表達(dá),葉柄中低表達(dá);? RpMYB3在種子中表達(dá)量最低,但在葉片中表達(dá)量為根中的4.39倍;相反的,? RpMYB5在葉片中的表達(dá)量?jī)H為根中的0.36倍,而在葉柄中表達(dá)量高;? RpMYB6在根中表達(dá)量最高,在其余器官中均較低。

    2.6? 激素應(yīng)答分析

    分析MYB基因于MeJA、SA、ABA激素處理下不同時(shí)間點(diǎn) (1、3、6、12、24 h) 的應(yīng)答情況 (圖6),以0 h(CK)為對(duì)照:? RpMYB2受SA誘導(dǎo),相對(duì)表達(dá)量隨時(shí)間推移呈現(xiàn)先下降后上升的變化趨勢(shì),同Mock變化趨勢(shì)一致,且于24 h最高,為CK的3.57倍;MeJA處理1 h、3 h顯著上調(diào)基因表達(dá);ABA處理1 h顯著上調(diào),12 h、24 h受到顯著抑制。

    RpMYB3受MeJA抑制表達(dá),相對(duì)表達(dá)量隨時(shí)間推移呈現(xiàn)先下降后上升的變化趨勢(shì),其中1、3、6 h與Mock相比表達(dá)量無差異,但于12 h最低,為CK的0.08倍;同樣,SA處理12 h顯著抑制基因表達(dá);ABA處理24 h顯著上調(diào),6 h下調(diào)表達(dá)。

    RpMYB5受SA誘導(dǎo),相對(duì)表達(dá)量呈先上升后下降再上升的變化趨勢(shì),24 h最高且與Mock相比無明顯變化;ABA處理也于24 h顯著上調(diào)基因表達(dá);MeJA處理1 h顯著誘導(dǎo)。

    RpMYB6受MeJA誘導(dǎo),相對(duì)表達(dá)量隨時(shí)間推移呈逐漸上升趨勢(shì),且于24 h最高,為CK的10.30倍;SA處理3 h顯著誘導(dǎo)表達(dá)且為CK的21.84倍,ABA處理6 h、12 h、24 h誘導(dǎo)表達(dá),1 h抑制表達(dá)。

    2.7? 非生物脅迫響應(yīng)

    分析掌葉大黃無菌苗于不同脅迫處理下 RpMYBs 基因在不同時(shí)間點(diǎn) (1、3、6、12、24 h) 的響應(yīng)情況 (圖7),以0 h為對(duì)照:

    RpMYB2受高溫脅迫誘導(dǎo),相對(duì)表達(dá)量于24 h達(dá)峰值,且為CK的5.19倍;損傷脅迫3 h顯著誘導(dǎo),為CK的2.44倍,1 h顯著抑制;干旱、低溫、鹽脅迫表達(dá)水平未見明顯變化。

    RpMYB3受損傷脅迫抑制表達(dá),相對(duì)表達(dá)量于3 h達(dá)谷值,為CK的0.30倍;干旱和低溫處理均在24 h受到顯著抑制;相反,鹽和高溫處理均在24 h顯著誘導(dǎo)。

    RpMYB5受損傷脅迫抑制表達(dá),相對(duì)表達(dá)量于24 h達(dá)谷值,為CK的0.14倍;干旱和高溫脅迫均在1 h、24 h誘導(dǎo)表達(dá);鹽脅迫于12 h顯著抑制;低溫脅迫于3 h、6 h抑制表達(dá)。

    RpMYB6受高溫脅迫誘導(dǎo),相對(duì)表達(dá)量隨時(shí)間推移呈先上升后下降再上升的變化趨勢(shì),且于24 h達(dá)峰值,為CK的4.68倍;受鹽、損傷脅迫抑制,均在24 h達(dá)谷值,且都為CK的0.09倍;干旱和低溫脅迫均在3 h上調(diào)表達(dá)。

    3? 討? 論

    MYB蛋白是植物生長(zhǎng)發(fā)育、代謝、防御和應(yīng)激反應(yīng)調(diào)控網(wǎng)絡(luò)中的關(guān)鍵因素[24]。鑒于此,對(duì)該類基因的研究不斷從模式植物向主要農(nóng)作物、藥用植物、園藝植物等方向擴(kuò)展,包括擬南芥[7]、小麥[25]、丹參[26]、棗樹[27]等。本研究利用RT-PCR技術(shù)從掌葉大黃中克隆到? RpMYB2、? RpMYB3、? RpMYB5和? RpMYB6的ORF序列,其中? RpMYB2編碼跨膜蛋白,可能介導(dǎo)細(xì)胞與外界之間的信號(hào)轉(zhuǎn)導(dǎo),并參與細(xì)胞膜內(nèi)外物質(zhì)、能量及信號(hào)的交換[28],其余3個(gè)屬于非跨膜蛋白。4個(gè)目標(biāo)序列均具有熱穩(wěn)定性和親水性,對(duì)其進(jìn)行亞細(xì)胞定位預(yù)測(cè)均為核蛋白,推測(cè)基因主要在細(xì)胞核內(nèi)行使功能。

    進(jìn)化樹分析主要用于基因家族成員分類及進(jìn)化關(guān)系研究,有助于推測(cè)親緣關(guān)系密切相關(guān)的蛋白家族成員的生物學(xué)功能。 RpMYBs 與其他7種植物MYBs的分子進(jìn)化分析表明,這些MYBs可分為兩個(gè)亞族,即R1-MYB和R2R3-MYB。其中? RpMYB2和? RpMYB5具有同源性,歸屬于R1-MYB亞族。? RpMYB2與? GmMYB176關(guān)系最近,R1-MYB轉(zhuǎn)錄因子? GmMYB176可調(diào)節(jié)大豆中異黃酮積累增加的關(guān)鍵基因? CHS8基因表達(dá),從而促進(jìn)異黃酮生物合成[29],推測(cè)? RpMYB2可能作用于苯丙烷類生物代謝途徑。同時(shí)? RpMYB2與? LcMYB2聚在一大支,猜測(cè)其可能具有干旱和ABA處理下促進(jìn)種子萌發(fā)和根系生長(zhǎng)的功能[30]。?? RpMYB5與? OsMYBS1親緣關(guān)系密切,而? OsMYBS1可抑制葡萄糖傳導(dǎo)因子HXK1的表達(dá),使ABA活力下降,繼而促進(jìn)種子發(fā)芽及幼苗的早期發(fā)育[31];類似地,? RpMYB3和? RpMYB6具有同源性,并且聚集在R2R3-MYB亞家族。?? RpMYB3與? FtMYB15親緣關(guān)系最近,推測(cè)其可能一定程度上參與植物花青素合成的調(diào)控[32]。? RpMYB6與擬南芥S20亞家族成員? AtMYB62和? AtMYB116聚在一支,可能通過赤霉素信號(hào)通路調(diào)節(jié)磷脅迫反應(yīng)[33]。然而,4個(gè)RpMYBs是否參與調(diào)控大黃蒽醌類化合物合成仍有待深入研究。

    基因結(jié)構(gòu)是決定基因功能的前提。分析RpMYBs蛋白結(jié)構(gòu)域和保守基序發(fā)現(xiàn),同一亞家族含有相似的結(jié)構(gòu)域和保守基序。R1-MYB類蛋白中所含Motif數(shù)量與分布不同,表明R1-MYB類轉(zhuǎn)錄因子的空間位置保守性相對(duì)較低。而R2R3-MYB類蛋白中Motif類型和分布情況均較為相似,表明其在R2R3-MYB中的分布不具有隨機(jī)性,但具有較高的保守性,推測(cè)其在植物體內(nèi)發(fā)揮核心調(diào)控功能。MYB轉(zhuǎn)錄因子家族雖在序列特征上相對(duì)保守,但在不同物種不同器官中功能呈現(xiàn)多樣化。擬南芥中? AtMYB33、? AtMYB65和? AtMYB101可促進(jìn)花藥和花粉發(fā)育[34]。? MoMYB53于葉中高表達(dá),參與調(diào)控巴戟天子葉中黃酮醇的合成[35]。本研究采用SYBRGreen實(shí)時(shí)熒光定量PCR方法,測(cè)定 RpMYBs基因在1 a生掌葉大黃不同組織中的相對(duì)表達(dá)量,揭示? RpMYB2、? RpMYB6在根中高表達(dá)特性,推測(cè)可能通過調(diào)控大黃藥用部位根的生長(zhǎng)發(fā)育及次生代謝進(jìn)而影響大黃藥材品質(zhì);? RpMYB3、? RpMYB5分別在葉片、葉柄中高表達(dá),可能在地上部分的生長(zhǎng)發(fā)育中發(fā)揮著重要的調(diào)控作用。葉片直接暴露于大氣中,受到陽光和溫差變化的影響大,根深埋地下,受到土壤微生態(tài)顯著影響,這勢(shì)必導(dǎo)致MYB基因在表達(dá)水平做出響應(yīng),從而可能在植物適應(yīng)環(huán)境和介導(dǎo)次生代謝產(chǎn)物合成方面發(fā)揮重要作用。

    MYB類轉(zhuǎn)錄因子廣泛參與植物對(duì)非生物脅迫的響應(yīng),如調(diào)控生長(zhǎng)發(fā)育,影響代謝產(chǎn)物的合成和影響激素信號(hào)等[36]。擬南芥? AtMYB41通過促進(jìn)ABA積累參與植物脅迫反應(yīng)[37]。? CaMYB306通過調(diào)控MDA和H2O2含量、ROS清除酶活性及抗性相關(guān)基因的表達(dá),負(fù)調(diào)控辣椒對(duì)低溫的響應(yīng)[38]。? MYB96既能通過CBFs途徑參與低溫反應(yīng),又能通過ABA的信號(hào)通路調(diào)節(jié)其抗旱性[39]。本研究結(jié)果表明4個(gè)大黃 RpMYBs響應(yīng)激素和非生物脅迫ABA處理導(dǎo)致? RpMYB2表達(dá)下調(diào),而高溫脅迫下24 h時(shí)顯著誘導(dǎo)其表達(dá),說明? RpMYB2可能通過ABA信號(hào)通路參與高溫脅迫應(yīng)答。MeJA處理? RpMYB6顯著上調(diào)表達(dá),且受低溫、高溫脅迫誘導(dǎo)表達(dá)。MeJA是植物CBF/DREB1冷脅迫響應(yīng)途徑的早期關(guān)鍵信號(hào)[24],而 RpMYB4[20]在MeJA和低溫脅迫處理下均下調(diào)表達(dá),可能參與MeJA負(fù)調(diào)控凍害脅迫的耐受過程。? RpMYB2和? RpMYB4可能響應(yīng)冷脅迫轉(zhuǎn)錄調(diào)控網(wǎng)絡(luò),更好地維持大黃生存,同時(shí)也為激素和冷信號(hào)通路互相協(xié)調(diào)傳遞和響應(yīng)寒冷信號(hào)的作用機(jī)制提供了新的理論基礎(chǔ)。? RpMYB3受MeJA應(yīng)答和高溫脅迫響應(yīng)表達(dá)趨勢(shì)一致,均在12 h達(dá)谷值后緩慢恢復(fù),可能通過MeJA信號(hào)傳遞正調(diào)控高溫脅迫應(yīng)答。? RpMYB5可能通過ABA和SA信號(hào)途徑共同負(fù)調(diào)控大黃的機(jī)械損傷脅迫防御。 RpMYBs 基因響應(yīng)激素和非生物脅迫差異表達(dá),表明它們可能通過參與多種生物學(xué)過程的調(diào)節(jié)影響掌葉大黃的生理? 代謝。

    4? 結(jié)? 論

    基于掌葉大黃轉(zhuǎn)錄組數(shù)據(jù),對(duì)掌葉大黃MYB轉(zhuǎn)錄因子進(jìn)行基因克隆、生物信息學(xué)研究,并合理預(yù)測(cè)了其在掌葉大黃中的部分功能,為進(jìn)一步揭示 RpMYBs 基因功能探究其生理機(jī)制奠定基礎(chǔ),也為掌葉大黃蒽醌類成分合成調(diào)控及遺傳改良提供科學(xué)支撐。

    參考文獻(xiàn)? Reference:

    [1]? 張椿雨,龍? 艷,馮? 吉,等.植物基因在轉(zhuǎn)錄水平上的調(diào)控及其生物學(xué)意義 [J].遺? 傳,2007(7):793-799.

    ZHANG CH Y,LONG Y,F(xiàn)ENG J,et al.Transcriptional regulation of plant genes and its signifi-cance biology [J].Proceedings of the National Academy of SciencesProceedings of the National Academy of Sciences,2007(7):793-799.

    [2]? CONE K C,BURR F A,BURR B.Molecular analysis of the maize anthocyanin regulatory locus C1 [J]. Proceedings of the National Academy of Sciences USA,1986,83:9631-5.

    [3]? ABDULLAH-ZAWAWI M R,AHMAD-NIZAMMUDDIN N F,Govender N,et al.Comparative genome-wide analysis of WRKY,MADS-box and MYB transcription factor families in Arabidopsis and rice [J]. Scientific Reports,2021,11:19678.

    [4]? 楊敬天,黎? 華,趙? 歡,等.過路黃MYB轉(zhuǎn)錄因子家族成員的挖掘及鑒定 [J].中草藥,2022,53(16):5149-5158.

    YANG J T,LI H,ZHAO H,et al.Excavation and identification of members of MYB transcription factor family in Lysimachia christinae [J]. Chinese Traditional and Herbal Drugs, 2022,53(16):5149-5158.

    [5]? 居利香,雷? 欣,趙成志,等.辣椒MYB基因家族的鑒定及與辣味關(guān)系分析 [J].園藝學(xué)報(bào),2020,47(5):875-892.

    JU L X,LEI X,ZHAO CH ZH,et al.Identification of MYB family genes and its relationship with pungency of? pepper[J]. Acta Horticulturae Sinica ,2022,47(5):875-892.

    [6]? OGATA K,KANEI-ISHII C,SASAKI M,et al.The cavity in the hydrophobic core of Myb DNA-binding domain is reserved for DNA recognition and trans-activation [J]. Nature Structural and Molecular Biology,1996,3:178-187.

    [7]? DUBOS C,STRACKE R,GROTEWORD E,et al.MYB transcription factors in Arabidopsis[J].Trends in Plant Science,2010,15:573-581.

    [8]? ROSINSKI J A,ATCHLEY W R.Molecular evolution of the Myb family of transcription factors:evidence for polyphyletic origin [J]. Journal of Molecular Evolution,1998,46:74-83.

    [9]? 位欣欣,蘭海燕.植物MYB轉(zhuǎn)錄因子調(diào)控次生代謝及逆境響應(yīng)的研究進(jìn)展 [J].生物技術(shù)通報(bào),2022,38(8):12-23.

    WEI X X,LAN H Y.Advances in the regulation of plant MYB transcription factors in secondary metabolism and stress response [J]. Biotechnology Bulletin ,2022,38(8):12-23.

    [10]? STRACKE R,ISHIHARA H,HUEP G,et al.Differential regulation of closely related R2R3-MYB transcription factors controls flavonol accumulation in different parts of the Arabidopsis thaliana seedling [J]. Plant Journal, 2007,50:660-677.

    [11]? RAMALINGAM A,KUDAPA H,PAZHAMALA L T,et al.Gene expression and yeast two-hybrid studies of 1R-MYB transcription factor mediating drought stress response in chickpea(Cicer arietinum L.) [J]. Frontiers in Plant Science,2015,6:1117.

    [12]? HAO S X,LU Y F,PENG Z,et al.? McMYB4 improves temperature adaptation by regulating phenylpropanoid metabolism and hormone signaling in? apple[J]. Horticulture Research,2021,8:182.

    [13]? YANG L,ZHAO X,RAN L Y,et al.? PtoMYB156 is involved in negative regulation of phenylpropanoid metabolism and secondary cell wall biosynthesis during wood formation in? poplar[J]. Scientific Reports, 2017,7:41209.

    [14]? PRESTON J,WHEELER J,HEAZLEWOOD J,et al. AtMYB32 is required for normal pollen development in Arabidopsis thaliana [J]. Plant Journal,2004,40:979-995.

    [15]? 國(guó)家藥典委員會(huì). 中華人民共和國(guó)藥典:一部[M]. 北京:中國(guó)醫(yī)藥科技出版社, 2020.

    Chinese Pharmacopoeia Commission.Pharmacopoeia of the Peoples Republic of China[M]. Beijing:China Pharmaceutical Science and Technology Press,2020.

    [16]? 陳艷琰,唐于平,陳嘉倩,等.大黃資源化學(xué)研究進(jìn)展與利用策略 [J].中草藥,2018,49(21):5170-5178.

    CHEN Y Y,TANG Y P,CHEN J Q,et al.Research progress and utilization strategy on resource chemistry of Rhei Radix et Rhizoma[J]. Chinese Traditional and Herbal Drugs, 2018,49:5170-5178.

    [17]? DONG N Q,LIN H X,Contribution of phenylpropanoid metabolism to plant development and plant-environment interactions [J]. Journal of Integrative Plant Biology,2021,63:180-209.

    [18]? PERASSOLO M,QUEVEDO C,BUSTO V,et al.Enhance of anthraquinone production by effect of praline and aminoindan-2-phosphonic acid in Rubia linctorum suspension cultures [J]. Enzyme and Microbial Technology,2007,41:181-185.

    [19]? 李元敏.掌葉大黃MYB轉(zhuǎn)錄因子基因克隆與特征研究 [D].陜西咸陽:陜西中醫(yī)藥大學(xué),2021.

    LI Y M.Cloning and characterization of MYB transcription factor gene in? Rheum palmatum? L [D].Xianyang Shaanxi:Shaanxi University of Chinese Medicine,2021.

    [20]? 黑小斌,李依民,李? 歡,等.掌葉大黃種子特性研究及無菌培養(yǎng)體系的構(gòu)建[J].中草藥,2019,50(18):4430-4437.

    HEI X B,LI Y M,LI H,et al.Investigation on seeds characteristics and establishment of aseptic cultural system for? Rheum palmatum? [J]. Chinese Traditional and Herbal Drugs ,2019,50(18):4430-4437.

    [21]? 李元敏,張? 娜,李? 慧,等.掌葉大黃轉(zhuǎn)錄因子? RpMYB4基因克隆及分子表達(dá)特性 [J].藥學(xué)學(xué)報(bào),2021,56(4):1170-1177.

    LI Y M,ZHANG N,LI H,et al.Cloning and expression characteristics of a transcription factor gene? RpMYB4? in? Rheum palmatum? L.[J]. Acta Pharmaceutica Sinica, 2021,56(4):1170-1177.

    [22]? 王寧杭,陸丹迎,常鵬杰,等.望春玉蘭實(shí)時(shí)定量PCR分析中內(nèi)參基因的篩選 [J].分子植物育種,2019,17(11):3674-3680.

    WANG N H,LU D Y,CHANG P J,et al.Selection of internal reference genes in real-time quantitative PCR analysis of Magnolia biondii[J]. Molecular Plant Breeding,2019,17(11):3674-3680.

    [23]? LIVAK K J,SCHMITTGEN T D.Analysis of relative gene expression data using realtime quantitative PCR and the 2-CT method [J]. Methods,2001,25:402-408.

    [24]? 邱文怡,王詩雨,李曉芳,等.MYB轉(zhuǎn)錄因子參與植物非生物脅迫響應(yīng)與植物激素應(yīng)答的研究進(jìn)展 [J].浙江農(nóng)業(yè)學(xué)報(bào),2020,32(7):1317-1328.

    QIU W Y,WANG SH Y,LI X F,et al.Functions of plant MYB transcription factors in response to abiotic stress and plant hormones [J]. Acta Agriculturae Zhejiangensis, 2020,32(7):1317-1328.

    [25]? 王延謙,張? 波,陳文杰,等.小麥R2R3-MYB轉(zhuǎn)錄因子 TaMYB3-4D的克隆及功能分析 [J].西北植物學(xué)報(bào),2015,35(4):646-652.

    WANG Y Q,ZHANG B,CHEN W J,et al.Cloning and functional verification of R2R3-MYB transcription factor? TaMYB3-4D? in wheat[J]. Acta Botanica Boreali-Occidentalia Sinica,2015,35(4):646-652.

    [26]? 彭晶晶,李? 凱,楊? 靖,等.丹參轉(zhuǎn)錄因子SmMYB52的基因克隆、表達(dá)分析和亞細(xì)胞定位 [J].基因組學(xué)與應(yīng)用生物學(xué),2021,40(2):802-808.

    PENG J J,LI K,YANG J,et al.Gene clone, expression?? analysis and subcellular localization of transcription factor?? SmMYB52 in Salvia miltiorrhiza [J]. Genomics and Applied Biology,2021,40(2):802-808.

    [27]? 王軍曉,李? 婷,馬志博,等.棗樹MYB基因家族鑒定、果實(shí)表達(dá)和鹽脅迫響應(yīng)研究 [J].西北林學(xué)院學(xué)報(bào),2022,? 37(2):20-28,36.

    WANG J X,LI T,MA ZH B,et al.Identification,fruitexpression and response to salt stress of MYB gene family in Ziziphus jujuba [J]. Journal of Northwest Forestry?? University, 2022,37(2):20-28,36.

    [28]? 宋江華,張立新.植物跨膜蛋白研究進(jìn)展 [J].生物學(xué)雜志,2009,26(6):62-64.

    SONG J H,ZHANG L X.Progress on the transmembrane protein in plants [J]. Journal of Biology, 2009,26(6):62-64.

    [29]? LI X Y,DHAUBHADEL S.14-3-3 proteins act as scaffolds for? GmMYB62 and? GmMYB176 and regulate their intracellular localization in soybean [J]. Plant Signaling and Behavior,2012,7:965-968.

    [30] ?ZHAO P C,HOU S L,GUO X F,et al.A MYB-related transcription factor from? sheepgrass, LcMYB2,promotes seed germination and root growth under drought stress [J]. BMC Plant Biology,2019,19:564.

    [31]? CHEN Y S,CHAO Y C,TSENG T W,et al.Two MYB-related transcription factors play opposite roles in sugar signaling in Arabidopsis [J]. Plant Molecular Biology,2017,93:299-311.

    [32]? ZHANG K X,LOGACHEVA M D,MENG Y,et al.Jasmonate-responsive MYB factors spatially repress rutin biosynthesis in Fagopyrum tataricum [J]. Journal of Experimental Botany,2018,69:1955-1966.

    [33]? DEVAIAH B N,MADHUVANTHI R,KARTHIKEYAN A S,et al.Phosphate starvation responses and gibberellic acid biosynthesis are regulated by the? MYB62? transcription factor in Arabidopsis [J]. Molecular Plant,2009,2:43-58.

    [34]? COSTACHE A,BERGHI ON,CERGAN R,et al.Respiratory allergies:Salicaceae sensitization (Review) [J]. Experimental and Therapeutic Medicine,2021,21(6):609.

    [35]? 謝德金,葉友杰,楊德明,等.基于巴戟天轉(zhuǎn)錄組數(shù)據(jù)的R2R3-MYB轉(zhuǎn)錄因子的鑒定和分析 [J].藥學(xué)學(xué)報(bào),2020,55(1):160-167.

    XIE D J,YE Y J,YANG D M,et al.Identification and?? analysis of the R2R3-MYB transcription factor based on a transcriptome database in Morinda officinalis [J]. Acta Pharmaceutica Sinica ,2020,55(1):160-167.

    [36]? 杜恬恬,李會(huì)萍,王博雅,等.梁山慈竹? DfMYB3基因克隆及啟動(dòng)子分析 [J].植物研究,2021,41(5):729-737.

    DU T T,LI H P,WANG B Y,et al.Cloning and promotor analysis of?? DfMYB3? from Dendrocalamus farinosus [J]. Bulletin of Botanical Research, 2021,41(5):729-737.

    [37]? LIPPOLD F,SANCHEZ D H,MUSIALAK M,et al.? AtMYB41? regulates transcriptional and metabolic responses to osmotic stress in Arabidopsis [J]. Plant Physiology,2009,149:1761-72.

    [38]? MA X,YU Y N,JIA J H,et al.The pepper? MYB transcription factor? CaMYB306? accelerates fruit coloration and negatively regulates cold resistance [J]. Scientia Horticulturae,2022,295:110892.

    [39]? LEE H G,SEO P J. MYB96 recruits the HDA15 protein to suppress negative regulators of ABA signaling in Arabidopsis [J]. Nature Communications,2019,10:1713.

    Molecular Cloning and Stress Expression of Four MYB

    Transcription Factors in? Rheum palmatum

    ZHAO Xia1,2, LI Yuanmin1,3, LI Yimin1,2, HU Xiaochen1,2, GAO Jing1,2, YAN Yonggang1,2 and? ZHANG Gang1,2

    (1.College of Pharmacy and Shaanxi Qinling Application Development and Engineering Center of Chinese Herbal

    Medicine,Shaanxi University of Chinese Medicine, Xian? 712046, China; 2.Key Laboratory for Research and

    Development of? “Qin Medicine”?? of Shaanxi Administration of Traditional Chinese Medicine, Shaanxi

    University of Chinese Medicine, Xian? 712046,China; 3.The Second Affiliated Hospital

    of Shaanxi University of Chinese Medicine, Xian? 712046, China)

    Abstract? In this study,we searched the transcriptome database to obtain the sequence information of four MYB gene family members,?? RpMYB2,?? RpMYB3,?? RpMYB5 and?? RpMYB6 genes. The corresponding physicochemical properties, phylogenetic relationships, tissue-specific expression, and induced expression under hormonal and abiotic stresses were analyzed. The results showed that?? RpMYB2,?? RpMYB3,?? RpMYB5 and?? RpMYB6 encoded 333, 275, 293 and 291 amino acids with molecular masses of 34.781, 30.525, 31.243 and 33.313 ku, respectively, and the four RpMYBs proteins were all hydrophilic proteins with high thermal stability, except for?? RpMYB2. The remaining three gene structures were mainly α-helix and randomly coiled, and all of them were localized in the nucleus. Analysis of the protein structural domains and conserved motifs revealed that?? RpMYB2 and?? RpMYB5 contained one HTH-MYB structural domain, and?? RpMYB3 and?? RpMYB6 contained two HTH-MYB structural domains, and there were differences in the number and types of conserved elements contained in different MYB transcription factors. Phylogenetic analysis showed that?? RpMYB2 and?? RpMYB5 belonged to the R1-MYB subfamily, and?? RpMYB3 and?? RpMYB6 belonged to the R2R3-MYB subfamily, consistent with their structural domain classification. The protein structural domains and phylogenetic trees consistently indicated that?? RpMYB2 and?? RpMYB5 belonged to the R1-MYB subfamily, and?? RpMYB3 and?? RpMYB6 belonged to the R2R3-MYB subfamily. Real-time fluorescence quantification revealed that?? RpMYB2 and? ?RpMYB6 were the most highly expressed in roots, and?? RpMYB3 and?? RpMYB5 were highly expressed in leaves and petioles. Respectively,?? RpMYB5 peaked at 24 h when treated with 200 μmol·L-1 abscisic acid (ABA). The highest expression of?? RpMYB6 was induced by 200 μmol·L-1 methyl jasmonate (MeJA), and the expression level increased gradually with time and was 10.30 times higher at 24 h than at 0 h, and inhibited?? RpMYB3 expression. 200 μmol·L-1 salicylic acid (SA) induced the most pronounced expression of ??RpMYB6 at 21.84 fold. The expression of? RpMYBs? also varied under abiotic stress treatments.?? RpMYB2 was induced by high temperature.?? RpMYB3 was inhibited by damage and high temperature.?? RpMYB5 was inhibited by damage stress.?? RpMYB6 responded to high temperature and salt stress response. This study provides a foundation for the functional identification of RpMYBs , which will be useful for the development of germplasm innovation and molecular breeding of R. palmatum L.

    Key words? ?Rheum palmatum; MYB; Expression analysis; Stress; qRT-PCR

    Received ??2023-01-05??? Returned? 2023-03-07

    Foundation item? National Natural Science Foundation of China(No.82104334,No.81973430); Middle Youth Science and Technology Innovation Leaders of Xianyang Municipality (No.L2022CXNLRC009);Project for Discipline Innovation Team of Shaanxi University of Traditional Chinese Medicine(No.2019-QN01).

    First author? ZHAO Xia, female,master student.Research area:molecular biopharmacology.?? E-mail:2443284316@qq.com

    LI Yuanmin, female,master student.Research area:molecular biopharmacology.E-mail:1272151264@qq.com

    Corresponding?? author? ZHANG Gang,male,professor.Research area:molecular biopharmacology.? E-mail:jay_gumling2003@aliyun.com

    LI? Yimin,female,associate professor.Research area:molecular biopharmacology.E-mail:2051058@sntcm.edu.cn

    (責(zé)任編輯:郭柏壽? Responsible editor:GUO Baishou)

    猜你喜歡
    脅迫實(shí)時(shí)熒光定量PCR表達(dá)分析
    丹參蛋白激酶SmSnRK2.4的克隆及表達(dá)分析
    雞卵清蛋白啟動(dòng)子慢病毒載體構(gòu)建及熒光定量PCR檢測(cè)重組慢病毒滴度
    1—脫氧野尻霉素對(duì)家蠶中腸BmSuc1基因表達(dá)及其酶活性的影響
    鹽堿脅迫對(duì)知母種子萌發(fā)的影響
    αB—crystallin在喉癌組織中的表達(dá)及臨床意義
    紅花生育酚環(huán)化酶基因的克隆及表達(dá)分析
    民事脅迫構(gòu)成要件的法理分析
    商(2016年20期)2016-07-04 01:05:25
    玉米紋枯病病菌y—谷氨酰轉(zhuǎn)肽酶基因克隆與表達(dá)分析
    信號(hào)分子與葉銹菌誘導(dǎo)下小麥病程相關(guān)蛋白1基因的表達(dá)分析
    水稻OsCATA基因相互作用蛋白的篩選和表達(dá)分析
    h视频一区二区三区| 蜜桃在线观看..| 大陆偷拍与自拍| av播播在线观看一区| 又粗又硬又长又爽又黄的视频| 亚洲色图综合在线观看| 青春草国产在线视频| 十八禁高潮呻吟视频| 99国产精品免费福利视频| 欧美日韩一区二区视频在线观看视频在线| 亚洲美女视频黄频| 亚洲高清免费不卡视频| 自线自在国产av| 91精品三级在线观看| 国产永久视频网站| 中文精品一卡2卡3卡4更新| 国产免费一区二区三区四区乱码| 亚洲精品中文字幕在线视频| 黑人巨大精品欧美一区二区蜜桃 | 在线观看三级黄色| 久久国产精品大桥未久av| 99热这里只有精品一区| 丰满少妇做爰视频| 亚洲av.av天堂| 中文天堂在线官网| 你懂的网址亚洲精品在线观看| 色5月婷婷丁香| 久久亚洲国产成人精品v| 伦理电影大哥的女人| 91精品国产九色| 成年人午夜在线观看视频| 亚洲婷婷狠狠爱综合网| 一区二区三区四区激情视频| 一本大道久久a久久精品| 国产成人aa在线观看| 韩国av在线不卡| 飞空精品影院首页| 亚洲国产成人一精品久久久| videos熟女内射| 欧美精品高潮呻吟av久久| 久热这里只有精品99| 欧美激情 高清一区二区三区| 人成视频在线观看免费观看| 日韩三级伦理在线观看| 最后的刺客免费高清国语| 18禁观看日本| 人人妻人人澡人人爽人人夜夜| av播播在线观看一区| 搡老乐熟女国产| 人妻少妇偷人精品九色| 80岁老熟妇乱子伦牲交| 人妻夜夜爽99麻豆av| 亚洲av男天堂| 免费日韩欧美在线观看| 男女无遮挡免费网站观看| 欧美日韩一区二区视频在线观看视频在线| 午夜激情久久久久久久| 亚洲人与动物交配视频| 日韩成人伦理影院| 亚洲熟女精品中文字幕| 91精品伊人久久大香线蕉| 18禁动态无遮挡网站| 日韩,欧美,国产一区二区三区| 欧美日韩成人在线一区二区| 天堂8中文在线网| 麻豆乱淫一区二区| 欧美激情国产日韩精品一区| 哪个播放器可以免费观看大片| 国产亚洲午夜精品一区二区久久| 韩国av在线不卡| 18+在线观看网站| 天天躁夜夜躁狠狠久久av| 大片电影免费在线观看免费| 在线免费观看不下载黄p国产| 日韩精品免费视频一区二区三区 | 寂寞人妻少妇视频99o| 色吧在线观看| 女性被躁到高潮视频| 精品久久久噜噜| 亚洲欧美精品自产自拍| 一级a做视频免费观看| 男人爽女人下面视频在线观看| 欧美一级a爱片免费观看看| 狠狠婷婷综合久久久久久88av| 九九久久精品国产亚洲av麻豆| 毛片一级片免费看久久久久| 国产日韩欧美亚洲二区| 成人毛片60女人毛片免费| 国产成人精品婷婷| 韩国av在线不卡| 亚洲欧美成人精品一区二区| 日本午夜av视频| 十分钟在线观看高清视频www| 亚洲精品久久午夜乱码| 亚洲欧美日韩另类电影网站| av在线app专区| 国产av码专区亚洲av| 美女国产视频在线观看| 中文乱码字字幕精品一区二区三区| 亚洲三级黄色毛片| 日韩av不卡免费在线播放| √禁漫天堂资源中文www| 狂野欧美激情性bbbbbb| 久久久久久久久久成人| 亚洲av成人精品一二三区| 精品少妇内射三级| 国产高清三级在线| 精品人妻偷拍中文字幕| 美女xxoo啪啪120秒动态图| 国产精品国产三级国产专区5o| 97超碰精品成人国产| 一区二区三区精品91| 一级二级三级毛片免费看| 热re99久久精品国产66热6| 麻豆成人av视频| 成人手机av| 黄色怎么调成土黄色| 婷婷成人精品国产| 99精国产麻豆久久婷婷| 国产精品三级大全| 国产一级毛片在线| 久久精品人人爽人人爽视色| 777米奇影视久久| 男人爽女人下面视频在线观看| 一区二区日韩欧美中文字幕 | 欧美人与性动交α欧美精品济南到 | 美女cb高潮喷水在线观看| tube8黄色片| 国产免费福利视频在线观看| 精品一区在线观看国产| 中文字幕亚洲精品专区| 天堂中文最新版在线下载| 大码成人一级视频| 久久久精品区二区三区| 亚洲成人av在线免费| 各种免费的搞黄视频| 女性被躁到高潮视频| 久久精品熟女亚洲av麻豆精品| 美女大奶头黄色视频| 搡老乐熟女国产| 国产精品久久久久久精品古装| 亚洲国产日韩一区二区| 国产精品国产三级国产av玫瑰| 涩涩av久久男人的天堂| 这个男人来自地球电影免费观看 | 亚洲一级一片aⅴ在线观看| 一区在线观看完整版| 精品一区二区免费观看| 国产精品女同一区二区软件| 99久久综合免费| 久久人人爽人人片av| 亚洲欧美清纯卡通| 久久久久网色| 国产爽快片一区二区三区| 精品人妻在线不人妻| 国产熟女欧美一区二区| 91久久精品电影网| av免费在线看不卡| 亚洲精品av麻豆狂野| 国产亚洲一区二区精品| 国产精品麻豆人妻色哟哟久久| a级片在线免费高清观看视频| 人体艺术视频欧美日本| 国产一区二区三区av在线| 亚洲内射少妇av| av国产久精品久网站免费入址| 久久人人爽人人片av| 九九爱精品视频在线观看| 久久精品人人爽人人爽视色| 午夜激情久久久久久久| 午夜精品国产一区二区电影| 人妻夜夜爽99麻豆av| 亚洲国产精品一区三区| 精品亚洲成a人片在线观看| 高清黄色对白视频在线免费看| 国产片内射在线| 亚洲国产欧美日韩在线播放| 大话2 男鬼变身卡| 黄色毛片三级朝国网站| 午夜福利视频在线观看免费| 国产午夜精品久久久久久一区二区三区| 一区二区三区四区激情视频| 成年av动漫网址| 你懂的网址亚洲精品在线观看| 夫妻午夜视频| 极品少妇高潮喷水抽搐| 亚洲精品,欧美精品| 亚洲成人手机| 蜜桃在线观看..| 亚洲av成人精品一二三区| 曰老女人黄片| 最近中文字幕2019免费版| 一级a做视频免费观看| 日韩一区二区三区影片| 日本黄大片高清| 欧美另类一区| 亚洲欧美色中文字幕在线| 蜜臀久久99精品久久宅男| 在线观看免费高清a一片| 精品久久久噜噜| 制服诱惑二区| av国产精品久久久久影院| 亚洲精品456在线播放app| 最近中文字幕2019免费版| 国产亚洲最大av| 成年av动漫网址| 精品人妻偷拍中文字幕| 午夜福利影视在线免费观看| 久久久久国产网址| 国产片内射在线| 乱人伦中国视频| 久久热精品热| 欧美日韩国产mv在线观看视频| 80岁老熟妇乱子伦牲交| 各种免费的搞黄视频| 性色av一级| 国产熟女欧美一区二区| 久久 成人 亚洲| 精品一品国产午夜福利视频| 精品一区二区免费观看| 飞空精品影院首页| 日韩av不卡免费在线播放| 大片免费播放器 马上看| 久久婷婷青草| 亚洲av男天堂| 久久久久久久国产电影| 亚洲精品,欧美精品| 中文字幕精品免费在线观看视频 | 国产永久视频网站| 成人18禁高潮啪啪吃奶动态图 | 久久久久网色| 91成人精品电影| 亚洲性久久影院| 久久精品国产亚洲av涩爱| 日韩一本色道免费dvd| 97超视频在线观看视频| 欧美日韩国产mv在线观看视频| 女的被弄到高潮叫床怎么办| 久久精品国产亚洲av涩爱| 亚洲天堂av无毛| 日本爱情动作片www.在线观看| 蜜桃国产av成人99| 91精品三级在线观看| 久久精品国产自在天天线| 免费观看性生交大片5| 91精品国产九色| 国产白丝娇喘喷水9色精品| av福利片在线| 日本爱情动作片www.在线观看| 在线观看www视频免费| 人体艺术视频欧美日本| 久久 成人 亚洲| 另类亚洲欧美激情| 在线观看三级黄色| 高清在线视频一区二区三区| 在线看a的网站| 2018国产大陆天天弄谢| 久久国产亚洲av麻豆专区| 欧美性感艳星| 国内精品宾馆在线| 韩国高清视频一区二区三区| 日韩电影二区| 三级国产精品欧美在线观看| 亚洲美女黄色视频免费看| av天堂久久9| 超碰97精品在线观看| 久久这里有精品视频免费| 亚洲av中文av极速乱| 99久久中文字幕三级久久日本| 一本大道久久a久久精品| 超碰97精品在线观看| 免费看光身美女| 如日韩欧美国产精品一区二区三区 | 国产午夜精品久久久久久一区二区三区| 卡戴珊不雅视频在线播放| 女的被弄到高潮叫床怎么办| 久久久久久久久久久免费av| 亚洲av电影在线观看一区二区三区| 亚洲五月色婷婷综合| 99视频精品全部免费 在线| 亚洲一级一片aⅴ在线观看| 亚洲av不卡在线观看| 天天躁夜夜躁狠狠久久av| 久久影院123| 亚洲欧美色中文字幕在线| 免费大片黄手机在线观看| 一边亲一边摸免费视频| 久久精品国产鲁丝片午夜精品| 国产免费福利视频在线观看| 99久久精品国产国产毛片| 性高湖久久久久久久久免费观看| √禁漫天堂资源中文www| 国产亚洲av片在线观看秒播厂| 成人亚洲精品一区在线观看| 亚洲精品国产色婷婷电影| 国产色爽女视频免费观看| 纵有疾风起免费观看全集完整版| 国语对白做爰xxxⅹ性视频网站| 国产午夜精品一二区理论片| 国产精品秋霞免费鲁丝片| 成人二区视频| 国产精品成人在线| 日日摸夜夜添夜夜爱| 国产色爽女视频免费观看| 80岁老熟妇乱子伦牲交| 亚洲第一区二区三区不卡| 免费久久久久久久精品成人欧美视频 | 汤姆久久久久久久影院中文字幕| 精品一区二区三区视频在线| 满18在线观看网站| 久久人人爽av亚洲精品天堂| kizo精华| 黑人猛操日本美女一级片| 亚洲人成网站在线播| 只有这里有精品99| 在线免费观看不下载黄p国产| 日韩熟女老妇一区二区性免费视频| 黑人欧美特级aaaaaa片| 特大巨黑吊av在线直播| 日本与韩国留学比较| 丝袜脚勾引网站| 国产精品一区二区三区四区免费观看| xxxhd国产人妻xxx| 伦精品一区二区三区| 我的女老师完整版在线观看| 日本91视频免费播放| 99热国产这里只有精品6| 老司机影院成人| 熟女电影av网| 欧美精品人与动牲交sv欧美| 国产有黄有色有爽视频| 亚洲精品av麻豆狂野| 国产成人精品一,二区| 在线天堂最新版资源| 亚洲国产av新网站| 最黄视频免费看| 日韩一区二区三区影片| 久久久午夜欧美精品| 99久国产av精品国产电影| 欧美精品高潮呻吟av久久| 欧美日韩成人在线一区二区| 久久精品久久久久久久性| 日韩在线高清观看一区二区三区| 久久影院123| 80岁老熟妇乱子伦牲交| 极品人妻少妇av视频| 一本大道久久a久久精品| 欧美老熟妇乱子伦牲交| 又大又黄又爽视频免费| 久久久久久久国产电影| 免费高清在线观看视频在线观看| 韩国高清视频一区二区三区| 亚洲精品aⅴ在线观看| 老司机亚洲免费影院| 日日摸夜夜添夜夜添av毛片| 亚洲国产精品一区三区| 国产成人精品一,二区| 18禁裸乳无遮挡动漫免费视频| 人妻系列 视频| 久久久午夜欧美精品| 亚洲婷婷狠狠爱综合网| 人体艺术视频欧美日本| 99九九线精品视频在线观看视频| 麻豆成人av视频| 久久久久久久久大av| 日日啪夜夜爽| 少妇被粗大猛烈的视频| av线在线观看网站| 免费高清在线观看视频在线观看| 少妇熟女欧美另类| 大话2 男鬼变身卡| 国产精品99久久99久久久不卡 | 久久精品久久久久久噜噜老黄| 国产黄色视频一区二区在线观看| 日韩av不卡免费在线播放| 精品久久久精品久久久| 久久久久久久久久成人| av网站免费在线观看视频| 国产熟女欧美一区二区| 婷婷色综合大香蕉| 亚洲婷婷狠狠爱综合网| 国产片特级美女逼逼视频| 国产av一区二区精品久久| 欧美xxⅹ黑人| 色哟哟·www| 熟女人妻精品中文字幕| 亚洲精品一区蜜桃| freevideosex欧美| 人人妻人人澡人人看| 日韩成人av中文字幕在线观看| 丝袜美足系列| 丝袜喷水一区| 亚洲成色77777| 亚洲国产av影院在线观看| 我要看黄色一级片免费的| 精品人妻在线不人妻| 日本免费在线观看一区| 美女xxoo啪啪120秒动态图| 如何舔出高潮| 夜夜爽夜夜爽视频| 黄色一级大片看看| 日韩伦理黄色片| 免费观看性生交大片5| 特大巨黑吊av在线直播| 两个人的视频大全免费| 国产一区二区在线观看av| 制服人妻中文乱码| av不卡在线播放| 久久久欧美国产精品| 高清黄色对白视频在线免费看| 国产成人免费观看mmmm| 欧美精品亚洲一区二区| 午夜福利视频在线观看免费| 日本黄大片高清| 国产成人精品久久久久久| 久久人人爽av亚洲精品天堂| 国产乱来视频区| 夜夜爽夜夜爽视频| 男人操女人黄网站| 国产成人a∨麻豆精品| 国产在线一区二区三区精| 超色免费av| 欧美另类一区| 秋霞在线观看毛片| 麻豆精品久久久久久蜜桃| 精品亚洲成国产av| 免费av中文字幕在线| 亚洲av电影在线观看一区二区三区| 欧美xxⅹ黑人| 国产精品国产三级国产av玫瑰| 中文字幕制服av| 99国产精品免费福利视频| 色婷婷av一区二区三区视频| 色哟哟·www| 中国国产av一级| 老司机影院成人| 免费观看无遮挡的男女| 久久精品国产亚洲av天美| 精品酒店卫生间| 亚洲经典国产精华液单| 国产免费现黄频在线看| 久久久久国产网址| 亚洲精华国产精华液的使用体验| 婷婷色麻豆天堂久久| 精品久久久久久久久av| av不卡在线播放| a级毛色黄片| 免费黄频网站在线观看国产| 汤姆久久久久久久影院中文字幕| 中文字幕最新亚洲高清| 乱码一卡2卡4卡精品| 亚洲精品自拍成人| 久久人人爽av亚洲精品天堂| 啦啦啦啦在线视频资源| 亚洲精品aⅴ在线观看| 18禁在线无遮挡免费观看视频| 国产精品久久久久久久久免| 成人亚洲精品一区在线观看| 极品人妻少妇av视频| 久久精品夜色国产| 亚洲人成77777在线视频| 亚洲av.av天堂| 久久精品国产亚洲网站| 欧美xxⅹ黑人| 9色porny在线观看| 国产乱人偷精品视频| 国产精品无大码| 一区二区三区免费毛片| av视频免费观看在线观看| 在线观看三级黄色| 在线观看www视频免费| 亚洲精品一二三| 2018国产大陆天天弄谢| 91精品国产九色| 夫妻性生交免费视频一级片| www.色视频.com| 亚洲精品成人av观看孕妇| 亚洲在久久综合| 伊人亚洲综合成人网| 成人18禁高潮啪啪吃奶动态图 | 中文乱码字字幕精品一区二区三区| 熟女av电影| 久久99蜜桃精品久久| 黄色怎么调成土黄色| 国产高清三级在线| 中文欧美无线码| 王馨瑶露胸无遮挡在线观看| 啦啦啦啦在线视频资源| 人妻夜夜爽99麻豆av| 亚洲精品av麻豆狂野| 美女中出高潮动态图| 男人操女人黄网站| 国产亚洲精品久久久com| 成人亚洲欧美一区二区av| 国产一区亚洲一区在线观看| 成人二区视频| 中文字幕人妻丝袜制服| 多毛熟女@视频| 亚洲国产精品国产精品| 大话2 男鬼变身卡| 丰满迷人的少妇在线观看| 久久99热6这里只有精品| 久久国产精品大桥未久av| 日韩av不卡免费在线播放| 日韩成人伦理影院| 极品人妻少妇av视频| 青春草国产在线视频| 免费久久久久久久精品成人欧美视频 | 精品国产一区二区三区久久久樱花| 26uuu在线亚洲综合色| 色婷婷久久久亚洲欧美| 黄色欧美视频在线观看| 亚洲欧美一区二区三区黑人 | 久久久久久久国产电影| av有码第一页| 亚洲内射少妇av| 国产免费又黄又爽又色| 成人黄色视频免费在线看| 丰满少妇做爰视频| 国产白丝娇喘喷水9色精品| 91久久精品国产一区二区三区| 精品99又大又爽又粗少妇毛片| 免费观看a级毛片全部| 日本欧美国产在线视频| 日韩在线高清观看一区二区三区| 看非洲黑人一级黄片| 久久精品久久久久久久性| 一级毛片黄色毛片免费观看视频| 在线精品无人区一区二区三| 婷婷色麻豆天堂久久| 亚洲欧洲精品一区二区精品久久久 | 精品久久蜜臀av无| 中文乱码字字幕精品一区二区三区| 国产又色又爽无遮挡免| 日韩中字成人| 插阴视频在线观看视频| 三上悠亚av全集在线观看| 街头女战士在线观看网站| 国精品久久久久久国模美| 中文字幕久久专区| 精品国产乱码久久久久久小说| 日韩亚洲欧美综合| 亚洲国产精品999| 看非洲黑人一级黄片| 在现免费观看毛片| 日韩在线高清观看一区二区三区| 亚洲精品乱久久久久久| 国产亚洲精品第一综合不卡 | 18禁观看日本| 国产精品.久久久| 久久人人爽人人爽人人片va| 纯流量卡能插随身wifi吗| 久久久久久久久久人人人人人人| 色哟哟·www| 亚洲三级黄色毛片| 久久99蜜桃精品久久| 国产精品久久久久成人av| 在线精品无人区一区二区三| 丰满迷人的少妇在线观看| 国产精品麻豆人妻色哟哟久久| 国产极品粉嫩免费观看在线 | 王馨瑶露胸无遮挡在线观看| 日本欧美视频一区| 超色免费av| xxxhd国产人妻xxx| 亚洲久久久国产精品| 国产精品三级大全| 少妇的逼水好多| 人人妻人人澡人人看| 国产老妇伦熟女老妇高清| 亚洲第一区二区三区不卡| 亚洲精品久久午夜乱码| 亚洲丝袜综合中文字幕| 亚洲内射少妇av| 国产深夜福利视频在线观看| 国产成人精品久久久久久| 91精品国产国语对白视频| 麻豆成人av视频| 最近中文字幕高清免费大全6| 久久久欧美国产精品| 一级毛片 在线播放| 一级爰片在线观看| 国产极品粉嫩免费观看在线 | 欧美日韩av久久| 亚洲精品自拍成人| 久久午夜福利片| 国产伦理片在线播放av一区| 久久久久精品久久久久真实原创| 亚洲精品成人av观看孕妇| 欧美亚洲 丝袜 人妻 在线| 色婷婷久久久亚洲欧美| kizo精华| 日韩大片免费观看网站| 久久精品国产亚洲网站| 夜夜看夜夜爽夜夜摸| 国模一区二区三区四区视频| 99久久精品一区二区三区| 精品久久久久久久久亚洲| 成年人午夜在线观看视频| 大香蕉久久网| 国产男人的电影天堂91| 99九九在线精品视频| 午夜福利,免费看| 桃花免费在线播放| 少妇 在线观看| 伦理电影大哥的女人| 各种免费的搞黄视频| 丰满饥渴人妻一区二区三| 肉色欧美久久久久久久蜜桃| 高清不卡的av网站| 久久国产精品男人的天堂亚洲 | a级毛片免费高清观看在线播放| 美女国产高潮福利片在线看| 国产在线免费精品| 少妇的逼好多水|