Search

Background: Heat stress has deleterious effects on physiological and performance traits in livestock. Within this context, using tropically adapted cattle breeds in pure herds or terminal crossbreeding schemes to explore heterosis is attractive for increasing animal production in warmer climate regions. This study aimed to identify biological processes, pathways, and potential biomarkers related to thermotolerance in Caracu, a tropically adapted beef cattle breed, by proteomic analysis of blood plasma. To achieve this goal, 61 bulls had their thermotolerance evaluated through a heat tolerance index. A subset of 14 extreme animals, including the seven most thermotolerant (HIGH group) and the seven least thermotolerant (LOW group), had their blood plasma samples used for proteomic analysis by liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS). The differentially regulated proteins detected between HIGH and LOW groups were used to perform functional enrichment analysis and a protein-protein interaction network analysis., Results: A total of 217 proteins were detected only in the HIGH thermotolerant group and 51 only in the LOW thermotolerant group. In addition, 81 and 87 proteins had significantly higher and lower abundancies in the HIGH group, respectively. Regarding proteins with the highest absolute log-fold change values, we highlighted those encoded by DUSP5, IGFALS, ROCK2, RTN4, IRAG1, and NNT genes based on their functions. The functional enrichment analysis detected several biological processes, molecular functions, and pathways related to cellular responses to stress, immune system, complement system, and hemostasis in both HIGH and LOW groups, in addition to terms and pathways related to lipids and calcium only in the HIGH group. Protein-protein interaction (PPI) network revealed as important nodes many proteins with roles in response to stress, hemostasis, immune system, inflammation, and homeostasis. Additionally, proteins with high absolute log-fold change values and proteins detected as essential nodes by PPI analysis highlighted herein are potential biomarkers for thermotolerance, such as ADRA1A, APOA1, APOB, APOC3, C4BPA, CAT, CFB, CFH, CLU, CXADR, DNAJB1, DNAJC13, DUSP5, FGA, FGB, FGG, HBA, HBB, HP, HSPD1, IGFALS, IRAG1, KNG1, NNT, OSGIN1, PROC, PROS1, ROCK2, RTN4, RYR1, TGFB2, VLDLR, VTN, and VWF., Conclusions: Identifying potential biomarkers, molecular mechanisms and pathways that act in response to heat stress in tropically adapted beef cattle contributes to developing strategies to improve performance and welfare traits in livestock under tropical climates., Competing Interests: Declarations Ethics approval All the experimental procedures were approved by the Institutional Animal Care and Use Committee at the Animal Science Institute, Nova Odessa, São Paulo, Brazil (CEUA n°292 − 19, October 7th, 2019), following guidelines for animal welfare according to State Law No. 11.977 of the São Paulo state, Brazil. Consent for publication Not applicable. Competing interests The authors declare no competing interests., (© 2024. The Author(s).)

Background: Nelore cattle play a key role in tropical production systems due to their resilience to harsh conditions, such as heat stress and seasonally poor nutrition. Monitoring their genetic diversity is essential to manage the negative impacts of inbreeding. Traditionally, inbreeding and inbreeding depression are assessed by pedigree-based coefficients (F), but recently, genetic markers have been preferred for their precision in capturing the inbreeding level and identifying animals at risk of reduced productive and reproductive performance. Hence, we compared the inbreeding and inbreeding depression for productive and reproductive performance traits in Nelore cattle using different inbreeding coefficient estimation methods from pedigree information (F Ped ), the genomic relationship matrix (F GRM ), runs of homozygosity (F ROH ) of different lengths (> 1 Mb (genome), between 1 and 2 Mb - F ROH 1-2 ; 2-4 Mb F ROH 2-4 or > 8 Mb F ROH >8 ) and excess homozygosity (F SNP )., Results: The correlation between F Ped and F ROH was lower when the latter was based on shorter segments (r = 0.15 with F ROH 1-2 , r = 0.20 with F ROH 2-4 and r = 0.28 with F ROH 4-8 ). Meanwhile, the F Ped had a moderate correlation with F SNP (r = 0.47) and high correlation with F ROH >8 (r = 0.58) and F ROH-genome (r = 0.60). The F ROH-genome was highly correlated with inbreeding based on F ROH>8 (r = 0.93) and F SNP (r = 0.88). The F GRM exhibited a high correlation with F ROH-genome (r = 0.55) and F ROH >8 (r = 0.51) and a lower correlation with other inbreeding estimators varying from 0.30 for F ROH 2-4 to 0.37 for F ROH 1-2 . Increased levels of inbreeding had a negative impact on the productive and reproductive performance of Nelore cattle. The unfavorable inbreeding effect on productive and reproductive traits ranged from 0.12 to 0.51 for F Ped , 0.19-0.59 for F GRM , 0.21-0.58 for F ROH-genome , and 0.19-0.54 for F SNP per 1% of inbreeding scaled on the percentage of the mean. When scaling the linear regression coefficients on the standard deviation, the unfavorable inbreeding effect varied from 0.43 to 1.56% for F Ped , 0.49-1.97% for F GRM , 0.34-2.2% for F ROH-genome , and 0.50-1.62% for F SNP per 1% of inbreeding. The impact of the homozygous segments on reproductive and performance traits varied based on the chromosomes. This shows that specific homozygous chromosome segments can be signs of positive selection due to their beneficial effects on the traits., Conclusions: The low correlation observed between F Ped and genomic-based inbreeding estimates suggests that the presence of animals with one unknown parent (sire or dam) in the pedigree does not account for ancient inbreeding. The ROH hotspots surround genes related to reproduction, growth, meat quality, and adaptation to environmental stress. Inbreeding depression has adverse effects on productive and reproductive traits in Nelore cattle, particularly on age at puberty in young bulls and heifer calving at 30 months, as well as on scrotal circumference and body weight when scaled on the standard deviation of the trait., (© 2024. The Author(s).)

