曼哈頓圖 Manhattan Plot
曼哈頓圖本質上是一個散點圖，用於顯示大量非零大范圍波動數值，最早應用於全基因組關聯分析(GWAS)研究展示高度相關位點。它得名源於樣式與曼哈頓天際線相似(如下圖）。
近幾年，在宏基因組領域，尤其是差異OTU結合分類學結果，采用Manhattan plot展示有非常好的效果，倍受推崇。
曼哈頓圖優點
大數據中，即展示數據全貌，又能快速找到目標基因或OTU，同時可知目標的具體位置和分類、顯著程度等信息。絕對高端大氣，而且還有內涵。
數據坐標軸介紹
以上圖GWAS研究結果為例：
- X軸為染色體編號，且每個基因組SNP位點沿染色體序列排列；在16S擴增子或宏基因組中則為OTU按Taxonomy某一級別排序。
- Y軸為該位點相關的統計顯著性Pvalue值，由於pvalue值范圍是從0-1，且越小越好，直接展示非常密集於0附近，很難區分。如何使越近0的顯著數值變大，且而容易區分開，log10變換是非常好的方法，直接把關注的高顯著性(Pvalue趨近零)值高位顯示，遠離整體，目標一目了然。
- 圖中水平線一般為設定的不同顯著性水平閾值，方便讀出每個點的顯著性水平；或只添加一條顯示性閾值，高於則顯著。
曼哈頓圖繪制工具
散點圖，自然還是R語言，ggplot2可以畫的非常漂亮。
看圖實戰(Result)
示例1. 雙曼哈頓圖展示WT和mutant間差異富含OTU分布在那些菌目
這篇文章分析了百脈根根瘤的微生物組成，同時在根瘤缺失突變體條件下發現根和根際微生物組均有較大差異的變化
圖5.A/B 曼哈頓圖展示野生型，突變體根相對於根際土顯著差異的OTU類型
1. 圖中元素解釋
- X軸標簽“OTU… respect to rhizosphere”表示：根際土壤作為背景對照，計算富集的OTU；
- X軸OTU按分類學目水平(order)字母順序排列顯示，由於數量太多，不顯示OTU編號標簽反而更美觀；
- Y軸為-log10(Pvalue)；將pvalue轉換為越顯著越大，便於觀察；
- 主圖區的每個圓點或圈代表1個OTU，大小代表其相對豐度；其中存在顯著富集OTU的目中所有OTU用彩色實心圓點顯示，並添加灰度背景，且該目的名稱標注於圖頂部；目中內無顯著富集OTU的目為空心灰點，且背景為白色。
2. 圖表結果：兩個曼哈頓圖展示WT和mutant間差異富含OTU分布在那些菌目；而且與野生型相比，在突變體中許多顯著富集的菌目消失；
3. 經驗和技巧：單曼哈頓圖顯示顯著富集的OTU已經信息非常豐度；采用曼哈頓圖展示兩中組差異的OTU，讓讀者自己去比較差異，反而更突出結果的顯著差別。分類學注釋級別選擇目，找到了一些差別的類，要保證這些類即不能太多，也不能太少，才便於傳遞給讀者工作即全面、又細致的印像。
示例2
這是我自己畫的一個樣式，對上圖的樣式做了一些改進，展示一個基因敲除突變體(KO/mutant)與野生型(WT)細菌組的比較；
1. 圖中元素解釋
- X軸為OTU，按分類學門水平字母排序；
- Y軸兩組比較的Pvalue值，取loge(P)，即自然對數轉換；
- 圖中點的大小代表該OTU的相對豐度，取log2(CPM)對數，即2的對數；CPM為count per million的縮寫，和RPM類似，都是百萬分數；
- 圖中點顏色代表分類學門類型，便於從門水平找規律；
- 圖中點的形狀標注了其變化的類型，是上調enriched(正實心三角)，還是下調depleted(倒空心三角)，還是沒有顯著差異變化nosig(實心圓點)；
2. 圖表結果：展示了KO突變體基因型相較WT有較明顯的細菌組變化，尤其是放線菌門上調較多，變型菌門上調和下調都很多，但上調的更顯著；
3. 圖表經驗：從門水平先看整體規律，再一步步往綱、目、科、屬去找規律的具體細節；用形狀區分上調或下調，讓結果更清楚。