Клинические штаммы E. coli
Для иллюстрации использования ОНТ для секвенирования плазмид из клинических штаммов авторы секвенировали 10 штаммов — щедрый подарок доктора Стивена Салипанте из Университета Вашингтона. Эти штаммы были частью коллекции из 312 изолятов экстракишечного патогенного кишечного кишечника (ExPEC), полученных в ходе планового клинического лечения в Медицинском центре УниверситетаВашингтона 45.
Конъюгация и характеристика плазмид с помощью гелевого теста на мобильность
Исследователи также провели конъюгационный анализ с использованием штамма E. coli LMB100 в качестве реципиента. Конъюгация — это процесс, который приводит к переносу плазмиды, присутствующей в одной клетке (доноре), в другую клетку (называемую реципиентом)46. Перенос требует физического контакта между двумя клетками и обычно обнаруживается с помощью маркеров, указывающих на наличие конъюгативной плазмиды в клетках реципиента. Это позволило авторам исключить различия в штаммах хозяина как переменной и проиллюстрировать использование протокола для секвенирования плазмид, полученных с помощью захвата плазмид (метод профилирования плазмид из определённой экологической ниши). Все конъюгационные анализы производили транссоюганты, подтверждая наличие конъюгативной плазмиды у доноров. Протокол, которому следовали авторы, описан встатье 36. Группа идентифицировала плазмиды, присутствующие в донорских и реципиентных штаммах, и оценила их подвижность с помощью пульсно-полевой гель-электрофорез (PFGE). Список из десяти репрезентативных донорских штаммов приведён в таблице 8 , а предполагаемые размеры плазмид, присутствующих в них, приведены в третьем столбце. Плазмиды за пределами диапазона разрешения геля PFGE исключены из этой таблицы, но исследователям удалось собрать несколько меньших плазмид около ~4 кб (не показано).
| Донорский штамм | Плазмида | Оценочный размер PFGE | Длина сборки | Средняя глубина | bp Zero Depth | Процент покрытия | Изменения позиций | Процент точности до SR |
blood_08 _0081 | pblood_08_ 0081_1 | 170000 | 171829 | 100 | 99.941803 | 54 | 0.0314 | 99.9686 |
blood_08 _0081 | pBlood_08_ 0081_2 | 82000 | 85576 | 46.9 | 0 | 100 | 2 | 99.9977 |
blood_d-08 _0094 | pBlood_d-08 _0094 | 179000 | 173334 | 382.6 | 3121 | 98.19943 | 2,347 | 98.646 |
blood_08 _1447 | pBlood_08 _1447 | 97000 | 99624 | 36.2 | 0 | 100 | 0 | 100 |
blood_10 _0913 | pBlood_10 _0913_1 | 70000 | 71788 | 49 | 4 | 99.994428 | 188 | 99.7381 |
blood_10 _0913 | pBlood_10 _0913_2 | 112000 | 118727 | 41.2 | 126 | 99.893874 | 276 | 99.7675 |
blood_1 1_184 | pBlood_11 _184_1 | 150000 | 146583 | 65.1 | 287 | 99.804206 | 134 | 99.9086 |
| blood_11_184 | pBlood_11 _184_2 | 75000 | 74263 | 82.1 | 0 | 100 | 32 | 99.9569 |
blood_2011 _0238 | pBlood_ 2011_0238 | 170000 | 167286 | 110.9 | 0 | 100 | 0 | 100 |
| upec_108 | pUPEC_108 | 130500 | 131039 | 103.9 | 0 | 100 | 0 | 100 |
| upec_134 | pUPEC_134 | 160500 | 157829 | 169 | 0 | 100 | 0 | 100 |
| upec_271 | pUPEC_271 | 145500 | 144519 | 86 | 33311 | 76.950436 | 1,557 | 98.9226 |
| upec_90 | pUPEC_90 | 70000 | 74624 | 139.8 | 0 | 100 | 6 | 99.992 |
Таблица 8: Оценки размера плазмид и метрики точности секвенирования для панели из 10 донорских штаммов. Перечислены оценки размера плазмид, основанные на PFGE или собранных донорских плазмидных последовательностях ONT. Попарные сравнения плазмидных сборок только с ONT и SR-последовательностей также приводятся как индикатор точности. Метрики включают: длину сборки, среднюю глубину, области с нулевой глубиной, процент покрытия, изменённые позиции и процент точности до SR.
Выходные метрики секвенирования плазмид на основе ONT
Последовательность была запланирована на 24 часа. Всего было мультиплексировано 24 образца. Технические данные по каждой выборке, включая показатель качества (Q), длину и охват, приведены в таблице 9. Секвенирование нанопор обеспечивало достаточный выход считывания и качество для всех анализируемых образцов. Во всех библиотеках секвенирование генерировало в среднем 134 115 чтений на выборку, что соответствует среднему выходу 576,2 Мб. Средняя длина чтений составляла 8 797 bp со средним N50 3 693 bp, обеспечивая достаточную непрерывность чтения для сборки плазмид. Средний показатель качества чтения составлял Q16.0 с медианой Q16.6, а индивидуальные забеги варьировались от Q13.7 до Q18.0.
| Секвенированная деформация | Чтения (kb) | Базы (Mb) | Средняя длина | Максимальная длина | Длина мин. | N50 | Медианная длина | Средний Q | Медиана Q |
| Blood_08_0081 донор | 109.8 | 261 | 5713 | 508418 | 78 | 2377 | 1351 | 14.5 | 15 |
| Blood_08_0094 донор | 89.3 | 229 | 4177 | 114948 | 81 | 2565 | 1065 | 13.7 | 14 |
| Blood_08_1447 донор | 87.5 | 246 | 5041 | 135147 | 75 | 2816 | 1609 | 15.1 | 15.7 |
| Bblood_10_0913 донор | 247.7 | 1026 | 9260 | 174231 | 70 | 4143 | 1661 | 15.5 | 16.1 |
| Blood_11_154 донор | 51.5 | 110 | 4172 | 91148 | 73 | 2127 | 1387 | 16.1 | 16.5 |
| Blood_2011_0238 донор | 16.8 | 42 | 6707 | 104848 | 62 | 2534 | 845 | 15.3 | 15.2 |
| UPEC_108 донор | 50.6 | 77 | 1915 | 90234 | 67 | 1525 | 1103 | 15.8 | 16.2 |
| UPEC_134 донор | 50 | 143 | 7110 | 158393 | 72 | 2865 | 1133 | 16 | 16.2 |
| UPEC_217 донор | 64 | 236 | 7512 | 144275 | 78 | 3678 | 1895 | 17.4 | 17.9 |
| UPEC_90 донор | 52.9 | 157 | 5666 | 95944 | 72 | 2969 | 1593 | 16.8 | 17.2 |
| Blood_08_0081 LMB100 трансконъюгант | 53.5 | 119 | 5162 | 225856 | 1 | 2232 | 1009 | 16 | 16.8 |
| Blood_08_0094 LMB100 трансконъюгант | 160 | 442 | 7257 | 252176 | 1 | 2761 | 1083 | 17.3 | 18.2 |
| Blood_08_1447 LMB100 трансконъюгант | 238 | 1029 | 10249 | 209221 | 1 | 4320 | 1857 | 17.7 | 18.7 |
| Blood_10_0913 LMB100 трансконъюгант | 792.3 | 4341 | 11508 | 474325 | 1 | 5479 | 2731 | 17.7 | 18.8 |
| Blood_11_184 LMB100 трансконъюгант | 42.5 | 281 | 17931 | 146626 | 87 | 6615 | 2448 | 16 | 16.7 |
| Blood_2011_0238 LMB100 трансконъюгант | 147.7 | 585 | 9000 | 344203 | 1 | 3962 | 1781 | 18 | 18.9 |
| UPEC_108 LMB100 трансконъюгант | 173.4 | 778 | 10127 | 168756 | 80 | 4485 | 2079 | 15.3 | 16.1 |
| UPEC_134 Трансконъюгантный LMB100 | 105.4 | 667 | 17568 | 349230 | 71 | 6323 | 2232 | 15.5 | 16.3 |
| UPEC_271 LMB100 трансконъюгант | 74.2 | 336 | 14521 | 342195 | 83 | 4520 | 1306 | 14.6 | 14.8 |
| UPEC_90 LMB100 трансконъюгант | 75.2 | 418 | 15348 | 225661 | 83 | 5561 | 2016 | 15.7 | 16.4 |
Таблица 9: Отчёт NanoQ о выходных метриках секвенирования для панели из 10 донорских штамм. Показатели выхода секвенирования перечислены. К ним относятся значения чтения, оснований, длины (среднего, медианного, максимального и минимального), длин, N50 и Q (среднее и медианное) от доноров и транссопряжентов.
Сборка плазмид и точность
Плазмидные последовательности собирались с помощьюAutocycler 34. Размеры собранных плазмидных последовательностей для 10 образцов приведены в таблице 8, четвёртый столбец, и во всех случаях они совпадали с оценками, основанными на мобильности PFGE (перечисленные в третьем столбце). Точность оценивалась путём сравнения сборок длинного чтения, отполированных Medaka, с сборками Medaka (версия 2.2.0) плюс полиполь (v0.6.1). Идентификация последовательности была рассчитана с использованием dnadiff из MUMmer (v3.2.4) (Таблица 8). Во всех анализируемых плазмидах (n = 12) сборки имели среднюю длину 120,4 кб и поддерживали среднюю глубину секвенирования 109x, при среднем покрытии 97,9% по последовательностям плазмид. До полировки сборки имели среднюю идентичность последовательности 99,74% по сравнению с полированной сборкой SR-последовательности.
Оценка загрязнения хромосом
Для оценки покрытия плазмидами и загрязнения хромосом риды были сопоставлены с собранными плазмидными последовательностями и с хромосомной последовательностью E. coli LMB100. Чтения Nanopore с базовыми названиями выравнивались с помощью minimap2 (v2.2.6). Полученные выравнивания обрабатывались с помощью SAMtools (v1.1.7) для генерации отсортированных файлов BAM и расчёта статистики отображения. Покрытие плазмид и средняя глубина оценивались с помощью покрытия SAMtools. Глубина секвенирования плазмид варьировалась между образцами, но оставалась значительно выше минимального покрытия, необходимого для точной сборки (50–70 x). Небольшая часть ридов отображалась на хромосому, что указывает на перенос ДНК хозяина во время экстракции плазмиды (Таблица 10).
| Плазмиды | Общее количество прочитанных | Хромосомные чтения | % хромосомных показаний | Среднее покрытие хромосом | Среднее покрытие плазмидами | Соотношение P:C |
| pBlood_08_0081_2 | 53491 | 20116 | 37.61% | 5 | 1068 | 236 |
| pBlood_d-08_0094 | 160024 | 74827 | 46.76% | 37 | 1618 | 44 |
| pBlood_08_1447 | 238031 | 127967 | 53.76% | 62 | 7842 | 127 |
| pBlood_10_0913_1 | 792263 | 221610 | 27.97% | 87 | 7801 | 90 |
| pBlood_10_0913_2 | 792263 | 221610 | 27.97% | 87 | 7772 | 89 |
| pBlood_11_184_1 | 42494 | 15919 | 37.46% | 7 | 1637 | 250 |
| pBlood_11_184_2 | 42494 | 15919 | 37.46% | 7 | 2188 | 335 |
| pBlood_2011_0238 | 147680 | 80526 | 54.53% | 56 | 2965 | 53 |
| pUPEC_108 | 173409 | 107884 | 62.21% | 61 | 4008 | 65 |
| pUPEC_134 | 105418 | 68645 | 65.12% | 29 | 3640 | 127 |
| pUPEC_271 | 74232 | 53554 | 72.14% | 34 | 1202 | 35 |
| pUPEC_90 | 75193 | 30052 | 39.97% | 16 | 4488 | 278 |
Таблица 10: Оценки загрязнения хромосомной ДНК. Загрязнение хромосом в плазмидах, очищенных от реципиентных штамм, количественно оценивалось как доля ридов, сопоставленных с хромосомной последовательностью LMB100 по сравнению с плазмидными считываниями. Для каждой плазмиды перечислены считывания плазмиды, общее считывание, хромосомные показания, среднее покрытие хромосом, среднее покрытие плазмид, а также соотношение P:C. Обратите внимание, что pBlood_08_0081_1 отсутствует, потому что его не удалось восстановить у получателя, и эти оценки проводились на реципиентах из-за их однородного хромосомного фона.
Загрязнение хромосом количественно оценивалось как доля чтений, сопоставленных с собранной хромосомной последовательностью LMB100, по сравнению с плазмидными считываниями. Уровень загрязнения хромосом варьировался от 72,1% до 28,0%, со средним показателем 46,9% (Таблица 10). Этот метод, вероятно, переоценивает загрязнение хромосом, поскольку отображения ридов как на плазмиду, так и на хромосому считаются только хромосомными, но всё равно указывает на значительную хромосомную переносность. Несмотря на этот хромосомный перенос, процесс сборки успешно восстановил полные плазмидные последовательности, которые всё ещё соответствовали размеру, предсказанному по мобильности PFGE для всех образцов (Таблица 8, третий и четвёртый столбцы).
Репрезентативные пары конъюгативных плазмид донор-реципиент
Штамм upec271 был выделен из образца мочи в период с 2011 по 2013 год и отнесён к филогруппе B2 и типу последовательности ST131. Анализ PFGE выявил наличие полосы, соответствующей плазмиде примерно 145 кб [Рисунок 3A, (Донор)]. Анализ последовательностей с коротким чтением выявил наличие репликона IncF.

