$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Mycobacterium bovis Штамм БЦЖ (bacilli de Calmette et Guérin) 1173P2, претерпевающий экспоненциальный рост, подвергался временным рядам (0, 4, 10 и 20 дней) питательного голодания с последующей 6-дневной реанимацией в богатой питательными веществами среде, как это было представлено ранее в Hu et al.7. Малые РНК были выделены из бактериальной культуры с тремя биологическими репликами в каждой из пяти обозначенных временных точек. Библиотеки Illumina были сконструированы с использованием описанного выше рабочего процесса подготовки библиотеки AQRNA-seq (рис. 1) с последующим секвенированием на секвенсоре в Центре BioMicro Массачусетского технологического института. Затем данные секвенирования были обработаны с помощью конвейера анализа данных AQRNA-seq (рис. 2), настроенного для количественной оценки распространенности тРНК.
После ПЦР-амплификации библиотеки кДНК с помощью праймеров секвенирования во всех образцах наблюдалось присутствие продуктов ПЦР размером 175 пар оснований (.н.), что свидетельствует об образовании димеров праймеров. Чтобы уменьшить перенос димеров праймеров, продукты ПЦР размером более 195.о. вырезали из геля и очищали (рис. 3B).
Отфильтрованные по качеству и обрезанные последовательности были сопоставлены с пользовательской библиотекой референсных последовательностей, включающей 45 изоакцепторов тРНК, внутренний стандарт и контрольные последовательности (т.е. 23S рРНК, 16S рРНК, 5S рРНК, rnpB и ssr). Изоакцепторы тРНК составляли от 10,5% до 40,2% от общего числа картированных прочтений данного образца и показали гораздо более высокую распространенность, чем контрольные последовательности (рис. 4). Важно отметить, что относительно низкие доли изоакцепторов тРНК при считывании можно объяснить более высокой относительной распространенностью внутренних стандартов. Таким образом, считанные пропорции изоакцепторов тРНК относительно внутренних стандартов (рис. 4, розовые и зеленые цветовые блоки) могут контролироваться оператором путем тонкой настройки количества внутреннего стандарта, вносимого в реакцию.
Данные о распространенности необработанных тРНК нормализовали с использованием метода медианы отношений, реализованного в версии пакета DESeq2 (далее обозначается как v)1.36.0 16 в R Statistical Programming Environment (далее – R) v 4.2.117. После нормализации достигается количественный ландшафт изоакцепторов тРНК в БЦЖ Mycobacterium bovis во время нутриентного голодания и реанимации (рис. 5).
Чтобы выявить различные кластеры образцов с различными фенотипами на основе закономерностей в распространенности изоакцепторов тРНК, был проведен анализ главных компонент (PCA) нормализованных данных о распространенности тРНК с использованием пакета stats v 4.2.117 in R (рис. 6). В ходе анализа были выделены образцы 0-го дня голодания и 6-го дня реанимации из образцов 4-го, 10-го и 20-го дней голодания, что свидетельствует о значительной разнице в ландшафте тРНК Mycobacterium bovis BCG, выращенных в среде, лишенной питательных веществ, и в среде, богатой питательными веществами.
Для профилирования динамики распространенности каждого изоакцептора тРНК в пяти обозначенных временных точках был проведен дифференциальный анализ экспрессии нормализованных данных о распространенности тРНК с использованием пакета DESeq2 v 1.36.0 в R (рис. 7). Анализ показал, что 17 из 20 семейств изоакцепторов содержали изоакцепторы, которые дифференциально экспрессировались (т.е. значительно повышались или понижались) по крайней мере в одном из временных моментов, что предполагает потенциальную роль регуляции пула тРНК в персистирующем состоянии Mycobacterium bovis BCG при туберкулезе.

Рисунок 1: Схема рабочего процесса подготовки библиотеки AQRNA-seq. Ключевые шаги, описанные в рабочем процессе, перечислены в центре схемы и соединены с соответствующими графическими иллюстрациями пунктирными линиями. Подробное описание каждого шага можно найти в разделе Протокол. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этой цифры.

