吳祎程，冉 濤，周傳社*，譚支良

(1.中國(guó)科學(xué)院亞熱帶農(nóng)業(yè)生態(tài)研究所 亞熱帶農(nóng)業(yè)生態(tài)過(guò)程重點(diǎn)實(shí)驗(yàn)室 畜禽養(yǎng)殖污染控制與資源化技術(shù)國(guó)家工程實(shí)驗(yàn)室 湖南省動(dòng)物營(yíng)養(yǎng)與生理代謝實(shí)驗(yàn)室 農(nóng)業(yè)部中南動(dòng)物營(yíng)養(yǎng)與飼料科學(xué)觀測(cè)站，長(zhǎng)沙 410125； 2. 中國(guó)科學(xué)院大學(xué)，北京 100049；3.蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院，蘭州 730020)

病毒(virus)廣泛存在于地球各種環(huán)境中，是地球上分布最廣、數(shù)量最為巨大的生命實(shí)體[1]。研究發(fā)現(xiàn)，大部分病毒能特異性侵染細(xì)菌，這些病毒被稱作噬菌體(bacteriophage或phage)[2-3]。在噬菌體被發(fā)現(xiàn)后不久，研究人員便證實(shí)了反芻動(dòng)物瘤胃內(nèi)也存在噬菌體[4-5]。長(zhǎng)期以來(lái)，盡管研究人員意識(shí)到病毒可能在反芻動(dòng)物瘤胃微生態(tài)系統(tǒng)中發(fā)揮著重要作用，但受限于技術(shù)手段和瘤胃環(huán)境的復(fù)雜性，瘤胃環(huán)境中噬菌體研究進(jìn)展相對(duì)緩慢[6]。目前，對(duì)瘤胃噬菌體的認(rèn)知多來(lái)自于傳統(tǒng)的培養(yǎng)、顯微鏡觀察和單一噬菌體的基因組學(xué)研究[7-8]，對(duì)其在瘤胃微生態(tài)系統(tǒng)層面的認(rèn)識(shí)仍十分匱乏且片面，以致至今對(duì)其在瘤胃微生態(tài)系統(tǒng)中的組成和功能仍知之甚少。近年來(lái)，隨著宏基因組學(xué)技術(shù)的發(fā)展，極大地推動(dòng)了海洋生態(tài)環(huán)境中病毒生態(tài)學(xué)研究的發(fā)展，并取得了引人注目的成績(jī)[9]。宏基因組學(xué)也被逐步運(yùn)用于胃腸道病毒組(virome)的研究，同時(shí)刷新了研究人員對(duì)胃腸道噬菌體的復(fù)雜性和多樣性的認(rèn)知[10]。Manrique等[11]首次利用病毒宏基因組學(xué)技術(shù)發(fā)現(xiàn)腸道中81%～93%的噬菌體為新型噬菌體，由于病毒缺乏合適的標(biāo)記基因和分析基準(zhǔn)，它們無(wú)法被歸類或找到對(duì)應(yīng)宿主。瘤胃微生態(tài)環(huán)境復(fù)雜，瘤胃作為一個(gè)巨大的基因資源庫(kù)，具有巨大的新病毒發(fā)掘空間，病毒宏基因組學(xué)技術(shù)與傳統(tǒng)技術(shù)的聯(lián)用在這一方面具有獨(dú)特優(yōu)勢(shì)，進(jìn)而推動(dòng)瘤胃噬菌體相關(guān)研究。

本文主要聚焦反芻動(dòng)物瘤胃生態(tài)系統(tǒng)中的噬菌體群落，并討論瘤胃噬菌體研究方法及現(xiàn)狀，及從其他生境中的噬菌體研究所得到的啟發(fā)。本文還展望了病毒宏基因組學(xué)在瘤胃環(huán)境中的應(yīng)用前景，以期為后續(xù)瘤胃噬菌體組研究提供科學(xué)參考。

1 瘤胃噬菌體概述

