木薯PG基因家族鑒定及表達分析

摘""要：多聚半乳糖醛酸酶（polygalacturonase，"PG）在果膠的分解過程中起著關(guān)鍵作用，從而導(dǎo)致細胞壁結(jié)構(gòu)的變化，參與植物的生長發(fā)育和應(yīng)對逆境脅迫等多個階段。盡管PG基因在其他植物中已有研究，但在木薯中尚未見對該基因家族成員鑒定和功能研究的報道。本研究從木薯基因組中鑒定出89個MePGs，編碼的氨基酸數(shù)量在183～808之間，分子量介于19.75～87.07"kDa之間，理論等電點（pI）在4.64～9.71之間，大部分成員預(yù)測定位于細胞膜。染色體定位分析顯示MePGs不均勻地分布于17條染色體上。基于進化關(guān)系，MePGs被分為A～G共7個亞組，亞組成員基因結(jié)構(gòu)相似，且存在明顯的串聯(lián)重復(fù)現(xiàn)象；共線性分析結(jié)果表示，相比擬南芥，木薯PG家族基因的親緣性與橡膠更接近。MePGs啟動子區(qū)域富含光響應(yīng)、激素響應(yīng)與逆境響應(yīng)元件。MePGs基因具有組織特異性，與木薯生長發(fā)育有關(guān)。在木薯采后變質(zhì)過程中，MePGs表現(xiàn)出相似的表達模式，特別是MePG20、MePG21、MePG25、MePG64和MePG72，均呈現(xiàn)先上升后下降的表達趨勢，表明這些基因可能在木薯塊根采后變質(zhì)過程中先響應(yīng)逆境脅迫，高表達以分解果膠，隨后隨著細胞壁的水解完成而降低表達，這表明PG基因在逆境脅迫中發(fā)揮重要作用。綜上所述，MePG基因家族成員可能通過片段復(fù)制和內(nèi)含子減少的方式進化，并可能通過感知不同種類的信號而行使不同的功能，從而形成多樣化的表達模式。本研究為進一步探索MePGs在木薯塊根采后變質(zhì)機制中的作用提供基礎(chǔ)。

關(guān)鍵詞：木薯；PG基因家族；表達分析中圖分類號：S31""""""文獻標(biāo)志碼：A

Identification"and"Expression"Analysis"of"the"PG"Gene"Family"in"Cassava

TANG"Xiangning1，2，3，"WANG"Xiaotong1，2，3，"WANG"Luhua2，3，4，"Li"Ruimei2，3，"Yao"Yuan2，3，"WANG"Yajie2，3，"GUO"Jianchun2，3*，"LIU"Jiao2，3*

1."School"of"Life"and"Health"Sciences，"Hainan"University，"Haikou，"Hainan"570228，"China;"2."Institute"of"Tropical"Biotechnology，"Chinese"Academy"of"Tropical"Agriculturalnbsp;Sciences"/"National"Key"Laboratory"of"Tropical"Crop"Breeding，"Haikou，"Hainan"571101，"China;"3."Sanya"Research"Institute，"Chinese"Academy"of"Tropical"Agricultural"Sciences，"Sanya，"Hainan"572025，"China;"4."School"of"Tropical"Agriculture"and"Forestry，"Hainan"University，"Danzhou，"Hainan"571737，"China

Abstract："Polygalacturonase"（PG）"plays"a"crucial"role"in"plant"growth"and"development，"and"in"stress"responses，"by"participating"in"the"degradation"of"pectin"and"altering"cell"wall"structure."Although"studies"on"PG"genes"have"been"conducted"in"various"plants，"there"is"a"lack"of"research"on"the"identification"and"functional"characterization"of"the"PG"gene"family"in"cassava"（Manihot"esculenta"Crantz）."This"study"identified"a"total"of"89"members"of"the"MePG"family"in"the"cassava"genome，"encoding"proteins"with"183"to"808"amino"acids，"molecular"weights"ranging"from"19.75"to"87.07"kDa，"and"theoretical"pI"values"between"4.64"and"9.71."Most"of"the"family"members"are"predicted"to"be"localized"in"the"cell"membrane."Chromosome"mapping"analysis"revealed"that"MePG"family"members"were"unevenly"distributed"across"17"chromosomes."Based"on"evolutionary"relationships，"MePGs"were"classified"into"seven"subgroups"（A"to"G），"with"similar"gene"structures"within"subgroups"and"evidence"of"tandem"duplication."Comparative"analysis"indicated"that"the"PG"family"genes"in"cassava"were"more"closely"related"to"rubber"than"to"Arabidopsis."The"promoter"regions"of"MePG"genes"were"enriched"with"elements"responsive"to"light，"hormones，"and"stress."MePG"genes"exhibit"tissue"specificity"and"were"associated"with"cassava"growth"and"development."During"postharvest"deterioration"of"cassava，"MePGs"displayed"similar"expression"patterns，"particularly"MePG20，"MePG21，"MePG25，"MePG64"and"MePG72，"which"showed"an"initial"increase"followed"by"a"decrease"in"expression，"suggesting"that"the"genes"may"initially"respond"to"stress，"upregulate"to"break"down"pectin，"and"then"downregulate"as"cell"wall"hydrolysis"is"completed."This"implies"a"significant"role"for"PG"in"stress"responses."Collectively，"the"findings"suggest"that"MePG"gene"family"members"may"have"evolved"through"segmental"duplication"and"intron"loss"and"may"perform"different"functions"by"sensing"various"types"of"signals，"leading"to"diverse"expression"patterns."This"would"study"provide"a"foundation"for"further"exploration"of"the"role"of"MePGs"in"the"postharvest"deterioration"mechanism"of"cassava"roots.

