Genetic Architecture of Carcass and Meat Quality Traits in Montana Tropical® Composite Beef Cattle

Genomic Prediction Using Imputed Whole-Genome Sequence Data in Australian Angus Cattle

Whole-genome sequence (WGS) data was used to estimate genomic breeding values for growth and carcass traits in Australian Angus cattle. The study aimed to compare the accuracy and bias of genomic predictions with three marker densities, including 50K, high-density (HD) and WGS. The dataset used in this study consisted of animals born between 2013 and 2022. Body weight traits included birthweight, weight at 400 days and weight at 600 days of age. The carcass traits were carcass weight, carcass intramuscular fat and carcass marbling score. The accuracy and bias of prediction were assessed using the cross-validation. Further, for the growth traits, animals in the validation group were subdivided into two subgroups, which were moderately or highly related to the reference. Genomic best linear unbiased prediction (GBLUP) was used to compare genomic predictions with the three marker densities. The prediction accuracies were generally similar across the marker densities, ranging between 0.61 and 0.68 for the body weight traits and between 0.40 and 0.52 for the carcass traits. However, the accuracies marginally decreased as the marker density increased for all the traits studied. A similar lack of difference was found when considering the accuracy by the relatedness subgroups. The results indicated that no meaningful difference in prediction accuracy was estimated when comparing the three marker densities due to the population structure. In conclusion, there was no substantial improvement in genomic prediction when using the WGS in this study.

Genetic spatialization in Montana cattle breed

Brazil is among the world's top beef producers and exporters, and the Montana breed of livestock has gained popularity due to its exceptional meat quality and ability to adapt to the tropical climate. This differentiates it from other breeds, such as the Zebu and Taurine. A methodology that mapped sixteen farms across five Brazilian states was adopted to investigate the link between SNP genetic data and animal collection location. The data was georeferenced using the Geographic Information System and the ArcGIS 10.3 software, considering a geographic coordinate system and the WGS 84 datum. The Alleles in Space (AIS) software was used to study the relationship between the SNP genetic data and the animal collection location. Various analyses were conducted, including the Mantel correlation, genetic barrier identification, and spatial autocorrelation. The results suggest that there is a correlation between genetic and geographical distances, supporting the idea that geographical distance can affect genetic diversity. This information is vital to understand the genetic structure of populations concerning geographical space and is crucial for optimizing Brazilian agricultural production.

Identification of Heilongjiang crossbred beef cattle pedigrees and reveals functional genes related to economic traits based on whole-genome SNP data

Introduction: To enhance the beef cattle industry, Heilongjiang Province has developed a new Crossbred beef cattle variety through crossbreeding with exotic commercial breeds. This new variety exhibits relatively excellent meat quality, and efficient reproductive performance, catering to market demands. Method: This study employed whole genome resequencing technology to analyze the genetic pedigree and diversity of 19 Heilongjiang Crossbred beef cattle, alongside 59 published genomes from East Asian, Eurasian, and European taurine cattle as controls. In addition, genes related to production traits were also searched by identifying Runs of Homozygosity (ROH) islands and important fragments from ancestors. Results: A total of 14,427,729 biallelic SNPs were discovered, with the majority located in intergenic and intron regions and a small percentage in exon regions, impacting protein function. Population genetic analyses including Principal Component Analysis (PCA), Neighbor-Joining (NJ) tree, and ADMIXTURE identified Angus, Holstein, and Mishima as the main ancestors of Crossbred beef cattle. In genetic diversity analysis, nucleotide diversity, linkage disequilibrium, and inbreeding coefficient analysis reveal that the genetic diversity of Crossbred beef cattle is at a moderate level, and a higher inbreeding coefficient indicates the need for careful breeding management. In addition, some genes related to economic traits are identified through the identification of Runs of Homozygosity (ROH) islands and important fragments from ancestors. Conclusion: This comprehensive genomic characterization supports the targeted improvement of economically important traits in Crossbred beef cattle, facilitating advanced breeding strategies.

Genetic background of hematological parameters in Holstein cattle based on genome-wide association and RNA sequencing analyses

Hematological parameters refer to the assessment of changes in the number and distribution of blood cells, including leukocytes (LES), erythrocytes (ERS), and platelets (PLS), which are essential for the early diagnosis of hematological system disorders and other systemic diseases in livestock. In this context, the primary objectives of this study were to investigate the genomic background of 19 hematological parameters in Holstein cattle, focusing on LES, ERS, and PLS blood components. Genetic and phenotypic (co)variances of hematological parameters were calculated based on the average information restricted maximum likelihood method and 1,610 genotyped individuals and 5,499 hematological parameter records from 4,543 cows. Furthermore, we assessed the genetic relationship between these hematological parameters and other economically important traits in dairy cattle breeding programs. We also carried out genome-wide association studies and candidate gene analyses. Blood samples from 21 primiparous cows were used to identify candidate genes further through RNA sequencing (RNA-seq) analyses. Hematological parameters generally exhibited low-to-moderate heritabilities ranging from 0.01 to 0.29, with genetic correlations between them ranging from -0.88 ± 0.09 (between mononuclear cell ratio and lymphocyte cell ratio) to 0.99 ± 0.01 (between white blood cell count and granulocyte cell count). Furthermore, low-to-moderate approximate genetic correlations between hematological parameters with one longevity, 4 fertility, and 5 health traits were observed. One hundred ninety-nine significant SNP located primarily on the Bos taurus autosomes (BTA) BTA4, BTA6, and BTA8 were associated with 16 hematological parameters. Based on the RNA-seq analyses, 6,687 genes were significantly downregulated and 4,119 genes were upregulated when comparing 2 groups of cows with high and low phenotypic values. By integrating genome-wide association studies (GWAS), RNA-seq, and previously published results, the main candidate genes associated with hematological parameters in Holstein cattle were ACRBP, ADAMTS3, CANT1, CCM2L, CNN3, CPLANE1, GPAT3, GRIP2, PLAGL2, RTL6, SOX4, WDFY3, and ZNF614. Hematological parameters are heritable and moderately to highly genetically correlated among themselves. The large number of candidate genes identified based on GWAS and RNA-seq indicate the polygenic nature and complex genetic determinism of hematological parameters in Holstein cattle.

Selective signatures in composite MONTANA TROPICAL beef cattle reveal potential genomic regions for tropical adaptation

Genomic regions related to tropical adaptability are of paramount importance for animal breeding nowadays, especially in the context of global climate change. Moreover, understanding the genomic architecture of these regions may be very relevant for aiding breeding programs in choosing the best selection scheme for tropical adaptation and/or implementing a crossbreeding scheme. The composite MONTANA TROPICAL® population was developed by crossing cattle of four different biological types to improve production in harsh environments. Pedigree and genotype data (51962 SNPs) from 3215 MONTANA TROPICAL® cattle were used to i) characterize the population structure; ii) identify signatures of selection with complementary approaches, i.e. Integrated Haplotype Score (iHS) and Runs of Homozygosity (ROH); and iii) understand genes and traits related to each selected region. The population structure based on principal components had a weak relationship with the genetic contribution of the different biological types. Clustering analyses (ADMIXTURE) showed different clusters according to the number of generations within the composite population. Considering results of both selection signatures approaches, we identified only one consensus region on chromosome 20 (35399405-40329703 bp). Genes in this region are related to immune function, regulation of epithelial cell differentiation, and cell response to ionizing radiation. This region harbors the slick locus which is related to slick hair and epidermis anatomy, both of which are related to heat stress adaptation. Also, QTLs in this region were related to feed intake, milk yield, mastitis, reproduction, and slick hair coat. The signatures of selection detected here arose in a few generations after crossbreeding between contrasting breeds. Therefore, it shows how important this genomic region may be for these animals to thrive in tropical conditions. Further investigations on sequencing this region can identify candidate genes for animal breeding and/or gene editing to tackle the challenges of climate change.

