Genome-wide association and functional genomic analyses for various hoof health traits in North American Holstein cattle.

Hoof diseases are a major welfare and economic issue in the global dairy cattle production industry, which can be minimized through improved management and breeding practices. Optimal genetic improvement of hoof health could benefit from a deep understanding of the genetic background and biological underpinning of indicators of hoof health. Therefore, the primary objectives of this study were to perform genome-wide association studies, using imputed high-density genetic markers data from North American Holstein cattle, for 8 hoof-related traits: digital dermatitis, sole ulcer, sole hemorrhage, white line lesion, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, and toe ulcer, and a hoof health index. De-regressed estimated breeding values from 25,580 Holstein animals were used as pseudo-phenotypes for the association analyses. The genomic quality control, genotype phasing, and genotype imputation were performed using the PLINK (version 1.9), Eagle (version 2.4.1), and Minimac4 software, respectively. The functional genomic analyses were performed using the GALLO R package and the DAVID platform. We identified 22, 34, 14, 22, 28, 33, 24, 43, and 15 significant markers for digital dermatitis, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, sole hemorrhage, sole ulcer, toe ulcer, white line lesion disease, and the hoof health index, respectively. The significant markers were located across all autosomes, except BTA10, BTA12, BTA20, BTA26, BTA27, and BTA28. Moreover, the genomic regions identified overlap with various previously reported quantitative trait loci for exterior, health, meat and carcass, milk, production, and reproduction traits. The enrichment analyses identified 44 significant gene ontology terms. These enriched genomic regions harbor various candidate genes previously associated with bone development, metabolism, and infectious and immunological diseases. These findings indicate that hoof health traits are highly polygenic and influenced by a wide range of biological processes.

A systematic review with meta-analysis of heritability estimates for temperament-related traits in beef and dairy cattle populations.

Temperament (docility) is a key breeding goal in the cattle industry due to its direct relationship with animal welfare, cattle handler's safety and animal productivity. Over the past six decades, numerous studies have reported heritability estimates for temperament-related traits in cattle populations ranging from low to high values. Therefore, the primary objective of this study was to perform a comprehensive systematic review with meta-analysis to obtain weighted estimates of heritability for temperament-related traits in worldwide cattle populations. After data editing and quality control, 106 studies were included in the systematic review, of which 29.2% and 70.8% reported estimates of heritability for temperament-related traits in dairy and beef cattle populations, respectively. Meta-analyses were performed for 95 heritability estimates using a random model approach. The weighted heritability estimates were as follow: (a) flight score at weaning = 0.23 (95% CI: 0.15-0.32); (b) flight speed at weaning = 0.30 (95% CI: 0.26-0.33); (c) joint analysis of flight speed and flight score at weaning = 0.27 (95% CI: 0.22-0.31); (d) flight speed at yearling = 0.26 (95% CI: 0.21-0.30); (e) joint analysis of flight speed at weaning and yearling = 0.27 (95% CI: 0.24-0.30); (f) movement score = 0.12 (95% CI: 0.08-0.15); (g) crush score at weaning = 0.21 (95% CI: 0.17-0.25); (h) pen score at weaning = 0.27 (95% CI: 0.19-0.34); (i) pen score at yearling = 0.20 (95% CI: 0.17-0.23); (j) joint analysis of pen score at weaning and yearling = 0.22 (95% CI: 0.18-0.26); (k) cow's aggressiveness at calving = 0.10 (95% CI: 0.01-0.19); (l) general temperament = 0.13 (95% CI: 0.06-0.19); (m) milking temperament = 0.16 (95% CI: 0.11-0.21); and (n) joint analysis of general and milking temperament = 0.14 (95% CI: 0.11-0.18). The heterogeneity index ranged from 0% to 77%, and the Q-test was significant (p < 0.05) for four single-trait meta-analyses. In conclusion, temperament is moderately heritable in beef cattle populations, and flight speed at weaning had the highest weighted heritability estimate. Moreover, between-study heterogeneity was low or moderate in beef cattle traits, suggesting reasonable standardization across studies. On the other hand, low-weighted heritability and high between-study heterogeneity were estimated for temperament-related traits in dairy cattle, suggesting that more studies are needed to better understand the genetic inheritance of temperament in dairy cattle populations.

Characterization of runs of homozygosity, heterozygosity-enriched regions, and population structure in cattle populations selected for different breeding goals.

