橡膠樹橡膠生物合成調控相關基因的表達相關性分析

楊署光 楊秀光 趙悅 史敏晶 李言 鄧小敏 晁金泉 田維敏

摘?要：該文通過qPCR獲得了正常割膠條件下，9個茉莉酸信號途徑關鍵環節基因HbCOI1、HbJAZ1、HbJAZ2、HbJAZ3、HbMYC1、HbMYC2、HbMYC3、HbMYC4、HbMYC5和6個橡膠生物合成相關基因HbHRT2、HbSRPP、HbREF、HbHMGR1、HbHRT1、HbGAPDH在5個橡膠樹魏克漢種質和5個1981IRRDB種質膠乳中的表達數據;通過皮爾遜相關系數分析了這15個基因彼此間的表達相關性，分別獲得105對基因的雙變量相關系數（r12）和偏相關系數（r12.3），所有105對基因的雙變量相關系數（︱r12︱）和偏相關系數（︱r12.3︱）的平均值分別為0.486 ± 0.220和0.304 ± 0.211，達到0.05的差異顯著水平。其中，r12與r12.3方向相同的63對（60%），方向相反的42對（40%）;︱r12︱<︱r12.3︱的23對（21.905%），︱r12︱>︱r12.3︱的82對（78.095%）;雙變量相關顯著性P<0.05的76對（72.38%），P<0.01的59對（56.19%）;偏相關顯著性P<0.05的21對（20%），P<0.01的16對（15.24%）。結果顯示，這兩條途徑中的基因表達彼此相關，這為基因表達譜分析中普遍采用的前提假設“功能相關基因的表達相關”提供了進一步的實驗證據，可用于挖掘、篩選和預測與橡膠生物合成相關的未知基因，為研究橡膠樹產量形成的分子調控機理提供理論基礎。

關鍵詞：巴西橡膠樹， 割膠， 茉莉酸信號途徑， 橡膠生物合成， 基因表達， Pearson相關

中圖分類號：Q943.2

文獻標識碼：A

文章編號：1000-3142（2020）12-1790-10

Abstract：Tapping-enhanced rubber biosynthesis is closely related to the activation of jasmonic acid signaling in laticifer cells of rubber tree.The expression level of genes related to jasmonic acid signaling pathway and rubber biosynthesis were both positively correlated with the dry rubber yield.However， the exact relationship between the expression of genes involved in jasmonic acid signaling and rubber biosynthesis is not to be elucidated yet.In the present study， qPCR was used to analyze the expression of nine jasmonic acid signaling genes， HbCOI1， HbJAZ1， HbJAZ2， HbJAZ3， HbMYC1， HbMYC2， HbMYC3， HbMYC4， HbMYC5 and six rubber biosynthesis genes， HbHRT2， HbSRPP， HbREF， HbHMGR1， HbHRT1， HbGAPDH， in laticifer cells of five Wichham germplasms and five 1981IRRDB germplasms following tapping them with S/2D d3 tapping system， Hb18S was used as internal reference gene.The correlation between these 15 genes was analyzed through Pearson correlation coefficient.The bivariate correlation coefficient （r12） and partial correlation coefficient （r12.3） of 105 gene pairs were obtained with the mean ± standard deviation of 0.486 ± 0.220 and 0.304 ± 0.211， respectively.Among them， 63 gene pairs （60%） with r12 in the same direction as r12.3 and 42 genes pairs （40%） in the opposite direction; In the degree of correlation， 23 gene pairs （21.905%） with | r12 | less than | r12.3 |， and 82 gene pairs （78.095%） with | r12 | greater than | r12.3 |.There 76 gene pairs （72.38%） bivariate correlation coefficient （r12） were significant correlation at P<0.05 level and 59 gene pairs （56.19%） at P<0.01 level.In contrast， there less partial correlation coefficient （r12.3） were significant correlation at P<0.05 level （21pairs， 20%） and P<0.01 level （16 pairs， 15.24%）.These results suggested that the expression of genes involved in these two pathways were related to each other， which provides a further evidence for the assumption that “expression correlation of functional related genes”， which is widely adopted in gene expression profile analysis.It can be used for excavating， screening and predicting unknown genes related to rubber biosynthesis， as well as provides a theoretical basis for studying the molecular regulation mechanism of rubber yield formation in Hevea brasiliensis.

