羅布麻CesA基因家族的生物信息學分析

解盛 李國旗 宋立肖 謝博勛 王雅芳 劉星

摘?要：?植物體中纖維素是細胞壁形成的主要成分，不僅參與細胞形態(tài)的建成，調控細胞發(fā)育，還參與細胞內多種細胞信號轉導途徑，進而影響植物體的生長發(fā)育。纖維素合酶是植物體合成纖維素的主要酶類。為了探究CesA基因家族對羅布麻生長發(fā)育及纖維素合成的調控機理，該文通過生物信息學分析方法，從基因家族鑒定、結構分析、蛋白理化性質與多級結構預測、亞細胞定位、信號肽、進化關系和順式作用元件等方面，對羅布麻CesA基因家族進行系統(tǒng)鑒定和分子特征分析。結果表明：基于全基因組測序，羅布麻CesA基因家族鑒定含有15個成員，分布在羅布麻11條染色體中的8條上，其編碼蛋白的氨基酸數(shù)量為730～1 158，相對分子質量81 280.81～130 123.18 kDa，理論等電點6.18～8.83。除了AvCesA3、AvCesA5、AvCesA7、AvCesA10和AvCesA11蛋白為穩(wěn)定蛋白，其余成員均為不穩(wěn)定蛋白;除了AvCesA12蛋白為疏水性蛋白，其余成員均為親水性蛋白。該家族成員包含3～14個外顯子，8～15個保守基序。編碼蛋白主要分布于質膜與高爾基體上，無信號肽，二級結構以無規(guī)則卷曲與α-螺旋為主要構成元件。AvCesA15蛋白的跨膜結構域和三級結構與其他成員存在顯著不同。羅布麻CesA基因進化時主要受純化選擇作用。對上游1 500 bp區(qū)域順式作用元件分析，結果顯示羅布麻CesA基因受到光、溫度、水分、氧氣等環(huán)境因子及生長素、赤霉素、脫落酸、乙烯、水楊酸等植物激素調控。該研究為進一步探究羅布麻CesA基因家族的生物學功能、提高纖維品質與品種改良奠定理論基礎。

關鍵詞： 羅布麻， 纖維發(fā)育， 生物信息學分析， CesA基因家族， 生長發(fā)育

中圖分類號：?Q943

文獻標識碼：?A

文章編號：?1000-3142（2021）04-0522-13

Abstract：?Cellulose is the main component of the cell wall in plants. It is not only participates in cell morphology and development， but also participates in various cellular signaling transduction pathways in cells， so as to affect the growth and development of plants. Cellulose synthase is the main enzyme used to?synthesize cellulose in plants. Here， we explore the regulatory mechanism of the CesA gene family for cellulose synthesis， growth and development in Apocynum venetum. Using bioinformatic methods， we conducted systematic identification and molecular characterization of A. venetum CesA gene family which include gene family member identification， structural analysis， protein physicochemical properties and multi-level structure prediction， subcellular localization， signal peptides， evolutionary relationships， cis-acting elements. The results were as follows： Based on whole-genome sequencing， the A. venetum CesA gene family contained 15 members， which were distributed on 8 of the 11 A. venetum chromosomes. These family proteins were encoded by 730-1 158 amino acids， with the molecular weight at 81 280.81-130 123.18 kDa and theoretical isoelectric point at 6.18-8.83. Among these proteins， AvCesA3， AvCesA5， AvCesA7， AvCesA10 and AvCesA11 proteins were stable proteins， the rest of the members were unstable proteins. Except for AvCesA12 protein as a hydrophobic protein， the remaining members were hydropathicity proteins. Members of this family contained 3-14 exons and 8-15 conserved motifs. The encoded proteins were mainly distributed on the plasma membrane and the Golgi apparatus， without a clear signal peptide. The secondary structure was mainly composed of random coils and α-helices. The transmembrane domained and tertiary structure of AvCesA15 protein were significantly different from other members. Evolutionary selection of A. venetum CesA gene was mainly affected by purification selection. The analysis of cis-acting elements in the upstream 1 500 bp region showed that these gene might be regulated by environmental factors such as light， temperature， water， oxygen and plant hormones such as auxin， gibberellin， abscisic acid， ethylene， salicylic acid， and so on. This study lays a theoretical foundation for further exploring the biological function of CesA gene family and improving fiber quality and variety of A. venetum.

Key words： Apocynum venetum， fiber development， bioinformatics analysis， CesA gene family， growth and development

羅布麻（Apocynum venetum）屬夾竹桃科（Apocynaceae）羅布麻屬（Apocynum）多年生宿根草本植物，具有耐旱、耐寒、耐堿、適應性強等特點（王東清等，2012;劉克彪和姜生秀，2018）。隨著研究的深入，羅布麻的多種經(jīng)濟價值被挖掘，主要體現(xiàn)在纖維的紡織利用（徐宗昌等，2018），黃酮類、槲皮素等化學成分的藥用價值（侯晉軍等，2006），茶飲保健功能（錢學射等，2005）等。自20世紀50年代，對擁有“野生纖維之王”美譽的羅布麻的研究從未間斷，但其分子層次的研究還比較薄弱。纖維素合酶基因（CesA）編碼纖維素合酶亞基，是廣泛參與陸地植物纖維素生物合成的多基因家族（Yin et al.， 2014）。因此，對纖維素合酶的研究是一個有價值的研究課題。通過在羅布麻基因組中系統(tǒng)鑒定CesA基因家族成員，并對其基因結構、蛋白理化性質、結構預測、進化關系進行分析，對進一步認識羅布麻纖維發(fā)育過程及培育改良羅布麻品種具有重要意義。

