摘要：【目的】通過多組學綜合分析白木通幼果在自花授粉后的脫落機制，為白木通良種繁育提供理論依據(jù)?！痉椒ā客ㄟ^人工授粉技術(shù)對白木通栽培種BMTGKAN1進行自花授粉和栽培種間授粉，授粉后第13和18d取樣，利用苯胺藍染色法分析授粉后花粉管的萌發(fā)情況，石蠟切片技術(shù)觀察授粉后胚珠發(fā)育情況，統(tǒng)計不同授粉處理的幼果留存率，并結(jié)合轉(zhuǎn)錄組和代謝組分析授粉后幼果脫落機制?！窘Y(jié)果】自花授粉和栽培種授粉的花粉在BMTGKAN1雌蕊內(nèi)萌發(fā)生長情況相似，但在受精過程中，栽培種間授粉胚珠在授粉后6h發(fā)生核融合，自花授粉組的胚珠面積較小，未發(fā)生核融合，授粉后18d栽培種間授粉的胚珠發(fā)育出胚乳和胚，而自花授粉胚珠未發(fā)育出胚乳和胚，出現(xiàn)內(nèi)外珠被分離，且此時幼果留存率存在明顯差異，即自花授粉后第18d幼果留存率僅為5.92%，栽培種間授粉后第18d幼果留存率高達86.24%。在轉(zhuǎn)錄和代謝水平上，栽培種間授粉與自花授粉相比，授粉后第13d和第18d差異表達基因（DEGs）分別為444和3646條，差異代謝物（DEMs）分別為427和412個。在相同授粉方式的比較中，自花授粉后第13d與第18d相比，有3305個DEGs和457個DEMs，栽培種間授粉第13d與第18d相比，有1123個DEGs和420個DEMs。實時熒光定量PCR檢測結(jié)果與轉(zhuǎn)錄組數(shù)據(jù)變化趨勢基本相同，說明轉(zhuǎn)錄組數(shù)據(jù)分析結(jié)果相對準確。在各對比組DEGs和DEMs聯(lián)合分析中，至少在一個對比組中被富集到的KEGG信號通路共有162條，其中有48條在各對比組中均被富集到。植物激素信號轉(zhuǎn)導通路、亞油酸代謝通路、ABC轉(zhuǎn)運蛋白通路和苯丙烷生物合成通路極有可能參與調(diào)控白木通自花授粉后的幼果脫落?！窘Y(jié)論】BMTGKAN1屬于晚期自交不親和植物，其自花授粉后13～18d大量脫落幼果，受其體內(nèi)植物激素信號轉(zhuǎn)導通路、亞油酸代謝通路、ABC轉(zhuǎn)運蛋白通路、苯丙烷生物合成通路調(diào)控，在種植實踐中選擇2個及以上花期相近的白木通栽培品種隔行交叉種植以提高坐果率。

關(guān)鍵詞：白木通；自花授粉；幼果脫落；轉(zhuǎn)錄組；代謝組

中圖分類號：S567.190.353文獻標志碼：A文章編號：2095-1191（2024）01-0137-14

Multi-omics comprehensive analysis of fruitlet abscission mecha-nism after self-pollination of Akebia trifoliata ssp.australis

LI Xian-hang，LONG Feng，XIONG Yi，CHEN Song-shu，LIU Hong-chang

（College of Agriculture，Guizhou University/Guizhou Key Laboratory of Propagation and Cultivation on Medicinal Plants，Guiyang，Guizhou550025，China）

Abstract：[Objective]In this study，multi-omics comprehensive analysis methods were used to explore the mecha-nism of fuitlet abscission after self-pollination of Akebia trifoliata ssp.anstralis，and to provide atheoretical basis for the breeding of improved varieties of A.trifoliata ssp.australis.【Method]Self-pollination and interspecific polination of A trifoliata ssp.australis cultivated variety BMTGKANI were carried out by artificial pollination techniques.Samples were taken on the13and18#d after pollination.The pollen tube germination was analyzed by aniline blue staining method，and ovule development was observed by paraffin section technique，and the retention rate of fruitlet under different polli-nation treatments was calculated.The mechanism offruitlet shedding after pollination was analyzed by transcriptome and metabolome.The germination and growth of pollen in BMTGKANI pisticles were similarbetween self-pollinated and cul-tivated varieties.However，during fertilization，the ovule of cultivated interspecific pollination occurred nuclear fusion at6h after pollination.The ovule area of self-polination group was small，and nuclear fusion did not occur.[Result]At18d after pollination，the ovules of inter-specific pollination developed endosperm and embryo，while the ovules of self-pollination developed no endosperm and embryo，and the inner and outer integuments were separated，and the fruitlet re-tention rate was quite different at this time.On the18*d after pollination，the fruitlet retention rate of self-pollinated was only5.92%，while that of cultivated species was as high as86.24%.At the level of transcription and metabolism，com-pared with self-pollination，there were444and3646differential expressed genes（DEGs）and427and412differentialme-tabolites（DEMs）on the13and18d after pollination.In the comparison of the same pollination methods，there were3305DEGs and457DEMs on the13d after self-pollination compared with the18#d after self-polination.There were1123DEGs and420DEMs in inter-species pollination.The results of real-time fluorescence quantitative PCR showed the same trend as the transcriptome data，indicating that the transcriptome data analysis results were relatively accurate.In the joint analysis ofDEGs and DEMs in each comparison group，a total of162KEGG metabolic pathways were enriched in at least one comparison group，of which48were enriched in each comparison group.Metabolic pathways such as planthor-mone signal transduction，linoleic acid metabolism，ABC transporters，and phenylpropanoid biosynthesis were likely to be involved in regulating the abscission of fruitlets after self-pollination of A.trifoliata ssp.australis.【Conclusion】BMTGKANI belongs to late self-incompatibility plant，its fruitlets of A.trifoliata ssp.australis will ll fall off within13-18d after self-pollination，and regulated by pathways of plant hormone signal transduction，linoleic acid metabolism，ABC transporters，phenylpropanoid biosynthesis.It is suggested that at least two varieties of A.trifoliata ssp.anstralis with similar flowering period should be selected for interlaced cross planting to improve fruit seting rate.

Key words：Akebia trifoliatassp.australis;self-pollination;fruitlet abscission;transcriptome;metabolome

Foundation items：National Key Research and Development Program of China（2021YFD1601002）;Guizhou Scien-ceand Technology Plan Project（QKHPTRC[2020]2101）

0引言