反芻家畜瘤胃內(nèi)棲息著數(shù)以萬(wàn)億計(jì)的微生物[12-15](圖1)，包括細(xì)菌、真菌、原蟲、古生菌和病毒，其中病毒的數(shù)量(>108～109pfu·mL-1)僅次于細(xì)菌，且主要為噬菌體(bacteriophage或phage)[16]。瘤胃內(nèi)的噬菌體具有不同形態(tài)，現(xiàn)有研究表明,反芻動(dòng)物瘤胃中，有尾噬菌體目(Caudovirales)普遍存在，包括肌尾噬菌體科(Myoviridae)、長(zhǎng)尾噬菌體科(Siphoviridae)、短尾噬菌體科(Podpviridae)占主導(dǎo)地位[17-18]，以及最近被鑒定出來(lái)的阿克曼噬菌體科(Ackermannviridae)和赫雷爾噬菌體科(Herelleviridae)[19]。其實(shí)，早在20世紀(jì)60年代，“瘤胃噬菌體”一詞就已進(jìn)入人們視野，但對(duì)其在瘤胃微生態(tài)系統(tǒng)中的功能知之甚少，直到20世紀(jì)80年代末，研究人員才開(kāi)始對(duì)其進(jìn)行探索[20]。隨著測(cè)序技術(shù)的發(fā)展，宏病毒組學(xué)技術(shù)日臻完善，借助該技術(shù)手段能夠更好地探索噬菌體這個(gè)“暗物質(zhì)”。已有研究表明，噬菌體在塑造瘤胃菌群結(jié)構(gòu)、維持微生物多樣性和調(diào)節(jié)宿主菌代謝等方面發(fā)揮著重要作用[21-22]。對(duì)反芻動(dòng)物瘤胃噬菌體多樣性的研究表明，噬菌體可能是通過(guò)裂解宿主細(xì)菌從而影響瘤胃細(xì)菌種群動(dòng)態(tài)變化[7,23-24]。Anderson等[23]通過(guò)營(yíng)養(yǎng)調(diào)控試驗(yàn)首次報(bào)道了飼糧變化對(duì)反芻動(dòng)物瘤胃噬菌體群落的影響，并表明，噬菌體可通過(guò)其編碼的輔助代謝基因(auxiliary metabolic genes, AMGs)，經(jīng)過(guò)裂解宿主細(xì)菌、產(chǎn)能、復(fù)制和微生物代謝的重新編程等一系列過(guò)程來(lái)影響瘤胃微生態(tài)系統(tǒng)[17,21,23]，而其中宿主菌代謝與反芻動(dòng)物飼料轉(zhuǎn)化效率密切相關(guān)。眾所周知，瘤胃微生物群體結(jié)構(gòu)的穩(wěn)定性和多樣性對(duì)于反芻家畜的健康、營(yíng)養(yǎng)、免疫和生存至關(guān)重要，因此，需要重視瘤胃噬菌體在整個(gè)瘤胃微生態(tài)系統(tǒng)中的地位、功能及作用。

2 瘤胃噬菌體研究方法

研究噬菌體的傳統(tǒng)方法是顯微鏡觀察和體外分離培養(yǎng)，但在實(shí)驗(yàn)室環(huán)境中培養(yǎng)的宿主菌數(shù)量有限；噬菌體還具有形態(tài)小、遺傳進(jìn)化快等特點(diǎn)，缺乏細(xì)菌所擁有的相對(duì)保守的系統(tǒng)發(fā)育標(biāo)記，如16S rRNA基因；且感染不同宿主細(xì)菌的噬菌體具有高度特異性，不同噬菌體序列的重疊群非常小[11,25-26]。在瘤胃噬菌體的早期研究中，利用透射電子顯微鏡(transmission electron microscope, TEM)技術(shù)，研究人員發(fā)現(xiàn)，瘤胃液中存在著形態(tài)各異、高度多樣化的噬菌體種群，并猜測(cè)其可能影響瘤胃內(nèi)細(xì)菌的種群[27]。后來(lái)，隨著DNA分析技術(shù)、脈沖場(chǎng)凝膠電泳技術(shù)(pulsed field gel electrophoresis, PFGE)和微生物分離鑒定技術(shù)的發(fā)展，研究人員證實(shí)，反芻動(dòng)物瘤胃液中存在大量溶菌性噬菌體[7]，并陸續(xù)分離鑒定了多種，包括以壞死梭桿菌(Fusobacteriumnecrophorum)[28]、白色瘤胃球菌(Ruminococcusalbus)[29]、溶纖維丁酸弧菌(Butyrivibriofibrisolvens)[30]和瘤胃擬桿菌(Prevotellaruminicola)[31]等作為宿主菌的噬菌體。TEM技術(shù)也隨之被廣泛用于新型噬菌體分離株的形態(tài)學(xué)表征[25, 32]。如今，宏基因組學(xué)測(cè)序技術(shù)和生物信息學(xué)的快速發(fā)展，使研究人員能更深入地了解病毒群落的結(jié)構(gòu)和功能。宏基因組學(xué)(metagenomics)手段是以某一特定環(huán)境樣本中的微生物群體基因組為研究對(duì)象，以功能基因篩選和序列測(cè)定分析為研究手段，以微生物多樣性、種群結(jié)構(gòu)、進(jìn)化關(guān)系、功能活性、相互作用以及其與環(huán)境間的關(guān)系為研究目的的一種新的微生物研究方法[33]。而宏病毒組學(xué)(viral metagenomics)則是結(jié)合病毒自身特點(diǎn)，將宏基因組學(xué)方法應(yīng)用于病毒學(xué)領(lǐng)域[34]。

瘤胃噬菌體宏基因組學(xué)研究，主要由病毒樣顆粒(virus-like particles, VLPs)濃縮液制備與生物信息學(xué)分析兩部分組成(圖2)。具體步驟包括：1)瘤胃樣本采集，經(jīng)過(guò)濾、消化以及高速離心等技術(shù)，去除細(xì)菌等其他潛在宿主；2)對(duì)樣品進(jìn)行富集、純化，得到病毒濃縮液；3)根據(jù)病毒核酸類型，選擇合適的病毒核酸提取試劑盒或手工提取的方式獲取病毒核酸，構(gòu)建病毒基因組文庫(kù)，并通過(guò)高通量測(cè)序技術(shù)進(jìn)行宏病毒組測(cè)序；4)對(duì)測(cè)序所得原始數(shù)據(jù)進(jìn)行質(zhì)控，基于重疊區(qū)(overlap)將高質(zhì)量測(cè)序讀段(reads)拼接為重疊群(contigs)；5)通過(guò)病毒序列鑒定軟件在重疊群中鑒別、篩選病毒序列，從而了解該環(huán)境中病毒構(gòu)成；6)對(duì)病毒基因組開(kāi)放閱讀框(open reading frame, ORF)進(jìn)行預(yù)測(cè)，再通過(guò)注釋工具將 ORF 與多個(gè)數(shù)據(jù)庫(kù)(KEGG、COG、CAZy、CARD等)比對(duì)進(jìn)行功能基因注釋。

