劉卓毅 于文菲 蔡文麗 劉子珊 張雨 袁媛 伍炳華 呂美玲

摘 ?要：幾丁質酶（chitinase， EC 3.2.1.14）是一種降解幾丁質的糖苷酶，在植物應對非生物和生物脅迫中起重要作用。本研究對辣椒幾丁質酶類（Capsicum annuum chitinase-like， CaCTL）基因家族成員進行了全基因組注釋、進化分析和基因表達模式分析。從最新的辣椒基因組數(shù)據庫中鑒定出31個CaCTL成員。根據系統(tǒng)發(fā)育進化樹將這些成員分為GH18和GH19亞家族，并進一步分為5個亞類（Ⅰ～Ⅴ類）。保守基序分析結果表明，每個亞類中的成員存在功能相關性。對辣椒、矮牽牛以及番茄的CTL基因家族成員同源性分析表明，在矮牽牛和番茄中分別有7個和11個辣椒CaCTL的同源基因。通過qRT-PCR分析發(fā)現(xiàn)，在正常環(huán)境生長的辣椒植株中，64.2%的GH18亞家族的CaCTL成員表達水平很低，有64.7%的GH19亞家族成員在不同組織中存在差異表達。順式作用元件分析表明，CaCTL成員啟動子區(qū)域存在許多激素反應以及生物和非生物應激相關元件。當植物遭受不同的非生物脅迫時，qRT-PCR結果顯示，多數(shù)CaCTL成員表達上調。當植物遭受高溫干旱條件時，多個CaCTL成員的表達上調明顯。上述結果豐富了CTL基因家族進化的研究，為探索CaCTL成員的功能提供了數(shù)據基礎。

關鍵詞：辣椒;幾丁質酶類基因;非生物脅迫;進化分析;表達分析

中圖分類號：S641.3 ? ? ?文獻標識碼：A

Abstract： Chitinase （EC 3.2.1.14） is chitin-degrading glycosidase， which plays important roles in the abiotic and biotic defense of plants. In this study， we conducted whole-genome annotation， molecular evolution， and gene expression analyses on the Capsicum annuum chitinase-like （CaCTL） gene family members. Thirty-one CaCTL genes （CaCTLs） were identified from the latest Capsicum genome database. According to the phylogenetic analysis， the genes were divided into GH18 and GH19 subgroups and further subdivided into five classes （Class I to Class V）. Conserved motif arrangements indicated the functional relevance within each group. Homology analysis of CTL gene family members in pepper， petunia and tomato indicated that there were 7 and 11 homologous genes of CaCTL members in petunia and tomato， respectively. By qRT-PCR analysis， in the normal growing plants， we found that 64.2% GH18 subgroup CaCTLs showed a very low expression level， 64.7% GH19 subgroup CaCTLs were differentially expressed in different tissues. Cis-acting element analysis showed that many hormone responses， as well as biotic and abiotic stress elements were found in the promoters of CaCTLs. Correspondingly， qRT-PCR showed that many CaCTLs had upre gulated expression profiles when the plants suffered different abiotic stress conditions， especially when plants are subjected to high temperature and drought conditions. The results would enrich the evolutionary process of the CTL gene family and provide a foundation for the exploration of the function of CaCTLs.

Keywords： Capsicum annuum L.; chitinase-like gene; abiotic stress; evolution analysis; expression analysis

DOI： 10.3969/j.issn.1000-2561.2021.11.008

幾丁質酶（chitinase，EC.3.2.1.14）是一種由單基因編碼的降解幾丁質的糖苷酶，主要水解幾丁質多聚體中的β-1，4糖苷鍵，產生N-乙酰氨基葡萄糖寡聚體。幾丁質酶廣泛存在于自然界的各類生物中。細菌可以通過合成的幾丁質酶分解幾丁質獲得營養(yǎng)維持生存;真菌中的幾丁質酶在其細胞分裂和形態(tài)建成等方面發(fā)揮作用;昆蟲等節(jié)肢動物中的幾丁質酶在其蛻皮等生長發(fā)育中發(fā)揮重要的作用[1];人體內的類幾丁質酶大多與免疫相關[2]。植物中的幾丁質酶表達具有組織特異性，參與植物的營養(yǎng)生長[3]和有性生殖[4]的調控過程。幾丁質酶在植物抵御疾病以及環(huán)境脅迫方面也發(fā)揮重要作用。研究表明，當植物受到干旱[5]、低溫[6]、病蟲害[7-8]等非生物和生物脅迫時，會引起體內幾丁質酶含量的增加。