Genome-wide scans identify biological and metabolic pathways regulating carcass and meat quality traits in beef cattle

Genome association studies (GWAS) provides knowledge about the genetic architecture of beef-related traits that allow linking the target phenotype to genomic information aiding breeding decision. Thus, the present study aims to uncover the genetic mechanism involved in carcass (REA: rib eye area, BF: backfat thickness, and HCW: hot carcass weight) and meat quality traits (SF: shear-force, MARB: marbling score, and IMF: intramuscular fat content) in Nellore cattle. For this, 6910 young bulls with phenotypic information and 23,859 animals genotyped with 435 k markers were used to perform the weighted single-step GBLUP (WssGBLUP) approach, considering two iterations. The top 10 genomic regions explained 8.13, 11.81, and 9.58% of the additive genetic variance, harboring a total of 119, 143, and 95 positional candidate genes for REA, BF, and HCW, respectively. For meat quality traits, the top 10 windows explained a large proportion of the total genetic variance for SF (14.95%), MARB (17.56%), and IMF (21.41%) surrounding 92, 155, and 111 candidate genes, respectively. Relevant candidate genes (CAST, PLAG1, XKR4, PLAGL2, AQP3/AQP7, MYLK2, WWOX, CARTPT, and PLA2G16) are related to physiological aspects affecting growth, carcass, meat quality, feed intake, and reproductive traits by signaling pathways controlling muscle control, key signal metabolic molecules INS / IGF-1 pathway, lipid metabolism, and adipose tissue development. The GWAS results provided insights into the genetic control of the traits studied and the genes found are potential candidates to be used in the improvement of carcass and meat quality traits.

Genotype-by-gestational thermal environment interaction and its impact on the future performance of tropical composite beef cattle offspring

With global warming, there are growing challenges for raising taurine and composite beef cattle populations in tropical regions, including elevated temperatures, limited forage availability, parasite infestation, and infectious diseases. These environmental factors can trigger specific physiological responses in the developing fetus, which may have long-term implications on its performance. Therefore, the main objective of this study was to assess the influence of naturally induced thermal stress during the gestation period on the subsequent performance of tropical composite beef cattle progeny. Furthermore, we aimed to investigate the impact of genotype-by-gestational thermal environment interaction (G×Eg) on traits under selection pressure in the breeding population. A total of 157,414 animals from 58 farms located in various Brazilian states were recorded for birth weight (BW), preweaning weight gain (PWG), yearling weight (YW), hip height (HH), scrotal circumference (SC), and days to first calving (DFC). We first applied a linear regression model to the BW data, which revealed that the last 40 d of gestation were suitable for calculating the mean temperature humidity index (THIg). Subsequent regression analyses revealed that for every 10-unit increase in THIg, detrimental effects of approximately 1.13% to 16.34% are expected for all traits evaluated. Genetic parameters were estimated through a reaction norm model using THIg as the environmental descriptor. The posterior means of heritability estimates (SD) were 0.35 (0.07), 0.25 (0.03), 0.31 (0.03), 0.37 (0.01), 0.29 (0.07), and 0.20 (0.09) for the direct effect of BW, PWG, YW, HH, SC, and DFC, respectively. These estimates varied along the range of THIg values, suggesting a variable response to selection depending on the thermal environment during gestation. Genetic correlation estimates between more divergent THIg values were low or negative for YW, PWG, and DFC, indicating that the best-performing individuals at low THIg values may not perform as well at high THIg values and vice versa. Overall, thermal stress during gestation impacts the future performance of beef cattle offspring. Our results indicate the need for developing effective breeding strategies that take into account G×Eg effects and the re-ranking of breeding animals along the THIg scale, particularly for traits such as DFC that are highly sensitive to thermal stress.

Genetic Association Analysis of Copy Number Variations for Meat Quality in Beef Cattle

Meat quality is an economically important trait for global food production. Copy number variations (CNVs) have been previously implicated in elucidating the genetic basis of complex traits. In this article, we detected a total of 112,198 CNVs and 10,102 CNV regions (CNVRs) based on the Bovine HD SNP array. Next, we performed a CNV-based genome-wide association analysis (GWAS) of six meat quality traits and identified 12 significant CNV segments corresponding to eight candidate genes, including PCDH15, CSMD3, etc. Using region-based association analysis, we further identified six CNV segments relevant to meat quality in beef cattle. Among these, TRIM77 and TRIM64 within CNVR4 on BTA29 were detected as candidate genes for backfat thickness (BFT). Notably, we identified a 34 kb duplication for meat color (MC) which was supported by read-depth signals, and this duplication was embedded within the keratin gene family including KRT4, KRT78, and KRT79. Our findings will help to dissect the genetic architecture of meat quality traits from the aspects of CNVs, and subsequently improve the selection process in breeding programs.

Genome-wide association study for primal cut lean traits in Canadian beef cattle

This study identified genomic variants and underlying candidate genes related to the whole carcass and individual primal cut lean content in Canadian commercial crossbred beef cattle. Genotyping information of 1035 crossbred beef cattle were available alongside estimated and actual carcass lean meat yield and individual primal cut lean content in all carcasses. Significant fixed effects and covariates were identified and included in the animal model. Genome-wide association analysis were implemented using the weighted single-step genomic best linear unbiased prediction (WssGBLUP). A number of candidate genes identified linked to lean tissue production were unrelated to estimated lean meat yield and were specific to the actual lean traits. Among these, 41 genes were common for actual lean traits, on specific regions of BTA4, BTA13 and BTA25 indicating potential involvement in lean mass synthesis. Therefore, the results suggested the inclusion of primal cut lean traits as a selection objective in breeding programs with consideration of further functional studies of the identified genes could help in optimizing lean yield for maximal carcass value.

Estimation of genetic parameters for primal tissue component traits in commercial crossbred beef cattle

Knowledge of genetic parameters is required to select for optimal yield of primal cuts that may be used as the selection criteria for designing future breeding programs. This study aimed to estimate the heritability, as well as genetic and phenotypic correlations of primal cut lean and fat tissue components, and carcass traits in Canadian crossbred beef cattle. All tissue component traits presented a medium to high heritability (lean 0.41 to 0.61; fat 0.46 to 0.62; bone 0.22 to 0.48), which indicates a probable increase in their response to genetic selection. In addition, high genetic correlations were found among the primal cut lean trait group (0.63 to 0.94) and fat trait group (0.63 to 0.94), as well as strong negative correlations between lean and fat component traits (-0.63 to -1). Therefore, results suggested inclusion of primal cut tissue composition traits as a selection objective in breeding programs with consideration of correlations among the traits could help in optimizing lean yield for the highest carcass value.

