李凱航 王浩臣 程可心 楊艷 金一 何曉青

（北京林業(yè)大學(xué)生物科學(xué)與技術(shù)學(xué)院，北京 100083）

自然界中，正常生長的植物表面和植物內(nèi)部富集了數(shù)量龐大且種類繁多的微生物，這些微生物的集合被稱為植物微生物組［1］。植物為微生物組（包括細菌、真菌、原生生物、線蟲和病毒）的生長繁殖提供生態(tài)位點，在植物微生物組中數(shù)目最多的是細菌和真菌［2］。植物的健康與高度多樣化和動態(tài)的微生物組有關(guān)，微生物組對其宿主具有重要的功能［3-5］，在促進植物礦質(zhì)元素吸收、提高植物抗病性、協(xié)助植物抵抗非生物脅迫、直接調(diào)節(jié)植物生長發(fā)育等方面有重要影響［6-7］。植物與微生物的互作現(xiàn)象及其機制，是生命科學(xué)研究關(guān)注的熱點，也是利用微生物的關(guān)鍵問題。目前只有一小部分微生物與植物的相互作用得到了深入研究，其他大量微生物種類與功能缺乏系統(tǒng)性研究，它們與植物的互作機理尚不清楚［8-9］。隨著國內(nèi)外微生物組研究的開展，植物微生物組的整體功能更受關(guān)注，闡明植物-微生物組互作的遺傳機制及互作在促進植物健康中的作用，將為植物的優(yōu)良性狀培育提供新的思路，具有重要的科學(xué)意義和應(yīng)用價值。

隨著測序技術(shù)和統(tǒng)計學(xué)的發(fā)展，GWAS 已成為有效方法來識別與微生物組相關(guān)的植物遺傳變異［4］，并且已經(jīng)被應(yīng)用在多種植物上［10-14］。GWAS 是利用全基因組范圍內(nèi)群體中高密度的分子標記，鑒定與復(fù)雜性狀表型變異相關(guān)聯(lián)的分子標記，進而挖掘與表型相關(guān)基因的方法。宿主由于基因型的差異直接影響了表型的變化，因此，GWAS 已經(jīng)被廣泛應(yīng)用于動物和植物以及人體的研究中［15-17］。而由于宿主基因型差異而改變的植物微生物組，同樣影響著植物的健康，可以將微生物組的差異作為植物的“拓展表型”，并利用關(guān)聯(lián)分析找到影響植物微生物組的基因［4］，為分子育種等提供重要的理論基礎(chǔ)和遺傳資源，是未來研究方向的熱點。在植物與微生物組GWAS 分析中，如何表征微生物組差異并用于后續(xù)關(guān)聯(lián)是保證結(jié)果準確性的關(guān)鍵以及難點所在。

1 植物微生物組的結(jié)構(gòu)和功能

植物并非單獨存在，而是與多樣化的微生物共存，在環(huán)境、宿主植物與微生物以及微生物與微生物之間的相互作用下共同塑造植物相關(guān)的微生物群落［7］。在不同的生態(tài)位展現(xiàn)了一個多樣性的植物微生物組，主要包括根際微生物組、葉際微生物組以及內(nèi)生微生物組。受宿主和環(huán)境雙重影響，不同區(qū)系的微生物組在組成和功能上存在差異，是植物與土壤、大氣交互作用的媒介。

