摘 要：【目的】研究不同徑級的大花序桉形成層的差異表達基因和表達模式，從基因表達水平揭示大花序桉徑向生長的相關(guān)機制，為林木遺傳育種和培育提供有益參考?！痉椒ā坎捎肐llumina HiSeq 2000測序平臺對5個不同徑級的大花序桉形成層的轉(zhuǎn)錄組進行測序，使用featureCount軟件計算基因表達水平FPKM，通過edgeR和TCseq進行基因差異表達分析和表達模式分析?！窘Y(jié)果】轉(zhuǎn)錄組測序經(jīng)質(zhì)控后獲高質(zhì)量reads共1 364 837 824條，平均GC含量51.07%，與參考基因組的比對率平均為68.49%。對不同徑級的大花序桉轉(zhuǎn)錄組進行差異表達基因分析，共篩選出8 794個差異表達基因，其中超高徑級組Vs中徑級組篩選出3 292個（上調(diào)表達1 878個，下調(diào)表達1 414個），大徑級組Vs中徑級組篩選出3 876個（上調(diào)表達2 314個，下調(diào)表達1 562個），小徑級組Vs中徑級組篩選出3 604個（上調(diào)表達1 599個，下調(diào)表達2 005個），超低徑級組Vs中徑級組篩選出4 565個（上調(diào)表達2 850個，下調(diào)表達1 715個）。差異表達基因的GO和KEGG分析顯示，差異表達基因主要富集在糖基轉(zhuǎn)移酶活性、有機環(huán)化合物結(jié)合、氧化還原酶活性等功能中，參與細胞色素P450代謝、苯丙烷生物合成、苯丙氨酸代謝、黃酮類生物合成等與次生木質(zhì)部形成相關(guān)的代謝途徑?；虮磉_模式分析將所有基因聚類為10個模式，其中模式9的基因隨徑級的降低表達量呈下降趨勢，主要富集在Hippo信號通路、萜類骨架生物合成等代謝通路。【結(jié)論】大花序桉徑向生長與多個代謝途徑相關(guān)，其中與木質(zhì)素合成相關(guān)的苯丙氨酸代謝途徑可能是較為重要的途徑，Hippo信號通路可能參與樹干大小的調(diào)控。

關(guān)鍵詞：大花序桉；轉(zhuǎn)錄組；徑向生長；差異表達基因；基因表達趨勢

中圖分類號：S722.3 文獻標志碼：A 文章編號：1673-923X（2024）10-0149-08

基金項目：廣西優(yōu)良用材林資源培育重點實驗室項目（2020-A-01-01）；廣西林業(yè)科技項目（桂林科研[2021]5號）。

Analysis of gene expression related to radial growth of Eucalyptus cloeziana

CHEN Shengkan1， LU Youwei2， DENG Ziyu1， LI Changrong1

（1. Key Laboratory of Central South Fast-growing Timber Cultivation of National Forestry and Grassland Administration， Guangxi Key Laboratory of Superior Timber Trees Resource Cultivation， Guangxi Forestry Research Institute， Nanning 530002， Guangxi， China； 2. Guangxi Gaofeng State-owned Forestry Farm， Nanning 530025， Guangxi， China）

Abstract：【Objective】The study was intended to study the gene expression profile of different diameter grade of Eucalyptus cloeziana to reveal the mechanism of radial growth at the gene expression level.【Method】We sequenced the cambium transcriptomes of five different diameter grades E. cloeziana using the Illumina HiSeq 2000 sequencing platform. Gene expression level was calculated using featureCount software. Difference expression and expression pattern analysis were performed by edgeR and TCseq.【Result】A total of 1 364 837 824 high-quality reads were obtained after quality control， with an average GC content of 51.07% and an average alignment rate of 68.49% to the reference genome. Difference expression analysis of all comparison groups screened 8 794 different expressed genes （DEGs）. 3 292， 3 876， 3 604 and 4 565 DEGs were screened in the comparison ultra-high diameter group， large diameter group， small diameter group， ultra-low diameter group Vs medium diameter group， respectively， and with 1 878 up-regulated and 1 414 down-regulated， 2 314 up-regulated and 1 562 down-regulated， 1 599 up-regulated and 2 005 down-regulated， 2 850 up-regulated and 1 715 down-regulated in these comparison， respectively. GO and KEGG analysis showed that DEGs were mainly enriched in glycosyltransferase activity， organocyclic compound binding， oxidorereductase activity， and were reported to be involved in cytochrome P450 metabolism， phenylpropane biosynthesis， phenylalanine metabolism， flavonoid biosynthesis and other metabolic pathways related to secondary xylem formation. Gene expression pattern analysis clustered all genes into 10 patterns， among which the genes in pattern 9 showed a downward trend with the decrease of diameter and were mainly enriched in Hippo signaling pathway and terpenoid backbone biosynthesis.【Conclusion】The radial growth of E. cloeziana is related to multiple metabolic pathways， among which the phenylalanine metabolic pathway related to lignin synthesis may be an important one， and the Hippo signaling pathway may be involved in the regulation of diameter size.

