藜麥ZF-HD轉(zhuǎn)錄因子的全基因組鑒定及其對鹽脅迫的響應(yīng)分析

1、PAGE 16 -藜麥ZF-HD轉(zhuǎn)錄因子的全基因組鑒定及其對鹽脅迫的響應(yīng)分析摘要：以藜麥的高質(zhì)量基因組為參考，對藜麥鋅指同源異型結(jié)構(gòu)域（ZF-HD）轉(zhuǎn)錄因子基因CqZF-HD進行全基因組鑒定，并利用生物信息學(xué)方法對其編碼的蛋白質(zhì)理化性質(zhì)、二級結(jié)構(gòu)、亞細胞定位、保守結(jié)構(gòu)域以及基因結(jié)構(gòu)、系統(tǒng)進化關(guān)系進行了分析。同時，利用前期轉(zhuǎn)錄組測序結(jié)果分析了CqZF-HD家族基因在鹽脅迫下的表達模式。結(jié)果表明，在藜麥中共鑒定出21個ZF-HD轉(zhuǎn)錄因子，大部分定位在細胞質(zhì)，二級結(jié)構(gòu)以無規(guī)則卷曲為主，氨基酸序列長度68284aa，相對分子質(zhì)量7.5410330.20103，理論等電點6.019.54;每個ZF-HD

2、蛋白都含有ZF或HD保守結(jié)構(gòu)域，系統(tǒng)進化樹將其劃分為5個亞家族。鹽處理下表達分析結(jié)果顯示，11個CqZF-HD基因響應(yīng)藜麥鹽脅迫，其中8個基因呈上調(diào)表達，3個基因呈下調(diào)表達，說明ZF-HD家族基因在參與藜麥鹽脅迫中發(fā)揮著不同的作用。關(guān)鍵詞：藜麥;鋅指同源異型結(jié)構(gòu)域（ZF-HD）;全基因組;鹽脅迫;基因表達中圖分類號：S512.9文獻標識碼：A文章編號：1000-4440（2022）02-0304-08Genome-wideidentificationofZF-HDtranscriptionfactorsandexpressionanalysisofresponsetosaltstressinq

3、uinoaSHIPi-biao，WANGDe-ling，JIANGRun-zhi，LIBin，YANJun，F(xiàn)ANGDi，YANKai-xuan，GUMin-fengAbstract：Zincfinger-homeodomain（ZF-HD）proteinsareplant-specifictranscriptionfactorsthatplayimportantrolesinplantgrowth，developmentandvariousstressresponses.Inthisstudy，theCqZF-HDgenewasidentifiedatgenome-widelevelwith

4、thehigh-qualityquinoagenomeasareference.Thephysicochemicalproperties，secondarystructure，subcellularlocalizationandconserveddomainsofCqZF-HDs，phylogeneticrelationshipandgenestructurewereanalyzedbybioinformaticsmethod.Meanwhile，theexpressionpatternofCqZF-HDfamilygenesundersaltstresswasalsoanalyzedbypr

5、eliminarytranscriptomesequencingresults.Theresultsshowedthat21ZF-HDtranscriptionfactorswereidentifiedinquinoa，mostofwhichwerelocatedinthecytoplasm，thesecondarystructurewasdominatedbyrandomcoils，aminoacidsequencelengthwas68-284aa，relativemolecularweightwas7.54103-30.20103，andtheoreticalisoelectricpoi

6、ntwas6.01-9.54.AlloftheZF-HDsweredividedintofivesubfamiliesbyphylogenetictreeanalysis，andeachmembercontainedconserveddomainsZForHD.Theexpressionanalysisundersalttreatmentshowedthat11CqZF-HDsrespondedtosaltstressinquinoa，eightgeneswereup-regulated，andthreegenesweredown-regulated.Inconclusion，ZF-HDfam

