連 輝， 江淑珍，b， 熊遠(yuǎn)芳，b， 沈 軍， 陳世品，①

(福建農(nóng)林大學(xué)： a. 林學(xué)院， b. 園林學(xué)院 蘭科植物保護(hù)與利用國家林業(yè)和草原局重點(diǎn)實(shí)驗(yàn)室， 福建 福州 350002)

在被子植物中，花是重要的生殖器官，也是進(jìn)化過程中最具多樣性的結(jié)構(gòu)[1]。同一個(gè)體(雌雄同株)或不同個(gè)體(雌雄異株)的單性雄花和雌花可促進(jìn)異花授粉和基因交流，從而減少近親繁殖[2]。殼斗科(Fagaceae)含900多種，是重要的經(jīng)濟(jì)資源，對森林生態(tài)系統(tǒng)的動態(tài)平衡有重要作用[3]。單性花是殼斗科植物的主要繁殖策略[4]。在殼斗科植物中，不同屬植物有著不同類型花序，如栗屬(CastaneaMill.)和錐屬(CastanopsisSpach)植物，前者具有雄性和兩性的柔荑花序，后者具有雄性和雌性的花序[4]。在栗屬植物中，板栗(CastaneamollissimaBlume)雌花中有12枚發(fā)育不完全的退化雄蕊，具有短的花絲和可見的花藥[5]。在錐栗〔Castaneahenryi(Skan) Rehd. et Wils.〕雌花中也存在退化雄蕊，在雌蕊原基的發(fā)育過程中雄蕊原基也發(fā)育，隨著雌蕊原基的分化完成雄蕊停止發(fā)育[6]。而在其他殼斗科種類的研究中，單性花由花分生組織產(chǎn)生，不能同時(shí)形成雌性和雄性的性器官原基[4]。因此，在殼斗科植物中，栗屬植物由于雄蕊停止生長而形成單性花，而同科其他屬種類的單性花不存在雌蕊或雄蕊敗育的現(xiàn)象，花發(fā)育起始階段就是單性的。在這些花中，性別決定基因通過控制花分生組織的發(fā)生和花器官的形成來調(diào)控性別分化[7]。目前，雌雄同株植物的花器官發(fā)育已被廣泛研究，在大多數(shù)情況下，花器官分化主要由MADS-box轉(zhuǎn)錄因子控制[8]。

米櫧〔Castanopsiscarlesii(Hemsl.) Hayata〕為殼斗科錐屬植物，廣泛分布于中亞熱帶海拔1 300 m以下的山地丘陵，為亞熱帶常綠闊葉林頂級群落之一[9]。米櫧具有2種類型花序：雄性柔荑花序在當(dāng)年生的葉軸上發(fā)育，雌花序在枝稍附近發(fā)育。目前，在米櫧的雄花中尚未發(fā)現(xiàn)雌蕊敗育的跡象，表明米櫧雄花在花發(fā)育的起始階段就是單性的。在殼斗科植物中，花發(fā)育起始階段就是單性的現(xiàn)象較為普遍，如歐洲栓皮櫟(QuercussuberLinn.)[10]。然而，關(guān)于殼斗科植物從起始階段就是單性的雌、雄花發(fā)育機(jī)制的研究還不夠深入。本文基于轉(zhuǎn)錄組測序技術(shù)，對米櫧雄花的花芽期、半開期和盛花期3個(gè)時(shí)期進(jìn)行轉(zhuǎn)錄組分析，對差異表達(dá)unigenes進(jìn)行功能注釋。通過系統(tǒng)發(fā)育分析，鑒定出米櫧MADS-box家族基因，并對其在不同發(fā)育階段雄花的表達(dá)量進(jìn)行分析。以期進(jìn)一步認(rèn)識MADS-box家族基因在米櫧雄花發(fā)育中的作用，為深入探究殼斗科植物單性花發(fā)育的分子機(jī)制提供參考。

1 材料和方法