Key words：Hevea brasiliensis， tapping， jasmonate signaling， rubber biosynthesis， gene expression， Pearson correlation

茉莉酸是植物逆境響應的關鍵激素信號（Qi et al.， 2011），COI1、JAZ和MYC是茉莉酸信號途徑的關鍵環節（Chini et al.， 2009）;JA-Ile、COI1、JAZ和MYC通過精細的互作（Donnell et al.， 1996; Chini et al.， 2007; Fonseca et al.， 2009; Qi et al.， 2011）參與植物次生代謝的調節（Qi et al.， 2011; Schweizer et al.， 2013）。巴西橡膠樹中的天然橡膠生物合成是一種典型的植物類異戊二烯代謝，受茉莉酸信號途徑的調節（Deng et al.， 2018）;橡膠生物合成途徑包括一系列橡膠生物合成關鍵酶（Koyama & Tanaka， 1996），如REF（Dennis & Light， 1989）、HMGR（Chye et al.， 1992）、HRT（Asawatreratanakul et al.， 2003）、SRPP（Collins， 2009）等;茉莉酸信號途徑對這些酶的精細調節機制尚不完全清楚。橡膠樹種質PR107、RRIM600、熱墾628、熱墾525和熱墾523的干膠產量普遍高于RO/CM/10 44/160、MT/IT/13 29/8、RO/C/8 24/104、RO/I/103 107和RO/CM/10 44/454，并且茉莉酸信號途徑關鍵環節基因HbCOI1、HbJAZ1、HbJAZ2、HbJAZ3、HbMYC1、HbMYC2、HbMYC3、HbMYC4、HbMYC5的表達和橡膠生物合成酶基因HbHRT2、HbSRPP、HbREF、HbHMGR1、HbHRT1、HbGAPDH的表達與他們的橡膠產量正相關（楊署光等，2019a，b）。但是，這些基因表達的株間差異以及這些基因間的表達相關性尚不清楚。本研究通過分析這些基因在這10份種質中的表達相關性，可為研究橡膠樹產量形成的分子調控機理提供理論基礎。

1?材料與方法

1.1 材料

10份巴西橡膠樹（Hevea brasiliensis）種質、試劑、耗材與先前的報道相同（楊署光等，2019a，b）。其中，RO/C/8 24/104為2株、RO/I/103 107為4株，其他種質各3株，共30株樹。

1.2 方法

1.2.1 材料處理?生產中按S/2D d3正常割膠，在8月份的某次割膠時，分別收集前10 min流出的膠乳，用于提取膠乳總RNA。

1.2.2 總RNA的提取與cDNA的合成

用曾日中等（2003）的方法提取膠乳總RNA。cDNA第一鏈的合成與先前的報道相同（楊署光等，2019a，b）。

1.2.3 基因表達分析

1.2.3.1 qPCR引物合成?引物為序列來自文獻（Tian et al.， 2013;楊署光等，2019a，b），委托Invitrogen公司合成。

1.2.3.2 qPCR反應?qPCR反應按先前報道的方法完成（楊署光等，2019a，b）。

1.3 數據處理

以Hb18S作內參，按先前報道的方法（楊署光等，2019a，b）計算目的基因的相對表達值。用SPSS軟件分析基因表達值的Pearson 相關性，P<0.01表示極顯著相關，P<0.05表示顯著相關;相關性的密切程度參照文獻報道（朱婉麗，2019）的標準，根據相關系數（r）的數值大小，分為四個等級：0<|r︱≤0.3為微弱相關;0.3<︱r︱≤0.5為低度相關;0.5<︱r︱≤0.8為中度相關;0.8<︱r︱≤1為高度相關。

分別用SPSS軟件的Duncan檢驗和Excel TTEST（Array 1，Array 2，Tails 2，Type 2）完成多重比較分析和雙樣本比較分析：P<0.01表示組間差異極顯著，P<0.05表示組間差異顯著。用GraphPad Prism 5作圖。

2?結果與分析

2.1 橡膠樹種質膠乳中橡膠生物合成調控相關基因的表達分析

