Artículo de método

AQRNA-seq para cuantificar ARNs pequeños

DOI:

10.3791/66335

2 de febrero de 2024

En este artículo

Resumen

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

La secuenciación de ARN de cuantificación absoluta (AQRNA-seq) es una tecnología desarrollada para cuantificar el panorama de todos los ARN pequeños en mezclas biológicas. Aquí, se demuestran los pasos de preparación de la biblioteca y procesamiento de datos de AQRNA-seq, cuantificando los cambios en el grupo de ARN de transferencia (ARNt) en Mycobacterium bovis BCG durante la latencia inducida por inanición.

Resumen

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

AQRNA-seq proporciona una relación lineal directa entre los recuentos de lecturas de secuenciación y el número pequeño de copias de ARN en una muestra biológica, lo que permite una cuantificación precisa del grupo de ARN pequeños. El procedimiento de preparación de la biblioteca AQRNA-seq descrito aquí implica el uso de enlazadores de secuenciación diseñados a medida y un paso para reducir las modificaciones del ARN de metilación que bloquean la procesividad de transcripción inversa, lo que da como resultado un mayor rendimiento de ADNc de longitud completa. Además, se presenta una implementación detallada de la cartera de productos bioinformáticos que la acompaña. Esta demostración de AQRNA-seq se llevó a cabo a través de un análisis cuantitativo de los 45 ARNt en Mycobacterium bovis BCG cosechados en 5 días seleccionados a lo largo de un curso de 20 días de privación de nutrientes y 6 días de reanimación. Los esfuerzos en curso para mejorar la eficiencia y el rigor de AQRNA-seq también se discutirán aquí. Esto incluye explorar métodos para obviar la purificación en gel para mitigar los problemas de dímero de cebador después de la amplificación de PCR y para aumentar la proporción de lecturas completas para permitir un mapeo de lectura más preciso. Las futuras mejoras de AQRNA-seq se centrarán en facilitar la automatización y la implementación de alto rendimiento de esta tecnología para cuantificar todas las especies pequeñas de ARN en muestras de células y tejidos de diversos organismos.

Introducción

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

La secuenciación de próxima generación (NGS), también conocida como secuenciación paralela masiva, es una tecnología de secuenciación de ADN que implica la fragmentación del ADN, la ligadura de oligonucleótidos adaptadores, la amplificación basada en la reacción en cadena de la polimerasa (PCR), la secuenciación del ADN y el reensamblaje de las secuencias de fragmentos en un genoma. La adaptación de NGS a la secuencia de ARN (RNA-seq) es un enfoque poderoso para identificar y cuantificar las transcripciones de ARN y sus variantes1. Los desarrollos innovadores en los flujos de trabajo de preparación de bibliotecas de ARN y las tuberías de análisis bioinformático, junto con los avances en la instrumentación de laboratorio, han ampliado el repertorio de aplicaciones de secuenciación de ARN, progresando más allá de la secuenciación del exoma hacia ómicas funcionales avanzadas como el perfil de ARN no codificante2, el análisis de células individuales3, la transcriptómica espacial 4,5, el análisis de empalme alternativo6entre otros. Estos métodos avanzados de secuenciación de ARN revelan funciones complejas de ARN a través del análisis cuantitativo del transcriptoma en células y tejidos normales y enfermos.

A pesar de estos avances en la secuenciación de ARN, varias características técnicas clave limitan el poder cuantitativo del método. Si bien la mayoría de los métodos de secuenciación de ARN permiten una cuantificación precisa y precisa de los cambios en los niveles de ARN entre variables experimentales (es decir, muestras biológicas y/o estados fisiológicos), no pueden proporcionar comparaciones cuantitativas de los niveles de moléculas de ARN dentro de una muestra. Por ejemplo, la mayoría de los métodos de secuenciación de ARN no pueden cuantificar con precisión el número relativo de copias de moléculas isoaceptoras de ARNt individuales en un grupo celular de ARNt expresados. Como se destaca en la publicación complementaria7, esta limitación a la secuenciación de ARN surge de varias características de la estructura del ARN y la bioquímica de la preparación de la biblioteca. Por ejemplo, la actividad de las enzimas de ligadura utilizadas para unir los enlazadores de secuenciación de los extremos 3' y 5' a las moléculas de ARN está fuertemente influenciada por la identidad de los nucleótidos terminales del ARN y los enlazadores de secuenciación. Esto conduce a grandes variaciones en la eficiencia de las ligaduras de enlazadores y profundos aumentos artifactuales en las lecturas de secuenciación 8,9,10.

Un segundo conjunto de limitaciones surgen de las propiedades estructurales inherentes de las moléculas de ARN. Específicamente, la formación de la estructura secundaria del ARN y los cambios dinámicos en las docenas de modificaciones postranscripcionales del ARN del epitranscriptoma pueden causar la caída de la polimerasa o la mutación durante la transcripción inversa. Estos errores dan lugar a una síntesis de ADNc incompleta o truncada o a una secuencia de ARN alterada. Si bien ambos fenómenos pueden explotarse para mapear estructuras secundarias o algunas modificaciones, degradan la precisión cuantitativa de la secuenciación de ARN si los pasos posteriores de preparación de la biblioteca no logran capturar ADNc truncados o si el procesamiento de datos arroja secuencias mutadas que no coinciden con un conjunto de datos de referencia11,12. Además, la inmensa diversidad química, estructural y de longitud de los transcritos de ARN, así como la falta de herramientas para fragmentar uniformemente los ARN largos, disminuye la aplicabilidad de la mayoría de los métodos de secuenciación de ARN a todas las especies de ARN13.

El método AQRNA-seq (secuenciación de ARN de cuantificación absoluta) se ha desarrollado para eliminar varias de estas restricciones técnicas y biológicas que limitan la precisión cuantitativa7. Al minimizar los sesgos dependientes de la secuencia en la captura, ligadura y amplificación durante la preparación de la biblioteca de secuenciación de ARN, AQRNA-seq logra una linealidad superior en comparación con otros métodos, cuantificando con precisión el 75% de una biblioteca de referencia de 963 miARN con una precisión de 2 veces. Esta correlación lineal entre el recuento de lecturas de secuenciación y la abundancia de ARN también se observa en un análisis de un grupo de longitud variable de estándares de oligonucleótidos de ARN y en referencia a métodos ortogonales como el Northern blot. El establecimiento de la linealidad entre el recuento de lecturas de secuenciación y la abundancia de ARN permite a AQRNA-seq lograr una cuantificación precisa y absoluta de todas las especies de ARN dentro de una muestra.

A continuación, se muestra una descripción del protocolo para el flujo de trabajo de preparación de la biblioteca AQRNA-seq y la canalización de análisis de datos descendente que lo acompaña. El método se aplicó para dilucidar la dinámica de la abundancia de ARNt durante la latencia inducida por inanición y la posterior reanimación en el modelo de tuberculosis Mycobacterium bovis bacilli de Calmette et Guérin (BCG). Se presentaron los resultados para la visualización exploratoria de los datos de secuenciación, junto con los posteriores análisis de agrupamiento y expresión diferencial que revelaron patrones discernibles en la abundancia de ARNt asociados con varios fenotipos.

