Аннотация результатов, полученных в 2014 году
На базе лаборатории в КФУ оптимизированы условия поддержания клеточной линии Pv11 способной к индуцируемому ангидробиозу. Был опробованы ряд сред и сделан вывод, что среда IPL является наиболее подходящей для поддержания культуры клеток и инициации в них ангидробиоза. Осуществлено секвенирование генома клеточной линии Pv11 хирономиды с использованием двух подходов спаренных прочтений: paired ends (PE) и mate pairs (MP). С использованием комбинации специализированного ПО и различных методик олучена более точная и целостная сборка генома клеточной культуры хирономиды (значение N50 увеличено более чем в 2 раза). Установлено, что геном клеточной линии отличается от генома личинок на 1-2% на нуклеотидном уровне. С использованием новых данных mRNASeq было идентифицировано 27 новых кодирующих белки генов не описанных в изначальной версии расшифровки генома хирономиды. На основе новых данных о структуре генома клеточной культуры Pv11 и анализа экспрессии мРНК (методом SE mRNA-Seq) были получены полногеномные профили экспрессии для каждой стадии цикла ангидробиоза. Результаты оказались достаточно неожиданными: детальный анализ экспрессии 132 генов, активность которых непосредственно в личинках связана с процессом ангидробиоза (LEA белки, белки теплового шока, ферменты синтеза трегалозы, антиоксиданты и др) показал, что сходная по динамике и интенсивности экспрессия этих генов в клетках Pv11 инициируется уже на стадии преинкубации с трегалозой еще до начала обезвоживания.Анализ данных экспрессии генов расположенных в ARIds (участки генома содержащие кластеры генов-паралогов, активность которых определяет успешность ангидробиоза, см. Gusev et.al, 2014) показал что лишь часть из этих генов экспрессируется в клетках Pv11 в гидратированном состоянии и повышает экспрессию в ответ на обезвоживание. Таким образом, можно сделать вывод, что одной из характеристик ARIds и генов находящихся в них является ткане-специфичность. И, вероятно, активная паралогизация генов связанных с ангидробиозом связана с их “специализацией” на определенном типе клеток или тканей. Произведена подготовка библиотек для CAGE анализа из общей РНК клеток Pv11 на разных этапах ангидробиоза в трех биологических повторностях. Данные (по 30-45 миллионов прочтений в каждой библиотеке) картированы на геном Pv11 и интегрированы в браузер ZENBU для выделения специфичных для цикла регуляторных элементов, включая энхансерные РНК и промотеры.К моменту отчета выделено 17 типов (многокопийных, высоко-консервативных, 500-800 п.н.) некодирующих РНК активность которых сопряжена со стадиями ангидробиоза. Предварительный анализ промотерной активности генов в ARIds участках генома свидетельствует о дифференциальных сайтах инициации транскрипции в зависимости от стадии цикла ангидробиоза.

