UV-C處理穿心蓮轉錄組分析及穿心蓮內酯合成相關基因挖掘

孫銘陽　徐世強　李靜宇　顧艷　梅瑜　張聞婷　周芳　王繼華

摘要：【目的】對紫外線C（UV-C）處理的穿心蓮葉片進行轉錄組測序分析，挖掘穿心蓮內酯合成相關基因，為深入探究穿心蓮內酯合成相關基因的調控機制和生物學功能提供理論參考。【方法】選用生長60 d的穿心蓮幼苗為材料，經UV-C照射2 h后，利用Illumina HiSeq二代測序平臺進行轉錄組測序，并結合生物信息學軟件和實時熒光定量PCR（qRT-PCR）對所得序列進行功能注釋和相對表達水平驗證，挖掘響應UV-C的穿心蓮內酯合成相關基因。【結果】共注釋到21545個穿心蓮基因，有2137個基因呈顯著差異表達模式，其中，1147個基因顯著上調表達，990個基因顯著下調表達。GO功能注釋結果顯示，UV-C處理造成了氧化脅迫，線粒體氧化磷酸化電子傳遞鏈的NADH脫氫酶基因上調表達。KEGG代謝通路富集結果顯示，二萜類物質穿心蓮內酯的前體（E，E，E）-香葉基香葉基二磷酸（GGPP）合成途徑[甲羥戊酸（MVA）途徑和2-C-甲基-D-赤蘚糖醇-4-磷酸（MEP）途徑]差異表達基因顯著富集。轉錄組測序結果和實時熒光定量PCR（qRT-PCR）檢測結果均顯示，MVA途徑的4種合成酶基因HMGS（CXN00016849）、HMGR（CXN00000058）、MVK（CXN00002900）和MVD（CXN00004398）及GGPP合成關鍵酶基因FPS的表達量顯著上調;MEP途徑中，僅DXS（CXN00013302）的表達量顯著下調。蛋白互作預測結果顯示，2個途徑中的差異表達基因編碼蛋白之間存在互作，HMGS與CYP71家族成員之間，以及FPS1與UGT71、UGT73和UGT76之間的互作網絡豐富。【結論】UV-C照射調控穿心蓮重要生命進程，其中，次生代謝、氧化磷酸化、光合作用等途徑基因轉錄本的相對表達量呈UV-C特異性響應模式。合成穿心蓮內酯前體GGPP的細胞質MVA途徑基因顯著上調表達，推測UV-C誘導穿心蓮內酯合成的主要場所為細胞質。穿心蓮MVA途徑基因及CYP71、UGT71、UGT73和UGT76可作為穿心蓮內酯合成途徑研究的候選分子靶標。

關鍵詞： 穿心蓮;紫外線C（UV-C）;轉錄組分析;穿心蓮內酯;基因挖掘;可變剪切;蛋白互作

中圖分類號： S567.239? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼： A 文章編號：2095-1191（2022）03-0618-10

Transcriptome analysis and discovery of genes related to

andrographolide synthesis of Andrographis paniculata

under UV-C treatment

SUN Ming-yang XU Shi-qiang LI Jing-yu GU Yan MEI Yu ZHANG Wen-ting ZHOU Fang WANG Ji-hua

（1Crops Research Institute，Guangdong Academy of Agricultural Science/Guangdong Provincial Key Laboratory of Crops Genetics Improvement， Guangzhou， Guangdong? 510640， China; 2Guangdong Provincial Engineering and Technology Research Center for Conservation and Utilization of the Genuine Southern Medicinal Resources， Guangzhou，

Guangdong? 510640， China）

