Bronquitis infecciosa en gallinas ponedoras

La industria avícola es el sector productivo de mayor crecimiento en el mundo, contribuyendo con la nutrición humana y generando ingresos para las economías familiares (Farrell, 2013). En Ecuador tiene gran proyección económica, debido a que la carne de pollo y los huevos de gallinas son fuentes de nutrientes con precios accesibles en el mercado (Aillón, 2012; Corporación Nacional de Avicultores del Ecuador - CONAVE-, 2020a). El consumo de huevos per cápita en el año 2018 fue de 218 unidades y ascendió a 226 por persona/año durante 2019 (Corporación Nacional de Avicultores del Ecuador -CONAVE-, 2020b). Según el censo avícola del año 2015, existen 310 granjas de aves ponedoras comerciales que se ubican mayoritariamente en la región del centro (Pastaza, Cotopaxi, Tungurahua y Chimborazo) y en la provincia de Manabí (Instituto Nacional de Estadística y Censos -INEC-, 2015). No obstante, la provincia de Tungurahua es la que cubre el 45% de la demanda de huevos de mesa en el país (Ministerio de Agricultura y Ganadería -MAG-, 2021).

El crecimiento de la avicultura ecuatoriana es resultado del trabajo colaborativo de los diferentes eslabones de esta cadena productiva (Agencia Ecuatoriana de Aseguramiento de la Calidad del Agro - AGROCALIDAD-, 2013). Sin embargo, el progreso de la avicultura puede verse afectado por enfermedades que influyen negativamente en la producción y, en consecuencia, la identificación y el control de estas resultan fundamentales para evitar pérdidas económicas (Bermudez, 2008).

La bronquitis infecciosa aviar es una de las enfermedades respiratorias más frecuentes y problemáticas dentro de las granjas de aves ponedoras comerciales en el mundo (Alexander & Senne, 2008; Cavanagh & Gelb, 2008). El agente causal es el IBV(Infectious bronchitis virus) perteneciente al género Gammacoronavirus de la subfamilia Orthocoronavirinae, familia Coronaviridae y orden Nidovirales (International Committee on Taxonomy of Viruses -ICTV-, 2020). El IBV es pleomórfico y posee un genoma de ARN de cadena simple con polaridad positiva y de longitud aproximada de 27,6 kilobases (Feng et al., 2017; Rafique et al., 2018). El genoma codifica para cuatro proteínas estructurales: la proteína de espícula (S), la glicoproteína de membrana (M), la proteína de la nucleocápside y la pequeña proteína de membrana (E) (Abro, 2013; Marandino, 2013). Además, incluye, en sus extremos, regiones no traducidas 5’ y 3’ UTR de baja variación, por lo que se suelen utilizar para determinar la presencia del IBV en muestras de granjas (Acevedo-Beiras, 2017). La entrada delIBVa la célula hospedadora se produce debido a la unión de la proteína S de la envoltura viral a los receptores específicos de la membrana celular (Shang et al., 2018). Esta proteína S se sintetiza como un precursor proteico (S0) de 1162 aminoácidos, que por modificación postraduccional genera dos subunidades funcionales, S1 y S2 (Abro, 2013; Liu et al., 2009a; Rafique et al., 2018). Las regiones hipervariables de S1 poseen una alta tasa de mutación, siendo no solo responsables de la antigenicidad viral, sino también de la mayoría de las variaciones genéticas del virus (AcevedoBeiras, 2010; Dolz et al., 2008; Montassier, 2010).

Debido a la existencia de diversas cepas del IBV, la instauración de un programa vacunal efectivo para la prevención y control de la bronquitis infecciosa constituye un verdadero desafío (de Wit et al., 2011). La información científica en Ecuador se limita a escasas publicaciones enfocadas al monitoreo en las aves de traspatio mediante la utilización de las técnicas de ELISA en el noroeste, hemoaglutinación en Galápagos y PCR en tiempo real en la parroquia de Puéllaro (García, 2016; Hernandez et al., 2006; Soos et al., 2008). Las investigaciones mencionadas y las vacunas empleadas en Ecuador demuestran la presencia del IBV; sin embargo, no existen publicaciones científicas que expongan de manera fehaciente la presencia de las variantes que pueden estar circulando. Por lo tanto, el objetivo de la investigación fue caracterizar, mediante técnicas moleculares, las cepas del IBV circulantes en las explotaciones de gallinas ponedoras de la provincia de Tungurahua.

