擬南芥花粉管與柱頭互作的乙醇代謝耦合模型

擬南芥花粉管與柱頭互作的乙醇代謝耦合模型

陶璐，岳訓(xùn)*

(山東農(nóng)業(yè)大學(xué)信息科學(xué)與工程學(xué)院，山東 泰安 271018)

摘要：依據(jù)擬南芥公開的代謝途徑數(shù)據(jù)庫，構(gòu)建了基于酶與酶的擬南芥代謝網(wǎng)絡(luò)模型。利用擬南芥花粉管與柱頭互作過程中的轉(zhuǎn)錄組數(shù)據(jù)，挖掘出花粉管與柱頭在互作過程中的特異表達(dá)基因，進(jìn)一步將特異表達(dá)的酶基因匹配到已構(gòu)建的擬南芥代謝網(wǎng)絡(luò)中，根據(jù)網(wǎng)絡(luò)拓?fù)淠Ｐ椭械墓?jié)點(diǎn)(酶)之間的共表達(dá)關(guān)聯(lián)性，最后給出了一個擬南芥花粉管與柱頭互作的乙醇代謝耦合模型。

關(guān)鍵詞：擬南芥；授粉；代謝網(wǎng)絡(luò)

中圖分類號：Q942文獻(xiàn)標(biāo)志碼：A

收稿日期：2014-11-26;修回日期：2015-01-09.

基金項(xiàng)目：國家社會科學(xué)

作者簡介：種樂熹，男，碩士研究生，研究方向：生物信息學(xué)；E-mail：chong_lx@163.com.

doi:10.3969/j.issn.1672-5565.2015.01.10

Ethanol metabolism coupling model betweenArabidopsispollen tube and stigma

TAO Lu, YUE Xun*

(CollegeofInformationSciencesandEngineering,ShandongAgriculturalUniversity,TaianShandong271018，China)

Abstract：According to the public database of Arabidopsis metabolism pathway, this paper constructed a metabolism correlation networks model based on enzymes correlation. Transcriptional expression data of the interaction between Arabidopsis pollen tub and stigma were used to mine significant differentially expressed genes. The significant differentially expressed genes that encoded enzymes were mapped to the enzyme correlation networks model. According to the co-regulation and inter-regulation model of the co-expressed enzymes in the consecutive steps for metabolic routes, we proposed an ethanol metabolism coupling model between Arabidopsis pollen tube and stigma.

Keywords：Arabidopsis; Pollination; Enzyme correlation networks model

傳粉與受精是植物產(chǎn)生種子的基礎(chǔ)，也是植物進(jìn)行有性生殖的關(guān)鍵所在。植物雌雄配子體互作是植物生殖與發(fā)育生物學(xué)研究領(lǐng)域重要和熱點(diǎn)問題之一[1]。傳粉與受精的效率和質(zhì)量直接關(guān)系到作物產(chǎn)量的高低。系統(tǒng)地闡明植物花粉管與柱頭互作的機(jī)制，對提高作物產(chǎn)量具有重要意義[2]。

花粉管在雌蕊組織內(nèi)的導(dǎo)向性生長是雌雄配子體相互作用的重要階段之一[3]。與植物其他器官不同，花粉管以極性而非擴(kuò)散性方式生長?；ǚ酃芸梢栽隗w外簡單培養(yǎng)基中生長[4]，而當(dāng)其在雌蕊體內(nèi)生長時，雌蕊提供其快速生長所需的物質(zhì)。在過去的20多年中，為了探究植物花粉管與柱頭在授粉過程中的調(diào)控機(jī)理，科學(xué)家們對此進(jìn)行了大量的研究[5-6]。但是，授粉過程中雌雄配子體互作的分子調(diào)控機(jī)制細(xì)節(jié)仍然不是很清楚[7]。

隨著擬南芥(Arabidopsis thaliana)等植物基因組全序列測定的完成，在后基因組學(xué)與蛋白質(zhì)組學(xué)時代，轉(zhuǎn)錄組、蛋白組等組學(xué)技術(shù)的應(yīng)用和從計(jì)算系統(tǒng)生物學(xué)水平上系統(tǒng)深入開展雌雄配子體互作調(diào)控網(wǎng)絡(luò)的研究，已成為植物生殖與發(fā)育生物學(xué)研究領(lǐng)域的重要發(fā)展趨勢。

