Этот протокол был разработан и протестирован на основе образцов, собранных в рамках исследований, одобренных Институциональным обзорным советом UCSF.

Рисунок 1: Обработка образцов и аналитический рабочий процесс для обнаружения патогенов в области метагеномного секвенирования мозговой жидкости (ЦСЖ) с помощью спинномозговой жидкости (ЦСЖ).Эта диаграмма демонстрирует основные этапы влажного лабораторного и биоинформатического путей мСП (CSF). Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.
1. Вопросы обращения с образцами CSF
ПРИМЕЧАНИЕ: результаты mNGS оптимальны, если подготовка библиотеки ликвора происходит во время сбора или если образцы СМЖ сохраняются методом мгновенного замораживания в жидком азоте с последующим длительное хранение при -80 °C 15. Альтернативно, использование консервантных растворов, таких как ДНК/РНК-щит, во время сбора образцов позволяет хранить их при комнатной температуре более длительноевремя — 16,17. Это может быть предпочтительнее в ситуациях, когда возникают сомнения по поводу надёжности холодовой цепи. Кроме того, предотвращение ненужных циклов замораживания и оттаивания также помогает сохранить целостность нуклеиновых кислот, и может быть разумно предварительно аликвотировать образцы при сборе, если для образцов СМЖ ожидается несколько аналитических процессов18,19.
2. Экстракция нуклеиновой кислоты ликвора
ВНИМАНИЕ: Выполняйте работу в соответствующем шкафу с уровнем биобезопасности в перчатках и лабораторном халате, чтобы снизить риск передачи инфекционных организмов. Реагенты вредны при глотке, вдыхании или контакте с кожей. Утилизировать использованные материалы как опасные химические отходы в соответствии с местными нормативами.
- Набор быстрой ДНК/РНК обычно хорошо работает для образцов с низким содержанием нуклеиновых кислот.
- Разморозить образцы СМФ на льду, затем перенести по 100 мкл–1 мл каждого образца в трубку на 1,5 мл и центрифугировать при 16 000 г при 4 °C в течение 10 минут.
- Аккуратно уберите пипетку от супернатанта, оставив после себя гранулу (которая часто невидима невооружённым глазом) и 100 мкл супернатанта.
- Добавьте 100 мкл ДНК/РНК-щита к оставшимся 100 мкл супернатанта и гранул, и тщательно перемешайте пипетку.
- Затем добавьте 600 мкл буфера лизиса ДНК/РНК и смесьте пипетку, пока раствор не станет однородным.
- Следуйте оставшейся части протокола экстракции ДНК/РНК, описанному производителем в дополнительных приложениях A, со следующими важными изменениями, специфичными для ликвордовой жидкости:
- Выполните все этапы центрифугации в течение 1 минуты на максимальной скорости (~21 000 x g), кроме финального этапа промывания ДНК/РНК, который должен длиться 3 минуты на максимальной скорости.
- Элюируйте образец, добавив 22 мкл (вместо 25 мкл) безнуклеазной воды непосредственно на матрицу столбца, затем инкубируйте 3 минуты при комнатной температуре. Больший объём обычно приводит к большему общему объёму нуклеиновой кислоты, элюируемой из столбца.
- После элюции пропустите проток обратно на матрицу колонки и снова центрифугуйте в течение одной минуты на максимальной скорости. Перенесите образец в новую накопительную трубку объёмом 1,5 мл и либо сразу приступите к подготовке библиотеки, либо храните при -80 °C.
3. Подготовка библиотеки РНК CSF
ВНИМАНИЕ: Выполняйте работу в соответствующем шкафу для биобезопасности, в перчатках и лабораторном халате. Реагенты вредны при глотке, вдыхании или контакте с кожей. Утилизировать использованные материалы как опасные химические отходы в соответствии с местными нормативами.
ПРИМЕЧАНИЕ: Важные соображения перед началом: истощение рРНК предпочтительнее полиаденилированного обогащения для выделения мРНК, поскольку истощение рРНК не требует дополнительных этапов очистки шариков (что может снизить доход полезной мРНК), а некоторые транскрипты нейроинвазивных вирусов РНК не имеют полиаденилирования. ERCC (консорциум контроля внешних РНК) РНК Spike-In может быть добавлен в качестве меры контроля качества, хотя это не является строго обязательным. Если РНК хранилась при -80 °C после извлечения, оттаивайте на льду до начала.
- Фрагментация и прайминг РНК: Пипетт смешивает следующее в пробирке полимеразной цепной реакции (ПЦР): образец РНК 3,5 мкл, 0,5 мкл 1:2500 разбавленный ERCC spike-in, 4 мкл буфера реакции первой цепи (5 раз), 1 мкл случайных праймеров и 1 мкл 1:100 разбавленной рРНК (таблица материалов). В термоциклере установите нагретую крышку на 105 °C и инкубуйте образцы при 75 °C в течение 2 минут, 70 °C 2 минуты, 65 °C 2 минуты, 60 °C 2 минуты, 55 °C 2 минуты, 37 °C 5 минут и 25 °C 5 минут.
- Синтез кДНК первой цепи: Pipette смешивает фрагментированную и праймированную РНК (10 мкл) с 8 мкл безнуклеазной водой и 2 мкл ферментов синтеза первой цепи. В термоциклере установите нагретую крышку на 105 °C и инкубируете образцы при 25 °C в течение 10 минут, 42 °C 15 минут и 70 °C 15 минут.
- Синтез кДНК второй цепи: Пипетт смешивает синтезированную ДНК первой цепи (20 мкл) с буфером реакции синтеза второй цепи 8 мкл, смесью ферментов синтеза второй цепочки 4 мкл и 48 мкл водой без нуклеазы. В термоциклере инкубировать образцы в течение 1 часа при 16 °C, с нагретой крышкой.
- Выполните очистку магнитных шариков (Дополнительные приложения B) с использованием соотношения 1,8x магнитных шариков SPRI (144 мкл). Элюировать в 53 мкл безнуклеазную воду и перенести 50 мкл последнего супернатанта (содержащего очищенную двухцепочечную кДНК) в чистую ПЦР-трубку без нуклеазы. В этот момент образцы можно безопасно заморозить при -20 °C на ночь, если это необходимо.
- Завершение подготовки библиотеки кДНК: если образец хранили при -20 °C ночью, разморозить на льду перед повторным началом. Пипетты смешивают очищенную ds-cDNA (50 мкл) с 7 мкл буфера для последней подготовки реакции и 3 мкл ферментной смеси для подготовки концов.
- В термоциклере инкубировать образцы при 20 °C в течение 30 минут и при 65 °C в течение 30 минут, при закрытой нагретой крышке.
- Перевязка адаптера (выполните этот этап на льду): Создайте разбавление адаптера в соотношении 1:100 в воде без нуклеаз. Пипетировать конечную реакционную смесь (60 мкл) в трубку с мастер-смесью 30 мкл лигации, 1 мкл усилителем лигации и 2,5 мкл разбавлённого адаптера 1:100. Убедитесь, что адаптер добавлен отдельно (например, не заранее смешан с мастер-миксом или усилителем лигации), чтобы избежать образования димеров адаптера.
- В термоциклере инкубировать образцы при 20 °C в течение 15 минут, с нагретой крышкой.
- Немедленно переходите к очистке магнитных шариков (Дополнительное приложение 2) с использованием соотношения 0,9x магнитных шариков SPRI (87 мкл).
- Элюировать в воду без нуклеазы 17 мкл и перевести 15 мкл последнего супернатанта в чистую ПЦР-трубку без нуклеазы
- Штрих-кодирование ПЦР: Пипетте очищенную кДНК с адаптером (15 мкл) в пробирку с 3 мкл ферментом USER, 25 мкл Q5 Master Mix и 10 мкл уникальных штрих-кодированных праймеров. В термоциклере установите нагретую крышку на 105 °C, инкубируете образцы при 37 °C в течение 15 минут, затем при 98 °C 30 секунд, а затем проведите 19 циклов при температуре 98 °C в течение 10 с и при 65 °C в течение 75 секунд. Завершите ПЦР, инкубировав при 65 °C в течение 5 минут.
- Выполните окончательную очистку магнитных шариков с помощью соотношения магнитных шариков 0,8x (43 мкл). Элюировать в воду без нуклеазы 23 мкл и после последнего этапа разделения шариков передать 20 мкл в чистую ПЦР-трубку без нуклеаз. В этот момент образцы можно безопасно заморозить при -20 °C на ночь, если это необходимо.
4. Подготовка библиотеки ДНК CSF
ВНИМАНИЕ: Проводите работу в шкафчике соответствующего уровня биобезопасности в перчатках и лабораторном халате. Реагенты вредны при глотке, вдыхании или контакте с кожей. Утилизировать использованные материалы как опасные химические отходы в соответствии с местными нормативами.
- Если извлеченная ДНК хранилась при -20 °C ночью, её разморозили на льду перед повторным началом.
- Фрагментация ДНК и подготовка на конце: Тщательно смешайте буфер реакции первой цепи путём вихря и пипетки, смешивая раствор для повторной суспензии осадка. Пипетная смесь содержит 3,5 мкл экстракции ДНК, 22,5 мкл воды без нуклеаз, 7 мкл буфера реакции первой цепи и 2 мкл смеси ферментов первой цепи.
- В термоциклере установите нагретую крышку на 105 °C и инкубируете образцы при 37 °C в течение 5 минут, затем при 65 °C в течение 30 минут.
- Процесс ПЦР с помощью адаптера и штрихкодирования идентичен для библиотек ДНК и РНК. Повторите разделы 7–10 «Подготовка библиотеки РНК по ЦВО» с фрагментированной и готовой на концах ДНК, чтобы завершить подготовку библиотеки ДНК.
5. Контроль качества библиотеки и пулирование
- Количественно определите концентрацию библиотек образцов с помощью комплекта для количественной оценки ДНК и соответствующего флюорометра. Водные контроли должны иметь значительно меньшие (или неизмеримые) концентрации по сравнению с остальными образцами.
- Оценивайте размеры библиотеки образцов, запуская образцы на автоматическом капиллярном электрофорезном аппарате.
- В качестве альтернативы, если автоматизированная капиллярная электрофорезная машина недоступна, дополнительно усилите 1-2 мкл образца, проведя ПЦР-реакцию с универсальными праймерами Illumina, а затем проведите электрофорез 10 мкл готового продукта на 2% агарозном геле. Утилизировать использованные материалы как опасные химические отходы в соответствии с местными нормативами.
ПРИМЕЧАНИЕ: Оптимальная длина библиотеки составляет примерно 400–500 пар оснований, поэтому будущие инкубации фрагментации следует корректировать, если размеры библиотек будут больше или короче желаемого. Кроме того, крайне важно оценить наличие адаптерных димеров, которые имеют примерно 150 пар оснований в длину и возникают, когда две адаптерные молекулы случайно связываются друг с другом без вставной последовательности. Даже если адаптерные димеры присутствуют в небольших количествах, они обычно эффективно кластерируются и секвенируются на ячейке потока, уменьшая долю пригодных считываний во время секвенирования. Таким образом, при наличии адаптерных димеров проводите дополнительную очистку магнитных шариков с соотношением 0,8x магнитных шариков SPRI. (Стоит отметить, что если >10% образца — адаптерный димер, может потребоваться 2-3 круга очистки магнитных шариков, чтобы удалить большинство адаптерных димеров).
- Для обеспечения одинаковой глубины секвенирования между образцами в одном и том же цикле необходимо объединить эквимолярные величины каждого образца. Если размеры библиотек сопоставимы для всех образцов, просто объедините эквивалентные массы каждой библиотеки РНК и повторите это отдельно для каждой библиотеки.
6. Вопросы по секвенированию
- Объединённые образцы теперь готовы к секвенированию на устройствах секвенирования Illumina. Если секвенирование проводится внутри компании (а не через специализированное ядро секвенирования), аккуратно денатурируйте и разбавляйте объединённые библиотеки в соответствии с инструкциями устройства секвенирования перед загрузкой пула на картридж реагента (типичная концентрация составляет 4 нМ для большинства платформ секвенирования Illumina).
- Рассчитайте предполагаемую глубину секвенирования, деля ожидаемый выход данных для ячейки потока Illumina на количество мультиплексированных образцов в пробеге.
ПРИМЕЧАНИЕ: Определение желаемой глубины секвенирования зависит от нескольких факторов, включая стоимость, используемое устройство секвенирования и цели эксперимента. Два примера различных глубин секвенирования приведены ниже в разделе «Репрезентативные результаты».
7. Анализ данных для обнаружения патогенов
- Используйте открытую веб-платформу Chan Zuckerberg ID (CZID, https://www.czid.org, Illumina mNGS Pipeline v8.3) для проведения метагеномного анализа с целью обнаружения патогенов. CZID также включает опции загрузки сэмплов через командный интерфейс и прямой передачи из аккаунта Illumina BaseSpace в CZID. Подробный обзор анализа данных CZID mNGS доступен по адресу https://chanzuckerberg.zendesk.com/hc/en-us/articles/13770737266196-Guide-to-mNGS-Data-Analysis
- Войдите в CZID, нажмите «загрузить», выберите «метагеномика» для типа анализа и выберите входные файлы fastq из секвенирования. Вводите необходимую информацию для каждого образца, включая название образца, тип образца (CSF) и нуклеотид (РНК или ДНК).
ПРИМЕЧАНИЕ: Загруженные файлы проходят автоматическую обработку через конвейер CZID, который включает несколько вычислительных инструментов, таких как STAR, Bowtie2, Trimmomatic, PriceSeq и GSNAP, среди прочих, для удаления низкокачественных чтений и последовательностей чтения человека. Конвейер затем выравнивается и собирает фильтрованные последовательности с помощью баз данных нуклеотидов (NT) и белков (NR) Национального центра биотехнологической информации (NCBI) через программы Minimap2 и DIAMOND соответственно. Итоговый результат CZID включает подсчёты таксонов NT/NR и контиг-счёты (смежные, перекрывающиеся сегменты, создаваемые при сборке) для каждого образца.
- После обработки нажмите на папку проекта с образцами и перейдите к панели сводки.
- Прокрутите вниз, чтобы изучить количество считываний на образец (и сравните с ожидаемым количеством считываний на основе используемой ячейки Illumina и количества образцов, мультиплексированных в ходе секвенирования), процент считываний, проходящих фильтрацию контроля качества, и коэффициент дублирования сжатия для оценки предвзятого переусиления ПЦР.
- Наличие загрязнения фрагментами нуклеиновых кислот окружающей среды — либо в процессе сбора образцов ликвора, либо при извлечении/подготовке в библиотеке — является распространённымявлением 20. Чтобы определить, являются ли бактериальные, грибковые и паразитические показания истинно положительными, можно рассмотреть потенциальные методы, использующие нормированные показатели чтений (чтения на миллион, или rPM):
- Создайте фоновую модель (например, используя образцы воды, извлечённые, подготовленные и секвенированные) непосредственно внутри CZID, чтобы отфильтровать таксоны с rPM, которые схожи или меньше по количеству по сравнению с фоновой моделью (и, следовательно, могут представлять загрязнение). В папке проекта проверьте каждый из образцов контроля воды и нажмите «фоновая модель» для создания. При анализе каждого образца выберите фоновую модель и нажмите «threshold filters», чтобы применить «NTZ Score» > 1 (или выше) для фильтрации вероятных загрязнителей.
- Используйте заранее заданные пороги rPM, такие каквыборка rpM > 10* (средний rPM контрольный) или log10 (выборка rPM) как минимум на 1 лог больше log 10 (средний rPMпо всей когорте). Сначала скачайте специфичный для каждого образца rPM для каждого образца, проверив каждый образец на странице проекта, нажав «Скачать и выбрав «Результаты комбинированных выборочных таксонов». Полученный файл .csv можно анализировать с помощью статистического программного обеспечения, такого как R или Stata, по желанию.
ПРИМЕЧАНИЕ: Учитывая низкую относительную численность вирусных патогенов, любой вирус с известным нейроинвазивным потенциалом и хотя бы одно считывание, совпадающее с вирусным геномом, должны считаться положительным результатом (и желательно подтверждены повторным mNGS или специфическим клиническим тестированием патогена). Кроме того, аналитический конвейер требует , чтобы каждое чтение файла fastq было вычислительно согласовано с наиболее вероятным эталонным геномом базы данных NCBI. В результате некоторые варианты могут быть относительно неспецифичными для конкретного патогена и теоретически соответствовать широкому кругу организмов, а не конкретному таксону, указанному в CZID. Учитывая этот риск, ручное подтверждение положительных вызовов следует выполнять путем подтверждения выравнивания через веб-инструмент NCBI BLAST (конкретные считывания таксона можно напрямую отправить из CZID в NCBI BLAST, нажав на значок взрыва рядом с соответствующим таксоном).