• <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分離和鑒定
    我要看黄色一级片免费的| 国产亚洲5aaaaa淫片| 亚洲伊人久久精品综合| 综合色丁香网| 狂野欧美白嫩少妇大欣赏| 中文字幕人妻熟人妻熟丝袜美| 日本黄色日本黄色录像| av播播在线观看一区| 啦啦啦在线观看免费高清www| 久久午夜福利片| 国产精品蜜桃在线观看| 国产欧美日韩一区二区三区在线 | 少妇 在线观看| 亚洲av成人精品一区久久| kizo精华| 亚洲欧美精品专区久久| 精品少妇久久久久久888优播| 女性被躁到高潮视频| 国产在线一区二区三区精| 欧美激情国产日韩精品一区| 久久国产乱子免费精品| 亚洲av电影在线观看一区二区三区| 中文乱码字字幕精品一区二区三区| 新久久久久国产一级毛片| 少妇高潮的动态图| 成人二区视频| 国产精品嫩草影院av在线观看| 国产有黄有色有爽视频| 高清视频免费观看一区二区| 亚洲国产精品999| kizo精华| 女人久久www免费人成看片| av天堂中文字幕网| 91狼人影院| 日日摸夜夜添夜夜添av毛片| 亚洲性久久影院| 看十八女毛片水多多多| 亚洲国产高清在线一区二区三| 国产极品天堂在线| 亚洲美女视频黄频| 亚洲精品乱久久久久久| 伦理电影免费视频| 色婷婷久久久亚洲欧美| 97在线人人人人妻| 欧美人与善性xxx| 久久久久久久久久人人人人人人| 少妇人妻久久综合中文| 午夜老司机福利剧场| 亚洲国产成人一精品久久久| 成人高潮视频无遮挡免费网站| 国产亚洲5aaaaa淫片| 成人18禁高潮啪啪吃奶动态图 | 国产成人91sexporn| 天天躁日日操中文字幕| 精品亚洲成国产av| av.在线天堂| 我的女老师完整版在线观看| 国产淫片久久久久久久久| 国产成人a区在线观看| 一级毛片黄色毛片免费观看视频| 2018国产大陆天天弄谢| 久久精品国产亚洲网站| 蜜桃在线观看..| 亚洲精品乱码久久久久久按摩| 一级a做视频免费观看| av女优亚洲男人天堂| 精品一区二区三卡| 日日啪夜夜撸| 久久97久久精品| 国产精品久久久久久久久免| 看免费成人av毛片| 黄片wwwwww| 精品少妇久久久久久888优播| 香蕉精品网在线| 日韩国内少妇激情av| 成人黄色视频免费在线看| 色视频在线一区二区三区| 亚洲av综合色区一区| 人妻一区二区av| 亚洲精品,欧美精品| 一级二级三级毛片免费看| 国产精品一二三区在线看| 国产伦理片在线播放av一区| 亚洲欧美成人综合另类久久久| 卡戴珊不雅视频在线播放| 欧美高清性xxxxhd video| 亚洲成色77777| 日日撸夜夜添| 亚洲精品第二区| 97超碰精品成人国产| 国产午夜精品久久久久久一区二区三区| 欧美一区二区亚洲| av线在线观看网站| 国产免费一级a男人的天堂| 国产精品国产三级国产av玫瑰| 日本爱情动作片www.在线观看| 日韩成人av中文字幕在线观看| 丰满迷人的少妇在线观看| 一二三四中文在线观看免费高清| 高清欧美精品videossex| 国模一区二区三区四区视频| 麻豆成人午夜福利视频| 香蕉精品网在线| 免费人成在线观看视频色| 欧美激情国产日韩精品一区| 在线天堂最新版资源| 亚洲av中文av极速乱| 精品久久久久久电影网| av一本久久久久| 日本一二三区视频观看| h视频一区二区三区| 亚洲av欧美aⅴ国产| 免费av中文字幕在线| 韩国av在线不卡| 高清视频免费观看一区二区| 亚洲欧洲国产日韩| 亚洲国产色片| 在线播放无遮挡| 美女主播在线视频| 一本一本综合久久| 中国美白少妇内射xxxbb| 国产一区二区三区av在线| 久久99热这里只频精品6学生| 18禁动态无遮挡网站| 国产精品一区二区在线观看99| 高清日韩中文字幕在线| 熟女人妻精品中文字幕| 国产白丝娇喘喷水9色精品| 色婷婷av一区二区三区视频| 老司机影院成人| 人妻 亚洲 视频| 亚洲欧美日韩无卡精品| 午夜福利视频精品| 99热这里只有精品一区| 少妇精品久久久久久久| 成年女人在线观看亚洲视频| 成人免费观看视频高清| 精品酒店卫生间| 91午夜精品亚洲一区二区三区| 久久久精品94久久精品| 日本vs欧美在线观看视频 | 精品少妇久久久久久888优播| 韩国av在线不卡| 18+在线观看网站| 欧美精品人与动牲交sv欧美| 国产亚洲91精品色在线| 久久久久网色| 欧美一级a爱片免费观看看| 午夜福利网站1000一区二区三区| 中国三级夫妇交换| 午夜福利影视在线免费观看| 蜜桃亚洲精品一区二区三区| 免费少妇av软件| 国产精品一区二区三区四区免费观看| 蜜臀久久99精品久久宅男| 国国产精品蜜臀av免费| 中文乱码字字幕精品一区二区三区| 成人一区二区视频在线观看| 寂寞人妻少妇视频99o| 在线观看三级黄色| 日韩 亚洲 欧美在线| 免费人成在线观看视频色| 天堂中文最新版在线下载| 亚洲国产欧美在线一区| 日韩免费高清中文字幕av| 晚上一个人看的免费电影| 免费黄频网站在线观看国产| 国产在线免费精品| 亚洲精品日韩在线中文字幕| 国产高清国产精品国产三级 | 免费观看av网站的网址| 色婷婷av一区二区三区视频| 91在线精品国自产拍蜜月| 热re99久久精品国产66热6| 少妇的逼水好多| 91午夜精品亚洲一区二区三区| 日韩av免费高清视频| 18禁裸乳无遮挡免费网站照片| 国产91av在线免费观看| 看免费成人av毛片| 一区在线观看完整版| 99久久综合免费| 久久精品国产a三级三级三级| 熟女电影av网| av又黄又爽大尺度在线免费看| 中文乱码字字幕精品一区二区三区| 伊人久久国产一区二区| 国产男女超爽视频在线观看| 亚洲av成人精品一二三区| 免费观看无遮挡的男女| 九九爱精品视频在线观看| 国产av码专区亚洲av| 我的女老师完整版在线观看| 日本欧美国产在线视频| 久久精品国产亚洲av天美| 一个人看视频在线观看www免费| 大话2 男鬼变身卡| 日日撸夜夜添| 久久韩国三级中文字幕| 成人二区视频| 久久亚洲国产成人精品v| 成人毛片a级毛片在线播放| 成人国产av品久久久| 久久97久久精品| 亚洲精品国产av蜜桃| 日韩国内少妇激情av| 日韩中文字幕视频在线看片 | 日本wwww免费看| 成年人午夜在线观看视频| 久久鲁丝午夜福利片| 国产一区亚洲一区在线观看| 在线亚洲精品国产二区图片欧美 | 欧美日韩综合久久久久久| 男女国产视频网站| 夜夜爽夜夜爽视频| 在线观看一区二区三区| 国产精品国产三级国产专区5o| 亚洲成人av在线免费| a级毛片免费高清观看在线播放| 中文字幕亚洲精品专区| 亚洲精品自拍成人| 国内精品宾馆在线| 国产欧美亚洲国产| 一级二级三级毛片免费看| 久久久a久久爽久久v久久| 一区二区三区精品91| 高清欧美精品videossex| 久久精品国产亚洲av涩爱| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品国产av蜜桃| 22中文网久久字幕| 夜夜骑夜夜射夜夜干| 久久精品国产亚洲av涩爱| 2018国产大陆天天弄谢| 国产色爽女视频免费观看| 国产 一区精品| 在现免费观看毛片| 波野结衣二区三区在线| 校园人妻丝袜中文字幕| 另类亚洲欧美激情| 2018国产大陆天天弄谢| 亚洲国产毛片av蜜桃av| 国产成人91sexporn| 免费看av在线观看网站| 日韩av在线免费看完整版不卡| 97热精品久久久久久| 97超碰精品成人国产| 爱豆传媒免费全集在线观看| 蜜桃久久精品国产亚洲av| 国产乱人视频| 中文乱码字字幕精品一区二区三区| 女性生殖器流出的白浆| 天堂8中文在线网| 成人美女网站在线观看视频| 99久久精品热视频| 极品教师在线视频| 日产精品乱码卡一卡2卡三| 成人无遮挡网站| 91狼人影院| 高清av免费在线| 国产在线一区二区三区精| 亚洲av国产av综合av卡| 精品一区二区三区视频在线| 精品熟女少妇av免费看| 免费观看无遮挡的男女| 免费看不卡的av| 日韩大片免费观看网站| 国产深夜福利视频在线观看| 自拍偷自拍亚洲精品老妇| 国产精品嫩草影院av在线观看| 久久久久久久久久久免费av| 色视频在线一区二区三区| 久久这里有精品视频免费| 99久久精品热视频| 亚洲国产精品成人久久小说| 大陆偷拍与自拍| 久久精品久久精品一区二区三区| 亚洲,欧美,日韩| 国产成人精品婷婷| 日本欧美视频一区| 爱豆传媒免费全集在线观看| av在线app专区| 国产黄频视频在线观看| 韩国av在线不卡| 1000部很黄的大片| 免费观看的影片在线观看| 人妻制服诱惑在线中文字幕| 欧美一区二区亚洲| 在线观看免费视频网站a站| 精品一区二区免费观看| 久久精品国产亚洲av天美| 91aial.com中文字幕在线观看| 亚洲国产精品成人久久小说| 高清黄色对白视频在线免费看 | 亚洲av在线观看美女高潮| 日本与韩国留学比较| 精品一区二区三区视频在线| 国产熟女欧美一区二区| 好男人视频免费观看在线| 亚洲欧美日韩东京热| 久久影院123| 成人美女网站在线观看视频| 亚洲自偷自拍三级| 一级爰片在线观看| 亚洲av成人精品一二三区| 免费看日本二区| 日韩av在线免费看完整版不卡| 热99国产精品久久久久久7| 国产高清不卡午夜福利| av国产久精品久网站免费入址| 欧美激情国产日韩精品一区| 国产精品女同一区二区软件| 日韩一区二区视频免费看| 国产一区亚洲一区在线观看| 一级爰片在线观看| 国产一区亚洲一区在线观看| 涩涩av久久男人的天堂| 十分钟在线观看高清视频www | 亚洲第一av免费看| 一本色道久久久久久精品综合| 欧美成人精品欧美一级黄| 国产黄色免费在线视频| 水蜜桃什么品种好| 亚洲精品乱码久久久久久按摩| 中文乱码字字幕精品一区二区三区| 九九在线视频观看精品| 欧美日韩国产mv在线观看视频 | 热99国产精品久久久久久7| 久久久久久久久久人人人人人人| 国模一区二区三区四区视频| 美女脱内裤让男人舔精品视频| 少妇猛男粗大的猛烈进出视频| 久久久久久久久久成人| 久久人人爽av亚洲精品天堂 | 欧美日韩视频高清一区二区三区二| 亚洲av中文字字幕乱码综合| 有码 亚洲区| 国产一区有黄有色的免费视频| 国产成人91sexporn| 欧美成人午夜免费资源| 精品一区二区免费观看| videos熟女内射| 国内揄拍国产精品人妻在线| 亚洲av不卡在线观看| 久久99蜜桃精品久久| a级毛色黄片| 国产极品天堂在线| 热re99久久精品国产66热6| 国产亚洲5aaaaa淫片| 亚洲精品成人av观看孕妇| 1000部很黄的大片| 国产精品一区www在线观看| 午夜福利在线观看免费完整高清在| 亚洲综合精品二区| 你懂的网址亚洲精品在线观看| 欧美成人a在线观看| 久久久久国产网址| 99久久综合免费| 黄色日韩在线| 国模一区二区三区四区视频| 国产精品久久久久久久电影| 在现免费观看毛片| 国内少妇人妻偷人精品xxx网站| 日本欧美国产在线视频| 成人毛片a级毛片在线播放| 精品一区二区三区视频在线| 久久久亚洲精品成人影院| 在线 av 中文字幕| 成人毛片a级毛片在线播放| 亚洲国产色片| 亚洲av国产av综合av卡| 韩国av在线不卡| 国产亚洲5aaaaa淫片| 三级国产精品片| 国产亚洲5aaaaa淫片| 三级国产精品片| 日本与韩国留学比较| 免费看av在线观看网站| 精品国产乱码久久久久久小说| 亚洲成人一二三区av| 国产一区二区三区av在线| 亚洲精品456在线播放app| 亚洲美女黄色视频免费看| 男人狂女人下面高潮的视频| 国内揄拍国产精品人妻在线| videos熟女内射| 97超视频在线观看视频| 在线观看免费视频网站a站| 国产亚洲午夜精品一区二区久久| 另类亚洲欧美激情| 观看美女的网站| 亚洲最大成人中文| 午夜福利视频精品| 亚洲国产高清在线一区二区三| 2021少妇久久久久久久久久久| 国产成人精品久久久久久| 大香蕉97超碰在线| 91久久精品电影网| 亚洲精品第二区| 日本与韩国留学比较| 一级av片app| 一区在线观看完整版| 亚洲人成网站在线播| 国内揄拍国产精品人妻在线| 国产亚洲av片在线观看秒播厂| 日韩中文字幕视频在线看片 | 18禁在线播放成人免费| 不卡视频在线观看欧美| 国产免费视频播放在线视频| 久久久久久久久久人人人人人人| 在现免费观看毛片| 亚洲av中文av极速乱| 色视频www国产| 国产精品无大码| 午夜精品国产一区二区电影| 99久久精品热视频| 国产成人a区在线观看| 精品久久久久久久久av| kizo精华| 免费观看性生交大片5| www.av在线官网国产| 亚洲高清免费不卡视频| 免费观看在线日韩| 久久精品夜色国产| 一区在线观看完整版| 亚洲不卡免费看| 色吧在线观看| 免费人成在线观看视频色| 国产视频首页在线观看| 内射极品少妇av片p| 久久国产亚洲av麻豆专区| 国产av精品麻豆| 春色校园在线视频观看| 日韩av不卡免费在线播放| 麻豆乱淫一区二区| 美女主播在线视频| 国产 一区 欧美 日韩| 少妇的逼好多水| 在线播放无遮挡| 91在线精品国自产拍蜜月| 国产在线视频一区二区| 超碰av人人做人人爽久久| 国产精品人妻久久久久久| 丰满人妻一区二区三区视频av| 国产午夜精品久久久久久一区二区三区| 国产黄频视频在线观看| 高清视频免费观看一区二区| 一级毛片 在线播放| 又黄又爽又刺激的免费视频.| 日韩一区二区三区影片| 1000部很黄的大片| 亚洲精品国产av成人精品| 免费看不卡的av| av一本久久久久| 又黄又爽又刺激的免费视频.| 国产一区亚洲一区在线观看| 欧美一级a爱片免费观看看| 美女主播在线视频| 秋霞在线观看毛片| 免费人妻精品一区二区三区视频| 国产乱来视频区| 久久99蜜桃精品久久| 亚洲成人一二三区av| 日本午夜av视频| 国产精品一二三区在线看| 在线观看av片永久免费下载| 五月天丁香电影| 高清黄色对白视频在线免费看 | 丰满乱子伦码专区| 秋霞伦理黄片| 黑人高潮一二区| 大又大粗又爽又黄少妇毛片口| 亚洲av福利一区| 秋霞伦理黄片| 各种免费的搞黄视频| 国产高清三级在线| 国产精品秋霞免费鲁丝片| 久久精品国产鲁丝片午夜精品| 亚洲怡红院男人天堂| 极品少妇高潮喷水抽搐| 看非洲黑人一级黄片| 亚洲成人一二三区av| 黄片无遮挡物在线观看| 亚洲精品第二区| 在线观看av片永久免费下载| 在线免费十八禁| 国产熟女欧美一区二区| 丰满少妇做爰视频| 黄色日韩在线| 免费观看无遮挡的男女| 黄色一级大片看看| 亚洲欧美日韩卡通动漫| 亚洲国产日韩一区二区| 久久97久久精品| 黄色配什么色好看| 国产高清有码在线观看视频| 国产探花极品一区二区| 一二三四中文在线观看免费高清| 亚洲国产日韩一区二区| 久久97久久精品| 国产精品一区二区三区四区免费观看| 六月丁香七月| 国产av码专区亚洲av| 国语对白做爰xxxⅹ性视频网站| 日本免费在线观看一区| 国内精品宾馆在线| 永久免费av网站大全| 在线免费十八禁| 这个男人来自地球电影免费观看 | 日本-黄色视频高清免费观看| 免费观看的影片在线观看| 国产精品一区二区三区四区免费观看| 国产一区亚洲一区在线观看| 伊人久久国产一区二区| 日韩中字成人| 97精品久久久久久久久久精品| 女人久久www免费人成看片| 欧美少妇被猛烈插入视频| 精品亚洲成国产av| 国产精品国产三级国产av玫瑰| 亚洲人成网站在线播| 亚洲欧美成人综合另类久久久| 如何舔出高潮| 国产精品无大码| 国产爽快片一区二区三区| 少妇丰满av| 91精品国产九色| 免费播放大片免费观看视频在线观看| 国产一级毛片在线| 国产极品天堂在线| 久久人人爽人人片av| 日韩人妻高清精品专区| a 毛片基地| 日韩人妻高清精品专区| 久久精品国产亚洲网站| 国产午夜精品久久久久久一区二区三区| 国产黄色视频一区二区在线观看| 国产高潮美女av| 久久精品夜色国产| 久久精品人妻少妇| 97精品久久久久久久久久精品| 国内少妇人妻偷人精品xxx网站| 国产精品国产三级国产专区5o| av免费观看日本| 嘟嘟电影网在线观看| 精品国产三级普通话版| 99热这里只有是精品50| 大码成人一级视频| videos熟女内射| 爱豆传媒免费全集在线观看| 青春草视频在线免费观看| 中国三级夫妇交换| 精品一区二区三卡| 91久久精品国产一区二区三区| a级一级毛片免费在线观看| 精品一区在线观看国产| 国产欧美日韩精品一区二区| 日本猛色少妇xxxxx猛交久久| 欧美日本视频| 色视频在线一区二区三区| 三级经典国产精品| 成人无遮挡网站| 午夜福利在线观看免费完整高清在| 国产老妇伦熟女老妇高清| 精品人妻偷拍中文字幕| 日韩一区二区三区影片| 欧美bdsm另类| 亚洲国产精品国产精品| 久久久久久伊人网av| 精品久久久噜噜| 一边亲一边摸免费视频| 精品视频人人做人人爽| 国产亚洲欧美精品永久| 秋霞伦理黄片| 在线观看国产h片| 国产成人一区二区在线| 一级毛片久久久久久久久女| 成人亚洲精品一区在线观看 | 免费看日本二区| 少妇被粗大猛烈的视频| 91精品国产九色| 精品国产乱码久久久久久小说| 内地一区二区视频在线| 亚洲综合色惰| tube8黄色片| 国产国拍精品亚洲av在线观看| 日韩电影二区| 国产亚洲最大av| 久久毛片免费看一区二区三区| 在线观看美女被高潮喷水网站| 亚洲中文av在线| 欧美日韩国产mv在线观看视频 | 老司机影院成人| 一区二区三区精品91| 内地一区二区视频在线| 国产精品人妻久久久影院| 蜜桃在线观看..| 亚洲欧美日韩卡通动漫| 亚洲国产高清在线一区二区三| 狂野欧美激情性bbbbbb| 黄片无遮挡物在线观看| 亚洲图色成人| 自拍欧美九色日韩亚洲蝌蚪91 | 男人添女人高潮全过程视频| 综合色丁香网| 全区人妻精品视频| 少妇熟女欧美另类| 日韩av不卡免费在线播放| 蜜桃亚洲精品一区二区三区| 日日撸夜夜添| 18+在线观看网站| 欧美成人a在线观看| 日韩欧美 国产精品| 嘟嘟电影网在线观看| 少妇人妻久久综合中文|