999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

蛋白核小球藻Fesod基因的克隆及表達分析

2015-10-24 09:19:22沈佳江靈芝孫雪
生物技術通報 2015年5期
關鍵詞:植物

沈佳 江靈芝 孫雪

(寧波大學海洋學院 浙江省海洋生物工程重點實驗室,寧波 315211)

蛋白核小球藻Fesod基因的克隆及表達分析

沈佳 江靈芝 孫雪

(寧波大學海洋學院浙江省海洋生物工程重點實驗室,寧波315211)

利用RACE技術克隆了蛋白核小球藻Fesod全基因序列,得到1 215 bp序列,其5'非翻譯區長52 bp,3'非翻譯區長452 bp,共編碼236個氨基酸。生物信息學分析結果表明,該Fesod基因編碼蛋白的分子量為26.42 kD,等電點為6.98;位于線粒體中的概率是73.9%;該蛋白沒有信號肽序列;蛋白質二級結構預測結果表明FeSOD中α-螺旋和無規卷曲分別占53.39%和39.41%。再利用實時熒光定量PCR技術檢測了不同鹽度和水楊酸濃度對蛋白核小球藻Fesod基因表達的影響,結果表明Fesod表達量隨著鹽度升高而增加,45‰鹽度表達量為15‰鹽度的3.61倍;而植物激素水楊酸的添加一定程度抑制了45‰鹽度培養藻Fesod基因的表達,并且隨其濃度增加呈現先升后降的趨勢。表明蛋白核小球藻Fesod受鹽度誘導,而水楊酸對其影響不明顯。

蛋白核小球藻;鐵超氧化物歧化酶;轉錄表達;鹽度;水楊酸

超氧化物歧化酶(SOD)是一種普遍存在于植物界中的含金屬原子的抗氧化酶,常被作為研究植物抗逆性的指標[1]。據SOD所結合金屬原子的不同,植物SOD可分為3種類型:MnSOD、Cu/ZnSOD和FeSOD,在低等植物中以FeSOD和MnSOD為主[2]。SOD作為生物體內第一個參與活性氧清除反應的重要保護酶[3],可將超氧陰離子(O2·-)轉化為H2O2,H2O2再在過氧化物酶(POD)和過氧化氫酶(CAT)等作用進一步還原為H2O和O2,從而清除植物體內過多的自由基,以減少自由基對植物的毒害[4]。因此植物可以通過提高SOD等抗氧化酶活性來增強其抗逆性,使植物體在一定程度上能耐受或抵抗高溫、干旱、高鹽等逆境脅迫帶來的危害[1,5]。但逆境脅迫對SOD的影響并不完全一致,如杜氏藻D. tertiolecta在3 mol/L NaCl脅迫下的SOD和CAT等抗氧化酶活性則沒有明顯變化[6]。

此外,植物激素如水楊酸(SA)和脫落酸(ABA)等也能誘導sod基因的表達[7]。水楊酸是一種植物內源信號物質和新型植物激素,大量試驗證明水楊酸具有提高植物抗病害等生物脅迫和鹽度、高低溫等非生物脅迫的作用[8]。其中減輕逆境脅迫引起的抗氧化損傷是水楊酸提高植物抗逆性的重要方式[9]。SA可以提高鹽脅迫下的SOD、POD和CAT等抗氧化酶的活性,從而逆轉鹽度壓力對抗氧化酶的擾亂[10-12]。

小球藻(Chlorella)是一類廣泛應用于食品、醫療、化工等方面的重要經濟微藻,同時小球藻具有較廣泛的耐受性,因此經常用于重金屬、污染物及其它逆境脅迫的研究[13-15]。目前,在低等藻類中關于FeSOD的研究不多,本研究以蛋白核小球藻(C. pyrenoidosa)為材料,克隆該藻的鐵超氧化物歧化酶全基因序列,研究不同鹽度以及水楊酸濃度對Fesod表達變化的影響,旨在豐富低等藻類植物中Fesod基因序列信息,并為FeSOD酶在轉錄水平上響應逆境脅迫研究提供參考。