7、ilygenesplaydifferentrolesinresponsetosaltstressinquinoa.Keywords：quinoa;zincfinger-homeodomain（ZF-HD）;genome-wide;saltstress;geneexpression鋅指同源異型結(jié)構(gòu)域（ZF-HD）蛋白是植物所特有的一類轉(zhuǎn)錄因子，在植物生長發(fā)育和各種脅迫反應(yīng)中發(fā)揮著重要作用1-2。ZF-HD轉(zhuǎn)錄因子的氨基酸序列中具有保守的HD結(jié)構(gòu)域（Homeodomain）和C2H2型ZF結(jié)構(gòu)域（C2H2-typezincfingerdomain）3-4。HD結(jié)構(gòu)域是一種特征明確的DNA結(jié)合結(jié)構(gòu)域

8、，通過與DNA結(jié)合調(diào)控靶基因的表達水平，從而參與植物的生長發(fā)育5。ZF結(jié)構(gòu)域廣泛存在于調(diào)控蛋白中，并與DNA結(jié)合和蛋白質(zhì)互作，在植物抵御脅迫反應(yīng)中起著重要作用6-7。ZF-HD蛋白通常以一種核心序列ATTA與特異性DNA序列結(jié)合，形成同型二聚體或異型二聚體1。ZF-HD家族可分為ZHD和MIF（Minizincfinger）2個亞家族，其中MIF基因編碼ZF-HD蛋白的ZF結(jié)構(gòu)域，并不編碼HD結(jié)構(gòu)域3。在黃頂菊中首次發(fā)現(xiàn)ZF-HD蛋白，可能是編碼C4磷酸烯醇式丙酮酸羧化酶基因的調(diào)控因子8。隨后，在擬南芥1、水稻9、番茄10、大白菜4、棉花11和苦蕎12等植物中也發(fā)現(xiàn)了ZF-HD基因，并對其進行

9、了基因家族鑒定及功能分析。擬南芥中有17個ZF-HD成員，它們作為轉(zhuǎn)錄因子，具有獨特的生理特性，在花的發(fā)育中起著非常重要的作用3。AtZHD1能夠與ERD1（EARLYRESPONSETODEHYDRATIONSTRESS1）啟動子特異性結(jié)合，并且其表達受干旱、鹽和脫落酸（ABA）處理誘導(dǎo)，NAC和AtZHD1的過表達增強了擬南芥的耐旱性13。AtHB33受ARF2（Auxinresponsefactor2）負調(diào)控，是種子萌發(fā)和主根生長所必需的14。OsZHD1和OsZHD2通過調(diào)控鱗狀細胞的數(shù)量和排列方式而引起水稻葉片卷曲15;OsZHD2的過表達不僅能促進根的生長，還能增加水稻產(chǎn)量16。干

10、旱、NaCl和低溫處理對番茄SlZHD18基因具有明顯的誘導(dǎo)作用，而熱脅迫對其沒有明顯誘導(dǎo)作用10。在響應(yīng)干旱脅迫和高鹽脅迫時，葡萄VvZHD基因呈上調(diào)表達17。因此，ZF-HD基因不僅在植物的各個生長發(fā)育階段發(fā)揮著一定的作用，還參與了植物對鹽分等逆境脅迫的響應(yīng)。藜麥是一種新型全營養(yǎng)作物，也是未來最具潛力的農(nóng)作物之一。然而，關(guān)于藜麥ZF-HD基因的研究至今尚未見報道。隨著藜麥基因組的公布18，使得在全基因組范圍內(nèi)分析藜麥ZF-HD基因家族成為可能。本研究以現(xiàn)有的藜麥基因組信息為基礎(chǔ)，借助生物信息學(xué)手段對藜麥ZF-HD基因家族進行全基因組鑒定，并對其基因結(jié)構(gòu)、系統(tǒng)進化關(guān)系以及ZF-HD蛋白的理化

