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

    調(diào)控萼脊蘭類(lèi)黃酮生物合成的MYB轉(zhuǎn)錄因子挖掘分析

    2021-05-19 15:57:38蔣素華牛蘇燕袁秀云梁芳王喜蒙崔波
    熱帶作物學(xué)報(bào) 2021年3期

    蔣素華 牛蘇燕 袁秀云 梁芳 王喜蒙 崔波

    摘 ?要:本研究通過(guò)生物信息學(xué)方法,挖掘萼脊蘭中可能調(diào)控類(lèi)黃酮生物合成的MYB轉(zhuǎn)錄因子。從萼脊蘭轉(zhuǎn)錄組數(shù)據(jù)庫(kù)中篩選MYB轉(zhuǎn)錄因子,獲得MYB的完整ORF,并對(duì)其蛋白理化性質(zhì)、基序、功能注釋、系統(tǒng)進(jìn)化、表達(dá)特性等進(jìn)行分析,并預(yù)測(cè)其功能。結(jié)果表明:在萼脊蘭轉(zhuǎn)錄組水平上共鑒定出27個(gè)具有完整閱讀框的MYB轉(zhuǎn)錄因子基因,萼脊蘭MYB轉(zhuǎn)錄因子TRINITY_DN38485_c0_g4和TRINITY_DN38485_c0_g1基因均在花朵和葉片中表達(dá)上調(diào),再根據(jù)與蝴蝶蘭的系統(tǒng)進(jìn)化關(guān)系與功能分析,預(yù)測(cè)這2個(gè)轉(zhuǎn)錄因子可能通過(guò)黃酮途徑完成類(lèi)黃酮生物合成。因此,MYB轉(zhuǎn)錄因子的研究對(duì)從分子水平上研究和調(diào)節(jié)類(lèi)黃酮的合成具有重要意義,轉(zhuǎn)錄因子的應(yīng)用是類(lèi)黃酮生物合成基因工程中的新途徑。

    關(guān)鍵詞:萼脊蘭;MYB轉(zhuǎn)錄因子;次生代謝

    中圖分類(lèi)號(hào):Q946.8 ? ? ?文獻(xiàn)標(biāo)識(shí)碼:A

    Mining Analysis of MYB Transcription Factors Related to Flavonoid Biosynthesis of Sedirea japonica

    JIANG Suhua1, NIU Suyan1, YUAN Xiuyun1, LIANG Fang1, WANG Ximeng2, CUI Bo1*

    1. Biological Engineering Research Center, Zhengzhou Normal University, Zhengzhou, Henan 450044, China; 2. College of Life Sciences, Zhengzhou University, Zhengzhou, Henan 450002, China

    Abstract: MYB transcription factors that may regulate the biosynthesis of the flavonoids in Sedirea japonica were detected using bioinformatics methods. The complete ORF of MYB was obtained, and the physicochemical properties, motif, functional annotation, phylogenetic evolution, and expression characteristics of its proteins were analyzed, and its functions were also predicted. A total of 27 MYB transcription factor genes with complete reading frames were identified at the transcriptome level of S. japonica. MYB transcription factors TRINITY_DN38485_c0_g4 and TRINITY_ DN38485_c0_g1 were up-regulated in flowers and leaves. Then according to the phylogenetic relationship with phalaenopsis and functional analysis,it is predicted that the two transcription factors were involved in flavonoid biosynthesis through the flavonoid pathway. Therefore, the study of MYB transcription factor is of great significance to study and regulate the synthesis of flavonoids at the molecular level. The application of transcription factors is a new approach in the genetic engineering of flavonoid biosynthesis, the application of transcription factors is a new approach in the genetic engineering of flavonoid biosynthesis.

    Keywords: S. japonica; MYB transcription factors; secondary metabolism

    DOI: 10.3969/j.issn.1000-2561.2021.03.003

    轉(zhuǎn)錄因子(transcription factor)是一群能與基因5? 端上游特定序列專(zhuān)一性結(jié)合,從而保證目的基因以特定強(qiáng)度在特定時(shí)間與空間表達(dá)的蛋白質(zhì)分子。根據(jù)轉(zhuǎn)錄因子的作用特點(diǎn)可分為兩類(lèi):第一類(lèi)為普遍轉(zhuǎn)錄因子;第二類(lèi)轉(zhuǎn)錄因子為組織細(xì)胞特異性轉(zhuǎn)錄因子[1-2]。人們?cè)谥参镏幸呀?jīng)找到了很多調(diào)控黃酮類(lèi)次生代謝的特異性轉(zhuǎn)錄因子基因,主要包括編碼MYB、MYC、bZIP、WD40蛋白和鋅指蛋白等基因家族[3-4]。MYB基因家族廣泛存在于植物中,是植物中最大的轉(zhuǎn)錄因子基因家族之一。在類(lèi)黃酮生物合成過(guò)程中,MYB轉(zhuǎn)錄因子扮演著重要角色,可調(diào)控與類(lèi)黃酮合成相關(guān)的酶基因的表達(dá),從而有效地調(diào)控類(lèi)黃酮物質(zhì)的生物合成。根據(jù)MYB的DNA結(jié)合結(jié)構(gòu)域中不完全重復(fù)序列的數(shù)目將其分為4個(gè)亞類(lèi):1R-MYB、R2R3-MYB、3R-MYB和4R-MYB。其中,R2R3-MYB是MYB轉(zhuǎn)錄因子家族中最大的亞類(lèi),其在類(lèi)黃酮合成中的作用已被廣泛研究[5-6]。目前擬南芥、水稻和小蘭嶼蝴蝶蘭的MYB轉(zhuǎn)錄因子的數(shù)據(jù)庫(kù)已經(jīng)建成,為其他蘭科植物的研究打下了良好的基礎(chǔ)。

    類(lèi)黃酮在植物中廣泛存在,并廣泛參與植物中花色素的形成、花香的調(diào)節(jié)、抗病、抗紫外線損傷等多種生物學(xué)過(guò)程[7]。李崇暉等[8]對(duì)不同顏色的蝴蝶石斛蘭品種花朵中的類(lèi)黃酮組成進(jìn)行分析,推測(cè)了蝴蝶石斛蘭花中的類(lèi)黃酮代謝途徑。鐘淮欽等[9]克隆了類(lèi)黃酮-3?5?-羥基化酶基因(F3?5?H),并對(duì)該基因進(jìn)行了表達(dá)分析。張加強(qiáng)等[10]鑒定了蝴蝶蘭花色苷合成的相關(guān)基因,從轉(zhuǎn)錄組水平揭示蝴蝶蘭花色苷生物合成的分子調(diào)控機(jī)制,結(jié)果顯示類(lèi)黃酮生物合成中差異基因?yàn)?個(gè),綜上所述,關(guān)于類(lèi)黃酮代謝途徑及相關(guān)基因研究的較多,而類(lèi)黃酮合成相關(guān)的MYB轉(zhuǎn)錄因子在植物優(yōu)良性狀形成及微量次生代謝物的合成中的調(diào)控研究還較少。

    萼脊蘭(Sedirea japonica)為蘭科(Orchidaceae)萼脊蘭屬(Sedirea Garay & H. R. Sweet)多年生植物,總狀花序從葉腋中發(fā)出,花具橘子香氣,花多且花期長(zhǎng),因此研究類(lèi)黃酮對(duì)萼脊蘭花香的調(diào)節(jié)具有重要意義。目前還未見(jiàn)萼脊蘭MYB轉(zhuǎn)錄因子家族基因的研究報(bào)道,本研究基于測(cè)序的萼脊蘭轉(zhuǎn)錄組數(shù)據(jù)庫(kù)信息,篩選并鑒定出萼脊蘭所有MYB轉(zhuǎn)錄因子,并對(duì)其進(jìn)行理化性質(zhì)分析、蛋白基序分析、系統(tǒng)發(fā)育樹(shù)構(gòu)建以及基因表達(dá)分析,這些結(jié)果為進(jìn)一步探索萼脊蘭MYB基因提供理論資料。

    1 ?材料與方法

    1.1 ?材料

    萼脊蘭采自鄭州師范學(xué)院觀賞、藥用蘭花河南省重點(diǎn)實(shí)驗(yàn)室蘭花種植基地。萼脊蘭MYB轉(zhuǎn)錄因子是本實(shí)驗(yàn)室構(gòu)建的轉(zhuǎn)錄組數(shù)據(jù)庫(kù)。蝴蝶蘭MYB轉(zhuǎn)錄因子家族氨基酸序列下載于蝴蝶蘭信息資源(PlantTFDB)數(shù)據(jù)庫(kù)(http://planttfdb. cbi.pku.edu.cn/quick_search_result.php)。

    1.2 ?萼脊蘭MYB轉(zhuǎn)錄本篩選

    高通量測(cè)序由上海生工生物科技有限公司完成,拼接組裝后獲得萼脊蘭轉(zhuǎn)錄組,采用轉(zhuǎn)錄組基因ID搜索獲取氨基酸序列,使用ORF Finder在線軟件分析預(yù)測(cè)開(kāi)放閱讀框(ORF),獲得具有全長(zhǎng)氨基酸序列的萼脊蘭MYB轉(zhuǎn)錄因子家族蛋白基因序列。

    1.3 ?蛋白基序分析、結(jié)構(gòu)預(yù)測(cè)、GO注釋

    運(yùn)用The MEME Suite 4.12.0程序分析萼脊蘭MYB家族蛋白的基序;應(yīng)用結(jié)構(gòu)域在線分析軟件PROSITE和SMART對(duì)保守結(jié)構(gòu)域進(jìn)行分析。選擇具有完整ORF的萼脊蘭MYB轉(zhuǎn)錄因子氨基酸序列,利用在線工具ProtParam對(duì)蛋白進(jìn)行理化性質(zhì)分析;利用在線軟件SOPMA對(duì)蛋白二級(jí)結(jié)構(gòu)進(jìn)行預(yù)測(cè)分析。使用BLAST2GO軟件對(duì)篩選獲得的萼脊蘭MYB轉(zhuǎn)錄因子進(jìn)行功能注釋分析。

    1.4 ?系統(tǒng)進(jìn)化樹(shù)構(gòu)建及功能分析

    蝴蝶蘭MYB轉(zhuǎn)錄因子的蛋白序列從PlantTFDB數(shù)據(jù)庫(kù)中下載,完整的氨基酸序列用于進(jìn)化分析。利用Clustal X程序和MEGA4對(duì)萼脊蘭MYB家族蛋白的氨基酸序列和下載的蝴蝶蘭所有MYB轉(zhuǎn)錄因子的蛋白序列進(jìn)行比對(duì)分析,將比對(duì)結(jié)果采用鄰接法構(gòu)建系統(tǒng)發(fā)育樹(shù)。

    通過(guò)萼脊蘭MYB轉(zhuǎn)錄因子與蝴蝶蘭所有MYB轉(zhuǎn)錄因子構(gòu)建的系統(tǒng)進(jìn)化樹(shù),篩選獲得與各個(gè)萼脊蘭MYB轉(zhuǎn)錄因子相似性最高的蝴蝶蘭MYB轉(zhuǎn)錄因子,然后通過(guò)分析預(yù)測(cè)MYB轉(zhuǎn)錄因子的分子功能。

    1.5 ?萼脊蘭MYB轉(zhuǎn)錄因子表達(dá)分析

    基于轉(zhuǎn)錄組數(shù)據(jù),獲取萼脊蘭MYB轉(zhuǎn)錄因子在葉、花蕾和花朵的表達(dá)量,采用HemI 1.0.3.7軟件對(duì)MYB轉(zhuǎn)錄因子在各組織表達(dá)數(shù)據(jù)進(jìn)行分層聚類(lèi)分析。

    2 ?結(jié)果與分析

    2.1 ?萼脊蘭MYB家族成員的獲得

    從萼脊蘭轉(zhuǎn)錄組數(shù)據(jù)庫(kù)中共篩選出55條MYB基因序列,片段大小為205~3982 bp。通過(guò)ORF預(yù)測(cè)全長(zhǎng)氨基酸序列,發(fā)現(xiàn)27條序列具有完整的ORF,片段最短的為499 bp(TRINITY_ DN34588_c1_g2_i3),編碼的154個(gè)氨基酸;片段最長(zhǎng)的為2461 bp(TRINITY_DN29682_ c0_g1_i1),編碼535個(gè)氨基酸。

    2.2 ?萼脊蘭MYB轉(zhuǎn)錄因子基序分析

    利用MEME在線軟件分析了27個(gè)萼脊蘭MYB轉(zhuǎn)錄因子家族的10個(gè)基序,基序1由50個(gè)氨基酸殘基組成,其中包含3個(gè)高度保守的色氨

    酸(W)殘基,每個(gè)色氨酸之間間隔19個(gè)氨基酸,屬于R2MYB結(jié)構(gòu)域(圖1C);基序2由41個(gè)氨基酸殘基組成,其中第1個(gè)色氨酸被苯丙氨酸(F)所代替,每個(gè)色氨酸之間間隔18個(gè)氨基酸,其他2個(gè)色氨酸殘基高度保守,屬于R3 MYB結(jié)構(gòu)域(圖1D);萼脊蘭21個(gè)轉(zhuǎn)錄因子包含基序1和基序2,這21條序列均包含R2R3結(jié)構(gòu)域,該結(jié)構(gòu)域是高度保守的,在該區(qū)域2(甘氨酸G)、4(色氨酸W)、8~9(谷氨酸E、天冬氨酸D)、12(亮氨酸L)、20(G)和38~44(賴(lài)氨酸K、絲氨酸S、半胱氨酸C、精氨酸R、亮氨酸L、精氨酸R、色氨酸W)等30個(gè)位點(diǎn)保守不變(圖1)。進(jìn)一步分析保守的motif序列特征發(fā)現(xiàn),27個(gè)MYB轉(zhuǎn)錄因子均具有典型的HTH結(jié)構(gòu)域。利用分析軟件SMART分析發(fā)現(xiàn),27個(gè)蛋白均具有DNA結(jié)合區(qū)-Myb結(jié)構(gòu)域。

    2.3 ?萼脊蘭MYB家族蛋白理化性質(zhì)分析

    萼脊蘭MYB家族蛋白的理化性質(zhì)分析結(jié)果表明(表1),27個(gè)MYB家族蛋白的平均相對(duì)分子質(zhì)量為32 651.65,TRINITY_DN67618_c0_g1_i1編碼蛋白的分子質(zhì)量最小,為17 527.05;TRINITY_ DN29682_c0_g1_i1編碼蛋白的分子質(zhì)量最大,為58 852.81。等電點(diǎn)(pI)平均值為7.37,TRINITY_ DN43325_c1_g2_i2編碼蛋白的pI最小,為4.67;TRINITY_DN30733_c0_g1_i1編碼蛋白的pI最大,為9.93。其中有13個(gè)MYB家族蛋白pI小于7(偏酸性),14個(gè)蛋白pI大于7(偏堿性),說(shuō)明大多數(shù)萼脊蘭MYB家族蛋白表現(xiàn)為偏堿性。MYB家族蛋白的平均疏水性值均小于0,說(shuō)明27個(gè)萼脊蘭MYB家族蛋白均屬于疏水蛋白。MYB家族蛋白二級(jí)結(jié)構(gòu)特征預(yù)測(cè)結(jié)果顯示,27個(gè)萼脊蘭MYB蛋白均具有α螺旋、無(wú)規(guī)卷曲、延伸鏈和β轉(zhuǎn)角4種構(gòu)型,主要由α螺旋和無(wú)規(guī)卷曲構(gòu)成(所占比平均值分別為32.28%和54.01%),延伸鏈和β轉(zhuǎn)角所占比例較?。ㄋ急绕骄捣謩e為8.45%和5.25%),TRINITY_DN43325_c1_g2_i2和TRI NITY_DN38485_c0_g4_i2基因編碼蛋白α螺旋所占比例最大,其余25個(gè)MYB蛋白均表現(xiàn)為無(wú)規(guī)卷曲所占比例最大。

    2.4 ?萼脊蘭MYB轉(zhuǎn)錄因子的GO注釋分析

    GO注釋分析表明(圖2),萼脊蘭27條MYB轉(zhuǎn)錄因子序列中分別有22、27、27條被注釋到生物過(guò)程、分子功能和細(xì)胞組分,所有的序列被進(jìn)一步富集為27個(gè)功能類(lèi)別,其中DNA結(jié)合功能(27個(gè))、細(xì)胞核功能(27個(gè))、轉(zhuǎn)錄作用(24個(gè))、轉(zhuǎn)錄調(diào)節(jié)(17個(gè))、轉(zhuǎn)錄因子活性(13個(gè))、防御響應(yīng)(9)、茉莉酸響應(yīng)(6個(gè))、花器官發(fā)育(6個(gè))等功能組中包含的序列較多,水楊酸響應(yīng)(5個(gè))、細(xì)胞分化(4個(gè))、鹽脅迫響應(yīng)(4個(gè))、乙烯響應(yīng)(4個(gè))、脫落酸響應(yīng)(4個(gè))、赤霉素調(diào)控(3個(gè))、鋅離子結(jié)合(2個(gè))、細(xì)胞壁生物發(fā)生(2個(gè))、缺水響應(yīng)(2個(gè))、生長(zhǎng)素響應(yīng)(2個(gè))、脫落酸激活信號(hào)通路調(diào)控(2個(gè))、幾丁質(zhì)響應(yīng)(2個(gè))、對(duì)鎘離子的響應(yīng)(2個(gè))、苯丙醇代謝過(guò)程的調(diào)控(2個(gè))、茉莉酸介導(dǎo)的信號(hào)通路的調(diào)控(2個(gè))、鋅離子結(jié)合(2個(gè))、水楊酸介導(dǎo)的信號(hào)通路調(diào)控(1個(gè))、花藥壁絨氈層發(fā)育(1個(gè))、花粉精細(xì)胞分化(1個(gè))等功能組中包含的序列較少。其中生物過(guò)程主要富集在細(xì)胞合成、生物響應(yīng)和調(diào)控上,分子功能集中富集在結(jié)合和轉(zhuǎn)錄因子活性上,細(xì)胞組分主要富集在細(xì)胞核部分。綜上所述,萼脊蘭MYB轉(zhuǎn)錄因子通過(guò)結(jié)合到DNA區(qū)域發(fā)揮相應(yīng)的調(diào)控功能,主要涉及植物的細(xì)胞過(guò)程、發(fā)育過(guò)程和物質(zhì)代謝等過(guò)程。

    2.5 ?萼脊蘭MYB蛋白進(jìn)化分析

    為了研究萼脊蘭MYB蛋白的進(jìn)化關(guān)系并預(yù)測(cè)其功能,從PlantTFDB數(shù)據(jù)庫(kù)中下載蝴蝶蘭已有的11個(gè)MYB家族蛋白氨基酸序列,采用軟件MEGA4構(gòu)建萼脊蘭MYB家族蛋白和蝴蝶蘭MYB家族蛋白進(jìn)化樹(shù)。由圖3可知,4個(gè)萼脊蘭MYB家族蛋白分別與蝴蝶蘭的MYB蛋白聚集在一起,5個(gè)蝴蝶蘭的MYB蛋白聚在一起,且比較集中,說(shuō)明萼脊蘭MYB家族蛋白與蝴蝶蘭MYB蛋白的保守性不是很高。

    2.6 ?萼脊蘭次生代謝相關(guān)MYB轉(zhuǎn)錄因子功能分析

    通過(guò)構(gòu)建萼脊蘭MYB轉(zhuǎn)錄因子與蝴蝶蘭所有MYB轉(zhuǎn)錄因子系統(tǒng)進(jìn)化樹(shù),篩選同源性較高的蝴蝶蘭MYB轉(zhuǎn)錄因子,然后分析MYB轉(zhuǎn)錄因子的分子功能。由表2可知,MYB轉(zhuǎn)錄因子基因TRINITY_DN38485_c0_g4_i2編碼蛋白同源性較高的蝴蝶蘭轉(zhuǎn)錄因子PEQU 10866,編碼的蛋白為MYB8,推測(cè)其具有調(diào)控花色素沉積和參與調(diào)節(jié)類(lèi)黃酮的生物合成的功能;TRINITY DN38485_ c0_g1_i1編碼蛋白同源性較高的蝴蝶蘭轉(zhuǎn)錄因子PEQU 23041,編碼的蛋白為MYB6,推測(cè)其可能參與調(diào)節(jié)類(lèi)黃酮的生物合成。TRINITY DN40180_ c3_g1_i1編碼蛋白同源性較高的蝴蝶蘭轉(zhuǎn)錄因子PEQU 05119,編碼的蛋白為MYB10,推測(cè)其可能具有協(xié)調(diào)植物生長(zhǎng)發(fā)育的功能;TRINITY DN43325 c1 g2 i2編碼蛋白同源性較高的蝴蝶蘭轉(zhuǎn)錄因子PEQU 30611,編碼的蛋白MYB3,推測(cè)其可能具有抑制C4H基因表達(dá)的功能。

    2.7 ?萼脊蘭MYB轉(zhuǎn)錄因子基因表達(dá)分析

    基于萼脊蘭轉(zhuǎn)錄組數(shù)據(jù),獲得萼脊蘭MYB轉(zhuǎn)錄因子在葉、花蕾、花朵各個(gè)組織的基因表達(dá)量,并繪制熱圖(圖4)。從圖4可見(jiàn),MYB轉(zhuǎn)錄因子在組織表達(dá)模式上存在差異,大多數(shù)成員在葉、花蕾、花朵中均有表達(dá),表達(dá)模式呈現(xiàn)多樣化。萼脊蘭MYB轉(zhuǎn)錄因子中10個(gè)MYB轉(zhuǎn)錄因子基因在花朵中的表達(dá)上調(diào),5個(gè)MYB轉(zhuǎn)錄因子基因在花蕾中表達(dá)上調(diào),12個(gè)MYB轉(zhuǎn)錄因子基因在葉片中的表達(dá)上調(diào),而萼脊蘭類(lèi)黃酮的代謝產(chǎn)物主要在花朵和葉片中積累。萼脊蘭MYB 轉(zhuǎn)錄因子TRINITY_ DN38485_c0_g4和TRINITY_DN38485_c0_g1基因均在花朵和葉片中表達(dá)量上調(diào),根據(jù)表達(dá)模式分析:萼脊蘭MYB轉(zhuǎn)錄因子TRINITY_DN38485_ c0_g4和TRINITY_DN38485_c0_g1可能調(diào)節(jié)類(lèi)黃酮的生物合成,這也與功能預(yù)測(cè)一致。

    3 ?討論

    通過(guò)萼脊蘭轉(zhuǎn)錄組數(shù)據(jù)分析,共篩選出55條MYB基因序列,片段大小為205~3982 bp。通過(guò)ORF預(yù)測(cè)全長(zhǎng)氨基酸序列,發(fā)現(xiàn)27條序列具有完整的ORF。GO注釋分析表明,萼脊蘭27條MYB轉(zhuǎn)錄因子序列中分別有22、27、27條被注釋到生物過(guò)程、分子功能和細(xì)胞組分。系統(tǒng)進(jìn)化樹(shù)分析發(fā)現(xiàn),萼脊蘭MYB家族蛋白與蝴蝶蘭MYB蛋白的保守性不是很高。

    類(lèi)黃酮廣泛存在于高等植物果實(shí)、花瓣和葉片等器官中,萼脊蘭MYB 轉(zhuǎn)錄因子在組織表達(dá)模式上存在差異,在葉片、花蕾和花瓣中均有表達(dá),表達(dá)模式呈現(xiàn)多樣化。類(lèi)黃酮的合成主要包括7個(gè)分支:花青素(anthocyanins)途徑、原花青素(proanthocyanidins)途徑、黃酮(flavones)途徑、黃酮醇(flavonols)途徑、異黃酮(isoflavone)和異類(lèi)黃酮(isoflavonoids)途徑、鞣酐(phlobaphenes)途徑以及橙酮(aurones)途徑[11-12]。預(yù)測(cè)萼脊蘭MYB轉(zhuǎn)錄因子TRINITY_DN38485_c0_g4和TRINITY_DN38485_c0_g1可能通過(guò)黃酮途徑調(diào)節(jié)類(lèi)黃酮的生物合成。

    除了對(duì)MYB蛋白功能的發(fā)掘,一些研究者還對(duì)其結(jié)構(gòu)進(jìn)行了深入研究,發(fā)現(xiàn)調(diào)控類(lèi)黃酮合成的MYB轉(zhuǎn)錄因子大部分都是R2R3MYB轉(zhuǎn)錄因子[13]。蔣卉等[14]研究了蝴蝶蘭響應(yīng)病毒關(guān)鍵MYB家族基因,發(fā)現(xiàn)152個(gè)MYB家族基因中R2R3型有138個(gè)。本研究通過(guò)對(duì)萼脊蘭蛋白基序分析發(fā)現(xiàn),27個(gè)MYB家族基因種21個(gè)基因?qū)儆赗2R3-MYB類(lèi)基因,大部分也都是R2R3型轉(zhuǎn)錄因子,兩者研究結(jié)果一致。R2R3-MYB類(lèi)4個(gè)轉(zhuǎn)錄因子AtMYB3、AtMYB4、AtMYB7和AtMYB32均為類(lèi)黃酮的生物合成途徑的抑制因子[15-17]。日本龍膽GtMYBP3和GtMYBP4基因在擬南芥中的異源表達(dá)激活了內(nèi)源黃酮醇生物合成基因的表達(dá),導(dǎo)致幼苗中黃酮醇積累增加[18]。在擬南芥中,較早發(fā)現(xiàn)并證實(shí)的調(diào)控類(lèi)黃酮合成的MYB轉(zhuǎn)錄因子基因?yàn)锳tMYB12[19],同時(shí),在其他物種中也發(fā)現(xiàn)了調(diào)控類(lèi)黃酮合成MYB12的同源基因,如番茄SlMYB12[20]、龍膽花GtMYBP3和GtMYBP4[18]。萼脊蘭RINITY_DN38485_c0_g4和TRINITY_DN38485_c0_g1同源的蝴蝶蘭MYB 蛋白有MYB7,TRINITY_DN38485_c0_g1同源的蝴蝶蘭MYB蛋白為MYB6和MYB4,同日本龍膽、龍膽、番茄和擬南芥的研究結(jié)果一致,說(shuō)明萼脊蘭的這3個(gè)預(yù)測(cè)的基因與類(lèi)黃酮的生物合成密切相關(guān)。王雪霽等[21]研究發(fā)現(xiàn)39個(gè)PeMYBs基因在4種器官(根、莖、葉、花)中均表達(dá),48個(gè)基因在不同器官中有特異性不表達(dá)現(xiàn)象,一些基因呈現(xiàn)器官特異性表達(dá)特點(diǎn),與該研究萼脊蘭MYB轉(zhuǎn)錄因子在組織表達(dá)模式上存在差異,大多數(shù)成員在葉、花蕾、花朵中均有表達(dá),表達(dá)模式呈現(xiàn)多樣化的結(jié)果一致。

    參考文獻(xiàn)

    [1] Liu J Y, Osbourn A, Ma P D. MYB transcription factors as regulators of phenylpropanoid metabolism in plants[J]. Molecular Plant, 2015, 8(5): 689-708.

    [2] Xu F, Ning Y J, Zhang W W, et al. An R2R3-MYB transcription factor as a negative regulator of the flavonoid biosynthesis pathway in Ginkgo biloba[J]. Functional & Integrative Genomics, 2014, 14(1): 177-189.

    [3] Wang N, Xu H F, Jiang S H, et al. MYB12 and MYB22 play essential roles in proanthocyanidin and flavonol synthesis in red-fleshed apple (Malus sieversii f. niedzwetzkyana)[J]. The Plant Journal, 2017, 90(2): 276-292.

    [4] Liu C Y, Long J M, Zhu K J, et al. Characterization of a Citrus R2R3-MYB transcription factor that regulates the flavonol and hydroxycinnamic acid biosynthesis[J]. Scientific Reports, 2016, 6: 25352.

    [5] Huang W J, Khaldun A B M, Chen J J, et al. A R2R3-MYB transcription factor regulates the flavonol biosynthetic pathway in a traditional chinese medicinal plant, Epimedium sagittatum[J]. Frontiers in Plant Science, 2016, 7: 1089.

    [6] 宋 ?楊, 劉紅弟, 張紅軍, 等. 越橘果實(shí)轉(zhuǎn)錄組及R2R3-MYB轉(zhuǎn)錄因子分析[J]. 園藝學(xué)報(bào), 2015, 42(12): 2383-2394.

    [7] Winkel-shirley B. Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology[J]. Plant Physiology, 2001, 126(2): 485-493.

    [8] 李崇暉, 任 ?羽, 黃素榮, 等. 蝴蝶石斛蘭花色表型及類(lèi)黃酮成分分析[J]. 園藝學(xué)報(bào), 2013, 40(1): 107-116.

    [9] 鐘淮欽, 黃敏玲, 吳建設(shè), 等. 蝴蝶蘭花青素苷合成途徑關(guān)鍵酶基因的克隆與表達(dá)[J]. 福建農(nóng)業(yè)學(xué)報(bào), 2013, 28(9): 854-858.

    [10] 張加強(qiáng), 史小華, 劉慧春, 等. 基于轉(zhuǎn)錄組學(xué)的不同色系蝴蝶蘭花色苷差異積累分析[J]. 分子植物育種, 2018, 16(14): 4530-4542.

    [11] Lai Y S, Li H X, Yamagishi M. A review of target gene specificity of flavonoid R2R3-MYB transcription factors and a discussion of factors contributing to the target gene selectivity[J]. Frontiers in Biology, 2013, 8(6): 577-598.

    [12] 伍小方, 高國(guó)應(yīng), 左 ?倩, 等. FtMYB1轉(zhuǎn)錄因子調(diào)控苦蕎毛狀根黃酮醇合成的機(jī)理研究[J]. 植物遺傳資源學(xué)報(bào), 2020, 21(5): 1270-1278.

    [13] 邢 ?文, 金曉玲. 調(diào)控植物類(lèi)黃酮生物合成的MYB轉(zhuǎn)錄因子研究進(jìn)展[J]. 分子植物育種, 2015, 13(3): 689-696.

    [14] 蔣 ?卉, 張 ?晶, 袁 ?欣, 等. 蝴蝶蘭MYB家族基因鑒定及應(yīng)對(duì)病毒的表達(dá)響應(yīng)[J]. 分子植物育種, 2018, 16(21): 6921-6930.

    [15] Dubos C, Stracke R, Grotewold E, ?et al. MYB transcription factors in Arabidopsis[J]. Trends in Plant Science, 2010, 15(10): 573-581.

    [16] Jin H, Cominelli E, Bailey P, et al. Transcriptional repression by AtMYB4 controls production of UV-protecting sunscreens in Arabidopsis[J]. The EMBO Journal, 2000, 19(22): 6150-6161.

    [17] Fornalé S, Lopez E, Salazar-Henao J E, et al. AtMYB7, a new player in the regulation of UV-sunscreens in Arabidopsis thaliana[J]. Plant & Cell Physiology, 2014, 55(3): 507-516.

    [18] Nakatsuka T, Saito M, Yamada E, et al. Isolation and characterization of GtMYBP3 and GtMYBP4, orthologues of R2R3-MYB transcription factors that regulate early flavonoid biosynthesis, in gentian flowers[J]. Journal of Experimental Botany, 2012, 63(18): 6505-6517.

    [19] Mehrtens F, Kranz H, Bednarek P, et al. The Arabidopsis transcription factor MYB12 is a flavonol-specific regulator of phenylpropanoid biosynthesis[J]. Plant Physiology, 2005, 138(2): 1083-1096.

    [20] Ballester A R, Molthoff J, De vos R, et al. Biochemical and molecular analysis of pink tomatoes: deregulated ex-pression of the gene encoding transcription factor SlMYB12 leads to pink tomato fruit color[J]. Plant Physiology, 2010, 152(1): 71-84.

    [21] 王雪霽, 梁立雄, 李潞濱, 等. 小蘭嶼蝴蝶蘭R2R3-MYB轉(zhuǎn)錄因子分析[J]. 林業(yè)科學(xué)研究, 2018, 31(3): 104-113.

    責(zé)任編輯:黃東杰

    收稿日期 ?2020-05-05;修回日期 ?2020-05-26

    基金項(xiàng)目 ?河南省科技廳科技發(fā)展計(jì)劃項(xiàng)目(No. 172102110227);河南省科技攻關(guān)項(xiàng)目(No. 182102110357);河南省重大科技專(zhuān)項(xiàng)(No. 701290)。

    作者簡(jiǎn)介 ?蔣素華(1983—),女,碩士,副教授,研究方向:花卉分子生物學(xué)。*通信作者(Corresponding author):崔 ?波(CUI Bo),E-mail:laocuibo@163.com。

    白带黄色成豆腐渣| 乱人视频在线观看| 国产精品三级大全| 一个人观看的视频www高清免费观看| 国产高清有码在线观看视频| 亚洲婷婷狠狠爱综合网| 欧美日韩乱码在线| 一区二区三区免费毛片| 午夜激情福利司机影院| 女的被弄到高潮叫床怎么办| 成人综合一区亚洲| 国产精品野战在线观看| 国产高清三级在线| 国产精品久久久久久久久免| 国产淫片久久久久久久久| 欧美xxxx性猛交bbbb| 国产视频内射| 久久精品久久久久久久性| 99热这里只有是精品在线观看| 久久99热6这里只有精品| 特大巨黑吊av在线直播| 中文字幕免费在线视频6| 中文字幕人妻熟人妻熟丝袜美| 最近最新中文字幕大全电影3| 亚洲精华国产精华液的使用体验 | 国产精品av视频在线免费观看| 淫秽高清视频在线观看| 成人午夜高清在线视频| 午夜福利视频1000在线观看| 热99re8久久精品国产| 久久久成人免费电影| 国产不卡一卡二| 婷婷六月久久综合丁香| 亚洲18禁久久av| av在线天堂中文字幕| 午夜老司机福利剧场| 又爽又黄无遮挡网站| 亚洲欧美日韩高清在线视频| 国产成人午夜福利电影在线观看| 欧美不卡视频在线免费观看| 大型黄色视频在线免费观看| 国产精品av视频在线免费观看| 久久久久久国产a免费观看| 欧美性猛交╳xxx乱大交人| 成年女人永久免费观看视频| 免费看美女性在线毛片视频| 六月丁香七月| 国产人妻一区二区三区在| 少妇熟女欧美另类| 国产精品一二三区在线看| 九九爱精品视频在线观看| 亚洲在线自拍视频| 国产成人午夜福利电影在线观看| 岛国毛片在线播放| av.在线天堂| 日韩一本色道免费dvd| 床上黄色一级片| 亚洲国产精品成人久久小说 | 久久99热这里只有精品18| 欧美成人a在线观看| 亚洲美女搞黄在线观看| 在线a可以看的网站| 国产成人91sexporn| 亚洲七黄色美女视频| 午夜爱爱视频在线播放| 国内久久婷婷六月综合欲色啪| 欧美成人a在线观看| 久久这里只有精品中国| 国产伦在线观看视频一区| 一个人看视频在线观看www免费| 少妇熟女aⅴ在线视频| 成年版毛片免费区| 寂寞人妻少妇视频99o| 日本与韩国留学比较| 99热网站在线观看| 亚洲欧洲日产国产| 亚洲欧美日韩卡通动漫| 一本一本综合久久| 国产69精品久久久久777片| 国产综合懂色| 精品少妇黑人巨大在线播放 | 久久鲁丝午夜福利片| 一进一出抽搐gif免费好疼| 国产不卡一卡二| 国产精品国产三级国产av玫瑰| 欧美精品国产亚洲| 少妇裸体淫交视频免费看高清| 免费人成视频x8x8入口观看| 日本欧美国产在线视频| 国产日本99.免费观看| 成年av动漫网址| 国产在线男女| 免费无遮挡裸体视频| 成人毛片a级毛片在线播放| 亚洲成人中文字幕在线播放| 亚洲国产欧洲综合997久久,| 91在线精品国自产拍蜜月| 国产精品嫩草影院av在线观看| 免费看日本二区| 国产在视频线在精品| 麻豆成人av视频| 看片在线看免费视频| 欧美激情久久久久久爽电影| 成人毛片60女人毛片免费| 日韩欧美精品v在线| 亚洲综合色惰| 欧美性感艳星| 一级毛片我不卡| 三级毛片av免费| 久久久久久久久久久丰满| 最近最新中文字幕大全电影3| 成人一区二区视频在线观看| 国产成人freesex在线| 久久久久久九九精品二区国产| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品亚洲一区二区| 听说在线观看完整版免费高清| 久久精品国产99精品国产亚洲性色| 伊人久久精品亚洲午夜| 99热网站在线观看| 久久精品夜色国产| 中文字幕久久专区| 熟妇人妻久久中文字幕3abv| 免费av不卡在线播放| 又粗又爽又猛毛片免费看| 国产一级毛片七仙女欲春2| 久久久午夜欧美精品| 三级经典国产精品| 国产精品人妻久久久影院| 亚洲国产色片| 国产成人一区二区在线| 美女大奶头视频| 国产精品久久久久久av不卡| 成人一区二区视频在线观看| 黄片wwwwww| 毛片一级片免费看久久久久| 又黄又爽又刺激的免费视频.| 夜夜看夜夜爽夜夜摸| 九色成人免费人妻av| 精品无人区乱码1区二区| 亚洲性久久影院| 亚洲成人中文字幕在线播放| 午夜精品在线福利| 国产成人一区二区在线| 在线免费观看的www视频| 亚洲欧美中文字幕日韩二区| 成人高潮视频无遮挡免费网站| 菩萨蛮人人尽说江南好唐韦庄 | 午夜激情福利司机影院| 人人妻人人澡人人爽人人夜夜 | h日本视频在线播放| 少妇人妻一区二区三区视频| 少妇裸体淫交视频免费看高清| 插阴视频在线观看视频| 婷婷六月久久综合丁香| 久久精品影院6| 一级毛片我不卡| 国产高清三级在线| 黄色一级大片看看| 男人舔女人下体高潮全视频| 国产精品爽爽va在线观看网站| 欧美色欧美亚洲另类二区| 国产精品不卡视频一区二区| 午夜福利在线在线| 少妇裸体淫交视频免费看高清| 久久久国产成人精品二区| 亚洲欧美日韩东京热| 男女下面进入的视频免费午夜| 男女边吃奶边做爰视频| 欧美色欧美亚洲另类二区| 国产探花极品一区二区| 婷婷色综合大香蕉| 成人午夜高清在线视频| 美女 人体艺术 gogo| 国产精品电影一区二区三区| 91精品国产九色| 国产精品1区2区在线观看.| 久久精品国产清高在天天线| 深爱激情五月婷婷| 午夜爱爱视频在线播放| 18+在线观看网站| 亚洲乱码一区二区免费版| 岛国在线免费视频观看| 欧美极品一区二区三区四区| 国产人妻一区二区三区在| 青春草视频在线免费观看| 国产高潮美女av| 亚洲欧洲国产日韩| 18禁在线无遮挡免费观看视频| 国产在视频线在精品| 国产一级毛片七仙女欲春2| 一边亲一边摸免费视频| 美女 人体艺术 gogo| 午夜免费男女啪啪视频观看| 久久人人精品亚洲av| av免费在线看不卡| 青青草视频在线视频观看| 99热全是精品| 中文在线观看免费www的网站| 免费无遮挡裸体视频| 99久久中文字幕三级久久日本| 毛片女人毛片| 可以在线观看毛片的网站| 日韩一本色道免费dvd| 国产成人福利小说| 日韩亚洲欧美综合| 亚洲av电影不卡..在线观看| 免费观看精品视频网站| 熟女电影av网| 久久国产乱子免费精品| 在线免费观看的www视频| av卡一久久| 色视频www国产| 看黄色毛片网站| 人妻夜夜爽99麻豆av| 国产精品久久久久久精品电影| 岛国在线免费视频观看| 久久精品国产99精品国产亚洲性色| 亚洲无线观看免费| 国产高潮美女av| 亚洲国产欧美人成| 久久精品夜色国产| 日本-黄色视频高清免费观看| 欧美一区二区亚洲| 看片在线看免费视频| 国产精品av视频在线免费观看| 日韩 亚洲 欧美在线| 黄色日韩在线| 久久久精品94久久精品| 国产日韩欧美在线精品| 午夜a级毛片| 美女 人体艺术 gogo| 亚洲成人av在线免费| 99九九线精品视频在线观看视频| 成人一区二区视频在线观看| 日韩精品有码人妻一区| 久久99热6这里只有精品| 国产精品久久久久久久电影| 女人十人毛片免费观看3o分钟| 亚洲三级黄色毛片| 亚洲综合色惰| 亚洲精品456在线播放app| 日本-黄色视频高清免费观看| 日本免费a在线| 成人二区视频| 精品不卡国产一区二区三区| 一区二区三区四区激情视频 | 久久精品国产自在天天线| 乱系列少妇在线播放| 看十八女毛片水多多多| 偷拍熟女少妇极品色| 国产色婷婷99| 免费电影在线观看免费观看| 高清毛片免费看| 只有这里有精品99| 国产精品美女特级片免费视频播放器| 国内精品一区二区在线观看| 深爱激情五月婷婷| 国产精品人妻久久久久久| 黑人高潮一二区| 成年av动漫网址| 久久99精品国语久久久| 人妻少妇偷人精品九色| av在线蜜桃| 久久韩国三级中文字幕| 国产探花极品一区二区| 亚洲国产欧美人成| 国产精品电影一区二区三区| av又黄又爽大尺度在线免费看 | 天堂网av新在线| 男人舔女人下体高潮全视频| 青春草国产在线视频 | 能在线免费看毛片的网站| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲欧洲日产国产| 一级毛片aaaaaa免费看小| 国产蜜桃级精品一区二区三区| 免费看av在线观看网站| 2021天堂中文幕一二区在线观| 99久国产av精品国产电影| 国产毛片a区久久久久| 激情 狠狠 欧美| 免费人成在线观看视频色| 国产私拍福利视频在线观看| 高清在线视频一区二区三区 | 91麻豆精品激情在线观看国产| 国产精品伦人一区二区| 91aial.com中文字幕在线观看| 国产片特级美女逼逼视频| 小说图片视频综合网站| 亚洲av熟女| 欧美日韩国产亚洲二区| 国产成人福利小说| 成人一区二区视频在线观看| 五月玫瑰六月丁香| 久久人妻av系列| 看非洲黑人一级黄片| 久久精品久久久久久久性| 日韩av不卡免费在线播放| 久久午夜亚洲精品久久| 夜夜看夜夜爽夜夜摸| 欧美成人a在线观看| 天天一区二区日本电影三级| 亚洲av第一区精品v没综合| 97超碰精品成人国产| 哪里可以看免费的av片| 日韩视频在线欧美| 亚洲18禁久久av| 直男gayav资源| 五月伊人婷婷丁香| 晚上一个人看的免费电影| 99国产极品粉嫩在线观看| 国产精品一区www在线观看| 永久网站在线| 久久午夜亚洲精品久久| 国内精品美女久久久久久| 欧美激情在线99| 国产免费男女视频| 欧美一区二区亚洲| 欧美激情国产日韩精品一区| 偷拍熟女少妇极品色| 人人妻人人看人人澡| 色哟哟·www| 国产成人精品一,二区 | 麻豆成人av视频| a级毛片免费高清观看在线播放| 亚洲av成人精品一区久久| 国产精品久久久久久精品电影| 亚洲久久久久久中文字幕| 亚洲精品日韩av片在线观看| 亚洲欧美精品综合久久99| 国产亚洲5aaaaa淫片| 能在线免费观看的黄片| 久久精品人妻少妇| 国产中年淑女户外野战色| 久久热精品热| 少妇人妻精品综合一区二区 | 插阴视频在线观看视频| 欧美三级亚洲精品| 99riav亚洲国产免费| 一边摸一边抽搐一进一小说| 亚洲av二区三区四区| 一个人免费在线观看电影| 尾随美女入室| 日韩成人av中文字幕在线观看| 亚洲欧美成人精品一区二区| 久久综合国产亚洲精品| 伊人久久精品亚洲午夜| 日韩av在线大香蕉| 夫妻性生交免费视频一级片| 久久午夜福利片| 久久久久久久久久久免费av| 国模一区二区三区四区视频| 夜夜看夜夜爽夜夜摸| 亚洲va在线va天堂va国产| 精品不卡国产一区二区三区| 99久久精品国产国产毛片| 欧美区成人在线视频| 午夜老司机福利剧场| 欧美区成人在线视频| 亚洲精品久久国产高清桃花| 全区人妻精品视频| 老司机影院成人| 寂寞人妻少妇视频99o| 午夜精品在线福利| 国内少妇人妻偷人精品xxx网站| 国产真实乱freesex| 一进一出抽搐gif免费好疼| av女优亚洲男人天堂| 成人毛片a级毛片在线播放| 亚洲欧美日韩高清专用| 22中文网久久字幕| 麻豆精品久久久久久蜜桃| 日日啪夜夜撸| 国产精品伦人一区二区| 国内揄拍国产精品人妻在线| 别揉我奶头 嗯啊视频| 草草在线视频免费看| 晚上一个人看的免费电影| 3wmmmm亚洲av在线观看| 亚洲第一电影网av| 亚州av有码| 国产高清三级在线| 国产一区二区在线av高清观看| 可以在线观看毛片的网站| 国产精品三级大全| 免费av观看视频| 国产精品精品国产色婷婷| 亚洲,欧美,日韩| 日本免费a在线| 少妇熟女aⅴ在线视频| 久久6这里有精品| 九草在线视频观看| 欧美xxxx黑人xx丫x性爽| 美女脱内裤让男人舔精品视频 | www.av在线官网国产| 日日摸夜夜添夜夜爱| 不卡视频在线观看欧美| 亚洲aⅴ乱码一区二区在线播放| 欧美激情久久久久久爽电影| videossex国产| 欧美日韩一区二区视频在线观看视频在线 | 日本撒尿小便嘘嘘汇集6| 久久99精品国语久久久| 成人二区视频| 婷婷六月久久综合丁香| 国产成人精品一,二区 | 白带黄色成豆腐渣| 亚洲精品久久国产高清桃花| 悠悠久久av| 日韩 亚洲 欧美在线| 国产极品天堂在线| 嘟嘟电影网在线观看| 最新中文字幕久久久久| 国产在视频线在精品| 国产一区亚洲一区在线观看| 国产av麻豆久久久久久久| a级毛片免费高清观看在线播放| 日韩成人伦理影院| 69av精品久久久久久| 三级经典国产精品| 春色校园在线视频观看| 欧美一区二区精品小视频在线| 日韩 亚洲 欧美在线| 最好的美女福利视频网| 免费黄网站久久成人精品| 国产私拍福利视频在线观看| 美女被艹到高潮喷水动态| 欧美性猛交黑人性爽| 麻豆成人av视频| 日日摸夜夜添夜夜爱| 女同久久另类99精品国产91| 午夜爱爱视频在线播放| 国产午夜精品久久久久久一区二区三区| 晚上一个人看的免费电影| 菩萨蛮人人尽说江南好唐韦庄 | 久久久久久久久久成人| 18禁在线无遮挡免费观看视频| 欧美三级亚洲精品| 精品无人区乱码1区二区| 国产 一区精品| 亚洲图色成人| 又粗又硬又长又爽又黄的视频 | АⅤ资源中文在线天堂| 中国美白少妇内射xxxbb| а√天堂www在线а√下载| h日本视频在线播放| 亚洲四区av| 久久热精品热| 特大巨黑吊av在线直播| 在线免费观看的www视频| 日韩一区二区视频免费看| 97在线视频观看| 日本一本二区三区精品| h日本视频在线播放| 69av精品久久久久久| 大香蕉久久网| 免费看美女性在线毛片视频| 久久久久免费精品人妻一区二区| 成年女人永久免费观看视频| 大又大粗又爽又黄少妇毛片口| 亚洲内射少妇av| 日韩欧美国产在线观看| 少妇被粗大猛烈的视频| 蜜桃久久精品国产亚洲av| 最新中文字幕久久久久| 亚洲av不卡在线观看| 国产三级在线视频| 欧美3d第一页| 日韩中字成人| 午夜精品国产一区二区电影 | 日韩一本色道免费dvd| 在现免费观看毛片| 国产成人精品久久久久久| 高清毛片免费观看视频网站| 免费观看在线日韩| 久久久久久九九精品二区国产| 最近中文字幕高清免费大全6| 国产成人一区二区在线| 久久久久久国产a免费观看| 天堂av国产一区二区熟女人妻| 亚洲精品色激情综合| 成人亚洲精品av一区二区| 国产午夜精品久久久久久一区二区三区| 亚洲成人久久爱视频| 不卡视频在线观看欧美| 我要看日韩黄色一级片| 日韩一区二区视频免费看| 99热网站在线观看| 亚洲国产精品久久男人天堂| 国产 一区精品| 国产中年淑女户外野战色| 成人毛片60女人毛片免费| 给我免费播放毛片高清在线观看| 99视频精品全部免费 在线| av卡一久久| 看十八女毛片水多多多| 精品免费久久久久久久清纯| 国产免费男女视频| 天美传媒精品一区二区| a级一级毛片免费在线观看| 狠狠狠狠99中文字幕| 国产极品精品免费视频能看的| 久久精品国产亚洲网站| 久久精品国产亚洲av香蕉五月| 精品国内亚洲2022精品成人| а√天堂www在线а√下载| 久久精品91蜜桃| 夜夜夜夜夜久久久久| 菩萨蛮人人尽说江南好唐韦庄 | 国产视频首页在线观看| 看免费成人av毛片| av女优亚洲男人天堂| 亚洲成人中文字幕在线播放| 黄色欧美视频在线观看| 一夜夜www| 日本黄色视频三级网站网址| 99热这里只有是精品在线观看| 国产伦精品一区二区三区视频9| 亚洲丝袜综合中文字幕| 免费人成视频x8x8入口观看| 国内少妇人妻偷人精品xxx网站| 中文字幕制服av| 精品午夜福利在线看| 亚洲,欧美,日韩| 禁无遮挡网站| 少妇人妻一区二区三区视频| 久久久精品大字幕| 嫩草影院新地址| av在线蜜桃| 老司机福利观看| 哪里可以看免费的av片| 婷婷亚洲欧美| 精品人妻熟女av久视频| 国产精品av视频在线免费观看| 国产在线男女| 亚洲五月天丁香| 男女啪啪激烈高潮av片| 亚洲国产精品国产精品| 最近视频中文字幕2019在线8| 亚洲成人av在线免费| 久久精品国产亚洲网站| 久久中文看片网| av在线亚洲专区| 免费观看的影片在线观看| 欧美日本亚洲视频在线播放| 国产大屁股一区二区在线视频| 欧美性感艳星| 成人特级黄色片久久久久久久| 国产视频内射| av又黄又爽大尺度在线免费看 | 小蜜桃在线观看免费完整版高清| 嫩草影院精品99| 少妇熟女aⅴ在线视频| 亚洲欧美清纯卡通| 亚洲激情五月婷婷啪啪| 久久草成人影院| 成人国产麻豆网| 高清在线视频一区二区三区 | 一区二区三区免费毛片| 久久99精品国语久久久| 99热这里只有是精品50| 亚洲七黄色美女视频| 男的添女的下面高潮视频| 欧美另类亚洲清纯唯美| 欧美一区二区国产精品久久精品| 波野结衣二区三区在线| 久久人人爽人人片av| 美女国产视频在线观看| 伦精品一区二区三区| 好男人在线观看高清免费视频| 婷婷色综合大香蕉| 看片在线看免费视频| 国模一区二区三区四区视频| 亚洲乱码一区二区免费版| 亚洲在线自拍视频| 岛国在线免费视频观看| 午夜福利高清视频| 亚洲国产精品国产精品| 少妇裸体淫交视频免费看高清| 日本黄大片高清| 蜜桃久久精品国产亚洲av| 久久久精品欧美日韩精品| 全区人妻精品视频| 女的被弄到高潮叫床怎么办| 一级毛片电影观看 | 亚洲精品乱码久久久久久按摩| 亚洲天堂国产精品一区在线| 国产黄色小视频在线观看| 三级国产精品欧美在线观看| 内地一区二区视频在线| 舔av片在线| 成年女人看的毛片在线观看| 91精品国产九色| 成人特级av手机在线观看| 你懂的网址亚洲精品在线观看 | 亚洲av二区三区四区| 久久99热6这里只有精品| 中文精品一卡2卡3卡4更新| 网址你懂的国产日韩在线| 桃色一区二区三区在线观看| 天堂√8在线中文| av福利片在线观看| 桃色一区二区三区在线观看| 天堂√8在线中文| 在线观看66精品国产| 啦啦啦啦在线视频资源| 亚洲欧洲日产国产| 成人特级黄色片久久久久久久| 一卡2卡三卡四卡精品乱码亚洲| 中文字幕熟女人妻在线| 三级国产精品欧美在线观看| 日本在线视频免费播放| 青青草视频在线视频观看| 在线观看av片永久免费下载|