Genome-Wide Detection of Copy Number Variations and Their Potential Association with Carcass and Meat Quality Traits in Pingliang Red Cattle.

Copy number variation (CNV) serves as a significant source of genetic diversity in mammals and exerts substantial effects on various complex traits. Pingliang red cattle, an outstanding indigenous resource in China, possess remarkable breeding value attributed to their tender meat and superior marbling quality. However, the genetic mechanisms influencing carcass and meat quality traits in Pingliang red cattle are not well understood. We generated a comprehensive genome-wide CNV map for Pingliang red cattle using the GGP Bovine 100K SNP chip. A total of 755 copy number variable regions (CNVRs) spanning 81.03 Mb were identified, accounting for approximately 3.24% of the bovine autosomal genome. Among these, we discovered 270 potentially breed-specific CNVRs in Pingliang red cattle, including 143 gains, 73 losses, and 54 mixed events. Functional annotation analysis revealed significant associations between these specific CNVRs and important traits such as carcass and meat quality, reproduction, exterior traits, growth traits, and health traits. Additionally, our network and transcriptome analysis highlighted CACNA2D1, CYLD, UBXN2B, TG, NADK, and ITGA9 as promising candidate genes associated with carcass weight and intramuscular fat deposition. The current study presents a genome-wide CNV map in Pingliang red cattle, highlighting breed-specific CNVRs, and transcriptome findings provide valuable insights into the underlying genetic characteristics of Pingliang red cattle. These results offer potential avenues for enhancing meat quality through a targeted breeding program.

Improving Genomic Predictions in Multi-Breed Cattle Populations: A Comparative Analysis of BayesR and GBLUP Models.

Numerous studies have shown that combining populations from similar or closely related genetic breeds improves the accuracy of genomic predictions (GP). Extensive experimentation with diverse Bayesian and genomic best linear unbiased prediction (GBLUP) models have been developed to explore multi-breed genomic selection (GS) in livestock, ultimately establishing them as successful approaches for predicting genomic estimated breeding value (GEBV). This study aimed to assess the effectiveness of using BayesR and GBLUP models with linkage disequilibrium (LD)-weighted genomic relationship matrices (GRMs) for genomic prediction in three different beef cattle breeds to identify the best approach for enhancing the accuracy of multi-breed genomic selection in beef cattle. Additionally, a comparison was conducted to evaluate the predictive precision of different marker densities and genetic correlations among the three breeds of beef cattle. The GRM between Yunling cattle (YL) and other breeds demonstrated modest affinity and highlighted a notable genetic concordance of 0.87 between Chinese Wagyu (WG) and Huaxi (HX) cattle. In the within-breed GS, BayesR demonstrated an advantage over GBLUP. The prediction accuracies for HX cattle using the BayesR model were 0.52 with BovineHD BeadChip data (HD) and 0.46 with whole-genome sequencing data (WGS). In comparison to the GBLUP model, the accuracy increased by 26.8% for HD data and 9.5% for WGS data. For WG and YL, BayesR doubled the within-breed prediction accuracy to 14.3% from 7.1%, outperforming GBLUP across both HD and WGS datasets. Moreover, analyzing multiple breeds using genomic selection showed that BayesR consistently outperformed GBLUP in terms of predictive accuracy, especially when using WGS. For instance, in a mixed reference population of HX and WG, BayesR achieved a significant accuracy of 0.53 using WGS for HX, which was a substantial enhancement over the accuracies obtained with GBLUP models. The research further highlights the benefit of including various breeds in the reference group, leading to enhanced accuracy in predictions and emphasizing the importance of comprehensive genomic selection methods. Our research findings indicate that BayesR exhibits superior performance compared to GBLUP in multi-breed genomic prediction accuracy, achieving a maximum improvement of 33.3%, especially in genetically diverse breeds. The improvement can be attributed to the effective utilization of higher single nucleotide polymorphism (SNP) marker density by BayesR, resulting in enhanced prediction accuracy. This evidence conclusively demonstrates the significant impact of BayesR on enhancing genomic predictions in diverse cattle populations, underscoring the crucial role of genetic relatedness in selection methodologies. In parallel, subsequent studies should focus on refining GRM and exploring alternative models for GP.

Genetic Origin and Introgression Pattern of Pingliang Red Cattle Revealed Using Genome-Wide SNP Analyses.

