9
REVISTA CON-CIENCIA N°1/VOL.6: 107-124. JUNIO 2018. ISSN: 2310-0265
Diseño y evaluación de primers in silico del gen E1 del virus de chikungunya para Real-Time PCR (qPCR)
Design and Evaluation of primers in silico of the E1 gene of chukungunya virus for Real -Time PCR (q PCR)
DE LA FUENTE, ARIEL1
ROMERO CALLE, DANITZA XIOMARA1 CÁRDENAS ALEGRÍA, OSCAR1 ÁLVAREZ ALIAGA, MARÍA TERESA1
FECHA DE RECEPCIÓN: 13 DE MARZO DE 2018 FECHA DE ACEPTACIÓN: 9 DE MAYO DE 2018
Abstract
Primer design is a substantial part of PCR assays, which is used in characteriza- tion of microorganisms, sequencing, quanti- fication, etc. In recent years, new tools have been added for this purpose, new variants and approaches have been created, many of these tools are commercial and others are free including open source, giving the pos- sibility to improve or modify as required. In this work, we compared the basic charac- teristics of the primers design, using 2 types of free and open source web tools, Primer- BLAST and Primer3Plus, in this sense the in silico design, analysis of secondary struc- tures and verification of primers designed. The results showed that the use of addi- tional variables or use of penalties using the software Primer3Plus affect in specificity, secondary structures formation in primers and in the products of amplification. In ad- dition, default parameters, that many soft-
Instituto de Investigaciones Fármaco Bioquímicas (IIFB), Facultad de Ciencias Farmacéuticas y Bioquímicas, UMSA.
ware have, do not imply that suitable and / or specific primers can be obtained, mean- ing the importance of the use of penalties. Therefore, the use of Primer3Plus has al- lowed to decrease the formation of second- ary structures without affecting the ampli- fication specificity of the selected marker.
Diseño de primers, estructuras secundarias, qPCR y software.
KEY WORDS
Design of primers, secondary structures, qPCR and software.
La tecnología genómica moderna permite estudiar millones de regiones de cada muestra de ADN. La mayoría de estas tecnologías requiere la amplifi- cación de la reacción en cadena de la polimerasa (PCR), técnica que requiere el uso de primers específicos los mismos no deben formar estructuras secun- darias que afecten el rendimiento de la PCR.
El amplio uso de primers universales o publicados en artículos de inves- tigación de referencia son generalmente el resultado del diseño primers, el cual requiere experiencia y destreza por estar sujeto a posibles errores (Ca- milo et al., 2016), en ausencia de primers de referencia es imprescindible el diseño de los mismos. En este sentido, existe una demanda reciente de méto- dos o algoritmos para su diseño, debido también a la proliferación de nuevas tecnologías dentro la genómica y otros campos relacionados (Kim, Kang, An, Koo, & Kim, 2016), además se ha reportado que estos métodos presentan re- sultados favorables (Chuang, Cheng, & Yang, 2013) en los diferentes ensayos e investigaciones incluyendo aplicaciones clínicas (Thierry, y otros, 2014) para la identificación de diferentes microorganismos emergentes. No obstante, va- rios de estos primers han sido diseñados para marcadores específicos o par- ticulares, por lo tanto, la amplificación de estos marcadores no podría ser las mismas o semejantes, afectando de esta forma el éxito del PCR.
Aunque no todos los aspectos del proceso de la qPCR están completa- mente descritos y modelados, se conocen ciertos factores que afectan la tasa de éxito de la qPCR, en sistemas que usan el SYBER1 como intercalan- te Green, la detección es más simple y económica de qPCR pero la unión del primer a regiones repetidas de la secuencia de DNA puede ser errónea y ge- nerar productos inespecíficos, otro factor determinante es la temperatura de fusión (Tm) (Pachón, 2017). Los falsos positivos resultan en la formación de primers dimers y amplificación no especifica. El éxito de la PCR depende de la calidad del templado, Taq polimerasa, solución buffer y el diseño de pri- mers en los cuales debe haber un balance entre especificidad y eficiencia.
El desarrollo informático ha permitido la creación de diversos softwares que facilitan el diseño de primers incrementando la tasa de éxito del qPCR. Los cuales permiten tanto diseño automático como diseño particular o espe- cifico, de tal forma que nuevos usuarios puedan realizar el diseño manual, y un usuario con experiencia puede realizar un diseño más avanzado (Wang et al., 2016; Lee et al., 2015; Yang et al., 2015), en el que se emplea el uso de pa- rámetros adicionales que son un conjunto de variables que intervienen en su características (Presti et al., 2014).
Existen una variedad de software para el diseño de primers muchos de ac- ceso libre como PrimerBlast y Primer3Plus empleado. Los autores del diseño y la creación digital de Primer3Plus, han examinado las versiones más recien- tes de navegadores tales como Firefox, Internet Explorer, Safari and Konque- ror, permitiendo su uso en cualquier navegador. Por otro lado, Primer3Plus funciona con JavaScript, Microsoft Windows y en cualquier sistema operati- vo basado en Unix/Linux.
Cuando se diseñan y analizan primers se deben tomar en cuenta princi- palmente 3 aspectos descritos por Thornton & Basu: 1) Los dímeros y horqui- llas de los primers que se produzcan no deberán tener un ΔG muy negativo (< -3.5 kcal/mol), 2) No debe existir horquillas en el extremo 3’, y 3) Las es- tructuras secundarias que resulten, no deben tener un Tm mayor a la tempe- ratura de alineamiento.
El diseño de primers está siendo extensamente implementado en mu- chos laboratorios, sean clínicos o de investigación, permitiendo reducir cos- tos, también permite tener un control personalizado y más completo sobre el experimento de la PCR.
El presente estudio tiene por finalidad diseñar primers para el gen E1 del virus Chikungunya mediante 2 tipos de software, ambos de acceso libre, Pri- merBlast y Primer3Plus para qPCR. La identificación de este gen (E1 codifica una proteína de membrana conocida por incrementar su habilidad de infec- ción posterior a su mutación) podría ser útil para emplearlo como marcador genético si inicialmente se realiza una retrotranscripción puesto que el virus de tipo RNA de cadena sencilla.
El diseño de primers inicialmente se seleccionó el marcador o región ge- nética del blanco, en gen E1 del virus de Chikunguya (Rodríguez et al., 2015), se utilizó dos programas diferentes para diseño de oligonucleótidos denomi- nados PrimerBlast y Primer3Plus y posteriormente se realizó el análisis de los oligos diseñados con el programa Beacon Designer Free Edition y del ampli- cón resultante mediante UNAFold, figura1.
La secuencia nucleotídica (563 bp) del gen E1 del virus Chikungunya (CHKV) (Evia, 2015) fue obtenida mediante una búsqueda en la base de datos GenBank, número de acceso JX142137.1. (Benson et al., 2000) para el diseño de los primers orientados al ensayo de PCR cuantitativo (qPCR) mediante pa- rámetros estándar descritos por Rodríguez et al. (2015).
Se ha diseñado primers que amplifiquen la región del gen E1 del virus CHKV mediante dos sitios web que proporcionan software online: 1) Primer- BLAST (Ye et al., 2012) y 2) Primer3Plus (Untergasser, et al., 2007).
El diseño de primers en el software Primer-BLAST del NCBI (https://www.ncbi.nlm.nih.gov/tools/primer-blast/index.cgi?LINK_ LOC=BlastHome) fue realizado de acuerdo a Rodríguez et al. (2015) y Untergasser et al. (2015), para el resto de los parámetros se utilizó los valores predeterminados por el software, figura2 y tabla 1 (Rodriguez et al., 2015; Untergasser, et al., 2007).
PARÁMETRO | VALOR | DESCRIPCIÓN |
Secuencia | JX142137.1 | Nombre de la secuencia |
Sequence ID | Chikungunya virus CHIKV/IRL/2009 E1 gene | Código de la secuencia |
Product Size Ranges | 80-150 100-200 | Rangos de tamaño para los amplicones |
Primer Size | 20 25 28 | Longitudes mínimas, optimas y máximas para los primers diseñados |
Primer Tm | 60 64 70 | Tm mínimo, optimo y máximo para los primers diseñados |
Max Tm Difference | 2 | Diferencia máxima aceptable entre los Tm de los primers |
Product Tm | 50 | Tm mínimo, optimo y máximo para el amplicón |
Primer GC% | 35 65 80 | Porcentajes mínimos, óptimos y máximos de G y C en cualquier primer |
GC Clamp | 0 | Numero especifico de Gs y Cs consecutivos en el extremo 3’ de ambos primers |
Concentration of divalent cations | 3.5 | Concentración milimolar de cationes divalentes |
Concentration of dNTPs | 0.20 | Concentración milimolar de trifosfatos deoxiribonucleotidos |
Penaltis para los primers (for Primers) | ||
Tm | Lt=0 - 4 Gt=0 - 4 | Describen los penaltis que permiten modificar el criterio que usa Primer3 para seleccionar los mejores primers |
Size | Lt=0 y1 Gt=0 y1 | |
Self-Complementarity | 0 y 3 | |
3' Self Complementarity | 0 y 3 | |
Penaltis para el par de primers (for Primer Pairs) | ||
Product tm | Lt=0 y1 Gt=0 y1 | Describen los penaltis que permiten modificar y ajustar el criterio que usa Primer3 para seleccionar los mejores primers |
Tm Difference | 0 y 2 | |
Any Complementarity | 0 y 3 | |
3' Self Complementarity | 0 y 3 |
Primer3Plus (Bioinformatics, http://www.bioinformatics.nl/cgi- bin/primer3plus/primer3plus.cgi) fue empleado para el diseño de primers(figura 3), inicialmente la secuencia del gen E1 en formato FAS- TA fue ingresada en Principal (Main), posteriormente se asignaron va- lores fisicoquímicos óptimos para el diseño de primers (tabla 1) en la sección de configuraciones avanzadas (Advanced Settings), consecuti- vamente los penaltis (Penalty weights) fueron ajustados en un rango de 0 a 4 en mediante la variación de la temperatura de fusión (Tm) uno de los parámetros más influyentes para la obtención de primers específi- cos (Untergasser et al., 2007; Rodríguez et al., 2015), de acuerdo a estas modificaciones se obtuvieron 6 grupos de primers (A, B, C, D, E, y F), en el grupo A se asignó el valor de 0 a todos los parámetros en la sección de Penaltis y para el resto de los grupos se asignó los valores de penal- tis predeterminados por el software, el resto de los parámetros fueron estándar (Thornton & Basu, 2011).
La calidad y especificidad de los primers obtenidos fueron analizados em- pleando Beacon Designer Free Edition (PRIMER Biosoft, http://www.pre- mierbiosoft.com/qOligo/Oligo.jsp?PID=1), se ingresó el par de primers den- tro su interface, seleccionando la opción SYBR®Green y dejando el resto de los campos de forma predeterminada (Figura 4). Asimismo, se identificó la formación de posibles estructuras secundarias (Apte, 2007).
Por otro lado, se analizaron los amplicones resultantes mediante el soft- ware UNAFold (IDT, Integrated DNA Technologies, Inc., https://www.idtdna. com/UNAFold), se ingresó la secuencia del amplicon resultante y asigno el valor de 60 ºC a la temperatura y realizó click en Submit, consecutivamen- te las variables termodinámicas de los amplicones y la formación de dímeros fueron evaluados.
Finalmente, se determinó la especificidad de los primers de acuerdo a los valores de identidad (Ident) y E-Value presentados en Nucleotide-BLAST (Na- tional Center for Biotechnology Information, U.S. National Library of Medi- cine 8600 Rockville Pike, Bethesda MD, 20894 USA https://blast.ncbi.nlm.nih. gov/Blast.cgi?PAGE_TYPE=BlastSearch).
Inicialmente, se utilizó la secuencia obtenida con número de acce- so JX142137.1 (Figura 5), se ingresó esta secuencia en la interfaz de Primer- BLAST (Figura 2).
El tiempo de procesamiento del programa para la búsqueda de los pri- mers fue aproximadamente 72 segundos, al terminar el programa reporto una lista de los primers encontrados (Figura 6), en esta lista el primer par de pri- mers obtenidos se muestran en la tabla 2.
Primer | Secuencia |
Forward | AATTCCGCGTACAAGTCTGTTCTAC |
Reverse | GTAGTTGACTATGTGGTCCTTCGGG |
Posteriormente, el par de primers obtenidos fueron analizados de forma
In silico con los programas correspondientes.
Inicialmente, se analizó el par de primers con Beacon Designer Free Edi- tion, diseñado para la evaluación de oligos por sus propiedades termodiná- micas y estructuras secundarias. Posteriormente, se analizó la secuencia am- plificada con UNAFold, que permite predecir los plegamientos o hibridación de ácidos nucleicos (Markham & Zuker, 2008).
Los valores fisicoquímicos, los dímeros y horquillas que podrían formar los primers en el experimento de PCR fueron determinados mediante el soft- ware Beacon Designer (Figura 7).
De acuerdo a las condiciones de reacción estándar, los resultados obteni- dos presentan un vínculo de color rojo que al activarlos se observa la sección donde se encuentran las estructuras secundarias identificadas, agrupadas se- gún el tipo de estructura secundaria, Cross Dimer son las formaciones en- tre Forward y Reverse, Self Dimer son las formaciones entre dos primers del mismo tipo. Para este experimento en la Figura 7 se observan los resultados con ΔG negativos, lo que implica la formación de varios dímeros y horquillas.
Para el primer forward (Sense primer) el valor de ΔG fue de -5,2 kcal/ mol en categorizado un valor bajo de acuerdo a los valores de ΔG re- portados por Thornton &Basu, 2011, los mismos no deben ser meno- res a -3,5 kcal/mol. Por lo tanto, existe la probabilidad que se formen estructuras secundarias de tipo ‘Self Dimer’, la energía necesaria para romper los enlaces de estos dímeros tendrá que ser mayor. Sin embar- go, se debe observar si los dímeros tienen interacciones en los extre- mos 3'.
Para el primer reverse (Anti-sense primer) la posible formación del ‘Self Dimer’ presento un ΔG de -1,5 kcal/mol.’.
Asimismo, se observó la presencia de una horquilla (ΔG bajo -1 kcal/mol) y un dímero en el primer Forward con interacción en el extremo 3’ (Lorenz, 2012). Asimismo, en el parámetro ‘Cross Dimer’ se observó valores inferiores a -1 kcal/mol lo cual indica la formación dímeros con interacciones en el ex- tremo 3’ con (Thornton & Basu, 2011).
El amplicon obtenido de 80 bp fue analizado mediante el software UNA- Fold, se identificó una estructura secundaria (Figura 4). Asimismo, se determi-
nó que la temperatura fusión era superior a la temperatura de alineamiento, por tanto, la formación de estructuras secundarias podría interferir en el de- sarrollo de la qPCR.
Por lo tanto, el Grupo de primers diseñados mediante Primer-BLAST para amplificar el gen E1 del virus CHKV deben ser descartados, y se deben eva- luar otros primers.
Se obtuvo 6 grupos de primers: A, B, C, D, E y F de acuerdo a la variación de penaltis, los grupos C y D no formaron estructuras secundarias (ΔG nega- tivas) de acuerdo a el software Beacon Designer, el primer reverse del Grupo B (Tabla 3), el ΔG se encuentra dentro de los rangos de referencia (Thornton & Basu, 2011). Adicionalmente, se observa que los otros grupos si bien presen- tan una temperatura dentro del rango establecido, forman estructuras secun- darias. La temperatura de fusión no debe superar los 60 ºC (Tm) (Thornton & Basu, 2011), ya que podría existir la formación de estructuras secundarias.
118
De la Fuente, A. y col.
GRUPO | PENALTI (Tm) | PRIMERS | Longitud | TM (°C) | GC% | GC CLAMP | CROSS DIMER | SELF DIMER | HORQUI- LLAS | ESTRUCTURAS SECUNDARIAS | CROSS DIMER | ||
(∆G) | (∆G) | (∆G) | Dímero | Horquilla | |||||||||
A | 0 | F | GTTCAGCAGGCACTAACTTGAC | 22 | 57.28 | 50 | 1 | -2.4 | -0.9 | -0.9 | 2 | 2 | 6 |
R | GGTTGCGTAGCCCTTTGATCT | 21 | 57.89 | 52.38 | 1 | -2.4 | -2.0 | 0.0 | 1 | 0 | |||
B | 0 | F | GTGTGGGACTGGTTGTTGCTGTTG | 24 | 62.62 | 54.17 | 1 | -1.1 | 0.0 | 0.0 | 0 | 0 | 9 |
R | TGTCAAGTTAGTGCCTGCTGAACGA | 25 | 62.68 | 48 | 1 | -1.1 | -0.9 | -0.9 | 2 | 2 | |||
C | 1 | F | CGGCGTCACATACCACCCTC | 20 | 59.9 | 65 | 1 | 0.0 | 0.0 | 0.0 | 0 | 0 | 0 |
R | AACAGCAACAACCAGTCCCACA | 22 | 60.26 | 50 | 1 | 0.0 | 0.0 | 0.0 | 0 | 0 | |||
D | 2 | F | CGGCGTCACATACCACCCTC | 20 | 59.9 | 65 | 1 | 0.0 | 0.0 | 0.0 | 0 | 0 | 0 |
R | AACAGCAACAACCAGTCCCACA | 22 | 60.26 | 50 | 1 | 0.0 | 0.0 | 0.0 | 0 | 0 | |||
E | 3 | F | CTACCCGGCGTCACATACCA | 20 | 58.63 | 60 | 2 | 0.0 | -4.3 | 0.0 | 1 | 0 | 0 |
R | CAACAGCAACAACCAGTCCCA | 21 | 58.5 | 52.38 | 2 | 0.0 | 0.0 | 0.0 | 0 | 0 | |||
F | 4 | F | CTACCCGGCGTCACATACCA | 20 | 58.63 | 60 | 2 | 0.0 | -4.3 | 0.0 | 1 | 0 | 0 |
R | CAACAGCAACAACCAGTCCCA | 21 | 58.5 | 52.38 | 2 | 0.0 | 0.0 | 0.0 | 0 | 0 |
Se observaron estructuras secundarias en los grupos de primers A, B, E y F, no obstante, los grupos A y B están dentro de los rangos aceptables (> -3.5 kcal/mol y las interacciones entre las bases están distantes al extremo 3’). Sin embargo, los grupos E y F deben descartarse porque presentan ΔG muy ne- gativas y muestran interacciones en más de 2 pares de bases, estas interaccio- nes pondrían afectar el rendimiento de la PCR.
GRUPO A | |||
Dímeros Secuencia Forward: | ∆G | ||
5' GTTCAGCAGGCACTAACTTGAC 3' | -0.9 | 5'GTTCAGCAGGCACTAACTTGAC 3' | -0.7 |
||| ¦¦¦ | ||| ¦¦¦ | ||
3' CAGTTCAATCACGGACGACTTG 5' | 3' CAGTTCAATCACGGACGACTTG 5' | ||
Horquillas Secuencia Forward: | |||
/CGGACGACTTG 5' | -0.9 | /ACGACTTG 5' | -0.7 |
A ||| | G ||| | ||
\CTAACTTGAC 3' | \GCACTAACTTGAC 3' | ||
Dímeros Secuencia Reverse: | |||
5' GGTTGCGTAGCCCTTTGATCT 3' | -2.0 | ||
|||| | |||
3' TCTAGTTTCCCGATGCGTTGG 5' | |||
Dímeros entre Primers Forward y Reverse: | |||
5' GTTCAGCAGGCACTAACTTGAC 3' | -2.4 | 5' GTTCAGCAGGCACTAACTTGAC 3' | -2.0 |
¦ ¦¦¦ ||| | ¦ ¦ ||| ¦ | ||
3' TCTAGTTTCCCGATGCGTTGG 5' | 3' TCTAGTTTCCCGATGCGTTGG 5' | ||
5' GTTCAGCAGGCACTAACTTGAC 3' | -2.0 | 5' GTTCAGCAGGCACTAACTTGAC 3' | -1.3 |
¦ ||| | ||| ¦ | ||
3' TCTAGTTTCCCGATGCGTTGG 5' | 3' TCTAGTTTCCCGATGCGTTGG 5' | ||
5' GTTCAGCAGGCACTAACTTGAC 3' | -0.7 | 5' GTTCAGCAGGCACTAACTTGAC 3' | -0.1 |
¦¦ ||| | ¦ ¦ ||| ¦ | ||
3' TCTAGTTTCCCGATGCGTTGG 5' | 3' TCTAGTTTCCCGATGCGTTGG 5' | ||
GRUPO B | |||
Dímeros Secuencia Reverse: | |||
5' TGTCAAGTTAGTGCCTGCTGAACGA 3' | -0.9 | 5' TGTCAAGTTAGTGCCTGCTGAACGA 3' | -0.7 |
||| ¦¦¦ | ||| ¦¦¦ | ||
3' AGCAAGTCGTCCGTGATTGAACTGT 5' | 3' AGCAAGTCGTCCGTGATTGAACTGT 5' |
Horquillas Secuencia Reverse: | |||
/GATTGAACTGT 5' | -0.9 | \GCCTGCTGAACGA 3' | -0.7 |
T ||| | C ||| | ||
/CGTGATTGAACTGT 5' | \TGCTGAACGA 3' | ||
Dímeros entre Primers Forward y Reverse: | |||
5' GTGTGGGACTGGTTGTTGCTGTTG 3' | -1.1 | 5' GTGTGGGACTGGTTGTTGCTGTTG 3' | -0.8 |
||| | ¦¦ ||| ¦ | ||
3' AGCAAGTCGTCCGTGATTGAACTGT 5' | 3' AGCAAGTCGTCCGTGATTGAACTGT 5' | ||
5'GTGTGGGACTGGTTGTTGCTGTTG 3' | -0.8 | 5'GTGTGGGACTGGTTGTTGCTGTTG3' | -0.7 |
¦ ||| ¦ | ¦ ¦ ||| | ||
3' AGCAAGTCGTCCGTGATTGAACTGT 5' | 3' AGCAAGTCGTCCGTGATTGAACTGT 5' | ||
5' GTGTGGGACTGGTTGTTGCTGTTG 3' | -0.7 | 5' GTGTGGGACTGGTTGTTGCTGTTG 3' | -0.7 |
||| | ¦ ¦ ||| | ||
3' AGCAAGTCGTCCGTGATTGAACTGT 5' | 3' AGCAAGTCGTCCGTGATTGAACTGT 5' | ||
5' GTGTGGGACTGGTTGTTGCTGTTG 3' | -0.7 | 5' GTGTGGGACTGGTTGTTGCTGTTG 3' | |
¦ ||| ¦ | ¦ ||| | ||
3' AGCAAGTCGTCCGTGATTGAACTGT 5' | 3' AGCAAGTCGTCCGTGATTGAACTGT 5' | ||
5' GTGTGGGACTGGTTGTTGCTGTTG 3' | -0.7 | ||
¦ ||| ¦ | |||
3' AGCAAGTCGTCCGTGATTGAACTGT 5' | |||
GRUPO E y F | |||
Dímeros Secuencia Forward: | |||
5' CTACCCGGCGTCACATACCA 3' | -4.3 | ||
|||| | |||
3' ACCATACACTGCGGCCCATC 5' |
Se determinó que al asignar valores 1 y 2 como penaltis, los primers dise- ñados son de mayor estabilidad, y al analizarlos en Beacon Designer no for- maron estructuras secundarias, lo que implica que la PCR que se realice con estos primers tendría mejores resultados. Por otro lado, al asignar valores de penaltis mayores a 2, se observa la formación de estructuras secundarias.
Los amplicones por cada Grupo de primers diseñados en Primer3Plus fue- ron analizados en UNAFold. En la Figura 5, los primers del Grupo C y D pre- sentaban un ΔG positivo o nulo, no obstante los valores de Tm están dentro de los valores de referencia, por tanto se determinó que los mismos son ade- cuados para realizar un ensayo de qPCR en el cual la reacción no excede la temperatura de alineamiento de 60ºC, por tanto, esta reacción es efectiva.
A
B
C
Con respecto a los resultados obtenidos, se observó que el uso de penal- tis en el software Primers3plus permitió una búsqueda más adecuada de pri- mers, que poseen un equilibrio entre especificidad y eficiencia para un mejor rendimiento en el ensayo de PCR en tiempo real. A diferencia del diseño de primers mediante el uso del software Primers3 en el cual se observó la forma- ción de estructuras secundarias.
Para diseñar primers se debe tener en cuenta no solo los tiempos de aná- lisis, sino el tiempo de búsqueda que dependerá principalmente de la lon- gitud de la secuencia. En este sentido, diferentes software tienen en sus pa- rámetros o variables, valores predeterminados automatizados y avanzados, dirigidos específicamente al diseño de primers de acuerdo a los tipos de PCR.
La especificad de los primers fueron evaluados mediante el sitio web de Nucleotide-BLAST, los resultados obtenidos se observan en la Tabla 5, todos los primers presentaron un porcentaje de cobertura de 100%, lo cual indica alta similitud de la secuencia con el gen E1 de virus de la Chikungunya.
Se determinó que el E-valor más cercano a la nulidad corresponde al gru- po B con 0.002, seguido del grupo A con 0.024, pero este grupo no coincide con el target que se busca; y por último los grupos C, D, E y F coinciden con el target que se busca, esta medida es más significativa cuanto más cerca de
0.0 esté, (Sánchez, 2012), la cual es la medida más significativa para evaluar la especificidad de primers en el portal Blast del NCBI.
Por lo tanto, todos los grupos, excepto el grupo A, han tenido resultados significativos para el target de búsqueda, óptimos para su uso en la identifica- ción del virus de la Chikungunya mediante PCR en tiempo real.
Grupo | Descripción | Puntaje Max | Puntaje Total | Cobertura | Valor E | Ident | Nº Acceso |
A | Chikungunya virus isolate CHIKV/ Homo apiens/ Brazil/2016/35AP | 44.1 | 44.1 | 100% | 0.024 | 100% | KY704955.1 |
B | Chikungunya virus isolate 13DX-16-0491 E1 gene, partial cds | 48.1 | 48.1 | 100% | 0.002 | 100% | KY131440.1 |
C y D | Chikungunya virus isolate 13DX-16-0491 E1 gene, partial cds | 40.1 | 40.1 | 100% | 0.25 | 100% | KY131440.1 |
E y F | Chikungunya virus isolate 13DX-16-0491 E1 gene, partial cds | 40.1 | 40.1 | 100% | 0.25 | 100% | KY131440.1 |
Existen varias guías para el diseño de primers, en este trabajo se han dise- ñado primers In silico empleando dos tipos de software libres, Primer-BLAST y Primer3Plus, para amplificar la región del gen E1 del virus de la Chikungun- ya empleando los mismos parámetros de entrada para ambos softwares y se determinó la factibilidad, eficacia que presenta el uso de penaltis y la eficien- cia de Primer3Plus que permitió ajustar los parámetros establecidos y el ajus- te de penaltis a un solo parámetro (Tm), para obtener primers más específi- cos sin alterar los parámetros requeridos para llevar a cabo el ensayo de PCR. Los primers obtenidos mediante Primer3Plus son más estables y específicos. Por lo tanto, emplear el uso de penaltis a los diferentes parámetros permitiría mejorar la calidad de los primers y por ende efectivizar un ensayo de qPCR.
Este trabajo ha sido realizado gracias al apoyo académico del Dr. Oscar Cárdenas.
De la Fuente, A. y col.
DORES (Primers). Obtenido de Centro Thierry, A. R., Mouliere, F., Messaoudi, S. E., de Bioinformática del Instituto de Bio- Mollevi, C., Lopez-Crapez, E., Rolet, F., . tecnología Universidad Nacional de Co- . . Nouaille, M. (2014). Clinical validation lombia: http://bioinf.ibun.unal.edu.co/ of the detection of KRAS and BRAF mu- documentos/primers/primer.php tations from circulating tumor DNA. Na-
ne. 7 (12): 925-932. Wang, Y., Tiwari, V. K., Rawat, N., Gill, B. S.,
(1) :31-56. Yang, Y., Kung, T., Hu, C., & Lin, S. (2015). De-