Understanding the phenotypic and genetic background of the lactose content in Sarda dairy sheep

Lactose, the principal carbohydrate found in milk, plays an important role in the physiological processes of milk production because it is related to milk volume, and it is responsible for the osmotic equilibrium between blood and milk in the mammary gland. In this study, factors affecting lactose content (LC) in sheep milk are investigated. For this purpose, 2,358 test-day records were sampled from 509 ewes (3-7 records per animal). The LC and other main milk traits were analyzed using a mixed linear model that included days in milk (DIM) class, parity, lambing month, and type of lambing as fixed effects and animal, permanent environment, and flock test day as random effects. The pedigree-based approach was used to estimate the heritability and repeatability of LC. Moreover, the genomic background of LC was investigated through a GWAS. The LC was affected by all tested factors (i.e., DIM class, parity, lambing month, and type of lambing). Low heritability (0.10 ± 0.05) and moderate repeatability (0.42 ± 0.02) were estimated for LC. High negative genetic correlations were estimated between LC and NaCl (-0.99 ± 0.01) and between LC and somatic cell count (-0.94 ± 0.05). Only 2 markers passed the chromosome-wide Bonferroni threshold. Results of the present study, although obtained on a relatively small sample, suggest the possibility to include LC in the breeding programs, particularly because of its strong relationship with NaCl and somatic cell count.

A look under the hood of genomic-estimated breed compositions for brangus cattle: What have we learned?

The Brangus cattle were developed to utilize the superior traits of Angus and Brahman cattle. Their genetic compositions are expected to be stabilized at 3/8 Brahman and 5/8 Angus. Previous studies have shown more than expected Angus lineage with Brangus cattle, and the reasons are yet to be investigated. In this study, we revisited the breed compositions for 3,605 Brangus cattle from three perspectives: genome-wise (GBC), per chromosomes (CBC), and per chromosome segments (SBC). The former (GBC) depicted an overall picture of the “mosaic” genome of the Brangus attributable to their ancestors, whereas the latter two criteria (CBC and SBC) corresponded to local ancestral contributions. The average GBC for the 3,605 Brangus cattle were 70.2% Angus and 29.8% Brahman. The K-means clustering supported the postulation of the mixture of 1/2 Ultrablack (UB) animals in Brangus. For the non-UB Brangus animals, the average GBC were estimated to be 67.4% Angus and 32.6% Brahman. The 95% confidence intervals of their overall GBC were 60.4%–73.5% Angus and 26.5%–39.6% Brahman. Possibly, genetic selection and drifting have resulted in an approximately 5% average deviation toward Angus lineage. The estimated ancestral contributions by chromosomes were heavily distributed toward Angus, with 27 chromosomes having an average Angus CBC greater than 62.5% but only two chromosomes (5 and 20) having Brahman CBC greater than 37.5%. The chromosomal regions with high Angus breed proportions were prevalent, tending to form larger blocks on most chromosomes. In contrast, chromosome segments with high Brahman breed proportion were relatively few and isolated, presenting only on seven chromosomes. Hence, genomic hitchhiking effects were strong where Angus favorable alleles resided but weak where Brahman favorable alleles were present. The functions of genes identified in the chromosomal regions with high (≥75%) Angus compositions were diverse yet may were related to growth and body development. In contrast, the genes identified in the regions with high (≥37.5%) Brahman compositions were primarily responsible for disease resistance. In conclusion, we have addressed the questions concerning the Brangus genetic make-ups. The results can help form a dynamic picture of the Brangus breed formation and the genomic reshaping.

Genomic diversity and relationship analyses of endangered German Black Pied cattle (DSN) to 68 other taurine breeds based on whole-genome sequencing

German Black Pied cattle (Deutsches Schwarzbuntes Niederungsrind, DSN) are an endangered dual-purpose cattle breed originating from the North Sea region. The population comprises about 2,500 cattle and is considered one of the ancestral populations of the modern Holstein breed. The current study aimed at defining the breeds closest related to DSN cattle, characterizing their genomic diversity and inbreeding. In addition, the detection of selection signatures between DSN and Holstein was a goal. Relationship analyses using fixation index (F_ST), phylogenetic, and admixture analyses were performed between DSN and 68 other breeds from the 1000 Bull Genomes Project. Nucleotide diversity, observed heterozygosity, and expected heterozygosity were calculated as metrics for genomic diversity. Inbreeding was measured as excess of homozygosity (F_Hom) and genomic inbreeding (F_RoH) through runs of homozygosity (RoHs). Region-wide F_ST and cross-population-extended haplotype homozygosity (XP-EHH) between DSN and Holstein were used to detect selection signatures between the two breeds, and RoH islands were used to detect selection signatures within DSN and Holstein. DSN showed a close genetic relationship with breeds from the Netherlands, Belgium, Northern Germany, and Scandinavia, such as Dutch Friesian Red, Dutch Improved Red, Belgian Red White Campine, Red White Dual Purpose, Modern Angler, Modern Danish Red, and Holstein. The nucleotide diversity in DSN (0.151%) was higher than in Holstein (0.147%) and other breeds, e.g., Norwegian Red (0.149%), Red White Dual Purpose (0.149%), Swedish Red (0.149%), Hereford (0.145%), Angus (0.143%), and Jersey (0.136%). The F_Hom and F_RoH values in DSN were among the lowest. Regions with high F_ST between DSN and Holstein, significant XP-EHH regions, and RoH islands detected in both breeds harbor candidate genes that were previously reported for milk, meat, fertility, production, and health traits, including one QTL detected in DSN for endoparasite infection resistance. The selection signatures between DSN and Holstein provide evidence of regions responsible for the dual-purpose properties of DSN and the milk type of Holstein. Despite the small population size, DSN has a high level of diversity and low inbreeding. F_ST supports its relatedness to breeds from the same geographic origin and provides information on potential gene pools that could be used to maintain diversity in DSN.

Feedlot growth performance and carcass characteristics of steers backgrounded on buffel grass or buffel-Desmanthus mixed pastures