【研究意義】白木通（Akebia trifoliata ssp.austra-lis）是木通科（Lardizabalaceae）木通屬（Akebia）多年生木質(zhì)藤本植物，具有一定的觀賞價值，在歐洲一些國家常作為觀賞植物種植，其果實具有較高的藥用價值，可加工成為中藥材預知子（中國科學院中國植物志編輯委員會，2001;Liet al.，2010;國家藥典委員會，2020;Maciag et al.，2021）。隨著市場對白木通果實的需求和相關(guān)研究的推進，白木通已在廣西、貴州、四川、湖南等地進行人工種植，但野生狀態(tài)下木通屬植物的坐果率較低，通過人工異花授粉可大副提高其坐果率（熊大勝，1996;Zou et al.，2022）。目前對白木通自花授粉的親和性、自花授粉后影響白木通幼果脫落的關(guān)鍵原因以及自花授粉后白木通幼果中的轉(zhuǎn)錄和代謝活動變化情況均尚不明確，使得白木通無論是在良種繁育還是在高產(chǎn)栽培上的科學實踐較困難，也是導致目前栽培上品種混亂的原因之一。因此，研究白木通自花授粉后幼果脫落機制對了解其育性和開發(fā)利用具有重要意義。【前人研究進展】花器官發(fā)育不良和受精不良均會導致幼果脫落，受精不良主要指自交不親和（Mesejo et al.，2014;Garner and Lovatt，2016;Zhao et al.，2020）。自交不親和現(xiàn)象是開花植物在遺傳演化中為了阻止自花受精演變而來的一種精確的遺傳機制，可分為孢子體自交不親和、配子體自交不親和及晚期自交不親和三大類（Franklin-Tong and Franklin，2003;Goldberg et al.，2010;Waligórski and Szaleniec，2010;Zhou et al.，2020）。在晚期自交不親和植物中，自花授粉后花粉與柱頭能成功識別，且花粉管在花柱中能成功生長，胚珠抑制發(fā)生在合子前或合子后（Gibbs，2014）。胚珠抑制發(fā)生后進一步引發(fā)幼果脫落。晚期自交不親和為植物的遺傳進化提供幫助，但該遺傳機制會導致幼果大量脫落而降低坐果率，成為生產(chǎn)上制約作物產(chǎn)量的關(guān)鍵因素（Tatari et al.，2018）。在細胞學水平上，植物器官的脫落區(qū)細胞分化后獲得在脫落過程中響應環(huán)境刺激和生長發(fā)育信號的能力，受到脫落信號刺激時脫落區(qū)的細胞降解形成離層（Shi et al.，2023）。在授粉失敗、幼果受脅迫或果實成熟時，脫落區(qū)的附著能力就會減弱致使果實或幼果脫落（Nakano et al.，2013;Ito and Nakano，2015）。隨著現(xiàn)代組學技術(shù)發(fā)展，植物自花授粉后幼果脫落在轉(zhuǎn)錄水平和代謝水平上的調(diào)控機制逐步被揭示。有研究表明，植物激素在調(diào)控幼果脫落中發(fā)揮著重要作用，如乙烯、生長素、脫落酸、赤霉素等激素相互協(xié)調(diào)作用促進或延緩幼果脫落（Mahouachi et al.，2009;Nakano et al.，2014;Pathar-kar and Walker，2017;Mornya and Cheng，2018）。此外，亞油酸代謝、ABC轉(zhuǎn)運蛋白和苯丙烷生物合成等通路常在脫落幼果的轉(zhuǎn)錄組和代謝組富集分析中出現(xiàn)（Rees et al.，2009;Zhou et al.，2020;Qiu et al.，2021）。因此，這些通路在幼果脫落調(diào)控中發(fā)揮主要作用?！颈狙芯壳腥朦c】目前關(guān)于白木通的研究主要集中在藥理和資源開發(fā)，對其繁育特性研究較少，尤其是自花授粉后在細胞、轉(zhuǎn)錄和代謝層次的幼果脫落機制研究鮮見報道?！緮M解決的關(guān)鍵問題】通過人工授粉技術(shù)進行白木通栽培種BMTGKAN1自花授粉和栽培種間授粉，利用苯胺藍染色法分析授粉后花粉管的萌發(fā)情況，石蠟切片技術(shù)觀察授粉后胚珠發(fā)育情況，統(tǒng)計各授粉處理的幼果留存率，并結(jié)合轉(zhuǎn)錄組和代謝組分析授粉后幼果脫落機制，以期為白木通的良種繁育提供理論依據(jù)。

1材料與方法

1.1試驗材料

供試材料來源于貴州省開陽縣南江鄉(xiāng)的2個白木通栽培種BMTGKAN1和BMTGKAN2，經(jīng)貴州大學植物鑒定中心鑒定為白木通Akebia trifoliata（Thunb）Koidz.ssp.australis（Diels）Rehd。供試材料株齡5年，株行距1.5m×1.5m。選擇長勢、分支數(shù)和開花量相當?shù)闹仓曜鳛樵囼炛仓?。主要試劑：氯仿購自重慶川東化工（集團）有限公司；石蠟、番紅和固綠購自北京索菜寶科技有限公司；透明液購自南昌雨露實驗器材有限公司；RNAprep Pure多糖多酚植物總RNA提取試劑盒購自天根生化科技（北京）有限公司；乙醇、冰醋酸、甲醇（純度≥99.0%）、2-氯-L-苯丙氨酸（內(nèi)標標準物質(zhì)）（純度98%）、乙腈（純度≥99.9%）、甲酸（LC-MS grade）和甲酸銨（純度≥99.9%）購自天津市富宇精細化工有限公司；Bio-marker One Step SYBR Green RT-qPCR Kit試劑盒和無核酸酶水購自北京百邁客生物科技有限公司。主要儀器設備：BX61正置熒光相差顯微鏡（Olympus，日本）、LeiCA RM2016病理切片機[徠卡顯微系統(tǒng)（上海）貿(mào)易有限公司]、Pannoramic MIDI高分辨玻片掃描系統(tǒng)（上海乾思生物科技有限公司）、Nano-Drop2000/2000c超微量紫外分光光度計（Thermo-Fisher，美國）、DYY-6C電泳儀（北京六一生物科技有限公司）、GelDoc2000凝膠成像系統(tǒng)（Bio-Rad，美國）、NovaSeq6000測序儀（IIlumina，美國）、Vanquish液相色譜儀（ThermoFisher，美國）、ACQUITYUPLCHSST3色譜柱（2.1mm×150mm，1.8μm）（Waters，美國）、Q Exactive HF-X質(zhì)譜儀（ThermoFisher，美國）、SB-5200DT超聲波清洗器（寧波新芝生物科技股份有限公司）、QL-866混勻儀（海門市其林貝爾儀器制造有限公司）和qTower2.2型熒光定量PCR（Analy-tik Jena，德國）

