Genome-wide association study identified genomic regions and putative candidate genes affecting meat color traits in Nellore cattle

Assessment of runs of homozygosity, heterozygosity-rich regions and genomic inbreeding estimates in a subpopulation of Guzerá (Bos indicus) dual-purpose cattle

For decades, inbreeding in cattle has been evaluated using pedigree information. Nowadays, inbreeding coefficients can be obtained using genomic information such as runs of homozygosity (ROH). The aims of this study were to quantify ROH and heterozygosity-rich regions (HRR) in a subpopulation of Guzerá dual-purpose cattle, to examine ROH and HRR islands, and to compare inbreeding coefficients obtained by ROH with alternative genomic inbreeding coefficients. A subpopulation of 1733 Guzerá animals genotyped for 50k SNPs was used to obtain the ROH and HRR segments. Inbreeding coefficients by ROH (FROH ), by genomic relationship matrix based on VanRaden's method 1 using reference allele frequency in the population (FGRM ), by genomic relationship matrix based on VanRaden's method 1 using allele frequency fixed in 0.5 (FGRM_0.5 ), and by the proportion of homozygous loci (FHOM ) were calculated. A total of 15,660 ROH were identified, and the chromosome with the highest number of ROH was BTA6. A total of 4843 HRRs were identified, and the chromosome with the highest number of HRRs was BTA23. No ROH and HRR islands were identified according to established criteria, but the regions closest to the definition of an island were examined from 64 to 67 Mb of BTA6, from 36 to 37 Mb of BTA2 and from 0.50 to 1.25 Mb of BTA23. The genes identified in ROH islands have previously been associated with dairy and beef traits, while genes identified on HRR islands have previously been associated with reproductive traits and disease resistance. FROH was equal to 0.095 ± 0.084, and its Spearman correlation with FGRM was low (0.44) and moderate-high with FHOM (0.79) and with FGRM_0.5 (0.80). The inbreeding coefficients determined by ROH were higher than other cattle breeds' and higher than pedigree-based inbreeding in the Guzerá breed obtained in previous studies. It is recommended that future studies investigate the effects of inbreeding determined by ROH on the traits under selection in the subpopulation studied.

A whole genome scan reveals distinct features of selection in Zhaotong cattle of Yunnan province

Over the years, indigenous cattle have not only played an essential role in securing primary food sources but have also been utilized for labor by humans, making them invaluable genetic resources. The Zhaotong cattle, a native Chinese breed from the Yunnan province, possess excellent meat quality and resistance to heat and humidity. Here we used whole genome sequencing data of 104 animals to delve into the population structure, genomic diversity and potential positive selection signals in Zhaotong cattle. The findings of this study demonstrate that the genetic composition of Zhaotong cattle was primarily derived from Chinese indicine cattle and East Asian cattle. The nucleotide diversity of Zhaotong cattle was only lower than that of Chinese indicine cattle, which was much higher than that of other taurine cattle. Genome-wide selection scans detected a series of positive candidate regions containing multiple key genes related to bone development and metabolism (CA10, GABRG3, GLDN and NOTUM), meat quality traits (ALG8, LINGO2, MYO5B, PRKG1 and GABRB1), immune response (ADA2, BMF, LEF1 and PAK6) and heat resistance (EIF2AK4 and LEF1). In summary, this study supplies essential genetic insights into the genome diversity within Zhaotong cattle and provides a foundational framework for comprehending the genetic basis of indigenous cattle breeds.

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.

Genome-wide association study reveals genetic loci and candidate genes for meat quality traits in a four-way crossbred pig population

Meat quality traits (MQTs) have gained more attention from breeders due to their increasing economic value in the commercial pig industry. In this genome-wide association study (GWAS), 223 four-way intercross pigs were genotyped using the specific-locus amplified fragment sequencing (SLAF-seq) and phenotyped for PH at 45 min post mortem (PH45), meat color score (MC), marbling score (MA), water loss rate (WL), drip loss (DL) in the longissimus muscle, and cooking loss (CL) in the psoas major muscle. A total of 227, 921 filtered single nucleotide polymorphisms (SNPs) evenly distributed across the entire genome were detected to perform GWAS. A total of 64 SNPs were identified for six meat quality traits using the mixed linear model (MLM), of which 24 SNPs were located in previously reported QTL regions. The phenotypic variation explained (PVE) by the significant SNPs was from 2.43% to 16.32%. The genomic heritability estimates based on SNP for six meat-quality traits were low to moderate (0.07-0.47) being the lowest for CL and the highest for DL. A total of 30 genes located within 10 kb upstream or downstream of these significant SNPs were found. Furthermore, several candidate genes for MQTs were detected, including pH45 (GRM8), MC (ANKRD6), MA (MACROD2 and ABCG1), WL (TMEM50A), CL (PIP4K2A) and DL (CDYL2, CHL1, ABCA4, ZAG and SLC1A2). This study provided substantial new evidence for several candidate genes to participate in different pork quality traits. The identification of these SNPs and candidate genes provided a basis for molecular marker-assisted breeding and improvement of pork quality traits.

