Русская Википедия:Транскриптомные технологии

Материал из Онлайн справочника
Версия от 19:39, 20 сентября 2023; EducationBot (обсуждение | вклад) (Новая страница: «{{Русская Википедия/Панель перехода}} '''Транскрипто́мные техноло́гии''' ({{lang-en|transcriptomics technologies}}) — методы, разработанные для изучения транскриптома (то есть совокупности всех РНК-транскриптов) организма. В состав транскрип...»)
(разн.) ← Предыдущая версия | Текущая версия (разн.) | Следующая версия → (разн.)
Перейти к навигацииПерейти к поиску

Транскрипто́мные техноло́гии (Шаблон:Lang-en) — методы, разработанные для изучения транскриптома (то есть совокупности всех РНК-транскриптов) организма. В состав транскриптома входят все транскрипты, которые присутствовали в клетке на момент Шаблон:Нп5. Исследуя транскриптом, можно установить, какие клеточные процессы были активны в тот или иной момент времени.

Первые попытки изучения транскриптома были предприняты в начале 1990-х годов. Благодаря развитию новых технологий в конце 1990-х транскриптомика стала важной биологической наукой. В настоящий момент в транскриптомике есть два основополагающих метода: микрочипы, позволяющие выявить наличие и количество определённых транскриптов, и секвенирование РНК (РНК-Seq), в котором используются методы секвенирования нового поколения для получения последовательностей всех транскриптов. С улучшением методик количество данных, получаемых в ходе одного транскриптомного эксперимента, увеличивалось. В связи с этим методы анализа данных также совершенствовались, чтобы обеспечить точный и эффективный анализ возрастающего объёма данных. Транскриптомные базы данных постоянно растут и становятся всё более полезными для исследователей. Это связано с тем, что правильная интерпретация данных, полученных в ходе транскриптомного эксперимента, практически невозможна без опоры на предшествующие исследования.

Измерение уровня экспрессии определённых генов в клетках разных тканей и при разных условиях или же в разные моменты времени даёт информацию о регуляторных механизмах, связанных с экспрессией генов. С помощью этих данных могут быть определены функции ранее Шаблон:Нп5 генов. Анализ транскриптомов позволяет выявить различия в экспрессии определённых генов у разных организмов, что может быть особенно полезно для понимания молекулярных основ заболеваний человека.

История

Файл:Transcriptomics technique publications over time-ru.svg
Количество публикаций, касающихся РНК-Seq (чёрный), микрочипов (красный), экспрессируемых меток последовательностей (синий) и сериального/кэпового анализа экспрессии генов (жёлтый) с 1990 года по 2016 год[1]

Первая попытка получения части человеческого транскриптома была предпринята в 1991 году; в ходе этого исследования были получены последовательности 609 мРНК из мозга человека[2]. В 2008 году были опубликованы два человеческих транскриптома, состоящие из миллионов последовательностей, происходящих от транскриптов 16 тысяч генов[3][4]. К 2015 году были опубликованы транскриптомы сотен людей[5][6]. Получение транскриптомов индивидуумов с тем или иным заболеванием, разных тканей и даже одиночных клеток в настоящее время является рутинной процедурой[6][7]. Бурное развитие транскриптомики было возможно благодаря быстрому развитию новых экономичных технологий с повышенной чувствительностью[8][9][10][11].

До транскриптомики

Исследования отдельных транскриптов проводились ещё за несколько десятилетий до того, как методы транскриптомики стали общедоступны. В конце 1970-х были получены Шаблон:Нп5, которые были конвертированы в комплементарную ДНК (кДНК) с помощью обратной транскриптазы для бабочки Шаблон:Нп5[12]. В 1980-х с помощью низкопроизводительного секвенирования по Сэнгеру были получены последовательности случайных транскриптов; так появились так называемые Шаблон:Нп5 (Шаблон:Lang-en)[2][13][14]. Метод секвенирования по Сэнгеру доминировал до появления технологий высокопроизводительного секвенирования, например, секвенирования синтезом (Solexa/Illumina). EST стали активно использоваться в 1990-х как эффективный метод определения генного состава организма без Шаблон:Нп5[14]. Количество отдельных транскриптов оценивалось с помощью нозерн-блотов, Шаблон:Нп5 и количественной ПЦР с обратной транскрипцией (RT-qPCR)[15][16]. Однако эти методы очень трудоёмки и охватывают лишь крошечную долю целого транскриптома[11].

Первые попытки

Слово «транскриптом» (Шаблон:Lang-en) было введено в употребление в 1990-х годах[17][18]. В 1995 году появился первый транскриптомный метод, основанный на секвенировании, — Шаблон:Нп5 (Шаблон:Lang-en), который заключался в секвенировании по Сэнгеру соединённых фрагментов случайных транскриптов. Количество транскриптов оценивалось по числу совпадений с фрагментами известных генов[19]. Вскоре появился вариант SAGE, использующий вместо секвенирования по Сэнгеру технологии секвенирования нового поколения — цифровой анализ экспрессии генов (Шаблон:Lang-en)[8][20]. Однако эти методы практически полностью были вытеснены методами высокопроизводительного секвенирования целых транскриптов, которые давали дополнительную информацию о транскрипте, например, сведения о сплайсинговых вариантах[8].

Развитие современных методов

Сравнение современных методов[21][22][9]
РНК-Seq Микрочипы
Шаблон:Нп5 От 1 дня до 1 недели на эксперимент[9] 1—2 дня на эксперимент[9]
Необходимое количество РНК Низкое ~ 1 нг тотальной РНК[23] Высокое ~ 1 мкг РНК[24]
Трудоёмкость Высокая (пробоподготовка и анализ данных)[9][21] Низкая[9][21]
Предшествующая информация Не требуется, хотя наличие Шаблон:Нп5/транскриптома упрощает работу[21] Для создания Шаблон:Нп5 необходим референсный геном/транскриптом[21]
Точность Шаблон:Нп5 ~ 90 % (ограничена Шаблон:Нп5 последовательностей)[25] > 90 % (ограничена точностью детекции флуоресценции)[25]
Разрешение последовательностей Может детектировать однонуклеотидные полиморфизмы и сплайсинговые варианты (ограничение — точность секвенирования (~ 99 %))[25] Специализированные микрочипы могут детектировать сплайсинговые варианты (ограничения — создание проб и кросс-гибридизация)[25]
Чувствительность 1 транскрипт на миллион (приблизительно, ограничение — покрытие последовательности)[25] 1 транскрипт на тысячу (приблизительно, ограничено детекцией флуоресценции)[25]
Динамический диапазон 100000 : 1 (ограничен покрытием последовательности)[26] 1000 : 1 (ограничен насыщением флуоресценции)[26]
Техническая воспроизводимость > 99 %[27][28] > 99 %[29][30]