The Pingliang red cattle, an outstanding indigenous resource in China, possesses an exceptional breeding value attributed to its tender meat and superior marbling quality. Currently, research efforts have predominantly concentrated on exploring its maternal origin and conducting conventional phenotypic studies. However, there remains a lack of comprehensive understanding regarding its genetic basis. To address this gap, we conducted a thorough whole-genome analysis to investigate the population structure, phylogenetic relationships, and gene flows of this breed using genomic SNP chip data from 17 bovine breeds. The results demonstrate that Pingliang red cattle have evolved distinct genetic characteristics unique to this breed, clearly distinguishing it from other breeds. Based on the analysis of the population structure and phylogenetic tree, it can be classified as a hybrid lineage between Bos taurus and Bos indicus. Furthermore, Pingliang red cattle display a more prominent B. taurus pedigree in comparison with Jinnan, Qinchuan, Zaosheng, Nanyang, and Luxi cattle. Moreover, this study also revealed closer genetic proximity within the Chinese indigenous cattle breed, particularly Qinchuan cattle, which shares the longest identical by descent (IBD) fragment with Pingliang red cattle. Gene introgression analysis shows that Pingliang red cattle have undergone gene exchange with South Devon and Red Angus cattle from Europe. Admixture analysis revealed that the proportions of East Asian taurine and Chinese indicine in the ancestry of Pingliang red cattle are approximately 52.44% and 21.00%, respectively, while Eurasian taurine, European taurine, and Indian indicine account for approximately 17.55%, 7.27%, and 1.74%. Our findings unveil distinct genetic characteristics in Pingliang red cattle and attribute their origin to B. taurus and B. indicus ancestry, as well as contributions from Qinchuan cattle, South Devon, and Red Angus.

Characterization of Chromatin Accessibility in Fetal Bovine Chondrocytes.

Despite significant advances of the bovine epigenome investigation, new evidence for the epigenetic basis of fetal cartilage development remains lacking. In this study, the chondrocytes were isolated from long bone tissues of bovine fetuses at 90 days. The Assay for Transposase-Accessible Chromatin with high throughput sequencing (ATAC-seq) and transcriptome sequencing (RNA-seq) were used to characterize gene expression and chromatin accessibility profile in bovine chondrocytes. A total of 9686 open chromatin regions in bovine fetal chondrocytes were identified and 45% of the peaks were enriched in the promoter regions. Then, all peaks were annotated to the nearest gene for Gene Ontology (GO) and Kyoto Encylopaedia of Genes and Genomes (KEGG) analysis. Growth and development-related processes such as amide biosynthesis process (GO: 0043604) and translation regulation (GO: 006417) were enriched in the GO analysis. The KEGG analysis enriched endoplasmic reticulum protein processing signal pathway, TGF-ß signaling pathway and cell cycle pathway, which are closely related to protein synthesis and processing during cell proliferation. Active transcription factors (TFs) were enriched by ATAC-seq, and were fully verified with gene expression levels obtained by RNA-seq. Among the top50 TFs from footprint analysis, known or potential cartilage development-related transcription factors FOS, FOSL2 and NFY were found. Overall, our data provide a theoretical basis for further determining the regulatory mechanism of cartilage development in bovine.

The eQTL colocalization and transcriptome-wide association study identify potentially causal genes responsible for economic traits in Simmental beef cattle.

BACKGROUND: A detailed understanding of genetic variants that affect beef merit helps maximize the efficiency of breeding for improved production merit in beef cattle. To prioritize the putative variants and genes, we ran a comprehensive genome-wide association studies (GWAS) analysis for 21 agronomic traits using imputed whole-genome variants in Simmental beef cattle. Then, we applied expression quantitative trait loci (eQTL) mapping between the genotype variants and transcriptome of three tissues (longissimus dorsi muscle, backfat, and liver) in 120 cattle. RESULTS: We identified 1,580 association signals for 21 beef agronomic traits using GWAS. We then illuminated 854,498 cis-eQTLs for 6,017 genes and 46,970 trans-eQTLs for 1,903 genes in three tissues and built a synergistic network by integrating transcriptomics with agronomic traits. These cis-eQTLs were preferentially close to the transcription start site and enriched in functional regulatory regions. We observed an average of 43.5% improvement in cis-eQTL discovery using multi-tissue eQTL mapping. Fine-mapping analysis revealed that 111, 192, and 194 variants were most likely to be causative to regulate gene expression in backfat, liver, and muscle, respectively. The transcriptome-wide association studies identified 722 genes significantly associated with 11 agronomic traits. Via the colocalization and Mendelian randomization analyses, we found that eQTLs of several genes were associated with the GWAS signals of agronomic traits in three tissues, which included genes, such as NADSYN1, NDUFS3, LTF and KIFC2 in liver, GRAMD1C, TMTC2 and ZNF613 in backfat, as well as TIGAR, NDUFS3 and L3HYPDH in muscle that could serve as the candidate genes for economic traits. CONCLUSIONS: The extensive atlas of GWAS, eQTL, fine-mapping, and transcriptome-wide association studies aid in the suggestion of potentially functional variants and genes in cattle agronomic traits and will be an invaluable source for genomics and breeding in beef cattle.

