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

    棉花EXO70基因家族全基因組的鑒定及種間比較

    2023-12-28 02:59:24董琰鈺徐碧玉董澤宇汪露瑤陳錦文方磊
    中國農(nóng)業(yè)科學(xué) 2023年23期
    關(guān)鍵詞:海島棉同源分支

    董琰鈺,徐碧玉,董澤宇,汪露瑤,陳錦文,方磊

    棉花EXO70基因家族全基因組的鑒定及種間比較

    董琰鈺,徐碧玉,董澤宇,汪露瑤,陳錦文,方磊

    浙江大學(xué)農(nóng)業(yè)與生物技術(shù)學(xué)院,杭州 310058

    【目的】當(dāng)前四倍體棉花的主栽品種是陸地棉與海島棉,陸地棉產(chǎn)量高且適應(yīng)性好,海島棉纖維品質(zhì)優(yōu)異但產(chǎn)量一般。EXO70是胞外分泌體(exocyst complex)中的關(guān)鍵亞基,在植物生長發(fā)育、脅迫響應(yīng)等方面發(fā)揮重要作用。通過全基因組水平鑒定和分析陸地棉與海島棉EXO70基因家族成員,并研究其在纖維發(fā)育、環(huán)境適應(yīng)等方面的功能,有助于揭示陸地棉與海島棉種間性狀差異形成的分子基礎(chǔ)。【方法】從擬南芥數(shù)據(jù)庫(TAIR)獲取EXO70家族蛋白序列作為參考序列,采用HMMER、ExPASy、MEME、TBtools等分析工具對(duì)陸地棉TM-1和海島棉Hai7124基因組中EXO70家族成員進(jìn)行鑒定與分析,系統(tǒng)比較該家族在基因表達(dá)模式、重要農(nóng)藝性狀的相關(guān)性、逆境響應(yīng)等方面的異同點(diǎn)?!窘Y(jié)果】通過對(duì)陸地棉和海島棉的基因組水平分析,均鑒定出54個(gè)EXO70家族成員,根據(jù)系統(tǒng)進(jìn)化發(fā)育可將其分為8個(gè)分支。陸地棉和海島棉同源基因可一一對(duì)應(yīng),分布在陸地棉與海島棉的20條染色體上,絕大多數(shù)成員為單外顯子,但是兩者間有12對(duì)同源基因在閱讀框序列上存在較大差異。陸地棉與海島棉大部分同源基因表達(dá)模式類似,但在同一纖維發(fā)育時(shí)期表達(dá)量存在差異,如分支A中的與同源基因。重要農(nóng)藝性狀關(guān)聯(lián)分析中發(fā)現(xiàn),海島棉有16個(gè)EXO70基因與纖維品質(zhì)性狀顯著相關(guān),多于陸地棉;陸地棉有13個(gè)基因與產(chǎn)量性狀衣分顯著相關(guān),多于海島棉。發(fā)現(xiàn)部分性狀關(guān)聯(lián)的差異可能由海陸種間的基因結(jié)構(gòu)差異造成,如和。在不同逆境脅迫下,陸地棉中有25個(gè)基因受誘導(dǎo),顯著多于海島棉中10個(gè)受誘導(dǎo)基因。【結(jié)論】在棉花四倍體形成、分化與馴化過程中,陸地棉與海島棉EXO70家族基因序列結(jié)構(gòu)相對(duì)保守,基因表達(dá)模式也保持較高的一致性。海島棉有更多的EXO70基因與纖維品質(zhì)性狀相關(guān),陸地棉有更多的基因與產(chǎn)量性狀相關(guān)且對(duì)環(huán)境逆境更為敏感。

    陸地棉;海島棉;EXO70基因家族;表達(dá)模式;性狀關(guān)聯(lián)分析;逆境響應(yīng)

    0 引言

    【研究意義】棉花是世界上最重要的經(jīng)濟(jì)作物之一,陸地棉()與海島棉()是目前最主要的2個(gè)棉花栽培種。陸地棉產(chǎn)量高、適應(yīng)性強(qiáng),而海島棉產(chǎn)量低但纖維品質(zhì)較好,因此,研究陸地棉適應(yīng)性強(qiáng)和海島棉纖維品質(zhì)好的相關(guān)基因?qū)τ谂嘤呖箖?yōu)質(zhì)的棉花新品種具有重要意義。EXO70(exocyst component of 70 kDa)是控制囊泡運(yùn)輸?shù)闹匾騕1],參與細(xì)胞壁形成、細(xì)胞分泌和抗逆應(yīng)答等重要細(xì)胞學(xué)過程。鑒定陸地棉和海島棉的EXO70家族基因,可為揭示EXO70在棉花生長發(fā)育尤其是纖維發(fā)育中的作用提供基礎(chǔ)數(shù)據(jù),進(jìn)行該類基因與重要農(nóng)藝性狀的關(guān)聯(lián)分析為挖掘關(guān)鍵EXO70基因找到切入點(diǎn),分析和對(duì)比該類基因在陸地棉、海島棉中的表達(dá)模式,將為深入探究EXO70調(diào)控機(jī)制奠定參考基礎(chǔ)?!厩叭搜芯窟M(jìn)展】囊泡運(yùn)輸是真核細(xì)胞中一種關(guān)鍵的跨膜運(yùn)輸和物質(zhì)交流方式,用于轉(zhuǎn)運(yùn)大分子等無法直接穿越細(xì)胞膜和質(zhì)膜的物質(zhì)[2],為細(xì)胞的分裂、生長、形態(tài)建成和極性建立等基本的生理活動(dòng)運(yùn)送原料或傳遞信號(hào)[3]。經(jīng)典的囊泡運(yùn)輸可劃分為4個(gè)過程:出芽(小泡形成)、轉(zhuǎn)運(yùn)(沿細(xì)胞骨架運(yùn)輸)、拴系(兩膜最初接觸)、膜融合(運(yùn)輸完成),由拴系因子介導(dǎo)的拴系過程是膜融合發(fā)生的必要條件[4]。Exocyst復(fù)合體是重要的拴系因子之一,是由SEC3(Secretory3)、SEC5、SEC6、SEC8、SEC10、SEC15、EXO70和EXO84(exocyst component of 84 kDa)組成的八亞基復(fù)合體[5-8]。EXO70是Exocyst復(fù)合體的關(guān)鍵亞基[9],決定胞吐時(shí)的靶標(biāo)膜位點(diǎn)。近年來,越來越多植物中EXO70基因家族成員的功能研究被報(bào)道,主要集中在調(diào)控花粉管伸長和極化、根毛生長、細(xì)胞壁物質(zhì)沉積、細(xì)胞板激活和成熟、防御和自噬方面[10-15],廣泛參與植物生長發(fā)育和免疫反應(yīng)。在植物生長方面,擬南芥影響細(xì)胞板的組裝、導(dǎo)管組織分化及運(yùn)輸、凱氏帶及次生壁的形成[11, 16-17];和參與調(diào)節(jié)花粉管的極性生長與成熟[1, 18];突變會(huì)影響水稻對(duì)礦質(zhì)營養(yǎng)素的吸收及維管束的分化[19]。在免疫調(diào)節(jié)方面,擬南芥EXO70B1和EXO70B2蛋白可以相互作用調(diào)節(jié)鞭毛受體蛋白(FLS2)向質(zhì)膜的運(yùn)輸,從而參與植物的免疫調(diào)節(jié)[12]。水稻與AVR-Pii互作,進(jìn)而參與調(diào)節(jié)免疫過程[20];沉默煙草中的2個(gè)EXO70B基因可以增加植株對(duì)晚疫病菌的感病性[21]。在不同逆境脅迫條件下,葡萄中14個(gè)EXO70基因的表達(dá)有顯著差異,這些基因在抗逆調(diào)控中起重要作用[22]。目前,對(duì)于棉花EXO70家族的報(bào)道主要集中在抗逆與免疫調(diào)節(jié)方面。在棉花應(yīng)對(duì)鹽脅迫過程中發(fā)揮重要作用,并且還與自噬相關(guān)[23]。2021年,Zhu等[24]在亞洲棉、雷蒙德氏棉、陸地棉和海島棉共鑒定到165個(gè)EXO70基因,系統(tǒng)發(fā)育分析表明擬南芥與棉花這兩種雙子葉植物中EXO70基因相似。通過VIGS技術(shù)沉默陸地棉驗(yàn)證基因功能,揭示了棉花EXO70基因可能廣泛參與代謝,也可能影響植物激素的生物合成和分布[24]。【本研究切入點(diǎn)】雖然關(guān)于棉花EXO70基因家族及部分棉花EXO70基因的功能已有相關(guān)報(bào)道,但還未見針對(duì)陸地棉與海島棉EXO70家族成員的組織表達(dá)及關(guān)聯(lián)重要農(nóng)藝性狀的海陸種間比較研究,且EXO70家族成員對(duì)棉花纖維發(fā)育的影響尚不明確?!緮M解決的關(guān)鍵問題】本研究基于陸地棉和海島棉基因組數(shù)據(jù)和二者不同組織及逆境脅迫的轉(zhuǎn)錄組數(shù)據(jù),利用生物信息學(xué)方法進(jìn)行海島棉和陸地棉的EXO70基因家族全基因組鑒定和種間比較分析等工作,旨在揭示EXO70基因家族在生長發(fā)育過程和抗逆境脅迫過程中的作用,結(jié)合陸地棉與海島棉EXO70基因各自優(yōu)勢可進(jìn)一步應(yīng)用于棉花選育與改良方面。

    1 材料與方法

    1.1 棉花EXO70基因家族成員的全基因組鑒定

    為了鑒定EXO70基因家族成員,從Cottonomics數(shù)據(jù)庫(http://cotton.zju.edu.cn/)[25]中下載異源四倍體陸地棉TM-1(v2.1)、異源四倍體海島棉Hai7124(v1.1)基因組序列及基因注釋文件(general feature format,GFF),從Pfam[26]數(shù)據(jù)庫(http://pfam.xfam.org/)中下載EXO70蛋白結(jié)構(gòu)域的隱馬爾可夫模型(編號(hào)PF3081),從擬南芥基因組網(wǎng)站(https://www.arabidopsis. org/)下載已報(bào)道的擬南芥EXO70蛋白的氨基酸序列。利用HMMER 3.0軟件的hmmsearch搜索程序鑒定棉花全基因組中含有EXO70保守結(jié)構(gòu)域的蛋白序列,同時(shí),利用BLAST+v2.6.0[27]搜索棉花全基因組中與已報(bào)道的擬南芥EXO70蛋白的氨基酸序列最相近的蛋白序列(值<1×10-5),將所獲得的EXO70蛋白序列提交至NCBI網(wǎng)站的保守結(jié)構(gòu)域數(shù)據(jù)庫(conserved domain database,CDD)進(jìn)行保守結(jié)構(gòu)域驗(yàn)證。

    1.2 棉花EXO70基因家族成員蛋白序列的理化性質(zhì)分析

    根據(jù)陸地棉與海島棉注釋文件,通過ExPASy在線網(wǎng)站(http://cn.expasy.org/tools)對(duì)GhEXO70和GbEXO70蛋白序列的等電點(diǎn)(isoelectric point,)、相對(duì)分子質(zhì)量(relative molecular mass,Mr)、親水性平均系數(shù)(grand average of hydropathicity,GRAVY)等理化性質(zhì)進(jìn)行預(yù)測。

    1.3 棉花EXO70基因系統(tǒng)發(fā)育樹的構(gòu)建及保守基序分析

    利用MEGA 11.0軟件中的MUSCLE對(duì)擬南芥、海島棉、陸地棉的EXO70家族成員的氨基酸序列進(jìn)行多重序列比對(duì),采用最大似然法自動(dòng)選擇最佳替代模型,自展值(bootstrap value)設(shè)為1 000構(gòu)建系統(tǒng)進(jìn)化樹,利用iTOL(https://itol.embl.de/)對(duì)結(jié)果進(jìn)行可視化分析。通過MEME程序[28]對(duì)家族成員進(jìn)行基序分析(最大發(fā)現(xiàn)數(shù)設(shè)為20)。通過陸地棉與海島棉基因組注釋文件,利用TBtools[29]軟件對(duì)EXO70基因的基因結(jié)構(gòu)及保守基序進(jìn)行數(shù)據(jù)可視化處理。

    1.4 棉花EXO70基因家族成員種內(nèi)共線性分析

    通過MCScanX[30]軟件對(duì)海島棉和陸地棉中EXO70基因進(jìn)行共線性分析,鑒定基因復(fù)制事件,用Circos[31]對(duì)共線性分析結(jié)果進(jìn)行可視化,參數(shù)均使用軟件默認(rèn)參數(shù)。利用TBTools軟件分別計(jì)算陸地棉和海島棉復(fù)制基因?qū)Φ姆峭x突變率(non-synonymous mutation rate,Ka)、同義突變率(synonymous mutation rate,Ks)和二者的比值(Ka/Ks)。

    1.5 棉花EXO70基因家族不同組織中的表達(dá)模式分析

    從Cottonomics數(shù)據(jù)庫和NCBI數(shù)據(jù)庫中下載海島棉與陸地棉不同組織、不同纖維發(fā)育時(shí)期的轉(zhuǎn)錄組原始數(shù)據(jù),數(shù)據(jù)編號(hào)為PRJNA490626。使用軟件fastp[32]過濾以獲得高質(zhì)量讀長,使用HISAT2(v2.1)[33]將高質(zhì)量讀長比對(duì)到TM-1參考基因組或Hai7124參考基因組上,使用featureCounts 2.0.0[34]對(duì)基因表達(dá)進(jìn)行定量,計(jì)算EXO70基因的表達(dá)分度FPKM值以評(píng)估基因的表達(dá)水平。取陸地棉和海島棉葉、根、莖、花瓣及以開花后天數(shù)(day post-anthesis,DPA)計(jì)算的纖維發(fā)育時(shí)期(0、1、3、5、10、20和25 DPA)的EXO70家族基因FPKM數(shù)據(jù),利用R語言做組織表達(dá)模式可視化分析。

    1.6 棉花EXO70基因家族主要纖維發(fā)育時(shí)期基因表達(dá)分析

    采用陌凡生物公司EASYspin PLUS植物RNA快速提取試劑盒,提取2021年種植在杭州市永安基地的陸地棉TM-1的0 DPA、3 DPA、5 DPA和10 DPA的RNA,分別取3 μl總RNA用諾唯贊公司的HiScript Ⅱ 1st Strand cDNA Synthesis Kit(+gDNA wiper)進(jìn)行逆轉(zhuǎn)錄,將反轉(zhuǎn)后的cDNA稀釋為100—200 ng·μl-1。選擇作為內(nèi)參基因,隨機(jī)挑選4個(gè)基因(、、和)進(jìn)行qRT-PCR。反應(yīng)體系選用ChamQ SYBR qPCR Master Mix(諾唯贊)體系。每個(gè)樣品進(jìn)行3次重復(fù)試驗(yàn),用2-ΔΔCt方法計(jì)算該基因的相對(duì)表達(dá)量。

    1.7 棉花EXO70基因家族成員關(guān)于重要農(nóng)藝性狀的關(guān)聯(lián)分析

    下載陸地棉群體和海島棉群體遺傳變異信息和對(duì)應(yīng)的表型數(shù)據(jù)[35-37],利用ANNOVAR軟件[38]分別對(duì)遺傳變異位點(diǎn)進(jìn)行注釋,隨后基于非同義突變的分型對(duì)EXO70基因進(jìn)行單倍型分析。EXO70基因的不同單倍型的纖維品質(zhì)和產(chǎn)量相關(guān)的重要農(nóng)藝性狀(即纖維長度、纖維強(qiáng)度和衣分)進(jìn)行差異分析(<0.05),利用R語言可視化。

    1.8 棉花EXO70基因家族逆境脅迫下表達(dá)模式分析

    從Cottonomics數(shù)據(jù)庫和NCBI數(shù)據(jù)庫中下載海島棉與陸地棉不同逆境脅迫處理的轉(zhuǎn)錄組原始數(shù)據(jù)(編號(hào)PRJNA490626)。原始數(shù)據(jù)的處理和基因表達(dá)量的計(jì)算方法同1.5,以FPKM=1作為條件進(jìn)行過濾,即FPKM>1的基因視為表達(dá),若所有組織中FPKM<1的基因視為不表達(dá),不表達(dá)的基因過濾不進(jìn)行差異基因分析。基于基因轉(zhuǎn)移格式(gene transfer format,GTF)的注釋文件和DESeq2軟件包對(duì)陸地棉和海島棉應(yīng)對(duì)逆境脅迫差異表達(dá)基因(differentially expressed genes,DEGs)進(jìn)行鑒定并篩選,差異基因的篩選標(biāo)準(zhǔn)為取偽發(fā)現(xiàn)率(false discovery rate,F(xiàn)DR)矯正后的<0.05。對(duì)應(yīng)逆境脅迫條件下篩選到的差異基因也即受到脅迫誘導(dǎo)的基因記為“1”,其他記為“0”,并用R語言可視化標(biāo)記EXO70家族基因是否受到脅迫誘導(dǎo)。

    2 結(jié)果

    2.1 棉花EXO70基因家族成員的鑒定和分布

    通過序列相似性與保守結(jié)構(gòu)域搜索在陸地棉TM-1和海島棉Hai7124基因組中均鑒定到了54個(gè)EXO70家族基因。以擬南芥、陸地棉和海島棉EXO70氨基酸序列構(gòu)建棉花EXO70家族成員系統(tǒng)發(fā)育樹(圖1),所有EXO70可分為8個(gè)分支,每個(gè)分支的成員數(shù)量不等但各分支內(nèi)海島棉和陸地棉EXO70基因數(shù)量相同,分支A、B、C、D、E、F、G和H中的陸地棉成員數(shù)量分別為4、2、9、4、9、4、7和15個(gè)。從鑒定家族成員與系統(tǒng)進(jìn)化發(fā)育來看,陸地棉與海島棉在EXO70家族的每個(gè)分支數(shù)量保持一致,并且陸地棉與海島棉同源基因能夠一一對(duì)應(yīng),表明EXO70基因家族在海島棉與陸地棉分化過程中相對(duì)保守。

    圖1 擬南芥(At)、陸地棉(Gh)和海島棉(Gb)EXO70基因家族的系統(tǒng)發(fā)育樹

    2.2 棉花EXO70家族成員的蛋白質(zhì)理化性質(zhì)和結(jié)構(gòu)特征

    通過對(duì)陸地棉EXO70家族成員編碼的蛋白質(zhì)進(jìn)行理化性質(zhì)分析(電子附表1)。EXO70家族成員編碼的肽鏈由310—743個(gè)氨基酸組成,相對(duì)分子質(zhì)量介于34 871.98—86 021.92 Da,平均值為71 029.997 Da;理論等電點(diǎn)介于4.84—8.94,平均值為6.22,多數(shù)偏酸性;蛋白不穩(wěn)定指數(shù)介于37.42—61.76,平均值為46.32,其中8個(gè)EXO70成員為穩(wěn)定蛋白(蛋白不穩(wěn)定指數(shù)<40),其余的46個(gè)基因?yàn)椴环€(wěn)定蛋白(蛋白不穩(wěn)定指數(shù)>40)。亞細(xì)胞定位預(yù)測顯示,EXO70家族基因主要定位于葉綠體、細(xì)胞核、細(xì)胞質(zhì)以及線粒體中,其中,有48個(gè)基因定位在細(xì)胞核中,22個(gè)定位在葉綠體上,9個(gè)定位在細(xì)胞質(zhì)中,4個(gè)定位在線粒體上,以細(xì)胞核中分布最廣泛。海島棉家族成員理化性質(zhì)(電子附表2)與陸地棉類似,無明顯差異。

    2.3 棉花EXO70基因家族保守基序和基因結(jié)構(gòu)分析及比較

    通過對(duì)陸地棉與海島棉EXO70家族的保守基序、基因結(jié)構(gòu)進(jìn)行查詢和比較后發(fā)現(xiàn),同一分支中的海陸同源基因含有相似的保守基序和外顯子結(jié)構(gòu)。保守基序分析(圖2)將預(yù)測到的20種基序依次命名為motif 1—motif 20。同一分支的蛋白序列呈現(xiàn)相似的保守基序排列方式,不同EXO70蛋白所包含的保守基序數(shù)目及種類存在差異。陸地棉與海島棉EXO70蛋白中,基序種類最少的蛋白僅含8種,而最多的含有19種motif,說明家族中的一些成員在保守基序上存在差異。陸地棉和海島棉除了分支A的全部成員與分支G的個(gè)別成員,其他成員都為單外顯子基因。從基因長度來看,不同的基因具有相似長度的編碼區(qū)序列,說明EXO70基因在進(jìn)化過程中內(nèi)含子長度存在較大的變異,這可能是導(dǎo)致其功能多樣化的原因之一。

    通過對(duì)比陸地棉和海島棉同源基因全長序列與編碼序列,發(fā)現(xiàn)有12對(duì)陸地棉海島棉同源基因編碼序列有差異(電子附表3)。其中,編碼序列差異較大的有6對(duì)同源基因集中在分支H中,如的長度小于陸地棉同源基因,且氨基酸數(shù)目比后者少83個(gè)。此外,分支A中的海島棉的編碼序列比它在陸地棉中的同源基因長1 035 bp,相應(yīng)氨基酸數(shù)目多345個(gè)。分析對(duì)比外顯子個(gè)數(shù)可知,在54對(duì)陸地棉海島棉同源基因中,45對(duì)同源基因?yàn)閱瓮怙@子,6對(duì)多外顯子,其余3對(duì)同源基因外顯子個(gè)數(shù)存在差異。如有2個(gè)外顯子,而海島棉同源基因?yàn)閱瓮怙@子基因。

    2.4 棉花EXO70基因家族成員的共線性分析

    對(duì)陸地棉和海島棉中EXO70基因分別進(jìn)行種內(nèi)共線性分析,對(duì)該家族成員進(jìn)行復(fù)制基因比對(duì)。在陸地棉中共檢測出108個(gè)同源基因?qū)Γ▓D3),其中,A、D亞組內(nèi)分別有23對(duì)、20對(duì)直系同源基因,A、D亞組間有65對(duì)旁系同源基因,這些同源基因?qū)Φ腒a/Ks值范圍在0.04—1.05(電子附表4)。而同源基因?qū)χ兄挥信c的Ka/Ks值大于1,為1.043。有6個(gè)同源基因?qū)Πl(fā)生同義突變,有22個(gè)基因沒有找到復(fù)制的同源基因?qū)?。在海島棉中共檢測出104個(gè)同源基因?qū)Γ▓D3),其中,A亞組內(nèi)有24對(duì)同源基因,D亞組內(nèi)有18對(duì),A、D亞組間有62對(duì)旁系同源基因。這些同源基因?qū)Φ腒a/Ks值范圍在0.03—1.05(電子附表5),也只有與基因?qū)Φ腒a/Ks值>1,為1.047。其中,有2個(gè)同源基因?qū)Πl(fā)生同義突變,有21個(gè)基因沒有找到復(fù)制的同源基因?qū)?。從圖3可知,陸地棉除了A亞組的第2、6和8染色體與D亞組的第2、6和8染色體沒有EXO70家族成員外,其他每條染色體上都至少含有一個(gè)EXO70基因,其中,A亞組第5染色體上最多,包含7個(gè)EXO70基因,海島棉EXO70基因在染色體上的分布情況與陸地棉一致。

    但這個(gè)趨勢越來越明顯。像教育、家裝、交通,醫(yī)療,都在發(fā)生很大的變化,都出現(xiàn)了有一定代表性的網(wǎng)絡(luò)協(xié)同或智能商業(yè)的企業(yè)。

    2.5 棉花EXO70家族基因組織表達(dá)模式分析及比較

    根據(jù)已發(fā)表的陸地棉和海島棉不同組織和不同發(fā)育階段的轉(zhuǎn)錄組數(shù)據(jù)(電子附表5和電子附表6),比較分析陸地棉和海島棉EXO70基因的表達(dá)模式(圖4)??傮w來看,不同EXO70家族分支呈現(xiàn)不同的表達(dá)規(guī)律,陸地棉和海島棉在同一分支中的表達(dá)趨勢又大致相似。在分支A中,陸地棉與海島棉存在4對(duì)同源基因在纖維發(fā)育時(shí)期較其他時(shí)期呈現(xiàn)優(yōu)勢表達(dá),分支F中的4對(duì)同源基因與分支A情況一致。分支B中的2對(duì)同源基因都呈現(xiàn)組成型表達(dá),在主要組織和纖維發(fā)育時(shí)期表達(dá)量都較高。陸地棉與海島棉在分支C與分支G中的成員在組織和胚珠纖維發(fā)育中表達(dá)量極低。由圖4可知,存在多數(shù)基因呈組成型表達(dá),陸地棉中共有18個(gè),海島棉中共有17個(gè)。陸地棉在纖維發(fā)育時(shí)期表達(dá)基因有34個(gè),海島棉中有36個(gè),多于陸地棉。

    圖2 陸地棉(A)和海島棉(B)EXO70基因家族系統(tǒng)進(jìn)化樹、蛋白基序、超家族和基因結(jié)構(gòu)

    圖3 陸地棉(A)與海島棉(B)EXO70家族種內(nèi)共線性分析

    海島棉和陸地棉同源基因表達(dá)模式大致相似,但在纖維發(fā)育時(shí)期中,海島棉整體表達(dá)量高于陸地棉。如,分支C中的海島棉在纖維發(fā)育時(shí)期20 DPA時(shí)表達(dá)量高于同源基因(圖4)。從纖維整體發(fā)育時(shí)期來看,陸地棉與海島棉同源基因的表達(dá)峰值時(shí)期也有所差異。分支A中陸地棉與同源基因都在纖維發(fā)育后期表達(dá)較高(圖4),前者集中表達(dá)在10 DPA,而海島棉在20 DPA表達(dá)達(dá)到高峰并高于陸地棉,表達(dá)高峰較陸地棉延后。為了進(jìn)一步驗(yàn)證轉(zhuǎn)錄組數(shù)據(jù)的準(zhǔn)確性,隨機(jī)選取4個(gè)基因設(shè)計(jì)引物(電子附表8)對(duì)其在纖維發(fā)育時(shí)期中的0、3、5和10 DPA進(jìn)行定量表達(dá)分析驗(yàn)證,發(fā)現(xiàn)表達(dá)趨勢均與轉(zhuǎn)錄組數(shù)據(jù)一致(圖4和圖5)。

    2.6 棉花EXO70家族與重要農(nóng)藝性狀的關(guān)聯(lián)分析

    通過對(duì)EXO70家族基因的不同單倍型的表型特征進(jìn)行差異分析,發(fā)現(xiàn)15個(gè)陸地棉和18個(gè)海島棉的EXO70基因的不同單倍型存在顯著性差異(<0.05)(電子附表9)。經(jīng)過分析對(duì)比可知,海島棉EXO70基因成員大部分與纖維長度和纖維強(qiáng)度這兩個(gè)纖維品質(zhì)性狀相關(guān),陸地棉EXO70基因成員更多是與產(chǎn)量性狀衣分相關(guān)(電子附表9)。如圖6-A可知,與纖維長度性狀呈顯著相關(guān)的海島棉基因個(gè)數(shù)有15個(gè),而陸地棉只有7個(gè),與纖維強(qiáng)度性狀顯著相關(guān)的海島棉基因有16個(gè),而陸地棉只有7個(gè)。在產(chǎn)量相關(guān)性狀——衣分方面,陸地棉與之顯著相關(guān)基因個(gè)數(shù)多于海島棉,分別為13和8個(gè)。分析比對(duì)同源基因可知,陸地棉與海島棉只有2對(duì)基因與纖維長度性狀相關(guān)聯(lián),有3對(duì)同源基因與纖維強(qiáng)度相關(guān)聯(lián),有2對(duì)同源基因與衣分相關(guān)聯(lián)。其中,只有分支B中的陸地棉和海島棉同源基因可找到與3個(gè)性狀同時(shí)關(guān)聯(lián)的位點(diǎn)。

    結(jié)合EXO70家族成員基因結(jié)構(gòu)分析與基因關(guān)聯(lián)分析可知,基因結(jié)構(gòu)存在差異的陸地棉與海島棉同源基因,在關(guān)聯(lián)性狀時(shí)也有差異。如,在分支H中的陸地棉與同源基因海島棉編碼序列存在較大差異(圖6-B),前者相對(duì)于后者存在提前終止情況,而海島棉在A11:115905472位點(diǎn)處可顯著關(guān)聯(lián)纖維長度、纖維強(qiáng)度性狀但陸地棉同源基因在對(duì)應(yīng)位點(diǎn)無關(guān)聯(lián)性狀。分支A中存在編碼序列差異的與(電子附表3),海島棉在D09:42925528與D09:42925541 2個(gè)位點(diǎn)與纖維長度、纖維強(qiáng)度和衣分3個(gè)性狀顯著關(guān)聯(lián)(電子附表8),而在陸地棉同源基因?qū)?yīng)位點(diǎn)處無性狀關(guān)聯(lián)。

    圖5 4個(gè)基因在纖維發(fā)育時(shí)期的定量表達(dá)分析

    結(jié)合組織表達(dá)模式與關(guān)聯(lián)分析結(jié)果可知,纖維發(fā)育時(shí)期的表達(dá)差異與基因關(guān)聯(lián)表型的差異相互印證。如,陸地棉在纖維發(fā)育時(shí)期基本不表達(dá),而海島棉同源基因有表達(dá),并且在20 DPA時(shí)達(dá)到表達(dá)高峰(圖4),群體數(shù)據(jù)中顯示部分海島棉的 EXO基因可關(guān)聯(lián)纖維長度、纖維強(qiáng)度等纖維品質(zhì)性狀,而陸地棉并無EXO70基因與此類纖維相關(guān)性狀關(guān)聯(lián),且只能關(guān)聯(lián)到衣分這種產(chǎn)量性狀(電子附表9)。

    2.7 棉花EXO70家族基因脅迫表達(dá)模式分析及比較

    如圖7所示,在相同算法條件下,將脅迫處理中是否受誘導(dǎo)基因區(qū)分,陸地棉和海島棉的同源基因在應(yīng)對(duì)低溫(4 ℃)、高溫(37 ℃)、鹽脅迫(NaCl)、干旱脅迫(PEG)逆境條件下響應(yīng)表達(dá)有所不同。4種非生物脅迫條件下,陸地棉受逆境脅迫誘導(dǎo)的EXO70基因有25個(gè),而海島棉只存在10個(gè)。整體數(shù)量上看,陸地棉要多于海島棉:在低溫4 ℃脅迫條件下,陸地棉受誘導(dǎo)基因有11個(gè),海島棉中只有1個(gè);在高溫37 ℃脅迫條件下,受誘導(dǎo)陸地棉基因個(gè)數(shù)為11,而海島棉為6個(gè);在鹽(NaCl)脅迫條件下,陸地棉有14個(gè)基因受誘導(dǎo),而海島棉中有2個(gè);在干旱(PEG)脅迫條件下,陸地棉有4個(gè)基因受誘導(dǎo),海島棉中也有4個(gè)。

    不同進(jìn)化分支中,陸地棉與海島棉在4種非生物脅迫條件下受誘導(dǎo)情況各自都具有一定的偏好性。陸地棉受誘導(dǎo)基因分布在分支C、D、E、F、G和H中,集中分布于E和H這2個(gè)分支中。而海島棉受誘導(dǎo)基因存在于A、C和H中,在分支H中受誘導(dǎo)基因個(gè)數(shù)最多,為5個(gè),占總海島棉EXO70整體受誘導(dǎo)基因個(gè)數(shù)的50%,其余分支不存在海島棉基因受誘導(dǎo)的情況。

    對(duì)于同一分支中,陸地棉與海島棉同源基因在不同脅迫條件下受誘導(dǎo)情況也不同。例如,在分支H中的陸地棉在4 ℃與NaCl條件下都會(huì)受到誘導(dǎo),而海島棉同源基因卻只在干旱條件(PEG)下受誘導(dǎo)。存在于分支C的陸地棉在高溫、低溫條件下都會(huì)受誘導(dǎo),而同源基因只在高溫中存在誘導(dǎo)情況。

    A:陸地棉(Gh)和海島棉(Gb)與纖維長度(FL)、纖維強(qiáng)度(FS)、衣分(LP)性狀呈顯著相關(guān)基因個(gè)數(shù);B:GH_D09G1462與同源基因GB_D09G1473的基因結(jié)構(gòu),和位于外顯子的非同義SNP(上部);GB_D09G1473的2個(gè)非同義SNP的不同單倍型在纖維長度、纖維強(qiáng)度和衣分的差異;****表示在P<0.05時(shí)單倍型間差異極顯著(t檢驗(yàn))

    3 討論

    3.1 棉花EXO70基因家族的鑒定和表達(dá)分析為后續(xù)深入研究提供借鑒

    陸地棉和海島棉是目前世界上最主要的2個(gè)棉花栽培種,且均已完成全基因組序列的測序和組裝工作[39-41],這為本研究進(jìn)行海島棉和陸地棉中EXO70基因家族的鑒定和分析提供極有利的先決條件。本文從陸地棉(TM-1)、海島棉(Hai7124)2個(gè)四倍體栽培棉種中都分別鑒定出54個(gè)EXO70基因,說明EXO70這類多拷貝基因在四倍體棉花的分別進(jìn)化階段在數(shù)量上較為保守。二倍體擬南芥EXO70基因家族成員數(shù)是23個(gè),本文鑒定出異源四倍體棉花EXO70家族成員數(shù)量是擬南芥的2.35倍,可能是因?yàn)槊藁ㄟM(jìn)化過程中發(fā)生了基因組的加倍和復(fù)制事件。近幾年,隨著測序技術(shù)的不斷更新和改進(jìn),基因組測序和組裝的水平得到了大幅提高。這樣的技術(shù)進(jìn)展導(dǎo)致不同的基因組版本中基因注釋等信息存在一些差異。未來,整合多個(gè)版本的基因組信息和分子試驗(yàn),可以更加有效和準(zhǔn)確地鑒定EXO70基因家族成員。酵母與大多數(shù)動(dòng)物的EXO70基因?yàn)閱慰截?,而大部分植物尤其是陸生植物中EXO70基因?yàn)槎嗫截?,這可能是因?yàn)橹参锞哂卸喾N多樣的細(xì)胞壁及復(fù)雜的液泡結(jié)構(gòu),需要多種胞泌復(fù)合體滿足其適應(yīng)環(huán)境,進(jìn)而演化出不同基因特有的調(diào)控機(jī)制[42-43],相關(guān)報(bào)道也證明不同的EXO70蛋白會(huì)參與細(xì)胞壁形成、細(xì)胞分泌和細(xì)胞免疫等不同的重要細(xì)胞學(xué)過程[44]。本研究鑒定的EXO70家族成員中既有組織特異型基因,也有少數(shù)組成型基因,這說明不同的棉花EXO70基因通過組織特異的表達(dá)模式參與不同組織的生長發(fā)育進(jìn)程。已有報(bào)道證明,在非生物脅迫下,囊泡運(yùn)輸相關(guān)蛋白可通過參與質(zhì)膜修復(fù)來增加植物對(duì)溫度等非生物脅迫的耐受性[28],而在本研究的另一個(gè)表達(dá)模式——逆境脅迫誘導(dǎo)表達(dá)模式研究中,不同分支的EXO70成員存在脅迫誘導(dǎo)多樣化、不同基因特異逆境誘導(dǎo)的情況,說明EXO70基因在棉花抗逆境脅迫過程中具有一定的功能。

    熱圖表示基因在應(yīng)對(duì)非生物脅迫過程中,是否顯著差異表達(dá)(紅色,差異表達(dá);白色,未差異表達(dá))

    3.2 EXO70家族基因在棉花進(jìn)化過程中具有保守性

    陸地棉與海島棉EXO70家族同源基因結(jié)構(gòu)有所差異,并且外顯子的個(gè)數(shù)也有所不同,而內(nèi)含子-外顯子的基因結(jié)構(gòu)通常被認(rèn)為是進(jìn)化的印記[45]。陸地棉與海島棉同源基因存在的基因結(jié)構(gòu)差異,可能會(huì)導(dǎo)致二者基因表達(dá)的差異從而引起基因功能的不同。從種內(nèi)共線性分析來看,陸地棉與海島棉中存在發(fā)生同義突變的基因?qū)?,說明這些序列分歧度太大,進(jìn)化距離很遠(yuǎn)。此外,陸地棉與海島棉都只有1個(gè)基因?qū)Φ腒a/Ks值大于1,其余都小于1,說明這些基因是受到正向選擇,EXO70家族基因在進(jìn)化過程中經(jīng)歷了較為強(qiáng)烈的凈化選擇[46],暗示大部分EXO70基因進(jìn)化保守,功能穩(wěn)定。

    3.3 EXO70家族可能參與棉花纖維發(fā)育過程

    棉花纖維的各個(gè)發(fā)育時(shí)期對(duì)纖維品質(zhì)與產(chǎn)量都至關(guān)重要。棉花纖維是從棉花胚珠表皮單個(gè)細(xì)胞突起進(jìn)而發(fā)育形成,與擬南芥表皮毛形成過程類似,并且當(dāng)今學(xué)者普遍認(rèn)為二者表皮毛發(fā)育相關(guān)基因的調(diào)控網(wǎng)絡(luò)類似[47]。在擬南芥表皮毛無胼胝質(zhì)的突變體中,功能缺失可導(dǎo)致表皮毛中胼胝質(zhì)合成酶PMR4分泌受阻,在毛狀體成熟階段細(xì)胞壁成分積累產(chǎn)生缺陷,而野生型毛狀體細(xì)胞壁可發(fā)生大量二次增厚[42, 48]。還有研究表明是在擬南芥氣管元件(tracheary element)發(fā)育過程中表達(dá)的基因,并在TE的囊泡運(yùn)輸中起作用從而調(diào)節(jié)細(xì)胞壁增厚[16]。本研究通過對(duì)比陸地棉與海島棉EXO70基因家族中同源基因在編碼序列、組織表達(dá)及重要農(nóng)藝性狀的基因關(guān)聯(lián)分析,可為后續(xù)棉花選育與改良方面結(jié)合陸地棉與海島棉優(yōu)勢基因,提供基礎(chǔ)。

    4 結(jié)論

    陸地棉和海島棉均存在54個(gè)EXO70家族基因,在進(jìn)化過程中基因家族數(shù)目相對(duì)保守,但基因結(jié)構(gòu)序列存在部分差異,同源基因的分子功能也發(fā)生了分化。海島棉存在更多EXO70基因與纖維長度、強(qiáng)度等品質(zhì)性狀顯著關(guān)聯(lián),可能與其優(yōu)異的纖維品質(zhì)特性相關(guān);陸地棉存在更多EXO70基因與衣分、籽指等產(chǎn)量性狀顯著關(guān)聯(lián),并且具有更強(qiáng)的非生物脅迫響應(yīng)能力,可能與其高產(chǎn)、廣適性等特性相關(guān)。

    [1] LI S P, VAN OS G M A, REN S C, YU D L, KETELAAR T, EMONS A M C, LIU C M. Expression and functional analyses ofgenes inimplicate their roles in regulating cell type-specific exocytosis. Plant Physiology, 2010, 154(4): 1819-1830.

    [2] WU L G, HAMID E, SHIN W, CHIANG H C. Exocytosis and endocytosis: modes, functions, and coupling mechanisms. Annual Review of Physiology, 2014, 76: 301-331.

    [3] ZHANG L, XING J J, LIN J X. At the intersection of exocytosis and endocytosis in plants. The New Phytologist, 2019, 224(4): 1479-1489.

    [4] BONIFACINO J S, GLICK B S. The mechanisms of vesicle budding and fusion. Cell, 2004, 116(2): 153-166.

    [5] PLESKOT R, CWIKLIK L, JUNGWIRTH P, ?áRSKY V, POTOCKY M. Membrane targeting of the yeast exocyst complex. Biochimica et Biophysica Acta, 2015, 1848(7): 1481-1489.

    [6] YU I M, HUGHSON F M. Tethering factors as organizers of intracellular vesicular traffic. Annual Review of Cell and Developmental Biology, 2010, 26: 137-156.

    [7] TERBUSH D R, MAURICE T, ROTH D, NOVICK P. The Exocyst is a multiprotein complex required for exocytosis in. The EMBO Journal, 1996, 15(23): 6483-6494.

    [8] TERBUSH D R, NOVICK P. Sec6, Sec8, and Sec15 are components of a multisubunit complex which localizes to small bud tips in. The Journal of Cell Biology, 1995, 130(2): 299-312.

    [9] HE B, GUO W. The exocyst complex in polarized exocytosis. Current Opinion in Cell Biology, 2009, 21(4): 537-542.

    [10] SYNEK L, SCHLAGER N, ELIáS M, QUENTIN M, HAUSER M-T, ZáRSKy V. AtEXO70A1, a member of a family of putative exocyst subunits specifically expanded in land plants, is important for polar growth and plant development. The Plant Journal, 2006, 48(1): 54-72.

    [11] FENDRYCH M, SYNEK L, PE?ENKOVá T, TOUPALOVá H, COLE R, DRDOVá E, NEBESá?OVá J, ?EDINOVá M, HáLA M, FOWLER J E, ?áRSKY V. Theexocyst complex is involved in cytokinesis and cell plate maturation. The Plant Cell, 2010, 22(9): 3053-3065.

    [12] WANG W, LIU N, GAO C Y, CAI H R, ROMEIS T, TANG D Z. Theexocyst subunits EXO70B1 and EXO70B2 regulate FLS2 homeostasis at the plasma membrane. The New Phytologist, 2020, 227(2): 529-544.

    [13] KULICH I, PE?ENKOVá T, SEKERE? J, SMETANA O, FENDRYCH M, FOISSNER I, H?FTBERGER M, ZáRSKy V.exocyst sub complex containing subunit EXO70B1 is involved in autophagy-related transport to the vacuole. Traffic (Copenhagen, Denmark), 2013, 14(11): 1155-1165.

    [14] ZáRSKy V, KULICH I, FENDRYCH M, PE?ENKOVá T. Exocyst complexes multiple functions in plant cells secretory pathways. Current Opinion in Plant Biology, 2013, 16(6): 726-733.

    [15] LIU N, HAKE K, WANG W, ZHAO T, ROMEIS T, TANG D Z. Calcium-dependent protein kinase5 associates with the truncated NLR protein TIR-NBS2 to contribute to-mediated immunity. The Plant Cell, 2017, 29(4): 746-759.

    [16] LI S P, CHEN M, YU D L, REN S C, SUN S F, LIU L D, KETELAAR T, EMONS A M C, LIU C M. EXO70A1-mediated vesicle trafficking is critical for tracheary element development in. The Plant Cell, 2013, 25(5): 1774-1786.

    [17] KALMBACH L, HéMATY K, DE BELLIS D, BARBERON M, FUJITA S, URSACHE R, DARASPE J, GELDNER N. Transient cell-specific EXO70A1 activity in the CASP domain and Casparian strip localization. Nature Plants, 2017, 3: 17058.

    [18] SYNEK L, VUKA?INOVI? N, KULICH I, HáLA M, ALDORFOVá K, FENDRYCH M, ?áRSKY V. EXO70C2 is a key regulatory factor for optimal tip growth of pollen. Plant Physiology, 2017, 174(1): 223-240.

    [19] TU B, HU L, CHEN W L, LI T, HU B H, ZHENG L, LV Z, YOU S J, WANG Y P, MA B T, CHEN X W, QIN P, LI S G. Disruption ofcauses irregular vascular bundles and perturbs mineral nutrient assimilation in rice. Scientific Reports, 2015, 5: 18609.

    [20] FUJISAKI K, ABE Y, ITO A, SAITOH H, YOSHIDA K, KANZAKI H, KANZAKI E, UTSUSHI H, YAMASHITA T, KAMOUN S, TERAUCHI R. Rice Exo70 interacts with a fungal effector, AVR-Pii, and is required for AVR-Pii-triggered immunity. The Plant Journal, 2015, 83(5): 875-887.

    [21] DU Y, OVERDIJK E J R, BERG J A, GOVERS F, BOUWMEESTER K. Solanaceous exocyst subunits are involved in immunity to diverse plant pathogens. Journal of Experimental Botany, 2018, 69(3): 655-666.

    [22] WANG H, MA Z H, MAO J, CHEN B H. Genome-wide identification and expression analysis of the EXO70 gene family in grape (l.). PeerJ, 2021, 9: e11176.

    [23] 丁顏朋, 葛曉陽, 王鵬, 吳潔, 王省芬, 李付廣. 棉花鹽脅迫應(yīng)答基因功能分析. 棉花學(xué)報(bào), 2018, 30(6): 423-434.

    DING Y P, GE X Y, WANG P, WU J, WANG X F, LI F G. Functional analysis of a salt stress response gene,, in upland cotton (L.). Cotton Science, 2018, 30(6): 423-434. (in Chinese)

    [24] ZHU Y Q, QIU L, LIU L L, LUO L, HAN X P, ZHAI Y H, WANG W J, REN M Z, XING Y D. Identification and comprehensive structural and functional analyses of thegene family in cotton. Genes, 2021, 12(10): 1594.

    [25] DAI F, CHEN J D, ZHANG Z Q, LIU F J, LI J, ZHAO T, HU Y, ZHANG T Z, FANG L. COTTONOMICS: A comprehensive cotton multi-omics database. Database, 2022, 2022: baac080.

    [26] EL-GEBALI S, MISTRY J, BATEMAN A, EDDY S R, LUCIANI A, POTTER S C, QURESHI M, RICHARDSON L J, SALAZAR G A, SMART A, SONNHAMMER E L, HIRSH L, PALADIN L, PIOVESAN D, TOSATTO S C, FINN R D. The Pfam protein families database in 2019. Nucleic Acids Research, 2019, 47(D1): D427-D432.

    [27] CAMACHO C, COULOURIS G, AVAGYAN V, MA N, PAPADOPOULOSJ, BEALER K, MADDEN T L. BLAST+: architecture and applications. BMC Bioinformatics, 2009, 10: 421.

    [28] BAILEY T L, BODEN M, BUSKE F A, FRITH M, GRANT C E, CLEMENTI L, REN J Y, LI W W, NOBLE W S. MEME Suite: tools for motif discovery and searching. Nucleic Acids Research, 2009, 37(supp1_2): W202- W208.

    [29] CHEN C J, CHEN H, ZHANG Y, THOMAS H R, FRANK M H, HE Y H, XIA R. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Molecular Plant, 2020, 13(8): 1194-1202.

    [30] WANG Y P, TANG H B, DEBARRY J D, TAN X, LI J P, WANG X Y, LEE T H, JIN H Z, MARLER B, GUO H, KISSINGER J C, PATERSON A H.: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Research, 2012, 40(7): e49.

    [31] KRZYWINSKI M, SCHEIN J, BIROL I, CONNORS J, GASCOYNE R, HORSMAN D, JONES S J, MARRA M A. Circos: An information aesthetic for comparative genomics. Genome Research, 2009, 19(9): 1639-1645.

    [32] CHEN S F, ZHOU Y Q, CHEN Y R, GU J. Fastp: an ultra-fast all-in- one FASTQ preprocessor. Bioinformatics, 2018, 34(17): i884-i890.

    [33] KIM D, LANGMEAD B, SALZBERG S L. HISAT: a fast spliced aligner with low memory requirements. Nature Methods, 2015, 12(4): 357-360.

    [34] LIAO Y, SMYTH G K, SHI W. FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics, 2014, 30(7): 923-930.

    [35] HU Y, CHEN J D, FANG L, ZHANG Z Y, MA W, NIU Y C, JU L Z, DENG J Q, ZHAO T, LIAN J M, BARUCH K, FANG D, LIU X, RUAN Y L, RAHMAN M U, HAN J L, WANG K, WANG Q, WU H T, MEI G F, ZANG Y H, HAN Z G, XU C Y, SHEN W J, YANG D F, SI Z F, DAI F, ZOU L F, HUANG F, BAI Y L, ZHANG Y G, BRODT A, BEN-HAMO H, ZHU X F, ZHOU B L, GUAN X Y, ZHU S J, CHEN X Y, ZHANG T Z.andgenomes provide insights into the origin and evolution of allotetraploid cotton. Nature Genetics, 2019, 51(4): 739-748.

    [36] FANG L, GONG H, HU Y, LIU C X, ZHOU B L, HUANG T, WANG Y K, CHEN S Q, FANG D D, DU X M, CHEN H, CHEN J D, WANG S, WANG Q, WAN Q, LIU B L, PAN M Q, CHANG L J, WU H T, MEI G F, XIANG D, LI X H, CAI C P, ZHU X F, CHEN Z J, HAN B, CHEN X Y, GUO W Z, ZHANG T Z, HUANG X H. Genomic insights into divergence and dual domestication of cultivated allotetraploid cottons. Genome Biology, 2017, 18(1): 33.

    [37] FANG L, ZHAO T, HU Y, SI Z F, ZHU X F, HAN Z G, LIU G Z, WANG S, JU L Z, GUO M L, MEI H, WANG L Y, QI B W, WANG H, GUAN X Y, ZHANG T Z. Divergent improvement of two cultivated allotetraploid cotton species. Plant Biotechnology Journal, 2021, 19(7): 1325-1336.

    [38] WANG K, LI M Y, HAKONARSON H. Annovar: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Research, 2010, 38(16): e164.

    [39] LI F G, FAN G Y, LU C R, XIAO G H, ZOU C S, KOHEL R J, MA Z Y, SHANG H H, MA X F, WU J Y, LIANG X M, HUANG G, PERCY R G, LIU K, YANG W H, CHEN W B, DU X M, SHI C C, YUAN Y L, YE W W, LIU X, ZHANG X Y, LIU W Q, WEI H L, WEI S J, HUANG G D, ZHANG X L, ZHU S J, ZHANG H, SUN F M, WANG X F, LIANG J, WANG J H, HE Q, HUANG L H, WANG J, CUI J J, SONG G L, WANG K B, XU X, YU J Z, ZHU Y X, YU S X. Genome sequence of cultivated Upland cotton (TM-1) provides insights into genome evolution. Nature Biotechnology, 2015, 33(5): 524-530.

    [40] ZHU T, LIANG C Z, MENG Z G, SUN G Q, MENG Z, GUO S D, ZHANG R. CottonFGD: an integrated functional genomics database for cotton. BMC Plant Biology, 2017, 17(1): 101.

    [41] WANG M J, TU L L, YUAN D J, ZHU D, SHEN C, LI J Y, LIU F Y, PEI L L, WANG P C, ZHAO G N, YE Z X, HUANG H, YAN F L, MA Y Z, ZHANG L, LIU M, YOU J Q, YANG Y C, LIU Z P, HUANG F, LI B Q, QIU P, ZHANG Q H, ZHU L F, JIN S X, YANG X Y, MIN L, LI G L, CHEN L L, ZHENG H K, LINDSEY K, LIN Z X, UDALL J A, ZHANG X L. Reference genome sequences of two cultivated allotetraploid cottons,and. Nature Genetics, 2019, 51(2): 224-229.

    [42] KULICH I, VOJTíKOVá Z, SABOL P, ORTMANNOVá J, NEDěLA V, TIHLA?íKOVá E, ?áRSKY V. Exocyst subunit EXO70H4 has a specific role in callose synthase secretion and silica accumulation. Plant Physiology, 2018, 176(3): 2040-2051.

    [43] YAMAZAKI T, KAWAMURA Y, MINAMI A, UEMURA M. Calcium-dependent freezing tolerance ininvolves membrane resealing via synaptotagmin SYT1. The Plant Cell, 2009, 20(12): 3389-3404.

    [44] SEKERE? J, PEJCHAR P, ?ANTR??EK J, VUKA?INOVI? N, ?áRSKy V, POTOCKy M. Analysis of exocyst subunit EXO70 family reveals distinct membrane polar domains in tobacco pollen tubes. Plant Physiology, 2017, 173(3): 1659-1675.

    [45] YU H H, LI M, SANDHU J, SUN G C, SCHNABLE J C, WALIA H, XIE W B, YU B, MOWER J P, ZHANG C. Pervasive misannotation of microexons that are evolutionarily conserved and crucial for gene function in plants. Nature Communications, 2022, 13: 820.

    [46] HURST L D. Theratio: Diagnosing the form of sequence evolution. Trends in Genetics, 2002, 18(9): 486.

    [47] GUAN X Y, YU N, SHANGGUAN X X, WANG S, LU S, WANG L J, CHEN X Y.trichome research sheds light on cotton fiber development mechanisms. Chinese Science Bulletin, 2007, 52(13): 1734-1741.

    [48] KULICH I, VOJTíKOVá Z, GLANC M, ORTMANNOVá J, RASMANN S, ?áRSKy V. Cell wall maturation oftrichomes is dependent on exocyst subunit EXO70H4 and involves callose deposition. Plant Physiology, 2015, 168(1): 120-131.

    Genome-wide identification and interspecific comparative analysis of the EXO70 gene family in cotton

    DONG YanYu, XU BiYu, DONG ZeYu, WANG LuYao, CHEN JinWen, FANG Lei

    College of Agriculture & Biotechnology, Zhejiang University, Hangzhou 310058

    【Objective】The two main types of allotetraploid cotton that are currently cultivated are upland cotton, known for its high yield and good adaptability, and island cotton, which boasts excellent fiber quality but lower yield. EXO70, a vital subunit of the exocyst complex, plays a significant role in plant growth, development, and stress response. By identifying and analyzing members of the EXO70 gene family in upland and island cotton at the whole-genome level, and studying their functions in fiber development and environmental adaptation, we can shed light on the molecular basis for the differences in traits between these two varieties. 【Method】The reference sequences of EXO70 protein inwere obtained from the TAIR database. HMMER, ExPASy, MEME, TBtools, and other analysis tools were used to identify and analyze the members of EXO70 gene family in the genomes of upland cotton TM-1 and island cotton Hai7124. The similarities and differences in gene expression patterns, correlations with crucial agronomic traits, and stress responses of this family were systematically compared. 【Result】Through genome-level analysis of upland cotton and sea-island cotton, 54 EXO70 family members were identified in both upland and sea-island cotton, which could be divided into eight subgroups based on phylogenetic analysis. Orthologous genes between upland cotton and sea-island cotton can be paired one-to-one and are distributed across the 20 chromosomes of both species. The majority of the members have single exons, while 12 pairs of homologous genes displayed significant differences in the reading frame sequences. Most orthologous genes inanddisplay similar expression patterns, but differences in expression levels are observed during the same fiber development stage, such asand its orthologous geneand. Under different stresses, upland cotton showed a significantly higher number of induced genes compared to sea-island cotton. 【Conclusion】The sequence structure and gene expression patterns of the EXO70 family were found to be relatively conserved in both upland and island cotton during the formation, differentiation, and domestication of tetraploid cotton. However, in terms of EXO70 family members, island cotton had more genes related to fiber quality traits, while upland cotton had more genes related to yield traits and exhibited greater sensitivity to environmental stress.

    ;; EXO70 gene family; expression pattern; trait association analysis; abiotic stress response

    10.3864/j.issn.0578-1752.2023.23.005

    2023-02-14;

    2023-05-04

    國家自然科學(xué)基金面上基金(32172008)、中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金(226-2022-00153)

    董琰鈺,E-mail:dongyy@zju.edu.cn。通信作者方磊,E-mail:fangl@zju.edu.cn

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

    猜你喜歡
    海島棉同源分支
    藥食同源
    ——紫 蘇
    兩岸年味連根同源
    以同源詞看《詩經(jīng)》的訓(xùn)釋三則
    巧分支與枝
    一類擬齊次多項(xiàng)式中心的極限環(huán)分支
    轉(zhuǎn)Bt基因海島棉分子檢測及抗蟲性鑒定
    棉花陸海雜交親本及子代光合葉綠素?zé)晒鈪?shù)分析
    虔誠書畫乃同源
    新疆海島棉生產(chǎn)現(xiàn)狀與發(fā)展建議
    海島棉染色體片段代換系BC4F4產(chǎn)量及纖維品質(zhì)主成分分析
    欧美午夜高清在线| 在线天堂最新版资源| 尤物成人国产欧美一区二区三区| 午夜福利在线观看吧| 国产精品免费一区二区三区在线| 亚洲中文字幕一区二区三区有码在线看| 日本三级黄在线观看| 国产高清有码在线观看视频| 国产精品av视频在线免费观看| 亚洲国产日韩欧美精品在线观看| 欧美黄色淫秽网站| 免费观看的影片在线观看| 亚洲无线在线观看| 乱人视频在线观看| 波多野结衣高清无吗| 18+在线观看网站| 美女高潮的动态| 国产精品亚洲一级av第二区| 亚洲av电影不卡..在线观看| 性色av乱码一区二区三区2| 高清日韩中文字幕在线| 日韩成人在线观看一区二区三区| 国产私拍福利视频在线观看| 亚洲av日韩精品久久久久久密| 美女黄网站色视频| 国产乱人伦免费视频| 久久热精品热| 亚洲第一区二区三区不卡| 中文字幕高清在线视频| 久久久国产成人精品二区| 最近中文字幕高清免费大全6 | 国产三级中文精品| 99热精品在线国产| 亚洲精品乱码久久久v下载方式| 免费黄网站久久成人精品 | 69人妻影院| 日韩欧美精品免费久久 | 国产人妻一区二区三区在| 国内揄拍国产精品人妻在线| 欧美三级亚洲精品| 国产午夜福利久久久久久| 亚洲无线观看免费| 性欧美人与动物交配| 亚洲av第一区精品v没综合| 国产一区二区亚洲精品在线观看| 久久精品影院6| 亚洲综合色惰| 色哟哟哟哟哟哟| 久久久久久久亚洲中文字幕 | 看十八女毛片水多多多| 人妻制服诱惑在线中文字幕| 啦啦啦韩国在线观看视频| 亚洲人成网站在线播| 一级黄片播放器| 男人的好看免费观看在线视频| 久久99热这里只有精品18| 又粗又爽又猛毛片免费看| 每晚都被弄得嗷嗷叫到高潮| 在线a可以看的网站| 熟女电影av网| 午夜两性在线视频| 在线十欧美十亚洲十日本专区| 99久久精品一区二区三区| 色尼玛亚洲综合影院| 久久草成人影院| 国产69精品久久久久777片| 日本黄色片子视频| 简卡轻食公司| 国模一区二区三区四区视频| 桃红色精品国产亚洲av| 久久久久国内视频| 国产中年淑女户外野战色| 精品无人区乱码1区二区| 久久精品综合一区二区三区| 无遮挡黄片免费观看| 桃色一区二区三区在线观看| 亚洲最大成人av| 国产成人福利小说| 美女高潮喷水抽搐中文字幕| 男女床上黄色一级片免费看| 欧美黑人欧美精品刺激| or卡值多少钱| 99久久无色码亚洲精品果冻| 夜夜夜夜夜久久久久| 校园春色视频在线观看| 麻豆成人午夜福利视频| 美女cb高潮喷水在线观看| 成人无遮挡网站| 色精品久久人妻99蜜桃| 色精品久久人妻99蜜桃| 国产精品嫩草影院av在线观看 | 成人特级黄色片久久久久久久| 欧美+日韩+精品| 在线十欧美十亚洲十日本专区| 精品福利观看| 国产精品一区二区性色av| 婷婷精品国产亚洲av在线| 国产精品国产高清国产av| 国产亚洲精品av在线| 麻豆一二三区av精品| 国产在线男女| 亚洲人成网站在线播放欧美日韩| 男女床上黄色一级片免费看| 男人和女人高潮做爰伦理| 国产日本99.免费观看| 一本一本综合久久| 美女cb高潮喷水在线观看| 一个人免费在线观看电影| 18禁在线播放成人免费| 精品免费久久久久久久清纯| 久久性视频一级片| 嫁个100分男人电影在线观看| 免费搜索国产男女视频| 欧美在线一区亚洲| 一边摸一边抽搐一进一小说| 嫁个100分男人电影在线观看| 久久久久精品国产欧美久久久| 成人特级黄色片久久久久久久| 男人舔女人下体高潮全视频| 精品国产三级普通话版| 无遮挡黄片免费观看| 国产高潮美女av| av在线蜜桃| 午夜福利18| 亚洲av免费在线观看| 久久久久久久午夜电影| 久久婷婷人人爽人人干人人爱| 亚洲国产精品999在线| 99热只有精品国产| 成人特级av手机在线观看| 老鸭窝网址在线观看| 九色国产91popny在线| 夜夜看夜夜爽夜夜摸| 精品人妻1区二区| av中文乱码字幕在线| 午夜福利在线在线| 亚洲精品一区av在线观看| 日韩成人在线观看一区二区三区| www.色视频.com| 成人特级av手机在线观看| 一个人免费在线观看电影| 国内毛片毛片毛片毛片毛片| 51午夜福利影视在线观看| 亚洲精品久久国产高清桃花| 狂野欧美白嫩少妇大欣赏| 国产v大片淫在线免费观看| 婷婷亚洲欧美| 久久久久久国产a免费观看| 亚洲精品久久国产高清桃花| x7x7x7水蜜桃| 久久99热这里只有精品18| 国产不卡一卡二| 国产一区二区三区在线臀色熟女| 精品久久久久久久末码| 久久九九热精品免费| 久99久视频精品免费| 免费在线观看亚洲国产| 天堂影院成人在线观看| 免费av观看视频| 一个人免费在线观看的高清视频| 波多野结衣巨乳人妻| 中文资源天堂在线| 永久网站在线| 午夜免费成人在线视频| 深夜a级毛片| 深夜a级毛片| 宅男免费午夜| 我要搜黄色片| 久久草成人影院| 欧美精品啪啪一区二区三区| 国产一区二区三区视频了| 两个人的视频大全免费| 国产av在哪里看| 18禁黄网站禁片免费观看直播| 天堂动漫精品| 久久久久九九精品影院| 男女床上黄色一级片免费看| 成人三级黄色视频| 成人午夜高清在线视频| 午夜亚洲福利在线播放| 国产单亲对白刺激| 三级男女做爰猛烈吃奶摸视频| 久久久久亚洲av毛片大全| 精品久久久久久久久av| 欧美3d第一页| 网址你懂的国产日韩在线| 黄片小视频在线播放| 日韩 亚洲 欧美在线| 国产亚洲精品av在线| 婷婷色综合大香蕉| 麻豆国产av国片精品| 亚洲中文日韩欧美视频| 在线观看一区二区三区| 男人和女人高潮做爰伦理| 看黄色毛片网站| 丁香欧美五月| 欧美乱妇无乱码| 欧美潮喷喷水| 久久精品国产自在天天线| a在线观看视频网站| av专区在线播放| 日韩欧美国产一区二区入口| 久久精品影院6| 免费观看的影片在线观看| 国产成人啪精品午夜网站| 午夜福利在线观看免费完整高清在 | 日韩欧美精品v在线| 亚洲精华国产精华精| 淫妇啪啪啪对白视频| 国产av在哪里看| 一区二区三区四区激情视频 | 嫩草影院精品99| 一级毛片久久久久久久久女| 色视频www国产| 午夜免费男女啪啪视频观看 | 美女cb高潮喷水在线观看| 日本 欧美在线| 免费看日本二区| 中文在线观看免费www的网站| 色5月婷婷丁香| 色哟哟哟哟哟哟| 久久久精品欧美日韩精品| 精品久久久久久久久久久久久| 国产乱人视频| 91在线精品国自产拍蜜月| 亚洲成a人片在线一区二区| 最近最新中文字幕大全电影3| 蜜桃久久精品国产亚洲av| 成人av一区二区三区在线看| 久久天躁狠狠躁夜夜2o2o| 麻豆一二三区av精品| avwww免费| 国产成人av教育| 国产乱人伦免费视频| АⅤ资源中文在线天堂| 亚洲第一电影网av| 人人妻人人澡欧美一区二区| 青草久久国产| 国语自产精品视频在线第100页| 亚洲国产欧洲综合997久久,| 国产久久久一区二区三区| 日韩欧美 国产精品| 一个人免费在线观看电影| 国产在视频线在精品| 日日啪夜夜撸| 中文乱码字字幕精品一区二区三区| 国产成人精品福利久久| 久久精品综合一区二区三区| 久久精品国产亚洲av天美| 国产乱人视频| 国产亚洲av嫩草精品影院| 内射极品少妇av片p| 国产精品久久久久久精品电影| 亚洲精品成人av观看孕妇| 国产欧美亚洲国产| a级毛片免费高清观看在线播放| 亚洲图色成人| 黄片wwwwww| 精品人妻熟女av久视频| 大话2 男鬼变身卡| av国产久精品久网站免费入址| 欧美高清成人免费视频www| av免费在线看不卡| 亚洲人成网站在线播| 男女啪啪激烈高潮av片| 日本欧美国产在线视频| 国产高清不卡午夜福利| 最近中文字幕高清免费大全6| 男人添女人高潮全过程视频| 亚洲av.av天堂| 干丝袜人妻中文字幕| 国产成人精品久久久久久| 国产高清有码在线观看视频| 欧美老熟妇乱子伦牲交| 综合色av麻豆| 亚洲丝袜综合中文字幕| 免费不卡的大黄色大毛片视频在线观看| av在线天堂中文字幕| 亚洲av成人精品一区久久| 精华霜和精华液先用哪个| 日本一二三区视频观看| 亚洲色图av天堂| 噜噜噜噜噜久久久久久91| 国产精品麻豆人妻色哟哟久久| av国产免费在线观看| 国产 精品1| 久久久久久久久久人人人人人人| 91在线精品国自产拍蜜月| 人人妻人人看人人澡| 亚洲一级一片aⅴ在线观看| 国产欧美日韩一区二区三区在线 | 三级国产精品片| 欧美xxⅹ黑人| 国语对白做爰xxxⅹ性视频网站| 内射极品少妇av片p| 男女国产视频网站| 国产日韩欧美在线精品| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品国产av蜜桃| 神马国产精品三级电影在线观看| 97精品久久久久久久久久精品| 你懂的网址亚洲精品在线观看| 久久97久久精品| 看黄色毛片网站| 啦啦啦在线观看免费高清www| 内地一区二区视频在线| 精华霜和精华液先用哪个| 韩国av在线不卡| 麻豆久久精品国产亚洲av| 精品久久久久久久末码| 白带黄色成豆腐渣| 国产av国产精品国产| 97超碰精品成人国产| 欧美精品国产亚洲| 亚洲人成网站高清观看| 黑人高潮一二区| 又爽又黄a免费视频| 五月伊人婷婷丁香| 又爽又黄无遮挡网站| 国产乱人视频| 少妇人妻精品综合一区二区| 久久热精品热| 亚洲欧美中文字幕日韩二区| 国产免费一级a男人的天堂| 一区二区三区免费毛片| h日本视频在线播放| 一区二区三区四区激情视频| 久久久久国产网址| 69人妻影院| 亚洲欧美成人精品一区二区| 亚洲精品国产色婷婷电影| 夜夜爽夜夜爽视频| 一本色道久久久久久精品综合| 成年av动漫网址| 深夜a级毛片| 久久精品夜色国产| 日韩一本色道免费dvd| 一本久久精品| 熟女电影av网| 王馨瑶露胸无遮挡在线观看| 五月开心婷婷网| 超碰av人人做人人爽久久| 国产在线男女| 激情五月婷婷亚洲| 熟女av电影| 九色成人免费人妻av| 免费黄频网站在线观看国产| 最近最新中文字幕免费大全7| 欧美高清性xxxxhd video| 国产 一区精品| 干丝袜人妻中文字幕| 自拍偷自拍亚洲精品老妇| 99热全是精品| 汤姆久久久久久久影院中文字幕| 国产亚洲一区二区精品| av播播在线观看一区| 青春草国产在线视频| 天天躁夜夜躁狠狠久久av| 午夜爱爱视频在线播放| 男人舔奶头视频| 黄色视频在线播放观看不卡| 神马国产精品三级电影在线观看| 成年免费大片在线观看| 中文欧美无线码| 少妇熟女欧美另类| 国产精品国产三级专区第一集| 亚洲欧美一区二区三区国产| 日韩av在线免费看完整版不卡| 下体分泌物呈黄色| 精品人妻偷拍中文字幕| 久久精品国产亚洲网站| 少妇人妻一区二区三区视频| 极品教师在线视频| 日韩伦理黄色片| 嘟嘟电影网在线观看| 日韩欧美精品v在线| 一级av片app| 亚洲精品日韩av片在线观看| 蜜桃久久精品国产亚洲av| 亚洲欧美成人综合另类久久久| 国产一级毛片在线| 免费观看在线日韩| 国产精品精品国产色婷婷| 在线亚洲精品国产二区图片欧美 | 日韩强制内射视频| av在线老鸭窝| 午夜福利高清视频| av在线天堂中文字幕| 亚洲精品中文字幕在线视频 | 五月天丁香电影| 永久免费av网站大全| 五月天丁香电影| 亚洲av欧美aⅴ国产| 下体分泌物呈黄色| 2022亚洲国产成人精品| 啦啦啦啦在线视频资源| 国产精品无大码| 特大巨黑吊av在线直播| 少妇人妻 视频| 亚洲精品乱码久久久v下载方式| 日韩制服骚丝袜av| 最近最新中文字幕大全电影3| 国产欧美另类精品又又久久亚洲欧美| 国产欧美日韩精品一区二区| 国产精品爽爽va在线观看网站| 国产又色又爽无遮挡免| 久久人人爽av亚洲精品天堂 | 美女脱内裤让男人舔精品视频| 亚洲欧美清纯卡通| 在线看a的网站| 国产老妇伦熟女老妇高清| 亚洲欧美精品自产自拍| 精品人妻熟女av久视频| 日产精品乱码卡一卡2卡三| 国产精品一二三区在线看| av播播在线观看一区| 欧美日韩亚洲高清精品| 久久99热这里只有精品18| 日日啪夜夜爽| 亚洲精品乱久久久久久| 可以在线观看毛片的网站| 欧美国产精品一级二级三级 | 热99国产精品久久久久久7| 国产精品麻豆人妻色哟哟久久| 热99国产精品久久久久久7| 国产又色又爽无遮挡免| 深夜a级毛片| 成人二区视频| 亚洲精品乱码久久久久久按摩| 91久久精品国产一区二区成人| 18+在线观看网站| 国产成人精品婷婷| 伊人久久国产一区二区| eeuss影院久久| 91久久精品国产一区二区成人| 亚洲在久久综合| 久久97久久精品| 日日摸夜夜添夜夜添av毛片| 久久久国产一区二区| 国产在线一区二区三区精| 五月伊人婷婷丁香| 99久国产av精品国产电影| 黄色一级大片看看| 精品久久国产蜜桃| 一级爰片在线观看| 亚洲怡红院男人天堂| 嘟嘟电影网在线观看| 2021天堂中文幕一二区在线观| 丝袜脚勾引网站| 美女高潮的动态| 五月开心婷婷网| 深爱激情五月婷婷| a级一级毛片免费在线观看| 777米奇影视久久| 观看免费一级毛片| 亚洲av成人精品一二三区| 久久久久久久大尺度免费视频| 午夜福利在线在线| 水蜜桃什么品种好| 永久免费av网站大全| 最近手机中文字幕大全| 网址你懂的国产日韩在线| 久久97久久精品| 久久精品国产亚洲av涩爱| 观看免费一级毛片| 少妇的逼水好多| 高清av免费在线| 天堂网av新在线| 各种免费的搞黄视频| 欧美高清性xxxxhd video| 精品久久久久久久久亚洲| 一级av片app| 自拍偷自拍亚洲精品老妇| 99热6这里只有精品| 亚洲电影在线观看av| 性插视频无遮挡在线免费观看| 日本熟妇午夜| 欧美性感艳星| av在线亚洲专区| 深爱激情五月婷婷| 欧美变态另类bdsm刘玥| 成人国产麻豆网| 亚洲国产成人一精品久久久| 亚洲精品日韩av片在线观看| 成人鲁丝片一二三区免费| 国产极品天堂在线| 亚洲成人久久爱视频| 制服丝袜香蕉在线| 国语对白做爰xxxⅹ性视频网站| 国产精品一区二区在线观看99| 午夜福利视频精品| 白带黄色成豆腐渣| 亚洲aⅴ乱码一区二区在线播放| 特级一级黄色大片| 午夜亚洲福利在线播放| 日韩在线高清观看一区二区三区| 听说在线观看完整版免费高清| 另类亚洲欧美激情| 一级a做视频免费观看| 精品一区二区三卡| 国产亚洲91精品色在线| 久久久a久久爽久久v久久| 搡女人真爽免费视频火全软件| 国产精品一二三区在线看| 亚洲怡红院男人天堂| 成年女人在线观看亚洲视频 | 王馨瑶露胸无遮挡在线观看| 亚洲av日韩在线播放| 日韩 亚洲 欧美在线| 男女啪啪激烈高潮av片| 久久99热6这里只有精品| 国产成人91sexporn| 看黄色毛片网站| 毛片一级片免费看久久久久| 精品人妻偷拍中文字幕| 亚洲美女搞黄在线观看| 久久久久久久午夜电影| 少妇高潮的动态图| 男人爽女人下面视频在线观看| 免费观看a级毛片全部| 久久久a久久爽久久v久久| 成人漫画全彩无遮挡| 欧美最新免费一区二区三区| 亚洲欧美一区二区三区国产| 亚洲美女视频黄频| 看黄色毛片网站| 亚洲婷婷狠狠爱综合网| 亚洲av国产av综合av卡| 久久精品国产亚洲av天美| 久久久精品94久久精品| 久久久久久久久久久免费av| 最近2019中文字幕mv第一页| 亚洲国产av新网站| 免费av观看视频| 一级毛片我不卡| 午夜日本视频在线| 久久精品国产a三级三级三级| 丰满乱子伦码专区| 欧美日韩综合久久久久久| 久久精品国产亚洲网站| 日韩大片免费观看网站| 亚洲无线观看免费| 亚洲四区av| 亚洲自拍偷在线| 国产成人精品一,二区| 日韩成人av中文字幕在线观看| 制服丝袜香蕉在线| 久久久久久久大尺度免费视频| 国产av不卡久久| 搞女人的毛片| 亚洲综合精品二区| 汤姆久久久久久久影院中文字幕| 中文在线观看免费www的网站| 亚洲最大成人av| 国产高清不卡午夜福利| 韩国高清视频一区二区三区| 五月玫瑰六月丁香| 不卡视频在线观看欧美| 好男人在线观看高清免费视频| 日本色播在线视频| 成人综合一区亚洲| 久久6这里有精品| 欧美精品一区二区大全| 在线亚洲精品国产二区图片欧美 | av在线亚洲专区| 毛片女人毛片| 成人黄色视频免费在线看| 九九爱精品视频在线观看| 一本色道久久久久久精品综合| 亚洲丝袜综合中文字幕| 免费不卡的大黄色大毛片视频在线观看| 天美传媒精品一区二区| 人妻制服诱惑在线中文字幕| 成人无遮挡网站| 在线播放无遮挡| 国产免费视频播放在线视频| 国内精品美女久久久久久| 成年免费大片在线观看| 精华霜和精华液先用哪个| 男女那种视频在线观看| 成人漫画全彩无遮挡| 黄色一级大片看看| 久久97久久精品| 亚洲精华国产精华液的使用体验| 91在线精品国自产拍蜜月| 熟女电影av网| 久久久色成人| 久久久久久九九精品二区国产| 超碰97精品在线观看| 日产精品乱码卡一卡2卡三| 高清毛片免费看| 在线观看三级黄色| 男女无遮挡免费网站观看| 青春草亚洲视频在线观看| 伊人久久精品亚洲午夜| av国产精品久久久久影院| 日本-黄色视频高清免费观看| 简卡轻食公司| 亚洲精品日韩在线中文字幕| 极品教师在线视频| 国产精品久久久久久精品电影| 国产探花在线观看一区二区| 嫩草影院精品99| 国产一区二区三区综合在线观看 | 日本午夜av视频| 乱码一卡2卡4卡精品| 国产成人福利小说| 中国美白少妇内射xxxbb| 天天躁日日操中文字幕| 午夜福利网站1000一区二区三区| 韩国av在线不卡| 亚洲欧美成人综合另类久久久| 青春草视频在线免费观看| 综合色丁香网| 夜夜看夜夜爽夜夜摸|