農(nóng)作物病原卵菌基因組數(shù)據(jù)庫(kù)資源概述

高亞梅， 王相晶， 向文勝＊

（1.東北農(nóng)業(yè)大學(xué)生命科學(xué)學(xué)院，哈爾濱 150030；2.黑龍江八一農(nóng)墾大學(xué)生命科學(xué)技術(shù)學(xué)院，大慶 163319）

卵菌（Oomycetes）屬于色藻界（Chromista），包括腐生菌及植物、動(dòng)物和許多微生物的病原菌。病原卵菌導(dǎo)致許多作物、花卉等發(fā)生災(zāi)難性病害，例如引起大豆疫霉根腐病的大豆疫霉（Phytophthora sojae Kaufmann et Gerdemann）；引起的馬鈴薯晚疫病的 致 病 疫 霉 ［P.infestans （Montagne）de Bary］［1－2］。在進(jìn)化地位上，卵菌與褐藻和硅藻有較近的親緣關(guān)系，形成一個(gè)獨(dú)特的二倍體微生物類群［3］，其分子遺傳學(xué)的研究相對(duì)落后。隨著分子生物學(xué)技術(shù)的發(fā)展，卵菌基因組學(xué)及功能基因組學(xué)的研究逐漸受到研究者的重視，迅速積累了大量基因組數(shù)據(jù)資源，這些資源又成為卵菌功能基因組學(xué)研究的基礎(chǔ)之一，同時(shí)也是基因組數(shù)據(jù)挖掘的重要來(lái)源，加快了卵菌分子生物學(xué)的各項(xiàng)研究。本文對(duì)農(nóng)作物病原卵菌基因組數(shù)據(jù)庫(kù)資源及其檢索方法做一概述。

1 卵菌結(jié)構(gòu)基因組學(xué)數(shù)據(jù)庫(kù)資源及檢索方法

隨著測(cè)序技術(shù)的不斷發(fā)展，高通量測(cè)序使得測(cè)序成本不斷降低，越來(lái)越多的基因組計(jì)劃啟動(dòng)，甚至可以用宏基因組學(xué)（metagenomics）技術(shù)測(cè)定一個(gè)生物群體的序列。因此，在各類數(shù)據(jù)庫(kù)中積累了海量的生物序列信息，供研究者檢索、下載和分析?；蚪M序列及結(jié)構(gòu)、基因功能和調(diào)控元件的相關(guān)注釋可為深入研究病原菌毒力、致病性、寄主特異性調(diào)控、病原菌生化與生理機(jī)制、生態(tài)位適應(yīng)性等提供大量有用信息。目前，共有7個(gè)卵菌物種［大豆疫霉P.sojae、分枝疫霉P.ramorumWerres，De Cock et Man、致病疫霉 P.infestans、辣椒疫霉P.capsici Leonian、終極腐霉 Pythium ultimum Trow、Hyaloperonospora parasitica （G?um.）G?ker，Riethm.，Voglmayr，Weiss ＆ Oberw.和古巴假霜霉 Pseudoperonospora cubensis （Berkeley et Curtis）Rostovzev］的全基因組序列正式對(duì)外公布（詳見表1）［4－9］；另外2個(gè)卵菌物種（Peronospora parasitica Tul.及Saprolegnia parasitica Coker）的基因組測(cè)序工作正在進(jìn)行中。通過(guò)對(duì)這些基因組序列和結(jié)構(gòu)的分析揭示了卵菌基因組在大小、致病相關(guān)基因、專性活體寄生等方面存在多樣性［10］，并為病原菌遺傳與致病機(jī)理、宿主－病原菌相互作用的特征（宿主的特異性，毒性策略等）以及宿主－病原菌相互作用的進(jìn)化等各方面研究提供了大量的信息?；蚪M分析揭示的卵菌基因組進(jìn)化特征包括重復(fù)序列在基因組內(nèi)的擴(kuò)增、缺失、基因融合和基因垂直轉(zhuǎn)移等［11－13］。這些無(wú)疑都加深了我們對(duì)于卵菌基因組結(jié)構(gòu)及進(jìn)化的認(rèn)識(shí)。