Keywords： Eucalyptus cloeziana； transcriptome； radial growth； different expression gene； gene expression trend

大花序桉Eucalyptus cloeziana是桃金娘科Myrtaceae桉屬Eucalyptus昆士蘭桉亞屬Idiogenes樹種，天然分布于澳大利亞昆士蘭州，具有4個不連續(xù)的地理分布區(qū)，其木材堅硬耐久、紋理通直、結(jié)構(gòu)均勻，材色呈黃褐色，具有黑金條紋，廣泛用于家具、建筑、礦柱和坑木等，是較好的鋸材樹種，在我國已作為大徑材進行培育[1-4]。大花序桉作為大經(jīng)級鋸材培育的周期在15～20 a，與尾葉桉E. urophylla、巨桉E. grandis及其雜種無性系等纖維材相比，其生長速度較慢，而生長量一直以來都是林木育種和培育的重要指標，包括高生長和徑向生長，林木徑向生長速度的提高能夠極大地促進林木蓄積量的增長，因此，解析林木徑向生長機制對林木育種和培育均具有重要意義。

對林木徑向生長的研究大多數(shù)集中在對氣候等環(huán)境因子的響應(yīng)，如賈存等[5]研究了華北落葉松Larix principis-rupprechtii人工林和天然林徑向生長對氣溫變化的響應(yīng)；楊保國等[6]研究了米老排Mytilaria laosensis徑向生長的變化規(guī)律及其對輻射總量、土壤含水量、相對空氣濕度、空氣溫度等環(huán)境因子的響應(yīng)；曹受金等[7]研究了氣候變化對不同緯度馬尾松Pinus massoniana徑向生長的影響。對林木徑向生長差異的分子調(diào)控機制的研究較少，主要在楊樹Populus[8]、巨桉[9]等用材樹種中開展，林木徑向生長的遺傳調(diào)控機制尚未清晰，仍需進一步研究。

林木的徑向生長（即次生木質(zhì)部生成）主要依賴于樹干側(cè)面分生組織維管束形成層細胞的分裂和分化。維管束形成層細胞向外分化形成次生韌皮部，向內(nèi)分化形成次生木質(zhì)部，形成層通過不斷分裂和分化增加木質(zhì)部的細胞數(shù)量，使林木徑向增粗[10-11]。次生木質(zhì)部的形成過程中，形成層向內(nèi)分裂形成的木質(zhì)部母細胞經(jīng)過初生細胞壁擴張而使細胞增大，隨后一系列纖維素和半纖維素等多糖沉積，次生細胞壁開始形成并增厚，同時開始木質(zhì)化。隨著木質(zhì)化次生細胞壁的成熟，細胞進入程序化死亡階段，最終完成木質(zhì)部細胞的分化，這一過程的各個階段均受到內(nèi)部遺傳信息的嚴格調(diào)控[12-13]。

