Бесплатный автореферат и диссертация по биологии на тему

Разнообразие репертуаров Т-клеточных рецепторов человека и его изменения в ходе старения

ВАК РФ 03.01.03, Молекулярная биология

Автореферат диссертации по теме "Разнообразие репертуаров Т-клеточных рецепторов человека и его изменения в ходе старения"

Федеральное государственное бюджетное учреждение науки ИНСТИТУТ БИООРГАНИЧЕСКОЙ ХИМИИ имени академиков М.М. Шемякина и Ю.А. Овчинникова Российской академии наук

На правах рукописи

Путинцева Екатерина Викторовна Разнообразие репертуаров Т-клеточных рецепторов человека и его изменения в ходе старения

Специальность - 03.01.03 - молекулярная биология

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата биологических наук

Москва 2014

1 8 СЕН 2014

005552608

Работа выполнена в Лаборатории геномики адаптивного иммунитета в Федеральном государственном бюджетном учреждении науки Институте биоорганической химии им. академиков М.М. Шемякина и Ю.А. Овчинникова (ИБХ РАН).

Научный руководитель:

Чудаков Дмитрий Михайлович, доктор биологических наук. Официальные оппоненты:

Атауллаханов Равшан Иноятович, доктор медицинских наук, профессор, заведующий Лабораторией активации иммунитета (№25) Федерального государственного бюджетного учреждения "ГНЦ институт иммунологии" ФМБА России.

Логачёва Мария Дмитриевна, кандидат биологических наук, ведущий научный сотрудник Лаборатории эволюционной геномики Факультета биоинженерии и биоинформатики Московского государственного университета им. М.В. Ломоносова.

Ведущая организация:

Федеральное государственное бюджетное учреждение науки Институт общей генетики им. Н.И. Вавилова РАН.

Защита диссертации состоится «8» октября 2014 года в 10:00 часов на заседании совета Д 002.019.01 при ИБХ РАН по адресу: 117997, г. Москва, ул. Миклухо-Маклая, д. 16/10.

С диссертацией можно ознакомиться в библиотеке Института биоорганической химии им. академиков М.М. Шемякина и Ю.А. Овчинникова РАН и на сайте www.ibch.ru. Автореферат разослан « » 2014 года.

Ученый секретарь диссертационного совета.

доктор физико-математических наук В.А. Олейников

Общая характеристика работы

Актуальность проблемы

Эффективность и сбалансированность адаптивной иммунной защиты организма от инфекционных и онкологических заболеваний в значительной степени определяется функциональным и рецепторным разнообразием Т-лимфоцитов. Каждый функциональный клон Т-лимфоцитов несет уникальный Т-клеточный рецептор (ТСИ), распознающий чужеродный антиген в составе молекулы главного комплекса гистосовместимости (МНС) и определяющий специфичность адаптивного иммунного ответа. Индивидуальное разнообразие наивного репертуара ТС!* определяет вероятность того, что для новой инфекции или онкологического заболевания найдутся клоны Т-лимфоцитов, специфично распознающие характерные для данного патогена или патологических клеток антигены. После распознавания антигена антиген-специфичные Т-лимфоциты активно размножаются, при этом их концентрация в организме и периферической крови значительно возрастает.

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

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

Настоящая работа посвящена анализу индивидуальных репертуаров ТС1*р человека методом глубокого массированного секвенирования.

1 Например, методов массированного секвенирования репертуаров ТСЯ, иммунофенотипирования Т-клеточных субпопуляций, транскриптомики, в том числе с единичных клеток, статистического анализа репертуаров ТС1?, и их потенциальных эпитопов в составе МНС, а также методов рентгеноструктурного анализа комплексов ТО!/пептид/МНС и масс-спектрометрического анализа репертуаров презентируемых пептидов.

