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

    油菜鹽脅迫響應(yīng)miRNA及其靶基因鑒定

    2022-03-16 10:31:58王越周婷岳彩鵬黃進(jìn)勇華營(yíng)鵬
    關(guān)鍵詞:文庫(kù)擬南芥抗性

    王越,周婷,岳彩鵬,黃進(jìn)勇,華營(yíng)鵬

    (鄭州大學(xué)農(nóng)學(xué)院,河南 鄭州,450001)

    土壤鹽堿化是全球重要的環(huán)境和農(nóng)業(yè)問題。據(jù)估計(jì),鹽脅迫可能影響全球一半的可耕地。土壤鹽脅迫嚴(yán)重抑制了作物生長(zhǎng)發(fā)育并顯著降低了作物的產(chǎn)量[1],其主要原因是鹽脅迫會(huì)導(dǎo)致高滲透脅迫、氧化應(yīng)激、營(yíng)養(yǎng)失衡、離子毒性、生物膜紊亂,并嚴(yán)重破壞細(xì)胞分裂及關(guān)鍵代謝等過程,導(dǎo)致水勢(shì)下降,從而限制水分的吸收和細(xì)胞的擴(kuò)張[2]。為了在鹽脅迫條件下生存下來,植物已經(jīng)進(jìn)化出各種復(fù)雜的機(jī)制來緩解鹽脅迫對(duì)自身造成的損傷[3,4]。因此,了解鹽脅迫的分子機(jī)制對(duì)于提高作物耐鹽性具有重要意義[5]。

    MicroRNA(miRNA)是一類由內(nèi)源基因編碼的長(zhǎng)度約為18~24 nt 的非編碼單鏈RNA 分子,在轉(zhuǎn)錄后和翻譯水平負(fù)調(diào)控基因表達(dá)[6]。越來越多的證據(jù)表明,miRNA 參與了植物對(duì)非生物脅迫的響應(yīng)[6~10],如溫度[11]、干旱[12]和鹽[13,14]脅迫。在擬南芥鹽脅迫反應(yīng)中,miR156、miR167、miR168 和miR396 等多個(gè)miRNA 存在差異表達(dá)[15]。過表達(dá)水稻miR319a的匍匐禾草(Agrostis stolonifera)表現(xiàn)出更強(qiáng)的耐鹽性[16],其原因是鹽脅迫誘導(dǎo)miR319 表達(dá)升高,導(dǎo)致至少4個(gè)miR319 的靶基因(AsPCF5、AsPCF6、AsPCF8和AsTCP14)和AsNAC60的表達(dá)下調(diào),從而提高了植株的鹽脅迫抗性。過表達(dá)miR396和miR394的棉花和擬南芥表現(xiàn)出高度的鹽脅迫敏感性[17~19],其原因是miR396 和miR394 可能作為一種負(fù)調(diào)控因子,分別以GRF 和LCR 等調(diào)控蛋白為靶點(diǎn)。其中GRF 抑制擬南芥的應(yīng)激反應(yīng)基因AtDREB2A的表達(dá),增加了植物對(duì)滲透脅迫的耐受性,LCR 則以ABA 依賴性方式參與擬南芥鹽脅迫抗性反應(yīng)。鹽脅迫激活WRKY并調(diào)節(jié)植物中脅迫耐受性相關(guān)途徑越來越受關(guān)注[15]。前人研究表明,擬南芥RPD3-like 組蛋白去乙酰化酶HDA9 通過調(diào)節(jié)WRKY53 的DNA 結(jié)合能力和轉(zhuǎn)錄活性來降低擬南芥的鹽脅迫抗性[20]。近期研究發(fā)現(xiàn)蘋果中MdSPL13 與MdWRKY100啟動(dòng)子區(qū)結(jié)合,并正向調(diào)節(jié)MdWRKY100的表達(dá),通過促進(jìn)基因轉(zhuǎn)錄,如抗氧化劑生物合成,糖代謝和脯氨酸生物合成有關(guān)基因的轉(zhuǎn)錄,進(jìn)一步維持ROS 和滲透平衡,從而增強(qiáng)鹽脅迫的耐受性[21]。植物在鹽脅迫下的生存需要對(duì)脅迫條件的快速感知和適應(yīng)。前人研究[22]表明RACK1A誘導(dǎo)的miR393 生物發(fā)生,下調(diào)了生長(zhǎng)素信號(hào),可能以此調(diào)節(jié)擬南芥對(duì)鹽度的適應(yīng)。最近研究表明[23]通過表達(dá)TIR1,增加植物生長(zhǎng)素信號(hào)傳導(dǎo),可能會(huì)觸發(fā)植物生長(zhǎng)素介導(dǎo)的下游途徑,以通過滲透壓調(diào)節(jié)和增加Na+外排來增強(qiáng)植物對(duì)鹽脅迫的抗性。這一發(fā)現(xiàn)使我們更好地理解miR393 介導(dǎo)生長(zhǎng)素調(diào)節(jié)的鹽脅迫反應(yīng)。上述研究結(jié)果均表明,miRNA 在植物的鹽脅迫抗性反應(yīng)中發(fā)揮重要作用。

    油菜(AnAnCnCn,2n=4x=38)是世界上重要的油料作物之一,起源于其二倍體祖先Brassica rapa(ArAr, 2n=2x=20)和Brassica oleracea(CoCo, 2n=2x=18)的自發(fā)雜交[24]。目前,在擬南芥以及甘藍(lán)型油菜等作物中,已有一些耐鹽性相關(guān)的基因被克隆與鑒定。如脯氨酸合成的關(guān)鍵基因P5CS1[25],通過緩解滲透脅迫來參與油菜鹽脅迫響應(yīng);控制油菜離子穩(wěn)態(tài)的Na+/H+逆向轉(zhuǎn)運(yùn)蛋白基因BnNHX1[26],可在液泡中隔離Na+和K+,并通過維持Na+/K+穩(wěn)態(tài)并增強(qiáng)滲透和抗氧化調(diào)節(jié)能力,賦予轉(zhuǎn)基因油菜耐鹽性。還有諸多轉(zhuǎn)錄因子參與調(diào)控耐鹽性,如WRKY家族[27],NAC家族[28]等。盡管前人已經(jīng)對(duì)油菜的鹽脅迫抗性展開了系列研究,但由于油菜基因組龐大復(fù)雜,耐鹽機(jī)制中許多問題仍有待探索。

    本研究的目的是鑒定油菜生長(zhǎng)發(fā)育過程中對(duì)響應(yīng)鹽脅迫的miRNA 及其核心靶基因并研究其在鹽脅迫中發(fā)揮的作用。本研究在200 mmol·L-1NaCl(T)和無鹽處理對(duì)照(C)條件下培養(yǎng)油菜,分別從地上部(S)和根(R)中取樣,提取總RNA 后構(gòu)建了12個(gè)miRNA 文庫(kù)。通過使用高通量測(cè)序技術(shù)分析了12個(gè)miRNA測(cè)序文庫(kù),最終鑒定到了穩(wěn)定差異表達(dá)的26個(gè)miRNA,并根據(jù)與相應(yīng)的轉(zhuǎn)錄組表達(dá)譜的比較預(yù)測(cè)了它們的靶基因。這些鑒定到的差異表達(dá)miRNA 及其核心靶基因有助于深入了解植物鹽脅迫抗性的分子機(jī)制,并為油菜鹽脅迫抗性的遺傳改良提供了優(yōu)異的基因資源。

    1 材料與方法

    1.1 植物材料及試驗(yàn)設(shè)置

    油菜種子為中雙11,由鄭州大學(xué)黃進(jìn)勇教授提供,在25℃黑暗條件下催芽約24 h,挑選露白一致的油菜種子均勻擺放在鋪有濾紙的育苗盤中,放于溫室25℃±3℃培養(yǎng)架上培養(yǎng),幼苗長(zhǎng)到5 cm 左右選取生長(zhǎng)一致的健壯幼苗,移栽到盛有適量霍格蘭營(yíng)養(yǎng)液的黑色塑料盆中水培生長(zhǎng),每5 d 更換一次營(yíng)養(yǎng)液[29]。試驗(yàn)在光照培養(yǎng)室中進(jìn)行,生長(zhǎng)環(huán)境設(shè)定如下:光照強(qiáng)度150 μmol/(m2·s),室溫24℃(白天)/22℃(夜晚),光照周期16 h(光照)/8 h(黑暗),相對(duì)濕度60%。

    全轉(zhuǎn)錄組測(cè)序是在種子萌發(fā)7 d 后,選取長(zhǎng)勢(shì)一致的油菜幼苗,將其培養(yǎng)在不含NaCl的條件下生長(zhǎng)10 d,隨后轉(zhuǎn)移到含200 mmol·L-1NaCl 的溶液中培養(yǎng)12 h 直至取樣。處理組和對(duì)照組油菜幼苗分別取其地上部和根部,每個(gè)部位均包含3 個(gè)生物學(xué)重復(fù)。

    1.2 mRNA文庫(kù)構(gòu)建和測(cè)序

    mRNA 文庫(kù)構(gòu)建、測(cè)序詳情以及本文所用的轉(zhuǎn)錄組數(shù)據(jù)參考本課題組已經(jīng)發(fā)表的文章[30]。本文中重點(diǎn)分析油菜鹽脅迫響應(yīng)的miRNA及其靶基因。

    1.3 sRNA文庫(kù)的構(gòu)建和測(cè)序

    使用TRIzol法提取油菜總RNA,用Nanophotometer?分光光度計(jì)(IMPLEN, CA, USA)檢測(cè)RNA 純度,使用Qubit?RNA 檢測(cè)試劑盒在Qubit?2.0 熒光計(jì)(Life Technologies, CA, USA)中檢測(cè)其濃度。為了評(píng)估RNA 的完整性,使用Agilent Bioanalyzer 2000 system(Agilent Technologies,CA,USA)的RNA Nano 6000 檢測(cè)試劑盒測(cè)定RNA 的完整性值。隨后使用Illumina?的NEBNext?Multiplex 小RNA 文庫(kù)制備試劑盒(NEB,USA)構(gòu)建小RNA 文庫(kù),將索引代碼添加到每個(gè)樣本的屬性序列中,每個(gè)樣品共使用約3.0 μg RNA。在miRNA 的3'和5'端添加特定的索引代碼,然后使用M-MuLV 逆轉(zhuǎn)錄酶(RNase H-)作為催化劑合成第一鏈cDNA。使用LongAmp Taq 2× Master Mix、Illumina 小RNA 引物和index 接頭引物進(jìn)行PCR 擴(kuò)增。PCR 產(chǎn)物經(jīng)8%聚丙烯酰胺凝膠電泳純化,回收140~160 bp 的DNA 片段。使用Agilent Bioanalyzer 2100 system(Agilent Technologies,CA,USA)對(duì)文庫(kù)質(zhì)量進(jìn)行評(píng)估。隨后使用TruSeq SR Kit v3-cBOt-HS(Illumina)在cBot 集群生成系統(tǒng)上對(duì)索引代碼樣本進(jìn)行聚類。集群生成后,在Illumina Hiseq 2500平臺(tái)上對(duì)庫(kù)進(jìn)行測(cè)序,生成50 bp的單端reads。

    1.4 數(shù)據(jù)過濾和讀取

    測(cè)序完成后,首先對(duì)原始reads 進(jìn)行過濾,去除reads中的3’接頭序列、由于接頭自連等原因?qū)е聸]有插入片段的reads、剪切3’端測(cè)序質(zhì)量較低的堿基(質(zhì)量值小于20)、含未知堿基N 的reads、長(zhǎng)度過短的reads(<18 nt)和長(zhǎng)度過長(zhǎng)的reads(>32 nt)以及含有poly A 尾的reads。計(jì)算數(shù)據(jù)的Q20、Q30和GC 含量,以便進(jìn)一步分析。然后,將高質(zhì)量的sRNA 標(biāo)記通過Bowtie 無錯(cuò)配的方式對(duì)比到油菜參考序列上,分析其在油菜參考基因組(中雙11)上的表達(dá)和分布。

    1.5 miRNA分類,靶基因預(yù)測(cè)及其差異表達(dá)分析

    對(duì)過濾后的reads 進(jìn)行分類,分析sRNA 數(shù)據(jù)的組成,選擇篩選18~32 nt的reads 作為后續(xù)分析的有效性數(shù)據(jù),并進(jìn)行該區(qū)間sRNA 長(zhǎng)度分布的統(tǒng)計(jì)[31]。參照miRBase 20.0,使用mirdeep2 和srna-tolls-cli獲得潛在的miRNA 并繪制二級(jí)結(jié)構(gòu)。使用自定義腳本獲得miRNA 計(jì)數(shù)和堿基偏置。整合miREvo 和miRdeep2,通過探索不加注釋的小RNA 標(biāo)簽的二級(jí)結(jié)構(gòu)、Dicer 裂解位點(diǎn)和最小自由能,預(yù)測(cè)新的miRNA。psRobot 中通過psRobot_tar 對(duì)miRNA 的靶基因進(jìn)行預(yù)測(cè)。miRNA 表達(dá)水平通過每百萬轉(zhuǎn)錄量(TPM)來估計(jì)。使用DESeq R 軟件包(1.8.3)將校正后的p值設(shè)為0.05 作為鑒別差異表達(dá)miRNA 的閾值。差異表達(dá)的miRNA 及其靶基因通過Wallenius non-central hyper-geometric distribution 進(jìn)行GO分析。利用KEGG(http://www. kegg. jp/)數(shù)據(jù)庫(kù)對(duì)miRNA 及其靶基因的代謝途徑富集進(jìn)行分析[30,32],并使用KOBAS 來檢測(cè)KEGG 通路中miRNA 及其靶基因的富集情況。

    2 結(jié)果與分析

    2.1 sRNA序列分析

    為了探究油菜鹽脅迫條件下miRNA 的調(diào)控機(jī)制,提取了在200 mmol·L-1鹽濃度下生長(zhǎng)的油菜總RNA 并進(jìn)行了小RNA 測(cè)序。使用Illumina 測(cè)序技術(shù),在對(duì)照組地上部(control shoot,CS)、處理組地上部(treatment shoot,TS)、對(duì)照組根部(control root,CR)和處理組根部(treatment root, TR)小RNA 文庫(kù)中分別獲得約35 704 556、40 977 538、45 958 866 和41 322 651 條raw reads(表1)。對(duì)raw reads 進(jìn)行過濾后,單個(gè)樣品的clean reads 比率、Q20、Q30以及GC含量分別在90%、98%以及39%以上??偣搏@得了164 070 484個(gè)clean reads(表2)。

    表1 油菜文庫(kù)小RNA測(cè)序結(jié)果Table 1 Results of sRNA sequencing of B.napus library

    在本研究的轉(zhuǎn)錄組中,小RNA 序列長(zhǎng)度在18~32 nt 范圍內(nèi),最豐富的是21~24 nt(圖1A)。這些小RNA 由已知的miRNA、預(yù)測(cè)的miRNA、核糖體RNA(rRNAs)、轉(zhuǎn)移核糖核酸(tRNAs)、核小RNA(sn-RNAs)、核仁小RNA(snoRNAs)、反作用小干擾RNA(TAS)和未知的片段組成(表2)。本研究發(fā)現(xiàn),在總的小RNA 中有6.68%~14.34%是已知的miRNA,而73.91%~80.86%是未知的miRNA。未知的miRNA數(shù)量如此之多,表明在油菜中仍有很多的miRNA 有待鑒定。從地上部和根中分別鑒定出2897 和2435個(gè)miRNA(圖1B),其中,不論在地上部還是根中,已知的miRNA 數(shù)量總是少于預(yù)測(cè)的新的miRNA。最豐富的miRNA 序列長(zhǎng)度為24 nt(12.91%),其次是21 nt(8.40%)和23 nt(8.85%)(圖1A),這些長(zhǎng)度與被子植物中典型的miRNA 分布相似。樣品中長(zhǎng)度為24 nt的miRNA 首位堿基以“U”為主,這符合miRNA首位堿基偏向“U”的特征(圖1C,D)。

    表2 小RNA注釋分類Table 2 Annotation classification of small RNA

    圖1 不同處理?xiàng)l件下油菜地上部和根中miRNA堿基長(zhǎng)度分布(A)、數(shù)量(B)、首位堿基偏好性(C)、各位點(diǎn)堿基偏好性(D)Fig.1 Base length distribution of miRNAs(A),numbers of miRNAs(B),preference distribution of the first base of different length miRNAs(C),base preference distribution of miRNAs(D)in shoots and roots of B.napus under different treatments

    2.2 油菜中保守miRNA及新miRNA的鑒定

    為了鑒定油菜地上部和根中的保守miRNA,將小RNA 序列與miRBase 數(shù)據(jù)庫(kù)中已知的植物miRNA序列進(jìn)行比較。通過BLASTN搜索和進(jìn)一步的序列分析,共發(fā)現(xiàn)了26 個(gè)miRNA 家族中的91 個(gè)miRNA。最大的miRNA 家族miR169 有10 個(gè)成員(圖2)。從文庫(kù)中可以看出,這些miRNA 家族的豐度存在明顯差異。為了比較它們表達(dá)水平,計(jì)算了12 個(gè)文庫(kù)中每個(gè)miRNA 的每百萬讀次的平均轉(zhuǎn)錄量(TPM)[33,34]。TPM 標(biāo)準(zhǔn)化后,bna-miR395f 的表達(dá)水平最低,而bna-miR171f的相對(duì)表達(dá)水平最高,其次是bna-miR166b 和bna-miR169j。為了識(shí)別新的miRNA,使用miREvo 和miRDeep2 程序預(yù)測(cè)二級(jí)結(jié)構(gòu)和Dicer 裂解位點(diǎn),并測(cè)量最小自由能??偣矎?2 個(gè)文庫(kù)中預(yù)測(cè)了2116 個(gè)新miRNA。新的成熟miRNA的長(zhǎng)度為18~24 nt,約60%為24 nt。

    圖2 油菜地上部和根中每個(gè)miRNA家族的保守miRNA數(shù)量Fig.2 Number of conserved miRNAs of each miRNA family in shoots and roots of B.napus

    2.3 不同鹽處理下油菜miRNA的差異表達(dá)分析

    為了確定miRNA 的表達(dá)與鹽脅迫的關(guān)系,對(duì)miRNA 的表達(dá)量進(jìn)行了分析,篩選出了表達(dá)量相對(duì)較高的15個(gè)miRNA,并分析它們?cè)谟筒说厣喜亢透械谋磉_(dá)情況(圖3A)。使用fold-change 進(jìn)行計(jì)算,將處理組地上部和根部中miRNA 的表達(dá)與對(duì)照組地上部和根部中miRNA 表達(dá)進(jìn)行比較,規(guī)定|fold change|≥2 且p<0.05 條件下的miRNA 在鹽脅迫條件下存在表達(dá)差異[35,36]。

    與CS文庫(kù)中的表達(dá)相比,TS文庫(kù)中有9個(gè)miRNA 表達(dá)下調(diào),6 個(gè)miRNA 表達(dá)上調(diào)(圖3B),與CR文庫(kù)中的表達(dá)相比,TR 文庫(kù)中有25 個(gè)miRNA 表達(dá)下調(diào),13 個(gè)miRNA 表達(dá)上調(diào)(圖3C)。這些差異表明,鹽脅迫影響著這些miRNA的表達(dá)。

    圖3 不同處理下差異表達(dá)量前15的miRNA在地上部和根中所占的比例(A),地上部處理組和對(duì)照組之間差異表達(dá)miRNA的表達(dá)量(B)和根部處理組和對(duì)照組之間差異表達(dá)miRNA的表達(dá)量(C)Fig.3 Proportion of miRNA in shoots and roots of the top 15 with differential expression under different treatments(A),differential expressions of miRNA between treatment shoot and control shoot(B),and differential expressions of miRNA between control root and treatment root(C)

    2.4 差異表達(dá)miRNA的靶基因數(shù)量分析

    為了進(jìn)一步解釋miRNA 在鹽脅迫反應(yīng)中的作用,我們對(duì)miRNA 的靶基因進(jìn)行了分析,并統(tǒng)計(jì)了每個(gè)差異表達(dá)miRNA 的靶基因數(shù)目,其中靶基因數(shù)目最多的是bna-miR156d,bna-miR156e和bnamiR156f,分別都含有85個(gè)受調(diào)控的靶基因(圖4)。

    圖4 差異表達(dá)miRNA的靶基因數(shù)量Fig.4 Numbers of target genes of miRNAs

    2.5 miRNA靶基因GO注釋和KEGG分析

    在CS 和TS 小RNA 文庫(kù)中,鑒定出17 個(gè)差異表達(dá)的miRNA,共含有93 個(gè)靶基因;其中已知miRNA的靶基因有89 個(gè),新miRNA 的靶基因有4 個(gè)。在CR 和TR 小RNA 文庫(kù)中,鑒定出38 個(gè)差異表達(dá)的miRNA,共含有684個(gè)預(yù)測(cè)靶基因,其中已知miRNA的靶基因有544 個(gè),新miRNA 的靶基因有140 個(gè)。大多數(shù)的預(yù)測(cè)靶基因編碼的一些與脅迫相關(guān)的轉(zhuǎn)錄因子,其中包括SPL、LAC、UBC、NFY以及其它蛋白質(zhì)如AGO2、超氧化物歧化酶等。為了了解miRNA 的生物學(xué)功能,所有推測(cè)的靶基因都通過Blast2GO 軟件進(jìn)行了GO 功能分類,以進(jìn)一步探討miRNA在響應(yīng)鹽脅迫中的作用。

    已知miRNA 和新miRNA 的靶基因被分為三個(gè)GO類別:細(xì)胞組分(CC)、生物學(xué)過程(BP)和分子功能(MF)。圖5(A)是在CS 和TS 小RNA 文庫(kù)中鑒定到的miRNA 靶基因GO 富集圖,圖5(B)是在CR 和TR 小RNA 文庫(kù)中鑒定到的miRNA 靶基因GO 富集圖。CS 和TS 小RNA 文庫(kù)中共有93 個(gè)靶基因參與了1835 個(gè)GO 條目(圖5A)。其中BP 中高度富集的GO 條目是細(xì)胞過程(cellular process);CC 中高度富集的GO 條目是細(xì)胞(cell);MF 中高度富集的GO 條目是結(jié)合活性(binding)。CR 和TR 兩個(gè)小RNA 庫(kù)共有684 個(gè)靶基因參與了4941 生物過程(圖5B)。其中BP 中高度富集的GO 條目是細(xì)胞過程(cellular process);CC 中高度富集的GO 條目是細(xì)胞(cell);MF 中高度富集的GO 條目是結(jié)合活性(binding)。四個(gè)數(shù)據(jù)庫(kù)中都發(fā)現(xiàn)了刺激反應(yīng)(response to stimulus),表明這些靶基因可以在應(yīng)對(duì)鹽脅迫中發(fā)揮作用。已知KEGG 代謝通路被分為7 大類:代謝(Metabolism)、遺傳信息處理(genetic information processing)、環(huán)境信息處理(environmental in-formation processing)、細(xì)胞過程(cellular processes)、生物體系統(tǒng)(organismal systems)、人類疾?。╤uman diseases)、藥物開發(fā)(drug development)。圖6是差異表達(dá)miRNA靶基因參與的pathway通路,其中A圖為CS和TS小RNA 文庫(kù)中鑒定到的差異表達(dá)miRNA 靶基因參與的pathway 通路,B 圖為CR 和TR 小RNA 文庫(kù)中鑒定到的差異表達(dá)miRNA 靶基因參與的pathway 通路。在CS和TS小RNA文庫(kù)中的生物體系統(tǒng)(organismal systems)中有18 個(gè)靶基因參與了免疫防御系統(tǒng)(immune system)。在CR 和TR 兩個(gè)小RNA 庫(kù)在CS 和TS 小RNA 文庫(kù)中的生物體系統(tǒng)(organismal systems)中有21 個(gè)靶基因參與了免疫防御系統(tǒng)(immune system)。上述結(jié)果表明這些靶基因可以積極參與油菜鹽脅迫抗性反應(yīng)。

    圖5 差異表達(dá)miRNA靶基因功能分類Fig.5 Functional classification of different expression miRNA targets

    圖6 差異表達(dá)miRNA靶基因參與的pathway通路Fig.6 Pathway involved in differentially expressed miRNAs target genes

    2.6 油菜地上部和根中差異表達(dá)miRNA的分析

    通過分析油菜地上部和根中miRNA 的差異表達(dá)譜,圖7A 和圖7B分別是鹽脅迫條件下,油菜地上部中上調(diào)表達(dá)和下調(diào)表達(dá)的miRNA;圖7C 和圖7D是鹽脅迫條件下,油菜根中上調(diào)和下調(diào)表達(dá)的miRNA。在地上部中,發(fā)現(xiàn)NaCl處理后bna-miR397a和bna-miR397b表達(dá)顯著上調(diào),bna-miR169a、bnamiR169b和bna-miR169n的表達(dá)顯著下調(diào)。在根中,bna-miR156a、bna-miR156d、bna-miR156e、bnamiR156f以及bna-miR824的表達(dá)顯著上調(diào),而bnamiR169j、bna-miR169i、bna-miR169h、bna-miR169k、bna-miR169g、bna-miR169a和bna-miR169b的表達(dá)并沒有顯著變化,bna-miR169m、bna-miR169n、bnamiR169c、bna-miR169e、bna-miR169f、bna-miR169d和bna-miR169b表達(dá)顯著下調(diào)。

    圖7 鹽脅迫條件下油菜的差異表達(dá)miRNAFig.7 Differentially expressed miRNAs in B.napus under salt stress

    2.7 靶基因表達(dá)分析

    前人研究表明有多個(gè)miRNA 參與植物鹽脅迫反應(yīng)[23,37~41],分別為miR397-LAC(bna-miR397a、bnamiR397b) 、 miR156-SPL(bna-miR156a、bnamiR156d、bna-miR156f、bna-miR156e)、miR169-NFY(bna-miR169m、bna-miR169g、bna-miR169n)以及miR399-UBC(bna-miR399b)。其中轉(zhuǎn)錄因子miR156 的靶基因是SPL[37],在本研究中(圖8)發(fā)現(xiàn)miR156在根中表達(dá)量上調(diào),而在本課題組前期轉(zhuǎn)錄組測(cè)序結(jié)果中發(fā)現(xiàn)SPL15表達(dá)量下調(diào)[30],因此miR156 可能通過調(diào)控SPL15來響應(yīng)鹽脅迫,從而緩解鹽脅迫給植物帶來的損傷。LAC是miR397 的靶基因,其可以通過調(diào)控木質(zhì)素含量來響應(yīng)外界非生物脅迫,差異表達(dá)分析顯示miR397在地上部中表達(dá)上調(diào),LAC12在本課題組前期轉(zhuǎn)錄組測(cè)序結(jié)果中表達(dá)量下調(diào)[30],miRNA397 可能通過調(diào)控LAC12的表達(dá),促進(jìn)原生木質(zhì)部細(xì)胞壁增厚,從而響應(yīng)鹽脅迫。轉(zhuǎn)錄因子NFY是miR169的靶基因[42],它是非生物脅迫反應(yīng)的重要調(diào)節(jié)因子,在本研究中發(fā)現(xiàn)miR169在地上部和根中的表達(dá)量都下調(diào),而在本課題組前期轉(zhuǎn)錄組測(cè)序結(jié)果中NFYA5的表達(dá)量顯著上調(diào)[30],miR169 可能通過調(diào)控NFYA5的表達(dá)參與鹽脅迫反應(yīng)。UBC是miR399 的靶基因,差異表達(dá)分析顯示miR399在根中表達(dá)量下調(diào),而在本課題組前期轉(zhuǎn)錄組測(cè)序結(jié)果中UBC29的表達(dá)量顯著性上調(diào)[30],miR399 可能通過調(diào)控UBC29的表達(dá)從而參與異常蛋白質(zhì)降解過程,提高作物的耐鹽性。

    圖8 鹽脅迫條件下對(duì)miRNA397靶基因LAC的影響(A)、miRNA156靶基因SPL的影響(B)、miRNA169靶基因NFY的影響(C)、miRNA399靶基因UBC的影響(D)Fig.8 Impact on miRNA397 target gene LAC(A),miRNA156 target gene SPL(B),miRNA169 target gene NFY(C),and miRNA399 target gene UBC(D)under salt stress

    3 討論

    鹽漬土使植物暴露在滲透脅迫下,滲透脅迫對(duì)植物最重要的后果之一是產(chǎn)生大量活性氧,隨后是氧化損傷,例如蛋白質(zhì)、脂質(zhì)、色素和DNA 的降解。生長(zhǎng)在含鹽條件下的植物吸收有害離子,尤其是鈉離子和氯離子,其大量積累對(duì)細(xì)胞有毒,并加劇滲透脅迫。這些離子破壞膜完整性、細(xì)胞代謝、酶結(jié)構(gòu)、細(xì)胞生長(zhǎng)和光合作用[43]。近年來,miRNA 在基因表達(dá)調(diào)控中的作用越來越清晰,許多miRNA 如miR156、miR169、miR393、miR396、miR397 及其靶基因[44]等在擬南芥[13]、大豆[45]、大麥[46]等多種植物的鹽脅迫響應(yīng)中發(fā)揮重要作用。

    前人研究表明油菜中有56 個(gè)基因編碼假定的轉(zhuǎn)錄因子在非生物脅迫下會(huì)發(fā)生改變[47]。在這些基因中,已經(jīng)顯示在鹽脅迫下上調(diào)超過5 倍的基因來自AP2-EREBP家族(ATERF11,CBF4/DREB1D,CBF1/DREB1B, ATERF4/RAP2.5),堿性-螺旋-環(huán)-螺旋(bHLH)家族(AtbHLH17),堿性區(qū)亮氨酸拉鏈(bZIP)家族(AtbZIP55/GBF3),C2H2家族(ZAT10、ZAT12、RHL41、ZAT6和ZAT102/RHL41) ,NAC 家 族(ANAC036、ANAC029/ATNAP、ANAC055/ATNAC3、ANAC047、ANAC072/RD26、ANAC002/ATAF1、ANAC019和ANAC032),WRKY家族(ATWRKY53,ATWRKY40和ATWRKY33),熱激家族(ATHSFA1E)以及Homeobox家族(ATHB-7)。此后,由轉(zhuǎn)錄因子和其他蛋白質(zhì)組成的復(fù)雜基因調(diào)控網(wǎng)絡(luò)控制著許多基因的表達(dá)。但是在本研究中,發(fā)現(xiàn)SPL15在鹽脅迫條件下表達(dá)量下調(diào),而LAC12、NFYA5和UBC29在鹽脅迫條件下表達(dá)量顯著上調(diào)。遺憾的是,目前對(duì)油菜鹽度脅迫下miRNA及其作用靶點(diǎn)的研究還不夠全面。

    為了鑒定miRNA 及其靶基因在油菜鹽脅迫響應(yīng)中發(fā)揮的重要作用,在本研究中,鑒定了來自油菜的數(shù)百萬個(gè)小RNA 序列,以了解鹽脅迫條件下miRNA 對(duì)植物生長(zhǎng)發(fā)育的影響。分別在200 mmol·L-1NaCl(T)和無鹽處理對(duì)照(C)條件下培養(yǎng)油菜品種中雙11,取其地上部(S)和根(R)構(gòu)建了12 個(gè)小RNA 文庫(kù),并用Illumina 基因組分析儀對(duì)其進(jìn)行測(cè)序,從CS、TS、CR 和TR 庫(kù)中共得到198 831 780 個(gè)raw reads。進(jìn)一步對(duì)全基因組的miRNA 進(jìn)行差異表達(dá)分析,一共鑒定到了26 個(gè)差異表達(dá)的miRNA。這些差異表達(dá)的miRNA 導(dǎo)致了下游多個(gè)基因的表達(dá)發(fā)生顯著變化,并且植物對(duì)鹽脅迫的響應(yīng)伴隨著一系列廣泛的細(xì)胞內(nèi)過程,包括信號(hào)感知、信號(hào)轉(zhuǎn)導(dǎo)、轉(zhuǎn)錄和蛋白質(zhì)生物合成等。這說明不同的植物通過多種miRNA 介導(dǎo)的調(diào)控策略來響應(yīng)脅迫[48],例如miR159、miR393、miR396 和miR398 等與鹽脅迫有關(guān)[49~52]。

    在本研究中鑒定出miR156-SPL15、miR397-LAC12、miR169-NFYA5和miR399-UBC29在鹽脅迫條件下發(fā)揮作用。有研究表明miR156 及其靶基因SPL在植物非生物脅迫的響應(yīng)中發(fā)揮著重要的作用[53],在本研究中,miR156 特異性的在根中上調(diào)表達(dá),這一結(jié)果與前人的研究結(jié)果形成對(duì)比,即擬南芥在鹽脅迫下miR156 的表達(dá)下調(diào)[21],并通過調(diào)控其靶基因SPL15參與到鹽脅迫響應(yīng)中;蘋果在鹽脅迫條件下miR156的表達(dá)下調(diào),其靶基因SPL13通過激活WRKY100,促進(jìn)與抗氧化劑生物合成,糖代謝和脯氨酸生物合成有關(guān)基因的轉(zhuǎn)錄,進(jìn)一步維持ROS和滲透平衡,從而增強(qiáng)鹽脅迫的耐受性[21]。這表明不同植物對(duì)鹽脅迫的響應(yīng)存在差異[54]。miR397 及其靶基因LAC可以通過調(diào)控木質(zhì)素含量來參與非生物脅迫(58),在本研究中,鑒定到LAC12是miR397a 的靶基因,該基因可以導(dǎo)致木質(zhì)素含量升高,進(jìn)一步改變Na+遷移途徑從而提高植物對(duì)鹽脅迫的抗性。已有研究表明,番茄和玉米根系中高濃度的NaCl均能提高LAC的轉(zhuǎn)錄水平,并指出鹽脅迫下根系轉(zhuǎn)錄水平的提高可能是植物的普遍響應(yīng)。因此,LAC可能在植物根系適應(yīng)鹽脅迫過程中發(fā)揮重要作用[55,56]。越來越多的證據(jù)表明NFY是非生物脅迫反應(yīng)的重要調(diào)節(jié)因子[16,17,43,57~61],并受到miR169 的調(diào)控。在本研究中miR169 家族成員miR169g 和miR169n 在根中的表達(dá)量下調(diào),這和前人之前的研究一致,但是過表達(dá)miR396c 的水稻和擬南芥對(duì)鹽脅迫抗性都降低,而過表達(dá)miR396a/b 轉(zhuǎn)基因植株的耐鹽性提高。這表明同一個(gè)miRNA 家族的不同成員在不同的植物中采用不同的方式參與鹽脅迫響應(yīng)[60]。然而,鹽脅迫條件下油菜中NFY的作用還沒有被詳細(xì)研究。前人研究表明UBC[42,62]定位在內(nèi)質(zhì)網(wǎng)上,是miR399 的靶基因,其表達(dá)受高鹽度或ABA 誘導(dǎo)并參與異常蛋白質(zhì)的降解過程[42],且在鹽脅迫條件下UBC轉(zhuǎn)錄水平升高,通過調(diào)控泛素化途徑參與植物鹽脅迫抗性反應(yīng)。這表明,植物在鹽脅迫反應(yīng)中可能存在著共同的調(diào)控機(jī)制。這些結(jié)果為進(jìn)一步研究作物抵抗鹽脅迫的分子機(jī)制提供了理論依據(jù)。

    4 結(jié)論

    植物的耐鹽過程是復(fù)雜的,由多種因素決定,miRNA 的參與進(jìn)一步增加了這一過程的復(fù)雜性。本研究通過油菜miRNA 測(cè)序以及轉(zhuǎn)錄組組測(cè)序技術(shù)發(fā)掘鹽脅迫響應(yīng)相關(guān)的miRNA及其核心靶基因,豐富了油菜地上部和根部miRNA 的情況。共篩選出26 個(gè)與鹽脅迫應(yīng)激反應(yīng)相關(guān)的候選miRNA,并且同一家族的不同miRNA 表達(dá)模式有所不同。差異表達(dá)的miRNA 的靶基因顯著富集的通路不同。這些miRNA 的鑒定及其功能意義的闡明拓寬了我們對(duì)鹽脅迫反應(yīng)中轉(zhuǎn)錄后基因調(diào)控的理解,為探索油菜耐鹽機(jī)制提供重要的參考。

    猜你喜歡
    文庫(kù)擬南芥抗性
    擬南芥:活得粗糙,才讓我有了上太空的資格
    一個(gè)控制超強(qiáng)電離輻射抗性開關(guān)基因的研究進(jìn)展
    專家文庫(kù)
    優(yōu)秀傳統(tǒng)文化啟蒙文庫(kù)
    幽默大師(2020年10期)2020-11-10 09:07:22
    關(guān)于推薦《當(dāng)代詩(shī)壇百家文庫(kù)》入選詩(shī)家的啟事
    尿黑酸對(duì)擬南芥酪氨酸降解缺陷突變體sscd1的影響
    專家文庫(kù)
    甲基對(duì)硫磷抗性菌的篩選及特性研究
    兩種LED光源作為擬南芥生長(zhǎng)光源的應(yīng)用探究
    擬南芥干旱敏感突變體篩選及其干旱脅迫響應(yīng)機(jī)制探究
    av网站免费在线观看视频| 欧美日本视频| 亚洲国产高清在线一区二区三 | 国产野战对白在线观看| 视频区欧美日本亚洲| 十八禁人妻一区二区| 亚洲第一欧美日韩一区二区三区| 黑人操中国人逼视频| ponron亚洲| 欧美激情 高清一区二区三区| 日本免费a在线| 99久久国产精品久久久| 久久久久久久久中文| 国产精品亚洲一级av第二区| 午夜免费观看网址| 午夜精品在线福利| 日日爽夜夜爽网站| 成人国产一区最新在线观看| 视频区欧美日本亚洲| 国产熟女午夜一区二区三区| 黄色毛片三级朝国网站| 黑人巨大精品欧美一区二区蜜桃| 久久中文字幕人妻熟女| 成年女人毛片免费观看观看9| 亚洲午夜理论影院| 亚洲国产毛片av蜜桃av| 搡老妇女老女人老熟妇| 午夜福利成人在线免费观看| 777久久人妻少妇嫩草av网站| 97人妻精品一区二区三区麻豆 | 久久久精品欧美日韩精品| 99精品久久久久人妻精品| 精品第一国产精品| 每晚都被弄得嗷嗷叫到高潮| 国产精品日韩av在线免费观看 | 亚洲国产毛片av蜜桃av| 亚洲精华国产精华精| 成熟少妇高潮喷水视频| 亚洲欧美日韩另类电影网站| 免费av毛片视频| 免费搜索国产男女视频| 日本三级黄在线观看| 国产精品国产高清国产av| 亚洲第一av免费看| 国产一区二区在线av高清观看| 欧美日本中文国产一区发布| 国产精品1区2区在线观看.| 大型黄色视频在线免费观看| 免费在线观看完整版高清| 婷婷六月久久综合丁香| 黄色视频,在线免费观看| 嫁个100分男人电影在线观看| 久久精品人人爽人人爽视色| 黄片小视频在线播放| 免费看美女性在线毛片视频| 99riav亚洲国产免费| 女人精品久久久久毛片| 国产aⅴ精品一区二区三区波| 免费观看精品视频网站| 中文字幕高清在线视频| 午夜福利影视在线免费观看| 变态另类成人亚洲欧美熟女 | 极品教师在线免费播放| 久久精品成人免费网站| 国产一卡二卡三卡精品| 成人免费观看视频高清| 午夜福利在线观看吧| 大码成人一级视频| 国产亚洲精品av在线| 999久久久精品免费观看国产| 亚洲国产毛片av蜜桃av| 成人免费观看视频高清| 国产欧美日韩一区二区精品| 精品久久久久久成人av| 精品国产国语对白av| 国产高清视频在线播放一区| 他把我摸到了高潮在线观看| 亚洲精品粉嫩美女一区| 国产精品久久久久久人妻精品电影| 男人的好看免费观看在线视频 | 韩国精品一区二区三区| 久久亚洲精品不卡| 成人免费观看视频高清| 99在线视频只有这里精品首页| 女警被强在线播放| 两人在一起打扑克的视频| 精品国内亚洲2022精品成人| 国产精品99久久99久久久不卡| 国产精品久久久av美女十八| 少妇被粗大的猛进出69影院| 成在线人永久免费视频| 久久这里只有精品19| 亚洲avbb在线观看| 久久青草综合色| 国产视频一区二区在线看| 好男人电影高清在线观看| 国产伦人伦偷精品视频| 亚洲七黄色美女视频| 三级毛片av免费| 最新美女视频免费是黄的| 搡老熟女国产l中国老女人| 亚洲欧美激情在线| 男人舔女人的私密视频| 韩国精品一区二区三区| 久久精品国产亚洲av香蕉五月| 国产精品爽爽va在线观看网站 | 国产精品久久久人人做人人爽| 欧美激情高清一区二区三区| 久9热在线精品视频| 亚洲国产精品久久男人天堂| 午夜免费鲁丝| 成年版毛片免费区| 成人三级黄色视频| АⅤ资源中文在线天堂| 欧美日韩乱码在线| 国产亚洲av嫩草精品影院| 国产成人啪精品午夜网站| 91精品国产国语对白视频| 国产亚洲av高清不卡| 亚洲成人精品中文字幕电影| 国产一区二区激情短视频| 老鸭窝网址在线观看| 久久久久国内视频| 精品欧美国产一区二区三| 亚洲成国产人片在线观看| 丰满人妻熟妇乱又伦精品不卡| 日日爽夜夜爽网站| 又紧又爽又黄一区二区| 韩国精品一区二区三区| 久久精品国产亚洲av香蕉五月| 国产欧美日韩精品亚洲av| 免费高清在线观看日韩| 午夜免费激情av| 中文字幕色久视频| 99久久99久久久精品蜜桃| 午夜视频精品福利| 久久精品aⅴ一区二区三区四区| 国内毛片毛片毛片毛片毛片| 久久午夜亚洲精品久久| 老司机福利观看| 成人三级做爰电影| 国产成人av激情在线播放| 午夜免费鲁丝| 亚洲自偷自拍图片 自拍| 欧美精品亚洲一区二区| 久久国产精品男人的天堂亚洲| 国产亚洲欧美在线一区二区| 可以在线观看的亚洲视频| 精品久久蜜臀av无| 1024视频免费在线观看| 最新在线观看一区二区三区| 久久人人爽av亚洲精品天堂| 久久精品91蜜桃| 国产一区在线观看成人免费| 成人免费观看视频高清| 给我免费播放毛片高清在线观看| 欧美成人免费av一区二区三区| 在线国产一区二区在线| ponron亚洲| 91精品三级在线观看| 男女之事视频高清在线观看| 香蕉国产在线看| 成熟少妇高潮喷水视频| 看片在线看免费视频| 成人特级黄色片久久久久久久| 村上凉子中文字幕在线| 欧美日韩中文字幕国产精品一区二区三区 | 动漫黄色视频在线观看| 日本欧美视频一区| 午夜福利免费观看在线| 成年人黄色毛片网站| 国产精品永久免费网站| 亚洲国产精品合色在线| 久久欧美精品欧美久久欧美| 国产三级黄色录像| 1024视频免费在线观看| aaaaa片日本免费| 久久青草综合色| 精品一区二区三区视频在线观看免费| 97人妻精品一区二区三区麻豆 | 国产精品久久久av美女十八| 女生性感内裤真人,穿戴方法视频| 91国产中文字幕| 亚洲成国产人片在线观看| 香蕉久久夜色| 久久久久亚洲av毛片大全| 国产精品二区激情视频| 免费人成视频x8x8入口观看| 国产区一区二久久| 亚洲色图综合在线观看| 日韩精品青青久久久久久| 亚洲精品国产精品久久久不卡| 亚洲五月婷婷丁香| 法律面前人人平等表现在哪些方面| 又紧又爽又黄一区二区| 男女做爰动态图高潮gif福利片 | 丝袜人妻中文字幕| 色综合欧美亚洲国产小说| 国产亚洲精品综合一区在线观看 | 欧美黄色淫秽网站| 国产成人精品久久二区二区免费| 韩国精品一区二区三区| 深夜精品福利| 日韩高清综合在线| 久久伊人香网站| 日韩精品中文字幕看吧| 人人妻人人澡欧美一区二区 | 法律面前人人平等表现在哪些方面| 精品高清国产在线一区| 伊人久久大香线蕉亚洲五| 日韩欧美三级三区| 国产精品久久久久久人妻精品电影| 18禁黄网站禁片午夜丰满| 午夜成年电影在线免费观看| 激情在线观看视频在线高清| 午夜影院日韩av| 国产精品久久久久久精品电影 | 日韩视频一区二区在线观看| 精品欧美国产一区二区三| 九色国产91popny在线| 成人特级黄色片久久久久久久| 久久精品国产亚洲av香蕉五月| a级毛片在线看网站| 亚洲国产中文字幕在线视频| 亚洲专区字幕在线| 国产高清视频在线播放一区| 国产一区二区在线av高清观看| 99香蕉大伊视频| 色播在线永久视频| 久久久久亚洲av毛片大全| 国产成人精品在线电影| 首页视频小说图片口味搜索| av电影中文网址| 亚洲成人精品中文字幕电影| 黑丝袜美女国产一区| 免费在线观看黄色视频的| 国产精品一区二区免费欧美| 日韩免费av在线播放| 日日干狠狠操夜夜爽| 午夜久久久在线观看| 亚洲专区字幕在线| av视频免费观看在线观看| 18禁裸乳无遮挡免费网站照片 | 精品无人区乱码1区二区| 啦啦啦韩国在线观看视频| 久久人人97超碰香蕉20202| av在线天堂中文字幕| 亚洲第一欧美日韩一区二区三区| 美国免费a级毛片| av在线播放免费不卡| 久久性视频一级片| 欧美久久黑人一区二区| 一级黄色大片毛片| 久久婷婷人人爽人人干人人爱 | 欧美国产精品va在线观看不卡| 人人妻,人人澡人人爽秒播| 国产一区二区在线av高清观看| 亚洲人成伊人成综合网2020| 日日夜夜操网爽| 老汉色∧v一级毛片| 亚洲精品美女久久久久99蜜臀| 高清毛片免费观看视频网站| 亚洲色图综合在线观看| 男人舔女人的私密视频| 99国产综合亚洲精品| 日本在线视频免费播放| 亚洲熟女毛片儿| 两个人免费观看高清视频| 视频在线观看一区二区三区| 国产精品久久电影中文字幕| 欧美乱色亚洲激情| 日本欧美视频一区| 午夜福利在线观看吧| 精品一区二区三区av网在线观看| 免费在线观看影片大全网站| 国产成人精品在线电影| 美女扒开内裤让男人捅视频| 久久久久国产精品人妻aⅴ院| 少妇被粗大的猛进出69影院| 嫁个100分男人电影在线观看| 长腿黑丝高跟| 母亲3免费完整高清在线观看| 国内久久婷婷六月综合欲色啪| 免费在线观看影片大全网站| 女性生殖器流出的白浆| 日本撒尿小便嘘嘘汇集6| 美女国产高潮福利片在线看| 丰满的人妻完整版| 欧美乱妇无乱码| 黄色毛片三级朝国网站| 最好的美女福利视频网| 久久人妻福利社区极品人妻图片| 久久草成人影院| 亚洲专区国产一区二区| 日韩欧美一区视频在线观看| 丝袜人妻中文字幕| 免费在线观看日本一区| 婷婷精品国产亚洲av在线| 亚洲欧美一区二区三区黑人| 国产成人欧美在线观看| 国产精品98久久久久久宅男小说| 免费av毛片视频| 日本 欧美在线| 50天的宝宝边吃奶边哭怎么回事| 国产av在哪里看| 国产精品 国内视频| 在线国产一区二区在线| 国产又色又爽无遮挡免费看| av超薄肉色丝袜交足视频| 成人18禁高潮啪啪吃奶动态图| 999久久久精品免费观看国产| 国产一区二区三区在线臀色熟女| 大码成人一级视频| 亚洲中文字幕一区二区三区有码在线看 | 久久人人97超碰香蕉20202| 黑丝袜美女国产一区| 国产片内射在线| 午夜两性在线视频| 天堂√8在线中文| 国产区一区二久久| 亚洲人成网站在线播放欧美日韩| 国产亚洲欧美在线一区二区| 性欧美人与动物交配| 国产免费男女视频| 99re在线观看精品视频| 啪啪无遮挡十八禁网站| 亚洲精品国产色婷婷电影| 久久久久久人人人人人| 国产精品自产拍在线观看55亚洲| 在线播放国产精品三级| 色在线成人网| 一区二区三区精品91| 9191精品国产免费久久| 国产私拍福利视频在线观看| 欧美人与性动交α欧美精品济南到| 国产亚洲精品一区二区www| 亚洲美女黄片视频| 亚洲精华国产精华精| 国产亚洲精品综合一区在线观看 | 男男h啪啪无遮挡| 久久天堂一区二区三区四区| 国产麻豆成人av免费视频| 曰老女人黄片| 十八禁网站免费在线| 亚洲欧美激情综合另类| 亚洲欧美精品综合一区二区三区| 亚洲男人天堂网一区| 女人精品久久久久毛片| 中文字幕高清在线视频| 国产主播在线观看一区二区| 午夜精品在线福利| 天天一区二区日本电影三级 | 级片在线观看| 麻豆久久精品国产亚洲av| 香蕉丝袜av| 久久青草综合色| 国产私拍福利视频在线观看| 亚洲人成电影免费在线| 国产亚洲精品一区二区www| av在线播放免费不卡| 欧美人与性动交α欧美精品济南到| 欧美精品啪啪一区二区三区| 国产91精品成人一区二区三区| 黄色片一级片一级黄色片| 欧美不卡视频在线免费观看 | 亚洲一区二区三区不卡视频| 国产私拍福利视频在线观看| 日本一区二区免费在线视频| 欧美激情久久久久久爽电影 | 亚洲欧美精品综合久久99| 久久香蕉精品热| 欧美激情高清一区二区三区| 欧美大码av| 成人国产一区最新在线观看| 99riav亚洲国产免费| 两性午夜刺激爽爽歪歪视频在线观看 | 久久天躁狠狠躁夜夜2o2o| 精品高清国产在线一区| 午夜a级毛片| 99香蕉大伊视频| 精品不卡国产一区二区三区| 亚洲九九香蕉| 国产精品久久电影中文字幕| 韩国精品一区二区三区| 禁无遮挡网站| 咕卡用的链子| 亚洲av成人av| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品中文字幕在线视频| 日本免费a在线| 亚洲国产精品久久男人天堂| 丁香欧美五月| 美女 人体艺术 gogo| 黄网站色视频无遮挡免费观看| 19禁男女啪啪无遮挡网站| 日本在线视频免费播放| 国产精品免费一区二区三区在线| 欧美一级a爱片免费观看看 | 免费av毛片视频| xxx96com| 美女高潮到喷水免费观看| 中文字幕高清在线视频| 国产亚洲精品久久久久5区| 亚洲最大成人中文| 久久精品成人免费网站| 老司机靠b影院| 日本一区二区免费在线视频| 午夜福利18| 欧美成人免费av一区二区三区| 日本 av在线| 在线播放国产精品三级| 一二三四社区在线视频社区8| 后天国语完整版免费观看| 亚洲欧美激情在线| 亚洲人成电影免费在线| 色在线成人网| 国产精品 欧美亚洲| 18禁国产床啪视频网站| 久久青草综合色| 91成年电影在线观看| 国产精品亚洲美女久久久| 最近最新中文字幕大全免费视频| 变态另类成人亚洲欧美熟女 | 久久人人精品亚洲av| 色综合欧美亚洲国产小说| 两个人免费观看高清视频| 久久热在线av| 久久久久久大精品| 亚洲va日本ⅴa欧美va伊人久久| cao死你这个sao货| 亚洲熟女毛片儿| 久久香蕉激情| 男女午夜视频在线观看| 日韩欧美在线二视频| 精品乱码久久久久久99久播| 一区二区三区国产精品乱码| 日日干狠狠操夜夜爽| 久久国产精品影院| 欧美色视频一区免费| 久久欧美精品欧美久久欧美| 成人亚洲精品av一区二区| 亚洲激情在线av| 国产99久久九九免费精品| 露出奶头的视频| 欧美乱色亚洲激情| 这个男人来自地球电影免费观看| 老熟妇仑乱视频hdxx| 亚洲免费av在线视频| 日韩精品免费视频一区二区三区| 久久久久久久久免费视频了| videosex国产| 天天一区二区日本电影三级 | 嫩草影视91久久| 少妇的丰满在线观看| 国产不卡一卡二| 正在播放国产对白刺激| 亚洲精品国产精品久久久不卡| 香蕉久久夜色| 国产精品久久视频播放| 纯流量卡能插随身wifi吗| 欧美丝袜亚洲另类 | 久久久久精品国产欧美久久久| 日本精品一区二区三区蜜桃| 欧美乱码精品一区二区三区| 一个人观看的视频www高清免费观看 | 亚洲av美国av| 久久精品国产综合久久久| cao死你这个sao货| 老汉色∧v一级毛片| 欧美色欧美亚洲另类二区 | 一级毛片女人18水好多| 天天添夜夜摸| 在线观看66精品国产| 一级作爱视频免费观看| 美女免费视频网站| 黄片播放在线免费| av在线天堂中文字幕| 777久久人妻少妇嫩草av网站| 国产av又大| 久久久久亚洲av毛片大全| 亚洲avbb在线观看| 欧美日本中文国产一区发布| 两性午夜刺激爽爽歪歪视频在线观看 | av视频免费观看在线观看| 脱女人内裤的视频| 欧美激情高清一区二区三区| 午夜两性在线视频| 欧美色欧美亚洲另类二区 | 欧美在线黄色| 9色porny在线观看| 丝袜美腿诱惑在线| 又紧又爽又黄一区二区| 精品国产一区二区久久| 男人舔女人下体高潮全视频| 脱女人内裤的视频| 满18在线观看网站| 亚洲人成网站在线播放欧美日韩| 精品一区二区三区四区五区乱码| 精品人妻在线不人妻| 成人欧美大片| 久久久久久久精品吃奶| 一区二区三区高清视频在线| 久久精品国产99精品国产亚洲性色 | 色精品久久人妻99蜜桃| 久久人人97超碰香蕉20202| 色av中文字幕| 日韩三级视频一区二区三区| 又黄又爽又免费观看的视频| 真人一进一出gif抽搐免费| 波多野结衣av一区二区av| 国产aⅴ精品一区二区三区波| 一个人免费在线观看的高清视频| 国产精品免费视频内射| 精品无人区乱码1区二区| 日韩欧美一区视频在线观看| 真人做人爱边吃奶动态| 一个人观看的视频www高清免费观看 | 亚洲视频免费观看视频| 丝袜美足系列| 亚洲精品久久国产高清桃花| 成人三级黄色视频| 精品一区二区三区av网在线观看| 欧美日韩乱码在线| 国产一级毛片七仙女欲春2 | 国产精品国产高清国产av| 午夜福利在线观看吧| 一区二区三区激情视频| 国产色视频综合| 精品国产乱子伦一区二区三区| 成熟少妇高潮喷水视频| 久久中文看片网| 看片在线看免费视频| 午夜亚洲福利在线播放| 成年版毛片免费区| 亚洲五月婷婷丁香| 999久久久国产精品视频| 精品无人区乱码1区二区| 久久久久亚洲av毛片大全| 国产成人一区二区三区免费视频网站| 91字幕亚洲| 亚洲成a人片在线一区二区| 午夜福利18| 久久久久国产一级毛片高清牌| 巨乳人妻的诱惑在线观看| 久久热在线av| or卡值多少钱| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品成人综合色| 一边摸一边做爽爽视频免费| 成人手机av| 精品国产超薄肉色丝袜足j| 两性夫妻黄色片| 91字幕亚洲| 人妻久久中文字幕网| 在线观看一区二区三区| 亚洲av电影不卡..在线观看| 国产精品久久久久久精品电影 | 后天国语完整版免费观看| 午夜精品国产一区二区电影| 国产欧美日韩综合在线一区二区| 国产精品爽爽va在线观看网站 | 两性夫妻黄色片| 精品少妇一区二区三区视频日本电影| 亚洲人成电影观看| 久久人人97超碰香蕉20202| 老司机午夜福利在线观看视频| 久久精品aⅴ一区二区三区四区| 久久久久久久久免费视频了| 久久精品人人爽人人爽视色| 欧美激情久久久久久爽电影 | 不卡一级毛片| 三级毛片av免费| 一a级毛片在线观看| 黄色丝袜av网址大全| 在线免费观看的www视频| 91精品国产国语对白视频| 搡老妇女老女人老熟妇| 麻豆成人av在线观看| 欧美日韩一级在线毛片| 久久国产精品男人的天堂亚洲| 九色国产91popny在线| 国产精品日韩av在线免费观看 | 亚洲性夜色夜夜综合| 国产精品亚洲美女久久久| 欧美成人一区二区免费高清观看 | 99国产精品免费福利视频| 看黄色毛片网站| 免费高清视频大片| 99国产精品一区二区三区| 欧美成狂野欧美在线观看| 大型黄色视频在线免费观看| 两人在一起打扑克的视频| 黄色女人牲交| 高潮久久久久久久久久久不卡| 久久久久久久午夜电影| 欧美性长视频在线观看| 给我免费播放毛片高清在线观看| 午夜影院日韩av| 多毛熟女@视频| 深夜精品福利| 亚洲午夜理论影院| 人成视频在线观看免费观看| 国产精品亚洲一级av第二区| 亚洲人成电影免费在线| 国产精品久久久久久亚洲av鲁大| 女人被狂操c到高潮| 他把我摸到了高潮在线观看| 亚洲精品av麻豆狂野| 国产亚洲精品久久久久5区| 国产精品亚洲一级av第二区| 无人区码免费观看不卡| 国产成人系列免费观看| 久久影院123| 91成年电影在线观看| 757午夜福利合集在线观看|