En numerosos países con clima tropical, el cruzamiento de Cebú con razas Bos taurus (Holstein, Jersey, Suizo Pardo) ha sido una práctica común en hatos de doble propósito, para mejorar rasgos de composición y producción de leche, salud, supervivencia y fertilidad. Entre estas características, se ha observado que la fertilidad tiene el mayor impacto en la eficiencia del hato.

El cruzamiento permite la introducción de genes favorables y obtiene ventaja de la complementariedad racial y la heterósis. La complementariedad racial permite a los criadores capitalizar las fortalezas de las diferentes razas, ya que ninguna raza es superior en todas las características que afectan la rentabilidad. En México, el sistema de producción de doble propósito está constituido principalmente por cruzas entre Cebú y Suizo Pardo, Holstein Friesian y Simmental (1,2,3).

Varios estudios realizados en diferentes regiones del mundo han evaluado el desempeño reproductivo de vacas Holstein Friesian x Cebú (4,5,6,7,8,9) o Jersey x Cebú (10,11), o han comparado el desempeño reproductivo de vacas puras indígenas y cruzadas Holstein Friesian x Cebú (12,13), de vacas Holstein Friesian x Cebú y Jersey x Cebú (12,14,15), de vacas Cebú y 3/8 Simmental x 5/8 Cebú (16), de vacas Suizo Pardo x Cebú y Holstein Friesian x Cebú (17), y de vacas Holstein Friesian x Cebú y Simmental x Cebú (18,19). Sin embargo, son muy escasos los artículos científicos relacionados con la evaluación reproductiva de vacas Suizo Pardo x Cebú y Simmental x Cebú bajo las mismas condiciones tropicales o subtropicales; aparentemente, solo se ha publicado un artículo relacionado con la comparación de estos dos últimos genotipos (1). Un estudio mexicano previo mostró que vacas Simmental x Cebú y Suizo Pardo x Cebú produjeron niveles similares de leche (3).

Con base en estos antecedentes, el objetivo del presente estudio fue comparar la fertilidad de vacas Suizo Pardo x Cebú y Simmental x Cebú criadas en un ambiente tropical.

Sitio de estudio. El estudio se implementó en un hato bovino de doble propósito en la estación experimental Playa Vicente, que pertenece al Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias (INIFAP), y se localiza a 95 msnm, en los 17° 19’ latitud norte y 95° 41’ longitud oeste, en Veracruz, México. La región cuenta con clima tropical húmedo, temperatura media de 26.8°C y precipitación de 2200 mm al año (20).

Animales. En la tabla 1 se presenta el número de vacas evaluadas y de los padres y madres que se utilizaron para producirlas. Se evaluaron características reproductivas de 185 vacas cruzadas Suizo Pardo x Cebú y Simmental x Cebú con diversos porcentajes de raza europea (50.0, 62.5 o 75.0%).

Tabla 1 Tabla 1 Número de vacas, padres y madres, por grupo genético .

Grupo genético	Vacas	Padres	Madres
Cebú	-	-	65
Simmental	-	21	-
Suizo Pardo	-	24	-
½ Simmental x ½ Cebú	36	2	21
½ Suizo Pardo x ½ Cebú	47	1	28
¾ Simmental x ¼ Cebú	35	2	4
¾ Suizo Pardo x ¼ Cebú	45	6	7
⅝ Simmental x ⅜ Cebú	7	-	-
⅝ Suizo Pardo x ⅜ Cebú	15	4	4
Total	185	60	129

Las vacas evaluadas nacieron de 1981 a 2003 y se produjeron con 60 toros puros e híbridos con diversos porcentajes de raza Suizo Pardo o Simmental y 129 madres puras e híbridas con diversos porcentajes de raza Cebú. La raza Cebú de las madres puras fue Indubrasil. Las 129 madres se aparearon con los 60 toros mediante IA (en la gran mayoría de los casos) y monta natural. Las vacas de ambos grupos genéticos fueron manejadas juntas, de la misma manera, y pastorearon en los mismos potreros.