qPCR分析表明，不同種質中的基因表達差異明顯，相同種質的生物學重復間也能達到極顯著差異，并且高產種質中的基因表達普遍高于低產種質（圖1、圖2）。在高低產種質組間，COI1（圖1：A）、JAZ1（圖1：B）、JAZ2（圖1：C）、JAZ3（圖1：D）、MYC1（圖1：E）、MYC2（圖1：F）、MYC3（圖1：G）、MYC5（圖2：I）、GAPDH（圖2：J）、HMGR1（圖2：K）、REF（圖2：M）、HRT2（圖2：O）差異極顯著，HRT1（圖2：N）差異顯著，SRPP（圖2：L）差異不大，MYC4（圖1：H）差異不顯著;高低產之比分別為1.74、1.84、2.17、3.73、2.58、2.38、8.72、0.34、2.58、1.36、2.28、3.02、1.84、1.36、1.01。15個基因在高低產種質組間的平均表達差異極顯著，兩者之比為2.08（圖2：P）。

2.2 橡膠樹種質膠乳中橡膠生物合成調控相關基因表達的相關性分析

共分析了15個基因間的表達相關性，分別獲得105對基因的雙變量相關系數（r12）和偏相關系數（r12.3）（圖3、圖4）。其中，r12與r12.3方向相同的占60%，相反的占40%;微弱相關、低度相關、中度相關、高度相關的r12分別為20.952%、 33.333%、37.143%、8.571%，在r12.3中分別為54.286%、22.667%、18.095%、0.952%;︱r12︱<︱r12.3︱的為21.905%，︱r12︱>︱r12.3︱的為78.095%;顯著相關和極顯著相關的r12分別為72.38%和56.19%，在r12.3中分別為20%和15.24%;在30個植株中，9個茉莉酸信號途徑相關基因表達的平均值與6個橡膠生物合成相關基因表達的平均值間的r12為0.892，極顯著正相關。

分析了單個基因分別與其余14個基因間Pearson 相關系數的平均值（圖4：P）。所有105對基因的︱r12︱和︱r12.3︱的平均值分別為0.486和0.304，兩者差異顯著;HRT2、REF、COI1、JAZ1、JAZ3、MYC1、GAPDH、MYC2、MYC3的︱r12︱與︱r12.3︱差異顯著，而HMGR1、SRPP、HRT1、JAZ2、MYC4、MYC5的︱r12︱與︱r12.3︱差異不顯著;各組中︱r12︱普遍大于︱r12.3︱，但彼此間的差異顯著性具有廣泛的交集;5個MYC家族

成員中，MYC1，2，3與其他14個基因的相關性相當，且極顯著高于MYC4，5;3個JAZ家族成員中，JAZ1，2，3與其他14個基因的相關性差異不大，僅JAZ1，2間的︱r12︱達到0.05的顯著水平;相對而言，REF、GAPDH、MYC1、MYC2和MYC3與其他14個基因的相關性高于MYC4和MYC5。

3?討論與結論

橡膠樹的橡膠生物合成是一種典型的植物類異戊二烯代謝，茉莉酸信號途徑與類異戊二烯代因的表達與橡膠產量的相關性;同時，該研究結果顯示，這兩條途徑中的基因表達彼此間具有一定的相關性，說明功能相關的蛋白，基因表達傾向相關;這為基因表達譜分析中普遍采用的前提假設“功能相關基因的表達相關”提供了理論依據，可用于挖掘、篩選和預測與橡膠生物合成相關的未知基因，為研究橡膠樹產量形成的分子調控機理提供理論基礎。

從多個角度研究了“功能相關基因的表達相關”這一現象。在線蟲（Roy et al.， 2002; Lercher et al.， 2003）、酵母（Cohen et al.， 2000; 王晨光等，2006）、果蠅（Spellman & Rubin，2002）和人類（Lercher et al.， 2002;郭政等，2003）等不同物種上陸續發現了染色體上鄰近基因表達相關這一現象。酵母中互作蛋白質的編碼基因（王磊等，2004;王晨光等，2006;歐陽玉梅，2008）、蛋白質復合物的編碼基因（Jansen et al.， 2002;王磊等，2004）表達顯著相關，并且功能聯系越強基因表達相關性越高（王磊等，2004），相互作用蛋白質更傾向于對具有相同的亞細胞定位（王晨光等，2006;歐陽玉梅，2008）;對基因數據庫的分析表明，同一代謝通路上的基因表達傾向于高度相關（Harris et al.， 2004; Kanehisa et al.， 2004;Curtis et al.， 2005），在酵母（李傳星等，2004;王晨光等，2006）和狗（華琳等，2008）的研究中進一步證實了這種相關性。以上結果表明真核生物中的基因存在模塊化的共表達趨勢。

