劉思琪，趙明明

(1.中國(guó)中藥公司，北京 100195；2.中國(guó)醫(yī)學(xué)科學(xué)院 北京協(xié)和醫(yī)院藥用植物研究所，北京 100193)

·綜述·

基于高通量測(cè)序技術(shù)的菌根共生研究概況

劉思琪1，趙明明2*

(1.中國(guó)中藥公司，北京100195；2.中國(guó)醫(yī)學(xué)科學(xué)院 北京協(xié)和醫(yī)院藥用植物研究所，北京100193)

菌根是自然界中普遍存在的一種共生現(xiàn)象，它是由土壤中的菌根真菌與高等植物根系形成的一種共生體。鑒于其在自然界中的重要作用，菌根研究日益引起世界各國(guó)學(xué)者的普遍關(guān)注。本文歸納總結(jié)了主要菌根互作類型叢枝菌根、外生菌根和蘭科菌根利用高通量測(cè)序技術(shù)的研究進(jìn)展。

高通量測(cè)序；叢枝菌根；外生菌根；蘭科菌根

大約有90%的陸地植物可與土壤菌根真菌共生，形成菌根。菌根的形成在土壤結(jié)構(gòu)、植物養(yǎng)分吸收與生長(zhǎng)、生物多樣性及農(nóng)業(yè)和自然生態(tài)系統(tǒng)的生產(chǎn)力等方面發(fā)揮重要作用[1]。根據(jù)形態(tài)結(jié)構(gòu)的不同，菌根可分為叢枝菌根(Arbuscular mycorrhiza，AM)、漿果鵑類菌根(Arbutoid mycorrhiza)、外生菌根(Ectomycorrhiza，EM)、歐石南類菌根(Ericoid)、石晶蘭類菌根(Monotropoid mycorrhiza)和蘭科菌根(Orchid mycorrhiza)。其中最普遍且最有經(jīng)濟(jì)價(jià)值的兩大類當(dāng)屬與農(nóng)業(yè)中絕大部分的谷類作物、蔬菜和水果等共生的叢枝菌根以及與樹(shù)木、灌木共生的外生菌根。

1 叢枝菌根

叢枝菌根(AM)真菌是一類以特定的形式定植于植物體根細(xì)胞，從而協(xié)調(diào)植物生長(zhǎng)的真菌種類。AM真菌從特定的接口入侵植物，從而建立共生關(guān)系[2]?；诜肿由飳W(xué)的叢枝菌根共生研究已開(kāi)展多年，并取得了一定成果。在植物方面，已經(jīng)建立有關(guān)AM共生cDNA的文庫(kù)，通過(guò)抑制性消減雜交技術(shù)獲得大量的克隆，測(cè)序得到AM共生相關(guān)表達(dá)序列標(biāo)簽EST[3]等方面的工作。所獲得的ESTs可以轉(zhuǎn)移到固相載體(陣列)，從而使不同來(lái)源獲得的cDNA雜交RNA以積累AM共生過(guò)程中的基因表達(dá)概況；在真菌方面，EST庫(kù)也正在建立，同樣采用的是直接克隆及抑制性消減雜交技術(shù)。雙方的互作關(guān)系主要表現(xiàn)為一些相關(guān)基因的調(diào)節(jié)和功能驗(yàn)證，以及各自在共生體系中扮演的角色，這已成為科學(xué)界關(guān)注的熱點(diǎn)[4]。

