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

    油棕EgGRF轉(zhuǎn)錄因子家族的全基因組鑒定與表達分析

    2021-08-26 08:58:41金龍飛尹欣幸曹紅星
    熱帶作物學(xué)報 2021年7期
    關(guān)鍵詞:油棕基因表達生物信息學(xué)

    金龍飛 尹欣幸 曹紅星

    摘? 要:生長調(diào)控因子(Growth regulating factor,GRF)是植物中重要的轉(zhuǎn)錄因子,參與植物生長、發(fā)育和逆境脅迫響應(yīng)等多種生物學(xué)過程。本研究從油棕(Elaeis guineensis)基因組中鑒定出15個EgGRF基因家族成員,并對其理化性質(zhì)、染色體定位、基因結(jié)構(gòu)、保守功能域、進化關(guān)系、啟動子順式作用元件、組織表達模式和果肉不同時期表達模式進行分析。結(jié)果表明:EgGRF基因家族成員編碼的肽鏈平均為409個氨基酸,分子量為27.62~65.51 kDa,等電點為6.24~9.38,蛋白不穩(wěn)定指數(shù)為47.24~68.37,脂溶指數(shù)為47.52~67.69,總平均親水性為–0.945~–0.400。EgGRF基因家族成員含有3~5個外顯子,均含有特征結(jié)構(gòu)域QLQ(Gln、Leu、Gln)和WRC(Trp、Arg、Cys),基于系統(tǒng)進化關(guān)系將EgGRF家族分為5個亞族,油棕EgGRF的親緣性與擬南芥較近。啟動子上鑒定出大量植物激素響應(yīng)、逆境脅迫響應(yīng)、光響應(yīng)和分生組織特異表達順式作用元件。不同組織的轉(zhuǎn)錄組數(shù)據(jù)分析結(jié)果表明,EgGRF基因家族在莖尖和花中表達量較高,8個EgGRF在果肉的不同時期特異表達。本研究為進一步探索EgGRF調(diào)控油棕生長發(fā)育過程的機制奠定基礎(chǔ)。

    關(guān)鍵詞:EgGRF;油棕;生物信息學(xué);基因表達

    中圖分類號:S565.9????? 文獻標識碼:A

    Genome-wide Identification and Expression of EgGRF Gene Family in Oil Palm

    JIN Longfei, YIN Xinxing, CAO Hongxing*

    Coconut Research Institute, Chinese Academy of Tropical Agricultural Sciences / Hainan Key Laboratory of Tropical Oil Crops Biology, Wenchang, Hainan 571339, China

    Abstract: Growth regulating factor (GRF) is an important transcriptional factor in plants which is involved in various biological processes, including regulating plant growth, development, and stress response. In this study, 15 EgGRF transcription factor family members were identified from the oil palm genome using bioinformatics methods, and the physicochemical properties, chromosomal location, gene structure, conserved functional domains, evolutionary relationships, promoter cis-acting elements, tissue expression patterns and mesocarp expression patterns at different stages were analyzed. The results showed that the EgGRF family members encoded an average of 409 aa in the peptide. The molecular weight was between 27.62 and 65.51 kDa, the isoelectric point was between 6.24 and 9.38, the instability index was between 47.24 and 68.37, the aliphatic index was between 47.52 and 67.69, and the gravy was between –0.945 and –0.400. Each EgGRF members contained 3-5 exons, and contained QLQ (Gln, Leu, Gln) and WRC (Trp, Arg, Cys) domain. EgGRF family was divided into 5 groups based on phylogenetic relationships, and was closely related to the Arabidopsis. Analysis of transcriptome data from different tissues showed that the EgGRF family was significantly high expressed in the shoot and flower. Eight EgGRFs were specifically expressed in different stages of mesocarp. The results would provide useful information and resource for the further investigation of EgGRF on oil palm growth and development.

    Keywords: EgGRF; oil palm; bioinformatics; gene expression

    DOI: 10.3969/j.issn.1000-2561.2021.07.002

    轉(zhuǎn)錄因子是一類重要的調(diào)控基因,參與植物生長、發(fā)育、代謝、繁殖、分化、逆境響應(yīng)等多種生物學(xué)過程[1-3]。目前在植物中已經(jīng)發(fā)現(xiàn)了60多個轉(zhuǎn)錄因子家族,其中生長調(diào)控因子(Growth regulating factor,GRF)主要在植物生長發(fā)育中起重要的調(diào)控作用[4-5]。GRF轉(zhuǎn)錄因子的N端有2個保守結(jié)構(gòu)域,分別是QLQ(Gln、Leu、Gln)和WRC(Trp、Arg、Cys)保守結(jié)構(gòu)域[6]。QLQ結(jié)構(gòu)域可通過與GRF互作因子(GRF-interacting factor1,GIF1)結(jié)合,形成轉(zhuǎn)錄共激活因子,激活下游的靶基因表達;WRC結(jié)構(gòu)域包含一個核定位信號NLS區(qū)域和一個與DNA結(jié)合的鋅指基序,均與GRF轉(zhuǎn)錄因子生物學(xué)功能密切相關(guān)[7]。植物中第一個編碼GRF轉(zhuǎn)錄因子的基因OsGRF1是在水稻中發(fā)現(xiàn)的,該基因通過調(diào)控赤霉素代謝誘導(dǎo)莖伸長[8]。OsGRF4通過調(diào)控2個細胞分裂素脫氫酶前體基因(OsCKX5和OsCKX1)的表達,增加細胞分裂素的含量,進而促進稻穗的伸長[9]。擬南芥AtGRF5通過調(diào)控細胞增殖,刺激葉綠體分裂,增加葉綠素含量,促進葉片生長[10]。AtGRF1、AtGRF4、AtGRF8、AtGRF9通過調(diào)控PTL基因控

    制體胚發(fā)生過程中內(nèi)源生長素的合成[11]。OsGRF4還能與miR396c、OsGIF1互作調(diào)控水稻籽粒大小[12],對提高水稻產(chǎn)量具有重要意義。在油菜中的研究也發(fā)現(xiàn),超表達BnGRF2能夠調(diào)控細胞分裂、葉綠素含量和光合作用效率,進而增加油菜的產(chǎn)油量[13]。隨著越來越多植物基因組測序完成,GRF基因家族在擬南芥[14]、水稻[15]、番茄[16]、藜麥[17]、大豆[18]、茶樹[19]、棉花[20]、龍眼[21]等多個物種中被鑒定出來,而在油棕中的研究未見報道。

    油棕(Elaeis guineensis)是世界上產(chǎn)油效率最高的油料作物,油量高達4.27 t/hm2,是花生的7~8倍、大豆的9~10倍[22]。油棕果壓榨的棕櫚油,廣泛應(yīng)用于食品加工、日用化工、機械潤滑和生物柴油等領(lǐng)域[23]。近年來市場對棕櫚油的需求逐年增加,油棕種植業(yè)的快速發(fā)展毀壞了大量的熱帶雨林。因此,提高油棕單產(chǎn)是維持熱帶雨林保護和食用油短缺平衡的重要途徑,而挖掘控制油棕重要農(nóng)藝性狀的關(guān)鍵基因并采用分子育種加快高產(chǎn)油棕品種的選育是油棕育種研究的重要工

    作[24]。本研究基于油棕基因組數(shù)據(jù)庫挖掘和鑒定調(diào)控油棕生長的轉(zhuǎn)錄因子EgGRF,分析其理化性質(zhì)、染色體定位、基因結(jié)構(gòu)、保守結(jié)構(gòu)域、進化關(guān)系、啟動子順式作用元件,同時分析其在不同組織和果肉不同發(fā)育時期的表達模式,為進一步探索EgGRF調(diào)控油棕生長發(fā)育過程的機制奠定基礎(chǔ)。

    1? 材料與方法

    1.1? 油棕EgGRF基因家族的挖掘與鑒定

    油棕全基因組數(shù)據(jù)從NCBI數(shù)據(jù)庫(https://www.ncbi.nlm.nih.gov/ ID:2669)中下載,擬南芥AtGRF蛋白質(zhì)序列從TAIR數(shù)據(jù)庫(https://www.arabidopsis.org/)中下載,水稻OsGRF蛋白序列從NCBI數(shù)據(jù)庫中下載。以AtGRF和OsGRF蛋白質(zhì)序列作為參考序列,通過Blsatp在油棕基因組數(shù)據(jù)庫中進行比對。將獲得的蛋白序列提交到NCBI的保守結(jié)構(gòu)域數(shù)據(jù)庫,利用Batch Web CD-search Tool軟件[25]進行保守結(jié)構(gòu)域分析。利用在線工具ExPAsy(http://www.expasy.org)分析油棕EgGRF基因家族的蛋白質(zhì)分子量、等電點、蛋白不穩(wěn)定指數(shù)、脂溶指數(shù)和總平均親水性等理化性質(zhì)。

    1.2? 油棕EgGRF基因家族的基因結(jié)構(gòu)、保守結(jié)構(gòu)域、染色體定位及進化分析

    油棕EgGRF基因家族的基因結(jié)構(gòu)、保守結(jié)構(gòu)域和染色體定位均采用Tbtool軟件[26]進行可視化分析。采用ClustalW軟件[27]對油棕、擬南芥和水稻的GRF蛋白質(zhì)序列進行多重比對,采用MEGA 6.0軟件的Neighbor-joining法進行建樹,分析進化關(guān)系,設(shè)置校驗值Bootstrap為1000。

    1.3? 油棕EgGRF基因家族的啟動子順式作用元件分析

    從NCBI數(shù)據(jù)庫中下載油棕EgGRF基因CDS上游2000 bp的序列。利用Plant Care軟件[28]進行順式作用元件分析,采用Tbtool軟件進行可視化分析[26]。

    1. 4? 油棕EgGRF基因家族的表達分析

    油棕的根(SRR851071)、莖尖(SRR851103)、葉(SRR851096)、花(SRR851108)和花后15周(SRR190698)、17周(SRR190699)、21周(SRR190701)和23周(SRR190702)果肉的轉(zhuǎn)錄組數(shù)據(jù)從NCBI的SRA數(shù)據(jù)庫(https://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?)中下載。利用RPKM值對轉(zhuǎn)錄組數(shù)據(jù)進行歸一化處理,利用基迪奧云平臺在線工具(https://www.omicshare.com/tools/)繪制熱圖。

    2? 結(jié)果與分析

    2.1? 油棕EgGRF基因家族的全基因組鑒定和理化性質(zhì)分析

    通過Blastp比對分析獲得油棕EgGRF的候選蛋白質(zhì)序列,預(yù)測保守結(jié)構(gòu)域,去除不含QLQ和WRC保守結(jié)構(gòu)域的序列,共獲得15個EgGRF家族成員。根據(jù)其與擬南芥的同源性順序進行命名,依次命名為EgGRF1~EgGRF15。油棕EgGRF多肽鏈氨基酸數(shù)目為260~611,平均為409;蛋白質(zhì)分子量為27.62~65.51 kDa,平均為44.71 kDa;等電點為6.24~9.38,平均為7.86;蛋白不穩(wěn)定指數(shù)為47.24~68.37,平均為59.48,表明EgGRF家族成員編碼的均為不穩(wěn)定蛋白;脂溶系數(shù)為47.52~67.69,平均為59.23,脂溶性一般;總平均親水性為–0.945~–0.400,平均為–0.633,都為負值,表明EgGRF均為親水性蛋白(表1)。

    2.2? 油棕EgGRF基因家族成員在染色體上的定位

    對油棕EgGRF基因家族在染色體上位置進行分析(圖1),發(fā)現(xiàn)15個EgGRF基因家族成員分布 在7條染色體和5個Scaffolds上。EgGRF10在Chr3上,EgGRF2、EgGRF7在Chr4上,EgGRF5在Chr10上,EgGRF4在Chr11上,EgGRF9、EgGRF11在Chr13上,EgGRF6、EgGRF12在Chr15上,EgGRF14在Chr16,EgGRF1、EgGRF3、EgGRF8、EgGRF13、EgGRF15分別在5個Scaffolds上。EgGRF4和EgGRF7,EgGRF6和EgGRF9存在基因共線性現(xiàn)象。

    2.3? 油棕EgGRF的基因結(jié)構(gòu)、保守結(jié)構(gòu)域和進化關(guān)系

    15個油棕EgGRF基因家族具有3~5個外顯子,其中EgGRF4、EgGRF7有3個外顯子,EgGRF8、EgGRF13有5個外顯子,其余EgGRF成員均為4個外顯子(圖2)。15個EgGRF基因家族成員均含有GRF轉(zhuǎn)錄因子特有的WRC和QLQ保守結(jié)構(gòu)域(圖3)。為了分析油棕EgGRF基因家族的進化發(fā)育關(guān)系,將油棕EgGRF和AtGRF、OsGRF蛋白質(zhì)序列進行多重比對構(gòu)建進化樹。根據(jù)進化樹分支(圖4),可將15個EgGRF分成5個亞族,EgGRF5、EgGRF10、EgGRF11、EgGRF12、EgGRF13、EgGRF14、EgGRF15為Ⅰ亞族,EgGRF4、EgGRF7為Ⅱ亞族,EgGRF3為Ⅲ亞族,EgGRF2、EgGRF8為Ⅳ亞族,EgGRF1、EgGRF6、EgGRF9為V亞族。其中第Ⅲ和第Ⅳ亞族EgGRF僅與AtGRF聚在一起,表明油棕的EgGRF與擬南芥的親緣關(guān)系更近。

    2.4? 油棕EgGRF基因家族的啟動子順式作用元件的鑒定

    在油棕EgGRF基因家族的啟動子中鑒定了大量的順式作用元件(EgGRF8的啟動子中有大量的堿基未測序,故順式作用元件較少)(圖5),包括生長素響應(yīng)(10個)、脫落酸響應(yīng)(21個)、赤霉素響應(yīng)(19個)、茉莉酸甲酯響應(yīng)(21個)、水楊酸響應(yīng)(10個)、低溫響應(yīng)(3個)、干旱響應(yīng)(7個)、光響應(yīng)(53個)、防衛(wèi)和逆境響應(yīng)(4個)、分生組織特異表達(12個)等順式作用元件,表明EgGRF參與植物激素和環(huán)境刺激的應(yīng)答。

    2.5? 油棕EgGRF基因家族的表達分析

    油棕EgGRF基因家族在不同組織中表達量差異較大(圖6),其中EgGRF1、EgGRF2、EgGRF3、EgGRF5、EgGRF6、EgGRF7、EgGRF9、EgGRF10、EgGRF12、EgGRF13、EgGRF14在莖尖中高表達,EgGRF4、EgGRF11、EgGRF15在花中高表達;EgGRF基因家族在根和葉中的表達量均較低,僅EgGRF8在葉片中表達量較高,這表明GRF在油棕的莖尖和花發(fā)育過程中起重要調(diào)控作用。8個EgGRF基因在果肉發(fā)育的不同時期中特異表達不同(圖7),其中EgGRF2、EgGRF9、EgGRF14在花后15周的果肉中高表達,EgGRF6、EgGRF8在花后17周的果肉中高表達,EgGRF3在花后21周的果肉中高表達,EgGRF11、EgGRF15在花后23周的果肉中高表達,這表明不同的EgGRF分別在油棕發(fā)育的不同時期起調(diào)控作用。

    3? 討論

    隨著基因組測序技術(shù)的發(fā)展和生物信息學(xué)研

    究的不斷深入,大量植物完成了全基因組測序,這為植物基因家族的挖掘、鑒定和功能分析提供了有利條件。GRF轉(zhuǎn)錄因子是植物中特有的一類轉(zhuǎn)錄因子,在葉片生長、子葉生長、莖伸長、花發(fā)育、籽粒發(fā)育、胚胎發(fā)生等生長發(fā)育過程中起重要的調(diào)控作用[5, 8, 11, 14]。油棕作為世界上產(chǎn)油效率最高的油料作物,其GRF基因家族的研究尚未見報道,而在擬南芥和水稻等模式植物中研究較為深入,這些研究給油棕EgGRF的研究提供參考。

    本研究利用生物信息學(xué)分析,從油棕基因組中鑒定了15個EgGRF基因家族成員,比擬南芥(9個)和水稻(12個)成員數(shù)目更多,EgGRF4和EgGRF7、EgGRF6和EgGRF9存在基因共線性,這可能是油棕在進化過程中發(fā)生了基因串聯(lián)復(fù)制而導(dǎo)致EgGRF基因數(shù)目增加?;驈?fù)制事件能夠?qū)е轮参锘蚪M中形成大量的重復(fù)基因,重復(fù)基因的存在能夠促進基因新功能的進化,增強植物對環(huán)境變化的適應(yīng)性[29]。油棕EgGRF轉(zhuǎn)錄因子的蛋白長度、分子量、等電點、蛋白不穩(wěn)定指數(shù)、脂溶指數(shù)和總平均親水性等基本特性存在較大差異,與擬南芥和水稻中的研究相似?;虻慕Y(jié)構(gòu)分析發(fā)現(xiàn),絕大部分EgGRF基因具有4個外顯子,與擬南芥AtGRF家族的基因結(jié)構(gòu)分布相似,其中EgGRF8的第4個和第5個外顯子之間有個極短的內(nèi)含子,植物進化過程中轉(zhuǎn)座子的插入可能導(dǎo)致內(nèi)含子數(shù)目的增加[30],EgGRF8內(nèi)含子的產(chǎn)生可能是轉(zhuǎn)座子插入導(dǎo)致的。進化分析發(fā)現(xiàn)油棕EgGRF優(yōu)先與AtGRF聚在一起,表明油棕EgGRF與擬南芥的親緣關(guān)系更近。

    EgGRF主要在油棕的莖尖、花和果肉中表達,可能是由于EgGRF在這些組織或器官的生長發(fā)育中發(fā)揮著重要調(diào)控作用。研究發(fā)現(xiàn)GRF在幼嫩組織中的表達較高,在成熟的組織和器官中表達量低[31-33]。在本研究中EgGRF1、EgGRF2、EgGRF3、EgGRF5、EgGRF6、EgGRF7、EgGRF9、EgGRF10、EgGRF12、EgGRF13、EgGRF14在莖尖中的表達顯著高于其他組織。莖尖是油棕的頂端生長點,分生能力極強,EgGRF的高表達可能與油棕莖尖旺盛的分生能力相關(guān),同時EgGRF啟動子中也鑒定了大量分生組織特異表達響應(yīng)元件。在楊樹中的研究發(fā)現(xiàn)PtGRF在莖尖生長中起著重要的調(diào)控作用[34],在玉米中的研究發(fā)現(xiàn)ZmGRF與ZmGIF1基因互作調(diào)控莖尖的發(fā)育[35]。GRF還參與植物根的發(fā)育[36],而本研究發(fā)現(xiàn)在生長極為旺盛的根尖部分EgGRF整體表達都很低,這表明油棕根發(fā)育可能主要受其他轉(zhuǎn)錄因子的調(diào)控,例如NAC1、TIR1等。EgGRF4、EgGRF11、EgGRF15在花中的表達量顯著高于其他組織,表明EgGRF在油棕花發(fā)育中起重要的調(diào)控作用,在擬南芥和毛果楊中的研究也發(fā)現(xiàn)GRF對花的發(fā)育起重要的調(diào)控作用[37-38]。GRF在果實發(fā)育中起重要調(diào)控作用,在柑橘中的研究發(fā)現(xiàn)CsGRF3和CsGRF4在果實發(fā)育前期高表達,且受赤霉素的誘導(dǎo)[39]。花后15周,油棕完成細胞分裂和膨脹,進入成熟階段;花后17周,油棕果皮開始轉(zhuǎn)色;花后21周,油棕進入內(nèi)含物的快速積累期,乙烯含量達到峰值;花后23周,油棕進入完熟期,內(nèi)含物積累達到峰值[40]。為探索EgGRF在油棕果肉發(fā)育中的作用,選擇4個時期進行基因表達分析,結(jié)果發(fā)現(xiàn)EgGRF2、EgGRF6、EgGRF8、EgGRF9、EgGRF14在油棕果肉發(fā)育前期高表達,EgGRF3、EgGRF5、EgGRF11在油棕果肉發(fā)育后期高表達,這表明不同EgGRF在油棕果肉發(fā)育的不同時期都起調(diào)控作用,同時在EgGRF啟動子上也發(fā)現(xiàn)大量與果實發(fā)育密切相關(guān)的植物激素應(yīng)答元件,如生長素、赤霉素和脫落酸應(yīng)答元件。油棕果肉發(fā)育的不同階段受不同激素的調(diào)控,前期主要受生長素、赤霉素和細胞分裂素的調(diào)控,后期主要受乙烯和脫落酸的調(diào)控[40],而EgGRF的表達受赤霉素、脫落酸和水楊酸等植物激素的誘導(dǎo)[16],這表明EgGRF在激素控制果肉發(fā)育中起著重要的調(diào)控作用。果肉是油棕主要的產(chǎn)油器官,因此明確EgGRF在果肉發(fā)育中的調(diào)控機制對提高油棕的產(chǎn)油量具有重要的意義。在后續(xù)的研究中還需要利用轉(zhuǎn)基因技術(shù)或基因編輯技術(shù)對這些基因的功能和作用機制進一步驗證。

    參考文獻

    [1]? Amorim L L B, Da Fonseca Dos Santos R, Neto J P B, et al. Transcription factors involved in plant resistance to pathogens[J]. Current Protein and Peptide Science, 2017, 18(4): 335-351.

    [2]? Ng D W, Abeysinghe J K, Kamali M. Regulating the regulators: the control of transcription factors in plant defense signaling[J]. International Journal of Molecular Sciences, 2018, 19(12): 3737.

    [3]? Gaudinier A, Rodriguez-Medina J, Zhang L, et al. Transcriptional regulation of nitrogen-associated metabolism and growth[J]. Nature, 2018, 563(7730): 259-264.

    [4]? Debernardi J M, Mecchia M A, Vercruyssen L, et al. Post-transcriptional control of GRF transcription factors by microRNA miR396 and GIF co-activator affects leaf size and longevity[J]. Plant Journal, 2014, 79(3): 413-426.

    [5]? Baucher M, Moussawi J, Vandeputte O M, et al. A role for the miR396/GRF network in specification of organ type during flower development, as supported by ectopic expression of Populus trichocarpa miR396c in transgenic tobacco[J]. Plant Biology, 2013, 15(5): 892-898.

    [6]? Fonini L S, Lazzarotto F, Barros P M, et al. Molecular evolution and diversification of the GRF transcription factor family[J]. Genetics and Molecular Biology, 2020, 43(3): e20200080.

    [7]? Kim J H. Biological roles and an evolutionary sketch of the GRF-GIF transcriptional complex in plants[J]. BMB Reports, 2019, 52(4): 227-238.

    [8]? Van Der Knaap E, Kim J H, Kende H. A novel gibberellin-induced gene from rice and its potential regulatory role in stem growth[J]. Plant Physiology, 2000, 122(3): 695-704.

    [9]? Sun P, Zhang W, Wang Y, et al. OsGRF4 controls grain shape, panicle length and seed shattering in rice[J]. Journal of Integrative Plant Biology, 2016, 58(10): 836-847.

    [10]????? Vercruyssen L, Tognetti V B, Gonzalez N, et al. GROWTH REGULATING FACTOR5 stimulates Arabidopsis chloroplast division, photosynthesis, and leaf longevity[J]. Plant Physiology, 2015, 167(3): 817-832.

    [11]????? Szczygiel-Sommer A, Gaj M D. The miR396-GRF regulatory module controls the embryogenic response in Arabidopsis via an auxin-related pathway[J]. International Journal of Molecular Sciences, 2019, 20(20): 1-18.

    [12]????? Li S, Gao F, Xie K, et al. The OsmiR396c-OsGRF4-OsGIF1 regulatory module determines grain size and yield in rice[J]. Plant Biotechnology Journal, 2016, 14(11): 2134-2146.

    [13]????? Liu J, Hua W, Yang H L, et al. The BnGRF2 gene (GRF2-like gene from Brassica napus) enhances seed oil production through regulating cell number and plant photosynthesis[J]. Journal of Experimental Botany, 2012, 63(10): 3727-3740.

    [14]????? 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.

    [15]????? Choi D, Kim J H, Kende H. Whole genome analysis of the OsGRF gene family encoding plant-specific putative transcription activators in rice (Oryza sativa L.)[J]. Plant and Cell Physiology, 2004, 45(7): 897-904.

    [16]????? Khatun K, Robin A H K, Park J, et al. Molecular characterization and expression profiling of tomato GRF transcription factor family genes in response to abiotic stresses and phytohormones[J]. International Journal of Molecular Sciences, 2017, 18(5): 1056.

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

    [18]????? Chen F, Yang Y, Luo X, et al. Genome-wide identification of GRF transcription factors in soybean and expression analysis of GmGRF family under shade stress[J]. BMC Plant Biology, 2019, 19(1): 269.

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

    [20]????? Cao J F, Huang J Q, Liu X, et al. Genome-wide characterization of the GRF family and their roles in response to salt stress in Gossypium[J]. BMC Genomics, 2020, 21(1): 575.

    [21]????? 劉蒲東, 張舒婷, 陳曉慧, 等. 龍眼GRF家族全基因組鑒定及表達模式[J]. 應(yīng)用與環(huán)境生物學(xué)報, 2020, 26(2): 236-245.

    [22]????? 雷新濤, 曹紅星, 馮美利, 等. 熱帶木本生物質(zhì)能源樹種——油棕[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報, 2012, 17(6): 185-190.

    [23]????? Mahlia T M I, Ismail N, Hossain N, et al. Palm oil and its wastes as bioenergy sources: a comprehensive review[J]. Environmental Science and Pollution Research International, 2019, 26(15): 14849-14866.

    [24]????? Cros D, Tchounke B, Nkague-Nkamba L. Training genomic selection models across several breeding cycles increases genetic gain in oil palm in silico study[J]. Molecular Breeding, 2018, 38(7): 1-12.

    [25]????? Lu S, Wang J, Chitsaz F, et al. CDD/SPARCLE: the conserved domain database in 2020[J]. Nucleic Acids Research, 2020, 48(d1): D265-D268.

    [26]????? Chen C, Chen H, Zhang Y, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data[J]. Molecular Plant, 2020, 13(8): 1194-1202.

    [27]????? Edgar R C, Batzoglou S. Multiple sequence alignment[J]. Current Opinion in Structural Biology, 2006, 16(3): 368-373.

    [28]????? Lescot M, Déhais P, Thijs G, et al. Plant CARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences[J]. Nucleic Acids Research, 2002, 30(1): 325-327.

    [29]????? Panchy N, Lehti-Shiu M, Shiu S H. Evolution of gene duplication in plants[J]. Plant Physiology, 2016, 171(4): 2294-2316.

    [30]????? Huff J T, Zilberman D, Roy S W. Mechanism for DNA transposons to generate introns on genomic scales[J]. Nature, 2016, 538(7626): 533-536.

    [31]????? Rodriguez R E, Mecchia M A, Debernardi J M, et al. Control of cell proliferation in Arabidopsis thaliana by microRNA miR396[J]. Development, 2010, 137(1): 103-112.

    [32]????? Ercoli M F, Rojas A M, Debernardi J M, et al. Control of cell proliferation and elongation by miR396[J]. Plant Signaling & Behavior, 2016, 11(6): e1184809.

    [33]????? Kim J H, Lee B H. GROWTH-REGULATING FACTOR4 of Arabidopsis thaliana is required for development of leaves, cotyledons, and shoot apical meristem[J]. Journal of Plant Biology, 2006, 49(6): 463-468.

    [34]????? Zhang D, Sun W, Singh R, et al. GRF-interacting factor1 regulates shoot architecture and meristem determinacy in maize[J]. Plant Cell, 2018, 30(2): 360-374.

    [35]????? Chao Q, Gao Z F, Zhang D, et al. The developmental dynamics of the Populus stem transcriptome[J]. Plant Biotechnology Journal, 2019, 17(1): 206-219.

    [36]????? Samad A F A, Sajad M, Nazaruddin N, et al. MicroRNA and transcription factor: key players in plant regulatory network[J]. Front Plant Science, 2017, 8: 565.

    [37]????? Yang F X, Liang G, Liu D M, et al. Arabidopsis MiR396 mediates the development of leaves and flowers in transgenic tobacco[J]. Journal of Plant Biology, 2009, 52(5): 475-481.

    [38]????? Baucher M, Moussawi J, Vandeputte O M, et al. A role for the miR396/GRF network in specification of organ type during flower development, as supported by ectopic expression of Populus trichocarpa miR396c in transgenic tobacco[J]. Plant Biology, 2013, 15(5): 892-898.

    [39]????? Liu X, Guo L X, Jin L F, et al. Identification and transcript profiles of Citrus growth-regulating factor genes involved in the regulation of leaf and fruit development[J]. Molecular Biology Reports, 2016, 43(10): 1059-1067.

    [40]????? Tranbarger T J, Stéphane D, Thierry J, et al. Regulatory mechanisms underlying oil palm fruit mesocarp maturation, ripening, and functional specialization in lipid and carotenoid metabolism[J]. Plant Physiology, 2011, 156(2): 564-584.

    責(zé)任編輯:黃東杰

    猜你喜歡
    油棕基因表達生物信息學(xué)
    棕櫚栽培管理初探
    抗菌肽對細菌作用機制的研究
    基因芯片在胃癌及腫瘤球細胞差異表達基因篩選中的應(yīng)用
    “PBL+E—learning”教學(xué)模式探索
    移動教學(xué)在生物信息學(xué)課程改革中的應(yīng)用
    今傳媒(2016年11期)2016-12-19 11:35:50
    中醫(yī)大數(shù)據(jù)下生物信息學(xué)的發(fā)展及教育模式淺析
    數(shù)據(jù)挖掘技術(shù)在生物信息學(xué)中的應(yīng)用
    美洲大蠊提取液對大鼠難愈合創(chuàng)面VEGF表達影響的研究
    二甲基砷酸毒理學(xué)的研究進展
    泰國油棕果市價連跌2個月
    亚洲av中文字字幕乱码综合| 国产乱人偷精品视频| 狠狠精品人妻久久久久久综合| 日本三级黄在线观看| 中文资源天堂在线| 亚洲激情五月婷婷啪啪| 欧美精品国产亚洲| 禁无遮挡网站| 伊人久久国产一区二区| 亚洲色图av天堂| 日本色播在线视频| 搡老乐熟女国产| 精品一区二区三区视频在线| 欧美一级a爱片免费观看看| 波野结衣二区三区在线| 日韩伦理黄色片| 久久久精品欧美日韩精品| 人妻少妇偷人精品九色| 男女无遮挡免费网站观看| 国产精品久久久久久久电影| 亚洲国产精品专区欧美| 18禁在线无遮挡免费观看视频| 欧美 日韩 精品 国产| kizo精华| 亚洲精品久久久久久婷婷小说| 国产黄a三级三级三级人| 日本一二三区视频观看| 日本一本二区三区精品| 成年女人看的毛片在线观看| 久久久久久久久大av| 91久久精品国产一区二区成人| 国产精品一二三区在线看| 精品久久久久久久人妻蜜臀av| 美女脱内裤让男人舔精品视频| 久久久久精品久久久久真实原创| 久久97久久精品| 超碰97精品在线观看| 91精品一卡2卡3卡4卡| 国产乱人视频| 97在线视频观看| 久久久久精品性色| 大又大粗又爽又黄少妇毛片口| 日本-黄色视频高清免费观看| 直男gayav资源| 亚洲欧美成人精品一区二区| 国产黄色免费在线视频| 搡老乐熟女国产| 秋霞在线观看毛片| 免费av观看视频| 久久精品人妻少妇| 激情五月婷婷亚洲| 最近2019中文字幕mv第一页| 国产又色又爽无遮挡免| 久久久亚洲精品成人影院| 一边亲一边摸免费视频| 最近2019中文字幕mv第一页| 97超碰精品成人国产| 亚洲精品视频女| 欧美少妇被猛烈插入视频| 狂野欧美激情性bbbbbb| 欧美老熟妇乱子伦牲交| 2021天堂中文幕一二区在线观| 嫩草影院精品99| 久久精品国产亚洲网站| 国产成人一区二区在线| 精品人妻视频免费看| 色视频在线一区二区三区| 免费黄色在线免费观看| av在线亚洲专区| 亚洲欧美精品专区久久| 91在线精品国自产拍蜜月| 欧美 日韩 精品 国产| 九草在线视频观看| 精品国产三级普通话版| 丰满少妇做爰视频| 成人美女网站在线观看视频| 交换朋友夫妻互换小说| 亚洲精品国产色婷婷电影| av女优亚洲男人天堂| 国产探花在线观看一区二区| 亚洲激情五月婷婷啪啪| 亚洲av一区综合| 亚洲成人一二三区av| 不卡视频在线观看欧美| 久久国产乱子免费精品| 91久久精品国产一区二区成人| 欧美成人午夜免费资源| 精品国产露脸久久av麻豆| 丰满人妻一区二区三区视频av| 午夜福利视频1000在线观看| 亚洲丝袜综合中文字幕| 免费看a级黄色片| 欧美日韩视频高清一区二区三区二| 亚洲精品中文字幕在线视频 | 国产精品久久久久久av不卡| 亚洲天堂国产精品一区在线| 日本黄大片高清| 色哟哟·www| 午夜激情福利司机影院| 欧美老熟妇乱子伦牲交| 精品国产乱码久久久久久小说| 91在线精品国自产拍蜜月| 成年av动漫网址| 黄色日韩在线| 中国三级夫妇交换| 熟女人妻精品中文字幕| 人妻一区二区av| 2022亚洲国产成人精品| 一本色道久久久久久精品综合| 麻豆国产97在线/欧美| 亚洲内射少妇av| 大片电影免费在线观看免费| 中国三级夫妇交换| 国产精品熟女久久久久浪| 大又大粗又爽又黄少妇毛片口| 91午夜精品亚洲一区二区三区| 听说在线观看完整版免费高清| 又爽又黄无遮挡网站| 国产成人精品福利久久| 中国三级夫妇交换| 下体分泌物呈黄色| 欧美+日韩+精品| 在线观看一区二区三区| 最近手机中文字幕大全| av专区在线播放| 国内精品美女久久久久久| 一级a做视频免费观看| 韩国av在线不卡| 国产成人freesex在线| 少妇人妻 视频| 大片免费播放器 马上看| 一级爰片在线观看| 国产黄色免费在线视频| 卡戴珊不雅视频在线播放| 涩涩av久久男人的天堂| 日韩av在线免费看完整版不卡| 久久人人爽av亚洲精品天堂 | 亚洲在久久综合| 又大又黄又爽视频免费| 99久国产av精品国产电影| 久久久精品免费免费高清| 国产精品成人在线| 97热精品久久久久久| 天堂网av新在线| 亚洲精品自拍成人| av国产精品久久久久影院| 亚洲天堂国产精品一区在线| 亚洲av中文av极速乱| 欧美性猛交╳xxx乱大交人| 欧美成人精品欧美一级黄| 国产高清不卡午夜福利| 夫妻午夜视频| 香蕉精品网在线| 最近中文字幕2019免费版| 亚洲国产欧美在线一区| 热re99久久精品国产66热6| 看十八女毛片水多多多| 2018国产大陆天天弄谢| 蜜桃亚洲精品一区二区三区| 日本免费在线观看一区| 偷拍熟女少妇极品色| 九九爱精品视频在线观看| 国内精品美女久久久久久| 中国国产av一级| 亚洲人成网站在线播| 日本免费在线观看一区| 在线观看国产h片| 中文字幕制服av| 在线天堂最新版资源| 夫妻午夜视频| 欧美高清性xxxxhd video| 日本黄大片高清| 一级二级三级毛片免费看| 99热这里只有精品一区| 天天躁夜夜躁狠狠久久av| 日韩成人av中文字幕在线观看| 新久久久久国产一级毛片| av在线蜜桃| 久久综合国产亚洲精品| 免费黄频网站在线观看国产| 国产大屁股一区二区在线视频| 中文天堂在线官网| 69av精品久久久久久| 亚洲人成网站高清观看| 黄片无遮挡物在线观看| 免费观看的影片在线观看| 国产精品不卡视频一区二区| 国产精品久久久久久精品电影| 乱系列少妇在线播放| 国产欧美亚洲国产| 色哟哟·www| 国产欧美日韩一区二区三区在线 | 久久久久性生活片| 99视频精品全部免费 在线| 欧美少妇被猛烈插入视频| 狠狠精品人妻久久久久久综合| 麻豆久久精品国产亚洲av| 午夜亚洲福利在线播放| 欧美97在线视频| 国产爱豆传媒在线观看| 亚洲内射少妇av| 久久97久久精品| 色视频www国产| 亚洲三级黄色毛片| 我的老师免费观看完整版| 国产一区二区三区av在线| 国国产精品蜜臀av免费| 男人狂女人下面高潮的视频| 日本黄大片高清| 天美传媒精品一区二区| 黄色欧美视频在线观看| 黄片wwwwww| 精品一区二区三区视频在线| 成人免费观看视频高清| 国产高清国产精品国产三级 | 成年版毛片免费区| 国产高清国产精品国产三级 | 你懂的网址亚洲精品在线观看| 国产精品精品国产色婷婷| 伊人久久国产一区二区| 中文欧美无线码| 男插女下体视频免费在线播放| 日韩成人伦理影院| 亚州av有码| 极品教师在线视频| 久久热精品热| 婷婷色综合www| 欧美激情国产日韩精品一区| 一边亲一边摸免费视频| 伊人久久国产一区二区| 亚洲熟女精品中文字幕| 狂野欧美白嫩少妇大欣赏| 国产黄片美女视频| 亚洲av电影在线观看一区二区三区 | 午夜精品一区二区三区免费看| 国产一区二区三区av在线| 国产精品嫩草影院av在线观看| 十八禁网站网址无遮挡 | 亚洲怡红院男人天堂| 久久人人爽人人片av| 老司机影院毛片| 中国国产av一级| 国产一区有黄有色的免费视频| 国产 一区精品| 亚洲国产色片| 大香蕉久久网| 久久影院123| 男女那种视频在线观看| 在线观看一区二区三区| 日本黄色片子视频| 成人毛片60女人毛片免费| 一二三四中文在线观看免费高清| 97超碰精品成人国产| 日韩三级伦理在线观看| 别揉我奶头 嗯啊视频| 一区二区三区精品91| 在现免费观看毛片| 成人免费观看视频高清| 日韩成人av中文字幕在线观看| 一级毛片aaaaaa免费看小| 搡女人真爽免费视频火全软件| 亚洲欧美日韩另类电影网站 | 日日啪夜夜撸| 亚洲精品日本国产第一区| 少妇高潮的动态图| 99九九线精品视频在线观看视频| 精品人妻一区二区三区麻豆| 欧美精品人与动牲交sv欧美| av在线观看视频网站免费| 国产毛片a区久久久久| 国产亚洲5aaaaa淫片| 精品一区二区免费观看| 插阴视频在线观看视频| 天堂俺去俺来也www色官网| 国产成人一区二区在线| 亚洲成人中文字幕在线播放| 少妇猛男粗大的猛烈进出视频 | 在线观看av片永久免费下载| 少妇裸体淫交视频免费看高清| av卡一久久| 人妻 亚洲 视频| 国产精品av视频在线免费观看| 欧美成人一区二区免费高清观看| 小蜜桃在线观看免费完整版高清| 欧美bdsm另类| 熟女人妻精品中文字幕| av天堂中文字幕网| 国产男人的电影天堂91| 18禁裸乳无遮挡免费网站照片| 蜜桃亚洲精品一区二区三区| 久久久久久九九精品二区国产| 亚洲精品国产色婷婷电影| av黄色大香蕉| 久久99蜜桃精品久久| 亚洲婷婷狠狠爱综合网| 亚洲天堂国产精品一区在线| 亚洲精品自拍成人| 成人无遮挡网站| 午夜免费男女啪啪视频观看| 亚洲欧美日韩无卡精品| 丰满少妇做爰视频| 精品一区二区三卡| 少妇人妻 视频| 狂野欧美激情性bbbbbb| 在线a可以看的网站| 久久人人爽av亚洲精品天堂 | 久久这里有精品视频免费| 少妇人妻一区二区三区视频| 日本欧美国产在线视频| 99re6热这里在线精品视频| 一级毛片 在线播放| 韩国av在线不卡| 成年女人在线观看亚洲视频 | 三级国产精品片| 熟女电影av网| 美女xxoo啪啪120秒动态图| 亚洲不卡免费看| 国产成人精品久久久久久| 狂野欧美激情性xxxx在线观看| 免费看日本二区| 国产av码专区亚洲av| 黄色怎么调成土黄色| 一级毛片aaaaaa免费看小| 大又大粗又爽又黄少妇毛片口| 91午夜精品亚洲一区二区三区| 一区二区三区乱码不卡18| 伊人久久精品亚洲午夜| www.色视频.com| 青春草国产在线视频| 欧美三级亚洲精品| 国产综合精华液| 亚洲人成网站高清观看| 只有这里有精品99| 国产免费一区二区三区四区乱码| 亚洲丝袜综合中文字幕| 哪个播放器可以免费观看大片| 一区二区av电影网| 成人毛片60女人毛片免费| 久久ye,这里只有精品| 国产精品蜜桃在线观看| 国产午夜福利久久久久久| 午夜激情久久久久久久| 成人亚洲精品av一区二区| 亚洲av中文av极速乱| 日本欧美国产在线视频| 女人久久www免费人成看片| 成人漫画全彩无遮挡| 国产老妇伦熟女老妇高清| 午夜福利视频1000在线观看| 欧美日本视频| 一本色道久久久久久精品综合| 欧美丝袜亚洲另类| 国产男女超爽视频在线观看| 亚洲最大成人av| 欧美成人精品欧美一级黄| av天堂中文字幕网| 亚洲国产欧美在线一区| 亚洲精品日本国产第一区| 国产又色又爽无遮挡免| 哪个播放器可以免费观看大片| 色播亚洲综合网| 国产成人a区在线观看| 2022亚洲国产成人精品| 久久97久久精品| 欧美97在线视频| 日产精品乱码卡一卡2卡三| 91精品伊人久久大香线蕉| 联通29元200g的流量卡| 2022亚洲国产成人精品| 一级二级三级毛片免费看| 久久精品夜色国产| 三级男女做爰猛烈吃奶摸视频| 色哟哟·www| 内地一区二区视频在线| 99热网站在线观看| 观看免费一级毛片| 精品午夜福利在线看| 国产免费又黄又爽又色| 黑人高潮一二区| 亚洲av不卡在线观看| 中文字幕av成人在线电影| 亚洲av免费高清在线观看| 人体艺术视频欧美日本| 欧美zozozo另类| 搞女人的毛片| 欧美性猛交╳xxx乱大交人| 美女被艹到高潮喷水动态| 久久久久久伊人网av| 18禁动态无遮挡网站| 亚洲av中文av极速乱| 一区二区三区乱码不卡18| 亚洲无线观看免费| 国产精品国产三级国产av玫瑰| 欧美丝袜亚洲另类| 国产午夜福利久久久久久| 一二三四中文在线观看免费高清| 人妻 亚洲 视频| 亚洲伊人久久精品综合| a级毛色黄片| 天堂中文最新版在线下载 | 日本欧美国产在线视频| 看十八女毛片水多多多| 中国国产av一级| 国产成人a区在线观看| 久久久久精品久久久久真实原创| 日韩一区二区三区影片| 国产免费福利视频在线观看| 美女内射精品一级片tv| 日韩国内少妇激情av| 久久精品综合一区二区三区| 美女被艹到高潮喷水动态| 成人特级av手机在线观看| 成人一区二区视频在线观看| 国产亚洲最大av| 色5月婷婷丁香| 免费观看性生交大片5| 成年版毛片免费区| 99久久精品国产国产毛片| 永久免费av网站大全| 亚洲成人久久爱视频| 免费高清在线观看视频在线观看| 亚洲国产成人一精品久久久| 我的女老师完整版在线观看| 精品一区二区三卡| 欧美性猛交╳xxx乱大交人| 不卡视频在线观看欧美| 欧美日韩国产mv在线观看视频 | 免费不卡的大黄色大毛片视频在线观看| 免费少妇av软件| 色哟哟·www| 尾随美女入室| 一二三四中文在线观看免费高清| 国产精品久久久久久久电影| 亚洲欧美一区二区三区国产| 少妇人妻一区二区三区视频| 伊人久久国产一区二区| 2018国产大陆天天弄谢| 人妻一区二区av| 丝袜美腿在线中文| 亚洲电影在线观看av| 麻豆成人午夜福利视频| 国产精品国产三级专区第一集| 国产精品久久久久久久久免| 高清日韩中文字幕在线| 亚洲av中文字字幕乱码综合| 久久久久久久精品精品| 美女脱内裤让男人舔精品视频| 中文乱码字字幕精品一区二区三区| 亚洲一区二区三区欧美精品 | 免费不卡的大黄色大毛片视频在线观看| 免费人成在线观看视频色| 一边亲一边摸免费视频| 欧美97在线视频| 欧美精品国产亚洲| 国产成人午夜福利电影在线观看| 亚洲欧美日韩卡通动漫| 日本免费在线观看一区| 一级黄片播放器| 美女主播在线视频| 亚洲国产成人一精品久久久| 成人黄色视频免费在线看| 天美传媒精品一区二区| 插阴视频在线观看视频| 特级一级黄色大片| 狂野欧美激情性xxxx在线观看| 午夜免费观看性视频| 2018国产大陆天天弄谢| 一级爰片在线观看| 国产成人免费观看mmmm| 九九爱精品视频在线观看| 亚洲精品,欧美精品| 美女cb高潮喷水在线观看| 亚洲欧美日韩东京热| 久久ye,这里只有精品| 成人午夜精彩视频在线观看| 欧美精品国产亚洲| 国产色爽女视频免费观看| 亚洲国产精品专区欧美| 内地一区二区视频在线| 少妇被粗大猛烈的视频| 18禁动态无遮挡网站| 深夜a级毛片| 2021少妇久久久久久久久久久| 亚洲成人一二三区av| 免费高清在线观看视频在线观看| 色婷婷久久久亚洲欧美| 熟女电影av网| 亚洲精品久久午夜乱码| 91狼人影院| 久久影院123| 看黄色毛片网站| 精品一区二区免费观看| 日韩av免费高清视频| 国产精品麻豆人妻色哟哟久久| av国产免费在线观看| 男人狂女人下面高潮的视频| 黄片无遮挡物在线观看| 日韩一区二区三区影片| 久久久午夜欧美精品| 丝袜脚勾引网站| 永久免费av网站大全| 国产淫片久久久久久久久| 日韩一区二区视频免费看| 成人免费观看视频高清| 观看美女的网站| 91狼人影院| 一级毛片久久久久久久久女| 国产午夜精品一二区理论片| 精品一区在线观看国产| 看黄色毛片网站| 国产高清国产精品国产三级 | 精品人妻一区二区三区麻豆| 美女国产视频在线观看| av在线app专区| 99久久中文字幕三级久久日本| 在线观看人妻少妇| 一本色道久久久久久精品综合| 男女那种视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 国产高清三级在线| 亚洲自偷自拍三级| 成人一区二区视频在线观看| 男女高潮啪啪啪动态图| 欧美激情 高清一区二区三区| 成人影院久久| 国产一区二区三区综合在线观看| 一级毛片我不卡| 国产黄色免费在线视频| 亚洲精品一区蜜桃| 成人影院久久| 国产熟女欧美一区二区| 热re99久久精品国产66热6| 最近最新中文字幕免费大全7| 亚洲欧美成人精品一区二区| 美国免费a级毛片| 欧美日韩成人在线一区二区| 国产亚洲精品第一综合不卡| 久久人人爽av亚洲精品天堂| 国产在线免费精品| 9191精品国产免费久久| 亚洲在久久综合| 香蕉丝袜av| 欧美日韩国产mv在线观看视频| 国产在线视频一区二区| 黑人欧美特级aaaaaa片| 久久精品aⅴ一区二区三区四区| 国产精品 欧美亚洲| 国产毛片在线视频| 91老司机精品| 亚洲成av片中文字幕在线观看| 激情五月婷婷亚洲| 日韩欧美一区视频在线观看| 国产片内射在线| 午夜福利乱码中文字幕| 成人影院久久| 精品人妻熟女毛片av久久网站| 国产成人免费无遮挡视频| 老司机影院毛片| 中文字幕精品免费在线观看视频| 成人毛片60女人毛片免费| 亚洲精品第二区| 如日韩欧美国产精品一区二区三区| 久久久久久人人人人人| 日本wwww免费看| 亚洲欧美色中文字幕在线| 国产精品三级大全| 尾随美女入室| 少妇猛男粗大的猛烈进出视频| 美女脱内裤让男人舔精品视频| 中文乱码字字幕精品一区二区三区| 国产精品 国内视频| 久久久久久免费高清国产稀缺| 操出白浆在线播放| 国产亚洲av片在线观看秒播厂| 一边亲一边摸免费视频| 亚洲成人国产一区在线观看 | 亚洲,一卡二卡三卡| 亚洲成人av在线免费| 亚洲四区av| 久久性视频一级片| 亚洲国产看品久久| 美女福利国产在线| 国产精品国产三级专区第一集| 国产一区有黄有色的免费视频| 亚洲色图综合在线观看| 色精品久久人妻99蜜桃| 精品亚洲成国产av| 校园人妻丝袜中文字幕| 国产福利在线免费观看视频| 免费黄网站久久成人精品| 搡老乐熟女国产| 51午夜福利影视在线观看| 在线亚洲精品国产二区图片欧美| 激情五月婷婷亚洲| av.在线天堂| 亚洲精品国产av蜜桃| 国产人伦9x9x在线观看| 欧美xxⅹ黑人| 嫩草影视91久久| 9色porny在线观看| 天天操日日干夜夜撸| 极品少妇高潮喷水抽搐| 久久99热这里只频精品6学生| av网站免费在线观看视频| 亚洲少妇的诱惑av| 免费观看性生交大片5| 麻豆av在线久日| 国产成人系列免费观看| 国精品久久久久久国模美| 这个男人来自地球电影免费观看 | 国产日韩一区二区三区精品不卡| 亚洲av日韩在线播放| 99久久人妻综合| 亚洲精品久久午夜乱码| 日韩一本色道免费dvd|