Efecto de la raza del macho sobre la densidad de empadre e indicadores reproductivos de cuyes (Cavia porcellus) mejoradas de Chota-Cajamarca

INTRODUCCIÓN

El cuy (Cavia porcellus) constituye una especie pecuaria de gran importancia socioeconómica en la región andina, particularmente en el Perú, destacando la explotación de la raza Perú por su alta prolificidad, precocidad y eficiencia en la conversión de forraje en proteína animal de alta calidad (Chauca, 2023), características que influyen en las técnicas de manejo aplicado por los productores a otras razas esperando encontrar los mismos o mejores resultados ignorando que cada una tiene una potencialidad diferente como la raza Kuri, liberada por el INIA en 2021, caracterizada por tener mayor cantidad de crías por parto con respecto a la Perú y mayor porcentaje de celo post parto en las hembras, lo cual favorece la productividad por jaula o poza durante la vida útil de los reproductores.

La relación macho-hembra o densidad de empadre es un factor de manejo determinante del desempeño productivo y de acuerdo y comúnmente se maneja la relación 1:7 para raza Perú, la cual es la más

utilizada en la provincia de Chota, región Cajamarca, lo cual limita el refrescamiento de sangre en los productores y buscan otras razas como alternativa a esta situación.

En este contexto, la raza Kuri se presenta como una alternativa para empadrarse con hembras mejoradas con raza Perú, porque existen antecedentes que demuestran un incremento en el tamaño de camada de hembras de la región Lambayeque (Bustamante, 2015) y La Libertad (Carrasco, 2019).

Ante esta situación, el presente estudio se justifica por la necesidad de evaluar alternativas de manejo específicos para las razas de cuyes, buscando contribuir al conocimiento técnico sobre estrategias de manejo reproductivo eficiente de las razas disponibles para la producción de cuyes.

Por tanto, el objetivo del estudio fue evaluar el efecto de la raza sobre la densidad de empadre e indicadores reproductivos de cuyes mejoradas de Chota-Cajamarca.

MATERIALES Y MÉTODOS

Zona de estudio

El estudio se realizó entre mayo de 2022 y mayo de 2024 en el Caserío de Llasavilca Centro, provincia y distrito de Chota, Región Cajamarca, Perú. La zona se encuentra ubicada a una altitud de 2 270 m.s.n.m., registrando una temperatura media anual que oscila entre los 10 °C como mínima y 25 °C como máxima, con una precipitación fluvial promedio de 916 mm (Senamhi, 2024).

Manejo experimental

Los animales se manejaron bajo un sistema de empadre continuo en jaulas. Se aplicó una estrategia de alimentación mixta basada en el suministro de alfalfa fresca y alimento balanceado comercial peletizado con un valor nutricional de 18 % de proteína cruda (PC), 2,9 Mcal/Kg de energía digestible (ED) y un mínimo de 6 % de fibra cruda (FC), disponiendo de agua ad libitum.

Se emplearon 42 hembras jóvenes mejoradas con base genética de la raza Perú procedentes de la provincia de Chota, las cuales registraron un peso promedio inicial de 771,07 ± 8,03 g. Las hembras fueron distribuidas de forma aleatoria para ser empadradas con seis machos reproductores: tres de la raza Kuri provenientes del INIA Lima (peso promedio: 1 138,67 ± 15,28 g) y tres de la raza Perú procedentes del INIA Cajamarca (peso promedio: 1 215,33 ± 41,63 g). La evaluación del peso inicial de las hembras demostró normalidad y homogeneidad de varianzas (p > 0,05) sin diferencias estadísticas entre grupos, mientras que el análisis del peso de los machos mediante la prueba de Welch tampoco determinó diferencias significativas (p = 0,072).

Diseño experimental y tratamientos