11、性質(zhì)和保守結(jié)構(gòu)域等進行分析，同時利用轉(zhuǎn)錄組數(shù)據(jù)探究其在鹽脅迫下的表達情況，為進一步解析ZF-HD基因家族與藜麥耐鹽性的關(guān)系及其抗性遺傳改良奠定基礎(chǔ)。1材料與方法1.1藜麥ZF-HD基因家族的鑒定和藜麥ZF-HD蛋白理化性質(zhì)分析從植物轉(zhuǎn)錄因子數(shù)據(jù)庫PlantTFDB（http：/）下載擬南芥AtZF-HD轉(zhuǎn)錄因子的氨基酸序列，并在藜麥基因組數(shù)據(jù)庫Chenopodiumquinoav1.0（https：/pz/portal.html#！info？alias=Org_Cquinoa_er）進行Blastp比對，設(shè)置E-value為1e-5。在Pfam（http：/）和CDD（https：/./cdd

12、）上在線分析候選CqZF-HD蛋白的保守結(jié)構(gòu)域，去除編碼不完整蛋白質(zhì)和不含ZF或HD結(jié)構(gòu)域的蛋白質(zhì)，最終得到藜麥ZF-HD轉(zhuǎn)錄因子家族所有成員。利用ExPasy網(wǎng)站（http：/tool.html）分析藜麥ZF-HD蛋白家族的氨基酸序列長度、相對分子質(zhì)量和理論等電點。1.2藜麥ZF-HD蛋白的二級結(jié)構(gòu)和亞細胞定位利用在線網(wǎng)站GORIV（https：/npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl？page=npsa_gor4.html）分析藜麥ZF-HD蛋白的二級結(jié)構(gòu)特征。使用PSORT（http：/psort1.hgc.jp/form.html）對ZF-

13、HD蛋白進行亞細胞定位預(yù)測。1.3藜麥ZF-HD基因家族的基因結(jié)構(gòu)和藜麥ZF-HD蛋白保守結(jié)構(gòu)域分析根據(jù)藜麥ZF-HD基因的基因組序列和編碼序列，確定基因的外顯子和內(nèi)含子起止位置，并利用在線軟件GSDS2.0（http：/）繪制基因結(jié)構(gòu)圖。經(jīng)CDD分析得到藜麥ZF-HD蛋白保守結(jié)構(gòu)域的氨基酸序列起始位點，利用IBS1.0軟件19將保守結(jié)構(gòu)域位置進行可視化顯示，再使用ClustalX2.0軟件20對保守結(jié)構(gòu)域氨基酸序列進行多重比對。1.4ZF-HD基因家族的系統(tǒng)進化分析利用MEGA5.021內(nèi)置軟件ClustalW對藜麥（21個）和擬南芥（17個）ZF-HD蛋白的氨基酸序列進行多重比對，再利用M

14、EGA5.0軟件對比對結(jié)果進行系統(tǒng)進化關(guān)系分析，用鄰接法（Neighbor-Joining，NJ）構(gòu)建無根進化樹，bootstrap重復(fù)次數(shù)設(shè)置為1000次。1.5鹽脅迫下藜麥ZF-HD基因的表達分析利用本課題組已經(jīng)完成的藜麥高通量測序轉(zhuǎn)錄組數(shù)據(jù)，分析藜麥ZF-HD基因家族在鹽脅迫下的表達模式。用3%的雙氧水對供試藜麥品種37TES的種子進行表面消毒1min，然后用無菌水將種子沖洗干凈，在光周期為16h光照/8h黑暗的28恒溫光照培養(yǎng)箱進行發(fā)芽試驗。培養(yǎng)7d后，將長勢一致的幼苗轉(zhuǎn)移到Hoagland營養(yǎng)液22中進行水培，繼續(xù)培養(yǎng)28d。然后進行300mmol/LNaCl處理，分別在0h、0.5