Sobre la base de los datos del censo avícola disponible, se incluyeron, al azar, 47 granjas. Para realizar un muestreo homogéneo se generaron clústeres con la georreferenciación de las granjas, tomando en cuenta zonas de tolerancia de 2 km. La creación de clústeres se realizó con el software QGIS 3.18.1 (Figura 1). El muestreo se realizó entre noviembre de 2020 y agosto de 2021.

Figura 1. Clústeres para muestreo del virus de la bronquitis infecciosa aviar (IBV) en la provincia de Tungurahua. a. Mapa político del Ecuador (https://www.prolipa.com.ec/actividades/mapapolitico). b. Imagen generada en QGIS 3.28.1 a partir de las coordenadas geográficas de las granjas avícolas de la provincia de Tungurahua que se identifican con círculos. Se observan siete clústeres identificados con letras de la I a la O.

De cada granja se seleccionó un galpón con aves de 15 a 100 semanas de edad y de cada galpón se tomaron cinco aves al azar, las que fueron transportadas hasta la Facultad de Medicina Veterinaria y Zootecnia de la Universidad Central del Ecuador.

Se tomaron hisopados cloacales y laringotraqueales de todas las aves y luego se sacrificaron mediante electrocución de una sola aplicación (FCVUNLP-CICUAL 124-5-22T). Se tomaron, además, muestras de tráquea, pulmón, oviducto y tonsilas cecales, las que fueron colocadas en frascos individuales con solución tamponada de fosfato estéril con antibióticos (penicilina 10.000 U/ml, estreptomicina 10.000 µg/ml, Eurobio SA) y posteriormente almacenadas a -80 °C. Se escogieron, además, 17 vacunas a virus vivo, usadas en la zona de estudio, para comparar sus secuencias con aquellas obtenidas a partir de las muestras de campo.

Las muestras de hisopados se descongelaron de manera rápida a 37 °C en baño María, se centrifugaron a 4500 rpm durante 30 minutos a 4 °C y se tomaron 1,5 ml de cada uno de los sobrenadantes para ser conservados a -80 °C. Las muestras de órganos se descongelaron como fue descripto, se cortaron en trozos pequeños y se molieron y homogeneizaron con medio mínimo esencial estéril (MEM, MP Biomedicals) para lograr suspensiones al 20% p/v. Estas suspensiones se centrifugaron a 4500 rpm durante 30 minutos a 4 °C y 0,5 ml de cada uno de los sobrenadantes obtenidos se agruparon de acuerdo con el siguiente criterio:

a) muestras de aparato respiratorio (hisopado laringotraqueal, suspensión de tráquea y pulmón) b) muestras de aparato reproductivo y sistema inmune (suspensión de oviducto y tonsilas cecales) c) hisopado cloacal

Se realizó, además, la reconstitución de las vacunas mencionadas en la Tabla 1, con el diluyente correspondiente y según las indicaciones del fabricante.

El ARN de cada grupo de muestras (a, b y c), así como el de las vacunas, fue extraído por separado mediante el uso del kit comercial QIAmp® Viral RNA Mini Kit (QIAGEN). Para la extracción se utilizaron 140 µl del líquido sobrenadante y se siguieron los pasos mencionados en el instructivo del fabricante. A partir de los 80 µl de ARN obtenido se procedió con la retrotranscripción a ADNc utilizando para ello el OneScript Plus cADN Synthesis Kit (abm, G236). Se siguió el protocolo establecido en las indicaciones, empleando 25 ng/µl como concentración final mínima y oligonucleótidos al azar.

Tabla 1. Cepas del virus de la bronquitis infecciosa aviar IBV usadas en las vacunas que se comercializan en la zona de muestreo.

Se realizó un control interno de extracción y retrotranscripción mediante PCR para una región del exón 3 de la β-actina aviar. Se utilizaron los cebadores FW 685 (5’-GAGAAATTGTGCGTGACATCA-3’) y RV 818 (5’- CCTGAACCTCTCATTGCCA- 3’) para amplificar un segmento de 152 pb (Marandino, 2013). La prueba para β-actina se llevó a cabo con cada tipo de muestra (a, b y c).

Para amplificar un segmento de 143 pb de la región 5’UTR del IBV se utilizaron los siguientes cebadores: IBV5’ GU391 (5’- GCTTTTGAGCCTAGCGTT-3’) y IBV5’ GL533 (5’- GCCATGTTGTCACTGTCTATTG-3’) (Callison et al., 2006).