根據氨基酸序列的差異以及結構域的有無，植物幾丁質酶可以分為糖苷水解酶家族18（glycoside hydrolase family 18，GH18）和糖苷水解酶家族19（glycoside hydrolase family 19，GH19）2個亞家族。GH18亞家族成員廣泛存在于細菌、真菌、植物、哺乳動物、昆蟲和病毒等生物體中，而GH19亞家族成員僅出現(xiàn)在植物以及鏈霉菌屬中[9]。Cao等[10]在番茄中鑒定了43個SlCTL（Solanum lycopersicum chitinase-like）成員，有11個屬于GH18亞家族，有32個屬于GH19亞家族，發(fā)現(xiàn)脅迫處理以及外源激素處理會引起多數(shù)成員的響應。GH19亞族成員具有Chitin-bind domain（CBD）結構。CBD可以與幾丁質結合，加速或加劇植株對侵入病蟲菌的識別與反應，增強植物體對環(huán)境脅迫的抵抗能力[11-12]。Ali等[13]從辣椒中鑒定了16個具有CBD結構的幾丁質酶類基因（Capsicum annuum chitinase-like，CaCTL），包括CM334數(shù)據庫中的8個和Zunla-1數(shù)據庫中的8個，發(fā)現(xiàn)這些CaCTL成員參與逆境防御反應和信號轉導，在辣椒的生長調節(jié)以及脅迫響應中發(fā)揮至關重要的作用。

辣椒是一種重要的經濟作物，具有極高的生產價值。在生長過程中極易遭受炭疽病、灰霉病、軟腐病等各類生物病害和異常溫度、水量、pH、環(huán)境等非生物因素的傷害，影響其產量與品質，導致經濟價值下降。本研究通過序列檢索和比對，在辣椒CM334數(shù)據庫中檢索到了31個CaCTL成員（包括Ali等[13]從辣椒CM334數(shù)據庫中檢索到的8個含有CBD結構的CaCTL成員），分析了CaCTL的進化關系、基因結構、組織特異表達以及響應環(huán)境脅迫的表達特征等。

1 ?材料與方法

1.1 ?CaCTL基因家族成員的序列鑒定

根據已發(fā)表文獻[10， 14]，在擬南芥（Arabidopsis thaliana）基因組數(shù)據庫（TAIR，http：//arabidopsis. org/）和茄科基因組數(shù)據庫（Sol genomics network，https：//solgenomics.net/）中分別下載到25條擬南芥和43條番茄的CTL成員的氨基酸序列。將上述序列在茄科基因組網站的Capsicum annuum L. v CM344數(shù)據庫中進行蛋白BLAST比對，收集比對結果中有“Chitinase”或者“Chitinase-like”注釋的氨基酸序列。為了防止遺漏，同時在Capsicum annuum L. v CM344蛋白注釋數(shù)據庫中以“Chitinase”為關鍵詞直接檢索，收集檢索結果中的氨基酸序列。將上述2次獲得的序列結果進行整合，去除冗余序列。使用NCBI數(shù)據庫中的“CD-Search”功能分析非冗余蛋白序列中的功能域。根據上述分析結果，確定最終的CaCTL基因家族成員。

1.2 ?CaCTL基因家族成員的序列比對、進化樹構建和理化性質預測

在NCBI中搜索獲得擬南芥（A. thaliana）、人類（Homo sapiens）、歐洲云杉（Picea abies）、玉米（Zea mays）、羅晃子（Tamarindus indica）和釀酒酵母（Saccharomyces cerevisiae）6個物種中鑒定的CTL的氨基酸序列。利用ClustalW軟件進行CaCTL成員與上述蛋白的氨基酸多重序列比對，使用MEGA 7軟件的NJ（Neighbor Joining）法構建系統(tǒng)發(fā)育樹，并使用EvolView（https：//www. evolgenius.info）在線美化。通過ExPASY網站（https： //web.expasy.org/protparam/）中的檢索計算功能獲得CaCTL成員的分子量、等電點信息。利用CELLO在線工具（CELLO v.2.5，http：//cello.life. nctu.edu.tw/）預測蛋白的亞細胞定位。

1.3 ?基因結構、保守基序、順式作用元件和同源性分析

利用MEME 5.0.5在線工具（http：//meme-suite. org/tools/meme）分析CaCTL基因家族成員的保守基序，最大檢索參數(shù)設置為16，其他數(shù)值默認。為了增加可信度，運用保守域數(shù)據庫（CDD，http：//www.ncbi.nlm.nih.gov/cdd/）進一步檢查基因家族成員序列的保守結構域。進而利用TBtools工具繪制CaCTL家族成員的基因結構（外顯子-內含子）圖。

