土壤宏病毒組的研究方法與進展*

祁慧鹓，鄭曉璇，孫明明，王金鋒，馬迎飛，朱 冬，王風賀，蔣 新，葉 茂?

（1. 中國科學院南京土壤研究所，土壤環境與污染修復重點實驗室，南京 210008；2. 中國科學院大學，北京 100049；3. 南京農業大學資源與環境科學學院，土壤生態學實驗室，南京 210095；4. 中國科學院北京生命科學研究院，北京 100101；5. 中國科學院深圳先進技術研究院，深圳 518055；6. 中國科學院生態環境研究中心，北京 100085；7. 南京師范大學環境學院，南京 210023）

病毒是地球上數量最多的生命體，廣泛存在于地球各種環境中[1]。土壤病毒在調控微生物群落組成、影響元素循環利用、促進生物進化等方面發揮重要作用[2-4]。近年來，海洋環境中病毒生態學研究進展迅速，雖然學術界逐漸意識到病毒在土壤環境中發揮著重要作用，但由于病毒缺少通用的標記基因，以及受到土壤異質性、多樣性的限制，陸地系統，特別是土壤環境中病毒研究進展相對緩慢[5-6]，土壤病毒學研究遠落后于海洋等水體環境[2]。噬菌體是寄生在細菌、古菌等原核生物里的病毒，也是土壤中最主要的病毒類群。單株噬菌體的全基因組測序有助于了解噬菌體的基因組結構及功能特征；而隨著測序技術的發展，宏病毒組學研究日益受到關注，它擺脫了以往微生物分離純培養的限制，以環境樣品為研究對象，直接從土壤樣品中富集和提取病毒基因組后進行測序和生物信息學分析[7-8]，為土壤病毒研究提供了新的技術手段[9]。

由此，本文主要介紹了土壤研究方法及病毒組研究進展，指出了土壤病毒研究當前面臨的困境，并對未來的研究方向和發展趨勢進行了展望，以期為后續土壤宏病毒研究提供科學參考。

1 土壤宏病毒組的研究方法

1.1 土壤病毒核酸提取

大部分土壤病毒被土壤顆粒吸附固定，對土壤病毒研究的關鍵是建立精準、專性、高效的提取方法，獲得吸附于土壤顆粒表面及內部的病毒粒子[10]。但土壤環境復雜且微生物多樣性高，在提取病毒宏基因組核酸時易混入細菌、真菌等其他微生物的核酸序列，從而影響后續基因序列分析、功能基因注釋的準確性。因此，土壤病毒提取方法的選擇（如濾膜、提取劑等）對后續病毒核酸純度、功能基因注釋的準確性具有較大影響。G?ller等[6]分別使用0.22 μm和0.45 μm濾膜過濾病毒，發現0.22 μm濾膜可去除更多的細菌DNA，且不影響土壤病毒多樣性，而通過0.45 μm孔徑過濾的病毒體中不僅細菌污染程度增加，病毒多樣性也有所下降。Williamson等[11]比較了10%（w/v）的牛肉膏、250 mmol·L-1甘氨酸溶液、10 mmol·L-1焦磷酸鈉和1%（w/v）檸檬酸鉀溶液四種提取劑對粉砂壤土和砂壤土中病毒的提取效果。通過活菌計數法（也稱間接計數法，僅測活菌總數）計算提取率，發現250 mmol·L-1甘氨酸溶液提取率可達28.9%；其次是10%（w/v）牛肉膏，提取率為26.0%；而1%（w/v）檸檬酸鉀溶液和10 mmol·L-1焦磷酸鈉提取率較低，分別為16.9%和15.0%。后續熒光顯微鏡計數時，發現牛肉膏提取物、粉砂壤土中焦磷酸鈉提取物和砂壤土中的甘氨酸提取物均無法計數類病毒顆粒（virus-like particles，VLPs），故該研究認為1%（w/v）檸檬酸鉀溶液是最佳的土壤病毒提取劑。與此同時，土壤pH、土壤陽離子交換量、溫度、含水量等環境因子均會影響土壤病毒粒子的提取效率[10]。此外，當病毒遺傳物質含量較低時，會影響后續宏病毒組分析，雖可經多重置換擴增（multiple displacement amplification，MDA）或引物延伸預擴增（primer extension preamplification，PEP）等全基因組擴增技術增加病毒核酸濃度，但這些技術在擴增產物上可能存在偏好性。為保證結果的準確性，建議提高原始土壤病毒樣品富集量，再進行宏病毒組基因測序[2]。

