マンハッタンプロット（Manhattan plot）は、数百万件に及ぶ関連解析テストの結果を視覚的に要約する強力なツールですが、ゲノムワイド関連解析（GWAS）自体の不備を補正することはできません。集団構造化（Population structure）、表現型の定義、インピュテーション（Imputation）の精度、血縁関係、アレル頻度（Allele frequency）によるフィルタリング、そして多重比較の補正といった重要な課題は、プロット作成段階ではなく、事前の解析段階で解決しておく必要があります。図の作成は、十分に吟味された解析結果テーブルのみに基づいて行ってください。
SciDrawのManhattan Plot Generatorは、染色体上の位置と -log10(p) をマッピングし、作成した図をSciVisにシームレスに引き渡して微調整することができます。用意されているテンプレート例には、基本レイアウト、ゲノムワイド表示、アノテーション付き、特定領域のハイライト、示唆的なしきい値（Suggestive threshold）付き、そしてマイアミプロット（Miami plot）風レイアウトが含まれています。
マンハッタンプロット作成における一般的な誤り
- 未検証のp値のプロット: 数値以外の値、ゼロ、負の値、または1を超える値が含まれていると、対数変換が正常に機能しないか、上流の解析プロセスに問題があることを示しています。
- 染色体のアルファベット順ソート: 単純にソートすると、染色体10が染色体2の前に配置されたり、性染色体の順序が乱れたりする原因になります。
- 一律のしきい値の盲信:
5 × 10⁻⁸はヒトGWASにおける一般的な慣習的しきい値ですが、すべての生物種、マーカーセット、表現型、あるいは解析手法に適用できる万能のルールではありません。 - すべての有意な変異へのラベル付与: 同一遺伝子座内で連鎖不平衡にある複数の変異すべてにラベルを付けると、文字が重なり合って判読不能になります。
- ピークを「原因遺伝子」と直ちに解釈すること: 関連解析で特定されるのは、さらなる検証が必要な「領域」であり、それ自体が直接的な因果変異、原因遺伝子、または生物学的メカニズムを証明するものではありません。
マンハッタンプロットに必要なデータ構造
解析対象となった変異ごとに1行のデータを作成します。
| 列名 | 役割 | 例 |
|---|---|---|
snp | 変異の識別子（ID） | rs100102 |
chromosome | 染色体番号またはコンティグ名 | 3 |
position_bp | 塩基対（bp）位置 | 18420351 |
p_value | 変換前の関連解析p値 | 2.4e-9 |
trait | 表現型またはシリーズのラベル（任意） | Trait_A |
annotated_locus | 遺伝子座のラベル（任意） | Locus 3q |
具体的な構成については、基本的なマンハッタンプロット用CSVまたはアノテーション付きピーク用CSVを参照してください。p値は必ず数値形式とし、強いシグナルを正しく識別できるよう十分な精度（桁数）を維持してください。
ステップ 1: 検証済み関連解析テーブルの出力
適切な遺伝統計学ワークフローを用いて、品質管理（QC）と関連解析モデリングを完了させてください。プロット用テーブルには、表示対象となる最終的な変異のみを含め、ゲノムビルド、染色体表記、および塩基対座標に不整合がないことを確認します。使用したゲノムビルドは、図のキャプションや論文のMethodsセクションに必ず明記してください。
ステップ 2: CSVまたはExcelデータのアップロード
Manhattan Plot Generatorを開き、テンプレート例を選択するか、ご自身のCSV、TSV、TXT、XLSX、XLSファイルをアップロードします。Excelファイルの場合は、最初のワークシートが読み込まれます。染色体と位置情報が各ゲノム座標を正しく示していること、およびすべてのp値が 0 < p ≤ 1 の範囲に収まっていることを確認してください。
ステップ 3: 変換処理の実行と検証
y軸は -log10(p) を表すため、p値が小さい（有意性が高い）ほどプロットの上方に配置されます。例えば、p = 10⁻⁸ は 8 に変換されます。いくつかのデータポイントを手動で計算してプロットと照合し、染色体がゲノム上の順序通りに並んでいること、および隣接する染色体を区別するために交互の配色が正しく適用されていることを確認します。
ゲノムワイド有意水準として頻繁に用いられる 5 × 10⁻⁸ のラインは、-log10(p) 軸上では約 7.30 に相当します。実際の研究においては、テスト回数や独立性、研究固有のプロトコル、あるいは分野ごとの慣習に基づいて、適切なしきい値を設定する必要があります。
ステップ 4: すべての点ではなく、代表的な遺伝子座に絞ってラベルを付与する
センチネル変異（Sentinel variant）や独立して定義された遺伝子座にアノテーションを付与し、その選択基準を明記します。連鎖不平衡（LD）によって近接する変異群は一つのピークを形成しますが、これらを個別の独立した発見として提示すべきではありません。
探索的な文脈を示すために「示唆的なしきい値（Suggestive threshold）」を設ける場合は、主要な有意水準のラインとは異なるスタイル（破線など）で描き、凡例やキャプションで両者の違いを明確に説明してください。
その後、SciVisに移行して、しきい値ライン、染色体ラベル、アノテーションの引き出し線、配色、および全体的なレイアウトを微調整します。見栄えを良くするという理由だけで、データポイントの位置を動かしたり数値を書き換えたりしてはなりません。
マンハッタンプロットの解釈方法
プロット上の各点は、テストされた個々の変異を表します。横軸はゲノム上の位置を示し、縦軸は統計的証拠の強さを -log10(p) で示しています。垂直にそびえ立つ高いドットの集まり（ピーク）は、互いに関連（多くの場合、連鎖不平衡による相関）を持つ変異が集まる領域を示しています。
ある遺伝子座に注目する前に、事前に設定したしきい値、効果量（Effect size）とその不確実性、アレル頻度、インピュテーションやジェノタイピングの品質、ゲノムインフレーション因子（Genomic inflation factor）、局所的な関連パターン、再現性の証拠、そして解析プロトコル全体を再確認してください。プロット単体では、これらの品質管理プロセスを証明することはできません。
マンハッタンプロット、マイアミプロット、局所関連プロットの使い分け
- マンハッタンプロット: 単一の関連解析結果について、ゲノムワイドな全体像を俯瞰する際に使用します。
- マイアミプロット: 共通の染色体軸を挟んで、上下に対称となるように2つのゲノムワイド解析シリーズを並べて比較する際に使用します。
- 局所関連プロット（Regional association plot）: 特定の遺伝子座に焦点を当て、局所的な連鎖不平衡や周辺遺伝子との位置関係を詳細に確認する際に使用します。
- ボルケーノプロット: 発現変動解析（Differential expression）などで、効果量と有意水準の関係を可視化する際にVolcano Plot Generatorを使用します。こちらはゲノム上の位置情報は保持されません。
よくある質問（FAQ）
SciDraw上でGWASの解析自体を実行できますか？
いいえ。SciDrawは、事前に作成された関連解析結果テーブルを可視化するツールです。品質管理、モデリング、集団構造化や血縁関係の補正、および統計的推論に関する意思決定は、すべてご自身の解析ワークフロー内で行ってください。
すべての研究で 5 × 10⁻⁸ をしきい値として使用すべきですか？
いいえ。これはヒトのGWASで広く用いられている慣習的なしきい値ですが、妥当なしきい値は、研究デザイン、対象生物、マーカーセット、テスト間の独立性、および解析計画によって異なります。
p値がゼロ（0）と出力されたデータはどう処理すべきですか？
元の解析出力に戻り、対数表記（log p-value）をサポートしている場合はそれを使用するなど、十分な数値精度を保った形式でエクスポートし直してください。ゼロを任意の値に置き換えることは、統計的証拠を歪める可能性があります。
変異IDの代わりに遺伝子名をラベルとして表示できますか？
マッピング方法を明記した上で、近傍遺伝子や候補遺伝子のアノテーションを付与することは可能ですが、物理的な近接性が因果関係を証明するわけではありません。代表的な変異（Sentinel variant）や遺伝子座の定義は、キャプションや補足資料に必ず残しておいてください。
図を作成する
まずは Manhattan Plot Generator にアクセスし、解析記録と照らし合わせながら入力データとしきい値を検証してください。関連する図のタイプについては、科学データ可視化ガイドもあわせてご参照ください。
