贛南早’臍橙在干旱脅迫下的生理及轉(zhuǎn)錄組研究

劉林芝　歐陽歡　李興濤　陳健美

‘

摘 ?要：我國臍橙產(chǎn)區(qū)的季節(jié)性干旱對臍橙產(chǎn)量和品質(zhì)影響較大。‘贛南早臍橙作為一個新品種臍橙，目前在我國臍橙產(chǎn)區(qū)已大面積推廣。為深入了解這個新品種的耐旱性，探究早熟品種‘贛南早臍橙應(yīng)對干旱脅迫的調(diào)控機制，以‘贛南早臍橙與‘紐荷爾臍橙（對照）為材料，測定比較不同干旱脅迫程度下二者光合作用、干旱相關(guān)生理指標等差異，并通過RNA-Seq分析比較轉(zhuǎn)錄水平差異及抗氧化物酶基因表達調(diào)控。結(jié)果表明：干旱脅迫下‘贛南早臍橙凈光合速率（P_n）、蒸騰速率（T_r）、氣孔導(dǎo)度（G_s）均顯著高于‘紐荷爾臍橙;隨著干旱脅迫程度增加，‘贛南早臍橙葉片較‘紐荷爾臍橙更舒展，‘贛南早臍橙相對電導(dǎo)率和丙二醛含量顯著低于‘紐荷爾臍橙，而保護酶超氧化物歧化酶（SOD）和過氧化物歧化酶（POD）活性變化幅度更大;‘贛南早臍橙和‘紐荷爾臍橙葉片可溶性糖含量無顯著性差異，復(fù)水后，‘紐荷爾臍橙葉片中可溶性糖含量極顯著高于‘贛南早臍橙。轉(zhuǎn)錄組測序分析表明，干旱脅迫0、10、20 d時，‘贛南早臍橙和‘紐荷爾臍橙間DEGs數(shù)量分別為1266、683、658個。GO富集分析顯示，在干旱脅迫過程中‘贛南早臍橙差異基因主要集中在細胞蛋白修飾過程、高分子修飾作用、含磷化合物代謝過程、蛋白修飾過程等通路，而‘紐荷爾臍橙未見明顯富集。KEGG富集分析顯示，除了富集于淀粉及蔗糖通路和氨基酸及核苷酸糖代謝途徑，‘贛南早臍橙其他差異基因富集途徑與‘紐荷爾臍橙基本一致。差異基因轉(zhuǎn)錄因子分析顯示二者在ERF家族、MYB家族、NAC家族、MYB_related家族、WRKY家族、bHLH家族、HB-other家族、HSF家族、B3家族和bZIP家族均有分布，此外，‘贛南早臍橙特異分布于GRAS家族。根據(jù)轉(zhuǎn)錄組分析篩選出抗氧化酶相關(guān)基因30個，其中上調(diào)表達48%，下調(diào)表達52%。本研究結(jié)果為‘贛南早響應(yīng)干旱脅迫的生理變化提供理論依據(jù)，并為其抗旱性研究提供分子基礎(chǔ)。

關(guān)鍵詞：臍橙;干旱脅迫;光合作用;酶活性;丙二醛;轉(zhuǎn)錄組分析中圖分類號：S666.4??????文獻標識碼：A

Physiological and Transcriptome Analysis of ‘Gannan Zao Navel Orange under Drought Stress

Abstract： The seasonal drought in the navel orange producing areas in China has a great impact on the yield and quality of navel orange. ‘Gannan zao navel orange， a new variety， is widely popularized in navel orange producing areas in China. In order to deeply understand the drought tolerance of ‘Gannan zao and explore its regulation mechanism on drought stress， the differences of photosynthesis and drought related physiological indexes between ‘Gannan zao and ‘Newhall under different drought stress were measured， and the transcriptional differences and the regulation of antioxidant enzyme gene expression were compared by RNA-seq analysis. The results showed that the net photosynthetic rate （P_n）， stomatal conductance （G_s） and transpiration rate （T_r） of ‘Gannan zao were significantly higher than those of ‘Newhall under drought stress. With the increase of drought stress， the leaves of ‘Gannan zao were more stretched than that of ‘Newhall， and the relative conductivity and malondialdehyde content of ‘Gannan zao were significantly lower than that of ‘Newhall， while the activities of protective enzymes superoxide dismutase and peroxide dismutase changed more. There was no significant difference in soluble sugar content between ‘Gannan zao and ‘Newhall. After rehydration， the soluble sugar content of ‘Newhall was significantly higher than that of ‘Gannan zao . Transcriptome sequencing analysis showed that the number of DEGs between ‘Gannan Zao and ‘Newhall was 1266， 683 and 658 respectively at 0， 10 and 20 days of drought stress. Go enrichment analysis showed that the differential genes of ‘Gannan zao were mainly concentrated in the process of cell protein modification， protein modification， polymer modification and phosphorus compound metabolism， while ‘Newhall was not significantly enriched. KEGG enrichment analysis showed that the differential gene concentration pathways of ‘Gannan Zao and ‘Newhall were basically the same， while ‘Gannan zao was also enriched in starch and sucrose pathway and amino acid and nucleotide sugar metabolism pathway. ‘Gannan zao was also enriched in starch and sucrose pathway and amino acid and nucleotide sugar metabolism pathway. Both transcription factors were in ERF family， MYB family， NAC family， MYB_Related family， WRKY family， bHLH family， HB-other family， HSF family， B3 family and bZIP family were distributed. In addition， ‘Gannan zao was specifically distributed in GRAS family. According to the transcriptome molecular results， 30 antioxidant enzyme related genes were screened， of which 48% were up-regulated and 52% were down regulated. The results of this study would provide a theoretical basis for the physiological changes of gannanzao in response to drought stress， and provide a molecular basis for the study of its drought resistance.9A12BED8-D238-43C5-A21B-7259BF963ED4