本文從計(jì)算系統(tǒng)生物學(xué)水平上對擬南芥花粉管在柱頭引導(dǎo)組織中導(dǎo)向生長過程中的互作機(jī)制進(jìn)行探究，依據(jù)擬南芥公開的代謝途徑數(shù)據(jù)庫，構(gòu)建了基于酶與酶的擬南芥代謝網(wǎng)絡(luò)模型。然后利用擬南芥花粉管和柱頭互作過程中的轉(zhuǎn)錄表達(dá)數(shù)據(jù)，挖掘互作過程中的特異表達(dá)的基因，進(jìn)一步將特異表達(dá)基因編碼的酶匹配到已構(gòu)建的代謝互作網(wǎng)絡(luò)中，根據(jù)網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)中的節(jié)點(diǎn)關(guān)聯(lián)性，發(fā)現(xiàn)柱頭組織可以通過產(chǎn)生的乙醇參與到花粉管TCA循環(huán)中，間接為花粉管在柱頭組織中的快速生長提供能量和營養(yǎng)物質(zhì)。最后給出了一個擬南芥花粉管與柱頭互作的乙醇代謝耦合模型。

1材料、數(shù)據(jù)及方法

1.1擬南芥柱頭數(shù)據(jù)

Boavida利用擬南芥ATH1的全基因組陣列, 比較了野生型未授粉柱頭與授粉后0.5H，3.5H和8H柱頭的表達(dá)譜，得到了擬南芥Affymetrix ATH1寡核苷酸微陣列芯片數(shù)據(jù)[8]。本文以這組芯片數(shù)據(jù)為基礎(chǔ)，通過數(shù)據(jù)分析等方法對擬南芥花粉管與柱頭的互作機(jī)制進(jìn)行分析(NCBI基因綜合數(shù)據(jù)庫編號：GSE27281)。

1.2擬南芥花粉管數(shù)據(jù)

采用Qin et al.研究獲得的包括干花粉粒、未萌發(fā)花粉、離體培養(yǎng)0.5H和4H的花粉、萌發(fā)的花粉和在柱頭中生長的花粉管等的芯片數(shù)據(jù)[9](NCBI基因綜合數(shù)據(jù)庫編號：GSE17343)。

1.3擬南芥代謝網(wǎng)絡(luò)的構(gòu)建

構(gòu)建擬南芥代謝網(wǎng)絡(luò)需按照以下七步進(jìn)行：

(1)從擬南芥代謝途徑數(shù)據(jù)庫(http://www. Arabidopsis. org/biocyc/)中下載源數(shù)據(jù)文件(TXT格式)。文件中包含有540個代謝途徑，7 127個酶，3 418個生化反應(yīng)，3 323個化合物等內(nèi)容。

(2)源數(shù)據(jù)文件中的內(nèi)容并不是都對構(gòu)建代謝網(wǎng)絡(luò)有效，所有先將擬南芥代謝途徑數(shù)據(jù)文件(包括代謝途徑、反應(yīng)酶、相關(guān)基因、反應(yīng)式和各種化合物)導(dǎo)入Oracle數(shù)據(jù)庫，做數(shù)據(jù)預(yù)處理。

(3)代謝網(wǎng)絡(luò)模型構(gòu)建規(guī)則的確定。如圖1所示，如果在酶1的催化下反應(yīng)生成了底物A和B，底物A或者底物B又經(jīng)過酶2的催化反應(yīng)生成了底物C，那么酶1和酶2就建立起了關(guān)聯(lián)。代謝物反應(yīng)式是確定酶與酶是否存在聯(lián)系的依據(jù)。在代謝網(wǎng)絡(luò)模型中，網(wǎng)絡(luò)拓?fù)渲械囊粋€節(jié)點(diǎn)就代表一個酶，一條邊代表兩個節(jié)點(diǎn)(酶)之間存在關(guān)聯(lián)。

圖1　代謝網(wǎng)絡(luò)模型構(gòu)建規(guī)則

