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

    萱草根莖低溫脅迫轉(zhuǎn)錄組分析

    2021-08-26 05:51:04黃長(zhǎng)兵程培蕾楊紹宗張煥朝姜正之金立敏
    關(guān)鍵詞:萱草差異基因低溫

    黃長(zhǎng)兵,程培蕾,楊紹宗,張煥朝,姜正之,金立敏

    (1.南京林業(yè)大學(xué) 南方現(xiàn)代林業(yè)協(xié)同創(chuàng)新中心/林學(xué)院,江蘇 南京 210037; 2.蘇州農(nóng)業(yè)職業(yè)技術(shù)學(xué)院 園藝科技學(xué)院,江蘇 蘇州 215008; 3.蘇州市華冠園創(chuàng)園藝科技有限公司,江蘇 蘇州 215557)

    萱草(Hemerocallisfulva)又被稱為蘐草、忘憂等,屬于百合科萱草屬多年生植物,其根莖較短但其紡錘形肉質(zhì)根較厚[1-3],在中國(guó)有著悠久的栽培歷史。萱草既可以作為園林植物作觀賞之用,也可以作為蔬菜食用,亦可以作為藥物治??;由于其品種多樣,色彩豐富,抗逆性強(qiáng)等因素,成為道路綠化、公園綠地、屋頂花園等生態(tài)城市建設(shè)廣泛使用的園林綠化材料[4]。作為萱草屬中較少幾種可食用植物,萱草擁有較高的經(jīng)濟(jì)價(jià)值。此外,萱草中富含大量次級(jí)代謝產(chǎn)物,這些代謝物對(duì)于人體的焦慮、腫脹等癥狀具有良好的治療效果[5-6]。目前,研究者們已經(jīng)對(duì)萱草的生理生化、分類鑒定、組織培養(yǎng)和育種等領(lǐng)域進(jìn)行了較多研究[7-11],但有關(guān)于萱草逆境脅迫的報(bào)道和研究相對(duì)較少。低溫脅迫作為最常見(jiàn)的非生物脅迫因素之一,能夠?qū)е轮参锏纳L(zhǎng)和發(fā)育出現(xiàn)延遲,同時(shí)也會(huì)影響作物的產(chǎn)量[12-14]。冷害是影響萱草觀賞價(jià)值和經(jīng)濟(jì)價(jià)值的重要因素,因此,研究萱草在低溫脅迫下基因的表達(dá)變化對(duì)于解決萱草的冷害具有重要意義。

    轉(zhuǎn)錄組測(cè)序與生物信息學(xué)分析技術(shù)相結(jié)合可以大大促進(jìn)功能基因組學(xué)的研究,通過(guò)測(cè)序得到一定條件下細(xì)胞的整體轉(zhuǎn)錄信息,為從轉(zhuǎn)錄水平來(lái)發(fā)掘生物體的生長(zhǎng)發(fā)育、次級(jí)代謝與調(diào)控機(jī)制[15-16]提供了便利。近些年,研究者們用轉(zhuǎn)錄組分析策略來(lái)研究植物抗低溫脅迫基因組學(xué),已經(jīng)獲得擬南芥[17]、番茄[18]、紫花苜蓿[19]、黃瓜[20]和水稻[21]等眾多植物轉(zhuǎn)錄組數(shù)據(jù),這些都為闡明植物抗低溫脅迫機(jī)制研究奠定了基礎(chǔ)。本研究利用高通量測(cè)序平臺(tái)Illumina HiSeq對(duì)低溫脅迫條件下2種萱草的根莖進(jìn)行轉(zhuǎn)錄組測(cè)序分析,期望揭示萱草根莖轉(zhuǎn)錄組在低溫脅迫條件下發(fā)生的表達(dá)差異性特征,為萱草冷應(yīng)激功能基因的鑒定、次生代謝途徑的解析及其低溫脅迫應(yīng)激調(diào)控機(jī)制的研究提供參考。

    1 材料與方法

    1.1 材料

    萱草Jinyan 7(耐冷型品種,冬季地上部常綠)和Lucretius 4(冷敏感型品種,冬季地上部枯萎)由蘇州農(nóng)業(yè)職業(yè)技術(shù)學(xué)院相城科技園提供。在通風(fēng)溫室中(20 ℃,光周期:16 h光照8 h黑暗),植物于塑料盆中生長(zhǎng)至6片葉子,然后移入生長(zhǎng)室進(jìn)行處理。將溫度從20 ℃逐漸降低至-4 ℃開(kāi)始低溫處理,具體方法:讓植株20 ℃生長(zhǎng)48 h,然后分別在10 ℃和4 ℃誘導(dǎo)48 h,最后將它們?cè)?4 ℃處理8 h。收集在20 ℃(作為對(duì)照組,CK)和-4 ℃(作為處理組,T)生長(zhǎng)的萱草根莖,在液氮中冷凍用于RNA提取和生理分析。收集的4個(gè)樣品分別命名為4CK(20 ℃時(shí)Lucretius 4的根莖)、4T(-4 ℃處理的Lucretius 4根莖)、7CK(20 ℃時(shí)Jinyan 7根莖)、7T(-4 ℃處理的Jinyan 7根莖)。每3株植物根莖作為1個(gè)生物學(xué)重復(fù),共設(shè)置2組重復(fù)。

    1.2 RNA提取與文庫(kù)構(gòu)建

    采用TRNzol Universal Reagent植物總RNA提取試劑盒,分別提取冷處理前后耐冷型和冷敏感型萱草根莖總RNA,Agilent 2100 Bioanalyzer檢測(cè)總RNA的質(zhì)量,NanoDrop TM2000檢測(cè)總RNA的濃度。利用Oligo(dT)磁珠富集每個(gè)萱草文庫(kù)的mRNA,加入適量打斷劑高溫下片斷化mRNA,并以片段化的mRNA為模板合成cDNA,選擇合適的片段進(jìn)行PCR擴(kuò)增,構(gòu)建文庫(kù)。

    1.3 轉(zhuǎn)錄組測(cè)序與組裝

    將構(gòu)建的萱草根莖轉(zhuǎn)錄組文庫(kù)提交Illumina HiSeq平臺(tái)進(jìn)行高通量測(cè)序。首先對(duì)Illumina測(cè)序的原始reads進(jìn)行處理,去除接頭污染、未知堿基N含量偏高(>5%)和低質(zhì)量reads,然后將獲得的高質(zhì)量clean reads利用Trinity[22]進(jìn)行denovo組裝分析,最后通過(guò)Tgicl對(duì)其進(jìn)行聚類和去冗余處理,得到unigenes。

    1.4 轉(zhuǎn)錄組功能注釋

    利用Blastn對(duì)得到的unigene序列進(jìn)行Nt數(shù)據(jù)庫(kù)注釋,通過(guò)Blastx將unigenes序列與以下蛋白數(shù)據(jù)庫(kù)進(jìn)行比對(duì)(E-value<1×10-5,E-value為比對(duì)的期望值):Nr(NCBI nonredundant protein)、SwissPort、KOG(clusters of orthologous groups)和KEGG(kyoto encyclopedia of genes and genomes),得到和相應(yīng)unigenes序列相似性最高的蛋白序列,從而得到unigenes的對(duì)應(yīng)注釋信息。使用Blast2GO軟件,同時(shí)結(jié)合Nr注釋得到unigenes的GO(gene ontology)功能注釋。

    1.5 蛋白編碼框和轉(zhuǎn)錄因子預(yù)測(cè)

    通過(guò)軟件Transdecoder(https://transdecoder.github.io)來(lái)識(shí)別unigenes的候選編碼區(qū)域,具體過(guò)程主要包括通過(guò)TransDecoder.LongOrfs軟件從unigenes中將最長(zhǎng)的開(kāi)放閱讀框(open reading frame,ORF)提取出來(lái),接著將其與SwissProt數(shù)據(jù)庫(kù)進(jìn)行Blast比對(duì),其結(jié)果經(jīng)Hmmscan搜索得到Pfam蛋白的同源序列,從而實(shí)現(xiàn)蛋白編碼框(coding sequence,CDS)的預(yù)測(cè)。通過(guò)getorf軟件來(lái)檢測(cè)unigenes的ORF,接著利用Hmmsearch檢索將ORF與PlantTFDB數(shù)據(jù)庫(kù)中轉(zhuǎn)錄因子蛋白結(jié)構(gòu)域進(jìn)行比對(duì),最后根據(jù)PlantTFDB中關(guān)于轉(zhuǎn)錄因子家族特征的描述對(duì)unigenes進(jìn)行鑒定。

    1.6 Unigenes基因表達(dá)水平分析

    通過(guò)MISA(http://pgrc.ipk-gatersleben.de/misa)軟件來(lái)檢測(cè)萱草轉(zhuǎn)錄組的unigenes,搜索和驗(yàn)證簡(jiǎn)單重復(fù)序列(simple sequence repeats, SSR),并進(jìn)行相關(guān)統(tǒng)計(jì)分析。利用軟件Bowtie2將clean reads與unigenes進(jìn)行比對(duì),然后通過(guò)RSEM軟件計(jì)算得到各個(gè)樣品的基因表達(dá)水平,并參照Saddhe等[23]的差異基因檢測(cè)方法,通過(guò)軟件PossionDis [fold change≥2.00和錯(cuò)誤發(fā)現(xiàn)率(FDR)≤0.001]得到不同樣本之間的差異基因。

    1.7 差異表達(dá)基因聚類、GO功能、Pathway功能分析

    差異基因的層次聚類分析采用R軟件中的pheatmap函數(shù)分析法。依據(jù)GO功能的官方分類對(duì)差異基因進(jìn)行功能注釋,根據(jù)KEGG的官方分類對(duì)差異基因的生物通路進(jìn)行注釋分析,差異基因富集分析則利用R軟件中的phyper函數(shù)分析法。

    1.8 差異表達(dá)基因的qRT-PCR驗(yàn)證

    通過(guò)對(duì)轉(zhuǎn)錄組中差異表達(dá)基因的進(jìn)一步篩選,隨機(jī)選擇8個(gè)顯著差異表達(dá)的基因,進(jìn)行qRT-PCR驗(yàn)證。qRT-PCR檢測(cè)使用Takara公司的SYBR?Premix Ex TaqTM試劑,在熒光定量?jī)x中進(jìn)行。內(nèi)參基因?yàn)?8S rRNA,正反向引物序列分別為:5′-CCATAAACGATGCCGACC-3′;5′-TTGCGACCATACTCCCCC-3′。PCR反應(yīng)條件為:95 ℃ 5 min;95 ℃ 5 s,60 ℃ 30 s,72 ℃ 10 s,40個(gè)循環(huán)。qRT-PCR中使用的引物序列如表1所示。對(duì)每個(gè)基因和樣本進(jìn)行3個(gè)生物重復(fù)和3個(gè)技術(shù)重復(fù)。用2-ΔΔCt法計(jì)算差異基因相對(duì)表達(dá)水平。

    2 結(jié)果與分析

    2.1 萱草轉(zhuǎn)錄組組裝與質(zhì)量分析

    從表2可以看出,4個(gè)樣品4CK、4T、7CK、7T分別得到47 738 570、50 110 130、50 918 552、38 540 808條raw reads,過(guò)濾后分別產(chǎn)生46 171 660、47 791 962、47 411 290、35 880 734條高質(zhì)量clean reads,各包含6 858 685 119、7 099 131 160、7 047 333 003、5 342 671 114個(gè)核苷酸信息,Q20(質(zhì)量值≥20的堿基所占百分比)分別為98.61%、98.63%、98.72%、98.82%,Q30(質(zhì)量值≥30的堿基所占百分比)分別為95.65%、95.7%、95.97%、96.23%,GC含量分別為48.83%、48.89%、48.57%、48.48%,表明clean reads的質(zhì)量符合實(shí)驗(yàn)要求,可用于后續(xù)分析。樣品4CK、4T、7CK、7T經(jīng)Trinity組裝后各獲得55 348、51 400、5 6962、49 915個(gè)unigene,平均長(zhǎng)度分別為796、824、763、753 nt,N50(按轉(zhuǎn)錄本長(zhǎng)度從大到小排序后逐個(gè)累加至所有轉(zhuǎn)錄本總長(zhǎng)度的50%)分別為1 294、1 337、1 225、1 197 nt。

    表2 有效數(shù)據(jù)評(píng)估統(tǒng)計(jì)Table 2 Valid data evaluation statistics

    Unigenes的長(zhǎng)度分布圖(圖1)顯示,Unigenes長(zhǎng)度在300~>3 000 nt;長(zhǎng)度為300 nt的unigenes數(shù)量最多,共25 011條;400 nt的unigenes數(shù)量次之,共11 341條;其他長(zhǎng)度的unigenes均低于10 000條,大于3 000 nt的unigenes共2 194條。

    2.2 萱草轉(zhuǎn)錄組unigenes的功能注釋

    為獲得組裝好的unigenes的表達(dá)和功能注釋,將unigene序列與Nr、Nt、Swissprot、KEGG、KOG等數(shù)據(jù)庫(kù)進(jìn)行比對(duì)。結(jié)果表明,通過(guò)Nt數(shù)據(jù)庫(kù)可對(duì)10 453條unigenes(12.00%)進(jìn)行成功比對(duì)并進(jìn)行注釋,Nr數(shù)據(jù)庫(kù)可對(duì)38 447條unigenes(44.15%)進(jìn)行成功比對(duì)并進(jìn)行注釋,在SwissPort、KOG、KEGG、Pfam、GO數(shù)據(jù)庫(kù)中分別有28 095(32.27%)、41 418(47.57%)、27 835(31.97%)、36 767(42.23%)、15 933(18.30%)條unigenes獲得成功比對(duì)和注釋。

    對(duì)Nr數(shù)據(jù)庫(kù)的比對(duì)注釋結(jié)果進(jìn)行分析,圖2是unigenes注釋同源基因的物種分布圖,可以發(fā)現(xiàn),在相似序列較高的物種中,蘆筍(Asparagusofficinalis)所占比例最高,其數(shù)量可以達(dá)到24 115條(45.71%);其次為油棕(Elaeisguineensis)、海棗(Phoenixdactylifera)和鳳梨(Ananascomosus),其數(shù)量分別為5 695條(10.79%)、3 673條(6.96%)和1 638條(3.1%);注釋到其他物種的unigenes共17 633條,占比為33.43%。

    2.3 CDS和轉(zhuǎn)錄因子分析

    通過(guò)BLAST比對(duì)分析,從4組萱草樣品轉(zhuǎn)錄組的全部unigenes中共獲得42 133個(gè)CDS 序列,CDS的長(zhǎng)度全部大于200 nt,其中有16 337條CDS長(zhǎng)度超過(guò)1 000 nt,3 923條CDS長(zhǎng)度超過(guò)2 000 nt(圖3)。通過(guò)轉(zhuǎn)錄因子預(yù)測(cè)分析可以發(fā)現(xiàn),這些unigenes分屬于56個(gè)家族(圖4),而其中Homeobox、zf-C2H2、bHLH、HMG、NGFIB-like、MYB、zf-GATA、ZBTB、TF_bZIP與CG-1類占主體,這說(shuō)明眾多轉(zhuǎn)錄因子參與萱草根莖的生理代謝。

    圖3 CDS長(zhǎng)度分布圖Fig.3 Distribution diagram of coding sequence length

    圖4 轉(zhuǎn)錄因子家族分類圖Fig.4 Classification of transcription factor families

    2.4 差異表達(dá)基因統(tǒng)計(jì)、GO和KEGG分析

    表3 Unigenes注釋結(jié)果統(tǒng)計(jì)Table 3 Annotation results of unigene sequences

    2.4.1 差異表達(dá)基因統(tǒng)計(jì)

    對(duì)冷敏感型和耐冷型萱草冷處理后的差異基因數(shù)量進(jìn)行統(tǒng)計(jì)發(fā)現(xiàn),與對(duì)照相比,冷敏感型萱草根莖在低溫脅迫下,總共檢測(cè)到6 153個(gè)基因表達(dá)變化,其中,上調(diào)表達(dá)基因有1 849個(gè),下調(diào)表達(dá)基因共檢測(cè)到4 304個(gè)。而耐冷萱草根莖在低溫脅迫下,共有14 813個(gè)差異表達(dá)基因,其中,有6786個(gè)基因上調(diào)表達(dá),8 027個(gè)基因表達(dá)量下調(diào)。對(duì)照組中,與冷敏感型萱草相比,耐冷萱草有10 282個(gè)基因上調(diào)表達(dá),12 234個(gè)基因下調(diào)表達(dá);而在低溫處理?xiàng)l件下,與冷敏感型萱草相比,耐冷萱草上調(diào)表達(dá)基因有9 224個(gè),下調(diào)表達(dá)基因有10 614個(gè)(圖5)。

    圖5 上調(diào)與下調(diào)差異基因數(shù)量統(tǒng)計(jì)Fig.5 Numbers of up-regulated and down-regulated differentially expressed genes

    2.4.2 差異表達(dá)基因GO富集分析

    為了進(jìn)一步了解低溫脅迫條件下萱草的抗逆機(jī)制,對(duì)差異表達(dá)基因進(jìn)行GO富集分析。GO分析結(jié)果表明(圖6),這些差異表達(dá)基因分為3大類:分子功能(molecular function)、細(xì)胞組分(cellular component)和生物學(xué)過(guò)程(biological process),每個(gè)大類對(duì)應(yīng)的功能組數(shù)量分別為16、19、29。差異基因大部分富集在生物學(xué)過(guò)程中,而參與分子功能的差異基因相對(duì)較少。在生物學(xué)過(guò)程中,富集較明顯的包括細(xì)胞過(guò)程(7 216)、代謝過(guò)程(6 652)、生物調(diào)節(jié)(2 138)、生物過(guò)程調(diào)節(jié)(1 973)和定位(1 700);參與細(xì)胞組分的差異基因主要富集在細(xì)胞(6 705)、細(xì)胞組分(6 614)、膜(5 386)、細(xì)胞器(4 947)和膜部分(4 934);而參與分子功能的差異基因富集類別主要有催化活性(8 314)、結(jié)合(7 916)和轉(zhuǎn)運(yùn)活性(1 916)。

    圖6 差異基因GO功能分類圖Fig.6 Gene ontology classification of differentially expressed genes

    2.4.3 差異表達(dá)基因KEGG分析

    根據(jù)KEGG代謝通路對(duì)差異基因進(jìn)行功能注釋,發(fā)現(xiàn)參與的代謝通路包含6個(gè)分支:細(xì)胞過(guò)程(cellular processes)、環(huán)境信息處理(environmental information processing)、遺傳信息處理(genetic information processing)、人類疾病(human diseases)、代謝(metabolism)和有機(jī)系統(tǒng)(organismal systems),具體的差異基因數(shù)量見(jiàn)圖7。在代謝途徑中,占據(jù)前5的類別分別為全局概述圖(5 527)、碳水化合物代謝(1 694)、核苷酸代謝(1 673)、氨基酸代謝(1 311)和脂質(zhì)代謝(1 183);在遺傳與信息處理類別中,翻譯(2 612)、轉(zhuǎn)錄(2 165)和折疊分揀與降解(1 781)差異基因顯著富集,而信號(hào)轉(zhuǎn)導(dǎo)(5 416)在環(huán)境信息處理通路中差異基因富集較多。

    圖7 差異表達(dá)基因KEGG分類Fig.7 KEGG classification of differentially expressed genes

    2.4.4 萱草低溫脅迫相關(guān)代謝通路分析

    對(duì)差異基因的代謝通路進(jìn)行分析發(fā)現(xiàn),差異基因富集程度比較明顯的途徑主要包括:異黃酮生物合成(isoflavonoid biosynthesis)、?;撬岷蛠喤;撬岽x(taurine and hypotaurine metabolism)、光合天線作用蛋白(photosynthesis-antenna proteins)、二苯基庚酮和姜酚生物合成(stilbenoid, diarylheptanoid and gingerol biosynthesis)、植物晝夜節(jié)律(circadian rhythm-plant)、不飽和脂肪酸的生物合成(biosynthesis of unsaturated fatty acids)等,具體的代謝通路與差異基因數(shù)量見(jiàn)表4。

    表4 與萱草低溫脅迫作用相關(guān)的16條代謝通路Table 4 Sixteen metabolic pathways related to cold stress of Hemerocallis fulva

    2.5 差異基因參與對(duì)低溫的反應(yīng)

    通過(guò)對(duì)差異基因的進(jìn)一步篩選發(fā)現(xiàn),Ca2+信號(hào)通路相關(guān)基因表達(dá)發(fā)生明顯差異(表5)。與對(duì)照組相比,低溫處理后,冷敏感型萱草中CBL互作絲氨酸/蘇氨酸蛋白激酶3基因顯著上調(diào);

    表5 鈣離子信號(hào)通路相關(guān)差異基因Table 5 Differentially expressed genes involved in Ca2+ signalling pathway in each comparison

    在耐冷型萱草中多種Ca2+信號(hào)通路相關(guān)基因明顯差異表達(dá),鈣依賴性蛋白激酶基因CAMK2A、鈣調(diào)蛋白基因(Calm5、Calm3)表達(dá)顯著上調(diào),而鈣調(diào)蛋白和肌聯(lián)蛋白結(jié)合RhoGEF基因、CBL互作絲氨酸/蘇氨酸蛋白激酶14基因則顯著下調(diào)。對(duì)照組中,與耐冷型萱草相比,冷敏感型萱草鈣調(diào)蛋白基因(Calm5、Calm3)、CBL互作絲氨酸/蘇氨酸蛋白激酶14基因顯著上調(diào),鈣依賴性蛋白激酶基因CAMK2A顯著下調(diào)。在低溫脅迫下,與耐冷型萱草相比,冷敏感型萱草鈣調(diào)蛋白基因(Calm5、Calm3)、鈣依賴性蛋白激酶基因CAMK2A與對(duì)照組表達(dá)趨勢(shì)相同。CL7725.Contig5_All基因

    編碼了鈣調(diào)蛋白和肌聯(lián)蛋白結(jié)合RhoGEF(calmodulin and titin-interacting RhoGEF)。低溫脅迫下,與對(duì)照組相比,耐冷型萱草中CL7725.Contig5_All基因顯著上調(diào);與耐冷型萱草相比,冷敏感型萱草的CL7725.Contig5_All基因顯著下調(diào);與對(duì)照組相比,冷敏感型萱草的CL7725.Contig5_All基因表達(dá)無(wú)顯著差異。對(duì)照組中,耐冷型萱草與冷敏感型萱草的CL7725.Contig5_All基因表達(dá)也沒(méi)有顯著差異。Unigene6945_All基因編碼鈣調(diào)素結(jié)合轉(zhuǎn)錄激活子3(calmodulin-binding transcription activator 3),低溫脅迫抑制了其在耐冷萱草中的表達(dá),在對(duì)照組中,該基因在2種萱草中的表達(dá)沒(méi)有差異。

    此外,也篩選到MAPK信號(hào)通路相關(guān)差異基因(表6),低溫處理前后比對(duì)結(jié)果發(fā)現(xiàn),在冷敏感型萱草(4CK-VS-4T)中,促分裂原活化的蛋白激酶基因MAPKKKK3顯著上調(diào),MAPK1-like、MAPKKKK5顯著下調(diào)。在耐冷型萱草中(7CK-VS-7T),促分裂原活化的蛋白激酶基因MLK4-like、MAPKKKK3、MAPK1、MAPK4顯著上調(diào),MAPK1-like、MAPK7-like、MAPKKKK5顯著下調(diào)。對(duì)照組中,與冷敏感型萱草相比,MAPK8互作蛋白2基因(MAPK8IP2)、MAPKKKK3基因、MAPK7-like基因上調(diào)表達(dá),MAPK4、MAPKKKK5基因下調(diào)表達(dá)。低溫處理后,耐冷型萱草與冷敏感型萱草的MAPKKKK3、MAPK8IP2基因表達(dá)上調(diào),MAPK7-like、MLK4-like、MAPKKKK5表達(dá)下調(diào)。進(jìn)一步分析發(fā)現(xiàn),低溫可以抑制MAPK1基因(CL4304.Contig1_All)的表達(dá),常溫下,2種萱草中沒(méi)有檢測(cè)到該基因的表達(dá)差異,而在低溫下,冷敏感萱草中,其下調(diào)倍數(shù)為-1.36,在耐冷型萱草中,其表達(dá)被抑制的更加明顯,下調(diào)倍數(shù)為-3.90。

    表6 MAPK信號(hào)通路相關(guān)差異基因Table 6 Differentially expressed genes involved in MAPK cascades pathway in each comparison

    對(duì)編碼熱激蛋白(表7)相關(guān)的基因進(jìn)行篩選發(fā)現(xiàn),冷敏感型萱草(4CK-VS-4T)中,熱激蛋白基因HSP70-2、DNAJC6、HSP90顯著上調(diào),而在耐冷型萱草中(7CK-VS-7T),僅熱激蛋白基因HSP86上調(diào)表達(dá),其他熱激蛋白基因HSP90-5、HSP90均顯著下調(diào)。不同的溫度條件下,對(duì)2種萱草進(jìn)行對(duì)比發(fā)現(xiàn),對(duì)照條件下,與冷敏感型萱草相比,Hspa8基因上調(diào)表達(dá),HSP70-2、HSP90-5基因下調(diào)表達(dá),低溫處理后,2種萱草的對(duì)比中篩選到的熱激蛋白基因均下調(diào)表達(dá)。HSP90基因(Unigene6560_All)在2種萱草中表現(xiàn)出相反的表達(dá)趨勢(shì),冷敏感萱草中,低溫處理后上調(diào),而在耐冷型萱草中,與對(duì)照條件相比下調(diào)-1.61,與冷敏感萱草相比下調(diào)-2.61。

    表7 編碼熱激蛋白的差異表達(dá)基因Table 7 Differentially expressed genes encoding HSPs protein in each comparison

    對(duì)差異表達(dá)的轉(zhuǎn)錄因子(表8)進(jìn)行分析發(fā)現(xiàn),在冷敏感型萱草(4CK-VS-4T)中,轉(zhuǎn)錄因子WRKY11、MYB102顯著下調(diào)表達(dá),而在耐冷型萱草(7CK-VS-7T)中,轉(zhuǎn)錄因子MYB20、bHLH2上調(diào)表達(dá),MYB102下調(diào)表達(dá)。與冷敏感萱草相比,對(duì)照組耐冷型萱草WRKY11上調(diào)表達(dá),bHLH2、4CK-VS-4T,Lucretius 4低溫與對(duì)照相比;7CK-VS-7T,Jinyan 7低溫和對(duì)照相比;4CK-VS-7CK,對(duì)照條件下Jinyan 7與Lucretius 4相比;4T-VS-7T,低溫處理下Jinyan 7與Lucretius 4相比。下同。

    表8 差異表達(dá)的轉(zhuǎn)錄因子Table 8 Differential expression of transcription factors

    2.6 熒光定量PCR驗(yàn)證

    為了驗(yàn)證轉(zhuǎn)錄組測(cè)序結(jié)果,隨機(jī)選擇8個(gè)表達(dá)水平差異顯著的基因進(jìn)行qRT-PCR試驗(yàn),分別為:類硫氧還蛋白基因(thioredoxin-like)、絲氨酸/精氨酸重復(fù)基質(zhì)蛋白基因(serine/arginine repetitive matrix protein 3-like)、溶磷酰膽堿酰轉(zhuǎn)移酶1基因(Lysophos phatidylcholine acyltransferase 1)、1-磷脂酰肌醇 3-磷酸 5-激酶基因(1-phosphatidylinositol 3-phosphate 5-kinase)、粘液素-2基因(Mucin-2)、肌球蛋白VIIa基因(myosin-VIIa-like)、

    含有LOB結(jié)構(gòu)域的蛋白38基因(LOB domain-containing protein 38-like)、多聚泛素-A基因(polyubiquitin-a)。結(jié)果顯示:低溫脅迫下,8個(gè)基因的qRT-PCR結(jié)果與轉(zhuǎn)錄組分析的趨勢(shì)一致(圖8),只是具體表達(dá)倍數(shù)上存在差異,這些差異可能是由于轉(zhuǎn)錄組測(cè)序和qRT-PCR實(shí)驗(yàn)的檢測(cè)靈敏度,以及實(shí)驗(yàn)結(jié)果分析方法不同造成的,表明轉(zhuǎn)錄組測(cè)序結(jié)果可靠。

    3 結(jié)論與討論

    低溫是許多植物所面臨的常見(jiàn)威脅,其不僅可誘導(dǎo)細(xì)胞脫水和膜損傷[23-24],還會(huì)改變植物細(xì)胞的正常生理代謝過(guò)程,主要表現(xiàn)在2個(gè)方面:呼吸作用障礙[25]、活性氧的積累和去除[26-28]。萱草是一種在許多地區(qū)廣泛種植的營(yíng)養(yǎng)植物,耐寒品種的選育對(duì)西北地區(qū)的營(yíng)養(yǎng)植物生長(zhǎng)具有重要意義。低溫脅迫通過(guò)抑制生長(zhǎng)和發(fā)育嚴(yán)重限制了植物物種的地理分布,植物可以通過(guò)調(diào)節(jié)基因表達(dá)來(lái)適應(yīng)低溫脅迫??沟蜏孛{迫的植物可以通過(guò)防御機(jī)制調(diào)節(jié)自身,包括生理反應(yīng)、代謝和信號(hào)轉(zhuǎn)導(dǎo)途徑的適應(yīng),以及基因表達(dá)變化[29]。本研究對(duì)2種不同抗性萱草在低溫脅迫條件下的轉(zhuǎn)錄組進(jìn)行分析發(fā)現(xiàn),在7個(gè)數(shù)據(jù)庫(kù)同時(shí)注釋的unigenes共有6 727條,至少有1種數(shù)據(jù)庫(kù)注釋成功的unigenes共有56 679條,其

    qRT-PCR結(jié)果中的誤差線表示平均值的標(biāo)準(zhǔn)差(n=3)。The error line in results represents the standard deviation of the mean value (n=3).圖8 差異表達(dá)基因的qRT-PCR驗(yàn)證Fig.8 qRT-PCR validation for differentially expressed genes identified by RNA-Seq

    比例達(dá)到65.09%;還有30 394條unigenes未獲得注釋,這一比例和多種植物的測(cè)定結(jié)果相類似[30-32],說(shuō)明萱草轉(zhuǎn)錄組中存在著大量特征序列,以及大量unigenes功能有待開(kāi)發(fā)。與對(duì)照組相比,冷敏感型萱草根莖在低溫脅迫下,上調(diào)表達(dá)基因有1 849條,下調(diào)表達(dá)基因有4 304條;而耐冷萱草根莖在低溫脅迫下,有6 786條基因上調(diào)表達(dá),8 027條基因下調(diào)表達(dá),差異表達(dá)基因數(shù)量顯著多于冷敏感型萱草品種,猜測(cè)其抗凍性可能與迅速調(diào)控多種基因差異表達(dá)有關(guān)。差異表達(dá)基因的GO分類結(jié)果顯示,2種萱草根莖的DEG功能特性與生物過(guò)程、分子功能和細(xì)胞組分相關(guān)。2種萱草低溫脅迫和非低溫脅迫條件下差異表達(dá)基因均涉及KEGG代謝通路的五大分支。

    次生代謝物在植物逆境脅迫下的合成和積累具有重要意義。已有研究發(fā)現(xiàn),低溫等逆境脅迫能夠誘導(dǎo)類黃酮化合物的積累,提高植物抗逆性[33];一些植物在低溫脅迫下體內(nèi)類黃酮的生物合成會(huì)發(fā)生顯著變化[34]。本研究也發(fā)現(xiàn),低溫脅迫下萱草根莖的類黃酮合成基因表達(dá)水平顯著高于對(duì)照組,這說(shuō)明低溫脅迫在一定程度上可能誘導(dǎo)了類黃酮的積累?;ㄇ嗨?、吲哚生物堿等物質(zhì)的合成也會(huì)受到低溫脅迫的影響。

    低溫刺激植物體內(nèi)Ca2+信號(hào)通路的變化,Ca2+信號(hào)通路被三大類蛋白質(zhì)感知:CDPKs、CaMs和CBLs[35]。在水稻中,低溫脅迫下,鈣調(diào)蛋白OsMSR2基因在不同發(fā)育階段不同組織均顯著上調(diào)表達(dá);擬南芥中外源表達(dá)該基因可以影響其耐鹽性和耐旱性[36]。在耐冷型萱草中,2個(gè)鈣調(diào)蛋白基因Calm5/Calm3均顯著上調(diào)表達(dá),而在冷敏感型萱草中,沒(méi)有發(fā)現(xiàn)相應(yīng)基因的變化。在低溫處理的茶樹(shù)中,5個(gè)鈣調(diào)蛋白基因、2個(gè)CDPK基因和1個(gè)CBL基因參與了信號(hào)轉(zhuǎn)導(dǎo)[37]。本研究中,CDPK3基因表達(dá)在冷敏感型萱草中顯著上調(diào),耐冷型萱草中CDPK14基因表達(dá)顯著下調(diào)。Ca2+信號(hào)通路相關(guān)基因表達(dá)在2種萱草中差異較顯著,其中calmodulin and titin-interacting RhoGEF(CL7725.Contig5_All)與calmodulin-binding transcription activator 3(Unigene6945_All)基因均在低溫脅迫的耐冷萱草中被顯著抑制,而在冷敏感萱草中無(wú)明顯差異,猜測(cè)這2個(gè)基因在Ca2+信號(hào)通路對(duì)萱草耐冷調(diào)控中發(fā)揮重要作用。MAPKs信號(hào)通路也影響植物非生物脅迫響應(yīng),許多受體蛋白激酶(RLKs)負(fù)責(zé)感知外部環(huán)境的變化和相應(yīng)的信號(hào)轉(zhuǎn)導(dǎo)[38]。在冷敏感型萱草中,MAPKKKK3上調(diào)表達(dá),MAPKKKK5下調(diào)表達(dá);在耐冷型萱草中也發(fā)現(xiàn)這2個(gè)基因有相同的表達(dá)趨勢(shì),且上下調(diào)程度更顯著。冷敏感型萱草中僅發(fā)現(xiàn)1個(gè)MAPK1基因下調(diào)表達(dá),而耐冷型萱草中共檢測(cè)到3個(gè)差異表達(dá)的MAPK基因(MAPK1、MAPK4、MAPK7-like)和1個(gè)差異表達(dá)的MLK基因。這些發(fā)現(xiàn)為Ca2+在萱草低溫響應(yīng)中的關(guān)鍵作用提供了證據(jù),同時(shí)進(jìn)一步證明MAPK通路影響萱草對(duì)低溫脅迫的響應(yīng)。MAPKs信號(hào)通路基因MAPK1(CL4304.Contig1_All)的表達(dá)受低溫抑制,在耐冷型萱草中的抑制程度明顯高于冷敏感萱草,表明MAPK1基因可能作為MAPKs信號(hào)通路的關(guān)鍵基因發(fā)揮作用。

    轉(zhuǎn)錄因子在植物非生物脅迫應(yīng)答中發(fā)揮重要作用,在擬南芥中過(guò)量表達(dá)水稻MYB4(R2R3-型MYB轉(zhuǎn)錄因子)和OsMYB3R-2(R1R2R3型MYB轉(zhuǎn)錄因子)基因,可顯著提高擬南芥抗凍性[39]。冷敏感型萱草中,MYB家族轉(zhuǎn)錄因子MYB102顯著下調(diào)表達(dá),該基因在耐冷型萱草中下調(diào)倍數(shù)較低。此外,另一個(gè)MYB家族轉(zhuǎn)錄因子MYB20僅在耐冷型萱草中上調(diào)表達(dá),在冷敏感型萱草未發(fā)現(xiàn)該基因表達(dá)變化。研究表明,WRKY和bHLH轉(zhuǎn)錄因子家族與植物低溫脅迫密切相關(guān)[40-41]。bHLH轉(zhuǎn)錄因子家族基因bHLH2僅在耐冷型萱草低溫處理后上調(diào)表達(dá),而WRKY轉(zhuǎn)錄因子家族基因WRKY11在冷敏感型萱草中下調(diào)表達(dá);與冷敏感型萱草相比,耐冷萱草中的WRKY11在常溫和低溫條件表達(dá)水平均較高,WRKY11在2種萱草中的表達(dá)差異及其在應(yīng)對(duì)低溫脅迫中的作用機(jī)理需要進(jìn)一步探討。以上結(jié)果表明,MYB、WRKY、bHLH轉(zhuǎn)錄因子家族在冷敏感型萱草和耐冷型萱草低溫脅迫響應(yīng)中發(fā)揮了重要而不同的作用。除了蛋白激酶和轉(zhuǎn)錄因子,一些功能蛋白也參與了低溫脅迫過(guò)程。熱休克蛋白(HSPs)最初是通過(guò)參與對(duì)高溫的反應(yīng)而被識(shí)別的,越來(lái)越多的研究發(fā)現(xiàn),許多熱休克蛋白也對(duì)低溫有反應(yīng)。在冷敏感型萱草和耐冷型萱草中均發(fā)現(xiàn)一些HSP蛋白基因呈現(xiàn)差異表達(dá),但在2種萱草中篩選到的HSP家族成員不同。有趣的是,在2種萱草中均發(fā)生變化的HSP90,在冷敏感型萱草中上調(diào)表達(dá),而在耐冷型萱草中則下調(diào)表達(dá),推測(cè)HSP90可能是萱草響應(yīng)低溫脅迫的重要候選基因,其在萱草低溫脅迫中的作用還需要進(jìn)一步研究。

    綜上所述,本研究利用高通量測(cè)序鑒定和比較了萱草Jinyan 7(耐冷型品種)和Lucretius 4(冷敏感型品種)的低溫響應(yīng)基因,發(fā)現(xiàn)在低溫脅迫下,Ca2+信號(hào)通路基因、MAPKs信號(hào)通路基因、轉(zhuǎn)錄因子和熱激蛋白基因在2種萱草中表現(xiàn)出不同的變化趨勢(shì),這些關(guān)鍵路徑基因表達(dá)的差異可能導(dǎo)致了2種萱草在低溫脅迫中的不同反應(yīng)。本研究為萱草低溫響應(yīng)的分子機(jī)制提供了思路,同時(shí),篩選出的差異基因也可能成為培育耐冷型萱草品種的候選基因。

    猜你喜歡
    萱草差異基因低溫
    低溫也能“燙傷”嗎
    ICR鼠肝和腎毒性損傷生物標(biāo)志物的篩選
    基于低溫等離子體修飾的PET/PVC浮選分離
    四月節(jié)令百草之萱草 誰(shuí)言四月芳菲盡 萱草葵花覺(jué)昨非
    紫禁城(2020年5期)2021-01-07 02:13:34
    基于RNA 測(cè)序研究人參二醇對(duì)大鼠心血管內(nèi)皮細(xì)胞基因表達(dá)的影響 (正文見(jiàn)第26 頁(yè))
    零下低溫引發(fā)的火災(zāi)
    北堂種萱草
    海峽姐妹(2016年7期)2016-02-27 15:21:17
    秋葉有情,萱草有愛(ài)
    人間(2015年19期)2016-01-04 12:46:44
    不同品種大花萱草抗寒性比較
    低溫休眠不是夢(mèng)
    天堂网av新在线| 色综合婷婷激情| 欧美黑人巨大hd| 欧美午夜高清在线| 国产不卡一卡二| av国产免费在线观看| 国产成人系列免费观看| 久久婷婷人人爽人人干人人爱| 90打野战视频偷拍视频| 一级作爱视频免费观看| 99国产精品99久久久久| 日本免费一区二区三区高清不卡| 色综合站精品国产| 久久这里只有精品19| 日本一本二区三区精品| 欧美三级亚洲精品| 长腿黑丝高跟| 可以在线观看的亚洲视频| 午夜久久久久精精品| 久久精品国产综合久久久| 999久久久国产精品视频| 亚洲18禁久久av| 又大又爽又粗| 夜夜躁狠狠躁天天躁| 亚洲精品在线美女| 美女高潮喷水抽搐中文字幕| 一本一本综合久久| 日日夜夜操网爽| 女同久久另类99精品国产91| 亚洲片人在线观看| 日韩av在线大香蕉| 精品久久久久久久久久久久久| 欧美激情久久久久久爽电影| 国产免费av片在线观看野外av| 亚洲国产精品成人综合色| 狂野欧美激情性xxxx| 欧美乱码精品一区二区三区| 18禁黄网站禁片免费观看直播| 国内毛片毛片毛片毛片毛片| 成年免费大片在线观看| 亚洲欧美日韩东京热| 亚洲国产精品sss在线观看| 俄罗斯特黄特色一大片| www国产在线视频色| 久久久久免费精品人妻一区二区| 国产一区二区在线观看日韩 | 制服丝袜大香蕉在线| 在线看三级毛片| 99精品在免费线老司机午夜| 亚洲午夜理论影院| 色综合欧美亚洲国产小说| 嫩草影院精品99| 亚洲av美国av| 精品国产三级普通话版| 日韩人妻高清精品专区| 亚洲 欧美 日韩 在线 免费| 白带黄色成豆腐渣| 国产精品一区二区精品视频观看| 亚洲人成网站在线播放欧美日韩| 亚洲中文字幕一区二区三区有码在线看 | 欧美日韩国产亚洲二区| av天堂在线播放| 淫秽高清视频在线观看| 亚洲成a人片在线一区二区| 亚洲av成人一区二区三| 欧美在线黄色| 欧美成人性av电影在线观看| 色综合欧美亚洲国产小说| 又黄又爽又免费观看的视频| 国内揄拍国产精品人妻在线| 日本一二三区视频观看| 亚洲精华国产精华精| e午夜精品久久久久久久| 欧美乱码精品一区二区三区| 男女下面进入的视频免费午夜| 两个人视频免费观看高清| 久久精品国产99精品国产亚洲性色| 欧美日韩国产亚洲二区| 亚洲乱码一区二区免费版| www.999成人在线观看| 欧美日韩瑟瑟在线播放| 好看av亚洲va欧美ⅴa在| 真人做人爱边吃奶动态| 日本一本二区三区精品| 老司机深夜福利视频在线观看| 一级黄色大片毛片| 在线视频色国产色| 丰满人妻一区二区三区视频av | 成人永久免费在线观看视频| 99久久精品一区二区三区| 中文字幕久久专区| 欧美高清成人免费视频www| 中文字幕av在线有码专区| netflix在线观看网站| 欧美日韩瑟瑟在线播放| 全区人妻精品视频| av天堂在线播放| 国产私拍福利视频在线观看| 噜噜噜噜噜久久久久久91| 国模一区二区三区四区视频 | 变态另类丝袜制服| 老司机午夜福利在线观看视频| www日本在线高清视频| 一区福利在线观看| 久久久色成人| 国产伦一二天堂av在线观看| 国产精品一及| 很黄的视频免费| 国产精品亚洲美女久久久| 一级a爱片免费观看的视频| 老鸭窝网址在线观看| 午夜久久久久精精品| 老汉色∧v一级毛片| 国产精品日韩av在线免费观看| 免费在线观看亚洲国产| 一二三四在线观看免费中文在| 午夜影院日韩av| 亚洲精品美女久久久久99蜜臀| 国产91精品成人一区二区三区| 久久99热这里只有精品18| 色吧在线观看| 国产精品自产拍在线观看55亚洲| 久久亚洲真实| 岛国在线免费视频观看| 日本一二三区视频观看| 很黄的视频免费| 久久久国产成人免费| 国产成人av教育| 一卡2卡三卡四卡精品乱码亚洲| 蜜桃久久精品国产亚洲av| 亚洲男人的天堂狠狠| 国产成人欧美在线观看| 97超级碰碰碰精品色视频在线观看| 精品无人区乱码1区二区| 少妇的丰满在线观看| 欧美极品一区二区三区四区| 在线播放国产精品三级| 国产精品 欧美亚洲| 怎么达到女性高潮| 99热这里只有是精品50| 亚洲av成人一区二区三| 久久精品亚洲精品国产色婷小说| 白带黄色成豆腐渣| 在线观看66精品国产| av在线蜜桃| av欧美777| 亚洲九九香蕉| 欧美日韩一级在线毛片| 国产人伦9x9x在线观看| 又黄又粗又硬又大视频| 此物有八面人人有两片| 日本 欧美在线| 日韩大尺度精品在线看网址| 国产野战对白在线观看| 男插女下体视频免费在线播放| 窝窝影院91人妻| 亚洲欧美日韩无卡精品| 1024手机看黄色片| 手机成人av网站| 99精品欧美一区二区三区四区| 很黄的视频免费| 亚洲av五月六月丁香网| 视频区欧美日本亚洲| 成熟少妇高潮喷水视频| 亚洲 国产 在线| 国产高清有码在线观看视频| www.熟女人妻精品国产| 亚洲色图av天堂| 国产麻豆成人av免费视频| xxxwww97欧美| 人妻丰满熟妇av一区二区三区| 国产精品久久久久久久电影 | 欧美日韩瑟瑟在线播放| 国产乱人视频| 悠悠久久av| 久久精品国产清高在天天线| 十八禁人妻一区二区| 久久国产精品人妻蜜桃| 亚洲电影在线观看av| 久久精品影院6| 国产精品女同一区二区软件 | 制服人妻中文乱码| 999久久久精品免费观看国产| 国产欧美日韩精品亚洲av| 一级作爱视频免费观看| 国产精品电影一区二区三区| 国产亚洲av高清不卡| 国产麻豆成人av免费视频| 亚洲av电影不卡..在线观看| 18禁观看日本| 成人无遮挡网站| 国产av麻豆久久久久久久| 啦啦啦免费观看视频1| 男插女下体视频免费在线播放| 国产亚洲精品av在线| 男女做爰动态图高潮gif福利片| 久9热在线精品视频| 熟妇人妻久久中文字幕3abv| 哪里可以看免费的av片| 欧美国产日韩亚洲一区| 成人三级黄色视频| 亚洲人与动物交配视频| 全区人妻精品视频| 亚洲狠狠婷婷综合久久图片| xxxwww97欧美| 久久午夜综合久久蜜桃| 男人的好看免费观看在线视频| 国产亚洲av嫩草精品影院| 熟女少妇亚洲综合色aaa.| 一本综合久久免费| 免费看a级黄色片| 成年女人永久免费观看视频| 麻豆成人午夜福利视频| 国产精品一区二区免费欧美| 国产精品98久久久久久宅男小说| 色综合亚洲欧美另类图片| 午夜激情福利司机影院| 麻豆成人av在线观看| 噜噜噜噜噜久久久久久91| 亚洲欧洲精品一区二区精品久久久| 嫩草影视91久久| 欧美xxxx黑人xx丫x性爽| 久久久久久久久中文| 国产免费av片在线观看野外av| 日韩有码中文字幕| 变态另类成人亚洲欧美熟女| 国产伦人伦偷精品视频| 一区二区三区国产精品乱码| 一二三四社区在线视频社区8| 99久国产av精品| 一区二区三区高清视频在线| 91久久精品国产一区二区成人 | 黄色丝袜av网址大全| 中文字幕av在线有码专区| 亚洲一区二区三区色噜噜| 九九久久精品国产亚洲av麻豆 | 久久久国产成人精品二区| 国产成人av激情在线播放| 久久国产精品影院| 色在线成人网| 亚洲国产欧美人成| 制服丝袜大香蕉在线| 日本撒尿小便嘘嘘汇集6| 亚洲欧美精品综合久久99| 91在线观看av| ponron亚洲| 一进一出好大好爽视频| 性色av乱码一区二区三区2| x7x7x7水蜜桃| 淫妇啪啪啪对白视频| 国产精品1区2区在线观看.| 国产真人三级小视频在线观看| 免费观看人在逋| 久久久国产成人精品二区| 亚洲自偷自拍图片 自拍| 人人妻,人人澡人人爽秒播| 黄色视频,在线免费观看| 观看免费一级毛片| 在线永久观看黄色视频| 亚洲av电影不卡..在线观看| 色在线成人网| 操出白浆在线播放| av在线蜜桃| 色精品久久人妻99蜜桃| 欧美激情久久久久久爽电影| 午夜亚洲福利在线播放| 黑人操中国人逼视频| 亚洲av免费在线观看| 欧美另类亚洲清纯唯美| 色在线成人网| 午夜激情福利司机影院| 在线十欧美十亚洲十日本专区| www日本黄色视频网| 叶爱在线成人免费视频播放| 嫁个100分男人电影在线观看| 色av中文字幕| 久久精品亚洲精品国产色婷小说| 日韩成人在线观看一区二区三区| 欧美一区二区国产精品久久精品| 一边摸一边抽搐一进一小说| 精品午夜福利视频在线观看一区| 99国产精品一区二区三区| 亚洲精品在线观看二区| 噜噜噜噜噜久久久久久91| av片东京热男人的天堂| 黑人操中国人逼视频| 天天躁狠狠躁夜夜躁狠狠躁| 性色av乱码一区二区三区2| 亚洲精品乱码久久久v下载方式 | 丰满人妻一区二区三区视频av | 18禁黄网站禁片免费观看直播| 性色av乱码一区二区三区2| 熟女少妇亚洲综合色aaa.| 国语自产精品视频在线第100页| 久久久国产欧美日韩av| 国产黄色小视频在线观看| 最新在线观看一区二区三区| 99久久99久久久精品蜜桃| 无遮挡黄片免费观看| 男人和女人高潮做爰伦理| 麻豆成人av在线观看| 亚洲av熟女| 国产伦一二天堂av在线观看| 免费看美女性在线毛片视频| 亚洲午夜理论影院| 舔av片在线| 国产1区2区3区精品| 国产成年人精品一区二区| 国产精品永久免费网站| 久9热在线精品视频| 亚洲人成网站在线播放欧美日韩| 国产精品久久视频播放| 少妇熟女aⅴ在线视频| av在线天堂中文字幕| 国产激情偷乱视频一区二区| 国产午夜精品久久久久久| 国产久久久一区二区三区| 男人舔女人的私密视频| 久久久久久国产a免费观看| 国产免费av片在线观看野外av| 可以在线观看的亚洲视频| 免费无遮挡裸体视频| 亚洲av成人av| 夜夜夜夜夜久久久久| 麻豆成人av在线观看| 午夜福利免费观看在线| 91麻豆av在线| 全区人妻精品视频| 琪琪午夜伦伦电影理论片6080| 老司机午夜福利在线观看视频| 国产精品一区二区精品视频观看| 国产精品免费一区二区三区在线| 亚洲成人久久爱视频| 亚洲va日本ⅴa欧美va伊人久久| 男女下面进入的视频免费午夜| 国产亚洲精品久久久com| 天天一区二区日本电影三级| 亚洲,欧美精品.| 在线看三级毛片| 免费av不卡在线播放| 国产精品一区二区三区四区久久| 久久午夜亚洲精品久久| 成人欧美大片| 一a级毛片在线观看| 男人和女人高潮做爰伦理| 午夜日韩欧美国产| 国产一区二区在线观看日韩 | 啦啦啦免费观看视频1| 男人和女人高潮做爰伦理| 国产伦在线观看视频一区| 色综合亚洲欧美另类图片| 麻豆国产av国片精品| 91九色精品人成在线观看| 国产精品久久久人人做人人爽| 别揉我奶头~嗯~啊~动态视频| 一个人看视频在线观看www免费 | 18禁国产床啪视频网站| 别揉我奶头~嗯~啊~动态视频| 国产免费男女视频| 精华霜和精华液先用哪个| 午夜福利在线观看吧| 国产激情偷乱视频一区二区| 欧美av亚洲av综合av国产av| a在线观看视频网站| 波多野结衣巨乳人妻| 久久精品夜夜夜夜夜久久蜜豆| 欧美日韩一级在线毛片| 久久香蕉精品热| 国产精品久久视频播放| 亚洲国产高清在线一区二区三| 欧美zozozo另类| 国产av麻豆久久久久久久| 国产一区二区在线观看日韩 | 婷婷六月久久综合丁香| 熟女人妻精品中文字幕| 欧美精品啪啪一区二区三区| 亚洲,欧美精品.| 国产亚洲欧美98| 成年女人看的毛片在线观看| 中国美女看黄片| 欧美日韩综合久久久久久 | 久久久精品大字幕| 久久欧美精品欧美久久欧美| 国产精品 国内视频| 免费一级毛片在线播放高清视频| 国产精品久久久av美女十八| 欧美成狂野欧美在线观看| 日本黄色视频三级网站网址| 久久久国产成人免费| 久久久久久九九精品二区国产| 欧美日韩一级在线毛片| 成人国产一区最新在线观看| 国产69精品久久久久777片 | 99热这里只有是精品50| 女生性感内裤真人,穿戴方法视频| 亚洲av中文字字幕乱码综合| 90打野战视频偷拍视频| 制服人妻中文乱码| 国产男靠女视频免费网站| 男女视频在线观看网站免费| 啦啦啦韩国在线观看视频| 久久精品夜夜夜夜夜久久蜜豆| 欧美日韩福利视频一区二区| 亚洲av五月六月丁香网| a在线观看视频网站| 精品一区二区三区视频在线观看免费| 99精品久久久久人妻精品| 男人舔女人的私密视频| 首页视频小说图片口味搜索| 在线观看一区二区三区| 九色国产91popny在线| 亚洲精品美女久久久久99蜜臀| 美女被艹到高潮喷水动态| 一个人免费在线观看的高清视频| 国产激情偷乱视频一区二区| 国产私拍福利视频在线观看| 五月伊人婷婷丁香| 精品免费久久久久久久清纯| 亚洲第一电影网av| 国产高潮美女av| 国产精品精品国产色婷婷| 日韩欧美国产在线观看| 两个人的视频大全免费| 制服丝袜大香蕉在线| 99精品在免费线老司机午夜| 久久亚洲精品不卡| 桃色一区二区三区在线观看| 日本 av在线| 亚洲成人久久爱视频| 成人国产一区最新在线观看| 99久国产av精品| 男女做爰动态图高潮gif福利片| 欧美大码av| 变态另类丝袜制服| 国产成+人综合+亚洲专区| 91av网站免费观看| 亚洲精品一卡2卡三卡4卡5卡| 国产精品 国内视频| 中文字幕av在线有码专区| 国产伦精品一区二区三区视频9 | 熟女人妻精品中文字幕| 天天一区二区日本电影三级| 国产一区二区在线av高清观看| 日韩欧美在线二视频| 国产野战对白在线观看| 黄色片一级片一级黄色片| 不卡一级毛片| 国产成人影院久久av| 国产黄色小视频在线观看| 国产伦精品一区二区三区四那| 亚洲成人中文字幕在线播放| 亚洲色图 男人天堂 中文字幕| 日本熟妇午夜| 最近视频中文字幕2019在线8| 成人精品一区二区免费| 男女之事视频高清在线观看| 在线a可以看的网站| 麻豆国产av国片精品| 亚洲真实伦在线观看| 老司机福利观看| 18禁国产床啪视频网站| 久久久久久久久中文| 日日干狠狠操夜夜爽| 国产亚洲精品久久久久久毛片| 在线观看免费午夜福利视频| 久久国产乱子伦精品免费另类| 国产久久久一区二区三区| 波多野结衣高清作品| 看黄色毛片网站| 岛国在线免费视频观看| 国产午夜福利久久久久久| 热99re8久久精品国产| 亚洲成人中文字幕在线播放| 亚洲avbb在线观看| 一二三四社区在线视频社区8| 亚洲成av人片在线播放无| 亚洲aⅴ乱码一区二区在线播放| 一区二区三区激情视频| 很黄的视频免费| 午夜a级毛片| 夜夜躁狠狠躁天天躁| 亚洲国产欧洲综合997久久,| 亚洲av成人精品一区久久| 国产乱人伦免费视频| 亚洲成人精品中文字幕电影| 操出白浆在线播放| 一区二区三区激情视频| 亚洲精品久久国产高清桃花| 精品福利观看| 亚洲av熟女| 在线观看日韩欧美| 国产激情偷乱视频一区二区| 老汉色av国产亚洲站长工具| 久久人妻av系列| 99久久精品热视频| 欧洲精品卡2卡3卡4卡5卡区| 久久国产乱子伦精品免费另类| e午夜精品久久久久久久| 亚洲精品一区av在线观看| 国产精品美女特级片免费视频播放器 | 成人特级黄色片久久久久久久| bbb黄色大片| 国产不卡一卡二| 日韩大尺度精品在线看网址| 精品免费久久久久久久清纯| 999久久久精品免费观看国产| 国产亚洲欧美98| 亚洲激情在线av| 看黄色毛片网站| 老鸭窝网址在线观看| 国产高潮美女av| 欧美成人免费av一区二区三区| 国产精品一区二区三区四区久久| 免费在线观看视频国产中文字幕亚洲| 午夜福利视频1000在线观看| 国产av在哪里看| 色尼玛亚洲综合影院| 亚洲精品国产精品久久久不卡| 岛国在线观看网站| av女优亚洲男人天堂 | 国产精品一区二区免费欧美| 国内毛片毛片毛片毛片毛片| 欧美三级亚洲精品| 欧美日韩一级在线毛片| 神马国产精品三级电影在线观看| 禁无遮挡网站| 欧美中文综合在线视频| 悠悠久久av| 在线观看日韩欧美| 黄色 视频免费看| 精品一区二区三区视频在线 | 手机成人av网站| 老鸭窝网址在线观看| 国产高潮美女av| 校园春色视频在线观看| 亚洲国产欧美人成| 免费av不卡在线播放| x7x7x7水蜜桃| 美女黄网站色视频| 免费在线观看视频国产中文字幕亚洲| 久久久色成人| 一a级毛片在线观看| 丁香六月欧美| 午夜福利免费观看在线| 香蕉久久夜色| 在线观看美女被高潮喷水网站 | 亚洲精品美女久久久久99蜜臀| 久久久久久久久免费视频了| 99在线人妻在线中文字幕| 国产精品野战在线观看| 欧美一区二区国产精品久久精品| 国产激情欧美一区二区| 亚洲中文av在线| 9191精品国产免费久久| 精品久久久久久久末码| 婷婷精品国产亚洲av在线| 级片在线观看| 制服丝袜大香蕉在线| 黄色 视频免费看| 精品久久久久久成人av| www.999成人在线观看| 99久久99久久久精品蜜桃| 国产高清有码在线观看视频| 麻豆成人av在线观看| 香蕉国产在线看| 精品电影一区二区在线| 免费在线观看成人毛片| 午夜影院日韩av| 真实男女啪啪啪动态图| 免费观看人在逋| 黄色 视频免费看| www.精华液| 天天添夜夜摸| 欧美日韩黄片免| 色视频www国产| 国产又色又爽无遮挡免费看| 69av精品久久久久久| 中文在线观看免费www的网站| 1024香蕉在线观看| 淫妇啪啪啪对白视频| 国产av不卡久久| 亚洲中文日韩欧美视频| 免费观看人在逋| 国产伦一二天堂av在线观看| 99国产精品一区二区蜜桃av| 国产精品98久久久久久宅男小说| a在线观看视频网站| 午夜视频精品福利| x7x7x7水蜜桃| 18禁黄网站禁片免费观看直播| 国内久久婷婷六月综合欲色啪| 长腿黑丝高跟| 一二三四在线观看免费中文在| 最新在线观看一区二区三区| АⅤ资源中文在线天堂| 变态另类成人亚洲欧美熟女| 又爽又黄无遮挡网站| 亚洲av免费在线观看| 色综合婷婷激情| 精品免费久久久久久久清纯| av天堂中文字幕网| 午夜福利视频1000在线观看| 午夜福利在线在线| 亚洲在线观看片| 嫩草影院入口| 色在线成人网| 国产伦精品一区二区三区四那| 18禁裸乳无遮挡免费网站照片| 国产三级中文精品| 亚洲avbb在线观看| 老汉色av国产亚洲站长工具| 黑人操中国人逼视频| 国产成人精品无人区| 久久精品91蜜桃| 久久久久性生活片|