15、h、2.0h和24.0h時取根部，用于轉(zhuǎn)錄組測序和qRT-PCR驗證，每個處理設(shè)3次重復(fù)。藜麥ZF-HD基因家族成員的轉(zhuǎn)錄組表達量經(jīng)logx2均一化處理，使用HemI軟件23繪制鹽脅迫下轉(zhuǎn)錄組基因表達圖譜。2結(jié)果與分析2.1藜麥ZF-HD基因家族的全基因組鑒定與特征分析以擬南芥AtZF-HD基因編碼的蛋白質(zhì)氨基酸序列為查詢對象，在藜麥全基因組數(shù)據(jù)庫進行Blastp比對，得到23個具有編碼ZF或HD結(jié)構(gòu)域的CqZF-HD基因，其中有2個CqZF-HD編碼不完整蛋白質(zhì)，最終鑒定出21個藜麥ZF-HD基因家族成員，分布在16條scaffold上（表1）。理化性質(zhì)分析結(jié)果表明，不同ZF-HD蛋白氨基酸

16、序列的長度差異較大，其中AUR62022438最小，為68aa，AUR62022022最大，為284aa。除AUR62027182和AUR62044191因蛋白質(zhì)氨基酸序列含有幾個連續(xù)未定義的氨基酸而無法計算其相對分子質(zhì)量和等電點（pI）外，其余19個ZF-HD蛋白的相對分子質(zhì)量為7.5410330.20103，pI為6.019.54。pI小于7的只有AUR62022438（6.68）和AUR62039796（6.01），為酸性蛋白質(zhì)，其他ZF-HD的pI均大于7，為堿性蛋白質(zhì)。2.2藜麥ZF-HD蛋白的二級結(jié)構(gòu)和亞細胞定位預(yù)測藜麥ZF-HD基因編碼的蛋白質(zhì)二級結(jié)構(gòu)主要包括-螺旋、延伸鏈和無規(guī)

17、則卷曲，均不含-折疊，其中AUR62022682、AUR62022414、AUR62022438、AUR62039796、AUR62028772和AUR62022926也沒有-螺旋結(jié)構(gòu)（表2）。大部分ZF-HD成員的二級結(jié)構(gòu)成分中，無規(guī)則卷曲占比最高，其次是延伸鏈，-螺旋占比相對較少。亞細胞定位預(yù)測結(jié)果顯示，14個ZF-HD成員定位于細胞質(zhì)，3個成員定位于細胞核，3個成員定位于線粒體基質(zhì)，1個成員定位于葉綠體基質(zhì)，說明它們可能調(diào)控著植物不同的生物學(xué)過程。2.3ZF-HD基因家族的系統(tǒng)進化關(guān)系分析為更好地了解藜麥ZF-HD基因家族成員之間的同源性關(guān)系，將17個擬南芥AtZF-HD基因與21個藜麥

18、CqZF-HD基因編碼的蛋白質(zhì)氨基酸序列進行多重比對并構(gòu)建系統(tǒng)進化樹（圖1）。聚類分析結(jié)果顯示，38個ZF-HD轉(zhuǎn)錄因子可以劃分為5個亞家族：MIF、ZHD、ZHD、ZHD和ZHD，位于不同亞家族的蛋白質(zhì)可能發(fā)揮著不同的生物學(xué)作用。藜麥的ZF-HD基因在MIF亞家族數(shù)目較多，為12個，在ZHD存在4個，在ZHD和ZHD均含2個，在ZHD只存在1個，即AUR62000606。藜麥的AUR62022575和AUR62022062與擬南芥的AT2G18350和AT1G75240處于同一分支，說明它們的親緣關(guān)系較近，可能具有類似的功能。2.4藜麥ZF-HD基因結(jié)構(gòu)特征分析基因結(jié)構(gòu)分析結(jié)果表明，藜麥ZF

