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

    睡蓮葉片胎生發(fā)育轉(zhuǎn)錄組分析

    2021-01-13 05:05:15蘇群田敏李春牛李先民盧家仕黃展文李杰梅卜朝陽(yáng)王虹妍
    熱帶作物學(xué)報(bào) 2021年12期
    關(guān)鍵詞:差異表達(dá)基因睡蓮

    蘇群 田敏 李春牛 李先民 盧家仕 黃展文 李杰梅 卜朝陽(yáng) 王虹妍

    摘 ?要:睡蓮葉片胎生現(xiàn)象是繁育途徑的重要補(bǔ)充,對(duì)種群的傳播、擴(kuò)散和生態(tài)環(huán)境的適應(yīng)性有重要作用。通過(guò)轉(zhuǎn)錄組測(cè)序技術(shù)篩選和分析睡蓮葉片胎生現(xiàn)象相關(guān)的代謝路徑和調(diào)控基因,為深入認(rèn)識(shí)睡蓮葉片胎生發(fā)育的分子機(jī)制提供參考。以葉片具有胎生現(xiàn)象的‘小花睡蓮’(X)和葉片無(wú)胎生現(xiàn)象的‘藍(lán)星睡蓮’(L)為材料,利用RNA-Seq技術(shù)對(duì)葉片4個(gè)發(fā)育階段的葉臍部分進(jìn)行生物信息學(xué)分析。分析對(duì)照(L)和樣品(X)葉片不同發(fā)育階段測(cè)序結(jié)果:篩選出的差異表達(dá)基因(DEGs)中,34 909個(gè)基因(48.65%)表達(dá)上調(diào),36 850個(gè)基因(51.35%)表達(dá)下調(diào)。DEGs分析顯示,隨著葉片的發(fā)育,X和L上調(diào)基因和下調(diào)基因數(shù)均呈增加趨勢(shì)。對(duì)L1-vs-X1、L4-vs-X4階段的GO和KEGG功能富集分析表明,DEGs主要富集在質(zhì)膜和膜相關(guān)成分、胞外區(qū)域、細(xì)胞壁等相關(guān)的細(xì)胞組分中,涉及到代謝過(guò)程、生物合成和應(yīng)急響應(yīng)等;Pathway代謝通路表明,DEGs主要參與到植物激素信號(hào)轉(zhuǎn)導(dǎo)、苯丙烷類生物合成、氨基酸類代謝、類黃酮生物合成、甘油磷脂類代謝以及細(xì)胞周期相關(guān)等過(guò)程。對(duì)DEGs進(jìn)一步分析,克隆出了4個(gè)可能參與睡蓮葉片胎生發(fā)育的轉(zhuǎn)錄因子。

    關(guān)鍵詞:睡蓮;葉片胎生;轉(zhuǎn)錄組分析;差異表達(dá)基因;代謝通路

    中圖分類號(hào):Q949.746.1 ? ? ?文獻(xiàn)標(biāo)識(shí)碼:A

    Abstract: Leaf vivipary in water lily (Nymphaea) is an important supplement to breeding, which plays an important role in the propagation and adaptability of species. Transcriptome sequencing technology was used to screen and explore genome-wide analysis of regulatory genes and metabolic pathways involved in leaf vivipary in water lily. This study would lay a foundation for further understanding the molecular mechanism of leaf viviparous development. Viviparous N. micrantha (X) and non-viviparous N. colorata (L) were selected as the experimental materials. By using RNA-Seq technology, four stages of leaf development of the leaf stalk and stem and stem join were tested and analyzed using a series of bioinformatics analysis. The result showed that 34 909 (48.65%) DEGs were up-regulated and 36 850 (51.35%) DEGs were down-regulated. DEGs analysis showed that up-regulated genes and down-regulated genes both increased with the leaves development of N. micrantha amd N. colorata. GO and KEGG enrichment analysis of L1-vs-X1 and L4-vs-X4 stage showed that DEGs were mainly enriched in plasmalemma, cytomembrane and extracellular domain cytoderm associated with ?metabolic process, biosynthesis and stimulus response. Pathway metabolism pathways indicated that DEGs were mainly involved in plant hormone signal transduction, phenylpropanoid biosynthesis, amino acids metabolism, flavonoid biosynthesis, glycerolipid metabolism and cell cycle and other processes. Four transcription factors potentially involved in leaf viviparous development of water lily were cloned based on DEGs analysis.

    Keywords: water lily; leaf vivipary; transcriptome sequencing; DEGs; metabolic pathway

    DOI: 10.3969/j.issn.1000-2561.2021.12.009

    睡蓮為睡蓮科(Nymphaeaceae)睡蓮屬(Nymphaea L.)多年生草本植物,花色豐富,花期長(zhǎng),適應(yīng)性與抗逆性強(qiáng),栽培容易且分布廣泛[1]。睡蓮屬為睡蓮科中最大的屬,約50余種(含變種),可分為subg. Nymphaea、subg. Anecphya、subg. Brachyceras、subg. Hydrocallis和subg. Lotos 5個(gè)亞屬[2-3]。睡蓮因其極高的觀賞價(jià)值和在植物進(jìn)化、分類中的重要地位而越來(lái)越受到愛好者、育種家和植物學(xué)家的追捧[4-7]。部分熱帶睡蓮的葉片除進(jìn)行光合作用外,還具有葉片胎生(vivipary)的特性,在其葉臍部位(葉片與葉柄連接處)長(zhǎng)出新的植株[1]。胎生是植物另一種繁殖途徑,可使其種群在較短時(shí)間內(nèi)快速傳播和擴(kuò)散,更易適應(yīng)多變而復(fù)雜的自然環(huán)境,此外胎生對(duì)生物多樣性保護(hù),生態(tài)平衡的發(fā)展具有重要意義[8]。與睡蓮葉片胎生現(xiàn)象相關(guān)的內(nèi)部代謝路徑和調(diào)控基因的研究則鮮見報(bào)道。

    本研究利用Illumina測(cè)序平臺(tái),以葉片具有胎生現(xiàn)象的‘小花睡蓮 ’(Nymphaea micrantha, X)和葉片無(wú)胎生現(xiàn)象的‘藍(lán)星睡蓮’(Nymphaea col-orata, L)為材料,對(duì)葉片4個(gè)發(fā)育階段的葉臍部分進(jìn)行轉(zhuǎn)錄組測(cè)序,并重點(diǎn)比對(duì)分析了L1-vs-X1和L4-vs-X4階段的差異表達(dá)基因和代謝通路富集,以期篩選睡蓮葉片不同發(fā)育模式下的相關(guān)代謝路徑和參與葉片胎生發(fā)育的調(diào)控基因,為睡蓮葉片胎生發(fā)育的分子遺傳機(jī)制研究提供參考。

    1 ?材料與方法

    1.1 ?材料

    ‘小花睡蓮’(X)和‘藍(lán)星睡蓮’(L)均種植于廣西壯族自治區(qū)農(nóng)業(yè)科學(xué)院花卉研究所睡蓮資源圃內(nèi),采用缸栽形式,水缸尺寸為口徑85 cm,深度50 cm,水深常年維持在25 cm左右。取‘X’(樣品)葉片和‘L’(對(duì)照)葉片4個(gè)不同發(fā)育階段的葉臍部位為材料(圖1),將相同基因型的材料6株混合取樣,并重復(fù)取樣1次。取樣后迅速放入液氮中保存。

    1.2 ?方法

    1.2.1 ?RNA提取、文庫(kù)構(gòu)建和轉(zhuǎn)錄組測(cè)序 ?委托攸歸(上海)生物科技有限公司完成RNA提取、文庫(kù)構(gòu)建和轉(zhuǎn)錄組測(cè)序工作。采用CTAB法提取睡蓮總RNA并采用Agilent 2100進(jìn)行質(zhì)檢。RNA樣品檢測(cè)合格后,富集并打斷成短片段,用六堿基隨機(jī)引物(random hexamers)合成一鏈cDNA,然后加入緩沖液、dNTPs、DNA polymerase I和RNase H合成二鏈cDNA,隨后利用AMPure XP beads純化雙鏈cDNA。純化的雙鏈cDNA先進(jìn)行末端修復(fù)、加A尾并連接測(cè)序接頭,再用AMPure XP beads進(jìn)行片段大小選擇。最后進(jìn)行PCR擴(kuò)增,并用AMPure XP beads純化PCR產(chǎn)物,得到最終的文庫(kù)。采用Illumina NovaSeq 6000平臺(tái)進(jìn)行雙末端測(cè)序。

    使用NGS QC Toolkit[9]軟件對(duì)測(cè)序得到的原始數(shù)據(jù)進(jìn)行質(zhì)控以得到高質(zhì)量的clean reads。使用hisat2[10]將clean reads與‘藍(lán)星睡蓮’的參考基因組進(jìn)行比對(duì)(ftp://download.big.ac.cn/gwh/Plants/ Nymphaea_colorata_Nym_GWHAAYW00000000/ GWHAAYW00000000.genome.fasta.gz),比對(duì)結(jié)

    果以二進(jìn)制binary文件即bam文件進(jìn)行儲(chǔ)存。之后使用Cufflinks[11]對(duì)基因定量獲取FPKM值。采用htseq-count[12]軟件獲取落到各個(gè)樣本中基因的reads數(shù)目。

    1.2.2 ?差異表達(dá)基因(DEGs)檢測(cè)、功能注釋和富集分析 ?使用STEM (Short Time-series Ex-pression Miner)軟件對(duì)不同發(fā)育階段的睡蓮葉片中基因的表達(dá)模式進(jìn)行分析。并采用R中的pheatmap函數(shù)進(jìn)行熱圖聚類分析。使用DESeq (2012) Rpackage的estimateSizeFactors函數(shù)對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化,并使用nbinomTest函數(shù)計(jì)算差異比較的P value和fold change值。以P<0.05且fold change值2倍以上(上調(diào)或下調(diào))作為篩選差異基因的標(biāo)準(zhǔn)。利用Blast2GO[12]軟件和KAAS軟件(https://www.genome.jp/tools/kaas/)分別進(jìn)行GO和KEGG[13]注釋和富集分析。重點(diǎn)選擇比對(duì)了L1-vs-X1、L4-vs-X4階段基因的差異表達(dá)情況,統(tǒng)計(jì)同種材料不同時(shí)期和不同材料同一時(shí)期轉(zhuǎn)錄組的DEGs,并進(jìn)行GO和KEGG富集分析,以判定差異基因主要影響的代謝途徑和信號(hào)通路。

    2 ?結(jié)果與分析

    2.1 ?總RNA質(zhì)量檢測(cè)和組裝結(jié)果分析

    ‘小花睡蓮’葉片和‘藍(lán)星睡蓮’葉片4個(gè)不同發(fā)育階段總計(jì)16份樣本,經(jīng)檢測(cè),樣品總RNA質(zhì)量達(dá)到建庫(kù)要求,其中濃度為70~1230 ng/μL,A260/280為2.1~2.2,28S/18S為1.1~1.9,RIN≥7。

    轉(zhuǎn)錄組測(cè)序共獲得134.97 Gb原始數(shù)據(jù)(表1),各樣本Q30值達(dá)到94%以上,GC≥49%,說(shuō)明測(cè)序結(jié)果準(zhǔn)確度較好,可用于后續(xù)分析。

    2.2 ?轉(zhuǎn)錄組基因表達(dá)數(shù)據(jù)分析

    為了分析‘小花睡蓮’和‘藍(lán)星睡蓮’睡蓮葉片轉(zhuǎn)錄組中基因表達(dá)差異情況,做了全轉(zhuǎn)錄組表達(dá)基因熱圖。從圖2可看出,‘小花睡蓮’和‘藍(lán)星睡蓮’中不同發(fā)育階段的葉片中基因表達(dá)差異較為明顯,L1-vs-X1、L4-vs-X4階段有著較為顯著的差異表達(dá)基因。同時(shí)進(jìn)一步分析‘小花睡蓮’和‘藍(lán)星睡蓮’葉片中不同表達(dá)基因的時(shí)空表達(dá)模式,利用STEM軟件對(duì)‘小花睡蓮’和‘藍(lán)星睡蓮’8個(gè)組織進(jìn)行了基因表達(dá)趨勢(shì)(Profile)分析。結(jié)果顯示共檢測(cè)到60個(gè)表達(dá)趨勢(shì),其中‘小花睡蓮’和‘藍(lán)星睡蓮’各30個(gè)(圖3)。在‘小花睡蓮’4個(gè)發(fā)育階段葉片中,表達(dá)趨勢(shì)8、12、4、18、0、14、25、13、21是顯著富集(P<0.05),分別包含2337、1302、1254、1157、1008、953、913、662、643個(gè)表達(dá)基因;在‘藍(lán)星睡蓮’4個(gè)發(fā)育階段葉片中,表達(dá)趨勢(shì)4、8、25、0、21、12是顯著富集(P<0.05),分別包含了3773、1919、1491、1153、838、675、530個(gè)表達(dá)基因。這些基因的大量表達(dá),提示它們?cè)诖税l(fā)育階段發(fā)揮重要作用。

    2.3 ?差異表達(dá)基因的數(shù)據(jù)分析

    如表2所示,在檢測(cè)到的DEGs中,上調(diào)差異基因34 909個(gè)(48.65%),下調(diào)差異基因36 850個(gè)(51.35%)。隨著葉片的發(fā)育,‘小花睡蓮’和‘藍(lán)星睡蓮’上調(diào)基因和下調(diào)基因數(shù)均呈增加趨勢(shì),各階段比對(duì)分析顯示,L1-vs-X1階段的DEGs為5559個(gè),其中2033個(gè)基因表現(xiàn)為上調(diào),3526個(gè)基因表現(xiàn)下調(diào),下調(diào)基因數(shù)明顯多于上調(diào)基因數(shù);L4-vs-X4階段的DEGs最多,高達(dá)7391個(gè),其中2979個(gè)基因上調(diào),4412個(gè)基因下調(diào),同樣也是下調(diào)基因多于上調(diào)基因。L1-vs-X1和L4-vs-X4的DEGs韋恩圖分析表明,3507個(gè)DEGs在2個(gè)比對(duì)階段皆有表達(dá),3884個(gè)DEGs在L4-vs-X4中差異表達(dá),2052個(gè)DEGs在L1-vs-X1中表達(dá)(圖4)。根據(jù)試驗(yàn)設(shè)計(jì)及數(shù)據(jù)結(jié)果初步推測(cè)3884個(gè)基因中可能存在睡蓮葉片胎生發(fā)育相關(guān)的基因。

    2.4 ?差異表達(dá)基因的GO功能分析

    為了進(jìn)一步了解睡蓮葉片胎生發(fā)育分子機(jī)制,重點(diǎn)分析了L1-vs-X1階段的5559個(gè)DEGs以及L4-vs-X4階段的7391個(gè)DEGs的GO富集分析top30(篩選3種分類中對(duì)應(yīng)基因數(shù)目大于2的GO條目,按照每個(gè)條目對(duì)應(yīng)的-log10 P value由大到小排序的各前10條)。從圖5可知,DEGs在L1-vs-X1及L4-vs-X4階段均主要富集在細(xì)胞組分中,具體到質(zhì)膜和膜相關(guān)成分、核小體、類囊體、胞外區(qū)域、細(xì)胞壁等的細(xì)胞組分中,涉及到代謝過(guò)程、生物合成等;在生物過(guò)程中,較多的DEGs歸類為次生代謝、光合呼吸反應(yīng)、防衛(wèi)反應(yīng)、鐵離子運(yùn)輸、應(yīng)急響應(yīng)及DNA復(fù)制和修復(fù)等;在分子功能分類中,參與結(jié)合、轉(zhuǎn)運(yùn)和催化活性的DEGs最多。

    2.5 ?差異表達(dá)基因的KEGG功能分析

    KEGG是有關(guān)Pathway的主要公共數(shù)據(jù)庫(kù),可以利用KEGG數(shù)據(jù)庫(kù)對(duì)差異Unigene進(jìn)行Pathway富集分析(結(jié)合KEGG注釋結(jié)果)。L1-vs-X1階段以及L4-vs-X4階段的KEGG富集分析top20(以P≤0.01為篩選標(biāo)準(zhǔn),同時(shí)過(guò)濾掉差異基因數(shù)目小于3的條目)如表3~表4所示。Pathway代謝通路表明DEGs主要參與到植物激素信號(hào)轉(zhuǎn)導(dǎo)、苯丙烷類生物合成、氨基酸類代謝、細(xì)胞周期、光合作用、類黃酮生物合成以及磷脂類代謝等相關(guān)過(guò)程。

    L1-vs-X1階段植物激素信號(hào)轉(zhuǎn)導(dǎo)通路富集DEGs最多達(dá)到17個(gè),7個(gè)基因表達(dá)上調(diào),10個(gè)基因下調(diào)表達(dá);其次,光合作用固碳路徑中有8個(gè)基因上調(diào),2個(gè)基因下調(diào);氨基酸代謝路徑富集9個(gè)DEGs,6個(gè)基因上調(diào),3個(gè)基因表達(dá)下調(diào);磷脂類代謝及葉綠色代謝類路徑各有8個(gè)DEGs,此外囊泡轉(zhuǎn)運(yùn)互作路徑、DNA復(fù)制、類胡蘿卜素生物合成等路徑DEGs均達(dá)到了5個(gè)。L4-vs-X4階段植物激素信號(hào)轉(zhuǎn)導(dǎo)通路富集DEGs最多達(dá)18個(gè),其中7個(gè)基因表達(dá)上調(diào),11個(gè)基因表達(dá)下調(diào);苯丙烷生物合成富集DEGs 15個(gè),12個(gè)基因表達(dá)上調(diào),3個(gè)基因表達(dá)下調(diào);重點(diǎn)指出的是參與細(xì)胞周期-減數(shù)分裂相關(guān)的DEGs多達(dá)35個(gè),其中22個(gè)基因上調(diào),13個(gè)基因下調(diào)。此外,甘油磷脂類代謝路徑DEGs有10個(gè),上調(diào)和下調(diào)表達(dá)基因各有5個(gè);富集參與到脂肪酸代謝、類黃酮代謝、類固醇激素代謝及葉綠素代謝通路的DEGs也較多。綜上可知,L1-vs-X1階段及L4-vs-X4階段的植物激素信號(hào)轉(zhuǎn)導(dǎo)通路富集DEGs都是最多,且下調(diào)基因數(shù)多于上調(diào)基因數(shù)。

    2.6 ?激素信號(hào)轉(zhuǎn)導(dǎo)途徑中顯著差異表達(dá)基因的篩選

    為了探尋參與睡蓮葉片胎生的內(nèi)部調(diào)控機(jī)制的基因,通過(guò)重點(diǎn)分析‘小花睡蓮’和‘藍(lán)星睡蓮’轉(zhuǎn)錄組數(shù)據(jù),篩選‘小花睡蓮’顯著上調(diào)轉(zhuǎn)錄因子,并通過(guò)轉(zhuǎn)錄組reads組裝獲得基因全長(zhǎng)。

    以此為基礎(chǔ),設(shè)計(jì)相關(guān)基因擴(kuò)增引物,通過(guò)PCR克隆方法,得到4個(gè)轉(zhuǎn)錄因子基因,分別為ERF1B(NCBI登錄號(hào):MT542693)、ERF105(NCBI登錄號(hào):MT974544)、RAP2-3(NCBI登錄號(hào):MT799791)和WRKY22(NCBI登錄號(hào):MT 799792)。

    3 ?討論

    葉片通常被認(rèn)為是植物進(jìn)行光合作用、呼吸作用和蒸騰作用的主要場(chǎng)所,部分葉片也具有觀賞價(jià)值和營(yíng)養(yǎng)價(jià)值[14]。一些植物葉片還能夠從葉緣或者葉片中央部位長(zhǎng)出新的完整植株,如伽藍(lán)菜屬的大葉落地生根(Kalanchoe daigremontiana)[15]和睡蓮屬(Nymphaea spp.)的一些種(品種)[1, 8]等。葉片、花朵等植物營(yíng)養(yǎng)器官的胎生是植物另一種繁殖途徑,可使其種群在較短時(shí)間內(nèi)快速傳播和擴(kuò)散,更易適應(yīng)復(fù)雜而多變的自然環(huán)境。此外胎生對(duì)生物多樣性保護(hù)及生態(tài)平衡的發(fā)展具有重要意義,這種胎生能力受內(nèi)部基因的調(diào)控,具有可遺傳性,植物通過(guò)這種方式可將進(jìn)化的優(yōu)良性狀一代代遺傳固定下來(lái)[8]。部分熱帶睡蓮也多以葉胎生的形式繁殖,但其形態(tài)發(fā)育和分子水平的內(nèi)部調(diào)控的機(jī)理尚不清楚,有待于進(jìn)一步研究。

    本研究中L1-vs-X1和L4-vs-X4階段的植物激素信號(hào)轉(zhuǎn)導(dǎo)通路富集的DEGs都最多,說(shuō)明在調(diào)控睡蓮葉片胎生苗發(fā)育過(guò)程中,植物激素信號(hào)一直發(fā)揮著重要作用,這與Kane等[16]探討植物激素對(duì)睡蓮胎生苗的研究的結(jié)果一致。在L1-vs-X1階段上調(diào)基因7個(gè),下調(diào)基因10個(gè);L4-vs-X4階段上調(diào)基因7個(gè),下調(diào)基因11個(gè),2個(gè)比對(duì)階段均有較多的上下調(diào)基因,維持著動(dòng)態(tài)平衡。此外通過(guò)胎生和非胎生睡蓮種間轉(zhuǎn)錄組的比較試驗(yàn),發(fā)現(xiàn)生長(zhǎng)素(AUX/IAA)、生長(zhǎng)素應(yīng)答因子(ARF)、植物生長(zhǎng)素酰胺合成酶(GH3)和生長(zhǎng)素上調(diào)小RNA(SAUR)的顯著差異調(diào)節(jié)。Aux/IAA是早期反應(yīng)基因,能夠精確快速地觸發(fā)基因重編程[17]。ARF轉(zhuǎn)錄因子在生長(zhǎng)素感知時(shí)被激活,并啟動(dòng)下游信號(hào)通路,包括SAUR基因。SAURs調(diào)節(jié)許多生長(zhǎng)素介導(dǎo)的反應(yīng),特別是通過(guò)細(xì)胞伸長(zhǎng)的組織生長(zhǎng)[18]。這些基因可能通過(guò)參與細(xì)胞分裂、擴(kuò)大和分化而直接參與睡蓮胎生外部生長(zhǎng)的發(fā)育[17]。生長(zhǎng)素與細(xì)胞因子的結(jié)合也有助于細(xì)胞分化,細(xì)胞分裂素(CRE)的高表達(dá)可能是胎生細(xì)胞分化的原因之一,因?yàn)槠鋵?duì)莖頂端分生組織的活性是已知的[19]。此外,有研究表明其他激素如乙烯等和IAA的互作來(lái)調(diào)控胎生苗根系的發(fā)育[20]。郭書磊等[21]利用RNA-Seq技術(shù)研究玉米葉片形態(tài)建成也得出多種植物激素相互作用通過(guò)參與調(diào)控玉米葉片寬窄的形態(tài)建成研究相同。綜上所述,植物激素間的相互作用和動(dòng)態(tài)平衡對(duì)睡蓮葉片胎生發(fā)育有著重要的影響。

    隨著近年鄉(xiāng)村振興和休閑旅游的興起,睡蓮因其迷人的外表而得到快速的推廣應(yīng)用,但性狀優(yōu)良的睡蓮品種往往不具有胎生能力,種苗繁育困難,導(dǎo)致價(jià)格居高不下,嚴(yán)重制約了其進(jìn)一步的推廣應(yīng)用。雖然目前在胎生睡蓮的形態(tài)解剖、分子水平上取得了一些進(jìn)展,但引起胎生現(xiàn)象的內(nèi)部調(diào)控機(jī)制并沒(méi)有得到系統(tǒng)解答,仍需進(jìn)一步在分子水平上探究形成機(jī)理和挖掘重要功能基因。轉(zhuǎn)錄因子在植物的生長(zhǎng)發(fā)育及其對(duì)外界環(huán)境的反應(yīng)中起著重要調(diào)控作用,其序列和結(jié)構(gòu)的多樣性決定了其功能的多樣性。通過(guò)睡蓮種間胎生和非胎生轉(zhuǎn)錄組的比較研究,挖掘出4個(gè)顯著差異表達(dá)的轉(zhuǎn)錄因子,分別為ERF1B、ERF105、RAP2-3和WRKY22。通過(guò)前人的研究得知AP2/ERF家族有多個(gè)成員能夠調(diào)節(jié)體細(xì)胞胚胎發(fā)生[22],其中被研究得最多的是BABY BOOM (BBM)基因[23-27]和WOUND INDUCEDDEDIFFERENTI ATION 1 (WIND1)或RAP2.4基因[28]。也有研究表明WRKY家族轉(zhuǎn)錄因子通過(guò)參與其他激素介導(dǎo)的信號(hào)轉(zhuǎn)導(dǎo)來(lái)間接調(diào)節(jié)植物體細(xì)胞的發(fā)育[29-30]。AP2/ERF基因家族和WRKY等基因家族成員多而復(fù)雜,成員之間的功能差異很大,需具體到特定物種中進(jìn)行功能驗(yàn)證,如在蘋果中ERF1B基因能夠通過(guò)調(diào)控LOX途徑進(jìn)行香氣的合成[31]。而在棉花中則可能參與抗黃萎病的調(diào)控[32]。

    本研究通過(guò)轉(zhuǎn)錄組角度分析出植物激素間的相互作用和動(dòng)態(tài)平衡對(duì)睡蓮葉片胎生發(fā)育有著重要的影響,并結(jié)合RT-PCR技術(shù)克隆了ERF1B、ERF105、RAP2-3和WRKY22等4個(gè)可能參與葉片胎生過(guò)程調(diào)控的轉(zhuǎn)錄因子基因,以期為睡蓮葉片胎生的內(nèi)部代謝和基因調(diào)控提供一定理論基礎(chǔ),今后還需對(duì)克隆出的4個(gè)轉(zhuǎn)錄因子作進(jìn)一步的研究探討。

    參考文獻(xiàn)

    [1] 黃國(guó)振, 鄧惠勤, 李祖修. 睡蓮[M]. 北京: 中國(guó)林業(yè)出版社, 2009: 1-22.

    [2] Borsch T, L?hne C, Wiersema J H, et al. Towards a com-plete species tree of Nymphaea: shedding further light on subg[J]. Telopea, 2011, 13(1-2): 193-217.

    [3] 李淑娟, 尉 ?倩, 陳 ?塵, 等. 中國(guó)睡蓮屬植物育種研究進(jìn)展[J]. 植物遺傳資源學(xué)報(bào), 2019, 20(4): 829-835.

    [4] Kirie S, Iwasaki H, Noshita K, et al. A theoretical morpho-logical model for quantitative description of the three- di-mensional floral morphology in water lily (Nymphaea)[J]. PLoS One, 2020, 15(10): e0239781.

    [5] Povilus R A, Dacosta J M, Grassa C, et al. Water lily (Nym-phaea thermarum) genome reveals variable genomic signa-tures of ancient vascular cambium losses[J]. Proceedings of the National Academy of Sciences of the United States of America, 2020, 117(15): 8649-8656.

    [6] Zhang L, Chen F, Zhang X, et al. The water lily genome and the early evolution of flowering plants[J]. Nature, 2020, 577(7788): 79-84.

    [7] 蘇 ?群, 楊亞涵, 田 ?敏, 等. 睡蓮種質(zhì)資源遺傳多樣性分析及DNA指紋圖譜構(gòu)建[J]. 熱帶作物學(xué)報(bào), 2020, 41(2): 258-266.

    [8] 艾星梅, 陳龍清, 李宇航, 等. 植物營(yíng)養(yǎng)體胎生研究進(jìn)展[J]. 熱帶亞熱帶植物學(xué)報(bào), 2020, 28(2): 209-216.

    [9] Patel R K, Jain M. NGS QC Toolkit: a toolkit for quality control of next generation sequencing data[J]. PLoS One, 2012, 7(2): e30619.

    [10] Langmead B, Salzberg S L. Fast gapped-read alignment with Bowtie 2[J]. Nature Methods, 2012, 9(4): 357-359.

    [11] Trapnell C, Roberts A, Goff L, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks[J]. Nature Protocols, 2012, 7(3): 562- 578.

    [12] Conesa A, G?tz S, García-Gómez J M, et al. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research[J]. Bioinformatics, 2005, 21(18): 3674-3676.

    [13] Kanehisa M, Araki M, Goto S, et al. KEGG for linking genomes to life and the environment[J]. Nucleic Acids Re-search, 2008, 36(database issue): D480-D484.

    [14] 嚴(yán) ?松, 嚴(yán)長(zhǎng)杰, 顧銘洪. 植物葉發(fā)育的分子機(jī)理[J]. 遺傳, 2008, 30(9): 1127-1135.

    [15] 鐘天秀. 大葉落地生根胎生苗差減cDNA文庫(kù)的構(gòu)建及相關(guān)基因的克隆[D]. 北京: 北京林業(yè)大學(xué), 2015.

    [16] Kane M, Kauth P. Use of plant growth regulators to enhance viviparous waterlily propagation[J]. Water Garden Journal, 2006, 21(4): 7-10.

    [17] Luo J, Zhou J J, Zhang J Z. Aux/IAA gene family in plants: molecular structure, regulation, and function[J]. International Journal of Molecular Sciences, 2018, 19(1): 259.

    [18] Charles C. SAUR15 connects auxin perception to lateral and adventitious root formation[J], Plant Physiology, 2020, 184(2): 558-559.

    [19] Inoue T, Higuchi M, Hashimoto Y, et al. Identification of CRE1 as a cytokinin receptor from Arabidopsis[J]. Nature, 2001, 409(6823): 1060-1063.

    [20] Kulka R G. Hormonal control of root development on epi-phyllous plantlets of Bryophyllum (Kalanchoe) marnieria-num: role of auxin and ethylene[J]. Journal of Experimental Botany, 2008, 59(9): 2361-2370.

    [21] 郭書磊, 魯曉民, 齊建雙, 等. 利用RNA-Seq發(fā)掘玉米葉片形態(tài)建成相關(guān)的調(diào)控基因[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2020, 53(1): 1-17.

    [22] Mantiri F R, Kurdyukov S, Chen S K, et al. The transcription factor MtSERF1 may function as a Nexus between stress and development in somatic embryogenesis in Medicago truncatula[J]. Plant Signal Behav, 2008, 3(7): 498-500.

    [23] Boutilier K, Offringa R, Sharma V K, et al. Ectopic expres-sion of BABY BOOM triggers a conversion from vegetative to embryonic growth[J]. Plant Cell, 2002, 14(8): 1737-1749.

    [24] Horstman A, Li M, Heidmann I, et al. The BABY BOOM transcription factor activates the LEC1-ABI3-FUS3-LEC2 network to induce somatic embryogenesis[J]. Plant Physiology, 2017, 175(2): 848-857.

    [25] Liu D, Chen X, Liu J, et al. The rice ERF transcription factor OsERF922 negatively regulates resistance to Magnaporthe oryzae and salt tolerance[J]. Journal of Experimental Botany, 2012, 63(10): 3899-3911.

    [26] Lowe K, Wu E, Wang N, et al. Morphogenic regulators Baby boom and Wuschel improve monocot transformation[J]. Plant Cell, 2016, 28(9): 1998-2015.

    [27] Zhang Z, Zhao H, Li W, et al. Genome-wide association study of callus induction variation to explore the callus formation mechanism of rice[J]. Journal of Integrative Plant Biology, 2019, 61(11): 1134-1150.

    [28] Ikeuchi M, Sugimoto K, Iwase A. Plant callus: mechanisms of induction and repression[J]. Plant Cell, 2013, 25(9): 3159- 3173.

    [29] Shang Y, Yan L, Liu Z Q, et al. The Mg-chelatase H subunit of Arabidopsis antagonizes a group of WRKY transcription repressors to relieve ABA-responsive genes of inhibition[J]. Plant Cell, 2010, 22(6): 1909-1935.

    [30] Zou X, Seemann J R, Neuman D, et al. A WRKY gene from creosote bush encodes an activator of the abscisic acid sig-naling pathway[J]. The Journal of Biological Chemistry, 2004, 279(53): 55770-55779.

    [31] 岳璇璇. 蘋果乙烯響應(yīng)因子ERF1B調(diào)控LOX途徑香氣合成的機(jī)理研究[D]. 泰安: 山東農(nóng)業(yè)大學(xué), 2020.

    [32] 劉 ?坤. 海島棉ERF族B3和B1亞組轉(zhuǎn)錄因子基因的克隆和特征研究[D]. 北京: 中國(guó)農(nóng)業(yè)科學(xué)院, 2011.

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

    猜你喜歡
    差異表達(dá)基因睡蓮
    你好,睡蓮
    睡蓮盛放
    睡蓮
    睡蓮
    鹽脅迫對(duì)花生硝酸鹽積累及信號(hào)轉(zhuǎn)導(dǎo)的影響
    生物信息學(xué)分析患有乳腺癌的乳腺球樣本中與自我更新相關(guān)的關(guān)鍵基因
    睡蓮
    豆腐睡蓮
    條斑紫菜優(yōu)良品系的基因芯片表達(dá)譜分析
    高溫脅迫下草坪草高羊茅差異表達(dá)基因的分子研究
    狂野欧美白嫩少妇大欣赏| 有码 亚洲区| 噜噜噜噜噜久久久久久91| 一本一本综合久久| 成人鲁丝片一二三区免费| 精品日产1卡2卡| 亚洲精品久久国产高清桃花| 国内久久婷婷六月综合欲色啪| 亚洲人成网站高清观看| 亚洲国产精品合色在线| 欧美一区二区精品小视频在线| 久久精品人妻少妇| 九九久久精品国产亚洲av麻豆| 性色av乱码一区二区三区2| 最近最新免费中文字幕在线| 又爽又黄无遮挡网站| xxx96com| 国产精品国产高清国产av| 午夜影院日韩av| 午夜免费观看网址| 精品免费久久久久久久清纯| 欧美成狂野欧美在线观看| 久久久久精品国产欧美久久久| 亚洲人成网站高清观看| 亚洲精品一卡2卡三卡4卡5卡| 国产成人系列免费观看| 国产精品三级大全| 性色avwww在线观看| 男人的好看免费观看在线视频| 亚洲人与动物交配视频| 国产 一区 欧美 日韩| 在线国产一区二区在线| 12—13女人毛片做爰片一| 亚洲成人久久性| 又紧又爽又黄一区二区| 99热精品在线国产| 级片在线观看| av黄色大香蕉| 精品午夜福利视频在线观看一区| 美女cb高潮喷水在线观看| 中文字幕精品亚洲无线码一区| 校园春色视频在线观看| 香蕉久久夜色| 午夜福利在线观看免费完整高清在 | 欧美一区二区亚洲| 午夜老司机福利剧场| 午夜激情福利司机影院| 夜夜夜夜夜久久久久| 老汉色av国产亚洲站长工具| 白带黄色成豆腐渣| 内射极品少妇av片p| 亚洲va日本ⅴa欧美va伊人久久| 国产在视频线在精品| 很黄的视频免费| 一级毛片高清免费大全| 精品久久久久久久久久免费视频| 给我免费播放毛片高清在线观看| 69人妻影院| 国产私拍福利视频在线观看| 国产一区在线观看成人免费| 好男人电影高清在线观看| 动漫黄色视频在线观看| 久久久久性生活片| 男人的好看免费观看在线视频| 久久天躁狠狠躁夜夜2o2o| 噜噜噜噜噜久久久久久91| 国产亚洲欧美在线一区二区| 桃色一区二区三区在线观看| 久久久久精品国产欧美久久久| 欧美日韩国产亚洲二区| 欧美日韩黄片免| 亚洲av成人av| 母亲3免费完整高清在线观看| 中文字幕精品亚洲无线码一区| 中文字幕人妻熟人妻熟丝袜美 | 亚洲成av人片在线播放无| 亚洲成人免费电影在线观看| 欧美av亚洲av综合av国产av| 舔av片在线| 国产精华一区二区三区| 国产精品久久久久久久久免 | 欧美绝顶高潮抽搐喷水| 亚洲最大成人手机在线| 久久久久国内视频| 天美传媒精品一区二区| 岛国在线免费视频观看| 欧美一区二区国产精品久久精品| 琪琪午夜伦伦电影理论片6080| 免费av毛片视频| 中文字幕人妻丝袜一区二区| 五月玫瑰六月丁香| 精品久久久久久,| 内射极品少妇av片p| 99精品在免费线老司机午夜| bbb黄色大片| 一边摸一边抽搐一进一小说| 日日干狠狠操夜夜爽| 亚洲精品粉嫩美女一区| 97超视频在线观看视频| 亚洲国产欧美人成| 国产精品日韩av在线免费观看| 琪琪午夜伦伦电影理论片6080| 亚洲熟妇中文字幕五十中出| 久久久久久久久久黄片| 国产精品一区二区三区四区久久| 成年女人永久免费观看视频| 国产三级中文精品| 女同久久另类99精品国产91| 不卡一级毛片| 国产视频一区二区在线看| 免费人成在线观看视频色| 欧美成人免费av一区二区三区| 国产精华一区二区三区| 精品久久久久久久毛片微露脸| 变态另类丝袜制服| 国产亚洲欧美在线一区二区| 色综合亚洲欧美另类图片| 国内揄拍国产精品人妻在线| 草草在线视频免费看| 日韩成人在线观看一区二区三区| 国产精品影院久久| 大型黄色视频在线免费观看| 2021天堂中文幕一二区在线观| 欧美日本亚洲视频在线播放| 亚洲成人中文字幕在线播放| 婷婷精品国产亚洲av在线| 亚洲精品亚洲一区二区| 十八禁网站免费在线| 最近最新免费中文字幕在线| 黄色丝袜av网址大全| 国内精品美女久久久久久| 日韩欧美三级三区| 欧美色视频一区免费| 国产精品亚洲av一区麻豆| 岛国在线观看网站| 久久这里只有精品中国| 欧美黑人巨大hd| 99在线视频只有这里精品首页| 真人做人爱边吃奶动态| 久久精品人妻少妇| 午夜福利免费观看在线| 中国美女看黄片| 亚洲av电影不卡..在线观看| 网址你懂的国产日韩在线| 人人妻人人澡欧美一区二区| 午夜福利免费观看在线| 国产麻豆成人av免费视频| 人妻丰满熟妇av一区二区三区| 小说图片视频综合网站| 亚洲,欧美精品.| 男插女下体视频免费在线播放| 国产欧美日韩精品一区二区| 午夜福利高清视频| 国产av不卡久久| 可以在线观看的亚洲视频| 高潮久久久久久久久久久不卡| 亚洲专区国产一区二区| 欧美成人a在线观看| 日本黄色片子视频| 观看免费一级毛片| 9191精品国产免费久久| av片东京热男人的天堂| 两个人视频免费观看高清| 欧美日韩中文字幕国产精品一区二区三区| 波多野结衣高清作品| 岛国在线免费视频观看| 国产成人福利小说| 国产亚洲欧美在线一区二区| 中文字幕人妻丝袜一区二区| 国产亚洲欧美在线一区二区| 深爱激情五月婷婷| 制服丝袜大香蕉在线| 国产精品 国内视频| 欧美日韩乱码在线| 免费观看人在逋| 国产淫片久久久久久久久 | 亚洲成av人片在线播放无| 午夜精品在线福利| 欧美日韩福利视频一区二区| 亚洲五月婷婷丁香| 久久性视频一级片| 午夜激情欧美在线| 亚洲人成伊人成综合网2020| 首页视频小说图片口味搜索| 在线观看日韩欧美| 欧美日韩亚洲国产一区二区在线观看| 首页视频小说图片口味搜索| 日本在线视频免费播放| 国产久久久一区二区三区| 夜夜躁狠狠躁天天躁| 国产色爽女视频免费观看| 国产激情偷乱视频一区二区| 中文资源天堂在线| 日韩成人在线观看一区二区三区| 久久九九热精品免费| 看免费av毛片| 久久久久亚洲av毛片大全| 亚洲精品美女久久久久99蜜臀| 亚洲成a人片在线一区二区| 免费av观看视频| 人人妻人人澡欧美一区二区| 男女之事视频高清在线观看| 男女那种视频在线观看| 欧美+亚洲+日韩+国产| 99精品欧美一区二区三区四区| 欧美色欧美亚洲另类二区| 人妻丰满熟妇av一区二区三区| 久久精品国产综合久久久| 网址你懂的国产日韩在线| 色哟哟哟哟哟哟| 中文字幕av成人在线电影| 国内揄拍国产精品人妻在线| 国产亚洲精品久久久com| 又爽又黄无遮挡网站| 在线观看66精品国产| 18禁国产床啪视频网站| 禁无遮挡网站| 国产激情偷乱视频一区二区| 国产蜜桃级精品一区二区三区| 一边摸一边抽搐一进一小说| 国产探花在线观看一区二区| 久久精品夜夜夜夜夜久久蜜豆| 嫩草影院入口| 男女那种视频在线观看| 亚洲美女视频黄频| 91麻豆av在线| 欧美日本亚洲视频在线播放| 熟女人妻精品中文字幕| 午夜福利在线观看免费完整高清在 | 有码 亚洲区| 内地一区二区视频在线| 久久欧美精品欧美久久欧美| 一级作爱视频免费观看| 国产乱人伦免费视频| 国产三级黄色录像| 日本三级黄在线观看| 国产三级中文精品| 别揉我奶头~嗯~啊~动态视频| 久久久久国产精品人妻aⅴ院| 国产成人福利小说| 最近最新中文字幕大全电影3| 免费看a级黄色片| 欧美乱妇无乱码| 日本熟妇午夜| 欧美一区二区亚洲| a级一级毛片免费在线观看| 亚洲精品亚洲一区二区| 露出奶头的视频| 亚洲精品456在线播放app | 最新美女视频免费是黄的| 久久精品91蜜桃| 88av欧美| 婷婷丁香在线五月| 成年女人看的毛片在线观看| 精品久久久久久久毛片微露脸| 日韩欧美国产在线观看| 男女午夜视频在线观看| 亚洲精品久久国产高清桃花| 国产黄片美女视频| 男女床上黄色一级片免费看| aaaaa片日本免费| 两个人看的免费小视频| 亚洲专区中文字幕在线| 亚洲黑人精品在线| 国产毛片a区久久久久| 日韩欧美一区二区三区在线观看| 午夜福利视频1000在线观看| 欧美激情久久久久久爽电影| 99热只有精品国产| 免费在线观看成人毛片| 午夜两性在线视频| 免费电影在线观看免费观看| 亚洲一区高清亚洲精品| 国产精品亚洲一级av第二区| 男女床上黄色一级片免费看| 美女 人体艺术 gogo| 欧美精品啪啪一区二区三区| 一进一出好大好爽视频| 一卡2卡三卡四卡精品乱码亚洲| 成人特级黄色片久久久久久久| 九色成人免费人妻av| 欧美中文综合在线视频| 嫩草影视91久久| 人妻丰满熟妇av一区二区三区| 精品国产超薄肉色丝袜足j| 90打野战视频偷拍视频| 99久久成人亚洲精品观看| 中文亚洲av片在线观看爽| 久久精品影院6| 亚洲欧美激情综合另类| 在线a可以看的网站| 中文字幕人妻丝袜一区二区| 国产成人欧美在线观看| 亚洲一区二区三区色噜噜| 亚洲精品国产精品久久久不卡| 国产97色在线日韩免费| 亚洲va日本ⅴa欧美va伊人久久| 国产乱人视频| 99久久精品热视频| 制服丝袜大香蕉在线| 国产69精品久久久久777片| 天堂网av新在线| 搞女人的毛片| 最新美女视频免费是黄的| 色综合婷婷激情| 日本成人三级电影网站| 99国产极品粉嫩在线观看| 午夜福利高清视频| 亚洲熟妇中文字幕五十中出| 国产成人aa在线观看| 91九色精品人成在线观看| 长腿黑丝高跟| 夜夜夜夜夜久久久久| av天堂中文字幕网| 日本精品一区二区三区蜜桃| 成人18禁在线播放| 久久久久久人人人人人| 最近最新中文字幕大全免费视频| 黑人欧美特级aaaaaa片| 丁香欧美五月| 精品熟女少妇八av免费久了| 欧美大码av| 国产精品99久久99久久久不卡| 亚洲熟妇中文字幕五十中出| 国产成人欧美在线观看| 九九在线视频观看精品| 无遮挡黄片免费观看| 精品熟女少妇八av免费久了| 99视频精品全部免费 在线| 亚洲精品国产精品久久久不卡| 国产色爽女视频免费观看| 无遮挡黄片免费观看| 午夜福利视频1000在线观看| 亚洲一区二区三区不卡视频| 丰满乱子伦码专区| 中文字幕精品亚洲无线码一区| 变态另类成人亚洲欧美熟女| 国产精品免费一区二区三区在线| 欧美3d第一页| 国产精品综合久久久久久久免费| 久久性视频一级片| www日本在线高清视频| 日本精品一区二区三区蜜桃| 美女高潮的动态| 国产av在哪里看| 亚洲精品影视一区二区三区av| 看片在线看免费视频| 午夜福利高清视频| 可以在线观看的亚洲视频| 亚洲色图av天堂| 亚洲熟妇熟女久久| 欧美大码av| 久久久久亚洲av毛片大全| 欧美成人一区二区免费高清观看| netflix在线观看网站| 啦啦啦免费观看视频1| 一区二区三区免费毛片| 欧美日韩国产亚洲二区| 岛国视频午夜一区免费看| a在线观看视频网站| 欧美精品啪啪一区二区三区| 亚洲性夜色夜夜综合| 最近最新中文字幕大全电影3| 亚洲成人久久性| 免费观看的影片在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品色激情综合| 最新美女视频免费是黄的| 色在线成人网| 国产美女午夜福利| 日韩人妻高清精品专区| 禁无遮挡网站| 国产黄片美女视频| 亚洲精品久久国产高清桃花| 婷婷亚洲欧美| 国产精品女同一区二区软件 | 国产精品久久久久久人妻精品电影| 成人欧美大片| av中文乱码字幕在线| 少妇人妻精品综合一区二区 | 无遮挡黄片免费观看| 亚洲性夜色夜夜综合| 在线免费观看不下载黄p国产 | 亚洲不卡免费看| 国产爱豆传媒在线观看| 桃红色精品国产亚洲av| 熟妇人妻久久中文字幕3abv| 亚洲国产欧美人成| 日韩欧美国产在线观看| 国产亚洲精品久久久久久毛片| 俄罗斯特黄特色一大片| 99久久综合精品五月天人人| 亚洲专区中文字幕在线| 亚洲aⅴ乱码一区二区在线播放| 久久久久精品国产欧美久久久| 日韩成人在线观看一区二区三区| 国产欧美日韩一区二区精品| 欧美区成人在线视频| 搡老妇女老女人老熟妇| 天堂网av新在线| 免费高清视频大片| 怎么达到女性高潮| 国产伦精品一区二区三区视频9 | 国产精品一区二区三区四区免费观看 | 1000部很黄的大片| 麻豆成人午夜福利视频| 在线天堂最新版资源| 两个人视频免费观看高清| 国内精品久久久久精免费| 天堂√8在线中文| www.熟女人妻精品国产| 啦啦啦免费观看视频1| 老鸭窝网址在线观看| svipshipincom国产片| 亚洲真实伦在线观看| 毛片女人毛片| 欧美最新免费一区二区三区 | 婷婷六月久久综合丁香| 女人被狂操c到高潮| av天堂中文字幕网| 日韩高清综合在线| 婷婷精品国产亚洲av| 欧美日韩中文字幕国产精品一区二区三区| 国产亚洲精品av在线| 日本黄色片子视频| 老汉色∧v一级毛片| 国模一区二区三区四区视频| 亚洲成av人片在线播放无| a级毛片a级免费在线| 久久久精品大字幕| 热99在线观看视频| 中文字幕人妻熟人妻熟丝袜美 | 国产日本99.免费观看| 国产精品一区二区三区四区免费观看 | 19禁男女啪啪无遮挡网站| 精品人妻一区二区三区麻豆 | 国产精品野战在线观看| 亚洲,欧美精品.| 熟女少妇亚洲综合色aaa.| 日韩 欧美 亚洲 中文字幕| 很黄的视频免费| 91九色精品人成在线观看| 久久6这里有精品| 日韩欧美在线乱码| 久久精品人妻少妇| 欧美成人a在线观看| 日本黄色视频三级网站网址| 美女 人体艺术 gogo| 97人妻精品一区二区三区麻豆| 校园春色视频在线观看| 男女床上黄色一级片免费看| 9191精品国产免费久久| 免费人成视频x8x8入口观看| 亚洲一区二区三区色噜噜| 在线免费观看不下载黄p国产 | 夜夜看夜夜爽夜夜摸| 熟女人妻精品中文字幕| 18禁黄网站禁片免费观看直播| 国产av一区在线观看免费| 欧美黑人欧美精品刺激| 18+在线观看网站| 国产高清有码在线观看视频| 亚洲电影在线观看av| 少妇人妻精品综合一区二区 | 亚洲国产欧洲综合997久久,| 国产熟女xx| 嫩草影院入口| 欧美一区二区精品小视频在线| 亚洲一区高清亚洲精品| 高清日韩中文字幕在线| 国产伦精品一区二区三区四那| 午夜亚洲福利在线播放| 日韩中文字幕欧美一区二区| x7x7x7水蜜桃| 草草在线视频免费看| 人妻丰满熟妇av一区二区三区| 久久久久久大精品| 老司机在亚洲福利影院| 国产欧美日韩一区二区精品| 少妇人妻一区二区三区视频| 级片在线观看| 中文字幕熟女人妻在线| 国产精品98久久久久久宅男小说| 99在线视频只有这里精品首页| 欧美日本亚洲视频在线播放| 精品人妻一区二区三区麻豆 | 国产成人a区在线观看| 久久精品综合一区二区三区| 国产精品一区二区三区四区免费观看 | 成人国产综合亚洲| 国产亚洲精品久久久com| 精品人妻1区二区| 亚洲不卡免费看| 99在线人妻在线中文字幕| 特级一级黄色大片| 免费观看的影片在线观看| 99国产极品粉嫩在线观看| 日韩欧美在线二视频| 哪里可以看免费的av片| 又黄又粗又硬又大视频| 毛片女人毛片| 国产真实伦视频高清在线观看 | 老司机福利观看| 日本精品一区二区三区蜜桃| 熟妇人妻久久中文字幕3abv| 免费观看的影片在线观看| 亚洲精品一区av在线观看| 在线播放国产精品三级| 90打野战视频偷拍视频| av福利片在线观看| 国产精品美女特级片免费视频播放器| 欧美bdsm另类| 免费观看精品视频网站| 最近最新中文字幕大全免费视频| 国产免费一级a男人的天堂| 老司机午夜福利在线观看视频| 他把我摸到了高潮在线观看| 成人精品一区二区免费| 亚洲精华国产精华精| 91字幕亚洲| 国产成人啪精品午夜网站| 天堂√8在线中文| 少妇熟女aⅴ在线视频| 国产精品嫩草影院av在线观看 | 成人国产综合亚洲| 成人性生交大片免费视频hd| 欧美在线黄色| 国产av麻豆久久久久久久| 99国产精品一区二区蜜桃av| 精品一区二区三区人妻视频| 亚洲欧美精品综合久久99| 亚洲国产欧洲综合997久久,| 欧美精品啪啪一区二区三区| 别揉我奶头~嗯~啊~动态视频| 日韩 欧美 亚洲 中文字幕| 精品国内亚洲2022精品成人| 69人妻影院| 日韩欧美 国产精品| 99热6这里只有精品| 村上凉子中文字幕在线| 欧美日韩亚洲国产一区二区在线观看| 亚洲内射少妇av| 老司机深夜福利视频在线观看| 国产精品日韩av在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美一区二区三区在线观看| 欧美国产日韩亚洲一区| 国产精品免费一区二区三区在线| 人妻丰满熟妇av一区二区三区| 黄片小视频在线播放| 国产色婷婷99| 淫秽高清视频在线观看| 国产伦一二天堂av在线观看| 极品教师在线免费播放| 美女被艹到高潮喷水动态| 有码 亚洲区| 亚洲成av人片在线播放无| 99国产精品一区二区蜜桃av| 91字幕亚洲| 99久久成人亚洲精品观看| 99国产综合亚洲精品| 精品一区二区三区视频在线观看免费| 欧美一级毛片孕妇| 成年女人永久免费观看视频| 一级黄色大片毛片| 首页视频小说图片口味搜索| 天堂√8在线中文| 国产黄片美女视频| 99久久精品国产亚洲精品| 国产精品一及| 亚洲va日本ⅴa欧美va伊人久久| 麻豆成人午夜福利视频| 草草在线视频免费看| 国产亚洲精品av在线| 国产亚洲av嫩草精品影院| 国产一区在线观看成人免费| 国产激情欧美一区二区| 一个人观看的视频www高清免费观看| 亚洲av成人不卡在线观看播放网| 最新美女视频免费是黄的| 亚洲欧美一区二区三区黑人| 99久久久亚洲精品蜜臀av| 欧美国产日韩亚洲一区| 香蕉丝袜av| 国产高清视频在线播放一区| 久久久久精品国产欧美久久久| 免费大片18禁| 91久久精品电影网| 亚洲不卡免费看| 偷拍熟女少妇极品色| 精品一区二区三区人妻视频| 亚洲精品乱码久久久v下载方式 | 在线观看66精品国产| 国产成人欧美在线观看| 国产精品三级大全| 熟女电影av网| 免费看日本二区| 国产一区二区在线av高清观看| 一级黄色大片毛片| 老司机午夜福利在线观看视频| 久久精品综合一区二区三区| 日韩中文字幕欧美一区二区| 成熟少妇高潮喷水视频| 中文字幕熟女人妻在线| 国内精品美女久久久久久| 日本三级黄在线观看| 我要搜黄色片| 精品人妻偷拍中文字幕| 制服人妻中文乱码| 久久久久久久精品吃奶| 亚洲精品一区av在线观看| 午夜免费激情av| 啦啦啦免费观看视频1| 欧美高清成人免费视频www| 精品欧美国产一区二区三|