• <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
    漫畫集錦
    精品欧美一区二区三区在线| 国产女主播在线喷水免费视频网站| 国产野战对白在线观看| 日韩免费高清中文字幕av| 精品久久蜜臀av无| 母亲3免费完整高清在线观看| 午夜免费观看性视频| 自线自在国产av| 真人做人爱边吃奶动态| 在线av久久热| 51午夜福利影视在线观看| 国产主播在线观看一区二区 | a级片在线免费高清观看视频| 国产精品一二三区在线看| 久久午夜综合久久蜜桃| 男人操女人黄网站| 国产精品av久久久久免费| 在线观看免费视频网站a站| 一级黄片播放器| 欧美大码av| 中文字幕人妻熟女乱码| 中文字幕另类日韩欧美亚洲嫩草| 国产xxxxx性猛交| 亚洲欧美日韩另类电影网站| 高潮久久久久久久久久久不卡| 18禁裸乳无遮挡动漫免费视频| 狠狠婷婷综合久久久久久88av| 80岁老熟妇乱子伦牲交| 国产精品一二三区在线看| 亚洲色图综合在线观看| 久久久久国产一级毛片高清牌| 少妇人妻久久综合中文| 精品免费久久久久久久清纯 | 亚洲精品一区蜜桃| 免费在线观看影片大全网站 | 午夜两性在线视频| 精品一区二区三区av网在线观看 | 午夜老司机福利片| 香蕉丝袜av| 欧美97在线视频| 又紧又爽又黄一区二区| av天堂久久9| 国产成人一区二区三区免费视频网站 | 精品一区二区三区av网在线观看 | 国产精品亚洲av一区麻豆| 亚洲欧美一区二区三区黑人| 赤兔流量卡办理| 亚洲中文字幕日韩| 国语对白做爰xxxⅹ性视频网站| 亚洲国产精品成人久久小说| 97人妻天天添夜夜摸| 咕卡用的链子| 校园人妻丝袜中文字幕| 啦啦啦在线观看免费高清www| 久久久亚洲精品成人影院| 久久人人97超碰香蕉20202| 女人爽到高潮嗷嗷叫在线视频| 我的亚洲天堂| 国产高清videossex| 中文字幕色久视频| 精品国产乱码久久久久久小说| 久久久久久久久久久久大奶| a级片在线免费高清观看视频| 狠狠精品人妻久久久久久综合| 欧美日韩av久久| 黄色毛片三级朝国网站| 一本—道久久a久久精品蜜桃钙片| 国产精品国产三级专区第一集| 黄色视频不卡| 国产成人系列免费观看| 黄网站色视频无遮挡免费观看| 国产免费视频播放在线视频| 最新在线观看一区二区三区 | 欧美大码av| 久久这里只有精品19| 中国美女看黄片| 色94色欧美一区二区| 天堂中文最新版在线下载| 精品一区二区三卡| 欧美激情 高清一区二区三区| 成年美女黄网站色视频大全免费| 国产精品 国内视频| 亚洲欧洲国产日韩| av又黄又爽大尺度在线免费看| 在线观看免费视频网站a站| 久久av网站| 亚洲一码二码三码区别大吗| 国产福利在线免费观看视频| 肉色欧美久久久久久久蜜桃| 成人亚洲欧美一区二区av| 深夜精品福利| 99re6热这里在线精品视频| 你懂的网址亚洲精品在线观看| 成人亚洲欧美一区二区av| 美女中出高潮动态图| 99国产精品一区二区三区| 观看av在线不卡| 18禁观看日本| 国产日韩欧美视频二区| 亚洲精品一区蜜桃| avwww免费| 欧美黄色片欧美黄色片| 欧美黑人精品巨大| 中文字幕人妻丝袜一区二区| 国产在线视频一区二区| 黄色一级大片看看| 成人免费观看视频高清| 国产成人精品在线电影| 欧美激情高清一区二区三区| 亚洲av欧美aⅴ国产| 在线观看免费日韩欧美大片| 又大又黄又爽视频免费| 18禁黄网站禁片午夜丰满| 色精品久久人妻99蜜桃| 久久ye,这里只有精品| 人妻 亚洲 视频| 巨乳人妻的诱惑在线观看| 欧美精品高潮呻吟av久久| 美女午夜性视频免费| av天堂在线播放| 日韩精品免费视频一区二区三区| 人人妻人人爽人人添夜夜欢视频| 捣出白浆h1v1| av有码第一页| 最近中文字幕2019免费版| 肉色欧美久久久久久久蜜桃| 国产成人精品无人区| 伊人亚洲综合成人网| 一级片免费观看大全| 91国产中文字幕| 美女主播在线视频| 最新在线观看一区二区三区 | 欧美日韩成人在线一区二区| 亚洲精品美女久久av网站| 麻豆av在线久日| 欧美日韩福利视频一区二区| 十分钟在线观看高清视频www| 视频在线观看一区二区三区| 久久中文字幕一级| 久久人人97超碰香蕉20202| 欧美精品啪啪一区二区三区 | 最近中文字幕2019免费版| 久久鲁丝午夜福利片| 99国产精品免费福利视频| 日韩精品免费视频一区二区三区| 国产又色又爽无遮挡免| 老汉色av国产亚洲站长工具| 成人亚洲欧美一区二区av| 色精品久久人妻99蜜桃| 欧美日韩视频精品一区| 日韩一卡2卡3卡4卡2021年| av天堂久久9| av天堂久久9| 国产男女内射视频| 在线观看国产h片| 女性被躁到高潮视频| 悠悠久久av| 免费黄频网站在线观看国产| 新久久久久国产一级毛片| 久久国产精品男人的天堂亚洲| 尾随美女入室| 在线精品无人区一区二区三| 亚洲av美国av| 欧美黄色片欧美黄色片| 91麻豆精品激情在线观看国产 | 欧美精品一区二区免费开放| 汤姆久久久久久久影院中文字幕| 久久免费观看电影| 波多野结衣一区麻豆| 亚洲欧美成人综合另类久久久| av又黄又爽大尺度在线免费看| svipshipincom国产片| 欧美中文综合在线视频| 丰满少妇做爰视频| 欧美 亚洲 国产 日韩一| 日韩一区二区三区影片| videos熟女内射| 亚洲精品美女久久久久99蜜臀 | 黑人猛操日本美女一级片| av一本久久久久| 一边亲一边摸免费视频| 亚洲精品久久午夜乱码| 亚洲国产最新在线播放| 欧美精品啪啪一区二区三区 | 成人影院久久| 国产人伦9x9x在线观看| 国产高清videossex| 久久青草综合色| 黄色毛片三级朝国网站| 亚洲国产精品一区三区| 精品卡一卡二卡四卡免费| 可以免费在线观看a视频的电影网站| 最黄视频免费看| 亚洲国产欧美在线一区| 久久人人爽av亚洲精品天堂| 咕卡用的链子| 一级黄片播放器| 午夜视频精品福利| 日韩电影二区| www.999成人在线观看| 香蕉丝袜av| 亚洲熟女精品中文字幕| 亚洲国产av影院在线观看| 国产精品国产三级专区第一集| 免费观看人在逋| 国产成人精品久久久久久| 国产在线免费精品| 宅男免费午夜| 欧美黄色片欧美黄色片| 国产主播在线观看一区二区 | 成在线人永久免费视频| 亚洲九九香蕉| 99国产精品一区二区蜜桃av | 丝袜人妻中文字幕| 中文欧美无线码| 欧美黄色淫秽网站| 久久热在线av| 99国产精品一区二区三区| 自线自在国产av| 女人被躁到高潮嗷嗷叫费观| 亚洲成人免费电影在线观看 | 成人免费观看视频高清| 久久精品国产亚洲av高清一级| 少妇裸体淫交视频免费看高清 | 久久久国产欧美日韩av| 搡老岳熟女国产| 欧美黑人精品巨大| 深夜精品福利| 大陆偷拍与自拍| 伊人亚洲综合成人网| 久久精品亚洲av国产电影网| 国产精品秋霞免费鲁丝片| 免费久久久久久久精品成人欧美视频| 校园人妻丝袜中文字幕| 老司机在亚洲福利影院| 成人国产一区最新在线观看 | 欧美精品一区二区大全| 国产国语露脸激情在线看| 2018国产大陆天天弄谢| 一级毛片女人18水好多 | 捣出白浆h1v1| 精品福利永久在线观看| 精品一区二区三区四区五区乱码 | 日韩av免费高清视频| 久久久久久久精品精品| 一本大道久久a久久精品| 超碰97精品在线观看| 亚洲第一av免费看| 亚洲久久久国产精品| 国产av精品麻豆| 制服诱惑二区| 国产激情久久老熟女| 欧美日韩国产mv在线观看视频| 在线观看www视频免费| 巨乳人妻的诱惑在线观看| 国产亚洲欧美在线一区二区| 亚洲av国产av综合av卡| 777米奇影视久久| 国产成人a∨麻豆精品| 秋霞在线观看毛片| 国产免费又黄又爽又色| 新久久久久国产一级毛片| 一区二区日韩欧美中文字幕| 精品视频人人做人人爽| 久久久精品免费免费高清| 亚洲人成网站在线观看播放| 国产主播在线观看一区二区 | 久久精品亚洲av国产电影网| 超色免费av| 18禁裸乳无遮挡动漫免费视频| av网站在线播放免费| 欧美另类一区| 两个人看的免费小视频| 一边摸一边抽搐一进一出视频| 亚洲av成人不卡在线观看播放网 | 高清欧美精品videossex| 欧美精品啪啪一区二区三区 | www.精华液| 国产av国产精品国产| 国产日韩欧美视频二区| 在线观看免费视频网站a站| 男女床上黄色一级片免费看| 国产精品.久久久| 国产伦人伦偷精品视频| 观看av在线不卡| 啦啦啦在线观看免费高清www| 黄色a级毛片大全视频| 亚洲自偷自拍图片 自拍| 亚洲精品第二区| 日本五十路高清| 黄片小视频在线播放| 香蕉丝袜av| 69精品国产乱码久久久| 亚洲国产av影院在线观看| 久久热在线av| 亚洲成人手机| 纯流量卡能插随身wifi吗| 男人添女人高潮全过程视频| 99国产精品免费福利视频| 大陆偷拍与自拍| 久久精品成人免费网站| 国产97色在线日韩免费| 日日爽夜夜爽网站| 悠悠久久av| 99re6热这里在线精品视频| 国产精品国产三级专区第一集| 国产成人啪精品午夜网站| 免费观看av网站的网址| 欧美精品一区二区大全| 亚洲成国产人片在线观看| 午夜久久久在线观看| 久久精品国产综合久久久| 国产精品一二三区在线看| 男女无遮挡免费网站观看| 国产三级黄色录像| 制服诱惑二区| 久久国产精品人妻蜜桃| 久久精品亚洲熟妇少妇任你| 狠狠精品人妻久久久久久综合| 99国产精品免费福利视频| 欧美亚洲 丝袜 人妻 在线| 国产精品 欧美亚洲| 午夜激情av网站| 国产真人三级小视频在线观看| 中文字幕人妻丝袜制服| 十分钟在线观看高清视频www| 亚洲欧美一区二区三区久久| 国产在线一区二区三区精| 大香蕉久久网| 精品国产乱码久久久久久小说| 在线观看免费高清a一片| 欧美 亚洲 国产 日韩一| 一级黄色大片毛片| 91麻豆精品激情在线观看国产 | 亚洲国产欧美一区二区综合| 黄色片一级片一级黄色片| 老司机午夜十八禁免费视频| 日本wwww免费看| 国产精品欧美亚洲77777| 一级毛片我不卡| 国产野战对白在线观看| 叶爱在线成人免费视频播放| 久久这里只有精品19| 国产成人精品久久二区二区免费| 自线自在国产av| 色94色欧美一区二区| 99热网站在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲av成人不卡在线观看播放网 | 蜜桃在线观看..| 亚洲九九香蕉| 午夜久久久在线观看| 亚洲一区中文字幕在线| 精品亚洲成国产av| 国产一区二区三区综合在线观看| 国产精品 国内视频| 亚洲免费av在线视频| 久久人人爽av亚洲精品天堂| 秋霞在线观看毛片| 美女午夜性视频免费| 叶爱在线成人免费视频播放| 操出白浆在线播放| 精品一区二区三区av网在线观看 | 女人高潮潮喷娇喘18禁视频| 国产三级黄色录像| 在线观看免费午夜福利视频| 人人澡人人妻人| 777久久人妻少妇嫩草av网站| 成年动漫av网址| 天天躁夜夜躁狠狠久久av| 爱豆传媒免费全集在线观看| 高清不卡的av网站| 亚洲精品中文字幕在线视频| 老司机影院毛片| 色婷婷av一区二区三区视频| 夜夜骑夜夜射夜夜干| 18禁国产床啪视频网站| 婷婷色综合大香蕉| 9191精品国产免费久久| av在线老鸭窝| 亚洲欧美一区二区三区黑人| 国产精品欧美亚洲77777| 国产视频一区二区在线看| 欧美人与善性xxx| 亚洲欧洲国产日韩| 我要看黄色一级片免费的| 国产成人系列免费观看| 国产一区二区三区av在线| 欧美日韩亚洲国产一区二区在线观看 | 老熟女久久久| 嫁个100分男人电影在线观看 | 精品少妇黑人巨大在线播放| 国产成人一区二区三区免费视频网站 | 色综合欧美亚洲国产小说| 欧美人与性动交α欧美软件| 午夜影院在线不卡| 欧美精品人与动牲交sv欧美| 国产男人的电影天堂91| 国产片特级美女逼逼视频| 男女无遮挡免费网站观看| 国产精品久久久久久人妻精品电影 | 亚洲精品国产一区二区精华液| 在线 av 中文字幕| av欧美777| 亚洲精品自拍成人| 亚洲精品国产色婷婷电影| 侵犯人妻中文字幕一二三四区| 久久精品aⅴ一区二区三区四区| 天天躁夜夜躁狠狠久久av| 国产视频首页在线观看| 国产精品一二三区在线看| 亚洲国产日韩一区二区| 18禁黄网站禁片午夜丰满| 精品国产国语对白av| 国产日韩欧美亚洲二区| 久久中文字幕一级| 一级片'在线观看视频| 97人妻天天添夜夜摸| 亚洲精品一二三| 亚洲美女黄色视频免费看| 中文字幕av电影在线播放| 一区二区三区乱码不卡18| 午夜免费成人在线视频| 叶爱在线成人免费视频播放| 777米奇影视久久| 一级a爱视频在线免费观看| 免费在线观看完整版高清| 69精品国产乱码久久久| 国产成人精品久久二区二区免费| 一个人免费看片子| 亚洲av国产av综合av卡| 黄色a级毛片大全视频| 18禁观看日本| 欧美激情高清一区二区三区| 2018国产大陆天天弄谢| 国产视频一区二区在线看| 亚洲国产日韩一区二区| 精品久久久久久电影网| 久久精品亚洲熟妇少妇任你| 少妇被粗大的猛进出69影院| 欧美精品一区二区免费开放| 丁香六月欧美| 丝袜喷水一区| 亚洲精品久久久久久婷婷小说| 久久精品久久精品一区二区三区| 亚洲图色成人| 亚洲av欧美aⅴ国产| 我的亚洲天堂| 亚洲av在线观看美女高潮| 一级片免费观看大全| 久久国产亚洲av麻豆专区| www.自偷自拍.com| 青草久久国产| 一本久久精品| 久久精品aⅴ一区二区三区四区| 免费人妻精品一区二区三区视频| 又粗又硬又长又爽又黄的视频| 午夜激情av网站| 亚洲精品美女久久av网站| 观看av在线不卡| 日韩免费高清中文字幕av| 午夜福利一区二区在线看| 亚洲精品自拍成人| 久久 成人 亚洲| 巨乳人妻的诱惑在线观看| 国产精品成人在线| 丝袜美足系列| 日韩欧美一区视频在线观看| 免费观看a级毛片全部| 亚洲国产精品一区三区| 国产在线一区二区三区精| 满18在线观看网站| 午夜免费鲁丝| 国产一区有黄有色的免费视频| 久久免费观看电影| 高潮久久久久久久久久久不卡| 国产日韩欧美亚洲二区| 日韩av在线免费看完整版不卡| 99国产精品99久久久久| 久久精品熟女亚洲av麻豆精品| 精品国产乱码久久久久久小说| www日本在线高清视频| 久久av网站| 在线av久久热| 超碰97精品在线观看| 成在线人永久免费视频| 亚洲精品国产色婷婷电影| 精品人妻1区二区| a级毛片在线看网站| 久久午夜综合久久蜜桃| 精品亚洲成a人片在线观看| 最新的欧美精品一区二区| 亚洲国产最新在线播放| 91麻豆精品激情在线观看国产 | 午夜激情av网站| 亚洲中文av在线| 午夜激情久久久久久久| 天天躁夜夜躁狠狠躁躁| 99国产精品一区二区蜜桃av | 精品一区在线观看国产| 精品亚洲乱码少妇综合久久| 免费高清在线观看视频在线观看| 少妇猛男粗大的猛烈进出视频| 午夜福利一区二区在线看| 亚洲国产av影院在线观看| 久久人人爽人人片av| 99国产综合亚洲精品| 我的亚洲天堂| 一区二区日韩欧美中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 一二三四在线观看免费中文在| 色网站视频免费| 免费在线观看日本一区| 国产淫语在线视频| 18禁黄网站禁片午夜丰满| 又紧又爽又黄一区二区| 国产高清不卡午夜福利| 欧美精品av麻豆av| 日本91视频免费播放| 欧美人与性动交α欧美精品济南到| 水蜜桃什么品种好| 久久中文字幕一级| 国产成人a∨麻豆精品| 中文字幕人妻丝袜一区二区| 一边亲一边摸免费视频| 韩国精品一区二区三区| 免费人妻精品一区二区三区视频| 在现免费观看毛片| 老司机影院毛片| 人人妻人人添人人爽欧美一区卜| 精品国产一区二区三区四区第35| av福利片在线| 性色av一级| 久久免费观看电影| 超碰成人久久| 国产精品一区二区在线观看99| 亚洲精品国产色婷婷电影| 国产精品久久久av美女十八| 日本午夜av视频| 国产精品99久久99久久久不卡| 丝袜美腿诱惑在线| 免费一级毛片在线播放高清视频 | 久久天躁狠狠躁夜夜2o2o | 日韩电影二区| 亚洲av男天堂| 亚洲专区中文字幕在线| 中文欧美无线码| 97精品久久久久久久久久精品| 成年美女黄网站色视频大全免费| 啦啦啦在线免费观看视频4| 国产91精品成人一区二区三区 | 乱人伦中国视频| 麻豆国产av国片精品| 国产精品国产三级国产专区5o| 国产女主播在线喷水免费视频网站| 亚洲欧美一区二区三区久久| 国产欧美亚洲国产| 国产高清videossex| 国产精品香港三级国产av潘金莲 | 另类精品久久| 婷婷色麻豆天堂久久| av不卡在线播放| 国产av精品麻豆| 国产免费视频播放在线视频| 1024香蕉在线观看| 国产精品 欧美亚洲| 日韩视频在线欧美| 久久久久久亚洲精品国产蜜桃av| 好男人电影高清在线观看| 欧美日韩一级在线毛片| 精品国产超薄肉色丝袜足j| 亚洲精品自拍成人| 韩国高清视频一区二区三区| 丝袜人妻中文字幕| 一级黄片播放器| 日韩一本色道免费dvd| 美女国产高潮福利片在线看| 肉色欧美久久久久久久蜜桃| 午夜91福利影院| 色94色欧美一区二区| 国产成人a∨麻豆精品| 国产在线视频一区二区| 久久久久久久久免费视频了| 亚洲成人国产一区在线观看 | 国产野战对白在线观看| 老司机午夜十八禁免费视频| 男人舔女人的私密视频| 欧美日韩福利视频一区二区| 一区二区三区激情视频| 国产成人影院久久av| 97人妻天天添夜夜摸| 亚洲黑人精品在线| 手机成人av网站| 国产福利在线免费观看视频| www日本在线高清视频| 欧美xxⅹ黑人| 中文乱码字字幕精品一区二区三区| 性色av乱码一区二区三区2| 老司机深夜福利视频在线观看 | 成年人免费黄色播放视频| 精品少妇久久久久久888优播| 99re6热这里在线精品视频| 50天的宝宝边吃奶边哭怎么回事| 91精品三级在线观看| 亚洲精品久久午夜乱码| 老汉色av国产亚洲站长工具| 看免费av毛片| 捣出白浆h1v1| 久久国产亚洲av麻豆专区| 成人亚洲精品一区在线观看| 少妇猛男粗大的猛烈进出视频| videosex国产| 热re99久久国产66热| 国产精品av久久久久免费| 一本综合久久免费|