“功能相關基因的表達相關”具有實驗條件依賴性：不同的實驗條件對這種相關性有一定的影響（Bader et al.， 2004;王晨光等，2006）。橡膠生物合成途徑中也存在類似的實驗條件依賴性，如在該研究的S/2D d3的割膠制度下，MYC1和MYC2顯著正偏相關，而在S/2D d1的割膠制度下相反（Zhao et al.， 2011）。研究基因組內和跨基因組的基因表達譜與蛋白質相互作用之間的相關性，發現與人類、小鼠等物種相比，酵母的蛋白質相互作用與基因表達譜的相關性最弱（Bhardwaj & Lu， 2005），推測“功能相關基因的表達相關”具有進化程度依賴性。

“HbCOI1-HbJAZs-HbMYCs”彼此間、家族成員間以及他們與FPS、SRPP、REF、HRT2等橡膠生物合成關鍵酶的相互作用（馬弗明， 2010;趙悅，2011;劉偉，2011;何鑫， 2013;包杰，2014;肖華，2015;王靖等，2016;姚笛，2016;鄧小敏等， 2018;Deng et al， 2018）表明他們在調控橡膠生物合成過程中具有協同作用，為這些基因的表達相關性提供了理化基礎。

橡膠樹的產量性狀是一種典型的次生產量性狀，由反復合理的收獲脅迫（割膠）形成，是一個“形成-收獲-再形成-再收獲……”的過程（林位夫，2012）;2次割膠之間的橡膠再生以及相關基因的表達是一個動態過程;由于基因表達的動態性、瞬時性、多樣性和復合性，不同mRNA 的穩定性和降解速率不同，可能會影響特定實驗條件下的基因表達水平（Kruglyak & Tang， 2000）。因此，系統追蹤這些相關基因在排膠過程中的表達相關性是今后的一個研究方向。

參考文獻：

ASAWATRERATANAKUL K， ZHANG YW， WITITSUWANNAKUL 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.

BADER JS， CHAUDHURI A， ROTHBERG JM， et al.， 2004.Gaining confidence in high-throughput protein interaction networks[J].Nat Biotechnol， 22 （1）：78-85.

BAO J， 2014.Studies on the transcriptional regulation of key enzymes of natural rubber biosynthesis by MYC and Myb transcription factors[D].Haikou：Hainan University：1-70.[包杰， 2014.MYC和Myb轉錄因子對橡膠生物合成關鍵酶轉錄調節的研究[D].海口：海南大學：1-70.]

BHARDWAJ N， LU H， 2005.Correlation between gene expression profiles and protein-protein interactions within and across genomes[J].Bioinformatics， 21（11）：2730-2738.

CHINI A， BOTER M， SOLANO R， 2009.Plant oxylipins：COI1/JAZs/MYC2 as the core jasmonic acid-signalling module[J].Febs J， 276（17）：4682-4692.

CHINI A， FONSECA S， FERN NDEZ G， et al.， 2007.The JAZ family of repressors is the missing link in jasmonate signalling[J].Nature， 448（7154）：666-671.

CHYE ML， TAN CT， CHUA NH， 1992.Three genes encode 3-hydroxy-3-methylglutaryl-coenzyme A reductase in Hevea brasiliensis：hmgr1 and hmgr3 are differentially expressed[J].Plant Mol Biol， 19（3）：473-484.

COHEN BA， MITRA RD， HUGHES JD， et al.， 2000.A computational analysis of whole-genome expression data reveals chromosomal domains of gene expression[J].Nat Genet， 26（2）： 183-186.

COLLINS JE， 2009.The role of small rubber particle proteins in rubber biosynthesis[D].Reno， USA：University of Nevada：1-112.

CURTIS RK， ORESIC M， VIDAL-PUIG A， 2005.Pathways to the analysis of microarray data[J].Trends Biotechnol， 23（8）： 429-435.

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.

DENG XM， WANG J， CHEN DK， et al.， 2018.Preliminary screening of candidate interactors of HblJAZ3 from Hevea brasiliensis laticifer by using yeast two-hybridization methods[J].Chin J Trop Crop， 39（2）：287-292.[鄧小敏， 王靖， 陳多坤， 等， 2018.橡膠樹乳管細胞HblJAZ3互作蛋白的分離鑒定[J].熱帶作物學報， 39（2）：287-292.]

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.

