$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
La disponibilidad de secuenciación profunda rápida, confiable y rentable ha revolucionado muchos aspectos en el campo de Ciencias de la vida, permitiendo a gran profundidad en los análisis de secuenciación. Un desafío restante radica en el diseño innovador y la creación del representante de bibliotecas de la secuencia. Aquí se describe un protocolo para capturar las moléculas de ADNc viral naciente, específicamente los intermedios del proceso de transcripción inversa de VIH-1.
El paso más crítico de esta estrategia es la ligadura de un adaptador para el open 3'-termini de manera cuantitativa y objetiva. Eficiencias de trompas entre dos termini de ssDNA, tanto inter- y intramolecular, investigaron y optimizado para diversas aplicaciones11,26,27,28,29. La opción de utilizar un adaptador de horquilla con T4 ADN ligasa en las condiciones descritas en el paso 3.3 es el resultado de la optimización empírica en la que se evaluaron diferentes ligasas, adaptadores y reactivos para la ligación de oligonucleótidos sintéticos que representan Secuencias de HIV-1 (tabla 2) (datos no mostrados). En estas reacciones de la prueba en vitro , confirmamos que el ligase de la DNA de T4 mediada por la ligadura del adaptador de la horquilla, según lo descrito por Kwok et al. 11, tiene un sesgo muy bajo y logra cerca ligadura completa de moléculas aceptoras cuando el adaptador se usa en exceso. La eficacia de la ligadura era inafectada por la adición de la secuencia de nucleótidos para representar el adaptador compatible para el sistema de primer multiplex (ver figura 4). En comparación, encontramos que un termoestable 5' ADN/ARN ligasa («Ligasa A», véase Tabla de materiales para exactos ligasas comparado aquí), que es una ingeniería ligasa de RNA que se desarrolló en parte para mejorar la eficacia de la ligadura con ssDNA como aceptor 27, fue más efectivo que ligan dos moléculas de ssDNA de ligasa de RNA "(B de la ligasa del) pero tenía un sesgo importante, con fuertes diferencias en la eficacia de la ligadura entre oligonucleótidos con diferencias de longitud de base individual [tabla 2 ; HTP con media G (a) y (b)]. Además, encontramos sólo un sesgo mínimo en reacciones con "Ligasa C" combinado con un adaptador que lleva al azar 5'-termini (una estrategia utilizada para compensar el sesgo de nucleótidos conocido "C ligasa"; véase por ejemplo Ding et al. 30). sin embargo, la "C de ligasa"-mediadas trompas intermoleculares estaban incompletos, haciendo el sistema de T4 ADN ligasa la opción superior.
Varias medidas de control de calidad sobre el curso de protocolo y la inclusión de controles positivos y negativos permiten para la detección de problemas potenciales antes de continuación de ensayo y proporcionan orientación para la solución de problemas de esfuerzos. Las cuantificaciones de qPCR en pasos 2.2.2 2.3.12 aseguran de que la cantidad de material de entrada es suficiente. Números de copia cDNA típica en el rango de elución (del paso 2.1) 200 μL de alrededor de 10.000 a 300.000 por μl. El paso de captura de híbrido puede resultar en una pérdida general de VIH-1 cantidad de cDNA, pero dará lugar a un fuerte enriquecimiento de cDNA específico del VIH-1 sobre el ADN celular, que puede determinarse mediante el uso de primers adecuados para cuantificar ADN genómico antes y después de enriquecimiento por qPCR o midiendo la concentración de DNA total. CDNA de VIH-1 recuperado después de que el híbrido captura pasos debe ser al menos el 10% de la entrada. Baja a partir de material de lo contrario puede explicar un control positivo de éxito del oligonucleótido (ver paso 3.3.2) pero limitado sólo lee en las muestras. Bajo Lee números general también podrían explicarse por la sobreestimación de la concentración de la biblioteca debido a la presencia de especies de ADN irrelevantes sin adaptadores MiSeq. Esto daría lugar a densidad baja del cluster y puede mejorarse mediante la determinación de la concentración de secuencias de HIV-1 en la biblioteca por qPCR además de la cantidad de ADN total por análisis fluorométrico. Debido a la naturaleza altamente sensible del método, se debe tener especial cuidado para evitar la contaminación incluso de bajo nivel, tanto de otras muestras (en particular, de las existencias de oligonucleótidos de alta concentración control) así como de equipos de laboratorio. Trabajando en la estación de trabajo PCR de esterilización UV es beneficioso en este sentido. La electroforesis automatizada de la biblioteca final (paso 6.1.2) es otra medida de control de calidad. El rango de tamaño de ácidos nucleicos observado por lo general es entre 150 a 500 nt. imprimaciones que pueden detectarse en el control opcional después de la polimerización en cadena y antes de purificación (ver nota en el paso 5.2) ahora esté ausente. En un resultado representativo, la curva de intensidad de la muestra tiene un pico alrededor de 160 a 170 nt y un segundo más agudo pico alrededor de 320 a 350 nt. Esto probablemente refleja la abundancia superior ve a menudo en la transcripción inversa de relativamente corto (longitud de inserción de 1 a 20 nt) y larga duración fuerte parada (longitud de inserción de 180 a 182 nt) (figura 3b).
Mientras que el protocolo presentado y las cartillas son específicas para construcciones de transcripción reversa de HIV-1 tempranas, el método es generalmente aplicable a cualquier estudio con el objetivo de determinar abierto 3' termini de ADN. Las principales modificaciones en otros contextos será el método de captura híbrida y la estrategia de diseño de la cartilla. Por ejemplo, si el objetivo es adaptarse a las transcripciones finales de VIH-1, un mayor número de diferentes oligonucleótidos biotinilados captura recocido a lo largo del cDNA sería recomendable y probablemente disminuirá la pérdida en el paso de captura de híbridos. Como se mencionó en la introducción, es importante considerar las limitaciones en el diseño de la gama que 3'-termini está a detectar para evitar diversas fuentes de sesgo. En primer lugar, puede haber un sesgo en las reacciones de polimerización en cadena si las plantillas con el adaptador son de longitud muy variable. En segundo lugar, la plataforma la secuencia utilizada aquí (por ejemplo, MiSeq) tiene un rango de longitud de inserción preferido para el agrupamiento óptimo y significativamente productos más cortos y más largo no pueden ser ordenados con la misma eficiencia. En parte, esto puede abordarse computacionalmente, como se hizo mediante el cálculo de un factor de corrección por sesgo de longitud (ver figura 4, gráfico inferior). Sin embargo, si la región de donde se desea la 3'-termini mapeo es larga (> 1000 nt), es más aconsejable dividir las reacciones con las transcripciones ligadas y usar múltiples cebadores aguas arriba para evaluar 3'-termini en secciones.
El programa de análisis fue escrito en la empresa con el propósito específico de analizar ambos el último nucleótido de la secuencia del VIH-1 adyacente a la secuencia del adaptador fijo, así como la variación de la base de todas las bases para identificar cualquier mutaciones. Los pasos individuales comprenden las siguientes: en primer lugar, las secuencias de adaptador se recortan utilizando el toolkit de fastx-0.0.13; Luego, se eliminan cualquier secuencia que está duplicado (es decir, idénticas secuencias incluyendo el código de barras). Todas las lecturas únicas restantes entonces se alinean con la secuencia del VIH-1 con Bowtie (http://bowtie-bio.sourceforge.net/index.shtml) el desajuste máximo en tres bases. La secuencia de la plantilla está conformada por la primera 635 nt del cDNA de VIH-1 (cepa NL4.3), que incluye la secuencia - sss y el primer producto de transferencia de filamento hasta la pista de polypurine (U5-R-U3-PPT; ver figura 1). Por lo tanto, el software suministrado y plantillas aptas sólo directamente si el método se utiliza para la misma aplicación (detección de tempranas transcripciones inversas de la de VIH-1NL4.3). Ajustes tendrá que hacer para otras secuencias diana. Las posiciones de la 3'-termini para cada lectura se determinaron por la posición en la alineación. Se registran llamadas base para cada posición y se calculan tasas de mutación de la cobertura total de cada base, que varía, como Lee es de diferentes longitudes e insertos largos podrían no estar cubiertos completamente por la secuencia de la base de 125 en Read2.
Para concluir, creemos que el método descrito para ser una herramienta valiosa para muchos tipos de estudios. Aplicaciones obvias incluyen las investigaciones de los mecanismos de inhibición de reversa de la transcripción a través de medicamentos antirretrovirales o factores de restricción celular. Sin embargo, sólo relativamente pequeños ajustes deben ser necesarios adaptar el sistema a 3'-termini mapeo dentro de otros intermediarios virales monocatenario DNA, que están presentes, por ejemplo, en la replicación del parvovirus. Además, el principio del método, particularmente su paso ligadura optimizado, puede proporcionar una parte fundamental del diseño de preparación de biblioteca para la caracterización de las extensiones de ADN 3', incluyendo elongaciones catalizadas por ADN de doble cadena celular polimerasas.