橡膠樹橡膠生物合成調控相關基因表達豐度比較

楊署光 楊秀光 史敏晶 鄧小敏 晁金泉 李言 張世鑫 田維敏

摘?要：?蛋白質是構成生命系統的基本元件之一，是大部分生物學功能的執行者。蛋白質豐度與其生物學功能息息相關，其豐度受基因表達過程中各環節嚴格精密的調控。其中，蛋白質豐度與其相應mRNA豐度存在較強的相關性，蛋白質豐度差異的40%可由mRNA豐度來解釋。茉莉酸信號途徑調節巴西橡膠樹中的天然橡膠生物合成，但相關基因彼此間的表達豐度差異尚待闡明。該文比較了S/2D d3割膠制度下，15個橡膠生物合成調控相關基因COI1、JAZ1、JAZ2、JAZ3、MYC1、MYC2、MYC3、MYC4、MYC5、GAPDH、HMGR1、SRPP、REF、HRT1、HRT2以及2個常用內參基因18S、ACTIN1在10個橡膠樹種質膠乳中的表達豐度差異;將ACTIN1的表達豐度設定為1，以此為標準計算出樣品中其他基因的表達豐度。結果表明：相同個體中不同基因的轉錄豐度差異明顯，不同個體中相同基因集的豐度大小排序存在一定差異;同一基因在不同個體中的轉錄豐度差異明顯，這16個基因的最大豐度分別是最低豐度的9.43、6.04、10.02、12.29、18.82、9.22、38.46、112.83、121.36、15.34、19.09、13.54、10.05、19.80、24.83、11.82倍，他們的變異系數分別為73.05%、55.19%、69.09%、67.37%、66.59%、53.87%、83.25%、122.02%、166.34%、59.89%、70.59%、75.67%、74.20%、68.34%、84.23%、78.59%;總的來說，在群體水平上，16個基因的轉錄豐度從高到低依次為18S>SRPP>HMGR1>REF>MYC2/HRT1>COI1>MYC1/MYC4>GAPDH/JAZ1/MYC5>JAZ2>HRT2/MYC3/JAZ3，他們的群體平均豐度依次為ACTIN1的28 382.26、43.64、11.39、7.16、5.47、5.10、1.07、0.75、0.74、0.45、0.42、0.33、0.12、0.06、0.06、0.04倍。值得注意的是，無論在個體水平還是群體水平上，18S的豐度毫無疑問是最大的，在mRNA中，SRPP的豐度最大，JAZ1大于JAZ2和JAZ3，MYC2大于MYC1、MYC3、MYC4、MYC5，HRT1大于HRT2。綜上結果表明，結構基因和功能基因的豐度高于調控基因。在基因相對表達分析中，常對目的基因和內參基因作均一化處理，從而掩蓋了不同基因間的真實豐度差異，因此，在基因表達分析中，既要關注基因的相對表達量，也要關注基因間的豐度差異，這有助于更全面地理解基因的功能。

關鍵詞： 巴西橡膠樹， 橡膠生物合成調控， 基因豐度， 比較