纖維素由24～36個經(jīng)β-1，4糖苷鍵連接的葡聚糖殘基（500～14 000個單體）鏈組成，這些殘基通過氫鍵相互連接，形成規(guī)則的晶體結構，變成穩(wěn)定體。纖維素是植物細胞壁的主要組成成分，也是地球上最豐富的生物聚合體和可再生碳源（Andrew & Specht， 2011）。Cantarel et al.（2009）研究顯示，2-型糖基轉移酶（GT2）廣泛分布在植物界，而纖維素合成酶屬于2-型糖基轉移酶合成酶家族中的一個超基因家族。纖維素合成酶CesA基因編碼的蛋白能夠在質膜上以復合體的形式合成纖維素（ Taylor et al.， 2000;Robyn， 2001）。纖維素做為植物細胞壁的重要組成之一，不僅參與細胞形態(tài)的建成，控制細胞生長，而且還能參與到多種生物、非生物脅迫的信號應答調控網(wǎng)絡（Burton et al.， 2010）。因此，CesA基因在植物的生長發(fā)育過程中起著非常重要的作用。第一個在植物中報道的CesA基因是從棉花（陸地棉）纖維中根據(jù)其與細菌CesA基因序列的高相似性鑒定出來的（Pear et al.， 1996），后續(xù)學者在模式生物擬南芥中鑒定出了10個CesA基因，命名為CesA1、CesA2……CesA10（Farrokhi et al.， 2006），其中CesA1、CesA3和CesA6-like（CesA2，CesA5，CesA6或CesA9）成員之一的蛋白質產(chǎn)物形成具有18～24個CesA異構體的異三聚體構型的CesA復合物（CSCs），并合成初生細胞壁纖維素的成分（Desprez et al.， 2007;Persson et al.， 2007;Mcfarlane et al.， 2014）。通過免疫沉淀和突變體分析，CesA4、CesA7和CesA8在次生細胞壁形成過程中直接相互作用，是合成纖維素所必須的（Somerville， 1997;Taylor et al.， 2003）。此外，在植物生長過程中，每個重要基因的突變，造成的結果都是致命性的（例如CesA1、CesA3基因）（Persson et al.， 2007）。有研究結果顯示，每個CesA6-like基因之間存在部分功能冗余，在此基礎上有學者提議將CesA6-like基因列為一個單獨的亞家族，雖然一個CesA基因可以替代另一個，但又不能完全取代，所以推斷每一個CesA6-like基因可能具有專門的功能（Ruprecht et al.， 2017; Hu et al.， 2018）。除模式生物外，CesA基因家族也在其他植物中被鑒定，如大麥（Hordeum vulgare）（Burton et al.， 2004）、火炬松（Pinus taeda）（Naim & Haselkorn， 2005）、玉米（Zea mays）（Appenzeller et al.， 2004;張曉榕等，2019）、水稻（Oryza sativa）（Yin et al.， 2009）等。CesA基因家族成員在不同植物中具有相似功能，但由于環(huán)境和物種的差異，不同植物在進化過程中的不同選擇，CesA基因家族成員在不同物種中均有不同程度的分化。本實驗室通過對羅布麻進行全基因組測序（結果未發(fā)表），為其進化和研究相關基因的功能關系奠定理論基礎，而CesA基因家族在羅布麻基因組中并未進行深入分析和比較分析。本研究利用生物信息學技術鑒定羅布麻CesA基因家族，并系統(tǒng)分析其在基因組中的基因結構、分布、進化特征，為闡述羅布麻CesA基因家族分化歷程及生物學功能奠定基礎。

1?材料與方法

材料為采自寧夏石嘴山市平羅縣羅布麻實驗基地的羅布麻。實驗于2018—2019年在寧夏大學西北土地退化與生態(tài)恢復國家重點實驗室培育基地實驗室進行。

1.1 羅布麻CesA基因家族成員的鑒定

本研究主要基于羅布麻全基因組測序進行系統(tǒng)分析，同時利用引自美國華盛頓州艾倫斯堡（Washington Ellensburg）的大麻狀羅布麻基因組測序結果進行比較分析。研究中使用的羅布麻基因組數(shù)據(jù)來自實驗室前期測序所得。利用現(xiàn)有羅布麻注釋基因的蛋白序列，與從擬南芥數(shù)據(jù)庫（TAIR）（https：//www.arabidopsis.org/）中檢索得到的10個CesA基因家族成員蛋白序列進行BLASTP比對，E-value<1e-10，Identity>40%，輸出最優(yōu)比對結果，篩選羅布麻候選CesA基因。利用Pfam（http：//pfam.sanger.ac.uk/search）及 SMART（http：//smart.emblheidelberg. de/）進行驗證分析，刪除缺失結構域的序列，最終獲得羅布麻CesA基因家族候選基因。

