Визуализируйте данные NGS Используя приложение Genomics Viewer

Приложение Genomics Viewer позволяет вам просмотреть и исследовать интегрированные геномные данные со встроенной версией Integrative Genomics Viewer (IGV) [1][2]. Геномные данные включают выравнивания чтения NGS, варианты генома и сегментированные данные о номере копии.

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

Откройте приложение

В командной строке ввести genomicsViewer. В качестве альтернативы кликните по значку приложения на вкладке Apps. Приложение требует интернет-соединения.

По умолчанию, Человек загрузок приложения (GRCh38/hg38) как ссылочная последовательность и Гены Refseq как файл аннотации. В приложении существует две основных панели. Левая панель является панелью Tracks, и правильная панель является встроенным веб-приложением IGV. Панель Tracks является областью только для чтения, отображающей имена дорожки, имена исходного файла и типы дорожки. Панель Tracks обновляется соответственно, когда вы конфигурируете дорожки во встроенном приложении IGV.

Default view of the Genomics Viewer app. The toolstrip is at the top. The Tracks Panel is on the left. The embedded integrative genomics viewer IGV is on the right.

Кнопка Reset восстанавливает приложение к представлению по умолчанию с двумя дорожками (HG38 с Генами Refseq) и удаляет любые другие существующие дорожки. Перед сбросом можно сохранить текущее представление как сеанс (.json) файл и восстановление это позже.

Добавьте дорожки путем импортирования данных

Импортируйте ссылочную последовательность

Можно импортировать одну ссылочную последовательность. Ссылочная последовательность должна быть в файле FASTA. Выберите Import Reference на вкладке Home. Можно также импортировать соответствующий cytoband файл, который содержит цитогенетические данные G-соединения. Можно добавить локальные файлы или задать внешние URL. URL Должен запуститься или с https или с gs. Другие протоколы передачи файлов, такие как ftp, не поддерживаются.

Импортируйте данные о выравнивании чтения последовательности

Можно импортировать несколько наборов данных считанных данных о выравнивании последовательности. Данные о выравнивании должны быть файл CRAM или BAM. Не требуется, что у вас есть соответствующий индексный файл (.BAI или .CRAI) в том же месте как ваш BAM или файл CRAM. Однако отсутствие индексного файла сделает приложение медленнее.

Можно добавить файлы выравнивания чтения с помощью Add tracks from file и опций Add tracks from URL от кнопки Add tracks. Если вы задаете URL, URL должен запуститься или с https или с gs. Другие протоколы передачи файлов, такие как ftp, не поддерживаются.

Импортируйте аннотации функции и другие геномные данные

Можно импортировать несколько наборов аннотаций функции из нескольких файлов, которые содержат данные для одной ссылочной последовательности. Поддерживаемые файлы аннотации: .BED, .GFF, .GFF3, и .GTF.

Можно также импортировать структурные варианты (.VCF) и визуализировать генетические изменения, такие как вставки и удаления.

Можно просмотреть сегментированные данные о номере копии (.SEG) и количественные геномные данные (.WIG, .BIGWIG, и .BEDGRAPH), такие как peaks ChIP и покрытие выравнивания.

Можно добавить аннотацию и геномные файлы данных с помощью Add tracks from file и опций Add tracks from URL от кнопки Add tracks. Если вы задаете URL, URL должен запуститься или с https или с gs. Другие протоколы передачи файлов, такие как FTP, не поддерживаются.

Визуализируйте одно изменение нуклеотида цитохрома P450

Ген CYP2C19 является членом цитохрома семейство генов P450. Ферменты, произведенные из цитохрома гены P450, вовлечены в метаболизм различных молекул и химикатов в ячейках. Фермент CYP2C19 играет роль в усваивании по крайней мере 10 процентов обычно прописываемых лекарств [3]. Полиморфизмы в цитохроме p450 семейство могут вызвать неблагоприятные ответы препарата в индивидуумах. Одним примером одного изменения нуклеотида является rs4986893 в положении chr10:94,780,653 где G заменяется A. Этот аллельный вариант также известен как CYP2C19*3. Следующие шаги показывают, как визуализировать такое изменение приложения с помощью и низкого покрытия и высоких данных о покрытии.