Background: Potential functional variants (PFVs) can be defined as genetic variants responsible for a given phenotype. Ultimately, these are the best DNA markers for animal breeding and selection, especially for polygenic and complex phenotypes. Herein, we described the identification of PFVs for complex phenotypes (in this case, Feed Efficiency in beef cattle) using a systems-biology driven approach based on RNA-seq data from physiologically relevant organs., Results: The systems-biology coupled with deep molecular phenotyping by RNA-seq of liver, muscle, hypothalamus, pituitary, and adrenal glands of animals with high and low feed efficiency (FE) measured by residual feed intake (RFI) identified 2,000,936 uniquely variants. Among them, 9986 variants were significantly associated with FE and only 78 had a high impact on protein expression and were considered as PFVs. A set of 169 significant uniquely variants were expressed in all five organs, however, only 27 variants had a moderate impact and none of them a had high impact on protein expression. These results provide evidence of tissue-specific effects of high-impact PFVs. The PFVs were enriched (FDR < 0.05) for processing and presentation of MHC Class I and II mediated antigens, which are an important part of the adaptive immune response. The experimental validation of these PFVs was demonstrated by the increased prediction accuracy for RFI using the weighted G matrix (ssGBLUP+wG; Acc = 0.10 and b = 0.48) obtained in the ssGWAS in comparison to the unweighted G matrix (ssGBLUP; Acc = 0.29 and b = 1.10)., Conclusion: Here we identified PFVs for FE in beef cattle using a strategy based on systems-biology and deep molecular phenotyping. This approach has great potential to be used in genetic prediction programs, especially for polygenic phenotypes., (© 2022. The Author(s).)

Background: The cattle introduced by European conquerors during the Brazilian colonization period were exposed to a process of natural selection in different types of biomes throughout the country, leading to the development of locally adapted cattle breeds. In this study, whole-genome re-sequencing data from indicine and Brazilian locally adapted taurine cattle breeds were used to detect genomic regions under selective pressure. Within-population and cross-population statistics were combined separately in a single score using the de-correlated composite of multiple signals (DCMS) method. Putative sweep regions were revealed by assessing the top 1% of the empirical distribution generated by the DCMS statistics., Results: A total of 33,328,447 biallelic SNPs with an average read depth of 12.4X passed the hard filtering process and were used to access putative sweep regions. Admixture has occurred in some locally adapted taurine populations due to the introgression of exotic breeds. The genomic inbreeding coefficient based on runs of homozygosity (ROH) concurred with the populations' historical background. Signatures of selection retrieved from the DCMS statistics provided a comprehensive set of putative candidate genes and revealed QTLs disclosing cattle production traits and adaptation to the challenging environments. Additionally, several candidate regions overlapped with previous regions under selection described in the literature for other cattle breeds., Conclusion: The current study reported putative sweep regions that can provide important insights to better understand the selective forces shaping the genome of the indicine and Brazilian locally adapted taurine cattle breeds. Such regions likely harbor traces of natural selection pressures by which these populations have been exposed and may elucidate footprints for adaptation to the challenging climatic conditions.

