A small proportion of Zebu genetic background in crossbred calves may not be enough to improve resistance against natural bovine Babesia spp. infections.

The main objective of cattle breeders in tropical and subtropical regions is to acquire animals with taurine-productive traits adapted to the broad weather range of these regions. However, one of the main challenges on using taurine genetics in these areas is the high susceptibility of these animals to tick-borne diseases. Consequently, the present study evaluated from 10 November 2021-19 April 2022, the over 13 assessments, the Babesia bovis and Babesia bigemina DNA loads and the IgG anti-B. bovis and anti-B. bigemina levels in Angus (n = 17, 100% Taurine) and Ultrablack (n = 14, ∼82% taurine and 18% Zebu) calves. Data were analyzed using a multivariate mixed model with repeated measures of the same animal including the fixed effects of evaluation, genetic group, sex, Babesia spp., and their interactions. The repeatability values were estimated from the (co)variances matrix and expressed for each species. The correlations between the DNA loads (CNlog) and IgG titers (S/P) values for the two species were also estimated using the same model. Regarding the specific IgG antibody titers for both Babesia spp., no significant differences were observed between the two genetic groups. However, for B. bovis and B. bigemina DNA loads, Ultrablack calves presented significantly higher values than Angus calves. Under the conditions evaluated in this study, our findings suggest that the low percentage of Zebu genetic in the Ultrablack breed was insufficient to improve resistance against babesiosis. Further studies must demonstrate if the low percentages of Zebu genetics in Taurine breeds can modify the susceptibility to babesiosis infections.

Estimation of genetic parameters for the tick and hemoparasite burden in Angus cattle.

The study was conducted with the objective of estimating genetic and phenotypic parameters for tick (CRM) and Babesia bigemina (IBBi), Babesia bovis (IBBo), and Anaplasma marginale (IAM) burden in Angus female breed in Brazil. The sample group was composed of Angus females raised in herds located in a region of endemic instability for cattle tick fever in the state of Rio Grande Sul (RS), Brazil. The variance components were estimated using Bayesian inference and Gibbs sampling algorithm, considering a multi-trait animal model. Heritability estimates showed values of low magnitude, ranging from 0.03 (IBBo) to 0.16 (CRM), while repeatability estimates ranged between 0.07 (IBBo) and 0.21 (CRM). Regarding the genetic correlation estimates, the values showed low (-0.01 for IBBo × IAM) to moderate (0.55 between IBBi × IAM) magnitudes. The results indicate that it is possible to use tick count and hemoparasite infection levels as selection criteria, with small genetic gains.

Identification of candidate lethal haplotypes and genomic association with post-natal mortality and reproductive traits in Nellore cattle.

The wide use of genomic information has enabled the identification of lethal recessive alleles that are the major genetic causes of reduced conception rates, longer calving intervals, or lower survival for live-born animals. This study was carried out to screen the Nellore cattle genome for lethal recessive haplotypes based on deviation from the expected population homozygosity, and to test SNP markers surrounding the lethal haplotypes region for association with heifer rebreeding (HR), post-natal mortality (PNM) and stayability (STAY). This approach requires genotypes only from apparently normal individuals and not from affected embryos. A total of 62,022 animals were genotyped and imputed to a high-density panel (777,962 SNP markers). Expected numbers of homozygous individuals were calculated, and the probabilities of observing 0 homozygotes was obtained. Deregressed genomic breeding values [(G)EBVs] were used in a GWAS to identify candidate genes and biological mechanisms affecting HR, STAY and PNM. In the functional analyses, genes within 100 kb down and upstream of each significant SNP marker, were researched. Thirty haplotypes had high expected frequency, while no homozygotes were observed. Most of the alleles present in these haplotypes had a negative mean effect for PNM, HR and STAY. The GWAS revealed significant SNP markers involved in different physiological mechanisms, leading to harmful effect on the three traits. The functional analysis revealed 26 genes enriched for 19 GO terms. Most of the GO terms found for biological processes, molecular functions and pathways were related to tissue development and the immune system. More phenotypes underlying these putative regions in this population could be the subject of future investigation. Tests to find putative lethal haplotype carriers could help breeders to eliminate them from the population or manage matings in order to avoid homozygous.

