Whole-genome QTL scan for ultrasound and carcass merit traits in beef cattle using Bayesian shrinkage method

Genome-wide association study for meat tenderness in beef cattle identifies patterns of the genetic contribution in different post-mortem stages

The beef tenderization process during the post-mortem period is one of the most important sensorial attributes and it is well-established. The aim of this study was to identify the genetic contribution pattern to meat tenderness at 7-(LMD7), 14-(LMD14), and 21-(LMD21) days post-mortem. The heritabilities for LMD7 (0.194), LMD14 (0.142) and LMD21 (0.048) are well established in the population evaluated here. However, its genetic contribution in terms of genomic candidate regions is still poorly understood. Tenderness was measured in the Longissiums thoracis using Warner-Bratzler shear force in the three post-mortem periods. A total of 4323 crossbred beef cattle were phenotyped and genotyped using the Illumina BovineSNP50K. The percentage of the total genetic variance was estimated using the weighted single-step genomic best linear unbiased prediction method. The main candidate windows for LMD7 were associated with proteolysis of myofibrillar structures and the weakening endomysium and perimysium. Candidate windows for LMD14 and LMD21 were mapped in bovine QTLs for body composition, height and growth. Results presented herein highlight, the largest contribution of proteolysis related processes before 14-days post-mortem and body composition characteristics in later stages for meat tenderness.

Genomic Approaches Reveal Pleiotropic Effects in Crossbred Beef Cattle

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

Identification of Genomic Regions for Carcass Quality Traits within the American Simmental Association Carcass Merit Program

USDA quality and yield grade are primary driving forces for carcass value in the United States. Carcass improvements can be achieved by making selection decisions based on the results of genetic evaluations in the form of expected progeny differences (EPD), real-time ultrasound imaging, and physical evaluation of candidate breeding animals. In an effort to advance their ability to accurately predict the breeding value of potential sires for carcass traits, the American Simmental Association launched the Carcass Merit Program as a means to collect progeny sire group carcass information. All records were extracted from the American Simmental Association database. Progeny data were organized by sire family and progeny performance phenotypes were constructed. Sire genotypes were filtered, and a multi-locus mixed linear model was used to perform an association analysis on the genotype data, while correcting for cryptic relatedness and pedigree structure. Three chromosomes were found to have genome-wide significance and this conservative approach identified putative QTL in those regions. Three hundred ninety-three novel regions were identified across all traits, as well as 290 novel positional candidate genes. Correlations between carcass characteristics and maternal traits were less unfavorable than those previously reported.

Transcriptome profiling of longissimus lumborum in Holstein bulls and steers with different beef qualities

Previous research regarding Holstein cows has mainly focused on increasing milk yield. However, in order to maximize the economical profits of Holstein cattle farming, it is necessary to fully take advantage of Holstein bulls to produce high-grade beef. The present study aims to investigate different transcriptomic profiling of Holstein bulls and steers, via high-throughput RNA-sequencing (RNA-seq). The growth and beef quality traits of Holstein steers and bulls were characterized via assessment of weight, rib eye area, marbling score, shear force and intramuscular fat percentage of the longissimus lumborum (LL) muscle. The results indicated that castration improved the meat quality, yet reduced the meat yield. Subsequently, RNA-seq of the LL muscle from Holstein steers and bulls revealed a total of 56 differentially expressed genes (DEGs). We performed the functional enrichment analysis in Gene Ontology (GO) annotations of the DEGs using GOseq R package software and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using KOBAS tool. Through the integrated analysis of DEGs with reported QTLs and SNPs, seven promising candidate genes potentially affecting the beef quality of LL muscle following castration were discovered, including muscle structural protein coding genes (MYH1, MYH4, MYH10) and functional protein coding genes (GADL1, CYP2R1, EEPD1, SHISA3). Among them, MYH10, GADL1, CYP2R1, EEPD1 and SHISA3 were novel candidate genes associated with beef quality traits. Notably, EEPD1 was associated with both meat quality and reproduction traits, thus indicating its overlapping role in responding to hormone change, and subsequently inducing beef quality improvement. Our findings provide a complete dataset of gene expression profile of LL in Holstein bulls and steers, and will aid in understanding how castration influence meat yield and quality.