1 材料與方法

1.1 材料

本實驗所用藻種為蛋白核小球藻820,來自寧波大學海洋生物工程重點實驗室藻種室。

1.2 方法

1.2.1 材料培養與處理 蛋白核小球藻820培養溫度為25℃,光照強度約為40 μmol/(m2·s),光暗周期為12 L/12 D。培養基為添加f/2的人工海水培養基,通過控制NaCl的量來調節培養基鹽度分別為15‰、30‰和45‰;在鹽度為45‰時,添加水楊酸母液至水楊酸終濃度分別為0、0.1、0.5、1.0和2.0 mmol/L。

1.2.2 Fesod部分序列的克隆 根據GenBank數據庫中的蛋白核小球藻(DQ183067)、菜菌衣藻(Chlamydomonas reinhardtii,GQ413964、XM_001690539)和鹽生杜氏藻(Dunaliella salina,AY847684)的Fesod序列,運用MEGA5.1軟件中的ClustalW程序進行同源性比對,設計簡并引物sodF和sodR(表1)。將離心收集的對數生長中期的藻細胞加液氮研磨后,按照Trizol試劑(Invitrogen)說明進行RNA的提取。RNA經過電泳檢測與紫外微量分光光度計(Nanodrop 1000)測定總RNA濃度和純度,再按照PrimeScript RT reagent kit with gDNA eraser反轉錄試劑盒(TaKaRa,大連)的說明轉錄成cDNA,并以該cDNA作為模板進行PCR擴增。擴增反應條件為:94℃預變性5 min;94℃ 50 s,55℃ 1 min,72℃ 1 min,共35個循環;72℃延伸8 min。將電泳檢測擴增大小正確的產物進行回收,用pMD18-T 載體(TaKaRa,大連)進行連接與轉化,最后將經過挑菌檢測為陽性的重組克隆送去測序。所有引物序列的合成和序列測定均由Invitrogen(上海)公司完成。1.2.3 Fesod全基因序列擴增 測序得到的Fesod基因部分序列,利用Primer Premier 5.0設計RACE擴增引物(表1)。用上述相同方法進行RNA提取,以其反轉錄的cDNA為模板,進行5'-RACE和3'-RACE(RACE試劑盒來自TaKaRa,大連)擴增。將擴增產物進行克隆后送去測序。最后將獲得的Fesod部分序列、5'-RACE和3'-RACE序列一起拼接后獲得Fesod全基因序列,并將序列遞交到NCBI數據庫。

表1 Fesod基因的擴增引物