Превалирующие современные методы — микрочипы и РНК-Seq — появились в середине 1990-х и 2000-х соответственно[8][31]. Публикации по микрочипам, которые измеряли относительное содержание определённых транскриптов за счёт их гибридизации с комплементарными пробами, нанесёнными на микрочип, появились в 1995 году[32][33]. Метод микрочипов позволял одновременно исследовать тысячи транскриптов, и за счёт этого позволял снизить стоимость исследования транскриптома в расчёте на ген и сэкономить усилия[34]. До конца 2000-х лучшими методами транскрипционного профилинга были пятновые олигонуклеотидные чипы (Шаблон:Lang-en) и микрочипы Шаблон:Нп5 с высокой плотностью[11][31]. В течение этого периода было создано множество чипов, покрывающих известные гены модельных и экономически важных организмов. Улучшения технологий создания микрочипов привели к увеличению специфичности проб и количества генов, которые можно проанализировать с помощью одного чипа. Благодаря новым методам детекции флуоресценции стало возможным точно определять наличие и количество даже транскриптов, синтезируемых на низком уровне[33][35].

РНК-Seq подразумевает секвенирование кДНК, соответствующей транскриптам, причём численность отдельных фрагментов кДНК определяется численностью соответствующих транскриптов. Огромное влияние на РНК-Seq оказало развитие методов секвенирования нового поколения[8][10]. Первым транскриптомным методом стало Шаблон:Нп5 (Шаблон:Lang-en), в основе которого лежало образование коротких последовательностей длиной от 16 до 20 пар оснований (п. о.) в ходе сложной последовательности гибридизаций[36]. В 2004 году с помощью этого метода была оценена экспрессия 10 тысяч генов растения Arabidopsis thaliana[37]. Первая работа, посвящённая РНК-Seq, была опубликована в 2006 году. В ходе этого исследования с помощью технологии Шаблон:Нп5 была определена последовательность ста тысяч транскриптов[38]. Полученного покрытия было достаточно для оценки относительного количества отдельных транскриптов. Популярность РНК-Seq значительно повысилась после 2008 года, когда технологии Illumina/Solexa позволили секвенировать один миллиард транскриптов[4][9][39][40]. Благодаря этим данным сейчас возможно количественно оценивать и сравнивать транскриптомы разных людей[4].

Получение данных

Получение данных о транскриптах возможно двумя принципиально различающимися путями: секвенированием отдельных транскриптов (EST или РНК-Seq) или гибридизацией транскриптов на упорядоченный чип нуклеотидных последовательностей (микрочип)[21].

Выделение РНК

Для всех транскриптомных методов необходимо выделить РНК из исследуемого организма. Несмотря на огромное разнообразие биологических систем, методика выделения РНК во всех случаях примерно одна и та же. Она включает разрушение клеток и тканей, разрушение РНКаз при помощи Шаблон:Нп5 солей[41], разрушение макромолекул и комплексов, содержащих нуклеотиды, отделение РНК от ненужных биомолекул, включая ДНК, концентрирование РНК при помощи Шаблон:Нп5 из раствора и очищение с помощью специальных Шаблон:Нп5[41][42]. Выделенную РНК также можно дополнительно обработать ДНКазой, чтобы разрушить остатки ДНК[43]. Обычно необходимо концентрирование мРНК, поскольку 98 % выделенной РНК приходится на рРНК[44]. Концентрирование можно произвести с помощью методов, использующих наличие у мРНК поли(А)-хвоста, или путём удаления рРНК с помощью специфических проб[45]. На результаты эксперимента может повлиять разрушенная РНК. Например, если отбирать мРНК из повреждённых РНК, то отобранные молекулы могут быть лишены 5'-концов и привести к искажению данных. Чтобы избежать разрушения РНК, перед её выделением образец обычно подвергают Шаблон:Нп5[42].

Экспрессируемые метки последовательностей

Экспрессируемые метки последовательностей (EST) — это короткие нуклеотидные последовательности, полученные из целого транскрипта. Поскольку EST можно получить без какой-либо специфики относительно организма, из которого выделена РНК, их можно получить из смеси организмов или образцов, взятых из окружающей среды[14]. Хотя в настоящее время чаще всего используется высокопроизводительное секвенирование, Шаблон:Нп5 EST активно использовали при разработке первых микрочипов. Например, микрочип ячменя был получен из 350 тысяч предварительно секвенированных EST[46].

Сериальный и кэповый анализ экспрессии генов (SAGE/CAGE)

Файл:Summary of SAGE-ru.svg
Схема SAGE. Из организма выделяется мРНК, переводится в двуцепочечную ДНК (синий) посредством обратной транскрипции. Далее ДНК расщепляется рестриктазами (в точках ‘X’ и ‘X’+11) с образованием 11-нуклеотидных меток. Далее они конкатенируются и секвенируются с помощью метода Сэнгера с длинными прочтениями (метки разных оттенков синего соответствуют разным генам). Далее последовательности подвергают деконволюции и определяют частоту встречаемости каждой метки. Частота метки свидетельствует об интенсивности экспрессии соответствующего гена[47]

Сериальный анализ экспрессии генов является дальнейшим развитием технологии EST с большей выработкой меток. Он также позволяет провести некоторый количественный анализ численности транскриптов. РНК сначала переводится в кДНК, затем она разрезается на метки длиной 11 нуклеотидов с помощью рестриктаз, которые вносят разрывы в определённые последовательности ДНК. Полученные метки сшивают по типу «голова к хвосту» в длинные фрагменты длиной более 500 нуклеотидов, которые секвенируют с помощью низкопроизводительных, но дающих длинные Шаблон:Нп5 методов, таких как секвенирование по Сэнгеру. Далее последовательности снова делят на 11-нуклеотидные кусочки с помощью специальных компьютерных программ (деконволюция)[19]. Если референсный геном недоступен, то полученные метки можно непосредственно использовать в качестве диагностических маркеров, которые в случае болезни экспрессируются иначе, чем в здоровом организме[19].

Кэповый анализ экспрессии генов (Шаблон:Lang-en) представляет собой вариант SAGE, при котором в качестве меток берутся только 5'-концевые последовательности мРНК. Поэтому, когда метки выравниваются на референсный геном, можно идентифицировать точки начала транскрипции генов. Этот метод активно используется для анализа промоторов и для клонирования полноразмерных кДНК[48].

SAGE и CAGE дают информацию о большем количестве генов, чем секвенирование отдельных EST, однако пробоподготовка и анализ данных в этих методах существенно сложнее[48].

Микрочипы