1.2樣品處理及采集

試驗共設2個處理，即栽培種間授粉（IP）和自花授粉（SP）（表1）。在雄花散粉前2d將雄蕊剪除，并將幼果套袋，散粉時于9：00—11：00進行授粉，授粉后繼續(xù)套袋直至花期結(jié)束取下。在授粉后第13和18d每處理分別選取9株白木通植株，每3株4～6個幼果作為一個生物學重復，用IP1和IP2分別代表栽培種間授粉第13和18d樣品，SP1和SP2分別代表自花授粉第13和18d樣品。取樣后置于液氮中速凍，使用干冰轉(zhuǎn)運至蘇州帕諾米克生物醫(yī)藥科技有限公司進行轉(zhuǎn)錄組和代謝組相關(guān)檢測。幼果留存率統(tǒng)計方法：幼果留存率（%）=現(xiàn)存幼果數(shù)/幼果總數(shù)×100。

1.3苯胺藍染色方法

將采集的幼果用卡諾固定液（60%乙醇：30%氯仿：10%冰醋酸）固定24h后轉(zhuǎn)移至70%的酒精保存?zhèn)溆?。?mol/L的NaOH于65℃水浴軟化處理15～30min，再用0.1%苯胺藍溶液（苯胺藍0.05g，磷酸氫二鉀0.871g，雙蒸水50mL）染色。壓片后于正置熒光相差顯微鏡下觀測。

1.4石蠟切片方法

將SP和IP處理后6和18d的BMTGKAN1幼果用卡諾固定液（60%乙醇：30%氯仿：10%冰醋酸）固定24h后轉(zhuǎn)移至70%的酒精保存?zhèn)溆?。將試驗材料依次?5%、85%、95%、100%和100%乙醇梯度脫水。經(jīng)透明液透明后，依次經(jīng)過透明液和石蠟混合液（1：1）、純石蠟中進行滲蠟，滲蠟溫度65℃。然后將蠟塊于石蠟切片機中切片，切片厚度為8μm。將切片經(jīng)過酒精和透明液脫蠟后用番紅—固綠染色。將染色后載玻片于高分辨玻片掃描系統(tǒng)掃描后觀測。

1.5轉(zhuǎn)錄組測序及分析

使用RNAprep Pure多糖多酚植物總RNA提取試劑盒提取樣品總RNA，并用NanoDrop2000/2000c超微量紫外分光光度計和1%瓊脂糖凝膠電泳檢測樣品質(zhì)量。合格樣品經(jīng)過RNA抽提、純化、建庫之后，使用NovaSeq6000測序儀對文庫進行雙末端測序。