La tabla 2 presenta la forma en que se realizaron los cruzamientos para producir las 185 vacas evaluadas. Las vacas F1 se generaron con sementales Simmental o Suizo Pardo. Las vacas 75% europeo-25% Cebú se generaron con sementales puros de raza europea y sementales 75% europeo-25% Cebú; las vacas 62.5% europeo-37.5% Cebú se produjeron con sementales F1 y 75% europeo-25% Cebú. Además, las vacas 62.5% Suizo Pardo-37.5% Bos indicus se generaron con sementales 62.5% europeo-37.5% Cebú.

Tabla 2 Tabla 2 Apareamientos realizados para producir las vacas Bos taurus x Bos indicusa

Padre	Madre
	C	1S1C	1P1C	3S1C	3P1C	5P3C
S	50%S-50%C	75%S-25%C
P	50%P-50%C		75%P-25%C
1S1C				62.5%S-37.5%C
1P1C					62.5%P-37.5%C
3S1C		62.5%S-37.5%C		75%S-25%C
3P1C			62.5%P-37.5%C		75%P-25%C
5P3C						62.5%P-37.5%C
aC= Cebú, S= Simmental, P= Suizo Pardo, 1S1C= ½ Simmental x ½ Cebú, 1P1C= ½ Suizo Pardo x ½ Cebú, 3S1C= ¾ Simmental x ¼ Cebú, 3P1C= ¾ Suizo Pardo x ¼ Cebú, 5P3C= ⅝ Suizo Pardo x ⅜ Cebú

Alimentación. Las vacas se mantuvieron en pastoreo rotacional en praderas con Guinea (Panicum maximum) y Estrella Africana (Cynodon plectostachyus). Además, consumieron 2 kg/animal/d de un suplemento comercial con 70% de TND y 16% de PC, 30 días antes de parir y en cada ordeño hasta el secado. En el periodo de estiaje, las vacas se suplementaron con 15-20 kg/animal/d de caña japonesa picada (Saccharum sinense) o ensilaje de maíz (Zea mayz).

Inseminación artificial. Las vacas se inseminaron por primera vez cuando alcanzaron un peso corporal de 350 kg. La observación de celos se efectuó durante 1 h en la mañana (de 06:00 a.m. a 07:00 a.m.) y 1 h en la tarde (de 05:00 p.m. a 06:00 p.m.), con la ayuda de un toro celador. La eficiencia de detección del celo fue 60%. La inseminación de las vacas fue de la siguiente manera: aquellas que presentaron celo en la mañana se sirvieron en la tarde, y aquellas que lo presentaron en la tarde se inseminaron al siguiente día por la mañana, 12 h después de observado el celo, aproximadamente. La gestación se diagnosticó por palpación rectal 45 d después de la inseminación. Las vacas se desecharon principalmente por razones de pobre fertilidad.

Ordeño. El ordeño fue mecánico, dos veces al día, después de un amamantamiento ligero realizado por el becerro, para estimular la eyección de leche. Los becerros se mantuvieron atados, a un costado de las vacas, mientras ellas fueron ordeñadas. La producción de leche de cada vaca se registró en cada ordeño. Durante los tres primeros meses de lactancia, se ordeñaron tres cuartos, dejando uno para consumo de leche del becerro, más la leche residual de los tres cuartos ordeñados.

A partir de los 91 días de lactancia hasta que el becerro se destetó, se ordeñaron los cuatro pezones, dejando para la alimentación del becerro solamente la leche residual. Después del ordeño, se permitió a los becerros mamar 1 h, aproximadamente. Después de ello, los becerros se separaron de las vacas. Posterior al destete, los becerros solo se usaron para estimular la salida de leche hasta el secado de las vacas, que se realizó a los 7 meses de gestación o cuando produjeron menos de 3 kg de leche por día.

Variables de respuesta. Se analizaron registros de edad a primer parto, periodo parto-primer servicio, días abiertos, periodo interparto, servicios por concepción, tasa de gestación y peso al parto. Días abiertos se definió como el intervalo (días) del parto a la concepción; esta característica refleja la tasa de parto y la capacidad de la hembra para ciclar y manifestar estro. Tasa de gestación a primer servicio se definió como una variable binaria; por lo tanto, si una hembra quedó gestante después del servicio, se asignó un valor de 1; de lo contrario, se asignó un valor de 0.

