新冠病毒突變?yōu)椤案鼜?qiáng)”亞型？最新結(jié)論并不如此

防務(wù)

      新京報(bào)快訊（記者 李玉坤）《國家科學(xué)評論》3月3日在線發(fā)布了北大生命科學(xué)院研究員陸劍等人的論文，媒體對此論文解讀稱，新冠病毒已經(jīng)突變，演化為兩個(gè)亞型，其中一個(gè)亞型“傳染性”更強(qiáng)（aggressive）。
: c& S6 ]/ l3 {+ q/ y. H      3月11日晚，記者發(fā)現(xiàn)陸劍等人已經(jīng)重新提交了手稿，在文末解釋稱，“aggressive”（具有侵略性）用詞具有誤導(dǎo)性，他們沒有關(guān)于病毒致病性的結(jié)論。
6 {0 M4 }* o' _# S  d. r0 x: i4 B新冠病毒變異后會更有“侵略性”？
. V9 N% H0 U4 i1 J7 x9 _) W6 h      論文作者補(bǔ)充聲明
9 B' K; l- W  P6 k7 s0 Q      陸劍等人研究了103個(gè)新冠病毒基因組發(fā)現(xiàn)，由于一個(gè)蛋白編碼位點(diǎn)的突變產(chǎn)生的差別，將新冠病毒毒株清晰分為兩組，被媒體稱為“亞型”。但有學(xué)者認(rèn)為，目前變異造成的差異不足以叫“亞型”，叫“分枝”（clade）比較合適。
      陸劍等人在重新提交的論文中，做了補(bǔ)充性聲明：
/ w7 b9 H7 f& _" N2 W  d% t# |( f      “在我們最近發(fā)表的文章，分析結(jié)果顯示，103個(gè)SARS-CoV-2病毒基因組存在兩種不同的譜系；分別稱之為‘L’和‘S’譜系。我們用來定義L和S譜系的氨基酸位點(diǎn)位于ORF8（開放閱讀框8。所謂開放閱讀框是指是DNA序列中具有編碼蛋白質(zhì)潛能的序列——記者注。）基因，這個(gè)基因還沒有發(fā)現(xiàn)具有任何已知的重要功能。基于‘L’譜系的頻率高于譜系S的發(fā)現(xiàn)，我們將L譜系描述為‘a(chǎn)ggressive’（具有侵略性）。我們現(xiàn)在認(rèn)識到，在本研究闡述的內(nèi)容中，‘侵略性’一詞會具有誤導(dǎo)性，應(yīng)該用更精確的術(shù)語‘更高的頻率’代替。簡而言之，盡管我們已經(jīng)發(fā)現(xiàn)這兩個(gè)譜系自然并存，但我們沒有提供任何證據(jù)支持關(guān)于SARS-CoV-2毒力或致病性的任何流行病學(xué)結(jié)論。因此，我們將在本文的印刷版本中進(jìn)行更正，以避免產(chǎn)生誤導(dǎo)?！?/font>
      在此前的報(bào)道中，媒體多用“亞型”這一表述。從聲明中可以看出，陸劍等人并不認(rèn)為突變產(chǎn)生了兩個(gè)“亞型”，而是用了“譜系”這一說法。