Файл:Summary of RNA Microarray-ru.svg
Схема работы микрочипов. Из организма выделяется зрелая мРНК (красный) и путём обратной транскрипции переводится в двуцепочечную ДНК (синий). Далее ДНК фрагментируется и флуоресцентно метится (оранжевый). Меченые фрагменты связываются с комплементарными олигонуклеотидными пробами на микрочипе, и интенсивность флуоресценции в каждой ячейке свидетельствует о численности соответствующего гена[47]

Принципы и преимущества

Микрочип состоит из коротких олигонуклеотидов (проб), которые прикреплены в ячейках сетки на стеклянной подложке[49]. Многочисленность транскриптов определяется на основании гибридизации Шаблон:Нп5 транскриптов с этими пробами[50]. Шаблон:Нп5 в каждой ячейке свидетельствует о численности транскрипта, гибридизующегося с данной пробой[50].

Для создания микрочипа необходимо знать, хотя бы частично, геном исследуемого организма, например, в виде аннотированной последовательности или библиотеки EST; это необходимо для создания проб[34].

Методы

Микрочипы, использующиеся в транскриптомике, можно подразделить на два типа: пятновые чипы с низкой плотностью и чипы высокой плотности с короткими пробами[34]. Пятновые чипы низкой плотности обычно представляют собой стеклянную основу, на которую нанесены пиколитровые капли, содержащие разные фрагменты очищенной кДНК[51]. Эти пробы длиннее, чем в чипах с короткими пробами, и с их помощью нельзя выявить события альтернативного сплайсинга. В пятновых чипах используются два типа флуорофоров, которыми метят экспериментальные и контрольные образцы, а относительная многочисленность высчитывается из интенсивности флуоресценции[52]. Чипы высокой плотности используют только одну флуоресцентную метку, и каждый образец гибридизуется и детектируется отдельно[53]. Чипы высокой плотности распространялись компанией Affymetrix GeneChip. В этих чипах каждому транскрипту соответствует несколько 25-нуклеотидных проб[54]. Компания NimbleGen производит чипы высокой плотности при помощи Шаблон:Нп5, которая позволяла получать чипы разного строения. Один чип содержит сотни тысяч проб длиной от 45 до 85 нуклеотидов, которые гибридизуются с образцом, меченным флуоресцентной меткой одного вида[55].

РНК-Seq

Файл:Summary of RNA-Seq-ru.svg
Схема РНК-Seq. Из организма выделяется зрелая мРНК (красный) и путём обратной транскрипции переводится в двуцепочечную ДНК (синий). ДНК секвенируется при помощи высокопроизводительных методов с короткими прочтениями. Далее прочтения выравниваются на референсный геном, благодаря чему выявляются транскрибирующиеся участки генома. С помощью полученных данных можно установить, какие гены экспрессируются, какова интенсивность их экспрессии, а также наличие альтернативного сплайсинга[47]

Принципы и преимущества

РНК-Seq представляет собой сочетание высокопроизводительного секвенирования с вычислительными методами оценки численности отдельных транскриптов в экстракте РНК[9]. Обычно получаются последовательности длиной около 100 пар оснований (п. о.), однако в зависимости от метода секвенирования их длина может составлять от 30 п. о. до 10 тысяч п. о. РНК-Seq обеспечивает глубокое покрытие транскриптома множеством коротких фрагментов, благодаря которому возможно при помощи вычислительных методов реконструировать исходные транскрипты, выравнивая прочтения на референсный геном или друг на друга (сборка de novo)[8]. С помощью РНК-Seq можно рассчитать количество как многочисленных, так и малочисленных РНК, так как динамический диапазон метода составляет 5 порядков. В этом заключается главное преимущество РНК-Seq перед микрочипами. Кроме того, для РНК-Seq требуется очень мало исследуемой РНК по сравнению с микрочипами — нанограммы против микрограммов. Благодаря этому РНК-Seq в сочетании с линейной амплификацией кДНК позволяет исследовать очень небольшие клеточные структуры вплоть до отдельных клеток[23][56]. Теоретически верхнего предела количественной оценки в РНК-Seq не существует, и для прочтений длиной 100 п. о. фоновый шум в неповторяющихся участках очень низок[9].

С помощью РНК-Seq можно идентифицировать гены в геноме или установить, какие гены активны в данный момент времени. На основании количества прочтений можно точно установить относительный уровень экспрессии генов. Методология РНК-Seq постоянно совершенствуется, преимущественно за счёт улучшения технологий секвенирования, которые повышают производительность и точность метода, а также выдают прочтения всё большей длины[57]. Со времени первых публикаций в 2006 и 2008 годах[38][58] РНК-Seq интенсивно внедрялся в исследования и к 2015 году догнал микрочипы, став вторым доминирующим транскриптомным методом[59].

Попытки получения транскриптомных данных для отдельных клеток стимулировали совершенствование методов приготовления библиотек для РНК-Seq, что значительно увеличило чувствительность технологии. На данный момент получен ряд транскриптомов единичных клеток и даже появились методы РНК-Seq in situ, в которых транскриптомы отдельных клеток были получены непосредственно в фиксированных тканях[60].

Методы

РНК-Seq появился вместе с бурным развитием нескольких методов высокопроизводительного секвенирования[61]. Однако стадии секвенирования выделенных РНК предшествуют несколько этапов пробоподготовки, которые различаются в разных методах. Методы различаются способами концентрирования транскриптов, фрагментации, амплификации, способом секвенирования (одноконцевое или парноконцевое), а также тем, сохраняется ли информация об исходной цепи[61].

Чувствительность РНК-Seq в конкретном эксперименте можно повысить за счёт концентрирования интересующих классов РНК и удаления остальных. мРНК можно отделить с помощью олигонуклеотидных проб, которые связываются с их поли(А)-хвостами. Удалить неинформативные и чрезвычайно многочисленные рРНК можно с помощью гибридизующихся проб, созданных специально для рРНК данного таксона (например, млекопитающих или растений). Однако вместе с рРНК с помощью такого подхода можно удалить и другие РНК, что может исказить картину эксперимента. Малые РНК, например, микроРНК, можно выделить на основании их размера из агарозного геля после Шаблон:Нп5[62].

Поскольку мРНК, как правило, длиннее, чем единичные прочтения в большинстве методов высокопроизводительного секвенирования, обычно перед секвенированием транскрипты Шаблон:Нп5[63]. Метод фрагментации лежит в основе создания библиотеки для секвенирования. Фрагментацию можно проводить путём химического гидролиза, распыления, обработки ультразвуком (Шаблон:Нп5) или обратной транскрипции с использованием терминирующих нуклеотидов[63]. Кроме того, фрагментацию и добавление меток к кДНК можно осуществить одновременно при помощи транспозаз[64].