Análisis estadísticos. Edad a primer parto se analizó con el procedimiento MIXED de SAS (21) con un modelo simple que incluyó vaca anidada en grupo genético de la vaca, y padre de la vaca anidado en grupo genético del padre como efectos aleatorios, y época de parto, grupo genético, y año de parto como efectos fijos.

Las características restantes se analizaron con un modelo de mediciones repetidas que incluyó vaca anidada en grupo genético de la vaca, y padre de la vaca anidado en grupo genético del padre como efectos aleatorios (excepto para tasa de gestación y servicios por concepción), y año de parto, número de lactancia, grupo genético de la vaca y época de parto como efectos fijos. Además, para tasa de gestación a primer servicio, el modelo incluyó fase de lactancia (Fase 1: de 1 a 50 d; Fase 2: de 51 a 100 d; Fase 3: de 101 a 150 d; y Fase 4: ≥151 d postparto).

Días abiertos, peso al parto, periodo parto-primer servicio y periodo interparto se analizaron con PROC MIXED de SAS (21); tasa de gestación y servicios por concepción se analizaron con PROC GENMOD del mismo programa. Para servicios por concepción, se especificó una distribución Poisson en la declaración del modelo; en el análisis estadístico de la tasa de gestación se especificó una distribución binomial y se usó una función liga logit.

El modelo estadístico para analizar peso al parto, días abiertos, periodo parto-primer servicio y periodo interparto se ajustó preliminarmente probando diferentes estructuras de covarianzas (ante-dependencia, auto-regresiva de primer orden, auto-regresiva heterogénea, simetría compuesta, simetría compuesta heterogénea, simple, Toeplitz, Toeplitz heterogénea y sin estructura) para proporcionar el mejor ajuste a los datos.

Las estructuras de covarianzas que se probaron para analizar tasa de gestación y servicios por concepción fueron auto-regresiva de primer orden, simetría compuesta, independiente, Toeplitz y sin estructura. La selección de la estructura de covarianzas apropiada para días abiertos, periodo parto-primer servicio, peso al parto y periodo interparto se basó en los criterios de información de ajuste de Akaike, de segundo orden y Bayesiano de Schwarz. Para tasa de gestación y servicios por concepción, la estructura de covarianzas apropiada se seleccionó con base en el criterio de información de ajuste de cuasi-verosimilitud.

La tabla 3 resume las características de los datos para todas las variables de respuesta. Las medias no ajustadas de periodo interparto, servicios por concepción, edad a primer parto, periodo parto-primer servicio, días abiertos, tasa de gestación a primer servicio y peso al parto fueron: 447.5 d, 2.1 servicios, 36.1 meses, 119.5 d, 162.1 d, 43.9%, y 482.4 kg, respectivamente.

Tabla 3 Tabla 3 Estadísticas descriptivas.

Variablea	n	Media	Dev Est	CV (%)	M í nimo	M á ximo
EPP (meses)	176	36.1	7.2	19.9	23.4	56.0
PPS (días)	641	119.5	78.7	65.9	7	906
SXC	595	2.1	1.3	61.9	1	8
DA (días)	596	162.1	104.9	64.7	7	926
TG (%)	595	43.9	49.7	113.2	0	100
PIP (días)	596	447.5	100.1	22.4	298	975
PP (kg)	655	482.4	74.0	15.3	285	770
aEPP= edad a primer parto; PPS= periodo parto-primer servicio; SXC= servicios por concepción; DA= días abiertos; TG= tasa de gestación a primer servicio; PIP= periodo interparto; PP= peso al parto.

Los criterios de información de ajuste de Akaike, de segundo orden, Bayesiano de Schwarz y de cuasi-verosimilitud se presentan en la tabla 4.

Table 4 Table 4 Criterios de información de ajuste de Akaike (AIC), segundo orden (AICC), Bayesiano de Schwarz (BIC) y cuasi-verosimilitud (QIC).

Variable/Estructura de covarianzas	Criterio de informacióna
	AIC	AICC	BIC	QIC