Keywords： navel orange; drought stress; photosynthesis; enzymatic activity; malondialdehyde; transcriptome analysis

DOI： 10.3969/j.issn.1000-2561.2022.05.003

柑橘（Citrus L.）是蕓香科、柑橘屬植物，是世界上種植面積和產(chǎn)量最大的水果。我國是柑橘的主要原產(chǎn)國之一，正逐漸向區(qū)域化、商品化、良種化、專業(yè)化、集約化等方向發(fā)展^[1]。以贛南臍橙為主的江西優(yōu)勢柑橘產(chǎn)業(yè)帶，季節(jié)性干旱頻率較高，危害嚴重，引起柑橘產(chǎn)量和品質(zhì)下降，是柑橘生產(chǎn)上亟待解決的問題之一^[2]。水分是柑橘生長過程中的必要條件^[3]。當植物體處于干旱脅迫環(huán)境條件下時，可以通過氣孔關(guān)閉減少水分蒸發(fā)，調(diào)節(jié)可溶性糖含量，以保證植株正常代謝，抵抗逆境帶來的危害^[4]。隨著干旱脅迫程度增加，植物細胞膜過氧化程度增加，丙二醛（malonaldehyde， MDA）、超氧化物歧化酶（superoxide dismutase， SOD）和過氧化物歧化酶（peroxidase， POD）等保護酶可以反應(yīng)植物對干旱脅迫的耐受程度^[5]。不同的植物其生理物質(zhì)的活性和含量在干旱脅迫條件下存在差異，有關(guān)‘贛南早臍橙在此方面的研究尚未見報道?！M南早臍橙（‘Gannan Zaonavel orange， GNZ）是從‘紐荷爾臍橙（‘Newhallnavel orange， NHE）中選育的芽變品種，有其獨特的品質(zhì)，具有轉(zhuǎn)色快、糖積累快、低酸等特性^[6]，比目前主栽品種‘紐荷爾臍橙早熟30?d以上^[7]，‘贛南早臍橙是臍橙主產(chǎn)區(qū)用于品種結(jié)構(gòu)調(diào)整的一個優(yōu)良早熟品種，在臍橙產(chǎn)區(qū)種植面積已超6666.67 hm²。本研究以‘贛南早臍橙與其對照‘紐荷爾臍橙為材料，研究二者在干旱脅迫下的主要生理指標響應(yīng)差異，并進一步分析其抗旱能力及基因表達差異。本研究有助于深入了解‘贛南早臍橙抗旱性，為新品種‘贛南早臍橙的栽培技術(shù)與大面積推廣提供參考，同時為柑橘抗旱遺傳改良提供重要的基因資源。

1 ?材料與方法

1.1 材料

1.1.1 ?植物材料??自然干旱試驗材料‘贛南早（Gannan Zao， GNZ）臍橙和‘紐荷爾（Newhall， NHE）臍橙種植于江西省贛州市信豐縣正平鎮(zhèn)‘贛南早臍橙基地，為4年生枳殼砧臍橙樹，同排種植，每個品種各選3株。盆栽控水干旱試驗材料來自贛南師范大學國家臍橙工程技術(shù)研究中心柑橘種質(zhì)資源圃，是‘贛南早臍橙和‘紐荷爾臍橙枳殼砧一年半生苗，每個品種各20株苗。對自然干旱材料進行光合指標測定，生理及轉(zhuǎn)錄組分析材料均取自盆栽控水干旱試驗材料，生理指標測定材料于干旱處理0、5、10、15、20?d和復(fù)水5?d（F5）時取‘贛南早和‘紐荷爾臍橙成熟枝條葉片，轉(zhuǎn)錄組測序分析材料取干旱處理0、10、20?d的葉片。樣品液氮速凍后，迅速帶回實驗室–80℃保存?zhèn)溆?，以上取樣均設(shè)置3次生物學重復(fù)。

1.1.2 ?干旱處理??自然干旱處理為2019年9—11月田間自然干旱60 d進行光合指標測定，期間贛南地區(qū)連續(xù)無超10 mm降水且試驗用樹未人工施水，復(fù)水5?d后再次測定。盆栽控水干旱處理材料生長環(huán)境一致（盆高38 cm，直徑25 cm，基質(zhì)為園土∶甜葉菊渣∶沙土=1∶1∶1），正常管理6個月后，在處理前正常澆水3?d，最后一次澆透水后進行控水處理，分別處理0、5、10、15、20?d，干旱處理20 d后復(fù)水5 d，并拍照記錄葉片形態(tài)。