Co-expression network analysis predicts a key role of microRNAs in the adaptation of the porcine skeletal muscle to nutrient supply

BACKGROUND

The role of non-coding RNAs in the porcine muscle metabolism is poorly understood, with few studies investigating their expression patterns in response to nutrient supply. Therefore, we aimed to investigate the changes in microRNAs (miRNAs), long intergenic non-coding RNAs (lincRNAs) and mRNAs muscle expression before and after food intake.

RESULTS

We measured the miRNA, lincRNA and mRNA expression levels in the gluteus medius muscle of 12 gilts in a fasting condition (AL-T0) and 24 gilts fed ad libitum during either 5 h. (AL-T1, N = 12) or 7 h. (AL-T2, N = 12) prior to slaughter. The small RNA fraction was extracted from muscle samples retrieved from the 36 gilts and sequenced, whereas lincRNA and mRNA expression data were already available. In terms of mean and variance, the expression profiles of miRNAs and lincRNAs in the porcine muscle were quite different than those of mRNAs. Food intake induced the differential expression of 149 (AL-T0/AL-T1) and 435 (AL-T0/AL-T2) mRNAs, 6 (AL-T0/AL-T1) and 28 (AL-T0/AL-T2) miRNAs and none lincRNAs, while the number of differentially dispersed genes was much lower. Among the set of differentially expressed miRNAs, we identified ssc-miR-148a-3p, ssc-miR-22-3p and ssc-miR-1, which play key roles in the regulation of glucose and lipid metabolism. Besides, co-expression network analyses revealed several miRNAs that putatively interact with mRNAs playing key metabolic roles and that also showed differential expression before and after feeding. One case example was represented by seven miRNAs (ssc-miR-148a-3p, ssc-miR-151-3p, ssc-miR-30a-3p, ssc-miR-30e-3p, ssc-miR-421-5p, ssc-miR-493-5p and ssc-miR-503) which putatively interact with the PDK4 mRNA, one of the master regulators of glucose utilization and fatty acid oxidation.

CONCLUSIONS

As a whole, our results evidence that microRNAs are likely to play an important role in the porcine skeletal muscle metabolic adaptation to nutrient availability.

GWAS and fine-mapping of livability and six disease traits in Holstein cattle

BACKGROUND

Health traits are of significant economic importance to the dairy industry due to their effects on milk production and associated treatment costs. Genome-wide association studies (GWAS) provide a means to identify associated genomic variants and thus reveal insights into the genetic architecture of complex traits and diseases. The objective of this study is to investigate the genetic basis of seven health traits in dairy cattle and to identify potential candidate genes associated with cattle health using GWAS, fine mapping, and analyses of multi-tissue transcriptome data.

RESULTS

We studied cow livability and six direct disease traits, mastitis, ketosis, hypocalcemia, displaced abomasum, metritis, and retained placenta, using de-regressed breeding values and more than three million imputed DNA sequence variants. After data edits and filtering on reliability, the number of bulls included in the analyses ranged from 11,880 (hypocalcemia) to 24,699 (livability). GWAS was performed using a mixed-model association test, and a Bayesian fine-mapping procedure was conducted to calculate a posterior probability of causality to each variant and gene in the candidate regions. The GWAS detected a total of eight genome-wide significant associations for three traits, cow livability, ketosis, and hypocalcemia, including the bovine Major Histocompatibility Complex (MHC) region associated with livability. Our fine-mapping of associated regions reported 20 candidate genes with the highest posterior probabilities of causality for cattle health. Combined with transcriptome data across multiple tissues in cattle, we further exploited these candidate genes to identify specific expression patterns in disease-related tissues and relevant biological explanations such as the expression of Group-specific Component (GC) in the liver and association with mastitis as well as the Coiled-Coil Domain Containing 88C (CCDC88C) expression in CD8 cells and association with cow livability.

CONCLUSIONS

Collectively, our analyses report six significant associations and 20 candidate genes of cattle health. With the integration of multi-tissue transcriptome data, our results provide useful information for future functional studies and better understanding of the biological relationship between genetics and disease susceptibility in cattle.

Identifying Loci Under Positive Selection in Yellow Korean Cattle (Hanwoo)