以往對植物徑向生長發(fā)育及調(diào)控過程的研究發(fā)現(xiàn)，許多轉(zhuǎn)錄因子和植物激素在木質(zhì)部生成的各個階段均具有重要的調(diào)控作用。如形成層分裂及分化為木質(zhì)部母細胞的過程，受HD-ZIP III和WOX4等轉(zhuǎn)錄因子的調(diào)控[14-15]，楊樹中IAA3基因過表達會影響生長素信號，降低形成層細胞的分裂活動，抑制木材形成[16]；在木質(zhì)部細胞擴增和伸長過程中，兩種具有生物活性的赤霉素GA1和GA4主要聚集在擴增木質(zhì)部細胞區(qū)域參與擴增過程[17]；次生細胞壁的合成過程主要受纖維素合酶的影響，楊樹基因組中發(fā)現(xiàn)18個纖維素合酶基因[18]，其他一些蛋白如蔗糖酶也與纖維素合成有關(guān)[19]，此外，次生細胞壁的合成還受許多含NAC結(jié)構(gòu)域的轉(zhuǎn)錄因子、MYB轉(zhuǎn)錄因子以及赤霉素、水楊酸和油菜素內(nèi)酯的調(diào)控[20-22]。大量的研究表明，木本植物徑向生長的分子調(diào)控十分復(fù)雜，雖取得了較大的進展，但仍有許多未知亟待解析。

本研究對相同生長環(huán)境和林齡的大花序桉形成層進行轉(zhuǎn)錄組測序，開展不同徑級大花序桉形成層的差異表達基因和表達模式分析，挖掘大花序桉徑向生長相關(guān)的基因和代謝通路，為桉樹鋸材樹種的木材產(chǎn)量提高和新型優(yōu)良品種的培育提供理論基礎(chǔ)。

1 材料與方法

1.1 試驗材料

試驗點位于廣西壯族自治區(qū)玉林市容縣容西鎮(zhèn)（110°09′E，22°39′N），試驗材料為2004年5月份種植的大花序桉種源試驗林，試驗林采取統(tǒng)一的撫育措施。根據(jù)胸徑測定數(shù)據(jù)，設(shè)置5個不同的徑級，以總體均值μ為基礎(chǔ)徑級，大于和小于總體均值1倍和2倍標準差σ各設(shè)置1個徑級，每個徑級選擇6個單株取其形成層樣品，共30個樣品。取樣時間為2019年7月（大花序桉生長旺季）上午9：00—11：00，在胸徑處向南一側(cè)用小刀將5 cm×5 cm的樹皮割除，用解剖刀輕輕向內(nèi)刮取木質(zhì)部外的幾層水狀細胞，迅速將其用錫箔紙包好置于干冰中帶回實驗室提取RNA。

1.2 轉(zhuǎn)錄組測序及數(shù)據(jù)處理

RNA的提取利用EASY spin Plus植物RNA快速提取試劑盒（Aidlab）進行，濃度和純度采用NanoDrop2000（Thermo Scientific）和Qubit 3.0（Thermo Scientific）進行測定。采用Illumina TruseqTM RNA sample prep Kit（Illumina）試劑盒構(gòu)建文庫，在HiSeq 2000（Illumina）測序平臺上進行測序。使用Trimmomatic軟件對原始測序數(shù)據(jù)過濾，去除reads中的接頭序列，修剪掉3′端質(zhì)量值小于20的堿基，若剩余序列中仍含有質(zhì)量值小于10的堿基，則剔除整條序列，去除含N比率超過10%的reads，舍棄接頭及質(zhì)量修剪后長度小于75 bp的reads。使用Hisat2軟件將質(zhì)控后得到的高質(zhì)量reads與巨桉參考基因組（Eucalyptus grandis v2.0）比對，采用QualiMap軟件進行基因覆蓋度分析。

1.3 基因表達水平和表達模式分析

利用R程序包edgeR對不同徑級的樣本進行基因表達差異分析，差異基因的篩選標準為樣本組間差異倍數(shù)（Fold Change）≥2且錯誤發(fā)現(xiàn)率（FDR）<0.05，將得到的差異基因進行GO和 KEGG富集分析。通過TCseq程序包進行基因表達趨勢分析。

2 結(jié)果與分析

2.1 樣品胸徑顯著性檢驗

從5個不同的徑級中各選取6個單株，各組的胸徑平均值如表1所示，各組胸徑均值分別為36.0、31.0、27.2、23.3和18.4 cm，最大組和最小組之間胸徑差異較大。對各組胸徑進行方差分析，結(jié)果顯示不同經(jīng)級樣本的胸徑在0.001水平上差異顯著，多重比較顯示5個不同組別之間均存在顯著差異。大花序桉樣本徑向生長上的顯著差異表明對其進行形成層轉(zhuǎn)錄組測序分析可能獲得與徑向生長相關(guān)的差異表達基因。