Natural levels of Rhipicephalus microplus infestation and Anaplasma marginale infection in Angus and Ultrablack calves.

Infections by Anaplasma marginale and infestations by Rhipicephalus microplus occur endemically in Brazil, representing an obstacle to expanding the use of taurine breeds, which are more susceptible. In this study, the levels of infection by A. marginale and infestation by R. microplus were monitored in 31 calves that were either purebred or had a high degree of taurine blood: 17 Angus (100% taurine) and 14 Ultrablack (ca. 82% taurine and 18% Zebu). The animals were evaluated on 13 occasions at 12-day intervals. The levels of A. marginale infection were determined by quantification of DNA copy number (CN) by qPCR, and ticks were monitored by two methods: counting adult females (≥ 4.5 mm) and scoring the level of tick infestation considering all visible instars in the animals' bodies. No significant effects were observed between the means of CN of A. marginale, tick counts and scores among Angus and Ultrablack animals. The repeatability estimates for CN of A. marginale, tick counts and tick scores were 0.53, 0.12 and 0.16, respectively. The correlations between CN and tick counts and scores were close to zero, whereas the correlations between tick assessment methods were 0.57. The absence of differences between the two genetic groups indicates, under the conditions of the present study, that the low degree of Zebu blood did not influence the levels of infection by A. marginale or infestation by R. microplus. The results also suggest that the evaluation of the levels of infestation by ticks using scores can provide information closer to the real infestation rate considering that it uses all the visible instars of the parasites.

Selection effect for growth traits on energy requirements in beef Nelore steers.

Growth data of 77,372 Nelore steers were used to estimate the selection effect on energy requirements considering two beef production systems: cow-calf and slaughter cycles. All the animals had measures from 120 days to 7 years old. The parameters necessary to evaluate the selection effect on energy requirements were obtained by random regression analysis using Legendre polynomials. The models included additive direct and maternal effects, and animal and maternal permanent environmental effects as random terms. Contemporary group and dam age at calving (linear and quadratic effect) were included as fixed effects, and orthogonal Legendre polynomials of animal age (cubic regression) were considered random covariables. The coefficients from the model M3353_5 were used to calculate the genetic gains necessary to predict the increase in phenotypes. The selection was simulated for body weight (BW) and weight gain (WG) at different ages and energy requirements were calculated using NRC equations. The cost of feed was calculated for a cow-calf and slaughter cycle of production considering a system of Brachiaria decumbens pasture without supplementation. In slaughter system, the selection for weight of 365 days of age is the best option. In cow-calf systems, the selection W120 is the best choice.

Signatures of selection in Nelore cattle revealed by whole-genome sequencing data.

Nelore cattle breed was farmed worldwide due to its economic importance in the beef market and adaptation to the tropics. In Brazil, purebred Nelore animals (PO) receive a certificate from the breeders' association based on the animal's genealogy and morphological characterization. The top 20 to 30% of the superior animals are eligible to receive the Special Certificate of Identification and Production (CEIP), meaning animals from this category were selected and evaluated in a breeding program to improve economically important traits. We used whole-genome sequencing and approaches based on haplotype differentiation and allelic differentiation to detect regions of selection signatures in Nelore cattle by comparing animals from PO and CEIP categories. From a total of 150 animals, a hierarchical clustering analysis was performed to choose the more unrelated animals from each category (16 PO and 40 CEIP). The hapFLK statistic was performed, and extensions of hapFLK values were investigated considering continuous regions with significant q-values. The Weir and Cockerham's Fst estimator (wcFst) was computed using the GPAT++ software library. The total of 82,326 SNPs with hapFLK values passed the FDR control (q-value<0.05), and 718 segments were target as signatures of selection. A total of 1713 highly differentiated genomic regions were identified based on the segmentFst approach. The signatures of selection were spread across the genome. Annotation of overlapping selection signature regions between the two methods revealed 118 genes in common. A variant located within the 3' region of the BOLA-DRB3 gene was found as a promising candidate polymorphism. Within genomic regions that deserves attention, we found genes previously associated with adaptation to tropical environments (HELB), growth and navel size (HMGA2), fat deposition and domestication (IRAK3), and feed efficiency and postmortem carcass traits (GABRG3). The genes BOLA-DQA2, BOLA-DQB, BOLA-DQA5, BOLA-DQA1, BOLA-DRB3, ENSBTAG00000038397 on chromosome 23 are part of the Bovine Major Histocompatibility Complex (MHC) Class II gene family, representing good candidates for immune response and adaptation to tropical conditions. The BoLA family genes and the interaction of ROBO1 with SLIT genes appeared in the enrichment results. Genomic regions located in intronic regions were also identified and might play a regulatory role in traits under selection in PO and CEIP subpopulations. The regions here identified contribute to our knowledge regarding genes and variants that have an important role in complex traits selected in this breed.