利用TBtools軟件在茄科基因組網站（Sol genomics network，https：//solgenomics.net/）提取辣椒基因組文件中各家族成員起始密碼子上游2 kb的序列。將獲得的序列提交到PlantCare網站（http：//bioinformatics.psb.ugent.be/webtools/plantcare/html/）進行順勢作用元件預測。對預測的順勢元件信息進行整理篩選后，使用TBtools進行可視化作圖。

為分析CTL基因家族在茄科植物中的分布及同源關系，從茄科基因組數(shù)據庫中下載矮牽牛（ftp：//ftp.solgenomics.net/genomes/Petunia_axillaris）、辣椒（ftp：//ftp.solgenomics.net/genomes/Capsi-cum_annuum/C.annuum_cvCM334）和番茄（ftp：// ftp.solgenomics.net/tomato_genome/Fla.8924）的基因組序列文件和注釋文件（矮牽牛未組裝到染色體水平），使用TBtools軟件進行基因同源性分析并對結果進行可視化。

1.4 ?植物材料與脅迫處理

植物材料辣椒種植于福建農林大學園藝學院人工氣候室。培養(yǎng)光周期為14 h光照/10 h黑暗，溫度為（25±1）℃，濕度為70%。在生殖生長期取辣椒的根、莖、葉、花蕾、開放花和漿果，保存于–80 ℃冰箱。

植物材料的脅迫處理試驗：首先用75%的乙醇清洗辣椒種子30 s，進而用無菌水沖洗3次，然后用2%次氯酸鈉溶液消毒5 min，再次用無菌水沖洗3次。將上述處理過的種子播種于MS液體培養(yǎng)基中，置于在14 h光照/10 h黑暗環(huán)境下震蕩培養(yǎng)，生長溫度設置25 ℃。待辣椒幼苗子葉展開后，將幼苗分為3個處理組，第1組置于（4±1）℃條件下3 h進行低溫脅迫處理，第2組置于（42±1）℃下3 h進行高溫脅迫處理，第3組置于無塵濾紙上24 h進行干旱脅迫處理。以置于正常生長條件下培養(yǎng)的幼苗作為對照。收集上述處理后的植物材料，保存于–80 ℃冰箱。

1.5 ?qRT-PCR分析

取辣椒根、莖、葉、花蕾、開放花、漿果和受脅迫處理的幼苗材料各100 mg，使用Trizol試劑盒（全式金生物，北京）提取總RNA。分別取上述組織的1 μg總RNA，使用PrimeScript RT試劑盒（全式金生物，北京）反轉錄成一鏈cDNA。

利用NCBI網站的引物設計工具（http：// www.ncbi.nlm.nih.gov/tools/primer-blast/）設計辣椒CaCTL基因家族成員的qRT-PCR引物。根據文獻[15]，選擇CaUBI為內參基因（F：5-TGTCC ATCTGCTCTCTGTTG-3，R：5-CACCCCAAGCA CAATAAGAC-3）。用TransStart Top Green qPCR SuperMix試劑盒（全式金生物，北京），在CFX96實時系統(tǒng)（Bio-Rad，美國）實時熒光定量PCR儀上進行各基因表達水平的檢測。每一反應總體積為15 μL。qRT-PCR反應條件如下：94 ℃反應30 s后，進行42個循環(huán)。每個循環(huán)為94 ℃反應5 s，60 ℃反應15 s和72 ℃反應10 s。試驗設置3次生物學重復，每次生物學重復開展3次技術重復，采用2-CT方法計算相對表達水平。隨后利用TBtools軟件生成熱圖，利用GraphPad Prism 8軟件生成柱狀圖。

2 ?結果與分析

2.1 ?CaCTL基因家族成員的鑒定與系統(tǒng)發(fā)育分析

通過序列比對和功能域分析，發(fā)現(xiàn)辣椒基因組中存在31個被注釋的有水解幾丁質功能或具有幾丁質酶類功能的CaCTL成員（表1）。序列分析表明，CaCTL成員氨基酸的長度范圍為55 AA（CA10g09860）～436 AA（CA05g11510），蛋白質分子量為6.342～49.598 kDa，等電點（pI）為4.64～9.24。通過亞細胞定位分析，發(fā)現(xiàn)CaCTL成員大多定位于胞外基質中，少數(shù)定位在溶酶體、線粒體或質膜上。

