Meta-analysis across Nellore cattle populations identifies common metabolic mechanisms that regulate feed efficiency-related traits

Differential gene expression in neonatal calf muscle tissues from Hanwoo cows overfed during mid to late pregnancy period

Maternal nutrition significantly influences fetal development and postnatal outcomes. This study investigates the impact of maternal overfeeding during mid to late pregnancy on gene expression in the round and sirloin muscles of Hanwoo neonatal calves. Eight cows were assigned to either a control group receiving standard nutrition (100%) or a treated group receiving overnutrition (150%). After birth, tissue samples from the round and sirloin muscles of neonatal calves were collected and subjected to RNA sequencing to assess differentially expressed genes (DEGs). RNA sequencing identified 43 DEGs in round muscle and 15 in sirloin muscle, involving genes related to myogenesis, adipogenesis, and energy regulation. Key genes, including PPARGC1A, THBS1, CD44, JUND, CNN1, ENAH, and RUNX1, were predominantly downregulated. Gene ontology (GO) enrichment analyses revealed terms associated with muscle development, such as "biological regulation," "cellular process," and "response to stimulus." Protein-protein interaction networks highlighted complex interactions among DEGs. Random Forest analysis identified ARC, SLC1A5, and GNPTAB as influential genes for distinguishing between control and treated groups. Overall, maternal overnutrition during mid-to-late pregnancy results in the downregulation of genes involved in muscle development and energy metabolism in neonatal Hanwoo calves. These findings provide insights into the molecular effects of maternal nutrition on muscle development.

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.

Meta-analysis of six dairy cattle breeds reveals biologically relevant candidate genes for mastitis resistance

BACKGROUND

Mastitis is a disease that incurs significant costs in the dairy industry. A promising approach to mitigate its negative effects is to genetically improve the resistance of dairy cattle to mastitis. A meta-analysis of genome-wide association studies (GWAS) across multiple breeds for clinical mastitis (CM) and its indicator trait, somatic cell score (SCS), is a powerful method to identify functional genetic variants that impact mastitis resistance.

RESULTS

We conducted meta-analyses of eight and fourteen GWAS on CM and SCS, respectively, using 30,689 and 119,438 animals from six dairy cattle breeds. Methods for the meta-analyses were selected to properly account for the multi-breed structure of the GWAS data. Our study revealed 58 lead markers that were associated with mastitis incidence, including 16 loci that did not overlap with previously identified quantitative trait loci (QTL), as curated at the Animal QTLdb. Post-GWAS analysis techniques such as gene-based analysis and genomic feature enrichment analysis enabled prioritization of 31 candidate genes and 14 credible candidate causal variants that affect mastitis.

CONCLUSIONS

Our list of candidate genes can help to elucidate the genetic architecture underlying mastitis resistance and provide better tools for the prevention or treatment of mastitis, ultimately contributing to more sustainable animal production.

Genome-Wide Association Study for Weight Loss at the End of Dry-Curing of Hams Produced from Purebred Heavy Pigs

Dissecting the genetics of production traits in livestock is of outmost importance, both to understand biological mechanisms underlying those traits and to facilitate the design of selection programs incorporating that information. For the pig industry, traits related to curing are key for protected designation of origin productions. In particular, appropriate ham weight loss after dry-curing ensures high quality of the final product and avoids economic losses. In this study, we analyzed data (N = 410) of ham weight loss after approximately 20 months of dry-curing. The animals used for ham production were purebred pigs belonging to a commercial line. A genome-wide association study (GWAS) of 29,844 SNP markers revealed the polygenic nature of the trait: 221 loci explaining a small percentage of the variance (0.3-1.65%) were identified on almost all Sus scrofa chromosomes. Post-GWAS analyses revealed 32 windows located within regulatory regions and 94 windows located in intronic regions of specific genes. In total, 30 candidate genes encoding receptors and enzymes associated with ham weight loss (MTHFD1L, DUSP8), proteolysis (SPARCL1, MYH8), drip loss (TNNI2), growth (CDCA3, LSP1, CSMD1, AP2A2, TSPAN4), and fat metabolism (AGPAT4, IGF2R, PTDSS2, HRAS, TALDO1, BRSK2, TNNI2, SYT8, GTF2I, GTF2IRD1, LPCAT3, ATN1, GNB3, CMIP, SORCS2, CCSER1, SPP1) were detected.