Protocolo

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

NOTA: La Figura 1 proporciona una ilustración gráfica de los procedimientos involucrados en la preparación de la biblioteca AQRNA-seq. La información detallada sobre los reactivos, productos químicos y columnas/kits utilizados en el procedimiento se puede encontrar en la Tabla de Materiales. Se recomienda realizar una evaluación exhaustiva de la pureza, integridad y cantidad de las muestras de ARN de entrada utilizando (i) electroforesis en gel de agarosa al 3%, (ii) herramientas de electroforesis automatizadas para el control de la calidad de las muestras de biomoléculas (véase la Tabla de Materiales), y (iii) espectrofotometría UV-Visible y/o cuantificación fluorométrica. Es obligatorio mantener todas las reacciones y mezclas maestras en hielo, a menos que se especifique lo contrario. Transporte reactivos (por ejemplo, enzimas) en neveras portátiles hasta y desde el almacenamiento a -20 °C para preservar su vida útil y evitar múltiples ciclos de congelación y descongelación de productos intermedios de biblioteca.

1. Desfosforilación de ARN

NOTA: La eliminación del 5'-fosfato (P; donante) impide la autoligadura al 3'-hidroxilo (OH; aceptor) de los ARN. Los enlazadores no se autoligarán, ya que su extremo 3' se modifica para incorporar una didesoxicitidina (Linker 1) o un espaciador (Linker 2). Los enlazadores solo se pueden ligar uniendo su 5'-P al 3'-OH de los ARN o ADNc.

  1. Preparar la reacción de desfosforilación en un tubo de PCR estéril (por ejemplo, tubo de 200 μL o 500 μL) añadiendo hasta 2 μL de la muestra de ARN, 0,5 μL de 40 U/μL de inhibidor de la ARNasa, 1 μL de patrón interno de 0,5 μM (Tabla 1), 0,5 μL de tampón de reacción de ARN ligasa T4 10x, 1 μL de 1 U/μL de fosfatasa alcalina de camarones, y suficiente agua libre de RNasa para llevar el volumen total a 5 μL.
    NOTA: Para muestras típicas, 75 ng de ARN (o aproximadamente 2 pmol para un ARN de 80 nt) se consideran suficientes para la cuantificación de ARN pequeños utilizando este protocolo.
  2. Incubar a 37 °C durante 30 min para desfosforilar los ARN y luego a 65 °C durante 5 min para inactivar la enzima y desnaturalizar los ARN. Mantenga las muestras a 4 °C durante al menos 10 minutos para evitar la renaturalización.

2. Ligadura del enlazador 1 al extremo 3' de los ARN

  1. Prepare la reacción de ligadura del Linker 1 en un tubo de PCR estéril añadiendo 5 μL de los ARN desfosforilados (productos del paso 1), 0,5 μL de 40 U/μL de inhibidor de la ARNasa, 1 μL de 100 μM de Linker 1 (Tabla 1), 3 μL de 10 mM de ATP, 2,5 μL de tampón de reacción de ARN ligasa T4 10x, 15 μL de PEG8000 (solución al 50%), 2 μL de 30 U/μL de ARN ligasa T4 1 y 1 μL de agua libre de ARNasa.
    NOTA: Los reactivos se pueden convertir en una mezcla maestra para facilitar el procesamiento de las muestras. No incluya PEG8000 y ARN ligasa 1 T4 en la mezcla maestra.
  2. Incubar a 25 °C durante 2 h y luego a 16 °C durante 16 h para ligar el Linker 1 a los ARN.
  3. Purificar en columna los ARN ligados a Linker-1. Utilice un kit para la recuperación y limpieza de ADN/ARN (consulte la tabla de materiales).
    NOTA: Este protocolo de purificación de columnas se aplica a todas las purificaciones de columnas posteriores que utilicen el mismo kit. Los kits que utilizan la tecnología de filtración en gel para eliminar los terminadores de tinte de las reacciones de secuenciación (ver Tabla de Materiales) no se pueden utilizar aquí, ya PEG8000 es incompatible con los filtros de gel de dichos kits. Se recomienda guardar una alícuota (1,2 μL) de cada muestra después de la purificación. Cuando sea necesario, estas alícuotas se pueden utilizar para comprobar la eficacia de la ligadura mediante la ejecución de un analizador de ácido nucleico comercial. Mantenga las muestras purificadas restantes en hielo o a -20 °C hasta que se realicen nuevos pasos.
    1. Para un volumen de muestra < 50 μL, añada agua sin RNasa hasta llegar a 50 μL. Añada 2 volúmenes del tampón de unión a oligonucleótidos (incluido en el kit) a 1 volumen de la muestra. Agregue 8 volúmenes de etanol al 100% a 1 volumen de la muestra.
    2. Cargue la muestra (hasta 750 μL a la vez) en una columna colocada dentro de un tubo de recolección de 2 mL (incluido en el kit). Para unir los ARN a la columna, centrifugar a 10.000 x g durante 30 s y desechar el flujo (es decir, el líquido en el tubo de recolección). Vuelva a colocar la columna en el tubo de recolección.
    3. Para lavar las impurezas de la columna, añada 750 μL del tampón de lavado de ADN (incluido en el kit) a la columna. Centrifugar a 10.000 x g durante 30 s y desechar el flujo. Vuelva a colocar la columna en el tubo de recolección.
    4. Centrifugar a velocidad máxima (por ejemplo, 16.000 x g para una microcentrífuga de sobremesa) durante 1 minuto adicional para eliminar el tampón de lavado de ADN residual.
    5. Para eluir los ARN, mueva con cuidado la columna a un tubo estéril de 1,5 mL, agregue agua libre de ARNasa a la columna. Utilice un volumen mayor del necesario al eluir los ARN para tener en cuenta la posible pérdida de volumen durante el proceso de elución. Por ejemplo, si se requieren 15 μL para el siguiente paso, agregue 17 μL de agua libre de RNasa para la elución. Centrífuga a 10.000 x g durante 30 s.

3. Eliminación de metilaciones postranscripcionales por AlkB demetilasa