1.2 羅布麻CesA基因家族成員染色體定位

由羅布麻全基因組注釋文件中，獲取CesA基因在染色體上的位置信息，利用在線工具MG2C（http：//mg2c.iask.in/mg2c_v2.0/）繪制染色體定位圖。

1.3 羅布麻CesA基因家族成員分析

利用 ExPASy中的ProtParam（https：//web.expasy.org/protparam/）工具對羅布麻CesA蛋白基本理化性質進行分析，包含氨基酸長度、理論等電點、相對分子質量、親水性總平均值、不穩(wěn)定系數(shù)和脂溶指數(shù)等指標。

1.4 羅布麻CesA基因家族成員結構和保守結構域分析

CesA基因的外顯子和內含子位置信息參考羅布麻基因組注釋信息GFF3文件，利用在線工具GSDS（Hu et al.， 2015）（http：//gsds.cbi.pku.edu.cn/）對羅布麻CesA基因家族成員進行結構分析。應用在線工具MEME（http：//meme-suite.org/tools/meme）預測分析CesA蛋白序列的結構域，搜尋 motif 值設置為15，其他設定為默認參數(shù)。通過 TBtools軟件繪制出 MEME 結構。

1.5 羅布麻CesA基因家族成員蛋白結構與亞細胞定位分析

利用在線平臺SOPMA（https：//npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl？page=/NPSA/ npsa_seccons.html）分析AvCesA蛋白二級結構。同時利用SignalP-5.0 Server（http：//www.cbs.dtu.dk/services/SignalP/）對羅布麻CesA蛋白進行信號肽分析。應用在線工具ProtComp Version 9.0（http：//linux1.softberry.com/berry.phtml？group=programs&subgroup= proloc&topic=protcomppl）對羅布麻CesA蛋白進行亞細胞定位分析。

1.6 羅布麻CesA基因家族成員蛋白跨膜螺旋與三級結構預測

利用TMHMM Server v. 2.0（http：//www.cbs.dtu.dk/services/TMHMM/）工具對羅布麻CesA蛋白跨膜螺旋進行預測。利用Phyre2（Kelley et al.， 2015）（http：//www.sbg.bio.ic.ac.uk/ phyre2/ html/page.cgi？id=index）對羅布麻CesA蛋白進行三級結構預測，并通過Jmol進行展示。

1.7 羅布麻CesA基因家族成員共線性與進化選擇壓力分析

利用MCScanX軟件對羅布麻CesA基因家族與雙子葉模式生物擬南芥和單子葉模式生物水稻進行共線性分析。通過BLAST建庫和KaKs_ Calculator 2.0（Wang et al.， 2010）對羅布麻CesA基因核苷酸的非同義替換率（Ka）與同義替換率（Ks）進行計算，獲取基因的Ka/Ks比率，進行選擇壓力分析。

1.8 羅布麻CesA基因家族成員系統(tǒng)進化分析

將檢索得到的擬南芥（Arabidopsis thaliana）、亞麻（Linum usitatissimum）、水稻（Oryza sativa）與羅布麻（Apocynum venetum）的CesA蛋白的氨基酸序列通過Clustal X進行蛋白序列比對，并通過鄰接法（Neighbour-Joining，NJ）在MEGA-X（Kumar et al.， 2018）中構建CesA基因家族系統(tǒng)進化樹，自展值（Bootstrap）設定為1 000，其他參數(shù)為系統(tǒng)默認值。使用EvolView（https：// evolgenius.info//evolview-v2/#login）工具展示系統(tǒng)進化樹。

1.9 羅布麻CesA基因家族成員順式作用元件分析

利用TBtools軟件（Chen et al.， 2018）獲得CesA基因上游1 500 bp序列，利用在線網(wǎng)站PlantCARE（http：//bioinformatics. psb.ugent.be/webtools/plantcare/html/）對獲得的基因上游序列進行預測，利用在線工具GSDS進行繪圖展示，并對主要的順勢作用元件進行探討。

2?結果與分析

2.1 羅布麻CesA基因家族成員的鑒定

從羅布麻全基因組中共篩選到35個擬南芥CesA的同源蛋白序列，將篩選結果輸入到Pfam 數(shù)據(jù)庫與SMART數(shù)據(jù)庫進行比對，刪除Cellulose-synt保守結構域缺失的序列，最終獲得15條羅布麻CesA序列，依次命名為AvCesA1、AvCesA2……AvCesA15。

2.2 羅布麻CesA基因家族成員染色體定位與基因復制分析

根據(jù)CesA基因染色體定位分析，15條CesA基因不均勻地分布在8條染色體上（圖1）。其中，11號染色體上分布的基因最多，達到4條，1、9和10號染色體上只分布1條基因，其余染色體上各含有2條基因?；蜷g無串聯(lián)重復現(xiàn)象，推斷羅布麻CesA基因家族成員間存在功能分化。

2.3 羅布麻CesA基因家族成員蛋白理化性質分析

