exponenta event banner

chromosomeplot

Сюжетная хромосомная идеограмма с рисунком G-banding

Синтаксис

chromosomeplot(CytoData)
chromosomeplot(CytoData, ChromNum)
chromosomeplot(CytoData, ChromNum, ...,'Orientation', OrientationValue, ...)
chromosomeplot(CytoData, ChromNum, ...,'ShowBandLabel', ShowBandLabelValue, ...)
chromosomeplot(CytoData, ChromNum, ...,'AddToPlot', AddToPlotValue, ...)
chromosomeplot(..., 'Unit', UnitValue, ...)
chromosomeplot(..., 'CNV', CNVValue, ...)

Аргументы

CytoDataОдно из следующих действий:
  • Символьный вектор или строка, задающая файл, содержащий цитогенетические данные G-banding (в единицах bp), такой как текстовый файл идеограммы NCBI или текстовый файл цитобанда браузера генома UCSC.

  • Структура, содержащая цитогенетические данные G-banding (в п.в.) в следующих областях:

    • ChromLabels

    • BandStartBPs

    • BandEndBPs

    • BandLabels

    • GieStains

Совет

Используйте cytobandread для создания структуры, используемой для CytoData.

ChromNumСкалярный или символьный вектор или строка, задающая одну хромосому для построения графика. Допустимыми записями являются целые числа, 'X', и 'Y'.

Примечание

Настройка ChromNum кому 0 строит идеограммы для всех хромосом.

OrientationValueСимвольный вектор или строка или число, указывающее ориентацию идеограммы одной хромосомы, указанной ChromNum. Варианты: 'Vertical' или 1 (по умолчанию) и 'Horizontal' или 2.
ShowBandLabelValueУправление отображением меток областей данных (например, q25.3) при построении единой хромосомной идеограммы, указанной ChromNum. Варианты: true (по умолчанию) или false.
AddToPlotValueИмя переменной оси фигуры, к которой добавляется одиночная хромосомная идеограмма, указанная ChromNum.

Примечание

Если это свойство используется для добавления идеограммы к графику геномных данных в единицах, отличных от bp, используйте 'Unit' для преобразования данных идеограммы в соответствующие единицы.

Совет

Перед печатью фигуры, содержащей добавленную хромосомную идеограмму, измените фон на белый, выдав следующую команду:

set(gcf,'color','w') 
UnitValueЦелое число, указывающее единицы измерения (пары оснований, пары оснований килограммов или пары оснований мега) для начального и конечного геномных положений. Эта единица используется в подсказке данных, отображаемой при наведении курсора на хромосомы на идеограмме. Этот блок также может использоваться при использовании 'AddToPlot' для добавления идеограммы к графику в единицах, отличных от bp. Варианты: 1 (bp), 2 (kb), или 3 (мб). По умолчанию: 1 (п.н.).
CNVValueУправляет отображением данных отклонения номера копии (CNV), предоставленных CNVValue, выровненный по хромосомной идеограмме. Выигрыши отображаются зеленым цветом справа или выше идеограммы, а потери - красным цветом слева или ниже идеограммы. CNVValue - структурный массив, содержащий четыре поля, описанные в таблице ниже.

Описание

chromosomeplot(CytoData) строит график идеограммы всех хромосом, используя информацию из CytoDataструктура, содержащая цитогенетические данные G-banding (в единицах bp) или символьный вектор или строку, задающую файл, содержащий цитогенетические данные G-banding (в единицах bp), такой как текстовый файл идеограммы NCBI или текстовый файл цитобанда UCSC Genome Browser. Полосы G различают различные области хромосомы. Например, для идеограммы Homo sapiens возможными диапазонами G являются:

  • gneg - белый

  • gpos25 - светло-серый

  • gpos50 - средняя серая

  • gpos75 - темно-серый

  • gpos100 - черный

  • acen - красный (центромер)

  • stalk - область с отступом (область с повторами)

  • gvar - светло-синий

Более темные полосы богаты AT, в то время как более светлые полосы богаты GC.

chromosomeplot(CytoData, ChromNum) строит график идеограммы одной хромосомы, указанной ChromNum.

chromosomeplot(..., 'PropertyName', PropertyValue, ...) требования chromosomeplot с необязательными свойствами, использующими пары имя/значение свойства. Можно указать одно или несколько свойств в любом порядке. Каждый PropertyName должен быть заключен в одинарные кавычки и не учитывать регистр. Эти пары имя/значение свойства следующие:

chromosomeplot(CytoData, ChromNum, ...,'Orientation', OrientationValue, ...) определяет ориентацию идеограммы одной хромосомы, указанной ChromNum. Варианты: 'Vertical' или 1 (по умолчанию) и 'Horizontal' или 2.

Примечание

При построении идеограммы всех хромосом ориентация всегда вертикальная.

chromosomeplot(CytoData, ChromNum, ...,'ShowBandLabel', ShowBandLabelValue, ...) отображает метки полосы (например, q25.3) при построении одной хромосомной идеограммы, указанной ChromNum. Варианты: true (по умолчанию) или false.

chromosomeplot(CytoData, ChromNum, ...,'AddToPlot', AddToPlotValue, ...) добавляет однохромосомную идеограмму, заданную ChromNum, к оси фигуры, указанной AddToPlotValue.

Примечание