1.2 方法

1.2.1 ?生理指標測定??光合作用參數(shù)凈光合速率（P_n）、氣孔導(dǎo)度（G_s）、胞間CO₂濃度（C_i）和蒸騰速率（T_r）使用Li-6400XT便攜式光合儀測定，設(shè)定空氣流速為500?mol/s，CO₂濃度400?μmol/mol。相對電導(dǎo)率、丙二醛及葉片可溶性糖含量測定操作參考《植物生理學實驗指導(dǎo)》^[8]，用電導(dǎo)率儀（雷磁DDS-307，上海）測定相對電導(dǎo)率值。分別用硫代巴比妥酸法、蒽酮比色法檢測并計算丙二醛和可溶性糖含量^[9]，超氧化物歧化酶、過氧化物酶活性測定^[9]采用氮藍四唑法和愈創(chuàng)木酚法。

1.2.2 ?轉(zhuǎn)錄組測序??由美吉生物科技有限公司（上海，中國）對‘贛南早臍橙和‘紐荷爾臍橙干旱處理0、10、20 d葉片，共18個樣品進行轉(zhuǎn)錄組分析，以甜橙基因組為參考基因組進行序列比對分析。參考基因組來源：Citrus_sinensis， http：//citrus.hzau.edu.cn/orange/download/index.php，參考基因組版本：hzau。Illumina平臺通過將測序圖像信號經(jīng)CASAVA堿基識別（base calling）轉(zhuǎn)換成文字信號，并將其以fastq格式儲存起來作為原始數(shù)據(jù)。根據(jù)index序列區(qū)分各個樣本的數(shù)據(jù)，以便進行后續(xù)分析。原始數(shù)據(jù)已上傳NCBI數(shù)據(jù)庫，SubmissionID：SUB10857158，BioProjectID：?PRJNA792482。9A12BED8-D238-43C5-A21B-7259BF963ED4

1.2.3??轉(zhuǎn)錄數(shù)據(jù)質(zhì)量控制??為保證后續(xù)分析的準確性，使用SeqPrep和Sickle軟件處理原始數(shù)據(jù)，去除reads中的接頭序列，將序列末端（3'端）低質(zhì)量（質(zhì)量值小于20）的堿基修剪掉，如剩余序列中仍然有質(zhì)量值小于10的堿基，則將整條序列剔除，去除含N（模塊堿基）的reads，舍棄去adapter及質(zhì)量修剪后長度小于30?bp的序列，并對質(zhì)控后的數(shù)據(jù)進行質(zhì)量評估。

1.2.4??轉(zhuǎn)錄數(shù)據(jù)分析??利用表達定量軟件RSEM分別對轉(zhuǎn)錄本的表達水平進行定量分析，差異基因表達量及篩選，使用DESeq2軟件對樣本進行分析，參數(shù)：P-adjust<0.05 & |log2FC|≥1。采用軟件Goatools對DEGs進行GO富集分析，F(xiàn)isher精確檢驗，用Bonferroni方法對P值進行校正，P值（FDR）<0.05。采用R腳本對差異基因進行KEGG PATHWAY富集分析，F(xiàn)isher精確檢驗，P值（correctedP-value）以0.05為閾值，滿足此條件的KEGG通路定義為在基因集中顯著富集的KEGG通路。采用PlantTFDB 4.0（http：//planttfdb.?cbi.pku.edu.cn/）對差異基因進行轉(zhuǎn)錄因子分析。

1.3??數(shù)據(jù)處理

利用Microsoft Excel 2007軟件對實驗數(shù)據(jù)進行整理、分析和作圖，并用SPSS Statistical 23.0軟件分析差異顯著性。

2 ?結(jié)果與分析

2.1 ?自然干旱脅迫下‘贛南早和‘紐荷爾臍橙葉片形態(tài)及光合指標差異

葉片卷曲程度是植物受干旱脅迫最直觀的表現(xiàn)，從表型觀察來看（圖1），干旱脅迫15 d時，‘紐荷爾臍橙葉片已出現(xiàn)卷曲，但‘贛南早臍橙葉片仍舒展，在干旱脅迫20 d時，‘贛南早臍橙葉片出現(xiàn)輕微卷曲，‘紐荷爾臍橙葉片卷曲程度更加嚴重，復(fù)水后葉片均恢復(fù)舒展。

在自然干旱60 d時測定田間‘贛南早臍橙和‘紐荷爾臍橙葉片凈光合速率（P_n）、氣孔導(dǎo)度（G_s）、胞間CO₂濃度（C_i）和蒸騰速率（T_r），與復(fù)水后光合作用參數(shù)，結(jié)果表明，自然干旱下‘贛南早臍橙的G_s和T_r值極顯著高于‘紐荷爾臍橙，P_n顯著高于‘紐荷爾臍橙。復(fù)水后，‘贛南早臍橙的P_n值極顯著高于‘紐荷爾臍橙，‘贛南早臍橙的T_r值顯著高于‘紐荷爾臍橙，而G_s、C_i值無顯著性差異（表1）。