Genome-wide association study considering genotype-by-environment interaction for productive and reproductive traits using whole-genome sequencing in Nellore cattle

BACKGROUND

The genotype-by-environment interaction (GxE) in beef cattle can be investigated using reaction norm models to assess environmental sensitivity and, combined with genome-wide association studies (GWAS), to map genomic regions related to animal adaptation. Including genetic markers from whole-genome sequencing in reaction norm (RN) models allows us to identify high-resolution candidate genes across environmental gradients through GWAS. Hence, we performed a GWAS via the RN approach using whole-genome sequencing data, focusing on mapping candidate genes associated with the expression of reproductive and growth traits in Nellore cattle. For this purpose, we used phenotypic data for age at first calving (AFC), scrotal circumference (SC), post-weaning weight gain (PWG), and yearling weight (YW). A total of 20,000 males and 7,159 females genotyped with 770k were imputed to the whole sequence (29 M). After quality control and linkage disequilibrium (LD) pruning, there remained ∼ 2.41 M SNPs for SC, PWG, and YW and ∼ 5.06 M SNPs for AFC.

RESULTS

Significant SNPs were identified on Bos taurus autosomes (BTA) 10, 11, 14, 18, 19, 20, 21, 24, 25 and 27 for AFC and on BTA 4, 5 and 8 for SC. For growth traits, significant SNP markers were identified on BTA 3, 5 and 20 for YW and PWG. A total of 56 positional candidate genes were identified for AFC, 9 for SC, 3 for PWG, and 24 for YW. The significant SNPs detected for the reaction norm coefficients in Nellore cattle were found to be associated with growth, adaptative, and reproductive traits. These candidate genes are involved in biological mechanisms related to lipid metabolism, immune response, mitogen-activated protein kinase (MAPK) signaling pathway, and energy and phosphate metabolism.

CONCLUSIONS

GWAS results highlighted differences in the physiological processes linked to lipid metabolism, immune response, MAPK signaling pathway, and energy and phosphate metabolism, providing insights into how different environmental conditions interact with specific genes affecting animal adaptation, productivity, and reproductive performance. The shared genomic regions between the intercept and slope are directly implicated in the regulation of growth and reproductive traits in Nellore cattle raised under different environmental conditions.

Combining genetic markers, on-farm information and infrared data for the in-line prediction of blood biomarkers of metabolic disorders in Holstein cattle

BACKGROUND

Various blood metabolites are known to be useful indicators of health status in dairy cattle, but their routine assessment is time-consuming, expensive, and stressful for the cows at the herd level. Thus, we evaluated the effectiveness of combining in-line near infrared (NIR) milk spectra with on-farm (days in milk [DIM] and parity) and genetic markers for predicting blood metabolites in Holstein cattle. Data were obtained from 388 Holstein cows from a farm with an AfiLab system. NIR spectra, on-farm information, and single nucleotide polymorphisms (SNP) markers were blended to develop calibration equations for blood metabolites using the elastic net (ENet) approach, considering 3 models: (1) Model 1 (M1) including only NIR information, (2) Model 2 (M2) with both NIR and on-farm information, and (3) Model 3 (M3) combining NIR, on-farm and genomic information. Dimension reduction was considered for M3 by preselecting SNP markers from genome-wide association study (GWAS) results.

RESULTS