Keywords："cassava;"PG"gene"family;"expression"analysis

DOI："10.3969/j.issn.1000-2561.2025.06.004

植物細胞壁是細胞的外部支撐結(jié)構(gòu)，主要由纖維素、半纖維素和果膠構(gòu)成，為細胞提供結(jié)構(gòu)支持和保護功能[1]。果膠在植物細胞壁中扮演多重角色，包括維持細胞結(jié)構(gòu)，促進細胞間粘連，增強細胞壁的延展性，提供生長和發(fā)育信號以及參與植物的防御響應(yīng)[2]。

果膠的生物合成涉及多種糖基轉(zhuǎn)移酶和蛋白質(zhì)復(fù)合體[3]。果膠由同源半乳糖醛酸（homog alac tu ronan，"HG）、鼠李半乳糖醛酸聚糖Ⅰ（rhamn o g ala cturonan"I，"RG-I）、鼠李半乳糖醛酸聚糖Ⅱ（rhamnogalacturonan"II，"RG-II）、木糖半乳糖醛酸聚糖（xylogalacturonan，"XGA）等組成[4-6]，其中HG約占果膠總量的65%[7]。

多聚半乳糖醛酸酶（polygalacturonase，"PG）是一類催化果膠降解的酶，是HG的主要酶之一，參與細胞壁中果膠的降解，影響細胞壁穩(wěn)定性，在植物各生長階段以及不同組織中扮演著關(guān)鍵角色[8]。

植物中PG蛋白通常包含4個保守的結(jié)構(gòu)域。其中，Ⅰ域和Ⅱ域分別包含核心序列SPNTDG和GDDC，這2個序列中包含3個天冬氨酸殘基（D），這些殘基參與構(gòu)成酶的催化活性位點。Ⅲ域由序列CGPGHG組成，其中的組氨酸（H）在催化過程中發(fā)揮重要作用。Ⅳ域的序列為RIK，該結(jié)構(gòu)域涉及與底物羧基端的離子相互作用。這些結(jié)構(gòu)域的保守性和功能多樣性是PG蛋白在植物細胞壁降解過程中發(fā)揮作用的基礎(chǔ)[9]。在這4個保守結(jié)構(gòu)域中，半胱氨酸殘基（C）均有較好的保守性，推測其可能參與維持蛋白質(zhì)的三級結(jié)構(gòu)[10]。

木薯（Manihot"esculenta"Crantz）是以富含淀粉的塊根為主要利用形式的重要糧食和能源作物，其塊根在采收后極易發(fā)生生理性變質(zhì)（post h ar vest"physiological"deterioration，PPD）[11]。木薯變質(zhì)分為初發(fā)性變質(zhì)和次發(fā)性變質(zhì)2個階段，初發(fā)性變質(zhì)取決于內(nèi)源酶作用，而次發(fā)性變質(zhì)主要由微生物感染引起。在次發(fā)性變質(zhì)過程中，作為植物細胞第一道屏障的細胞壁結(jié)構(gòu)發(fā)生變化，微生物隨之侵入引發(fā)次級變質(zhì)，繼而發(fā)生嚴(yán)重的氧化損傷和細胞死亡，最終導(dǎo)致整個塊根腐爛[12]。在這個過程中，細胞壁結(jié)構(gòu)的變化可能與PG家族基因的活性密切相關(guān)。

目前PG基因家族在許多物種中都進行了全基因組鑒定及分析，如榴蓮[13]、甘薯[14]和蘋果[15]等，但在木薯中暫未見報道。本研究旨在木薯基因組數(shù)據(jù)庫中篩選鑒定木薯PG家族成員，對其理化性質(zhì)、基因結(jié)構(gòu)和啟動子序列等進行分析，并且分析其在木薯不同組織部位和塊根采后變質(zhì)過程中的表達情況，以探討PG基因在木薯塊根采后變質(zhì)過程中的作用。

1""材料與方法

1.1""材料

木薯種植基地位于海南省三亞市南濱農(nóng)場，年降水量1400～1600"mm，年平均氣溫24～26℃。選擇數(shù)棵生長成熟（約150"d）的華南8號（SC8）木薯作為試驗材料。將未損傷的木薯塊根完整挖出，置于實驗室室溫存放。每次選擇完整的塊根在0、1、3、5、7、10、14、21"d進行切片取樣，并存放在?80"℃冰箱中保存?zhèn)溆谩?/p>

