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

    福州野生蕉隱花色素和光敏色素基因家族的克隆及其表達(dá)分析

    2017-05-30 06:59:08劉轉(zhuǎn)霞余運(yùn)康陳裕坤馮新劉煒?gòu)O張梓浩程春振林玉玲賴鐘雄
    熱帶作物學(xué)報(bào) 2017年11期
    關(guān)鍵詞:表達(dá)分析

    劉轉(zhuǎn)霞 余運(yùn)康 陳裕坤 馮新 劉煒?gòu)O 張梓浩 程春振 林玉玲 賴鐘雄

    摘 要 采用RT-PCR結(jié)合RACE法從福州野生蕉試管苗中克隆隱花色素和光敏色素基因家族成員的cDNA序列,命名為MuCry1、MuCry2a、MuCry2b、MuPhyB、MuPhyC1。MuCrys和MuPhys的cDNA序列依次為2 330、2 859、2 752、3 272、3 912 bp,編碼698、666、669、1 089、1 003個(gè)氨基酸。生物信息學(xué)分析表明:MuCrys和MuPhys均屬不穩(wěn)定、親水蛋白,具跨膜結(jié)構(gòu)域和卷曲螺旋結(jié)構(gòu),都不具信號(hào)肽,均定位于細(xì)胞核。系統(tǒng)進(jìn)化樹分析結(jié)果表明:Cry1和Cry2聚為兩大類,其中MuCry1與小果野蕉、海棗和油棕Cry1的親緣關(guān)系最近,與雙子葉植物大豆Cry1的親緣關(guān)系最遠(yuǎn);MuCry2a和MuCry2b與小果野蕉、水稻和小麥Cry2的親緣關(guān)系最近,與雙子葉植物碧桃Cry2的親緣關(guān)系最遠(yuǎn);MuPhyB與小果野蕉PhyB的親緣關(guān)系最近,與雙子葉植物毛果楊PhyB的親緣關(guān)系最遠(yuǎn);MuPhyC1與小果野蕉和馬來蘭花蕉PhyC的親緣關(guān)系最近,與雙子葉甜橙PhyC親緣關(guān)系最遠(yuǎn)。qRT-PCR結(jié)果表明:MuCry的3個(gè)成員和MuPhyC1對不同光質(zhì)響應(yīng)存在差異,藍(lán)光都可促進(jìn)MuCry 3個(gè)成員 mRNA的轉(zhuǎn)錄,紅光促進(jìn)MuPhyC1 mRNA的轉(zhuǎn)錄,在黑暗、紅光、綠光和黃光下MuCry1的表達(dá)受到抑制;黃光和紅光對MuCry2a的轉(zhuǎn)錄水平影響不明顯,而黑暗、暖白光和綠光可抑制MuCry2a的表達(dá)水平;MuCry2b的轉(zhuǎn)錄水平受藍(lán)光和紅光的正調(diào)控,而受綠光的負(fù)調(diào)控;MuPhyC1在紅光處理下的表達(dá)量最高,藍(lán)光和暖白光下的表達(dá)量相近,白光、綠光、黃光和黑暗處理下的表達(dá)量差異不大,表明MuPhyC1積極響應(yīng)紅光刺激。

    關(guān)鍵詞 野生蕉;隱花色素;光敏色素;光質(zhì)處理;表達(dá)分析

    中圖分類號(hào) S668.1 文獻(xiàn)標(biāo)識(shí)碼 A

    Abstract The RT-PCR combined with RACE method was used to obtain Cryptochrome and Phytochrome gene family in the wild banana of Fuzhou, named MuCry1, MuCry2a, MuCry2b, MuPhyB and MuPhyC1. respectively. The cDNA sequences of MuCrys and MuPhys were 2 330, 2 859, 2 725, 3 272, 3 912 bp, encoding 698, 666, 669, 1 089 and 1 003 amino acids. Bioinformatics analysis showed that All MuCrys and MuPhys belong to the instability and hydrophilous proteins, and had transmembrane structures, and had no signal peptide, subcellularlocation predictions of them locating were in the nucleus. Phylogenetic anglicizing of Cry and Phy in plants indicated that not only Cry1 and Cry2, but also PhyB and PhyC belonged to the different branches, but they all had farthest genetic relationship with dicotyledons; MuCry1 have a close genetic relationship with MaCry1(Musa acuminata), EgCry1(Elaeis guineensis), and PdCry1(Phoenix dactylifera), MuCry2a and MuCry2b belongs to the same branch with a close genetic relationship with MuCry2, and TaCry2(Triticum aestivum), and OsCry2(Oryza sativaJaponica), MuPhyB have a close genetic relationship with MuPhyB(Musa acuminata), MuPhyC1 have a close genetic relationship with MuPhyC(Musa acuminata)and OmPhyC(Orchidantha maxillarioides). qRT-PCR results revealed that three members of Cryptochrome gene family and MuPhyC1 had different responses to different light quality, and blue light could promote the transcription of MuCrys mRNA, and red light could promote the transcription of MuPhyC1. The expression of MuCry1 was inhibited in dark, red, green and yellow light; The effect of light on the transcription level of MuCry2a was not obvious, while dark, warm white and green light could suppress the expression level of MuCry2a. The transcriptional level of MuCry2b was positively regulated by blue and red light, and was negatively controled by green light; MuPhyC1 had the highest expression in red light, The expression level of white red, green, yellow and dark light treatment was not significant, indicating that MuPhyC1 positively responded to red light stimulation.

    Key words Wild banana(Musa spp); Crypthchrome; Phytochrome; light quality treatment; gene expression

    doi 10.3969/j.issn.1000-2561.2017.11.016

    隱花色素(Cryptochrome,Cry)是重要的光受體之一,是一類通過吸收藍(lán)光和近紫外光[1],進(jìn)而影響植物生長發(fā)育形態(tài)、新陳代謝和其自身對向光性反應(yīng)的光敏受體[1-2],大多數(shù)植物的隱花色素都是黃素類蛋白,主要包括氨基酸PHR(photolyase-related)和羧基端CCT(Cryptochrome C Terminal)兩個(gè)明顯的功能區(qū)[1-3]。研究表明隱花色素與DNA光裂解酶結(jié)構(gòu)高度相似但功能各異,隱花色素?zé)o光裂解酶活性[1,4-5],其廣泛存在于細(xì)菌、植物、動(dòng)物和人體內(nèi)[6-7]。最近新發(fā)現(xiàn)的隱花色素家族成員Cry-DASH[8],僅在葉綠體和線粒體中檢測到,其缺少Cry蛋白典型的C末端延伸,但可能具有單鏈DNA光裂解酶活性[8-9]。植物隱花家族包括光裂解酶,Cry-DASH和植物隱花色素三類成員[1,7,9]。隱花色素作為藍(lán)光受體,其基因家族成員中Cry1和Cry2的C端包括一個(gè)可變區(qū)域,N端的相似性很高,包含碟呤和FAD兩種分子,Cry-DASH具有修復(fù)單鏈DNA的功能[8-9],隱花色素出現(xiàn)于真核和原核分化之前,其DNA結(jié)合和依賴氧化還原的功能與光解酶一致,是一種轉(zhuǎn)錄抑制劑[1,7,9,10]。光敏色素(Phytochrome,Phy)是由4個(gè)生色團(tuán)與脫輔機(jī)蛋白質(zhì)共價(jià)結(jié)合形成的一種易溶于水的調(diào)節(jié)生長發(fā)育的色素蛋白復(fù)合體,廣泛存在于植物界,在光合細(xì)菌中的含量較光合植物中高[10],光合植物中Phy的表達(dá)幾乎不受環(huán)境光照的影響。研究發(fā)現(xiàn)植物中光敏色素家族基因包括PhyA、PhyB、PhyC、PhyD和PhyE五個(gè)成員[4,11]。近年來植物中關(guān)于隱花色素的研究主要包括幼苗的去黃化[12-15]、植物開花調(diào)節(jié)[16-17]、植物生物鐘[18]和氣孔開放的調(diào)控[19]、引發(fā)晝夜節(jié)律性[20]和調(diào)控基因的表達(dá)模式[21]、并表現(xiàn)出趨光性彎曲[22]的特性等方面,其中Cry1蛋白主要介導(dǎo)生物鐘調(diào)節(jié)以及抑制幼苗下胚軸伸長[14],Cry2調(diào)控去黃化[12,14]等。近來還發(fā)現(xiàn)隱花色素參與植物對磁場的感應(yīng)[23]和細(xì)胞的程序性死亡[24],同時(shí)林辰濤等[25]證明了植物隱花色素的光誘導(dǎo)蛋白質(zhì)二聚化反應(yīng)為其原初光反應(yīng)的關(guān)鍵步驟,光敏色素對植物的生長發(fā)育的作用主要表現(xiàn)在調(diào)控植物花期[16-17]、介導(dǎo)植物光受體信號(hào)轉(zhuǎn)導(dǎo)[26]、促進(jìn)種子萌發(fā)[13-14,27]和改變植物抗逆性[28],通過調(diào)節(jié)植物自身的表達(dá)進(jìn)而調(diào)控作物的產(chǎn)量[21,29],其中PhyA介導(dǎo)遠(yuǎn)紅光信號(hào)下幼苗的光形態(tài)建成,是紅光下種子萌發(fā)和去黃化的主要光受體[26],同時(shí)PhyB也參與抑制胚軸的伸長和促進(jìn)種子萌發(fā)[14]以及紅光信號(hào)下幼苗的光形態(tài)建成,其缺失則影響種子葉柄延長[14]、發(fā)育以及早花現(xiàn)象[16,27],對于紅光感應(yīng)較弱的PhyC,與PhyB、PhyD、PhyE共同參與避蔭反應(yīng)和調(diào)節(jié)光下植物生長[13,19,28],PhyD基因啟動(dòng)子在植物不同組織的活性存在明顯差異[12,14]。此外隱花色素和光敏色素的磷酸化可改變其活性,其蛋白互作共同參與并調(diào)控紅藍(lán)光作用[10,13-14,29]。

    光是影響植物生長發(fā)育最重要的環(huán)境因子之一,隱花色素和光敏色素分別作為藍(lán)光和紅光受體,當(dāng)前研究已較為深入,但大多主要集中于擬南芥、水稻、小麥等模式植物幼苗去黃化、開花調(diào)節(jié)以及光形態(tài)建成等方面。目前關(guān)于隱花色素(MuCry)和光敏色素(MuPhy)在福州野生蕉(Musa spp.)[30]光形態(tài)建成方面的研究還未見報(bào)道,通過香蕉全基因組分析,香蕉基因組中隱花色素家族具有Cry1,Cry2a、Cry2b和Cry-DASH 4個(gè)成員,光敏色素有PhyA、PhyB、PhyC1和PhyC2 4個(gè)成員。香蕉作為世界4大水果之一,不僅在全世界種植廣泛,更是重要的經(jīng)濟(jì)作物和糧食作物,但由于栽培過程中遭受光照不足、病害、低溫等脅迫使得產(chǎn)量下降和品質(zhì)嚴(yán)重不佳。本研究以福州野生蕉為材料,克隆得到隱花色素和光敏色素基因家族不同成員的基因序列,進(jìn)行基因結(jié)構(gòu)以及蛋白的生物信息學(xué)分析,并借助qRT-PCR技術(shù)檢測各成員在不同光質(zhì)處理下的表達(dá)情況,以期為探討隱花色素響應(yīng)藍(lán)光和光敏色素響應(yīng)紅光的機(jī)制提供理論依據(jù),為隱花色素和光敏色素參與香蕉光質(zhì)應(yīng)答和光形態(tài)建成等方面的功能研究奠定理論基礎(chǔ),同時(shí)為解決生產(chǎn)中香蕉光照不足導(dǎo)致產(chǎn)量和品質(zhì)問題提供解決辦法的科學(xué)依據(jù)。

    1 材料與方法

    1.1 供試材料與光質(zhì)處理

    以福州野生蕉(Musa spp.)組培苗為供試材料提取總RNA;并選取大小和長勢一致的組培苗進(jìn)行光質(zhì)(紅、綠、黃、藍(lán)、白、暖白、黑暗)處理24 h,5個(gè)重復(fù)。取不同光質(zhì)(功率18 W)處理后的香蕉葉片于液氮中速凍,并保存于-80 ℃用于后續(xù)基因表達(dá)分析,以上材料均由福建農(nóng)林大學(xué)園藝植物生物工程研究所提供。

    1.2 方法

    1.2.1 RNA提取與cDNA合成 采用Column Plant RNAOUT 2.0試劑盒(TIANDZ,China)進(jìn)行香蕉總RNA的提取。經(jīng)1%瓊脂糖凝膠電泳檢測其完整性,后用紫外分光光度計(jì)檢測OD260/280在1.9~2.1之間且完整性良好的核酸備用。采用Thermo Scientific Rever-tAid First Strand cDNA Synthesis Kit(Fermentas,EU)將總RNA逆轉(zhuǎn)錄為帶AP(GGCCACGCGTC-GACTAGTACTTTTTTTTTTTTTTTTTT)接頭的cDNA,作為ORF和3′-UTR擴(kuò)增的模板。采用SMARTTM RACE cDNA Amplification kit(Takara,Japan)進(jìn)行cDNA逆轉(zhuǎn)錄,用于5′-UTR的擴(kuò)增。

    1.2.2 引物設(shè)計(jì)及PCR擴(kuò)增 保守區(qū)擴(kuò)增:根據(jù)GenBank中馬來西亞小果野蕉(Musa.acuminata,AA group)全基因組(http://banana-genome-hub.southgreen.fr/)結(jié)合GenBank中已經(jīng)登錄的多種植物相應(yīng)基因序列多重比對結(jié)果,選擇同源性比較高的基因片段進(jìn)行福州野生蕉MuCrys和MuPhys基因保守區(qū)的引物設(shè)計(jì)與PCR擴(kuò)增,選取與預(yù)期片段大小一致的條帶進(jìn)行回收,以pMD18-T為載體、DH5α為感受態(tài)對回收片段進(jìn)行TA克隆。菌液PCR后選取陽性克隆子測序,篩選與目的片段大小相符的結(jié)果,經(jīng)BLAST進(jìn)行同源比對確定保守區(qū)序列。

    3′-RACE和5′-RACE:采用RACE法,在得到保守序列的基礎(chǔ)上分別設(shè)計(jì)2條3′-RACE和5′-RACE特異引物,以福州野生蕉試管苗cDNA為模板,進(jìn)行巢式PCR反應(yīng)進(jìn)行3′未端序列和5′未端序列的擴(kuò)增。

    拼接驗(yàn)證:對克隆所得保守序列、3′未端序列和5′未端序列的測序結(jié)果進(jìn)行全長拼接,設(shè)計(jì)拼接驗(yàn)證引物,PCR反應(yīng)驗(yàn)證拼接結(jié)果。以上引物均由北京六合華大基因科技股份有限公司合成。所有引物的序列、擴(kuò)增的目的片段大小、退火溫度及擴(kuò)增用途見表1。

    PCR反應(yīng)體系和擴(kuò)增程序參照張銳[31]的方法,根據(jù)擴(kuò)增片段的不同,對PCR擴(kuò)增程序進(jìn)行相應(yīng)的調(diào)整。獲得目的片段后切膠回收,TA克隆后挑取陽性克隆子的菌液進(jìn)行PCR擴(kuò)增,將有目的條帶的菌液送至華大基因公司測序。

    1.2.3 基因序列分析 采用DNAMAN6.0進(jìn)行核苷酸序列比對和拼接。采用NCBI對福州野生蕉隱花色素和光敏色素家族基因及其推導(dǎo)的氨基酸序列進(jìn)行同源性比對;以ExPASy Protparam軟件預(yù)測編碼蛋白的理化性質(zhì);應(yīng)用SignalP 3.0 Server進(jìn)行編碼蛋白的信號(hào)肽預(yù)測;以PSORT亞細(xì)胞定位;TMpered軟件進(jìn)行蛋白質(zhì)跨膜結(jié)構(gòu)預(yù)測;應(yīng)用NetPhos2.0軟件進(jìn)行編碼蛋白的磷酸化位點(diǎn)預(yù)測,蛋白卷曲螺旋結(jié)構(gòu)則以EMBnet COILS軟件進(jìn)行預(yù)測;以NCBI blastp預(yù)測蛋白質(zhì)保守結(jié)構(gòu)域;經(jīng)PSIPRED在線軟件預(yù)測蛋白質(zhì)二級(jí)結(jié)構(gòu);蛋白質(zhì)三級(jí)結(jié)構(gòu)預(yù)測以SWISS-MODEL進(jìn)行;最后采用Mega5.05軟件的Neighbor-Joining(鄰位相連法,NJ法)構(gòu)建核苷酸序列的分子系統(tǒng)進(jìn)化樹(P-distance法),并用bootstrap法(重復(fù)1 000次)評(píng)估系統(tǒng)進(jìn)化樹。

    1.2.4 福州野生蕉隱花色素和光敏色素家族基因在不同光質(zhì)處理下的qRT-PCR分析 以Chen[32]等篩選的在非生物脅迫下表達(dá)較穩(wěn)定的CAC基因作為本研究定量分析的內(nèi)參基因,引物序列為CAC-qF(AACTCCTATGTTGCTCGCTTATG)和CAC-qR(GGC

    TACTACTTCGGTTCTTTCAC),采用SYBR ExScript試劑盒(Takara,Japan)和羅氏LightCycler480儀器(Thermo Fisher,USA),以經(jīng)不同光質(zhì)處理的福州野生蕉組培苗的cDNA為模板,根據(jù)熒光定量引物設(shè)計(jì)的原則在MuCrys和MuPhys各基因的特異位置設(shè)計(jì)上下游引物進(jìn)行qPCR擴(kuò)增。qRT-PCR反應(yīng)體系和擴(kuò)增程序參照馮[33]等的方法。待反應(yīng)結(jié)束后進(jìn)行擴(kuò)增曲線、溶解曲線(60~95 ℃)和凝膠電泳分析,檢測引物的特異性;每個(gè)反應(yīng)包括3個(gè)重復(fù),在樣品擴(kuò)增的同時(shí),將不同光質(zhì)處理的cDNA模板的混合樣進(jìn)行5倍梯度稀釋制作標(biāo)準(zhǔn)曲線。從擴(kuò)增曲線圖中得到Ct值,并根據(jù)標(biāo)準(zhǔn)曲線獲得不同處理下福州野生蕉隱花色素Crys和光敏色素Phys各基因的mRNA相對含量,通過內(nèi)參基因的校正最終得到目的基因的相對表達(dá)量,進(jìn)行福州野生蕉隱花色素Crys和光敏色素Phys各基因在不同光質(zhì)處理下的表達(dá)情況分析。數(shù)據(jù)分析采用Ex-cel和geNORM(version3.5)[34]軟件進(jìn)行。

    2 結(jié)果與分析

    2.1 福州野生蕉MuCrys和MuPhys基因家族成員cDNA全長序列克隆及分析

    本研究通過同源克隆的方法分別得到了MuCrys和MuPhys基因家族的保守區(qū),并采用RACE法擴(kuò)增得到MuCrys和MuPhys的3′末端序列和5′末端序列,通過DNAMAN6.0拼接驗(yàn)證,得到MuCry1的cDNA序列2 330 bp,MuCry2a的cDNA全長2 859 bp,MuCry2b的cDNA全長2 725 bp,MuPhyB的cDNA序列3 272 bp,MuPhyC1的cDNA全長3 912 bp(具體見表2)。MuCry1和MuCry2a的終止密碼子是TAG,MuCry2b和MuPhyC1的終止密碼子是TAA,MuPhyB的終止密碼子是TGA。將以上序列的核苷酸和推導(dǎo)的氨基酸序列分別在NCBI上進(jìn)行Blast分析,結(jié)果顯示MuCrys和MuPhys與數(shù)據(jù)庫中已知的毛果楊(Populus tremula)、海棗(Phoneix dactylifera)、油棕(Elaeis guineensis)等的Crys、Phys基因序列高度同源。因此推斷已經(jīng)成功克隆得到福州野生蕉隱花色素Crys基因家族和光敏色素Phys基因家族的序列,分別命名為MuCry1(登錄號(hào)KX236155)、MuCry2a(登錄號(hào)KX236156)、MuCry2b(登錄號(hào)KX236157)、MuPhyB(登錄號(hào)KX247367)、MuPhyC1(登錄號(hào)KX247368)。

    2.2 福州野生蕉MuCrys和MuPhys的生物信息學(xué)分析

    利用ExPASy Protparam預(yù)測MuCry1、MuCry2a、MuCry2b、MuPhyB和MuPhyC1的理化性質(zhì),具體見表3。其信號(hào)肽、跨膜結(jié)構(gòu)及其磷酸化位點(diǎn)預(yù)測等見表2,其中磷酸化位點(diǎn)MuCry1(Ser:27,Thr:6,Tyr:7)、MuCry2a(Ser:30,Thr:3,Tyr:7)、MuCry2b(Ser:33,Thr:3,Tyr:6)、MuPhyB (Ser:40,Thr:6,Tyr:12)和MuPhyC1(Ser:32,Thr:8,Tyr:6)。以PSORT預(yù)測可知MuCry1定位于微體,MuCry2a定位于細(xì)胞核和微體,MuCry2b定位于微體和質(zhì)膜,MuPhyB定位于線粒體類囊體膜、線粒體基質(zhì)和質(zhì)膜,MuPhyC1定位于線粒體基質(zhì)和質(zhì)膜。經(jīng)NCBI blastp預(yù)測,MuCrys都含有crypt_chrom_pln保守結(jié)構(gòu)域(20~484 aa),其中MuCry1的N端含有與DNA光裂合酶相關(guān)(photolyase-related,PHR)的保守結(jié)構(gòu)域(19~181 aa),中間為FAD綁定功能域(221~498 aa),C端為Cryptochrome-C超家族(524~637 aa);MuCry2a和MuCry2b的N端分別含有與DNA光裂合酶相關(guān)(photolyase-related,PHR)的保守結(jié)構(gòu)域(7~161 aa)和(7~162 aa),中間為FAD綁定功能域(212~489 aa),但C端無Cryptochrome-C超家族(圖1)。光敏色素包括有兩個(gè)中心區(qū)域,即N末端光感受區(qū)和C末端光調(diào)節(jié)區(qū),其中MuPhyB的N末端光感受區(qū)包括一個(gè)PAS光感應(yīng)區(qū)域(122~240 aa),GAF(cGMP phosphodiesterase-adenylyl cyclas-FhlA,GAF)綁定功能域(273~411 aa)和光敏色素結(jié)構(gòu)域PHY(418~588 aa),C末端光調(diào)節(jié)區(qū)包括PAS光感應(yīng)區(qū)(518~632 aa),HATPase區(qū)域(868~981 aa),HisKA區(qū)域(761-819 aa),MuPhyC1的N末端光感受區(qū)包括一個(gè)PAS光感應(yīng)區(qū)域(57~141 aa),GAF綁定功能域(174~308 aa)和光敏色素結(jié)構(gòu)域PHY(315~487 aa),C末端光調(diào)節(jié)區(qū)包括PAS光感應(yīng)區(qū)(624~681 aa),PAS-8區(qū)域(725~833 aa),HATPase區(qū)域(967~1 079 aa),HisKA區(qū)域(855~910 aa),見圖2。經(jīng)PSIPRED在線軟件預(yù)測蛋白質(zhì)二級(jí)結(jié)構(gòu),MuCry1、MuCry2a和MuCry2b主要由無規(guī)則卷曲和α-螺旋組成,MuCry1中α-螺旋31.8%、無規(guī)則卷曲61.9%、β-折疊6.3%;MuCry2a中α-螺旋34.8%、無規(guī)則卷曲59.5%、β-折疊5.7%;MuCry2b中α-螺旋35.7%、無規(guī)則卷曲57.8%、β-折疊6.5%;MuPhyB主要以α螺旋結(jié)構(gòu)為主,所占比例為45.91%;無規(guī)則卷曲為40.22%;β折疊13.87%;而MuPhyC11蛋白二級(jí)結(jié)構(gòu)中無規(guī)則卷曲占43.37%;α螺旋41.87%,β折疊14.76%。采用SWISS-MODEL軟件進(jìn)行三級(jí)結(jié)構(gòu)預(yù)測,可見MuCrys蛋白的三級(jí)結(jié)構(gòu)比較相近,N端區(qū)域主要由多個(gè)α-螺旋和β-折疊共同構(gòu)成復(fù)雜的結(jié)構(gòu),而C端區(qū)域主要由無規(guī)則卷曲組成,與二級(jí)結(jié)構(gòu)預(yù)測的結(jié)果相符(圖3),MuPhyB和MuPhyC1的三級(jí)結(jié)構(gòu)差異較大。最后采用Mega5.05 軟件的Neighbor-Joining(鄰位相連法,NJ法)構(gòu)建核苷酸序列的分子系統(tǒng)進(jìn)化樹(P-distance法),并用bootstrap法(重復(fù)1 000次)評(píng)估系統(tǒng)進(jìn)化樹(圖4、圖5)。

    2.3 氨基酸序列比對及系統(tǒng)進(jìn)化分析

    本研究通過同源克隆的方法分別得到了福州野生蕉隱花色素基因家族MuCry1基因的保守區(qū)和5′末端、MuCry2a和MuCry2b的cDNA全長,光敏色素基因家族MuPhyB的保守序列以及MuPhyC1的cDNA全長。并通過NCBI進(jìn)行blastn比對發(fā)現(xiàn),MuCry1與油棕(Elaeis guineensis,XP_010939647.1)Cry1基因的同源性高達(dá)81%,MuCry2a與海棗(Phoenix dactylifera,XP_008790851.1)Cry2a的同源性達(dá)75%,MuCry2b與菠蘿(Ananas comosus,OAY67638.1)Cry2b基因的同源性高達(dá)74%,MuPhyB與海棗(Phoenix dactylifera,XP_008790851.1)PhyB的同源性高達(dá)77%,MuPhyC1與毛果楊(Populus tremula,XP_002318913.1)PhyC的同源性達(dá)75%。為了進(jìn)一步研究MuCrys和MuPhys各基因的保守性,經(jīng)NCBI blastp預(yù)測,MuCrys都含有保守結(jié)構(gòu)域crypt_chrom_pln,MuCrys的N端含有與DNA光裂合酶相關(guān)(photolyase-related,PHR)的保守結(jié)構(gòu)域,中間為FAD綁定功能域,且MuCry1C端具Cryptochrome-C超家族。NCBI blastn預(yù)測表明MuPhyB和MuPhyC1各含有一段可變剪接體,可能會(huì)造成MuPhyB和MuPhyC1在香蕉表達(dá)中的功能缺失。為了研究植物MuCry和MuPhy的進(jìn)化關(guān)系,根據(jù)NCBI數(shù)據(jù)庫中油棕(Elaeis guineensis)、海棗(Phoenix dactylifera)、玉米(Zea mays)、小果野蕉(Musa acuminata)、菠蘿(Ananas comosus)、番茄(Solanum lycopersicum)和煙草(Nicotiana sylvestris)等的40多條Cry和Phy的氨基酸序列,利用Mega5.05軟件的鄰近相鄰法(NJ法)構(gòu)建植物Cry和Phy的系統(tǒng)進(jìn)化樹(圖4、5),結(jié)果表明:植物Cry1和Cry2根據(jù)種屬聚為兩大類,其中MuCry1與小果野蕉、海棗和油棕Cry1的親緣關(guān)系最近,與雙子葉植物大豆Cry1的親緣關(guān)系最遠(yuǎn);MuCry2a和MuCry2b與小果野蕉、水稻和小麥Cry2的親緣關(guān)系最近,與雙子葉植物碧桃Cry2的親緣關(guān)系最遠(yuǎn);MuPhyB與小果野蕉PhyB的親緣關(guān)系最近,與雙子葉毛果楊PhyB的親緣關(guān)系最遠(yuǎn);MuPhyC1與馬來蘭花蕉PhyC的親緣關(guān)系最近,與雙子葉甜橙PhyC親緣關(guān)系最遠(yuǎn)。

    2.4 福州野生蕉隱花色素MuCrys和光敏色素MuPhyC1在不同光質(zhì)條件下的表達(dá)分析

    采用qRT-PCR法研究福州野生蕉隱花色素MuCrys和光敏色素MuPhyC1基因家族共4個(gè)成員在不同光質(zhì)處理(黑暗、紅光、綠光、黃光、藍(lán)光、白光和暖白光)條件下的表達(dá)模式(圖6),結(jié)果表明不同光質(zhì)處理下,MuCry1對不同光質(zhì)響應(yīng)存在差異,其中藍(lán)光和暖白光下的表達(dá)量高于對照(白光),尤其是藍(lán)光下的相對表達(dá)量最高,說明藍(lán)光刺激可促進(jìn)MuCry1 mRNA的轉(zhuǎn)錄,而在黑暗,紅光,綠光和黃光下表達(dá)量低于對照,說明這些光質(zhì)處理抑制了MuCry1的表達(dá)。MuCry2a在藍(lán)光處理下的表達(dá)量高于對照,黃光和紅光處理與對照的相對表達(dá)水平相近,而黑暗、暖白光和綠光處理下表達(dá)量低于對照,說明藍(lán)光可提高M(jìn)uCry2a的表達(dá)水平。MuCry2b也差異性的響應(yīng)不同光質(zhì)的處理,在藍(lán)光下的表達(dá)量最高,依次分別為紅光,黃光,白光(CK),暖白,綠光和黑暗,說明MuCry2b的轉(zhuǎn)錄水平受藍(lán)光和紅光的正調(diào)控,而受綠光的負(fù)調(diào)控。MuPhyC1對于不同的光處理也顯示出差異性表達(dá),在紅光下的表達(dá)量最高,說明MuPhyC1 mRNA的轉(zhuǎn)錄受紅光的影響,其在藍(lán)光和暖白光下的表達(dá)量相近,白光、綠光、黃光和黑暗處理下的表達(dá)量差異不大,說明MuPhyC1積極響應(yīng)紅光刺激。實(shí)時(shí)熒光定量PCR的結(jié)果說明了福州野生蕉隱花色素MuCry基因家族3個(gè)成員的光響應(yīng)表達(dá)模式相似,它們都能正向響應(yīng)藍(lán)光的刺激,綠光和黑暗處理則對3個(gè)成員的表達(dá)均起抑制作用。而光敏色素家族基因MuPhyC1作為紅光的受體基因之一,其轉(zhuǎn)錄水平受紅光的正調(diào)控,受綠光的負(fù)調(diào)控。結(jié)果表明藍(lán)光和綠光能調(diào)控隱花色素基因的表達(dá),紅光和綠光調(diào)控光敏色素的表達(dá),兩者通過在光形態(tài)建成中起作用乃至相互作用,進(jìn)而影響植物的生長發(fā)育過程,參與光機(jī)制的建成。

    3 討論

    隱花色素和光敏色素家族基因在不同物種中的成員數(shù)量不同,如在擬南芥[8]中存在3個(gè)Cry成員,水稻、番茄和燕麥[3,6]中也存在至少3個(gè)Cry成員,在蘋果中也存在2個(gè)Cry成員,龍眼中存在3個(gè)Cry成員,蕨類和苔蘚[35]中亦存在至少5個(gè)Cry基因;光敏色素在蘋果、蘿卜[15]等雙子葉植物中有至少五個(gè)成員,而在香蕉、小麥、玉米[15]等單子葉植物中,目前只發(fā)現(xiàn)了三個(gè)成員。經(jīng)過香蕉全基因組分析檢索獲得MuCry1、MuCry2a、MuCry2b和MuCry-DASH等4個(gè)具有完整CDS的Cryptochrome家族基因,分別位于6號(hào)、5號(hào)(MuCry2a、MuCry2b)和1號(hào)染色體上;同時(shí)獲得MuPhyA、MuPhyB、MuPhyC1和MuPhyC2等4個(gè)具有完整CDS的Phytochrome家族基因,分別位于1號(hào)、3號(hào)、6號(hào)和4號(hào)染色體上。生物信息學(xué)分析表明MuCrys和MuPhys均屬親水蛋白,具跨膜結(jié)構(gòu),不具信號(hào)肽,均定位于細(xì)胞核,MuCrys都含有crypt_chrom_pln保守結(jié)構(gòu)域,但MuCry2a和MuCry2b無Cryptochrome-C超家族。MuCrys氨基酸末端區(qū)域PHR與光解酶的序列同源,但光解酶結(jié)合在DNA損傷部位和催化DNA修復(fù)的氨基酸殘基在隱花色素中不存在,因此除Cry-DASH外其余隱花色素雖為感光受體,但不具有光解酶的活性[2]。MuCrys和MuPhys的生物信息學(xué)分析結(jié)果與前人研究相符,故初步推斷MuCrys和MuPhys為福州野生蕉隱花色素和光敏色素,其可能參與并調(diào)控福州野生蕉光形態(tài)建成、幼苗去黃化以及蛋白互作等功能。

    實(shí)時(shí)熒光定量PCR分析MuCrys和MuPhyC1在不同光質(zhì)下的表達(dá)調(diào)控模式,結(jié)果表明MuCrys和MuPhyC1對不同光質(zhì)響應(yīng)存在差異,藍(lán)光刺激可促進(jìn)MuCrys mRNA的轉(zhuǎn)錄,MuCrys都表現(xiàn)出積極響應(yīng)藍(lán)光的特性,同時(shí)MuPhyC1正向調(diào)控紅光刺激,其mRNA的轉(zhuǎn)錄受紅光控制,作為編碼藍(lán)光受體的蛋白,MuCry2b在紅光條件下的相對表達(dá)量也較高,說明MuCry2b和MuPhyB蛋白互作進(jìn)而調(diào)控光形態(tài)建成。擬南芥突變體研究表明,Cry1介導(dǎo)抑制擬南芥下胚軸伸長[2];植物去黃化反應(yīng)的研究發(fā)現(xiàn),Cry1與PhyA和PhyB互作可以抑制種子下胚軸伸長和類黃酮生物合成[10,13,29],促使子葉開放和改變基因表達(dá)[20],除隱花色素外,光敏色素對藍(lán)光條件下的去黃化也起調(diào)節(jié)作用,擬南芥和番茄的研究表明藍(lán)光下PhyA和PhyB共同調(diào)節(jié)去黃化[4,8,19],香蕉生長過程中易出現(xiàn)幼苗黃化現(xiàn)象,通過調(diào)控隱花色素和光敏色素可控制其黃化。同時(shí)藍(lán)光的磷酸化作用會(huì)因磁場變化而影響隱花色素的活性[9,27],Ahmad等[7]認(rèn)為Cry2需要在紅光下被Phy磷酸化才能完全激活,MuCrys和MuPhys均含有極其豐富的蛋白磷酸化位點(diǎn),除此之外MuCrys含有其它潛在功能位點(diǎn),如N-糖基化位點(diǎn),蛋白酶C磷酸化位點(diǎn)。因此隱花色素還可能通過信號(hào)轉(zhuǎn)導(dǎo)、蛋白磷酸化或其他潛在功能位點(diǎn)等方式調(diào)控目標(biāo)蛋白發(fā)揮其功能[38]。不同光質(zhì)通過觸發(fā)光受體刺激其感知光信號(hào),進(jìn)而影響植物的光合特性、生長發(fā)育以及抗逆和衰老等[26,39]。香蕉生長過程中光照不足導(dǎo)致光合速率減緩,香蕉抗逆性降低,容易遭受一系列脅迫影響其生長發(fā)育,同時(shí)光照不足導(dǎo)致香蕉果實(shí)著色差,風(fēng)味不佳,嚴(yán)重影響其產(chǎn)量和品質(zhì),因而研究隱花色素和光敏色素在香蕉光質(zhì)應(yīng)答以及光形態(tài)建成方面等的作用和機(jī)理顯得至關(guān)重要,其研究可以為香蕉產(chǎn)業(yè)增產(chǎn)增收增質(zhì)奠定理論基礎(chǔ)。

    參考文獻(xiàn)

    [1] Lin C T. Blue light receptors and signal transduction[J]. Plant Cell, 2002, 14(1): S207-S225.

    [2] ZHU C L(朱春利), ZHANG G R(張桂榮), CAI A J(蔡愛軍), et al. Advances in the Structure and Function of Cryptochrome in Plant[J]. Genomics and Applied Biology(基因組學(xué)與應(yīng)用生物學(xué)), 2009, 28(01): 174-178.

    [3] Li Q H, Yang H Q. Cryptochrome Signaling in Plants[J]. Photochemistry and Photobiology, 2007, 83(1): 94-101.

    [4] Todo T. Functional diversity of the DNA photolyase/blue light receptor family[J]. Mutation Research, 1999, 434(2): 89-97.

    [5] Kleine T, Lockhart P, Batschauer A. An Arabidopsis protein closely related to Synechocystis cryptochrome is targeted to organelles[J]. Plant, 2003, 35(1): 93-103.

    [6] Cashmore A R, Jarillo J A, Wu Y J, et al. Cryptochromes: blue light receptors for plants and animals[J]. Science, 1999, 284(5 415): 760-765.

    [7] Ahmad M, Cashmore A R. HY4 gene of Arabidopsis thaliana encodes a protein with characteristics of a blue-light photoreceptor[J]. Nature, 1993, 366: 162-166.

    [8] Castrillo M, Bernhardt A, Batschauer A, et al. Biochemical Characterization of the DASH-Type Cryptochrome CryD from Fusarium fujikuroi[J]. Photochemistry and Photobiology, 2015, 91(6): 1 356-1 367.

    [9] YAN H(閆 海), ZHOU B(周 波), LI Y H(李玉花), et al. Plant Cryptochrome and Its Light Signal Transduction[J]. Plant Physiology Communications(植物生理學(xué)通訊), 2005, 41(6): 841-846.

    [10] KhudairA K A, Rboleda O P. Phytochrome-mediated carotenoid biosynthesis and its influence by plant hormones[J]. Physiology Plant, 1971, 10(1 111): 8-22.

    [11] Alba R, K elmenson P M, Cordonnier-Pratt M M, et al. The phytochrome gene family in tomato and rapid different evolution of this family in angiosperms[J]. Molecular Biology and Evolution, 2000, 17(3): 362-373

    [12] Lin C, Yang H, Cashmore A R, et al. Enhancement of blue-light sensitivity of Arabidopsis seedlings by a blue light receptor cryptochrome2[J]. Science, 1998, 95(5): 2 686-2 690.

    [13] Quail P H. Phytochromes: photosensory perception and signal transduction[J]. Science, 1995, 268(5 211): 675-680.

    [14] Zeng Jianxin, Wang Qiming, Lin Jianzhong, et al. Arabidopsis cryptochrome-1 restrains lateral roots growth by inhibiting auxin transport[J]. Journal of Plant Physiology, 2010, 167(8): 670.

    [15] Fankhauser C, Chory J. Light control of seeding development[J]. Plant Physiology and Plant Molecular Biology, 1997, 47(7): 215-243.

    [16] Hirose F, Shinomura T, Tanabata T, et al. Involvement of rice cryptochromes in de-etiolation responses and flowering[J]. Plant and Cell Physiology, 2006, 47(7): 915-925.

    [17] Guo H W, Yang H Y, Mockler T C, et al. Regulation of flowering time by Arabidopsis photoreceptors[J]. Science, 1998, 279(5 355): 1 360-1 363.

    [18] Somers D E, Devlin P F, Kay S A. Phytochromes and cryptochromes in the entrainment of the Arabidopsis circadian clock[J], Science, 1998, 282: 1 488-1 490

    [19] Kang C Y, Lian H L, Wang F F, et al. Cryptochromes, Phytochromes and COP1 regulate light-controlled stomatal development in Arabidopsis[J]. Plant Cell, 2009, 21(9): 2 624-2 641

    [20] Xu Chun xiao, Yan lu. Blue light-dependent phosphorylations of cryptochromes are affected by magnetic fields in Arabidopsis[J]. Advances in Space Research, 2014, 53(7): 1 118-1 124.

    [21] Imaizumi T, Kadota M, Wada M. Cryptochrome light signals control development to suppress auxin sensitivity in the moss Physcomitrella patens[J]. Plant Cell, 2002, 14: 373-386.

    [22] Zhao X(趙 翔), Zhao Q P(趙青平), Yang X(楊 煦), et al. Specificity and Crosstalk of Phototropin with Cryptochrome and Phytochrome in Regulating Hypocotyl Phototropism[J]. Chinese Bulletin of Botany(植物學(xué)報(bào)), 2015, 50(1): 122-132.

    [23] Johnsen S, Mattern E, Ritz T. Light-dependent magnetoreception: quantum catches and opponency mechanisms of possible photosensitive molecules[J]. Experimental Biology, 2007, 210: 3 171-3 178.

    [24] Danon A., Coll N. S and Apel K. Cryptochrome-1-dependent execution of programmed cell death induced by singlet oxygen in Arabidopsis thaliana[J]. Proceeding of the National Academy of Sciences, USA, 2006, 103(45): 17 036-17 041.

    [25] Wang Q, Zuo Z C, Lin C T, et al. Photoactivation and inactivation of Arabidopsis cryptochrome2[J]. Science, 354(6310): 343-347.

    [26] Briggs W R, Olney M A. Photoreceptors in plant photomorpho genesis to date five phytochromes, two cryptochromes, one phototropin, and one superchrome[J]. Plant Physiology, 2001, 1(125): 85-88.

    [27] Casal J J, Miguel L D, Burgin M J. Antagonistic action of low-fluence and high-irradiance modes of response of phytochrome on germination and β-mannanase activity in Datura ferox seed[J]. Journal of Experimental Botany, 2000, 51(347): 1 127-1 133.

    [28] Sharrock R A, Lissemore J L, Quail P H. Nucleotide and amino acid sequence of a Cucurbita phytochrome cDNA clone, identification of conserved features by comparision with Aveda phytochrome[J]. Gene, 1986, 47(2-3): 287-95.

    [29] Yang H Q, Wu Y J, Tang R H, et al. The C termini of Arabidopsis cryptochromes mediate a constitutive light response[J]. Cell, 2000, 103(5): 815-827.

    [30] Lai Z X(賴鐘雄), Chen Y(陳 源), Lin Y L(林玉玲), et al. Discovery and taxonomy of wild banana (Musa spp., ‘AA Group)in Fuzhou[J]. Subtropical Agriculture Research(亞熱帶農(nóng)業(yè)研究), 2007, 3(1): 1-5(in China).

    [31] 張 銳.福建野生蕉資源RAPD分析、離體培養(yǎng)與SOD基因克隆[D].福州: 福建農(nóng)林大學(xué), 2011.

    [32] Chen L, Zhong H Y, Kuang J F, et al. Validation of reference genes for quantitative RT-qPCR studies of gene expression in banana fruit under different experimental conditions[J]. BMC Research Notes, 2011, 234(2): 377-390.

    [33] Feng X(馮 新), Xu L(徐 蕾), Lai G T(賴恭梯), et al. Clone and Expression Analysis under Cold Stress of FeSOD Family Genes from the Wild Banana in Fuzhou[J]. Acta Botanica Boreali-Occidentalia Sinica(西北植物學(xué)報(bào)), 2015, 35(5): 853-864.

    [34] Vandesompele J, D E Preter K, Pattyn F, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes[J]. Genome Biology, 2002, 3(7): 31-34.

    [35] Small D G, Min B, Lefebvre P A. Characterization of a Chlamydomonas reingardtii gene encoding a protein of the DNA photolyase/blue light photoreceptor family[J]. Plant Molecular Biology, 1995, 28(3): 443-454.

    [36] Jiao Y, Lau S O, Deng X W. Light-regulated transcriptional networks in higher plants[J]. Nature Reviews Genetics, 2007, 8: 217-230.

    [37] Toth R, Kevei E, Hall A, et al. Circadian clock-regulated expression of phytochrome and cryptochrome genes in Arabidopsis[J]. Plant Physiology, 2001, 127(4): 1 607-1 616.

    [38] Wang X Q, Huang W D. Effects of Weak Light on the Ultmstructural Variatiom of Phloem Tissues in Source Leaves of Three-Year-Old Nectarine Trees[J]. Journal of Integrative Plant Biology(植物學(xué)報(bào)), 2003, 45(6): 688-697.

    [39] Franklin K A, Davis S J. Mutant analyses define multiple roles for phytochromes C in Arabidopsis thaliana photo-morphogenesis[J]. Plant Cell, 2003, 21(15): 981-1 989.

    猜你喜歡
    表達(dá)分析
    巴西橡膠樹CCaMK基因家族成員的鑒定和表達(dá)分析
    太子參玉米黃質(zhì)環(huán)氧化酶基因的克隆與表達(dá)分析
    木薯PHOR1基因的序列分析及其乙烯和茉莉酸甲酯誘導(dǎo)表達(dá)特性
    雷公藤貝殼杉烯酸氧化酶基因的全長cDNA克隆與表達(dá)分析
    紅花生育酚環(huán)化酶基因的克隆及表達(dá)分析
    三色堇DFR基因的克隆及表達(dá)分析
    香蕉抗壞血酸過氧化物酶基因的克隆及表達(dá)分析
    膠孢炭疽菌漆酶基因Lac2的序列特征與表達(dá)分析
    玉米紋枯病病菌y—谷氨酰轉(zhuǎn)肽酶基因克隆與表達(dá)分析
    信號(hào)分子與葉銹菌誘導(dǎo)下小麥病程相關(guān)蛋白1基因的表達(dá)分析
    一边摸一边抽搐一进一小说| 黄色 视频免费看| 国产三级在线视频| 国产午夜精品久久久久久| 国产精品久久视频播放| 91麻豆av在线| 精品卡一卡二卡四卡免费| 国产午夜福利久久久久久| 美女大奶头视频| 老汉色av国产亚洲站长工具| 国产亚洲av嫩草精品影院| 亚洲一区高清亚洲精品| 精品久久久久久久人妻蜜臀av| 欧美三级亚洲精品| 18禁国产床啪视频网站| 一区二区三区精品91| 美女扒开内裤让男人捅视频| 婷婷精品国产亚洲av| 日韩精品青青久久久久久| 亚洲精品美女久久久久99蜜臀| 一级片免费观看大全| 国内精品久久久久久久电影| 男女那种视频在线观看| 午夜两性在线视频| 亚洲男人天堂网一区| 色精品久久人妻99蜜桃| 久久人妻av系列| 久热爱精品视频在线9| 美女午夜性视频免费| 香蕉丝袜av| 亚洲欧美激情综合另类| 国产精品亚洲美女久久久| 韩国精品一区二区三区| 女生性感内裤真人,穿戴方法视频| 中国美女看黄片| av片东京热男人的天堂| 日日摸夜夜添夜夜添小说| 黑人巨大精品欧美一区二区mp4| 午夜福利一区二区在线看| 亚洲aⅴ乱码一区二区在线播放 | 亚洲男人的天堂狠狠| 成人永久免费在线观看视频| 99re在线观看精品视频| 欧美在线一区亚洲| 国产片内射在线| 精品福利观看| 18禁黄网站禁片免费观看直播| 免费观看人在逋| 久久久久久亚洲精品国产蜜桃av| 欧美不卡视频在线免费观看 | 黄色视频不卡| 国产麻豆成人av免费视频| 91大片在线观看| 在线十欧美十亚洲十日本专区| 婷婷精品国产亚洲av| 国产精品久久久av美女十八| 国产精品自产拍在线观看55亚洲| 女警被强在线播放| 黑人欧美特级aaaaaa片| 亚洲专区国产一区二区| 日韩精品青青久久久久久| 成人亚洲精品一区在线观看| 国产熟女午夜一区二区三区| 国产在线精品亚洲第一网站| 亚洲人成网站在线播放欧美日韩| 国产av一区二区精品久久| av在线天堂中文字幕| 欧美色视频一区免费| 国产精品久久电影中文字幕| 国产伦在线观看视频一区| 亚洲精品国产区一区二| 免费av毛片视频| 日韩欧美免费精品| 成在线人永久免费视频| 一级作爱视频免费观看| 午夜免费鲁丝| 中文亚洲av片在线观看爽| www日本黄色视频网| 亚洲狠狠婷婷综合久久图片| 可以在线观看的亚洲视频| 精品一区二区三区av网在线观看| 精品国产亚洲在线| 欧美午夜高清在线| 免费搜索国产男女视频| 国内久久婷婷六月综合欲色啪| 欧美乱码精品一区二区三区| 成人av一区二区三区在线看| 久久中文看片网| 日日干狠狠操夜夜爽| 国产成人一区二区三区免费视频网站| 真人做人爱边吃奶动态| 亚洲国产精品sss在线观看| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久亚洲精品国产蜜桃av| 韩国av一区二区三区四区| 亚洲成人久久性| 国产av一区二区精品久久| 满18在线观看网站| 国产v大片淫在线免费观看| 亚洲精品在线美女| 丰满的人妻完整版| 亚洲国产欧美一区二区综合| 国内精品久久久久精免费| 亚洲欧美精品综合一区二区三区| 男人舔女人下体高潮全视频| 国产私拍福利视频在线观看| 午夜福利视频1000在线观看| 欧美一级a爱片免费观看看 | 巨乳人妻的诱惑在线观看| 免费av毛片视频| 精品免费久久久久久久清纯| 特大巨黑吊av在线直播 | 久久久久久国产a免费观看| 久久性视频一级片| 淫妇啪啪啪对白视频| 女性被躁到高潮视频| 男人舔奶头视频| www.999成人在线观看| 亚洲五月天丁香| 91麻豆av在线| 国产免费男女视频| 国产av在哪里看| 亚洲,欧美精品.| 日韩高清综合在线| www.精华液| 亚洲一区二区三区不卡视频| 美女大奶头视频| 免费看美女性在线毛片视频| 色综合欧美亚洲国产小说| 好看av亚洲va欧美ⅴa在| 男人舔奶头视频| 国产欧美日韩精品亚洲av| 久久性视频一级片| 亚洲精品av麻豆狂野| 妹子高潮喷水视频| www.999成人在线观看| 男人的好看免费观看在线视频 | 日韩大尺度精品在线看网址| 妹子高潮喷水视频| 啦啦啦 在线观看视频| 91国产中文字幕| 国内久久婷婷六月综合欲色啪| 欧美精品啪啪一区二区三区| 久久人人精品亚洲av| videosex国产| bbb黄色大片| 亚洲 欧美一区二区三区| 久久久久久免费高清国产稀缺| 他把我摸到了高潮在线观看| 午夜福利视频1000在线观看| 免费在线观看日本一区| 亚洲国产毛片av蜜桃av| 亚洲一码二码三码区别大吗| 啦啦啦观看免费观看视频高清| 999久久久精品免费观看国产| 国产一卡二卡三卡精品| 在线观看免费午夜福利视频| 欧美日韩瑟瑟在线播放| 欧美激情 高清一区二区三区| 免费人成视频x8x8入口观看| 淫妇啪啪啪对白视频| 麻豆久久精品国产亚洲av| 男女床上黄色一级片免费看| 国产精华一区二区三区| 国产免费男女视频| 久久99热这里只有精品18| 亚洲七黄色美女视频| 日韩欧美国产在线观看| 老司机靠b影院| 亚洲av成人av| 好男人在线观看高清免费视频 | 日韩精品免费视频一区二区三区| 日韩大码丰满熟妇| 两个人看的免费小视频| 亚洲国产欧美日韩在线播放| 成在线人永久免费视频| 午夜成年电影在线免费观看| 黄色成人免费大全| 免费高清在线观看日韩| 日韩欧美在线二视频| 身体一侧抽搐| 亚洲狠狠婷婷综合久久图片| 两个人看的免费小视频| 91字幕亚洲| 视频区欧美日本亚洲| aaaaa片日本免费| 免费一级毛片在线播放高清视频| 国产三级在线视频| 国产成人精品久久二区二区91| 国产激情久久老熟女| 悠悠久久av| 女生性感内裤真人,穿戴方法视频| 国产国语露脸激情在线看| 亚洲在线自拍视频| 国产又爽黄色视频| 俺也久久电影网| 国产精品一区二区免费欧美| 18禁黄网站禁片免费观看直播| 日本 欧美在线| 成人18禁高潮啪啪吃奶动态图| 亚洲人成网站在线播放欧美日韩| 狠狠狠狠99中文字幕| 成人亚洲精品av一区二区| 久久国产精品影院| 亚洲熟妇熟女久久| 香蕉国产在线看| 看黄色毛片网站| 免费看a级黄色片| 国产伦人伦偷精品视频| 伦理电影免费视频| 一级a爱片免费观看的视频| 色播亚洲综合网| 色尼玛亚洲综合影院| 精品少妇一区二区三区视频日本电影| 岛国视频午夜一区免费看| 好男人在线观看高清免费视频 | 极品教师在线免费播放| 亚洲成国产人片在线观看| 757午夜福利合集在线观看| 男女下面进入的视频免费午夜 | 757午夜福利合集在线观看| 男女下面进入的视频免费午夜 | 精品第一国产精品| 久久国产精品影院| videosex国产| 特大巨黑吊av在线直播 | 一区二区三区精品91| 国内精品久久久久精免费| 天天躁夜夜躁狠狠躁躁| 99国产综合亚洲精品| 精品欧美一区二区三区在线| 国产成+人综合+亚洲专区| 亚洲国产精品sss在线观看| 夜夜爽天天搞| 法律面前人人平等表现在哪些方面| 免费在线观看影片大全网站| 久久国产乱子伦精品免费另类| 天天添夜夜摸| 亚洲色图 男人天堂 中文字幕| 中文在线观看免费www的网站 | 伊人久久大香线蕉亚洲五| 国内毛片毛片毛片毛片毛片| 亚洲成国产人片在线观看| 亚洲午夜精品一区,二区,三区| 国产精品久久久人人做人人爽| 嫁个100分男人电影在线观看| 淫妇啪啪啪对白视频| 亚洲成人久久性| a级毛片在线看网站| 国产精品一区二区免费欧美| 日韩成人在线观看一区二区三区| 18美女黄网站色大片免费观看| 精品一区二区三区av网在线观看| 波多野结衣av一区二区av| 久久亚洲真实| 热re99久久国产66热| 一进一出抽搐gif免费好疼| 久久久久久国产a免费观看| 最新在线观看一区二区三区| 日本黄色视频三级网站网址| av有码第一页| 成人亚洲精品av一区二区| 亚洲 欧美 日韩 在线 免费| 国内毛片毛片毛片毛片毛片| 亚洲在线自拍视频| av天堂在线播放| 大香蕉久久成人网| 亚洲人成网站高清观看| 亚洲av成人不卡在线观看播放网| 18美女黄网站色大片免费观看| 脱女人内裤的视频| 欧美日本亚洲视频在线播放| 一边摸一边抽搐一进一小说| 国产亚洲精品综合一区在线观看 | 国产爱豆传媒在线观看 | 亚洲五月婷婷丁香| 亚洲国产欧美网| 亚洲片人在线观看| 一进一出抽搐动态| 女人爽到高潮嗷嗷叫在线视频| 男女床上黄色一级片免费看| 久久草成人影院| 亚洲电影在线观看av| 久久久精品欧美日韩精品| 亚洲成国产人片在线观看| 国产亚洲av高清不卡| 亚洲av日韩精品久久久久久密| 色综合欧美亚洲国产小说| 亚洲熟女毛片儿| 久久久国产精品麻豆| 特大巨黑吊av在线直播 | 性色av乱码一区二区三区2| 色尼玛亚洲综合影院| 18禁观看日本| 18美女黄网站色大片免费观看| 国产精品一区二区免费欧美| 巨乳人妻的诱惑在线观看| 最近最新中文字幕大全免费视频| 欧美另类亚洲清纯唯美| 亚洲一区二区三区色噜噜| 丝袜美腿诱惑在线| 九色国产91popny在线| 日韩欧美国产一区二区入口| 男女视频在线观看网站免费 | 后天国语完整版免费观看| 亚洲av中文字字幕乱码综合 | 变态另类丝袜制服| 国产视频内射| 久热爱精品视频在线9| 国产亚洲欧美在线一区二区| 久久久久精品国产欧美久久久| 老熟妇仑乱视频hdxx| 最近最新中文字幕大全免费视频| 亚洲中文日韩欧美视频| 午夜福利18| 99国产精品一区二区蜜桃av| 99国产精品99久久久久| 亚洲精品久久成人aⅴ小说| 亚洲成人久久爱视频| 色婷婷久久久亚洲欧美| 国产精品九九99| 亚洲国产精品久久男人天堂| 无遮挡黄片免费观看| 国产极品粉嫩免费观看在线| 国产黄片美女视频| 日韩三级视频一区二区三区| 国产av一区在线观看免费| 午夜免费成人在线视频| 久久久久国内视频| 好看av亚洲va欧美ⅴa在| 亚洲成人国产一区在线观看| 精品日产1卡2卡| 12—13女人毛片做爰片一| videosex国产| 久久伊人香网站| 女人高潮潮喷娇喘18禁视频| 真人一进一出gif抽搐免费| 国产成+人综合+亚洲专区| 亚洲成av片中文字幕在线观看| 成年女人毛片免费观看观看9| 精品国产亚洲在线| 亚洲aⅴ乱码一区二区在线播放 | 搞女人的毛片| 欧美日韩中文字幕国产精品一区二区三区| 级片在线观看| 精品久久久久久,| 久久久久九九精品影院| 成人精品一区二区免费| 久久久精品欧美日韩精品| 成人18禁高潮啪啪吃奶动态图| 精品久久久久久,| 麻豆国产av国片精品| 色老头精品视频在线观看| 国产一区二区激情短视频| 亚洲国产精品久久男人天堂| 在线视频色国产色| 99在线视频只有这里精品首页| 亚洲av五月六月丁香网| 精品久久久久久久毛片微露脸| 99国产精品99久久久久| 亚洲久久久国产精品| 少妇的丰满在线观看| 亚洲电影在线观看av| 少妇裸体淫交视频免费看高清 | 久久伊人香网站| 久久欧美精品欧美久久欧美| 操出白浆在线播放| 亚洲精品国产区一区二| 亚洲色图 男人天堂 中文字幕| 国产黄片美女视频| а√天堂www在线а√下载| 国产成+人综合+亚洲专区| 欧美激情极品国产一区二区三区| 夜夜躁狠狠躁天天躁| 国产又色又爽无遮挡免费看| 国产激情久久老熟女| 国产亚洲精品第一综合不卡| 久久久精品国产亚洲av高清涩受| 亚洲专区中文字幕在线| 亚洲五月天丁香| 日本熟妇午夜| 亚洲性夜色夜夜综合| 国产在线精品亚洲第一网站| 99在线人妻在线中文字幕| 女人被狂操c到高潮| 美女国产高潮福利片在线看| 丝袜美腿诱惑在线| 国产成人精品久久二区二区免费| 亚洲成人精品中文字幕电影| 国产熟女午夜一区二区三区| 国产99白浆流出| 午夜精品在线福利| 久久国产亚洲av麻豆专区| 女生性感内裤真人,穿戴方法视频| 最近最新中文字幕大全电影3 | tocl精华| 在线观看免费视频日本深夜| 国产三级黄色录像| netflix在线观看网站| 亚洲成a人片在线一区二区| 亚洲 欧美 日韩 在线 免费| 久久天躁狠狠躁夜夜2o2o| 无限看片的www在线观看| 香蕉丝袜av| 人人妻人人看人人澡| 免费高清视频大片| 中文字幕久久专区| 国产私拍福利视频在线观看| √禁漫天堂资源中文www| 久久精品国产亚洲av香蕉五月| 久久久精品国产亚洲av高清涩受| 久久久久久九九精品二区国产 | 国内精品久久久久久久电影| av天堂在线播放| 日韩有码中文字幕| 国产午夜福利久久久久久| cao死你这个sao货| 免费观看精品视频网站| 久久久久久国产a免费观看| 成人国语在线视频| 日本五十路高清| 久9热在线精品视频| 麻豆成人av在线观看| 一夜夜www| av电影中文网址| 久久香蕉国产精品| 怎么达到女性高潮| 国产一区二区在线av高清观看| 亚洲人成网站高清观看| 国产精品久久久久久精品电影 | 老司机在亚洲福利影院| videosex国产| 夜夜夜夜夜久久久久| 嫩草影视91久久| 欧美人与性动交α欧美精品济南到| 在线观看www视频免费| 精品久久久久久久末码| 欧美一级a爱片免费观看看 | 婷婷六月久久综合丁香| 久久精品成人免费网站| 丰满人妻熟妇乱又伦精品不卡| 亚洲 欧美 日韩 在线 免费| 丝袜美腿诱惑在线| 日本一本二区三区精品| 国产爱豆传媒在线观看 | 黄色女人牲交| 麻豆成人午夜福利视频| 色婷婷久久久亚洲欧美| 成人三级黄色视频| 国产精品精品国产色婷婷| 在线av久久热| 欧美精品亚洲一区二区| 欧美激情 高清一区二区三区| 国产成人av教育| 欧美精品啪啪一区二区三区| 波多野结衣高清无吗| 19禁男女啪啪无遮挡网站| 最新在线观看一区二区三区| 亚洲专区国产一区二区| 久久久久久久午夜电影| 精品久久蜜臀av无| av免费在线观看网站| 国产av不卡久久| 制服诱惑二区| 国产熟女xx| av在线播放免费不卡| 91麻豆精品激情在线观看国产| 国产成人精品久久二区二区91| 中文字幕久久专区| 成人三级黄色视频| 欧美色视频一区免费| 成人永久免费在线观看视频| 99国产极品粉嫩在线观看| 国产欧美日韩精品亚洲av| 久久这里只有精品19| 91九色精品人成在线观看| 成人av一区二区三区在线看| 国产精品爽爽va在线观看网站 | 在线观看一区二区三区| 国产午夜精品久久久久久| 久久草成人影院| 日本黄色视频三级网站网址| 亚洲七黄色美女视频| 亚洲精品国产一区二区精华液| 999久久久精品免费观看国产| 黑丝袜美女国产一区| 免费看日本二区| 黄色成人免费大全| 国内毛片毛片毛片毛片毛片| svipshipincom国产片| 午夜视频精品福利| 男女下面进入的视频免费午夜 | 国产成年人精品一区二区| 这个男人来自地球电影免费观看| 男男h啪啪无遮挡| 欧美午夜高清在线| e午夜精品久久久久久久| 一本大道久久a久久精品| 久久精品国产清高在天天线| 国产在线观看jvid| 久久精品国产亚洲av高清一级| 人妻丰满熟妇av一区二区三区| 日韩欧美国产一区二区入口| 久久精品成人免费网站| 女人爽到高潮嗷嗷叫在线视频| 久久久久久久久久黄片| 99久久精品国产亚洲精品| 9191精品国产免费久久| 国产亚洲精品久久久久久毛片| 极品教师在线免费播放| 欧美一级a爱片免费观看看 | 一级毛片精品| 久久婷婷人人爽人人干人人爱| 后天国语完整版免费观看| 亚洲avbb在线观看| 天天躁夜夜躁狠狠躁躁| 老鸭窝网址在线观看| 亚洲国产中文字幕在线视频| 男女午夜视频在线观看| 午夜免费激情av| 国产精品亚洲美女久久久| 男女之事视频高清在线观看| 好看av亚洲va欧美ⅴa在| 国产亚洲精品久久久久5区| 一本综合久久免费| 亚洲av美国av| 精品午夜福利视频在线观看一区| 99国产极品粉嫩在线观看| 免费在线观看黄色视频的| 亚洲午夜精品一区,二区,三区| 久久人妻福利社区极品人妻图片| 久久精品国产综合久久久| 日韩欧美一区视频在线观看| 日本 欧美在线| 波多野结衣巨乳人妻| 18禁国产床啪视频网站| 成年女人毛片免费观看观看9| 麻豆成人av在线观看| 91av网站免费观看| 脱女人内裤的视频| 18禁国产床啪视频网站| 人妻久久中文字幕网| 男人舔奶头视频| 欧美成狂野欧美在线观看| 曰老女人黄片| av在线天堂中文字幕| 国内精品久久久久精免费| 在线观看一区二区三区| 国产区一区二久久| 午夜两性在线视频| 51午夜福利影视在线观看| 一本精品99久久精品77| 亚洲自拍偷在线| 1024手机看黄色片| 亚洲精品一卡2卡三卡4卡5卡| 18禁国产床啪视频网站| 亚洲在线自拍视频| 久久久久国产精品人妻aⅴ院| 国产91精品成人一区二区三区| 亚洲第一欧美日韩一区二区三区| 日韩免费av在线播放| 日韩一卡2卡3卡4卡2021年| 欧美三级亚洲精品| 香蕉丝袜av| 我的亚洲天堂| 亚洲av美国av| 日日爽夜夜爽网站| 黄片小视频在线播放| 在线av久久热| 久久精品人妻少妇| 一本久久中文字幕| 免费在线观看视频国产中文字幕亚洲| 欧美性猛交╳xxx乱大交人| 欧美乱码精品一区二区三区| av有码第一页| 脱女人内裤的视频| 亚洲国产欧洲综合997久久, | 亚洲国产高清在线一区二区三 | 18禁裸乳无遮挡免费网站照片 | 日本成人三级电影网站| 国产精品九九99| 免费观看精品视频网站| 日本精品一区二区三区蜜桃| 成人18禁高潮啪啪吃奶动态图| 草草在线视频免费看| 欧美一级毛片孕妇| 男人舔女人的私密视频| 香蕉丝袜av| 99久久99久久久精品蜜桃| 1024香蕉在线观看| 亚洲欧美日韩高清在线视频| 免费电影在线观看免费观看| 日本在线视频免费播放| 草草在线视频免费看| 久久国产精品人妻蜜桃| 精品国产美女av久久久久小说| 日韩欧美在线二视频| 欧美在线黄色| 亚洲免费av在线视频| 精品电影一区二区在线| 欧美成人一区二区免费高清观看 | 伊人久久大香线蕉亚洲五| 欧美日韩福利视频一区二区| 亚洲 欧美 日韩 在线 免费| 亚洲精品久久国产高清桃花| 一级a爱视频在线免费观看| 亚洲精品中文字幕一二三四区| 中文字幕久久专区| 在线免费观看的www视频| 国产精品久久电影中文字幕| www.自偷自拍.com| 午夜福利视频1000在线观看| 亚洲专区国产一区二区| 啪啪无遮挡十八禁网站| 婷婷亚洲欧美|