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

    玉米自交系響應高溫、干旱脅迫的關(guān)鍵基因及通路

    2021-03-25 13:50:00姚啟倫霍仕平張俊軍
    江蘇農(nóng)業(yè)學報 2021年1期
    關(guān)鍵詞:差異表達基因轉(zhuǎn)錄組干旱脅迫

    姚啟倫 霍仕平 張俊軍

    摘要: 以4個不同的玉米自交系為材料,對高溫、干旱處理后的苗期植株進行轉(zhuǎn)錄組測序。玉米自交系響應高溫和干旱脅迫的差異表達基因(DEGs)分別為6 966和6 272個,在高溫和干旱脅迫下4個玉米自交系相同的DEGs分別是705和871個。同時響應高溫和干旱的DEGs有100個。在耐旱、耐熱性強的玉米自交系中鑒定出18個特異的DEGs,其中鋅指轉(zhuǎn)錄因子、WRKY轉(zhuǎn)錄因子、GT轉(zhuǎn)錄因子和B2熱激轉(zhuǎn)錄因子在脅迫響應中發(fā)揮關(guān)鍵的調(diào)控作用。KEGG通路分析結(jié)果表明,耐旱、耐熱性強玉米自交響應高溫干旱脅迫的DEGs富集在生物學過程、分子功能、代謝過程、遍在蛋白代謝和氮代謝途徑5條通路。熱帶、亞熱帶玉米種質(zhì)的耐旱、耐熱性強于溫帶玉米種質(zhì),可在熱帶、亞熱帶玉米種質(zhì)中有效篩選耐旱、耐熱基因。

    關(guān)鍵詞: 玉米自交系;高溫、干旱脅迫;轉(zhuǎn)錄組;差異表達基因;通路

    中圖分類號: S513 文獻標識碼: A 文章編號: 1000-4440(2021)01-0029-09

    Key genes and pathways of maize inbred lines responding to heat and drought stress

    YAO Qi-lun1, HUO Shi-ping2, ZHANG Jun-jun3

    (1.School of Advanced Agriculture and Bioengineering, Yangtze Normal University, Chongqing 408100, China;2.Chongqing Three Gorges Acodemy of Agricultural Science, Chongqing 404100, China;3.Heilongjiang Hetian Fengze Agriculture Science and Technology Development Co., Ltd., Harbin? 067500, China)

    Abstract: Transcriptome sequencing of the four different maize inbred lines subjected to high temperatures and water deficits, was performed at the stage of seedling. In total, 6 966 heat-responsive and 6 272 drought-responsive differentially expressed genes(DEGs) were identified in the four maize lines. In addition, 705 and 871 DEGs were identified as being commonly associated with heat and drought stress, respectively. There were 100 DEGs responding to both heat and drought stress. Eighteen DEGs were identified in maize inbred lines with strong drought tolerance and heat tolerance. The zinc finger transcription factor, WRKY transcription factor, GT transcription factor and B2 heat shock transcription factor played key regulatory roles in stress response. Analyses of KEGG pathway enrichment showed that the pathways such as biological process, molecular function, metabolic process, protein ubiquitination metabolism and nitrogen metabolism were the most highly enriched in the maize lines tolerant to heat and drought. It can be concluded that tropical and subtropical maize is tolerant to heat and drought, implying that it is plausible to screen goal genes in the tropical and subtropical germplasm.

    Key words: maize inbred lines;heat and drought stress;transcriptome;differentially expressed genes(DEGs);pathways

    高溫和干旱是抑制玉米生長發(fā)育、制約玉米產(chǎn)量的主要非生物脅迫因子[1]。在干旱、半干旱區(qū),玉米往往受到高溫干旱復合脅迫的不利影響[2-3]。全球有43 %的干旱、半干旱耕地,中國干旱、半干旱土地面積約占全國土地面積的近50 %[4-5]。隨著溫室效應加劇,預計到本世紀末地表氣溫將上升1~5 ℃,1980-2008年,高溫脅迫導致全球玉米產(chǎn)量下降3.8%[6-7]。因此,篩選玉米耐旱、耐熱基因,培育耐旱、耐熱品種一直是玉米育種界面臨的重大課題。迄今關(guān)于高溫、干旱脅迫對玉米植株形態(tài)、生理生化和產(chǎn)量品質(zhì)不利影響的研究已有較多文獻報道。研究結(jié)果表明,高溫、干旱脅迫下,玉米葉片伸長速率減慢,葉面積減小,莖稈變細[8-9];高溫、干旱會減少玉米葉綠素含量、擾亂光合系統(tǒng),抑制光合作用[10-12];高溫、干旱會降低ADP-葡萄糖焦磷酸化酶、淀粉合成酶及分支酶活性,從而抑制碳水化合物合成[13-14];玉米灌漿期,高溫、干旱打破植株體內(nèi)激素間的平衡,減緩細胞分裂速率,縮小籽粒庫容,導致產(chǎn)量和品質(zhì)下降[15]。利用高通量轉(zhuǎn)錄組測序(RNA-seq)技術(shù)在分子水平解析植物對逆境脅迫的響應是生物技術(shù)研究的熱點。通過分析生物細胞中轉(zhuǎn)錄組的變化規(guī)律,可在分子水平探究特定生物學過程中基因組轉(zhuǎn)錄調(diào)控機制[16-17]。玉米在傳播過程中逐漸分化形成了適應其栽培環(huán)境的熱帶、亞熱帶和溫帶種質(zhì),在分子水平研究不同玉米種質(zhì)的逆境響應機制,可有效指導耐逆玉米基因的發(fā)掘與利用[18]。盡管轉(zhuǎn)錄組測序(RNA-seq)技術(shù)已成功應用于玉米逆境脅迫研究[19-20],但在逆境脅迫研究中往往只涉及單個脅迫因子以及單一玉米自交系和雜交種。本研究選取熱帶、亞熱帶和溫帶種質(zhì),且具有耐旱、耐熱性差異的玉米自交系為供試材料,設(shè)置高溫、干旱2個脅迫處理,結(jié)合高溫、干旱脅迫下不同玉米種質(zhì)的生理特性分析,在轉(zhuǎn)錄組水平分析不同種質(zhì)對高溫、干旱脅迫的分子響應機制,以豐富植物逆境脅迫研究的數(shù)據(jù)資料,并為玉米耐旱、耐熱基因的篩選提供理論支撐。

    1 材料與方法

    1.1 試驗材料

    以長江師范學院選育的4個玉米自交系Suwan-3、Suwan-10、Cim-5和Cim-16為供試材料,Suwan-3和Suwan-10分別以熱帶和亞熱帶玉米地方品種為基礎(chǔ)材料選育而成, Cim-5和Cim-16以國際玉米小麥改良中心(CIMMYT)的溫帶玉米群體為基礎(chǔ)材料選育而成。2016-2018年田間鑒定玉米自交系的耐旱、耐熱性,參照張麗梅等[21]的方法,分別計算各自交系經(jīng)濟產(chǎn)量的耐旱指數(shù)(DTI)和耐熱指數(shù)(HTI),即干旱(或高溫)脅迫下經(jīng)濟產(chǎn)量與正常條件(CK)下經(jīng)濟產(chǎn)量的比值,根據(jù)經(jīng)濟產(chǎn)量的DTI和HTI篩選出具有明顯耐旱耐熱性差異的熱帶、亞熱帶和溫帶玉米自交系Suwan-3、Suwan-10、Cim-5和Cim-16。 熱帶、亞熱帶玉米自交系Suwan-3、Suwan-10耐旱耐熱性強,經(jīng)濟產(chǎn)量的DTI分別為0.86和0.83,HTI分別為0.91和0.93;溫帶玉米自交系Cim-5和Cim-16耐旱耐熱性弱,經(jīng)濟產(chǎn)量的DTI分別為0.57和0.51,HTI分別為0.55和0.49。Suwan-3和Suwan-10作為耐旱、耐熱性強的玉米自交系(基因型)試驗材料,Cim-5和Cim-16作為耐旱、耐熱性弱的玉米自交系(基因型)試驗材料。

    1.2 脅迫處理

    高溫、干旱脅迫處理在人工氣候箱中進行,4個玉米自交系種子浸種催芽后在控制條件下培養(yǎng),培養(yǎng)溫度29 ℃(晝)/23 ℃(夜),相對濕度70%,每天光照12 h。待幼苗生長至5葉1心,4個自交系各種植18株,其中6株作對照、6株作高溫脅迫處理、6株作干旱處理,共72株(每株1個營養(yǎng)缽)。作對照的4株自交系幼苗置于上述人工氣候箱內(nèi)繼續(xù)生長;作高溫處理的4個自交系幼苗在40 ℃(晝)/34 ℃(夜)、相對濕度70%、每天光照12 h的高溫條件下生長7 d;作干旱處理的4個自交系幼苗在溫度29 ℃(晝)/23 ℃(夜)、相對濕度20%、每天光照12 h的控水條件下生長7 d。取每個自交系各個處理3株植株的葉片,葉片剪碎后分成2份分別用作生理指標測定和mRNA轉(zhuǎn)錄組分析。

    1.3 生理指標測定

    利用便攜式光合儀(LI-6400XT)室外檢測待測植株旗葉的葉綠素熒光參數(shù)值:光反應最大熒光(F′m)、光反應最小熒光(F′o)、暗適應最大熒光(Fm)、暗適應最小熒光(Fo)和穩(wěn)態(tài)熒光(Fs);基于F′m、F′o、Fm、Fo和Fs,計算PS II潛在光量子產(chǎn)量(Fv/Fo)、PS II最大光化學量子產(chǎn)量(Fv/Fm)、PS II 有效光化學量子產(chǎn)量(ΦPSII)和光化學熒光淬滅系數(shù)(QP) [22-27]。采用李小芳等[28]的方法室內(nèi)檢測葉綠素含量、丙二醛(MDA)含量、過氧化物酶(POD)活性和超氧化物歧化酶(SOD)活性。數(shù)據(jù)分析采用Microsoft Excel 2010和DPS軟件。

    1.4 葉片RNA提取

    采用TRIzol(Invitrogen, Carisbad, CA, USA)法提取對照、高溫處理和干旱處理下的4個自交系的葉片RNA,用無RNA的DNA酶Ⅰ去除樣品中的DNA,同時檢測樣品中RNA濃度及其完整性(要求RNA完整值大于6.5),共構(gòu)建12個mRNA文庫。

    1.5 轉(zhuǎn)錄組文庫構(gòu)建與測序

    經(jīng)片段化后的mRNA反轉(zhuǎn)錄合成cDNA,篩選300~500 bp大小的cDNA片段進行PCR擴增,轉(zhuǎn)錄組文庫質(zhì)量經(jīng)Agilent 2100 Bioanalyzer(Agilent Technologies Santa Clata, CA, USA)檢測合格后,在Illumina HiSeq 2500測序平臺(海元莘生物醫(yī)藥技術(shù)有限公司)上進行。測序的基本原理是邊合成邊測序,采用Illumina novaseq 6000測序,雙端測序讀數(shù)150 bp、讀長300 bp。在測序的流式細胞中加入4種熒光標記的 dNTP、DNA 聚合酶以及接頭引物進行擴增,在每一個測序簇延伸互補鏈時,每加入一個被熒光標記的 dNTP 就能釋放出相對應的熒光,測序儀通過捕獲熒光信號,并通過計算機軟件將光信號轉(zhuǎn)化為測序峰,從而獲得堿基序列信息。

    1.6 轉(zhuǎn)錄組序列分析

    過濾掉接頭污染、低質(zhì)量(Q<30)和位置堿基氮含量> 5%的初始序列,獲取過濾后序列作后續(xù)分析。過濾后序列比對到玉米B73參考基因組(https://download.maizegdb.rg/Zm-B73-REFERENCE-GRAMENE-4.0/Zm00001d.2.genomic. fa),利用RSEM(0.9.3版本)軟件作基因表達量估算。以P<0.05和|log2FC|≥1為標準,用edgeR軟件作差異表達基因(DEGs)分析 [29]。在NCBI的玉米基因數(shù)據(jù)庫中獲取DEGs的功能信息,基于內(nèi)部perl腳本在maizeGDB網(wǎng)站上獲取DEGs的GO注釋。

    1.7 差異表達基因(DEGs)的定量驗證

    以持家基因Actin為標準,選擇8個差異表達基因進行實時定量(qRT-PCR)分析,利用Primer 5.0軟件在線設(shè)計qRT-PCR反應引物(表1),qRT-PCR反應程序按照ABI2100反應體系(Applied Biosystems, CA, USA)的操作要求進行。各個試驗設(shè)置3次重復,采用相對定量(2-△△Ct)法[30]估算基因表達量。

    2 結(jié)果與分析

    2.1 玉米自交系對高溫、干旱的生理響應

    由表2可知,供試材料在單一脅迫因子高溫或干旱條件下,玉米自交系Cim-5和Cim-16 的過氧化物酶、超氧化物歧化酶活性以及丙二醛含量較對照顯著增加, 4個葉綠素熒光參數(shù)值(Fv/Fo、Fv/Fm、ΦPS II和QP)較對照顯著下降;熱帶、亞熱帶玉米自交系Suwan-3和Suwan-10 除丙二醛含量較對照顯著增加外,過氧化物酶、超氧化物歧化酶活性及4個葉綠素熒光參數(shù)值則無顯著變化,這表明在生理上熱帶、亞熱帶玉米種質(zhì)的耐旱、耐熱性強于溫帶種質(zhì)。

    2.2 玉米自交系響應高溫、干旱脅迫的轉(zhuǎn)錄組分析

    由表3可知,測序得到的初始序列經(jīng)過除去接頭污染、低質(zhì)量(Q<30)和位置堿基氮含量> 5%的過濾處理后,4個自交系在對照、高溫和干旱處理條件下的過濾后序列變幅為44 560 000~54 040 000條,其過濾后序列占初始序列的比例在97.5 %以上,說明轉(zhuǎn)錄組文庫質(zhì)量良好。

    在4個玉米自交系轉(zhuǎn)錄組中共檢測到6 966個響應高溫脅迫的差異表達基因,自交系Suwan-3、Suwan-10、Cim-5和Cim-16分別有1 444、1 281、2 509和1 732個差異表達基因(圖1A);在4個玉米自交系轉(zhuǎn)錄組中共檢測到6 272個響應干旱脅迫的DEGs,自交系Suwan-3、Suwan-10、Cim-5和Cim-16分別有1 843、1 841、1 414和1 174 個差異表達基因(圖1B);4個玉米自交系在高溫、干旱和高溫干旱脅迫下相同的DEGs分別是705、871和100個, 分別占差異表達基因總數(shù)的42%、52%和6%(圖1C)??傮w上,溫帶玉米自交系響應高溫、干旱脅迫的差異表達基因數(shù)(6 829)高于熱帶、亞熱帶玉米自交系(6 409)。這表明在分子水平,高溫、干旱脅迫對溫帶玉米自交系的影響更大。

    2.3 熱帶、亞熱帶玉米自交系在高溫、干旱脅迫下特異的差異表達基因

    為探明玉米自交系中的耐旱、耐熱關(guān)鍵基因,用在FPKM(1 000 000條序列中單個基因的數(shù)目)中比對上的Reads數(shù)目來度量基因表達水平,在P<0.05和|log2FC|≥1標準下,在耐旱、耐熱性強的熱帶、亞熱帶玉米自交系Suwan-3和Suwan-10 中篩選特異的差異表達基因。由表4可知,在耐旱、耐熱性強的玉米自交系Suwan-3和Suwan-10中有18個特異的差異表達基因同時響應高溫和干旱脅迫。這18個差異表達基因極顯著差異表達,其中鋅指轉(zhuǎn)錄因子、WRKY轉(zhuǎn)錄因子、GT轉(zhuǎn)錄因子和B2熱激轉(zhuǎn)錄因子一致上調(diào),這4個轉(zhuǎn)錄因子在玉米響應高溫、干旱脅迫中可能發(fā)揮關(guān)鍵的遺傳調(diào)控作用。

    采用qRT-PCR技術(shù),選擇8個響應高溫、干旱的關(guān)鍵基因Zm00001d053913、Zm00001d018803、Zm00001d035139、Zm00001d002801、Zm00001d012325、Zm00001d015743、Zm00001d017597和 Zm00001d020614,驗證轉(zhuǎn)錄組中差異表達基因的表達結(jié)果,盡管基因Zm00001d012325、Zm00001d015743、Zm00001d017597和 Zm00001d020614在轉(zhuǎn)錄組測序中差異表達不顯著,但這4個基因在植物響應高溫、干旱中發(fā)揮重要作用。由圖2可見,基因Zm00001d053913、Zm00001d002801、Zm00001d015743和Zm00001d017597上調(diào)表達,而基因Zm00001d018803、Zm00001d035139、Zm00001d012325和 Zm00001d020614下調(diào)表達,實時定量試驗數(shù)據(jù)與轉(zhuǎn)錄組測序數(shù)據(jù)(RNA-Seq prolifing)結(jié)果一致?;騔m00001d035139 (編碼MA3結(jié)構(gòu)域蛋白)在轉(zhuǎn)錄組測序中上調(diào)表達,而在差異表達基因的定量驗證中下調(diào)表達,這可能是在基因的定量驗證中以持家基因Actin作參照的結(jié)果。

    2.4 玉米自交系差異表達基因的GO注釋

    通過對差異表達基因的GO(Gene Ontology)分析,可確定與差異表達基因相關(guān)的基因本體項(GO terms),追溯差異表達基因的GO分類條目,揭示差異表達基因與其功能調(diào)節(jié)關(guān)系。根據(jù)Kim等[29]的方法,顯示4個玉米自交系差異表達基因的GO分類注釋結(jié)果,由圖3中的A2 Vs A1和B2 Vs B1可見,熱帶、亞熱帶玉米自交系Suwan-3和Suwan-10響應高溫脅迫的差異表達基因在腺嘌呤核苷結(jié)合、腺嘌呤脫氧核苷結(jié)合、脫氧核苷結(jié)合、核苷結(jié)合、細胞大分子代謝、細胞內(nèi)細胞器部分(Intracellur organelle part)、細胞器部分(Organelle part)、初級代謝過程、分子功能、催化活性、大分子代謝途徑和生物學過程12個基因本體項同時上調(diào)(A2 Vs A1和B2 Vs B1);響應干旱脅迫的差異表達基因在生物學過程、大分子代謝途徑、分子功能、代謝途徑、質(zhì)膜、脅迫響應、蛋白質(zhì)代謝途徑以及堿基、核苷、核苷酸和核酸代謝途徑8個基因本體項共富集(A3 Vs A1和B3 Vs B1);熱帶、亞熱帶和溫帶玉米自交系響應高溫、干旱脅迫的差異表達基因僅在分子功能和生物學過程2個基因本體項同時上調(diào),這表明高溫、干旱脅迫主要影響玉米分子功能和生物學過程。

    2.5 熱帶、亞熱帶玉米自交系差異表達基因的KEGG通路

    KEGG通路分析結(jié)果(圖4)表明,高溫、干旱脅迫下自交系Suwan-3和Suwan-10有30個通路受到調(diào)控,涉及16個生物學途徑、9個細胞組分和5個生物學功能通路,生物學過程、分子功能、代謝過程、遍在蛋白和氮代謝途徑5條通路上有大量差異表達基因富集。顯著性檢驗閥值(P值)顯示,大量極顯著的差異表達基因富集在分子功能、代謝過程、遍在蛋白代謝、氮代謝途徑和生物學調(diào)節(jié)通路。由此推測這些通路在玉米響應高溫、干旱脅迫中發(fā)揮主要作用。

    3 討論

    3.1 不同玉米種質(zhì)響應高溫、干旱脅迫的關(guān)鍵基因

    生產(chǎn)上高溫往往伴隨干旱,不同基因型的玉米種質(zhì)存在耐旱、耐熱性差異[4,7]。本研究轉(zhuǎn)錄組測序結(jié)果顯示,高溫、干旱脅迫下不同種質(zhì)來源的4個玉米基因型存在大量差異表達基因,4個玉米基因型均有100個相同的同時響應高溫和干旱脅迫的差異表達基因,在耐旱、耐熱性強的玉米自交系中有18個特異的差異表達基因,由此推測這18個差異表達基因調(diào)節(jié)玉米對高溫、干旱脅迫的適應性。植物對高溫、干旱脅迫的分子響應是一個包括眾多轉(zhuǎn)錄調(diào)控因子的調(diào)控網(wǎng)絡[31]。熱激轉(zhuǎn)錄因子(HSFs)在玉米響應高溫、干旱脅迫中發(fā)揮重要作用,植物中的HSFs包括A、B和C 3類,除響應高溫脅迫外,HSFs還參與其他逆境脅迫的響應[32]。本研究發(fā)現(xiàn),高溫、干旱脅迫下,一個B類HSFs在耐旱、耐熱性強的玉米自交系中顯著上調(diào)。GT轉(zhuǎn)錄因子、WRKY轉(zhuǎn)錄因子和鋅指轉(zhuǎn)錄因子也是玉米響應非生物脅迫的主要調(diào)控因子[33]。這3個轉(zhuǎn)錄因子同樣存在于本研究18個特異的差異表達基因編碼中。轉(zhuǎn)錄組分析結(jié)果表明,除轉(zhuǎn)錄調(diào)控因子外,還有一些基因參與玉米響應高溫、干旱脅迫,這些基因編碼重要的結(jié)構(gòu)蛋白和功能蛋白,其中Zm00001d035139和Zm00001d035475基因分別編碼MA3結(jié)構(gòu)域蛋白和異戊烯輔酶A脫氫酶。

    3.2 不同玉米種質(zhì)響應高溫、干旱脅迫的通路

    通路分析是脅迫生物學中常用的方法。通路即代謝途徑,通過調(diào)控植物體內(nèi)物質(zhì)合成、生理活性和光合產(chǎn)物轉(zhuǎn)運,進而響應逆境脅迫[34]。生物學過程、代謝過程和生物學調(diào)節(jié)是植物體內(nèi)重要的代謝通路,參與能量物質(zhì)和生理活性物質(zhì)的合成與運輸[35]。Zenda等[36]研究發(fā)現(xiàn),耐旱玉米自交系在干旱脅迫下生物學過程、代謝過程和生物學調(diào)節(jié)通路上的差異表達基因高度富集,這與本研究KEGG通路分析結(jié)果一致。氮代謝途徑是植物生長發(fā)育過程中的一條主要的生理代謝通路,直接參與植物細胞成分合成與細胞活性調(diào)節(jié)[33]。氮代謝途徑是旱生植物響應干旱脅迫的關(guān)鍵通路,遍在蛋白作為調(diào)控因子參與植物逆境脅迫響應,在高溫、干旱等逆境脅迫下遍在蛋白代謝通路上的差異表達基因高度富集[37]。本研究KEGG通路分析結(jié)果發(fā)現(xiàn),耐旱、耐熱性強的玉米基因型在高溫、干旱脅迫下,氮代謝和遍在蛋白代謝通路上的差異表達基因顯著富集。

    綜上,在高溫、干旱脅迫下,不同玉米種質(zhì)耐旱、耐熱基因通過差異表達調(diào)控其代謝通路,從而顯現(xiàn)生理上的耐旱、耐熱性差異,使玉米耐旱、耐熱基因的篩選成為可能?;谇捌谔镩g鑒定及本研究結(jié)果,熱帶、亞熱帶玉米種質(zhì)耐旱、耐熱性強于溫帶玉米種質(zhì),在熱帶、亞熱帶玉米種質(zhì)中篩選耐旱、耐熱基因更有效。

    參考文獻:

    [1] KADAM N N, XIAO G, MELGAR R J, et al. Agronomic and physiological responses to high temperature, drought, and elevated CO2 interactions in cereals[J]. Advances in Agronomy, 2014, 127: 111-156.

    [2] LOBELL D B, SCHLENKER W, COSTA-ROBERTS J. Climate trends and global crop production since 1980[J]. Science, 2011, 333: 616-620.

    [3] GONG F, YANG L, TAI F, et al. “Omics” of maize stress response for sustainable food production: opportunities and challenges[J]. OMICS, 2014, 18: 714-732.

    [4] 王曉琴,袁繼超,熊慶娥. 玉米抗旱性的研究及展望[J]. 玉米科學,2002,10(1):57-60.

    [5] 齊 健,宋鳳斌,劉勝群. 苗期玉米根葉對干旱脅迫的生理響應[J]. 生態(tài)環(huán)境,2006,15(6):1264-1268.

    [6] 吳 昊,李燕敏,謝傳曉. 作物耐熱生理基礎(chǔ)與基因發(fā)掘研究進展[J]. 作物雜志,2018(5):1-9.

    [7] 劉克祿,陳衛(wèi)國. 植物耐熱相關(guān)基因研究進展[J]. 植物遺傳資源學報,2015,16(1):127-132.

    [8] 付 景,孫寧寧,劉天學,等. 高溫脅迫對玉米形態(tài)、葉片結(jié)構(gòu)及其產(chǎn)量的影響[J].玉米科學,2019,27(1):46-53.

    [9] KARIM M A, FRACHEBOUD Y, STAMP P. Effect of high temperature on seedling growth and photosynthesis of tropical maize genotypes[J]. Agronomy and Crop Sci, 2000, 184: 217-223.

    [10]HASANUZZAMAN M, NAHAR K, ALAM M M, et al. Physiological, biochemical, and molecular mechanisms of heat stress tolerance in plants[J]. Int J Mol Sci, 2013, 14: 9643-9684.

    [11]BIAT C A, GERATS T. Plant tolerance to high temperature in a changing environment: scientific fundamentals and production of heat stress-tolerant crops[J]. Front Plant Sci, 2013, 4: 273.

    [12]KARIM M A, FRACHEBOUD Y, STAMP P. Photosynthetic activity of developing leaves of Zea mays less affected by heat stress than that of developed leaves[J]. Plant Physiol, 1999, 105: 685-691.

    [13]BERRY J, BJOKMAN O. Photosynthetic response and adaptation to temperature in higher plants[J]. Annu Rev Plant Physiol, 1980, 98(673): 491-543.

    [14]WATERS A J, MAKAREVITCH I, NOSHAY J. Natural variation for gene expression responses to abiotic stress in maize[J]. The Plant Journal, 2017, 89: 706-717.

    [15]裴志超,張偉強,周繼華,等. 干旱脅迫對不同基因型玉米產(chǎn)量及其構(gòu)成因素的影響[J]. 玉米科學,2019,27(4):115-124.

    [16]MCGETTIGAN P A. Transcriptomics in the RNA-seq era[J]. Curr Opin Chem Biol, 2013, 17: 4-11.

    [17]CASARETTO J A, El-KEREAMY A, ZENG B, et al. Expression of OsMYB55 in maize activates stress-responsive genes and enhances heat and drought tolerance[J]. BMC Genomics, 2016, 17: 312-327.

    [18]MEYER R S, PURUGGANAN M D. Evolution of crop species: genetics of domestication and diversification[J]. Nat Rev Genet, 2013, 14: 840-852.

    [19]張 健,唐 露,張雅潔,等. 轉(zhuǎn)錄組測序技術(shù)在植物水淹脅迫研究中的應用[J]. 分子植物育種,2019,17(4):1191-1202.

    [20]LI P C, CAO W, FANG H M, et al. Transcriptic profiling maize (Zay mays L.) leaf response to abiotic stresses at the seedling stage[J]. Front Plant Sci, 2017, 8: 290-305.

    [21]張麗梅,賀立源,李建生. 玉米自交系耐低磷材料苗期篩選研究 [J]. 中國農(nóng)業(yè)科學,2004,37(12): 1955-1959.

    [22]HORTENSTEINER S. Chlorophyll degradation during senescence[J]. Annu Rev Plant Biol, 2006, 57:55-77.

    [23]CHEN J P, BURKE J J, XIN Z G. Chlorophyll fluorescence analysis revealed essential roles of FtsH11 protease in regulation of the adaptive responses of photosynthetic systems to high temperature[J]. BMC Plant Biology, 2018, 18: 11-24.

    [24]DABROWSKI P, BACZEWSKA A H, PAWLUSKIEWICZ B, et al. Prompt chlorophyll a fluorescence as a rapid tool for diagnostic changes in PSII structure inhibited by salt stress in perennial ryegrass[J]. J Photochem Photobiol B, 2016, 157: 22-31.

    [25]SU L Y, DAI Z W, LI S H. A novel system for evaluating drought-cold tolerance of grapevines using chlorophyll fluorescence[J]. BMC Plant Biology, 2015, 15: 82-94.

    [26]SCHREIBER U, SCHLIWA U, BILGER W. Continuous recording of photochemical and non-photochemical chlorophyll fluorescence quenching with a new type of modulation fluorometer[J]. Photosynth Res, 1986, 10: 51-62.

    [27]YU Q, SHEN Y, WANG Q, et al. Light deficiency and waterlogging affect chlorophyll metabolism and photosynthesis in Magnolia sinostellata[J]. Trees, 2019, 33: 11-22.

    [28]李小芳,張志良. 植物生理學實驗指導 [M]. 北京:高等教育出版社,2015.

    [29]KIM D, LANGMEAD B, SALZBERG S L. HISAT: A fast spliced aligner with low memory requirements[J]. Nature Methods, 2015, 12(4): 357-360.

    [30]WU T D, NACU S. Fast and SNP-tolerant detection of complex variants and splicing in short reads[J]. Bioinformatics, 2010, 26: 873-881.

    [31]BIANCHI V J, RUBIO M, TRAINOTTI L, et al. Prunus transcription factors: breeding perspectives[J]. Front Plant Sci, 2015, 6: 443-449.

    [32]ZENDA T S, LIU S T, WANG X, et al. Comparative proteomic and physiological analyses of two divergent maize inbred lines provide more insights into drought-stress tolerance mechanisms[J]. Int J Mol Sci, 2018, 19: 3225.

    [33]WANG H, WANG H, SHAO H, et al. Recent advances in utilizing transcription factors to improve plant abiotic stress tolerance by transgenic technology[J]. Front Plant Sci, 2016,7:67-73.

    [34]LUO M, ZHAO Y, WANG Y, et al. Comparative proteomics of contrasting maize genotypes provides insights into salt-stress tolerance mechanisms[J]. J. Proteome Res, 2018, 17: 141-153.

    [35]SHINDE H, TANAKA K, DUDHATE A, et al. Comparative de novo transcriptomic profiling of the salinity stress responsiveness in contrasting pearl millet lines[J]. Environ Exp Bot, 2018, 155: 619-627.

    [36]ZENDA T S, LIU S T, WANG X, et al. Key maize drought-responsive genes and pathways revealed by comparative transcriptome and physiological analyses of contrasting inbred lines[J]. Int J Mol Sci, 2019,20:1268.

    [37]AHMADI A, EMAM Y, PESSARAKLI M. Biochemical changes in maize seedlings exposed to drought stress conditions at different nitrogen levels[J]. J Plant Nutr, 2010, 33: 541-556.

    (責任編輯:陳海霞)

    收稿日期:2020-03-17

    基金項目:國家自然科學基金項目(31371633);重慶市教委優(yōu)秀成果轉(zhuǎn)化項目(KJZH17133)

    作者簡介:姚啟倫(1964-),男,重慶萬州人,博士,教授,主要從事玉米遺傳育種研究。(Tel)13908253365;(E-mail)yql641@aliyun.com

    猜你喜歡
    差異表達基因轉(zhuǎn)錄組干旱脅迫
    生物信息學分析患有乳腺癌的乳腺球樣本中與自我更新相關(guān)的關(guān)鍵基因
    基于轉(zhuǎn)錄組測序的山茱萸次生代謝生物合成相關(guān)基因的挖掘
    金釵石斛轉(zhuǎn)錄組SSR位點信息分析
    一氧化氮參與水楊酸對玉米幼苗根系抗旱性的調(diào)控
    人參屬藥用植物轉(zhuǎn)錄組研究進展
    干旱脅迫對扁豆生長與生理特性的影響
    不同水分條件下硫肥對玉米幼苗葉片光合特性的影響
    干旱脅迫對金花茶幼苗光合生理特性的影響
    大豆轉(zhuǎn)錄組測序研究進展綜述
    條斑紫菜優(yōu)良品系的基因芯片表達譜分析
    性少妇av在线| 国产亚洲精品一区二区www | 99re6热这里在线精品视频| 亚洲性夜色夜夜综合| 国产精品久久久久久人妻精品电影 | av有码第一页| 免费高清在线观看视频在线观看| 欧美中文综合在线视频| av不卡在线播放| 十八禁高潮呻吟视频| 欧美日本中文国产一区发布| 亚洲人成电影免费在线| 中文字幕最新亚洲高清| 一区二区三区激情视频| 黄片小视频在线播放| 国产精品久久久久久人妻精品电影 | 亚洲精品国产一区二区精华液| 国产1区2区3区精品| 性色av一级| 久久毛片免费看一区二区三区| 永久免费av网站大全| 女人精品久久久久毛片| 超碰成人久久| 欧美在线黄色| 成年美女黄网站色视频大全免费| 成人国产av品久久久| 日韩 亚洲 欧美在线| 国产精品麻豆人妻色哟哟久久| 午夜福利免费观看在线| 日韩一卡2卡3卡4卡2021年| 啦啦啦中文免费视频观看日本| 狂野欧美激情性bbbbbb| 国产av一区二区精品久久| 大型av网站在线播放| 欧美黄色淫秽网站| 人妻人人澡人人爽人人| 久久国产精品男人的天堂亚洲| 俄罗斯特黄特色一大片| 欧美一级毛片孕妇| av不卡在线播放| 国产在线观看jvid| 亚洲七黄色美女视频| 欧美黄色淫秽网站| a 毛片基地| 美女午夜性视频免费| 免费女性裸体啪啪无遮挡网站| 国产区一区二久久| 亚洲一区中文字幕在线| 久久天堂一区二区三区四区| 久久免费观看电影| 成人18禁高潮啪啪吃奶动态图| 国产精品1区2区在线观看. | 9色porny在线观看| 老司机靠b影院| 国产无遮挡羞羞视频在线观看| 韩国精品一区二区三区| 人人妻,人人澡人人爽秒播| 人人妻,人人澡人人爽秒播| 悠悠久久av| 成人国产一区最新在线观看| 一级黄色大片毛片| 婷婷色av中文字幕| 飞空精品影院首页| 亚洲少妇的诱惑av| 国产一区有黄有色的免费视频| 国产一区有黄有色的免费视频| 亚洲视频免费观看视频| 日本91视频免费播放| 美女午夜性视频免费| 亚洲精品中文字幕在线视频| 一级a爱视频在线免费观看| 国产精品久久久久久精品古装| 丰满饥渴人妻一区二区三| 精品福利观看| 国产精品国产三级国产专区5o| 午夜视频精品福利| 国产成人系列免费观看| 国产av一区二区精品久久| av有码第一页| 欧美精品一区二区大全| 欧美精品一区二区免费开放| 国产一区二区在线观看av| 亚洲精品国产一区二区精华液| 天堂8中文在线网| 99九九在线精品视频| 午夜激情av网站| 黑人巨大精品欧美一区二区蜜桃| 亚洲成人国产一区在线观看| 久久精品亚洲av国产电影网| 热re99久久国产66热| 精品亚洲乱码少妇综合久久| 亚洲黑人精品在线| 精品一区二区三区av网在线观看 | 亚洲 国产 在线| 成在线人永久免费视频| 国产亚洲午夜精品一区二区久久| 久热爱精品视频在线9| 久久毛片免费看一区二区三区| 黄片大片在线免费观看| 黑人欧美特级aaaaaa片| 亚洲国产欧美日韩在线播放| 亚洲精品粉嫩美女一区| 777米奇影视久久| 成年人午夜在线观看视频| 久久精品国产亚洲av香蕉五月 | 亚洲国产精品一区二区三区在线| 一级毛片电影观看| 国产高清视频在线播放一区 | 女性被躁到高潮视频| 日韩精品免费视频一区二区三区| 18在线观看网站| 人人澡人人妻人| 欧美 日韩 精品 国产| 一个人免费在线观看的高清视频 | 丝袜在线中文字幕| 精品国产乱码久久久久久男人| 婷婷丁香在线五月| 精品国产超薄肉色丝袜足j| 欧美日本中文国产一区发布| 精品欧美一区二区三区在线| 蜜桃国产av成人99| 精品国产乱码久久久久久男人| 午夜激情久久久久久久| 亚洲精品国产区一区二| 大陆偷拍与自拍| 最新在线观看一区二区三区| 欧美成人午夜精品| 国产精品自产拍在线观看55亚洲 | 性少妇av在线| 99热全是精品| 最近中文字幕2019免费版| 久久久久视频综合| 热re99久久精品国产66热6| 亚洲精品国产精品久久久不卡| 亚洲伊人久久精品综合| 午夜激情av网站| 人人妻人人爽人人添夜夜欢视频| 女人精品久久久久毛片| 建设人人有责人人尽责人人享有的| 国产欧美亚洲国产| 亚洲欧洲日产国产| 日韩欧美国产一区二区入口| 久久久久国内视频| 在线观看舔阴道视频| 欧美97在线视频| 日韩大码丰满熟妇| 91九色精品人成在线观看| 美女午夜性视频免费| a在线观看视频网站| av国产精品久久久久影院| 亚洲少妇的诱惑av| 国产精品成人在线| 99久久综合免费| 精品亚洲乱码少妇综合久久| 国产亚洲欧美在线一区二区| 日韩欧美免费精品| 亚洲黑人精品在线| 欧美日韩视频精品一区| 精品久久久久久电影网| 亚洲中文日韩欧美视频| 两个人免费观看高清视频| 亚洲色图 男人天堂 中文字幕| 亚洲av片天天在线观看| av欧美777| 国产精品偷伦视频观看了| 免费一级毛片在线播放高清视频 | 日韩视频一区二区在线观看| 国产区一区二久久| 伦理电影免费视频| 亚洲欧美一区二区三区黑人| 亚洲欧美一区二区三区黑人| 久久免费观看电影| 人人妻人人澡人人看| 性色av乱码一区二区三区2| 天天添夜夜摸| tube8黄色片| 久久国产精品男人的天堂亚洲| 成年人免费黄色播放视频| 香蕉丝袜av| 久久久久久久国产电影| 91九色精品人成在线观看| 女人被躁到高潮嗷嗷叫费观| 女人高潮潮喷娇喘18禁视频| 国产亚洲欧美在线一区二区| 久久香蕉激情| 精品少妇久久久久久888优播| 另类精品久久| 亚洲,欧美精品.| 欧美日韩成人在线一区二区| 精品第一国产精品| 久久久久久免费高清国产稀缺| a级片在线免费高清观看视频| 黄片大片在线免费观看| 久久久久久久久免费视频了| 国产亚洲一区二区精品| 中文字幕另类日韩欧美亚洲嫩草| 啦啦啦在线免费观看视频4| 精品亚洲成a人片在线观看| 久久久久久久久久久久大奶| 久久精品国产综合久久久| 久久亚洲精品不卡| bbb黄色大片| 不卡av一区二区三区| 黄片小视频在线播放| 国产区一区二久久| 在线亚洲精品国产二区图片欧美| 久久天躁狠狠躁夜夜2o2o| 一区二区三区四区激情视频| 久久精品国产亚洲av香蕉五月 | www.999成人在线观看| 在线av久久热| 久久久久久人人人人人| 18在线观看网站| 天天躁日日躁夜夜躁夜夜| 欧美午夜高清在线| 成年女人毛片免费观看观看9 | 我的亚洲天堂| 久热这里只有精品99| 久9热在线精品视频| 一本综合久久免费| 亚洲精品粉嫩美女一区| 中文精品一卡2卡3卡4更新| 一区二区三区激情视频| 伊人亚洲综合成人网| av国产精品久久久久影院| 美女中出高潮动态图| 后天国语完整版免费观看| 欧美性长视频在线观看| 亚洲,欧美精品.| 嫁个100分男人电影在线观看| 成人影院久久| 国产精品熟女久久久久浪| 蜜桃国产av成人99| 精品亚洲乱码少妇综合久久| 精品亚洲成国产av| 国产在线免费精品| 亚洲av电影在线进入| 亚洲av国产av综合av卡| 美女扒开内裤让男人捅视频| 中文字幕av电影在线播放| 12—13女人毛片做爰片一| 蜜桃在线观看..| 国产男女超爽视频在线观看| 欧美精品人与动牲交sv欧美| 人妻久久中文字幕网| 亚洲午夜精品一区,二区,三区| 日韩精品免费视频一区二区三区| 伊人亚洲综合成人网| av天堂久久9| 亚洲精品国产区一区二| 日韩视频一区二区在线观看| 亚洲avbb在线观看| 亚洲五月色婷婷综合| 亚洲成av片中文字幕在线观看| 别揉我奶头~嗯~啊~动态视频 | 99精品欧美一区二区三区四区| tocl精华| 黄色毛片三级朝国网站| 91字幕亚洲| 久久国产精品大桥未久av| 各种免费的搞黄视频| 欧美精品亚洲一区二区| 亚洲av美国av| 如日韩欧美国产精品一区二区三区| 日韩免费高清中文字幕av| 在线永久观看黄色视频| 2018国产大陆天天弄谢| 亚洲欧美色中文字幕在线| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品在线美女| 亚洲人成电影观看| 国产精品一区二区精品视频观看| 亚洲欧美精品自产自拍| 国产亚洲精品一区二区www | www.999成人在线观看| 日本猛色少妇xxxxx猛交久久| 人人妻,人人澡人人爽秒播| 香蕉国产在线看| 欧美精品亚洲一区二区| 人成视频在线观看免费观看| 五月天丁香电影| 女人久久www免费人成看片| 国产日韩一区二区三区精品不卡| 另类亚洲欧美激情| 亚洲欧洲日产国产| 一区福利在线观看| 亚洲国产欧美在线一区| 波多野结衣av一区二区av| 午夜精品久久久久久毛片777| 久久久久国产精品人妻一区二区| cao死你这个sao货| 国产精品久久久久成人av| 日韩熟女老妇一区二区性免费视频| 制服人妻中文乱码| 欧美+亚洲+日韩+国产| 天堂中文最新版在线下载| 中文字幕av电影在线播放| 日本五十路高清| 满18在线观看网站| 老司机深夜福利视频在线观看 | 亚洲一区中文字幕在线| 国产福利在线免费观看视频| 美国免费a级毛片| 成人免费观看视频高清| 亚洲精品第二区| 大香蕉久久网| 黄色毛片三级朝国网站| 国产精品二区激情视频| 日本vs欧美在线观看视频| 九色亚洲精品在线播放| 国产精品亚洲av一区麻豆| 精品第一国产精品| 亚洲一区中文字幕在线| 狠狠精品人妻久久久久久综合| 国产主播在线观看一区二区| 亚洲色图 男人天堂 中文字幕| 中文字幕精品免费在线观看视频| 精品亚洲乱码少妇综合久久| 69av精品久久久久久 | 免费女性裸体啪啪无遮挡网站| 夫妻午夜视频| 国产精品久久久久久精品古装| 国产一区二区三区在线臀色熟女 | 亚洲色图综合在线观看| 亚洲全国av大片| 精品免费久久久久久久清纯 | 国产免费现黄频在线看| 午夜激情久久久久久久| 午夜两性在线视频| 人人妻人人澡人人爽人人夜夜| 亚洲成国产人片在线观看| 一二三四在线观看免费中文在| 免费观看av网站的网址| 免费在线观看视频国产中文字幕亚洲 | 亚洲熟女精品中文字幕| 亚洲成人免费电影在线观看| 桃红色精品国产亚洲av| 亚洲三区欧美一区| 久久九九热精品免费| 9热在线视频观看99| 国产有黄有色有爽视频| 一本久久精品| 中文字幕人妻丝袜制服| 久久精品亚洲熟妇少妇任你| 69av精品久久久久久 | 午夜影院在线不卡| 老司机福利观看| tube8黄色片| 最新在线观看一区二区三区| 欧美av亚洲av综合av国产av| 18禁裸乳无遮挡动漫免费视频| 新久久久久国产一级毛片| 欧美日韩一级在线毛片| 每晚都被弄得嗷嗷叫到高潮| 久久香蕉激情| 免费在线观看影片大全网站| 国产av又大| 精品国产超薄肉色丝袜足j| 免费高清在线观看视频在线观看| 亚洲熟女毛片儿| 欧美日韩成人在线一区二区| 国产无遮挡羞羞视频在线观看| 亚洲av国产av综合av卡| 热re99久久精品国产66热6| 国产精品偷伦视频观看了| 黄网站色视频无遮挡免费观看| 久久免费观看电影| 亚洲精品第二区| 久久精品成人免费网站| 建设人人有责人人尽责人人享有的| 欧美亚洲 丝袜 人妻 在线| 97精品久久久久久久久久精品| 母亲3免费完整高清在线观看| 在线十欧美十亚洲十日本专区| 精品免费久久久久久久清纯 | 亚洲av片天天在线观看| 俄罗斯特黄特色一大片| 色婷婷久久久亚洲欧美| 最近最新免费中文字幕在线| 侵犯人妻中文字幕一二三四区| 久久人妻熟女aⅴ| 亚洲va日本ⅴa欧美va伊人久久 | 国产免费一区二区三区四区乱码| 人人妻人人澡人人爽人人夜夜| 99久久综合免费| 成在线人永久免费视频| 在线十欧美十亚洲十日本专区| 久久人妻熟女aⅴ| 啦啦啦视频在线资源免费观看| 19禁男女啪啪无遮挡网站| 久久人人爽av亚洲精品天堂| 日本黄色日本黄色录像| 少妇猛男粗大的猛烈进出视频| 国产免费现黄频在线看| 夜夜骑夜夜射夜夜干| 免费观看人在逋| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久亚洲精品国产蜜桃av| 无限看片的www在线观看| 亚洲视频免费观看视频| 99热网站在线观看| 动漫黄色视频在线观看| 最近最新中文字幕大全免费视频| 少妇猛男粗大的猛烈进出视频| 欧美午夜高清在线| 可以免费在线观看a视频的电影网站| 纯流量卡能插随身wifi吗| 亚洲精品日韩在线中文字幕| 高清欧美精品videossex| 啦啦啦 在线观看视频| 亚洲精品国产色婷婷电影| 男女午夜视频在线观看| 国产欧美日韩一区二区精品| 国产精品免费视频内射| 精品国产一区二区三区久久久樱花| 国产一区二区 视频在线| 午夜影院在线不卡| 99久久国产精品久久久| 亚洲精品在线美女| 亚洲精品自拍成人| 建设人人有责人人尽责人人享有的| 亚洲色图 男人天堂 中文字幕| 国产欧美日韩一区二区三区在线| 亚洲欧美成人综合另类久久久| 国产亚洲精品第一综合不卡| 久久久久国产一级毛片高清牌| 久久中文看片网| 久久影院123| 99re6热这里在线精品视频| 精品国内亚洲2022精品成人 | 嫩草影视91久久| 国产成人免费观看mmmm| 男女无遮挡免费网站观看| av一本久久久久| 黄片小视频在线播放| 亚洲 欧美一区二区三区| 国产精品一区二区在线观看99| 精品久久久精品久久久| 69av精品久久久久久 | 午夜福利,免费看| 丁香六月欧美| 丝袜美腿诱惑在线| 国产精品一区二区精品视频观看| tocl精华| 国产日韩欧美亚洲二区| √禁漫天堂资源中文www| 久久精品aⅴ一区二区三区四区| 久久国产精品人妻蜜桃| av有码第一页| 老司机亚洲免费影院| 久久天躁狠狠躁夜夜2o2o| 超碰97精品在线观看| 最新在线观看一区二区三区| 国产精品香港三级国产av潘金莲| 人人澡人人妻人| 欧美日韩中文字幕国产精品一区二区三区 | 色婷婷av一区二区三区视频| 午夜福利视频在线观看免费| 女人精品久久久久毛片| 考比视频在线观看| 日本av手机在线免费观看| 午夜福利在线免费观看网站| av电影中文网址| 纵有疾风起免费观看全集完整版| 另类精品久久| 亚洲av男天堂| 在线永久观看黄色视频| 免费一级毛片在线播放高清视频 | 久久久国产一区二区| 宅男免费午夜| 亚洲国产欧美网| 男女下面插进去视频免费观看| 午夜两性在线视频| 亚洲国产av新网站| 久久亚洲精品不卡| 啦啦啦中文免费视频观看日本| 国产精品av久久久久免费| 国产精品久久久久久精品电影小说| 免费人妻精品一区二区三区视频| 91字幕亚洲| 中文精品一卡2卡3卡4更新| a级片在线免费高清观看视频| 91精品国产国语对白视频| 亚洲欧美精品自产自拍| 精品国产超薄肉色丝袜足j| 精品一区二区三区四区五区乱码| 国产主播在线观看一区二区| 啦啦啦中文免费视频观看日本| 久久中文字幕一级| 侵犯人妻中文字幕一二三四区| 日韩有码中文字幕| 国产xxxxx性猛交| 国产一区二区 视频在线| 90打野战视频偷拍视频| 九色亚洲精品在线播放| 国产欧美日韩综合在线一区二区| 精品一区二区三卡| 我的亚洲天堂| 黄频高清免费视频| 热re99久久精品国产66热6| 99久久精品国产亚洲精品| 国产精品香港三级国产av潘金莲| 精品少妇久久久久久888优播| 香蕉丝袜av| 欧美成人午夜精品| 性少妇av在线| 久久久欧美国产精品| 亚洲中文日韩欧美视频| 国产精品 国内视频| 成人av一区二区三区在线看 | 深夜精品福利| 中文字幕精品免费在线观看视频| 在线观看免费日韩欧美大片| 午夜免费成人在线视频| svipshipincom国产片| 日韩三级视频一区二区三区| 超碰成人久久| 国产一区二区三区av在线| av在线app专区| 成人影院久久| 最近中文字幕2019免费版| 久久人人爽人人片av| bbb黄色大片| 午夜免费观看性视频| 黄色毛片三级朝国网站| 老鸭窝网址在线观看| 女人被躁到高潮嗷嗷叫费观| 一级片'在线观看视频| 国产在线视频一区二区| 久久久久网色| 99热国产这里只有精品6| 欧美97在线视频| 国产99久久九九免费精品| 亚洲欧美精品自产自拍| 丝袜在线中文字幕| 欧美日韩亚洲国产一区二区在线观看 | netflix在线观看网站| a级毛片在线看网站| 韩国精品一区二区三区| 伦理电影免费视频| 黄色视频,在线免费观看| 日韩 亚洲 欧美在线| 国产精品秋霞免费鲁丝片| 80岁老熟妇乱子伦牲交| 9191精品国产免费久久| 欧美日韩成人在线一区二区| 无遮挡黄片免费观看| 久久人人爽av亚洲精品天堂| 亚洲欧美精品自产自拍| 女人被躁到高潮嗷嗷叫费观| www.精华液| 成人免费观看视频高清| 少妇 在线观看| 国产精品秋霞免费鲁丝片| 他把我摸到了高潮在线观看 | 无遮挡黄片免费观看| 人人澡人人妻人| 日本a在线网址| 免费在线观看影片大全网站| 欧美另类亚洲清纯唯美| 国产免费现黄频在线看| 久久国产精品影院| e午夜精品久久久久久久| 午夜福利,免费看| 欧美性长视频在线观看| 精品熟女少妇八av免费久了| 国产主播在线观看一区二区| 99九九在线精品视频| 久久精品亚洲av国产电影网| 久久国产精品男人的天堂亚洲| 日日爽夜夜爽网站| 一级片免费观看大全| av天堂在线播放| 人人澡人人妻人| 国产成人精品无人区| 韩国精品一区二区三区| 亚洲av电影在线观看一区二区三区| 精品国产一区二区三区久久久樱花| 精品国产超薄肉色丝袜足j| 国产男人的电影天堂91| 日韩熟女老妇一区二区性免费视频| 丰满少妇做爰视频| 91成年电影在线观看| 99久久综合免费| 天堂中文最新版在线下载| 五月天丁香电影| 又大又爽又粗| 成年动漫av网址| 女人精品久久久久毛片| xxxhd国产人妻xxx| 18禁国产床啪视频网站| 91精品国产国语对白视频| 男女下面插进去视频免费观看| 最近最新中文字幕大全免费视频| 黑人巨大精品欧美一区二区mp4| 国产精品成人在线| 男人添女人高潮全过程视频| 免费在线观看日本一区| 亚洲av日韩在线播放| 纯流量卡能插随身wifi吗| 多毛熟女@视频| 老司机影院毛片| 亚洲五月色婷婷综合| 操美女的视频在线观看| 国产主播在线观看一区二区| 中文字幕高清在线视频| 欧美精品一区二区大全| 人人妻,人人澡人人爽秒播| 丰满迷人的少妇在线观看| 男女之事视频高清在线观看| 亚洲国产欧美日韩在线播放| 精品欧美一区二区三区在线| 一二三四社区在线视频社区8| 久久国产亚洲av麻豆专区| a在线观看视频网站|