Genomic prediction of carcass traits using different haplotype block partitioning methods in beef cattle.

Genomic prediction (GP) based on haplotype alleles can capture quantitative trait loci (QTL) effects and increase predictive ability because the haplotypes are expected to be in linkage disequilibrium (LD) with QTL. In this study, we constructed haploblocks using LD-based and the fixed number of single nucleotide polymorphisms (fixed-SNP) methods with Illumina BovineHD chip in beef cattle. To evaluate the performance of different haplotype block partitioning methods, we constructed haploblocks based on LD thresholds (from r 2 > 0.2 to r 2 > 0.8) and the number of fixed-SNPs (5, 10, 20). The performance of predictive methods for three carcass traits including liveweight (LW), dressing percentage (DP), and longissimus dorsi muscle weight (LDMW) was evaluated using three approaches (GBLUP and BayesB model based on the SNP, GHBLUP, and BayesBH models based on the haploblock, and GHBLUP+GBLUP and BayesBH+BayesB models based on the combined haploblock and the nonblocked SNPs, which were located between blocks). In this study, we found the accuracies of LD-based and fixed-SNP haplotype Bayesian methods outperformed the Bayesian models (up to 8.54 ± 7.44% and 5.74 ± 2.95%, respectively). GHBLUP showed a high improvement (up to 11.29 ± 9.87%) compared with GBLUP. The Bayesian models have higher accuracies than BLUP models in most scenarios. The average computing time of the BayesBH+BayesB model can reduce by 29.3% compared with the BayesB model. The prediction accuracies using the LD-based haplotype method showed higher improvements than the fixed-SNP haplotype method. In addition, to avoid the influence of rare haplotypes generated from haplotype construction, we compared the performance of GP by filtering four types of minor haplotype allele frequency (MHAF) (0.01, 0.025, 0.05, and 0.1) under different conditions (LD levels were set at r 2 > 0.3, and the fixed number of SNPs was 5). We found the optimal MHAF threshold for LW was 0.01, and the optimal MHAF threshold for DP and LDMW was 0.025.

Transcriptional states and chromatin accessibility during bovine myoblasts proliferation and myogenic differentiation.

OBJECTIVES: Although major advances have been made in bovine epigenome study, the epigenetic basis for fetal skeletal muscle development still remains poorly understood. The aim is to recapitulated the time course of fetal skeletal muscle development in vitro, and explore the dynamic changes of chromatin accessibility and gene expression during bovine myoblasts proliferation and differentiation. METHODS: PDGFR- cells were isolated from bovine fetal skeletal muscle, then cultured and induced myogenic differentiation in vitro in a time-course study (P, D0, D2,and D4). The assay for transposase-accessible chromatin sequencing (ATAC-seq) and RNA sequencing (RNA-seq) were performed. RESULTS: Among the enriched transcriptional factors with high variability, we determined the effects of MAFF, ZNF384, and KLF6 in myogenesis using RNA interference (RNAi). In addition, we identified both stage-specific genes and chromatin accessibility regions to reveal the sequential order of gene expression, transcriptional regulatory, and signal pathways involved in bovine skeletal muscle development. Further investigation integrating chromatin accessibility and transcriptome data was conducted to explore cis-regulatory regions in line with gene expression. Moreover, we combined bovine GWAS results of growth traits with regulatory regions defined by chromatin accessibility, providing a suggestive means for a more precise annotation of genetic variants of bovine growth traits. CONCLUSION: Overall, these findings provide valuable information for understanding the stepwise regulatory mechanisms in skeletal muscle development and conducting beef cattle genetic improvement programs.

