Как построить график Манхэттен (Manhattan Plot) для данных GWAS

• Визуализация непроверенных значений p-value: наличие нечисловых, нулевых, отрицательных значений или значений больше единицы нарушает математическое преобразование и указывает на ошибки в исходных расчетах.
• Алфавитная сортировка хромосом: при таком подходе хромосома 10 может оказаться перед хромосомой 2, а половые хромосомы займут некорректные позиции.
• Слепое применение универсального порога: общепринятый порог значимости для GWAS, равный 5 × 10⁻⁸, не является абсолютным правилом для любого организма, набора маркеров, фенотипа или модели анализа.
• Добавление подписей ко всем значимым вариантам: высокая плотность сцепленных вариантов в одном локусе при избыточном аннотировании превращает график в нечитаемое нагромождение текста.
• Интерпретация пика как причинного гена: обнаружение ассоциации лишь указывает на перспективный геномный район для дальнейшего изучения, но само по себе не определяет причинный вариант, ген или биологический механизм.

Какие данные необходимы для построения графика?

Для визуализации требуется таблица, где каждому протестированному варианту соответствует одна строка:

Вы можете ознакомиться со структурой данных в базовом CSV-файле или в файле с аннотированными пиками. Следите за тем, чтобы значения p-value оставались строго числовыми и сохраняли достаточную точность для корректного отображения сильных сигналов.

Шаг 1: Экспорт проверенной таблицы ассоциаций

Выполните контроль качества данных и статистическое моделирование ассоциаций в специализированном биоинформатическом конвейере. Итоговая таблица должна содержать только те варианты, которые планируется отобразить на графике, с единой сборкой генома, унифицированными обозначениями хромосом и точными физическими координатами. Обязательно укажите использованную сборку генома в подрисуночной подписи или в разделе методов вашей работы.

Шаг 2: Загрузка данных в формате CSV или Excel

Перейдите в Manhattan Plot Generator, выберите подходящий демонстрационный шаблон или загрузите собственный файл в формате CSV, TSV, TXT, XLSX или XLS. При импорте файлов Excel система считывает первый рабочий лист. Убедитесь, что столбцы хромосом и координат однозначно определяют положение каждого маркера в геноме, а все значения p-value строго укладываются в диапазон 0 < p ≤ 1.

Шаг 3: Генерация и проверка преобразования данных

Ось Y на графике отображает величину -log10(p), благодаря чему наиболее статистически значимые точки (с наименьшими p-value) оказываются на самом верху. Например, значение p = 10⁻⁸ преобразуется в 8. Проверьте несколько точек вручную и убедитесь, что хромосомы расположены в правильном геномном порядке, а чередующиеся цвета используются исключительно для визуального разделения соседних хромосом.

Стандартный порог полногеномной значимости 5 × 10⁻⁸ соответствует значению примерно 7.30 по шкале -log10(p). В зависимости от специфики вашего исследования, количества тестов и их взаимной зависимости может потребоваться расчет индивидуального порога значимости, обусловленного протоколом анализа или принятыми в конкретной области стандартами.

Шаг 4: Аннотирование локусов, а не отдельных точек

Наносите подписи только для ключевых (sentinel) вариантов или независимо определенных локусов, обязательно поясняя критерии их отбора. Из-за неравновесия по сцеплению близкорасположенные варианты часто формируют единый пик; их не следует представлять как независимые открытия.

Если для демонстрации поискового этапа анализа вам необходимо показать порог предварительной значимости (suggestive threshold), визуально отделите его от основного порога значимости и опишите оба критерия в легенде или подрисуночной подписи.

Перейдите в редактор SciVis, чтобы настроить положение линий порогов, подписи хромосом, линии-выноски для аннотаций, цветовую палитру и общую композицию. Изменять положение точек или редактировать их исходные значения в угоду визуальной привлекательности категорически не допускается.

Как интерпретировать график Манхэттен

Каждая точка на графике представляет собой протестированный генетический вариант. Ее положение по горизонтали указывает на локализацию в геноме, а по вертикали — на статистическую силу ассоциации, выраженную как -log10(p). Высокие вертикальные кластеры точек (пики) указывают на геномные области, содержащие ассоциированные и, как правило, сильно сцепленные между собой маркеры.

Прежде чем делать выводы о значимости конкретного локуса, детально проанализируйте установленный порог значимости, величину эффекта и доверительные интервалы, частоту аллелей, качество импутации или генотипирования, показатель геномной инфляции (lambda), характер локальной ассоциации, результаты репликации в независимых выборках и общий протокол исследования. Сам по себе график не заменяет этих обязательных этапов верификации.

График Манхэттен, график Майами или локальный график ассоциаций?

• Используйте график Манхэттен для общего полногеномного обзора результатов одного анализа ассоциаций.
• Используйте график Майами (Miami plot) для сопоставления двух полногеномных серий данных, расположенных зеркально (сверху и снизу) относительно общей оси хромосом.
• Используйте локальный график ассоциаций (regional association plot) для детального анализа конкретного локуса с учетом структуры неравновесия по сцеплению и расположения генов.
• Используйте Volcano Plot Generator для сопоставления величины эффекта и статистической значимости в исследованиях дифференциальной экспрессии генов или аналогичных анализах; этот тип графиков не сохраняет информацию о геномных координатах.

Часто задаваемые вопросы (FAQ)

Выполняет ли SciDraw расчеты GWAS?

Нет. Инструмент предназначен исключительно для визуализации предварительно рассчитанной таблицы результатов ассоциативного анализа. Контроль качества, статистическое моделирование, учет популяционной структуры или родства, а также формулирование статистических выводов осуществляются исследователем на этапе первичного анализа данных.

Обязательно ли использовать порог 5 × 10⁻⁸ в любом исследовании?

Нет. Этот порог является общепринятым стандартом для GWAS человека, однако корректный и научно обоснованный порог всегда зависит от дизайна эксперимента, исследуемого организма, плотности маркеров, характера сцепления и общего плана статистического анализа.

Что делать со значениями p-value, которые в выводе программ указаны как нулевые?

Вернитесь к исходным результатам анализа и экспортируйте данные с более высокой точностью представления чисел с плавающей запятой (например, в виде логарифмов p-value, если это поддерживается ПО). Замена нулевых значений произвольными числами может исказить реальную статистическую картину.

Можно ли подписывать названия генов вместо идентификаторов вариантов?

Вы можете указать близлежащий ген или ген-кандидат, если в работе четко описан метод картирования. Однако помните, что физическая близость гена к маркеру не доказывает причинно-следственную связь. Обязательно приведите информацию о ключевом варианте (sentinel variant) или границах локуса в подписи к рисунку или в дополнительных материалах.

Создать график

Перейдите в Manhattan Plot Generator, загрузите данные, сверьте пороги значимости с протоколом вашего исследования и ознакомьтесь с нашим руководством по научной визуализации данных для выбора других типов графиков.