Genomic study of the resilience of buffalo cows to a negative energy balance.

The research article was carried out with the objective of studying the genetic variation on the resilience of buffaloes to negative energy balance-NEB (measured by changes in body weight in early lactation)-as well as investigating genomic regions of interest for this trait. A model of reaction norms was used, considering milk production as the trait to be analyzed and solutions of the contemporary groups to weight changes as environmental gradient. In this methodology, the genetic value of the slope represents the measure of resilience of the animals. After the estimation step, a genome-wide association analysis was performed for the slope of the reaction norms model, to obtain a list of windows and associated genes. The heritability estimates for milk production over the resilience gradient ranged from 0.13 to 0.28, with lower values in the intermediate environmental groups. Regarding the productive resilience of dairy buffalo cows to NEB, the genomic windows with the highest contribution to the genetic variance were detected on chromosomes BBU 1, 2, 3, 4, 9, 12, 19, and 21. A functional analysis of the genes described in the selected windows indicated association with metabolic routes related to growth and immunity of the animals, with an emphasis on the STAT6 gene. The results presented indicate that there is for this trait genetic variation to be used as selection criteria, in addition to genomic regions that can increase the precision of the selection.

Development of genomic predictions for Angus cattle in Brazil incorporating genotypes from related American sires.

Genomic prediction has become the new standard for genetic improvement programs, and currently, there is a desire to implement this technology for the evaluation of Angus cattle in Brazil. Thus, the main objective of this study was to assess the feasibility of evaluating young Brazilian Angus (BA) bulls and heifers for 12 routinely recorded traits using single-step genomic BLUP (ssGBLUP) with and without genotypes from American Angus (AA) sires. The second objective was to obtain estimates of effective population size (Ne) and linkage disequilibrium (LD) in the Brazilian Angus population. The dataset contained phenotypic information for up to 277,661 animals belonging to the Promebo breeding program, pedigree for 362,900, of which 1,386 were genotyped for 50k, 77k, and 150k single nucleotide polymorphism (SNP) panels. After imputation and quality control, 61,666 SNPs were available for the analyses. In addition, genotypes from 332 American Angus (AA) sires widely used in Brazil were retrieved from the AA Association database to be used for genomic predictions. Bivariate animal models were used to estimate variance components, traditional EBV, and genomic EBV (GEBV). Validation was carried out with the linear regression method (LR) using young-genotyped animals born between 2013 and 2015 without phenotypes in the reduced dataset and with records in the complete dataset. Validation animals were further split into progeny of BA and AA sires to evaluate if their progenies would benefit by including genotypes from AA sires. The Ne was 254 based on pedigree and 197 based on LD, and the average LD (±SD) and distance between adjacent single nucleotide polymorphisms (SNPs) across all chromosomes were 0.27 (±0.27) and 40743.68 bp, respectively. Prediction accuracies with ssGBLUP outperformed BLUP for all traits, improving accuracies by, on average, 16% for BA young bulls and heifers. The GEBV prediction accuracies ranged from 0.37 (total maternal for weaning weight and tick count) to 0.54 (yearling precocity) across all traits, and dispersion (LR coefficients) fluctuated between 0.92 and 1.06. Inclusion of genotyped sires from the AA improved GEBV accuracies by 2%, on average, compared to using only the BA reference population. Our study indicated that genomic information could help us to improve GEBV accuracies and hence genetic progress in the Brazilian Angus population. The inclusion of genotypes from American Angus sires heavily used in Brazil just marginally increased the GEBV accuracies for selection candidates.

