• <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级毛片在线播放| 国产熟女欧美一区二区| 久久av网站| 久久久久久久久久成人| 在线天堂最新版资源| 久久99一区二区三区| 亚洲精品视频女| 人成视频在线观看免费观看| 久久热精品热| 久久鲁丝午夜福利片| 国产不卡av网站在线观看| 丝瓜视频免费看黄片| 人妻系列 视频| 亚洲美女搞黄在线观看| av在线app专区| 欧美日韩亚洲高清精品| 少妇精品久久久久久久| 麻豆成人av视频| 成年美女黄网站色视频大全免费 | 免费高清在线观看日韩| 免费人成在线观看视频色| 狠狠精品人妻久久久久久综合| 国产老妇伦熟女老妇高清| 精品一区在线观看国产| 亚洲av.av天堂| 26uuu在线亚洲综合色| 国产精品国产三级专区第一集| 最近中文字幕2019免费版| 国产不卡av网站在线观看| 国产精品国产av在线观看| 一区二区三区免费毛片| 天美传媒精品一区二区| 亚洲欧美日韩卡通动漫| 日韩欧美精品免费久久| 好男人视频免费观看在线| 男女边吃奶边做爰视频| 高清av免费在线| 一级a做视频免费观看| 午夜激情久久久久久久| 亚洲精品国产色婷婷电影| xxx大片免费视频| 亚洲国产欧美在线一区| 最近最新中文字幕免费大全7| 五月玫瑰六月丁香| 天堂8中文在线网| 成年人午夜在线观看视频| 中文字幕免费在线视频6| 久久久久久久久久久丰满| 99久久人妻综合| 国产日韩欧美亚洲二区| 亚洲av成人精品一区久久| 少妇精品久久久久久久| 亚洲丝袜综合中文字幕| 另类亚洲欧美激情| 伦理电影大哥的女人| 女人久久www免费人成看片| 18+在线观看网站| 国产在线视频一区二区| 国产熟女欧美一区二区| 欧美人与性动交α欧美精品济南到 | 丰满饥渴人妻一区二区三| 性色avwww在线观看| 视频中文字幕在线观看| 美女内射精品一级片tv| 日韩亚洲欧美综合| 一级黄片播放器| 亚洲精品一区蜜桃| 日本91视频免费播放| 日本黄色日本黄色录像| 国产av精品麻豆| 九九在线视频观看精品| 九九在线视频观看精品| 亚洲精品一二三| 久热久热在线精品观看| 亚洲精品av麻豆狂野| 精品一区二区免费观看| 亚洲欧美成人综合另类久久久| 欧美精品一区二区免费开放| 欧美日韩视频精品一区| 亚洲精品日韩av片在线观看| 91国产中文字幕| 亚洲精品自拍成人| 欧美日韩综合久久久久久| 黄色视频在线播放观看不卡| 国产片特级美女逼逼视频| 老司机影院成人| 老司机影院成人| 黄色视频在线播放观看不卡| 中文天堂在线官网| 少妇高潮的动态图| 波野结衣二区三区在线| 免费观看在线日韩| 欧美激情极品国产一区二区三区 | 国产av码专区亚洲av| 如日韩欧美国产精品一区二区三区 | 一级a做视频免费观看| av卡一久久| 热re99久久精品国产66热6| 日韩精品免费视频一区二区三区 | 少妇高潮的动态图| 国产午夜精品久久久久久一区二区三区| 夫妻午夜视频| 亚洲国产色片| 日本欧美视频一区| 成年美女黄网站色视频大全免费 | 高清在线视频一区二区三区| 性高湖久久久久久久久免费观看| 美女主播在线视频| 一区二区三区四区激情视频| 国产乱来视频区| 18+在线观看网站| 亚洲精品视频女| 国国产精品蜜臀av免费| 一本一本综合久久| 午夜福利视频在线观看免费| 欧美最新免费一区二区三区| 国产极品粉嫩免费观看在线 | 午夜福利视频在线观看免费| 免费av不卡在线播放| 日韩不卡一区二区三区视频在线| 久久精品久久精品一区二区三区| 国产乱来视频区| 国产毛片在线视频| 亚洲精品,欧美精品| 欧美另类一区| 狂野欧美白嫩少妇大欣赏| 成人午夜精彩视频在线观看| 寂寞人妻少妇视频99o| 最近最新中文字幕免费大全7| 插阴视频在线观看视频| 久久精品国产亚洲av天美| 亚洲无线观看免费| 国产精品无大码| 中文字幕最新亚洲高清| 一级黄片播放器| 在线播放无遮挡| 婷婷色麻豆天堂久久| 国产亚洲欧美精品永久| 久久青草综合色| 欧美亚洲日本最大视频资源| 国产探花极品一区二区| 亚洲精品久久久久久婷婷小说| 我的女老师完整版在线观看| 中文字幕av电影在线播放| 一边摸一边做爽爽视频免费| 少妇精品久久久久久久| 少妇高潮的动态图| 日韩中字成人| xxxhd国产人妻xxx| 天美传媒精品一区二区| 国产亚洲午夜精品一区二区久久| 亚洲综合精品二区| 香蕉精品网在线| 午夜免费观看性视频| 午夜福利在线观看免费完整高清在| 欧美精品人与动牲交sv欧美| 国产一区二区在线观看日韩| 大码成人一级视频| 黄色毛片三级朝国网站| 黄色视频在线播放观看不卡| 亚洲人成网站在线播| 国产淫语在线视频| 国产av码专区亚洲av| 日韩av在线免费看完整版不卡| 一区二区三区精品91| 国产 一区精品| 国产国语露脸激情在线看| 桃花免费在线播放| 色哟哟·www| 久久精品久久精品一区二区三区| 好男人视频免费观看在线| 亚洲欧美一区二区三区国产| 天堂俺去俺来也www色官网| 欧美亚洲日本最大视频资源| 精品人妻一区二区三区麻豆| 男女无遮挡免费网站观看| 天天影视国产精品| 内地一区二区视频在线| 国产色爽女视频免费观看| 午夜免费男女啪啪视频观看| 日韩熟女老妇一区二区性免费视频| 日本免费在线观看一区| 2022亚洲国产成人精品| 少妇人妻久久综合中文| 男女高潮啪啪啪动态图| 最近手机中文字幕大全| 丰满饥渴人妻一区二区三| 亚洲精品国产av蜜桃| 在线观看一区二区三区激情| 在线精品无人区一区二区三| 高清午夜精品一区二区三区| 欧美精品一区二区大全| 亚洲美女黄色视频免费看| 99热网站在线观看| a级毛片免费高清观看在线播放| 国产黄片视频在线免费观看| 人妻人人澡人人爽人人| 熟女人妻精品中文字幕| 欧美成人精品欧美一级黄| 丝袜喷水一区| 亚洲人成网站在线观看播放| 91久久精品电影网| av播播在线观看一区| 有码 亚洲区| 美女cb高潮喷水在线观看| 18禁动态无遮挡网站| 伦精品一区二区三区| 91午夜精品亚洲一区二区三区| 亚洲av.av天堂| 午夜福利视频精品| 久久久精品区二区三区| 亚洲三级黄色毛片| 精品国产露脸久久av麻豆| 69精品国产乱码久久久| 久久午夜综合久久蜜桃| 热99国产精品久久久久久7| 18禁观看日本| 中国美白少妇内射xxxbb| 成年美女黄网站色视频大全免费 | 卡戴珊不雅视频在线播放| 国产精品久久久久久精品古装| 丁香六月天网| 色5月婷婷丁香| 男女免费视频国产| 最黄视频免费看| 激情五月婷婷亚洲| 黑人欧美特级aaaaaa片| 一级黄片播放器| 成人综合一区亚洲| 精品国产露脸久久av麻豆| 精品国产乱码久久久久久小说| 有码 亚洲区| 最近手机中文字幕大全| 成人二区视频| 在线观看免费视频网站a站| 国产av精品麻豆| 成人国产av品久久久| 老女人水多毛片| 中国美白少妇内射xxxbb| 最近中文字幕2019免费版| 亚洲精品久久久久久婷婷小说| 最近最新中文字幕免费大全7| 91在线精品国自产拍蜜月| 飞空精品影院首页| 另类亚洲欧美激情| 2022亚洲国产成人精品| 一级二级三级毛片免费看| 国产69精品久久久久777片| 欧美精品高潮呻吟av久久| 国产有黄有色有爽视频| 日韩三级伦理在线观看| 国模一区二区三区四区视频| 国产无遮挡羞羞视频在线观看| 少妇的逼水好多| 男女啪啪激烈高潮av片| 精品一区在线观看国产| 丝袜脚勾引网站| 久久精品国产鲁丝片午夜精品| 亚洲精品一二三| videosex国产| 街头女战士在线观看网站| 国产精品三级大全| 亚洲国产精品成人久久小说| 涩涩av久久男人的天堂| 亚洲欧美成人精品一区二区| 伊人亚洲综合成人网| 在线观看国产h片| 久久久精品区二区三区| 久久亚洲国产成人精品v| 亚洲欧美精品自产自拍| 亚洲av.av天堂| 久久午夜福利片| 高清毛片免费看| 看十八女毛片水多多多| 一区二区av电影网| 久久 成人 亚洲| 久久久久久久国产电影| 亚洲国产日韩一区二区| 69精品国产乱码久久久| 王馨瑶露胸无遮挡在线观看| 久久久久久久久久久丰满| 成人免费观看视频高清| 日韩精品有码人妻一区| 国产成人av激情在线播放 | 久久婷婷青草| 成人毛片a级毛片在线播放| 丰满乱子伦码专区| a 毛片基地| 精品久久久噜噜| 纵有疾风起免费观看全集完整版| 免费久久久久久久精品成人欧美视频 | 国产亚洲最大av| 99九九线精品视频在线观看视频| 校园人妻丝袜中文字幕| 色婷婷久久久亚洲欧美| 99热全是精品| 免费少妇av软件| 国产成人精品久久久久久| 各种免费的搞黄视频| 成人无遮挡网站| 久久久久久久久久久丰满| 国产成人精品福利久久| 亚洲精品日韩av片在线观看| 丰满迷人的少妇在线观看| 久久精品国产亚洲网站| 亚洲性久久影院| 在线观看免费日韩欧美大片 | tube8黄色片| 色网站视频免费| 亚洲伊人久久精品综合| 啦啦啦视频在线资源免费观看| 久久人人爽人人片av| 日韩电影二区| 成年美女黄网站色视频大全免费 | 免费看av在线观看网站| 亚洲av福利一区| 欧美日韩视频高清一区二区三区二| 99热全是精品| 赤兔流量卡办理| 交换朋友夫妻互换小说| 亚州av有码| 午夜91福利影院| 免费黄网站久久成人精品| 欧美激情国产日韩精品一区| 免费高清在线观看视频在线观看| 国产精品欧美亚洲77777| 国产亚洲精品久久久com| 亚洲图色成人| 久久精品人人爽人人爽视色| 3wmmmm亚洲av在线观看| 人人妻人人添人人爽欧美一区卜| 一级黄片播放器| 80岁老熟妇乱子伦牲交| 91久久精品电影网| 久久人人爽人人爽人人片va| 一本久久精品| 91成人精品电影| .国产精品久久| 最新中文字幕久久久久| 91在线精品国自产拍蜜月| 性色av一级| 熟女av电影| 欧美日韩成人在线一区二区| 中文字幕最新亚洲高清| 秋霞在线观看毛片| 亚洲精华国产精华液的使用体验| 日韩av在线免费看完整版不卡| 夫妻午夜视频| 久久精品夜色国产| 久久久久久久久久成人| 国产极品粉嫩免费观看在线 | 国产成人aa在线观看| 欧美97在线视频| 高清黄色对白视频在线免费看| 欧美精品国产亚洲| 亚洲国产精品成人久久小说| 80岁老熟妇乱子伦牲交| 伊人亚洲综合成人网| 久久久精品94久久精品| 99热全是精品| 欧美人与善性xxx| 又大又黄又爽视频免费| 国产成人免费无遮挡视频| 亚洲国产成人一精品久久久| 国产精品.久久久| 一级毛片 在线播放| 2021少妇久久久久久久久久久| 天天影视国产精品| 久久久久久久久久久久大奶| 亚洲内射少妇av| 免费黄色在线免费观看| 中国美白少妇内射xxxbb| 亚洲美女黄色视频免费看| 免费观看性生交大片5| 国产欧美日韩一区二区三区在线 | 成人漫画全彩无遮挡| 亚洲成人一二三区av| 三级国产精品欧美在线观看| 中文字幕av电影在线播放| 国产在线一区二区三区精| 高清黄色对白视频在线免费看| 国产精品女同一区二区软件| 自线自在国产av| 日本黄大片高清| 国产黄片视频在线免费观看| 免费黄频网站在线观看国产| 天天操日日干夜夜撸| 欧美日韩视频高清一区二区三区二| 在线观看人妻少妇| 18在线观看网站| 亚洲欧美日韩卡通动漫| 国产亚洲精品久久久com| 在线播放无遮挡| 爱豆传媒免费全集在线观看| 成人二区视频| 免费观看在线日韩| 狂野欧美激情性xxxx在线观看| 在线亚洲精品国产二区图片欧美 | 观看美女的网站| 老女人水多毛片| 国产精品国产三级专区第一集| 日韩三级伦理在线观看| 黄片播放在线免费| a级毛片黄视频| 国产探花极品一区二区| 又大又黄又爽视频免费| 欧美 亚洲 国产 日韩一| 伊人久久精品亚洲午夜| 大片电影免费在线观看免费| 99热全是精品| 亚洲成色77777| 久久人人爽人人爽人人片va| 天堂俺去俺来也www色官网| 亚洲精品美女久久av网站| 亚洲美女视频黄频| 亚洲综合精品二区| 国产色婷婷99| 亚洲av中文av极速乱| 精品国产国语对白av| 国产在线视频一区二区| 国产精品欧美亚洲77777| 免费看av在线观看网站| 狠狠婷婷综合久久久久久88av| 国产男女内射视频| 国产精品久久久久久精品古装| 99视频精品全部免费 在线| 91久久精品国产一区二区三区| 丁香六月天网| 超色免费av| 最后的刺客免费高清国语| h视频一区二区三区| 亚洲精品,欧美精品| 午夜老司机福利剧场| 亚洲国产精品999| 一区二区三区精品91| 精品少妇内射三级| 青春草视频在线免费观看| 黄片无遮挡物在线观看| 插阴视频在线观看视频| 99国产精品免费福利视频| 97超视频在线观看视频| 韩国高清视频一区二区三区| 人体艺术视频欧美日本| 97在线视频观看| 中文字幕av电影在线播放| 九九爱精品视频在线观看| 国产免费福利视频在线观看| 国产日韩一区二区三区精品不卡 | 日韩一区二区视频免费看| 国产一区二区三区av在线| 下体分泌物呈黄色| 国产高清不卡午夜福利| 18禁在线播放成人免费| 欧美精品高潮呻吟av久久| av免费观看日本| 国产视频内射| 亚洲av国产av综合av卡| 精品人妻熟女av久视频| 91午夜精品亚洲一区二区三区| 精品久久久噜噜| 国产午夜精品久久久久久一区二区三区| 国产精品熟女久久久久浪| 亚洲人成77777在线视频| 99九九线精品视频在线观看视频| 91成人精品电影| 大香蕉久久成人网| 亚洲精品国产av成人精品| 伊人久久精品亚洲午夜| 久久久精品94久久精品| 少妇精品久久久久久久| 亚洲国产精品专区欧美| 少妇熟女欧美另类| 在线看a的网站| 91精品伊人久久大香线蕉| 午夜福利,免费看| 国产男人的电影天堂91| 丝瓜视频免费看黄片| 亚洲一区二区三区欧美精品| 日韩免费高清中文字幕av| 精品人妻熟女毛片av久久网站| 自拍欧美九色日韩亚洲蝌蚪91| 国产一级毛片在线| 亚洲精品aⅴ在线观看| 少妇猛男粗大的猛烈进出视频| 午夜免费男女啪啪视频观看| www.色视频.com| 少妇熟女欧美另类| 全区人妻精品视频| 在线观看免费高清a一片| 日韩在线高清观看一区二区三区| 美女国产视频在线观看| 999精品在线视频| 国产永久视频网站| 18禁在线无遮挡免费观看视频| 日韩在线高清观看一区二区三区| 国产精品三级大全| 夜夜爽夜夜爽视频| 成人国产av品久久久| 国产色婷婷99| 国产国语露脸激情在线看| 韩国av在线不卡| 三级国产精品欧美在线观看| 18禁在线无遮挡免费观看视频| 久久av网站| 男女边吃奶边做爰视频| 中文字幕人妻熟人妻熟丝袜美| 中文欧美无线码| 亚洲欧美中文字幕日韩二区| 久久久久久久久久人人人人人人| 免费大片18禁| 一级毛片aaaaaa免费看小| 亚洲伊人久久精品综合| 妹子高潮喷水视频| 91久久精品电影网| 午夜老司机福利剧场| 91精品三级在线观看| 国产男女内射视频| 丝袜喷水一区| 国产国语露脸激情在线看| 国产亚洲一区二区精品| 伊人亚洲综合成人网| 国产精品久久久久成人av| 日韩熟女老妇一区二区性免费视频| 亚洲欧洲精品一区二区精品久久久 | 午夜福利网站1000一区二区三区| 国产精品久久久久久av不卡| 欧美老熟妇乱子伦牲交| 高清av免费在线| 精品亚洲成a人片在线观看| 国产 一区精品| 国产 精品1| 亚洲欧美精品自产自拍| av卡一久久| 大又大粗又爽又黄少妇毛片口| 超色免费av| 日韩 亚洲 欧美在线| 免费观看av网站的网址| av国产久精品久网站免费入址| 黄色配什么色好看| 国产综合精华液| 日韩成人伦理影院| 国产亚洲午夜精品一区二区久久| 亚洲欧美色中文字幕在线| 亚洲国产av影院在线观看| 中文字幕制服av| 日本午夜av视频| 999精品在线视频| 一区二区三区精品91| 亚洲人与动物交配视频| 我的老师免费观看完整版| 99久久精品国产国产毛片| 99国产综合亚洲精品| 好男人视频免费观看在线| 国产精品.久久久| av国产精品久久久久影院| 国产成人一区二区在线| 国产视频首页在线观看| 秋霞伦理黄片| 免费久久久久久久精品成人欧美视频 | 在线观看三级黄色| 亚洲人成网站在线观看播放| 人人澡人人妻人| 色94色欧美一区二区| 99久国产av精品国产电影| 22中文网久久字幕| 久久久精品区二区三区| av播播在线观看一区| 麻豆精品久久久久久蜜桃| 午夜日本视频在线| 亚洲精品第二区| 18在线观看网站| 久久人人爽人人片av| 日韩av不卡免费在线播放| 人妻少妇偷人精品九色| 欧美精品一区二区大全| 麻豆精品久久久久久蜜桃| 精品人妻偷拍中文字幕| 国产精品一区www在线观看| 人成视频在线观看免费观看| 国产老妇伦熟女老妇高清| 国产免费一级a男人的天堂| 女性被躁到高潮视频| 在线看a的网站| a 毛片基地| 男女啪啪激烈高潮av片| 少妇 在线观看| freevideosex欧美| 99热6这里只有精品| 亚洲精品久久久久久婷婷小说| 青春草国产在线视频| 只有这里有精品99| 国产无遮挡羞羞视频在线观看| 国产 一区精品| 18禁动态无遮挡网站| 少妇人妻 视频| 新久久久久国产一级毛片| 国产精品国产三级国产av玫瑰| 亚洲第一av免费看| 亚洲内射少妇av| 少妇被粗大的猛进出69影院 | 久久国产精品大桥未久av| av国产久精品久网站免费入址| 国产av精品麻豆| 我要看黄色一级片免费的| 亚洲av日韩在线播放| a 毛片基地| 最新的欧美精品一区二区|