MAK: a machine learning framework improved genomic prediction via multi-target ensemble regressor chains and automatic selection of assistant traits.

Incorporating the genotypic and phenotypic of the correlated traits into the multi-trait model can significantly improve the prediction accuracy of the target trait in animal and plant breeding, as well as human genetics. However, in most cases, the phenotypic information of the correlated and target trait of the individual to be evaluated was null simultaneously, particularly for the newborn. Therefore, we propose a machine learning framework, MAK, to improve the prediction accuracy of the target trait by constructing the multi-target ensemble regression chains and selecting the assistant trait automatically, which predicted the genomic estimated breeding values of the target trait using genotypic information only. The prediction ability of MAK was significantly more robust than the genomic best linear unbiased prediction, BayesB, BayesRR and the multi trait Bayesian method in the four real animal and plant datasets, and the computational efficiency of MAK was roughly 100 times faster than BayesB and BayesRR.

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.

Transcriptome Analysis of Compensatory Growth and Meat Quality Alteration after Varied Restricted Feeding Conditions in Beef Cattle.

Compensatory growth (CG) is a physiological response that accelerates growth following a period of nutrient limitation, with the potential to improve growth efficiency and meat quality in cattle. However, the underlying molecular mechanisms remain poorly understood. In this study, 60 Huaxi cattle were divided into one ad libitum feeding (ALF) group and two restricted feeding groups (75% restricted, RF75; 50% restricted, RF50) undergoing a short-term restriction period followed by evaluation of CG. Detailed comparisons of growth performance during the experimental period, as well as carcass and meat quality traits, were conducted, complemented by a comprehensive transcriptome analysis of the longissimus dorsi muscle using differential expression analysis, gene set enrichment analysis (GSEA), gene set variation analysis (GSVA), and weighted correlation network analysis (WGCNA). The results showed that irrespective of the restriction degree, the restricted animals exhibited CG, achieving final body weights comparable to the ALF group. Compensating animals showed differences in meat quality traits, such as pH, cooking loss, and fat content, compared to the ALF group. Transcriptomic analysis revealed 57 genes and 31 pathways differentially regulated during CG, covering immune response, acid-lipid metabolism, and protein synthesis. Notably, complement-coagulation-fibrinolytic system synergy was identified as potentially responsible for meat quality optimization in RF75. This study provides novel and valuable genetic insights into the regulatory mechanisms of CG in beef cattle.

KCRR: a nonlinear machine learning with a modified genomic similarity matrix improved the genomic prediction efficiency.

Nowadays, advances in high-throughput sequencing benefit the increasing application of genomic prediction (GP) in breeding programs. In this research, we designed a Cosine kernel-based KRR named KCRR to perform GP. This paper assessed the prediction accuracies of 12 traits with various heritability and genetic architectures from four populations using the genomic best linear unbiased prediction (GBLUP), BayesB, support vector regression (SVR), and KCRR. On the whole, KCRR performed stably for all traits of multiple species, indicating that the hypothesis of KCRR had the potential to be adapted to a wide range of genetic architectures. Moreover, we defined a modified genomic similarity matrix named Cosine similarity matrix (CS matrix). The results indicated that the accuracies between GBLUP_kinship and GBLUP_CS almost unanimously for all traits, but the computing efficiency has increased by an average of 20 times. Our research will be a significant promising strategy in future GP.

Transcriptional atlas analysis from multiple tissues reveals the expression specificity patterns in beef cattle.

BACKGROUND: A comprehensive analysis of gene expression profiling across tissues can provide necessary information for an in-depth understanding of their biological functions. We performed a large-scale gene expression analysis and generated a high-resolution atlas of the transcriptome in beef cattle. RESULTS: Our transcriptome atlas was generated from 135 bovine tissues in adult beef cattle, covering 51 tissue types of major organ systems (e.g., muscular system, digestive system, immune system, reproductive system). Approximately 94.76% of sequencing reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We detected a total of 60,488 transcripts, and 32% of them were not reported before. We identified 2654 housekeeping genes (HKGs) and 477 tissue-specific genes (TSGs) across tissues. Using weighted gene co-expression network analysis, we obtained 24 modules with 237 hub genes (HUBGs). Functional enrichment analysis showed that HKGs mainly maintain the basic biological activities of cells, while TSGs were involved in tissue differentiation and specific physiological processes. HKGs in bovine tissues were more conserved in terms of expression pattern as compared to TSGs and HUBGs among multiple species. Finally, we obtained a subset of tissue-specific differentially expressed genes (DEGs) between beef and dairy cattle and several functional pathways, which may be involved in production and health traits. CONCLUSIONS: We generated a large-scale gene expression atlas across the major tissues in beef cattle, providing valuable information for enhancing genome assembly and annotation. HKGs, TSGs, and HUBGs further contribute to better understanding the biology and evolution of multiple tissues in cattle. DEGs between beef and dairy cattle also fill in the knowledge gaps about differential transcriptome regulation of bovine tissues underlying economically important traits.

