Longitudinal genetic dynamics of weaning index and implications for cow-calf production efficiency.

In beef cattle, the selection for higher weights at young ages has been questioned with the argument that this criterion may increase the adult weight of cows, resulting in higher costs. Therefore, selection criteria should be employed to increase weights at young ages with minimal impact on the adult weight of cows. Additionally, the relationship between measures of cow production efficiency and other well-established selection criteria in breeding programs remains poorly understood. The objective of this study was to longitudinally evaluate the relationship between the weaning index (WIndex) as a measure of efficiency and growth traits of the cows. Possible changes over time in WIndex due to selection applied for yearling weight (YW) were also investigated. The WIndex was proposed to maximize genetic response in the weaning weight of the calf while maintaining genetic gain in BW of the cow at zero. A random regression model was adopted to estimate correlations between WIndex, BW, hip height (HH), and body condition score (BCS) using records of Nelore cows from three lines. Genetic trends were calculated for the control line (NeC) and lines selected for greater YW (NeS and NeT). The age of 3 years was the most critical for the weaning efficiency of the cows. At this stage, young cows are still growing and wean lighter calves than their adult counterparts. The genetic correlation estimates between WIndex and BW (-0.58 to 0.04), HH (-0.05 to -0.34), and BCS (-0.51 to -0.17) were close to zero or negative. BW and HH were strongly correlated genetically across all ages (0.73-0.76). In general, HH exhibited a weak and negative genetic relationship with BCS. The genetic correlation between BW and BCS was stronger for advanced ages (0.45-0.68). In lines selected for YW, important increases in WIndex were observed. However, NeS has been selected since the 1980s until the present for YW, and thus, it showed a more pronounced trend of increasing BW and, consequently, a more modest trend of increasing WIndex compared to NeT. In contrast, WIndex exhibited a trend close to zero for NeC. In this context, monitoring HH and BCS can be useful to avoid losses in the weaning efficiency of cows. Furthermore, we suggest that one way to mitigate efficiency losses in calf production could involve stabilizing the BW of cows and increasing the weaning weight of calves using the WIndex.

Assessment of runs of homozygosity, heterozygosity-rich regions and genomic inbreeding estimates in a subpopulation of Guzerá (Bos indicus) dual-purpose cattle.

For decades, inbreeding in cattle has been evaluated using pedigree information. Nowadays, inbreeding coefficients can be obtained using genomic information such as runs of homozygosity (ROH). The aims of this study were to quantify ROH and heterozygosity-rich regions (HRR) in a subpopulation of Guzerá dual-purpose cattle, to examine ROH and HRR islands, and to compare inbreeding coefficients obtained by ROH with alternative genomic inbreeding coefficients. A subpopulation of 1733 Guzerá animals genotyped for 50k SNPs was used to obtain the ROH and HRR segments. Inbreeding coefficients by ROH (FROH ), by genomic relationship matrix based on VanRaden's method 1 using reference allele frequency in the population (FGRM ), by genomic relationship matrix based on VanRaden's method 1 using allele frequency fixed in 0.5 (FGRM_0.5 ), and by the proportion of homozygous loci (FHOM ) were calculated. A total of 15,660 ROH were identified, and the chromosome with the highest number of ROH was BTA6. A total of 4843 HRRs were identified, and the chromosome with the highest number of HRRs was BTA23. No ROH and HRR islands were identified according to established criteria, but the regions closest to the definition of an island were examined from 64 to 67 Mb of BTA6, from 36 to 37 Mb of BTA2 and from 0.50 to 1.25 Mb of BTA23. The genes identified in ROH islands have previously been associated with dairy and beef traits, while genes identified on HRR islands have previously been associated with reproductive traits and disease resistance. FROH was equal to 0.095 ± 0.084, and its Spearman correlation with FGRM was low (0.44) and moderate-high with FHOM (0.79) and with FGRM_0.5 (0.80). The inbreeding coefficients determined by ROH were higher than other cattle breeds' and higher than pedigree-based inbreeding in the Guzerá breed obtained in previous studies. It is recommended that future studies investigate the effects of inbreeding determined by ROH on the traits under selection in the subpopulation studied.

Insights on sperm assays and cryopreservation in six Neotropical pit vipers.