Рисунок 3: Характеристика pUPEC271. (A) PFGE донора и трансконъюганта. Показаны полосы PFGE для донора (D) и трансконъюгантного (T); зелёные стрелки обозначают полосу плазмид, соответствующую pUPEC271. Протокол, используемый для ПФГЕ, описан висследовании 36. (B) Карта pUPEC271. Трансконъюгант был выстрелен против донора, и консенсус (который на 100% идентичен) показан как зелёное кольцо. Это представление было сгенерировано с помощьюProksee 48. Внешнее кольцо показывает местоположение CDS, два внутренних кольца — содержимое GC и смещение, а внутреннее — расположение генов тРНК. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.
Перенос этой плазмиды, которую авторы назвали pUPEC271, была дополнительно подтверждена анализом PFGE с использованием протокола, ранее описанного встатье 36 [Рисунок 3A, T(трансконъюгант)]. Плазмидная ДНК, очищенная из донорского штамма upec271 (pUPEC271_donor), была секвенирована с помощью нанопоры. В соответствии с размером полоски на геле PFGE, последовательность pUPEC271_donor выявила плазмиду размером 144 519 bp (рисунок 3B). Эта плазмида имела репликон типа IncF и содержала гены устойчивости к шести различным классам антибиотиков (см. дополнительную таблицу 1). Авторы также очистили ДНК одного из трансконъюгантов (pUPEC271_TC) и получили последовательность, идентичную донорской. Итоговая сборка, сравниваемая с pUPEC271_donor, показана на рисунке 3B. Это подтвердило, что конъюгативная плазмида, присутствующая в доноре, была перенесена в трансконъюгантный штамм. Технические данные приведены в таблице 9.
Штамм Blood_10_0913 — пример штамма с несколькими плазмидами. Этот образец был собран из крови в период с 2008 по 2013 год, и короткое чтение показало наличие двух репликонов — IncF и IncY репликон. Анализ PFGE этого штамма выявил две плазмиды примерно 70 и 110 кб соответственно (рисунок 4A). Это было подтверждено секвенированием нанопор, которое дало две сборки размером 71 182 и 111 226 пп соответственно. Авторы назвали эти две плазмиды pBlood_10_0913_1 и pBlood_10_0913_2 (рисунок 4B). Однако изначально сборки не сталкивались с циркулярностью из-за обширной гомологии последовательностей в срезе плазмид, несущих механизм сопряжения (раздел обширной гомологии выделен на рисунке 4B). Схема неудачной сборки, созданная с помощьюBandage 47 , показана на рисунке 4C. Полная сборка двух плазмид стала возможной только после агрессивного отбора последовательностей с более коротким считыванием (< 3000 bp).

Рисунок 4: Идентификация и сборка плазмид Blood_10_0913 штамма. (A) PFGE. Зелёные стрелки обозначают плазмидные полосы. Протокол, используемый для PFGE, описан в сборкепоследовательностей 36 (B). Сборка двух меньших плазмид, присутствующих в pBlood_10_0913 образце — pBlood_10_0913_2 и pBlood_10_0913_2. (C) Неудачная сборка. Представление неудачи сборки из-за наличия общей последовательности между двумя различными сборками. Это представление было сгенерировано с помощьюBandage 47. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.
Дополнительная таблица 1: Гены устойчивости, переносимые pUPEC271. Показаны размер плазмиды, репликоны и гены устойчивости, связанные с pUPEC271. Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительные команды. Пожалуйста, нажмите здесь, чтобы скачать этот файл.