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

植物響應(yīng)溫度脅迫的代謝組學(xué)研究進(jìn)展

2017-04-04 09:21:32周連玉李園媛王文妮鐘松
山西農(nóng)業(yè)科學(xué) 2017年2期
關(guān)鍵詞:植物研究

周連玉,李園媛,王文妮,鐘松

(青藏高原資源與環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,青海師范大學(xué)生命與地理科學(xué)學(xué)院,青海西寧810008)

植物響應(yīng)溫度脅迫的代謝組學(xué)研究進(jìn)展

周連玉,李園媛,王文妮,鐘松

(青藏高原資源與環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,青海師范大學(xué)生命與地理科學(xué)學(xué)院,青海西寧810008)

在極端溫度脅迫下,植物體內(nèi)會發(fā)生一系列復(fù)雜的生理、生化變化來響應(yīng)外界環(huán)境。代謝組學(xué)已成為系統(tǒng)生物學(xué)的重要組成部分之一,是對某一生物、組織或細(xì)胞中小分子量的代謝產(chǎn)物進(jìn)行定性和定量分析,以揭示生命現(xiàn)象與過程。綜述了植物響應(yīng)溫度脅迫的代謝譜分析,并展望了代謝組學(xué)急需解決的問題。這將為今后引種、馴化、改良和培育抗寒或耐熱植物提供充分的理論依據(jù)。

植物;代謝組學(xué);溫度脅迫

極端溫度、干旱、鹽堿、病原菌等非生物和生物因子嚴(yán)重影響植物的生長和生產(chǎn)力。溫度脅迫是重要的非生物脅迫因子之一,影響植物的生長發(fā)育、生產(chǎn)量和地理分布范圍,誘導(dǎo)植物在分子、細(xì)胞、生理和生化水平上發(fā)生顯著改變[1-2],從而使其克服溫度脅迫環(huán)境。

代謝組學(xué)、基因組學(xué)、轉(zhuǎn)錄組學(xué)和蛋白質(zhì)組學(xué)是系統(tǒng)生物學(xué)研究的重要組成部分。代謝組學(xué)(metabolomics/metabonomics)是繼基因組學(xué)、轉(zhuǎn)錄組學(xué)和蛋白質(zhì)組學(xué)后又一個新發(fā)展起來的組學(xué)。“metabolomics”是由德國科學(xué)家FIEHN在21世紀(jì)初提出的,即“定性、定量分析某一生物或細(xì)胞中所有分子量在1 000以內(nèi)的代謝物”,主要利用色譜-質(zhì)譜聯(lián)用技術(shù)對植物和微生物進(jìn)行相關(guān)研究。并且依據(jù)研究的對象和目的不同,F(xiàn)IEHN[3]將生物代謝分析(metabolite analysis)歸納為4類:目標(biāo)分析(target analysis)、代謝譜分析(metabolic profiling)、代謝物組分析(metabolomics analysis)和代謝指紋分析(metabolic fingerprinting)。“metabonomics”是由1999年英國NICHOLSON定義為“定量測定生命體系中受生理、病理或基因改造而引起的多種代謝物變化”,主要側(cè)重于運(yùn)用核磁共振法研究人或動物在疾病作用下總代謝物的變化情況。代謝組學(xué)的研究也已成為系統(tǒng)生物學(xué)的重要組成部分之一,其研究結(jié)果可揭示出基因組、轉(zhuǎn)錄組和蛋白組受遺傳或環(huán)境因素影響后的變化規(guī)律,更接近細(xì)胞或生物的表型,由此越來越廣泛地被應(yīng)用[3]。代謝組學(xué)以大規(guī)模的實(shí)驗(yàn)方法學(xué)和統(tǒng)計計算分析相結(jié)合為特征,涉及對生物樣品中的代謝物提取和預(yù)處理、代謝組分的定性和定量分析以及對海量數(shù)據(jù)的解析等方面的研究。通過代謝產(chǎn)物變化作為標(biāo)記,進(jìn)行代謝組的綜合研究,有助于找到調(diào)節(jié)代謝的關(guān)鍵因素,增強(qiáng)人們對植物如何響應(yīng)各種逆境脅迫的理解。筆者綜述了近年來代謝譜分析在植物對溫度逆境響應(yīng)中的研究進(jìn)展,并討論了代謝組研究的現(xiàn)狀和今后的發(fā)展方向。