Breeding snakes in captivity has become more and more relevant due not only to the growing interest on their venoms but also to the increasing number of endangered species worldwide. Unfortunately, studies on the formation of germplasm banks for these reptiles do not follow the same pace, and literature on sperm cryopreservation remains in its infancy when compared to other taxa. Herein, we first validated a sperm-egg binding assay (using chicken egg perivitelline membrane - EPM) and some nonfluorescent staining techniques for semen analysis of two pit viper genera (Bothrops and Crotalus), and then we investigated the protective effects of dimethylacetamide (DMA), dimethylformamide (DMF), and dimethylsulfoxide (DMSO) at different concentrations (3, 6 and 12%) throughout the freezing process in five species of lancehead and one of rattlesnake. Our validation process showed high correlations among sperm functional tests (including sperm-binding to EPM) and motion parameters. A total of 166 fresh ejaculates were acquired from 233 collection attempts, and 63.9% of these samples exhibited minimal motility for freezing (≥20%). During cryopreservation we observed that post-thaw motility and quality was improved by higher levels of cryoprotectants (CPA), regardless the CPA type. Lower concentrations of CPA were less harmful to sperm motility and progressive motility following the equilibrium phase, but were ineffective in protecting these cells from the freeze-thaw cycle. Likewise, higher CPA concentrations increased post-thaw integrity of the acrosome and plasma membrane for most species, except for rattlesnakes in which only 12% DMSO produced better outcomes.

Insights on sperm assays and cryopreservation in six Neotropical pit vipers

Breeding snakes in captivity has become more and more relevant due not only to the growing interest on their venoms but also to the increasing number of endangered species worldwide. Unfortunately, studies on the formation of germplasm banks for these reptiles do not follow the same pace, and literature on sperm cryopreservation remains in its infancy when compared to other taxa. Herein, we first validated a sperm-egg binding assay (using chicken egg perivitelline membrane – EPM) and some nonfluorescent staining techniques for semen analysis of two pit viper genera (Bothrops and Crotalus), and then we investigated the protective effects of dimethylacetamide (DMA), dimethylformamide (DMF), and dimethylsulfoxide (DMSO) at different concentrations (3, 6 and 12%) throughout the freezing process in five species of lancehead and one of rattlesnake. Our validation process showed high correlations among sperm functional tests (including sperm-binding to EPM) and motion parameters. A total of 166 fresh ejaculates were acquired from 233 collection attempts, and 63.9% of these samples exhibited minimal motility for freezing (≥20%). During cryopreservation we observed that post-thaw motility and quality was improved by higher levels of cryoprotectants (CPA), regardless the CPA type. Lower concentrations of CPA were less harmful to sperm motility and progressive motility following the equilibrium phase, but were ineffective in protecting these cells from the freeze-thaw cycle. Likewise, higher CPA concentrations increased post-thaw integrity of the acrosome and plasma membrane for most species, except for rattlesnakes in which only 12% DMSO produced better outcomes.

Physical and ecological isolation contribute to maintain genetic differentiation between fire salamander subspecies.

Landscape features shape patterns of gene flow among populations, ultimately determining where taxa lay along the continuum between panmixia to complete reproductive isolation. Gene flow can be restricted, leading to population differentiation in two non-exclusive ways: "physical isolation", in which geographic distance in combination with the landscape features restricts movement of individuals promoting genetic drift, and "ecological isolation", in which adaptive mechanisms constrain gene flow between different environments via divergent natural selection. In central Iberia, two fire salamander subspecies occur in parapatry across elevation gradients along the Iberian Central System mountains, while in the adjacent Montes de Toledo Region only one of them occurs. By integrating population and landscape genetic analyses, we show a ubiquitous role of physical isolation between and within mountain ranges, with unsuitable landscapes increasing differentiation between populations. However, across the Iberian Central System, we found strong support for a significant contribution of ecological isolation, with low genetic differentiation in environmentally homogeneous areas, but high differentiation across sharp transitions in precipitation seasonality. These patterns are consistent with a significant contribution of ecological isolation in restricting gene flow among subspecies. Overall, our results suggest that ecological divergence contributes to reduce genetic admixture, creating an opportunity for lineages to follow distinct evolutionary trajectories.

