2019新型冠狀病毒基因組的生物信息學(xué)分析

陳嘉源，施勁松，丘棟安，劉 暢，李 鑫，趙 強(qiáng)，阮吉壽，高 山*

(1.南開大學(xué) 生命科學(xué)學(xué)院, 天津300071; 2.東部戰(zhàn)區(qū)總醫(yī)院, 南京 210016;3.英國諾丁漢特倫特大學(xué) 生物科學(xué)系, 諾丁漢 NG11 8NS;4.南開大學(xué) 醫(yī)學(xué)院, 天津300071;5.南開大學(xué) 數(shù)學(xué)科學(xué)學(xué)院, 天津 300071)

冠狀病毒屬(Coronavirus)的病毒是具有包膜的正鏈RNA病毒。由于2003年嚴(yán)重急性呼吸綜合征(Severe Acute Respiratory Syndrome, SARS)和2012年中東呼吸綜合征(Middle East Respiratory Syndrome, MERS)的爆發(fā)，冠狀病毒逐漸成為病毒學(xué)領(lǐng)域的一個(gè)研究熱點(diǎn)。2019年12月，中國武漢報(bào)道了冠狀病毒引起的肺炎，其臨床癥狀與2003年爆發(fā)的SARS不同，因此推斷該病毒可能是冠狀病毒的一個(gè)新變種。當(dāng)前，冠狀病毒溯源相對(duì)比較清晰：2003年，Enserink等在國際上首次提出了SARS冠狀病毒(SARS Coronavirus, SARS-CoV)來自蝙蝠，而后通過果子貍傳播給人類[1]；再后來的病毒溯源研究都是根據(jù)蝙蝠-果子貍-人途徑回溯，并重點(diǎn)研究各種蝙蝠體內(nèi)的SARS樣(SARS-like)冠狀病毒與已知SARS冠狀病毒的進(jìn)化關(guān)系；特別是2013年，SARS冠狀病毒(GenBank: AY274119.3)被溯源到中華菊頭蝠(Rhinolophussinicus)[2]。盡管已有足夠的冠狀病毒基因組數(shù)據(jù)公開，有關(guān)冠狀病毒跨物種傳播和宿主適應(yīng)性方面的研究卻一直難以突破，特別是無法解釋冠狀病毒變異快、宿主多且具有較強(qiáng)的宿主適應(yīng)性等特點(diǎn)。其中一個(gè)最主要原因就是冠狀病毒基因組較大，多個(gè)基因功能未知。例如，SARS冠狀病毒基因組(GenBank: AY274119.3)的全長是29 751 bp，比HIV-1病毒全長三倍還要多。另外，冠狀病毒基因組注釋中基因命名不統(tǒng)一、且大部分注釋的蛋白質(zhì)序列都來自預(yù)測(cè)，無實(shí)驗(yàn)驗(yàn)證，難以保證正確性。高通量測(cè)序的迅速發(fā)展，使大量冠狀病毒基因組序列得以公開。然而，當(dāng)前大部分的冠狀病毒基因組研究都是簡單使用全基因組或某個(gè)病毒結(jié)構(gòu)基因的序列，后一種做法主觀性較大。簡單使用全基因組序列進(jìn)行進(jìn)化分析和重要功能或致病位點(diǎn)分析，不僅工作量巨大，而且大量與研究目的無關(guān)的基因組區(qū)域會(huì)對(duì)數(shù)據(jù)分析產(chǎn)生干擾。因此，不同于簡單使用全基因組序列的其它研究，我們于2018年在國際上首次提出分子功能與進(jìn)化分析相結(jié)合的研究思想，并應(yīng)用于冠狀病毒基因組的研究。分子功能與進(jìn)化分析相結(jié)合，即僅使用與特定功能相關(guān)的基因組特征對(duì)一定變異范圍內(nèi)的序列進(jìn)行進(jìn)化分析[3]。這樣不僅可以簡化進(jìn)化分析、提高信噪比，使分析結(jié)果更可靠，而且可以用于指導(dǎo)功能實(shí)驗(yàn)。2018年, 南開大學(xué)高山等在國際上首次報(bào)道了SARS冠狀病毒中存在互補(bǔ)回文小RNA(complemented palindromic Small RNAs, cpsRNAs)的現(xiàn)象，并確定了Beta冠狀病毒B亞群(BB冠狀病毒)與功能相關(guān)的兩個(gè)重要的基因組特征(Genomic feature)[3]，即一段22 bp的互補(bǔ)回文序列(命名為Nankai complemented palindrome)與其所在的編碼區(qū)(命名為Nankai CDS)。在上述研究中，對(duì)Nankai complemented palindrome進(jìn)行突變分析的結(jié)果與使用Nankai CDS進(jìn)行進(jìn)化分析的結(jié)果一致支持SARS冠狀病毒的跨物種傳播途徑是蝙蝠-果子貍-人。在分子功能與進(jìn)化分析相結(jié)合的研究思想的指導(dǎo)下，本研究使用Nankai complemented palindrome和Nankai CDS對(duì)新發(fā)布的2019新型冠狀病毒基因組(GenBank: MN908947)進(jìn)行分析以期準(zhǔn)確溯源，并對(duì)BB冠狀病毒跨物種傳播和宿主適應(yīng)性進(jìn)行初步研究。