" W5 H  G0 q) {+ z      在論文剛公布的時(shí)候，媒體報(bào)道，“新變異出的L亞型更具傳染性”。對于這個(gè)觀點(diǎn)，陸劍等人做了重點(diǎn)更正，認(rèn)為此前用詞具有誤導(dǎo)性。
" D* b0 ~, n; b5 R" Z- H% x6 D8 Y      新冠病毒有多愛突變？
      已記錄481個(gè)變異
) X5 B4 t1 `+ K" d2 ]3 b6 n6 n      湖北省疾控中心衛(wèi)生檢測檢驗(yàn)研究所流感參比實(shí)驗(yàn)室主任劉琳琳等人在另一篇論文中提到，根據(jù)數(shù)據(jù)庫，新冠病毒的基因組大小從29,409bp到29,911bp不等，包括183個(gè)可變位點(diǎn)。
      記者發(fā)現(xiàn)，截至3月12日5時(shí)，國家生物信息中心2019新型冠狀病毒信息庫中已經(jīng)記錄了482個(gè)新冠病毒的變異。
新型冠狀病毒信息庫記錄4種變異類型，分別是SNP（單核苷酸多性態(tài)，就是單個(gè)堿基的改變）、deletation（缺失，一個(gè)或者多個(gè)核苷酸丟失）、insertion（插入，多了一個(gè)或多個(gè)核苷酸）、indel（插入或缺失）。其中SNP最多，有461條。其次是deletation，有20條。因?yàn)樾鹿诓《臼荝NA單鏈，這種病毒SNP變異特別常見。
# J, g$ ~4 G" M: A* e

  \4 E0 S: X9 |! E  c$ l  P& P      新冠病毒各基因的變異情況。
      陸劍等人用來定義L和S譜系的氨基酸位點(diǎn)位于ORF8（開放閱讀框8）。從新型冠狀病毒信息庫的數(shù)據(jù)來看，ORF8的變異不是最多的。并且，陸劍等人表示這個(gè)基因還沒有發(fā)現(xiàn)具有任何已知的重要功能。
      記者了解到，S基因具有重要功能，這個(gè)基因用來生產(chǎn)S蛋白（spikeprotein，翻譯為棘突蛋白、長釘?shù)鞍?、穗蛋白），棘突蛋白會和病毒宿主?xì)胞受體ACE2結(jié)合，使得病毒能夠侵入到細(xì)胞內(nèi)。
6 y* B$ y' l6 s2 ?) G      在S基因中，有一處重要突變一直是各國科學(xué)家研究的重點(diǎn)。南開大學(xué)高山等人發(fā)現(xiàn)，新冠病毒的受體結(jié)合基序中多了一個(gè)Furin蛋白酶切位點(diǎn)，這在蝙蝠冠狀病毒和穿山甲冠狀病毒中都沒有發(fā)現(xiàn)。不同于SARS等其他大部分Beta冠狀病毒，這個(gè)變異使得新冠病毒可能采用了HIV等其他病毒的類似的包裝機(jī)制，它的棘突蛋白獲得了更高的侵染細(xì)胞的效率，更具感染性。
3 Z' h9 L8 q1 A; z5 y      兩個(gè)譜系有何差別？
/ Q; a0 h5 d) @! g      L型由S型衍化而來
      陸劍等人用完全連鎖的SNP區(qū)分了新冠病毒的兩種譜系，在103株新冠病毒病毒株中，有101株表現(xiàn)出兩個(gè)SNP之間的完全連鎖：72株表現(xiàn)出“CT”單倍型（定義為“L”型，因?yàn)槭橇涟彼岬拿艽a子），29株在這兩個(gè)位點(diǎn)表現(xiàn)出“TC”單倍型（定義為“S”型，因?yàn)槭墙z氨酸的密碼子）。由此可見，L為主要類型（占70％），S為次要類型（占30％）。
      他們將新冠病毒和其他高度相關(guān)的病毒基因組比對，發(fā)現(xiàn)S譜系的核苷酸與最密切相關(guān)的病毒的同源位點(diǎn)相同，L型比S型更多，但S型實(shí)際上更“老”。不過，盡管L型是從古代S型新進(jìn)化而來的，但它在人類群體中的傳播或復(fù)制速度更快，從而導(dǎo)致其積累的突變比S型更多。
      從武漢分離出的27種病毒中，有26種（96.3％）為L型，只有1種（3.7％）為S型。但是，在武漢以外地區(qū)分離出的其他73種病毒中，L型為45種（61.6％），S型為28種（38.4％）。這種比較表明，L型在武漢比其他地方普遍得多。

: U2 D8 j+ G- Z" Z4 L0 z
! h$ N8 H2 @( y9 h5 f/ H2 `      新冠病毒分化的兩個(gè)譜系。
1 E8 u. n. v9 @2 L- A; m      劉琳琳等人做的研究，也顯示出了類似的進(jìn)化和分布關(guān)系。
      他們通過分析146個(gè)基因組，發(fā)現(xiàn)了183個(gè)替換位點(diǎn)，將新冠病毒分為80個(gè)單倍型。其中，與其祖先單倍型相比，一個(gè)被稱為H3的單倍型及其衍生的單倍型數(shù)量很大，并在湖北（尤其武漢）廣泛分布。除中國外，其他15個(gè)國家提供的65個(gè)病毒樣本被分配給43個(gè)單倍型。其中，H3衍生出來的單倍型有27種，它的兩個(gè)祖先單倍型H14衍生的單倍型11種，H15衍生的單倍型3種。