В ходе пробоподготовки для секвенирования фрагменты кДНК, соответствующие транскриптам, можно размножить с помощью ПЦР так, чтобы повысить численность молекул, содержащих необходимые 3'- и 5'-концевые Шаблон:Нп5[65]. Стадия амплификации также необходима перед секвенированием образцов с очень низким содержанием РНК. Нижним пределом количества РНК, которое пригодно для секвенирования, является 50 пикограмм[66]. Чтобы оценить качество библиотеки и секвенирования (GC-состав, длина фрагментов, предпочтение фрагментов с определённым положением в транскрипте), можно использовать контрольные Шаблон:Нп5[67]. Шаблон:Нп5 (Шаблон:Lang-en) — это короткие случайные последовательности, которые используются для индивидуального мечения фрагментов при приготовлении библиотеки таким образом, чтобы после добавления идентификатора каждый фрагмент был уникален[68]. С помощью UMI можно измерить численность транскриптов в абсолютной шкале, чтобы скорректировать отклонения, возникшие при создании библиотеки, перед амплификацией, а также чтобы точно оценить количество ДНК в исходном образце. UMI особенно удобны для РНК-Seq одиночных клеток, в которых исходное количество РНК очень мало и требует неспецифичной амплификации[69][70][71].

После пробоподготовки молекулы транскриптов (точнее, соответствующих им кДНК) можно секвенировать в одном направлении (одноконцевое прочтение) или в обоих (парноконцевое прочтение). Одноконцевое секвенирование, как правило, быстрее и дешевле, и в большинстве случаев его достаточно для количественной оценки уровней экспрессии генов. Парноконцевое секвенирование позволяет получить более точные выравнивания и сборки, что очень важно для аннотации генов и описания изоформ транскриптов[9]. Методы РНК-Seq, специфичные по отношению к цепи, сохраняют сведения о цепи ДНК, с которой был транскрибирован каждый транскрипт. Без этой информации прочтения можно выровнять на локус, однако будет неясно, в каком направлении происходит транскрипция гена. Одноцепочечный РНК-Seq удобен для определения направления транскрипции Шаблон:Нп5, расположенных на разных цепях, что позволяет сделать предсказание генов у немодельных организмов более точным[72].

Технологии секвенирования, используемые в РНК-Seq[73][74]
Платформа Коммерческий релиз Типичная длина прочтения Максимальная производительность одного запуска Точность единичного прочтения Запуски РНК-Seq, размещённые в базе NCBI SRA по состоянию на октябрь 2016 года. RNA-Seq runs deposited in the NCBI SRA (Oct 2016)[75]
454 Life Sciences 2005 700 п. о. 0,7 миллиарда п. о. 99,9 % 3548
Illumina 2006 50-300 п. о. 900 миллиардов п. о. 99,9 % 362903
SOLiD 2008 50 п. о. 320 миллиардов п. о. 99,9 % 7032
Ion Torrent 2010 400 п. о. 30 миллиардов п. о. 98 % 1953
Шаблон:Нп5 2011 10000 п. о. 2 миллиарда п. о. 87 % 160

NCBI SRA — Шаблон:Нп5 Национального центра биотехнологической информации (США)

Поскольку в настоящий момент РНК-Seq включает перевод РНК в кДНК в ходе обратной транскрипции, платформы для последующего секвенирования одинаковы и для транскриптомных, и для геномных данных. По этой причине развитие РНК-Seq в значительной мере определяется усовершенствованием методов секвенирования ДНК[74][76][77]. Однако всё большее распространение получает непосредственное секвенирование РНК при помощи нанопор[78][79]. С помощью нанопорового секвенирования в РНК можно обнаружить модифицированные основания, которые нельзя было выявить при секвенировании кДНК, кроме того, для этого метода не нужна амплификация, которая вносит дополнительные искажения[10][80].

Чувствительность и точность РНК-Seq определяются числом прочтений, полученных с каждого образца[81][82]. Для достаточного покрытия транскриптома необходимо очень много прочтений, что позволяет обнаружить даже малочисленные транскрипты. Дополнительные сложности создаёт этап секвенирования, дающего прочтения ограниченной длины, различной точности и качества. Более того, организмы каждого вида имеют разное количество генов, поэтому для эффективной сборки транскриптома для каждого вида необходимо разное число прочтений. На ранних этапах это количество определялось эмпирически, однако с развитием технологий необходимое покрытие стало возможным предсказывать вычислительными методами. Наиболее эффективным способом улучшения точности выявления дифференциальной экспрессии низкоэкспрессируемых генов является не увеличение числа прочтений, а увеличение Шаблон:Нп5[83]. В настоящий момент Энциклопедия элементов ДНК рекомендует 70-кратное покрытие экзома для обычного РНК-Seq и до 500-кратного покрытия для обнаружения редких транскриптов и изоформ[84][85][86].

Анализ данных

Методы транскриптомики позволяют проводить параллельные эксперименты со множеством образцов, поэтому для получения результатов как с помощью РНК-Seq, так и с помощью микрочипов необходима серьёзная Шаблон:Нп5 вычислительными методами[87][88][89][90]. Данные с микрочипов представляют собой изображения с высоким разрешением, поэтому обработка данных включает Шаблон:Нп5 (Шаблон:Lang-en) и спектральный анализ[91]. Изображения, полученные с микрочипов, достигают размера 750 Мб, в то время как обработанные данные занимают 60 Мб. Множество коротких проб, соответствующих одному и тому же транскрипту, могут позволить определить экзон-интронную структуру гена, поэтому для определения достоверности финального сигнала необходимы статистические модели. В ходе экспериментов РНК-Seq получаются миллиарды коротких последовательностей ДНК, которые необходимо выровнять на референсный геном, включающий миллионы или миллиарды п. о. Сборка транскриптома de novo требует построения весьма сложных Шаблон:Нп5[92]. Операции обработки данных РНК-Seq требуют многократного повторения, поэтому для них могут удобны параллелизованные вычисления, однако с использованием современных алгоритмов обработку данных простых транскриптомных экспериментов, не требующих сборки de novo, можно осуществить даже на обычном персональном компьютере[93]. Человеческий транскриптом можно достаточно точно собрать из 300 миллионов 100-нуклеотидных прочтений, полученных с помощью РНК-Seq[81][82]. Для хранения такого объёма данных в сжатом Шаблон:Нп5 требуется 1,8 Гб дискового пространства на образец. Обработанные численные значения для каждого гена занимают ещё меньше памяти, сопоставимо с обработанными данными с микрочипов. Данные о последовательностях можно хранить в публичных данных, таких как SRA (от Шаблон:Lang-en — архив прочтений последовательностей)[94]. Набор данных РНК-Seq можно загрузить в базу данных Gene Expression Omnibus[95].

Обработка изображений