Integrating genomics and transcriptomics to identify candidate genes for subcutaneous fat deposition in beef cattle.

Fat deposition is a complex economic trait regulated by polygenic genetic basis and environmental factors. Therefore, integrating multi-omics data to uncover its internal regulatory mechanism has attracted extensive attention. Here, we performed genomics and transcriptomics analysis to detect candidates affecting subcutaneous fat (SCF) deposition in beef cattle. The association of 770K SNPs with the backfat thickness captured nine significant SNPs within or near 11 genes. Additionally, 13 overlapping genes regarding fat deposition were determined via the analysis of differentially expressed genes and weighted gene co-expression network analysis (WGCNA). We then calculated the correlations of these genes with BFT and constructed their interaction network. Finally, seven biomarkers including ACACA, SCD, FASN, ACOX1, ELOVL5, HACD2, and HSD17B12 were screened. Notably, ACACA, identified by the integration of genomics and transcriptomics, was more likely to exert profound effects on SCF deposition. These findings provided novel insights into the regulation mechanism underlying bovine fat accumulation.

Comparative transcriptomic analysis reveals region-specific expression patterns in different beef cuts.

BACKGROUND: Beef cuts in different regions of the carcass have different meat quality due to their distinct physiological function. The objective of this study was to characterize the region-specific expression differences using comparative transcriptomics analysis among five representative beef cuts (tenderloin, longissimus lumborum, rump, neck, chuck). RESULTS: We obtained 15,701 expressed genes in 30 muscle samples across five regions from carcass meat. We identified a total of 80 region-specific genes (RSGs), ranging from three (identified in the rump cut) to thirty (identified in the longissimus lumborum cut), and detected 25 transcription factors (TFs) for RSGs. Using a co-expression network analysis, we detected seven region-specific modules, including three positively correlated modules and four negatively correlated modules. We finally obtained 91 candidate genes related to meat quality, and the functional enrichment analyses showed that these genes were mainly involved in muscle fiber structure (e.g., TNNI1, TNNT1), fatty acids (e.g., SCD, LPL), amino acids (ALDH2, IVD, ACADS), ion channel binding (PHPT1, SNTA1, SUMO1, CNBP), protein processing (e.g., CDC37, GAPDH, NRBP1), as well as energy production and conversion (e.g., ATP8, COX8B, NDUFB6). Moreover, four candidate genes (ALDH2, CANX, IVD, PHPT1) were validated using RT-qPCR analyses which further supported our RNA-seq results. CONCLUSIONS: Our results provide valuable insights into understanding the transcriptome regulation of meat quality in different beef cuts, and these findings may further help to improve the selection for health-beneficial meat in beef cattle.

Towards the detection of copy number variation from single sperm sequencing in cattle.

BACKGROUND: Copy number variation (CNV) has been routinely studied using bulk-cell sequencing. However, CNV is not well studied on the single-cell level except for humans and a few model organisms. RESULTS: We sequenced 143 single sperms of two Holstein bulls, from which we predicted CNV events using 14 single sperms with deep sequencing. We then compared the CNV results derived from single sperms with the bulk-cell sequencing of one bull's family trio of diploid genomes. As a known CNV hotspot, segmental duplications were also predicted using the bovine ARS-UCD1.2 genome. Although the trio CNVs validated only some single sperm CNVs, they still showed a distal chromosomal distribution pattern and significant associations with segmental duplications and satellite repeats. CONCLUSION: Our preliminary results pointed out future research directions and highlighted the importance of uniform whole genome amplification, deep sequence coverage, and dedicated software pipelines for CNV detection using single cell sequencing data.

Genome-wide recombination map construction from single sperm sequencing in cattle.