Para realizar los ensayos de PCR se utilizó la enzima polimerasa GoTaq® Flexi DNA Polymerase (PROMEGA, M9295). Se trabajó inicialmente con las condiciones de termociclado recomendadas por el fabricante de la polimerasa, en un volumen total de reacción de 15 µl.

Como control positivo se usó ADNc de IBV facilitado por la Agencia de Regulación y control Fito y Zoosanitario del Ecuador (AGROCALIDAD). También se procedió a detectar la presencia del IBV en cada una de las 17 vacunas a virus vivo.

Se consideraron como positivas las granjas con, al menos, una muestra que amplificó el segmento de la región 5’UTR. Solo una muestra positiva a IBV de cada granja fue evaluada mediante PCR para la región S1. Del mismo modo, se amplificaron las muestras de vacunas, con cebadores previamente diseñados para caracterizar cepas del IBV en Argentina y Uruguay (Marandino, 2013). El par de cebadores utilizados fue: IBV-S1 5’-ACTGAACAAAAGACAGACTT-3’ e IBV-S2 5’- CCATCTGAAAAATTRCCAGT-3’, que amplifica un fragmento de 783 pb (Marandino, 2013; Marandino et al., 2015).

Para la reacción de PCR se utilizó una desnaturalización inicial a 94 °C durante 5 minutos, seguida de 40 ciclos de: desnaturalización a 94 °C durante 1 minuto, hibridación a 53 °C durante 1 minuto, extensión a 72 °C durante 1 minuto, y una extensión final a 72 °C durante 10 minutos. Se utilizó GoTaq® Flexi DNA Polymerase (PROMEGA M9295) a diferentes concentraciones (1,25U, 1,50U y 1,75U) y dNTP Mix (PROMEGA, U1515). Se ensayaron, además, diferentes concentraciones finales de cebadores (0,2 mM y 0,4 mM). El volumen final de reacción fue de 45 µl con 6 µl de muestra, 1,25 U de GoTaq y 0,2 mM de cebadores. Se verificó la especificidad de la técnica ensayando la PCR con dos vacunas a virus vivo contra la enfermedad de Newcastle (EN), ADNc del virus de la EN cedido por AGROCALIDAD y ADN del virus de Laringotraqueitis infecciosa aviar, también entregado por AGROCALIDAD. Las reacciones de PCR se realizaron en un termociclador SIMPLI AMPTM TERMAL CYCLER (APPLIED BIOSYSTEMS).

La electroforesis de los amplicones se realizó en una cuba horizontal ENDURO GelXL (LABNET), utilizando buffer TAE (40mM, EDTA 0,5M pH 8,0) a 120V y 400mA durante 60 minutos. El resultado de la migración electroforética en geles de agarosa (PROMEGA) al 2 % se observó en un transiluminador UV (LABNET) y se registró con el uso de una cámara fotográfica digital.

Los amplicones obtenidos para el fragmento del gen S1 del IBV fueron enviados para secuenciación a Macrogen Inc. (Corea). Las secuencias de nucleótidos (nt) fueron alineadas con los programas bioinformáticos MEGA-X 11.0.8. y Geneious Prime 2021.2.2. El alineamiento de las secuencias se realizó con el algoritmo de Clustal Omega. Para el análisis filogenético se emplearon secuencias parciales de la subunidad S1 (desde el codón de inicio ATG hasta el sitio 664 pb). Se incluyeron, además de las secuencias de las muestras y las 17 secuencias vacunales logradas durante la investigación, 32 cepas prototipos de la clasificación genotípica de Valastro et al. (2016) (Tabla 2) y 32 secuencias consideradas como representativas para América, Europa y China (Tabla 3).