蒺藜苜蓿Medicagotruncatula和百脈根Lotusjaponicas作為AM菌根互作研究的模式植物[5-6]，已分別進(jìn)行了基因組測(cè)序和突變體篩選[7-9]等方面的研究。Frenzel等[10]已成功構(gòu)建的隨機(jī)cDNA文庫(kù)MtAMP，來(lái)源于蒺藜苜蓿與叢枝菌根真菌Glomusintraradices共生后的主要發(fā)育階段的組織樣品，文庫(kù)中共包含了5’端ESTs 3448條。Journet等[11]也在2002年構(gòu)建了含有8567條ESTs的蒺藜苜蓿和叢枝菌根真菌G.intraradices的隨機(jī)cDNA文庫(kù)MtBC。隨后，Liu等[12]在2003年建立了包含9030條ESTs的蒺藜苜蓿和叢枝菌根真菌G.versiforme的隨機(jī)cDNA文庫(kù)MHAMHE MHAM2。由于從隨機(jī)文庫(kù)中發(fā)現(xiàn)高表達(dá)的基因是比較困難的，研究趨勢(shì)轉(zhuǎn)為用差減測(cè)序技術(shù)替代隨機(jī)cDNA文庫(kù)技術(shù)[13]。在采用抑制差減雜交(SSH)cDNA文庫(kù)方法建立的蒺藜苜蓿和叢枝菌根真菌根的差減雜交EST文庫(kù)MtGIM中，對(duì)照組分別是未被侵染的根、磷酸鹽處理的根、苜蓿根瘤菌Sinorhizobiummeliloti侵染的根和豌豆褐斑菌Aphanomyceseuteiches侵染的根，共獲得了1686條ESTs，這些序列中很大一部分被證明是菌根誘導(dǎo)相關(guān)基因[14-15]。AM菌根真菌G.mosseae建立了附著胞生長(zhǎng)相關(guān)的SSH cDNA文庫(kù)，獲得658條ESTs，這些序列被用于研究共生最早期的真菌基因表達(dá)[16]。Ouziad等[17]對(duì)AM菌根真菌G.intraradices進(jìn)行不同外源因子刺激后，菌絲生長(zhǎng)的SSH cDNA文庫(kù)中共獲得2059條ESTs。Hildebrandt等[18]以G.intraradices的萌發(fā)孢子和向外緣生長(zhǎng)的菌絲為材料構(gòu)建的2個(gè)SSH cDNA文庫(kù)中，共獲得3708條ESTs。由菌根組織獲得的ESTs中，真菌序列低于10%，而來(lái)源于AM根的序列主要用于發(fā)現(xiàn)與菌根共生的建立和功能相關(guān)的植物基因[19-20]。

微列陣技術(shù)也是近年來(lái)研究有多少基因能相互作用以及一個(gè)細(xì)胞網(wǎng)如何同時(shí)控制大量基因的新方法。微陣列分為cDNA微陣列和寡聚核苷酸微陣列。微陣列技術(shù)就是利用分子雜交原理，使同時(shí)被比較的標(biāo)本(用同位素或熒光素標(biāo)記)與微陣列雜交，通過(guò)檢測(cè)雜交信號(hào)強(qiáng)度及數(shù)據(jù)處理，把他們轉(zhuǎn)化成不同標(biāo)本中特異基因的豐度，從而全面比較不同標(biāo)本的基因表達(dá)水平的差異[21]。微陣列技術(shù)是一種探索基因組功能的有力手段[22]。在蒺藜苜蓿M.truncatula中首先建立了含有6359個(gè)EST簇的根的微陣列Mt6k-RIT[23]。這些序列代表了來(lái)源于氮缺乏根的cDNA文庫(kù)、幼小根瘤cDNA文庫(kù)和AM菌根cDNA文庫(kù)的21 473條ESTs[24]?；诨ê投骨v隨機(jī)cDNA文庫(kù)中的2516條ESTs、1776個(gè)cDNA克隆經(jīng)聚合酶鏈?zhǔn)椒磻?yīng)(PCR)擴(kuò)增后補(bǔ)充到Mt6k-RIT中，建立了新的Mt8k cDNA芯片，該芯片上大約有6300個(gè)M.truncatula的基因探針[25]。關(guān)于豆科模式植物M.truncatula的其他基因組表達(dá)分析還包括AM菌根的2.5 k cDNA宏陣列和6 k微陣列[26-27]。近幾年來(lái)cDNA探針逐漸被70mer寡聚核苷酸代替[28]。