系統(tǒng)發(fā)育分析結果表明（圖1），CaCTL基因家族成員分別有14個（45.1%）和17個（54.9%）屬于GH18和GH19亞家族。本研究將CaCTL家族成員進一步分為5類。其中Ⅲ類和Ⅴ類屬于GH18亞家族，分別包含9個和5個CaCTL成員，都包含有GH18結構域;Ⅰ類、Ⅱ類和Ⅳ類屬于GH19亞家族，分別包含2個、12個和3個CaCTL成員，都包含有GH19結構域。

2.2 ?CaCTL基因家族成員的基因結構、保守結構域以及順式元件分析

外顯子-內含子結構分析結果顯示（圖2），GH18和GH19亞家族間基因的結構和內含子的數(shù)目存在明顯差異。GH18亞家族成員多數(shù)（85.7%）不含有內含子，而GH19亞家族成員多數(shù)（58.9%）含有內含子。在31個CaCTL成員中有5個具有較長的內含子結構，包括CA12g20880、CA00g51890、CA08g10220、CA5g11510、CA04g08220;有19個成員（61.2%）不具備內含子，有3個成員（9.7%）具有2個內含子，1個成員（3.2%）具有3個內含子，1個成員CA05g11510（3.2%）具有7個內含子。

通過NCBI網站中的Pfam工具，檢索分析16個預測可能與GH18或者GH19相關的氨基酸保守基序。分析結果表明（圖3），同一亞家族的CaCTL成員間存在功能的相關性。在GH19亞家族中，Ⅰ類和Ⅳ類成員都含有Glyco_hydro_19 motif，Ⅱ類中只有3個成員含有Glyco_hydro_19 motif，其余成員均含有Chitin_bind motif （除CA08g10220），并且Chitin_bind motif僅分布于GH19亞家族中。在包含Ⅲ類和Ⅴ類的GH18亞家族中，除了CA05g11510沒有被檢索到motif外，其他的CaCTL成員廣泛分布有Glyco_hydro_18 motif。但是利用NCBI的CD-Search工具進行檢索時，發(fā)現(xiàn)CA05g11510具有GH18結構域。

啟動子區(qū)域往往存在許多防御相關的響應元件。對起始密碼子上游2 kb區(qū)域的序列進行預測，發(fā)現(xiàn)CaCTL成員啟動子區(qū)域含有的順式元件涉及響應植物激素和非生物脅迫等多個方面（圖4）。主要包括5種激素類響應元件，分別為茉莉酸響應元件（CGTCA/TGACG-motif）、赤霉素響應元件（TATC/P/GARE-motif）、脫落酸響應元件（ABRE）、水楊酸響應元件（TC A-element）和生長素響應元件（TGA-element/AuxRR-core）。CGTCA/TGACG-motif存在于22個CaCTL成員的啟動子區(qū)域，其中CA07g02470的啟動子區(qū)域中最多，含有12個;TATC/P/GARE-motif存在于15個CaCTL成員的啟動子區(qū)域，其中CA01g20040的啟動子區(qū)域中最多，含有4個;ABRE存在于23個CaCTL成員的啟動子區(qū)域，CA04g08220和CA07g02460的啟動子區(qū)域中最多，分別含有5個;TCA-element存在于10個CaCTL成員的啟動子區(qū)域，除了CA08g07860和CA11g13780的啟動子區(qū)域分別含有3個，其余基因分別含有1個;TGA-element/AuxRR-core存在于14個CaCTL成員的啟動子區(qū)域，CA08g10220的啟動子區(qū)域中最多，含有3個。干旱誘導響應元件（MBS）、低溫響應元件（LTR）、防御與脅迫響應元件（TC-rich repeats）以及晝夜節(jié)律響應元件（circadian）也分布在該家族的許多成員的啟動子區(qū)域。MBS存在于15個CaCTL成員的啟動子區(qū)域，CA03g30180的啟動子區(qū)域中存在5個;LTR存在于9個CaCTL成員的啟動子區(qū)域，除了CA11g13780的啟動子區(qū)域中有2個以外，其他8個基因中都只含有1個該元件;TC-rich repeats存在于12個CaCTL成員的啟動子區(qū)域，CA03g02990、CA03g30180、CA12g08860的啟動子區(qū)域中數(shù)量最多，各含有3個。其他順式元件，比如創(chuàng)傷響應元件（WUN- motif）僅存在與CA00g97030的啟動子區(qū)域，參與玉米蛋白代謝調控的元件（O2-site）在CA08g15980的啟動子區(qū)域中最多，含有3個。MYBHv1結合位點（CCAAT-box）以及AT-rich DNA結合蛋白（ATBP-1）的結合位點（AT-rich element）分別存在于9個和5個CaCTL成員的啟動子區(qū)域中。