1 數(shù)據(jù)與方法

在前期研究中，我們通過互補(bǔ)回文小RNA的發(fā)現(xiàn)從SARS冠狀病毒MA-15株(GenBank: DQ497008)中確定了BB冠狀病毒中與功能相關(guān)的兩個(gè)重要的基因組特征，即一段22 bp的互補(bǔ)回文序列與其所在的編碼區(qū)。這個(gè)編碼區(qū)(DQ497008: 25689-26153)是BB冠狀病毒基因組中一段高度保守的序列，在已知SARS冠狀病毒基因組(如DQ497008)中只對(duì)應(yīng)一個(gè)讀碼框ORF 3b；但在其它BB冠狀病毒基因組中普遍對(duì)應(yīng)多個(gè)讀碼框，因此ORF 3b無法代表其它BB冠狀病毒基因組中的這段序列。本研究將22 bp的互補(bǔ)回文序列與其所在的編碼區(qū)(可能是一個(gè)或多個(gè)讀碼框)分別命名為Nankai complemented palindrome和Nankai CDS。DQ497008和2003年爆發(fā)的SARS冠狀病毒基因組(GenBank: AY274119.3)含有相同的Nankai complemented palindrome和Nankai CDS，因此在本研究中唯一代表SARS冠狀病毒基因組。在前期研究中[3]，我們的研究對(duì)象僅限于BB冠狀病毒(本研究亦是如此)，并進(jìn)行如下操作[3]：提取NCBI GenBank數(shù)據(jù)庫中全部BB冠狀病毒的完整基因組(Complete genome)序列(以下簡稱病毒序列或序列)；去除與Nankai CDS相似度在0.999以上的序列后，留下11條病毒序列(GenBank: JX993987、JX993988、GQ153539、GQ153540、GQ153542、DQ071615、DQ412042、DQ412043、AY515512、AY572034和DQ497008)用于進(jìn)一步研究。本研究增加一條新發(fā)布的2019新型冠狀病毒基因組的序列(GenBank: MN908947)和一條來自蝙蝠的序列(GenBank: MG772934)，其Nankai CDS與2019新型冠狀病毒Nankai CDS相似度最高(92.52%)。另外，由于本論文接收時(shí)2019新型冠狀病毒基因組序列(GenBank: MN908947)尚未公開，實(shí)際使用的是來自互聯(lián)網(wǎng)公開的序列(見致謝)。在本研究中，13條序列根據(jù)其宿主分為五組用于進(jìn)一步研究，這五組命名為SARS(DQ497008)、果子貍Civet(AY515512和AY572034)、2019新型冠狀病毒2019-nCoV(MN908947)、蝙蝠群體bat2(MG772934)和蝙蝠群體bat1(MG772934之外8條來自蝙蝠的序列)。在本研究中，進(jìn)化分析使用軟件PHYLIP v3.69，統(tǒng)計(jì)與作圖使用軟件R v2.15.3[4]。其它數(shù)據(jù)處理包括：病毒蛋白質(zhì)預(yù)測(cè)和核苷酸三連子分析只考慮一個(gè)起始密碼子ATG、三個(gè)終止密碼子(TAA、TAG和TGA)和正向三個(gè)讀碼框；蛋白質(zhì)預(yù)測(cè)不考慮小于20個(gè)氨基酸的蛋白質(zhì)；最后，考慮四個(gè)起始密碼子(ATG、GTG、CTG和TTG)和三個(gè)終止密碼子重復(fù)蛋白質(zhì)預(yù)測(cè)和核苷酸三連子分析以檢驗(yàn)結(jié)果的可靠性(見結(jié)果)。