NOTA: AlkB es una enzima bacteriana que elimina los grupos metilo de algunos, pero no todos, los nucleótidos metilados en el ADN y el ARN. La eliminación de varios tipos de ribonucleótidos metilados en el ARN evita la caída de la transcriptasa inversa para permitir lecturas más completas e identificación de sitios de modificación. Este paso debe controlarse a un pH bajo para evitar la degradación inesperada de los ARN.

  1. Prepare una solución madre de 1 M 2-cetoglutarato disolviendo 1,4611 g de 2-cetoglutarato (146,11 g por M) en 10 mL de agua libre de RNasa. Filtro: esterilice la solución a través de un filtro de jeringa de 0,2 μm. Aliquotar la solución madre en tubos estériles de 2 ml y almacenar a -20 °C.
  2. Prepare una solución madre de 0,5 M de ácido L-ascórbico disolviendo 0,88 g de ácido L-ascórbico (176,12 g por M) en 10 ml de agua libre de RNasa. Filtro: esterilice la solución a través de un filtro de jeringa de 0,2 μm. Aliquotar la solución madre en tubos estériles de 2 ml y almacenar a -20 °C.
  3. Preparar una solución madre de 0,25 M de sulfato ferroso amónico hexahidratado disolviendo 0,9835 g de sulfato ferroso amónico hexahidratado (392,14 g por M) en 10 mL de agua libre de RNasa. Filtro: esterilice la solución a través de un filtro de jeringa de 0,2 μm. Aliquotar la solución madre en tubos estériles de 2 ml y almacenar a -20 °C.
  4. Prepare una solución madre de 1 M de HEPES disolviendo 2,383 g de HEPES (238,30 g por M) en 10 mL de agua libre de RNasa. Ajuste el pH de la solución a 8 usando NaOH y esterilice la solución a través de un filtro de jeringa de 0,2 μm. Aliquotar la solución madre en tubos estériles de 2 ml y almacenar a -20 °C.
  5. Prepare un tampón de reacción AlkB 2x. Para producir 10 mL del tampón, combine 1,5 μL de 1 M de 2-cetoglutarato (elaborado en el paso 3.1), 80 μL de ácido L-ascórbico 0,5 M (elaborado en el paso 3.2), 6 μL de 0,25 M de sulfato ferroso de amonio hexahidratado (elaborado en el paso 3.3), 100 μL de 10 mg/mL de BSA, 1000 μL de 1 M HEPES (fabricado en el paso 3.4; añadir el último), y 8812,5 μL de agua libre de RNasa. Filtro: esterilice el tampón a través de un filtro de jeringa de 0,2 μm.
    NOTA: El tampón de reacción 2x AlkB debe prepararse fresco inmediatamente antes de cada experimento debido a la labilidad química de los componentes.
  6. Prepare la reacción de digestión de AlkB en un tubo de PCR estéril añadiendo 20 μL de los ARN ligados a Linker-1 (productos del paso 2), 50 μL del tampón de reacción 2x AlkB (fabricado en el paso 3.5), 2 μL de AlkB desmetilasa, 1 μL de inhibidor de RNasa y 27 μL de agua libre de RNasa.
  7. Incubar a temperatura ambiente durante 2 h para eliminar las metilaciones postranscripcionales de los ARN.
  8. Para eliminar AlkB de la reacción, siga los pasos que se describen a continuación.
    1. Para una separación de fases limpia, agregue 50 μL de agua libre de RNasa en la reacción AlkB y luego agregue 100 μL de fenol: cloroformo: alcohol isoamílico 25:24:1 (pH = 5,2).
    2. Agitar con la mano durante 10 s y, a continuación, centrifugar a 16.000 x g durante 10 min. Asegúrese de que el rotor de la centrífuga de sobremesa sea compatible con los tubos de PCR. Utilice adaptadores si es necesario.
    3. Transfiera los ARN (es decir, la capa acuosa en la parte superior; aproximadamente 140 μL) a un tubo estéril de 1,5 mL. Si el cloroformo (es decir, la capa inferior) se mezcla con la capa acuosa, centrifugar nuevamente con los mismos ajustes.
    4. Agregue 100 μL de cloroformo a los ARN extraídos para eliminar el fenol residual. Agitar con la mano durante 10 s y, a continuación, centrifugar a 16.000 x g durante 10 min.
    5. Transfiera los ARN (es decir, la capa acuosa en la parte superior; aproximadamente 120 μL) a un tubo estéril de 1,5 mL.
  9. Purificar en columna los ARN extraídos. Utilice un kit para la recuperación y limpieza de ADN/ARN (consulte la tabla de materiales). Siga el protocolo detallado en el paso 2.3 para realizar la purificación de la columna.

4. Eliminación del exceso de Linker 1

NOTA: Se recomienda guardar una alícuota (1,2 μL) de cada muestra después de la purificación. Cuando sea necesario, estas alícuotas se pueden utilizar para comprobar la eficiencia de la digestión de RecJf mediante el funcionamiento de un analizador de ácidos nucleicos comercial. Proceda con las muestras purificadas inmediatamente para revertir la transcripción.

  1. Prepare la reacción de indeadenilación en un tubo de PCR estéril añadiendo 15 μL de ARN ligados a Linker-1 (productos del paso 3), 1 μL de inhibidor de ARNasa de 40 U/μL, 2 μL de tampón 2 del kit 10x (consulte la Tabla de materiales) y 2 μL de 50 U/μL de 5'-deadenilasa.
  2. Incubar a 30 °C durante 1 h para eliminar la adenina en el extremo 5' del Linker 1. Añadir 2 μL de 30 U/μL de RecJf en la reacción de deadenilación.
  3. Incubar a 37 °C durante 30 min para digerir el exceso de Linker 1. Agregue otros 2 μL de 30 U/μL de RecJf a la reacción.
  4. Incubar a 37 °C durante 30 min para continuar la digestión del exceso de Linker 1 y luego a 65 °C durante 20 min para desnaturalizar la enzima.
  5. Purificar en columna los ARN ligados a Linker-1. Utilice un kit que utilice tecnología de filtración en gel para eliminar los terminadores de colorante de las reacciones de secuenciación (consulte la Tabla de materiales), ya que es eficaz para eliminar restos cortos (por ejemplo, oligonucleótidos con una longitud de 2 a 10 pb). Siga los pasos que se describen a continuación para la purificación.
    1. Prepare las columnas de filtro de gel de acuerdo con el protocolo del fabricante. Coloque una columna en un tubo estéril de 1,5 ml y cargue la columna con 24 μl de la muestra.
    2. Para purificar los ARN, centrifugar a 800 x g durante 3 min y desechar la columna. Los ARN purificados se encuentran en el eluyente.

5. Reacción de transcripción inversa (RT)

NOTA: La siguiente configuración de reacción RT (ver Tabla de Materiales) sigue el protocolo del fabricante, con pequeñas modificaciones para permitir la compatibilidad con AQRNA-seq.

  1. Prepare la reacción de recocido del cebador RT en un tubo de PCR estéril añadiendo 24 μL de ARN molde (productos del paso 4), 1 μL de cebador RT de 2 μM (Tabla 1) y 1 μL de dNTP (10 mM de cada tipo de nucleótidos).
  2. Incubar a 80 °C durante 2 min para recocer los cebadores RT a los ARN molde y, a continuación, enfriar en hielo inmediatamente durante 2 min.
  3. Prepare la reacción RT añadiendo 6 μL de tampón de reacción RT 5x, 1 μL de inhibidor de RNasa 40 U/μL y 1 μL de transcriptasa inversa en el tubo de reacción de recocido.
  4. Incubar a 50 °C durante 2 h para realizar la transcripción inversa de las plantillas de ARN, y luego a 70 °C durante 15 min para inactivar la enzima. Los productos RT (es decir, híbridos de ARN y ADNc) pueden almacenarse a 4 °C o -20 °C durante la noche.