2 外生菌根

以歐洲山楊Populustremula和毒蠅鵝膏Amanitamuscaria為外生菌根互作研究的模式共生體系中，已建立了含有6669條ESTs的隨機(jī)cDNA文庫(kù)[29]。以外生菌根真菌A.muscaria在低氮源的葡萄糖培養(yǎng)基上生長(zhǎng)的菌絲和無(wú)葡萄糖高氮源培養(yǎng)基上生長(zhǎng)的菌絲為對(duì)照組，構(gòu)建的SSH cDNA文庫(kù)中，共獲得582條ESTs。外生菌根真菌Hebelomacylindrosporum建立了不同氮源條件下菌絲生長(zhǎng)的隨機(jī)cDNA文庫(kù)，共獲得459條ESTs[30]。未來(lái)研究可能更多采用差異表達(dá)分析策略，鑒定菌根共生互作過(guò)程的活性基因，這些策略基于高通量EST測(cè)序、EST聚類和注釋、菌根相關(guān)基因表達(dá)的生物信息學(xué)分析及基于基因芯片技術(shù)的轉(zhuǎn)錄組分析。

3 蘭科菌根

蘭科菌根是一種內(nèi)生菌根，主要寄生于蘭科(Orchidaceae)植物的種子及根系上。目前已開(kāi)展對(duì)蘭科菌根真菌的分類及真菌資源多樣性、蘭科菌根的形態(tài)和菌根對(duì)蘭科植物的效應(yīng)等最新的研究[31]。研究表明，感染蘭科植物根部并能與之共生的真菌絕大多數(shù)屬于擔(dān)子菌門(mén)(Basidiomycota)和半知菌門(mén)(Deuteromycotina)，也有部分屬于子囊菌門(mén)(Ascomycota)；蘭科菌根的形成可分為兩種情況：一是對(duì)蘭科植物種子的侵染；二是對(duì)成長(zhǎng)新根的侵染[32]。菌根真菌對(duì)蘭科植物的種子萌發(fā)和植株生長(zhǎng)發(fā)育均有一定影響。

遺傳信息的豐富不僅有利于發(fā)現(xiàn)功能基因，也對(duì)未來(lái)基因組注釋有幫助，對(duì)觀賞性蘭花黃花杓蘭Cypripediumflavum菌根形成的差異基因表達(dá)研究中，首次獲得蘭科菌根共生相關(guān)差異表達(dá)ESTs[33]。以11種蝴蝶蘭屬蘭花的cDNA文庫(kù)中的37 979 342條序列[34]構(gòu)建了蘭科數(shù)據(jù)庫(kù)OrchidBase，其中41 310條ESTs是由Sanger法獲得的，37 908 032條ESTs由第二代測(cè)序技術(shù)獲得，包括Roche 454測(cè)序技術(shù)和Solexa Illumina測(cè)序技術(shù)所獲得的ESTs經(jīng)過(guò)聚類分析，分為8501個(gè)重疊群(contigs)和76 116單一序列(singletons)，平均片段長(zhǎng)度為459 bp。該數(shù)據(jù)庫(kù)是基于網(wǎng)絡(luò)的蘭科EST數(shù)據(jù)庫(kù)，不僅能查詢到相關(guān)的EST數(shù)據(jù)信息，還提供了聚類信息、功能注釋、基因分析(gene ontology)和代謝途徑分析[35-36]。蘭科數(shù)據(jù)庫(kù)是第一個(gè)在線的蘭科EST序列整合信息庫(kù)，可以自由地上傳或下載相關(guān)數(shù)據(jù)信息[37]。從基因的分子水平和功能分析水平上都為蘭科植物研究提供了豐富的遺傳數(shù)據(jù)信息。