Se utilizó un Diseño Completamente al Azar (DCA) con arreglo factorial 2 × 3, evaluándose el efecto de dos factores principales: Raza del macho (Kuri y Perú) y Densidad de empadre (6, 7 y 8 hembras por macho), lo que generó un total de seis tratamientos independientes evaluados de forma longitudinal durante dos partos consecutivos (considerados como repeticiones en el tiempo). Debido a la naturaleza intrínseca de las densidades por poza, el diseño contempló un número desigual de unidades experimentales (repeticiones) por tratamiento. Los tratamientos se estructuraron de la siguiente manera:

• T1: Macho Kuri × densidad de 6 hembras mejoradas.

• T2: Macho Kuri × densidad de 7 hembras mejoradas.

• T3: Macho Kuri × densidad de 8 hembras mejoradas.

• T4: Macho Perú × densidad de 6 hembras mejoradas.

• T5: Macho Perú × densidad de 7 hembras mejoradas.

• T6: Macho Perú × densidad de 8 hembras mejoradas.

Variables evaluadas

Las variables de respuesta zootécnica registradas durante el primer y segundo parto correspondieron a:

• Tasa de fertilidad (TF): Proporción de hembras preñadas y paridas con respecto al total de hembras expuestas al macho (%).

• Tamaño de camada al nacimiento (TCn): Número total de crías nacidas por parto.

• Peso al nacimiento por sexo (PN): Masa corporal individual de los gazapos regis-

trada en gramos dentro de las primeras 24 horas post-parto, clasificada por sexo (machos y hembras).

• Frecuencia de ocurrencia del TCn (FOTCn): Distribución porcentual y probabilística del número de crías obtenidas por evento de parto.

Modelo estadístico y análisis de datos

Para contrastar el efecto de los tratamientos sobre las variables continuas, se planteó la hipótesis nula (H0 :μ1 =μ2 =μ3 =μ4 =μ5 =μ6) de igualdad de medias frente a la hipótesis alternativa (Ha) de diferencia en al menos uno de los tratamientos. Los datos fueron analizados mediante el siguiente modelo lineal aditivo:

Y =μ+R +D +P +(RD) +ε ijk i j k ij ijk

Donde:

• Y : Variable respuesta observada en la ijk unidad experimental.

• μ: Media general del experimento.

• R: Efecto fijo de la i-ésima raza del mai cho reproductor (i = Kuri, Perú).

• D: Efecto fijo de la j-ésima densidad de j empadre (j = 6, 7, 8).

• P : Efecto fijo del k-ésimo orden de k parto o repetición temporal (k = I, II).

• (RD) : Efecto de la interacción entre la ij i-ésima raza y la j-ésima densidad de empadre.

• ε : Efecto del error experimental resiijk dual asociado a las observaciones.

Los datos cuantitativos se procesaron mediante un Análisis de Varianza (ANOVA) factorial. En los escenarios donde la interacción Raza × Densidad no resultó estadísticamente significativa, el modelo estadístico se redujo de forma estricta a sus efectos aditivos con la finalidad

de suprimir los sesgos de estimación que induce la suma de cuadrados Tipo III ante la presencia de interacciones irrelevantes. Las comparaciones múltiples de medias se ejecutaron mediante la prueba de Tukey (p < 0,05). Para validar la consistencia del efecto del factor raza sobre la variable discreta de conteo (TCn), se aplicaron de forma complementaria la prueba t de Student, la prueba no paramétrica U de Mann-Whitney y un modelo de regresión de Poisson. El procesamiento estadístico de la información se realizó con el software Infostat versión 2020

RESULTADOS Y DISCUSIÓN

Tasa de fertilidad

De las 42 hembras evaluadas, la tasa de fertilidad (TF) general alcanzó el 90,5 % en el primer parto (38/42) y el 88,1 % en el segundo (37/42). El 83,3 % (35/42) de las hembras concibió exitosamente en ambos ciclos, mientras que el 95,2% (40/42) parió al menos en una oportunidad. Al discriminar por la variable del reproductor, el 100,0 % de las hembras empadradas con machos Kuri parieron en el primer parto, frente al 81,0 % registrado para la raza Perú. En el segundo parto, los índices de fertilidad se situaron en 95,2 % para Kuri y 81,0 % para Perú (Tabla 1).