對羅布麻CesA家族蛋白的理化性質進行分析，結果表明羅布麻CesA蛋白長度在730～1 158 aa之間，相對分子質量在81 280.81～130 123.18 kDa之間，理論等電點在6.18～8.83之間，編碼氨基酸序列最長的是AvCesA4，其相對分子質量最大，為130 123.18 kDa;編碼氨基酸序列最小的是AvCesA15，其相對分子質量也最小，為81 280.81 kDa。當?shù)鞍撞环€(wěn)定系數(shù)>40判斷其為不穩(wěn)定蛋白，AvCesA3、AvCesA5、AvCesA7、AvCesA10和AvCesA11蛋白不穩(wěn)定系數(shù)均小于40，為穩(wěn)定蛋白，AvCesA3蛋白的不穩(wěn)定系數(shù)最小，為37.84;其余家族成員蛋白不穩(wěn)定系數(shù)均大于40，為不穩(wěn)定蛋白，AvCesA15蛋白的不穩(wěn)定系數(shù)最大，為50.6。羅布麻CesA蛋白脂肪族氨基酸指數(shù)為66.96～93.32，親水性總平均值在-0.556～0.016之間。除AvCesA12蛋白親水性總平均值為0.016（大于0），為疏水性蛋白;其余成員蛋白親水性總平均值均小于0，為親水性蛋白。

2.4 羅布麻CesA基因家族成員基因結構和保守結構域分析

利用在線工具GSDS對15條羅布麻CesA基因家族成員進行基因結構分析，結果見圖2。從圖2可以看出，CesA基因家族外顯子數(shù)量主要分布在3～14之間。利用MEME對15條羅布麻CesA蛋白序列進行分析，結果見圖3。由圖3可知， 共鑒定獲得15個保守基序（motif）， 除AvCesA12與AvCesA15外，其他家族成員均含有15個motif，且分布相對比較均勻，其分布的數(shù)量與位置也基本相同，說明具有高度的保守性，由此推斷這部分羅布麻CesA基因家族成員的結構域和功能單位基本相似。AvCesA15所含基序最少，僅含有8個，AvCesA12所含基序為12個，說明AvCesA12和AvCesA15與其他羅布麻CesA家族成員可能存在功能上的差異。

2.5 羅布麻CesA基因家族成員蛋白結構特點、亞細胞定位與信號肽分析

對15條羅布麻CesA蛋白的二級結構的預測結果見表2。由表2可知，CesA蛋白均由α-螺旋、延伸鏈、β-轉角、無規(guī)則卷曲四部分組成。其中CesA蛋白的無規(guī)則卷曲比例最高，α-螺旋次之，而β-轉角比例最小。根據(jù)亞細胞定位分析，結果顯示羅布麻CesA基因家族主要定位在質膜與高爾基體上，這可能與基因的生物學功能相關聯(lián)。信號肽分析結果顯示，羅布麻CesA基因家族成員均不含有信號肽。

2.6 羅布麻CesA基因家族成員蛋白跨膜結構域與三級結構預測

跨膜結構一般由20個左右的疏水氨基酸殘基形成α-螺旋，是蛋白質與膜內蛋白的靜電相互作用和氫鍵鍵合作用與膜結合的一段氨基酸片段?？缒そY構域是膜中蛋白與生物膜脂質雙分子相結合的主要部位， 固著于細胞膜上起“錨定”作用。通過對跨膜結構域的預測和分析，有助于了解蛋白質的結構、功能以及在細胞中的作用位點等。通過前人研究結果可知，當前的基因組數(shù)據(jù)中，有20%～30%的基因產(chǎn)物被預測為膜蛋白（Taylor et al.， 2000），具備多種生物學功能。利用TMHMM對羅布麻CesA基因家族成員蛋白的跨膜結構進行預測，結果見圖4。從圖4可以看出，除AvCesA12與AvCesA15蛋白外，其余成員蛋白在C端均含有6個完整的跨膜結構域，且每個組成跨膜結構域的氨基酸數(shù)量與位置均相似，推斷其可能有相同的結構和功能。AvCesA12與AvCesA15可能因為發(fā)生突變，相應功能發(fā)生變化，從而導致跨膜結構缺失。除AvCesA2、AvCesA3、AvCesA4、AvCesA10、AvCesA13與AvCesA14蛋白外，其余成員蛋白在近N端均含有2個完整的跨膜結構域。

通過Phyre2對羅布麻CesA蛋白進行三級結構預測，結果見圖5。由圖5可知，除AvCesA15蛋白外，其余成員蛋白三級結構在空間結構上相似度很高，推斷其可能具有相似的生物學功能，而AvCesA15結構差異較大，可能發(fā)生變異，導致其結構與功能發(fā)生改變。

2.7 羅布麻CesA基因家族成員共線性與進化選擇壓力分析

