桑樹(shù)GASA基因家族全基因組鑒定及表達(dá)模式分析

摘""要：GASA（gibberellic"acid-stimulated"Arabidopsis）是植物體內(nèi)廣泛存在的小分子多肽，在植物生長(zhǎng)發(fā)育等過(guò)程中起重要作用，尤其在桑樹(shù)芽休眠調(diào)控方面具有潛在價(jià)值，但目前對(duì)桑樹(shù)中MaGASA基因在芽休眠過(guò)程中的功能了解有限。本研究旨在鑒定桑樹(shù)GASA基因家族成員，解析其在芽休眠中的作用。采用生物信息學(xué)方法鑒定GASA成員，并分析其理化性質(zhì)等特征，同時(shí)利用實(shí)時(shí)熒光定量PCR和轉(zhuǎn)錄組數(shù)據(jù)分析其表達(dá)模式。結(jié)果表明：在桑樹(shù)中存在12個(gè)MaGASA家族成員（MaGASA1～MaGASA12），分布在8條染色體上；成員間酸堿存在差異、親水性平均系數(shù)有正有負(fù)，但均含GASA保守結(jié)構(gòu)域；進(jìn)化樹(shù)顯示其與桃具有較近的親緣關(guān)系；共線性分析顯示，在MaGASA基因家族中存在3條共線性關(guān)系；MaGASAs啟動(dòng)子區(qū)域含響應(yīng)光照、低溫及多種激素的順式作用元件；MaGASA家族成員具有組織和品種特異性，其中MaGASA12在桑芽中的表達(dá)量顯著高于其他部位；轉(zhuǎn)錄組分析表明，在解除休眠過(guò)程中，MaGASA12表達(dá)量顯著降低，推測(cè)該基因是桑樹(shù)芽解除休眠的負(fù)調(diào)控基因；而MaGASA10表達(dá)量急劇上調(diào)，推測(cè)該基因可能是桑芽解除休眠的關(guān)鍵基因。通過(guò)實(shí)時(shí)熒光定量PCR技術(shù)對(duì)MaGASA10和MaGASA12在不同休眠階段的樣本表達(dá)量進(jìn)行檢測(cè)，結(jié)果進(jìn)一步驗(yàn)證了轉(zhuǎn)錄組分析的結(jié)論，表明MaGASA10和MaGASA12可能與桑屬冬芽休眠過(guò)程相關(guān)。研究結(jié)果為桑樹(shù)MaGASA生物功能的深入研究及遺傳育種提供理論參考。

關(guān)鍵詞：桑樹(shù)；GASA基因家族；休眠解除；生物信息學(xué)；轉(zhuǎn)錄組中圖分類號(hào)：S888.2""""""文獻(xiàn)標(biāo)志碼：A

Genome-wide"Identification"and"Expression"Pattern"Analysis"of"the"GASA"Gene"Family"in"Mulberry

HUANG"Jingjing1，"WANG"Shuchang2，3，"LIN"Peiqun2，3，"LOU"Dezhao2，3，"GENG"Tao2，3，"WU"Huazhou1，2，3*，"YU"Wengang1*

1."School"of"Tropical"Agriculture"and"Forestry，"Hainan"University，"Haikou，"Hainan"570228，"China;"2."National"Key"Laboratory"for"Tropical"Crop"Bio-breeding，"Sanya，"Hainan"572024，"China;"3."Environment"and"Plant"Protection"Institute，"Chinese"Academy"of"Tropical"Agricultural"Sciences"/"Key"Laboratory"of"Integrated"Pest"Management"of"Tropical"Crops，"Ministry"of"Agriculture"and"Rural"Affairs"/"Hainan"Key"Laboratory"of"Monitoring"and"Control"of"Tropical"Agricultural"Pests，"Haikou，"Hainan"571101，"China

Abstract："GASA"（gibberellic"acid-stimulated"Arabidopsis）"is"a"small"molecule"polypeptide"widely"presented"in"plants，"playing"an"important"role"in"plant"growth"and"development"processes."It"has"potential"value"especially"in"the"regulation"of"mulberry"bud"dormancy."However，"currently，"the"understanding"of"the"functions"of"MaGASA"genes"in"mulberry"during"the"bud"dormancy"process"is"limited."The"study"aimed"to"identify"the"members"of"the"mulberry"GASA"gene"family"and"analyze"the"roles"in"bud"dormancy."Bioinformatics"methods"were"used"to"identify"the"members"and"analyze"the"physicochemical"properties"and"other"characteristics."Meanwhile，"real-time"fluorescence"quantitative"PCR"and"transcriptome"data"were"utilized"to"analyze"the"expression"patterns."The"results"showed"that"there"were"12"MaGASA"family"members"（MaGASA1-MaGASA12）"in"mulberry，"which"were"distributed"on"eight"chromosomes."There"were"differences"in"acidity"and"alkalinity"among"the"members，"and"the"average"hydrophilicity"coefficients"were"positive"or"negative，"but"all"contained"the"GASA"conserved"domain."The"phylogenetic"tree"showed"that"they"had"a"close"genetic"relationship"with"peach."Collinearity"analysis"revealed"that"there"were"three"collinear"relationships"in"the"MaGASA"gene"family."The"promoter"regions"of"MaGASAs"contained"cis-acting"elements"that"respond"to"light，"low"temperature，"and"various"hormones."The"MaGASA"family"members"showed"tissue-variety"specific"expression."Among"them，"the"expression"level"of"MaGASA12"in"mulberry"buds"was"significantly"higher"than"that"in"other"parts."Transcriptome"analysis"indicated"that"the"expression"level"of"MaGASA12"decreased"significantly"during"the"dormancy"breaking"process，"while"the"expression"level"of"MaGASA10"increased"sharply"during"the"dormancy"breaking"process."The"expression"levels"of"MaGASA10"and"MaGASA12"at"different"dormancy"stages"were"detected"by"real-time"fluorescence"quantitative"PCR"technology，"and"the"results"further"verified"the"conclusions"of"the"transcriptome"analysis，"indicating"that"they"may"be"related"to"the"winter"bud"dormancy"process"of"Morus."The"research"results"would"provide"a"theoretical"reference"for"the"in-depth"study"of"the"biological"functions"of"MaGASA"in"mulberry"and"genetic"breeding.

Keywords："Morus"sp.;"GASA"gene"family;"dormancy"release;"bioinformatics;"transcriptome

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

