Sr. Editor. En los últimos años en el Perú se presentó un alarmante aumento de casos por varicela complicada en población pediátrica. El laboratorio de referencia nacional de Sarampión y Rubéola del Instituto Nacional de Salud (INS), como parte de sus funciones de diagnóstico especializado, conservó muestras clínicas provenientes de pacientes diagnosticados con varicela. Las muestras correspondieron a suero sanguíneo e hisopado de lesión vesicular obtenidas de casos clínicos de un brote epidémico ocurrido en las regiones Lima, La Libertad, Loreto, Cusco y Tumbes. Las muestras fueron conservadas en refrigeración entre dos y ocho °C, hasta su posterior procesamiento por ensayo por inmunoabsorción ligado a enzimas (ELISA) para determinar la respuesta de anticuerpos IgM, IgG; y por reacción en cadena de la polimerasa (PCR) en tiempo real para la detección del virus. Las muestras positivas fueron sujetas a un posterior análisis de caracterización molecular con el objetivo de genotipificar los virus de varicela zoster (VZV).

Entre el 2016 y 2017 fueron remitidas al laboratorio de Sarampión y Rubéola un total de 198 muestras de suero y 177 de hisopado de lesión vesicular. Del total de muestras 126 resultaron IgM positivos, y 104 fueron positivas en el PCR en tiempo real, las cuales fueron inoculadas en cultivos puros de células de Vero y Neuroblastoma. Las células experimentaron un efecto citopático como consecuencia de la infección producida por VZV. Posteriormente los aislamientos fueron cosechados y almacenados a -20°C; a partir de estos se realizaron las extracciones de ADN usando el kit comercial QIAamp DNA Mini (QIAGEN), seguidamente se realizaron las amplificaciones por PCR convencional con el uso de cebadores (primers) desarrollados por Li & Zhu para VZV y para las siguientes regiones: ORF1, ORF 21, ORF 22, ORF 38, ORF 54, ORF 62 y R5 ¹1. Li Y, Zhu B. Genotyping of clinical varicella-zoster virus isolates collected from Yunnan in Southwestern China. Biomed Rep. 2016;4(2):209-14. DOI: 10.3892/br.2015.562.
https://doi.org/10.3892/br.2015.562... .

Los resultados fueron observados en geles de agarosa al 2%, los productos de PCR se purificaron usando un kit comercial de ADN genómico QIAquick Gel Extraction (QIAGEN) y fueron secuenciados por el método de Sanger y procesados en un analizador genético Applied Biosystems 3500. Se obtuvieron cromatogramas de ambas hebras (forward y reverse) los cuales fueron editados y ensamblados con el programa SeqTrace versión 0.9.0 ²2. Stucky BJ. SeqTrace: A Graphical Tool for Rapidly Processing DNA Sequencing Chromatograms. J Biomol Tech JBT. 2012; 23(3):90-93. DOI: 10.7171/jbt.12-2303-004.
https://doi.org/10.7171/jbt.12-2303-004... , las secuencias obtenidas fueron luego contrastadas con secuencias de genomas completos o parciales depositadas en la base de datos GenBanK³3. Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Ostell J, Pruitt KD, et al. GenBank. Nucleic Acids Res. 2018; 46(D1):D41-D47. DOI: 10.1093/nar/gkx1094.
https://doi.org/10.1093/nar/gkx1094... , de la que se descargaron 165 para cada región estudiada (ORFs). Se realizaron alineamientos de secuencias con el programa MAFFT versión 7 ⁴4. Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002; 30(14):3059-66. y la concatenación de todas las regiones con el programa Sequence Matrix 1.8 ⁵5. Vaidya G, Lohman DJ, Meier R. Sequence Matrix: concatenation software for the fast assembly of multi-gene datasets with character set and codon information. Cladistics. 2011; 27(2):171-80.. El alineamiento total estuvo conformado por una matriz de 171 secuencias y 3456 pares de bases (pb). Se identificaron los cambios a nivel nucleotídico y aminoacídico de las muestras peruanas con respecto de la vacuna VariVax. Se realizó una reconstrucción filogenética con la metodología Neighbor-Joining (NJ) con la opción pairwise deleción, modelo de sustitución Kimura 2 parámetros (K2P), y 1000 réplicas bootstrap. Estas dos últimas operaciones fueron realizadas con el programa MEGA 7 ⁶6. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol. 2016; 33(7):1870-4..

El análisis de las secuencias nucleotídicas obtenidas a partir de las extracciones de ADN de los cultivos celulares, permitió determinar la presencia del VZV en las muestras, confirmándose los resultados previos de PCR en tiempo real. Como parte de la genotipificación de VZV se obtuvo un árbol filogenético NJ en el que se distinguen los clados a los que pertenecen las muestras analizadas (Figura 1). Se registraron 29 cambios nucleotídicos puntuales (SNPs) que se traducen en 13 cambios no sinónimos, es decir, cambios aminoácidicos con respecto de la vacuna circulante VariVax, resaltada como triangulo negro en el clado 2 de la Figura 1.

Resaltamos que la topología del árbol NJ se asemeja (aunque no es igual) a la obtenida por otros autores con datos de genomas completos, dado que se distinguen los clados 1,2,3,4,5,VI,VIII y IX como grupos independientes, asimismo, los soportes bootstrap de cada clado superan el 50 % considerado como límite ⁷7. Norberg P, Depledge DP, Kundu S, Atkinson C, Brown J, Haque T, et al. Recombination of Globally Circulating Varicella-Zoster Virus. J Virol. 2015; 89(14):7133-46.^,88. Zell R, Taudien S, Pfaff F, Wutzler P, Platzer M, Sauerbrei A. Sequencing of 21 Varicella-Zoster Virus Genomes Reveals Two Novel Genotypes and Evidence of Recombination. J Virol. 2012; 86(3):1608-22.. La mayoría de las cepas colectadas entre 2016 y 2017 se agrupan en el clado 1 (en rojo en la Figura 1), el mismo que está conformado por cepas circulantes en Alemania, Suecia, Reino Unido, EE.UU. y forman un sub-clado independiente; en tanto, solo una cepa pertenece al clado 6 (en verde en la Figura 1).

Las muestras peruanas aisladas entre 2016 y 2017 pertenecen a dos grupos genéticamente diferenciados (clados 1 y VI). Determinamos que las muestras peruanas con respecto a la vacuna VariVax presentan 29 cambios nucleotídicos y 13 a nivel aminoacídico. Las siete regiones empleadas para el análisis filogenético permiten caracterizar molecularmente las muestras peruanas almacenadas en el laboratorio de Sarampión y Rubeola; y podemos afirmar que corresponden a VZV. Sin embargo, mediante la secuenciación de alto rendimiento de todo el genoma, en el contexto de grandes virus de ADN, los datos recientes sugieren que los VZV también podrían considerarse candidatos para el modelado de la dinámica evolutiva de las poblaciones de estas especies y la diversidad dentro del huésped.

Agradecimientos

A Sara Martha Gordillo Vilchez e Irma Atalaya Atalaya.

REFERENCIAS BIBLIOGRÁFICAS

Fuente de financiamiento

Citar como

Fechas de Publicación

Histórico