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

    黑果枸杞SPL轉(zhuǎn)錄因子家族成員鑒定與生物信息學(xué)分析

    2024-12-05 00:00:00劉筠郭毅李效雄胡曉桐馬瑞杜雨賈西貝柳迪馬彥軍

    摘要: SPL(SQUAMOSA promoter-binding protein-like)家族成員參與植物的生長(zhǎng)發(fā)育和脅迫響應(yīng),為探究黑果枸杞SPL家族基因在鹽脅迫下的耐鹽分子機(jī)制,本研究對(duì)其基因家族成員進(jìn)行了鑒定和生物信息學(xué)分析。結(jié)果表明,從黑果枸杞轉(zhuǎn)錄組數(shù)據(jù)中共篩選出20個(gè)SPL基因家族成員,分為8個(gè)亞家族,蛋白質(zhì)亞細(xì)胞定位預(yù)測(cè)結(jié)果顯示SPL蛋白主要定位于細(xì)胞核,20個(gè)SPL蛋白均具有SBP保守結(jié)構(gòu)域。黑果枸杞SPL家族蛋白質(zhì)相對(duì)分子量為15 296.08~130 310.8;等電點(diǎn)為5.64~9.65;總原子數(shù)為2 106~18 214;氨基酸數(shù)量為133~1 176 aa;脂溶指數(shù)為46.35~87.53;均為親水性蛋白質(zhì)和不穩(wěn)定蛋白質(zhì)。黑果枸杞SPL家族蛋白質(zhì)結(jié)構(gòu)分析結(jié)果表明,無(wú)規(guī)則卷曲和α-螺旋占比較大。表達(dá)模式分析結(jié)果表明,在NaCl脅迫下不同器官中黑果枸杞SPL家族基因具有不同的表達(dá)趨勢(shì),LrSPL9在根和葉中均為高表達(dá)。該研究結(jié)果可為黑果枸杞SPL基因家族成員鹽脅迫響應(yīng)方面的進(jìn)一步研究提供參考。

    關(guān)鍵詞: 黑果枸杞;NaCl脅迫;SPL轉(zhuǎn)錄因子;生物信息學(xué)分析

    中圖分類號(hào): S722.3+6 文獻(xiàn)標(biāo)識(shí)碼: A 文章編號(hào): 1000-4440(2024)11-2001-12

    Identification and bioinformatics analysis of SPL transcription factor family members in Lycium ruthenicum

    LIU Yun1, GUO Yi2, LI Xiaoxiong3, HU Xiaotong1, MA Rui1, DU Yu1, JIA Xibei1, LIU Di1, MA Yanjun1

    (1.College of Forestry, Gansu Agricultural University, Lanzhou 730070, China;2.Gansu Forestry Polytechnic, Tianshui 741020, China;3.Lanzhou Resources amp; Environment Voc-Tech University, Lanzhou 730030, China)

    Abstract: The SQUAMOSA promoter-binding protein-like (SPL) family members are involved in plant growth and development and stress response. To investigate the salt tolerance molecular mechanism of Lycium ruthenicum SPL family genes under salt stress, the gene family members were identified and analyzed by bioinformatics. The results showed that 20 SPL family members were screened from transcriptome data of Lycium ruthenicum, and they were divided into eight subfamilies. The results of protein subcellular localization prediction showed that SPL proteins were mainly located in the nucleus, and 20 SPL proteins had SBP conserved domains. The relative molecular weight of SPL family proteins in Lycium ruthenicum was 15 296.08-130 310.8. The isoelectric point was 5.64-9.65, the total number of atoms was 2 106-18 214, the number of amino acids was 133-1 176 aa, and the lipid solubility index was 46.35-87.53. All of them were hydrophilic proteins and unstable proteins. The results of protein structure analysis of SPL family in Lycium ruthenicum showed that random coil and α-helix accounted for a large proportion. The results of expression pattern analysis showed that the SPL family genes of Lycium ruthenicum had different expression trends in different organs under NaCl stress, and LrSPL9 was highly expressed in roots and leaves. The results of this study can provide a reference for further research on salt stress response of SPL gene family members in Lycium ruthenicum.

    Key words: Lycium ruthenicum;NaCl stress;SPL transcription factor;bioinformatics analysis

    SPL(SQUAMOSA promoter-binding protein-like)是一類廣泛存在于植物中的重要的特異性轉(zhuǎn)錄因子,其蛋白質(zhì)均具有SBP結(jié)構(gòu)域[1-2]。SPL轉(zhuǎn)錄因子最初是在金魚(yú)草(Antirrhinum majus L.)中發(fā)現(xiàn)并由此得名[3],SPL在植物生長(zhǎng)發(fā)育過(guò)程中參與花的發(fā)育及調(diào)控、根的發(fā)育、非生物與生物的脅迫應(yīng)答[4-5]。例如在鹽脅迫或干旱脅迫下SPL家族成員會(huì)通過(guò)調(diào)節(jié)基因參與信號(hào)轉(zhuǎn)導(dǎo)的豐度、脯氨酸合成和花青素代謝等來(lái)響應(yīng)鹽脅迫或干旱脅迫[2]。研究人員在對(duì)擬南芥(Arabidopsis thaliana)的研究中發(fā)現(xiàn),沉默miR156靶向的SPL2、SPL9、SPL11基因,可增強(qiáng)擬南芥遭受高溫脅迫、鹽脅迫和干旱脅迫后的恢復(fù)能力[6];水稻(Oryza sativa)中miR156靶向調(diào)節(jié)SPL基因,可增強(qiáng)水稻遭受鹽脅迫、干旱脅迫后的恢復(fù)能力。白樺(Betula platyphylla)中BpSPL6基因轉(zhuǎn)至擬南芥后進(jìn)行鹽脅迫和干旱脅迫試驗(yàn)[7],結(jié)果發(fā)現(xiàn)脅迫發(fā)生后,BpSPL6基因啟動(dòng)子驅(qū)動(dòng)的GUS基因表達(dá)量有所下降。鐵皮石斛(Dendrobium officinale)中DoSTM3基因在NaCl脅迫下其相對(duì)表達(dá)水平與對(duì)照組相比顯著降低[8],推測(cè)SPL基因可能參與植物鹽脅迫和干旱脅迫的響應(yīng)。

    黑果枸杞(Lycium ruthenicum)為茄科(Solanaceae)枸杞屬(Lycium)多年生落葉小灌木。因果實(shí)中含有豐富的花青素、維生素、蔗糖和脂肪等[9],并且可廣泛應(yīng)用于醫(yī)藥、紡織、食品等領(lǐng)域,在植物中有“軟黃金”的美稱[10]。它還是中國(guó)鹽堿地區(qū)和荒漠地區(qū)的重要建群樹(shù)[11-12]。但是現(xiàn)在關(guān)于黑果枸杞耐鹽基因的相關(guān)報(bào)道相對(duì)較少,所以研究其耐鹽分子機(jī)理對(duì)于人工栽培和良種選育具有現(xiàn)實(shí)意義。

    目前已從擬南芥[6]、水稻[13]、小麥(Triticum aestivum)[14]、玉米(Zea mays)[15]、番茄(Solanum lycopersicum)[16]、圓齒野鴉椿(Euscaphis konishii)[17]、紅花(Carthamus tinctorius)[18]、芥藍(lán)(Brassica oleracea var.)[19]、馬尾松(Pinus massoniana)[20]和山羊草(Aegilops tauschii)[21]等多種植物中鑒定出SPL轉(zhuǎn)錄因子家族基因,但在黑果枸杞中尚未有SPL基因的相關(guān)研究報(bào)道。本研究基于不同濃度NaCl在不同脅迫時(shí)間處理下的黑果枸杞轉(zhuǎn)錄組數(shù)據(jù),通過(guò)生物信息學(xué)方法系統(tǒng)鑒定黑果枸杞SPL基因家族成員,并分析其編碼的蛋白質(zhì)理化性質(zhì)、亞細(xì)胞定位、保守結(jié)構(gòu)域、蛋白質(zhì)結(jié)構(gòu)及系統(tǒng)進(jìn)化樹(shù)及表達(dá)模式,為深入研究SPL基因在黑果枸杞抗鹽脅迫方面的作用提供參考依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)材料

    本試驗(yàn)所用植物材料是保存于甘肅農(nóng)業(yè)大學(xué)林學(xué)院組培實(shí)驗(yàn)室的黑果枸杞組培苗,組培苗培養(yǎng)條件及培養(yǎng)基配方見(jiàn)表1。在配置好的培養(yǎng)基中快速繁殖黑果枸杞組培苗,14 d后選取生長(zhǎng)健壯、長(zhǎng)勢(shì)相對(duì)一致的黑果枸杞組培苗煉苗5 d,然后將其移栽到配置好的1/2 Hoagland營(yíng)養(yǎng)液(pH為5.7)中再培養(yǎng)14 d,之后將其移入1/2 Hoagland營(yíng)養(yǎng)液中進(jìn)行鹽脅迫(NaCl)處理[22-23]。鹽脅迫處理方式有3種,分別是不加NaCl、加入50 mmol/L NaCl、加入250 mmol/L NaCl,分別在鹽脅迫處理0 h、1 h和12 h進(jìn)行取樣[24]。對(duì)相同的NaCl處理方式的樣品分別進(jìn)行3個(gè)生物學(xué)重復(fù)。經(jīng)過(guò)NaCl處理后,選擇不同方位生長(zhǎng)良好的主根和完整的葉片取樣,選取的測(cè)試樣品重量約0.1 g,采集樣品后用液氮速凍并保存于-80 ℃冰箱中備用。樣品密封后送至上海歐易生物醫(yī)學(xué)科技有限公司進(jìn)行轉(zhuǎn)錄組測(cè)序。

    1.2 生物信息學(xué)分析

    1.2.1 黑果枸杞SPL基因家族成員的鑒定 在Pfam數(shù)據(jù)庫(kù)(http://pfam.xfam.org/)下載SPL家族(PF03110)[21,25]的隱馬爾可夫模型文件,使用TBtools軟件進(jìn)行數(shù)據(jù)初篩,然后根據(jù)SPL家族蛋白質(zhì)保守結(jié)構(gòu)域SBP,用NCBI(https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi)在線軟件和SMART(https://smart.embl.de/)在線軟件再次篩選,最終得到黑果枸杞SPL基因家族成員。

    1.2.2 黑果枸杞SPL家族蛋白質(zhì)亞細(xì)胞定位預(yù)測(cè)及理化性質(zhì)分析 應(yīng)用WoLF PSORT(https://wolfpsort.hgc.jp/)在線軟件預(yù)測(cè)黑果枸杞SPL家族蛋白質(zhì)亞細(xì)胞定位,采用Expasy-ProtParam(https://web.expasy.org/protparam/)在線軟件對(duì)黑果枸杞SPL家族蛋白質(zhì)理化性質(zhì)和親疏水性進(jìn)行分析,使用SignalP-6.0(https://services.healthtech.dtu.dk/services/SignalP-6.0/)在線軟件預(yù)測(cè)黑果枸杞SPL家族蛋白質(zhì)信號(hào)肽,使用ChiPlot(https://www.chiplot.online/)在線軟件繪制黑果枸杞SPL家族蛋白質(zhì)亞細(xì)胞定位預(yù)測(cè)結(jié)果。

    1.2.3 黑果枸杞SPL家族蛋白質(zhì)結(jié)構(gòu)和保守結(jié)構(gòu)域分析 利用SOPMA(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html)在線軟件進(jìn)行黑果枸杞SPL家族蛋白質(zhì)二級(jí)結(jié)構(gòu)預(yù)測(cè)和分析,采用MEME(https://meme-suite.org/meme/tools/meme)在線軟件分析黑果枸杞SPL家族蛋白質(zhì)保守基序,使用TBtools軟件對(duì)黑果枸杞SPL蛋白的系統(tǒng)進(jìn)化、保守結(jié)構(gòu)域和保守基序進(jìn)行數(shù)據(jù)可視化。

    1.2.4 黑果枸杞SPL家族蛋白質(zhì)系統(tǒng)進(jìn)化分析 在TAIR(https://www.arabidopsis.org/)數(shù)據(jù)庫(kù)下載擬南芥SPL家族蛋白質(zhì)數(shù)據(jù),在NCBI(https://www.ncbi.nlm.nih.gov/)數(shù)據(jù)庫(kù)下載番茄SPL家族蛋白質(zhì)數(shù)據(jù)和煙草(Nicotiana tabacum)SPL家族蛋白質(zhì)數(shù)據(jù),利用MEGA-X軟件的MUSCLE工具進(jìn)行黑果枸杞、番茄、煙草和擬南芥SPL蛋白的多序列比對(duì),然后使用MEGA-X軟件最大似然法(Maximum likelihood,ML)構(gòu)建系統(tǒng)發(fā)育樹(shù),并把Bootstrap參數(shù)設(shè)置為1 000。使用ChiPlot TVBOT(https://www.chiplot.online/tvbot.html)[26]在線軟件繪制黑果枸杞、番茄、煙草和擬南芥SPL家族蛋白質(zhì)的系統(tǒng)進(jìn)化樹(shù)。

    1.2.5 黑果枸杞SPL基因家族成員的表達(dá)模式分析 根據(jù)軟件[27-28]計(jì)算得到的FPKM(Fragments per kilobase per million reads sequenced)數(shù)據(jù),使用ChiPlot在線軟件繪制黑果枸杞SPL基因家族成員在其根與葉中的表達(dá)模式以及黑果枸杞SPL基因響應(yīng)鹽脅迫的表達(dá)分析。

    2 結(jié)果與分析

    2.1 黑果枸杞SPL基因家族成員鑒定

    黑果枸杞轉(zhuǎn)錄組數(shù)據(jù)在NCBI數(shù)據(jù)庫(kù)(https://www.ncbi.nlm.nih.gov/sra)的登錄號(hào)為SRR7700825[29-31]。通過(guò)SPL家族的隱馬爾可夫模型文件和TBtools軟件初篩黑果枸杞SPL基因家族成員后,在NCBI Conserved Domains和SMART上進(jìn)一步去除無(wú)SBP保守結(jié)構(gòu)域的氨基酸序列,最終得到黑果枸杞SPL基因家族成員20個(gè),分別命名為L(zhǎng)rSPL1~LrSPL20(表2)。

    2.2 黑果枸杞SPL家族蛋白質(zhì)亞細(xì)胞定位及理化性質(zhì)

    從圖1可知,黑果枸杞SPL家族蛋白質(zhì)中除LrSPL10蛋白亞細(xì)胞定位于質(zhì)膜,其余19個(gè)蛋白質(zhì)亞細(xì)胞定位在細(xì)胞核,并且20條蛋白質(zhì)氨基酸序列都不含有信號(hào)肽。從表2可見(jiàn),黑果枸杞SPL蛋白相對(duì)分子量為15 296.08~130 310.8;等電點(diǎn)為5.64~9.65,除LrSPL9、LrSPL10、LrSPL11、LrSPL12、LrSPL15蛋白之外,其余15個(gè)LrSPL蛋白的等電點(diǎn)都大于7.0,說(shuō)明黑果枸杞SPL家族蛋白質(zhì)大多數(shù)富含堿性氨基酸;不穩(wěn)定系數(shù)為43.97~78.27,均為不穩(wěn)定蛋白質(zhì);脂溶指數(shù)為46.35~87.53;親水指數(shù)為-1.199~-0.164,全部是親水性蛋白質(zhì);總原子數(shù)為2 106~18 214;氨基酸個(gè)數(shù)為133~1 176 aa。

    Expasy-ProtScale在線軟件分析黑果枸杞SPL家族蛋白質(zhì)的親水性和疏水性,結(jié)果如圖2所示,LrSPL1最高值為2.456,在113位;最低值為-3.322,在152位。LrSPL2最高值為1.767,在192位;最低值為-3.989,在232位。LrSPL3最高值為1.233,在362位;最低值為-3.600,在226位。LrSPL4最高值為2.244,在150位;最低值為-3.356,在244位。LrSPL5最高值為1.789,在23位;最低值為-3.744,在96位。LrSPL6最高值為1.467,在243位;最低值為-3.600,在205位。LrSPL7最高值為1.667,在103位;最低值為-3.422,在142位。LrSPL8最高值為1.678,在186位;最低值為-3.356,在227位。LrSPL9最高值為3.311,在961位;最低值為-3.578,在210位。LrSPL10最高值為3.344,在1149位;最低值為-3.544,在215位。LrSPL11最高值為2.522,在344位;最低值為-4.067,在193位。LrSPL12最高值為3.344,在991位;最低值為-3.544,在215位。LrSPL13最高值為3.311,在951位;最低值為-3.611,在223位。LrSPL14最高值為2.544,在87位;最低值為-3.544,在127位。LrSPL15最高值為2.522,在374位;最低值為-4.067,在193位。LrSPL16最高值為1.667,在38位;最低值為-3.422,在77位。LrSPL17最高值為1.967,在32位;最低值為-3.544,在71位。LrSPL18最高值為1.778,在92位;最低值為-3.633,在278位。LrSPL19最高值為2.122,在242位;最低值為-3.633,在281。LrSPL20最高值為1.678,在268位;最低值為-3.633,在176位。由此可見(jiàn),20條LrSPL蛋白氨基酸序列均存在明確的親水區(qū)域和疏水區(qū)域,并且親水區(qū)域多于疏水區(qū)域,表明LrSPL家族蛋白質(zhì)均為親水性蛋白質(zhì)。

    2.3 黑果枸杞SPL家族蛋白質(zhì)結(jié)構(gòu)和保守結(jié)構(gòu)域

    LrSPL蛋白質(zhì)二級(jí)結(jié)構(gòu)預(yù)測(cè)和分析結(jié)果(表3、圖3)表明,20個(gè)蛋白質(zhì)均含有α-螺旋、β-轉(zhuǎn)角、無(wú)規(guī)則卷曲和延伸鏈。在LrSPL家族蛋白質(zhì)中α-螺旋分布占比11.00%~48.91%,β-轉(zhuǎn)角分布占比0.75%~8.43%,無(wú)規(guī)則卷曲分布占比37.23%~70.72%,延伸鏈分布占比7.30%~18.21%,其中無(wú)規(guī)則卷曲>α-螺旋>延伸鏈>β-轉(zhuǎn)角的LrSPL蛋白最多,有17個(gè),所占比例為85.00%;2個(gè)LrSPL蛋白二級(jí)結(jié)構(gòu)表現(xiàn)為無(wú)規(guī)則卷曲>延伸鏈>α-螺旋>β-轉(zhuǎn)角,比例是10.00%;1個(gè)LrSPL蛋白二級(jí)結(jié)構(gòu)表現(xiàn)為α-螺旋>無(wú)規(guī)則卷曲>延伸鏈>β-轉(zhuǎn)角,比例是5.00%。LrSPL蛋白空間結(jié)構(gòu)的不同決定了蛋白質(zhì)功能的差異,為后續(xù)黑果枸杞SPL家族的蛋白質(zhì)研究提供參考。

    黑果枸杞20個(gè)LrSPL蛋白氨基酸序列中均含有1個(gè)SBP結(jié)構(gòu)域;其中LrSPL9、LrSPL10、LrSPL12和LrSPL13還分別含有1個(gè)Ank_2 superfamily結(jié)構(gòu)域(圖4B)。MEME在線軟件對(duì)LrSPL家族蛋白質(zhì)的分析結(jié)果顯示,LrSPL蛋白共包含5個(gè)保守基序,分別為motif 1、 motif 2、 motif 3、 motif4和 motif 5,其中含5個(gè)motif的LrSPL蛋白有4個(gè),含4個(gè)motif的LrSPL蛋白有2個(gè),含3個(gè)motif的LrSPL蛋白有13個(gè),含2個(gè)motif的LrSPL蛋白有1個(gè),除LrSPL17,其余蛋白質(zhì)均含motif 1、motif 3和motif 4(圖4C)。由此推測(cè),motif 1、motif 3和motif 4相比其他基序的保守性較高,與LrSPL家族蛋白質(zhì)結(jié)構(gòu)域的保守性基本一致。

    2.4 黑果枸杞SPL家族蛋白質(zhì)的系統(tǒng)進(jìn)化

    為了更全面地了解黑果枸杞SPL家族蛋白質(zhì)的系統(tǒng)進(jìn)化關(guān)系,將其20個(gè)SPL蛋白,番茄的29個(gè)SPL蛋白,煙草的13個(gè)SPL蛋白和擬南芥的16個(gè)SPL蛋白,共 78個(gè)SPL蛋白構(gòu)建系統(tǒng)進(jìn)化樹(shù)(圖5)。根據(jù)4種植物的SPL家族蛋白質(zhì)系統(tǒng)進(jìn)化樹(shù)可知,黑果枸杞SPL家族成員可分為8個(gè)亞家族,黑果枸杞中存在3個(gè)同源基因?qū)ΑT赟PL家族蛋白質(zhì)系統(tǒng)進(jìn)化樹(shù)的78個(gè)成員中,第Ⅰ亞家族LrSPL4與煙草Nt59156的親緣關(guān)系最近,LrSPL8與番茄Sl010321004的親緣關(guān)系最近;第Ⅱ亞家族LrSPL7與番茄Sl004249164的親緣關(guān)系最近,LrSPL16與煙草Nt59145的親緣關(guān)系最近;第Ⅲ亞家族LrSPL3、LrSPL6與番茄Sl004252249、Sl025884439的親緣關(guān)系近,LrSPL18和LrSPL20與擬南芥At1G69170的親緣關(guān)系最近,LrSPL19與番茄Sl004239031的親緣關(guān)系最近,LrSPL5與煙草Nt59148的親緣關(guān)系近;第Ⅳ亞家族LrSPL10、LrSPL12、LrSPL13與番茄Sl004239889的親緣關(guān)系近,LrSPL9與煙草Nt59149的親緣關(guān)系近;第V亞家族,LrSPL2與番茄Sl004229775的親緣關(guān)系最近;第Ⅵ亞家族LrSPL1與煙草Nt59154、Nt59153的親緣關(guān)系近,LrSPL17與番茄Sl004243367的親緣關(guān)系最近;第Ⅶ亞家族LrSPL14與煙草Nt59157的親緣關(guān)系最近; 第Ⅷ亞家族LrSPL11、LrSPL15與番茄Sl019070792、Sl019070794、Sl010323003的親緣關(guān)系近。根據(jù)以上結(jié)果,黑果枸杞SPL家族蛋白質(zhì)與番茄SPL家族蛋白質(zhì)親緣關(guān)系最近,其次是煙草。因此推測(cè)LrSPL家族蛋白質(zhì)具有多樣性,而親緣關(guān)系較近的蛋白質(zhì)可能還具有相似或相近的生物學(xué)功能。

    2.5 黑果枸杞SPL基因家族成員在鹽脅迫下的表達(dá)模式

    為初步探索黑果枸杞SPL基因的功能,通過(guò)不同濃度NaCl、不同脅迫時(shí)間處理,研究20個(gè)LrSPL基因在黑果枸杞根與葉中的表達(dá)(圖6)。結(jié)果表明,在鹽脅迫下,不同LrSPL基因在不同器官中的表達(dá)量存在差異,并且SPL基因在黑果枸杞根中比在葉中更活躍。其中LrSPL20在根部和葉片的表達(dá)量較低,而LrSPL9在根部和葉片的表達(dá)量都較高,LrSPL1、LrSPL2、LrSPL10、LrSPL12、LrSPL14和LrSPL17在根部下調(diào),葉片上調(diào);LrSPL3、LrSPL4、LrSPL5、LrSPL6、LrSPL8、LrSPL9、LrSPL11、LrSPL15和LrSPL19在根部上調(diào),葉片下調(diào)。LrSPL1、LrSPL2和LrSPL5在NaCl濃度為50 mmol/L,脅迫時(shí)間為1 h和12 h時(shí)(L50_1、L50_12),在葉片的表達(dá)量與對(duì)照(L0)相比均呈上升趨勢(shì);NaCl濃度為250 mmol/L,脅迫時(shí)間為1 h和12 h時(shí)(L250_1、L250_12),在葉片的表達(dá)量呈下降趨勢(shì);LrSPL9、LrSPL19和LrSPL20在NaCl濃度為50 mmol/L,脅迫時(shí)間為1 h和12 h時(shí)(R50_1、R50_12),在根部的表達(dá)量與對(duì)照(R0)相比均呈上升趨勢(shì);NaCl濃度為250 mmol/L,脅迫時(shí)間為1 h和12 h時(shí)(R250_1、R250_12),在根部的表達(dá)量與對(duì)照相比也呈上升趨勢(shì)。由此推測(cè)LrSPL1、LrSPL2、LrSPL5、LrSPL9、LrSPL19、LrSPL20在響應(yīng)鹽脅迫的過(guò)程中發(fā)揮調(diào)控作用。

    3 討論

    SPL轉(zhuǎn)錄因子是植物中特有的重要轉(zhuǎn)錄因子[32],但是在植物基因家族的鑒定與生物信息學(xué)分析中SPL家族的研究相對(duì)較少。目前已在擬南芥、水稻、玉米、番茄、圓齒野鴉椿、紅花、芥藍(lán)、馬尾松、山羊草、小麥[25]、桂花(Osmanthus fragrans)[33]、甜橙(Citrus sinensis)[34]和甘薯(Dioscorea esculenta)[35]等不同植物中鑒定出SPL基因家族成員,但在不同植物中SPL基因家族的成員數(shù)量存在差異。本研究通過(guò)生物信息學(xué)方法鑒定得到20個(gè)黑果枸杞SPL基因家族成員,而擬南芥有16個(gè)SPL基因,水稻有19個(gè)SPL基因[36],小麥有56個(gè)SPL基因,番茄有29個(gè)SPL基因,圓齒野鴉椿有29個(gè)SPL基因,紅花有21個(gè)SPL基因,芥藍(lán)有31個(gè)SPL基因,馬尾松有11個(gè)SPL基因,山羊草有18個(gè)SPL基因,桂花有29個(gè)SPL基因,甜橙有15個(gè)SPL基因,甘薯有30個(gè)SPL基因。SPL基因家族成員在已完成測(cè)序的植物與黑果枸杞的鑒定結(jié)果中數(shù)量未出現(xiàn)較大差別。

    L50_1表示NaCl濃度為50 mmol/L、脅迫時(shí)間為1 h時(shí)LrSPL基因在葉片的表達(dá);L250_1表示NaCl濃度為250 mmol/L、脅迫時(shí)間為1 h時(shí)LrSPL基因在葉片的表達(dá);L0表示NaCl濃度為0 mmol/L、脅迫時(shí)間為0 h時(shí)LrSPL基因在葉片的表達(dá);L50_12表示NaCl濃度為50 mmol/L、脅迫時(shí)間為12 h時(shí)LrSPL基因在葉片的表達(dá);L250_12表示NaCl濃度為250 mmol/L、脅迫時(shí)間為12 h時(shí)LrSPL基因在葉片的表達(dá)。R50_1表示NaCl濃度為50 mmol/L、脅迫時(shí)間為1 h時(shí)LrSPL基因在根部的表達(dá);R250_1表示NaCl濃度為250 mmol/L、脅迫時(shí)間為1 h時(shí)LrSPL基因在根部的表達(dá);R0表示NaCl濃度為0 mmol/L、脅迫時(shí)間為0 h時(shí)LrSPL基因在根部的表達(dá);R50_12表示NaCl濃度為50 mmol/L、脅迫時(shí)間為12 h時(shí)LrSPL基因在根部的表達(dá);R250_12表示NaCl濃度為250 mmol/L、脅迫時(shí)間為12 h時(shí)LrSPL基因在根部的表達(dá)。

    本研究中LrSPL蛋白家族中多數(shù)成員的等電點(diǎn)大于7.0,這與大部分植物SPL蛋白家族成員的等電點(diǎn)大于7.0的情況基本一致[37-39]。通過(guò)黑果枸杞SPL蛋白保守基序分析結(jié)果可知,除了LrSPL17外,其余19個(gè)蛋白質(zhì)均含motif 1、motif 3和motif 4。這與圓齒野鴉椿、芥藍(lán)、馬尾松、山羊草、甜橙、擬南芥與番茄[35]等SPL蛋白均是motif 1在SPL家族蛋白中占多數(shù)的情況類似。根據(jù)蛋白質(zhì)系統(tǒng)進(jìn)化樹(shù)分析結(jié)果可知,LrSPL家族蛋白質(zhì)可以分為8個(gè)亞族。LrSPL18和LrSPL20聚類于第Ⅲ亞家族的同一分支中,LrSPL10和LrSPL12聚類于第Ⅳ亞家族的同一分支中,LrSPL11和LrSPL15聚類于第Ⅷ亞家族的同一分支中。聚類在同一分支的蛋白質(zhì)其保守結(jié)構(gòu)域的組成、排列順序及位置基本相同,說(shuō)明LrSPL的分類結(jié)果進(jìn)一步得到了保守基序的支持,也意味著LrSPL亞族成員間序列的高度保守性。所有亞族中的黑果枸杞SPL家族成員與番茄SPL家族成員相聚較近,表明黑果枸杞與番茄進(jìn)化距離比擬南芥、煙草更近。

    本研究結(jié)果表明,SPL轉(zhuǎn)錄因子在植物的生長(zhǎng)發(fā)育及脅迫應(yīng)答等過(guò)程中具備調(diào)控作用[32]。本研究鑒定到20個(gè)黑果枸杞SPL家族蛋白質(zhì)均為親水性蛋白質(zhì),在馬尾松、圓齒野鴉椿、鐵皮石斛等植物中,它們的SPL蛋白也均是親水性蛋白質(zhì),而親水性蛋白質(zhì)對(duì)植物抵抗非生物脅迫有利[40]。鹽脅迫下,黑果枸杞根與葉的LrSPL表達(dá)模式分析結(jié)果顯示,其在不同器官中的表達(dá)量存在明顯差異,這與其他植物SPL基因家族成員的鹽脅迫研究結(jié)果有相似之處[35]。在不同濃度鹽脅迫下LrSPL1、LrSPL2和LrSPL5在葉片的表達(dá)量,LrSPL9、LrSPL19和LrSPL20在根部的表達(dá)量,與榮譽(yù)磊等[35]對(duì)擬南芥SPL基因家族成員在鹽脅迫下的研究結(jié)果,楊樂(lè)等[8]對(duì)鐵皮石斛SPL基因家族成員在鹽脅迫下的研究結(jié)果相似,推測(cè)它們可能在響應(yīng)鹽脅迫的過(guò)程中發(fā)揮著調(diào)控作用。LrSPL基因的鑒定與分析結(jié)果為進(jìn)一步研究黑果枸杞在鹽脅迫方面的基因表達(dá)提供理論依據(jù),為進(jìn)一步研究此類基因是否具有鹽脅迫應(yīng)答能力以及轉(zhuǎn)基因品種的獲取提供參考。

    4 結(jié)論

    在黑果枸杞轉(zhuǎn)錄組中共篩選出20個(gè)SPL基因,分為8個(gè)亞家族。蛋白質(zhì)亞細(xì)胞定位預(yù)測(cè)中其家族成員主要定位于細(xì)胞核,蛋白質(zhì)結(jié)構(gòu)預(yù)測(cè)中20個(gè)蛋白質(zhì)均含有α-螺旋、β-轉(zhuǎn)角、無(wú)規(guī)則卷曲和延伸鏈,蛋白質(zhì)系統(tǒng)進(jìn)化樹(shù)顯示在擬南芥、煙草和番茄這3種植物中,黑果枸杞SPL家族蛋白質(zhì)與番茄SPL家族蛋白質(zhì)親緣關(guān)系最近,其次是煙草。表達(dá)模式分析結(jié)果表明,LrSPL20在根部和葉片中均為低表達(dá),LrSPL9在根部和葉片中均為高表達(dá)。LrSPL9、LrSPL19和LrSPL20在NaCl不同濃度,不同脅迫時(shí)間下,與對(duì)照相比在根部的表達(dá)量均呈上升趨勢(shì);LrSPL1、LrSPL2和LrSPL5在低濃度鹽脅迫,不同脅迫時(shí)間下,與對(duì)照相比葉片表達(dá)量呈上升趨勢(shì),而高濃度鹽脅迫下葉片表達(dá)量與對(duì)照相比呈下降趨勢(shì)。由此推測(cè)LrSPL基因的功能可能與鹽脅迫響應(yīng)相關(guān)。本研究結(jié)果為黑果枸杞SPL基因家族成員的功能鑒定及后續(xù)黑果枸杞耐鹽基因的進(jìn)一步篩選提供了參考。

    參考文獻(xiàn):

    [1] YAMASAKI K, KIGAWA T, INOUE M, et al. A novel zinc-binding motif revealed by solution structures of DNA-binding domains of arabidopsis SBP-family transcription factors[J]. Journal of Molecular Biology,2004,337(1):49-63.

    [2] 張 磊. SPL轉(zhuǎn)錄因子研究進(jìn)展[J]. 農(nóng)業(yè)與技術(shù),2022,42(8):25-27.

    [3] KLEIN J, SAEDLER H, HUIJSER P. A new family of DNA binding proteins includes putative transcriptional regulators of the Antirrhinum majus floral meristem identity gene SQUAMOSA[J]. Molecular amp; General Genetics,1996,250(1):7-16.

    [4] 張 麗. 煙草NtabSPL6-2和NtabSPL6-3基因的功能研究[D].西安:西北大學(xué),2019.

    [5] 吳 艷,侯智紅,程 群,等. SPL轉(zhuǎn)錄因子的研究進(jìn)展[J]. 大豆科學(xué),2019,38(2):304-310.

    [6] CARDON G H, HHMANN S, NETTESHEIM K, et al. Functional analysis of the Arabidopsis thaliana SBP-box gene SPL3:a novel gene involved in the floral transition[J]. The Plant Journal,1997,12(2):367-377.

    [7] 李 豆,蘇功博,胡曉晴,等. 白樺BpSPL6基因啟動(dòng)子的克隆及表達(dá)分析[J]. 北京林業(yè)大學(xué)學(xué)報(bào),2022,44(2):1-10.

    [8] 楊 樂(lè),聶 聰,龍小琴,等. 鐵皮石斛SPL膜結(jié)合(STM)轉(zhuǎn)錄因子的全基因組鑒定及表達(dá)分析[J]. 廣西植物,2023,43(6):1041-1050.

    [9] 趙建華,李浩霞,尹 躍, 等. 枸杞酸性轉(zhuǎn)化酶基因的克隆與表達(dá)[J]. 江蘇農(nóng)業(yè)學(xué)報(bào),2015,31(5):1140-1148.

    [10]王 琴,劉鳳蘭,毛金梅,等. ‘黑杞一號(hào)’與野生黑果枸杞的酚類物質(zhì)分析[J]. 西北植物學(xué)報(bào),2019,39(10):1835-1841.

    [11]郭有燕,劉宏軍,孔東升,等. 干旱脅迫對(duì)黑果枸杞幼苗光合特性的影響[J]. 西北植物學(xué)報(bào),2016,36(1):124-130.

    [12]張培通,郭文琦,李春宏. 枸杞在江蘇沿海灘涂鹽堿地的應(yīng)用前景及實(shí)用栽培技術(shù)[J]. 江蘇農(nóng)業(yè)科學(xué),2014,42(3):197-199.

    [13]MIURA K, IKEDA M, MATSUBARA A, et al. OsSPL14 promotes panicle branching and higher grain productivity in rice[J]. Nature Genetics,2010,42(6):545-549.

    [14]劉 霞,張 斌,毛新國(guó),等. 小麥tae-MIR156前體基因的克隆及其靶基因TaSPL17多態(tài)性分析[J]. 遺傳,2014,36(6):592-602.

    [15]HULTQUIST J F, DORWEILER J E. Feminized tassels of maize mop1 and ts1 mutants exhibit altered level of miR156 and specific SBP-box genes[J]. Planta,2008,229(1):99-113.

    [16]SALINAS M, XING S P, HHMANN S, et al. Genomic organization, phylogenetic comparison and differential expression of the SBP-box family of transcription factors in tomato[J]. Planta,2012,35(6):1171-1184.

    [17]孫維紅,盧 玉,劉鈺園,等. 圓齒野鴉椿SPL全基因組的鑒定及表達(dá)[J]. 亞熱帶農(nóng)業(yè)研究,2024,20(2):117-126.

    [18]宋夢(mèng)瑤,趙 樂(lè),董誠(chéng)明,等. 紅花SPL基因家族全基因組鑒定及表達(dá)分析[J]. 江蘇農(nóng)業(yè)科學(xué),2022,50(19):36-43.

    [19]曾東琳,唐為玲,曾家晶,等. 芥藍(lán)SPL基因家族鑒定及表達(dá)分析[J]. 西北植物學(xué)報(bào),2022,42(3):361-370.

    [20]王春昱,范付華. 馬尾松SPL基因家族鑒定及其響應(yīng)低磷脅迫的表達(dá)分析[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào),2023,31(3):509-517.

    [21]楊紫貽,米福貴,唐 芳,等. 山羊草SPL轉(zhuǎn)錄因子基因家族分析[J]. 中國(guó)草地學(xué)報(bào),2021,43(10):9-17.

    [22]杜 雨,劉 筠,胡曉桐,等. 黑果枸杞GRF轉(zhuǎn)錄因子家族鑒定與生物信息學(xué)分析[J]. 草地學(xué)報(bào),2022,30(6):1403-1413.

    [23]賈西貝,李效雄,杜 雨,等. 黑果枸杞XTH轉(zhuǎn)錄因子家族鑒定與生物信息學(xué)分析[J]. 草地學(xué)報(bào),2023,31(10):2916-2924.

    [24]馬彥軍,段慧榮,魏 佳,等. NaCl脅迫下黑果枸杞轉(zhuǎn)錄組測(cè)序分析[J]. 生物技術(shù)通報(bào),2020,36(2):100-109.

    [25]賀金秋,苗敬南,馬 超,等. 小麥SPL家族基因在非生物脅迫下的表達(dá)分析[J]. 麥類作物學(xué)報(bào),2023,43(10):1335-1343.

    [26]XIE J M, CHEN Y R, CAI G J, et al. Tree visualization by one table (tvBOT):a web application for visualizing, modifying and annotating phylogenetic trees[J]. Nucleic Acids Research,2023,51(W1):587-592.

    [27]LANGMEAD B, SALZBERG S L. Fast gapped-read alignment with Bowtie 2[J]. Nat Methods,2012,9(4):357-359.

    [28]ROBERTS A. Ambiguous fragment assignment for high-throughput sequencing experiments[M]. Berkeley:University of California,2013.

    [29]MA Y J, DUAN H R, ZHANG F, et al. Transcriptomic analysis of Lycium ruthenicum Murr. during fruit ripening provides insight into structural and regulatory genes in the anthocyanin biosynthetic pathway[J]. PLoS One,2018,13(12):e0208627.

    [30]TIIKA R J, WEI J, MA R, et al. Identification and expression analysis of the WRKY gene family during different developmental stages in Lycium ruthenicum Murr. fruit[J]. Peer J,2020,8:e10207.

    [31]WEI J, TIIKA R J, MA Y J, et al. Transcriptome-wide identification and analysis of the KT/HAK/KUP family in black goji under NaCl stress[J]. Agronomy Journal,2022,114(4):2069-2080.

    [32]葛 奇,席會(huì)鵬. 植物SPL轉(zhuǎn)錄因子研究進(jìn)展[J]. 安徽農(nóng)業(yè)科學(xué),2023,51(23):25-29,53.

    [33]吳秀怡,田清尹,楊秀蓮,等. 參與調(diào)控桂花花朵開(kāi)放的SPL基因鑒定及表達(dá)分析[J]. 西北植物學(xué)報(bào),2021,41(11):1801-1809.

    [34]楊 杰,陳 蓉,胡文娟,等. 甜橙SPL基因家族的鑒定及其在成花誘導(dǎo)中的表達(dá)分析[J]. 西北植物學(xué)報(bào),2022,42(1):48-56.

    [35]榮譽(yù)磊,周志林,趙冬蘭,等. 甘薯、番茄、擬南芥中SPL轉(zhuǎn)錄因子的生物信息學(xué)分析[J]. 江蘇農(nóng)業(yè)科學(xué),2021,49(20):74-83.

    [36]胡 麗,楊范敏,陳薇蘭,等. 水稻SPL家族轉(zhuǎn)錄因子的生物學(xué)功能研究進(jìn)展[J]. 中國(guó)水稻科學(xué),2024,38(3):223-232.

    [37]張曉紅. 陸地棉開(kāi)花相關(guān)基因的功能研究及調(diào)控分析[D]. 楊凌:西北農(nóng)林科技大學(xué),2016.

    [38]劉 闖. 18個(gè)白樺SPL基因的鑒定及BpSPL8基因的功能分析[D]. 哈爾濱:東北林業(yè)大學(xué),2017.

    [39]祁香寧. 獼猴桃SBP-box轉(zhuǎn)錄因子基因的鑒定、進(jìn)化及其表達(dá)分析[D]. 楊凌:西北農(nóng)林科技大學(xué),2018.

    [40]駱 洋. 茶樹(shù)花青素還原酶基因的原核表達(dá)[D]. 合肥:安徽農(nóng)業(yè)大學(xué),2011.

    (責(zé)任編輯:黃克玲)

    久久亚洲精品不卡| 中文字幕精品亚洲无线码一区| 亚洲电影在线观看av| 亚洲五月天丁香| 日韩欧美在线乱码| 听说在线观看完整版免费高清| 亚洲中文日韩欧美视频| 久久精品影院6| 亚洲欧美日韩高清在线视频| 国产91精品成人一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 久久婷婷人人爽人人干人人爱| 99久久成人亚洲精品观看| 中文字幕人成人乱码亚洲影| 又爽又黄无遮挡网站| 精品人妻1区二区| 99久久精品国产亚洲精品| 男女午夜视频在线观看| 国产精品亚洲av一区麻豆| 搡老岳熟女国产| 国产激情久久老熟女| 欧美午夜高清在线| 日本一本二区三区精品| 免费高清视频大片| 国产激情欧美一区二区| 哪里可以看免费的av片| 成人欧美大片| 最近视频中文字幕2019在线8| 好看av亚洲va欧美ⅴa在| 亚洲欧美一区二区三区黑人| 成年版毛片免费区| www.www免费av| 悠悠久久av| 久久久久性生活片| 精华霜和精华液先用哪个| 欧洲精品卡2卡3卡4卡5卡区| 亚洲,欧美精品.| 精品久久久久久久人妻蜜臀av| 一个人免费在线观看的高清视频| 国产精品影院久久| 国产日本99.免费观看| 18禁黄网站禁片免费观看直播| 国内久久婷婷六月综合欲色啪| 亚洲18禁久久av| 91在线观看av| 久久久久久久精品吃奶| 亚洲成人免费电影在线观看| 麻豆成人午夜福利视频| 久久人人精品亚洲av| 亚洲成人久久性| 老熟妇乱子伦视频在线观看| 国产成人啪精品午夜网站| 亚洲人成电影免费在线| 一夜夜www| 成人三级黄色视频| 中文字幕久久专区| 18禁美女被吸乳视频| 亚洲国产高清在线一区二区三| 亚洲人成伊人成综合网2020| 天天躁日日操中文字幕| 成人午夜高清在线视频| 夜夜爽天天搞| 国产精品 国内视频| 欧美日韩一级在线毛片| 午夜成年电影在线免费观看| 国内精品久久久久精免费| 99热只有精品国产| 亚洲一区二区三区不卡视频| 青草久久国产| 国产精品香港三级国产av潘金莲| 国产精品av视频在线免费观看| 男女下面进入的视频免费午夜| 日韩人妻高清精品专区| 在线a可以看的网站| 久久精品亚洲精品国产色婷小说| 男女那种视频在线观看| 午夜福利成人在线免费观看| 国产综合懂色| 久久性视频一级片| 日韩 欧美 亚洲 中文字幕| 亚洲专区字幕在线| 亚洲欧美一区二区三区黑人| 99国产极品粉嫩在线观看| 婷婷精品国产亚洲av在线| 身体一侧抽搐| 欧美日本亚洲视频在线播放| 国产综合懂色| 久久久久久久精品吃奶| 此物有八面人人有两片| 午夜久久久久精精品| 制服丝袜大香蕉在线| 我的老师免费观看完整版| 国产v大片淫在线免费观看| 亚洲精品美女久久av网站| 国产精品国产高清国产av| 美女 人体艺术 gogo| 1024香蕉在线观看| 黄色 视频免费看| 国产三级在线视频| 国产主播在线观看一区二区| 国产成人福利小说| 午夜视频精品福利| 丰满人妻熟妇乱又伦精品不卡| 亚洲性夜色夜夜综合| 午夜福利在线观看免费完整高清在 | 91老司机精品| 性色av乱码一区二区三区2| 亚洲一区二区三区色噜噜| 亚洲欧美精品综合久久99| 久久久久久久午夜电影| 午夜福利视频1000在线观看| 五月玫瑰六月丁香| 亚洲色图av天堂| 国产日本99.免费观看| 久久精品影院6| 国产精品亚洲一级av第二区| 老司机午夜福利在线观看视频| 女生性感内裤真人,穿戴方法视频| 熟女电影av网| 亚洲avbb在线观看| 叶爱在线成人免费视频播放| 日韩欧美 国产精品| 免费在线观看亚洲国产| 99精品欧美一区二区三区四区| 国产精品98久久久久久宅男小说| 两人在一起打扑克的视频| 国产成人福利小说| 999久久久精品免费观看国产| 成人18禁在线播放| 国产成人精品无人区| 午夜福利18| 免费一级毛片在线播放高清视频| 欧美日韩亚洲国产一区二区在线观看| 成在线人永久免费视频| 老熟妇仑乱视频hdxx| 午夜精品一区二区三区免费看| 免费看十八禁软件| 久久亚洲精品不卡| 国产精品一区二区三区四区免费观看 | 亚洲欧美日韩卡通动漫| 国产精品自产拍在线观看55亚洲| 黄色日韩在线| 精品午夜福利视频在线观看一区| 日本免费a在线| 久久人人精品亚洲av| 黑人操中国人逼视频| 美女 人体艺术 gogo| a级毛片在线看网站| x7x7x7水蜜桃| 最新美女视频免费是黄的| 久久久久国内视频| 久久天堂一区二区三区四区| 一二三四在线观看免费中文在| 欧美日韩亚洲国产一区二区在线观看| 日本免费a在线| 亚洲国产精品sss在线观看| 搡老熟女国产l中国老女人| 99久久精品国产亚洲精品| 舔av片在线| 人人妻人人澡欧美一区二区| 在线国产一区二区在线| 99热6这里只有精品| 国产精品自产拍在线观看55亚洲| 99国产精品一区二区蜜桃av| 两性夫妻黄色片| 午夜激情福利司机影院| 在线观看免费午夜福利视频| 熟女人妻精品中文字幕| 亚洲人与动物交配视频| 国产麻豆成人av免费视频| 男女之事视频高清在线观看| 最新中文字幕久久久久 | 亚洲 欧美 日韩 在线 免费| 免费在线观看日本一区| 亚洲人成伊人成综合网2020| 久久中文字幕一级| 少妇丰满av| 国产亚洲精品久久久久久毛片| 亚洲在线自拍视频| 国产精品美女特级片免费视频播放器 | 久久人妻av系列| 日韩国内少妇激情av| 亚洲专区字幕在线| 国产精品99久久久久久久久| 天天一区二区日本电影三级| 曰老女人黄片| 在线a可以看的网站| 香蕉久久夜色| 欧美国产日韩亚洲一区| 久9热在线精品视频| av在线天堂中文字幕| 国产不卡一卡二| 婷婷六月久久综合丁香| 一二三四社区在线视频社区8| 色播亚洲综合网| 国产成人一区二区三区免费视频网站| 国内毛片毛片毛片毛片毛片| 啦啦啦免费观看视频1| 国产私拍福利视频在线观看| 欧美激情在线99| 人妻久久中文字幕网| 青草久久国产| 国产一区二区三区在线臀色熟女| 搡老岳熟女国产| 国产视频一区二区在线看| 一级毛片精品| 美女高潮的动态| 成人无遮挡网站| 国产一区二区激情短视频| 少妇人妻一区二区三区视频| 久久久久久久午夜电影| 在线观看午夜福利视频| 亚洲精品美女久久久久99蜜臀| 日本成人三级电影网站| 一区二区三区国产精品乱码| 午夜福利成人在线免费观看| 国产精华一区二区三区| 一个人免费在线观看的高清视频| av国产免费在线观看| 精品一区二区三区视频在线观看免费| 成人精品一区二区免费| 欧美极品一区二区三区四区| 毛片女人毛片| 香蕉丝袜av| 日韩欧美三级三区| 亚洲精品一区av在线观看| 精品不卡国产一区二区三区| 在线永久观看黄色视频| 午夜两性在线视频| 日韩欧美在线乱码| 欧美日韩亚洲国产一区二区在线观看| 嫁个100分男人电影在线观看| 国产成人福利小说| 国产一区在线观看成人免费| 日韩中文字幕欧美一区二区| 一级毛片高清免费大全| 色av中文字幕| 最近在线观看免费完整版| 国产精品电影一区二区三区| 久久精品91蜜桃| 久久亚洲精品不卡| 日韩精品青青久久久久久| 一二三四在线观看免费中文在| 日韩欧美在线乱码| 一本精品99久久精品77| 两人在一起打扑克的视频| 俺也久久电影网| 成人鲁丝片一二三区免费| 搞女人的毛片| 欧美乱色亚洲激情| 久久久久久久久免费视频了| 久久久国产成人精品二区| 香蕉国产在线看| 亚洲av熟女| 亚洲欧美激情综合另类| 精品久久久久久久久久久久久| 国产单亲对白刺激| 亚洲国产色片| 国内精品久久久久久久电影| 久久久国产精品麻豆| 国产成人aa在线观看| 色尼玛亚洲综合影院| 亚洲在线观看片| 99久久综合精品五月天人人| av女优亚洲男人天堂 | 99久久综合精品五月天人人| 亚洲人与动物交配视频| 99riav亚洲国产免费| av视频在线观看入口| 成人特级av手机在线观看| 黄片小视频在线播放| av国产免费在线观看| 欧美又色又爽又黄视频| av在线蜜桃| 亚洲成人免费电影在线观看| 亚洲最大成人中文| 99精品在免费线老司机午夜| 亚洲天堂国产精品一区在线| 欧美一级毛片孕妇| 亚洲欧美激情综合另类| 别揉我奶头~嗯~啊~动态视频| 成人无遮挡网站| 日韩av在线大香蕉| 性色av乱码一区二区三区2| 国产精品爽爽va在线观看网站| 美女高潮喷水抽搐中文字幕| 少妇丰满av| 亚洲性夜色夜夜综合| 性欧美人与动物交配| 亚洲一区二区三区不卡视频| 国内揄拍国产精品人妻在线| 久久亚洲精品不卡| av国产免费在线观看| 久久人人精品亚洲av| 午夜福利高清视频| 午夜福利在线观看免费完整高清在 | 岛国视频午夜一区免费看| 草草在线视频免费看| 一区福利在线观看| 色噜噜av男人的天堂激情| 午夜福利免费观看在线| 亚洲国产色片| 99热这里只有是精品50| 成人高潮视频无遮挡免费网站| 国产高清视频在线播放一区| 日韩欧美一区二区三区在线观看| 午夜亚洲福利在线播放| 狂野欧美白嫩少妇大欣赏| 国产亚洲欧美98| 99久久综合精品五月天人人| 色综合欧美亚洲国产小说| 欧美另类亚洲清纯唯美| 亚洲avbb在线观看| 亚洲男人的天堂狠狠| 蜜桃久久精品国产亚洲av| 国产久久久一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 欧洲精品卡2卡3卡4卡5卡区| 国产精品永久免费网站| 国产精品永久免费网站| 精品久久久久久久久久久久久| 亚洲欧洲精品一区二区精品久久久| 久久中文字幕一级| aaaaa片日本免费| 中文亚洲av片在线观看爽| 舔av片在线| 久久久精品大字幕| www日本黄色视频网| 最近在线观看免费完整版| 又黄又爽又免费观看的视频| 亚洲国产欧美人成| 久久久久性生活片| 婷婷精品国产亚洲av| 美女高潮喷水抽搐中文字幕| 精品久久蜜臀av无| ponron亚洲| 久久久久九九精品影院| 国产亚洲精品一区二区www| 国产高清视频在线观看网站| 最新在线观看一区二区三区| 国产极品精品免费视频能看的| 免费看美女性在线毛片视频| 国产精品一区二区免费欧美| 成人av一区二区三区在线看| 91久久精品国产一区二区成人 | www.www免费av| 亚洲av片天天在线观看| av国产免费在线观看| 后天国语完整版免费观看| 一个人看视频在线观看www免费 | 亚洲av电影在线进入| 人妻丰满熟妇av一区二区三区| 久久人人精品亚洲av| 可以在线观看的亚洲视频| 99视频精品全部免费 在线 | 日韩欧美在线乱码| 亚洲激情在线av| 99久国产av精品| 欧美又色又爽又黄视频| 在线观看午夜福利视频| 岛国在线观看网站| 午夜福利免费观看在线| 老熟妇仑乱视频hdxx| 一二三四社区在线视频社区8| 熟女少妇亚洲综合色aaa.| 午夜福利在线在线| 国产高清有码在线观看视频| 亚洲av片天天在线观看| 精品国产美女av久久久久小说| 日韩有码中文字幕| 国产亚洲精品久久久久久毛片| 51午夜福利影视在线观看| 免费观看精品视频网站| 国产精品免费一区二区三区在线| 在线观看日韩欧美| 精品不卡国产一区二区三区| 久久久精品欧美日韩精品| 久久人妻av系列| 午夜福利在线观看免费完整高清在 | 美女大奶头视频| a在线观看视频网站| 搞女人的毛片| 最近视频中文字幕2019在线8| 国产aⅴ精品一区二区三区波| 国产精品综合久久久久久久免费| 嫩草影院精品99| 国产69精品久久久久777片 | 欧美xxxx黑人xx丫x性爽| 成人18禁在线播放| 91av网一区二区| 两人在一起打扑克的视频| 波多野结衣高清无吗| 亚洲欧美日韩高清专用| 免费无遮挡裸体视频| av在线蜜桃| 琪琪午夜伦伦电影理论片6080| 日本一本二区三区精品| 久久这里只有精品中国| 非洲黑人性xxxx精品又粗又长| 精品国产亚洲在线| 国产精品久久久久久人妻精品电影| 欧美一级a爱片免费观看看| 久久精品91蜜桃| 最近最新中文字幕大全免费视频| 99热精品在线国产| 美女午夜性视频免费| 欧美激情在线99| 在线观看舔阴道视频| 亚洲国产看品久久| 一本久久中文字幕| 超碰成人久久| 免费高清视频大片| 精品无人区乱码1区二区| h日本视频在线播放| 国产精品女同一区二区软件 | 欧美绝顶高潮抽搐喷水| 欧美午夜高清在线| 国产97色在线日韩免费| 成人国产一区最新在线观看| a在线观看视频网站| 手机成人av网站| 免费观看人在逋| 脱女人内裤的视频| 啦啦啦免费观看视频1| 两人在一起打扑克的视频| 日韩人妻高清精品专区| 日韩av在线大香蕉| 两性午夜刺激爽爽歪歪视频在线观看| 色综合站精品国产| 国产伦精品一区二区三区四那| 一a级毛片在线观看| 欧美成狂野欧美在线观看| 免费看光身美女| 亚洲九九香蕉| 国产伦精品一区二区三区四那| 岛国视频午夜一区免费看| 久久久久国产一级毛片高清牌| 97碰自拍视频| 欧美xxxx黑人xx丫x性爽| 老司机深夜福利视频在线观看| 琪琪午夜伦伦电影理论片6080| 美女扒开内裤让男人捅视频| 国内少妇人妻偷人精品xxx网站 | 国产精品久久电影中文字幕| 18禁观看日本| 18禁裸乳无遮挡免费网站照片| 天天一区二区日本电影三级| 国产黄a三级三级三级人| 日本在线视频免费播放| 国产爱豆传媒在线观看| 亚洲国产欧美一区二区综合| 中文字幕av在线有码专区| 88av欧美| 国产亚洲欧美98| 久久精品人妻少妇| tocl精华| 老汉色∧v一级毛片| 久久精品91无色码中文字幕| 亚洲色图av天堂| 国产三级黄色录像| 欧美一级a爱片免费观看看| 99久久精品热视频| 成人亚洲精品av一区二区| 狂野欧美白嫩少妇大欣赏| 国产一区二区三区在线臀色熟女| 久久人妻av系列| 久久久成人免费电影| 特大巨黑吊av在线直播| 女人高潮潮喷娇喘18禁视频| 一进一出抽搐动态| 国产99白浆流出| 男人和女人高潮做爰伦理| 小说图片视频综合网站| 在线观看午夜福利视频| 精品不卡国产一区二区三区| 嫁个100分男人电影在线观看| 麻豆成人av在线观看| 天天躁日日操中文字幕| 久久国产精品影院| www.www免费av| 亚洲中文字幕日韩| 黄片大片在线免费观看| 黑人巨大精品欧美一区二区mp4| 亚洲精华国产精华精| 国内精品美女久久久久久| 国产成人av教育| 色噜噜av男人的天堂激情| 99热这里只有是精品50| 国产精品女同一区二区软件 | 91久久精品国产一区二区成人 | 天天躁狠狠躁夜夜躁狠狠躁| 国产久久久一区二区三区| 亚洲18禁久久av| 香蕉久久夜色| 99riav亚洲国产免费| 亚洲 国产 在线| 欧美日韩亚洲国产一区二区在线观看| 十八禁网站免费在线| 一本综合久久免费| 男插女下体视频免费在线播放| 久久久久亚洲av毛片大全| www.精华液| 国产av不卡久久| 欧美日本视频| 又大又爽又粗| 午夜免费成人在线视频| av黄色大香蕉| 亚洲性夜色夜夜综合| 又紧又爽又黄一区二区| 18禁黄网站禁片午夜丰满| 少妇丰满av| 亚洲 欧美一区二区三区| а√天堂www在线а√下载| 国产乱人伦免费视频| 欧美成人一区二区免费高清观看 | 日本 欧美在线| 欧美不卡视频在线免费观看| 亚洲国产精品成人综合色| av福利片在线观看| 日本 欧美在线| 又紧又爽又黄一区二区| 中国美女看黄片| 激情在线观看视频在线高清| 69av精品久久久久久| 少妇熟女aⅴ在线视频| 69av精品久久久久久| 免费av不卡在线播放| 欧美xxxx黑人xx丫x性爽| 九九在线视频观看精品| 丝袜人妻中文字幕| 动漫黄色视频在线观看| 久久午夜亚洲精品久久| 成人av一区二区三区在线看| 女警被强在线播放| 精品人妻1区二区| 色视频www国产| www日本黄色视频网| 国产av一区在线观看免费| 在线免费观看的www视频| 精品国产乱子伦一区二区三区| 女生性感内裤真人,穿戴方法视频| 少妇的逼水好多| 99国产综合亚洲精品| 一个人看视频在线观看www免费 | 嫁个100分男人电影在线观看| 日本黄大片高清| 美女 人体艺术 gogo| 亚洲片人在线观看| 天堂av国产一区二区熟女人妻| 天天躁日日操中文字幕| 午夜福利在线在线| 观看免费一级毛片| 一个人免费在线观看电影 | 色视频www国产| 亚洲av免费在线观看| 免费看光身美女| 99精品久久久久人妻精品| 在线观看一区二区三区| 窝窝影院91人妻| 嫩草影院入口| av天堂在线播放| 男女那种视频在线观看| 制服丝袜大香蕉在线| 我的老师免费观看完整版| 99久久综合精品五月天人人| 亚洲精品国产精品久久久不卡| 国产乱人伦免费视频| 国内精品久久久久久久电影| 国产精品,欧美在线| 一个人看视频在线观看www免费 | 看片在线看免费视频| 麻豆av在线久日| 人人妻人人澡欧美一区二区| 最新美女视频免费是黄的| 18禁美女被吸乳视频| 亚洲成人久久性| 午夜成年电影在线免费观看| 黄色视频,在线免费观看| 国产一区二区三区视频了| 色av中文字幕| 亚洲无线在线观看| 在线观看66精品国产| 首页视频小说图片口味搜索| 精品乱码久久久久久99久播| 亚洲专区字幕在线| 好男人在线观看高清免费视频| 久久亚洲精品不卡| 好看av亚洲va欧美ⅴa在| 十八禁人妻一区二区| 18禁黄网站禁片免费观看直播| 国产激情偷乱视频一区二区| cao死你这个sao货| 日本五十路高清| 男女那种视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 日韩免费av在线播放| 欧美中文综合在线视频| 久久中文字幕人妻熟女| 麻豆久久精品国产亚洲av| 欧美国产日韩亚洲一区| 少妇裸体淫交视频免费看高清| 国产高潮美女av| 亚洲天堂国产精品一区在线| 51午夜福利影视在线观看| 欧美+亚洲+日韩+国产| 最近最新中文字幕大全电影3| 亚洲av熟女| 国产午夜福利久久久久久| 欧美一级毛片孕妇| 国产精品九九99| 亚洲人成网站在线播放欧美日韩| 国产单亲对白刺激| 麻豆国产97在线/欧美| 全区人妻精品视频| 美女扒开内裤让男人捅视频|