2.2 轉(zhuǎn)錄組測序數(shù)據(jù)統(tǒng)計

對大花序桉群體各個徑級樣本形成層進行轉(zhuǎn)錄組測序，共獲得204 GB原始測序數(shù)據(jù)。經(jīng)質(zhì)控后，獲得高質(zhì)量reads共1 364 837 824條，平均每個樣品45 494 594條；獲得高質(zhì)量堿基數(shù)202 857 185 891 bp，GC含量為48.53%～53.03%，平均GC含量為51.07%，Q30%（質(zhì)量值≥30的堿基所占比例）為93.42%～95.44%，均值為94.48%。以巨桉基因組為參考基因組進行比對，比對率為59.05%～72.18%，平均比對率為68.49%，比對上的序列中86.42%～93.25%比對到基因組的外顯子區(qū)域。測序統(tǒng)計結(jié)果表明，轉(zhuǎn)錄組測序數(shù)據(jù)質(zhì)量良好，可用于下游分析。

2.3 不同徑級大花序桉差異基因分析

對所有樣品表達基因數(shù)量進行統(tǒng)計，可檢測到的表達基因數(shù)量平均為20 656個，F(xiàn)PKM大于10的基因數(shù)量平均為7 858個。以FC≥2且FDR<0.05為顯著性標準閾值，以中徑級組（μ組）為對照，對不同徑級的大花序桉樣品進行差異表達基因分析（圖1）。超高徑級組（μ+2σ）樣品相對中徑級組篩選出3 292個顯著差異表達基因，其中上調(diào)表達基因1 878個，下調(diào)表達基因1 414個；大徑級組（μ+σ）樣品相對中徑級組篩選出3 876個顯著差異表達基因，其中上調(diào)表達基因2 314個，下調(diào)表達基因1 562個；小徑級組（μ-σ）樣品相對中徑級組篩選出3 604個顯著差異表達基因，其中上調(diào)表達基因1 599個，下調(diào)表達基因2 005個；超低徑級組（μ-2σ）樣品相對中徑級組篩選出4 565個顯著差異表達基因，其中上調(diào)表達基因2 850個，下調(diào)表達基因1 715個。所有差異表達基因中，4組比較共有的差異表達基因共456個，超高徑級組和大徑級組上調(diào)表達基因與小徑級組和超低徑級組下調(diào)表達基因的共有基因共166個，超高徑級組和大徑級組下調(diào)表達基因與小徑級組和超低徑級組上調(diào)表達基因的共有基因共174個。

4個比較組共獲得8 794個差異表達基因，對差異表達基因進行GO功能注釋，主要集中在分子功能和生物過程兩大類，較少差異基因富集在細胞成分。GO富集分析顯示（圖2a），差異基因的主要功能富集于ADP結(jié)合、腺苷酸結(jié)合、核苷酸結(jié)合、碳水化合物衍生物結(jié)合、糖基轉(zhuǎn)移酶活性、有機環(huán)化合物結(jié)合、氧化還原酶活性等，與次生木質(zhì)部的形成密切相關(guān)，包括細胞的分裂分化以及纖維素和半纖維素等多糖沉積、木質(zhì)素單體合成和聚合等次生細胞壁生物合成的重要過程。差異基因的KEGG富集分析顯示（圖2b），差異基因顯著富集于細胞色素P450代謝、苯丙烷生物合成、苯丙氨酸代謝、黃酮類生物合成、半乳糖代謝等通路，其中苯丙氨酸代謝途徑在植物生長發(fā)育中具有重要作用，其下游主要有黃酮類化合物的合成分支途徑和木質(zhì)素的合成分支途徑。

2.4 不同徑級大花序桉形成層基因表達模式分析

對大花序桉形成層的表達基因進行表達模式分析，將相似表達模式的基因劃分聚類，結(jié)果顯示（圖3），所有表達的基因共被聚類為10種表達模式，其中模式1和模式10隨大花序桉徑級的降低基因的表達量呈上升趨勢，模式6和模式9隨徑級的降低基因的表達量呈下降趨勢，其余模式隨徑級的降低基因的表達量呈不規(guī)則趨勢。模式1和模式9隨徑級降低呈上升或者下降的趨勢較為明顯，所有表達基因中共2 417個基因被聚類到模式1，其中780個基因為差異表達基因，被聚類到模式9的基因共有3 060個，其中711個基因為差異表達基因。