1.1 材料

供試米櫧雄花來自于福建省福州市鼓山風(fēng)景區(qū)(東經(jīng)119°23′04″、北緯26°03′13″，海拔402 m)米櫧野生群落。樹種鑒定參照《中國植物志》[11]和《福建植物志》[12]。于2019年3月至5月，持續(xù)觀察米櫧雄花發(fā)育時(shí)期，分別在米櫧雄花的花芽期(花被片和雄蕊均未分化，F(xiàn)1期)、半開期(花被片和雄蕊分化完成，花被片微微張開，F(xiàn)2期)及盛花期(花被片張開并開始散粉，F(xiàn)3期)采集雄花花序，每個(gè)時(shí)期分別從3株植株上采集，即為3個(gè)生物學(xué)重復(fù)，共采集9個(gè)樣品。取樣后立即用液氮速凍，保存于-80 ℃冰箱中，用于后續(xù)轉(zhuǎn)錄組測序分析。

1.2 方法

1.2.1 轉(zhuǎn)錄組測序與組裝 使用PureLinkTMPlant RNA試劑盒(美國Invitrogen公司)分別提取米櫧3個(gè)時(shí)期9個(gè)樣品的總RNA，再使用RNase-FreeDNase (美國Promega公司)去除DNA污染。采用質(zhì)量體積分?jǐn)?shù)1.25%瓊脂糖凝膠電泳檢測總RNA的完整性，使用NanoDrop 1000超微量分光光度計(jì)(美國Thermo Fisher Scientific公司)測定總RNA的濃度和純度。使用寡聚體(dT)貼附磁珠對總RNA進(jìn)行純化，將純化后的總RNA用片段緩沖液破碎成小塊，利用隨機(jī)的六聚體引物反轉(zhuǎn)錄產(chǎn)生第1鏈cDNA，接著合成第2鏈cDNA。雙鏈DNA末端補(bǔ)平，5′端磷酸化，3′端進(jìn)行接頭。使用PCR擴(kuò)增雙鏈cDNA片段，用AMPure XP磁珠(美國Thermo Fisher Scientific公司)純化擴(kuò)增產(chǎn)物，純化后的雙鏈PCR產(chǎn)物經(jīng)熱變性，溶解于EB溶液中，經(jīng)安捷倫2100生物分析儀(美國Agilent公司)質(zhì)量控制驗(yàn)證后，溶解于EB溶液中的雙鏈PCR產(chǎn)物經(jīng)95 ℃熱變性成單鏈，經(jīng)寡核苷酸序列循環(huán)，最終得到單鏈環(huán)狀DNA作為文庫。采用第二代測序技術(shù)，并基于Illumina HiSeq測序平臺，對這些文庫進(jìn)行雙末端測序。使用Trinity軟件[13]進(jìn)行轉(zhuǎn)錄組從頭組裝，為確保RNA測序數(shù)據(jù)的準(zhǔn)確性和可靠性，從raw reads中剔除低質(zhì)量、接頭污染以及未知堿基N含量過高的reads，僅對剩下的高質(zhì)量測序數(shù)據(jù)(clean reads)進(jìn)行統(tǒng)計(jì)分析。

1.2.2 轉(zhuǎn)錄組差異表達(dá)unigenes注釋 采用FPKM值表示轉(zhuǎn)錄組差異表達(dá)unigenes的表達(dá)豐度。采用Benjamini-Hochberg校正法對原有假設(shè)檢驗(yàn)得到的P值進(jìn)行校正，最終采用FDR(錯(cuò)誤發(fā)現(xiàn)率)作為差異表達(dá)unigenes篩選的關(guān)鍵指標(biāo)。利用DESeq軟件，以FDR<0.05、FC(差異倍數(shù))≥2作為差異表達(dá)unigenes的篩選標(biāo)準(zhǔn)，差異表達(dá)unigenes的功能注釋信息通過KEGG數(shù)據(jù)庫(http：∥www.genome.jp/kegg/)和GO數(shù)據(jù)庫(http：∥www. geneontology.org/)獲得。使用BLAST軟件將差異表達(dá)unigenes序列與GO數(shù)據(jù)庫比對，進(jìn)行注釋，在獲得每個(gè)差異表達(dá)unigene的GO注釋后，使用WEGO軟件(http：∥wego.genomics.org.cn/cgi-bin/wego/index.pl)對所有差異表達(dá)unigenes進(jìn)行GO功能分類統(tǒng)計(jì)，通過KEGG數(shù)據(jù)庫的注釋信息得到基因的path-way注釋。