Transcriptome Analysis of Bovine Rumen Tissue in Three Developmental Stages.

Rumen development is a crucial physiological challenge for ruminants. However, the molecular mechanism regulating rumen development has not been clearly elucidated. In this study, we investigated genes involved in rumen development in 13 rumen tissues from three developmental stages (birth, youth, and adult) using RNA sequencing. We identified that 6,048 genes were differentially expressed among three developmental stages. Using weighted correlation network analysis, we found that 12 modules were significantly associated with developmental stages. Functional annotation and protein-protein interaction (PPI) network analysis revealed that CCNB1, CCNB2, IGF1, IGF2, HMGCL, BDH1, ACAT1, HMGCS2, and CREBBP involved in rumen development. Integrated transcriptome with GWAS information of carcass weight (CW), stomach weight (SW), marbling score (MS), backfat thickness (BFT), ribeye area (REA), and lean meat weight (LMW), we found that upregulated DEGs (fold change 0∼1) in birth-youth comparison were significantly enriched with GWAS signals of MS, downregulated DEGs (fold change >3) were significantly enriched with GWAS signals of SW, and fold change 0∼1 up/downregulated DEGs in birth-adult comparison were significantly enriched with GWAS signals of CW, LMW, REA, and BFT. Furthermore, we found that GWAS signals for CW, LMW, and REA were enriched in turquoise module, and GWAS signals for CW was enriched in lightgreen module. Our study provides novel insights into the molecular mechanism underlying rumen development in cattle and highlights an integrative analysis for illustrating the genetic architecture of beef complex traits.

Identification of Candidate Variants Associated With Bone Weight Using Whole Genome Sequence in Beef Cattle.

Bone weight is critical to affect body conformation and stature in cattle. In this study, we conducted a genome-wide association study for bone weight in Chinese Simmental beef cattle based on the imputed sequence variants. We identified 364 variants associated with bone weight, while 350 of them were not included in the Illumina BovineHD SNP array, and several candidate genes and GO terms were captured to be associated with bone weight. Remarkably, we identified four potential variants in a candidate region on BTA6 using Bayesian fine-mapping. Several important candidate genes were captured, including LAP3, MED28, NCAPG, LCORL, SLIT2, and IBSP, which have been previously reported to be associated with carcass traits, body measurements, and growth traits. Notably, we found that the transcription factors related to MED28 and LCORL showed high conservation across multiple species. Our findings provide some valuable information for understanding the genetic basis of body stature in beef cattle.

Assessing the Genetic Background and Selection Signatures of Huaxi Cattle Using High-Density SNP Array.

Huaxi cattle, a specialized beef cattle breed in China, has the characteristics of fast growth, high slaughter rate, and net meat rate, good reproductive performance, strong stress resistance, and wide adaptability. In this study, we evaluated the genetic diversity, population structure, and genetic relationships of Huaxi cattle and its ancestor populations at the genome-wide level, as well as detecting the selection signatures of Huaxi cattle. Principal component analysis (PCA) and phylogenetic analysis revealed that Huaxi cattle were obviously separated from other cattle populations. The admixture analysis showed that Huaxi cattle has distinct genetic structures among all populations at K = 4. It can be concluded that Huaxi cattle has formed its own unique genetic features. Using integrated haplotype score (iHS) and composite likelihood ratio (CLR) methods, we identified 143 and 199 potentially selected genes in Huaxi cattle, respectively, among which nine selected genes (KCNK1, PDLIM5, CPXM2, CAPN14, MIR2285D, MYOF, PKDCC, FOXN3, and EHD3) related to ion binding, muscle growth and differentiation, and immunity were detected by both methods. Our study sheds light on the unique genetic feature and phylogenetic relationship of Huaxi cattle, provides a basis for the genetic mechanism analysis of important economic traits, and guides further intensive breeding improvement of Huaxi cattle.

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.

Genomic Prediction Using LD-Based Haplotypes Inferred From High-Density Chip and Imputed Sequence Variants in Chinese Simmental Beef Cattle.