DONNELL PJ， CALVERT C， ATZORN R， et al.， 1996.Ethylene as a signal mediating the wound response of tomato plants[J].Science， 274（5294）：1914-1917.

FONSECA S， CHINI A， HAMBERG M， et al.， 2009.（+）-7-iso-Jasmonoyl-L-isoleucine is the endogenous bioactive jasmonate[J].Nat Chem Biol， 5（5）：344-350.

GUO Z， ZHANG TW， LI X， et al.， 2003.Analysis of expression correlation within adjacent genes along the human chromosomes by GeneHub[J].Chin J Bioinform， 1（1）：11-14.[郭政， 張田文， 李霞， 等， 2003.利用GeneHub軟件分析人類染色體相鄰基因的表達相關性[J].生物信息學， 1（1）：11-14.]

HARRIS MA， CLARK J， IRELAND A， et al.， 2004.The gene ontology （GO） database and informatics resource[J].Nucl Acids Res， 32：D258-D261.

HE X， 2013.Research on the expression of several members of JAZ and MYC genes family correlated with natural rubber yield in Hevea brasilensis Muell.Arg.[D].Haikou：Hainan University：1-77.[何鑫， 2013.巴西橡膠樹JAZ和MYC家族幾個成員基因表達和產量相關性的研究[D].海口：海南大學：1-77.]

HUA L， ZHENG WY， LIU H， et al.， 2008.Relativity of gene expression and Co-regulated gene patterns in feature KEGG pathways[J].Chin J Biotechnol， 24（9）：1643-1648.[華琳， 鄭衛英， 劉紅， 等， 2008.特征代謝通路上的基因表達相關性及共調控表達模式[J].生物工程學報， 24（9）：1643-1648.]

JANSEN R， GREENBAUM D， GERSTEIN M， 2002.Relating whole-genome expression data with protein-protein interactions[J].Genome Res， 12（1）：37-46.

KANEHISA M， GOTO S， KAWASHIMA S， et al.， 2004.The KEGG resource for deciphering the genome[J].Nucl Acids Res， 32：D277-D280.

KOYAMA T， TANAKA Y， 1996.Isopentenyl diphosphate isomerase in rubber latex[J].Phytochemistry， 43（4）：769-772.

KRUGLYAK S， TANG H， 2000.Regulation of adjacent yeast genes[J].Trends Genet， 16 （3）：109-111.

LERCHER MJ， BLUMENTHAL T， HURST LD， 2003.Coexpression of neighboring genes in Caenorhabditis elegans is mostly due to operons and duplicate genes[J].Genome Res， 13（2）：238-243.

LERCHER MJ， URRUTIA AO， HURST LD， 2002.Clustering of housekeeping genes provides a unified model of gene order in the human genome[J].Nat Genet， 31（2）：180-183.

LI CX， LI X， GUO Z， et al.， 2004.Analysis of the mRNA expression similarity of genes in the same gene expression regulatory pathway[J].Heredias， 26 （6）：929-933.[李傳星， 李霞， 郭政， 等， 2004.調控通路內基因表達的相關性分析[J].遺傳， 26 （6）：929-933.]

LIN WF， 2012.Secondary yield characters of crops and their exploitation[J].Chin J Trop Crop， 33（4）：602-604.[林位夫， 2012.作物次生產量特性及其利用對策探討[J].熱帶作物學報， 33（4）：602-604.]

LIU W， 2011.Screen the proteins interacting with HbJAZ of Hevea brasiliensis[D].Haikou：Hainan University：1-55.[劉偉， 2011.橡膠樹HbJAZ1基因的互作蛋白篩選研究[D].海口：海南大學：1-55.]

LU SX， 2010.Structural and HblMYC1 gene correlated with natural rubber yield in Hevea brasilensis Muell.Arg.[D].Haikou：Hainan University：1-49.[盧世香， 2010.樹皮結構和HblMYC1基因與橡膠樹產量相關性研究[D].海口：海南大學：1-49.]

MA FM， 2010.Cloning and expression of transcription factor for cis-isopentenyl transferase gene from Hevea brasiliensis[D].Haikou：Hainan University：1-80.[馬弗明， 2010.橡膠樹順式-異戊烯基轉移酶基因轉錄因子的克隆及表達[D].海口：海南大學：1-80.]