(4)在確定了代謝網(wǎng)絡(luò)構(gòu)建規(guī)則的基礎(chǔ)上，利用Oracle數(shù)據(jù)庫和SQL(結(jié)構(gòu)化查詢語言)對源文件處理。首先從源數(shù)據(jù)文件中提取出存在的反應(yīng)關(guān)系式以及相對應(yīng)催化反應(yīng)關(guān)系式的酶。然后利用SQL語句，將反應(yīng)關(guān)系式中的反應(yīng)物與催化該反應(yīng)的酶進(jìn)行一一對應(yīng)，同時將反應(yīng)關(guān)系式中的生成物與催化該反應(yīng)的酶進(jìn)行一一對應(yīng)。

(5)剔除上述兩組數(shù)據(jù)中含有小分子(H+, NADH, NADP, NADPH, NH3, ATP, ADP, AMP, NAD, CoA, O2, CO2, Glu 和pyrophosphate)的反應(yīng)物或生成物記錄。由于這些小分子通常是作為細(xì)胞電子轉(zhuǎn)移的媒介，與代謝過程作用不大。如不剔除，會增加所構(gòu)建網(wǎng)絡(luò)的負(fù)擔(dān)，對進(jìn)一步的網(wǎng)絡(luò)拓?fù)浞治鲈斐梢欢ǖ母蓴_。

(6)利用Oracle數(shù)據(jù)庫平臺和SQL語句，在代謝網(wǎng)絡(luò)構(gòu)建規(guī)則的基礎(chǔ)上構(gòu)建酶與酶一一對應(yīng)的關(guān)系。完成構(gòu)建網(wǎng)絡(luò)前的必需數(shù)據(jù)預(yù)處理準(zhǔn)備。

(7)將上一步中已經(jīng)處理準(zhǔn)備好的酶與酶一一對應(yīng)的數(shù)據(jù)導(dǎo)入生物分子網(wǎng)絡(luò)分析軟件Cytoscape。最后，擬南芥代謝網(wǎng)絡(luò)模型在Cytoscape平臺上直觀呈現(xiàn)出來[10]。如圖2(a)所示。

1.4基因本體分析

基因本體(Gene ontology, GO)是一個在生物信息學(xué)領(lǐng)域中廣泛使用的本體。它主要包括三個分支：細(xì)胞組件、分子功能和生物過程[11]。本文使用的生物網(wǎng)絡(luò)基因本體分析工具是軟件Cytoscape中的一個插件(BiNGO 2.3)[12]，用來識別過表達(dá)基因的P值在5×10-5水平上的功能團(tuán)。

1.5差異基因

那些在至少一個生物學(xué)過程中特異表達(dá)，并且野生型和實(shí)驗(yàn)體差異倍數(shù)大于等于2，錯誤反現(xiàn)率小于等于1×10-10的基因，被認(rèn)為是差異基因。利用差異基因進(jìn)行GO分析，可以大致確定差異基因的生物功能。

1.6擬南芥花粉管和柱頭的差異基因編碼酶的子網(wǎng)絡(luò)的構(gòu)建

從擬南芥代謝途徑數(shù)據(jù)庫(http://www.arabidopsis.org/biocyc/)中下載代謝途徑(pathway)的源數(shù)據(jù)文件。文件中包含有酶編號和基因編號等內(nèi)容。將文件導(dǎo)入到Oracle數(shù)據(jù)庫中，利用SQL語句找到擬南芥花粉管和柱頭中差異基因編碼的酶。在擬南芥代謝網(wǎng)絡(luò)模型中將擬南芥花粉管和柱頭在授粉過程中特異表達(dá)基因編碼的酶表示出來。如圖2(b, c)，分別是擬南芥花粉管和柱頭在授粉過程中特意表達(dá)的基因編碼酶在已構(gòu)建的擬南芥代謝網(wǎng)絡(luò)中的標(biāo)記，圖2(d, e)分別是擬南芥花粉管和柱頭在授粉過程中特意表達(dá)的基因編碼的酶所構(gòu)成的子網(wǎng)絡(luò)模型。

圖2　擬南芥花粉管和柱頭代謝網(wǎng)絡(luò)拓?fù)?/p>

