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

    茶樹NUDIX基因家族的鑒定及表達分析

    2023-05-08 06:49:20陳振艷張相琴陳蘭謝思藝劉碩謙田娜
    茶葉科學(xué) 2023年2期
    關(guān)鍵詞:基序內(nèi)含子逆境

    陳振艷,張相琴,陳蘭,謝思藝,劉碩謙,2*,田娜,2*

    茶樹基因家族的鑒定及表達分析

    陳振艷1,張相琴1,陳蘭1,謝思藝1,劉碩謙1,2*,田娜1,2*

    1. 湖南農(nóng)業(yè)大學(xué)園藝學(xué)院茶學(xué)教育部重點實驗室,湖南 長沙 410128;2. 湖南農(nóng)業(yè)大學(xué)園藝學(xué)院茶學(xué)系,湖南 長沙 410128

    NUDIX水解酶屬于焦磷酸酶,在信息傳遞、植物生長協(xié)調(diào)和響應(yīng)逆境脅迫等方面起著重要作用?;诓铇洌ǎ┗蚪M鑒定出43個家族基因,并運用生物信息學(xué)和實時熒光定量PCR等方法對其分析。結(jié)果表明,43個CsNUDXs蛋白質(zhì)分子量介于11.8~89.2?kDa,氨基酸長度為102~342個,理論等電點為4.49~9.26,18.6%的蛋白為穩(wěn)定性蛋白;根據(jù)進化關(guān)系將分為6個亞家族;啟動子順式作用元件分析發(fā)現(xiàn)具有許多與激素響應(yīng)、逆境脅迫和生長發(fā)育等調(diào)控相關(guān)的作用元件;對家族基因在不同器官表達模式進行分析,發(fā)現(xiàn)、在果實中表達量較高,而、在果實中表達量極低,在老葉中表達量極低;不同處理組的實時熒光定量PCR結(jié)果表明,1?mmol·L-1茉莉酸甲酯(MeJA)處理下、和等表達量呈現(xiàn)先升后降再升的趨勢,而1?mmol·L-1水楊酸(SA)處理下、和等表達量則呈現(xiàn)先降后升再降的趨勢,300?mmol·L-1NaCl處理下和等表達量呈現(xiàn)先升后降的趨勢。本研究應(yīng)用生物信息學(xué)技術(shù)初步解析了的基本特征和功能,發(fā)現(xiàn)可響應(yīng)高鹽脅迫和MeJA、SA處理。

    茶樹;基因家族;非生物脅迫;表達分析

    NUDIX水解酶(Nucleoside diphosphate linked to x hydrolases)是一類能夠催化水解核苷糖、二核苷多聚磷酸鹽和三磷酸核苷等多類有機焦磷酸鹽的焦磷酸酶[1],廣泛分布在生物體內(nèi)。大腸桿菌中MutT(NudA)是最早被發(fā)現(xiàn)的NUDIX水解酶,故NUDIX水解酶家族又被稱為MutT家族[2]。該家族特征是具有一段高度保守的氨基酸序列——GX5EX7REUXEEXGU(其中U表示疏水基團,X表示任意氨基酸),當其核心序列REUXEEX中的谷氨酰胺殘基與二價陽離子結(jié)合時,NUDIX蛋白發(fā)揮作用[2-3]。

    NUDIX水解酶在植物逆境脅迫響應(yīng)過程中發(fā)揮著重要作用。干旱脅迫下,過表達山桃()可提高、和等抗性基因的表達,有效減少煙草的干旱損傷[4];敲除擬南芥()中的基因會降低其抗氧化脅迫能力[5],當被沉默后會誘導(dǎo)防御相關(guān)基因的表達,并提高植物對丁香假胞菌的抗性[6];過表達可維持植物體內(nèi)正常的煙酰胺腺嘌呤二核苷酸(Nicotinamide adenine dinucleotide,NAD)和腺苷三磷酸(Adenosine triphosphate,ATP)含量,增強植物抗氧化脅迫的能力[5,7]。此外,NUDIX水解酶也廣泛參與了植物激素信號的轉(zhuǎn)導(dǎo)過程。對山桃苗進行脫落酸(Abscisicacid,ABA)處理后,會受到明顯誘導(dǎo),且轉(zhuǎn)基因煙草中的和等ABA信號通路中關(guān)鍵基因的轉(zhuǎn)錄水平明顯高于野生型(WT),而的轉(zhuǎn)錄水平明顯低于WT,表明可調(diào)控ABA信號轉(zhuǎn)導(dǎo)通路[4];是NPR1依水楊酸(Salicylic acid,SA)信號通路的陽性調(diào)節(jié)因子,敲除或過表達會使參與SA誘導(dǎo)NPR1激活過程中相關(guān)基因和的表達量增加或降低[4,8];與野生型相比,和基因共缺失型擬南芥的差異基因?qū)岳蛩峒柞ィ∕ethyl jasmonate,MeJA)的響應(yīng)顯著[9]。目前,家族已經(jīng)在擬南芥[10]、陸地棉()[11]、桃樹()[4]、玫瑰()[12]、甘菊()[13]、大豆()[14]和葡萄()[15]等植物中被鑒定,其在擬南芥中響應(yīng)植物激素和逆境脅迫的作用機制已逐漸被揭開,但在茶樹中的研究尚淺[16]。

    因此本研究基于茶樹基因組數(shù)據(jù)庫,對茶樹家族基因進行全基因組鑒定,并分析茶樹基因家族成員在不同器官及高鹽脅迫和激素處理(MeJA和SA)下的表達模式,為進一步探究家族基因在茶樹中的作用奠定基礎(chǔ)。

    1 材料與方法

    1.1 試驗材料

    以碧香早茶樹枝梢為試驗材料,采自湖南農(nóng)業(yè)大學(xué)茶葉研究所茶園(113° E,28° N)。

    高鹽處理:將長勢一致且生長狀況良好的茶樹枝梢插入300?mmol·L-1NaCl溶液中。激素處理:向茶樹枝梢分別噴灑1?mmol·L-1MeJA溶液和1?mmol·L-1SA溶液,以液體布滿葉片為標準;并將插入ddH2O中的茶樹枝梢作為對照。將以上處理材料放入溫度為(25±2)℃、相對濕度為70%~75%、光照周期為14?h光/10?h暗的氣候室中進行培養(yǎng)。在處理0、3、9、24?h和48?h后隨機摘取3~5片嫩葉,取樣后進行液氮冷凍固樣,置于–80?℃保存?zhèn)溆?,每個處理3次生物學(xué)重復(fù)。

    1.2 試驗方法

    1.2.1基因家族的篩選

    從茶樹基因組數(shù)據(jù)庫(http://teaplant.org/index.html)中下載茶樹的全基因組數(shù)據(jù),從擬南芥基因組數(shù)據(jù)庫(https://www. arabidopsis.org)中下載擬南芥NUDIX蛋白序列,用于后續(xù)分析。以擬南芥的25個NUDIX蛋白序列為參考,在茶樹基因組數(shù)據(jù)庫中進行BlastP比對,設(shè)定閾值E<1E-6,獲得潛在的基因家族序列;從PFAM(http://pfam. xfam.org)數(shù)據(jù)庫中下載NUDIX基因特有的保守結(jié)構(gòu)域模型(PF00293),并通過HMMER 3.0軟件中的hmmersearch搜索鑒定含有NUDIX保守結(jié)構(gòu)域的茶樹蛋白序列。合并已獲得的兩個推定的蛋白序列文件,提交至SMART(http://smart.embl-heidelberg.de)、PFAM和NCBI-CDD(http://www.ncbi.nlm.nih.gov)等網(wǎng)站進行驗證,整理并剔除冗余序列,獲得最終的茶樹NUDIX基因家族成員。

    1.2.2 CsNUDXs蛋白理化性質(zhì)及亞細胞定位預(yù)測

    利用ExPASy Compute pI/Mw(http://web. expasy.org/compute_pi)對NUDIX蛋白的理化性質(zhì)進行分析,用WoLF PSORT(https://wolfpsort. hgc.jp)工具預(yù)測蛋白質(zhì)亞細胞定位。

    1.2.3系統(tǒng)發(fā)育關(guān)系、基因結(jié)構(gòu)和基序分析

    通過MEGA 7.0軟件對茶樹和擬南芥的NUDIX蛋白序列進行多序列比對,并采用鄰近連接法(Neighbour-joining method)構(gòu)建系統(tǒng)進化樹,校驗參數(shù)設(shè)置為1?000次[17]。將家族基因的編碼序列(Coding sequence,CDS)提交至MEME網(wǎng)站(https:// meme-suite.org)進行基因基序分析,基序數(shù)目設(shè)置為10個,氨基酸長度設(shè)置為6~55個,其余參數(shù)為默認。通過TBtools[18]軟件進行可視化分析。

    1.2.4染色體定位和啟動子順式元件分析

    利用HMMER 3.0軟件中的perl腳本從茶樹基因組注釋信息中獲取在染色體上的起始位置,并通過MG2C(http://mg2c. iask.in/mg2c_v2.1)進行可視化。運用TBtools軟件提取上游2?000?bp的序列,通過PlantCare(http://bioinformatics. psb. ugent. be/webtools/plantcare)網(wǎng)站進行啟動子順式作用元件預(yù)測,并通過TBtools[18]軟件進行可視化。

    1.2.5表達模式分析

    從茶樹基因組數(shù)據(jù)庫中下載茶樹8個不同器官(頂芽、嫩葉、成熟葉、老葉、莖、根、花和果實)中的表達量數(shù)據(jù),篩選高表達的家族成員,并通過TBtools[18]軟件進行可視化。

    1.2.6 實時熒光定量PCR分析

    利用FastPure?Universal Plant Total RNA Isolation Kit(諾唯贊,南京)分別提取對照組和3個處理組(NaCl、MeJA和SA)的茶樹嫩葉總RNA,并通過超微量紫外分光光度計(賽默飛,美國)和1.0%瓊脂糖凝膠電泳檢測其完整性、純度和濃度。使用Fastking gDNA Dispelling RT SuperMix Kit(天根,北京)將RNA反轉(zhuǎn)錄為cDNA。

    2 結(jié)果與分析

    2.1 CsNUDXs基因家族的篩選及理化性質(zhì)

    從茶樹基因組中共鑒定到43個基因,根據(jù)其在染色體上的位置進行命名,并進一步對其編碼的蛋白質(zhì)進行理化性質(zhì)分析。結(jié)果表明,茶樹NUDIX含有氨基酸殘基為102~342個;分子量為11.8~89.2?kDa;脂肪族指數(shù)為50.67~103.67;理論等電點(pI)介于4.49~9.26,其中pI<7的CsNUDXs蛋白達到70%以上。43個CsNUDXs蛋白中僅有CsNUDX1、CsNUDX3、CsNUDX6、CsNUDX10、CsNUDX21、CsNUDX22、CsNUDX30和CsNUDX31的不穩(wěn)定指數(shù)低于40,其余均為不穩(wěn)定蛋白。除CsNUDX36和CsNUDX29外,其余CsNUDIXs蛋白均為親水性蛋白。亞細胞定位預(yù)測結(jié)果顯示(表2),主要分布在細胞核、葉綠體、線粒體中。

    2.2 CsNUDXs系統(tǒng)進化

    將茶樹和擬南芥的基因家族成員構(gòu)建系統(tǒng)發(fā)育進化樹。結(jié)果表明(圖1),茶樹基因家族可分為6個亞家族,分別含有10、3、10、8、3個和9個。其中(CSS0008608.1)和(CSS0042662.1)、(CSS0021705.1)和(CSS0015340.1)等26個基因在進化樹上分別成對出現(xiàn),說明它們進化關(guān)系較近;(CSS0000911.1)、(CSS0028900.1)、(CSS0035220.1)、(CSS0026873.1)和(CSS0026618.1)等5個與進化關(guān)系較近。

    表1 CsNUDXs表達分析的實時熒光定量PCR引物

    表2 CsNUDXs的亞細胞定位

    2.3 CsNUDXs結(jié)構(gòu)和基序分析

    對43個的保守基序進行分析,共鑒定出10個保守基序,依次命名為motif 1 ~motif 10(圖2)。所包含的基序數(shù)量介于1~6,其中motif 1出現(xiàn)最頻繁,其次是motif 10、motif 7和motif 6。均含有motif 1,motif 10為亞家族Ⅰ和亞家族Ⅵ共有;此外,motif 2、motif 5和motif 8僅存在于亞家族Ⅲ中,motif 3、motif 4和motif 9僅存在于亞家族Ⅵ中。

    對6個亞家族的的結(jié)構(gòu)分析發(fā)現(xiàn)(圖3),亞家族Ⅰ中內(nèi)含子數(shù)目為1~7,其中和內(nèi)含子數(shù)目最多,為7個;亞家族Ⅱ中的含有20個內(nèi)含子,和各有5個內(nèi)含子;亞家族Ⅲ中除(4個內(nèi)含子)和(5個內(nèi)含子)外,其余基因均有3個內(nèi)含子;亞家族Ⅳ的內(nèi)含子數(shù)目為3~5,、和均有3個內(nèi)含子;亞家族Ⅴ中和有7個內(nèi)含子,僅有4個內(nèi)含子;亞家族Ⅵ的內(nèi)含子數(shù)目為5~8。綜上所述,各亞家族中內(nèi)含子數(shù)量差異較大,但同一亞家族中內(nèi)含子個數(shù)相同的基因結(jié)構(gòu)更為相似。

    2.4 CsNUDXs染色體分布和啟動子順式作用元件分析

    染色體分布結(jié)果顯示(圖4),有39個不均勻的分布在Chr1、Chr2、Chr4、Chr5、Chr6、Chr7、Chr8、Chr10、Chr12、Chr14和Chr15號染色體上。其中Chr12號染色體有6個,Chr6和Chr8號染色體各有5個,Chr1、Chr2、Chr5和Chr7號染色體各有4個,Chr14號染色體有3個,Chr10號染色體有2個,Chr4和Chr15號染色體各有1個;、、、這4個基因未定位在染色體上。

    圖1 茶樹和擬南芥NUDIX基因家族的系統(tǒng)進化樹

    圖2 茶樹NUDIX基因的motif分析

    圖3 茶樹NUDIX基因結(jié)構(gòu)分析

    對上游2?000?bp的啟動子區(qū)域進行順式作用元件分析,結(jié)果顯示(圖5),除了含有基礎(chǔ)的作用元件外(如和轉(zhuǎn)錄相關(guān)的TATA-box、CAAT-box等),還含有與光響應(yīng)、激素響應(yīng)、逆境脅迫和生長發(fā)育調(diào)控等相關(guān)的順式作用元件。

    2.5 CsNUDXs在茶樹不同器官的表達分析

    在不同器官表達模式分析結(jié)果顯示(圖6),在頂芽中表達量較高的基因有和,嫩葉中表達量較高的基因有、和,成熟葉中表達量較高的基因有、、和,老葉中表達量較高的基因有和,莖中表達量較高的基因有、和,根中表達量較高的基因有、、和等,花中表達量較高的基因有、和,、和在果實中表達量較高;此外,和在果實中表達量極低,在老葉中表達量極低,這表明在茶樹不同器官中具有表達特異性。

    圖4 茶樹NUDIX基因在茶樹染色體上的定位

    注:A表示光響應(yīng),B表示脫落酸響應(yīng),C表示損傷響應(yīng),D表示低溫響應(yīng),E表示赤霉素響應(yīng),F(xiàn)表示胚乳特異性負表達,G表示水楊酸響應(yīng),H表示干旱反應(yīng),I表示激素響應(yīng),J表示類黃酮合成基因調(diào)控,K表示玉米醇溶蛋白代謝調(diào)節(jié),L表示種子特異性調(diào)節(jié),M表示分生組織表達,N表示逆境脅迫,O表示缺氧反應(yīng),P表示厭氧反應(yīng),Q表示茉莉酸甲酯響應(yīng)

    2.6 CsNUDXs在MeJA、SA及NaCl處理下的表達模式

    篩選出、、和等9個在茶樹不同器官中表達量較高且含有大量與激素響應(yīng)或非生物脅迫相關(guān)元件的基因做進一步的驗證。

    由圖7可知,在MeJA處理下,、、和的表達量較高,的表達量持續(xù)上升,其余基因在處理后的表達量均呈先上升后下降再上升的趨勢。

    在SA處理下(圖8),和的表達量在0~9?h下降,9~24?h上升,然后繼續(xù)呈下降趨勢;和的表達量在處理0~3?h后下降,3~24?h上升,然后呈下降趨勢;和的表達量在處理0~3?h上升,3~9?h下降,9~24?h上升,然后繼續(xù)下降趨勢,其中表達量較高。

    在NaCl處理下(圖9),和的表達量較高。其中,和的表達量在處理后0~3?h增加,之后呈下降趨勢;的表達量在處理后0~3?h和9~24?h上升,3~9?h和24?h之后呈下降趨勢;在處理后0~3?h增加,3~9?h下降,之后呈上升趨勢。

    3 討論

    NUDIX水解酶可以降解植物體內(nèi)積累過多的有害代謝物,在信息傳遞、生長協(xié)調(diào)和機體保護等方面起著重要作用[20]。本研究基于茶樹全基因組數(shù)據(jù)共鑒定出43個茶樹基因。在分析43個茶樹基因保守基序時發(fā)現(xiàn),每個基因中包含的基序數(shù)量在1~6,且進化關(guān)系較近的在基序的種類和長度上都較相似,這可能是由于同一亞家族的基因在功能選擇和長期進化下形成的[21];其中,motif 1為各共有,表明它可能承擔著重要且相似的功能;motif 10為亞家族Ⅰ和亞家族Ⅵ共有,推測亞家族Ⅰ和亞家族Ⅵ的行使的功能可能有交叉。真核生物基因組的復(fù)雜程度與內(nèi)含子數(shù)量有關(guān)[22],基因間內(nèi)含子數(shù)目差異易導(dǎo)致形成不同的可變剪切體,基因組中的可變剪切體可增強茶樹對非生物脅迫的抗逆能力[23];本研究發(fā)現(xiàn)大部分家族基因的內(nèi)含子數(shù)目介于1~20,且6個亞家族間存在差異(Ⅰ類:1~7,Ⅱ類:5~20,Ⅲ、Ⅳ類:3~5,Ⅴ類:4~7,Ⅵ類:5~8),家族基因的內(nèi)含子數(shù)目差異可能導(dǎo)致形成不同的可變剪切體進而增強茶樹的抗逆能力。在鹽脅迫環(huán)境中,過量的Na+取代植物生長必需的K+,且Na+、Cl-與蛋白質(zhì)或氨基酸結(jié)合時會降低蛋白質(zhì)或氨基酸的活性,鹽脅迫引起的滲透脅迫和營養(yǎng)不協(xié)調(diào)進一步引起氧化脅迫,使植物的光合作用和呼吸作用的電子鏈發(fā)生斷裂,最終導(dǎo)致植物死亡[24]。研究發(fā)現(xiàn),甘菊在NaCl澆灌處理后,在0~48?h的表達量先上升后下降,表明可響應(yīng)鹽脅迫,可能通過降低超氧化物歧化酶(Superoxide dismutase,SOD)和過氧化氫酶(Catalase,CAT)活性[25]、增加過氧化物酶(Peroxidase,POD)活性以提升甘菊清除活性氧的能力,進而提升了甘菊適應(yīng)鹽脅迫的能力[13]。茶樹在響應(yīng)鹽脅迫時,表達量也呈先升后降的趨勢,因此推測可能通過清除茶樹體內(nèi)的活性氧參與了茶樹響應(yīng)鹽脅迫的過程。

    注:AB表示頂芽,YL表示嫩葉,ML表示成熟葉,OL表示老葉,S表示莖,R表示根,F(xiàn)L表示花,F(xiàn)R表示果實

    注:不同大、小寫字母分別表示在0.01和0.05水平差異顯著,下同

    本研究還發(fā)現(xiàn),在啟動子區(qū)富含激素響應(yīng)相關(guān)的作用元件,表明可能通過多個順式作用元件參與茶樹對激素的響應(yīng)過程。在SA處理下呈現(xiàn)先上升后下降的表達趨勢,敲除或者過表達會使得參與到SA誘導(dǎo)NPR1激活過程中的相關(guān)基因和的表達量上升或下降,表明可正向調(diào)節(jié)SA誘導(dǎo)NPR1激活過程進而影響植物免疫反應(yīng)[8];甘菊的表達量在SA的誘導(dǎo)下也先上升后下降[13]。本研究發(fā)現(xiàn)和的表達量在SA處理后也呈現(xiàn)出先上升后下降的趨勢,這與和表達模式相同,說明和可能通過正向調(diào)節(jié)SA誘導(dǎo)NPR1激活過程參與茶樹對SA的響應(yīng)。MeJA是植物體內(nèi)重要的脂類衍生化合物,可作為信號分子和生長調(diào)節(jié)物質(zhì)參與到植物的生長發(fā)育和對逆境脅迫的響應(yīng)等過程[26]。MeJA可通過加深植物因鹽脅迫導(dǎo)致的生長抑制,降低6-芐氨基腺嘌呤(6-Benzylaminopurine,6-BA)和赤霉素(Gibberellin,GA)等植物激素的轉(zhuǎn)錄水平,使植株能更有效地進行離子穩(wěn)態(tài)、滲透調(diào)節(jié)和抗氧化防御等反應(yīng)[27]。本研究發(fā)現(xiàn),在MeJA處理下大部分的表達模式有差異,表明響應(yīng)MeJA刺激的作用機制不同。

    圖8 外源SA處理下部分CsNUDXs的相對表達量

    圖9 外源NaCl處理下部分CsNUDXs的相對表達量

    此外,NUDIX水解酶在植物生長過程中起著重要的協(xié)調(diào)作用[15],亞細胞定位預(yù)測結(jié)果顯示主要分布在細胞核、葉綠體、線粒體中,這與大豆()的結(jié)果一致[14],說明可能在細胞核、葉綠體和線粒體中發(fā)揮作用。玫瑰在花瓣中表達量最高,且花發(fā)育后期香氣散發(fā)達到最大時表達量持續(xù)增加,而抑制的表達后玫瑰芳香物質(zhì)合成減少[28]。玫瑰品種唐紅中的表達量與花朵開放程度呈正相關(guān),且主要芳香成分也隨花朵開放程度增加而增加,在盛開期達到最高;將轉(zhuǎn)入矮牽牛(Vilm)中可顯著增加矮牽牛的主要香氣成分含量,且感官上可聞到濃郁香氣,過表達也可顯著增強矮牽?;ㄏ愕臐庥舫潭萚29]。以上結(jié)果表明,家族基因的組織表達特異性與其功能具有較高的相關(guān)性。本研究發(fā)現(xiàn),和在花中表達量較高而在其他組織部位中表達量極低,推測其主要在花中發(fā)揮作用。

    本研究首次對全基因組進行了鑒定和表達模式分析,在43個的啟動子區(qū)發(fā)現(xiàn)了許多與光響應(yīng)、激素響應(yīng)和逆境脅迫相關(guān)的作用元件,且在茶樹各組織中均有表達,表明參與茶樹的生長發(fā)育調(diào)控。此外,利用實時熒光定量PCR技術(shù)檢測了9個在SA、MeJA和NaCl處理下的表達模式,發(fā)現(xiàn)其均可響應(yīng)SA、MeJA和NaCl處理。綜上所述,可以響應(yīng)多種逆境脅迫,但在協(xié)調(diào)茶樹生長發(fā)育和響應(yīng)逆境脅迫過程中的具體機制仍需進一步探清。

    [1] Mildvan A S, Xia Z, Azurmendi H F, et al. Structures and mechanisms of Nudix hydrolases [J]. Archives of Biochemistry and Biophysics, 2005, 433(1): 129-143.

    [2] Bessman M J, Frick D N, O'Handley S F. The MutT proteins or "Nudix" hydrolases, a family of versatile, widely distributed, "housecleaning" enzymes [J]. The Journal of Biological Chemistry, 1996, 271(41): 25059-25062.

    [3] McLennan A G. The Nudix hydrolase superfamily [J]. Cellular and Molecular Life Sciences, 2006, 63(2): 123-143.

    [4] 和華杰. 桃基因響應(yīng)干旱脅迫的分子機制研究[D]. 泰安: 山東農(nóng)業(yè)大學(xué), 2022. He H J. Molecular mechanism ofgene response to drought stress in peach [D]. Tai'an: Shandong Agricultural University, 2022.

    [5] Ogawa T, Yoshimura K, Miyake H, et al. Molecular characterization of organelle-type Nudix hydrolases in Arabidopsis [J]. Plant Physiology, 2008, 148(3): 1412-1424.

    [6] Dong S M, Wang Y C. Correction: Nudix effectors: a common weapon in the arsenal of plant pathogens [J]. Plos Pathogens, 2016, 12(10): e1005970. doi: 10.1371/journal. ppat.1005970.

    [7] Ogawa T, Ishikawa K, Harada K, et al. Overexpression of an ADP-ribose pyrophosphatase,, confers enhanced tolerance to oxidative stress in Arabidopsis plants [J]. The Plant Journal, 2009, 57(2): 289-301.

    [8] Ishikawa K, Yoshimura K, Harada K, et al. AtNUDX6, an ADP-ribose/NADH pyrophosphohydrolase in Arabidopsis, positively regulates NPR1-dependent salicylic acid signaling [J]. Plant Physiology, 2010, 152(4): 2000-2012.

    [9] 尹虹.基因缺失型擬南芥基因芯片表達譜的生物信息學(xué)分析[J]. 錦州醫(yī)科大學(xué)學(xué)報, 2021, 42(5): 30-35. Yin H. Bioinformatic analysis of expression profiles of gene chip of arabidopsis withgene deletion [J]. Journal of Jinzhou Medical University, 2022, 42(5): 30-35.

    [10] Ogawa T, Uedo Y, Yoshimura K, et al. Comprehensive analysis of cytosolic Nudix hydrolases in[J]. The Journal of Biological Chemistry, 2005, 280(26): 25277-25283.

    [11] 竇玲玲, 孫亞如, 趙琴, 等. 陸地棉基因家族的全基因組鑒定及表達分析[J]. 棉花學(xué)報, 2021, 33(2): 112-123. Dou L L, Sun Y R, Zhao Q, et al. Genome-wide identification and expression analysis ofgene family in upland cotton [J]. Cotton Science, 2021, 33(2): 112-123.

    [12] Bergman M E, Bhardwaj M, Phillips M A. Cytosolic geraniol and citronellol biosynthesis require a Nudix hydrolase in rose scented geranium () [J]. The Plant Journal, 2021, 107(2): 493-510.

    [13] 黃河. 甘菊響應(yīng)鹽誘導(dǎo)的分子機理研究[D]. 北京: 北京林業(yè)大學(xué), 2012. Huang H. The molecular mechanism of salt response in(Fisch. ex Trautv.) Makino [D]. Beijing: Beijing Forestry University, 2012.

    [14] ?,|, 王娟, 于洋, 等. 大豆基因家族全基因組分析[J]. 西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版), 2018, 46(9): 27-34. Chang W, Wang J, Yu Y, et al. Genome-wide analysis of soybeangene family [J]. Journal of Northwest A&F University (Natural Science Edition), 2018, 46(9): 27-34.

    [15] Wang P P, Wang Z K, Guan L, et al. Versatile physiological functions of the Nudix hydrolase family in berry development and stress response in grapevine [J]. Journal of Integrative Agriculture, 2022, 21(1): 91-112.

    [16] 龐靜怡. 茶樹核苷水解酶基因與香葉醇生物合成的研究[D]. 合肥: 安徽農(nóng)業(yè)大學(xué), 2018. Pang J Y. Studies on tea (var.) nucleoside hydrolasegene and geraniol biosynthesis [D]. Hefei: Anhui Agricultural University, 2018.

    [17] Newman L, Duffus A L J, Lee C. Using the free program MEGA to build phylogenetic trees from molecular data [J]. The American Biology Teacher, 2016, 78(7): 608-612.

    [18] Chen C J, Chen H, Zhang Y, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data [J]. Molecular Plant, 2020, 13(8): 1194-1202.

    [20] Yoshimura K, Shigeoka S. Versatile physiological functions of the Nudix hydrolase family in Arabidopsis [J]. Bioscience, Biotechnology, and Biochemistry, 2015, 79(3): 354-366.

    [21] 練珊珊, 張寶會, 姚新轉(zhuǎn), 等. 茶樹CsIQD基因家族的鑒定和表達模式分析[J]. 特產(chǎn)研究, 2021, 43(5): 26-34. Lian S S, Zhang B H, Yao X Z, et al. Identification and expression pattern analysis of CsIDQ gene family inL. [J]. Special Wild Economic Animal and Plant Research, 2021, 43(5): 26-34.

    [22] Xu G X, Guo C C, Shan H Y, et al. Divergence of duplicate genes in exon-intron structure [J]. PNAS, 2012, 109(4): 1187-1192.

    [23] Ding Y Q, Wang Y, Qiu C, et al. Alternative splicing in tea plants was extensively triggered by drought, heat and their combined stresses [J]. PeerJ, 2020, 8(15): e8258. doi: 10.7717/peerj.8258.

    [24] Sun J, Chen S L, Dai S X, et al. Ion flux profiles and plant ion homeostasis control under salt stress [J]. Plant Signaling & Behavior, 2009, 4(4): 261-264.

    [25] 彭玉華. 甘菊耐鹽性分析及基因啟動子轉(zhuǎn)化菊花試驗[D]. 北京: 北京林業(yè)大學(xué), 2011. Peng Y H. Salt-tolerance analysis ofand transformation ofpromoters to chrysanthemum [D]. Beijing: Beijing Forestry University, 2011.

    [26] 汪新文. 茉莉酸參與植物逆境脅迫的研究進展[J]. 安徽農(nóng)學(xué)通報, 2008, 14(6): 29-35, 24. Wang X W. Research progress of jasminic acid participating in plant adversity stress [J]. Anhui Agricultural Science Bulletin, 2008, 14(6): 29-35, 24.

    [27] 邢勇翔. 外源茉莉酸甲酯處理影響玫瑰抗旱性與花香合成的機理研究[D]. 揚州: 揚州大學(xué), 2021. Xing Y X. Study on the affecting mechanism of drought resistance and floral fragrance synthesis by treatment of exogenous methyl jasmonate in[D]. Yangzhou: Yangzhou University, 2021.

    [28] Magnard J L, Roccia A, Caissard J C, et al. Biosynthesis of monoterpene scent compounds in roses [J]. Science, 2015, 349(6243): 81-83.

    [29] 臧姝. 玫瑰花香成分生物合成相關(guān)基因和的功能驗證研究[D]. 揚州: 揚州大學(xué), 2017. Zang S. Functional verification ofandgene related to biosynthesis of floral constituents in[D]. Yangzhou: Yangzhou University, 2017.

    Identification and Expression Pattern Analysis ofGene Family in

    CHEN Zhenyan1, ZHANG Xiangqin1, CHEN Lan1, XIE Siyi1, LIU Shuoqian1,2*, TIAN Na1,2*

    1. Key Laboratory of Tea Science of Ministry of Education, Changsha 410128, China; 2. Department of tea science, College of horticulture, Hunan Agricultural university, Changsha 410128, China

    NUDIX hydrolase belongs to pyrophosphatase, which plays an important role in information transmission, plant growth coordination and responses to adversity stresses. In this study, 43family genes were identified based on the tea () genome, and analyzed by bioinformatics and fluorescence quantitative PCR. The results show that the protein molecular weight of 43 CsNUDXs rang from 11.8~89.2?kDa, with 102 to 342 amino acids. The theoretical isoelectric point was from 4.49 to 9.26, and 18.6% of them were stable proteins. According to the evolutionary relationship,is divided into six subfamilies.-acting element analysis of promoter shows thathave many functional elements related to hormone response, adversity stress, growth and development. The expression patterns ofin different organs were analyzed. It was found that the expression levels ofandwere high in fruits, while the expression levels ofandwere extremely low in fruits, and the expression level ofwas extremely low in old leaves. In addition, the results of real-time fluorescence quantitative PCR show that the expression levels of, such as,andincreased first, then decreased and then increased under the treatment of 1?mmol·L-1MeJA. However, under the treatment of 1?mmol·L-1SA, the expression levels of,anddecreased first, then increased and then decreased. While under the treatment of 300?mmol·L-1NaCl, the expression levels of,andincreased first and then decreased. In summary, the basic characteristics and functions ofwere preliminarily analyzed by bioinformatics technology, and it was found thatcould respond to high salt stress, MeJA and SA treatments.

    ,gene family, abiotic stress, expression analysis

    S571.1;Q786

    A

    1000-369X(2023)02-159-14

    2023-01-17

    2023-03-04

    國家自然科學(xué)基金(32172629)、湖南省種業(yè)創(chuàng)新項目(2021NK1008)

    陳振艷,女,碩士研究生,主要從事茶樹栽培育種及分子生物學(xué)方面研究。*通信作者:shuoqianliu@hunau.net;tianna5678@163.com

    猜你喜歡
    基序內(nèi)含子逆境
    EPIYA 基序與幽門螺桿菌感染相關(guān)胃病關(guān)系的研究進展
    超越逆境
    做人與處世(2022年6期)2022-05-26 10:26:35
    How adversity makes you stronger逆境如何讓你更強大
    帶TRS基序突變的新型冠狀病毒威脅更大
    芥藍Aux/IAA家族基因生物信息學(xué)與表達分析
    線粒體核糖體蛋白基因中內(nèi)含子序列間匹配特性分析
    不同方向內(nèi)含子對重組CHO細胞中神經(jīng)生長因子表達的影響
    更 正
    內(nèi)含子的特異性識別與選擇性剪切*
    完形填空Ⅳ
    久久精品夜夜夜夜夜久久蜜豆| 日韩人妻高清精品专区| 亚洲va日本ⅴa欧美va伊人久久| 精品久久久久久久久久免费视频| 国产精品久久久久久久久免| 欧美xxxx性猛交bbbb| 午夜精品在线福利| 日本a在线网址| 欧美日韩中文字幕国产精品一区二区三区| 国产私拍福利视频在线观看| 久久欧美精品欧美久久欧美| 亚洲专区国产一区二区| 国产视频内射| 亚洲av免费高清在线观看| 久久午夜福利片| 69人妻影院| 校园人妻丝袜中文字幕| 我要搜黄色片| 中文字幕熟女人妻在线| АⅤ资源中文在线天堂| 国产精品一区二区三区四区久久| 亚洲 国产 在线| 我的老师免费观看完整版| 成年女人看的毛片在线观看| 亚洲色图av天堂| 国产探花极品一区二区| 小蜜桃在线观看免费完整版高清| 欧美极品一区二区三区四区| 中文字幕精品亚洲无线码一区| 好男人在线观看高清免费视频| 日本 欧美在线| 狠狠狠狠99中文字幕| 国产一区二区激情短视频| 日韩欧美免费精品| 69人妻影院| 免费人成视频x8x8入口观看| 国产一区二区在线av高清观看| 美女黄网站色视频| 久久婷婷人人爽人人干人人爱| 精品久久久久久,| 欧美日本亚洲视频在线播放| 99国产极品粉嫩在线观看| 国产精品综合久久久久久久免费| 精品一区二区三区视频在线观看免费| 国产精品亚洲一级av第二区| 看黄色毛片网站| 欧美日韩亚洲国产一区二区在线观看| 成年人黄色毛片网站| 69av精品久久久久久| 亚洲无线在线观看| 少妇被粗大猛烈的视频| 婷婷色综合大香蕉| 3wmmmm亚洲av在线观看| eeuss影院久久| 3wmmmm亚洲av在线观看| 一级a爱片免费观看的视频| 日韩国内少妇激情av| av在线亚洲专区| 欧美+日韩+精品| 国产精品亚洲美女久久久| 欧美+日韩+精品| 欧美不卡视频在线免费观看| av福利片在线观看| 欧美最新免费一区二区三区| 欧美日韩综合久久久久久 | 亚洲avbb在线观看| 久久久国产成人精品二区| 日本撒尿小便嘘嘘汇集6| 欧美xxxx性猛交bbbb| 成人二区视频| 国产精品久久久久久久久免| 高清在线国产一区| 国产aⅴ精品一区二区三区波| 久久久久精品国产欧美久久久| 国产精品永久免费网站| 在线播放国产精品三级| av黄色大香蕉| 久久久精品大字幕| 黄色欧美视频在线观看| 人妻丰满熟妇av一区二区三区| 男女那种视频在线观看| 亚洲av.av天堂| 国产精品,欧美在线| 嫩草影院精品99| 露出奶头的视频| 少妇的逼好多水| 欧美性猛交黑人性爽| 深爱激情五月婷婷| 午夜亚洲福利在线播放| a级毛片免费高清观看在线播放| 久久久久久国产a免费观看| 一级毛片久久久久久久久女| 色5月婷婷丁香| 日韩大尺度精品在线看网址| 中文字幕高清在线视频| 欧美绝顶高潮抽搐喷水| 亚洲国产精品合色在线| 看黄色毛片网站| 国产高清有码在线观看视频| 噜噜噜噜噜久久久久久91| 国产亚洲精品久久久com| 日韩欧美精品免费久久| 给我免费播放毛片高清在线观看| 夜夜爽天天搞| 99riav亚洲国产免费| 欧美色视频一区免费| 乱系列少妇在线播放| 精品午夜福利视频在线观看一区| 最新在线观看一区二区三区| 日本黄色片子视频| 久久精品国产亚洲av天美| 淫妇啪啪啪对白视频| 欧美最黄视频在线播放免费| 麻豆一二三区av精品| av福利片在线观看| 12—13女人毛片做爰片一| ponron亚洲| 午夜福利在线在线| 丰满的人妻完整版| 亚洲av.av天堂| 哪里可以看免费的av片| 免费大片18禁| 国产亚洲av嫩草精品影院| 一进一出抽搐动态| 国产真实乱freesex| 韩国av在线不卡| 内地一区二区视频在线| 国产aⅴ精品一区二区三区波| 99精品在免费线老司机午夜| 国产精品野战在线观看| 熟妇人妻久久中文字幕3abv| 欧美最新免费一区二区三区| 国产精品综合久久久久久久免费| 久久久午夜欧美精品| 国内揄拍国产精品人妻在线| 国产精品人妻久久久影院| 国内毛片毛片毛片毛片毛片| 男女视频在线观看网站免费| 嫁个100分男人电影在线观看| 日本一二三区视频观看| 少妇人妻精品综合一区二区 | 免费看美女性在线毛片视频| 免费人成视频x8x8入口观看| 精品一区二区三区视频在线| 成年女人毛片免费观看观看9| 国产久久久一区二区三区| 我的女老师完整版在线观看| 亚洲国产精品sss在线观看| 日日撸夜夜添| 俄罗斯特黄特色一大片| 国产高清有码在线观看视频| 亚洲四区av| 校园春色视频在线观看| 色尼玛亚洲综合影院| 欧美一级a爱片免费观看看| 国产欧美日韩精品一区二区| 黄色配什么色好看| 国产高清视频在线播放一区| 人人妻,人人澡人人爽秒播| 国内精品久久久久久久电影| 亚洲成a人片在线一区二区| 99热这里只有精品一区| 国产精品人妻久久久久久| 黄色日韩在线| 国产爱豆传媒在线观看| 欧美bdsm另类| 欧美性猛交黑人性爽| 久久亚洲真实| 真人一进一出gif抽搐免费| 精品99又大又爽又粗少妇毛片 | 国产精品久久久久久精品电影| 国产亚洲91精品色在线| av视频在线观看入口| 精品一区二区三区视频在线| 别揉我奶头 嗯啊视频| 精品人妻1区二区| 99国产极品粉嫩在线观看| 国产私拍福利视频在线观看| 国产一区二区三区在线臀色熟女| 亚洲精品影视一区二区三区av| 日韩亚洲欧美综合| 一卡2卡三卡四卡精品乱码亚洲| 在线播放国产精品三级| 久久久久久久久久成人| 不卡一级毛片| 亚洲美女黄片视频| 久久国产乱子免费精品| 亚洲熟妇熟女久久| 精品午夜福利在线看| 亚洲狠狠婷婷综合久久图片| 俺也久久电影网| 女的被弄到高潮叫床怎么办 | 亚洲精品一区av在线观看| 日韩精品有码人妻一区| 综合色av麻豆| 国产亚洲精品久久久com| 亚洲avbb在线观看| 国产色爽女视频免费观看| 国产精品亚洲美女久久久| 中文亚洲av片在线观看爽| 美女免费视频网站| 亚洲精品色激情综合| av视频在线观看入口| 人妻少妇偷人精品九色| 国产在视频线在精品| 又黄又爽又免费观看的视频| 婷婷精品国产亚洲av| 午夜福利在线在线| 亚洲18禁久久av| videossex国产| 国产单亲对白刺激| 欧美极品一区二区三区四区| 午夜视频国产福利| 成人特级黄色片久久久久久久| 不卡一级毛片| 国产男靠女视频免费网站| 88av欧美| 男女视频在线观看网站免费| 久久精品国产亚洲网站| 老女人水多毛片| 久久精品国产亚洲av涩爱 | 国产精品女同一区二区软件 | 一个人看的www免费观看视频| 中文字幕av成人在线电影| 一a级毛片在线观看| 成人国产麻豆网| 狂野欧美激情性xxxx在线观看| 久久久久久久精品吃奶| 亚洲精品影视一区二区三区av| 亚洲欧美清纯卡通| 成人亚洲精品av一区二区| 看片在线看免费视频| 伦理电影大哥的女人| 亚洲成人久久性| 日韩中文字幕欧美一区二区| 在线观看av片永久免费下载| 国产一级毛片七仙女欲春2| 免费av毛片视频| 久久久久久久精品吃奶| 亚洲欧美日韩高清专用| 日韩欧美在线二视频| 老熟妇乱子伦视频在线观看| 国产精品综合久久久久久久免费| 久久久久久伊人网av| 亚洲aⅴ乱码一区二区在线播放| 桃色一区二区三区在线观看| 在线播放无遮挡| 男人舔奶头视频| 欧美一区二区国产精品久久精品| 桃红色精品国产亚洲av| 色av中文字幕| 久久精品国产亚洲av涩爱 | 欧美成人一区二区免费高清观看| 又爽又黄无遮挡网站| 亚洲性久久影院| 嫩草影视91久久| 国产av不卡久久| 午夜精品久久久久久毛片777| 久久精品国产鲁丝片午夜精品 | 成年女人永久免费观看视频| 麻豆国产av国片精品| 在线播放国产精品三级| 老司机深夜福利视频在线观看| 国产精品综合久久久久久久免费| 色综合亚洲欧美另类图片| 日本a在线网址| 麻豆精品久久久久久蜜桃| 日日摸夜夜添夜夜添av毛片 | 高清毛片免费观看视频网站| 直男gayav资源| 欧美最新免费一区二区三区| 99久久精品热视频| 成人无遮挡网站| 欧洲精品卡2卡3卡4卡5卡区| 老司机深夜福利视频在线观看| 成人国产一区最新在线观看| 亚洲精品粉嫩美女一区| 国产亚洲精品久久久com| 国产色婷婷99| 国产女主播在线喷水免费视频网站 | 日日撸夜夜添| 18禁裸乳无遮挡免费网站照片| 色尼玛亚洲综合影院| 久久人妻av系列| 91在线观看av| 999久久久精品免费观看国产| 国产成年人精品一区二区| 给我免费播放毛片高清在线观看| 欧美极品一区二区三区四区| 婷婷亚洲欧美| 五月伊人婷婷丁香| 69av精品久久久久久| 特大巨黑吊av在线直播| 赤兔流量卡办理| 网址你懂的国产日韩在线| 99久久成人亚洲精品观看| 亚洲av一区综合| 97人妻精品一区二区三区麻豆| 成年女人毛片免费观看观看9| 成人三级黄色视频| 88av欧美| 日本欧美国产在线视频| 天堂网av新在线| 1000部很黄的大片| 亚洲自偷自拍三级| 国产精品一区www在线观看 | 草草在线视频免费看| 毛片一级片免费看久久久久 | 国产免费男女视频| 村上凉子中文字幕在线| 超碰av人人做人人爽久久| 成人亚洲精品av一区二区| 老熟妇乱子伦视频在线观看| 精品人妻视频免费看| 亚洲av第一区精品v没综合| 色视频www国产| 免费看光身美女| 一进一出好大好爽视频| 久久久久国内视频| 国产在视频线在精品| 国产一区二区三区视频了| 丰满人妻一区二区三区视频av| www.www免费av| 淫妇啪啪啪对白视频| 12—13女人毛片做爰片一| 97超视频在线观看视频| 可以在线观看的亚洲视频| 波多野结衣巨乳人妻| 国产精品无大码| 夜夜爽天天搞| 亚洲人成网站在线播| 一区二区三区激情视频| 国产亚洲91精品色在线| 成人av在线播放网站| 赤兔流量卡办理| 成人三级黄色视频| 女人十人毛片免费观看3o分钟| 国产精品久久久久久久久免| 黄色日韩在线| 精品乱码久久久久久99久播| 久久久色成人| 午夜日韩欧美国产| 午夜影院日韩av| 亚洲av熟女| 乱码一卡2卡4卡精品| 精品一区二区三区人妻视频| 欧美日韩中文字幕国产精品一区二区三区| 国产麻豆成人av免费视频| 亚洲国产精品合色在线| 国产主播在线观看一区二区| 天天躁日日操中文字幕| av在线老鸭窝| 大型黄色视频在线免费观看| 亚洲av日韩精品久久久久久密| 一级黄片播放器| 在线国产一区二区在线| 97人妻精品一区二区三区麻豆| 性插视频无遮挡在线免费观看| 人妻久久中文字幕网| 欧美最新免费一区二区三区| 波多野结衣高清作品| 狂野欧美白嫩少妇大欣赏| 欧美绝顶高潮抽搐喷水| 亚洲av.av天堂| 99精品久久久久人妻精品| 欧美+日韩+精品| 我的女老师完整版在线观看| 很黄的视频免费| 欧美日本视频| 久久精品久久久久久噜噜老黄 | 禁无遮挡网站| 午夜视频国产福利| 亚洲成人中文字幕在线播放| 禁无遮挡网站| 18+在线观看网站| 日本 av在线| 亚洲四区av| 中文资源天堂在线| 日日撸夜夜添| 99热6这里只有精品| 国语自产精品视频在线第100页| 日韩中字成人| 久久久色成人| 成人永久免费在线观看视频| 午夜a级毛片| 日韩亚洲欧美综合| 欧美日韩国产亚洲二区| 丰满的人妻完整版| 午夜老司机福利剧场| 午夜日韩欧美国产| 男女做爰动态图高潮gif福利片| 高清毛片免费观看视频网站| 日日夜夜操网爽| av在线观看视频网站免费| 欧美xxxx黑人xx丫x性爽| 中文字幕久久专区| 午夜福利在线在线| 亚洲精华国产精华精| 精品国内亚洲2022精品成人| 国产 一区精品| 国产精华一区二区三区| 国产精品野战在线观看| 精品久久久久久久久亚洲 | 中文资源天堂在线| 99久久精品热视频| 热99re8久久精品国产| 老司机深夜福利视频在线观看| 久久婷婷人人爽人人干人人爱| 毛片女人毛片| 久久久久久久久久成人| 免费不卡的大黄色大毛片视频在线观看 | 国产精品野战在线观看| 欧美最新免费一区二区三区| 欧美性猛交黑人性爽| 99久久成人亚洲精品观看| 99国产极品粉嫩在线观看| 一夜夜www| 精品久久久久久久久久久久久| 亚洲成人久久性| 校园人妻丝袜中文字幕| 国产成人a区在线观看| 日本色播在线视频| 国产真实乱freesex| 婷婷精品国产亚洲av| 中文字幕熟女人妻在线| 国产 一区 欧美 日韩| 五月伊人婷婷丁香| 一本一本综合久久| 人妻制服诱惑在线中文字幕| 变态另类丝袜制服| 亚洲精品乱码久久久v下载方式| 熟女电影av网| 国产一区二区三区在线臀色熟女| 成熟少妇高潮喷水视频| 精品一区二区三区视频在线观看免费| 舔av片在线| 国产色爽女视频免费观看| 婷婷色综合大香蕉| 午夜激情福利司机影院| 在线观看免费视频日本深夜| 国产精品爽爽va在线观看网站| 亚洲av电影不卡..在线观看| 久久欧美精品欧美久久欧美| 美女cb高潮喷水在线观看| 亚洲va日本ⅴa欧美va伊人久久| 美女高潮喷水抽搐中文字幕| 成年女人毛片免费观看观看9| 国产老妇女一区| 我要搜黄色片| 免费在线观看日本一区| 五月伊人婷婷丁香| 午夜激情欧美在线| 亚洲中文字幕日韩| 亚洲av一区综合| 成人欧美大片| 久久久久久伊人网av| 成人二区视频| 嫩草影视91久久| 午夜精品在线福利| 亚洲最大成人中文| 91久久精品国产一区二区三区| 欧美绝顶高潮抽搐喷水| 日韩精品有码人妻一区| bbb黄色大片| 一个人看的www免费观看视频| 久久人人爽人人爽人人片va| 亚洲人成网站在线播| 亚洲美女黄片视频| 亚洲最大成人手机在线| 亚洲av五月六月丁香网| 亚洲av一区综合| 亚洲中文日韩欧美视频| 国产高清激情床上av| 国产日本99.免费观看| 成人av一区二区三区在线看| 国产精品一区二区免费欧美| 精品不卡国产一区二区三区| 在线观看美女被高潮喷水网站| av在线蜜桃| 国产一区二区三区av在线 | 国产精品一区二区三区四区免费观看 | 色综合色国产| а√天堂www在线а√下载| 可以在线观看毛片的网站| 97超视频在线观看视频| 校园春色视频在线观看| 老熟妇仑乱视频hdxx| 毛片一级片免费看久久久久 | 成人国产综合亚洲| 舔av片在线| 看片在线看免费视频| 精品一区二区免费观看| 成人亚洲精品av一区二区| 久久人人精品亚洲av| bbb黄色大片| 精品不卡国产一区二区三区| 亚洲国产精品成人综合色| 欧美精品啪啪一区二区三区| 精品乱码久久久久久99久播| 成年免费大片在线观看| 男人舔奶头视频| 一进一出抽搐动态| 精品免费久久久久久久清纯| 免费观看精品视频网站| 久久精品国产亚洲av香蕉五月| 国内精品久久久久久久电影| 1000部很黄的大片| 最近最新免费中文字幕在线| 联通29元200g的流量卡| 搡女人真爽免费视频火全软件 | 五月玫瑰六月丁香| 中文在线观看免费www的网站| 日本与韩国留学比较| 日本撒尿小便嘘嘘汇集6| 亚洲专区国产一区二区| 欧美一区二区亚洲| 色综合婷婷激情| 在线免费观看的www视频| 亚洲美女搞黄在线观看 | 成人国产一区最新在线观看| 嫁个100分男人电影在线观看| 成人国产综合亚洲| 欧美激情在线99| 少妇高潮的动态图| 久久久久久久精品吃奶| 中文字幕高清在线视频| 琪琪午夜伦伦电影理论片6080| 欧美最黄视频在线播放免费| 国产成年人精品一区二区| 久久国产精品人妻蜜桃| 久久精品夜夜夜夜夜久久蜜豆| 热99re8久久精品国产| 亚洲国产高清在线一区二区三| 两个人视频免费观看高清| 日韩在线高清观看一区二区三区 | 男人和女人高潮做爰伦理| 成年女人看的毛片在线观看| 一本久久中文字幕| 夜夜看夜夜爽夜夜摸| 久久国内精品自在自线图片| 亚洲电影在线观看av| 亚洲最大成人手机在线| 亚洲中文字幕一区二区三区有码在线看| 亚洲av熟女| 搞女人的毛片| а√天堂www在线а√下载| 亚洲国产精品合色在线| 天天一区二区日本电影三级| 一夜夜www| 日韩国内少妇激情av| 欧美zozozo另类| 变态另类成人亚洲欧美熟女| 成人无遮挡网站| 小说图片视频综合网站| 性色avwww在线观看| 国产极品精品免费视频能看的| 国产成人影院久久av| 国产精品1区2区在线观看.| 日韩高清综合在线| 国产一级毛片七仙女欲春2| 午夜精品久久久久久毛片777| 国产熟女欧美一区二区| 久久久色成人| 日本免费一区二区三区高清不卡| 日韩高清综合在线| 不卡一级毛片| 春色校园在线视频观看| av在线蜜桃| 中文字幕精品亚洲无线码一区| 大型黄色视频在线免费观看| 国产美女午夜福利| 国产三级在线视频| 免费看日本二区| 亚洲美女视频黄频| av在线老鸭窝| 日本免费a在线| 91久久精品电影网| 亚洲美女搞黄在线观看 | 别揉我奶头~嗯~啊~动态视频| 日韩欧美精品v在线| 免费人成在线观看视频色| 日本 欧美在线| 国产精品久久久久久亚洲av鲁大| 麻豆av噜噜一区二区三区| 成人国产麻豆网| 99久久精品一区二区三区| 成人无遮挡网站| av在线观看视频网站免费| 内地一区二区视频在线| 99热只有精品国产| 看黄色毛片网站| 在线观看午夜福利视频| 免费看a级黄色片| 观看免费一级毛片| 亚洲av熟女| 22中文网久久字幕| 国产大屁股一区二区在线视频| 亚洲狠狠婷婷综合久久图片| 国产午夜福利久久久久久| 精品99又大又爽又粗少妇毛片 | 久久久久久久久久成人| 简卡轻食公司| 免费av毛片视频| 亚洲va日本ⅴa欧美va伊人久久| 男人和女人高潮做爰伦理| 热99在线观看视频| 国产伦在线观看视频一区| 久久久国产成人免费| 日本 av在线| 伦精品一区二区三区| 国产真实乱freesex| 国产精品一及| 成人国产综合亚洲| 两人在一起打扑克的视频| 免费不卡的大黄色大毛片视频在线观看 |