• <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手机在线免费观看| 麻豆成人午夜福利视频| 免费人成在线观看视频色| 看十八女毛片水多多多| 亚洲欧美精品自产自拍| 国产永久视频网站| 国产免费一级a男人的天堂| 免费大片黄手机在线观看| 特级一级黄色大片| 人妻制服诱惑在线中文字幕| 精品国产一区二区三区久久久樱花 | 日产精品乱码卡一卡2卡三| 五月伊人婷婷丁香| 久久久久久久久大av| 夫妻午夜视频| 麻豆乱淫一区二区| 久久精品国产亚洲av涩爱| 精品人妻偷拍中文字幕| 国产亚洲午夜精品一区二区久久 | 男人舔奶头视频| 丰满人妻一区二区三区视频av| 亚洲国产av新网站| 欧美日韩国产mv在线观看视频 | 不卡视频在线观看欧美| 久久精品国产自在天天线| 午夜免费鲁丝| 99热这里只有精品一区| av天堂中文字幕网| 国产欧美另类精品又又久久亚洲欧美| 成年人午夜在线观看视频| 成人美女网站在线观看视频| 国内少妇人妻偷人精品xxx网站| 国产精品人妻久久久影院| 18禁裸乳无遮挡免费网站照片| 亚洲图色成人| 高清午夜精品一区二区三区| 亚洲一级一片aⅴ在线观看| 亚洲高清免费不卡视频| 亚洲人成网站在线播| 看免费成人av毛片| 草草在线视频免费看| 色播亚洲综合网| 狠狠精品人妻久久久久久综合| 内射极品少妇av片p| 亚洲精品久久久久久婷婷小说| 欧美少妇被猛烈插入视频| 在线免费十八禁| 日韩视频在线欧美| 国产免费一级a男人的天堂| 啦啦啦在线观看免费高清www| 一本一本综合久久| 成人一区二区视频在线观看| 交换朋友夫妻互换小说| 26uuu在线亚洲综合色| 亚洲国产欧美人成| 2021少妇久久久久久久久久久| 国产精品一二三区在线看| 日韩成人伦理影院| 久久精品久久久久久噜噜老黄| 女人十人毛片免费观看3o分钟| 晚上一个人看的免费电影| 国产探花极品一区二区| 在线观看av片永久免费下载| 2021天堂中文幕一二区在线观| 成人综合一区亚洲| 夜夜爽夜夜爽视频| 亚洲性久久影院| 精品国产露脸久久av麻豆| 国产女主播在线喷水免费视频网站| 一区二区av电影网| 一级毛片电影观看| 插阴视频在线观看视频| 国产精品久久久久久精品古装| 国产色婷婷99| 白带黄色成豆腐渣| 免费高清在线观看视频在线观看| 亚洲av免费在线观看| 欧美性感艳星| 亚洲av免费高清在线观看| 亚洲人成网站高清观看| 男人舔奶头视频| 国产精品女同一区二区软件| 成年av动漫网址| 亚洲av二区三区四区| 久久久久久久久大av| 午夜福利视频1000在线观看| 国产午夜精品一二区理论片| 街头女战士在线观看网站| 最近中文字幕高清免费大全6| 超碰97精品在线观看| 日日摸夜夜添夜夜添av毛片| 少妇人妻 视频| 一个人观看的视频www高清免费观看| 丝瓜视频免费看黄片| 欧美最新免费一区二区三区| 国产成人精品婷婷| 三级国产精品片| 亚洲在久久综合| 久久精品久久久久久噜噜老黄| 国产免费一区二区三区四区乱码| 亚洲丝袜综合中文字幕| 久久国产乱子免费精品| 能在线免费看毛片的网站| 最近最新中文字幕大全电影3| 国产精品一区二区性色av| 色哟哟·www| 麻豆精品久久久久久蜜桃| 国产成人一区二区在线| 日韩不卡一区二区三区视频在线| 毛片女人毛片| 色5月婷婷丁香| av免费观看日本| 好男人视频免费观看在线| 国产高清三级在线| 国产黄频视频在线观看| 国产精品久久久久久av不卡| tube8黄色片| 国产高潮美女av| 我要看日韩黄色一级片| 国产乱人偷精品视频| 国产免费福利视频在线观看| 美女xxoo啪啪120秒动态图| 极品教师在线视频| 亚洲第一区二区三区不卡| 国产一区二区亚洲精品在线观看| 亚洲人成网站在线播| 夜夜看夜夜爽夜夜摸| 一级a做视频免费观看| 人妻系列 视频| 国产毛片a区久久久久| 我的女老师完整版在线观看| 一二三四中文在线观看免费高清| 久久午夜福利片| 黄片wwwwww| 国产成年人精品一区二区| 18禁在线播放成人免费| 老司机影院成人| 别揉我奶头 嗯啊视频| 亚洲国产精品999| 亚洲激情五月婷婷啪啪| 国产亚洲午夜精品一区二区久久 | 国产日韩欧美亚洲二区| 国产成人freesex在线| 男女边摸边吃奶| 大片电影免费在线观看免费| av国产久精品久网站免费入址| 国产一级毛片在线| 久热这里只有精品99| 欧美xxⅹ黑人| 免费av毛片视频| 国产片特级美女逼逼视频| 国产伦理片在线播放av一区| 亚洲国产欧美在线一区| 天堂俺去俺来也www色官网| 久久久久国产网址| 色吧在线观看| 久久精品久久久久久久性| 伦理电影大哥的女人| 国产精品国产三级国产av玫瑰| 97精品久久久久久久久久精品| 高清视频免费观看一区二区| 日韩三级伦理在线观看| 女的被弄到高潮叫床怎么办| 男女边摸边吃奶| 免费观看无遮挡的男女| 久久久欧美国产精品| 免费播放大片免费观看视频在线观看| 国产淫片久久久久久久久| 乱码一卡2卡4卡精品| 热re99久久精品国产66热6| 久久午夜福利片| 99热6这里只有精品| 男人添女人高潮全过程视频| 国产熟女欧美一区二区| 国内揄拍国产精品人妻在线| 久久综合国产亚洲精品| 国产中年淑女户外野战色| 国产淫语在线视频| 国产老妇女一区| 久久久久精品久久久久真实原创| 综合色丁香网| 国产成人aa在线观看| 成人一区二区视频在线观看| 你懂的网址亚洲精品在线观看| 99久久精品热视频| 亚洲国产精品成人久久小说| 成人免费观看视频高清| 亚洲最大成人av| 亚洲精品一区蜜桃| 男女边摸边吃奶| 欧美成人a在线观看| 波多野结衣巨乳人妻| 亚洲av成人精品一二三区| 亚洲成人久久爱视频| 国产精品.久久久| 丰满人妻一区二区三区视频av| 99热这里只有是精品在线观看| 美女脱内裤让男人舔精品视频| 搡老乐熟女国产| 又爽又黄a免费视频| 久久精品人妻少妇| 国产亚洲av嫩草精品影院| 欧美变态另类bdsm刘玥| 18禁裸乳无遮挡动漫免费视频 | 男的添女的下面高潮视频| 久久6这里有精品| 欧美97在线视频| 精品人妻视频免费看| 性色avwww在线观看| 欧美高清成人免费视频www| 久久久成人免费电影| 99热6这里只有精品| 欧美最新免费一区二区三区| 国产高清国产精品国产三级 | 99re6热这里在线精品视频| 各种免费的搞黄视频| 麻豆乱淫一区二区| 在线看a的网站| 卡戴珊不雅视频在线播放| 精品久久久久久久人妻蜜臀av| 日韩制服骚丝袜av| 国产精品伦人一区二区| 久久久亚洲精品成人影院| 91精品伊人久久大香线蕉| 精品视频人人做人人爽| 国产成人91sexporn| 水蜜桃什么品种好| 国产有黄有色有爽视频| 69人妻影院| 联通29元200g的流量卡| 精品久久久久久久久亚洲| 黄色视频在线播放观看不卡| 一级毛片黄色毛片免费观看视频| 69av精品久久久久久| 国产伦精品一区二区三区四那| 亚洲av成人精品一二三区| 99热全是精品| 国产高潮美女av| 亚洲婷婷狠狠爱综合网| kizo精华| 精品久久久久久久久av| 久久精品国产亚洲av天美| 中文字幕av成人在线电影| 精品久久久精品久久久| 日韩,欧美,国产一区二区三区| 欧美精品国产亚洲| 一级av片app| 成人毛片a级毛片在线播放| 卡戴珊不雅视频在线播放| 在线播放无遮挡| 各种免费的搞黄视频| 久热久热在线精品观看| 日产精品乱码卡一卡2卡三| 日本一二三区视频观看| 色综合色国产| av在线播放精品| 成人午夜精彩视频在线观看| 777米奇影视久久| 香蕉精品网在线| 午夜福利网站1000一区二区三区| 毛片女人毛片| 你懂的网址亚洲精品在线观看| 久久精品久久精品一区二区三区| 国产精品不卡视频一区二区| av一本久久久久| 国产精品麻豆人妻色哟哟久久| 男女下面进入的视频免费午夜| 久久精品国产亚洲网站| 永久网站在线| 插阴视频在线观看视频| 九草在线视频观看| 亚洲av电影在线观看一区二区三区 | 日产精品乱码卡一卡2卡三| 色婷婷久久久亚洲欧美| 亚洲第一区二区三区不卡| 毛片一级片免费看久久久久| 免费大片黄手机在线观看| 激情 狠狠 欧美| 中文在线观看免费www的网站| 国产免费一级a男人的天堂| 高清日韩中文字幕在线| 亚洲av中文av极速乱| 在线观看三级黄色| 精品人妻熟女av久视频| 欧美变态另类bdsm刘玥| 22中文网久久字幕| 一级毛片电影观看| 欧美丝袜亚洲另类| 日日摸夜夜添夜夜添av毛片| 日本午夜av视频| 欧美xxxx黑人xx丫x性爽| 国产大屁股一区二区在线视频| 1000部很黄的大片| 激情五月婷婷亚洲| 在线亚洲精品国产二区图片欧美 | 亚洲av中文av极速乱| 久久精品久久精品一区二区三区| 国产 一区 欧美 日韩| 中文欧美无线码| 大香蕉97超碰在线| 免费av观看视频| 久久ye,这里只有精品| 国产在线一区二区三区精| 久久国产乱子免费精品| 美女视频免费永久观看网站| 亚洲成人中文字幕在线播放| 国产精品伦人一区二区| 亚洲成人中文字幕在线播放| 色视频www国产| 亚洲第一区二区三区不卡| videossex国产| 亚洲国产日韩一区二区| 深爱激情五月婷婷| 亚洲无线观看免费| 18禁在线播放成人免费| 国产精品人妻久久久久久| 国产黄a三级三级三级人| 真实男女啪啪啪动态图| 国产精品人妻久久久影院| 国产综合懂色| 又大又黄又爽视频免费| 在线 av 中文字幕| 国产免费又黄又爽又色| 欧美极品一区二区三区四区| 国产伦理片在线播放av一区| 国产爱豆传媒在线观看| 18禁在线播放成人免费| 秋霞在线观看毛片| 国产高清不卡午夜福利| 亚洲自偷自拍三级| 97精品久久久久久久久久精品| 国产成人精品福利久久| 午夜激情久久久久久久| 女人被狂操c到高潮| 最近中文字幕2019免费版| 久久久久久久午夜电影| 成人亚洲精品av一区二区| 边亲边吃奶的免费视频| 水蜜桃什么品种好| 国产黄色免费在线视频| 久久精品国产亚洲网站| 男人舔奶头视频| 欧美xxxx性猛交bbbb| 我的老师免费观看完整版| 最近手机中文字幕大全| 免费观看在线日韩| 日本欧美国产在线视频| 男的添女的下面高潮视频| 亚洲国产色片| 久久热精品热| 国产精品国产三级国产专区5o| 日日啪夜夜撸| 国产精品国产三级国产专区5o| 神马国产精品三级电影在线观看| 韩国av在线不卡| 日韩av在线免费看完整版不卡| 人妻少妇偷人精品九色| freevideosex欧美| 干丝袜人妻中文字幕| 国产午夜精品久久久久久一区二区三区| 日韩一区二区三区影片| 91aial.com中文字幕在线观看| 三级国产精品片| 国产亚洲5aaaaa淫片| 久久精品国产亚洲网站| 久久国产乱子免费精品| 熟妇人妻不卡中文字幕| 一级毛片 在线播放| 乱码一卡2卡4卡精品| 99久久中文字幕三级久久日本| videos熟女内射| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 高清午夜精品一区二区三区| 五月天丁香电影| 91久久精品电影网| 韩国av在线不卡| 久久99蜜桃精品久久| 国产乱来视频区| av播播在线观看一区| 干丝袜人妻中文字幕| 免费观看a级毛片全部| 久久久久久久大尺度免费视频| 在线观看一区二区三区| 中文在线观看免费www的网站| 我的女老师完整版在线观看| 亚洲三级黄色毛片| 亚洲最大成人av| 国产毛片在线视频| 色视频在线一区二区三区| 亚洲国产日韩一区二区| videossex国产| 亚洲国产高清在线一区二区三| 国内精品宾馆在线| 91午夜精品亚洲一区二区三区| 午夜激情福利司机影院| 精品人妻视频免费看| 国产欧美日韩一区二区三区在线 | 男人添女人高潮全过程视频| 综合色av麻豆| 免费看光身美女| 欧美zozozo另类| 成人免费观看视频高清| 最近手机中文字幕大全| 久久热精品热| 亚洲成人一二三区av| 日韩伦理黄色片| 亚洲美女搞黄在线观看| 一级毛片 在线播放| 国产免费又黄又爽又色| 嫩草影院精品99| 天天一区二区日本电影三级| 欧美日韩国产mv在线观看视频 | 日韩av免费高清视频| 国产成人午夜福利电影在线观看| 亚洲精品久久午夜乱码| 黄色日韩在线| 一本一本综合久久| 免费高清在线观看视频在线观看| 狠狠精品人妻久久久久久综合| 亚洲综合色惰| 久久这里有精品视频免费| 在线观看国产h片| 一个人看视频在线观看www免费| 少妇猛男粗大的猛烈进出视频 | 久久精品国产亚洲网站| 久久精品国产自在天天线| 97在线人人人人妻| 大码成人一级视频| 赤兔流量卡办理| 视频中文字幕在线观看| 亚洲色图av天堂| 免费观看的影片在线观看| 精品人妻偷拍中文字幕| 干丝袜人妻中文字幕| 日韩不卡一区二区三区视频在线| 成人美女网站在线观看视频| 舔av片在线| 黄色一级大片看看| 毛片女人毛片| 欧美亚洲 丝袜 人妻 在线| 亚洲精品456在线播放app| 亚洲美女视频黄频| 免费少妇av软件| 亚洲最大成人av| 免费大片黄手机在线观看| 日韩欧美精品免费久久| 亚洲国产精品成人久久小说| 日日撸夜夜添| 麻豆精品久久久久久蜜桃| 在线观看av片永久免费下载| 国产亚洲av嫩草精品影院| 成人亚洲精品一区在线观看 | 中文欧美无线码| 亚洲精品成人av观看孕妇| 日韩人妻高清精品专区| 丰满少妇做爰视频| 国产精品一区二区在线观看99| 午夜亚洲福利在线播放| kizo精华| 激情 狠狠 欧美| 看黄色毛片网站| 欧美少妇被猛烈插入视频| 亚洲丝袜综合中文字幕| 国产黄a三级三级三级人| 日韩视频在线欧美| 亚洲国产欧美在线一区| 日本一本二区三区精品| 最近最新中文字幕免费大全7| 精品一区二区三卡| 日韩成人av中文字幕在线观看| 国产综合懂色| 综合色av麻豆| 三级经典国产精品| 简卡轻食公司| 少妇被粗大猛烈的视频| 国产免费又黄又爽又色| 综合色丁香网| 欧美三级亚洲精品| 色视频www国产| 国产精品久久久久久精品电影| 成人欧美大片| 天美传媒精品一区二区| 99视频精品全部免费 在线| 国产精品一区www在线观看| 亚洲欧美精品专区久久| 国产大屁股一区二区在线视频| 亚洲精品国产成人久久av| 国产精品一区二区三区四区免费观看| 五月天丁香电影| 国产一区二区三区av在线| 精品亚洲乱码少妇综合久久| 麻豆国产97在线/欧美| 久久精品国产亚洲av涩爱| 一本久久精品| 亚洲av欧美aⅴ国产| 日韩 亚洲 欧美在线| 男的添女的下面高潮视频| 久久久久久伊人网av| 丰满少妇做爰视频| 国产精品一区二区三区四区免费观看| 精品少妇久久久久久888优播| 亚洲国产精品999| 国产高清国产精品国产三级 | 中文字幕久久专区| 亚州av有码| 在线观看一区二区三区激情| 亚洲伊人久久精品综合| 高清在线视频一区二区三区| 九九久久精品国产亚洲av麻豆| 亚洲国产精品国产精品| 亚洲伊人久久精品综合| 亚洲国产精品国产精品| 我的女老师完整版在线观看| 精品人妻一区二区三区麻豆| 亚洲伊人久久精品综合| 老女人水多毛片| 丰满人妻一区二区三区视频av| 国产淫语在线视频| 亚洲av电影在线观看一区二区三区 | 国产精品秋霞免费鲁丝片| av专区在线播放| 亚洲欧美日韩卡通动漫| 国产爱豆传媒在线观看| 亚洲美女搞黄在线观看| 亚洲欧洲日产国产| 国产白丝娇喘喷水9色精品| 亚洲国产精品成人久久小说| 99热国产这里只有精品6| 可以在线观看毛片的网站| 啦啦啦中文免费视频观看日本| 国产免费一级a男人的天堂| 欧美日韩亚洲高清精品| 国产男女内射视频| 大香蕉久久网| 最近的中文字幕免费完整| 最新中文字幕久久久久| 亚洲av中文字字幕乱码综合| 精品久久久久久久末码| 欧美亚洲 丝袜 人妻 在线| 中文在线观看免费www的网站| 一级二级三级毛片免费看| 最近中文字幕2019免费版| 天堂俺去俺来也www色官网| 97超视频在线观看视频| 国产成人freesex在线| 在线免费十八禁| 日韩不卡一区二区三区视频在线| 人人妻人人看人人澡| 一级毛片我不卡| 亚洲成人精品中文字幕电影| 国产免费又黄又爽又色| 亚洲国产精品成人综合色| 精品人妻偷拍中文字幕| 久久精品国产a三级三级三级| 国产高清国产精品国产三级 | 亚洲国产欧美在线一区| 国产毛片a区久久久久| 51国产日韩欧美| 精品一区二区免费观看| 国产精品不卡视频一区二区| 大香蕉97超碰在线| 精品国产乱码久久久久久小说| 精品久久久久久久人妻蜜臀av| 丰满少妇做爰视频| 国产真实伦视频高清在线观看| 青春草视频在线免费观看| 国产精品久久久久久av不卡| 免费高清在线观看视频在线观看| 国产一区有黄有色的免费视频| 午夜免费鲁丝| 亚洲av欧美aⅴ国产| 伊人久久国产一区二区| 视频中文字幕在线观看| 国内少妇人妻偷人精品xxx网站| 亚洲国产精品国产精品| 欧美日韩视频高清一区二区三区二| 身体一侧抽搐| 日本三级黄在线观看| 搡老乐熟女国产| 国产成人精品久久久久久| 国产爽快片一区二区三区| 嫩草影院精品99| 大话2 男鬼变身卡| 亚洲精品乱码久久久v下载方式| 听说在线观看完整版免费高清| 寂寞人妻少妇视频99o| 欧美成人午夜免费资源| 久久综合国产亚洲精品| 欧美精品一区二区大全| 新久久久久国产一级毛片| 青青草视频在线视频观看| av在线蜜桃| 亚洲美女搞黄在线观看| 久久综合国产亚洲精品| av网站免费在线观看视频| 极品少妇高潮喷水抽搐| 视频中文字幕在线观看| 国产乱人视频| 亚洲人成网站高清观看| 国产精品久久久久久av不卡| 国产精品福利在线免费观看| 一级毛片aaaaaa免费看小| 成人无遮挡网站| 国产 一区精品| 亚洲av中文字字幕乱码综合| 午夜福利视频1000在线观看|