1.2 測序與病毒序列鑒定

第二代測序技術具有通量高、速度快、成本低的優點[12-13]，使得宏病毒組研究發生巨大變化。其中Illumina目前擁有MiSeq、HiSeq、NovaSeq等多個技術平臺，并有多種讀長和通量模式供選擇。但總體而言，第二代測序技術仍存在序列讀長短的局限性。近年來，第三代測序技術的快速發展在宏病毒組研究中顯示出巨大潛力[14]，它無需進行PCR擴增，且具有讀長長的優點。隨著第三代測序技術的不斷成熟及測序成本的降低，預計它的普及程度將不斷提升。

在含有病毒和宿主的混合基因組數據中，鑒定病毒序列是解析病毒信息的關鍵[15]，會直接影響后續病毒功能基因注釋的效果。目前，識別完整微生物基因組中原噬菌體的工具主要有Phage_Finder[16]，Prophinder[17]和PHAST。PHAST速度快、準確性高的優勢使其成為極具吸引力的工具，后續由于序列數據庫規模的擴大及用戶數量的增加，PHAST發行了新版本PHASTER[18]。與PHAST相比，新版本的突出優勢在于可以識別宏基因組拼接產物的原噬菌體[19]。盡管PHASTER與PHAST是在細菌基因組中鑒別原噬菌體的兩個廣泛使用的工具，但值得注意的是，對預測的原噬菌體精準程度，如融合位點的位置仍存在一些不確定性[19]。此外，這些軟件大多沒有設計專門的算法用于相對較短的重疊群或支架，且無法在短時間內同時處理大量的序列，因此并不適合從宏基因組數據中鑒別病毒序列[15]。2015年Roux等[20]開發了一款工具—VirSorter，不僅可識別完整微生物基因組中的原噬菌體，還可用于檢測拼接宏基因組數據中的病毒序列。VirSorter雖然在很大程度上依賴于對已有病毒基因組的相似性搜索，但它卻使用了一個自定義的病毒參考基因組數據庫，增加了從淡水、海水、人體腸道、肺和唾液中取樣的宏病毒組序列[15]。2017年Ren等[15]開發了一款基于k-mer的工具VirFinder，用于從宏基因組數據中識別原核病毒序列。與基于基因的病毒分類工具VirSorter相比，VirFinder在識別病毒序列方面明顯優于VirSorter。在模擬數據集中，VirFinder從1、3和5 kb序列中識別病毒序列的真陽性率分別較VirSorter高78倍、2.4倍和1.8倍；其假陽性率與VirSorter相同，表明VirFinder真陽性更高、對短片段序列的識別效果更好[15]。但這兩款軟件也有其局限性，VirSorter和VirFinder為檢測細菌和古細菌病毒而優化設計，不能很好地檢測真核病毒；且兩者在微生物群落中對病毒分析的功能相對有限，在鑒定出病毒序列后，沒有進一步分析病毒與宿主的對應關系。相比之下，VirMiner是一款病毒預測與分析宏病毒組樣本的綜合工具[21]，能夠捕獲到高豐度的噬菌體序列，這些噬菌體在感染細菌和影響微生物群落動態方面起著關鍵作用；更重要的是，VirMiner提供了更全面的噬菌體分析流程，包括宏基因組原始讀段處理、功能注釋、噬菌體序列鑒定、噬菌體-宿主侵染關系預測；此外，當宏基因組序列包括不同條件的數據（如處理組和對照組）時，還可支持不同組之間的統計比較。

1.3 病毒功能基因注釋