Abstract：【Objective】Transcriptome sequencing was performed on the leaves of Andrographis paniculata treated with ultraviolet C （UV-C）， and the genes related to andrographolide synthesis were excavated to provide reference for further exploring the regulatory mechanism and biological functions of andrographolide synthesis-related genes. 【Method】A.paniculata seedlings grown for 60 days were selected as experimental materials. After 2 h of UV-C irradiation， their transcriptome was sequenced using the 2nd-generation sequencing platform， andwith bioinformatics software and qRT-PCR，the sequence were evaluated by gene expression analysis and functional annotation to identify genes related to andrographolide synthesis in response to UV-C. 【Result】A total of 21545 A.paniculata genes were annotated， and 2137 genes showed significant differential expression （SDE） patterns. Among them， 1147 genes were significantly up-regulated and 990 genes were significantly down-regulated. GO annotation results showed that UV-C treatment caused oxidative stress， and the expression of NADH dehydrogenase genes in mitochondrial oxidative phosphorylation electron transport chain were significantly up-regulated. KEGG metabolic pathway enrichment analysis showed that the precursor of andrographolide（E，E，E）-geranylgeranyl diphosphate （GGPP） synthesis pathway （Mevalonate， MVA pathway and 2-C-methyl-D-erythritol-4-phosphate，MEP pathway） SDE genes were significantly enriched. According to the results of RNA-seq and real-time quantitative PCR（qRT-PCR）， the expression of four synthase genes in the MVA pathway， HMGS （CXN00016849）， HMGR （CXN00000058）， MVK （CXN00002900） and MVD （CXN00004398）， as well as GGPP synthesis key enzyme FPS were significantly up-regulated. In the MEP pathway， only the expression of DXS （CXN00013302） was significantly down-regulated. The type of alternative splicing （AS） that commonly occurs in the genes of the two pathways was exon skipping （SE）. The protein interaction prediction showed that there were interactions between the differentially expressed genes （DEGs） of the two pathways， and the interaction networks between HMGS and CYP71 family members， and the ones among FPS1， UGT71， UGT73 and UGT76 were enriched. 【Conclusion】 UV-C irradiation regulates important life processes of A. paniculata. The relative expression of gene transcripts in secondary metabolism， oxidative phosphorylation， photosynthesis and other pathways show UV-C-specific response pattern. The expressions of cytoplasmic MVA pathway genes are significantly up-regulated. It is speculated that the main site of UV-C-induced andrographolide synthesis is the cytoplasm. A. paniculata MVA pathway genes，CYP71， UGT71， UGT73 and UGT76 can be used as candidate molecular targets for andrographolide synthesis pathway research.938D37B0-0BF6-4583-A58F-FB85E89CBAAF

Key words： Andrographis paniculata; UV-C; transcriptome analysis; andrographolide; gene mining; alternative splicing; protein interaction

Foundation items：Guangdong Special Funds for Science and Technology Innovation Strategy（R2019PY-JX003）;Guangdong Province Key Field Research and Development Project（2020B020221001）; Crop Research Institute， Guangdong Academy of Agricultural Sciences/Guangdong Provincial Key Laboratory of Crop Genetic Improvement Open Research Fund （202101）

0 引言

【研究意義】穿心蓮（Andrographis paniculata）為唇形目爵床科穿心蓮屬植物，是傳統中藥材，廣泛分布于東南亞地區（Okhuarobo et al.，2014），自20世紀50年代引入我國后，已被多種復方所收錄（張曉等，2018）。穿心蓮的藥效強勁，有效地抑制了印度1919年爆發的大流感（Sun et al.，2018）。最新研究表明，穿心蓮可緩解新冠肺炎患者的發熱和咳嗽（Jalali et al.，2021）。穿心蓮內酯（C₂₀H₃₀O₅）是穿心蓮發揮藥效的主要次生代謝活性成分，屬于二萜類內酯，食味發苦（Abu-Ghefreh et al.，2009）。然而，穿心蓮內酯的含量受不同品種、產地、栽培技術、收獲時間和加工技術等因素影響較大，很難滿足穿心蓮產品的市場化需求。因此，開展穿心蓮內酯合成途徑相關基因的分子機制研究，對高效挖掘調控穿心蓮次生代謝的主效基因具有重要意義。【前人研究進展】2019年之前，由于穿心蓮缺少基因組信息，相關研究重點主要為臨床藥理活性和代謝產物分離等方面。2019年起，已陸續公開2個穿心蓮基因組數據庫，但仍存在基因功能注釋不全面和部分基因組位置組裝質量不高等問題，導致穿心蓮二萜類物質合成相關功能基因研究仍處于較淺層面。Wang等（2019）利用克隆技術獲得4個假定的穿心蓮二萜合成關鍵酶基因（ApGGPPS）的cDNA序列，并通過原核表達驗證了其中3個基因能夠表達目的蛋白。Srinath等（2020）對1個假定的穿心蓮MVA途徑限速酶基因（ApHMGR）的cDNA進行分離，結合qRT-PCR和高效液相色譜技術證實了該基因的表達量與穿心蓮內酯的含量呈正相關。目前，大量研究證實，紫外線C（UV-C，波長介于100～280 nm）照射可增強植物次生代謝產物的積累。該手段簡便且見效較快，在草藥、蔬菜和水果產業中得到廣泛應用（Ma et al.，2019）。Xu等（2018）研究發現，草莓（Fragaria ana-nassa Duchesne）經UV-C照射后其萜類物質含量上升26%。Alsoufi等（2019）研究發現，UV-C處理的金盞花（Calendula officinalis）毛狀根中萜類物質含量是未處理的8.5倍。然而，植物細胞過度暴露于UV-C下會誘導其產生活性氧（Reactive oxygen species，ROS），從而促進光氧化，損害光系統（Photosystem II，PSII），導致光合活性降低（Darras et al.，2020）。轉錄組測序技術能在RNA水平反映基因表達的時空特異性和脅迫應激特異性，現已成為探索藥用植物次生代謝分子機制的重要手段（李鉑等，2020）。目前已超過50種傳統藥用植物實現轉錄組測序數據解析，包括穿心蓮、人參（Panax ginseng）和冬蟲夏草（Cordyceps）等（李慧等，2018）。Garg等（2015）利用轉錄組測序技術對穿心蓮不同組織的二萜類物質合成相關基因表達水平進行解析。Gao等（2019）利用茉莉酸甲酯處理的穿心蓮轉錄組測序數據評估了穿心蓮內酯合成關鍵基因轉錄本的可變剪切（Alternative splicing，AS）情況。【本研究切入點】目前，鮮見有關UV-C處理穿心蓮的轉錄組測序的研究報道。【擬解決的關鍵問題】對UV-C處理穿心蓮葉片進行轉錄組測序，并結合生物信息學軟件和實時熒光定量PCR（qRT-PCR）對所得序列進行功能注釋和相對表達水平驗證，重點挖掘響應UV-C的穿心蓮內酯合成相關基因，以期解析UV-C對穿心蓮內酯合成的影響，為探究相關基因的調控機制提供參考。