Abstract：?Protein is one of the basic components of life system and the executor of most biological functions. Protein abundance is closely related to its biological function， and its abundance is strictly and precisely regulated by each link in the process of gene expression. Among them，?there is a strong correlation between protein abundance and its corresponding mRNA abundance， about 40% of the difference in protein abundance can be explained by mRNA abundance. Jasmonic acid signaling pathway regulates the biosynthesis of natural rubber in Hevea brasiliensis， but the difference of expression abundance among related genes needs to be elucidated. In the present study， the expression abundance differences of 15 rubber biosynthesis regulatory genes COI1， JAZ1， JAZ2， JAZ3， MYC1， MYC2， MYC3， MYC4， MYC5， GAPDH， HMGR1， SRPP， REF， HRT1， HRT2， and 2 common internal reference genes 18S， ACTIN1 in 10 rubber tree germplasms latex following tapping them with S/2D d3 tapping system were compared. The expression abundance of ACTIN1 in each sample is set to 1， and the expression abundance of other genes in the sample is calculated according to the standard. The results were as follows： The transcriptional abundance of different genes in the same individual was significantly different， and the abundance order of the same gene set was different in different individuals; The transcription abundance of the same gene was significantly different in different individuals， the maximum abundance of the 16 genes were 9.43， 6.04， 10.02， 12.29， 18.82， 9.22， 38.46， 112.83， 121.36， 15.34， 19.09， 13.54， 10.05， 19.80， 24.83， 11.82 times of the lowest abundance， and the coefficient of variation were 73.05%， 55.19%， 69.09%， 67.37%， 66.59%， 53.87%， 83.25%， 122.02%， 166.34%， 59.89%， 70.59%， 75.67%， 74.20%， 68.34%， 84.23%， 78.59%， respectively; Overall， at the population level， the transcription abundance of the 16 genes from high to low was 18S>SRPP>HMGR1>REF>MYC2/HRT1>COI1>MYC1/MYC4>GAPDH/JAZ1/MYC5>JAZ2>HRT2/MYC3/JAZ3， correspondingly， the average abundance were 28 382.26， 43.64， 11.39， 7.16， 5.47， 5.10， 1.07， 0.75， 0.74， 0.45， 0.42， 0.33， 0.12， 0.06， 0.06， 0.04 times than that of ACTIN1， respectively. It is worth noting that， the abundance of 18S is undoubtedly the highest， and in mRNA， SRPP is the largest， JAZ1 is greater than that of JAZ2 and JAZ3， MYC2 is greater than that of MYC1， MYC3， MYC4 and MYC5， HRT1 is greater than HRT2 at both the individual and population levels. The results showed that， the abundance of structural genes and functional genes is higher than that of regulatory genes. In the analysis of gene relative expression， the target gene and the internal reference gene are usually homogenized， thus masking the real abundance difference between different genes， therefore， in the gene expression analysis， we should pay attention not only to the relative expression of genes， but also to the abundance difference between genes， which is helpful for understanding the function of genes in a more comprehensive way.

Key words： Hevea brasiliensis， rubber biosynthesis regulation， gene abundance， comparison

蛋白質是構成生命系統的基本元件之一，是大部分生物學功能的執行者。蛋白質豐度與其生物學功能息息相關，其豐度受基因表達過程中各環節嚴格精密的調控。其中，mRNA豐度可以解釋蛋白質豐度差異的主要部分。

蛋白質豐度與其相應的mRNA豐度有一定的相關性。Lu et al.（2007）對大腸桿菌和酵母細胞蛋白質組進行定量時發現，細胞內蛋白質豐度與相應的mRNA豐度具有較高相關性（大腸桿菌R2=0.47，酵母R2=0.73）。Laurent et al.（2010）研究了7個代表物種的蛋白質和mRNA豐度，發現蛋白質豐度與其mRNA豐度存在0.36～0.70的正相關性。Schwanhusser et al.（2011）在研究小鼠NIH3T3細胞內定量蛋白質組時發現蛋白質豐度與相應mRNA豐度存在較強的相關性（R2=0.41）。Marquerat et al. （2012）對裂殖酵母的增殖細胞和靜止細胞進行定量蛋白質組研究，發現細胞內蛋白質豐度與其相應mRNA豐度具有一定相關性（R2=0.55）。綜上所述，在細胞群體水平上，蛋白質豐度與其相應mRNA豐度存在較強的相關性，mRNA豐度可以解釋蛋白質豐度差異的主要部分（約40%）。

茉莉酸信號途徑調控巴西橡膠樹的橡膠生物合成（Deng et al.， 2018），橡膠產量與茉莉酸信號途徑關鍵環節基因和橡膠生物合成酶基因的表達正相關（楊署光等，2019a，b）。目前的qPCR基因表達分析主要關注基因間、樣本組織間或處理條件下基因的相對表達量，很少關注功能相關基因的豐度差異。在基因相對表達分析中，常對目的基因和內參基因作均一化處理，從而掩蓋了不同基因間的真實豐度差異，因此，在基因表達分析中，既要關注基因的相對表達量，也要關注基因間的豐度差異，這有助于更全面地理解基因的功能。該研究將每個樣品中ACTIN1的基因豐度設為1，分析15個橡膠生物合成調控相關基因以及常用作內參基因的基因18S的豐度差異，以期能更全面地理解這些相關基因在橡膠生物合成調控中的地位和作用。

