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

    普通小麥 TaCAT 9基因的克隆與表達(dá)分析

    2018-07-02 08:40:28褚瑩瑩王晨陽陳雨露劉衛(wèi)星胡陽陽李莎莎
    麥類作物學(xué)報(bào) 2018年6期
    關(guān)鍵詞:同源克隆氨基酸

    褚瑩瑩,王晨陽,2,陳雨露,康 娟,馬 耕,劉衛(wèi)星,胡陽陽,李莎莎

    (1.河南農(nóng)業(yè)大學(xué)農(nóng)學(xué)院,河南鄭州 450002; 2.國家小麥工程技術(shù)研究中心,河南鄭州 450002)

    氨基酸是植物體內(nèi)氮素存在的主要形式,能夠作為含氮通貨參與到各種代謝途徑中,而這個(gè)過程有賴于氨基酸轉(zhuǎn)運(yùn)蛋白的運(yùn)輸功能,因此研究氨基酸轉(zhuǎn)運(yùn)蛋白的功能對(duì)于植物氮代謝的研究具有重要意義[1-2]。植物中的氨基酸轉(zhuǎn)運(yùn)蛋白可分為五大超家族,但研究較多的主要有兩個(gè)超家族:一個(gè)是氨基酸轉(zhuǎn)運(yùn)蛋白超家族(ATFs),約有46個(gè)成員,至少又可分為5個(gè)亞家族;另一個(gè)是氨基酸-多胺-膽堿轉(zhuǎn)運(yùn)蛋白超家族(APCs),有14個(gè)成員,又分為L型氨基酸轉(zhuǎn)運(yùn)蛋白(LATs)和陽離子氨基酸轉(zhuǎn)運(yùn)蛋白(CATs)兩個(gè)亞家族[3-5]。研究表明,除了脯氨酸轉(zhuǎn)運(yùn)蛋白,所有的氨基酸轉(zhuǎn)運(yùn)蛋白都表現(xiàn)出顯著的底物專一性[6]。許多研究證明氨基酸轉(zhuǎn)運(yùn)蛋白的表達(dá)具有組織及器官特異性,且受多種因素的影響,如光、線蟲感染、蔗糖、滲透勢、銨鹽、硝酸鹽及谷氨酸等[7-9]。目前對(duì)ATFs超家族的研究較為透徹,而對(duì)APCs超家族的研究相對(duì)較少,本試驗(yàn)的研究對(duì)象CAT9是CATs亞家族中的一員[9]。

    對(duì)植物CATs亞家族的研究報(bào)道主要集中在模式植物中,擬南芥CATs亞家族有9個(gè)成員,大都定位于質(zhì)膜或液泡膜上,具有11~14個(gè)推測的跨膜結(jié)構(gòu)域[9-10]。AtCAT1啟動(dòng)子與GUS融合表達(dá)分析顯示,該基因主要高表達(dá)于衰老的葉片,在初生根的根尖及花中也能檢測到少量表達(dá),猜測該基因可能參與氨基酸向種子的運(yùn)輸[11]。AtCAT2主要定位于液泡膜,可能通過參與氨基酸的跨液泡膜轉(zhuǎn)運(yùn)來調(diào)節(jié)葉片可溶性氨基酸濃度;AtCAT3和AtCAT4分別定位于高爾基體膜和液泡膜,具體功能有待研究[12]。AtCAT5主要在籽粒中表達(dá),可能在籽粒形成或種子萌發(fā)過程中裝載種柄中的氨基酸,又有研究表明該基因在真葉邊緣部位的表達(dá)量也比較高,可能參與重新吸收葉緣部位滲漏的氨基酸[13-14]。AtCAT6只在根原基、花及種子等作為庫的組織器官中表達(dá),而且在根部的表達(dá)水平能夠受到根結(jié)線蟲感染的誘導(dǎo)[9]。AtCAT7在正常環(huán)境下的表達(dá)水平非常的低[13],AtCAT8在根和植株上部的分生組織中有較高的表達(dá)豐度,其功能可能是介導(dǎo)氨基酸特別是谷氨酸和谷氨酰胺向植株幼嫩分生組織分配[14]。研究表明,擬南芥中AtCAT9基因主要定位于囊泡膜上,具有比較廣泛的表達(dá)特性,在大多數(shù)組織器官中都有表達(dá),可能參與調(diào)控液泡內(nèi)外氨基酸濃度的平衡,過表達(dá)AtCAT9的擬南芥植株在低氮脅迫下具有生長發(fā)育延遲、葉片衰老緩慢的現(xiàn)象[15]。番茄中的研究表明,CAT9可能充當(dāng)膜兩側(cè)GLU/ASP與GABA之間交換的媒介,導(dǎo)致番茄成熟過程中大液泡內(nèi)GLU/ASP增加,GABA減少,從而影響番茄的口感與營養(yǎng)品質(zhì)[16]。

    小麥?zhǔn)侨蜃钪匾募Z食作物之一,其產(chǎn)量高低直接關(guān)系著世界糧食安全。研究氮代謝相關(guān)基因的功能對(duì)提高氮素利用率及產(chǎn)量具有重要意義。本實(shí)驗(yàn)室通過前期的蛋白組學(xué)研究發(fā)現(xiàn)小麥幼苗葉片中的TaCAT9蛋白參與低氮脅迫的應(yīng)答,說明其可能在氮素脅迫條件下起著某種調(diào)控作用(數(shù)據(jù)未發(fā)表)。通過查閱文獻(xiàn)發(fā)現(xiàn),目前有關(guān)小麥TaCAT9基因及CATs基因家族的研究尚未見報(bào)道。鑒于此,本研究采用同源克隆的方法,從普通小麥品種豫麥49-198中克隆TaCAT9兩個(gè)部分同源基因的cDNA序列,對(duì)其進(jìn)行生物信息學(xué)分析,同時(shí)研究了TaCAT9在不同組織器官中及氮素脅迫下的表達(dá)特性,以期為進(jìn)一步探索TaCAT9的生物學(xué)功能奠定基礎(chǔ)。

    1 材料與方法

    1.1 材料及處理

    供試材料為河南省大面積種植的小麥品種豫麥49-198,該品種具有豐產(chǎn)、優(yōu)質(zhì)、抗逆性好的特點(diǎn),由河南平安種業(yè)有限公司惠贈(zèng)。

    選取籽粒飽滿的豫麥49-198種子,經(jīng)75%酒精浸泡10 min后,用蒸餾水清洗干凈,置于干凈的濕毛巾上,在黑暗條件下放置48 h,期間保持毛巾濕潤狀態(tài)。待種子露白后,選取發(fā)芽一致的種子置于滅菌培養(yǎng)皿中,在光暗周期16/8 h、光照強(qiáng)度250 μmol·m-2·s-1、晝夜溫度25/15 ℃、晝夜相對(duì)濕度60%和70%的培養(yǎng)箱中培養(yǎng)。培養(yǎng)液為Hoagland全營養(yǎng)液,每天更換一次。待幼苗長至兩葉一心,將幼苗做兩種處理:一種為對(duì)照組,在Hoagland全營養(yǎng)液中生長;另一種為處理組,生長于缺氮的Hoagland營養(yǎng)液(用KCl和CaCl2代替Hoagland營養(yǎng)液中的KNO3和Ca(NO3)2)中。分別在0、0.25、0.5、1和7 d時(shí)取對(duì)照組與處理組的小麥幼苗葉片用于低氮脅迫下目的基因的表達(dá)模式分析。剪取正常處理兩葉一心期的葉片用于目的基因的克隆。所有材料用液氮速凍后保存于-80 ℃冰箱備用。

    目的基因組織特異性表達(dá)分析:于2016年10月將豫麥49-198種植于河南農(nóng)業(yè)大學(xué)科教園區(qū)試驗(yàn)田,正常水肥管理,即施純氮240 kg·hm-2、磷肥(P2O5)150 kg·hm-2、鉀肥(K2O)120 kg·hm-2,分別于拔節(jié)期和揚(yáng)花期灌水,水表計(jì)量,每次灌水量750 m3·hm-2。小區(qū)面積為21 m2(7 m×3 m)。由于從2010年起,連續(xù)六年對(duì)該小區(qū)采取了精細(xì)化的栽培管理措施,因此土壤養(yǎng)分等條件基本一致。在小麥開花期選取花期與長勢一致的植株掛牌標(biāo)記,于花后15 d取標(biāo)記植株的根、莖、葉和籽粒,液氮速凍后保存于-80 ℃ 冰箱備用。

    1.2 RNA的提取及cDNA第一鏈的合成

    小麥根、莖、葉和籽??俁NA的提取采用Trizol法,具體提取步驟參照Trizol Reagent(北京索萊寶科技有限公司)說明書進(jìn)行。用1%的瓊脂糖凝膠電泳檢測提取總RNA的完整性,并采用ND-1000微量紫外分光光度計(jì)(基因有限公司)測定RNA的濃度和純度。總RNA經(jīng)DNaseⅠ(普洛麥格(北京)生物技術(shù)有限公司)處理后作為模板,參照PrimeScriptTMⅡ 1st Strand cDNA Synthesis Kit(寶生物工程(大連)有限公司)合成cDNA第一條鏈,-20 ℃保存?zhèn)溆谩?/p>

    1.3 引物設(shè)計(jì)及TaCAT9基因的克隆

    以二穗短柄草CAT9(XM_003569968.3)的cDNA序列為探針在Ensembl Plants六倍體小麥數(shù)據(jù)庫(http://plants.ensembl.org/Triticum_aestivum/Info/Index)進(jìn)行比對(duì),得到三條高度相似且分別定位在6 A、6 B和6 D長臂上的序列,對(duì)應(yīng)的序列號(hào)分別為TRIAE_CS42_6AL_TGACv1_471594_AA1511590.1,RIAE_CS42_6BL_TGACv1_501182_AA1614650.1和TRIAE_CS42_6DL_TGACv1_526398_AA1682300.1;長度分別為2 462、2 572和2 365 bp,均包含一個(gè)長度為1 818 bp的開放閱讀框。比對(duì)這三條序列開放閱讀框的上下游,利用Primer5.0軟件在三條序列相似的部位設(shè)計(jì)兩對(duì)引物TaCAT9-F-1/TaCAT9-R和TaCAT9-F-2/TaCAT9-R用于三個(gè)部分同源基因cDNA序列的擴(kuò)增。由于上游引物5′端第五個(gè)堿基在三條序列的結(jié)合位點(diǎn)上存在一個(gè)由T到C的替換,所以設(shè)計(jì)兩條上游引物TaCAT9-F-1和TaCAT9-F-2,它們只在的5′端第五個(gè)堿基存在差異,分別為T和C。引物由北京奧科鼎盛生物科技有限公司合成,序列信息見表1。以豫麥49-198兩葉一心時(shí)期葉片的cDNA為模板進(jìn)行PCR擴(kuò)增,所用的酶為2×Phanta Max Master Mix(南京諾唯贊生物科技有限公司)。擴(kuò)增體系與程序參照說明書,其中循環(huán)數(shù)為32,退火溫度為56 ℃,延伸時(shí)間為2 min 30 s。擴(kuò)增結(jié)束后,用1%的瓊脂糖對(duì)PCR產(chǎn)物進(jìn)行電泳檢測,回收純化目的片段后,將目的片段連接至克隆載體pTOPO-Blunt simple(北京愛德萊生物科技有限公司),轉(zhuǎn)化大腸桿菌感受態(tài)細(xì)胞DH10B(北京中美泰和生物技術(shù)有限公司),從兩對(duì)引物組合的擴(kuò)增產(chǎn)物中挑選經(jīng)PCR鑒定正確的陽性克隆送至北京奧科鼎盛生物科技有限公司進(jìn)行測序。

    1.4 TaCAT9基因編碼蛋白的生物信息學(xué)分析

    運(yùn)用ExPASy-ProtParam(https://web.expasy.org/protparam/)對(duì)TaCAT9基因編碼蛋白的理化性質(zhì)進(jìn)行預(yù)測;運(yùn)用WOLF PSORT(https:// psort.hgc.jp/)和TMHMM(http://www.cbs.dtu.dk/services/TMHMM/)對(duì)編碼蛋白的信號(hào)肽、亞細(xì)胞定位及跨膜區(qū)進(jìn)行預(yù)測;運(yùn)用InterProScan(http://www.ebi.ac.uk/interpro/)對(duì)編碼蛋白保守結(jié)構(gòu)域進(jìn)行預(yù)測。利用SOPMA(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=/NPSA/npsa_sopma.html)預(yù)測TaCAT9蛋白質(zhì)的二級(jí)結(jié)構(gòu),用Phyre 2(http://www.sbg.bio.ic.ac.uk/Phyre2/html/page.cgi?id=index)預(yù)測蛋白質(zhì)的三級(jí)結(jié)構(gòu)。以TaCAT9蛋白序列為探針,在NCBI(https://www.ncbi.nlm.nih.gov/)搜索不同物種的同源蛋白,并利用DNAMAN5.2.2進(jìn)行蛋白序列的多重比對(duì),利用MEGA5.0構(gòu)建基于neigbor-joining(bootstrap值為1 000)的系統(tǒng)發(fā)育樹。

    1.5 TaCAT9基因的表達(dá)分析

    由于序列高度相似,未能設(shè)計(jì)出分別針對(duì)TaCAT9-A和TaCAT9-B的特異引物,只設(shè)計(jì)了一對(duì)針對(duì)兩序列的特異性引物。利用NCBI在線引物設(shè)計(jì)工具primer-BLAST(https://www.ncbi.nlm.nih.gov/tools/primer-blast/index.cgi?LINK_LOC=BlastHome)在TaCAT9-A和TaCAT9-B序列的保守區(qū)(避開核苷酸突變點(diǎn))設(shè)計(jì)引物對(duì)TaCAT9-qPCR-F/TaCAT9-qPCR-R(表1),用于檢測TaCAT9基因的相對(duì)表達(dá)量;內(nèi)參基因采用小麥看家基因β-Actin(AB181991),所用引物Actin-F/Actin-R序列見表1[17]。qRT-PCR體系與程序參照SYBR Premix ExTaqTMⅡ(寶生物工程(大連)有限公司)說明書,其中循環(huán)數(shù)為40,退火溫度為60 ℃,延伸時(shí)間為30 s。所用反應(yīng)平臺(tái)為CFX Connect Real-time Systerm(伯樂生命醫(yī)學(xué)產(chǎn)品(上海)有限公司),每個(gè)樣品設(shè)置3次生物學(xué)重復(fù)。目的基因相對(duì)表達(dá)量的計(jì)算采用2-△△CT法,用SPSS 19.0對(duì)數(shù)據(jù)進(jìn)行顯著性分析,用Excel 2003作圖。

    表1 本研究中所用引物的基本信息Table 1 Primers used in this study

    2 結(jié)果與分析

    2.1 小麥TaCAT9基因克隆結(jié)果

    利用引物TaCAT9-F1/TaCAT9-R和TaCAT9-F2/TaCAT9-R分別對(duì)豫麥49-198的cDNA進(jìn)行擴(kuò)增,分別得到兩條大小約為2 000 bp的特異性條帶(圖1)。將目的條帶回收純化后,經(jīng)連接、轉(zhuǎn)化,各篩選30個(gè)陽性克隆進(jìn)行測序。測序結(jié)果與中國春的三條序列(見材料方法)進(jìn)行比對(duì),結(jié)果表明:TaCAT9-F1/TaCAT9-R擴(kuò)增產(chǎn)物連接轉(zhuǎn)化的30個(gè)陽性單克隆的序列信息都與TRIAE_CS42_6AL_TGACv1_471594_AA1511590.1相同,而TaCAT9-F2/ TaCAT9-R有3個(gè)單克隆的序列與TRIAE_CS42_6BL_TGACv1_501182_AA1614650.1相同,其余都與TRIAE_CS42_6AL_TGACv1_471594_AA1511590.1相同。所以本試驗(yàn)共得到兩個(gè)TaCAT9基因的cDNA序列,因其分別位于小麥A和B亞基因組上,故將其分別命名為TaCAT9-A和TaCAT9-B。對(duì)于TaCAT9-D,可能由于表達(dá)量過低或者沒有表達(dá),本研究沒能克隆到其cDNA序列。通過比對(duì)分析發(fā)現(xiàn),兩序列間存在28個(gè)單核苷酸突變位點(diǎn),其中8個(gè)為有義突變,導(dǎo)致7個(gè)氨基酸的改變,分別為第116、136、154、652、1 420、1 546、1 547和1 711位堿基,其中第1 546位與第1 547位共同導(dǎo)致一個(gè)氨基酸的改變;其余20個(gè)為沉默突變 ,沒有導(dǎo)致氨基酸的變化(圖2)。

    1:引物TaCAT9-F1/TaCAT9-R 的擴(kuò)增產(chǎn)物;2:引物TaCAT9-F2/TaCAT9-R 的擴(kuò)增產(chǎn)物;M:Trans2k Plus DNA marker。

    1:Amplified products of primer TaCAT9-F1/TaCAT9-R;2:Amplified products of primer TaCAT9-F2/TaCAT9-R;M:Trans2k Plus DNA marker.

    圖1TaCAT9基因的RT-PCR擴(kuò)增產(chǎn)物

    Fig.1AmplifiedproductsofTaCAT9genebyRT-PCR

    圖2 TaCAT9-A和TaCAT9-B的cDNA序列比對(duì)Fig.2 Alignment of cDNA sequences between TaCAT9-A and TaCAT9-B

    2.2 TaCAT9基因編碼蛋白的生物信息學(xué)分析

    2.2.1TaCAT9基因編碼蛋白的理化性質(zhì)

    由于序列相似度極高,protparam在線預(yù)測分析顯示TaCAT9-A和TaCAT9-B兩蛋白基本理化性質(zhì)非常接近,分子式分別為C2959H4619N743O802S19和C2962H4626N746O799S19,分子量分別為64.04和64.08 kDa,理論等電點(diǎn)(pI)分別為8.23和8.27,屬于弱堿性蛋白。脂肪族氨基酸指數(shù)分別為114.91和114.76,負(fù)電荷的氨基酸殘基數(shù)(Asp+Glu)分別為31和30;正電荷的氨基酸殘基數(shù)(Arg+Lys)分別為34和35??偲骄H水性值(GRAVY)分別為0.745和0.739,不穩(wěn)定系數(shù)分別為35.62和35.68,兩者都屬于疏水穩(wěn)定蛋白。

    2.2.2TaCAT9基因編碼蛋白的亞細(xì)胞定位、跨膜區(qū)及保守結(jié)構(gòu)域

    WOLF PSORT預(yù)測分析顯示,兩蛋白均定位于細(xì)胞膜上,沒有分泌通路信號(hào)肽,即兩蛋白都不屬于分泌型蛋白。TMHMM在線預(yù)測表明TaCAT9-A和TaCAT9-B蛋白都存在13個(gè)跨膜區(qū),屬于膜蛋白,這和前人對(duì)CAT亞家族蛋白跨膜區(qū)的描述相吻合。InterProScan保守結(jié)構(gòu)域分析顯示(圖3),兩蛋白都屬于氨基酸透性酶家族,且都含有一個(gè)典型的陽離子氨基酸轉(zhuǎn)運(yùn)蛋白的保守C末端。

    2.2.3TaCAT9基因編碼蛋白的氨基酸序列比對(duì)及進(jìn)化分析

    利用DNAMAN將TaCAT9-A和TaCAT9-B蛋白的氨基酸序列與來自玉米(NP_001168279.1)、水稻(XP_015627412.1)和高粱(XP_002454544.1)的CAT9氨基酸序列進(jìn)行多重比對(duì)。結(jié)果(圖4)表明,4個(gè)物種間的CAT9蛋白氨基酸序列長度差別不大,都約有600個(gè)氨基酸,且同源性較高。TaCAT9-A和TaCAT9-B兩者的同源性達(dá)到98.84%,4個(gè)物種間的同源性達(dá)到91.26% ,氨基酸序列高度保守,都包含13個(gè)跨膜區(qū)和1個(gè)保守的陽離子氨基酸轉(zhuǎn)運(yùn)蛋白的C末端。

    圖3 TaCAT9-A和TaCAT9-B蛋白的保守結(jié)構(gòu)域分析Fig.3 Conserved domains of TaCAT9-A protein and TaCAT9-B protein

    方框指示陽離子氨基酸轉(zhuǎn)運(yùn)蛋白的C末端,橫線指示13個(gè)跨膜區(qū)。

    The box indicates the C-terminal of cationic amino acid transporter and the horizontal lines indicate 13 transmembrane regions.

    圖4TaCAT9與其他物種同源蛋白的多序列比對(duì)

    Fig.4Multi-alignmentofTaCAT9aminoacidsequencewithotherCAT9s

    進(jìn)化分析(圖5)表明,所有參比的CAT9蛋白可以被分為兩大類,小麥等單子葉植物為一類,擬南芥等雙子葉植物為一類。其中與TaCAT9-A蛋白親緣關(guān)系最近的為烏拉爾圖小麥的CAT9蛋白,而山羊草CAT9蛋白與前二者聚類在一起,TaCAT9-B蛋白與上述三者聚類在一起。

    圖5 小麥與其他物種CAT9蛋白的進(jìn)化分析Fig.5 Phylogene tree of CAT9 proteins from wheat and other species表2 TaCAT9蛋白的二級(jí)結(jié)構(gòu)預(yù)測及所占比例Table 2 Prediction proportion of secondary structure of TaCAT9 protein %

    2.2.4TaCAT9基因編碼蛋白的二級(jí)及三級(jí)結(jié)構(gòu)預(yù)測結(jié)果

    通過SOPMA在線預(yù)測分析了TaCAT9-A和TACAT9-B的二級(jí)結(jié)構(gòu)組成(表2),兩者都是由α-螺旋、延伸直鏈、β-轉(zhuǎn)角和無規(guī)則卷曲四種二級(jí)結(jié)構(gòu)組成,兩個(gè)蛋白所含的各種二級(jí)結(jié)構(gòu)比例相近,其中α-螺旋所占比例最高,無規(guī)則卷曲次之。通過Phyre 2構(gòu)建了蛋白的3D模型,因TaCAT9-A 與TACAT9-B氨基酸序列高度相似,所以同源建模選用的蛋白模型相同,3D結(jié)構(gòu)圖也相同。該模型中有73%的氨基酸序列結(jié)構(gòu)預(yù)測置信度達(dá)到100%,因此該蛋白質(zhì)結(jié)構(gòu)模型可信度較高。由圖6可知,兩蛋白主要由13個(gè)α-螺旋組成,這與TMHMM預(yù)測的兩個(gè)蛋白都有13個(gè)跨膜區(qū)一致。

    2.3 TaCAT9基因的表達(dá)特性分析

    在根、莖、葉和籽粒中均可檢測到TaCAT9基因的表達(dá),但在根和籽粒中的表達(dá)量最低,莖中次之,葉片中的表達(dá)量最高,推測TaCAT9基因主要在葉片中表達(dá)(圖7)。進(jìn)一步的分析顯示,在氮脅迫下,小麥葉片TaCAT9基因的表達(dá)呈上升趨勢,在脅迫初期變化不明顯,隨著脅迫時(shí)間的增加表達(dá)量逐漸上調(diào)。

    圖6 TaCAT9蛋白的三級(jí)結(jié)構(gòu)預(yù)測Fig.6 Prediction of three dimensional structure of TaCAT9 protein

    圖柱上的小寫字母不同表示在0.05水平差異顯著。

    The different lower-case letters above the bar indicate significant difference at 0.05 level.

    圖7TaCAT9基因的組織器官表達(dá)模式(A)和氮脅迫下TaCAT9基因在葉片中的表達(dá)模式(B)

    Fig.7ExpressionpatternsofTaCAT9indifferenttissuesofwheatandexpressionpatternsofTaCAT9geneinleavesundernitrogenstress

    3 討 論

    普通小麥?zhǔn)且粋€(gè)典型的異源六倍體物種,由三個(gè)序列高度同源的亞基因組A、B和D組成[18]。因此,小麥中大部分基因在A、B和D三個(gè)亞基因組上都有一個(gè)基因位點(diǎn),這三個(gè)基因稱為部分同源基因。研究表明,在六倍體小麥的產(chǎn)生以及進(jìn)化過程中,部分同源基因的表達(dá)可能產(chǎn)生分化[1,19]。如小麥膨脹素基因(TaEXPA1),它的三個(gè)部分同源基因TaEXPA1-A/B/D在幼苗的根中都是沉默的,但TaEXPA1-A和TaEXPA1-D在幼葉中表達(dá),而TaEXPA1-B在該組織中則是沉默的[20]。小麥TaWLHSI-A/B/D在不同的組織器官中具有不同的表達(dá)模式:TaWLHSI-B幾乎在所有組織器官中的表達(dá)量都低于另外兩個(gè)基因,TaWLHSI-A在幼穗中表達(dá)量最高,TaWLHSI-D則在稃中高表達(dá)[21]。以上研究表明,三個(gè)部分同源基因可能通過表達(dá)的分化來實(shí)現(xiàn)功能的相互協(xié)調(diào),從而來調(diào)控小麥的生長發(fā)育。這一調(diào)控網(wǎng)絡(luò)的復(fù)雜性,大大增加了小麥遺傳和功能分析的難度。在研究小麥的某個(gè)基因時(shí),研究者需要面對(duì)三個(gè)序列高度同源的基因,由于缺乏小麥基因組信息,使得區(qū)分部分同源基因的難度會(huì)很大,所以很多研究者并沒有對(duì)它們加以區(qū)分。2014年,IWGSC(國際小麥基因組測序組織)采用先分離染色體臂,再對(duì)染色體臂分別測序的方法,發(fā)表了基于染色體測序的六倍體小麥中國春的序列草圖,這為研究小麥部分同源基因的遺傳和功能分析提供了一個(gè)很好的參考[22]。本研究所用的Ensemble plants網(wǎng)站中的小麥數(shù)據(jù)庫是基于中國春基因組信息而形成的,其中很多基因的轉(zhuǎn)錄本信息都是通過分析軟件預(yù)測的,而且不同品種間的序列信息也可能有所不同,比如存在SNP位點(diǎn),存在不同的剪切體,這些都需要實(shí)驗(yàn)加以驗(yàn)證[23-24]。本研究通過二穗短柄草CAT9基因的 cDNA序列在Ensembl plants小麥數(shù)據(jù)庫中比對(duì)到三條高度同源的轉(zhuǎn)錄本,分別位于小麥6A、6B和6D染色體上,可以初步判定為TaCAT9的三個(gè)部分同源基因。以中國春的這三條轉(zhuǎn)錄本為參考,設(shè)計(jì)引物,從豫麥49-198克隆到位于6A和6B上的兩條序列,分別將其命名為TaCAT9-A和TaCAT9-B。比對(duì)分析發(fā)現(xiàn)克隆的兩條序列與中國春的序列相似性達(dá)到100%,猜測該基因序列在不同品種間較為保守??赡苡捎贒基因組上的序列表達(dá)量過低或者沒有表達(dá),所以本研究沒有克隆到。另一方面從篩選單克隆測序的結(jié)果來看,本研究所采用的組織材料中,A基因組上的序列表達(dá)量最高,其他兩種遠(yuǎn)遠(yuǎn)低于該序列,推測三個(gè)部分同源基因中TaCAT9-A基因發(fā)揮主要功能。當(dāng)然不排除在其他組織器官或受到內(nèi)外因素的刺激后,三個(gè)基因的表達(dá)量會(huì)發(fā)生變化,這需要進(jìn)一步的研究。

    本研究從豫麥49-198中克隆到的兩個(gè)TaCAT9基因的cDNA序列編碼區(qū)長度都為1 818 bp,編碼605個(gè)氨基酸,它們所編碼的蛋白無論從基本理化性質(zhì)還是結(jié)構(gòu)特征上都高度相似。保守結(jié)構(gòu)域分析表明,兩蛋白都有一個(gè)陽離子氨基酸轉(zhuǎn)運(yùn)蛋白家族的C末端,且都有13個(gè)跨膜結(jié)構(gòu),這都符合模式作物中CATs亞家族的結(jié)構(gòu)特點(diǎn)。由于尚未發(fā)現(xiàn)與小麥TaCAT9基因相關(guān)的研究報(bào)道,因此了解該基因的功能只能通過參考其他物種中已有的報(bào)道。擬南芥中的研究表明氨基酸轉(zhuǎn)運(yùn)蛋白是個(gè)龐大的多基因家族,而CATs作為APCs家族中的一個(gè)亞族,主要參與中性氨基酸的運(yùn)輸[14]。2015年,Christopher等[16]通過定量蛋白組學(xué)的方法分析了番茄由綠變紅過程中果肉細(xì)胞液泡膜蛋白的變化,發(fā)現(xiàn)CAT9蛋白占液泡膜總蛋白的比例由0.02% 升至0.12%,亞細(xì)胞定位分析表明該蛋白不僅存在于液泡膜,而且廣泛存在于細(xì)胞的各種膜結(jié)構(gòu)中;進(jìn)一步研究表明,該蛋白可能在番茄果實(shí)成熟過程中充當(dāng)液泡膜內(nèi)外GLU/ASP與GABA交換的載體。同一年,Yang等[15]在擬南芥中研究發(fā)現(xiàn),在氮素虧缺的生長條件下,AtCAT9缺失會(huì)造成植株萎黃失綠,早衰死亡,而過表達(dá)AtCAT9則完全相反,植株發(fā)育遲緩,存活率升高,葉片的持綠性明顯增強(qiáng);進(jìn)一步分析發(fā)現(xiàn)相對(duì)于AtCAT9缺失植株,轉(zhuǎn)基因植株葉片的可溶性氨基酸含量顯著降低,這說明過表達(dá)AtCAT9有利于提高氮素利用率,緩解氮素缺失造成的不利影響。本研究中,WOLF PSORT預(yù)測分析顯示TaCAT9-A和TaCAT9-B均定位于細(xì)胞膜,說明TaCAT9為膜蛋白,但軟件預(yù)測的結(jié)果較為模糊,該蛋白也有可能位于細(xì)胞的多種膜結(jié)構(gòu)中,需要進(jìn)一步的實(shí)驗(yàn)去明確小麥中CAT9的亞細(xì)胞定位。另外,采用同源建模的方法構(gòu)建TaCAT9蛋白的 3D模型時(shí),Phyre2服務(wù)器預(yù)測該蛋白可能為GLU與GABA的逆向轉(zhuǎn)運(yùn)蛋白,這和番茄中的研究結(jié)果較為相似,所以猜測GLU與GABA這兩種氨基酸可能為TaCAT9轉(zhuǎn)運(yùn)的底物。qRT-PCR結(jié)果顯示,TaCAT9基因主要在小麥葉片中表達(dá),莖中次之,根和籽粒中最少;與正常供氮的植株相比,缺氮植株葉片TaCAT9基因在脅迫初期無變化,隨著脅迫時(shí)間的延長該基因的表達(dá)量顯著上調(diào),說明TaCAT9參與低氮脅迫應(yīng)答。而脅迫響應(yīng)表現(xiàn)出一定的滯后可能是因?yàn)闊o機(jī)氮不是氨基酸轉(zhuǎn)運(yùn)蛋白的直接作用底物,其在根部吸收后,需要同化為氨基酸及經(jīng)過一系列反應(yīng)才能影響到氨基酸轉(zhuǎn)運(yùn)蛋白的表達(dá),所以缺氮不會(huì)立即影響到TaCAT9的表達(dá)。

    參考文獻(xiàn):

    [1] 盧永恩.水稻谷氨酸合酶基因和胞質(zhì)異檸檬酸脫氫酶基因的功能研究以及氨基酸轉(zhuǎn)運(yùn)蛋白基因家族分析[D].武漢:華中農(nóng)業(yè)大學(xué),2014:31.

    LU Y E.Functional analysis of glutamate synthase genes and cytosolic isocitrate dehydrogenase genes and genome-wide analysis of amino acid transporter gene family in rice [D].Wuhan:Huazhong Agricultural University,2014:31.

    [2] RENTSCH D,SCHMIDT S,TEGEDER M.Transporters for uptake and allocation of organic nitrogen compounds in plants [J].FebsLetters,2007,581(12):2281.

    [3] 馬 清,管 超,夏曾潤,等.高等植物氮素轉(zhuǎn)運(yùn)蛋白研究進(jìn)展[J].蘭州大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,51(2):222.

    MA Q,GUAN C,XIA Z R,etal.Membrane proteins mediating nitrogen transport in higher plants [J].JournalofLanzhouUniversity(NaturalSciences),2015,51(2):222.

    [4] FISCHER W N,ANDR B,RENTSCH D,etal.Plant amino acid transport [J].TrendsinPlantScience,1998,3(5):189.

    [5] 宋奇超,曹鳳秋,鞏元勇,等.高等植物氨基酸吸收與轉(zhuǎn)運(yùn)及生物學(xué)功能的研究進(jìn)展[J].植物營養(yǎng)與肥料學(xué)報(bào),2012,18(6):1509.

    SONG Q C,CAO F Q,GONG Y Y,etal.Current research progresses of amino acids uptake transport and their biological roles in higher plants [J].PlantNutritionandFertilizerScience,2012,18(6):1509.

    [6] RENTSCH D,HIRNER B,SCHMELZER E,etal.Salt stress-induced proline transporters and salt stress-repressed broad specificity amino acid permeases identified by suppression of a yeast amino acid permease-targeting mutant [J].ThePlantCell,1996,8(8):1437.

    [7] OKUMOTO S,KOCH W,TEGEDER M,etal.Root phloem-specific expression of the plasma membrane amino acid proton co-transporter AAP3 [J].JournalofExperimentalBotany,2004,55(406):2155.

    [8] HIMER A,LADWIG F,STRANSKY H,etal.ArabidopsisLHT1 is a high-affinity transporter for cellular amino acid uptake in both root epidermis and leaf mesophyll [J].ThePlantCell,2006,18(8):1931.

    [9] HAMMES U Z,NIELSEN E,HONAAS L A,etal.AtCAT6,a sink-tissue-localized transporter for essential amino acids inArabidopsis[J].PlantJournal,2006,48(3):415.

    [10] WIPF D,LUDEWIG U,TEGEDER M,etal.Conservation of amino acid transporters in fungi,plants and animals [J].TrendsinBiochemicalSciences,2002,27(3):141.

    [11] YANG H,POSTEL S,KEMMERLING B,etal.Altered growth and improved resistance ofArabidopsis,against pseudomonas syringae by overexpression of the basic amino acid transporterAtCAT1 [J].PlantCell&Environment,2014,37(6):1404.

    [12] YANG H,KREBS M,STIERHOF Y D,etal.Characterization of the putative amino acid transporter genesAtCAT2,3&4:the tonoplast localized AtCAT2 regulates soluble leaf amino acids [J].JournalofPlantPhysiology,2014,171(8):594.

    [13] LIU X,BUSH D R.Expression and transcriptional regulation of amino acid transporters in plants [J].AminoAcids,2006,30(2):115.

    [14] SU Y H,FROMMER W B,LUDEWIG U.Molecular and functional characterization of a family of amino acid transporters fromArabidopsis[J].PlantPhysiology,2004,136(2):3104.

    [15] YANG H,STIERHOF Y D,LUDEWIG U.The putativeCationicAminoAcidTransporter9 is targeted to vesicles and may be involved in plant amino acid homeostasis [J].FrontiersinPlantScience,2015,6:212.

    [16] SNOWDEN C J,BENJAMIN T,BAXTER C J,etal.A tonoplast Glu/Asp/GABA exchanger that affects tomato fruit amino acid composition [J].PlantJournal,2015,81(5):651.

    [17] 軒紅梅,王永華,魏利婷,等.小麥幼苗葉片中硝酸鹽轉(zhuǎn)運(yùn)蛋白NRT1和NRT2家族基因?qū)Φ囸I響應(yīng)的表達(dá)分析[J].麥類作物學(xué)報(bào),2014,34(8):1022.

    XUAN H M,WANG Y H,WEI L T,etal.Transcription analysis of the genes encoding nitrate transporter NRT1 and NRT2 families in response to nitrogen starvation in wheat seedlings leaves [J].JournalofTriticeaeCrops,2014,34(8):1022.

    [18] MARCUSSEN T,SANDVE S R,HEIER L,etal.Ancient hybridizations among the ancestral genomes of bread wheat [J].Science,2014,345(6194):1250092.

    [19] LIU Z L,ADAMS K L.Expression partitioning between genes duplicated by polyploidy under abiotic stress and during organ development [J].CurrentBiology,2007,17(19):1669.

    [20] HU Z,HAN Z,SONG N,etal.Epigenetic modification contributes to the expression divergence of threeTaEXPA1 homoeologs in hexaploid wheat(Triticumaestivum) [J].NewPhytologist,2013,197(4):1344.

    [21] SHITSUKAWA N,TAHIRA C,KASSAI K,etal.Genetic and epigenetic alteration among three homoeologous genes of a classEMADSbox gene in hexaploid wheat [J].ThePlantCell,2007,19(6):1728.

    [22] MAYER K F X,ROGERS J,DOLEEL J,etal.A chromsome-based draft sequence of the hexaploid bread wheat(Triticumaestivum) genome [J].Science,2014,345(6194):1251788.

    [23] 韓俊杰,王昊龍,李淼淼,等.小麥淀粉合成酶SSⅡa基因克隆及生物學(xué)分析[J].麥類作物學(xué)報(bào),2015,35(7):908.

    HAN J J,WANG H L,LI M M,etal.Cloning and biological analysis of wheatSSⅡagene [J].JournalofTriticeaeCrops,2015,35(7) :908.

    [24] 郎丹瑩,呂庚寅,梁廣旺,等.小麥UROS基因的可變剪接[J].麥類作物學(xué)報(bào),2016,36(12):1560.

    LANG D Y,Lü G Y,LIANG G W,etal,Alternative splicing of uroporphyrinogen Ⅲ synthase genes in wheat [J].JournalofTriticeaeCrops,2016,36(12):1560.

    猜你喜歡
    同源克隆氨基酸
    藥食同源
    ——紫 蘇
    兩岸年味連根同源
    克隆狼
    以同源詞看《詩經(jīng)》的訓(xùn)釋三則
    浙江:誕生首批體細(xì)胞克隆豬
    月桂酰丙氨基酸鈉的抑菌性能研究
    UFLC-QTRAP-MS/MS法同時(shí)測定絞股藍(lán)中11種氨基酸
    中成藥(2018年1期)2018-02-02 07:20:05
    抗BP5-KLH多克隆抗體的制備及鑒定
    虔誠書畫乃同源
    一株Nsp2蛋白自然缺失123個(gè)氨基酸的PRRSV分離和鑒定
    只有这里有精品99| 伦理电影免费视频| a 毛片基地| 涩涩av久久男人的天堂| 国产精品.久久久| 欧美日韩视频精品一区| 校园人妻丝袜中文字幕| 超碰97精品在线观看| 精品一区二区三卡| 大话2 男鬼变身卡| 一级a做视频免费观看| videossex国产| 特大巨黑吊av在线直播| 一个人看视频在线观看www免费| 男女啪啪激烈高潮av片| 啦啦啦在线观看免费高清www| 好男人视频免费观看在线| 亚洲精品中文字幕在线视频 | 在线观看人妻少妇| 久久国产精品大桥未久av | 少妇猛男粗大的猛烈进出视频| 51国产日韩欧美| 夜夜爽夜夜爽视频| av网站免费在线观看视频| 18禁裸乳无遮挡动漫免费视频| 欧美日韩国产mv在线观看视频 | 777米奇影视久久| 亚洲精品国产av成人精品| 亚洲精品一区蜜桃| 夜夜看夜夜爽夜夜摸| www.色视频.com| 国产免费视频播放在线视频| 九九爱精品视频在线观看| av视频免费观看在线观看| 亚洲成色77777| 日日啪夜夜撸| 国产精品国产av在线观看| 人妻 亚洲 视频| 国产成人午夜福利电影在线观看| 亚洲国产毛片av蜜桃av| 日韩,欧美,国产一区二区三区| 精品国产一区二区三区久久久樱花 | 蜜桃亚洲精品一区二区三区| 日韩欧美精品免费久久| 国产精品熟女久久久久浪| 免费黄网站久久成人精品| 欧美高清成人免费视频www| 国产在线视频一区二区| 成人亚洲欧美一区二区av| 欧美精品一区二区免费开放| 久久人人爽人人爽人人片va| 日本爱情动作片www.在线观看| 欧美xxxx黑人xx丫x性爽| 在现免费观看毛片| 久久久久国产精品人妻一区二区| 亚洲精品456在线播放app| 亚洲国产成人一精品久久久| 国产日韩欧美亚洲二区| 国产色婷婷99| 欧美老熟妇乱子伦牲交| 久久鲁丝午夜福利片| 高清午夜精品一区二区三区| 午夜福利网站1000一区二区三区| 大话2 男鬼变身卡| 亚洲真实伦在线观看| 18禁裸乳无遮挡动漫免费视频| 亚洲综合精品二区| 国产欧美日韩精品一区二区| 国产精品成人在线| 在线观看免费高清a一片| 午夜精品国产一区二区电影| 国产有黄有色有爽视频| 一级a做视频免费观看| 狂野欧美激情性bbbbbb| 国产探花极品一区二区| 日韩精品有码人妻一区| 一本一本综合久久| 在线播放无遮挡| 国产精品国产三级专区第一集| 噜噜噜噜噜久久久久久91| 麻豆精品久久久久久蜜桃| 国产伦理片在线播放av一区| 黄片wwwwww| 啦啦啦视频在线资源免费观看| 久久久久久人妻| 黄色日韩在线| 大片电影免费在线观看免费| 麻豆精品久久久久久蜜桃| 国产伦在线观看视频一区| 亚洲熟女精品中文字幕| videossex国产| 最近最新中文字幕大全电影3| 国产有黄有色有爽视频| 久久久国产一区二区| 毛片女人毛片| 男人爽女人下面视频在线观看| 精品人妻视频免费看| av国产免费在线观看| av又黄又爽大尺度在线免费看| 久久综合国产亚洲精品| 啦啦啦视频在线资源免费观看| 最近2019中文字幕mv第一页| 啦啦啦啦在线视频资源| 色吧在线观看| 久久人人爽av亚洲精品天堂 | 十八禁网站网址无遮挡 | 国精品久久久久久国模美| 精品一区二区三区视频在线| 亚洲精品国产av蜜桃| 成人高潮视频无遮挡免费网站| 久久久久精品久久久久真实原创| 精品少妇久久久久久888优播| 亚洲欧美清纯卡通| 亚洲人成网站在线观看播放| 欧美xxⅹ黑人| 2018国产大陆天天弄谢| 久久久久精品性色| 欧美人与善性xxx| 99热国产这里只有精品6| 王馨瑶露胸无遮挡在线观看| 大又大粗又爽又黄少妇毛片口| 天天躁夜夜躁狠狠久久av| 舔av片在线| 嘟嘟电影网在线观看| 亚洲欧美日韩东京热| 91午夜精品亚洲一区二区三区| 日韩av在线免费看完整版不卡| 九九久久精品国产亚洲av麻豆| 六月丁香七月| 午夜视频国产福利| 26uuu在线亚洲综合色| 大片免费播放器 马上看| 老女人水多毛片| 亚洲av日韩在线播放| 黑人猛操日本美女一级片| 久久精品国产亚洲网站| 亚洲av二区三区四区| 亚州av有码| 国产黄色视频一区二区在线观看| 国产精品.久久久| 久久久久视频综合| 欧美丝袜亚洲另类| 色网站视频免费| 高清av免费在线| 天堂俺去俺来也www色官网| 日本黄色日本黄色录像| 最近最新中文字幕免费大全7| 热99国产精品久久久久久7| 午夜免费男女啪啪视频观看| 久久精品久久久久久噜噜老黄| 性色avwww在线观看| 一级毛片aaaaaa免费看小| 久久99精品国语久久久| 最近手机中文字幕大全| 亚洲无线观看免费| 久久久久视频综合| 嘟嘟电影网在线观看| 精品国产一区二区三区久久久樱花 | 色视频在线一区二区三区| 熟妇人妻不卡中文字幕| 午夜福利网站1000一区二区三区| 国产av国产精品国产| videos熟女内射| 纵有疾风起免费观看全集完整版| 校园人妻丝袜中文字幕| av黄色大香蕉| 欧美亚洲 丝袜 人妻 在线| 国产亚洲午夜精品一区二区久久| 成人影院久久| 国产有黄有色有爽视频| 国产 一区精品| 五月玫瑰六月丁香| 久久久久网色| 伊人久久国产一区二区| 老司机影院毛片| 高清日韩中文字幕在线| 久久国产亚洲av麻豆专区| 国产精品久久久久成人av| 建设人人有责人人尽责人人享有的 | 熟女人妻精品中文字幕| 精品99又大又爽又粗少妇毛片| 久久婷婷青草| 欧美日韩综合久久久久久| 一区二区三区精品91| 一级片'在线观看视频| 一级毛片我不卡| 亚洲精品456在线播放app| 秋霞伦理黄片| 久久热精品热| 一区二区三区四区激情视频| av专区在线播放| 最黄视频免费看| 日本欧美国产在线视频| 亚洲色图av天堂| 免费久久久久久久精品成人欧美视频 | 亚洲欧美成人精品一区二区| 美女国产视频在线观看| 亚洲精品中文字幕在线视频 | 久久久久久九九精品二区国产| 国产综合精华液| 韩国高清视频一区二区三区| 久久亚洲国产成人精品v| 日韩,欧美,国产一区二区三区| 丝瓜视频免费看黄片| 国产伦精品一区二区三区四那| 欧美日韩一区二区视频在线观看视频在线| 国产午夜精品久久久久久一区二区三区| 视频区图区小说| 一区二区三区四区激情视频| 如何舔出高潮| 久久久久久久久久成人| 欧美精品一区二区免费开放| 国产高清国产精品国产三级 | 99re6热这里在线精品视频| 欧美亚洲 丝袜 人妻 在线| 国产精品一区二区在线不卡| 一级黄片播放器| 久久婷婷青草| 国产乱人视频| 激情五月婷婷亚洲| 国产精品一二三区在线看| 国产精品人妻久久久久久| 亚洲精品,欧美精品| 久久久国产一区二区| 国产在视频线精品| 免费黄频网站在线观看国产| 日韩人妻高清精品专区| 婷婷色综合www| 校园人妻丝袜中文字幕| 欧美97在线视频| 精品久久久久久久久亚洲| 国产老妇伦熟女老妇高清| 亚洲国产日韩一区二区| 久久ye,这里只有精品| 内地一区二区视频在线| 日本av免费视频播放| 老司机影院成人| 精品久久久噜噜| 亚洲av欧美aⅴ国产| 91久久精品国产一区二区成人| 国产淫语在线视频| 亚洲av.av天堂| 久久99热6这里只有精品| 日本-黄色视频高清免费观看| 午夜免费鲁丝| 免费看不卡的av| 国产精品国产三级国产av玫瑰| 啦啦啦啦在线视频资源| 亚洲一区二区三区欧美精品| 麻豆成人午夜福利视频| 亚洲第一区二区三区不卡| 免费观看无遮挡的男女| 久久毛片免费看一区二区三区| 日韩制服骚丝袜av| 狠狠精品人妻久久久久久综合| 久久精品人妻少妇| 亚洲欧美中文字幕日韩二区| 欧美三级亚洲精品| 啦啦啦视频在线资源免费观看| 久久毛片免费看一区二区三区| 国产精品国产三级国产av玫瑰| 特大巨黑吊av在线直播| 午夜日本视频在线| 中文乱码字字幕精品一区二区三区| 三级经典国产精品| 久久ye,这里只有精品| 国产成人a区在线观看| 久久久久精品久久久久真实原创| 国产v大片淫在线免费观看| 中文字幕av成人在线电影| 最近手机中文字幕大全| 大陆偷拍与自拍| 天堂8中文在线网| 在线观看免费视频网站a站| 亚洲国产精品专区欧美| 久久久久久人妻| 亚洲美女黄色视频免费看| 国产片特级美女逼逼视频| 成人免费观看视频高清| 亚洲色图综合在线观看| 国产精品不卡视频一区二区| 九色成人免费人妻av| 中文在线观看免费www的网站| 精品久久久久久久久亚洲| 一级黄片播放器| 我的老师免费观看完整版| 丰满少妇做爰视频| 黄色日韩在线| 久久久久久久精品精品| 国产成人一区二区在线| 91久久精品国产一区二区三区| 丝瓜视频免费看黄片| 久久国产乱子免费精品| 免费看av在线观看网站| 亚洲美女黄色视频免费看| 精品人妻熟女av久视频| 欧美精品人与动牲交sv欧美| 简卡轻食公司| 亚洲精品aⅴ在线观看| 一本久久精品| 一级毛片黄色毛片免费观看视频| 成人无遮挡网站| 亚洲久久久国产精品| 黑人猛操日本美女一级片| 成人综合一区亚洲| 中文精品一卡2卡3卡4更新| 欧美成人午夜免费资源| 在线观看一区二区三区| 亚洲精品国产av蜜桃| 亚洲欧美日韩无卡精品| 免费看av在线观看网站| 少妇人妻一区二区三区视频| 国产一级毛片在线| 久久国产乱子免费精品| 如何舔出高潮| 99热国产这里只有精品6| 成人二区视频| 黄片无遮挡物在线观看| 久久99热6这里只有精品| 男女啪啪激烈高潮av片| 亚洲精品成人av观看孕妇| 1000部很黄的大片| 男人和女人高潮做爰伦理| 黄色日韩在线| 一级毛片我不卡| 久久女婷五月综合色啪小说| 国产免费一区二区三区四区乱码| 亚洲成色77777| 少妇裸体淫交视频免费看高清| 国产午夜精品一二区理论片| 欧美少妇被猛烈插入视频| 亚洲三级黄色毛片| 日本猛色少妇xxxxx猛交久久| 精品久久久噜噜| 爱豆传媒免费全集在线观看| 欧美一级a爱片免费观看看| 久久久久久伊人网av| 直男gayav资源| 成人国产麻豆网| 晚上一个人看的免费电影| 亚洲美女黄色视频免费看| 99热这里只有是精品在线观看| 久久精品人妻少妇| 国产乱来视频区| 亚洲av中文av极速乱| 51国产日韩欧美| 天天躁夜夜躁狠狠久久av| 久久久久精品性色| 亚洲一区二区三区欧美精品| 青春草国产在线视频| 国产黄色免费在线视频| 天堂俺去俺来也www色官网| 午夜福利在线在线| 久久99热6这里只有精品| 国产精品国产三级专区第一集| 99热6这里只有精品| 午夜日本视频在线| 久久精品久久久久久久性| 婷婷色麻豆天堂久久| 狂野欧美白嫩少妇大欣赏| 亚洲精品aⅴ在线观看| 在现免费观看毛片| 日韩欧美 国产精品| 99久国产av精品国产电影| 国产av码专区亚洲av| 午夜老司机福利剧场| 欧美日韩综合久久久久久| 国产免费一区二区三区四区乱码| 欧美日韩综合久久久久久| 天天躁日日操中文字幕| av天堂中文字幕网| 99久久综合免费| 99re6热这里在线精品视频| 午夜免费鲁丝| 亚洲色图av天堂| 嫩草影院新地址| 久久久久久久久久成人| 国产免费视频播放在线视频| 大陆偷拍与自拍| 美女福利国产在线 | 亚洲一区二区三区欧美精品| 男男h啪啪无遮挡| 国产亚洲av片在线观看秒播厂| 久久女婷五月综合色啪小说| 亚洲欧洲国产日韩| 干丝袜人妻中文字幕| 在线观看av片永久免费下载| 免费大片18禁| videos熟女内射| 街头女战士在线观看网站| 国产综合精华液| 国产又色又爽无遮挡免| 亚洲国产欧美在线一区| 99热6这里只有精品| 九九在线视频观看精品| 中国国产av一级| 高清毛片免费看| 日韩精品有码人妻一区| 99热全是精品| 制服丝袜香蕉在线| 女人十人毛片免费观看3o分钟| 男人添女人高潮全过程视频| 国产精品国产三级专区第一集| 国产精品熟女久久久久浪| 不卡视频在线观看欧美| 2018国产大陆天天弄谢| 国产91av在线免费观看| 国产高潮美女av| 国产免费视频播放在线视频| 老女人水多毛片| av播播在线观看一区| 日韩免费高清中文字幕av| 内射极品少妇av片p| 国产伦精品一区二区三区视频9| 免费在线观看成人毛片| 亚洲第一区二区三区不卡| 少妇猛男粗大的猛烈进出视频| 欧美老熟妇乱子伦牲交| 国产男人的电影天堂91| 亚洲精品乱码久久久v下载方式| 免费高清在线观看视频在线观看| 国产片特级美女逼逼视频| 一个人看的www免费观看视频| 国产高清国产精品国产三级 | 亚洲国产精品专区欧美| 免费看日本二区| 中文字幕av成人在线电影| 久久午夜福利片| 不卡视频在线观看欧美| 久久久精品94久久精品| 又黄又爽又刺激的免费视频.| 日韩伦理黄色片| 99热全是精品| 国产高潮美女av| 亚洲精品国产av成人精品| 精品酒店卫生间| 中国国产av一级| 久久婷婷青草| 自拍欧美九色日韩亚洲蝌蚪91 | 国产黄色免费在线视频| 欧美成人精品欧美一级黄| 国产精品秋霞免费鲁丝片| 国产色爽女视频免费观看| 一级a做视频免费观看| 亚洲国产精品999| 91在线精品国自产拍蜜月| 网址你懂的国产日韩在线| 最后的刺客免费高清国语| 午夜福利在线观看免费完整高清在| 51国产日韩欧美| 亚洲欧美中文字幕日韩二区| 中国美白少妇内射xxxbb| 少妇的逼好多水| 少妇精品久久久久久久| av国产免费在线观看| 免费看av在线观看网站| 欧美精品人与动牲交sv欧美| 一级爰片在线观看| 亚洲精品色激情综合| av线在线观看网站| 一区在线观看完整版| 成年免费大片在线观看| 久久久亚洲精品成人影院| 色婷婷久久久亚洲欧美| 交换朋友夫妻互换小说| 美女xxoo啪啪120秒动态图| 国产午夜精品一二区理论片| 国产欧美另类精品又又久久亚洲欧美| 日本黄色片子视频| 黄片无遮挡物在线观看| 99久久精品国产国产毛片| 精品国产露脸久久av麻豆| 婷婷色综合www| 1000部很黄的大片| 国产老妇伦熟女老妇高清| 欧美xxⅹ黑人| 国产亚洲精品久久久com| 亚洲欧美精品自产自拍| 五月天丁香电影| 午夜福利网站1000一区二区三区| 熟女人妻精品中文字幕| 这个男人来自地球电影免费观看 | 国产精品99久久久久久久久| 夫妻午夜视频| 久久久久久久亚洲中文字幕| 亚洲自偷自拍三级| av专区在线播放| 搡老乐熟女国产| 高清视频免费观看一区二区| 久久97久久精品| 大陆偷拍与自拍| 最后的刺客免费高清国语| 成人免费观看视频高清| 一区二区三区免费毛片| 国产av精品麻豆| 女人十人毛片免费观看3o分钟| 国产在线视频一区二区| 狂野欧美激情性xxxx在线观看| 秋霞伦理黄片| 国产精品免费大片| 一级毛片aaaaaa免费看小| 九色成人免费人妻av| 免费人成在线观看视频色| 久久人人爽人人爽人人片va| 大又大粗又爽又黄少妇毛片口| 欧美精品国产亚洲| 国产色婷婷99| 国产国拍精品亚洲av在线观看| 联通29元200g的流量卡| 国产精品免费大片| 久久久精品免费免费高清| 2018国产大陆天天弄谢| 久久精品熟女亚洲av麻豆精品| 亚洲自偷自拍三级| 人妻一区二区av| 一级毛片我不卡| freevideosex欧美| 国产乱人偷精品视频| 丰满乱子伦码专区| 国产精品.久久久| 亚洲欧美中文字幕日韩二区| av国产免费在线观看| 日本欧美视频一区| 在线观看美女被高潮喷水网站| 熟女电影av网| 99热这里只有是精品在线观看| 亚洲国产精品专区欧美| kizo精华| 一本久久精品| 亚洲美女黄色视频免费看| 国产亚洲av片在线观看秒播厂| 久久久久人妻精品一区果冻| 大陆偷拍与自拍| 日韩一区二区三区影片| 成人国产av品久久久| 亚洲精品456在线播放app| 国产一区亚洲一区在线观看| 一个人看视频在线观看www免费| 亚洲av免费高清在线观看| 精品人妻视频免费看| 国产免费又黄又爽又色| 日本av手机在线免费观看| 免费观看av网站的网址| 免费看日本二区| 又粗又硬又长又爽又黄的视频| 久久鲁丝午夜福利片| 在线免费十八禁| 婷婷色综合大香蕉| 国产色爽女视频免费观看| 视频中文字幕在线观看| 高清黄色对白视频在线免费看 | 女人久久www免费人成看片| 国产精品久久久久久精品电影小说 | 日本欧美国产在线视频| av免费观看日本| 久久精品熟女亚洲av麻豆精品| 亚洲aⅴ乱码一区二区在线播放| 国产爱豆传媒在线观看| 九九爱精品视频在线观看| 精品一区在线观看国产| 插逼视频在线观看| 亚洲精品中文字幕在线视频 | 少妇裸体淫交视频免费看高清| 国产成人精品久久久久久| 国产v大片淫在线免费观看| 少妇的逼好多水| 亚洲av二区三区四区| 麻豆乱淫一区二区| 18禁裸乳无遮挡动漫免费视频| 欧美+日韩+精品| 女性被躁到高潮视频| 又大又黄又爽视频免费| 国产成人午夜福利电影在线观看| 一本久久精品| 亚洲成人中文字幕在线播放| 日韩欧美一区视频在线观看 | 高清在线视频一区二区三区| 偷拍熟女少妇极品色| 三级国产精品欧美在线观看| 成年av动漫网址| 我要看日韩黄色一级片| 亚洲欧美精品自产自拍| 狠狠精品人妻久久久久久综合| 91久久精品国产一区二区三区| 全区人妻精品视频| 国产一级毛片在线| 欧美国产精品一级二级三级 | 国产欧美另类精品又又久久亚洲欧美| 久久鲁丝午夜福利片| 最近的中文字幕免费完整| 最近最新中文字幕免费大全7| 在线精品无人区一区二区三 | 26uuu在线亚洲综合色| 午夜免费男女啪啪视频观看| 美女cb高潮喷水在线观看| av国产免费在线观看| 亚洲人成网站在线观看播放| 日本一二三区视频观看| 日韩伦理黄色片| 91精品一卡2卡3卡4卡| 国产精品久久久久成人av| 少妇丰满av| 日韩成人伦理影院| 在现免费观看毛片| 日本wwww免费看| 欧美+日韩+精品| 最近2019中文字幕mv第一页| 中国三级夫妇交换| 精品酒店卫生间| 国产精品一二三区在线看| 高清av免费在线| 欧美区成人在线视频| 美女视频免费永久观看网站| 一级毛片电影观看| 精品熟女少妇av免费看|