1 材料與方法

1. 1 試驗材料

選用福建漳州生產用穿心蓮種，將其種植并培養于廣東省農業科學院作物研究所農作物遺傳改良重點實驗室的光照培養箱（光照/黑暗：16 h/8 h;溫度：28 ℃）中，土壤為草炭土（HAWITA），培養期間保持土壤濕潤。幼苗培養60 d，生長時間為2020年4月25日—6月24日。主要試劑：RNA提取試劑盒購自寶日醫生物技術（北京）有限公司;RNA反轉錄試劑盒購自諾唯贊生物科技股份有限公司;實時熒光定量PCR（qRT-PCR）擴增試劑盒購自伯樂生命醫學產品有限公司;引物由湖北擎科生物科技有限公司合成。主要設備儀器：GXZ智能型人工氣候箱（寧波江南儀器廠）、Illumina HiSeq二代轉錄組測序儀（上海仁科生物科技有限公司）、實時熒光定量PCR儀CFX96TM System（伯樂生命醫學產品有限公司）。

1. 2 樣品處理

將幼苗從培養箱中移出，分為2組，分別用白光（對照組，CK）和UV-C（處理組）照射2 h。利用UV-C照度計測得處理組幼苗接收的UV-C照度為3.80 W/m2。輻照2 h后，UV-C的總照射劑量為27.36 kJ/m2。隨后剪取處理組和對照組幼苗的倒三葉，迅速置于液氮速凍，用于RNA提取。對照組和處理組設3次獨立的生物學重復，共計6份樣品，分別表示為CK1、CK2、CK3和UV1、UV2、UV3。938D37B0-0BF6-4583-A58F-FB85E89CBAAF

1. 3 轉錄組測序數據分析

利用RNAiso Plus試劑提取樣品的總RNA，利用Nano Photometer Spectrophotometer、Qubit 2.0 Fluorometer和Agilent 2100 Bioanalyzer精細檢測RNA的純度和完整性。高質量RNA經過去除核糖體RNA（rRNA）、反轉錄合成雙鏈cDNA、cDNA末端修復及cDNA加接頭處理，所得產物于Illumina HiSeq二代測序平臺進行文庫構建。利用fastp對序列進行質控，得到Clean reads，并計算Q30值和GC含量（%）。利用bowtie2將Clean reads比對到穿心蓮rRNA數據庫，去除比對上的序列，剩余序列用于后續分析。利用HISAT2將序列進行穿心蓮基因組比對，并計算基因比對區域。從NCBI數據庫（https：//www.ncbi.nlm.nih.gov/genome/？term=Andrographis+paniculata）下載穿心蓮基因組數據。利用Stringtie重構轉錄本后，篩選出在本次測序結果中發現但未被參考基因組收錄的基因，定義為新基因。

利用RSEM計算每個樣本中所有基因的表達量。利用DESeq分析組間差異表達基因（Differentially expressed genes，DEGs）。基于差異分析結果，篩選FDR<0.05，|log₂FC|>1的基因為顯著差異表達基因。利用NCBI中的BLASTx將轉錄本序列對應的多肽序列比對到蛋白數據庫GO（Gene ontology）和KEGG（Kyoto encyclopedia of genes and genomes）中。利用rMATS檢測外顯子跳躍型（Skipping exon，SE）、外顯子互斥型（Mutually exons exclusively，MXE）、5′端或3′端可變型（Alternative 5′donor or 3′ acceptor sites，A5/A3SS）及內含子保留型（Intron retention，IR）AS轉錄本。利用STRING蛋白質互作數據庫中的互作關系選項進行差異表達基因編碼蛋白的互作網絡分析，互作指數（0～1000）越高，蛋白間存在互作的可能性越大。本研究挑選互作指數超過500的差異表達基因，利用Cytoscape 3制作蛋白互作圖。利用TBtools將基因的FPKM值進行log2后繪制表達水平熱圖。