1.2.3 MADS-box家族基因鑒定與表達(dá)量分析 在TAIR在線數(shù)據(jù)庫(https：∥www.arabidopsis.org/index.jsp)下載擬南芥〔Arabidopsisthaliana(Linn.) Heynh.〕MADS-box家族蛋白序列。以擬南芥MADS-box家庭蛋白序列作為索引在米櫧基因組中進(jìn)行本地BLASTp檢索。將獲得的候選序列提交到NCBI-CDD網(wǎng)站(https：∥www.ncbi.nlm.nih.gov/cdd)和HMMER網(wǎng)站(https：∥www.ebi.ac.uk/Tools/hmmer/search/hmmsearch)，根據(jù)Pfam(PF00319)在米櫧基因組數(shù)據(jù)庫中篩選出具有完整保守結(jié)構(gòu)域的成員作為米櫧MADS-box家族基因。由于米櫧type-Ⅰ型MADS-box家族基因在花發(fā)育過程中幾乎不表達(dá)，且該家族基因在花發(fā)育過程中調(diào)控作用的研究較少，因此，未開展該家族基因的系統(tǒng)進(jìn)化樹及表達(dá)量分析。采用MEGA 7.0軟件對米櫧和擬南芥的type-Ⅱ型MADS-box家族基因的蛋白序列進(jìn)行多序列比對后，采用鄰接法(neighbor-joining method)構(gòu)建系統(tǒng)進(jìn)化樹，對建成的系統(tǒng)進(jìn)化樹進(jìn)行自展法系數(shù)(bootstrap)校正，設(shè)置1 000次重復(fù)檢驗(yàn)，最后優(yōu)化進(jìn)化樹。從轉(zhuǎn)錄組數(shù)據(jù)中獲取米櫧雄花3個(gè)時(shí)期MADS-box家族基因的表達(dá)豐度值(FPKM)，通過log2(FPKM+1)計(jì)算表達(dá)量，并使用TBtools軟件中的HeatMap程序繪制基因表達(dá)量熱圖。

2 結(jié)果和分析

2.1 測序質(zhì)量分析和差異表達(dá)unigenes篩選

2.1.1 測序質(zhì)量分析 測序結(jié)果(表1)顯示：米櫧雄花3個(gè)時(shí)期樣品的clean reads為628 800 000～643 700 000，高質(zhì)量總堿基量為6.29～6.44 Gb，Q30堿基百分比均不小于90.67%，GC含量為45.09%～45.39%。

2.1.2 差異表達(dá)unigenes篩選 對花芽期、半開期和盛花期3個(gè)時(shí)期米櫧雄花的轉(zhuǎn)錄組結(jié)果進(jìn)行比對分析，結(jié)果(表2)顯示：花芽期與半開期間有15 086個(gè)差異表達(dá)unigenes，其中，與花芽期相比，半開期表達(dá)量上調(diào)的unigenes有5 283個(gè)，表達(dá)量下調(diào)的unigenes有9 803個(gè)。半開期與盛花期間有6 770個(gè)差異表達(dá)unigenes，其中，與半開期相比，盛花期表達(dá)量上調(diào)的unigenes有3 536個(gè)，表達(dá)量下調(diào)的unigenes有3 234個(gè)?；ㄑ科谂c盛花期間有16 150個(gè)差異表達(dá)unigenes，其中，與花芽期相比，盛花期表達(dá)量上調(diào)的unigenes有6 588個(gè)，表達(dá)量下調(diào)的unigenes有9 562個(gè)。

