• <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
    免费观看的影片在线观看| 永久网站在线| 久久精品人妻少妇| 国产高清视频在线观看网站| 高清毛片免费观看视频网站| 老司机午夜福利在线观看视频| 国产精品野战在线观看| 在线观看免费视频日本深夜| 午夜a级毛片| 制服丝袜大香蕉在线| 三级毛片av免费| 国内精品久久久久精免费| 一级黄色大片毛片| 在线观看一区二区三区| 超碰av人人做人人爽久久| 99精品久久久久人妻精品| 国产aⅴ精品一区二区三区波| 老熟妇乱子伦视频在线观看| 91av网一区二区| 哪里可以看免费的av片| 给我免费播放毛片高清在线观看| 最近最新中文字幕大全电影3| 欧美不卡视频在线免费观看| 美女大奶头视频| 乱人视频在线观看| 日本撒尿小便嘘嘘汇集6| 国产av麻豆久久久久久久| av.在线天堂| 国产私拍福利视频在线观看| 性色avwww在线观看| 亚洲欧美精品综合久久99| 国产精品不卡视频一区二区| 日本五十路高清| 国产一区二区三区在线臀色熟女| 国产成人aa在线观看| 一本精品99久久精品77| 全区人妻精品视频| 美女高潮喷水抽搐中文字幕| 给我免费播放毛片高清在线观看| 中文字幕久久专区| 身体一侧抽搐| 哪里可以看免费的av片| 51国产日韩欧美| 啦啦啦观看免费观看视频高清| 免费观看在线日韩| 精品乱码久久久久久99久播| 一本久久中文字幕| 又黄又爽又刺激的免费视频.| 精品不卡国产一区二区三区| 我的女老师完整版在线观看| 午夜免费激情av| av黄色大香蕉| 国产熟女欧美一区二区| 亚洲av.av天堂| 国产精品久久久久久av不卡| 欧美日韩亚洲国产一区二区在线观看| 老熟妇乱子伦视频在线观看| 亚洲三级黄色毛片| 我的老师免费观看完整版| 亚洲精品色激情综合| 一本久久中文字幕| 亚洲人成网站在线播| 国产探花极品一区二区| 亚洲精品成人久久久久久| 欧美性猛交╳xxx乱大交人| 三级毛片av免费| 色综合婷婷激情| 97碰自拍视频| 听说在线观看完整版免费高清| 乱系列少妇在线播放| 麻豆一二三区av精品| 亚洲成人久久性| 国产av不卡久久| 99在线视频只有这里精品首页| 国产在线精品亚洲第一网站| av在线蜜桃| 国产亚洲91精品色在线| 久久午夜亚洲精品久久| 免费观看精品视频网站| 波多野结衣高清无吗| 最近视频中文字幕2019在线8| x7x7x7水蜜桃| 老司机深夜福利视频在线观看| 精品人妻视频免费看| 搡老妇女老女人老熟妇| 亚洲av不卡在线观看| 久久久久国产精品人妻aⅴ院| 久久精品国产自在天天线| 一个人观看的视频www高清免费观看| 欧美极品一区二区三区四区| 日本成人三级电影网站| 亚洲国产欧洲综合997久久,| 一级黄片播放器| 国产精品不卡视频一区二区| 欧美bdsm另类| 男插女下体视频免费在线播放| 黄色欧美视频在线观看| 国产三级中文精品| 两个人视频免费观看高清| 亚洲av电影不卡..在线观看| 黄色日韩在线| 日本精品一区二区三区蜜桃| 国产精品美女特级片免费视频播放器| 搡老熟女国产l中国老女人| 深夜a级毛片| 啪啪无遮挡十八禁网站| 日本一本二区三区精品| 在线观看av片永久免费下载| 一级毛片久久久久久久久女| 内地一区二区视频在线| 一区二区三区激情视频| 久久精品综合一区二区三区| 国产v大片淫在线免费观看| 国产毛片a区久久久久| 一区福利在线观看| 嫁个100分男人电影在线观看| 最新在线观看一区二区三区| 看黄色毛片网站| 国产免费一级a男人的天堂| 久久人人精品亚洲av| 日本爱情动作片www.在线观看 | 久久久久久久久久成人| 精品国产三级普通话版| 亚洲四区av| 国产精品野战在线观看| 亚洲精品成人久久久久久| 久久午夜福利片| 观看免费一级毛片| 久久精品人妻少妇| 欧美在线一区亚洲| 久久久精品大字幕| 国产男靠女视频免费网站| 日本免费a在线| 亚洲人成网站在线播| 成人特级黄色片久久久久久久| 天堂动漫精品| 国产日本99.免费观看| 一个人看视频在线观看www免费| 一a级毛片在线观看| 18禁裸乳无遮挡免费网站照片| 成人一区二区视频在线观看| 日韩欧美精品v在线| 欧美成人a在线观看| 久久久精品大字幕| 亚洲欧美日韩无卡精品| 国产成人影院久久av| 联通29元200g的流量卡| 日韩欧美精品v在线| 国产精品久久视频播放| 免费一级毛片在线播放高清视频| 97人妻精品一区二区三区麻豆| 天天躁日日操中文字幕| 淫妇啪啪啪对白视频| 又粗又爽又猛毛片免费看| 在线观看美女被高潮喷水网站| 亚洲图色成人| 91午夜精品亚洲一区二区三区 | 亚洲电影在线观看av| 亚洲av日韩精品久久久久久密| 国产v大片淫在线免费观看| 国产视频一区二区在线看| 成年女人永久免费观看视频| 欧美成人免费av一区二区三区| 国产高清不卡午夜福利| 亚洲av日韩精品久久久久久密| 91麻豆av在线| 夜夜看夜夜爽夜夜摸| 噜噜噜噜噜久久久久久91| 婷婷丁香在线五月| 国产一级毛片七仙女欲春2| 国产一区二区三区在线臀色熟女| 99久久成人亚洲精品观看| 国产精品一区二区性色av| 嫩草影院精品99| 亚洲avbb在线观看| 91午夜精品亚洲一区二区三区 | av在线蜜桃| 一级毛片久久久久久久久女| 99热这里只有是精品50| 村上凉子中文字幕在线| 禁无遮挡网站| 亚洲欧美日韩东京热| av免费在线看不卡| 天堂俺去俺来也www色官网| 成人亚洲精品一区在线观看 | 极品少妇高潮喷水抽搐| 午夜福利在线在线| 美女cb高潮喷水在线观看| 另类亚洲欧美激情| 亚洲国产精品专区欧美| 欧美xxxx性猛交bbbb| 九草在线视频观看| 国产一区亚洲一区在线观看| 2021少妇久久久久久久久久久| av专区在线播放| 身体一侧抽搐| 欧美另类一区| 极品少妇高潮喷水抽搐| 只有这里有精品99| videossex国产| 国产精品麻豆人妻色哟哟久久| 亚洲激情五月婷婷啪啪| 亚洲欧美清纯卡通| 久久精品国产亚洲av涩爱| 国内揄拍国产精品人妻在线| 狂野欧美激情性bbbbbb| 夫妻性生交免费视频一级片| 在线观看国产h片| 涩涩av久久男人的天堂| 亚洲精品久久久久久婷婷小说| 美女内射精品一级片tv| 欧美老熟妇乱子伦牲交| 97超视频在线观看视频| 免费看av在线观看网站| 日本-黄色视频高清免费观看| 免费大片18禁| 国产免费又黄又爽又色| 美女主播在线视频| 最近2019中文字幕mv第一页| 成人亚洲精品一区在线观看 | 美女中出高潮动态图| 最黄视频免费看| 亚洲欧美中文字幕日韩二区| 人体艺术视频欧美日本| 校园人妻丝袜中文字幕| 一级毛片aaaaaa免费看小| 婷婷色综合大香蕉| 亚洲国产色片| 国产熟女欧美一区二区| 一二三四中文在线观看免费高清| 久热这里只有精品99| 又大又黄又爽视频免费| 日本黄色片子视频| 身体一侧抽搐| 丰满少妇做爰视频| 伊人久久精品亚洲午夜| 日本猛色少妇xxxxx猛交久久| 欧美高清性xxxxhd video| 天天躁夜夜躁狠狠久久av| 午夜福利视频精品| 亚洲人成网站在线观看播放| 亚洲伊人久久精品综合| 国产精品国产三级国产专区5o| 欧美日韩一区二区视频在线观看视频在线| 91久久精品国产一区二区成人| 成年女人在线观看亚洲视频| 日本vs欧美在线观看视频 | 国产毛片在线视频| 又粗又硬又长又爽又黄的视频| 亚洲成人一二三区av| 大话2 男鬼变身卡| 精品久久久久久久末码| 久久精品国产亚洲av涩爱| 国产男女内射视频| 精品亚洲成a人片在线观看 | 免费黄频网站在线观看国产| 日韩伦理黄色片| 国产中年淑女户外野战色| 国产熟女欧美一区二区| 免费观看性生交大片5| h视频一区二区三区| av免费观看日本| 在线 av 中文字幕| 国产精品一区二区在线观看99| 成人漫画全彩无遮挡| 国产 一区 欧美 日韩| 亚洲精品亚洲一区二区| 久久久色成人| 国产av精品麻豆| 最近中文字幕高清免费大全6| 亚洲av中文字字幕乱码综合| 99久久中文字幕三级久久日本| 18禁动态无遮挡网站| 午夜福利网站1000一区二区三区| 高清不卡的av网站| 91精品一卡2卡3卡4卡| 久久人人爽人人爽人人片va| 中文字幕亚洲精品专区| 亚洲欧美精品专区久久| 国产精品国产三级国产专区5o| 精品熟女少妇av免费看| 18禁在线无遮挡免费观看视频| 尾随美女入室| 97在线视频观看| 又黄又爽又刺激的免费视频.| 日韩伦理黄色片| 久久99热6这里只有精品| 交换朋友夫妻互换小说| 2021少妇久久久久久久久久久| 国产黄片美女视频| 91久久精品国产一区二区成人| 熟女av电影| 各种免费的搞黄视频| 中文字幕免费在线视频6| 亚洲欧美精品专区久久| 久久99热6这里只有精品| 亚洲av欧美aⅴ国产| 亚洲国产最新在线播放| 亚洲成人av在线免费| 亚洲精品久久午夜乱码| 精品少妇黑人巨大在线播放| 国产综合精华液| 国产欧美日韩精品一区二区| 一本—道久久a久久精品蜜桃钙片| 免费黄网站久久成人精品| 一区二区三区精品91| 久久久久精品性色| 高清黄色对白视频在线免费看 | 国产 一区 欧美 日韩| 国产 一区精品| 少妇人妻精品综合一区二区| 久久热精品热| 一本久久精品| 免费黄频网站在线观看国产| 国产高潮美女av| 国产精品人妻久久久影院| 最近2019中文字幕mv第一页| 欧美激情国产日韩精品一区| 欧美日韩视频高清一区二区三区二| 少妇的逼水好多| 777米奇影视久久| 在线观看人妻少妇| 亚洲精品乱久久久久久| av在线老鸭窝| 亚州av有码| 蜜桃久久精品国产亚洲av| 成人黄色视频免费在线看| 国产中年淑女户外野战色| av网站免费在线观看视频| 最近手机中文字幕大全| 成年美女黄网站色视频大全免费 | 在线观看美女被高潮喷水网站| 免费看不卡的av| 少妇的逼好多水| 国产免费又黄又爽又色| 国产精品av视频在线免费观看| 97精品久久久久久久久久精品| 麻豆乱淫一区二区| 97在线人人人人妻| 国产av精品麻豆| 美女xxoo啪啪120秒动态图| 亚洲人成网站在线播| 国产高潮美女av| 18禁在线无遮挡免费观看视频| 午夜视频国产福利| 亚洲精品日韩av片在线观看| 国产在视频线精品| 国产色婷婷99| 只有这里有精品99| 美女内射精品一级片tv| 人人妻人人看人人澡| 少妇人妻一区二区三区视频| 亚洲精品久久午夜乱码| 99热这里只有精品一区| av福利片在线观看| 18禁在线播放成人免费| 亚洲精品乱码久久久v下载方式| 亚洲激情五月婷婷啪啪| 欧美性感艳星| 国产又色又爽无遮挡免| 亚洲一级一片aⅴ在线观看| 青春草国产在线视频| 成人亚洲精品一区在线观看 | 干丝袜人妻中文字幕| tube8黄色片| 99热国产这里只有精品6| 黄色视频在线播放观看不卡| 看十八女毛片水多多多| 久久ye,这里只有精品| 久久久a久久爽久久v久久| 中文欧美无线码| 高清欧美精品videossex| 亚洲av二区三区四区| 国产亚洲欧美精品永久| 高清毛片免费看| 美女脱内裤让男人舔精品视频| 久久影院123| 99热这里只有精品一区| 日本午夜av视频| 国产成人aa在线观看| 波野结衣二区三区在线| 一边亲一边摸免费视频| 国产久久久一区二区三区| 免费观看av网站的网址| 精品久久久久久久久av| 国产高清三级在线| 国内揄拍国产精品人妻在线| 一级黄片播放器| 大话2 男鬼变身卡| 久久 成人 亚洲| 国产精品国产av在线观看| 男女边摸边吃奶| 日韩 亚洲 欧美在线| 国产av码专区亚洲av| 欧美日韩视频高清一区二区三区二| 联通29元200g的流量卡| 亚洲综合色惰| 国产精品久久久久久精品古装| 大话2 男鬼变身卡| 午夜老司机福利剧场| 欧美最新免费一区二区三区| 黑人猛操日本美女一级片| 97在线人人人人妻| 3wmmmm亚洲av在线观看| 国产老妇伦熟女老妇高清| 五月开心婷婷网| 一二三四中文在线观看免费高清| 欧美日本视频| 久久久久久伊人网av| 欧美成人精品欧美一级黄| 日韩欧美 国产精品| 一区二区三区乱码不卡18| 国产av精品麻豆| 免费观看性生交大片5| 夜夜爽夜夜爽视频| 国产伦精品一区二区三区视频9| 国产成人aa在线观看| 欧美三级亚洲精品| av国产久精品久网站免费入址| 国产精品一区二区在线观看99| 国产精品欧美亚洲77777| 精品亚洲乱码少妇综合久久| 亚洲欧美精品专区久久| 老司机影院成人| 久久久久久久久久人人人人人人| a级毛片免费高清观看在线播放| 国产成人a区在线观看| 男的添女的下面高潮视频| 精品久久久精品久久久| 91精品一卡2卡3卡4卡| 高清av免费在线| 国产精品麻豆人妻色哟哟久久| 国产一区二区三区av在线| 最近的中文字幕免费完整| 插逼视频在线观看| 国模一区二区三区四区视频| 中文资源天堂在线| 久久人人爽人人爽人人片va| tube8黄色片| 少妇人妻 视频| 久久精品久久久久久久性| 欧美区成人在线视频| 久久ye,这里只有精品| 一个人免费看片子| 国产一区亚洲一区在线观看| 国产又色又爽无遮挡免| 一区二区三区四区激情视频| 国产有黄有色有爽视频| 亚洲人与动物交配视频| av国产免费在线观看| 老熟女久久久| 麻豆国产97在线/欧美| 日日摸夜夜添夜夜添av毛片| 午夜福利视频精品| 两个人的视频大全免费| 天堂俺去俺来也www色官网| 在线观看一区二区三区| 久久国产精品男人的天堂亚洲 | a级一级毛片免费在线观看| 成人毛片60女人毛片免费| 亚洲国产精品成人久久小说| 日日啪夜夜撸| 成年av动漫网址| av免费在线看不卡| 国产精品国产av在线观看| 最近中文字幕高清免费大全6| 亚洲成人av在线免费| 国产精品久久久久久精品电影小说 | av在线播放精品| 免费av不卡在线播放| 男人爽女人下面视频在线观看| freevideosex欧美| 国产深夜福利视频在线观看| 日本猛色少妇xxxxx猛交久久| 麻豆成人av视频| 卡戴珊不雅视频在线播放| 夜夜看夜夜爽夜夜摸| 亚洲精品第二区| 26uuu在线亚洲综合色| 国产精品麻豆人妻色哟哟久久| 亚洲精品国产色婷婷电影| 国产成人精品福利久久| 久久久久网色| 街头女战士在线观看网站| 亚洲国产精品国产精品| 欧美少妇被猛烈插入视频| 日韩,欧美,国产一区二区三区| 国模一区二区三区四区视频| 日韩av在线免费看完整版不卡| 免费看日本二区| 成人国产av品久久久| 中文字幕免费在线视频6| 国产一区亚洲一区在线观看| 人妻一区二区av| videossex国产| 日韩人妻高清精品专区| 成年av动漫网址| 高清av免费在线| 日韩国内少妇激情av| 一级毛片我不卡| 男人舔奶头视频| 高清毛片免费看| 日本猛色少妇xxxxx猛交久久| 日韩国内少妇激情av| 免费观看的影片在线观看| 久久99精品国语久久久| 中文在线观看免费www的网站| 伦理电影免费视频| 女性生殖器流出的白浆| 肉色欧美久久久久久久蜜桃| 国产精品99久久99久久久不卡 | 日本av免费视频播放| 乱码一卡2卡4卡精品| 噜噜噜噜噜久久久久久91| 黄色日韩在线| 免费黄色在线免费观看| 日韩制服骚丝袜av| a级一级毛片免费在线观看| 91久久精品国产一区二区成人| 国产免费一级a男人的天堂| 蜜桃亚洲精品一区二区三区| 一本—道久久a久久精品蜜桃钙片| 欧美老熟妇乱子伦牲交| 亚洲va在线va天堂va国产| 亚洲真实伦在线观看| 亚洲不卡免费看| 有码 亚洲区| 女人久久www免费人成看片| 丝袜脚勾引网站| 最黄视频免费看| 联通29元200g的流量卡| 国产成人免费无遮挡视频| 亚洲av欧美aⅴ国产| 狠狠精品人妻久久久久久综合| 九色成人免费人妻av| 大香蕉久久网| 下体分泌物呈黄色| 久久精品久久久久久噜噜老黄| 秋霞伦理黄片| 亚洲成人手机| 亚洲美女搞黄在线观看| 成人影院久久| 在线免费十八禁| av在线观看视频网站免费| 日本黄色片子视频| 成年美女黄网站色视频大全免费 | 亚洲av不卡在线观看| 久久久久久久久大av| 日韩视频在线欧美| 亚洲色图综合在线观看| 色网站视频免费| 国产黄色视频一区二区在线观看| 日韩伦理黄色片| 亚洲av成人精品一区久久| 深爱激情五月婷婷| 最近中文字幕2019免费版| 精品一区二区三区视频在线| 亚洲熟女精品中文字幕| 国产免费又黄又爽又色| 久久青草综合色| 午夜福利在线观看免费完整高清在| 欧美丝袜亚洲另类| www.色视频.com| 18+在线观看网站| videossex国产| 亚洲精品久久久久久婷婷小说| 伊人久久国产一区二区| av一本久久久久| 国产有黄有色有爽视频| 国产69精品久久久久777片| 中国三级夫妇交换| 99re6热这里在线精品视频| 亚洲va在线va天堂va国产| 国产精品偷伦视频观看了| av在线老鸭窝| 最近中文字幕高清免费大全6| 国产精品一区二区性色av| 18禁动态无遮挡网站| 国产男女内射视频| 国产黄片美女视频| 一级黄片播放器| 99久久中文字幕三级久久日本| 日日摸夜夜添夜夜添av毛片| 美女cb高潮喷水在线观看| 亚洲欧美日韩无卡精品| 中文字幕久久专区| 在线观看免费高清a一片| 男女边摸边吃奶| av在线app专区| 一级毛片久久久久久久久女| 亚洲欧美一区二区三区黑人 | 国产白丝娇喘喷水9色精品| 久久精品夜色国产| 99热网站在线观看| 伦精品一区二区三区| 五月天丁香电影| 亚洲怡红院男人天堂| 国产精品不卡视频一区二区| 中文字幕精品免费在线观看视频 | 成人综合一区亚洲| 91久久精品国产一区二区成人| 久久精品人妻少妇| 国产日韩欧美亚洲二区| 观看美女的网站| 91久久精品电影网| 直男gayav资源| 免费大片18禁| 国产女主播在线喷水免费视频网站| 免费黄网站久久成人精品| 激情 狠狠 欧美| 欧美xxⅹ黑人| 99热这里只有精品一区| 搡老乐熟女国产| 永久免费av网站大全| 直男gayav资源|