Servicios por concepciónb
Auto-regresiva de primer orden	-	-	-	879.33
Simetría compuesta	-	-	-	881.01
Independiente	-	-	-	880.28
Toeplitz	-	-	-	881.54
Tasa de gestación a primer serviciob
Auto-regresiva de primer orden	-	-	-	825.51
Simetría compuesta	-	-	-	825.37
Independiente	-	-	-	825.46
Toeplitz	-	-	-	825.52
Periodo parto-primer servicioc
Simple	6584.2	6584.2	6593.9	-
Auto-regresiva de primer orden	6585.5	6585.6	6598.6	-
Días abiertosc
Simple	6375.6	6375.7	6385.4	-
Auto-regresiva de primer orden	6377.6	6377.7	6390.6	-
Periodo interpartoc
Simple	6305.4	6305.5	6315.2	-
Auto-regresiva de primer orden	6304.0	6304.1	6317.0	-
Simetría compuesta	6317.4	6317.5	6327.2	-
Peso al partoc
Simple	6288.2	6288.3	6298.0	-
Auto-regresiva de primer orden	6284.1	6284.1	6297.1	-
aValores menores indican mejor ajuste. bLa covarianza sin estructura no fue estimable. cLas covarianzas ante-dependencia, sin estructura, auto-regresiva heterogénea, simetría compuesta, simetría compuesta heterogénea, Toeplitz y Toeplitz heterogénea no fueron estimables.

Las estructuras de covarianzas apropiadas usadas en los modelos definitivos fueron: auto-regresiva de primer orden para peso al parto y servicios por concepción; simetría compuesta para tasa de gestación; y simple para periodo interparto, periodo parto-primer servicio y días abiertos.

La significación estadística de los efectos genéticos y ambientales se muestra en la tabla 5. La vaca anidada en grupo genético de la vaca fue una fuente de variación significativa para peso al parto (p<0.0001), días abiertos (p<0.0075), edad a primer parto (p<0.0001), periodo parto-primer servicio (p<0.0244) y periodo interparto (p<0.0072). Padre de la vaca anidado en grupo genético del padre fue significativo para peso al parto (p<0.0206), días abiertos (p<0.0107), periodo parto-primer servicio (p<0.0278) y periodo interparto (p<0.0125). El grupo genético explicó la variación en días abiertos (p<0.0023), periodo interparto (p<0.0018), periodo parto-primer servicio (p<0.0002) y peso al parto (p<0.0008). Número de lactancia fue una fuente de variación significativa para días abiertos (p<0.0001), peso al parto (p<0.0001), periodo interparto (p<0.0001), periodo parto-primer servicio (p<0.0001) y tasa de gestación (p<0.0130).

Tabla 5. Tabla 5. Niveles de significancia estadística de los efectos fijos y aleatorios para las variables de respuesta.

Efecto	Variable de respuestaa
	EPP	PPS	SXC	DA	TG	PIP	PP
Vacab	<0.0001	0.0244	---	0.0075	---	0.0072	<0.0001
Padrec	0.0623	0.0278	---	0.0107	---	0.0125	0.0206
Grupo genético	0.1315	0.0002	0.8939	0.0023	0.3617	0.0018	0.0008
Año (A)	<0.0001	<.0001	0.0042	<0.0001	0.0358	<0.0001	<0.0001
Época (E)	0.1155	0.6925	0.3036	0.7028	0.8084	0.9605	0.2618
Lactanciad	---	<.0001	0.5058	<0.0001	0.0130	<0.0001	<0.0001
A x E	---	0.0392	---	0.0002	---	0.0001	0.0010
Fasee	---	---	---	---	0.4262	---	---
aEPP= edad a primer parto; PPS= periodo parto-primer servicio; SXC= servicios por concepción; DA= días abiertos; TG= tasa de gestación a primer servicio; PIP= periodo interparto; PP= peso al parto; bVaca anidada en grupo genético de la vaca; cPadre de la vaca anidado en grupo genético del padre; dLactancia= número de lactancia; eFase= Fase de la lactancia.

Las medias ajustadas para las variables de respuesta analizadas se muestran en la taba 6. Las vacas Simmental x Cebú se sirvieron después del parto 39 días antes (p<0.05), tuvieron 47 días abiertos menos (p<0.05) y periodo interparto 45 días más corto (p<0.05) y fueron 34 kg más pesadas (p<0.05) al parto que las vacas Suizo Pardo x Cebú.