鐵皮石斛DendrobiumofficinaleKimura et Migo為蘭科(Orchidaceae)石斛屬Dendrobium名貴藥材，其種子細(xì)小如粉，自然條件下需共生真菌侵染以提供萌發(fā)所需的碳、氮源等營(yíng)養(yǎng)元素。真菌侵染蘭科植物種子，形成共生關(guān)系是蘭科種子萌發(fā)的前提，但共生機(jī)理尚不清楚。在分子水平上揭示蘭科種子共生萌發(fā)成為近來(lái)科學(xué)界關(guān)注的熱點(diǎn)。趙明明等(2013)以鐵皮石斛種子接種蠟殼菌屬真菌Sebacinasp.共生萌發(fā)至第三階段的組織樣品成功構(gòu)建抑制性差減雜交cDNA文庫(kù)。測(cè)序分析共獲得1437個(gè)EST標(biāo)簽，聚類拼接得到1074個(gè)unigenes，包括172個(gè)重疊群(contigs)和902個(gè)單一序列(singletons)。真核生物蛋白相鄰類聚簇族(KOG)數(shù)據(jù)庫(kù)近似蛋白功能分類、基因本體(GO)注釋及京都基因與基因組百科全書(shū)(KEGG)途徑分析發(fā)現(xiàn)，植物來(lái)源蛋白涉及23個(gè)功能類群，主要包含在信號(hào)轉(zhuǎn)導(dǎo)、抗逆境脅迫及代謝等途徑中[38]。

我國(guó)叢枝菌根、外生菌根和蘭科菌根資源均非常豐富，隨著對(duì)菌根資源的藥用價(jià)值和食用價(jià)值認(rèn)識(shí)的不斷深入，越來(lái)越多的學(xué)者更加關(guān)注菌根互作分子生物學(xué)領(lǐng)域的研究，通過(guò)高通量測(cè)序技術(shù)、微列陣技術(shù)、差減文庫(kù)技術(shù)深層挖掘菌根互作過(guò)程中分子水平的基因表達(dá)信息，這些研究成果為今后菌根資源的合理開(kāi)發(fā)、利用和保護(hù)提供了有力的理論依據(jù)。

[1] Smith S E，Read D J.Mycorrhizal Symbiosis[M].London：Academic Press，1997.

[2] Harrison M J.Signaling in the arbuscular mycorrhizal symbiosis[J].Annual Review of Microbiology，2005，59(1)：19-42.

[3] Gianinazzipearson V，Brechenmacher L.Functional genomics of arbuscular mycorrhiza：decoding the symbiotic c.[J].Canadian Journal of Botany，2004，82(8)：1228-1234(7).

[4] Franken P，Requena N.Analysis of gene expression in arbuscular mycorrhizas：new approaches and challenges.[J].New Phytologist，2001，150(3)：517-523.

[5] Barker D G，Bianchi S，Blondon F，et al.Medicago truncatula，a model plant for studying the molecular genetics of theRhizobium-legume symbiosis[J].Plant Molecular Biology Reporter，1990，8(1)：40-49.

[6] Handberg K，Stougaard J.Lotus japonicus，an autogamous，diploid legume species for classical and molecular genetics[J].Plant Journal，1992，2(4)：487-496.

[7] Tadege M，Ratet P，Mysore K S.Insertional mutagenesis：a Swiss Army knife for functional genomics of Medicago truncatula[J].Trends in Plant Science，2005，10(5)：229-235.

[8] Udvardi M K，Tabata S，Parniske M，et al.Lotus japonicus：legume research in the fast lane.[J].Trends in Plant Science，2005，10(5)：222.

[9] Town C D.Annotating the genome of Medicago truncatula[J].Current Opinion in Plant Biology，2006，9(2)：122.