表1 已公布的測(cè)序卵菌基因組信息資源Table 1 Oomycete sequenced genome database resources

聯(lián)合基因組研究所JGI網(wǎng)站（DOE Joint Genome Institute）提供細(xì)菌、古細(xì)菌、真核生物和宏基因組等眾多物種的基因組序列信息和分析服務(wù)，在其主頁(yè)的搜索工具欄選擇相應(yīng)的物種即可進(jìn)入相應(yīng)物種的數(shù)據(jù)庫(kù)，其中卵菌物種包括P.sojae、P.ramorum和P.cinnamomi［14］。基因組序列信息可通過(guò)該數(shù)據(jù)庫(kù)SEARCH功能進(jìn)行關(guān)鍵詞檢索，利用BLAST進(jìn)行序列相似性搜索，利用DOWNLOAD下載，利用BROWSE以圖形化方式瀏覽（見圖1）。另外 GO、KEGG、KOG、CLUSTER、SYNTENY提供基因功能、代謝圖譜、共線性、簇等注釋信息的查詢檢索。

圖1 JGI網(wǎng)站P.cinnamomi var.cinnamomi基因組瀏覽界面Fig.1 The genome browse interface of P.cinnamomi var.cinnamomi in JGI

VMD數(shù)據(jù)庫(kù)（VBI Microbial Database）［15］是另一重要的卵菌基因組數(shù)據(jù)庫(kù)，包括Pythium ultimum（V1.0），Phytophthora infestans（V4.0），P.capsici（V11.0），P.sojae（V1.0，V4.0，V5.0），P.ramorum（V1.0），H.arabidopsidis（V3.0，V6.0，V8.3）基因組及注釋信息。每個(gè)數(shù)據(jù)庫(kù)均提供圖形化的瀏覽方式，顯示每個(gè)基因的位置、模式和相關(guān)基因組的BLAST信息等。VMD數(shù)據(jù)庫(kù)還提供檢索、下載、注釋和工具箱等服務(wù)。檢索可以通過(guò)不同方式進(jìn)行，包括 Gene ID、GO ID、Scaffold Number、protein domains／motifs／functions、primary annotation、affymetrixID 等。檢索結(jié)果頁(yè)面提供基因的序列、基因組位置、功能注釋及軟件預(yù)測(cè)結(jié)果等信息。BLAST工具提供多個(gè)可選數(shù)據(jù)庫(kù)，包括已測(cè)序卵菌基因組、表達(dá)數(shù)據(jù)庫(kù)、大豆和擬南芥數(shù)據(jù)庫(kù)等，進(jìn)行有針對(duì)性的序列相似性比對(duì)分析。P.infestans T30－4基因組序列可以在Phytophthora infestans Database數(shù)據(jù)庫(kù)下載，該數(shù)據(jù)庫(kù)同樣提供各類分析工具和可視化的瀏覽頁(yè)面。

Pythium Genome Database是專門存儲(chǔ)Pythium ultimum序列信息的數(shù)據(jù)庫(kù)，包括其代表性菌株P(guān)ythium ultimum DAOM BR144的基因組序列和不同營(yíng)養(yǎng)條件下的菌絲體的EST序列。在Pythium Genome Database可以下載 Genome Assembly、Transcripts、Proteins和Gene Model GFF3的信息。通過(guò)Genome Browser可視化瀏覽基因組，可以使用序列名、基因名、遺傳位點(diǎn)或區(qū)域標(biāo)記進(jìn)行檢索。BLAST工具可以選擇卵菌基因組、表達(dá)數(shù)據(jù)等12個(gè)數(shù)據(jù)庫(kù)進(jìn)行序列比對(duì)。該數(shù)據(jù)庫(kù)還將儲(chǔ)存Pythium ultimum另一菌株高通量測(cè)序結(jié)果，開展對(duì)兩菌株及它們與其他已測(cè)序卵菌的基因組比較研究。

卵菌基因組測(cè)序多采用第二代測(cè)序技術(shù)，隨著測(cè)序成本的下降，將會(huì)有更多的卵菌基因組數(shù)據(jù)產(chǎn)生，包括致病菌和非致病菌，為功能基因組學(xué)和比較基因組學(xué)研究提供更多的數(shù)據(jù)。