Tabla 6 Tabla 6 Medias de cuadrados mínimos y errores estándar para edad a primer parto (EPP), periodo parto-primer servicio (PPS), servicios por concepción (SXC), días abiertos (DA), tasa de gestación a primer servicio (TG), periodo interparto (PIP) y peso al parto (PP), por grupo genético, época de parto, número de lactancia y fase de la lactancia.

	Variable de respuesta
	EPP	PPS	SXC	DA	TG	PIP	PP
Grupo genéticod
P x C	36.5±.62a	143±6.2a	2.1±.08a	192±9.6a	41±3a	476±9.1a	451±6.3a
S x C	38.0±.82a	104±7.0b	2.1±.09a	145±10.8b	45±4a	431±10.3b	485±7.2b
Época de parto
Fría	38.1±.79a	120±6.3a	2.1±.10a	167±9.5a	43±5a	453±8.9a	470±5.5a
Seca	37.5±.75a	126±6.7a	2.1±.12a	173±9.8a	41±4a	455±9.3a	464±5.6a
Lluviosa	36.3±.72a	124±5.7a	2.0±.09a	165±8.4a	45±4a	452±7.9a	470±5.2a
Número lactancia
1		144±6.3a	2.2±.10a	200±9.2a	34±4a	487±8.7a	428±5.3a
2		125±6.9b	2.0±.12a	161±9.9b	52±5b	445±9.4b	469±5.5b
≥3		102±5.8c	2.0±.08a	144±8.7b	44±4ab	428±8.3b	506±5.5c
Fase de lactanciae
1					37±6a
2					45±4a
3					43±4a
4					48±5a
a,b ,cMedias con diferente superíndice dentro de la misma columna en cada factor son diferentes (P<0.05). dP x C= Suizo Pardo x Cebú; S x C= Simmental x Cebú. e1= de 1 a 50 d; 2= de 51 a 100 d; 3= de 101 a 150 d; 4= ≥151 d postparto.

Las vacas con tres lactancias tuvieron menos (p<0.05) días a primer servicio después del parto que las de una o dos lactancias. Las vacas con dos lactancias tuvieron menos (p<0.05) días a primer servicio después del parto que las de una lactancia. Las vacas con dos o tres lactancias tuvieron menos días abiertos (p<0.05) e intervalos entre partos más cortos (p<0.05) que las de una lactancia. Las vacas con dos lactancias tuvieron mayor (p<0.05) tasa de gestación a primer servicio que las de una lactancia; la tasa de gestación de las vacas con tres lactancias fue intermedia. Las vacas con tres lactancias fueron más pesadas al parto (p<0.05) que las de una o dos lactancias. Las vacas con dos lactancias fueron más pesadas al parto (p<0.05) que las de una lactancia.

Artículos científicos que comparan la capacidad reproductiva de vacas Simmental x Cebú y Suizo Pardo x Cebú son muy escasos en la literatura, por lo tanto, en la mayoría de los casos, vacas Simmental x Cebú fueron comparadas con vacas Holstein Friesian x Cebú; sin embargo, artículos científicos que comparan el desempeño reproductivo de estos dos últimos genotipos son también limitados.

Las medias de cuadrados mínimos para edad a primer parto y periodo interparto reportadas para los grupos genéticos Simmental x Cebú y Suizo Pardo x Cebú del presente estudio son similares a las reportadas para el grupo genético Suizo Pardo x Cebú en condiciones tropicales de Yucatán, México (22), y para el grupo genético Holstein Friesian x Cebú en granjas lecheras de pequeños agricultores de Etiopía (5).

En un estudio realizado en el trópico húmedo de México, los genotipos Suizo Pardo x Cebú y Simmental x Cebú tuvieron similares edades a primer parto y periodos interparto (1). Estos resultados concuerdan parcialmente con los del presente estudio, en el cual los genotipos Simmental x Cebú y Suizo Pardo x Cebú tuvieron similar edad a primer parto (p>0.05), pero el genotipo Simmental x Cebú tuvo periodos interparto más cortos (p<0.05) que el genotipo Suizo Pardo x Cebú.