В первой части работы проводится сравнительный анализ состава и структуры индивидуальных репертуаров ТС[*Р родственных и неродственных пар мама-ребенок.

Вторая часть работы посвящена разработке методики нормализации библиотек генов ТСГ! для корректного и воспроизводимого сравнения разнообразия репертуаров.

В третьей части работы с помощью разработанной методики исследуются изменения, происходящие в репертуаре ТСИР человека при старении.

Цель работы

Цель настоящей работы заключалась в изучении индивидуального репертуара 3-цепей Т-клеточных рецепторов периферической крови человека и влияния на него наследственных факторов и старения. Были поставлены следующие задачи:

1. Подготовить для массированного секвенирования библиотеки генов Р-цепей Т-клеточных рецепторов из образцов периферической крови группы матерей и их детей;

2. Провести сравнительный анализ состава и структуры репертуаров генов Р-цепей Т-клеточных рецепторов родственных и неродственных пар мама-ребенок;

3. Разработать метод на основе молекулярного баркодирования, позволяющий проводить корректное сравнение разнообразия репертуаров С01?3 участков библиотек генов ТС1?, полученных от различных доноров, в условиях разного количества клеток на старте, разного качества РНК, и разной глубины массированного секвенирования;

4. Подготовить для массированного секвенирования молекулярно-баркодированные библиотеки генов р-цепей Т-клеточных рецепторов из периферической крови системно здоровых доноров разных возрастов;

5. Провести сравнительный анализ нормированного разнообразия репертуаров генов Р-цепей Т-клеточных рецепторов периферической крови для доноров различного возраста.

Научная новизна и практическая ценность работы

Массированное секвенирование открыло новые возможности для исследования адаптивного иммунитета на уровне индивидуальных репертуаров антител и Т-клеточных рецепторов. В настоящей работе мы применили этот метод для анализа индивидуальных репертуаров СРКЗ ТИф различных когорт системно здоровых доноров. Достигнутая глубина анализа (порядка 0,5-1 миллиона индивидуальных клонотипов для каждого донора), а также применение разработанной нами технологии молекулярного

баркодирования для количественной нормализации сравниваемых библиотек позволило нам на качественно новом уровне охарактеризовать индивидуальные репертуары TCRp.

В данной работе мы изучали процесс формирования индивидуальных репертуаров TCR(5 и влияние наследственных факторов на этот процесс, исследовали структуру индивидуальных репертуаров и степень их пересечения между различными индивидуумами, а также оценивали изменения в разнообразии и общей структуре репертуара TCRp в ходе старения.

Мы уверены, что полученные нами данные интересны как с фундаментальной точки зрения, так и с практической. В частности, мы надеемся, что они будут полезны для разработки вакцин и адаптации их к иммунитету пожилых людей, в планировании трансплантаций, а также в разработке препаратов, поддерживающих иммунитет человека в старости. Кроме того, разработанный нами способ сравнительного анализа репертуаров TCR с помощью метода молекулярного баркодирования позволяет проводить точное и аккуратное сравнение образцов, и мы надеемся, что этот метод будет полезен в дальнейших исследованиях репертуаров TCR и антител.

Структура диссертации

Диссертационная работа изложена на 114 страницах и состоит из введения, обзора литературы, описания материалов и методов исследования, результатов и их обсуждения, выводов, заключения, списка сокращений, а также списка литературы, включающего 153 источника. Диссертация содержит 27 рисунков и 4 таблицы.

Апробация работы

Основные результаты диссертации были доложены на конференции «ZIBI summer symposium», 2013, Берлин, Германия.

Публикации

По результатам работы опубликовано 7 статей в рецензируемых зарубежных журналах.

Содержание работы

1. Сравнительный анализ репертуаров ТСРф матерей и детей

Мы использовали массированное секвенирование для сравнения индивидуальных репертуаров ТС1?р матерей и их детей. Мы стремились выявить специфические характеристики, свойственные репертуарам ТСЯр родственных пар, чтобы определить степень влияния наследственности на формирование индивидуальных репертуаров. Кроме того, мы анализировали влияние отбора в тимусе на формирование индивидуального репертуара ТСРф человека.

Взаимодействие иммунитетов матери и плода активно изучается в контексте осложнений в ходе беременности, аутоиммунных заболеваний, а также гаплоидентичных трансплантаций органов и гематопоэтических стволовых клеток. Достигнутая глубина секвенирования позволила нам исследовать наличие зрелых микрохимерных Т-лимфоцитов матери или плода (то есть Т-клеток, мигрировавших через плаценту в качестве зрелых а/р Т-лимфоцитов и сохранившихся в периферической крови обоих доноров).

В исследовании разнообразия репертуаров ТСЯр принимали участие три мамы и их дети (по два ребенка каждой из мам). Мы получили более 1х107 С01*3-содержащих прочтений секвенирования (ридов) для каждой мамы и около ЗхЮ6 ридов - для каждого ребенка. Из полученных массивов последовательностей с помощью нашего программного обеспечения МГГС1* было идентифицировано 500 ООО - 2 ООО ООО различных клонотипов ТСВР в образцах каждого из доноров. Полученные индивидуальные репертуары ТСРР мы использовали для сравнительного анализа (Рисунок 1).

I

I

Ребе ной А1 Ребенок А2

I

I

I

Ребенок Б1 Ребенок Б2 Мама В РебенокВ1 Ребенок В2

Тип Клона

■ Редкие

■ Малые Средние

■ Большие

■ Очень большие

Размер Клока Менее 0,001% От 0.001% до 0.01% От 0.01% до 0.1% От 0.1% до 1% Более 1%

Рисунок 1. Представленность клонов Т-лимфоцитов разной частоты встречаемости в индивидуальных репертуарах тенр.

Разноцветные столбцы представляют собой Т-клеточное гомеостатическое пространство каждого из девяти доноров. В них цветами показана частота встречаемости клонов Т-лимфоцитов.

Светло-зеленые столбцы соответствуют доле наивных С027Ыб'1С045РАЫг'1 Т-клеток, измеренную с помощью проточной цитофлуориметрии.

Темные горизонтальные линии отсекают долю низкочастотных ТС(ф клонотипов, эквивалентную содержанию наивных Т-лимфоцитов в периферической крови каждого донора. Квадратные скобки показывают долю высокочастотных клонотипов (более 0.001% от общего числа прочтений секвенирования).

1.1. Распределение частот использования \/р сегментов

Существует 54 варианта генов УР-сегмента у человека, вносящих большой вклад в общее индивидуальное разнообразие ТС(*р. Гены этого сегмента кодируют СОИ и СО(<2 участки ТСИР, взаимодействующие с молекулой МНС в составе комплекса пептид-МНС. Мы проанализировали относительное распределение использования УР сегментов в

репертуарах ТСЯР пар мама-ребенок на трех уровнях: на уровне репертуаров нефункциональных ТСЯР, репертуаров ТСЯР низкочастотных клонотипов и репертуаров ТСК|3 высокочастотных клонотипов.

1) Распределение частот использования Vв сегментов в репертуарах нефункциональных ТСИв

Нефункциональные варианты ТСЯР - гены со сбитой рамкой считывания или со стоп-кодоном в области СО?в - находятся на гомологичной хромосоме примерно в 40% Т-лимфоцитов. Такие варианты присутствуют наряду с функциональным вариантом Т-клеточного рецептора в качестве пассивного балласта, не экспрессируются, и, соответственно, не подвергаются отбору в тимусе. Поэтому по их репертуару можно судить об особенностях исходной сборки ТСВ в ходе рекомбинационных событий.

КЛОНОТИПЫ со сбитой рамкой клонотипы считывании

1 Высокочастотные Все клонотипы клонотипы с правильной

рамкой считыванин

Клонотипы НИ1КОЧаСТО!НМ1 со сбитой рамкой клонотипы считывания

высокочастотный Все клонотипы клонотипы с правильной

рамкой считыванин

Рисунок 2. Средние значения дивергенции Дженсена-Шеннона представленности V)}-

сегментов в репертуарах ТСЯ родственных {Р, серые столбцы) и неродственных (НР, синие

столбцы) пар мама-ребенок (с указанием стандартного отклонения). Р-значения с

поправкой Бонферрони представлены с учетом множественных выборок.

(А) Сравнение представленности \/р-сегментов среди общего числа ТС(ф клонотипов (но не

ридов - то есть влияние относительных размеров исследуемых клонотипов ТСВР не

учитывается).

(Б) Сравнение представленности \/р-сегментов среди общего числа ридов -то есть с учетом относительных размеров исследуемых клонотипов ТСКР.

Чем меньше значение дивергенции Дженсена-Шеннона, тем ближе распределение использования УР-сегментов в репертуарах ТС[< исследуемых пар. * - различия статистически достоверны.

Сравнительный анализ репертуаров ТС^Р мам и детей на этом уровне показал, что распределение использования Ур-сегментов сопоставимо для родственных и неродственных пар мама-ребенок. У родственных пар мама-ребенок наблюдается некоторое (статистически недостоверное) уменьшение значения дивергенции Дженсена-Шеннона2 (Рисунок 2, два первых столбца), что может свидетельствовать о большем сходстве распределения частот использования Ур-сегментов в ходе рекомбинации.

2) Распределение частот использования Ув сегментов в репертуарах ТСКв низкочастотных клонотипов

Мы предположили, что влияние селекции в тимусе на распределение представленности \/р-сегментов можно оценить, сравнив репертуар нефункциональных ТС(*р с репертуаром ТС1*Р наивных Т-клеток.

Долю наивных Т-лимфоцитов (С034С027Ы8ЬС045ВАЫ811) в образцах периферической крови каждого донора мы определяли с помощью метода проточной цитофлуориметрии.

Мы предположили, что фракция низкочастотных клонотипов, составляющая ту же долю общего количества ридов (полученного для каждого из доноров), что и определенная с помощью метода проточной цитофлуориметрии доля наивных Т-клеток в периферической крови этого донора, включает в себя преимущественно наивные Т-клетки (Рисунок 1).

Распределение использования УР-сегментов на этом уровне рассмотрения было более сходным для всех доноров, чем на уровне нефункциональных клонотипов (Рисунок 2). По всей вероятности, это является следствием дальнейшего формирования репертуара наивных Т-лимфоцитов под действием их отбора в тимусе. В таком случае его влияние приводит к выраженной конвергенции частот использования различных Ур-сегментов как в родственных (Р = 0.0006), так и в неродственных (Р = 0.0015) парах.

Представленность УР-сегментов в репертуарах функциональных ТС(*Р на уровне низкочастотных Т-лимфоцитов (то есть после позитивной и негативной селекции в тимусе) оказалась достоверно (Р = 0.02) более сходной у родственных, чем у неродственных пар (Рисунок 2, столбцы 3 и 4).

В соответствии с анализом Дженсена-Шеннона, в репертуарах ТСИР родственных групп доноров заметна тенденция к выравниванию распределения использования УР-

2 Дивергенция Дженсена-Шеннона - статистический метод для определения близости двух

распределений вероятностей.

сегментов в низкочастотных клонотипах по сравнению с нефункциональными клонотипами. Интересно, что наблюдаемые изменения в целом схожи и для неродственных доноров. Зависимость степени сближения частоты использования УЗ-сегментов от количества общих аллелей ША оказалась недостоверной (Я = 0.12, Р = 0.63). То есть, по всей вероятности, тимус оказывает значительное и схожее для разных индивидуумов давление на формирование наивного репертуара ТСИР в отношении предпочтительно отбираемых У-сегментов.

3) Распределение частот использования 1/6 сегментов в репертуарах ТСЯв высокочастотных клонотипов

Мы предположили, что влияние встречи Т-лимфоцитов со специфичным антигеном на разнообразие используемых \/-сегментов может быть отслежено при сравнении репертуаров ТС Я наивных и активированных Т-лимфоцитов. В популяцию высокочастотных клонотипов (преимущественно состоящую из активированных Т-лимфоцитов) мы включили те, доля прочтений секвенирования которых превышает 0.001% от всех полученных СОЯЗ-содержащих сиквенсов. Таким образом, нижняя граница этой группы приблизительно на порядок превосходит верхнюю границу группы низкочастотных клонотипов каждого из доноров (Рисунок 1). Это позволяет снизить количество "засоров" наивными Т-клонотипами группы высокочастотных клонотипов.

Мы показали, что высокочастотные клонотипы характеризуется повышенной гетерогенностью используемых ур-сегментов. При этом достоверно большего сходства распределения частот использования ур-сегментов в репертуарах родственных пар по сравнению с неродственными обнаружено не было.

Повышенная гетерогенность используемых ур-сегментов можетявляться следствием увеличения разброса за счет сравнительно небольшого количества высокочастотных клонотипов, входящих в анализ. Также его причиной может быть случайный характер встречи специфичных Т-клонотипов и антигенов, приводящий к экспансии разных клонотипов в разных индивидуумах - даже при иммунном ответе против одних и тех же антигенов.

1.2. Пересечения репертуаров ТСЯр СОИЗ в родственных и неродственных парах мама-ребенок

Исследования последних лет выявили значительные межиндивидуальные пересечения репертуаров ТСЖр. Тем не менее, до сих пор неизвестно, влияет ли гаплоидентичность доноров (то есть идентичность половины аллелей ША) на степень пересечения их репертуаров ТСй.

Чтобы ответить на этот вопрос, мы провели сравнительный анализ репертуаров ТСРр родственных и неродственных пар мама-ребенок на трех уровнях:

1) на уровне аминокислотной последовательности СРРЗ,

2) на уровне нуклеотидной последовательности СЭЯЗ,

3) на уровне нуклеотидной последовательности СОРЗ при совпадающих V- и .1-сегментах (т.е. полностью идентичных ТСИр цепей).

Мы оценили такие пересечения отдельно для низкочастотных и высокочастотных клонотипов, а также для всех клонотипов с правильной рамкой считывания. Для корректного сравнения степени перекрывания репертуаров ТСЖР в разных субпопуляциях, мы нормализовали количество совпадающих СОЯЗ вариантов на объемы сравниваемых выборок по формуле:

[нормализованное количество общих ТС1*р С01*3 клонотипов выборок А и Б] = [количество общих ТСКР С01?3 клонотипов выборок А и Б] / ([количество ТСЯР СЭЯЗ клонотипов в выборке А] х [количество ТС[?Р С01*3 клонотипов в выборке Б])

Нормализованные результаты пересечений представлены на Рисунке 3.

На всех уровнях сравнения пересечений репертуаров СОИЗ относительное количество общих клонотипов было выше для родственных пар, по сравнению с неродственными. Тем не менее, эти различия не достигали статистически значимых величин. Вероятно, для достижения статистически значимых различий необходим анализ большего количества пар доноров.

Наибольшая степень пересечения репертуаров ТСЯР СР1*3 у всех рассматриваемых пар доноров наблюдается в субпопуляции высокочастотных клонотипов. Это может быть следствием сочетания двух причин.

Во-первых, наивные кпонотипы с упрощенной структурой СОЙЗ продуцируются в тимусе с повышенной частотой - как внутри одного индивидуума, так и в рамках всей

популяции. Это повышает вероятность встречи таких клонотипов со специфичными антигенами и, как следствие, их экспансии.

Во-вторых, поскольку все изучаемые нами доноры живут в схожей среде, их иммунная система сталкивается в основном с одними и теми же антигенами. Отбор на уровне специфичности к этим антигенам повышает вероятность клональной экспансии схожих клонов в гомеостатическом пространстве Т-лимфоцитов разных доноров.

1СИ

§ 10"1

а

\о о о

ю-6

,-7.

10

£ 10"8-

10"!

р НР Р нр р НР

низко- высоко- все частотные частотные клонотипы

Аминокислотные клонотипы С0!)3

Р НР Р НР Р НР

низко- высоко- все частотные частотные клонотипы

Нуклеотидные клонотипы СО!»

Р НР Р НР Р НР

низко- высоко- все частотные частотные клонотипы

Нуклеотидные клонотипы СО!», идентичные по У|5- и сегментам

Рисунок 3. Относительное количество общих ТСРр СЭИЗ клонотипов пар родственных и неродственных доноров.

Пересечения показаны на уровне аминокислотных последовательностей СЭНЗ, нуклеотидных последовательностей СЭРЗ и нуклеотидных последовательностей СЭИЗ с идентичными Ур-и ,1Р-сегментами; для родственных (серые круги) и неродственных (синие круги) пар мама-ребенок; для низкочастотных, высокочастотных и всех клонотипов с правильной рамкой считывания. На графике представлены нормализованные значения пересечений (см. в тексте). Р - родственные пары, НР - неродственные пары.

1.3. Одинаковые аминокислотные последовательности CDR3 чаще сочетаются с совпадающим Vß-сегментом в высокочастотных клонотипах родственных пар доноров

Известно, что участок CDR3 преимущественно взаимодействует с пептидом антигена, в то время как участки CDR1 и CDR2, последовательность которых закодирована в V-сегменте, преимущественно взаимодействуют с молекулой МНС. Участки некоторых Vß-сегментов, участвующие в формировании CDR3, очень похожи друг на друга. Поэтому разные Vß-сегменты могут формировать идентичные аминокислотные последовательности CDR3.

В анализ вошло по 10 ООО наиболее часто встречающихся аминокислотных CDR3 клонотипов каждого из доноров. Мы вычисляли пересечения репертуаровэтихклонотипов в родственных и неродственных парах, а также их пересечения при условии совпадения Vß-сегментов (то есть при условии совпадения CDR1, CDR2, и CDR3 участков TCRß).

Мы определяли соотношение количества аминокислотных клонотипов с общими CDR1, CDR2, CDR3 участками к количеству аминокислотных клонотипов с общим CDR3 участком (то есть соотношение Vß-CDR3/CDR3). Оно всегда было выше в родственных парах.

Более того, мы наблюдали достоверную положительную корреляцию этого соотношения и количества общих аллелей MHCI у анализируемых доноров. Несмотря на то, что сравнение проводилось на несортированных Т-лимфоцитах, вклад CD8* клеток был достаточным для обеспечения корреляции.

Это наблюдение, по всей видимости, показывает, что для узнавания определенного комплекса пептид-МНС CD8* клетками часто необходима полная конвергенция TCRß цепи.

1.4. Селекция в тимусе снижает среднюю длину CDR3

Сравнение репертуаров функциональных и нефункциональных3 CDR3 показало, что у всех анализируемых доноров последние характеризуются большей средней длиной CDR3 участка (45,6+ 0.4 против 43,3+-0.2) и большим количеством добавленных нуклеотидов (8,6+-0.2 против 7,4+-0.1). То есть репертуар TCRß CDR3, возникающий в результате соматической рекомбинации (параметры которого отражены в репертуаре нефункциональных CDR3), в среднем характеризуется большей длиной. Последующая же

3 Гены TCRß со сбитой рамкой считывания или стоп-кодоном в области CDR3.

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

1.5. Поиск зрелых микрохимерных клонов Т-клеток

На сегодняшний день считается, что в ходе беременности материнские клетки и клетки плода (в частности, предшественники Т-лимфоцитов и зрелые Т-лимфоциты) активно переходят через плаценту и могут сохраняться в организме нового хозяина в течение десятков лет после рождения. Тем не менее, до сих пор неизвестно, могут ли Т-лимфоциты, перешедшие через плаценту в зрелом состоянии, сохраняться в крови нового индивидуума и работать, как функциональные клоны.

Мы предположили, что проводимый нами анализ методом глубокого секвенирования значительной порции репертуаров ТСР?3 матери и ребенка может выявить присутствие клонов Т-клеток, перешедших через плаценту в зрелом состоянии и давших начало функциональной популяции Т-клеток в новом хозяине.

В качестве маркеров таких клеток мы выбрали:

1) Полное совпадение \/р- и .1Р- сегментов в сочетании с идентичными нуклеотидными последовательностями Сй(?3;

2) Количество добавленных нуклеотидов как характеристику вероятности независимого возникновения данной последовательности СО(?3 в разных людях4.

Тем не менее, чувствительность выбранного нами метода оказалась ограничена естественной дисперсией кривых распределения количества добавленных в СЭКЗ нуклеотидов, характеризующих совпадение части репертуара ТС1^Р любых двух людей.

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

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

1.6. Оценка полного пересечения индивидуальных TCRP репертуаров

Мы предприняли попытку оценить полное пересечение репертуаров TCRp двух индивидуумов. Для этого мы экстраполировали данные по пересечениям репертуаров 12 неродственных пар мама-ребенок.

Увеличивая количество входящих в анализ уникальных клонотипов, выбираемых случайным образом из каждого образца, мы анализировали число идентичных вариантов CDR3, встречающихся в этих выборках. При увеличении объемов выборок в арифметической прогрессии количество общих клонотипов в каждой паре росло в геометрической прогрессии (Рисунок 4 А-В).

При максимальной глубине анализа (порядка одного миллиона клонотипов для каждого из доноров) было обнаружено 72 ООО клонотипов, идентичных на аминокислотном уровне; 68 ООО клонотипов, идентичных на аминокислотном уровне при условии несовпадения нуклеотидных последовательностей; и б ООО клонотипов, идентичных на нуклеотидном уровне. Такое количество общих клонотипов в несколько раз превосходит предыдущие оценки пересечения разнообразия репертуаров человека.

Экстраполируя эти данные на оцениваемое общее индивидуальное разнообразие TCRP, мы провели приблизительную оценку полного перекрывания CDR3 репертуаров TCRP двух индивидуумов, составившую порядка 2 200 ООО аминокислотных последовательностей CDR3; 2 060 ООО аминокислотных последовательностей CDR3 с несовпадающими нуклеотидными последовательностями и 180 ООО нуклеотидных последовательностей CDR3 (что представляет собой 44.1%, 41.3% и 3.6% индивидуальных репертуаров, соответственно).

Эти наблюдения вносят существенный вклад в понимание адаптивной иммунной системы человека. По всей видимости, частотные предпочтения в ходе рекомбинации и селекция в тимусе настолько сильно влияют на формирование репертуара человека, что большинство CDR3 вариантов TCRP, экспрессируемых в наивных Т-лимфоцитах, покидающих тимус, выбираются из некоторого "шорт-листа", состоящего из менее, чем 108 вариантов аминокислотных последовательностей.

50.000 500,000 5,000,000

Количество уникальных АК покледовательностей в образце

50,000 500,000 5,000,000

Количество уникальных НТ последовательностей в образце

50,000 500,000 5,000,000

Количество уникальных АК последовательностей в образце

Количество образцов, содержащих клонотип

Рисунок 4. Перекрывание индивидуальных ТС(ф СОИЗ репертуаров растет геометрически при увеличении объема анализируемых массивов данных.

На графиках представлено количество общих ТСНР СОЯЗ клонотипов 12 неродственных пар доноров в зависимости от объема анализируемых образцов. Анализ пересечения репертуаров был произведен для: (А) аминокислотных последовательностей (АК), (Б) аминокислотных последовательностей с несовпадающими нуклеотидными последовательностями (АК, не НТ) и (В) нуклеотидных последовательностей (НТ). Каждая из 12 линий разных цветов представляет собой наблюдаемое количество общих ТС1?Р СРКЗ клонотипов каждой пары неродственных доноров. Пунктирной прямой показана степенная зависимость (У = аХЬ), которая использовалась для экстраполяции величины перекрытия индивидуальных репертуаров ТС1^р при анализе полных индивидуальных репертуаров.

(Г) На данном графике представлено распределение количества общих клонотипов для 9 анализируемых людей. То есть показано, сколько общих клонотипов встречается одновременно у разного количества доноров. Из графика видно, что частота, с которой ТСЯр клонотипы встречаются в репертуарах человека, распределена по степенному закону.

2. Разработка метода нормированного анализа репертуаров ТСЯ с использованием молекулярного баркодирования

Любая из стадий подготовки библиотеки генов антител или Т-клеточных рецепторов (ТС1*) к анализу (забор образца клеток, выделение нуклеиновых кислот, ПЦР-амплификация, секвенирование) может приводить к искажению количественной информации о реальном репертуаре.

В нашей работе мы разработали метод для решения этой проблемы. Наш метод объединил в себе принципы разработанной нами ранее технологии подготовки библиотек генов ТСК и технологии введения молекулярных баркодов (Рисунок 5). Молекулярный баркод - это последовательность из случайных нуклеотидов, уникальная для каждой молекулы кДНК.

В методе используется эффект смены матрицы ревертазой. В синтезируемую в процессе обратной транскрипции кДНК происходит встраивание 5'-адаптера, содержащего участок из 12 случайных нуклеотидов (молекулярный баркод. Рисунок 5). Таким образом, каждая синтезированная молекула кДНК оказывается специфически мечена одним из 412 (>16.7 миллионов) вариантов уникальных молекулярных баркодов.

синтез первой цепи кДНК

УДГ. Э7*С 30 мин

первая амплификация ПЦР, 18 циклов

I «иутретиД праймер I M1St бериллом ■ обри^цл

вторая амплификация ПЦР.

М1 циклов

Illumina л пап тес

Illumina адаптер

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

I >-1 V NON J

покрытие 100НТ Двусторомиим I |

прочтением сеявеииромми* CDfiJ

I_»

необяодимый участок дчя надежного определения (DR).V- и сегмемпм

I_С '

покрытие 100нт прочтение и с праАмерэ Illumina

Рисунок 5. Схема подготовки библиотеки, созданной на основе мРНК, для массированного секвенирования.

(dUNNNNb = молекулярный баркод; ХХХХХ = баркод образца; (N)2-4 - случайные нуклеотиды, добавленные к 5' концу библиотеки с целью создания разнообразия для лучшего распознавания кластеров на секвенаторе Illumina.

При дальнейшем анализе данных массированного секвенирования библиотеки такие молекулярные баркоды позволяют провести точную оценку стартового количества молекул «ДНК в образце, а также значительно снизить уровень искажения количественной информации в ходе подготовки библиотеки генов TCR. В результате снижается разброс данных для независимо полученных реплик и повышается детектируемое разнообразие репертуара ß цепей TCR при фиксированном количестве анализируемых прочтенных последовательностей (Рисунок 6). Кроме того, применение данного метода нормализации позволяет снизить количественные искажения представленности вариантов внутри каждого конкретного массива.

Разработанный метод мы применили в следующей части работы для нормированного сравнительного анализа разнообразия репертуаров ТС1ф в периферической крови здоровых доноров различного возраста.

А

Б

СО. 600.000 п ОС

и

СУМ.7Х

СУ=12.5%

СУ=25.4Ч

2x10* случайным

прочтений секвенирования

2*10* прочтений секвенирования с уникальными кДНК метками

1x10* случайных

прочтений секвенирования

1x10* прочтении секвенирования с уникальными кДНК метками

Рисунок 6. Нормализация с использованием молекулярного баркодирования увеличивает разнообразие ТС[?0 внутри образца и уменьшает разброс между образцами.

СУ - коэффициент вариации.

(А) Библиотека ТС КЗ, полученная из трех независимо полученных образцов крови одного и того же донора объемом 2 мл. Левы й столбец соответствует количеству уникальных ТС(ф клонотипов, обнаруженных в 2 х 105 случайно выбранных, СОЯЗ-содержащих прочтений секвенирования (без учета молекулярных баркодов). Правый столбец соответствует аналогичному количеству прочтений секвенирования, характеризующихся уникальными молекулярными баркодами.

(Б) Библиотека ТОф, полученная из двух образцов с разным количеством РВМС - 5 х 106 и 1 х 107 клеток. Оба образца были получены в результате одного забора крови у донора. Левый столбец соответствует количеству уникальных ТС1*3 клонотипов, обнаруженных среди 1 х 106 случайно выбранных, СРИЗ-содержащих прочтений (без учета молекулярных баркодов). Правый столбец соответствует аналогичному количеству прочтений с уникальными молекулярными баркодами.

3. Исследование изменений репертуара ТСК|3 в ходе старения человека

В нашем исследовании влияния старения на индивидуальный репертуар ТСИР человека принимало участие 39 системно здоровых донора, разделенных на 4 возрастные группы:

Группа 1 - молодые доноры (в среднем -16 лет),

Группа 2 - доноры среднего возраста (в среднем - 39 лет),

Группа 3 - пожилые доноры (в среднем - 62 года).

Группа 4 - долгоживущие доноры (в среднем - 82 года).

3.1. Наблюдаемое индивидуальное разнообразие ТС1ф СйРЗ падает в течение всей жизни человека

Мы показали практически линейное убывание индивидуального разнообразия ТСНР (наблюдаемого на 1 миллион взятых Т лимфоцитов) с возрастом (Рисунок 7) - начиная от примерно 600 ООО клонотипов в детском возрасте, с кинетикой снижения разнообразия примерно 5 ООО клонотипов в год. Наблюдается статистически значимое уменьшение числа уникальных клонотипов в группе доноров среднего возраста по сравнению с группой молодых доноров (Р<0.01, падение в 1.3 раза), а также в группе пожилых доноров относительно группы доноров среднего возраста (Р<0.01, падение в 1.6 раз). Однако в группе долгоживущих доноров - по сравнению с группой пожилых доноров, -напротив, наблюдается небольшое увеличение разнообразия ТС(ф, которое, впрочем, не является статистически достоверным.

3.2. Определение нижней границы общего индивидуального разнообразия ТС1ф

СОЯЗ

Для оценки нижней границы общего индивидуального разнообразия ТСЯР в периферической крови изучаемых доноров, нами была применена модель ненаблюдаемых видов5. Были получены следующие - последовательно снижающиеся -значения оценочной нижней границы общего разнообразия ТСРР для разных возрастных групп:

Группа 1 - ~7х106 клонотипов;

Группа 2 - ~4х10б клонотипов;

5 Модель, позволяющая на основании доступной части видов оценить их полное

разнообразие.

Группа 3 - ~2.4х106 клонотипов.

Нижняя оценочная граница для доноров Группы 4 оказалась чуть выше, чем у доноров Группы 3, и составила ~2.6х10б клонотипов.

Следует отметить, что мы получили лишь нижнюю оценку разнообразия, которая, вероятно, будет возрастать при увеличении глубины анализа. Чтобы проверить это предположение, мы использовали разделенные на 2 реплики (на уровне периферических мононуклеарных клеток крови) образцы двух доноров возрастом 25 и 87 лет. Доля наивных Т-клеток (CD45RAhigh/CD27hiEh) в образцах периферической крови этих доноров составляла 44% и 15%, соответственно.

В анализе использовалось 106 уникальных стартовых молекул кДНК каждой из реплик. Сначала мы оценивали детектируемое разнообразие и нижнюю границу разнообразия СЭКЗ ТСКР клонотипов, основываясь на одной из двух реплик. Затем в анализ добавляли вторые реплики. Их добавление увеличивало напрямую детектируемое разнообразие клонотипов в 1.91 и в 1.74 раза, а оценку нижней границы разнообразия — в 2.3±0.4 и в 2.4±0.1 раза для молодого и долгоживущего доноров, соответственно.

Таким образом, оценка нижней границы общего индивидуального разнообразия ТСКР, основанная на одном миллионе Т-клеток, действительно, ниже реального разнообразия ТСЯР человека. Скорее всего, с дальнейшим увеличением глубины секвенирования эта оценка будет и далее увеличиваться. Однако, учитывая клональносгь выходящих из тимуса наивных Т-лимфоцитов, воспроизводимость наиболее частых вариантов рекомбинации, а также конечный размер общего пула Т-лимфоцитов человека (до 1012 клеток), общее индивидуальное разнообразие ТСЯР не должно превышать 108-109 вариантов.

3.3. Заполнение гомеостатического пространства активированными клонами Т-клеток

При сравнительном анализе репертуаров ТСКР групп молодого, среднего и пожилого возраста (Группы 1, 2 и 3) прослеживается возрастная зависимость объема Т-клеточного гомеостатического пространства, занимаемого высоко-представленными клонами. Напротив, процентное содержание единично представленных клонотипов (то есть клонотипов, кДНК которых представлена 1 молекулой в анализируемом миллионе молекул кДНК образца), падает с возрастом, составляя в среднем 55% в Группе 1 и 24% - в Группе 3 (Рисунок 8).

А

В

Б

3 _ '

Уникальные молекулы кДНК ТС1*0

Рисунок 7. Возрастные изменения разнообразия репертуаров ТС1ЧЗ. Р - коэффициент корреляции.

(A) На графике изображена зависимость количества уникальных ТСРР СОРЗ клонотипов на 106 Т-клеток от возраста донора. Наблюдаемое разнообразие падает с увеличением возраста практически линейно и достоверно = -0.75; Р < 0.0001, критерий Краскела-Уоллиса). ** - Р < 0.01, двусторонний Ьтест для зависимых выборок; нд - недостоверно. (Б) Кривые видового накопления ТС(*Р клонотипов для четырех возрастных групп, с приращением 1 х 105. В анализ включены 106 молекул кДНК каждого образца. Графики показывают среднее количество уникальных ТСРр СЭКЗ клонотипов, полученных методом случайной выборки (п=5) молекул кДНК из образца каждого донора.

(B) Оценка нижней границы суммарного разнообразия ТС1*3 доноров в зависимости от возраста.

Группа 2

Молодые доноры Доноры среднего Пожилые доноры Долгоживущие (6-25) возраста (34-48) (61-66) доноры (71-90)

Тип Клона Размер Клона Корреляция (R) с возрастом (<70 лет)

■ Редкий единичные клоны -0.8, ***

■ Маленький больше, чем редкие, но меньше, чем 0.001% -0.2, нд

Средний от 0.001% до 0.01% 0.6,***

Большой от0.01%до1% 0.6,***

■ Очень большой более 1% 0.5,**

Рисунок 8. Распределение Т-клеточного гомеопатического пространства доноров разных возрастов между клонами разной представленности.

На графике показана доля гомеостатического пространства, занимаемого клонами разной представленности для четырех возрастных групп. Представленность клонов указана в процентах от общего количества анализируемых Т-клеток (106). *** - Р < 0.0001, ** - Р = 0.001, нд - недостоверно. Я - коэффициент корреляции.

Учитывая относительно стабильное общее количество Т-лимфоцитов (практически не изменяющееся с возрастом), разрастание субпопуляции активированных Т-клеток неизбежно приводит к снижению количества наивных Т-клеток. Это приводит к снижению вероятности эффективного иммунного ответа против новых патогенов или патогенных клеток.

Интересно, что эти закономерности не распространяются на состоящую из относительно долгоживущих доноров Группу 4. В последней, напротив, наблюдается снижение объема, занимаемого высоко-представленными клонами и увеличение количества редко встречающихся клонов, по сравнению с Группой 3 (Рисунок 8).

3.4. Корреляция между разнообразием CDR3 TCRP и долей наивных Т-клеток в периферической крови

Чтобы оценить корреляцию между процентом наивных Т-клеток в периферической крови, возрастом донора и наблюдаемым разнообразием TCRP, мы измерили долю наивных Т-клеток (CD45RAh'gh/CD27high) каждого из исследуемых образцов методом проточной цитофлуориметрии.

Как и предполагалось, была обнаружена обратная зависимость доли наивных Т-клеток в периферической крови донора от его возраста (Рисунок 9). Мы также показали наличие выраженной корреляции между долей наивных Т-клеток и наблюдаемым разнообразием TCRP (Рисунок 10).

Эта корреляция подтверждает как высокую точность методики молекулярного баркодирования для анализа репертуаров TCRp, так и верность определения CD45RAhigh/CD27high субпопуляции Т-лимфоцитов человека, как наивных клеток.

3.5. Образцы периферической крови долгоживущих доноров характеризуются высокой долей CD4* наивных Т-клеток

При анализе доли наивных Т-лимфоцитов (CD45RAWgh/CD27high) в CD8* субпопуляции периферической крови доноров оказалось, что она падает линейно с возрастом донора (R = -0.88; Р < 1х10"8; Рисунок 9). Доля наивных Т-клеток в CD4+ субпопуляции также падает линейно (R = -0.75) и со схожей кинетикой, но лишь у доноров моложе 70лет. В группе же долгоживущих доноров происходит увеличение процента содержания наивных Т-лимфоцитов в CD4* субпопуляции в 1.7 раз по сравнению с донорами Группы 3 - различие близкое к достоверному (Р = 0.054).

Кроме того, в отличие от других рассматриваемых в исследовании возрастных групп, в Группе 4 доля наивных Т-клеток в CD4* субполуляции в 3.6 раз превышает их долю в CD8+ субпопуляции. Более высокая доля наивных CD4* Т-клеток сочетается с большим количеством низкочастотных клонотипов и сравнительно высоким разнообразием TCRP 8 Группе 4 (см. раздел 3.3).

20 40 60 60

Возраст, годы

Возраст, годы

Рисунок 9. Возрастное изменение содержания наивных Т-клеток в периферической крови.

(A) Способ определения наивных Т-клеток методом проточной цитофлуориметрии (FACS). Периферические мононуклеарные клетки крови были проанализированы с помощью проточной цитофлуориметрии. Т-клетки определяли с помощью антител к CD3. Среди CD3* клеток в дальнейшем выделяли наивные Т-клетки, используя антитела к CD45RA и CD27. На графиках представлены типичные данные анализа доли наивных Т-клеток в пуле CD3+ клеток для Групп 1-4.

(Б) Процентное содержание наивных Т-клеток (CD45RAhlgh/CD27high) 8 пуле CD3+ клеток показано в зависимости от возраста донора. Наблюдаемая зависимость линейна и статистически достоверна (R = -0.85, Р = 1 х 10 s) для Групп 1-3 (синие круги), но не для Группы 4 (серые круги). ** - Р = 0.0012, нд- недостоверно.

(B) Процентное содержание наивных Т-клеток в пуле CD8+ клеток (серые круги, пунктирная линия) и в пуле CD4+ клеток (синие круги, сплошная линия для Групп 1-3 и для среднего Группы 4). Доля наивных Т-клеток в CD8' субпопуляции линейно убывает с возрастом (R = -0.88). Однако в CD4+ субпопуляции такое убывание характерно только для Групп 1-3 (R = -0.75). В Группе 4 в CD4+ субпопуляции доля наивных клеток статистически достоверно превосходит их долю в CD8* субпопуляции. ** - Р = 0.004, двусторонний t-тест для зависимых выборок.

900,000

800,000

^ 700,000 -о

g 600,000

£ 500,000 -

§ 400,000 -

OJ

" 300,000 -200,000 -100,000 -

R = 0.9

20 40

Доля наивных Т-клеток, i

60

80

Рисунок 10. Наблюдаемое разнообразие тсяр пропорционально доле наивных Т-клеток в периферической крови.

Показана зависимость разнообразия ТСРр, наблюдаемого в 106 Т-клеток периферической крови, и доли наивных Т-клеток (CD45RAh'eh/CD27high) в пуле СОЗ* клеток (К = 0.90. Р = 1.4 х 10~14). Р - коэффициент корреляции.

Вероятным объяснением этих наблюдений является тот факт, что доноры Группы 4, дожив до глубокой старости, тем самым уже прошли некоторый отбор. То есть их долголетие является следствием (а не причиной) отклонений, наблюдаемых в их периферической крови.

В этом случае можно предположить связь повышенной доли наивных Т-клеток в CD4+ субпопуляции с увеличенной продолжительностью жизни. Такая зависимость может объясняться тем, что более разнообразная популяция CD4+ Т-лимфоцитов в целом лучше регулирует работу иммунной системы, ограничивая развитие аутоиммунных реакций, частых среди пожилых людей, и обеспечивая более сбалансированный выбор типа и интенсивности иммунного ответа.

3.6. Общественные клонотипы6 и старение

Сравнительный анализ репертуаров ТСЯР на 1 миллион Т лимфоцитов для каждых двух доноров оказал значительное пересечение этих репертуаров, состоящее из 8 ООО - 55 ООО идентичных аминокислотных вариантов СРЯЗ и 500 - 7 ООО идентичных нуклеотидных вариантов С01}3.

Чтобы получить список наиболее распространенных клонотипов ТС1?Р в популяции, был проведен сравнительный анализ репертуаров исследуемых доноров. В анализ вошли по 100 ООО наиболее часто встречающихся клонотипов каждого из образцов. Из них были выбраны аминокислотные варианты СОЯЗ, найденные как минимум в 6 из 39 доноров (т.е. в 15% всех доноров). Так, был получен список из примерно 10 ООО самых часто встречающихся клонотипов в популяции.

Клонотипы, входящие в этот список, характеризуются короткой длиной СОИЗ с медианой в 39 нт (в сравнение с медианой длин С01?3 всех проанализированных клонотипов из всех доноров, составляющей 45 нт). Кроме того, среднее количество добавленных к СЭКЗ букв также ниже у этой группы клонотипов.

Интересно, что количество таких наиболее распространенных общественных клонотипов в изучаемой популяции, наблюдаемых среди 100 ООО самых представленных клонотипов, уменьшается от Группы 1 к Группе 3 и коррелирует с наблюдаемым разнообразием репертуара ТСИР (И = 0.48, Р<0.01).

6 Клонотипы, встречающиеся в организмах многих индивидуумов.

Выводы

1. Был проведен сравнительный анализ структуры репертуаров ТС1*Р по частотам использования УР-сегментов для родственных и не родственных пар мама-ребенок. Было показано, что:

- в исходно генерируемом в ходе рекомбинации репертуаре ТС1?р относительная частота использования ур-сегментов сопоставима для родственных и не родственных пар;

- давление селекции в тимусе значительно и сонаправленно для всех доноров изменяет частоты использования УР-сегментов при формировании наивного репертуара ТСВР;

- в репертуаре наивных Т-лимфоцитов (после селекции в тимусе) частоты использования Ур-сегментов достоверно ближе для родственных пар.

2. Впервые был проведен глубокий сравнительный анализ репертуаров последовательностей СОЯЗ ТС(ф для родственных и неродственных пар мама-ребенок на аминокислотном и нуклеотидном уровнях. Было показано, что:

- достоверных различий в пересечении репертуаров для родственных и неродственных пар не наблюдается;

- наибольшая степень пересечения репертуаров ТСИР СОЯЗ для любых двух индивидуумов наблюдается на уровне высокочастотных клонотипов;

- на уровне высокочастотных клонотипов идентичные аминокислотные последовательности СОЯЗ в родственных парах мама-ребенок чаще сочетаются с идентичными УР-сегментами, и эта величина хорошо коррелирует с числом общих аллелей МНС1;

- селекция в тимусе снижает среднюю длину СОЯЗ участка репертуара ТСЯР человека;

- микрохимерных клонов Т-лимфоцитов, перешедших через плаценту в зрелом состоянии, не наблюдается;

- степень пересечения аминокислотных репертуаров ТСЙр СОЯЗ для любых двух людей может составлять до 40% от общего индивидуального разнообразия.

3. Был разработан метод на основе молекулярного баркодирования библиотек, который позволяет корректно и глубоко сравнивать разнообразие репертуаров ТСЯЗ С01*3 образцов, полученных от различных пациентов, в условиях разного количества клеток на старте, разного качества РНК, и разной глубины секвенирования.

28

4. С использованием разработанного метода было показано, что относительное разнообразие репертуара СОЯЗ ТС1*Р периферической крови (на 1 миллион Т-лимфоцитов) практически линейно убывает с возрастом, начиная от ~6х105 клонотипов в детском возрасте, с кинетикой снижения разнообразия ~5х103 клонотипов в год. Снижение разнообразия коррелирует с долей наивных Т-лимфоцитов, которая также линейно убывает с возрастом.

5. На основе данных по разнообразию последовательностей СО(*3 ТСРР на 1 миллион Т-лимфоцитов была проведена оценка нижней границы общего разнообразия ТСЯР в периферической крови человека в зависимости от возраста. Эта величина составила ~7.0х106 вариантов у детей и ~2.5х10б у пожилых людей.

6. Было показано, что образцы периферической крови доноров в возрасте >70 лет характеризуются:

- повышенным процентным содержанием наивных клеток в популяции С04+ Т-лимфоцитов;

- повышенным содержанием низкочастотных клонотипов;

- относительно высоким разнообразием СййЗ ТСИр.

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

7. Были идентифицированы последовательности около 10 ООО наиболее часто встречающихся в популяции клонотипов ТС1*р. Было показано, что:

- такие клонотипы характеризуются малой длиной СйИЗ и малым количеством добавленных нуклеотидов;

- содержание таких клонотипов в периферической крови снижается с возрастом.

Список работ, опубликованных по теме диссертации

Публикации

1. Putintseva EV, Britanova OV, Staroverov DB, Merzlyak EM, Turchaninova MA, Shugay M, Bofotin DA, Pogorelyy MV, Mamedov IZ, Bobrynina V, Maschan M, Lebedev YB, Chudakov DM. Mother and child T cell receptor repertoires: deep profiling study. Front Immunol. 2013; 4:463.

2. Britanova OV*, Putintseva EV*, Shugay M*, Merzlyak EM, Turchaninova MA,Staroverov DB, Bolotin DA, Lukyanov S, Bogdanova EA, Mamedov IZ, Lebedev YB, Chudakov DM. Age-related decrease in TCR repertoire diversity measured with deep and normalized sequence profiling. J Immunol. 2014;192(6):2689-98.

3. Shugay M, Bolotin DA, Putintseva EV, Pogorelyy MV, Mamedov IZ, Chudakov DM. Huge Overlap of Individual TCR Beta Repertoires. Front Immunol. 2013;4:466.

4. Shugay M., Britanova O.V., Merzlyak E.M., Turchaninova M.A., Mamedov I.Z., Tuganbaev T.R., Bolotin D.A., Staroverov D.B., Putintseva E.V., Plevova K., Linnemann C., Shagin D., Pospisilova S., Lukyanov S.„ Schumacher T.N., Chudakov D.M. Towards error-free profiling of immune repertoires. Nat Methods. 2014;ll(6):653-5.

5. Bolotin DA, Shugay M, Mamedov IZ, Putintseva EV, Turchaninova MA, Zvyagin IV, Britanova OV, Chudakov DM. MiTCR: software for T-cell receptor sequencing data analysis. Nat Methods. 2013;10(9):813-4.

6. Mamedov IZ, Britanova OV, Zvyagin IV, Turchaninova MA, Bolotin DA, Putintseva EV, Lebedev YB, Chudakov DM. Preparing unbiased T-cell receptor and antibody cDNA libraries for the deep next generation sequencing profiling. Front Immunol. 2013;4:456.

7. Turchaninova MA, Britanova OV, Bolotin DA, Shugay M, Putintseva EV, Staroverov DB, Sharonov G, Shcherbo D, Zvyagin IV, Mamedov IZ, Linnemann C, Schumacher TN, Chudakov DM. Pairing of T-cell receptor chains via emulsion PCR. Eur J Immunol. 2013;43(9):2507-15.

Тезисы докладов на конференциях

1. Ekaterina Putintseva, Dmitry Bolotin, Miklail Pogorely, Dmitry Staroverov, Mikhail Shugay, Olga Britanova, Dmitry Chudakov "Mother and child TCR repertoires: deep NGS-based profiling study", летний симпозиум ZIBI (Центр инфекционной биологии и иммунологии), июнь 2013, Берлин, Германия.

* Вклад в работу данных авторов равнозначен.

Подписано в печать 4.08.2014 г. Формат А5 Бумага офсетная. Печать цифровая. Тираж 70 Экз. Заказ № 1559-8-14 Типография ООО "Ай-клуб" (Печатный салон МДМ) 119146, г. Москва, Комсомольский пр-кт, д.28 Тел. 8-495-782-88-39