2.3 ?CaCTL基因家族成員染色體定位以及基因復制和進化分析

染色體定位分析結果表明（圖5），在辣椒的12條染色體中，除第2號染色體外，其他染色體上均有CaCTL基因的分布。有29個CaCTL基因定位在11條染色體上，2個CaCTL基因分別定位在2條scaffold上（CA00g51890定位在scaffold5575上，CA00g97030定位在scaffold1231上）。第3號和第7號染色體上分布最多，各含有6個CaCTL成員，其次是第8號染色體，含有4個CaCTL成員。在第7號、第8號和第11號染色體上分別含有1個串聯(lián)重復基因對（CA07g02440和CA07g02460、CA08g07860和CA08g07870、CA11g13780和CA11g13790）。

為了分析不同物種中CTL的進化關系，本研究對茄科植物辣椒、矮牽牛和番茄的CTL基因家族成員的同源性進行了分析。結果表明（圖6），矮牽牛和番茄基因組中有多個CTLs與辣椒的CaCTL成員存在同源性。辣椒基因組中分別有7個CaCTLs（CA08g07860、CA06g08030、CA07g 02290、CA07g02440、CA11g13780、CA12g02880、CA03g02990）和矮牽牛CTL成員存在同源性;有11個CaCTLs（CA08g07860、CA04g14160、CA05g11 510、CA06g08030、CA11g13780、CA12g20880、CA03g02990、CA03g31390、CA03g30170、CA10g09 860和CA10g09850）和番茄CTL成員存在同源性。

圖中灰色連線代表2個基因組具有同源關系的成員，藍色連線表示2個基因組中具有同源關系的CTL基因家族成員。

The gray lines in the figure represent the members of the two genomes that have a homologous relationship， and the blue. lines indicate the members of the CTL gene family that have a homologous relationship in the two genomes.

2.4 ?CaCTL基因家族成員的表達特征

通過qRT-PCR試驗檢測CaCTL成員在辣椒6個組織（根、莖、葉、花蕾、開放花和漿果）中的表達特征。結果表明，GH18和GH19兩個亞家族CaCTL成員之間的表達差異非常明顯（圖7）。在GH19亞家族的17個CaCTL成員中，有6個基因（CA06g08030、CA00g97030、CA10g09850、CA03g30180、CA03g31390和CA03g30190）在所有組織中低表達或者不表達，而其他11個基因（64.7%）在不同組織中呈現(xiàn)差異表達。CA10g 09860、CA07g09480、CA00g51890和CA03g02990分別在根、葉片、漿果和花苞優(yōu)勢表達。在GH18亞家族的CaCTL成員中，有9個成員（64.2%）呈現(xiàn)低表達態(tài)勢，優(yōu)勢表達基因主要集中在根和花組織中。綜上所述，有超過半數(shù)的CaCTL成員在正常生長的植株中顯示出較低的表達水平。

通過qRT-PCR分析，檢測CaCTL基因家族成員在受低溫（4 ℃）、高溫（42 ℃）以及干旱環(huán)境脅迫處理的幼苗中的表達特征。結果表明，與正常生長的幼苗相比，在受低溫脅迫的幼苗中，CA03g30190、CA03g30170、CA07g09480、CA07g02470、CA12g08860等5個成員相對表達量顯著提高，而CA12g20880、CA00g97030、CA03g31430、CA08g07860、CA11g13780和CA01g 20040等6個成員的表達明顯受到抑制（圖8）。其中CA07g02470的相對表達水平最多，被上調6倍;CA12g20880的相對表達水平最少，被下調3.3倍。在受42 ℃高溫脅迫的幼苗中，相對表達水平被上調成員數(shù)量少但幅度大，被下調的成員數(shù)量多但幅度小。其中CA12g20880、CA11g13790、CA03g 31390、CA03g30190和CA12g08860等5個成員的相對表達量顯著升高，與對照中的表達水平存在極顯著差異，并以CA12g20880的相對表達水平最高，被上調21.1倍，而CA00g97030、CA10g 09850、CA10g09860、CA03g31430、CA05g 11510、CA08g15980、CA01g20040、CA07g02470、CA04g 14160、CA11g13780、CA03g30170和CA07g02460共計12個成員的表達水平明顯下調。在受干旱脅迫的幼苗中，絕大多數(shù)的CaCTL成員呈現(xiàn)相對表達量被上調，CA06g08030、CA10g09860、CA03g 30180、CA03g31390、CA07g 09480、CA08g10220、CA01g15470、CA11g13790、CA07g02290、CA04g 08220、CA07g02470、CA12g 08860、CA07g02440、CA09g07340和CA07g02460共計15個成員響應強烈，尤其是CA07g02470的相對表達水平比對照幼苗植株高出114.7倍。僅有CA03g02990、CA05g11510、CA08g07870、CA08g07860、CA01g 20040和CA11g137806等6個成員顯示出明顯的下調。綜合來看，多數(shù)CaCTL成員的相對表達量在2個或者2個以上的逆境環(huán)境中都有上調。其中CA03g31390、CA11g13790和CA12g08860在3個脅迫處理中響應強烈，CA03g30190在低溫和高溫脅迫條件下均顯著上調，CA08g10220、CA07g 09480和CA07g02470在低溫以及干旱兩種環(huán)境中上調。CA03g31430、CA05g11510、CA08g07870、CA08g07860、CA01g20040和CA11g13780在3種處理中相對表達量都顯著低于對照中的表達水平。