6. Hidrólisis del ARN

  1. Añadir 1 μL de NaOH 5 M en el híbrido ARN-ADNc (productos del paso 5). Incubar a 93 °C durante 3 min para hidrolizar la cadena de ARN del híbrido ARN-ADNc.
  2. Añadir 0,77 μL de HCl 5 M para neutralizar la reacción. Después de agregar HCl, mueva para mezclar y gire por el tubo. La neutralización es instantánea.
    NOTA: Se recomienda probar la cantidad precisa de 5 M de HCl que se necesita para neutralizar 1 μL de NaOH (por ejemplo, usando tiras de pH) en condiciones tamponadas.
  3. Purificar en columna los ADNc monocatenarios. Utilice un kit para la recuperación y limpieza de ADN/ARN (consulte la tabla de materiales). Siga el protocolo detallado en el paso 2.3 para realizar la purificación de la columna.
    NOTA: Los kits que utilizan tecnología de filtración en gel para eliminar los terminadores de tinte de las reacciones de secuenciación (consulte la Tabla de materiales) no se pueden usar aquí debido a la variación del pH en los pasos anteriores.
  4. Vacíe rápidamente los ADNc purificados a < 5 μL y luego agregue agua libre de RNasa para que el volumen vuelva a 5 μL. Tenga cuidado de no vaciar rápidamente los ADNc hasta que se sequen por completo.
  5. Transfiera los ADNc purificados a un tubo de PCR estéril. Los ADNc purificados pueden almacenarse a -20 °C durante un máximo de 1 semana.

7. Ligadura del enlazador 2 hasta el extremo 3' de los ADNc

  1. Prepare la reacción de ligadura de Linker 2 en un tubo de PCR estéril agregando 5 μL de ADNc (productos del paso 6), 1 μL de 50 μM de Linker 2 (Tabla 1), 2 μL de tampón de reacción de ADN ligasa T4 10x, 1 μL de ATP 10 m, 9 μL de PEG8000 (solución al 50%) y 2 μL de 400 U/μL de ADN ligasa T4.
    NOTA: Los reactivos se pueden convertir en una mezcla maestra para facilitar el procesamiento de las muestras. No incluya PEG8000 y T4 ADN ligasa en la mezcla maestra.
  2. Incubar a 16 °C durante 16 h para ligar Linker 2 a los ADNc.
  3. Purificar en columna los ADNc ligados a Linker-2. Utilice un kit para la recuperación y limpieza de ADN/ARN (consulte la tabla de materiales). Siga el protocolo detallado en el paso 2.3 para realizar la purificación de la columna.

8. Eliminación del exceso de Linker 2

  1. Prepare la reacción de indeadenilación en un tubo de PCR estéril añadiendo 16 μL de ADNc ligados a Linker-2, 2 μL de tampón 2 del kit 10x y 2 μL de 50 U/μL de 5'-deadenilasa. Incubar a 30 °C durante 1 h para eliminar la adenina en el extremo 5' del Linker 2.
  2. Añadir 2 μL de 30 U/μL de RecJf en la reacción de deadenilación. Incubar a 37 °C durante 30 min para digerir el exceso de Linker 2. Agregue otros 2 μL de 30 U/μL de RecJf a la reacción.
  3. Incubar a 37 °C durante 30 min para continuar la digestión del exceso de Linker 2 y luego a 65 °C durante 20 min para desnaturalizar la enzima.

9. Amplificación por PCR de los ADNc con cebadores de secuenciación

  1. Asigne cebadores de PCR a las muestras. Cada muestra necesita una combinación única de cebadores directos e inversos (Tabla 1) para una multiplexación eficaz.
  2. Agregue agua sin RNasa para llevar el volumen de la muestra a 25 μL. Guarde 5 μL de cada muestra en un tubo de PCR estéril como respaldo en caso de que sea necesario repetir la PCR.
  3. Prepare la reacción de PCR (consulte la Tabla de materiales para el kit de PCR) agregando 20 μL de ADNc (productos del paso 8), 1 μL de cebador directo de 2,5 μM, 1 μL de cebador inverso de 2,5 μM, 25 μL de tampón de ADN polimerasa 2x, 2 μL de agua libre de ARNasa y 1 μL de ADN polimerasa.
    NOTA: Los reactivos se pueden convertir en una mezcla maestra para facilitar el procesamiento de las muestras. No agregue ADN polimerasa a la mezcla maestra.
  4. Realice la PCR con una desnaturalización inicial a 94 °C durante 1 min, seguida de 18 ciclos de desnaturalización a 98 °C durante 20 s - recocido a 58 °C durante 20 s - extensión a 68 °C durante 1 min.
    NOTA: No amplifique por PCR más allá del rango lineal. Un total de 18 ciclos será óptimo para la mayoría de los experimentos, pero esto puede depender del contexto.
  5. Vacíe los productos de PCR a menos de 25 μL y, a continuación, añada agua sin RNasa para que el volumen vuelva a ser de 25 μL. Transfiera 5 μL de los productos de PCR a un tubo estéril de 0,5 mL para comprobar la distribución del tamaño (consulte el paso 9.6). Almacene los 20 μL restantes de los productos PCR a -20 °C hasta que se realicen los siguientes pasos.
  6. Compruebe la distribución del tamaño de los productos PCR como se describe a continuación.
    1. Prepare gel de agarosa al 3% en tampón TAE.
      NOTA: En este protocolo, el bromuro de etidio (EtBr) se utiliza para la tinción en gel posterior a la electroforesis (consulte el paso 9.6.5). Las tinciones de ADN apropiadas se pueden agregar a la solución de gel en este paso o se pueden usar para la tinción en gel más adelante.
    2. Mezcle 1 μL de colorante de carga 6x en 5 μL de productos PCR (del paso 9.5) y cargue el gel con la mezcla.
    3. Cargue 5 μL de escaleras de ADN en el pocillo antes de la primera muestra y en el pocillo después de la última muestra. Utilice escaleras de ADN de 50 pb o 100 pb para permitir una mejor discriminación del tamaño de los productos de PCR entre 150 pb y 300 pb.
    4. Ejecute una electroforesis en gel para localizar los productos de PCR. La condición de ejecución adecuada puede depender del contexto. En este caso, funciona a 120 V, 400 mA durante 75 minutos para una losa de gel de 17,78 cm (ancho) x 10,16 cm (alto) x 1 cm (grosor).
    5. Coloque el gel en una caja y llénela con agua desionizada (DI) hasta que el gel esté completamente sumergido. Agregue 10 μL de EtBr en el agua desionizada empapando el gel. Envuelve la caja con papel de aluminio y colócala en una coctelera. Tiñe el gel durante 30 minutos mientras agitas.
    6. Deseche los residuos que contengan EtBr en una botella de residuos colocada en una campana extractora. Enjuague el gel con agua desionizada una vez y deseche el agua que contiene EtBr en la botella de desecho.
    7. Llene la caja con agua desionizada hasta que el gel esté completamente sumergido. Envuelve la caja con papel de aluminio y colócala en una coctelera. Lavar el gel durante 10 min mientras se agita.
    8. Deseche el agua que contiene EtBr en la botella de residuos en la campana extractora. Utilice un generador de imágenes en gel para visualizar las bandas. Adquiere una imagen de alta resolución del gel.

