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

    辣椒ASR基因的克隆與表達分析

    2021-03-05 11:56:03蔡小桃張穎韋麗麗徐小萬謝炳春莊澤岳吳智明
    南方農(nóng)業(yè)學報 2021年11期
    關鍵詞:非生物脅迫基因表達辣椒

    蔡小桃 張穎 韋麗麗 徐小萬 謝炳春 莊澤岳 吳智明

    辣椒ASR基因的克隆與表達分析

    蔡小桃1,張 穎1,韋麗麗1,徐小萬2,謝炳春1,莊澤岳1,吳智明1*

    (1仲愷農(nóng)業(yè)工程學院園藝園林學院,廣州 510225;2廣東省農(nóng)業(yè)科學院蔬菜研究所/廣東省蔬菜新技術(shù)研究重點實驗室,

    廣州 510640)

    摘要:【目的】克隆辣椒ASR(ABA-stress-ripening)基因進行生物信息學分析并探究其在不同組織和非生物脅迫處理下的表達模式,為深入研究ASR基因家族在辣椒果實發(fā)育和逆境脅迫中的作用機制提供參考。【方法】以遵辣1號辣椒為材料克隆其ASR基因,對該基因家族成員進行生物信息學分析和系統(tǒng)進化分析,運用實時熒光定量PCR檢測其在不同組織及在ABA、高溫、高鹽和干旱脅迫處理下的表達模式?!窘Y(jié)果】克隆的2個辣椒ASR基因分別命名為CaASR2和CaASR3,基因全長分別為987和938 bp,最長開放閱讀框(ORF)分別為333和339 bp,分別編碼110和112個氨基酸,蛋白分子量分別為12.45和12.73 kD,理論等電點(pI)分別為7.47和6.50,均為親水性蛋白。結(jié)合前人克隆的辣椒CaASR1基因分析,結(jié)果顯示,CaASRs預測位于細胞核,均含有ABA/WDS特征結(jié)構(gòu)域;基因結(jié)構(gòu)分析結(jié)果顯示,辣椒ASR基因均含有2個外顯子和1個內(nèi)含子,啟動子區(qū)域均包含ABA響應元件ABRE。CaASR1、CaASR2和CaASR3氨基酸序列分別與番茄(Solanum lycopersicum)基因組中的ASR蛋白氨基酸序列NP_001269248.1(80.20%)、NP_001307920.1(91.30%)和NP_001234137.1(96.43%)有非常高的相似性,表明其進化的相對保守性。3個ASR基因在辣椒花中表達量最低,在根、莖、葉和種子中表達量中等,但均隨著果實的發(fā)育表達量逐漸升高,特別是在轉(zhuǎn)色和成熟果實中。ABA、高溫、高鹽和干旱脅迫處理下,基因的表達量均顯著升高(P<0.05),其中以CaASR3響應最快,且上調(diào)表達量高于其他同源基因?!窘Y(jié)論】辣椒基因組中3個ASR基因具有ASR基因家族的典型特征,并在辣椒不同組織和非生物脅迫響應中發(fā)揮重要作用。

    關鍵詞: 辣椒;ASR;果實發(fā)育;非生物脅迫;基因表達

    中圖分類號: S641.3? ? ? ? ? ? ? ? ? ? ? ? ?文獻標志碼: A 文章編號:2095-1191(2021)11-3049-10

    Cloning and expression analysis of ASR genes in pepper

    (Capsicum annuum L.)

    CAI Xiao-tao1, ZHANG Ying1, WEI Li-li1, XU Xiao-wan2, XIE Bing-chun1,

    ZHUANG Ze-yue1, WU Zhi-ming1*

    (1College of Horticulture and Landscape Architecture, Zhongkai University of Agriculture and Engineering, Guangzhou

    510225, China; 2Vegetable Research Institute/Guangdong Key Laboratory for New Technology Research

    of Vegetables, Guangdong Academy of Agricultural Sciences, Guangzhou? 510640, China)

    Abstract:【Objective】ASR(ABA-stress-ripening) genes were cloned in this study. By bioinformatics analysis and expression analysis under different tissues and abiotic stress treatments, to provide a reference for further study on the mecha-nism of ASR gene family in pepper fruit development and abiotic stress. 【Method】ASR genes were cloned from Zunla-1 cultivar. A series of bioinformatics and phylogenetic analysis were carried out to evaluate the family members. The tissue specific expression patterns and the expression patterns under ABA, high temperature, high salt, and drought stress were detected by real-time fluorescence quantitative PCR. 【Result】Two ASR genes were cloned, named as CaASR2 and CaASR3. The full lengths of CaASR2 and CaASR3 were 987 and 938 bp, and the open reading frame(ORF) were 333 and 339 bp in length, which encoded 110 and 112 amino acids, respectively. Their molecular weights were 12.45 and 12.73 kD, and their theoretical isoelectric points(pI) were 7.47 and 6.50, respectively, and they were hydrophilic protein. The results showed that CaASRs were predicted to be in the nucleus and contained ABA/WDS(ABA/water de?cit stress) domain. All CaASRs genes contained two exons and one intron, and the promoter region contained ABA response element ABRE. The amino acid sequences of CaASR1, CaASR2 and CaASR3 were highly similar to? those of ASR protein amino acid NP_001269248.1(80.20%), NP_001307920.1(91.30%) and NP_001234137.1(96.43%) in tomato(Solanum lycopersicum), which indicated the relative conservation of evolution. The results of tissue specific expression analysis showed that they were expressed the lowest in flower, and the moderate in root, stem, leaf, and seed, gradually increased with fruit grew, and especially highly expressed in color-changing fruit and mature fruit. The expression levels were significantly increased under ABA stress, high temperature, high salt and drought(P<0.05), and CaASR3 had the fastest response, whose up regulated expression levels were higher than that of other homologous genes. 【Conclusion】Three ASR genes in pepper genome have typical characteristics of the ASR family and play important roles in different tissue and abiotic stress resistance in pepper.

    Key words: pepper(Capsicum annuum L.); ASR(ABA-stress-ripening); fruit development; abiotic stress; gene expression

    Foundation item: National Natural Science Foundation of China(32072598); The Characteristic Innovation Project of Colleges and Universities in Guangdong(2018KTSCX099)

    0 引言

    【研究意義】辣椒(Capsicum annuum L.)富含維生素C和辣椒素,營養(yǎng)價值極高,品種資源豐富,栽培適應性廣,在保證我國蔬菜周年均衡供應、農(nóng)民增收與脫貧致富中發(fā)揮了重要作用。2018年我國辣椒播種面積已達213.3萬ha,年產(chǎn)值約2500億元,是我國栽培面積最大和經(jīng)濟效益最高的蔬菜作物之一(鄒學校等,2020)。由于栽培范圍廣,加上近年來氣候的劇烈變化,高溫、干旱及土壤鹽漬化等逆境嚴重影響了辣椒栽培生產(chǎn),帶來巨大損失。ASR(Abscisic acid-,stress- and ripening-induced)蛋白是一類植物特有的親水性蛋白,在植物的生長、發(fā)育與逆境脅迫響應中發(fā)揮重要作用(González and Iusem,2014)。克隆并深入研究ASR基因的特征、進化和表達模式,對辣椒豐產(chǎn)和抗逆性改良具有重要意義?!厩叭搜芯窟M展】ASR基因家族是從種子植物進化而來,首次在番茄果實中被鑒定(Iusem et al.,1993),隨后在單子葉和雙子葉植物均有發(fā)現(xiàn)ASR蛋白,海岸生長的鹽生植物海蓬子(Salicornia brachiata)也鑒定出含ASR基因(Jha et al.,2012),這些同源基因主要受水分虧缺、鹽脅迫、激素脫落酸(ABA)、低溫和光強的調(diào)控。在不同物種中,ASR基因的數(shù)量不同,如水稻(Oryza sativa)ASR基因家族有6個已知成員(Philippe et al.,2010),番茄(Solanum lycopersicum)有5個成員(Frankel et al.,2006),葡萄(Vitis vinifera)僅1個已知成員(?akir et al.,2003),玉米(Zea mays)則多達9個成員(Virlouvet et al.,2011);而在十字花科植物如擬南芥(Arabidopsis thaliana)中缺失(Carrari et al.,2004)。研究發(fā)現(xiàn),ASR蛋白大多包含2個高度保守的區(qū)域:N端富含6~7個組氨酸殘基區(qū)域,C端約70個氨基酸殘基和核定位信號肽區(qū)域,核定位信號NLS已被證明無功能(Ricardi et al.,2012),亞細胞定位研究證實大多數(shù)植物ASR蛋白主要存在于細胞核和細胞質(zhì)中。ASR蛋白的另一個顯著特征是包含1個特有的保守ABA/WDS(ABA/water-de?cit stress)氨基酸結(jié)構(gòu)域,其以基因家族的形式存在。大多數(shù)植物ASR基因很短,編碼低分子量蛋白,通常作為糖、脫落酸和脅迫信號通路的共同成分,并參與果實發(fā)育和糖代謝等植物生命活動的許多過程,但ASR基因的生理作用和調(diào)控機制尚未明確(程維舜等,2013)。ASR蛋白還作為轉(zhuǎn)錄因子影響植物代謝,如水稻中OsASR5蛋白調(diào)節(jié)保護細胞免受Al誘導的應激反應(Arenhart et al.,2013)。Huang等(2020)發(fā)現(xiàn)CaASR1(Capana04g001289)作為一個正調(diào)控因子能與CabZIP63相互作用,加強辣椒對青枯病的抗性。對番茄ASR基因進行的表觀遺傳試驗證明,水分脅迫會引起表觀遺傳標記的變化(González et al.,2011)?!颈狙芯壳腥朦c】研究辣椒ASR基因家族的生物學功能,對于闡述辣椒果實成熟和對不利環(huán)境的耐受性機制具有重要意義。辣椒中有關ASR基因的研究鮮見報道,因此本課題組前期基于測序的辣椒基因組信息,通過同源性比對方法,從遵辣1號辣椒基因組中鑒定到3個ASR基因成員(基因ID分別為Capana04g001287、Capana04g001288和Capana04g00 1289),并初步分析了其在果實發(fā)育中的表達(張圣旭,2016)?!緮M解決的關鍵問題】利用生物信息學手段全面分析ASR基因家族的特征,檢測其在不同組織及在ABA處理和多種非生物脅迫下的表達模式,為后續(xù)研究ASR基因家族在辣椒生長發(fā)育及響應非生物脅迫中的功能機制提供理論依據(jù)。

    1 材料與方法

    1. 1 試驗材料

    辣椒自交系ZY35和遵辣1號品種分別由廣東省農(nóng)業(yè)科學院蔬菜研究所和遵義市農(nóng)業(yè)科學研究院辣椒研究所提供。植物RNA提取試劑盒購自北京華越洋生物科技有限公司;2×Hieff Robust PCR Master Mix,Hifair Ⅲ 1st Strand cDNA Synthesis SuperMix for qPCR反轉(zhuǎn)錄試劑盒,2×Hieff Canace Gold PCR Master Mix高保真酶預混液,以及Hieff clone zero TOPO-Blunt Cloning Kit等試劑盒購自上海翊圣生物科技有限公司;大腸桿菌DH5α感受態(tài)細胞購自上海唯地生物技術(shù)有限公司;膠回收試劑盒E.Z.N.A. Gel Extraction Kit和質(zhì)粒提取試劑盒E.Z.N.A. Plasmid DNA Mini KitⅠprotocal-spin protocol購自美國Omega Bio-Tek公司;熒光定量試劑盒SYBR Premix Ex TaqTM購自日本TaKaRa公司;其他生化試劑均為國產(chǎn)分析純。主要儀器設備:5810 R高速冷凍離心機(Eppendorf,德國);PCR儀(Bio-Rad,美國);4100凝膠成像儀(上海天能科技有限公司);熒光定量PCR儀CFX ConnectTM Real-Time System(Bio-Rad,美國)。

    1. 2 樣品采集及脅迫處理

    辣椒自交系ZY35于2020年8月進行播種,9月露地定植于廣東省農(nóng)業(yè)科學院白云試驗基地,盛果期從5個單株上采取根、莖、葉、花、幼果、轉(zhuǎn)色期果實、成熟果和種子等組織,每個組織樣品分成3份,用于熒光定量分析。于2021年1月將遵辣1號種子浸泡過夜,在28 ℃恒溫培養(yǎng)箱內(nèi)催芽,種子發(fā)芽后播種于塑料營養(yǎng)缽,在塑料大棚內(nèi)養(yǎng)護,待長至6片真葉時(約40 d),選取生長勢良好、長勢基本一致的幼苗進行脅迫處理。試驗設鹽脅迫(200 mmol/L NaCl)、干旱脅迫(400 mmol/L甘露醇)、高溫脅迫(42 ℃)和ABA脅迫(30 μmol/L ABA)4個處理,以正常栽培下清水澆灌處理作為對照(CK)。每組處理10棵幼苗,每個處理設3次重復。處理0、3、6、12和24 h后取辣椒第3~4片葉用于非生物脅迫后續(xù)基因定量分析。辣椒組織樣品經(jīng)液氮速凍后立即放入-80 ℃超低溫冰箱保存?zhèn)溆谩?/p>

    1. 3 基因克隆與測序

    自交系ZY35辣椒果實總DNA的提取參照CTAB法(王關林和方宏筠,2014)。參照試劑盒說明書提取不同組織和逆境處理下葉片的總RNA,經(jīng)1%瓊脂糖凝膠電泳檢測其質(zhì)量后利用試劑盒進行反轉(zhuǎn)錄,合成的cDNA保存在-20 ℃冰箱備用。根據(jù)遵辣1號基因組網(wǎng)站(https://www.solgenomics.net/)下載已有的辣椒ASR基因(登錄號Capana04g001287和Capana04g001288)序列信息(Qin et al.,2014),利用Primer Premier 5.0設計合成基因全長擴增的上、下游引物(表1),分別以自交系ZY35辣椒果實的DNA和cDNA為模板進行PCR擴增。PCR反應體系:模板2.0 μL,上、下游引物各2.5 μL,Gold Mix 25.0 μL;ddH2O補至50.0 μL。擴增程序:98 ℃預變性3 min;98 ℃變性10 s,68 ℃退火30 s,進行33個循環(huán);72 ℃延伸5 min。PCR擴增產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢測,確定條帶大小正確,進行膠回收。根據(jù)試劑盒操作說明,將純化產(chǎn)物連接到pESI-Blunt載體并轉(zhuǎn)化大腸桿菌DH5α感受態(tài)細胞,挑取單克隆在LB液體培養(yǎng)基中培養(yǎng),進行菌液PCR鑒定,挑選陽性克隆送至生工生物工程(上海)股份有限公司測序。測序結(jié)果使用DNAMAN 6.0進行比對。

    1. 4 基因表達分析

    利用實時熒光定量PCR對基因進行表達分析,以Actin為內(nèi)參,引物見表1。反應體系:TB Green Premix Ex Taq(Tli RNaseH Plus)10.0 μL,模板cDNA 2.0 μL,上、下游引物各0.5 μL,ddH2O補至20.0 μL。擴增程序:95 ℃ 30 s;95 ℃ 5 s,55 ℃ 30 s,進行40個循環(huán);從65 ℃到95 ℃每0.5 ℃依次讀板,繪制溶解曲線。每個樣品在在同一板中重復3次,作為技術(shù)重復。

    1. 5 辣椒ASR基因家族生物信息學分析

    利用ExPASy(https://www.expasy.org/protparam)在線分析工具ProtParam統(tǒng)計辣椒ASR家族蛋白的分子量、理論等電點(pI)、穩(wěn)定性及親/疏水性等理化性質(zhì)。利用NCBI數(shù)據(jù)庫CDD(Conserved domain database,https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi)預測蛋白質(zhì)保守結(jié)構(gòu)域。利用GSDS網(wǎng)站(http://gsds.cbi.pku.edu.cn/)繪制基因結(jié)構(gòu)圖,在BUSCA(busca.biocomp.unibo.it)在線網(wǎng)站進行亞細胞定位預測。運用PRABI和SWISS-MODEL預測ASR蛋白的二、三級結(jié)構(gòu)。提取辣椒ASR基因家族成員起始密碼子上游2000 bp基因組序列,利用Plant CARE(http://bioinformatics.psb.ugent.be/web-tools/plantcare)網(wǎng)站在線分析啟動子區(qū)域順式作用元件,之后在GSDS網(wǎng)站進行啟動子順式作用元件可視化。

    為分析辣椒ASR蛋白與不同物種間的親緣關系及系統(tǒng)進化關系,從NCBI數(shù)據(jù)庫中獲取水稻、玉米、番茄和葡萄的ASR家族蛋白質(zhì)序列。利用MEGA-X軟件進行系統(tǒng)進化分析,采取鄰接法(Neighbor joining,NJ)構(gòu)建系統(tǒng)發(fā)育進化樹,校驗參數(shù)Bootstrap method設置為1000。應用在線網(wǎng)站MEME(http://meme-suite.org)預測分析所有蛋白序列的結(jié)構(gòu)域,搜尋motif值為7,結(jié)構(gòu)域?qū)挾仍O定最小為10、最大為100,其他設定為默認參數(shù)。最后通過TBtools V1.09854(Chen et al.,2020)進行可視化拼接。

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

    采集的數(shù)據(jù)采用2-ΔΔCt法計算基因的相對表達量(Livak and Schmittgen,2002)。使用SPSS 23.0對數(shù)據(jù)進行顯著性分析,采用Origin 2019制圖。

    2 結(jié)果與分析

    2. 1 辣椒ASR基因的克隆與測序

    根據(jù)辣椒ASR基因的CDS序列和基因組序列,設計定量PCR特異引物。分別以cDNA和DNA為模板,對目標基因進行PCR擴增,獲得包含辣椒ASR基因完整開放閱讀框(ORF)的全長基因(圖1)。將擴增的片段克隆后測序,將測序序列與辣椒數(shù)據(jù)庫中ORF參考序列進行比對,結(jié)果表明,該基因的ORF序列和全長序列與數(shù)據(jù)庫中的序列完全一致。將登錄號為Capana04g001287的基因命名為CaASR2,Capana04g 001288命名為CaASR3。CaASR2基因全長為987 bp,最大ORF長333 bp,CaASR3全長938 bp,ORF長339 bp。

    2. 2 ASR編碼蛋白的序列特征

    通過ExPASy工具進行理化性質(zhì)分析,結(jié)果(表2)顯示,CaASR1、CaASR2和CaASR3分別編碼279、110和112個氨基酸,分子量分別為30.49、12.45和12.73 kD,pI分別為5.17、7.47和6.5,親水性平均系數(shù)均為負值,均為親水性蛋白。亞細胞定位預測結(jié)果顯示辣椒3個ASR蛋白均位于細胞核內(nèi)。利用NCBI數(shù)據(jù)庫中的保守結(jié)構(gòu)域預測工具CDD MbACO2蛋白保守結(jié)構(gòu)域,結(jié)果(圖2)發(fā)現(xiàn)其含保守結(jié)構(gòu)域PLN02299,故推測MbA-CO2屬于ACO家族成員。多重比對結(jié)果(圖3)顯示,辣椒ASR基因C端均含有保守ABA/WDS結(jié)構(gòu)域。

    蛋白質(zhì)一級結(jié)構(gòu)分析結(jié)果顯示,CaASR1蛋白占比在前5位的氨基酸依次為Gly(G)、Glu(E)、Thr(T)、Lys(K)和Ala(A),分別占15.8%、14.3%、12.5%、11.1%和7.2%。CaASR2蛋白占比在前5位的氨基酸依次為His(H)、Lys(K)、Ala(A)、Glu(E)和Gly(G),分別占22.7%、15.5%、15.5%、12.7%和8.2%。CaASR3蛋白占比在前5位的氨基酸依次為His(H)、Glu(E)、Ala(A)、Lys(K)和Gly(G),分別占19.6%、17.0%、16.1%、15.2%和8.0%。二級結(jié)構(gòu)預測分析結(jié)果(圖4)表明,CaASR1蛋白包含:67個α-螺旋,占24.01%;20個延伸鏈,占7.17%;12個β-轉(zhuǎn)角,占4.30%;180個無規(guī)則卷曲,占64.52%。CaASR2蛋白包含:64個α-螺旋,占58.18%;8個延伸鏈,占7.27%;10個β-轉(zhuǎn)角,占9.09%;28個無規(guī)則卷曲,占25.45%。CaASR3蛋白包含:65個α-螺旋,占56.25%;10個延伸鏈,占8.93%;8個β-轉(zhuǎn)角,占7.14%;31個無規(guī)則卷曲,占27.68%。利用SWISS-MODEL對ASR蛋白進行同源建模,獲得其三級結(jié)構(gòu)模型如圖5。

    2. 3 辣椒ASR基因結(jié)構(gòu)與啟動子順式作用元件預測

    利用GSDS網(wǎng)絡對ASR基因成員的結(jié)構(gòu)進行分析,發(fā)現(xiàn)辣椒中的3個ASR基因均含有2個外顯子和1個內(nèi)含子(圖6-A)。對辣椒ASR基因家族啟動子區(qū)域的順式作用元件進行預測,結(jié)果發(fā)現(xiàn),3個辣椒ASR基因除包含真核生物共有的順式調(diào)控元件(TATA框和CAAT框)外,還包括一些與脅迫響應相關的順式調(diào)控元件,即ABA響應元件ABRE及厭氧誘導(ARE)、低溫(LTR)、干旱誘導相關MYB結(jié)合位點(MBS)等元件(圖6-B),這些順式作用元件決定了該基因?qū)δ婢趁{迫的響應機制。

    2. 4 辣椒ASR蛋白質(zhì)進化和保守基序分析

    為了解辣椒ASR基因家族的進化關系,通過NCBI數(shù)據(jù)庫下載玉米、水稻、番茄和葡萄ASR家族蛋白序列,通過MEGA-X軟件構(gòu)建辣椒與玉米、水稻、番茄和葡萄ASR家族蛋白的系統(tǒng)進化樹(圖7-A)。結(jié)果顯示,24個ASR蛋白序列相似,與辣椒CaASRs基因遺傳背景相近的是番茄和葡萄2個物種,親緣關系較遠的是水稻和玉米。同時CaASR2和CaASR3的親緣關系較近。CaASR1、CaASR2和CaASR3氨基酸序列分別與番茄基因組中的ASR蛋白氨基酸序列NP_001269248.1(80.20%)、NP_001307920.1(91.3%)和NP_001234137.1(96.43%)有非常高的相似性。保守基序統(tǒng)計分析,共識別到7個motif(保守基序),并依次命名為motif1~motif7(圖7-B)。其中motif1~ motif4最為保守,在所有辣椒ASR基因中均存在。CaASR2和CaASR3功能序列的數(shù)量和位點基本一致,具有高度保守性,推測其功能單位基本相似。進一步分析發(fā)現(xiàn),3個基因均含有motif1~motif4,CaASR1還含有motif5和motif6。表明同一進化簇成員間具有相似的保守結(jié)構(gòu)域組成,而不同簇具有各自特異的保守結(jié)構(gòu)域,可能與其具體的生物學功能相關。

    2. 5 CaASRs基因的組織特異性表達分析

    為充分了解CaASR基因在辣椒不同組織中的表達模式,分別以自交系ZY35辣椒的根、莖、葉、花、幼果、轉(zhuǎn)色期果實、成熟果實和種子的cDNA為模板,對CaASRs基因表達量進行實時熒光定量PCR分析,結(jié)果(圖8)顯示,3個ASR基因呈現(xiàn)相對一致的組織特異性表達,即在辣椒花中表達量最低,其次是在根中表達量相對較低,在莖、葉和種子中表達量中等,在果實中最高。CaASR1隨果實發(fā)育表達量逐漸升高,在成熟果實中表達量比幼果增加30倍。CaASR2和CaASR3均在轉(zhuǎn)色期果實中表達量最高,在成熟果實中的表達量反而顯著降低(P<0.05,下同)。說明CaASR2和CaASR3在果實發(fā)育早期被誘導表達,CaASR1在果實成熟后期表現(xiàn)高表達量。綜上所述,推測ASR基因可能參與辣椒的生長發(fā)育過程,特別是在辣椒果實發(fā)育和成熟過程中可能發(fā)揮重要作用。

    2. 6 CaASRs基因在非生物脅迫下的表達分析

    對辣椒ASR基因在42 ℃高溫、甘露醇誘導的干旱脅迫、ABA噴施和NaCl脅迫處理下葉片中表達量進行分析,結(jié)果(圖9)表明,辣椒在4組非生物脅迫下,3個ASR基因均表達量上調(diào),表達模式存在差異,其中,CaASR3的表達量相對于其他2個基因上調(diào)更多。CaASR1在處理后表達水平先上調(diào),隨后恢復至對照水平,表達量均在12 h時顯著高于其他時間段,在4組不同處理中呈現(xiàn)相對一致的表達模式;CaASR2在ABA和干旱處理下,表達量在6 h時顯著高于其他時間段,在鹽脅迫和高溫處理下,24 h時仍具有相對較高的表達量。CaASR3基因在處理3 h后就開始表達量上調(diào),響應快于其他同源基因;12 h后已經(jīng)下降,恢復至對照水平,在4組處理中表達量趨勢也呈現(xiàn)相對一致。推測CaASRs在辣椒抵御鹽脅迫、干旱脅迫及氧化脅迫中可能發(fā)揮正向調(diào)控功能。

    3 討論

    ASR基因作為一個小的基因家族,最早從成熟的番茄果實中分離,其在果實的發(fā)育與成熟、非生物脅迫和ABA應答中發(fā)揮了重要作用,在馬鈴薯(Dóczi et al.,2005)、草莓(Chen et al.,2011)、番茄(Ricardi et al.,2014)、玉米(Li et al.,2018)和水稻(Park et al.,2020)等多種作物中均有報道。辣椒基因組中包含3個ASR基因,均含1個內(nèi)含子和2個外顯子的典型結(jié)構(gòu)和保守的ABA/WDS氨基酸結(jié)構(gòu)域,基因結(jié)構(gòu)簡單且相似,并呈串聯(lián)重復分布于辣椒第4號染色體上。進化分析結(jié)果表明,辣椒ASR基因與番茄基因組中的ASR蛋白氨基酸序列同源性最近。3個ASR蛋白預測均位于細胞核內(nèi),蛋白質(zhì)理化性質(zhì)相似,結(jié)構(gòu)簡單,該結(jié)果與Golan等(2014)研究所得的ASR序列特征相似,辣椒中該基因家族具有ASR蛋白的主要特征。

    ASR基因在植物中存在組織特異性表達,不同成員可能在不同的器官、環(huán)境條件和表達模式中差異表達(Wang et al.,2016;Li et al.,2020)。組織特異性表達分析結(jié)果表明,CaASRs基因在莖、葉和種子中表達量高,在根中表達量相對較低,在花中表達量最低,該基因家族3個成員在正常生理條件下具有廣泛的表達模式。ASR基因在不同組織中的表達呈現(xiàn)出不同的表達模式,暗示了ASR基因家族功能的復雜性。然而CaASR1隨果實發(fā)育表達量逐漸升高,在成熟果實中表達量達最高值;CaASR2和CaASR3均在轉(zhuǎn)色期果實中表達量最高,高于其他組織。?akir等(2003)通過酵母單雜交試驗表明,葡萄ASR蛋白與己糖轉(zhuǎn)運基因的啟動子結(jié)合,作為參與糖代謝的轉(zhuǎn)錄因子發(fā)揮作用。Chen等(2011)發(fā)現(xiàn)草莓果實中增加內(nèi)源ABA含量,F(xiàn)aASR在轉(zhuǎn)錄和翻譯水平的表達增加,可能對草莓果實成熟有一定的促進作用。辣椒是非呼吸躍變型果實,已有研究證明ABA含量隨糖分積累和果實成熟逐漸積累,導致果實成熟過程中ASR基因大量表達。ASR基因在辣椒果實特異高表達,推測其在辣椒果實發(fā)育與成熟過程中可能發(fā)揮重要作用。

    關于ASR基因在逆境脅迫中的研究較多。最早從番茄中分離的ASR基因顯示,該基因的表達受脫水處理、ABA、茉莉酸甲酯和水楊酸的誘導,導致植物體對水和滲透脅迫有更高的耐受性。在煙草中過量表達番茄ASR1增強了煙草對高鹽脅迫的耐受性(Fischer et al.,2011)。水稻ASR1(OsAsr1)基因受鹽、甘露醇和ABA誘導,OsAsr1轉(zhuǎn)入酵母中能增強酵母對NaCl、雙氧水、硫酸、水楊酸和熱等脅迫的耐受性(Kim et al.,2009)。香蕉MaASR1在擬南芥中表達后使植株變小、葉片變小且變厚、抗旱性增強,同時能增強擬南芥對ABA信號的敏感性(苗紅霞等,2014)。在擬南芥中過表達玉米ZmASR3基因能顯著提高植株抗?jié)B透脅迫能力,增加ABA含量(Liang et al.,2019)。通過分析CaASRs在ABA、干旱、高鹽和高溫脅迫下的表達模式,發(fā)現(xiàn)多種處理時,ASR基因在辣椒葉片中的表達量整體呈升高趨勢,表達受非生物脅迫的共同調(diào)控,與辣椒的逆境響應密切相關。已有研究證明ABA處理可在轉(zhuǎn)錄和翻譯水平上誘導ASR的表達,上游順式作用元件分析和RT-qPCR結(jié)果表明ABA參與非生物脅迫下辣椒ASR基因表達的信號轉(zhuǎn)導。因此,推斷CaASRs基因可提高辣椒對非生物脅迫的耐受性,但有待進一步通過轉(zhuǎn)CaASRs驗證其功能。

    4 結(jié)論

    辣椒基因組中含3個ASR基因,是一個結(jié)構(gòu)簡單的小基因家族,進化非常保守。辣椒ASR基因的表達具有組織特異性,在果實中表達量最高。在ABA、高鹽、干旱和高溫等非生物逆境脅迫下,CaASRs基因被誘導表達,推測ASR基因家族在辣椒果實成熟和抗逆過程起重要作用。

    參考文獻:

    程維舜,孫玉宏,曾紅霞,杜念華,施先鋒,蔡新忠. 2013. ASR蛋白與植物的抗逆性研究進展[J]. 園藝學報,40(10):2049-2057. [Cheng W S,Sun Y H,Zeng H X,Du N H,Shi X F,Cai X Z. 2013. ASR protein and plant stress to-lerance[J]. Acta Horticulturae Sinica,40(10):2049-2057.] doi:10.16420/j.issn.0513-353x.2013.10.004.

    苗紅霞,王園,徐碧玉,劉菊華,賈彩紅,張建斌,王卓,孫佩光,金志強. 2014. 香蕉MaASR1基因的抗干旱作用[J]. 植物學報,49(5):548-559. [Miao H X,Wang Y,Xu B Y,Liu J H,Jia C H,Zhang J B,Wang Z,Sun P G,Jin Z Q. 2014. The role of banana MaASR1 in drought stress tolerance[J]. Chinese Bulletin of Botany,49(5):548-559.] doi:10.3724/SP.J.1259.2014.00548.

    王關林,方宏筠. 2014. 植物基因工程[M]. 第二版. 北京:科學出版社. [Wang G L,F(xiàn)ang H J. 2014. Plant genetic engineering[M]. The 2nd Edition. Beijing:Science Press.]

    張圣旭. 2016. ABA信號在‘遵辣1號’果實發(fā)育及成熟過程中的作用研究[D]. 廣州:仲愷農(nóng)業(yè)工程學院. [Zhang S X. 2016. Study on the effect of ABA signal in the fruit development and maturation process in ‘Zunla-1’[D]. Guangzhou:Zhongkai University of Agriculture and Enginee-ring.]

    鄒學校,馬艷青,戴雄澤,李雪峰,楊莎. 2020. 辣椒在中國的傳播與產(chǎn)業(yè)發(fā)展[J]. 園藝學報,47(9):1715-1726. [Zou X X,Ma Y Q,Dai X Z,Li X F,Yang S. 2020. Spread and industry development of pepper in China[J]. Acta Horticulturae Sinica,47(9):1715-1726.] doi:10.16420/j.issn.0513-353x.2020-0103.

    Arenhart R A,Lima J,Pedron M,Carvalho F E L,Silveira J A G,Rosa S B,Caverzan A,Andrade C M B,Schune-mann M,Margis R,Margis-Pinheiro M. 2013. Involvement of ASR genes in aluminium tolerance mechanisms in rice[J]. Plant,Cell & Environment,36(1):52-67. doi:10.1111/j.1365-3040.2012.02553.x.

    ?akir B,Agasse A,Gaillard C,Saumonneau A,Delrot S,Atanassova R. 2003. A grape ASR protein involved in sugar and abscisic acid signaling[J]. Plant Cell,15(9):2165-2180. doi:10.1105/tpc.013854.

    Carrari F,F(xiàn)ernie A R,Iusem N D. 2004. Heard it through the grapevine?ABA and sugar crosstalk:The ASR story[J]. Trends in Plant Science,9(2):57-59. doi:10.1016/j.tplants.2003.12.004.

    Chen J Y,Liu D J,Jiang Y M,Zhao M L,Shan W,Kuang J F,Lu W J. 2011. Molecular characterization of a strawberry FaASR gene in relation to fruit ripening[J]. PLoS ONE,6(9):e24649. doi:10.1371/journal.pone.0024649.

    Chen C,Chen H,Zhang Y,Thomas H R,F(xiàn)rank M H,He Y H,Xia R. 2020. TBtools:An integrative toolkit developed for interactive analyses of big biological data[J]. Molecular Plant,13(8):1194-1202. doi:10.1016/j.molp.2020.06. 009.

    Dóczi R,Kondrák M,Kovács G,Beczner F,Bánfalvi Z. 2005. Conservation of the drought inducible DS2 genes and divergences from their ASR paralogues in solanaceous species[J]. Plant Physiology and Biochemistry,43(3):269-276. doi:10.1016/j.plaphy.2005.02.002.

    Fischer I,Camus-Kulandaivelu L,Allal F,Stephan W. 2011. Adaptation to drought in two wild tomato species:The evolution of the Asr gene family[J]. New Phytologist,190(4):1032-1044. doi:10.1111/j.1469-8137.2011.03648.x.

    Frankel N,Carrari F,Hasson E,Allal F,Iusem N D. 2006. Evolutionary history of the Asr gene family[J]. Gene,378:74-83. doi:10.1016/j.gene.2006.05.010.

    Golan I,Dominguez P G,Konrad Z,Shkolnik-Inbar D,Carrari F,Bar-Zvi D. 2014. Tomato ABSCISIC ACID STRESS RIPENING(ASR) gene family revisited[J]. PLoS ONE,9(10):e107117. doi:10.1371/journal.pone.0107117.

    González R M,Iusem N D. 2014. Twenty years of research on Asr(ABA-stress-ripening) genes and proteins[J]. Planta,239(5):941-949. doi:10.1007/s00425-014-2039-9.

    González R M,Ricardi M M,Iusem N D. 2011. Atypical epigenetic mark in an atypical location:Cytosine methylation at asymmetric(CNN) sites within the body of a non-repetitive tomato gene[J]. BMC Plant Biology,11(1):94. doi:10.1186/1471-2229-11-94.

    Huang J F,Shen L,Sheng Y,Guan D Y,He S L. 2020. CaASR1 promotes salicylic acid- but represses jasmonic acid-dependent signaling to enhance the resistance of Capsicum annuum to bacterial wilt by modulating CabZIP63[J]. Journal of Experimental Botany t,71(20):6538-6554. doi:10.1093/jxb/eraa350.

    Iusem N D,Bartholomew D M,Hitz W D,Scolnik P A. 1993. Tomato(Lycopersicon esculentum) transcript Induced by water deficit and ripening[J]. Plant Physiology,102(4):1353-1354. doi:10.1104/pp.102.4.1353.

    Jha B,Lal S,Tiwari V,Yadav S K,Agarwal P K. 2012. The SbASR-1 gene cloned from an extreme halophyte Salicornia brachiata enhances salt tolerance in transgenic tobacco[J]. Marine Biotechnology,14(6):782-792. doi:10.1007/ s10126-012-9442-7.

    Kim S J,Lee S C,Hong S K,An K,An G,Kim S R. 2009. Ectopic expression of a cold-responsive OsAsr1 cDNA gives enhanced cold tolerance in transgenic rice plants[J]. Molecules and Cells,27(4):449-458. doi:10.1007/s10059-009-0055-6.

    Livak K J,Schmittgen T D. 2002. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCt method[J]. Methods,25(4):402-408. doi:10.1006/meth.2001.1262.

    Li H W,Guan H Y,Zhuo Q C,Wang Z S,Li S D,Si J S,Zhang B,F(xiàn)eng B,Kong L A,Wang F H,Wang Z,Zhang L S. 2020. Genome-wide characterization of the abscisic acid-,stress- and ripening-induced(ASR) gene family in wheat(Triticum aestivum L.)[J]. Biological Research,53(1):23. doi:10.1186/s40659-020-00291-6.

    Li X Y,Li L J,Zuo S Y,Li J,Wei S. 2018. Differentially expressed ZmASR genes associated with chilling tolerance in maize(Zea mays) varieties[J]. Functional Plant Biology,45(12):1173. doi:10.1071/FP17356.

    Liang Y,Jiang Y L,Du M,Li B Y,Chen L,Chen M C,Jin D M,Wu J D. 2019. ZmASR3 from the maize ASR gene family positively regulates drought tolerance in transgenic Arabidopsis[J]. International Journal of Molecular Scien-ces,20(9):2278. doi:10.3390/ijms20092278.

    Park S I,Kim J J,Shin S Y,Kim Y S,Yoon H S. 2020. ASR enhances environmental stress tolerance and improves grain yield by modulating stomatal slosure in rice[J]. Frontiers in Plant Science,10:1752. doi:10.3389/fpls.2019. 01752.

    Philippe R,Courtois B,McNally K L,Mournet P,El-Malki R,Le Paslier M C,F(xiàn)abre D,Billot C,Brunel D,Glaszmann J,This D. 2010. Structure,allelic diversity and selection of Asr genes,candidate for drought tolerance,in Oryza sativa L. and wild relatives[J]. Theoretical and Applied Genetics,121(4):769-787. doi:10.1007/s00122-010-1348-z.

    Qin C,Yu C S,Shen Y U,F(xiàn)ang X D,Chen L,Min J M,Cheng J W,Zhao S C,Xu M,Luo Y,Yang Y L,Wu Z M,Mao L K,Wu H Y,Ling-Hu C Y,Zhou H K,Lin H J,Gonzalez-Morales S I,Trejo-Saavedra D L,Tian H,Tang X,Zhao M J,Huang Z Y,Zhou A W,Yao X M,Cui J J,Li W Q,Chen Z,F(xiàn)eng Y Q,Niu Y C,Bi S M,Yang X W,Li W P,Cai H M,Luo X R,Montes-Hernandez S,Leyva-Gonzalez M A,Xiong Z Q,He X J,Bai L J,Tan S,Tang X Q,Liu D,Liu J W,Zhang S X,Chen M S,Zhang L,Zhang L,Zhang Y H,Liao W Q,Zhang Y,Wang M,Lv X D,Wen B,Liu H J,Luan H M,Zhang Y G,Yang S,Wang X D,Batley J,Li X Q,Li S C,Wang J Y,Palloix A,Bosland P,Li Y R,Krogh A,Rivera-Bustamante R F,Herrera-Estrella L,Yin Y,Yu J P,Hu K L,Zhang Z M. 2014. Whole-genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization[J]. Proceedings of the National Academy of Sciences,111(14):5135-5140. doi:10.1073/ pnas.1400975111.

    Ricardi M M,Guaimas F F,González R M,Burrieza H P,López-Fernández M P,Jares-Erijman E A,Estévez J M,Iusem N D. 2012. Nuclear import and dimerization of tomato ASR1,a water stress-inducible protein exclusive to plants[J]. PLoS One,7(8):e41008. doi:10.1371/journal.pone.0041008.

    Ricardi M M,González R M,Zhong S,Domínguez P G,Duffy T,Turjanski P G,Salgado Salter J D,Alleva K,Carrari F,Giovannoni J J,Estévez J M,Iusem N D. 2014. Genome-wide data(ChIP-seq) enabled identification of cell wall-related and aquaporin genes as targets of tomato ASR1,a drought stress-responsive transcription factor[J]. BMC Plant Biology,14(1):29. doi:10.1186/1471-2229-14-29.

    Virlouvet L,Jacquemot M,Gerentes D,Corti H,Bouton S,Gilard F,Valot B,Trouverie J,Tcherkez G,F(xiàn)alque M,Damerval C,Rogowsky P,Perez P,Noctor G,Zivy M,Coursol S. 2011. The ZmASR1 protein influences bran-chedchain amino acid biosynthesis and maintains kernel yield in maize under water-limited conditions[J]. Plant Physiology,157(2):917-936. doi:10.1104/pp.111.176818.

    Wang L,Hu W,F(xiàn)eng J,Yang X,Huang Q,Xiao J,Liu Y,Yang G,He G. 2016. Identification of the ASR gene family from Brachypodium distachyon and functional characterization of BdASR1 in response to drought stress[J]. Plant Cell Reports,35(6):1221-1234. doi:10.1007/s00299-016-1954-6.

    收稿日期:2021-05-06

    基金項目:國家自然科學基金項目(32072598);廣東省普通高校特色創(chuàng)新項目(2018KTSCX099)

    通訊作者:吳智明(1981-),https://orcid.org/0000-0001-8988-6060,博士,教授,主要從事蔬菜生物技術(shù)與遺傳育種研究工作,E-mail:wuzm2012@zhku.edu.cn

    第一作者:蔡小桃(1996-),https://orcid.org/0000-0001-7165-7075,研究方向為蔬菜生物技術(shù),E-mail:2392023499@qq.com

    猜你喜歡
    非生物脅迫基因表達辣椒
    原來,你還是這樣的辣椒
    你的辣椒結(jié)出果實了嗎?
    辣椒也瘋狂
    揀辣椒
    中外文摘(2020年9期)2020-06-01 13:47:56
    白木香轉(zhuǎn)錄因子AsMYB1和AsMYB2克隆及表達分析
    抗菌肽對細菌作用機制的研究
    基因芯片在胃癌及腫瘤球細胞差異表達基因篩選中的應用
    美洲大蠊提取液對大鼠難愈合創(chuàng)面VEGF表達影響的研究
    二甲基砷酸毒理學的研究進展
    非生物脅迫對擬南芥IQM4基因表達的影響
    科技視界(2016年16期)2016-06-29 11:55:38
    一级毛片久久久久久久久女| 欧美成人免费av一区二区三区| 男女下面进入的视频免费午夜| 国产亚洲精品av在线| 极品教师在线视频| 九九爱精品视频在线观看| 在线免费十八禁| 日韩亚洲欧美综合| 免费人成视频x8x8入口观看| 日韩强制内射视频| 九色成人免费人妻av| 亚洲熟妇熟女久久| 国产伦精品一区二区三区视频9| 日本五十路高清| 国产午夜精品久久久久久一区二区三区 | 成年免费大片在线观看| 啦啦啦啦在线视频资源| 亚洲一区高清亚洲精品| 一级黄片播放器| 国产麻豆成人av免费视频| 非洲黑人性xxxx精品又粗又长| 日本撒尿小便嘘嘘汇集6| 免费在线观看成人毛片| 真人一进一出gif抽搐免费| 午夜免费激情av| 国产精品久久视频播放| a在线观看视频网站| 日韩高清综合在线| 69av精品久久久久久| 欧美成人一区二区免费高清观看| 午夜福利视频1000在线观看| 久9热在线精品视频| 国产中年淑女户外野战色| 久久精品国产99精品国产亚洲性色| 日本在线视频免费播放| 乱码一卡2卡4卡精品| 两人在一起打扑克的视频| 精品日产1卡2卡| 高清在线国产一区| 国内精品久久久久精免费| 国产色婷婷99| 国产激情偷乱视频一区二区| 亚洲无线观看免费| 亚洲色图av天堂| 免费观看在线日韩| 欧美+亚洲+日韩+国产| 亚洲精品一区av在线观看| 久久国产乱子免费精品| 久久99热这里只有精品18| 99久久中文字幕三级久久日本| 亚洲熟妇熟女久久| 日本精品一区二区三区蜜桃| 88av欧美| 久久久久久久午夜电影| 日本a在线网址| 日本一二三区视频观看| 俺也久久电影网| 国产麻豆成人av免费视频| 成人二区视频| 日本在线视频免费播放| 国产综合懂色| 午夜日韩欧美国产| 少妇高潮的动态图| 999久久久精品免费观看国产| 国产精品久久久久久亚洲av鲁大| 精品久久久久久久久亚洲 | 91麻豆av在线| 国产高清三级在线| 最近最新中文字幕大全电影3| 看片在线看免费视频| 两性午夜刺激爽爽歪歪视频在线观看| 全区人妻精品视频| 日日啪夜夜撸| 色噜噜av男人的天堂激情| 国产免费一级a男人的天堂| 美女大奶头视频| 国产精品国产三级国产av玫瑰| 五月玫瑰六月丁香| 九九热线精品视视频播放| 日本a在线网址| 春色校园在线视频观看| 少妇人妻精品综合一区二区 | 九色国产91popny在线| 日日摸夜夜添夜夜添小说| 国产精品不卡视频一区二区| 久久精品国产亚洲av涩爱 | 国产精品自产拍在线观看55亚洲| 国产伦人伦偷精品视频| 亚洲欧美清纯卡通| 国产视频内射| 特大巨黑吊av在线直播| 久久久精品大字幕| 国产精品亚洲美女久久久| 国产精品亚洲一级av第二区| 国产精品综合久久久久久久免费| 亚洲精品一区av在线观看| 国产蜜桃级精品一区二区三区| 成年女人永久免费观看视频| 亚洲av.av天堂| 我的女老师完整版在线观看| 亚洲一级一片aⅴ在线观看| 久久久精品欧美日韩精品| 国产精品精品国产色婷婷| 岛国在线免费视频观看| 91在线精品国自产拍蜜月| 深爱激情五月婷婷| 国产黄色小视频在线观看| 春色校园在线视频观看| 99热6这里只有精品| 久久久久精品国产欧美久久久| 国产av麻豆久久久久久久| 亚洲av五月六月丁香网| 日韩一本色道免费dvd| 少妇被粗大猛烈的视频| 久久人人爽人人爽人人片va| 少妇裸体淫交视频免费看高清| 亚洲国产精品成人综合色| 欧美色视频一区免费| 午夜福利成人在线免费观看| 国产精品电影一区二区三区| 国产精品乱码一区二三区的特点| 永久网站在线| 亚洲真实伦在线观看| 国产免费av片在线观看野外av| 免费高清视频大片| 欧美日韩乱码在线| 1024手机看黄色片| 女人被狂操c到高潮| 尤物成人国产欧美一区二区三区| 欧美+日韩+精品| 国产视频一区二区在线看| av中文乱码字幕在线| 欧美精品国产亚洲| 欧美日韩乱码在线| 真人一进一出gif抽搐免费| 国产单亲对白刺激| 黄色丝袜av网址大全| a级一级毛片免费在线观看| 特级一级黄色大片| 亚洲中文字幕一区二区三区有码在线看| 亚洲乱码一区二区免费版| 长腿黑丝高跟| 丝袜美腿在线中文| 成年免费大片在线观看| 精品一区二区三区av网在线观看| 在线观看舔阴道视频| 一区二区三区激情视频| 国产真实伦视频高清在线观看 | 极品教师在线视频| 国内久久婷婷六月综合欲色啪| 成人欧美大片| 少妇高潮的动态图| 自拍偷自拍亚洲精品老妇| 亚洲熟妇中文字幕五十中出| 亚洲国产精品成人综合色| 三级毛片av免费| 亚洲第一电影网av| 久久久久免费精品人妻一区二区| 三级国产精品欧美在线观看| 美女大奶头视频| 九九热线精品视视频播放| 日本一二三区视频观看| 欧美成人免费av一区二区三区| 久久久久久久久中文| 一个人看的www免费观看视频| 精品人妻视频免费看| 中文字幕熟女人妻在线| 国产美女午夜福利| 性插视频无遮挡在线免费观看| xxxwww97欧美| 久久久久久久亚洲中文字幕| 又黄又爽又免费观看的视频| 最好的美女福利视频网| 国产精品人妻久久久影院| 日本a在线网址| 非洲黑人性xxxx精品又粗又长| a级一级毛片免费在线观看| 又紧又爽又黄一区二区| 2021天堂中文幕一二区在线观| www.色视频.com| 欧美xxxx黑人xx丫x性爽| 国产v大片淫在线免费观看| 看黄色毛片网站| 亚洲经典国产精华液单| av在线观看视频网站免费| 亚州av有码| 午夜亚洲福利在线播放| 一进一出抽搐动态| 又黄又爽又刺激的免费视频.| 国产精品国产高清国产av| 精品一区二区三区人妻视频| 搡女人真爽免费视频火全软件 | 欧美日本亚洲视频在线播放| 国产亚洲91精品色在线| 精品久久久久久久末码| 18禁黄网站禁片午夜丰满| 啦啦啦韩国在线观看视频| 免费av不卡在线播放| 国产单亲对白刺激| 欧美xxxx性猛交bbbb| 人妻丰满熟妇av一区二区三区| 欧美最新免费一区二区三区| 成人国产麻豆网| 欧美激情在线99| 中文字幕精品亚洲无线码一区| 波多野结衣巨乳人妻| 三级国产精品欧美在线观看| 天堂动漫精品| 性插视频无遮挡在线免费观看| 亚洲男人的天堂狠狠| 黄色配什么色好看| 亚洲精品一区av在线观看| 国产精品,欧美在线| 嫩草影院入口| 亚洲国产精品久久男人天堂| 日韩国内少妇激情av| 可以在线观看毛片的网站| 亚洲国产高清在线一区二区三| 午夜福利18| 露出奶头的视频| 99久久精品国产国产毛片| 干丝袜人妻中文字幕| 国产精品99久久久久久久久| av在线蜜桃| 亚洲av美国av| 日本 欧美在线| 男人狂女人下面高潮的视频| 国产精品亚洲美女久久久| 99国产极品粉嫩在线观看| 久久久久久国产a免费观看| 好男人在线观看高清免费视频| 国产成人a区在线观看| 身体一侧抽搐| 十八禁网站免费在线| 免费观看的影片在线观看| 99久久成人亚洲精品观看| 网址你懂的国产日韩在线| 精品久久久久久,| 日本a在线网址| 国产精品久久电影中文字幕| 色5月婷婷丁香| 神马国产精品三级电影在线观看| 哪里可以看免费的av片| 99国产精品一区二区蜜桃av| 久久精品国产99精品国产亚洲性色| 国产在视频线在精品| 亚洲中文字幕一区二区三区有码在线看| 亚洲内射少妇av| 少妇人妻精品综合一区二区 | 欧美+亚洲+日韩+国产| 极品教师在线免费播放| 两人在一起打扑克的视频| 国产伦一二天堂av在线观看| 一a级毛片在线观看| 啪啪无遮挡十八禁网站| 蜜桃久久精品国产亚洲av| 亚洲第一电影网av| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区在线观看日韩| 日韩av在线大香蕉| 尤物成人国产欧美一区二区三区| 精品不卡国产一区二区三区| 日本 av在线| 日韩欧美精品免费久久| 国产精品一区二区免费欧美| 日韩欧美 国产精品| 九色成人免费人妻av| 一进一出抽搐动态| 99久久精品一区二区三区| 国产精品久久久久久亚洲av鲁大| 九九热线精品视视频播放| 久久国内精品自在自线图片| 九色国产91popny在线| 成年女人看的毛片在线观看| 免费不卡的大黄色大毛片视频在线观看 | 国产精品人妻久久久久久| 亚洲美女黄片视频| 99九九线精品视频在线观看视频| 久久久国产成人精品二区| 熟妇人妻久久中文字幕3abv| 搡老岳熟女国产| 两个人视频免费观看高清| 日韩欧美在线二视频| 又紧又爽又黄一区二区| 亚洲国产欧洲综合997久久,| 国产精品美女特级片免费视频播放器| 国产精品野战在线观看| 亚洲成人久久性| 国产av不卡久久| 99热6这里只有精品| 国产成人一区二区在线| 91av网一区二区| 夜夜爽天天搞| 国产aⅴ精品一区二区三区波| 1000部很黄的大片| 久久精品综合一区二区三区| 亚洲人与动物交配视频| 赤兔流量卡办理| 午夜福利高清视频| 久久中文看片网| 给我免费播放毛片高清在线观看| 人人妻人人看人人澡| av天堂在线播放| 超碰av人人做人人爽久久| 女的被弄到高潮叫床怎么办 | 在线观看美女被高潮喷水网站| 欧美中文日本在线观看视频| 18+在线观看网站| 久久人人精品亚洲av| 男女之事视频高清在线观看| 蜜桃久久精品国产亚洲av| 亚洲人成网站在线播放欧美日韩| 国产精品一及| 久9热在线精品视频| 韩国av一区二区三区四区| 2021天堂中文幕一二区在线观| 国产高清有码在线观看视频| 亚洲不卡免费看| 亚洲国产高清在线一区二区三| 男女边吃奶边做爰视频| 成人鲁丝片一二三区免费| 国产精品国产高清国产av| 97超视频在线观看视频| 99精品久久久久人妻精品| 一进一出好大好爽视频| 色播亚洲综合网| 国产高清三级在线| 此物有八面人人有两片| 搡老妇女老女人老熟妇| 亚洲人成网站在线播| 免费观看在线日韩| 露出奶头的视频| 亚洲欧美精品综合久久99| 变态另类成人亚洲欧美熟女| 国产精品三级大全| 国产日本99.免费观看| 不卡视频在线观看欧美| 色综合亚洲欧美另类图片| 精品一区二区三区视频在线| 亚洲午夜理论影院| 大型黄色视频在线免费观看| 舔av片在线| 免费看日本二区| 九九在线视频观看精品| 欧美精品国产亚洲| 韩国av在线不卡| 欧美不卡视频在线免费观看| 99久久中文字幕三级久久日本| 国产精品美女特级片免费视频播放器| netflix在线观看网站| 丰满的人妻完整版| 久99久视频精品免费| 日韩,欧美,国产一区二区三区 | 99久久精品热视频| 国产精品嫩草影院av在线观看 | 人人妻,人人澡人人爽秒播| а√天堂www在线а√下载| 亚洲精品乱码久久久v下载方式| 99视频精品全部免费 在线| 国产男人的电影天堂91| 桃色一区二区三区在线观看| 国产精品99久久久久久久久| 国产精品久久久久久精品电影| 欧美xxxx性猛交bbbb| 国语自产精品视频在线第100页| 久久人人爽人人爽人人片va| 国产精品伦人一区二区| 成年女人看的毛片在线观看| a级毛片a级免费在线| 日韩欧美在线二视频| 91午夜精品亚洲一区二区三区 | 天天一区二区日本电影三级| АⅤ资源中文在线天堂| 一级av片app| 亚洲成人免费电影在线观看| 乱人视频在线观看| 高清在线国产一区| 久久久久久九九精品二区国产| 91在线精品国自产拍蜜月| 联通29元200g的流量卡| 久久国产乱子免费精品| 亚洲最大成人手机在线| 国产精品乱码一区二三区的特点| 欧美黑人欧美精品刺激| 特大巨黑吊av在线直播| 欧美一区二区精品小视频在线| 岛国在线免费视频观看| 国产激情偷乱视频一区二区| 五月伊人婷婷丁香| 成年人黄色毛片网站| 欧美+日韩+精品| 国产伦精品一区二区三区视频9| 在线免费十八禁| 国产黄片美女视频| 久久久久久久午夜电影| 九九爱精品视频在线观看| 国产视频一区二区在线看| 舔av片在线| 国产 一区 欧美 日韩| 久久国产乱子免费精品| 国产淫片久久久久久久久| .国产精品久久| 国产精品女同一区二区软件 | 99久久中文字幕三级久久日本| 亚洲成a人片在线一区二区| 亚洲中文字幕日韩| 日本a在线网址| 99热精品在线国产| 国产伦精品一区二区三区视频9| 无人区码免费观看不卡| 欧美区成人在线视频| 嫩草影院精品99| 亚洲天堂国产精品一区在线| 亚洲乱码一区二区免费版| 亚洲最大成人手机在线| 天天一区二区日本电影三级| 亚洲av日韩精品久久久久久密| 国产av一区在线观看免费| 亚洲精品在线观看二区| 久久国内精品自在自线图片| 一个人看的www免费观看视频| 嫩草影院精品99| 淫妇啪啪啪对白视频| 国产三级在线视频| 欧美色视频一区免费| 白带黄色成豆腐渣| 亚洲最大成人av| 乱码一卡2卡4卡精品| 午夜福利18| 国产真实乱freesex| 变态另类成人亚洲欧美熟女| 国产精品精品国产色婷婷| 亚洲在线观看片| 欧美日本亚洲视频在线播放| 1024手机看黄色片| 一级黄色大片毛片| 日韩欧美免费精品| 免费一级毛片在线播放高清视频| 国产高清激情床上av| 干丝袜人妻中文字幕| 此物有八面人人有两片| 午夜精品在线福利| 波野结衣二区三区在线| 天美传媒精品一区二区| 九色成人免费人妻av| 日韩亚洲欧美综合| 亚洲经典国产精华液单| 极品教师在线免费播放| 国产色婷婷99| 丰满的人妻完整版| 亚洲精品456在线播放app | 国内精品宾馆在线| 成年版毛片免费区| av在线天堂中文字幕| 精品久久久噜噜| 99久久精品热视频| 老司机午夜福利在线观看视频| 一区二区三区激情视频| 日韩av在线大香蕉| 内射极品少妇av片p| 欧美国产日韩亚洲一区| 国产三级在线视频| 简卡轻食公司| 亚洲电影在线观看av| 91久久精品国产一区二区成人| 欧美高清成人免费视频www| 亚洲国产欧洲综合997久久,| 亚洲精品乱码久久久v下载方式| 欧美日本亚洲视频在线播放| 无人区码免费观看不卡| 黄色配什么色好看| 亚洲av.av天堂| 97热精品久久久久久| 老师上课跳d突然被开到最大视频| 搡老熟女国产l中国老女人| 欧美日韩黄片免| 香蕉av资源在线| 三级毛片av免费| av在线老鸭窝| 亚洲在线自拍视频| 男人舔奶头视频| 永久网站在线| 天美传媒精品一区二区| 99热这里只有精品一区| 一区福利在线观看| 性欧美人与动物交配| 亚洲专区中文字幕在线| 99在线视频只有这里精品首页| 麻豆成人av在线观看| 亚洲精华国产精华精| 网址你懂的国产日韩在线| 亚洲av熟女| 欧美绝顶高潮抽搐喷水| 天美传媒精品一区二区| 国产探花极品一区二区| 制服丝袜大香蕉在线| 在线国产一区二区在线| 久久久久久大精品| 久久久久久久久中文| 国产精品日韩av在线免费观看| 人妻丰满熟妇av一区二区三区| 亚洲18禁久久av| 成人综合一区亚洲| 国产成年人精品一区二区| 精品人妻一区二区三区麻豆 | 在线a可以看的网站| 在线观看午夜福利视频| 变态另类成人亚洲欧美熟女| 一边摸一边抽搐一进一小说| av天堂中文字幕网| 69av精品久久久久久| 亚洲最大成人中文| 69人妻影院| 国产探花极品一区二区| 亚洲av一区综合| 欧美+日韩+精品| 国产黄色小视频在线观看| 搞女人的毛片| 国产伦人伦偷精品视频| 国产免费男女视频| 欧美一区二区亚洲| 又黄又爽又刺激的免费视频.| 高清日韩中文字幕在线| 欧美高清性xxxxhd video| 在线观看一区二区三区| 欧美另类亚洲清纯唯美| 18禁黄网站禁片免费观看直播| 国产单亲对白刺激| 亚洲国产精品sss在线观看| 久久久国产成人精品二区| 日韩欧美精品v在线| 久久久国产成人精品二区| 国产精品福利在线免费观看| 麻豆精品久久久久久蜜桃| 欧美+亚洲+日韩+国产| 国产真实乱freesex| 国产精品综合久久久久久久免费| 精品一区二区三区视频在线| 日韩欧美 国产精品| 能在线免费观看的黄片| 成人国产一区最新在线观看| 国产人妻一区二区三区在| 久久热精品热| 国内精品久久久久久久电影| av中文乱码字幕在线| 成人鲁丝片一二三区免费| 亚洲国产精品sss在线观看| 女人被狂操c到高潮| 村上凉子中文字幕在线| xxxwww97欧美| 精品久久久久久,| 最新在线观看一区二区三区| 美女免费视频网站| 欧美日韩亚洲国产一区二区在线观看| 日韩精品有码人妻一区| 日本三级黄在线观看| 日韩欧美三级三区| 国内精品久久久久精免费| 亚洲精品色激情综合| 精品国内亚洲2022精品成人| 国产精品三级大全| 韩国av一区二区三区四区| 最近中文字幕高清免费大全6 | 国产av在哪里看| 97超级碰碰碰精品色视频在线观看| 亚洲精品在线观看二区| 精品无人区乱码1区二区| 国产精品综合久久久久久久免费| 久99久视频精品免费| 精品人妻1区二区| 国产单亲对白刺激| 亚洲美女视频黄频| 老师上课跳d突然被开到最大视频| 久久中文看片网| 国产成人a区在线观看| 国产黄a三级三级三级人| 精品国内亚洲2022精品成人| 成人特级av手机在线观看| 麻豆一二三区av精品| 国国产精品蜜臀av免费| 中文资源天堂在线| 成人午夜高清在线视频| 一级黄色大片毛片| 亚洲va在线va天堂va国产| 国产高潮美女av| 久久99热6这里只有精品| 婷婷色综合大香蕉| 欧美中文日本在线观看视频| 最新中文字幕久久久久| 国产又黄又爽又无遮挡在线| 欧美黑人欧美精品刺激| 亚洲中文日韩欧美视频| 狠狠狠狠99中文字幕| 亚洲av二区三区四区| 禁无遮挡网站| 在线观看66精品国产| 在现免费观看毛片| 在线观看av片永久免费下载| 国产视频内射| 搡老岳熟女国产| 99久久成人亚洲精品观看| 久久精品影院6| 在线国产一区二区在线| 欧美性猛交╳xxx乱大交人| 色综合色国产| 日韩欧美精品免费久久| 国产精品一区二区性色av| bbb黄色大片| av国产免费在线观看| 成人亚洲精品av一区二区| 在线免费观看的www视频| 免费人成在线观看视频色| 偷拍熟女少妇极品色| 天天一区二区日本电影三级| 欧美xxxx黑人xx丫x性爽|