木薯（Manihotnbsp;esculenta"Crantz）基因組及注釋數(shù)據(jù)下載于NCBI數(shù)據(jù)中心（https：//www.ncbi."nlm.nih.gov/），擬南芥[Arabidopsis"thaliana"（L.）"H eynh.]和橡膠[Hevea"brasiliensis"（Willd."ex"A."Juss.）"Müll."Arg.]PG基因家族序列下載于NC BI（https：//www.ncbi.nlm.nih.gov/）和Phytozome（ht tps：//phytozome-next.jgi.doe.gov/）數(shù)據(jù)庫。

1.2""方法

1.2.1""木薯PG基因家族成員鑒定""通過TBtools軟件獲取木薯蛋白序列，在擬南芥的蛋白序列中搜索與木薯PG蛋白相似的蛋白序列，將這些蛋白放到swissprot（http：//www.gpmaw.com/html/s w issprot.html）數(shù)據(jù)庫中進行驗證，篩選出木薯多聚半乳糖醛酸酶（MePG）家族基因的候選序列。通過InterPro（https：//www.ebi.ac.uk/interpro/）在線軟件和TBtools軟件的Advanced"Hmmer"Search插件，將候選序列與已知的PG保守結(jié)構(gòu)域（PF00295）進行篩選，以確定最終的MePG基因家族成員，并對篩選出的基因進行重命名。

1.2.2""蛋白質(zhì)特性分析""利用TBtools軟件的Protein"Parameter"Calc插件對MePG蛋白的特性進行預(yù)測，包括氨基酸數(shù)量、分子量（MW）和等電點（pI）。通過PSORT"Prediction（https：//psort.hgc."jp/form.html）在線軟件進行亞細胞定位的預(yù)測。

1.2.3""木薯PG蛋白多序列比對""利用DNAM A N軟件對PG蛋白家族進行蛋白多序列比較。

1.2.4""木薯PG蛋白系統(tǒng)發(fā)育樹的構(gòu)建""使用TBtools軟件的One"Step"Build"a"ML"Tree插件創(chuàng)建MePG蛋白序列的系統(tǒng)發(fā)育樹。

1.2.5""木薯PG家族基因結(jié)構(gòu)和motif分析""利用NCBI"batch"wab"CD-search"tool（https：//www.ncbi."nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi）在線軟

件進行保守結(jié)構(gòu)域分析，通過MEME（https：//"me me-suite.org/meme/tools/meme）在線軟件分析MePG蛋白的保守基序；利用TBtools軟件進行可視化處理。

1.2.6""染色體分布和物種間共線性""使用TBtools軟件確定89個MePG基因的染色體位置。使用TBtools軟件的One"Step"MCScanX插件進行木薯與擬南芥和橡膠中的基因同源性分析。

1.2.7""MePG基因家族啟動子順式作用元件分析及功能預(yù)測""利用TBtools軟件提取MePG基因家族序列上游長度為2000"bp的啟動子序列，利用PlantCARE（http：//bioinformatics.psb.ugent.be/"w e b tools/plantcare/html/）軟件進行預(yù)測分析，并通過TBtools軟件進行可視化分析。

1.2.8""MePG家族成員的表達模式分析""根據(jù)本實驗室的SC8木薯轉(zhuǎn)錄組數(shù)據(jù)，分析MePG基因在各組織部位以及將SC8木薯采后變質(zhì)21"d的FPKM值，利用ChiPlot（https：//www.chiplot."on l ine/）在線軟件繪制MePG基因的表達熱圖。選取MePG20/21/25/64/72基因，設(shè)計qRT-PCR特異性引物（表1），由北京擎科生物科技有限公司（海口）合成引物。使用M5"HiPer"Plant"RNeasy"Complex"Mini"Kit“繁可簡”多糖多酚植物RNA提取試劑盒（北京聚合美生物科技有限公司）提取塊根RNA，使用MonScript?"RTIII"All-in-One"Mix"with"dsDNase（莫納生物科技有限公司）進行逆轉(zhuǎn)錄。使用Taq-HS"SYBR?"Green"qPCR"Premix"（Universal）（江蘇愚公生命科技有限公司）和上海力康HealForce的CG-05儀器進行qRT-PCR驗證。使用2?ΔΔCt法計算基因的相對表達量。

2""結(jié)果與分析

2.1""MePG基因家族全基因組鑒定與理化分析

利用生物信息學(xué)方法分析木薯基因組數(shù)據(jù)庫，共鑒定到89個PG基因。根據(jù)染色體位置分

布，將89個MePG基因依次命名為MePG1～"MePG89。氨基酸序列分析表明，MePG蛋白質(zhì)的氨基酸數(shù)目為183（MePG75）～808（MePG22），平均長度為442"aa；分子量為19.75（MePG75）～

87.08（MePG22）kDa，平均值為47.79"kDa；理論等電點為4.64（MePG21）～9.71（MePG32）（表2）。亞細胞定位預(yù)測結(jié)果顯示82個MePG家族成員定位于細胞膜，其余成員可能分布在細胞壁、葉綠體、細胞質(zhì)或細胞核中。