注：(a) 擬南芥代謝網(wǎng)絡(luò)模型；(b) 突出顯示的是擬南芥花粉管在授粉過程中特意表達(dá)的基因所編碼的酶在已構(gòu)建的擬南芥代謝網(wǎng)絡(luò)中的標(biāo)記；(c) 突出顯示的是擬南芥柱頭在授粉過程中特意表達(dá)的基因所編碼的酶在已構(gòu)建的擬南芥代謝網(wǎng)絡(luò)中的標(biāo)記；(d) 擬南芥花粉管在授粉過程中特意表達(dá)的基因所編碼的酶所構(gòu)成的子網(wǎng)絡(luò)模型；(e) 擬南芥花粉管在授粉過程中特意表達(dá)的基因所編碼的酶所構(gòu)成的子網(wǎng)絡(luò)模型.

Notes:(a) The metabolic network model of Arabidopsis; (b) The enzymes encoded by differentially expressed genes in pollen tube during the process of pollination were highlighted in the metabolic network model of Arabidopsis; (c) The enzymes encoded by differentially expressed genes in stigma during the process of pollination were highlighted in the metabolic network model of Arabidopsis; (d) The sub-network constructed by the enzymes which encoded by differentially expressed genes in pollen tube; (e) The sub-network constructed by the enzymes which encoded by differentially expressed genes in stigma.

2結(jié)果與分析

2.1擬南芥花粉管和柱頭的差異基因

在擬南芥花粉管中找到1 121個差異基因，在擬南芥柱頭中找到了1 468個差異基因，這些差異基因在授粉過程中的表達(dá)水平都比較顯著。而后，利用生物網(wǎng)絡(luò)基因本體分析工具(BiNGO 2.3)[12]，分析擬南芥花粉管和柱頭的差異基因，發(fā)現(xiàn)大部分的差異基因都富集在幾個分類類別中，如“代謝過程”，“生物合成過程”，“細(xì)胞壁修復(fù)”，“株型細(xì)胞壁松弛”和“細(xì)菌應(yīng)對防御”。

2.3擬南芥花粉管中乙醇二級代謝和TCA代謝途徑上酶基因的特異表達(dá)

通過對擬南芥花粉管和柱頭數(shù)據(jù)的分析，發(fā)現(xiàn)了一些編碼催化連續(xù)反應(yīng)的共表達(dá)酶的特意表達(dá)基因。

擬南芥授粉花粉管中乙醇二級代謝(如圖3)：(1)AT5G63620 編碼乙醇脫氫酶(alcohol dehydrogenase)(EC:1.1.1.1)，催化乙醇(ethanol)轉(zhuǎn)化為乙醛(acetaldehyde)；(2)AT1G54100 編碼乙醛脫氫酶(aldehyde dehydrogenase)(EC:1.2.1.3)，催化乙醛(acetaldehyde)轉(zhuǎn)化為醋酸鹽(acetate)；(3)AT5G36880 編碼乙酰輔酶A合成酶(acetyl-CoA synthetase)(EC:6.2.1.1)，催化醋酸鹽(acetate)轉(zhuǎn)化為乙酰輔酶( acetyl-CoA)。

擬南芥授粉花粉管中TCA循環(huán)代謝(如圖4)：(1)AT2G17130 編碼異檸檬酸脫氫酶(isocitrate dehydrogenase)(NAD+)(EC:1.1.1.41)，催化檸檬酸脫氫酶(D-threo-isocitrate)轉(zhuǎn)化為2-氧戊二酸鹽(2-oxoglutarate)；(2)AT5G65750和AT3G55410編碼2-氧戊二酸鹽脫氫酶E1 (2-oxoglutarate dehydrogenase E1)(EC 1.2.4.2)，催化2-氧戊二酸鹽(2-oxoglutarate)轉(zhuǎn)化為琥珀酰輔酶A(succinyl-CoA)；(3)AT5G08300 編碼琥珀酰輔酶A合成酶 (Succinyl-CoA synthetase)(EC 6.2.1.5)，催化琥珀酰輔酶A (Succinyl-CoA)轉(zhuǎn)化為琥珀酸鹽(succinate)；(4)AT5G66760編碼琥珀酸脫氫酶(succinate dehydrogenase (ubiquinone) flavoprotein subunit)(EC 1.3.5.1)，催化琥珀酸鹽(Succinate)轉(zhuǎn)化為延胡索酸鹽(Fumarate)。

