Просмотр кода страницы Поиск пиков в данных ChiP-Seq

'''Поиск пиков''' ({{lang-en|Peak calling}}) — компьютерный метод поиска областей [[геном]]а, обогащённых выровненными ридами из данных [[ChIP-seq|ChiP-Seq эксперимента]]. Такие области называют пиками, поскольку они являются локальными максимумами покрытия генома ридами (то есть количества прочтений, выровненных в данной позиции). Определённые таким образом регионы считаются местами связывания белка, с которым проводился опыт. Если речь идёт о [[Факторы транскрипции|транскрипционном факторе]], то обогащённые области по сути являются сайтами связывания транскрипционного фактора. Определение областей пиков является важным шагом любого Chip-Seq эксперимента, поскольку именно на этом этапе происходит интерпретация данных — нужно отличить шумы от сигнала<ref>{{статья |заглавие=Genome-wide analysis of transcription factor binding sites based on ChIP-seq data |издание=[[Nature Methods]] |том=6 |номер=5 |страницы=829—834 |pmc=2917543 |ссылка=http://www.nature.com/nmeth/journal/v6/n11s/full/nmeth.1371.html |pmid=19160518 |doi=10.1038/nmeth.1246 |язык=en |тип=journal |автор=Valouev A., etal |месяц=9 |год=2008 |archivedate=2011-08-10 |archiveurl=https://web.archive.org/web/20110810065702/http://www.nature.com/nmeth/journal/v6/n11s/full/nmeth.1371.html }}</ref>.

== Обобщённый алгоритм ==
Операция нахождения пиков содержит в себе несколько этапов, необходимых для получения итогового списка пиков<ref>{{Cite web |url=http://epigenie.com/guide-peak-calling-for-chip-seq/ |title=Guide: Peak Calling for ChIP-Seq |access-date=2016-05-02 |archive-date=2016-05-01 |archive-url=https://web.archive.org/web/20160501140254/http://epigenie.com/guide-peak-calling-for-chip-seq/ |deadlink=no }}</ref>:

# '''Смещение ридов''' (важно в случае секвенирования с одного конца — single-end sequencing). Каждый рид, полученный после секвенирования, имеют длину 50-300 нуклеотидов и маркирует только один из концов фрагмента ДНК, очищенного [[Иммунопреципитация|иммунопреципитацией]]. Соответственно, он может не содержать в себе сайта связывания белка, если длина ДНК фрагментов больше длины рида, а белок был связан с центром фрагмента. На начальном этапе поиска идентифицируют области, обогащённые ридами. После этого находят отдельно максимумы обогащения для «+» и «-» цепей. Расстояние между максимумами считают средней длиной фрагментов ДНК. После этого все риды смещают на половину такой длины в сторону центра пика (то есть в сторону их 3’-конца), и определяют истинную область связывания. В случае секвенирования парных ридов (с обоих концов каждого фрагмента), эта операция не требуется, поскольку пара ридов точно маркирует каждый фрагмент, а значит сайт связывания заведомо находится между ними.
# '''Оценка уровня шума'''. Чаще всего для оценки уровня шума проводят контрольные Chip-Seq эксперименты, в которых стадия иммунопреципитации пропущена или используются неспецифические антитела (IgG-контроли). Однако существуют программы, способные оценить шум по данным одного эксперимента.
# '''Идентификация пиков.''' Чаще всего используется скользящее окно, которое сканирует геном и определяет регионы, в которых покрытие ридами превышает некий заданный порог, или отношение покрытия условного пика превышает окружающий его шум в заданное число раз.
# '''Оценка качества определённых пиков.''' Многие программы подсчитывают [[P-значение|P-value]] для каждого пика, его высоту и обогащения ридами относительно фона. Исходя из статистической значимости пиков можно сделать оценку ожидаемого уровня ложных предсказаний ([[:en:False discovery rate|FDR]]).
# '''Удаление артефактов.''' В данных Chip-Seq встречаются артефакты [[Полимеразная цепная реакция|ПЦР]], возникающие при пробоподготовке библиотек для секвенирования. При фильтрации пиков на предмет артефактов, удаляются такие из них, которые содержат только дуплицированные (то есть идентичные) риды, поскольку считается, что такие пики возникли из-за ПЦР. Также критически относятся к ситуациям, когда риды картированы резко асимметрично (только на 3’ или 5’ цепь).

== Параметры, влияющие на качество определения пиков ==
На качество определения пиков влияет множество параметров, начиная от способа картирования ридов на геном. Так, количество разрешённых при картировании несовпадений влияет на долю ридов, которые будут успешно картированы, а значит и на профиль покрытия генома. А операции по удалению ридов, способных выравниваться с несколькими местами в геноме, могут привести к потере пиков, расположенных в областях с повторами. Кроме того, в каждом из пяти перечисленных пунктов обобщённого алгоритма идентификации пиков существуют свои параметры, аккуратный подбор которых позволяет добиться максимальной селективности и специфичности<ref>{{cite doi|10.1371/journal.pone.0096303}}</ref><ref>{{cite doi|10.1093/bioinformatics/btu722}}</ref>.

== Программы для поиска пиков ==
К настоящему времени существуют десятки программ для определения пиков, отличающиеся по своей применимости и используемым алгоритмам. В данной статье будет разобрана самая популярная программа — [https://web.archive.org/web/20140921172022/http://liulab.dfci.harvard.edu/MACS/ MACS], которая хорошо подходит для поиска относительно узких пиков, характерных для сайтов связывания транскрипционных факторов. Также будет рассмотрена программа [https://web.archive.org/web/20160514012212/http://omictools.com/sicer-tool SICER], целью которой является поиск слабых и диффузных сигналов (то есть широких и слабых пиков), характерных для распределения модификаций [[Гистоны|гистонов]] или отражающих работающие [[Полимераза|полимеразы]]<ref>{{статья |заглавие=Identifying ChIP-seq enrichment using MACS |издание={{Нп3|Nature Protocols}} |том=7 |номер=9 |страницы=1728—1740 |doi=10.1038/nprot.2012.101 |язык=en |автор=Feng Jianxing, Liu Tao, Qin Bo, Zhang Yong, Liu Xiaole Shirley |число=29 |месяц=8 |год=2012 |тип=journal}}</ref>.

== MACS ==
'''MACS''' ({{lang-en|Model-based Analysis of Chip-Seq}}) — программа анализа данных ChiP-Seq эксперимента. Используемый программой метод хорошо подходит для поиска относительно узких пиков, характерных для сайтов связывания транскрипционных факторов. MACS является наиболее часто используемым алгоритмом для поиска пиков.

=== Стадии алгоритма ===
1. '''Моделирование размера сдвига ChiP-Seq ридов'''
[[Файл:Рис. 1. Модель MACS для ChiP-Seq.png|мини|'''Рисунок 1.''' ''Модель MACS для ChiP-Seq.Риды из случайной выборки 1000 модельных пиков, выровненные относительно центров распределений пиков «+» и «-» цепей. d — расстояние между максимумами «+» и «-» цепей в выравнивании.'']]
ChiP-Seq риды представляют собой концы фрагментов в ChiP-DNA библиотеке. По описанным выше причинам (см. Обобщённый алгоритм), они часто сдвинуты в сторону 5'-концов и размер этого сдвига неизвестен экспериментатору. Так как фрагменты ChiP-DNA равновероятно секвенированы с обоих концов, то плотность ридов вокруг истинного сайта связывания будет бимодальным: на «+» цепи плотность возрастает до сайта связывания, на «-» цепи — после (рис.1). MACS использует этот двойной паттерн для построения эмпирической модели размера сдвига ридов в сторону их 3’-конца для лучшего определения сайта связывания. Для первичного поиска регионов, богатых ридами, данный алгоритм использует 2 параметра: '''bandwidth''' (средняя длина фрагментов ДНК, получаемая при пробоподготовке, например, при облучении ультразвуком) и '''mfold''' (пороговое обогащение относительно уровня фона, используемого моделью, выше которого заданная область считается статистически обогащённой ридами). Скользя по геному «окном» размером 2*bandwidth алгоритм ищет регионы, mfold которых выше, чем mfold при случайном распределении ридов по геному. Затем MACS выбирает 1000 случайных, из определённых на предыдущем шаге областей, и отделяет в них риды «+» и «-» цепей. Если максимум плотности ридов «+» цепи лежит левее максимума для «-» цепи, то алгоритм выравнивает их относительно центров распределений. Расстояние между максимумами «+» и «-» цепей в выравнивании будет определяться как '''d''' (рис.1). MACS сдвигает все теги на расстояние d/2 в сторону 3'-конца для определения наиболее вероятного сайта взаимодействия ДНК с белком.
 
2. '''Поиск пиков'''

В экспериментах с контролем (см. Обобщённый алгоритм), MACS линейно увеличивает число контрольных ридов до общего количества ридов, полученных в ChiP-Seq эксперименте. В случае дуплицированных ридов, MACS способен эффективно удалять идентичные последовательности для снижения уровня шума.

Для статистической оценки значимости пиков программа использует предположение, что количество ридов, выровненных в данной позиции генома, распределено по Пуассону<ref>{{статья |заглавие=Genome-wide maps of chromatin state in pluripotent and lineage-committed cells |издание=Nature |номер=448 |страницы=553—560 |pmc=2921165 |ссылка=http://www.nature.com/nature/journal/v448/n7153/full/nature06008.html |pmid=17603471 |doi=10.1038/nature06008 |язык=en |тип=journal |автор=Mikkelsen T. S., etal |месяц=7 |год=2007 |archivedate=2011-08-13 |archiveurl=https://web.archive.org/web/20110813101643/http://www.nature.com/nature/journal/v448/n7153/full/nature06008.html }}</ref>. Преимуществом данной модели является то, что один параметр — '''λ<sub>BG</sub>''', оценённый для всего генома по контрольным точкам, определяет среднее значение и дисперсию распределения. После того, как MACS сдвинул каждый относящийся к обогащённому региону рид на d/2, алгоритм ещё раз сканирует геном, теперь уже «окном» шириной 2d и находит статистически значимые пики (используя приближение о Пуассоновском распределении числа ридов с параметром λ<sub>BG</sub> и уровнем значимости по умолчанию 10<sup>−5</sup>). Пересекающиеся пики сливаются, а участок их слияния называют '''summit'''. Он определяет локализацию сайта связывания.

В контрольных образцах часто наблюдаются флуктуации в распределении ридов. На это могут влиять многие факторы, такие как локальные структуры хроматина, шум, возникающий при амплификации и секвенировании и вариации в числе копий генома. Поэтому, MACS может использовать не только λ<sub>BG</sub>, оценённый по всему геному, но и динамический параметр λ<sub>local</sub>, который вычисляют для каждого пика как:

<math>\textstyle \lambda_{local} = \max\bigl(\lambda_{BG} , \lambda_{1k} , \lambda_{5k} , \lambda_{10k} \bigr) </math>

где λ<sub>1k</sub>, λ<sub>5k</sub> и λ<sub>10k</sub> — это параметры, определённые для «окон» шириной в 1 кб, 5 кб и 10 кб соответственно, центрированных относительно рассматриваемого пика . Таким образом, λ<sub>local</sub> отражает влияние местных флуктуаций.
[[Файл:Рисунок 2. Cравнение алгоритма MACS с алгоритмами ChIPSeq Peak Finder, FindPeaks и QuEST..png|мини|224x224пкс|'''Рисунок 2.''' ''Сравнение алгоритма MACS с алгоритмами ChIPSeq Peak Finder, FindPeaks и QuEST. Показаны подсчитанные FDR для таких факторов как FoxA1 и NRSF.'']]
Далее MACS использует λ<sub>local</sub> для подсчёта p-value каждого потенциального пика и удаляет ложно положительные результаты, вызванные шумом в покрытии ридами. Потенциальные пики с p-value ниже заданного (по умолчанию 10<sup>−5</sup>) отмечаются, а отношение между количеством ридов ChiP-Seq, выровненных в данной позиции и λ<sub>local</sub> называют '''fold enrichment'''. Для Chip-Seq эксперимента с контролем, MACS вычисляет ожидаемый уровень ложных предсказаний (FDR) для каждого определённого пика, используя процедуру, которая применялась и в ранних алгоритмах поиска, таких как MAT<ref>{{статья |заглавие=Model-based analysis of tiling-arrays for ChIP-chip |издание=[[Proceedings of the National Academy of Sciences|Proceedings of the National Academy of Sciences of the United States of America]] |номер=103 |страницы=12457—12462 |pmc=1567901 |pmid=16895995 |doi=10.1073/pnas.0601180103 |язык=en |тип=journal |автор=Johnson W. E., etal |месяц=8 |год=2006}}</ref> и MA2C<ref>{{статья |заглавие=Model-based analysis of two-color arrays (MA2C) |издание={{Нп3|Genome Biology|Genome Biol.||Genome Biology}} |номер=8 |страницы=R178 |pmc=2375008 |pmid=17727723 |doi=10.1186/gb-2007-8-8-r178 |язык=en |автор=Song J. S., etal |год=2007 |тип=journal}}</ref>.

=== Преимущества MACS ===
MACS сравнивали с другими алгоритмами по поиску пиков, например — ChIPSeq Peak Finder<ref>{{статья |заглавие=Genome-wide mapping of in vivo protein-DNA interactions |издание=Science |номер=316 |страницы=1497—1502 |ссылка=http://science.sciencemag.org/content/316/5830/1497.long |pmid=17540862 |doi=10.1126/science.1141319 |язык=en |тип=journal |автор=Johnson D. S., etal |год=2007 |archivedate=2016-05-30 |archiveurl=https://web.archive.org/web/20160530211131/http://science.sciencemag.org/content/316/5830/1497.long }}</ref>, [http://www.bcgsc.ca/platform/bioinfo/software/findpeaks/ FindPeaks] и [http://mendel.stanford.edu/SidowLab/downloads/quest/ QuEST]. Для того чтобы сравнить специфичность их предсказаний, ChiP-Seq данные были поменяны местами с соответствующими им контролями, а по результатам программной обработки вычислили FDR каждого алгоритма. MACS каждый раз давал меньше ложноположительных результатов, нежели другие методы (рис.2).

== SICER ==
'''SICER''' ({{lang-en|spatial identification of ChIP-enriched regions}}) — это компьютерный метод, разработанный для идентификации кластеров, обогащённых выровненными ридами на основании данных ChiP-Seq эксперимента с учётом оценки уровня шума. Хотя данный метод создан для анализа слабых и диффузных сигналов, характерных для модификаций [[Гистоны|гистонов]], при правильном выборе параметров он может быть применён для поиска узких пиков, соответствующих сайтам связывания, например, [[Факторы транскрипции|транскрипционных факторов]].

Программа представляет собой пакет, использующий в качестве интерфейса командную строку, подходит для операционных систем [[UNIX|Unix]]/[[Linux]], язык программирования: [[Python]].

Главной особенностью метода является объединение сигналов от соседних [[Нуклеосома|нуклеосом]], несущих одинаковую модификацию и расположенных близко друг к другу. Данный подход улучшает отношение сигнал-шум и особенно полезен при диффузных модификациях [[хроматин]]а, которые могут иметь большую (килобазы или даже мегабазы) протяжённость в геноме. C помощью математического подхода идентифицируются статистически значимые области генома, обогащённые выровненными ридами. Далее проводится оценка уровня шума с использованием контрольной библиотеки (создаётся на основе контрольных ChIP-Seq экспериментов на тотальной клеточной ДНК без стадии [[Иммунопреципитация|иммунопреципитации]])<ref name=":0">{{статья |заглавие=A clustering approach for identification of enriched domains from histone modification ChIP-Seq data |язык=en |издание=Bioinformatics |том=25 |номер=15 |страницы=1952—1958 |doi=10.1093/bioinformatics/btp340 |тип=journal |автор=Chongzhi Zang, Dustin E. Schones, Chen Zeng, Kairong Cui, Keji Zhao |число=1 |месяц=8 |год=2009}}</ref>.
[[Файл:Острова SICER.jpg|мини|'''Рисунок 3.''' ''Схематическая иллюстрация определения островов. По оси х откладываются координаты генома, где каждый интервал представляет собой «окно». По оси y откладываются покрытия ридами: каждая чёрная вертикальная линия представляет собой степень обогащения ридами относительно фона конкретного пика. Участки, подчеркнутые фиолетовыми горизонтальными линиями, — идентифицированные «острова». Внутри первого острова есть локальные отсутствия покрытия, которые называются гэпы (отмечены жёлтым)'']]
На основании полученных данных идентифицируется набор локальных максимумов покрытия ридами (пиков), для которых отношения покрытия ридами превышает фоновый шум в определённое пороговое число раз. Геном делится на неперекрывающиеся участки фиксированной длины — «окна», статистически значимые области генома называются «острова». Статистически значимыми областями генома принято считать кластеры окон, для которых отношение их покрытия ридами к фоновому покрытию превышает некий заданный порог. Окна внутри островов, для которых велик фоновый шум или маленькое покрытие ридами — «[[гэп]]ы» (рис. 3).

Острова представляют особый интерес. Гэпы разрешены в островах для учёта:
* Неточности ридов и неоднозначности данных [[Секвенирование|секвенирования]]
* Повторяющихся геномных участков и дуплицированных ридов
* Немодифицированных нуклеосом

Размер гэпа может быть скорректирован в зависимости от природы модификации хроматина. [[Статистический вес]] острова высчитывается, основываясь на степени обогащения ридами на всем идентифицированном участке, а не на степени обогащения ридами каждого отдельного пика<ref>{{Статья |автор=Joel Rozowsky, Ghia Euskirchen, Raymond K. Auerbach, Zhengdong D. Zhang, Theodore Gibson |заглавие=PeakSeq enables systematic scoring of ChIP-seq experiments relative to controls |ссылка=http://www.nature.com/nbt/journal/v27/n1/full/nbt.1518.html |язык=en |издание=[[Nature Biotechnology]] |год=2009-01-01 |том=27 |выпуск=1 |страницы=66—75 |issn=1087-0156 |doi=10.1038/nbt.1518 |издательство=[[Nature Publishing Group]] |archivedate=2016-05-11 |archiveurl=https://web.archive.org/web/20160511214059/http://www.nature.com/nbt/journal/v27/n1/full/nbt.1518.html }}</ref>.
[[Файл:Вес островов и гэпы.jpg|мини|'''Рисунок 4.''' ''Суммарный статистический вес островов как функция от размера гэпа (измеряется в окнах) при фиксированном размере окна (в данном случае — 200 пн)'']]
На вход программе подаётся библиотека выровненных ридов из данных ChIP-Seq эксперимента, выбирается размер окна, эффективная длина генома (наибольшая длина генома, покрытая выровненными ридами) и размер гэпа. Размер гэпа должен быть кратным выбранному размеру окна. Как правильно, чем большую протяжённость имеют домены, тем большого размера следует выбирать гэп. Для точной оценки размера гэпа можно представить суммарный статистический вес всех значимых островов как функцию от размера гэпа. Таким образом, лучшим размером гэпа будет тот, который соответствует наибольшему весу островов. На рис. 4 рассмотрен пример такого графика, в данном случае размер гэпа следуют выбирать равным шести окнам.
[[Файл:Схема SICER.jpg|мини|'''Рисунок 5.''' ''Общая схема алгоритма SICER'']]
При наличии контрольной библиотеки идентифицированный набор островов проверяется на способность кластеризоваться с мягким порогом Е-value (основан на предварительном анализе обогащения ридами относительно фона в контрольной библиотеке). Далее SICER вычисляет p-value и q-value для каждого идентифицированного острова с использованием [[Распределение Пуассона|распределения Пуассона]] и FDR уже с учётом удаления артефактов. При отсутствии контрольной библиотеки тоже можно получить статистически значимые острова, используя различные пороги E-value. На выходе программы получается набор островов, высоко обогащённых выровненными ридами<ref>{{Статья|автор=Shiliyang Xu, Sean Grullon, Kai Ge, Weiqun Peng|заглавие=Spatial Clustering for Identification of ChIP-Enriched Regions (SICER) to Map Regions of Histone Methylation Patterns in Embryonic Stem Cells|ссылка=https://dx.doi.org/10.1007/978-1-4939-0512-6_5|язык=en|ответственный=Benjamin L. Kidder|издание=Stem Cell Transcriptional Networks|издательство=Springer New York|год=2014-01-01|страницы=97—111|isbn=9781493905119|doi=10.1007/978-1-4939-0512-6_5}}</ref> (рис. 5).

Особенность данной программы состоит в том, что для её алгоритма локальная степень обогащения сильно зависима от контекста. Иными словами, если окно является обогащённым, этого ещё недостаточно, чтобы оно стало статистически значимым. Однако в присутствии соседних обогащённых окон оно (наряду с другими членами кластера) становится таковым. Это и отличает SICER от подходов, основанных на локальной статистике, которые бы идентифицировали бы его как не значимый в обоих случаях. Данный метод, благодаря нивелированию общего шума, подходит для таких задач, как статистический анализ данных ChIP-Seq экспериментов, нормализация данных и масштабный анализ больших последовательностей с диффузным обогащением<ref name=":0" />.

'''Математика'''

1. '''Схема подсчёта статистических весов'''

Геном эффективной длины <math>L
</math> делится на неперекрывающиеся окна размером <math>w
</math>. Статистический вес <math>s
</math> для окна с <math>l
</math> ридами идентифицируется как <math>s(l)=-\log{P(l,\lambda)}
</math>, где <math>P(l,\lambda)
</math> — распределение Пуассона с параметром <math>\lambda=wN/L
</math>'','' учитывающим среднее число ридов в окне; <math>N
</math> ''-'' общее число ридов в библиотеке ChIP-Seq. Статистические веса окон представляют собой отрицательный логарифм вероятности нахождения <math>l
</math> ридов в окне при условии, что риды могут охватывать любой участок в геноме с равной вероятностью (модель случайного фона). Статистические веса кластеров окон являются аддитивными, представляющими собой отрицательный логарифм совместной вероятности нахождения наблюдаемой конфигурации в модели случайного фона. Чем выше статистический вес, тем меньше вероятность того, что наблюдаемый профиль является случайным.

2. '''Идентификация острова'''

Для числа ридов в каждом окне устанавливается порог <math>l_0
</math>''<sub>,</sub>'' полученный из значения P-value, основанного на распределении Пуассона:

<math>\sum_{l=l_0}^\infty{P(l,\lambda)\leq{p_0}} 
</math>

Более того, <math>l_0
</math> зависит от размера ChIP-Seq библиотеки.

3. '''Рекуррентное соотношение для вероятности острова с данным статистическим весом на случайном фоне'''

Сначала высчитывается вероятность распределения статистических весов отдельного окна:

<math>p(s)=\sum_{l\geq{l_0}}{\delta(s-s(l))P(b,\lambda)}
</math>

где <math>\delta()
</math> — дельта-функция Дирака. Вероятность <math>t
</math> того, что окно является «неподходящим» равна:

<math>t=P(0,\lambda)+P(1,\lambda)+...+P(l_0-1,\lambda)
</math>

Далее вводится параметр <math>g
</math>'','' отражающий число гэпов в окне. Острова идентифицируются как кластеры окон, разделённые гэпами, не большего размера, чем <math>g
</math>. На основании этого высчитывается гэп-фактор <math>G
</math>. Статистический вес острова является совокупным статистическим весом всех окон на данном острове.

<math>G=1+t+t^2+...+t^g</math>

<math>M(s)</math> — вероятность нахождения острова со статистическим весом <math>s</math>, начиная от заданной позиции вдоль генома; зависит от <math>\lambda</math>, <math>l_0</math> и <math>g</math> через <math>p(s)</math>'','' <math>t</math> и <math>G</math>''.'' Так как остров должен быть окружён гэпами размеров не меньше, чем <math>g+1</math>, <math>M(s)</math> можно записать как сумму:

<math>M(s)=t^{g+1}\Mu(s)t^{g+1}</math>

Тогда [[Рекуррентная формула|рекуррентное соотношение]]:

<math>\Mu(s)=G(\lambda,l_0,g)\int\limits_{s_0}^{s} ds'\Mu(s-s')p(s')</math>

Граничные условия: <math>\Mu(0)=G(\lambda,l_0,g)^{-1}</math> и <math>s_0=-ln{P(l_0,\lambda)}</math>.

Число островов со статистическим весом: <math>s-LM(s)</math>.

4. '''Асимптотика для распределения статистических весов островов на случайном фоне'''

Подставим анзац в рекуррентное соотношение:

<math>\Mu(s)=\alpha{exp(-\beta{s})}</math>

Тогда

<math>G(\lambda,l_0,g)\sum_{l\geq{l_0}}P(l,\lambda)^{1-\beta}=1</math>

5. '''Нахождение значимых островов без использования контрольной библиотеки'''

Распределение статистических весов островов на случайном фоне позволяет выявить пороговое значение статистического веса <math>s_T</math>''<sub>,.</sub>'' Данный параметр используется для поиска статистически значимых островов. <math>s_T</math> определяется c учётом того, что ожидаемое число островов со статистическим весом выше, чем порог <math>s_T</math>''<sub>,</sub>'' должно быть не больше порогового значения ''E-''value <math>e</math>'':''

<math>\sum_{s\geq{s_T}}LM(s)\leq{e}</math>

6.'''Нахождение значимых островов c использованием контрольной библиотеки'''

Сначала используется мягкий порог ''E-''value для идентификации набора островов-кандидатов, которые демонстрируют кластеризацию при модели случайного фона. Затем для каждого кандидата высчитывается число картированных ридов <math>n_s</math> и ридов контрольной библиотеки <math>n_c</math>. ''P''-value вычисляется как:

<math>\textstyle \sum_{n=n_s}^\infty \displaystyle P(n_s,cn_c)</math>

где ''с'' равно соотношению размера библиотеки ChIP к размеру контрольной библиотеки (<math>c=N_s/N_c</math>). Острова-кандидаты с <math>n_s\leq{cn_c}</math> отбрасываются, потому что необходимо рассматривать только случаи обогащения. Значимые острова могут быть идентифицированы с помощью порога ''P''-value с использованием [[Коррекция на множественное тестирование|поправки Бонферрони]]<ref name=":0" />.

== Примечания ==
{{примечания}}

[[Категория:Биоинформатика]]
[[Категория:Методы молекулярной биологии]]