最新的一項(xiàng)研究顯示，去年在云南勐臘縣獲得的蝙蝠樣本中發(fā)現(xiàn)一種新冠狀病毒，該病毒在最長(zhǎng)編碼基因區(qū)1ab與新冠病毒（SARS-CoV-2）一致性達(dá)到97.2%。同時(shí)，研究團(tuán)隊(duì)在這一新冠狀病毒中S蛋白（刺突糖蛋白）S1和S2交界處也發(fā)現(xiàn)插入了三個(gè)殘基PAA，他們由此認(rèn)為SARS-CoV-2在S1和S2的交界處的四氨基酸插入并非“人工痕跡”，有力證明了其起源于自然。

當(dāng)?shù)貢r(shí)間5月10日，《細(xì)胞》（Cell）旗下子刊《現(xiàn)代生物學(xué)》（Current Biology）在線發(fā)表了“A novel bat coronavirus closely related to SARS-CoV-2 contains natural insertions at the S1/S2 cleavage site of the spike protein ”，來(lái)自山東第一醫(yī)科大學(xué)、山東省高等學(xué)校新發(fā)傳染病病因流行病學(xué)實(shí)驗(yàn)室、中國(guó)科學(xué)院西雙版納熱帶植物園、中國(guó)科學(xué)院北京生命科學(xué)研究所、中國(guó)科學(xué)院武漢病毒研究所、中國(guó)科學(xué)院微生物研究所、澳大利亞悉尼大學(xué)聯(lián)合團(tuán)隊(duì)的研究人員給出了上述結(jié)論。

在這篇論文中，研究團(tuán)隊(duì)報(bào)告了一種新的蝙蝠來(lái)源冠狀病毒，命名為RmYN02，是在2019年5月至10月期間從中國(guó)云南省收集的227只蝙蝠的基因組分析中鑒定出來(lái)的。

研究顯示，RmYN02和新型冠狀病毒SARS-CoV-2在全病毒基因組中的同源性為93.3%，在與SARS-CoV-2最接近的1ab基因中的同源性為97.2%。相比之下，RmYN02受體結(jié)合域（RBD）與SARS-CoV-2之間呈現(xiàn)低序列同源性（61.3%），可能不與血管緊張素轉(zhuǎn)換酶2（ACE2）結(jié)合。

然而，與SARS-CoV-2相似，RmYN02的特點(diǎn)是在S蛋白的S1和S2亞單位的剪切位點(diǎn)處有多個(gè)氨基酸插入。研究團(tuán)隊(duì)認(rèn)為，這就有力證明了這種插入事件在自然界中是可以發(fā)生的。

這些數(shù)據(jù)表明，SARS-CoV-2源于蝙蝠和其他野生動(dòng)物中存在的病毒之間的多重自然重組。

該論文通訊作者為山東第一醫(yī)科大學(xué)基礎(chǔ)醫(yī)學(xué)院病原生物學(xué)研究所所長(zhǎng)、山東省高等學(xué)校新發(fā)傳染病病因流行病學(xué)實(shí)驗(yàn)室主任史衛(wèi)峰，中國(guó)科學(xué)院微生物研究所病毒傳播預(yù)警與致病機(jī)制研究組組長(zhǎng)、項(xiàng)目研究員畢玉海，中國(guó)科學(xué)院西雙版納熱帶植物園綜合保護(hù)中心景觀生態(tài)組組長(zhǎng)、副教授Alice Catherine Hughes。該研究此前于當(dāng)?shù)貢r(shí)間3月5日發(fā)表在預(yù)印本網(wǎng)站bioRxiv上。

新冠病毒起源尚不明確，野生動(dòng)物中仍有大量冠狀病毒

SARS-CoV-2在中國(guó)及其他地區(qū)造成了前所未有的肺炎流行，已引起全球范圍的公共衛(wèi)生關(guān)注。盡管蝙蝠被認(rèn)為是SARS-CoV-2 最有可能的自然宿主，但病毒的起源仍不清楚。

系統(tǒng)發(fā)育分析表明，SARS-CoV-2是一種不同于SARS-CoV和MERS-CoV的新型冠狀病毒。到目前為止，與SARS-CoV-2最密切相關(guān)的病毒是RaTG13，它是中科院武漢病毒所石正麗團(tuán)隊(duì)于2013年在云南采集的一只中華菊頭蝠標(biāo)本中分離發(fā)現(xiàn)。RaTG13病毒株和新冠病毒具有96.1%的核苷酸同源性、92.9%的S基因同源性。這些數(shù)據(jù)再次說(shuō)明蝙蝠是冠狀病毒的重要宿主。

然而，值得注意的是，香港大學(xué)公共衛(wèi)生學(xué)院新發(fā)傳染病國(guó)家重點(diǎn)實(shí)驗(yàn)室管軼教授、廣西醫(yī)科大學(xué)胡艷玲教授團(tuán)隊(duì)，以及華南農(nóng)業(yè)大學(xué)、嶺南現(xiàn)代農(nóng)業(yè)科學(xué)與技術(shù)廣東省實(shí)驗(yàn)室沈永義教授、肖立華教授團(tuán)隊(duì)兩個(gè)研究小組此前均報(bào)告了存在于馬來(lái)穿山甲的SARS-CoV-2相關(guān)冠狀病毒，這些穿山甲通過(guò)非法走私進(jìn)入廣西和廣東。