GASA（gibberellic"acid-stimulated"Arabidopsis）基因家族最初在擬南芥（Arabidopsis"thaliana）中被發(fā)現(xiàn)，現(xiàn)已在多種植物中被鑒定，包括番茄（Lycopersicon"esculentum）、馬鈴薯（Solanum"tuberosum）、葡萄（Vitis"vinifera）和桃（Prunus"persica）等[1-4]。該家族編碼的富含半胱氨酸的小分子蛋白質(zhì)主要定位于細(xì)胞壁，在植物的生長(zhǎng)發(fā)育、激素信號(hào)轉(zhuǎn)導(dǎo)以及生物和非生物脅迫響應(yīng)等過(guò)程中發(fā)揮重要作用[5]。

GASA基因家族參與作物生長(zhǎng)和發(fā)育過(guò)程的調(diào)節(jié)，包括莖伸長(zhǎng)[6]、成花誘導(dǎo)[5]、種子萌發(fā)[7]、側(cè)根形成[8]和果實(shí)發(fā)育[9]。還與多種植物激素關(guān)系極為密切，尤其是與赤霉素（GA）、脫落酸（ABA）共同構(gòu)成了復(fù)雜的激素調(diào)控網(wǎng)絡(luò)[10]。例如，擬南芥中的GASA4和GASA6通過(guò)調(diào)節(jié)活性氧（ROS）水平和氧化還原狀態(tài)來(lái)促進(jìn)種子萌發(fā)和莖伸長(zhǎng)[6]。在葡萄和桃樹(shù)中，GASA基因家族成員對(duì)激素的響應(yīng)會(huì)因植物的發(fā)育階段以及組織類型的不同而表現(xiàn)出差異[4，"7]。

赤霉素是一類四環(huán)二萜類植物激素，在種子萌發(fā)、細(xì)胞伸長(zhǎng)和細(xì)胞分裂等多個(gè)植物生長(zhǎng)發(fā)育過(guò)程及逆境脅迫響應(yīng)中發(fā)揮重要功能[11-12]。外源GA4處理能夠加速青梅（Prunus"mume）休眠的解除，而在葡萄藤的芽休眠解除期間檢測(cè)到GA1的增加[13-14]。因此，GA參與果樹(shù)冬芽休眠解除，具有正向作用[13-14]。而GASA基因家族作為DELLA蛋白下游的重要靶基因，與赤霉素信號(hào)傳導(dǎo)路徑緊密相連，極有可能在植物的休眠起始、維持以及解除等過(guò)程中發(fā)揮著至關(guān)重要的功能[11]。

桑樹(shù)（Morus"sp.）不僅是傳統(tǒng)養(yǎng)蠶業(yè)的主要飼料來(lái)源，同時(shí)也是藥食同源的重要植物資源，具有極高的經(jīng)濟(jì)價(jià)值和生態(tài)價(jià)值。隨著“北桑南移”“立桑為業(yè)，多元發(fā)展”的產(chǎn)業(yè)結(jié)構(gòu)調(diào)整，果桑在鄉(xiāng)村振興和人們健康食品來(lái)源中發(fā)揮重要作用，但隨著全球溫度升高和產(chǎn)業(yè)南移，低溫不足引起的冬芽萌芽率低和不整齊成為果桑產(chǎn)業(yè)發(fā)展的主要障礙之一。前期研究發(fā)現(xiàn)，引種至熱帶地區(qū)的桑樹(shù)無(wú)法正常萌發(fā)，利用單氰胺涂抹休眠芽可以使其在同一時(shí)間萌芽并結(jié)果[15]。已有研究表明，桑樹(shù)的DAM基因、ABA合成與分解相關(guān)基因等具有調(diào)控桑樹(shù)休眠的功能[16-17]。然而，可能參與休眠調(diào)控的GASA基因在桑樹(shù)中仍未被鑒定和分析。這限制了我們?cè)谏?shù)品種改良中充分利用GASA基因的潛力——即通過(guò)調(diào)控其表達(dá)、挖掘功能等方式，來(lái)提升桑樹(shù)抗逆性、改善桑葉品質(zhì)及優(yōu)化休眠特性等重要農(nóng)藝性狀。因此，本研究在全基因組范圍內(nèi)對(duì)桑樹(shù)MaGASA家族成員進(jìn)行鑒定，對(duì)其基因和蛋白結(jié)構(gòu)、蛋白特性、系統(tǒng)進(jìn)化關(guān)系、啟動(dòng)子順式作用元件等進(jìn)行分析，并探究MaGASA基因家族在不同組織、不同品種以及單氰胺處理桑芽解除休眠過(guò)程中的表達(dá)模式，探討桑樹(shù)MaGASA基因家族可能的生物學(xué)功能，旨在為后續(xù)深入研究桑樹(shù)MaGASA基因家族在冬芽休眠中的功能提供參考，同時(shí)為桑樹(shù)分子育種提供候選基因及分子理論基礎(chǔ)。

1""材料與方法

1.1""桑樹(shù)GASA家族基因成員鑒定

基于NCBI官網(wǎng)（https：//www.ncbi.nlm.nih."gov/）獲取的桑樹(shù)基因組數(shù)據(jù)，使用hmmer軟件鑒定基因家族。在InterPro數(shù)據(jù)庫(kù)（https：//www."ebi.ac.uk/interpro/）直接搜索關(guān)鍵詞，得到相應(yīng)的pfam號(hào)（PF03106），下載hmm文件。用hmmsearch工具鑒定基因家族，提取結(jié)構(gòu)域序列后截取結(jié)構(gòu)域，使用muscle工具對(duì)結(jié)構(gòu)域進(jìn)行比對(duì)后建立物種特異性hmm模型，再次搜索。篩選hmmsearch結(jié)果中E-valuelt;0.001的記錄。獲得初步基因家族的基因列表，利用seqkit軟件得到對(duì)應(yīng)基因的蛋白序列，通過(guò)在線數(shù)據(jù)庫(kù)確認(rèn)結(jié)構(gòu)域是否存在，同時(shí)確認(rèn)結(jié)構(gòu)域長(zhǎng)度。最后共鑒定出12個(gè)桑樹(shù)GASA基因。

1.2""桑樹(shù)總RNA提取及基因全長(zhǎng)的獲得

取桑樹(shù)品種四季果72的莖、葉、果、芽不同部位，四季果72、大十、拐桑等12個(gè)桑樹(shù)品種的休眠芽，以及單氰胺處理0、72、120、168、216"h的無(wú)籽大10的休眠芽，每個(gè)樣品設(shè)置3個(gè)重復(fù)，液氮研磨后使用CTAB法分別提取RNA，利用HiScript"Ⅲ"RT"SuperMix"for"qPCR"（+gDNA"wiper）反轉(zhuǎn)錄得到cDNA。

1.3""桑樹(shù)GASA家族成員理化性質(zhì)分析

利用TBtools軟件對(duì)得到的蛋白序列文件進(jìn)行理化性質(zhì)分析，得到GASA基因家族中12個(gè)基因的氨基酸、理論等電點(diǎn)及分子量等信息。

1.4""桑樹(shù)GASA家族成員的染色體定位

從數(shù)據(jù)庫(kù)中得到GASA家族成員的染色體位置信息，利用MapChart軟件（https：//www.wur.nl/"en/show/Mapchart.htm）繪制GASA在染色體中的位置。

1.5""桑樹(shù)GASA蛋白進(jìn)化樹(shù)分析

采用鄰接法（neighbor-joining，"NJ），使用MEGA軟件構(gòu)建桑樹(shù)、擬南芥、楊樹(shù)、桃、葡萄GASA蛋白家族的系統(tǒng)進(jìn)化樹(shù)，檢驗(yàn)參數(shù)bootstrap重復(fù)為1000次，其他參數(shù)均為默認(rèn)值。

1.6""桑樹(shù)GASA基因啟動(dòng)子順式作用元件分析

利用NCBI網(wǎng)站獲得GASA家族各成員起始密碼子上游2000"bp的序列，將獲得的序列提交至Plant"Care網(wǎng)站（http：//bioinformatics.psb.ugent."be/webtools/plantcare/html/）進(jìn)行順式作用元件預(yù)測(cè)。

1.7""桑樹(shù)GASA家族結(jié)構(gòu)域及保守motif分析

利用MEME、GSDS"2.0（http：//gsds.gao-lab."org/）在線工具分析桑樹(shù)GASA家族結(jié)構(gòu)域及保守motif。

1.8""桑樹(shù)GASA蛋白理化性質(zhì)及結(jié)構(gòu)分析

分別利用SOPMA（https：//npsa.lyon.inserm.fr/"cgi-bin/npsa_automat.pl？page=/NPSA/npsa_sopma.html）及swiss-model在線工具（https：//swissmodel."expasy.org/）預(yù)測(cè)GASA蛋白的二級(jí)結(jié)構(gòu)和三級(jí)結(jié)構(gòu)。

1.9""轉(zhuǎn)錄組數(shù)據(jù)的獲得

試驗(yàn)采樣前1"d，對(duì)桑樹(shù)田進(jìn)行漫灌，直至土壤充分濕潤(rùn)，表面可見(jiàn)少量積水。選擇桑樹(shù)品種大十樹(shù)齡一致的12株植株作為試材，每株植株選取6根生長(zhǎng)狀況相近的枝條，先取各枝條上的休眠芽作為0"h（CK0h）的樣本，采后迅速置于液氮中保存；然后將每株植株的6根枝條平均分為2組，其中3根枝條使用2%單氰胺溶液均勻涂抹休眠芽，其余3根枝條使用等量的清水進(jìn)行涂抹。涂抹處理后，分別在72、120、216"h，從各處理組的枝條上采集休眠芽樣本。用清水處理分別為CK72h、CK120h、CK216h，單氰胺溶液處理分別為HC72h、HC120h、HC216h。樣本采集后，在液氮環(huán)境下迅速進(jìn)行研磨，樣品送至武漢邁維代謝生物科技股份有限公司進(jìn)行轉(zhuǎn)錄組檢測(cè)。

1.10""桑樹(shù)GASA基因的表達(dá)分析

根據(jù)GASA各基因序列，設(shè)計(jì)特異定量引物（表1），利用南京諾唯贊生物科技股份有限公司的定量試劑盒（ChamQ"Universal"SYBR"qPCR"Master"Mix），每個(gè)樣品設(shè)置3個(gè)重復(fù)，參照說(shuō)明書(shū)進(jìn)行GASA的特異擴(kuò)增，以桑樹(shù)的Actin基因作為內(nèi)參，使用熒光定量PCR儀（Thermo"Fisher"Scientific）擴(kuò)增。反應(yīng)程序：95"℃預(yù)變性15"min；95"℃變性10"s，60"℃退火30"s，72"℃延伸30"s，40個(gè)循環(huán)。采用2–??CT進(jìn)行數(shù)據(jù)分析，并使用IBM"SPSS"Statistics"27軟件進(jìn)行差異顯著性分析，利用origin"2021軟件作圖。

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

2.1""桑樹(shù)GASA家族成員理化性質(zhì)分析以及染色體定位

在InterPro網(wǎng)站下載hmm模型，使用hmmsearch工具鑒定基因家族，通過(guò)2次比對(duì)，篩選到桑樹(shù)12個(gè)GASA基因家族成員，使用TBtools軟件進(jìn)行理化性質(zhì)分析。根據(jù)chr1～chr14染色體從上到下的位置，將其命名為MaGASA1～"MaGASA12。從表2中可以看出，GASA基因家族成員間氨基酸數(shù)量與分子量存在明顯差異，氨基酸數(shù)從87～272不等，分子量相應(yīng)在9232.87～28"043.82"kD之間變化，暗示其功能多樣性。理論等電點(diǎn)范圍為5.02～10.10，MaGASA1等電點(diǎn)為5.02，呈酸性，MaGASA9等電點(diǎn)為10.10，呈堿性，其余成員的等電點(diǎn)分布在二者之間，顯示出成員酸堿性質(zhì)的差異。不穩(wěn)定指數(shù)高低有別，高指數(shù)成員MaGASA1（60.15）可能更易降解，低指數(shù)成員則相對(duì)穩(wěn)定，不穩(wěn)定指數(shù)較高的成員可能在體內(nèi)更容易降解，其半衰期較短，可能需要更頻繁地合成來(lái)維持其功能；而不穩(wěn)定指數(shù)低的成員相對(duì)更穩(wěn)定，可能在細(xì)胞內(nèi)長(zhǎng)期存在并持續(xù)發(fā)揮作用。脂肪族指數(shù)的差異意味著熱穩(wěn)定性和柔韌性不同。親水性平均系數(shù)反映了成員的親疏水特性，有正有負(fù)，MaGASA3的親水性平均系數(shù)為0.141，相對(duì)疏水，MaGASA2為–0.447，相對(duì)親水。親水性強(qiáng)的可能更多存在于細(xì)胞質(zhì)等水環(huán)境，疏水性稍強(qiáng)的可能與細(xì)胞膜等疏水結(jié)構(gòu)關(guān)聯(lián)更緊密。

通過(guò)NCBI檢索得到GASA基因家族成員在染色體中的起始位置和基因長(zhǎng)度，并利用Map"Chart軟件繪制染色體分布圖（圖1）。與蘋(píng)果[18]、大豆[19]及萵苣[20]類似，GASA基因在桑樹(shù)14條染色體上分布不均勻。發(fā)現(xiàn)GASA基因家族成員分別分布于2號(hào)、3號(hào)、5號(hào)、6號(hào)、7號(hào)、8號(hào)、10號(hào)、14號(hào)等8條染色體上；MaGASA1、MaGASA7、MaGASA3、MaGASA4位于染色體中位區(qū)，其中MaGASA3、MaGASA4可能屬于大片段復(fù)制基因?qū)?；MaGASA2、MaGASA8、MaGASA6、MaGASA9、MaGASA10、MaGASA12位于染色體上端，其中MaGASA9和MaGASA10基因不僅處于同一染色體上，且二者距離很近，這種緊密的毗鄰關(guān)系強(qiáng)烈暗示它們構(gòu)成了基因簇；MaGASA5、MaGASA11位于染色體下端。

2.2""不同物種GASA基因家族進(jìn)化樹(shù)分析以及基因復(fù)制

對(duì)擬南芥、桃、葡萄、毛果楊和桑樹(shù)共計(jì)77個(gè)GASA蛋白進(jìn)行進(jìn)化樹(shù)分析。根據(jù)序列相似性聚類，系統(tǒng)發(fā)育關(guān)系可分為3類：Ⅰ類（18個(gè)基因）、Ⅱ類（25個(gè)基因）、Ⅲ類（34個(gè)基因）。系統(tǒng)進(jìn)化樹(shù)顯示，桑樹(shù)與桃出現(xiàn)了相對(duì)頻繁的基因共支特征。具體表現(xiàn)為：MaGASA11和桃的PRUPE"3G161500處于同一分支，MaGASA9與PRUPE"3G279500在另一分支上匯聚，MaGASA3與PRUPE"1G365600也在另一分支上匯聚（圖2A）。鑒于這些基因在進(jìn)化樹(shù)上的緊密聚類，從進(jìn)化角度推測(cè)，這些處于同一分支的基因很可能具有相似的功能。與本研究所選的其他物種相比，桑樹(shù)與桃相關(guān)基因的共支情況表明，這可能是由于二者都是薔薇目植物，也暗示了二者在這些基因所代表的遺傳層面上擁有比本研究所選其他物種更近的親緣關(guān)系以及更為相似的親緣演變軌跡。這意味著在漫長(zhǎng)的進(jìn)化歷程中，這2個(gè)物種在這些基因相關(guān)的遺傳信息傳遞與演化上，可能經(jīng)歷了較為相似的選擇壓力和遺傳漂變過(guò)程。

基因復(fù)制驅(qū)動(dòng)著基因組和遺傳系統(tǒng)的進(jìn)化，并促成了基因家族中基因結(jié)構(gòu)和功能的多樣性。MaGASA基因的3個(gè)全基因組重復(fù)（WGD）對(duì)（MaGASA5/7、MaGASA8/10和MaGASA9/12）分布在5條染色體上（圖2B）。

通過(guò)創(chuàng)建多個(gè)物種共線性圖揭示MaGASA基因與擬南芥、毛果楊、桃和葡萄的進(jìn)化關(guān)系（圖2C）。共線性分析顯示，MaGASA基因的多個(gè)拷貝存在于其他物種中，其中在桑樹(shù)基因組中發(fā)現(xiàn)了6個(gè)拷貝，在擬南芥基因組中發(fā)現(xiàn)了7個(gè)拷貝，在毛果楊基因組中發(fā)現(xiàn)了11個(gè)拷貝，桃基因組中發(fā)現(xiàn)7個(gè)拷貝。這表明不同物種間MaGASA基因拷貝數(shù)的差異可能反映了它們?cè)谶M(jìn)化歷程中經(jīng)歷了不同的選擇壓力和基因復(fù)制事件。桑樹(shù)、擬南芥、毛果楊、桃以及葡萄中拷貝數(shù)的不同或許暗示著該基因在各物種中的功能分化和適應(yīng)性變化。特別是桑樹(shù)的6個(gè)拷貝中有3個(gè)在5個(gè)物種中均存在同源基因，這可能代表著這些拷貝在物種進(jìn)化過(guò)程中具有相對(duì)保守的功能和重要的生物學(xué)意義。