Background: The aim of this study was to use transcriptome RNA-Seq data from longissimus thoracis muscle of uncastrated Nelore males to identify hub genes based on co-expression network obtained from differentially expressed genes (DEGs) associated with intramuscular fat content., Results: A total of 30 transcriptomics datasets (RNA-Seq) obtained from longissimus thoracis muscle were selected based on the phenotypic value of divergent intramuscular fat content: 15 with the highest intramuscular fat content (HIF) and 15 with the lowest intramuscular fat content (LIF). The transcriptomics datasets were aligned with a reference genome and 65 differentially expressed genes (DEGs) were identified, including 21 upregulated and 44 downregulated genes in HIF animals. The normalized count data from DEGs was then used for co-expression network construction. From the co-expression network, four modules were identified. The topological properties of the network were analyzed; those genes engaging in the most interactions (maximal clique centrality method) with other DEGs were predicted to be hub genes (PDE4D, KLHL30 and IL1RAP), which consequently may play a role in cellular and/or systemic lipid biology in Nelore cattle. Top modules screened from the gene co-expression network were identify. The two candidate modules had clear associated biological pathways related to fat development, cell adhesion, and muscle differentiation, immune system, among others. The hub genes belonged in top modules and were downregulated in HIF animals. PDE4D and IL1RAP have known effects on lipid metabolism and the immune system through the regulation of cAMP signaling. Given that cAMP is known to play a role in lipid systems, PDE4D and IL1RAP downregulation may contribute to increased levels of intracellular cAMP and thus may have effects on IF content differences in Nelore cattle. KLHL30 may have effects on muscle metabolism. Klhl protein families play a role in protein degradation. However, the downregulation of this gene and its role in lipid metabolism has not yet been clarified., Conclusions: The results reported in this study indicate candidate genes and molecular mechanisms involved in IF content difference in Nelore cattle.

Background: Genome-wide association studies (GWAS) are utilized in cattle to identify regions or genetic variants associated with phenotypes of interest, and thus, to identify design strategies that allow for the increase of the frequency of favorable alleles. Visual scores are important traits of cattle production in Brazil because they are utilized as selection criteria, helping to choose more harmonious animals. Despite its importance, there are still no studies on the genome association for these traits. This study aimed to identify genome regions associated with the traits of conformation, precocity and muscling, based on a visual score measured at weaning., Results: Bayesian approaches with BayesC and Bayesian LASSO were utilized with 2873 phenotypes of Nellore cattle for a GWAS. The animals were genotyped with Illumina BovineHD BeadChip, and a total of 309,865 SNPs were utilized after quality control. In the analyses, phenotype and deregressed breeding values were utilized as dependent variables; a threshold model was utilized for the former and a linear model for the latter. The association criterion was the percentage of genetic variance explained by SNPs found in 1 Mb-long windows. The Bayesian approach BayesC was better adjusted to the data because it could explain a larger phenotypic variance for both dependent variables., Conclusions: There were no large effects for the visual scores, indicating that they have a polygenic nature; however, regions in chromosomes 1, 3, 5, 7, 14, 15, 16, 19, 20 and 23 were identified and explained a large part of the genetic variance.

Background: The aim of this study was to assess genome-wide autozygosity in a Nellore cattle population and to characterize ROH patterns and autozygosity islands that may have occurred due to selection within its lineages. It attempts also to compare estimates of inbreeding calculated from ROH (F ROH ), genomic relationship matrix (F GRM ), and pedigree-based coefficient (F PED )., Results: The average number of ROH per animal was 55.15 ± 13.01 with an average size of 3.24 Mb. The Nellore genome is composed mostly by a high number of shorter segments accounting for 78% of all ROH, although the proportion of the genome covered by them was relatively small. The genome autozygosity proportion indicates moderate to high inbreeding levels for classical standards, with an average value of 7.15% (178.70 Mb). The average of F PED and F ROH , and their correlations (- 0.05 to 0.26) were low. Estimates of correlation between F GRM -F PED was zero, while the correlation (- 0.01 to - 0.07) between F GRM -F ROH decreased as a function of ROH length, except for F ROH > 8Mb (- 0.03). Overall, inbreeding coefficients were not high for the genotyped animals. Autozygosity islands were evident across the genome (n = 62) and their genomic location did not largely differ within lineages. Enriched terms (p < 0.01) associated with defense response to bacteria (GO:0042742), immune complex reaction (GO:0045647), pregnancy-associated glycoproteins genes (GO:0030163), and organism growth (GO:0040014) were described within the autozygotic islands., Conclusions: Low F PED- F ROH correlation estimates indicate that F PED is not the most suitable method for capturing ancient inbreeding when the pedigree does not extend back many generations and F ROH should be used instead. Enriched terms (p < 0.01) suggest a strong selection for immune response. Non-overlapping islands within the lineages greatly explain the mechanism underlying selection for functionally important traits in Nellore cattle.

