Dual-purpose Guzerá cattle exhibit high dairy performance under heat stress.

The present study evaluated the heat stress response pattern of dual-purpose Guzerá cattle for test-day (TD) milk yield records of first lactation and estimated genetic parameters and trends related to heat stress. A total of 31,435 TD records from 4,486 first lactations of Guzerá cows, collected between 1986 and 2012, were analysed. Two random regression models considered days in milk (DIM) and/or temperature × humidity-dependent (THI) covariate. Impacts of -0.037, -0.019 and -0.006 kg/day/THI for initial and intermediate stages of lactation were observed when considering the mean maximum daily temperature and humidity to calculate THI. Heritability estimates ranged from 0.16 to 0.35 throughout lactation and THI values, suggesting the possibility to expect gains from selection for such trait. The variable trajectory of breeding values for dual-purpose Guzerá sires in response to changes in THI values confirms that the genotype × environment interaction due to heat stress can have some effect on TD milk yield. Despite the high dairy performance of Guzerá cattle under heat stress, estimated genetic trends showed a progressive reduction in heat tolerance. Therefore, new strategies should be adopted to prevent negative impacts of heat stress over milk production in Guzerá animals in future.

Reducing cryptic relatedness in genomic data sets via a central node exclusion algorithm.

Cryptic relatedness is a confounding factor in genetic diversity and genetic association studies. Development of strategies to reduce cryptic relatedness in a sample is a crucial step for downstream genetic analyses. This study uses a node selection algorithm, based on network degrees of centrality, to evaluate its applicability and impact on evaluation of genetic diversity and population stratification. 1,036 Guzerá (Bos indicus) females were genotyped using Illumina Bovine SNP50 v2 BeadChip. Four strategies were compared. The first and second strategies consist on a iterative exclusion of most related individuals based on PLINK kinship coefficient (φij) and VanRaden's φij, respectively. The third and fourth strategies were based on a node selection algorithm. The fourth strategy, Network G matrix, preserved the larger number of individuals with a better diversity and representation from the initial sample. Determining the most probable number of populations was directly affected by the kinship metric. Network G matrix was the better strategy for reducing relatedness due to producing a larger sample, with more distant individuals, a more similar distribution when compared with the full data set in the MDS plots and keeping a better representation of the population structure. Resampling strategies using VanRaden's φij as a relationship metric was better to infer the relationships among individuals. Moreover, the resampling strategies directly impact the genomic inflation values in genomewide association studies. The use of the node selection algorithm also implies better selection of the most central individuals to be removed, providing a more representative sample.