過濾原始下機數(shù)據(jù)[NGDC數(shù)據(jù)庫（https：//ngdc.cncb.ac.cn/），Bioproject為PRJCA012651]，去除平均質(zhì)量分數(shù)低于Q20的Reads，Cutadapt去除3'端帶接頭的序列，得到Clean data。使用HISAT2v2.0.5將配對末端Clean reads與白木通參照基因組（Biopro-ject為PRJCA012651）比對。使用FeatureCounts比對到每一條基因上Read Count值，采用TPM（每千個堿基的轉(zhuǎn)錄每百萬映射讀取的Transcripts）進行表達定量標準化，使用DESeq2進行表達分析，篩選出|log?Fold Change|gt;1且Plt;0.05的差異表達基因（DEGs）。

1.6實時熒光定量PCR檢測

挑選參與自花授粉后白木通幼果脫落代謝通路的關(guān)鍵基因，對其進行實時熒光定量PCR檢測，遵照實時熒光定量PCR的最小公開信息（Minimum information for publication of quantitative real-time PCR experiments，MIQE）指南進行試驗。使用Primer3Input（https：/primer3.ut.ee/）設計定量引物（表2），選用ACT7作為內(nèi)參基因（Vashisth et al.，2011）。所有基因引物序列由通用生物（安徽）股份有限公司合成。采用qTower2.2型熒光定量PCR儀進行PCR擴增，設3次重復。反應體系參照Bio-marker One Step SYBR Green RT-qPCR Kit說明書配制。擴增程序：42℃5min，95℃60s，進行40個循環(huán)；95℃5s，47℃32s。采用2-AAQ法計算目的基因的相對表達量，使用Rv4.3.0進行數(shù)據(jù)分析。

1.7代謝物檢測及分析

使用甲醇（含2-氯-L-苯丙氨酸）提取樣品，采用Vanquish液相色譜儀進行色譜檢測[ACQUITYUPLCHSST3色譜柱（2.1mm×150.0mm，1.8μm），流速0.25mL/min，柱溫40℃，進樣量2μL]，色譜條件參照Zelena等（2009）的方法。使用QExactiveHF-X質(zhì)譜檢測器、電噴霧離子源（ESI）及正負離子模式進行質(zhì)譜檢測，質(zhì)譜條件參照Want等（2013）的方法。

使用ProteoWizard對原始下機數(shù)據(jù)[NGDC數(shù)據(jù)庫（https：/ngdc.cncb.ac.cn/），Bioproject為PRJCA012651]進行處理，使用R XCMS并設定bw=2，ppm=15，peakwidth=c（5，30），mzwid=0.015，mzdiff=-0.01，method='centWave'對峰進行檢測、對齊、過濾處理。設定ppmlt;30，在HMDB、MassBank、LipidMaps、mzCloud、KEGG數(shù)據(jù)庫及蘇州帕諾米克生物醫(yī)藥科技有限公司自建物質(zhì)庫中進行物質(zhì)鑒定，并過濾掉QC樣本中相對標準偏差（RSD）gt;30%的物質(zhì)。根據(jù)統(tǒng)計檢驗計算P值、OPLS-DA降維方法計算變量投影重要度（VIP）和差異倍數(shù)（Fold Change），輔助篩選標志代謝物，當Plt;0.05和VIPgt;1時，判定代謝物具有統(tǒng)計學意義。采用MetaboAnalyst對篩選差異代謝分子進行功能通路富集。

2結(jié)果與分析

2.1授粉后花粉萌發(fā)和胚珠發(fā)育情況

SP組第13d的幼果留存率為21.62%，第18d為5.92%;IP組第13d的幼果留存率為98.72%，第18d為86.24%（圖1-A）。幼果內(nèi)花粉管生長試驗結(jié)果顯示，SP與IP組的花粉萌發(fā)和花粉管生長情況相似，均表現(xiàn)為授粉后3h萌發(fā)，12h聚集成簇向子房內(nèi)生長（圖2-B～圖2-G）。通過石蠟切片技術(shù)對授粉后胚珠進行觀察，結(jié)果發(fā)現(xiàn)IP組授粉后6h的胚珠極核附近被染紅面積較大，發(fā)生了核融合，而SP組的胚珠面積較小，未發(fā)生核融合，此后并未觀察到核融合現(xiàn)象（圖2-H和圖2-I）。授粉后18d，IP組的胚珠發(fā)育出胚乳和胚；而SP組的胚珠內(nèi)外珠被分離，未出現(xiàn)胚乳和胚，且中央細胞極核消失，表現(xiàn)出非正常雙子葉植物種子形成的特征（圖2-J和圖2-K）。綜上所述，SP和IP處理后花粉對白木通柱頭均具有親和性，均能產(chǎn)生正常生長的花粉管，但SP組的胚珠卻不能成功受精致使種子不能正常發(fā)育，表現(xiàn)為在授粉后第13和18d，IP組的幼果留存率均遠高于SP處理，最終SP處理的柱頭全部脫落，表明在BMTGKAN1中存在晚期自交不親和。

2.2轉(zhuǎn)錄組測序分析結(jié)果

去除在IP1 vs IP2與SP1vs SP2對比組中表達趨勢相同的DEGs，將其余4632個DEGs可分為四大類群（圖3），第I和第Ⅲ類DEGs在各樣本間表達情況差異最大，推測其與BMTGKAN1自花授粉導致

幼果脫落生理調(diào)控密切相關(guān)。使用皮爾遜相關(guān)系數(shù)表示樣品間基因的表達水平相關(guān)性，結(jié)果如圖4所示。相同生物學重復間具有極強線性相關(guān)性；授粉后13d SP與IP兩組相關(guān)性較高，授粉后第18d SP與IP組相關(guān)性較低，授粉后第13和18d2個時期相關(guān)性較低。同時，聚類結(jié)果表明SP2的轉(zhuǎn)錄情況異于SP1、IP1和IP2。

從SP與IP組授粉后第13d轉(zhuǎn)錄數(shù)據(jù)中共鑒定到444個DEGs，其中有127個下調(diào)表達，有317個上調(diào)表達；授粉后第18d轉(zhuǎn)錄組數(shù)據(jù)中共鑒定到3646個DEGs，其中2062個下調(diào)表達，1584個上調(diào)表達。從IP組授粉后第13和18d轉(zhuǎn)錄數(shù)據(jù)中共鑒定到DEGs有1123個，其中有493個上調(diào)表達，有630個下調(diào)表達；從SP組授粉后第13～18d轉(zhuǎn)錄數(shù)據(jù)中共鑒定到3305個DEGs，其中1247個上調(diào)表達，2058個下調(diào)表達（圖5），表明授粉后第18d時SP組基因表達與IP組存在明顯差異。SP組授粉后第13和18d，在所有分組中共有的DEGs為51個；在相同時期不同處理間共有的DEGs為220個，授粉后第13d IP組與SP組特有的DEGs為196個，授粉后第18d IP組與SP組特有的DEGs為3397個；在不同時期相同處理中，授粉后第13和18d時共有的DEGs為567個，特有的DEGs在SP組和IP組中分別為2739和557個（圖6）。

綜上所述，從對比組DEGs數(shù)量來看，授粉后第13d的DEGs數(shù)目遠少于授粉后第18d，IP組的DEGs基因數(shù)目遠少于SP組。這可能是在授粉后第13和18d時，SP組生理活動與IP組存在較大差異的原因，該結(jié)果與幼果在13～18d內(nèi)出現(xiàn)大量脫落的現(xiàn)象保持一致，表明試驗設計和取樣合理。

2.3代謝組分析結(jié)果

為揭示白木通自花授粉后幼果脫落機制，使用非靶向代謝技術(shù)鑒定了所有樣本中初生代謝物和次生代謝物。代謝模式相似的代謝物具有相似的功能，或共同參與同一代謝過程或細胞通路。本研究檢測到二級代謝物可聚類為三大類，其中第Ⅲ的表達量最高但在樣本間變化較少，第I和第Ⅱ類則相反（圖7）。由圖8可知，代謝檢測樣本在第一主成分（PC1）上表現(xiàn)出組內(nèi)樣本聚集，組間樣本分散，表明分析結(jié)果較可靠。本研究共檢測到958種代謝物（MS/MS），其中110種羧酸及其衍生物（Carboxylic acids and derivatives）、89種苯及其取代衍生物（Ben-zene and substituted derivatives）、75種脂肪?；‵atty acyls）、61種有機氧化合物（Organooxygen com-pounds）、48種甾體及其衍生物（Steroids and steroid derivatives）、42種異戊醇脂質(zhì)（Prenol lipids）、28種黃酮類（Flavonoids）、19種有機氧化物（Organic oxides）、18種氮雜環(huán)化合物（Azacyclic compounds）、18種酚類（Phenols）和其他類453種（圖9）。

IP與SP組在授粉后第13d時共有427種差異代謝物（DEMs），其中235種上調(diào)，192種下調(diào)；在授粉后第18d時，共有412種DEMs，其中166種上調(diào)，246種下調(diào)；IP組授粉后第13和18d間共有420種DEMs，其中254種上調(diào)，166種下調(diào)；SP組授粉后第13和18d間共有DEMs457種，其中221種上調(diào)，236種下調(diào)（圖10）?？傮w來說，各對比組中DEMs數(shù)目保持在同一水平。

2.4木通幼果脫落相關(guān)基因?qū)崟r熒光定量PCR驗證結(jié)果

為了驗證轉(zhuǎn)錄組測序數(shù)據(jù)可靠性，從轉(zhuǎn)錄組數(shù)據(jù)中選取9個ABC轉(zhuǎn)運蛋白代謝通路和激素信號轉(zhuǎn)導途徑關(guān)鍵基因（KEGG基因編號見表2）進行實時熒光定量PCR檢測，其中ABCBI、ABCC10和ABCC4屬于ABC轉(zhuǎn)運蛋白代謝通路調(diào)控基因；ARF18、AUX22、AUX28、GH3.11、GH3.6和IAA12屬于激素信號轉(zhuǎn)導代謝通路調(diào)控基因，ACT7作為內(nèi)參基因。由圖11可知，9個基因在不同時期和處理轉(zhuǎn)錄組數(shù)據(jù)與實時熒光定量PCR檢測結(jié)果變化趨勢基本一致，證明轉(zhuǎn)錄組數(shù)據(jù)分析結(jié)果相對準確。