La TF global obtenida (95,2 % con al menos un parto) se alinea con el comportamiento reproductivo estándar de 95 % reportado por Cahui (2019). Sin embargo, se sitúa ligeramente por debajo del 98 % para la raza Perú y del 100 % para la raza Kuri estimado por Chauca (2022). Estas variaciones son atribuidas comúnmente a factores fisiológicos y ambientales como el retraso en el primer celo, la manifestación de periodos de anestro o condiciones críticas de alimentación y variabilidad climáti-

ca que restringen el óptimo funcional del aparato reproductor (Rosel, 2001).

Por otra parte, se determinó una asociación estadística altamente significativa entre la raza del macho y la fertilidad (χ² = 4,48; p = 0,034; Cochran-Mantel-Haenszel = 6,16; p = 0,013; OR = 11,18). Este comportamiento demuestra que las hembras asignadas a reproductores Kuri poseen una probabilidad significativa-

Tabla 1. Porcentaje de hembras fértiles por raza

Raza n % P

Kuri 10

Perú 8

Tamaño de camada al nacimiento

El análisis de varianza factorial (ANAVA) determinó la ausencia de un efecto de interacción significativo entre los factores raza y densidad de empadre sobre el tamaño de camada al nacimiento (TCn) (F = 0,41; p = 0,664). Debido a esto, el modelo estadístico se redujo a sus efectos aditivos para corregir el sesgo inherente a la suma de cuadrados tipo III ante interacciones no significativas.

El modelo aditivo reveló diferencias estadísticas significativas para el factor principal raza (F = 10,70; p = 0,0016), donde la variedad Kuri superó en 0,93 crías por parto a la raza Perú. El promedio para Kuri se fijó en 3,45 ± 1,21 crías/parto frente a 2,52 ± 1,52 crías/partos obtenidos en Perú, lo que representa una ventaja productiva del +37,0% (Tabla 2).

La variable densidad de empadre no evidenció modificaciones estadísticas significativas sobre el TCn (F = 2,43; p = 0,095). En contraposición, el orden de parto ejerció un efecto significativo (F = 7,66;

mente mayor de concretar un parto exitoso. Es fundamental destacar que la TF no está supeditada de forma exclusiva al potencial de la hembra. Factores del macho como la calidad seminal, el vigor sexual (libido) y la capacidad de servicio alteran los indicadores bajo esquemas de parición continua, donde limitantes adicionales como el estrés térmico pueden comprometer la viabilidad reproductiva (Benavides et al., 2020; Quispe-Ccasa et al., 2023).

el macho y número de parto

rto I % Parto II

,0 95,2

0 81,0

p = 0,0070), constatándose un incremento real de 0,78 crías durante el segundo parto (3,38 ± 1,55 crías/parto) en comparación con el primero (2,60 ± 1,23 crías/parto).

El incremento genotípico en el tamaño de camada al emplear reproductores machos Kuri concuerda con la base conceptual y las evaluaciones de prolificidad descritas por Chauca (2023). Pese a ello, la diferencia del 37,0 % detectada en la provincia de Chota resultó ligeramente inferior al parámetro de ventaja del 41 % difundido originalmente por el MIDAGRI/INIA (2021).

Respecto al orden cronológico de los partos, la mayor productividad cuantitativa del segundo ciclo coincide con las dinámicas documentadas por Rodríguez et al. (2015) y Moya (2024). Este fenómeno responde a la madurez fisiológica uterina de la madre y al desarrollo de sus capacidades homeostáticas y gestacionales tras superar el primer ciclo reproductivo.

Tabla 2. Tamaño de camada de cuyes según la raza del repro ± DE)

Factor evaluado Nivel Interacción Raza ×Densidad Kuri —1 Kuri —1 Kuri —1 Perú —1 Perú —1 Perú —1 Efecto Fijo: Raza Kuri Perú Efecto Fijo: Parto I II

Nota: Medias con letras distintas dentro de una misma sección d (p < 0,05).