El árbol filogenético se construyó empleando el programa IQ-Tree (https://www.hiv.lanl.gov/content/sequence/IQTREE/iqtree.html). IQTree realizó automáticamente los cálculos y estimó el mejor modelo a través de ModelFinder (Kalyaanamoorthy et al., 2017; Nguyen et al., 2015). Además, el programa construyó un árbol consenso de máxima similitud con 1000 ultra-bootstrap para el conjunto de datos (Hoang et al., 2017).

Tabla 2. Cepas prototipo de cada linaje según clasificación genotípica del virus de la bronquitis infecciosa aviar (IBV) Tomado de Valastro et al (2016).

Tabla 3. Cepas representativas del virus de la bronquitis infecciosa aviar IBV procedentes de América Europa y China.

El control interno de extracción y retrotranscripción funcionó correctamente y todas las muestras evaluadas resultaron positivas para ß-actina aviar (Figura 2). Por otro lado, 16 granjas (34,04%) resultaron positivas a IBV al registrar, al menos, una de las tres muestras empleadas que amplificara para la región 5’UTR (Tabla 4) (Figura 2). Del mismo modo, las 16 granjas resultaron positivas a la amplificación del segmento S1 (Tabla 4) (Figura 3).

Tabla 4. Resultados de pruebas de PCR para la región 5’UTR y S1 del virus de la bronquitis infecciosa aviar (IBV).

Figura 2. Técnica de RT PCR realizadas con muestras de cuatro granjas de la provincia de Tungurahua, Ecuador. Electroforesis en gel de agarosa al 2% teñidas con SYBR Safe DNA gel stain (INVITROGEN). a. Amplificación de un segmento de 152 pb de la β-actina aviar. b. Amplificación de un segmento de 143 pb de la región 5’UTR del virus de la bronquitis infecciosa aviar (IBV). Referencias: M: marcador de pares de bases, C+ control positivo, C- control negativo.

Con el servidor en línea gratuito IQ-TREE el mejor modelo que se obtuvo para el conjunto de secuencias fue: TPM3u+F+I+G4. A partir de este árbol consenso generado con IQ-TREE (Figura 4) se pudo determinar que cinco (U5250, U5261, U5272, U5281, U5284) de las 16 secuencias obtenidas a partir de las muestras se ubicaron en el mismo clado donde estaban presentes las cepas de tipo Massachusetts del linaje GI-1 (Figura 4). Por otro lado, siete de las secuencias recuperadas (U5240, U5241, U5263, U5264, U5268, U5271, U5280) se agruparon filogenéticamente en el clado del linaje GI-13 de las cepas tipo 793B. Finalmente, cuatro secuencias (U5247, U5256, U5262, U5270) de las granjas estudiadas se agruparon en un clado donde se ubicaron tanto cepas asiáticas como latinoamericanas del tipo Q1 del linaje G1-16 (Figura 4). Con relación a las vacunas comerciales analizadas, se observó que 15 estaban agrupadas en el linaje GI-1, en el que se ubican las cepas tipo Massachusetts, y solo dos vacunas (2B y 4B) se situaron en el mismo clado de las cepas vacunales de tipo 793B (Figura 4).

Figura 3. Gel de agarosa al 2% teñido con SYBR Safe DNA gel stain (INVITROGEN). Se observa la amplificación de un segmento de 783 pb del gen codificante de la proteína S1 del virus de la bronquitis infecciosa aviar (IBV). La muestra analizada correspondió a la codificada como U5280 y que resultó positiva a la región 5’UTR. Referencia. M: marcador de pares de bases, C+ control positivo, C- control negativo.

El análisis filogenético realizado mostró que cinco de las cepas encontradas en la provincia de Tungurahua correspondían a cepas vacunales de tipo Massachusetts. Es importante mencionar que las cepas vacunales de tipo Massachusetts son las recomendadas para utilizar en Ecuador (Organización Mundial de Sanidad Animal -OIE-, 2018). En un análisis filogenético realizado con aislamientos obtenidos en Uruguay entre 2009 y 2012 y en Argentina entre 2009 y 2011 se concluyó que las cepas no agrupaban con la cepa vacunal Massachusetts, la que también era la única autorizada para vacunación en estos países durante el desarrollo del trabajo (Marandino, 2013). Se debe mencionar que las muestras con las que trabajó esta investigadora provenían de brotes de infección por IBV y las muestras del presente estudio procedieron de aves que no manifestaron signología relacionada con IBV, pero que, al igual que en Argentina y Uruguay, habían sido vacunadas con la cepa Massachussets. Por su parte, Callison et al. (2006), también usando muestras de brotes de bronquitis infecciosa para amplificar la subunidad S1, concluyeron que los aislados estaban altamente relacionados con la cepa Massachusetts y que existía la probabilidad de detección y aislamiento de cepas vacunales. Del mismo modo, estudios de epizootiología molecular realizados en Suecia revelaron una fuerte probabilidad de que la vacunación con la cepa Massachussets había llevado a la propagación del virus vacunal, inclusive hacia una zona libre de vacunación (Farsang et al., 2002). En consecuencia, resulta razonable considerar que el virus circulante identificado en la provincia de Tungurahua sea el virus vacunal o cepas homólogas.

Siete de las secuencias obtenidas se encontraron en el clado del linaje GI13 (las cepas tipo 793B) (Valastro et al., 2016). Según Valastro et al. (2016) la cepa Morroccan-G/83, aislada en Marruecos en 1983, sería la cepa progenitora del linaje, a pesar de que la primera cepa conocida fue aislada en Europa. Actualmente, las cepas de este linaje se encuentran presentes en muchas partes del mundo (Picault et al., 1995). En particular en América del Sur, las cepas de este tipo fueron introducidas como vacunas con la cepa 4/91 (Eldemery et al., 2017). Las cepas analizadas en el presente estudio podrían tener un origen vacunal, dado que dos de las vacunas analizadas (2B y 4B) se encuentran en este clado. Aunque la OIE (2018) sugiere para América del Sur la utilización de la cepa vacunal Massachusetts, en Ecuador están registradas para uso otras cepas vacunales como Arkansas, 4/91 y QX (Agencia Ecuatoriana de Aseguramiento de la Calidad del Agro -AGROCALIDAD-, 2021). De esta manera, se podría explicar el origen de las cepas variantes tipo 793 o 4/91 encontradas en el estudio desarrollado, tomando en cuenta que las cepas emergentes podrían ser el producto de la recombinación entre cepas de campo y cepas utilizadas para la vacunación (Quinteros et al., 2016).

Por otro lado, se estableció que cuatro de las muestras obtenidas en la provincia de Tungurahua estaban asociadas con cepas tipo Q1. Estas cepas tipo Q1 fueron aisladas inicialmente en Asia (Liu et al., 2009b; Yu et al., 2001) y luego se pudo evidenciar su presencia en América del Sur (de Wit et al., 2017; Marandino et al., 2015). Por ejemplo, en el análisis filogenético basado en la región S1 y realizado a partir de aislados chilenos de los años 2008 y 2009 se encontraron cepas de genotipo Q1 (de Wit et al., 2017). En consecuencia, la introducción de la cepa asiática podría haber tomado una ruta indirecta desde mercados latinoamericanos, dado que Ecuador no tiene comercio de importación directo con Asia. Además, debe recordarse que nuevas cepas de IBV pueden surgir también como producto del proceso replicativo característico de los coronavirus (Jackwood, 2012; de Wit et al., 2011).

Este trabajo constituye el primer estudio de caracterización molecular y análisis filogenético realizado con cepas encontradas en gallinas de postura comercial de Ecuador. Se concluyó que hay tres tipos de cepas circulando en la provincia de Tungurahua, dos de tipo vacunal o sus homólogas y una de origen desconocido. Se requieren estudios moleculares en los que se incluya la amplificación del gen S completo, a fin de que, con base en los estudios filogenéticos actuales que incluyen toda esta proteína, se puedan obtener datos más concluyentes.

Figura 4. Árbol filogenético consenso con modelo TPM3u+F+I+G4. Los colores de letras hacen referencia a la procedencia de las cepas: rojo (América), amarillo (Asia), verde (Europa), celeste (Oceanía) y negro (África). Las cepas muestrales se observan en color negro y destacadas en recuadros rosa y las cepas de vacunas comerciales utilizadas en la zona de muestreo se muestran en color negro y ubicadas en recuadros lila. Los números a lo largo de las ramas hacen referencia a los valores de ultrabootstrap.

El presente trabajo se financió con fondos otorgados como parte del programa de becas doctorales de la Universidad Central del Ecuador. Agradecemos a los avicultores de la provincia del Tungurahua, al Instituto Ecuatoriano del Huevo (INEH), a los médicos veterinarios Juan Valle y Richard Tuapanta, a los técnicos José Luis Medina y Sofía De Janón de la Unidad de Investigación de Enfermedades Transmitidas por Alimentos y Resistencia a los Antimicrobianos (UNIETAR) de Facultad de Medicina Veterinaria y Zootecnia de la Universidad Central Ecuador. También se agradece al personal del Laboratorio de Virología (LAVIR) de la Facultad de Ciencias Veterinarias de la Universidad Nacional de La Plata. El aporte de todos ellos fue de gran utilidad para llevar a cabo el estudio.

Declaración de conflicto de intereses: Todos los autores declaran que no existe conflicto de intereses, incluyendo las relaciones financieras, personales o de otro tipo que pudieran influir en el trabajo.