對篩選出的DEGs和DEMs進行KEGG信號通路富集分析，結(jié)果顯示SP與IP組相同時期和不同時期DEGs共9010個，被富集到KEGG數(shù)據(jù)庫的新陳代謝（Metabolism）、遺傳信息處理（Genetic informa-tion processing）、環(huán)境信息處理（Environmental infor-mation processing）、細胞過程（Cellular processes）和生物系統(tǒng)（Organismal systems）5個分支（圖12-A）。對代謝組檢測到各對比組中DEMs進行KEGG信號通路富集分析，并挑選差異顯著性排名前10的通路進行分析，結(jié)果顯示在SP1vs SP2對比組中富集到植物激素信號轉(zhuǎn)導（Plant hormone signal trans-duction）、苯丙烷生物合成（Phenylpropanoid biosyn-thesis）等通路，IP1vs SP1和IP1vs IP2對比組中均富集到苯丙烷生物合成和ABC轉(zhuǎn)運蛋白（ABC transporters）等通路（圖12-B）。SP和IP組在不同時期差異主要是由授粉時間不同導致，相同時期差異主要是授粉方式不同所導致。

為進一步了解由白木通自花授粉引起幼果內(nèi)的轉(zhuǎn)錄和代謝變化機制，本研究將同一個比較組的DEGs和DEMs進行關(guān)聯(lián)分析，結(jié)果顯示有162條通路在轉(zhuǎn)錄組和代謝組中均被富集到，其中在4個比較組中均被富集到的代謝通路有48條（圖13）。結(jié)合Rees等（2009）、Zhou等（2020）、Qiu等（2021）和Shi等（2023）的研究報道進行分析，篩選出4條與自花授粉導致幼果脫落反應密切相關(guān)代謝通路，依次為植物激素信號轉(zhuǎn)導通路、亞油酸代謝通路、ATP結(jié)合盒轉(zhuǎn)運蛋白（ATP-binding cassette transporter pro-teins，ABC轉(zhuǎn)運蛋白）通路和苯丙烷生物合成通路，其中富集到植物激素信號轉(zhuǎn)導通路中的5種植物激素和77個相關(guān)基因，富集到苯丙烷生物合成通路中的19種代謝物和19個相關(guān)基因；富集到ABC轉(zhuǎn)運蛋白通路中的30種代謝物和12個相關(guān)基因。在亞油酸代謝途徑中共檢測到7種代謝物和8個相關(guān)基因（圖14）。

在植物激素信號轉(zhuǎn)導通路中，赤霉素A4（Gib-berellinA4）、（S）-脫落酸[（S）-abscisic acid]、順式玉米素（Cis-zeatin）、茉莉酸（Jasmonic acid）和水楊酸（Salicylic acid）5種植物激素參與調(diào)控白木通幼果脫落。在授粉后第13d至第18d過程中，除順式玉米素在IP1vs IP2對比組中相對表達量上調(diào)外，其余4種植物激素在SP1vs SP2和IP1vs IP2對比組中均下調(diào)。但在IP1vs SP1和IP2vs SP2對比組中表達趨勢出現(xiàn)差異，具體表現(xiàn)為水楊酸和赤霉素A4上調(diào)，但水楊酸不是DEM;順式玉米素下調(diào)；茉莉酸和（S）—脫落酸在IP1vs SP1對比組中上調(diào)，在IP2vs SP2對比組中下調(diào)。乙烯響應因子基因（ERFIB）在所有對比組中均上調(diào)表達，SP組的ERFIB基因表達量較IP組顯著上調(diào)。乙烯受體基因（ETR1、ETR2、CTRI、EIN2和EIN3）在SPI vs SP2對比組中均表達上調(diào)，在IP1vs IP2和IP1vs SP1對比組中均不是DEGs，但均上調(diào)表達，在IP2vs SP2對比組中除ETR2不是DEG外均上調(diào)表達。富集到生長素轉(zhuǎn)導相關(guān)通路的DEGs共37個，分別為ARF家族成員5個、GH3家族成員5個、AUX/IAA家族成員10個、AUX/LAXs家族成員2個、SAUR家族成員13個、TIR1和AFB2。其中，GH3家族5個DEGs雖然同屬一個家族，但在各對比組中表達趨勢各不相同，LAX家族成員LAX5與LAX2基因表達模式相反。茉莉酸調(diào)控基因TIFY9在SP1vs SP2對比組中下調(diào)表達，在其他對比組中差異不明顯。脫落酸調(diào)控基因AFB2在SP1vs SP2對比組中上調(diào)表達，在IP1vs IP2處理中表達呈上調(diào)趨勢，但AFB2基因在該對比組中差異未達到顯著水平（Pgt;0.05，下同）（圖15）。

亞油酸代謝途徑調(diào)控基因LOX5、LOX1.5和SDPI在IP2vs SP2和SP1vs SP2對比組中顯著上調(diào)表達，其他對比組中則差異不顯著;代謝物9-OxoODE、9，12，13-TriHOME和13-OxoODE在IP1vs IP2和SP1vs SP2對比組呈顯著上調(diào)，IP1vs SP1和IP2vs SP2對比組中顯著下調(diào)，而13S-hydroxyoctadecadienoicacid則相反（圖16-A）。ABC轉(zhuǎn)運蛋白途徑7個功能模塊共涉及30個代謝物，寡糖、多元醇和脂質(zhì)轉(zhuǎn)運蛋白模塊代謝物最多，共16個。磷酸鹽和氨基酸轉(zhuǎn)運蛋白模塊共有7個。ABC轉(zhuǎn)運蛋白途徑涉及4個亞族12個基因，分別為ABCB（6個）、ABCC（4個）、ABCC10（1個）和ABCG2（1個）（圖16-B）。從苯丙烷生物合成途徑中篩選到最多的基因是編碼過氧化物酶的基因，共19個，PER家族基因有15個（圖16-C）。

3討論

3.1晚期自交不親和或早期近交衰退導致白木通自花授粉后幼果脫落

