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

    小黑楊HD-Zip轉(zhuǎn)錄因子家族生物信息學及應答鹽脅迫分析

    2017-11-10 01:47:18張雪梅姚文靜姜廷波周博如
    植物研究 2017年5期
    關鍵詞:信息學楊樹結構域

    張雪梅 姚文靜 趙 凱 姜廷波 周博如

    (東北林業(yè)大學林木遺傳育種國家重點實驗室,哈爾濱 150040)

    小黑楊HD-Zip轉(zhuǎn)錄因子家族生物信息學及應答鹽脅迫分析

    張雪梅 姚文靜 趙 凱 姜廷波 周博如*

    (東北林業(yè)大學林木遺傳育種國家重點實驗室,哈爾濱 150040)

    HD-Zip轉(zhuǎn)錄因子基因是植物中特有的一類蛋白家族,在植物生長發(fā)育和逆境應答脅迫過程中發(fā)揮重要作用。HD-Zip轉(zhuǎn)錄因子基因是由高度保守的同源異型結構域(HD)和亮氨酸拉鏈域(LZ)結構域構成的特殊結構模型。楊樹HD-Zip轉(zhuǎn)錄因子家族共有63個基因,可被分為HD-ZipⅠ、HD-ZipⅡ、HD-ZipⅢ和HD-ZipⅣ四個亞家族。本文利用RNA-Seq分析了鹽脅迫條件下HD-Zip基因家族在小黑楊根、莖、葉等不同組織的基因表達差異,從轉(zhuǎn)錄組水平揭示其應答脅迫環(huán)境的分子機制,結果表明,鹽脅迫下在葉中有25個HD-Zip基因下調(diào)表達,21個基因上調(diào)表達;莖中有42個基因下調(diào)表達,11個基因上調(diào)表達;根中有26個基因下調(diào)表達,24個基因上調(diào)表達。另外,本文根據(jù)擬南芥HD-Zip轉(zhuǎn)錄因子家族基因的已知功能,預測了楊樹HD-Zip轉(zhuǎn)錄因子同源基因的功能,并利用生物信息學方法分析了楊樹HD-Zip轉(zhuǎn)錄因子蛋白序列的保守結構域、氨基酸組成和理化性質(zhì)等,為進一步研究楊樹HD-Zip轉(zhuǎn)錄因子基因功能提供參考。

    HD-Zip轉(zhuǎn)錄因子;非生物脅迫;生長發(fā)育;小黑楊;生物信息學

    楊樹生長速度快、適應性強、分布范圍廣、材質(zhì)優(yōu)良,是重要的綠化、造林和用材樹種,主要應用于四旁綠化、農(nóng)林防護、民用建筑、板料和紙漿制造方面。隨著楊樹全基因組測序的完成,楊樹成為研究植物分子生物學和遺傳學上的木本模式植物,在木本植物生長發(fā)育和植物抗旱、耐鹽、抗低溫等非生物脅迫的研究和應用中非常廣泛。

    轉(zhuǎn)錄因子(Transcription factor)也稱反式作用因子,是能與真核生物啟動子特定DNA序列結合的蛋白質(zhì)分子,或者是具有已知DNA結合域結構特征的蛋白質(zhì)分子[1]。轉(zhuǎn)錄因子能與基因5′端上游特定序列進行專一性結合,保證目的基因在特定的時間與空間表達。植物中的轉(zhuǎn)錄調(diào)控能嚴格控制生物表達過程,從而應答環(huán)境脅迫,調(diào)節(jié)對環(huán)境變化的反應。

    同源異型域—亮氨酸拉鏈蛋白(Homeodomain-leucine Zipper,HD-Zip)是植物特有的一類轉(zhuǎn)錄因子,屬于Homeobox蛋白家族,具有高度保守的HD(Homeodomain,HD)結構域,由61個氨基酸組成[2],可折疊成典型的螺旋—環(huán)—螺旋—轉(zhuǎn)角—螺旋的結構,能與DNA特異結合[3],其末端與亮氨酸拉鏈域(LZ)緊密相連,亮氨酸拉鏈結構域主要是介導蛋白二聚體的形成[4]。根據(jù)氨基酸序列同源性的不同,植物HD-Zip轉(zhuǎn)錄因子基因可分為HD-ZipⅠ、HD-ZipⅡ、HD-ZipⅢ和HD-ZipⅣ 4類[5],并且四大亞類均有一定的結構特征和相應的功能特點,HD-ZipⅠ主要包括HD和LZ兩個保守結構域,在應答非生物脅迫方面發(fā)揮重要作用[6];HD-ZipⅡ亞家族基因除了含有高度保守的HD和LZ序列外,在LZ下游緊挨著LZ結構域還有一個以五個保守氨基酸半胱氨酸(Cys)、脯氨酸(Pro)、絲氨酸(Ser)、半胱氨酸(Cys)和谷氨酸(Glu)命名的N端共有序列,主要負責細胞氧化還原傳感[7],該亞家族的基因負責感應光信號并能對植物激素產(chǎn)生應答效應[8];HD-ZipⅢ亞家族成員大多含有START(steroidogenic acute regulatory protein-related lipid transfer)結構域和與START相鄰的SAD(START-adjacent domain),除此以外,幾乎所有成員在C端均有一個可以結合信號分子來改變HD-ZipⅢ蛋白構象的MEKHLA結構域,該亞家族主要在胚胎形態(tài)發(fā)生、分生組織調(diào)控、側(cè)生器官形成、葉片的極性運輸以及生長素的運輸?shù)确矫姘l(fā)揮重要作用[9~10];HD-ZipⅣ亞家族成員相較于HD-ZipⅢ亞家族缺少了MEKHLA結構域,該亞家族成員主要參與表皮細胞的分化、花青素的積累、根的發(fā)育和毛狀體的形成[11~12]等。當前已對擬南芥、水稻[13]、蘋果、紫花苜蓿等植物的HD-Zip轉(zhuǎn)錄因子基因家族有過相關研究,但對楊樹的HD-Zip轉(zhuǎn)錄因子基因家族介紹較少。

    楊樹共有63個HD-Zip轉(zhuǎn)錄因子基因[14],本文利用轉(zhuǎn)錄組測序(RNA-Seq)法比較了楊樹在高鹽脅迫下分別在根、莖、葉中的相對表達量,從而分析了HD-Zip轉(zhuǎn)錄因子基因在不同組織中的差異性表達;并利用生物信息學的方法對楊樹的HD-Zip轉(zhuǎn)錄因子基因家族成員進行分析,為深度研究楊樹轉(zhuǎn)錄因子的功能提供依據(jù)。

    1 材料與方法

    1.1 非生物脅迫處理楊樹幼苗

    從同一株楊樹上采集嫩枝,截成15 cm,25℃室溫14 h·d-1光照水培兩個月,新長出根、莖(頂尖嫩莖處)、葉的楊樹苗即為本實驗材料。然后分別用水(對照)和0.15 mol·L-1NaCl溶液處理24 h,收集根、莖、葉樣品(每個處理有3個生物學重復),液氮速凍、-70℃保存,然后送生物公司提取RNA并利用IlluminaHiSeq 2000基因分析系統(tǒng)進行轉(zhuǎn)錄組高通量測序,最后進行RNA測序數(shù)據(jù)處理[15]。

    1.2 楊樹HD-Zip轉(zhuǎn)錄因子基因的生物信息學分析

    從楊樹基因組數(shù)據(jù)庫(http://planttfdb.cbi.pku.edu.cn/)中獲得了63條轉(zhuǎn)錄因子氨基酸序列,從擬南芥基因組數(shù)據(jù)庫(http://www.arabidopsis.org/)中得到48條HD-Zip轉(zhuǎn)錄因子氨基酸序列。利用BioEdit對楊樹HD-Zip轉(zhuǎn)錄因子蛋白序列進行比對分析,利用ClustalX2.1和MEME(Version 4.10.0)對楊樹HD-Zip轉(zhuǎn)錄因子蛋白序列進行基因保守元件(conserved motifs)的多元比較分析[16]。

    將獲得的楊樹HD-Zip轉(zhuǎn)錄因子蛋白氨基酸序列放入ExPaSy在線Protparam軟件(http://web/expasy/org/protparam/)進行分析,獲得HD-Zip基因的蛋白分子量和等電點等理化屬性。用PBIL在線網(wǎng)站(https://npsa-prabi.ibcp.fr/)對63條楊樹HD-Zip轉(zhuǎn)錄因子蛋白序列的二級結構進行預測分析,并用Swiss-Model程序從四個亞家族中隨機抽取3條蛋白序列進行三級結構同源建模。

    2 結果與分析

    2.1 鹽脅迫對楊樹HD-Zip轉(zhuǎn)錄因子基因表達影響

    根據(jù)本實驗室之前的鹽脅迫處理實驗,用0.15 mol·L-1鹽溶液處理楊樹苗24 h后葉片出現(xiàn)干枯現(xiàn)象,故以0.15 mol·L-1NaCl溶液作為脅迫處理楊樹的最大濃度。將楊樹苗浸泡在水和0.15 mol·L-1NaCl溶液中24 h,分別收集根、莖、葉樣品,經(jīng)轉(zhuǎn)錄組測序分析發(fā)現(xiàn),在葉中除了17個基因在鹽脅迫條件下沒有表達量變化以外,有25個HD-Zip基因下調(diào)表達,21個基因上調(diào)表達(表1);莖中有10個基因沒有表達量變化,42個基因下調(diào)表達,11個基因上調(diào)表達;根中有13個基因沒有表達量變化,26個基因下調(diào)表達,24個基因上調(diào)表達(圖1A),其中,在葉中上調(diào)表達2倍以上的有5個,莖中上調(diào)表達2倍以上的有5個,最高為5.14倍,根中上調(diào)表達2倍以上的有7個,最高為8.25倍。除此以外,Potri.014G103000.1、Potri.001G083700.1和Potri.007G008200.1等5個基因在鹽脅迫條件下均在根、莖、葉中顯著上調(diào)表達,Potri.017G081700.1、Potri.001G372300.1和Potri.008G148200.1等14個基因在根、莖、葉中均顯著下調(diào)表達(圖1B),而在葉中上調(diào)4.55倍的Potri.007G029500.1基因在莖和根中表現(xiàn)下調(diào)趨勢,還有的基因如Potri.009G023600.1、Potri.012G139300.1和Potri.001G229700.1在葉和根中表達量均顯著增加,在莖中卻呈現(xiàn)下調(diào)趨勢。結果表明了,楊樹HD-Zip轉(zhuǎn)錄因子家族各成員在不同的組織中參與植物的抗逆過程,從而實現(xiàn)高鹽應答脅迫。

    2.2 楊樹HD-Zip轉(zhuǎn)錄因子家族生物信息學分析

    利用BioEdit對63條楊樹HD-Zip轉(zhuǎn)錄因子蛋白結構域和48條擬南芥HD-Zip轉(zhuǎn)錄因子蛋白結構域進行序列比對,根據(jù)擬南芥HD-Zip基因家族的分類[17], 將楊樹的63個HD-Zip基因同樣分成四個亞家族HD-ZipⅠ、HD-ZipⅡ、HD-ZipⅢ和HD-ZipⅣ(表2),從表中可以看出,楊樹有22個HD-Zip基因?qū)儆贖D-ZipⅠ亞家族,17個基因?qū)儆贖D-ZipⅡ亞家族,8個基因?qū)儆贖D-ZipⅢ亞家族,16個基因?qū)儆贖D-ZipⅣ亞家族。

    表1HD-Zip家族基因應答鹽脅迫的相對表達量

    Table1RelativeexpressionofHD-Zipfamilygenesundersaltstress

    基因IDGeneID葉中相對表達量Relativeexpressioninleaves莖中相對表達量Relativeexpressioninstems根中相對表達量RelativeexpressioninrootsPotri.006G071600.10.022±0.001-1.383±0.0048.246±0.686Potri.004G020400.10.048±0.007-0.438±0.0170.048±0.003Potri.008G129500.10.108±0.022-0.403±0.0110.912±0.037Potri.018G045100.10.228±0.0010-1.826±0.144Potri.003G052400.10.250±0.0201.045±0.017-0.282±0.025Potri.002G113400.10.336±0.040-1.502±0.069-1.185±0.064Potri.010G112600.10.403±0.0142.912±0.2060.699±0.101Potri.014G045100.10.445±0.016-0.716±0.0321.212±0.026Potri.015G065400.10.477±0.0960.350±0.013-2.969±0.098Potri.007G135100.10.481±0.021-2.785±0.036-1.293±0.132Potri.001G229700.10.606±0.129-0.103±0.0331.791±0.709Potri.007G008200.10.733±0.0480.365±0.0380.357±0.025Potri.012G139300.10.794±0.014-0.355±0.0244.411±0.315Potri.012G070900.11.078±0.0380.845±0.107-1.408±0.521Potri.009G023600.11.420±0.269-0.022±0.0012.389±0.510Potri.002G136400.11.488±0.0330.993±0.0211.488±0.016Potri.012G023700.12.199±0.1892.547±0.0113.805±0.102Potri.002G176300.12.870±0.1052.089±0.195-1.191±0.976Potri.007G029500.14.555±0.138-0.020±0.004-0.144±0.061Potri.001G083700.14.642±0.1313.847±0.5321.209±0.121Potri.014G103000.14.843±0.0595.143±0.3821.116±0.143

    圖1 鹽脅迫條件下楊樹HD-Zip家族基因表達 A.HD-Zip家族基因在不同組織中的表達;B.上調(diào)表達基因和下調(diào)表達基因Fig.1 Expression of HD-Zip family genes in poplar by salt stress A.Expression of HD-Zip family genes in different tissues; B.Up-regulated and down-regulated genes

    Phytozome登錄號PhytozomeID類型Type染色體Chromosome大小Size(aa)分子量Molecularweight等電點TheoreticalPI脂肪系數(shù)Aliphaticindex疏水性平均系數(shù)Grandaverageofhydropathicityα?螺旋Alphahelix(%)二級結構預測Secondarystructureprediction延伸鏈Extendedstrand(%)無規(guī)則鏈Randomcoil(%)Potri.007G097100.1Ⅰ732836845.64.8665.46-0.76742.389.4548.17Potri.010G093400.1Ⅰ1031135250.56.2269.07-0.78136.0111.5852.41Potri.005G071900.1Ⅰ532936971.74.8266.41-0.78536.4713.9849.54Potri.008G148200.1Ⅰ828932985.85.9464.19-0.86233.918.3057.79Potri.001G155100.1Ⅱ130233642.98.7464.3-0.78227.4811.5960.93Potri.017G081700.1Ⅰ1730634748.86.2764.12-0.84338.246.5455.23Potri.015G065400.1Ⅰ15314354934.6766.82-0.81840.4510.8348.73Potri.012G070900.1Ⅰ12317354944.762.81-0.79840.387.5752.05Potri.016G059000.1Ⅱ1634438311.98.9861.89-0.79722.3814.8362.79Potri.009G023600.1Ⅱ938342234.56.0356.37-0.8423.5024.8051.70Potri.001G229700.1Ⅱ137441439.77.6551.23-0.90122.1920.3257.49Potri.006G193700.1Ⅱ632636355.99.2863.19-0.7924.8515.6459.51Potri.007G008200.1Ⅱ727330530.58.6865.19-0.70934.8015.3849.82Potri.014G045100.1Ⅱ14298331458.463.89-0.82830.5410.4059.06Potri.002G136400.1Ⅱ232035667.18.8865.59-0.83530.639.3860.00Potri.005G147100.1Ⅱ529032480.77.5370.31-0.6828.9713.4557.59Potri.002G100600.1Ⅰ226530733.16.2368.45-0.95549.817.1743.02Potri.008G194400.1Ⅰ830235004.25.0463.61-0.828.4814.9056.62Potri.012G023700.1Ⅰ12275313985.6163.16-0.79153.8210.9135.27Potri.018G133100.1Ⅱ18207239258.7959.37-1.10532.8511.1156.04Potri.014G103000.1Ⅰ1423827160.24.9759.87-0.91637.399.6652.94Potri.016G058600.1Ⅱ1620323012.19.7762.07-0.88735.4714.7849.75Potri.002G176300.1Ⅰ223827282.45.0459.87-0.93836.978.4054.62Potri.006G071600.1Ⅱ621424587.58.3254.72-1.18632.2411.2156.54Potri.002G113400.1Ⅱ2269300706.1675.39-0.64228.2517.1054.65Potri.005G161500.1Ⅰ526530781.15.6769.58-0.91148.3010.1941.51Potri.005G126100.1Ⅰ520423701.99.0270.74-0.87533.8214.2251.96Potri.003G079800.1Ⅱ335439484.56.6868.08-0.66723.7313.5662.71Potri.006G117700.1Ⅰ623227132.48.4563.92-0.9547.846.9045.26Potri.007G029500.1Ⅰ720423834.16.0469.8-0.89850.988.3340.69Potri.008G129500.1Ⅱ8232255738.4867.28-0.62128.4517.6753.88Potri.010G112600.1Ⅱ1022525185.49.1665.91-0.7332.8916.0051.11Potri.001G083700.1Ⅰ119623014.48.9458.21-1.24552.553.0644.39Potri.011G114500.1Ⅰ1121324526.36.4569.15-0.86846.484.6948.83Potri.001G395100.1Ⅰ120523745.78.7776.59-0.64947.809.7642.44Potri.007G135100.1Ⅰ717920311.78.5576.87-0.74549.169.5041.34Potri.017G016900.1Ⅰ1717920296.97.8376.87-0.70453.073.9143.02Potri.016G058400.1Ⅱ1614416280.69.9262.36-0.75131.9420.8347.22Potri.003G147200.1Ⅰ39511456.29.3492.42-0.39466.3214.7418.95Potri.003G052400.1Ⅳ375984384.36.1174.91-0.50837.0217.9245.06Potri.001G184100.1Ⅳ176485053.26.0374.54-0.47938.0916.1045.81Potri.002G230200.1Ⅳ276583840.85.677.88-0.31830.4617.6551.90Potri.014G152000.1Ⅳ1476283224.45.5679.23-0.29237.2714.9647.77Potri.004G020400.1Ⅳ472579870.55.5479.89-0.31834.7617.1048.14Potri.011G025000.1Ⅳ1172078975.55.5680.29-0.29732.6418.4748.89Potri.014G075200.1Ⅳ1482389393.45.8981.41-0.27434.5116.7748.72Potri.002G154700.1Ⅳ282089297.25.9381.35-0.26935.8518.0546.10Potri.015G034100.1Ⅳ1572579739.56.3381.46-0.25730.0719.1750.76Potri.012G038500.1Ⅳ1271178026.46.0883.5-0.2935.5815.6148.80Potri.006G237500.1Ⅲ684492883.35.9783.58-0.13142.5416.1141.35Potri.015G141800.1Ⅳ1575383020.35.7183.84-0.2935.8618.0646.08Potri.018G045100.1Ⅲ1884392805.26.0785.41-0.14240.3317.9141.76Potri.011G098300.1Ⅲ1184392739.85.9585.99-0.1442.7017.0840.21Potri.012G139300.1Ⅳ1274882431.45.5684.81-0.32433.4218.8547.73Potri.001G188800.1Ⅲ185293855.36.0687.14-0.13740.4917.0242.49Potri.001G372300.1Ⅲ184492675.75.9684.62-0.15242.7716.1141.11Potri.003G050100.1Ⅲ383792230.46.0386.71-0.16142.1715.4142.41Potri.003G096000.1Ⅳ380889791.85.4380.82-0.31128.7119.9351.36Potri.001G137800.1Ⅳ180489190.15.3580.63-0.31528.1119.1552.74Potri.009G014500.1Ⅲ985593554.55.990.51-0.07343.9816.3739.65Potri.004G211300.1Ⅲ4846928415.8691.82-0.04648.2314.8936.88Potri.017G144600.1Ⅳ1764875924.36.3882.54-0.28130.2619.0150.73Potri.T120700.1ⅣScaffold?20164170604.86.4681.86-0.39935.1016.85 48.05

    分析楊樹HD-Zip轉(zhuǎn)錄因子基因家族的理化性質(zhì)可知,楊樹HD-Zip轉(zhuǎn)錄因子基因家族的氨基酸個數(shù)在95~855,分布在HD-ZipⅠ和HD-ZipⅡ亞家族的基因氨基酸個數(shù)較少,一般在100~400個,而HD-ZipⅠ亞家族中的Potri.003G147200.1僅有95個氨基酸。分布在HD-ZipⅢ和HD-ZipⅣ亞家族中的基因氨基酸個數(shù)在641~855,氨基酸個數(shù)最多的是HD-ZipⅢ中的Potri.009G014500.1,多達855個。楊樹HD-Zip轉(zhuǎn)錄因子基因家族蛋白的等電點大多在酸性范圍內(nèi),脂肪系數(shù)均在51.23以上,最高達到92.42,這表明了HD-Zip轉(zhuǎn)錄因子家族的蛋白熱穩(wěn)定性較高。其疏水性平均系數(shù)均為負數(shù),說明HD-Zip蛋白均表現(xiàn)為親水性。63個楊樹HD-Zip蛋白序列中主要有α-螺旋、延伸鏈和無規(guī)則卷曲結構且三種結構交錯排列。利用planttfdb數(shù)據(jù)庫找到63個楊樹HD-Zip轉(zhuǎn)錄因子氨基酸序列,根據(jù)popgenie在線網(wǎng)站上Chromosome Diagram功能,進行染色體定位發(fā)現(xiàn),HD-Zip基因在楊樹染色體上呈不均勻分布,在第1、2、3號染色體上分布較多,分別有8、6、5條基因,在其它染色體上分布較少,第4、9、10、16和18號染色體上均只存在2條HD-Zip基因,并且在第13和第19號染色體上不存在HD-Zip基因。HD-Zip基因在染色體上的這種不均衡分布可能跟楊樹進化過程中基因的功能特性相關。

    利用MEME在線軟件確定了20個保守序列(表3),根據(jù)20個氨基酸保守序列可將63個楊樹HD-Zip基因分成四類,與系統(tǒng)發(fā)育樹的分類一致。對63條HD-Zip基因分析得到,基序1和基序2編碼高度保守的HD結構域,基序17編碼亮氨酸拉鏈結構域(LZ),而基序3、7、8、9、11、12、13、14、16、18可能與HD-ZipⅢ和HD-ZipⅣ中的START結構域的形成有關,基序4、5、6、10、13、15、19、20則可能與僅存在與HD-ZipⅢ中的MEKHLA結構域的編碼有關。不同的亞家族包含有不同的保守結構域,這可能與HD-Zip基因蛋白特殊的生物學功能相關。

    利用Swiss-Model分別HD-Zip 4個亞家族中的各個成員進行三級結構的同源建模,結果表明,HD-ZipⅠ、HD-ZipⅡ和HD-ZipⅣ亞家族中蛋白的三級結構大體相同,而HD-ZipⅢ亞家族中的蛋白三級結構的構象之間略有不同,這可能與HD-ZipⅢ中基因的特殊保守域MEKHLA相關,從而使HD-ZipⅢ中不同的成員具有不同的功能[18],而且四個亞家族之間的蛋白構型也存在一定的差異。分析其二級結構發(fā)現(xiàn),HD-ZipⅠ亞家族中蛋白的α-螺旋所占比例偏大,HD-ZipⅡ亞家族中蛋白的無規(guī)則卷曲結構所占比例偏大,而HD-ZipⅢ和HD-ZipⅣ亞家族中蛋白的延伸鏈的含量較HD-ZipⅠ和HD-ZipⅡ偏大,因此不同亞家族之間蛋白三級結構的差異可能與蛋白中的α-螺旋、延伸鏈和無規(guī)則卷曲結構的含量有關。

    表3HD-Zip基序的長度和保守氨基酸序列

    Table3ThelengthandconservedaminoacidsequenceofHD-Zipmotifs

    基序Motif長度Length保守氨基酸序列ConservedaminoacidsequenceMotif150LGLRPRQVEVWFQNRRARTKLKQTEVDCEFLKRCCENLTAENRRLQKEVQMotif229RYHRHTARQIQEMEAMFKECPHPDDKQRMMotif329GVAGNYNGALQVMTAEFQVPSPLVPTRENYFVRYCKQMotif450IEPKQIKVWFQNRRCREKQRKEASRLQTVNRKLTAMNKLLMEENDRLQKQMotif541KATGTAVEWVQMPGMKPGPDSIGIVAISHGCTGVAARACGLMotif650PVFTFANQAGLDMLETTLVALQDITLEKIFDDHGRKTLCSEFPQIMQQGFMotif724THQLYRPLISSGMGFGAQRWIATLMotif841TSFWLPVPPKRVFDFLRDENTRNEWDILSNGGVVQEMAHIAMotif941QESCTDSTGSYVIYAPVDISAMNIVLSGGDPDYVALLPSGFMotif1050DVFLLQLCCGMDENAVGTCAELIFAPIDATFADDAPLLPSGFRIIPLDSGMotif1115LPSGCLIQDMPNGYSMotif1229RMVISFCAGVSASTAHTWTTLSGTGADDVMotif1332SLLTVAFQILVNSLPTAKLTVESVETVNNLISMotif1441RSEASRESAVVIMNHVNLVEILMDANQWSTIFCGIVSRAMTMotif1550VLCAKASMLLQNVPPAILLRFLREHRSEWADSGIDAYAAAAVKAGPCSLPMotif1630DKPMIIELAVAAMEELVRMAQMDEPLWMGSMotif1721PATTLTMCPSCERVATTTTSSMotif1841FDEQKLRIENAQLKEELDRVSSIAAKYIGRPISQLPPVQPVMotif1929CLPGGICLSSMGRPVSYERAVAWKVLNEEMotif2029VGLEPTRVAEILKDRPSWFRDCRAVDVLN

    3 討論

    HD-Zip轉(zhuǎn)錄因子在植物的生長發(fā)育過程中發(fā)揮重要功能,不同的亞家族轉(zhuǎn)錄因子發(fā)揮不同的功能,如HD-ZipⅠ蛋白在番茄的生長發(fā)育過程中參與果實成熟的調(diào)控[19],HD-ZipⅡ亞家族中基因的過表達誘導擬南芥的早花表型[20],HD-ZipⅢ中的成員調(diào)控木質(zhì)部維管束形成[21]、頂端分生組織的發(fā)生[22],HD-ZipⅣ亞家族中的基因成員在玉米中過表達后抑制表皮毛的形成[23]等。本文通過轉(zhuǎn)錄組測序比較了鹽脅迫條件下楊樹HD-Zip轉(zhuǎn)錄因子基因在根、莖、葉中的差異性表達,找到了鹽脅迫下在葉中高表達的21個HD-Zip基因,莖中11個高表達的HD-Zip基因以及根中24個高表達的HD-Zip基因,分析發(fā)現(xiàn)Potri.014G103000.1等5個基因經(jīng)過高鹽處理后均顯著上調(diào)表達,Potri.017G081700.1等14個基因在根、莖、葉中均下調(diào)表達,而有的基因則選擇性地在不同的組織中高表達,表明了楊樹HD-Zip轉(zhuǎn)錄因子家族各成員在不同的組織中參與植物的抗逆過程,從而發(fā)揮鹽脅迫應答作用。

    本文利用生物信息學的方法分析了63個楊樹HD-Zip轉(zhuǎn)錄因子基因的結構及其相關信息,通過與擬南芥HD-Zip轉(zhuǎn)錄因子家族中的48個基因進行比較以及對其包含的保守結構域進行分析將其分為四個亞家族,并分析了整個家族中63個基因的理化性質(zhì)和染色體定位,得到分布在17條染色體上的63個HD-Zip轉(zhuǎn)錄因子蛋白大多處于酸性范圍內(nèi),而等電點值較大的基因蛋白大多含有較多的亮氨酸等堿性氨基酸,有助于形成亮氨酸拉鏈保守結構域,從而介導蛋白二聚體的形成;通過蛋白二級結構和三級結構的預測分析發(fā)現(xiàn)HD-ZipⅠ、HD-ZipⅡ和HD-ZipⅣ亞家族中各成員蛋白結構相似,而HD-ZipⅢ中各成員蛋白結構略有不同,可能與其特殊的結構域有關。統(tǒng)計得到,鹽脅迫條件下表現(xiàn)出上調(diào)表達趨勢的21個HD-Zip轉(zhuǎn)錄因子基因除Potri.018G045100.1屬于HD-ZipⅢ亞家族外,有3個基因?qū)儆贖D-ZipⅣ亞家族,而屬于HD-ZipⅠ和HD-ZipⅡ亞家族的分別有8個和9個;在根、莖、葉中應答鹽脅迫均呈上調(diào)表達趨勢的基因除Potri.007G008200.1屬于HD-ZipⅡ亞家族外,其余都屬于HD-ZipⅠ亞家族。這些基因在染色體上呈不均勻分布,其氨基酸個數(shù)與其分子量不大,除了Potri.001G083700.1基因的等電點處于堿性范圍外,其它基因均處于酸性范圍內(nèi),其脂肪系數(shù)均在58.21以上,蛋白質(zhì)的熱穩(wěn)定性較高,疏水性平均值在-0.791~-1.245,是親水性極強的蛋白,這些基因所表現(xiàn)出的較強的脅迫應答能力可能與其特殊的理化性質(zhì)以及HD-ZipⅠ亞家族中轉(zhuǎn)錄因子的HD和LZ結構域有關,而HD-ZipⅢ和HD-ZipⅣ中的START和MEKHLA結構域在脅迫處理環(huán)境中可能抑制基因的表達??傊?,通過RNA-Seq和生物信息學方法的分析,將進一步為研究楊樹HD-Zip轉(zhuǎn)錄因子家族蛋白在應答逆境脅迫方面提供有效依據(jù)。

    1.Yanagisawas.Dof domain proteins:plant-specific transcription factors associated with diverse phenomena unique to plants[J].Plant and Cell Physiology,2004,45(4):386-391.

    2.Schena M,Davis R W.HD-Zip proteins:members of anArabidopsishomeodomain protein superfamily[J].Proceedings of the National Academy of Sciences of the United States of America,1992,89(9):3894-3898.

    3.Bharathan G,Janssen B J,Kellogg E A,et al.Did homeodomain proteins duplicate before the origin of angiosperms,fungi,andmetazoa?[J].Proceedings of the National Academy of Sciences of the United States of America,1997,94(25):13749-13753.

    4.Jain M,Tyagi A K,Khurana J P.Genome-wide identification,classification,evolutionary expansion and expression analyses of homeobox genes in rice[J].The FEBS Journal,2008,275(11):2845-2861.

    5.Sessa G,Steindler C,Morelli G,et al.TheArabidopsisAthb-8,-9 and -14 genes are members of a small gene family coding for highly related HD-ZIP proteins[J].Plant Molecular Biology,1998,38(4):609-622.

    6.Olsson A S,Engstr?m P,S?derman E.The homeobox genesATHB12 andATHB7 encode potential regulators of growth in response to water deficit inArabidopsis[J].Plant Molecular Biology,2004,55(5):663-677.

    7.Tron A E,Bertoncini C W,Chan R L,et al.Redox regulation of plant homeodomain transcription factors[J].Journal of Biological Chemistry,2002,277(38):34800-34807.

    8.Sawa S,Ohgishi M,Goda H,et al.TheHAT2 gene,a member of the HD-Zip gene family,isolated as an auxin inducible gene by DNA microarray screening,affects auxin response inArabidopsis[J].The Plant Journal,2002,32(6):1011-1022.

    9.Prigge M J,Otsuga D,Alonso J M,et al.ClassⅢ homeodomain-leucine zipper gene family members have overlapping,antagonistic,and distinct roles inArabidopsisdevelopment[J].The Plant Cell,2005,17(1):61-76.

    10.Otsuga D,Deguzman B,Prigge M J,et al.REVOLUTAregulates meristem initiation at lateral positions[J].The Plant Journal,2001,25(2):223-236.

    11.Nakamura M,Katsumata H,Abe M,et al.Characterization of the class Ⅳ homeodomain-leucine zipper gene family inArabidopsis[J].Plant Physiology,2006,141(4):1363-1375.

    12.Luo D L,Oppenheimer D G.Genetic control of trichome branch number inArabidopsis:the roles of theFURCAloci[J].Development,1999,126(24):5547-5557.

    13.Agalou A,Purwantomo S,?vern?se,et al.A genome-wide survey of HD-Zip genes in rice and analysis of drought-responsive family members[J].Plant Molecular Biology,2008,66(1-2):87-103.

    14.Hu R B,Chi X Y,Chai GH,et al.Genome-wide identification,evolutionary expansion,and expression profile of homeodomain-leucine zipper gene family in poplar(Populustrichocarpa)[J].PLoS One,2012,7(2):e31149.

    15.Yao W J,Wang S J,Zhou B R,et al.Transgenic poplar overexpressing the endogenous transcription factor ERF76 gene improves salinity tolerance[J].Tree Physiology,2016,36(7):896-908.

    16.Bailey T L,Elkan C.Fitting a mixture model by expectation maximization to discover motifs in biopolymers[J].Proceeding International Conference on Intelligent Systems for Molecular Biology,1994,2:28-36.

    17.Aso K,Kato M,Banks JA,et al.Characterization of homeodomain-leucine zipper genes in the fernCeratopterisrichardiiand the evolution of the homeodomain-leucine zipper gene family in vascular plants[J].Molecular Biology and Evolution,1999,16(4):544-552.

    19.Lin Z F,Hong Y G,Yin M G,et al.A tomato HD-Zip homeobox protein,LeHB-1,plays an important role in floral organogenesis and ripening[J].The Plant Journal,2008,55(2):301-310.

    20.Dezar C A,Giacomelli J I,Manavella P A,et al.HAHB10,a sunflower HD-ZipⅡ transcription factor,participates in the induction of flowering and in the control of phytohormone-mediated responses to biotic stress[J].Journal of Experimental Botany,2011,62(3):1061-1076.

    21.Itoh J I,Hibara K I,Sato Y,et al.Developmental role and auxin responsiveness of classⅢ homeodomain leucine zipper gene family members in rice[J].Plant Physiology,2008,147(4):1960-1975.

    22.Kim J,Jung J H,Reyes J L,et al.microRNA-directed cleavage ofATHB15mRNA regulates vascular development inArabidopsisinflorescence stems[J].The Plant Journal,2005,42(1):84-94.

    23.Vernoud V,Laigle G,Rozier F,et al.The HD-ZIPⅣ transcription factor OCL4 is necessary for trichome patterning and anther development in maize[J].The Plant Journal,2009,59(6):883-894.

    National 863 project(2013AA102701)

    introduction:ZHANG Xue-Mei(1991—),female,master,mainly engaged in the forestry research of stress resistance.

    date:2017-07-14

    BioinformaticsandResponsetoSaltStressAnalysisoftheHD-ZipTranscriptionFactorFamilyinPopulussimonii×P.nigra

    ZHANG Xue-Mei YAO Wen-Jing ZHAO Kai JIANG Ting-Bo ZHOU Bo-Ru*

    (State Key Laboratory of Tree Genetics and Breeding,Northeast Forestry University,Harbin 150040)

    The HD-Zip transcription factor family is a kind of special protein in plant, playing crucial roles in plant development processes and abiotic stresses response. The protein family includes a highly conserved homeodomain(HD) and a leucine zipper domain(LZ), forming a special structural model. The HD-Zip transcription factor family, including 63 members, can be classified into 4 subfamilies HD-Zip Ⅰ, HD-Zip Ⅱ, HD-Zip Ⅲ and HD-Zip Ⅳ. RNA-Seq analysis of poplar roots, stems and leaves under salt stress was conducted for the differentially expressed genes in different tissues, and reveal the molecular mechanism of the response to stress environment. Therewere 25 HD-Zip genes down-regulated and 21 HD-Zip genes up-regulated in leaves under salt treatment, 42 genes down-regulated and 11 genes up-regulated in stems, and 26 genes down-regulated and 24 genes up-regulated in roots. According to gene annotations in Arabidopsis thaliana, the functions of HD-Zip transcription factor in poplarwere predicted.The conserved domains, structural characteristics and chemical properties of HD-Zip family members were also studied by bioinformatics, providing reference for gene functions of HD-Zip family in poplar.

    HD-Zip transcription factor;abiotic stresses;growth and development;Populussimonii×P.nigra;bioinformatics

    國家863課題資助(2013AA102701)

    張雪梅(1991—),女,碩士研究生,主要從事林木抗逆機理方面研究。

    * 通信作者:E-mail:boruzhou@yahoo.com

    2017-07-14

    * Corresponding author:E-mail:boruzhou@yahoo.com

    S722.3+6

    A

    10.7525/j.issn.1673-5102.2017.05.011

    猜你喜歡
    信息學楊樹結構域
    雞NRF1基因啟動子區(qū)生物信息學分析
    蛋白質(zhì)結構域劃分方法及在線服務綜述
    跟蹤導練(五)
    初論博物館信息學的形成
    中國博物館(2018年2期)2018-12-05 05:28:50
    口算本失蹤記
    楊樹山漫畫作品欣賞
    喜劇世界(2017年10期)2017-06-01 12:39:29
    重組綠豆BBI(6-33)結構域的抗腫瘤作用分析
    組蛋白甲基化酶Set2片段調(diào)控SET結構域催化活性的探討
    miRNA-148a在膀胱癌組織中的表達及生物信息學分析
    泛素結合結構域與泛素化信號的識別
    高清av免费在线| 亚洲全国av大片| 欧美精品啪啪一区二区三区| 久久影院123| 成人黄色视频免费在线看| 国产精品久久久人人做人人爽| 欧美日韩av久久| 成年版毛片免费区| 999久久久精品免费观看国产| 精品一区二区三区视频在线观看免费 | 亚洲色图av天堂| 精品一区二区三区视频在线观看免费 | 97碰自拍视频| 亚洲精品在线观看二区| 色在线成人网| 国产精品自产拍在线观看55亚洲| 精品国产一区二区久久| 精品卡一卡二卡四卡免费| 亚洲熟女毛片儿| 国产av一区二区精品久久| 精品国产乱码久久久久久男人| 一进一出抽搐动态| 中文字幕av电影在线播放| 国产精品亚洲一级av第二区| 欧美黑人精品巨大| 18禁裸乳无遮挡免费网站照片 | 久久久久久久久久久久大奶| 狂野欧美激情性xxxx| 91在线观看av| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人精品无人区| 国产区一区二久久| 成人永久免费在线观看视频| 日韩成人在线观看一区二区三区| 免费高清在线观看日韩| 久久草成人影院| 韩国av一区二区三区四区| 伊人久久大香线蕉亚洲五| 亚洲熟妇熟女久久| 波多野结衣一区麻豆| 丰满迷人的少妇在线观看| 又紧又爽又黄一区二区| 国产精品久久久人人做人人爽| 久久亚洲真实| 亚洲avbb在线观看| 又紧又爽又黄一区二区| 亚洲av成人av| 国产国语露脸激情在线看| 脱女人内裤的视频| 老熟妇仑乱视频hdxx| 老汉色av国产亚洲站长工具| 久久午夜亚洲精品久久| 久久 成人 亚洲| 纯流量卡能插随身wifi吗| 99在线人妻在线中文字幕| 波多野结衣一区麻豆| 亚洲精品久久午夜乱码| 丝袜在线中文字幕| 久久久久久大精品| 桃色一区二区三区在线观看| 国产伦一二天堂av在线观看| 亚洲av成人不卡在线观看播放网| 午夜亚洲福利在线播放| svipshipincom国产片| 国产精品99久久99久久久不卡| 午夜福利欧美成人| www.www免费av| 黄色丝袜av网址大全| 亚洲欧美日韩另类电影网站| 老汉色av国产亚洲站长工具| 欧美日韩福利视频一区二区| 好看av亚洲va欧美ⅴa在| 日本五十路高清| 亚洲性夜色夜夜综合| av片东京热男人的天堂| 欧美激情久久久久久爽电影 | 大型av网站在线播放| 亚洲第一欧美日韩一区二区三区| a级毛片黄视频| 男女下面进入的视频免费午夜 | 在线观看免费视频网站a站| 久久久久久久精品吃奶| 丰满迷人的少妇在线观看| 久久久久久免费高清国产稀缺| www.999成人在线观看| 99在线视频只有这里精品首页| 亚洲欧美日韩高清在线视频| 国产av在哪里看| 俄罗斯特黄特色一大片| 正在播放国产对白刺激| 国产99白浆流出| 长腿黑丝高跟| e午夜精品久久久久久久| 国产亚洲精品一区二区www| 91成人精品电影| 51午夜福利影视在线观看| 成人亚洲精品av一区二区 | 在线观看日韩欧美| 久久精品亚洲精品国产色婷小说| 午夜福利一区二区在线看| 这个男人来自地球电影免费观看| 伦理电影免费视频| 身体一侧抽搐| 国产一区二区在线av高清观看| 国产精品美女特级片免费视频播放器 | netflix在线观看网站| 国产三级在线视频| 一级片免费观看大全| 国产人伦9x9x在线观看| 18禁观看日本| 琪琪午夜伦伦电影理论片6080| 日本wwww免费看| 好看av亚洲va欧美ⅴa在| 巨乳人妻的诱惑在线观看| 国产精品久久久久久人妻精品电影| 日本三级黄在线观看| 亚洲av五月六月丁香网| 亚洲精品粉嫩美女一区| 热99re8久久精品国产| 香蕉久久夜色| 国产欧美日韩一区二区三| 大香蕉久久成人网| 精品国内亚洲2022精品成人| 亚洲国产毛片av蜜桃av| 久久精品亚洲av国产电影网| 窝窝影院91人妻| 999久久久国产精品视频| 免费在线观看完整版高清| 一边摸一边做爽爽视频免费| 伊人久久大香线蕉亚洲五| 18禁国产床啪视频网站| 在线观看一区二区三区激情| videosex国产| 精品无人区乱码1区二区| 在线播放国产精品三级| 香蕉久久夜色| 大型黄色视频在线免费观看| 巨乳人妻的诱惑在线观看| 深夜精品福利| 一区二区日韩欧美中文字幕| 91成年电影在线观看| 成人手机av| 国内久久婷婷六月综合欲色啪| 久久精品亚洲精品国产色婷小说| 欧美一区二区精品小视频在线| 精品国产美女av久久久久小说| 久久精品国产亚洲av高清一级| 久久这里只有精品19| 黄色a级毛片大全视频| 国产精品一区二区精品视频观看| 国产av在哪里看| 国产欧美日韩精品亚洲av| 天堂中文最新版在线下载| 欧美日韩亚洲综合一区二区三区_| 国产精品久久久久久人妻精品电影| 国产极品粉嫩免费观看在线| 午夜成年电影在线免费观看| 亚洲精品国产一区二区精华液| 日韩大尺度精品在线看网址 | 亚洲国产欧美网| 亚洲精品在线美女| 精品人妻在线不人妻| 久久精品国产亚洲av高清一级| 女人被狂操c到高潮| 看黄色毛片网站| 老熟妇仑乱视频hdxx| 国产欧美日韩精品亚洲av| 18禁裸乳无遮挡免费网站照片 | 国产成人一区二区三区免费视频网站| 黄色成人免费大全| 俄罗斯特黄特色一大片| 我的亚洲天堂| 99re在线观看精品视频| 国产成人一区二区三区免费视频网站| 亚洲欧美精品综合一区二区三区| 极品教师在线免费播放| 正在播放国产对白刺激| 日韩成人在线观看一区二区三区| 国产精华一区二区三区| 日韩国内少妇激情av| 在线免费观看的www视频| 黄色片一级片一级黄色片| 99国产精品一区二区三区| 欧美激情 高清一区二区三区| 亚洲人成伊人成综合网2020| 国产99白浆流出| 久久中文字幕人妻熟女| 亚洲在线自拍视频| 在线观看免费视频网站a站| 久久久国产成人免费| 不卡av一区二区三区| 欧美激情 高清一区二区三区| 久久久精品国产亚洲av高清涩受| 一区二区三区激情视频| 51午夜福利影视在线观看| 人人澡人人妻人| 久久久国产欧美日韩av| 亚洲人成电影免费在线| 欧美成人午夜精品| 18禁美女被吸乳视频| 91老司机精品| 亚洲国产精品sss在线观看 | 亚洲欧美激情在线| 高清黄色对白视频在线免费看| 高潮久久久久久久久久久不卡| 亚洲自拍偷在线| 久久精品亚洲av国产电影网| 色婷婷久久久亚洲欧美| 亚洲精品美女久久av网站| 精品欧美一区二区三区在线| 男人舔女人的私密视频| 两性夫妻黄色片| 亚洲自偷自拍图片 自拍| 免费在线观看完整版高清| 最新在线观看一区二区三区| 成人av一区二区三区在线看| 免费看十八禁软件| 亚洲一区二区三区不卡视频| 成年人免费黄色播放视频| 亚洲avbb在线观看| 久久人人精品亚洲av| 国产精品偷伦视频观看了| 美国免费a级毛片| 午夜免费观看网址| 亚洲男人的天堂狠狠| 久久人人97超碰香蕉20202| 夜夜爽天天搞| 午夜福利影视在线免费观看| 亚洲成av片中文字幕在线观看| 少妇的丰满在线观看| 在线观看www视频免费| 18禁美女被吸乳视频| 久久青草综合色| 在线观看免费视频网站a站| 热99re8久久精品国产| 精品少妇一区二区三区视频日本电影| 俄罗斯特黄特色一大片| 国产精品 国内视频| 黄色怎么调成土黄色| 国产精品久久久av美女十八| 欧美日韩一级在线毛片| av免费在线观看网站| 91成年电影在线观看| 国产aⅴ精品一区二区三区波| 日韩成人在线观看一区二区三区| 日本wwww免费看| 久久精品91蜜桃| 高清黄色对白视频在线免费看| 黄色毛片三级朝国网站| 国产成人系列免费观看| 在线观看免费午夜福利视频| 侵犯人妻中文字幕一二三四区| 久久中文字幕人妻熟女| 一区二区三区国产精品乱码| 国产精品一区二区三区四区久久 | 交换朋友夫妻互换小说| 亚洲自偷自拍图片 自拍| 国产一区二区三区综合在线观看| 搡老乐熟女国产| 国产又爽黄色视频| 国产精品亚洲av一区麻豆| 身体一侧抽搐| 91国产中文字幕| 欧美日韩瑟瑟在线播放| 免费观看人在逋| 日韩大尺度精品在线看网址 | 久久久久国内视频| 国产乱人伦免费视频| 久久中文字幕人妻熟女| 精品卡一卡二卡四卡免费| 88av欧美| 精品第一国产精品| 在线永久观看黄色视频| 亚洲avbb在线观看| 欧美在线黄色| 中文字幕色久视频| 成人国产一区最新在线观看| 真人一进一出gif抽搐免费| 国产精品国产高清国产av| 亚洲av日韩精品久久久久久密| 欧美大码av| 国产一区在线观看成人免费| 老司机午夜福利在线观看视频| 女人被躁到高潮嗷嗷叫费观| 波多野结衣一区麻豆| 侵犯人妻中文字幕一二三四区| 18禁美女被吸乳视频| 热re99久久国产66热| 欧美黄色片欧美黄色片| 久久人妻av系列| 美女午夜性视频免费| 色婷婷av一区二区三区视频| 亚洲伊人色综图| 国产av又大| 国产精品98久久久久久宅男小说| 满18在线观看网站| 中文字幕精品免费在线观看视频| 亚洲av日韩精品久久久久久密| 亚洲情色 制服丝袜| 91成人精品电影| 99riav亚洲国产免费| 精品国内亚洲2022精品成人| 免费日韩欧美在线观看| 午夜免费成人在线视频| 欧美日韩乱码在线| 亚洲精品国产色婷婷电影| 超碰97精品在线观看| 午夜福利免费观看在线| 99久久综合精品五月天人人| 水蜜桃什么品种好| 91精品三级在线观看| 桃色一区二区三区在线观看| 精品福利观看| 99国产精品免费福利视频| 美女 人体艺术 gogo| 美女国产高潮福利片在线看| 亚洲男人天堂网一区| 99久久精品国产亚洲精品| 欧美精品亚洲一区二区| 满18在线观看网站| a级毛片黄视频| 久久人妻熟女aⅴ| 久久精品国产综合久久久| 欧美激情极品国产一区二区三区| 可以在线观看毛片的网站| 欧美老熟妇乱子伦牲交| 最新美女视频免费是黄的| 80岁老熟妇乱子伦牲交| 国产激情久久老熟女| 无遮挡黄片免费观看| 热99re8久久精品国产| 国产高清视频在线播放一区| 中文字幕人妻熟女乱码| 999久久久国产精品视频| 国产熟女午夜一区二区三区| 亚洲熟女毛片儿| 中国美女看黄片| 国产乱人伦免费视频| 狠狠狠狠99中文字幕| 久久精品国产综合久久久| 久久久国产精品麻豆| 丁香六月欧美| 成人免费观看视频高清| 精品卡一卡二卡四卡免费| 色综合站精品国产| 91在线观看av| 欧美日韩乱码在线| 色婷婷久久久亚洲欧美| 国产精品久久久久久人妻精品电影| 免费日韩欧美在线观看| 精品福利永久在线观看| 欧美日韩精品网址| 一级片免费观看大全| 亚洲av五月六月丁香网| 亚洲一区二区三区欧美精品| 黄网站色视频无遮挡免费观看| 一级,二级,三级黄色视频| 九色亚洲精品在线播放| a在线观看视频网站| 国产亚洲av高清不卡| 亚洲五月天丁香| 欧美不卡视频在线免费观看 | 妹子高潮喷水视频| 别揉我奶头~嗯~啊~动态视频| 久久婷婷成人综合色麻豆| 在线观看免费日韩欧美大片| 亚洲一区二区三区色噜噜 | 成年版毛片免费区| 免费人成视频x8x8入口观看| 在线观看66精品国产| 1024视频免费在线观看| 成人特级黄色片久久久久久久| 亚洲五月婷婷丁香| 国产99白浆流出| tocl精华| 亚洲专区中文字幕在线| 日本一区二区免费在线视频| av中文乱码字幕在线| 桃色一区二区三区在线观看| 久久国产精品男人的天堂亚洲| 老司机午夜福利在线观看视频| 搡老岳熟女国产| 久久草成人影院| 久久国产精品男人的天堂亚洲| 欧美成人免费av一区二区三区| 午夜福利欧美成人| 国产高清videossex| 国产精品久久久人人做人人爽| 国产片内射在线| 国产精品国产av在线观看| 国产av一区在线观看免费| 1024香蕉在线观看| 国产成人av激情在线播放| 男女床上黄色一级片免费看| 免费一级毛片在线播放高清视频 | 中文字幕人妻熟女乱码| 婷婷精品国产亚洲av在线| 国产亚洲欧美98| 亚洲美女黄片视频| 亚洲国产中文字幕在线视频| 欧美午夜高清在线| 91老司机精品| 成人亚洲精品一区在线观看| 高清欧美精品videossex| 亚洲av熟女| 亚洲午夜精品一区,二区,三区| 啪啪无遮挡十八禁网站| 亚洲三区欧美一区| 亚洲人成电影免费在线| 午夜影院日韩av| 亚洲第一欧美日韩一区二区三区| av视频免费观看在线观看| 亚洲欧洲精品一区二区精品久久久| 在线永久观看黄色视频| 亚洲全国av大片| 久久九九热精品免费| 亚洲成人久久性| 久久久精品国产亚洲av高清涩受| 国产精品久久久av美女十八| www.熟女人妻精品国产| 十八禁网站免费在线| 在线观看免费高清a一片| 琪琪午夜伦伦电影理论片6080| 免费av中文字幕在线| 国产三级在线视频| 大码成人一级视频| 黄色女人牲交| 亚洲国产中文字幕在线视频| 国产又爽黄色视频| 国产精品 欧美亚洲| 日韩欧美三级三区| 国产精品野战在线观看 | av免费在线观看网站| 18禁裸乳无遮挡免费网站照片 | 国产激情欧美一区二区| 免费高清在线观看日韩| 美女国产高潮福利片在线看| 成人手机av| 91字幕亚洲| 热99re8久久精品国产| 精品久久久久久,| www.精华液| 在线免费观看的www视频| 可以在线观看毛片的网站| 久久久精品欧美日韩精品| 色播在线永久视频| 男女之事视频高清在线观看| 精品少妇一区二区三区视频日本电影| 69精品国产乱码久久久| 波多野结衣av一区二区av| 欧洲精品卡2卡3卡4卡5卡区| 岛国视频午夜一区免费看| 真人做人爱边吃奶动态| 亚洲精品中文字幕一二三四区| 久久国产精品影院| 午夜精品久久久久久毛片777| 神马国产精品三级电影在线观看 | 99热只有精品国产| 亚洲av片天天在线观看| 亚洲精品中文字幕一二三四区| 变态另类成人亚洲欧美熟女 | 欧美最黄视频在线播放免费 | 国产精品野战在线观看 | 日本黄色视频三级网站网址| 精品国产一区二区久久| 日韩精品免费视频一区二区三区| 91老司机精品| 丝袜美腿诱惑在线| 一级a爱视频在线免费观看| 青草久久国产| 亚洲伊人色综图| 校园春色视频在线观看| 日本精品一区二区三区蜜桃| 国产成人精品久久二区二区免费| 亚洲七黄色美女视频| 午夜两性在线视频| 美女 人体艺术 gogo| 日韩 欧美 亚洲 中文字幕| 国产精品一区二区在线不卡| 91成年电影在线观看| 成人精品一区二区免费| 女生性感内裤真人,穿戴方法视频| 50天的宝宝边吃奶边哭怎么回事| 精品一品国产午夜福利视频| 亚洲人成伊人成综合网2020| 悠悠久久av| 国产精品成人在线| 午夜两性在线视频| 久久精品成人免费网站| 国产精品99久久99久久久不卡| 欧美在线一区亚洲| 久久久久国产精品人妻aⅴ院| 少妇裸体淫交视频免费看高清 | 99国产精品一区二区三区| 叶爱在线成人免费视频播放| 91九色精品人成在线观看| av在线天堂中文字幕 | 国产xxxxx性猛交| www国产在线视频色| 一级a爱片免费观看的视频| 欧美成人免费av一区二区三区| 久久中文看片网| 国产亚洲精品综合一区在线观看 | 精品一区二区三区视频在线观看免费 | 啦啦啦在线免费观看视频4| 黄频高清免费视频| 精品一区二区三区视频在线观看免费 | 看黄色毛片网站| 交换朋友夫妻互换小说| 欧美一区二区精品小视频在线| 大码成人一级视频| 免费高清视频大片| 色婷婷久久久亚洲欧美| 欧美人与性动交α欧美精品济南到| 精品久久久久久,| 久久影院123| 黑丝袜美女国产一区| 麻豆成人av在线观看| 精品国产一区二区久久| 久久精品国产亚洲av香蕉五月| 国产成人免费无遮挡视频| 69av精品久久久久久| 老司机靠b影院| 久久久久国产一级毛片高清牌| 欧美黑人欧美精品刺激| 超碰成人久久| 亚洲色图 男人天堂 中文字幕| 欧美av亚洲av综合av国产av| 淫妇啪啪啪对白视频| 少妇裸体淫交视频免费看高清 | 免费人成视频x8x8入口观看| ponron亚洲| 老司机福利观看| 久久久水蜜桃国产精品网| 岛国在线观看网站| 国产熟女午夜一区二区三区| 亚洲国产看品久久| av国产精品久久久久影院| 国产有黄有色有爽视频| 夜夜看夜夜爽夜夜摸 | 国产一区二区三区视频了| 国产单亲对白刺激| 99国产精品免费福利视频| 热99国产精品久久久久久7| 99在线视频只有这里精品首页| 啦啦啦免费观看视频1| 亚洲 欧美一区二区三区| 69av精品久久久久久| 日韩免费av在线播放| 法律面前人人平等表现在哪些方面| 国产亚洲欧美在线一区二区| 一级a爱视频在线免费观看| 亚洲人成伊人成综合网2020| 亚洲成人国产一区在线观看| 高清欧美精品videossex| tocl精华| 在线观看一区二区三区| 国产成人精品久久二区二区免费| 日韩高清综合在线| 久久久国产欧美日韩av| 日韩大尺度精品在线看网址 | 国产精品 欧美亚洲| 一个人免费在线观看的高清视频| 国产视频一区二区在线看| 欧美丝袜亚洲另类 | 亚洲五月婷婷丁香| 国产在线观看jvid| 中文字幕av电影在线播放| 日韩人妻精品一区2区三区| 男人舔女人的私密视频| 亚洲欧洲精品一区二区精品久久久| av在线天堂中文字幕 | 午夜福利在线观看吧| 亚洲精品一区av在线观看| 亚洲在线自拍视频| 欧美在线一区亚洲| 99久久人妻综合| xxx96com| 亚洲一码二码三码区别大吗| 嫩草影视91久久| 伦理电影免费视频| 免费女性裸体啪啪无遮挡网站| 99久久久亚洲精品蜜臀av| 亚洲国产中文字幕在线视频| 天天添夜夜摸| 国产一区二区三区视频了| 美女扒开内裤让男人捅视频| 亚洲av日韩精品久久久久久密| 亚洲专区国产一区二区| 首页视频小说图片口味搜索| 久久久水蜜桃国产精品网| 变态另类成人亚洲欧美熟女 | 国产精品乱码一区二三区的特点 | 99久久99久久久精品蜜桃| 91国产中文字幕| 国产精品香港三级国产av潘金莲| 亚洲国产欧美网| 午夜精品在线福利| 久久精品aⅴ一区二区三区四区| 欧美日韩国产mv在线观看视频| 香蕉国产在线看| 日韩欧美一区视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 久9热在线精品视频| 亚洲片人在线观看| 亚洲午夜理论影院| 真人一进一出gif抽搐免费| 男女做爰动态图高潮gif福利片 | 天天躁夜夜躁狠狠躁躁| 伦理电影免费视频| 18禁美女被吸乳视频| 97碰自拍视频| 国产成人av激情在线播放| 夜夜爽天天搞| 色在线成人网|