Файл:Microarray and sequencing flow cell-ru.svg
Микрочип и проточная ячейка для секвенирования (РНК-Seq). Обработка изображений при использовании микрочипов и РНК-Seq различается. В микрочипе каждое пятно соответствует определённой олигонуклеотидной пробе, и интенсивность флуоресценции соответствует численности определённого транскрипта (Affymetrix). В ячейке для высокопроизводительного секвенирования пятно соответствует одному нуклеотиду, секвенированному в данный момент времени, и цвет флуоресценции соответствует виду нуклеотида (Illumina Hiseq). В различных видах этих методов может быть использовано большее или меньшее количество цветов[47][96]

Обработка изображений, полученных с микрочипов, должна сохранять Шаблон:Нп5 ячеек изображения и независимо количественно оценивать интенсивность флуоресценции в каждой ячейке. Необходимо также выявлять Шаблон:Нп5 и исключать их из финального анализа. Интенсивность флуоресценции свидетельствует о представленности каждой последовательности, поскольку последовательность пробы в каждой ячейке известна[96].

Первые этапы РНК-Seq также включают схожую обработку изображений, однако перевод изображений в данные о последовательностях производится автоматически специальными программами. Результатом секвенирования синтезом по технологии Illumina является совокупность кластеров, расположенных на поверхности проточной ячейки[97]. В ходе каждого цикла секвенирования изображение каждой проточной ячейки снимается до четырёх раз, причём один запуск включает десятки или сотни циклов. Кластеры проточных ячеек аналогичны пятнам в микрочипах и на ранних этапах секвенирования должны быть правильно определены. При пиросеквенировании (компания Шаблон:Нп5) интенсивность испускаемого света соответствует количеству одинаковых нуклеотидов в гомополимерном участке. Существует множество вариаций перечисленных методов, и каждая предполагает использование разных профилей ошибок для получающихся данных[98].

Анализ данных РНК-Seq

В ходе экспериментов РНК-Seq получается огромный объём прочтений, которые необходимо обработать для получения полезной информации. Анализ данных, как правило, включает использование комбинаций различных биоинформатических программ, которые необходимо подбирать сообразно эксперименту и целям. Процесс обработки данных можно разбить на четыре этапа: контроль качества, выравнивание, количественный анализ и дифференциальная экспрессия[99]. Наиболее популярные программы для обработки данных РНК-Seq запускаются из командной строки в среде Unix или R/Bioconductor[89].

Контроль качества

Прочтения небезупречны, поэтому необходимо определить точность прочтения каждого основания в последовательности. Прочтения, прошедшие контроль качества, гарантированно имеют высокую точность определения каждого основания, их GC-состав соответствует ожидаемому распределению, в них нет перепредставленности коротких мотивов и редки дупликации[82] . Существует несколько программ для анализа качества, например, FastQC и FaQCs. Низкокачественные прочтения или удаляются, или специальным образом помечаются, что учитывается при дальнейшем анализе[100][101].

Выравнивание

Чтобы связать численность прочтений с определённым геном, прочтения необходимо выровнять на референсный геном или друг на друга, если референсный геном неизвестен (сборка транскриптома de novo)[102][103]. Основными требованиями, которым должны удовлетворять программы для выравнивания, являются скорость, позволяющая за приемлемое время выровнять миллиарды коротких прочтений, определённая гибкость, чтобы выявлять случаи сплайсинга эукариотических мРНК, а также правильно выбирать расположение прочтений, соответствующих нескольким местам в геноме. Программы постоянно совершенствуются в соответствии с перечисленными требованиями, а увеличение длины прочтений снижает вероятность неоднозначного выравнивания. Шаблон:Нп5 (EBI) поддерживает список доступных в данный момент инструментов для выравнивания прочтений, полученных при высокопроизводительном секвенировании[104].

Выравнивание первичных эукариотических транскриптов на референсный геном требует особое обращение с интронами, которых нет в зрелых мРНК[105]. Программы для выравнивания коротких прочтений могут создавать особые выравнивания, предназначенные специально для идентификации сайтов сплайсинга на основе канонических последовательностей сайтов сплайсинга. Идентификация сайтов сплайсинга предотвращает их ошибочное выравнивание или отбраковывание, позволяя выровнять на референсный геном большее количество прочтений и увеличивая качество количественных оценок экспрессии генов. Поскольку Шаблон:Нп5 может осуществляться на уровне изоформ мРНК, выравнивания, учитывающие сплайсинг, позволяют обнаружить изменения в численности тех или иных изоформ, что было бы невозможно с помощью обычного анализа[106].

Для сборки транскриптома de novo производится выравнивание прочтений друг на друга, что позволяет реконструировать полноразмерные транскрипты без использования референсного генома[107]. Сложностями сборки de novo являются необходимость больших вычислительных мощностей, чем для сборки на основе референсного генома, дополнительной проверки вариантов и фрагментов генов и дополнительной аннотации собранных транскриптов. Первые метрики, предназначенные для оценки качества сборки транскриптома, такие как N50, были признаны ошибочными[108], и в данный момент доступны усовершенствованные методы оценки. Метрики, основанные на аннотации, хорошо подходят для оценки степени сборки генома. Транскриптом, собранный de novo, может быть использован в качестве референсного при выравнивании последовательностей и количественном анализе экспрессии генов[109][110].

Программы для сборки транскриптома de novo
Программа Дата релиза Дата последнего обновления Вычислительная эффективность Сильные и слабые стороны
Velvet-Oases[111][112] 2008 2011 Низкая, один поток выполнения, необходимо много памяти с произвольным доступом Первый сборщик коротких прочтений. В настоящее время почти не используется.
SOAPdenovo-trans[103] 2011 2014 Средняя, несколько потоков выполнения, умеренная необходимость в памяти с произвольным доступом Один из первых сборщиков коротких прочтений. Адаптирован для сборки транскриптомов.
Trans-ABySS[113] 2010 2016 Средняя, несколько потоков выполнения, умеренная необходимость в памяти с произвольным доступом Предназначена для коротких прочтений, но может быть использована и для сложных транскриптомов. Доступна MPI-параллельная версия для вычислительных кластеров.
Trinity[114][92] 2011 2017 Средняя, несколько потоков выполнения, умеренная необходимость в памяти с произвольным доступом Предназначена для коротких прочтений. Можно использовать для сложных транскриптомов, но требуется много памяти.
miraEST[115] 1999 2016 Средняя, несколько потоков выполнения, умеренная необходимость в памяти с произвольным доступом Может обрабатывать повторяющиеся последовательности, комбинирует несколько форматов данных секвенирования, совместима с большим числом платформ секвенирования.
Newbler[116] 2004 2012 Низкая, один поток выполнения, необходимо много памяти с произвольным доступом Специализирована на устранении ошибок секвенаторов 454 Roche, связанных с гомополимерными последовательностями.
CLC genomics workbench[117] 2008 2014 Высокая, несколько потоков выполнения, низкая необходимость памяти с произвольным доступом Имеет графический интерфейс, может комбинировать различные технологии секвенирования. Не специализирована для транскриптомов, перед использованием необходимо приобрести лицензию.
SPAdes[118] 2012 2017 Высокая, несколько потоков выполнения, низкая необходимость памяти с произвольным доступом Предназначена для транскриптомных экспериментов с единичными клетками.
RSEM[119] 2011 2017 Высокая, несколько потоков выполнения, низкая необходимость памяти с произвольным доступом Может оценивать частоту альтернативно сплайсированных транскриптов. Удобна в использовании.
StringTie[93] 2015 2018 Высокая, несколько потоков выполнения, низкая необходимость памяти с произвольным доступом Может использовать комбинацию методов сборки на основе референсного генома и de novo для идентификации транскриптов.