Загрузите файл сеанса

В целях этого примера начните с файла сеанса, который имеет некоторые предварительно загруженные дорожки. Чтобы загрузить файл, нажмите Open. Перейдите к matlabroot\examples\bioinfo\, где matlabroot является папкой, где вы установили MATLAB®. Выберите rs4986893.json.

The tracks panel has three tracks, namely, hg38.fa sequence, NA18564 alignment data, and refseq genes annotation. The IGV shows the aligned reads graphically.

Сеанс содержит три дорожки:

  • Human (GRCh38/hg38) как ссылка

  • NA18564 как низкие данные о выравнивании покрытия

  • Гены Refseq

Низкие данные о выравнивании покрытия прибывают из розеточные ханьцы из Пекина, Китай. Демонстрационным ID является NA18564, и выборка была идентифицирована с мутацией CYP2C19*3 [4].

Исследуйте низкие данные о покрытии

В этом файле сеанса данные о выравнивании были сосредоточены вокруг местоположения мутации на гене CYP2C19.

  1. Кликните по оранжевой панели в зоне охвата, чтобы посмотреть на положение и информацию о распределении аллели.

    Image of aligned reads with overlapping context menu in IGV. The context menu shows the counts of A, C, G, T, N.

    Это показывает, что 71% чтений имеет G, в то время как 29% имеют в местоположении chr10:94,780,653. Эти данные являются низкими данными о покрытии и не могут показать все случаи этой мутации. Высокие данные о покрытии будут исследованы позже в примере.

    Закройте окно всплывающей подсказки.

  2. Можно настроить различные аспекты отображения данных в приложении. Например, можно изменить высоту дорожки, чтобы сделать больше комнаты для более поздних дорожек. Кликните по второму значку механизма. Выберите Set track height. Войдите 200.

    Image of the context menu. It has various options to change the appearance of the track, such as track color, track name, and track height.

    Для получения дополнительной информации на встроенном приложении IGV и его доступных параметрах, посетите здесь.

Исследуйте высокие данные о покрытии

Можно посмотреть на высокие данные о покрытии из той же выборки, чтобы видеть случаи этой мутации.

  1. Перейдите к Международному Демонстрационному веб-сайту Ресурса Генома.

  2. Ищите демонстрационный NA18564.

  3. Загрузите файл выравнивания Exome, который находится в .CRAM формат.

  4. Также загрузите соответствующий индексный файл, который находится в .CRAI формат. Сохраните файл в том же месте как источник .CRAM файл.

  5. Щелкните (+) значок на вкладке Home. Выберите загруженный .CRAM файл и нажимает Open.

    The tracks panel now shows the fourth track for the alignment data that was loaded. IGV also shows the additional track for the added alignment data.

    Высокие данные о покрытии появляются как track3. Можно теперь видеть много случаев мутации в нескольких чтениях.

  6. Кликните по оранжевой панели в зоне охвата, чтобы видеть распределение аллели. Это показывает, что G заменяется почти в 50% времени.

    Image of the context menu which shows the counts for A, C, G, T, and N. The counts for A is 79 (49%) and the counts for G is 82 (51%). Other counts are zero.

Ссылки

[1] Робинсон, J., Х. Торволдсдоттир, В. Винклер, М. Гуттман, E. Высаживающийся на берег, Г. Гетц, Й. Месиров. 2011. Интегральный Genomics Viewer. Биотехнология природы. 29:24–26.

[2] Thorvaldsdóttir, H., Дж. Робинсон, Й. Месиров. 2013. Интегральный Genomics Viewer (IGV): высокоэффективная визуализация данных о геномике и исследование. Брифинги в Биоинформатике. 14:178–192.

Смотрите также

| |

Для просмотра документации необходимо авторизоваться на сайте