A haplotype is defined as a combination of alleles at adjacent loci belonging to the same chromosome that can be transmitted as a unit. In this study, we used both the Illumina BovineHD chip (HD chip) and imputed whole-genome sequence (WGS) data to explore haploblocks and assess haplotype effects, and the haploblocks were defined based on the different LD thresholds. The accuracies of genomic prediction (GP) for dressing percentage (DP), meat percentage (MP), and rib eye roll weight (RERW) based on haplotype were investigated and compared for both data sets in Chinese Simmental beef cattle. The accuracies of GP using the entire imputed WGS data were lower than those using the HD chip data in all cases. For DP and MP, the accuracy of GP using haploblock approaches outperformed the individual single nucleotide polymorphism (SNP) approach (GBLUP_In_Block) at specific LD levels. Hotelling's test confirmed that GP using LD-based haplotypes from WGS data can significantly increase the accuracies of GP for RERW, compared with the individual SNP approach (∼1.4 and 1.9% for GHBLUP and GHBLUP+GBLUP, respectively). We found that the accuracies using haploblock approach varied with different LD thresholds. The LD thresholds (r 2 ≥ 0.5) were optimal for most scenarios. Our results suggested that LD-based haploblock approach can improve accuracy of genomic prediction for carcass traits using both HD chip and imputed WGS data under the optimal LD thresholds in Chinese Simmental beef cattle.

Validation of the Prediction Accuracy for 13 Traits in Chinese Simmental Beef Cattle Using a Preselected Low-Density SNP Panel.

Chinese Simmental beef cattle play a key role in the Chinese beef industry due to their great adaptability and marketability. To achieve efficient genetic gain at a low breeding cost, it is crucial to develop a customized cost-effective low-density SNP panel for this cattle population. Thirteen growth, carcass, and meat quality traits and a BovineHD Beadchip genotyping of 1346 individuals were used to select trait-associated variants and variants contributing to great genetic variance. In addition, highly informative SNPs with high MAF in each 500 kb sliding window and in each genic region were also included separately. A low-density SNP panel consisting of 30,684 SNPs was developed, with an imputation accuracy of 97.4% when imputed to the 770 K level. Among 13 traits, the average prediction accuracy levels evaluated by genomic best linear unbiased prediction (GBLUP) and BayesA/B/Cπ were 0.22-0.47 and 0.18-0.60 for the ~30 K array and BovineHD Beadchip, respectively. Generally, the predictive performance of the ~30 K array was trait-dependent, with reduced prediction accuracies for seven traits. While differences in terms of prediction accuracy were observed among the 13 traits, the low-density SNP panel achieved moderate to high accuracies for most of the traits and even improved the accuracies for some traits.

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.

Genome-Wide Assessment of Runs of Homozygosity in Chinese Wagyu Beef Cattle.

Runs of homozygosity (ROH) are continuous homozygous regions that generally exist in the DNA sequence of diploid organisms. Identifications of ROH leading to reduction in performance can provide valuable insight into the genetic architecture of complex traits. Here, we evaluated genome-wide patterns of homozygosity and their association with important traits in Chinese Wagyu beef cattle. We identified a total of 29,271 ROH segments from 462 animals. Within each animal, an average number of ROH was 63.36 while an average length was 62.19 Mb. To evaluate the enrichment of ROH across genomes, we initially identified 280 ROH regions by merging ROH events across all individuals. Of these, nine regions containing 154 candidate genes, were significantly associated with six traits (body height, chest circumference, fat coverage, backfat thickness, ribeye area, and carcass length; p < 0.01). Moreover, we found 26 consensus ROH regions with frequencies exceeding 10%, and several regions overlapped with QTLs, which are associated with body weight, calving ease, and stillbirth. Among them, we observed 41 candidate genes, including BCKDHB, MAB21L1, SLC2A13, FGFR3, FGFRL1, CPLX1, CTNNA1, CORT, CTNNBIP1, and NMNAT1, which have been previously reported to be related to body conformation, meat quality, susceptibility, and reproductive traits. In summary, we assessed genome-wide autozygosity patterns and inbreeding levels in Chinese Wagyu beef cattle. Our study identified many candidate regions and genes overlapped with ROH for several important traits, which could be unitized to assist the design of a selection mating strategy in beef cattle.

Incorporating Genome Annotation Into Genomic Prediction for Carcass Traits in Chinese Simmental Beef Cattle.