OUYANG YM， 2008.Correlation between protein-protein interactions in yeast and gene expression profiles and subcellular localization[J].Chin J Biochem Mol Biol， 24（10）：950-956.[歐陽玉梅， 2008.酵母蛋白質相互作用與基因表達譜和亞細胞定位的相關性[J].中國生物化學與分子生物學報， 24（10）：950-956.]

QI T， SONG S， REN Q， et al.， 2011.The Jasmonate-ZIM-domain proteins interact with the WD-Repeat/bHLH/MYB complexes to regulate jasmonate-mediated anthocyanin accumulation and trichome initiation in Arabidopsis thaliana[J].Plant Cell， 23（5）：1795-1814.

ROY PJ， STUART JM， LUND J， et al.， 2002.Chromosomal clustering of muscle-expressed genes in caenorhabditis elegans[J].Nature， 418（6901）：975-979.

SCHWEIZER F， FERN NDEZ-CALVO P， ZANDER M， et al.， 2013.Arabidopsis basic helix-loop-helix transcription factors MYC2， MYC3， and MYC4 regulate glucosinolate biosynthesis， insect performance， and feeding behavior[J].Plant Cell， 25（8）：3117-3132.

SPELLMAN PT， RUBIN GM， 2002.Evidence for large domains of similarly expressed genes in the drosophila genome[J].J Biol， 1（1）：5.

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 CG， GUO Z， LI X， et al.， 2006.Multiple analysis of the relationship between correlation of gene function and expression by GeneHub software[J].J Biomed Eng， 23（6）：1157-1161.[王晨光， 郭政， 李霞， 等， 2006.利用 GeneHub軟件多角度分析基因功能相關與表達相關的聯系[J].生物醫學工程學雜志， 23（6）：1157-1161.]

WANG J， DENG XM， TIAN WM， 2016.Screening and in vitro functional analysis of the HblMYC3 interacted proteins in laticifer cell of Hevea brasiliensis[J].Chin J Trop Crop， 37（7）：1312-1317.[王靖， 鄧小敏， 田維敏， 2016.巴西橡膠樹乳管細胞HblMYC3互作蛋白篩選與功能分析[J].熱帶作物學報， 37（7）：1312-1317.]

WANG L， GUO Z， LI X， et al.， 2004.Comparative analysis of gene expression correlation between interactive proteins and among subunits of protein complexes[J].Chin J Bioinform， （4）：1-4.[王磊， 郭政， 李霞， 等， 2004.互作蛋白質與復合物亞基的基因表達相關性比較分析[J].生物信息學， 2（4）：1-4.]

XIAO H， 2015.Validation and characterization of HbJAZ1 interacted proteins of Hevea brasiliensis[D].Haikou：Hainan University：1-70.[肖華， 2015.巴西橡膠樹HbJAZ1互作蛋白的驗證及功能鑒定[D].海口：海南大學：1-70.]

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.]

YAO D， 2016.Cloning and interactions of HbJAZs of Hevea brasiliensis[D].Haikou：Hainan University：1-55.[姚笛， 2016.巴西橡膠樹HbJAZs基因的克隆與互作蛋白鑒定[D].海口：海南大學：1-55.]

ZENG RZ， DUAN CF， LI Y， et al.， 2003.Construction of cDNA library and analyses of cDNA sequences in latex from JA-stimulated rubber trees[J].Chin J Trop Crop， 24（3）：1-6.[曾日中， 段翠芳， 黎瑜， 等， 2003.茉莉酸刺激的橡膠樹膠乳cDNA消減文庫的構建及其序列分析[J].熱帶作物學報， 24（3）：1-6.]

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.]

ZHAO Y， ZHOU LM， CHEN YY， et al.， 2011.MYC genes with differential responses to tapping， mechanical wounding， ethrel and methyl jasmonate in laticifers of rubber tree （Hevea brasiliensis Muell.Arg.）[J].J Plant Physiol， 168（14）： 1649-1658.

ZHU WL， 2019.Study on application of correlation analysis principle in library and information analysis[J].Jiangsu Sci Technol Inform， （1）：9-12.[朱婉麗， 2019.相關性分析原理在圖書情報分析中的運用[J].江蘇科技信息， （1）：9-12.]

（責任編輯?李?莉）