2.3""桑樹(shù)GASA基因結(jié)構(gòu)域及蛋白保守motif分析

保守motif的存在可以進(jìn)一步用于對(duì)基因家族進(jìn)行分類。在桑樹(shù)MaGASA基因結(jié)構(gòu)域及蛋白保守motif分析中，由于motif"1普遍存在于所有MaGASA基因中，推測(cè)motif"1所對(duì)應(yīng)的氨基酸序列可能執(zhí)行GASA基因家族最基礎(chǔ)和保守的功能；motif"2和motif"4在除MaGASA1之外的所有基因中發(fā)現(xiàn)；而motif"9僅在MaGASA2和MaGASA7中存在（圖3），說(shuō)明這2個(gè)基因在功能上可能存在緊密聯(lián)系，擁有一些區(qū)別于其他成員的特殊功能，或許與特定的生理過(guò)程或細(xì)胞環(huán)境響應(yīng)有關(guān)。含相同保守motif的蛋白，或形成相似結(jié)構(gòu)域，具備相似功能，參與類似生化反應(yīng)或信號(hào)通路[21]。例如，motif"6和motif"8在MaGASA1和MaGASA2中被發(fā)現(xiàn)，可能表明這2個(gè)基因在結(jié)構(gòu)和功能上具有相似性。根據(jù)motif的分布和差異，進(jìn)一步推測(cè)基因家族內(nèi)不同亞組的motif差異可能導(dǎo)致它們?cè)谏飳W(xué)功能上的多樣性。這與已有研究一致，已有研究表明，具有相同motif的基因在同一亞組內(nèi)可能執(zhí)行相似的功能，而不同亞組之間的motif差異可能使其具有不同的功能[22]。motif"3僅存在于MaGASA9，這暗示MaGASA9可能在進(jìn)化過(guò)程中獲得了獨(dú)特的功能，與其他家族成員有所差異。值得注意的是，4個(gè)MaGASA家族蛋白成員同時(shí)含有4個(gè)motif，這些motif可能共同構(gòu)成特定的蛋白結(jié)構(gòu)域。結(jié)構(gòu)域往往具有特定的空間結(jié)構(gòu)和功能，推測(cè)這些蛋白在結(jié)構(gòu)上具有相似性，可能在細(xì)胞內(nèi)參與類似的生化反應(yīng)或信號(hào)通路[23]。不同的motif組合可能決定了MaGASA蛋白與其他蛋白之間的相互作用特異性，不同motif組合的蛋白與另一組蛋白相互作用，從而構(gòu)建起復(fù)雜的細(xì)胞內(nèi)相互作用網(wǎng)絡(luò)，共同調(diào)控生物學(xué)過(guò)程。這類似于擬南芥中fruitfull（ful）亞家族的MADS結(jié)構(gòu)域轉(zhuǎn)錄因子通過(guò)一個(gè)關(guān)鍵的氨基酸motif決定其與agamous和sepallata蛋白的相互作用特異性[24]。

