В настоящее время актуальным является изучение не только фенотипических, но и генетических свойств промышленно-перспективных штаммов лактобацилл с использованием различных молекулярно-генетических методов, в том числе основанных на секвенировании нуклеиновых кислот.
Современные технологии массивного параллельного секвенирования (next-generation sequencing – NGS) позволяют получать информацию о полногеномных последовательностях широкого спектра бактерий и вирусов, и проводить молекулярно-генетические исследования на качественно новом уровне [1].
Цель исследования: изучение структуры генома штамма Lactobacillus fermentum 630 с использованием метода полногеномного секвенирования (whole-genome shotgun sequencing).
Материалы и методы.
Объектом исследования являлся штамм L. fermentum 630 из Государственной коллекции лактобацилл Нижегородского НИИЭМ им. академика И.Н. Блохиной, выделенный из кишечника здорового человека.
Видовая принадлежность исследуемого штамма подтверждена путем изучения его биохимических свойств с использованием тест-систем API 50 CHL (Biomerueux, Франция) и масс-профилей рибосомальных белков, полученных с помощью времяпролетного MALDI масс-спектрометра Autoflex (Bruker Daltonics, Германия). Клетки исследуемого микроорганизма имеют вид неподвижных полиморфных грамположительных палочек, не образующих спор и капсул.
Выделение и очистку ДНК исследуемого штамма проводили с использованием набора «Ампли-Сенс ДНКсорб В» (ЦНИИЭ, Москва). Концентрацию ДНК в образцах определяли с помощью флуориметра Qubit (Invitrogen, Австрия). Подготовку библиотеки ДНК для секвенирования осуществляли с использованием набора Nextera XT (Illumina, США) на 500 циклов согласно инструкции производителя. Секвенирование генома проводили на платформе NGS MiSeq (Illumina, США) в режиме ресеквенирования. В качестве референса служила полногеномная последовательность штамма L. fermentum F6 (номер GenBank NC_021235.1). Для выравнивания нуклеотидных последовательностей использовали программу Burrows-Wheeler Aligner (BWA), для поиска однонуклеотидных полиморфизмов, инсерций и делеций – программное обеспечение Genome Analysis Toolkit (GATK). Визуализацию и анализ полученных данных проводили с помощью программного обеспечения UGENE [6]. Аннотацию генома проводили с использованием сервера RAST (Rapid Annotation using Subsystem Technology) [7]. Филогенетическое дерево строили с помощью программы Mega 6.0 [8] методом ближайших соседей (neighbor-joining), предварительно проведя множественное выравнивание нуклеотидных последовательностей 19 геномов L. fermentum, задепонированных в базе данных NCBI (National Center for Biotechnology Information) и L. fermentum 630 с помощью сервера RALPHY 1.10 (Reference sequence Alignment based Phylogeny builder) [5].
Результаты и обсуждение
В результате секвенирования было получено 2 237 480 коротких чтений длиной до 250 н.о. Вторичный анализ позволил получить 390 контигов, которые были объединены в 250 скаффолдов. Общая длина последовательности составила 1,830,322 пар нуклеотидов, средний уровень покрытия - 54,6. С помощью программного обеспечения GATK было обнаружено 15076 однонуклеотидных полиморфизмов, 213 инсерций и 253 делеции относительно референсного штамма L. fermentum F6.
При аннотировании с помощью сервера RAST в геноме было выявлено 1668 белок кодирующих последовательностей, 58 последовательностей тРНК и 19 рРНК.
В базе данных NCBI представлены нуклеотидные последовательности геномов 15 штаммов L. fermentum в виде контигов и скаффолдов (на 12 ноября 2015 года). В таблицедана сравнительная характеристика структур генома L. fermentum 630 и неполных геномов штаммов L. fermentum, задепонированных в базе данных Whole Genome Shotgun contigs (WGS) NCBI.
Сравнительная характеристика структур геномов L. fermentum 630 и задепонированных в базе данных WGS NCBI
Штамм L. fermentum | Размер, Mb | % ГЦ | Число белок кодирующих последовательностей | Кол-во тРНК | Кол-во рРНК |
630 | 1,83032 | 52,45 | 1668 | 58 | 19 |
ATCC 14931 | 1,86701 | 52,60 | 1656 | 56 | 4 |
28-3-CHN | 2,02652 | 52,00 | 1806 | 53 | 4 |
FTDC8312 | 1,96655 | 51,80 | 1718 | 16 | 2 |
NB-22 | 2,01131 | 51,80 | 1826 | 31 | 4 |
LfQi6 | 2,1998 | 51,50 | 1847 | 56 | 8 |
DSM 20055 | 1,90005 | 52,4 | 1713 | 47 | 3 |
Lf1 | 1,81565 | 52,50 | 1633 | 57 | 15 |
MTCC8711 | 2,5663 | 46,65 | 1911 | 57 | 15 |
39 | 1,82966 | 51,60 | 1658 | 51 | 13 |
90 TC-4 | 1,82248 | 51,90 | 1571 | 54 | 11 |
L930BB | 1,97184 | 52,10 | 1767 | 56 | 7 |
779_LFER | 1,93581 | 52,10 | 1729 | 49 | 15 |
HFB3 | 2,04336 | 51,80 | 1382 | 58 | 20 |
UCO-979C | 2,01183 | 51,90 | 1492 | 57 | 13 |
222 | 1,95041 | 52,10 | 1864 | 54 | 4 |
Следует отметить, что данные, полученные в нашем исследовании, совпадают с результатами аналогичных проектов, представленных в таблице. Высокопроизводительное секвенирование генома штамма L. fermentum 630 позволило определить 80-90% нуклеотидной последовательности.
При проведении филогенетического анализа исследуемого штамма с последовательностями геномов L. fermentum, представленных в базе данных NCBI, было установлено что штамм L. fermentum 630 находится на одной ветке со штаммом L. fermentum 39. Филогенетически близкородственными являются также штаммы L. fermentum NB-22 и L. fermentum UCO-979C (рисунок).
Филогенетическое дерево штаммов L. fermentum, построенное на основании анализа нуклеотидных последовательностей геномов методом ближайших соседей. Цифрами указана статистическая значимость порядка ветвления (в %), определенная с помощью бутстрап-анализа 1000 альтернативных деревьев
В структуре генома с помощью сервера RAST обнаружены нуклеотидные последовательности 21 мобильного элемента, 2 генов, кодирующих фактор элонгации G, обеспечивающего устойчивость к тетрациклину, и 2 генов, кодирующих белки семейства бета-лактамаз, гена хлорамфеникол ацетилтрансферазы, детерминирующего устойчивость к хлорамфениколу, а также гены белков семейства фаговых интеграз.
Выводы:
- С использованием высокопроизводительного секвенирования на секвенаторе MiSeq получена молекулярно-генетическая характеристика штамма L. fermentum 630.
- В геноме L. fermentum 630 выявлено значительное количество однонуклеотидных полиморфизмов, а также инсерций и делеций относительно референсного штамма L. fermentum F6
- Наиболее близкородственным, по результатам филогенетического анализа является штамм L. fermentum 39.
- Обнаруженные последовательности мобильных элементов, фаговых интеграз, генетических детерминант антибиотикорезистентности послужат основой для создания генетической карты генома и обнаружения генетических детерминант, обусловливающих уникальные фенотипические свойства.
Рецензенты:
Смирнов В.Ф., д.б.н., профессор кафедры физиологии и биохимии человека и животных Федерального государственного автономного образовательного учреждения высшего образования «Национальный исследовательский государственный Нижегородский университет им. Н.И. Лобачевского» Министерства образования России, г.Нижний Новгород;
Заславская М.И., д.б.н., профессор кафедры микробиологии и иммунологии Государственного бюджетного образовательного учреждения высшего профессионального образования «Нижегородская государственная медицинская академия» Министерства образования и социального развития России, г.Нижний Новгород.
Библиографическая ссылка
Алексеева А.Е., Точилина А.Г., Бруснигина Н.Ф., Солнцев Л.А., Белова И.В., Ефимов Е.И., Соловьева И.В. ПОЛНОГЕНОМНOE СЕКВЕНИРОВАНИE ГЕНОМА ШТАММА LACTOBACILLUS FERMENTUM 630 // Современные проблемы науки и образования. – 2015. – № 6. ;URL: https://science-education.ru/ru/article/view?id=24029 (дата обращения: 10.02.2025).