Распознавание генов Факультет биоинженерии и биоинформатики МГУ, второй набор, второй курс. Апрель 2005 М. Гельфанд (лекции) А. Неверов (задания) Е. Ермакова, Р.Нуртдинов (занятия) А. Казаков (примеры)
Распознавание генов Поиск открытых рамок считывания Использование статистики (отличия белок-кодирующих и некодирующих областей) Идентификация начал генов – участки связывания рибосом (прокариоты) Экзон-интронная структура (эукариоты) Сравнения с известными генами Геномные сравнения
Вероятность открытой рамки считывания длины не меньше данной
ORFы в геноме K. pneumoniae
Перепред- сказание (E. coli)
Сильное перепред- сказание (Aeropyrum pernix)
Поиск открытых рамок в заданной последова- тельности
Генетический код
Codon usage (статистика употребления кодонов) частоты кодонов отличаются от частот триплетов в некодирующих областях –различия в частотах аминокислот в белках –различия в частотах синонимичных кодонов частоты синонимичных кодонов –специфичны для генома –коррелируют с концентрациями тРНК
Ещё про codon usage различается у высоко- и низко- экспрессируемых генов (у высокоэкспрессируемых генов больше доля «оптимальных» кодонов) – прокариоты, дрожжи нестандартный у горизонтально перенесенных генов у фага T4 – близок к хозяйскому (E. coli) у ранних генов, специфический (соответствует своему набору тРНК) – у поздних
Кодирующий потенциал Функция, измеряющая, насколько участок генома похож на белок-кодирующий (и отличается от некодирующего) с точки зрения статистики Можно вычислять кодирующий потенциал –скользящего окна (не слишком маленького!) –открытой рамки считывания Нужна обучающая выборка генов (и межгенных промежутков) из данного организма
E. coli. Окно 48 нт
E. coli. Окно 96 нт
Сравнение предсказаний при разной длине окон
Gene- Mark
Сигналы на границах генов dnaN ACATTATCCGTTAGGAGGATAAAAATG gyrA GTGATACTTCAGGGAGGTTTTTTAATG serS TCAATAAAAAAAGGAGTGTTTCGCATG bofA CAAGCGAAGGAGATGAGAAGATTCATG csfB GCTAACTGTACGGAGGTGGAGAAGATG xpaC ATAGACACAGGAGTCGATTATCTCATG metS ACATTCTGATTAGGAGGTTTCAAGATG gcaD AAAAGGGATATTGGAGGCCAATAAATG spoVC TATGTGACTAAGGGAGGATTCGCCATG ftsH GCTTACTGTGGGAGGAGGTAAGGAATG pabB AAAGAAAATAGAGGAATGATACAAATG rplJ CAAGAATCTACAGGAGGTGTAACCATG tufA AAAGCTCTTAAGGAGGATTTTAGAATG rpsJ TGTAGGCGAAAAGGAGGGAAAATAATG rpoA CGTTTTGAAGGAGGGTTTTAAGTAATG rplM AGATCATTTAGGAGGGGAAATTCAATG
… после выравнивания dnaN ACATTATCCGTTAGGAGGATAAAAATG gyrA GTGATACTTCAGGGAGGTTTTTTAATG serS TCAATAAAAAAAGGAGTGTTTCGCATG bofA CAAGCGAAGGAGATGAGAAGATTCATG csfB GCTAACTGTACGGAGGTGGAGAAGATG xpaC ATAGACACAGGAGTCGATTATCTCATG metS ACATTCTGATTAGGAGGTTTCAAGATG gcaD AAAAGGGATATTGGAGGCCAATAAATG spoVC TATGTGACTAAGGGAGGATTCGCCATG ftsH GCTTACTGTGGGAGGAGGTAAGGAATG pabB AAAGAAAATAGAGGAATGATACAAATG rplJ CAAGAATCTACAGGAGGTGTAACCATG tufA AAAGCTCTTAAGGAGGATTTTAGAATG rpsJ TGTAGGCGAAAAGGAGGGAAAATAATG rpoA CGTTTTGAAGGAGGGTTTTAAGTAATG rplM AGATCATTTAGGAGGGGAAATTCAATG cons. tacataaaggaggtttaaaaat num
Участки связывания рибосом
rbsD в E. coli Eco AGGATTAAACTGTGGGTCAGCGAAACGTTTCGCTGATGGAGAAAAAAATGAAAAAAGGC Eco ACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTG
rbsD в энтеробактериях Sty AGGGTTACACTGCGGC-CAGCGAAACGTTTCGCTAGTGGAGCAGAAAAATGAAGAAAGGC Sen AGGGTTACACTGCGGC-CAGCGAAACGTTTCGCTAGTGGAGCAGAAAAATGAAGAAAGGC Stm GGGGTTACACTGCGGC-CAGCGAAACGTTTCGCTAGTGGAGCAGAAAAATGAAGAAAGGC Eco AGGATTAAACTGTGGGTCAGCGAAACGTTTCGCTGATGGAGAA-AAAAATGAAAAAAGGC Ype TTTTCTAAACTCCTTGTTAGCGAAACGTTTCGCTCTTGGAGTA-GATCATGAAAAAAGGT ** *** **************** ***** * * ***** ***** Sty ACCGTACTCAACTCTGAAATCTCGTCGGTCATTTCCCGTCTGGGGCATACTGATACTCTG Sen ACCGTACTCAACTCTGAAATCTCGTCGGTCATTTCCCGTCTGGGGCATACTGATACTCTG Stm ACCGTACTCAACTCTGAAATCTCGTCGGTCATTTCCCGTCTGGGGCATACTGATACTCTG Eco ACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTG Ype GTATTACTGAACGCTGATATTTCCGCGGTTATCTCCCGTCTGGGCCATACCGATCAGATT * ** ** **** ** ** **** ** *********** ***** *** *
rbsD в энтеробактериях: ответ Sty AGGGTTACACTGCGGC-CAGCGAAACGTTTCGCTAGTGGAGCAGAAAAATGAAGAAAGGC Sen AGGGTTACACTGCGGC-CAGCGAAACGTTTCGCTAGTGGAGCAGAAAAATGAAGAAAGGC Stm GGGGTTACACTGCGGC-CAGCGAAACGTTTCGCTAGTGGAGCAGAAAAATGAAGAAAGGC Eco AGGATTAAACTGTGGGTCAGCGAAACGTTTCGCTGATGGAGAA-AAAAATGAAAAAAGGC Ype TTTTCTAAACTCCTTGTTAGCGAAACGTTTCGCTCTTGGAGTA-GATCATGAAAAAAGGT ** *** **************** ***** * * ***** ***** Sty ACCGTACTCAACTCTGAAATCTCGTCGGTCATTTCCCGTCTGGGGCATACTGATACTCTG Sen ACCGTACTCAACTCTGAAATCTCGTCGGTCATTTCCCGTCTGGGGCATACTGATACTCTG Stm ACCGTACTCAACTCTGAAATCTCGTCGGTCATTTCCCGTCTGGGGCATACTGATACTCTG Eco ACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTG Ype GTATTACTGAACGCTGATATTTCCGCGGTTATCTCCCGTCTGGGCCATACCGATCAGATT * ** ** **** ** ** **** ** *********** ***** *** *
Паттерн нуклеотидных замен в белок-кодирующих областях: pdxB в энтеробактериях Sty TCGCTCG--CAGCGGAAAGAGGATTACGCCCTTCGCCTGGAGGCTGTGCAGGGGC---GCCGGAGATGGGATGCATAATT Stm TCGCTCG--CAGCGGAAAGAGGATTACGCCCTTCGCCTGGAGGCTGTGCAGGGGC---GCCGGAGATGGGATGCATAATT Sen TCGCTCG--CAGCGGAAAGAGGATTACGCCCTTCGCCTGGAGGCTGTGCAGGGGC---GCCGGAGATGGGATGCATAATT Eco TTGCCCG--TGCCAGACGGCAGATTATCTCCCTGACCTGGTGGTTGCCCAGGAGGAGGGCCGGAAATAGGTTGTATCATT Kpn ----CGG--TGGCGCAGTGCCTGATGGG-CCTCGCCCTGGAGGACGGTCTGGCAT---ATCAGCAAGGGGGTGCGTCATG Ype TTGTTAGAACAGGGGAAAACGGTAAACAGTGTGGCATTAGATGTCGGTTATAGCT-----CCGCCTCTGCTTTTATCGCC * * * * * * * * * * * Sty AATTATCCTTTAAC CATAAATCTGAGCAATA-TATGCTTGGCGGCCAGATTATGGC--ACACTTGTCCGG Stm AATTATCCTTTAAC CATAAATCTGAGCAATA-TATGCCTGGCGGCCAGATTATGGC--ACACTTGTCCGG Sen AATTATCCTTTAAC CATAAATCTGAGCAATA-TATGCCTGGCGGCCAGATTATGGC--ACACTTGTCCGG Eco ACGTATCCTTATAC CTGAAATCTTCGCAAG--TATGCCTGGCCGCGAGATTATGGC--ACACTTGTCCGG Kpn ATTCATCCTTTCGATATCGCGGTGCTGGAACCAGGTGATGAGTATGCCTGGCGGCCAGATTATGGC--ACACTTCCCCAG Ype ATGTTTCAGCAAATAT CGGGTACCA-CGCCTGAGCGTTTCCGGCGGGGCAATAGTGGCTTATACTAAGCCCC * ** * * * * *** * ** **** * *** ** Sty TTAACTCTCGTT-CTCAAACAG------GTACGACAGTC--GTGAAAATTCTCGTTGATGAAAATATGCCTTACGCCCGC Stm TTAACTCTCGTT-CTCAAACAG------GTACGACAGTC--GTGAAAATTCTCGTTGATGAAAATATGCCTTACGCCCGC Sen TTAACTCTCGTT-CTCAAACAG------GTACGACAGTC--GTGAAAATTCTCGTTGATGAAAATATGCCTTACGCCCGC Eco TTAACTCTCGT--CTCATACAG------GTAACACAAAC--GTGAAAATCCTTGTTGATGAAAATATGCCTTATGCCCGC Kpn TTAACTCTCGTT-CTCAGACAG------GTACTGAACT---GTGAAAATCCTCGTTGATGAAAATATGCCCTATGCCCGT Ype CTGTTTTTCATCTGTATGGCAGTTCGCTGTCGGAGAGTAAAGTGAAAATTCTGGTTGATGAAAATATGCCGTACGCTGAG * * ** * * *** ** * ******** ** ***************** ** **
Белковое выравнивание (ribD) Eco V_____QDEYYMARALKLAQRGRFTTHPNPNVGCVIVKDGEIVGEGYHQRAGEPHAEVHA QD +M RAL LA +G +TT PNP VGCV VK+GEIVGEG+H +AG+PHAE A Hin MLEFSSQDCVFMQRALDLAAKGQYTTTPNPSVGCVLVKNGEIVGEGFHFKAGQPHAERVA Eco GCGCGCCTGGAGGACTAA----G CCGTGCAGGAC-GAGTATTACATGGCGCGGGCGCTAA Hin GAAAAATTAAAGGATTAATTATGCTTGAATTTTCCTCACAAGATTGCGTATTT-ATGCAACGTGCCTTAG * * **** *** * ** ** ** * ***** *** ** ** **
Множественное выравнивание REC tttttatttcaggcaatcggggtgaat gtggcgcaggcggaagtgttgaat RECO tttttatttcaggcaatcggggtgaat gtggcgcaggcggaagtgttgaat RECS tttttatttcaggcaatcggggtgaat gtggcgcaggcggaagtgttgaat RTY tagcgcctgttttgatttatggtgaacggggttaatgtggcgcaggcggaagtgttgaat RSY tagcgcctgttttgatttatggtgaacggggttaatgtggcgcaggcggaagtgttgaat REO atagcgcctgtttgatttcattgaattggggaaggcgtgtctacggcggaagtattgaat RYPK gccggcctgtgcagatctaatagttgggggaaaagtgtgtcgaccgcagcagtgataaac RYP gccggcctgtgcagatctaatagttgggggaaaagtgtgtcgaccgcagcagtgataaac RYE aaccggcctgtgcagatctcatagttggggaatagtgtgtcaaccgcagcagtgataaat RVFI tattattgatgagttttttatgtccagcatgatcgcagagcaaccaatggaa REC06584 f l f q a i g v n = = = V A Q A E V L N RECO04717 f l f q a i g v n = = = V A Q A E V L N RECS04752 f l f q a i g v n = = = V A Q A E V L N RTY01088 * r l f * f m v n g v n V A Q A E V L N RSY05814 * r l f * f m v n g v n V A Q A E V L N REO01497 i a p v * f h * i g e g V S T A E V L N RYPK00397 a g l c r s n s w g k s V S T A A V I N RYP04048 a g l c r s n s w g k s V S T A A V I N RYE04903 n r p v q i s * l g n s V S T A A V I N RVFI i i d e f f m s s M I A E Q P M E
Распознавание генов в отсутствие обучающей выборки «псевдообучающая выборка»: протяженные рамки считывания гены, предсказанные по сходству
Репликация и статистика ДНК GC-сдвиг (G-C)/(G+C) Направление транскрипции DnaA сайты
Эукариоты (человек) В среднем 9-10 экзонов (кодирующих) на ген Средняя длина (внутреннего) экзона нуклеотидов Часто очень длинные интроны
Длины экзонов: человек, нематода C. elegans, дрозофила
Длины интронов
Бета-глобин человека
Хемотрипсин крысы
… ничего … (28S рРНК человека)
Статистические методы Скользящее окно не работает! (~ 1990) Статистика кодирующих и некодирующих областей + сайты сплайсинга – ещё одна вариация на тему динамического программирования
Сайты сплайсинга Donor sites gtgggatgatgtaagtattggggcggcccg tcaaaacaaggtaagaaatgaggtatgcct agctcccaaggtaggaggttgagtgttgtg agtggccaaggtatggtggatggaaattgc tggaaaaagcgtaagtcactctaattttat ctctcaaaaagtaagctttgtgagcatttc atcttcaagggtgagcatgtgtgttatgct tttcagaattgtaagagtacacattttaag gccagaaaaggtcagtactttctttcacac tacctcacaggtatgaattttctagttctt atctttcaaggtagagtatatgaatgttac atgtggattcgtaagtattcaacacattca aaaatatccagtaagcagttctgatgtttg ccaggagccggtgaggggctggtgggctct aatggatgaggtgggtacttagggcttctg atttcaaaaagtaagttttccctggagaaa aatttgtagagtatccttgatttgacgaat cagacaatgggtaagtacatgcttgttccc gtctgttaaggtaggtataccccatcacaa gttcaaaaaggttggtcacatgttcttgat attcggccaggtatgggtagtgtgctgaga acatatgcaggtaaacaacttaactcaaat aaagaaagaggtgagagggtgttttaattt ccagctccaggtaagccatctggaaagagc gtcttaacaggtaaatgccaccctttcccc Acceptor sites gtttcttcttacatttctaggactcaacta ttcacgtttttgccttccaggagacagagc tttcaatatttattacccaggaccccaaat gtgttatttacatttttcaggaatggacaa tttttctgcttctccaacagctatactaaa ttgttgtgttcacttcacagcatatatcgc tccgttgttttatttcccagaatgattcaa tggtttttcattgtttttagtggtgcaaaa tctaacttcatttcctccaggacaaatatc gttttgttggtgttttatagctggccaact acatgtgttctcatttttaggaagtgatag ctgttcttgttctcccttagcccaaagcag atgcctttcatttctattagctggaatctg ctgttattaaaatttgacaggagaagctga ttttttattcctacttccaggggactgctg tttgttgttgcttaactcagaaagaaataa tacttaacatgatggtccagatataacaaa cttgtgtttttgatactcagacctggctat ttgatttattgattttctagattatttcag gtccttaatgtcctttgtaggtggttcttc gcattattctcaccttccaggctatcacta aatatctcttccctatttagatgtcatcga aaggatatttataattttaggctgatcctg ttttatcttttatattacaggttctgtaaa ttcatattcatttgttgcagaagtggaagc
Распознавание сайтов сплайсинга
Список потенциальных экзонов
Граф динамического программирования
Путь = экзон-интронная структура
Gen- Scan
Сравнительные методы BLASTN: ESTs и альтернативный сплайсинг BLASTX BLASTX+статистика Сравнение с известными белками Геномные сравнения –выравнивание ДНК –выравнивание белков All of the above and more…
Семей- ство про- грамм BLAST
ESTs: короткие фрагменты (клонированной) мРНК Характерная длина ~300 нт Ошибки секвенирования Ошибки клонирования –несплайсированный транскрипты –геномная ДНК Обогащение к 3-концу (PolyA-праймеры) Альтернативный сплайсинг: 30-50% генов
Human Genome Browser – поиск по имени гена
Результат
MAGE-C1
Ещё о сравнении предсказаний
Альтернативный сплайсинг генов человека
MAGE- A2
GenomeScan=GenScan+BLASTX
Сплайсированное выравнивание Сравнение (формально транслированной) ДНК с аминокислотной последовательностью родственного белка. Динамическое программирование, дополнительная операция – интрон –Только на потенциальных сайтах сплайсинга –Небольшой штраф –Учёт особенностей экзон-интронной структуры – минимальная длина интрона (зависит от генома)
VISTA (human-dog-mouse)
Сплайсированое выравнивание геномных последовательностей
Другие возможности
Человек- мышь (мульти- генное семейство)
Ткане- специфич- ная экспрессия