Количественный анализ

Файл:Transcriptomes heatmap example-ru.svg
Тепловая карта, отображающая наличие одинаковых паттернов экспрессии в различных образцах. Каждая колонка отображает данные измерений изменений экспрессии генов в единичном образце. Относительная экспрессия генов показана цветом: высокая (красный), средняя (белый) и низкая (синий). Гены и образцы с одинаковой экспрессией можно автоматически сгруппировать (левое и верхнее деревья). Образцы могут соответствовать разным индивидуумам, тканям, условиям среды и состоянию здоровья. В данном примере в образцах 1, 2 и 3 экспрессия набора генов 1 высокая, а набора 2 — низкая[47][120]

Количественный анализ выравниваний прочтений может быть произведён на уровне гена, экзона и транскрипта. Типичным результатом анализа является количество прочтений для каждого элемента анализа (гена, экзона или транскрипта). Например, для генов она выдаётся в формате general feature format (GFF). Количество прочтений для генов и экзонов можно определить при помощи разных программ, например, HTSeq[121]. Анализ на уровне транскрипта более сложен и требует привлечения вероятностных методов для оценки численности транскрипта на основании коротких прочтений; например, это может выполнить программа cufflinks[106]. Прочтения, одинаково хорошо ложащиеся на разные места в геноме, должны быть идентифицированы и удалены либо выровнены на одно из возможных мест либо на наиболее вероятное из них. Некоторые методы оценки вовсе не предполагают выравнивания прочтения на референсный геном. Например, метод, используемый в программе kallisto, объединяет псевдовыравнивание и количественный анализ в один шаг, который идёт на два порядка быстрее, чем методы программ tophat и cufflinks, и требует меньших вычислительных затрат[122].

Дифференциальная экспрессия

Когда для каждого транскрипта получены количественные данные, с помощью их статистического анализа, моделирования и нормализации анализируется дифференциальная экспрессия генов[102]. Большинство программ, которые её анализируют, берут на вход таблицу из названий генов и числа транскриптов для каждого из них, но некоторые программы, например, cuffdiff, в качестве входных данных получают выравнивание прочтений в Шаблон:Нп5 (от Шаблон:Lang-en — карта попарных выравниваний). На выходе программы выдают список генов с результатами попарных статистических тестов, проверяющих значимость различий в экспрессии между экспериментальными и контрольными данными[123].

Программы для анализа дифференциальной экспрессии генов в экспериментах РНК-Seq
Программа Среда Специализация
Cuffdiff2[102] Основанные на Unix Анализ транскриптов, направленный на обнаружение событий альтернативного сплайсинга мРНК
EdgeR[88] R/Bioconductor Любые количественные геномные данные
DEseq2[124] R/Bioconductor Разные типы данных
Limma/Voom[87] R/Bioconductor Microarray or RNA-Seq data, flexible experiment design
Ballgown[125] R/Bioconductor Эффективное и чувствительное отыскание транскриптов

Подтверждение

Результаты транскриптомного анализа можно подтвердить с помощью других методов, например, количественной ПЦР (qPCR)[126]. Экспрессия генов измеряется относительно стандартной экспрессии изучаемого гена и контрольных генов. Принцип измерения в qPCR такой же, что и в РНК-Seq, а именно, значение для данного гена рассчитывается на основе концентрации участка-мишени в исследуемом образце. Однако qPCR подходит только для Шаблон:Нп5, имеющих менее 300 п. о. и расположенных вблизи 3'-конца кодирующей области[127]. Если необходимо проверить данные по изоформам транскрипта, с помощью тщательного анализа выравниваний прочтений РНК-Seq можно определить, каким участкам должны соответствовать праймеры для qPCR, чтобы сделать различия наиболее явными[128]. Измерение экспрессии контрольных генов наряду с исследуемыми даёт стабильные референсные данные. Проверка данных РНК-Seq с помощью контрольной ПЦР показала, что разные варианты РНК-Seq в целом дают схожие данные[58][129][130].

Для анализа транскриптомных данных очень важна информация о функциях изучаемых генов. Наблюдаемые паттерны экспрессии генов можно связать с определённым фенотипом при помощи экспериментов по нокдауну изучаемых генов и Шаблон:Нп5[131].

Применение

Диагностика и профилирование заболеваний

Транскриптомные технологии нашли применение в различных областях биомедицины, в частности, в диагностике и профилировании заболеваний[9]. С помощью РНК-Seq стало возможным обнаружение сайтов начала транскрипции, использования альтернативных промоторов и новых вариантов альтернативного сплайсинга. Поскольку геномные Шаблон:Нп5 играют важную роль в патогенезе многих заболеваний, определение их вариантов чрезвычайно важно для интерпретации данных полногеномного поиска ассоциаций[132]. С помощью РНК-Seq можно обнаруживать однонуклеотидные полиморфизмы, связанные с болезнями, случаи аллелеспецифичной экспрессии, Шаблон:Нп5, которые позволяют пролить свет на генетические основы развития заболеваний[133].

С помощью РНК-Seq можно получить информацию о транскрипции эндогенных ретротранспозонов, которые могут влиять на транскрипцию соседних генов посредством разнообразных эпигенетических механизмов, что может привести к развитию болезней[134]. Важной потенциальной сферой применения РНК-Seq является исследование молекулярных основ нарушений иммунной системы, потому что этот метод позволяет разделить популяции иммунных клеток разных типов и секвенировать репертуары T- и B-клеточных рецепторов пациентов[135][136].

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

С помощью РНК-Seq можно находить изменения в экспрессии генов у человеческих патогенов, что может помочь идентифицировать новые Шаблон:Нп5, предсказывать устойчивость к антибиотикам, а также понять детали Шаблон:Нп5[137][138]. С помощью РНК-Seq можно разработать оптимизированные меры Шаблон:Нп5, а также направленные индивидуальные стратегии лечения[136].