白木通在自花授粉后幼果約在1周內(nèi)表現(xiàn)同栽培種間授粉，而且自花授粉和栽培種間授粉的花粉管萌發(fā)現(xiàn)象一致，但在授粉后第13～18d時幼果大量脫落。本研究發(fā)現(xiàn)，IP組胚珠在授粉后6h發(fā)生核融合，SP組則無該現(xiàn)象，授粉后18d IP組胚珠發(fā)育出胚乳和胚，而SP組胚珠未發(fā)育出胚乳和胚，并出現(xiàn)內(nèi)外珠被分離。晚期自交不親和的植物花粉能正常萌發(fā)，甚至到達胚珠、穿透胚珠或形成合子（Gibbs，2014）?；谏鲜霈F(xiàn)象，推測BMTGKAN1屬于自交不親和型。在前一周內(nèi)SP組幼果的表現(xiàn)與IP組一致，其原因可能是受到花粉刺激使幼果開始發(fā)育，但由于胚珠不能受精，進而不能正常發(fā)育，最終導致幼果脫落。

3.2植物激素信號轉(zhuǎn)導參與調(diào)控白木通自花授粉后幼果脫落

乙烯信號轉(zhuǎn)導在植物器官脫落中具有不可替代的作用（Shi et al.，2023）。ERF1B是乙烯響應因子（ERFs）家族成員。本研究檢測到ERFIB基因在所有對比組中表達上調(diào)，SP組的ERFIB基因表達較IP組顯著上調(diào)。前人研究也發(fā)現(xiàn)，經(jīng)芐基氨基嘌呤（BA）處理的蘋果脫落的幼果中，乙烯受體ETR1和ETR2表達上調(diào)（Eccher et al.，2015）;經(jīng)乙烯誘導后芒果MiERF3和MiERF113基因在幼果中顯著上調(diào)，引起芒果脫落，且乙烯下游信號調(diào)控基因CTR1和EIN3在脫落區(qū)顯著上調(diào)表達（Rai et al.，2021）。可見，乙烯信號轉(zhuǎn)導在蘋果和芒果幼果脫落中具有重要作用。本研究中ETR1、ETR2、CTRI、EIN2和EIN3基因在SP1vs SP2對比組中均上調(diào)表達，在IP1vs IP2和IP1vs SP1對比組中均不是DEGs，但表現(xiàn)出上調(diào)趨勢；在IP2vs SP2對比組中除ETR2基因不是DEG外，其余均為DEGs，表現(xiàn)為上調(diào)表達，與上述前人研究結(jié)果基本一致，推測該轉(zhuǎn)導途徑參與調(diào)控BMTGKAN1種間授粉后幼果脫落。生長素參與植物的生長發(fā)育，包括植物器官的衰老和脫落（Pathar-kar and Walker，2017）。本研究轉(zhuǎn)錄組數(shù)據(jù)中大量有關(guān)生長素信號轉(zhuǎn)導的轉(zhuǎn)錄本被檢測到，其中，LAX5與LAX2基因表達模式剛好相反，再次驗證了生長素內(nèi)流載體家族成員在果實脫落中發(fā)揮不同作用的結(jié)論（Denisov et al.，2017）。本研究共檢測GH3家族的DEGs共有5個，該家族成員很可能介導花和果的脫落（Cheng et al.，2015）。生長素信號通路負抑制因子（AUX/IAAs）、生長素響應因子（ARFs和SAURs）在生長素信號轉(zhuǎn)導中至關(guān)重要（Shi et al.，2023）。本研究發(fā)現(xiàn)，富集到生長素轉(zhuǎn)導相關(guān)途徑的DEGs共37個，包含ARF、GH3、AUX/IAA、AUX/LAX、SAUR、TIR1等家族基因。其中AUX/IAA、GH3、SAUR、ARF等家族基因與荔枝（Litchi chinen-sis）幼果脫落相關(guān)（Kuang et al.，2012）。本研究從BMTGKAN1SP和IP組幼果中檢測到許多與生長素信號轉(zhuǎn)導途徑相關(guān)的DEGs，推測該途徑參與調(diào)控白木通幼果脫落。綜上所述，植物激素信號轉(zhuǎn)導途徑極有可能參與白木通自花授粉后幼果脫落，尤其是乙烯信號轉(zhuǎn)導途徑和生長素信號轉(zhuǎn)導途徑。白木通自花授粉后幼果脫落的調(diào)控是復雜的，具體調(diào)控網(wǎng)絡還有待深入研究。

3.3白木通自花授粉后幼果脫落的相關(guān)代謝途徑分析

本研究中除植物激素信號轉(zhuǎn)導通路外，亞油酸代謝通路、ABC轉(zhuǎn)運蛋白通路和苯丙烷生物合成通路在IP1vs SP1、IP2vs SP2、SPI vs SP2對比組中存在明顯差異，推測這3個通路與BMTGKAN1幼果的脫落緊密相關(guān)。前人研究發(fā)現(xiàn)，亞油酸代謝途徑與小麥雄性育性、玉米灌漿有關(guān)，而胚敗育往往會導致果實脫落（Malik et al.，2003;Yu et al.，2023;Zhang et al.，2023）。本研究LOX5、LOX?.5和SDP1基因在IP2vs SP2和SP1vs SP2對比組中顯著上調(diào)表達，其他對比組中差異不顯著;代謝物9-OxoODE、9，12，13-TriHOME和13-OxoODE在IP1vs IP2和SP1vs SP2對比組表達上調(diào)，而13S-hydroxyoctadecadienoic acid則相反，且BMTGKAN1SP組胚珠同樣表現(xiàn)出敗育特征，即在自花授粉后18d出現(xiàn)內(nèi)外珠被分離的現(xiàn)象。本研究中檢測到較多ABC轉(zhuǎn)運蛋白通路的基因和代謝物，如ABCC10、ABCG5、L-Lysine和Nopa-line等。ABC轉(zhuǎn)運蛋白通路主要負責利用ATP水解產(chǎn)生能量將底物逆濃度梯度跨膜運輸，該通路與油茶晚期自交不親和相關(guān)（Rees et al.，2009;Zhou et al.，2020）。還有研究發(fā)現(xiàn)，糖、多胺和多肽與植物落花落果相關(guān)，而已鑒定的ABCB亞族成員中大多與激素（生長素）調(diào)控和運輸密切相關(guān)，ABCC亞族成員具有轉(zhuǎn)運代謝物功能（Francisco et al.，2013;Oforiet al.，2018;Shi et al.，2023）。據(jù)此推測，ABC轉(zhuǎn)運蛋白通路與植物激素信號轉(zhuǎn)導通路共同調(diào)控BMTGKAN1自花授粉后幼果脫落。研究表明，苯丙烷生物合成代謝通路與細胞木質(zhì)化有關(guān)，可促進細胞分離過程中的機械性細胞壁斷裂（Liljegren et al.，2000）。此外，Qiu等（2021）在甜櫻桃（Prunus avium）脫落的幼果中也富集到了該通路。本研究在該途徑中篩選到最多的基因是編碼過氧化物酶的基因，共19個，PER家族基因有15個。綜上所述，亞油酸代謝、ABC轉(zhuǎn)運蛋白、苯丙烷生物合成同樣極有可能參與白木通自花授粉后幼果脫落，其具體調(diào)控機制還有待進一步研究驗證。此外，BMTGKAN1由自花授粉引起的幼果脫落是一個高度程序化且較為復雜的生理過程，自花授粉后受精障礙有待進一步研究。植物激素信號轉(zhuǎn)導、ABC轉(zhuǎn)運蛋白、亞油酸代謝、苯丙烷的生物合成及其他代謝通路協(xié)同調(diào)控白木通幼果脫落的調(diào)控網(wǎng)絡還需進一步探究。