Various methods have been proposed for genomic prediction (GP) in livestock. These methods have mainly focused on statistical considerations and did not include genome annotation information. In this study, to improve the predictive performance of carcass traits in Chinese Simmental beef cattle, we incorporated the genome annotation information into GP. Single nucleotide polymorphisms (SNPs) were annotated to five genomic classes: intergenic, gene, exon, protein coding sequences, and 3'/5' untranslated region. Haploblocks were constructed for all markers and these five genomic classes by defining a biologically functional unit, and haplotype effects were modeled in both numerical dosage and categorical coding strategies. The first-order epistatic effects among SNPs and haplotypes were modeled using a categorical epistasis model. For all makers, the extension from the SNP-based model to a haplotype-based model improved the accuracy by 5.4-9.8% for carcass weight (CW), live weight (LW), and striploin (SI). For the five genomic classes using the haplotype-based prediction model, the incorporation of gene class information into the model improved the accuracies by an average of 1.4, 2.1, and 1.3% for CW, LW, and SI, respectively, compared with their corresponding results for all markers. Including the first-order epistatic effects into the prediction models improved the accuracies in some traits and genomic classes. Therefore, for traits with moderate-to-high heritability, incorporating genome annotation information of gene class into haplotype-based prediction models could be considered as a promising tool for GP in Chinese Simmental beef cattle, and modeling epistasis in prediction can further increase the accuracy to some degree.

Genomic Prediction and Association Analysis with Models Including Dominance Effects for Important Traits in Chinese Simmental Beef Cattle.

Non-additive effects play important roles in determining genetic changes with regard to complex traits; however, such effects are usually ignored in genetic evaluation and quantitative trait locus (QTL) mapping analysis. In this study, a two-component genome-based restricted maximum likelihood (GREML) was applied to obtain the additive genetic variance and dominance variance for carcass weight (CW), dressing percentage (DP), meat percentage (MP), average daily gain (ADG), and chuck roll (CR) in 1233 Simmental beef cattle. We estimated predictive abilities using additive models (genomic best linear unbiased prediction (GBLUP) and BayesA) and dominance models (GBLUP-D and BayesAD). Moreover, genome-wide association studies (GWAS) considering both additive and dominance effects were performed using a multi-locus mixed-model (MLMM) approach. We found that the estimated dominance variances accounted for 15.8%, 16.1%, 5.1%, 4.2%, and 9.7% of the total phenotypic variance for CW, DP, MP, ADG, and CR, respectively. Compared with BayesA and GBLUP, we observed 0.5-1.1% increases in predictive abilities of BayesAD and 0.5-0.9% increases in predictive abilities of GBLUP-D, respectively. Notably, we identified a dominance association signal for carcass weight within RIMS2, a candidate gene that has been associated with carcass weight in beef cattle. Our results suggest that dominance effects yield variable degrees of contribution to the total genetic variance of the studied traits in Simmental beef cattle. BayesAD and GBLUP-D are convenient models for the improvement of genomic prediction, and the detection of QTLs using a dominance model shows promise for use in GWAS in cattle.

Theoretical Evaluation of Multi-Breed Genomic Prediction in Chinese Indigenous Cattle.

Genomic selection (GS) has been widely considered as a valuable strategy for enhancing the rate of genetic gain in farm animals. However, the construction of a large reference population is a big challenge for small populations like indigenous cattle. In order to evaluate the potential application of GS for Chinese indigenous cattle, we assessed the influence of combining multiple populations on the reliability of genomic predictions for 10 indigenous breeds of Chinese cattle using simulated data. Also, we examined the effect of different genetic architecture on prediction accuracy. In this study, we simulated a set of genotype data by a resampling approach which can reflect the realistic linkage disequilibrium pattern for multiple populations. We found within-breed evaluations yielded the highest accuracies ranged from 0.64 to 0.68 for four different simulated genetic architectures. For scenarios using multiple breeds as reference, the predictive accuracies were higher when the reference was comprised of breeds with a close relationship, while the accuracies were low when prediction were carried out among breeds. In addition, the accuracy increased in all scenarios with the heritability increased. Our results suggested that using meta-population as reference can increase accuracy of genomic predictions for small populations. Moreover, multi-breed genomic selection was feasible for Chinese indigenous populations with genetic relationships.

Identification of muscle-specific candidate genes in Simmental beef cattle using imputed next generation sequencing.