Аннотация результатов, полученных в 2015 году
На 2015 год были запланированы ключевые эксперименты по профилированию экспрессионной активности регуляторных элементов в тканях хирономид и клеточной культуре Pv11 полученной из эмбриональной клеточной массы хирономид. На конец 2015 были выполнены все работы по проекту, а также проведен ряд дополнительных экспериментов, что позволило создать более целостную картину процессов связанных с ангидробиозом в клетках и тканях хирономиды. Были подобраны наиболее перспективные в качестве модельных ткани хирономиды для полногеномного и протеомного анализа. Это – жировое тело, кишечник, нервная трубка. Были получены полногеномные профили экспрессии мРНК в выбранных органах хирономиды на активной стадии и после 24 часов медленного обезвоживания (это точка, когда генетические процессы связанные с ангидробиозом полностью активирутся, но в то же время в теле личинок присутствует еще значительное количество воды, так что извлечение органов не сопряжено с повреждением последних). Показано, что значительная часть из 17154 генов в геноме хирономиды является тканеспецифичными. Среди тканей, наиболее большое количество дифференциально экспрессирующихся генов (P<0.05 среднее в трех повторностях) в ответ на обезвоживание было обнаружено в нервной системе (780 различных транскриптов). Анализ экспрессии 27 генов кодирующих Lea белки в хирономиде показал, что часть из них ( в том числе PvLea4, PvLea3, PvLea6, PvLea11) экспрессируются на высоком уровне во всех исследованных тканях, и, в то же время их экспрессия дополнительно индуцируется обезвоживанием. В то же время, экспрессия ряда других PvLea генов носила более тканеспецифичный характер. Схожие паттерны наблюдались и в группах генов кодирующих антиоксиданты, в частности тиоредоксины. Был продолжен анализ экспрессионной активности (с учетом новых данных по методике CAGE) в клеточной культуре Pv11. Основная массаа дифференциально экспрессирующихся CAGE кластеров приходилось на период преинкубации клеток с трегалозой, и регидратации после полностью безводного состояния, в то время как период непосредственно процесса обезвоживания характеризовался значительно меньшим изменением экспрессии генов. Мы отметили факт того что бОльшая часть генов показывающая дифференциальную экспрессию хотя бы на одной стадии ангидробиоза характеризуется двумя пиками сайтов инициации транскрипции и на разных стадиях ангидробиоза “соотношение” интенсивности пиков значительно изменяется. Это дает основание предположить существование особых транскрипционных факторов, активность которых связана с изменением концентрации трегалозы, либо c процессом обезвоживания. Кроме того мы выделили ряд некодирующих РНК экспрессирующихся на мультикопийных участках генома хирономиы и обладающих признаками регуляторных (в том числе и энхансерных РНК), - экспрессия которых прямо или обратно коррелирует с активностью генов известных участников ангидробиза, а также экспрессирующихся в участках обладающими свойствами энхансерных зон в геноме хирономиды. Эти РНК будут являться кандидатными целями для in vitro экспериментов. Cравнительный анализ размера геномов хирономид и структуры хромосом P.vanderplanki из различных популяций показал признаки значительных изменений в структуре генома (а также появления ряда новых микро-морфологических признаков) в разных популяциях хирономид. Одним из возможных катализаторов этого процесса может являться процесс повреждения ДНК при ангидробиозе, в том числе генерирующих полиморфные участки в регуляторных областях генома. Был получен профиль распространения белков в клеточной культуры Pv11 на разных этапах цикла ангидробиоза: лабораторный контроль, клетки после 48 часов инкубации с трегалозой, полностью обезвоженные клетки и клетки через 24 часа после регидратации (только живые клетки). В результате был идентифицирован 19001 пептид предствляющий 408 белков, динамика распространенности которых были проанализированы во четырех группах образцов. По результатам соотнесения белкового и РНК профилей наблюдалась значительная корреляция высокоэкспрессирующихся РНК и их белковых продуктов. Были обнаружено присутствие высоких концентрация PvLEA4 и PvLea11 белков, а также ряда тиоредоксинов и других антиоксидантов, что предоставило уникальную информацию о фактическом распространении белковых компонентов в клетках с формирующимся “молекулярным щитом”. На основании полученных данных готовится рукопись статьи о динамике изменения транскрипционного, метаболомного и белкового профилей в клеточной культуре Pv11 в цикле ангидробиоза В соответствии с условиями РНФ проведены два международных семинара с элементами молодежной школы (июль 2015, Москва : http://kpfu.ru/biology-medicine/struktura-instituta/openlab/extreme-biology-lab/novosti/seminar-s-elementami-shkoly-39funkcionalnaya_137652.html и август 2015, Иркутск: http://kpfu.ru/news/seminar-s-elementami-shkoly-39funkcionalnaya_141092.html) ориентированные на полногеномный анализ в экстремофильных организмах, в частности с применением различных методов анализа транскрипционной активности генома. По ключевым результатам опубликованы две статьи и еще четыре рукописи рассматриваются редакциями профильных журналов для публикации.

Аннотация результатов, полученных в 2016 году
На данном этапе проекта, нам удалось показать, что, как и предполагалось по данным биоинформатического анализа, фактор теплового шока является одним из ключевых регуляторов активации транскрипции в ответ на обезвоживание в клетках хирономиды. Результаты направленного подавления активности транскрипционного фактора теплового шока (HSFT) показали значительное снижение выживаемости и общую низкую пролиферацию клеток Pv11 с подавленной экспрессией HSTF после обезвоживания, что также может быть объяснено тем что большинство клеток обезвоженных в условиях сниженной экспрессии HSFT погибали в первый же день регидратации. Также, используя РНК извлеченную из клеток после воздействия siRNA мы провели оценку динамики экспрессии набора из 29 генов, активирующихся на ранних стадиях ангидробиоза (уже на стадии инкубации с трегалозой), сохраняющих высокую экспрессию на всем протяжении ангидробиоза (“ангидробиозные гены”). Для максимальной точности, выбор в оценке экспрессии был сделан в пользу ПЦР в реальном времени в трех независимых биологических повторностях. Изучение экспрессии в клетках инкубирующихся в присутствии трегалозы и подверженных воздействию siРНК, показало что 18 генов из выбранного набора, значительно изменяли экспрессию (t-test, BH adjusted p < 0.05) в клетках с направленной сниженной экспрессией HSFT. Картирование данных CAGE на обновленную сборку генома позволило выявить новые закономерности транскрипционной активности в тканях хирономид в цикле ангидробиоза. Для каждой пары тканей появилась возможность провести анализ обогащенных GO тегов, среди которых 10 являются общими для всех тканей: cell redox homeostasis, protein disulfide oxidoreductase activity, protein-L-isoaspartate (D-aspartate) O-methyltransferase activity, ATPase regulator activity, protein homodimerization activity и другие. Некоторые теги из этой группы соответствуют определенным “Ангидробиотическим генам” – генам напрямую связанным с ангидробиозом - таким как метилтрансферазы PIMT или оксидоредуктазы Trx. Такой результат подчеркивает универсальную роль этих ферментов при переходе к состоянию ангидробиоза. Большая часть “ангидробиозных генов” активировалась при высыхании тканей (75%), при этом ~39% из них проявляли активность в разных тканях схожим образом. С другой стороны, были обнаружены некоторые особенности, связанные с ткане-специфичной активацией этой группы генов: для жирового тела характерны глобины , в кишечнике только один ген тиоредоксина, в мозге это набор из оксидоредуктаз, метилтрансферазы, LEA и др. Среди 31 “ангидробиозного гена”, не повышающих экспрессию в экспериментеприсутствуют стандартные ортологи, либо гены активные в культуре клеток Pv11, но не в тканях, как, например, ген трегалазы. При общем рассмотрении экспрессии “ангидробиозных генов” со всеми остальными активными генами виден общий тренд их активации в экспериментах с обезвоживанием, тогда как все прочие гены в большинстве инактивируются. Кроме того, “ангидробиозные гены” имеют относительно высокую активность и в нормальных условиях. Вместе с достаточно высокой корреляцией между образцами сухих и мокрых тканей мозга и жирового тела, можно предположить, что эти ткани независимо от других частей личинки готовы к переходу в состояние ангидробиоза в любое время (что подкрепляет данные полученные с использованием изолированных нервных тканей личинки). Частичное секвенирование генома и транскриптома хирономиди из других африканских популяций (в Нигерии и Малави) позволил выявить значительную изменчивость в нуклеотидных последовательностях генов, особенно в таковых кодирующих Lea белки. Так, например, мы обнаружили присутствие в геномных последовательностях хирономид из региона Zaria гетерозиготный вариант гена кодирующего белок LEA-5, который нес в себе дополнительную “вставку” кодирующую экстра- канонический повтор из 12 аминокислот, что может быть признаком процесса микроэволюции генов этой группы. В свою очередь, несмотря на успешную идентификацию генов-ортологов в геномах малавийских популяций хирономид, мы обнаружили значительные отличия в первичных последовательностях генов. Так, в среднем на уровне нуклеотидных последовательностей между кодирующими участками нигерийских и малавийских хирономид наблюдалось сходство на уровне 90-95%. В то же время, последовательности кодирующие LEA белки отличались значительность вариабельностью и здесь уровень сходства между ортологами не превышал 85%.В то же время, ген-ортологи LEA отличались схожей динамикой экспрессии, т.е. наиболее представленными транскриптами в обезвоженных личинках как в Малави так и в Нигерии (включая лабораторную линию хирономид), являлись Lea4, Lea1, Lea3, Lea5, Lea7. Полученные данные могут свидетельствовать о значительном влиянии ангидробиоза на процессы радиации в популяциях хирономид. Стоит отметить, что совокупность морфологических и генетических признаков хирономид позволила нам сделать вывод что малавийские популяции можно с достаточной степенью уверенность выделить в отдельный вид.