Рисунок 2: Схема конвейера анализа данных AQRNA-seq. Основные этапы, описанные в конвейере, перечислены в центре схемы и соединены с соответствующими графическими иллюстрациями пунктирными линиями. Подробное описание каждого шага доступно на GitHub (https://github.com/Chenrx9293/AQRNA-seq-JoVE.git). Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этой цифры.

Рисунок 3: Электрофорез фрагментов кДНК в агарозном геле после амплификации ПЦР с секвенированием праймеров. (A) Изображение геля до экстракции и очистки геля. Дорожки 7 и 14 с левой стороны содержат 5 μL из 50.о. лестницы ДНК, в то время как другие дорожки содержат по 20 μL каждого из 15 образцов. Размерная локализация продуктов ПЦР указывает на их максимальную концентрацию в диапазоне от 175.о. (димеры праймеров) до 300.о. (два праймера + 120.о. 5S рРНК). (Б) Изображение геля после экстракции и очистки геля. Для каждого образца гелевый блок в диапазоне от 200 до 400.о. был исключён, чтобы свести к минимуму загрязнение димерами праймеров в библиотеке секвенирования. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этой цифры.

Рисунок 4: Количество успешно сопоставленных операций чтения последовательностей с библиотекой эталонных последовательностей. По оси X отображаются названия образцов (например, D18-69XX), сгруппированных по временным точкам (например, День голодания 0). Для каждой выборки количество прочтений, связанных с различными категориями целевых предметов, представлено с помощью цветных блоков, расположенных друг над другом. Числа, расположенные в центре цветовых блоков, представляют собой пропорции прочтений, соответствующих соответствующим целевым объектам в данной выборке. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этой цифры.

Рисунок 5: Количественный ландшафт изоакцепторов тРНК Mycobacterium bovis BCG в различные моменты времени голодания и реанимации. Данные о распространенности необработанных тРНК нормализовались с помощью метода медианы отношений. Здесь каждая строка изображает нормализованное содержание тРНК (ось y) в виде среднего ± стандартной ошибки для 3 биологических репликатов в каждый момент времени. По оси x изоакцепторы из одного семейства были сгруппированы вместе и помечены соответствующей аминокислотой. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этой цифры.

Рисунок 6: Квадратичный косинусный график образцов, полученных в результате анализа главных компонент (PCA). ПЦА проводили на основе нормализованного содержания тРНК. Квадрат косинуса указывает на важность главных компонент для образцов, и образцы были построены относительно квадрата косинуса первых двух главных компонентов. Образцы были помечены с использованием идентификаторов образцов и цветовой кодировкой по временным точкам. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этой цифры.

Рисунок 7: Дифференциальная экспрессия изоакцепторов тРНК в разные моменты времени. Нормализованное содержание тРНК было обобщено в виде средних значений (линейных узлов) ± стандартной ошибки (полос погрешностей) по 3 биологическим репликам. Из-за ограниченности места условия были сокращены следующим образом: S0-S20 = голодные дни 0-20; R6 = 6-й день реанимации. Дифференциальный анализ экспрессии был проведен для каждого изоакцептора тРНК, сравнивая различные временные точки попарно с использованием теста отношения правдоподобия и теста Вальда. Компактные буквы использовались для представления статистической значимости, когда распространенность данного изоакцептора тРНК в моменты времени, имеющие по крайней мере одну общую букву, существенно не отличались друг от друга. Например, содержание тРНК-Lys-CTT-1-1 (в лизиновой панели) было значительно снижено от S0 до S4 и от S4 до S10, но не от S10 до S20. Затем он был значительно отрегулирован с S20 до R6. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этой цифры.
Таблица 1: Олигонуклеотиды, участвующие в процессе подготовки библиотеки AQRNA-seq. Внутренним стандартом является РНК, в то время как все остальные олигонуклеотиды являются ДНК. Перечисленные праймеры для ПЦР и пользовательские праймеры для секвенирования специфичны для платформ секвенирования. Дополнительные праймеры для ПЦР могут быть разработаны с помощью новых индексных последовательностей. Пожалуйста, нажмите здесь, чтобы скачать эту таблицу.