Genotype by environment interaction due to heat stress during gestation and postpartum for milk production of Holstein cattle.

Remarkable increases in the production of dairy animals have negatively impacted their tolerance to heat stress (HS). The evaluation of the effect of HS on milk yield is based on the direct impact of HS on performance. However, in practical terms, HS also exerts its influence during gestation (indirect effect). The main purpose of this study was to identify and characterize the genotype by environment interaction (G × E) due to HS during the last 60 days of gestation (THI_g) and also the HS postpartum (THI_m) over first lactation milk production of Brazilian Holstein cattle. A total of 389 127 test day milk yield (TD) records from 1572 first lactation Holstein cows born in Brazil (daughters of 1248 dams and 70 sires) and the corresponding temperature-humidity index (THI) obtained between December 2007 and January 2013 were analyzed using different random regression models. Cows in the cold environment (THI_g = 64 to 73) during the last 60 days of gestation produced more milk than those cows in a hot environment (THI_g = 74 to 84), particularly during the first 150 days of lactation (DIM). The heritabilities (h2) of TD were similar throughout DIM for cows in THI_g hot (0.11 to 0.20) or (0.10 to 0.22), while the genetic correlations (rg) for TD between these two environments ranged from 0.11 to 0.52 along the first 250 DIM. The h2 estimates for TD across THI_m were similar for cows in THI_g hot (0.07 to 0.25) and THI_g cold (0.08 to 0.19). The rg estimates ranged from 0.17 to 0.42 along THI_m between TD of cows in cold and hot THI_g. The results were consistent in demonstrating the existence of an additional source of G × E for TD due to THI_g and THI_m. The present study is probably the first to provide evidence of this source of G × E; further research is needed because of its importance when the breeding objective is to select animals that are more tolerant to HS.

The dietary supplementation of canthaxanthin in combination with 25OHD3 results in reproductive, performance, and progeny quality gains in broiler breeders.

The objective of this study was to evaluate the effects of the dietary supplementation of canthaxanthin and cholecalciferol (25OHD3) to broiler breeders and their progenies on their performance. Eighty 25-wk-old Cobb 500 broiler breeders were distributed according to a completely randomized design in 2 experimental treatments, with 10 replicates of 4 birds each. The treatments consisted of the supplementation or not of a corn-soybean meal diet with canthaxanthin + 25OHD3 to broiler breeders fed from 25 to 62 wk of age. Egg production, fertility, hatchability, hatchability of fertile eggs, embryonic mortality, and egg specific gravity were evaluated. Breeders were inseminated at 35, 45, and 62 wk of age, the eggs were incubated, and the performance of the progenies was evaluated. From the progeny of each breeder age, 300 male chicks were distributed according to a completely randomized design in a 2 × 2 factorial arrangement (chicks from breeders supplemented or not with canthaxanthin + 25OHD3 and chicks supplemented or not with canthaxanthin + 25OHD3), totaling 4 treatments with 5 replicates of 15 birds each. Canthaxanthin + 25OHD3 were supplied to the chicks until 21 D of age. The combination of canthaxanthin and 25OHD3, containing 69 mg of 25OHD3 and 6 g of canthaxanthin, was supplemented at 1 kg/t of feed. Breeders supplemented with canthaxanthin + 25OHD3 showed higher egg production, total hatchability, hatchability of fertile eggs, and lower early embryo mortality compared with those fed the control diet. Broilers from breeders fed canthaxanthin + 25OHD3 and supplemented with this additive up to 21 D of age presented a better feed conversion ratio and higher carcass and breast yields than those derived from nonsupplemented breeders.

Genetic parameters for stayability to consecutive calvings in Zebu cattle.