Results indicate that M2 improved the predictive ability by an average of 19% for energy-related metabolites (glucose, cholesterol, NEFA, BHB, urea, and creatinine), 20% for liver function/hepatic damage, 7% for inflammation/innate immunity, 24% for oxidative stress metabolites, and 23% for minerals compared to M1. Meanwhile, M3 further enhanced the predictive ability by 34% for energy-related metabolites, 32% for liver function/hepatic damage, 22% for inflammation/innate immunity, 42.1% for oxidative stress metabolites, and 41% for minerals, compared to M1. We found improved predictive ability of M3 using selected SNP markers from GWAS results using a threshold of > 2.0 by 5% for energy-related metabolites, 9% for liver function/hepatic damage, 8% for inflammation/innate immunity, 22% for oxidative stress metabolites, and 9% for minerals. Slight reductions were observed for phosphorus (2%), ferric-reducing antioxidant power (1%), and glucose (3%). Furthermore, it was found that prediction accuracies are influenced by using more restrictive thresholds (-log10(P-value) > 2.5 and 3.0), with a lower increase in the predictive ability.

CONCLUSION

Our results highlighted the potential of combining several sources of information, such as genetic markers, on-farm information, and in-line NIR infrared data improves the predictive ability of blood metabolites in dairy cattle, representing an effective strategy for large-scale in-line health monitoring in commercial herds.

Benchmarking machine learning and parametric methods for genomic prediction of feed efficiency-related traits in Nellore cattle

Genomic selection (GS) offers a promising opportunity for selecting more efficient animals to use consumed energy for maintenance and growth functions, impacting profitability and environmental sustainability. Here, we compared the prediction accuracy of multi-layer neural network (MLNN) and support vector regression (SVR) against single-trait (STGBLUP), multi-trait genomic best linear unbiased prediction (MTGBLUP), and Bayesian regression (BayesA, BayesB, BayesC, BRR, and BLasso) for feed efficiency (FE) traits. FE-related traits were measured in 1156 Nellore cattle from an experimental breeding program genotyped for ~ 300 K markers after quality control. Prediction accuracy (Acc) was evaluated using a forward validation splitting the dataset based on birth year, considering the phenotypes adjusted for the fixed effects and covariates as pseudo-phenotypes. The MLNN and SVR approaches were trained by randomly splitting the training population into fivefold to select the best hyperparameters. The results show that the machine learning methods (MLNN and SVR) and MTGBLUP outperformed STGBLUP and the Bayesian regression approaches, increasing the Acc by approximately 8.9%, 14.6%, and 13.7% using MLNN, SVR, and MTGBLUP, respectively. Acc for SVR and MTGBLUP were slightly different, ranging from 0.62 to 0.69 and 0.62 to 0.68, respectively, with empirically unbiased for both models (0.97 and 1.09). Our results indicated that SVR and MTGBLUBP approaches were more accurate in predicting FE-related traits than Bayesian regression and STGBLUP and seemed competitive for GS of complex phenotypes with various degrees of inheritance.

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.

Characterization of heterozygosity-rich regions in Italian and worldwide goat breeds

Heterozygosity-rich regions (HRR) are genomic regions of high heterozygosity, which may harbor loci related to key functional traits such as immune response, survival rate, fertility, and other fitness traits. This study considered 30 Italian and 19 worldwide goat breeds genotyped with the Illumina GoatSNP50k BeadChip. The aim of the work was to study inter-breed relationships and HRR patterns using Sliding Window (SW) and Consecutive Runs (CR) detection methods. Genetic relationships highlighted a clear separation between non-European and European breeds, as well as the north-south geographic cline within the latter. The Pearson correlation coefficients between the descriptive HRR parameters obtained with the SW and CR methods were higher than 0.9. A total of 166 HRR islands were detected. CHI1, CHI11, CHI12 and CHI18 were the chromosomes harboring the highest number of HRR islands. The genes annotated in the islands were linked to various factors such as productive, reproductive, immune, and environmental adaptation mechanisms. Notably, the Montecristo feral goat showed the highest number of HRR islands despite the high level of inbreeding, underlining potential balancing selection events characterizing its evolutionary history. Identifying a species-specific HRR pattern could provide a clearer view of the mechanisms regulating the genome modelling following anthropogenic selection combined with environmental interaction.