Genome-wide association studies (GWAS) have commonly been used to identify candidate genes that control economically important traits in livestock. Our objective was to detect potential candidate genes associated mainly with muscle development traits related to dimension of hindquarter in cattle. A next generation sequencing (NGS) dataset to imputed to 12 million single nucleotide polymorphisms (SNPs) (from 1252 Simmental beef cattle) were used to search for genes affecting hindquarter traits using a linear, mixed model approach. We also used haplotype and linkage disequilibrium blocks to further support our identifications. We identified 202 significant SNPs in the bovine BTA4 chromosome region associated with width of hind leg, based on a stringent statistical threshold (p = 0.05/ effective number of SNPs identified). After exploring the region around these SNPs, we found candidate genes that were potentially related to the associated markers. More importantly, we identified a region of approximately 280 Kb on the BTA4 chromosome that harbored several muscle specific candidate genes, genes to be in a potential region for muscle development. However, we also found candidate gene SLC13A1 on BTA4, which seems to be associated with bone disorders (such as chondrodysplasia) in Simmental beef cattle.

Genome-Wide Scan Identifies Selection Signatures in Chinese Wagyu Cattle Using a High-Density SNP Array.

Selective breeding can lead to genetic diversity and diverse phenotypes in farm animals. Analysis of the genomic regions under selection can provide important insights into the genetic basis of complex traits. In this study, a high-density SNP array was used for analysis of genome selection signatures in Chinese Wagyu cattle. In total, we obtained 478,903 SNPs and 24,820 no-overlap regions for |iHS| (integrated haplotype score) estimations. Under the threshold of the top 1%, 239 regions were finally identified as candidate selected regions and 162 candidate genes were found based on the UMD3.1 genome assembly. These genes were reported to be associated with fatty acids, such as Bos taurus nitric oxide synthase 1 adaptor protein (NOS1AP), Bos taurus hydroxysteroid 17-beta dehydrogenase 7 (HSD17B7), Bos taurus WD repeat domain 7 (WDR7), Bos taurus ELOVL fatty acid elongase 2 (ELOVL2), Bos taurus calpain 1 (CAPN1), Bos taurus parkin RBR E3 ubiquitin protein ligase (PRKN, also known as PARK2), Bos taurus mitogen-activated protein kinase kinase 6 (MAP2K6), meat quality, including Bos taurus ADAM metallopeptidase domain 12 (ADAM12), Bos taurus 5'-aminolevulinate synthase 1 (ALAS1), Bos taurus small integral membrane protein 13 (SMIM13) and Bos taurus potassium two pore domain channel subfamily K member 2 (KCNK2), growth, and developmental traits, such as Bos taurus insulin like growth factor 2 receptor (IGF2R), Bos taurus RAR related orphan receptor A (RORA), Bos taurus fibroblast growth factor 14 (FGF14), Bos taurus paired box 6 (PAX6) and Bos taurus LIM homeobox 6 (LHX6). In addition, we identified several genes that are associated with body size and weight, including Bos taurus sorting nexin 29 (SNX29), Bos taurus zinc finger imprinted 2 (ZIM2), Bos taurus family with sequence similarity 110 member A (FAM110A), immune system, including Bos taurus toll like receptor 9 (TLR9), Bos taurus TAFA chemokine like family member 1 (TAFA1), Bos taurus glutathione peroxidase 8 (putative) (GPX8), Bos taurus interleukin 5 (IL5), Bos taurus PR domain containing 9 (PRDM9), Bos taurus glutamate ionotropic receptor kainate type subunit 2 (GRIK2) and feed intake efficiency, Bos taurus sodium voltage-gated channel alpha subunit 9 (SCN9A), Bos taurus relaxin family peptide/INSL5 receptor 4 (RXFP4), Bos taurus RNA polymerase II associated protein 3 (RPAP3). Moreover, four GO terms of biological regulation (GO:0009987, GO:0008152) and metabolic process (GO:0003824, GO:0005488) were found based on these genes. In addition, we found that 232 candidate regions (~18 Mb) overlapped with the Quantitative trait loci (QTL)regions extracted from cattle QTLdb. Our findings imply that many genes were selected for important traits in Chinese Wagyu cattle. Moreover, these results can contribute to the understanding of the genetic basis of the studied traits during the formation of this population.