Ко-эволюция транскриптома и транслятома у млекопитающих
Комментарий к статье: Wang, ZY., Leushkin, E., Liechti, A. et al. Transcriptome and translatome co-evolution in mammals. Nature (2020). https://doi.org/10.1038/s41586-020-2899-z
![]() ![]() © Александр Финошин, аспирант лаборатории биохимии процессов онтогенеза. Область научных интересов: клеточная пластичность, протеомика, гомеостаз клеточного протеома, регенерация, металлопротеины. © Ким Адамейко, младший научный сотрудник лаборатории биохимии процессов онтогенеза. Научные интересы: биоинформатика, транскриптомика, хроматин, беспозвоночные. Junior research scientist Laboratory of Biochemistry. Interests: bioinformatics, transcriptomics, chromatin, invertebrates. Благодаря новым методам на базе высокопроизводительного секвенирования и достижениям масс-спектрометрии мы узнаём всё больше о тонкой взаимосвязи различных уровней экспрессии генетической информации. Хорошо известно, что транскриптом клеток (совокупность мРНК) не определяет автоматически состав клеточного протеома (совокупности белков). Дифференциальная трансляции мРНК на рибосомах вносит существенный вклад в формирование протеома. Поэтому транслятом (совокупность мРНК, транслируемых на рибосомах в данный момент), является важным промежуточным этапом в реализации генетической информации. В статье Wang и соавторов, опубликованной в журнале Nature, впервые проведен сравнительный анализ транскриптомов, протеомов и транслятомов в клетках млекопитающих. Транслятом количественно характеризуют с помощью метода рибосомного профайлинга (Ribo-seq), основанного на изоляции защищённых рибосомами фрагментов мРНК от остальных видов клеточной РНК с последующим их секвенированием. Были получены данные RNA-seq и Ribo-seq для трёх органов (мозг, печень, семенники), формирующихся из разных зародышевых листков, пяти видов млекопитающих (человек, макака, мышь, опоссум и утконос) и курицы. Протеомные данные для человека и мыши взяли из предыдущего исследования. Анализ трёх уровней экспрессии генов у разных организмов позволил сделать интересные эволюционные выводы. Транскрипция не всегда напрямую сопряжена с трансляцией [1]; анализ ранговой корреляции экспрессии на трёх уровнях (транскриптом, транслятом и протеом) подтвердил более тесную связь между «соседними» уровнями, чем между «крайними» - транскриптомом и протеомом (рис. 1A), что согласуется с данными о том, что эффективность синтеза белка является лучшим показателем содержания белка, чем мРНК [2]. При этом наиболее высокая корреляция наблюдалась между транскриптомным и транслятомным уровнями (0.78-0.90). Дисперсия экспрессии возрастала на 12-32% в транслятоме по сравнению с транскриптомом у всех организмов в соматических органах и падала в репродуктивном (рис. 1В), что говорит о тканеспецифичной трансляционной эффективности (ТЭ). Это наблюдается и на внутритканевом уровне. При анализе ТЭ в четырех типах клеток в семенниках мыши авторы установили, что кластер генов, избегающих общей трансляционной репрессии в сперматоцитах, действительно обогащен генами, связанными со сперматогенезом. Исследуя ко-эволюцию экспрессионных уровней, авторы построили филогенетические деревья на основе экспрессии в транскриптоме и транслятоме клеток трёх тканей для примерно 5000 стабильно экспрессирующихся ортологов. В качестве попарного межвидового расстояния взята дисперсия изменений экспрессии (у близких видов она невелика), отнесённая к общей дисперсии экспрессии по всем видам. Таким образом, длины ветвей в деревьях отражают долю изменений экспрессии, связанную с эволюционными процессами. Интересно, что ветви транслятомных деревьев оказались на 20-22% короче, чем у транскриптомных (рис. 1С), что говорит о большей консервативности транслятомного уровня регуляции.
Рисунок 1. А – Ранговая корреляция Спирмана (Spearman’s p) между тремя уровнями экспрессии, рассчитанная по 9642 генам, стабильно экспрессирующихся в человеческих семенниках, мозге и печени. В – Дисперсия экспрессии на транкриптомном и транслятомном уровнях у шести видов в трех органах из разных зародышевых листков. C – Деревья, в которых длина ветвей соответствует доле дисперсии экспрессии, обусловленной эволюционными процессами. Для интерпретации результатов авторы вводят понятие буферизации (“buffering”), или сглаживания изменений, происходящих на разных уровнях (примером могут служить нетранслирующиеся сплайс-изоформы, находящиеся под меньшим давлением отбора и вносящим вклад в большее разнообразие транскриптома). Они также исследуют скорости ко-эволюции экспрессии индивидуальных генов Δ, выводимые из межвидовых дисперсий. Генов, предположительно находящихся под стабилизирующим отбором (Δ<0, FDR < 0.1), оказалось в несколько раз больше, чем тех, которые находятся под положительным (Δ>0). В целом бо́льшая скорость эволюции в транскриптоме может ещё раз подтверждать бо́льшую гибкость регуляции этого уровня экспрессии генов. Подробный анализ некоторых характеристик исследуемых ортологов показал, что гены, более чувствительные к мутациям, проявляют более низкую дивергенцию экспрессии на обоих уровнях вместе с более сильной буферизацией по сравнению с генами, которые относительно толерантны к мутациям (рис. 2). Эффект «важности» (essentiality) выше в соматических тканях по сравнению с семенниками и особенно силен в мозге. Для housekeeping генов (группа Broadly expressed) наблюдались очень низкие и в целом сходные скорости эволюции экспрессии. В отличие от housekeeping генов, тканеспецифичные гены эволюционируют с разной скоростью в разных органах. Гены, специфичные для мозга эволюционируют значительно медленнее, чем гены, специфичные для печени или семенников. Также установлено, что более эволюционно древние гены имеют более низкую дивергенцию экспрессии генов и более сильную буферизацию, чем у эволюционно молодых генов [6-7].
Рисунок 2. Экспрессионная дивергенция на двух экспрессионных уровнях была рассчитана для всех 8,109 стабильно экспрессирующихся ортологов (FPKM>1) для макаки, мыши и опоссума (референс) и для специфических групп генов: Broadly Expressed, Tissue specific, Mutation tolerant/intolerant, Haplo-insufcient/sufcient, old/young. Анализ был ограничен тремя видами, с целью увеличения количества 1:1 ортологов. Соотношения между скоростями дивергенции экспрессии генов (от трансляционной к транскрипционной) показаны справа от каждого набора гистограмм. Планки погрешностей рассчитанные на основе 1000 повторений (bootstrap). Для понимания, в какой степени буферизация посредством активации трансляции могла ослабить эффект наличия лишь одной X хромосомы у самцов и/или эффект мейотического сайленсинга в половых хромосомах, сравнили экспрессионные уровни ортологов на современной Х хромосоме млекопитающих с более древними ортологами на “предковой” половой хромосоме (условная “предковая” половая хромосома была составлена из аутосомных ортологов курицы). Было установлено, что: 1. Эффективность трансляции для генов, сцепленных с Х-хромосомой плацентарных животных, значительно выше, чем для аутосомных генов, особенно в семенниках, тогда как для ортологов Х-сцепленных генов у яйцекладущих (курица и утконос) эффективность трансляции не выше по сравнению с другими аутосомными генами; 2. Более высокие уровни экспрессии X-сцепленных генов на транслятомном и протеомных уровнях по сравнению со транскриптомным уровнем наблюдаются во всех трех органах, с наиболее сильной активацией в семенниках; 3. Повышение уровня экспрессии от транслятома к протеому для X-сцепленных генов является относительно незначительным. Из этого следует, что повышенная экспрессия генов на X-хромосоме в основном обеспечивается транслятомным уровнем (например образованием полирибосом). Авторы статьи предоставляют удобный веб-ресурс для изучения ко-эволюционных характеристик исследованных 5000 ортологов, который доступен по адресу https://ex2plorer.kaessmannlab.org/. Полученные авторами результаты показывают важную роль транслятомного уровня регуляции экспрессии в реализации генетической информации и формировании клеточного протеома в разных органах эволюционно удаленных видов животных. Эволюционный процесс идет за счет более быстрого изменения тканеспецифичных генов и эволюционно молодых генов. При помощи веб-сервиса https://ex2plorer.kaessmannlab.org/ каждый может посмотреть данные по трансляционной эффективности интересующего гена, что, вероятно, поможет более точной интерпретации данных РНК-секвенирования. Результаты, полученные в данной статье, могут стать мощным стимулом в развитии и использовании метода рибосомного профайлинга.
Новость подготовили © Финошин А. © Адамейко К. |