2.4""桑樹(shù)GASA基因啟動(dòng)子順式作用元件分析

基因的轉(zhuǎn)錄表達(dá)高度依賴于該基因的上游啟動(dòng)子區(qū)，選取MaGASA上游2000"bp序列作為啟動(dòng)子，分析其順式作用元件。結(jié)果顯示，12個(gè)MaGASA基因啟動(dòng)子區(qū)包含多種順式作用元件（表3）。選取與桑樹(shù)休眠現(xiàn)象相關(guān)的激素類順式作用元件作為研究對(duì)象，發(fā)現(xiàn)其中含有脫落酸響應(yīng)元件（ABRE、ABRE2、ABRE3a、ABRE4）、赤霉素響應(yīng)元件（GARE-motif、P-box）、乙烯響應(yīng)元件（ERE），以及其他參與植物激素信號(hào)轉(zhuǎn)導(dǎo)的一類順式作用元件（TGA-box、TGACG-motif）。此外，還發(fā)現(xiàn)了與桑樹(shù)休眠相關(guān)的低溫響應(yīng)元件（LTR）、光響應(yīng)元件（TCT-motif），這和其他物種中GASA基因啟動(dòng)子中順式作用元件基本一致[25]。GASA基因啟動(dòng)子上的順式作用元件分析表明，桑樹(shù)可能經(jīng)過(guò)低溫春化后，通過(guò)調(diào)控激素相關(guān)作用元件，特別是赤霉素響應(yīng)元件，使桑芽萌動(dòng)伸長(zhǎng)。通過(guò)進(jìn)一步研究，如使用CRISPR/Cas基因組編輯系統(tǒng)和過(guò)表達(dá)等手段，可以闡明這些順式調(diào)節(jié)元件的作用[26]。