Транскриптомный анализ можно проводить и для хозяина, и для патогена. С помощью двойного РНК-Seq можно одновременно построить профили экспрессии генов и хозяина, и патогена в ходе всего инфекционного процесса. Такой подход позволяет изучать динамический иммунный ответ и межвидовые Шаблон:Нп5 для обоих взаимодействующих организмов с момента начального контакта до инвазии и финальной персистенции патогена или его разрушения иммунной системой хозяина[139][140].

Ответы на условия окружающей среды

Транскриптомика позволяет идентифицировать гены и метаболические пути, ответственные за ответ и противодействие стрессам, связанным с биотическими и абиотическими факторами внешней среды[141][131]. Благодаря неспецифичным методам транскриптомики с её помощью можно находить новые генные сети даже в сложных системах. Например, сравнительный анализ нескольких линий нута на различных стадиях развития позволил идентифицировать транскрипционные профили, связанные со стрессами, вызванными засухой и повышенной солёностью; в частности, была показана роль изоформ транскриптов Шаблон:Нп5-Шаблон:Нп5[141]. Изучение экспрессии генов при формировании биоплёнок патогенными дрожжами Candida albicans позволило выявить набор совместно регулируемых генов, критически важных для образования и поддержания биоплёнки[142].

Транскриптомное профилирование даёт ценнейшую информацию о механизмах лекарственной устойчивости. Анализ более тысячи Шаблон:Нп5 малярийного плазмодия Plasmodium falciparum[143] показал, что устойчивость к артемизинину изолятов из Юго-Восточной Азии связана с повышенной активностью Шаблон:Нп5 и с более медленным прохождением внутриэритроцитарной стадии жизненного цикла[144].

Аннотация функций генов

Одно из применений транскриптомных технологий заключается в определении функций генов, а также аллелей, ответственных за конкретный фенотип. Транскриптомика экотипов растения Arabidopsis, которые Шаблон:Нп5 металлы, показала связь с этим фенотипом генов, отвечающих за проникновение металлов в организм, толерантность и гомеостаз[145]. Объединение данных РНК-Seq, полученных для разных тканей, позволило усовершенствовать аннотацию функций генов у коммерчески важных организмов, например, огурца[146], или видов, находящихся под угрозой исчезновения, таких как коала[147].

Сборка прочтений РНК-Seq не зависит от референсного генома[114], поэтому этот метод идеально подходит для изучения экспрессии генов у немодельных организмов, для которых ещё нет готовых геномных данных. Например, база данных однонуклеотидных полиморфизмов, которая использовалась в программах по размножению псевдотсуги Мензиса, была создана при транскриптомном анализе de novo при отсутствии секвенированного генома[148]. Схожим образом гены, участвующие в развитии сердечной, мышечной и нервной тканей у омара, были идентифицированы путём сравнения транскриптомов разных тканей без использования последовательности генома. РНК-Seq также можно использовать для обнаружения ранее неизвестных белоккодирующих областей в уже секвенированных геномах[149].

Некодирующие РНК

Обычно транскриптомика рассматривает только мРНК клетки. Однако те же методы можно применить и к некодирующим РНК, которые участвуют в трансляции, репликации геномной ДНК, сплайсинге и регуляции транскрипции[150][151][152][153]. Многие из этих некодирующих тРНК связаны с развитием болезней, в том числе рака, сердечно-сосудистых заболеваний и Шаблон:Нп5[154].

Базы данных транскриптомов

При изучении транскриптомов создаются огромные объёмы данных, которые потенциально могут быть использованы и в других проектах. Поэтому сырые или обработанные данные размещают в публичных базах данных, чтобы сделать их доступными для всего научного сообщества. Например, по состоянию на 2018 год, база данных Gene Expression Omnibus содержит данные миллионов экспериментов[155].

Базы данных транскриптомов
Название Владелец Данные Описание
Gene Expression Omnibus[95] NCBI Микрочипы, РНК-Seq Первая база данных транскриптомов, полученных из разных источников. Первая ввела стандарты Шаблон:Нп5 и MINSEQE, которые регламентируют необходимые метаданные для эксперимента, чтобы он был хорошо интерпретируем и воспроизводим[156][157].
ArrayExpress[158] Шаблон:Нп5 Микрочипы Импортирует наборы данных из Gene Expression Omnibus и подчиняется ей. Обработанные данные и метаданные экспериментов хранятся в ArrayExpress, а сырые прочтения — в ENA. Соответствует стандартам MIAME и MINSEQE[156][157].
Шаблон:Нп5[159] EBI Микрочипы, РНК-Seq Содержит данные о тканеспецифичной экспрессии генов у животных и растений. Содержит данные вторичных анализов и их визуализацию, использует термины Gene Ontology, домены Шаблон:Нп5 и метаболические пути. Содержит ссылки на данные по многочисленности белков, если они имеются.
Шаблон:Нп5[160] Частное курирование Микрочипы, РНК-Seq Содержит справочные пояснения к публично доступным транскриптомным данным, в основном касающиеся медицины и биологии растений. Данные отдельных экспериментов нормализованы, что позволяет сравнивать экспрессию генов в различных экспериментах. Для полного доступа необходимо приобрести лицензию, бесплатно доступна лишь часть базы.
RefEx[161] DDBJ Все Транскриптомы, полученные из 40 различных органов человека, мыши и крысы. Данные по экспрессии генов визуализированы в виде тепловой карты, наложенной на трёхмерную модель анатомической структуры.
NONCODE[162] noncode.org РНК-Seq Некодирующие РНК (кроме тРНК и рРНК)

Примечания

Шаблон:Примечания