分別對模式1和模式9的基因進行GO和KEGG富集分析，GO結(jié)果顯示，模式1中基因的主要功能富集在UDP-糖基轉(zhuǎn)移酶活性、有機物運輸和基于微管過程等，模式9中基因的主要功能富集在大分子生物合成過程、細胞酰胺代謝過程和非膜結(jié)構(gòu)細胞器等。KEGG富集分析結(jié)果顯示，模式1基因主要富集在縫隙鏈接、磷酸肌醇代謝、內(nèi)質(zhì)網(wǎng)蛋白質(zhì)加工等通路，模式9基因主要富集在Hippo信號通路和萜類骨架生物合成等代謝通路。

3 討 論

本研究比較分析了超高徑級、大徑級、中徑級、小徑級和超小徑級5個不同徑級大花序桉形成層的轉(zhuǎn)錄組，通過與中徑級轉(zhuǎn)錄組進行比較，共篩選出8 794個差異表達基因。進一步分析發(fā)現(xiàn)，差異基因主要功能富集在碳水化合物衍生物結(jié)合、糖基轉(zhuǎn)移酶活性、有機環(huán)化合物結(jié)合、氧化還原酶活性等，這些功能均與次生木質(zhì)部的形成密切相關(guān)。其中，糖基轉(zhuǎn)移酶在次生細胞壁的生物合成中較為重要，GO分析顯示基因表達模式分析中模式1基因的主要功能同樣富集在UDP-糖基轉(zhuǎn)移酶活性，糖基轉(zhuǎn)移酶的活性在大花序桉徑向生長中具有重要作用。糖基轉(zhuǎn)移酶廣泛存在于植物中，是專門催化糖基從活化的供體分子轉(zhuǎn)移到糖類或非糖類受體分子上形成糖苷的酶，其供體分子通常為尿苷二磷酸葡萄糖（UDP-glucose）[23-24]。較多研究表明糖基轉(zhuǎn)移酶參與木本植物次生細胞壁中碳水化合物的合成和架構(gòu)，直接影響木質(zhì)部的發(fā)育過程，如楊樹中已鑒定出25個與木材相關(guān)的糖基轉(zhuǎn)移酶[25]；若干糖基轉(zhuǎn)移酶家族成員參與葡糖木聚糖的生物合成[26]；通過對特定組織進行基因表達轉(zhuǎn)錄組學(xué)和生物信息學(xué)等分析，發(fā)現(xiàn)楊樹的纖維素合酶基因?qū)儆谔腔D(zhuǎn)移酶家族[27]。

木本植物的徑向生長與木質(zhì)素的生物合成密切相關(guān)，木質(zhì)素的合成和積累速度越快，徑向生長速度越快。木質(zhì)素是植物次生細胞壁所特有的一種苯丙烷類次生代謝物，主要由3種單體ρ-香豆醇（ρ-coumaryl alcohol）、松柏醇（coniferyl alcohol）和芥子醇（sinapyl alcohol）通過自由基偶聯(lián)聚合而成[28]。3種木質(zhì)素單體是苯丙氨酸代謝途徑的終產(chǎn)物，苯丙氨酸經(jīng)過脫氨反應(yīng)后，在芳香環(huán)上發(fā)生一系列的羥基化和甲基化修飾，側(cè)鏈發(fā)生氧化還原反應(yīng)，最終形成3種單體[29]。木質(zhì)素單體合成的一系列反應(yīng)由多個酶介導(dǎo)，包括一個解氨酶（PAL）、3個不同的細胞色素P450酶（C4H、C3′H和F5H）、2個甲基轉(zhuǎn)移酶（CCoAOMT和COMT）和2個氧化還原酶（CCR和CAD）等。陳沫等[30]挖掘和鑒定了不同種植密度下尾巨桉木質(zhì)部生成的關(guān)鍵基因，獲得了一批響應(yīng)種植密度的木質(zhì)部生成重要轉(zhuǎn)錄因子。本研究在對不同徑級差異基因的KEGG富集分析中發(fā)現(xiàn)，差異基因顯著富集于細胞色素P450代謝、苯丙烷生物合成、苯丙氨酸代謝、黃酮類生物合成等代謝途徑，與木質(zhì)素的生物合成密切相關(guān)，表明這些差異基因?qū)Υ蠡ㄐ蜩竦膹较蛏L有重要作用。