CK代表正常環(huán)境條件，4 ℃、42 ℃和干旱為脅迫處理條件。橫坐標大于0為上調表達，小于0為下調表達?；疑煮w的基因為GH19亞家族成員，黑色字體的基因為GH18亞家族成員。

CK replicates the normal environmental conditions; 4 ℃， 42 ℃ and drought replicate the stress treatment conditions. The abscissa data greater than 0 replicate the upregulated expression level， less than 0 replicate the downregulated expression level. Genes with gray font color are members of the GH19 subfamily， while genes with black font color are members of the GH18 subfamily.

3 ?討論

幾丁質酶在植物的生長發(fā)育中扮演重要角色。在水稻中，幾丁質酶類基因OsCLP的T-DNA插入突變體（osclp）表現(xiàn)出根和莖生長的顯著阻滯[16]。在擬南芥中，幾丁質酶類基因AtCTL1和AtCTL2被敲除時，ctl1單突變體和ctl1ctl2雙突變體均表現(xiàn)出植株變矮，開花推遲和角果變小的表型[17]。在矮牽牛中，幾乎所有的花發(fā)育相關幾丁質酶活性都位于柱頭上，并且在花藥開裂后活性增加約5倍[18]，說明植物幾丁質酶在高等植物的有性生殖特別是授粉授精過程中也發(fā)揮功能。在本研究中，CaCTL成員在辣椒的根、莖、葉、花、漿果中的表達特征存在明顯差異，31個被鑒定的CaCTL成員可以區(qū)分為優(yōu)勢表達、微量表達和不表達3種模式。優(yōu)勢表達的成員在不同的組織中的表達水平存在差異，在同一器官的不同階段也呈現(xiàn)顯著變化。在花器官中，CA03g02990、CA12g20880、CA04g14160、CA00g51890、CA03g 31430、CA05g11510和CA01g20040在花苞到開放花的發(fā)育階段相對表達均呈下調模式，而CA07g13240、CA10g09860、CA07g09480、CA08g 10220、CA04g08220、CA07g02460和CA07g02440表現(xiàn)出上調表達，推測它們可能在花器官的發(fā)育中發(fā)揮作用。另外，CA03g02990、CA07g13240、CA00g51890、CA10g09860、CA07g09480、CA08g 10220和CA04g08220在營養(yǎng)器官（根、莖、葉）中優(yōu)勢表達，推測這些CaCTL成員或許在辣椒植株營養(yǎng)生長時期發(fā)揮作用。這些推測還需要進一步的試驗驗證。