[10] Frenzel A，Manthey K，Perlick A M，et al.Combined transcriptome profiling reveals a novel family of arbuscular mycorrhizal-specific Medicago truncatula lectin genes[J].Molecular plant-microbe interactions：MPMI，2005，18(8)：771-782.

[11] Journet E P，Van T D，Gouzy J，et al.Exploring root symbiotic programs in the model legume Medicago truncatula using EST analysis.[J].Nucleic Acids Research，2002，30(24)：5579-5592.

[12] Liu J，Blaylock L A，Endre G，et al.Transcript Profiling Coupled with Spatial Expression Analyses Reveals Genes Involved in Distinct Developmental Stages of an Arbuscular Mycorrhizal Symbiosis[J].Plant Cell，2003，15(9)：2106-2123.

[13] Diatchenko L，Lau Y F，Campbell A P，et al.Suppression subtractive hybridization：a method for generating differentially regulated or tissue-specific cDNA probes and libraries[J].Proceedings of the National Academy of Sciences of the United States of America，1996，93(12)：6025-6030.

[14] Wulf A，Manthey K，Doll J，et al.Transcriptional changes in response to arbuscular mycorrhiza development in the model plant Medicago truncatula[J].Molecular plant-microbe interactions：MPMI，2003，16(4)：306-314.

[15] Brechenmacher L，Weidmann S，Tuinen D V，et al.Expression profiling of up-regulated plant and fungal genes in early and late stages of Medicago truncatula-Glomus mosseae interactions[J].Mycorrhiza，2004，14(4)：253-262.

[16] Breuninger M，Requena N.Recognition events in AM symbiosis：analysis of fungal gene expression at the early appressorium stage[J].Fungal Genetics & Biology Fg & B，2004，41(8)：794-804.

[17] Ouziad F，Hildebrandt U，Schmelzer E，et al.Differential gene expressions in arbuscular mycorrhizal-colonized tomato grown under heavy metal stress[J].Journal of Plant Physiology，2005，162(6)：634-649.

[18] Hildebrandt U，Ouziad F，Marner F J，et al.The bacterium Paenibacillus validus stimulates growth of the arbuscular mycorrhizal fungus Glomus intraradices up to the formation of fertile spores[J].Fems Microbiology Letters，2006，254(2)：258-267.

[19] Martin F，Podila G K.Symbiotic Sequencing for the Populus Mesocosm[J].New Phytologist，2004，161(2)：330-335.

[20] Lammers P，Tuskan G A，Difazio S P，et al.Mycorrhizal symbionts of Populus to be sequenced by the United States Department of Energy's Joint Genome Institute.[J].Mycorrhiza，2004，14(1)：63-64.

[21] Lam S H,Mathavan S,Gong Z. Zebrafish Spotted-Microarray for Genome-Wide Expression Profiling Experiments.Part I:Array Printing and Hybridization[J].Methods in Molecular Biology,2009,546(546):175-195.

[22] Quackenbush J.Microarray data normalization and transformation[C]// Nat Genet.2002.

[23] Küster H，Hohnjec N，Krajinski F，et al.Construction and validation of cDNA-based Mt6k-RIT macro- and microarrays to explore root endosymbioses in the model legume Medicago truncatula.[J].Journal of Biotechnology，2004，108(2)：95-113.

[24] Yahyaoui F E，ster，Amor B B，et al.Expression Profiling in Medicago truncatula Identifies More Than 750 Genes Differentially Expressed during Nodulation，Including Many Potential Regulators of the Symbiotic Program[J].Plant Physiology，2004，136(2)：3159-3176.

[25] Firnhaber C，Küster H.EST sequencing and time course microarray hybridizations identify more than 700 Medicago truncatula genes with developmental expression regulation in flowers and pods.[J].Planta，2005，222(2)：269-283.

[26] Lohar D P，Sharopova N，Endre G，et al.Transcript analysis of early nodulation events in Medicago truncatula.[J].Plant Physiology，2006，140(1)：221-234.