2 結(jié) 果

在經(jīng)典定義中，DNA回文(DNA palindrome)和DNA互補(bǔ)回文(DNA complemented palindrome)序列都叫做DNA回文序列，而且大多數(shù)情況下，DNA回文序列的經(jīng)典定義與我們新定義中的互補(bǔ)回文序列相同(例如ATCGCGAT)。根據(jù)我們的新定義[3]，回文序列(DNA或RNA)要求序列(一條鏈)從5’到3’方向與從3’到5’方向的讀取結(jié)果相同(例如ATCGGCTA)?；匚暮突パa(bǔ)回文序列應(yīng)該加以區(qū)別，主要基于三個(gè)原因：(1)我們的前期研究發(fā)現(xiàn)了回文和互補(bǔ)回文序列可能具有不同的生物學(xué)意義；(2)在動(dòng)物病毒基因組中，回文和互補(bǔ)回文序列具有不同的分布特征；(3)經(jīng)典定義是從DNA水平定義的，我們的定義考慮了單鏈DNA和RNA(特別針對(duì)病毒)情況。2004年，Chew等在國際上首次報(bào)道了SARS冠狀病毒基因組中較短互補(bǔ)回文序列(新定義)的分布特征[5]，其主要發(fā)現(xiàn)包括兩點(diǎn)：(1)在所有分析的冠狀病毒基因組中，長度為4 bp的互補(bǔ)回文序列出現(xiàn)頻率顯著較低；(2)長度為6 bp的互補(bǔ)回文序列僅在SARS冠狀病毒(而不是所有冠狀病毒)基因組中出現(xiàn)頻率顯著較低。因此，該研究的結(jié)論是SARS冠狀病毒基因組含有較少的6 bp互補(bǔ)回文序列有可能利于該病毒有效規(guī)避宿主細(xì)胞內(nèi)的某些防御機(jī)制而有利于病毒存活；該研究同時(shí)還發(fā)現(xiàn)了兩個(gè)較長(14 bp以上)的互補(bǔ)回文序列分別是TCTTTAACAAGCTTGTTAAAGA和TAAAATTAATTTTA。但是，僅僅根據(jù)概率模型計(jì)算得到的此類回文或互補(bǔ)回文序列數(shù)量巨大，絕大部分無法與分子功能關(guān)聯(lián)，因此沒有進(jìn)一步研究的價(jià)值。直到2018年，我們偶然發(fā)現(xiàn)了互補(bǔ)回文小RNA恰好來自這個(gè)22 bp的互補(bǔ)回文序列(命名為Nankai complemented palindrome)，特別是發(fā)現(xiàn)了根據(jù)這個(gè)序列預(yù)測(cè)的RNA二級(jí)結(jié)構(gòu)在BB冠狀病毒基因組中高度保守(見圖1A)。RNA水平的證據(jù)是將這個(gè)序列與分子功能建立關(guān)聯(lián)的關(guān)鍵[3]。在前期研究中，我們對(duì)11個(gè)BB冠狀病毒(見數(shù)據(jù)與方法)的Nankai complemented palindrome進(jìn)行突變分析的結(jié)果與使用其所在的編碼區(qū)(命名為Nankai CDS)進(jìn)行進(jìn)化分析的結(jié)果一致支持SARS冠狀病毒的跨物種傳播途徑是蝙蝠-果子貍-人(見圖1B紅色方框內(nèi))[3]。

圖1 2019新型冠狀病毒的起源與可變翻譯Fig.1 Origin and alternative translation of the 2019 novel coronavirus genome