BACKGROUND: Meiotic recombination is one of the important phenomena contributing to gamete genome diversity. However, except for human and a few model organisms, it is not well studied in livestock, including cattle. RESULTS: To investigate their distributions in the cattle sperm genome, we sequenced 143 single sperms from two Holstein bulls. We mapped meiotic recombination events at high resolution based on phased heterozygous single nucleotide polymorphism (SNP). In the absence of evolutionary selection pressure in fertilization and survival, recombination events in sperm are enriched near distal chromosomal ends, revealing that such a pattern is intrinsic to the molecular mechanism of meiosis. Furthermore, we further validated these findings in single sperms with results derived from sequencing its family trio of diploid genomes and our previous studies of recombination in cattle. CONCLUSIONS: To our knowledge, this is the first large-scale single sperm whole-genome sequencing effort in livestock, which provided useful information for future studies of recombination, genome instability, and male infertility.

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.

Genomic sequencing analysis reveals copy number variations and their associations with economically important traits in beef cattle.

Copy number variation (CNV) represents a major source of genetic variation, which may have potentially large effects, including alternating gene regulation and dosage, as well as contributing to gene expression and risk for normal phenotypic variability. We carried out a comprehensive analysis of CNV based on whole genome sequencing in Chinese Simmental beef cattle. Totally, we found 9313 deletion and 234 duplication events, covering 147.5 Mb autosomal regions. Within them, 257 deletion events of high frequency overlapped with 193 known RefGenes. Among these genes, we observed several genes were related to economically important traits, like residual feed intake, immune responding, pregnancy rate and muscle differentiation. Using a locus-based analysis, we identified 11 deletions and 1 duplication, which were significantly associated with three traits including carcass weight, tenderloin and longissimus muscle area. Our sequencing-based study provided important insights into investigating the association of CNVs with important traits in beef cattle.

Cis-eQTL Analysis and Functional Validation of Candidate Genes for Carcass Yield Traits in Beef Cattle.

Carcass yield traits are of considerable economic importance for farm animals, which act as a major contributor to the world's food supply. Genome-wide association studies (GWASs) have identified many genetic variants associated with carcass yield traits in beef cattle. However, their functions are not effectively illustrated. In this study, we performed an integrative analysis of gene-based GWAS with expression quantitative trait locus (eQTL) analysis to detect candidate genes for carcass yield traits and validate their effects on bovine skeletal muscle satellite cells (BSCs). The gene-based GWAS and cis-eQTL analysis revealed 1780 GWAS and 1538 cis-expression genes. Among them, we identified 153 shared genes that may play important roles in carcass yield traits. Notably, the identified cis-eQTLs of PON3 and PRIM2 were significantly (p < 0.001) enriched in previous GWAS loci for carcass traits. Furthermore, overexpression of PON3 and PRIM2 promoted the BSCs' proliferation, increased the expression of MYOD and downregulated the expression of MYOG, which indicated that these genes may inhibit myogenic differentiation. In contrast, PON3 and PRIM2 were significantly downregulated during the differentiation of BSCs. These findings suggested that PON3 and PRIM2 may promote the proliferation of BSCs and inhibit them in the pre-differentiation stage. Our results further contribute to the understanding of the molecular mechanisms of carcass yield traits in beef 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.

Application of ensemble learning to genomic selection in chinese simmental beef cattle.

Genomic selection (GS) using the whole-genome molecular makers to predict genomic estimated breeding values (GEBVs) is revolutionizing the livestock and plant breeding. Seeking out novel strategies with higher prediction accuracy for GS has been the ultimate goal of breeders. With the rapid development of artificial intelligence, machine learning algorithms were applied to estimate the GEBVs increasingly. Although some machine learning methods have better performance in phenotype prediction, there is still considerable room for improvement. In this study, we applied an ensemble-learning algorithm, Adaboost.RT, which integrated support vector regression (SVR), kernel ridge regression (KRR) and random forest (RF), to predict genomic breeding values of three economic traits (carcass weight, live weight, and eye muscle area) in Chinese Simmental beef cattle. Predictive accuracy measured as the Pearson correlation between the corrected phenotypes and predicted GEBVs. Moreover, we compared the reliability of SVR, KRR, RF, Adaboost.RT and GBLUP methods. The result showed that machine learning methods outperformed GBLUP, and the average improvement of four machine learning methods over the GBLUP was 12.8%, 14.9%, 5.4% and 14.4%, respectively. Among the four machine learning methods, the reliability of Adaboost.RT was comparable to KRR with higher stability. We therefore believe that the Adaboost.RT algorithm is a reliable and efficient method for GS.

Systematic analyses reveal RNA editing events involved in skeletal muscle development of goat (Capra hircus).