2.2 干旱脅迫下‘贛南早和‘紐荷爾臍橙葉片電導(dǎo)率、丙二醛及可溶性固形物含量

相對電導(dǎo)率、丙二醛含量等可判斷柑橘干旱脅迫下細胞膜受損傷程度。結(jié)果表明，干旱處理前‘贛南早和‘紐荷爾臍橙葉片相對電導(dǎo)率差異不顯著，干旱處理中二者相對電導(dǎo)率值均呈上升趨勢，且在復(fù)水后下降（圖2A）。在干旱脅迫10 d，‘紐荷爾的相對電導(dǎo)率顯著高于‘贛南早臍橙，干旱脅迫15、20 d和復(fù)水5 d后極顯著高于‘贛南早臍橙。

如圖2B所示，干旱處理下，‘紐荷爾臍橙葉片丙二醛含量呈上升趨勢，而‘贛南早臍橙在干旱脅迫5～10 d時呈下降趨勢，后又迅速積累。干旱脅迫0、10、15?d時‘紐荷爾臍橙的丙二醛含量極顯著高于‘贛南早臍橙，干旱脅迫5、20 d時‘贛南早臍橙丙二醛含量顯著高于‘紐荷爾臍橙，復(fù)水后均下降。

可溶性糖含量檢測結(jié)果表明（圖2C），二者葉片可溶性糖含量在干旱脅迫0～5 d時呈下降趨勢，且‘紐荷爾臍橙在干旱處理15 d時顯著高于‘贛南早臍橙，二者在干旱脅迫15、20 d時含量均增加，復(fù)水后均下降，但‘紐荷爾臍橙葉片可溶性糖含量極顯著高于‘贛南早臍橙。

2.3 干旱脅迫下‘贛南早和‘紐荷爾臍橙葉片保護酶活性

保護酶活性檢測結(jié)果表明（圖3A），在干旱處理中，‘贛南早和‘紐荷爾臍橙葉片SOD活性均呈先降后升趨勢，‘贛南早葉片SOD活性在干旱處理5 d后呈上升趨勢，而‘紐荷爾臍橙在10 d后呈上升趨勢，且均在復(fù)水后下降。

如圖3B所示，在干旱處理0 d時，‘紐荷爾臍橙葉片POD活性顯著高于‘贛南早臍橙，0～10 d時POD活性均下降，10～20 d時POD活性均增強，且干旱處理10、15?d時‘贛南早臍橙葉片POD活性極顯著高于‘紐荷爾臍橙，復(fù)水后，‘贛南早臍橙葉片POD活性仍增強，而‘紐荷爾臍橙葉片POD活性下降。

2.4 測序數(shù)據(jù)質(zhì)量控制

對‘贛南早和‘紐荷爾臍橙在干旱脅迫0、10、20 d下的18個轉(zhuǎn)錄組樣本進行測序（表2），共獲得130.96 Gb clean reads，各樣本過濾后堿基數(shù)均在6.12 Gb及以上，GC含量均在44.50%左右，較穩(wěn)定，Q₂₀≥97.71%，Q₃₀>93.19%，錯誤率在0.025%左右，表明此測序結(jié)果及質(zhì)控效果較好，可用于后續(xù)轉(zhuǎn)錄組分析。9A12BED8-D238-43C5-A21B-7259BF963ED4

2.5 干旱脅迫下‘贛南早臍橙和‘紐荷爾臍橙的差異表達基因

將干旱脅迫0、10、20 d的轉(zhuǎn)錄組測序數(shù)據(jù)進行分析（圖4），在相同脅迫時間點進行比較，‘贛南早和‘紐荷爾臍橙間獲得了大量差異表達基因（DEGs），分別為1266、683、658個，

其中上調(diào)基因數(shù)為350、231、410個，下調(diào)基因數(shù)為916、452、248個。以干旱脅迫0 d為對照，‘贛南早臍橙在脅迫10、20 d時差異基因數(shù)為3942、3304個，‘紐荷爾臍橙在脅迫10、20 d時分別有236、2041個差異表達基因。干旱脅迫10 d和20 d比較，‘贛南早臍橙鑒定出3752個差異表達基因，而‘紐荷爾臍橙差異基因為2108個，以上篩選出的差異基因均為下調(diào)表達更多。

2.6差異表達基因GO富集分析及KEGG富集分析

為闡明差異表達基因行使的基因功能，采用軟件Goatools對差異基因進行GO富集分析，按P-adjust<0.05的條件選取top20的富集結(jié)果，干旱脅迫0 d和10 d差異基因?qū)Ρ确治霭l(fā)現(xiàn)，‘贛南早臍橙2組間的差異基因主要集中在細胞蛋白修飾過程（cellular protein modification process）、蛋白修飾過程（protein modification process）、高分子修飾作用（macromolecule modification）、含磷化合物代謝過程（phosphate-con taining compound metabolic process）、磷代謝過程（phosphorus metabolic process）和小分子代謝過程（small molecule metabolic process）（圖5A），而‘紐荷爾臍橙未見明顯富集。干旱脅迫0?d與20?d差異基因?qū)Ρ劝l(fā)現(xiàn)，‘贛南早和‘紐荷爾臍橙差異基因均集中于生物學過程（biological?process）。此外，‘贛南早臍橙還特異性富集于氧化還原酶（oxidoreductase activity）、蛋白質(zhì)絲氨酸/蘇氨酸激酶活性（protein serine/threonine kinase activity）、磷酸化作用（phosphorylation）和氧己酸代謝過程（oxoacid metabolic process）等通路中（圖5B、5C）。