Investigadores colombianos (19) observaron que el intervalo del parto al primer calor de vacas Simmental x Cebú fue 48.8 días más corto que el de vacas Holstein x Cebú; sin embargo, estos investigadores también encontraron que estos grupos genéticos no difirieron significativamente en la edad a primer calor. En un estudio realizado en Etiopía (18) se encontró que el intervalo entre partos de vacas ¾ Simmental x ¼ Horro (Cebú) fue 178 días más corto que el de vacas ¾ Friesian x ¼ Horro (Cebú). En contraste, en un estudio realizado en Brasil (23) vacas Simmental x Cebú y Holstein Friesian x Cebú tuvieron similares intervalos entre partos.

La menor fertilidad de las vacas Suizo Pardo x Cebú encontrada en el presente estudio se puede deber a una mayor intensidad de selección para producción de leche en la raza Suizo Pardo que en la raza Simmental. Se ha demostrado que la producción de leche tiene una correlación genética desfavorable con la fertilidad, lo que indica que la selección para mayor producción de leche disminuye la fertilidad (24).

En varios estudios, la evaluación reproductiva de vacas cruzadas Bos taurus x Bos indicus ha revelado que vacas con dos o tres lactancias tienen intervalos entre partos más cortos que vacas con una lactancia (1,13,15,22), lo cual concuerda con los presentes hallazgos. Esto se explica, en parte, por el hecho de que las vacas primerizas no han completado su desarrollo corporal, por lo que parte del alimento que consumen lo destinan a este fin, sacrificando su producción de leche y fertilidad (22). En contraste, en un estudio realizado en Tabasco, México (9), con vacas Holstein Friesian x Cebú, vacas con tres lactancias tuvieron intervalos entre partos más largos que vacas con una lactancia (440 vs 414 días), hallazgo que no concuerda con el presente resultado; también en discordancia están los resultados obtenidos con vacas Jersey x Red Sindhi en la India, donde vacas con una, dos o tres lactancias tuvieron similares intervalos entre partos (11). La discrepancia con este último estudio puede deberse a que las hembras de reemplazo tuvieron una excelente alimentación, lo cual resultó en un buen desarrollo corporal antes de iniciar su vida reproductiva.

En la presente investigación, la tasa de gestación a primer servicio no fue afectada por la fase de lactancia. Por el contrario, en un cruce dialélico con Holstein Friesian y Suizo Pardo (25), vacas en Fase 4 de la lactancia (≥151 d postparto) tuvieron mayor tasa de gestación a primer servicio que vacas en Fase 1 (de 1 a 50 d postparto) y 2 (de 51 a 100 d postparto) de la lactancia (63% vs 44 y 50%, respectivamente), y vacas en Fase 3 (de 101 a 150 d postparto) fueron superiores a vacas en Fase 1 (56% vs 44%); sin embargo, la discrepancia entre estudios pudo deberse a diferencias en los genotipos evaluados y prácticas de manejo del hato. En hatos lecheros, los becerros son separados de sus madres 4-5 d después del parto, pero en hatos de doble propósito los becerros se mantienen atados, a un costado de sus madres, mientras ellas son ordeñadas. Además, durante los primeros meses de lactancia, solo se ordeñan tres pezones, dejando un pezón para consumo de leche del becerro, más la leche residual de los tres pezones ordeñados. La presencia del becerro y el amamantamiento alteran la interacción entre hipotálamo, hipófisis y ovarios, disminuyendo la liberación de GnRH, lo cual conduce a insuficientes pulsos de LH, evitando que se presente la ovulación (26).

En conclusión, las vacas Simmental x Cebú fueron más pesadas al parto y tuvieron intervalos del parto al primer servicio, del parto a la gestación, y del parto al próximo parto más cortos que las Suizo Pardo x Cebú. Esto sugiere que vacas Simmental x Cebú son una mejor alternativa para la producción pecuaria en las regiones tropicales que vacas Suizo Pardo x Cebú.

Conflicto de interés

Los autores no tienen ningún conflicto de intereses.