2.5""MaGASA蛋白的二級(jí)和三級(jí)結(jié)構(gòu)預(yù)測(cè)

蛋白質(zhì)的三級(jí)結(jié)構(gòu)有助于蛋白質(zhì)功能的準(zhǔn)確表征[27]?；趯?duì)MaGASA基因家族的二級(jí)空間結(jié)構(gòu)的預(yù)測(cè)（表4），發(fā)現(xiàn)MaGASA基因家族的蛋白質(zhì)結(jié)構(gòu)主要由無(wú)規(guī)卷曲和α-螺旋組成，但也存在延伸鏈。在葡萄[4]、蘋(píng)果[18]、白菜[25]和谷子[28]的GASA蛋白中也發(fā)現(xiàn)了類似的結(jié)構(gòu)特征，并且由于α-螺旋的存在，GASA基因的所有蛋白質(zhì)均具有柔性結(jié)構(gòu)，MaGASA的三級(jí)結(jié)構(gòu)如圖4所示。

2.6""GASA基因家族在桑樹(shù)不同部位的表達(dá)分析

從圖5中可以看出，MaGASA基因家族成員的組織表達(dá)模式存在差異。MaGASA1、MaGASA2、MaGASA4、MaGASA6在葉片中表達(dá)量最高，其中MaGASA1、MaGASA2、MaGASA6在葉片的相對(duì)表達(dá)量顯著高于其他部位，MaGASA4顯著高于除芽外的其他部位，值得注意的是，這幾個(gè)基因均含有光響應(yīng)順式作用元件，這暗示它們可能在桑樹(shù)葉片

的發(fā)育和光合作用中發(fā)揮重要功能。MaGASA3、MaGASA5、MaGASA7、MaGASA8、MaGASA10在莖中的表達(dá)量最高，其中MaGASA3、MaGASA7、MaGASA8、MaGASA10的相對(duì)表達(dá)量顯著高于其他部位，MaGASA5顯著高于除葉片外的其他部位。MaGASA11與MaGASA12在芽中的相對(duì)表達(dá)量顯著高于其他部位，結(jié)合轉(zhuǎn)錄組結(jié)果表明，MaGASA12可能在調(diào)控芽活動(dòng)狀態(tài)中發(fā)揮重要作用，也暗示MaGASA11與MaGASA12在芽發(fā)育或生理活動(dòng)中具有重要意義。MaGASA基因家族成員在桑樹(shù)不同組織部位呈現(xiàn)顯著的差異表達(dá)模式，這種組織特異性表達(dá)為深入了解MaGASA基因家族在桑樹(shù)不同器官生長(zhǎng)發(fā)育中的功能提供線索。

2.7""GASA基因家族在桑樹(shù)不同品種芽中的表達(dá)分析

通過(guò)對(duì)四季果72、大十、拐桑等12個(gè)具有代表性的桑樹(shù)品種展開(kāi)GASA基因家族成員的表達(dá)量分析，結(jié)果顯示，MaGASA基因家族成員在各品種中呈現(xiàn)出豐富多樣的表達(dá)模式（圖6）。MaGASA1在強(qiáng)桑5號(hào)中的表達(dá)量顯著高于其他品種；MaGASA2在強(qiáng)桑5號(hào)和四季果72中的表達(dá)量顯著高于其他品種；MaGASA3在冀桑和紅瑪瑙中的表達(dá)量顯著高于其他品種；MaGASA4在黃金桑中的表達(dá)量最高，顯著高于其他品種；MaGASA5在川桑和紅瑪瑙中的表達(dá)量最高，顯著高于其他品種；MaGASA6在黑珍珠中的基因表達(dá)量顯著高于其他品種；MaGASA7在川桑中的表達(dá)量顯著高于除大十外的其他品種；MaGASA8在大十和拐桑中的表達(dá)量最高；MaGASA9在紅瑪瑙中的表達(dá)量最高；MaGASA10在黑珍珠和拐桑中的表達(dá)量顯著高于其他品種，其他品種間則無(wú)顯著差異；MaGASA11在冀桑中的表達(dá)量顯著高于其他品種，除冀桑外，其他品種之間無(wú)顯著差異；MaGASA12在川桑中的表達(dá)量顯著高于其他品種。

2.8""GASA基因家族在桑樹(shù)休眠解除中的表達(dá)分析

GASA基因家族各成員轉(zhuǎn)錄組熱圖表明，在休眠芽經(jīng)單氰胺處理不同時(shí)間后，其表達(dá)量發(fā)生了明顯變化（圖7）。與對(duì)照（CK）相比，經(jīng)單氰胺處理（HC）不同時(shí)間后其整體表達(dá)量呈現(xiàn)下調(diào)趨勢(shì)，MaGASA12和MaGASA10明顯響應(yīng)休眠信號(hào)的解除。通過(guò)RT-qPCR分析（圖8），單氰胺處