1. 4 qRT-PCR檢測

提取處理組和對照組的總RNA，反轉錄合成cDNA，以cDNA為模板對穿心蓮內酯合成相關基因進行qRT-PCR檢測。 反應體系：TB Green Fast qPCR Mix（2×）10.0 μL，10 μmol/L上、下游引物0.8 μL，150 ng/μL cDNA模板 2 μL ，ddH₂O補充至20.0 μL。擴增程序：95 ℃? 30 s;95 ℃? 5 s，60 ℃ 10 s，進行40個循環。以Actin（GenBank登錄號JX444056）基因作為內參（Wang et al.，2019）。試驗所用引物如表1所示。

1. 5 統計分析

處理組和對照組均設3次獨立的生物學重復，每個生物學重復設3次技術重復，根據2^?ΔΔCt算法確定基因的相對表達量，數據用獨立生物學實驗的平均值±樣本標準誤差（SE）表示。利用Excel 2019進行計算及數據整理。

2 結果與分析

2. 1 UV-C處理下穿心蓮轉錄組的測序質量

原始數據經過濾后得到Clean reads，其中，從CK1、CK2和CK3分別得到36156632、47112300和37869802條Clean reads，從UV1、UV2和UV3分別得到39754768、53187418和46009804條Clean reads。堿基質量分析結果（表2）顯示，6份樣品的Q30值均超過93.00%，質量較高。Clean reads的堿基組成越平衡，代表測序質量越好。6份測序樣品的平均GC含量為51.28%，堿基含量趨于平均，說明測序質量好，可用于后續研究。

雖然總RNA樣品中的rRNA已被去除，但為避免殘留，利用短reads比對工具將Clean reads與穿心蓮rRNA數據庫進行比對。未比對到rRNA數據庫中的Unmapped reads用于后續功能分析。本研究中，UV-C處理組樣品的平均Unmapped reads數量為46317330條，占平均Clean reads數量的96.81%;對照組的平均Unmapped reads數量為40379578，占平均Clean reads數量的99.54%。將去除rRNA的Clean reads與穿心蓮基因組進行比對，結果如表2所示。6份樣品的平均基因組比對率高達95.29%。將比對到基因組上的序列進行基因結構分析，結果發現UV-C處理使穿心蓮RNA來源呈現出由基因區向基因間區發展的趨勢，編碼蛋白的RNA含量降低，非編碼蛋白的RNA含量平均增加22.64%（表3）。

2. 2 差異基因篩選及表達分析結果

穿心蓮基因組共包含25428個基因的注釋信息（Sun et al.，2018）。本研究映射到其中的21545個（84.73%）基因，其原因可能為基因組組裝存在部分偏差、基因組和轉錄組送測樣品的生長時期、組織不同等，使轉錄組測序結果的部分信息無法比對到基因組中。本研究對轉錄組序列進行重構分析，結果發現678個新基因，其中有476個得到GO或KEGG數據庫的功能注釋。

對處理組和對照組數據進行基因表達量分析，結果顯示，共5193個基因上調表達，16989個基因下調表達。經UV-C處理后，將表達量FDR<0.05且|log2FC|>1的基因作為顯著差異表達基因，結果如圖1所示。UV-C照射處理后共產生2137個顯著差異表達基因，顯著上調表達的基因有1147個，顯著下調表達的基因有990個。938D37B0-0BF6-4583-A58F-FB85E89CBAAF

2. 3 差異表達基因的GO功能注釋分析結果

GO是國際標準化的基因功能分類體系。將差異表達基因進行GO功能注釋，結果顯示，分別有1039、1069和1242個差異表達基因歸屬于細胞組成、分子功能和生物進程三大類。由表4可知，UV-C照射穿心蓮后，差異表達基因主要參與的生物進程為響應有機氮化合物、響應內源性刺激及信號的轉導和傳遞;主要分子功能為氧化磷酸化途徑的NADH脫氫酶活性及電子載流活性;功能發揮的主要位置為生物膜。