圖2 瘤胃病毒組分析流程圖Fig.2 Flowcharts of rumen virome analysis

2.1 樣品前處理及噬菌體核酸提取

瘤胃噬菌體即為那些以瘤胃細(xì)菌為宿主菌的噬菌體，或是那些存在于反芻動(dòng)物瘤胃體內(nèi)的噬菌體[31]。由于采樣后，細(xì)菌和噬菌體仍然彼此接觸，噬菌體的侵染作用仍持續(xù)存在，長(zhǎng)時(shí)間孵育會(huì)影響噬菌體與微生物比例[35]。因此，樣品采集后要立即超低溫凍存或立即對(duì)樣品進(jìn)行處理，將病毒樣顆粒從樣品中分離并純化，用于病毒遺傳物質(zhì)提取?，F(xiàn)有的提取瘤胃噬菌體核酸的方法大多參考人或小鼠腸道噬菌體及土壤噬菌體分離方法[17,23,31]，截至目前，僅有一篇文獻(xiàn)以山羊和綿羊瘤胃液作為瘤胃噬菌體來(lái)源，建立了瘤胃噬菌體富集及DNA提取方法(包括過(guò)濾、離心、沉淀及核酸提取一系列過(guò)程)[18]。

瘤胃環(huán)境異常復(fù)雜且微生物多樣性高，瘤胃樣品中噬菌體的分離質(zhì)量直接關(guān)系到病毒宏基因組數(shù)據(jù)的質(zhì)量。Hungate[36]研究表明，大約75%的瘤胃細(xì)菌緊密附著在飼料顆粒上，以細(xì)菌為宿主的噬菌體也會(huì)附著在飼料上，因此，研究瘤胃噬菌體時(shí)，需同時(shí)采集固相和液相樣品，以便充分囊括瘤胃內(nèi)的所有噬菌體。在從樣品中分離噬菌體時(shí)，常利用噬菌體與其他微生物顆粒大小不同這一特征，借助微孔濾膜、離心等方式來(lái)實(shí)現(xiàn)噬菌體與其他微生物的分離。但是，實(shí)際操作中噬菌體樣極易混入細(xì)菌、真菌等微生物，尤其是那些與噬菌體顆粒大小相近的微生物，使得提取的病毒宏基因組中?；煊屑?xì)菌、真菌等其他微生物的核酸[37]。于是，在制備噬菌體粗制樣時(shí)，微孔濾膜、化學(xué)試劑等的選擇對(duì)去除細(xì)菌等其他微生物的污染和保證噬菌體群體的代表性具有重要影響，進(jìn)而對(duì)后續(xù)病毒核酸純度、基因功能注釋的準(zhǔn)確性產(chǎn)生影響。Concei??o-Neto等[38]研究發(fā)現(xiàn)，常見(jiàn)代表性病毒的大小由17到1 000 nm不等，其中，常見(jiàn)噬菌體的大小為125 nm左右。胃腸道病毒顆粒中絕大部分是噬菌體，為了避免細(xì)菌的污染，在研究胃腸道噬菌體組時(shí)，常選擇孔徑為0.22或0.45 μm的微孔濾膜，當(dāng)樣本環(huán)境復(fù)雜時(shí)，為避免遺漏巨型病毒信息，常選擇0.45 μm的微孔濾膜，但為了更有效地防止噬菌體宿主菌污染，0.22 μm的微孔濾膜仍被廣泛結(jié)合使用[39]。Friedersdorff等[30]選擇0.45和0.22 μm兩種孔徑的微孔濾膜依次過(guò)濾瘤胃液樣品，最終分離得到5種瘤胃溶菌性噬菌體，并對(duì)其進(jìn)行了全基因組學(xué)測(cè)序。同樣，Namonyo等[18]在建立瘤胃病毒DNA提取方法時(shí)也選用了0.45和0.22 μm兩種孔徑的微孔濾膜依次進(jìn)行過(guò)濾。而Klieve等[31]只選擇了孔徑為0.45 μm的微孔濾膜對(duì)瘤胃相關(guān)樣本進(jìn)行過(guò)濾。雖然現(xiàn)有文獻(xiàn)表明暫未鑒定出巨型瘤胃噬菌體，但并不排除技術(shù)手段造成的誤差，為了更全面地研究瘤胃環(huán)境噬菌體群落，研究人員仍需根據(jù)研究目的謹(jǐn)慎選擇微孔濾膜的孔徑。