1?材料與方法

1.1 材料

實驗材料（表1）和實驗設計與先前的報道相同，先前報道了9個茉莉酸信號途徑關鍵環節基因HbCOI1、HbJAZ1、HbJAZ2、HbJAZ3、HbMYC1、HbMYC2、HbMYC3、HbMYC4、HbMYC5和6個橡膠生物合成酶基因HbHRT2、HbSRPP、HbREF、HbHMGR1、HbHRT1、HbGAPDH在這10個橡膠樹種質中的相對表達差異，以及這些基因間的表達相關性（楊署光等，2020）;該研究進一步利用這些基因的qPCR 結果，分析這些基因間的豐度差異。

1.2 方法

根據GenBank（JF775488）的全長cDNA序列設計ACTIN1的qPCR引物，F：GTTCTACAAGTGCTTTGATGGCGA，R：GCAGCCATAACATGAAACGCAATAG;引物的擴增效率為92.7%，qPCR產物為190 bp。

根據“A=2△Cq=2Actin1 Cq-Gene Cq”計算每個樣品中各基因的豐度（abundance，A），即將每個樣品中ACTIN1的豐度均設為1（A=2△Cq=2Actin1 Cq- Actin1 Cq=20=1）。

1.3 數據處理

采用GraphPad Prism 5作圖;采用SPSS軟件的Duncan檢驗進行多重比較分析。

2?結果與分析

2.1 基因轉錄的平均豐度

總的來說，在群體水平上，16個基因的轉錄豐度差異很大（圖1）。從高到低排序為18S>SRPP>HMGR1>REF>MYC2/HRT1>COI1>MYC1/MYC4>GAPDH/JAZ1/MYC5>JAZ2>HRT2/MYC3/JAZ3。其中MYC2/HRT1、MYC1/MYC4、GAPDH/JAZ1/MYC5、HRT2/MYC3/JAZ3之間差異不顯著（P>0.05）;SRPP基因的豐度極顯著（P<0.01）高于REF;HRT1基因的豐度極顯著（P<0.01）高于HRT2;5個MYCs家族成員中，MYC2>MYC1/MYC4>MYC5>MYC3（P<0.01）;3個JAZs家族成員中，JAZ1>JAZ2>JAZ3（P<0.01）。

2.2 相同樣品中不同基因的表達豐度

在個體水平上，相同個體中不同基因的轉錄豐度差異明顯（圖2-圖6）：如樣品23和28中，16個基因的豐度彼此間差異均到達P<0.01或P<0.05的顯著水平。不同樣品中相同基因集的豐度大小排序存在一定差異（圖2-圖6，表1）：所有樣品中，排名一致的是18S和SRPP，分別為第1和第2;排名第3的主要是HMGR1（67%）;排名第4的主要是REF（53%）;排名第5的主要是 MYC2（60%）;排名第6的主要是 HRT1（57%）;排名第7、第8的主要是 COI1（57%，40%）;排名第9的主要是 MYC4（30%）和GAPDH（30%）;排名第10、第11的主要是 JAZ1（47%，30%）;排名第12的主要是 MYC4（23%）、MYC5（17%）、GAPDH（17%）、JAZ2（17%）、MYC3（17%）;排名第13的主要是 JAZ2（50%）;排名第14、第15、第16的主要是 JAZ3（27%，30%，37%）;排名第16的主要是 JAZ3（37%）和MYC3（33%）。同一基因在不同樣品中的豐度排序范圍不一樣，如MYC2在30個樣品中有3種位次，而MYC5則多達9種位次。SRPP和REF是構成橡膠樹橡膠粒子的2種主要蛋白，本研究的所有樣品中，SRPP基因的豐度均極顯著（P<0.01）高于REF。HRT1和HRT2是橡膠樹橡膠轉移酶家族的2個成員，本研究的所有樣品中，HRT1基因的豐度均極顯著（P<0.01）高于HRT2。該研究的所有樣品中，MYCs家族的5個成員中，MYC2的豐度最高，MYC2基因的豐度均極顯著（P<0.01）高于MYC1、MYC3、MYC4、MYC5。該研究的所有樣品中，JAZs家族的3個成員中，JAZ1的豐度最高，JAZ1基因的豐度均極顯著（P<0.01）高于JAZ2，JAZ3;絕大部分樣品（25個，83.3%）的JAZ2豐度大于JAZ3，僅有5個樣品（7，8，9，15，21;16.7%）的JAZ3豐度大于（7，8，9;10%）或相當于（15，21;6.7%）JAZ2。