There was a desire to implement genomic selection for Angus cattle in Brazil since the technology has been proved to increase genetic gain in animal breeding programs. Single-step genomic best linear unbiased prediction (ssGBLUP), which simultaneously combines pedigree and genomic information, was used to estimate individuals' genomic breeding values (GEBV) or genetic merit. Genomic selection can accelerate genetic progress by increasing accuracy, especially in young animals without progeny. The accuracy of GEBV can also be improved by combing data from other countries to increase the reference population (i.e., genotyped and phenotyped animals) in small, genotyped populations. Thus, the main objective of this study was to evaluate the accuracy of GEBV for young Brazilian Angus (BA) bulls and heifers with ssGBLUP, including or not the genotypes from American Angus sires. The accuracies with ssGBLUP were higher than those from traditional BLUP (EBV calculated from pedigree), improving accuracies by, on average, 16% for young bulls and heifers. Including genotypes from American Angus sires heavily used in Brazil just marginally increased the GEBV accuracies for selection candidates.

How long does the mRNA remains stable in untreated whole bovine blood?

BACKGROUND: High quality and quantity of messenger RNA (mRNA) are required for accuracy of gene expression studies and other RNA-based downstream applications. Since RNA is considered a labile macromolecular prone to degradation, which may result in falsely altered gene expression patterns, several commercial stabilizing reagents have been developed aiming to keep RNA stable for long period. However, for studies involving large number of experimental samples, the high costs related to these specific reagents may constitute a barrier. METHODS AND RESULTS: In this context the present study was designed aiming to evaluate the stability of mRNA in whole bovine blood collected in EDTA tubes during storage at common fridge (4 °C). Whole blood samples were collected from six Holstein calves and submitted to RNA extraction in each different interval: immediately after blood sampling (< 2 h), at 1-day post-sampling (dps), 2 dps, 3 dps, 7 dps and 14dps intervals. RNA integrity and purity were evaluated, and RT-qPCR assays were run using seven different genes (B2M, ACTB, PPIA, GAPDH, YWHAZ, CD4 and IFN-γ) aiming to evaluate the presence of altered gene transcription during storage. All extracted RNA samples presented high purity, while optimal integrity and unaltered gene expression were observed in whole experimental group up to 3 days of storage. CONCLUSION: Bovine blood RNA remained stable in K3EDTA tubes for 3 days stored at common fridge and can be successfully and accurately used for gene expression studies.

Semi-quantitative evaluation of Babesia bovis and B. bigemina infection levels estimated by HRM analysis.

Bovine babesiosis is economically the most important arthropod-borne disease of cattle worldwide. The most significant damage caused by bovine babesiosis is attributed to Babesia bovis due to its higher pathogenicity. This study aimed to develop a real-time PCR method followed by HRM (high-resolution melting) analysis for the simultaneous detection of B. bovis and B. bigemina, enabling a semi-quantitative analysis of Babesia levels using a single-tube reaction. The HRM was compared with real-time PCR using species-specific hydrolysis probes. The HRM analysis allowed to differentiate both Babesia species and was sensitive in the detection and differentiation of 10% for each Babesia species in the sample. Our results suggest the use of this method to estimate the prevalence of infections by B. bovis or B. bigemina as an alternative to the methods of absolute quantification by real-time PCR since it neither requires precise estimates of the number of DNA loads nor the construction of calibration curves. The simultaneous detection of the two Babesia species can be used to characterise the infection levels in cattle populations from different geographical regions, allowing a better control of these diseases.