2 卵菌功能基因組學(xué)數(shù)據(jù)庫(kù)資源及檢索方法

目前，在各種數(shù)據(jù)庫(kù)中均存儲(chǔ)了大量的各種類型、代表不同信息的數(shù)據(jù)，如何對(duì)這些數(shù)據(jù)進(jìn)行分析，獲得蘊(yùn)含其中的生物學(xué)意義是當(dāng)前生物信息學(xué)的重要研究領(lǐng)域。利用生物信息學(xué)方法對(duì)卵菌基因組內(nèi)大量未知功能基因的分析比較是基因功能研究和數(shù)據(jù)注釋的重要手段。以下介紹的一些常用的卵菌EST數(shù)據(jù)庫(kù)、綜合性功能基因組資源平臺(tái)和專門數(shù)據(jù)庫(kù)是開展卵菌功能基因組學(xué)研究的重要資源。

2.1 卵菌EST數(shù)據(jù)庫(kù)資源

卵菌表達(dá)序列標(biāo)簽EST序列（Expressed Sequence Tags，ESTs）的大量公布為進(jìn)行全面的基因結(jié)構(gòu)和功能研究提供了重要的數(shù)據(jù)，必將促進(jìn)對(duì)卵菌致病分子機(jī)制的認(rèn)識(shí)，從而為防治卵菌病害設(shè)計(jì)合理有效的策略。目前數(shù)據(jù)庫(kù)中的EST序列集中在已正式公布的基因組序列的卵菌物種中，例如已有大豆疫霉（P.sojae）的 EST約33350條，而致病疫霉（P.infestans）有99320條（http：∥phytophthora.vbi.vt.edu／EST），還包括一些正在或已完成測(cè)序的卵菌物種P.capsici、Phytophthora brassicae De Cock et al、Albugo candida （Pers.）Roussel等［16］。這些 EST序列大部分為病原卵菌侵染寄主植物過(guò)程中的轉(zhuǎn)錄組，代表疫霉侵染與致病過(guò)程中表達(dá)的基因，還包括卵菌不同發(fā)育階段及不同環(huán)境條件下的EST。卵菌EST存放于多個(gè)數(shù)據(jù)庫(kù)中，如基因組的數(shù)據(jù)庫(kù)、NCBI、轉(zhuǎn)錄組數(shù)據(jù)庫(kù)OTD等，提供多種途徑可以獲取和檢索。

OTD卵菌轉(zhuǎn)錄組學(xué)數(shù)據(jù)庫(kù)（Oomycetes Transcriptomics Database，V4.0，http：∥vmd.vbi.vt.edu／transcripts／）是一個(gè)整合的轉(zhuǎn)錄組和EST數(shù)據(jù)資源的數(shù)據(jù)庫(kù)［17］，儲(chǔ)存卵菌侵染寄主和不同生長(zhǎng)條件下的轉(zhuǎn)錄組數(shù)據(jù)，包括EST序列、NGS（next generation sequencing）轉(zhuǎn)錄組序列、數(shù)字基因表達(dá)譜序列 （Digital Gene Expression Profile，DGE）等。OTD瀏覽器以可視化方式顯示轉(zhuǎn)錄本的組裝，轉(zhuǎn)錄本在基因組上的定位、表達(dá)譜及轉(zhuǎn)錄本在基因組特定位置覆蓋度等信息。該數(shù)據(jù)庫(kù)提供多種檢索方式，可利用EST序列和重疊群ID號(hào)、關(guān)鍵詞、表達(dá)量等進(jìn)行檢索。來(lái)自H.arabidopsidis的37492 ESTs通過(guò)聚類形成Unigenes，并利用BLAST、InterProScan、TMHMM、TargetP和SignalP進(jìn)行了注釋。EST序列可從該數(shù)據(jù)庫(kù)下載。OTD提供與VMD、PTD數(shù)據(jù)庫(kù)的鏈接界面。疫霉轉(zhuǎn)錄組學(xué)數(shù)據(jù)庫(kù)（Phytophthora transcriptomics database，PTD http：∥phy.njau.edu.cn／ptd）提供大豆疫霉的數(shù)字基因表達(dá)譜信息［18］，可通過(guò)基因ID、關(guān)鍵詞進(jìn)行檢索，利用序列比對(duì)BLAST搜索同源基因，可查看不同發(fā)育階段的差異表達(dá)基因。基于轉(zhuǎn)錄信息的基因研究是功能基因組學(xué)研究的重要策略，例如EST數(shù)據(jù)不僅為基因的功能研究提供有價(jià)值的表達(dá)信息，而且基于EST的電子基因克隆技術(shù)也可用來(lái)發(fā)現(xiàn)新基因。近年來(lái)，二代測(cè)序技術(shù)產(chǎn)生的轉(zhuǎn)錄組序列為基因功能研究提供了更為有利的工具，例如RNA－Seq產(chǎn)生的轉(zhuǎn)錄本信息可以更為有效地發(fā)現(xiàn)未知轉(zhuǎn)錄本和稀有轉(zhuǎn)錄本，提供轉(zhuǎn)錄本的結(jié)構(gòu)和表達(dá)水平信息，用于精確地識(shí)別可變剪切位點(diǎn)以及cSNP（編碼序列單核苷酸多態(tài)性）［19］。