表1 米櫧雄花轉(zhuǎn)錄組數(shù)據(jù)統(tǒng)計(jì)

表2 不同時(shí)期間米櫧雄花轉(zhuǎn)錄組差異表達(dá)unigenes數(shù)

2.2 功能注釋分析

2.2.1 GO功能注釋分析 在GO數(shù)據(jù)庫對米櫧雄花的花芽期與半開期、半開期與盛花期以及花芽期與盛花期3組的差異表達(dá)unigenes進(jìn)行比對和功能注釋，結(jié)果見圖1。

由圖1可以看出：米櫧雄花的花芽期與半開期、半開期與盛花期以及花芽期與盛花期3組分別有31 888、13 228和35 203個(gè)差異表達(dá)unigenes注釋在分子功能、細(xì)胞組分和生物過程3大類47個(gè)亞類中。在分子功能的13個(gè)亞類中，3組差異表達(dá)unigenes主要注釋在催化活性和連接2個(gè)亞類中，其中，注釋在催化活性亞類的差異表達(dá)unigenes分別有4 411、2 001和4 876個(gè)；注釋在連接亞類的差異表達(dá)unigenes分別有4 437、1 911和4 896個(gè)。在細(xì)胞組分的14個(gè)亞類中，3組差異表達(dá)unigenes主要注釋在膜、膜要素、細(xì)胞和細(xì)胞器4個(gè)亞類中，其中，注釋在膜亞類的差異表達(dá)unigenes分別有3 215、1 493和3 464個(gè)；注釋在膜要素亞類的差異表達(dá)unigenes分別有2 997、1 399和3 206個(gè)；注釋在細(xì)胞亞類的差異表達(dá)unigenes分別有2 594、997和2 915個(gè)；注釋在細(xì)胞器亞類的差異表達(dá)unigenes分別有1 862、674和2 111個(gè)。在生物過程的20個(gè)亞類中，3組的差異表達(dá)unigenes主要注釋在細(xì)胞過程、代謝過程、生物調(diào)節(jié)和應(yīng)激反應(yīng)4個(gè)亞類中，其中，注釋在細(xì)胞過程亞類的差異表達(dá)unigenes分別有2 738、983和3 071個(gè)；注釋在代謝過程亞類的差異表達(dá)unigenes分別有2 347、690和2 403個(gè)；注釋在生物調(diào)節(jié)亞類的差異表達(dá)unigenes分別有1 033、398和1 121個(gè)；注釋在應(yīng)激反應(yīng)亞類的差異表達(dá)unigenes分別有981、393和1 037個(gè)。

2.2.2 KEGG代謝通路分析 在KEGG數(shù)據(jù)庫中對米櫧雄花的花芽期與半開期、半開期與盛花期以及花芽期與盛花期3組的差異表達(dá)unigenes進(jìn)行比對和功能注釋，結(jié)果見圖2。