19、-HD基因家族的基因結(jié)構(gòu)比較簡單，均不含非翻譯區(qū)，但不同成員間也存在一定差異（圖2）。其中，13個基因只有1個外顯子，沒有內(nèi)含子;7個基因具有2個外顯子和1個內(nèi)含子;1個基因（AUR62027182）具有4個外顯子和3個內(nèi)含子。AUR62022022和AUR62027848同屬ZHD亞家族，有著相同的基因結(jié)構(gòu)和相似的CDS長度，但因內(nèi)含子長度差別較大使其基因全長存在顯著差異。為進一步研究藜麥ZF-HD蛋白的特征區(qū)域，在CDD網(wǎng)站上對其進行保守結(jié)構(gòu)域分析，結(jié)果表明21個ZF-HD蛋白均含有ZF或HD結(jié)構(gòu)域，但這2個結(jié)構(gòu)域在各個亞家族成員的分布有所不同（圖3）。MIF亞家族所有成員HD結(jié)構(gòu)域缺失，

20、均只含有1個ZF結(jié)構(gòu)域，并且同源基因?qū)幋a蛋白質(zhì)的氨基酸數(shù)量和結(jié)構(gòu)域位置都具有較高的相似性;ZHD亞家族中除AUR62044191只含有1個ZF結(jié)構(gòu)域外，其他3個成員都含有ZF和HD結(jié)構(gòu)域;ZHD和ZHD的成員都含有ZF和HD結(jié)構(gòu)域;ZHD的2個成員的保守結(jié)構(gòu)域均位于蛋白質(zhì)的N端，AUR62039796只含有1個HD結(jié)構(gòu)域，而AUR62028772只含有1個ZF結(jié)構(gòu)域。對藜麥ZF-HD蛋白保守結(jié)構(gòu)域氨基酸序列進行多重比對，發(fā)現(xiàn)ZF結(jié)構(gòu)域由3155個氨基酸組成，HD結(jié)構(gòu)域由3158個氨基酸組成，并且它們都有多個保守的氨基酸位點，HD結(jié)構(gòu)域在藜麥ZF-HD家族中的保守性比ZF結(jié)構(gòu)域更強（圖4）。保

21、守區(qū)域可能通過與其他DNA結(jié)合或互作，從而導(dǎo)致ZF-HD轉(zhuǎn)錄因子的功能多樣性。2.5藜麥ZF-HD家族基因在鹽脅迫下的表達模式分析為挖掘藜麥耐鹽相關(guān)ZF-HD基因，提供藜麥耐鹽遺傳改良基因資源，利用課題組前期完成的轉(zhuǎn)錄組數(shù)據(jù)，對ZF-HD家族基因在鹽脅迫下的表達模式進行分析。結(jié)果如圖5所示，藜麥ZF-HD家族基因在鹽脅迫處理后的表達模式存在差異，其中10個基因表達量未發(fā)生變化，即不參與對鹽脅迫的響應(yīng)，11個基因表達量在不同處理時間發(fā)生不同程度的變化。與對照（0h）相比，AUR62027848、AUR62022414和AUR62028772基因在各個時間點都呈上調(diào)表達，而AUR62022062、

22、AUR62000606和AUR62023917則呈下調(diào)表達。鹽脅迫0.5h時，AUR62022022和AUR62022575的表達量上調(diào)至最高值，在2.0h和24.0h時表達量有所下降并明顯低于對照，說明這2個基因可能迅速響應(yīng)早期鹽脅迫。鹽脅迫下AUR62022682的表達量呈先降低后升高再降低的趨勢，在2.0h時達到最高峰;AUR62039798和AUR62039796的表達量呈先降低后上升的趨勢，在24.0h達到最大值。為驗證轉(zhuǎn)錄組數(shù)據(jù)的可靠性，課題組已對用于測序的樣品進行了qRT-PCR分析（圖6），結(jié)果表明兩者具有很強的相關(guān)性（R2=0.7973），證明RNA-seq表達數(shù)據(jù)是可靠的2