4結(jié)論

BMTGKAN1屬于晚期自交不親和植物，其自花授粉后13～18d大量脫落幼果，受其體內(nèi)植物激素信號轉(zhuǎn)導通路、亞油酸代謝通路、ABC轉(zhuǎn)運蛋白通路和苯丙烷生物合成通路調(diào)控，在種植實踐中選擇2個及以上花期相近的白木通栽培品種隔行交叉種植以提高坐果率。

參考文獻：

國家藥典委員會.2020.中華人民共和國藥典[M].北京：中國醫(yī)藥科技出版社.[Pharmacopoeia Commission of the People's Republic of China.2020.Pharmacopoeia of the People's Republic ofChina[M].Beijing：China Medical Science Press.]

熊大勝，曹庸，朱金桃.1996.三葉木通生物學特性研究[J].西南農(nóng)業(yè)大學學報，18（1）：85-90.[Xiong DS，Cao Y，Zhu JT.1996.Studieson biological characteristics ofAke-bia trifoliata[J].Journal of Southwest Agricultural Univer-sity，18（1）：85-90.]

中國科學院中國植物志編輯委員會.2001.中國植物志[M].北京：科學出版社.[Editorial Committee of Flora of China，ChineseAcademy of Sciences.2001.Flora of China[M]Beijing：Science Press.]

Cheng CZ，Zhang LY，Yang XL，ZhongG Y.2015.Profiling gene expression in citrus fruit calyx abscission zone（AZ-C）treated with ethylene[J].Molecular Genetics and Geno-mics，290：1991-2006.doi：10.1007/s00438-015-1054-2.

Denisov Y，Glick S，Zviran T，Ish-ShalomM，LevinA，F(xiàn)aigen-boim A，Cohen Y，Irihimovitch V.2017.Distinct organ-specific and temporal expression profiles of auxin-related genes during mango fruitlet drop[J].Plant Physiology and Biochemistry，115：439-448.doi：10.1016j.plaphy.2017.04.021.

Eccher G，Begheldo M，Boschetti A，Ruperti B，Botton A.2015.Roles of ethylene production and ethylene receptor expression in regulating apple fruitlet abscission[J].Plant Physiology，169（1）：125-137.doi：10.1104/pp.15.00358.

Francisco RM，Regalado A，Ageorges A，BurlaBJ，Bassin B，Eisenach C，Zarrouk O，Vialet S，Marlin T，Chaves MM，Martinoia E，Nagy R.2013.ABCC1，an ATP binding cas-sette protein from grape berry，transports anthocyanidin3-O-glucosides[J].The Plant Cell，25（5）：1840-1854.doi：10.1105/tpc.112.102152.

Franklin-Tong NVE，F(xiàn)ranklin FCH.2003.Gametophyticself-incompatibility inhibits pollen tube growth using different mechanisms[J].Trends in Plant Science，8（12）：598-605.doi：10.1016/j.tplants.2003.10.008.

Garmer LC，Lovatt CJ.2016.Physiological factors affecting flower and fruit abscission of‘Hass'avocado[J].Scientia Horticulturae，199：32-40.doi：10.1016/j.scienta.2015.12.009.

Gibbs PE.2014.Late-acting self-incompatibility-the pariah breeding system in flowering plants[J].New Phytologist，203：717-734.doi：10.1111/nph.12874.

Goldberg EE，Kohn JR，Lande R，Robertson KA，Smith SA，Igié B.2010.Species selection maintains self-incompa-tibility[J].Science，330：493-495.doi：10.1126/science.1194513.Ito Y，Nakano T.2015.Development and regulation of pedicel abscission in tomato[J].Frontiers in Plant Science，6：422.doi：10.3389/fpls.2015.00442.

Kuang JF，WuJY，Zhong HY，LiCQ，ChenJY，Lu WJ，LiJG.2012.Carbohydrate stress affecting fruitlet abscission and expression of genes related to auxin signal transduc-tion pathway in litchi[J].International Journal of Molecu-larSciences，13（12）：16084-16103.doi：10.3390/ijms131216084.

LiL，YaoXH，Zhong CH，Chen XZ，Huang HW.2010.Ake-bia：A potential new fruitcrop in China[J].HortScience，45（1）：4-10.doi：10.21273/hortsci.45.1.4

Liljegren SJ，Ditta GS，Eshed Y，Savidge B，Bowman JL，Yanofsky MF.2000.SHATTERPROOFMADS-box genes control seed dispersal in Arabidopsis[J].Nature，404（6779）：766-770.doi：10.1038/35008089.

Maciag D，Dobrowolska E，Sharafan M，Ekiert H，Tomczyk M，"Szopa A.2021.Akebia quinata and Akebia trifoliata—a"review of phytochemical composition，ethnopharmacologi-cal approaches and biological studies[J].Journal ofEthno-pharmacology，280：114486.doi：10.1016/jjep.2021.114486.

Mahouachi J，Iglesias DJ，Agustí M，Talon M.2009.Delay of early fruitlet abscission by branch girdling in citrus coin-cides with previous increases in carbohydrate and gibberel-lin concentrations[J].Plant Growth Regulation，58：15-23.doi：10.1007/s10725-008-9348-6.