RNA editing is a posttranscriptional molecular process involved with specific nucleic modification, which can enhance the diversity of gene products. Adenosine-to-inosine (A-to-I, I is read as guanosine by both splicing and translation machinery) is the main type of RNA editing in mammals, which manifested as AG (adenosine-to-guanosine) in sequence data. Here, we aimed to explore patterns of RNA editing using RNA sequencing data from skeletal muscle at four developmental stages (three fetal periods and one postnatal period) in goat. We found the occurrences of RNA editing events raised at fetal periods and declined at the postnatal period. Also, we observed distinct editing levels of AG editing across stages, and significant difference was found between postnatal period and fetal periods. AG editing patterns in newborn goats are similar to those of 45-day embryo compared with embryo at 105 days and embryo at 60 days. In this study, we found a total of 1415 significantly differential edited AG sites among four groups. Moreover, 420 sites were obviously clustered into six time-series profiles, and one profile had significant association between editing level and gene expression. Our findings provided some novel insights into understanding the molecular mechanism of muscle development in mammals.

Multiple association analysis of loci and candidate genes that regulate body size at three growth stages in Simmental beef cattle.

BACKGROUND: Body size traits as one of the main breeding selection criteria was widely used to monitor cattle growth and to evaluate the selection response. In this study, body size was defined as body height (BH), body length (BL), hip height (HH), heart size (HS), abdominal size (AS), and cannon bone size (CS). We performed genome-wide association studies (GWAS) of these traits over the course of three growth stages (6, 12 and 18 months after birth) using three statistical models, single-trait GWAS, multi-trait GWAS and LONG-GWAS. The Illumina Bovine HD 770 K BeadChip was used to identify genomic single nucleotide polymorphisms (SNPs) in 1217 individuals. RESULTS: In total, 19, 29, and 10 significant SNPs were identified by the three models, respectively. Among these, 21 genes were promising candidate genes, including SOX2, SNRPD1, RASGEF1B, EFNA5, PTBP1, SNX9, SV2C, PKDCC, SYNDIG1, AKR1E2, and PRIM2 identified by single-trait analysis; SLC37A1, LAP3, PCDH7, MANEA, and LHCGR identified by multi-trait analysis; and P2RY1, MPZL1, LINGO2, CMIP, and WSCD1 identified by LONG-GWAS. CONCLUSIONS: Multiple association analysis was performed for six growth traits at each growth stage. These findings offer valuable insights for the further investigation of potential genetic mechanism of growth traits in Simmental beef cattle.

Genome-wide association studies using binned genotypes.

Linear mixed models (LMM) that tests trait association one marker at a time have been the most popular methods for genome-wide association studies. However, this approach has potential pitfalls: over conservativeness after Bonferroni correction, ignorance of linkage disequilibrium (LD) between neighboring markers, and power reduction due to overfitting SNP effects. So, multiple locus models that can simultaneously estimate and test all markers in the genome are more appropriate. Based on the multiple locus models, we proposed a bin model that combines markers into bins based on their LD relationships. A bin is treated as a new synthetic marker and we detect the associations between bins and traits. Since the number of bins can be substantially smaller than the number of markers, a penalized multiple regression method can be adopted by fitting all bins to a single model. We developed an innovative method to bin the neighboring markers and used the least absolute shrinkage and selection operator (LASSO) method. We compared BIN-Lasso with SNP-Lasso and Q + K-LMM in a simulation experiment, and showed that the new method is more powerful with less Type I error than the other two methods. We also applied the bin model to a Chinese Simmental beef cattle population for bone weight association study. The new method identified more significant associations than the classical LMM. The bin model is a new dimension reduction technique that takes advantage of biological information (i.e., LD). The new method will be a significant breakthrough in associative genomics in the big data era.

Comparative Transcriptome Analyses of Drug-sensitive and Drug-resistant Strains of Eimeria tenella by RNA-sequencing.