對林木徑向生長分子調(diào)控機制的研究多采用幼苗作為研究材料，如Du等[31]對銀腺楊“84K”P. alba×P. glandulosa幼苗6個不同發(fā)育階段的莖段，通過單細胞空間轉(zhuǎn)錄組研究其從初生生長到次生徑向生長的連續(xù)發(fā)育過程；Zhou等[9]對巨桉幼苗第3節(jié)間到第11節(jié)間莖段的發(fā)育梯度進行了轉(zhuǎn)錄組分析，解析了巨桉莖段從初生生長向次生生長轉(zhuǎn)變的差異表達基因，得到了可能對次生生長重要的基因和轉(zhuǎn)錄因子。然而，林木在主莖干上的徑向生長是林木蓄積增長的重要來源，主莖干上形成層的分裂分化活動可能與幼苗頂芽莖段存在差異。本研究以15年生的大花序桉為材料，對其胸徑處的形成層進行轉(zhuǎn)錄組分析，解析徑向生長相關(guān)的基因表達，對林業(yè)生產(chǎn)更具意義。

基因表達模式分析中，模式9隨徑級的降低基因的表達量呈下降趨勢，對模式9基因的KEGG富集分析顯示，模式9基因主要富集在Hippo信號通路、萜類骨架生物合成等代謝通路。其中Hippo信號通路是調(diào)控動物細胞分裂、器官大小的信號通路，植物中同樣存在與Hippo信號通路核心元件Hippo/MST1/2和MOB1/Mats同源的SIK1和MOB1蛋白[32]。植物Hippo信號通路的首要功能可能是促進細胞分裂退出與啟動核內(nèi)復(fù)制周期，從而調(diào)控植物器官大小。對擬南芥sik1突變體的研究發(fā)現(xiàn)，突變體各器官大小、成熟器官細胞數(shù)量、細胞大小、核內(nèi)復(fù)制水平均顯著小于野生型，葉肉細胞退出有絲分裂、啟動核內(nèi)復(fù)制的時機晚于野生型[33]。模式9中富集在Hippo信號通路的基因，其表達趨勢與擬南芥sik1突變體的研究結(jié)果相符，徑級即為器官大小，隨徑級的降低基因的表達量呈下降趨勢，模式9的基因與大花序桉徑向生長相關(guān)。

本研究以大花序桉為材料開展徑向生長相關(guān)的轉(zhuǎn)錄組研究，分析不同徑級形成層的差異表達基因和表達趨勢，挖掘徑向生長相關(guān)的基因和代謝通路，有助于理解桉樹徑向生長的遺傳調(diào)控機制，促進桉樹分子育種，提高桉樹鋸材樹種的木材產(chǎn)量和質(zhì)量，為新型優(yōu)良品種的培育提供理論基礎(chǔ)。然而，本研究只初步解析了徑向生長相關(guān)的差異表達基因及相關(guān)代謝通路，后續(xù)仍需對徑向生長起關(guān)鍵作用的基因或轉(zhuǎn)錄因子進一步挖掘、解析和驗證。

4 結(jié) 論

為了解大花序桉徑向生長的機制，本研究對不同徑級的大花序桉形成層轉(zhuǎn)錄組進行了測序，并開展了基因差異表達分析和表達模式分析。基因差異表達分析結(jié)果顯示，不同徑級的大花序桉存在大量的差異表達基因，主要富集于細胞色素P450代謝、苯丙烷生物合成、苯丙氨酸代謝和黃酮類生物合成等與次生木質(zhì)部密切相關(guān)的重要代謝途徑?；虮磉_模式分析結(jié)果顯示，所有表達基因聚類為10個不同的表達模式，其中模式1和模式9分別隨著徑級的降低表達量呈現(xiàn)上升和下降趨勢，模式1基因的主要功能富集在UDP-糖基轉(zhuǎn)移酶活性，模式9基因主要富集在Hippo信號通路、萜類骨架生物合成等代謝通路。本研究從轉(zhuǎn)錄組水平初步揭示了大花序桉徑向生長相關(guān)的差異表達基因及相關(guān)的代謝通路，有關(guān)結(jié)果對深入研究林木徑向生長機制，創(chuàng)制桉樹優(yōu)良新品種具有一定的參考價值。