病毒功能基因注釋是將預測出的編碼基因通過與相關數據庫的參考序列進行比對，在與現有病毒進行同源性搜索的基礎上，獲取該基因的功能信息。通過對病毒功能基因注釋，不僅為深入認識病毒個體生命過程提供理論基礎；還有助于了解病毒群落的生態過程及與宿主群落的生態互作關系，從而闡釋病毒與宿主、環境間復雜的相互作用機制。本文重點從單個噬菌體的全基因組測序及土壤宏病毒組兩方面闡釋了土壤病毒的功能基因注釋。

單個土壤噬菌體的全基因組注釋流程，首先從土壤篩選分離得到噬菌體，經純化、濃縮后采用透射電子顯微鏡觀察其形態特征；隨后提取噬菌體核酸，進行全基因組測序，在過濾掉低質量序列后進行全基因組序列組裝，并通過注釋軟件等在線工具，對噬菌體全基因組序列進行功能注釋[22]。

土壤宏病毒組的功能基因注釋由土壤樣品制備與宏病毒組分析兩個主要部分組成（圖1），具體步驟包括：（1）根據研究目的從相應的土壤中采集樣本，并根據病毒類型、土壤理化性質等選擇合適的緩沖液進行土壤病毒提取，隨后對提取液進行過濾，去除細菌等其他潛在宿主，進而濃縮富集、純化病毒；（2）通過病毒核酸提取試劑盒或手工提取的方式獲取病毒核酸、構建測序文庫、并通過測序平臺進行宏病毒組測序；（3）對測序得到的原始數據進行質量控制，基于重疊區（overlap）將高質量測序讀段（reads）拼接為重疊群（contigs），進一步組裝成支架（scaffolds）；（4）通過病毒序列識別軟件在重疊群或支架中判別、篩選出病毒序列；（5）使用基因預測軟件對病毒基因組ORF進行預測，再通過注釋工具將ORF與多個數據庫比對進行功能基因注釋；（6）自動注釋后，為保證結果準確性，可手動修正自動注釋結果并進行適當補充。

目前，普遍采用與各種數據庫進行蛋白質序列比對的方法，對樣本中的基因功能進行注釋分析。在注釋過程中，研究人員通常根據自身需求選擇合適的功能數據庫。其中Pfam數據庫是一個基于多序列比對和隱馬爾可夫模型的蛋白質結構域和家族數據庫，可提供蛋白質家族和結構域的完整準確分類，被廣泛用于查詢蛋白家族或蛋白結構域的注釋，結構及多序列比對信息，在基因功能注釋上可用性較強[23-24]。它有A和B 兩個質量級別的家族數據庫，Pfam A通過比對人工校正過的種子序列，并使用隱馬爾可夫模型進行選擇，數據質量較高；Pfam B為算法自動生成，雖可靠性降低，但也可以被用于鑒別功能保守的區域。GO（gene ontology）數據庫分別從細胞學組件、分子功能、生物學途徑對基因表達產物進行簡單注釋。經GO數據庫注釋后，可得到基因在不同類別中注釋的具體情況。而KEGG（kyoto encyclopedia of genes and genomes）是一個系統分析基因功能的知識庫，核心為KEGG PATHWAY數據庫。利用KEGG進行注釋后，能清楚地反映出基因與相關代謝的關系[25]。COG（clusters of orthologous groups）數據庫是由NCBI創建并維護的蛋白數據庫[26]。隨著測序技術的發展，COG陸續在不同物種中建立相關的同源蛋白簇。基于完整噬菌體基因組中的編碼蛋白系統進化關系，構建而成的POG（phage orthologous groups）數據庫便是其中的一個分支[27]。通過比對，可將某個特定蛋白序列注釋到一個由直系同源序列構成的POG中，從而推測該序列的功能。此外，POG數據庫包含了進化過程中基因得失信息，還可用于系統發育的統計推斷和祖先基因組的重建[27]。CAZy（carbohydrate-active enzyme）數據庫則針對性較強，是一類與合成或分解復雜碳水化合物和糖復合物酶類有關的數據庫，可提供碳水化合物酶類物種來源、酶功能分類、基因序列、蛋白質序列及其結構等信息[28]。還有一些小眾數據庫，如抗性基因數據庫CARD（comprehensive antibiotic resistance database）在細菌耐藥性的分子基礎上，提供了參考DNA和蛋白質序列、檢測模型和生物信息學工具[29]。ARO（antibiotic resistance ontology）是該數據庫的核心，包含了與抗生素抗性基因、抗性機制、抗生素相關的條目。通過與該數據庫進行比對，可得到與耐藥基因相關的注釋信息。