10. Purificación en gel

  1. Prepare gel de agarosa al 3% en tampón TAE. Haga un gel de 1 cm de espesor con peines anchos (1 mm de grosor; 5 mm de ancho; 15 mm de profundidad), de modo que cada pocillo pueda contener al menos 25 μL de mezcla de colorante de carga de muestra.
  2. Mezcle 4 μL de colorante de carga 6x con 20 μL de productos PCR (del paso 9.5) y cargue el gel con la mezcla. Deje carriles vacíos entre las muestras para minimizar la contaminación cruzada durante la escisión en gel.
  3. Cargue escaleras de ADN, ejecute electroforesis en gel, tiña y lave el gel, y tome imágenes de gel como se describe en el paso 9.6.
  4. Elimine los bloques de gel que contienen productos de PCR dentro del rango de tamaño objetivo. Para minimizar la contaminación con dímeros de cebadores (enlazadores de 175 pb sin insertos), extraiga productos de PCR con un tamaño superior a 195 pb (enlazadores de 175 pb + miRNAs de 20 pb).
  5. Purifique los productos de PCR mediante extracción en gel. Utilice un kit de extracción en gel (consulte la tabla de materiales). El protocolo de purificación se basa en el protocolo del fabricante, con pequeñas modificaciones para la compatibilidad con AQRNA-seq. Todos los pasos de centrifugación deben realizarse a 17.900 x g durante 1 min utilizando una centrífuga de sobremesa a temperatura ambiente, a menos que se especifique lo contrario. Siga los pasos que se describen a continuación.
    1. Mida el peso de los bloques de gel dentro de los tubos. Agregue 6 volúmenes de tampón QG (incluido en el kit) a 1 volumen de bloque de gel (1 mg de gel es aproximadamente 1 μL).
    2. Incubar a 50 °C durante 10 min o hasta la disolución completa de los bloques de gel. Tubos de vórtice cada 2 min para facilitar la disolución del gel. Después de disolver el gel, la mezcla debe parecerse al color del Buffer QG sin el gel disuelto. Si el color es naranja o violeta, agregue 10 μL de acetato de sodio 3 M (pH = 5.0) y mezcle bien.
    3. Agregue 1 volumen de gel de isopropanol a la mezcla y mezcle bien. Coloque una columna de centrifugado en un tubo de recolección de 2 ml (incluido en el kit).
    4. Para unir el ADN, aplique la muestra (hasta 750 μL cada vez) a la columna y centrifuga. Deseche el flujo y vuelva a colocar la columna en el mismo tubo de recolección. La cantidad máxima de gel por columna de centrifugado es de 400 mg.
    5. Añada 500 μL de tampón QG a la columna y centrifuga. Deseche el flujo y vuelva a colocar la columna en el mismo tubo de recolección.
    6. Para lavar las impurezas, agregue 750 μL de tampón PE (incluido en el kit) a la columna, deje reposar la columna durante 5 minutos y centrifuga. Deseche el flujo y vuelva a colocar la columna en el mismo tubo de recolección. Centrifugar una vez más para eliminar el tampón de lavado residual.
    7. Coloque la columna en un tubo estéril de 1,5 ml. Para eluir el ADN, agregue 30 μL de tampón EB (incluido en el kit) al centro de la membrana de la columna, deje reposar la columna durante 4 minutos y centrifuga.
    8. Acelere y vuelva a suspender los productos de PCR purificados en gel en 12 μL de Buffer EB.
  6. Mida la concentración de las bibliotecas construidas mediante espectrofotometría UV-Visible y/o cuantificación fluorométrica.

11. Secuenciación de bibliotecas

  1. Envíe las bibliotecas construidas a un centro de secuenciación externo para la evaluación de la calidad y la secuenciación de Illumina. Para garantizar una sensibilidad suficiente en el mapeo cuantitativo de pequeños paisajes de ARN, opte por la secuenciación de extremos emparejados con lecturas de 75 pb desde cada dirección (es decir, PE75), con el objetivo de obtener al menos 1,5 M de lecturas de secuencia bruta en cada dirección para cada muestra. Utilice cebadores personalizados (Tabla 1) para la secuenciación de NextSeq, pero esto es opcional para la secuenciación en MiSeq.
    NOTA: La secuenciación se puede realizar utilizando las plataformas MiSeq o NextSeq500. La elección de la plataforma puede depender de la naturaleza de las muestras y del recuento total de muestras.

12. Canalización de análisis de datos

