R E S U M E N

El género Enterovirus es un grupo viral que afecta a un amplio rango de hospederos, entre ellos los humanos (especies A, B, C, y D), causan enfermedades respiratorias, gastrointestinales, neurológicas, y otras, y son altamente contagiosos. Los síntomas pueden ser leves o graves. El objetivo del trabajo fue analizar la variación nucleotídica, filogenética y de presión evolutiva de secuencias nucleotídicas del gen VP4 de las cuatro especies que afectan a los humanos. Se emplearon 92 secuencias nucleotídicas disponibles en la base de datos GenBank; éstas se editaron con el software BioEdit y se alinearon con Clustal W; las relaciones filogenéticas se determinaron con MEGA6, y las presiones evolutivas con los algoritmos SNAP y SLAC. Se encontró que la identidad nucleotídica mínima intra-especie fue de 43,2% (especie B) a 72,6% (especie D). Los genotipos más variables por especie fueron EV-71 (A), Echovirus 2 (B), EV-118 (C), y EV-94 (D). El análisis de presión evolutiva mostró que el gen VP4 en las cuatro especies evoluciona bajo presión selectiva negativa. Esto indicaría que la alta tasa mutacional y eventos de recombinación no tienen un rol significativo en la evolución de este gen, debido probablemente a la localización interna de la proteína VP4.

Palabras claves: Enterovirus humanos, secuencias nucleotídicas, gen VP4, análisis filogenético, presión selectiva.

A B S T R A C T

The Enterovirus genus is a viral group that affects a wide host range, including humans (species A, B, C and D), cause respiratory, gastrointestinal, and neurologic disease, among others, and are highly contagious. The symptoms range from mild to severe. The objective of this study was to perform a nucleotidic variation, phylogenetic and selective pressure analyses of the VP4 gene from the four enterovirus species that affect humans. Ninety-two nucleotide sequences (available in the GenBank database) were employed; they were edited with BioEdit software and aligned with Clustal W; the phylogenetic relationships were determined with MEGA6, and the evolutive pressures with SNAP and SLAC algorithms. It was found an intra-species nucleotide identity of at least 43,2% (species B) to 72,6% (species D). The more variable genotypes by species were EV-71 (A), Echovirus 2 (B), EV-118 (C), and EV-94 (D). The selective pressure analysis showed that VP4 gene of the four species evolves by negative pressure. This would indicate that the high mutation rate and recombination events do not have a significant role in the evolution of this gene, probably due to the internal localization of the VP4 protein.

Key words: Human enteroviruses, nucleotide sequences, VP4 gene, phylogenetic analysis, selective pressure.

INTRODUCCCIÓN

Los Enterovirus, pertenecientes a la familia Picornaviridae, comprenden un grupo de virus de tamaño pequeño y RNA como material genético (1). Los viriones son esferoidales, con un diámetro de aproximadamente 30 nanómetros. La cápside es simétrica y posee un arreglo de 60 subunidades llamadas protómeros, dispuestas en una estructura icosaédrica altamente empaquetada; cada protómero está formado por 4 proteínas estructurales denominadas proteínas virales (VP). En la superficie, cada bloque está constituido por VP1, VP2 y VP3, y en la parte interna están relacionadas con VP4 (1). El RNA viral es infeccioso porque una vez dentro de la célula huésped genera todas las proteínas virales requeridas para la replicación viral (1).

Se clasifican en 10 especies: Human enterovirus A, B, C y D; Simian enterovirus, Bovine enterovirus, Porcine enterovirus B; y Human rinhovirus A, B y C. La clasificación se realiza en función a la identidad aminoacídica de proteínas externas compartida por los miembros de una especie; por lo tanto, cada especie enteroviral comprende características en común compartidas por los miembros que los conforman (1).

Los Enterovirus humanos no se limitan a afectar exclusivamente el tracto gastrointestinal o respiratorio; también pueden producir infecciones del sistema nervioso central (meninges y encéfalo), páncreas, corazón y también se los han relacionado con la aparición de diabetes mellitus (2).