在獲得含噬菌體的粗制樣后，常利用聚乙二醇(polyethylene glycol，PEG)配合不同濃度的NaCl來(lái)沉淀噬菌體顆粒以實(shí)現(xiàn)噬菌體的濃縮。較為常用的PEG和NaCl用量為25% PEG6000(w/v)和1.0 mol·L-1NaCl[18]；10% PEG8000和0.5 mol·L-1NaCl[31]及20% PEG6000和25 mol·L-1NaCl[30]。添加PEG后，常需在低溫下(4 ℃)孵育過(guò)夜，以使噬菌體充分沉淀，經(jīng)過(guò)高速離心(13 000×g, 30 min, 4 ℃； 12 000×g, 30 min, 4 ℃或52 350×g, 10 min, 4 ℃)后，倒掉上清液，獲得的沉淀即為噬菌體樣品。根據(jù)試驗(yàn)需求，該樣品可以直接用于后續(xù)試驗(yàn)，也可經(jīng)超速離心進(jìn)一步純化后用于后續(xù)試驗(yàn)。較為常用的是氯化銫(CsCl)密度梯度離心等，此法可對(duì)特定密度范圍內(nèi)的噬菌體進(jìn)行純化，并根據(jù)密度將VLPs分層。Carroll-Portillo等[40]的測(cè)序結(jié)果表明，用氯化銫密度梯度離心法富集的VLPs樣品宿主細(xì)菌核酸去除率比未經(jīng)富集的高，但氯化銫密度梯度離心法在富集樣品時(shí)，對(duì)不同密度的噬菌體有明顯的偏好性。Cordova等[41]也表明，某些噬菌體會(huì)由于氯化銫產(chǎn)生的滲透壓脅迫(osmotic shock)而導(dǎo)致噬菌體核酸丟失，因此，研究人員會(huì)選擇性采用氯化銫密度梯度離心法。

與瘤胃微生物組研究類似，病毒組研究也需先從分離得到的樣本中提取遺傳物質(zhì)。鑒于DNA病毒在瘤胃環(huán)境中的優(yōu)勢(shì)地位，目前，對(duì)于瘤胃噬菌體的研究主要集中于DNA病毒的核酸提取[20]。噬菌體核酸提取主要有兩種方式：第一種是試劑盒法，如Roche high pure viral nucleic acid large volume kit (Roche, Switzerland)[18]、QIAamp Ultra Sens Virus Kit (Qiagen)[23]和Fast DNATMSpin Kit for Soil (MP Biomedicals, Solon, OH, United States)[31]等，這類試劑盒可用于DNA病毒核酸的提取，試劑盒具有操作方便、提取純度高且不需要直接接觸有害化學(xué)試劑等優(yōu)點(diǎn)；第二種是傳統(tǒng)的苯酚-氯仿提取法，利用核酸、蛋白質(zhì)等雜質(zhì)在水相和有機(jī)相中溶解性不同而重新分配的性質(zhì)來(lái)實(shí)現(xiàn)核酸提取，該方法價(jià)格低廉，但耗時(shí)長(zhǎng)且提取濃度低，常需進(jìn)一步純化后才能達(dá)到測(cè)序要求[42]。試驗(yàn)中，若遇提取的噬菌體DNA濃度過(guò)低而不能滿足測(cè)序要求的情況，可利用多重置換擴(kuò)增(multiple displacement amplification, MDA)對(duì)核酸進(jìn)行擴(kuò)增，得到的DNA可用于構(gòu)建基因文庫(kù)或高通量測(cè)序[43]，但擴(kuò)增技術(shù)所得到的產(chǎn)物可能存在偏好性，還是建議提高原始噬菌體樣品富集量以保證結(jié)果的準(zhǔn)確性[3]。

2.2 高通量測(cè)序與病毒序列鑒定

高通量測(cè)序技術(shù)(high-throughput sequencing, HTS) 是對(duì)傳統(tǒng)Sanger測(cè)序技術(shù)革命性的變革，可對(duì)數(shù)百萬(wàn)個(gè)DNA分子進(jìn)行同時(shí)測(cè)序，并可深入地對(duì)一個(gè)物種的基因組和轉(zhuǎn)錄組進(jìn)行整體分析，因此，也稱其為下一代測(cè)序技術(shù) (next generation sequencing, NGS)[44]。HTS具有通量高、速度快、成本低等優(yōu)點(diǎn)，已廣泛應(yīng)用于宏病毒組學(xué)的研究中。近年來(lái)，牛津納米孔技術(shù)(Oxford nanopore technology, ONT)作為新興的單分子實(shí)時(shí)測(cè)序技術(shù)(single molecule real-time sequencing technology, SMRT)之一，具有快速制備文庫(kù)，超長(zhǎng)讀取和實(shí)時(shí)數(shù)據(jù)采集等優(yōu)勢(shì)[45]。該技術(shù)的核心是蛋白質(zhì)納米孔，通過(guò)產(chǎn)生覆蓋單個(gè)病毒顆粒內(nèi)所有突變的基因組長(zhǎng)度讀數(shù)來(lái)獲得病毒基因組。在此基礎(chǔ)上，中國(guó)科學(xué)院微生物研究所王軍研究員團(tuán)隊(duì)開(kāi)發(fā)了一種新的工作流程(包括病毒顆粒富集、核酸的逆轉(zhuǎn)錄和擴(kuò)增以及生物信息學(xué)分析)，并首次使用ONT PromethlON平臺(tái)對(duì)人類病毒組進(jìn)行表征[46]。相較于傳統(tǒng)的二代測(cè)序，ONT具有讀長(zhǎng)長(zhǎng)的優(yōu)點(diǎn)，因此，無(wú)需進(jìn)行PCR擴(kuò)增，便可以在更短時(shí)間內(nèi)、更準(zhǔn)確地從樣品中生成數(shù)據(jù)。