由圖2可以看出：花芽期與半開期、半開期與盛花期以及花芽期與盛花期3組分別有10 826、5 241和11 896個(gè)差異表達(dá)unigenes注釋在KEGG代謝通路，共6大類21亞類，且主要注釋在代謝和遺傳信息處理2大類中。3組注釋在代謝相關(guān)通路的差異表達(dá)unigenes最多，分別有6 580、3 414和7 130個(gè)，包括全局及概要圖、碳水化合物代謝、其他次生代謝產(chǎn)物的生物合成、氨基酸代謝和脂質(zhì)代謝等11個(gè)亞類，其中，注釋在全局及概要圖亞類的差異表達(dá)unigenes最多，分別有2 473、1 266和2 673個(gè)。3組注釋在遺傳信息處理相關(guān)通路的差異表達(dá)unigenes分別有2 060、843和2 379個(gè)，包括翻譯、轉(zhuǎn)錄、復(fù)制和修復(fù)以及折疊、分類和降解4個(gè)亞類。3組注釋在環(huán)境信息處理相關(guān)通路的差異表達(dá)unigenes分別有846、384和943個(gè)，包括信號轉(zhuǎn)錄和膜運(yùn)輸2個(gè)亞類。在生物系統(tǒng)相關(guān)通路中，3組分別有846、384和943個(gè)差異表達(dá)unigenes注釋在環(huán)境適應(yīng)性亞類中。在細(xì)胞過程相關(guān)通路中，3組分別有372、171和443個(gè)差異表達(dá)unigenes注釋在運(yùn)輸和分解代謝亞類中。在人類疾病相關(guān)通路中，3組共有32、8和40個(gè)差異表達(dá)unigenes注釋在內(nèi)分泌和代謝性疾病以及耐藥性：抗菌2個(gè)亞類中。

2.3 type-Ⅱ型MADS-box家族基因鑒定及表達(dá)量分析

利用MEGA 7.0對米櫧和擬南芥type-Ⅱ型MADS-box蛋白構(gòu)建系統(tǒng)進(jìn)化樹(圖3)。結(jié)果顯示：米櫧type-Ⅱ型MADS-box家族基因有53個(gè)成員，可分為13類：C/D類(6個(gè)成員)、ANR1類(4個(gè)成員)、AGL15類(4個(gè)成員)、AP3/PI類(1個(gè)成員)、Bs類(1個(gè)成員)、MIKC*類(7個(gè)成員)、SVP類(14個(gè)成員)、FLC類(2個(gè)成員)、SOC1類(3個(gè)成員)、AP1/FUL類(3個(gè)成員)、AGL12類(2個(gè)成員)、AGL6類(2個(gè)成員)和SEP類(4個(gè)成員)。

： 花芽期與半開期的差異表達(dá)unigenes的代謝通路Metabolic pathway of differentially expressed unigenes between flower bud stage and half flowering stage； ： 半開期與盛花期的差異表達(dá)unigenes的代謝通路Metabolic pathway of differentially expressed unigenes between half flowering stage and flowering stage； ： 花芽期與盛花期的差異表達(dá)unigenes的代謝通路Metabolic pathway of differentially expressed unigenes between flower bud stage and flowering stage.

基于花芽期、半開期和盛花期3個(gè)時(shí)期米櫧雄花的轉(zhuǎn)錄組數(shù)據(jù)對米櫧type-Ⅱ型MADS-box家族基因進(jìn)行表達(dá)分析，發(fā)現(xiàn)其在雄花不同發(fā)育階段的表達(dá)具有差異性。3個(gè)時(shí)期米櫧雄花type-Ⅱ型MADS-box家族基因的表達(dá)量熱圖(圖4)顯示：在C/D類同源基因中，SHP基因(Maker00006492)和AG基因(Maker00069044)表達(dá)趨勢類似，在花芽期表達(dá)量較低，半開期和盛花期表達(dá)量升高；其余4個(gè)C/D類同源基因在雄花發(fā)育過程中幾乎不表達(dá)。在AP3/PI類同源基因中，AP3基因(Maker00011188)隨雄花發(fā)育表達(dá)量逐漸升高。在MIKC*類同源基因中， Maker00016147和Maker00054870在雄花發(fā)育過程中表達(dá)量較高，其余基因表達(dá)量較低或不表達(dá)；其中，Maker00016147在花芽期的表達(dá)量較低，在半開期的表達(dá)量升至最高，在盛花期的表達(dá)量降至半開期的一半；Maker00054870的表達(dá)趨勢類似，在半開期的表達(dá)量最高，在盛花期的表達(dá)量略下降。在AP1/FUL類同源基因中，Maker00031743和Maker00060929為AP1基因，Maker00031743在花芽期的表達(dá)量最高，在半開期的表達(dá)量降至最低，在盛花期的表達(dá)量較半開期升高；Maker00060929在花芽期、半開期和盛開期的表達(dá)量均較高，且?guī)缀鯖]有變化；FUL基因(Maker00039576)在花芽期和盛花期的表達(dá)量較高，但在半開期的表達(dá)量卻較低。在AGL6類同源基因中，AG6基因(Maker00067599)在雄花發(fā)育3個(gè)時(shí)期的表達(dá)量均較高。在SEP類同源基因中， Maker00050389、Maker00039517、Maker00060937和Maker00031685均在雄花發(fā)育過程中高表達(dá)；其中，SEP3基因(Maker00050389)的表達(dá)量最高，且在雄花發(fā)育3個(gè)時(shí)期的波動不大；Maker00039517和Maker00031685的表達(dá)趨勢相似，在花芽期的表達(dá)量較高，在半開期和盛花期的表達(dá)量逐漸升高；Maker00060937在花芽期的表達(dá)量較低，在半開期和盛花期的表達(dá)量逐漸升高。其余6類同源基因在米櫧雄花發(fā)育過程中幾乎不表達(dá)。