BACKGROUND: A decline in the level of genetic diversity in livestock can result in reduced response to selection, greater incidence of genetic defects, and inbreeding depression. In this context, various metrics have been proposed to assess the level of genetic diversity in selected populations. Therefore, the main goals of this study were to: 1) investigate the population structure of 16 cattle populations from 15 different pure breeds or composite populations, which have been selected for different breeds goals; and, 2) identify and compare runs of homozygosity (ROH) and heterozygosity-enriched regions (HER) based on different single nucleotide polymorphism (SNP) panels and whole-genome sequence data (WGS), followed by functional genomic analyses. RESULTS: A total of 24,187 ROH were found across all cattle populations, with 55% classified in the 2-4 Mb size group. Fourteen homozygosity islands were found in five populations, where four ROH islands located on BTA1, BTA5, BTA16, and BTA19 overlapped between the Brahman (BRM) and Gyr (GIR) breeds. A functional analysis of the genes found in these islands revealed candidate genes known to play a role in the melanogenesis, prolactin signaling, and calcium signaling pathways. The correlations between inbreeding metrics ranged from 0.02 to 0.95, where the methods based on homozygous genotypes (FHOM), uniting of gametes (FUNI), and genotype additive variance (FGRM) showed strong correlations among them. All methods yielded low to moderate correlations with the inbreeding coefficients based on runs of homozygosity (FROH). For the HER, 3576 runs and 26 islands, distributed across all autosomal chromosomes, were found in regions containing genes mainly related to the immune system, indicating potential balancing selection. Although the analyses with WGS did not enable detection of the same island patterns, it unraveled novel regions not captured when using SNP panel data. CONCLUSIONS: The cattle populations that showed the largest amount of ROH and HER were Senepol (SEN) and Montana (MON), respectively. Overlapping ROH islands were identified between GIR and BRM breeds, indicating a possible historical connection between the populations. The distribution and pattern of ROH and HER are population specific, indicating that different breeds have experienced divergent selection processes or different genetic processes.

Variants in myostatin and MyoD family genes are associated with meat quality traits in Santa Inês sheep.

Myostatin and MyoD family genes play vital roles in myogenesis and this study aimed to identify association of variants in MyoD1, MyoG, MyF5, MyF6, and MSTN genes with meat traits in Santa Inês sheep. A dataset with 44 variants and records of seven meat traits in 192 lambs (pH0, pH24, a*, b*, L*, tenderness assessed by shear force, and water-holding capacity) was used. Single-locus and haplotype association analyses were performed, and the significance threshold was established according to Bonferroni's method. Single-locus analysis revealed two associations at a Bonferroni level, where the variant c.935-185C > G in MyoD1 had an additive effect (-4.31 ± 1.08 N) on tenderness, while the variant c.464 + 185G > A in MyoG had an additive effect (-2.86 ± 0.64) on a*. Additionally, the haplotype replacement GT>AC in MSTN was associated with pH0 (1.26 ± 0.31), pH24 (1.07 ± 0.27), a* (-1.40 ± 0.51), and tenderness (3.83 ± 1.22 N), while the replacement GT > AG in MyoD1 was associated with pH0 (1.43 ± 0.26), pH24 (1.25 ± 0.22), b* (-1.06 ± 0.39), and tenderness (-4.13 ± 1.16 N). Our results have demonstrated that some variants in MyoG, MyF6, MyoD1, and MSTN can be associated with physicochemical meat traits in Santa Inês sheep.

CNV detection and their association with growth, efficiency and carcass traits in Santa Inês sheep.

Copy number variations (CNV) are an important source of genetic variation. CNV has been increasingly studied and frequently associated with diseases and productive traits in livestock animals. However, CNV-based genome-wide association studies (GWAS) in Santa Inês sheep, one of the principal sheep breeds in Brazil, have not yet been reported. Thus, the aim of this study was to investigate the association between CNV and growth, efficiency and carcass traits in sheep. The Illumina OvineSNP50 BeadChip array was used to detect CNV in 491 Santa Inês individuals. Then, CNV-based GWAS was performed with a linear mixed model approach considering a genomic relationship matrix, for ten traits: (1) growth: body weight at three (W3) and six (W6) months of age; (2) efficiency: residual feed intake (RFI) and feed efficiency (FE) and (3) carcass: external carcass length (ECL), leg length (LL), carcass yield (CY), commercial cuts weight (CCW), loin eye area (LEA) and subcutaneous fat thickness (SFT). We identified 1,167 autosomal CNV in 438 sheep, with 294 non-redundant CNV, ranging from 21.8 to 861.9 kb, merged into 216 distinct copy number variation regions (CNVRs). One significant CNV segment (pFDR -value<0.05) in OAR3 was associated with CY, while another significant CNV in OAR6 was associated with RFI. Additionally, another 5 CNV segments were considered relevant for investigation in the future studies. The significant segments overlapped 4 QTLs and spanned 8 genes, including the SPAST,TGFA and ADGRL3 genes, involved in cell differentiation and energy metabolism. Therefore, the results of the present study increase knowledge about CNV in sheep, their possible impacts on productive traits, and provide information for future investigations, being especially useful for those interested in structural variations in the sheep genome.