Comparison of methods for predicting genomic breeding values for growth traits in Nellore cattle.

The objective of this study was to evaluate the accuracy of genomic predictions of growth traits in Nellore cattle. Data from 5064 animals belonging to farms that participate in the Conexão DeltaGen and PAINT breeding programs were used. Genotyping was performed with the Illumina BovineHD BeadChip (777,962 SNPs). After quality control of the genomic data, 412,993 SNPs were used. Deregressed EBVs (DEBVs) were calculated using the estimated breeding values (EBVs) and accuracies of birth weight (BW), weight gain from birth to weaning (GBW), postweaning weight gain (PWG), yearling height (YH), and cow weight (CW) provided by GenSys. Three models were used to estimate marker effects: genomic best linear unbiased prediction (GBLUP), BayesCπ, and improved Bayesian least absolute shrinkage and selection operator (IBLASSO). The prediction ability of genomic estimated breeding value (GEBVs) was estimated by the average Pearson correlation between DEBVs and GEBVs, predicted with the different methodologies in the validation populations. The regression coefficients of DEBVs on GEBVs in the validation population were calculated and used as indicators of prediction bias of GEBV. In general, the Bayesian methods provided slightly more accurate predictions of genomic breeding values than GBLUP. The BayesCπ and IBLASSO were similar for all traits (BW, GBW, PWG, and YH), except for CW. Thus, there does not seem to be a more suitable method for the estimation of SNP effects and genomic breeding values. Bayesian regression models are of interest for future applications of genomic selection in this population, but further improvements are needed to reduce deflation of their predictions.

Genomic Study of Babesia bovis Infection Level and Its Association With Tick Count in Hereford and Braford Cattle.

Bovine babesiosis is a tick-borne disease caused by intraerythrocytic protozoa and leads to substantial economic losses for the livestock industry throughout the world. Babesia bovis is considered the most pathogenic species, which causes bovine babesiosis in Brazil. Genomic data could be used to evaluate the viability of improving resistance against B. bovis infection level (IB) through genomic selection, and, for that, knowledge of genetic parameters is needed. Furthermore, genome-wide association studies (GWAS) could be conducted to provide a better understanding of the genetic basis of the host response to B. bovis infection. No previous work in quantitative genetics of B. bovis infection was found. Thus, the objective of this study was to estimate the genetic correlation between IB and tick count (TC), evaluate predictive ability and applicability of genomic selection, and perform GWAS in Hereford and Braford cattle. The single-step genomic best linear unbiased prediction method was used, which allows the estimation of both breeding values and marker effects. Standard phenotyping was conducted for both traits. IB quantifications from the blood of 1,858 animals were carried using quantitative PCR assays. For TC, one to three subsequent tick counts were performed by manually counting adult female ticks on one side of each animal's body that was naturally exposed to ticks. Animals were genotyped using the Illumina BovineSNP50 panel. The posterior mean of IB heritability, estimated by the Bayesian animal model in a bivariate analysis, was low (0.10), and the estimations of genetic correlation between IB and TC were also low (0.15). The cross-validation genomic prediction accuracy for IB ranged from 0.18 to 0.35 and from 0.29 to 0.32 using k-means and random clustering, respectively, suggesting that genomic predictions could be used as a tool to improve genetics for IB, especially if a larger training population is developed. The top 10 single nucleotide polymorphisms from the GWAS explained 5.04% of total genetic variance for IB, which were located on chromosomes 1, 2, 5, 6, 12, 17, 18, 16, 24, and 26. Some candidate genes participate in immunity system pathways indicating that those genes are involved in resistance to B. bovis in cattle. Although the genetic correlation between IB and TC was weak, some candidate genes for IB were also reported in tick infestation studies, and they were also involved in biological resistance processes. This study contributes to improving genetic knowledge regarding infection by B. bovis in cattle.

A polymorphic CD4 epitope related to increased susceptibility to Babesia bovis in Canchim calves.

