INTRODUCCIÓN
Las enfermedades caninas transmitidas por vectores constituyen un peligro para la salud animal. Algunos patógenos causantes de estas enfermedades son de gran preocupación zoonótica y constituyen un grave peligro para la salud humana en todo el mundo. 1
Ehrlichia canis es el agente de la ehrlichiosis monocítica canina en perros, lobos y chacales y se ha reportado infectando a humanos. 2 La infección por Anaplasma platys causa trombocitopenia cíclica infecciosa en perros. El patógeno también se ha identificado en un amplio abanico de hospedadores, que incluye a los humanos, además de los perros. 3 Las infecciones por A. platys pueden llegar a ser graves o mortales, sobre todo cuando hay coinfecciones con otros patógenos transmitidos por garrapatas, como E. canis. La presencia de A. platys y E. canis ha sido descrita previamente en Cuba, 4,5 pero la información relativa a su prevalencia y diversidad genética sigue siendo escasa.
La borreliosis de Lyme es la enfermedad zoonótica transmitida por garrapatas más prevalente en el hemisferio norte, causada por Borrelia burgdorferi sensu lato.6 La anaplasmosis granulocítica humana es otra enfermedad transmitida por garrapatas con importancia para la salud pública, causada por Anaplasma phagocytophilum.
La rickettsiosis es una enfermedad ampliamente distribuida por todo el mundo; varias especies son patógenos zoonóticos emergentes o reemergentes transmitidos por artrópodos hematófagos. 7
La hepatozoonosis canina es una enfermedad protozoaria transmitida por garrapatas y causada por Hepatozoon canis y Hepatozoon americanum. Hasta la fecha, en Cuba solo hay un breve informe sobre la detección por PCR de Hepatozoon spp. en muestras de sangre de perros, y no se ha informado de la confirmación de la especie, ni de su caracterización molecular. 8
En Cuba, el género Babesia se ha informado en perros, pero su diagnóstico se ha basado en la observación de signos clínicos, hallazgos hematológicos, observación directa mediante frotis de sangre y la detección de anticuerpos a través de métodos serológicos. González et al.9 realizaron la primera detección molecular de B. vogeli infectando garrapatas Rhipicephalus sanguineus s.l. colectadas de perros con dueño y plantearon la necesidad de la detección de las especies de Babesia infectando perros. El presente estudio tuvo como objetivo determinar la prevalencia de patógenos zoonóticos (Anaplasma platys, Anaplasma phagocytophilum, Borrelia burgdorferi sensu lato, Ehrlichia canis y Rickettsia spp.) y no zoonóticos (Hepatozoon canis y Babesia spp.) transmitidos por garrapatas en perros sin dueño de La Habana.
MÉTODOS
Diagnóstico molecular, prevalencia e importancia de los patógenos zoonóticos transmitidos por garrapatas
Recogida de muestras y extracción de ADN
Se tomaron muestras de sangre de un total de 100 perros seleccionados al azar de 10 municipios de la provincia La Habana. La extracción del ADN se realizó utilizando el juego de reactivos Wizard Genomic DNA Purification Kit (Promega, Madison, WI, EE. UU.). Se extrajo una muestra de hasta 10 garrapatas por perro infestado. La clasificación taxonómica se realizó según la clave taxonómica descrita por Estrada-Peña et al. 10
Amplificación y secuenciación por PCR
Para verificar la presencia de ADN amplificable en las muestras, se realizó un ensayo de qPCR en tiempo real para el gen de la gliceraldehído-3-fosfato deshidrogenasa (gapdh) 11) Todas las muestras de ADN se analizaron mediante qPCR en tiempo real utilizando los cebadores y sondas descritos para A. phagocytophilum,12A. platys,13B. burgdorferi s.l., 14) E. canis15 y Rickettsia spp. 16).Se realizó la secuenciación del ADN en muestras seleccionadas al azar y se utilizaron en ensayos de PCR con cebadores específicos de género y especie para Ehrlichia/Anaplasma spp. (gen ARNr 16S), 17A. platys (gen groEL), 18E. canis (gen gltA), 19 y Rickettsia spp. (genes ompA, ompB y htrA) 20
Análisis de secuencias
Los productos de la PCR se purificaron y se clonaron utilizando los reactivos pCR 2.1 de clonación TOPO TA (Thermo Fisher Scientific, Dreieich, Alemania) seguido de la transformación en células competentes Escherichia coli Top 10F'. El ADN del plásmido recombinante se envió para secuenciación con cebadores universales del gen M13 (Microsynth, Balgach, Suiza). Las secuencias obtenidas se analizaron mediante BLAST: Basic Local Alignment Search Tool (http://blast.ncbi.nlm.nih.gov/Blast.cgi). La traducción teórica de las secuencias de nucleótidos a secuencias de aminoácidos se realizó con la herramienta ExPASy translate (http://www.expasy.org) y las secuencias de proteínas se alinearon con ClustalW, incluido en el paquete BioEdit v.7.0.0 (Ibis Biosciences, Carlsbad, CA, EE. UU.).
Análisis filogenético
El análisis filogenético se realizó con el paquete de software Molecular Evolutionary Genetics Analysis versión 7.0 (MEGA7), 21 utilizando el método neighbor-joining. Las secuencias se alinearon utilizando MAFFT configurado para la mayor precisión y se identificaron las regiones conservadas. 22 Tras la alineación, las regiones ambiguas se eliminaron con la versión 0.91 b de Gblocks 23 Para la construcción de los árboles filogenéticos, se seleccionó el modelo que mejor se ajustaba a la evolución de la secuencia basándose en el criterio de información de Akaike Corregido (cAIC) y el criterio de información bayesiano (BIC) implementados en MEGA7.
Análisis de datos
Los datos obtenidos se compilaron y analizaron con el software Excel 2016 (Microsoft Corporation, WA, EE. UU.), y el análisis estadístico se realizó con el software R. 24 Las tasas de prevalencia de A. platys, E. canis, Rickettsia spp. y coinfecciones con intervalos de confianza (IC) del 95 % se calcularon mediante un enfoque bayesiano basado en la distribución beta (s + 1; n-s + 1), donde s = positivos; n = animales analizados.
Detección y caracterización molecular de Hepatozoon canis
Toma de muestras de sangre y examen microscópico
Se tomaron muestras de sangre total de 80 perros sin dueños que vivían en diferentes lugares de La Habana. Se prepararon 2 frotis sanguíneos finos, se tiñeron con Giemsa 25 y se examinaron para detectar la presencia de gametos de Hepatozoon utilizando un microscopio de luz (Carl Zeiss Microscopy GmbH, Jena, Alemania), con lente de inmersión en aceite a un aumento final 1000 X.
Extracción de ADN y ensayos de PCR en tiempo real
El ADN se extrajo de sangre total utilizando el juego de reactivos Wizard Genomic DNA Purification Kit (Promega, Madison, WI, EE.UU.). Se realizó un ensayo de PCR en tiempo real para el gen gapdh.26 La presencia de ADN de H. canis se evaluó mediante un ensayo de PCR en tiempo real con SBYR Green. 27 Se utilizó como control positivo una porción sintética de ADN de doble cadena (600 pb) del gen del ARNr 18 S de H. canis (GeneArt String DNA, Thermo Fisher Scientific Inc., Waltham, MA, EE. UU.).
Secuenciación y análisis filogenético
Las muestras positivas a la PCR de H. canis con valores de ciclo umbral (Ct) bajos se utilizaron para amplificar el gen ARNr 18S (~1750 pb). 28 Se enviaron para la secuenciación automatizada (Eurofins, Toronto, Canadá).
Las secuencias se enviaron a GenBank con los números de acceso MN393910-MN393913 y se compararon con secuencias del gen ARNr 18 S depositadas en el Gen Bank utilizando la herramienta nBLAST. Para el análisis filogenético, las secuencias se alinearon con las secuencias conocidas de Hepatozoon sp. utilizando MAFFT v.7 configurado para la máxima actividad (MAFFT con ajustes por defecto). 29 El modelo de evolución de la secuencia que mejor se ajustaba se seleccionó basándose en el criterio de información de Akaike corregido (cAIC) y el criterio de información bayesiano (BIC), tal y como se implementó en el análisis de genética evolutiva molecular (MEGA) v.7.0.14. 30 Se eligió para la reconstrucción del árbol el modelo de 3 parámetros de Tamura. 31 El árbol filogenético se reconstruyó utilizando el método neighbor-joining (NJ) implementado en MEGA.
Análisis estadístico
Todos los datos obtenidos se compilaron en el software Excel 2016 (Microsoft Corporation, Redmond, WA). La concordancia entre los resultados de frotis de sangre y la PCR en tiempo real del gen ARNr 18 S se comprobó mediante la estadística kappa (κ) basada en el esquema propuesto por Altman et al. 32 La prevalencia observada, los estadísticos kappa, los intervalos de confianza (IC) de 95 % y el análisis estadístico se realizaron utilizando el software R. 24 La prueba U de Mann-Whitney se realizó utilizando la función mwu del paquete sjstats v0.17.5 y el análisis de concordancia se realizó utilizando el paquete fmsb v0.6.3 del software R.
Ocurrencia de Babesia spp. en perros sin dueño
Toma de muestras y extracción de la sangre
Se colectaron muestras de sangre completa de 60 perros sin dueño procedentes de 11 municipios de la provincia La Habana. La sangre se obtuvo mediante punción de la vena yugular, por sistema de vacío (BD Vacutainer®), con 7,2 mg de EDTA K2 como anticoagulante.
Colecta de garrapatas
De cada animal se extrajo, al menos, 10 de las garrapatas presentes y se agruparon por especie, estadio y sexo. Para la identificación se empleó la clave taxonómica publicada por Nava et al.33
Diagnóstico por frotis sanguíneo
Las formaciones compatibles con Babesia spp. se identificaron como inclusiones intracitoplasmáticas, ubicadas en los eritrocitos, en forma de peras, simples, pareadas o múltiples; según lo descrito por Lempereur et al.34 extracelulares, según Coralic et al.35
Extracción y calidad del ADN
Para la extracción del ADN se utilizó el juego de reactivos Genomic DNA Purification Kit (Promega, Madison, EE. UU.). Para detectar la presencia de posibles inhibidores en el ADN extraído, se amplificó un fragmento del gen gapdh (gliceraldehido 3 fosfato deshidrogenasa). 36
Diagnóstico de Babesia spp.
Se seleccionaron los cebadores descritos por Birkenheuer et al.36 con los que se amplifica en un primer ensayo de PCR el ARNr 18 S de Babesia y en un segundo ensayo de PCR se amplifica una región de aproximadamente 340 pb dentro de este gen. Los resultados del PCR se aplicaron en geles de agarosa (Sigma-Aldrich, St. Louis, MI, EE. UU.) al 2 % en tampón TBE 0,5X, se tiñeron con Bromuro de Etidio (0,5 μg/mL) y se visualizaron en un transiluminador de luz ultravioleta Macro Vue (Pharmacia Biotech Inc., EE. UU.). En todos los casos se utilizó un marcador de peso molecular de 100 pb (Promega, Madison, EE. UU.).
Análisis estadístico
Los datos se agruparon conformando una base de datos en Microsoft Excel 2016 (Microsoft Corporation, Redmond, WA). Se utilizó la prueba exacta de Fisher para determinar si existía asociación entre las variables infestación por garrapatas y el sexo del animal, con la infección por Babesia spp. El análisis de los datos obtenidos se realizó con el empleo del software R, Package stats version 3.6.0. 24 con un 95 % de intervalo de confianza (IC).
RESULTADOS
Diagnóstico molecular, prevalencia e importancia de los patógenos zoonóticos transmitidos por garrapatas
De las 100 muestras de sangre, 85 (85 %; IC del 95 %: 77, 88-92, 12) resultaron positivas para al menos un patógeno. E. canis fue el patógeno más prevalente 62 % (IC 95 %: 52,32-71,68), seguido de A. platys 40 % (IC 95 %: 30,23-49,77) y Rickettsia spp. 27 % (IC 95 %: 18,15-35,85). Los perros estaban coinfectados más frecuentemente por E. canis y A. platys en 28 (28 %; IC 95 %: 19,05-36,95), seguidos por E. canis y Rickettsia spp. en 13 (13 %; IC 95 %: 6,29-19,71), A. platys y Rickettsia spp. en 11 (11 %; IC 95 %: 4,76-17,24), y se detectaron infecciones mixtas triples en 8 perros (8 %; IC 95 %: 2,59-13,41). Se observó infestación por garrapatas en 57 de los 100 perros muestreados. Todas se identificaron morfológicamente como R. sanguineus (111 hembras, 131 machos y 6 ninfas). La presencia de garrapatas no se asoció estadísticamente con la aparición de infecciones positivas a la PCR (p = 0,115): E. canis (p = 0,269), A. platys (p = 0,268), Rickettsia spp. (p = 0,814) y coinfecciones (p = 0,411).
El análisis de las secuencias del gen del ARNr 16 S obtenidas de los aislados cubanos de E. canis (1434 pb) y A. platys (1431 pb) reveló identidades > 99 % con varias secuencias de E. canis (LC269822, EF139459) y A. platys (EU106856, CP000107) disponibles en el Gen Bank. Las secuencias de nucleótidos obtenidas en el presente estudio no eran 100 % idénticas entre sí y pueden representar variantes locales que existen dentro de la región estudiada.
Las secuencias obtenidas de E. canis-gltA fueron 100 % idénticas entre sí y a las secuencias notificadas de China (CP025749), Filipinas (LC428206) y Zambia (LC373038), mientras que, para A. platys-groEL, las secuencias fueron 100 % idénticas entre sí y a las notificadas de la República Democrática del Congo (AF478129), Japón (AY077621) y Venezuela (AF399916) Las secuencias de Rickettsia spp.-htrA (434 pb) fueron 100 % idénticas entre sí y mostraron > 99 % de identidad con las secuencias del gen del antígeno de superficie de 17 kDa de los aislados de referencia de R. felis notificados en México (GU447234) y Estados Unidos (CP000053). Para las secuencias de los genes E. canis-gltA, A. platys-groEL y R. felis-htrA no se observó ninguna variación nucleotídica entre los amplicones de PCR secuenciados.
El análisis filogenético basado en las secuencias del gen ARNr 16S se agrupó en 2 clados principales de Anaplasma spp. y Ehrlichia spp. El árbol filogenético resultante reveló que los aislados cubanos de A. platys y E. canis se agrupaban estrechamente con otras cepas de A. platys y E. canis notificadas en todo el mundo. El genotipo cubano de A. platys-groEL se clasificó estrechamente en el clúster de A. platys agrupado con otras cepas aisladas de diferentes especies hospedadoras de todo el mundo, con un valor bootstrap del 100 %. La cepa cubana de E. canis-gltA se agrupó dentro del clado de E. canis formado por cepas aisladas de diferentes especies hospedadoras de todo el mundo, con un apoyo del 100 % del valor bootstrap.
Detección y caracterización molecular de Hepatozoon canis
El examen microscópico de los frotis de sangre reveló la presencia de gametos de Hepatozoon spp. en 8/80 (10,0 %; IC del 95 %: de 4,80 % a 18,0 %) perros. Los análisis por PCR detectaron la presencia de ADN de H. canis en 38/80 (47,5 %; IC 95 %: de 36,8 % a 58,4 %) muestras. Todas las muestras positivas por microscopía se confirmaron como infectadas por H. canis mediante PCR y secuenciación de ADN. Se observó una concordancia "pobre" entre los resultados obtenidos por el examen microscópico y el PCR, con un valor κ de 0,20 (IC 95 %: de 0,07 a 0,33). Las secuencias obtenidas para los tres aislados de H. canis fueron muy similares entre sí (de 99,3 % a 100 %) y > 99 % idénticas a las de los aislados de H. canis de República Checa (KX712124), Brasil (AY461375 y España (AY150067). El análisis filogenético produjo un árbol en el que los aislados cubanos de H. canis se situaron dentro de un clado con un fuerte apoyo de bootstrap (93 %) que contenía todos los demás aislados de H. canis con exclusión de todas las demás especies del género.
Ocurrencia de Babesia spp. en perros sin dueño
A la exploración clínica el 51,6 % animales (31/60) estaban infestados con garrapatas. Se colectaron 69 garrapatas adultas y 3 ninfas de Rhipicephalus sanguineus s.l. En 10 de los perros donde se encontró R. sanguineus s.l., también se colectaron 15 garrapatas adultas y 2 ninfas de Rhipicephalus microplus.
Se observó que 11,6 % animales (7/60) presentaron formaciones intraeritrocíticas compatibles con Babesia spp. Se visualizaron microorganismos intraeritrocíticos en forma de pera en pares o simples. En un animal (1,6 %) se observaron merozoitos extracelulares.
Como resultado del diagnóstico molecular mediante PCR, 20 % (12/60) de los animales resultaron positivos a la presencia de Babesia spp, con una banda de 340 pb, lo que se corresponde con lo descrito para este ensayo por Birkenheuer et al.24 El estudio de asociación entre el sexo de los animales y la presencia de garrapatas con la infección por Babesia spp, mostró que no existe asociación entre las variables evaluadas para un valor de p = 1.
DISCUSIÓN
Hasta donde los autores conocen, este es el primer estudio dirigido a investigar los patógenos transmitidos por garrapatas en perros sin dueño de La Habana. El 85 % de los perros resultaron positivos con los ensayos de PCR para al menos 1 de los patógenos en estudio. Las infecciones detectadas tienen un potencial zoonótico relevante, ya que se han notificado infecciones humanas por E. canis, A. platys y Rickettsia spp. en Venezuela, 37 Granada 38 y España, 39 respectivamente.
En este estudio se observó una alta prevalencia de infecciones mixtas; sin embargo, los perros no mostraron signos clínicos consistentes con las infecciones por E. canis, A. platys y Rickettsia spp. La mayoría de los perros eran asintomáticos, lo que puede reflejar la existencia de infecciones crónicas, subclínicas o leves que dificultan el diagnóstico clínico de los perros infectados.
El análisis filogenético, basado en las secuencias de ADNr 16 S, reveló que las cepas de A. platys y E. canis estaban estrechamente agrupadas con otras cepas de estos patógenos aisladas de perros de todo el mundo. Se observó un perfil genético altamente conservado para las cepas de A. platys y E. canis basado en el análisis de las secuencias parciales de los genes groEL y gltA, respectivamente. El análisis de las secuencias de los genes E. canis-gltA y A. platys-groEL realizado en 3 cepas cubanas reveló una identidad del 100 %, a pesar de que las muestras analizadas se obtuvieron de zonas diferentes. Los alineamientos de las secuencias y los árboles filogenéticos sugirieron una escasa diversidad genética y una evolución homogénea dentro de las cepas de A. platys y E. canis, basándose en la estrecha identidad entre sus secuencias de ADNr 16 S, groEL y gltA. Los resultados obtenidos concuerdan con resultados anteriores sobre una ligera variación genética entre los genes secuenciados de diferentes cepas de A. platys y E. canis.40
No se identificó la presencia de ADN de A. phagocytophilum y B. burgdorferi s.l. en ninguna de las muestras de sangre analizadas. Estos resultados están en consonancia con la ausencia del principal vector de estos patógenos, que son garrapatas duras distintas de R. sanguineus, que nunca han sido reportadas en Cuba. 41
De las 100 muestras, 27 resultaron positivas, en la qPCR de detección de especies de Rickettsia, y el posterior análisis de la secuencia identificó la presencia de R. felis, un miembro del grupo de las Rickettsias de la fiebre manchada (SFGR), por primera vez en perros de La Habana. Un estudio anterior realizado por Noda et al. 42 describió la detección de "Candidatus Rickettsia amblyommii" en la garrapata Amblyomma mixtum mediante PCR, lo que constituye el primer informe de un miembro de la SFGR en Cuba. La alta prevalencia de R. felis encontrada en los perros analizados pone de relieve la importancia sustancial de este patógeno para la salud humana, ya que la Rickettsiosis se ha convertido en un problema reemergente en todo el mundo. Estos hallazgos indican la necesidad de realizar más estudios sobre la presencia de Rickettsia spp. en R. sanguineus de Cuba.
En el presente estudio se informa por primera vez la deteción de H. canis en perros sin dueños de La Habana. El examen microscópico de los frotis sanguíneos reveló una baja incidencia de infecciones por Hepatozoon spp. (10 %) en la población de perros examinados. El ensayo de PCR en tiempo real del gen ARNr 18 S, reveló que el 47,5 % de los perros estaban infectados por H. canis, demostrándose una escasa concordancia entre ambas pruebas. Este resultado demuestra que el examen microscópico de frotis sanguíneos es un método ineficaz, comparado con la PCR para la detección de Hepatozoon spp. en animales infectados, conclusión inferida por estudios anteriores realizados en otras regiones geográficas. 43
La prevalencia de la infección por H. canis en perros sin dueños en La Habana basada en la PCR fue similar a la de informes anteriores basados en el análisis molecular en perros domésticos de Brasil (53,3 %), 44 República Checa (50 %), 45 e Italia (57,8 %). 46 Sin embargo, la prevalencia de la infección reportada para La Habana (47,5 %) fue mayor comparada con la descrita en estudios realizados en Haití (9,3 %) 47 y San Cristóbal y Nieves (6 %). 48 Las diferencias en la prevalencia de H. canis pueden deberse a varios factores, como la distribución geográfica y la densidad de población del vector, las condiciones climáticas, la metodología de muestreo utilizada, el estado inmunitario del hospedador y la población canina a la que se dirige. 49
Los resultados de los análisis filogenéticos revelaron que H. canis de los perros de La Habana pertenecía al mismo clado que todos los demás aislados de H. canis, y con un apoyo estadístico muy fuerte para la separación de este clado de todas las demás especies de Hepatozoon. El presente estudio proporciona la primera caracterización molecular de H. canis en perros sin dueños de La Habana. Los resultados confirman mediante PCR en tiempo real una alta prevalencia de H. canis en perros de zonas urbanas de La Habana, Cuba. Son necesarios más estudios para determinar la prevalencia de H. canis en diferentes poblaciones caninas y en vectores potenciales de diferentes zonas geográficas de Cuba.
En la detección de Babesia spp. el frotis sanguíneo mostró un bajo desempeño con respecto al PCR. Resultados similares fueron obtenidos por O’Dwyer et al.,50 que reportaron el 8 % de las muestras positivas por PCR y solo en el 2 % se observaron los merozoitos de Babesia spp.
Singh et al. 51 no encontraron diferencias significativas (p ≤ 0,05) entre el sexo y la infección por Babesia spp. en perros de la India. Estos autores plantean que existen otros factores de riesgo como la edad (animales jóvenes) y la época del año (verano) que si influyen de forma directa en la presencia de este patógeno en perros.
En el caso de la infestación por garrapatas, Galay et al. 52 reportaron que no existía asociación entre esta variable y la infección por Babesia spp.; sin embargo, los perros sin dueño pueden ser vehículos para las garrapatas, contribuyendo así a la propagación de las hemoparasitosis, por lo que la prevención de las enfermedades transmitida por garrapatas se logra eliminando la posibilidad de exposición al vector mediante un control activo del ambiente. Hasta donde los autores conocen este es el primer estudio de detección Babesia spp. en perros sin dueño de La Habana y demuestra la necesidad de investigaciones futuras dirigidas a determinar la especie, con estudios adicionales sobre los vectores que la pueden transmitir.
Conclusiones
En el presente estudio se determinó la prevalencia de patógenos zoonóticos y no zoonóticos transmitidos por garrapatas en perros sin dueño de La Habana. Estos resultados constituyen el primer reporte de Rickettsia felis en perros sin dueño y demuestran la alta prevalencia de patógenos transmitidos por garrapatas, con potencial zoonótico, teniendo en cuenta, además las coinfecciones. Por primera vez se tiene evidencia molecular de Hepatozoon canis y Babesia spp. en perros sin dueño de La Habana. Estos resultados son de gran importancia para la vigilancia de las enfermedades transmitidas por vectores en los animales de compañía, y demuestran la necesidad de estudios posteriores sobre la prevención de la transmisión y expansión de las enfermedades que provocan, siendo de gran valor para orientar medidas de prevención y control dirigidas a evitar su diseminación.