23、4。同時，ZF-HD家族基因在鹽脅迫下呈現(xiàn)出的不同表達模式也說明它們在藜麥鹽響應(yīng)過程中可能發(fā)揮著不同的作用。3討論土壤鹽分是干旱半干旱地區(qū)現(xiàn)代農(nóng)業(yè)的主要非生物脅迫之一，嚴重影響種子萌發(fā)、植物生長發(fā)育和作物產(chǎn)量25-26。據(jù)統(tǒng)計，約有8108hm2的土地（占世界土地總面積的6.5%）受到鹽脅迫的危害27。為了適應(yīng)鹽脅迫，植物進化出了一系列的響應(yīng)機制來提高其耐鹽性，包括復(fù)雜的生理特性、代謝和信號通路以及分子調(diào)控網(wǎng)絡(luò)等。隨著植物基因組測序的完成及全基因組數(shù)據(jù)庫的公布，越來越多的耐鹽基因被挖掘出來，其中鑒定和分析植物基因家族就是挖掘關(guān)鍵功能基因的一種有效手段。植物特異性ZF-HD基因編碼一個轉(zhuǎn)錄因子家

24、族，該家族存在于陸生植物的主要類群中，包括維管植物和非維管植物，但在原核生物、綠藻和真菌中沒有發(fā)現(xiàn)。ZF-HD基因家族參與多種過程，包括植物發(fā)育和生理過程，以及對生物脅迫和非生物脅迫的抗性。目前，在擬南芥、水稻、玉米和棉花等多種植物中均已鑒定出ZF-HD基因家族，但關(guān)于藜麥ZF-HD基因家族的研究至今尚未見報道。在本研究中，我們鑒定并分析了藜麥ZF-HD基因的全基因組及其在鹽脅迫反應(yīng)中的表達。在藜麥ZF-HD基因家族中共鑒定出21個成員，其數(shù)量要多于模式植物擬南芥（17）和水稻（15），可能是因為基因組復(fù)制事件導(dǎo)致了異源四倍體物種藜麥ZF-HD基因的擴增。21個ZF-HD基因不均勻地分布在藜麥

25、16條scaffold上，其編碼蛋白質(zhì)的氨基酸數(shù)量和相對分子質(zhì)量存在差異，并以堿性蛋白質(zhì)居多?；蛲怙@子和內(nèi)含子結(jié)構(gòu)的差異主要有3種機制，即外顯子或內(nèi)含子的獲得或丟失，外顯子化或假外顯子化和插入或缺失28-29。藜麥所有ZF-HD成員中，13個基因沒有內(nèi)含子，7個基因只有1個內(nèi)含子，1個基因具有3個內(nèi)含子，這種差異可能是上述3種機制參與所致。有研究結(jié)果表明，植物基因內(nèi)含子的數(shù)量越少，其對外界環(huán)境的適應(yīng)性越強30。將擬南芥和藜麥的ZF-HD蛋白構(gòu)建系統(tǒng)進化樹，38個ZF-HD轉(zhuǎn)錄因子劃分為5個亞家族，2種植物ZF-HD在每個亞家族均有分布，說明基因復(fù)制是植物基因組進化中基因家族擴增的主要驅(qū)動力。

26、藜麥ZF-HD蛋白在結(jié)構(gòu)上與其他物種是相似的，含有ZF和HD結(jié)構(gòu)域，但在亞家族分布上也有不同之處，比如與擬南芥1和小麥2等植物相比，藜麥的MIF亞家族成員同樣只含有ZF結(jié)構(gòu)域，而ZHD亞家族并非都同時含有ZF和HD結(jié)構(gòu)域。氨基酸序列分析結(jié)果顯示，藜麥ZF和HD結(jié)構(gòu)域都具有高度保守性，并且HD結(jié)構(gòu)域的保守性更強一些，結(jié)構(gòu)域中保守區(qū)域可能通過與其他DNA結(jié)合或互作，從而導(dǎo)致ZF-HD轉(zhuǎn)錄因子的功能多樣性。二級結(jié)構(gòu)和亞細胞定位分析結(jié)果顯示，大部分ZF-HD成員中無規(guī)則卷曲占比最高，其次是延伸鏈，-螺旋占比相對較少，其中3對同源基因不含-螺旋;14個ZF-HD屬于胞質(zhì)蛋白，分別有3個成員定位于細胞核和