Jeju Black cattle is one of the aboriginal Korean cattle breeds that has been isolated in Jeju island for a long time, while Yellow Hanwoo cattle has been extensively selected for beef production traits for the last several decades. Aiming to investigate broader patterns of selection, we genotyped 352 Yellow Hanwoo and 169 Jeju Black cattle using a customized 150 K bovine chip. Our composite selection signals’ analysis to identify selection signatures (cross-population extended haplotype homozygosity [XP-EHH], ΔSAF, and F_ST) identified recent and strong signature of selection near many loci with mutations affecting the traits under strong selection as outlier in Yellow Hanwoo, including SCP2 (P = 8.41 × 10⁻¹⁰) that may be involved in the meat quality. We found nine candidate regions with significant clusters of selection signals, and further bioinformatics analyses of the genes located within these regions revealed mainly genes involved in G-protein coupled receptor signaling pathway (GO:0007186) or olfactory transduction (bta04740), which may be due to adaptation to natural environments in Jeju island. Based on the stronger correlation of Ne₁₀/Ne₁₀₀ ratio between Yellow Hanwoo (0.61) and Jeju Black (0.66) cattle, our results suggest that the difference of chromosomal regions of selection signature between the 2 cattle breeds was due to a consequence of selection processes to adapt to environmental differences between Jeju island and the main inland, Korean peninsula.

Genomewide association study for production and meat quality traits in Canchim beef cattle

The commercial value of the bovine carcass is determined by a set of traits, such as weight, yield, back fat thickness, and marbling; therefore, the genetic improvement of growth, meat, and carcass quality traits is an important tool to add value to the supply chain. Genomewide association studies (GWAS) enable the identification of loci that control phenotypic expression of quantitative traits (QTL). Therefore, the objective of this work was to perform a GWAS to identify genomic regions and genes associated with growth, carcass traits, and meat quality in Canchim beef cattle. These traits were yearling weight (YW), rib eye area (REA), back fat thickness (BFT), and marbling (MARB). To increase sample size and marker density, genotype imputation was performed, and only markers imputed with greater than 95% accuracy were used. Genomewide association study was performed using a Bayesian approach, by the Bayes B statistical method, incorporating genotypes and phenotypes from 614 animals from both the Canchim breed and the MA genetic group (offspring of Charolais bulls and one-half Canchim + one-half Zebu cows). This investigation identified 1 and 4 genomic regions explaining 0.23 and 7.35% of the genetic variance for REA and YW, respectively. These regions harbor a total of 19 genes, 7 of which were classified for biological functions by functional analysis. Significant associations were not observed for BFT and MARB. The identification of QTL that had been previously described in the literature reinforces associations found in this study.

TRIENNIAL GROWTH AND DEVELOPMENT SYMPOSIUM: Factors influencing bovine intramuscular adipose tissue development and cellularity

Appearance, distribution, and amount of intramuscular fat (IMF), often referred to as marbling, are highly variable and depend on environmental and genetic factors. On the molecular level, the concerted action of several drivers, including hormones, receptors, transcription factors, etc., determines where clusters of adipocytes arise. Therefore, the aim of future studies remains to identify such factors as biological markers of IMF to increase the ability to identify animals that deposit IMF early in age to increase efficiency of high-quality meat production. In an attempt to unravel the cellular development of marbling, we investigated the abundance of markers for adipogenic differentiation during fattening of cattle and the transcriptome of muscle and dissected IMF. Markers of different stages of adipogenic differentiation are well known from cell culture experiments. They are usually transiently expressed, such as delta-like homolog 1 (DLK1) that is abundant in preadipocytes and absent during differentiation to mature adipocytes. It is even a greater challenge to detect those markers in live animals. Within skeletal muscles, hyperplasia and hypertrophy of adipocytes can be observed throughout life. Therefore, development of marbling requires, on the cellular level, recruitment, proliferation, and differentiation of adipogenic cells to store excess energy in the form of lipids in new cells. In a recent study, we investigated the localization and abundance of early markers of adipogenic differentiation, such as DLK1, in bovine muscle tissue. An inverse relationship between IMF content and number of DLK1-positive cells in bovine muscle was demonstrated. Considering the cellular environment of differentiating adipocytes in muscle and the secretory action of adipocytes and myocytes, it becomes obvious that cross talk between cells via adipokines and myokines may be important for IMF development. Secreted proteins can act on other cells, inhibiting or stimulating their function via autocrine and paracrine actions. Such factors with potential influence on IMF, among them, agouti signaling protein and thrombospondin 4, were identified in transcriptome analyses and further investigated. Furthermore, results from transcriptome analysis indicate involvement of genes that are not directly related to adipogenesis and lipid metabolism, providing new candidates for future research.