2.2 卵菌綜合性功能基因組資源

綜合性疫霉病原菌基因組學(xué)資源（Comprehensive Phytopathogen Genomics Resource，CPGR）包括4個(gè)數(shù)據(jù)庫(kù)：基因組倉(cāng)庫(kù)、注釋數(shù)據(jù)庫(kù)、轉(zhuǎn)錄物組裝和rDNA數(shù)據(jù)庫(kù)［20］。基因組倉(cāng)庫(kù)儲(chǔ)存包括植物病原菌基因組及EST測(cè)序計(jì)劃信息，每年至少更新2次。注釋數(shù)據(jù)庫(kù)包含基因組注釋信息，包括物種名、所屬類群、引起的疾病、NCBI序列登錄號(hào)、基因組大小、測(cè)序計(jì)劃狀態(tài)及機(jī)構(gòu)、特征、相關(guān)注釋以及文獻(xiàn)信息等。數(shù)據(jù)庫(kù)頁(yè)面包含依據(jù)分類單元和測(cè)序狀態(tài)的過(guò)濾器功能。轉(zhuǎn)錄物組裝數(shù)據(jù)來(lái)自NCBI數(shù)據(jù)庫(kù)的EST和mRNA，去除低重復(fù)和污染序列后，經(jīng)拼接注釋而成。CPGR的轉(zhuǎn)錄物組裝和來(lái)自GenBank的序列可在其FTP下載。rDNA數(shù)據(jù)庫(kù)包括植物病原菌的名稱和分類單元標(biāo)識(shí)符，另外該數(shù)據(jù)庫(kù)在病原菌和植物病害之間做了超鏈接。CPGR的搜索工具允許用戶使用序列和功能注釋搜索該數(shù)據(jù)庫(kù)。卵菌注釋可以通過(guò)位點(diǎn)標(biāo)識(shí)符、功能、Pfam功能域、Interpro功能域搜索。對(duì)于已測(cè)序的卵菌，用戶可在CPGR網(wǎng)站利用其提供的一系列界面實(shí)現(xiàn)對(duì)基因組信息（序列和注釋信息）的瀏覽、下載和檢索。目前只有P.infestans T30－4和Pythium ultimum DAOM BR144可利 用 CPGR 的 genome browser進(jìn)行可視化瀏覽，包括等位基因、基因模式（gene models）、rRNA 和tRNA 基因、推測(cè)的簡(jiǎn)單重復(fù)序列SSRs（putative simple sequence repeats）、GC 含 量 （GC content）和 六 框 翻 譯 （sixframe translation）。位點(diǎn)信息和基因模式與其他內(nèi)容形成超鏈接?；蚪M瀏覽器還包括基因列表（gene list）、RNA 基因列表（RNA gene list）、Pfam功能域列表、Interpro功能域列表、推測(cè)SSRs，縮放工具允許用戶在各個(gè)分辨率水平瀏覽基因組。主頁(yè)的resource項(xiàng)下包含的簡(jiǎn)單重復(fù)序列標(biāo)記搜索工具（SSR Candidate Marker Search Tool）可以用于對(duì)提交的序列進(jìn)行SSR搜索，并利用PRIMER3搜索合適的擴(kuò)增引物。用戶可自行設(shè)定SSR類型和重復(fù)數(shù)目。SSRs廣泛存在于各類生物基因組中，是應(yīng)用較為廣泛的遺傳標(biāo)記［21－22］。在基因組序列水平上對(duì)SSR的系統(tǒng)分析是SSR分子標(biāo)記的基礎(chǔ)性工作［23－24］。我 們 利 用 CPGR 數(shù) 據(jù) 庫(kù) 的 SSR 數(shù) 據(jù) 對(duì)P.infestans與Pythium ultimum 基因組內(nèi)的SSR進(jìn)行了比較分析。P.infestans包含1958個(gè)SSR，平均每122.5kb有一個(gè)SSR，Pythium ultimum 包含1911個(gè)SSR，平均每22.39kb有一個(gè)SSR，兩基因組內(nèi)均以二、三、六堿基重復(fù)為最豐富，在基序的重復(fù)次數(shù)上也具有一定的相似性。在二、三、六堿基重復(fù)中，每種都存在優(yōu)勢(shì)基序，詳見表2。在P.infestans基因組中，二、三、六堿基重復(fù)分別占總SSR的26.3%，21.81%和35.44%，在 Pythium ultimum中它們所占的百分比分別為27.68%，26.84%和27.42%，其他堿基重復(fù)所占的比例很少。在兩基因組二、三堿基重復(fù)的優(yōu)勢(shì)基序中，大部分是相同的優(yōu)勢(shì)基序類型，但是六堿基重復(fù)的優(yōu)勢(shì)基序相同的比較少，而且，六堿基重復(fù)的基序類型分布也比較分散。但是兩基因組的重復(fù)序列的長(zhǎng)度都在20bp左右，長(zhǎng)序列較少。而且有些基序類型沒有出現(xiàn)。由以上分析可以看出，兩個(gè)卵菌基因組內(nèi)的SSR數(shù)目相差不多，但由于基因組大小差異導(dǎo)致SSR在基因組內(nèi)的密度存在差異。優(yōu)勢(shì)堿基重復(fù)單元集中在二、三、六，說(shuō)明其遺傳變異速度中等。SSR長(zhǎng)度比較集中?；蚪M水平上的SSR分析為在以后的研究中應(yīng)用SSR標(biāo)記提供了有用的信息。由于SSR標(biāo)記在已測(cè)序卵菌基因組內(nèi)的位置及側(cè)翼序列是已知的，所以可以快速、準(zhǔn)確地找到連鎖的基因位點(diǎn)，從而完成候選基因的篩查、鑒定、基因注釋及功能研究，SSR標(biāo)記將在卵菌功能基因組研究中發(fā)揮重要作用。