Genomic Analysis of Heterosis in an Angus × Hereford Cattle Population

Crossbreeding is widely used in the beef cattle industry to exploit its several benefits. This study evaluated the effects of heterozygosity on growth traits in an Angus × Hereford cross-population. Moreover, a genome wide association study was conducted to detect regions in the genome with significant dominance effects on growth traits contributing to heterosis. A total of 1530 animals comprising of pure Line 1 Hereford, Angus and Angus × Line 1 Hereford cross. Growth phenotypes included birth weight, weaning weight and yearling weight. All animals were genotyped with GeneSeek GGP LD 50k. Significant effects of genomic heterozygosity on growth traits were detected. These effects were 0.03 kg (0.006), 5.13 kg (0.04), 6.02 kg (0.08) on birth weight, weaning weight and yearling weight, respectively. Genome wide association study revealed several SNP markers with significant heterotic effects associated with birth weight, weaning weight and yearling weight. These SNP markers were located on chromosomes 1, 2, 6, 21, 14, 19, 13 and 12. Genes in these regions were reported to be involved in growth and other important physiological mechanisms. Our study revealed several regions associated with dominance effects and contributing to heterosis. These results could be beneficial in optimizing crossbreeding.

Runs of Homozygosity Analysis Reveals Genomic Diversity and Population Structure of an Indigenous Cattle Breed in Southwest China

In aiming to achieve sustainable development goals in the livestock industry, it is becoming increasingly necessary and important for the effective conservation of genetic resources. There are some indigenous cattle breeds in Sichuan, southwest China, for which, however, the genetic diversity and population structures still remain unknown because of the unavailability of systematic breeding programs and pedigree information. Xieka cattle are an indigenous breed locally distributed in southeastern Sichuan and have a long-term evolutionary adaptation to local environments and climates. In this study, we obtained 796,828 single nucleotide polymorphisms (SNPs) through sequencing the genomes of 30 Xieka cattle and used them for analyzing the genetic diversity and runs of homozygosity (ROH). The mean nucleotide diversity was 0.28 and 72% of SNPs were found to be in the heterozygous states. A total of 4377 ROH were detected with even distribution among all autosomes, and 74% of them were lower than 1 Mb in length. Meanwhile, only five ROH were found longer than 5 Mb. We further determined 19 significant genomic regions that were obviously enriched by ROH, in which 35 positional candidate genes were found. Some of these genes have been previously reported to be significantly associated with various production traits in cattle, such as meat quality, carcass performances, and diseases. In conclusion, the relatively high degree of genetic diversity of Xieka cattle was revealed using the genomic information, and the proposed candidate genes will help us optimize the breeding programs regarding this indigenous breed.

Peroxisome Proliferator-Activated Receptor Activation in Precision-Cut Bovine Liver Slices Reveals Novel Putative PPAR Targets in Periparturient Dairy Cows