Different allelic forms of bovine CD4 were previously described in cattle and were also observed in Canchim calves examined in the present experiment. However, the functional relevance of these different CD4 phenotypes has not yet been investigated. CD4 + T helper cells are known to play a central role in immune control against Babesia bovis infection. Thus, our study aimed to compare the profiles of immune cells, specific antibody titers and blood infection levels measured by qPCR (quantitative polymerase chain reaction) in calves naturally infected with B. bovis, phenotyped as CD4- (absence of anti-CD4 staining), CD4 + (intermediate staining) or CD4 ++ (high staining). The CD4 mRNA precursor was also measured in these animals. Calves with the CD4- phenotype showed higher amounts of B. bovis DNA in blood samples, compared to the other CD4 phenotypes. It was also observed that these calves with higher levels of infection had lower amounts of natural killer cells and higher expression of the CD4 gene, which can be interpreted as a compensation for the failure of the altered CD4 receptor to recognize relevant B. bovis epitopes.

Publisher Correction: Uncovering Sub-Structure and Genomic Profiles in Across-Countries Subpopulations of Angus Cattle.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Identification of genomic regions related to age at first calving and first calving interval in water buffalo using single-step GBLUP.

In Brazil, water buffaloes have been used to produce milk for mozzarella cheese production. Consequently, the main selection criterion applied for the buffalo genetic improvement is the estimated mozzarella yield as a function of milk, fat and protein production. However, given the importance of reproductive traits in production systems, this study aimed to use techniques for identifying genomic regions that affect the age at first calving (AFC) and first calving interval (FCI) in buffalo cows and to select candidate genes for the identification of QTL and gene expression studies. The single-step GBLUP method was used for the identification of genomic regions. Windows of 1 Mb containing single-nucleotide polymorphisms were constructed and the 10 windows that explained the greatest proportion of genetic variance were considered candidate regions for each trait. Genes present into the selected windows were identified using the UOA_WB_1 assembly as the reference, and their ontology was defined with the Panther tool. Candidate regions for both traits were identified on BBU 3, 12, 21 and 22; for AFC, candidates were detected on BBU 6, 7, 8, 9 and 15 and for first calving interval on BBU 4, 14 and 19. This study identified regions with great contribution to the additive genetic variance of age at first calving and first calving interval in the population of buffalo cows studied. The ROCK2, PMVK, ADCY2, MAP2K6, BMP10 and GFPT1 genes are main candidates for reproductive traits in water dairy buffaloes, and these results may have future applications in animal breeding programs or in gene expression studies of the species.

Correlations and repeatability between Babesia spp. infection levels using two dairy cattle breeding systems.

Babesia bovis and Babesia bigemina are tick-transmitted piroplasms that cause severe damage to the livestock industry in tropical regions of the world. Recent studies demonstrated differences in infection levels of these haemoparasites among bovine breeds and variation between individual cows regarding resistance to these diseases. This study aimed to estimate the repeatability and correlations between B. bovis and B. bigemina using two cattle breeding systems, an individual system (IS) and a collective paddock system (CPS). All animals were Holstein breed, and the levels of B. bovis and B. bigemina in blood samples were estimated by quantitative polymerase chain reaction (qPCR). The estimated correlations for the B. bigemina and B. bovis DNA copy number for IS and CPS were moderate and high, respectively, whereas repeatability estimates for both systems and both Babesia species were moderate. Although we cannot infer that the type of rearing system directly influenced the correlation and repeatability coefficients, it appears that the bovine parasitemia burden may be dependent on (or determine) the parasitemia burden on ticks because the bovines remained in the same place for a longer time in both systems. Thus, the babesiosis infection levels of the ticks may have been uniform, a phenomenon that also ensures greater uniformity in cattle infection. This factor may have favored the occurrence of infected ticks leading to higher repeatability estimates and correlations. Our study confirms high variability in resistance/susceptibility between breeds, and the high correlations found may be linked to this characteristic and the most intensive breeding type of dairy cattle. Besides, under the present study conditions, the estimated correlations suggest that measuring an infection level of one Babesia species can predict the level of infection of the other.