Malik AU，Agrez V，Singh Z2003.Fruitlet abscission of mango in relation to ethylene[J].The Journal of Horticul-tural Science and Biotechnology，78：458-462.doi：10.1080/14620316.2003.11511650.

Mesejo C，Munoz-Fambuena N，Reig C，Martínez-Fuentes A，Agustí M.2014.Cell division interference in newly ferti-lized ovules inducestenospermocarpy in cross-pollinated citrus fruit[J].Plant Science：An International Journal of Experimental Plant Biology，225：86-94.doi：10.1016/j.plantsci.2014.05.019

Mornya PMP，Cheng F.2018.Effect of combined chilling and GA?treatment on bud abortion in forced*Luoyanghong’treepeony（Paeonia sufizuticosa Andr.）[J].Horticultura"Plant Journal，4（6）：250-256.doi：10.1016j.hpj.2018.09.003.

Nakano T，F(xiàn)ujisawa M，Shima Y，Ito Y.2013.Expression pro-filing of tomato pre-abscission pedicels provides insights into abscission zone properties including competence to respond to abscission signals[J].BMC Plant Biology，13：40.doi：10.1186/1471-2229-13-40.

Nakano T，F(xiàn)ujisawa M，Shima Y，Ito Y.2014.The AP2/ERF transcription factor SIERF52functions in flower pedicel abscission in tomato[J].Journalof Experimental Botany，65（12）：3111-3119.doi：10.1093/jxb/eru154.

Ofori PA，Mizuno A，Suzuki M，Martinoia E，Reuscher S，Aoki K，Shibata D，OtagakiS，Matsumoto S，Shiratake K.2018.Genome-wide analysis of ATP binding cassette（ABC）transporters intomato[J].PLoS One，13（7）：e0200854.doi：10.1371/journal.pone.0200854.

Patharkar OR，Walker JC.2017.Advances in abscission sig-naling[J].Journal of Experimental Botany，69（4）：733-740.doi：10.1093/jxb/erx256.

Qiu ZL，Wen Z，Hou QD，QiaoG，YangK，Hong Y，Wen XP.2021.Cross-talk between transcriptome，phytohormone and HD-ZIP gene family analysis illuminates the molecu-lar mechanism underlying fruitlet abscission in sweet cherry（Prumus avium L）[J].BMC Plant Biology，21（1）：173.doi：10.1186/s12870-021-02940-8.

Rai AC，Halon E，Zemach H，Zviran T，Sisai I，Philosoph-Hadas S，Meir S，Cohen Y，Irihimovitch V.2021.Charac-terization of two ethephon-induced IDA-like genes from mango，and elucidation of their involvement in regulating organ abscission[J].Genes（Basel），12（3）：439.doi：10.3390/genes12030439

Rees DC，JohnsonE，Lewinson0.2009.ABC transporters：the power to change[J].Nature Reviews.Molecular Cell Bio-logy，10（3）：218-227.doi：10.1038/nrm2646.

Shi Y，SongB Q，Liang Q，Su DD，LuW，Liu YD，Li ZG.2023.Molecular regulatory events of flower and fruit abscission in horticultural plants[J].Horticultural Plant Journal，9（5）：867-883.doi：10.1016/j.hpj.2023.03.008.

Tatari M，Abdollahi H，Mousavi A.2018.Effect ofpollination on dropping offlowers and fruits in new quince（Cydomia oblonga Mill.）cultivar and promising genotypes[J].Scien-tia Horticulturae，231：126-132.doi：10.1016/j.scienta.2017.10.045.

Vashisth T，Johnson LK，Malladi A.2011.An efficient RNA isolation procedure and identification of reference genes for normalization of gene expression in blueberry[J].Plant Cell Reports，30：2167-2176.doi：10.1007/s00299-011-1121-z.

Waligórski P，Szaleniec M.2010.Prediction of white cabbage"（Brassica oleraceavar.capitata）self-incompatibility based"on neural network and discriminant analysis of complex electrophoretic patterns[J].Computational Biology Che-mistry，34（2）：115-121.doi：10.1016/j.compbiolchem.2010.03.002.

Want EJ，Masson P，Michopoulos F，Wilson ID，Theodoridis G，Plumb RS，Shockcor J，Loftus N，Holmes E，Nicholson JK.2013.Gobal metabolic profiling of animal and human tissues via UPLC-MS[J].Nature Protocols，8（1）：17-32.doi：10.1038/nprot.2012.135.

Yu HL，Wu YQ，Zhang Y，Wang YB.2023.Analysis of key biological pathways for sugaraccumulation during kernel development in different genotypes of sweetcorn[J].Paki-stan Journal of Botany，55：1001-1012.doi：10.30848/PJB2023-3（17）.

Zelena E，Dunn WB，Broadhurst D，F(xiàn)rancis-McIntyre S，Car-roll KM，Begley P，O'Hagan S，Knowles JD，Halsall A，Consortium H，Wilson ID，Kell DB.2009.Development of arobust and repeatable UPLC-MS method for the long term metabolomic study of human serum[J].Analytical"Chemistry，81（4）：1357-1364.doi：10.1021/ac8019366.

Zhang LL，Hou MG，Zhang XG，Cao YY，Sun S，Zhu ZD，Han SB，Chen YH，Ku LX，Duan CX.2023.Integrative transcriptome and proteome analysis reveals maize res-ponses to Fusarium verticillioides infection inside the stalks[J].Molecular Plant Pathology，24（7）：693-710.doi：10.1111/mpp.13317.

Zhao ML，LiCQ，Ma XS，Xia R，Chen JY，Liu XC，Ying PY，Peng MJ，Wang J，ShiC L，LiJ G.2020.KNOX pro-tein KNAT1regulates fruitlet abscission in litchi by rep-ressing ethylene biosyntheticgenes[J].Journal of Experi-mental Botany，71（14）：4069-4082.doi：10.1093/jxb/eraa162.

Zhou JQ，Lu MQ，Yu SS，Liu YY，Yang J，Tan XF.2020.In-depth understanding of Camellia oleifera self-incompati-bility by comparative transcriptome，proteome and metabo-lome[J].International Journal of Molecular Sciences，21（5）：1600.doi：10.3390/ijms21051600.

Zou SY，F(xiàn)eng C，Gao PX，LiTJ，Jia TJ，Huang HW.2022.Germplasm resources and genetic improvement of Akebia：A new fruitcrop in China[J].Plant Diversity，45（6）：712-721.doi：10.1016j.pld.2022.12.001.

（責任編輯 陳燕）