注：A：Nankai CDS是BB冠狀病毒基因組中一段高度保守的序列，包括一段22 bp的互補(bǔ)回文序列(用*表示)。Nankai CDS還包括一個(gè)8或11 bp的可變區(qū)，分別屬于兩種不同長度(是465或468 bp，用#表示)的Nankai CDS；B：進(jìn)化樹構(gòu)建使用13條去除可變區(qū)的Nankai CDS。使用多種方法進(jìn)行進(jìn)化分析的結(jié)果一致，這里只展示鄰接(Neighbor Joining，NJ)法的結(jié)果：2019新型冠狀病毒源自蝙蝠某個(gè)群體(藍(lán)色方框內(nèi))，但與SARS冠狀病毒差異巨大(紅色方框內(nèi))；用P1到P17(見圖1C)可以對(duì)13個(gè)BB冠狀病毒進(jìn)行基因分型(小括號(hào)內(nèi))，其結(jié)果與進(jìn)化分析的結(jié)果一致。C：可變翻譯導(dǎo)致從Nankai CDS可以預(yù)測(cè)出至少17個(gè)蛋白質(zhì)，分別命名為P1到P17。這里展示的是使用一個(gè)起始密碼子ATG和三個(gè)終止密碼子進(jìn)行蛋白質(zhì)預(yù)測(cè)的結(jié)果，使用四個(gè)起始密碼子(ATG、GTG、CTG和TTG)和三個(gè)終止密碼子會(huì)預(yù)測(cè)更多的蛋白質(zhì)，但不改變可變翻譯相關(guān)結(jié)論。

Notes:A: The 22-bp complemented palindrome (named Nankai complemented palindrome) and the CDS (named Nankai CDS) have a high degree of evolutionary conservation in genomes of betacoronavirus subgroup B(BB coronavirus). An 8- or 11-bp long variable region in Nankai CDS resides in two types of Nankai CDS. The first type has a length of 465 bp, while the second type has a length of 468 bp.B: Consistent results were obtained by different phylogenetic analysis methods. Here, we only show the result using the NJ (Neighbor Joining) method. The 2019 novel coronavirus (in blue box) with large differences from the SARS coronavirus (in red box), may originate from BB coronaviruses in Chinese horseshoe bats. Phylogenetic analysis and genotyping (in parentheses) of 13 viruses using 17 putative proteins (Fig.1C) achieved consistent results. C: The putative proteins (named as P1 to P17) were predicted from the Nankai CDS using ATG as the start codon, and TAA, TAG and TGA as stop codons. The prediction using ATG, GTG, CTG and TTG as start codons, and TAA, TAG and TGA as stop codons did not change all conclusions in the present study.

在本研究中，所有13條Nankai CDS(見數(shù)據(jù)與方法)都含有一個(gè)完整的Nankai complemented palindrome和一個(gè)緊鄰其3'端的長度為8或11 bp的可變區(qū)(Variable region)。使用MN908947(2019新型冠狀病毒)中Nankai complemented palindrome預(yù)測(cè)的RNA二級(jí)結(jié)構(gòu)仍然能形成穩(wěn)定的莖環(huán)(Stem-loop)結(jié)構(gòu)，再次驗(yàn)證了Nankai complemented palindrome在BB冠狀病毒基因組中的保守性[3]。Nankai CDS中的可變區(qū)導(dǎo)致了長度變異(見圖1A)：MG772934(來自蝙蝠)和MN908947(2019新型冠狀病毒)的Nankai CDS全長為468 bp(含11 bp可變區(qū))；其余11條序列的Nankai CDS全長為465 bp(含8 bp可變區(qū))。由于無法確定每條Nankai CDS的可變區(qū)的變異方式，進(jìn)行進(jìn)化分析前，13條Nankai CDS去除可變區(qū)后即可對(duì)齊，長度統(tǒng)一到457 bp。使用多種方法(UPGMA、鄰接法、最大簡約法、最大似然法和貝葉斯法)對(duì)13條Nankai CDS進(jìn)行進(jìn)化分析的結(jié)果(見圖1B)一致支持：(1)MN908947(2019新型冠狀病毒)與MG772934(來自蝙蝠)在進(jìn)化樹的一個(gè)分支內(nèi)(見圖1B藍(lán)色方框內(nèi))，說明2019新型冠狀病毒源自蝙蝠；(2)2019新型冠狀病毒(見圖1B藍(lán)色方框內(nèi))與SARS冠狀病毒(見圖1B紅色方框內(nèi))在不同分支，而且距離很大，這一結(jié)果與兩者臨床癥狀差異一致；(3)SARS冠狀病毒的跨物種傳播途徑“蝙蝠-果子貍-人”(見圖1B紅色方框內(nèi))與2018年分析結(jié)果完全相同[3]，證明了我們方法的可靠性。

Nankai CDS是BB冠狀病毒基因組中一段高度保守的序列(見數(shù)據(jù)與方法)，靠近其3'端的99個(gè)堿基(編碼33個(gè)氨基酸)極端保守，不僅沒有起始或終止密碼子突變，也極少有其它突變，具有重要的研究和應(yīng)用(如引物設(shè)計(jì))價(jià)值。Nankai CDS在來自SARS和果子貍的beta冠狀病毒基因組中是一個(gè)完整的讀碼框，編碼蛋白P16(見圖1C)，但在來自其它宿主的冠狀病毒基因組中(不包括3'端的99個(gè)堿基)出現(xiàn)了大量的起始或終止密碼子突變(見表1)，并導(dǎo)致Nankai CDS出現(xiàn)讀碼框的多樣性，即可變翻譯。由于去冗余后的完整基因組數(shù)據(jù)較少，需要從兩個(gè)方面排除這些突變是來自測(cè)序錯(cuò)誤等導(dǎo)致的假陽性。一方面，用13條Nankai CDS比對(duì)到NCBI NT數(shù)據(jù)庫，看是否有其它冠狀病毒序列支持這些突變。比對(duì)結(jié)果顯示，大部分序列上的突變都能得到NT數(shù)據(jù)庫中其它數(shù)據(jù)支持，只有三個(gè)例外：JX993988在Nankai CDS上的三個(gè)突變(C80A、T80A和T254A)，MG772934(來自蝙蝠)和MN908947(2019新型冠狀病毒)在Nankai CDS上共有的四個(gè)突變(T182A、C182A、T237A和T236G)沒有其它數(shù)據(jù)支持。JX993988代表的病毒類型比較特殊，與其它來自蝙蝠的病毒差異較大(見圖1B)，這里不做進(jìn)一步研究。另外一方面，統(tǒng)計(jì)13條Nankai CDS中核苷酸三連子出現(xiàn)次數(shù)，次數(shù)最高的前14位依次是：ATG(123次)、ACT(106次)、AAG(98次)、TTG(83次)、CAG(77次)、CAA(71次)、GTG(68次)、AAA(66次)、CTA(65次)、ACA(60次)、TTA(59次)、ATT(58次)、CGA(57次)和CTG(55次)。核苷酸三連子統(tǒng)計(jì)結(jié)果顯示：前14位三連子出現(xiàn)次數(shù)之和超過三連子總數(shù)一半以上(1 046次/2 016次)；6個(gè)三連子與起始密碼子ATG只有0或1個(gè)核苷酸的差異；8個(gè)三連子與終止密碼子(TAA、TAG和TGA)只有1個(gè)核苷酸的差異。統(tǒng)計(jì)結(jié)果說明了通過1個(gè)核苷酸突變(點(diǎn)突變)導(dǎo)致起始或終止密碼子的丟失和獲得的概率非常高。根據(jù)以上兩個(gè)方面的分析結(jié)果，基本可以排除Nankai CDS中的大量起始或終止密碼子突變是假陽性(見表1)。