CFGP（Comparative Fungal Genomics Platform）是一個(gè)綜合性的真菌比較基因組學(xué)平臺(tái)，包含283個(gè)真菌基因組的數(shù)據(jù)信息（6個(gè)卵菌）、專門數(shù)據(jù)庫(kù)的超鏈接和各類比較基因組分析工具，包括BLAST、ClustalW、InterProScan、SignalP等［25］。平臺(tái)采用數(shù)據(jù)驅(qū)動(dòng)的用戶界面（Data－driven User Interface DUI），用戶可實(shí)現(xiàn)數(shù)據(jù)瀏覽，選擇和分析的連續(xù)進(jìn)行，節(jié)省收集數(shù)據(jù)、轉(zhuǎn)化數(shù)據(jù)格式、提交數(shù)據(jù)、不同數(shù)據(jù)庫(kù)間瀏覽查詢的時(shí)間，極大地方便了大量數(shù)據(jù)的分析。該平臺(tái)的數(shù)據(jù)查詢主要是通過(guò)頁(yè)面的SEQUENCE菜單實(shí)現(xiàn)，有10種查詢方式：Taxonomy、Contig、Genome reference、MyGene、SWISSprot、MSIPI、NR、PDB、PDBchain、SequenceSet browser。數(shù)據(jù)通過(guò) SNUGB（Seoul National University Genome Browser）以可視化方式顯示，并可以選擇顯示不同的數(shù)據(jù)庫(kù)信息。查詢到的序列通過(guò)勾選進(jìn)入到Favorites功能框架。Favorites作為生物信息分析的工作界面，包括Edit、Function、Anal－ysis和Download 4項(xiàng)功能，每個(gè)功能分若干項(xiàng)（見圖2）。其中的Function包括BLAST、BLASTmatrix等27種分析工具，是主要的分析界面，可用于序列比對(duì)分析、功能結(jié)構(gòu)域分析查詢、系統(tǒng)進(jìn)化分析、分泌蛋白分析、亞細(xì)胞定位、跨膜螺旋預(yù)測(cè)、RNA二級(jí)結(jié)構(gòu)預(yù)測(cè)、翻譯后修飾和保守結(jié)構(gòu)域搜索。選擇一種工具后，CFGP即通過(guò)PERL模塊實(shí)現(xiàn)與相應(yīng)程序的鏈接，設(shè)定參數(shù)即可運(yùn)行獲得結(jié)果。分析結(jié)果保存在History中。除了查看自己的分析結(jié)果，在History中也可查看其他人在CFGP分析的結(jié)果。在分析工具中，值得一提的是BLASTMatrix工具，它用于在多個(gè)不同真菌物種中搜索潛在的同源基因。BLASTMatrix結(jié)果提供提交序列在每個(gè)物種最佳命中序列的列表，該列表以每個(gè)物種的分類地位組織在一起，同時(shí)會(huì)提供在不同分類群組內(nèi)同源基因的分布模式圖。輸出結(jié)果還包括InterPro或GO術(shù)語(yǔ)，可以幫助用戶預(yù)測(cè)可能的基因功能，進(jìn)一步分析確定其直系同源關(guān)系。

