text in 
English (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Similars in
SciELO 
Permalink
Introducción
Los sistemas de producción de leche a pequeña escala mejoran la seguridad alimentaria y proporcionan ingresos económicos en zonas rurales a nivel mundial1. En México, este sistema de producción posee aproximadamente el 23 % del inventario ganadero2, contribuye con el 30 % de la producción nacional de leche3 y representa el 73 % de los establos lecheros4. Las unidades de producción en este sistema se caracterizan por incorporar mano de obra familiar, utilizan en mayor medida razas especializadas en la producción de leche, tienen pocos vientres en producción y niveles medios-bajos de tecnificación5-7. La mejora productiva de estos establos lecheros a pequeña escala contribuye a disminuir la pobreza en zonas rurales(1), y promover el desarrollo de las comunidades8,9.
La rentabilidad de los establos lecheros está fuertemente asociada con el desempeño reproductivo eficiente10,11. La identificación de factores de riesgo que comprometen el desempeño reproductivo es primordial para diseñar e implementar estrategias que mejoren la productividad. Existen estudios en hatos de producción de leche a pequeña escala que indican que la IA puede afectar el intervalo parto a primer servicio postparto y la tasa de concepción por servicio en comparación con la monta directa12,13. Además, en sistemas intensivos de producción se ha observado que el tamaño de hato influye sobre el desempeño reproductivo14-16, al igual que el nivel de seroprevalencia de enfermedades infecciosas reproductivas como neosporosis, rinotraqueitis infecciosa bovina (IBR) y diarrea viral bovina (BVD)17,18. Por lo tanto, el objetivo del presente trabajo fue determinar el impacto, como factores potenciales de riesgo a nivel de establo, del uso de IA, tamaño del hato y prevalencia de enfermedades infecciosas reproductivas sobre el desempeño reproductivo en establos productores de leche a pequeña escala en México. La hipótesis de trabajo fue que estos factores están asociados con el desempeño reproductivo en vacas lecheras.
Material y métodos
Selección de establos y registro de la información
Se llevó a cabo un estudio de cohorte prospectivo observacional (959 registros) en seis estados de México donde el sistema de producción de leche a pequeña escala tiene presencia importante. Se incluyeron al estudio 52 establos con la siguiente distribución por estado: Jalisco (23), Estado de México (10), Tlaxcala (9), Guanajuato (4), Puebla (3), y Querétaro (3). Los criterios de selección utilizados fueron los siguientes: soporte primario de mano de obra familiar en la unidad de producción, que contaran con un número de vientres en producción entre 10 y 100, que la producción de leche fuera el objetivo primario del establo y que el nivel de incorporación de tecnología fuera medio-bajo. El porcentaje de vacas Holstein y el número promedio de vacas en los establos fueron 91.3 % y 30.3 ± 2.4, respectivamente. La tasa estimada de desecho y la producción de leche por vaca fueron 26.4 % y 17.10 ± 0.5 kg/día, respectivamente. Los establos incorporados al estudio cumplen características de establos de producción a pequeña escala en México descritos en otras publicaciones2,7,19. El periodo de captura de información a nivel de campo fue durante 18 meses, periodo en el que se registraron los siguientes eventos reproductivos: fechas de parto, fechas y tipo de servicios (inseminación artificial o monta directa), ocurrencia de asistencias al parto o retenciones de placenta, y resultado de diagnóstico de gestación alrededor de 50 días post-servicio.
Eventos de interés y clasificación de factores potenciales de riesgo relacionados con el desempeño reproductivo
Los eventos de interés considerados fueron: asistencia al parto (asistencia menor y asistencia mayor en la misma categoría), retención de placenta (>12 h), días a primer servicio mayor a 70 días en leche (DPS>70), días abiertos mayores a 110 días en leche (DA>110) y vacas no gestantes al primer servicio (NG1S). Con base a lo que se ha sugerido para el sistema de producción de leche a pequeña escala se establecieron los valores límite >70 DPS y >110 DA como indicadores de falla reproductiva(19).
Los factores potenciales de riesgo incluyeron al uso de la IA, hatos grandes y seroprevalencias altas de neosporosis, IBR y BVD. Para clasificar los establos por tipo de servicio, se consideraron en la categoría IA a aquellas unidades de producción donde al menos 75 % de los servicios fueron brindados utilizando esta tecnología y en la categoría de monta natural cuando al menos el 75 % de los servicios se proporcionó por monta natural. Los límites para tamaño del hato y seroprevalencias se establecieron de acuerdo con la distribución por cuartiles en la muestra de estudio20. La clasificación del tamaño de hato se hizo con base en el número de vientres en producción promedio por establo durante el periodo
de captura de información a nivel de campo; el tercer cuartil correspondió a 33 vacas (clasificación <33 o ≥33). El tercer cuartil se estableció como límite para clasificar a los establos con seroprevalencia alta; neosporosis (≥84 %), IBR (≥38 %) y BVD (=100 %) (Cuadro 1).
Cuadro 1 Seroprevalencias para neosporosis, rinotraqueítis infecciosa bovina (IBR) y diarrea viral bovina (BVD) en 52 establos de producción de leche a pequeña escala
| Media ± EE | Mínimo | Cuartil 1 | Mediana | Cuartil 3 |
Máximo | |
|---|---|---|---|---|---|---|
| Neosporosis | 52.7±4.5 | 0 | 33 | 50 | 84 | 100 |
| IBR | 23.3±1.8 | 0 | 0 | 23.5 | 38 | 75 |
| BVD | 59.7±3.2 | 0 | 28 | 75 | 100 | 100 |
Identificación de animales seropositivos a neosporosis, IBR y BVD
Se colectaron muestras de sangre por punción de vena coccígea (sistema vacutainer), de forma aleatoria en el 10 % de las vacas en producción en cada hato en estudio. Las muestras de sangre fueron mantenidas a 4 °C por 24 h, posteriormente fueron centrifugadas (2,500 xg por 10 min a 4 °C) para separar el suero, el cual fue congelado a -20° C hasta su análisis.
La detección de anticuerpos contra Neospora caninum se realizó mediante la prueba de ELISA utilizando un kit comercial (Laboratorios IDEXX), de acuerdo con las instrucciones del fabricante. El análisis de los sueros para DVB se realizó con un kit comercial de ELISA por bloqueo (CIVTEST bovis BVD/Bd P80, Laboratorios Hipra), siguiendo las instrucciones del fabricante. El análisis de IBR se realizó mediante la técnica de seroneutralización en placa, empleando la línea celular MDBK (células de riñón de bovino) y el virus de Referencia IBR758, con título de 105.6 TCID 50% a una dilución entre 500-1000 dosis infectante/ml. Los sueros se diluyeron de 1:2 hasta 1:128 y se observó el efecto citopático producido por el virus para determinar la muestra como positiva(18). No se contó con historial de vacunación en cada unidad de producción, sin embargo, en las regiones estudiadas esta práctica es común(13,18).
Análisis estadístico
Todos los análisis fueron llevados a cabo utilizando el paquete estadístico SAS 9.3 (SAS Institute Inc., Cary, NC). Para identificar a los factores de riesgo se utilizaron análisis de regresión logística múltiple (PROC LOGISTIC) siguiendo la metodología implementada por Potter et al21. Para elaborar estos modelos, el primer paso fue realizar pruebas de regresión logística simple entre los eventos de interés y los factores potenciales de riesgo. Los factores
con un valor de P<0.35 fueron retenidos y analizados posteriormente para colinealidad21. Para prevenir colinealidad en los modelos múltiples, se obtuvieron coeficientes de correlación y se aplicaron pruebas de χ2 en pruebas pareadas de factores retenidos utilizando el procedimiento FREQ opción CHISQ. Cuando en un par de factores, el límite de confianza del coeficiente de correlación no incluyó al 0 y el valor de P de χ2 fue <0.05, ambas variables no fueron incluidas en el mismo modelo múltiple. Finalmente, con el objetivo de obtener modelos múltiples parsimoniosos, se utilizó la opción BACKWARD para retener variables significativas a un valor de P<0.121. Los modelos múltiples finales incluyeron únicamente a los efectos principales y la razón de momios (OR) fue utilizada como una medida de asociación entre los factores de riesgo y las variables de interés.
Resultados
Eventos de interés y factores potenciales de riesgo a nivel de establo relacionados con el desempeño reproductivo
Estadísticas descriptivas para eventos de interés. Las prevalencias de los eventos de interés asistencia al parto y retención de placenta fueron 13.2 y 11.7 % respectivamente. La proporción de vacas con DPS>70, DA>110 y NG1S fueron 64.9, 46.4 y 50.5 % respectivamente.
Factores potenciales de riesgo. Los establos con el factor potencial de riesgo IA representaron el 73.9 %, mientras que los clasificados con tamaño ≥33 vacas representaron el 41.3 %. En el Cuadro 1 se muestra la distribución de seroprevalencias para neosporosis, IBR y BVD.
Factores de riesgo asociados a la falla en el desempeño reproductivo
El nivel de significancia y las razones de momios obtenidos en los análisis de regresión logística simple se muestran en el Cuadro 2. En el Cuadro 3 se muestran los modelos múltiples para cada evento de interés. Para la variable de interés asistencia al parto, los factores de riesgo identificados (P<0.10) fueron establos con ≥33 vacas y seroprevalencia alta de neosporosis (Cuadro 4). Con respecto a la variable de interés retención de placenta, no se detectaron factores de riesgo (OR>1). Sin embargo, IA y seroprevalencia alta de IBR fueron factores significativos (OR<1, P<0.10; Cuadro 4). Para la variable de interés DPS>70 los factores de riesgo identificados fueron IA y seroprevalencias altas de IBR y BVD (P<0.10; Cuadro 5). Con relación a la variable de interés DA>110, la IA fue el único factor de riesgo identificado (P<0.10; Cuadro 5). Para la variable de interés NG1S, establos ≥33 vacas fue el único factor de riesgo identificado (P<0.10).
Cuadro 2 Valor de probabilidad (P) y razón de momios (OR) para factores potenciales de riesgo considerando diferentes eventos de interés; análisis de regresión logística simple
| Eventos de interés | |||||
|---|---|---|---|---|---|
| Factores | AP (P;OR) | RP (P;OR) | DPS>70 (P;OR) |
DA>110 (P;OR) |
NG1S (P;OR) |
| IA | 0.048; 0.67 | 0.072; 0.68 | <0.001; 2.39 | 0.054; 1.32 | 0.162; 0.82 |
| Establos ≥33 vacas | 0.090; 1.52 | 0.821; NC | 0.853; NC | 0.823; NC | 0.001; 1.69 |
| Neosporosis alta | <0.001; 2.29 | 0.099; 0.64 | 0.414; NC | 0.357; NC | 0.718; NC |
| IBR alta | 0.202; 1.32 | 0.005; 0.42 | 0.169; 1.25 | 0.484; NC | 0.916; NC |
| BVD alta | 0.010; 0.46 | 0.793; NC | <0.001; 1.86 | 0.553; NC | 0.049; 0.72 |
IA= inseminación artificial; AP=asistencia ak parto; RP0retención placentaria; DPS>70= días a primer servicio mayor a 70 días en leche; DA>110= días abiertos mayores a 110 días en leche; NG1S= vacas no gestantes al primer servicio; NC= indica OR no calculado por no ser significativo.
Cuadro 3 Factores potenciales de riesgo no colineales para eventos de interés incluidos en modelos múltiples
| Evento de interés | Modelo | Factores potenciales de riesgo |
|---|---|---|
| Asistencia al parto | 1 | IA + IBR Alta |
| 2 | Establos ≥33 vacas + Neosporosis alta | |
| 3 | BVD alta | |
| Retención placentaria | 1 | IA + Neosporosis alta |
| 2 | IA + IBR alta | |
| DPS>70 días | 1 | IA + IBR alta |
| 2 | BVD alta + IBR alta | |
| DA>110 | 1 | IA |
| NG1S | 1 | Establos ≥33 vacas |
| 2 | BVD alta | |
| 3 | IA |
DPS>70= días a primer servicio mayor a 70 días en leche; DA>110= días abiertos mayores a 110
días en leche; NG1S= vacas no gestantes al primer servicio; IA= inseminación artificial; BVD= diarrea viral bovina; IBR= rinotraqueitis infecciosa bovina.
Cuadro 4 Efecto de variables de estudio sobre asistencia al parto y retención de placenta en modelos múltiples
| Variable de interés | Efectos | OR | IC 95% | P |
|---|---|---|---|---|
| Asistencia al parto Modelo 1 | Tipo de servicio: MN | Ref. | N/A | N/A |
| MNTipo de servicio: IA | 0.67 | 0.67 0.45-0.99 | 0.048 | |
| Modelo 2 | Tamaño del hato: <33 vacas | Ref. | N/A | N/A |
| Tamaño del hato: ≥33 vacas | 1.51 | 0.90-2.45 | 0.090 | |
| Neosporosis: Resto | Ref. | N/A | N/A | |
| Neosporosis: Alta | 2.28 | 1.52-3.40 | <0.001 | |
| Modelo 3 | BVD: Resto | Ref. | N/A | N/A |
| BVD: Alta | 0.46 | 0.26-0.83 | 0.010 | |
| Retención placentaria | ||||
| Modelo 1 | Tipo de servicio: MN | Ref. | N/A | N/A |
| Tipo de servicio: IA | 0.68 | 0.45-1.04 | 0.072 | |
| Modelo 2 | Tipo de servicio: MN | Ref. | N/A | N/A |
| Tipo de servicio: IA | 0.66 | 0.43-1.01 | 0.055 | |
| IBR: Resto | Ref. | N/A | N/A | |
| IBR: Alta | 0.41 | 0.23-0.75 | 0.004 |
P= valor de probabilidad; MN= monta natural; IA= inseminación artificial; BVD= diarrea viral bovina; IBR= rinotraqueítis infecciosa bovina; IC= intervalo de confianza de razón de momios; OR= razón de momios; N/A= no aplicable.
Cuadro 5 Efecto de variables de estudio sobre días a primer servicio (DPS>70), días abiertos (DA>110) y no gestante al primer servicio (NG1S) con diferentes modelos múltiples
| Variable de interés | Efectos | OR | IC 95% | P |
|---|---|---|---|---|
| DPS>70 Modelo 1 | Tipo de servicio: MN Tipo de servicio: IA IBR: Resto IBR: Alta Resto | Ref. 2.42 Ref. 1.32 Ref. | N/A 1.81-3.2 N/A 0.95-1.84 N/A | N/A <0.001 N/A 0.097 N/A |
| Modelo 2 | BVD: Resto | Ref. | N/A | N/A |
| BVD: Alta | 1.86 | 1.31-2.64 | <0.001 | |
| DA>110 | ||||
| Modelo 1 | Tipoo de servicio: MN | Ref. | N/A | N/A |
| Tipo de servicio: IA | 1.32 | 0.99-1.8 | 0.054 | |
| NG1S | ||||
| Modelo 1 | Tamaño del hato:<33 vacas | Ref. | N/A | N/A |
| Tamaño del Hato: ≥33 vacas | 1.69 | 1.23-2.32 | 0.001 | |
| Modelo 2 | BVD: Resto | Ref. | N/A | N/A |
| BVD: Alta | 0.72 | 0.25-1.00 | 0.049 |
P= valor de probabilidad; IC= intervalo de confianza de razón de momios; MN= monta natural; IA= inseminación artificial; BVD= diarrea viral bovina; IBR= rinotraqueitis infecciosa bovina; OR= razón de momios; N/A= no aplicable.
Discusión
El sistema de producción de leche a pequeña escala en México es altamente heterogéneo en cuanto a manejo productivo, reproductivo y estatus sanitario13,19. Desde un punto de vista sanitario, neosporosis, IBR y BVD son enfermedades asociadas a desórdenes reproductivos22,23. En el presente estudio, la media en seroprevalencias para neosporosis por establo representó valores similares a lo reportado previamente para el sistema de producción de leche a pequeña escala (51.7 %)24. Sin embargo, es mayor a lo reportado para el sistema intensivo en México (~43 %)(18,25). Es probable que las medidas de bioseguridad sean menos estrictas en los establos de producción de leche a pequeña escala, lo cual pudiera estar incrementando factores de riesgo asociados a la presencia de Neospora, como la presencia de perros en las unidades de producción26,27. Por otra parte, la seroprevalencia de IBR en el presente estudio fue similar a lo reportado previamente en el sistema de producción de leche a pequeña escala en México con tasas que rondan el 22 %24,28. Sin embargo, algunos estudios han reportado hasta el 69 % de seroprevalencia en el sistema de producción a pequeña escala18. Referente a BVD, la seroprevalencia promedio obtenida por establo se encuentra en un punto medio entre lo reportado para este sistema de producción (52-81%)18,24. Anteriormente se ha propuesto que estas enfermedades son endémicas y con alta
prevalencia en ganado lechero en los sistemas intensivo, doble propósito y a pequeña escala en México18,29,30. A pesar de que las altas seroprevalencias encontradas pudieran deberse a anticuerpos vacunales para IBR y BVD, estos resultados confirman la importancia de estas enfermedades reproductivas en establos de producción a pequeña escala debido a su diseminación sobre una gran área o distancia, y alta prevalencia de estas enfermedades.
Hasta donde tenemos conocimiento, este es el primer trabajo que reporta factores de riesgo a nivel de establo asociados con el desempeño reproductivo en el sistema de producción de leche a pequeña escala en México. Los establos con seroprevalencia alta de neosporosis y establos con 33 o más vacas tuvieron 128 y 51 % respectivamente, más probabilidad de requerir asistencia al parto. Estos resultados son consistentes con reportes previos sobre una asociación significativa entre asistencia al parto y animales seropositivos a Neospora31. Sin embargo, otros estudios no han encontrado esta asociación32, por lo tanto, son necesarios más estudios para clarificar esta asociación potencial. La asistencia al parto fue más prevalente en establos ≥33 vacas. Una posible explicación es que en establos más grandes utilicen sementales (IA o monta natural) de mayor talla, o bien que existan problemas asociados a la condición corporal al parto, hipótesis pendiente de ser desafiada. Otra explicación radica en que, en este sistema de producción, el productor y su familia se hacen cargo personalmente de toda la operación del establo y del campo6. Por otra parte, se ha reportado que el manejo en el periparto de la vaca y la atención obstétrica correcta son algunos de los principales factores a considerar para controlar los problemas de distocia33. Es razonable pensar que, debido a las múltiples tareas del productor, él y los miembros de su familia son incapaces de proveer adecuado cuidado durante el parto en establos ≥33 vacas.
Por otra parte, los animales en los establos con mayor prevalencia de BVD y los establos que utilizan IA tienen un menor riesgo de asistencia al parto. Este resultado era esperado, en otros estudios se ha reportado que el peso al nacimiento de becerros positivos a antígenos de BVD son 7 kg más ligeros que su contraparte de becerros negativos34. Referente al impacto de la IA sobre la ocurrencia de asistencia al parto, actualmente existen en el mercado una gran variedad de sementales con diferentes características de conformación, incluyendo la facilidad de parto35,36. Estos resultados sugieren que los productores en las regiones en estudio pudieran estar seleccionando sementales con base a la facilidad de parto (comunicación personal, MC. Fernando Villaseñor). Aunado a esto, los animales en establos que utilizan la IA tienen un menor riesgo de sufrir Retención de Placenta (Cuadro 4). La asistencia al parto es uno de los factores de riesgo a nivel de individuo más importantes asociados con retención de placenta37. Por lo tanto, es razonable pensar que con la IA se reduzca la prevalencia de asistencias al parto y en consecuencia las retenciones de placenta. Por otra parte, las vacas en los establos con seroprevalencia alta de IBR tienen un menor riesgo de sufrir retención de placenta. Este resultado pudiera parecer contradictorio porque IBR ha sido asociado a aborto y retención de placenta38. Sin embargo, una posible explicación radica en que las vacas incluidas en nuestro estudio tengan anticuerpos de origen vacunal(18); reduciendo la incidencia de abortos por IBR y en consecuencia disminuyendo el riesgo de sufrir retención de placenta 39.
Las vacas en los establos que utilizan IA tuvieron 142 % más probabilidad de tener DPS>70. Desde un punto de vista reproductivo, el éxito al implementar esta tecnología depende de la eficiencia en la detección de estros40. Aunque en el presente estudio no fue determinado, en los establos de producción de leche a pequeña escala en México(12) y en establos intensivos es común que exista una deficiente tasa de detección de estros41. Las vacas en establos con seroprevalencias altas de IBR y BVD tuvieron 32 y 86 % más probabilidad de tener DPS>70; pudiendo ir en IBR desde una reducción en la probabilidad del 5 % hasta un incremento del 80 %. Las consecuencias reproductivas de la infección con BVD han sido documentadas42,43. De acuerdo con nuestro conocimiento, no existen reportes que hayan asociado seroprevalencias altas de estas enfermedades con días a primer servicio. Además, desde un punto de vista patológico, no está claro cómo estas enfermedades pudieran incrementar DPS>70. Una posible explicación radica en que la seroprevalencia alta de estas enfermedades pudiera impactar de forma indirecta a este indicador44. Por ejemplo, seroprevalencias altas de IBR y BVD pudieran estar correlacionadas con otros factores de riesgo como alto grado de estabulación o mal manejo de desechos biológicos45.
Por otra parte, las vacas en establos que brindan IA tuvieron 32 % más probabilidad de tener DA>110; este cambio puede ir desde una reducción en la probabilidad de 1 % hasta incremento del 80 %. Desde un punto de vista de desempeño reproductivo, esta tecnología para el mejoramiento genético también repercute negativamente en indicadores como días a primer servicio, sugiriendo la existencia de deficiencias en la detección de estros12,46. Para contrarrestar esto, en los establos intensivos de producción de leche se han implementado protocolos de sincronización de estros47. Considerando la buena fertilidad a primer servicio en este sistema (49.5 %), la implementación de protocolos de inseminación a tiempo fijo adaptados específicamente para este sistema de producción pudiera ser una estrategia para mejorar el desempeño reproductivo48.
Las vacas en establos con 33 o más vacas tuvieron 69 % más probabilidad para NG1S. En otros sistemas de producción, se ha reportado que conforme los establos incrementan su tamaño, la capacidad para manejar reproductivamente al hato disminuye14-16. En el sistema de producción de leche a pequeña escala, aunque son de menor tamaño comparados con los establos intensivos, el efecto de tamaño de hato es también evidente para este indicador. Por otra parte, se ha reportado que la BVD provoca muertes embrionarias tempranas y subfertilidad en ganado lechero(49-51). Los resultados indican que hatos con seroprevalencia alta de BVD tienen menor riesgo de NG1S, un resultado que pudiera parecer contradictorio. Sin embargo, una posible explicación es que la alta seroprevalencia en estos establos se deba a anticuerpos vacunales18.
La IA es un factor significativo común para la mayoría de los eventos de interés. El uso de esta tecnología para el mejoramiento genético reduce la prevalencia de problemas posteriores
al parto sin afectar la fertilidad a primer servicio, sin embargo, incrementa los días a primer servicio y los días abiertos (probablemente por deficiencias en la detección de estros). Una posible estrategia para aprovechar las ventajas de la IA pudiera ser la implementación de protocolos de sincronización a tiempo fijo para el primer servicio con semen sexado considerando la buena fertilidad observada en los hatos pertenecientes a este sistema de producción. Sin embargo, son necesarios estudios de factibilidad financiera para poder realizar recomendaciones a gran escala en el sector.
Finalmente, una limitación del presente estudio es la incertidumbre acerca del origen de los anticuerpos en las pruebas serológicas. En la base de datos no se contó con antecedentes de vacunación confiables por la falta de disciplina que existe en este tipo de establos para llevar registros. A pesar de que en las regiones estudiadas es común la vacunación13,18, esta incertidumbre limita la posibilidad de hacer inferencias más precisas con base en los resultados de este estudio. Sin embargo, se considera que la interpretación de los resultados es conservadora y aún fundacional para estudios posteriores en términos epidemiológicos y patológicos.
Conclusiones e implicaciones
En conclusión, la IA, el tamaño del hato y las seroprevalencias altas de neosporosis, IBR y BVD son factores asociados al desempeño reproductivo en establos pertenecientes al sistema de producción de leche a pequeña escala en México. Los factores de riesgo identificados para asistencia al parto fueron hatos con 33 o más vacas y neosporosis alta; para DPS>70 fueron IA y seroprevalencias altas de IBR y BVD; para DA>110 fue IA; para NG1S fueron hatos con 33 o más vacas; para retención de placenta no se identificaron factores de riesgo. Además, el presente estudio es un recordatorio sobre la vital importancia de prevenir neosporosis, IBR y BVD debido a que estas enfermedades están ampliamente distribuidas en establos lecheros en el país.