2 土壤噬菌體基因組研究進展

噬菌體是土壤中最主要的病毒類群，對單個噬菌體進行全基因組測序和功能基因注釋，有助于探明該噬菌體的基因信息和功能特征，進而挖掘該噬菌體在環境修復、疾病治療等實際應用中的巨大潛力。如近年來，噬菌體療法作為一種可以高效靶向追蹤滅活土壤體系中致病細菌的有效手段而日益受到關注[30]，故分離篩選出新型烈性噬菌體菌株資源對于噬菌體療法具有重要意義。但由于一些噬菌體攜帶毒力基因等原因，導致噬菌體防治存在“雙刃劍”風險[31]。因此，從土壤中分離單株噬菌體，并對其進行全基因組測序、功能基因注釋有助于我們更好的了解噬菌體的基因組結構及功能特征。

2.1 全基因組測序下的土壤噬菌體研究趨勢

借助Web of Science核心合集以“complete genome sequence”和“phage”為關鍵詞進行檢索，在2000年-2020年時間段內共有1 508篇相關文獻。發現近20年該領域發文量呈現整體上升趨勢（圖2a）。隨后利用VOSviewer[32]可視化分析軟件對關鍵詞進行聚類分析（圖2b）。圖中節點大小表示關鍵詞出現的頻率，顏色反映不同時間的研究熱點，其中黃色部分代表了較為前沿的關鍵詞。總體而言，“開放閱讀框 open reading frame”、“原噬菌體prophage”、“病原體pathogen”、“表征 characterization”、“作用 role”是出現頻率較高的詞，表明對噬菌體的研究主要聚焦于其形態和功能。值得注意的是，“土壤soil”、“抗生素抗性antibiotic resistance”、“噬菌體療法phage therapy”是近幾年出現的關鍵詞。“土壤 soil”的出現表明土壤噬菌體領域的研究將日益受到重視；此外，由于近年來土壤中抗生素抗性基因的增多、超級細菌的出現，噬菌體療法可能成為未來土壤中病原細菌滅活的重要手段。因此，研發對病原細菌具有廣譜性的混合噬菌體雞尾酒制劑，探明噬菌體療法對于土壤微生物群落結構、功能及養分循環的影響可能是未來土壤噬菌體領域的研究熱點。

2.2 全基因組測序下的土壤噬菌體研究進展

蘇靖芳等人[31]以引起煙草以及多種茄科植物產生萎蔫癥的青枯雷爾氏菌（Ralstonia solanacearum）為宿主，采用雙層平板法從煙田土壤中分離出一株烈性噬菌體RS-PII-1。隨后對噬菌體RS-PII-1進行全基因組測序，并通過RAST在線軟件對噬菌體全基因組序列進行功能注釋，發現該基因組含有46個ORF，其中25個ORF可被注釋為相應的功能蛋白、結構蛋白，但仍存在21個功能未知的ORF假想蛋白（hypothetical protein）。通過與已知噬菌體全基因組序列比較分析以及噬菌體RS-PII-1系統進化樹的構建，發現噬菌體RSJ2、RSB1與RS-PII-1相似度最高，但噬菌體RS-PII-1基因組中功能蛋白和假想蛋白區域與兩株噬菌體仍存在明顯差異，表明RS-PII-1是一株全新的青枯雷爾氏菌烈性噬菌體，這對防控由青枯雷爾氏菌引起的細菌性病害具有較強的應用意義。蠟狀芽孢桿菌（Bacillus cereus）可導致人體腹瀉、嘔吐等食源性中毒，Hock等[33]從土壤中分離出一株感染蠟狀芽孢桿菌的長尾病毒科噬菌體Deep-Purple，隨后對該噬菌體進行全基因組測序、預測潛在編碼序列（coding sequence，CDS），并對CDS進行功能注釋。結果表明，CDS可分為結構相關蛋白（structural related proteins）、DNA復制轉錄（DNA replication and transcription）、DNA包裝（DNA packaging）及宿主裂解（host lysis）四個功能組。同時發現，該噬菌體不存在編碼潛在毒力因子的基因，且具有熱穩定性、pH穩定性的特點，研究人員認為噬菌體Deep-Purple可作為治療由蠟狀芽孢桿菌引起的食物中毒的潛在藥物。綜上所述，全基因組測序下的土壤噬菌體研究不僅有助于了解土壤中新型噬菌體的結構和功能，在致病菌引起的細菌性病害治療方面也具有積極意義。