[27] Colebatch G，Kloska S，Trevaskis B，et al.Novel aspects of symbiotic nitrogen fixation uncovered by transcript profiling with cDNA arrays[J].Mol Plant Microbe Interact.2002，15(5)：411-420.

[28] Alba R，F(xiàn)ei Z，Payton P，et al.ESTs，cDNA microarrays，and gene expression profiling：tools for dissecting plant physiology and development[J].Plant Journal，2004，39(5)：697-714.

[29] Grunze N，Willmann M，Nehls U.The Impact of Ectomycorrhiza Formation on Monosaccharide Transporter Gene Expression in Poplar Roots[J].New Phytologist，2004，164(1)：147-155.

[30] Wipf D，Benjdia M，Rikirsch E，et al.An expression cDNA library for suppression cloning in yeast mutants，complementation of a yeast his4 mutant and EST analysis from the symbiotic basidiomycete Hebeloma cylindrosporum[J].Genome.2003，46(2)：177-181.

[31] Wipf D，Benjdia M，Rikirsch E，et al.An expression cDNA library for suppression cloning in yeast mutants，complementation of a yeast his4 mutant，and EST analysis from the symbiotic basidiomycete Hebeloma cylindrosporum[J].Genome / National Research Council Canada = Génome / Conseil national de recherches Canada，2003，46(2)：177.

[32] Delseny M，Han B，Yueie H.High throughput DNA sequencing：the new sequencing revolution.[J].Plant Science An International Journal of Experimental Plant Biology，2010，179(5)：407-422.

[33] Watkinson J I，Welbaum G E.Characterization of gene expression in roots of Cypripedium parviflorum var.pubescens incubated with a mycorrhizal fungus[J].Acta Horticulturae，2003,624(64):463-470.

[34] Hsiao Y Y，Tsai W C，Kuoh C S，et al.Comparison of transcripts in Phalaenopsis bellina and Phalaenopsis equestris (Orchidaceae) flowers to deduce monoterpene biosynthesis pathway[J].BMC Plant Biology，2006，6(1)：14.

[35] Otero J T，F(xiàn)lanagan N S.Orchid diversity-beyond deception[J].Trends in Ecology & Evolution，2006，21(2)：65-66.

[36] Peakall R.Speciation in the Orchidaceae：confronting the challenges[J].Mol Ecol.2007，16(14)：2834-2837.

[37] Tsai W C，Hsiao Y Y，Lee S H，et al.Expression analysis of the ESTs derived from the flower buds of Phalaenopsis equestris[J].Plant Sci.2006，170(3)：426-432.

[38] Yu H，Goh C J.Molecular genetics of reproductive biology in orchids[J].Plant Physiology，2001，127(127)：1390-1393.

ProgressofMycorrhizalSymbiosisAnalysisUsingHighThroughputSequencing

LIU Siqi1，ZHAO Mingming2*

(1.ChinaNationalTraditionalChineseMedicineCorporation,Beijing100195,China;2.InstituteofMedicinalPlantDevelopment,ChineseAcademyofMedicalSciencesPekingUnionMedicalCollage,Beijing100193,China)

Mycorrhiza is a common existence in nature, and it is a symbiosis formed by mycorrhizal fungi in soil and roots of higher plants. In view of its important role in nature, mycorrhizal research has attracted more and more attention from scholars all over the world. In this paper, we summarized recent progresses in the research of arbuscular mycorrhiza，ectotrophic mycorrhiza and orchidaceous mycorrhiza symbiosis using high throughput sequencing.

High throughput sequencing；Arbuscular mycorrhiza；Ectotrophic mycorrhiza；Orchidaceous mycorrhiza

] 趙明明，博士后，研究方向：中藥資源學(xué)，E-Mail：mmzhao@icmm.ac.cn

10.13313/j.issn.1673-4890.2017.5.029

2017-05-11)

*[