土壤作為微生物種子庫，為植物提供了豐富的微生物資源［18］。根際微生物主要來自于土壤，并召集了特定的微生物［19］。Wang 等［20］研究發(fā)現(xiàn)微生物的絕對數(shù)量從土壤到根際是逐漸增加的，提出根際微生物組“擴增-選擇”組裝新模型。有研究表明，相比于土壤，根際中的變形菌門、放線菌門和厚壁菌門細菌略有增加［2,21］，根際微生物中的復(fù)雜微生物群落具有活化根區(qū)養(yǎng)分、促進植物生長、增強植物抗逆、抑制土傳病害等功能［22-27］。除了根瘤中的固氮細菌外，與自然界植物相關(guān)的其他根際微生物也已成為植物氮素營養(yǎng)的巨大來源［28-29］。植物定殖細菌甚至可以調(diào)節(jié)植物開花的時間［30-31］，還能幫助宿主避免滅絕的風(fēng)險［32］。最近的研究表明，植物根際微生物群可緩解植物干旱脅迫，使作物的生產(chǎn)可持續(xù)進行［26,33］。葉際能夠過濾來自更廣闊環(huán)境的微生物，并允許特定類型的微生物在植物定殖，進而服務(wù)于宿主和葉際生態(tài)系統(tǒng)［34-36］。葉際微生物組主要由變形菌門、擬桿菌門、厚壁菌門和放線菌門組成，其中變形菌門占群落組成的近50%［2］。葉際微生物組也已被證明會影響宿主健康和生長、對非生物脅迫的恢復(fù)力和對病原體的抵抗力，具體包括固氮作用、通過產(chǎn)生VOCs 提高抗逆并且阻止病原菌侵染、介導(dǎo)植物防御、抗干旱和紫外線以及促進植物代謝等［37-39］。植物內(nèi)生菌群經(jīng)常富含有變形菌門和厚壁菌門，而擬桿菌門的含量則較低［2］，內(nèi)生菌在促進植物生長、提高抗病性和緩解脅迫耐受性也發(fā)揮了重要作用［2］。根際內(nèi)生菌可以幫助抵御病原菌的入侵，水稻種子內(nèi)生菌可以幫助抗病等［40-41］。

2 植物對微生物組的影響

植物相關(guān)微生物組的巨大功能潛力已經(jīng)得到證實，宿主因素對微生物組的影響也進行了大量的研究［42-43］。在特定的土壤類型及環(huán)境下，植物微生物組差異主要由宿主調(diào)控［44-45］。植物在募集微生物組、篩選有益微生物并且抵御有害菌以及調(diào)控微生物之間的互作等方面發(fā)揮重要作用。不同植物品種及基因型在對微生物的募集上表現(xiàn)出明顯的差異［21,26,46-47］。植物基因型影響了根系代謝物、免疫系統(tǒng)功能及根系分泌物的組成，進而影響了微生物組的活動與結(jié)構(gòu)［6,48］。在面對生物或者非生物脅迫時，植物會通過遺傳因子整合脅迫信號并參與主動重塑植物微生物群，利用“呼救”（cry for help）策略招募有益菌［46,49-50］，例如，小麥植株招募嗜麥芽寡養(yǎng)單孢菌（Stenotrophomonas rhizophilaSR80）來抑制土傳病原體［51］。甜菜根在根內(nèi)層富集噬幾丁質(zhì)菌科（Chitinophagaceae）和黃桿菌科（Flavobacteriaceae），這促進了植物抗真菌［40］。在最近的一項研究中，通過比較已感染受柑橘黑點病菌Diaporthe citri和未感染受柑橘黑點病菌D. citri柑橘葉之間的葉際微生物組的變化，發(fā)現(xiàn)柑橘會主動招募有益菌（包括泛菌Pantoeaasv90 和甲基桿菌Methylobacteriumasv41等）［39］。對優(yōu)先考慮的番茄數(shù)量性狀位點的基因內(nèi)容分析表明，各種根際菌系的不同招募有其遺傳基礎(chǔ)［52］。通過對143 個玉米自交系材料進行根系細菌群落組成研究發(fā)現(xiàn)，玉米根系細菌存在差異［53-54］。也有研究發(fā)現(xiàn)，由于物種繼承了祖先的遺傳物質(zhì)和性狀，親緣關(guān)系越密切的物種比進化距離遙遠的物種微生物組更相似［55］，說明宿主遺傳因素影響了微生物組組成。Bergelson 等［42］總結(jié)了植物影響微生物組的主要通路，包括生物鐘、生長發(fā)育、物理屏障、免疫系統(tǒng)、植物激素、初級代謝、次級代謝以及植物與微生物的共生。

2.1 植物基因型差異對微生物組網(wǎng)絡(luò)的影響