理后MaGASA12的表達(dá)量明顯下調(diào)，這種變化模式與番茄成熟過(guò)程中SiGASA1表達(dá)量的下調(diào)相似[29]，表明該基因可能使桑芽維持休眠狀態(tài)。但MaGASA10的表達(dá)量在經(jīng)過(guò)單氰胺處理后急劇上調(diào)，對(duì)休眠解除信號(hào)的響應(yīng)最為強(qiáng)烈，這與梅花的PmGASA16在花芽休眠解除后的高表達(dá)情況一致[30]，這充分表明MaGASA10極有可能是桑樹(shù)休眠解除過(guò)程的重要基因。

3""討論"

赤霉素主要通過(guò)降解DELLA蛋白，釋放PIFs轉(zhuǎn)錄因子與GASA基因家族成員的啟動(dòng)子位點(diǎn)相結(jié)合，從而推進(jìn)細(xì)胞伸長(zhǎng)進(jìn)程[11-12]。GASA蛋白在植物生長(zhǎng)發(fā)育、逆境脅迫、激素信號(hào)轉(zhuǎn)導(dǎo)等方面起著非常重要的作用[12，"25，"28]。

本研究在桑樹(shù)基因組中鑒定了12個(gè)GASA基因，并對(duì)這些基因的生物信息學(xué)和表達(dá)模式進(jìn)行了分析。與其他物種相比，可能由于桑樹(shù)與桃均為薔薇目植物，具有更近的親緣關(guān)系，出現(xiàn)了相對(duì)頻繁的基因共支特征。部分學(xué)者認(rèn)為果樹(shù)休眠受低溫調(diào)控，需要滿足一定的冷量才可解除其休眠狀態(tài)，在解除休眠過(guò)程中激素發(fā)揮著重要作用[31]。同時(shí)鈣網(wǎng)蛋白被認(rèn)為是GA信號(hào)通路中的關(guān)鍵蛋白，因?yàn)槠湓诩に匦盘?hào)通路中具有重要意義，而激素信號(hào)通路又可能參與山毛櫸種子休眠的解除及其萌發(fā)[32]?；趯?duì)激素信號(hào)基因表達(dá)模式的研究，其被認(rèn)為是休眠調(diào)控的關(guān)鍵因素[33]。在桑樹(shù)的MaGASA啟動(dòng)子區(qū)域篩選出與激素相關(guān)以及低溫和光響應(yīng)順式作用元件，其中赤霉素相關(guān)順式作用元件同樣在其啟動(dòng)子區(qū)域中被發(fā)現(xiàn)，這和其他物種GASA基因啟動(dòng)子的順式作用元件基本一致[25]。通過(guò)研究GASA基因啟動(dòng)子的順式作用元件表明，桑樹(shù)可能經(jīng)過(guò)低溫春化后，通過(guò)調(diào)控激素相關(guān)作用元件，特別是赤霉素響應(yīng)元件，使芽萌動(dòng)伸長(zhǎng)。

單氰胺在解除果樹(shù)休眠中被廣泛應(yīng)用，在打破桑樹(shù)休眠上同樣具有良好效果。其中MaGASA10表達(dá)變化劇烈，暗示其在單氰胺誘導(dǎo)的桑芽休眠解除中起關(guān)鍵作用，PmGASA16在梅花花芽休眠解除后具有較強(qiáng)的誘導(dǎo)作用，表明PmGASA16在梅花花芽休眠解除中可能具有促進(jìn)作用[30]。這或許表明單氰胺激活了這些基因，推動(dòng)休眠解除相關(guān)生理進(jìn)程。MaGASA12在單氰胺處理后顯著下調(diào)，伴隨休眠芽的萌動(dòng)其表達(dá)量持續(xù)下降，這表明它極有可能是單氰胺打破桑樹(shù)休眠的關(guān)鍵調(diào)控基因；在ZHANG等[34]的研究中，AtGASA4促進(jìn)擬南芥開(kāi)花，而AtGASA5延遲或抑制擬南芥開(kāi)花的情況相符。表明GASA基因家族成員在休眠解除過(guò)程中的引導(dǎo)通路較為復(fù)雜，而非單一途徑，可能引發(fā)多種信號(hào)通路互作，從而調(diào)整基因表達(dá)網(wǎng)絡(luò)，啟動(dòng)休眠解除關(guān)鍵生理過(guò)程。

4""結(jié)論

本研究系統(tǒng)鑒定了桑樹(shù)MaGASA基因家族成員，分析了其理化性質(zhì)、系統(tǒng)發(fā)育、基因復(fù)制和共線性，并對(duì)不同組織、品種和解除休眠過(guò)程的表達(dá)模式等進(jìn)行了分析。結(jié)果顯示，12個(gè)MaGASA基因成員均包含保守的GASA結(jié)構(gòu)域，大多數(shù)基因表現(xiàn)出相似的理化特性，且在系統(tǒng)發(fā)育分析中分為3個(gè)類群，表明該基因家族在桑樹(shù)中具有較高的多樣性?；驈?fù)制和共線性分析揭示了MaGASA基因在基因組中的擴(kuò)增，可能通過(guò)基因重復(fù)事件推動(dòng)了基因家族的多樣性。此外，組織特異性表達(dá)模式的分析表明，MaGASA12在芽中表達(dá)量顯著高于其他部位，提示其在芽發(fā)育過(guò)程中可能具有重要作用。同時(shí)，結(jié)合轉(zhuǎn)錄組數(shù)據(jù)，顯示該基因的表達(dá)可能與桑芽休眠有關(guān)。因此，今后將考慮通過(guò)敲除或沉默該基因以探究其作用。此外，MaGASA10在單氰胺處理后的表達(dá)量急劇上調(diào)，該基因可能是促進(jìn)桑樹(shù)萌芽的關(guān)鍵基因之一。未來(lái)的研究可以通過(guò)過(guò)表達(dá)MaGASA10來(lái)驗(yàn)證其功能，以進(jìn)一步闡明該基因在休眠解除中的作用。

綜上所述，本研究為桑樹(shù)MaGASA基因家族的功能研究提供了系統(tǒng)的生物信息學(xué)分析，揭示了基因家族成員在基因組中的擴(kuò)增、進(jìn)化以及組織特異性表達(dá)的特點(diǎn)，并重點(diǎn)關(guān)注了MaGASA10和MaGASA12基因表達(dá)量的變化，為今后的研究提供研究思路。

參考文獻(xiàn)