利用MCScanX將羅布麻分別與雙子葉模式生物擬南芥和單子葉模式生物水稻建立2組CesA基因家族的共線性比較圖譜（圖 6）。共線性關系（同源基因對）在擬南芥CesA中有11對，水稻CesA中有4對。其中AvCesA4與AvCesA8在兩個圖譜中都存在共線性關系。相較于雙子葉植物，羅布麻CesA基因家族與單子葉植物存在共線性關系的基因較少，其原因可能是植物起源基因組不同而造成的。根據(jù)Abrouk et al.（2010）的研究，雙子葉植物基因組可能源于7條古染色體，單子葉植物基因組可能源于5條古染色體。由此，我們推斷AvCesA4與AvCesA8在進化過程中高度保守，其功能有待發(fā)掘。

非同義替換率（Ka）與同義替換率（Ks）的比值可以判斷是否有選擇壓力作用于這個蛋白質編碼基因，從而反映生物進化的過程中所受到的選擇壓力。對來源于羅布麻CesA基因家族成員間親緣關系最近的8條基因序列進行Ka/Ks分析，結果見表3。由表3可知，基因對之間的Ka/Ks的比值均小于1，顯示羅布麻CesA基因在進化時主要受純化選擇作用。

2.8 羅布麻CesA基因家族成員系統(tǒng)進化分析

為進一步了解羅布麻CesA基因家族進化關系，通過MEGA-X對羅布麻（15條）、擬南芥（10條）、亞麻（14條）、水稻（11條）共50條CesA蛋白序列構建NJ系統(tǒng)進化樹（圖7）。由圖7可知，聚類樹分為5個亞族，羅布麻CesA基因在每一亞族中均有分布，表明各物種間CesA基因家族具有明顯的同源關系。

2.9 羅布麻CesA基因家族成員順式作用元件分析

通過對羅布麻CesA基因上游1 500 bp序列進行分析，預測基因的順式作用元件，結果見圖8。從圖8可以看出，將CesA基因家族順式作用元件劃分為基礎元件、植物激素調控元件、逆境脅迫相關元件、光反應元件、生長發(fā)育調控元件等5大類?；A元件包括TATA-box和CAAT-box等，數(shù)量最多。光反應元件包括Sp1、MRE、GT1-motif、G-Box、G-box、3-AF1 binding site、ACE、ATC-motif、ACA-motif、TCT-motif等，數(shù)量排第2位。植物激素調控元件包含生長素響應元件TGA-element、AuxRR-core，赤霉素響應元件P-box、F-Box、GARE-motif、CARE，乙烯應答元件ERE，脫落酸響應元件ABRE、ABRE4，水楊酸響應元件TCA-element，參與到茉莉酸甲酯信號轉導途徑的TGACG-motif、CGTCA-motif，茉莉酮酸酯和激發(fā)子應答元件JERE。逆境脅迫相關元件包含損傷誘導元件WRE3、WUN-motif，損傷與病原響應元件W-box，防御與應激反應元件TC-rich repeats，低溫調控元件TCA、LTR，脅迫響應元件STRE，干旱響應元件MBS，缺氧特異性誘導GC-motif，干旱和滲透脅迫誘導元件DRE1，厭氧誘導調控元件ARE。生長發(fā)育調控元件包含胚乳表達所需元件GCN4_motif，玉米醇溶蛋白代謝調控中的調控元件O2-site，參與種子特異性調控的順式作用調控元件RY-element，葉肉細胞分化的相關元件HD-Zip1，參與細胞周期調控的順式作用元件MSA-like，與分生組織特異性激活相關的順式調控元件CCGTCC-box，與分生組織表達相關的順式作用調控元件CCGTCC motif 、CAT-box，參與晝夜節(jié)律調控的順式調控元件circadian，根特異表達順式作用調控元件as-1。羅布麻CesA基因可能受到光、溫度、水分、氧氣等環(huán)境因子及內源激素調控，并且這些順式作用元件有可能參與到羅布麻抗逆機制中，在植物生長發(fā)育過程中起到防御與保護作用。除此之外，該段序列富含參與植株生長發(fā)育的順式作用元件，在植株的生長與纖維素合成中有重大作用。

3?討論與結論

對羅布麻、亞麻、擬南芥、水稻CesA基因家族構建進化樹，由聚類分析可知， AvCesA2與AvCesA3與擬南芥中參與初生細胞壁纖維素合成基因AtCesA1、AtCesA3（Desprez et al.， 2007）聚類在同一分支，推斷這些羅布麻CesA基因可能參與初生細胞壁纖維素的合成。AvCesA5、AvCesA8、AvCesA10分別與擬南芥中參與次生細胞壁纖維素合成基因AtCesA7、AtCesA4、AtCesA8（Arioli et al.， 1998）的親緣關系較近，因此推斷這些AvCesA基因也可能參與次生細胞壁纖維素的合成。AvCesA1、AvCesA4、AvCesA6、AvCesA9、AvCesA11、AvCesA12、AvCesA15聚類在同一分支，同緣關系顯著，可知這些基因可能具有類似的生物功能。