3 土壤宏病毒基因組研究進展

3.1 土壤環境中的宏病毒組研究

宏病毒組主要對種群結構、基因功能活性、病毒與宿主的協作關系以及與環境之間的聯系進行探究，這為土壤環境微生物群落的研究提供了有力支撐[9，34]。病毒識別及其功能基因注釋是其中一個重要環節，在完成數據質控、基因組裝、基因預測等過程后，將預測的病毒編碼基因與COG、eggNOG[35]等數據庫比對，獲得病毒功能信息。表1對近年來土壤病毒組研究中具有代表性的功能基因進行總結，以揭示病毒與宿主群落的相互作用機制及其在元素生物地球化學循環中的作用機理。

表1 典型土壤病毒基因組功能基因及作用Table 1 Typical functional genes of soil viral genome

Jin等[36]在中國廣西和海南三個紅樹林生境采集樣本，研究了紅樹林土壤病毒功能多樣性。通過將預測的病毒ORF與eggNOG數據庫進行比對，發現大多數ORF無法獲得注釋信息，但有138個ORF可以注釋到與碳水化合物活性酶（carbohydrate-active enzyme，CAZyme）相關的基因。隨后，經過CAZy數據庫的進一步注釋，顯示糖苷水解酶類（glycoside hydrolases）ORF最為豐富，其次是糖基轉移酶類（glycosyl transferases）、碳水化合物結合模塊（carbohydrate-binding modules）等。這表明病毒可通過復雜多糖的生物分解直接調控碳循環，揭示了病毒在有機碳分解中的重要作用。Bi等[37]在中國西南部地區采集了4個玉米根際土壤及4個非根際土壤樣本，探究農業生態系統中病毒的多樣性及其對潛在生物地球化學循環的影響機制。通過對土壤病毒組進行測序，將測序讀段組裝成237個長度為10 kb以上的重疊群，并通過VirSorter軟件鑒定病毒序列。對這些土壤病毒進行功能基因注釋，共有40個基因被鑒定為溶菌酶或幾丁質酶，可用于降解宿主細胞壁。此外，共鑒定出了48個ORF與碳水化合物活性酶相關，包括碳水化合物結合模塊、碳水化合物酯酶（carbohydrate esterases）及糖苷水解酶。表明病毒可能編碼改變宿主活性的輔助代謝基因，間接參與土壤碳元素的生物地球化學循環。值得注意的是，研究人員發現該農業土壤中病毒的大多數CAZyme基因與紅樹林土壤中的病毒并不同[36]，認為病毒編碼的酶可能具有環境特異性。Segobola等[38]通過宏病毒組技術對灌木土壤的病毒群落進行探究。經過病毒基因組測序、組裝及功能基因注釋，發現組裝后最長的重疊群近乎一個完整的噬菌體基因組。該基因組的基因15和基因16分別對應末端酶大亞基和小亞基基因，參與了噬菌體雙鏈DNA的裂解和包裝；基因34的翻譯產物被識別為假定的ERF超家族蛋白，可參與噬菌體基因組的重組；基因41的翻譯產物被鑒定為gp77蛋白，與分枝桿菌（Mycobacterium）噬菌體Che9d編碼的同源物有95%的相似性，該蛋白在噬菌體復制的早期起到關閉基因（shut-off genes）的作用。通過KEGG Orthology（KO）數據庫對宏病毒組進行分析，發現相關代謝蛋白（如碳水化合物代謝、氨基酸代謝和核苷酸代謝）識別率最高，表明土壤病毒很可能干擾宿主的新陳代謝。Liang等[39]從美國東南部農業區老成土（美國土壤分類系統的一個土綱）中提取病毒。通過對土壤表層（0～16 cm）和亞表層（55～92 cm）病毒組基因序列進行組裝、分類，隨后采用VIROME對ORF進行功能基因注釋，發現大量病毒組序列功能未知，50.3%的預測蛋白在所比對的數據庫中沒有顯著同源性；有35.4%～38.7%的ORF被注釋為與宿主代謝途徑相關，如細胞信號（cell signaling）、氧化磷酸化（oxidative phosphorylation）、遺傳信息處理（genetic information processing）和磷、蛋白質及碳水化合物的代謝（metabolisms of phosphorous，protein，and carbohydrates）等。此外，研究人員還在病毒組中發現了豐富的碳水化合物代謝（carbohydrate metabolism）基因，表明病毒可能參與土壤碳循環的調節。值得注意的是，在亞表層土壤病毒組中檢測到的功能蛋白（除參與氧化磷酸化的蛋白質）編碼基因豐度比表層土壤病毒體高16倍，表明亞表層土壤病毒雖密度較低，但可能與微生物介導的過程密切相關；與噬菌體感染循環（即噬菌體溶源、裂解循環和原噬菌體誘導）和噬菌體結構成分（如噬菌體衣殼）相關的蛋白標準化豐度（normalized abundances）在亞表層土壤病毒組中也更高，這可能與病毒宿主在貧瘠營養條件下的協同進化有關。由此可得，病毒在調控土壤中營養元素生物地球化學循環，調節宿主新陳代謝及微生物群落結構方面發揮重要作用。

