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

    巴西橡膠樹HbGRF基因的克隆、物理定位及表達(dá)分析

    2020-10-29 07:35張宇航潘冉冉李飛陶志強(qiáng)王英高和瓊莊南生
    熱帶作物學(xué)報 2020年9期
    關(guān)鍵詞:表達(dá)分析

    張宇航 潘冉冉 李飛 陶志強(qiáng) 王英 高和瓊 莊南生

    摘 ?要:生長調(diào)控因子(growth-regulating factors,GRF)是植物特有的一類轉(zhuǎn)錄因子,該轉(zhuǎn)錄因子家族成員數(shù)目較多,在植物生長發(fā)育過程中發(fā)揮著重要作用。以巴西橡膠樹優(yōu)良品種‘熱研7-33-97為材料,通過RT-PCR方法進(jìn)行HbGRF基因的克隆;利用生物信息學(xué)的方法對其蛋白序列、理化性質(zhì)、基因結(jié)構(gòu)、進(jìn)化關(guān)系進(jìn)行分析;利用IS-PCR技術(shù)和FISH技術(shù)對其進(jìn)行物理定位分析;采用qRT-PCR對橡膠樹中HbGRF基因的表達(dá)模式進(jìn)行研究。結(jié)果表明:在橡膠樹雌花中克隆到3個GRF基因,分別命名為HbGRF1、HbGRF2和HbGRF3,分子量分別為65.663、41.188、52.858 kDa,其編碼的蛋白長度分別為609、384、494 aa,均為不穩(wěn)定蛋白;3個基因均具有GRF完整的特征結(jié)構(gòu)域WRC和QLQ,分別屬于3個不同亞族。基因物理定位結(jié)果表明,HbGRF1、HbGRF2和HbGRF3基因分別位于橡膠樹染色體的第11號長臂、第5號短臂和第9號長臂上,基因到著絲粒的平均百分距離是77.65、42.66和65.27。表達(dá)結(jié)果顯示,橡膠樹3個HbGRF基因在莖尖、雌花等生長旺盛的組織中表達(dá)量較高,用外源赤霉素和脫落酸處理后,發(fā)現(xiàn)表達(dá)量呈上升趨勢,經(jīng)過48 h后表達(dá)量和初始值基本持平,3個基因在莖尖、雌花中有明顯的表達(dá)特異性,可能在橡膠樹花芽分化及雌花的發(fā)育過程中起到了重要作用。該結(jié)果為研究橡膠樹HbGRF基因的功能及作用機(jī)制奠定了理論基礎(chǔ),為橡膠樹精準(zhǔn)育種提供了分子細(xì)胞學(xué)依據(jù)。

    關(guān)鍵詞:巴西橡膠樹;GRF基因家族;基因定位;表達(dá)分析

    中圖分類號:S794.1 ? ? ?文獻(xiàn)標(biāo)識碼:A

    Abstract: Growth-regulating factor (GRF) is a plant-specific transcription factor. The transcription factor family has a large number of members and plays an important role in plant growth and development. In this study, the excellent clone ‘Reyan 7-33-97 was used as the experimental material, and the HbGRF genes were cloned by RT-PCR. The bioinformatics method was used to study its gene length, physicochemical properties, basic structure and evolutionary relationship. The physical localization analysis was carried out by in situ PCR and fluorescence in situ hybridization. The expression pattern of HbGRF genes in rubber trees was analyzed by qRT-PCR. Three GRF genes were cloned into the female flowers of rubber trees, named HbGRF1, HbGRF2 and HbGRF3, respectively. The protein length was 609 aa, 384 aa and 494 aa with molecular weight of 65.663 kDa, 41.188 kDa and 52.858 kDa, respectively, and they were all unstable protein. All the three members had the complete characteristic domains WRC and QLQ of GRF gene family. Evolutionary analysis showed that the three genes belonged to three different groups. HbGRF1, HbGRF2 and HbGRF3 were located on the long arm of chromosome 11, the short arm of chromosome 5 and the long arm of chromosome 9, and the average percentage distance between the signal site and the centromere was 77.65, 42.66 and 65.27, respectively. The expression results showed that the three HbGRF genes in rubber trees were highly expressed in the vigorous growth tissues like stem tips and female flowers. Under the treatment of exogenous GA and ABA, the expression levels of the three genes increased first, and the expression levels at 48 h were basically the same as those under the initial conditions. The three genes had obvious expression specificity in stem tips and female flowers, which may play important roles in flower bud differentiation and female flower development. This study would provide a theoretical basis and materials for further exploration of the function of HbGRF, and provide molecular cytology basis for accurate breeding of rubber trees.

    Keywords: Hevea brasiliensis; GRF gene family; gene mapping; expression analysis

    DOI: 10.3969/j.issn.1000-2561.2020.09.001

    巴西橡膠樹(Hevea brasiliensis)屬于大戟科橡膠樹屬[1],其產(chǎn)生的膠乳在工業(yè)、國防等領(lǐng)域有著十分重要的作用。橡膠樹喜高溫、高濕的生長環(huán)境,易受極端天氣的影響,微弱的生長環(huán)境變化就會導(dǎo)致膠乳品質(zhì)和產(chǎn)量的下降[2]。生長調(diào)控因子(growth-regulating factors, GRF)是一類成員數(shù)目較多的、植物特有的轉(zhuǎn)錄因子家族,參與了植物生長發(fā)育等過程的調(diào)控[3-5]。本研究通過克隆橡膠樹HbGRF基因,并進(jìn)行了基因定位及表達(dá)分析,為研究橡膠樹HbGRF基因的功能及作用機(jī)制奠定了理論基礎(chǔ),也為橡膠樹精準(zhǔn)育種提供了理論指導(dǎo)。植物里最先被發(fā)現(xiàn)的GRF基因是水稻的OsGRF1,它編碼的蛋白通過調(diào)控赤霉素影響莖的伸長[6]。目前,已報道的GRF家族基因在擬南芥和水稻中分別有9個和12個成員[7-8],楊樹中有19個成員[9],番茄中有13個成員[10],陸地棉中有33個成員[11],油菜中有34個成員[12],煙草中有25個成員[13],桃樹中有10個成員[14],茶樹中有11個成員[15],藜麥中有18個成員[16]。該基因家族編碼的蛋白序列的N末端一般含有QLQ(Gln、Leu、Gln)和WRC(Trp、Arg、Cys)2個保守結(jié)構(gòu)域,QLQ結(jié)構(gòu)域與GRF互作因子(GRF-interacting factor, GIF)相互作用形成轉(zhuǎn)錄激活因子[17],WRC中包含一個與DNA相結(jié)合的鋅指結(jié)構(gòu)和核定位信號區(qū)域(NLS),它們在基因行使其生物學(xué)功能方面起到了重要作用[7]。在植物中,GRF基因家族的成員在發(fā)育活躍的部位如花芽、根尖及幼葉中表達(dá)量較高,在成熟的部位表達(dá)量相對較低[18-19],多種GRF基因在葉和莖的生長、花的發(fā)育、胚珠形成等方面發(fā)揮了重要的作用,GRF基因還可通過負(fù)反饋調(diào)節(jié)機(jī)制在植物滲透脅迫下發(fā)揮作用[7]。目前,在楊樹、桃樹、蘋果等許多植物中已對GRF基因家族進(jìn)行了鑒定及功能分析,上述研究完善了我們對GRF基因家族的認(rèn)識,但該家族在橡膠樹中的研究及功能分析尚未見報道。本研究經(jīng)克隆得到HbGRF基因,運(yùn)用生物信息學(xué)和基因表達(dá)等方法分析HbGRF及其蛋白的基本信息,初步了解HbGRF在橡膠樹生長發(fā)育過程中的生物學(xué)功能及作用機(jī)制;利用原位PCR和熒光原位雜交技術(shù)探究HbGRF與其他功能基因之間的位置關(guān)系,這可為橡膠樹精準(zhǔn)育種設(shè)計提供理論指導(dǎo)。

    1 ?材料與方法

    1.1 ?材料

    本研究以正常生長開割10 a的巴西橡膠樹優(yōu)良品種‘熱研7-33-97為材料,采集根、莖、莖尖、葉片、樹皮、果實(shí)、雌花以及不同發(fā)育時期的花序。實(shí)驗(yàn)材料由中國熱帶農(nóng)業(yè)科學(xué)院國家橡膠樹種質(zhì)資源圃提供。

    外源赤霉素(GA)和脫落酸(ABA)處理所用材料為1 a的‘熱研7-33-97組培苗,每個處理組設(shè)置5棵幼苗。GA處理:對植株噴施濃度為200 ?mol/L的GA,使用過程中以無水乙醇溶解GA后噴施處理植株,并以噴施等濃度乙醇處理的植株為對照,分別采集處理0、0.5、2、6、12、24、48 h植株與對照植株的葉片,ABA處理與GA處理方法相同。以上所有材料取樣后用液氮速凍,保存至-80 ℃冰箱中備用。

    1.2 ?方法

    1.2.1 ?橡膠樹HbGRF基因的克隆 ?根據(jù)GRF家族基因的保守序列,在橡膠樹基因組數(shù)據(jù)庫(http://hevea.catas.cn)中做同源序列搜索,通過CDS編碼區(qū)預(yù)測和序列拼接得到橡膠樹HbGRF1、HbGRF2、HbGRF3基因的cDNA序列。運(yùn)用Primer Primer 5.0軟件進(jìn)行引物設(shè)計(表1)。用反轉(zhuǎn)錄得到的橡膠樹雌花cDNA為模板,對橡膠樹HbGRF基因家族的HbGRF1、HbGRF2和HbGRF3 3個成員的編碼區(qū)進(jìn)行擴(kuò)增,PCR反應(yīng)體系為20 ?L,其中包括Prime STAR Max Premix(2×)10 ?L,cDNA模板1 ?L,上、下游引物各0.4 ?L,滅菌的ddH2O 8.2 ?L。PCR擴(kuò)增程序?yàn)椋?4 ℃預(yù)變性5 min;94 ℃變性30 s,56 ℃退火30 s,72 ℃延伸90 s,共擴(kuò)增35個循環(huán);最后于72 ℃延伸10 min。使用OMEGA公司的凝膠回收試劑盒對產(chǎn)物進(jìn)行回收純化后與pMD19-T載體進(jìn)行連接,轉(zhuǎn)化至大腸桿菌DH5α中,挑選陽性克隆重組子交華大基因公司測序。

    1.2.2 ?橡膠樹HbGRF基因的生物信息學(xué)分析 ?HbGRF基因的編碼區(qū)序列及氨基酸序列的比對和預(yù)測分析是通過NCBI中的BLAST軟件和ORF Finder軟件完成,其所編碼蛋白的理化特性通過ExPASy ProtParam在線分析[20]。HbGRF蛋白的亞細(xì)胞定位、信號肽、編碼蛋白的跨膜區(qū)段及功能結(jié)構(gòu)域預(yù)測,分別采用PSORT、SignalP、TMpred及SMART軟件進(jìn)行分析。使用PredictProtein軟件對HbGRF基因家族編碼蛋白的二級結(jié)構(gòu)進(jìn)行預(yù)測分析。使用DNAMAN 6.0軟件對近源物種的HbGRF蛋白的氨基酸序列進(jìn)行多重比對分析。使用MEGA 7.0軟件對HbGRF蛋白進(jìn)行系統(tǒng)進(jìn)化樹分析。

    1.2.3 ?橡膠樹HbGRF基因在染色體上的物理定位

    (1)染色體標(biāo)本的制備:參照高和瓊等[21]、李懋學(xué)等[22]和張志丹等[23]的實(shí)驗(yàn)方法,選取古銅期幼葉進(jìn)行染色體標(biāo)本的制備,于-20 ℃冰箱中保存?zhèn)溆谩?/p>

    (2)特異性引物的合成與篩選:依據(jù)克隆得到的測序結(jié)果,對序列片段進(jìn)行分析,參照代素珍等[24]的方法對HbGRF基因進(jìn)行合成與篩選,將PCR的退火溫度調(diào)整為45~65 ℃,通過凝膠電泳對PCR產(chǎn)物進(jìn)行檢測,選取條帶單一、明亮的產(chǎn)物交華大基因公司測序,將測序結(jié)果與基因序列比對,共篩選出3對符合條件的特異性引物(表1)。

    (3)原位PCR:參照官錦燕等[25]的原位PCR方法進(jìn)行改進(jìn),對HbGRF1和HbGRF2基因進(jìn)行檢測。

    (4)熒光原位雜交:將PCR擴(kuò)增得到的HbGRF3基因的擴(kuò)增產(chǎn)物回收純化后,用地高辛切口平移試劑盒DIG-Nick Translation Mix將其標(biāo)記成探針,信號位點(diǎn)呈現(xiàn)綠色;標(biāo)記好的探針溶于100%去離子甲酰胺溶液中,且終濃度為50 ng/?L;將標(biāo)記好的探針置于-20 ℃冰箱中避光保存?zhèn)溆?。熒光原位雜交的方法參照本課題組邱海燕[26]報道的方法進(jìn)行檢測。

    (5)鏡檢與分析:將40 ?L含有DAPI染料的抗褪色劑滴加至染色體標(biāo)本上進(jìn)行封片,在黑暗環(huán)境中放置60 min后,使用Leica DM5000顯微鏡進(jìn)行觀察與拍照;用Photoshop軟件對所拍照片進(jìn)行處理,使用Image J軟件對染色體長度及信號位點(diǎn)距著絲粒的距離進(jìn)行測量,并利用擴(kuò)增信號位點(diǎn)百分距離的計算方法[27],結(jié)合‘熱研7-33-97的核型分析參數(shù)[28],分析和確定信號位點(diǎn)在染色體上的位置。

    擴(kuò)增信號位點(diǎn)的百分距離=?×100

    1.2.4 ?橡膠樹HbGRF基因的表達(dá)分析 ?依據(jù)橡膠樹HbGRF基因的ORF區(qū)設(shè)計熒光定量引物,使用TIANGEN試劑盒提取RNA,用TaKaRa的反轉(zhuǎn)錄試劑盒合成cDNA第1條鏈,以各樣品的cDNA為模板,濃度統(tǒng)一稀釋為100 ng/?L,進(jìn)行定量PCR擴(kuò)增。設(shè)計熒光定量特異性引物見表1,對HbGRF基因在不同組織、不同發(fā)育時期及不同外源激素處理后的表達(dá)情況進(jìn)行qRT-PCR分析,以巴西橡膠樹18S rRNA為內(nèi)參,采用2ΔΔCT法進(jìn)行定量分析,每個試驗(yàn)設(shè)置3次重復(fù),定量數(shù)據(jù)使用SAS 9.1.3軟件進(jìn)行統(tǒng)計分析[29]。

    2 ?結(jié)果與分析

    2.1 ?橡膠樹HbGRF基因的克隆和蛋白序列分析

    以巴西橡膠樹‘熱研7-33-97雌花cDNA為模板,使用特異性引物克隆得到3個片段大小分別為1830、1155、1485 bp的基因。測序后發(fā)現(xiàn)這3個基因片段和轉(zhuǎn)錄組測序后的序列一致,將這3個基因分別命名為HbGRF1、HbGRF2和HbGRF3,其登錄號分別為:MN396751、MN396752和MN396753。

    運(yùn)用生物信息學(xué)方法對巴西橡膠樹3個HbGRF基因進(jìn)行分析,分析結(jié)果表明其編碼384~609個氨基酸,分子量大小范圍為41.188~ 65.663 kDa,等電點(diǎn)為6.09~8.77,蛋白質(zhì)不穩(wěn)定指數(shù)為51.08~56.06,均為不穩(wěn)定蛋白,脂肪系數(shù)為57.98~69.03,疏水性平均系數(shù)范圍為0.541~ 0.582(表2)。對HbGRF蛋白進(jìn)行亞細(xì)胞定位分析,發(fā)現(xiàn)GRF蛋白均處在細(xì)胞核內(nèi),與轉(zhuǎn)錄因子的特征相符合。

    通過Predict Protein在線軟件對HbGRF1、HbGRF2、HbGRF3蛋白的二級結(jié)構(gòu)進(jìn)行預(yù)測,結(jié)果顯示這3個蛋白都含有DNA結(jié)合區(qū)和多個蛋白結(jié)合區(qū)域(圖1),說明它們可能在轉(zhuǎn)錄過程中發(fā)揮作用。

    紅色條帶:螺旋結(jié)構(gòu);藍(lán)色條帶:鏈;黃色條帶:隱藏結(jié)構(gòu);綠色條帶:無規(guī)律結(jié)構(gòu);紅色小方塊:蛋白結(jié)合區(qū);黃色小圓圈:DNA結(jié)合區(qū);紫色小圓圈:RNA結(jié)合區(qū);灰色小圓圈:核苷酸結(jié)合區(qū)。

    2.2 ?橡膠樹HbGRF蛋白的多重比對及進(jìn)化分析

    對HbGRF的蛋白序列和近緣物種如木薯(Manihot esculenta)、可可樹(Theobroma cacao)以及麻風(fēng)樹(Jatropha curcas)的GRF蛋白序列進(jìn)行比對。結(jié)果表明,橡膠中的3個HbGRF基因家族成員均含有QLQ和WRC保守結(jié)構(gòu)域,QLQ結(jié)構(gòu)域的序列是QX3LX2Q,WRC結(jié)構(gòu)域的序列是CX9CX10CX2H,在橡膠HbGRF家族的3個成員中,WRC結(jié)構(gòu)域比QLQ結(jié)構(gòu)域更為保守(圖2),表明這3個橡膠樹HbGRF蛋白均含有特有的GRF結(jié)構(gòu)域。

    進(jìn)一步選取木薯、可可樹、水稻、擬南芥、麻風(fēng)樹等已知的GRF家族基因使用MEGA 6.1軟件構(gòu)建系統(tǒng)進(jìn)化樹。結(jié)果顯示,在植物中GRF家族基因可以分成5個亞族(Ⅰ~Ⅴ)(圖3)。其中橡膠樹的HbGRF1、HbGRF2、HbGRF3基因分別屬于Ⅰ亞族、Ⅱ亞族、Ⅳ亞族,且HbGRF1與HbGRF3基因親緣關(guān)系略近。

    2.3 ?HbGRF1、HbGRF2和HbGRF3基因的物理定位分析

    分別用HbGRF1和HbGRF2基因定位的特異引物與Digoxigenin-11-dUTP等配制成擴(kuò)增液,在制備好的染色體標(biāo)本上進(jìn)行原位PCR反應(yīng)和熒光檢測。在有絲分裂的間、前、中期的染色體上均能觀察到2個綠色熒光信號點(diǎn)(圖4,圖5)。陰性對照組中不加入Taq酶,其細(xì)胞染色體上均未出現(xiàn)信號(圖4D,圖5D)。對信號位點(diǎn)所在的中期細(xì)胞進(jìn)行核型分析,得到信號位點(diǎn)與染色體的對應(yīng)關(guān)系,以及信號位點(diǎn)在核型模式圖上的位置關(guān)系。結(jié)果表明,HbGRF1和HbGRF2基因分別位于第11號染色體長臂和第5號染色體短臂上,基因到著絲粒的平均百分距離分別是77.65和42.66。

    用地高辛對HbGRF3基因的擴(kuò)增純化產(chǎn)物標(biāo)記為探針,隨后采用熒光原位雜交進(jìn)行檢測,在有絲分裂的間、前、中期的染色體上均能觀察到2個綠色熒光信號點(diǎn)(圖6),陰性對照組中不加入探針,其細(xì)胞染色體上均未出現(xiàn)信號(圖6D)。對信號位點(diǎn)所在的中期細(xì)胞進(jìn)行核型分析,得到信號位點(diǎn)與染色體的對應(yīng)關(guān)系,以及信號位點(diǎn)在核型模式圖上的位置關(guān)系(圖7)。結(jié)果表明,HbGRF3位于第9號染色體長臂上,基因到著絲粒的平均百分距離是65.27。

    2.4 ?橡膠樹HbGRF基因的表達(dá)分析

    2.4.1 ?橡膠樹HbGRF基因在不同組織中的表達(dá)分析 ?為進(jìn)一步研究橡膠樹HbGRF基因在生長發(fā)育過程中的作用,本研究使用qRT-PCR技術(shù)對橡膠樹根、莖、莖尖等8個組織或器官中HbGRF基因的相對表達(dá)量進(jìn)行分析。由圖8可知,HbGRF1、HbGRF2和HbGRF3基因在植株8個不同部位中表達(dá)量有差異。其中,HbGRF1基因在莖尖、雌花中高調(diào)表達(dá),在其他部位表達(dá)量較低。HbGRF2和HbGRF3基因在莖尖中顯著高調(diào)表達(dá),在雌花中的表達(dá)量也比別的部位高。由此推測HbGRF1、HbGRF2和HbGRF3基因可能參與橡膠樹花芽分化及雌花的發(fā)育過程。

    2.4.2 ?橡膠樹HbGRF基因在不同發(fā)育時期花序中的表達(dá)分析 ?對橡膠樹HbGRF基因在花蕾期(花序長度約為3 cm)、伸展期(花序長度約為6 cm)、盛花期(花序長度約為9 cm)和開花末期(花序長度約為12 cm)等4個花序發(fā)育階段的表達(dá)情況進(jìn)行分析,結(jié)果表明,在伸展期,HbGRF1、HbGRF2和HbGRF3基因的表達(dá)量最高,由此推測HbGRF基因在生長發(fā)育較為活躍的部位表達(dá)較強(qiáng)(圖9)。

    2.4.3 ?橡膠樹HbGRF基因在外源激素處理下的表達(dá)分析 ?使用外源赤霉素對橡膠樹幼葉進(jìn)行處理后,HbGRF1和HbGRF2在處理后0.5 h、HbGRF3在處理后6 h時,其表達(dá)量最高;使用外源脫落酸處理時發(fā)現(xiàn),3個基因均在處理后0.5 h時,表達(dá)量達(dá)到了最高值。當(dāng)2種外源激素處理后48 h時,3個基因的表達(dá)量和初始條件下的表達(dá)量基本相同(圖10)。

    3 ?討論

    GRF基因家族是一類成員數(shù)目較多的轉(zhuǎn)錄因子家族,該轉(zhuǎn)錄因子對植物的生長發(fā)育過程具有重要的調(diào)控作用。目前,GRF基因家族在水稻、擬南芥、楊樹、番茄、陸地棉、油菜、煙草、桃樹、茶樹、藜麥等多種植物中已有研究,但對于具有重要戰(zhàn)略意義的橡膠樹GRF基因的相關(guān)研究仍缺乏進(jìn)一步的報道。本研究從橡膠樹中克隆到3個HbGRF基因,首先對其進(jìn)行生物信息學(xué)與表達(dá)模式等分析,再使用原位PCR技術(shù)和熒光原位雜交技術(shù)對基因進(jìn)行了物理定位分析,試圖從分子生物學(xué)及細(xì)胞學(xué)的角度對功能基因進(jìn)行較為全面的研究。

    伴隨基因組測序技術(shù)的發(fā)展與完善,Livak等[29]對橡膠樹‘熱研7-33-97的全基因組進(jìn)行了測序,這對橡膠樹基因家族的分析與鑒定、功能基因的挖掘與研究都提供了很多便利,預(yù)測了43 792個蛋白編碼基因,其中包括7個預(yù)測的橡膠樹GRF基因,較為符合一般植物中GRF基因家族含有8~20個成員的規(guī)律。本研究根據(jù)擬南芥GRF的2個保守序列QLQ和WRC,在橡膠數(shù)據(jù)庫搜索到7個基團(tuán),然后設(shè)計引物,以橡膠樹總cDNA為模板進(jìn)行擴(kuò)增,結(jié)果顯示,只擴(kuò)增出3個基因,推測可能是其他基因的豐度比較低等原因造成。將本研究中擴(kuò)增的3個HbGRF基因和預(yù)測的7個橡膠樹GRF基因的序列比對,結(jié)果發(fā)現(xiàn)HbGRF1配對序列的一致性為99.84%,其余均完全比對,說明該比對結(jié)果與Tang等[30]的預(yù)測結(jié)果一致。

    對3個HbGRF基因的氨基酸序列進(jìn)行分析,結(jié)果發(fā)現(xiàn)3個成員均含有QLQ和WRC 2個保守結(jié)構(gòu)域,且WRC結(jié)構(gòu)域的特異性較高,這與藜麥、茶樹等植物中的研究相一致。對橡膠樹與其他植物GRF的系統(tǒng)進(jìn)化分析表明,可將GRF分為5個亞族,其中HbGRF1與MeGRF1、HbGRF2與JcGRF3、HbGRF3與AtGRF8親緣關(guān)系較近,基因之間的親緣關(guān)系越接近,其功能可能也越類似。如擬南芥中的AtGRF8基因在花的發(fā)育中具有重要的調(diào)控作用,尤其是在花器官發(fā)育后期。HbGRF3基因在橡膠樹雌花中的表達(dá)量相對較高,由此推測橡膠樹HbGRF3基因在花的發(fā)育過程中起重要的調(diào)控作用。

    前人研究結(jié)果表明,GRF基因家族的成員在發(fā)育活躍的器官或組織中表達(dá)量較高,在成熟的組織或器官中表達(dá)量相對較低。在茶樹中CsGRF在生長比較旺盛的組織尤其是芽葉中起到了重要的調(diào)控作用[15]。在油菜中BnGRF可調(diào)控花器官的發(fā)育[31]。本研究對3個HbGRF基因在不同部位的表達(dá)情況分析發(fā)現(xiàn),它們在莖尖及雌花中的表達(dá)量相對較高,推測這3個基因參與了調(diào)控橡膠樹花芽分化以及雌花發(fā)育的過程。對HbGRF在花序不同發(fā)育時期的表達(dá)量進(jìn)行分析發(fā)現(xiàn),其在伸展期表達(dá)量最高,符合GRF轉(zhuǎn)錄因子的特點(diǎn),但在開花末期,3個基因的表達(dá)量也相對較高,出現(xiàn)這一現(xiàn)象的原因及機(jī)制尚不明確,需要進(jìn)行更深入的研究。在實(shí)際生產(chǎn)實(shí)踐過程中,橡膠樹座果率較低,一般為5%~8%,最高不超過10%,這是橡膠樹常規(guī)育種工作效率較低的原因之一。赤霉素信號途徑在調(diào)控植物株高、開花方面有著重要的作用,脫落酸信號途徑在植物多種生理反應(yīng)及基因的脅迫應(yīng)答方面具有重要的調(diào)控作用。前人的研究發(fā)現(xiàn),蘋果樹MdGRF基因可能在外源激素作用下調(diào)控蘋果樹的生長,在調(diào)控嫁接到矮化砧木上接穗的生長以及2個品種節(jié)間伸長等過程中有著重要的作用[32]。本研究用外源赤霉素和外源脫落酸對橡膠樹進(jìn)行處理后發(fā)現(xiàn),其表達(dá)量在處理后0.5~6 h內(nèi)達(dá)到最高,在外源激素處理后48 h時,3個基因的表達(dá)量回落到初始水平,表明HbGRF1、HbGRF2和HbGRF3受赤霉素和脫落酸處理的調(diào)控。

    基因定位的結(jié)果表明,HbGRF1、HbGRF2和HbGRF3定位在3條不同的染色體上,它們之間可能相互獨(dú)立,結(jié)合本課題組已報道的橡膠樹其他功能基因的定位結(jié)果可以得知,HbGRF1與OVP-10390[33]、Hevb7[34]同時位于第11號染色體上;HbGRF2與HbMyb1[35]、HbN1N1[36]、HbNAC1[37]、HbPEPRK2[24]、HblMAC4[38]同時位于第5號染色體上;HbGRF3與HbCPT[39]、HbJAZ5[40]、HbNAC3[37]、HblMYC2、HblMYC5[38]、HbWRKY1、HbWRKY4、HbWRKY5[41]同時位于第9號染色體上。這些定位在相同染色體上的基因間可能存在互為連鎖的關(guān)系。研究功能基因之間的相互關(guān)系對了解這些基因間可能存在的遺傳關(guān)系提供了幫助,這些結(jié)果可以為橡膠樹精準(zhǔn)育種設(shè)計提供理論指導(dǎo)。

    參考文獻(xiàn)

    [1] 劉少軍, 周廣勝, 房世波.中國橡膠樹種植氣候適宜性區(qū)劃[J]. 中國農(nóng)業(yè)科學(xué), 2015, 48(12): 2335-2345.

    [2] 劉金河. 巴西橡膠樹的水分狀況與生長和產(chǎn)膠量的關(guān)系[J]. 生態(tài)學(xué)報, 1982(3): 217-224.

    [3] Chen Y Z, Cao J. Comparative analysis of Dof transcription factor family in maize[J]. Plant Molecular Biology Reporter, 2015, 33(5): 1245-1258.

    [4] Hoe Kim J, Tsukaya H. Regulation of plant growth and development by the GROWTH-REGULATING FACTOR and GRF-INTERACTING FACTOR duo[J]. Journal of Experimental Botany, 2015, 66(20): 6093-6107.

    [5] Schwechheimer C, Bevan M. The regulation of transcription factor activity in plants[J]. Trends in Plant Science, 1998, 3(10): 378-383.

    [6] Knaap E V D, Kende K H. A novel gibberellin-induced gene from rice and its potential regulatory role in stem growth[J]. Plant Physiology, 2000, 122(3): 695-704.

    [7] Kim J H, Choi D, Kende H. The AtGRF family of putative transcription factors is involved in leaf and cotyledon growth in Arabidopsis[J]. The Plant Journal, 2003, 36(1): 94-104.

    [8] Dongsu C, Hoe K J, Hans K. Whole genome analysis of the OsGRF gene family encoding plant-specific putative tran-scription activators in rice (Oryza sativa L.)[J]. Plant and Cell Physiology, 2004, 45(7): 897-904.

    [9] 周厚君. 楊樹GRF基因家族分析及PtGRF1/2d功能研究[D]. 北京: 中國林業(yè)科學(xué)研究院, 2016: 26-35.

    [10] 袁 ?岐, 張春利, 趙婷婷, 等. 番茄GRF轉(zhuǎn)錄因子家族的生物信息學(xué)分析[J]. 分子植物育種, 2017, 15(8): 2949-2956.

    [11] 張書芹, 樂 ?愉, 武 ?斐. 陸地棉GRF基因家族的鑒定和生物信息學(xué)分析[J]. 分子植物育種, 2019, 17(12): 3817-3824.

    [12] 阮先樂, 王俊生, 劉紅占, 等. 油菜GRF基因家族的鑒定和基本特征分析[J]. 分子植物育種, 2018, 16(8): 2420-2428.

    [13] Zhang J, Li Z, Jin J, et al. Genome-wide identification and analysis of the growth-regulating factor family in tobacco (Nicotiana tabacum)[J]. Gene, 2018, 639: 117-127.

    [14] 曹 ?珂, 薛靈姿, 王 ?蛟, 等. 桃GRF基因家族的序列及其組織特異性表達(dá)分析[J]. 植物遺傳資源學(xué)報, 2018, 19(3): 578-586.

    [15] 王鵬杰, 鄭玉成, 林 ?浥, 等. 茶樹GRF基因家族的全基因組鑒定及表達(dá)分析[J]. 西北植物學(xué)報, 2019, 39(3): 413-421.

    [16] 時丕彪, 何 ?冰, 費(fèi)月躍, 等. 藜麥GRF轉(zhuǎn)錄因子家族的鑒定及表達(dá)分析[J]. 作物學(xué)報, 2019, 45(12): 1841-1850.

    [17] Kim J H, Kende H. A transcriptional coactivator, AtGIF1, is involved in regulating leaf growth and morphology in Arabidopsis[J]. Proceedings of the National Academy of Sciences of the United States of America, 2004, 101(36): 13374-13379.

    [18] Ma J Q, Jian H J, Yang B, et al. Genome-wide analysis and expression profiling of the GRF gene family in oilseed rape (Brassica napus L.)[J]. Gene, 2017, 620: 36-45.

    [19] Zheng L W, Zhao C D, Mao J P, et al. Genome-wide identification and expression analysis of brassinosteroid biosynthesis and metabolism genes regulating apple tree shoot and lateral root growth[J]. Journal of Plant Physiology, 2018, 231: 68-85.

    [20] Elisabeth G, Alexandre G, Christine H, et al. ExPASy: The proteomics server for in-depth protein knowledge and analysis[J]. Nucleic Acids Research, 2003, 31(13): 3784-3788.

    [21] 高和瓊, 王 ?英, 金 ?鴿, 等. 橡膠樹葉片染色體制片方法的優(yōu)化[J]. 熱帶作物學(xué)報, 2009, 30(5): 565-569.

    [22] 李懋學(xué), 張贊平. 作物染色體及其研究技術(shù)[M]. 北京:中國農(nóng)業(yè)出版社, 1996: 30-34.

    [23] 張志丹, 王 ?英, 莊南生, 等. 巴西橡膠樹染色體制片方法的改良及FISH信號檢測[J]. 基因組學(xué)與應(yīng)用生物學(xué), 2019, 38(6): 2713-2717.

    [24] 代素珍, 莊南生, 王 ?英, 等. 巴西橡膠樹HbPEPRK基因家族5個成員在染色體上的定位[J]. 基因組學(xué)與應(yīng)用生物學(xué), 2019, 38(7): 3155-3160.

    [25] 官錦燕, 王 ?英, 高和瓊, 等. 巴西橡膠樹4個環(huán)鋅指蛋白基因(HbRZF)的物理定位[J]. 基因組學(xué)與應(yīng)用生物學(xué), 2014, 33(3): 610-616.

    [26] 邱海燕. 巴西橡膠樹REF、RT和SRPP基因物理定位的研究[D]. ??冢?海南大學(xué), 2010.

    [27] Song Y C, Gustafson J P J T, Genetics A. The physical location of fourteen RFLP markers in rice (Oryza sativa L.) [J]. Theoretical and Applied Genetics, 1995, 90(1): 113-119.

    [28] 高和瓊, 莊南生, 王 ?英, 等. 橡膠樹兩個品種的核型分析[J]. 武漢植物學(xué)研究, 2009, 27(5): 537-540.

    [29] Livak K J, Schmittgen T D. Analysis of relative gene ex-pression data using real-time quantitative PCR and the 2–CT method[J]. Methods, 2000, 25(4): 402-408.

    [30] Tang C R, Yang M, Fang Y J, et al. The rubber tree genome reveals new insights into rubber production and species ad-aptation[J]. Nature Plants, 2016, 2(6): 16073

    [31] Ma J Q, Jian H J, Yang B, et al. Genome-wide analysis and expression profiling of the, GRF, gene family in oilseed rape (Brassica napus L.) [J]. Gene, 2017, 620: 36-45.

    [32] Zheng L W, Ma J J, Song C H, et al. Genome-wide identification and expression analysis of GRF genes regulating apple tree architecture[J]. Tree Genetics & Genomes, 2018, 14(4): 54-60.

    [33] 高和瓊. 巴西橡膠樹HbMyb1基因和OPV-10390連鎖標(biāo)記原位PCR定位的研究[D]. ??冢?海南大學(xué), 2008: 26-44.

    [34] 高佳佳. 巴西橡膠樹幾個與橡膠合成相關(guān)基因的物理定位研究[D]. ??冢?海南大學(xué), 2014: 26-28.

    [35] 高和瓊, 莊南生, 王 ?英, 等. 巴西橡膠樹HbMyb1基因的原位PCR物理定位[J]. 熱帶亞熱帶植物學(xué)報, 2012, 20(4): 365-368.

    [36] 高佳佳, 王 ?英, 高和瓊, 等. 巴西橡膠樹膠乳轉(zhuǎn)化酶HbNIN基因家族物理定位的研究[J]. 熱帶作物學(xué)報, 2014, 35(9): 1704-1709.

    [37] 楊光涌, 鄭 ?菲, 王 ?英, 等. 巴西橡膠樹NAC基因家族5個成員的熒光原位雜交物理定位[J]. 分子植物育種, 2018, 16(2): 512-517.

    [38] 高 ?豫, 莊南生, 王 ?英, 等. MYC基因家族5成員在巴西橡膠樹染色體上的物理定位[J]. 熱帶生物學(xué)報, 2018, 9(2): 163-169.

    [39] 張新新. 巴西橡膠樹4個膠乳生物合成相關(guān)基因的物理定位研究[D]. 海口: 海南大學(xué), 2013: 33-45.

    [40] 李曉燕, 蘇莉莉, 高佳佳, 等. JAZ基因家族6個成員在橡膠樹上的物理定位[J]. 分子植物育種, 2017, 15(12): 4992-4999.

    [41] 劉正林, 莊南生, 王 ?英, 等. 巴西橡膠樹HbWRKY基因家族10個成員在染色體上的定位[J]. 植物遺傳資源學(xué)報, 2018, 19(6): 1170-1179.

    猜你喜歡
    表達(dá)分析
    巴西橡膠樹CCaMK基因家族成員的鑒定和表達(dá)分析
    太子參玉米黃質(zhì)環(huán)氧化酶基因的克隆與表達(dá)分析
    木薯PHOR1基因的序列分析及其乙烯和茉莉酸甲酯誘導(dǎo)表達(dá)特性
    雷公藤貝殼杉烯酸氧化酶基因的全長cDNA克隆與表達(dá)分析
    紅花生育酚環(huán)化酶基因的克隆及表達(dá)分析
    三色堇DFR基因的克隆及表達(dá)分析
    香蕉抗壞血酸過氧化物酶基因的克隆及表達(dá)分析
    膠孢炭疽菌漆酶基因Lac2的序列特征與表達(dá)分析
    玉米紋枯病病菌y—谷氨酰轉(zhuǎn)肽酶基因克隆與表達(dá)分析
    信號分子與葉銹菌誘導(dǎo)下小麥病程相關(guān)蛋白1基因的表達(dá)分析
    亚洲久久久久久中文字幕| 我要看日韩黄色一级片| 国产精品久久久久久亚洲av鲁大| 国产欧美日韩精品一区二区| 成年免费大片在线观看| 制服丝袜大香蕉在线| 香蕉av资源在线| 欧美一区二区精品小视频在线| 亚洲一区二区三区不卡视频| 免费看日本二区| 又爽又黄无遮挡网站| 一个人观看的视频www高清免费观看| 国产91精品成人一区二区三区| 亚洲一区二区三区色噜噜| 人人妻人人看人人澡| 校园春色视频在线观看| 国产亚洲av嫩草精品影院| a级毛片a级免费在线| 最新中文字幕久久久久| 亚洲国产高清在线一区二区三| 91九色精品人成在线观看| 校园春色视频在线观看| 1024手机看黄色片| 深爱激情五月婷婷| 他把我摸到了高潮在线观看| 日韩免费av在线播放| 国产毛片a区久久久久| 国产在线精品亚洲第一网站| 三级毛片av免费| av视频在线观看入口| 88av欧美| 18禁黄网站禁片免费观看直播| 欧美中文日本在线观看视频| 日韩av在线大香蕉| 亚洲熟妇熟女久久| 一级a爱片免费观看的视频| 女同久久另类99精品国产91| 成人一区二区视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 久久九九热精品免费| 一个人免费在线观看的高清视频| 国产高清激情床上av| 久久久久久九九精品二区国产| 久久欧美精品欧美久久欧美| 伦理电影大哥的女人| 欧美国产日韩亚洲一区| 国产淫片久久久久久久久 | 国产私拍福利视频在线观看| 成人国产一区最新在线观看| 一个人看的www免费观看视频| 精品久久久久久久久久免费视频| 老司机午夜福利在线观看视频| 久久久国产成人免费| 免费看美女性在线毛片视频| 五月玫瑰六月丁香| 黄色女人牲交| 亚洲 国产 在线| 亚洲人成网站在线播放欧美日韩| 人人妻人人看人人澡| 国产美女午夜福利| 18禁在线播放成人免费| 一进一出好大好爽视频| 午夜免费成人在线视频| 国产激情偷乱视频一区二区| 久久精品国产99精品国产亚洲性色| 久久久久九九精品影院| av在线天堂中文字幕| 国产精品久久久久久精品电影| 欧美一区二区精品小视频在线| h日本视频在线播放| 国产精品野战在线观看| 国产激情偷乱视频一区二区| 久久久久久久亚洲中文字幕 | 亚洲专区中文字幕在线| 又黄又爽又免费观看的视频| 日日摸夜夜添夜夜添小说| 精品99又大又爽又粗少妇毛片 | 日本免费a在线| 久久这里只有精品中国| 亚洲 国产 在线| 婷婷亚洲欧美| 国产大屁股一区二区在线视频| 好男人在线观看高清免费视频| 亚洲综合色惰| 日本黄色视频三级网站网址| 赤兔流量卡办理| 女人被狂操c到高潮| 久久午夜福利片| 欧美性猛交╳xxx乱大交人| 久久精品影院6| 村上凉子中文字幕在线| 91九色精品人成在线观看| 欧美极品一区二区三区四区| 蜜桃亚洲精品一区二区三区| 夜夜躁狠狠躁天天躁| 成人国产综合亚洲| 两个人的视频大全免费| 国产午夜福利久久久久久| 国产男靠女视频免费网站| 欧美+亚洲+日韩+国产| 日韩欧美精品v在线| 久久久成人免费电影| 国产私拍福利视频在线观看| 国产一区二区激情短视频| 亚洲乱码一区二区免费版| 国产色婷婷99| 最近最新免费中文字幕在线| 五月伊人婷婷丁香| 国产av麻豆久久久久久久| 1024手机看黄色片| 深爱激情五月婷婷| 亚洲精品一区av在线观看| 老司机深夜福利视频在线观看| 久久久久免费精品人妻一区二区| 人人妻人人看人人澡| 国产精品野战在线观看| 免费在线观看日本一区| 中文字幕av在线有码专区| 99国产精品一区二区三区| 成年版毛片免费区| 在线国产一区二区在线| 国产精品乱码一区二三区的特点| 亚洲成人免费电影在线观看| 免费看光身美女| 欧美乱色亚洲激情| 91九色精品人成在线观看| 中文字幕av在线有码专区| 日韩欧美在线乱码| 婷婷亚洲欧美| .国产精品久久| 欧美不卡视频在线免费观看| 欧美不卡视频在线免费观看| 色哟哟哟哟哟哟| 波多野结衣高清无吗| 搡女人真爽免费视频火全软件 | 露出奶头的视频| 制服丝袜大香蕉在线| 久久久久九九精品影院| 国内精品久久久久久久电影| 午夜福利视频1000在线观看| 免费看美女性在线毛片视频| 国产精品av视频在线免费观看| 五月伊人婷婷丁香| 精品人妻偷拍中文字幕| 亚洲色图av天堂| 久久香蕉精品热| 欧美日本视频| 午夜免费男女啪啪视频观看 | 黄色女人牲交| 亚洲熟妇中文字幕五十中出| 免费看美女性在线毛片视频| 在线观看午夜福利视频| 日本撒尿小便嘘嘘汇集6| 国产av一区在线观看免费| 亚洲中文日韩欧美视频| 少妇高潮的动态图| 在线a可以看的网站| 看十八女毛片水多多多| 日韩欧美 国产精品| 日韩精品青青久久久久久| 久久亚洲精品不卡| 天堂动漫精品| 午夜福利免费观看在线| 亚洲国产欧美人成| 宅男免费午夜| 一a级毛片在线观看| 亚洲片人在线观看| 搡女人真爽免费视频火全软件 | 欧美黑人巨大hd| 免费观看的影片在线观看| 在现免费观看毛片| bbb黄色大片| 免费看日本二区| 观看免费一级毛片| 淫妇啪啪啪对白视频| 国内精品久久久久精免费| 欧美最黄视频在线播放免费| 国产欧美日韩精品一区二区| 免费在线观看日本一区| 天堂av国产一区二区熟女人妻| 九色国产91popny在线| 夜夜躁狠狠躁天天躁| 99热这里只有是精品50| 天堂av国产一区二区熟女人妻| 9191精品国产免费久久| 欧美绝顶高潮抽搐喷水| 亚洲av二区三区四区| 成年女人看的毛片在线观看| 欧美又色又爽又黄视频| 久久久久久久久大av| 久久人人爽人人爽人人片va | 毛片一级片免费看久久久久 | 精品人妻视频免费看| 波多野结衣高清无吗| 他把我摸到了高潮在线观看| 国产一区二区激情短视频| 免费人成视频x8x8入口观看| 美女高潮的动态| 99久久精品一区二区三区| 亚洲专区国产一区二区| 国内少妇人妻偷人精品xxx网站| 一a级毛片在线观看| avwww免费| 久久精品国产亚洲av香蕉五月| 亚洲国产欧洲综合997久久,| 欧美成人性av电影在线观看| 无人区码免费观看不卡| 内地一区二区视频在线| 久久伊人香网站| 久久久精品大字幕| 韩国av一区二区三区四区| 国产白丝娇喘喷水9色精品| 色哟哟哟哟哟哟| 国产精品久久久久久久电影| 亚洲 欧美 日韩 在线 免费| 91字幕亚洲| 91av网一区二区| 精品久久久久久久人妻蜜臀av| 尤物成人国产欧美一区二区三区| 成人美女网站在线观看视频| 亚洲成av人片在线播放无| 日韩精品中文字幕看吧| 亚洲 欧美 日韩 在线 免费| 亚洲精品亚洲一区二区| www日本黄色视频网| 亚州av有码| 欧美性猛交黑人性爽| 性色av乱码一区二区三区2| 一区二区三区免费毛片| 欧美日韩福利视频一区二区| 成年女人毛片免费观看观看9| 黄色丝袜av网址大全| 成人欧美大片| 欧美绝顶高潮抽搐喷水| 午夜两性在线视频| 搡老岳熟女国产| 一进一出好大好爽视频| 黄色视频,在线免费观看| 99久久无色码亚洲精品果冻| 亚洲欧美激情综合另类| 啦啦啦观看免费观看视频高清| 精品久久久久久久久久免费视频| 99久久精品一区二区三区| 99久久九九国产精品国产免费| 熟女电影av网| 国产蜜桃级精品一区二区三区| 丰满乱子伦码专区| 日本a在线网址| 亚洲av一区综合| 女生性感内裤真人,穿戴方法视频| 欧美一区二区亚洲| 麻豆成人午夜福利视频| 国产在视频线在精品| 欧美+日韩+精品| 日本熟妇午夜| 国产午夜福利久久久久久| 搡女人真爽免费视频火全软件 | 亚洲成人久久性| 一a级毛片在线观看| 一夜夜www| 亚洲欧美日韩高清专用| 18禁黄网站禁片午夜丰满| 欧美色视频一区免费| 色视频www国产| 内地一区二区视频在线| 国产 一区 欧美 日韩| 日日夜夜操网爽| 国产伦一二天堂av在线观看| 一区二区三区激情视频| 两个人的视频大全免费| 18禁黄网站禁片免费观看直播| 久久这里只有精品中国| 少妇裸体淫交视频免费看高清| 最近在线观看免费完整版| 欧美zozozo另类| 最近中文字幕高清免费大全6 | 99久久成人亚洲精品观看| 中文字幕免费在线视频6| 国产aⅴ精品一区二区三区波| 日韩大尺度精品在线看网址| 欧美日韩国产亚洲二区| 欧美日韩综合久久久久久 | 1024手机看黄色片| 国产精品伦人一区二区| 最近最新中文字幕大全电影3| 日日摸夜夜添夜夜添av毛片 | 脱女人内裤的视频| 国产熟女xx| 在线观看免费视频日本深夜| 亚洲真实伦在线观看| www.色视频.com| а√天堂www在线а√下载| 国产精品美女特级片免费视频播放器| 首页视频小说图片口味搜索| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品久久电影中文字幕| 欧美最黄视频在线播放免费| 久久久久国产精品人妻aⅴ院| 国产高潮美女av| 免费看a级黄色片| 国产av不卡久久| 性插视频无遮挡在线免费观看| 99热这里只有是精品在线观看 | 国产高清视频在线播放一区| 日韩欧美 国产精品| 亚洲自拍偷在线| 丰满乱子伦码专区| 亚洲国产色片| 亚洲在线观看片| 又黄又爽又免费观看的视频| 亚洲精品在线观看二区| 舔av片在线| 又爽又黄无遮挡网站| 国产精品乱码一区二三区的特点| 精品久久久久久,| a级毛片免费高清观看在线播放| 禁无遮挡网站| 国产精品永久免费网站| 黄色日韩在线| 久久久久久久久久成人| 亚洲国产日韩欧美精品在线观看| 精品人妻偷拍中文字幕| 天堂动漫精品| 男女视频在线观看网站免费| 99热这里只有精品一区| 熟女电影av网| 最新中文字幕久久久久| 制服丝袜大香蕉在线| 桃色一区二区三区在线观看| a级一级毛片免费在线观看| 日韩欧美 国产精品| 亚洲精品粉嫩美女一区| 欧美xxxx性猛交bbbb| 真人做人爱边吃奶动态| 精华霜和精华液先用哪个| 久久九九热精品免费| 欧美+日韩+精品| 国产私拍福利视频在线观看| 国产成人av教育| 韩国av一区二区三区四区| 最好的美女福利视频网| 午夜福利免费观看在线| 免费观看人在逋| 亚洲av不卡在线观看| 欧美最新免费一区二区三区 | 欧美日韩中文字幕国产精品一区二区三区| bbb黄色大片| 婷婷精品国产亚洲av| 国产私拍福利视频在线观看| 简卡轻食公司| 在线观看舔阴道视频| 午夜a级毛片| 1024手机看黄色片| 久久久久国产精品人妻aⅴ院| 国产一区二区在线av高清观看| 两人在一起打扑克的视频| 精品福利观看| 国产伦精品一区二区三区视频9| 岛国在线免费视频观看| 精品午夜福利视频在线观看一区| 搞女人的毛片| 成人三级黄色视频| 69人妻影院| 男人狂女人下面高潮的视频| 内射极品少妇av片p| 欧美日本亚洲视频在线播放| 九九热线精品视视频播放| 日韩 亚洲 欧美在线| 韩国av一区二区三区四区| 日本五十路高清| 色综合站精品国产| 中文字幕精品亚洲无线码一区| 脱女人内裤的视频| 99视频精品全部免费 在线| 美女大奶头视频| 少妇裸体淫交视频免费看高清| x7x7x7水蜜桃| 91在线观看av| 亚洲乱码一区二区免费版| 欧美日韩乱码在线| 国产久久久一区二区三区| 亚洲av二区三区四区| 一a级毛片在线观看| 欧美黑人巨大hd| 日本三级黄在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲熟妇熟女久久| 一进一出抽搐动态| 午夜日韩欧美国产| 国产野战对白在线观看| 免费大片18禁| 麻豆av噜噜一区二区三区| 久久精品国产亚洲av涩爱 | 色精品久久人妻99蜜桃| 直男gayav资源| 久久香蕉精品热| 亚洲自偷自拍三级| 精品一区二区三区视频在线| 啦啦啦韩国在线观看视频| 久久欧美精品欧美久久欧美| 国产精品女同一区二区软件 | 99热6这里只有精品| 欧美+日韩+精品| 最新中文字幕久久久久| 国模一区二区三区四区视频| 老女人水多毛片| 12—13女人毛片做爰片一| 丰满乱子伦码专区| 欧美黄色淫秽网站| 欧美又色又爽又黄视频| 中国美女看黄片| 中文字幕av在线有码专区| 99国产极品粉嫩在线观看| 亚洲av一区综合| 亚洲精品乱码久久久v下载方式| 精品一区二区三区人妻视频| 女生性感内裤真人,穿戴方法视频| 无人区码免费观看不卡| 一本综合久久免费| 国产蜜桃级精品一区二区三区| 午夜福利视频1000在线观看| 美女被艹到高潮喷水动态| 精品人妻视频免费看| 99久久无色码亚洲精品果冻| 国产大屁股一区二区在线视频| 国产爱豆传媒在线观看| 欧美性猛交╳xxx乱大交人| 精品99又大又爽又粗少妇毛片 | 伦理电影大哥的女人| 宅男免费午夜| 国产中年淑女户外野战色| 欧美日韩瑟瑟在线播放| 一进一出好大好爽视频| 亚洲欧美日韩高清专用| 亚洲,欧美,日韩| 亚洲综合色惰| 国产美女午夜福利| 免费看光身美女| 欧美黄色片欧美黄色片| 18禁黄网站禁片午夜丰满| 国产免费男女视频| 久久婷婷人人爽人人干人人爱| 一区二区三区免费毛片| 窝窝影院91人妻| 日韩欧美精品v在线| a级毛片免费高清观看在线播放| 99热这里只有精品一区| 色哟哟哟哟哟哟| 国内毛片毛片毛片毛片毛片| 亚洲人成网站在线播| 999久久久精品免费观看国产| 综合色av麻豆| 99热这里只有精品一区| 中文字幕av成人在线电影| 亚洲激情在线av| 国产亚洲av嫩草精品影院| 天堂av国产一区二区熟女人妻| eeuss影院久久| 99精品久久久久人妻精品| 久久久久国产精品人妻aⅴ院| www.999成人在线观看| 欧美3d第一页| 日韩中文字幕欧美一区二区| 舔av片在线| 级片在线观看| av天堂在线播放| 十八禁人妻一区二区| 久久久久久久久久成人| 亚洲自偷自拍三级| 一级av片app| 亚洲一区二区三区色噜噜| 性色av乱码一区二区三区2| 国产色爽女视频免费观看| 人人妻人人澡欧美一区二区| 丝袜美腿在线中文| 在线国产一区二区在线| 日韩欧美一区二区三区在线观看| 热99re8久久精品国产| 少妇高潮的动态图| 亚洲国产欧美人成| 亚洲美女搞黄在线观看 | 麻豆国产av国片精品| 757午夜福利合集在线观看| 久久人人爽人人爽人人片va | 亚洲成av人片在线播放无| 丰满的人妻完整版| 久久99热6这里只有精品| 亚洲中文字幕日韩| 观看美女的网站| 亚洲三级黄色毛片| 亚洲片人在线观看| 中文字幕精品亚洲无线码一区| 在线观看66精品国产| 欧美性猛交黑人性爽| 一本精品99久久精品77| 国产av在哪里看| 午夜精品在线福利| 欧美激情久久久久久爽电影| 久久精品综合一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 99国产极品粉嫩在线观看| 亚洲欧美清纯卡通| 亚洲专区国产一区二区| 18禁在线播放成人免费| 亚洲av成人精品一区久久| 97人妻精品一区二区三区麻豆| av女优亚洲男人天堂| 国语自产精品视频在线第100页| 嫩草影院精品99| 99热这里只有是精品50| 欧美在线黄色| 99热6这里只有精品| 永久网站在线| 免费人成视频x8x8入口观看| 在线播放无遮挡| 十八禁网站免费在线| 日本撒尿小便嘘嘘汇集6| 欧美日本视频| 国产蜜桃级精品一区二区三区| 最新在线观看一区二区三区| 国产一级毛片七仙女欲春2| av在线天堂中文字幕| 97超级碰碰碰精品色视频在线观看| 一个人看视频在线观看www免费| 国产一区二区在线观看日韩| 日韩av在线大香蕉| 国产亚洲精品综合一区在线观看| 天堂动漫精品| 亚洲第一电影网av| 免费看a级黄色片| 国产真实乱freesex| 欧美三级亚洲精品| 女人十人毛片免费观看3o分钟| 看片在线看免费视频| 久久中文看片网| 亚洲精品456在线播放app | 日本成人三级电影网站| 制服丝袜大香蕉在线| 天堂av国产一区二区熟女人妻| 日本精品一区二区三区蜜桃| 精品人妻一区二区三区麻豆 | 在现免费观看毛片| 久久国产精品影院| 男女之事视频高清在线观看| 波多野结衣巨乳人妻| 免费高清视频大片| 国产大屁股一区二区在线视频| av欧美777| 91麻豆av在线| 午夜福利18| 久久99热6这里只有精品| 91字幕亚洲| 国产真实乱freesex| 亚洲 欧美 日韩 在线 免费| 亚洲国产精品成人综合色| 欧美乱色亚洲激情| 噜噜噜噜噜久久久久久91| 日韩大尺度精品在线看网址| 宅男免费午夜| www.熟女人妻精品国产| 日本 av在线| 国产v大片淫在线免费观看| 午夜免费成人在线视频| 一进一出抽搐动态| 国产高潮美女av| 1024手机看黄色片| 国产午夜福利久久久久久| 99riav亚洲国产免费| .国产精品久久| 在线观看午夜福利视频| 国产精品1区2区在线观看.| 亚洲欧美日韩卡通动漫| 一个人看视频在线观看www免费| 欧美性感艳星| 国产熟女xx| 欧美性猛交╳xxx乱大交人| 欧美高清成人免费视频www| 久久久久久大精品| 少妇的逼水好多| 欧美日韩中文字幕国产精品一区二区三区| 窝窝影院91人妻| 九色成人免费人妻av| 免费在线观看日本一区| 一二三四社区在线视频社区8| 国产探花在线观看一区二区| www.999成人在线观看| 欧美最黄视频在线播放免费| 人人妻,人人澡人人爽秒播| 亚洲成人久久爱视频| 国产精品综合久久久久久久免费| 麻豆国产av国片精品| 国产高清视频在线播放一区| 一级作爱视频免费观看| 欧美一区二区精品小视频在线| 能在线免费观看的黄片| 国产三级中文精品| 国产视频内射| 村上凉子中文字幕在线| 麻豆av噜噜一区二区三区| 最好的美女福利视频网| 久久99热6这里只有精品| 午夜福利18| 亚洲第一欧美日韩一区二区三区| 成人三级黄色视频| 精品日产1卡2卡| 久久这里只有精品中国| 欧美高清性xxxxhd video| 黄色女人牲交| 九色国产91popny在线| 天堂动漫精品| 国产av麻豆久久久久久久| 欧美日韩中文字幕国产精品一区二区三区| 麻豆一二三区av精品| 国产一区二区三区在线臀色熟女| 国产午夜精品论理片|