MePG的氨基酸序列比對結(jié)果如圖1所示，89個MePGs中有33個MePGs的氨基酸序列完整包含了4個典型的PG保守結(jié)構(gòu)域Ⅰ（序列為SPNTDGI）、Ⅱ（序列為GDDC）、Ⅲ（序列為CGPGHG）和Ⅳ（序列為RIK），如MePG1、MePG2、MePG3、MePG4等。部分MePGs的氨基酸序列存在丟失或缺少的現(xiàn)象，其中42個MePGs的氨基酸序列缺少結(jié)構(gòu)域Ⅰ，19個MePGs的氨基酸序列缺少結(jié)構(gòu)域Ⅱ，38個MePGs的氨基酸序列缺少結(jié)構(gòu)域Ⅲ，10個MePGs的氨基酸序列缺少結(jié)構(gòu)域Ⅳ。MePG17、MePG18、MePG20的氨基酸序列同時缺少4個保守結(jié)構(gòu)域，并伴有氨基酸替換或缺失，推測這些PG基因可能具有更高的進化程度。

2.2""MePG基因家族進化樹分析

為研究木薯PG基因的進化關(guān)系，本研究采用鄰接法（NJ），構(gòu)建木薯和擬南芥（Arabidopsis"thaliana）的PG基因家族成員蛋白序列（68個AtPG，89個MePG）進化樹。如圖2所示，2個物種的157個PG蛋白被分為7個亞族（Group"A～Group"G）。A至G亞族分別包含9、6、21、27、16、5、5個木薯PG。2個物種在A亞族、B亞族和E亞族中的基因數(shù)量大致相同，2個物種的PG蛋白在每個亞族中不均勻分布。亞族的劃分基于蛋白序列的相似度和進化關(guān)系，不同亞族的基因通常具有不同的功能。

2.3""MePG基因保守基序及基因結(jié)構(gòu)分析

在基因家族進化過程中，通常伴有外顯子的缺失和增加，為了理解木薯PG基因的特性，通過TBtools軟件對木薯89個MePG基因的外顯子和內(nèi)含子結(jié)構(gòu)進行分析，并繪制基因結(jié)構(gòu)圖。結(jié)果表明（圖3A），同一亞族中的PG基因具有相似的基因結(jié)構(gòu)，其內(nèi)含子的數(shù)量接近，而不同亞族（如Group"F和Group"G）包含基因之間內(nèi)含子數(shù)量差異較大。所有MePG基因中外顯子數(shù)量最少為3個，而位于Group"F中的MePG60外顯子的數(shù)量最多（18個）。以上結(jié)果表明，在木薯PG基因進化過程中，外顯子的丟失、獲得與PG基因的功能多樣密切相關(guān)，這是因為基因家族內(nèi)部的基因在分子功能上是冗余的，它們可能經(jīng)歷了進化選擇過程，基因中斷由內(nèi)含子插入是真核基因的特征之一。

木薯PG的保守基序分析結(jié)果如圖3B所示，親緣關(guān)系越近的成員其基序結(jié)構(gòu)越相似。大部分MePGs都包含有Motif"10和Motif"6，且均以N→C端按序排列。其中Motif"1的序列包含“SPNTDGI”（保守結(jié)構(gòu)域Ⅰ）與“GDDC”（保守結(jié)構(gòu)域Ⅱ），Motif"3和Motif"5的序列分別包含“CGPGHG”（保守結(jié)構(gòu)域Ⅲ）、“RIK”（保守結(jié)構(gòu)域Ⅳ）。保守結(jié)構(gòu)域Ⅰ包含1個天冬氨酸殘基（D），該殘基參與酶的催化活性位點，直接影響PG酶的催化效率，缺少該結(jié)構(gòu)域的基因可能催化活性較低或功能不同。保守結(jié)構(gòu)域Ⅱ同樣包含天冬氨酸殘基，參與底物的結(jié)合和催化過程，缺少該結(jié)構(gòu)域的基因可能無法有效結(jié)合果膠底物，導(dǎo)致酶活性降低。保守結(jié)構(gòu)域Ⅲ中的組氨酸（H）在催化過程中發(fā)揮重要作用，尤其是在果膠的水解過程中，缺少該結(jié)構(gòu)域的基因可能在果膠降解過程中效率較低。保守結(jié)構(gòu)域Ⅳ可能與底物羧基端的離子相互作用有關(guān)，影響PG酶與果膠的結(jié)合能力，缺少該結(jié)構(gòu)域的基因可能無法有效結(jié)合果膠，導(dǎo)致酶活性下降。

從進化關(guān)系上看，親緣關(guān)系較近的同一亞族的PG基因在基因結(jié)構(gòu)組成和蛋白結(jié)構(gòu)域分布上也具有較高保守性。如圖3C所示，68個（76.4%）木薯PG家族成員包含1個共同的結(jié)構(gòu)域PL-6"superfamily，16個（17.9%）木薯PG家族成員包含1個共同的結(jié)構(gòu)域Pgu1"superfamily，4個（4.4%）木薯PG家族成員包含結(jié)構(gòu)域Pectate_"lyase_3"superfamily，但MePG17并不包含這3種結(jié)構(gòu)域中的任何一種。