微生物在植物的作用下，通過競爭或者合作等方式，形成網(wǎng)絡(luò)結(jié)構(gòu)來維持微生物組內(nèi)部的穩(wěn)定［6］。Xiong 等［56］研究發(fā)現(xiàn)，由于植物的選擇效應(yīng)，影響了植物微生物組群落構(gòu)建及網(wǎng)絡(luò)復(fù)雜度，網(wǎng)絡(luò)復(fù)雜度從土壤到根際到植物內(nèi)部逐漸降低，并以葉內(nèi)最低。不同生態(tài)位的微生物組也因植物基因型因素的差異而有區(qū)別，包括根際和葉際等（圖1）［29,34］。盡管研究人員已在多種植物中開展植物基因型影響微生物組結(jié)構(gòu)的研究，但對植物如何控制微生物組平衡來影響其健康的探究相對不足［55］，參與招募并塑造微生物組的主要基因有待挖掘。

圖1 植物基因型差異影響根際和葉際的微生物互作網(wǎng)絡(luò)Fig. 1 Plant genotype differences affecting the microbial interaction network at the rhizosphere and phyllosphere

2.2 遺傳力表征植物對微生物以及微生物組的影響

為了解釋基因型和數(shù)量性狀之間復(fù)雜的關(guān)系，育種學(xué)家和遺傳學(xué)家已通過計算遺傳力來表示基因型對表型的影響大小。某一表型數(shù)據(jù)遺傳力越大，表示其與宿主基因型關(guān)系越大。植物微生物組差異作為一種數(shù)量性狀，也將遺傳力計算應(yīng)用在植物基因?qū)ξ⑸锝M的影響上。諸多學(xué)者將植物微生物組豐度、α 和β 多樣性等參數(shù)視為數(shù)量性狀，估計它們的遺傳力，辨別受宿主遺傳因素調(diào)控的微生物群落，并進一步通過GWAS 分析挖掘?qū)е轮参镂⑸锝M可遺傳的宿主遺傳變異［11-12］。研究微生物組遺傳力，一方面可以加強我們對植物微生物組互作的影響，另一方面也可以作為幫助我們篩選益生菌的途徑，加強對有益微生物的利用。

針對于植物微生物組的遺傳力，其大小并非表示后代傳遞微生物組的能力，更應(yīng)該準確的表示為描述最終由宿主遺傳引起的群落組裝的決定性成分［57］。通過對玉米根際微生物研究發(fā)現(xiàn)，單個微生物遺傳力大小范圍為0.15-0.25，雖然相比于開花時間的和其他傳統(tǒng)的農(nóng)藝性狀來說遺傳力較低，但仍能體現(xiàn)出植物基因?qū)μ囟ǚ诸惾旱挠绊?。其中，遺傳力較高的微生物普遍為α 變形菌綱，并且發(fā)現(xiàn)7 種核心微生物（多個樣品中普遍出現(xiàn)的微生物）中有5 種屬于高遺傳力的微生物［54］。Brachi等［58］以擬南芥葉際微生物為研究對象，篩選出了既是可遺傳的并且在微生物互作網(wǎng)絡(luò)中為hub 微生物，屬于原霉菌（Protomyces inouyei）。Li 等［59］計算了單個微生物的遺傳力，發(fā)現(xiàn)有15 個屬于網(wǎng)絡(luò)中hub 微生物同樣具有相對較高的遺傳力（H2>0.10），包括變形菌門（Proteobacteria）、擔(dān)子菌門（Basidiomycota）、子囊菌門（Ascomycota）以及放線菌門（Actinobacteria）。Deng 等［12］計算了高粱根際微生物組的遺傳力值，以降維后的主成分值數(shù)據(jù)作為不同樣品的微生物組結(jié)構(gòu)，研究發(fā)現(xiàn)，第一主成分解釋了微生物群落21%的方差，并且第一主成分遺傳力為0.35，證明了在高粱中微生物組具有較高的遺傳力。Wallace 等［11］對玉米葉際微生物組研究，計算出某些單個的微生物遺傳力較高，包括甲基桿菌屬（Methylobacterium，H2=0.457）等，利用PICRUST 預(yù)測的方法，提出大量微生物功能也具有遺傳性。

