$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Mycobacterium bovis La cepa BCG (bacilos de Calmette et Guérin) 1173P2 en crecimiento exponencial se sometió a una serie temporal (0, 4, 10 y 20 días) de inanición de nutrientes, seguida de una reanimación de 6 días en un medio rico en nutrientes, como se presentó previamente en Hu et al.7. Se aislaron pequeños ARN de cultivo bacteriano, con tres réplicas biológicas, en cada uno de los cinco puntos de tiempo designados. Las bibliotecas de Illumina se construyeron utilizando el flujo de trabajo de preparación de bibliotecas AQRNA-seq descrito anteriormente (Figura 1), seguido de la secuenciación en un secuenciador en el Centro BioMicro del Instituto Tecnológico de Massachusetts. A continuación, los datos de secuenciación se procesaron utilizando la canalización de análisis de datos AQRNA-seq (Figura 2) personalizada para la cuantificación de la abundancia de ARNt.
Después de la amplificación por PCR de la biblioteca de ADNc con los cebadores de secuenciación, se observó la presencia de productos de PCR con un tamaño de 175 pares de bases (pb) en todas las muestras (Figura 3A), lo que sugiere la formación de dímeros de cebadores. Para mitigar el arrastre de dímeros de cebador, los productos de PCR con un tamaño superior a 195 pb se extirparon del gel y se purificaron (Figura 3B).
Las lecturas de secuencias filtradas y recortadas por calidad se asignaron a una biblioteca de secuencias de referencia personalizada que incluía los 45 isoaceptores de ARNt, el estándar interno y las secuencias de control (es decir, ARNr 23S, ARNr 16S, ARNr 5S, rnpB y ssr). Los isoaceptores de ARNt representaron entre el 10,5% y el 40,2% del total de lecturas mapeadas de una muestra dada y mostraron una abundancia mucho mayor que las secuencias de control (Figura 4). Es importante destacar que las proporciones de lectura relativamente bajas de los isoaceptores de ARNt podrían atribuirse a la mayor abundancia relativa de los patrones internos. Por lo tanto, el operador puede controlar las proporciones de lectura de los isoaceptores de ARNt en relación con los estándares internos (Figura 4, bloques de color rosa frente a verde) mediante el ajuste fino de la cantidad de patrón interno introducido en la reacción.
Los datos brutos de abundancia de ARNt se normalizaron utilizando el método de mediana de proporciones implementado con la versión del paquete DESeq2 (en adelante, v) 1.36.016 en R Statistical Programming Environment (en adelante, R) v 4.2.117. Después de la normalización, se logra un panorama cuantitativo de los isoaceptores de ARNt en Mycobacterium bovis BCG durante un curso de tiempo de inanición de nutrientes y reanimación (Figura 5).
Para revelar distintos grupos de muestras con diferentes fenotipos basados en patrones en la abundancia de isoaceptores de ARNt, se realizó un Análisis de Componentes Principal (PCA) sobre los datos de abundancia de ARNt normalizados utilizando el paquete de estadísticas v 4.2.117 en R (Figura 6). El análisis distinguió las muestras del día 0 de inanición y el día 6 de reanimación de las muestras de los días 4, 10 y 20 de inanición, lo que sugiere una diferencia considerable en el panorama del ARNt de Mycobacterium bovis BCG cultivado en medio privado de nutrientes y medio rico en nutrientes.
Para perfilar la dinámica de la abundancia de cada isoaceptor de ARNt en los cinco puntos de tiempo designados, se realizó un análisis de expresión diferencial en los datos de abundancia de ARNt normalizados utilizando el paquete DESeq2 v 1.36.0 en R (Figura 7). El análisis reveló que 17 de las 20 familias de isoaceptores contenían isoaceptores que se expresaban diferencialmente (es decir, significativamente regulados al alza o a la baja) en al menos uno de los puntos temporales, lo que sugiere un papel potencial de la regulación del grupo de ARNt en el estado persistente de Mycobacterium bovis BCG durante la tuberculosis.

Figura 1: Esquema del flujo de trabajo de preparación de la biblioteca AQRNA-seq. Los pasos clave descritos en el flujo de trabajo se enumeran en el centro del esquema y se conectan a sus respectivas ilustraciones gráficas mediante líneas de puntos. La descripción detallada de cada paso se puede encontrar en la sección Protocolo. Haga clic aquí para ver una versión más grande de esta figura.