Genetic origin of donkeys in Brazil.

INTRODUCTION: The genetic groups of native donkeys in Brazil are characterized by adaptation to the local environment. However, the donkey population in the country is declining, mainly because of agricultural mechanization and transportation that has led to the abandonment and the consequent indiscriminate slaughter of these animals. There are three local genetic groups of distinct geographic and temporal formation. However, analyses of their origin, phylogenetic relationship, and population structure are scarce. Within this context, molecular markers such as the mitochondrial control region (D-loop) are useful for these analyses. MATERIAL AND METHODS: This study aimed to evaluate the variation and origin of maternal lineages of groups of naturalized donkeys in Brazil (Brazilian, Nordestino, and Pêga). RESULTS AND DISCUSSION: We detected five mitochondrial haplotypes with 19 polymorphic sites, two of them exclusively found in the Nordestino donkey; this group is in fact more distant from the others. Phylogenetic analysis indicates maternal contributions of two clades (Nubian and Somali) to the formation of the genetic groups of donkeys, a fact that explains the high diversity, structure, and distances of the groups, reported here for the first time. CONCLUSION : This analysis contributes production and conservation of native donkey breeds. It also gives clues about the formation of the Iberian breeds from which Brazilian donkeys originated.

Flaxseed added to the diet of Alpine goats affects the nutrients intake and blood parameters.

The flaxseed is a nutraceutical food used as a source of α-linolenic acid, which can bring benefits to the health of mammals. This study was carried out to examine the effect of flaxseed inclusion in the diets on the intake of nutrients, body weight, and blood parameters of Alpine goats. Twenty-one adult females with an initial average weight of 41.06 ± 1.84 kg were used in a completely randomized design, with four experimental treatments (0, 5, 10, and 15% of flaxseed in the total diet) and five replications per treatment. The intake of the dry matter presented a decreasing linear effect (P < 0.001), with a reduction of 53.5% between the control diet (0% of flaxseed) and 15% of flaxseed, but no effect was observed on weight gain (P > 0.05). Inclusion of flaxseed from 0 to 15% linearly decreased the intakes of organic matter, crude protein, non-fibrous carbohydrates, and neutral detergent fiber, but increased the ether extract intake (P < 0.001). Regarding plasma concentration traits, increasing the flaxseed levels from 0 to 15% had linear positive effect on LDL, VLDL, and triglycerides (P < 0.05), but no effect on cholesterol, HDL, glucose, creatinine, and urea (P > 0.05). The inclusion of flaxseed in the diet for 80 days changes the nutrients intake and blood parameters but shows no impact on body weight. However, further studies are required to determine the impact of flaxseed on the goat's health in the long term since there has been an increase in the concentration of triglycerides, LDL, and VLDL.

Genotype by environment interaction for fat and protein yields via reaction norms in Holstein cattle of southern Brazil.

Our objective was to evaluate the genetic merit of Holstein cattle population in southern Brazil in response to variations in the regional temperature by analyzing the genotype by environment interaction using reaction norms. Fat yield (FY) and protein yield (PY) data of 67 360 primiparous cows were obtained from the database of the Paraná Holstein Breeders Association, Brazil (APCBRH). The regional average annual temperature was used as the environmental variable. A random regression model was adopted applying mixed models with Restricted Maximum Likelihood (REML) algorithm using WOMBAT software. The genetic merit of the 15 most representative bulls, depending on the temperature gradient, was evaluated. Heritability ranged from 0.21 to 0.27 for FY and from 0.14 to 0.20 for PY. The genetic correlation observed among the environmental gradients proved to be higher than 0.80 for both traits. Slight reranking of bulls for both traits was detected, demonstrating that non-relevant genotype by environment interaction for FY and PY were observed. Consequently, no inclusion of the temperature effect in the model of genetic evaluation in southern Brazilian Holstein breed is required.