2.3 相同基因在不同樣品中的表達豐度

在個體水平上，同一基因在不同個體中的轉錄豐度差異明顯，變異系數在53.87%～166.34%之間（圖7-圖10）。COI1在3.9×10-1～3.7×100之間，變異系數73.05%，個體間差異可達1個數量級;JAZ1在1.9×10-1～1.15×100之間，變異系數55.19%，個體間差異可達1個數量級;JAZ2在3.0×10-2～3.0×10-1之間，變異系數69.09%，個體間差異可達1個數量級;JAZ3在1.0×10-2～1.3×10-1之間，變異系數67.37%，個體間差異可達1個數量級;MYC1在1.5×10-1～2.8×100之間，變異系數66.59%，個體間差異可達1個數量級;MYC2在1.4×100～1.24×101之間，變異系數53.87%，個體間差異可達1個數量級;MYC3在4.5×10-3～1.7×10-1之間，變異系數83.25%，個體間差異可達2個數量級; MYC4在4.0×10-2～4.4×100之間， 變異系數122.02%，個體間差異可達2個數量級;MYC5在2.0×10-2～4.4×100之間，變異系數166.34%，個體間差異可達2個數量級;GAPDH在8.0×10-2～1.2×100之間，變異系數59.89%，個體間差異可達2個數量級;HMGR1在1.5×100～2.9×101之間，變異系數70.59%，個體間差異可達1個數量級;SRPP在1.2×101～1.6×102之間，變異系數75.67%，個體間差異可達1個數量級;REF在2.8×100～2.8×101之間，變異系數74.20%，個體間差異可達1個數量級;HRT1在6.5×10-1～1.3×101之間，變異系數68.34%，個體間差異可達2個數量級;HRT2在1.1×10-2～2.8×10-1之間，變異系數84.23%，個體間差異可達1個數量級;18S在6.5×103～7.7×104之間，變異系數78.59%，個體間差異可達1個數量級。

3?討論與結論

核糖體（ribosome）是最古老、精細復雜的細胞器， 其結構和組成從原核到真核高度保守（靳聰聰等，2018）。核糖體RNA（ribosomal RNA，rRNA）是細胞中最為豐富的RNA，在真核細胞中約有50%的RNA是rRNA（蘇志寧和洪勵上，2009）;rRNA是核糖體的主要結構組分之一，占原核生物核糖體相對分子量的64.78%，占真核生物核糖體相對分子量的58.72%（劉望夷，2009）。本研究中，18S的基因豐度最高，與結構蛋白相對調控蛋白豐度較高（Sato et al.， 1999; Giegé & Brennicke， 1999; Lin et al.， 1999;王曦光等，2017），蛋白質表達水平較高的基因傾向于進化保守（Drummond et al.， 2005; Drummond & Wilke， 2008），與主要執行細胞核心功能（Beck et al.， 2011）的一般結果一致。