UV-C照射會刺激植物線粒體和葉綠體釋放ROS，破壞生物膜的功能及線粒體膜上的電子傳遞鏈，導致氧化脅迫（Darras et al.，2020）。GO功能注釋結果表明，UV-C照射處理對穿心蓮造成了氧化脅迫。NADH脫氫酶作為線粒體電子傳遞鏈中的第一個復合體，是線粒體氧化磷酸化的“入口酶”。細胞色素c（Cytochrome c，Cyt c）氧化酶是電子傳遞鏈的末端酶，是將H+由基質抽提到膜間隙的質子泵（Gao et al.，2008）。本研究發現，線粒體氧化磷酸化途徑的19個NADH脫氫酶基因、4個Cyt c氧化酶基因和5個F型ATP酶基因經UV-C處理后顯著上調表達（圖2）。此外，生物膜的主要成分為不飽和脂肪酸，可維持生物膜的流動性，保證細胞正常的生理功能。本研究中有3個不飽和脂肪酸合成關鍵基因的積累量受明顯抑制，分別為CXN00015046（超長鏈3-氧代酰基輔酶A還原酶1，log₂FC=-1.03）、CXN00009378（超長鏈3R-3-羥酰基輔酶A脫水酶PASTICCINO2，log2FC=-1.45）和CXN00013503（超長鏈烯酰輔酶A還原酶，log2FC=-1.69）。

2. 4 差異表達基因的KEGG代謝通路富集分析結果

KEGG是生物代謝通路分類的數據庫。本研究差異表達基因共富集到111個通路中，顯著富集的前20個代謝通路如圖3所示。經UV-C處理后，響應環境變化的基礎生命進程，如絲裂原活化蛋白激酶（Mitogen-activated protein kinase，MAPK）級聯途徑、氧化磷酸化、光合作用和生物膜流動性等途徑的差異基因得到富集。同時，穿心蓮主要藥效次生代謝物為二萜類化合物。由圖3（標紅處）可知，UV-C處理使穿心蓮二萜類物質前體骨架合成相關基因得到顯著富集。

MAPK是一組能被不同的胞外刺激激活的絲氨酸-蘇氨酸蛋白激酶，是信號從細胞表面傳導到細胞核內部的重要傳遞者。本研究發現，MAPK通路誘導過氧化氫、乙烯和植物抗毒素的合成，從而增強病原體抗性、促進氣孔發育及維持活性氧穩態。此外，光合作用產生的有機物是植物營養物質的重要來源。KEGG代謝通路注釋結果顯示，UV-C照射使PS復合物中的4個反應中心基因和5個葉綠素脫輔基基因上調表達;接收光信號的光敏色素B6/f復合物中，2個光敏色素基因被誘導表達，但3個葉綠素合成途徑基因下調表達。

2. 5 穿心蓮二萜前體合成途徑基因的表達模式及AS分析結果

二萜類物質的前體（E，E，E）-香葉基香葉基二磷酸[（E，E，E）-geranyl geranyl diphosphate，GGPP]的合成途徑為位于細胞質中的甲羥戊酸（Mevalonate，MVA）途徑和位于質體中的2-C-甲基-D-赤蘚糖醇-4-磷酸（2-C-methyl-D-erythritol-4-phosphate，MEP）途徑。由圖4可知，經UV-C處理后，MVA途徑的4種合成酶基因HMGS（CXN00016849）、HMGR（CXN00 000058）、MVK（CXN00002900）和MVD（CXN0000 4398）表達呈顯著上調。同時，GGPP合成關鍵酶基因FPS（CXN00005067）表達也顯著上調。而MEP途徑中，僅1個DXS（CXN00013302）基因呈顯著差異表達，且顯著下調。

AS是真核生物體內普遍存在的基因調控形式，主要分為SE、MXE、A5/A3SS和IR。本研究共檢測到8881個基因發生SE，705個基因發生MXE，1429個基因發生A5SS，1829個基因發生A3SS，1532個基因發生IR。由圖4可知，經UV-C處理后，12種萜類骨架合成途徑酶基因發生了AS事件，其中，SE發生的頻率最高（75%），其次為A3SS（58%）。

2. 6 二萜前體合成通路相關基因編碼的蛋白互作預測結果

為進一步明確UV-C照射對穿心蓮內酯合成途徑基因的影響，利用蛋白互作數據庫對差異表達基因進行互作網絡分析（互作指數>500），結果如圖5所示。UV-C處理后，MVA途徑中顯著差異表達基因編碼的蛋白HMGS、HMGR、MVK、MVD、DXS和FPS之間存在互作。此外，二萜類物質的前體GGPP形成后，經一系列特異性催化修飾形成具有獨特藥理活性的終產物。UDP依賴型糖基轉移酶（UDP-dependent glycosyltransferases，UGT）家族的UGT71、UGT73、UGT76、UGT85和UGT88成員參與二萜類物質的修飾（Sun et al.，2018）。本研究UV-C處理后發現，HMGS和MVK均與UGT71D13互作，GGPP合成關鍵酶FPS1與UGT71、UGT73、UGT76成員互作。細胞色素P450（Cytochrome P450，CYP450）基因家族的CYP71和CYP76成員參與唇形目植物二萜類物質的合成（Banerjee and Hamberger，2018）。UV-C處理下，穿心蓮HMGS與7種CYP71家族成員互作，MVK和DXS分別與CYP71AU50和CYP710A11互作。參與二萜生物合成的轉錄因子（Transcription factor，TF）家族有AP2/ERF、bHLH和WRKY（Sun et al.，2018）。UV-C處理使穿心蓮HMGR與WRKY53互作。938D37B0-0BF6-4583-A58F-FB85E89CBAAF