Если это свойство используется для добавления идеограммы к графику геномных данных в единицах, отличных от bp, используйте 'Unit' для преобразования данных идеограммы в соответствующие единицы.

Совет

Перед печатью фигуры, содержащей добавленную хромосомную идеограмму, измените фон на белый, выдав следующую команду:

set(gcf,'color','w') 

chromosomeplot(..., 'Unit', UnitValue, ...) определяет единицы измерения (пары оснований, пары оснований в килограммах или мегапараты оснований) для начального и конечного геномных положений. Эта единица используется в подсказке данных, отображаемой при наведении курсора на хромосомы на идеограмме. Этот блок также может использоваться при использовании 'AddToPlot' для добавления идеограммы к графику в единицах, отличных от bp. Варианты: 1 (bp), 2 (kb), или 3 (мб). По умолчанию: 1 (п.н.).

chromosomeplot(..., 'CNV', CNVValue, ...) отображает данные отклонения количества копий (CNV), предоставленные CNVValue, выровненный по хромосомной идеограмме. Выигрыши отображаются зеленым цветом справа или выше идеограммы, а потери - красным цветом слева или ниже идеограммы. CNVValue является массивом структуры, содержащим следующие поля. Каждое поле должно содержать одинаковое количество элементов.

ОбластьОписание
Chromosome

Одно из следующих действий:

  • Числовой вектор, содержащий номер хромосомы, на которой расположен каждый CNV.

    Примечание

    Для половой хромосомы, X, используйте N, где N = число аутосом + 1. Для половой хромосомы, Y, используйте М, где М = число аутосом + 2. Например, для Homo sapiens используется 23 для X и 24 для Yи для мускулов Mus (лабораторная мышь), используйте 20 для X и 21 для Y.

  • Матрица символов, содержащая номер хромосомы, на которой расположен каждый CNV.

    Примечание

    Использование символьного массива позволяет использовать символы X и Y (вместо цифр) для половых хромосом. Однако все элементы массива должны иметь одинаковую ширину, что может потребовать добавления пробелов в некоторые векторы символов. Например:

    [' 1'; ' 2'; '10'; ' X']

    Или вы можете использовать char функция с клеточным массивом для создания символьного массива хромосомных чисел и букв. Например:

    char({'1', '2', '10', 'X'})

CNVType

Числовой вектор, содержащий тип каждого CNV. 1 (потеря) или 2 (выигрыш).

Start

Числовой вектор, содержащий исходное геномное положение каждого CNV. Единицы измерения должны состоять из базовых пар.

End

Числовой вектор, содержащий конечное геномное положение каждого CNV. Единицы измерения должны состоять из базовых пар.

Примеры

свернуть все

Прочитайте информацию о цитогенетической связке для Homo sapiens в структуру.

hs_cytobands = cytobandread('hs_cytoBand.txt')
hs_cytobands = struct with fields:
     ChromLabels: {862x1 cell}
    BandStartBPs: [862x1 int32]
      BandEndBPs: [862x1 int32]
      BandLabels: {862x1 cell}
       GieStains: {862x1 cell}

Постройте всю хромосомную идеограмму.

chromosomeplot(hs_cytobands);
title('Human Karyogram')

Можно отобразить идеограмму определенной хромосомы, щелкнув ее правой кнопкой мыши на графике, а затем выбрав команду «Показать на новой фигуре» > «Вертикальная» или «Горизонтальная».

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

chromosomeplot(hs_cytobands, 15, 'Orientation', 2, 'Unit', 2);

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

Загрузить данные CGH (aCGH) на основе массива из исследования клеточной линии Кориелла (Snijders, A. et al., 2001).

load coriell_baccgh

Используйте cghcbs функция анализа хромосомы 10 образца 3 (GM05296) данных aCGH и данных дисперсии числа возвращаемых копий (CNV) в структуре, S. График сегмента означает, что исходные данные относятся только к хромосоме 10 образца 3.

S = cghcbs(coriell_data,'sampleindex',3,'chromosome',10,...
           'showplot',10);
Analyzing: GM05296. Current chromosome 10

Figure contains an axes. The axes with title GM05296 - Chr 10 contains 4 objects of type line.

Используйте функцию хромосомеплота с опцией «addtoplot», чтобы добавить на график идеограмму хромосомы 10 для Homo sapiens. Поскольку график данных CNV из исследования клеточной линии Кориелла находится в КБ, используйте свойство «Unit» для преобразования данных идеограммы в КБ.

set(gcf,'color','w'); % Set the background of the figure to white.
chromosomeplot('hs_cytoBand.txt', 10, 'addtoplot', gca,...
               'Unit', 2);

Figure contains an axes. The axes with title GM05296 - Chr 10 contains 4 objects of type line.

Ссылки

[1] Снайдерс, А.М., Новак, Н., Сегрейвс, Р., Блэквуд, С., Браун, Н., Конрой, Дж., Гамильтон, Г., Хиндл, А.К., Хьюи, Б., Кимура, К., Ло, С., Мямбо, К., Палмер, Дж. Грей, Джей У., Джейн, А.Н., Пинкел, Д. и Альбертсон, Д.Г. (2001). Сборка микрочипов для измерения числа копий ДНК по всему геному. Генетика природы 29, 263-264.

См. также

|

Представлен в R2007b