表2 P.infestans與Pythium ultimum基因組內(nèi)SSR二、三、六堿基重復(fù)優(yōu)勢(shì)基序Tabel 2 Major motifs of di－，tri－and hexa－nucleotide repeats in genome of P.infestans and Pythium ultimum

圖2 CFGP網(wǎng)站的Favorites工作界面Fig.2 The favorites interface in CFGP

FungiDB真菌功能基因組學(xué)資源，包括33種真菌和6種卵菌的基因組信息（V 2.1），整合了基因組序列和注釋信息，可開展比較基因組學(xué)研究、基因表達(dá)分析、生物信息學(xué)分析和數(shù)據(jù)挖掘等［26］。FungiDB提供利用多種策略在選定的單個(gè)或多個(gè)物種中進(jìn)行基因、ESTs、ORF、基因組序列等的快速檢索，可實(shí)現(xiàn)多步檢索策略的交叉、合并和刪除，創(chuàng)建用戶自己的檢索策略進(jìn)行數(shù)據(jù)分析。同時(shí)提供包括BLAST、序列檢索（Sequence Retrieval）、文獻(xiàn)檢索（PubMed and Entrez）等分析工具。

2.3 卵菌專門數(shù)據(jù)庫(kù)資源

除了綜合性資源外，一些專門數(shù)據(jù)庫(kù)也是卵菌功能基因組學(xué)研究的重要資源。專門數(shù)據(jù)庫(kù)是針對(duì)特定目標(biāo)由一級(jí)數(shù)據(jù)庫(kù)衍生而來(lái)，對(duì)生物學(xué)知識(shí)和信息進(jìn)一步整理形成的二級(jí)數(shù)據(jù)庫(kù)，例如真菌分泌體數(shù)據(jù)庫(kù)FSD，專門存儲(chǔ)那些含有信號(hào)肽、在高爾基體或內(nèi)質(zhì)網(wǎng)進(jìn)行加工能夠分泌到細(xì)胞外的蛋白，這些蛋白通過(guò)在全基因組水平進(jìn)行生物信息學(xué)分析確定，該數(shù)據(jù)庫(kù)中可以查詢到6種已測(cè)序卵菌的所有預(yù)測(cè)獲得的分泌體蛋白［27］。另外，真菌的轉(zhuǎn)錄因子數(shù)據(jù)庫(kù)FTFD［28］、真菌細(xì)胞色素 P450數(shù)據(jù)庫(kù) FCPD［29］等均存儲(chǔ)有卵菌的相應(yīng)蛋白的信息。這些專門數(shù)據(jù)庫(kù)可用于對(duì)于特定基因或領(lǐng)域的研究。