& W3 E+ k3 e+ ~* u( q2 q

      各單倍型之間的衍化關(guān)系。
  x; ]1 _  U& {3 o. O+ k8 R      值得一提的是，陸劍等人發(fā)現(xiàn)，有兩個(gè)樣本包含了S和L兩種類型，這兩個(gè)樣本都來自國外——美國和澳大利亞。他們猜測，患者可能感染了至少兩種不同的新冠病毒。
* J7 W4 \! L1 L0 s+ p6 r  S      直接來源到底是哪種動(dòng)物？
6 F  J: W, l) h1 S# s1 l      穿山甲不太可能
5 x+ Y1 m- n( I8 A7 Q      陸劍等人除了研究新冠病毒的進(jìn)化外，還研究了新冠病毒的起源，著重分析了是來自蝙蝠還是穿山甲。
      全基因組系統(tǒng)發(fā)生樹表明新冠病毒最接近蝙蝠冠狀病毒（RaTG13），其次是穿山甲相關(guān)冠狀病毒，然后是人的SARS病毒。劉琳琳等人也得出了相同的結(jié)論，基于系統(tǒng)發(fā)生樹，蝙蝠冠狀病毒（bat-RaTG13-CoV）與新冠病毒和兩個(gè)從廣東采集的穿山甲冠狀病毒的姊妹病毒的遺傳距離相對較大，穿山甲可能不是新冠病毒中間宿主。

& e( ~3 `1 t4 c  I
      蝙蝠、穿山甲和人類新冠病毒的遺傳距離。
3 Q7 x) e! Z5 K. F9 s1 b      不過，陸劍等人發(fā)現(xiàn)，盡管新冠病毒與蝙蝠SARS相關(guān)冠狀病毒（SARSr-CoV-RaTG13）之間的基因組核苷酸差異僅4％，但中性位點(diǎn)的差異為17％，表明兩種病毒之間的差異比以前估計(jì)的要大要多。作者做了一個(gè)對比，這個(gè)差異是人與黑猩猩之間差距的14倍，是人與恒河猴之間差距的2倍。
      此外，有一些研究表明穿山甲可能為新冠病毒提供了部分棘突基因（上文說的S基因）。棘突蛋白受體結(jié)合域中的6個(gè)關(guān)鍵氨基酸殘基，和SARS病毒相比，新冠病毒有5個(gè)殘基都不同，并且新冠病毒與受體ACE2的結(jié)合親和力更高。蝙蝠冠狀病毒（RaTG13）也是6個(gè)中只有1個(gè)相同，但有趣的是，穿山甲冠狀病毒和新冠病毒的6個(gè)關(guān)鍵氨基酸殘基是相同的。
      有人提出，新冠病毒的棘突蛋白受體結(jié)合域可能是由穿山甲中最近的重組事件引起的。陸劍等人認(rèn)為，盡管棘突蛋白產(chǎn)生過重組，但新冠病毒和穿山甲冠狀病毒中相同的功能位點(diǎn)也可能是巧合的趨同進(jìn)化的結(jié)果，香港大學(xué)管軼等人在此前的論文中表達(dá)了相同的觀點(diǎn)。
      【主要參考論文】
Genome-widedata inferring the evolution and population demography of the novel pneumonia coronavirus (SARS-CoV-2). bioRxiv2020.03.04.976662
On the origin and continuing evolution of SARS-CoV-2. National Science Review, nwaa036
  ]6 c- ^- D+ d; p- C& N      新京報(bào)記者 李玉坤


6 H0 b3 I+ H6 Y& T  Y# j

2 N0 U* }( E6 l. Y- _4 T, Q
, w- F  n1 q$ Y3 R8 K) w8 ^


& ?* M) b& R; K- j! n

6 y8 H) v3 u8 B' `. z& r


9 h' p1 s( ?1 n0 E9 w- v& |3 ~( N3 C

$ O, Y; V% ?# Y/ n4 M, c* ~

桃李無言