3 全基因組關(guān)聯(lián)分析在植物微生物組中的應(yīng)用

研究植物與微生物組的互作遺傳機制需要構(gòu)建系統(tǒng)的方法［4,19］。GWAS 是應(yīng)對這一挑戰(zhàn)的一個強大而靈活的方法，由于GWAS 分析只需要表型數(shù)據(jù)和基因型數(shù)據(jù)，通過鑒定對宿主基因型敏感的微生物或者微生物組，并將它們與影響其定殖的基因位點聯(lián)系起來，幫助我們挖掘各種性狀顯著相關(guān)的數(shù)量性狀座位（SNP）及其候選基因位點。因此，全基因組關(guān)聯(lián)分析作為一種有效的手段幫助我們理解植物關(guān)鍵調(diào)控因子和相關(guān)信號途徑如何主動重塑相關(guān)微生物組組裝以及調(diào)控微生物組的穩(wěn)定。

3.1 植物與微生物組關(guān)聯(lián)分析的數(shù)量遺傳模型的應(yīng)用

GWAS 的基礎(chǔ)是連鎖不平衡（linkage disequilibrium, LD）［60］，LD 的大小主要受群體遺傳多樣性的影響, 在不同物種和群體中差異很大。GWAS 需要考慮的問題包括群體的選取、群體結(jié)構(gòu)分析、表型鑒定、全基因組關(guān)聯(lián)分析方法選擇及結(jié)果矯正［61］。在選擇微生物組作為植物的表型進行關(guān)聯(lián)分析時，首先是考慮對樣品的選擇，樣本量是保證結(jié)果可靠性的一個重要因素，一般選擇200 個以上［62］。然后進行群體結(jié)構(gòu)分析，群體結(jié)構(gòu)分析主要分為兩類，包括主成分分析（principal component analysis, PCA）和顯示生成模型的方法［61］。在獲得基因型數(shù)據(jù)和表型數(shù)據(jù)后進行關(guān)聯(lián)分析時，由于宿主群體結(jié)構(gòu)以及個體間的關(guān)系，會導(dǎo)致關(guān)聯(lián)結(jié)果出現(xiàn)假陽性，在做關(guān)聯(lián)分析時，我們需要考慮群體結(jié)構(gòu)（G）和親緣關(guān)系（K）。并將其作為協(xié)變量，通過引入一般線性模型（general linear model, GLM）和混合線性模型（mixed linear model, MLM）來矯正群落結(jié)構(gòu)和親緣關(guān)系的影響，計算膨脹系數(shù)并結(jié)合Quantile-Quantile 圖選擇合適的模型篩選顯著性位點［61］。目前可用于GWAS的軟件包越來越多，Yin 等［63］利用R 語言，創(chuàng)建了rMVP 包，通過引入GLM、MLM 以及FarmCPU 三種模型，采用分塊矩陣計算、并行加速、計算流程優(yōu)化等策略，加快了計算。同時，更多的GWAS 模型包括Meta-GWAS、多表型GWAS 等，其用于植物微生物組關(guān)聯(lián)分析上仍有待探索［4］。

3.2 微生物組差異作為植物表型的選擇

在應(yīng)用關(guān)聯(lián)分析揭示植物對微生物組的調(diào)控中，如何表征以及量化植物體內(nèi)的微生物組是我們需要思考的問題。傳統(tǒng)的研究主要依賴于在實驗室條件下開展植物與單一微生物的互作關(guān)系研究，很少在自然狀態(tài)下研究微生物組與宿主植物共存的機制［1］，例如病原菌的豐度會影響植物健康，早期人們?yōu)榱搜芯坎≡谒拗鞯年P(guān)系，以病原菌豐度為表型數(shù)據(jù)調(diào)查植物對病原菌的調(diào)控［64］；也有研究者對共生微生物進行研究，比如豆科植物和根瘤菌之間的關(guān)系［65］。隨著測序技術(shù)的發(fā)展，擴增子測序使得人們可以了解到植物體內(nèi)的微生物群落。因此也有不少研究以相對豐度最高的微生物為表型數(shù)據(jù)進行關(guān)聯(lián)分析［10］。Bergelson 等［42］提出只考慮相對豐度高的微生物會忽略遺傳力高的微生物，確定可遺傳的微生物可能有助于建立宿主遺傳變異和微生物生態(tài)網(wǎng)絡(luò)之間的聯(lián)系。