27、線粒體基質(zhì)，1個定位于葉綠體基質(zhì)，說明它們在參與藜麥生物學(xué)過程中可能發(fā)揮著不同的作用。植物以一種復(fù)雜的方式響應(yīng)鹽脅迫，涉及許多基因、蛋白質(zhì)以及代謝和信號通路的協(xié)同功能與相互作用31-32?；虮磉_調(diào)控是植物應(yīng)對和適應(yīng)鹽脅迫的關(guān)鍵現(xiàn)象之一33-35。相關(guān)研究結(jié)果表明，ZF-HD在植物耐鹽脅迫中扮演著重要的角色10，17。本研究利用課題組已經(jīng)完成的藜麥在300mmol/LNaCl處理下不同時間點的高通量測序轉(zhuǎn)錄組數(shù)據(jù)，分析了所有ZF-HD基因在響應(yīng)鹽脅迫的表達模式，發(fā)現(xiàn)10個成員不參與響應(yīng)鹽脅迫，另外11個成員的表達量在不同處理時間呈現(xiàn)出不同變化，說明它們在藜麥響應(yīng)鹽脅迫過程中發(fā)揮著不同的作用。A

28、UR62027848、AUR62022414和AUR62028772基因在鹽處理0.524.0h時表達量受到明顯誘導(dǎo)，這些基因可能在藜麥抵抗鹽脅迫過程中發(fā)揮著重要的作用。而AUR62022062、AUR62000606和AUR62023917為鹽敏感基因，鹽脅迫明顯抑制了其表達。AUR62022022和AUR62022575在鹽脅迫0.5h時表達量較高，在脅迫后期顯著下降，說明這2個基因可在早期迅速響應(yīng)鹽脅迫，可能直接受到鹽脅迫信號的誘導(dǎo)。關(guān)于ZF-HD在藜麥抵抗鹽脅迫中的功能及作用機理仍需要借助遺傳轉(zhuǎn)化手段進行驗證。4結(jié)論本研究對藜麥ZF-HD轉(zhuǎn)錄因子進行了全基因組鑒定，共獲得21個家族成員

29、，劃分為5個亞家族，都具有典型的ZF或HD保守結(jié)構(gòu)域。21個ZF-HD基因在鹽脅迫處理下呈現(xiàn)出不同的表達模式，推測不同ZF-HD家族基因在藜麥響應(yīng)鹽脅迫過程中可能發(fā)揮著不同的作用。參考文獻：1TANQK，IRISHVF.TheArabidopsiszincfinger-homeodomaingenesencodeproteinswithuniquebiochemicalpropertiesthatarecoordinatelyexpressedduringfloraldevelopmentJ.PlantPhysiology，2022，140：1095-1108.2LIUH，YANGY，ZHAN

30、GLS.Zincfinger-homeodomaintranscriptionalfactors（ZF-HDs）inwheat（TriticumaestivumL.）：identification，evolution，expressionanalysisandresponsetoabioticstressesJ.Plants，2022，10：593.3HUW，DEPAMPHILISCW，MAH.Phylogeneticanalysisoftheplant-specificzincfinger-homeoboxandminizincfingergenefamiliesJ.JournalofInt

31、egrativePlantBiology，2022，50：1031-1045.4WANGW，WUP，LIY，etal.Genome-wideanalysisandexpressionpatternsofZF-HDtranscriptionfactorsunderdifferentdevelopmentaltissuesandabioticstressesinChinesecabbageJ.MolecularGeneticsandGenomics，2022，291：1451-1464.5BHATTACHARJEEA，GHANGALR，GARGR，etal.Genome-wideanalysiso