3.2 特殊及極端陸地環境中宏病毒組研究

Gao等[40]于廣東某鉛鋅礦尾礦庫采集樣本，研究高度分層硫化尾礦中病毒群落組成和功能特征。通過將預測的病毒蛋白與eggNOG數據庫比對，進行了病毒基因組直系同源基因簇COG的注釋分析。發現地表尾礦由于存在大量古菌和古菌病毒，導致大多數COG注釋困難；反之，以細菌為主的深層病毒群落存在著大量與同化硫酸鹽還原、轉座酶、DNA復制、噬菌體整合酶和重組酶相關的COG。隨后研究者試圖確定病毒編碼的輔助代謝基因，發現深層病毒群落含有豐富的同化硫酸鹽還原輔助代謝基因，這有利于宿主利用硫酸鹽，進而促進病毒的復制和繁殖。Daly等[41]對水力壓裂井中病毒與宿主的相互作用動力學進行探究，施加應激源對菌株（Halanaerobium. congolenseWG8）進行原噬菌體誘導實驗，并對病毒核酸進行純化、測序、基因注釋。結果顯示，病毒基因組中存在與整合酶（integrase）、切除酶（excisionase）、假想蛋白以及轉座酶（transposase）相關的基因；同時發現，其中一個基因被注釋為重疊感染排斥蛋白（superinfection exclusion protein）基因，推測該基因的存在可能有助于維持宿主容納原噬菌體與細胞裂解的平衡。Bezuidt等[42]在南極土壤群落中使用VirSorter工具從宏基因組序列數據中組裝了793個重疊群。通過病毒基因組數據庫對其進行分類注釋，有645個contigs被定義為病毒，且560個被進一步劃分為有尾噬菌體目。隨后使用eggNOG數據庫進行功能分析，發現噬菌體具有促進宿主感染的基因，如編碼幾丁質酶的基因，可參與宿主生物膜的降解。此外，eggNOG 功能分析也揭示了增強噬菌體毒性基因的存在，其中編碼甲基轉移酶的基因最為豐富，該基因編碼的酶有利于噬菌體規避宿主的限制修飾（restriction-modification，RM）系統。這一結果表明，噬菌體可能在面臨進化壓力時發展出對宿主的規避機制。Adriaenssens等[43]對納米比亞沙漠巖石下生物宏病毒組進行探究，發現大多數病毒為有尾噬菌體目，其中長尾病毒科是最常見的病毒型。通過功能基因注釋，有3%的基因被MG-RAST分類為“毒力、疾病、防御”（virulence，disease and defense）子系統，相應的編碼序列被鑒定為來自致病細菌的假想蛋白或噬菌體相關蛋白，如整合酶和復制蛋白（replication proteins）。此外，研究人員還發現噬菌體中存在與磷酸鹽調節相關的輔助代謝基因phoH，并通過MetaVir鑒定出18條完整的phoH基因和23條部分phoH序列。而在用極大或然法構建的系統發生樹卻顯示大部分沙漠巖石下病毒phoH氨基酸序列與NCBI數據庫中的完整噬菌體基因組序列關系較遠，且海洋和沙漠巖石下病毒phoH基因分布在不同的進化枝上，表明納米比亞沙漠巖石中存在獨特的噬菌體phoH基因序列。值得注意的是，該病毒樣本中沒有發現與光合作用或營養脅迫相關的宿主衍生基因，推測此生境中phoH基因較其他輔助代謝基因發揮更重要的作用。Emerson等[44]在瑞典斯托達倫沼澤地的泥炭巖芯處，采集了三個不同棲息地的病毒樣本，從中獲得了53個vOTUs（viral operational taxonomic units），發現僅有約30%的基因可被注釋，該結果可以佐證土壤是大量未知病毒遺傳多樣性的儲存庫[45]。此外，在13個vOTUs中鑒定出多個參與多糖結合（polysaccharide binding）、多糖降解（polysaccharide degradation）、中心碳代謝（central C metabolism）及孢子形成（sporulation）的輔助代謝基因。病毒的中心碳代謝基因可能在感染宿主期間增加核苷酸和能量產生；調節內孢子形成過程中的兩個輔助代謝基因spoVS和whiB分別有助于形成隔膜和外套，從而提高孢子的耐熱性。這揭示了病毒在介導碳代謝、土壤有機質降解、多糖結合和孢子形成過程中的調控作用。