Metabolic challenges experienced by dairy cows during the transition between pregnancy and lactation (also known as peripartum), are of considerable interest from a nutrigenomic perspective. The mobilization of large amounts of non-esterified fatty acids (NEFA) leads to an increase in NEFA uptake in the liver, the excess of which can cause hepatic accumulation of lipids and ultimately fatty liver. Interestingly, peripartum NEFA activate the Peroxisome Proliferator-activated Receptor (PPAR), a transcriptional regulator with known nutrigenomic properties. The study of PPAR activation in the liver of periparturient dairy cows is thus crucial; however, current in vitro models of the bovine liver are inadequate, and the isolation of primary hepatocytes is time consuming, resource intensive, and prone to errors, with the resulting cells losing characteristic phenotypical traits within hours. The objective of the current study was to evaluate the use of precision-cut liver slices (PCLS) from liver biopsies as a model for PPAR activation in periparturient dairy cows. Three primiparous Jersey cows were enrolled in the experiment, and PCLS from each were prepared prepartum (−8.0 ± 3.6 DIM) and postpartum (+7.7± 1.2 DIM) and treated independently with a variety of PPAR agonists and antagonists: the PPARα agonist WY-14643 and antagonist GW-6471; the PPARδ agonist GW-50156 and antagonist GSK-3787; and the PPARγ agonist rosiglitazone and antagonist GW-9662. Gene expression was assayed through RT-qPCR and RNAseq, and intracellular triacylglycerol (TAG) concentration was measured. PCLS obtained from postpartum cows and treated with a PPARγ agonist displayed upregulation of ACADVL and LIPC while those treated with PPARδ agonist had increased expression of LIPC, PPARD, and PDK4. In PCLS from prepartum cows, transcription of LIPC was increased by all PPAR agonists and NEFA. TAG concentration tended to be larger in tissue slices treated with PPARδ agonist compared to CTR. Use of PPAR isotype-specific antagonists in PCLS cultivated in autologous blood serum failed to decrease expression of PPAR targets, except for PDK4, which was confirmed to be a PPARδ target. Transcriptome sequencing revealed considerable differences in response to PPAR agonists at a false discovery rate-adjusted p-value of 0.2, with the most notable effects exerted by the PPARδ and PPARγ agonists. Differentially expressed genes were mainly related to pathways involved with lipid metabolism and the immune response. Among differentially expressed genes, a subset of 91 genes were identified as novel putative PPAR targets in the bovine liver, by cross-referencing our results with a publicly available dataset of predicted PPAR target genes, and supplementing our findings with prior literature. Our results provide important insights on the use of PCLS as a model for assaying PPAR activation in the periparturient dairy cow.

Within- and between-Breed Selection Signatures in the Original and Improved Valachian Sheep

This study explored the genomic diversity and selection signatures in two Slovakian national breeds, the Original Valachian and the Improved Valachian sheep. As they are an important animal genetic resource within the country, but with decreasing population size, our aim is to identify potentially valuable genomic regions. A total of 97 sheep (18 male and 79 female) from the Original Valachian, and 69 sheep (25 male and 44 female) from the Improved Valachian populations were genotyped using the GeneSeek GGP Ovine 50 K chip. The inbreeding levels were assessed with runs of homozygosity (ROH). The selection signatures within breeds were identified based on the top 1% of most homozygous regions within the breed, the so-called ROH islands. The selection signatures between breeds were assessed based on variance in linkage disequilibrium. Overall, we have identified selection signatures with quantitative trait loci (QTL) and genes pointing towards all three production purposes of the Valachian sheep, milk, meat, and wool, including their quality characteristics. Another group with apparent large importance was the various traits related to health and resistance to parasites, which is well in line with the sturdy nature of this breed.

An evaluation of maternal serum dynamic thiol-disulfide homeostasis and ischemia modified albumin changes in pregnant women with COVID-19

Objective:

It is thought that oxidative stress, free radicals, reactive oxygen species and reactive nitrogen species affect the pathophysiology of coronavirus disease-2019 (COVID-19). This study aimed to evaluate the oxidative status in pregnant patients with COVID-19 infection according to the changes seen in the levels of maternal serum thiol-disulfide and ischemia-modified albumin (IMA).

Materials and Methods:

A study group was formed of 40 pregnant women with confirmed COVID-19 infection (study group) and a control group of 40 healthy pregnant women with no risk factors determined. In this prospective, case-controlled study, analyses were made of the maternal serum native thiol, total thiol, disulfide, IMA, and disulfide/native thiol concentrations.

Results:

The maternal serum native thiol and total thiol concentrations in the study group were determined to be statistically significantly lower (p=0.007 and p=0.006, respectively), and the disulfide/native thiol ratio was higher but not to a level of statistical significance (p=0.473). There was no difference between the two groups regarding IMA levels (p=0.731).

Conclusion:

The thiol-disulfide balance was seen to shift in the oxidant direction in pregnancies with COVID-19, which might support the view that ischemic processes play a role in the etiopathogenesis of this novel disease.

Small genetic variation affecting mRNA isoforms associated with marbling and meat color in beef cattle

The aim of this study was to identify mRNA isoforms and small genetic variants that may be affecting marbling and beef color in Nellore cattle. Longissimus thoracis muscle samples from 20 bulls with different phenotypes (out of 80 bulls set) for marbling (moderate (n = 10) and low (n = 10) groups) and beef color (desirable (n = 10) and undesirable (n = 9) group) traits were used to perform transcriptomic analysis using RNA sequencing. Fourteen and 15 mRNA isoforms were detected as differentially expressed (DE) (P-value ≤ 0.001) between divergent groups for marbling and meat color traits, respectively. Some of those DE mRNA isoforms have shown sites of splicing modified by small structural variants as single nucleotide variant (SNV), insertion, and/or deletion. Enrichment analysis identified metabolic pathways, such as O2/CO2 exchange in erythrocytes, tyrosine biosynthesis, and phenylalanine degradation. The results obtained suggest potential key regulatory genes associated with these economically important traits for the beef industry and for the consumer.

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.