研究團(tuán)隊(duì)提到，盡管在全基因組水平上，這些穿山甲中檢測(cè)到的冠狀病毒與SARS-CoV-2的距離要遠(yuǎn)于RaTG13與SARS-CoV-2的距離，但它們?cè)赟蛋白的受體結(jié)合域(RBD)上與SARS-CoV-2非常相似。

因此，盡管目前還不清楚穿山甲是否是SARS-CoV-2傳播到人過(guò)程中的中間宿主，但它們可能在冠狀病毒的生態(tài)學(xué)和進(jìn)化中發(fā)揮重要作用。

他們認(rèn)為，這些在穿山甲中發(fā)現(xiàn)的病毒可以表明，在野生動(dòng)物中仍有大量的冠狀病毒樣本，其中一些可能直接參與了SARS-CoV-2的出現(xiàn)。

SARS-CoV-2和幾種代表性蝙蝠來(lái)源冠狀病毒序列比較。

一種新的蝙蝠來(lái)源冠狀病毒，命名RmYN02

論文中提到，2019年5月至10月,研究團(tuán)隊(duì)從云南勐臘縣的227只蝙蝠中一共收集了302個(gè)樣本。這些蝙蝠屬于20個(gè)不同種類，大多數(shù)樣本從馬來(lái)菊頭蝠Rhinolophus malayanus（n=48, 21.1%）、中蹄蝠Hipposideros larvatus （n=41，18.1%）和小褐菊頭蝠Rhinolophus stheno（n=39,17.2%)中獲得。樣本來(lái)源于多種組織，包括翼膜（219）、肺（2）、肝（3）和糞便（78）。

除了3只蝙蝠外，其余所有蝙蝠均為在活著時(shí)取樣并被釋放。

利用新一代宏基因組測(cè)序技術(shù)，研究團(tuán)隊(duì)首先鎖定了2個(gè)初步的一致序列。這些序列產(chǎn)生的樣本來(lái)自于2019年5月6日至7月30日期間的11份馬來(lái)菊頭蝠糞便。經(jīng)過(guò)一系列驗(yàn)證步驟，研究團(tuán)隊(duì)得到一個(gè)部分（23395bp）和一個(gè)完整（29671bp）的蝙蝠冠狀病毒基因組序列，分別命名為BetaCoV/Rm/Yunnan/YN01/2019 （RmYN01）和BetaCoV/Rm/Yunnan/YN02/2019 （RmYN02）。

相比之下，RmYN02與SARS-CoV-2密切相關(guān)，表現(xiàn)出93.3%的核苷酸序列一致性，但從全長(zhǎng)基因組層面來(lái)說(shuō)，RaTG13和SARS-CoV-2的一致性更高(96.1%)。RmYN02和SARS-CoV-2在大多數(shù)基因組區(qū)域(如1ab、3a、E、6、7a、N和10)非常相似（>96%序列一致性）。特別是，RmYN02在最長(zhǎng)編碼基因區(qū)1ab (n=21285)與SARS-CoV-2一致性達(dá)到97.2%。

不過(guò)，RmYN02和SARS-CoV-2在S基因中的序列一致性（核苷酸71.8%，氨基酸72.9%）遠(yuǎn)低于RaTG13和SARS-CoV-2之間的97.4%。另外值得注意的是，RmYN02和SARS-CoV-2在RBD中的氨基酸同源性僅為62.4%。而來(lái)自廣東的穿山甲冠狀病毒和SARS-CoV-2在RBD中的氨基酸同源性為97.4%，也是在RBD區(qū)域目前和SARS-CoV-2最接近的。

圍繞SARS-CoV-2的同源模型、體外實(shí)驗(yàn)和S蛋白的三維結(jié)構(gòu)分析結(jié)果都表明，SARS-CoV-2和SARS-CoV一樣，也可以利用ACE2作為細(xì)胞受體。研究團(tuán)隊(duì)也同樣使用同源模型分析了RmYN02、RaTG13和兩個(gè)穿山甲CoVs的RBD。

RmYN02和代表性冠狀病毒RBD結(jié)構(gòu)的同源建模和結(jié)構(gòu)比較。

研究發(fā)現(xiàn)，RmYN02 RBD中的氨基酸缺失在受體結(jié)合位點(diǎn)附近形成了兩個(gè)比SARS-CoV-2 RBD短的環(huán)。重要的是，在SARS-CoV、SARS-CoV-2、RaTG13、穿山甲/MP789/2019、穿山甲/GX/P5L/2017的外子域（external subdomain ）保守的二硫鍵在RmYN02中缺失了。研究團(tuán)隊(duì)推測(cè)，這些缺失可能導(dǎo)致構(gòu)象變化，從而減少RmYN02 RBD與ACE2的結(jié)合，甚至導(dǎo)致不結(jié)合。