Longer-lived cows tend to be more profitable and the stayability trait is a selection criterion correlated to longevity. An alternative to the traditional approach to evaluate stayability is its definition based on consecutive calvings, whose main advantage is the more accurate evaluation of young bulls. However, no study using this alternative approach has been conducted for Zebu breeds. Therefore, the objective of this study was to compare linear random regression models to fit stayability to consecutive calvings of Guzerá, Nelore and Tabapuã cows and to estimate genetic parameters for this trait in the respective breeds. Data up to the eighth calving were used. The models included the fixed effects of age at first calving and year-season of birth of the cow and the random effects of contemporary group, additive genetic, permanent environmental and residual. Random regressions were modeled by orthogonal Legendre polynomials of order 1 to 4 (2 to 5 coefficients) for contemporary group, additive genetic and permanent environmental effects. Using Deviance Information Criterion as the selection criterion, the model with 4 regression coefficients for each effect was the most adequate for the Nelore and Tabapuã breeds and the model with 5 coefficients is recommended for the Guzerá breed. For Guzerá, heritabilities ranged from 0.05 to 0.08, showing a quadratic trend with a peak between the fourth and sixth calving. For the Nelore and Tabapuã breeds, the estimates ranged from 0.03 to 0.07 and from 0.03 to 0.08, respectively, and increased with increasing calving number. The additive genetic correlations exhibited a similar trend among breeds and were higher for stayability between closer calvings. Even between more distant calvings (second v. eighth), stayability showed a moderate to high genetic correlation, which was 0.77, 0.57 and 0.79 for the Guzerá, Nelore and Tabapuã breeds, respectively. For Guzerá, when the models with 4 or 5 regression coefficients were compared, the rank correlations between predicted breeding values for the intercept were always higher than 0.99, indicating the possibility of practical application of the least parameterized model. In conclusion, the model with 4 random regression coefficients is recommended for the genetic evaluation of stayability to consecutive calvings in Zebu cattle.

Genetics of heat tolerance for milk yield and quality in Holsteins.

Tropical and sub-tropical climates are characterized by high temperature and humidity, during at least part of the year. Consequently, heat stress is common in Holstein cattle and productive and reproductive losses are frequent. Our objectives were as follows: (1) to quantify losses in production and quality of milk due to heat stress; (2) to estimate genetic correlations within and between milk yield (MY) and milk quality traits; and (3) to evaluate the trends of genetic components of tolerance to heat stress in multiple lactations of Brazilian Holstein cows. Thus, nine analyses using two-trait random regression animal models were carried out to estimate variance components and genetic parameters over temperature-humidity index (THI) values for MY and milk quality traits (three lactations: MY×fat percentage (F%), MY×protein percentage (P%) and MY×somatic cell score (SCS)) of Brazilian Holstein cattle. It was demonstrated that the effects of heat stress can be harmful for traits related to milk production and milk quality of Holstein cattle even though most herds were maintained in a modified environment, for example, with fans and sprinklers. For MY, the effect of heat stress was more detrimental in advanced lactations (-0.22 to -0.52 kg/day per increase of 1 THI unit). In general, the mean heritability estimates were higher for lower THI values and longer days in milk for all traits. In contrast, the heritability estimates for SCS increased with increasing THI values in the second and third lactation. For each trait studied, lower genetic correlations (different from unity) were observed between opposite extremes of THI (THI 47 v. THI 80) and in advanced lactations. The genetic correlations between MY and milk quality trait varied across the THI scale and lactations. The genotype×environment interaction due to heat stress was more important for MY and SCS, particularly in advanced lactations, and can affect the genetic relationship between MY and milk quality traits. Selection for higher MY, F% or P% may result in a poor response of the animals to heat stress, as a genetic antagonism was observed between the general production level and specific ability to respond to heat stress for these traits. Genetic trends confirm the adverse responses in the genetic components of heat stress over the years for milk production and quality. Consequently, the selection of Holstein cattle raised in modified environments in both tropical and sub-tropical regions should take into consideration the genetic variation in heat stress.

Inbreeding depression in Zebu cattle traits.

The productivity of herds may be negatively affected by inbreeding depression, and it is important to know how intense is this effect on the livestock performance. We performed a comprehensive analysis involving five Zebu breeds reared in Brazil to estimate inbreeding depression in productive and reproductive traits. Inbreeding depression was estimated for 13 traits by including the individual inbreeding rate as a linear covariate in the standard genetic evaluation models. For all breeds and for almost all traits (no effect was observed on gestation length), the performance of the animals was compromised by an increase in inbreeding. The average inbreeding depression was -0.222% and -0.859% per 1% of inbreeding for linear regression coefficients scaled on the percentage of mean (ßm ) and standard deviation (ßσ ), respectively. The means for ßm (and ßσ ) were -0.269% (-1.202%) for weight/growth traits and -0.174% (-0.546%) for reproductive traits. Hence, inbreeding depression is more pronounced in weight/growth traits than in reproductive traits. These findings highlight the need for the management of inbreeding in the respective breeding programmes of the breeds studied here.

