作者naieeliao (♪~~)
看板nCoV2019
標題[整理] 病毒演化系統發生樹
時間Sun Apr 19 18:37:27 2020
這裡只是來宣傳一下這兩週做的小project
https://github.com/PinguPinga/hCoV-2019-phylogenetic-analysis
從 GISAID 的資料庫抓病毒定序資料
Project1 => 篩出高品質的序列
Project2 => 分析並畫出樹狀圖
Git裡面有pdf檔可以參考
https://i.imgur.com/rELDNzQ.png
當初看 nextstrain.org 的資料
自己簡單玩了一下發現有出入 (C241T是假的，定序誤差)
所以就把 GISAID 原始的序列抓來分析了
https://i.imgur.com/WhUWCF9.png
概念圖
樹根的原點是使用大家統一定義的參考序列，採集地點在武漢，這並不重要
真要追溯到蝙蝠，說真的，沒完沒了
接下來要講故事了
https://imgur.com/kHlMZD5
一、鑽石公主號
node #450 (G11083T) EPI_ISL_408481 chongqing 2020-01-18
雖然在重慶第一次檢驗到，但實際上很多助攻來自第一艘爆炸的鑽石公主號
如果回過頭仔細分析不完整的序列，G11083T 貢獻全球 1/6 的上傳序列 (持續降低中)
還記得鑽石公主號乘客返回以色列之後被驗出來的事嗎？ 他有 G11083T
但這是不是日本該背的鍋呢？
1. 目前最盛行這個突變的地域是歐洲 (合併 G26144T)
G26144T 偏偏也有獨立的點突變
2. 日本目前很少報告有這個突變 (可是日本最近爆炸，卻很少上傳序列)
3. 早在鑽石公主號爆炸之前，重慶、香港、新加坡就有這個突變的個案了
結論是...很難講啊 我也不清楚
https://i.imgur.com/VPX47cI.png
二、惡名昭彰的Italy
node #136 EPI_ISL_412973 italy 2020-02-20 (C3037T C14408T A23403G)
其實在他前面有出現過一個
EPI_ISL_406862 germany 2020-01-28 (C3037T A23403G)
但是德國早期很快就遏止一間公司的群聚感染了
序列重複上傳筆數只有1，所以被程式篩掉了
但是義大利真的害人不淺阿，太多觀光客了，一下就傳到世界各處
這個分支占了 60% 以上的筆數，歐美遍地都是他的子孫
node #375 India 就是他的子孫之一
印度上傳的序列過半 (9/12) 都是他，在人口世界第二的國家，應該會爆炸吧 (?)
node #1159 Russia 也是 (11/22)
看來俄羅斯要起飛了
https://i.imgur.com/no8EZUk.png
三、Make America Great Again
事實上跟義大利同時發生的，是二月下旬到三月上旬的美國分支
美國現在最大的流派還是三月初暴漲的義大利系 (美國鎖國前的飛機超毒)
node #1 EPI_ISL_413456 usa 2020-02-20
C8782T C17747T A17858G C18060T T28144C
這個分支主要在北美跟澳洲出現，普遍度還是遠遠不及義大利的株系
https://i.imgur.com/ZMI0Qc5.png
四、日漸壯大的西班牙無敵艦隊
聽說最近西班牙死比較多人？ 其實西班牙的資料只侷限在兩個支系
很抱歉，這株系比較難從大量數據中找到源頭
node #409 EPI_ISL_418248 spain 2020-03-01
C8782T T9477A C14805T G25979T T28144C C28657T C28863T
智利就中獎了... 南美洲會不會都是他呢？
澳洲有一個小分支，接下來討論
五、世界的大贏家，澳洲
當你以為紐澳離亞洲、歐洲、美洲都隔著汪洋大海時
澳洲真的各個病毒株系都有！一二月亞洲的株系、三四月的歐美株系
雖然說澳洲表面上好像趨緩，但現在他可是集大成者！
實際上我們還不清楚究竟對不同株系會不會有共同的免疫力
聽說最近紐澳要形成共同的 travel bubble (往來的班機復航)
拭 目 以 待
六、蓋牌專家，中國
自從三月後上傳的有效資料數.... 手指數得出來
聽說最近有班機，我們可以來驗驗看有沒有新的演化支出現
結論：地理位置對於病毒的演化有很大的關係
冰島就有很多專屬的病毒株系，大英帝國也有少部分
以上是我個人的意見，解說到這邊，謝謝大家 (我不是專門研究病毒的，很抱歉)
※禁止轉載※
最近沒投2330，失望啊....
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 118.169.10.133 (臺灣)
※ 文章網址: https://webptt.com/m.aspx?n=bbs/nCoV2019/M.1587292651.A.060.html
※ 編輯: naieeliao (118.169.10.133 臺灣), 04/19/2020 18:40:37
1F：推 warrenchen: 感謝彙整101.137.108.236 04/19 18:42
2F：→ KunioFAN: 推整理~ 125.224.99.3 04/19 18:44
3F：→ soulgem: 想到上次看到的 graphen.ai 還在畫 0.0a 140.112.30.49 04/19 18:44
當我看到某些突變點都被抓出來的時候，表示他根本沒有篩選掉定序的誤差
當誤差為 0.01% 的時候，這隻病毒序列約 29k，就有3個誤差囉
以誤差 0.01% 拿去分析是不行的，需要用量去彌補誤差
而不是把誤差呈現在一張圖上，看起來眼花撩亂
5F：推 james732: 推 39.10.62.205 04/19 18:48
6F：推 mesenchymal: 感謝分析 67.87.236.128 04/19 18:54
7F：推 ilccc: 感謝整理！202.133.193.204 04/19 18:56
8F：推 volkyrie: 感謝整理最近就想看看這個，那請問有台 1.160.47.69 04/19 19:02
9F：→ volkyrie: 灣病株也會像澳洲那樣各地來嗎? 1.160.47.69 04/19 19:02
Git裡面有pdf跟txt，可以直接從 complete 搜尋 "taiwan"
長庚有一個小小的株系，至於突變出現四個我猜測是定序的誤差 (primer/機器)
其他大多來自國外
10F：推 snorker: push 61.223.168.12 04/19 19:09
11F：推 silentence: 要多大變異才能排除定序的偏誤呢@@?180.204.227.120 04/19 19:12
簡單說，當你把特定突變全部隱藏，而資料的亂度顯著下降
之前沒有仔細篩的時候，發現病毒突變好快，就跟graphen那張圖一樣亂
可是graphen連完整度 < 2% 的序列都不放過，所謂 GIGO
12F：推 aprilsheep: 公司名稱都打出來了，還甩鍋，就是中 223.140.197.11 04/19 19:12
13F：→ aprilsheep: 資呀，歐洲國家買到出問題的公司都會 223.140.197.11 04/19 19:12
14F：→ aprilsheep: 報名字推推 223.140.197.11 04/19 19:12
15F：→ aprilsheep: 抱歉上方為舊推錯 223.140.197.11 04/19 19:13
16F：推 becseger: 感謝專家解說！ 67.148.60.28 04/19 19:19
17F：→ SRadiant: 中國從很早就開始刻意不上傳序列了，不 111.82.201.234 04/19 19:29
18F：→ SRadiant: 然應該會更精采 111.82.201.234 04/19 19:29
19F：→ muskdeer: 看嘸...0.0? 1.160.137.252 04/19 19:32
20F：推 k879311: 這個要看出共同免疫並不是太容易，不過這 111.243.93.128 04/19 19:43
21F：→ k879311: 個可以提供很好的防疫方向。像是台灣找 111.243.93.128 04/19 19:43
22F：→ k879311: 不出感染源的，是否已經有台灣株了？ 111.243.93.128 04/19 19:43
DNA/RNA的序列不適合拿去推論蛋白質層面的東西，免疫是這部分無法分析的地方
23F：推 bunjie: 感謝分析和解說182.155.179.146 04/19 19:44
24F：→ kuma660224: 這跟免疫沒有關係啦 1.200.218.175 04/19 21:07
25F：→ kuma660224: 是看病毒傳往那個方向 1.200.218.175 04/19 21:07
26F：→ kuma660224: 例如台灣前期流行的都中國來源 1.200.218.175 04/19 21:08
27F：→ kuma660224: 後期來自歐美 而美國的來自歐洲 1.200.218.175 04/19 21:08
28F：推 kuma660224: 峈첑u你沒守好哪個方向 1.200.218.175 04/19 21:10
29F：→ kuma660224: 源頭武漢 但多數不是武漢直接去 1.200.218.175 04/19 21:11
30F：→ kuma660224: 源頭還好守 豬隊友才很難守 1.200.218.175 04/19 21:12
31F：推 abow0704: 好強的整理 223.138.97.30 04/19 21:15
32F：推 kennetharies: 推 111.82.149.74 04/19 21:21
33F：推 silentence: 長知識了 謝謝180.204.227.120 04/19 21:32
34F：推 clouse1234: 推 122.117.56.229 04/19 21:35
35F：推 arthurwang: 推整理 59.115.114.196 04/19 21:36
36F：推 ioupoiu: 垃圾支那 1.200.207.127 04/19 21:45
37F：推 MosDonalds: 推 49.216.169.230 04/19 21:57
38F：推 avans: 推研究 謝謝說明! 1.161.135.174 04/19 22:05
39F：推 fromwilda: 推111.242.148.127 04/19 22:39
40F：推 none024: 謝謝整理 73.43.106.240 04/19 22:44
※ 編輯: naieeliao (118.169.10.133 臺灣), 04/19/2020 23:03:54
41F：推 ThomasHSNU: 推114.136.195.235 04/19 23:29
42F：推 hotanger: 有問過張教授嗎？ 220.136.3.41 04/20 07:33
43F：推 soulgem: 推 (果然還是有很需要領域知識的地方) 140.112.30.49 04/20 07:50
44F：推 skycat2216: 知識推 101.136.131.89 04/20 08:02
45F：推 kisscat1023: 推專業 101.10.22.248 04/20 08:28
46F：推 windfolg: 推 223.141.30.65 04/20 10:00
47F：推 joann63: 推專業! 175.180.66.131 04/20 12:38
48F：推 ausagi: 推223.136.143.163 04/20 14:51