當前人們越來越重視微生物組，希望可以深入了解植物與微生物組的互作機制。已有多種表征微生物組差異的類型，包括α 多樣性指數(shù)、β 多樣性指數(shù)、微生物組功能差異、hub 微生物的相對豐度及網(wǎng)絡(luò)特征值［4,13-14］。α 多樣性主要關(guān)注局域均勻生境下的物種數(shù)目，因此也被稱為生境內(nèi)的多樣性，包括了Chao1、Ace、Shannon、Simpson、Richness 等指數(shù)。β 多樣性不僅可以反映樣本之間的多樣性距離關(guān)系，而且還可以反映生物群落之間的分化程度。其使用統(tǒng)計算法Euclidean，Bray-Curtis，Unweighted_unifrac，weighted_unifrac 等，計算兩兩樣品間距離，獲得距離矩陣，通過降維方式展示微生物組差異，包括主成分分析（principal component analysis，PCA）、主坐標分析（principal co-ordinates analysis，PCoA）、非度量多維尺度分析（non-metric multidimensional scaling，NMDS 分 析）、非 加 權(quán) 組平均聚類分析（unweighted pair-group method with arithmetic means，UPGMA）等分析方法，從中發(fā)現(xiàn)不同樣品之間微生物組的差異。還可以利用一些功能數(shù)據(jù)來表示微生物組之間的差異，包括利用宏基因組展示微生物組功能差異、基于擴增子測序數(shù)據(jù)預(yù)測的微生物群落功能（利用PICRUSt、Tax4Fun、FAPROTAX 及BugBase 等方法）［66-67］。

植物微生物組之間通過互作方式形成復(fù)雜的網(wǎng)絡(luò)結(jié)構(gòu)，特別是細菌-真菌之間的相互作用，在植物生長中發(fā)揮重要作用［68］。在陸地生態(tài)系統(tǒng)的養(yǎng)分循環(huán)及植物群落動態(tài)等過程有重要調(diào)節(jié)作用［69-70］。微生物組差異也可以通過一些網(wǎng)絡(luò)參數(shù)體現(xiàn)出來，其中，一些hub 微生物在網(wǎng)絡(luò)中發(fā)揮重要作用，宿主也會通過影響hub 微生物來調(diào)控微生物互作網(wǎng)絡(luò)［71］。因此，有研究利用不同樣品間微生物互作網(wǎng)絡(luò)中的hub 微生物為表型數(shù)據(jù)來代表不同樣品之間的微生物組差異。通過了解對核心微生物的分析加強對微生物組的理解［13］。Jiang 等［72］用數(shù)學(xué)模型將微生物的競爭與合作系統(tǒng)分解為共生、拮抗、攻擊和利他網(wǎng)絡(luò)，He 等［14］整合行為生態(tài)學(xué)與基因作圖理論，將網(wǎng)絡(luò)嵌入到基因定位的框架，構(gòu)建了網(wǎng)絡(luò)作圖模型（network mapping）（圖2）。利用數(shù)學(xué)模型計算了不同互作關(guān)系下的六種微生物互作網(wǎng)絡(luò)特征值（Connectivity、Closeness、Betweenness、Eccentricity、Eigenvector、PageRank），并以網(wǎng)絡(luò)特征值來表征不同樣品之間微生物互作網(wǎng)絡(luò)的差異。

圖2 網(wǎng)絡(luò)作圖研究植物與微生物組互作Fig. 2 Studying plant-microbiome interactions by network mapping

3.3 利用GWAS進行的植物微生物組研究