La variabilidad genética es una determinante importante en la generación de nuevos tipos o variantes virales. Entre los principales mecanismos generadores de variabilidad en virus RNA se encuentran las mutaciones y recombinaciones. Las mutaciones nucleotídicas pueden darse debido a la incapacidad de la RNA polimerasa de corregir errores durante la replicación, lo que ocasiona elevadas tasas de mutaciones por sustitución (3). La recombinación, originada por intercambio genético entre dos cepas virales al infectar una misma célula, puede originar cepas con características muy distintas a las predecesoras; este mecanismo, aparte de contribuir a la variabilidad genética y la generación de nuevos virus, reduce la carga mutacional (4, 5).

La presión selectiva se define como el grado con que los organismos son favorecidos o desfavorecidos por las mutaciones genéticas (6), pudiendo ocasionar dos escenarios posibles: 1) el cambio en el codón no produce cambio en el aminoácido codificado (denominado sustitución sinónima), y 2) el cambio en el codón produce cambio del aminoácido (denominado sustitución no sinónima). Una presión selectiva positiva implica que la variación aminoacídica que sufrió el gen, favoreció a la especie en su supervivencia. La presión negativa, en cambio, se refiere a que las condiciones de selección negativa para ese gen han permanecido iguales, eliminando los cambios que favoreciesen a la especie (7).

Las reconstrucciones filogenéticas son utilizadas como una herramienta útil para el análisis evolutivo de secuencias génicas, permitiendo la clasificación de los genes estudiados de manera más sencilla y relacionándolos con factores evolutivos de la especie y la historia demográfica de las poblaciones en estudio. Además, se puede obtener información sobre la dinámica de los virus a partir de la estructura de las ramas en un árbol filogenético (8).

En América Latina, existen datos de circulación de las cuatro especies de enterovirus humanos, en mayor porcentaje los genotipos de la especie B (9). En Paraguay no existen datos precisos sobre especies ni genotipos enterovirales asociados a manifestaciones clínicas, por lo que este estudio podría utilizarse como punto de partida para análisis moleculares de enterovirus de importancia en Salud Pública en nuestro país.

MATERIALES Y MÉTODOS

Obtención de secuencias nucleotídicas

Se obtuvieron en total 92 secuencias nucleotídicas de referencia del gen de VP4 de enterovirus humanos (especies A, B, C, y D), de la base de datos GenBank (http://www.ncbi.nlm.nih.gov); cada gen estaba conformado por 67 codones (201 nucleótidos); cada secuencia correspondió a un genotipo representativo por especie, según asignación taxonómica en la base de datos de los picornavirus (http://www.picornaviridae.com/). El número de secuencias obtenidas por especie fue de 17 para la especie A, 55 para la B, 17 para la C y 3 para la D.

Análisis filogenético y de identidad nucleotídica

Las secuencias nucleotídicas obtenidas fueron manualmente editadas con BioEdit v.7.0.5 (10). El alineamiento múltiple de secuencias fue realizado con Clustal W (11). La reconstrucción filogenética fue realizada con el método de neighbor-joining, utilizando el método de Kimura 2-parámetros como modelo de sustitución nucleotídica, y análisis de bootstrap de 1.000 réplicas con el programa MEGA v6 (12). Para el cálculo de identidad nucleotídica, se realizó un análisis de distancia entre pares de secuencias (pairwise distance) dentro de una misma especie, utilizando el programa MEGA v6 (12); a partir de las distancias, se dedujeron las identidades (expresadas en porcentajes).

Análisis de las presiones selectivas

Para el análisis de presión selectiva, se utilizó el método de Nei y Gojobori (7). Los datos fueron analizados utilizando el programa SNAP (13), basado en el lenguaje de programación Perl; para ello, se partió de un alineamiento múltiple de secuencias, y se contabilizó el número de cambios sinónimos y no sinónimos para cada codón. Al hallar la relación entre las sustituciones no sinónimas (dn) y sinónimas (ds), i.e. dn/ds, si el valor es 1, la presión de selección es neutra; si el valor es menor a 1 implica una presión de selección negativa, y si es mayor a 1, una presión selectiva positiva (7).

La determinación del número de codones bajo presión selectiva positiva, negativa o neutra se estimó con el algoritmo “Single likelihood ancestor counting” (SLAC) de la interfase Datamonkey (14).