采用R腳本對樣本差異基因進行KEGG PATHWAY富集分析，功能分類及PATHWAY注釋顯示（圖6），按P-adjust<0.05的條件選取top20的富集結(jié)果，干旱脅迫0 d和10 d差異基因?qū)Ρ确治霭l(fā)現(xiàn)，‘贛南早臍橙差異基因主要集中在植物MAPK信號通路（MAPK signaling pathway-?plant）、植物激素信號轉(zhuǎn)導(dǎo)（plant hormone signal transduction）和苯丙烷生物合成（phenylpropanoid biosynthesis）（圖6A），而‘紐荷爾臍橙主要集中在苯丙烷生物合成、倍半萜和三萜生物合成（sesquiterpenoid and triterpenoid biosynthesis）和內(nèi)質(zhì)網(wǎng)蛋白加工（protein processing in endoplasmic reticulum）通路（圖6C）。干旱脅迫0?d和20?d差異基因?qū)Ρ确治霭l(fā)現(xiàn)，‘贛南早臍橙還富集在

植物-病原相互作用（plant-pathogen interaction）通路中（圖6B），而‘紐荷爾富集在光合作用（photosynthesis）、光合生物的固碳作用（carbon fixation in photosynthetic organisms）和糖酵解/糖異生（glycolysis/glucon eogenesis）途徑（圖6D）。此外，‘贛南早臍橙特異性富集于淀粉及蔗糖通路（starch and sucrose metabolism）和氨基酸及核苷酸糖代謝（amio sugar and nucleotide sugar metabolism）途徑。

2.7轉(zhuǎn)錄因子分析及抗氧化保護酶相關(guān)基因的轉(zhuǎn)錄調(diào)控

植物響應(yīng)干旱脅迫過程中，轉(zhuǎn)錄因子（transcription factor， TFs）發(fā)揮著重要的調(diào)控作用。本研究對4組比較組獲得的差異基因進行了

轉(zhuǎn)錄因子分析，結(jié)果顯示，‘贛南早臍橙和‘紐荷爾臍橙轉(zhuǎn)錄因子在ERF家族、MYB家族、NAC家族、MYB_related家族、WRKY家族、bHLH家族、HB-other家族、HSF家族、B3家族和bZIP家族均有分布。二者比較發(fā)現(xiàn)，‘紐荷爾臍橙特異性分布于LBD（AS2/LOB）家族、ARF家族和AP2家族，‘贛南早臍橙特異性分布于GRAS家族。干旱脅迫0 d和10 d對比分析發(fā)現(xiàn)，‘贛南早臍橙差異基因主要分布于ERF家族（29個）（圖7A），‘紐荷爾臍橙為WRKY家族（7個）（圖7C），干旱脅迫0 d和20 d對比分析表明，‘贛南早臍橙差異基因主要分布于MYB家族（23個）（圖7B），‘紐荷爾臍橙為bHLH家族（17個）（圖7D）。

根據(jù)差異基因富集分析及轉(zhuǎn)錄因子分析結(jié)果，從轉(zhuǎn)錄組差異表達基因中挖掘出30個編碼抗氧化保護酶基因（表3），過氧化物酶（POD）基因Cs8g17370、Cs7g05450、Cs2g25450和Cs2g03110等上調(diào)表達，Cs3g21730、Cs6g20170、Cs7g19270、Cs1g18600和Cs9g05140等下調(diào)表達。編碼抗壞血酸過氧化物酶（aseorbateperoxidase， APX）基因Cs8g17370和Cs7g05450上調(diào)表達，而Cs7g05430下調(diào)表達。編碼超氧化物歧化酶（SOD）基因Cs8g15520和過氧化氫酶（catalase，?CAT）基因Cs3g27280均下調(diào)表達。9A12BED8-D238-43C5-A21B-7259BF963ED4

3 ?討論

凈光合速率反應(yīng)植物的光合能力的大小^[10]，隨著凈光合速率增大，植物消耗CO₂量增大，胞間CO₂減少，植物就通過增大氣孔導(dǎo)度獲取更多的CO₂^[11]。在自然干旱條件下，‘贛南早臍橙凈光合速率及氣孔導(dǎo)度均高于‘紐荷爾臍橙，說明‘贛南早臍橙在干旱脅迫下的光合作用能力強于‘紐荷爾臍橙。通常細胞膜能維持細胞內(nèi)穩(wěn)態(tài)，但干旱脅迫會導(dǎo)致植物細胞膜破壞導(dǎo)致大量電解質(zhì)外滲，因此可以根據(jù)植物在干旱脅迫下的相對電導(dǎo)率來評估其抗旱性。同時，當植物在逆境條件下，往往發(fā)生膜脂過氧化作用，丙二醛是其產(chǎn)物之一，通常利用它表示植物對逆境反應(yīng)的強弱^[12]。本研究結(jié)果表明，隨著干旱脅迫天數(shù)增加，‘贛南早和‘紐荷爾臍橙葉片相對電導(dǎo)率值均上升，且‘贛南早臍橙的相對電導(dǎo)率值

