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

    重金屬脅迫下白骨壤數(shù)字基因表達譜分析

    2017-04-19 08:44:55袁長春莫健新袁柳嬌余如鳳鐘軍弟劉鍇棟
    林業(yè)科學研究 2017年2期
    關鍵詞:白骨紅樹測序

    袁長春,莫健新,袁柳嬌,余如鳳,劉 倩,鐘軍弟,劉鍇棟

    (1.嶺南師范學院生命科學與技術學院,廣東 湛江 524048; 2.廣州基迪奧生物科技有限公司,廣東 廣州 510006)

    重金屬脅迫下白骨壤數(shù)字基因表達譜分析

    袁長春1*,莫健新1,袁柳嬌1,余如鳳1,劉 倩2,鐘軍弟1,劉鍇棟1*

    (1.嶺南師范學院生命科學與技術學院,廣東 湛江 524048; 2.廣州基迪奧生物科技有限公司,廣東 廣州 510006)

    [目的]從轉錄組水平研究重金屬污水脅迫處理下白骨壤的基因表達變化,以揭示白骨壤響應重金屬污染脅迫的機制。[方法]采用基于高通量測序的數(shù)字基因表達譜(DGE)技術對重金屬污水脅迫處理組和對照組白骨壤材料進行轉錄組測序分析。[結果]與對照組相比,重金屬處理下白骨壤共有10 459個差異表達基因,其中,8 685個表達量上升,1 774個表達量下降。MeV聚類表明:大部分基因在重金屬處理樣本中的表達量受到極大的誘導。GO功能注釋發(fā)現(xiàn),差異表達基因主要定位于葉綠體以及細胞內各種膜結構,參與“細胞代謝”、“細胞組分的組合和生物合成”、“對刺激的響應”等過程。KEGG 代謝通路分析顯示,差異表達基因分布于126條Pathways,涉及“核糖體”、“光合作用”、“乙醛酸鹽代謝”、“丙酮酸代謝”等途徑。重金屬脅迫還促進萜類、黃酮類化合物生物合成的關鍵酶基因(牻牛兒基焦磷酸合酶(0009238)、黃酮醇合成酶(0001833))的表達,進而促進白骨壤有效成分的積累;顯著誘導細胞分裂素水解酶編碼基因CKX7(0057124)、油菜素內脂信號轉導組分基因BSK(0043741)等的表達,進而提高白骨壤對重金屬脅迫的適應能力。此外,轉錄因子分析發(fā)現(xiàn),bHLH、NAC、MYB-related和WRKY在重金屬脅迫中發(fā)揮重要作用。實驗選取8個與環(huán)境刺激響應密切相關的基因,通過qRT-PCR驗證了它們在重金屬污水脅迫處理下的表達差異,結果與數(shù)字基因表達譜分析的結果較一致,證實了差異表達基因數(shù)據(jù)的有效性。[結論]重金屬處理影響了白骨壤大量的新陳代謝、光合作用、小分子酸合成、激素合成與信號轉導及次生代謝產(chǎn)物合成相關基因的表達。

    白骨壤;重金屬脅迫;數(shù)字基因表達譜(DEG);差異表達基因;機制

    隨著我國沿海地區(qū)工農(nóng)業(yè)的快速發(fā)展以及城市化區(qū)域的不斷擴大,工農(nóng)廢水和生活污水的排放量越來越大。大量污水涌入江河海灣造成了很大的環(huán)境壓力[1]。紅樹濕地系統(tǒng)處在陸地向海洋的過度帶,對各種污染物有很強的吸附能力,具有凈化污水、沉積污染物的能力[2-4];然而,污水排放對紅樹植物的正反影響效應存在較大的爭議。一方面,污水排放能為紅樹植物帶來豐富的營養(yǎng)物質,促進紅樹的生長;另一方面,污水中重金屬污染物的脅迫作用對紅樹植物存在一定的不利影響,減緩紅樹植物的生長發(fā)育。因此,探究紅樹植物響應污染脅迫的機制是當前研究的熱點[5]。有研究表明,紅樹植物能夠通過多種途徑應對重金屬脅迫,如:根部發(fā)達的凱氏帶可以防止重金屬離子向地上部分轉運;植物體內的小分子有機酸等物質可以將重金屬沉淀于細胞壁或局限于液泡等細胞器;通過增強植物體自身抗氧化系統(tǒng)活性或促進各種抗氧化物質合成來清除重金屬脅迫下產(chǎn)生的大量活性氧自由基[6-8];然而,這些研究大多停留在生理層面,從分子層面分析重金屬脅迫處理對紅樹植物生長影響的研究鮮有報道[9-10]。

    數(shù)字基因表達譜(DGE)是一種高效的高通量測序技術,DEG能用于分析某一物種特定組織在特定條件下基因的表達情況。通過對差異基因涉及到的生物學過程進行搜索、比較和分析,可預測基因的表達模式和蛋白質的功能等信息。目前,已有學者利用該技術在轉錄水平分析了美洲黑楊(PopulusdeltoidesMarsh)、泡泡刺(NitrariasphaerocarpaMaxim.)、花生(ArachishypogaeaLinn.)等植物基因表達譜[11-13]。如,賴曉娟等[14]通過分析壇紫菜(PyropiahaitanensisRhodophyta)的數(shù)字基因表達譜,篩選出了256個上調表達基因和3 820個下調表達基因,證明了高溫脅迫對壇紫菜的生長發(fā)育有一定的抑制作用,但同時也誘導了補救途徑。

    作為常見的紅樹林優(yōu)勢樹種白骨壤(Avicenniamarina(Forssk?l) Vierhapper)廣泛分布于陸地向海洋過渡的潮間帶,在我國分布于海南、廣東、廣西等省的濱海地區(qū)。白骨壤能耐受一定濃度的污染脅迫,可在污水污染的潮間帶生長[15]。有研究表明,白骨壤模擬濕地系統(tǒng)可對污水起到很好的凈化效應,在栽種白骨壤植物的土壤子系統(tǒng)中,重金屬積累含量顯著降低[16]。進一步研究表明,白骨壤正是通過調控濕地系統(tǒng)中氮、磷的分配、循環(huán)達到凈化作用[17-18]。周涵韜等[19]通過mRNA差異顯示技術,分離得到了白骨壤耐鹽相關cDNA,為揭示白骨壤鹽脅迫響應機制打下了基礎。本文采用基于高通量測序的數(shù)字基因表達譜技術研究污水脅迫處理下白骨壤的基因表達變化,從分子層面分析白骨壤對重金屬污染脅迫的耐受機制,為深入研究重金屬污水脅迫對紅樹植物生長的影響機制提供參考。

    1 材料與方法

    1.1 實驗試劑與儀器

    Hoagland營養(yǎng)液;人工污水;Trizol試劑盒(TaKaRa, 大連,中國);瓊脂糖;寡聚(dT)-纖維素;片段化緩沖液(fragmentation buffer);dNTPs;六堿基隨機引物;RNase H;DNA polymerase I;QiaQuick PCR試劑盒;EB緩沖液;PrimeScriptTMRT Master Mix (TaKaRa,大連,中國);SYBR Green supermix;PCR儀型:Bio-Rad CFX96 touchTM。

    1.2 白骨壤材料處理和收集

    1.3 總RNA提取和基因表達譜測序

    采用Trizol法提取污水脅迫處理組和對照組白骨壤葉片樣品的總RNA,對照和處理各3份生物學重復;然后,將對照和處理的各自3份RNA樣品等量混合,由廣州基迪奧公司構建測序文庫,采用Illumina HiSeqTM 2000平臺進行轉錄組測序。

    1.4 差異表達基因的篩選

    基因的表達量由比對到基因表達區(qū)域序列的readcount頻率來估算,通過edgeR軟件對得到的readcount數(shù)據(jù)進行差異分析,以篩選差異表達基因。本實驗差異基因篩選標準為P-value<0.05,|log2Fold Change|>1,其中,F(xiàn)old Change 為污水脅迫處理樣本與對照樣本中各基因表達量的比值。對于差異表達基因,其log2Fold Change>0時,認為該差異表達基因是上調的;反之,若log2Fold Change <0,則認為該差異表達基因是下調的。

    1.5 Gene Ontology (GO) 基因功能注釋和KEGG代謝通路分析

    采用GO數(shù)據(jù)庫(http://www.geneontology.org/)對污水脅迫處理下白骨壤的差異表達基因進行GO功能富集顯著性分析。采用GOseq方法計算得到一個特定的P-value。一般校正后的P-value≤0.05,表示差異表達基因在該GO中出現(xiàn)了富集。同時,使用KOBAS(2.0)軟件對差異基因進行KEGG 代謝通路富集性分析,當FDR(錯誤發(fā)現(xiàn)率)≤0.05時,表示該通路中差異表達基因富集顯著。

    1.6 qRT-PCR驗證

    表1 引物序列

    2 結果與分析

    2.1 測序質量分析

    分析白骨壤數(shù)字基因表達譜測序數(shù)據(jù)發(fā)現(xiàn),重金屬污水脅迫處理樣品和對照樣品測序所得數(shù)據(jù),經(jīng)過濾篩除含有接頭和低質量部分后,得到高質量Clean reads的比例分別為98.26%、98.23%,證明測序質量較高。將得到的Clean reads與參考基因比對,對照樣本和重金屬污水脅迫處理樣本分別有89.20%和87.45%的Clean reads可比對上參考基因,這為白骨壤轉錄組序列的注釋打下了結實的基礎。測序結果提交到Genebank,登錄號為SRA487295。

    飽和度是用來衡量單一樣本測序得到的基因數(shù)量是否覆蓋全基因組的標準,檢測所得的基因數(shù)量隨測序量的增加而增加,當測序量接近基因組的大小時,其檢測到的基因總數(shù)達到峰值,表明檢測到的基因數(shù)已基本覆蓋全基因組。本研究中,對照樣本和污水脅迫處理樣本的Reads數(shù)量超過10 M時,檢測到的基因數(shù)量趨于飽和,而實驗中2個樣品的測序量分別為25.26、25.64 M,證明絕大部分基因已被檢測到。隨機性分析表明:各樣本的Reads呈正態(tài)分布,具有很高的隨機性。上述結果表明:本次測序結果能夠較好地反映各樣本中基因表達的真實情況,可用于分析重金屬污水脅迫處理下白骨壤基因的表達變化。

    2.2 差異表達基因篩選

    MeV聚類分析顯示:聚在同一類的基因表達量在對照樣本和重金屬污水脅迫處理樣本間差異極大,大部分基因表現(xiàn)為在污水脅迫處理樣本中的表達量比對照樣本中的高。

    通過比較對照樣品與重金屬污水脅迫處理樣品基因表達譜,共篩選出10 459個差異表達基因,其中,有8 685個基因表達量增加,1 774個基因表達量減少,表明重金屬污水脅迫處理極大的誘導了白骨壤基因的表達。

    2.3 差異表達基因GO功能顯著性富集分析

    將白骨壤差異表達基因比對到GO數(shù)據(jù)庫,共得到2 041條terms,其中,“分子功能”514條,占25.18%;“細胞組分”202條,占9.90%;“生物學過程” 1 325條,占64.92%。以P-value≤0.01為標準篩選極顯著富集條目,共得到230條,其中,“分子功能”39條,“細胞組分”46條,“生物學過程”145條。對差異表達基因GO功能分類發(fā)現(xiàn),大部分類別中,表達量上調的基因數(shù)比表達量下調的基因數(shù)多,其中,與“生物學過程”相關的差異表達基因主要涉及“代謝過程”、“對刺激因素的響應”、“區(qū)域化”、“生物調節(jié)”等過程,與“分子功能”相關的差異表達基因主要與“結構分子活性”、“轉運體活性”、“催化活性”、“結合活性”等分子活性相關,與“細胞組分”相關的差異表達基因主要與“細胞器組分”、“細胞組分”、“大分子復合物”、“膜組分”等結構相關。此外,重金屬污水脅迫處理還提高了白骨壤的“抗氧化因子”、“蛋白質結合轉錄因子”、“核酸結合轉錄因子”、“電子載體”、“鳥嘌呤核苷酸交換因子”等關鍵因子的活性,進而影響“生物粘附”、“ 信號傳導”、“ 多細胞調節(jié)”等過程。

    2.4 KEGG 代謝通路顯著性富集分析

    通過KEGG 代謝通路富集分析,可以揭示差異表達基因參與的主要生理活動途徑。結果表明:重金屬污水脅迫處理下,白骨壤差異表達基因輻射到126條Pathway中。從表2可以看出:在P≤0.05顯著性富集的20條Pathway中,富集在“核糖體途徑”的基因數(shù)量最多,共有409個差異表達基因;其次是“碳代謝途徑”,富集了303個差異表達基因。“乙醛酸鹽代謝”途徑、“光合生物固碳途徑”和“丙酮酸代謝途徑”也是3個顯著性富集的代謝途徑,分別富集了107、103、113個差異表達基因,占差異表達基因總數(shù)的6.32%、6.08%、6.67%。此外,在“苯丙環(huán)素的生物合成”、“氮素代謝”、“光合作用”、“脂肪酸代謝”、“氨基酰-tRNA的生物合成”、“半乳糖代謝”等途徑也富集了一定數(shù)量的差異表達基因。

    表2 KEGG代謝通路富集分析

    2.5 污水脅迫處理對白骨壤激素合成、運輸和信號途徑基因的表達調控

    植物激素是一類參與植物生長發(fā)育、代謝和環(huán)境響應的重要化學物質[20]。通過對白骨壤DEG分析,得到了生長素、乙烯、細胞分裂素、脫落酸、赤霉素和油菜素內脂等一系列激素相關基因的表達差異情況,并分析了污水脅迫對上述基因的影響。結果(表3)表明:在污水脅迫處理下,一個生長素上調小RNA基因SAUR(0052024)顯著上調。多個乙烯信號途徑基因也得到了分析,其中,乙烯不敏感3基因EIN3(0016997)的表達保持穩(wěn)定,而乙烯響應轉錄因子1基因ERF1(0013995)的表達受到污水脅迫的明顯誘導。眾多細胞分裂素基因的表達受到污水脅迫的顯著影響,如細胞分裂素脫氫酶7基因CKX7(0057124)的表達顯著上調,而一個富含組氨酸的磷酸轉移蛋白1基因AHP1(0029198)的表達受到極大的抑制。另一個內源細胞分裂素相關基因玉米黃質環(huán)氧化酶基因ZEP(0021189)在污水脅迫處理下表達水平也顯著上調。在脫落酸途徑中,一個脫落酸不敏感5基因ABI5(0061511)的表達顯著上調。最后筆者發(fā)現(xiàn),油菜素內脂的一個信號轉導組分基因——油菜素內酯信號激酶基因BSK(0043741)的表達受到極大的誘導。

    表3 激素相關基因在污水脅迫下的表達變化

    2.6 污水脅迫處理影響白骨壤刺激代謝關鍵酶編碼基因的表達

    通過KEGG分析發(fā)現(xiàn):污水處理可顯著影響白骨壤中萜類和黃酮類次生代謝相關基因的表達(表4),其中,萜類物質生物合成的限速酶牻牛兒基焦磷酸合酶(0009238)、甲基赤蘚醇磷酸胞苷酰轉移酶(0057458)和1-脫氧-D-木酮糖-5-磷酸還原異構酶(0025369)的表達都受到污水脅迫的顯著誘導。另外,黃酮類生物合成的幾個關鍵酶,如:黃酮醇合成酶(0001833)和查爾酮合成酶(0000917)的編碼基因表達水平也受到明顯的誘導。大量次生代謝相關基因的上調表達,暗示污水脅迫可能加速白骨壤內次生代謝產(chǎn)物的積累。

    表4 次生代謝相關基因在污水脅迫下的表達變化

    2.7 重金屬污水脅迫下的轉錄因子分析

    一系列轉錄因子家族成員在植物環(huán)境脅迫響應中起到了重要的作用[21]?;谀J街参飻M南芥(Arabidopsisthaliana(L.) Heynh.)轉錄因子數(shù)據(jù)庫,筆者對白骨壤轉錄組進行比對和注釋,一共鑒定到了分屬于56個常見家族的6 662個轉錄因子,轉錄因子家族及其包含的家族成員數(shù)據(jù)見表5。研究發(fā)現(xiàn):bHLH、NAC、MYB_related、WRKY等家族在白骨壤中擁有大量的成員。此外,一批與激素相關的轉錄因子家族也得到了鑒定,如ARF、ERF、ARR-B等;一些與植物環(huán)境脅迫密切相關的轉錄因子家族也得到了不同程度的分析,如bZIP、GATA家族。鑒定轉錄因子家族,有利于形成一個完整的調控網(wǎng)絡,揭示白骨壤污水脅迫響應的分子機制。

    表5 轉錄因子分析

    2.8 差異表達基因qRT-PCR驗證

    為了驗證測序數(shù)據(jù)的可靠性,筆者選擇顯著富集的GO分類條目“GO:0016491”、“GO:0005198”、“GO:0016866”和“GO:0003723”中的8個差異表達基因進行qRT-PCR分析,這8個基因分別是:0003922、0048831、0026527、0053069、0020253、0055711、0016422和0049777,它們都是白骨壤中參與基礎代謝的關鍵基因。實驗結果表明:這8個基因的qRT-PCR結果與數(shù)字基因表達譜中檢測到的表達趨勢有較高的一致性。

    3 討論

    重金屬污染是影響植物生長發(fā)育較為嚴重的逆境。植物在長期進化過程中,響應重金屬脅迫而形成了特定的生長習性和形態(tài)、生理特征。作為一種有代表性的紅樹植物,白骨壤對污染脅迫表現(xiàn)出極強的耐受性[9],展開對白骨壤的轉錄組分析,有利于從基因層面揭示紅樹植物耐受污染脅迫的機制。在本研究中,污水脅迫處理下,白骨壤的差異表達基因大部分上調表達,說明白骨壤可能通過提高自身轉錄水平應對污染脅迫,并增強適應能力。GO顯著性分析發(fā)現(xiàn),“代謝過程”共富集了1 792個上調表達基因,178個下調表達基因;KEGG代謝通路顯著性分析表明,“核糖體”、“氮素代謝”、“碳代謝”、“丙酮酸代謝”、“脂肪酸代謝”等途徑中富集了大量的差異表達基因(大部分為上調基因),說明白骨壤可能通過上調表達與基本代謝相關的基因來加快有害物質的新陳代謝,從而增強自身對污染脅迫的耐受能力。

    光合作用的強弱能反映植物體對外界環(huán)境變化的適應情況。微量的重金屬能刺激紅樹植物葉綠素的合成,使光合作用增強,以提供更多的能量促進植物體的生長[22];然而,紅樹植物對污染脅迫也有一定的耐受范圍,當脅迫濃度超過其耐受的上限時,會抑制葉片的光合作用。用大于30 mg·kg-1的Cd處理桐花樹(Aegicerascorniculatum(Linn.) Blanco)發(fā)現(xiàn),葉片細胞內膜結構被破壞,葉綠素含量降低[22];另有報道,在高濃度Hg2+(50 mg·L-1)處理下,白骨壤幼苗葉片凈光合作用的速率降低[23]。此外,黃國勇等[24]報道,增加復合重金屬的脅迫濃度,秋茄(Kandeliacandel(L.) Druce)葉片葉綠素的合成量下降加劇。原因可能是重金屬離子與葉綠素合成酶的功能結構域結合,阻礙它們發(fā)揮作用。另外,重金屬離子也能加快葉綠素的分解速率。在本研究中,KEGG 代謝通路分析發(fā)現(xiàn),“光合作用途徑”中富集了83個差異表達基因,說明在污染脅迫下白骨壤可能通過調節(jié)與光合作用相關基因的表達,以維持葉綠素的正常合成和葉綠體的結構完整性,從而使光合作用得以正常進行。

    大量研究表明,激素作為植物體內一類重要的化學成分參與植物的抗逆反應[25]。在白骨壤的DEGs(差異表達基因)中,鑒定得到一批與激素相關的基因,它們的表達水平受到污水脅迫處理的調控。在擬南芥(Arabidopsisthaliana(L.) Heynh.)和甘草(GlycyrrhizauralensisFisch.)中,鎘、銅或者汞脅迫顯著誘導了ERF家族基因的表達[26-27]。本研究中,一個乙烯響應轉錄因子ERF1(0013995)受到誘導,這與模式植物的表達情況一致。脫落酸是已知的逆境激素,張春榮等[27]發(fā)現(xiàn),干旱脅迫顯著誘導ABA生物合成關鍵酶ZEP基因和信號轉導組分ABI5基因的表達。在本研究中,細胞分裂素相關基因玉米黃質環(huán)氧化酶ZEP(0021189)和脫落酸不敏感基因ABI5(0061511)都受到重金屬污水脅迫的顯著誘導。重金屬脅迫顯著促進了白骨壤乙烯、細胞分裂素的生物合成和信號轉導,這有可能進而促進根的發(fā)生和細胞增殖,從而提高脅迫下的植株適應能力。

    環(huán)境脅迫往往會導致一些次生代謝產(chǎn)物的積累,研究白骨壤次生代謝相關基因的表達,有助于揭示次生代謝途徑參與白骨壤環(huán)境脅迫的分子機制[28]。在眾多的次生代謝成分中,黃酮類和萜類物質與環(huán)境脅迫的聯(lián)系最緊密。有報道稱,多種環(huán)境脅迫可以誘導黃酮類植物在植物體內的積累[29]。通過對差異表達基因的分析,筆者發(fā)現(xiàn),與黃酮和萜類化合物次生代謝相關的關鍵基因的表達都受到污水脅迫處理的一致上調,其中,黃酮合成的2個上游酶編碼基因:黃酮醇合成酶(0001833)和查爾酮合成酶(0000917)都受到明顯的誘導。這進一步說明,積累的次生代謝產(chǎn)物有助于白骨壤污水脅迫耐受性的形成。

    前期的研究表明,一系列轉錄因子參與了植物對環(huán)境脅迫的響應。作為植物中最大的家族之一, MYB類轉錄因子廣泛參與植物次生代謝調控、細胞形態(tài)發(fā)生、脅迫應答、分生組織形成及細胞周期控制等[30]。在白骨壤中,筆者一共鑒定得到219個MYB轉錄因子和431個MYB-related轉錄因子。NAC家族轉錄因子也被認為是一類重要的脅迫調控因子。OsMAC5和OsNAC6都是水稻(OryzasativaL.)響應生物及非生物脅迫的關鍵因子,它們與植物的脅迫耐受性密切相關[31-32]。在白骨壤中,可能的NAC轉錄因子一共有487個。以上結果為進一步分析轉錄因子家族在白骨壤污水脅迫中的響應機制奠定基礎。

    4 結論

    本實驗以重金屬處理下的白骨壤差異基因表達譜為基礎,得到了大量重金屬響應基因,并預測了這些基因可能的生物學功能。筆者分析表明,重金屬處理影響了白骨壤大量的新陳代謝、光合作用、小分子酸合成、激素合成與信號轉導及次生代謝產(chǎn)物合成相關基因的表達。白骨壤對重金屬脅迫響應機制具有極高的復雜性,深入研究一些關鍵基因將有助于全面了解其分子機制。

    [1] 丘耀文,余克服. 海南紅樹林濕地沉積物中重金屬的累積[J]. 熱帶海洋學報, 2011, 30(2): 102-108.

    [2] 萬 璐,康麗華,廖寶文,等.紅樹林根際解磷菌分離、培養(yǎng)及解磷能力的研究[J]. 林業(yè)科學研究,2004,17(1):89-94.

    [3] Li R L, Chai M W, Qiu G Y. Distribution, fraction, and ecological assessment of heavy metals in sediment-plant system in Mangrove forest, South China Sea[J]. PLoS One,2016, 11(1): e0147308.

    [4] 白保勛,樊 巍,卞新民,等.生活污水脅迫對Ⅰ-72楊樹人工林的影響[J].林業(yè)科學研究,2010,23(4):602-606.

    [5] 陳桂葵, 陳桂珠, 黃玉山, 等. 人工污水對白骨壤生長的影響[J]. 海洋環(huán)境科學, 2003, 22(3): 39-43.

    [6] 程 皓, 陳桂珠, 葉志鴻. 紅樹林重金屬污染生態(tài)學研究進展[J]. 生態(tài)學報, 2009, 29(7): 3893-3900.

    [7] MacFarlane G R, Burchett M D. Zinc distribution and excretion in the leaves of the grey mangrove,Avicenniamarina(Forsk) Vierh[J]. Environmental and Experimental Botany, 1999, 41: 167-175.

    [8] Zhang Z W, Xu X R, Sun Y X,etal. Heavy metal and organic contaminants in mangrove ecosystems of China: a review[J]. Environ Sci Pollut Res Int, 2014, 21(20):11938-11950.

    [9] 陳昌徐, 陸志強, 章耕耘, 等. 鎘與萘復合脅迫對紅樹植物白骨壤幼苗萌芽及生長的影響[J]. 熱帶亞熱帶植物學報, 2016, 24(1):80-86.

    [10] 陸俊錕,康麗華,陳 俊,等.華南三地紅樹林土壤微生物及其與土壤化學性質的相關性研究[J].林業(yè)科學研究,2008,21(4):523-527.

    [11] 丁昌俊, 張偉溪, 高 暝, 等. 不同生長勢美洲黑楊轉錄組差異分析[J]. 林業(yè)科學, 2016, 52(3):47-58.

    [12] 馬 婧, 鄧 楠, 褚建民, 等. 泡泡刺高通量轉錄組鑒定及其黃酮類代謝途徑初步分析[J]. 林業(yè)科學研究, 2016, 29(1):61-66.

    [13] 孫愛清, 張杰道, 萬勇善, 等. 花生干旱脅迫響應基因的數(shù)字表達譜分析[J]. 作物學報, 2013, 39(6): 1045-1053.

    [14] 賴曉娟, 嚴小軍, 楊 銳, 等.高溫脅迫下壇紫菜的數(shù)字基因表達譜研究[J]. 海洋學報, 2014, 36(6):104-111.

    [15] Wang L Y, Wang Y S, Zhang J P,etal. Molecular cloning of class III chitinase gene fromAvicenniamarinaand its expression analysis in response to cadmium and lead stress[J]. Ecotoxicology, 2015, 24:1697-1704.

    [16] 陳桂珠,陳桂葵,譚鳳儀,等. 白骨壤模擬濕地系統(tǒng)對污水的凈化效應[J]. 海洋環(huán)境科學, 2000, 19(4):23-26.

    [17] 陳桂葵,陳桂珠. 白骨壤模擬濕地系統(tǒng)中磷的分配循環(huán)及其凈化效應[J]. 生態(tài)學報, 2005, 25(3):627-632.

    [18] 陳桂葵,陳桂珠. 模擬分析白骨壤濕地系統(tǒng)中Ni的分配循環(huán)及其凈化效果[J]. 海洋環(huán)境科學, 2005, 24(4):16-19.

    [19] 周涵韜,林 鵬. 利用mRNA差別顯示技術分離鹽脅迫下紅樹植物白骨壤耐鹽相關cDNA[J]. 生物工程學報, 2002, 18(1):51-54.

    [20] 齊飛艷,彭鎮(zhèn)華,胡 陶,等.毛竹花期不同器官內源激素含量的變化[J].林業(yè)科學研究,2013,26(3):332-336.

    [21] 王計平, 史華平, 毛 雪, 等. 轉錄因子網(wǎng)絡與植物對環(huán)境脅迫的響應[J]. 應用生態(tài)學報, 2006, 17(9):1740-1746.

    [22] 吳桂容, 嚴重玲. 鎘對桐花樹幼苗生長及滲透調節(jié)的影響[J]. 生態(tài)環(huán)境, 2006, 15(5): 1003-1008.

    [23] 袁彥婷, 丁振華, 張 玲. 汞脅迫對白骨壤(Avicenniamarina)幼苗生理生態(tài)的影響[J]. 生態(tài)學雜志, 2011, 30(5): 1013-1017.

    [24] 黃國勇, 王友紹, 孫翠慈, 等. 秋茄葉對復合重金屬的脅迫反應及其積累能力研究[J]. 熱帶海洋學報, 2010, 29(6): 104-109.

    [25] Montero-Palmero M B, Martin-Barranco A, Escobar C,etal. Early transcriptional responses to mercury: a role for ethylene in mercury-induced stress[J]. New Phytologist, 2014, 201(1):116-130.

    [26] Schellingen K, Van Der Straeten D, Vandenbussche F,etal. Cadmium-induced ethylene production and responses inArabidopsisthalianarely onACS2 andACS6 gene expression[J]. BMC Plant Biology,2014,14:214.

    [27] 張春榮, 桑雪雨, 渠 萌, 等. 基于轉錄組測序揭示適度干旱脅迫對甘草根基因表達的調控[J]. 中國中藥雜志, 2015, 40(24):4817-4823.

    [28] 覃柳燕, 李 力, 繆劍華. 環(huán)境脅迫對藥用植物次生代謝產(chǎn)物的影響[J]. 湖北農(nóng)業(yè)科學, 2012, 51(8):1528-1532.

    [29] 梁 濱, 季平揚, 唐宋輝, 等. La對UV-B脅迫下大豆幼苗類黃酮含量影響:I對類黃酮總量影響[J]. 農(nóng)業(yè)環(huán)境科學學報, 2006, (3):584-586.

    [30] 喬 孟, 于延沖, 向鳳寧. 擬南芥R2R3-MYB類轉錄因子在環(huán)境脅迫中的作用[J]. 生命科學, 2009, 21(1):145-150.

    [31] Nakashima K, Tran LSP, Dong VN,etal. Functional analysis of a NAC-type transcription factor OsNAC6 involved in abiotic and biotic stress-responsive gene expression in rice[J]. Plant Journal, 2007, 51(4):617-630.

    [32] Takasaki H, Maruyama K, Kidokoro S,etal. The abiotic stress-responsive NAC-type transcription factor OsNAC5 regulates stress-inducible genes and stress tolerance in rice[J]. Mgg Molecular & General Genetics, 2010, 284(3):173-183.

    (責任編輯:徐玉秀)

    Digital Gene Expression Profiling Analysis ofAvicenniamarinaunder Heavy Metal Stress Treatment

    YUANChang-chun1,MOJian-xin1,YUANGLiu-jiao1,YURu-feng1,LIUQian2,ZHONGJun-di1,LIUKai-dong1

    (1.Life Science and Technology School, Lingnan Normal University, Zhanjiang 524048, Guangdong, China; 2.Guangzhou Gene Denovo Company, Guangzhou 510006, Guangdong, China)

    [Objective]Transcriptome sequencing was perform to study the molecular mechanism in response to heavy metal stress inAvicenniamarina.[Method]Digital Gene Expression Profiling (DGE) technique was used to analyze the differences in gene expressions ofA.marinaunder artificial heavy metal stress.[Result]Compared to the controls, 10459 differentially expressed genes (DEGs) were detected, including 8685 up-regulated genes and 1774 down-regulated genes. MeV cluster analysis indicated that the expression level of most genes was largely induced in artificial heavy metal treatment sample than that in control sample. Gene Ontology analysis showed that most DEGs locate in “organelle envelope”, “plastid envelope” and “chloroplast part”, and were enriched in the biological process of “response to stimulus”, and “l(fā)ocalization”, etc. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed that these genes are distributed in 126 pathways involved in “ribosome”, “photosynthesis”, and various metabolism pathways, such as pyruvate and glyoxylate. Several terpenoid and flavonoid biosynthesis-related enzymes, such as geranylgeranyl pyrophosphate synthase (Unigene0009238) and flavonol synthase (Unigene0001833), were up-regulated by heavy metal stress, which might be the reason for the enhancement for the active ingredients accumulation inA.marina. Furthermore, heavy metal stress significantly induced the cytokinin dehydrogenase 7(CKX7) (Unigene0057124) and serine/threonine-protein kinase (BSK) (Unigene0043741). Besides, the transcription factor analysis demonstrated that bHLH, NAC, MYB-related and WRKY play vital roles in response to heavy metal stress.Furthermore, eight randomly selected genes were used to confirm the accuracy of DGE by qRT-PCR method.[Conclusion]Heavy metal stresses have effects on the expression of related genes involved in the metabolism, photosynthesis, small molecule acid synthesis, hormone synthesis and signal transduction, and secondary metabolites synthesis ofA.marina

    Avicenniamarina; heavy metal stress; digital gene expression profiling (DGE); differentially expressed genes (DEGs); mechanism

    2016-06-12 基金項目: 廣東省林業(yè)科技創(chuàng)新項目(2013KJCX011-03; 2015KJCX025);國家級大學生創(chuàng)新創(chuàng)業(yè)訓練計劃項目(201510579277);湛江市熱帶特色資源植物技術開發(fā)重點實驗室項目(2014A06008);嶺南師范學院自然科學研究項目(LZL1507) 作者簡介: 袁長春,博士,教授. 主要研究方向:植物分子生物學. E-mail:yuanchangchun@163.com * 通訊作者:劉鍇棟,碩士,副研究員. 研究方向:植物分子生物學. E-mail: liukaidong2001@126.com

    10.13275/j.cnki.lykxyj.2017.02.004

    S718.46

    A

    1001-1498(2017)02-0206-08

    猜你喜歡
    白骨紅樹測序
    杰 Sir 帶你認識宏基因二代測序(mNGS)
    新民周刊(2022年27期)2022-08-01 07:04:49
    紅樹林誕生記
    二代測序協(xié)助診斷AIDS合并馬爾尼菲籃狀菌腦膜炎1例
    傳染病信息(2021年6期)2021-02-12 01:52:58
    紅樹植物無瓣海桑中重金屬元素的分布與富集特征
    海洋通報(2020年5期)2021-01-14 09:27:04
    無題
    滇池(2019年1期)2019-02-14 02:36:54
    紅樹植物老鼠簕內生真菌Aspergillus flavipes AIL8化學成分的研究
    媽媽的眼睛
    藝術評鑒(2016年13期)2016-05-30 16:19:05
    執(zhí)子之手
    紅豆(2016年6期)2016-05-14 12:12:30
    漫畫集錦
    日韩强制内射视频| 午夜福利18| 亚洲欧美日韩高清专用| 久久精品影院6| 日韩精品青青久久久久久| 熟女电影av网| 不卡视频在线观看欧美| 欧美一区二区亚洲| 色在线成人网| 久久久久国产精品人妻aⅴ院| 精品99又大又爽又粗少妇毛片 | 精品久久国产蜜桃| 国产高潮美女av| 午夜久久久久精精品| 一级毛片久久久久久久久女| 国产白丝娇喘喷水9色精品| 极品教师在线视频| 日日干狠狠操夜夜爽| 色在线成人网| 欧美日韩瑟瑟在线播放| 日韩欧美在线二视频| 亚洲欧美日韩卡通动漫| www日本黄色视频网| 国产乱人伦免费视频| 观看免费一级毛片| 美女cb高潮喷水在线观看| 国产亚洲91精品色在线| 欧美+亚洲+日韩+国产| 可以在线观看的亚洲视频| 欧美性猛交╳xxx乱大交人| 国产av在哪里看| 中国美白少妇内射xxxbb| 国产亚洲精品综合一区在线观看| 给我免费播放毛片高清在线观看| 日韩欧美三级三区| 亚洲av中文av极速乱 | 欧美高清性xxxxhd video| 99久久精品国产国产毛片| 国产成年人精品一区二区| 精品一区二区三区人妻视频| 国产探花在线观看一区二区| 久久国内精品自在自线图片| 丰满人妻一区二区三区视频av| 欧美三级亚洲精品| 麻豆一二三区av精品| 三级毛片av免费| 精品免费久久久久久久清纯| 婷婷色综合大香蕉| 亚洲人成网站在线播放欧美日韩| 蜜桃久久精品国产亚洲av| 一本久久中文字幕| 欧美一区二区国产精品久久精品| 日韩欧美国产一区二区入口| 亚洲在线自拍视频| 午夜激情欧美在线| 最近最新免费中文字幕在线| avwww免费| 国产精品电影一区二区三区| av在线老鸭窝| 亚洲欧美清纯卡通| 国产精品日韩av在线免费观看| 2021天堂中文幕一二区在线观| 夜夜看夜夜爽夜夜摸| 色综合婷婷激情| 成人鲁丝片一二三区免费| 久久久久九九精品影院| 国产国拍精品亚洲av在线观看| 国产爱豆传媒在线观看| 久久精品国产亚洲av涩爱 | 免费人成视频x8x8入口观看| 老熟妇仑乱视频hdxx| 欧美高清成人免费视频www| 亚洲av不卡在线观看| 性色avwww在线观看| 91久久精品国产一区二区成人| 国产黄a三级三级三级人| 成人av一区二区三区在线看| 欧美+日韩+精品| 特大巨黑吊av在线直播| 精品久久久久久久末码| 九色成人免费人妻av| 久久久久久久久大av| 久久久久久九九精品二区国产| 18+在线观看网站| 内地一区二区视频在线| 美女cb高潮喷水在线观看| 国产精品爽爽va在线观看网站| 99热只有精品国产| 色在线成人网| 欧美三级亚洲精品| 黄色配什么色好看| 日韩一区二区视频免费看| 亚洲一区高清亚洲精品| 国模一区二区三区四区视频| 天堂√8在线中文| 午夜福利欧美成人| 在线天堂最新版资源| 亚洲欧美日韩东京热| 九色国产91popny在线| 亚洲人成伊人成综合网2020| 99九九线精品视频在线观看视频| 久久亚洲精品不卡| 国产精品综合久久久久久久免费| 成人欧美大片| 国产亚洲av嫩草精品影院| 久久精品国产亚洲av香蕉五月| 久久热精品热| 神马国产精品三级电影在线观看| 亚洲av五月六月丁香网| 午夜日韩欧美国产| 欧美+日韩+精品| 欧美日韩瑟瑟在线播放| 国产高清不卡午夜福利| 国产精品亚洲一级av第二区| 亚洲精品粉嫩美女一区| 国产一级毛片七仙女欲春2| 三级国产精品欧美在线观看| 亚洲中文字幕一区二区三区有码在线看| 国产黄片美女视频| 亚洲av二区三区四区| 99视频精品全部免费 在线| 国产精品亚洲一级av第二区| 在线观看免费视频日本深夜| 免费无遮挡裸体视频| 黄色日韩在线| 国产老妇女一区| 精品一区二区三区视频在线观看免费| 少妇熟女aⅴ在线视频| 欧美绝顶高潮抽搐喷水| 色5月婷婷丁香| 99久久精品一区二区三区| 精品一区二区三区av网在线观看| 精品一区二区三区av网在线观看| 午夜久久久久精精品| 人妻丰满熟妇av一区二区三区| 国产精品电影一区二区三区| 中文字幕久久专区| 午夜精品一区二区三区免费看| 日韩欧美国产在线观看| 久久久久久久久久久丰满 | 欧美中文日本在线观看视频| 99久久成人亚洲精品观看| 亚洲av中文字字幕乱码综合| 伊人久久精品亚洲午夜| 国产成人一区二区在线| 亚洲无线观看免费| 尾随美女入室| 久久精品国产亚洲av涩爱 | 日韩欧美三级三区| 国产欧美日韩精品亚洲av| 国产蜜桃级精品一区二区三区| 成年女人看的毛片在线观看| 国产视频一区二区在线看| 国产熟女欧美一区二区| 日本三级黄在线观看| 久久精品国产亚洲av天美| 国产精品嫩草影院av在线观看 | 97超级碰碰碰精品色视频在线观看| x7x7x7水蜜桃| 一级黄片播放器| 午夜福利成人在线免费观看| 日本在线视频免费播放| 欧美高清性xxxxhd video| 天堂影院成人在线观看| 欧美国产日韩亚洲一区| 国产精品电影一区二区三区| 国产精品人妻久久久久久| 日本熟妇午夜| 婷婷亚洲欧美| 变态另类成人亚洲欧美熟女| 午夜福利高清视频| 欧美又色又爽又黄视频| 亚洲成人中文字幕在线播放| 色综合亚洲欧美另类图片| 精品久久国产蜜桃| 最近中文字幕高清免费大全6 | 桃红色精品国产亚洲av| 国产成人影院久久av| 麻豆成人午夜福利视频| 免费在线观看成人毛片| 内地一区二区视频在线| 欧美一级a爱片免费观看看| 日日夜夜操网爽| 成人美女网站在线观看视频| 亚洲午夜理论影院| 亚洲国产精品合色在线| 悠悠久久av| 久久精品久久久久久噜噜老黄 | 亚洲成人中文字幕在线播放| av国产免费在线观看| 欧美日韩精品成人综合77777| 18禁黄网站禁片免费观看直播| 搡老岳熟女国产| 在线观看一区二区三区| 国产精品一及| 99在线人妻在线中文字幕| 我要看日韩黄色一级片| 91久久精品国产一区二区成人| a在线观看视频网站| 午夜爱爱视频在线播放| 人妻少妇偷人精品九色| 女的被弄到高潮叫床怎么办 | 亚洲精品日韩av片在线观看| 国产精品免费一区二区三区在线| 国产精品久久久久久精品电影| 他把我摸到了高潮在线观看| 男女啪啪激烈高潮av片| 国产免费男女视频| 又爽又黄a免费视频| 日本 av在线| av.在线天堂| 十八禁国产超污无遮挡网站| 一本一本综合久久| 色哟哟·www| 国产一区二区三区av在线 | 午夜免费男女啪啪视频观看 | 亚洲精品粉嫩美女一区| 国产男靠女视频免费网站| 国产白丝娇喘喷水9色精品| 又黄又爽又刺激的免费视频.| 欧美黑人巨大hd| 长腿黑丝高跟| 亚洲aⅴ乱码一区二区在线播放| 亚洲va在线va天堂va国产| 亚洲成人久久爱视频| 亚洲男人的天堂狠狠| 联通29元200g的流量卡| 色综合色国产| 亚洲成人精品中文字幕电影| 国产蜜桃级精品一区二区三区| АⅤ资源中文在线天堂| 舔av片在线| 欧美xxxx黑人xx丫x性爽| 国产69精品久久久久777片| 成人综合一区亚洲| 搡老熟女国产l中国老女人| 国产精品久久久久久久电影| 婷婷精品国产亚洲av在线| 亚洲精品日韩av片在线观看| 美女黄网站色视频| 俄罗斯特黄特色一大片| 春色校园在线视频观看| 亚洲成a人片在线一区二区| 久久国产精品人妻蜜桃| 老熟妇乱子伦视频在线观看| 国产精品一区二区性色av| av国产免费在线观看| 国产亚洲91精品色在线| 日本爱情动作片www.在线观看 | av在线天堂中文字幕| 久久久久国内视频| 很黄的视频免费| 别揉我奶头 嗯啊视频| 成人精品一区二区免费| 免费无遮挡裸体视频| 国内精品美女久久久久久| 成人三级黄色视频| 亚洲av免费在线观看| 亚洲内射少妇av| 国产黄色小视频在线观看| 亚洲av免费高清在线观看| 51国产日韩欧美| 国产精品永久免费网站| 18+在线观看网站| 夜夜看夜夜爽夜夜摸| 无遮挡黄片免费观看| 国产91精品成人一区二区三区| 国国产精品蜜臀av免费| 亚洲七黄色美女视频| 人妻少妇偷人精品九色| 尤物成人国产欧美一区二区三区| 国产白丝娇喘喷水9色精品| 成人三级黄色视频| 久久亚洲真实| 久久久久久久久久久丰满 | 91午夜精品亚洲一区二区三区 | 精品人妻一区二区三区麻豆 | 久久久久性生活片| 国产欧美日韩精品一区二区| 999久久久精品免费观看国产| 亚洲精品色激情综合| 99久久无色码亚洲精品果冻| 99热6这里只有精品| 无人区码免费观看不卡| 99热这里只有是精品在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产色爽女视频免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 我的老师免费观看完整版| 午夜福利欧美成人| 国产成人福利小说| 一个人观看的视频www高清免费观看| 国产精品嫩草影院av在线观看 | 淫妇啪啪啪对白视频| 精品久久久久久久人妻蜜臀av| 亚洲av五月六月丁香网| 日韩欧美国产一区二区入口| 99热网站在线观看| 国产精品女同一区二区软件 | 国产亚洲精品av在线| 九九在线视频观看精品| 人人妻人人澡欧美一区二区| 久久精品国产清高在天天线| av在线天堂中文字幕| 日本一二三区视频观看| 久久精品国产99精品国产亚洲性色| 精品人妻一区二区三区麻豆 | 狠狠狠狠99中文字幕| 丝袜美腿在线中文| 三级毛片av免费| 尤物成人国产欧美一区二区三区| 伦精品一区二区三区| 亚洲一区二区三区色噜噜| 国产精品av视频在线免费观看| 九色成人免费人妻av| 色精品久久人妻99蜜桃| 日本撒尿小便嘘嘘汇集6| 国产精品一区二区免费欧美| 久久精品国产亚洲av涩爱 | 露出奶头的视频| 97热精品久久久久久| 成人鲁丝片一二三区免费| 成人欧美大片| 亚洲国产高清在线一区二区三| 国产真实伦视频高清在线观看 | 国产精品久久电影中文字幕| 亚洲成av人片在线播放无| 成年女人看的毛片在线观看| 午夜福利欧美成人| 露出奶头的视频| 亚洲四区av| 亚洲男人的天堂狠狠| 日韩欧美在线乱码| 国产老妇女一区| 最新中文字幕久久久久| 偷拍熟女少妇极品色| 国产成人影院久久av| 精品日产1卡2卡| 久久久久免费精品人妻一区二区| 麻豆久久精品国产亚洲av| 久久久久久国产a免费观看| 真人一进一出gif抽搐免费| 内地一区二区视频在线| 搡老熟女国产l中国老女人| 日韩欧美 国产精品| 久久中文看片网| 午夜老司机福利剧场| 蜜桃久久精品国产亚洲av| 亚洲国产精品久久男人天堂| 亚洲av不卡在线观看| 免费黄网站久久成人精品| 一卡2卡三卡四卡精品乱码亚洲| 人妻久久中文字幕网| 亚洲av不卡在线观看| 一卡2卡三卡四卡精品乱码亚洲| 成人亚洲精品av一区二区| 精品99又大又爽又粗少妇毛片 | 国产在线男女| netflix在线观看网站| 欧美人与善性xxx| 一个人看视频在线观看www免费| 特级一级黄色大片| 看免费成人av毛片| 国产高清不卡午夜福利| 亚洲国产高清在线一区二区三| 国产爱豆传媒在线观看| 国产成人福利小说| 午夜精品久久久久久毛片777| 亚洲av成人精品一区久久| 亚洲人成网站在线播放欧美日韩| 国产蜜桃级精品一区二区三区| 一个人观看的视频www高清免费观看| 成年女人毛片免费观看观看9| 色5月婷婷丁香| 午夜福利18| 国产伦人伦偷精品视频| 女人被狂操c到高潮| 久久久久免费精品人妻一区二区| 欧美又色又爽又黄视频| 国产一区二区三区在线臀色熟女| 国产中年淑女户外野战色| 国产伦精品一区二区三区四那| 成人av一区二区三区在线看| 久久精品影院6| 可以在线观看的亚洲视频| 精品久久久久久久人妻蜜臀av| 中文字幕av成人在线电影| 麻豆国产97在线/欧美| 国产麻豆成人av免费视频| 午夜免费男女啪啪视频观看 | 国产极品精品免费视频能看的| 88av欧美| 精品一区二区三区视频在线| 99精品在免费线老司机午夜| ponron亚洲| 九九热线精品视视频播放| aaaaa片日本免费| 搞女人的毛片| 日韩精品有码人妻一区| 色精品久久人妻99蜜桃| 日韩欧美国产在线观看| 国产精品一区二区性色av| 色哟哟·www| 免费在线观看成人毛片| 国产精品99久久久久久久久| 少妇被粗大猛烈的视频| 亚洲五月天丁香| 真实男女啪啪啪动态图| 永久网站在线| 久久精品影院6| 少妇人妻精品综合一区二区 | 成人国产一区最新在线观看| 久久久精品欧美日韩精品| 在线免费观看不下载黄p国产 | 国产成人福利小说| xxxwww97欧美| 白带黄色成豆腐渣| 久久这里只有精品中国| 国产精华一区二区三区| 老女人水多毛片| 欧美人与善性xxx| 搡老岳熟女国产| 99久久无色码亚洲精品果冻| 亚洲天堂国产精品一区在线| 老熟妇仑乱视频hdxx| 亚洲精品国产成人久久av| 少妇人妻精品综合一区二区 | 国产精品一区二区三区四区久久| 一级av片app| 22中文网久久字幕| 亚洲欧美清纯卡通| 狂野欧美白嫩少妇大欣赏| 亚洲,欧美,日韩| 最新中文字幕久久久久| 美女cb高潮喷水在线观看| 可以在线观看毛片的网站| 日韩精品青青久久久久久| 我要搜黄色片| 日本三级黄在线观看| 免费在线观看影片大全网站| 久久久久久大精品| 亚洲精品456在线播放app | 嫩草影院精品99| 九色国产91popny在线| 亚洲欧美日韩无卡精品| 国产人妻一区二区三区在| 国产在线精品亚洲第一网站| 成人国产一区最新在线观看| 日本爱情动作片www.在线观看 | 日本黄色片子视频| 免费av不卡在线播放| 99在线人妻在线中文字幕| 国产亚洲欧美98| 亚洲av二区三区四区| 亚洲国产欧洲综合997久久,| 亚洲一区高清亚洲精品| 欧美最新免费一区二区三区| 国产高清视频在线观看网站| 无人区码免费观看不卡| 亚洲国产精品合色在线| 天堂动漫精品| 一级黄片播放器| 日韩精品中文字幕看吧| 亚洲精品色激情综合| 成人鲁丝片一二三区免费| av在线亚洲专区| eeuss影院久久| 性插视频无遮挡在线免费观看| 亚洲欧美日韩高清在线视频| 日本与韩国留学比较| 一区二区三区免费毛片| 日日摸夜夜添夜夜添小说| 日韩欧美精品v在线| 男人的好看免费观看在线视频| 看黄色毛片网站| 国内揄拍国产精品人妻在线| 在线免费十八禁| 国产精品野战在线观看| 色噜噜av男人的天堂激情| 麻豆成人午夜福利视频| 国产av麻豆久久久久久久| 少妇被粗大猛烈的视频| 别揉我奶头~嗯~啊~动态视频| 欧美日韩瑟瑟在线播放| 亚洲av免费在线观看| 老司机午夜福利在线观看视频| 永久网站在线| 狠狠狠狠99中文字幕| 黄色欧美视频在线观看| 美女免费视频网站| 久久国内精品自在自线图片| 欧美绝顶高潮抽搐喷水| 精品午夜福利在线看| 成熟少妇高潮喷水视频| 一区二区三区高清视频在线| 嫩草影院新地址| 精品久久久噜噜| 在线观看一区二区三区| 少妇人妻精品综合一区二区 | 村上凉子中文字幕在线| 国模一区二区三区四区视频| 亚洲精华国产精华液的使用体验 | 中出人妻视频一区二区| 国内精品美女久久久久久| 岛国在线免费视频观看| 国产激情偷乱视频一区二区| 女同久久另类99精品国产91| 我要搜黄色片| 真实男女啪啪啪动态图| 成人三级黄色视频| 久久久国产成人免费| 又黄又爽又免费观看的视频| 久久久久久大精品| 精品人妻一区二区三区麻豆 | 国产伦精品一区二区三区四那| 久久这里只有精品中国| 麻豆av噜噜一区二区三区| 中文亚洲av片在线观看爽| 91在线精品国自产拍蜜月| 日韩一区二区视频免费看| www.色视频.com| 免费看日本二区| 久久精品人妻少妇| 俺也久久电影网| 看免费成人av毛片| 嫩草影院精品99| 亚洲欧美日韩高清在线视频| avwww免费| 色视频www国产| 人妻夜夜爽99麻豆av| 最近最新免费中文字幕在线| 九九热线精品视视频播放| 亚洲欧美日韩高清在线视频| 黄色女人牲交| 如何舔出高潮| 自拍偷自拍亚洲精品老妇| 99精品久久久久人妻精品| 中文亚洲av片在线观看爽| 亚洲黑人精品在线| 欧美绝顶高潮抽搐喷水| АⅤ资源中文在线天堂| 高清在线国产一区| 日韩中字成人| 日本黄色片子视频| 22中文网久久字幕| 真人一进一出gif抽搐免费| 又紧又爽又黄一区二区| 狠狠狠狠99中文字幕| 韩国av一区二区三区四区| 变态另类丝袜制服| 亚洲第一电影网av| av国产免费在线观看| 桃色一区二区三区在线观看| 国产乱人视频| 午夜福利视频1000在线观看| 婷婷亚洲欧美| 深爱激情五月婷婷| 少妇的逼水好多| 亚洲内射少妇av| 成人美女网站在线观看视频| 国产精品,欧美在线| 亚洲欧美精品综合久久99| 亚洲美女黄片视频| 一级黄片播放器| 91午夜精品亚洲一区二区三区 | 国产精品国产高清国产av| 俄罗斯特黄特色一大片| 亚洲av中文字字幕乱码综合| 天堂动漫精品| 亚洲人成网站高清观看| 免费av毛片视频| 亚洲久久久久久中文字幕| 国产三级在线视频| 麻豆精品久久久久久蜜桃| 国产高清视频在线播放一区| 给我免费播放毛片高清在线观看| 一级黄色大片毛片| 少妇猛男粗大的猛烈进出视频 | 国产乱人伦免费视频| 久久久久久久精品吃奶| 久久精品国产亚洲av涩爱 | 日本在线视频免费播放| 91久久精品电影网| 一个人免费在线观看电影| 黄色丝袜av网址大全| 国产爱豆传媒在线观看| 亚洲国产精品成人综合色| 日韩大尺度精品在线看网址| 欧美日本亚洲视频在线播放| 欧美xxxx黑人xx丫x性爽| 国产伦精品一区二区三区四那| 美女免费视频网站| 99久久精品国产国产毛片| 国产亚洲精品综合一区在线观看| 最新中文字幕久久久久| av女优亚洲男人天堂| 天天一区二区日本电影三级| 亚洲熟妇中文字幕五十中出| 国产v大片淫在线免费观看| 国产午夜福利久久久久久| 制服丝袜大香蕉在线| 国产精品一区二区免费欧美| 精品久久久久久,| 国产美女午夜福利| 欧美日韩黄片免| 亚洲性夜色夜夜综合| 精品一区二区三区av网在线观看| 日本色播在线视频| 国产精品美女特级片免费视频播放器| 亚洲精品一卡2卡三卡4卡5卡| 在线观看美女被高潮喷水网站| 直男gayav资源| 此物有八面人人有两片| 美女cb高潮喷水在线观看| 人妻制服诱惑在线中文字幕|