Background: Runs of homozygosity (ROH) are continuous homozygous segments of the DNA sequence. They have been applied to quantify individual autozygosity and used as a potential inbreeding measure in livestock species. The aim of the present study was (i) to investigate genome-wide autozygosity to identify and characterize ROH patterns in Gyr dairy cattle genome; (ii) identify ROH islands for gene content and enrichment in segments shared by more than 50% of the samples, and (iii) compare estimates of molecular inbreeding calculated from ROH (F ROH ), genomic relationship matrix approach (F GRM ) and based on the observed versus expected number of homozygous genotypes (F HOM ), and from pedigree-based coefficient (F PED )., Results: ROH were identified in all animals, with an average number of 55.12 ± 10.37 segments and a mean length of 3.17 Mb. Short segments (ROH 1-2 Mb ) were abundant through the genomes, which accounted for 60% of all segments identified, even though the proportion of the genome covered by them was relatively small. The findings obtained in this study suggest that on average 7.01% (175.28 Mb) of the genome of this population is autozygous. Overlapping ROH were evident across the genomes and 14 regions were identified with ROH frequencies exceeding 50% of the whole population. Genes associated with lactation (TRAPPC9), milk yield and composition (IRS2 and ANG), and heat adaptation (HSF1, HSPB1, and HSPE1), were identified. Inbreeding coefficients were estimated through the application of F ROH , F GRM , F HOM , and F PED approaches. F PED estimates ranged from 0.00 to 0.327 and F ROH from 0.001 to 0.201. Low to moderate correlations were observed between F PED -F ROH and F GRM -F ROH , with values ranging from -0.11 to 0.51. Low to high correlations were observed between F ROH -F HOM and moderate between F PED -F HOM and F GRM -F HOM . Correlations between F ROH from different lengths and F PED gradually increased with ROH length., Conclusions: Genes inside ROH islands suggest a strong selection for dairy traits and enrichment for Gyr cattle environmental adaptation. Furthermore, low F PED- F ROH correlations for small segments indicate that F PED estimates are not the most suitable method to capture ancient inbreeding. The existence of a moderate correlation between larger ROH indicates that F ROH can be used as an alternative to inbreeding estimates in the absence of pedigree records.

Background: Meat tenderness is the consumer's most preferred sensory attribute. This trait is affected by a number of factors, including genotype, age, animal sex, and pre- and post-slaughter management. In view of the high percentage of Zebu genes in the Brazilian cattle population, mainly Nellore cattle, the improvement of meat tenderness is important since the increasing proportion of Zebu genes in the population reduces meat tenderness. However, the measurement of this trait is difficult once it can only be made after animal slaughtering. New technologies such as RNA-Seq have been used to increase our understanding of the genetic processes regulating quantitative traits phenotypes. The objective of this study was to identify differentially expressed genes related to meat tenderness, in Nellore cattle in order to elucidate the genetic factors associated with meat quality. Samples were collected 24 h postmortem and the meat was not aged., Results: We found 40 differentially expressed genes related to meat tenderness, 17 with known functions. Fourteen genes were up-regulated and 3 were down-regulated in the tender meat group. Genes related to ubiquitin metabolism, transport of molecules such as calcium and oxygen, acid-base balance, collagen production, actin, myosin, and fat were identified. The PCP4L1 (Purkinje cell protein 4 like 1) and BoLA-DQB (major histocompatibility complex, class II, DQ beta) genes were validated by qRT-PCR. The results showed relative expression values similar to those obtained by RNA-Seq, with the same direction of expression (i.e., the two techniques revealed higher expression of PCP4L1 in tender meat samples and of BoLA-DQB in tough meat samples)., Conclusions: This study revealed the differential expression of genes and functions in Nellore cattle muscle tissue, which may contain potential biomarkers involved in meat tenderness.