均低于‘紐荷爾臍橙，說明干旱脅對‘紐荷爾臍橙葉片細胞損傷程度大于‘贛南早臍橙。二者葉片丙二醛含量也均呈上升趨勢，在脅迫10、15?d時‘紐荷爾臍橙葉片丙二醛含量高于‘贛南早臍橙，說明‘紐荷爾臍橙葉片膜脂過氧化程度比‘贛南早臍橙高，細胞膜受損更嚴重。

SOD、POD是植物中防止膜脂過氧化的重要酶，因此SOD和POD活性能較好地反映植物對逆境的耐受能力^[13-14]，植物為避免細胞遭到氧化損傷，通過自我調(diào)節(jié)啟動一系列抗氧化酶清除過量的活性氧物質(zhì)，超氧化物歧化酶（SOD）是第一個參與植物細胞內(nèi)活性氧清除的抗氧化酶，當SOD催化超氧根陰離子發(fā)生歧化反應(yīng)后將產(chǎn)生H₂O₂和O₂，由POD和CAT等保護酶負責進一步將H₂O₂分解成H₂O和O₂^[15]。本研究發(fā)現(xiàn)，在干旱脅迫5 d時‘贛南早臍橙葉片SOD活性增強，而‘紐荷爾臍橙葉片SOD酶活性在干旱脅迫10 d時增強，且干旱脅迫10 d后‘贛南早臍橙葉片SOD活性顯著高于‘紐荷爾臍橙，而二者葉片POD活性變化較為一致，干旱脅迫10 d后活性均上升，且‘贛南早臍橙葉片POD活性極顯著高于‘紐荷爾臍橙，說明與‘紐荷爾臍橙相比，‘贛南早臍橙對干旱脅迫環(huán)境刺激更敏感，干旱脅迫加劇時，通過提高保護酶活性來抵御干旱逆境能力更強。

轉(zhuǎn)錄因子通過調(diào)控靶基因表達，在植物抵抗生物及非生物脅迫中起重要作用^[16]。‘贛南早臍橙和‘紐荷爾臍橙轉(zhuǎn)錄因子分析比較發(fā)現(xiàn)，‘贛南早臍橙特異性分布于GRAS家族，GRAS家族基因具有促進分生組織發(fā)育，光敏信號傳導(dǎo)和GA代謝調(diào)節(jié)等作用^[17]，說明‘贛南早臍橙抗干旱能力可能還與植物激素類代謝有關(guān)。干旱脅迫0 d和10 d對比發(fā)現(xiàn)，‘贛南早臍橙有29個差異基因分布于ERF家族，‘紐荷爾臍橙有7個差異基因分布于WRKY家族。干旱脅迫0 d和20 d對比分析表明，‘贛南早臍橙有23個差異基因分布于MYB家族，‘紐荷爾臍橙有17個差異基因分布于bHLH家族。ERF家族和MYB家族均為復(fù)雜的轉(zhuǎn)錄因子家族，不僅在種子萌發(fā)等植物發(fā)育過程中起重要作用，而且一些基因已被證明在干旱、高溫等非生物脅迫中起重要作用^[18-21]，而WRKY家族通過調(diào)節(jié)激素類物質(zhì)代謝幫助植物抵抗逆境^[22]。bHLH家族則以結(jié)合作用，與其他蛋白互作響應(yīng)干旱脅迫^[23]，說明干旱脅迫下ERF家族、MYB家族、WRKY家族和bHLH家族在調(diào)控網(wǎng)絡(luò)中起重要作用，‘贛南早臍橙和‘紐荷爾臍橙對干旱脅迫耐受能力不同可能與轉(zhuǎn)錄因子調(diào)控有關(guān)，有待進一步發(fā)掘與考證。

參考文獻

[1]?馬小煥， 賴九江. 贛南臍橙產(chǎn)業(yè)現(xiàn)狀淺析及對策研究[J]. 現(xiàn)代園藝， 2021， 44（17）：?49-51.MA X H， LAI J J. Analysis and countermeasures of navel orange industry in southern Jiangxi[J]. Modern Horticulture， 2021，?44（17）： 49-51. （in Chinese）

[2]?樊衛(wèi)國， 李慶宏， 吳素芳. 長期干旱環(huán)境對柑橘生長及養(yǎng)分吸收和相關(guān)生理的影響[J]. 中國生態(tài)農(nóng)業(yè)學報， 2012（11）： 76-85.FAN W G， LI Q H， WU S F. Effects of long-term arid environment on citrus growth， nutrient absorption and related physiology[J]. Chinese Journal of Ecological Agriculture， 2012（11）： 76-85. （in Chinese）