： 米櫧Castanopsis carlesii (Hemsl.) Hayata； ： 擬南芥Arabidopsis thaliana (Linn.) Heynh.

F1： 花芽期Flower bud stage； F2： 半開期Half flowering stage； F3： 盛花期Flowering stage.

3 討 論

功能注釋結(jié)果表明：通過與GO數(shù)據(jù)庫比對，花芽期與半開期、半開期與盛花期以及花芽期與盛花期3組差異表達(dá)unigenes注釋在分子功能、細(xì)胞組分和生物過程的47個(gè)亞類中。在分子功能中，注釋在催化活性和連接亞類的差異表達(dá)unigenes最多，可見催化酶在米櫧雄花發(fā)育過程中發(fā)揮著重要作用；在細(xì)胞組分中，差異表達(dá)unigenes主要注釋在膜、膜要素、細(xì)胞和細(xì)胞器4個(gè)亞類中，在雄花發(fā)育過程中，花粉粒的成熟伴隨著新的細(xì)胞膜、細(xì)胞壁以及細(xì)胞器的生成，推測差異表達(dá)unigenes與花粉粒的成熟相關(guān)。通過KEGG數(shù)據(jù)庫比對，3組差異表達(dá)unigenes主要注釋在代謝和遺傳信息處理2大類中。注釋在代謝相關(guān)通路的差異表達(dá)unigenes最多，主要注釋在全局及概要圖、碳水化合物代謝和其他次生代謝產(chǎn)物的生物合成3個(gè)亞類中。相關(guān)研究表明：植物成花過程中需要消耗大量的碳水化合物和蛋白質(zhì)等[14]，因此，注釋在代謝相關(guān)通路的差異表達(dá)unigenes對米櫧雄花成花有重要作用。此外，KEGG代謝通路注釋結(jié)果還顯示：注釋在遺傳信息處理相關(guān)通路的轉(zhuǎn)錄亞類和環(huán)境信息處理相關(guān)通路的信號轉(zhuǎn)錄亞類的差異表達(dá)unigenes較多，表明轉(zhuǎn)錄因子對米櫧雄花的基因表達(dá)有重要作用。