高通量測(cè)序后首先要對(duì)原始數(shù)據(jù)進(jìn)行質(zhì)量控制和過(guò)濾,去除宿主核酸，再進(jìn)行序列的拼接，并運(yùn)用專業(yè)軟件來(lái)鑒定病毒序列，確定該序列的生物群落來(lái)源，篩選出有用的基因信息。對(duì)噬菌體基因組進(jìn)行組裝時(shí)，一般挑選1 000 bp以上組裝序列，再應(yīng)用一系列軟件對(duì)病毒contig進(jìn)行鑒定。表1列出了常用于噬菌體鑒定的生物信息學(xué)分析工具。

從同時(shí)包含噬菌體及其宿主的混合基因組數(shù)據(jù)集中準(zhǔn)確識(shí)別噬菌體，是分析樣品中噬菌體組成的關(guān)鍵步驟。目前，有非常多基于序列或基于重疊群對(duì)病毒進(jìn)行鑒定的生物信息學(xué)分析軟件。VirSorter(https://github.com/simroux/VirSorter)可利用概率模型和大量病毒數(shù)據(jù)最大限度地檢測(cè)新病毒，雖然該工具在很大程度上依賴對(duì)已有病毒基因組的相似性搜索，但它有一個(gè)顯著優(yōu)勢(shì)，即可手動(dòng)管理病毒參考基因組數(shù)據(jù)庫(kù)，并補(bǔ)充了以海水、人體腸道、肺部組織及唾液樣本為來(lái)源的宏病毒組序列[47]。2017年，Ren等[48]開(kāi)發(fā)的VirFinder(https://github.com/jessieren/VirFinder)是通過(guò)利用細(xì)菌和病毒在K-mer上的差異，將病毒從宏基因組數(shù)據(jù)集中識(shí)別出來(lái)，與VirSorter相比，它能鑒定出更多潛在病毒，尤其是大小在1 000 bp以下的病毒contig。但VirSorter和VirFinder進(jìn)行噬菌體分析的功能仍相對(duì)有限，在鑒定噬菌體重疊群后，無(wú)法預(yù)測(cè)噬菌體與宿主間的相互作用。MARVEL(https://github.com/LaboratorioBioinformatica/MARVEL)是2018年由Amgarten等[49]開(kāi)發(fā)的工具，以其對(duì)Caudovirales目dsDNA病毒的有效注釋能力而為人所知，但該工具鑒定RNA病毒的錯(cuò)報(bào)率非常高。2019年開(kāi)發(fā)的VirMiner(https://github.com/TingtZHENG/VirMiner)則針對(duì)噬菌體分析功能做出了改進(jìn)，可通過(guò)預(yù)測(cè)模型從宏基因組數(shù)據(jù)集中預(yù)測(cè)噬菌體重疊群并進(jìn)行下游分析，包括功能基因注釋及噬菌體與宿主間關(guān)系預(yù)測(cè)等[50]。VIBRANT(https://github.com/AnantharamanLab/VIBRANT)，是2020年開(kāi)發(fā)的利用迭代注釋進(jìn)行病毒識(shí)別的工具，它能高效識(shí)別以細(xì)菌為宿主的多種病毒(包括dsDNA、ssDNA、dsRNA和ssRNA病毒)，并從宿主序列中提取前病毒區(qū)域，從而預(yù)測(cè)到其他分析工具無(wú)法預(yù)測(cè)到的前病毒(prophage)，且能解析環(huán)境中不同病毒之間新陳代謝的能力[51]。2021年，Guo等[52]改進(jìn)并發(fā)布了VirSorter2(https://github.com/jiarong/VirSorter2)，但VirSorter仍是最被廣泛使用的工具?，F(xiàn)有的每種工具都有其局限性，應(yīng)根據(jù)實(shí)際需求選擇使用或聯(lián)合使用不同的軟件。

表1 病毒鑒定工具

2.3 病毒功能基因注釋

使用BLAST軟件進(jìn)行基因預(yù)測(cè)，并將獲得的非冗余基因組與RefSeq病毒數(shù)據(jù)庫(kù)的參考序列進(jìn)行比對(duì)，獲取該基因的功能信息[53]。通過(guò)對(duì)病毒功能基因注釋，能夠?qū)ι钊胝J(rèn)識(shí)病毒個(gè)體生命過(guò)程提供理論基礎(chǔ)，還有助于了解病毒群落的生態(tài)過(guò)程及與宿主群落的互作網(wǎng)絡(luò)，從而闡釋病毒與宿主間復(fù)雜的相互作用機(jī)制。