Random regression analysis of test-day milk yields in the first and second lactations of Brazilian Gyr cows.

The objective of the present study was to estimate the genetic parameters for test-day milk yields (TDMY) in the first and second lactations using random regression models (RRM) in order to contribute to the application of these models in genetic evaluation of milk yield in Gyr cattle. A total of 53,328 TDMY records from 7118 lactations of 5853 Gyr cows were analyzed. The model included the direct additive, permanent environmental, and residual random effects. In addition, contemporary group and linear and quadratic effects of the age of cows at calving were included as fixed effects. A random regression model fitting fourth-order Legendre polynomials for additive genetic and permanent environmental effects, with five classes of residual variance, was applied. In the first lactation, the heritabilities increased from early lactation (0.26) until TDMY3 (0.38), followed by a decrease until the end of lactation. In the second lactation, the estimates increased from the first (0.29) to the fifth test day (0.36), with a slight decrease thereafter, and again increased on the last two test days (0.34 and 0.41). There were positive and high genetic correlations estimated between first-lactation TDMY and the remaining TDMY of the two lactations. The moderate heritability estimates, as well as the high genetic correlations between half the first-lactation TDMY and all TDMY of the two lactations, suggest that the selection based only on first lactation TDMY is the best selection strategy to increase milk production across first and second lactations of Gyr cows.

Annual changes in seminal variables of golden lanchead pitvipers (Bothrops insularis) maintained in captivity.

Bothrops insularis is an endemic and critically endangered snake with an estimated population of 2000 individuals restricted to Queimada Grande Island, in southeastern Brazil. Brazilian researchers established a captive breeding program for the species that includes the application of assisted reproductive technologies. The present study, therefore, aimed to evaluate semen samples from captive B. insularis throughout the year to ascertain seasonal differences in semen traits as well as correlations with body size and weight. Eighteen males with snout-vent length (SVL) ranging from 43.5 to 73.7 cm were collected at quarterly basis between August 2012 and May 2013. Macroscopic analysis revealed semen volumes ranging from 0.5 to 6.0 µL with samples featuring whitish to yellowish color and creamy and thick consistency. Viable sperm was obtained from all males indicating that individuals with SVL equal to or greater than 43.5 cm are sexually developed. However, adult and immature males (estimated by SVL) exhibited different seasonal profiles for motility and progressive motility. Adult males had a decrease in sperm motility and progressive motility during summer and spring, respectively, whereas the same variables did not vary throughout the year in immature snakes. Sperm concentration in all individuals was less (0.5 × 10(9) µL) during the winter, but no seasonal fluctuations were detected in semen volume. These findings are of particular importance to the development of reproductive tools such as male selection, artificial insemination and sperm freezing for the genetic management of this critically endangered snake.

Detrimental effect of selection for milk yield on genetic tolerance to heat stress in purebred Zebu cattle: Genetic parameters and trends.

In an attempt to determine the possible detrimental effects of continuous selection for milk yield on the genetic tolerance of Zebu cattle to heat stress, genetic parameters and trends of the response to heat stress for 86,950 test-day (TD) milk yield records from 14,670 first lactations of purebred dairy Gir cows were estimated. A random regression model with regression on days in milk (DIM) and temperature-humidity index (THI) values was applied to the data. The most detrimental effect of THI on milk yield was observed in the stage of lactation with higher milk production, DIM 61 to 120 (-0.099kg/d per THI). Although modest variations were observed for the THI scale, a reduction in additive genetic variance as well as in permanent environmental and residual variance was observed with increasing THI values. The heritability estimates showed a slight increase with increasing THI values for any DIM. The correlations between additive genetic effects across the THI scale showed that, for most of the THI values, genotype by environment interactions due to heat stress were less important for the ranking of bulls. However, for extreme THI values, this type of genotype by environment interaction may lead to an important error in selection. As a result of the selection for milk yield practiced in the dairy Gir population for 3 decades, the genetic trend of cumulative milk yield was significantly positive for production in both high (51.81kg/yr) and low THI values (78.48kg/yr). However, the difference between the breeding values of animals at high and low THI may be considered alarming (355kg in 2011). The genetic trends observed for the regression coefficients related to general production level (intercept of the reaction norm) and specific ability to respond to heat stress (slope of the reaction norm) indicate that the dairy Gir population is heading toward a higher production level at the expense of lower tolerance to heat stress. These trends reflect the genetic antagonism between production and tolerance to heat stress demonstrated by the negative genetic correlation between these components (-0.23). Monitoring trends of the genetic component of heat stress would be a reasonable measure to avoid deterioration in one of the main traits of Zebu cattle (i.e., high tolerance to heat stress). On the basis of current genetic trends, the need for future genetic evaluation of dairy Zebu animals for tolerance to heat stress cannot be ruled out.