2. 7 二萜前體合成通路基因的表達模式分析結果

利用qRT-PCR檢測2.5中顯著上調或顯著下調的基因HMGR（CXN00000058）、MVK（CXN0000 2900）、MVD（CXN00004398）、FPS（CXN00005067）、DXS（CXN00013302）和HMGS（CXN00016849）的相對表達量，再與轉錄組測序結果進行比較，結果如圖6所示。6個基因相對表達量的qRT-PCR檢測結果與轉錄組測序結果一致，表明UV-C處理后的轉錄測序數據可信度較高。

3 討論

本研究利用UV-C對穿心蓮幼苗進行照射處理，并進行轉錄組測序分析，測序結果堿基含量高度平衡，數據質量高。進一步分析發現，UV-C處理使穿心蓮轉錄組更龐大，產生更多rRNA，其原因可能是植株響應環境光變化，導致蛋白翻譯加速。轉錄本成分分析結果顯示，非編碼RNA在UV-C處理組中的占比大幅增加（平均含量增加22.64%），是穿心蓮響應UV-C的主要RNA成分。同時，本研究UV-C處理使穿心蓮二萜合成途徑差異表達基因產生豐富的AS亞型，其中SE事件最為普遍，說明其可通過單個基因的AS事件作出高度協調UV-C照射的分子響應。穿心蓮經UV-C處理后，共產生1147個顯著上調基因和990個顯著下調基因，差異表達基因的數量明顯高于番茄（Lycopersicon esculentum）（劉長虹等，2012）、煙草（Nicotiana tabacum）（易小艷，2015）等植物，為穿心蓮次生代謝途徑基因研究打下基礎。

穿心蓮內酯前體GGPP合成途徑分為胞質MVA途徑和質體MEP途徑。轉錄組測序結果和qRT-PCR結果顯示，UV-C處理下，MVA途徑基因HMGS、HMGR、MVK和MVD以及GGPP合成必需基因FPS的轉錄本積累被顯著誘導，但大多數MEP途徑基因未呈差異表達，僅DXS基因呈顯著下調表達模式。表明UV-C處理后，穿心蓮主要依靠MVA途徑而非MEP途徑進行二萜前體的合成，其原因可能為葉綠體是光合作用的主要場所，感應環境光的變化，當高劑量UV-C照射葉片時，以葉綠體為主的質體完整性受到一定程度破壞，MEP途徑合成GGPP的效率受到抑制（Darras et al.，2020）。此外，MVA和MEP途徑的差異表達基因編碼的蛋白之間存在復雜的互作網絡，主要為HMGS、HMGR、MVK、MVD、DXS和FPS之間的互作。Zheng等（2021）研究發現，葡萄VvHMGS、VvHMGR、VvMVK、VvMVD和VvFPS/GGPPS蛋白之間存在互作，與本研究結果相符。同時，CYPs和UGTs家族成員在調節萜類生物合成過程中發揮重要作用，其代表性成員為CYP71、 UGT71和UGT73（Liu et al.，2016;Cheng et al.，2020）。本研究UV-C處理后的穿心蓮CYP71與HMGS之間以及UGT71、UGT73和UGT76與FPS1之間也存在緊密的互作關系。

此外，當植株長期暴露于UV-C下，過量產生的 ROS會使線粒體跨膜電位丟失，導致Cyt c釋放，造成氧化脅迫并觸發細胞程序性凋亡（Programmed cell death，PCD）（Gao et al.，2008）。本研究穿心蓮細胞通過上調Cyt c氧化酶基因和NADH脫氫酶基因表達抑制UV-C引起的Cyt c釋放并促進膜電位重建，減緩PCD;穿心蓮MAPK級聯通路基因經UV-C處理后顯著上調表達，調控細胞生長、分化和應激等多種重要生理過程相關基因的表達水平。此外，UV-C會損害植物細胞PSII的光合效率和CO2 的吸收能力，加速葉綠素a/b的降解（Rusaczonek et al.，2015）。本研究發現，UV-C使穿心蓮葉綠素合成相關基因的表達水平降低后，細胞試圖恢復光合效率而提高了光捕獲PS復合體基因的表達水平。本研究結果為穿心蓮生長發育和次生代謝機制的研究提供新思路，后續將通過高效液相色譜分析、基因克隆、轉基因植株培育等手段深入研究穿心蓮內酯合成關鍵基因的分子功能。

4 結論