NOTA: La Figura 2 proporciona una ilustración gráfica de los procedimientos simplificados involucrados en la canalización de análisis de datos, que toma lecturas de secuencia sin procesar (en formato FASTQ) como entrada y genera una matriz de abundancia con filas que representan miembros de pequeñas especies de ARN de interés y columnas que representan muestras. Para la secuenciación de extremos emparejados, cada muestra corresponde a dos archivos FASTQ, uno para las lecturas directas y otro para las lecturas inversas. El pipeline completo de analítica de datos, con todos los scripts asociados y un manual con extensas anotaciones para cada paso, está disponible en GitHub (https://github.com/Chenrx9293/AQRNA-seq-JoVE.git).

  1. Recupere lecturas de secuencias sin procesar del centro de secuenciación externo y evalúe la calidad de la secuenciación mediante programas de código abierto como FastQC14 o fastp15.
  2. Cree una biblioteca de secuencias de referencia en formato FASTA.
    NOTA: La clave para la adaptabilidad de la tubería a diversas clases pequeñas de ARN es una biblioteca de secuencias de referencia adecuada. Para lograr estimaciones precisas de la abundancia de miembros de clases específicas de ARN de interés (por ejemplo, miARN), se prevé que los usuarios seleccionen meticulosamente una biblioteca de secuencias de referencia para su uso con la tubería. Todas las demás instrucciones para la ejecución de la tubería permanecen consistentes en diferentes clases de ARN pequeñas.
  3. Cree un directorio denominado AQRNA-seq para implementar la canalización de análisis de datos y coloque todos los scripts, las lecturas de secuencias filtradas por calidad y la biblioteca de secuencias de referencia en este directorio. Siga las instrucciones detalladas en GitHub (https://github.com/Chenrx9293/AQRNA-seq-JoVE.git) para preparar subdirectorios y realizar modificaciones esenciales en los archivos para dirigirse a diferentes organismos y/o pequeñas especies de ARN, así como la compatibilidad del sistema operativo y/o el programador de trabajos.
  4. Implemente la canalización de análisis de datos siguiendo los pasos descritos en el manual en GitHub (https://github.com/Chenrx9293/AQRNA-seq-JoVE.git), que contiene descripciones, archivos de entrada y salida, así como líneas de comandos para cada paso. En resumen, la canalización abarca (i) el recorte de las secuencias enlazadoras y los nucleótidos aleatorios de las lecturas, (ii) el filtrado de las lecturas en función de su longitud, (iii) el mapeo de las lecturas a las secuencias de referencia, (iv) la resolución de asignaciones ambiguas y (v) la generación de la matriz de abundancia.

Resultados

Loading...
$$\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.

figure-results-1
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.

figure-results-2
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.

figure-results-3
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.

figure-results-4
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.

figure-results-5
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.

figure-results-6
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.

figure-results-7
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.

Discusión

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

El flujo de trabajo de preparación de la biblioteca AQRNA-seq está diseñado para maximizar la captura de ARN dentro de una muestra y minimizar la caída de la polimerasa durante la transcripción inversa7. A través de una ligadura de enlazador de dos pasos, los nuevos oligonucleótidos de ADN (Linker 1 y Linker 2) se ligan en exceso para complementar completamente el ARN dentro de la muestra. El exceso de enlazadores se puede eliminar de manera eficiente con RecJf, una exonucleasa de 5' a 3' específica para ADN monocatenario, dejando intactos los productos ligados. Además, el tratamiento con AlkB reduce las modificaciones de metilo en los ARN18, que pueden causar la caída de la polimerasa durante la transcripción inversa y, por lo tanto, alivia el problema de los productos de ADNc truncados artefactuales. AQRNA-seq representa el primer método reportado para optimizar sistemáticamente la eficiencia de la ligadura y la amplificación. AQRNA-seq ha sido validado por su precisión cuantitativa y ausencia de artefactos de sesgo, y se ha demostrado que es la aplicación más precisa para los estudios de perfiles de miARN.

AQRNA-seq es notablemente extenso, y se recomienda el cuidado operativo para los pasos sensibles a la temperatura y/o al tiempo durante la preparación de la biblioteca. En concreto, los pasos previos a la síntesis de ADNc (Figura 1; pasos 1-4), en los que las muestras contienen ARN monocatenario (ARNs) inherentemente inestable o ARNs ligado a un enlazador de ADN, deben llevarse a cabo sin interrupciones innecesarias. La hidrólisis del ARN (Figura 1; paso 6) utilizando una base fuerte a alta temperatura es muy eficiente en la digestión de plantillas de ARN después de la transcripción inversa19. Sin embargo, se debe tener cuidado para evitar la exposición prolongada de la muestra a la base fuerte y la neutralización tardía con ácido fuerte, especialmente dado un gran recuento de muestras. En general, se recomienda encarecidamente preparar todas las reacciones en hielo y evitar ciclos repetidos de congelación y descongelación para los productos intermedios de la biblioteca.

Se recomienda el muestreo de los intermedios de la biblioteca en puntos clave a lo largo del flujo de trabajo para las comprobaciones de calidad y la resolución oportuna de problemas del proceso. El análisis de muestras en herramientas de electroforesis automatizadas puede determinar el éxito y la eficiencia de la ligadura del enlazador (Figura 1; pasos 2 y 7). La visualización de la localización de los productos de ADNc indexados después de la adición de cebadores de secuenciación por PCR (Figura 1; paso 9) a través de la electroforesis en gel de agarosa puede ayudar a evaluar la idoneidad de los ciclos de PCR, así como a identificar la posible degradación de las muestras y/o la presencia de dímeros de cebadores. Los cebadores de PCR y secuenciación pueden rediseñarse para que sean compatibles con la plataforma de secuenciación de interés del operador.

Debido en parte a los distintos diseños involucrados en la preparación de la biblioteca, la canalización de análisis de datos AQRNA-seq se construyó principalmente con scripts personalizados, incorporando solo un puñado de programas bioinformáticos existentes. Sin embargo, para alinear las lecturas de secuencia con la biblioteca de secuencias de referencia, se empleó la Herramienta Básica de Búsqueda de Alineación Local (BLAST)20 debido a sus diversas fortalezas que satisfacen las necesidades específicas del análisis de datos AQRNA-seq. Estos incluyen (i) la capacidad de alinear lecturas de menos de 20 pb y realizar alineaciones de subcadenas, (ii) el manejo sensato de bases ambiguas (N), (iii) la precisión de alineación ligeramente mayor y (iv) la salida informativa que proporciona detalles sustanciales como discrepancias y brechas, lo que es propicio para diversos análisis posteriores.

En particular, siguen existiendo varios desafíos de diseño tanto en el flujo de trabajo de preparación de bibliotecas como en la canalización de análisis de datos de AQRNA-seq, que justifican optimizaciones adicionales y representan el enfoque del trabajo actual en curso. En primer lugar, la precisión cuantitativa superior de AQRNA-seq en comparación con los métodos convencionales de RNA-seq se logró con 50-75 ng de entrada de ARN pequeño, mientras que el flujo de trabajo se prueba actualmente con cantidades reducidas de ARN de entrada para mejorar su sensibilidad y con ARN más largos para ampliar su utilidad. En segundo lugar, aunque la metilación postranscripcional se elimina enzimáticamente con AlkB, existen otras modificaciones del ARN que pueden causar la caída de la polimerasa durante la transcripción inversa10. Esto puede conducir a productos de ADNc truncados artifactuales que no podrían diferenciarse de manera efectiva de la degradación 5' biológicamente significativa de ARNt de longitud completa. Por lo tanto, actualmente se están evaluando varias formas posibles de mejorar aún más la procesividad de la transcriptasa inversa para determinar su eficacia y compatibilidad con AQRNA-seq. En tercer lugar, la eliminación de los dímeros de cebadores a través de la extracción y purificación en gel es ineficiente e invariablemente resulta en una pérdida de los fragmentos objetivo y una mayor variabilidad entre las bibliotecas. Por lo tanto, mitigar el arrastre potencial de dímeros de cebadores sin perder fragmentos objetivo representa una oportunidad para mejorar aún más la sensibilidad y la precisión cuantitativa de AQRNA-seq. Es importante destacar que tanto el tratamiento con AlkB como la extracción y purificación del gel requieren mucha mano de obra, lo que los hace subóptimos para procesar un gran número de muestras. Con las mejoras logradas en las áreas mencionadas, la automatización y el procesamiento de alto rendimiento del método AQRNA-seq se habilitarán en el futuro. El refinamiento continuo de la canalización de análisis de datos prioriza las mejoras en la eficiencia, la precisión cuantitativa y la flexibilidad. Si bien la tubería actual recorta efectivamente los adaptadores 3' de las lecturas de secuencia, un subconjunto de lecturas también posee secuencias de adaptadores 5', probablemente debido al bloqueo inadecuado del extremo 3' de los enlazadores durante la ligadura. La presencia de secuencias de adaptadores 5' puede comprometer la precisión cuantitativa al interferir con la alineación de la lectura de secuencias. Por lo tanto, se garantizará una eliminación adecuada de las secuencias de adaptadores de 3' y 5' mediante el empleo de herramientas de recorte de adaptadores adicionales como Cutadapt21. Para reforzar la eficiencia, la cartera revisada integrará estrategias racionales para reducir el tiempo de cómputo. Esto implica un paso de ensamblado de lectura final emparejado para fusionar lecturas hacia adelante e inversas mediante la identificación de una superposición significativa y la extracción de secuencias únicas junto con sus ocurrencias en cada biblioteca. Se prevé que estas estrategias reduzcan considerablemente el tiempo de cálculo descendente al evitar el análisis redundante de (i) lecturas desde ambas direcciones y (ii) secuencias duplicadas que pueden dominar el conjunto de lecturas de secuencias de las bibliotecas AQRNA-seq. En el futuro, la optimización sinérgica tanto del flujo de trabajo de preparación de bibliotecas como de la canalización de análisis de datos de AQRNA-seq probablemente abrirá diversas oportunidades en la investigación cuantitativa rigurosa de todas las formas de ARN (por ejemplo, transcriptoma, fragmentos de ARNt y especies raras de ARN como los ARN tumorales circulantes), el mapeo de las modificaciones del ARNt y otras áreas fundamentales de la biología de biomarcadores y descubrimientos.

Divulgaciones

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

P.C.D. es inventor de dos patentes (PCT/US2019/013714, US 2019/0284624 A1) relacionadas con la obra publicada.

Agradecimientos

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Los autores del presente trabajo agradecen a los autores del artículo original que describe la tecnología AQRNA-seq7. Este trabajo fue apoyado por subvenciones de los Institutos Nacionales de Salud (ES002109, AG063341, ES031576, ES031529, ES026856) y la Fundación Nacional de Investigación de Singapur a través de la Alianza Singapur-MIT para la Investigación y Tecnología de la Resistencia a los Antimicrobianos IRG.

Materiales

Lista de materiales utilizados en este artículo
NombreEmpresaNúmero de catálogoComentarios
2-cetoglutarato Sigma-Aldrich75890Prepare una solución de trabajo (1 M) y guárdela a -20º C
2100 Instrumento bioanalizadorAgilentG2938C
Adenosina 5'-trifosfato (ATP) New England BiolabsM0437M (componente #: N0437AVIAL)NEB M0437M contiene ARN ligasa T4 1 (30 U/μ L), tampón de reacción de ARN ligasa T4 (10X), PEG 8000 (1X) y ATP (100 mM); preparar una solución de trabajo (10 mM) y almacenarla a -20º C
AGAROSA GPG/LEAmericanBioAB00972-00500Almacenar a temperatura ambiente
Sulfato de hierro (II) hexahidratadoSigma-AldrichF2262Prepare una solución de trabajo (0,25 M) y guárdela a -20 °C; C
Bioanalyzer Análisis de ARN pequeñoAgilent5067-1548 El análisis de ARN pequeño se utiliza para comprobar la calidad de los ARN de entrada y la eficiencia de las reacciones enzimáticas (por ejemplo, ligadura de Linker 1)
Albúmina sérica bovina (BSA; 10 mg/mL) New England BiolabsB9000Este producto se suspendió el 15/12/2022 y se reemplaza con albúmina recombinante, grado de biología molecular (NEB B9200).
Cloroformo Macron Fine Chemicals4441-10
Desmetilasa ArrayStarAS-FS-004Demethylase viene con el pretratamiento de ARNt rtStar & Kit de síntesis de ADNc de primera cadena (AS-FS-004)
Mezcla de solución de desoxinucleótido (dNTP)New England BiolabsN0447L (componente #: N0447LVIAL)Esta mezcla de soluciones dNTP contiene concentraciones equimolares de dATP, dCTP, dGTP y dTTP (10 mM cada una)
Bloque de calor doble digitalVWR Scientific Products13259-052El bloque calefactor se utiliza con el QIAquick  Kit de extracción de gel 
DyeEx  2.0 Kit de giro Qiagen63204 Eficaz para eliminar restos cortos (por ejemplo, oligonucleótidos de menos de 10 pb de longitud)
Fuente de alimentación para electroforesisBio-Rad LabrotoriesPowerPac 300
Eppendorf Tubos PCR (0,5 mL)Eppendorf0030124537
Eppendorf Tubos Safe-Lock (0,5 mL)Eppendorf022363611
Eppendorf Tubos Safe-Lock (1,5 mL)Eppendorf022363204
Eppendorf Tubos Safe-Lock (2 mL)Eppendorf022363352
Alcohol etílico (etanol),Sigma-AldrichE7023 El etanol puro se utiliza con el kit de limpieza y concentrador Oligo de Zymo Research
Sistema de imágenes en gelAlpha InnotechFluorChem 8900
Tinte de carga en gel, púrpura (6X), sin SDSNew England BiolabsN0556S (componente #: B7025SVIAL)NEB N0556S contiene Escalera de ADN púrpura de carga rápida de 50 pb y tinte de carga en gel, púrpura (6X), sin SDS
Espectrofotómetro GENESYS 180 UV-VisThermo Fisher Scientific840-309000El espectrofotómetro se utiliza para medir las concentraciones de oligonucleótidos utilizando la ley de Beer
HEPESSigma-AldrichH4034Prepare una solución de trabajo (1 M; pH = 8 con NaOH) y guárdela a -20 ° C
Ácido clorhídrico (HCl) Productos científicos de VWRBDH3028 Prepare una solución de trabajo (5 M) y guárdela a temperatura ambiente
Alcohol isopropílico (isopropanol), puroMacron Fine Chemicals3032-16El isopropanol se utiliza con el QIAquick  Kit de extracción de gel 
Ácido L-ascórbicoSigma-AldrichA5960Prepare una solución de trabajo (0,5 M) y guárdela a -20º
C MicrocentrífugaEppendorf5415D
NanoDrop 2000 EspectrofotómetroThermo Fisher ScientificND-2000
NEBuffer 2 (10X)New England BiolabsM0264L (componente #: B7002SVIAL)NEB M0264L contiene RecJf (30 U/μ L) y NEBuffer 2 (10X); tienda a -20 ° C
Agua sin nucleasas (no tratada con DEPC)Thermo Fisher ScientificAM9938 
Oligo Limpio & Kit concentrador Investigación ZymoD4061 Almacenar a temperatura ambiente
PEG 8000 (solución al 50%) New England BiolabsM0437M (componente #: B1004SVIAL)NEB M0437M contiene ARN ligasa T4 1 (30 U/μ L), tampón de reacción de ARN ligasa T4 (10X), PEG 8000 (1X) y ATP (100 mM); preparar una solución de trabajo (10 mM) y almacenarla a -20º C
Termociclador PeltierMJ ResearchPTC-200
Fenol: cloroformo: isoamilo   alcohol 25:24:1 pH = 5.2 Thermo Fisher ScientificJ62336 
Búfer PrimeScript (5X)TaKaRa2680A 
Transcriptasa inversa PrimeScriptTaKaRa2680A 
QIAquick  Kit de extracción de gel Qiagen28704Este kit requiere un bloque calefactor e isopropanol para trabajar con
Quick-Load Purple 100 bp DNA LadderNew England BiolabsN0551S (componente #: N0551SVIAL)
Quick-Load Purple 50 bp DNA LadderNew England BiolabsN0556S (componente #: N0556SVIAL)NEB N0556S contiene Quick-Load Purple 50 bp DNA Ladder y Gel Loading Dye, Purple (6X), sin SDS
RecJf  (30 U/μ L) New England BiolabsM0264L (componente #: M0264LVIAL)NEB M0264L contiene RecJf (30 U/μ L) y NEBuffer 2 (10X); tienda a -20 °
C Inhibidor de la ARNasa (murino; 40 U/μ L) New England BiolabsM0314L ( componente #: M0314LVIAL)Almacenar a -20 ° C
SeqAMP  ADN polimerasa TaKaRa638509 TaKaRa 638509 contiene SeqAMP  ADN polimerasa  y SeqAMP  Tampón de PCR (2X)
SeqAMP  Tampón de PCR (2X) TaKaRa638509 TaKaRa 638509 contiene SeqAMP  ADN polimerasa  y SeqAMP  Tampón PCR (2X)
Fosfatasa alcalina de camarones (1 U/μ L) New England BiolabsM0371L ( componente #: M0371LVIAL)
Hidróxido de sodio (NaOH)Sigma-Aldrich  S5881Prepare una solución de trabajo (5 M) y guárdela a temperatura ambiente
T4 ADN Ligasa (400 U/μ L) New England BiolabsM0202L ( componente #: M0202LVIAL)NEB M0202L contiene ADN ligasa T4 (400 U/μ L) y Tampón de reacción de ADN ligasa T4 (10X)
Tampón de reacción de ADN ligasa T4 (10X) New England BiolabsM0202L ( componente #: B0202SVIAL)NEB M0202L contiene ADN ligasa T4 (400 U/μ L) y tampón de reacción de ADN ligasa T4 (10X)
T4 ARN ligasa 1 (30 U/μ L) New England BiolabsM0437M (componente #: M0437MVIAL)NEB M0437M contiene T4 ARN ligasa 1 (30 U/μ L), tampón de reacción de ARN ligasa T4 (10X), PEG 8000 (1X) y ATP (100 mM)
Tampón de reacción de ARN ligasa T4 (10X) New England BiolabsM0437M (componente #: B0216SVIAL)NEB M0437M contiene ARN ligasa T4 1 (30 U/μ L), tampón de reacción de ARN ligasa T4 (10X), PEG 8000 (1X) y ATP (100 mM)
puro

Referencias

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,
  1. Byron, S. A., Van Keuren-Jensen, K. R., Engelthaler, D. M., Carpten, J. D., Craig, D. W. Translating RNA sequencing into clinical diagnostics: opportunities and challenges. Nat Rev Gen. 17 (5), 257-271 (2016).
  2. Grillone, K., et al. Non-coding RNAs in cancer: platforms and strategies for investigating the genomic "dark matter.". J Exp Clin Cancer Res. 39 (1), 117(2020).
  3. Hwang, B., Lee, J. H., Bang, D. Single-cell RNA sequencing technologies and bioinformatics pipelines. Exp Mol Med. 50 (8), 1-14 (2018).
  4. Goh, J. J. L., et al. Highly specific multiplexed RNA imaging in tissues with split-FISH. Nat Methods. 17 (7), 689-693 (2020).
  5. Moses, L., Pachter, L. Museum of spatial transcriptomics. Nat Methods. 19 (5), 534-546 (2022).
  6. Cummings, B. B., et al. Improving genetic diagnosis in Mendelian disease with transcriptome sequencing. Sci Transl Med. 9 (386), 5209(2017).
  7. Hu, J. F., et al. Quantitative mapping of the cellular small RNA landscape with AQRNA-seq. Nat Biotech. 39 (8), 978-988 (2021).
  8. Alon, S., et al. Barcoding bias in high-throughput multiplex sequencing of miRNA. Genome Res. 21 (9), 1506-1511 (2011).
  9. Fuchs, R. T., Sun, Z., Zhuang, F., Robb, G. B. Bias in ligation-based small RNA sequencing library construction is determined by adaptor and RNA structure. PLoS One. 10 (5), e0126049(2015).
  10. Pang, Y. L. J., Abo, R., Levine, S. S., Dedon, P. C. Diverse cell stresses induce unique patterns of tRNA up- and down-regulation: tRNA-seq for quantifying changes in tRNA copy number. Nuc Acids Res. 42 (22), e170(2014).
  11. Machnicka, M. A., Olchowik, A., Grosjean, H., Bujnicki, J. M. Distribution and frequencies of post-transcriptional modifications in tRNAs. RNA Biol. 11 (12), 1619-1629 (2014).
  12. Li, F., et al. Regulatory impact of RNA secondary structure across the Arabidopsis transcriptome. Plant Cell. 24 (11), 4346-4359 (2012).
  13. García-Nieto, P. E., Wang, B., Fraser, H. B. Transcriptome diversity is a systematic source of variation in RNA-sequencing data. PLOS Comput Biol. 18 (3), e1009939(2022).
  14. FASTQC: a quality control tool for high throughput sequence data. , Available from: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (2010).
  15. Chen, S., Zhou, Y., Chen, Y., Gu, J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 34 (17), i884-i890 (2018).
  16. Love, M. I., Huber, W., Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15 (12), 550(2014).
  17. R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing. , Vienna, Austria. (2022).
  18. Ougland, R., et al. AlkB restores the biological function of mRNA and tRNA inactivated by chemical methylation. Mol Cell. 16 (1), 107-116 (2004).
  19. Bernhardt, H. S., Tate, W. P. Primordial soup or vinaigrette: did the RNA world evolve at acidic pH. Biol Direct. 7, 4(2012).
  20. Altschul, S. F., Gish, W., Miller, W., Myers, E. W., Lipman, D. J. Basic local alignment search tool. J Mol Biol. 215 (3), 403-410 (1990).
  21. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 17 (1), 10-12 (2011).

Reimpresiones y permisos

Solicitar permiso para reutilizar el texto o las figuras de este artículo de JoVE

Solicitar permiso

Etiquetas

AQRNA SeqSmall RNA QuantificationtRNA ModificationsLibrary PreparationReverse TranscriptionGel ElectrophoresisBioinformatics PipelineProtein TranslationMicroRNA MappingPhenol Chloroform Extraction

Artículos relacionados