如上數(shù)據(jù)都表明：在擬南芥花粉管中，乙醇通過二級代謝可參與到TCA循環(huán)代謝途徑中，促進(jìn)花粉管中ATP產(chǎn)量的提高，為花粉管在柱頭引導(dǎo)組織中極性、快速生長提供能量和物質(zhì)保障。

2.4擬南芥柱頭中糖酵解代謝途徑上酶基因的特異表達(dá)

花粉管是目前所知的生長的最快的植物細(xì)胞。而這種快速的生長模式，就需要大量的能量和物質(zhì)做基礎(chǔ)[13]。為了滿足快速生長花粉管的細(xì)胞需求，柱頭中大量的蛋白質(zhì)合成、生化反應(yīng)等代謝途徑被激活，為花粉管的快速生長提供充足的物質(zhì)和能量[14]。利用多組學(xué)網(wǎng)絡(luò)分析的方法同樣發(fā)現(xiàn)，在擬南芥授粉的柱頭中，所編碼共表達(dá)酶的基因、糖酵解的連續(xù)反應(yīng)以及與能產(chǎn)生高能營養(yǎng)(乙醇)的關(guān)鍵酶的轉(zhuǎn)錄表達(dá)水平都呈現(xiàn)出明顯升高趨勢。如圖5所示。

圖3　擬南芥花粉管中乙醇二級代謝

圖4　擬南芥花粉管中 TCA循環(huán)代謝途徑

(1)兩個擬南芥基因(AT1G20950和AT4G04040)編碼6-磷酸果糖轉(zhuǎn)磷酸酶(fructose-6-phosphate 1-phosphotransferase) (EC 2.7.1.90)。催化6-磷酸果糖轉(zhuǎn)磷酸酶(fructose-6-phosphate 1-phosphotransferase)轉(zhuǎn)化為果糖(fructose-1)和磷酸(6-bisphosphate)，作為糖酵解(glycolysis IV)代謝途徑的最開始反應(yīng)。

(2)在糖酵解代謝的第二步中，多組學(xué)網(wǎng)絡(luò)分析表明在擬南芥中沒有編碼果糖二磷酸醛縮酶(fructose-bisphosphate aldolase) (EC 4.1.2.13)的差異基因，果糖二磷酸醛縮酶可以催化果糖(fructose-1)、二磷酸(fructose-1)轉(zhuǎn)化為3-磷酸甘油醛(D-glyceraldehyde 3-phosphate)。本文推斷果糖二磷酸醛縮酶(fructose-bisphosphate aldolase)可能是起因于被稱為“自噬行為”的細(xì)胞降解過程。

(3)兩個擬南芥基因(AT1G16300 和 AT1G79530)編碼甘油醛-3-磷酸脫氫酶(glyceraldehyde-3-phosphate dehydrogenase) (GAPDH) (EC 1.2.1.12)，催化D-3-磷酸甘油醛(D-glyceraldehyde 3-phosphate)轉(zhuǎn)化為1,3-二磷酸-D-甘油酸(1,3-bisphospho-D-glycerate)。

(4)磷酸甘油酸激酶(Phosphoglycerate kinase) (EC 2.7.2.3)，催化1,3-二磷酸-D-甘油酸(1,3-bisphospho-D-glycerate)轉(zhuǎn)化為3-磷酸-D-甘油酸(3-phospho-D-glycerate)。重要的是，擬南芥的柱頭中并沒有發(fā)現(xiàn)編碼磷酸甘油酸激酶的基因。而后，通過重新查找這些基因的源數(shù)據(jù)庫，結(jié)果與分析發(fā)現(xiàn)結(jié)論一致。

(5)AT1G09780編碼2,3-二磷酸甘油酸依賴性磷酸甘油酸變位酶(2,3-bisphosphoglycerate-dependent phosphoglycerate mutase) (EC 5.4.2.1)，催化3-磷酸-D-甘油酸(3-phospho-D-glycerate)轉(zhuǎn)化為2-磷酸-D-甘油酸(2-phospho-D-glycerate)。

