Изучение генетического разнообразия вируса мозаики турнепса и создание удвоенных гаплоидных линий brassica – источников устойчивости к нему
На правах рукописи
Зубарева Ирина Александровна ИЗУЧЕНИЕ ГЕНЕТИЧЕСКОГО РАЗНООБРАЗИЯ ВИРУСА МОЗАИКИ ТУРНЕПСА И СОЗДАНИЕ УДВОЕННЫХ ГАПЛОИДНЫХ ЛИНИЙ BRASSICA – ИСТОЧНИКОВ УСТОЙЧИВОСТИ К НЕМУ Специальности: 03.01.06 – биотехнология (в том числе бионанотехнологии), 03.02.07 – генетика
АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата биологических наук
Москва - 2013
Работа выполнена в Федеральном государственном бюджетном учреждении науки Центре «Биоинженерия» Российской академии наук
Научный консультант:
доктор биологических наук Игнатов Александр Николаевич
Официальные оппоненты:
Карлов Геннадий Ильич, доктор биологических наук, профессор, Российский государственный аграрный университет МСХА имени К.А.
Тимирязева, руководитель научно-образовательного Центра молекулярной биотехнологии Супрунова Татьяна Павловна, кандидат сельскохозяйственных наук, Всероссийский научно-исследовательский институт селекции и семеноводства овощных культур, лаборатория биотехнологии, старший научный сотрудник
Ведущая организация: Всероссийский научно-исследовательский институт овощеводства Российской академии сельскохозяйственных наук
Защита состоится 23 октября 2013 г. в 16-00 часов на заседании диссертационного совета Д220.043.10 при Российском государственном аграрном университете – МСХА имени К.А.Тимирязева по адресу 127550, г. Москва, ул. Прянишникова, д. 15, тел/факс (499) 976-24-92, e-mail:
С диссертацией можно ознакомиться в Центральной научной библиотеке им. Н.И. Железнова РГАУ-МСХА имени К.А. Тимирязева.
Автореферат разослан 23 сентября 2013 г.
Ученый секретарь диссертационного совета Л.С. Большакова
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность проблемы Вирус мозаики турнепса (Turnip mosaic virus – TuMV) относится к крупнейшей группе фитопатогенных вирусов семейства Potyviridae (Shukla et al., 1994). Это наиболее известный представитель рода Роtyvirus, поражающий растения семейства Brassicaceae. При поражении растений происходит задержка в росте, ранняя дефолиация, снижение товарного вида кочанных и салатных форм, что приводит к заметному снижению урожайности (до 100 % в зависимости от культуры) и экономическим потерям. Во многих странах TuMV по экономической вредоносности стоит на втором месте, уступая лишь вирусу мозаики огурца (Cucumber Mosaic Virus – CMV) (Walsh et al., 2002).
Вредоносность TuMV резко усиливается при совместном поражении растений с другими вирусами (Camele et al., 1991).
Один из наиболее опасных способов распространения и сохранения вирусной инфекции – передача вируса с семенами, что имеет серьезные последствия и создает очаги первичной инфекции при посеве/посадке растений в поле. Вирусы могут сохраняться в семени в течение длительного времени и распространяться на большие расстояния (Matthews, 2009). Естественная устойчивость растений является единственно возможным, эффективным и экологически безопасным способом борьбы с TuMV.
В настоящее время в России лишь отдельные исследования связаны с TuMV (Микрюков, 2010) и нет данных о его патотипическом составе, притом что опасность распространения инфекции нарастает, а количество сортов Brassica, устойчивых к патогену, невелико. Известные сорта эффективны только против какого-либо одного патотипа TuMV, а гены устойчивости малоизученны. В качестве объекта исследований при изучении генов устойчивости наиболее часто используют удвоенные гаплоидные линии (Rusholme et al., 2007;
Jenner et al., 2010;
Qian et al., 2013). Поэтому для изучения генов устойчивости создание дигаплоидных линий (DH-линий), одновременно являющихся источниками генов устойчивости к TuMV, представляет собой важную научную и практическую задачу.
Цель и задачи работы Целью настоящих исследований было изучение генетического разнообразия TuMV и создание удвоенных гаплоидных линий рода Brassica – источников устойчивости к патогену.
Для достижения поставленной цели решали следующие задачи:
1. Провести молекулярно-генетический и филогенетический анализ изучаемых изолятов TuMV.
2. Оценить коллекции Brassica на устойчивость к TuMV и получить удвоенные гаплоидные линии – источники устойчивости к нему.
3. Оценить полиморфизм локусов eIF4E и eIF(iso)4E, связанных с устойчивостью к TuMV, у представителей Brassica.
4. Изучить возможность передачи изолята TuMV I2 через семена Brassica, получить и проанализировать полную нуклеотидную последовательность кодирующей части его генома.
Научная новизна работы Идентифицированы 6 новых изолятов TuMV, принадлежащих к двум различным патотипам, один из которых впервые обнаружен в России.
Впервые проведен анализ наличия локусов эукариотических фактров инициации трансляции у различных видов Brassica: B. napus, B. juncea, B. nigra, B. oleracea и B. carinata. В результате сравнения аминокислотных последовательностей локусов BraA.eIF(iso)4E.а и BraA.eIF4E.а выявлены характерные для всех устойчивых образцов замены, вероятно, ведущие к неспецифической устойчивости образцов к TuMV. Кроме того, показано наличие трех аллельных вариантов BraA.eIF4E.а у некоторых устойчивых к TuMV образцов B. rapa var. oleifera, один из которых оказался псевдогеном.
Впервые доказана возможность передачи TuMV (изолят I2) семенами растений рода Brassica через заражение зародыша семени.
Получена полная последовательность кодирующей части генома изолята TuMV I2, проведен ее биоинформатический анализ, выявлены уникальные аминокислотные замены, вероятно, позволяющие изоляту передаваться семенами растений рода Brassica.
Практическая значимость Получены удвоенные гаплоидные линии 13 генотипов видов B. napus и B. oleracea, обладающих устойчивостью к 6 изолятам TuMV 1 и 5 патотипов.
Описаны локусы генов eIF(iso)4E и eIF4E, представляющие собой новый и потенциально прочный источник устойчивости к TuMV, который может быть передан восприимчивым сортам и культурам рода Brassica.
Подобраны праймеры и проба для разработки набора диагностики TuMV в семенах и инфицированных растениях методом ПЦР в реальном времени.
Все нуклеотидные последовательности были опубликованы в международной базе данных GenBank и находятся в свободном доступе.
Апробация работы Результаты работы были представлены на 2-ой Международной школе конференции молодых ученых «Генетика и селекция растений, основанная на современных генетических знаниях и технологиях» (Звенигород, 2012);
XXIV и XXV зимних молодежных научных школах «Перспективы направления физико-химической биологии и биотехнологии» (Москва, 2012 и 2013);
II(X) Международной Ботанической Конференции молодых ученых (Санкт Петербург, 2012);
Осеннем финале «У.М.Н.И.К.» РАН (Москва, 2012);
Всероссийской Научной Конференции «Научное Наследие Н.И. Вавилова и Современность»;
I Международной Интернет-конференции «Современные тенденции в сельском хозяйстве» (Казань, 2012);
Международной научно практической конференции «Фармацевтические биотехнологии» (Москва, 2013);
XII Молодежной научной конференции «Биотехнология в растениеводстве, животноводстве и ветеринарии» (Москва, 2013);
17-ой Международной Пущинской школе-конференции молодых ученых (Пущино, 2013).
Личный вклад автора заключается в проведении экспериментальных и теоретических исследований. Основные результаты работы получены лично автором при его непосредственном участии в планировании и проведении экспериментов. Имена соавторов указаны в соответствующих публикациях.
Публикации По материалам диссертации опубликовано 15 научных работ, в том числе 6 статей в журналах, входящих в перечень ВАК.
Объем и структура диссертации. Материалы диссертации изложены на 144 страницах машинописного текста и включают 34 рисунка и 24 таблицы.
Диссертация состоит из введения, обзора литературы, описания материалов и методов исследования, результатов и их обсуждения, выводов, списка литературы из 230 источников и 1 приложения.
РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЯ И ИХ ОБСУЖДЕНИЕ Определение TuMV с использованием растений-индикаторов и ИФА Растительный материал с симптомами вирусной мозаики собирали в четырех районах Москвы и Московской области (Москва, Одинцовский р-н, Подольский р-н, Наро-Фоминский р-н) в ходе обследований посевов и посадок овощных и масличных капустных культур в открытом и закрытом грунте.
Всего было собрано 87 образцов, которые классифицировали по признакам поражения и растениям-хозяевам на 11 групп (I1, I2, I3a, I3b, I4, I5, I6, I7, I8, I и I10).
Результаты заражения растений-индикаторов (Nicotiana tabacum L., Chenopodium amaranticolor Coste, Lactuca sativa L., Brassica rapa subsp rapa) и ИФА (DAS-ELISA) с использованием поликлональных антител к TuMV в составе тест-системы NEOGEN показали, что не более 15 % из 87 образцов были инфицированы возбудителем мозаики турнепса. Отобранные изоляты были получены из пораженных растений B. oleracea (I2), B. rapa subsp rapa (I3a), B. napus (I3b) из закрытого грунта и из посадок B. rapa subsp chinensis в поле (I7, I8, I10).
Для выделения и очистки изолятов TuMV проводили серию пассажей, используя C. amaranticolor или C. quinoa, а затем через N. tabacum изоляты переносили и поддерживали на B. rapa или B. napus.
Молекулярно-генетический и филогенетический анализ изолятов TuMV Для молекулярно-генетического анализа из изолятов и контрольного образца (растение B. napus, проверенное ИФА на отсутствие TuMV) выделили тотальную РНК и с помощью ОТ-ПЦР с использованием специфических пар праймеров, амплифицировали и секвенировали фрагменты кДНК генов Р1, NIb и CP. Области секвенирования охватывали три «горячие точки» рекомбинации в геноме вируса. Нуклеотидные последовательности участков трех генов сравнили с последовательностями базы данных GenBank при помощи программы BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi), в результате чего выделили последовательности четырех наиболее близких изолятов TuMV (UK1, KEN1, RUS1 и CAR51) cо сходством 97-99 % к отдельным изучаемым изолятам и построили филогенетические деревья как для отдельных генов, так и для суммарной последовательности кДНК (рисунок 1). Секвенирование участков трех генов показало, что изучаемые изоляты отличались не только от всех известных мировых штаммов TuMV, но и между собой.
Рисунок 1 – Филогенетическое дерево, построенное с использованием алгоритма минимальной эволюции на основе сравнения нуклеотидных последовательностей объединенных фрагментов генов P1, NIb и CP изучаемых изолятов TuMV общей длинной 1650 п.н. На шкале представлено генетическое расстояние между последовательностями кДНК, определенное методом наибольшего правдоподобия Анализ деревьев, полученных методами минимальной эволюции, ближайших соседей и максимального правдоподобия, показал высокое сходство их топологии и продемонстрировал для последовательностей генов NIb, CP и P1 наибольшую близость группы изолятов I2, I3a и I3b к изолятам, выделенным из растений B. napus и B. oleracea в Великобритании и относящимся к 1 патотипу (UK1, KEN1). В то же время группа изолятов I7, I8 и I10 была близка к последовательностям изолятов, выделенных из Armoracia rusticana и Cochlearia armoracia в России и Польше и принадлежащих к патотипу TuMV (RUS1, CAR51). Генетическое сходство между двумя группами изолятов было менее 95 %, что сопоставимо с ранее описанным сходством групп изолятов вируса внутри вида в пределах 85-95 % (Rzhetsky et al., 1992).
Сравнение нуклеотидных и аминокислотных последовательностей участков трех генов изучаемых изолятов с наиболее близкими к ним изолятами TuMV UK1 и RUS1 показало, что наиболее полиморфным является фрагмент гена P1, содержащий от 1,6 до 3,0 % замен пар нуклеотидов, из которых почти половина ведет к заменам аминокислот. Наиболее консервативным оказался участок гена СР, имеющий всего 0,9-1,3 % нуклеотидных замен, практически не ведущих к заменам аминокислот в сравнении с UK1 и RUS1.
Рассмотрение последовательностей с помощью методов, реализованных в программе RDP, не выявило фактов рекомбинации в изучаемых изолятах по последовательности трех анализируемых генов (NIb, CP, P1). Таким образом, популяции вируса, поражающие растения капусты, рапса и турнепса в защищенном грунте и растения пекинской капусты в полевых условиях, принадлежат к разным патотипам и имеют независимое происхождение.
Анализ вирулентности изолятов TuMV Для анализа вирулентности и агрессивности изолятов TuMV, а также с целью отбора образцов для создания линий – источников устойчивости к вирусу, проводили механическую инокуляцию 136 образцов растений семейства Brassicaceae из двух коллекций. Коллекция № 1 состояла из образцов и включала в себя представителей B. napus, B. oleracea, B. carinata, B. juncea и B. nigra. Коллекция № 2 включала 72 образца, представляющих стержневую коллекцию Brassica rapa ВИРа.
Симптомы заболевания, такие как хлоротичная пятнистость и деформация листовой пластинки, оценивали визуально по 10-бальным шкалам с семи дневным интервалом в течение четырех недель после инокуляции, после чего со всех образцов собирали материал для проведения ИФА. Результаты, полученные ИФА, подтвердили результаты визуальной оценки и выявили дополнительно три толерантных (бессимптомных) к TuMV образца растений.
По результатам оценки устойчивости коллекций растений рода Brassica было выявлено, что изоляты TuMV отличаются между собой по вирулентности (таблица 1). Наиболее вирулентными оказались изоляты I3b и I2, поражающие 63 и 51 образец из 136 соответственно. Однако изолят I2 был умеренно агрессивным и поражал растения на уровне 0-5 баллов. Наименее вирулентный изолят I10 поражал всего 13 образцов.
Кроме того, характер проявления симптомов мозаики листьев при инфицировании растений изолятами I7, I8 и I10 отличался от симптомов, наблюдаемых при поражении изолятами I2, I3a и I3b.
Таблица 1 – Оценка вирулентности изолятов TuMV на коллекции из образцов рода Brassica Количество устойчивых образцов Количество восприимчивых образцов, шт.
Изолят (коллекция № 1 и № 2) TuMV коллекция № 1 коллекция № 2 шт. % I2 31 20 85 I3a 23 5 108 I3b 49 14 73 I7 1 28 107 I8 0 31 104 I10 1 12 123 В результате анализа вирулентности изолятов было отобрано 13 образцов (№ 444, Center-1, № 1, № 417, Jp-8, Ханна, Jp-4, Griffin, Сenter-2.1, Галант, Луговской, Cobra SR B. napus и ISA454 B. oleracea) с устойчивостью ко всем шести изолятам TuMV и еще 12 с устойчивостью к пяти из шести изолятов (Oro, Jp-2, Cobra Winter, Сenter -2.2, № 417.1.1, Jp-1, Griffin DR, Griffin CR, Сobra linc 1.8, STS К0 B. napus, FH1 и ISA454.1.36 B. oleracea). Данные образцы являются наиболее перспективными в селекции на устойчивость к TuMV.
Устойчивые ко всем шести изучаемым изолятам TuMV образцы Brassica в дальнейшем использовали для получения удвоенных гаплоидных линий – источников устойчивости к TuMV и изучения генов устойчивости к данному вирусу.
Получение удвоенных гаплоидных линий Brassica – источников устойчивости к TuMV Гаплоидные растения получали с применением изолированной культуры микроспор. Для выделения микроспор у B. napus использовали бутоны размером от 2,0 до 3,5 мм в зависимости от сорта, а у B. oleracea образец ISA – 5,0-6,0 мм. Температурный диапазон теплового шока, при котором были получены эмбриоиды, лежал в пределах от +32 до +35 °С, а продолжительность воздействия варьировала от 24 до 72 часов. Оптимальное для эмбриогенеза каждого образца сочетание температуры и продолжительности теплового шока приведено в таблице 2.
Таблица 2 – Оптимальное сочетание температуры и продолжительности теплового шока для эмбриогенеза у исследуемых образцов Brassica Образец Образец Параметры Параметры (количество (количество теплового шока теплового шока эмбриоидов, шт.) эмбриоидов, шт.) № 1 (20 шт.) № 417 (74 шт.) + 32 оС, 48 ч + 33 оС, 72 ч ISA454 (60 шт.) Луговской (9 шт.) Галант (12 шт.) Griffin (38 шт.) № 444 (32 шт.) Cobra SR (4 шт.) Jp-8 (24 шт.) о о + 32 С, 72 ч + 35 С, 24 ч Center-1 (53 шт.) Хана (5 шт.) Center-2.1 (6 шт.) Jp-4 (25 шт.) После тепловой обработки культуру микроспор переносили на +25 оС, и из отдельных микроспор формировались проэмбриональные структуры, которые развивались в эмбриоиды, дающие начало гаплоидным растениям.
Через 21 день образовавшиеся эмбриоиды переносили на твердую питательную среду MS (Murashige and Skoog, 1962) или B5 (Gamborg, 1968), приготовленную с уменьшением концентрации компонентов в 2 раза, где проходило формирование семядольных листьев и корневой системы.
Полученные растения-регенеранты адаптировали к почвенным условиям на стадии формирования 3-4 пары листьев. При оценке плоидности полученных из микроспор растений использовали прямой (подсчет количества хромосом в меристеме корня) и косвенный (подсчет числа хлоропластов в замыкающих клетках устьиц) методы. Гаплоидные растения отбирали и удваивали количество их хромосом впрыскиванием раствора колхицина в верхушечную и пазушные почки с целью получения дигаплоидных линий (таблица 3).
Таблица 3 Количество полученных удвоенных гаплоидных линий образцов Brassica Количество Количество Образец Brassica дигаплоидных Образец Brassica дигаплоидных линий, шт. линий, шт.
Луговской ISA454 50 № 444 Галант 22 Jp-8 20 Griffin Ханна 3 Cobra SR Jp-4 18 Center-1 №1 15 Center-2.1 № 417 Семена удвоенных гаплоидов вновь высевали для подтверждения устойчивости к TuMV. При достижении фазы 3-4 настоящих листьев проводили механическую инокуляцию шестью изолятами TuMV, по 5- растений потомства каждой дигаплоидной линии каждого образца. Симптомы заболевания оценивали визуально, а затем в растениях определяли содержание TuMV с помощью ИФА.
По результатам опыта отбирали по одной наилучшей устойчивой ко всем шести изолятам TuMV линии от каждого образца. Таким образом, было получено 13 удвоенных гаплоидных линий: DH № 444, DH Center-1, DH № 1, DH № 417, DH Jp-8, DH Ханна, DH Jp-4, DH Griffin, DH Сenter-2.1, DH Галант, DH Луговской, DH Cobra SR и DH ISA454.
Анализ полиморфизма локусов eIF(iso)4E и eIF4E у представителей Brassica При анализе вирулентности и агрессивности изолятов TuMV были отобраны устойчивые ко всем шести изолятам образцы Brassica: представителей B. napus, 1 B. oleracea и 22 B. rapa. Поскольку изучаемые изоляты относятся к двум различным патотипам (1 и 5 патотип) и в то же время не известно ни одного гена устойчивости к 5 патотипу, то было сделано предположение, что устойчивость образцов является неспецифической.
Неспецифическую устойчивость растений к TuMV связывают с факторами инициации трансляции eIF(iso)4E и eIF4E (Rusholme et al., 2007;
Jenner et al., 2010;
Qian et al., 2013). По литературным данным, в геноме В. rapa ssp.
trilocularis было обнаружено 6 таких локусов (BraA.eIF(iso)4E.a, и BraA.eIF(iso)4E.b, BraA.eIF(iso)4E.c, BraA.eIF4E.a, BraA.eIF4E.b BraA.eIF4E.c.), один из которых (BraA.eIF4E.b) был псевдогеном (Jenner et al., 2010). Наличие или отсутствие данных локусов в геномах остальных представителей Brassica до сих пор не было показано.
Для обнаружения и оценки полиморфизма данных локусов, а также установления их отличий по нуклеотидным последовательностям между устойчивыми и восприимчивыми к TuMV образцами использовали различных представителей Brassica: устойчивые к TuMV 12 DH-линии B. napus, 1 DH линия B. oleracea, 22 образца B. rapa;
по 1 образцу восприимчивых к TuMV B. juncea, B. nigra и B. carinata, B. napus (Галакси) и B. rapa (Kousaitai).
В результате проведенных исследований наличие локусов BraA.eIF(iso)4E.a и BraA.eIF(iso)4E.c с помощью ПЦР показано во всех образцах B. napus, B. oleracea (рисунок 2), B. nigra, B. carinata и в 11 из образцах (в т.ч. восприимчивый образец) B. rapa.
Локус BraA.eIF4E.a был обнаружен у всех представителей Brassica (в т.ч.
B. juncea) и исследуемых устойчивых к TuMV образцов. Кроме того, у образцов B. rapa наблюдалась его вариация по длине и наличие у некоторых образцов до 2-3 аллельных вариантов. Локусы BraA.eIF(iso)4E.b, и BraA.eIF4E.b BraA.eIF4E.c не были обнаружены ни у одного образца.
Рисунок 2 – Электрофореграмма продуктов ПЦР: а – BraA.eIF(iso)4E.a и BraA.eIF(iso)4E.c;
б – BraA.eIF4E.a. 1 – маркер молекулярной массы (Fermentas, Gene Ruler Dna Laders);
2 – ПЦР без ДНК-матрицы;
3 – ПЦР на препарате ДНК восприимчивого к TuMV образца B. napus Галакси;
4-16 - ПЦР на препаратах ДНК устойчивых к TuMV DH-линий:
№ 444, Center-1, № 1, № 417, Jp-8, Ханна, Jp-4, Center-2.1, Галант, Луговской, Griffin, Cobra SR и ISA 454 соответственно Для установления факта наличия или отсутствия локусов BraA.eIF(iso)4E.a и BraA.eIF(iso)4E.c в геноме 12 из 23 образцов B. rapa, не дающих продукт амплификации, проводили дот-блот анализ препаратов геномной ДНК. Продукты ПЦР локусов BraA.eIF(iso)4E.a и BraA.eIF(iso)4E.c метили с использованием дигоксигенин-11-дУТФ. Дот-блот гибридизация показала наличие искомых локусов в образцах B. rapa (рисунок 3). Вероятно, нуклеотидные последовательности локусов сильно изменены в местах посадки праймеров, что и приводит к отсутствию продукта амплификации.
Рисунок 3 – Дот-блот гибридизация ДНК образцов B. rapa с зондами на локусы: а – BraA.eIF(iso)4E.a;
б – BraA.eIF(iso)4E.c. +К – ДНК восприимчивого к TuMV образца Kousaitai (B. rapa);
R2–R67 – ДНК образцов B. rapa, не дающих продукты амплификации;
-К – ДНК плазмиды pBR328;
0 – ТЕ-буфер без ДНК-матрицы Продукты ПЦР локусов и BraA.eIF4E.a, BraA.eIF(iso)4E.a BraA.eIF(iso)4E.c пяти восприимчивых к TuMV: Галакси (B. napus), Kousaitai (B. rapa), горчица сарептская краснолистная (B. juncea), № 12 (B. nigra) и № (B. carinata);
и восьми устойчивых образцов: DH № 417 (B. napus), DH Jp- (B. napus), DH ISA454 (B. oleracea), IT 103369 (B. rapa), cурепица Palton sarson и К-374 88/47 (B. rapa), капуста китайская Dark green leaf и Местный К- (B. rapa);
секвенировали. Каждый ген состоит из пяти экзонов и четырех интронов (Jenner et al., 2010). После выравнивания нуклеотидные последовательности сравнивали между собой и с последовательностями, взятыми из GenBank, на наличие полиморфизма и отличий в нуклеотидных последовательностях между устойчивыми и восприимчивыми к TuMV образцами.
Полиморфизм белков оценивали по показателям энтропии (Hx). Наиболее полиморфным оказался локус BraA.eIF4E.а, в котором замены аминокислот присутствуют во всех пяти экзонах. Аллельные варианты данного локуса 1 и различались между собой длиной 1 интрона, а 3 аллельный вариант (образцы B. rapa var.oleifera R2, R4, R19, R40, R41 и R61) оказался псевдогеном, т.к. в нем отсутствовала часть 3 экзона, полностью 3 интрон и 4 экзон. Меньше полиморфизма наблюдалось в BraA.eIF(iso)4E.c и BraA.eIF(iso)4E.а, где аминокислотные замены встречались в 1, 2 и 3 экзонах и ни одной замены в четвертом и пятом.
В результате выравнивания аминокислотных последовательностей локуса BraA.eIF(iso)4E.а была выявлена замена аминокислоты 152Gly Asp (глицин на аспарагиновую кислоту), характерная для всех устойчивых образцов (рисунок 4). Данная аминокислотная замена также была обнаружена в B. rapa ssp. chinensis и отмечалась, как мутация, ведущая к устойчивости к 4 патотипу TuMV (Qian et al., 2013).
Рисунок 4 Выравнивание последовательностей аминокислот в коротком сегменте экзона 3 BraA.eIF(iso)4E.а среди пяти восприимчивых и восьми устойчивых к TuMV образцов Brassica. 80425 – аминокислотная последовательность локуса BraA.eIF(iso)4E.а восприимчивой к TuMV линии B. rapa ssp. chinensis Замены в локусе BraA.eIF(iso)4E.c преимущественно были локализованы в третьем экзоне, который, возможно, принимает участие в связывании с VPg белком TuMV, но не у всех устойчивых образцов.
В связи с тем, что локус BraA.eIF4E.а оказался наиболее полиморфным и присутствовал у всех представителей Brassica, секвенировали дополнительные образцы. В результате выравнивания аминокислотных последовательностей выявлена замена 40Thr Ile (треонин на изолейцин) у всех устойчивых к TuMV образцов (рисунок 5). У образца К-131 Tyre Индия данная замена наблюдалась в обоих аллельных вариантах локуса, а у К-13 Аргентина в одном. Однако нефункциональная для репликации TuMV аллель данного локуса имеет доминантный характер, что и позволяет образцам проявлять устойчивую фенотипическую реакцию к TuMV.
Рисунок 5 – Выравнивание последовательностей аминокислот в коротком сегменте экзона 1 BraA.eIF4E.а среди шести восприимчивых и девятнадцати устойчивых к TuMV образцов Brassica. eIF4E.a – аминокислотная последовательность локуса BraA.eIF4E.а восприимчивой к TuMV линии B. rapa ssp. trilocularis R-o- Неспецифическая устойчивость у Brassica обеспечивается двумя взаимодействующими генами retr01 и ConTR01 (Rusholme et al., 2007).
Рецессивная аллель гена retr01, совпадающая с локусом фактора инициации трансляции BraA.eIF(iso)4E.a не обеспечивает репликацию TuMV и усиливает проявление устойчивости, контролируемой вторым геном ConTR01, который наоборот имеет доминантную аллель, совпадающую с BraA.eIF4E.a и не функциональную для TuMV. Оба данных локуса (BraA.eIF(iso)4E.a и BraA.eIF4E.a) имеют мутации, характерные для устойчивых образцов и, вероятно, являются нефункциональными для присоединения VPg-белка TuMV, тем самым, обеспечивая растению устойчивость к вирусу.
Локусы eIF(iso)4E и eIF4E представляет собой новый и потенциально прочный источник неспецифической устойчивости к TuMV, который может быть передан восприимчивым сортам и культурам рода Brassica.
Анализ возможности передачи TuMV семенами растений рода Brassica Подбор растения-индикатора для диагностики TuMV по методике Холмса Методика Холмса была предложена для диагностики ВТМ на семядолях огурца (Holmes, 1929). По литературным данным Cucumis sativus не инфицируется TuMV. Поэтому индикаторное растение подбирали экспериментально из представителей семейства Cucurbitaceae (Тыквенные):
Cucurbita maxima L. (тыква гигантская), Cucurbita pepo L. (кабачок) Ecballium elaterium L. (бешеный огурец обыкновенный), Cucurbita melo L. (дыня).
Механическую инокуляцию семядольных листьев проводили соком неинфицированного растения B. napus и двумя изолятами, относящимися к разным патотипам TuMV (I2 и I7).
Через 5-7 дней после заражения инокулированные семядольные листья срезали, часть смотрели на наличие крахмальных паранекрозов при помощи йодного окрашивания, а часть проверяли на наличие TuMV ИФА. Результаты ИФА подтвердили результаты визуальной оценки.
Крахмальные паранекрозы на семядольных листьях, как результат механической инокуляции соком растения Brassica с вирусной инфекцией, были отмечены на C. pepo, E. elaterium и C. melo. Семядольные листья E. elaterium и C. melo очень нежные и небольших размеров, что делает работу более трудоемкой. Поэтому в качестве растения-индикатора TuMV по методике Холмса был выбран C. pepo, имеющий оптимальный размер и толщину семядольных листьев.
Диагностика TuMV в семенах, зародышах и взрослых растениях Растения, инфицированные исследуемыми штаммами TuMV, как правило, не достигали фазы цветения или не завязывали семена. Лишь при поражении TuMV I2 с толерантных и слабовосприимчивых образцов удалось собрать семена. Слабовосприимчивые образцы: яровой рапс (B. napus) Cobra Spring (1), Галакси (5), CrGC5.1 (10), № 423 (24), Aomory-3 (30), Aomory-1 (39), Ратник (45);
горчица индийская (B. juncea) Лада (3а), FBLM (8a);
и толерантный к вирусу образец B. napus STS K0 (47) после механической инокуляции штаммом TuMV I2 довели до цветения, самоопылили и собрали семена.
Перед тестированием семян на наличие в них вирусной инфекции проводили поверхностное обеззараживание тремя способами: 70 % этанолом + NaOCl (pH 12,6) – спирт и щелочная среда инактивирует поверхностную вирусную инфекцию;
мирамистином – противовирусный коммерческий препарат;
прогревом семян в стерильной дистиллированной воде при 50 оС в течение 20 мин, а затем перемещение на 5 мин в холодную стерильную дистиллированную воду – термическая обработка семян с целью инактивации вирусных частиц. После предварительного обеззараживания семена были рендомизированно высажены в изолированном боксе теплицы при соблюдении мер предосторожности против случайной контаминации вирусом.
При благоприятных для накопления вирусной инфекции в растениях условиях (+ 24 оС) через 20-25 дней растения шести образцов (1, 10, 24, 39, 45 и 3а) из девяти показали симптомы. При последующем снижении температуры до +18 оС симптомы исчезали в течение двух недель, что подтверждало толерантную и слабовосприимчивую реакцию отобранных растений. Со всех образцов был собран материал для проверки наличия вируса методами ИФА, ОТ-ПЦР, методом Холмса с модификациями.
ИФА и метод инокуляции семядольных листьев C. pepo по Холмсу (рисунок 6) и ОТ-ПЦР на участок гена Р1 (рисунок 7) показали наличие TuMV в восьми образцах (1, 5, 10, 24, 39, 45, 47Т и 3а) из девяти.
1 2 3 Рисунок 6 – Пример крахмальных паранекрозов на семядольных листьях C. pepo, инокулированных: 1 – образцом, не содержащим по результатам оценки TuMV;
2 – соком здорового растения B. napus;
3 –инокулюмом TuMV I2;
4 – образцом, содержащим TuMV, переданный семенами Рисунок 7 – Электрофореграмма продуктов ПЦР с праймерами на ген Р1: 1 – маркер молекулярной массы (100bp Fermentas);
2 –ПЦР без кДНК-матрицы;
3 – ПЦР на препарате кДНК здорового растения;
4 – ПЦР на препарате кДНК TuMV I2;
5-13 – ПЦР на препаратах кДНК исследуемых образцов (1, 5, 10, 24, 39, 45, 47Т, 3а, 8а соответственно) Для обнаружения TuMV в семенах изучаемых образцов пробы по 20 шт.
каждого образца (1, 5, 10, 24, 30, 39, 45, 3а и 8а) были взяты для ИФА и ОТ ПЦР. ОТ-ПЦР показал наличие TuMV в пробах семян для семи из девяти изучаемых образцов. ИФА выявил TuMV только в двух пробах № 423 (24) и Ратник (45), что говорит о большой концентрации вируса в них. Причем, вероятно, накопление вирусной инфекции в семенах находится в обратной корреляции с восприимчивостью к вирусу образцов, так как № 423 и Ратник поражались TuMV I2 на уровне 2-3 баллов.
Для изучения механизма передачи TuMV семенами отбирали пробы по семян каждого из четырех образцов (Aomory-3, Ратник, Лада и FBLM) для индивидуального анализа зародышей семян на присутствие в них TuMV методом ОТ-ПЦР, методикой Холмса с модификациями и ИФА. ОТ-ПЦР диагностировал TuMV только в четырех и трех из десяти зародышей образцов Ратник и Лада соответственно. Методика Холмса с модификациями показала наличие TuMV в трех зародышах из десяти в образце Ратник и в четырех из десяти зародышей образца Лада, а в зародышах образцов Aomory-3 и FBLM вирус не был обнаружен. Таким образом, TuMV в зародышах семян образцов Aomory-3 и FBLM не был обнаружен ни одним из методов. Однако в зародышах семян образцов Ратник и Лада реакция на вирус была положительной для всех используемых методов.
Для идентификации передаваемого семенами вируса из зараженных растений выделили тотальную РНК и с помощью ОТ-ПЦР с использованием специфических пар праймеров амплифицировали и секвенировали фрагмент кДНК гена Р1. Нуклеотидные последовательности были идентичны исходному изоляту TuMV I2, а сравнение их с последовательностями базы данных, выявило наиболее близкий изолят TuMV UK1 (Великобритания), с идентичностью 99 %.
Детекция TuMV методом ПЦР в реальном времени Передача вируса с семенами имеет важное экономическое значение, так как обеспечивает наиболее опасный способ распространения и сохранения вирусной инфекции в течение длительного времени. Для ограничения распространения TuMV на территории РФ необходимо проводить проверку зараженности этим вирусом семян овощных и масличных крестоцветных культур. Поэтому с целью своевременной детекции вируса при малых его концентрациях необходимо применение более чувствительных методов, например, метода ПЦР в реальном времени.
Для проведения ПЦР в реальном времени последовательности праймеров и зонда подбирали на наиболее консервативный и специфичный для TuMV участок гена P1. Для регистрации накопления продуктов ПЦР использовали линейные разрушаемые пробы (TaqMan), позволяющие детектировать накопление фрагмента ДНК строго определенной последовательности.
В результате анализа данных ПЦР в реальном времени была показана эффективность реакции на уровне 1,7-1,8. Пороговое значение флуоресценции для определения величины порогового цикла (Ct) было установлено на уровне 30 во время экспоненциального возрастания кривой амплификации. Реакции со значениями Ct менее 45 и экспоненциальными амплификационными участками считали положительными. Если значения Ct не были получены из-за отсутствия точки пересечения кривой амплификации с пороговым значением, реакции оценивали как отрицательные (таблица 4).
Таблица 4 – Значения порогового цикла (Ct) ПЦР в реальном времени для изолятов TuMV, проб семян и растений, выросших из них Наличие Наличие Образец FAM Сt Образец FAM Сt TuMV TuMV Изоляты TuMV -К (ПЦР без кДНК-матрицы) – 0,00 I3b 18,97 + -К (кДНК здорового растения) – 0,00 I7 18,36 + I2 17,93 + I8 19,42 + I3a 17,96 + I10 16,99 + Пробы семян -К (ПЦР без кДНК-матрицы) – 0,00 24 29,92 + -К (кДНК здорового растения) – 0,00 30 30,39 + +К (кДНК изолята I2 TuMV) 15,23 + 39 23,11 + 1 29,88 + 45 29,29 + 3а 5 29,87 + 21,33 + 8а 10 17,31 + 30,90 + Растения, выросшие из инфицированных семян -К (ПЦР без кДНК-матрицы) – 0,00 24 23,53 + -К (кДНК здорового растения) – 0,00 30 26,17 + +К (кДНК изолята I2 TuMV) 20,73 + 39 24,81 + 1 22,93 + 45 23,26 + 3а 5 24,91 + 21,94 + 8а 10 21,68 + 27,11 + Примечание: FAM Ct величины порогового цикла по флуоресцентному красителю FAM По результатам ПЦР в реальном времени было подтверждено аличие TuMV в образцах. Достижение пороговой величины наблюдали на препаратах кДНК изолятов TuMV на 17-19 циклах, на препаратах кДНК проб семян исследуемых образцов – 22-31 циклах и на препаратах кДНК растений, выросших из инфицированных семян, на 20-27 циклах. В отрицательных контролях (без кДНК матрицы и кДНК неинфицированного растения) экспоненциального возрастания кривой амплификации не наблюдали на протяжении всех 45 циклов. В пробах, приготовленных из семян и выросших из них растений, образцов Aomory-3 (30) и FBLM (8а) по данным ПЦР в реальном времени было показано наличие TuMV с величиной порогового около 30 и соответственно, что несколько меньше, чем в других образцах.
Полученные данные позволяют рекомендовать подобранные праймеры и пробу для разработки набора для диагностики TuMV методом ПЦР в реальном времени в семенах и инфицированных растениях.
Генетический и филогенетический анализ изолята I2 TuMV По литературным данным TuMV не передается семенами, однако в данной работе показана передача TuMV I2 с семенами растений рода Brassica.
С целью выявления мутаций, которые позволили изоляту TuMV I передаваться через семена, была определена полная кодирующая последовательность его генома (GenBank № KC297103).
Возможные генетические связи и филогенетическая группировка полных последовательностей геномов всех известных изолятов TuMV при исследовании изучали тремя методами: максимального правдоподобия, минимальной эволюции и ближайших соседей. Филогенетические деревья рассчитывали для геномных последовательностей 112 изолятов TuMV, включая все рекомбинантные штаммы. По патогенности и происхождению все изоляты TuMV делят на четыре группы: базальная-B, базальная-BR, азиатская-BR и мировая-B. Изолят I2 принадлежит к мировой-В группе, и наиболее близкими к нему изолятами являются UK1 (AF169561.2), GBR50 (AB252114.1), GBR (EU861593.1) и GBR36 (AB252113.1).
Таблица 5 Показатели изменчивости генов для мировой-B группы TuMV Гены TuMV Показатель P1 HC-Pro P3 6K1 CI 6K2 VPg Pro NIb CP GC состав, % 48,56 45,56 44,69 45,63 45,01 46,79 47,15 44,59 45,38 47, Ts / Tv 8,99 11,29 6,55 5,95 9,73 12,23 9,86 12,13 15,99 5, N/S 0,254 0,036 0,229 0,016 0,039 0,057 0,024 0,025 0,048 0, Число полиморфных 495 451 398 42 550 49 201 207 425 сайтов Нуклеотидное 0,079 0,058 0,069 0,049 0,041 0,040 0,056 0,047 0,034 0, разнообразие -1,37 -1,34 -1,22 -1,24 -1,76 -1,82 -1,48 -1,35 -1,92 -1, Tajimas D Р0,1 Р0,1 Р0,1 Р0,1 Р0,1 Р0,05 Р0,1 Р0,1 Р0,05 Р0, Примечание: Ts – транзиция;
Tv – трансверсия;
S – синонимичная замена аминокислоты;
N – несинонимичная замена аминокислоты По показателям изменчивости генов для мировой-В группы видно, что GC состав всех 10 генов TuMV в среднем находится на одном уровне (таблица 5). Отношение Тs/Тv далеко от единицы, что говорит об отсутствии повторных замен нуклеотидов (реверсий) в генах TuMV. Число несинонимичных замен меньше числа синонимичных, а их отношение отлично от 1, следовательно, в мировой-В группе TuMV имеет место стабилизирующий отбор. Нуклеотидное разнообразие генов показывает, что наиболее полиморфными генами TuMV являются P1 и P3, а самыми стабильными – CP и NIb. Тест Tajimas D говорит об отборе, т.к. его значения находятся в пределах -2…0, то отбора нет и существенной адаптации генов не выявлено.
Для детекции возможных рекомбинационных событий в изоляте TuMV I использовали методы, реализованные в программе RDP. Результаты показали, что изолят TuMV I2 имеет два рекомбинационных события, которые были идентичны с изолятом UK1. TuMV I2 и UK1 являются рекомбинантами между FKD001 (95,7 %) и Tu-3 (96,3 %), принадлежащими к мировой-В группе.
Нуклеотиды 1-328 генома I2 были из изолята FKD001J (R. sativus, Япония), а 465-9367 Tu-3 (B. oleracea, Япония). Сайты рекомбинации были расположены в районе середины гена P1 и 3'-участка гена CP. Рекомбинационные события были обнаружены по меньшей мере шестью методами (RDP, GENECONV, MaxChi, Chimaera, SiScan и 3Seq) и имели P-значения от 6,204 10-05 до 2,516 10-85.
Поскольку изоляты TuMV I2 и UK1 имеют общие рекомбинационные события, то, вероятно, изолят I2 произошел от UK1 в результате мутаций.
Результаты оценки времени наиболее вероятного расхождения генотипов I2 и UK1, полученные с помощью эволюционного анализа с использованием модели р-расстояния (Nei and Kumar, 2000), согласуется с информацией об эпифитотии заболевания в конце 1990-х годов после испытания коллекции сортов пекинской капусты, завезенной в Россию из Великобритании, на Селекционной станции РГАУ-МСХА имени К.А. Тимирязева (Г.Ф. Монахос, личное сообщение).
При сравнении нуклеотидной последовательности кодирующей части генома TuMV I2 с наиболее близкими изолятами (UK1, GBR50, GBR98, GBR36) было обнаружено 99 нуклеотидных замен (92 транзиции (40 С/U и G/A) и 7 трансверсий), распределенных по всем 10 генам вируса (таблица 7).
Таблица 7 – Количество и частота нуклеотидных замен в геноме TuMV I в сравнении с наиболее близкими изолятами (UK1, GBR50, GBR98, GBR36) Ген Количество замен Частота нуклеотидных замен, % название п.н. нуклеотидов аминокислот P1 1086 12 4 1, – HC-Pro 1374 6 0, – P3 1065 8 0, – 6K1 156 1 0, – CI 1932 33 1, – 6K2 159 3 1, – VPg 576 2 0, – Pro 729 4 0, – NIb 1551 21 1, CP 862 9 1 1, Пять замен были несинонимичными и уникальными в сравнении со всеми известными штаммами TuMV и привели к изменениям в двух из десяти закодированных вирусных белков. Четыре из них были в P1 и одна в СР (рисунок 8). Несмотря на большое количество нуклеотидных замен в геноме TuMV I2 в сравнении с наиболее близкими к нему изолятами, в среднем GC состав остался на уровне изолята UK1 и составил 46,38 %.
Рисунок 8 – Схематическое изображение структуры генома с TuMV I расположением уникальных аминокислотных замен Белок P1 относится к числу наиболее полиморфных потивирусных белков, участвует в репликации вируса и вместе с HC-Pro действует как усилитель патогенности, подавляя сайленсинг генов. Функции СР заключаются в раздевании и амплификации вируса, его транспорте от клетки к клетке, по всему растению и в передаче тлями (Urcuqui-Inchima et al., 2001). Вероятно, обнаруженные мутации усиливают амплификацию генома и транспорт вируса по растению, одновременно ослабляя его агрессивность, что позволяет растению достигать фазы цветения и завязывать семена, в которых накапливается вирусная инфекция.
ВЫВОДЫ 1. В результате молекулярно-генетического анализа были определены новых изолятов TuMV, принадлежащих к 1 и 5 патотипам вируса, имеющих различное происхождение, вирулентность и характер заражения растений.
2. Получены удвоенные гаплоидные линии 13 генотипов видов B. napus и B. oleracea, обладающих устойчивостью к 6 изолятам TuMV, относящихся к 1 и 5 патотипам.
3. Впервые показано наличие локуса эукариотического фактора инициации трансляции BraA.eIF4E.a у изученных представителей всех видов рода Brassica и отсутствие его изоформ BraA.eIF(iso)4E.a и BraA.eIF(iso)4E.c в геноме B. juncea.
4. В результате сравнения аминокислотных последовательностей локусов BraA.eIF(iso)4E.а и BraA.eIF4E.а у устойчивых к TuMV образцов Brassica выявлены аминокислотные замены 152Gly Asp и 40Thr Ile, соответственно.
5. Впервые показана возможность передачи TuMV (изолят I2) семенами растений рода Brassica через зародыш семени.
6. В результате анализа полной последовательности кодирующей части генома изолята TuMV I2 показано, что способность передаваться семенами может быть связана с уникальными аминокислотными заменами в последовательностях генов P1 (122, 134, 261, 301) и CP (2896).
Список работ, опубликованных по теме диссертации:
Зубарева, И.А. Вирус мозаики турнепса передается семенами / И.А.
1.
Зубарева, А.Н. Игнатов, Т.Н. Грибова, С.Г. Монахос // Защита и карантин растений. Москва. – 2013. – № 1. – С. 37-40.
Зубарева, И.А. Молекулярная характеристика изолятов Turnip 2.
mosaic virus из растений рода Brassica / И.А. Зубарева, Т.Н. Грибова, С.Г.
Монахос, А.Н. Игнатов // Известия ТСХА. Москва. – 2013. – № 1. – С. 29-35.
Зубарева, И.А. Передача штамма I2 Turnip mosaic virus c семенами 3.
растений рода Brassica / И.А. Зубарева, А.Н. Игнатов, Т.Н. Грибова, С.Г.
Монахос // Известия ТСХА. Москва. – 2013. – № 2. – С. 34-41.
Зубарева, И.А. Генетическое разнообразие вируса мозаики турнепса 4.
и механизм его передачи семенами растений рода Brassica / И.А. Зубарева, С.В.
Виноградова, Т.Н. Грибова, С.Г. Монахос, академик К.Г. Скрябин, А.Н.
Игнатов // Доклады академии наук. Москва. – 2013. – Том 450, № 1. – С. 105 108.
Зубарева, И.А. Создание дигаплоидных линий Brassica napus – 5.
доноров устойчивости к вирусу мозаики турнепса / И.А. Зубарева, Е.Н.
Головешкина, С.В. Виноградова, Т.Н. Грибова, С.Г. Монахос, А.Н. Игнатов // Сельскохозяйственная биология. Москва. – 2013. – № 5. – С. 122-125.
Зубарева, И.А. Бактериальные и вирусные болезни 6.
сельскохозяйственных культур: распространение и диагностика / А.Н. Игнатов, С.В. Виноградова, Е.Н. Головешкина, И.А. Зубарева // Овощи России. – 2013. – № 2. – С. 67-68.
Зубарева, И.А. Изучение коллекции исходного материала растений 7.
рода Brassica как источника генов устойчивости к вирусу мозаики турнепса (TuMV) / И.А. Зубарева, Т.Н. Грибова, А.Н. Игнатов // Тезисы докладов II международной школы-конференции молодых ученых "Генетика и селекция растений, основанная на современных генетических знаниях и технологиях".
Москва: Звенигород. – 2011. – С. 41.
Зубарева, И.А. Физиологическая и молекулярная характеристика 8.
изолятов вируса мозаики турнепса (TuMV) из растений рода Brassica / И.А.
Зубарева, Т.Н. Грибова, А.Н. Игнатов // Тезисы докладов и стендовых сообщений XXIV зимней молодежной научной школы «Перспективные направления физико-химической биологии и биотехнологии». Москва. – 2012.
– С. 90.
Зубарева, И.А. Передача изолята I2 Turnip mosaic virus с семенами 9.
растений рода Brassica / И.А. Зубарева, Т.Н. Грибова, А.Н. Игнатов // Тезисы докладов II (X) Международной ботанической конференции молодцх ученых в Санкт-Петербурге. Санкт-Петербург. – 2012. – С. 62.
10. Зубарева, И.А. Создание дигаплоидных линий масличного рапса – доноров устойчивости к вирусу мозаики турнепса / И.А. Зубарева, Т.Н.
Грибова, А.Н. Игнатов // Сборник тезисов докладов Осеннего финала по программе «У.М.Н.И.К.» РАН. Москва, Президиум РАН. – 2012. – С. 16.
11. Зубарева, И.А. Механизм передачи изолята I2 Turnip mosaic virus с семенами растений рода Brassica / И.А. Зубарева, С.В. Виноградова, Т.Н.
Грибова, А.Н. Игнатов // сборник трудов I международной интернет конференции «Современные тенденции в сельском хозяйстве». Казань. – 2012.
– С. 96-97.
12. Зубарева, И.А. Создание устойчивых к TuMV удвоенных гаплоидных линий Brassica / И.А. Зубарева, Е.Н. Головешкина, Т.Н. Грибова, А.Н. Игнатов // сборник трудов I международной интернет-конференции «Современные тенденции в сельском хозяйстве». Казань. – 2012. – С. 98-99.
13. Зубарева, И.А. Генетическое разнообразие TuMV и создание удвоенных гаплоидных линий – доноров устойчивости к нему / И.А. Зубарева, Т.Н. Грибова, А.Н. Игнатов // Тезисы докладов и стендовых сообщений XXV международной зимней молодежной научной школы «Перспективные направления физико-химической биологии и биотехнологии». Том 2. Конкурс молодых ученых. Москва. – 2013. – С. 49.
14. Зубарева, И.А. Молекулярная характеристика изолятов TuMV и создание дигаплоидных линий Brassica – доноров устойчивости к нему / И.А.
Зубарева, А.Н. Игнатов // XII Молодежная научная конференция «Биотехнология в растениеводстве, животноводстве и ветеринарии». – 2013. – С. 24-25.
15. Зубарева, И.А. Молекулярная характеристика изолятов TuMV и создание дигаплоидных линий Brassica napus – доноров устойчивости к нему / И.А. Зубарева, Т.Н. Грибова, А.Н. Игнатов // 17 Международная Пущинская школа-конференция молодых ученых. – 2013. – С. 335-336.