Genetic parameters and investigation of genotype × environment interactions in Nellore × Hereford crossbred for resistance to cattle ticks in different regions of Brazil.

Data from 6,156 Nellore × Hereford crossbred cattle, distributed in 18 herds located in the Brazilian states of Mato Grosso do Sul (MS), São Paulo (SP), Paraná (PR) and Rio Grande do Sul (RS), were analysed in order to investigate genetic variation for resistance to the cattle tick through the estimation of variance components and genetic parameters for counting ticks (Rhipicephalus (Boophilus) microplus) in natural infestation, and also a possible genotype × environment interaction. The tick count data (CC) were transformed to log10(CC + 1) and grouped into two regions, defined by cluster analysis and analysed using the method of restricted maximum likelihood. The statistical model included the additive genetic effect as random and fixed effects of the contemporary group (CG) and genetic group (GG) as classificatory and the age of the animal at the time of counting (linear effect) and individual heterozygosity (linear effect, ranging from 0 to 1) as covariates. In the studied regions, the effect of heterosis found was inversely proportional; in other words, the count of ticks decreased with the increase of heterozygosity. The observed heritability estimation for resistance to tick infestation were 0.12 ± 0.04 and 0.11 ± 0.04 for groups A (RS and south PR) and B (MS, SP and north PR), respectively. The results of this study suggest that selection for animals resistant to ticks would be possible using the tick count to estimate the genetic value of animals, but with a slow genetic progress. The genetic correlation for tick count between the two groups (A and B) was 0.84 ± 0.27 and genotype × environment interaction for this trait was not observed in the studied population.

Utilisation of sperm-binding assay combined with computer-assisted sperm analysis to evaluate frozen-thawed bull semen.

Due to homologies between the chicken egg perivitelline membrane with mammalian zona pellucida proteins, spermatozoa of several species are able to bind to this membrane. However, adequate standardisation is required to attest possible applications of this technique for semen evaluation of a given species. Therefore, we thawed and divided cryopreserved semen samples into two aliquotes, one kept in water bath at 37 °C (thawed) and the other submitted to snap-freezing to damage sperm cells (dead spermatozoa). Aliquotes were mixed into different ratios of thawed:dead cells and analysed for motility, membrane and acrosomal integrity, and mitochondrial activity. In parallel, chicken egg perivitelline membranes were inseminated with these ratios, and the number of spermatozoa bound per mm(2) of membrane was assessed by conventional microscopy (CM) and computer-assisted sperm analysis (CASA). Linear regression showed high correlation between thawed:dead sperm ratio and number of spermatozoa bound to the membrane (CM: r(2) = 0.91 and CASA: r(2) = 0.92 respectively). Additionally, positive correlations were found between the number of spermatozoa bound to the membrane and acrosomal integrity, membrane integrity, mitochondrial activity and motility. These findings indicate that sperm-egg-binding assay associated with CASA is a reliable, practical and inexpensive method for examining the fertilising capacity of cryopreserved bull semen.

Genotype by environment interaction for post-weaning weight gain, scrotal circumference, and muscling score of composite beef cattle in different regions of Brazil.