當(dāng)然也有可能的一種情況是，包括RmYN02、ZXC21和ZC45在內(nèi)的環(huán)缺失的SARS相關(guān)冠狀病毒，使用了一種我們目前未知的受體。

值得一提的是，此前有研究認(rèn)為，RBD中的6個(gè)氨基酸殘基(L455、F486、Q493、S494、N501和Y505)是SARS-CoV-2與ACE2有效受體結(jié)合的主要決定因素。與同源建模一致，穿山甲/MP789/2019在所有6個(gè)位置上都具有與SARS-CoV-2相同的氨基酸殘基。相比之下，RaTG13、RmYN02和RmYN01與SARS-CoV-2均只在1個(gè)位置上有相同的氨基酸殘基。

研究團(tuán)隊(duì)認(rèn)為，這種進(jìn)化模式是重組和自然選擇的復(fù)雜結(jié)合的表現(xiàn)。

系統(tǒng)發(fā)育樹：A 全長(zhǎng)基因組；B、S基因；C、RBD ；D、RdRp(RNA依賴的RNA聚合酶

研究團(tuán)隊(duì)還對(duì)RmYN02、RaTG13、SARS-CoV-2和穿山甲中的蝙蝠冠狀病毒進(jìn)行了系統(tǒng)發(fā)育分析。與先前的研究相符，穿山甲beta-CoVs形成兩個(gè)亞型。然而，論文中提到，穿山甲是否是這些病毒的天然蓄水池，或者它們從蝙蝠或其他野生動(dòng)物中獨(dú)立獲得，這都需要進(jìn)一步驗(yàn)證。

更值得注意的是，在大多數(shù)病毒基因組中，RmYN02與SARS-CoV-2的親緣關(guān)系最近，盡管這兩種病毒之間仍然有一段較長(zhǎng)的分支距離。S基因樹顯示，SARS-CoV-2離RaTG13較近，和RmYN02較遠(yuǎn)，這表明后者在S基因經(jīng)歷了重組。在RBD的系統(tǒng)發(fā)育樹上，SARS-CoV-2和pangolin-CoV/GD最密切相關(guān)，與蝙蝠病毒都較遠(yuǎn)，再次表明重組發(fā)生。最后，完整的RNA依賴的RNA聚合酶(RdRp)基因(在RNA病毒的系統(tǒng)發(fā)育分析中經(jīng)常被使用)系統(tǒng)發(fā)育分析顯示，RmYN02、RaTG13和SARS-CoV-2形成了一個(gè)與穿山甲病毒完全不同的亞群。

SARS-CoV-2是自然起源，可能通過(guò)重組獲得

類似禽流感病毒(AIVs)血凝素(HA)蛋白的方式，冠狀病毒的S蛋白在功能上分裂成兩個(gè)亞基S1和S2。而在某些AIV亞型的剪切位點(diǎn)上插入多堿基氨基酸被認(rèn)為與增強(qiáng)致病性有關(guān)。

值得注意的是，SARS-CoV-2的特征之一即為在S1和S2的交界處有一個(gè)四氨基酸的插入，這在其他β冠狀病毒的其他譜系中沒(méi)有觀察到過(guò)。這種被稱為弗林蛋白酶剪切位點(diǎn)的插入是SARS-CoV-2獨(dú)有的，目前在所有檢測(cè)的SARS-CoV-2序列中都發(fā)現(xiàn)存在。

研究團(tuán)隊(duì)此番在RmYN02中S1和S2交界處也發(fā)現(xiàn)插入了三個(gè)殘基PAA，他們認(rèn)為這是非常重要的?！半m然SARS-CoV-2插入的殘基(以及由此產(chǎn)生的核苷酸)與RmYN02中的殘基不同，但可以表明它們是一個(gè)獨(dú)立的插入事件，它們?cè)谝吧鷦?dòng)物(蝙蝠)中出現(xiàn)強(qiáng)烈表明它們是自然起源的，可能是通過(guò)重組獲得的?！?

因此，這些數(shù)據(jù)有力地表明了SARS-CoV-2的自然起源。

此外，研究團(tuán)隊(duì)再次確定了RmYN02的蝙蝠宿主為馬來(lái)菊頭蝠（Rhinolophus malayanus），和一種馬來(lái)菊頭蝠標(biāo)本中獲得的序列（GenBank accession MK900703 ）100%一致。

馬來(lái)菊頭蝠和中華菊頭蝠都廣泛分布于中國(guó)西南部和東南亞地區(qū)。論文中提到，一般來(lái)說(shuō)，這些蝙蝠不會(huì)長(zhǎng)距離遷移，群居性很強(qiáng)，很可能生活在同一個(gè)洞穴中，這可能會(huì)促進(jìn)它們之間的病毒交換和重組的發(fā)生。

值得注意的是，科學(xué)家從肛門拭子中發(fā)現(xiàn)了RaTG13，從糞便中發(fā)現(xiàn)了RmYN02。因此，糞便是蝙蝠將病毒傳播給其他動(dòng)物，特別是能夠利用洞穴環(huán)境的物種的一種簡(jiǎn)單但可行的方法。