SRPP、HMGR1、REF、HRT1是橡膠生物合成的關鍵酶，SRPP和REF也是橡膠粒子的2個主要構成蛋白（Dennis & Light， 1989; Oh et al.， 1999; Berthelot et al.， 2014）;COI1、JAZ1、JAZ2、JAZ3、MYC1、MYC2、MYC3、MYC4、MYC5是茉莉酸信號途徑的關鍵蛋白，參與橡膠生物合成調控（Deng et al.， 2018）。本研究表明，SRPP、HMGR1、REF、HRT1的基因豐度顯著高于JAZs和MYCs家族成員，與功能蛋白/結構蛋白的豐度普遍高于調控蛋白（Ishihama et al.， 2008; Beck et al.， 2011; Nagaraj et al.， 2011）的結果相符。HMGR1位于多個代謝途徑的上游，細胞質中的GAPDH是參與糖酵解反應的關鍵酶，其基因豐度高于JAZ2、JAZ3和MYC3，與功能上參與基礎“物質流”的蛋白質豐度高于調控精細“信息流”的蛋白質豐度（Zhong et al.， 2012）的結果相似。

體外分析表明，HRT2而非HRT1與橡膠合成有關（Asawatreratanakul et al.， 2003），這與不同產量的橡膠樹種質中HRT1的基因表達比HRT2保守（楊署光等，2019a）的結果一致;然而，HRT1的基因豐度顯著高于HRT2，表明HRT1可能起著一種更為保守的功能。

基因組序列分析草圖表明，橡膠樹中分別存在10個REF和12個SRPP基因成員（Rahman et al.， 2013），這些基因在基因組中的定位相同（Oh et al.， 1999; Sookmark， 1999），進化分析表明REF和SRPP是同源蛋白，起源于共同的祖先基因，同屬于一個大的植物脅迫相關蛋白家族（Karine et al.， 2014）。REF和SRPP的mRNA在乳管細胞中高表達（Han et al.， 2000; Ko et al.， 2003; Chow et al.， 2007; Chotigeat et al.， 2010; Tan et al.， 2014），本研究獲得類似結果;在膠乳中，REFs家族的轉錄豐度最高（9.44%），其次是SRPPs家族（1.21%）（Chotigeat et al.， 2010），在個體成員水平上，本研究中SRPP的豐度大于REF。

MYC1、MYC2、MYC3在膠乳中特異表達，而MYC4、MYC5主要在花中表達（趙悅，2011），這和前者與橡膠產量的相關性高于后者（楊署光等，2019b）的結果一致;然而，膠乳中MYC4、MYC5的豐度還高于MYC3，MYC4的豐度與MYC1相當，說明MYC4、MYC5也在膠乳代謝過程中起作用。因此，在基因表達分析中，既要關注基因的相對表達量，也要關注基因間的豐度差異，這有助于更全面地理解基因的功能。

本研究結果表明，不同功能/環節的橡膠生物合成調控相關基因間的轉錄豐度差異顯著，并且這種差異在不同橡膠樹種質間基本一致;無論在種質內還是種質間，這些基因在個體內/間的轉錄水平均處于一定的變動狀態，相對而言，調控基因的變異程度大于功能基因和結構基因。本研究獲得的是割膠后3 d的基因轉錄豐度結果，由于割膠顯著促進橡膠樹的橡膠生物合成，并且割膠包含排膠和機械傷害2種效應，因此，研究傷害和排膠對相關基因轉錄豐度的影響以及系統的追蹤割膠后這些基因轉錄豐度的動態變化過程，有助于完善橡膠生物合成調控的理論機制。

參考文獻：

ASAWATRERATANAKUL K， ZHANG YW， WITITSUWAN?NAKUL D， et al.， 2003. Molecular cloning， expression and characterization of cDNA encoding cis-prenyltransferases from Hevea brasiliensis—A key factor participating in natural rubber biosynthesis[J]. Eur J Biochem， 270（23）： 4671-4680.

BECK M， SCHMIDT A， MALMSTROEM J， et al.， 2011. The quantitative proteome of a human cell line[J]. Mol Syst Biol， 7： 549.

BERTHELOT K， LECOMTE S， ESTEVEZ Y， et al.， 2014. Homologous Hevea brasiliensis REF （Hevb1） and SRPP （Hevb3） present different autoassembling[J]. Biochim Biophys Acta， 1844（2）： 473-485.