表1 Nankai CDS內(nèi)起始或終止密碼子突變Table 1 Mutations in start or stop codons in Nankai CDS

注：可變翻譯導(dǎo)致從Nankai CDS可以預(yù)測(cè)出至少17個(gè)蛋白質(zhì)，分別命名為P1到P17(見圖1C)。所有突變均來自17個(gè)預(yù)測(cè)的蛋白質(zhì)，PP表示預(yù)測(cè)的蛋白質(zhì)(Putative Protein)的名稱；Start一列記錄的是起始密碼子的第一個(gè)堿基的位置；Stop一列記錄的是終止密碼子的最后一個(gè)堿基位置；Allele一列記錄所在位置的密碼子的全部類型；Mutation一列只記錄導(dǎo)致獲得起始或終止密碼子的突變，突變格式為“突變前堿基-突變位置-突變后堿基”；AA表示PP的長度，單位是氨基酸(Amino Acid)。這里展示的是只使用一個(gè)起始密碼子ATG和三個(gè)終止密碼子進(jìn)行蛋白質(zhì)預(yù)測(cè)的結(jié)果；使用四個(gè)起始密碼子(ATG、GTG、CTG和TTG)和三個(gè)終止密碼子會(huì)預(yù)測(cè)更多的蛋白質(zhì)，但不改變可變翻譯相關(guān)結(jié)論。