3 總結(jié)

基因組學(xué)極大地提高了我們對(duì)于植物病原菌的理解，例如卵菌病原菌的大量測(cè)序揭示了調(diào)控宿主與寄生物相互作用的效應(yīng)物分子類型［4，30－31］。通過(guò)在全基因組水平上對(duì)數(shù)據(jù)信息的挖掘，尤其是比較基因組學(xué)研究，將對(duì)揭示卵菌的特異性、病原菌群體結(jié)構(gòu)、卵菌與病害關(guān)系、發(fā)現(xiàn)病原菌診斷標(biāo)記起到重要作用［32］。隨著測(cè)序方法的發(fā)展，數(shù)據(jù)產(chǎn)生規(guī)模的提高，對(duì)數(shù)據(jù)的處理和挖掘已經(jīng)成為研究者的新挑戰(zhàn)。

［1］ Erwin D C，Riberro O K.Phytophthora diseases worldwide［M］.St Paul：APS Press，1996：200－562.

［2］ 鄭小波.疫霉菌及其研究技術(shù)［M］.北京：中國(guó)農(nóng)業(yè)出版社，1997：1－29.

［3］ Baldauf S L，Roger A J，Wenk－Siefert I，et al.A kingdom－level phylogeny of eukaryotes based on combined protein data［J］.Science，2000，290（5493）：972－977.

［4］ Tyler B M，Tripathy S，Zhang X，et al.Phytophthora genome sequences uncover evolutionary origins and mechanisms of pathogenesis［J］.Science，2006，313（5791）：1261－1266.

［5］ Haas B J，Kamoun S，Zody M C，et al.Genome sequence and analysis of the Irish potato famine pathogen Phytophthora infestans［J］.Nature，2009，461（7262）：393－398.

［6］ Lamour K H，Mudge J，Gobena D，et al.Genome sequencing and mapping reveal loss of heterozygosity as a mechanism for rapid adaptation in the vegetable pathogen Phytophthora capsici［J］ Molecular Plant－Microbe Interactions，2012，25（10）：1350－1360.

［7］ Lévesque C A，Brouwer H，Cano L，et al.Genome sequence of the necrotrophic plant pathogen Pythium ultimumreveals original pathogenicity mechanisms and effector repertoire［J］.Genome Biology，2010，11（7）：R73－R94.

［8］ Baxter L，Tripathy S，Ishaque N，et al.Signatures of adaptation to obligate biotrophy in the Hyaloperonospora arabidopsidis genome［J］.Science，2010，330（6010）：1549－1551.

［9］ Tian M，Win J，Savory E，et al.454genome sequencing of Pseudoperonospora cubensis reveals effector proteins with a putative QXLR translocation motif［J］.Molecular Plant－Microbe Interactions，2011，24（5）：543－553.

［10］Lamour K H，Kamoun S.Oomycete genetics and genomics：Diversity，interactions and research tools［M］.Hoboken，USA：Wiley－Blackwell Press，2009：540－582.

［11］Richards T A，Dacks J B，Jenkinson J M，et al.Evolution of filamentous plant pathogens：Gene exchange across eukaryotic kingdoms［J］.Current Biology，2006，16（18）：1857－1864.

［12］Kemen E，Gardiner A，Schultz－Larsen T，et al.Gene gain and loss during evolution of obligate parasitism in the white rust pathogen of Arabidopsis thaliana［J］.PLoS Biology，2011，9（7）：e1001094.