1 植物響應(yīng)低溫脅迫的代謝組學(xué)

代謝可反映生物活體和狀態(tài),隨著不斷改進(jìn)代謝產(chǎn)物的檢測和鑒定技術(shù),植物應(yīng)答非生物脅迫的代謝調(diào)節(jié)得到廣泛的關(guān)注[4-5]。目前,植物逆境防御過程中的代謝組學(xué)已成為廣泛研究的熱點(diǎn)之一,主要是涉及一些模式植物和作物品種[6-9]。在低溫(4℃)脅迫下,擬南芥植株體內(nèi)天冬氨酸、鳥氨酸、精氨酸含量的增加表明了尿素合成途徑被上調(diào);α-酮戊二酸、延胡索酸、蘋果酸、檸檬酸等化合物在低溫脅迫下積累揭示了三羧酸循環(huán)途徑被激活[6]。在后續(xù)的研究中,這些代謝組數(shù)據(jù)和轉(zhuǎn)錄組學(xué)數(shù)據(jù)相互結(jié)合揭示擬南芥低溫適應(yīng)的機(jī)理[10]。外源物質(zhì)對植物響應(yīng)低溫脅迫的代謝譜的作用效果也受到一些學(xué)者關(guān)注,泡葉藻(Ascophyllum nodosum)提取物能促進(jìn)擬南芥在低溫脅迫下可溶性糖、糖醇、有機(jī)酸和脂肪酸的積累,這一結(jié)果也被轉(zhuǎn)錄組學(xué)數(shù)據(jù)分析所證實(shí)[11]。另外,SHI等[12]研究外源鈣對狗牙草(Cynodon dactylon)在正常溫度(28℃)和低溫(4℃)條件下產(chǎn)生的特異性差異蛋白和代謝產(chǎn)物的影響,結(jié)果表明,外源鈣可使51種蛋白質(zhì)的豐度發(fā)生改變,基于氣相-質(zhì)譜聯(lián)用技術(shù)(GC-MS)平臺鑒定出43種代謝物,有39種化合物受到外源鈣的顯著影響,其中,33種代謝物在外源鈣作用下增多。

有關(guān)不同物種間代謝譜的比較研究,SON等[13]應(yīng)用核磁共振技術(shù)(NMR)比較6個葡萄(Vitis vinifera)品種響應(yīng)10℃條件的代謝譜,發(fā)現(xiàn)不同品種間代謝譜不同,3個品種的脯氨酸水平高,另外3個則反之。BENINA等[14]基于GC-MS平臺比較3個物種于21,4℃生長下的代謝譜,正常溫度下巴爾干苣苔(Haberlea rhodopensis)中棉子糖、蜜二糖、海藻糖、鼠李糖、肌醇、山梨糖醇、半乳糖、蘇糖酸、α-酮戊二酸、檸檬酸和甘油等含量高于其他2種;腐胺和延胡索酸含量在擬南芥(Arabidopsis thaliana)中最高;鹽芥(Thellungiella halophyla)中具有高濃度的丙氨酸、β-丙氨酸、天冬酰胺、組氨酸、異亮氨酸、苯丙氨酸、纈氨酸、絲氨酸、蘇氨酸;3個物種對低溫反應(yīng)不一致,鹽芥合成較多的麥芽糖,擬南芥棉子糖合成量增多,而蔗糖含量在3個物種中均上升。VACLAVIK等[15]基于超高效液相-質(zhì)譜聯(lián)用技術(shù)(UPLC-MS)平臺分析8個擬南芥品系的低溫代謝譜,發(fā)現(xiàn)冷敏感型品系中3-丁烯基硫苷含量較高,而抗寒性品系中黃酮-3-醇糖苷合成量增加。總體上,不同物種響應(yīng)低溫脅迫的代謝譜表現(xiàn)不同,也暗示出一些代謝產(chǎn)物在低溫脅迫下發(fā)揮著協(xié)同作用,由此可借助代謝組學(xué)方法來篩選抗寒性品種。