2.4""MePG染色體定位分析

用TBtools軟件中的染色體定位工具對89個MePG基因在染色體上的位置分析發(fā)現(xiàn)，89個MePG基因不均勻地分布在17條染色體上，其中第1、5和9號染色體上的基因數(shù)量最多，各含有11個MePG基因，而第4、6號染色體上的基因數(shù)量最少，僅含有1個MePG基因（圖4）。這種不均勻的分布可能反映了基因復(fù)制事件和染色體重排。MePG88和MePG89由于基因組未完全組裝，僅定位到基因組的特定區(qū)域。MePG基因在染色體上的分布模式表明，這些PG基因可能通過片段復(fù)制和串聯(lián)復(fù)制的方式進化，且不同染色體上的基因可能具有不同的功能或表達模式。

2.5""MePG基因的物種間共線性分析

為了進一步了解木薯PG基因家族的系統(tǒng)發(fā)育機制，本研究分析了木薯與擬南芥、橡膠的種間共線性關(guān)系。如圖5所示，在木薯與擬南芥共有51對PG同源基因?qū)Γ臼砼c橡膠共有70對PG同源基因?qū)Γ砻髂臼砼c橡膠的親緣關(guān)系更近。

2.6""MePG基因啟動子順式作用元件分析

為了弄清MePG基因的潛在功能和調(diào)控機制，本研究截取89個木薯PG基因上游2000"bp的啟動子序列，利用PlantCARE軟件對MePGs啟動子區(qū)域的順式作用元件進行分析。結(jié)果（圖6）顯示，大部分基因都包含與植物激素和非生物脅迫響應(yīng)相關(guān)的作用元件，這些元件與木薯的生長發(fā)育息息相關(guān)。所含順式作用元件中，最多的是光響應(yīng)元件（671個），其次是茉莉酸甲酯響應(yīng)元件（185個）和脫落酸響應(yīng)元件（140個），其余作用元件除損傷響應(yīng)元件（6個）和柵欄葉肉細胞分化響應(yīng)元件（4個）外，數(shù)量分布范圍在20～68個之間。不同MePG基因家族成員啟動子區(qū)域所含的順式作用元件類型差異較大，大部分MePGs的啟動子區(qū)域含有光響應(yīng)元件，表明這些基因可能受到光信號的調(diào)控，尤其是在光合作用活躍的組織中。茉莉酸甲酯和脫落酸響應(yīng)元件藍色線條代表木薯和擬南芥、橡膠的共線性PG基因?qū)Γ疑€條表示木薯和擬南芥、橡膠的同源基因。

在MePGs的啟動子區(qū)域中非常豐富，表明這些基因可能在木薯逆境響應(yīng)中發(fā)揮重要作用，脫落酸響應(yīng)元件可能在干旱或鹽脅迫條件下調(diào)控MePG基因的表達。部分MePG基因（如MePG3、MePG20、MePG26等）的啟動子區(qū)域含有損傷響應(yīng)元件，表明這些基因可能在機械損傷或病原體侵染時被激活，參與細胞壁的修復(fù)和防御過程。如赤霉素、生長素等激素響應(yīng)元件也存在于部分MePG基因的啟動子區(qū)域，表明這些基因可能在木薯的生長發(fā)育過程中發(fā)揮調(diào)控作用。

2.7""MePG基因表達模式分析

本研究分析了69個MePG基因在木薯愈傷組織、體胚、嫩葉、成熟葉、莖、須根、塊根、塊根韌皮部和木質(zhì)部等9個組織器官中的表達模式，結(jié)果如圖7所示。根據(jù)MePGs的表達譜可將其分為3個組：A組，21個基因在大多數(shù)組織中顯示出高表達水平；B組，37個基因在特定組織中高度轉(zhuǎn)錄；C組，11個基因在大多數(shù)組織中顯示出低表達水平。B組中的MePGs基本只在愈傷組織和體胚中表達，而A組中的MePG20、MePG72、MePG25、MePG16、MePG21、MePG5等基因在塊根中表達量較高。此外，MePG15、MePG31為愈傷組織特異表達基因，MePG8、MePG9、MePG32為體胚特異表達基因。

目前研究表明PG參與細胞壁結(jié)構(gòu)分解，木薯塊根采后變質(zhì)中細胞壁結(jié)構(gòu)的變化可能與PG家族基因的活性密切相關(guān)。根據(jù)本實驗室前期對木薯果膠和PG活性的研究，發(fā)現(xiàn)同時期PG活性低的木薯細胞壁可以保留更完整的酯化果膠[16]。

本研究進一步對SC8木薯塊根采后變質(zhì)過程中MePGs的表達情況進行分析發(fā)現(xiàn)：一部分MePGs在采后變質(zhì)中幾乎不表達；一部分MePGs的表達量無明顯變化，如MePG16、MePG65、MePG13和MePG12等；MePG20、MePG21、MePG25、MePG64、MePG72在木薯塊根采后變質(zhì)21"d的表達量下調(diào)（圖8）。