Para validar estadísticamente que los efectos detectados correspondían a un comportamiento genético real y descartar sesgos derivados de los niveles de infertilidad observados en el grupo asociado a la raza Perú, se ejecutaron pruebas estadísticas de contraste independientes del tipo de suma de cuadrados lineal. Las metodologías analíticas aplicadas (Prueba t de Student, U de Mann-Whitney y Regresión de Poisson) ratificaron de manera unánime y robusta la significancia estadística crítica de los factores raza y orden de parto sobre el TCn (Tabla 3).

Peso al nacimiento de gazapos por sexo

El factor densidad de empadre (DE) modificó de forma estadísticamente significativa el peso al nacimiento (PN) tanto en ejemplares machos (F = 10,02; p = 0,0001) como en hembras (F = 13,23; p < 0,001). El factor principal raza no demostró efectos significativos sobre el peso de las crías macho (F = 2,45; p = 0,120), pero sí generó diferencias en las crías hembra (F = 10,29; p = 0,0017).

Se determinó una interacción estadística significativa entre los factores raza

uctor, densidad de empadre y orden de parto (Media

n / Σ crías Crías por parto — 3,58 ±1,16 — 3,93 ±1,21 — 2,94 ±1,12 — 2,75 ±1,22 — 2,64 ±1,45 — 2,25 ±1,81 42 / 145 3,45 ±1,21 a 42 / 106 2,52 ±1,52 b 42 / 109 2,60 ±1,23 a 42 / 142 3,38 ±1,55 b

factor indican diferencias significativas según la prueba de Tukey

y densidad de empadre para ambos sexos (machos: F = 4,08, p = 0,019; hembras: F = 3,42, p = 0,036). Esto indica rigurosamente que el impacto de la densidad de manejo pecuario sobre el PN está condicionado por el tipo de raza utilizado y viceversa.

Asimismo, el orden de parto ejerció una influencia altamente significativa en el peso de ambos sexos (machos: F = 69,16, p < 0,001; hembras: F = 44,34, p < 0,001). Se registraron pesos significativamente superiores en el primer parto (machos: 150,82 ± 28,28 g; hembras: 145,86 ± 25,93 g) respecto al segundo parto (machos: 118,12 ± 25,82 g; hembras: 127,23 ± 14,51 g).

Al evaluar las interacciones específicas (Tabla 4), los mayores Pesos al Nacimiento en los machos correspondieron a los tratamientos Kuri × 1:6 (155,00 ± 7,89 g) y Kuri × 1:8 (136,01 ± 5,42 g). En el caso de las hembras, las medias más altas se concentraron en los grupos Kuri × 1:6 (141,49 ± 5,88 g), Kuri × 1:8 (145,12 ± 3,59 g) y Perú × 1:8 (147,32 ± 4,26 g). El peso general de nacimiento integrado para la totalidad de gazapos evaluados (n = 251) se situó en 133,96 ± 1,71 g.

Tabla 3. Modelos estadísticos alternativos de verificación para el efecto de la raza y número de parto sobre el tamaño de camada

Valor de significancia Factor evaluado Metodología estadística Estadístico de prueba (p) ANOVA (Efectos Raza del macho F = 10,70 0,0016 aditivos) (Kuri vs. Perú) Prueba t de Student t = 3,10 0,0027 U de Mann-Whitney U = 1172,5 0,0076 Regresión de Poisson LR = 6,08 0,0136 ANOVA (Efectos Número de parto F = 7,66 0,0070 aditivos) (Parto I vs. Parto II) Prueba t de Student t = 2,58 0,0117 U de Mann-Whitney U = 540,5 0,0017 Regresión de Poisson LR = 4,35 0,0370

Tabla 4. Peso al nacimiento (g) clasificado por sexo bajo efectos de la interacción raza y densidad de empadre (Media ± EE)

Raza Densidad n Peso machos (g) Peso hembras (g)

Kuri 1:6 16/27 155,00 ±7,89 a 141,49 ±5,88 a

Kuri 1:7 31/24 120,81 ±6,94 b 122,75 ±2,83 b

Kuri 1:8 25/22 136,01 ±5,42 a 145,12 ±3,59 a

Perú 1:6 16/17 137,46 ±7,85 a 125,59 ±3,21 b

Perú 1:7 21/16 118,62 ±5,48 b 123,79 ±3,52 b

Perú 1:8 18/18 143,47 ±1,35 a 147,32 ±4,26 a

Nota: Medias con letras distintas dentro de una misma columna difieren significativamente según la prueba de Tukey (p < 0,05).

El notable descenso del PN evidenra transversal entre las razas Kuri y Perú ciado durante el segundo parto responde (Tabla 4) confirma que este parámetro está de forma directa al incremento simultáneo condicionado prioritariamente por la escaverificado en el tamaño de la camada. Esta la física del tamaño de la camada más que respuesta biológica coincide plenamente por determinantes genéticos exclusivos del con los principios de fisiología e interaccioprogenitor. Los pesos promedio obtenidos nes de espacio intrauterino reportados por para la raza Kuri guardan estrecha corres- Wuebker (2021), donde un mayor número pondencia con los estándares productivos de fetos en desarrollo disminuye la disponide líneas genéticas mejoradas reportados bilidad individual de nutrientes y restringe previamente en el país por Carrasco (2019) el peso final al parto. y Cruz et al. (2021).

La marcada homogeneidad en la La ausencia de diferencias signifireducción del PN observada de manecativas entre ambas razas al aislar los aná-

lisis dentro de un mismo número de parto ciclos productivos. Al analizar la distribufundamenta la postura de que el TCn consción de frecuencias (Tabla 5), la mayor tasa tituye el factor crítico determinante del de ocurrencia para el grupo Kuri se concenpeso fetal individual, tal como sostienen tró en camadas de 3 crías durante el parto I investigaciones previas en la especie (Ro- (42,9 %) y de 4 crías en el parto II (52,4 %). dríguez et al., 2015; Rubio et al., 2019 y En contraste, el grupo de empadre con raza Guerrero et al., 2020). Perú concentró sus mayores frecuencias en camadas de 2 y 3 crías para el parto I (33,3 Frecuencia de ocurrencia del tamaño de camada % y 23,8 %) y de 3 y 4 crías en el parto II (42,9 % y 19,0 %). La tasa de nuliparidad (ausencia de partos) se situó en un 4,8 % para las hem- Las frecuencias de camadas numebras empadradas con machos Kuri durante rosas (5 crías) alcanzaron niveles de 9,5 % el parto II, mientras que el grupo asociado a y 14,3 % (partos I y II) en el grupo Kuri, machos Perú reflejó una ausencia de partos frente a un 14,3 % exclusivo en el grupo del 19,0 % de manera constante en ambos de la raza Perú. Por último, la ocurrencia

Tabla 5. Frecuencia de número de crías por parto según tratamiento reproductivo e interacción en cuyes mejoradas

N°Crías K-1:6 K-1:7 K-1:8 P-1:6 P-1:7 P-1:8 Frecuencia Acumulada General 0 0 (0%) 0 (0%) 1 (6,3%) 1 (8,3%) 2 (14,3%) 5 (31,3%) 1 0 (0%) 0 (0%) 0 (0%) 1 (8,3%) 0 (0%) 0 (0%) 2 2 (16,7%) 2 (14,3%) 4 (25%) 1 (8,3%) 4 (28,6%) 3 (18,8%) 3 4 (33,3%) 3 (21,4%) 5 (31,3%) 6 (50%) 4 (28,6%) 4 (25%) 4 4 (33,3%) 4 (28,6%) 6 (37,5%) 3 (25%) 3 (21,4%) 2 (12,5%) 5 1 (8,3%) 4 (28,6%) 0 (0%) 0 (0%) 1 (7,1%) 2 (12,5%) 6 1 (8,3%) 1 (7,1%) 0 (0%) 0 (0%) 0 (0%) 0 (0%) Frecuencia en Parto I 0 0 (0%) 0 (0%) 0 (0%) 1 (16,7%) 0 (0%) 3 (37,5%) 1 0 (0%) 0 (0%) 0 (0%) 1 (16,7%) 0 (0%) 0 (0%) 2 1 (16,7%) 2 (28,6%) 4 (50%) 1 (16,7%) 3 (42,9%) 3 (37,5%) 3 4 (66,7%) 2 (28,6%) 3 (37,5%) 2 (33,3%) 1 (14,3%) 2 (25%) 4 1 (16,7%) 1 (14,3%) 1 (12,5%) 1 (16,7%) 3 (42,9%) 0 (0%) 5 0 (0%) 2 (28,6%) 0 (0%) 0 (0%) 0 (0%) 0 (0%) 6 0 (0%) 0 (0%) 0 (0%) 0 (0%) 0 (0%) 0 (0%) Frecuencia en Parto II 0 0 (0%) 0 (0%) 1 (12,5%) 0 (0%) 2 (28,6%) 2 (25%) 1 0 (0%) 0 (0%) 0 (0%) 0 (0%) 0 (0%) 0 (0%) 2 1 (16,7%) 0 (0%) 0 (0%) 0 (0%) 1 (14,3%) 0 (0%) 3 0 (0%) 1 (14,3%) 2 (25%) 4 (66,7%) 3 (42,9%) 2 (25%) 4 3 (50%) 3 (42,9%) 5 (62,5%) 2 (33,3%) 0 (0%) 2 (25%) 5 1 (16,7%) 2 (28,6%) 0 (0%) 0 (0%) 1 (14,3%) 2 (25%) 6 1 (16,7%) 1 (14,3%) 0 (0%) 0 (0%) 0 (0%) 0 (0%)

extrema de camadas de 6 crías se registró únicamente en el segundo parto de las hembras empadradas con Kuri (9,5 %) (Figura 1).

La mayor recurrencia de partos múltiples numerosos, específicamente en los rangos de 5 y 6 crías observados en el lote Kuri, se asocia directamente a la elevada tasa de prolificidad que caracteriza a este genotipo compuesto bajo condiciones controladas de densidad pecuaria. La integración analítica de los datos sugiere que los componentes de manejo zootécnico (densidad de empadre), la base genética (raza del macho reproductor) y las variaciones intrínsecas del tamaño de camada deben interactuar de manera coordinada para optimizar los pesos al nacimiento y elevar los techos de eficiencia en los sistemas comerciales de cría de la región.

CONCLUSIONES

El genotipo del reproductor macho constituye un factor determinante sobre el tamaño de la camada y la distribución del peso al nacimiento de los gazapos en la provincia de Chota, Región Cajamarca,

demostrando que la introducción de la raza Kuri incrementa significativamente la eficiencia productiva frente a la raza Perú.

La utilización de machos de la raza Kuri genera una ventaja zootécnica crítica del +37,0% en el tamaño de camada al nacimiento en comparación con los machos de la raza Perú, alcanzando un promedio de 3,45 ± 1,21 crías por parto.

Existe un efecto de interacción estadísticamente significativo entre la raza del reproductor y la densidad de empadre sobre el peso al nacimiento de las crías de ambos sexos, identificándose que los mayores pesos individuales se maximizan bajo relaciones de manejo pecuario de 1:6 y 1:8 empleando la raza Kuri en hembras mejoradas.

El orden de parto ejerce una influencia inversa entre la prolificidad y el desarrollo fetal; el segundo parto consolida una mayor productividad cuantitativa en el número de crías (+0,78 gazapos), pero induce un descenso consecuente en el peso individual al nacimiento debido a restricciones de espacio intrauterino, independientemente de la raza del macho.

Declaración de roles de autores: • Benito Bautista Espinoza: Conceptualización e implementación del proyecto; determinación de la metodología y administración; registro y tabulación de datos. • Julia Mercedes Ramírez Sánchez: Supervisión; apoyo en la determinación de la metodología; redacción, investigación, revisión y edición. • Napoleón Corrales Roriguez: Análisis estadístico de los datos, análisis formal, búsqueda de información e investigación, redacción y borrador original.