Notes:The putative proteins (named as P1 to P17) were predicted from the Nankai CDS (Fig.1C). This table only includes the mutations in start or stop codons in 17 putative proteins. “PP” is Putative Protein; “Start” records the positions of the first nucleotides in start codons; “Stop” records the positions of the last nucleotides in stop codons; “Allele” records all types of codons on the positions; “Mutation” only records mutations which cause the acquisition of start or stop codons in the format “nucleotide before mutaion-position-nucleotide after mutaion”; “AA” is Amino Acid. The putative proteins were predicted from the Nankai CDS using ATG as the start codon, and TAA, TAG and TGA as stop codons. The prediction using ATG, GTG, CTG and TTG as start codons, and TAA, TAG and TGA as stop codons did not change all conclusions in the present study.

可變翻譯導(dǎo)致從Nankai CDS可以預(yù)測(cè)出至少17個(gè)蛋白質(zhì)(Putative Protein)，分別命名為P1到P17(見圖1C)。用P1到P17可以對(duì)13個(gè)病毒進(jìn)行基因分型(基于可變翻譯的基因分型)，例如，MN908947(2019新型冠狀病毒)可以表示為P10P11P12P13P17。我們發(fā)現(xiàn)基因分型相同或相近的病毒都在進(jìn)化樹上的同一個(gè)分支內(nèi)(見圖1B)。例如，bat1分支的病毒的基因分型都由三個(gè)蛋白質(zhì)組成；SARS和果子貍分支都是P16；MN908947(2019新型冠狀病毒)和MG772934(來自蝙蝠)都有P10。根據(jù)以上事實(shí)，可以得到以下推論：(1)BB冠狀病毒(不限于Nankai CDS內(nèi))存在大量的可變翻譯；(2)BB冠狀病毒變異快、多樣性高，僅僅13個(gè)病毒即可得到多達(dá)9種基因分型(見圖1BC)；(3)在SARS冠狀病毒從蝙蝠到人傳播過程中，Nankai CDS內(nèi)多個(gè)點(diǎn)突變積累才形成P16；(4)在2019新型冠狀病毒從蝙蝠到人傳播過程中，P15只需一個(gè)點(diǎn)突變(C355T)即形成P12和P13，P14只需兩個(gè)點(diǎn)突變(T170A和T191C)即形成P11和P17；(5)可變翻譯可用于(但不限于)BB冠狀病毒的基因分型，用P1到P17進(jìn)行基因分型的結(jié)果與用Nankai CDS進(jìn)行進(jìn)化分析的結(jié)果一致。(6)BB冠狀病毒可能通過可變翻譯以適應(yīng)不同宿主。

3 結(jié) 論

本研究證明了我們提出的分子功能與進(jìn)化分析相結(jié)合的研究思想的可行性與可靠性，對(duì)生物大分子的功能或進(jìn)化研究具有一定的指導(dǎo)意義，并對(duì)病毒或細(xì)菌等病原體的研究提供了新的路線和方向。因此，本研究除了解決2019新型冠狀病毒溯源等具體問題，更主要貢獻(xiàn)是提供了研究思想和方法。

1)使用Nankai CDS溯源分析的結(jié)果支持2019新型冠狀病毒源自蝙蝠，但與SARS冠狀病毒差異巨大；

2)從BB冠狀病毒可變翻譯中獲取的信息可應(yīng)用于其快速檢測(cè)、基因分型、疫苗開發(fā)以及藥物設(shè)計(jì)；

3)BB冠狀病毒可能通過可變翻譯以適應(yīng)不同宿主；

4)對(duì)BB冠狀病毒可變翻譯的研究將有助于研究該病毒的變異、宿主適應(yīng)性、感染以及致病機(jī)制等問題。

4 討 論