如圖9所示，選擇8個時間段的木薯采后塊根進行取樣，在每顆塊根中部選擇大小一致位置作為取樣點，切片并迅速研磨提取RNA。可以發(fā)現(xiàn)，SC8木薯整根塊根在采后前3"d的薯塊顏色、硬度以及韌皮部形態(tài)變化不大，3～5"d開始發(fā)生明顯褐化、薯塊變軟等現(xiàn)象。而到14"d，薯塊褐變面積明顯增大。到21"d時，薯塊褐變顏色加深，薯塊木質(zhì)部出現(xiàn)開裂，薯塊進一步變軟。

進一步驗證MePG在塊根采后變質(zhì)過程中的表達模式。選擇MePG20、MePG21、MePG25、MePG64、MePG72進行qPCR分析，結(jié)果如圖10所示。在SC8變質(zhì)21"d的過程中，MePG20基因除第1天有短暫上升外，后面時間段基本呈下調(diào)的表達趨勢，在第2周（7～14"d）的表達量無顯著差異，第21天的表達量最低；MePG21基因在第5天的表達量升高至最高，之后隨著變質(zhì)時間的延長其表達量持續(xù)降低；MePG25基因表達量在前3"d短暫降低，在5～7"d升高，7"d后其表達量開始隨著變質(zhì)時間降低；MePG64基因表達量隨變質(zhì)時間的延長出現(xiàn)降低-升高-降低-升高-迅速降低的變化；MePG72基因在第1天的表達量升至最高，3"d后降低，然后在第5天升高，后續(xù)隨著變質(zhì)時間持續(xù)降低。5個MePG基因的相對表達量變化模式與轉(zhuǎn)錄組相符。

3""討論

作為果膠修飾酶的關(guān)鍵成員，植物PG基因家族的功能已受到廣泛關(guān)注。PG基因在植物不同

發(fā)育階段和組織中的表達模式表明其在細胞壁水解中發(fā)揮重要作用。本研究在木薯基因組中鑒定出89個MePG基因家族成員，這些成員在理化性質(zhì)如肽鏈長度、相對分子量和等電點等方面表現(xiàn)出顯著差異，與甘薯PG基因家族的研究結(jié)果相似[14]。多種高等植物中均已鑒定出PG基因家族成員，如擬南芥（68個）、榴蓮（57個）[13]、甘薯（103個）[14]、蘋果（85個）[15]、桃（84個）[17]、水稻（44個）[18]、甘藍（99個）[19]、小麥（113個）[20]、蓮霧（25個）[21]等。系統(tǒng)進化分析顯示，木薯與擬南芥、榴蓮、甘薯和桃等的PG基因家族均分為7個亞族，亞族內(nèi)基因結(jié)構(gòu)和保守基序相似。木薯中PG基因的數(shù)量與蘋果、桃接近，而與甘薯、水稻、小麥、蓮霧等有較大差異，這可能反映了不同植物物種在進化過程中基因復(fù)制事件和環(huán)境適應(yīng)性選擇的影響。

植物激素如水楊酸、脫落酸、赤霉素、生長素、獨腳金內(nèi)酯、茉莉酸甲酯等在植物逆境響應(yīng)中起到信號傳導(dǎo)作用，介導(dǎo)逆境脅迫下的代謝物質(zhì)合成[22]。木薯與榴蓮中PG基因家族啟動子含有相似的激素與逆境響應(yīng)元件[13]，推測多數(shù)PG蛋白能夠?qū)Ψ巧飸?yīng)激做出反應(yīng)，特別是脫落酸和茉莉酸甲酯，作為植物逆境生長和發(fā)育的重要調(diào)控激素，其順式作用元件在MePG基因家族中非常豐富，表明其可能在調(diào)節(jié)木薯生長和發(fā)育中發(fā)揮重要作用，并參與抗病過程[22]。此外，MePG3、MePG20、MePG26、MePG41、MePG55和MePG63的啟動子區(qū)域存在損傷響應(yīng)元件，這可能是研究植物機械損傷和病原體侵染的關(guān)鍵點。本研究關(guān)注的幾個在木薯采后過程中表達量顯著變化的基因如MePG20、MePG21和MePG25等，其啟動子區(qū)域富含光響應(yīng)元件、茉莉酸甲酯響應(yīng)元件、脫落酸響應(yīng)元件、低溫響應(yīng)元件以及損傷響應(yīng)元件。這些元件的存在表明木薯PG基因可能在光調(diào)控的生長發(fā)育過程中發(fā)揮作用，并在多種逆境脅迫條件下被激活，參與細胞壁降解和防御機制。后續(xù)研究應(yīng)該關(guān)注這些基因在逆境脅迫條件下的表達調(diào)控機制，以及其在木薯塊根采后變質(zhì)中的具體功能。

PG家族成員在不同組織部位具有不同的表達模式，部分基因在采后變質(zhì)過程中的表達量呈變化顯著，表明這些PG基因可能響應(yīng)其他生長機制的調(diào)控。MePG20、MePG21、MePG25、MePG64、MePG72等基因。在采后變質(zhì)過程中的表達量先上升后下降，表明PG酶活性可能在果實成熟軟化的某個階段達到高峰，隨后下降，這與葡萄成熟中的PG活性變化趨勢相似[23]。這一過程中涉及到復(fù)雜的激素和酶的相互作用，以及細胞壁結(jié)構(gòu)的變化。MePGs的表達變化說明其在植物細胞壁的降解中扮演重要角色，在木薯采后階段響應(yīng)逆境脅迫中起到水解果膠改變細胞壁的作用。