Whole-genome sequencing provides new insights into genetic mechanisms of tropical adaptation in Nellore (Bos primigenius indicus).

Most of the knowledge about genetic variants at the sequence level in cattle is for Bos primigenius taurus populations. Here, we presented a complete genomic characterization of 52 Nellore (Bos primigenius indicus) bulls, revealing specific zebu DNA variants with putative impact in tropical adaptation and productive traits. Single nucleotide polymorphisms (SNPs) and insertion/deletion (INDELs) mutations were identified using the newest bovine reference genome ARS_UCD1.2, and variant functional consequences were predicted using the Ensembl VEP software. A total of 35,753,707 SNPs and 4,492,636 INDELs were detected and annotated to their functional effects. We identified 400 genes that comprised both, a SNP and an INDEL, of high functional impact on proteins (i.e. variants that cause protein truncation, loss of function or triggering nonsense-mediated decay). Among these, we highlight the following genes: BoLA, associated with cattle immune response to infections and reproduction aspects; HSPA8, DNAJC27, and DNAJC28, involved with thermoregulatory protective mechanisms in mammals; and many olfactory signaling pathway related genes that are important genetic factors in the evolution of mammalian species. All these functional aspects are directly related to cattle adaptability to tropical environments.

Uncovering Sub-Structure and Genomic Profiles in Across-Countries Subpopulations of Angus Cattle.

Highlighting genomic profiles for geographically distinct subpopulations of the same breed may provide insights into adaptation mechanisms to different environments, reveal genomic regions divergently selected, and offer initial guidance to joint genomic analysis. Here, we characterized similarities and differences between the genomic patterns of Angus subpopulations, born and raised in Canada (N = 382) and Brazil (N = 566). Furthermore, we systematically scanned for selection signatures based on the detection of autozygosity islands common between the two subpopulations, and signals of divergent selection, via FST and varLD tests. The principal component analysis revealed a sub-structure with a close connection between the two subpopulations. The averages of genomic relationships, inbreeding coefficients, and linkage disequilibrium at varying genomic distances were rather similar across them, suggesting non-accentuated differences in overall genomic diversity. Autozygosity islands revealed selection signatures common to both subpopulations at chromosomes 13 (63.77-65.25 Mb) and 14 (22.81-23.57 Mb), which are notably known regions affecting growth traits. Nevertheless, further autozygosity islands along with FST and varLD tests unravel particular sites with accentuated population subdivision at BTAs 7 and 18 overlapping with known QTL and candidate genes of reproductive performance, thermoregulation, and resistance to infectious diseases. Our findings indicate overall genomic similarity between Angus subpopulations, with noticeable signals of divergent selection in genomic regions associated with the adaptation in different environments.

Genome-wide association studies for growth traits in buffaloes using the single step genomic BLUP.

Growth traits are important for the profitability of buffalo breeding systems, since in general, these animals are raised both for meat and milk. In this study, the single-step genomic BLUP method was employed to prospect the genomic regions' associated with weight at standard ages of 100, 210, 365, and 550 days in a buffalo population, aiming to identify genes with stronger expression for those characteristics. We found 6, 1, 2, and 5 SNPs significantly associated (p value < 10-5) with weight at 100, 210, 365, and 550 days of age, respectively, where those SNPs respectively explained 0.164, 0.040, 0.044, and 0.213% of the additive variance of each trait. SNP AX-85099682 (BBU24) was significant for weight at 100, 210, and 365 days, indicating the existence of a possible QTL affecting the initial growth rate of buffaloes. All told, eight genes (CBLB, TRNAG-UCC, GADD45B, LOC112583811, MGAT4C, KCNMA1, SLC5A2, and TGFB1I1) were identified as candidates for the growth traits of buffaloes. However, molecular and gene expression studies are necessary to validate these genes for subsequent use in programs for genetic improvement of the species.