Genotype-by-environment interactions for feed efficiency traits in Nellore cattle based on bi-trait reaction norm models

BACKGROUND

Selecting animals for feed efficiency directly impacts the profitability of the beef cattle industry, which contributes to minimizing the environmental footprint of beef production. Genetic and environmental factors influence animal feed efficiency, leading to phenotypic variability when exposed to different environmental conditions (i.e., temperature and nutritional level). Thus, our aim was to assess potential genotype-by-environment (G × E) interactions for dry matter intake (DMI) and residual feed intake (RFI) in Nellore cattle (Bos taurus indicus) based on bi-trait reaction norm models (RN) and evaluate the genetic association between RFI and DMI across different environmental gradient (EG) levels. For this, we used phenotypic information on 12,958 animals (young bulls and heifers) for DMI and RFI recorded during 158 feed efficiency trials.

RESULTS

The heritability estimates for DMI and RFI across EG ranged from 0.26 to 0.54 and from 0.07 to 0.41, respectively. The average genetic correlations (± standard deviation) across EG for DMI and RFI were 0.83 ± 0.19 and 0.81 ± 0.21, respectively, with the lowest genetic correlation estimates observed between extreme EG levels (low vs. high) i.e. 0.22 for RFI and 0.26 for DMI, indicating the presence of G × E interactions. The genetic correlation between RFI and DMI across EG levels decreased as the EG became more favorable and ranged from 0.79 (lowest EG) to 0.52 (highest EG). Based on the estimated breeding values from extreme EG levels (low vs. high), we observed a moderate Spearman correlation of 0.61 (RFI) and 0.55 (DMI) and a selection coincidence of 53.3% and 40.0% for RFI and DMI, respectively.

CONCLUSIONS

Our results show evidence of G × E interactions on feed efficiency traits in Nellore cattle, especially in feeding trials with an average daily gain (ADG) that is far from the expected of 1 kg/day, thus increasing reranking of animals.

Transcriptome analysis of divergent residual feed intake phenotypes in the M. longissimus thoracis et lumborum of Wannan Yellow rabbits

Introduction: Feed efficiency is an important economic trait in rabbit meat production. The identification of molecular mechanisms and candidate genes for feed efficiency may improve the economic and environmental benefits of the rabbit meat industry. As an alternative to the conventional feed conversion ratio, residual feed intake (RFI) can be used as an accurate indicator of feed efficiency. Methods: RNA sequencing was used to identify the differentially expressed genes (DEGs) in the M. longissimus thoracis et lumborum of eight Wannan Yellow rabbits with excessively high or low RFIs (HRFI or LRFI, respectively). Thereafter, Gene Ontology (GO) analysis, enrichment using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, gene set enrichment analysis (GSEA), and protein-protein interaction (PPI) network analysis was conducted. Results: In total, 445 DEGs were identified in the M. longissimus thoracis et lumborum of rabbits with high and low RFIs. The significantly enriched GO terms identified in these two groups were primarily involved in energy and mitochondrial metabolism and oxidation-reduction processes. KEGG analysis identified 11 significantly enriched pathways, including oxidative phosphorylation, PI3K-Akt signaling, and extracellular matrix-receptor interaction pathways. According to GSEA, the expressions of genes and pathways related to mitochondrial function were upregulated in HRFI rabbits, whereas genes with upregulated expressions in LRFI rabbits were related to immune response and energy metabolism. Additionally, PPI network analysis revealed five potential candidate genetic markers. Conclusion: Comparative analysis of the M. longissimus thoracis et lumborum transcriptomes in HRFI and LRFI rabbits revealed FOS, MYC, PRKACB, ITGA2, and FN1 as potential candidate genes that affect feed efficiency in rabbits. In addition, key signaling pathways involved in oxidative phosphorylation and PI3K-Akt and ECM-receptor interaction signaling impact rabbit feed efficiency. These findings will aid in breeding programs to improve feed efficiency and optimize RFI selection of rabbits for meat production.

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.