植物作為一種固著生物，其生長發(fā)育過程深受立地環(huán)境的影響，嚴重時甚至導致死亡，比如異溫[19-21]、干旱[22]、水淹、鹽堿危害[23-24]等非生物脅迫，以及由灰霉菌[25]、甜菜夜蛾[26]等造成的生物脅迫。為了應對脅迫傷害，植物會合成一系列“病程相關蛋白”（pathogenesis-related proteins，PRs）度過危機。幾丁質酶也屬于一類PRs，在增強植物抗性脅迫方面的作用被日漸關注。本研究發(fā)現(xiàn)，CaCTL基因家族成員對非生物脅迫環(huán)境也產生了廣泛的響應。比如CA07g02470在低溫和干旱脅迫環(huán)境中的上調倍數(shù)顯著，CA12g20880對高溫脅迫響應強烈，可能在相應脅迫環(huán)境的應對中發(fā)揮重要作用。部分成員與多種脅迫響應相關，如CA03g31390、CA11g13790、CA12g08860、CA03g30190、CA08g10220、CA07g09480和CA07g 02470。這種現(xiàn)象在楊樹[27]和大豆[28]的CTL家族成員研究中也有出現(xiàn)，即同一個家族成員可以參與到不同的逆境響應機制中。本研究發(fā)現(xiàn)脅迫環(huán)境中響應較為敏感的CaCTL成員都具備有相應的順式元件，比如CA03g31390含有1個低溫響應元件、CA11g13790含有1個干旱響應元件、CA12g08860含有1個低溫響應元件以及3個防御應激元件、CA03g30190還有1個防御應激元件、CA08g10220含有1個干旱響應元件、CA07g09480含有1個干旱響應元件以及1個防御應激元件?；蛟S還有一些脅迫響應元件不在檢索的2 kb長度之內，但這些元件基本都與CaCTL成員在非脅迫環(huán)境中的響應規(guī)律相一致。已有研究報道表明，茉莉酸（JA）、水楊酸（SA）和脫落酸（ABA）等植物激素和植物信號物質都參與了植物體內各種應激信號通路的運行[29-32]。辣椒中含有CBD結構的幾丁質酶對外源激素會產生響應[13]，番茄CTL基因家族成員也已經被證實積極參與到植物激素和植物信號的反應過程[10]。本研究中順式元件預測發(fā)現(xiàn)，許多CaCTL成員的啟動子區(qū)域也包含有多個激素響應元件，說明這些成員可能在植物激素響應過程中發(fā)揮作用。但具體的觸發(fā)過程以及響應模式還有待進一步的研究探討。該研究結果豐富了CaCTL基因家族進化的研究，為后期CaCTL基因家族成員的功能鑒定和辣椒的抗性育種提供數(shù)據基礎。

參考文獻

[1] Ding X， Gopalakrishnan B， Johnson L B， et al. Insect resis-tance of transgenic tobacco expressing an insect chitinase gene[J]. Transgenic Research， 1998， 7（2）： 77-84.

[2] Li H， Greene L H. Sequence and structural analysis of the chitinase insertion domain reveals two conserved motifs in-volved in chitin-binding[J]. PLoS One， 2010， 5（1）： e8654.

[3] Singh H R， Deka M， Das S. Enhanced resistance to blister blight in transgenic tea （Camellia sinensis [L.] O. Kuntze） by overexpression of class I chitinase gene from potato （Solanum tuberosum）[J]. Functional & Integrative Genomics， 2015， 15（4）： 461-480.

[4] Robinson S P， Jacobs A K， Dry I B. A class IV chitinase is highly expressed in grape berries during ripening[J]. Plant Physiology， 1997， 114（3）： 771-778.

[5] Hong J K， Hwang B K. Induction by pathogen， salt and drought of a basic class II chitinase mRNA and its in situ localization in pepper （Capsicum annuum）[J]. Physiologia Plantarum， 2002， 114（4）： 549-558.

[6] Yeh S， Moffatt B A， Griffith M， et al. Chitinase genes res-ponsive to cold encode antifreeze proteins in winter ce-reals[J]. Plant Physiology， 2000， 124（3）： 1251-1264.

[7] Hermans C， Porco S， Verbruggen N， et al. Chitinase-like protein CTL1 plays a role in altering root system architecture in response to multiple environmental conditions[J]. Plant Physiology， 2010， 152（2）： 904-917.

[8] Wu J N， Wang Y M， Kim S T， et al. Characterization of a newly identified rice chitinase-like protein （OsCLP） homo-logous to xylanase inhibitor[J]. BMC Biotechnology， 2013， 13（1）： 1-8.

[9] Cletus J， Balasubramanian V， Vashisht D， et al. Transgenic expression of plant chitinases to enhance disease resis-tance[J]. Biotechnology Letters， 2013， 35（11）： 1719-1732.

[10] Cao J， Tan X N. comprehensive analysis of the chitinase family genes in tomato （Solanum lycopersicum）[J]. Plants， 2019， 8（3）： 52.

[11] Yamagami T， Funatsu G. Identification of the aspartic acid residue located at or near substrate-binding site of rye seed chitinase-C[J]. Bioscience， Biotechnology， and Biochemistry， 1998， 62（2）： 383-385.

[12] Taira T， Toma N， Ishihara M. Purification， characterization， and antifungal activity of chitinases from pineapple （Ananas comosus） leaf[J]. Bioscience， Biotechnology， and Biochemistry， 2005， 69（1）： 189-196.