Genome-wide association study and pathway analysis for fat deposition traits in nellore cattle raised in pasture-based systems.

Genome-wide association study (GWAS) is a powerful tool to identify candidate genes and genomic regions underlying key biological mechanisms associated with economically important traits. In this context, the aim of this study was to identify genomic regions and metabolic pathways associated with backfat thickness (BFT) and rump fat thickness (RFT) in Nellore cattle, raised in pasture-based systems. Ultrasound-based measurements of BFT and RFT (adjusted to 18 months of age) were collected in 11,750 animals, with 39,903 animals in the pedigree file. Additionally, 1,440 animals were genotyped using the GGP-indicus 35K SNP chip, containing 33,623 SNPs after the quality control. The single-step GWAS analyses were performed using the BLUPF90 family programs. Candidate genes were identified through the Ensembl database incorporated in the BioMart tool, while PANTHER and REVIGO were used to identify the key metabolic pathways and gene networks. A total of 18 genomic regions located on 10 different chromosomes and harbouring 23 candidate genes were identified for BFT. For RFT, 22 genomic regions were found on 14 chromosomes, with a total of 29 candidate genes identified. The results of the pathway analyses showed important genes for BFT, including TBL1XR1, AHCYL2, SLC4A7, AADAT, VPS53, IDH2 and ETS1, which are involved in lipid metabolism, synthesis of cellular amino acids, transport of solutes, transport between Golgi Complex membranes, cell differentiation and cellular development. The main genes identified for RFT were GSK3ß, LRP1B, EXT1, GRB2, SORCS1 and SLMAP, which are involved in metabolic pathways such as glycogen synthesis, lipid transport and homeostasis, polysaccharide and carbohydrate metabolism. Polymorphisms located in these candidate genes can be incorporated in commercial genotyping platforms to improve the accuracy of imputation and genomic evaluations for carcass fatness. In addition to uncovering biological mechanisms associated with carcass quality, the key gene pathways identified can also be incorporated in biology-driven genomic prediction methods.

Water buffaloes (Bubalus bubalis) only have A2A2 genotype for beta-casein.

Beta-casein is a milk protein that has two variants: A1 and A2. Some individuals have difficulties digesting beta-casein A1, which can cause gastrointestinal disorders. A2 milk has emerged as an alternative. This milk only contains beta-casein A2 and is obtained from females carrying the A2A2 genotype of the gene. In cattle, allele and genotype frequencies vary according to breed and marker-assisted selection is performed to obtain A2A2 animals and the consequent production of A2 milk that is easier to digest. This study aimed to evaluate the alleles of beta-casein in buffaloes. A total of 657 buffaloes of four different breeds were genotyped and all animals carried the A2A2 genotype, i.e., allele A1 does not exist in the buffalo species. Thus, all milk products of buffaloes are naturally A2. This result adds value to products derived from buffalo milk.

Evaluation of genotype by environment interactions on milk production traits of Holstein cows in southern Brazil.

OBJECTIVE: This study assessed the possible existence of genotype by environment interactions for milk, fat and protein yields in Holstein cattle raised in one of the most important milk production basins in Brazil. METHODS: Changes in the genetic parameters and breeding values were evaluated for 57,967 animals from three distinct regions of southern Brazil, divided according to differences in climate. The genotype by environment interaction was determined by genetic correlations between regions, estimated by the restricted maximum likelihood, considering the animal model. Bull rankings were investigated to verify the ratio of coincident selected animals between regions for each trait. RESULTS: The estimates of heritability coefficients were similar between two regions, but were lower in the third evaluated area, for all traits. Genetic correlations between regions were high, ranging from 0.91 to 0.99 for milk, fat and protein yields, representing the absence of a genotype by environment interaction for productive traits. The percentage of selection error between regions for the top 10% of animals ranged from 0.88% to 2.07% for milk yield, 0.99% to 2.46% for fat yield and 0.59% to 3.15% for protein yield. CONCLUSION: A slight change in genotype between areas was expected since no significant genotype by environment interactions were identified, facilitating the process of selecting Holstein cattle in southern Brazil.

Effects of birth type and family on the variation of carcass and meat traits in Santa Ines sheep.

