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

    山藥塊莖發(fā)育階段基因共表達(dá)網(wǎng)絡(luò)構(gòu)建及階段特異性分析

    2022-03-16 01:18:35殷劍美張鉛蔣璐張培通
    關(guān)鍵詞:山藥

    殷劍美 張鉛 蔣璐 張培通

    摘要: 山藥是一種藥食同源作物,塊莖是其主要產(chǎn)品器官,其發(fā)育機(jī)制對(duì)于山藥產(chǎn)量和品質(zhì)提升相關(guān)育種有著重要的指導(dǎo)作用。以不同生長階段(塊莖膨大第1 d、第11 d、第21 d、第31 d、第41 d、第51 d,T1~T6)的山藥塊莖為材料進(jìn)行轉(zhuǎn)錄組測序,經(jīng)從頭(de novo)組裝后獲得42 042個(gè)unigenes,利用七大數(shù)據(jù)庫進(jìn)行注釋,鑒定后發(fā)現(xiàn)山藥塊莖形成過程中涉及大量基因,包括激素合成、信號(hào)轉(zhuǎn)導(dǎo)等相關(guān)基因。階段特異性分析結(jié)果表明,塊莖形成初始階段為塊莖發(fā)育的關(guān)鍵時(shí)期;加權(quán)基因共表達(dá)網(wǎng)絡(luò)分析結(jié)果表明,在山藥塊莖形成初期及后期,激素相關(guān)基因的表達(dá)特別活躍,脫落酸(ABA)、茉莉酸(JA)是山藥塊莖生長發(fā)育的重要激素。研究結(jié)果為進(jìn)一步研究山藥塊莖膨大分子機(jī)制奠定了基礎(chǔ)。

    關(guān)鍵詞: 山藥;轉(zhuǎn)錄組測序;階段特異性基因;基因共表達(dá)網(wǎng)絡(luò)

    中圖分類號(hào): S632.101 文獻(xiàn)標(biāo)識(shí)碼: A 文章編號(hào): 1000-4440(2022)01-0030-09

    Abstract: Chinese yam is a crop with hamology of medicine and food, and tuber is its main product organ. The development mechanism of tubers plays an important guiding role in the breeding related to the improvement of yield and quality of yam. Based on transcriptome sequencing of tubers at different growth stages (T1-T6), after de novo assembly, 42 042 unigenes were obtained and annotated using seven databases. After identification, it was found that a large number of genes were involved in the process of tuber formation, including hormone synthesis-related genes, signal transduction-related genes, and so on. The results of stage specificity analysis showed that the initial stage of tuber formation was a critical period. The results based on the weighted gene co-expression network analysis (WGCNA) showed that the expression of hormone-related genes was very active in the early and late stages of tuber growth. Abscisic acid (ABA) and jasmonic acid (JA) were important hormones for the growth and development of yam tubers. The results lay a foundation for further study on the molecular mechanism of tuber enlargement of yam.

    Key words: yam;transcriptome sequencing;stage-specific gene;weighted gene co-expression network analysis

    參薯山藥(Dioscorea alata)是薯蕷科薯蕷屬藥食同源的作物,其地下塊莖是參薯山藥的主要產(chǎn)品器官[1],具有高鉀、低纖維、高蛋白質(zhì)、營養(yǎng)豐富等特點(diǎn)[2]。參薯山藥塊莖在發(fā)育過程中邊膨大邊積累養(yǎng)分,關(guān)于其塊莖形成的分子機(jī)制研究較少。隨著高通量測序技術(shù)的發(fā)展,人們可以從分子水平研究植物的動(dòng)態(tài)發(fā)育[3-4],包括激素和次生代謝物形成的分子機(jī)制等[5]。目前,關(guān)于地下塊莖發(fā)育的研究主要集中在馬鈴薯、甘薯等作物上,且塊莖發(fā)育與激素密切相關(guān)。前人研究發(fā)現(xiàn),ABR1在芥菜塊莖中的表達(dá)量顯著增加,脫落酸(ABA)介導(dǎo)的信號(hào)參與了塊莖膨大[6];在馬鈴薯塊莖膨大早期,調(diào)控赤霉素(GA)降解的基因表達(dá)量上調(diào),GA含量迅速下降,從而促進(jìn)塊莖膨大[7]。Sarker等[8]認(rèn)為,茉莉酸(JA)和茉莉酸甲酯(MeJA)在馬鈴薯塊莖形成中起著重要作用。筆者之前對(duì)參薯山藥塊莖不同發(fā)育時(shí)期激素含量的研究發(fā)現(xiàn),GA在塊莖形成初期的含量較高,可以抑制塊莖形成;ABA和JA對(duì)塊莖膨大具有促進(jìn)作用[9]。

    本研究以不同生長時(shí)期參薯山藥塊莖為材料進(jìn)行轉(zhuǎn)錄組測序,采用階段特異性(Stage specific,SS)方法分析塊莖發(fā)育過程中相關(guān)激素基因的表達(dá),構(gòu)建其基因共表達(dá)網(wǎng)絡(luò)(Weighted gene co-expression network analysis, WGCNA),以期為揭示參薯山藥塊莖形成的分子調(diào)控機(jī)制提供理論依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)材料

    試驗(yàn)材料為米易山藥(Dioscorea alata L.),引自四川省米易縣,屬于國家地理標(biāo)志農(nóng)產(chǎn)品,于2018年4月26日種植于江蘇省農(nóng)業(yè)科學(xué)院六合動(dòng)物科學(xué)基地,當(dāng)年10月24日收獲。塊莖于9月4日(播種后132 d)開始膨大,分別于塊莖膨大的第11 d(9月14日)、21 d(9月24日)、31 d(10月4日)、41 d(10月14日)、51 d(10月24日)取地下塊莖,每次取3株混合,用液氮速凍后保存于-70 ℃冰柜中備用。

    1.2 RNA和cDNA文庫的制備

    采用TRIzol法(Invitrogen,USA)進(jìn)行總RNA的提取,用1%瓊脂糖凝膠電泳檢測RNA是否降解及污染,用Nano光度計(jì)(IMPLEN,USA)檢測RNA的純度(OD260/280),用QubitRNA試劑盒(Life Technologies,USA)測定RNA的濃度,用RNA Nano 6000試劑盒(Agilent Technologies,USA)評(píng)估RNA的完整性。

    取1.5 μg RNA,用poly-T寡聚物磁珠純化mRNA。利用隨機(jī)六聚體引物和M-MuLV逆轉(zhuǎn)錄酶(RNase H)合成第1鏈cDNA,用DNA聚合酶Ⅰ和RNase H合成第2鏈cDNA。用TruSeq PE Cluster Kit v3 cBot HS(Illumia)在cBot聚類生成系統(tǒng)上對(duì)索引編碼樣品進(jìn)行聚類。在Illumina Hiseq平臺(tái)上對(duì)文庫進(jìn)行測序,獲得大小為125 bp/150 bp的雙端reads。

    1.3 測序數(shù)據(jù)質(zhì)量的評(píng)估

    對(duì)測序所得原始數(shù)據(jù)進(jìn)行過濾,去除接頭(未知堿基比例>10%,Qphred值≤20的堿基數(shù)>50%)。

    1.4 拼接及功能注釋

    采用Grabherr[10]對(duì)clean reads進(jìn)行拼接,獲得unigene,隨后進(jìn)行七大數(shù)據(jù)庫(包括NR、NT、PFAM、KOG、Swiss-Prot、KEGG、GO)功能的注釋,同時(shí)基于斯皮爾曼相關(guān)系數(shù)(Spearman correlation coefficient,SCC)進(jìn)行主成分分析(Principal component analysis,PCA)。

    1.5 基因表達(dá)量的計(jì)算和差異基因的篩選

    采用基因表達(dá)水平估算方法[11]對(duì)readcount(read數(shù))進(jìn)行每百萬堿基對(duì)測序的每千堿基轉(zhuǎn)錄序列片段(fragments)數(shù)目(FPKM)估算,計(jì)算公式如下:

    FPKM=total exon fragmentsmapped reads × exon length

    式中,total exon fragments表示總外顯子片段數(shù),mapped reads表示定位read數(shù)(百萬);exon length表示外顯子長度(kB)。

    以P值<0.05且 |log2FoldChange|>1為閾值篩選差異表達(dá)基因(Differentially expressed genes,DEGs)。

    SS(i,j)=1-MaxEikEij,1≤k≤6,k≠j

    1.6 階段特異基因的篩選及共表達(dá)網(wǎng)絡(luò)的構(gòu)建

    通過階段特異性(Stage specific,SS)算法[12],篩選得分>0.5的unigene作為階段特異基因。對(duì)基因變異系數(shù)>0.5的基因進(jìn)行加權(quán)基因共表達(dá)網(wǎng)絡(luò)分析(WGCNA)[13]。

    1.7 實(shí)時(shí)熒光定量PCR驗(yàn)證

    采用PrimeScript反轉(zhuǎn)錄試劑盒(TaKaRa,大連)合成cDNA,根據(jù)轉(zhuǎn)錄組表達(dá)量的差異,隨機(jī)挑選8個(gè)基因,使用Primer 5.0(Premier BIosoft International,USA)設(shè)計(jì)引物,選擇DaActin作為內(nèi)參基因進(jìn)行RT-qPCR。反應(yīng)程序:95? ℃ 30 min;95 ℃ 5 s,58 ℃ 30 s,72 ℃ 30 s,45個(gè)循環(huán)。通過定量反應(yīng)的熔解曲線,驗(yàn)證引物的特異性是否滿足試驗(yàn)要求。RT-qRCR反應(yīng)在Light Cycle 480熒光定量PCR檢測系統(tǒng)(Roche, 瑞士)中進(jìn)行,依據(jù)2-△△Ct算法計(jì)算基因的相對(duì)表達(dá)量[14]。

    2 結(jié)果與分析

    2.1 測序數(shù)據(jù)分析與基因注釋

    將測序的原始數(shù)據(jù)過濾后,平均每個(gè)文庫含有50 266 758個(gè)clean reads,共獲得186 793個(gè)transcripts和42 042個(gè)unigenes。不同F(xiàn)PKM區(qū)間的unigene數(shù)量如圖1所示,可以看出,在不同時(shí)期,大部分基因的FPKM值為 0~0.10、0.31~3.57。

    T1:塊莖膨大第1 d;T1:塊莖膨大的第11 d;T3:塊莖膨大的第21 d;T4:塊莖膨大的第31 d;T5:塊莖膨大的第41 d;T6:塊莖膨大的第51 d。

    分別有46.48%、24.15%、14.62%、56.74%、35.30%、46.48%和44.01%的基因在GO、KEGG(KO)、KOG、NR、NT、PFAM和SwissProt數(shù)據(jù)庫中得到注釋。其中在NR數(shù)據(jù)庫中的注釋結(jié)果表明,與參薯山藥比對(duì)率較高的植物是棕櫚(23.6%)、海棗(18.7%)、菠蘿(7.1%)、小果野蕉(6.8%)和蘆筍(5.6%),它們與參薯山藥均屬于單子葉植物,具有較相似的基因信息。

    由圖2可以看出,在GO數(shù)據(jù)庫中注釋到55個(gè)富集條目,其中包括25個(gè)“生物學(xué)過程”、20個(gè)“細(xì)胞組成”和10個(gè)“分子功能”?!吧飳W(xué)過程”基因主要富集在“細(xì)胞過程”(GO:0009987)、“代謝過程”(GO:0008152)和“單一組織過程”(GO:0044699);“細(xì)胞組成”基因主要富集在“細(xì)胞”(GO:0005623)、“部分細(xì)胞”(GO:0044464)和“細(xì)胞器官”(GO:0043226);“分子功能”基因主要富集在“結(jié)合”(GO:0005488)和“催化活性”(GO:0003824)。

    共有10 154個(gè)unigenes注釋到130個(gè)KEGG通路中,可分為新陳代謝、遺傳信息處理、環(huán)境信息處理、細(xì)胞過程和組織系統(tǒng)5類,其中大部分基因注釋到新陳代謝的代謝通路中(圖3)。

    主成分分析結(jié)果顯示,6個(gè)樣本可明顯分為T1時(shí)期、T2~T6時(shí)期,T1時(shí)期和其他各時(shí)期的相關(guān)性較低,T2~T6時(shí)期較為聚集;塊莖在生長早期(T1時(shí)期)和后期(T2~T6時(shí)期)存在較大差異(圖4)。

    2.2 塊莖發(fā)育過程中不同表達(dá)基因和階段特異基因分析

    分析差異表達(dá)基因發(fā)現(xiàn),兩兩樣品之間的DEGs數(shù)量為442~6 774個(gè),其中與T5時(shí)期相比,T1時(shí)期的DEGs數(shù)量最多,有3 790個(gè)上調(diào)基因,2 984個(gè)下調(diào)基因(表1)。共獲得8 222個(gè)階段特異表達(dá)基因,數(shù)量范圍為599~3 621個(gè)(圖5)。GO分析發(fā)現(xiàn),T1時(shí)期的大多數(shù)基因與蛋白質(zhì)磷酸化、去磷酸化、細(xì)胞蛋白修飾過程、蛋白質(zhì)修飾過程、生長素應(yīng)答、脂質(zhì)轉(zhuǎn)運(yùn)、激素應(yīng)答等有關(guān)。

    柱旁數(shù)據(jù)表示注釋到相應(yīng)基因的數(shù)量。

    2.3 T1時(shí)期差異表達(dá)基因及特異基因分析

    主成分分析和不同發(fā)育階段特異基因分析結(jié)果顯示,T1時(shí)期和T2~T6時(shí)期之間存在明顯差異,在其有表達(dá)差異的基因中,共有211個(gè)基因編碼轉(zhuǎn)錄因子,其中36個(gè)基因編碼AP2/ERF,另外編碼WRKY、C2H2、MYB、NAC、bHLH、bZIP的基因均表現(xiàn)為高表達(dá)。

    通過比較T1時(shí)期外的其他時(shí)期與T1時(shí)期發(fā)現(xiàn),差異表達(dá)的基因達(dá)9 701個(gè),其中4 546個(gè)基因的表達(dá)量上調(diào),5 181個(gè)基因的表達(dá)量下調(diào)。GO分析結(jié)果顯示,上調(diào)表達(dá)的基因與細(xì)胞壁構(gòu)建、細(xì)胞骨架構(gòu)建、DNA復(fù)制、焦磷酸酶活性、染色體、DNA代謝過程調(diào)控有關(guān);下調(diào)表達(dá)的基因與核糖核酸結(jié)合、蛋白激酶活性、蛋白激酶磷酸化活性、大分子代謝、初級(jí)代謝及氮化合物過程的調(diào)控、細(xì)胞大分子和大分子調(diào)節(jié)及細(xì)胞生物合成過程有關(guān)。KEGG分析結(jié)果表明,上調(diào)表達(dá)的基因主要與次生代謝物生物合成、氨基酸合成、代謝和降解、檸檬酸循環(huán)有關(guān);下調(diào)表達(dá)的基因主要涉及植物激素信號(hào)轉(zhuǎn)導(dǎo)、MAPK信號(hào)轉(zhuǎn)導(dǎo)、苯丙氨酸生物合成、酪氨酸生物合成、脂肪酸降解、色氨酸生物合成、玉米素生物合成、生物合成酯類代謝。

    2.4 激素相關(guān)基因轉(zhuǎn)錄譜分析

    對(duì)所有時(shí)期的特異基因進(jìn)行分析發(fā)現(xiàn),與GA相關(guān)的3個(gè)關(guān)鍵基因在整個(gè)生長期的表達(dá)量均較低,0

    2.5 基因共表達(dá)網(wǎng)絡(luò)的構(gòu)建

    對(duì)20 988個(gè)基因進(jìn)行共表達(dá)網(wǎng)絡(luò)構(gòu)建,可將其分為圖6中的25個(gè)不同模塊。進(jìn)一步分析發(fā)現(xiàn),模塊1(綠黃模塊,181個(gè)基因)與T1時(shí)期高度相關(guān),其中有18個(gè)編碼轉(zhuǎn)錄因子,包括轉(zhuǎn)錄因子BHLH148、WRKY與乙烯應(yīng)答相關(guān)轉(zhuǎn)錄因子ERF036 like等。構(gòu)建T1時(shí)期特征模塊基因間聯(lián)通性的共表達(dá)網(wǎng)絡(luò),詳見圖7。

    T1-1、T1-2、T1-3:T1時(shí)期取樣塊莖的3個(gè)重復(fù),其他處理以此類推。共分為25個(gè)不同模塊,其中模塊1(綠黃模塊,181個(gè)基因)與T1時(shí)期高度相關(guān);圖例數(shù)據(jù)表示模塊與樣本間的相關(guān)性。

    正方形節(jié)點(diǎn):核心基因;菱形節(jié)點(diǎn):轉(zhuǎn)錄因子及調(diào)控激酶基因;正六邊形節(jié)點(diǎn):合成和代謝相關(guān)基因;圓形節(jié)點(diǎn):其他基因。

    2.6 RT-qPCR驗(yàn)證

    以Actin為內(nèi)參基因,隨機(jī)選取8個(gè)基因進(jìn)行RT-qPCR,測定其表達(dá)量。引物序列見表3,基因表達(dá)量見圖8?;驘晒舛勘磉_(dá)趨勢與轉(zhuǎn)錄組測序結(jié)果一致,說明轉(zhuǎn)錄組測序結(jié)果可靠。

    左縱坐標(biāo):基于熒光定量法的基因相對(duì)表達(dá)水平;右縱坐標(biāo):基于轉(zhuǎn)錄組測序的FPKM值。A:NCED基因(cluster-5213.7062);B:CYP707A基因(cluster-5213.14981);C:GID1C基因(cluster-5213.13061);D:OPR基因(cluster-5213.15137);E:AOS基因(cluster-5213.15182);F:LOG基因(cluster-5213.2615);G:GBSS基因(cluster-5213.17103);H:SSS基因(cluster-5213.22318)。T1、T2、T3、T4、T5、T6見圖1注。不同發(fā)育時(shí)期標(biāo)有不同小寫字母的表示差異顯著(P<0.05)。

    3 討論

    3.1 山藥塊莖發(fā)育關(guān)鍵時(shí)期分析

    DEGs主要發(fā)生在塊莖形成初期,T1時(shí)期的特異基因數(shù)量最多,推測塊莖形成初期是最重要的階段,這與馬鈴薯相關(guān)研究結(jié)果[15]一致。T1時(shí)期上調(diào)表達(dá)的DEGs在蛋白質(zhì)磷酸化、蛋白激酶途徑中顯著富集,而蛋白質(zhì)磷酸化作為細(xì)胞信號(hào)系統(tǒng)的重要組成部分,經(jīng)常參與許多生物過程,如植物生長發(fā)育、細(xì)胞分化、刺激反應(yīng)和抗應(yīng)激調(diào)節(jié)[16]。T1時(shí)期的特異基因編碼211個(gè)轉(zhuǎn)錄因子,揭示許多編碼TFs的基因參與塊莖形成。另外T1時(shí)期表達(dá)的大部分基因與激素有關(guān),說明激素與塊莖形成關(guān)系密切。

    3.2 參與山藥塊莖發(fā)育過程的激素相關(guān)基因分析

    植物激素的協(xié)同作用對(duì)塊莖的形成和生長調(diào)節(jié)具有重要意義。GA和ABA是1對(duì)經(jīng)典的植物激素組合,它們的前體都是甲戊龍酸,對(duì)種子休眠和萌發(fā)、根系發(fā)育和開花有拮抗作用,GA和ABA之間的平衡和相互調(diào)節(jié)在植物生長發(fā)育中起著重要作用[17]。關(guān)于GA對(duì)塊莖生長的影響有不同觀點(diǎn),但一般認(rèn)為,GA可以抑制或延緩塊莖的發(fā)育。GA的信號(hào)轉(zhuǎn)導(dǎo)通過環(huán)境和內(nèi)部條件整合GA受體GID1和抑制因子DELLA,形成GA-GID1-DELLA三聚體復(fù)合體,從而調(diào)節(jié)植物的生長發(fā)育[18]。已發(fā)現(xiàn)的DELLA蛋白包括擬南芥GAI、RGA、RGL1、RGL2和RGL3[19]以及小麥RHT[20]、大麥SLN1[21]和水稻SLR1[22]等。本研究發(fā)現(xiàn),編碼DELLA蛋白的GAI在T1時(shí)期高表達(dá),推測塊莖開始形成時(shí),GA就開始受到抑制,導(dǎo)致T2時(shí)期的GA含量迅速下降后,在以后生育期處于較低水平[9]。

    研究發(fā)現(xiàn),NCED3高表達(dá)顯著增加擬南芥的ABA含量[23];ZEP參與擬南芥種子發(fā)育中ABA的合成[24];CYP707A[25]編碼ABA8′羥化酶。本研究發(fā)現(xiàn),1個(gè)NCED基因在T1~T5時(shí)期下調(diào)表達(dá),但在T6期強(qiáng)烈上調(diào)表達(dá);3個(gè)ZEP基因中有2個(gè)基因在早期高表達(dá);3個(gè)CYP707A基因在T1時(shí)期上調(diào)表達(dá),在T2~T6期的表達(dá)量顯著降低。在塊莖形成過程中,調(diào)控ABA的基因的表達(dá)一直較活躍,含量保持上升[9]。對(duì)水稻的研究發(fā)現(xiàn),高ABA含量會(huì)誘導(dǎo)淀粉積累量增加[26],在山藥塊莖發(fā)育后期ABA含量持續(xù)增長,不僅有利于塊莖生長,而且有利于淀粉積累。

    丙二烯氧化物合酶(AOS)是JA合成途徑中的第1個(gè)特異性酶,可以催化氧化丙二烯的合成,12-氧-植物二烯酸還原酶(OPR)是JA合成的關(guān)鍵酶,AOS基因的過度表達(dá)可增加馬鈴薯中JA含量[27]。在本研究中,AOS、OPR在整個(gè)生育期一直高表達(dá),尤其在T1時(shí)期,推測在山藥塊莖的整個(gè)發(fā)育過程中,JA的合成和代謝相對(duì)較強(qiáng),使得塊莖中的JA含量較高[9]。

    共表達(dá)網(wǎng)絡(luò)分析結(jié)果顯示,在T1時(shí)期的特異性模塊中,乙烯應(yīng)答轉(zhuǎn)錄因子ERF036 like、真核細(xì)胞翻譯起始因子等基因的表達(dá)水平較高。ERF轉(zhuǎn)錄因子家族參與多種生物過程,包括激素信號(hào)轉(zhuǎn)導(dǎo)、生長發(fā)育、代謝及對(duì)生物和非生物脅迫的響應(yīng)[28-30],推測在參薯山藥塊莖發(fā)生初期涉及多種激素信號(hào)轉(zhuǎn)導(dǎo)過程。在T6時(shí)期的特異性模塊中,含有2個(gè)轉(zhuǎn)錄因子(WRKY71、WRKY79)編碼基因。在水稻胚、糊粉層細(xì)胞中,WRKY71表達(dá)由ABA誘導(dǎo)并被GA抑制,α-淀粉酶amy32B和赤霉素誘導(dǎo)轉(zhuǎn)錄激活因子GAMYB協(xié)同表達(dá)[31];WRKY79在玉米中過度表達(dá)會(huì)增加茉莉酸和乙烯途徑、活性氧清除中相關(guān)基因及其他相關(guān)基因的表達(dá)量[32],推測在塊莖發(fā)育后期,激素相關(guān)基因的表達(dá)特別活躍。

    4 結(jié)論

    本研究對(duì)山藥轉(zhuǎn)錄組數(shù)據(jù)的分析發(fā)現(xiàn),在山藥塊莖形成初期,尤其是T1時(shí)期為塊莖發(fā)育的重要時(shí)期,階段特異基因數(shù)目最多,上調(diào)的DEGs與細(xì)胞生長、激素均密切相關(guān);塊莖發(fā)育后期激素相關(guān)基因表達(dá)特別活躍,ABA和JA是塊莖生長發(fā)育的重要激素。研究結(jié)果與之前關(guān)于塊莖發(fā)育中激素含量的研究結(jié)果一致,轉(zhuǎn)錄組測序結(jié)合生理測定結(jié)果可以為進(jìn)一步分析山藥塊莖膨大的分子機(jī)制奠定基礎(chǔ)。

    參考文獻(xiàn):

    [1] 韋本輝.中國淮山藥栽培[M]. 北京:中國農(nóng)業(yè)出版社,2013:2-4.

    [2] OSELEBE H O, OKPORIE E O. Evaluation of water yam (Dioscorea alata L.) genotypes for yield and yield components in Abakaliki agro-ecological zone of Nigeria[J]. Agro-Science, 2008, 7(3): 179-185.

    [3] LE B H, CHENG C, BUI A Q, et al. Global analysis of gene activity during Arabidopsis seed development and identification of seed specific transcription factors[J]. Proceedings of the National Academy of Sciences of the United States of America, 2010, 107(18): 8063-8070.

    [4] BELMONTE M F, KIRKBRIDE R C, STONE S L, et al. Comprehensive developmental profiles of gene activity in regions and subregions of the Arabidopsis seed[J]. Proceedings of the National Academy of Sciences of the United States of America, 2013, 110(5): E435-E444.

    [5] ZHOU T, LUO X J, YU C N, et al. Transcriptome analyses provide insights into the expression pattern and sequence similarity of several taxol biosynthesis-related genes in three Taxus species[J]. BMC Plant Biology, 2019, 19(1): 33-42.

    [6] SUN Q, ZHOU G F, CAI Y F, et al. Transcriptome analysis of stem development in the tumourous stem mustard Brassica juncea var. tumida Tsen et Lee by RNA sequencing[J]. BMC Plant Biology, 2012, 12: 53.

    [7] ROUMELIOTIS E, VISSER R G F, BACHEM C W B. A crosstalk of auxin and GA during tuber development[J]. Plant Signaling & Behavior, 2012, 7(10): 1360-1363.

    [8] SARKER D, PANDEY S K, SHARMA S. Cytokinins antagonize the jasmonates action on the regulation of potato (Solanum tuberosum) tuber formation in vitro[J]. Plant Cell, Tissue and Organ Culture, 2006, 87(3): 285-295.

    [9] 張 鉛,蔣 璐,張培通,等. 長江中下游地區(qū)參薯塊莖發(fā)育動(dòng)態(tài)特征分析[J]. 植物生理學(xué)報(bào),2020,56(12):2736-2744.

    [10]GRABHERR M G, HAAS B J, YASSOUR M, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome[J]. Nature Biotechnology, 2011, 29: 644-652.

    [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]ZHAN J P, THAKARE D, MA C, et al. RNA sequencing of laser-capture microdissected compartments of the maize kernel identifies regulatory modules associated with endosperm cell differentiation[J]. The Plant Cell, 2015, 27(3): 513-531.

    [13]LANGFELDER P, HORVATH S. WGCNA: an R package for weighted correlation network analysis[J]. BMC Bioinformatics, 2008, 9(1): 559.

    [14]LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-△△Ctmethod[J]. Methods, 2001, 25(4): 402-408.

    [15]KLOOSTERMAN B, VORST O, HALL R D, et al. Tuber on a chip: differential gene expression during potato tuber development[J]. Plant Biotechnology Journal, 2005, 3(5): 505-519.

    [16]PAWSON T, SCOTT J D. Signaling through scaffold, anchoring, and adaptor proteins[J]. Science, 1997, 278(5346): 2075-2080.

    [17]SHU K, LIU X D, XIE Q, et al. Two faces of one seed: hormonal regulation of dormancy and germination[J]. Molecular Plant, 2016, 9(1): 34-45.

    [18]SUN T P. The molecular mechanism and evolution of the review GAGID1-DELLA signaling module in plants[J]. Current Biology, 2011, 21: 338-345.

    [19]TYLER L, THOMAS S G, HU J H, et al. DELLA proteins and gibberellin regulated seed germination and floral development in Arabidopsis[J]. Plant Physiology, 2004, 135: 1008-1019.

    [20]PENG J, RICHARDS D E, HARTLEY N M, et al. ‘Green Revolution genes encode mutant gibberellin response modu lators[J]. Nature, 1999, 400: 256-261.

    [21]FU X D, RICHARDS D E, AIT-AIT T, et al. Gibberellin-mediated proteasome-dependent degradation of the barley DELLA protein SLN1 repressor[J]. Plant Cell, 2002, 14(12): 3191-3200.

    [22]IKEDA A, UEGUCHI-TANAKA M, SONODA Y, et al. Slender rice, a constitutive gibberellin response mutant, is caused by a null mutation of the SLR1 gene, an ortholog of the height-regulating gene GAI/RGA/RHT/D8[J]. Plant Cell, 2001, 13: 999-1010.

    [23]IUCHI S, KOBAYASHI M, TAJI T, et al. Regulation of drought tolerance by gene manipulation of 9-cis-epoxycarotenoid dioxygenase, a key enzyme in abscisic acid biosynthesis in Arabidopsis[J]. Plant Journal for Cell & Molecular Biology, 2010, 27(4): 325-333.

    [24]SCHWARZ N, ARMBRUSTER U, LVEN T, et al. Tissue-specific accumula tion and regulation of zeaxanthin epoxidase in Arabidopsi reflect the multiple functions of the enzyme in plastids[J]. Plant Cell Physiology, 2015, 56(2): 346-357.

    [25]KROCHKO J E, ABRAMS G D, LOEWEN M K, et al. (+)-Abscisic acid 8-hydroxylase is a cytochrome P450 monooxygenase[J]. Plant Physiol, 1998, 118(3): 849-860.

    [26]WANG Z Q, XU Y J, CHEN T T, et al. Abscisic acid and the key enzymes and genes in sucrose-to-starch conversion in rice spikelets in response to soil drying during grain filling[J]. Planta, 2015, 241(5): 1091-1107.

    [27]PAN Z Q, CAMARA B, GARDNER H W, et al. Aspirin inhibition and acetylation of the plant cytochrome P450, allene oxide synthase, resembles that of animal prostaglandin endoperoxide H synthase[J]. The Journal of Biological Chemistry, 1998, 273(17): 18139-18145.

    [28]CAO F Y, DEFALCO T A, MOEDER W, et al. Arabidopsis ETHYLENE RESPONSE FACTOR 8(ERF8) has dual functions in ABA signaling and immunity[J]. BMC Plant Biol, 2018, 18 (1): 211.

    [29]LIU K, LI Y H, CHEN X N, et al. ERF72 interacts with ARF6 and BZR1 to regulate hypocotyl elongation in Arabidopsis[J]. Journal of Experimental Botany, 2018, 69 (16): 3933-3947.

    [30]JIN Y, PAN W Y, ZHENG X F, et al. OsERF101, an ERF family transcription factor, regulates drought stress response in reproductive tissues[J]. Plant Molecular Biology, 2018, 98(1/2): 51-65.

    [31]XIE Z, ZHANG Z L, ZOU X L, et al. Interactions of two abscisic-acid induced WRKY genes in repressing gibberellin signaling in aleurone cells[J]. The Plant Journal, 2006, 46(2): 231-242.

    [32]FU J Y, LIU Q, WANG C, et al. Zmwrky79 positively regulates maize phytoalexin biosynthetic gene expression and is involved in stress response[J]. Journal of Experimental Botany, 2018, 69(3): 497-510.

    (責(zé)任編輯:徐 艷)

    猜你喜歡
    山藥
    可以飽食的人參 晉臻長山藥
    山藥買粗的還是細(xì)的
    山藥的養(yǎng)生吃法
    治病山藥 美味山藥
    山藥記
    文苑(2020年6期)2020-06-22 08:41:50
    山藥提早上市栽培技術(shù)
    我的“魔法山藥”
    山藥被稱“長壽因子”
    飛人食譜山藥飯
    中國自行車(2018年2期)2018-05-09 07:02:50
    山藥貯藏留種技術(shù)
    午夜av观看不卡| 乱人伦中国视频| 在线观看免费午夜福利视频| 汤姆久久久久久久影院中文字幕| 日韩一区二区三区影片| 美女扒开内裤让男人捅视频| 久久久国产欧美日韩av| 精品亚洲成a人片在线观看| 天天操日日干夜夜撸| 日本五十路高清| 国产精品 欧美亚洲| 免费高清在线观看日韩| 久久久久久久大尺度免费视频| 欧美中文综合在线视频| 9热在线视频观看99| 免费女性裸体啪啪无遮挡网站| 国产高清videossex| 国产片内射在线| 精品少妇内射三级| 免费少妇av软件| 免费av中文字幕在线| 国产深夜福利视频在线观看| 久久久久国产一级毛片高清牌| 日韩 亚洲 欧美在线| 欧美精品人与动牲交sv欧美| 亚洲一区中文字幕在线| 国产成人精品久久久久久| 欧美日韩精品网址| 国产成人av教育| 久久久欧美国产精品| 亚洲,一卡二卡三卡| 这个男人来自地球电影免费观看| 精品一区二区三卡| 91精品伊人久久大香线蕉| 观看av在线不卡| 亚洲国产看品久久| 王馨瑶露胸无遮挡在线观看| 色精品久久人妻99蜜桃| 日本欧美视频一区| 欧美人与性动交α欧美软件| 国产成人av教育| 啦啦啦 在线观看视频| 亚洲国产毛片av蜜桃av| 免费在线观看影片大全网站 | 日韩大码丰满熟妇| 中文字幕制服av| 黄片播放在线免费| 亚洲人成电影观看| 男女无遮挡免费网站观看| 国产97色在线日韩免费| 另类精品久久| 日韩大码丰满熟妇| 我要看黄色一级片免费的| 午夜福利一区二区在线看| 一区二区av电影网| 久久精品国产综合久久久| 欧美黑人精品巨大| 黄色片一级片一级黄色片| 黄色视频在线播放观看不卡| 欧美xxⅹ黑人| 校园人妻丝袜中文字幕| 男女无遮挡免费网站观看| 亚洲欧美成人综合另类久久久| 夜夜骑夜夜射夜夜干| 人人妻人人澡人人爽人人夜夜| 亚洲精品国产av蜜桃| 女人爽到高潮嗷嗷叫在线视频| 国产片内射在线| bbb黄色大片| 色网站视频免费| www.自偷自拍.com| 亚洲成人手机| bbb黄色大片| 狠狠婷婷综合久久久久久88av| av在线老鸭窝| 久久久久网色| 可以免费在线观看a视频的电影网站| 99精国产麻豆久久婷婷| 嫩草影视91久久| 国产高清videossex| 亚洲精品久久久久久婷婷小说| 中文字幕人妻丝袜一区二区| 国产一区二区激情短视频 | 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成色77777| 国产深夜福利视频在线观看| 国产成人精品久久二区二区免费| 女人精品久久久久毛片| 一级毛片黄色毛片免费观看视频| 丁香六月天网| 美国免费a级毛片| 欧美精品av麻豆av| 欧美日韩成人在线一区二区| 91麻豆精品激情在线观看国产 | 少妇被粗大的猛进出69影院| 亚洲熟女毛片儿| 免费在线观看日本一区| 99久久精品国产亚洲精品| 日韩大码丰满熟妇| 丝袜美腿诱惑在线| 两个人免费观看高清视频| 久久人人97超碰香蕉20202| 欧美国产精品va在线观看不卡| 天堂8中文在线网| av又黄又爽大尺度在线免费看| 一区二区三区精品91| 90打野战视频偷拍视频| 蜜桃在线观看..| 无遮挡黄片免费观看| 丰满饥渴人妻一区二区三| 欧美人与性动交α欧美软件| av国产精品久久久久影院| √禁漫天堂资源中文www| 如日韩欧美国产精品一区二区三区| 狠狠婷婷综合久久久久久88av| 国产高清不卡午夜福利| 每晚都被弄得嗷嗷叫到高潮| 2021少妇久久久久久久久久久| 国产精品三级大全| 18在线观看网站| 日韩熟女老妇一区二区性免费视频| 成人午夜精彩视频在线观看| 欧美 日韩 精品 国产| 99久久99久久久精品蜜桃| 夫妻午夜视频| 亚洲欧美成人综合另类久久久| 久久天躁狠狠躁夜夜2o2o | 亚洲成国产人片在线观看| 免费在线观看黄色视频的| 成年人黄色毛片网站| 亚洲黑人精品在线| 亚洲欧美一区二区三区黑人| 成年人免费黄色播放视频| 国产不卡av网站在线观看| 丝袜喷水一区| 狠狠精品人妻久久久久久综合| 免费看av在线观看网站| 亚洲欧美一区二区三区国产| 99国产精品一区二区三区| 午夜免费鲁丝| 多毛熟女@视频| 亚洲国产欧美在线一区| 精品国产一区二区久久| 亚洲色图 男人天堂 中文字幕| 国产精品成人在线| 久久人妻熟女aⅴ| 日韩熟女老妇一区二区性免费视频| 成年人免费黄色播放视频| 99国产精品免费福利视频| 高清视频免费观看一区二区| 两个人免费观看高清视频| 国产xxxxx性猛交| 女人高潮潮喷娇喘18禁视频| 乱人伦中国视频| 天天躁夜夜躁狠狠躁躁| 狠狠精品人妻久久久久久综合| 精品国产超薄肉色丝袜足j| 丝袜脚勾引网站| 最近手机中文字幕大全| 夜夜骑夜夜射夜夜干| 国产精品久久久久成人av| 色视频在线一区二区三区| 久久久精品94久久精品| 国产一区亚洲一区在线观看| 国产国语露脸激情在线看| 亚洲av成人精品一二三区| 婷婷色综合大香蕉| 如日韩欧美国产精品一区二区三区| 如日韩欧美国产精品一区二区三区| 欧美精品人与动牲交sv欧美| 免费女性裸体啪啪无遮挡网站| 国产亚洲精品久久久久5区| 国产在线免费精品| 精品人妻熟女毛片av久久网站| 久久久欧美国产精品| 巨乳人妻的诱惑在线观看| 99久久人妻综合| 亚洲成人免费电影在线观看 | 国产欧美日韩一区二区三 | 色网站视频免费| 午夜免费观看性视频| 久久精品aⅴ一区二区三区四区| 久久狼人影院| 国产真人三级小视频在线观看| 丁香六月天网| 新久久久久国产一级毛片| 欧美av亚洲av综合av国产av| 亚洲图色成人| 99久久人妻综合| 极品少妇高潮喷水抽搐| 亚洲精品一区蜜桃| 老熟女久久久| 乱人伦中国视频| 无遮挡黄片免费观看| videos熟女内射| 日本午夜av视频| 国产成人欧美| 亚洲一卡2卡3卡4卡5卡精品中文| 中文字幕制服av| 亚洲精品av麻豆狂野| 美女福利国产在线| 国产精品99久久99久久久不卡| 悠悠久久av| 国产无遮挡羞羞视频在线观看| 国产亚洲一区二区精品| 亚洲精品国产区一区二| 一本综合久久免费| 国产片内射在线| 青春草亚洲视频在线观看| 国产一区亚洲一区在线观看| 国产精品国产三级专区第一集| 国产av精品麻豆| 国产精品.久久久| 老司机影院毛片| 最近最新中文字幕大全免费视频 | 国产成人免费无遮挡视频| 午夜福利视频精品| 精品少妇内射三级| 91九色精品人成在线观看| 一本大道久久a久久精品| 中文字幕精品免费在线观看视频| 9热在线视频观看99| 欧美日韩av久久| 热re99久久精品国产66热6| 欧美亚洲日本最大视频资源| 在线亚洲精品国产二区图片欧美| 777久久人妻少妇嫩草av网站| 欧美精品亚洲一区二区| 亚洲欧美激情在线| 老司机靠b影院| 亚洲少妇的诱惑av| 久久精品aⅴ一区二区三区四区| 一区二区日韩欧美中文字幕| 老汉色av国产亚洲站长工具| 黄色视频不卡| 久久国产精品大桥未久av| 国产免费福利视频在线观看| 免费高清在线观看视频在线观看| 欧美成狂野欧美在线观看| 亚洲人成电影免费在线| 日韩中文字幕视频在线看片| 青春草亚洲视频在线观看| 久久国产亚洲av麻豆专区| netflix在线观看网站| 99国产精品一区二区三区| 韩国高清视频一区二区三区| 乱人伦中国视频| avwww免费| 精品人妻熟女毛片av久久网站| 亚洲欧美一区二区三区黑人| 电影成人av| 国产欧美日韩综合在线一区二区| 成年女人毛片免费观看观看9 | 岛国毛片在线播放| 丝袜美腿诱惑在线| 大香蕉久久成人网| 欧美日韩一级在线毛片| 首页视频小说图片口味搜索 | 久久精品人人爽人人爽视色| 免费在线观看黄色视频的| 男男h啪啪无遮挡| 美女国产高潮福利片在线看| 成人午夜精彩视频在线观看| 亚洲精品自拍成人| 国产国语露脸激情在线看| 在线av久久热| 一级毛片电影观看| 国精品久久久久久国模美| av国产久精品久网站免费入址| 热99国产精品久久久久久7| 不卡av一区二区三区| videos熟女内射| 国产三级黄色录像| 国产激情久久老熟女| 一边亲一边摸免费视频| 欧美日韩亚洲国产一区二区在线观看 | 男女床上黄色一级片免费看| 看免费成人av毛片| 日韩熟女老妇一区二区性免费视频| 久久精品国产亚洲av高清一级| 一区二区三区精品91| av片东京热男人的天堂| 成人国语在线视频| 久久影院123| 久久久久网色| 国产午夜精品一二区理论片| 男女边摸边吃奶| 国产片内射在线| 亚洲av片天天在线观看| 不卡av一区二区三区| 免费在线观看影片大全网站 | videos熟女内射| 国产伦理片在线播放av一区| 女性生殖器流出的白浆| 成在线人永久免费视频| 只有这里有精品99| 亚洲国产av影院在线观看| 多毛熟女@视频| 无限看片的www在线观看| 9热在线视频观看99| 丝袜在线中文字幕| 日韩电影二区| 欧美激情高清一区二区三区| 伊人亚洲综合成人网| 日本色播在线视频| 七月丁香在线播放| 亚洲精品国产色婷婷电影| 国产成人一区二区三区免费视频网站 | 99久久人妻综合| 考比视频在线观看| 欧美激情高清一区二区三区| 成年人黄色毛片网站| 亚洲综合色网址| 亚洲欧美日韩高清在线视频 | 亚洲av在线观看美女高潮| 国产主播在线观看一区二区 | 国产成人一区二区三区免费视频网站 | 男女无遮挡免费网站观看| 亚洲av美国av| 99国产精品一区二区蜜桃av | 欧美日韩亚洲高清精品| 91成人精品电影| 欧美黑人欧美精品刺激| 狠狠婷婷综合久久久久久88av| a级毛片在线看网站| 欧美人与性动交α欧美精品济南到| 少妇的丰满在线观看| 午夜免费鲁丝| 久久人人爽人人片av| 成人18禁高潮啪啪吃奶动态图| 亚洲精品日本国产第一区| 男女午夜视频在线观看| 高清视频免费观看一区二区| 人成视频在线观看免费观看| 精品一区二区三区四区五区乱码 | 国产成人av激情在线播放| 曰老女人黄片| 少妇精品久久久久久久| 成人免费观看视频高清| 国产色视频综合| 亚洲欧洲精品一区二区精品久久久| 日本欧美国产在线视频| 午夜福利视频在线观看免费| 国产一区二区激情短视频 | 老汉色av国产亚洲站长工具| 欧美激情高清一区二区三区| 国产欧美日韩综合在线一区二区| 亚洲av电影在线观看一区二区三区| 午夜激情av网站| 午夜福利视频在线观看免费| 成在线人永久免费视频| 免费人妻精品一区二区三区视频| 日本猛色少妇xxxxx猛交久久| 亚洲自偷自拍图片 自拍| 欧美+亚洲+日韩+国产| 国产午夜精品一二区理论片| 欧美精品人与动牲交sv欧美| 一级片'在线观看视频| 免费在线观看视频国产中文字幕亚洲 | 亚洲成国产人片在线观看| 亚洲国产最新在线播放| 99久久人妻综合| 美女脱内裤让男人舔精品视频| 男女边摸边吃奶| 免费av中文字幕在线| 国产精品免费视频内射| 熟女av电影| 如日韩欧美国产精品一区二区三区| 日本av手机在线免费观看| 亚洲av综合色区一区| 亚洲久久久国产精品| 国产成人精品久久久久久| 色网站视频免费| 日韩一卡2卡3卡4卡2021年| 99国产精品99久久久久| 黑人巨大精品欧美一区二区蜜桃| 午夜福利,免费看| 国产色视频综合| 久久青草综合色| 亚洲国产成人一精品久久久| 中文字幕人妻丝袜一区二区| 国产一区二区三区av在线| 亚洲情色 制服丝袜| 国产日韩欧美视频二区| 亚洲久久久国产精品| 亚洲五月婷婷丁香| 亚洲av电影在线观看一区二区三区| 亚洲一区中文字幕在线| 亚洲国产看品久久| 亚洲国产欧美一区二区综合| 成年女人毛片免费观看观看9 | 免费女性裸体啪啪无遮挡网站| 色综合欧美亚洲国产小说| a级毛片黄视频| 80岁老熟妇乱子伦牲交| 亚洲精品久久久久久婷婷小说| 国产片内射在线| 97人妻天天添夜夜摸| 成年动漫av网址| 久久九九热精品免费| 日日摸夜夜添夜夜爱| 日韩制服丝袜自拍偷拍| 欧美激情 高清一区二区三区| 1024视频免费在线观看| www.精华液| 天堂俺去俺来也www色官网| 国产精品二区激情视频| 桃花免费在线播放| 成人手机av| 精品第一国产精品| 亚洲精品在线美女| 久久鲁丝午夜福利片| 欧美日韩一级在线毛片| 欧美乱码精品一区二区三区| 亚洲 国产 在线| 交换朋友夫妻互换小说| 两个人看的免费小视频| 美女视频免费永久观看网站| 国产高清不卡午夜福利| 久久久久视频综合| 国产老妇伦熟女老妇高清| 国产不卡av网站在线观看| 别揉我奶头~嗯~啊~动态视频 | 亚洲情色 制服丝袜| 久久久久久久大尺度免费视频| 国产精品一区二区在线观看99| 亚洲精品乱久久久久久| 手机成人av网站| 1024香蕉在线观看| 国产熟女午夜一区二区三区| 91精品伊人久久大香线蕉| 国产精品久久久av美女十八| 日韩大码丰满熟妇| 在线观看免费高清a一片| 国产野战对白在线观看| 亚洲国产欧美一区二区综合| 波野结衣二区三区在线| 99国产精品一区二区三区| 亚洲精品国产av蜜桃| 中文字幕另类日韩欧美亚洲嫩草| 操美女的视频在线观看| 免费在线观看视频国产中文字幕亚洲 | 一区二区日韩欧美中文字幕| 国产免费福利视频在线观看| 精品一区在线观看国产| 亚洲精品国产区一区二| 成人亚洲欧美一区二区av| 午夜激情久久久久久久| 国产精品99久久99久久久不卡| 欧美激情 高清一区二区三区| 久9热在线精品视频| 国产免费视频播放在线视频| 久久久国产欧美日韩av| 免费看十八禁软件| 国产成人av激情在线播放| 精品人妻1区二区| 亚洲av国产av综合av卡| 国产成人免费无遮挡视频| 老熟女久久久| 一级片免费观看大全| 免费在线观看日本一区| 亚洲第一av免费看| 亚洲国产欧美一区二区综合| 91老司机精品| 悠悠久久av| 亚洲国产毛片av蜜桃av| 亚洲精品美女久久av网站| 久久久欧美国产精品| 大香蕉久久网| 91麻豆精品激情在线观看国产 | 国产男女超爽视频在线观看| 免费不卡黄色视频| 老司机亚洲免费影院| 久久久久久久精品精品| 晚上一个人看的免费电影| 久久女婷五月综合色啪小说| 大片电影免费在线观看免费| 亚洲中文字幕日韩| 九草在线视频观看| 中文精品一卡2卡3卡4更新| 精品少妇一区二区三区视频日本电影| 九草在线视频观看| 中文字幕色久视频| 国产亚洲av高清不卡| 一二三四社区在线视频社区8| 各种免费的搞黄视频| 欧美日韩精品网址| 国产在线视频一区二区| 各种免费的搞黄视频| 啦啦啦啦在线视频资源| 国产亚洲精品久久久久5区| 国产精品三级大全| 欧美日韩亚洲高清精品| 成在线人永久免费视频| 一级片'在线观看视频| 亚洲精品中文字幕在线视频| 成人亚洲欧美一区二区av| 激情五月婷婷亚洲| 国产精品久久久久久精品古装| 成年美女黄网站色视频大全免费| 丰满迷人的少妇在线观看| 99久久人妻综合| a级毛片在线看网站| 国产精品一区二区免费欧美 | 亚洲午夜精品一区,二区,三区| 99国产精品免费福利视频| 亚洲av男天堂| 超碰97精品在线观看| 免费女性裸体啪啪无遮挡网站| 高清不卡的av网站| 在线观看一区二区三区激情| 只有这里有精品99| 日本欧美国产在线视频| www.999成人在线观看| 国产男女内射视频| 晚上一个人看的免费电影| 亚洲av欧美aⅴ国产| 午夜91福利影院| videosex国产| 久久99热这里只频精品6学生| 成年动漫av网址| 桃花免费在线播放| 精品少妇一区二区三区视频日本电影| 亚洲av日韩精品久久久久久密 | 亚洲中文日韩欧美视频| 久久精品aⅴ一区二区三区四区| 日本vs欧美在线观看视频| 国产亚洲一区二区精品| 久久精品人人爽人人爽视色| 少妇人妻 视频| 老司机影院成人| 国产一区二区激情短视频 | 99精国产麻豆久久婷婷| av天堂久久9| 999精品在线视频| 女性生殖器流出的白浆| √禁漫天堂资源中文www| 人人妻,人人澡人人爽秒播 | 欧美另类一区| 久久久国产欧美日韩av| 亚洲激情五月婷婷啪啪| 9191精品国产免费久久| 国产一卡二卡三卡精品| 日韩大片免费观看网站| 欧美日韩亚洲高清精品| 国产成人影院久久av| 久久久久网色| 美女脱内裤让男人舔精品视频| 激情五月婷婷亚洲| 精品一区二区三卡| 亚洲人成电影免费在线| 黄网站色视频无遮挡免费观看| 纵有疾风起免费观看全集完整版| 黄色毛片三级朝国网站| 人妻一区二区av| 亚洲国产欧美在线一区| 晚上一个人看的免费电影| 婷婷色麻豆天堂久久| 欧美精品啪啪一区二区三区 | 国产精品久久久久久精品电影小说| 夫妻性生交免费视频一级片| 久久99热这里只频精品6学生| 欧美黄色片欧美黄色片| 十分钟在线观看高清视频www| 自线自在国产av| 亚洲人成电影免费在线| 欧美日韩黄片免| 人成视频在线观看免费观看| 欧美 亚洲 国产 日韩一| 夫妻性生交免费视频一级片| 啦啦啦在线免费观看视频4| 人妻一区二区av| 在线观看国产h片| 精品视频人人做人人爽| 最新在线观看一区二区三区 | 热99久久久久精品小说推荐| 亚洲成国产人片在线观看| 伊人久久大香线蕉亚洲五| 国产黄色免费在线视频| 亚洲欧美清纯卡通| 国产男女超爽视频在线观看| 最近手机中文字幕大全| 免费观看人在逋| 手机成人av网站| 又大又爽又粗| 激情视频va一区二区三区| 又大又爽又粗| 亚洲成人免费电影在线观看 | 午夜日韩欧美国产| 亚洲国产日韩一区二区| 日本av手机在线免费观看| 精品人妻熟女毛片av久久网站| 国产视频首页在线观看| 亚洲视频免费观看视频| 777米奇影视久久| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产av新网站| 黄色 视频免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产无遮挡羞羞视频在线观看| 母亲3免费完整高清在线观看| 欧美国产精品va在线观看不卡| 国产免费福利视频在线观看| 久久精品亚洲熟妇少妇任你| 久久中文字幕一级| 国语对白做爰xxxⅹ性视频网站| 国产精品av久久久久免费| 丰满少妇做爰视频| 精品国产一区二区三区久久久樱花| 国语对白做爰xxxⅹ性视频网站| 大陆偷拍与自拍| 最近手机中文字幕大全| 国产伦人伦偷精品视频| 国语对白做爰xxxⅹ性视频网站| 国产精品av久久久久免费|