[3]?王鳳麗. 持續(xù)高溫干旱對果樹的影響及高產(chǎn)栽培管理措施[J]. 現(xiàn)代農(nóng)業(yè)科技， 2020（24）： 67-68.WANG F L. Effects of continuous high temperature and drought on fruit trees and high-yield cultivation and management measures[J]. Modern Agricultural Science and Technology， 2020（24）： 67-68. （in Chinese）

[4]?楊建瑩， 霍治國， 王培娟， 鄔定榮， 毛紅丹， 孔??瑞. 中國北方蘋果干旱等級指標構(gòu)建及危險性評價[J]. 應(yīng)用氣象學報， 2021， 32（1）： 25-37.YANG J Y， HUO Z G， WANG P J， WU D R， MAO H D， KONG R. Construction of drought grade index and risk assessment of apple in northern China[J]. Journal of Applied Meteorology， 2021，?32（1）： 25-37. （in Chinese）9A12BED8-D238-43C5-A21B-7259BF963ED4

[5]?姜露露， 于??坤， 劉冬冬， 王軍武， 包興成， 鄭??重. 干旱脅迫對無花果葉片生理特性的影響[J]. 中國農(nóng)學通報， 2021， 37（28）： 62-67.JIANG L L， YU K， LIU D D， WANG J W， BAO X C， ZHENG Z. Effects of drought stress on physiological characteristics of fig leaves[J]. Chinese Agronomy Bulletin， 2021，?37（28）： 62-67. （in Chinese）

[6]?陳健美， 謝麗紅， 周??娟， 江小美， 鐘八蓮， 李淑惠， 楊斌華， 喻芳琴.‘贛南早臍橙早熟性狀回復(fù)型突變體的生理與轉(zhuǎn)錄組分析[J]. 果樹學報， 2019， 36（4）： 401-409.CHEN J M， XIE L H， ZHOU J， JIANG X M， ZHONG B L， LI S H， YANG B H， YU F Q. Physiological and transcriptome analysis of revertant mutant of early maturing trait of ‘Gannan Zao Navel Orange[J]. Journal of fruit trees， 2019，?36（4）： 401-409. （in Chinese）

[7]?鐘八蓮， 賴曉樺， 楊斌華， 米蘭芳， 謝上海， 黃彩英， 楊文俠， 張??湟. 紐荷爾臍橙芽變早熟品種—贛南早臍橙[J].中國南方樹， 2013， 42（2）：48-51.ZHONG B L， LAI X H， YANG B H，?MI L F， XIE S H， HUANG C Y， YANG W X， ZHANG H. Early maturing variety of ‘Newhall Navel Orange-Gannan early navel orange[J]. Fruit Trees in Southern China， 2013，?42（2）： 48-51. （in Chinese）

[8]?高俊鳳. 植物生理學實驗指導(dǎo)[M]. 北京： 高等教育出版社， 2006：?74-77， 144-148， 210-219.GAO J F. Experimental guidance of plant physiology[M]. Beijing： Higher Education Press， 2006： 74-77， 144-148， 210-219. （in Chinese）

[9]?李合生. 植物生理生化實驗原理和技術(shù)[M]. 北京： 高等教育出版社， 2000： 260-261.LI H S. Principles and techniques of plant physiological and biochemical experiments[M]. Beijing： Higher Education Press， 2000： 260-261. （in Chinese）

• 劉懷年， 鄧曉建， 李??平. 水稻品種資源光合速率研究[J]. 四川農(nóng)業(yè)大學報，?2007， 25（4）：?379-383.LIU H N， DENG X J， LI P. Study on photosynthetic rate of rice variety resources[J]. Journal of Sichuan Agricultural University， 2007，?25（4）： 379-383. （in Chinese）
• 張效星， 樊??毅， 賈??悅， 崔寧博， 趙??璐， 胡笑濤， 龔道枝. 水分虧缺對滴灌柑橘光合和產(chǎn)量及水分利用效率的影響[J].農(nóng)業(yè)工程學報， 2018， 34（3）： 143-150.ZHANG X X， FAN Y， JIA Y ， CUI N B， ZHAO L， HU X T， GONG D Z. Effects of water deficit on photosynthesis， yield and water use efficiency of Citrus under drip irrigation[J]. Journal of Agricultural Engineering， 2018，?34（3）： 143-150. （in Chinese）
• 王海珍， 梁宗鎖， 韓蕊蓮， 韓 ?路. 土壤干旱對黃土高原鄉(xiāng)土樹種水分代謝與滲透調(diào)節(jié)物質(zhì)的影響[J]. 西北植物學報， 2004（10）： 1822-1827.WANG H Z， LIANG Z S， HAN R L， HAN L. Effects of soil drought on water metabolism and osmoregulation substances of native tree species in the Loess Plateau[J]. Journal of Northwest Botany， 2004（10）： 1822-1827. （in Chinese）
• 張珊珊， 康洪梅， 楊文忠. 干旱脅迫下水楊酸浸種對云南藍果樹幼苗生理響應(yīng)的影響[J]. 東北林業(yè)大學學報， 2016， 44（9）： 34-39.ZHANG S S， KANG H M， YANG W Z. Effects of Seed soaking with salicylic acid under drought stress on physiological response of Yunnan blue fruit tree seedlings[J]. Journal of Northeast Forestry University， 2016，?44（9）： 34-39. （in Chinese）
• 徐迎春， 陳尚武. 水分脅迫期間及脅迫解除后蘋果樹源葉碳同化物代謝規(guī)律的研究[J]. 果樹學報， 2001， 18（1）： 1-6.XU Y C， CHEN S W. Study on the metabolism of carbon assimilates in source leaves of apple trees during and after water stress[J]. Journal of Fruit Trees， 2001， 18（1）： 1-6. （in Chinese）
• 尹永強， 胡建斌， 鄧明軍. 植物葉片抗氧化系統(tǒng)及其對逆境脅迫的響應(yīng)研究進展[J]. 中國農(nóng)學通報， 2007（1）：?115-120.YIN Y Q， HU J B， DENG M J. Research progress on antioxidant system of plant leaves and its response to stress[J]. Chinese Agronomy Bulletin， 2007（1）： 115-120. （in Chinese）