擬南芥CesA基因家族蛋白在N端都含有一個保守的鋅指蛋白結構域，該保守域能在CesA蛋白形成蛋白復合體時進行蛋白之間的相互識別（周曉馥等，2002）。通過對羅布麻CesA基因家族蛋白保守基序分析顯示，多數(shù)成員包含的保守結構域基本一致，具有高度的保守性，這說明蛋白功能表達方面存在相似性。AvCesA12 與AvCesA15相比較，缺失部分保守結構域，可能是基因序列發(fā)生突變引起，從而使蛋白執(zhí)行功能發(fā)生變化，這一現(xiàn)象在跨膜結構預測和蛋白三級結構中均有體現(xiàn)。AvCesA12蛋白的跨膜結構預測中C端第六個跨膜結構缺失，AvCesA15蛋白的跨膜結構預測中C端無跨膜結構，并且三級構象發(fā)生明顯變化。在羅布麻CesA蛋白跨膜結構預測中發(fā)現(xiàn)最多包含8個跨膜結構，與Nawaz et al.（2019）的總結一致，在植物CesA蛋白中最多存在8個跨膜結構域。徐昌宗等（2017）的研究結果發(fā)現(xiàn)，普通煙草CesA蛋白跨膜結構數(shù)量差異與該基因參與初生細胞壁或者次生細胞壁纖維素合成間存在一定關系。羅布麻CesA蛋白跨膜結構數(shù)量存在差異，以此為思路有助于探究羅布麻CesA基因的功能。

順式作用元件能被特定的轉錄因子識別，對下游基因的特異性表達起調控作用。對羅布麻CesA基因上游1 500 bp區(qū)域進行順式作用元件分析，該區(qū)域包含多種順式作用元件，在植物生長發(fā)育過程中可能存在重要作用。水楊酸是參與植物光敏反應和系統(tǒng)獲得抗性反應的重要信號分子（陳沖等， 2018），茉莉酸甲酯在植物受到病原侵害、創(chuàng)傷或者脅迫時誘導防御基因表達的信號分子，且在植株生命周期的各個階段也發(fā)揮著調節(jié)作用（Titarenko et al.， 1997）。除此之外，多種激素與環(huán)境因素可以影響CesA基因的表達。魏凱莉等（2017）的研究顯示，赤霉素處理下楊樹次生壁CesA基因表達量明顯上調，而生長素與乙烯處理下楊樹次生壁CesA基因表達量下調，影響細胞壁中纖維素含量。油菜素甾醇和光照可以促進擬南芥下胚軸伸長和初生壁CesA基因轉錄，促使擬南芥株形增高（Xie et al.， 2011;Leivar & Quail， 2011）。在水稻中，赤霉素信號轉導可以促進纖維素的合成，影響水稻植株節(jié)間的發(fā)育（Huang et al.， 2015）。擬南芥AtCesA8突變體植株，抗旱性顯著增強（Chen et al.， 2005）。龍眼DiCesA基因被脫落酸處理表達量下調，可能參與脫落酸對干旱脅迫的調控，通過對MeJA處理下DiCesA基因表達量下調，證實DiCesA參與龍眼對干旱脅迫的調控（朱永靜等，2020）。玉米CesA基因在種子生長階段參與抵御干旱脅迫的調控（張曉榕等，2019）。在羅布麻AvCesA基因家族成員的啟動子區(qū)域都包含多種激素響應元件，并且80%以上包含防御與逆境響應元件，60%以上包含了干旱脅迫以及光響應MYB結合位點。由此可見，CesA基因家族包含的大量光響應、脅迫響應、激素響應等元件，在植株生長發(fā)育與抗逆防御方面發(fā)揮著巨大作用。

在羅布麻纖維發(fā)育有關的CesA類基因家族中，被報道的基因數(shù)目很少，并且大多數(shù)基因的功能還不清楚。本研究通過生物信息學方法分析了羅布麻基因組中的CesA基因家族的鑒定、基因結構、蛋白特性及順式作用元件分析等，對于發(fā)掘新的與植物纖維發(fā)育相關的基因，加快羅布麻纖維發(fā)育分子機制的研究具有十分重要的意義，也為羅布麻培育高質量纖維育種研究奠定了基礎。

參考文獻：

ABROUK M， MURAT F， PONT C， et al.，2010. Palaeogenomics of plants： synteny-based modelling of extinct ancestors[J]. Trend Plant Sci， 15（9）： 479-487.

ANDREW C， SPECHT CD， 2011. Understanding plant cellulose synthases through a comprehensive investigation of the cellulose synthase family sequences[J]. Front Plant Sci， 2（5）： 2-5.

APPENZELLER L， DOBLIN M， BARREIRO R， et al.，2004. Cellulose synthesis in maize： Isolation and expression analysis of the cellulose synthase （CesA） gene family[J]. Cellulose， 11（3-4）： 287-299.

ARIOLI T， PENG LC， BETZNER AS， et al.，1998. Molecular analysis of cellulose biosynthesis in Arabidopsis[J]. Science， 279（5351）： 717-720.

BURTON RA， 2004. The CesA gene family of Barley. Quantitative analysis of transcripts reveals two groups of co-expressed genes[J]. Plant Physiol， 134（1）： 224-236.

BURTON RA， GIDLEY MJ， FINCHER GB， 2010. Heterogeneity in the chemistry， structure and function of plant cell walls[J]. Nat Chem Biol， 6（10）： 724-732.