［13］Judelson H S.Dynamics and innovations within Oomycete genomes：Insights into biology，pathology，and evolution［J］.Eukaryot Cell，2012，11（11）：1304－1312.

［14］Grigoriev I V，Nordberg H，Shabalov I，et al.The genome portal of the department of Energy Joint Genome Institute［J］.Nucleic Acids Research，2012，40（Database issue）：D26－D32.

［15］Tripathy S，Pandey V N，F(xiàn)ang B，et al.VMD：A community annotation database for oomycetes and microbial genomes［J］.Nucleic Acids Research，2006，34（Database issue）：D379－D381.

［16］Links M G，Holub E，Jiang R H，et al.De novo sequence assembly of Albugo candida reveals a small genome relative to other biotrophic oomycetes［J］.BMC Genomics，2011，12：503－514.

［17］Tripathy S，Deo T，Tyler B M.Oomycete transcriptomics database：A resource for oomycete transcriptomes［J］.BMC Genomics，2012，13（1）：303－316.

［18］Ye W，Wang X，Tao K，et al.Digital gene expression profiling of the Phytophthora sojae transcriptome［J］.Molecular Plant－Microbe Interactions，2011，24（12）：1530－1539.

［19］祁云霞，劉永斌，榮威恒.轉(zhuǎn)錄組研究新技術(shù)：RNA－Seq及其應(yīng)用［J］.遺傳，2011，33（11）：1191－1202.

［20］Hamilton J P，Neeno－Eckwall E C，Adhikari B N，et al.The comprehensive phytopathogen genomics resource：A webbased resource for data－mining plant pathogen genomes［J］.Database：The Journal of Biological Databases and Curation，2011：bar53－65.

［21］Brurberg M B，Elameen A，Le V H，et al.Genetic analysis of Phytophthora infestans populations in the Nordic European countries reveals high genetic variability［J］.Fungal Biology，2011，115（4－5）：335－342.

［22］Vercauteren A，Larsen M，Goss E，et al.Identification of new polymorphic microsatellite markers in the NA1and NA2lineages of Phytophthora ramorum［J］.Mycologia，2011，103（6）：1245－1249.

［23］Schena L，Cardle L，Cooke D E.Use of genome sequence data in the design and testing of SSR markers for Phytophthoraspecies［J］.BMC Genomics，2008，9：620－642.

［24］Garnica D P，Pinzón A M，Quesada－Ocampo L M，et al.Survey and analysis of microsatellites from transcript sequences in Phytophthoraspecies：Frequency，distribution，and potential as markers for the genus［J］.BMC Genomics，2006，7：245－255.

［25］Park J，Park B，Jung K，et al.CFGP：A web－based，comparative fungal genomics platform［J］.Nucleic Acids Research，2008，36（Database issue）：D562－D571.

［26］Stajich J E，Harris T，Brunk B P，et al.FungiDB：An integrated functional genomics database for fungi［J］.Nucleic Acids Research，2012，40（Database issue）：D675－D681.

［27］Choi J，Park J，Kim D，et al.Fungal Secretome Database：Integrated platform for annotation of fungal secretomes［J］.BMC Genomics，2010，11：105－119.

［28］Park J，Park J，Jang S，et al.FTFD：An informatics pipeline supporting phylogenomic analysis of fungal transcription factors［J］.Bioinformatics，2008，24（7）：1024－1025.

［29］Park J，Lee S，Choi J，et al.Fungal cytochrome P450database［J］.BMC Genomics，2008，9：402－412.

［30］Win J，Kamoun S.Adaptive evolution has targeted the C－terminal domain of the RXLR effectors of plant pathogenic oomycetes［J］.Plant Signal Behavior，2008，3（4）：251－253.

［31］Savory E A，Zou C，Adhikari B N，et al.Alternative splicing of a multi－drug transporter from Pseudoperonospora cubensis generates an RXLR effector protein that elicits a rapid cell death［J］.PLoS One，2012，7（4）：e34701－e34713.

［32］Seidl M F，Van den Ackerveken G，Govers F，et al.A domain－centric analysis of oomycete plant pathogen genomes reveals unique protein organization［J］.Plant Physiology，2011，155（2）：628－644.