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

    黃藤2個(gè)NAC基因的分子特征及其SSR分子標(biāo)記開發(fā)*

    2017-10-14 08:02:15孫化雨李利超趙韓生楊意宏王思寧高志民
    林業(yè)科學(xué) 2017年8期
    關(guān)鍵詞:基因組位點(diǎn)引物

    孫化雨 李利超 趙韓生 楊意宏,2 王思寧 高志民

    (1. 國(guó)際竹藤中心 國(guó)家林業(yè)局竹藤科學(xué)與技術(shù)重點(diǎn)開放實(shí)驗(yàn)室 北京 100102; 2. 河北農(nóng)業(yè)大學(xué)園藝學(xué)院 保定 071001)

    黃藤2個(gè)NAC基因的分子特征及其SSR分子標(biāo)記開發(fā)*

    孫化雨1李利超1趙韓生1楊意宏1,2王思寧1高志民1

    (1. 國(guó)際竹藤中心 國(guó)家林業(yè)局竹藤科學(xué)與技術(shù)重點(diǎn)開放實(shí)驗(yàn)室 北京 100102; 2. 河北農(nóng)業(yè)大學(xué)園藝學(xué)院 保定 071001)

    【目的】 棕櫚藤是重要的森林植物,纖鞭是其重要的攀援器官,也是重要的分類依據(jù)。研究黃藤NAC(NAM,ATAF和CUC)轉(zhuǎn)錄因子基因的分子特征并開發(fā)SSR分子標(biāo)記,以期為棕櫚藤的分子育種和輔助分類提供參考依據(jù)?!痉椒ā?以黃藤為材料,借助轉(zhuǎn)錄組數(shù)據(jù),采用同源克隆的方法從黃藤中分離NAC基因,采用生物信息學(xué)方法進(jìn)行基因結(jié)構(gòu)、蛋白性質(zhì)與結(jié)構(gòu)分析和SSR位點(diǎn)預(yù)測(cè),采用實(shí)時(shí)定量PCR技術(shù)分析基因的組織表達(dá)特性,利用PAGE電泳和測(cè)序技術(shù)分析SSR分子標(biāo)記在不同棕櫚藤樣品中的通用性和多態(tài)性?!窘Y(jié)果】 從黃藤葉片中獲得了2個(gè)NAC同源基因DjNAC3 (GenBank登錄號(hào): KU556738)和DjNAC4(GenBank登錄號(hào): KX579750),二者的開放閱讀框長(zhǎng)度分別為729 bp和1 326 bp,對(duì)應(yīng)的基因組序列為850 bp和1 441 bp,均包含2個(gè)外顯子和1個(gè)內(nèi)含子。DjNAC3和DjNAC4編碼的蛋白分別為242 aa和441 aa。蛋白結(jié)構(gòu)分析表明,DjNAC3和DjNAC4具有典型的NAC轉(zhuǎn)錄因子結(jié)構(gòu)特征,屬于NAC家族的CUC亞家族,但二者之間的相似系數(shù)僅為23.6%,表明它們?cè)邳S藤生長(zhǎng)發(fā)育過程中可能具有不同的功能。DjNAC3和DjNAC4在不同組織中的表達(dá)模式存在明顯差異,DjNAC3在發(fā)育成熟纖鞭中的表達(dá)豐度最高,葉片中的表達(dá)豐度最低,而DjNAC4則在發(fā)育成熟的鉤刺中表達(dá)豐度最高,而發(fā)育初期的鉤刺中最低。在DjNAC3和DjNAC4的基因組序列中分別包含1個(gè)SSR位點(diǎn),其中前者的SSR位點(diǎn)位于內(nèi)含子區(qū)域,為(TA)6,后者的SSR位點(diǎn)位于第1個(gè)外顯子區(qū)域,為(GCA)5。根據(jù)DjNAC3和DjNAC4中SSR位點(diǎn)旁側(cè)序列設(shè)計(jì)引物,以黃藤和另外20個(gè)不同棕櫚藤樣品的基因組DNA為模板進(jìn)行擴(kuò)增,PAGE電泳結(jié)果分析表明,引物具有較高的通用性,且擴(kuò)增產(chǎn)物具有多態(tài)性。6個(gè)樣品的測(cè)序結(jié)果證實(shí),用DjNAC3設(shè)計(jì)引物的擴(kuò)增產(chǎn)物測(cè)序獲得的SSR位點(diǎn)序列存在著一定的差異,既包括SSR類型變異、重復(fù)次數(shù)變化等多態(tài)性,又有SSR位點(diǎn)缺失的現(xiàn)象; 而根據(jù)DjNAC4設(shè)計(jì)引物擴(kuò)增獲得的SSR位點(diǎn)序列差異主要為重復(fù)次數(shù)的變化?!窘Y(jié)論】 黃藤DjNAC3和DjNAC4基因的基因結(jié)構(gòu)、表達(dá)模式、SSR位點(diǎn)等均存在明顯的差異,這表明它們?cè)邳S藤生長(zhǎng)發(fā)育中可能具有不同的功能,二者所包含的SSR分子標(biāo)記具有通用性和多態(tài)性,可以作為分子標(biāo)記應(yīng)用于棕櫚藤的輔助分類和分子輔助育種。

    黃藤;NAC; 基因表達(dá); SSR開發(fā)

    Abstract: 【Objective】 Rattan is one of the important forest plants, cirri and flagella are important organs for its climbing habit, which are also the important basis for rattan classification. In order to provide a basis for molecular breeding and SSR-assisted classification of rattan, study on the molecular characteristics of NAC (NAM, ATAF and CUC) transcription factor genes inDaemonoropsjenkinsianaand the development of SSR markers was performed in this paper.【Method】 With the aid of transcriptome data, the sequences ofNAChomologue genes inD.jenkinsianawere isolated by PCR method. The analyses of gene structure, protein properties and structure as well as the prediction of SSR locus inNACgenes were conducted using bioinformatics method. The tissue specific expression ofNACgenes in different tissues were analyzed using real-time quantitative PCR (qPCR). To evaluate the universality and polymorphism of the developed SSR markers, PAGE electrophoresis and sequencing analyses were used on the basis of products amplified with different rattan samples.【Result】 Two homologous genes,DjNAC3 (GenBank No. KU556738) andDjNAC4 (GenBank No. KX579750), were obtained from leaves ofD.jenkinsiana, of which open reading frame (ORF) were 729 bp and 1 326 bp, respectively. The genomic sequence corresponding to the ORFs ofDjNAC3 andDjNAC4 were 850 bp and 1 441 bp, which all contained two exons and one intron. The proteins encoded byDjNAC3 andDjNAC4 were 242 aa and 441 aa respectively. Protein structure analysis showed that DjNAC3 and DjNAC4 had the typical structural features of NAC transcription factors, which were belonged to the CUC subfamily of NAC family. However, the similarity coefficient between DjNAC3 and DjNAC4 was only 23.6%, indicating that they might have different functions in the growth process ofD.jenkinsiana. The expression patterns ofDjNAC3 andDjNAC4 in different tissues were obviously different.DjNAC3 was expressed in developed cirri with the highest level, and the lowest in leaves, while that ofDjNAC4 was the highest abundance in developed barbs and the lowest in early developed barbs. The SSR loci were detected in the genomic sequences ofDjNAC3 andDjNAC4, the SSR locus of (TA)6was located in the intron region ofDjNAC3 and that of (GCA)5was in the first exon ofDjNAC4. SSR primers were designed according to the flanking sequences of SSR loci inDjNAC3 andDjNAC4. Genomic DNAs ofD.jenkinsianaand other 20 different rattan samples were selected as templates for amplification. PAGE electrophoresis analysis showed that the primers had universality and polymorphism in the samples. The sequencing result of amplification products from six different templates further confirmed the polymorphism such as variation of SSR type, the number change of repetitions, and SSR locus missing, which were found in the sequences generated by the primer pair designed based onDjNAC3. The main difference in the sequences generated by the primer pair designed according toDjNAC4 was the number change of repetitions.【Conclusion】 There are significant differences betweenDjNAC3 andDjNAC4, such as gene structures, gene expression patterns, and SSR locus, which means thatDjNAC3 andDjNAC4 might play different roles in the growth and development ofD.jenkinsiana. The universality and polymorphism of SSR markers developed fromDjNAC3 andDjNAC4 indicate that they can serve as molecular markers for rattan classification and molecular-assisted breeding.

    Keywords:Daemonoropsjenkinsiana;NAC; gene expression; SSR development

    NAC (NAM,ATAF和CUC)是植物所特有的轉(zhuǎn)錄因子家族,其成員數(shù)量眾多。NAC家族成員的N端相對(duì)保守,都由約150個(gè)氨基酸組成結(jié)構(gòu)域,且被分為5個(gè)亞結(jié)構(gòu)域(A、B、C、D和E),但其C端的序列呈現(xiàn)出高度多樣化,富含絲氨酸、谷氨酸等(Olsenetal., 2005)。迄今為止,對(duì)許多物種全基因組范圍內(nèi)的NAC轉(zhuǎn)錄因子進(jìn)行了鑒定,包括擬南芥(Arabidopsisthaliana)(117個(gè))、水稻(Oryzasativa)(151個(gè))(Nuruzzamanetal., 2010)、大豆(Glycinemax)(152個(gè))(Leetal., 2011)、白菜(Brassicarapa)(204個(gè))(Liuetal., 2014)、玉米(Zeamays)(152個(gè))(Shirigaetal., 2014)、鷹嘴豆(Cicerarietinum)(71個(gè))(Haetal., 2014)、番茄(Solanumlycopersicum)(104個(gè))(Suetal., 2015)和毛竹(Phyllostachysedulis)(125個(gè))(黎幫勇等, 2015)等。NAC轉(zhuǎn)錄因子的功能多樣,參與許多生物學(xué)過程,包括花發(fā)育、次生壁形成、細(xì)胞分裂、頂端分生組織的形成、葉片衰老,以及生物和非生物脅迫應(yīng)答(Olsenetal., 2005; Tranetal., 2010; Nakashimaetal., 2012; Nuruzzamanetal., 2013; Banerjeeetal., 2015)等,在植物的生長(zhǎng)發(fā)育過程中起著重要的調(diào)控作用。

    隨著研究的不斷深入,越來越多NAC在植物中的功能被揭示。NAM基因在矮牽牛(Petuniahybrida)的分生組織和原基邊界表達(dá),是其胚胎和花形成所需的(Soueretal., 1996); 擬南芥CUC2的表達(dá)決定其葉緣鋸齒的程度(Nikovicsetal., 2006); 玉米ZmCUC3在頂端分生組織外圍表達(dá),邊界明顯,部分細(xì)胞將形成一個(gè)新的葉原基(Zimmermannetal., 2005)。擬南芥NAC轉(zhuǎn)錄因子NTL8參與赤霉酸介導(dǎo)的鹽信號(hào)調(diào)控(Kimetal., 2008),而且通過FLOWERING LOCUST調(diào)節(jié)擬南芥鹽響應(yīng)開花(Kimetal., 2007),冷激活NAC轉(zhuǎn)錄因子調(diào)控誘導(dǎo)擬南芥病原體的抗性反應(yīng)(Seoetal., 2010),AtNAC2作為乙烯和生長(zhǎng)素信號(hào)轉(zhuǎn)導(dǎo)通路下游的轉(zhuǎn)錄因子,參與鹽應(yīng)激反應(yīng)和側(cè)根發(fā)育調(diào)控(Heetal., 2005)。擬南芥NAC轉(zhuǎn)錄因子AtNAP在葉片衰老中發(fā)揮著重要作用,AtNAP的誘導(dǎo)表達(dá)導(dǎo)致擬南芥早熟衰老(Guoetal., 2006),ABA-AtNAP通過控制氣孔的開閉來調(diào)節(jié)乙烯刺激呼吸,實(shí)現(xiàn)對(duì)果實(shí)的衰老調(diào)控(Kouetal., 2012)。在祖野麥(Triticumturgidumssp.dicoccoides)中NAC轉(zhuǎn)錄因子(NAM-B1)可以加速衰老,提高葉片的營(yíng)養(yǎng)轉(zhuǎn)運(yùn)形成籽粒(Uauyetal., 2006)。水稻OsNAC19是一種轉(zhuǎn)錄激活因子,參與水稻響應(yīng)稻瘟病菌感染,可能在茉莉酸甲酯介導(dǎo)的信號(hào)傳導(dǎo)途徑中起重要作用(Linetal., 2007)。歐洲油菜(Brassicanapus)BnaNAC55通過激活一些活性氧(Reactive Oxygen Species, ROS)以及與防御有關(guān)的基因的表達(dá),來調(diào)控ROS的積累和細(xì)胞死亡(Niuetal., 2016)。

    棕櫚藤屬于棕櫚科(Palmae)植物,是重要的熱帶森林資源,其藤莖是優(yōu)良的非木材森林產(chǎn)品之一,是編制各種高檔家具及工藝品的理想材料,具有很高的經(jīng)濟(jì)價(jià)值。由于棕櫚藤資源極其匱乏,其原料是國(guó)內(nèi)外市場(chǎng)上的緊缺物資。絕大多數(shù)棕櫚藤具有極強(qiáng)的攀緣性,總是和其他的樹木糾纏在一起,使其成為森林的伴生物種。黃藤(Daemonoropsjenkinsiana)為棕櫚科鱗果亞科(Lepidocaryoideae)省藤族(CALAMEAE)黃藤屬(Daemonorops)植物,分布于廣東東南部、香港、海南及廣西西南部,為我國(guó)棕櫚科單屬單種植物,葉羽狀全裂,莖初時(shí)直立,后攀援,其攀援主要借助于葉片頂端延伸而成的具爪狀刺的纖鞭。研究表明NAC參與葉片形態(tài)建成的調(diào)控(Zimmermannetal., 2005; Nikovicsetal., 2006),目前雖然對(duì)黃藤的生長(zhǎng)特性以及形態(tài)特征具有比較詳細(xì)的描述(江澤慧等, 2013),但關(guān)于黃藤的生長(zhǎng)發(fā)育的分子基礎(chǔ)尚屬空白。因此,本研究以黃藤為對(duì)象,從葉片中分離克隆了2個(gè)NAC基因,并對(duì)其分子特征進(jìn)行深入分析,研究其表達(dá)模式,并開發(fā)SSR分子標(biāo)記,以期為深入研究其功能和分子標(biāo)記輔助育種提供參考。

    1 材料與方法

    1.1 試驗(yàn)材料 黃藤樣品取自中國(guó)林業(yè)科學(xué)研究院熱帶林業(yè)研究所(廣州,龍洞),分別采集葉片、發(fā)育初期的纖鞭、發(fā)育成熟的纖鞭、發(fā)育初期的鉤刺、發(fā)育成熟的鉤刺,放入RNAlater? RNA Stabilization Reagent(QIAGEN,德國(guó))中,帶回實(shí)驗(yàn)室備用。另外從中國(guó)科學(xué)院西雙版納植物園分別采集的20個(gè)藤種(表1)葉片,用硅膠干燥處理后帶回實(shí)驗(yàn)室備用。

    1.2 DNA提取、RNA分離與cDNA的合成 采用改良CTAB法提取不同藤種的基因組DNA(高志民等, 2006),存-20 ℃?zhèn)溆?。采用改良Trizol法提取黃藤葉片、發(fā)育初期鉤刺、發(fā)育成熟鉤刺、發(fā)育初期纖鞭和發(fā)育成熟纖鞭的總RNA(Gaoetal., 2006),按照反轉(zhuǎn)錄試劑盒(Promega, 美國(guó))說明書分別合成cDNA。

    表1 本研究所用藤種Tab.1 The list of rattan accessions used in this study

    1.3 基因克隆與分析 根據(jù)黃藤轉(zhuǎn)錄組數(shù)據(jù)(NCBI SRA 注冊(cè)號(hào): SRR3089417、SRR3089429、SRR3089432、SRR3089433、SRR3089434、SRR3089435、SRR3089436和SRR3089437)中預(yù)測(cè)的NAC轉(zhuǎn)錄因子基因序列設(shè)計(jì)引物,Teng-NAC3-F和Teng-NAC3-R,以及Teng-NAC4-F和Teng-NAC4-R(表2),由生工生物工程股份有限公司合成。分別以黃藤cDNA和基因組DNA為模板進(jìn)行擴(kuò)增。PCR反應(yīng)體系為: 5×Prime STARTMBuffer 4 μL,dNTP Mixture (2.5 mmol·L-1)2 μL,正向引物(10 μmol·L-1)0.5 μL,反向引物(10 μmol·L-1)0.5 μL,cDNA/基因組DNA 2 μL,Prime STAR HS DNA polymerase 0.2 μL,DMSO 1 μL,加水補(bǔ)至總體積為20 μL。PCR反應(yīng)程序?yàn)椋?98 ℃ 4 min; 98 ℃ 30 s,60 ℃ 30 s,72 ℃ 1 min 50 s(基因組DNA擴(kuò)增延伸為2 min 30 s),反應(yīng)循環(huán)35個(gè); 72 ℃ 10 min,4 ℃保存。用膠回收試劑盒(Biomega, 美國(guó))回收PCR產(chǎn)物,進(jìn)行加A反應(yīng),反應(yīng)體系為: 10×LA PCR Buffer 1.0 μL,dATP(2.5 mmol·L-1)1.0 μL,回收產(chǎn)物7.9 μL,LATaq酶0.1 μL; 反應(yīng)程序: 70 ℃ 30 min。用加A產(chǎn)物連接 pGEM-T easy載體(Promega, 美國(guó)),轉(zhuǎn)化大腸桿菌(Escherichiacoli)DH5α,選取陽(yáng)性克隆,酶切檢測(cè)分析后,送生工生物工程股份有限公司測(cè)序。

    用DNAstar對(duì)測(cè)序獲得的序列進(jìn)行初步分析。用NCBI在線軟件BLAST (http:∥blast.ncbi.nlm.nih.gov/Blast.cgi)程序進(jìn)行序列同源性比較。蛋白結(jié)構(gòu)域預(yù)測(cè)用TMHMM 2.0 (http:∥myhits.isb-sib.ch/cgi-bin/motif_scan)來完成。對(duì)獲得基因編碼的蛋白采用Phyre2.0 (http:∥www.sbg.bio.ic.ac.uk/~phyre2/html/page.cgi?id=index)進(jìn)行三級(jí)結(jié)構(gòu)建模。

    表2 PCR擴(kuò)增所用引物Tab.2 List of PCR primer sequences

    1.4 基因表達(dá)分析 分別利用定量引物qNAC3-F和qNAC3-R、qNAC4-F和qNAC4-R(表2)對(duì)2個(gè)NAC基因在葉片、發(fā)育初期的纖鞭、發(fā)育成熟的纖鞭、發(fā)育初期的鉤刺、發(fā)育成熟的鉤刺中的表達(dá)變化情況進(jìn)行定量分析。反應(yīng)在耶拿qTOWER2.2儀器上進(jìn)行,總體系(10.0 μL): 5.0 μL LightCycler?480 SYBR Green I Master Mix (Roche,美國(guó)),0.8 μL cDNA,正、反向引物各0.2 μL,3.8 μL ddH2O。反應(yīng)程序: 95 ℃ 5 min; 95 ℃ 30 s,60 ℃ 10 s,40個(gè)循環(huán)。定量分析所用內(nèi)參基因?yàn)棣?tubulin(GenBank No. KX793703)。反應(yīng)重復(fù)3次,結(jié)果用2-ΔΔCt法(Livaketal., 2001)分析。

    1.5 SSR位點(diǎn)分析 利用SSRIT(Simple Sequence Repeat Identification Tool)在線軟件 (Temnykhetal., 2001)對(duì)獲得的黃藤NAC轉(zhuǎn)錄因子基因序列進(jìn)行SSR搜索,SSR類型最多核苷酸數(shù)量設(shè)置為4個(gè),重復(fù)次數(shù)等于或大于5。

    根據(jù)DjNAC3和DjNAC4基因組序列分別設(shè)計(jì)SSR引物(ssrNAC3-F和ssrNAC3-R、ssrNAC4-F和ssrNAC4-R),分別以21個(gè)不同藤種(表1)的基因組DNA為模板,進(jìn)行擴(kuò)增,PCR產(chǎn)物用8%的PAGE電泳進(jìn)行分析。回收以伸長(zhǎng)鉤葉藤、老撾鉤葉藤、狹葉黃藤、黃藤屬種1、勐臘鞭藤、斑嶺省藤6個(gè)藤種基因組DNA為模板擴(kuò)增獲得的目的片段,并連接到pGEM-T easy載體,轉(zhuǎn)化大腸桿菌,并送公司測(cè)序。

    2 結(jié)果與分析

    2.1DjNAC的克隆與基因結(jié)構(gòu)分析 對(duì)應(yīng)黃藤cDNA和基因組DNA模板,引物Teng-NAC3-F和Teng-NAC3-R的PCR擴(kuò)增產(chǎn)物電泳結(jié)果顯示,分別約在0.7 kb和0.8 kb有1條目的條帶,同時(shí)引物Teng-NAC4-F和Teng-NAC4-R分別約在1.3 kb和1.4 kb有1條目的條帶。分別回收目的條帶測(cè)序結(jié)果表明,引物Teng-NAC3-F和Teng-NAC3-R擴(kuò)增的cDNA為731 bp,包含1個(gè)729 bp的開放閱讀(ORF),ORF對(duì)應(yīng)的基因組序列為850 bp,包含2個(gè)外顯子和1個(gè)121 bp內(nèi)含子(454-574 bp),將該基因cDNA命名為DjNAC3 (GenBank登錄號(hào): KU556738)。引物Teng-NAC4-F和Teng-NAC4-R擴(kuò)增的cDNA為1 334 bp,包含1個(gè)1 326 bp的ORF,ORF對(duì)應(yīng)的基因組序列為1 441 bp,包含2個(gè)外顯子和1個(gè)107 bp的內(nèi)含子(875-981 bp),將該基因cDNA命名為DjNAC4 (GenBank登錄號(hào): KX579750)。DjNAC3和DjNAC4的開放閱讀(ORF)對(duì)應(yīng)的基因組所包含的內(nèi)含子均符合GT-AG剪切原則(Mooreetal., 1993),基因結(jié)構(gòu)如圖1所示。

    2.2 DjNAC蛋白性質(zhì)與結(jié)構(gòu)分析DjNAC3編碼1個(gè)242 aa的蛋白,預(yù)測(cè)的分子量為27.95 kDa,等電點(diǎn)為8.323。DjNAC4編碼1個(gè)441 aa的蛋白,預(yù)測(cè)的分子量為49.66 kDa,等電點(diǎn)為7.082。Blast分析表明,DjNAC3和DjNAC4具有典型的NAC轉(zhuǎn)錄因子結(jié)構(gòu)特征,N端保守性較強(qiáng),由大約150個(gè)高度保守的氨基酸殘基組成,包含A、B、C、D、E 5個(gè)亞結(jié)構(gòu)域,C端的氨基酸序列保守性較弱。DjNAC3和DjNAC4均屬于CUC (development-related NAC)亞家族(Zhuetal., 2015),但二者之間的相似系數(shù)僅為23.6%,以c3ulxA_為模板構(gòu)建DjNAC3和DjNAC4的三級(jí)結(jié)構(gòu)(圖2),結(jié)構(gòu)上的差異意味著它們?cè)邳S藤生長(zhǎng)發(fā)育過程中可能具有不同的功能。

    圖1 DjNAC3和DjNAC4的基因結(jié)構(gòu)Fig.1 Gene structures of DjNAC3 and DjNAC4E: 外顯子; I: 內(nèi)含子。E: Exon; I: Intron.

    圖2 DjNAC3和DjNAC4蛋白的三級(jí)結(jié)構(gòu)Fig.2 Tertiary structures of DjNAC3 and DjNAC4

    2.3DjNAC的表達(dá)分析 以黃藤葉片、發(fā)育初期鉤刺、發(fā)育成熟鉤刺、發(fā)育初期纖鞭和發(fā)育成熟纖鞭的cDNA為模板,進(jìn)行實(shí)時(shí)定量PCR分析,檢測(cè)DjNAC在不同組織中的表達(dá)模式。結(jié)果顯示,DjNAC3和DjNAC4的表達(dá)模式存在著明顯的差異。DjNAC3在葉片中的表達(dá)豐度最低,遠(yuǎn)低于鉤刺和纖鞭,而纖鞭中的表達(dá)又高于鉤刺,以發(fā)育成熟的纖鞭中的表達(dá)豐度最高;DjNAC4則在發(fā)育成熟的鉤刺中表達(dá)豐度最高,發(fā)育成熟的纖鞭次之,葉片和發(fā)育初期的纖鞭中均較多,而發(fā)育初期的鉤刺中最低(圖3)。由此表明,DjNAC3和DjNAC4在黃藤葉片發(fā)育過程中可能具有不同的功能。

    2.4 SSR位點(diǎn)分析與分子標(biāo)記開發(fā) SSR搜索結(jié)果表明,在DjNAC3和DjNAC4的基因組序列中分別找到1個(gè)SSR位點(diǎn),其中前者的SSR位點(diǎn)位于內(nèi)含子區(qū)域,SSR類型為TA,重復(fù)次數(shù)為6; 后者的SSR位點(diǎn)位于第1個(gè)外顯子區(qū)域,SSR類型為GCA,重復(fù)次數(shù)為5(表3)。

    圖3 DjNAC3和DjNAC4在黃藤不同組織中的表達(dá)分析Fig.3 Expression analysis of DjNAC3 and DjNAC4 in different tissues of Daemonorops jenkinsiana1: 葉片; 2: 發(fā)育初期鉤刺; 3: 發(fā)育成熟鉤刺; 4: 發(fā)育初期纖鞭; 5: 發(fā)育成熟纖鞭。1: Leaves; 2: Early developed barbs; 3: Developed barbs; 4: Early developed cirri; 5: Developed cirri.

    表3 SSR位點(diǎn)分析Tab.3 Analysis of SSR loci

    圖4 不同藤種樣品SSR擴(kuò)增產(chǎn)物的測(cè)序分析Fig.4 Analysis of SSR sequences of PCR products amplified from six rattan samplesA: ssrNAC3-F和ssrNAC3-R擴(kuò)增產(chǎn)物序列; B: ssrNAC4-F和ssrNAC4-R擴(kuò)增產(chǎn)物序列。編號(hào)與表1中藤種相對(duì)應(yīng)。A: Sequences of products amplified by ssrNAC3-F and ssrNAC3-R; B: Sequences of products amplified by ssrNAC4-F and ssrNAC4-R. The No. represents the species in Tab.1 correspondingly.

    為進(jìn)一步驗(yàn)證SSR位點(diǎn)的可靠性,分別回收了伸長(zhǎng)鉤葉藤、老撾鉤葉藤、狹葉黃藤、黃藤屬種1、勐臘鞭藤、斑嶺省藤6個(gè)藤種的PCR產(chǎn)物進(jìn)行測(cè)序。結(jié)果表明,引物對(duì)ssrNAC3-F和ssrNAC3-R的擴(kuò)增產(chǎn)物在205~220 bp之間,不同藤種樣品中的SSR序列存在著一定的差異,既包括SSR類型變異、重復(fù)次數(shù)變化等多態(tài)性,又有SSR位點(diǎn)缺失的現(xiàn)象。如狹葉黃藤、黃藤屬種1與黃藤的SSR位點(diǎn)均為(TA)6,而伸長(zhǎng)鉤葉藤、老撾鉤葉藤和勐臘鞭藤的SSR位點(diǎn)均為(TA)4,少了2個(gè)TA重復(fù),但在斑嶺省藤中SSR位點(diǎn)卻缺失了,同時(shí)在SSR位點(diǎn)之后的30個(gè)堿基各種之間存在著不同程度的多態(tài)性(圖4A)。引物對(duì)ssrNAC4-F和ssrNAC4-R的擴(kuò)增產(chǎn)物在331~346 bp之間,其中狹葉黃藤、黃藤屬種1與黃藤的SSR位點(diǎn)均為(GCA)5,而伸長(zhǎng)鉤葉藤、老撾鉤葉藤、勐臘鞭藤和斑嶺省藤的SSR位點(diǎn)分別為(GCA)8、(GCA)9、(GCA)7和(GCA)4,不同屬間存在明顯的多態(tài)性(圖4B)。

    2.5 SSR分子標(biāo)記的不同藤種通用性分析 為驗(yàn)證SSR位點(diǎn)的通用性和多態(tài)性,利用SSR引物(ssrNAC3-F和ssrNAC3-R、ssrNAC4-F和ssrNAC4-R)(表2),以黃藤和另外20個(gè)不同藤種的基因組DNA為模板進(jìn)行擴(kuò)增。PAGE電泳結(jié)果表明,所有PCR產(chǎn)物均有明顯的目的條帶,表明2個(gè)分子標(biāo)記在不同藤種中具有通用性(圖5)。

    引物對(duì)ssrNAC3-F和ssrNAC3-R的21個(gè)樣品擴(kuò)增產(chǎn)物約在220 bp均有一明顯特異條帶(圖5A),其中鉤葉藤屬的伸長(zhǎng)鉤葉藤和老撾鉤葉藤的目的條帶大小一致,而細(xì)鉤葉藤的明顯偏大; 黃藤屬的長(zhǎng)柄黃藤、狹葉黃藤和黃藤種1的帶型一致,而黃藤的目的條帶明顯偏??; 省藤屬中除了褐鱗省藤和勐捧省藤帶型一致、種2(C. sp. 2)和種3(C. sp. 3)的帶型較一致外,其他藤種的帶型差異明顯。引物對(duì)ssrNAC4-F和ssrNAC4-R的擴(kuò)增產(chǎn)物則約在350 bp均有一明顯特異條帶(圖5 B),其中鉤葉藤屬3個(gè)藤種的目的條帶大小一致; 黃藤屬的長(zhǎng)柄黃藤、狹葉黃藤和種1(D. sp. 1)的帶型一致,而黃藤的目的條帶明顯偏小; 省藤屬柳條省藤和褐鱗省藤帶型目的條帶差異明顯,而褐鱗省藤和勐捧省藤的相一致,其余藤種的帶型差異明顯,但元江省藤(C. sp.)與種1(C. sp. 1)、種2(C. sp. 2)和種3(C. sp. 3)的基本一致。由此表明,2個(gè)SSR分子標(biāo)記在不同藤種中具有多態(tài)性,尤其是不同屬的不同藤種之間更為明顯。

    圖5 不同藤種樣品SSR引物擴(kuò)增產(chǎn)物的PAGE電泳分析Fig.5 PAGE analysis of PCR products amplified from different rattan samples using SSR primersA: ssrNAC3-F和ssrNAC3-R的擴(kuò)增產(chǎn)物; B: ssrNAC4-F和ssrNAC4-R的擴(kuò)增產(chǎn)物。M: DNA分子量標(biāo)記; 1-21: 不同藤種樣品。A. PCR products of ssrNAC3-F and ssrNAC3-R; B. PCR products of ssrNAC4-F and ssrNAC4-R. M: DNA ladder; 1-21: Different rattan samples.

    3 討論

    NAC轉(zhuǎn)錄因子對(duì)植物不同組織器官的系統(tǒng)發(fā)育具有重要作用,是植物形態(tài)建成中不可或缺的調(diào)控因子。根據(jù)蛋白結(jié)構(gòu)特征將被子植物的NAC劃分為6個(gè)大的家族,每個(gè)家族包含不同的亞家族,即家族Ⅰ(VND和NST/BRN/SMB)、Ⅱ(NAC1、CUC和ORE)、Ⅲ(NTL、NAC2和TMM)、Ⅳ(FEZ/JUB和LOV1)、Ⅴ(XND1、脅迫相關(guān)家族和NARS)和家族Ⅵ(其他) (Pereira-Santanaetal., 2015)。本研究從黃藤葉片中分離得到的DjNAC3和DjNAC4基因均屬于NAC家族的CUC亞家族。基因結(jié)構(gòu)的差異,意味著它們具有不同的功能,DjNAC3和DjNAC4在葉片、發(fā)育初期鉤刺、發(fā)育成熟鉤刺、發(fā)育初期纖鞭和發(fā)育成熟纖鞭中的基因表達(dá)模式差異進(jìn)一步證實(shí)了這一點(diǎn)。纖鞭是棕櫚藤植物為適應(yīng)森林攀援生活形成的最為重要的器官,是由葉軸頂端延伸出去形成的,或是葉鞘發(fā)育形成。黃藤的葉片屬于大型羽狀復(fù)葉(全裂),其纖鞭由葉軸頂端延伸發(fā)育而成,纖鞭上具有輪生的鉤刺,利于絞纏攀爬。DjNAC3和DjNAC4在黃藤葉片的不同部位均有表達(dá),且存在明顯差異,但其功能如何,在黃藤葉片生長(zhǎng)發(fā)育中是怎么實(shí)現(xiàn)其調(diào)控作用的,有待于進(jìn)一步深入研究。

    SSR標(biāo)記作為一種有效的分子標(biāo)記已被廣泛應(yīng)用于遺傳多樣性分析、系統(tǒng)發(fā)育、輔助分類和育種研究等。DjNAC3和DjNAC4的基因組中均具有1個(gè)SSR位點(diǎn),并證實(shí)在21個(gè)不同棕櫚藤樣品中具有通用性,且在不同屬種間具有多態(tài)性,這為進(jìn)一步研究棕櫚藤的遺傳多樣性和輔助分類提供了分子標(biāo)記。另外研究表明,SSR位點(diǎn)在基因序列中所處位置不同,其功能也存在著差異,如位于5′非翻譯區(qū)的SSR 可以調(diào)控基因的轉(zhuǎn)錄和翻譯,位于基因內(nèi)部的SSR具有較強(qiáng)的選擇壓力,可以調(diào)控基因的表達(dá),位于內(nèi)含子區(qū)域的SSR可以影響基因的轉(zhuǎn)錄,位于3′非翻譯區(qū)的SSR可以引起mRNA的延伸 (Lietal., 2004)。DjNAC3和DjNAC4的SSR位點(diǎn),分別位于內(nèi)含子區(qū)域和外顯子區(qū)域,這對(duì)二者的基因表達(dá)可能具有不同的影響,具體情況需要進(jìn)一步試驗(yàn)驗(yàn)證。

    隨著現(xiàn)代生物技術(shù)的快速發(fā)展,利用基因工程手段開發(fā)利用植物中具有調(diào)控作用的關(guān)鍵NAC 轉(zhuǎn)錄因子,實(shí)現(xiàn)植物性狀的定向改良已成為現(xiàn)實(shí)。水稻SNAC1轉(zhuǎn)入小麥(Triticumaestivum)導(dǎo)致轉(zhuǎn)基因植物生長(zhǎng)遲緩(Saadetal., 2013),毛竹PeNAC在擬南芥中過量表達(dá)促進(jìn)轉(zhuǎn)基因植株的側(cè)根發(fā)育(Wangetal., 2016),在水稻中過量表達(dá)穇子(Eleusinecoracana)的EcNAC67能夠提高水稻的耐鹽和抗旱能力(Rahmanetal., 2016),大麥(Hordeumvulgare)中過量表達(dá)HvNAC005導(dǎo)致提早成熟(Christiansenetal., 2016)。棕櫚藤植物作為重要的森林資源之一,其潛在的價(jià)值遠(yuǎn)遠(yuǎn)沒有得到開發(fā)利用,隨著對(duì)其生長(zhǎng)發(fā)育分子基礎(chǔ)研究的不斷深入,諸如NAC家族的轉(zhuǎn)錄因子等具有重大經(jīng)濟(jì)價(jià)值的功能基因?qū)⒉粩啾煌诰虺鰜恚瑢?duì)其性狀的調(diào)控將更具有操作性。因此,未來棕櫚藤的分子育種有著廣闊的前景。

    4 結(jié)論

    本研究從黃藤中克隆了2個(gè)NAC轉(zhuǎn)錄因子家族成員基因DjNAC3和DjNAC4,二者在基因結(jié)構(gòu)、編碼蛋白特征等方面具有明顯不同的分子特征,在黃藤葉片形態(tài)建成中可能具有不同的調(diào)控功能。同時(shí),DjNAC3和DjNAC4的序列中均包含1個(gè)SSR位點(diǎn),分別為(TA)6和(GCA)5,據(jù)此開發(fā)的分子標(biāo)記在不同的棕櫚藤中具有較高的通用性和多態(tài)性,可應(yīng)用于棕櫚藤輔助分類和分子輔助育種。

    高志民, 范少輝, 高 健, 等. 2006. 基于CTAB法提取毛竹基因組DNA的探討. 林業(yè)科學(xué)研究, 19(6): 725-728.

    (Gao Z M, Fan S H, Gao J,etal. 2006. Extract genomic DNA fromPhyllostachysedulisby CTAB-based method. Forest Research, 19(6):725-728. [in Chinese])

    江澤慧,王康林. 2013. 中國(guó)棕櫚藤. 北京: 科學(xué)出版社, 90-92.

    (Jiang Z H, Wang K L. 2013. Rattan in China. Beijing: Science Press, 90-92.[in Chinese])

    黎幫勇, 胡尚連, 曹 穎, 等. 2015. 毛竹NAC 轉(zhuǎn)錄因子家族生物信息學(xué)分析. 基因組學(xué)與應(yīng)用生物學(xué), 34(8): 1769-1777.

    (Li B Y, Hu S L, Cao Y,etal. 2015. Bioinformatics analysis of NAC gene family in moso bamboo. Genomics Applied Biol, 34(8): 1769-1777. [in Chinese])

    Banerjee A, Roychoudhury A. 2015. WRKY proteins: signaling and regulation of expression during abiotic stress responses. Sci World J, 2015: 807560.

    Christiansen M W, Matthewman C,Podzimska-Sroka D,etal. 2016. Barley plants over-expressing the NAC transcription factor geneHvNAC005 show stunting and delay in development combined with early senescence. J Exp Bot, 67 (17): 5259-5273.

    Gao Z M, Li X P, Li L B,etal. 2006. An effective method for total RNA isolation from bamboo. Chinese For Sci Tech, 5(3): 52-54.

    Guo Y, Gan S. 2006. AtNAP, a NAC family transcription factor, has an important role in leaf senescence. Plant J, 46(4): 601-612.

    Ha C V, Esfahani M N, Watanabe Y,etal. 2014. Genome-wide identification and expression analysis of the CaNAC family members in chickpea during development, dehydration and ABA treatments. PLoS One, 9(12): e114107.

    He X J, Mu R L, Cao W H,etal. 2005. AtNAC2, a transcription factor downstream of ethylene and auxin signaling pathways, is involved in salt stress response and lateral root development. Plant J, 44(6): 903-916.

    Kim S G, Kim S Y, Park C M. 2007. A membrane-associated NAC transcription factor regulates salt-responsive flowering via FLOWERING LOCUST inArabidopsis. Planta, 6(3): 647-654.

    Kim S G, Lee A K, Yoon H K,etal. 2008. A membrane-bound NAC transcription factor NTL8 regulates gibberellic acid-mediated salt signaling inArabidopsisseed germination. Plant J, 55(1):77-88.

    Kou X, Watkins C B, Gan S S. 2012.ArabidopsisAtNAP regulates fruit senescence. J Exp Bot, 63(17): 6139-6147.

    Le D T, Nishiyama R, Watanabe Y,etal. 2011. Genome-wide survey and expression analysis of the plant-specific NAC transcription factor family in soybean during development and dehydration stress. DNA Res, 18(4): 263-276.

    Li Y C, Korol A B, Fahima T,etal. 2004. Microsatellites within genes: structure, function, and evolution. Mol Biol Evol, 21(6): 991-1007.

    Lin R M, Zhao W S, Meng X B,etal. 2007. Rice geneOsNAC19 encodes a novel NAC-domain transcription factor and responds to infection byMagnaporthegrisea. Plant Sci, 172(1): 120-130.

    Liu T K, Song X M, Duan W K,etal. 2014. Genome-wide analysis and expression patterns of NAC transcription factor family under different developmental stages and abiotic stresses in Chinese cabbage. Plant Mol Biol Report, 32(5): 1041-1056.

    Livak K J, Schmittgen T D. 2001. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCtmethod. Methods, 25: 402-408.

    Moore M J, Query C C, Sharp P A. 1993. Splicing of precursors to mRNA by the spliceosome∥ Gesteland R F, Atkins J F. The RNA world. Cold Spring Harbor: Cold Spring Harbor Laboratory Press, 303-357.

    Nakashima K, Takasaki H, Mizoi J,etal. 2012. NAC transcription factors in plant abiotic stress responses. Biochim Biophys Acta, 1819(2): 97-103.

    Nikovics K, Blein T, Peaucelle A,etal. 2006. The balance between theMIR164AandCUC2 genes controls leaf margin serration inArabidopsis. Plant Cell, 18(11): 2929-2945.

    Niu F, Wang C, Yan J,etal. 2016. Functional characterization of NAC55 transcription factor from oilseed rape (BrassicanapusL.) as a novel transcriptional activator modulating reactive oxygen species accumulation and cell death. Plant Mol Biol, 92(1/2): 89-104.

    Nuruzzaman M, Manimekalai R, Sharoni A M,etal. 2010. Genome-wide analysis of NAC transcription factor family in rice. Gene, 465(1/2): 30-44.

    Nuruzzaman M, Sharoni A M, Kikuchi S. 2013. Roles of NAC transcription factors in the regulation of biotic and abiotic stress responses in plants. Front Microbiol, 4: 248.

    Olsen A N, Ernst H A, Leggio L L,etal. 2005. NAC transcription factors: structurally distinct, functionally diverse. Trends Plant Sci, 10(2): 79-87.

    Pereira-Santana A, Alcaraz L D, Castao E,etal. 2015. Comparative genomics of NAC transcriptional factors in angiosperms: Implications for the adaptation and diversification of flowering plants. PLoS One, 10(11): e0141866.

    Rahman H, Ramanathan V, Nallathambi J,etal. 2016. Over-expression of a NAC 67 transcription factor from finger millet (EleusinecoracanaL.) confers tolerance against salinity and drought stress in rice. BMC Biotechnol, Suppl 1: 35.

    Saad A S, Li X, Li H P,etal. 2013. A rice stress-responsive NAC gene enhances tolerance of transgenic wheat to drought and salt stresses. Plant Sci, 203/204: 33-40.

    Seo P J, Kim M J, Park J Y,etal. 2010. Cold activation of a plasma membrane-tethered NAC transcription factor induces a pathogen resistance response inArabidopsis. Plant J, 61(4): 661-671.

    Shiriga K, Sharma R, Kumar K,etal. 2014. Genome-wide identification and expression pattern of drought-responsive members of the NAC family in maize. Meta Gene, 2: 407-417.

    Souer E, van Houwelingen A, Kloos D,etal. 1996. Thenoapicalmeristemgene of petunia is required for pattern formation in embryos and flowers and is expressed at meristem and primordia boundaries. Cell, 85(2): 159-170.

    Su H, Zhang S, Yin Y,etal. 2015. Genome-wide analysis of NAM-ATAF1,2-CUC2 transcription factor family inSolanumlycopersicum. J Plant Biochem Biotechnol, 24(2): 176-183.

    Temnykh S, DeClerck G, Lukashova A,etal. 2001. Computational and experimental analysis of microsatellites in rice (OryzasativaL.): frequency, length variation, transposon associations, and genetic marker potential. Genome Res, 11(8): 1441-1452.

    Tran L S, Nishiyama R,Yamaguchi-Shinozaki K,etal. 2010. Potential utilization of NAC transcription factors to enhance abiotic stress tolerance in plants by biotechnological approach. GM Crops, 1(1): 32-39.

    Uauy C, Distelfeld A, Fahima T,etal. 2006. A NAC gene regulating senescence improves grain protein, zinc, and iron content in wheat. Science, 314(5803): 1298-1301.

    Wang L, Zhao H, Chen D,etal. 2016. Characterization and primary functional analysis of a bambooNACgene targeted by miR164b. Plant Cell Rep, (6):1371-1383.

    Zhu G, Chen G, Zhu J,etal. 2015. Molecular characterization and expression profiling of NAC transcription factors inBrachypodiumdistachyonL. PLoS One, 10(10): e0139794.

    Zimmermann R, Werr W. 2005. Pattern formation in the monocot embryo as revealed by NAM and CUC3 orthologues fromZeamaysL. Plant Mol Biol, 58(5): 669-685.

    (責(zé)任編輯 徐 紅)

    MolecularCharacteristicsandSSRMarkerDevelopmentofTwoNACGenesfromDaemonoropsjenkinsiana

    Sun Huayu1Li Lichao1Zhao Hansheng1Yang Yihong1,2Wang Sining1Gao Zhimin1

    (1.KeyOpenLaboratoryontheScienceandTechnologyofBambooandRattanofStateForestryAdministrationInternationalCenterforBambooandRattanBeijing100102; 2.CollegeofHorticulture,HebeiAgriculturalUniversityBaoding071001)

    S718.46

    A

    1001-7488(2017)08-0132-09

    10.11707/j.1001-7488.20170815

    2016-09-12;

    2016-10-24。

    “十二五”農(nóng)村領(lǐng)域國(guó)家科技計(jì)劃項(xiàng)目“竹藤種質(zhì)資源創(chuàng)新利用研究”第一課題“竹藤優(yōu)異種質(zhì)創(chuàng)制創(chuàng)新與種苗培育標(biāo)準(zhǔn)化示范”(2015BAD04B01)。

    *高志民為通訊作者。

    猜你喜歡
    基因組位點(diǎn)引物
    DNA引物合成起始的分子基礎(chǔ)
    高中生物學(xué)PCR技術(shù)中“引物”相關(guān)問題歸類分析
    鎳基單晶高溫合金多組元置換的第一性原理研究
    上海金屬(2021年6期)2021-12-02 10:47:20
    SSR-based hybrid identification, genetic analyses and fingerprint development of hybridization progenies from sympodial bamboo (Bambusoideae, Poaceae)
    牛參考基因組中發(fā)現(xiàn)被忽視基因
    CLOCK基因rs4580704多態(tài)性位點(diǎn)與2型糖尿病和睡眠質(zhì)量的相關(guān)性
    二項(xiàng)式通項(xiàng)公式在遺傳學(xué)計(jì)算中的運(yùn)用*
    火炬松SSR-PCR反應(yīng)體系的建立及引物篩選
    基因組DNA甲基化及組蛋白甲基化
    遺傳(2014年3期)2014-02-28 20:58:49
    有趣的植物基因組
    a级一级毛片免费在线观看| 人妻夜夜爽99麻豆av| h视频一区二区三区| 99热6这里只有精品| 午夜福利视频精品| 国产精品一及| 在线天堂最新版资源| 免费黄频网站在线观看国产| 国产精品99久久99久久久不卡 | 国产成人91sexporn| 亚洲av日韩在线播放| 久久亚洲国产成人精品v| 晚上一个人看的免费电影| 久久青草综合色| 爱豆传媒免费全集在线观看| av播播在线观看一区| 国产免费福利视频在线观看| 国产精品不卡视频一区二区| 成人二区视频| 少妇猛男粗大的猛烈进出视频| 久久人人爽人人爽人人片va| 又爽又黄a免费视频| 亚洲欧美中文字幕日韩二区| 亚洲精品久久午夜乱码| 亚州av有码| av播播在线观看一区| 欧美三级亚洲精品| 国产伦在线观看视频一区| 老女人水多毛片| 国产69精品久久久久777片| 99久久综合免费| 精品久久久久久久久av| 亚洲成人中文字幕在线播放| 欧美xxⅹ黑人| 在线播放无遮挡| 极品少妇高潮喷水抽搐| 国产美女午夜福利| 高清不卡的av网站| 国产亚洲午夜精品一区二区久久| 免费播放大片免费观看视频在线观看| 日本黄大片高清| 日韩,欧美,国产一区二区三区| 尤物成人国产欧美一区二区三区| 免费黄色在线免费观看| 国模一区二区三区四区视频| 久久久久人妻精品一区果冻| 国产伦精品一区二区三区视频9| 亚洲精品久久久久久婷婷小说| 国产成人精品福利久久| 成人一区二区视频在线观看| 中文字幕免费在线视频6| 大香蕉97超碰在线| 欧美一区二区亚洲| 大码成人一级视频| 男女下面进入的视频免费午夜| 亚洲精品视频女| 精品国产露脸久久av麻豆| 欧美+日韩+精品| 在线观看人妻少妇| 中文字幕免费在线视频6| 亚洲av综合色区一区| 久久久色成人| 免费久久久久久久精品成人欧美视频 | 亚洲真实伦在线观看| 亚洲经典国产精华液单| 国产精品国产av在线观看| 老司机影院毛片| 边亲边吃奶的免费视频| 自拍偷自拍亚洲精品老妇| 日韩不卡一区二区三区视频在线| 国产精品人妻久久久影院| 欧美一区二区亚洲| 91精品伊人久久大香线蕉| 日韩免费高清中文字幕av| 99久久精品国产国产毛片| av不卡在线播放| 在现免费观看毛片| 秋霞在线观看毛片| 熟妇人妻不卡中文字幕| 又爽又黄a免费视频| 涩涩av久久男人的天堂| 人人妻人人看人人澡| 人人妻人人爽人人添夜夜欢视频 | 国产精品久久久久成人av| 麻豆成人午夜福利视频| 婷婷色av中文字幕| 国产 一区 欧美 日韩| 国产淫语在线视频| 亚洲国产精品999| 成人一区二区视频在线观看| 日本欧美视频一区| 直男gayav资源| 国产深夜福利视频在线观看| 欧美成人午夜免费资源| 99热这里只有精品一区| 国产乱人偷精品视频| 日韩视频在线欧美| 欧美老熟妇乱子伦牲交| 麻豆成人av视频| 精品少妇黑人巨大在线播放| 人体艺术视频欧美日本| 国产熟女欧美一区二区| 免费观看在线日韩| 最近中文字幕2019免费版| av在线蜜桃| 久久久久久九九精品二区国产| 免费高清在线观看视频在线观看| 欧美亚洲 丝袜 人妻 在线| 日韩人妻高清精品专区| 精品国产一区二区三区久久久樱花 | 高清在线视频一区二区三区| 中文字幕免费在线视频6| 久久精品国产亚洲网站| 成人二区视频| 久久久久久人妻| 国产精品麻豆人妻色哟哟久久| 观看免费一级毛片| 久热这里只有精品99| a级一级毛片免费在线观看| 成人高潮视频无遮挡免费网站| 欧美xxⅹ黑人| 欧美日韩视频精品一区| 亚洲美女视频黄频| 九色成人免费人妻av| 久久久久人妻精品一区果冻| 欧美日韩视频高清一区二区三区二| 最近中文字幕高清免费大全6| 欧美成人一区二区免费高清观看| 毛片一级片免费看久久久久| 好男人视频免费观看在线| 精品久久久精品久久久| 成年人午夜在线观看视频| 夜夜骑夜夜射夜夜干| 亚洲欧美一区二区三区国产| 女的被弄到高潮叫床怎么办| 乱码一卡2卡4卡精品| 中文字幕久久专区| 欧美高清成人免费视频www| 九草在线视频观看| 亚洲精品一二三| 亚洲欧美成人精品一区二区| 在线播放无遮挡| 精品亚洲成国产av| 人妻一区二区av| 亚洲经典国产精华液单| 1000部很黄的大片| 亚洲成人中文字幕在线播放| 亚洲综合色惰| 国产国拍精品亚洲av在线观看| 丰满迷人的少妇在线观看| 国产亚洲欧美精品永久| 十分钟在线观看高清视频www | .国产精品久久| 91久久精品电影网| 久久99精品国语久久久| 国产精品国产三级专区第一集| 日韩精品有码人妻一区| 在线观看一区二区三区| 一区在线观看完整版| 亚洲综合色惰| 97在线视频观看| 免费观看无遮挡的男女| 春色校园在线视频观看| 99热这里只有是精品在线观看| 国产在线一区二区三区精| 赤兔流量卡办理| h日本视频在线播放| 欧美xxⅹ黑人| 亚洲一区二区三区欧美精品| 在线观看免费日韩欧美大片 | 建设人人有责人人尽责人人享有的 | 男人和女人高潮做爰伦理| 久久久久网色| 妹子高潮喷水视频| 国产精品麻豆人妻色哟哟久久| 99re6热这里在线精品视频| 免费观看在线日韩| 一区在线观看完整版| 777米奇影视久久| 日韩成人av中文字幕在线观看| 亚洲av免费高清在线观看| 日本黄色片子视频| 新久久久久国产一级毛片| 国产av一区二区精品久久 | 在线观看一区二区三区| 最近最新中文字幕免费大全7| 青春草亚洲视频在线观看| 91在线精品国自产拍蜜月| 亚洲va在线va天堂va国产| 国产精品一二三区在线看| 又爽又黄a免费视频| 免费黄频网站在线观看国产| 在线观看一区二区三区激情| 在线精品无人区一区二区三 | 国产免费视频播放在线视频| 欧美3d第一页| 国产 一区 欧美 日韩| 中文资源天堂在线| 夜夜看夜夜爽夜夜摸| 一区二区三区精品91| 国产亚洲91精品色在线| 日本av手机在线免费观看| 亚洲欧美中文字幕日韩二区| 亚洲精品456在线播放app| 少妇精品久久久久久久| 一个人看视频在线观看www免费| av在线蜜桃| 国产成人91sexporn| 偷拍熟女少妇极品色| 夜夜骑夜夜射夜夜干| 特大巨黑吊av在线直播| 国产69精品久久久久777片| 尤物成人国产欧美一区二区三区| 日日摸夜夜添夜夜添av毛片| 九九在线视频观看精品| 肉色欧美久久久久久久蜜桃| 日本与韩国留学比较| 色婷婷av一区二区三区视频| 噜噜噜噜噜久久久久久91| 能在线免费看毛片的网站| 麻豆精品久久久久久蜜桃| 国产精品麻豆人妻色哟哟久久| 黄色一级大片看看| 久久99精品国语久久久| 国产男女内射视频| 精品久久久久久久末码| 久久人人爽人人片av| 亚洲av成人精品一区久久| 国产一区二区三区综合在线观看 | 自拍偷自拍亚洲精品老妇| 九色成人免费人妻av| 简卡轻食公司| 欧美 日韩 精品 国产| 国产色婷婷99| 不卡视频在线观看欧美| 中文天堂在线官网| 少妇人妻久久综合中文| 久久久久人妻精品一区果冻| 国产精品久久久久久精品电影小说 | 久热这里只有精品99| 久久国产精品大桥未久av | 欧美精品一区二区大全| 久久精品国产自在天天线| 五月玫瑰六月丁香| 激情五月婷婷亚洲| 国产亚洲最大av| 热99国产精品久久久久久7| a级毛色黄片| 观看免费一级毛片| 亚洲国产日韩一区二区| av网站免费在线观看视频| 亚洲国产av新网站| 一级a做视频免费观看| 大片免费播放器 马上看| 在线免费观看不下载黄p国产| 免费黄频网站在线观看国产| 国产精品.久久久| 中文天堂在线官网| 国内揄拍国产精品人妻在线| 亚洲av.av天堂| 视频中文字幕在线观看| 国产片特级美女逼逼视频| 在线亚洲精品国产二区图片欧美 | 精华霜和精华液先用哪个| 国产成人免费观看mmmm| 一级片'在线观看视频| 精品久久久久久久末码| 亚洲国产成人一精品久久久| 十八禁网站网址无遮挡 | 国产熟女欧美一区二区| 午夜福利在线在线| 国产欧美另类精品又又久久亚洲欧美| 精品少妇久久久久久888优播| 男人和女人高潮做爰伦理| 久久久久久久国产电影| 人人妻人人看人人澡| 永久网站在线| 搡女人真爽免费视频火全软件| 少妇裸体淫交视频免费看高清| 天堂俺去俺来也www色官网| a级一级毛片免费在线观看| 高清黄色对白视频在线免费看 | 国产精品无大码| 一级毛片电影观看| 一级黄片播放器| 午夜免费男女啪啪视频观看| 亚洲电影在线观看av| 99热国产这里只有精品6| 少妇裸体淫交视频免费看高清| 国产乱来视频区| av专区在线播放| 1000部很黄的大片| 亚洲精品乱久久久久久| 嫩草影院入口| 国产一级毛片在线| 人妻一区二区av| 哪个播放器可以免费观看大片| 青青草视频在线视频观看| 人妻系列 视频| 97超碰精品成人国产| 男女啪啪激烈高潮av片| 国产成人免费观看mmmm| 国国产精品蜜臀av免费| av不卡在线播放| videos熟女内射| 免费看光身美女| 日本与韩国留学比较| 狂野欧美激情性xxxx在线观看| 一本一本综合久久| 亚洲图色成人| 国产成人freesex在线| 欧美精品人与动牲交sv欧美| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品国产av成人精品| 丝瓜视频免费看黄片| 久久久久久久久久成人| 国产一区二区三区av在线| 青青草视频在线视频观看| 嘟嘟电影网在线观看| 91狼人影院| 在线看a的网站| 日韩不卡一区二区三区视频在线| 国产精品女同一区二区软件| 成人二区视频| 不卡视频在线观看欧美| av国产免费在线观看| 一区二区三区乱码不卡18| 少妇人妻精品综合一区二区| 一级毛片aaaaaa免费看小| 中文天堂在线官网| 欧美性感艳星| 美女视频免费永久观看网站| 大香蕉久久网| www.色视频.com| 精品久久久久久久末码| 欧美日韩亚洲高清精品| 一个人看视频在线观看www免费| 日韩中字成人| 天天躁日日操中文字幕| 精品亚洲成a人片在线观看 | 九色成人免费人妻av| 少妇人妻精品综合一区二区| 日韩成人伦理影院| 这个男人来自地球电影免费观看 | 啦啦啦视频在线资源免费观看| 国产精品久久久久久久久免| 国产在线视频一区二区| 国产一区二区在线观看日韩| 最近手机中文字幕大全| 国产伦精品一区二区三区四那| 蜜臀久久99精品久久宅男| 色视频在线一区二区三区| 亚洲美女搞黄在线观看| 高清黄色对白视频在线免费看 | 91久久精品国产一区二区三区| 久久人人爽人人片av| 高清毛片免费看| 亚洲国产精品999| 在线看a的网站| 少妇人妻 视频| 嫩草影院入口| 人妻 亚洲 视频| 蜜臀久久99精品久久宅男| 国产成人免费无遮挡视频| 亚洲国产精品国产精品| 插逼视频在线观看| 97在线视频观看| 亚洲av综合色区一区| 国产黄片美女视频| 美女主播在线视频| 啦啦啦中文免费视频观看日本| 99久久人妻综合| 日本黄大片高清| 激情五月婷婷亚洲| 天堂俺去俺来也www色官网| 不卡视频在线观看欧美| 精品久久久久久电影网| 免费观看的影片在线观看| 国产av一区二区精品久久 | 又粗又硬又长又爽又黄的视频| 男人狂女人下面高潮的视频| 亚洲人成网站在线播| av天堂中文字幕网| 一本色道久久久久久精品综合| 亚洲国产色片| 欧美一级a爱片免费观看看| 男人舔奶头视频| 亚洲,一卡二卡三卡| 国产 一区 欧美 日韩| 日韩一本色道免费dvd| 最近中文字幕高清免费大全6| 特大巨黑吊av在线直播| 熟女电影av网| 丰满乱子伦码专区| 欧美日韩亚洲高清精品| 美女脱内裤让男人舔精品视频| 日本猛色少妇xxxxx猛交久久| 国产亚洲一区二区精品| 伊人久久精品亚洲午夜| 久久99热这里只有精品18| 午夜福利影视在线免费观看| 久久这里有精品视频免费| 亚洲熟女精品中文字幕| videossex国产| a级一级毛片免费在线观看| 国产精品无大码| videossex国产| h视频一区二区三区| 亚洲成人中文字幕在线播放| 美女cb高潮喷水在线观看| 亚洲人成网站在线观看播放| 久久精品国产自在天天线| 插阴视频在线观看视频| 大香蕉久久网| 视频中文字幕在线观看| kizo精华| 国产亚洲欧美精品永久| 亚洲成人一二三区av| 寂寞人妻少妇视频99o| 99九九线精品视频在线观看视频| 久久精品国产a三级三级三级| 国产黄色免费在线视频| 国产成人a∨麻豆精品| 小蜜桃在线观看免费完整版高清| 精品国产乱码久久久久久小说| 简卡轻食公司| 国产精品.久久久| 91在线精品国自产拍蜜月| 99国产精品免费福利视频| 在线看a的网站| 国产高清有码在线观看视频| 成人一区二区视频在线观看| 亚洲国产av新网站| 久久99热6这里只有精品| 国产精品久久久久久久电影| 丰满迷人的少妇在线观看| 97在线人人人人妻| 精品人妻视频免费看| 大香蕉97超碰在线| 人人妻人人爽人人添夜夜欢视频 | 人妻 亚洲 视频| 亚洲av成人精品一二三区| 久久97久久精品| 夜夜爽夜夜爽视频| 久久久成人免费电影| 国产深夜福利视频在线观看| 亚洲精品,欧美精品| 日韩三级伦理在线观看| 最近2019中文字幕mv第一页| av在线蜜桃| 精品久久久精品久久久| 色综合色国产| 深爱激情五月婷婷| 亚洲av中文字字幕乱码综合| 精品午夜福利在线看| 99热这里只有精品一区| 少妇人妻一区二区三区视频| 人妻系列 视频| 男女啪啪激烈高潮av片| 麻豆成人午夜福利视频| 久久人妻熟女aⅴ| 成人无遮挡网站| 日韩国内少妇激情av| 偷拍熟女少妇极品色| 人人妻人人看人人澡| 国产亚洲午夜精品一区二区久久| 国产精品麻豆人妻色哟哟久久| 国产欧美另类精品又又久久亚洲欧美| 又爽又黄a免费视频| 草草在线视频免费看| 亚洲精品456在线播放app| av免费观看日本| 精品亚洲成a人片在线观看 | 国产成人午夜福利电影在线观看| 欧美zozozo另类| 黄色配什么色好看| 国精品久久久久久国模美| 亚洲在久久综合| 黑人猛操日本美女一级片| 夜夜爽夜夜爽视频| 五月玫瑰六月丁香| 少妇精品久久久久久久| 我的女老师完整版在线观看| 妹子高潮喷水视频| 午夜福利在线在线| 亚洲av免费高清在线观看| 久久久久久久久久久免费av| 久久综合国产亚洲精品| 小蜜桃在线观看免费完整版高清| 各种免费的搞黄视频| 尤物成人国产欧美一区二区三区| 免费播放大片免费观看视频在线观看| 美女视频免费永久观看网站| av播播在线观看一区| 99热这里只有是精品在线观看| 婷婷色麻豆天堂久久| 51国产日韩欧美| 中文天堂在线官网| 在线免费十八禁| 成年美女黄网站色视频大全免费 | 国产精品无大码| 婷婷色av中文字幕| 亚洲熟女精品中文字幕| 一级av片app| 亚洲欧美日韩无卡精品| 一级av片app| 亚洲成人中文字幕在线播放| 久久久久久久大尺度免费视频| av免费观看日本| 欧美 日韩 精品 国产| 日韩伦理黄色片| 蜜桃久久精品国产亚洲av| 精品国产三级普通话版| 日日啪夜夜爽| 中国三级夫妇交换| 久久国产精品男人的天堂亚洲 | 联通29元200g的流量卡| 国产黄片视频在线免费观看| 全区人妻精品视频| 色婷婷久久久亚洲欧美| 91aial.com中文字幕在线观看| 老女人水多毛片| 一级二级三级毛片免费看| 美女脱内裤让男人舔精品视频| 99热6这里只有精品| 菩萨蛮人人尽说江南好唐韦庄| 少妇的逼水好多| 大片免费播放器 马上看| 91午夜精品亚洲一区二区三区| 蜜桃久久精品国产亚洲av| 国产亚洲午夜精品一区二区久久| 女性生殖器流出的白浆| 欧美xxxx性猛交bbbb| 草草在线视频免费看| 少妇的逼水好多| 国产免费一区二区三区四区乱码| 亚洲精品久久久久久婷婷小说| 精品亚洲乱码少妇综合久久| 久久久久人妻精品一区果冻| 全区人妻精品视频| 男人狂女人下面高潮的视频| 一级片'在线观看视频| 看免费成人av毛片| 简卡轻食公司| 国产精品av视频在线免费观看| 国产大屁股一区二区在线视频| 久久人人爽人人爽人人片va| 欧美变态另类bdsm刘玥| 久久这里有精品视频免费| 国产 一区精品| 成年av动漫网址| 爱豆传媒免费全集在线观看| 免费人成在线观看视频色| 日韩制服骚丝袜av| 又黄又爽又刺激的免费视频.| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲在久久综合| 国产高清有码在线观看视频| 少妇高潮的动态图| 免费观看无遮挡的男女| 国产伦在线观看视频一区| 黑人猛操日本美女一级片| 中文资源天堂在线| 美女xxoo啪啪120秒动态图| 内射极品少妇av片p| 毛片一级片免费看久久久久| 国产免费一级a男人的天堂| 色视频www国产| 免费播放大片免费观看视频在线观看| 国产成人精品一,二区| av一本久久久久| 亚洲最大成人中文| 国产在视频线精品| 久久精品人妻少妇| 中文字幕亚洲精品专区| 五月玫瑰六月丁香| 亚洲无线观看免费| videossex国产| 久久人妻熟女aⅴ| 99热全是精品| 国产伦精品一区二区三区视频9| 成人综合一区亚洲| 久久影院123| 性高湖久久久久久久久免费观看| av国产久精品久网站免费入址| 国产精品久久久久久久电影| 一本一本综合久久| 精品熟女少妇av免费看| 亚洲婷婷狠狠爱综合网| 午夜福利在线观看免费完整高清在| 少妇 在线观看| 日韩av不卡免费在线播放| 成人二区视频| 99热网站在线观看| 国产在线男女| 毛片一级片免费看久久久久| 肉色欧美久久久久久久蜜桃| 丝袜脚勾引网站| 国产亚洲5aaaaa淫片| 性色av一级| 少妇被粗大猛烈的视频| 内地一区二区视频在线| 婷婷色综合大香蕉| 久久久久久久久久久丰满| 午夜激情久久久久久久| 高清视频免费观看一区二区| 免费大片18禁| 国产精品.久久久| 国产成人a∨麻豆精品| 一级毛片黄色毛片免费观看视频| 男人和女人高潮做爰伦理| 免费大片黄手机在线观看| 日韩欧美精品免费久久| 国产精品久久久久久久电影| 我要看黄色一级片免费的|