相關(guān)研究表明：轉(zhuǎn)錄因子家族在植物的生殖過程中發(fā)揮著重要作用，其中MADS-box家族是控制花器官形成的重要轉(zhuǎn)錄因子家族[15]。相關(guān)研究根據(jù)系統(tǒng)發(fā)育學(xué)和基因結(jié)構(gòu)的特點(diǎn)，將MADS-box家族分為type-Ⅰ型和type-Ⅱ型[16，17]。通過對米櫧雄花的花芽期、半開期和盛花期的基因表達(dá)量分析，發(fā)現(xiàn)type-Ⅰ型MADS-box家族基因在米櫧成花過程中幾乎不表達(dá)，印證了大多數(shù)植物中以type-Ⅱ型MADS-box家族基因?yàn)橹鞯耐茢郲18]。通過對米櫧type-Ⅱ型MADS-box家族基因進(jìn)行鑒定及表達(dá)量分析，發(fā)現(xiàn)其在米櫧雄花不同發(fā)育階段的表達(dá)具有差異性，其中，鑒定到AP1/FUL類同源基因3個(gè)，AP1基因2個(gè)，其中1個(gè)在花芽期和盛花期表達(dá)量較高，另外1個(gè)在3個(gè)時(shí)期的表達(dá)量均較高；FUL基因1個(gè)，在花芽期和盛花期表達(dá)量較高。研究發(fā)現(xiàn)，F(xiàn)UL基因在花序分生組織中表達(dá)，對花分生組織的起始分化有重要作用，可能參與調(diào)控開花時(shí)間，而AP1基因在花分生組織和萼片的形成中有重要作用[19]。因此，米櫧FUL基因在花芽期和盛花期表達(dá)量較高，推測對米櫧雄花花序分化和開花起重要作用，AP1基因在半開期和盛花期高表達(dá)，調(diào)控花被片的形成。鑒定到AP3/PI類同源基因1個(gè)，為AP3基因，表達(dá)量隨著雄花的發(fā)育逐漸升高。AP3/PI類基因通常在花瓣和雄蕊中表達(dá)，對控制花瓣和雄蕊的形成有重要作用[19]。因此，在米櫧雄花發(fā)育過程中，AP3/PI類同源基因?qū)π廴锖突ū黄陌l(fā)育有重要作用。鑒定到C/D類同源基因6個(gè)，其中米櫧AG基因表達(dá)量隨著雄花的發(fā)育逐漸升高，主要在雄花發(fā)育后期表達(dá)。在歐洲栗(CastaneasativaMill.)和歐洲栓皮櫟的研究中，AG基因主要在雄花發(fā)育后期表達(dá)，與花粉發(fā)育密切相關(guān)[4，10]。在米櫧雄花發(fā)育過程中，AG基因的表達(dá)趨勢與歐洲栗和歐洲栓皮櫟相似，因此，AG類同源基因可能參與調(diào)控米櫧的花粉發(fā)育。米櫧SHP基因在時(shí)間上的表達(dá)模式與AG基因相似，主要在雄花發(fā)育后期表達(dá)。在歐洲栗和歐洲栓皮櫟的研究中，SHP基因在雄花發(fā)育后期高表達(dá)，同時(shí)還在雌花中表達(dá)，表明殼斗科的SHP類同源基因可能是功能性的C類基因[4，10]。而STK基因主要參與調(diào)控胚珠的發(fā)育[20]，在米櫧雄花發(fā)育的3個(gè)時(shí)期均未表達(dá)，印證了米櫧雄花中不存在雌蕊的結(jié)構(gòu)。鑒定到SEP類同源基因4個(gè)，在米櫧雄花發(fā)育的3個(gè)時(shí)期有較高表達(dá)，表明SEP類同源基因在米櫧雄花發(fā)育過程中起協(xié)同作用。

綜上所述，在米櫧雄花發(fā)育過程中，AP1/FUL類、AP3/PI類、C/D類和SEP類同源基因有較高的表達(dá)，且在不同時(shí)期的表達(dá)具有差異性，而這些基因在其他物種中已被證實(shí)參與調(diào)控各個(gè)花器官的發(fā)育。米櫧是雌雄同株的單性花植物，AP3/PI類和C/D類同源基因?qū)π廴锖痛迫锏男纬捎兄匾淖饔?。因此，下一步可重點(diǎn)分析MADS-box家族基因分別在米櫧雄花和雌花不同花器官中的表達(dá)情況，以期闡明這些基因在米櫧性別分化和花器官發(fā)育中的功能和作用機(jī)制。