This study aimed to evaluate the effects of birth type and family on quantitative and qualitative meat traits in Santa Ines lambs and to do a descriptive study of variation inside of breed. A total of 50 carcass and meat traits were evaluated in up to 101 lambs slaughtered at 424 days of age. Random effect of family, the fixed effect of birth type (single or double) and year, and the covariates of dam's weight at weaning and age (in days) at slaughter were inserted in the model. This study shows that some carcass and meat traits in Santa Ines raised on pasture has large variability, specially for meat cut yields, several carcass non-components yields, and rib eye area, subcutaneous fat thickness, shear force, cooking losses, and b* color parameter. The birth type effect influenced (P < 0.05) rump width; reticulum, leg, heart and kidney yields; and the carcass pH at 24 h after slaughter. The family effect was significant (P < 0.05) for cold carcass weight; hot, cold and true carcass yields; external carcass and leg lengths; perimeter and width of the rump; kidney, liver and tail yields; and carcass pH at slaughter.

Performance of the Santa Ines breed raised on pasture in semiarid tropical regions and factors that explain trait variation.

This study aimed to evaluate sex, the number of lambs per birth, and the family effects on production traits in the Santa Ines breed of sheep by estimating the least square means and coefficient of variance for those traits. A total of 484 lambs were evaluated for the following traits: weight at birth, at weaning, and at 240 days of age; weight gain during the pre-weaning and post-weaning periods; height, width, and length of different body regions; and rib eye area and fat thickness between the 12th and 13th ribs. We observed coefficients of variation higher than 10 % for several traits. Generally, males were larger than females (P < 0.05), while lambs from single and double births were larger than lambs from triple births (P < 0.05). Family effect was significant (P < 0.05) for most traits and explained the highest percentage of residual variance. The results showed good development of Santa Ines sheep, especially during the pre-weaning period but no in post-weaning. Our study also showed that there is an effect of sex, birth type, and family, which must be included in any statistical model for the estimation of least square means and residual variance in ANOVA.

Genome-wide association and functional annotation analysis for the calving interval in Nellore cattle.

Calving interval (CI) measures the number of days between two consecutive calves of the same cow, and previous studies based on phenotype and pedigree data reported low heritability for this trait. However, the genetic architecture of CI in the Nellore breed was not evaluated based on genomic data. Thus, this study aimed to estimate the heritability based on genomic data and carry out a genome-wide association study (GWAS) for CI in the Nellore breed, using 12,599 pedigree records, 5078 CI records, and 3818 animals genotyped with 50k SNPchip panel. Both quality control and GWAS were performed in BLUPF90 family packages, which use the single-step genomic best linear unbiased predictor (ssGBLUP) method. The average CI was 427.6 days, with a standard deviation of 106.9 and a total range of 270-730 days. The heritability estimate was 0.04 ± 0.04. The p-values of GWAS analysis resulted in a genomic inflation factor (lambda) of 1.08. The only significant SNP (rs136725686) at the genome-wide level (p-value = 1.53E-06) was located on BTA13. Other 19 SNPs were significant at the chromosome-wide level, distributed on BTA1, 2, 3, 6, 10, 13, 14, 17, 18, 22, and 26. Functional annotation analysis found thirty-six protein-coding genes, including genes related to cell cycle (RAD21, BCAR3), oocyte function (LHX8, CLPX, UTP23), immune system (TXK, TEC, NFATC2), endocrine function (LRRFIP2, GPR158), estrous cycle (SLC38A7), and female fertility (CCK, LYZL4, TRAK1, FOXP1, STAC). Therefore, CI is a complex trait with small heritability in Nellore cattle, and various biological processes may be involved with the genetic architecture of CI in Nellore cattle.

Live weight, carcass ultrasound images, and visual scores in Angus cattle under feeding regimes in Brazil.

This study aimed to evaluate visual precocity, muscling, conformation, skeletal, and breed scores; live weights at birth, at 205, and at 550 days of age; and, besides, rib eye area and fat thickness between the 12th and 13th ribs obtained by ultrasound. Those traits were evaluated in 1,645 Angus cattle kept in five feeding conditions as follows: supplemented or non-supplemented, grazing native pasture or grazing cultivated pasture, and feedlot. Descriptive statistics, Pearson's correlations, and principal component analysis were carried out. Gender and feeding conditions were fixed effects, while animal's age and mother's weight at weaning were the covariates analyzed. Gender and feeding conditions were very significant for the studied traits, but visual scores were not influenced by gender. Animal's age and mother's weight at weaning influenced many traits and must be appropriately adjusted in the statistical models. An important correlation between visual scores, live weights, and carcass traits obtained by ultrasound was found, which can be analyzed by univariate procedure. However, the multivariate approach revealed some information that cannot be neglected in order to ensure a more detailed assessment.