1.2.4 Fesod基因的系統進化分析 在NCBI網站(http://www.ncbi.nlm.nih.gov)中的protein數據庫中查找幾種綠藻的Fesod基因,將它們與本實驗得到的蛋白核小球藻一起用MEGA5.1軟件中的ClustalW程序進行比對,再用鄰接法進行系統進化樹的構建,Bootstrap值為1 000。

1.2.5 Fesod的生物信息學分析 在NCBI網站,通過 ORF finder(http://www.ncbi.nlm.nih.gov/gorf/gorf. html),查找其可能的開放閱讀框。蛋白質分子量和等電點預測使用http://www.expasy.org,亞細胞定位用http://psort.hgc.jp,信號肽預測用http://www. cbs.dtu.dk/services/SignalP,蛋白質跨膜區預測用http://www.cbs.dtu.dk/services/TMHMM,蛋白質二級結構預測用http://npsa-pbil.ibcp.fr/ cgi-bin/npsa_ automat.pl?page=/NPSA/npsa_seccons.html。

1.2.6 Fesod的表達分析 利用Primer Premier5.0軟件設計Fesod熒光定量PCR引物(表1),用18S rDNA作為內標[16]。按照SYBR Premix Ex Taq 說明書進行熒光定量PCR反應,反應體系中SYBR Premix Ex Taq 10 μL,上下游引物各0.4 μL,cDNA 2 μL,加水補足總體積到20 μL。熒光定量PCR(Rotor-Gene 6000 PCR儀)反應條件為:94℃預變性3 min;94℃ 10 s,55℃ 15 s,72℃ 20 s,共40次循環。熒光定量PCR數據采用2-△△Ct法進行分析[17]。

1.2.7 數據處理 數據均采用Excel軟件作圖,用SPSS19.0軟件進行數據分析,用SPSS的獨立樣本T檢驗進行差異顯著性分析,用P<0.05表示差異顯著,P<0.01表示差異極顯著。

2 結果

2.1 蛋白核小球藻Fesod基因的克隆

以蛋白核小球藻cDNA為模板,經過PCR擴增得到約250 bp大小的cDNA片段,測序后其長度為269 bp(圖1-A)。利用GenBank中的BLASTn工具進行搜索發現,該序列與萊茵衣藻Fesod(GenBank登錄號為U22416)和鹽生杜氏藻(AY847684)的Fesod相似性均為75%,與蛋白核小球藻211-8bFesod(DQ183067)的相似性為73%。該比較結果表明獲得片段是蛋白核小球藻的Fesod基因。

根據獲得的269 bp的Fesod部分序列,設計RACE引物進行上下游未知序列的擴增。5'-RACE擴增后,電泳檢測(圖1-B)顯示大約在600 bp的位置有一條PCR條帶,3'-RACE擴增反應后,電泳檢測(圖1-C)顯示在700 bp附近有一條PCR條帶。克隆測序后得到5'-RACE和3'-RACE的擴增產物分別為604 bp和639 bp。拼接后得到長為1 215 bp的全Fesod基因序列。

圖1 蛋白核小球藻Fesod基因擴增結果

2.2 幾株綠藻Fesod基因編碼蛋白的進化分析

在NCBI中利用該蛋白核小球藻Fesod(GenBank登錄號為KF736945,蛋白序列號為AHD25899)進行BLASTx搜索,發現該蛋白核小球藻Fesod編碼氨基酸序列與小球藻C. variabilis(XP_005851580)相似性最高(為68%),與萊茵衣藻Fesod(XP_0016-90591)相似性為64%,而與蛋白核小球藻Fesod(ABA71697)相似性為65%。

我們又從NCBI數據庫中下載了11條綠藻和1條小麥的FeSOD氨基酸序列,與本研究得到的蛋白核小球藻FeSOD共13條序列一起構建系統發育樹。結果(圖2)表明,12條綠藻依次聚在一起,高等植物小麥(與其它綠藻的遺傳距離在0.529-0.548)在最外群。其中6株綠藻——兩株萊茵衣藻(綠藻門、綠藻綱、衣藻目和衣藻科)(XP_001690591和AAB04944),團藻目的雨生紅球藻(Haematococcus pluvialis)(AFC35159)、團藻Volvox carteri f. nagariensis(XP_002946030)和鹽生杜氏藻(AAX92665),以及衣藻科的Polytomella parva(ABH11031)最先聚類;外面依次是石莼綱的裂片石莼(Ulva fasciata)(ABR23159)、Trebouxiophyceae綱的蛋白核小球藻(ABA71697)、Coccomyxa subellipsoidea(XP_00564-7238)、本研究的蛋白核小球藻(AHD25899)和另外兩條小球藻C. variabilis(XP_005851580和XP_005850533)。可見,根據FeSOD氨基酸序列可以把綠藻門的不同綱分開,但是可能不太適合更低分類水平的區分。如蛋白核小球藻(ABA71697)并沒有與同屬的小球藻聚在一起,但該小球藻與裂片石莼的遺傳距離為0.310,與同屬其他3株小球藻的距離為0.329-0.394,差別較小。

圖2 根據FeSOD序列構建的系統進化樹

2.3 蛋白核小球藻Fesod的生物信息學分析

獲得的Fesod基因全長1 215 bp,其中5'-非翻譯區長52 bp,3'-非翻譯區長452 bp,開放閱讀框長711 bp,共編碼236個氨基酸。生物信息學分析結果表明,該Fesod編碼蛋白的理論分子量為26.42 kD,等電點為6.98。亞細胞定位結果表明,該Fesod基因定位于線粒體中的概率是73.9%,位于細胞核中的概率是13.0%,位于細胞質、胞外和細胞骨架的概率均是4.3%。

信號肽預測結果表明該FeSOD蛋白沒有信號肽。跨膜結構預測結果表明編碼的236個氨基酸全部位于胞外。蛋白質二級結構預測結果表明α-螺旋占53.39%,β-折疊占5.93%,無規卷曲所占比例為39.41%,不明確的氨基酸占1.27%(圖3)。

圖3 蛋白核小球藻FeSOD二級結構預測圖

2.4 蛋白核小球藻Fesod基因的轉錄表達分析

圖4顯示,蛋白核小球藻Fesod表達量隨著鹽度升高而增加。以15‰鹽度培養藻的Fesod表達量為1,30‰和45‰鹽度培養藻的Fesod表達量分別是15‰鹽度培養藻的Fesod表達量的1.89倍(P<0.05)和3.65倍(P<0.01)。結果表明,高鹽脅迫一定程度地促進了蛋白核小球藻Fesod的表達,推測高鹽脅迫時,藻細胞通過增加Fesod基因表達從而增加FeSOD酶的含量來對抗高鹽逆境傷害。

植物激素水楊酸的添加對45‰鹽度培養蛋白核小球藻Fesod表達的影響,見圖5(30‰鹽度結果與45‰類似,文中未顯示)。從0.1-2.0 mmol/L SA濃度變化,Fesod表達呈現先下降后上升再下降的趨勢。以未添加SA的Fesod基因表達量作為1,0.1 mmol/L SA組Fesod表達量為對照組的0.54倍(P<0.05),0.5和2.0 mmol/L SA組Fesod表達量分別為對照組的0.84倍和0.80倍,而1.0 mmol/L SA組Fesod表達量為對照的1.17倍。表明水楊酸誘導Fesod表達的作用不明顯,低濃度(0.1和0.5 mmol/L)和高濃度(2.0 mmol/L)水楊酸處理Fesod基因表達均低于對照組,只有1.0 mmol/L SA添加后Fesod表達水平比對照略高,但差異不顯著(P>0.05)。

圖4 鹽度對蛋白核小球藻Fesod轉錄表達的影響

圖5 水楊酸對蛋白核小球藻Fesod轉錄表達的影響

3 討論

SOD是植物抗氧化系統的重要酶,相對于高等植物中Cu/ZnSOD和MnSOD的研究,低等藻類植物中關于SOD的研究并不多,如大型紅藻壇紫菜(Pyropia haitanensis)和龍須菜(Gracilaria lemaneiformis)中Mnsod基因的克隆與熱脅迫下的表達[18,19];微藻中的鈍頂螺旋藻(Spirulina platensis)、雨生紅球藻和萊茵衣藻的Fesod的克隆等[20-22]。FeSOD一般位于植物葉綠體中[23],也有FeSOD位于細胞質中的報道,如苔蘚植物Pogonatum inflexum和大豆[24,25]。但本研究生物信息學分析結果表明該FeSOD蛋白位于線粒體中(其概率是73.9%)。已進行基因組測序的小球藻C. variabilis中有兩條Fesod基因,分別編碼236個和204個氨基酸。該小球藻C. variabilis中編碼236個氨基酸的Fesod(XP_005851580)與本研究小球藻Fesod編碼的氨基酸個數相同,對其序列進行亞細胞定位分析結果表明其FeSOD蛋白位于線粒體和細胞質中的概率分別是34.8%和30.4%。蛋白質二級結構預測結果表明該蛋白核小球藻中α-螺旋占53.39%,β-折疊占5.93%;鈍頂螺旋藻FeSOD的二級結構中α-螺旋所占比例很大,而β-折疊只有一個[20];而雨生紅球藻FeSOD中則含有3個β-折疊片和11個α-螺旋[21],這些結果都符合植物中FeSOD的α-螺旋所占比例遠遠大于β-折疊的序列結構特征[23]。并且,該蛋白核小球藻FeSOD序列也有高等單子葉植物FeSOD蛋白的WEHAYY保守序列[23]。

鹽脅迫是影響植物生命活動的重要環境因素。鹽脅迫會對植物體內活性氧的代謝系統平衡造成破壞,具體表現在活性氧(ROS)含量增加,SOD活性增加。而作為SOD家族的一員,Fesod表達量也受鹽度誘導。如杜氏鹽藻中3.0 mol/L NaCl培養的Fesod表達量明顯比1.5 mol/L NaCl增加[26];100 mmol/L NaCl 處理后銀杏葉中Fesod基因表達量提高為對照組的4.3倍[27]。這些結果與本研究一致。

此外,據報道其他植物激素對高等植物Fesod基因表達的影響不同,如ABA可以誘導銀杏Fesod基因的表達[27],但IAA則無此作用[27];而ABA和IAA不能誘導煙草Fesod的表達,并且ABA作用下煙草Fesod mRNA反而有所降低[28]。本研究結果顯示水楊酸對蛋白核小球藻Fesod表達作用不明顯。

4 結論

本研究從蛋白核小球藻中克隆到1 215 bp的Fesod序列,共編碼236個氨基酸,該蛋白很可能位于線粒體中,無信號肽序列,α-螺旋含量遠高于β-折疊。該Fesod基因表達受鹽度誘導,而植物激素水楊酸對其表達影響較小。

[1] 覃鵬, 劉飛虎, 梁雪妮. 超氧化物歧化酶與植物抗逆性[J].黑龍江農業科學, 2002, 24(1):31-34.

[2] 趙詠梅. 植物SOD在抵抗干旱脅迫中的作用[J]. 生物學教學,2011, 36(3):4-5.

[3] McCord JM, Fridovich I. Superoxide dismutase:An enzymic function for erythrocuprein(hemocuprein)[J]. J Biol Chem, 1969,224(22):6049-6055.

[4] Matés JM. Effects of antioxidant enzymes in the molecular control of reactive oxygen species toxicology[J]. Toxicology, 2000, 153(1-3):83-104.

[5]焉婷婷, 周亞維, 李朋富, 等. 廣鹽性硅藻披針舟形藻在高鹽和低鹽脅迫下的抗氧化響應[J]. 鹽業與化工, 2010, 39(5):15-19.

[6]Jahnke LS, White AL. Long-term hyposaline and hypersaline stresses produce distinct antioxidant responses in the marine alga Dunaliella tertiolecta[J]. Plant Physiol, 2003, 160(10):1193-1202.

[7] 馬旭俊, 朱大海. 植物超氧化物歧化酶(SOD)的研究進展[J].遺傳, 2003, 25(2):225-231.

[8]Hayat Q, Hayat S, Irfan M, et al. Effect of exogenous salicylic acid under changing environment:A review[J]. Environmental and Experimental Botany, 2010, 68(1):14-25.

[9]Parashar A, Yusuf M, Fariduddin Q, et al. Salicylic acid enhances antioxidant system in Brassica juncea grown under different levels of manganese[J]. Int J Biol Macromol, 2014, 70:551-558.

[10]Li T, Hu Y, Du X, et al. Salicylic acid alleviates the adverse effects of salt stress in Torreya grandis cv. Merrillii seedlings by activating photosynthesis and enhancing antioxidant systems[J]. PLoS One, 2014, 9(10):e109492.

[11]Mutlu S, Atici O. Alleviation of high salt toxicity-induced oxidative damage by salicylic acid pretreatment in two wheat cultivars[J]. Toxicol Ind Health, 2013, 29(1):89-96.

[12]Al-Whaibi MH, Siddiqui MH, Basalah MO. Salicylic acid and calcium-induced protection of wheat against salinity[J]. Protoplasma, 2012, 249(3):769-778.

[13]李建宏, 浩云濤, 翁永萍. Cd2+脅迫條件下橢圓小球藻的生理應答[J]. 水生生物學報, 2004, 28(6):659-663.

[14]王淑智, 李利, 張道勇, 等. NaCl與Cd對小球藻光系統Ⅱ(PSⅡ)活性的影響[J]. 應用與環境生物學報, 2011, 17(6):839-846.

[15]許超, 李梅梅, 曾偉, 等. 高效氯氰菊酯對蛋白核小球藻的毒性影響[J]. 環境科學與技術, 2013, 36(12M):6-10, 41.

[16]鄧意龍, 孫雪, 徐年軍, 等. 蛋白核小球藻(Chlorella pyrenoidosa)Rubisco活化酶基因的克隆與表達分析[J]. 海洋與湖沼, 2012, 43(1):41-46.

[17]Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCTmethod[J]. Method, 2001, 25(4):402-408.

[18]Yang R, Liu W, Zhang XL, et al. Sequences of Mn-sod gene from Pyropia haitanensis (Bangiales, Rhodophyta)and its expression under heat shock[J]. Botanica Marina, 2013, 56(3):249-259.

[19]Lu N, Zang X, Zhang X, et al. Gene cloning, expression and activity analysis of manganese superoxide dismutase from two strains of Gracilaria lemaneiformis(Gracilariaceae, Rhodophyta)under heat stress[J]. Molecules, 2012, 17(4):4522-4532.

[20]郭建軍, 龔興國, 魯宇文, 等. 鈍頂螺旋藻Fe-SOD基因的克隆及序列分析[J]. 浙江大學學報:理學版, 2004, 31(6):674-678.

[21]樊寅卯, 張蕾, 秦松, 等. 雨生紅球藻鐵型超氧化物酶基因的克隆與序列分析[J]. 海洋科學, 2012, 36(4):857-865.

[22] McKim SM, Durnford DG. Translational regulation of light-harvesting complex expression during photoacclimation to high-light in Chlamydomonas reinhardtii[J]. Plant Physiol Biochem, 2006,44(11-12):857-865.

[23]李娟, 楊立新, 鄭文寅, 等. 小麥FeSOD基因的克隆與序列分析[J]. 核農學報, 2013, 27(8):1111-1117.

[24] Kanematsu S, Okayasu M, Ueno S. Atypical cytosol-localized Fesuperoxide dismutase in the moss Pogonatum inflexum[J]. Bull Minamikyushu Univ, 2013, 43 A:23-31.

[25] Salin ML. Toxic oxygen species and protective systems of the chloroplast[J]. Physiologia Plantarum, 1988, 72(3):681-689.

[26] 盧雪景, 劉紅濤, 柴玉榮, 等. 外源一氧化氮緩解鹽脅迫下杜氏鹽藻氧化損傷的初步研究[J]. 河南農業科學, 2009, 37(5):47-52.

[27] 李琳玲, 程華, 許鋒, 等. 銀杏鐵型超氧化物歧化酶基因(GbFeSOD)的克隆與表達[J]. 果樹學報, 2009, 26(6):840-846.

[28] Kurepa J, Herouart D, Van MM. Differential expression of CuZnand Fe-superoxide dismutase genes of tobacco during development,oxidative stress, and hormonal treatments[J]. Plant Cell Physiol,1997, 38(4):463-470.

(責任編輯 馬鑫)

Gene Cloning and Expression Analysis of Fesod Gene from Chlorella pyrenoidosa

Shen Jia Jiang Lingzhi Sun Xue
(Key Laboratory of Marine Biotechnology of Zhejiang Province,School of Marine Sciences,Ningbo University,Ningbo315211)

The Fesod gene of Chlorella pyrenoidosa was cloned by RACE method. Total 1 215 bp Fesod sequence was obtained, containing 52 bp 5' untranslated region, 452 bp 3' untranslated region, and an open reading frame of 711 bp. The bioinformatics analysis showed that the molecular weight of FeSOD was 26.42 kD and isoelectric point was 6.98; the FeSOD protein located in the mitochondria with a probability of 73.9%; the FeSOD protein had no signal peptide; and the secondary structure prediction results showed that a-helix and random coil in FeSOD accounted for 53.39% and 39.41%, respectively. Then, the transcriptional expression of Fesod at different concentrations of salinity and salicylic acid was detected by real-time fluorescent quantitative PCR. The results showed that Fesod expression increased with the salinity varying from 15‰ to 45‰, and the expression quantity was 3.61 times under 45‰ salinity culture than that under 15‰ salinity. Moreover, salicylic acid could inhibit the mRNA accumulation of Fesod to a certain extent under 45‰ salinity, and its accumulation firstly increased and then decreased with the salicylic acid concentration from 0.1 to 2.0 mmol/L. In conclusion, the expression of Fesod was induced by salinity, while salicylic acid showed no significant effect on its expression.

Chlorella pyrenoidosa;Fe-superoxide dismutase;transcriptional expression;salinity;salicylic acid

10.13560/j.cnki.biotech.bull.1985.2015.05.027

2014-09-05

浙江省自然科學基金項目(LY13D060007),浙江省科技創新團隊(2010R50025-25)

沈佳,女,碩士研究生,研究方向:藻類生化與分子生物學;E-mail:sj677501@163.com

孫雪,女,博士,副研究員,研究方向:藻類生化與分子生物學;E-mail:sunxue@nbu.edu.cn

猜你喜歡
植物
誰是最好的植物?
為什么植物也要睡覺
長得最快的植物
各種有趣的植物
植物也會感到痛苦
會喝水的植物
植物的防身術
把植物做成藥
哦,不怕,不怕
將植物穿身上
主站蜘蛛池模板: 欧美国产综合色视频| 少妇人妻无码首页| 国产人碰人摸人爱免费视频| 亚洲成在线观看 | 亚洲一区二区三区在线视频| 无码精油按摩潮喷在线播放| 波多野结衣视频网站| 91亚洲视频下载| 人人艹人人爽| 欧美性爱精品一区二区三区| 在线精品亚洲一区二区古装| 亚洲日韩日本中文在线| 亚洲欧洲日韩久久狠狠爱| 亚洲欧美精品一中文字幕| 久久久久无码国产精品不卡| 久久综合九九亚洲一区| 日本道综合一本久久久88| 91精品国产一区自在线拍| 亚洲无码电影| 欧美午夜小视频| 亚洲国产亚洲综合在线尤物| 国产欧美日韩资源在线观看 | 国产青青操| 精品超清无码视频在线观看| 九九这里只有精品视频| 黄色福利在线| 日韩AV手机在线观看蜜芽| 毛片视频网| 久久久久夜色精品波多野结衣| 日本不卡在线| 不卡视频国产| 一区二区影院| 亚洲美女高潮久久久久久久| 国产乱人乱偷精品视频a人人澡| 亚洲综合狠狠| 成人精品免费视频| 国产精品不卡片视频免费观看| 狂欢视频在线观看不卡| 亚洲首页国产精品丝袜| 国产激情无码一区二区APP| 97在线免费| av手机版在线播放| 欧美成人精品一级在线观看| 91国内在线视频| 中文字幕精品一区二区三区视频 | 欧美亚洲激情| 97成人在线视频| 高清视频一区| 亚洲欧美自拍中文| 最新国产午夜精品视频成人| 久久九九热视频| 亚洲欧美一区二区三区图片| 中文字幕乱码二三区免费| 国产午夜精品一区二区三| 成人毛片免费在线观看| 成人亚洲国产| AV不卡在线永久免费观看| 亚洲美女一级毛片| 免费人成在线观看视频色| 四虎影院国产| 国产一级做美女做受视频| 欧美综合中文字幕久久| 综合久久久久久久综合网| 午夜成人在线视频| 97人妻精品专区久久久久| 欧美精品一区在线看| 欧美日韩第三页| 亚洲网综合| 国产亚洲欧美另类一区二区| 精品人妻一区二区三区蜜桃AⅤ| 欧美曰批视频免费播放免费| 亚洲va精品中文字幕| 国产免费高清无需播放器| 亚洲色图欧美视频| 国产成人精品18| 香蕉综合在线视频91| 最新日本中文字幕| 91福利在线观看视频| 欧美日韩福利| 亚洲丝袜第一页| 一本大道香蕉中文日本不卡高清二区| 成年人国产视频|