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

    桂花類胡蘿卜素異構(gòu)酶基因的克隆及表達(dá)分析

    2017-06-21 12:28:47張超劉玉成王藝光付建新趙宏波
    生物技術(shù)通報 2017年6期
    關(guān)鍵詞:金球花序胡蘿卜素

    張超 劉玉成 王藝光 付建新 趙宏波

    (浙江農(nóng)林大學(xué)風(fēng)景園林與建筑學(xué)院,臨安 311300)

    桂花類胡蘿卜素異構(gòu)酶基因的克隆及表達(dá)分析

    張超 劉玉成 王藝光 付建新 趙宏波

    (浙江農(nóng)林大學(xué)風(fēng)景園林與建筑學(xué)院,臨安 311300)

    利用已構(gòu)建的桂花‘堰虹桂’花瓣轉(zhuǎn)錄組數(shù)據(jù)庫中Unigene序列信息,結(jié)合RT-PCR技術(shù)對1個桂花類胡蘿卜素異構(gòu)酶基因(OfCRTISO)的最大閱讀框進(jìn)行擴(kuò)增,并測序驗證。生物信息學(xué)分析表明OfCRTISO基因含有長為1 842 bp的開放讀碼框,編碼613個氨基酸殘基。桂花OfCRTISO與已知胡麻科和茄科植物的CRTISO序列相似性最高,系統(tǒng)進(jìn)化樹分析發(fā)現(xiàn)OfCRTISO與芝麻CRTISO(XP_011077785)親緣關(guān)系最近,聚為同一小支。表達(dá)分析發(fā)現(xiàn),OfCRTISO基因在桂花‘堰虹桂’花蕾前期花序中的表達(dá)量極低,花蕾后期其表達(dá)量逐漸增加,初開期達(dá)到高峰,而盛開期時,其表達(dá)量顯著下降。對桂花不同花色品種花瓣中OfCRTISO基因的表達(dá)分析表明,桂花初開期OfCRTISO基因的表達(dá)量可能影響盛開期各品種花瓣中類胡蘿卜素的積累。

    桂花;花色;類胡蘿卜素;類胡蘿卜素異構(gòu)酶基因

    桂花(Osmanthus fragrans Lour.)是我國十大傳統(tǒng)名花之一,也是我國重要的園林植物之一。根據(jù)花期桂花可分為秋桂和四季桂兩大類,秋桂又可根據(jù)花色分為金桂(黃色至金黃色系)、銀桂(黃白色至黃色系)和丹桂(橙色至橙紅色系)三大品種群[1]。前人通過對秋桂不同花色品種花瓣呈色物質(zhì)成分進(jìn)行分析發(fā)現(xiàn),類胡蘿卜素(Carotenoid)的種類及其含量是影響秋桂花色呈現(xiàn)的最主要因素[2,3]。其中,葉黃素(Lutein)、α-胡蘿卜素(α-carotene)和β-胡蘿卜素(β-carotene)這3種類胡蘿卜素的含量與秋桂花瓣顏色的呈現(xiàn)緊密相關(guān)[3]。

    類胡蘿卜素屬于萜烯類化合物,呈現(xiàn)黃色、橙紅色和紅色,是觀賞植物花瓣的重要呈色物質(zhì)之一[4]。根據(jù)類胡蘿卜素的化學(xué)結(jié)構(gòu),主要分為胡蘿卜素和葉黃素類化合物兩大類[5]。到目前為止,通過各種生物學(xué)方法和手段已經(jīng)分離和鑒定了參與植物類胡蘿卜素生物合成的全部基因[4,6]。

    類胡蘿卜素異構(gòu)酶(Carotenoid isomerase,CRTISO)能催化前番茄紅素(7,9,7',9'-四順式-番茄紅素,7,9,7',9'-tetra-cis-lycopene)生成全反式番茄紅素(All-trans-lycopene)[7,8]。植物CRTISO基因首先在擬南芥(Arabidopsis thaliana)[9]和番茄(Lycopersicon esculentum)[10]中被分離和鑒定。Park等[9]在擬南芥中定位并克隆了CRTISO基因,將其轉(zhuǎn)化大腸桿菌,證實(shí)CRTISO具有催化多順式類胡蘿卜素向全反式類胡蘿卜素異構(gòu)化的功能。相對于野生型番茄,Isaacson等[10]發(fā)現(xiàn)CRTISO缺失突變體tangerinemic和CRTISO低豐度表達(dá)突變體tangerine3183果實(shí)中積累大量前番茄紅素和ζ-胡蘿卜素,而番茄紅素和β-胡蘿卜素含量較低。

    此后,除了模式植物和園藝作物[11-13]以外,許多觀賞植物的CRTISO基因也相繼被分離,如菊花(Chrysanthemum morifolium)[14]、 百 合(Lilium spp.)[15]、金盞菊(Calendula officinalis)[16]和香石竹(Dianthus caryophyllus)[17]。通過對菊花的研究發(fā)現(xiàn),富含葉黃素的黃花品種的CRTISO基因表達(dá)量顯著高于葉黃素含量低的白花品種[14]。另有研究發(fā)現(xiàn),不同牽牛屬(Ipomoea)植物和百合不同品種番茄紅素下游類胡蘿卜素的含量與CRTISO基因表達(dá)量緊密相關(guān)[15,18]。

    本研究利用已構(gòu)建的桂花‘堰虹桂’花瓣轉(zhuǎn)錄組數(shù)據(jù)庫[19]中Unigene序列信息,結(jié)合RT-PCR技術(shù)克隆得到桂花OfCRTISO基因cDNA全長,對其基因序列及其編碼氨基酸序列進(jìn)行分析,并采用實(shí)時熒光定量PCR技術(shù)對該基因在‘堰虹桂’不同發(fā)育階段花序和營養(yǎng)器官中以及不同花色品種花瓣中的表達(dá)模式進(jìn)行分析,旨在為進(jìn)一步研究桂花OfCRTISO基因的功能,探討桂花花色形成的分子機(jī)制提供理論依據(jù)。

    1 材料與方法

    1.1 材料

    以浙江農(nóng)林大學(xué)桂花資源圃地栽桂花丹桂品種‘堰虹桂’、銀桂品種‘小葉蘇桂’和金桂品種‘金球桂’為試驗材料。

    不同發(fā)育階段花序(圖1)取樣:2014年8月9日、8月15日、8月22日、8月29日、9月4日、9月10日、9月17日、9月23日、9月24日和9月25日對‘堰虹桂’不同發(fā)育狀態(tài)花序進(jìn)行取樣。其中,8月9日至9月17日花序狀態(tài)為花蕾期(Bud stage),簡寫為B1-B7;9月23日花序狀態(tài)為鈴梗期(Linggeng stage,L);9月24日花序狀態(tài)為初開期(Initial flowering stage,IF);9月25日花序狀態(tài)為盛開期(Full flowering stage,F(xiàn)F)。取樣時間均為上午10時,花序樣品離體后速凍于液氮中,于-80℃儲藏備用。

    營養(yǎng)器官取樣:2014年9月25日對‘堰虹桂’一年生莖(當(dāng)年生枝條先端第2至第3節(jié)之間的莖)、嫩葉(當(dāng)年生枝條先端第1對葉片)和成熟葉(當(dāng)年生枝條先端第3對葉片)進(jìn)行取樣,液氮速凍后,于-80℃儲藏備用。

    不同花色品種花瓣取樣:2014年9月下旬對‘堰虹桂’、‘小葉蘇桂’和‘金球桂’鈴梗期、初開期和盛開期花序的花瓣進(jìn)行取樣,取樣時間均為上午10時,離體后速凍于液氮中,于-80℃儲藏備用。

    圖1 桂花‘堰虹桂’不同發(fā)育階段花序

    1.2 方法

    1.2.1 總RNA提取與cDNA的第一鏈的合成 采用RNAprep pure Plant Kit(天根,北京)按照說明書步驟提取各樣品的總RNA。紫外分光光度計和1%瓊脂糖凝膠電泳檢測總RNA濃度和質(zhì)量。采用Reverse Transcriptase M-MLV(TaKaRa,大連)按照說明書合成cDNA的第一鏈后,儲存于-20℃?zhèn)溆谩?/p>

    1.2.2 OfCRTISO基因的克隆 根據(jù)前期課題組已構(gòu)建的桂花‘堰虹桂’花瓣轉(zhuǎn)錄組數(shù)據(jù)庫[19],從中篩選出與CRTISO基因相關(guān)序列Locus_7474,設(shè)計特異性引物F1(5'- TTACTTTTACTTGTTAACCACCG -3') 和 R1(5'- GCAATCACTATTTCTAATCCTCT -3'),以‘堰虹桂’盛開期花瓣RNA反轉(zhuǎn)的cDNA為模板,對桂花OfCRTISO基因最大閱讀框序列進(jìn)行PCR擴(kuò)增。PCR反應(yīng)體系:10×Buffer 2.5 μL,25 mmol/L MgCl21.5 μL,10 mmol/L dNTP 0.5 μL,10 μmol/L上、下游引物各0.5 μL,Taq DNA聚合酶1.25 U,cDNA 2 μL,加去離子水至25 μL。PCR反應(yīng)條件:94℃預(yù)變性5 min,94℃變性30 s、60℃退火30 s、72℃延伸2.5 min,35個循環(huán),最后72℃保溫10 min,4℃結(jié)束反應(yīng)。

    PCR反應(yīng)產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢測后回收、純化,與pMD18-T載體(TaKaRa)連接,轉(zhuǎn)化大腸桿菌DH5α感受態(tài)細(xì)胞,藍(lán)白斑篩選陽性克隆,經(jīng)PCR鑒定后送上海生工生物工程技術(shù)服務(wù)有限公司測序。

    1.2.3 序列分析 測序獲得的序列通過NCBI提供的ORF Finder進(jìn)行ORF查找,序列翻譯、氨基酸序列比對采用DNAMAN 6.0軟件進(jìn)行。通過ExPASy(http://www.expasy.org/tools/protparam.html)提供的Prot-Param在線軟件對蛋白質(zhì)的分子量和等電點(diǎn)進(jìn)行分析。通過NetPhos 2.0(http://www.cbs. dtu.dk/services/NetPhos/)對蛋白磷酸化位點(diǎn)進(jìn)行預(yù)測;使用Plant-PLoc軟件(http://www.csbio.sjtu.edu. cn/bioinf/plant/)對編碼蛋白的亞細(xì)胞定位進(jìn)行預(yù)測;使用MEGA 5.05軟件中的 Neighbor-Joining(鄰位相連法,NJ)法建立系統(tǒng)進(jìn)化樹。

    1.2.4 實(shí)時熒光定量PCR 按照Takara公司的SYBR?Premix Ex TaqTMII(Tli RNaseH Plus) 試 劑盒的說明書,采用7300 real-time PCR儀(Applied Biosystems)檢測桂花OfCRTISO基因的相對表達(dá)量。OfCRTISO基因表達(dá)引物為F2(5'- GAAAAACAAGGGATTCTCGGA -3')和R2(5'- GCAGATAGTAGGCAAGGGTCAA -3'),內(nèi)參基因OfACT表達(dá)引物為F3(5'- CCCAAGGCAAACAGAGAAAAAAT -3') 和R3(5'- ACCCCATCACCAGAATCAAGAA -3')[20]。所用反應(yīng)體系和程序如下,PCR反應(yīng)體系為20 μL,其 中 2×SYBR Premix Ex Taq II(Tli RNaseH Plus)10 μL,10 μmol/L上下游定量引物各0.8 μL,模板cDNA 2 μL,50×ROX Reference Dye 0.4 μL,用滅菌雙蒸水補(bǔ)齊至20 μL。熒光定量PCR擴(kuò)增采用兩步法,即在95℃預(yù)變性30 s之后,先運(yùn)行40個循環(huán)的95℃ 5 s,60℃ 31 s,再運(yùn)行熔解曲線階段的95℃15 s,60℃ 1 min,95℃ 30 s,60℃ 15 s。反應(yīng)后根據(jù)熔解曲線分析PCR產(chǎn)物的特異性,分別導(dǎo)出目的基因和內(nèi)參基因的Ct值,利用參照基因的ΔCt法計算目的基因的相對表達(dá)量。

    2 結(jié)果

    2.1 桂花OfCRTISO基因的克隆與序列分析

    提取所有樣品總RNA均有28S和18S兩條清晰條帶,RNA無降解,部分樣品總RNA經(jīng)1%瓊脂糖凝膠電泳檢測如圖2所示。紫外分光光度計檢測OD260/OD280均在1.9-2.0之間,濃度為210-420 ng/μL之間。以‘堰虹桂’盛開期花瓣總RNA反轉(zhuǎn)的cDNA為模板,擴(kuò)增出一條約2 200 bp的單一特異性條帶(圖3)。經(jīng)測序獲得一條長度為2 241 bp的cDNA序列(圖4),Blast比對發(fā)現(xiàn)該序列與CRTISO基因同源,命名為OfCRTISO。

    利用NCBI提供的ORF Finder進(jìn)行分析發(fā)現(xiàn),OfCRTISO包含一個長度為1 842 bp的開放讀碼框,5'非翻譯區(qū)長50 bp,3'非翻譯區(qū)長349 bp。其開放讀碼框編碼一個含有613個氨基酸殘基的蛋白質(zhì)(圖4)。

    2.2 桂花OfCRTISO基因編碼氨基酸的基本特征

    利用NCBI Conserved Domain Search進(jìn)行保守域分析表明,OfCRTISO編碼的蛋白含有1個carot_ isom保守域(第96-590位氨基酸殘基)。利用Prot-Param軟件分析桂花OfCRTISO氨基酸序列的基本特征發(fā)現(xiàn),該蛋白的分子式為C3079H4818N816O875S24,相對分子量為68.04 kD,理論等電點(diǎn)(pI)為8.58。該蛋白中亮氨酸(Leu)最多占10.4%;色氨酸(Trp)含量最低,只占0.8%。

    利用NetPhos 2.0軟件對OfCRTISO蛋白進(jìn)行磷酸化位點(diǎn)的分析(圖4)發(fā)現(xiàn),該蛋白含有14個絲氨酸激酶磷酸化位點(diǎn)(分別為第72、85、91、149、155、200、238、282、363、442、473、481、498和531位氨基酸),5個蘇氨酸激酶磷酸化位點(diǎn)(分別為第299、382、425、519和522位氨基酸)和7個酪氨酸激酶磷酸化位點(diǎn)(分別為第86、127、244、343、403、485和534位氨基酸)。經(jīng)Plant-PLoc軟件分析,推測OfCRTISO蛋白定位于葉綠體中。

    2.3 桂花OfCRTISO基因編碼氨基酸的同源性比較

    利用DNAMAN軟件對桂花OfCRTISO及其他11個物種的CRTISO氨基酸序列進(jìn)行了多重比對,發(fā)現(xiàn)桂花OfCRTISO與已知胡麻科和茄科植物的CRTISO氨基酸序列有很高的同源性(圖5),其中與芝麻(Sesamum indicum)CRTISO序列(XP_011077785)相似性最高,序列一致性為81%。

    為進(jìn)一步了解OfCRTISO與其他植物CRTISO之間的進(jìn)化關(guān)系,用上述11個物種的CRTISO氨基酸序列與桂花OfCRTISO構(gòu)建了系統(tǒng)進(jìn)化樹(圖6)發(fā)現(xiàn),桂花OfCRTISO與其他雙子葉植物CRTISO聚為一大支,而水稻(Oryza sativa)和玉米(Zea mays)等單子葉植物則聚在另一大支。馬鈴薯(Solanum tuberosum)、番茄(Solanum lycopersicum)和枸杞(Lycium barbarum)等茄科植物聚為一小支,表明CRTISO氨基酸的進(jìn)化與植物的進(jìn)化基本一致。桂花OfCRTISO與芝麻CRTISO聚為同一小支,表明其具有更近的親緣關(guān)系(圖6)。

    2.4 桂花OfCRTISO基因的表達(dá)分析

    桂花OfCRTISO基因在‘堰虹桂’營養(yǎng)器官(一年生莖、嫩葉和成熟葉)和不同發(fā)育階段花序中的表達(dá)模式如圖7所示,結(jié)果表明,OfCRTISO基因在花蕾前期(B1-B5)花序中的表達(dá)量極低,花蕾后期(B6-B7)時其表達(dá)量逐漸增加。隨著花瓣展開,花序中OfCRTISO基因的表達(dá)量急劇增大,于初開期達(dá)到最大值;而當(dāng)花序開放至盛開期時,其表達(dá)量顯著下降。OfCRTISO基因在一年生莖、嫩葉和成熟葉片中的表達(dá)量無顯著差異,與B6時期花序中表達(dá)量相近。

    隨著花序開放過程(鈴梗期至盛開期),OfCRTISO基因在銀桂品種‘小葉蘇桂’花瓣中的表達(dá)量不變,在金桂品種‘金球桂’花瓣中表達(dá)量逐漸增加,而在丹桂品種‘堰虹桂’花瓣中的表達(dá)量呈現(xiàn)先增加再減少的趨勢(圖8)。

    OfCRTISO基因在鈴梗期不同花色桂花品種花瓣中的表達(dá)量相近;初開期時,‘堰虹桂’花瓣中OfCRTISO基因的表達(dá)量顯著高于‘金球桂’和‘小葉蘇桂’,‘金球桂’花瓣中的OfCRTISO基因的表達(dá)量顯著高于‘小葉蘇桂’;盛開期‘金球桂’和‘堰虹桂’花瓣中OfCRTISO基因的表達(dá)量相近,均顯著高于‘小葉蘇桂’(圖8)。

    圖2 桂花部分樣品總RNA電泳圖

    圖3 桂花OfCRTISO基因的PCR擴(kuò)增電泳圖

    3 討論

    對于桂花而言,雖然已明確類胡蘿卜素的種類與含量是影響其花色呈現(xiàn)的最主要因素[21],但是,桂花類胡蘿卜素合成的分子機(jī)制仍不清晰[2,3,19]。本研究從桂花中獲得了CRTISO基因cDNA全長序列,并分析該基因在‘堰虹桂’不同發(fā)育階段花序和營養(yǎng)器官中以及不同花色品種花瓣中的表達(dá)模式,為探討桂花花色形成的分子機(jī)制奠定基礎(chǔ)。

    圖4 桂花OfCRTISO基因的核苷酸序列及推測的氨基酸序列

    前人研究證實(shí)CRTISO基因在擬南芥和番茄基因組中均以單拷貝形式存在[9,10],本研究利用已構(gòu)建的桂花轉(zhuǎn)錄組數(shù)據(jù)庫只獲得1個桂花CRTISO基因序列,推測CRTISO基因在桂花基因組中也可能以單拷貝形式存在。氨基酸序列對比和系統(tǒng)進(jìn)化樹分析表明桂花OfCRTISO與芝麻CRTISO(XP_011077785)序列相似性最高,兩者親緣關(guān)系最近(圖6)。

    圖5 桂花OfCRTISO與其他物種CRTISO氨基酸序列對比

    圖6 桂花OfCRTISO與其他物種CRTISO氨基酸序列的進(jìn)化樹分析

    CRTISO可以催化前番茄紅素生成全反式番茄紅素[7,8],是植物類胡蘿卜素生物合成過程中一個重要的催化酶。觀賞植物花瓣中類胡蘿卜素合成與CRTISO基因的表達(dá)量緊密相關(guān)[14,15,18]。與菊花品種‘Yellow Paragon’[14]和百合品種‘Connecticut King’[15]相似,桂花‘堰虹桂’花序發(fā)育過程中,類胡蘿卜素總量逐漸增加,至盛開期時達(dá)到最大[19]。本研究發(fā)現(xiàn),‘堰虹桂’花序和花瓣中OfCRTISO基因的表達(dá)量均于初開期達(dá)到最大,但是盛開期時,其表達(dá)量顯著下降。這與百合品種‘Connecticut King’花蕾發(fā)育過程中LhCRTISO基因的表達(dá)模式[15]一致,推測桂花‘堰虹桂’初開期OfCRTISO基因的高豐度表達(dá)對盛開期類胡蘿卜素的大量積累至關(guān)重要。

    百合、香石竹和菊花不同花色品種開花早期CRTISO基因的表達(dá)量均無差異[6,14,15],這與本研究結(jié)果相似,桂花不同花色品種鈴梗期時花瓣中OfCRTISO基因的表達(dá)量也相同。而開花中期,菊花黃色花品種‘Yellow Paragon’類胡蘿卜素總量和CRTISO基因的表達(dá)量則顯著高于白色花品種‘Paragon’[14]。同樣的,百合品種‘Connecticut king’ LhCRTISO基因的表達(dá)量于開花中期顯著高于百合品種‘Montreux’,而在該時期‘Connecticut king’類胡蘿卜素總量也顯著高于‘Montreux’[15],表明CRTISO基因的表達(dá)量與類胡蘿卜素的積累緊密相關(guān)。在本研究中,初開期時‘堰虹桂’和‘金球桂’花瓣中OfCRTISO基因的表達(dá)量急劇上升,均顯著高于‘小葉蘇桂’中的表達(dá)量,而‘堰虹桂’中OfCRTISO基因的表達(dá)量上升更為顯著。侯丹測定‘堰虹桂’、‘金球桂’和‘小葉蘇桂’花瓣中類胡蘿卜素含量發(fā)現(xiàn),‘堰虹桂’花瓣中類胡蘿卜素總量最高,‘金球桂’次之,而‘小葉蘇桂’最低[21],這可能與3個品種初開期OfCRTISO基因的表達(dá)量正相關(guān)。盛開期時,‘堰虹桂’和‘金球桂’花瓣中OfCRTISO基因的表達(dá)量相近,而類胡蘿卜素含量存在差異,表明桂花初開期OfCRTISO基因的表達(dá)高低對盛開期花瓣中類胡蘿卜素積累的影響更大。

    圖7 OfCRTISO基因在‘堰虹桂’營養(yǎng)器官和不同發(fā)育階段花序中的表達(dá)量

    圖8 OfCRTISO基因在不同花色桂花品種花瓣中的表達(dá)量

    4 結(jié)論

    桂花OfCRTISO基因含有長為1 842 bp的開放讀碼框,編碼613個氨基酸殘基。桂花OfCRTISO與已知胡麻科和茄科植物的CRTISO序列相似性最高。OfCRTISO基因在桂花‘堰虹桂’初開期表達(dá)量最高,初開期OfCRTISO基因的表達(dá)量影響盛開期各品種花瓣中類胡蘿卜素的積累。

    [1] 向其柏, 劉玉蓮. 中國桂花品種圖志[M]. 杭州:浙江科學(xué)技術(shù)出版社, 2008.

    [2] Han Y, Li L, Dong M, et al. cDNA cloning of the phytoene synthase(PSY)and expression analysis of PSY and carotenoid cleavage dioxygenase genes in Osmanthus fragrans[J]. Biologia, 2013, 68(2):258-263.

    [3] Han Y, Wang X, Chen W, et al. Differential expression of carotenoidrelated genes determines diversified carotenoid coloration in flowerpetal of Osmanthus fragrans[J]. Tree Genetics & Genomes, 2014, 10(2):329-338.

    [4] Zhu C, Bai C, Sanahuja G, et al. The regulation of carotenoid pigmentation in flowers[J]. Archives of Biochemistry and Biophysics, 2010, 504(1):132-141.

    [5] Tanaka Y, Ohmiya A. Seeing is believing:engineering anthocyanin and carotenoid biosynthetic pathways[J]. Current Opinion in Biotechnology, 2008, 19(2):190-197.

    [6] Ohmiya A. Qualitative and quantitative control of carotenoid accumulation in flower petals[J]. Scientia Horticulturae, 2013, 163(6):10-19.

    [7] Breitenbach J, Vioque A, Sandmann G. Gene sll0033 from Synechocystis 6803 encodes a carotene isomerase involved in the biosynthesis of all-E Lycopene[J]. Zeitschrift für Naturforschung Section C-A Journal of Biosciences, 2001, 56(9-10):915-917.

    [8] Masamoto K, Wada H, Kaneko T, et al. Identification of a gene required for cis-to-trans carotene isomerization in carotenogenesis of the cyanobacterium Synechocystis sp. PCC 6803[J]. Plant & Cell Physiology, 2001, 42(12):1398-1402.

    [9] Park H, Kreunen SS, Cuttriss AJ, et al. Identification of the carotenoid isomerase provides insight into carotenoid biosynthesis, prolamellar body formation, and photomorphogenesis[J]. Plant Cell, 2002, 14(2):321-332.

    [10] Isaacson T, Ronen G, Zamir D, et al. Cloning of tangerine from tomato reveals a carotenoid isomerase essential for the production of beta-carotene and xanthophylls in plants[J]. Plant Cell, 2002, 14(2):333-342.

    [11] Chai C, Fang J, Liu Y, et al. ZEBRA2, encoding a carotenoid isomerase, is involved in photoprotection in rice[J]. Plant Molecular Biology, 2011, 75(3):211-221.

    [12] Kato M, Ikoma Y, Matsumoto H, et al. Accumulation of carotenoids and expression of carotenoid biosynthetic genes during maturation in citrus fruit[J]. Plant Physiology, 2004, 134(2):824-837.

    [13] Liu Y, Zeng S, Sun W, et al. Comparative analysis of carotenoid accumulation in two goji(Lycium barbarum L. and L. ruthenicum Murr. )fruits[J]. BMC Plant Biology, 2014, 14(1):269.

    [14] Kishimoto S, Ohmiya A. Regulation of carotenoid biosynthesis in petals and leaves of chrysanthemum(Chrysanthemum morifolium)[J]. Physiologia Plantarum, 2006, 128(3):436-447.

    [15] Yamagishi M, Kishimoto S, Nakayama M. Carotenoid composition and changes in expression of carotenoid biosynthetic genes in tepals of Asiatic hybrid lily[J]. Plant Breeding, 2010, 129(1):100-107.

    [16] Sanae K, Akemi O. Carotenoid isomerase is key determinant of petal color of Calendula officinalis[J]. Journal of Biological Chemistry, 2012, 287(1):276-285.

    [17] Ohmiya A, Tanase K, Hirashima M, et al. Analysis of carotenogenic gene expression in petals and leaves of carnation(Dianthus caryophyllus L. )[J]. Plant Breeding, 2013, 132(4):423-429.

    [18] Yamamizo C, Kishimoto S, Ohmiya A. Carotenoid composition and carotenogenic gene expression during Ipomoea petal development[J]. Journal of Experimental Botany, 2010, 61(3):709-719.

    [19] Zhang C, Wang Y, Fu J, et al. Transcriptomic analysis and carotenogenic gene expression related to petal coloration in Osmanthus fragrans ‘Yanhong Gui’[J]. Trees-Structure and Function, 2016, 30(4):1207-1223.

    [20] Zhang C, Fu J, Wang Y, et al. Identification of suitable reference genes for gene expression normalization in the quantitative realtime PCR analysis of sweet osmanthus(Osmanthus fragrans Lour.)[J]. PLoS One, 2015, 10(8):e0136355.

    [21] 侯丹. 桂花主要品種花香和花色及其對溫度變化的響應(yīng)[D].臨安:浙江農(nóng)林大學(xué), 2014.

    (責(zé)任編輯 狄艷紅)

    Cloning and Expression Analysis of Carotenoid Isomerase Gene in Osmanthus fragrans

    ZHANG Chao LIU Yu-cheng WANG Yi-guang FU Jian-xin ZHAO Hong-bo
    (School of Landscape Architecture,Zhejiang Agriculture and Forestry University,Lin’an 311300)

    Based on the sequence information of Unigene of Osmanthus fragrans ‘Yanhong Gui’ from the petal transcriptome database,sequence of open reading frame in one carotenoid isomerase gene(OfCRTISO)from Osmanthus fragrans was amplified using RT-PCR technology and confirmed by sequencing. Bioinformatics analysis indicated that OfCRTISO contained a 1 842 bp open reading frame,encoding 613 amino acid residues. OfCRTISO from Osmanthus fragrans shared high sequence similarity with CRTISOs from Pedaliaceae and Solanaceae plants. Phylogenetic analysis showed that OfCRTISO and Sesamum indicum CRTISO(XP_011077785)gathered into one branch,revealing that these two proteins showed the closest genetic relationship. Gene expression analysis implicated that the expression level of OfCRTISO was extremely low at bud former stage,gradually increased at bud later stage,and reached the peak at initial flowering stage;however,it significantly decreased at full flowering stage. According to the expression analysis of OfCRTISO in the petals of sweet osmanthus cultivars with different flower color,it was demonstrated that the expression level of OfCRTISO at initial flowering stage might have an impact on the carotenoid accumulation in the petals of sweet osmanthus cultivars at full flowering stage.

    Osmanthus fragransl;petal color;carotenoids;carotenoid isomerase gene

    10.13560/j.cnki.biotech.bull.1985.2016-1163

    2017-01-05

    國家自然科學(xué)基金項目(31501790),浙江省自然科學(xué)基金(LQ15C160004和LQ16C160003),浙江農(nóng)林大學(xué)科研發(fā)展基金人才啟動項目(2016FR031)

    張超,男,博士,研究方向:觀賞植物分子生物學(xué);E-mail:zhangchao24804@gmail.com

    趙宏波,男,博士,研究方向:觀賞植物遺傳育種;E-mail:zhaohb@zafu.edu.cn

    猜你喜歡
    金球花序胡蘿卜素
    一些花序的新定義和一個新的花序分類系統(tǒng)
    ——植物學(xué)教材質(zhì)疑(六)
    2017年金球獎30人名單
    足球周刊(2017年22期)2018-04-03 07:05:21
    傀儡師
    南風(fēng)(2017年34期)2018-01-16 07:29:09
    花序最大的木本植物
    花序最大的草本植物
    β-胡蘿卜素微乳液的體外抗氧化性初探
    中國果菜(2016年9期)2016-03-01 01:28:39
    決賽無功仍金球亞洲代表遭群嘲
    RP-HPLC法測定螺旋藻中β-胡蘿卜素的含量
    論我國類胡蘿卜素保健食品的審評審批
    類胡蘿卜素與癌癥風(fēng)險相關(guān)性的研究
    人人妻,人人澡人人爽秒播| 三级毛片av免费| √禁漫天堂资源中文www| 国产成人免费无遮挡视频| 国产精品98久久久久久宅男小说| 亚洲国产毛片av蜜桃av| 国产有黄有色有爽视频| 国产精品亚洲av一区麻豆| 99精品在免费线老司机午夜| 老熟妇仑乱视频hdxx| 少妇粗大呻吟视频| 考比视频在线观看| 久久国产精品大桥未久av| 肉色欧美久久久久久久蜜桃| 亚洲成国产人片在线观看| 亚洲av日韩在线播放| 亚洲国产av影院在线观看| 亚洲精品久久午夜乱码| 久久婷婷成人综合色麻豆| 国产免费av片在线观看野外av| 久久人妻熟女aⅴ| 汤姆久久久久久久影院中文字幕| 亚洲伊人久久精品综合| 高清在线国产一区| 超碰成人久久| 婷婷成人精品国产| 美女午夜性视频免费| www.999成人在线观看| 久久精品国产亚洲av高清一级| 成人黄色视频免费在线看| 日韩大片免费观看网站| 亚洲 欧美一区二区三区| 欧美日韩视频精品一区| 亚洲伊人久久精品综合| 91大片在线观看| 亚洲精品久久午夜乱码| 亚洲午夜精品一区,二区,三区| 亚洲精品中文字幕在线视频| 亚洲 国产 在线| 国产真人三级小视频在线观看| 1024香蕉在线观看| 国产一区二区激情短视频| 亚洲精品中文字幕一二三四区 | 一本综合久久免费| 在线观看一区二区三区激情| 满18在线观看网站| 如日韩欧美国产精品一区二区三区| 国产精品秋霞免费鲁丝片| 日韩中文字幕视频在线看片| 国产男女内射视频| 一级毛片电影观看| 国产av国产精品国产| 一区二区三区激情视频| 国产在线视频一区二区| 欧美乱妇无乱码| 精品久久久久久电影网| 一本色道久久久久久精品综合| 久久99热这里只频精品6学生| 国产无遮挡羞羞视频在线观看| 国产91精品成人一区二区三区 | 不卡一级毛片| 欧美激情 高清一区二区三区| 黄片小视频在线播放| 亚洲国产欧美一区二区综合| 午夜久久久在线观看| 国产精品国产高清国产av | 满18在线观看网站| 国产人伦9x9x在线观看| 一边摸一边做爽爽视频免费| 91字幕亚洲| 色在线成人网| 国产有黄有色有爽视频| 免费av中文字幕在线| 亚洲专区国产一区二区| 欧美老熟妇乱子伦牲交| 日本av免费视频播放| 我的亚洲天堂| 999精品在线视频| 国产在线视频一区二区| 免费在线观看黄色视频的| 91精品国产国语对白视频| 三上悠亚av全集在线观看| 国产午夜精品久久久久久| 精品熟女少妇八av免费久了| 一本大道久久a久久精品| 亚洲欧美色中文字幕在线| 满18在线观看网站| 久久久久久久久免费视频了| 乱人伦中国视频| 99精国产麻豆久久婷婷| 多毛熟女@视频| 在线观看免费视频日本深夜| 五月开心婷婷网| 一进一出抽搐动态| 欧美国产精品va在线观看不卡| 啦啦啦免费观看视频1| cao死你这个sao货| 国产成人av激情在线播放| 妹子高潮喷水视频| 亚洲性夜色夜夜综合| 最新在线观看一区二区三区| 不卡av一区二区三区| 国产又色又爽无遮挡免费看| 后天国语完整版免费观看| 脱女人内裤的视频| 国产在线一区二区三区精| 俄罗斯特黄特色一大片| 日韩中文字幕视频在线看片| 一级黄色大片毛片| av有码第一页| 色综合欧美亚洲国产小说| 老熟妇乱子伦视频在线观看| 亚洲精品美女久久久久99蜜臀| 久久久国产成人免费| 中文字幕高清在线视频| 日韩 欧美 亚洲 中文字幕| 免费看十八禁软件| 午夜免费鲁丝| 亚洲美女黄片视频| 国产精品 国内视频| 在线播放国产精品三级| 自线自在国产av| 久久中文字幕人妻熟女| 在线播放国产精品三级| 国产深夜福利视频在线观看| 黄频高清免费视频| 青青草视频在线视频观看| 午夜精品国产一区二区电影| 老司机影院毛片| 中文字幕色久视频| 午夜激情av网站| 欧美黄色片欧美黄色片| aaaaa片日本免费| 老司机深夜福利视频在线观看| 国产淫语在线视频| 国产精品偷伦视频观看了| 国精品久久久久久国模美| 国产在线一区二区三区精| 亚洲七黄色美女视频| 国产精品欧美亚洲77777| 婷婷成人精品国产| 91av网站免费观看| av欧美777| 免费看十八禁软件| 91九色精品人成在线观看| 免费少妇av软件| 18在线观看网站| 淫妇啪啪啪对白视频| 男人舔女人的私密视频| 天天操日日干夜夜撸| 黑人欧美特级aaaaaa片| 国产在线一区二区三区精| 久久99一区二区三区| 国产日韩一区二区三区精品不卡| 色在线成人网| 又大又爽又粗| 日日夜夜操网爽| 色在线成人网| 亚洲色图 男人天堂 中文字幕| 王馨瑶露胸无遮挡在线观看| 久久天堂一区二区三区四区| 女人精品久久久久毛片| 日韩中文字幕欧美一区二区| 国产亚洲精品久久久久5区| 国产高清videossex| 国产高清视频在线播放一区| 精品福利永久在线观看| 丝袜美足系列| 久久中文看片网| 三级毛片av免费| 一边摸一边抽搐一进一出视频| 欧美黄色片欧美黄色片| 成年女人毛片免费观看观看9 | 大型黄色视频在线免费观看| 一级片免费观看大全| 黄片播放在线免费| 日韩精品免费视频一区二区三区| 99久久国产精品久久久| 日韩一区二区三区影片| 欧美激情高清一区二区三区| xxxhd国产人妻xxx| 午夜福利乱码中文字幕| 两人在一起打扑克的视频| 手机成人av网站| 人妻一区二区av| 国产区一区二久久| 不卡一级毛片| a级毛片在线看网站| 国精品久久久久久国模美| 亚洲第一av免费看| 操出白浆在线播放| 女同久久另类99精品国产91| 国产伦人伦偷精品视频| 国精品久久久久久国模美| 巨乳人妻的诱惑在线观看| 久久精品亚洲av国产电影网| 老司机亚洲免费影院| 国产伦理片在线播放av一区| 日韩中文字幕视频在线看片| 男人舔女人的私密视频| 超色免费av| 黑人操中国人逼视频| 午夜福利一区二区在线看| 99精品久久久久人妻精品| 在线 av 中文字幕| videosex国产| 男女高潮啪啪啪动态图| 一个人免费看片子| 日本一区二区免费在线视频| 色老头精品视频在线观看| 91精品三级在线观看| 亚洲精品乱久久久久久| 三上悠亚av全集在线观看| 欧美精品一区二区大全| av国产精品久久久久影院| 99re在线观看精品视频| 亚洲国产av新网站| 另类精品久久| 久久亚洲真实| 一本综合久久免费| 一二三四社区在线视频社区8| 在线av久久热| 丝袜在线中文字幕| 又大又爽又粗| 9热在线视频观看99| 亚洲精品粉嫩美女一区| 久久午夜综合久久蜜桃| 黄网站色视频无遮挡免费观看| 国产又爽黄色视频| 视频区欧美日本亚洲| 日本一区二区免费在线视频| 热re99久久国产66热| 91字幕亚洲| 男女之事视频高清在线观看| 欧美中文综合在线视频| 女同久久另类99精品国产91| 亚洲成av片中文字幕在线观看| 国产欧美日韩综合在线一区二区| 亚洲午夜理论影院| 考比视频在线观看| 后天国语完整版免费观看| 精品国产国语对白av| 免费黄频网站在线观看国产| 亚洲黑人精品在线| 男女高潮啪啪啪动态图| 在线观看66精品国产| 欧美精品高潮呻吟av久久| 日韩大码丰满熟妇| 午夜激情av网站| 亚洲av电影在线进入| 精品第一国产精品| 少妇被粗大的猛进出69影院| 亚洲国产欧美在线一区| 精品人妻1区二区| 久久性视频一级片| 欧美日韩精品网址| 久久人妻熟女aⅴ| 丁香六月欧美| 男女床上黄色一级片免费看| 久久久久精品国产欧美久久久| 女人高潮潮喷娇喘18禁视频| 国产色视频综合| 老司机亚洲免费影院| 色94色欧美一区二区| 9色porny在线观看| 考比视频在线观看| 国产亚洲欧美精品永久| tube8黄色片| 99热国产这里只有精品6| 亚洲欧美精品综合一区二区三区| 日本a在线网址| 国产成人免费无遮挡视频| 国产一卡二卡三卡精品| 午夜久久久在线观看| 涩涩av久久男人的天堂| 精品国产一区二区久久| 狠狠精品人妻久久久久久综合| 日韩大码丰满熟妇| 侵犯人妻中文字幕一二三四区| 免费女性裸体啪啪无遮挡网站| 日日摸夜夜添夜夜添小说| 久久久欧美国产精品| 精品久久蜜臀av无| 国产精品麻豆人妻色哟哟久久| 午夜精品久久久久久毛片777| 又紧又爽又黄一区二区| cao死你这个sao货| 五月天丁香电影| 成人亚洲精品一区在线观看| 亚洲色图综合在线观看| 精品人妻熟女毛片av久久网站| 免费久久久久久久精品成人欧美视频| 巨乳人妻的诱惑在线观看| 精品一区二区三区四区五区乱码| 欧美日韩黄片免| 久久精品91无色码中文字幕| 亚洲午夜精品一区,二区,三区| 精品国产一区二区三区久久久樱花| 精品国内亚洲2022精品成人 | 亚洲久久久国产精品| 欧美日韩亚洲国产一区二区在线观看 | 老熟妇乱子伦视频在线观看| 建设人人有责人人尽责人人享有的| 国产亚洲精品久久久久5区| 黄色毛片三级朝国网站| 国产欧美日韩精品亚洲av| 国产精品 欧美亚洲| 丁香六月天网| e午夜精品久久久久久久| 999精品在线视频| 视频在线观看一区二区三区| 丝袜美足系列| 一个人免费在线观看的高清视频| 丝袜美足系列| 国产男女超爽视频在线观看| 久久精品国产综合久久久| 国产精品九九99| 久久人妻av系列| 国产免费av片在线观看野外av| 精品欧美一区二区三区在线| 欧美精品高潮呻吟av久久| aaaaa片日本免费| 欧美+亚洲+日韩+国产| 69精品国产乱码久久久| 妹子高潮喷水视频| 十八禁人妻一区二区| 国产精品1区2区在线观看. | 女警被强在线播放| 好男人电影高清在线观看| 久久九九热精品免费| 黄色成人免费大全| 色精品久久人妻99蜜桃| 女同久久另类99精品国产91| 桃红色精品国产亚洲av| 中文亚洲av片在线观看爽 | 欧美日韩亚洲国产一区二区在线观看 | 悠悠久久av| 国产精品成人在线| 一二三四在线观看免费中文在| 亚洲一码二码三码区别大吗| 亚洲色图av天堂| 99国产精品一区二区蜜桃av | 精品一品国产午夜福利视频| 亚洲欧洲日产国产| 久久精品熟女亚洲av麻豆精品| 日韩中文字幕视频在线看片| 免费少妇av软件| 久久久久视频综合| 国产真人三级小视频在线观看| 亚洲成国产人片在线观看| 精品亚洲成a人片在线观看| 国产精品秋霞免费鲁丝片| 999久久久国产精品视频| 精品熟女少妇八av免费久了| 成人国产一区最新在线观看| 午夜老司机福利片| 涩涩av久久男人的天堂| 中文字幕色久视频| 精品视频人人做人人爽| 男女高潮啪啪啪动态图| 18禁国产床啪视频网站| 色视频在线一区二区三区| 人人妻人人爽人人添夜夜欢视频| 中文字幕最新亚洲高清| 人妻一区二区av| 国产精品一区二区免费欧美| 亚洲熟妇熟女久久| 两个人看的免费小视频| 国产日韩欧美在线精品| 一级片'在线观看视频| 十八禁高潮呻吟视频| 久久这里只有精品19| 纯流量卡能插随身wifi吗| 国产精品偷伦视频观看了| 国产精品亚洲一级av第二区| 久久午夜亚洲精品久久| 久久久久久久大尺度免费视频| 如日韩欧美国产精品一区二区三区| 国产精品久久久久久精品古装| 中文欧美无线码| 欧美日韩黄片免| 男人操女人黄网站| 亚洲全国av大片| 热99re8久久精品国产| 亚洲av成人不卡在线观看播放网| 日韩大片免费观看网站| 日韩精品免费视频一区二区三区| 国产高清videossex| 宅男免费午夜| 国产99久久九九免费精品| 国产一区二区激情短视频| 成人18禁高潮啪啪吃奶动态图| 欧美黄色片欧美黄色片| 亚洲精品av麻豆狂野| 五月开心婷婷网| 午夜日韩欧美国产| 91av网站免费观看| 老司机在亚洲福利影院| 亚洲五月婷婷丁香| 国产成人欧美| 国产国语露脸激情在线看| 纯流量卡能插随身wifi吗| 99久久国产精品久久久| 嫩草影视91久久| 99久久精品国产亚洲精品| 丝袜喷水一区| 欧美国产精品一级二级三级| 午夜福利,免费看| 最近最新中文字幕大全电影3 | 成年人免费黄色播放视频| 久久青草综合色| 国产精品.久久久| 午夜日韩欧美国产| 男女免费视频国产| 9热在线视频观看99| 视频在线观看一区二区三区| 午夜两性在线视频| 日本av免费视频播放| 午夜视频精品福利| 久久久久久亚洲精品国产蜜桃av| 超碰97精品在线观看| 亚洲国产欧美一区二区综合| 久久久久久免费高清国产稀缺| 最黄视频免费看| 十八禁高潮呻吟视频| av又黄又爽大尺度在线免费看| av在线播放免费不卡| 久久久久精品国产欧美久久久| 亚洲av美国av| 丝袜美足系列| 99国产精品一区二区蜜桃av | 人妻 亚洲 视频| 最新美女视频免费是黄的| 亚洲精品美女久久av网站| 色播在线永久视频| 纵有疾风起免费观看全集完整版| 麻豆乱淫一区二区| 别揉我奶头~嗯~啊~动态视频| 久久精品国产a三级三级三级| 亚洲精品乱久久久久久| 热99久久久久精品小说推荐| 中文字幕av电影在线播放| 久久国产精品大桥未久av| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产欧美一区二区综合| 蜜桃国产av成人99| 桃花免费在线播放| 欧美变态另类bdsm刘玥| 看免费av毛片| 中文字幕人妻丝袜制服| 在线观看免费视频日本深夜| 国产极品粉嫩免费观看在线| 热99久久久久精品小说推荐| 搡老乐熟女国产| 国产一区二区三区综合在线观看| 一本久久精品| 亚洲一区二区三区欧美精品| 在线亚洲精品国产二区图片欧美| 一本色道久久久久久精品综合| av欧美777| 国产在线精品亚洲第一网站| 亚洲av日韩精品久久久久久密| 无限看片的www在线观看| 国产精品久久久久久人妻精品电影 | 久久精品人人爽人人爽视色| 成年人午夜在线观看视频| 午夜福利在线观看吧| 免费在线观看影片大全网站| 日本vs欧美在线观看视频| 国产日韩欧美在线精品| 国产一区二区 视频在线| 午夜福利免费观看在线| 在线十欧美十亚洲十日本专区| 日韩人妻精品一区2区三区| 欧美亚洲日本最大视频资源| aaaaa片日本免费| h视频一区二区三区| 久久人人爽av亚洲精品天堂| 99国产综合亚洲精品| 九色亚洲精品在线播放| 操出白浆在线播放| 精品国产一区二区三区四区第35| 国产一区二区 视频在线| 国产成人精品在线电影| 水蜜桃什么品种好| 亚洲人成电影免费在线| 日韩 欧美 亚洲 中文字幕| 9色porny在线观看| 俄罗斯特黄特色一大片| 在线观看免费视频日本深夜| 国产aⅴ精品一区二区三区波| 黄色视频不卡| 女性被躁到高潮视频| 黄片大片在线免费观看| 亚洲精品在线观看二区| 一级毛片电影观看| 久久精品人人爽人人爽视色| 欧美日韩精品网址| 99热国产这里只有精品6| 巨乳人妻的诱惑在线观看| 9色porny在线观看| 香蕉丝袜av| 国产亚洲av高清不卡| 国产极品粉嫩免费观看在线| 午夜福利视频在线观看免费| 国产高清国产精品国产三级| 老司机影院毛片| 黄色成人免费大全| 国产成人一区二区三区免费视频网站| 国产主播在线观看一区二区| 黄色毛片三级朝国网站| 高清黄色对白视频在线免费看| 热re99久久精品国产66热6| 国产成人免费观看mmmm| 久久中文看片网| 久久精品国产综合久久久| 麻豆av在线久日| 嫩草影视91久久| 夫妻午夜视频| 国产精品.久久久| 成年版毛片免费区| 99精品久久久久人妻精品| 日韩视频在线欧美| 精品第一国产精品| 免费观看人在逋| 精品国产亚洲在线| 国产单亲对白刺激| www.999成人在线观看| 老司机深夜福利视频在线观看| 久久性视频一级片| 亚洲欧美色中文字幕在线| 精品一区二区三区av网在线观看 | 国产av国产精品国产| 国产福利在线免费观看视频| 久久久久久人人人人人| 国产成人啪精品午夜网站| 国产一区二区三区综合在线观看| www日本在线高清视频| 亚洲色图 男人天堂 中文字幕| 亚洲欧美色中文字幕在线| 色视频在线一区二区三区| 日本a在线网址| 老司机靠b影院| 亚洲第一欧美日韩一区二区三区 | 成人影院久久| 国产有黄有色有爽视频| 国产无遮挡羞羞视频在线观看| 中文字幕色久视频| 99在线人妻在线中文字幕 | 少妇被粗大的猛进出69影院| 大片免费播放器 马上看| 久久精品91无色码中文字幕| 日本五十路高清| 精品第一国产精品| 精品免费久久久久久久清纯 | 亚洲中文av在线| 欧美日韩av久久| 久久久精品国产亚洲av高清涩受| 亚洲视频免费观看视频| 国产成+人综合+亚洲专区| 制服诱惑二区| 啦啦啦免费观看视频1| 窝窝影院91人妻| 啪啪无遮挡十八禁网站| 老熟女久久久| 制服人妻中文乱码| 99国产综合亚洲精品| 国产无遮挡羞羞视频在线观看| 精品午夜福利视频在线观看一区 | 9色porny在线观看| 亚洲精品中文字幕一二三四区 | 成人三级做爰电影| 成人免费观看视频高清| 一边摸一边做爽爽视频免费| 久久久久视频综合| 岛国毛片在线播放| 97人妻天天添夜夜摸| 日本撒尿小便嘘嘘汇集6| 免费在线观看完整版高清| 国产精品熟女久久久久浪| 免费在线观看黄色视频的| 在线永久观看黄色视频| 丰满少妇做爰视频| 后天国语完整版免费观看| h视频一区二区三区| 午夜91福利影院| tocl精华| 久热这里只有精品99| 精品久久久精品久久久| 制服诱惑二区| 国产高清激情床上av| 黑人操中国人逼视频| 日本av免费视频播放| 在线观看免费视频网站a站| 黑人操中国人逼视频| 国产av又大| 久久久精品免费免费高清| 亚洲国产欧美网| 如日韩欧美国产精品一区二区三区| avwww免费| 99九九在线精品视频| 国产av又大| 欧美变态另类bdsm刘玥| 我要看黄色一级片免费的| 99香蕉大伊视频| 1024香蕉在线观看| 色婷婷久久久亚洲欧美| 亚洲国产av新网站| 天堂动漫精品| 一级a爱视频在线免费观看| 欧美黄色淫秽网站| 啦啦啦 在线观看视频| 国产亚洲av高清不卡| 成人三级做爰电影| 大陆偷拍与自拍|