A comprehensive characterization of longevity and culling reasons in Canadian Holstein cattle based on various systematic factors.

The decision of premature culling cows directly impacts the profitability of dairy farms. A comprehensive characterization of the primary causes of culling reasons would greatly improve both management and selection objectives in dairy cattle breeding programs. Therefore, this study aimed to analyze the temporal frequencies of 34 culling reasons in Canadian Holstein cows. After data editing and quality control, records from 3,096,872 cows culled from 9,683 herds spread across Canada were used for the analyses covering the periods from 1996 to 2020. Reproductive issues were the main culling reason accounting for 23.02%, followed by milk production (20.82%), health (20.39%), conformation problems (13.69%), economic factors (13.10%), accidents (5.67%), age-related causes (1.67%), and workability (1.63%). Nearly fifty-eight percent of cows were culled after 47 months of age. The observed frequencies of culling due to economic factors were lower than expected from 1996 to 2014 and higher than expected between 2015 and 2020. Reproduction issues had the highest culling frequencies during fall (24.54%), winter (24.02%), and spring (22.51%), while health issues were the most frequent (22.51%) culling reason in the summer season. Health issues (25.50%) and milk production (27.71%) were the most frequent culling reasons in the provinces of Quebec and Ontario, respectively. Reproductive issues showed the highest frequency across climates based on the Köppen climate classification, except for Csb (Dry-summer subtropical or Mediterranean climate) and Bsk (Middle latitude steppe climate), which correspond to small regions in Canada, where production was the most frequent culling reason (29.42% and 21.56%, respectively). Reproductive and milk performance issues were the two main culling reasons in most ecozones, except in Boreal Shield and Atlantic Marine, where health issues had the highest frequencies (25.12 and 23.75%, respectively). These results will contribute to improving management practices and selective decisions to reduce involuntary culling of Holstein cows.

Genetic and Genomic Characterization of a New Beef Cattle Composite Breed (Purunã) Developed for Production in Pasture-Based Systems.

Purunã is a composite beef cattle breed, developed in Southern Brazil by crossing the Angus, Charolais, Canchim, and Caracu breeds. The goal of this study was to perform the first genetic characterization of the Purunã breed, based on both pedigree and genomic information. For this, 100 randomly selected animals were genotyped, and 11,205 animals born from 1997 to 2019 had pedigree information. The genetic analyses performed were principal component analysis, admixture, phylogenic tree, pedigree and genomic inbreeding, linkage disequilibrium (LD), effective population size (Ne), consistency of the gametic phase, runs of homozygosity (ROH), heterozygosity-enriched regions (HERs), and functional analyses of the ROH and HER regions identified. Our findings indicate that Purunã is more genetically related to the Charolais, Canchim, and Angus breeds than Caracu or Nellore. The levels of inbreeding were shown to be small based on all the metrics evaluated and ranged from -0.009 to 0.029. A low (-0.12-0.31) correlation of the pedigree-based inbreeding compared to all the genomic inbreeding coefficients evaluated was observed. The LD average was 0.031 (±0.0517), and the consistency of the gametic phase was shown to be low for all the breed pairs, ranging from 0.42 to 0.27 to the distance of 20 Mb. The Ne values based on pedigree and genomic information were 158 and 115, respectively. A total of 1,839 ROHs were found, and the majority of them are of small length (<4 Mb). An important homozygous region was identified on BTA5 with pathways related to behavioral traits (sensory perception, detection of stimulus, and others), as well as candidate genes related to heat tolerance (MY O 1A), feed conversion rate (RDH5), and reproduction (AMDHD1). A total of 1,799 HERs were identified in the Purunã breed with 92.3% of them classified within the 0.5-1 Mb length group, and 19 HER islands were identified in the autosomal genome. These HER islands harbor genes involved in growth pathways, carcass weight (SDCBP), meat and carcass quality (MT2A), and marbling deposition (CISH). Despite the genetic relationship between Purunã and the founder breeds, a multi-breed genomic evaluation is likely not feasible due to their population structure and low consistency of the gametic phase among them.