京都基因與基因組百科全書(Kyoto encyclopedia of genes and genomes，KEGG)是一個(gè)系統(tǒng)分析基因功能的知識(shí)庫(kù)，KEGG具有強(qiáng)大的圖形功能，可利用圖形來(lái)介紹眾多代謝途徑(包括各種代謝通路、合成通路、膜轉(zhuǎn)運(yùn)、信號(hào)傳遞、細(xì)胞周期以及疾病相關(guān)通路等)以及各個(gè)途徑間的關(guān)系，這樣可以直觀地反映基因與相關(guān)代謝的關(guān)系[54]。蛋白質(zhì)直系同源簇(clusters of orthologous groups of proteins，COG)數(shù)據(jù)庫(kù)，可將不同物種中的直系同源基因進(jìn)行聚類，通過(guò)在不同物種中建立相關(guān)同源蛋白簇來(lái)預(yù)知位置蛋白質(zhì)的功能，同時(shí)也為分子系統(tǒng)發(fā)育分析提供數(shù)據(jù)基礎(chǔ)。2015年，Galperin等[55]基于完整噬菌體基因組中的編碼蛋白系統(tǒng)進(jìn)化關(guān)系構(gòu)建了 POG (phage orthologous groups)數(shù)據(jù)庫(kù)，可在雙鏈DNA噬菌體的全基因組序列中鑒定出保守的直系同源簇。Pfam 數(shù)據(jù)庫(kù)是蛋白質(zhì)家族的數(shù)據(jù)庫(kù)，根據(jù)多序列比對(duì)結(jié)果和隱馬爾可夫模型(hidden markov model, HMMs)，可查詢蛋白質(zhì)家族或蛋白結(jié)構(gòu)域的注釋、結(jié)構(gòu)及多序列比對(duì)信息，被廣泛用于基因功能注釋[56]?？股氐臑E用對(duì)畜牧業(yè)造成了深遠(yuǎn)影響，雖然抗生素對(duì)病毒并無(wú)直接作用，但有研究發(fā)現(xiàn)，瘤胃噬菌體攜帶抗生素耐藥性基因[23]，因此，對(duì)病毒功能基因序列進(jìn)行耐藥性檢測(cè)就顯得尤為重要?？股鼐C合研究數(shù)據(jù)庫(kù)(comprehensive antibiotic resistance database，CARD) 可在細(xì)菌耐藥性的分子基礎(chǔ)上，提供參考 DNA 和蛋白質(zhì)序列、檢測(cè)模型和生物信息學(xué)工具，通過(guò)與該數(shù)據(jù)庫(kù)進(jìn)行比對(duì)，可用于關(guān)聯(lián)抗生素模塊及其目標(biāo)、抗性機(jī)制、基因變異等信息與耐藥基因相關(guān)的注釋信息[57]。

3 瘤胃噬菌體基因組研究趨勢(shì)及進(jìn)展

借助Web of Science核心合集以“rumen bacteriophage”、“rumen phage”或“rumen virus”為關(guān)鍵詞進(jìn)行檢索，2000—2021年時(shí)間段內(nèi)共有108篇相關(guān)文獻(xiàn)，且發(fā)文量呈現(xiàn)整體上升趨勢(shì)(圖3a)，2021年發(fā)表的5篇文獻(xiàn)并未在圖中標(biāo)明。其中，美國(guó)、澳大利亞、印度、新西蘭和加拿大的發(fā)文量分別為33.66%、15.39%、9.62%、8.65%和7.69%，而我國(guó)在瘤胃噬菌體這一領(lǐng)域的研究嚴(yán)重滯后。所發(fā)表的關(guān)于瘤胃噬菌體的文章中，基于宏基因組學(xué)的文獻(xiàn)僅有7篇[17-18,21,23-24,58-59]，涉及研究動(dòng)物包括水牛、奶牛、綿羊、山羊和麋鹿，且很少涉及瘤胃噬菌體群落動(dòng)態(tài)變化規(guī)律及瘤胃噬菌體在整個(gè)瘤胃微生態(tài)系統(tǒng)中的作用。隨后，利用VOSviewer進(jìn)行關(guān)鍵詞密度聚類，并進(jìn)行可視化分析[60](圖3b)，圖中結(jié)點(diǎn)越接近黃色，表示此為研究熱點(diǎn)，可發(fā)現(xiàn)瘤胃環(huán)境病毒研究熱點(diǎn)為噬菌體及宏基因組層面。而以“gut phage”或“gut virus”作為關(guān)鍵詞檢索，在2000—2021年時(shí)間段內(nèi)共有4 449篇相關(guān)文獻(xiàn)，這說(shuō)明，瘤胃噬菌體是一個(gè)極具發(fā)掘空間的研究領(lǐng)域。

對(duì)病毒進(jìn)行全基因組測(cè)序(complete genome sequencing)，又叫單病毒基因組學(xué)(single-virus genomics, SVG)，可獲得特定病毒(噬菌體)的全部核酸序列。SVG即以純培養(yǎng)的病毒為研究對(duì)象，進(jìn)行核酸提取、文庫(kù)構(gòu)建并對(duì)單個(gè)病毒進(jìn)行全基因組測(cè)序，SVG從頭組裝可以規(guī)避病毒遺傳物質(zhì)多樣性的問(wèn)題，從而構(gòu)建更完整的病毒基因組，以獲得更全面的病毒數(shù)據(jù)庫(kù)的信息[61]。目前，僅有兩篇基于SVG來(lái)研究瘤胃噬菌體的文章，表2綜述了噬菌體及其宿主信息。