Avian coccidiosis is a widespread and economically significant disease in poultry. At present, treatment of coccidiosis mainly relies on drugs. Anticoccidial drugs can be divided into two categories: ionophorous compounds and synthetic drugs. However, the emergence of drug-resistant strains has become a challenge for coccidiosis control with anticoccidial drugs. To gain insights into the molecular mechanism governing the drug resistance of Eimeria tenella, two drug-resistant strains of E. tenella, one maduramicin-resistant (MRR) strain and one diclazuril-resistant (DZR) strain, were generated. We carried out comparative transcriptome analyses of a drug-sensitive strain (DS) and two drug-resistant MRR and DZR strains of E. tenella using RNA-sequencing. A total of 1,070 differentially expressed genes (DEGs), 672 upregulated and 398 downregulated, were identified in MRR vs. DS, and 379 DEGs, 330 upregulated and 49 downregulated, were detected in DZR vs. DS. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed to better understand the functions of these DEGs. In the comparison of DZR vs. DS, some DEGs were involved in peroxisome, biosynthesis of unsaturated fatty acids, and fatty acid metabolism. In the comparison of MRR vs. DS, some DEGs were involved in glycolysis/gluconeogenesis, regulation of actin cytoskeleton, and DNA replication. In addition, some DEGs coded for surface antigens that were downregulated in two drug-resistant strains involved invasion, pathogenesis, and host-parasite interactions. These results provided suggestions for further research toward unraveling the molecular mechanisms of drug resistance in Eimeria species and contribute to developing rapid molecular methods to detect resistance to these drugs in Eimeria species in poultry.

Genome-wide association study identifies the PLAG1-OXR1 region on BTA14 for carcass meat yield in cattle.

Carcass meat yield is an important carcass trait that contributes to the production efficiency and economic benefits in beef cattle. It is therefore critical to identify quantitative trait loci associated with carcass traits to enable selection. Our previous studies have identified several causal variants within the pleomorphic adenoma gene 1 (PLAG1) and coiled-coil-helix-coiled-coil-helix domain-containing 7 (CHCHD7) genes on BTA14 for carcass traits in Chinese Simmental. In the current study, we carried out a genome-wide association study for carcass meat yield in 472 Wagyu cattle with Bovine HD SNP array. Our results showed that 27 single nucleotide polymorphisms (SNPs) were identified for tenderloin weight (TDW), striploin weight (SPW), chuck roll weight (CRW), bicep weight (BPW), knuckle weight (KCW), and flank steak weight (FSW) in Wagyu cattle. Of these SNPs, 10 distinct SNPs were detected within the oxidation resistance 1 (OXR1), fatty acid binding protein 5 (FABP5), TNF receptor superfamily member 11b (TNFRSF11B), and zinc finger CCCH-type containing 3 (ZC3H3) genes on BTA14. Notably, three significant SNPs, BovineHD1400016738, BovineHD1400016743, and BovineHD1400016665 within OXR1, were shown strong linkage disequilibrium (r2 > 0.8) and significantly associated with CRW (P = 1.37 × 10-8 ~ 1.94 × 10-8). Moreover, Ingenuity Pathway Analysis showed that OXR1, FABP5, and CAP1A genes were involved in a single network and FABP5 may regulate the expression of OXR1 gene via node gene, peroxisome proliferator-activated receptor gamma (PPARG). Overall, this study suggests that OXR1 and FABP5 are candidate genes affecting carcass traits in Wagyu and the PLAG1-OXR1 region on BTA14 as a putative susceptibility locus for carcass meat yield for both Chinese Simmental and Wagyu.

Probe-based association analysis identifies several deletions associated with average daily gain in beef cattle.

BACKGROUND: Average daily gain (ADG) is an important trait that contributes to the production efficiency and economic benefits in the beef cattle industry. The molecular mechanisms of ADG have not yet been fully explored because most recent association studies for ADG are based on SNPs or haplotypes. We reported a systematic CNV discovery and association analysis for ADG in Chinese Simmental beef cattle. RESULTS: Our study identified 4912 nonredundant CNVRs with a total length of ~ 248.7 Mb, corresponding to ~ 8.9% of the cattle genome. Using probe-based CNV association, we identified 24 and 12 significant SNP probes within five deletions and two duplications for ADG, respectively. Among them, we found one common deletion with 89 kb imbedded in LHFPL Tetraspan Subfamily Member 6 (LHFPL6) at 22.9 Mb on BTA12, which has high frequency (12.9%) dispersing across population. CNV selection test using VST statistic suggested this common deletion may be under positive selection in Chinese Simmental cattle. Moreover, this deletion was not overlapped with any candidate SNP for ADG compared with previous SNPs-based association studies, suggesting its important role for ADG. In addition, we identified one rare deletion near gene Growth Factor Receptor-bound Protein 10 (GRB10) at 5.1 Mb on BTA4 for ADG using both probe-based association and region-based approaches. CONCLUSIONS: Our results provided some valuable insights to elucidate the genetic basis of ADG in beef cattle, and these findings offer an alternative perspective to understand the genetic mechanism of complex traits in terms of copy number variations in farm animals.