The objectives of this study were to characterize and define homogenous production environments of composite beef cattle in Brazil in terms of climatic and geographic variables by using multivariate exploratory techniques; to evaluate the presence of genotype by environment interaction (GxE) for post-weaning weight gain (PWG), yearling scrotal circumference (SC), and yearling muscling (MUS). Hierarchical and nonhierarchical cluster analysis was used to group farms located in regions with similar environmental variables into clusters. Six clusters of farms were formed. The effect of sire-cluster interaction was tested by single-trait analysis. Genetic parameters were estimated by multi-trait analysis considering the same trait to be different in each cluster. The effect of sire-cluster interaction was significant (P <0.01) for PWG and MUS. Estimates of genetic correlations among clusters ranged from 0.31 to 0.93 for PWG, 0.64 to 0.89 for SC, and 0.18 to 0.80 for MUS. These results indicate the need for a genetic analysis on a regional basis or inclusion of the GxE effect in the statistical model to permit appropriate evaluation of the animals.

Linear and Poisson models for genetic evaluation of tick resistance in cross-bred Hereford x Nellore cattle.

Cattle resistance to ticks is measured by the number of ticks infesting the animal. The model used for the genetic analysis of cattle resistance to ticks frequently requires logarithmic transformation of the observations. The objective of this study was to evaluate the predictive ability and goodness of fit of different models for the analysis of this trait in cross-bred Hereford x Nellore cattle. Three models were tested: a linear model using logarithmic transformation of the observations (MLOG); a linear model without transformation of the observations (MLIN); and a generalized linear Poisson model with residual term (MPOI). All models included the classificatory effects of contemporary group and genetic group and the covariates age of animal at the time of recording and individual heterozygosis, as well as additive genetic effects as random effects. Heritability estimates were 0.08 ± 0.02, 0.10 ± 0.02 and 0.14 ± 0.04 for MLIN, MLOG and MPOI models, respectively. The model fit quality, verified by deviance information criterion (DIC) and residual mean square, indicated fit superiority of MPOI model. The predictive ability of the models was compared by validation test in independent sample. The MPOI model was slightly superior in terms of goodness of fit and predictive ability, whereas the correlations between observed and predicted tick counts were practically the same for all models. A higher rank correlation between breeding values was observed between models MLOG and MPOI. Poisson model can be used for the selection of tick-resistant animals.

Comparison of factor-analytic and reduced rank models for test-day milk yield in Gyr dairy cattle (Bos indicus).

We analyzed 46,161 monthly test-day records of milk production from 7453 first lactations of crossbred dairy Gyr (Bos indicus) x Holstein cows. The following seven models were compared: standard multivariate model (M10), three reduced rank models fitting the first 2, 3, or 4 genetic principal components, and three models considering a 2-, 3-, or 4-factor structure for the genetic covariance matrix. Full rank residual covariance matrices were considered for all models. The model fitting the first two principal components (PC2) was the best according to the model selection criteria. Similar phenotypic, genetic, and residual variances were obtained with models M10 and PC2. The heritability estimates ranged from 0.14 to 0.21 and from 0.13 to 0.21 for models M10 and PC2, respectively. The genetic correlations obtained with model PC2 were slightly higher than those estimated with model M10. PC2 markedly reduced the number of parameters estimated and the time spent to reach convergence. We concluded that two principal components are sufficient to model the structure of genetic covariances between test-day milk yields.

Random regression models using Legendre polynomials or linear splines for test-day milk yield of dairy Gyr (Bos indicus) cattle.

Studies investigating the use of random regression models for genetic evaluation of milk production in Zebu cattle are scarce. In this study, 59,744 test-day milk yield records from 7,810 first lactations of purebred dairy Gyr (Bos indicus) and crossbred (dairy Gyr × Holstein) cows were used to compare random regression models in which additive genetic and permanent environmental effects were modeled using orthogonal Legendre polynomials or linear spline functions. Residual variances were modeled considering 1, 5, or 10 classes of days in milk. Five classes fitted the changes in residual variances over the lactation adequately and were used for model comparison. The model that fitted linear spline functions with 6 knots provided the lowest sum of residual variances across lactation. On the other hand, according to the deviance information criterion (DIC) and bayesian information criterion (BIC), a model using third-order and fourth-order Legendre polynomials for additive genetic and permanent environmental effects, respectively, provided the best fit. However, the high rank correlation (0.998) between this model and that applying third-order Legendre polynomials for additive genetic and permanent environmental effects, indicates that, in practice, the same bulls would be selected by both models. The last model, which is less parameterized, is a parsimonious option for fitting dairy Gyr breed test-day milk yield records.