Background: Fatty acid type in beef can be detrimental to human health and has received considerable attention in recent years. The aim of this study was to identify differentially expressed genes in longissimus thoracis muscle of 48 Nellore young bulls with extreme phenotypes for fatty acid composition of intramuscular fat by RNA-seq technique., Results: Differential expression analyses between animals with extreme phenotype for fatty acid composition showed a total of 13 differentially expressed genes for myristic (C14:0), 35 for palmitic (C16:0), 187 for stearic (C18:0), 371 for oleic (C18:1, cis-9), 24 for conjugated linoleic (C18:2 cis-9, trans11, CLA), 89 for linoleic (C18:2 cis-9,12 n6), and 110 genes for α-linolenic (C18:3 n3) fatty acids. For the respective sums of the individual fatty acids, 51 differentially expressed genes for saturated fatty acids (SFA), 336 for monounsaturated (MUFA), 131 for polyunsaturated (PUFA), 92 for PUFA/SFA ratio, 55 for ω3, 627 for ω6, and 22 for ω6/ω3 ratio were identified. Functional annotation analyses identified several genes associated with fatty acid metabolism, such as those involved in intra and extra-cellular transport of fatty acid synthesis precursors in intramuscular fat of longissimus thoracis muscle. Some of them must be highlighted, such as: ACSM3 and ACSS1 genes, which work as a precursor in fatty acid synthesis; DGAT2 gene that acts in the deposition of saturated fat in the adipose tissue; GPP and LPL genes that support the synthesis of insulin, stimulating both the glucose synthesis and the amino acids entry into the cells; and the BDH1 gene, which is responsible for the synthesis and degradation of ketone bodies used in the synthesis of ATP., Conclusion: Several genes related to lipid metabolism and fatty acid composition were identified. These findings must contribute to the elucidation of the genetic basis to improve Nellore meat quality traits, with emphasis on human health. Additionally, it can also contribute to improve the knowledge of fatty acid biosynthesis and the selection of animals with better nutritional quality.

Background: Saturated fatty acids can be detrimental to human health and have received considerable attention in recent years. Several studies using taurine breeds showed the existence of genetic variability and thus the possibility of genetic improvement of the fatty acid profile in beef. This study identified the regions of the genome associated with saturated, mono- and polyunsaturated fatty acids, and n-6 to n-3 ratios in the Longissimus thoracis of Nellore finished in feedlot, using the single-step method., Results: The results showed that 115 windows explain more than 1 % of the additive genetic variance for the 22 studied fatty acids. Thirty-one genomic regions that explain more than 1 % of the additive genetic variance were observed for total saturated fatty acids, C12:0, C14:0, C16:0 and C18:0. Nineteen genomic regions, distributed in sixteen different chromosomes accounted for more than 1 % of the additive genetic variance for the monounsaturated fatty acids, such as the sum of monounsaturated fatty acids, C14:1 cis-9, C18:1 trans-11, C18:1 cis-9, and C18:1 trans-9. Forty genomic regions explained more than 1 % of the additive variance for the polyunsaturated fatty acids group, which are related to the total polyunsaturated fatty acids, C20:4 n-6, C18:2 cis-9 cis12 n-6, C18:3 n-3, C18:3 n-6, C22:6 n-3 and C20:3 n-6 cis-8 cis-11 cis-14. Twenty-one genomic regions accounted for more than 1 % of the genetic variance for the group of omega-3, omega-6 and the n-6:n-3 ratio., Conclusions: The identification of such regions and the respective candidate genes, such as ELOVL5, ESSRG, PCYT1A and genes of the ABC group (ABC5, ABC6 and ABC10), should contribute to form a genetic basis of the fatty acid profile of Nellore (Bos indicus) beef, contributing to better selection of the traits associated with improving human health.

Background: Knowledge of the linkage disequilibrium (LD) between markers is important to establish the number of markers necessary for association studies and genomic selection. The objective of this study was to evaluate the extent of LD in Nellore cattle using a high density SNP panel and 795 genotyped steers., Results: After data editing, 446,986 SNPs were used for the estimation of LD, comprising 2508.4 Mb of the genome. The mean distance between adjacent markers was 4.90 ± 2.89 kb. The minor allele frequency (MAF) was less than 0.20 in a considerable proportion of SNPs. The overall mean LD between marker pairs measured by r(2) and |D'| was 0.17 and 0.52, respectively. The LD (r(2)) decreased with increasing physical distance between markers from 0.34 (1 kb) to 0.11 (100 kb). In contrast to this clear decrease of LD measured by r(2), the changes in |D'| indicated a less pronounced decline of LD. Chromosomes BTA1, BTA27, BTA28 and BTA29 showed lower levels of LD at any distance between markers. Except for these four chromosomes, the level of LD (r(2)) was higher than 0.20 for markers separated by less than 20 kb. At distances < 3 kb, the level of LD was higher than 0.30. The LD (r(2)) between markers was higher when the MAF threshold was high (0.15), especially when the distance between markers was short., Conclusions: The level of LD estimated for markers separated by less than 30 kb indicates that the High Density Bovine SNP BeadChip will likely be a suitable tool for prediction of genomic breeding values in Nellore cattle.