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

    橡膠樹產(chǎn)膠生物學研究進展

    2020-12-09 05:37唐朝榮
    熱帶作物學報 2020年10期
    關鍵詞:橡膠樹

    唐朝榮

    摘? 要:天然橡膠(順式-1,4-聚異戊二烯)是一種重要的工業(yè)原料,主要來自橡膠樹。以天然橡膠的生物合成與產(chǎn)量形成為主要內(nèi)容的產(chǎn)膠生物學研究為橡膠樹高產(chǎn)遺傳改良提供理論指導,近10年取得了重要進展。本文從橡膠樹基因組測序、橡膠轉(zhuǎn)移酶、轉(zhuǎn)基因、以及轉(zhuǎn)錄組和蛋白質(zhì)組等4個方面介紹產(chǎn)膠生物學研究主要進展,并討論了相關領域研究存在的問題,對未來5~10年需重點關注的研究內(nèi)容提出了建議。在介紹橡膠轉(zhuǎn)移酶時,同時概述其他幾種產(chǎn)膠植物的相關研究進展。

    關鍵詞:橡膠樹;產(chǎn)膠生物學;基因組測序;橡膠轉(zhuǎn)移酶

    中圖分類號:S794.1? ? ? 文獻標識碼:A

    Abstract: Natural rubber (cis-1, 4-polyisoprene) is an important industrial raw material, commercially harvested mainly from rubber tree (Hevea brasiliensis). Important progresses have been made in the past decade on Hevea biology of rubber production, focusing on rubber biosynthesis and latex production. The present review summarized the major progress in four research areas with relation to rubber production, i.e. Hevea genome sequencing, rubber transferase, transgenic research, and transcriptomics & proteomics. Moreover, the problems were discussed and a few of research emphasis in the next 5 to 10 years were proposed. When describing the studies of rubber transferase, the recent progress in several other rubber-producing plants were included.

    Keywords: Hevea brasiliensis (para rubber tree); biology of rubber production; genome sequencing; rubber transferase

    DOI: 10.3969/j.issn.1000-2561.2020.10.003

    天然橡膠是一種重要的工業(yè)原料,在交通運輸、醫(yī)療衛(wèi)生、國防軍工等領域應用廣泛,其優(yōu)異的綜合性能迄今仍無法被人工合成橡膠完全替代[1]。我國是世界最大的天然橡膠消費國,年消費量(約550萬t)約占世界天然橡膠總產(chǎn)量(約1400萬t)的40%,但自給率僅約15%[2]。目前,天然橡膠的商業(yè)來源幾乎全部來自巴西橡膠樹這一個熱帶樹種,無論是理論預測[3]或是生產(chǎn)上超高產(chǎn)橡膠樹單株的發(fā)現(xiàn)[4]都顯示橡膠樹的產(chǎn)膠潛力巨大,可在現(xiàn)有最高產(chǎn)橡膠樹品種單產(chǎn)水平(年產(chǎn)約2500 kg/hm2)的基礎上提升3~5倍。產(chǎn)膠生物學基礎研究,將為橡膠樹高產(chǎn)分子育種提供理論依據(jù)與靶標基因。

    近10年,隨著高通量、低成本的二代和三代測序技術以及蛋白質(zhì)組學技術的飛速發(fā)展,橡膠樹產(chǎn)膠生物學研究取得了顯著進展。本文從橡膠樹基因組、橡膠轉(zhuǎn)移酶、轉(zhuǎn)基因以及轉(zhuǎn)錄組和蛋白質(zhì)組學等4個方向介紹產(chǎn)膠生物研究領域的主要進展,并對該研究領域存在的問題和發(fā)展方向提出了思考,以期增進研究者對本領域發(fā)展現(xiàn)狀和未來工作重點的了解。

    1? 橡膠樹基因組測序研究

    2013年,馬來西亞研究者公布了第一個橡膠樹基因組[5],此后包含筆者團隊在內(nèi)的不同研究者又陸續(xù)報道了4個版本[6-9](表1)。從拼接質(zhì)量上看,筆者團隊在2016年發(fā)表的版本[6]僅次于最近以三代測序為主完成的版本[9],但明顯優(yōu)于其他3個同樣以二代測序技術為主的版本(表1)[10]。從對產(chǎn)膠生物學問題的闡述上看,筆者團隊發(fā)現(xiàn)了橡膠樹物種高產(chǎn)橡膠的遺傳線索、提出乙烯刺激產(chǎn)膠新機制等重要產(chǎn)膠生物學問題[6],是迄今最受關注和認可的橡膠樹基因組版本。

    1.1? 橡膠樹基因組的大小

    最早基于孚爾根染色微密度測定法(Feulgen microdensitometry)測定橡膠樹基因組的大小為2.15 Gb[11],而筆者團隊用二代測序k-mer分析法估算橡膠樹‘熱研7-33-97等6個品種的基因組大小在1.41~1.55 Gb之間[6],這與53份橡膠樹栽培種質(zhì)基因組大小的流式細胞分析結果(1.46~1.60 Gb, 平均1.53 Gb)[12]一致。結合幾個較高質(zhì)量橡膠樹基因組的實際拼接大小[6-9](表1),我們認為橡膠樹這個物種基因組的大小應為1.5 Gb左右。

    1.2? 橡膠生物合成單體IPP的來源

    在橡膠樹的產(chǎn)膠細胞乳管中存在2種異戊烯基焦磷酸(IPP)的合成途徑,即細胞質(zhì)胞漿的甲羥戊酸途徑(MVA)和質(zhì)體的2-C-甲基-D-赤蘚糖醇-4-磷酸途徑(MEP)。多數(shù)研究認為MVA途徑可能是橡膠生物合成單體IPP的主要來源[13],但MEP途徑是否參與橡膠生物合成長期存在爭議[14-15]。

    筆者團隊在橡膠樹基因組中鑒定到18個MVA和22個MEP途徑家族基因(圖1A)。組織表達分析發(fā)現(xiàn),MVA途徑的6個催化酶中都至少有1個在膠乳(乳管細胞的細胞質(zhì))中特異高表達的基因,而MEP途徑中僅DXS有2個編碼基因在膠乳中的表達量相對較高,其他6個酶的編碼基因在膠乳中表達量很低或明顯低于其他組織。這些結果表明,MVA途徑是橡膠樹中橡膠生物合成IPP的主要來源,MEP途徑可能對橡膠生物合成的貢獻很小。

    1.3? 橡膠樹高產(chǎn)橡膠的遺傳線索

    在已知的2500余種產(chǎn)膠植物中,為何橡膠樹具有優(yōu)異的產(chǎn)膠能力,筆者團隊從基因組研究上獲得了科學線索。

    通過與其他17種植物基因組比較,我們發(fā)現(xiàn)橡膠樹中橡膠延伸因子/小橡膠粒子蛋白(REF/ SRPP)基因家族發(fā)生了顯著擴增,共有18個基因,數(shù)目遠高于其他植物[6],類似發(fā)現(xiàn)在另外2個關于橡膠樹基因組的研究[7, 9]中也有報道。在橡膠樹基因組上,有12個REF/SRPP基因形成一個跨度約160 kb的基因簇,并且在膠乳中特異高表達的4個基因(REF1、SRPP1、REF3和REF7)都在這個基因簇上;從進化關系上看,膠乳中表達量最高的REF1基因可能是由表達量較低的REF/SRPP基因逐步進化而來(圖1)。從蛋白序

    列比對結果來看,REF1類似于一個C-端被截斷的SRPP1(膠乳表達量僅次于REF1基因)。有趣的是,僅在橡膠樹中存在與REF1編碼蛋白(138 aa)大小相當?shù)腞EF/SRPP基因,而其他植物中只有編碼大小與SRPP1(204 aa)相當或更大的REF/SRPP蛋白。

    按粒徑尺寸將橡膠粒子分為大(>400 nm)、小(<400 nm)兩類,雖然大橡膠粒子的數(shù)目只占膠乳中橡膠粒子總數(shù)的6%,卻貢獻了93%的橡膠產(chǎn)量[16]。SRPP和REF在橡膠粒子上的分布模式存在明顯差異,SRPP主要位于小橡膠粒子膜的表面,而REF主要鑲嵌在大橡膠粒子的膜中[17]。綜合上述結果,我們推測REF/SRPP基因家族的顯著擴增和乳管特異性功能分化,特別是乳管中超高表達特異REF1的出現(xiàn)可能是橡膠樹進化出優(yōu)異高產(chǎn)橡膠性狀的根本原因,這種性狀有利于橡膠樹在熱帶叢林中抵御蛀食性害蟲的侵害[18]。

    1.4? 乙烯刺激橡膠樹產(chǎn)膠的新機制

    在橡膠樹的樹干上涂抹乙烯利(2-氯乙基磷酸,乙烯釋放劑)或直接施用氣體乙烯可刺激橡膠樹排膠、顯著提高膠乳產(chǎn)量,是橡膠樹生產(chǎn)上的一項重要增產(chǎn)措施[18]。在過去40多年中,國內(nèi)外發(fā)表了大量關于乙烯刺激橡膠樹膠乳產(chǎn)生機制的論文,發(fā)現(xiàn)乙烯處理增強了乳管細胞中蔗糖的吸收和降解[19-22]、促進水分吸收[23]和能量代謝[24]、導致胞漿pH堿性化[24]、促進氮同化[25]和引起防御應答[26]等,這些過程與膠乳再生或膠乳排出直接或間接相關,但都未能從根本上回答乙烯處理能顯著刺激橡膠樹乳管產(chǎn)膠的原因。

    利用RNAseq深度測序,筆者團隊共鑒定出500多個受乙烯調(diào)控的膠乳差異表達基因。進一步分析發(fā)現(xiàn),乙烯合成關鍵酶ACO家族基因在膠乳中不表達或表達量極低[6],結合膠乳中氧含量低[19]和ACO催化ACC氧化產(chǎn)生乙烯時需要氧氣的事實,推測橡膠樹乳管細胞中的乙烯合成能力很弱;同時,膠乳中乙烯信號傳遞和應答通路中的多個關鍵基因,包括4個ETR、2個EIN2和1個EIL1基因受乙烯刺激后顯著上調(diào)表達[6],表明乳管細胞中存在活躍的乙烯信號傳遞與應答通路。這一發(fā)現(xiàn),從源頭上回答了在橡膠樹上使用外源乙烯刺激可顯著刺激膠乳增產(chǎn)的根本原因。

    2? 橡膠轉(zhuǎn)移酶研究

    1969年Archer等[27]將催化橡膠生物合成的酶稱為橡膠轉(zhuǎn)移酶(rubber transferase),其催化橡膠生物合成單體IPP逐個摻入到不斷延伸的橡膠烴(順式-1,4-聚異戊二烯)上,分布在膠乳的水相和橡膠粒子上。此后至今的50余年,人們對橡膠轉(zhuǎn)移酶的生化、分子特性和調(diào)控機制的認識不斷深入,特別是在2015—2016年出現(xiàn)了突破性進展。研究者分別以萵苣[28]、短角蒲公英[29]和橡膠樹[30]這3種產(chǎn)膠植物為研究對象,發(fā)現(xiàn)橡膠轉(zhuǎn)移酶不是一種單一蛋白,而是以順式異戊烯基轉(zhuǎn)移酶(cis-prenyl-transferase, CPT)為核心,并包含CPT結合蛋白(一種人NogoB受體類似蛋白,也稱CPT類似蛋白、橡膠轉(zhuǎn)移酶活化劑)、橡膠延伸因子(rubber elongation factor,REF)和小橡膠粒子蛋白(small rubber particle protein,SRPP)等多種蛋白組分,定位在橡膠粒子上催化橡膠烴大分子合成的一種復雜的蛋白復合體[31]。

    2.1? 順式異戊烯基轉(zhuǎn)移酶

    1989年,Light等[32]報道從橡膠樹膠乳中純化了橡膠轉(zhuǎn)移酶,在橡膠粒子存在時可催化橡膠烴的合成,在無橡膠粒子但有二甲基烯丙基焦磷酸(DMAPP)存在時可催化合成牻牛兒基焦磷酸(GPP)和法尼基焦磷酸(FPP)。但Cornish[33]隨后的研究發(fā)現(xiàn)Light等純化的不是橡膠轉(zhuǎn)移酶,而是在膠乳中催化橡膠合成起始分子FPP的可溶性反式異戊烯基轉(zhuǎn)移酶,真正的橡膠轉(zhuǎn)移酶是與橡膠粒子密切結合的順式異戊烯基轉(zhuǎn)移酶(CPT)。Asawatreratanakul等[34]首先克隆了2個在橡膠樹膠乳中特異表達的CPT基因(HRT1和HRT2),其中原核表達的HRT2重組蛋白在膠乳離心洗滌后的底層顆粒存在時可催化合成橡膠大分子。Takahashi等[35]在酵母和擬南芥細胞中表達HRT1和HRT2基因,但發(fā)現(xiàn)HRT重組蛋白并不能催化生成橡膠大分子,推測橡膠轉(zhuǎn)移酶活性需要一些特定的膠乳成分。Post等[36]證明短角蒲公英中乳管特異表達的CPT是橡膠生物合成所必需的,相關基因經(jīng)RNAi沉默后,膠乳中的橡膠含量減低90%以上。利用蛋白質(zhì)組學分析,Dai等[37]從橡膠樹膠乳的橡膠粒子上鑒定到6個CPT蛋白,證實了早期關于CPT與橡膠粒子緊密結合的報道[33]。Tang等[6]發(fā)現(xiàn)橡膠樹基因組上共有11個CPT基因,其中3個在膠乳中高豐度表達,這與Uthup等[38]認為橡膠樹中僅有3個與橡膠生物合成相關的CPT基因(RubCPT1、RubCPT2和RubCPT3)的結論一致。Uthup等[38]同時發(fā)現(xiàn)RubCPT1在不同橡膠品系中的單倍體類型與其基因表達水平和產(chǎn)膠性狀密切相關。最近,Ding等[39]利用轉(zhuǎn)錄組學和基因組學整合分析研究手段發(fā)現(xiàn)1個CPT基因(CPT2)是所鑒定的3個橡膠生物合成樞紐基因之一。

    2.2? 橡膠延伸因子和小橡膠粒子蛋白

    橡膠延伸因子(REF)和小橡膠粒子蛋白(SRPP)是橡膠樹膠乳中豐度最高的蛋白,同屬REF/SRPP蛋白家族。橡膠樹基因組中有18個REF/SRPP基因,是已報道植物基因組中家族基因數(shù)量最多的物種[6],其中一些膠乳特異高表達的家族基因在多個蛋白質(zhì)組研究中被鑒定到[37, 40-42]。關于橡膠樹REF/SRPP蛋白參與橡膠生物合成的研究主要集中在2個膠乳特異高表達家族成員上,即REF1(138 aa, 14.6 kDa)和SRPP1(204 aa, 24 kDa)[6],它們具有較高的氨基酸序列同源性,都主要分布在膠乳中的橡膠粒子上,但與橡膠粒子的結合方式和生化特性存在明顯差別[43-44]:SRPP1松散地結合在小橡膠粒子膜的表面,而REF1同時存在于大、小橡膠粒子上且與膜結合緊密;在溶液中,REF1聚合成富含β-折疊的淀粉樣蛋白,并很快形成微米級別的大聚合體,SRPP1則形成穩(wěn)定的納米級別的近球形多聚體;SRPP1與REF1可發(fā)生相互作用,SRPP1會抑制REF1的聚集。

    Dennis等[13]發(fā)現(xiàn)橡膠粒子上的REF蛋白與橡膠烴分子的比例約1∶1,膠乳中REF蛋白的含量與橡膠含量正相關,同時在體外橡膠生物合成體系中加入REF抗體或去除橡膠粒子上的REF蛋白均會顯著抑制橡膠合成。Priya等[45]發(fā)現(xiàn)REF基因在橡膠樹高產(chǎn)品系膠乳中的表達水平高于低產(chǎn)品系,割膠和乙烯利刺激均能誘導REF基因的表達,表明REF基因表達與膠乳產(chǎn)量正相關。Oh等[46]發(fā)現(xiàn)在體外橡膠生物合成體系中添加SRPP重組蛋白可顯著促進橡膠合成,推測SRPP和REF可能都是橡膠生物合成系統(tǒng)的重要組成部分。

    在2種產(chǎn)膠蒲公英——俄羅斯蒲公英和短角蒲公英中,SRPP/REF被RNA干擾下調(diào)表達后,根中的橡膠含量顯著降低,幅度可高達50%以上[47-49]。不同的是,在俄羅斯蒲公英干擾植株中,橡膠粒子的穩(wěn)定性和橡膠分子量都顯著下降[47];而在短角蒲公英干擾植株中,橡膠粒子的穩(wěn)定性、橡膠分子量和分子分散度則均不受影響[48-49]。Hillebrand等[48]推測,在短角蒲公英中SRPP是維持橡膠粒子穩(wěn)定性的重要蛋白,RNA干擾植株中橡膠含量下降的主要原因是橡膠粒子的穩(wěn)定性受到影響。奇怪的是,在另一種產(chǎn)膠植物萵苣中,RNA干擾膠乳中2個主要SRPP基因并不影響橡膠生物合成,干擾植株膠乳中的橡膠含量、橡膠分子量和分子分散度都不受影響[50]。

    需要指出的是,系統(tǒng)進化分析結果顯示橡膠樹膠乳中特異高表達的REF/SRPP基因單獨聚為一簇,其他產(chǎn)膠植物中一些REF/SRPP基因聚為另一簇,而橡膠樹和其他產(chǎn)膠植物其余的REF/SRPP基因則分散開來與非產(chǎn)膠植物聚為不同的簇[6],這表明其他產(chǎn)膠植物膠乳中關于特異高表達REF/ SRPP基因的研究結果可能難以真實反映橡膠樹中相應基因的功能。最近,Ding等[39]研究鑒定了3個參與橡膠生物合成的樞紐基因,其中包括膠乳中表達豐度最高的2個REF/SRPP基因(REF1和SRPP1)。

    2.3? 橡膠轉(zhuǎn)移酶復合體

    Qu等[28]從萵苣中鑒定到一個在膠乳中高表達的CPT-like 2(CPTL2)蛋白,該蛋白與人NogoB受體和典型的CPT都有一定的同源性,但缺乏CPT的保守motif。萵苣中CPTL2基因經(jīng)RNA干擾后,膠乳中的橡膠含量降至野生植株的5%;CPT3和CPTL2可直接互作,在煙草和酵母細胞中共表達CPT3和CPTL2基因時,原本在細胞溶質(zhì)中表達的CPT3定位到內(nèi)質(zhì)網(wǎng)上;在酵母微粒體中添加CPT3/CPTL2重組蛋白后低聚合度順式-聚異戊二烯的合成能力增強,卻不能合成橡膠大分子。根據(jù)以上結果,作者推測CPTL2可能作為一種腳手架蛋白將CPT3拉到內(nèi)質(zhì)網(wǎng)上,進而形成橡膠生

    物合成的細胞器——橡膠粒子。在短角蒲公英中,Epping等[29]從膠乳蛋白質(zhì)組上鑒定了一個與萵苣CPTL2功能近似的NogoB受體蛋白,該蛋白在橡膠粒子膜上與CPT互作,是橡膠轉(zhuǎn)移酶復合物的必需組分;在該蛋白基因的RNA干擾的植株中,橡膠生物合成能力幾乎完全喪失,膠乳中檢測不到CPT蛋白,但多萜醇(dolichol)的含量和蛋白質(zhì)的糖基化不受影響,因此將該蛋白命名為橡膠轉(zhuǎn)移酶活化劑(rubber transferase activator, RTA)。

    Yamashita等[30]在橡膠轉(zhuǎn)移酶復合體研究中更進了一步,發(fā)現(xiàn)REF1也是橡膠轉(zhuǎn)移酶復合體的一個關鍵組分。利用去污劑(8 mmol/L CHAPS)處理,獲得了去除大部分膜蛋白的橡膠樹橡膠粒子(WRP),結合無細胞蛋白翻譯系統(tǒng)建立了橡膠體外合成反應的研究體系。在WRP懸浮液中同時表達橡膠樹CPT(HRT1)、REF和HRT1-REF橋梁蛋白(HRBP,系萵苣CPTL2和短角蒲公英RTA的同源蛋白)這3種蛋白,發(fā)現(xiàn)橡膠轉(zhuǎn)移酶的活性遠高于單獨表達HRT1或同時表達HRT1和HRBP的活性。橡膠粒子蛋白質(zhì)組學和互作網(wǎng)絡的研究表明,HRT1、REF和HRBP三者可能在橡膠粒子上形成復合體,其中HRBP同時與HRT1和REF互作,在復合體形成中發(fā)揮橋梁作用。基于相關結果,提出了在橡膠樹膠乳中的橡膠粒子上進行橡膠合成以及橡膠粒子發(fā)生機制的模型圖(圖2)[30],指出在橡膠粒子上形成正確的HRT1蛋白復合體是進行橡膠大分子生物合成的關鍵。

    3? 橡膠樹轉(zhuǎn)基因研究

    橡膠樹遺傳轉(zhuǎn)化的效率低、周期長,轉(zhuǎn)基因研究進展緩慢。1994年Arokiaraj等[51]將GUS報告基因?qū)胂鹉z樹的基因組,但其后十幾年僅有幾篇進行橡膠樹轉(zhuǎn)基因體系優(yōu)化研究的報道[52-56]。2012年發(fā)表了第一篇對轉(zhuǎn)基因橡膠樹植株表型和所轉(zhuǎn)基因的功能進行系統(tǒng)研究的文章[57],研究者在橡膠樹中過表達了其自身的細胞質(zhì)胞漿CuZnSOD酶,轉(zhuǎn)基因植株在干旱脅迫處理時氣孔導度下降、脯氨酸含量增加,干旱耐受能力明顯增強,并發(fā)現(xiàn)這種增強與轉(zhuǎn)基因植株中活性氧清除能力的提高直接相關[57]。近幾年,橡膠樹轉(zhuǎn)基因研究取得可喜進展,有望通過轉(zhuǎn)基因手段明顯改進橡膠樹的抗逆和產(chǎn)量性狀[58-60]。

    Lestari等[59]在橡膠樹中過表達1個擬南芥ERF1(ethylene response factor 1)的橡膠樹同源基因HbERF-IXc5,發(fā)現(xiàn)轉(zhuǎn)基因植株的根系發(fā)達,莖圍和株高明顯增加,植株干重顯著增大,這與以往報道的ERF1過表達植株矮化嚴重的現(xiàn)象明顯不同,表明HbERF-IXc5的功能可能有別于經(jīng)典的植物ERF1;HbERF-IXc5轉(zhuǎn)基因植株對干旱、冷、鹽等非生物脅迫的耐受力有所增強,與相關植株生態(tài)生理指標的測定結果一致;轉(zhuǎn)基因植株主葉脈和綠色嫩莖中的初生乳管數(shù)量明顯增加,其中一個轉(zhuǎn)基因株系老化莖中的次生乳管數(shù)量也明顯增多,推測HbERF-IXc5基因可能通過參與乙烯和茉莉酸信號傳導間的“竄擾”(crosstalk),從而控制一些直接參與乳管分化的基因。Jayashree等[60]在橡膠樹中過表達IPP甲羥戊酸合成途徑的關鍵酶HMGR基因(hmgr1),結果令人振奮,所有轉(zhuǎn)基因植株的莖圍和膠乳產(chǎn)量均高于對照植株,膠乳產(chǎn)量最高可達對照的5倍,顯示通過轉(zhuǎn)基因手段培養(yǎng)高產(chǎn)甚至超高產(chǎn)橡膠樹的前景誘人。

    最近,中國熱帶農(nóng)業(yè)科學研究院橡膠研究所的科研人員[61]將體外組合的Cas9/sgRNA核蛋白導入橡膠樹原生質(zhì)體中,實現(xiàn)了對橡膠樹靶標基因FT和TFL1的有效編輯,結合該團隊前期建立的橡膠樹原生質(zhì)體植株再生體系[62],有望將無外源DNA導入的基因編輯技術應用到橡膠樹遺傳改良中。

    4? 橡膠樹轉(zhuǎn)錄組與蛋白質(zhì)組研究

    橡膠樹基因組測序推動了與產(chǎn)膠直接或間接相關的轉(zhuǎn)錄組和蛋白質(zhì)組研究。研究內(nèi)容涉及以下方面:膠乳蛋白質(zhì)組[63]與轉(zhuǎn)錄組的研究方法[10, 64-65];不同橡膠樹組織的轉(zhuǎn)錄組[66-69];膠乳[41, 70]及其不同亞細胞組分(橡膠粒子、c-乳清和黃色體)的蛋白組[37, 42, 63, 71-73];乙烯和茉莉酸刺激[68, 74],以及與死皮相關[75]的膠乳或樹皮轉(zhuǎn)錄組;乙烯刺激的膠乳蛋白質(zhì)組[41, 70];不同產(chǎn)膠水平單株或品系[4, 76]以及排膠相關的膠乳轉(zhuǎn)錄組[77];基于轉(zhuǎn)錄組的產(chǎn)膠相關基因家族發(fā)掘[78-79]、分子標記開發(fā)與遺傳連鎖圖譜構建[61, 80];橡膠樹基因組與轉(zhuǎn)錄組的整合數(shù)據(jù)庫[81]等。這些研究產(chǎn)生了海量的以產(chǎn)膠組織膠乳為主的轉(zhuǎn)錄組和蛋白質(zhì)組數(shù)據(jù),加深了人們對橡膠樹的橡膠生物合成、激素應答和割膠脅迫等生物學問題的認識。但從整體上看,多數(shù)研究只是對相關組學數(shù)據(jù)進行初步分析,未能對所發(fā)現(xiàn)的產(chǎn)膠相關候選基因或蛋白進行深入探究。

    最近幾年,在橡膠樹組學研究上有了一些新變化,取得一些新進展,這里簡要介紹幾項代表性的研究。Wang等[70]鑒定到143個表達豐度顯著受乙烯刺激調(diào)控的膠乳蛋白,進一步利用磷酸化蛋白質(zhì)組學技術鑒定到59個應答乙烯處理的磷酸化蛋白,其中包括一些REF和SRPP的同源異構體(isoform),推測蛋白翻譯后修飾和isoform特異性磷酸化修飾可能在乙烯刺激產(chǎn)膠中發(fā)揮重要作用。Makita等[65]建立了橡膠樹不同組織的全長cDNA文庫,利用Sanger和Illumina兩種測序手段結合更新了5500個基因結構,新注釋了9500個轉(zhuǎn)錄起始位點,結合橡膠樹不同品種和組織的RNAseq數(shù)據(jù),對橡膠生物合成、膠乳產(chǎn)量和抗病有了新認識。Chow等[10]利用PacBio三代測序技術對橡膠樹全長cDNA文庫進行轉(zhuǎn)錄組測序,共得到3.7萬余個平均長度約2 kb的全長轉(zhuǎn)錄本,這些轉(zhuǎn)錄本對應約1.5萬個基因座,與‘熱研7-33-97基因組比對后發(fā)現(xiàn),超過一半的轉(zhuǎn)錄本可能是新的基因isoform。Ding等[82]利用公共數(shù)據(jù)庫中的129個RNAseq數(shù)據(jù)包進行基因共表達分析,共鑒定到25個基因共表達模塊,其中1個模塊基因注釋富集為類異戊二烯代謝,在膠乳中高豐度表達并顯著應答乙烯、茉莉酸處理和死皮發(fā)生,分析認為模塊中的SRPP1、CPT2和REF1是橡膠生物合成的中心(hub)基因;基因組進化分析發(fā)現(xiàn)REF/SRPP基因家族在橡膠樹物種中進化產(chǎn)生2個分別包含SRPP1和REF1基因的特異基因簇,推測與橡膠生物合成密切相關,這些結果與基因組研究[6]發(fā)現(xiàn)“REF1基因可能是橡膠樹物種進化出高產(chǎn)橡膠性狀的關鍵事件”一致。

    5? 問題與展望

    過去10年,橡膠樹產(chǎn)膠生物學研究取得了顯著進展,加深了人們對橡膠樹的產(chǎn)膠機制以及與產(chǎn)膠相關的抗逆與激素應答等重要生物學問題的認識,也為橡膠樹高產(chǎn)分子改良提供了思路和技術儲備。但與水稻等重要糧食作物相比,橡膠樹橡膠產(chǎn)量性狀形成的分子機制研究尚處于比較初級的水平,相關研究成果還無法對橡膠樹高產(chǎn)遺傳改良提供切實支持。根據(jù)橡膠樹高產(chǎn)分子育種研究與實踐發(fā)展的需要,筆者認為在未來5~10年,應對以下幾個方面予以重點關注:

    (1)構建橡膠樹主要栽培品種和核心種質(zhì)的泛基因組及膠乳等主要組織的泛轉(zhuǎn)錄組,揭示橡膠樹基因組物種水平的基因結構、序列變異和可變剪輯模式;

    (2)利用各組學相結合的整合生物學研究手段以及豐富的橡膠樹栽培與種質(zhì)材料,揭示橡膠產(chǎn)量形成的關鍵代謝途徑與節(jié)點基因;

    (3)利用橡膠樹膠乳體外橡膠合成體系,以及酵母或植物懸浮細胞等表達體系,揭示橡膠轉(zhuǎn)移酶復合體的精準結構與作用機制;

    (4)利用橡膠樹懸浮培養(yǎng)細胞以及產(chǎn)膠模式植物——橡膠草的轉(zhuǎn)基因與基因編輯研究,揭示產(chǎn)膠細胞器(橡膠粒子)的發(fā)生與發(fā)育機制;

    (5)探索橡膠樹體胚發(fā)生與植株再生的分子調(diào)控機制,建立橡膠樹高效遺傳轉(zhuǎn)化和基因編輯技術平臺。

    參考文獻

    Van Beilen J B, Poirier Y. Establishment of new crops for the production of natural rubber[J]. Trends in Biotechnology, 2007, 25(11): 522-529.

    莫業(yè)勇. 天然橡膠供需形勢和風險分析[J]. 中國熱帶農(nóng)業(yè), 2019(2): 4-6, 10.

    Paardekooper E. Exploitation of the rubber tree[M]//Webster C, Baulkwill W. Rubber. New York: Longman Scientific and Technical, 1989.

    Tang C, Xiao X, Li H, et al. Comparative analysis of latex transcriptome reveals putative molecular mechanisms underlying super productivity of Hevea brasiliensis[J]. PLoS One, 2013, 8(9): e75307.

    Rahman A Y A, Usharraj A O, Misra B B, et al. Draft genome sequence of the rubber tree Hevea brasiliensis[J]. BMC Genomics, 2013, 14: 75.

    Tang C, Yang M, Fang Y, et al. The rubber tree genome reveals new insights into rubber production and species adaptation[J]. Nature Plants, 2016, 2(6): 16073.

    Lau N S, Makita Y, Kawashima M, et al. The rubber tree genome shows expansion of gene family associated with rubber biosynthesis[J]. Scientific Reports, 2016, 6(1): 28594.

    Pootakham W, Sonthirod C, Naktang C, et al. De novo hybrid assembly of the rubber tree genome reveals evidence of paleotetraploidy in Hevea species[J]. Scientific Reports, 2017, 7: 41457.

    Liu J, Shi C, Shi C C, et al. The Chromosome-based rubber tree genome provides new insights into spurge genome evolution and rubber biosynthesis[J]. Molecular Plant, 2020, 13(2): 336-350.

    Chow K S, Khoo J S, Mohd-Zainuddin Z, et al. Utility of PacBio Iso-Seq for transcript and gene discovery in Hevea latex[J]. Journal of Rubber Research, 2019, 22(4): 169-186.

    Bennett M D, Leitch I J. Nuclear DNA amounts in angiosperms-583 new estimates[J]. Annals of Botany, 1997, 80(2): 169-196.

    柳? 覲, 牛迎鳳, 吳? 裕, 等. 巴西橡膠樹栽培種質(zhì)基因組C值測定和變異分析[J]. 熱帶亞熱帶植物學報, 2018, 26 (5): 523-528.

    Dennis M S, Light D R. Rubber elongation factor from Hevea brasiliensis. Identification, characterization, and role in rubber biosynthesis[J]. The Journal of Biological Chemistry, 1989, 264(31): 18608-18617.

    Sando T, Takeno S, Watanabe N, et al. Cloning and characterization of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway genes of a natural-rubber producing plant, Hevea brasiliensis[J]. Bioscience, Biotechnology and Biochemistry, 2008, 72(11): 2903-2917.

    Chow K S, Matisa M N, Bahari A, et al. Metabolic routes affecting rubber biosynthesis in Hevea brasiliensis latex[J]. Journal of Experimental Botany, 2012, 63(5): 1863-1871.

    Yeang H Y, Yip E, Hamzah S. Characterisation of Zone 1 and Zone 2 rubber particles in Hevea brasiliensis latex[J]. Journal of Natural Rubber Ressearch, 1995, 10: 108-123.

    Berthelot K, Lecomte S, Estevez Y, et al. Rubber particle proteins, HbREF and HbSRPP, show different interactions with model membranes[J]. Biochimica et Biophysica Acta (BBA)-Biomembranes, 2014, 1838(1): 287-299.

    dAuzac J, Jacob J L, Prév?t J C, et al. The regulation of cis-polyisoprene production (natural rubber) from Hevea brasiliensis[M]//Pandalai S G. Recent research developments in plant physiology. Trivandrum: Research Singpost, 1997.

    Tupy J. Sucrose supply and utilization for latex production[M]//DAuzac J, Jacob J-L, Chrestin H. Physiology of rubber tree latex. Boca Raton: CRC Press, 1989.

    Tang C, Huang D, Yang J, et al. The sucrose transporter HbSUT3 plays an active role in sucrose loading to laticifer and rubber productivity in exploited trees of Hevea brasiliensis (para rubber tree)[J]. Plant, Cell and Environment, 2010, 33(10): 1708-1720.

    Dusotoit-Coucaud A, Kongsawadworakul P, Maurousset L, et al. Ethylene stimulation of latex yield depends on the expression of a sucrose transporter (HbSUT1B) in rubber tree (Hevea brasiliensis)[J]. Tree Physiology, 2010, 30(12): 1586-1598.

    Liu S, Lan J, Zhou B, et al. HbNIN2, a cytosolic alkaline/neutral-invertase, is responsible for sucrose catabolism in rubber-producing laticifers of Hevea brasiliensis (para rubber tree)[J]. New Phytologist, 2015, 206(2): 709-725.

    Tungngoen K, Kongsawadworakul P, Viboonjun U, et al. Involvement of HbPIP2; 1 and HbTIP1;1 aquaporins in ethylene stimulation of latex yield through regulation of water exchanges between inner liber and latex cells in Hevea brasiliensis[J]. Plant Physiology, 2009, 151: 843-856.

    Amalou Z, Bangratz J, Chrestin H. Ethrel (ethylene releaser)-induced increases in the adenylate pool and transtonoplast delta pH within Hevea latex cells[J]. Plant Physiology, 1992, 98(4): 1270-1276.

    Pujade-Renaud V, Clement A, Perrotrechenmann C, et al. Ethylene-Induced increase in glutamine synthetase activity and mRNA levels in Hevea brasiliensis latex cells[J]. Plant Physiology, 1994, 105(1): 127-132.

    Putranto R A, Duan C, Kuswanhadi, et al. Ethylene response factors are controlled by multiple harvesting stresses in Hevea brasiliensis[J]. PLoS One, 2015, 10(4): e0123618.

    Archer B L, Cockbain E G. Rubber transferase from Hevea brasiliensis latex[J]. Methods in Enzymology, 1969, 15: 476-480.

    Qu Y, Chakrabarty R, Tran H T, et al. A lettuce (Lactuca sativa) homolog of human Nogo-B receptor interacts with cis-prenyltransferase and is necessary for natural rubber biosynthesis[J]. Journal of Biological Chemistry, 2015, 290(4): 1898-1914.

    Epping J, Van Deenen N, Niephaus E, et al. A rubber transferase activator is necessary for natural rubber biosynthesis in dandelion[J]. Nature Plants, 2015, 1(5): 15048.

    Yamashita S, Yamaguchi H, Waki T, et al. Identification and reconstitution of the rubber biosynthetic machinery on rubber particles from Hevea brasiliensis[J]. eLife, 2016, 5: e19022.

    Cherian S, Ryu S B, Cornish K. Natural rubber biosynthesis in plants, the rubber transferase complex, and metabolic engineering progress and prospects[J]. Plant Biotechnology Journal, 2019, 17(11): 2041-2061.

    Light D R, Dennis M S. Purification of a prenyltransferase that elongates cis-isoprene rubber from latex of Hevea brasiliensis[J]. Journal of Biological Chemistry, 1989, 264(31): 18589-18597.

    Cornish K. The separate roles of plant cis and trans prenyl transferases in cis-1,4-polyisoprene biosynthesis[J]. European Journal of Biochemistry, 1993, 218(1): 267-271.

    Asawatreratanakul K, Zhang Y W, Wititsuwannakul D, et al. Molecular cloning, expression and characterization of cDNA encoding cis-prenyltransferases from Hevea brasiliensis: a key factor participating in natural rubber biosynthesis[J]. European Journal of Biochemistry, 2003, 270(23): 4671- 4680.

    Takahashi S, Lee H J, Yamashita S, et al. Characterization of cis-prenyltransferases from the rubber producing plant Hevea brasiliensis heterologously expressed in yeast and plant cells[J]. Plant Biotechnology, 2012, 29(4): 411-417.

    Post J, van Deenen N, Fricke J, et al. Laticifer-specific cis-prenyltransferase silencing affects the rubber, triterpene, and inulin content of Taraxacum brevicorniculatum[J]. Plant Physiology, 2012, 158(3): 1406-1417.

    Dai L, Kang G, Li Y, et al. In-depth proteome analysis of the rubber particle of Hevea brasiliensis (para rubber tree)[J]. Plant Molecular Biology, 2013, 82(1-2): 155-168.

    Uthup T K, Rajamani A, Ravindran M, et al. Distinguishing CPT gene family members and vetting the sequence structure of a putative rubber synthesizing variant in Hevea brasiliensis[J]. Gene, 2019, 689: 183-193.

    Ding Z, Fu L, Tan D, et al. An integrative transcriptomic and genomic analysis reveals novel insights into the hub genes and regulatory networks associated with rubber synthesis in H. brasiliensis[J]. Industrial Crops and Products, 2020, 153: 112562 .

    Dai L, Nie Z, Kang G, et al. Identification and subcellular localization analysis of two rubber elongation factor isoforms on Hevea brasiliensis rubber particles[J]. Plant Physiology and Biochemistry, 2017, 111: 97-106.

    Tong Z, Wang D, Sun Y, et al. Comparative proteomics of rubber latex revealed multiple protein species of REF/SRPP family respond diversely to ethylene stimulation among different rubber tree clones[J]. International Journal of Molecular Sciences, 2017, 18(5): 958.

    Wang D, Sun Y, Chang L L, et al. Subcellular proteome profiles of different latex fractions revealed washed solutions from rubber particles contain crucial enzymes for natural rubber biosynthesis[J]. Journal of Proteomics, 2018, 182: 53-64.

    Berthelot K, Lecomte S, Estevez Y, et al. Homologous Hevea brasiliensis REF (Hevb1) and SRPP (Hevb3) present different auto-assembling[J]. Biochimica et Biophysica Acta (BBA) - Proteins and Proteomics, 2014, 1844(2): 473-485.

    Berthelot K, Lecomte S, Estevez Y, et al. Hevea brasiliensis REF (Hev b 1) and SRPP (Hev b 3): An overview on rubber particle proteins[J]. Biochimie, 2014, 106: 1-9.

    Priya P, Venkatachalam P, Thulaseedharan A, et al. Differential expression pattern of rubber elongation factor (REF) mRNA transcripts from high and low yielding clones of rubber tree (Hevea brasiliensis Muell. Arg.) [J]. Plant Cell Reports, 2007, 26(10): 1833-1838.

    Oh S K, Kang H, Shin D H, et al. Isolation, characterization and functional analysis of a novel cDNA clone encoding a small rubber particle protein from Hevea brasiliensis[J]. Journal of Biological Chemistry, 1999, 274(24): 17132- 17138.

    Collinssilva J, Nural A T, Skaggs A, et al. Altered levels of the Taraxacum kok-saghyz (Russian dandelion) small rubber particle protein, TkSRPP3, result in qualitative and quantitative changes in rubber metabolism[J]. Phytochemistry, 2012, 79: 46-56.

    Hillebrand A, Post J, Wurbs D, et al. Down-regulation of small rubber particle protein expression affects integrity of rubber particles and rubber content in Taraxacum brevicorniculatum[J]. PLoS One, 2012, 7(7): e41874.

    Laibach N, Hillebrand A, Twyman R M, et al. Identification of a Taraxacum brevicorniculatum rubber elongation factor protein that is localized on rubber particles and promotes rubber biosynthesis[J]. The Plant Journal, 2015, 82(4): 609- 620.

    Chakrabarty R, Qu Y, Ro D K. Silencing the lettuce homologs of small rubber particle protein does not influence natural rubber biosynthesis in lettuce (Lactuca sativa)[J]. Phytochemistry, 2015, 113: 121-129.

    Arokiaraj P, Jones H, Cheong K F, et al. Gene insertion into Hevea brasiliensis[J]. Plant Cell Reports, 1994, 13(8): 425-431.

    Arokiaraj P, Yeang H Y, Cheong K F, et al. CaMV 35S promoter directs β-glucuronidas expression in the laticiferous system of transgenic Hevea brasiliensis (rubber tree)[J]. Plant Cell Reports, 1998, 17: 621-625.

    Montoro P, Rattana W, Pugade-Renaud V, et al. Production of Hevea brasiliensis transgenic embryogenic callus lines by Agrobacterium tumefaciens: roles of calcium[J]. Plant Cell Reports, 2003, 21: 1095-1102.

    Blanc G, Baptiste C, Oliver G, et al. Efficient Agrobacterium tumefaciens-mediated transformation of embryogenic calli and regeneration of Hevea brasiliensis Mull Arg. plants[J]. Plant Cell Reports, 2006, 24: 724-733.

    Montoro P, Lagier S, Baptiste C, et al. Expression of the HEV2.1 gene promoter in transgenic Hevea brasiliensis[J]. Plant Cell Tissue and Organ Culture, 2008, 94(1): 55-63.

    Leclercq J, Lardet L, Martin F, et al. The green fluorescent protein as an efficient selection marker for Agrobacterium tumefaciens-mediated transformation in Hevea brasiliensis (Mull. Arg)[J]. Plant Cell Reports, 2010, 29: 513-522.

    Leclercq J, Martin F, Sanier C, et al. Over-expression of a cytosolic isoform of the HbCuZnSOD gene in Hevea brasiliensis changes its response to a water deficit[J]. Plant Molecular Biology, 2012, 80: 255-272.

    Rekha K, Nazeem P A, Venkatachalam P, et al. Development of osmotin transgenics in Hevea brasiliensis Muell. Arg. using explants of zygotic origin[J]. Journal of Tropical Agriculture, 2014, 52(1): 7-20.

    Lestari R, Rio M, Martin F, et al. Overexpression of Hevea brasiliensis ethylene response factor HbERF‐IXc5 enhances growth and tolerance to abiotic stress and affects laticifer differentiation[J]. Plant Biotechnology Journal, 2018, 16(1): 322-336.

    Jayashree R, Nazeem P A, Rekha K, et al. Over-expression of 3-hydroxy-3- methylglutaryl-coenzyme A reductase 1 (hmgr1) gene under super-promoter for enhanced latex biosynthesis in rubber tree (Hevea brasiliensis Muell. Arg.)[J]. Plant Physiology and Biochemistry, 2018, 127: 414-424.

    Fan Y, Xin S, Dai X, et al. Efficient genome editing of rubber tree (Hevea brasiliensis) protoplasts using CRISPR/ Cas9 ribonucleoproteins[J]. Industrial Crops and Products, 2020, 146: 112146 .

    戴雪梅, 黃天帶, 李? 季, 等. 不同外植體對橡膠樹原生質(zhì)體分離和再生的影響[J]. 分子植物育種, 2014, 12(6): 1259-1264.

    Wang X, Shi M, Lu X, et al. A method for protein extraction from different subcellular fractions of laticifer latex in Hevea brasiliensis compatible with 2-DE and MS[J]. Proteome Science, 2010, 8: 35.

    Chow K, Ghazali A, Hoh C, et al. RNA sequencing read depth requirement for optimal transcriptome coverage in Hevea brasiliensis[J]. BMC Research Notes, 2014, 7(1): 69.

    Makita Y, Ng K K, Singham G V, et al. Large-scale collection of full-length cDNA and transcriptome analysis in Hevea brasiliensis[J]. DNA Research, 2017, 24(2): 159-167.

    Xia Z, Xu H, Zhai J, et al. RNA-Seq analysis and de novo transcriptome assembly of Hevea brasiliensis[J]. Plant Molecular Biology, 2011, 77(3): 299.

    Li D, Deng Z, Qin B, et al. De novo assembly and characterization of bark transcriptome using Illumina sequencing and development of EST-SSR markers in rubber tree (Hevea brasiliensis Muell. Arg.)[J]. BMC Genomics, 2012, 13(1): 192.

    Pirrello J, Leclercq J, Dessailly F, et al. Transcriptional and post-transcriptional regulation of the jasmonate signalling pathway in response to abiotic and harvesting stress in Hevea brasiliensis[J]. BMC Plant Biology, 2014, 14(1): 341.

    Fang Y, Mei H, Zhou B, et al. De novo transcriptome analysis reveals distinct defense mechanisms by young and mature leaves of Hevea brasiliensis (para rubber tree)[J]. Scientific Reports, 2016, 6: 33151.

    Wang X, Wang D, Sun Y, et al. Comprehensive proteomics analysis of laticifer latex reveals new insights into ethylene stimulation of natural rubber production[J]. Scientific Reports, 2015, 5: 13778.

    Xiang Q, Xia K, Dai L, et al. Proteome analysis of the large and the small rubber particles of Hevea brasiliensis using 2D-DIGE[J]. Plant Physiology and Biochemistry, 2012, 60: 207-213.

    Wang X, Shi M, Wang D, et al. Comparative proteomics of primary and secondary lutoids reveals that chitinase and glucanase play a crucial combined role in rubber particle aggregation in Hevea brasiliensis[J]. Journal of Proteome Research, 2013, 12(11): 5146-5159.

    Habib M A H, Gan C Y, Othman F, et al. Proteomics analysis of latex from Hevea brasiliensis (clone RRIM 600)[J]. Biochemistry and Cell Biology, 2017, 95(2): 232-242.

    Liu J, Zhuang Y, Guo X, et al. Molecular mechanism of ethylene stimulation of latex yield in rubber tree (Hevea brasiliensis) revealed by de novo sequencing and transcriptome analysis[J]. BMC Genomics, 2016, 17(1): 257.

    Li D, Wang X, Deng Z, et al. Transcriptome analyses reveal molecular mechanism underlying tapping panel dryness of rubber tree ( Hevea brasiliensis)[J]. Scientific Reports, 2016, 6: 23540.

    Chao J, Chen Y, Wu S, et al. Comparative transcriptome analysis of latex from rubber tree clone CATAS8-79 and PR107 reveals new cues for the regulation of latex regeneration and duration of latex flow[J]. BMC Plant Biology, 2015, 15: 104.

    Wei F, Luo S, Zheng Q, et al. Transcriptome sequencing and comparative analysis reveal long-term flowing mechanisms in Hevea brasiliensis latex[J]. Gene, 2015, 556(2): 153-162.

    Piyatrakul P, Yang M, Putranto R A, et al. Sequence and expression analyses of ethylene response factors highly expressed in latex cells from Hevea brasiliensis[J]. PLoS One, 2014, 9(6): e99367

    Nie Z, Kang G, Li Y, et al. Whole-transcriptome survey of the putative ATP-binding cassette (ABC) transporter family genes in the latex-producing laticifers of Hevea brasiliensis[J]. PLoS One, 2015, 10(1): e0116857.

    Shearman J R, Sangsrakru D, Jomchai N, et al. SNP identification from RNA sequencing and linkage map construction of rubber tree for anchoring the draft genome[J]. PLoS One, 2015, 10(4): e0121961.

    Makita Y, Kawashima M, Lau N S, et al. Construction of Pará rubber tree genome and multi-transcriptome database accelerates rubber researches[J]. BMC Genomics, 2018, 19(Suppl 1): 922.

    Ding Z, Fu L, Tan D, et al. An integrative transcriptomic and genomic analysis reveals novel insights into the hub genes and regulatory networks associated with rubber synthesis in H. brasiliensis[J]. Industrial Crops and Products, 2020, 153: 112562.

    猜你喜歡
    橡膠樹
    天然橡膠種植管理技術研究
    生如橡膠樹
    中國熱科院在橡膠樹低溫應答機制研究中取得重要進展
    橡膠樹開割季在5月已經(jīng)開始
    海南島橡膠樹風害坡度的影響特征
    不同懸空培養(yǎng)高度對橡膠樹小筒苗接穗生長的影響
    橡膠鏈格孢葉斑病病原菌生物學特性研究
    云南山地膠園橡膠樹專用肥肥效研究
    2013年海南省橡膠樹種苗基地調(diào)查報告
    4個不同橡膠樹品系的光合特性研究
    a级毛色黄片| 国产精品久久久久久精品古装| 亚洲av中文av极速乱| 99久久人妻综合| 中文精品一卡2卡3卡4更新| 免费黄网站久久成人精品| av又黄又爽大尺度在线免费看| 少妇高潮的动态图| 韩国av在线不卡| 精品一区二区免费观看| av女优亚洲男人天堂| 一级毛片 在线播放| 伦精品一区二区三区| 亚洲三级黄色毛片| 精品少妇久久久久久888优播| 网址你懂的国产日韩在线| 嘟嘟电影网在线观看| 欧美日韩国产mv在线观看视频 | 91久久精品国产一区二区三区| 国产91av在线免费观看| 日韩不卡一区二区三区视频在线| 在线观看免费视频网站a站| 日韩欧美精品免费久久| 看非洲黑人一级黄片| 国产精品不卡视频一区二区| 免费在线观看成人毛片| 男男h啪啪无遮挡| 久久久精品94久久精品| 中国国产av一级| 国产精品一及| 国产 一区精品| 久久久久久九九精品二区国产| 国产高潮美女av| 中文欧美无线码| 男女下面进入的视频免费午夜| 国产免费福利视频在线观看| 成人午夜精彩视频在线观看| 人妻系列 视频| 欧美zozozo另类| 草草在线视频免费看| 午夜福利高清视频| 亚洲性久久影院| 亚洲人成网站高清观看| 成年人午夜在线观看视频| 国国产精品蜜臀av免费| 久久久午夜欧美精品| 久久久久久久久久成人| 久久精品国产a三级三级三级| 一本久久精品| 各种免费的搞黄视频| 精华霜和精华液先用哪个| 久久精品熟女亚洲av麻豆精品| av国产久精品久网站免费入址| 免费看不卡的av| 亚洲av成人精品一二三区| 一级毛片黄色毛片免费观看视频| 少妇的逼水好多| 麻豆成人午夜福利视频| 日韩av免费高清视频| 哪个播放器可以免费观看大片| 亚洲三级黄色毛片| 国产片特级美女逼逼视频| 2018国产大陆天天弄谢| 婷婷色麻豆天堂久久| 一个人看视频在线观看www免费| 人妻少妇偷人精品九色| 国产精品三级大全| av免费观看日本| 在线观看av片永久免费下载| 国产乱人偷精品视频| 高清不卡的av网站| 99热网站在线观看| 一级毛片黄色毛片免费观看视频| 亚洲国产精品国产精品| 国产在视频线精品| 高清日韩中文字幕在线| 国产熟女欧美一区二区| 午夜免费男女啪啪视频观看| 国产在线男女| 成人高潮视频无遮挡免费网站| 亚洲成人av在线免费| 日本黄大片高清| 少妇人妻 视频| 激情五月婷婷亚洲| 欧美日韩国产mv在线观看视频 | 成人漫画全彩无遮挡| 两个人的视频大全免费| 少妇高潮的动态图| 欧美变态另类bdsm刘玥| 欧美成人午夜免费资源| 国内精品宾馆在线| 日韩av免费高清视频| 少妇人妻一区二区三区视频| 九九在线视频观看精品| 99久久精品一区二区三区| 亚洲va在线va天堂va国产| 久久久国产一区二区| 人人妻人人澡人人爽人人夜夜| 国产亚洲av片在线观看秒播厂| 国产亚洲精品久久久com| 国产综合精华液| 久久久久久伊人网av| 寂寞人妻少妇视频99o| 精品一区二区三卡| 啦啦啦啦在线视频资源| 国产精品久久久久成人av| 在线观看三级黄色| 久久久午夜欧美精品| 亚洲av二区三区四区| 精品视频人人做人人爽| 99热这里只有是精品在线观看| 久久精品久久精品一区二区三区| 日韩中文字幕视频在线看片 | 国精品久久久久久国模美| 久久久精品94久久精品| 国产真实伦视频高清在线观看| 国产精品一区二区三区四区免费观看| tube8黄色片| 狠狠精品人妻久久久久久综合| 乱系列少妇在线播放| 精品久久久精品久久久| 亚洲av福利一区| 寂寞人妻少妇视频99o| 大香蕉久久网| 亚洲av成人精品一区久久| 免费在线观看成人毛片| 有码 亚洲区| 搡女人真爽免费视频火全软件| 久久久久久久久久人人人人人人| 黄色欧美视频在线观看| 亚洲欧美中文字幕日韩二区| 在线观看美女被高潮喷水网站| 亚洲av男天堂| 亚洲av欧美aⅴ国产| 国产欧美亚洲国产| 精品亚洲成a人片在线观看 | 赤兔流量卡办理| 色视频www国产| 国产黄片视频在线免费观看| 日本猛色少妇xxxxx猛交久久| 国产 一区精品| 国产精品久久久久久久电影| 国产精品麻豆人妻色哟哟久久| 在线免费十八禁| 丰满迷人的少妇在线观看| 只有这里有精品99| 久久久久国产网址| 99热这里只有精品一区| 精品人妻熟女av久视频| 伦理电影免费视频| 国产精品福利在线免费观看| 国产精品一区二区在线观看99| 蜜桃在线观看..| 黄色欧美视频在线观看| 中文字幕亚洲精品专区| 亚洲av免费高清在线观看| 舔av片在线| 大片电影免费在线观看免费| 中文欧美无线码| 九九在线视频观看精品| 国产精品人妻久久久久久| 99热这里只有是精品50| 韩国高清视频一区二区三区| 最后的刺客免费高清国语| 久久影院123| 大片电影免费在线观看免费| 久久久久视频综合| 国产一区二区在线观看日韩| 纵有疾风起免费观看全集完整版| 看十八女毛片水多多多| 人妻制服诱惑在线中文字幕| 亚洲av中文字字幕乱码综合| 美女视频免费永久观看网站| 国产亚洲精品久久久com| 在线免费十八禁| 色视频www国产| 99久久人妻综合| 国产伦理片在线播放av一区| av.在线天堂| av在线蜜桃| www.av在线官网国产| 99热这里只有是精品50| 女的被弄到高潮叫床怎么办| av女优亚洲男人天堂| 午夜老司机福利剧场| 18禁在线无遮挡免费观看视频| 久久久久久久精品精品| 国产成人精品婷婷| 美女主播在线视频| 26uuu在线亚洲综合色| 男女边摸边吃奶| 99视频精品全部免费 在线| www.av在线官网国产| 18+在线观看网站| 国产亚洲av片在线观看秒播厂| 晚上一个人看的免费电影| 国产老妇伦熟女老妇高清| 精品久久久精品久久久| 毛片一级片免费看久久久久| av在线播放精品| 国国产精品蜜臀av免费| 亚洲美女视频黄频| 日本wwww免费看| 国产一区二区在线观看日韩| 精品国产露脸久久av麻豆| 女人久久www免费人成看片| 秋霞伦理黄片| 日韩欧美一区视频在线观看 | 婷婷色综合www| 国产av国产精品国产| 脱女人内裤的视频| 国产精品 欧美亚洲| www.精华液| 亚洲精品国产色婷婷电影| 精品少妇久久久久久888优播| 国产高清国产精品国产三级| 女人被躁到高潮嗷嗷叫费观| 亚洲精品第二区| 建设人人有责人人尽责人人享有的| 亚洲熟女毛片儿| 水蜜桃什么品种好| 国产精品久久久久成人av| 久久精品亚洲熟妇少妇任你| 午夜两性在线视频| 亚洲精品在线美女| 又紧又爽又黄一区二区| 亚洲自偷自拍图片 自拍| 久久精品亚洲熟妇少妇任你| 欧美国产精品va在线观看不卡| www日本在线高清视频| 一本久久精品| 91九色精品人成在线观看| 美国免费a级毛片| 老司机影院成人| 99久久99久久久精品蜜桃| 亚洲成人手机| 午夜日韩欧美国产| 国产精品久久久久成人av| 美女国产高潮福利片在线看| 老司机影院成人| 无限看片的www在线观看| 美女福利国产在线| 一级毛片我不卡| 美女国产高潮福利片在线看| 777久久人妻少妇嫩草av网站| 欧美日韩精品网址| 欧美性长视频在线观看| 国产精品一区二区精品视频观看| 亚洲熟女毛片儿| 国产亚洲一区二区精品| 久热爱精品视频在线9| a级毛片黄视频| 日韩中文字幕欧美一区二区 | 久久鲁丝午夜福利片| 国产欧美日韩精品亚洲av| 熟女av电影| 看免费成人av毛片| 久久精品人人爽人人爽视色| 中文字幕人妻丝袜制服| 在线天堂中文资源库| 亚洲人成网站在线观看播放| 1024香蕉在线观看| 精品人妻熟女毛片av久久网站| 国产亚洲欧美在线一区二区| 国产成人精品在线电影| 亚洲精品美女久久av网站| 国产在线免费精品| 欧美激情高清一区二区三区| 亚洲国产精品一区二区三区在线| 亚洲成人国产一区在线观看 | 久9热在线精品视频| 日韩一卡2卡3卡4卡2021年| 大码成人一级视频| 青春草亚洲视频在线观看| 亚洲精品久久午夜乱码| 人人妻人人添人人爽欧美一区卜| 久久久久久久精品精品| 我要看黄色一级片免费的| 超碰成人久久| 欧美日韩国产mv在线观看视频| 又黄又粗又硬又大视频| www日本在线高清视频| 老司机在亚洲福利影院| 一级黄色大片毛片| 国产片特级美女逼逼视频| 丰满饥渴人妻一区二区三| av天堂在线播放| 天天影视国产精品| 宅男免费午夜| www.av在线官网国产| 在线天堂中文资源库| 最近中文字幕2019免费版| 人成视频在线观看免费观看| 成人手机av| 国产成人一区二区三区免费视频网站 | 伊人亚洲综合成人网| 交换朋友夫妻互换小说| 欧美日韩成人在线一区二区| 麻豆国产av国片精品| 亚洲国产av新网站| 国产人伦9x9x在线观看| 色播在线永久视频| 在线 av 中文字幕| 一本大道久久a久久精品| 咕卡用的链子| 亚洲欧美日韩另类电影网站| 新久久久久国产一级毛片| 久久精品aⅴ一区二区三区四区| 岛国毛片在线播放| 欧美激情高清一区二区三区| 欧美激情高清一区二区三区| 一级,二级,三级黄色视频| 69精品国产乱码久久久| 成人三级做爰电影| 免费av中文字幕在线| netflix在线观看网站| 国产免费视频播放在线视频| 一本综合久久免费| 黄色 视频免费看| 青春草视频在线免费观看| 国产成人一区二区在线| 久久国产精品大桥未久av| 少妇裸体淫交视频免费看高清 | 女性被躁到高潮视频| 丰满饥渴人妻一区二区三| 国产一区有黄有色的免费视频| 亚洲成人国产一区在线观看 | 黑人猛操日本美女一级片| 免费一级毛片在线播放高清视频 | 亚洲一卡2卡3卡4卡5卡精品中文| 日韩一卡2卡3卡4卡2021年| 天堂8中文在线网| 曰老女人黄片| 男人舔女人的私密视频| av欧美777| 看十八女毛片水多多多| 国产黄色免费在线视频| 亚洲成国产人片在线观看| xxxhd国产人妻xxx| 欧美 亚洲 国产 日韩一| 国产片特级美女逼逼视频| 啦啦啦在线观看免费高清www| 欧美另类一区| 亚洲欧美日韩高清在线视频 | 伊人亚洲综合成人网| 欧美在线一区亚洲| 亚洲色图 男人天堂 中文字幕| 超色免费av| 日本一区二区免费在线视频| 亚洲国产av影院在线观看| 狂野欧美激情性xxxx| 国产精品一区二区免费欧美 | 国产成人啪精品午夜网站| 精品少妇内射三级| 久久精品久久久久久久性| 精品亚洲乱码少妇综合久久| 免费在线观看完整版高清| 国产精品久久久久久人妻精品电影 | 国产高清不卡午夜福利| 18在线观看网站| 只有这里有精品99| 欧美亚洲 丝袜 人妻 在线| 又粗又硬又长又爽又黄的视频| 看免费av毛片| 久久久久国产精品人妻一区二区| 欧美 日韩 精品 国产| 国产精品二区激情视频| 午夜福利免费观看在线| 黄频高清免费视频| 国产黄色免费在线视频| 日本wwww免费看| 欧美大码av| 亚洲熟女毛片儿| 波野结衣二区三区在线| 69精品国产乱码久久久| 91老司机精品| 性色av一级| 国产三级黄色录像| 欧美人与性动交α欧美精品济南到| 中国国产av一级| 飞空精品影院首页| 国产成人一区二区在线| 亚洲人成电影观看| 成人国产av品久久久| 久久人人97超碰香蕉20202| 韩国精品一区二区三区| 夫妻性生交免费视频一级片| 大码成人一级视频| 9191精品国产免费久久| 在线观看免费视频网站a站| 日本色播在线视频| 国产精品三级大全| 美女视频免费永久观看网站| 国产av精品麻豆| 精品久久久久久电影网| 久久久久久久精品精品| 不卡av一区二区三区| 免费日韩欧美在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲一码二码三码区别大吗| 日韩一卡2卡3卡4卡2021年| 亚洲欧美清纯卡通| 国产深夜福利视频在线观看| 一区二区日韩欧美中文字幕| 国产黄色免费在线视频| 国产精品久久久久成人av| 黄色视频在线播放观看不卡| 午夜免费成人在线视频| 欧美日韩综合久久久久久| 在线观看免费午夜福利视频| 校园人妻丝袜中文字幕| 国产淫语在线视频| 久久久亚洲精品成人影院| 在线观看免费高清a一片| 美女视频免费永久观看网站| 黄色视频不卡| 别揉我奶头~嗯~啊~动态视频 | 亚洲成av片中文字幕在线观看| 亚洲精品日韩在线中文字幕| 亚洲av男天堂| 国产免费视频播放在线视频| 久久久久久免费高清国产稀缺| 日韩制服丝袜自拍偷拍| 色94色欧美一区二区| 视频在线观看一区二区三区| av网站在线播放免费| 在线观看国产h片| 日韩 亚洲 欧美在线| 日本91视频免费播放| 国产精品二区激情视频| 国产有黄有色有爽视频| 亚洲色图 男人天堂 中文字幕| 可以免费在线观看a视频的电影网站| 久久99热这里只频精品6学生| 一级a爱视频在线免费观看| 精品人妻一区二区三区麻豆| 午夜福利,免费看| 欧美精品一区二区大全| 久久99精品国语久久久| 又粗又硬又长又爽又黄的视频| 亚洲色图综合在线观看| 狂野欧美激情性bbbbbb| videos熟女内射| 一本一本久久a久久精品综合妖精| 男女之事视频高清在线观看 | 最近最新中文字幕大全免费视频 | 亚洲精品国产av成人精品| 免费少妇av软件| 超色免费av| 一本综合久久免费| 国产精品.久久久| 日本vs欧美在线观看视频| 中文乱码字字幕精品一区二区三区| 国产在线观看jvid| 超色免费av| 一边摸一边做爽爽视频免费| kizo精华| 中文欧美无线码| 99久久综合免费| 99香蕉大伊视频| 少妇裸体淫交视频免费看高清 | 大香蕉久久成人网| 国产一区亚洲一区在线观看| 中文字幕av电影在线播放| 一本一本久久a久久精品综合妖精| 一区福利在线观看| 看十八女毛片水多多多| 国产精品国产三级专区第一集| 日韩制服骚丝袜av| 欧美中文综合在线视频| 久久热在线av| 欧美日韩亚洲国产一区二区在线观看 | 在线观看国产h片| 丰满少妇做爰视频| 午夜精品国产一区二区电影| 免费观看人在逋| 国产成人免费观看mmmm| 久久精品国产综合久久久| 999久久久国产精品视频| 欧美变态另类bdsm刘玥| 天天影视国产精品| 视频在线观看一区二区三区| svipshipincom国产片| 亚洲精品国产av成人精品| 丝瓜视频免费看黄片| 91麻豆精品激情在线观看国产 | 亚洲av国产av综合av卡| 日本色播在线视频| 精品第一国产精品| 秋霞在线观看毛片| 人妻一区二区av| 久久九九热精品免费| 国产精品 欧美亚洲| 免费在线观看视频国产中文字幕亚洲 | a 毛片基地| 美女中出高潮动态图| 精品国产一区二区三区久久久樱花| av网站在线播放免费| 一本大道久久a久久精品| 乱人伦中国视频| 成人亚洲欧美一区二区av| 欧美另类一区| a 毛片基地| 午夜福利一区二区在线看| 亚洲av日韩在线播放| 午夜免费鲁丝| 另类精品久久| 中文字幕人妻丝袜一区二区| 伊人亚洲综合成人网| 中文字幕最新亚洲高清| 黄色 视频免费看| 黄色视频不卡| 欧美成人精品欧美一级黄| 欧美变态另类bdsm刘玥| 国产黄色免费在线视频| 日韩av免费高清视频| 亚洲欧美日韩另类电影网站| 久久久久精品人妻al黑| 亚洲国产av影院在线观看| 91精品国产国语对白视频| 国产精品久久久久久精品古装| 国产伦理片在线播放av一区| 久久性视频一级片| 人人澡人人妻人| 黄色 视频免费看| 深夜精品福利| 人成视频在线观看免费观看| 国产亚洲欧美精品永久| 男人爽女人下面视频在线观看| 亚洲国产看品久久| 交换朋友夫妻互换小说| 久久久久精品国产欧美久久久 | 秋霞在线观看毛片| 美女午夜性视频免费| 国产成人精品无人区| 大片电影免费在线观看免费| 777久久人妻少妇嫩草av网站| 亚洲人成网站在线观看播放| 婷婷色麻豆天堂久久| 在线观看国产h片| 亚洲三区欧美一区| 国产精品二区激情视频| 免费在线观看黄色视频的| 国产1区2区3区精品| 18禁观看日本| av网站免费在线观看视频| 欧美日韩亚洲国产一区二区在线观看 | 人体艺术视频欧美日本| 一区二区av电影网| 日韩免费高清中文字幕av| 国产成人系列免费观看| 婷婷色麻豆天堂久久| 我要看黄色一级片免费的| 无遮挡黄片免费观看| 男女下面插进去视频免费观看| 91老司机精品| av天堂在线播放| 国产免费福利视频在线观看| √禁漫天堂资源中文www| 精品少妇一区二区三区视频日本电影| 久久av网站| 美女主播在线视频| netflix在线观看网站| 欧美av亚洲av综合av国产av| 高清黄色对白视频在线免费看| 国产男人的电影天堂91| 一级a爱视频在线免费观看| 亚洲男人天堂网一区| 亚洲欧美成人综合另类久久久| 亚洲色图综合在线观看| 老熟女久久久| 国产亚洲精品久久久久5区| 一级片'在线观看视频| 精品国产超薄肉色丝袜足j| 亚洲久久久国产精品| 人人妻人人澡人人看| 久久久久网色| 青春草视频在线免费观看| cao死你这个sao货| 丰满迷人的少妇在线观看| 999精品在线视频| 国产av精品麻豆| 男女国产视频网站| 好男人电影高清在线观看| 午夜av观看不卡| 91精品伊人久久大香线蕉| 美女午夜性视频免费| 操美女的视频在线观看| 欧美国产精品va在线观看不卡| 精品国产一区二区三区四区第35| 高清黄色对白视频在线免费看| 久久久久精品国产欧美久久久 | 免费日韩欧美在线观看| 久久精品国产亚洲av涩爱| 中国国产av一级| 丝袜脚勾引网站| 午夜免费观看性视频| 亚洲熟女精品中文字幕| 亚洲成人免费电影在线观看 | 18在线观看网站| 咕卡用的链子| 51午夜福利影视在线观看| 青青草视频在线视频观看| av在线老鸭窝| 亚洲欧洲日产国产| 午夜影院在线不卡| 欧美另类一区| 国产真人三级小视频在线观看| 亚洲欧洲日产国产| 嫩草影视91久久| 亚洲人成电影观看| 老司机影院成人| 亚洲中文av在线| 男女国产视频网站| 妹子高潮喷水视频|