植物是一個復(fù)雜而精致的有機(jī)體,體內(nèi)代謝物的合成與分解總是處于微妙的動態(tài)平衡狀態(tài)以感應(yīng)外部環(huán)境變化,維護(hù)機(jī)體的正常代謝。溫度引起許多代謝產(chǎn)物的變化,尤其是低溫影響更大,幼齡葉片中大約50%的代謝產(chǎn)物于低溫脅迫49 d時表現(xiàn)出顯著的變化,而成熟葉片在低溫脅迫28 d或49 d時代謝產(chǎn)物幾乎無變化[16]。對擬南芥響應(yīng)低溫脅迫過程(0,4,12,24,48,96 h)的代謝產(chǎn)物動態(tài)變化分析,229種化合物的含量呈現(xiàn)出持續(xù)性或瞬間性增加或降低,其中,92種化合物為低溫脅迫初期時增加或降低,66種化合物發(fā)生在低溫脅迫下中間時段(12~24 h),71種化合物在低溫處理后期(48~96 h)降低或增加[6]。ROHLOFF等[17]基于GCMS技術(shù)檢測3個野草莓品系(Fragaria vesca)響應(yīng)2℃下0,3,24,72,240 h時的代謝譜,發(fā)現(xiàn)在低溫脅迫下許多代謝產(chǎn)物于脅迫后期(72,240 h)合成量增加;而根系是植物最敏感的器官之一,其多數(shù)代謝物的水平高于植株地上部分,低溫對3個品系的代謝途徑的影響也不一致。由此可見,植物響應(yīng)低溫脅迫其代謝產(chǎn)物在時間和空間都具有高度的動態(tài)性。

2 植物響應(yīng)高溫脅迫代謝組學(xué)

從代謝組學(xué)方面對植物應(yīng)答高溫脅迫進(jìn)行研究的報道并不多見,KAPLAN等[6]發(fā)現(xiàn)高溫脅迫影響擬南芥體內(nèi)氨基酸、碳水化合物等代謝物的變化。許多氨基酸如天冬酰胺、亮氨酸、異亮氨酸、蘇氨酸、丙氨酸和纈氨酸,是由丙酮酸和草酰乙酸合成,這些氨基酸在高溫脅迫下增加;延胡索酸和蘋果酸是草酰乙酸的前體物,也在高溫脅迫下增加;此外,高溫會促進(jìn)少數(shù)含氮化合物如β-丙氨酸、γ-氨基丁酸和腐胺含量上升;碳水化合物如麥芽糖、蔗糖、半乳糖醇、肌醇、棉子糖和單糖在高溫脅迫下積累。RIZHSKY等[18]研究結(jié)果與上述相似,擬南芥體內(nèi)β-丙氨酸、甘油、麥芽糖、蔗糖和海藻糖在高溫6 h時增加。由此可見,植物對高溫脅迫的響應(yīng)方式是復(fù)雜多樣的,涉及多種氨基酸、碳水化合物及一些含氮化合物。植物在高溫脅迫時常常經(jīng)歷著干旱脅迫,RIZHSKY等[18]使用GC-MS技術(shù)研究擬南芥應(yīng)答高溫和干旱復(fù)合脅迫下的代謝物變化,復(fù)合脅迫下植物積累可溶性糖如蔗糖、麥芽糖、海藻糖、果糖和葡萄糖等,脯氨酸含量僅在干旱脅迫下積累,而在高溫和干旱復(fù)合脅迫或高溫脅迫下脯氨酸含量均不增加。