Genome-wide association study identifies loci and candidate genes for meat quality traits in Simmental beef cattle

Improving meat quality is the best way to enhance profitability and strengthen competitiveness in beef industry. Identification of genetic variants that control beef quality traits can help breeders design optimal breeding programs to achieve this goal. We carried out a genome-wide association study for meat quality traits in 1141 Simmental cattle using the Illumina Bovine HD 770K SNP array to identify the candidate genes and genomic regions associated with meat quality traits for beef cattle, including fat color, meat color, marbling score, longissimus muscle area, and shear force. In our study, we identified twenty significant single-nucleotide polymorphisms (SNPs) (p < 1.47 × 10(-6)) associated with these five meat quality traits. Notably, we observed several SNPs were in or near eleven genes which have been reported previously, including TMEM236, SORL1, TRDN, S100A10, AP2S1, KCTD16, LOC506594, DHX15, LAMA4, PREX1, and BRINP3. We identified a haplotype block on BTA13 containing five significant SNPs associated with fat color trait. We also found one of 19 SNPs was associated with multiple traits (shear force and longissimus muscle area) on BTA7. Our results offer valuable insights to further explore the potential mechanism of meat quality traits in Simmental beef cattle.

Whole-Genome Resequencing Analysis of Hanwoo and Yanbian Cattle to Identify Genome-Wide SNPs and Signatures of Selection

Over the last 30 years, Hanwoo has been selectively bred to improve economically important traits. Hanwoo is currently the representative Korean native beef cattle breed, and it is believed that it shared an ancestor with a Chinese breed, Yanbian cattle, until the last century. However, these two breeds have experienced different selection pressures during recent decades. Here, we whole-genome sequenced 10 animals each of Hanwoo and Yanbian cattle (20 total) using the Illumina HiSeq 2000 sequencer. A total of approximately 3.12 and 3.07 billion sequence reads were mapped to the bovine reference sequence assembly (UMD 3.1) at an average of approximately 10.71- and 10.53-fold coverage for Hanwoo and Yanbian cattle, respectively. A total of 17,936,399 single nucleotide polymorphisms (SNPs) were yielded, of which 22.3% were found to be novel. By annotating the SNPs, we further retrieved numerous nonsynonymous SNPs that may be associated with traits of interest in cattle. Furthermore, we performed whole-genome screening to detect signatures of selection throughout the genome. We located several promising selective sweeps that are potentially responsible for economically important traits in cattle; the PPP1R12A gene is an example of a gene that potentially affects intramuscular fat content. These discoveries provide valuable genomic information regarding potential genomic markers that could predict traits of interest for breeding programs of these cattle breeds.

Genome-wide association studies using haplotypes and individual SNPs in Simmental cattle

Recent advances in high-throughput genotyping technologies have provided the opportunity to map genes using associations between complex traits and markers. Genome-wide association studies (GWAS) based on either a single marker or haplotype have identified genetic variants and underlying genetic mechanisms of quantitative traits. Prompted by the achievements of studies examining economic traits in cattle and to verify the consistency of these two methods using real data, the current study was conducted to construct the haplotype structure in the bovine genome and to detect relevant genes genuinely affecting a carcass trait and a meat quality trait. Using the Illumina BovineHD BeadChip, 942 young bulls with genotyping data were introduced as a reference population to identify the genes in the beef cattle genome significantly associated with foreshank weight and triglyceride levels. In total, 92,553 haplotype blocks were detected in the genome. The regions of high linkage disequilibrium extended up to approximately 200 kb, and the size of haplotype blocks ranged from 22 bp to 199,266 bp. Additionally, the individual SNP analysis and the haplotype-based analysis detected similar regions and common SNPs for these two representative traits. A total of 12 and 7 SNPs in the bovine genome were significantly associated with foreshank weight and triglyceride levels, respectively. By comparison, 4 and 5 haplotype blocks containing the majority of significant SNPs were strongly associated with foreshank weight and triglyceride levels, respectively. In addition, 36 SNPs with high linkage disequilibrium were detected in the GNAQ gene, a potential hotspot that may play a crucial role for regulating carcass trait components.