• ZHANG"S"C，"WANG"X"J."One"new"kind"of"phytohormonal"signaling"integrator："up-and-coming"GASA"family"genes[J]."Plant"Signaling"amp;"Behavior，"2017，"12（2）："e1226453.
• SHI"L，"GAST"R"T，"GOPALRAJ"M，"OLSZEWSKI"N"E."Characterization"of"a"shoot-specific，"GA3"and"ABA-regulated"gene"from"tomato[J]."The"Plant"Journal，"1992，"2（2）："153-159.
• NAHIR?AK"V，"RIVAROLA"M，"GONZALEZ"DE"URRETA"M，"PANIEGO"N，"HOPP"H"E，"ALMASIA"N"I，"VAZQUEZ-"ROVERE"C."Genome-wide"analysis"of"the"Snakin/GASA"gene"family"in"Solanum"tuberosum"cv."Kennebec[J]."American"Journal"of"Potato"Research，"2016，"93："172-188.
• AHMAD"B，"YAO"J，"ZHANG"S，"LI"X，"ZHANG"X，"YADAV"V，"WANG"X."Genome-wide"characterization"and"expression"profiling"of"GASA"genes"during"different"stages"of"seed"development"in"grapevine"（Vitis"vinifera"L.）"predict"their"involvement"in"seed"development[J]."International"Journal"of"Molecular"Sciences，"2020，"21（3）："1088.
• QU"J，"KANG"S"G，"HAH"C，"JANG"J"C."Molecular"and"cellular"characterization"of"GA-stimulated"transcripts"GASA4"and"GASA6"in"Arabidopsis"thaliana[J]."Plant"Science，"2016，"246："1-10.
• BEN-NISSAN"G，"LEE"J"Y，"BOROHOV"A，"WEISS"D."GIP，"a"Petunia"hybrida"GA-induced"cysteine-rich"protein："a"possible"role"in"shoot"elongation"and"transition"to"flowering[J]."The"Plant"Journal，"2004，"37（2）："229-238.
• ROXRUD"I，"LID"S"E，"FLETCHER"J"C，"SCHMIDT"E"D，"OPSAHL-SORTEBERG"H"G."GASA4，"one"of"the"14-"member"Arabidopsis"GASA"family"of"small"polypeptides，"regulates"flowering"and"seed"development[J]."Plant"and"Cell"Physiology，"2007，"48（3）："471-483.
• ZIMMERMANN"R，"SAKAI"H，"HOCHHOLDINGER"F."The"gibberellic"acid"stimulated-like"gene"family"in"maize"and"its"role"in"lateral"root"development[J]."Plant"Physiology，"2010，"152（1）："356-365.
• MOYANO-CANETE"E，"BELLIDO"M"L，"GARCíA-"CAPARRóS"N，"MEDINA-PUCHE"L，"AMIL-RUIZ"F，"GONZáLEZ-REYES"J"A，"CABALLERO"J"L，"MUNOZ-"BLANCO"J，"BLANCO-PORTALES"R."FaGAST2，"a"strawberry"ripening-related"gene，"acts"together"with"FaGAST1"to"determine"cell"size"of"the"fruit"receptacle[J]."Plant"and"Cell"Physiology，"2013，"54（2）："218-236.
• MUHAMMAD"I，"LI"W"Q，"JING"X"Q，"ZHOU"M"R，"SHALMANI"A，"ALI"M，"WEI"X"Y，"SHARIF"R，"LIU"W"T，"CHEN"K"M."A"systematic"in"silico"prediction"of"gibberellic"acid"stimulated"GASA"family"members："a"novel"small"peptide"contributes"to"floral"architecture"and"transcriptomic"changes"induced"by"external"stimuli"in"rice[J]."Journal"of"Plant"Physiology，"2019，"234："117-132.
• 令狐前前，"徐剛，"李雪，"陳雨倩."麻瘋樹(shù)JcGASA6基因克隆及其原核表達(dá)分析[J]."植物生理學(xué)報(bào)，"2021，"57（12）："2346-2356.

LINGHU"Q"Q，"XU"G，"LI"X，"CHEN"Y"Q."Cloning"and"prokaryotic"expression"analysis"of"JcGASA6"gene"in"Jatropha"curcas[J]."Plant"Physiology"Journal，"2021，"57（12）："2346-2356."（in"Chinese）

• HAN"S，"JIAO"Z，"NIU"M"X，"YU"X，"HUANG"M，"LIU"C，"WANG"H"L，"ZHOU"Y，"MAO"W，"WANG"X，"YIN"W，"XIA"X."Genome-wide"comprehensive"analysis"of"the"GASA"gene"family"in"Populus[J]."International"Journal"of"Molecular"Sciences，"2021，"22（22）："12336.
• ZHENG"C，"ACHEAMPONG"A"K，"SHI"Z，"MUGZECH"A，"HALALY-BASHA"T，"SHAYA"F，"SUN"Y，"COLOVA"V，"MOSQUNA"A，"OPHIR"R."Abscisic"acid"catabolism"enhances"dormancy"release"of"grapevine"buds[J]."Plant，"Cell"amp;"Environment，"2018，"41（10）："2490-2503.
• ZHUANG"W"B，"GAO"Z"H，"WANG"L"J，"ZHONG"W"J，"NI"Z"J，"ZHANG"Z."Comparative"proteomic"and"transcriptomic"approaches"to"address"the"active"role"of"GA4"in"Japanese"apricot"flower"bud"dormancy"release[J]."Journal"of"Experimental"Botany，"2013，"64（16）："4953-4966.
• 武華周，"王樹(shù)昌，"耿濤，"盧芙萍，"婁德釗，"韋顯紅."一種解除熱帶地區(qū)果桑冬芽休眠的高效栽培方法及其破眠劑："CN202111264326.1[P]."2021-10-28.

WU"H"Z，"WANG"S"C，"GENG"T，"LU"F"P，"LOU"D"Z，"WEI"X"H."An"efficient"cultivation"method"for"breaking"winter"bud"dormancy"of"fruit"mulberry"in"tropical"regions"and"its"dormancy-breaking"agent："CN202111264326.1[P]."2021-10-28."（in"Chinese）