DU等[19]通過GC-MS技術(shù)對狗牙根和草地早熟禾響應(yīng)高溫脅迫的代謝譜進(jìn)行分析,發(fā)現(xiàn)二者的36種代謝物發(fā)生改變,狗牙根尤其在脅迫18 d時,大多數(shù)代謝物的水平比草地早熟禾高,這些代謝物有7種糖類(蔗糖、果糖、半乳糖、紅藻糖、蜜二糖、麥芽糖、木糖)、1種糖醇(肌醇)、6種有機(jī)酸(蘋果酸、檸檬酸、異檸檬酸、蘇糖酸、半乳糖醛酸、甲基丙二酸)和9種氨基酸(天冬酰胺、丙氨酸、纈氨酸、蘇氨酸、γ-氨基丁酸、異亮氨酸、甘氨酸、賴氨酸和蛋氨酸)。該研究表明,對高溫脅迫的響應(yīng),不同物種之間代謝譜不同,不同脅迫階段表現(xiàn)的代謝譜也不一樣,植物體內(nèi)往往積累多種氨基酸和糖類物質(zhì)以增強(qiáng)植株的耐熱性。

3 溫度脅迫代謝譜的比較

植物對高溫和低溫反應(yīng)的區(qū)別在于,低溫脅迫對代謝譜的影響更大于高溫脅迫,在低溫脅迫早期、中期、后期過程中,代謝產(chǎn)物相對均衡地發(fā)生變化;而在高溫脅迫過程時,最初30 min內(nèi)大多數(shù)代謝產(chǎn)物發(fā)生了明顯的變化;其次為高溫脅迫所產(chǎn)生的大部分代謝物也會在低溫脅迫下產(chǎn)生。KAPLAN等[6]通過GC-MS技術(shù)對響應(yīng)40,4℃的擬南芥植株間的代謝產(chǎn)物進(jìn)行比較,發(fā)現(xiàn)響應(yīng)高溫脅迫所產(chǎn)生的大部分代謝物與低溫脅迫產(chǎn)生的代謝物相同,有143種代謝產(chǎn)物在高溫脅迫下發(fā)生顯著變化;311種代謝物在低溫脅迫下發(fā)生顯著改變,而其中大約60%代謝產(chǎn)物不受高溫影響。WIENKOOP等[20]使用GCMS檢測擬南芥植株響應(yīng)熱脅迫和冷脅迫的代謝研究,也發(fā)現(xiàn)響應(yīng)高溫脅迫所產(chǎn)生的大部分代謝物與低溫脅迫產(chǎn)生的代謝物相同,此外,結(jié)果還表明,脯氨酸、單糖(葡萄糖、果糖)、肌醇半乳糖苷和棉子糖等物質(zhì)對溫度脅迫具有重要作用。

總體上,植物體內(nèi)的氨基酸、三羧酸循環(huán)的中間代謝物、許多碳水化合物含量均會受到溫度脅迫的影響[10,16,18,19-22]。

4 展望

植物對溫度脅迫產(chǎn)生耐受性的過程中涉及諸多代謝物的變化。近10 a內(nèi),從代謝水平上闡述生物響應(yīng)不良環(huán)境的機(jī)理研究受到關(guān)注。代謝組學(xué)是針對所有代謝物進(jìn)行分析,但是如今仍然沒有一種能適合所有代謝產(chǎn)物的提取和預(yù)處理方法,實(shí)際研究中只能根據(jù)不同的研究對象、目的和檢測手段而選擇不同的提取方法。代謝產(chǎn)物的檢測、分析和鑒定是代謝組學(xué)的核心技術(shù),但是目前檢測分析技術(shù)和計算技術(shù)均不夠成熟和完善。大多數(shù)檢測技術(shù)最多只能分離出數(shù)百種代謝物,在實(shí)際操作中往往需要采用多種方法進(jìn)行分析,然而不同檢測技術(shù)或不同儀器所獲得的數(shù)據(jù)之間的轉(zhuǎn)換、整合分析及誤差校驗(yàn)等還沒有準(zhǔn)確而通用的方法,因此,今后要不斷完善和提高代謝組學(xué)分析系統(tǒng)的硬件、軟件技術(shù),朝著整合化、自動化和高通量方向發(fā)展。代謝產(chǎn)物分析譜圖包含大量的信息,為了挖掘數(shù)據(jù)中潛在的信息,需要運(yùn)用多種數(shù)學(xué)和統(tǒng)計學(xué)的分析工具,如主成分分析法[23]、聚類分析[24]、偏最小二乘判別分析[25-26]和辨別功能分析[27]等方法,而多種數(shù)據(jù)分析技術(shù)的開發(fā)研究十分重要。