MePG20、MePG21、MePG25、MePG64、MePG72等基因在采后變質(zhì)過程中表達量變化顯著，表明其可能在木薯塊根采后變質(zhì)中起重要作用。qPCR結(jié)果進一步驗證了這些基因的表達模式，表明其在木薯塊根細胞壁軟化和變質(zhì)過程中發(fā)揮關(guān)鍵作用。

綜上所述，本研究從木薯全基因組中鑒定出89個MePG基因，并分析其理化性質(zhì)、保守結(jié)構(gòu)域、保守基序、順式作用元件等和其在不同組織部位及塊根采后變質(zhì)過程的表達模式。本研究初步認(rèn)定MePG可能在木薯塊根響應(yīng)采后變質(zhì)具有重要作用，為進一步研究木薯PG基因?qū)δ婢稠憫?yīng)機制奠定基礎(chǔ)。

參考文獻

[1]"COSGROVE"D"J."Structure"and"growth"of"plant"cell"walls[J]."Nature"Reviews"Moleculra"Cell"Biology，"2024，"25："340-358.

[2]"JOBERT"F，"YADAV"S，"ROBERT"S，"ROBERT"S."Auxin"as"an"architect"of"the"pectin"matrix[J]."Journal"of"Experimental"Botany，"2023，"74（22）："6933-6949.

[3]"ATMODJO"M"A，"HAO"Z，"MOHNEN"D."Evolving"views"of"pectin"biosynthesis[J]."Annual"Review"of"Plant"Biology，"2013，"64（1）："747-779.

[4]"MOHNEN"D."Pectin"structure"and"biosynthesis[J]."Current"Opinion"in"Plant"Biology，"2008，"11（3）："266-277.

[5]"LI"L，"GAO"X"L，"LIU"J"G，"CHITRAKAR"B，"WANG"B，"WA NG"Y"C."Hawthorn"pectin："extraction，"function"and"utilization[J]."Current"Research"in"Food"Science，"2021，"4："429-435.

[6]"LINARES-GARCíA"J"A，"RAMOS-RAMíREZ"E"G，"SALA ZAR- MONTOYA"J"A."Viscoelastic"properties"and"textural"chara cterisation"of"high"methoxyl"pectin"of"hawthorn"（Cra ta e gus"pubescens）"in"a"gelling"system[J]."International"Journal"of"Food"Science"amp;"Technology，"2015，"50（6）："1484-1493.

[7]"WOLF"S，"MOUILLE"G，"PELLOUX"J."Homogalacturonan"methyl-esterification"and"plant"development[J]."Molecular"Plant，"2009，"2（5）："851-860.

[8]"鄭秋萍，"林育釗，"李美玲，"王慧，"林河通."果實采后軟化的影響因素及抑制技術(shù)研究進展[J]."亞熱帶農(nóng)業(yè)研究，"2019，"15（4）："262-270.ZHENG"Q"P，"LIN"Y"Z，"LI"M"L，"WANG"H，"LIN"H"T."Research"progress"on"the"influencing"factors"and"inhibiting"technologies"of"postharvest"fruit"softening[J]."Subtropical"Agriculture"Research，"2019，"15（4）："262-270."（in"Chineses）

[9]"LU"L，"HOU"Q"C，"WANG"L"L，"WANG"L"L，"ZHANG"T"Y，"ZHAO"W，"YAN"T"W，"ZHAO"L"N，"LI"J"P，"WAN"X"Y."Genome-wide"identification"and"characterization"of"polygalacturonase"gene"family"in"maize"（Zea"mays"L.）[J]."International"Journal"of"Molecular"Sciences，"2021，"22（19）："10722.

[10]"葛廷."柑橘PG基因的鑒定及其與幼果脫落的相關(guān)性分析[D]."重慶："西南大學(xué)，"2020.GE"T."Identification"of"Citrus"PG"genes"and"their"roles"in"fruitlet"abscission[D]."Chongqing："Southwest"University，"2020."（in"Chineses）