CANTAREL BL， COUTINHO PM， RANCUREL C， et al.，2009. The carbohydrate-active enzymes database （CAZy）： an expert resource for glycogenomics[J]. Nucl Acid Res， 37（Suppl. 1）： D233-D238.

CHEN C， LIU S， WANG DD， et al.，2018. Identification of PCD induced by salicylic acid and expression analysis of genes related with the PCD in cucumber[J]. Acta Agric Boreal-Sin， 33（6）： 60-67.[陳沖，劉雙，王丹丹，等，2018. 水楊酸誘導黃瓜PCD的鑒定及相關基因的表達分析[J]. 華北農學報，33（6）：60-67.]

CHEN CJ， CHEN H， ZHANG Y， et al.，2020. TBtools： An integrative toolkit developed for interactive analyses of big biological data[J]. Mol Plant， 13（8）： 1194-1202.

CHEN ZZ， HONG XH， ZHANG HR， et al.， 2005. Disruption of the cellulose synthase gene， AtCesA8/IRX1， enhances drought and osmotic stress tolerance in Arabidopsis[J]. Plant J， 43（2）： 273-283.

DESPREZ T， JURANIEC M， CROWELL EF， et al.，2007. Organization of cellulose synthase complexes involved in primary cell wall synthesis in Arabidopsis thaliana[J]. Proc Natl Acad Sci， 104（39）： 15572-15577.

FARROKHI N，BURTON RA， BROWNFIELD L， et al.，2006. Plant cell wall biosynthesis： Genetic， biochemical and functional genomics approaches to the identification of key genes[J]. Plant Biotechnol J， 4（2）： 145-167.

HOU JJ， HAN LW， YANG GE， et al.，2006. Advances in studies on chemical constituents and pharmacological activities of Apocatharum venetum[J]. Chin Trad Herb Drugs， 37（10）： 1603-1605.[侯晉軍， 韓利文， 楊官娥， 等， 2006. 羅布麻葉化學成分和藥理活性研究進展[J]. 中草藥， 37（10）：1603-1605.]

HU B， JIN JP， GUO AY， et al.， 2015. GSDS 2.0： An upgraded gene feature visualization server. Bioinformatics， 31（8）： 1296-1297.

HU HZ， ZHANG R， FENG SQ， et al.， 2018. Three AtCesA6-like members enhance biomass production by distinctively promoting cell growth in Arabidopsis[J]. Plant Biotechnol J， 16（5）：976-988.

HUANG D， WANG S， ZHANG B， et al.， 2015. A gibberellin-mediated DELLA-NAC signaling cascade regulates cellulose synthesis in rice[J]. Plant Cell， 27（6）： 1681-1696.

KELLEY LA， MEZULIS S， YATES CM， et al.， 2015. The Phyre2 web portal for protein modeling， prediction and analysis[J]. Nat Protoc， 10（6）： 845-858.

KROGH A， LARSSON B， VON HG， et al.， 2001. Predicting transmembrane protein topology with a hidden Markov model： Application to complete genomes[J]. J Mol Biol， 305（3）：?567-580.

KUMAR S， STECHER G， LI M， et al.， 2018. MEGA X： Molecular evolutionary genetics analysis across computing platforms[J]. Mol Biol Evol， 35（6）： 1547-1549.

LEIVAR P， QUAIL PH， 2011. PIFs： Pivotal components in a cellular signaling hub[J]. Trends Plant Sci， 16（1）： 19-28.

LIU KB， JIANG SX， 2016. Responses of Apocynum venetum seed germination to drought and salt stress[J]. Acta Pratac Sin， 25（5）： 214-221.[劉克彪， 姜生秀， 2016. 干旱和鈉鹽脅迫對羅布麻種子萌發(fā)的影響[J]. 草業(yè)學報， 25（5）：214-221.]

MCFARLANE HE， DRING A， PERSSON S， 2014. The cell biology of cellulose synthesis[J]. Ann Rev Plant Biol， 65（1）： 69-94.

NAIRN CJ， HASELKORN T， 2005. Three loblolly pine CesA genes expressed in developing xylem are orthologous to secondary cell wall CesA genes of angiosperms[J]. New Phytol， 166（3）： 907-915.

NAWAZ MA， LIN XA， CHAN TF， et al.， 2019. Characterization of cellulose synthase A （CESA） gene family in eudicots[J]. Biochem Genet， 57（2）： 248-272.

PEAR JR， KAWAGOE Y， SCHRECKENGOST WE， et al.， 1996. Higher plants contain homologs of the bacterial cela genes encoding the catalytic subunit of cellulose synthase[J]. Proc Natl Acad Sci， 93（22）： 12637-12642.

PERSSON S， PAREDEZ A， CARROLL A， et al.， 2007. Genetic evidence for three unique components in primary cell-wall cellulose synthase complexes in Arabidopsis[J]. Proc Natl Acad Sci， 104（39）： 15566-15571.

QIAN XS， ZHANG WM， GU GP， et al.， 2005. Health tea development of Apocynum venetum L. and its medicinal diet[J]. Chin Wild Plant Resour， 24（6）： 21-25.[錢學射， 張衛(wèi)明， 顧龔平， 等， 2005. 羅布麻保健茶的開發(fā)與藥膳[J]. 中國野生植物資源， 24（6）：21-25.]