(6)AT1G74030編碼烯醇酶(enolase) (EC 4.2.1.11)，催化2-磷酸-D-甘油酸(2-phospho-D-glycerate)轉(zhuǎn)化為磷酸烯醇丙酮酸鹽(phosphoenolpyruvate)。

(7)丙酮酸激酶(pyruvate kinase) (PKM2) (EC 2.7.1.40)，催化磷酸烯醇丙酮酸鹽(phosphoenolpyruvate)轉(zhuǎn)化為丙酮酸鹽(pyruvate)，是糖酵解反應(yīng)最后的一個步驟。

有趣的是，在分析過程中發(fā)現(xiàn)了幾個能產(chǎn)生高能營養(yǎng)物質(zhì)(乙醇)的關(guān)鍵酶。首先三個擬南芥基因(AT1G04040, AT3G17790 和 AT3G52820)編碼了酸性磷酸酶(acid phosphatase) (EC 3.1.3.2)催化磷酸單酯(phosphate monoester)轉(zhuǎn)化為乙醇(ethanol)。其次，兩個擬南芥基因編碼了甘油磷酸二脂酶(glycerophosphodiester phosphodiesterase) (EC 3.1.4.46)，催化甘油磷酸二脂酶(glycerophosphodiester)轉(zhuǎn)化為乙醇(ethanol)。

圖5　擬南芥柱頭中的糖酵解代謝

3討論

通過上面的數(shù)據(jù)分析表明：當(dāng)授粉后的花粉管在柱頭分泌區(qū)的細(xì)胞外基質(zhì)中生長時，柱頭可能確實(shí)能被激活許多代謝活動，為花粉管生長提供能量和所需物質(zhì)。例如，糖酵解相關(guān)的代謝途徑被激活：糖分轉(zhuǎn)化為丙酮酸，丙酮酸被進(jìn)一步代謝為高能營養(yǎng)物質(zhì)乙醇，柱頭組織中的乙醇又被轉(zhuǎn)運(yùn)到極性、快速生長的花粉管中，參與到花粉管的TCA循環(huán)代謝途徑，產(chǎn)生能量和營養(yǎng)物質(zhì)。

根據(jù)上述分析，本文提出了擬南芥花粉管與柱頭互作的乙醇耦合模型(如圖6)。

圖6　擬南芥花粉管與柱頭互作的乙醇代謝耦合模型

花粉萌發(fā)和花粉管的快速生長是一個需要大量能量和物質(zhì)的生物過程。特別是某些代謝途徑，如細(xì)胞質(zhì)溶質(zhì)的糖酵解代謝等，這些代謝過程中產(chǎn)生的丙酮酸脫氫酶、TCA(檸檬酸)循環(huán)、乙醇等都對花粉管的生長有著非常重要的作用。

(1)當(dāng)花粉管在柱頭的引導(dǎo)組織中生長延長時，花粉管的延長觸發(fā)了柱頭引導(dǎo)組織中糖酵解產(chǎn)生的能量的轉(zhuǎn)移，激活了糖類物質(zhì)轉(zhuǎn)化為丙酮酸鹽。

(2)丙酮酸鹽進(jìn)一步代謝為能被花粉管吸收利用的高能營養(yǎng)物質(zhì)乙醇。

(3)在花粉管中乙醇通過II級降解，進(jìn)一步代謝為乙酰輔酶A(acetyl-CoA)。

(4)乙酰輔酶A(acetyl-CoA)參與到花粉管的TCA循環(huán)匯總，促進(jìn)了ATP產(chǎn)量的增加，從而為花粉管的快速生長提供了能量保證。

