
{"id":10,"date":"2013-05-29T20:28:51","date_gmt":"2013-05-29T18:28:51","guid":{"rendered":"https:\/\/bioinfo2.ugr.es\/bioinfo\/?page_id=10"},"modified":"2026-01-25T11:55:13","modified_gmt":"2026-01-25T09:55:13","slug":"analisis-de-secuencias","status":"publish","type":"page","link":"https:\/\/bioinfo2.ugr.es\/bioinfo\/analisis-de-secuencias\/","title":{"rendered":"An\u00e1lisis b\u00e1sico"},"content":{"rendered":"<p>A partir de una secuencia de ADN podemos calcular una serie de propiedades composicionales que nos pueden facilitar informaci\u00f3n valida acerca de la secuencia e incluso del organismo del que proviene. El contenido en G+C de una secuencia es simplemente la fracci\u00f3n de guaninas y citosinas sobre el total de las 4 bases (en ausencia de bases ambiguas, este n\u00famero es id\u00e9ntico a la longitud de la secuencia). Su an\u00e1lisis a lo largo de la secuencia nos puede dar una idea donde se localizan promotores y exones, regiones que se suelen caracterizar por destacar en su contenido en G+C sobre su entorno. Otras aplicaciones son la detecci\u00f3n de isocoras, islas CpG o eventos de transferencia horizontal entre procariotas.<\/p>\n<p>El\u00a0paquete de software EMBOSS dispone de diferentes programas para analizar propiedades b\u00e1sicas de secuencias de ADN y prote\u00ednas<\/p>\n<ul>\n<li><a href=\"http:\/\/emboss.sourceforge.net\/index.html\" target=\"_blank\" rel=\"noopener noreferrer\">Web site\u00a0<\/a> (<a href=\"https:\/\/emboss.sourceforge.net\/apps\/cvs\/emboss\/apps\/index.html\" target=\"_blank\" rel=\"noopener\">lista de programas<\/a>)<\/li>\n<li><a href=\"https:\/\/emboss.sourceforge.net\/docs\/emboss_tutorial\/emboss_tutorial.html\" target=\"_blank\" rel=\"noopener\">Un tutorial<\/a><\/li>\n<li>P<a href=\"https:\/\/www.ebi.ac.uk\/jdispatcher\/\" target=\"_blank\" rel=\"noopener\">rogramas disponibles en el EBI<\/a><\/li>\n<li><a href=\"https:\/\/emboss.toulouse.inra.fr\/cgi-bin\/emboss\/\" target=\"_blank\" rel=\"noopener\">servidor toulouse<\/a><\/li>\n<\/ul>\n<div>\n<hr \/>\n<h3>Secuencias an\u00f3nimas:<\/h3>\n<ul>\n<li><a href=\"https:\/\/bioinfo2.ugr.es\/bioinfo\/wp-content\/uploads\/2024\/01\/anon2.txt\" target=\"_blank\" rel=\"noopener\">Secuencia del genoma humano\u00a0<\/a><\/li>\n<li><a href=\"https:\/\/bioinfo2.ugr.es\/bioinfo\/wp-content\/uploads\/2024\/01\/anon3.txt\" target=\"_blank\" rel=\"noopener\">Genoma de un procariota<\/a><\/li>\n<\/ul>\n<table style=\"height: 174px; width: 100%; border-collapse: collapse; background-color: #e6dfdf;\">\n<tbody>\n<tr style=\"height: 174px;\">\n<td style=\"width: 100%; height: 174px;\"><strong>Necesitaremos el sistema operativo Linux:<\/strong><\/p>\n<ul>\n<li><a href=\"https:\/\/bioinfo2.ugr.es\/bioinfo\/linux\/\" target=\"_blank\" rel=\"noopener\">Informaci\u00f3n general<\/a><\/li>\n<li><a href=\"https:\/\/bioinfo2.ugr.es\/bioinfo\/ejercicios-linux\/\" target=\"_blank\" rel=\"noopener\">Ejercicios Linux<\/a><\/li>\n<\/ul>\n<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<h3><span id=\"Programas_de_EMBOSS\">Analisis b\u00e1sico mediante EMBOSS en linea de comando<\/span><\/h3>\n<h3><a href=\"http:\/\/emboss.bioinformatics.nl\/cgi-bin\/emboss\/infoseq\" target=\"_blank\" rel=\"noopener noreferrer\">infoseq<\/a><\/h3>\n<p>Este programa nos permite extraer anotaciones de una secuencia y determinar propiedades b\u00e1sicas como la longitud y el contenido en G+C<\/p>\n<p>Ejecutar el programa visualizando todos los parametros<\/p>\n<pre>infoseq --help<\/pre>\n<p>Ejecutar infoseq proporcionando el fichero que queremos analizar:<\/p>\n<pre>infoseq -sequence anon2.txt -outfile anon2_infoseq.txt<\/pre>\n<h3><a href=\"http:\/\/emboss.bioinformatics.nl\/cgi-bin\/emboss\/extractseq\" target=\"_blank\" rel=\"noopener noreferrer\">extractseq<\/a><\/h3>\n<p>Extraer una subsecuencia.<\/p>\n<ol>\n<li>Utilizando <i>extractseq<\/i>, haga el splicing del <a href=\"http:\/\/www.ncbi.nlm.nih.gov\/nuccore\/93007319?report=genbank\" target=\"_blank\" rel=\"noopener\">gen de la rodopsina<\/a> de <i>Xenopus<\/i> y extraiga la CDS completa.<\/li>\n<li>Use <strong><em>plotorf<\/em><\/strong> o <a href=\"https:\/\/www.ncbi.nlm.nih.gov\/orffinder\/\" target=\"_blank\" rel=\"noopener\">ORFfinder<\/a> para visualizar las ORFs en la CDS obtenida y confirmar que el splicing ha sido correcto.<\/li>\n<li>Tambi\u00e9n puede utilizar <em><strong>cusp<\/strong> <\/em>(v\u00e9ase m\u00e1s adelante)\u00a0para contar los codones y confirmar que esta CDS tiene un solo codon de stop al final de la secuencia.<\/li>\n<\/ol>\n<h3><a href=\"http:\/\/emboss.bioinformatics.nl\/cgi-bin\/emboss\/revseq\" target=\"_blank\" rel=\"noopener noreferrer\">revseq<\/a><\/h3>\n<p>Programa que permite determinar la secuencia complementaria inversa<\/p>\n<hr \/>\n<h2><span id=\"Analisis_composicional\">An\u00e1lisis composicional<\/span><\/h2>\n<div class=\"entry-content clearfix\">\n<h3>Fluctuaci\u00f3n espacial de la composici\u00f3n: <a href=\"http:\/\/emboss.bioinformatics.nl\/cgi-bin\/emboss\/freak\" target=\"_blank\" rel=\"noopener noreferrer\"><b><i>freak<\/i><\/b><\/a><\/h3>\n<p>Este programa nos permite determinar la variaci\u00f3n espacial del %GC a lo largo de una secuencia<\/p>\n<ul>\n<li>Determinar la fluctuaci\u00f3n del contenido en G+C a lo largo de esta secuencia: <a href=\"http:\/\/bioinfo2.ugr.es\/secuencias\/wp-content\/uploads\/2015\/11\/MHC_LH_clean.txt\" target=\"_blank\" rel=\"noopener\">secuencia humana<\/a><\/li>\n<li>Determinar la fluctuaci\u00f3n a lo largo de las secuencias anon2 y anon3.\n<ul>\n<li>\u00bfQu\u00e9 diferencias podemos observar entre las secuencias del genoma humano y de la Bacteria?<\/li>\n<li>\u00bfQue podr\u00edan representan las estructuras detectadas?<\/li>\n<li>\u00bfComo se puede verificar la identidad de estas estructuras?<\/li>\n<\/ul>\n<\/li>\n<\/ul>\n<p><strong>Ejemplo de uso:<\/strong><\/p>\n<\/div>\n<\/div>\n<div>\n<div class=\"entry-content clearfix\">\n<pre>freak -seqall anon3.txt -letters gc -outfile anon.freak -step 500 -window 1000\n<\/pre>\n<p><strong>Con salida gr\u00e1fica:<\/strong><\/p>\n<pre>freak -seqall anon3.txt -letters gc -outfile anon.freak -step 500 -window 1000 -graph png -plot\n\n<\/pre>\n<\/div>\n<div class=\"entry-content clearfix\">\n<h3>An\u00e1lisis de las frecuencias de k-meros: <a href=\"http:\/\/emboss.bioinformatics.nl\/cgi-bin\/emboss\/compseq\" target=\"_blank\" rel=\"noopener noreferrer\"><b><i>compseq<\/i><\/b><\/a><\/h3>\n<p>La composici\u00f3n de una secuencia de ADN se puede estudiar determinando las frecuencias de mononucle\u00f3tidos, dinucle\u00f3tidos, etc<\/p>\n<p>Para ello se usa una ventana movil que se va desplazando posici\u00f3n a posici\u00f3n a lo largo de la secuencia (Por ejemplo, para determinar la frecuencia de dinucle\u00f3tidos se usa una ventana de tama\u00f1o 2 y salto de 1:<\/p>\n<p>ATTCCGTGAACTG\u2026<\/p>\n<p>AT, TT, TC, CG, GT\u2026<\/p>\n<ul>\n<li>Determinar las frecuencias de dinucle\u00f3tidos de esta <a href=\"https:\/\/bioinfo2.ugr.es\/bioinfo\/wp-content\/uploads\/2022\/01\/secuencia_genoma.txt\">secuencia an\u00f3nima (hsa)<\/a> y de la bacteria <a href=\"http:\/\/www.ncbi.nlm.nih.gov\/nuccore\/NC_006461.1\" target=\"_blank\" rel=\"noopener\">Thermus thermophilus<\/a><\/li>\n<li>Determinar las frecuencias de dinucle\u00f3tidos de una secuencia de C. elegans mediante el <a href=\"http:\/\/genome.ucsc.edu\/index.html\" target=\"_blank\" rel=\"noopener\">navegador gen\u00f3mico de Santa Cruz<\/a><\/li>\n<li>Comparar los resultados con los de una secuencia aleatorizada mediante <a href=\"http:\/\/emboss.bioinformatics.nl\/cgi-bin\/emboss\/shuffleseq\" target=\"_blank\" rel=\"noopener noreferrer\">suffleseq<\/a>.<\/li>\n<\/ul>\n<p><em><strong>Cuestiones<\/strong><\/em><\/p>\n<ul>\n<li>\u00bfQu\u00e9 diferencias podemos observar entre la secuencia del genoma humano de de la Bacteria?\u00a0\u00bfA que se pueden deber estas diferencias?<\/li>\n<li>\u00bfQue observamos en <em>C. elegans<\/em>?\u00bfQu\u00e9 explicaci\u00f3n puede tener?<\/li>\n<\/ul>\n<hr \/>\n<h2><span id=\"Uso_de_codones\">Uso de codones<\/span><\/h2>\n<p>Es importante subrayar que los codones no se cuentan c\u00f3mo los trinucle\u00f3tidos. Estos \u00faltimos se cuentan mediante <i><b>compseq <\/b><\/i>con una ventana solapante. Los codones son tambi\u00e9n trinucl\u00e9otidos, pero NO solapantes. Otra cuesti\u00f3n a tener en cuenta es que los codones hay que contarlos en la fase correcta. El recuento de codones hay que hacerlo pues con el programa <i><b>cusp<\/b><\/i>. El uso de codones (o dialecto gen\u00e9tico) es especie-espec\u00edfico en organismos unicelulares y regi\u00f3n-espec\u00edfico en multicelulares.<\/p>\n<ul>\n<li><a href=\"http:\/\/www.ncbi.nlm.nih.gov\/Taxonomy\/Utils\/wprintgc.cgi?mode=c#SG1\" target=\"_blank\" rel=\"noopener noreferrer\">Tabla del c\u00f3digo gen\u00e9tico standard<\/a><\/li>\n<li><a href=\"http:\/\/www.ncbi.nlm.nih.gov\/nuccore\/NM_000558.3\" target=\"_blank\" rel=\"noopener noreferrer\">CDS de la alfa-globina humana<\/a><\/li>\n<li><a href=\"http:\/\/emboss.bioinformatics.nl\/cgi-bin\/emboss\/cusp\" target=\"_blank\" rel=\"noopener noreferrer\">cusp \u2013 EMBOSS<\/a><\/li>\n<li><a href=\"http:\/\/www.kazusa.or.jp\/codon\/\" target=\"_blank\" rel=\"noopener noreferrer\">Codon Usage Database<\/a><\/li>\n<li>Codon optimization: <a href=\"https:\/\/www.genscript.com\/codon-opt.html?src=google&amp;gclid=Cj0KCQiA7OnxBRCNARIsAIW53B81InatqrqI_WsMt4y6dvveZSXm2KoXGIM2DUYbAemQZhX--djYQHoaAr2kEALw_wcB\" target=\"_blank\" rel=\"noopener noreferrer\">GenScript<\/a> <a href=\"https:\/\/www.genewiz.com\/en-GB\/Public\/Services\/Gene-Synthesis\/Codon-Optimization\/\" target=\"_blank\" rel=\"noopener noreferrer\">GeneWiz<\/a><\/li>\n<\/ul>\n<p><strong>Ejercicios<\/strong><\/p>\n<ul>\n<li>Analizar el uso de codones de la bacteria <a href=\"http:\/\/www.ncbi.nlm.nih.gov\/nuccore\/NC_006461.1\" target=\"_blank\" rel=\"noopener\">Thermus thermophilus<\/a>\n<ol>\n<li>Obtener todas las secuencias de las CDS<\/li>\n<li>Obtener el uso de codones mediante el programa <a href=\"http:\/\/emboss.bioinformatics.nl\/cgi-bin\/emboss\/cusp\" target=\"_blank\" rel=\"noopener\">CUSP<\/a><\/li>\n<\/ol>\n<\/li>\n<li>Analizar el uso de codones en <a href=\"http:\/\/www.ncbi.nlm.nih.gov\/nuccore\/NC_007633.1\" target=\"_blank\" rel=\"noopener\">Mycoplasma capricolum<\/a><\/li>\n<li>Usar la base de datos del NCBI: <a href=\"http:\/\/www.ncbi.nlm.nih.gov\/genome\/browse\/\">http:\/\/www.ncbi.nlm.nih.gov\/genome\/browse\/<\/a>\n<ol>\n<li>Seleccionar una bacteria<\/li>\n<li>Determinar el uso de codones y el G+C del genoma<\/li>\n<\/ol>\n<\/li>\n<\/ul>\n<p><strong>Cuestiones<\/strong><\/p>\n<ul>\n<li>\u00bfA que se se pueden deber las diferencias en el contenido en G+C?<\/li>\n<li>\u00bfPor que el G+C es tan diferente en las distintas posiciones del cod\u00f3n?<\/li>\n<\/ul>\n<\/div>\n<div class=\"entry-content clearfix\">\n<h2><del><span id=\"Patrones_en_el_ADN\">Patrones en el ADN<\/span><\/del><\/h2>\n<div class=\"entry-content clearfix\">\n<p><del>Existen varias utilidades de EMBOSS para el descubrimiento de patrones (o motivos), tanto en el ADN (<a href=\"http:\/\/www.bioinformatics.nl\/cgi-bin\/emboss\/fuzznuc\" target=\"_blank\" rel=\"noopener\"><b>fuzznuc<\/b><\/a>) como en las prote\u00ednas (<b><a href=\"http:\/\/www.bioinformatics.nl\/cgi-bin\/emboss\/fuzzpro\" target=\"_blank\" rel=\"noopener\">fuzzpro<\/a>, antigenic<\/b>).<\/del><\/p>\n<p><del>Una manera visualizar estos patrones es mediante el llamado \u2018juego del caos\u2019. El programa de EMBOSS que permite obtener una representaci\u00f3n ca\u00f3tica de una secuencia de ADN es <a href=\"http:\/\/emboss.bioinformatics.nl\/cgi-bin\/emboss\/chaos\" target=\"_blank\" rel=\"noopener noreferrer\">chaos<\/a>. Puede usar la secuencia de la regi\u00f3n de la <a href=\"http:\/\/www.ncbi.nlm.nih.gov\/nuccore\/455025?report=fasta\" target=\"_blank\" rel=\"noopener noreferrer\">beta-globina humana<\/a> para obtener su representaci\u00f3n ca\u00f3tica, y compararla con la imagen que se obtiene tras aleatorizar la secuencia mediante <a href=\"http:\/\/emboss.bioinformatics.nl\/cgi-bin\/emboss\/shuffleseq\" target=\"_blank\" rel=\"noopener noreferrer\">suffleseq<\/a>.<\/del><\/p>\n<hr \/>\n<h2><del><span id=\"Mapas_de_restriccion\">Mapas de restricci\u00f3n<\/span><\/del><\/h2>\n<p><del>El programa de EMBOSS <em><strong>restrict\u00a0<\/strong><\/em>permite obtener un listado de las dianas que presenta una secuencia problema de ADN para distintas enzimas de restricci\u00f3n. Se pueden especificar una o dos enzimas de restricci\u00f3n, o bien todas las enzimas de REBASE. Existen opciones para controlar los sitios que se quieren obtener y para elegir el formato de salida.<\/del><\/p>\n<\/div>\n<\/div>\n<\/div>\n<p><!--more--><\/p>\n<h2><a href=\"https:\/\/bioinfo2.ugr.es\/bioinfo\/resultados_analisis_basico\/\" target=\"_blank\" rel=\"noopener\">Resultados<\/a><\/h2>\n","protected":false},"excerpt":{"rendered":"<p>A partir de una secuencia de ADN podemos calcular una serie de propiedades composicionales que nos pueden facilitar informaci\u00f3n valida acerca de la secuencia e incluso del organismo del que proviene. El contenido en G+C de una secuencia es simplemente la fracci\u00f3n de guaninas y citosinas sobre el total de las 4 bases (en ausencia &hellip; <a href=\"https:\/\/bioinfo2.ugr.es\/bioinfo\/analisis-de-secuencias\/\">Continue reading <span class=\"meta-nav\">&rarr;<\/span><\/a><\/p>\n","protected":false},"author":1,"featured_media":0,"parent":0,"menu_order":5,"comment_status":"closed","ping_status":"closed","template":"","meta":[],"_links":{"self":[{"href":"https:\/\/bioinfo2.ugr.es\/bioinfo\/wp-json\/wp\/v2\/pages\/10"}],"collection":[{"href":"https:\/\/bioinfo2.ugr.es\/bioinfo\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/bioinfo2.ugr.es\/bioinfo\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/bioinfo2.ugr.es\/bioinfo\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/bioinfo2.ugr.es\/bioinfo\/wp-json\/wp\/v2\/comments?post=10"}],"version-history":[{"count":36,"href":"https:\/\/bioinfo2.ugr.es\/bioinfo\/wp-json\/wp\/v2\/pages\/10\/revisions"}],"predecessor-version":[{"id":1966,"href":"https:\/\/bioinfo2.ugr.es\/bioinfo\/wp-json\/wp\/v2\/pages\/10\/revisions\/1966"}],"wp:attachment":[{"href":"https:\/\/bioinfo2.ugr.es\/bioinfo\/wp-json\/wp\/v2\/media?parent=10"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}