[16]?WANG Y， XU W， CHEN Z，?HAN B，?HAQUE M E，?LIU A Z. Gene structure， expression pattern and interaction of Nuclear Factor-Y family in castor bean （Ricinus communis）[J]. Planta，?2018， 247（3）： 559-572.

[17]?李桂英. 植物GRAS家族轉(zhuǎn)錄因子的研究現(xiàn)狀[J]. 安徽農(nóng)業(yè)科學， 2014（14）： 4207-4210.LI G Y. Research status of plant GRAS family transcription factors[J]. Anhui Agricultural Science， 2014（14）： 4207-4210. （in Chinese）

[18]?TSUBASA S J， LING Y. ERF gene clusters： working together to regulate metabolism[J]. Trends Plant Science， 2021， 26： 23–32.

[19]?崔喜艷， 張瑩瑩， 周??瑩. 植物響應(yīng)干旱脅迫轉(zhuǎn)錄因子研究進展[J/OL]. 吉林農(nóng)業(yè)大學學報： 1-7[2021-11-27].?http：//?kns.cnki.net/kcms/detail/22.1100.S.20210412.1425.006.html.CUI X Y， ZHANG Y Y， ZHOU Y. Research progress of plant transcription factors in response to drought stress[J/OL]. Journal of Jilin Agricultural University： 1-7[2021-11-27].?http：//kns.cnki.net/kcms/detail/22.1100.S.20210412.1425.006.html. （in Chinese）

[20]?李嘉哲， 田廣宇， 張敬鑫， 茶新有， 楊建坤， 董晨肖， 張榮沭. 山新楊MYB4基因在脅迫條件下的組織表達分析[J/OL]. 分子植物育種： 1-9[2021-12-04].?http：//kns.cnki.net/?kcms/detail/46.1068.S.20211111.1502.004.html.LI J Z， TIAN G Y， ZHANG J X， CHA X Y， YANG J K， DONG C X， ZHANG R S. Tissue expression analysis of MYB4 gene of Shanxin poplar under stress[J/OL]. Molecular Plant Breeding： 1-9[2021-12-04]. http：//kns.cnki.net/?kcms/detail/46.1068.S.20211111.1502.004.html. （in Chinese）

• 柯玉潔， 陳明堃， 馬山虎， 歐??悅， 王??藝， 鄭清冬， 劉仲健， 艾??葉. 蘭科植物MYB轉(zhuǎn)錄因子研究進展[J/OL]. 園藝學報： 1-10[2021-12-04]. https：//doi.org/10.16420/j.issn.?0513-353x.2020-0843.KE Y J， CHEN M K， MA S H， OU Y， WANG Y， ZHENG Q D， LIU Z J， AI Y. Research progress of MYB transcription factors in Orchidaceae[J/OL]. Journal of Horticulture： 1-10 [2021-12-04]. https：//doi.org/10.16420/j.issn.0513-353x.2020-?0843. （in Chinese）
• 陳林英， 李佳佳， 王??博， 杜婉清， 高夢雪， 劉??慧， 檀淑琴， 邱麗娟， 王曉波. WRKY轉(zhuǎn)錄因子在大豆響應(yīng)生物和非生物脅迫中的功能研究進展[J/OL]. 植物遺傳資源學報： 1-11[2021-12-04]. https：//doi.org/10.13430/j.cnki.jpgr.2021?1011001.CHEN L Y， LI J J， WANG B， DU W Q， GAO M X， LIU H， TAN S Q， QIU L J， WANG X B. Research progress on the function of WRKY transcription factor in soybean response to biological and abiotic stress[J/OL]. Journal of Plant Genetic Resources： 1-11[2021-12-04]. https：//doi.org/10.13430/?j.cnki.jpgr.20211011001. （in Chinese）
• 劉曉月， 王文生， 傅彬英. 植物bHLH轉(zhuǎn)錄因子家族的功能研究進展[J]. 生物技術(shù)進展， 2011， 1（6）： 391-397.LIU X Y， WANG W S， FU B Y. Research progress on the function of plant bHLH transcription factor family[J]. Biotechnology progress， 2011，?1（6）： 391-397. （in Chinese）