參考文獻：

[1] PRY L D， JOHNSON L A S. A classification of the eucalyptus[M]. Australia： Australian National University，1971.

[2] KEATING W G， BOLZA E. Characteristics， properties and uses of timbers[M]. Australia： Northern Australia and the Pacific，1982.

[3] 王豁然.桉樹生物學(xué)概論[M].北京：科學(xué)出版社，2010. WANG H R. A Chinese appreciation of eucalyptus[M]. Beijing： Science Press，2010.

[4] 鄧紫宇，陳健波，郭東強，等.大花序桉的遺傳多樣性分析[J].林業(yè)科學(xué)研究，2019，32（4）：41-46. DENG Z Y， CHEN J B， GUO D Q， et al. Genetic diversity of Eucalyptus cloeziana[J]. Forest Research，2019，32（4）：41-46.

[5] 賈存，郭明明，王倩，等.華北落葉松人工林和天然林徑向生長對氣候變化的響應(yīng)[J].中南林業(yè)科技大學(xué)學(xué)報，2022，42（1）： 120-128. JIA C， GUO M M， WANG Q， et al. Response of the radial growth of Larix principis-rupprechtii plantations and natural forests to climate change[J]. Journal of Central South University of Forestry Technology，2022，42（1）：120-128.

[6] 楊保國，劉士玲，龐圣江，等.大青山米老排徑向生長變化及其對環(huán)境因子的響應(yīng)[J].中南林業(yè)科技大學(xué)學(xué)報，2020，40（2）： 29-35，53. YANG B G， LIU S L， PANG S J， et al. The stem radial growth of Mytilaria laosensis and its response to environmental factors in Daqing mountain[J]. Journal of Central South University of Forestry Technology，2022，40（2）：29-35，53.

[7] 曹受金，潘百紅.氣候變化對不同緯度馬尾松徑向生長的影響[J].中南林業(yè)科技大學(xué)學(xué)報，2015，35（12）：10-19. CAO S J， PAN B H. Effects of climate changes on Pinus massoniana growth in different latitudes in Nanling regions of south China[J]. Journal of Central South University of Forestry Technology，2015，35（12）：10-19.

[8] CHAO Q， GAO Z F， ZHANG D， et al. The developmental dynamics of the Populus stem transcriptome[J]. Plant Biotechnology Journal，2019，17：206-219.

[9] ZHOU F P， ZHANG H N， CHEN S S， et al. Transcriptome analysis of the transition from primary to secondary growth of vertical stem in Eucalyptus grandis[J]. BMC Plant Biology，2024，24：96.

[10] LARSON P R. The vascular cambium： development and structure[M]. Berlin： Springer-Verlag Berlin Heidelberg，1994： 9-30.

[11] RATHGEBER C B， CUNY H E， FONTI P. Biological basis of tree-ring formation： a crash course[J]. Frontiers in Plant Science， 2016，7：734.

[12] PLOMION C， LEPROVOST G， STOKES A. Wood formation in trees[J]. Plant Physiology，2001，127（4）：1513-1523.

[13] SUNDELL D， STREET N R， KUMAR M， et al. AspWood： highspatial-resolution transcriptome profiles reveal uncharacterized modularity of wood formation in Populus tremula[J]. The Plant Cell， 2017，29：1585-1604.

[14] CAMPBELL L， TURNER S. Regulation of vascular cell division[J]. Journal of Experimental Botany，2017，68：27-43.

[15] BRACKMANN K， QI J， GEBERT M， et al. Spatial specificity of auxin responses coordinates wood formation[J]. Nature Communications，2018，9：875.

[16] NILSSON J， KARLBERG A， ANTTI H， et al. Dissecting the molecular basis of the regulation of wood formation by auxin in hybrid aspen[J]. The Plant Cell，2008，20：843.

[17] ISRAELSSON M， SUNDBERG B， MORITZ T. Tissue-specific localization of gibberellins and expression of gibberellinbiosynthetic and signaling genes in wood-forming tissues in aspen[J]. The Plant Journal： for Cell and Molecular Biology， 2005，44：494-504.