Genome-wide association study for feedlot average daily gain in Nellore cattle (Bos indicus)

The genome-wide association study (GWAS) results are presented for average daily gain (ADG) in Nellore cattle. Phenotype of 720 male Bos indicus animals with information of ADG in feedlots and 354,147 single-nucleotide polymorphisms (SNPs) obtained from a database added by information from Illumina Bovine HD (777,962 SNPs) and Illumina BovineSNP50 (54,609) by imputation were used. After quality control and imputation, 290,620 SNPs remained in the association analysis, using R package Genome-wide Rapid Association using Mixed Model and Regression method GRAMMAR-Gamma. A genomic region with six significant SNPs, at Bonferroni-corrected significance, was found on chromosome 3. The most significant SNP (rs42518459, BTA3: 85849977, p = 9.49 × 10(-8)) explained 5.62% of the phenotypic variance and had the allele substitution effect of -0.269 kg/day. Important genes such as PDE4B, LEPR, CYP2J2 and FGGY are located near this region, which is overlapped by 12 quantitative trait locus (QTLs) described for several production traits. Other regions with markers with suggestive effects were identified in BTA6 and BTA10. This study showed regions with major effects on ADG in Bos indicus in feedlots. This information may be useful to increase the efficiency of selecting this trait and to understand the physiological processes involved in its regulation.

MicroRNAs in bovine adipogenesis: genomic context, expression and function

Background

MicroRNAs (miRNAs) are small non-coding RNAs found to regulate several biological processes including adipogenesis. Understanding adipose tissue regulation is critical for beef cattle as fat is an important determinant of beef quality and nutrient value. This study analyzed the association between genomic context characteristics of miRNAs with their expression and function in bovine adipose tissue. Twenty-four subcutaneous adipose tissue biopsies were obtained from eight British-continental crossbred steers at 3 different time points. Total RNA was extracted and miRNAs were profiled using a miRNA microarray with expression further validated by qRT-PCR.

Results

A total of 224 miRNAs were detected of which 155 were expressed in all steers (n = 8), and defined as the core miRNAs of bovine subcutaneous adipose tissue. Core adipose miRNAs varied in terms of genomic location (59.5% intergenic, 38.7% intronic, 1.2% exonic, and 0.6% mirtron), organization (55.5% non-clustered and 44.5% clustered), and conservation (49% highly conserved, 14% conserved and 37% poorly conserved). Clustered miRNAs and highly conserved miRNAs were more highly expressed (p < 0.05) and had more predicted targets than non-clustered or less conserved miRNAs (p < 0.001). A total of 34 miRNAs were coordinately expressed, being part of six identified relevant networks. Two intronic miRNAs (miR-33a and miR-1281) were confirmed to have coordinated expression with their host genes, transcriptional factor SREBF2 and EP300 (a transcriptional co-activator of transcriptional factor C/EBPα), respectively which are involved in lipid metabolism, suggesting these miRNAs may also play a role in regulation of bovine lipid metabolism/adipogenesis. Furthermore, a total of 17 bovine specific miRNAs were predicted to be involved in the regulation of energy balance in adipose tissue.

Conclusions

These findings improve our understanding on the behavior of miRNAs in the regulation of bovine adipogenesis and fat metabolism as it reveals that miRNA expression patterns and functions are associated with miRNA genomic location, organization and conservation.

Expression variants of the lipogenic AGPAT6 gene affect diverse milk composition phenotypes in Bos taurus