選擇合適的微生物組表型差異數(shù)據(jù)后，可以將基因型數(shù)據(jù)與表型數(shù)據(jù)進行關(guān)聯(lián)，挖掘調(diào)控微生物組的植物候選基因。Wallace 等［11］計算了玉米葉際微生物組的α 和β 多樣性的多種參數(shù)（α 多樣性指數(shù)計算了除Michaelis-Menten 之外的其他所有參數(shù)，利用加權(quán)的UniFrac、未加權(quán)的UniFrac 和Bray-Curtis 距離度量計算了β 多樣性，同時利用PICRUST 的方法預(yù)測了微生物組的功能）并以此為表型數(shù)據(jù)進行關(guān)聯(lián)分析，發(fā)現(xiàn)了影響玉米葉際微生物組的管家基因（包括AGPv3等）。在最近的一項研究中，通過對高粱根際微生物組進行分析，以不同高粱根際微生物組PCA 結(jié)果的第一主成分來表征不同樣品間微生物組差異，選到RGA2、CHAF1A等基因，結(jié)果顯示這些基因表現(xiàn)出較強的根特異性活性［12］。Koprivova 等［73］以擬南芥根際土壤中微生物組硫酸鹽酶活性為表型數(shù)據(jù)，提出編碼色素P450 的基因影響根際硫酸鹽酶活性來反映出微生物組的差異。有研究以水稻葉際微生物互作網(wǎng)絡(luò)中的hub 微生物的相對豐度為表型數(shù)據(jù)，獲得了編碼氧化氫同工酶、ATP 依賴的蛋白酶等的基因（OsNippo02g019000、OsNippo11g101050）［13］；He 等［14］將 擬 南 芥 基 因型與根際微生物組網(wǎng)絡(luò)特征值進行關(guān)聯(lián)，篩選到一些影響擬南芥根系生長以及一些抗逆基因，例如TMK3，IBR1B等；Li 等［59］同樣利用網(wǎng)絡(luò)作圖發(fā)現(xiàn)一些在植物生長中發(fā)揮重要作用的基因也影響了擬南芥葉際微生物互作，包括ADA2B,HAL3A等。同時，我們總結(jié)了在不同植物以及不同生態(tài)位中與微生物組有關(guān)的植物基因（表1）

表1 與塑造微生物組有關(guān)的植物基因Table 1 Plant genes associated with the microbiome

4 展望

植物與自身微生物組作為一個整體，對于植物調(diào)控微生物組組裝的機制已有相關(guān)研究，微生物組在植物抗逆抗病等方面的作用也得到了廣泛證實，植物微生物組對植物生長發(fā)育的影響及其與植物的互作機制是國際前沿研究領(lǐng)域和重點研究方向［21］。微生物組差異作為植物的“拓展表型”應(yīng)用于GWAS 幫助我們尋找可以調(diào)控微生物組的顯著位點，但是仍然有很多未來需要關(guān)注的方向。

首先，隨著測序技術(shù)迅速發(fā)展，提高了我們對微生物組的認識。為了進一步研究植物對微生物組的調(diào)控，我們?nèi)匀恍枰獜亩嘟M學(xué)的角度出發(fā)，來量化微生物組功能差異，包括代謝組學(xué)以及蛋白質(zhì)組學(xué)等，并將其作為表型數(shù)據(jù)引入到GWAS 分析中成為未來發(fā)展的方向。其次，為了充分發(fā)揮植物與微生物組之間的有益協(xié)作關(guān)系，保證微生物組在可持續(xù)農(nóng)業(yè)中發(fā)揮更積極的作用，將篩選到的顯著基因應(yīng)用到植物中，通過改善植物基因表達情況來操控微生物群的穩(wěn)態(tài)或者優(yōu)化微生物群，進而幫助植物健康生長和抵抗脅迫環(huán)境。

對于在關(guān)聯(lián)分析中尋找到的益生菌，合成菌群（SynComs）作為一種有效的手段，已經(jīng)嘗試將微生物組利用在植物上［82］。有關(guān)微生物群落中樞紐菌株和關(guān)鍵物種的新信息為利用有益微生物提供了進一步的指導(dǎo)，通過培養(yǎng)和篩選核心微生物，進一步創(chuàng)制新型菌劑，作為挖掘利用微生物組的重要途經(jīng)，加強對有益微生物的利用。