綜上，在擬南芥花粉管和柱頭的互作生物學(xué)問題背景下，從計(jì)算系統(tǒng)生物學(xué)水平上，以擬南芥花粉管和柱頭代謝組學(xué)數(shù)據(jù)為基礎(chǔ)，運(yùn)用多組學(xué)網(wǎng)絡(luò)分析方法，構(gòu)建了擬南芥代謝互作網(wǎng)絡(luò)模型，進(jìn)一步挖掘擬南芥花粉管和柱頭在授粉過程中特異表達(dá)的基因，根據(jù)網(wǎng)絡(luò)拓?fù)渲械墓?jié)點(diǎn)關(guān)聯(lián)性和代謝互作網(wǎng)絡(luò)的實(shí)際生物學(xué)問題，發(fā)現(xiàn)擬南芥花粉管中的乙醇和柱頭中糖酵解代謝產(chǎn)生的被轉(zhuǎn)運(yùn)到花粉管中的乙醇可以參與到花粉管TCA循環(huán)代謝途徑中，為花粉管在柱頭引導(dǎo)組織中快速、極性生長所需要的能量和營養(yǎng)物質(zhì)提供了保障。最后，本文建立了擬南芥花粉管和柱頭互作的乙醇代謝耦合模型。

參考文獻(xiàn)(References)

[1]HISCOCK S J, ALLEN A M. Diverse cell signaling pathways regulate pollen-stigma interactions: the search for consensus [J]. New Phytologist, 2008, 179(2):286-317.

[2]SIVITZ A B,REINDERS A,WARD J M, et al. Arabidopsis sucrose transporter AtSUC1 is important for pollen germination and sucrose induced anthocyanin accumulation [J]. Plant Physiology, 2008, 147(1): 92-100.

[3]劉敬婧,侯英楠,瞿禮嘉.高等植物受精過程中雌-雄相互作用的分子調(diào)控機(jī)制[J].中國科學(xué)：生命科學(xué),2013,43(10):842-853.

LIU Jingjing, HOU Yingnan, QU Lijia. Molecular mechanisms behind the male-female interaction in fertilization in higher plants [J]. Scientia Sinica Vitae, 2013, 43(10):842-853.

[4]WANG Y X, MAGNARD J L, MCCORMICK S, et al. Progression through meiosis I and meiosis II in arabidopsis anthers is regulated by an A-type cyclin predominately expressed in prophase I [J].Plant Physiology, 2004, 136(4): 4127-4135.

[5]SWANSON R, EDLUND A F, PREUSS D. Species specificity in pollen-pistil interactions[J]. Annual Review of Genetics, 2004, 38: 793-818.

[6]GAO X Q, ZHU D Z, ZHANG X S. Stigma factors regulating self-compatible pollination[J].Frontiers in Biology, 2010, 5(2): 156-163.

[7]HIGASHIYAMA T. Peptide signaling in pollen-pistil interactions [J]. Plant Cell Physiology, 2010, 51(2): 177-189.

[8]BOAVIDA L C, BORGES F, BECKER J D, et al. Whole genome analysis of gene expression reveals coordinated activation of signaling and metabolic pathways during pollen-pistil interactions in arabidopsis [J]. Plant Physiology, 2011, 155(4): 2066-2080.

[9]QIN Y, LEYDON A R, MANZIELLO A, et al. Penetration of the stigma and style elicits a novel transcriptome in pollen tubes, pointing to genes critical for growth in a pistil [J]. Plos Genetics, 2009, 5(8): e1000621.

[10]LI M N, XU W Y, YANG W Q, et al. Genome-wide gene expression profiling reveals conserved and novel molecular functions of the stigma of rice [J]. Plant Physiology, 2007, 144(4):1797-1812.

[11]YOUNG M D, WAKEFIELD M J, SMYTH G K, et al. Gene ontology analysis for RNA-SEQ: accounting for selection bias [J]. Genome Biology, 2010, 11(2):R14.

[12]MAERE S, HEYMANS K, KUIPER M. Bingo: a cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks [J]. Bioinformatics, 2005, 21(16):3448-3449.

[13]TAYLOR L P, HEPLER P K. Pollen germination and tube growth[J]. Annual Review of Plant Physiology and Plant Molecular Biology, 1997, 48:461-491.

[14]OBERMEYER G, FRAGNER L, LANG V,et al. Dynamic adaption of metabolic pathways during germination and growth of lily pollen tubes after inhibition of the electron transport chain [J]. Plant Physiology, 2013, 162(4):1822-1833.

*通信作者：胡德華，男，教授，博士生導(dǎo)師，研究方向：生物信息學(xué);E-mail：hudehua2000@163.com.