[18] SUZUKI S， LI L， SUN Y H， et al. The cellulose synthase gene superfamily and biochemical functions of xylem-specific cellulose synthase-like genes in Populus trichocarpa[J]. Plant Physiology，2006，142：1233-1245.

[19] COLEMAN H D， YAN J， MANSFIELD S D. Sucrose synthase affects carbon partitioning to increase cellulose production and altered cell wall ultrastructure[J]. Proceedings of the National Academy of Sciences，2009，106：13118-13123.

[20] GUO H， WANG Y， LIU H， et al. Exogenous GA3 application enhances xylem development and induces the expression of secondary wall biosynthesis related genes in Betula platyphylla[J]. International Journal of Molecular Sciences，2015，16： 22960-22975.

[21] ZHONG R， YE Z H. Secondary cell walls： biosynthesis， patterned deposition and transcriptional regulation[J]. Plant Cell Physiology，2015，56：195-214.

[22] SANCHEZ-RODRIGUEZ C， KETELAAR K， SCHNEIDER R， et al. BRASSINOSTEROID INSENSITIVE2 negatively regulates cellulose synthesis in Arabidopsis by phosphorylating cellulose synthase 1[J]. Proceedings of the National Academy of Sciences of the United States of America，2017，114：3533-3538.

[23] LIM E K， BOWLES D J. A class of plant glycosyltransferases involved in cellular homeostasis[J]. The EMBO Journal，2004，23： 2915-2922.

[24] YIN R， HAN K， HELLER W， et al. Kaempferol 3-O-rhamnoside-7-O-rhamnoside is an endogenous flavonol inhibitor of polar auxin transport in Arabidopsis shoots[J]. New Phytologist， 2014，201：466-475.

[25] 楊少輝，王潔華，宋英今，等.楊樹中與次生細胞壁生物合成相關(guān)的糖基轉(zhuǎn)移酶研究進展[J].生物學(xué)雜志，2013，30（3）： 68-71. YANG S H， WANG J H， SONG Y J， et al. Involved in the secondary cell wall biosynthesis in popar[J]. Journal of Biology， 2013，30（3）：68-71.

[26] HIMMEL M E， DING S Y， JOHNSON D K， et al. Biomass recalcitrance： engineering plants and enzymes for biofuels production[J]. Science，2007，315：804-807.

[27] DJERBI S， ASPEBORG H， NILSSON P， et al. Identification and expression analysis of genes encoding putative cellulose synthases （CesA） in the hybrid aspen， Populus tremula （L.）×P. remuloides （Michx.）[J]. Cellulose，2004，11：301-312.

[28] YUKI T， CHEN F， JIN N， et al. Coexistence but independent biosynthesis of catechyl and guaiacyl/syringyl lignin polymers in seed coats[J]. The Plant Cell，2013，25（7）：2587-2600.

[29] MURO-VILLANUEVA F， MAO X Y， CHAPPLE C. Linking phenylpropanoid metabolism， lignin deposition， and plant growth inhibition[J]. Current Opinion in Biotechnology，2019，56： 202-208.

[30] 陳沫，何沙娥，陳少雄，等.基于轉(zhuǎn)錄水平解析尾巨桉徑向生長對種植密度的響應(yīng)[J].林業(yè)科學(xué)研究，2021，34（5）：1-12. CHEN M， HE S E， CHEN S X， et al. Response of stem radial growth of Eucalyptus urophylla×E. grandis to planting density based on transcriptome analysis[J]. Forest Research，2021，34（5）： 1-12.

[31] DU J， WANG Y C， CHEN W F， et al. High-resolution anatomical and spatial transcriptome analyses reveal two types of meristematic cell pools within the secondary vascular tissue of poplar stem[J]. Molecular Plant，2023，16（5）：809-828.

[32] 張平平，佟鑫，張?zhí)鞓?，?植物Hippo信號通路研究進展[J].遺傳，2017，39（7）：568-575. ZHANG P P， TONG X， ZHANG T L， et al. The emerging Hippo signaling pathway in plants[J]. Hereditas，2017，39（7）：568-575.

[33] XIONG J， CUI X， YUAN X， et al. The Hippo/STE20 homolog SIK1 interacts with MOB1 to regulate cell proliferation and cell expansion in Arabidopsis[J]. Journal of Experimental Botany， 2016，67（5）：1461-1475.

[本文編校：吳 彬]