上述案例揭示病毒不僅在農業土壤中與宿主、環境之間存在復雜的相互作用聯系，而且在極端、特殊的陸地環境中對調控微生物群落組成、影響生物地球化學循環、促進生物協同進化等方面也具有巨大潛能。

4 土壤宏病毒基因組的前沿與展望

現階段學術界對土壤病毒及其功能基因的科學認知依然十分有限。基于病毒組學領域的發展趨勢，今后土壤病毒組研究方向主要聚焦在以下方面：

1）土壤病毒主要通過微孔濾膜過濾的方式進行富集，該方法易將較大的病毒類型（如最近發現的巨型噬菌體）屏蔽在外[46-47]，從而缺乏對此類病毒功能基因的認識；此外，目前的研究主要聚焦于DNA 病毒，對RNA病毒研究較少。未來需研發針對巨型病毒的提取和富集技術、關注土壤RNA病毒的基因功能，這有助于學者探明土壤整體病毒的生態功能及作用機制。

2）現階段土壤病毒的提取、宏病毒組分析等缺乏統一技術規范。開發宏病毒組研究獨有的新方法，逐步規范技術流程，制定可以廣泛通用于土壤宏病毒組分析的技術導則和標準十分必要。

3）目前病毒功能基因注釋較為費時費力，即使通過自動注釋也因其準確性不足，需后期人工注釋進行修正，因此未來需開發更加高效、準確的生物信息學工具，識別病毒并注釋其基因組中的基因功能。

4）土壤病毒學的研究還處于起步階段，病毒基因組測序的數量遠落后于相應宿主細菌的基因組測序數量，通過同源蛋白對土壤病毒基因組上的蛋白功能進行注釋時，由于病毒ORF更短、進化更快，以及全球病毒取樣的有限性、數據庫較小等原因，導致部分ORF找不到匹配的功能注釋。故仍需大力發展全基因組擴增技術和測序技術，不斷完善土壤病毒資源庫，為病毒功能基因注釋提供有力支撐。

5）關注土壤病毒群落與宿主菌群的生態關系，進一步探究土壤病毒群落在元素生物地球化學循環中的直接與間接調控作用，深入探明土壤病毒與污染物的響應機制。

致 謝 衷心感謝美國Rice大學土木與環境工程系俞萍鋒博士在本文撰寫和修改過程中給予的學術指導和建議。