CHOTIGEAT W， DUANGCHU S， PHONGDARA A， 2010. cDNA library from the latex of Hevea brasiliensis[J]. Songklanakarin J Sci Technol， 32（6）： 555-559.

CHOW KS， WAN KL， ISA MN， et al.， 2007. Insights into rubber biosynthesis from transcriptome analysis of Hevea brasiliensis Latex[J]. J Exp Bot， 58（10）： 2429-2440.

DENG XM， GUO D， YANG SG， et al.， 2018. Jasmonate signalling in regulation of rubber biosynthesis in laticifer cells of rubber tree （Hevea brasiliensis Muell. Arg.）[J]. J Exp Bot， 69（15）： 3559-3571.

DENNIS MS， LIGHT DR， 1989. Rubber elongation factor from Hevea brasiliensis. Identification， characterization， and role in rubber biosynthesis[J]. J Biol Chem， 264（31）： 18608-18617.

DRUMMOND DA， BLOOM JD， ADAMI C， et al.， 2005. Why highly expressed proteins evolve slowly[J]. Proc Natl Acad Sci USA， 102（40）： 14338-14343.

DRUMMOND DA， WILKE CO， 2008. Mistranslation-induced protein misfolding as a dominant constraint on coding-sequence evolution[J]. Cell， 134（2）： 341-352.

GIEG P， BRENNICKE A， 1999. RNA editing in Arabidopsis mitochondria effects 441 C to U changes in ORFs[J]. Proc Natl Acad Sci USA， 96（26）： 15324-15329.

HAN KH， SHIN DH， YANG J， et al.， 2000. Genes expressed in the latex of Hevea brasiliensis[J]. Tree Physiol， 20（8）： 503-510.

ISHIHAMA Y， SCHMIDT T， RAPPSILBER J， et al.， 2008. Protein abundance profiling of the Escherichia coli cytosol[J]. BMC Genomics， 9： 102.

JIN CC， HOU MY， PAN YY， 2018. Research progress of ribosomal protein function in Arabidopsis thaliana[J]. J Plant Physiol， 54（2）： 203-212.[靳聰聰， 侯名語， 潘延云， 2018. 擬南芥核糖體蛋白生物學功能研究進展[J]. 植物生理學報， 54（2）： 203-212.]

KARINE B， SOPHIE L， YANNICK E， et al.， 2014. Hevea brasiliensis REF （Hev b 1） and SRPP （Hev b 3）： An overview on rubber particle proteins[J]. Biochimie， 106： 1-9.

KO JH， CHOW KS， HAN KH， 2003. Transcriptome analysis reveals novel features of the molecular events occurring in the laticifers of Hevea brasiliensis（para rubber tree）[J]. Plant Mol Biol， 53（4）： 479-492.

LAURENT JM， VOGEL C， KWON T， et al.， 2010. Protein abundances are more conserved than mRNA abundances across diverse taxa[J]. Proteomics， 10（23）： 4209-4212.

LIN XY， KAUL S， ROUNSLEY S， et al.， 1999. Sequence and analysis of chromosome 2 of the plant Arabidopsis thaliana[J]. Nature， 402（6763）： 761-768.

LIU WY， 2009. Structure and function of the bacterial ribosome[J]. Chin Bull Life Sci， 21（6）： 771-780.[劉望夷， 2009. 細菌核糖體的結構和功能[J]. 生命科學， 21（6）： 771-780.]

LU P， VOGEL C， WANG R， et al.， 2007. Absolute protein expression profiling estimates the relative contributions of transcriptional and translational regulation[J]. Nat Biotechnol， 25（1）： 117-124.

MARQUERAT S， SCHMIDT A， CODLIN S， et al.， 2012. Quantitative analysis of fission yeast transcriptomes and proteomes in proliferating and quiescent cells[J]. Cell， 151（3）： 671-683.

NAGARAJ N， WISNIEWSKI JR， GEIGER T， et al.， 2011. Deep proteome and transcriptome mapping of a human cancer cell line[J]. Mol Syst Biol， 7： 548.

