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

    中華獼猴桃WOX轉錄因子的生物信息學分析

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

    付春 黃春莉 江納 楊瑤君

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

    關鍵詞:中華獼猴桃,WOX,轉錄因子,基因家族,生物信息學分析

    中圖分類號:Q943.2

    文獻標識碼:A

    文章編號: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)的一種經濟樹木。中華獼猴桃原產自我國,果實富含維生素C,且維生素C在人體中的利用率高達94%(宋圃菊等,1984),一個成年人每天只攝食一個鮮獼猴桃就可滿足人體對維生素C的需求,因而被譽為“維C之王”。果實酸甜可口,除了鮮食,還可加工成各種食品和飲料。臨床試驗發(fā)現,獼猴桃汁能有效抑制麻風桿菌,從而明顯改善麻風病的臨床癥狀。除此之外,獼猴桃汁所含的一種抗突變成分谷胱甘肽,可抑制誘發(fā)致癌基因的突變,并有效阻斷致癌物質亞硝酸基嗎啉的合成,從而防止癌癥的發(fā)生(《湖北林業(yè)科技》編輯部,2007)。

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

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

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

    1 材料與方法

    1.1 材料

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

    1.2 方法

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

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

    1.2.3 磷酸化位點的預測分析 通過在線工具NetPhos 3.1 Server (http://www.cbs.dtu.dk/services/NetPhos/)對中華獼猴桃的WOX轉錄因子潛在磷酸化位點進行預測分析。

    1.2.4 二級結構、三級結構的分析預測 使用SOPMA程序(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html)對中華獼猴桃WOX轉錄因子的蛋白質二級結構進行預測分析;利用在線工具SWISS-MODEL(https://swissmodel.expasy.org/)對其蛋白質的三級結構進行預測分析。

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

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

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

    2 結果與分析

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

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

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

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

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

    2.2 保守結構域預測、信號肽的分析及跨膜域的預測分析

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

    SignalP4.1 Server信號肽預測分析結果顯示除Achn141001外,其余15個中華獼猴桃轉錄因子均無信號肽。因此,推測Achn141001為分泌蛋白,其余均為非分泌蛋白(圖4)。TMHMM Severv.2.0跨膜域分析結果顯示,除Achn362451、Achn141001有明顯的跨膜結構域外,其余14個WOX轉錄因子均無明顯的跨膜結構域。因此,推測Achn362451、Achn141001為跨膜蛋白(圖5)。

    2.3 磷酸化位點的預測分析

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

    2.4 二級結構和三級結構的預測分析

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

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

    2.5 保守基序分析

    利用在線工具MEME對中華獼猴桃的WOX家族進行保守基序分析,結果表明:motif1出現在中華獼猴桃的所有的WOX蛋白成員中,motif2出現在除Achn107771和Achn141001外的其他WOX蛋白成員中。根據不同成員中所含motif不同,將中華獼猴桃的16個WOX蛋白成員分為A、B、C三個組。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基因家族在不同時期的表達水平分析

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

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

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

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

    3 討論與結論

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

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

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

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

    參考文獻:

    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.[陳強,劉建,陳丹陽,等,2018. 中華獼猴桃EIN3/EIL轉錄因子家族的成員鑒定、系統(tǒng)進化和基因表達模式[J]. 應用與環(huán)境生物學報,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)價值與保健功能[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]. 生物技術通報,31(5):7-12.]

    GUO YJ,FU 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基因家族成員的鑒定及其表達譜分析[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)基因家族鑒定及表達分析[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基因在植物生長發(fā)育和抗逆脅迫中的調控作用研究[D]. 北京:中國林業(yè)科學研究院.]

    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轉錄因子家族的鑒定及其進化、表達分析[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轉錄因子基因的全基因組鑒定與分析[J]. 山西農業(yè)科學,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]. 山東:山東農業(yè)大學.]

    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轉錄因子家族基因的全基因組鑒定與分析[J]. 華北農學報,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.[單強強,陽文龍,李英潔,等,2019. 烏拉爾圖小麥WOX基因的克隆與表達分析[J]. 麥類作物學報,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)學報,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基因家族全基因組鑒定與表達分析[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]. 園藝學報,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.[王爽,張洋,任夢軒,等,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基因啟動子的克隆與序列分析[J]. 黃山學院學報,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家族成員的比較基因組學分析[J]. 應用與環(huán)境生物學報,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轉錄因子基因家族分析[M]. 哈爾濱:哈爾濱師范大學.]

    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,FENG MY,HAN JM,et al.,2019. Identification and analysis of WOX family genes in maize[J]. J Maize Sci,27(3):16-22.[鄭玲,馮夢園,韓建明,2019. 玉米WOX家族基因的鑒定與分析[J]. 玉米科學,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.

    (責任編輯 李 莉)

    国产精品一区二区性色av| 丰满的人妻完整版| 日韩三级伦理在线观看| 亚洲高清免费不卡视频| 亚洲欧美精品综合久久99| 成人欧美大片| 国产精品久久久久久久电影| 免费电影在线观看免费观看| 欧美激情在线99| 麻豆乱淫一区二区| 在线播放无遮挡| 悠悠久久av| 久久中文看片网| 两个人视频免费观看高清| 国产中年淑女户外野战色| 中文字幕人妻熟人妻熟丝袜美| 国产精品乱码一区二三区的特点| 午夜福利高清视频| 成人特级av手机在线观看| 日本一二三区视频观看| 中文字幕久久专区| 91久久精品电影网| 国产成人精品婷婷| 精品一区二区三区视频在线| 久久人妻av系列| 国产黄片美女视频| 国产精品久久久久久精品电影| 亚洲七黄色美女视频| 久久人人爽人人片av| 精品久久久噜噜| 日本成人三级电影网站| 亚洲激情五月婷婷啪啪| 深夜精品福利| 亚洲国产精品成人综合色| 高清毛片免费观看视频网站| 亚洲国产精品国产精品| 免费一级毛片在线播放高清视频| 国产淫片久久久久久久久| 精品免费久久久久久久清纯| 一区二区三区免费毛片| 男女视频在线观看网站免费| 插阴视频在线观看视频| 最近的中文字幕免费完整| 特大巨黑吊av在线直播| 赤兔流量卡办理| 欧美最新免费一区二区三区| 久久欧美精品欧美久久欧美| 国产男人的电影天堂91| 男女视频在线观看网站免费| 97热精品久久久久久| 日韩亚洲欧美综合| 色视频www国产| av黄色大香蕉| 欧美高清成人免费视频www| .国产精品久久| 男人舔奶头视频| 丰满人妻一区二区三区视频av| 舔av片在线| 成人无遮挡网站| 狂野欧美白嫩少妇大欣赏| 久久精品国产99精品国产亚洲性色| 免费在线观看成人毛片| 人妻久久中文字幕网| 亚洲人成网站高清观看| 免费看光身美女| 成人特级黄色片久久久久久久| 大型黄色视频在线免费观看| 日韩欧美一区二区三区在线观看| 深夜a级毛片| 久久精品91蜜桃| 国产亚洲5aaaaa淫片| 在线观看一区二区三区| 国产一区二区三区在线臀色熟女| 中国国产av一级| 日韩av在线大香蕉| 午夜老司机福利剧场| 乱系列少妇在线播放| 人妻久久中文字幕网| 国产探花在线观看一区二区| 国产午夜精品论理片| 九九在线视频观看精品| 亚洲一区二区三区色噜噜| 亚洲乱码一区二区免费版| 九九久久精品国产亚洲av麻豆| 色综合站精品国产| 一级毛片我不卡| 最近2019中文字幕mv第一页| 国产精品电影一区二区三区| 日本五十路高清| 欧美日韩乱码在线| 免费看日本二区| 人妻制服诱惑在线中文字幕| 久久久久久久久久黄片| 99热6这里只有精品| 国产成人午夜福利电影在线观看| 久久精品久久久久久噜噜老黄 | 久久久精品94久久精品| 久久久久性生活片| 白带黄色成豆腐渣| 又粗又硬又长又爽又黄的视频 | 高清午夜精品一区二区三区 | 国产伦理片在线播放av一区 | 久久人人爽人人片av| 欧美日韩国产亚洲二区| 婷婷六月久久综合丁香| 简卡轻食公司| 国产精品一区二区三区四区免费观看| 中文精品一卡2卡3卡4更新| 国产老妇伦熟女老妇高清| 啦啦啦韩国在线观看视频| 天天躁夜夜躁狠狠久久av| 国产精品精品国产色婷婷| 欧美3d第一页| 中文亚洲av片在线观看爽| АⅤ资源中文在线天堂| 一级av片app| 成年版毛片免费区| 九色成人免费人妻av| 韩国av在线不卡| 欧美人与善性xxx| 欧美日韩精品成人综合77777| 免费人成视频x8x8入口观看| 一级毛片aaaaaa免费看小| 国产精品久久久久久久电影| 五月伊人婷婷丁香| 99国产极品粉嫩在线观看| avwww免费| 国产精品,欧美在线| 舔av片在线| 男女视频在线观看网站免费| 久久精品国产亚洲网站| 国产久久久一区二区三区| kizo精华| 亚洲第一区二区三区不卡| 欧美日韩精品成人综合77777| a级毛片免费高清观看在线播放| 18禁裸乳无遮挡免费网站照片| 桃色一区二区三区在线观看| 久久久久久久亚洲中文字幕| 天天一区二区日本电影三级| 国产乱人视频| 免费不卡的大黄色大毛片视频在线观看 | 国产精品一区二区三区四区免费观看| 床上黄色一级片| 女人被狂操c到高潮| 岛国毛片在线播放| 国产成人精品婷婷| 精品99又大又爽又粗少妇毛片| 欧美日韩一区二区视频在线观看视频在线 | 性色avwww在线观看| 欧美最黄视频在线播放免费| 少妇被粗大猛烈的视频| 午夜爱爱视频在线播放| .国产精品久久| 床上黄色一级片| 夜夜夜夜夜久久久久| www.av在线官网国产| 91麻豆精品激情在线观看国产| 乱系列少妇在线播放| 一个人看视频在线观看www免费| 搡女人真爽免费视频火全软件| 禁无遮挡网站| 两性午夜刺激爽爽歪歪视频在线观看| 日韩制服骚丝袜av| 桃色一区二区三区在线观看| 日本-黄色视频高清免费观看| 久久久久久大精品| 精品久久久久久久久久久久久| 日本免费一区二区三区高清不卡| 夜夜爽天天搞| 国产成人午夜福利电影在线观看| 免费av毛片视频| 久久精品夜色国产| 一级黄色大片毛片| 成人漫画全彩无遮挡| 欧美高清性xxxxhd video| 国产精品一区二区性色av| 亚洲av一区综合| 国产乱人偷精品视频| 一进一出抽搐gif免费好疼| 国产黄色小视频在线观看| 亚洲国产精品国产精品| 自拍偷自拍亚洲精品老妇| 中国美女看黄片| 国产av麻豆久久久久久久| 国产 一区精品| 午夜福利视频1000在线观看| 亚洲欧美日韩东京热| 五月玫瑰六月丁香| 青春草国产在线视频 | 一个人免费在线观看电影| 国产精品永久免费网站| 99国产极品粉嫩在线观看| 国产精品久久视频播放| 内射极品少妇av片p| 国产私拍福利视频在线观看| 国产精品蜜桃在线观看 | 国产亚洲精品久久久久久毛片| 99久久中文字幕三级久久日本| 一区二区三区免费毛片| 免费av观看视频| 最新中文字幕久久久久| 可以在线观看的亚洲视频| 亚洲欧美日韩东京热| 又粗又硬又长又爽又黄的视频 | 亚洲精品久久久久久婷婷小说 | 美女 人体艺术 gogo| 亚洲经典国产精华液单| 欧美性感艳星| 日本一本二区三区精品| av在线观看视频网站免费| 99riav亚洲国产免费| 亚洲国产日韩欧美精品在线观看| 亚洲av成人av| 中文在线观看免费www的网站| 99国产极品粉嫩在线观看| 在线观看午夜福利视频| 亚洲无线观看免费| 国产一区二区在线av高清观看| 精品日产1卡2卡| 人妻夜夜爽99麻豆av| 日本爱情动作片www.在线观看| 日韩视频在线欧美| 国产精品久久久久久精品电影| 麻豆乱淫一区二区| 亚洲一区二区三区色噜噜| 91在线精品国自产拍蜜月| 国产探花在线观看一区二区| 国产精品,欧美在线| 国产单亲对白刺激| 亚洲精品久久久久久婷婷小说 | a级毛色黄片| 精品不卡国产一区二区三区| 午夜亚洲福利在线播放| 精品久久久久久久人妻蜜臀av| 久久久成人免费电影| 青春草国产在线视频 | 国产精品久久视频播放| 免费看光身美女| 国产一区二区在线av高清观看| 亚洲精品乱码久久久v下载方式| 日韩人妻高清精品专区| 日韩亚洲欧美综合| 高清午夜精品一区二区三区 | 国产色婷婷99| 久久韩国三级中文字幕| 秋霞在线观看毛片| 成人美女网站在线观看视频| 国产三级中文精品| 夜夜爽天天搞| 亚洲天堂国产精品一区在线| 内射极品少妇av片p| 久久草成人影院| 午夜福利成人在线免费观看| 免费看光身美女| 日韩欧美精品v在线| 丝袜美腿在线中文| 国内精品久久久久精免费| 国产高清不卡午夜福利| 色5月婷婷丁香| 日韩一本色道免费dvd| 亚洲欧美成人综合另类久久久 | 夜夜看夜夜爽夜夜摸| 亚洲精品粉嫩美女一区| 国产精品蜜桃在线观看 | 亚洲国产精品sss在线观看| 欧美一区二区亚洲| 亚洲一区高清亚洲精品| 搞女人的毛片| www日本黄色视频网| 搡老妇女老女人老熟妇| 精品久久久久久成人av| 小说图片视频综合网站| 国产亚洲91精品色在线| 秋霞在线观看毛片| 日韩国内少妇激情av| 婷婷六月久久综合丁香| 久久久色成人| 成人一区二区视频在线观看| 久久精品国产自在天天线| 午夜激情欧美在线| 大型黄色视频在线免费观看| 人妻少妇偷人精品九色| 成人毛片60女人毛片免费| 欧美成人一区二区免费高清观看| 亚洲经典国产精华液单| 99久久无色码亚洲精品果冻| 欧美激情国产日韩精品一区| 乱人视频在线观看| 美女被艹到高潮喷水动态| 亚洲欧洲国产日韩| 精品日产1卡2卡| 免费看a级黄色片| 1024手机看黄色片| 欧美三级亚洲精品| 高清毛片免费观看视频网站| 美女脱内裤让男人舔精品视频 | 成人午夜精彩视频在线观看| 床上黄色一级片| 最近最新中文字幕大全电影3| 日韩精品有码人妻一区| 性欧美人与动物交配| 亚洲精品亚洲一区二区| 午夜激情福利司机影院| 国内少妇人妻偷人精品xxx网站| 99久久精品国产国产毛片| .国产精品久久| 欧美zozozo另类| 91狼人影院| 成年女人永久免费观看视频| 美女被艹到高潮喷水动态| 国产成年人精品一区二区| 亚洲第一电影网av| 美女国产视频在线观看| 中文在线观看免费www的网站| 一区福利在线观看| 国产91av在线免费观看| eeuss影院久久| 欧美丝袜亚洲另类| 欧美+亚洲+日韩+国产| 给我免费播放毛片高清在线观看| 日本成人三级电影网站| 亚洲五月天丁香| 亚洲在久久综合| 18禁在线播放成人免费| 亚洲,欧美,日韩| 国产视频首页在线观看| 久久久欧美国产精品| 亚洲一级一片aⅴ在线观看| 国产乱人偷精品视频| 在线观看一区二区三区| 日本三级黄在线观看| 在线观看一区二区三区| 九九爱精品视频在线观看| 久久精品国产自在天天线| 91精品国产九色| 久久99精品国语久久久| 欧美xxxx性猛交bbbb| 成人毛片60女人毛片免费| 久久午夜亚洲精品久久| 国内精品美女久久久久久| 国产亚洲av嫩草精品影院| 在线天堂最新版资源| 日韩欧美精品免费久久| 亚洲国产高清在线一区二区三| 最后的刺客免费高清国语| 婷婷精品国产亚洲av| 久久6这里有精品| 特级一级黄色大片| 国产黄片视频在线免费观看| 婷婷亚洲欧美| av在线播放精品| 国产精品99久久久久久久久| 一个人免费在线观看电影| 欧美+日韩+精品| 成人漫画全彩无遮挡| 亚洲成av人片在线播放无| 亚洲最大成人av| 午夜视频国产福利| 欧美性猛交╳xxx乱大交人| 国产亚洲精品av在线| 最近视频中文字幕2019在线8| 亚洲av中文字字幕乱码综合| 国产一区二区三区在线臀色熟女| 中国美白少妇内射xxxbb| 听说在线观看完整版免费高清| kizo精华| 嫩草影院新地址| 精品少妇黑人巨大在线播放 | 晚上一个人看的免费电影| 少妇人妻一区二区三区视频| 久久久国产成人精品二区| 久久人人精品亚洲av| 日本黄色片子视频| 久久精品国产自在天天线| 精品久久久久久成人av| 日日撸夜夜添| 99热这里只有是精品50| 久久精品影院6| 欧美最黄视频在线播放免费| 亚洲无线在线观看| 国产伦一二天堂av在线观看| 搡老妇女老女人老熟妇| 三级经典国产精品| 美女cb高潮喷水在线观看| 国产午夜福利久久久久久| 亚洲最大成人手机在线| 日韩人妻高清精品专区| 日本黄色片子视频| 色噜噜av男人的天堂激情| 欧美高清性xxxxhd video| 久久精品国产亚洲av香蕉五月| 在线播放国产精品三级| 国产一区亚洲一区在线观看| 免费人成在线观看视频色| 国产男人的电影天堂91| 久久久久久久久中文| 国产精品精品国产色婷婷| 欧美+日韩+精品| 两个人视频免费观看高清| 99久久九九国产精品国产免费| 性插视频无遮挡在线免费观看| 欧美日韩综合久久久久久| 日韩欧美精品免费久久| 国产精品一区二区性色av| 亚洲成人久久爱视频| 亚洲精品久久久久久婷婷小说 | 国产亚洲av片在线观看秒播厂 | 91久久精品电影网| 国产日本99.免费观看| 只有这里有精品99| 特大巨黑吊av在线直播| 99热网站在线观看| 啦啦啦韩国在线观看视频| 欧美性猛交黑人性爽| 日韩成人伦理影院| 亚洲av熟女| 亚洲经典国产精华液单| 一本久久精品| 日韩欧美 国产精品| 久99久视频精品免费| 亚洲精品日韩在线中文字幕 | av免费在线看不卡| 亚洲人成网站高清观看| 国产一区二区激情短视频| 亚洲欧美成人精品一区二区| 日韩成人av中文字幕在线观看| 国产探花极品一区二区| 97热精品久久久久久| 精品一区二区三区人妻视频| 麻豆乱淫一区二区| avwww免费| 国产大屁股一区二区在线视频| 长腿黑丝高跟| 亚洲无线观看免费| 性插视频无遮挡在线免费观看| 日本黄色视频三级网站网址| 午夜爱爱视频在线播放| 少妇人妻一区二区三区视频| 精品人妻一区二区三区麻豆| 亚洲av成人av| 国产伦在线观看视频一区| 91在线精品国自产拍蜜月| 欧美最新免费一区二区三区| 麻豆成人午夜福利视频| 黄片无遮挡物在线观看| 看片在线看免费视频| 一级二级三级毛片免费看| 两性午夜刺激爽爽歪歪视频在线观看| 一边摸一边抽搐一进一小说| 久久热精品热| 欧美色欧美亚洲另类二区| 精品久久久久久久久亚洲| 欧洲精品卡2卡3卡4卡5卡区| 我的女老师完整版在线观看| 国产成人精品一,二区 | 网址你懂的国产日韩在线| 禁无遮挡网站| 国产极品天堂在线| 欧美高清性xxxxhd video| 国产成年人精品一区二区| 看黄色毛片网站| 一个人免费在线观看电影| 欧美最新免费一区二区三区| 在线观看免费视频日本深夜| 亚洲av一区综合| 婷婷亚洲欧美| 好男人在线观看高清免费视频| 2022亚洲国产成人精品| 亚洲精品乱码久久久v下载方式| 国产一区亚洲一区在线观看| 在线播放无遮挡| 欧美在线一区亚洲| 国产一级毛片七仙女欲春2| 久久人人爽人人片av| 欧美bdsm另类| 午夜精品一区二区三区免费看| 亚洲最大成人av| 99久久九九国产精品国产免费| 国产亚洲av嫩草精品影院| 亚洲av中文av极速乱| 嘟嘟电影网在线观看| 免费电影在线观看免费观看| 日韩欧美在线乱码| 91久久精品国产一区二区成人| 亚洲中文字幕日韩| 久久久久性生活片| 日韩欧美一区二区三区在线观看| 天堂av国产一区二区熟女人妻| 国产黄片视频在线免费观看| 天堂中文最新版在线下载 | 22中文网久久字幕| 国产精品伦人一区二区| 亚洲av免费在线观看| 少妇的逼水好多| 丝袜美腿在线中文| 边亲边吃奶的免费视频| 极品教师在线视频| 爱豆传媒免费全集在线观看| 欧美日韩国产亚洲二区| 久久久久九九精品影院| 亚洲精品乱码久久久久久按摩| 成熟少妇高潮喷水视频| 99在线人妻在线中文字幕| 97超碰精品成人国产| 国产视频首页在线观看| 日韩在线高清观看一区二区三区| 在线播放无遮挡| 男女边吃奶边做爰视频| 亚洲人成网站在线观看播放| 青青草视频在线视频观看| 国产精品99久久久久久久久| 狠狠狠狠99中文字幕| 人人妻人人澡人人爽人人夜夜 | 日韩一区二区视频免费看| 亚洲熟妇中文字幕五十中出| 久久人妻av系列| 97超视频在线观看视频| 男女那种视频在线观看| 99热网站在线观看| 久久久久久久久久久丰满| 日韩制服骚丝袜av| 国产大屁股一区二区在线视频| 亚洲自偷自拍三级| 天堂网av新在线| 亚洲av男天堂| 亚洲第一电影网av| 亚洲综合色惰| 真实男女啪啪啪动态图| 美女国产视频在线观看| 最近最新中文字幕大全电影3| 久久久久久久久久黄片| 国产 一区 欧美 日韩| 国产精品一及| av黄色大香蕉| 国产成年人精品一区二区| 天堂网av新在线| 乱人视频在线观看| 欧美区成人在线视频| 一级黄片播放器| 高清日韩中文字幕在线| 国产精品不卡视频一区二区| 一区二区三区免费毛片| 国产亚洲5aaaaa淫片| 国产乱人偷精品视频| 乱人视频在线观看| 亚洲第一区二区三区不卡| 嘟嘟电影网在线观看| 国产成年人精品一区二区| 国产伦一二天堂av在线观看| 欧美成人一区二区免费高清观看| 精品人妻偷拍中文字幕| 欧美高清性xxxxhd video| 亚洲无线在线观看| 黄色日韩在线| 非洲黑人性xxxx精品又粗又长| 99热网站在线观看| 伦精品一区二区三区| 美女cb高潮喷水在线观看| 亚洲乱码一区二区免费版| 久久人妻av系列| 免费看美女性在线毛片视频| 欧美三级亚洲精品| 亚洲av不卡在线观看| 麻豆国产97在线/欧美| 国内精品久久久久精免费| 麻豆国产97在线/欧美| 91精品一卡2卡3卡4卡| 日韩成人av中文字幕在线观看| 看十八女毛片水多多多| 国产亚洲精品久久久com| 男的添女的下面高潮视频| 国产精品久久久久久久久免| 久久久久久久久久黄片| 亚洲无线观看免费| 日韩欧美精品免费久久| 国产乱人视频| 老司机福利观看| 亚洲18禁久久av| 欧美高清成人免费视频www| 天美传媒精品一区二区| 美女脱内裤让男人舔精品视频 | 少妇的逼水好多| 亚洲欧美日韩卡通动漫| 亚洲中文字幕日韩| 午夜激情福利司机影院| 亚洲精品成人久久久久久| 亚洲人成网站高清观看| 欧美一区二区亚洲| 九九在线视频观看精品| 51国产日韩欧美| 黄片wwwwww| 亚洲在线自拍视频| 亚洲av一区综合| 在线天堂最新版资源| 国产日本99.免费观看| 国产精品久久久久久精品电影| 日韩一区二区视频免费看| 欧美zozozo另类| 国产女主播在线喷水免费视频网站 | 亚洲在久久综合| 久久久久九九精品影院| 国产三级中文精品| 人妻久久中文字幕网| 波多野结衣高清作品| 国产探花极品一区二区| 亚洲四区av| 久久热精品热| 欧美另类亚洲清纯唯美| 变态另类丝袜制服| 亚洲欧美日韩东京热| 日本欧美国产在线视频| 国产三级在线视频| 日韩欧美国产在线观看| 嫩草影院精品99|