生物適應(yīng)脅迫環(huán)境是一種許多代謝產(chǎn)物和代謝途徑組成的動態(tài)變化的復(fù)雜過程。不同物種積累不同的代謝產(chǎn)物,代謝途徑也常常涉及到多種細(xì)胞器的協(xié)作,代謝組學(xué)技術(shù)應(yīng)用到植物研究中還需要建立植物代謝的專門數(shù)據(jù)庫以及代謝模擬工具。為了更清晰地解釋抗逆性代謝物變化的機(jī)理,今后的研究需要整合代謝組學(xué)、轉(zhuǎn)錄組學(xué)和基因組學(xué)等數(shù)據(jù)[28]。

[1]KOTAK S,LARKINDALE J,LEE U,et al.Complexity of the heat stress response in plants[J].Current Opinion in Plant Biology,2007,10(3):310-316.

[2]THOMASHOWMF.Plant cold acclimation:freezingtolerance genes and regulatory mechanisms[J].Annual Reviewof Plant Physiology,1999,50:571-599.

[3]FIEHN O.Metabolomics-the link between genotypes and phenotypes[J].Plant Mol Biology,2002,48(1/2):155-171.

[4]SHULAE V,CORTES D,MILLER G,et al.Metabolomics for plant stress response[J].Physiologia Plantarum,2008,132(2):199-200.

[5]GUY C,KAPLAN F,KOPKA J,et al.Metabolomics of temperature stress[J].Physiologia Plantarum,2008,132(2):220-235.

[6]KAPLAN F,KOPKA J,HASKELL D W,et al.Exploring the temperature-stress metabolome of Arabidopsis[J].Plant Physiology,2004,136(4):4159-4168.

[7]COOK D,F(xiàn)OWLER S,F(xiàn)IEHN O,et al.A prominent role for the CBF cold response pathway in configuring the low-temperature metabolome of Arabidopsis[J].Proceedings of the National Academy of Sciences of the United States of Ameri ca,2004,101(42):15243-15248.

[8]WITT S,GALICIA L,LISEC J,et al.Metabolic and phenotypic responses of greenhouse-grown maize hybrids to experimentally controlled drought stress[J].Molecular Plant,2012,5(2):401-417.

[9]BOWNE J B,ERWIN T A,JUTTNER J,et al.Drought responses of leaf tissues from wheat cultivars of differing drought tolerance at the metabolite level[J].Molecular Plant,2012,5(2):418-429.

[10]KAPLANF,KOPKAJ,SUNGDY,et al.Transcript and metabolite profilingduringcold acclimation of Arabidopsis reveals an intricate relationship ofcold-regulated gene expression with modifications in metabolite content[J].The Plant Journal,2007,50(6):967-981.

[11]NAIR P,KANDASAMY S,ZHANG J Z,et al.Transcriptional and metabolomic analysis of Ascophyllum nodosum mediated freezing tolerance in Arabidopsis thaliana[J].BMC Genomics,2012,13:643.

[12]SHI H T,YE T T,ZHONG B,et al.Comparative proteomic and metabolomic analyses reveal mechanisms of improved cold stress tolerance in bermudagrass(Cynodon dactylon(L.)Pers.)by exogenous calcium[J].Journal of Integrative Plant Biology,2014,56(11):1064-1079.

[13]SON H S,LIM K S,CHUNG H J,et al.Metabolic phenotyping of berries in different sixgrape(Vitis vinifera)cultivars[J].Journal of the Korean Society for Applied Biological Chemistry,2014,57(4):491-502.