Friedersdorff等[30]以瘤胃擬桿菌(Prevotellaruminicola)、白色瘤胃球菌(Streptococcusbovis)及瘤胃鏈球菌(Ruminococcusalbus)為宿主，采用雙層平板法(plaque assay)從瘤胃相關(guān)樣品(廢水、牛糞及瘤胃液)中分離出5株烈性噬菌體，這是首次對(duì)瘤胃環(huán)境中以特定細(xì)菌為宿主的噬菌體的基因組研究。隨后對(duì)這5株噬菌體進(jìn)行全基因組測(cè)序，并通過(guò) Glimmer和Prodigal在線工具對(duì)噬菌體全基因組序列進(jìn)行潛在編碼序列(coding sequence, CDS)預(yù)測(cè)和功能注釋。發(fā)現(xiàn)在φBrb01和φBrb02基因組中，CDS 可分為DNA 復(fù)制轉(zhuǎn)錄(DNA replication and transcription)、DNA 包裝(DNA packaging)及宿主裂解(Host lysis)3個(gè)功能組。溶纖維丁酸弧菌(Butyrivibriofibrisolvens)是瘤胃內(nèi)的一種優(yōu)勢(shì)菌，對(duì)纖維降解和蛋白質(zhì)分解等有重要作用。Friedersdorff等[30]以其為宿主，從瘤胃液及糞便中分離了5株噬菌體，并進(jìn)行全基因組測(cè)序和利用 Glimmer和GeneMarkS在線工具進(jìn)行功能注釋。結(jié)果顯示，在科水平上，瘤胃噬菌體以有尾噬菌體為主，其中長(zhǎng)尾噬菌體科占主導(dǎo)地位，與前人研究結(jié)果一致[17-18]。Ceridwen屬噬菌體基因組中與其他噬菌體基因組同源的ORF中，其中20個(gè)與Siphoviridae科噬菌體有同源性，1個(gè)與Podoviridae科噬菌體有同源性。另外，用PHACTS軟件對(duì)噬菌體生活史進(jìn)行預(yù)測(cè)發(fā)現(xiàn)，這10株噬菌體中僅有Arawn和Idris屬噬菌體基因組中存在編碼整合酶的基因，表明這兩株噬菌體可能為溶源性噬菌體(lysogenic phage)，其他8株噬菌體為溶菌性噬菌體(lytic phage)。值得注意的一點(diǎn)是，從瘤胃液中分離出的Bo-Finn噬菌體與從牛糞便中分離出的Arian噬菌體基因組序列相似性可達(dá)98.6%，它們屬于同種噬菌體。

4 不同生境噬菌體研究進(jìn)展

自2002年美國(guó)學(xué)者Breitbart等[62]首次對(duì)海洋病毒群(virome)進(jìn)行宏基因組測(cè)序，到2020年，全球范圍內(nèi)的研究已報(bào)道了近20萬(wàn)個(gè)病毒種群[63]。目前，噬菌體的研究已然成為研究不同生境內(nèi)微生物群落結(jié)構(gòu)和功能的焦點(diǎn)，如哺乳動(dòng)物腸道噬菌體與機(jī)體健康和疾病發(fā)生已受到廣泛關(guān)注。關(guān)于不同生境噬菌體(如人體腸道、土壤等)的研究已證實(shí)，噬菌體是導(dǎo)致環(huán)境中宿主細(xì)菌死亡的主要原因[64]，且噬菌體可作為水平轉(zhuǎn)移基因(horizontal gene trans-fer, HGT)的重要載體，與宿主菌進(jìn)行遺傳物質(zhì)的交換，進(jìn)而調(diào)節(jié)細(xì)菌的進(jìn)化和多樣性；值得注意的是，HGT也是導(dǎo)致抗生素抗性基因(antibiotic resistance genes, ARGs)擴(kuò)散的重要因素[65]。鑒于瘤胃內(nèi)噬菌體研究相對(duì)滯后的現(xiàn)狀，可以適度借鑒其他環(huán)境樣本中的噬菌體研究方法并用于瘤胃噬菌體研究。

圖3 瘤胃噬菌體領(lǐng)域年發(fā)表論文數(shù)量(a)及關(guān)鍵詞密度圖譜(b)Fig.3 Number of publications(a) and keyword dendity network map(b) in the field of rumen virus (bacteriophage)

許多環(huán)境噬菌體的研究得益于微生物多組學(xué)方法聯(lián)用[66-67]，取得了喜人的成果。例如，Emerson等[66]利用宏基因組和宏轉(zhuǎn)錄組揭示了氣候變化中噬菌體對(duì)復(fù)雜碳降解的潛能。此外，Brum等[67]結(jié)合蛋白質(zhì)組學(xué)和宏基因組學(xué)，從海洋環(huán)境中鑒定出多樣性極高的病毒衣殼蛋白。類似的，蛋白質(zhì)組學(xué)這一新興領(lǐng)域在瘤胃微生態(tài)系統(tǒng)中的研究表明，瘤胃噬菌體介導(dǎo)的細(xì)胞裂解會(huì)釋放微生物胞內(nèi)酶，包括那些參與碳水化合物分解的微生物酶，這些酶的釋放會(huì)促進(jìn)瘤胃內(nèi)飼料的降解[59]。

除了以上能通過(guò)組學(xué)技術(shù)挖掘出的數(shù)據(jù)外，病毒宏基因組測(cè)序結(jié)果中大部分序列為未知序列，為了充分利用這些未知數(shù)據(jù)，Brum等[67]通過(guò)對(duì)預(yù)測(cè)的蛋白質(zhì)序列進(jìn)行全長(zhǎng)相似性聚類，揭示了七大洋海洋樣本中存在一個(gè)核心病毒基因集，這樣可以從未知序列中得到另一種“已知”結(jié)果。Swain等[22]考慮到瘤胃病毒在不同個(gè)體間存在的差異，運(yùn)用類似方法發(fā)現(xiàn)了核心瘤胃病毒基因集。此外，Allen等[68]研究海洋微生物群落時(shí)采用了一種基于噬菌體基因組標(biāo)簽的方法(phage genome signature-based recovery, PGSR)，從宏基因組中提取噬菌體序列，并結(jié)合病毒/細(xì)菌比(virus-to-bacteria ratio，VBR)來(lái)評(píng)估病毒對(duì)細(xì)菌群落的影響，該方法有助于研究人員了解噬菌體在塑造整個(gè)海洋微生物群落結(jié)構(gòu)中的作用。但目前未見(jiàn)此法運(yùn)用于瘤胃噬菌體研究。