Feedlot performance and carcass characteristics of tropical beef steers backgrounded on buffel grass (Cenchrus ciliaris) only or buffel grass oversown with desmanthus (Desmanthus spp. ; 11.5% initial sward botanical composition) were evaluated. It was hypothesized that tropical beef cattle steers backgrounded on buffel grass only or buffel grass oversown with desmanthus with similar backgrounding growth performance will not differ in feedlot growth performance and carcass quality. Three hundred and twelve Bos indicus × Bos taurus tropical composite steers, 20-23 months old and weighing 413 ± 24 kg, previously backgrounded on buffel grass only or buffel-desmanthus mixed pastures for 147 days were finished on a concentrate diet in the feedlot for 110 days before slaughter. Buffel-desmanthus backgrounded steers had a slightly higher average daily gain (ADG; 1.8 kg/day) than the buffel grass backgrounded steers that had 1.7 kg/day ADG (p < 0.01). However, the final live weight and dry matter intake were not different (p ≥ 0.59). All the carcass traits measured were not different (p ≥ 0.18). Only 4% buffel grass and 8% buffel-desmanthus backgrounded steers fell short of the Meat Standards Australia (MSA) index, a level that is within the 4-9% reported for cattle produced in Queensland and slaughtered between July 2019 and June 2020. These findings indicate that desmanthus can be used to background beef cattle in northern Australia vertosol soil regions, where there is a paucity of adapted pasture legumes, with no negative impact on feedlot performance and carcass quality. The hypothesis that tropical beef cattle steers backgrounded on buffel grass only pastures or buffel grass oversown with desmanthus with similar backgrounding growth performance will have similar feedlot growth performance and carcass quality was accepted.

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 (&lt;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.

Evaluation of the association between plasma glucose-dependent insulinotropic polypeptide, respiratory quotient, and intramuscular fat deposition in feedlot cattle fed different levels of dry matter intake

The objectives of this trial were to evaluate the association between different levels of dry matter intake (DMI) on gas exchange, plasma glucose-dependent insulinotropic polypeptide (GIP) concentration, and intramuscular (IM) fat deposition. We used 60 individually fed backgrounded Angus × SimAngus-crossbred steers (n = 30) in a randomized complete block design. Steers (paired by body weight [BW] and gain to feed ratio [G:F]) were randomly allocated to one of the following treatments: ad libitum intake (AI) or restricted intake (RI; the same diet fed at 85% of the AI) of a finishing diet. The diet contained 61% cracked corn, 9% corn silage, 15% distillers’ dried grains with solubles, 5% soyhulls, and 10% of a protein-mineral-vitamin premix. Measurements of CO₂ emission and consumption of O₂, and respiratory quotient (RQ) were taken using the GreenFeed system (n = 15/treatment). Plasma and gas samples were collected 10 d before slaughter, 1 h before and 2 h after feeding. Plasma glucose, non-esterified fatty acids, GIP, and insulin concentration and gasses (O_2, CO_2, and RQ) were analyzed using the MIXED procedure of SAS evaluating the fixed effect of treatment, time (repeated measurement) and their interaction, and the random effect of the block. Final BW and carcass characteristics were analyzed with a similar model, without the time statement and its interaction. Compared with RI, AI steers had greater (P < 0.01) DMI and average daily gain (ADG). Steers on AI had greater final BW (P = 0.02), tended to have a greater ribeye area (P = 0.09), and had lower plasma GIP concentration (P = 0.04). There was no treatment effect (P ≥ 0.11) on G:F, subcutaneous backfat (BF), and IM fat, O₂ consumption, CO₂ emission, and RQ. Plasma glucose concentration of AI steers was greater before and after feeding than RI (P < 0.05). In conclusion, feeding steers ad libitum increased DMI, ADG, and plasma glucose and GIP concentration but does not affect G:F, BF, IM fat, CO₂ emission, and O₂ consumption. Plasma GIP concentration and RQ are not associated with IM fat deposition.

Genome-Wide Selection Signatures and Human-Mediated Introgression Events in Bos taurus indicus-influenced Composite Beef Cattle

Genetic introgression from interbreeding hybridization of European Bos taurus taurus (EBT) and Indian Bos taurus indicus (IBI) cattle breeds have been widely used to combine the climatic resilience of the IBI cattle and the higher productivity of EBT when forming new composite beef cattle (CB) populations. The subsequent breeding strategies have shifted their initial genomic compositions. To uncover population structure, signatures of selection, and potential introgression events in CB populations, high-density genotypes [containing 492,954 single nucleotide polymorphisms (SNPs) after the quality control] of 486 individuals from 15 cattle breeds, including EBT, IBI, and CB populations, along with two Bos grunniens genotypes as outgroup were used in this study. Then, in-depth population genetics analyses were performed for three CB breeds of Beefmaster, Brangus, and Santa Gertrudis. Neighbor-joining, principal components, and admixture analyses confirmed the historical introgression of EBT and IBI haplotypes into CB breeds. The f_dM statistics revealed that only 12.9% of CB populations' genetic components are of IBI origin. The results of signatures of selection analysis indicated different patterns of selection signals in the three CB breeds with primary pressure on pathways involved in protein processing and stress response in Beefmaster, cell proliferation regulation and immune response in Brangus, and amino acids and glucose metabolisms in Santa Gertrudis. An average of >90% of genomic regions underlying selection signatures were of EBT origin in the studied CB populations. Investigating the CB breeds' genome allows the estimation of EBT and IBI ancestral proportions and the locations within the genome where either taurine or indicine origin alleles are under selective pressure. Such findings highlight various opportunities to control the selection process more efficiently and explore complementarity at the genomic level in CB populations.

Rhamnolipids Regulate Lipid Metabolism, Immune Response, and Gut Microbiota in Rats

Objectives

Gut microbes influence lipid metabolism and immune responses that are key features of metabolic disorders. This study examined effects of bacterial rhamnolipids (RLS) on lipid metabolism, immune response, and gut microbiota in rats.

Methods

Twenty-four Sprague-Dawley rats were randomly divided into three groups and gavage-fed for seven weeks with normal saline (NCO group), 50 mg/kg bw RLS (RLS1 group), and 100 mg/kg bw RLS (RLS2 group).

Results

Compared with those of the NCO group, the RLS1 and RLS2 groups showed significantly decreased fat weight, relative fat weight, and adipocyte size (P < 0.05). Furthermore, RLS1 and RLS2 significantly decreased concentrations of triglycerides, low-density lipoprotein-cholesterol, and non-esterified fatty acids and increased high-density lipoprotein-cholesterol levels (P < 0.05). However, the total cholesterol content among the three groups (P > 0.05) were not significantly different. Serum concentrations of interleukin-1β, interleukin-6, and tumor necrosis factor-α were significantly lower in the RLS2 group than those in the NCO group (P < 0.05). The relative mRNA expression of fatty acid synthase was significantly decreased, while those of carnitine palmitoyltransferase-1, carnitine palmitoyltransferase-2, and peroxisome proliferator-activated receptor-gamma coactivator-1α were significantly increased in the RLS2 group compared with those in the NCO group (P < 0.05). Moreover, the relative abundances of Lactobacillus, Roseburia, Ruminococcus-1, and Parabacteroides were significantly higher in the RLS2 group than those in the NCO group (P < 0.05).

Conclusion

Our findings suggest that RLS reduces fat deposition, inhibits inflammation, regulates intestinal flora, and promotes the proliferation of beneficial bacteria in rats.

CircUBE2Q2 promotes differentiation of cattle muscle stem cells and is a potential regulatory molecule of skeletal muscle development

Background

The growth and development of muscle stem cells (MuSCs) are significant events known to affect muscle plasticity, disease, meat production, and meat quality, which involves the types and functions of mRNA and non-coding RNA. Here, MuSCs were cultured from Guangxi fetal cattle. RNA sequencing was used to analyze the RNA expression of mRNA and non-coding RNAs during the cell proliferation and differentiation phases.

Results

Two thousand one hundred forty-eight mRNAs and 888 non-coding RNAs were differentially expressed between cell proliferation and differentiation phases, including 113 miRNAs, 662 lncRNAs, and 113 circRNAs. RT-qPCR verified the differential expression levels of mRNAs and non-coding RNAs, and the differentially expressed circUBE2Q2 was subsequently characterized. Expression profile analysis revealed that circUBE2Q2 was abundant in muscle tissues and intramuscular fat. The expression of cricUBE2Q2 was also significantly upregulated during MuSCs myogenic differentiation and SVFs adipogenic differentiation and decreased with age in cattle muscle tissue. Finally, the molecular mechanism of circUBE2Q2 regulating MuSCs function that affects skeletal muscle development was investigated. The results showed that circUBE2Q2 could serve as a sponge for miR-133a, significantly promoting differentiation and apoptosis of cultured MuSCs, and inhibiting proliferation of MuSCs.

Conclusions

CircUBE2Q2 is associated with muscle growth and development and induces MuSCs myogenic differentiation through sponging miR-133a. This study will provide new clues for the mechanisms by which mRNAs and non-coding RNAs regulate skeletal muscle growth and development, affecting muscle quality and diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08518-4.

Sustained Effects of Muscle Calpain System Genotypes on Tenderness Phenotypes of South African Beef Bulls during Ageing up to 20 Days

Simple Summary

When searching for genetic markers for the selection of more tender beef, it is important to maintain minimal environmental variation from pre-slaughter, right through to the ageing process, to ensure the accuracy of the obtained phenotypes. This is because beef quality traits have a large environmental component that can greatly alter the characteristics of the meat, which would not reflect a true genetic effect. We propose that variable ageing times are especially important in determining whether markers are associated with tenderization or not. Our analyses included candidate genes for the protein degrading enzyme system for calpains, because they contribute the most to tenderization. We were able to validate these markers in South African beef cattle, where they could be useful for selection. The timing of the collection of tenderness data was critical, as only a few (6/134) genetic markers sustained their association with tenderization over ageing to 20 days. A larger tenderization effect earlier in ageing, as shown here for the capn1_187 and capn1_4751 markers, would decrease the length of ageing. This would not only increase profits, but also decrease the energy needed during the storage and refrigeration of aged beef, decreasing the carbon footprint of beef production.

Abstract

The most important factor that determines beef tenderness is its proteolytic activity, and the balance between calpain-1 protease activity and calpastatin inhibition is especially important, while contributions can also arise from calpain-2 and, possibly, calpain-3. The meat ageing process itself affects these processes. To determine whether genotypes in the calpain–calpastatin system can enhance tenderness through a 20-day ageing period, South African purebred beef bulls (n = 166) were genotyped using the Illumina BovineHD SNP BeadChip through a gene-based association analysis targeting the cast, capn3, capn2 and capn1 genes. The Warner–Bratzler shear force (WBSF) and myofibril fragment length (MFL) of Longissimus thoracis et lumborum (LTL) steaks were evaluated between d 3 and d 20 of ageing, with protease enzyme activity in the first 20 h post-mortem. Although several of the 134 SNPs are associated with tenderness, only seven SNP in the cast, capn2 and capn1 genes sustained genetic associations, additive to the ageing-associated increases in tenderness for at least three of the four ageing periods. While most genomic associations were relatively stable over time, some genotypes within the SNP responded differently to ageing, resulting in altered genomic effects over time. The level of ageing at which genomic associations are performed is an important factor that determines whether SNPs affect tenderness phenotypes.

Genome Wide Association Study of Beef Traits in Local Alpine Breed Reveals the Diversity of the Pathways Involved and the Role of Time Stratification

Knowledge of the genetic architecture of key growth and beef traits in livestock species has greatly improved worldwide thanks to genome-wide association studies (GWAS), which allow to link target phenotypes to Single Nucleotide Polymorphisms (SNPs) across the genome. Local dual-purpose breeds have rarely been the focus of such studies; recently, however, their value as a possible alternative to intensively farmed breeds has become clear, especially for their greater adaptability to environmental change and potential for survival in less productive areas. We performed single-step GWAS and post-GWAS analysis for body weight (BW), average daily gain (ADG), carcass fleshiness (CF) and dressing percentage (DP) in 1,690 individuals of local alpine cattle breed, Rendena. This breed is typical of alpine pastures, with a marked dual-purpose attitude and good genetic diversity. Moreover, we considered two of the target phenotypes (BW and ADG) at different times in the individuals' life, a potentially important aspect in the study of the traits' genetic architecture. We identified 8 significant and 47 suggestively associated SNPs, located in 14 autosomal chromosomes (BTA). Among the strongest signals, 3 significant and 16 suggestive SNPs were associated with ADG and were located on BTA10 (50-60 Mb), while the hotspot associated with CF and DP was on BTA18 (55-62 MB). Among the significant SNPs some were mapped within genes, such as SLC12A1, CGNL1, PRTG (ADG), LOC513941 (CF), NLRP2 (CF and DP), CDC155 (DP). Pathway analysis showed great diversity in the biological pathways linked to the different traits; several were associated with neurogenesis and synaptic transmission, but actin-related and transmembrane transport pathways were also represented. Time-stratification highlighted how the genetic architectures of the same traits were markedly different between different ages. The results from our GWAS of beef traits in Rendena led to the detection of a variety of genes both well-known and novel. We argue that our results show that expanding genomic research to local breeds can reveal hitherto undetected genetic architectures in livestock worldwide. This could greatly help efforts to map genomic complexity of the traits of interest and to make appropriate breeding decisions.

Genetic parameters for traditional and novel ultrasound carcass traits in Nellore cattle

The objective of this study was to contribute to genetic improvement for carcass quality in Nellore cattle by estimating genetic parameters for carcass characteristics measured by traditional and novel ultrasound adjusted for body weight. Records of 19,910 Nellore breed animals born between 2001 and 2015 as part of a genetic improvement program were evaluated for this purpose. The (co)variance components and genetic parameters for yearling weight (YW), slaughter conformation at yearling age (SCY), ribeye area (REA), REA adjusted for 100 kg of body weight (REA100), REA adjusted for 450 kg of body weight (REA450), subcutaneous fat thickness over the loin (SFTL) and the rump (SFTR), intramuscular fat (IF), and height × width ratio of the sirloin steak (RATIO) were estimated using the restricted maximum likelihood (REML) method and the three-trait animal model. The estimated heritability values were 0.39 (YW), 0.33 (SCY), 0.31 (REA), 0.42 (REA100), 0.38 (REA450), 0.31 (RATIO), 0.17 (SFTL), 0.27 (SFTR), and 0.41 (MAR). The genetic correlations between REA and YW (0.43) and SCY and YW (0.78), respectively, suggested that selecting for REA or SCY could result in an increase in body weight. Because of its moderate and negative correlation with YW, REA100 was ineffective as a muscling indicator (- 0.38). The indicators REA450 and RATIO had a - 0.12 and 0.11 low correlation with YW, respectively. Thus, it is concluded that REA450 and RATIO were the best indicators of muscling, with heritability of medium magnitude, making the selection viable, as well as having a low correlation with body weight.

Genome-wide Association Study for Carcass Primal Cut Yields Using Single-step Bayesian Approach in Hanwoo Cattle

The importance of meat and carcass quality is growing in beef cattle production to meet both producer and consumer demands. Primal cut yields, which reflect the body compositions of carcass, could determine the carcass grade and, consequently, command premium prices. Despite its importance, there have been few genome-wide association studies on these traits. This study aimed to identify genomic regions and putative candidate genes related to 10 primal cut traits, including tenderloin, sirloin, striploin, chuck, brisket, top round, bottom round, shank, flank, and rib in Hanwoo cattle using a single-step Bayesian regression (ssBR) approach. After genomic data quality control, 43,987 SNPs from 3,745 genotyped animals were available, of which 3,467 had phenotypic records for the analyzed traits. A total of 16 significant genomic regions (1-Mb window) were identified, of which five large-effect quantitative trait loci (QTLs) located on chromosomes 6 at 38–39 Mb, 11 at 21–22 Mb, 14 at 6–7 Mb and 26–27 Mb, and 19 at 26–27 Mb were associated with more than one trait, while the remaining 11 QTLs were trait-specific. These significant regions were harbored by 154 genes, among which TOX, FAM184B, SPP1, IBSP, PKD2, SDCBP, PIGY, LCORL, NCAPG, and ABCG2 were noteworthy. Enrichment analysis revealed biological processes and functional terms involved in growth and lipid metabolism, such as growth (GO:0040007), muscle structure development (GO:0061061), skeletal system development (GO:0001501), animal organ development (GO:0048513), lipid metabolic process (GO:0006629), response to lipid (GO:0033993), metabolic pathways (bta01100), focal adhesion (bta04510), ECM–receptor interaction (bta04512), fat digestion and absorption (bta04975), and Rap1 signaling pathway (bta04015) being the most significant for the carcass primal cut traits. Thus, identification of quantitative trait loci regions and plausible candidate genes will aid in a better understanding of the genetic and biological mechanisms regulating carcass primal cut yields.

Genetic and Genomic Analyses of Service Sire Effect on Female Reproductive Traits in Holstein Cattle

Fertility and reproductive performance are key drivers of dairy farm profitability. Hence, reproduction traits have been included in a large majority of worldwide dairy cattle selection indexes. The reproductive traits are lowly heritable but can be improved through direct genetic selection. However, most scientific studies and dairy cattle breeding programs have focused solely on the genetic effects of the dam (GED) on reproductive performance and, therefore, ignored the contribution of the service sire in the phenotypic outcomes. This study aimed to investigate the service sire effects on female reproductive traits in Holstein cattle from a genomic perspective. Genetic parameter estimation and genome-wide association studies (GWAS) were performed for the genetic effect of service sire (GESS) on conception rate (CR), 56-day non-return rate (NRR56), calving ease (CE), stillbirth (SB), and gestation length (GL). Our findings indicate that the additive genetic effects of both sire and dam contribute to the phenotypic variance of reproductive traits measured in females (0.0196 vs. 0.0109, 0.0237 vs. 0.0133, 0.0040 vs. 0.0289, 0.0782 vs. 0.0083, and 0.1024 vs. 0.1020 for GESS and GED heritability estimates for CR, NRR56, CE, SB, and GL, respectively), and these two genetic effects are positively correlated for SB (0.1394) and GL (0.7871). Interestingly, the breeding values for GESS on insemination success traits (CR and NRR56) are unfavorably and significantly correlated with some production, health, and type breeding values (ranging from -0.449 to 0.274), while the GESS values on calving traits (CE, SB, and GL) are usually favorably associated with those traits (ranging from -0.493 to 0.313). One hundred sixty-two significant single-nucleotide polymorphisms (SNPs) and their surrounding protein-coding genes were identified as significantly associated with GESS and GED, respectively. Six genes overlapped between GESS and GED for calving traits and 10 genes overlapped between GESS for success traits and calving traits. Our findings indicate the importance of considering the GESS when genetically evaluating the female reproductive traits in Holstein cattle.

The Expression Profiles of mRNAs and lncRNAs in Buffalo Muscle Stem Cells Driving Myogenic Differentiation

Buffalo breeding has become an important branch of the beef cattle industry. Hence, it is of great significance to study buffalo meat production and meat quality. However, the expression profiles of mRNA and long non-coding RNAs (lncRNA) molecules in muscle stem cells (MuSCs) development in buffalo have not been explored fully. We, therefore, performed mRNA and lncRNA expression profiling analysis during the proliferation and differentiation phases of MuSCs in buffalo. The results showed that there were 4,820 differentially expressed genes as well as 12,227 mRNAs and 1,352 lncRNAs. These genes were shown to be enriched in essential biological processes such as cell cycle, p53 signaling pathway, RNA transport and calcium signaling pathway. We also identified a number of functionally important genes, such as MCMC4, SERDINE1, ISLR, LOC102394806, and LOC102403551, and found that interference with MYLPF expression significantly inhibited the differentiation of MuSCs. In conclusion, our research revealed the characteristics of mRNA and lncRNA expression during the differentiation of buffalo MuSCs. This study can be used as an important reference for the study of RNA regulation during muscle development in buffalo.

Integration of a single-step genome-wide association study with a multi-tissue transcriptome analysis provides novel insights into the genetic basis of wool and weight traits in sheep

BACKGROUND

Genetic improvement of wool and growth traits is a major goal in the sheep industry, but their underlying genetic architecture remains elusive. To improve our understanding of these mechanisms, we conducted a weighted single-step genome-wide association study (WssGWAS) and then integrated the results with large-scale transcriptome data for five wool traits and one growth trait in Merino sheep: mean fibre diameter (MFD), coefficient of variation of the fibre diameter (CVFD), crimp number (CN), mean staple length (MSL), greasy fleece weight (GFW), and live weight (LW).

RESULTS

Our dataset comprised 7135 individuals with phenotype data, among which 1217 had high-density (HD) genotype data (n = 372,534). The genotypes of 707 of these animals were imputed from the Illumina Ovine single nucleotide polymorphism (SNP) 54 BeadChip to the HD Array. The heritability of these traits ranged from 0.05 (CVFD) to 0.36 (MFD), and between-trait genetic correlations ranged from - 0.44 (CN vs. LW) to 0.77 (GFW vs. LW). By integrating the GWAS signals with RNA-seq data from 500 samples (representing 87 tissue types from 16 animals), we detected tissues that were relevant to each of the six traits, e.g. liver, muscle and the gastrointestinal (GI) tract were the most relevant tissues for LW, and leukocytes and macrophages were the most relevant cells for CN. For the six traits, 54 quantitative trait loci (QTL) were identified covering 81 candidate genes on 21 ovine autosomes. Multiple candidate genes showed strong tissue-specific expression, e.g. BNC1 (associated with MFD) and CHRNB1 (LW) were specifically expressed in skin and muscle, respectively. By conducting phenome-wide association studies (PheWAS) in humans, we found that orthologues of several of these candidate genes were significantly (FDR < 0.05) associated with similar traits in humans, e.g. BNC1 was significantly associated with MFD in sheep and with hair colour in humans, and CHRNB1 was significantly associated with LW in sheep and with body mass index in humans.

CONCLUSIONS

Our findings provide novel insights into the biological and genetic mechanisms underlying wool and growth traits, and thus will contribute to the genetic improvement and gene mapping of complex traits in sheep.

Genetic Analysis of Major Carcass Traits of Korean Hanwoo Males Raised for Thirty Months

Simple Summary

Generally, Korean Hanwoo males produced under a 24-month production cycle (PROD24) are evaluated as a part of the progeny test program. However, there is little information on other males outside the PROD24, such as those raised under a 30-month production cycle (PROD30) for higher profits. Therefore, we investigated PROD30 males for important carcass traits (carcass weight, eye muscle area, backfat thickness, and marbling score) using a reasonably large dataset to understand their genetic merit. To do so, we estimated the genetic parameters of traits using animal model. Our analysis revealed moderate to high heritability values for the studied traits. The marbling score was found to be highly heritable at 0.56. The genetic correlation between traits was mostly moderate to low, and the backfat thickness was poorly correlated with the marbling score. These results are consistent with many previous reports on PROD24. Our study suggests that PROD30 and PROD24 males might have somewhat similar genetic potential, as well as similar genetic backgrounds. Thus, it could be concluded that there is further scope for PROD30 males to improve Hanwoo males’ overall prediction accuracy, especially under a genomic selection program, together with PROD24 males.

Abstract

Understanding animals’ genetic potential for carcass traits is the key to genetic improvements in any beef cattle. In this study, we investigated the genetic merits of carcass traits using Hanwoo males raised in a 30-month production system (PROD30). We achieved this using a dataset comprising 6092 Hanwoo males born between 2005 and 2017 and measures of four carcass traits (carcass weight, CWT; eye muscle area, EMA; backfat thickness, BFT; and marbling score, MS). Genetic parameters were estimated using a multiple-trait animal model through the AIREMLF90 program. According to the multiple-trait model, the h² of CWT, EMA, BFT, and MS were 0.35 ± 0.04, 0.43 ± 0.05, 0.48 ± 0.05, and 0.56 ± 0.05, respectively. The strongest genetic correlation (r_g) was obtained between CWT and EMA (0.49 ± 0.07), whereas it was negligible between CWT and BFT. EMA and MS were also moderately correlated, whereas there was a relatively low negative correlation between EMA and BFT (−0.26 ± 0.08). Our study revealed a consistent indirect genetic improvement in animals from 2005 onwards. Although Hanwoo improvement has mainly focused on males under a 24-month production cycle, we observed PROD30 males to have somewhat similar genetic potential. Our results provide useful insights into the genetic merits of PROD30 males for the first time, which may facilitate future studies on them and their integration into the Hanwoo National Evaluation for genomic selection.

Genomic Approaches Reveal Pleiotropic Effects in Crossbred Beef Cattle

Carcass and meat quality are two important attributes for the beef industry because they drive profitability and consumer demand. These traits are of even greater importance in crossbred cattle used in subtropical and tropical regions for their superior adaptability because they tend to underperform compared to their purebred counterparts. Many of these traits are challenging and expensive to measure and unavailable until late in life or after the animal is harvested, hence unrealistic to improve through traditional phenotypic selection, but perfect candidates for genomic selection. Before genomic selection can be implemented in crossbred populations, it is important to explore if pleiotropic effects exist between carcass and meat quality traits. Therefore, the objective of this study was to identify genomic regions with pleiotropic effects on carcass and meat quality traits in a multibreed Angus-Brahman population that included purebred and crossbred animals. Data included phenotypes for 10 carcass and meat quality traits from 2,384 steers, of which 1,038 were genotyped with the GGP Bovine F-250. Single-trait genome-wide association studies were first used to investigate the relevance of direct additive genetic effects on each carcass, sensory and visual meat quality traits. A second analysis for each trait included all other phenotypes as covariates to correct for direct causal effects from identified genomic regions with pure direct effects on the trait under analysis. Five genomic windows on chromosomes BTA5, BTA7, BTA18, and BTA29 explained more than 1% of additive genetic variance of two or more traits. Moreover, three suggestive pleiotropic regions were identified on BTA10 and BTA19. The 317 genes uncovered in pleiotropic regions included anchoring and cytoskeletal proteins, key players in cell growth, muscle development, lipid metabolism and fat deposition, and important factors in muscle proteolysis. A functional analysis of these genes revealed GO terms directly related to carcass quality, meat quality, and tenderness in beef cattle, including calcium-related processes, cell signaling, and modulation of cell-cell adhesion. These results contribute with novel information about the complex genetic architecture and pleiotropic effects of carcass and meat quality traits in crossbred beef cattle.

Single-Nucleotide Polymorphism Variations Associated With Specific Genes Putatively Identified Enhanced Genetic Predisposition for 305-Day Milk Yield in the Girolando Crossbreed

Milk production phenotypes are the main focus of genetic selection in dairy herds, and although there are many genes identified as related to the biology of these traits in pure breeds, little is known about crossbreed animals. This study aimed to identify potential genes associated with the 305-day milk yield in 337 crossbreed Gir × Holstein (Girolando) animals. Milk production records were genotyped for 45,613 single-nucleotide polymorphisms (SNPs). This dataset was used for a genome-wide association study (GWAS) using the 305-day milk yield adjusted for the fixed effects of herd and year and linear and quadratic effects of age at calving (in days) and calving factor averaged per animal. Genes within the significant SNPs were retrieved from the Bos taurus ARS-UCD1.2 assembly (bosTau9) for gene ontology analysis. In summary, the GWAS identified 52 SNPs associated [p ≤ 10^–4, false discovery rate (FDR) = 8.77%] with milk production, including NUB1 and SLC24A2, which were previously described as related to milk production traits in cattle. The results suggest that SNPs associated mainly with NUB1 and SLC24A2 could be useful to understand milk production in Girolando and used as predictive markers for selecting genetic predisposition for milk yield in Girolando.

Signatures of Selection in Composite Vrindavani Cattle of India

Vrindavani is an Indian composite cattle breed developed by crossbreeding taurine dairy breeds with native indicine cattle. The constituent breeds were selected for higher milk production and adaptation to the tropical climate. However, the selection response for production and adaptation traits in the Vrindavani genome is not explored. In this study, we provide the first overview of the selection signatures in the Vrindavani genome. A total of 96 Vrindavani cattle were genotyped using the BovineSNP50 BeadChip and the SNP genotype data of its constituent breeds were collected from a public database. Within-breed selection signatures in Vrindavani were investigated using the integrated haplotype score (iHS). The Vrindavani breed was also compared to each of its parental breeds to discover between-population signatures of selection using two approaches, cross-population extended haplotype homozygosity (XP-EHH) and fixation index (F_ST). We identified 11 common regions detected by more than one method harboring genes such as LRP1B, TNNI3K, APOB, CACNA2D1, FAM110B, and SPATA17 associated with production and adaptation. Overall, our results suggested stronger selective pressure on regions responsible for adaptation compared to milk yield.

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.

Selection signatures in two oldest Russian native cattle breeds revealed using high-density single nucleotide polymorphism analysis

Native cattle breeds can carry specific signatures of selection reflecting their adaptation to the local environmental conditions and response to the breeding strategy used. In this study, we comprehensively analysed high-density single nucleotide polymorphism (SNP) genotypes to characterise the population structure and detect the selection signatures in Russian native Yaroslavl and Kholmogor dairy cattle breeds, which have been little influenced by introgression with transboundary breeds. Fifty-six samples of pedigree-recorded purebred animals, originating from different breeding farms and representing different sire lines, of the two studied breeds were genotyped using a genome-wide bovine genotyping array (Bovine HD BeadChip). Three statistical analyses—calculation of fixation index (F_ST) for each SNP for the comparison of the pairs of breeds, hapFLK analysis, and estimation of the runs of homozygosity (ROH) islands shared in more than 50% of animals—were combined for detecting the selection signatures in the genome of the studied cattle breeds. We confirmed nine and six known regions under putative selection in the genomes of Yaroslavl and Kholmogor cattle, respectively; the flanking positions of most of these regions were elucidated. Only two of the selected regions (localised on BTA 14 at 24.4–25.1 Mbp and on BTA 16 at 42.5–43.5 Mb) overlapped in Yaroslavl, Kholmogor and Holstein breeds. In addition, we detected three novel selection sweeps in the genome of Yaroslavl (BTA 4 at 4.74–5.36 Mbp, BTA 15 at 17.80–18.77 Mbp, and BTA 17 at 45.59–45.61 Mbp) and Kholmogor breeds (BTA 12 at 82.40–81.69 Mbp, BTA 15 at 16.04–16.62 Mbp, and BTA 18 at 0.19–1.46 Mbp) by using at least two of the above-mentioned methods. We expanded the list of candidate genes associated with the selected genomic regions and performed their functional annotation. We discussed the possible involvement of the identified candidate genes in artificial selection in connection with the origin and development of the breeds. Our findings on the Yaroslavl and Kholmogor breeds obtained using high-density SNP genotyping and three different statistical methods allowed the detection of novel putative genomic regions and candidate genes that might be under selection. These results might be useful for the sustainable development and conservation of these two oldest Russian native cattle breeds.

Genome-wide scan for selection signatures reveals novel insights into the adaptive capacity in local North African cattle

Natural-driven selection is supposed to have left detectable signatures on the genome of North African cattle which are often characterized by the fixation of genetic variants associated with traits under selection pressure and/or an outstanding genetic differentiation with other populations at particular loci. Here, we investigate the population genetic structure and we provide a first outline of potential selection signatures in North African cattle using single nucleotide polymorphism genotyping data. After comparing our data to African, European and indicine cattle populations, we identified 36 genomic regions using three extended haplotype homozygosity statistics and 92 outlier markers based on Bayescan test. The 13 outlier windows detected by at least two approaches, harboured genes (e.g. GH1, ACE, ASIC3, HSPH1, MVD, BCL2, HIGD2A, CBFA2T3) that may be involved in physiological adaptations required to cope with environmental stressors that are typical of the North African area such as infectious diseases, extended drought periods, scarce food supply, oxygen scarcity in the mountainous areas and high-intensity solar radiation. Our data also point to candidate genes involved in transcriptional regulation suggesting that regulatory elements had also a prominent role in North African cattle response to environmental constraints. Our study yields novel insights into the unique adaptive capacity in these endangered populations emphasizing the need for the use of whole genome sequence data to gain a better understanding of the underlying molecular mechanisms.

Transcriptional regulation of adipogenic marker genes for the improvement of intramuscular fat in Qinchuan beef cattle

The intramuscular fat content plays a crucial role in meat quality traits. Increasing the degree of adipogenesis in beef cattle leads to an increase in the content of intramuscular fat. Adipogenesis a complex biochemical process which is under firm genetic control. Over the last three decades, the Qinchuan beef cattle have been extensively studied for the improvement of meat production and quality traits. In this study, we reviewed the literature regarding adipogenesis and intramuscular fat deposition. Then, we summarized the research conducted on the transcriptional regulation of key adipogenic marker genes, and also reviewed the roles of adipogenic marker genes in adipogenesis of Qinchuan beef cattle. This review will elaborate our understanding regarding transcriptional regulation which is a vital physiological process regulated by a cascade of transcription factors (TFs), key target marker genes, and regulatory proteins. This synergistic action of TFs and target genes ensures the accurate and diverse transmission of the genetic information for the accomplishment of central physiological processes. This information will provide an insight into the transcriptional regulation of the adipogenic marker genes and its role in bovine adipogenesis for the breed improvement programs especially for the trait of intramuscular fat deposition.

Fifteen Shades of Grey: Combined Analysis of Genome-Wide SNP Data in Steppe and Mediterranean Grey Cattle Sheds New Light on the Molecular Basis of Coat Color

Coat color is among the most distinctive phenotypes in cattle. Worldwide, several breeds share peculiar coat color features such as the presence of a fawn pigmentation of the calf at birth, turning over time to grey, and sexual dichromatism. The aim of this study was to search for polymorphisms under differential selection by contrasting grey cattle breeds displaying the above phenotype with non-grey cattle breeds, and to identify the underlying genes. Using medium-density SNP array genotype data, a multi-cohort F_ST-outlier approach was adopted for a total of 60 pair-wise comparisons of the 15 grey with 4 non-grey cattle breeds (Angus, Limousin, Charolais, and Holstein), with the latter selected as representative of solid and piebald phenotypes, respectively. Overall, more than 50 candidate genes were detected; almost all were either directly or indirectly involved in pigmentation, and some of them were already known for their role in phenotypes related with hair graying in mammals. Notably, 17 relevant genes, including SDR16C5, MOS, SDCBP, and NSMAF, were located in a signal on BTA14 convergently observed in all the four considered scenarios. Overall, the key stages of pigmentation (melanocyte development, melanogenesis, and pigment trafficking/transfer) were all represented among the pleiotropic functions of the candidate genes, suggesting the complex nature of the grey phenotype in cattle.

Genetic Parameters and Genome-Wide Association Studies for Anti-Müllerian Hormone Levels and Antral Follicle Populations Measured After Estrus Synchronization in Nellore Cattle

Reproductive efficiency plays a major role in the long-term sustainability of livestock industries and can be improved through genetic and genomic selection. This study aimed to estimate genetic parameters (heritability and genetic correlation) and identify genomic regions and candidate genes associated with anti-Müllerian hormone levels (AMH) and antral follicle populations measured after estrous synchronization (AFP) in Nellore cattle. The datasets included phenotypic records for 1099 and 289 Nellore females for AFP and AMH, respectively, high-density single nucleotide polymorphism (SNP) genotypes for 944 animals, and 4129 individuals in the pedigree. The heritability estimates for AMH and AFP were 0.28 ± 0.07 and 0.30 ± 0.09, and the traits were highly and positively genetically correlated (rG = 0.81 ± 0.02). These findings indicated that these traits can be improved through selective breeding, and substantial indirect genetic gains are expected by selecting for only one of the two traits. A total of 31 genomic regions were shown to be associated with AMH or AFP, and two genomic regions located on BTA1 (64.9-65.0 Mb and 109.1-109.2 Mb) overlapped between the traits. Various candidate genes were identified to be potentially linked to important biological processes such as ovulation, tissue remodeling, and the immune system. Our findings support the use of AMH and AFP as indicator traits to genetically improve fertility rates in Nellore cattle and identify better oocyte donors.