[14]BENINA M,OBATA T,MEHTEROV N,et al.Comparative metabolic profiling of Haberlea rhodopensis,Thellungiella halophyla,and Arabidopsis thaliana exposed to low temperature[J]. Frontiers in Plant Science,2013,4:499.

[15]VACLAVIK L,MISHRA A,MISHRA K B,et al.Mass spectrometry-based metabolomic fingerprinting for screening cold tolerance in Arabidopsis thaliana accessions[J].Analytical and Bioanalytical Chemistry,2013,405(8):2671-2683.

[16]GRAY G R,HEATh D.A global reorganization of the metabolome in Arabidopsis during cold acclimation is revealed by metabolic fingerprinting[J].Physiologia Plantarum,2005,124(2):236-248.

[17]ROHLOFF J,KOPKA J,ERBAN A,et al.Metabolite profiling reveals novel multi-level cold responses in the diploid model Fragaria vesca(woodland strawberry)[J].Phytochemistry,2012,77(1):99-109.

[18]RIZHSKYL,LIANGH,SHUMANJ,et al.When defense pathways collide.The response of Arabidopsis to a combination of drought and heat stress[J].Plant Physiology,2004,134(4):1683-1696.

[19]DU H M,WANG Z L,YU W J,et al.Differential metabolic responses of perennial grass Cynodon transvaalensis×Cynodon dactylon(C4)and Poa Pratensis(C3)to heat stress[J].Physiologia Plantarum,2011,141(3):251-264.

[20]WIENKOOP S,MORGENTHAL K,WOLSCHIN F,et al.Intergration of metabolomic and proteomic phenotypes:analysis of data covariance dissects starch and RFO metabolism from low and high temperature compensation response in Arabidopsis thaliana[J]. Molecular and Cell Proteomics,2008,7(9):1725-l736.

[21]HANNAH MA,WIESE D,F(xiàn)REUNDS,et al.Natural genetic variation of freezing tolerance in Arabidopsis[J].Plant Physiology,2006,142(1):98-112.

[22]ANGELCHEVA L,MISHRA Y,ANTTI H,et al.Metabolomic analysis of extreme freezing tolerance in Siberian spruce(Picea obovata)[J].NewPhytologist,2014,204(3):545-555.

[23]RAFTERY D,GU H W,PAN Z Z,et al.Principal component directed partial least squares analysis for combiningnuclear magnetic resonance and mass spectrometry data in metablomics:application to the detection of breast cancer[J].Analytica Chimica Acta,2011,686(1/2):57-63.

[24]XU G W,LI X,LU X,et al.Application of fuzzy c-means clustering in data analysis of metabolomics[J].Analytical Chemistry,2009,81(11):4468-4475.

[25]KEMSLEY E K.Discriminant analysis of high-dimensional data:a comparison of principal components analysis and partial least squares data reduction methods[J].Chemometrics and Intelligent LaboratorySystems,1996,33(1):47-61.

[26]ALSBERG B K,KELL D B,GOODACRE R.Variable selection in discriminant partial least-squares analysis[J].Analytical Chemistry,1998,70(19):4126-4133.

[27]KAWAHARA F K,SANTNER J F,JULIAN E C.Characterization of heavy residual fuel oils and asphalts by infrared spectrophotometry using statistical discriminant function analysis[J].Analytical Chemistry,1974,46(2):266-273.

[28]OKSMAN-CALDENTEYKM,SAITOK.Integratinggenomics and metabolomics for engineeringplant metabolic pathways[J].Current Opinion in Biotechnology,2005,16(2):174-179.

Research Progress in the Metabolomics for Plants Response to Temperature Stress

ZHOULianyu,LI Yuanyuan,WANGWenni,ZHONGSong
(Educational MinistryLaboratoryofEnvironment and Resource in Qinghai-Tibet Platean,College ofLife Science and Geography,Qinghai Normal University,Xining810008,China)