UV-C照射調控穿心蓮重要生命進程，其中，次生代謝、氧化磷酸化、光合作用等途徑基因轉錄本的相對表達量呈UV-C特異性響應模式。合成穿心蓮內酯前體GGPP的細胞質MVA途徑基因顯著上調表達，推測UV-C誘導穿心蓮內酯合成的主要場所為細胞質。穿心蓮MVA途徑基因及CYP71、UGT71、UGT73和UGT76可作為穿心蓮內酯合成途徑研究的候選分子靶標。

參考文獻：

李鉑，彭亮，王楠，孫曉春，何懿菡，宋忠興，張永強，唐志書，黃文靜. 2020. 聚乙二醇模擬干旱脅迫下半夏莖段轉錄組分析[J]. 中草藥，51（21）：5579-5589. [Li B，Peng L，Wang N，Sun X C，He Y H，Song Z X，Zhang Y Q，Tang Z S，Huang W J. 2020. Transcriptome analysis of stems of Pinellia ternate under drought stress simulated by PEG[J]. Chinese Traditional and Herbal Drugs，51（21）：5579-5589.] doi：10.7501/j.issn.0253-2670.2020.21.023.

李慧，馬德志，姜明，劉溪，陳曉婷，劉振艷. 2018. 傳統藥用植物轉錄組研究進展[J]. 中醫藥信息，35（6）：114-120. [Li H，Ma D Z，Jiang M，Liu X，Chen X T，Liu Z Y. 2018. Research progress on transcriptome sequencing of traditional medicinal plants[J]. Information on Traditional Chinese Medicine，35（6）：114-120.] doi：10.19656/j.cnki. 1002-2406.180192.938D37B0-0BF6-4583-A58F-FB85E89CBAAF

劉長虹，蔡路昀，韓曉旭，應鐵進. 2012. UV-C照射后番茄差異表達基因的基因芯片檢測[J]. 農業機械學報，43（1）：134-140. [Liu C H，Cai L Y，Han X X，Ying T J. 2012. Microarray analysis of gene expression profile induced by postharvest UV-C irradiation in tomato fruit[J]. Tran-sactions of the Chinese Society for Agricultural Machi-nery，43（1）：134-140.] doi：10.6041/j.issn.1000-1298.2012. 01.024.

易小艷. 2015. UV-B和UV-C對煙苗形態生理及基因表達譜的影響 [D]. 成都：四川農業大學. [Yi X Y. 2015. Effects of UV-B and UV-C on tobacco seedlings[D]. Chengdu：Sichuan Agricultural University.]

張曉，唐力英，吳宏偉，許夢瑩，郭日新，于現闊，魯亞奇，羅寒燕，劉地發，王祝舉，李志勇. 2018. 穿心蓮現代研究進展[J]. 中國實驗方劑學雜志，24（18）：222-234. [Zhang X，Tang L Y，Wu H W，Xu M Y，Guo R X，Yu X K，Lu Y Q，Luo H Y，Liu D F，Wang Z J，Li Z Y. 2018. Development of modern research on Andrographis paniculate[J]. Chinese Journal of Experimental Traditional Medical Formulae，24（18）：222-234.] doi：10.13422/j.cnki.syfjx.2018 1516.

Abu-Ghefreh A A，Canatan H，Ezeamuzie C H. 2009. In vitro and in vivo anti-inflammatory effects of andrographolide[J]. International Immunopharmacology，9（3）：313-318. doi：10.1016/j.intimp.2008.12.002.

Alsoufi A S M，Paczkowski C，Dlugosz M，Szakiel A. 2019. Influence of selected abiotic factors on triterpenoid biosynthesis and saponin secretion in marigold（Calendula officinalis L.） in vitro hairy root cultures[J]. Molecules，24（16）：2907. doi：10.3390/molecules24162907.

Banerjee A，Hamberger B. 2018. P450s controlling metabolic bifurcations in plant terpene specialized metabolism[J]. Phytochemistry Reviews，17（1）：81-111. doi：10.1007/s11101-017-9530-4.

Cheng Y，Liu H B，Tong X J，Liu Z M，Zhang X，Li D L，Jiang X M，Yu X H. 2020. Identification and analysis of CYP450 and UGT supergene family members from the transcriptome of Aralia elata（Miq.） seem reveal candidate genes for triterpenoid saponin biosynthesis[J]. BMC Plant Biology，20（1）：214. doi：10.1186/s12870-020-02411-6.

Darras A，Tsikaloudakis G，Lycoskoufis I，Lycoskoufis I，Dimitrisdis C I，Karamousantas D. 2020. Low doses of UV-C irradiation affects growth，fruit yield and photosynthetic activity of tomato plants[J]. Scientia Horticulturae，267：109357. doi：10.1016/j.scienta.2020.109357.

Gao C J，Xing D，Li L L，Zhang L R. 2008. Implication of reactive oxygen species and mitochondrial dysfunction in the early stages of plant programmed cell death induced by ultraviolet-C overexposure[J]. Planta，227（4）：755-767. doi：10.1007/s00425-007-0654-4.