表2 已報(bào)道的獲得全基因組信息的瘤胃噬菌體

值得注意的是，不同生態(tài)環(huán)境中研究方法的適用性有限。即使兩個(gè)生態(tài)系統(tǒng)具有非常相似的病毒群落結(jié)構(gòu)，他們潛在的微生態(tài)關(guān)系也不盡相同。比如，極地水體環(huán)境噬菌體和胃腸道環(huán)境噬菌體中，溶源性噬菌體較溶菌性噬菌體在數(shù)量上占有優(yōu)勢(shì)[69-71]。但極地海洋生態(tài)系統(tǒng)中，隨著細(xì)菌豐度增加，溶源性噬菌體會(huì)從溶源循環(huán)(lysogenic cycle)轉(zhuǎn)變?yōu)槿芫h(huán)(lytic cycle)[69]；而在腸道中，宿主菌豐度增加時(shí)，溶源性噬菌體會(huì)持續(xù)占領(lǐng)主導(dǎo)地位[70-71]。因此，在將不同生境中噬菌體研究方法及生物學(xué)概念外推到瘤胃環(huán)境時(shí)仍需謹(jǐn)慎。

5 瘤胃宏病毒組的前景與展望

隨著宏基因組學(xué)和生物信息學(xué)的發(fā)展，研究人員可以不再需要體外分離培養(yǎng)噬菌體，就能獲得一些分散、豐度低的病毒的結(jié)構(gòu)和功能信息，也逐步增進(jìn)對(duì)瘤胃內(nèi)噬菌體的多樣性及潛在功能的認(rèn)知。但是，現(xiàn)階段學(xué)術(shù)界對(duì)瘤胃噬菌體及其功能基因的科學(xué)認(rèn)知仍然十分有限，筆者結(jié)合自身體會(huì)，認(rèn)為今后此領(lǐng)域的研究應(yīng)主要聚焦在以下幾個(gè)方面：

5.1宏基因組學(xué)研究中，瘤胃噬菌體樣品制備時(shí)主要通過(guò)微孔濾膜過(guò)濾的方式進(jìn)行純化富集，因此，難免會(huì)因?yàn)闉V膜孔徑選擇過(guò)小而造成顆粒較大的噬菌體的遺傳信息的遺漏?，F(xiàn)階段瘤胃噬菌體的提取、宏病毒組分析等缺乏統(tǒng)一技術(shù)規(guī)范，建立成熟的噬菌體分離純化實(shí)驗(yàn)室標(biāo)準(zhǔn)化操作流程以及針對(duì)巨型噬菌體的優(yōu)化提取和富集技術(shù)是當(dāng)務(wù)之急。

5.2目前的研究主要聚焦在瘤胃DNA病毒，對(duì)于RNA病毒的研究較少。未來(lái)要利用多組學(xué)手段(如宏基因組、宏蛋白組、宏轉(zhuǎn)錄組和宏代謝組等)結(jié)合生物信息學(xué)技術(shù)，同步關(guān)注瘤胃RNA病毒的基因功能，盡可能獲得全面的噬菌體信息，這有助于探明瘤胃整體病毒的生態(tài)功能和作用機(jī)制。

5.3鑒于噬菌體種群在個(gè)體間差異顯著且噬菌體種群存在于反芻動(dòng)物整個(gè)消化道，需進(jìn)一步研究不同生理階段、胃腸道不同部位噬菌體群落對(duì)微生物群落的影響。

5.4盡管病毒宏基因組學(xué)為研究人員提供了一種不需分離培養(yǎng)即可了解環(huán)境中全部病毒基因組信息的方法，但仍需要以瘤胃優(yōu)勢(shì)菌群為宿主菌，對(duì)噬菌體進(jìn)行分離培養(yǎng)，這將有助于全面了解噬菌體的生物學(xué)特性以及它們?cè)诰S持瘤胃動(dòng)態(tài)平衡中的作用。

5.5關(guān)于瘤胃病毒的研究在國(guó)內(nèi)處于亟待開(kāi)發(fā)的階段，瘤胃是一個(gè)復(fù)雜、多變的微生態(tài)系統(tǒng)，蘊(yùn)藏著豐富的基因資源，有大量新病毒等待研究人員去挖掘。故需要大力發(fā)展全基因組擴(kuò)增技術(shù)和測(cè)序技術(shù)，不斷完善動(dòng)物胃腸道病毒資源庫(kù)，豐富瘤胃噬菌體基因組數(shù)據(jù)庫(kù)，揭示噬菌體塑造瘤胃細(xì)菌群體結(jié)構(gòu)的機(jī)制。

5.6我國(guó)于2020年開(kāi)始對(duì)飼料全面禁抗，鑒于瘤胃噬菌體對(duì)瘤胃微生物的裂解作用及對(duì)瘤胃營(yíng)養(yǎng)物質(zhì)循環(huán)的影響，噬菌體制劑作為抗生素替代品將成為畜牧界的重大課題。