病毒或細(xì)菌等原核生物的可變翻譯現(xiàn)象很早就已經(jīng)發(fā)現(xiàn)，但是其生物學(xué)意義尚未清楚，特別是缺乏基于較大數(shù)據(jù)的實(shí)證研究。我們?cè)谇捌谘芯恐参锊《綶6]、昆蟲病毒[7]、人病毒[8]、哺乳動(dòng)物病毒(特別是非洲豬瘟病毒[9])和布魯氏菌[10]的過程中，都沒有明顯觀察到可變翻譯現(xiàn)象，特別是BB冠狀病毒Nankai CDS中這樣大量而集中存在的可變翻譯(僅僅13個(gè)病毒即可得到多達(dá)9種基因分型)極其罕見。此次實(shí)證研究結(jié)果達(dá)到了基于可變翻譯的基因分型與進(jìn)化分析的一致性，是一次較大突破，主要得益于以下幾點(diǎn)：(1)分子功能與進(jìn)化分析相結(jié)合的研究思想得以應(yīng)用；(2)有蝙蝠-果子貍-人這一跨物種傳播途徑作為方向性指導(dǎo)；(3)無意中發(fā)現(xiàn)Nankai complemented palindrome和Nankai CDS；(4)來自多種宿主(蝙蝠、果子貍和人)的病毒數(shù)據(jù)，特別是陽性樣本(如來自SARS和武漢肺炎患者)數(shù)據(jù)的積累。經(jīng)典理論認(rèn)為病毒與宿主共用一套翻譯系統(tǒng)，然而病毒的蛋白質(zhì)翻譯也有其特殊性，導(dǎo)致本研究從Nankai CDS中預(yù)測(cè)的17種蛋白質(zhì)不能提供準(zhǔn)確的序列。蛋白質(zhì)序列的準(zhǔn)確度不影響本研究的所有結(jié)論和推論，而且為下一步蛋白質(zhì)、細(xì)胞和動(dòng)物水平的實(shí)驗(yàn)指出方向。

基于可變翻譯的基因分型相同的病毒都在進(jìn)化樹上的同一個(gè)分支(同種宿主)內(nèi)預(yù)示了beta冠狀病毒可能通過可變翻譯以適應(yīng)不同宿主(見圖1B)。但是，由于目前能獲得的病毒數(shù)據(jù)僅僅來自少量宿主(蝙蝠、果子貍和人)，這一推斷的驗(yàn)證還需要積累來自更多種類宿主的病毒的數(shù)據(jù)。我們推測(cè)，病毒通過可變翻譯以適應(yīng)不同宿主作為一種普遍機(jī)制，并不限于冠狀病毒或原核生物；真核生物中某些基因的可變剪接可能也有相似或相近的生物學(xué)功能。

BB冠狀病毒通過可變翻譯產(chǎn)生的多種蛋白質(zhì)的致病性研究是一個(gè)更為重要的課題。當(dāng)前公開的基因組數(shù)據(jù)雖然總量很大，但是，實(shí)際有效的陽性樣本(有實(shí)驗(yàn)或臨床記錄證明其致病性)種類不足和缺乏陰性(臨床實(shí)驗(yàn)確認(rèn)的非SARS病毒)樣本。致病性最明確的2019新型冠狀病毒和SARS冠狀病毒在Nankai CDS上各產(chǎn)生五個(gè)小蛋白和一個(gè)大蛋白(見圖1C)，并沒有顯示出蛋白質(zhì)某個(gè)屬性(如大小)與致病性的關(guān)系。在預(yù)測(cè)的17個(gè)蛋白質(zhì)中，只有P16(長度為155個(gè)氨基酸)和P9(長度為115個(gè)氨基酸)是大蛋白質(zhì)(見表1)；P16(來自SARS冠狀病毒)的致病性已知，而P9(來自蝙蝠的冠狀病毒)未知。因此，來自蝙蝠的其它冠狀病毒產(chǎn)生的蛋白(如P9)未必沒有致病性，或許也具有導(dǎo)致2019肺炎或SARS的能力。本研究預(yù)測(cè)的17種蛋白質(zhì)大小適中而且差異明顯，可以很容易通過蛋白質(zhì)組或外源表達(dá)等技術(shù)進(jìn)一步驗(yàn)證其致病性。

致謝：感謝新型冠狀病毒相關(guān)工作人員(省疾控中心等)的基礎(chǔ)性工作，特別復(fù)旦大學(xué)張永振等公開的2019新型冠狀病毒基因組數(shù)據(jù)，感謝南開大學(xué)生命科學(xué)學(xué)院動(dòng)物系卜文俊、張濤、黃大衛(wèi)、劉燕強(qiáng)和賀秉軍等各位老師對(duì)我們生物信息學(xué)研究的長期支持。