Milk is composed of a complex mixture of lipids, proteins, carbohydrates and various vitamins and minerals as a source of nutrition for young mammals. The composition of milk varies between individuals, with lipid composition in particular being highly heritable. Recent reports have highlighted a region of bovine chromosome 27 harbouring variants affecting milk fat percentage and fatty acid content. We aimed to further investigate this locus in two independent cattle populations, consisting of a Holstein-Friesian x Jersey crossbreed pedigree of 711 F2 cows, and a collection of 32,530 mixed ancestry Bos taurus cows. Bayesian genome-wide association mapping using markers imputed from the Illumina BovineHD chip revealed a large quantitative trait locus (QTL) for milk fat percentage on chromosome 27, present in both populations. We also investigated a range of other milk composition phenotypes, and report additional associations at this locus for fat yield, protein percentage and yield, lactose percentage and yield, milk volume, and the proportions of numerous milk fatty acids. We then used mammary RNA sequence data from 212 lactating cows to assess the transcript abundance of genes located in the milk fat percentage QTL interval. This analysis revealed a strong eQTL for AGPAT6, demonstrating that high milk fat percentage genotype is also additively associated with increased expression of the AGPAT6 gene. Finally, we used whole genome sequence data from six F1 sires to target a panel of novel AGPAT6 locus variants for genotyping in the F2 crossbreed population. Association analysis of 58 of these variants revealed highly significant association for polymorphisms mapping to the 5′UTR exons and intron 1 of AGPAT6. Taken together, these data suggest that variants affecting the expression of AGPAT6 are causally involved in differential milk fat synthesis, with pleiotropic consequences for a diverse range of other milk components.

Association analysis for feed efficiency traits in beef cattle using preserved haplotypes

This study reports a genome wide scan for chromosome regions and their haplotypes that significantly associated with average daily gain (ADG), dry matter intake (DMI), and residual feed intake (RFI) in beef cattle. The study used data from 597 Angus, 450 Charolais, and 616 crossbred beef cattle, and the Illumina Bovine SNP50 beadchip. Extended haplotype homozygosity was used to identify chromosome regions that had been recently selected for in the three groups of animals. Such regions in the crossbreds were tested for association with ADG, DMI, and RFI. At false discovery rates of 5% and 10%, there were six and eight chromosome regions showing significant associations with the traits, respectively. At nominal significance levels (at least P < 0.05), 23 regions with a total number of 31 haplotypes were found significantly associated with at least one of the three traits. The proportion of phenotypic variance explained by these 23 regions varied depending on the trait; the highest proportion for ADG, DMI, and RFI was 13.50%, 9.92%, and 2.64%, respectively. Most of the haplotypes affected single traits, except for GAA (BTA4), GCG (BTA7), and TAGT (BTA12) that affected multiple traits. Thirty-six quantitative trait loci for 16 production traits, from the current literature, covered fully or in part the 23 chromosome regions. The findings from this study might be an important contribution to the current knowledge of the beef cattle genome and to the effective identification of causative genes associated with important traits in cattle.

Genome-wide analysis of zygotic linkage disequilibrium and its components in crossbred cattle

Background

Linkage disequilibrium (LD) between genes at linked or independent loci can occur at gametic and zygotic levels known asgametic LD and zygotic LD, respectively. Gametic LD is well known for its roles in fine-scale mapping of quantitative trait loci, genomic selection and evolutionary inference. The less-well studied is the zygotic LD and its components that can be also estimated directly from the unphased SNPs.

Results

This study was set up to investigate the genome-wide extent and patterns of zygotic LD and its components in a crossbred cattle population using the genomic data from the Illumina BovineSNP50 beadchip. The animal population arose from repeated crossbreeding of multiple breeds and selection for growth and cow reproduction. The study showed that similar genomic structures in gametic and zygotic LD were observed, with zygotic LD decaying faster than gametic LD over marker distance. The trigenic and quadrigenic disequilibria were generally two- to three-fold smaller than the usual digenic disequilibria (gametic or composite LD). There was less power of testing for these high-order genic disequilibria than for the digenic disequilibria. The power estimates decreased with the marker distance between markers though the decay trend is more obvious for the digenic disequilibria than for high-order disequilibria.

Conclusions

This study is the first major genome-wide survey of all non-allelic associations between pairs of SNPs in a cattle population. Such analysis allows us to assess the relative importance of gametic LD vs. all other non-allelic genic LDs regardless of whether or not the population is in HWE. The observed predominance of digenic LD (gametic or composite LD) coupled with insignificant high-order trigenic and quadrigenic disequilibria supports the current intensive focus on the use of high-density SNP markers for genome-wide association studies and genomic selection activities in the cattle population.