• 張鵬."能源木薯種質(zhì)資源面臨的問題與解決策略[J]."生物產(chǎn)業(yè)技術(shù)，"2008（5）："25-30.ZHANG"P."Problems"and"solutions"of"energy"cassava"germplasm"resources[J]."Biotechnology"amp;"Business，"2008（5）："25-30."（in"Chineses）
• LIU"S，"ZAINUDDIN"I"M，"VANDERSCHUREN"H，"JAMES"D，"BEECHING"J"R."RNAi"inhibition"of"feruloyl"CoA"6′-hyd r oxylase"reduces"scopoletin"biosynthesis"and"post-har vest"physiological"deterioration"in"cassava"（Manihot"esculenta"Cr antz）"storage"roots[J]."Plant"Molecular"Biology，"2017，"94（1/2）："185-195.
• 劉玲玲."榴蓮PG基因家族鑒定及其裂果的果膠變化分析[D]."合肥：nbsp;安徽農(nóng)業(yè)大學(xué)，"2022.LIU"L"L."Identification"of"PG"gene"family"in"durian"and"pectin"changes"during"fruit"dehiscence[D]."Heifei："Anhui"Agricultural"University，"2022."（in"Chineses）
• HE"P，"ZHANG"J，"LV"Z，"CUI"P，"XU"X"M，"GEORGE"M"S，"LU"G"Q."Genome-wide"identification"and"expression"analysis"of"the"polygalacturonase"gene"family"in"sweetpotato[J]."BMC"Plant"Biology，"2023，"23（1）："300.
• 陳鴻飛，"邵紅霞，"樊勝，"馬娟娟，"張東，"韓明玉."蘋果全基因組多聚半乳糖醛酸酶基因家族的鑒定及進化分析[J]."園藝學(xué)報，"2016，"43（10）："1863-1877.CHEN"H"F，"SHAO"H"X，"FAN"S，"MA"J"J，"ZHANG"D，"HAN"M"Y."Identification"and"phylogenetic"analysis"of"the"polygalacturonase"gene"family"in"apple[J]."Acta"Horticulturae"Sinica，"2016，"43（10）："1863-1877."（in"Chineses）
• 周楊驕."木薯MePMEI1對植物細胞壁果膠的影響及其響應(yīng)重金屬"Pb、Cd脅迫的功能研究[D]."海口："海南大學(xué)，"2023.ZHOU"Y"J."Effect"of"cassava"MePMEI1"on"plant"cell"wall"pectin"and"its"function"in"response"to"heavy"metal"Pb"and"Cd"stresses[D]."Haikou："Hainan"University，"2023."（in"Chineses）
• 霍如雪，"劉振寧，"楊青，"王光全."桃PG基因家族的鑒定與分析[J]."江蘇農(nóng)業(yè)科學(xué)，"2016，"44（6）："33-40.HUO"R"X，"LIU"Z"N，"YANG"Q，"WANG"G"Q."Identification"and"analysis"of"PG"gene"family"in"peach[J]."Jiangsu"Agricultural"Sciences，"2016，"44（6）："33-40."（in"Chineses）
• 鄭銀珍."水稻PLL和PG基因家族的鑒定及功能研究[D]."武漢："武漢大學(xué)，"2018.ZHENG"Y"Z."Genomic"identification"of"PLL"and"PG"gene"families"and"functional"analysis"in"rice[D]."Wuhan："Wuhan"University，"2018."（in"Chineses）
• LYU"M，"IFTIKHAR"J，"GUO"R，"WU"B"H，"CAO"J"S."Patterns"of"expansion"and"expression"divergence"of"the"polygalacturonase"gene"family"in"Brassica"oleracea[J]."International"Journal"of"Molecular"Sciences，"2020，"21（16）："5706.

[20]"YE"J，"YANG"X，"YANG"Z，"NIU"F"Q，"CHEN"Y"R，"ZHANG"L"L，"SONG"X"Y."Comprehensive"analysis"of"polygalacturonase"gene"family"highlights"candidate"genes"related"to"pollen"development"and"male"fertility"in"wheat"（Triticum"aestivum"L.）[J]."Planta，"2020，"252（2）："31-48.

[21]"陳迪飛，"魏秀清，"許玲，"許家輝，"曾黎輝."蓮霧PG基因家族全基因組鑒定及表達分析[J]."果樹學(xué)報，"2022，"39（4）："548-563.CHEN"D"F，"WEI"X"Q，"XU"L，"XU"J"H，"ZENG"L"H."Genome-wide"identification"and"expression"analysis"of"PG"gene"family"in"wax"apple"[Syzygium"samarangense"（Bl.）"Merr."et"Perry][J]."Journal"of"Fruit"Science，"2022，"39（4）："548-563."（in"Chineses）

[22]"CANNON"S"B，"MITRA"A，"BAUMGARTEN"A，"YOUNG"N"D，"MAY"G."The"roles"of"segmental"and"tandem"gene"duplication"in"the"evolution"of"large"gene"families"in"Arabidopsis"thaliana[J]."BMC"Plant"Biology，"2004，"4（1）："4.

[23]"海龍飛，"栗溫新，"李志謙，"李猛，"陳超陽，"魏紅麗，"鄒東方，"何暢，"馮建燦，"葉霞."軟/硬肉葡萄果實細胞壁結(jié)構(gòu)、組分及降解酶活性的變化[J]."果樹學(xué)報，"2023，"40（4）："690-698.HAI"L"F，"LI"W"X，"LI"Z"Q，"LI"M，"CHEN"C"Y，"WEI"H"L，"ZOU"D"F，"HE"C，"FENG"J"C，"YE"X."Variations"in"cell"wall"microstructure"and"components"and"activities"of"their"degradation"enzymes"in"grapes"with"soft"or"hard"textures[J]."Journal"of"Fruit"Science，"2023，"40（4）："690-698."（in"Chineses）