32、fhomeoboxgenefamilyinlegumes：Identification，geneduplicationandexpressionprofilingJ.PLoSOne，2022，10：e0119198.6KRISHNASS，MAJUMDARI，GRISHINNV.Structuralclassificationofzincfingers：SurveyandsummaryJ.NucleicAcidsResearch，2022，31：532-550.7MACKAYJP，CROSSLEYM.ZincfingersarestickingtogetherJ.TrendsinBiochemi

33、calSciences，1998，23：1-4.8WINDHOVELA，HEINI，DABROWAR，etal.CharacterizationofanovelclassofplanthomeodomainproteinsthatbindtotheC4phosphoenolpyruvatecarboxylasegeneofFlaveriatrinerviaJ.PlantMolecularBiology，2022，45：201-214.9JRGENSENJE，GRNLUNDM，PALLISGAARDN，etal.Anewclassofplanthomeoboxgenesisexpressedin

34、specificregionsofdeterminatesymbioticrootnodulesJ.PlantMolecularBiology，1999，40：65-77.10KHATUNK，NATHUK，ROBINA，etal.Genome-wideanalysisandexpressionprofilingofzincfingerhomeodomain（ZHD）familygenesreveallikelyrolesinorgandevelopmentandstressresponsesintomatoJ.BMCGenomics，2022，18：695.11ABDULLAHM，CHENGX

35、，CAOY，etal.Zincfinger-homeodomaintranscriptionalfactors（ZHDs）inuplandcotton（Gossypiumhirsutum）：Genome-wideidentificationandexpressionanalysisinfiberdevelopmentJ.FrontiersinGenetics，2022，9：357.12LIUM，WANGX，SUNW，etal.Genome-wideinvestigationoftheZF-HDgenefamilyinTartarybuckwheat（Fagopyrumtataricum）J.B

36、MCPlantBiology，2022，19：248.13TRANLS，NAKASHIMAK，SAKUMAY，etal.Co-expressionofthestress-induciblezincfingerhomeodomainZFHD1andNACtranscriptionfactorsenhancesexpressionoftheERD1geneinArabidopsisJ.PlantJournal，2022，49：46-63.14WANGL，HUAD，HEJ，etal.AuxinResponseFactor2（ARF2）anditsregulatedhomeodomaingeneHB3

37、3mediateabscisicacidresponseinArabidopsisJ.PLoSGenetics，2022，7：e1002172.15XUY，WANGY，LONGQ，etal.OverexpressionofOsZHD1，azincfingerhomeodomainclasshomeoboxtranscriptionfactor，inducesabaxiallycurledanddroopingleafinriceJ.Planta，2022，239：803-816.16YOONJ，CHOLH，YANGW，etal.HomeoboxtranscriptionfactorOsZHD2

38、promotesrootmeristemactivityinricebyinducingethylenebiosynthesisJ.JournalofExperimentalBotany，2022，71（18）：5348-5364.17WANGH，YINX，LIX，etal.Genome-wideidentification，evolutionandexpressionanalysisofthegrape（VitisviniferaL.）zincfinger-homeodomaingenefamilyJ.InternationalJournalofMolecularSciences，2022，

39、15：5730-5748.18JARVISDE，HOYS，LIGHTFOOTDJ，etal.ThegenomeofChenopodiumquinoaJ.Nature，2022，542：307-312.19LIUW，XIEY，MAJ，etal.IBS：anillustratorforthepresentationandvisualizationofbiologicalsequencesJ.Bioinformatics，2022，31（20）：3359-3361.20THOMPSONJD，GIBSONTJ，PLEWNIAKF，etal.TheCLUSTAL_Xwindowsinterface：flexiblestrategiesformultiplesequencealignmentaidedbyqualityanalysistoolsJ.NucleicAcidsResearch，1997，25（24）：4876-4882.21TAMURAK，PETERSOND，PETERSONN，etal.MEGA5：molecularevolutionarygeneticsanalysisusingmaximumlikelihood，evolutionarydistance，andmaximumparsimonymethodsJ.MolecularBiolog