Gao H，Li F H，Xu Z C，Huang C H，Xiong C，Jiang C H，Xie N，Leng L，Zhang Y，Yousaf? Z，Liu X，Sun W. 2019. Genome-wide analysis of methyl jasmonate-regulated isoform expression in the medicinal plant Andrographis paniculata[J]. Industrial Crops and Products，135（1）：39-48. doi：10.1016/j.indcrop.2019.04.023.938D37B0-0BF6-4583-A58F-FB85E89CBAAF

Garg A，Agrawal L，Misra R C，Sharma S，Ghosh S. 2015. Andrographis paniculata transcriptome provides molecular insights into tissue-specific accumulation of medicinal diterpenes [J]. BMC Genomics，16（1）：659. doi：10.1186/s12864-015-1864-y.

Jalali A，Dabaghian F，Akbrialiabad H，Foroughinia F，Zarshenas M M. 2021. A pharmacology-based comprehensive review on medicinal plants and phytoactive constituents possibly effective in the management of COVID-19[J]. Phytotherapy Research，35（4）：1925-1938. doi：10.1002/ptr.6936.

Liu X L，Cao F L，Cai J F，Wang H L. 2016. The molecular cloning and expression analysis of a CYP71 gene in Ginkgo biloba L[J]. Notulae Botanicae Horti agrobotanici Cluj-Napoca，44（1）：77-84. doi：10.15835/nbha44110214.

Ma M，Zhang H X，Xie Y J，Yang M，Tang J F，Wang P，Yang R Q，Gu Z X. 2019. Response of nutritional and functio-nal composition，anti-nutritional factors and antioxidant activity in germinated soybean under UV-B radiation[J]. LWT-Food Science and Technology，118（4）：108709. doi：10.1016/j.lwt.2019.108709.

Okhuarobo A，Falodun J E，Erharuyi O，Imieje V，Falodun A，Langer P. 2014. Harnessing the medicinal properties of Andrographis paniculata for diseases and beyond：A review of its phytochemistry and pharmacology[J]. Asian Pacific Journal of Tropical Disease，4（3）：213-222. doi：10.1016/S2222-1808（14）60509-0.

Rusaczonek A，Czarnocka W，Kacprzak S，Witoń D，?lesak I，Szechyńska-Hebda M，Gawroński P，Karpinski S. 2015. Role of phytochromes A and B in the regulation of cell death and acclimatory responses to UV stress in Arabidopsis thaliana[J]. Journal of Experimental Botany，66（21）：6679-6695. doi：10.1093/jxb/erv375.

Srinath M，Bindu B B V，Shailaja A，Giri C C. 2020. Isolation，characterization and in silico analysis of 3?hydroxy?3? methylglutaryl? coenzyme A reductase（HMGR） gene from Andrographis paniculata（Burm.f） Nees[J]. Molecular Biology Reports，47（10）：639-654. doi：10.1007/s11033-019- 05172-0.

Sun W，Leng L，Yin Q G，Xu M M，Huang M K，Xu Z C，Zhang Y J，Yao H，Wang C X，Xiong C，Chen S，Jiang C H，Xie N，Zheng X L，Wang Y，Song C，Peters R J，Chen S L. 2018. The medicinal plant Andrographis paniculata genome provides insight into biosynthesis of the bioactive diterpenoid neoandrographolide[J]. The Plant Journal，97（5）：841-857. doi：10.1111/tpj.14162.

Wang J，Lin H X，Su P，Chen T，Guo J，Gao W，Huang L Q. 2019. Molecular cloning and functional characterization of multiple geranylgeranyl pyrophosphate synthases（ApGGPPS） from Andrographis paniculata[J]. Plant Cell Reports，38（1）：117-128. doi：10.1007/s00299-018-2353-y.

Xu Y Q，Charles M T，Luo Z S，Mimee B，Tong Z C，Véronneau P Y，Roussel D，Rolland D. 2018. Ultraviolet-C priming of strawberry leaves against subsequent Mycosphaerella fragariae infection involves the action of reactive oxygen species，plant hormones，and terpenes[J]. Plant Cell and Environment，42（3）：815-831. doi：10.1111/ pce.13491.

Zheng T，Guan L B，Yu K，Salman H M，Nasim M，Liu Z J，Li T，Zhang K K，Jiu S T，Jia H F，Fang J G. 2021. Expressional diversity of grapevine 3-hydroxy-3-methylglutaryl-CoA reductase（VvHMGR） in different grapes genotypes[J]. BMC Plant Biology，21（1）：279. doi：10.1186/s12870-021-03073-8.

（責任編輯 陳 燕）938D37B0-0BF6-4583-A58F-FB85E89CBAAF