Next Generation Sequencing of Single Nucleotide Polymorphic DNA-Markers in Selecting for Intramuscular Fat, Fat Melting Point, Omega-3 Long-Chain Polyunsaturated Fatty Acids and Meat Eating Quality in Tattykeel Australian White MARGRA Lamb

Meat quality data can only be obtained after slaughter when selection decisions about the live animal are already too late. Carcass estimated breeding values present major precision problems due to low accuracy, and by the time an informed decision on the genetic merit for meat quality is made, the animal is already dead. We report for the first time, a targeted next-generation sequencing (NGS) of single nucleotide polymorphisms (SNP) of lipid metabolism genes in Tattykeel Australian White (TAW) sheep of the MARGRA lamb brand, utilizing an innovative and minimally invasive muscle biopsy sampling technique for directly quantifying the genetic worth of live lambs for health-beneficial omega-3 long-chain polyunsaturated fatty acids (n-3 LC-PUFA), intramuscular fat (IMF), and fat melting point (FMP). NGS of stearoyl-CoA desaturase (SCD), fatty acid binding protein-4 (FABP4), and fatty acid synthase (FASN) genes identified functional SNP with unique DNA marker signatures for TAW genetics. The SCD g.23881050T>C locus was significantly associated with IMF, C22:6n-3, and C22:5n-3; FASN g.12323864A>G locus with FMP, C18:3n-3, C18:1n-9, C18:0, C16:0, MUFA, and FABP4 g.62829478A>T locus with IMF. These add new knowledge, precision, and reliability in directly making early and informed decisions on live sheep selection and breeding for health-beneficial n-3 LC-PUFA, FMP, IMF and superior meat-eating quality at the farmgate level. The findings provide evidence that significant associations exist between SNP of lipid metabolism genes and n-3 LC-PUFA, IMF, and FMP, thus underpinning potential marker-assisted selection for meat-eating quality traits in TAW lambs.

Runs of homozygosity analysis reveals consensus homozygous regions affecting production traits in Chinese Simmental beef cattle

BACKGROUND

Genomic regions with a high frequency of runs of homozygosity (ROH) are related to important traits in farm animals. We carried out a comprehensive analysis of ROH and evaluated their association with production traits using the BovineHD (770 K) SNP array in Chinese Simmental beef cattle.

RESULTS

We detected a total of 116,953 homozygous segments with 2.47Gb across the genome in the studied population. The average number of ROH per individual was 99.03 and the average length was 117.29 Mb. Notably, we detected 42 regions with a frequency of more than 0.2. We obtained 17 candidate genes related to body size, meat quality, and reproductive traits. Furthermore, using Fisher's exact test, we found 101 regions were associated with production traits by comparing high groups with low groups in terms of production traits. Of those, we identified several significant regions for production traits (P < 0.05) by association analysis, within which candidate genes including ECT2, GABRA4, and GABRB1 have been previously reported for those traits in beef cattle.

CONCLUSIONS

Our study explored ROH patterns and their potential associations with production traits in beef cattle. These results may help to better understand the association between production traits and genome homozygosity and offer valuable insights into managing inbreeding by designing reasonable breeding programs in farm animals.

Integration of selection signatures and multi-trait GWAS reveals polygenic genetic architecture of carcass traits in beef cattle

Carcass merits are widely considered as economically important traits affecting beef production in the beef cattle industry. However, the genetic basis of carcass traits remains to be well understood. Here, we applied multiple methods, including the Composite of Likelihood Ratio (CLR) and Genome-wide Association Study (GWAS), to explore the selection signatures and candidate variants affecting carcass traits. We identified 11,600 selected regions overlapping with 2214 candidate genes, and most of those were enriched in binding and gene regulation. Notably, we identified 66 and 110 potential variants significantly associated with carcass traits using single-trait and multi-traits analyses, respectively. By integrating selection signatures with single and multi-traits associations, we identified 12 and 27 putative genes, respectively. Several highly conserved missense variants were identified in OR5M13D, NCAPG, and TEX2. Our study supported polygenic genetic architecture of carcass traits and provided novel insights into the genetic basis of complex traits in beef cattle.