Under acute temperature stress,a complexprocess involvinga number ofphysiological and biochemical changes occured within plants.One promise of systems biology was metabolomics that could be defined as the complete set of small-molecule metabolites within a biology,organ or cell for revealing the process and phenomenon of the life.Here,the paper highlighted a selection of recent publications in the temperature stress field,and the problems that need to be settled in time in metabolomics were prospected.These provided important information for introduction,acclimation,improvement,and breed ofplant cold hardeningor thermotolence.

plant;metabolomics;temperature stress

Q946

A文獻(xiàn)標(biāo)識碼:1002-2481(2017)02-0317-04

10.3969/j.issn.1002-2481.2017.02.42

2016-10-13

青海師范大學(xué)本科生科技創(chuàng)新項(xiàng)目(qhnuxskj2016019);青海師范大學(xué)博士科研啟動經(jīng)費(fèi)(080021323)

周連玉(1976-),女,湖南安仁人,副教授,博士,主要從事微生物學(xué)教學(xué)及科研工作。

猜你喜歡
植物研究
FMS與YBT相關(guān)性的實(shí)證研究
2020年國內(nèi)翻譯研究述評
遼代千人邑研究述論
視錯覺在平面設(shè)計中的應(yīng)用與研究
科技傳播(2019年22期)2020-01-14 03:06:54
EMA伺服控制系統(tǒng)研究
新版C-NCAP側(cè)面碰撞假人損傷研究
植物的防身術(shù)
把植物做成藥
哦,不怕,不怕
將植物穿身上
主站蜘蛛池模板: 亚洲色图欧美视频| 国产精品久线在线观看| 亚洲无限乱码| 国产成人无码播放| 日韩第九页| 欧美a级完整在线观看| 亚洲妓女综合网995久久| 精品国产成人高清在线| 人妻21p大胆| 亚洲最大福利网站| 91香蕉国产亚洲一二三区| 色综合国产| 免费国产小视频在线观看 | 中文国产成人久久精品小说| 亚洲无码视频一区二区三区 | a毛片免费看| 99久久国产自偷自偷免费一区| 欧美国产在线看| 欧美成人免费午夜全| 99re视频在线| 在线欧美日韩| 无码人妻热线精品视频| 色精品视频| 六月婷婷精品视频在线观看| 久久精品只有这里有| 久久先锋资源| 亚洲欧美成人网| 91口爆吞精国产对白第三集| 国产视频a| 无码高潮喷水在线观看| 午夜毛片免费看| 国产草草影院18成年视频| 亚洲男人的天堂在线观看| 欧美有码在线| 久久成人免费| 久久久久亚洲Av片无码观看| 美女一区二区在线观看| 国产成人精品高清不卡在线| 92午夜福利影院一区二区三区| 国产精品一区二区国产主播| 精品亚洲麻豆1区2区3区| 精品国产网站| 99热这里都是国产精品| 亚洲色婷婷一区二区| 91久久偷偷做嫩草影院电| 日韩精品久久久久久久电影蜜臀| 日韩一二三区视频精品| 久久狠狠色噜噜狠狠狠狠97视色| 欧美成人精品在线| 1769国产精品视频免费观看| 久久精品这里只有精99品| 欧美日韩导航| 91九色最新地址| www.日韩三级| 色婷婷综合在线| 亚洲狼网站狼狼鲁亚洲下载| 免费在线观看av| 国产拍揄自揄精品视频网站| 久久semm亚洲国产| 亚洲国产成人综合精品2020| 国产在线拍偷自揄观看视频网站| 国产毛片高清一级国语| 亚洲人成人伊人成综合网无码| 91福利片| 国产一级视频在线观看网站| 91po国产在线精品免费观看| 99尹人香蕉国产免费天天拍| 久青草网站| 亚洲成aⅴ人在线观看| 久久精品国产精品一区二区| 国产永久在线观看| 国产精品乱偷免费视频| 一边摸一边做爽的视频17国产 | 午夜一级做a爰片久久毛片| 四虎成人精品在永久免费| 99偷拍视频精品一区二区| 国产精品亚洲日韩AⅤ在线观看| 欧美人人干| 亚洲人成网站在线播放2019| 午夜欧美在线| 日韩午夜福利在线观看| 尤物在线观看乱码|