Figura 2: Esquema de la canalización de análisis de datos AQRNA-seq. Los pasos clave descritos en la canalización se enumeran en el centro del esquema y están conectados a sus respectivas ilustraciones gráficas mediante líneas punteadas. La descripción detallada de cada paso está disponible en GitHub (https://github.com/Chenrx9293/AQRNA-seq-JoVE.git). Haga clic aquí para ver una versión más grande de esta figura.

Figura 3: Electroforesis en gel de agarosa de los fragmentos de ADNc después de la amplificación por PCR con cebadores de secuenciación. (A) Imagen del gel antes de la extracción y purificación del gel. Los carriles 7 y 14 del lado izquierdo contienen 5 μL de la escalera de ADN de 50 pb, mientras que los otros carriles contienen 20 μL de cada una de las 15 muestras. La localización del tamaño de los productos de PCR indica su concentración más alta dentro del rango de 175 pb (dímeros de cebador) a 300 pb (dos cebadores + 120 pb ARNr 5S). (B) Imagen del gel después de la extracción y purificación del gel. Para cada muestra, se extirpó el bloque de gel entre 200 pb y 400 pb para minimizar la contaminación de los dímeros de cebadores en la biblioteca de secuenciación. Haga clic aquí para ver una versión más grande de esta figura.

Figura 4: Número de lecturas de secuencia asignadas con éxito a la biblioteca de secuencias de referencia. El eje x muestra los nombres de las muestras (por ejemplo, D18-69XX) agrupadas por punto de tiempo (por ejemplo, Día 0 de inanición). Para cada muestra, el recuento de lecturas asociado con varias categorías de temas objetivo se representa mediante bloques de color apilados uno encima del otro. Los números ubicados en el centro de los bloques de color representan las proporciones de lecturas correspondientes a los respectivos sujetos objetivo dentro de una muestra determinada. Haga clic aquí para ver una versión más grande de esta figura.

Figura 5: Panorama cuantitativo de los isoaceptores de ARNt de Mycobacterium bovis BCG en varios puntos temporales a lo largo del curso del tiempo de inanición y reanimación. Los datos brutos de abundancia de ARNt se normalizaron utilizando el método de la mediana de las proporciones. Aquí, cada fila representa las abundancias normalizadas de ARNt (eje Y) como media ± error estándar para 3 réplicas biológicas en cada punto de tiempo. En el eje x, los isoaceptores de la misma familia se agruparon y se marcaron con el aminoácido correspondiente. Haga clic aquí para ver una versión más grande de esta figura.

Figura 6: Diagrama de coseno cuadrado de muestras derivadas del análisis de componentes principales (PCA). La PCA se realizó en base a la abundancia normalizada de ARNt. El coseno al cuadrado indica la importancia de los componentes principales para las muestras, y las muestras se trazaron con respecto al coseno al cuadrado de los dos primeros componentes principales. Las muestras se etiquetaron con identificadores de muestra y se codificaron por colores según el punto de tiempo. Haga clic aquí para ver una versión más grande de esta figura.

Figura 7: Expresión diferencial de isoaceptores de ARNt en diferentes puntos temporales. Las abundancias normalizadas de ARNt se resumieron como medias (nudos de línea) ± error estándar (barras de error) en 3 réplicas biológicas. Debido a la limitación de espacio, las condiciones se abreviaron de la siguiente manera: S0-S20 = días de inanición 0-20; R6 = día de reanimación 6. Se realizó un análisis de expresión diferencial para cada isoaceptor de ARNt, comparando varios puntos temporales de forma pareada utilizando la prueba de razón de verosimilitud y la prueba de Wald. Se emplearon letras compactas para representar la significación estadística, donde las abundancias de un isoaceptor de ARNt dado en puntos de tiempo que compartían al menos una letra común no eran significativamente diferentes entre sí. Por ejemplo, la abundancia de tRNA-Lys-CTT-1-1 (en el panel de lisina) se reguló significativamente a la baja de S0 a S4 y de S4 a S10, pero no de S10 a S20. Luego se reguló significativamente de S20 a R6. Haga clic aquí para ver una versión más grande de esta figura.
Tabla 1: Oligonucleótidos implicados en el flujo de trabajo de preparación de la biblioteca AQRNA-seq. El estándar interno es el ARN, mientras que todos los demás oligonucleótidos son ADN. Los cebadores de PCR y los cebadores de secuenciación personalizados enumerados son específicos para las plataformas de secuenciación. Se pueden diseñar cebadores de PCR adicionales con secuencias de índice novedosas. Haga clic aquí para descargar esta tabla.