En genética, la secuenciación de escopeta o shotgun sequencing es un método utilizado para secuenciar cadenas de ADN aleatorias. Recibe su nombre por analogía con el patrón de disparo cuasialeatorio de rápida expansión de una escopeta.
El método de secuenciación de ADN por terminación de cadena ("secuenciación de Sanger") solo se puede utilizar para cadenas cortas de ADN de 100 a 1000 pares de bases. Debido a este límite de tamaño, las secuencias más largas se subdividen en fragmentos más pequeños que se pueden secuenciar por separado, y estas secuencias se ensamblan para dar la secuencia general.
Hay dos métodos principales para este proceso de fragmentación y secuenciación. Primer walking (o "chromosome walking") progresa a través de toda la hebra pieza por pieza, mientras que la secuenciación rápida es un proceso más rápido, pero más complejo que utiliza fragmentos aleatorios.
En la secuenciación de escopeta, ADN se divide al azar en numerosos segmentos pequeños, que se secuencian utilizando el método de terminación de cadena para obtener lecturas. Se obtienen múltiples lecturas superpuestas para el ADN diana realizando varias rondas de esta fragmentación y secuenciación. Los programas de computadora luego usan los extremos superpuestos de diferentes lecturas para ensamblarlos en una secuencia continua.
La secuenciación por escopeta fue una de las tecnologías precursoras que fue responsable de permitir la secuenciación completa del genoma.
Ejemplo Por ejemplo, considere las siguientes dos rondas de lecturas de escopeta: Hebra Secuencia Original AGCATGCTGCAGTCATGCTTAGGCTA Primera secuencia de escopeta AGCATGCTGCAGTCATGCT------- -------------------TAGGCTA Segunda secuencia de escopeta AGCATG-------------------- ------CTGCAGTCATGCTTAGGCTA Reconstrucción AGCATGCTGCAGTCATGCTTAGGCTA En este ejemplo extremadamente simplificado, ninguna de las lecturas cubre la longitud completa de la secuencia original, pero las cuatro lecturas se pueden ensamblar en la secuencia original usando la superposición de sus extremos para alinearlas y ordenarlas. En realidad, este proceso utiliza enormes cantidades de información plagada de ambigüedades y errores de secuenciación. El ensamblaje de genomas complejos también se complica por la gran abundancia de secuencias repetitivas, lo que significa que lecturas cortas similares podrían provenir de partes completamente diferentes de la secuencia.
Se necesitan muchas lecturas superpuestas para cada segmento del ADN original para superar estas dificultades y ensamblar con precisión la secuencia. Por ejemplo, para completar el Proyecto Genoma Humano, la mayor parte del genoma humano se secuenció a una cobertura 12X o mayor; es decir, cada base de la secuencia final estuvo presente en promedio en 12 lecturas diferentes. Aun así, los métodos actuales no han logrado aislar o ensamblar una secuencia confiable para aproximadamente el 1% del genoma humano (eucromático), hasta 2004.
Secuenciación escopeta del genoma completo
Historia La secuenciación rápida del genoma completo para genomas pequeños (4000 a 7000 pares de bases) se sugirió por primera vez en 1979. El primer genoma secuenciado por secuenciación de escopeta fue el del virus del mosaico de la coliflor, publicado en 1981.
Secuenciación paired-end La aplicación más amplia se benefició de la secuenciación de extremos por pares, conocida coloquialmente como secuenciación de escopeta de doble cañón (double-barrel shotgun sequencing). A medida que los proyectos de secuenciación comenzaron a asumir secuencias de ADN más largas y complicadas, varios grupos comenzaron a darse cuenta de que se podía obtener información útil secuenciando ambos extremos de un fragmento de ADN. Aunque secuenciar ambos extremos del mismo fragmento y realizar un seguimiento de los datos emparejados fue más engorroso que secuenciar un solo extremo de dos fragmentos distintos, saber que las dos secuencias estaban orientadas en direcciones opuestas y tenían aproximadamente la longitud de un fragmento aparte de cada uno fue valioso para reconstruir la secuencia del fragmento objetivo original.
Historia. La primera descripción publicada del uso de extremos emparejados fue en 1990 como parte de la secuenciación del locus HGPRT humano, aunque el uso de extremos emparejados se limitó a cerrar brechas después de la aplicación de un enfoque de secuenciación de escopeta tradicional. La primera descripción teórica de una estrategia de secuenciación de extremos por pares pura, asumiendo fragmentos de longitud constante, fue en 1991. En ese momento, existía un consenso de la comunidad de que la longitud óptima del fragmento para la secuenciación final por pares sería tres veces la longitud de lectura de la secuencia. En 1995, Roach et al. introdujo la innovación de usar fragmentos de diferentes tamaños y demostró que una estrategia pura de secuenciación final por pares sería posible en grandes objetivos. La estrategia fue posteriormente adoptada por el Instituto de Investigación Genómica (TIGR) para secuenciar el genoma de la bacteria Haemophilus influenzae en 1995, y luego por Celera Genomics para secuenciar el genoma de Drosophila melanogaster (mosca de la fruta) en 2000, y posteriormente el genoma humano.
Acercamiento Para aplicar la estrategia, una hebra de ADN de alto peso molecular se corta en fragmentos aleatorios, se selecciona por tamaño (generalmente 2, 10, 50 y 150 kb) y se clona en un vector apropiado. A continuación, los clones se secuencian desde ambos extremos usando el método de terminación de cadena produciendo dos secuencias cortas. Cada secuencia se denomina lectura final o lectura 1 y lectura 2 y dos lecturas del mismo clon se denominan pares de apareamiento. Dado que el método de terminación de la cadena generalmente solo puede producir lecturas de entre 500 y 1000 bases de largo, en todos los clones, excepto en los más pequeños, los pares de apareamiento rara vez se superponen.
Montaje La secuencia original se reconstruye a partir de las lecturas utilizando un software de ensamblaje de secuencias. Primero, las lecturas superpuestas se recopilan en secuencias compuestas más largas conocidas como contigs. Los contigs se pueden unir en andamios siguiendo las conexiones entre pares de aparejos. La distancia entre los contigs se puede inferir de las posiciones de los pares coincidentes si se conoce la longitud promedio del fragmento de la biblioteca y tiene una ventana de desviación estrecha. Dependiendo del tamaño del espacio entre contigs, se pueden usar diferentes técnicas para encontrar la secuencia en los espacios. Si el espacio es pequeño (5-20 kb), entonces se requiere el uso de la reacción en cadena de la polimerasa (PCR) para amplificar la región, seguida de la secuenciación. Si la brecha es grande (>20 kb), el fragmento grande se clona en vectores especiales como los cromosomas artificiales bacterianos (BAC) seguido de la secuenciación del vector.
Pros y contras Los defensores de este enfoque argumentan que es posible secuenciar todo el genoma a la vez utilizando grandes conjuntos de secuenciadores, lo que hace que todo el proceso sea mucho más eficiente que los enfoques más tradicionales. Los detractores argumentan que, aunque la técnica secuencia rápidamente grandes regiones de ADN, su capacidad para unir correctamente estas regiones es sospechosa, particularmente para genomas con regiones repetidas. A medida que los programas de ensamblaje de secuencias se vuelven más sofisticados y la potencia informática se vuelve más barata, puede ser posible superar esta limitación.
Cobertura
La cobertura (profundidad o profundidad de lectura) es el número medio de lecturas que representan un nucleótido dado en la secuencia reconstruida. Se puede calcular a partir de la longitud del genoma original (G), el número de lecturas (N) y la longitud media de lectura (L) como . Por ejemplo, un genoma hipotético con 2.000 pares de bases reconstruidos a partir de 8 lecturas con una longitud promedio de 500 nucleótidos tendrá una redundancia 2x. Este parámetro también permite estimar otras cantidades, como el porcentaje del genoma cubierto por lecturas (a veces también llamado cobertura). Se desea una alta cobertura en la secuenciación de escopeta porque puede superar errores en la llamada y el ensamblaje de bases. El tema de la teoría de la secuenciación del ADN aborda las relaciones de tales cantidades.
A veces se hace una distinción entre cobertura de secuencia y cobertura física. La cobertura de secuencia es el número promedio de veces que se lee una base. La cobertura física es el número medio de veces que se lee o se amplía una base mediante lecturas emparejadas.
Secuencia de escopeta jerárquica miniaturadeimagen| En la secuenciación de escopeta del genoma completo (arriba), todo el genoma se corta al azar en pequeños fragmentos (del tamaño adecuado para la secuenciación) y luego se vuelve a ensamblar. En la secuenciación jerárquica de escopeta (abajo), el genoma se divide primero en segmentos más grandes. Después de deducir el orden de estos segmentos, se cortan en fragmentos de tamaño apropiado para la secuenciación.|350x350px Aunque, en teoría, la secuenciación rápida se puede aplicar a un genoma de cualquier tamaño, su aplicación directa a la secuenciación de genomas grandes (por ejemplo, el genoma humano) fue limitada hasta finales de la década de 1990, cuando los avances tecnológicos hicieron práctico el manejo de grandes cantidades de datos complejos involucrados en el proceso. Históricamente, se creía que la secuenciación de escopeta de genoma completo estaba limitada tanto por el gran tamaño de los genomas grandes como por la complejidad agregada por el alto porcentaje de ADN repetitivo (más del 50% para el genoma humano) presente en genomas grandes. No fue ampliamente aceptado que una secuencia de escopeta de genoma completo de un genoma grande proporcionaría datos confiables. Por estas razones, se tuvieron que utilizar otras estrategias que redujeron la carga computacional del ensamblaje de secuencias antes de realizar la secuenciación de escopeta. En la secuenciación jerárquica, también conocida como secuenciación de arriba hacia abajo, se realiza un mapa físico de baja resolución del genoma antes de la secuenciación real. A partir de este mapa, se selecciona un número mínimo de fragmentos que cubren todo el cromosoma para secuenciar. De esta manera, se requiere la cantidad mínima de secuenciación y ensamblaje de alto rendimiento.
El genoma amplificado primero se corta en trozos más grandes (50-200 kb) y se clona en un huésped bacteriano utilizando BAC o cromosomas artificiales derivados de P1 (PAC). Debido a que se han cortado al azar múltiples copias del genoma, los fragmentos contenidos en estos clones tienen diferentes extremos, y con suficiente cobertura encontrar un andamio de contigs BAC que cubra todo el genoma es teóricamente posible. Este andamio se llama camino de mosaico (tiling path). miniaturadeimagen| Un contig BAC que cubre toda el área genómica de interés constituye el camino del mosaico.|350x350px Una vez que se ha encontrado una ruta de mosaico, los BAC que forman esta ruta se cortan al azar en fragmentos más pequeños y se pueden secuenciar utilizando el método de la escopeta en una escala más pequeña.