• ZHU"H，"CHEN"P"Y，"ZHONG"S"L，"DARDICK"C，"CALLAHAN"A，"AN"Y"Q，"KNOCKER"S"V，"YANG"Y"Z，"ZHONG"G"Y，"ABBOTT"A，"LIU"Z"R."Thermal-responsive"genetic"and"epigenetic"regulation"of"DAM"cluster"controlling"dormancy"and"chilling"requirement"in"peach"floral"buds[J]."Horticulture"Research，"2023（1）："1444-1457.
• ZHENG"C"L，"HALALY"T，"ACHEAMPONG"A"K，"TAKEBAYASHI"Y，"JIKUMARU"Y，"KAMIYA"Y，"OR"E."Abscisic"acid"（ABA）"regulates"grape"bud"dormancy，nbsp;and"dormancy"release"stimuli"may"act"through"modification"of"ABA"metabolism[J]."Journal"of"Experimental"Botany，"2015，"66（5）："1527-1542.
• FAN"S，"ZHANG"D，"ZHANG"L，"GAO"C，"XIN"M，"TAHIR"M"M，"LI"Y，"MA"J，"HAN"M."Comprehensive"analysis"of"GASA"family"members"in"the"Malus"domestica"genome："identification，"characterization，"and"their"expressions"in"response"to"apple"flower"induction[J]."BMC"Genomics，"2017，"18："1-19.
• AHMAD"M"Z，"SANA"A，"JAMIL"A，"NASIR"J"A，"AHMED"S，"HAMEED"M"U，"ABDULLAH."A"genome-wide"approach"to"the"comprehensive"analysis"of"GASA"gene"family"in"Glycine"max[J]."Plant"Molecular"Biology，"2019，"100："607-620.
• SUN"H"L，"JIN"S"Y，"JUNG"W"J，"KIM"D"Y，"SEO"Y"W."Genome-wide"identification"and"characterization"of"the"lettuce"GASA"family"in"response"to"abiotic"stresses[J]."BMC"Plant"Biology，"2023，"23（1）："106.
• HONG"Z"Z，"WANG"X"R，"YANG"A"Y，"YAN"G"C，"HE"Y，"ZHU"Z"J，"XU"Y"M."Tendril"morphogenesis"is"regulated"by"a"CsaTEN-"CsaUFO"module"in"cucumber[J]."New"Phytologist，"2023，"239："364-373.
• HU"M"Q，"LI"Z"Q，"LIN"X"J，"TANG"B"Q，"XING"M，"ZHU"H"B."Comparative"analysis"of"the"LEA"gene"family"in"seven"Ipomoea"species，"focuses"on"sweet"potato"（Ipomoea"batatas"L.）[J]."BMC"Plant"Biology，"2024，"24（1）："1256.
• 雷彪，"曲瑞芳，"任超，"盧成達(dá)，"李紅英，"韓彥卿，"馬芳芳，"韓淵懷，"邢國(guó)芳."谷子LBD基因家族的鑒定及其對(duì)逆境脅迫的響應(yīng)[J]."植物生理學(xué)報(bào)，"2023，"59（3）："527-542.

LEI"B，"QU"R"F，"REN"C，"LU"C"D，"LI"H"Y，"HAN"Y"Q，"MA"F"F，"HAN"Y"H，"XING"G"F."Identification"of"the"LBD"gene"family"in"foxtail"millet"and"its"response"to"abiotic"stresses[J]."Plant"Physiology"Journal，"2023，"59（3）："527-542."（in"Chinese）

• KAI"T，"MIGUEL"C"M，"MARTIJN"F，"LAI"X，"IRISE"Z，"JIANG"X，"MARK"M，"STEFAN"B，"STUART"J，"MAX"N."Uncoupling"FRUITFULL’s"functions"through"modification"of"a"protein"motif"identified"by"co-ortholog"analysis[J]."Nucleic"Acids"Research，"2024，"52（21）："13290-13304.
• SUN"B，"ZHAO"X，"GAO"J，"LI"J，"XIN"Y，"ZHAO"Y，"LIU"Z，"FENG"H，"TAN"C."Genome-wide"identification"and"expression"analysis"of"the"GASA"gene"family"in"Chinese"cabbage"（Brassica"rapa"L."ssp."pekinensis）[J]."BMC"Genomics，"2023，"24（1）："668.
• ABDULLAH，"FARAJI"S，"MEHMOOD"F，"MALIK"H"M"T，"AHMED"I，"HEIDARI"P，"POCZAI"P."The"GASA"gene"family"in"cacao"（Theobroma"cacao，"Malvaceae）："genome"wide"identification"and"expression"analysis[J]."Agronomy，"2021，"11（7）："1425.
• HAN"S，"JIAO"Z"Y，"NIU"M"X，"YU"X，"HUANG"M"B，"LIU"C，"WANG"H"L，"ZHOU"Y"Y，"MAO"W，"WANG"X"F，"YIN"W"L，"XIA"X"L."Genome-wide"comprehensive"analysis"of"the"GASA"gene"family"in"Populus[J]."International"Journal"of"Molecular"Sciences，"2021，"22（22）："12336.
• 王東勝，"張小朵，"尚炫辰，"王振華，"劉維仲."谷子GASA基因家族鑒定及其對(duì)非生物脅迫的響應(yīng)[J]."植物生理學(xué)報(bào)，"2024，"60（5）："859-869.

WANG"D"S，"ZHANG"X"D，"SHANG"X"C，"WANG"Z"H，"LIU"W"Z."Identification"of"the"GASA"gene"family"in"foxtail"millet"and"its"response"to"abiotic"stresses[J]."Plant"Physiology"Journal，"2024，"60（5）："859-869."（in"Chinese）

• SU"D，"LIU"K"D，"YU"Z"S，"LI"Y，"ZHANG"Y"X，"ZHU"Y"Q，"WU"Y，"HE"H"Y，"ZENG"X"D，"CHEN"H"L，"GRIERSON"D，"DENG"H，"LIU"M"C."Genome-wide"characterization"of"the"tomato"GASA"family"identifies"SlGASA1"as"a"repressor"of"fruit"ripening[J]."Horticulture"Research，"2023（1）："135-145.".
• ZHANG"M，"CHENG"W，"WANG"J，"CHENG"T，"ZHANG"Q."Genome-wide"identification，"evolution，"and"expression"analysis"of"GASA"gene"family"in"Prunus"mume[J].nbsp;International"Journal"of"Molecular"Sciences，"2022，"23（18）："10923.
• HEIDE"O"M，"PRESTRUD"A"K."Low"temperature，"but"not"photoperiod，"controls"growth"cessation"and"dormancy"induction"and"release"in"apple"and"pear[J]."Tree"Physiology，"2005，"25（1）："109-114.
• PAW?OWSKI"T"A."Proteomics"of"European"beech"（Fagus"sylvatica"L.）"seed"dormancy"breaking："influence"of"abscisic"and"gibberellic"acids[J]."Proteomics，"2007，"7（13）："2246-2257.
• RINNE"P"L，"WELLING"A，"VAHALA"J，"RIPEL"L，"RUONALA"R，"KANGASJ?RVI"J，"VAN"DER"SCHOOT"C."Chilling"of"dormant"buds"hyperinduces"FLOWERING"LOCUS"T"and"recruits"GA-inducible"1，3-β-glucanases"to"reopen"signal"conduits"and"release"dormancy"in"Populus[J]."The"Plant"Cell，"2011，"23（1）："130-146.
• ZHANG"S"C，"WANG"X"J."Overexpression"of"GASA5"increases"the"sensitivity"of"Arabidopsis"to"heat"stress[J]."Journal"of"Plant"Physiology，"2011，"168（17）："2093-2101.