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

    中華獼猴桃WOX轉(zhuǎn)錄因子的生物信息學(xué)分析

    2021-10-16 14:29:00付春黃春莉江納楊瑤君
    廣西植物 2021年9期

    付春 黃春莉 江納 楊瑤君

    摘 要:WOX基因家族是植物特有的一類轉(zhuǎn)錄因子,在植物發(fā)育關(guān)鍵時(shí)期,如胚的形成、維持干細(xì)胞穩(wěn)定性和器官形成過(guò)程中發(fā)揮重要調(diào)控作用。為探究WOX轉(zhuǎn)錄因子調(diào)控中華獼猴桃的生長(zhǎng)發(fā)育的分子理論基礎(chǔ),該文利用ProtParam、Cell-PLoc2.0、SignalP4.1 Server等在線軟件對(duì)中華獼猴桃的16個(gè)WOX轉(zhuǎn)錄因子進(jìn)行理化特性、亞細(xì)胞定位和信號(hào)肽等詳細(xì)的生物信息學(xué)分析。結(jié)果表明:(1)中華獼猴WOX轉(zhuǎn)錄因子的所有蛋白成員的氨基酸數(shù)目在138~371之間,分子量為16.222~42.185 kD,均定位于細(xì)胞核的不穩(wěn)定親水蛋白。(2)中華獼猴桃的16個(gè)WOX蛋白成員均屬于Homodomain超家族,僅Achn362451是分泌蛋白,且僅Achn362451和Achn141001具有跨膜區(qū)。(3)大部分WOX蛋白成員的潛在磷酸化位點(diǎn)都位于絲氨酸處。(4)其二級(jí)結(jié)構(gòu)主要以無(wú)規(guī)則卷曲為主,其次為α-螺旋。(5)保守基序分析結(jié)果表明中華獼猴桃的16個(gè)WOX蛋白成員均含有motif1。(6)系統(tǒng)進(jìn)化分析結(jié)果表明中華獼猴桃與菠蘿的親緣關(guān)系最近,與番茄的親緣關(guān)系最遠(yuǎn)。(7)在中華獼猴桃果實(shí)成熟過(guò)程中不同時(shí)期的轉(zhuǎn)錄組分析結(jié)果表明其WOX基因家族中4個(gè)成員Achn131681、Achn145561、Achn336591、Achn362451基因在授粉后20、120、127 d表達(dá)含量都較高,其他12個(gè)WOX基因成員表達(dá)含量低甚至不表達(dá)。該研究為進(jìn)一步研究中華獼猴桃WOX轉(zhuǎn)錄因子在其生長(zhǎng)發(fā)育過(guò)程中的分子機(jī)理提供了一定的參考依據(jù)。

    關(guān)鍵詞:中華獼猴桃,WOX,轉(zhuǎn)錄因子,基因家族,生物信息學(xué)分析

    中圖分類號(hào):Q943.2

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1000-3142(2021)09-1561-15

    Abstract:WOX gene family is a kind of transcription factor peculiar to plants,and plays an important regulatory role in the critical stages of plant development,such as embryo formation,stem cell stability and organ formation. In order to explore the molecular basis of WOX transcription factors regulating the growth and development of Actinidia chinensis,the physical and chemical characteristics,subcellular localization and signal peptide of 16 WOX transcription factors in A. chinensis were analyzed by ProtParam,Cell-PLoc2.0,SignalP4.1 Server online software in detail in this study. The results were as follows:(1)The amino acid number of all protein members of WOX transcription factors in A. chinensis ranged from 138 to 371,with molecular weight ranging from 16.222 to 42.185 kD,all of which were unstable hydrophilic proteins located in nucleus. (2)The 16 WOX protein members of A. chinensis belonged to Homodomain superfamily,only Achn362451 was secreted protein,and only Achn362451 and Achn141001 had transmembrane regions. (3)The potential phosphorylation sites of most WOX members were located in serine site. (4)The secondary structure of WOX proteins was mainly random coil,followed by α-helix. (5)Conservative motifs analysis results showed that all 16 WOX protein members of A. chinensis contained motif1. (6)Phylogenetic analysis results showed that the genetic relationship between A. chinensis and Ananas comosus was the closest,and that between Actinidia chinensis and Solanum lycopersicum was the farthest. (7)The results of transcriptome analysis in different stages of fruit ripening process of Actinidia chinensis showed that four members of the WOX gene family,Achn131681,Achn145561,Achn336591,Achn362451,were highly expressed on 20,120,127 d after pollination,while those of the other 12 WOX gene members were low or even not expressed. This study provides a reference basis for further study on the molecular mechanism of WOX transcription factors in the growth and development of A. chinensis.

    Key words:Actinidia chinensis,WOX,transcription factor,gene family,bioinformatics analysis

    中華獼猴桃(Actinidia chinensis),多年生藤本植物,是杜鵑花目(Ericales)獼猴桃科(Actinidiaceae)的一種經(jīng)濟(jì)樹木。中華獼猴桃原產(chǎn)自我國(guó),果實(shí)富含維生素C,且維生素C在人體中的利用率高達(dá)94%(宋圃菊等,1984),一個(gè)成年人每天只攝食一個(gè)鮮獼猴桃就可滿足人體對(duì)維生素C的需求,因而被譽(yù)為“維C之王”。果實(shí)酸甜可口,除了鮮食,還可加工成各種食品和飲料。臨床試驗(yàn)發(fā)現(xiàn),獼猴桃汁能有效抑制麻風(fēng)桿菌,從而明顯改善麻風(fēng)病的臨床癥狀。除此之外,獼猴桃汁所含的一種抗突變成分谷胱甘肽,可抑制誘發(fā)致癌基因的突變,并有效阻斷致癌物質(zhì)亞硝酸基嗎啉的合成,從而防止癌癥的發(fā)生(《湖北林業(yè)科技》編輯部,2007)。

    WUSCHEL(WUS)相關(guān)同源框(WOX)蛋白家族成員特異存在于植物中,在莖頂端分生組織、根頂端分生組織和形成層中干細(xì)胞的維持和增殖等生長(zhǎng)發(fā)育中發(fā)揮重要作用。根據(jù)系統(tǒng)進(jìn)化分析,模式植物擬南芥(Arabidopsis thaliana)中15個(gè)WOX轉(zhuǎn)錄因子被分為以ATWOX10,13,14在內(nèi)的遠(yuǎn)古進(jìn)化支,以ATWOX8,9,11,12在內(nèi)的中間進(jìn)化支,和以WUS,ATWOX1-7在內(nèi)的現(xiàn)代進(jìn)化支/WUS進(jìn)化支(高麗等,2015)。其中,遠(yuǎn)古進(jìn)化支起源于藻類植物,中間進(jìn)化支起源于維管植物,現(xiàn)代進(jìn)化支起源于種子植物(Lian et al.,2014)。隨著對(duì)WOX轉(zhuǎn)錄因子的深入研究,WOX蛋白家族在各個(gè)植物生長(zhǎng)發(fā)育中的功能及其作用機(jī)理逐漸明確,特別是在模式植物擬南芥中,大多數(shù)WOX蛋白成員的功能及作用機(jī)理都已明確。ATWOX1的過(guò)度表達(dá)導(dǎo)致擬南芥分生組織發(fā)育和多胺穩(wěn)態(tài)的缺陷(Zhang et al.,2011)。ATWOX2和ATWOX8在卵細(xì)胞和合子中共表達(dá),但由于其不對(duì)稱分裂分別僅局限于合子的頂端和基部子細(xì)胞(Lie et al.,2012)。ATWOX3調(diào)控旁側(cè)器官(葉和花)的發(fā)育,其突變體會(huì)出現(xiàn)花瓣缺失或變窄(Shimizu et al.,2008)。ATWOX4涉及形成層中維管束的干細(xì)胞的維持。ATWOX5在根尖靜止中心特異性表達(dá),主要維持根尖分生組織(RAM)中干細(xì)胞的活性(Haecker,2004)。ATWOX6主要在胚珠中表達(dá),在珠被和卵細(xì)胞的形成過(guò)程中可防止細(xì)胞提前分化,進(jìn)而調(diào)控胚珠的生長(zhǎng)發(fā)育(Park SO,2005);Zhu et al.(2004)指出ATWOX6可影響冷脅迫的應(yīng)答。ATWOX9可識(shí)別發(fā)育信號(hào)并調(diào)控細(xì)胞周期(Wu et al.,2005)。ATWOX10的具體作用雖未明確,但已有研究發(fā)現(xiàn)ATWOX10可能受到光信號(hào)、生長(zhǎng)素、赤霉素和脫落酸等植物激素以及脅迫等多種信號(hào)途徑的協(xié)同作用(王占軍等,2017)。ATWOX11和ATWOX12共同影響外植體根的發(fā)育(Liu et al.,2014)。ATWOX13和ATWOX14在根部及花藥中表達(dá),主要防止其提前分化(Deveaux et al.,2008)。

    關(guān)于植物基因組中的WOX基因家族鑒定和生物信息學(xué)分析的報(bào)道較多。利用隱馬爾科夫模型對(duì)番茄(Solanum lycopersicum)基因組進(jìn)行WOX基因家族鑒定并發(fā)現(xiàn)有10個(gè)WOX基因成員,同時(shí)分析其WOX蛋白結(jié)構(gòu)域和基因結(jié)構(gòu)及系統(tǒng)進(jìn)化關(guān)系(李曉旭等,2016)。蝴蝶蘭(Phalaenopsis aphrodite)和石斛蘭(Dendrobium catenatum)基因組中分別含有14個(gè)和10個(gè)WOX基因成員,其中,蝴蝶蘭有3個(gè)WOX基因成員出現(xiàn)基因重復(fù),而在石斛蘭基因組中沒(méi)有出現(xiàn)WOX基因重復(fù)(Ramkumar et al.,2018)。茶樹(Camellia sinensis)基因組中有18個(gè)WOX基因成員,且其系統(tǒng)進(jìn)化樹分為3組,分別是古老分支、中間分支和現(xiàn)代分支三類(Wang et al.,2019)。在甘藍(lán)型油菜、白菜和甘藍(lán)基因組序列中分別鑒定出52、25、29個(gè)WOX基因成員(Li et al.,2019)。大豆(Glycine max)基因組中含有33個(gè)WOX基因成員,并定位在19條染色體上,而16號(hào)染色體沒(méi)有WOX基因(Hao et al.,2019)。對(duì)毛白楊(Populus tomentosa)的PtoWUSa、PtoWOX4a和PtoWOX5a三個(gè)基因編碼蛋白質(zhì)的結(jié)構(gòu)域和在不同組織中的表達(dá)進(jìn)行了詳細(xì)分析(李建波,2017;王爽等,2019)。在玉米(Zea mays)基因組中鑒定出31個(gè)WOX基因成員,且分析了其編碼蛋白的理化特性和系統(tǒng)進(jìn)化關(guān)系(鄭玲等,2019)。毛竹(Phyllostachys heterocycla)基因組中有27個(gè)WOX基因成員(李晨曦等,2016),而在其他竹類植物基因組中,俄莉竹(Olyra latifolia)有12個(gè)WOX基因成員,Raddia guianensis含有9個(gè)WOX基因成員,瓜多竹(Guadua angustifolia)含有18個(gè)WOX基因成員以及蕓香竹(Bonia amplexicaulis)含有26個(gè)WOX基因成員(Li et al.,2020)。在葡萄(Vitis vinifera)基因組中鑒定出10個(gè)WOX基因成員,且其與擬南芥的進(jìn)化樹分為祖先分支、中間分支和現(xiàn)代分支三類(王鵬飛等,2018)。小麥(Triticum aestivum)基因組中含有80個(gè)WOX基因成員,其中定位在A、B、D基因組上的分別有23、26、29個(gè),另有2個(gè)WOX基因成員定位不明確(李仲卿,2018)。蘋果(Malus pumila)基因組中含有12個(gè)WOX基因成員,且在果實(shí)中的表達(dá)含量較高(王楚堃等,2020)。小桐子(Jatropha curcas)基因組中有12個(gè)WOX基因成員,且在不同組織中的表達(dá)呈顯著性差異(唐躍輝等,2019)。棉花(Gossypium hirsutum)基因組中有37個(gè)WOX基因成員,其中,26個(gè)WOX基因在棉花根、莖、葉、花等器官中都有表達(dá)(呂有軍等,2017)。簸箕柳(Salix suchowensis)基因組中被鑒定出15個(gè)WOX基因成員,且定位在20條染色體上(Wang et al.,2019)。在人參(Panax ginseng)基因組中鑒定出10個(gè)WOX基因成員(Liu et al.,2019)。雖然關(guān)于植物基因組中WOX基因家族的鑒定、功能分析的研究較多,但對(duì)于中華獼猴桃基因組中WOX基因家族的鑒定及其分子功能研究的文獻(xiàn)報(bào)道較少。

    關(guān)于中華獼猴桃基因組中轉(zhuǎn)錄因子的研究報(bào)道有R1R2R3-MYB轉(zhuǎn)錄因子增強(qiáng)其干旱脅迫反應(yīng)(Ya et al.,2019),EIN3/EIL轉(zhuǎn)錄因子的鑒定,系統(tǒng)進(jìn)化及其基因表達(dá)分析(陳強(qiáng)等,2018),MADS-box轉(zhuǎn)錄因子與中華獼猴桃乙烯促成熟相關(guān)(Mcatee et al.,2015),還有bZIP家族的相關(guān)報(bào)道(于豐源等,2020),但未見中華獼猴桃基因組中WOX轉(zhuǎn)錄因子鑒定和生物信息學(xué)研究的報(bào)道。為研究中華獼猴桃WOX基因家族編碼蛋白質(zhì)的分子功能及其在生長(zhǎng)發(fā)育過(guò)程中的作用如何和非生物脅迫條件下的應(yīng)激反應(yīng)如何。本研究擬對(duì)中華獼猴桃16個(gè)WOX基因成員進(jìn)行理化特性、亞細(xì)胞定位、蛋白高級(jí)結(jié)構(gòu)及系統(tǒng)進(jìn)化等詳細(xì)的生物信息學(xué)分析,為后續(xù)進(jìn)一步研究中華獼猴桃的生長(zhǎng)發(fā)育和分子功能等提供一定的參考依據(jù)和理論基礎(chǔ)。

    1 材料與方法

    1.1 材料

    中華獼猴桃的16個(gè)WOX轉(zhuǎn)錄因子的蛋白序列來(lái)自于PlantTFDB(http://planttfdb.cbi.pku.edu.cn/)數(shù)據(jù)庫(kù),并通過(guò)該數(shù)據(jù)庫(kù)下載與中華獼猴桃WOX轉(zhuǎn)錄因子成員Achn172381同源的番茄、矮牽?;ǎ≒etunia inflata)、辣椒(Capsicum annuum)、擬南芥、油菜(Brassica napus)、葡萄、菠菜(Spinacia oleracea)、大豆、菠蘿(Ananas comosus)、小麥、毛竹、玉米、水稻(Oryza rufipogon)、花生(Arachis hypogaea)、馬鈴薯(Solanum tuberosum)的WOX轉(zhuǎn)錄因子成員的蛋白序列用于不同物種間WOX轉(zhuǎn)錄因子的系統(tǒng)進(jìn)化關(guān)系分析研究。

    1.2 方法

    1.2.1 理化性質(zhì)分析 利用在線程序ProtParam (https://web.expasy.org/protparam/)對(duì)中華獼猴桃的WOX基因進(jìn)行理化性質(zhì)分析,包括分子量、氨基酸數(shù)、等電點(diǎn)、帶負(fù)電荷殘基數(shù)、帶正電荷殘基數(shù)、不穩(wěn)定指數(shù)、脂肪指數(shù),以及親疏水性的預(yù)測(cè)分析。并利用數(shù)據(jù)庫(kù)PlantTFDB(http://planttfdb.cbi.pku.edu.cn/blast.php)進(jìn)行Blastp,分別找出中華獼猴桃16個(gè)WOX轉(zhuǎn)錄因子在模式植物擬南芥對(duì)應(yīng)的基因。通過(guò)Cell-PLoc2.0 (http://www.csbio.sjtu.edu.cn/bioinf/ Cell-PLoc-2/)對(duì)該家族成員進(jìn)行亞細(xì)胞定位分析。

    1.2.2 保守結(jié)構(gòu)域、信號(hào)肽的分析及跨膜域的預(yù)測(cè) 利用NCBI(https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi)預(yù)測(cè)WOX蛋白的保守結(jié)構(gòu)域;利用SignalP4.1 Server (http://www.cbs.dtu.dk/services/SignalP/)分析WOX蛋白的信號(hào)肽;在TMHMM Severv.2.0 (http://www.cbs.dtu.dk/services/TMHMM/)在線程序輸入中華獼猴桃16個(gè)WOX轉(zhuǎn)錄因子蛋白序列進(jìn)行跨膜域預(yù)測(cè)分析。

    1.2.3 磷酸化位點(diǎn)的預(yù)測(cè)分析 通過(guò)在線工具NetPhos 3.1 Server (http://www.cbs.dtu.dk/services/NetPhos/)對(duì)中華獼猴桃的WOX轉(zhuǎn)錄因子潛在磷酸化位點(diǎn)進(jìn)行預(yù)測(cè)分析。

    1.2.4 二級(jí)結(jié)構(gòu)、三級(jí)結(jié)構(gòu)的分析預(yù)測(cè) 使用SOPMA程序(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html)對(duì)中華獼猴桃WOX轉(zhuǎn)錄因子的蛋白質(zhì)二級(jí)結(jié)構(gòu)進(jìn)行預(yù)測(cè)分析;利用在線工具SWISS-MODEL(https://swissmodel.expasy.org/)對(duì)其蛋白質(zhì)的三級(jí)結(jié)構(gòu)進(jìn)行預(yù)測(cè)分析。

    1.2.5 保守基序分析 利用在線工具M(jìn)EME(http://meme-suite.org/tools/meme)對(duì)中華獼猴桃的WOX蛋白家族進(jìn)行保守基序分析,參數(shù)設(shè)置:motifs的數(shù)值設(shè)為10,其余參數(shù)默認(rèn)。

    1.2.6 序列的多重比對(duì)及進(jìn)化樹的構(gòu)建 利用中華獼猴桃的Achn172381蛋白序列在PlantTFDB(http://planttfdb.cbi.pku.edu.cn/blast.php)數(shù)據(jù)庫(kù)進(jìn)行Blastp,分別找到番茄、辣椒、油菜、葡萄、玉米、擬南芥等其他15物種中的同源蛋白序列,并分別下載下來(lái)。通過(guò)軟件MEGA-X(Sudhir et al.,2018)中的ClutalW比對(duì)功能對(duì)中華獼猴桃的Achn172381蛋白序列與擬南芥等其他15種植物的同源蛋白序列進(jìn)行多重比對(duì),參數(shù)默認(rèn);并用鄰接法(nerghbor-join,NJ)構(gòu)建進(jìn)化樹,Bootstrap replications值設(shè)定為1 000。

    1.2.7 中華獼猴桃在不同生長(zhǎng)發(fā)育時(shí)期的WOX基因表達(dá)分析 中華獼猴桃授粉后20、120、127 d的轉(zhuǎn)錄組原始數(shù)據(jù)下載于NCBI(https://www.ncbi.nlm.nih.gov/)數(shù)據(jù)庫(kù)(SRA No.SRA065642)(Huang et al.,2013),利用Tang et al.(2016)對(duì)該轉(zhuǎn)錄組原始數(shù)據(jù)進(jìn)行注釋分析其授粉20、120、127 d后的表達(dá)譜結(jié)果,通過(guò)篩選中華獼猴桃WOX基因家族的16個(gè)基因成員表達(dá)水平FPKM值分析其WOX基因表達(dá)水平情況。

    2 結(jié)果與分析

    2.1 中華獼猴桃WOX蛋白的理化性質(zhì)分析

    利用ProtParam預(yù)測(cè)分析中華獼猴桃16個(gè)WOX轉(zhuǎn)錄因子的理化性質(zhì),分析結(jié)果顯示中華獼猴桃WOX轉(zhuǎn)錄因子的氨基酸范圍為138(Achn107771)~371(Achn363401);分子量大小范圍為16.222(Achn107771)~42.185 kD(Achn363401);等電點(diǎn)為5.14(Achn362341)~9.79(Achn141001),其中,Achn053291、Achn362341、Achn104531、Achn362451、Achn145561的等電點(diǎn)均小于7,它們可能編碼弱酸性蛋白,在酸性亞細(xì)胞環(huán)境中發(fā)揮作用;不穩(wěn)定指數(shù)為46.01~74.15,表明均為不穩(wěn)定蛋白(不穩(wěn)定指數(shù)>40為不穩(wěn)定蛋白);脂肪指數(shù)為44.53(Achn053291)~79.49(Achn141001),說(shuō)明蛋白成員間熱穩(wěn)定性差異較大。負(fù)電荷的氨基酸殘基數(shù)(Asp + Glu)為13(Achn107771)~37(Achn362451),正電荷的氨基酸殘基數(shù)(Arg + Lys)為14(Achn107771)~42(Achn363401),其中,除Achn053291、Achn362341、Achn104531、Achn362451、Achn145561這5個(gè)成員的負(fù)電荷的氨基酸殘基總數(shù)大于正電荷的氨基酸殘基總數(shù)外,其余13個(gè)WOX蛋白成員的正電荷的氨基酸殘基總數(shù)大于負(fù)電荷的氨基酸殘基總數(shù)。亞細(xì)胞定位分析表明,中華獼猴桃的16個(gè)WOX蛋白成員均位于細(xì)胞核(表1)。

    根據(jù)各個(gè)WOX蛋白成員的氨基酸組成不同,主要分為a、b、c、d、e 五類:a類以絲氨酸為主,包括Achn172381、Achn220151、Achn285781、Achn207451、

    Achn053291、Achn104531、Achn145561、Achn172571、Achn362341、Achn141001這10個(gè)成員,其中Achn053291絲氨酸的含量最高,為15.20%;b類中僅有Achn362451這1個(gè)成員,其氨基酸組成以絲氨酸和谷氨酰胺為主,絲氨酸和谷氨酰各占9.30%;c類以亮氨酸為主,包括Achn131681、Achn107771、Achn336591這3個(gè)成員,其中Achn107771亮氨酸最高,占10.90%;d類中以亮氨酸和蘇氨酸為主,僅含Achn299661這1個(gè)成員,亮氨酸和蘇氨酸各占10.40%;e類中以蘇氨酸為主,只含Achn363401這1個(gè)成員,蘇氨酸占為11.10%(圖1)。

    對(duì)中華獼猴桃WOX轉(zhuǎn)錄子進(jìn)行氨基酸親疏水性Protscole預(yù)測(cè)發(fā)現(xiàn)中華獼猴桃的16個(gè)WOX蛋白成員的GRAVY值均為負(fù),即這16個(gè)成員均為親水蛋白。其中總平均親水指數(shù)最大的為Achn053291,其值為-0.954,該成員的第8、第9、第10、第11、第12、第13、第14位氨基酸表現(xiàn)為最大親水性,其值為-3.5;第152位氨基酸表現(xiàn)為最大疏水性,其值為1.533。總平均親水指數(shù)最小的是Achn104531,其值為-0.402,該成員的第84、第85、第86位氨基酸表現(xiàn)為最大親水性,其值為-3.567;第250位氨基酸表現(xiàn)為最大疏水性,其值為2.5(表2,圖2)。

    2.2 保守結(jié)構(gòu)域預(yù)測(cè)、信號(hào)肽的分析及跨膜域的預(yù)測(cè)分析

    中華獼猴桃的16個(gè)WOX轉(zhuǎn)錄因子的保守結(jié)構(gòu)域預(yù)測(cè)分析表明中華獼猴桃的16個(gè)WOX轉(zhuǎn)錄因子中,僅Achn145561有Homodomain superfamily和PIG-H superfamily兩個(gè)保守結(jié)構(gòu)域,其余15個(gè)WOX蛋白成員只含Homodomain superfamily這一個(gè)保守結(jié)構(gòu)域。中華獼猴桃WOX蛋白成員的保守結(jié)構(gòu)域長(zhǎng)度為38~61,其中,Achn172381、Achn131681、Achn299661、Achn363401這4個(gè)成員的保守結(jié)構(gòu)域等長(zhǎng),均為59,但位置各異;Achn207451與Achn053291的保守結(jié)構(gòu)域等長(zhǎng),均為60,且位置吻合;Achn362341與Achn104531的保守結(jié)構(gòu)域長(zhǎng)度和位置均吻合;Achn362451與Achn145561中的Homodomain superfamily保守結(jié)構(gòu)域等長(zhǎng),均為58,但位置不同(表3,圖3)。

    SignalP4.1 Server信號(hào)肽預(yù)測(cè)分析結(jié)果顯示除Achn141001外,其余15個(gè)中華獼猴桃轉(zhuǎn)錄因子均無(wú)信號(hào)肽。因此,推測(cè)Achn141001為分泌蛋白,其余均為非分泌蛋白(圖4)。TMHMM Severv.2.0跨膜域分析結(jié)果顯示,除Achn362451、Achn141001有明顯的跨膜結(jié)構(gòu)域外,其余14個(gè)WOX轉(zhuǎn)錄因子均無(wú)明顯的跨膜結(jié)構(gòu)域。因此,推測(cè)Achn362451、Achn141001為跨膜蛋白(圖5)。

    2.3 磷酸化位點(diǎn)的預(yù)測(cè)分析

    NetPhos 3.1磷酸化位點(diǎn)分析表明16個(gè)WOX蛋白成員中,Achn172381、Achn336591、Achn285781、Achn207451、Achn053291等13個(gè)WOX蛋白成員中磷酸化位點(diǎn)的總數(shù):絲氨酸>蘇氨酸>酪氨酸;Achn131681和Achn299661中磷酸化位點(diǎn)總數(shù):蘇氨酸>絲氨酸>酪氨酸;Achn363401中磷酸化位點(diǎn)總數(shù):絲氨酸=蘇氨酸>酪氨酸。就單個(gè)成員而言,除Achn299661的潛在磷酸化位點(diǎn)為蘇氨酸外,其余15個(gè)WOX蛋白成員潛在磷酸化位點(diǎn)均為絲氨酸。在16個(gè)WOX蛋白成員中,極易在絲氨酸處磷酸化的成員為Achn207451、Achn053291、Achn172571(0.998),極易在蘇氨酸處磷酸化的成員為Achn362341和Achn104531(0.995),極易在酪氨酸處磷酸化的成員為Achn362451(0.978)(表4,圖6)。

    2.4 二級(jí)結(jié)構(gòu)和三級(jí)結(jié)構(gòu)的預(yù)測(cè)分析

    使用SOPMA程序?qū)χ腥A獼猴桃WOX轉(zhuǎn)錄因子的二級(jí)結(jié)構(gòu)進(jìn)行預(yù)測(cè)分析。結(jié)果顯示,大部分WOX蛋白成員的二級(jí)結(jié)構(gòu)都以無(wú)規(guī)則卷曲為主,其次為α-螺旋和延伸鏈,β-轉(zhuǎn)角含量最少。其中,Achn172381、Achn131681、Achn299661、Achn336591、Achn285781等12個(gè)成員的二級(jí)結(jié)構(gòu)含量為無(wú)規(guī)則卷曲>α-螺旋>延伸鏈>β-轉(zhuǎn)角;Achn220151、Achn220151和Achn104531三個(gè)成員的二級(jí)結(jié)構(gòu)含量為無(wú)規(guī)則卷曲>延伸鏈>α-螺旋>β-轉(zhuǎn)角;Achn141001的二級(jí)結(jié)構(gòu)含量為無(wú)規(guī)則卷曲=α-螺旋>延伸鏈>β-轉(zhuǎn)角。α-螺旋在成員中占比為12.79%(Achn362341)~42.00%(Achn362451);延伸鏈在成員中占比為4.00%(Achn362451)~21.18%(Achn104531);β-轉(zhuǎn)角在成員中占比為2.38%(Achn172381)~9.49%(Achn141001);無(wú)規(guī)則卷曲在成員中占比為69.38%(Achn362341)~37.34%(Achn141001)(表5,圖7)。

    通過(guò)在線工具SWISS-MODEL對(duì)中華獼猴桃WOX轉(zhuǎn)錄因子進(jìn)行蛋白質(zhì)的三級(jí)結(jié)構(gòu)預(yù)測(cè),結(jié)果顯示,16個(gè)WOX蛋白成員的蛋白質(zhì)三級(jí)結(jié)構(gòu)可分為11類。Achn172381與Achn220151為一類,其代表為Achn220151;Achn131681與Achn336591為一類,其代表為Achn131681;Achn207451、Achn053291、Achn172571為一類,其代表為Achn053291;Achn362341與Achn104531為一類,其代表為Achn362341;Achn299661、Achn285781、Achn363401、Achn362451、Achn145561、Achn107771、Achn141001各代表一類(圖8)。

    2.5 保守基序分析

    利用在線工具M(jìn)EME對(duì)中華獼猴桃的WOX家族進(jìn)行保守基序分析,結(jié)果表明:motif1出現(xiàn)在中華獼猴桃的所有的WOX蛋白成員中,motif2出現(xiàn)在除Achn107771和Achn141001外的其他WOX蛋白成員中。根據(jù)不同成員中所含motif不同,將中華獼猴桃的16個(gè)WOX蛋白成員分為A、B、C三個(gè)組。A組主要含motif1、motif2、motif3、motif5和motif10的蛋白成員有Achn104531與Achn362341,只含motif1、motif2、motif3和motif5的有Achn362451與Achn145561,只含motif1、motif2、motif3和motif10的有Achn220151;B組主要含motif1、motif2、motif4、motif7、motif8和motif9的有Achn131681與Achn336591,只含motif1、motif2、motif7和motif9的有Achn299661,只含motif1、motif2和motif4的有Achn172381與Achn363401;C組主要含motif1、motif2、motif4和motif6的有Achn285781與Achn172571,含有motif1、motif2、motif4和motif7的有Achn053291與Achn207451,只含motif1和motif4的有Achn141001,只含motif1的為Achn107771(圖9)。

    2.6 中華獼猴桃WOX基因家族在不同時(shí)期的表達(dá)水平分析

    通過(guò)對(duì)中華獼猴桃授粉后20、120、127 d的轉(zhuǎn)錄組表達(dá)分析發(fā)現(xiàn)其WOX基因家族中Achn131681、Achn145561、Achn336591和Achn362451成員在中華獼猴桃授粉后20 d的基因表達(dá)水平顯著提高,成員Achn362451在中華獼猴桃16個(gè)WOX基因成員中表達(dá)含量最高。在授粉后120 d,成員Achn131681、Achn145561和Achn362451基因表達(dá)水平顯著提高,Achn362451表達(dá)含量最高。在授粉后127 d,成員Achn145561和Achn362451的表達(dá)含量較高,且Achn362451表達(dá)含量在所有WOX基因中表達(dá)含量中最高。從中華獼猴桃授粉后20、120、127 d的生長(zhǎng)發(fā)育成熟過(guò)程來(lái)看,成員Achn131681基因表達(dá)水平是先上升后下降;成員Achn145561基因表達(dá)水平一直上升,在授粉后127 d達(dá)最高水平;成員Achn336591基因表達(dá)水平是先上升后下降,在授粉后127 d幾乎不表達(dá);成員Achn362451基因表達(dá)水平是先上升后下降,但在授粉127 d后表達(dá)含量也較高。而除以上4個(gè)基因成員外,其他12個(gè)WOX基因在其授粉后20、120、127 d表達(dá)含量很低甚至不表達(dá)(圖10)。

    2.7 系統(tǒng)發(fā)育樹構(gòu)建

    中華獼猴桃與擬南芥WOX轉(zhuǎn)錄因子構(gòu)建的系統(tǒng)進(jìn)化樹結(jié)果表明,其進(jìn)化樹分為三大分支,與擬南芥的三大進(jìn)化支即遠(yuǎn)古進(jìn)化支、中間進(jìn)化支和現(xiàn)代進(jìn)化支/WUS進(jìn)化支相吻合。因此,推測(cè)Achn145561和Achn362451為古代進(jìn)化支,Achn220151、Achn362341、Achn104531為中間進(jìn)化支,Achn172381、Achn107771、Achn285761、Achn172571、Acn299661、Achn131681、Achn336591、Achn363401、Achn207451、Achn053291、Achn141001為現(xiàn)代進(jìn)化支,即WUS進(jìn)化支(圖11)。

    對(duì)中華獼猴桃、擬南芥、番茄等16種植物的WOX基因家族的蛋白序列進(jìn)行多重比對(duì)并構(gòu)建進(jìn)化樹。其進(jìn)化樹可分為GroupⅠ、GroupⅡ、Group Ⅲ三組,分別包括番茄、辣椒、矮牽?;āM南芥、油菜、葡萄、菠菜、大豆在內(nèi)的Group I,包括中華獼猴桃、菠蘿、玉米、水稻、小麥以及毛竹在內(nèi)的Group Ⅱ,包括馬鈴薯及花生在內(nèi)的Group Ⅲ。分支的遠(yuǎn)近及距離長(zhǎng)短表明,中華獼猴與菠蘿的親緣關(guān)系最近,與番茄的親緣關(guān)系最遠(yuǎn)(圖12)。

    3 討論與結(jié)論

    本研究對(duì)中華獼猴桃的16個(gè)WOX轉(zhuǎn)錄因子進(jìn)行了蛋白質(zhì)理化性質(zhì)、保守結(jié)構(gòu)域預(yù)測(cè)、信號(hào)肽、跨膜結(jié)構(gòu)域預(yù)測(cè)、磷酸化位點(diǎn)、保守基序、二級(jí)結(jié)構(gòu)分析、三級(jí)結(jié)構(gòu)分析和進(jìn)化樹等詳細(xì)生物信息學(xué)分析。王楚堃等(2019)對(duì)蘋果WOX轉(zhuǎn)錄因子研究發(fā)現(xiàn)蘋果中WOX蛋白成員的氨基酸數(shù)量為176~381,分子量為20.122~45.337 kD。李雪垠等(2019)對(duì)谷子和狗尾草WOX轉(zhuǎn)錄因子研究發(fā)現(xiàn),谷子中WOX蛋白成員的氨基酸數(shù)量為212~533,分子量為20.011~54.998 kD;狗尾草中WOX蛋白成員的氨基酸數(shù)量為212~531,分子量為23.969~54.788 kD。郭亞靜等(2020)對(duì)星油藤WOX轉(zhuǎn)錄因子的研究中發(fā)現(xiàn)星油藤WOX蛋白成員氨基酸數(shù)量為220~420,分子量為24.95~45.37 kD。同時(shí),在對(duì)小桐子(唐躍輝等,2019)、番茄(李曉旭等,2016)、棉花(呂有軍等,2017)、烏拉爾圖小麥(單強(qiáng)強(qiáng)等,2019)、玉米(鄭玲等,2019)、簸箕柳(Wang et al.,2019)、蒺藜苜蓿和鷹嘴豆(張洪濤,2018)中WOX蛋白成員的研究中發(fā)現(xiàn),其WOX蛋白成員的氨基酸數(shù)量和分子量差異很大。其中,對(duì)棉花(呂有軍等,2017)、蘋果(王楚堃等,2019)和烏拉爾圖小麥(單強(qiáng)強(qiáng)等,2019)WOX蛋白成員進(jìn)行穩(wěn)定性和親疏水性預(yù)測(cè)分析發(fā)現(xiàn)其WOX蛋白成員均為不穩(wěn)定的親水蛋白,張洪濤等(2018)對(duì)蒺藜苜蓿和鷹嘴豆WOX蛋白成員的親疏水性分析發(fā)現(xiàn)蒺藜苜蓿和鷹嘴豆中WOX蛋白成員均為親水蛋白。本研究結(jié)果表明,中華獼猴桃的16個(gè)WOX蛋白成員氨基酸數(shù)量為138~371,分子量為16.222~42.185 kD的不穩(wěn)定親水蛋白。結(jié)合前人以及本研究中WOX蛋白質(zhì)理化特性分析結(jié)果可以推測(cè)植物大部分WOX蛋白成員的氨基酸數(shù)量和分子量差異性很大且為不穩(wěn)定的親水蛋白。對(duì)烏拉爾圖小麥(單強(qiáng)強(qiáng)等,2019)、玉米(鄭玲等,2019)、鷹嘴豆(張洪濤,2018)、簸箕柳(Wang et al.,2019)、谷子和狗尾草(李雪垠等,2019)基因組中的WOX蛋白理化特性分析結(jié)果發(fā)現(xiàn),其大部分WOX蛋白成員為定位在細(xì)胞核的堿性蛋白,而在蒺藜苜蓿(張洪濤,2018)、棉花(呂有軍等,2017)、蘋果(王楚堃等,2019)、星油藤(郭亞靜等,2020)中WOX蛋白成員以酸性蛋白為主,在番茄的WOX蛋白成員中酸性蛋白和堿性蛋白各占一半。棉花WOX轉(zhuǎn)錄因子的亞細(xì)胞預(yù)測(cè)結(jié)果顯示,棉花中的WOX蛋白成員大部分定位于細(xì)胞質(zhì)和微體,少部分定位于細(xì)胞核。中華獼猴桃的WOX蛋白成員大多為定位于細(xì)胞核的堿性蛋白,極少數(shù)為定位于細(xì)胞核的酸性蛋白。此結(jié)果與烏拉爾圖小麥、玉米鷹嘴豆、簸箕柳、谷子中WOX蛋白成員的分析結(jié)果有相似之處,但由于缺少具體實(shí)驗(yàn)數(shù)據(jù)的支撐,暫不能得出WOX蛋白成員中亞細(xì)胞定位與酸堿性的關(guān)系結(jié)論,但可以推測(cè)不同的物種中WOX蛋白成員由于發(fā)揮的功能不同造成其WOX蛋白成員的亞細(xì)胞定位和酸堿性也會(huì)有所差異。由氨基酸組成來(lái)看,中華獼猴桃中16個(gè)WOX蛋白成員主要以絲氨酸為主,其次為亮氨酸、蘇氨酸和谷氨酰胺,而在其他物種中對(duì)WOX蛋白成員的氨基酸組成分析尚未見報(bào)道。

    保守結(jié)構(gòu)域分析表明,中華獼猴桃的16個(gè)WOX轉(zhuǎn)錄因子均屬于Homodomain superfamily,此結(jié)果與玉米(鄭玲等,2019)、番茄(李曉旭等,2016)、葡萄(王鵬飛等,2018)等研究結(jié)果一致,在玉米(鄭玲等,2019)、番茄(李曉旭等,2016)、葡萄(王鵬飛等,2018)等研究中,所有的WOX蛋白成員都含有homeobox這一結(jié)構(gòu)域。信號(hào)肽和跨膜結(jié)構(gòu)域顯示,在中華獼猴桃中,大多數(shù)成員無(wú)跨膜區(qū)和信號(hào)肽,僅Achn362451是分泌蛋白,且僅Achn362451和Achn141001具有跨膜區(qū),為跨膜蛋白。而在陸地棉花(呂有軍等,2017)中,所有的WOX都不含跨膜區(qū)和信號(hào)肽。故推測(cè)大多數(shù)WOX蛋白成員為非跨膜蛋白或非分泌蛋白。

    中華獼猴桃大多數(shù)WOX蛋白二級(jí)結(jié)構(gòu)主要以無(wú)規(guī)則卷曲為主,其次為α-螺旋和延伸鏈,β-轉(zhuǎn)角含量最少。此結(jié)果與陸地棉中WOX轉(zhuǎn)錄因子的二級(jí)結(jié)構(gòu)含量預(yù)測(cè)結(jié)果相似(呂有軍等,2017;Hao et al.,2019)。保守基序分析結(jié)果顯示,中華獼猴桃中所有的WOX蛋白成員均含有motif1。李曉旭等(2016)對(duì)番茄的WOX基因家族研究中發(fā)現(xiàn)所有成員都含有motif1,并指出motif1是由60個(gè)氨基酸殘基組成的保守的同源異型結(jié)構(gòu)域homeodomain。王鵬飛等(2018)在葡萄的WOX基因家族研究中發(fā)現(xiàn)所有WOX蛋白成員均含有motif1和motif2,并進(jìn)一步分析指出這兩個(gè)結(jié)構(gòu)域序列是包含WOX蛋白的homeobox結(jié)構(gòu)域的序列。王楚堃等(2019)在蘋果WOX基因家族研究中發(fā)現(xiàn),其WOX所有蛋白成員均含motif1,且motif1結(jié)構(gòu)域包含WOX家族homeobox的結(jié)構(gòu)域序列。李晨曦等(2016)在毛竹WOX基因家族研究中發(fā)現(xiàn),所有WOX蛋白成員都含有motif1和motif3,且指出motif1和motif3共同組成了homeodamin。本研究的中華獼猴桃WOX蛋白基序分析結(jié)果與以上這些植物基因組中WOX蛋白的保守基序分析結(jié)果相似,故推測(cè)中華獼猴桃中motif1在所有基序中最為保守,可能包含有WOX蛋白的homeobox結(jié)構(gòu)域序列。

    系統(tǒng)進(jìn)化研究揭示了中華獼猴桃的16個(gè)WOX蛋白成員被分為遠(yuǎn)古進(jìn)化支、中間進(jìn)化支和現(xiàn)代進(jìn)化支。此結(jié)果與番茄、棉花、玉米、蘋果、葡萄、烏拉爾圖小麥、星油藤、歐洲油菜(Li et al.,2019)、人參(Liu et al.,2019)、玫瑰和野薔薇、蒺藜苜蓿和鷹嘴豆中構(gòu)建的WOX進(jìn)化樹的分支相同,但與毛竹構(gòu)建的WOX進(jìn)化樹不同。李晨曦等(2016)在構(gòu)建的毛竹WOX進(jìn)化樹中發(fā)現(xiàn),其WOX蛋白成員分為現(xiàn)代/WUS進(jìn)化分支和古老進(jìn)化支,并指出由于缺少與10個(gè)擬南芥對(duì)應(yīng)同源基因,且不含中間進(jìn)化支,故與一般種子植物含有的3個(gè)分支不同。在中華獼猴桃與其他15個(gè)物種構(gòu)建的進(jìn)化樹中,中華獼猴與菠蘿的親緣關(guān)系最近,與番茄的親緣關(guān)系最遠(yuǎn)。中華獼猴桃授粉后20、120、127 d的轉(zhuǎn)錄組表達(dá)分析結(jié)果表明,其Achn131681、Achn145561、Achn336591和Achn362451基因表達(dá)水平顯著提高,Achn131681、Achn336591和Achn362451基因表達(dá)水平呈先上升后下降的趨勢(shì),而Achn145561基因表達(dá)水平一直上升,除了這4個(gè)WOX基因外,中華獼猴桃其他12個(gè)WOX基因在其生長(zhǎng)發(fā)育及其成熟過(guò)程中表達(dá)含量很低甚至不表達(dá)。本研究為進(jìn)一步研究中華獼猴桃WOX轉(zhuǎn)錄因子的分子功能提供了一定的理論依據(jù)。

    參考文獻(xiàn):

    CHEN Q,LIU J,CHEN DY,et al.,2018. Genome-wide identification,phylogenetic evolution,and gene expression pattern of EIN3/EIL transcription factors in Actinidia chinensis [J]. Chin J Appl Environ Biol,24(2):315-321.[陳強(qiáng),劉建,陳丹陽(yáng),等,2018. 中華獼猴桃EIN3/EIL轉(zhuǎn)錄因子家族的成員鑒定、系統(tǒng)進(jìn)化和基因表達(dá)模式[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào),24(2):315-321.]

    DEVEAUX Y,TOFFANO-NIOCHE C,CLAISSE G,et al.,2008. Genes of the most conserved WOX clade in plants affect root and flower development in Arabidopsis[J]. BMC Evol Biol,8(1):291-0.

    Editorial Department of Hubei Forestry Science and Technology. 2007. Nutritional value and health function of kiwifruit[J]. Hubei For Sci Technol,(2):9.[《湖北林業(yè)科技》編輯部,2007. 獼猴桃的營(yíng)養(yǎng)價(jià)值與保健功能[J]. 湖北林業(yè)科技,(2):9.]

    GAO L,SUN YM,SHAO TM,et al.,2015. Studies on WUSCHEL-related homebox(WOX)protein family[J]. Biotechnol Bull,31(5):7-12.[高麗,孫祎敏,邵鐵梅,等,2015. 植物WUSCHEL-related homeobox (WOX)家族研究進(jìn)展[J]. 生物技術(shù)通報(bào),31(5):7-12.]

    GUO YJ,F(xiàn)U QT,HU XD,et al.,2020. Identification and expression analysis of WOX family genes in sacha inchi (Plukenetia volubilis)[J]. Mol Plant Breed,18(12):3907-3915.[郭亞靜,付乾堂,胡曉笛,等,2020. 星油藤WOX基因家族成員的鑒定及其表達(dá)譜分析[J]. 分子植物育種,18(12):3907-3915.]

    HAECKER A,2004. Expression dynamics of WOX genes mark cell fate decisions during early embryonic patterning in Arabidopsis thaliana[J]. Development,131(3):657-668.

    HAO QN,ZHANG L,YANG YY,et al.,2019. Genome-wide analysis of the WOX gene family and function exploration of GmWOX18 in soybean[J]. Plants (Basel,Switzerland),8(7).

    HUANG SX,DING J,DENG DT,et al.,2013. Draft genome of the kiwifruit Actinidia chinensis[J]. Nat Comm,4(4): 2640.

    LI CX,CHEN L,JIANG T,et al.,2016. Identification and expression analysis of WUSCHEL-related homeobox(WOX)gene family in moso bamboo[J]. Mol Plant Breed,14(8):2003-2011.[李晨曦,陳磊,江濤,等,2016. 毛竹WUSCHEL-related homeobox(WOX)基因家族鑒定及表達(dá)分析[J]. 分子植物育種,14(8):2003-2011. ]

    LI JB,2017. Functional analysis of WOX and HSP17.8 in poplar development and stress resistance[D]. Beijing:Chinese Academy of Forestry.[李建波,2017. 楊樹WOX和HSP17.8基因在植物生長(zhǎng)發(fā)育和抗逆脅迫中的調(diào)控作用研究[D]. 北京:中國(guó)林業(yè)科學(xué)研究院.]

    LI MD,WANG RH,LIU ZY,et al.,2019. Genome-wide identification and analysis of the WUSCHEL-related homeobox (WOX)gene family in allotetraploid Brassica napus reveals changes in WOX genes during polyploidization[J]. BMC Genomics,20(1).

    LI XX,LIU C,LI W,et al.,2016. Genome-wide identification,phylogenetic analysis and expression profiling of the WOX family genes in Solanum lycopersicum[J]. Hereditas,38(5):444-460.[李曉旭,劉成,李偉,等,2016. 番茄WOX轉(zhuǎn)錄因子家族的鑒定及其進(jìn)化、表達(dá)分析[J]. 遺傳,38(5):444-460.

    LI XY,LI J,CAI MM,et al.,2020. Identification and evolution of the WUSCHEL related homeobox protein family in Bambusoideae[J]. Biomolecules,10:739.

    LI XY,WU YM,ZHANG FJ,2019. Genome-wide identification and analysis of WOX transcription factor genes in Setaria italic and Setaria viridis[J]. J Shanxi Agric Sci,47(6):957-962.[李雪垠,武懿茂,張鳳潔,2019. 谷子和狗尾草WOX轉(zhuǎn)錄因子基因的全基因組鑒定與分析[J]. 山西農(nóng)業(yè)科學(xué),47(6):957-962. ]

    LI ZQ,2018. Genome-wideanalysis of WOX gene family and functional research of TaWUS in wheat[D]. Shandong:Shandong Agricultural University.[李仲卿,2018. 小麥WOX基因家族的全基因組分析及TaWUS的功能研究[D]. 山東:山東農(nóng)業(yè)大學(xué).]

    LIAN GB,DING ZW,WANG Q,et al.,2014. Origins and evolution of WUSCHEL-related homeobox protein family in plant Kingdom[J]. Sci World J,(1):1-12.

    LIE C,KELSOM C,WU XL,2012. WOX2 and STIMPY-LIKE/WOX8 promote cotyledon boundary formation in Arabidopsis[J]. Plant J,72(4):674-682.

    LIU J,JIANG C,CHEN T,et al.,2019. Identification and 3D gene expression patterns of WUSCEHEL-related homeobox (WOX)genes from Panax ginseng[J]. Plant Physiol Biochem,143(2019):257-264.

    LIU JC,SHENG LH,XU YQ,et al.,2014. WOX11 and 12 are involved in the first-step cell fate transition during de novo root organogenesis in Arabidopsis[J]. Plant Cell,26(3):1081-1093.

    L YJ,YANG WJ,ZHAO LJ,et al.,2017. Genome-wide identification and analysis of WOX family genes in Gossypium hirsutum L.[J]. Acta Agric Boreal-Sin,32(2):21-31.[呂有軍,楊衛(wèi)軍,趙蘭杰,等,2017. 棉花WOX轉(zhuǎn)錄因子家族基因的全基因組鑒定與分析[J]. 華北農(nóng)學(xué)報(bào),32(2): 21-31.]

    MCATEE PA,RICHARDSON AC,NIEUWENHUIZEN NJ,et al.,2015. The hybrid non-ethylene and ethylene ripening response in kiwifruit (Actinidia chinensis)is associated with differential regulation of MADS-box transcription factors[J]. BMC Plant Biol,15(1):304.

    PARK SO,2005. The PRETTY FEW SEEDS2 gene encodes an Arabidopsis homeodomain protein that regulates ovule development[J]. Development,132(4):841-849.

    RAMKUMAR TR,KANCHAN M,UPADHYAY SK,et al.,2018. Identification and characterization

    of WUSCHEL-related homeobox (WOX)gene family in economically important orchid species Phalaenopsis equestris and Dendrobium catenatum[J]. Plant Gene,14:37-45.

    SHAN QQ,YANG WL,LI YJ,et al.,2019. Cloning and exoression analysis of WOX genes in Triticum urartu[J]. J Tritic Crops,39(8):903-911.[單強(qiáng)強(qiáng),陽(yáng)文龍,李英潔,等,2019. 烏拉爾圖小麥WOX基因的克隆與表達(dá)分析[J]. 麥類作物學(xué)報(bào),39(8):903-911.]

    SHIMIZU R,JI J,KELSEY E,et al.,2008. Tissue specificity and evolution of meristematic WOX3 function[J]. Plant Physiol,149(2):841-850.

    SONG PJ,TANG Y,YU XX,1984. Utilization of ascorbic acid in concentrated kiwifruit juice in human body[J]. Acta Nutr Sin,6(1):35-40.[宋圃菊,唐儀,俞錫璇,1984. 獼猴桃濃縮果汁中抗壞血酸在人體的利用[J]. 營(yíng)養(yǎng)學(xué)報(bào),6(1): 35-40.]

    SUDHIR K,GLEN S,MICHAEL L,et al.,2018. MEGA X:Molecular evolutionary genetics analysis across computing platforms[J]. Mol Biol Evol,(6):6.

    TANG W,ZHENG Y,DONG J,et al.,2016. Comprehensive transcriptome profiling reveals long noncoding RNA expression and alternative splicing regulation during fruit development and ripening in kiwifruit (Actinidia chinensis)[J]. Front Plant Sci,7:335.

    TANG YH,BAO XX,WANG J,et al.,2019. Genome-wide identification and expression analysis of the WOX gene family in physicnut[J]. Mol Plant Breed,17(4):1154-1162.[唐躍輝,包欣欣,王健,等,2019. 小桐子WOX基因家族全基因組鑒定與表達(dá)分析[J]. 分子植物育種,17(4):1154-1162.]

    WANG CK,HAN PM,WANG YM,et al.,2019. Genome-wide identification and analysis of apple WUSCHEL-related homeobox(WOX)family genes[J]. Acta Hortic Sin,46(6):1021-1032.[王楚堃,韓朋良,王詠梅,等,2019. 蘋果WUSCHEL-related homeobox(WOX)家族基因的鑒定與分析[J]. 園藝學(xué)報(bào),46(6):1021-1032.]

    WANG P,GUO Y,CHEN X,et al.,2019. Genome-wide identification of WOX genes and their expression patterns under different hormone and abiotic stress treatments in tea plant (Camellia sinensis)[J]. Trees,4(33):1129-1142.

    WANG PF,WU XY,YANG LY,et al.,2018. Identification and ananlysis of grape WUSCHEL-related homebox(WOX)family genes[J]. Sino-Overseas Grapevine Wine,(3):16-23.[王鵬飛,吳新穎,楊立英,等,2018. 葡萄WUSCHEL-related homeobox(WOX)家族基因鑒定與分析[J]. 中外葡萄與葡萄酒,(3):16-23.]

    WANG S,ZHANG Y,REN MX,et al.,2019. Genome-wide analysis of the WOX family reveals their involvement in stem growth of Populus trichocarpa[J]. Bull Bot Res,(4):568-577.[王爽,張洋,任夢(mèng)軒,等,2019. 全基因組分析楊樹WOX基因家族在莖部發(fā)育中的作用[J]. 植物研究,(4):568-577.]

    WANG XL,BI CW,WANG CY,et al.,2019. Genome-wide identification and characterization of WUSCHEL-related homeobox(WOX)genes in Salix suchowensis[J]. J For Res,30:1811-1822.

    WANG ZJ,HUANG SP,OU ZL,et al.,2017. The cloning and sequence analysis of promoter of WOX10 gene family in Arabidopsis thaliana[J]. J Huangshan Univ,19(3):58-61.[王占軍,黃士平,歐祖蘭,等,2017. 擬南芥WOX10基因啟動(dòng)子的克隆與序列分析[J]. 黃山學(xué)院學(xué)報(bào),19(3): 58-61.]

    WU XL,DABI T,WEIGEL D,2005. Requirement of homeobox gene STIMPY/WOX9 for Arabidopsis meristem growth and maintenance[J]. Curr Biol,15(5):436-440.

    YA BZ,TANG W,LI HW,et al.,2020. Kiwifruit(Actinidia chinensis)R1R2R3-MYB transcription factor AcMYB3R enhances drought and salinity tolerance in Arabidopsis thaliana[J]. J Integr Agric,(2):417-427.

    YU FY,JIANG LL,JIA JQ,et al.,2020. Comparative genomics analyses of bZIP family members in Actinidia chinensis[J]. Chin J Appl Environ Biol,26(2):226-235.[于豐源,蔣禮玲,賈舉慶,等,2020. 中華獼猴桃bZIP家族成員的比較基因組學(xué)分析[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào),26(2):226-235.]

    ZHANG HT,2018. Analysis of Medicago and Cicer arietinum WOX transcription factor gene family[M]. Harbin:Harbin Normal University.[張洪濤,2018. 蒺藜苜蓿和鷹嘴豆WOX轉(zhuǎn)錄因子基因家族分析[M]. 哈爾濱:哈爾濱師范大學(xué).]

    ZHANG YX,WU RH,QIN GT,et al.,2011. Over-expression of WOX1 leads to defects in meristem development and polyamine homeostasis in Arabidopsis[J]. J Integr Plant Biol,53(6):493-506.

    ZHENG L,F(xiàn)ENG MY,HAN JM,et al.,2019. Identification and analysis of WOX family genes in maize[J]. J Maize Sci,27(3):16-22.[鄭玲,馮夢(mèng)園,韓建明,2019. 玉米WOX家族基因的鑒定與分析[J]. 玉米科學(xué),27(3):16-22.]

    ZHU JH,SHI HZ,LEE BH,et al.,2004. An Arabidopsis homeodomain transcription factor gene,HOS9,mediates cold tolerance through a CBF-independent pathway[J]. Proc Natl Acad Sci USA,101(26):9873-9878.

    (責(zé)任編輯 李 莉)

    激情视频va一区二区三区| 久久久久精品人妻al黑| 日韩一本色道免费dvd| 免费在线观看黄色视频的| 制服人妻中文乱码| 少妇人妻精品综合一区二区| 免费少妇av软件| 亚洲av日韩精品久久久久久密 | 欧美精品一区二区免费开放| 国产片内射在线| 亚洲欧美一区二区三区久久| 男女午夜视频在线观看| 国产男女超爽视频在线观看| 色婷婷av一区二区三区视频| 最近手机中文字幕大全| 18禁观看日本| 亚洲,欧美,日韩| 午夜日本视频在线| 狂野欧美激情性xxxx| 99久久99久久久精品蜜桃| 欧美 日韩 精品 国产| 久久亚洲国产成人精品v| 99国产综合亚洲精品| 久久精品久久久久久噜噜老黄| 久久久精品94久久精品| 国产黄色视频一区二区在线观看| avwww免费| 黄频高清免费视频| 成人三级做爰电影| 五月天丁香电影| 一个人免费看片子| 在线看a的网站| 亚洲一级一片aⅴ在线观看| 亚洲 欧美一区二区三区| 激情视频va一区二区三区| 欧美老熟妇乱子伦牲交| 午夜激情久久久久久久| 亚洲一区中文字幕在线| 精品人妻在线不人妻| 亚洲七黄色美女视频| 免费不卡黄色视频| netflix在线观看网站| 熟妇人妻不卡中文字幕| 大陆偷拍与自拍| 99香蕉大伊视频| 综合色丁香网| 97精品久久久久久久久久精品| 搡老岳熟女国产| 精品国产乱码久久久久久男人| 国产极品粉嫩免费观看在线| 精品久久久精品久久久| 亚洲国产毛片av蜜桃av| 国产精品国产av在线观看| 无遮挡黄片免费观看| 亚洲av电影在线观看一区二区三区| 免费观看a级毛片全部| 欧美人与性动交α欧美精品济南到| 久久久久精品久久久久真实原创| 中国三级夫妇交换| 国产亚洲最大av| 久久久久久久久免费视频了| 9191精品国产免费久久| 国产老妇伦熟女老妇高清| 国产精品一区二区精品视频观看| 亚洲精品在线美女| 国产成人系列免费观看| 午夜福利网站1000一区二区三区| 啦啦啦在线观看免费高清www| 人妻人人澡人人爽人人| tube8黄色片| 亚洲第一区二区三区不卡| 十分钟在线观看高清视频www| 午夜91福利影院| 久久久国产欧美日韩av| 久久免费观看电影| 精品一区二区三区av网在线观看 | 赤兔流量卡办理| 男女边摸边吃奶| √禁漫天堂资源中文www| 国产精品免费视频内射| 午夜福利乱码中文字幕| 91国产中文字幕| 成人18禁高潮啪啪吃奶动态图| 在线免费观看不下载黄p国产| 久热这里只有精品99| 亚洲精品自拍成人| 精品一区二区三区av网在线观看 | 男女之事视频高清在线观看 | 日韩av免费高清视频| 老司机在亚洲福利影院| 国产探花极品一区二区| 精品国产一区二区久久| 亚洲国产精品成人久久小说| 亚洲人成网站在线观看播放| 91国产中文字幕| 日本爱情动作片www.在线观看| 久久精品久久久久久久性| 国产成人一区二区在线| 精品一区二区三区四区五区乱码 | 美女国产高潮福利片在线看| 国产伦理片在线播放av一区| av福利片在线| 一本大道久久a久久精品| 国产免费又黄又爽又色| 天天躁夜夜躁狠狠久久av| 美女扒开内裤让男人捅视频| 妹子高潮喷水视频| 欧美日韩综合久久久久久| av片东京热男人的天堂| 最近2019中文字幕mv第一页| 80岁老熟妇乱子伦牲交| 欧美精品一区二区免费开放| 国产在视频线精品| 精品一品国产午夜福利视频| 只有这里有精品99| 99久久99久久久精品蜜桃| 少妇精品久久久久久久| 街头女战士在线观看网站| 亚洲情色 制服丝袜| 晚上一个人看的免费电影| 亚洲精品美女久久av网站| 亚洲美女视频黄频| 亚洲欧美精品综合一区二区三区| 欧美人与善性xxx| 日韩大码丰满熟妇| 丁香六月欧美| 秋霞伦理黄片| 国产精品一国产av| 我要看黄色一级片免费的| 男女无遮挡免费网站观看| 日韩大码丰满熟妇| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产毛片av蜜桃av| 久久久亚洲精品成人影院| 久久久精品94久久精品| 久久久久久免费高清国产稀缺| 91精品三级在线观看| 人人妻人人澡人人看| 亚洲欧美一区二区三区久久| 欧美激情 高清一区二区三区| 91aial.com中文字幕在线观看| av在线播放精品| 亚洲精品国产av蜜桃| 精品国产露脸久久av麻豆| 美女扒开内裤让男人捅视频| 国产亚洲一区二区精品| 欧美最新免费一区二区三区| 多毛熟女@视频| 久久久精品区二区三区| 香蕉丝袜av| 欧美中文综合在线视频| 免费黄网站久久成人精品| 久久女婷五月综合色啪小说| xxx大片免费视频| 一级a爱视频在线免费观看| 在线精品无人区一区二区三| 少妇人妻 视频| 不卡av一区二区三区| 1024香蕉在线观看| 日韩不卡一区二区三区视频在线| 一区二区三区乱码不卡18| 欧美在线黄色| 日韩大码丰满熟妇| 亚洲av电影在线观看一区二区三区| 9热在线视频观看99| 黄片小视频在线播放| 国产麻豆69| 亚洲欧美中文字幕日韩二区| 国产在线视频一区二区| 国产精品人妻久久久影院| 国产 精品1| 99国产精品免费福利视频| e午夜精品久久久久久久| 老司机在亚洲福利影院| 亚洲成人av在线免费| 在线观看三级黄色| 伊人久久大香线蕉亚洲五| 国产成人精品久久久久久| 国产精品 欧美亚洲| 欧美成人精品欧美一级黄| 国产高清不卡午夜福利| av一本久久久久| 亚洲精品国产区一区二| 国产伦理片在线播放av一区| 久久人妻熟女aⅴ| 亚洲精品aⅴ在线观看| 大香蕉久久成人网| 久久午夜综合久久蜜桃| 欧美精品亚洲一区二区| 精品国产一区二区三区四区第35| 女性生殖器流出的白浆| 看非洲黑人一级黄片| 久久综合国产亚洲精品| 日韩 亚洲 欧美在线| 久久精品熟女亚洲av麻豆精品| 中文字幕另类日韩欧美亚洲嫩草| 看非洲黑人一级黄片| 两个人看的免费小视频| 丁香六月天网| 只有这里有精品99| 久久精品亚洲熟妇少妇任你| videosex国产| 久久久久精品国产欧美久久久 | 国产日韩一区二区三区精品不卡| 亚洲av在线观看美女高潮| 国产精品 国内视频| 少妇猛男粗大的猛烈进出视频| 国产亚洲av片在线观看秒播厂| av网站免费在线观看视频| 看免费av毛片| 国产免费又黄又爽又色| 免费在线观看黄色视频的| 曰老女人黄片| 黄色一级大片看看| 欧美少妇被猛烈插入视频| 飞空精品影院首页| 十八禁人妻一区二区| 久久韩国三级中文字幕| 亚洲视频免费观看视频| 欧美日韩亚洲国产一区二区在线观看 | 国产成人啪精品午夜网站| 飞空精品影院首页| 亚洲人成电影观看| 国产片特级美女逼逼视频| 国产伦理片在线播放av一区| 亚洲国产欧美在线一区| 另类精品久久| 亚洲av中文av极速乱| 久久韩国三级中文字幕| 亚洲精华国产精华液的使用体验| 欧美精品一区二区免费开放| 国产精品免费视频内射| 亚洲av日韩在线播放| 人人妻人人添人人爽欧美一区卜| 日韩视频在线欧美| 在线观看免费午夜福利视频| 老司机在亚洲福利影院| 天堂中文最新版在线下载| 欧美日韩视频高清一区二区三区二| 国产深夜福利视频在线观看| 大香蕉久久成人网| 久热爱精品视频在线9| 久久女婷五月综合色啪小说| 啦啦啦中文免费视频观看日本| 午夜福利视频在线观看免费| 日本wwww免费看| 777久久人妻少妇嫩草av网站| 国产欧美亚洲国产| 久久午夜综合久久蜜桃| 精品人妻在线不人妻| 色婷婷av一区二区三区视频| 三上悠亚av全集在线观看| a级毛片在线看网站| 丝袜美足系列| www.av在线官网国产| 日韩av免费高清视频| 国产在线一区二区三区精| 中文欧美无线码| 中文字幕精品免费在线观看视频| 一本一本久久a久久精品综合妖精| 国产精品一二三区在线看| 日韩一区二区视频免费看| 69精品国产乱码久久久| 日韩一区二区三区影片| 永久免费av网站大全| 亚洲精品av麻豆狂野| 午夜免费男女啪啪视频观看| 色吧在线观看| www.自偷自拍.com| 国产精品一国产av| 精品国产超薄肉色丝袜足j| 国产女主播在线喷水免费视频网站| 久久99热这里只频精品6学生| 久久国产精品男人的天堂亚洲| 亚洲美女黄色视频免费看| www日本在线高清视频| 伦理电影大哥的女人| 久久影院123| 男女无遮挡免费网站观看| 国产成人91sexporn| 亚洲欧美日韩另类电影网站| 亚洲av欧美aⅴ国产| 国产免费一区二区三区四区乱码| 亚洲男人天堂网一区| 成人三级做爰电影| 午夜影院在线不卡| 一级毛片 在线播放| 纯流量卡能插随身wifi吗| 亚洲av中文av极速乱| 亚洲综合精品二区| 90打野战视频偷拍视频| 少妇猛男粗大的猛烈进出视频| 国产精品久久久人人做人人爽| 九九爱精品视频在线观看| 亚洲激情五月婷婷啪啪| 五月天丁香电影| 国产高清不卡午夜福利| 啦啦啦视频在线资源免费观看| 亚洲精品久久成人aⅴ小说| 99热网站在线观看| 久久久久久人人人人人| 欧美日韩一级在线毛片| 亚洲欧美一区二区三区国产| 狠狠婷婷综合久久久久久88av| 我的亚洲天堂| 在线精品无人区一区二区三| 综合色丁香网| 热99久久久久精品小说推荐| 男女高潮啪啪啪动态图| 国产av一区二区精品久久| 久久久久国产一级毛片高清牌| 三上悠亚av全集在线观看| 街头女战士在线观看网站| 侵犯人妻中文字幕一二三四区| 国产av码专区亚洲av| 免费观看性生交大片5| 久久精品aⅴ一区二区三区四区| 在线观看免费日韩欧美大片| 亚洲精品av麻豆狂野| 亚洲精品一区蜜桃| 欧美最新免费一区二区三区| 性少妇av在线| 国产精品久久久久久精品电影小说| 另类亚洲欧美激情| 在线观看免费午夜福利视频| 免费在线观看视频国产中文字幕亚洲 | 一级,二级,三级黄色视频| 国产99久久九九免费精品| 波多野结衣av一区二区av| 久久天躁狠狠躁夜夜2o2o | 久久综合国产亚洲精品| 女的被弄到高潮叫床怎么办| 99国产综合亚洲精品| 国产人伦9x9x在线观看| 另类亚洲欧美激情| 国产高清国产精品国产三级| 亚洲第一青青草原| 精品第一国产精品| 91精品三级在线观看| 久久热在线av| 99国产精品免费福利视频| 精品第一国产精品| 精品少妇黑人巨大在线播放| 国精品久久久久久国模美| 精品人妻熟女毛片av久久网站| 亚洲欧洲精品一区二区精品久久久 | 免费观看人在逋| 日韩欧美精品免费久久| 欧美激情极品国产一区二区三区| 久久99热这里只频精品6学生| 亚洲精华国产精华液的使用体验| 男女国产视频网站| 两个人免费观看高清视频| 一边亲一边摸免费视频| 美女午夜性视频免费| 一本色道久久久久久精品综合| 自线自在国产av| 操出白浆在线播放| 亚洲精品一区蜜桃| 极品人妻少妇av视频| 九九爱精品视频在线观看| 亚洲天堂av无毛| 啦啦啦啦在线视频资源| 国产精品成人在线| 18禁动态无遮挡网站| 91国产中文字幕| 青春草国产在线视频| av网站免费在线观看视频| 亚洲欧美中文字幕日韩二区| 国产女主播在线喷水免费视频网站| 无限看片的www在线观看| 你懂的网址亚洲精品在线观看| 久久韩国三级中文字幕| 看十八女毛片水多多多| 美国免费a级毛片| 欧美激情高清一区二区三区 | 纵有疾风起免费观看全集完整版| 中文乱码字字幕精品一区二区三区| 中文字幕人妻丝袜制服| av女优亚洲男人天堂| 亚洲国产精品一区二区三区在线| 最近手机中文字幕大全| av.在线天堂| 精品久久久精品久久久| 国产福利在线免费观看视频| 狂野欧美激情性bbbbbb| 十八禁人妻一区二区| 天堂俺去俺来也www色官网| 欧美日韩成人在线一区二区| 99热全是精品| 国产极品天堂在线| 99热全是精品| 超色免费av| 建设人人有责人人尽责人人享有的| 亚洲精品自拍成人| 国产日韩欧美视频二区| 亚洲色图 男人天堂 中文字幕| 国产亚洲一区二区精品| 超色免费av| 999精品在线视频| 精品国产乱码久久久久久小说| 欧美成人精品欧美一级黄| 日本色播在线视频| 国产欧美日韩一区二区三区在线| 美女主播在线视频| 新久久久久国产一级毛片| e午夜精品久久久久久久| 少妇精品久久久久久久| 国产精品 欧美亚洲| kizo精华| 国产一区二区 视频在线| 高清视频免费观看一区二区| 亚洲,一卡二卡三卡| 女人久久www免费人成看片| 精品国产国语对白av| 国产伦人伦偷精品视频| 亚洲欧美精品综合一区二区三区| 久久久亚洲精品成人影院| 日韩 亚洲 欧美在线| 国产精品国产三级专区第一集| 亚洲欧美激情在线| 欧美日韩一级在线毛片| 国产亚洲精品第一综合不卡| 成人三级做爰电影| a级片在线免费高清观看视频| 欧美中文综合在线视频| 日本91视频免费播放| 波野结衣二区三区在线| 这个男人来自地球电影免费观看 | 秋霞在线观看毛片| 十八禁网站网址无遮挡| 国产男人的电影天堂91| 免费黄色在线免费观看| 国产精品一区二区在线观看99| 精品一区在线观看国产| 国产色婷婷99| 性少妇av在线| 纯流量卡能插随身wifi吗| 91国产中文字幕| 水蜜桃什么品种好| 波野结衣二区三区在线| av国产精品久久久久影院| 亚洲男人天堂网一区| 国产成人免费无遮挡视频| 丝袜美腿诱惑在线| 国产精品一区二区精品视频观看| av在线观看视频网站免费| 久久久久久久大尺度免费视频| 国产欧美日韩综合在线一区二区| 日日啪夜夜爽| 国产亚洲午夜精品一区二区久久| 亚洲,欧美精品.| 看免费av毛片| 免费女性裸体啪啪无遮挡网站| 在线 av 中文字幕| 亚洲人成电影观看| 黑人巨大精品欧美一区二区蜜桃| 成人国产av品久久久| 精品一区二区免费观看| 免费在线观看完整版高清| 亚洲色图综合在线观看| 丝袜美腿诱惑在线| 少妇 在线观看| 亚洲美女黄色视频免费看| 亚洲美女搞黄在线观看| 99热网站在线观看| 青春草视频在线免费观看| 丰满乱子伦码专区| 国产在线免费精品| 五月天丁香电影| 一区二区三区精品91| 日韩中文字幕欧美一区二区 | 我的亚洲天堂| 欧美日韩亚洲综合一区二区三区_| 国产精品嫩草影院av在线观看| 电影成人av| 婷婷色av中文字幕| 精品一区二区三区四区五区乱码 | 国产一卡二卡三卡精品 | 男女无遮挡免费网站观看| 国产一区二区 视频在线| 90打野战视频偷拍视频| 国产 精品1| 欧美在线一区亚洲| av女优亚洲男人天堂| videos熟女内射| 欧美黑人精品巨大| 日本vs欧美在线观看视频| 飞空精品影院首页| 国产精品无大码| 侵犯人妻中文字幕一二三四区| 制服丝袜香蕉在线| 亚洲自偷自拍图片 自拍| 香蕉丝袜av| 我要看黄色一级片免费的| 香蕉丝袜av| 亚洲欧美清纯卡通| 欧美日韩一区二区视频在线观看视频在线| 91精品国产国语对白视频| 美女高潮到喷水免费观看| 黄色毛片三级朝国网站| 黄色一级大片看看| 最近中文字幕高清免费大全6| kizo精华| 国产精品 欧美亚洲| 赤兔流量卡办理| 在线观看国产h片| 黄片小视频在线播放| 久久综合国产亚洲精品| a级片在线免费高清观看视频| 久久久精品94久久精品| 久久精品熟女亚洲av麻豆精品| 中国国产av一级| 女性生殖器流出的白浆| 97人妻天天添夜夜摸| 欧美人与性动交α欧美精品济南到| av线在线观看网站| 成人三级做爰电影| 19禁男女啪啪无遮挡网站| 2021少妇久久久久久久久久久| 午夜日本视频在线| 狠狠精品人妻久久久久久综合| 狂野欧美激情性xxxx| 一本一本久久a久久精品综合妖精| 日日爽夜夜爽网站| 美女大奶头黄色视频| 熟妇人妻不卡中文字幕| 男女高潮啪啪啪动态图| 天天添夜夜摸| 国产日韩欧美视频二区| 一级毛片 在线播放| 亚洲,欧美精品.| 亚洲精品国产av成人精品| 日韩av在线免费看完整版不卡| 女人被躁到高潮嗷嗷叫费观| 国产男女内射视频| 国产在线免费精品| 最近2019中文字幕mv第一页| 精品国产乱码久久久久久小说| 亚洲欧美色中文字幕在线| 午夜免费鲁丝| 久久99一区二区三区| 美女视频免费永久观看网站| av又黄又爽大尺度在线免费看| 尾随美女入室| 女人爽到高潮嗷嗷叫在线视频| 久久天躁狠狠躁夜夜2o2o | 国产免费又黄又爽又色| 观看美女的网站| 亚洲国产精品一区二区三区在线| 超碰97精品在线观看| 色吧在线观看| 99re6热这里在线精品视频| 欧美国产精品va在线观看不卡| 久久性视频一级片| 巨乳人妻的诱惑在线观看| 日韩成人av中文字幕在线观看| 国产有黄有色有爽视频| 久久精品久久久久久久性| av片东京热男人的天堂| 亚洲,欧美,日韩| 五月天丁香电影| 精品国产乱码久久久久久小说| 日韩 欧美 亚洲 中文字幕| 国产1区2区3区精品| 麻豆av在线久日| 伊人久久国产一区二区| 青青草视频在线视频观看| 精品国产一区二区三区久久久樱花| 国产 精品1| 精品一区二区免费观看| 一本一本久久a久久精品综合妖精| 亚洲精品国产av成人精品| 一区二区三区精品91| 国产日韩欧美在线精品| 校园人妻丝袜中文字幕| 久久精品国产亚洲av高清一级| 老汉色∧v一级毛片| av网站在线播放免费| 色94色欧美一区二区| 亚洲三区欧美一区| 国产成人欧美在线观看 | 青草久久国产| 国产xxxxx性猛交| 91精品伊人久久大香线蕉| 国产精品久久久久久人妻精品电影 | 亚洲国产欧美一区二区综合| 国产乱人偷精品视频| 久久久久人妻精品一区果冻| 国产伦人伦偷精品视频| 亚洲专区中文字幕在线 | 国产成人一区二区在线| 欧美人与性动交α欧美软件| 免费在线观看黄色视频的| 亚洲视频免费观看视频| 热99国产精品久久久久久7| 成人国产麻豆网| 国产又色又爽无遮挡免| 女人爽到高潮嗷嗷叫在线视频| 国产精品偷伦视频观看了| 亚洲国产看品久久| 精品国产国语对白av| 成人黄色视频免费在线看| 美女中出高潮动态图| 亚洲欧美一区二区三区久久| 男人操女人黄网站| 80岁老熟妇乱子伦牲交| 婷婷色综合www| 婷婷色综合大香蕉| 亚洲视频免费观看视频| 在现免费观看毛片| 日韩制服骚丝袜av| 久久天堂一区二区三区四区| 一本大道久久a久久精品| bbb黄色大片|