Внешние ссылки

  1. Шаблон:Cite web
  2. 2,0 2,1 Шаблон:Cite pmid
  3. Шаблон:Cite pmid
  4. 4,0 4,1 4,2 Шаблон:Cite pmid
  5. Шаблон:Cite pmid
  6. 6,0 6,1 Шаблон:Cite pmid
  7. Шаблон:Cite pmid
  8. 8,0 8,1 8,2 8,3 8,4 8,5 Шаблон:Cite pmid
  9. 9,00 9,01 9,02 9,03 9,04 9,05 9,06 9,07 9,08 9,09 9,10 Шаблон:Cite pmid
  10. 10,0 10,1 10,2 Шаблон:Cite pmid
  11. 11,0 11,1 11,2 Шаблон:Cite pmid
  12. Шаблон:Cite pmid
  13. Шаблон:Cite pmid
  14. 14,0 14,1 14,2 Шаблон:Cite pmid
  15. Шаблон:Cite pmid
  16. Шаблон:Cite pmid
  17. Шаблон:Cite pmid
  18. Шаблон:Cite pmid
  19. 19,0 19,1 19,2 Шаблон:Cite pmid
  20. Шаблон:Cite pmid
  21. 21,0 21,1 21,2 21,3 21,4 21,5 Шаблон:Cite pmid
  22. Шаблон:Cite pmid
  23. 23,0 23,1 Шаблон:Cite pmid
  24. Шаблон:Cite pmid
  25. 25,0 25,1 25,2 25,3 25,4 25,5 Шаблон:Cite web
  26. 26,0 26,1 Шаблон:Cite pmid
  27. Шаблон:Cite pmid
  28. Шаблон:Cite pmid
  29. Шаблон:Cite pmid
  30. Шаблон:Cite pmid
  31. 31,0 31,1 Шаблон:Cite pmid
  32. Шаблон:Cite pmid
  33. 33,0 33,1 Шаблон:Cite pmid
  34. 34,0 34,1 34,2 Шаблон:Cite pmid
  35. Шаблон:Книга
  36. Шаблон:Cite pmid
  37. Шаблон:Cite pmid
  38. 38,0 38,1 Шаблон:Cite pmid
  39. Шаблон:Cite pmid
  40. Шаблон:Cite pmid
  41. 41,0 41,1 Шаблон:Cite pmid
  42. 42,0 42,1 Шаблон:Cite pmid
  43. Шаблон:Cite pmid
  44. Шаблон:Cite pmid
  45. Шаблон:Cite pmid
  46. Шаблон:Cite pmid
  47. 47,0 47,1 47,2 47,3 47,4 Шаблон:Cite pmid
  48. 48,0 48,1 Шаблон:Cite pmid
  49. Шаблон:Cite pmid
  50. 50,0 50,1 Шаблон:Cite pmid
  51. Шаблон:Cite pmid
  52. Шаблон:Cite pmid
  53. Шаблон:Cite pmid
  54. Шаблон:Cite pmid
  55. Шаблон:Cite pmid
  56. Шаблон:Cite pmid
  57. Шаблон:Cite doi
  58. 58,0 58,1 Шаблон:Cite pmid
  59. Шаблон:Cite pmid
  60. Шаблон:Cite pmid
  61. 61,0 61,1 Шаблон:Cite pmid
  62. Шаблон:Cite pmid
  63. 63,0 63,1 Шаблон:Cite pmid
  64. Шаблон:Cite pmid
  65. Шаблон:Cite pmid
  66. Шаблон:Cite pmid
  67. Шаблон:Cite pmid
  68. Шаблон:Cite pmid
  69. Шаблон:Cite pmid
  70. Шаблон:Cite pmid
  71. Шаблон:Cite pmid
  72. Шаблон:Cite pmid
  73. Шаблон:Cite pmid
  74. 74,0 74,1 Шаблон:Cite pmid
  75. Шаблон:Cite web
  76. Шаблон:Cite pmid
  77. Шаблон:Cite pmid
  78. Шаблон:Cite pmid
  79. Шаблон:Cite pmid
  80. Шаблон:Cite pmid
  81. 81,0 81,1 Шаблон:Cite pmid
  82. 82,0 82,1 82,2 Шаблон:Cite pmid
  83. Шаблон:Cite pmid
  84. Шаблон:Cite pmid
  85. Шаблон:Cite pmid
  86. Шаблон:Cite web
  87. 87,0 87,1 Шаблон:Cite pmid
  88. 88,0 88,1 Шаблон:Cite pmid
  89. 89,0 89,1 Шаблон:Cite pmid
  90. Шаблон:Книга
  91. Шаблон:Книга
  92. 92,0 92,1 Шаблон:Cite pmid
  93. 93,0 93,1 Шаблон:Cite pmid
  94. Шаблон:Cite pmid
  95. 95,0 95,1 Шаблон:Cite pmid
  96. 96,0 96,1 Шаблон:Cite doi
  97. Шаблон:Книга
  98. Шаблон:Cite pmid
  99. Шаблон:Cite pmid
  100. Шаблон:Cite web
  101. Шаблон:Cite pmid
  102. 102,0 102,1 102,2 Шаблон:Cite pmid
  103. 103,0 103,1 Шаблон:Cite pmid
  104. Шаблон:Cite pmid
  105. Шаблон:Cite pmid
  106. 106,0 106,1 Шаблон:Cite pmid
  107. Шаблон:Cite pmid
  108. Шаблон:Cite pmid
  109. Шаблон:Cite pmid
  110. Шаблон:Cite pmid
  111. Шаблон:Cite pmid
  112. Шаблон:Cite pmid
  113. Шаблон:Cite pmid
  114. 114,0 114,1 Шаблон:Cite pmid
  115. Шаблон:Cite pmid
  116. Шаблон:Cite pmid
  117. Шаблон:Cite pmid
  118. Шаблон:Cite pmid
  119. Шаблон:Cite pmid
  120. Шаблон:Cite pmid
  121. Шаблон:Cite pmid
  122. Шаблон:Cite pmid
  123. Шаблон:Cite pmid
  124. Шаблон:Cite pmid
  125. Шаблон:Cite pmid
  126. Шаблон:Cite pmid
  127. Шаблон:Cite pmid
  128. Шаблон:Cite pmid
  129. Шаблон:Cite pmid
  130. Шаблон:Cite pmid
  131. 131,0 131,1 Шаблон:Cite pmid
  132. Шаблон:Cite pmid
  133. Шаблон:Cite pmid
  134. Шаблон:Cite pmid
  135. Шаблон:Cite pmid
  136. 136,0 136,1 Шаблон:Cite pmid
  137. Шаблон:Cite pmid
  138. Шаблон:Cite pmid
  139. Шаблон:Cite pmid
  140. Шаблон:Cite pmid
  141. 141,0 141,1 Шаблон:Cite pmid
  142. Шаблон:Cite pmid
  143. Шаблон:Cite pmid
  144. Шаблон:Cite pmid
  145. Шаблон:Cite pmid
  146. Шаблон:Cite pmid
  147. Шаблон:Cite pmid
  148. Шаблон:Cite pmid
  149. Шаблон:Cite pmid
  150. Шаблон:Cite pmid
  151. Шаблон:Cite pmid
  152. Шаблон:Cite pmid
  153. Шаблон:Cite pmid
  154. Шаблон:Cite pmid
  155. Шаблон:Cite web
  156. 156,0 156,1 Шаблон:Cite pmid
  157. 157,0 157,1 Шаблон:Cite pmid
  158. Шаблон:Cite pmid
  159. Шаблон:Cite pmid
  160. Шаблон:Cite pmid
  161. Шаблон:Cite pmid
  162. Шаблон:Cite pmid

Шаблон:Выбор языка Шаблон:Избранная статья