OH SK， KANG H， SHIN DH， et al.， 1999. Isolation， characterization， and functional analysis of a novel cDNA clone encoding a small rubber particle protein from Hevea brasiliensis[J]. J Biol Chem， 274（24）： 17132-17138.

RAHMAN AY， USHARRAJ AO， MISRA BB， et al.， 2013. Draft genome sequence of the rubber tree Hevea brasiliensis[J]. BMC Genomics， 14： 75.

SATO S， NAKAMURA Y， KANEKO T， et al.， 1999. Complete structure of the chloroplast genome of Arabidopsis thaliana[J]. DNA Res， 6（5）： 283-290.

SCHWANHAUSSER B， BUSSE D， LI N， et al.， 2011. Global quantification of mammalian gene expression control[J]. Nature， 473（7347）： 337-342.

SOOKMARK U， PUJADE-RENAUD V， CHRESTIN H， et al.， 2002. Characterization of polypeptides accumulated in the latex cytosol of rubber trees affected by the tapping panel dryness syndrome[J]. Plant Cell Physiol， 43（11）： 1323-1333.

SU ZN， HONG LS， 2009. The transicription and regulation of ribosomal RNA[J]. J Langfang Teach Coll （Nat Sci Ed）， 9（4）： 74-77.[蘇志寧， 洪勵上， 2009. 核糖體RNA的轉錄與調控[J]. 廊坊師范學院學報（自然科學版）， 9（4）： 74-77.]

TAN DG， SUN XP， ZHANG JM， 2014. Age-dependent and jasmonic acid-induced laticifer cell differentiation in anther callus cultures of rubber tree[J]. Planta， 240（2）： 337-344.

TIAN WM， ZHANG H， YANG SG， et al.， 2013. Molecular and biochemical characterization of a cyanogenic β-glucosidase in the inner bark tissues of rubber tree （Hevea brasiliensis Muell. Arg.）[J]. J Plant Physiol， 170（8）： 723-730.

WANG XG， WANG J， ZHANG L， 2017. A. thaliana protein abundance analysis coresponding with elongation efficiency[J]. Chin Biotechnol， 37（2）： 40-47.[王曦光， 王娟， 張琳， 2017. 擬南芥蛋白質豐度與基因翻譯效率關聯分析[J]. 中國生物工程雜志， 37（2）： 40-47.]

YANG SG， CHEN YY， LI Y， et al.， 2019a. Correlation between the expression level of rubber biosynthesis genes and rubber yield[J]. Chin J Trop Crop， 40（3）： 475-482.[楊署光， 陳月異， 李言， 等， 2019a. 橡膠樹橡膠生物合成相關基因表達與橡膠產量的相關性[J]. 熱帶作物學報， 40（3）： 475-482.]

YANG SG， ZHAO Y， CHEN YY， et al.， 2019b. Correlation between the expression level of genes related to jasmonate signaling and rubber yield[J]. Guihaia， 39（5）： 641-649.[楊署光， 趙悅， 陳月異， 等， 2019b. 橡膠樹茉莉酸信號途徑相關基因表達與橡膠產量的相關性[J]. 廣西植物， 39（5）： 641-649.]

YANG SG， YANG XG， ZHAO Y， et al.， 2020. Correlation analysis on expression of genes related to rubber biosynthesis regulation of rubber tree[J]. Guihaia， 40（12）： 1790-1799.[楊署光， 楊秀光， 趙悅， 等， 2020. 橡膠樹橡膠生物合成調控相關基因的表達相關性分析[J].?廣西植物， 40（12）： 1790-1799.]

ZHAO Y， 2011. Involvement of jasmonate signaling pathway in regulating rubber biosynthesis in laticifer cells of Hevea brasilensis[D]. Haikou： Hainan University： 1-159.[趙悅， 2011. 巴西橡膠樹乳管細胞茉莉酸信號途徑對橡膠生物合成調節的研究[D]. 海口： 海南大學： 1-159.]

ZHONG F， YANG D， HAO YW， et al.， 2012. Regular patterns for proteome-wide distribution of protein abundance across species[J]. PLoS ONE， 7（3）： e32423.

（責任編輯?周翠鳴）