[13] Ali M， Luo D X， Khan A， et al. Classification and genome-wide analysis of chitin-binding proteins gene family in pepper （Capsicum annuum L.） and transcriptional regulation to Phytophthora capsici， abiotic stresses and hormonal applications[J]. International Journal of Molecular Sciences， 2018， 19（8）： 2216.

[14] Passarinho P A， De Vries S C. Arabidopsis chitinases： a genomic survey[J]. The Arabidopsis Book， 2002， 1： e0023.

[15] Wan H J， Yuan W， Ruan M Y， et al. Identification of reference genes for reverse transcription quantitative real- time PCR normalization in pepper （Capsicum annuum L.）[J]. Biochemical and Biophysical Research Communications， 2011， 416（1-2）： 24-30.

[16] Wu J N， Wang Y M， Kim S G， et al. A secreted chitinase-like protein （OsCLP） supports root growth through calcium signaling in Oryza sativa[J]. Physiologia Plantarum， 2017， 161（2）： 273-284.

[17] Aktar Hossain M， Noh H N， Kim K I， et al. Mutation of the chitinase-like protein-encoding AtCTL2 gene enhances lignin accumulation in dark-grown Arabidopsis seedlings[J]. Journal of Plant Physiology， 2010， 167（8）： 650-658.

[18] Leung D W M. Involvement of plant chitinase in sexual reproduction of higher plants[J]. Phytochemistry， 1992， 31（6）： 1899-1900.

[19] Orvar B L， Sangwan V， Omann F， et al. Early steps in cold sensing by plant cells： the role of actin cytoskeleton and membrane fluidity[J]. Plant Journal， 2000， 23（6）： 785-794.

[20] Pearce R. Plant freezing and damage[J]. 2001， 87（4）： 417-424.

[21] Haak D C， Fukao T， Grene R， et al. Multilevel regulation of abiotic stress responses in plants[J]. Frontiers in Plant Science， 2017， 8： 1564.

[22] Jaleel C A， Gopi R， Sankar B， et al. Differential responses in water use efficiency in two varieties of Catharanthus roseus under drought stress[J]. Comptes Rendus Biologies， 2008， 331（1）： 42-47.

[23] Chen W C， Cui P J， Sun H Y， et al. Comparative effects of salt and alkali stresses on organic acid accumulation and io-nic balance of seabuckthorn （Hippophae rhamnoides L.）[J]. Industrial Crops and Products， 2009， 30（3）： 351-358.

[24] Ismail A， Takeda S， Nick P. Life and death under salt stress： same players， different timing[J]. Journal of Experimental Botany， 2014， 65（12）： 2963-2979.

[25] El Oirdi M， El Rahman T A， Rigano L， et al. Botrytis cinerea manipulates the antagonistic effects between immune pathways to promote disease development in tomato[J]. Plant Cell， 2011， 23（6）： 2405-2421.

[26] Musser R O， Hum-Musser S M， Eichenseer H， et al. Cater-pillar saliva beats plant defences[J]. Nature， 2002， 416（6881）： 599-600.

[27] Jiang C， Song J L， Huang R F， et al. Molecular cloning and expression analysis of a class V chitinase gene from popu-lus[J]. Molecular Plant Breeding， 2013， 11（2）： 225-231.

[28] Mészáros P， Rybansky ?， Spieb N， et al. Plant chitinase responses to different metal-type stresses reveal specificity[J]. Plant Cell Reports， 2014， 33（11）： 1789-1799.

[29] Fujita M， Fujita Y， Noutoshi Y， et al. Crosstalk between abiotic and biotic stress responses： a current view from the points of convergence in the stress signaling networks[J]. Current Opinion in Plant Biology， 2006， 9（4）： 436-442.

[30] Ma X W， Ma F W， Mi Y F， et al. Morphological and physi-ological responses of two contrasting Malus species to ex-ogenous abscisic acid application[J]. Plant Growth Regula-tion， 2008， 56（1）： 77-87.

[31] Liu Y J， Jiang H F， Zhao Z G， et al. Abscisic acid is involved in brassinosteroids-induced chilling tolerance in the suspension cultured cells from Chorispora bungeana[J]. Journal of Plant Physiology， 2011， 168（9）： 853-862.

[32] Guo M， Zhai Y F， Lu J P， et al. Characterization of CaHsp70-1， a pepper heat-shock protein gene in response to heat stress and some regulation exogenous substances in Capsicum annuum L.[J]. International Journal of Molecular Sciences， 2014， 15（11）： 19741-19759.

責任編輯：黃東杰