ROBYN M， PERRIN， 2001. Cellulose： how many cellulose synthases to make a plant？[J]. Curr Biol， 11（6）： R213-R214.

RUPRECHT C， PROOST S， HERNANDEZ-CORONADO M， et al.， 2017. Phylogenomic analysis of gene co-expression networks reveals the evolution of functional modules[J]. Plant J， 90（3）： 447-465.

SOMERVILLE TCR， 1997. Collapsed xylem phenotype of Arabidopsis identifies mutants deficient in cellulose deposition in the secondary cell wall[J]. Plant Cell， 9（5）： 689-701.

TAYLOR NG， HOWELLS RM， HUTTLY AK， et al.， 2003. Interactions among three distinct CesA proteins essential for cellulose synthesis[J]. Proc Natl Acad Sci USA， 100： 1450-1455.

TAYLOR NG， LAURIE S， TURNER SR， 2000. Multiple cellulose synthase catalytic subunits are required for cellulose synthesis in Arabidopsis[J]. Plant cell， 12（12）： 2529-2540.

TITARENKO E， ROJO E， LEN J， et al.， 1997. Jasmonic acid-dependent and independent signaling pathways control wound-induced gene activation in Arabidopsis thaliana [J]. Plant Physiol， 115（2）： 817-826.

WANG DP， ZHANG YB， ZHANG Z， et al.， 2010. KaKs_Calculator 2.0： A toolkit incorporating gamma-series methods and sliding window strategies[J]. Genom Proteom Bioinform， 8（1）： 77-80.

WANG DQ， LI GQ， SU DX， 2012. Effect of drought stress on osmotic adjustment substances and activity of protective enzymes in two species of Apocynum[J]. J Arid Land Resour Environ， 26（12）： 177-181.[王東清， 李國旗， 蘇德喜， 2012. 干旱脅迫對兩種羅布麻滲透調節(jié)物質積累和保護酶活性的影響[J]. 干旱區(qū)資源與環(huán)境， 26（12）： 177-181.]

WEI KL， ZHOU HJ， JIANG C，et al.， 2017. Interaction and expression of secondary wall CESAs in Populus [J]. For Res， 30（2）： 245-253.[魏凱莉， 周厚君， 江成， 等， 2017. 楊樹次生壁纖維素合酶的表達與互作模式分析[J]. 林業(yè)科學研究， 30（2）： 245-253.]

XIE L， YANG C， WANG X， 2011. Brassinosteroids can regulate cellulose biosynthesis by controlling the expression of CESA genes in Arabidopsis[J]. J Exp Bot， 62（13）： 4495-4506.

XU ZC， KONG YZ， 2017. Genome-wide identification， subcellular localization and gene expression analysis of the members of CESA gene family in common tobacco[J]. Hereditas， 39（6）： 512-524.[徐宗昌， 孔英珍，

2017. 普通煙草CESA基因家族成員的鑒定、亞細胞定位及表達分析[J]. 遺傳， 39（6）： 512-524.]

XU ZC， ZHOU JH， ZHANG CS， et al.， 2018. Review of current research and utilization status of Apocynum venetum germplasm in China[J]. Bull Bot， 53（3）： 382-390.[徐宗昌， 周金輝， 張成省， 等， 2018. 我國羅布麻種質資源研究利用現(xiàn)狀[J]. 植物學報， 53（3）： 382-390.]

YIN Y， HUANG J， XU Y， 2009. The cellulose synthase superfamily in fully sequenced plants and algae[J]. Bmc Plant Biol， 9（1）： 99-101

YINY， JOHNS MA， CAO H， et al.， 2014. A survey of plant and algal genomes and transcriptomes reveals new insights into the evolution and function of the cellulose synthase superfamily[J]. BMC Genom， 15（1）： 260.

ZHOU XF， WANG JY， WANG XZ， 2002. Research progress of cellulose synthase genes in higher plant[J]. Hereditas， 24（3）： 376-378.[周曉馥， 王景余， 王興智， 2002. 植物纖維素合成酶基因的研究進展[J]. 遺傳， 24（3）： 376-378.]

ZHANG XR， TAN JF， WEN MQ， et al.， 2019. Systematic identification and functional study of CesA family in maize[J]. J NW A & F Univ （Nat Sci Ed）， 47（2）： 45-53.[張曉榕， 譚俊峰， 溫曼晴， 等， 2019. 玉米CesA家族的系統(tǒng)鑒定及功能研究[J]. 西北農林科技大學學報（自然科學版）， 47（2）： 45-53.]

ZHU YJ， LU BS， ZHANG ST， et al.，?2020. Identification and expression patterns of the cellulose synthase gene family in Dimocarpus longan ?LOUR[J]. Chin J Appl Environ Biol， 26（5）： 1235-1243.[朱永靜， 路保順， 張舒婷， 等， 2020. 龍眼纖維素合成酶基因家族成員鑒定及表達模式[J]. 應用與環(huán)境生物學報， 26（5）： 1235-1243.]

（責任編輯?何永艷）