Genome‐wide association studies for hematological traits in swine

Improving immune capacity may increase the profitability of animal production if it enables animals to better cope with infections. Hematological traits play pivotal roles in animal immune capacity and disease resistance. Thus far, few studies have been conducted using a high‐density swine SNP chip panel to unravel the genetic mechanism of the immune capability in domestic animals. In this study, using mixed model‐based single‐locus regression analyses, we carried out genome‐wide association studies, using the Porcine SNP60 BeadChip, for immune responses in piglets for 18 hematological traits (seven leukocyte traits, seven erythrocyte traits, and four platelet traits) after being immunized with classical swine fever vaccine. After adjusting for multiple testing based on permutations, 10, 24, and 77 chromosome‐wise significant SNPs were identified for the leukocyte traits, erythrocyte traits, and platelet traits respectively, of which 10 reached genome‐wise significance level. Among the 53 SNPs for mean platelet volume, 29 are located in a linkage disequilibrium block between 32.77 and 40.59 Mb on SSC6. Four genes of interest are located within the block, providing genetic evidence that this genomic segment may be considered a candidate region relevant to the platelet traits. Other candidate genes of interest for red blood cell, hemoglobin, and red blood cell volume distribution width also have been found near the significant SNPs. Our genome‐wide association study provides a list of significant SNPs and candidate genes that offer valuable information for future dissection of molecular mechanisms regulating hematological traits.     

A genome-wide association study for quantitative trait loci of show-jumping in Hanoverian warmblood horses

Show-jumping is an economically important breeding goal in Hanoverian warmblood horses. The aim of this study was a genome-wide association study (GWAS) for quantitative trait loci (QTL) for show-jumping in Hanoverian warmblood horses, employing the Illumina equine SNP50 Beadchip. For our analyses, we genotyped 115 stallions of the National State stud of Lower Saxony. The show-jumping talent of a horse includes style and ability in free-jumping. To control spurious associations based on population stratification, two different mixed linear animal model (MLM) approaches were employed, besides linear models with fixed effects only and adaptive permutations for correcting multiple testing. Population stratification was explained best in the MLM considering Hanoverian, Thoroughbred, Trakehner and Holsteiner genes and the marker identity-by-state relationship matrix. We identified six QTL for show-jumping on horse chromosomes (ECA) 1, 8, 9 and 26 (-log(10) P-value >5) and further putative QTL with -log(10) P-values of 3-5 on ECA1, 3, 11, 17 and 21. Within six QTL regions, we identified human performance-related genes including PAPSS2 on ECA1, MYL2 on ECA8, TRHR on ECA9 and GABPA on ECA26 and within the putative QTL regions NRAP on ECA1, and TBX4 on ECA11. The results of our GWAS suggest that genes involved in muscle structure, development and metabolism are crucial for elite show-jumping performance. Further studies are required to validate these QTL in larger data sets and further horse populations.     

A genome‐wide association analysis for carcass traits in a commercial Duroc pig population

We performed a genome‐wide association study to map the genetic determinants of carcass traits in 350 Duroc pigs typed with the Porcine SNP60 BeadChip. Association analyses were carried out using the gemma software. The proportion of phenotypic variance explained by the SNPs ranged between negligible to moderate (= 0.01–0.30) depending on the trait under consideration. At the genome‐wide level, we detected one significant association between backfat thickness between the 3^rd and 4^th ribs and six SNPs mapping to SSC12 (37–40 Mb). We also identified several chromosome‐wide significant associations for ham weight (SSC11: 51–53 Mb, three SNPs; 67–68 Mb, two SNPs), carcass weight (SSC11: 66–68 Mb, two SNPs), backfat thickness between the 3^rd and 4^th ribs (SSC12: 21 Mb, one SNP; 33–40 Mb, 17 SNPs; 51–58 Mb, two SNPs), backfat thickness in the last rib (SSC12: 37 Mb, one SNP; 40–41 Mb, nine SNPs) and lean meat content (SSC13: 34 Mb, three SNPs and SSC16: 45.1 Mb, one SNP; 62–63 Mb, 10 SNPs; 71–75 Mb, nine SNPs). The ham weight trait‐associated region on SSC11 contains two genes (UCHL3 and LMO7) related to muscle development. In addition, the ACACA gene, which encodes an enzyme for the catalysis of fatty acid synthesis, maps to the SSC12 (37–41 Mb) region harbouring trait‐associated regions for backfat thickness traits. Sequencing of these candidate genes may help to uncover the causal mutations responsible for the associations found in the present study.     

Genome‐wide association mapping for female fertility traits in Danish and Swedish Holstein cattle

A genome‐wide association study was conducted using a mixed model analysis for QTL for fertility traits in Danish and Swedish Holstein cattle. The analysis incorporated 2,531 progeny tested bulls, and a total of 36 387 SNP markers on 29 bovine autosomes were used. Eleven fertility traits were analyzed for SNP association. Furthermore, mixed model analysis was used for association analyses where a polygenic effect was fitted as a random effect, and genotypes at single SNPs were successively included as a fixed effect in the model. The Bonferroni correction for multiple testing was applied to adjust the significance threshold. Seventy‐four SNP‐trait combinations showed chromosome‐wide significance, and five of these were significant genome‐wide. Twenty‐four QTL regions on 14 chromosomes were detected. Strong evidence for the presence of QTL that affect fertility traits were observed on chromosomes 3, 5, 10, 13, 19, 20, and 24. The QTL intervals were generally smaller than those described in earlier linkage studies. The identification of fertility trait‐associated SNPs and mapping of the corresponding QTL in small chromosomal regions reported here will facilitate searches for candidate genes and candidate polymorphisms.     

Genome‐wide association study for jumping performances in French sport horses

A genome‐wide association study was performed to identify single nucleotide polymorphisms (SNPs) associated with jumping performances of warmbloods in France. The 999 horses included in the study for jumping performances were sport horses [mostly Selle Français (68%), Anglo‐Arabians (13%) and horses from the other European studbooks]. Horses were genotyped using the Illumina EquineSNP50 BeadChip. Of the 54 602 SNPs available on this chip, 44 424 were retained after quality testing. Phenotypes were obtained by deregressing official breeding values for jumping competitions to use all available information, that is, the performances of each horse as well as those of its relatives. Two models were used to test the effects of the genotypes on deregressed phenotypes: a single‐marker mixed model and a haplotype‐based mixed model (significant: P < 1E‐05; suggestive: P < 1E‐04). Both models included a polygenic effect to take into account familial structures. For jumping performances, one suggestive quantitative trait locus (QTL) located on chromosome 1 (BIEC2_31196 and BIEC2_31198) was detected with both models. This QTL explains 0.7% of the phenotypic variance. RYR2, a gene encoding a major calcium channel in cardiac muscle in humans and mice, is located 0.55 Mb from this potential QTL.     

A genome‐wide search for quantitative trait loci affecting the cortical surface area and thickness of Heschl's gyrus

Heschl's gyrus (HG) is a core region of the auditory cortex whose morphology is highly variable across individuals. This variability has been linked to sound perception ability in both speech and music domains. Previous studies show that variations in morphological features of HG, such as cortical surface area and thickness, are heritable. To identify genetic variants that affect HG morphology, we conducted a genome‐wide association scan (GWAS) meta‐analysis in 3054 healthy individuals using HG surface area and thickness as quantitative traits. None of the single nucleotide polymorphisms (SNPs) showed association P values that would survive correction for multiple testing over the genome. The most significant association was found between right HG area and SNP rs72932726 close to gene DCBLD2 (3q12.1; P = 2.77 × 10^?7). This SNP was also associated with other regions involved in speech processing. The SNP rs333332 within gene KALRN (3q21.2; P = 2.27 × 10^?6) and rs143000161 near gene COBLL1 (2q24.3; P = 2.40 × 10^?6) were associated with the area and thickness of left HG, respectively. Both genes are involved in the development of the nervous system. The SNP rs7062395 close to the X‐linked deafness gene POU3F4 was associated with right HG thickness (Xq21.1; P = 2.38 × 10^?6). This is the first molecular genetic analysis of variability in HG morphology.     

Genome‐wide association analysis reveals novel loci for hypoxia adaptability in Tibetan chicken

The Tibetan chicken (TBC), an indigenous chicken breed of the Tibetan Plateau, has adapted to its hypoxic, high‐altitude environment over hundreds of years. The objective of this study was to identify the polymorphisms and genes associated with adaptation to hypoxia in this chicken breed. In the present study, samples were collected during days 18–21 of the incubation period from both surviving chicks and dead embryos, all of which were hatched under hypoxic conditions. A genome‐wide association study was conducted using the Illumina iSelect 60K SNP array with a case–control design, in which the case group consisted of the dead chicken embryos (n = 54) and controls were the surviving chicks (n = 82). Four significant SNPs were detected at the genome‐wide level (P < 0.05), and the results indicated that fork head box G1 (FOXG1) was the main candidate gene. The lead SNP NC_006092.4:g.33368893T>C was confirmed with a polymerase chain reaction‐restriction fragment length polymorphism analysis of 122 cases and 212 controls. A chi‐square test showed a significant association between NC_006092.4:g.33368893T>C and hatchability under hypoxic conditions (P < 0.01). Our results revealed novel polymorphisms and a candidate gene associated with hypoxic adaptation, facilitating further study in this field.     

Susceptibility loci for metabolic syndrome and metabolic components identified in Han Chinese: a multi‐stage genome‐wide association study

Metabolic syndrome (MetS), a cluster of metabolic disturbances that increase the risk for cardiovascular disease and diabetes, was because of genetic susceptibility and environmental risk factors. To identify the genetic variants associated with MetS and metabolic components, we conducted a genome‐wide association study followed by replications in totally 12,720 participants from the north, north‐eastern and eastern China. In combined analyses, independent of the top known signal at rs651821 on APOA5, we newly identified a secondary triglyceride‐associated signal at rs180326 on BUD13 (P_combined = 2.4 × 10⁻⁸). Notably, by an integrated analysis of the genotypes and the serum levels of APOA5, BUD13 and triglyceride, we observed that BUD13 was another potential mediator, besides APOA5, of the association between rs651821 and serum triglyceride. rs671 (ALDH2), an east Asian‐specific common variant, was found to be associated with MetS (P_combined = 9.7 × 10⁻²²) in Han Chinese. The effects of rs671 on metabolic components were more prominent in drinkers than in non‐drinkers. The replicated loci provided information on the genetic basis and mechanisms of MetS and metabolic components in Han Chinese.     

A genome‐wide association study of copy number variations with umbilical hernia in swine

Umbilical hernia (UH) is one of the most common congenital defects in pigs, leading to considerable economic loss and serious animal welfare problems. To test whether copy number variations (CNVs) contribute to pig UH, we performed a case–control genome‐wide CNV association study on 905 pigs from the Duroc, Landrace and Yorkshire breeds using the Porcine SNP60 BeadChip and penncnv algorithm. We first constructed a genomic map comprising 6193 CNVs that pertain to 737 CNV regions. Then, we identified eight CNVs significantly associated with the risk for UH in the three pig breeds. Six of seven significantly associated CNVs were validated using quantitative real‐time PCR. Notably, a rare CNV (CNV14:13030843–13059455) encompassing the NUGGC gene was strongly associated with UH (permutation‐corrected P = 0.0015) in Duroc pigs. This CNV occurred exclusively in seven Duroc UH‐affected individuals. SNPs surrounding the CNV did not show association signals, indicating that rare CNVs may play an important role in complex pig diseases such as UH. The NUGGC gene has been implicated in human omphalocele and inguinal hernia. Our finding supports that CNVs, including the NUGGC CNV, contribute to the pathogenesis of pig UH.     

Meta-analysis of results from quantitative trait loci mapping studies on pig chromosome 4

Meta-analysis of results from multiple studies could lead to more precise quantitative trait loci (QTL) position estimates compared to the individual experiments. As the raw data from many different studies are not readily available, the use of results from published articles may be helpful. In this study, we performed a meta-analysis of QTL on chromosome 4 in pig, using data from 25 separate experiments. First, a meta-analysis was performed for individual traits: average daily gain and backfat thickness. Second, a meta-analysis was performed for the QTL of three traits affecting loin yield: loin eye area, carcass length and loin meat weight. Third, 78 QTL were selected from 20 traits that could be assigned to one of three broad categories: carcass, fatness or growth traits. For each analysis, the number of identified meta-QTL was smaller than the number of initial QTL. The reduction in the number of QTL ranged from 71% to 86% compared to the total number before the meta-analysis. In addition, the meta-analysis reduced the QTL confidence intervals by as much as 85% compared to individual QTL estimates. The reduction in the confidence interval was greater when a large number of independent QTL was included in the meta-analysis. Meta-QTL related to growth and fatness were found in the same region as the FAT1 region. Results indicate that the meta-analysis is an efficient strategy to estimate the number and refine the positions of QTL when QTL estimates are available from multiple populations and experiments. This strategy can be used to better target further studies such as the selection of candidate genes related to trait variation.     

A genome‐wide association study suggests new candidate genes for milk production traits in Chinese Holstein cattle

A genome‐wide association study (GWAS) was conducted on 15 milk production traits in Chinese Holstein. The experimental population consisted of 445 cattle, each genotyped by the GGP (GeneSeek genomic profiling)‐BovineLD V3 SNP chip, which had 26 151 public SNPs in its manifest file. After data cleaning, 20 326 SNPs were retained for the GWAS. The phenotypes were estimated breeding values of traits, provided by a public dairy herd improvement program center that had been collected once a month for 3 years. Two statistical models, a fixed‐effect linear regression model and a mixed‐effect linear model, were used to estimate the association effects of SNPs on each of the phenotypes. Genome‐wide significant and suggestive thresholds were set at 2.46E‐06 and 4.95E‐05 respectively. The two statistical models concurrently identified two genome‐wide significant (P < 0.05) SNPs on milk production traits in this Chinese Holstein population. The positional candidate genes, which were the ones closest to these two identified SNPs, were EEF2K (eukaryotic elongation factor 2 kinase) and KLHL1 (kelch like family member 1). These two genes could serve as new candidate genes for milk yield and lactation persistence, yet their roles need to be verified in further function studies.     

Genetic‐based dissection of arsenic accumulation in maize using a genome‐wide association analysis method

Understanding the mechanism of arsenic (As) accumulation in plants is important in reducing As's toxicity to plants and its potential risks to human health. Here, we performed a genome‐wide association study to dissect the genetic basis of the As contents of different maize tissues in Xixian, which was irrigated with As‐rich surface water, and Changge using an association population consisting of 230 representative maize inbred lines. Phenotypic data revealed a wide normal distribution and high repeatability for the As contents in maize tissues. The As concentrations in maize tissues followed the same trend in the two locations: kernels < axes < stems < bracts < leaves. In total, 15, 16 and 15 non‐redundant quantitative trait loci (QTL s) associated with As concentrations were identified (P  ≤ 2.04 × 10^?6) in five tissues from Xixian, Changge, and the combination of the locations, respectively, explaining 9.70%–24.65% of the phenotypic variation for each QTL , on average. Additionally, four QTL s [involving 15 single nucleotide polymorphisms (SNP s)] were detected in the single and the combined locations, indicating that these loci/SNP s might be stable across different environments. The candidate genes associated with these four loci were predicted. In addition, four non‐redundant QTL s (6 SNP s), including a QTL that was detected in multiple locations according to the genome‐wide association study, were found to co‐localize with four previously reported QTL intervals. These results are valuable to understand the genetic architecture of As mechanism in maize and facilitate the genetic improvement of varieties without As toxicity.     

A genome‐wide scan for quantitative trait loci affecting respiratory disease and immune capacity in Landrace pigs

Respiratory disease is the most important health concern for the swine industry. Genetic improvement for disease resistance is challenging because of the difficulty in obtaining good phenotypes related with disease resistance; however, identification of genes or markers associated with disease resistance can help in the genetic improvement of pig health. The purpose of our study was to investigate whether quantitative trait loci (QTL) associated with disease resistance were segregated in a purebred population of Landrace pigs that had been selected for meat production traits and mycoplasmal pneumonia of swine (MPS) scores over five generations. We analysed 1395 pigs from the base to the fifth generation of this population. Two respiratory disease traits [MPS scores and atrophic rhinitis (AR) scores] and 11 immune‐capacity traits were measured in 630–1332 animals at 7 weeks of age and when the animal's body weight reached 105 kg. Each of the pigs, except sires in the base population, was genotyped using 109 microsatellite markers, and then, QTL analysis of the full‐sib family population with a multi‐generational pedigree structure was performed. Variance component analysis was used to detect QTL associated with MPS or AR scores, and the logarithm of odds (LOD) score and genotypic heritability of the QTL were estimated. Five significant (LOD > 2.51) and 18 suggestive (LOD > 1.35) QTL for respiratory disease traits and immune‐capacity traits were detected. The significant QTL for Log‐MPS score, located on S. scrofa chromosome 2, could explain 87% of the genetic variance of this score in this analysis. This is the first report of QTL associated with respiratory disease lesions.     

Genome-wide association mapping in Norwegian Red cattle identifies quantitative trait loci for fertility and milk production on BTA12

Reproductive performance is a critical trait in dairy cattle. Poor reproductive performance leads to prolonged calving intervals, higher culling rates and extra expenses related to multiple inseminations, veterinary treatments and replacements. Genetic gain for improved reproduction through traditional selection is often slow because of low heritability and negative correlations with production traits. Detection of DNA markers associated with improved reproductive performance through genome-wide association studies could lead to genetic gain that is more balanced between fertility and production. Norwegian Red cattle are well suited for such studies, as very large numbers of detailed reproduction records are available. We conducted a genome-wide association study for non-return rate, fertility treatments and retained placenta using almost 1 million records on these traits and 17 343 genome-wide single-nucleotide polymorphisms. Genotyping costs were minimized by genotyping the sires of the cows recorded and by using daughter averages as phenotypes. The genotyped sires were assigned to either a discovery or a validation population. Associations were only considered to be validated if they were significant in both groups. Strong associations were found and validated on chromosomes 1, 5, 8, 9, 11 and 12. Several of these were highly supported by findings in other studies. The most important result was an association for non-return rate in heifers in a region of BTA12 where several associations for milk production traits have previously been found. Subsequent fine-mapping verified the presence of a quantitative trait loci (QTL) having opposing effects on non-return rate and milk production at 18 Mb. The other reproduction QTL did not have pleiotropic effects on milk production, and these are therefore of considerable interest for use in marker-assisted selection.     

Genome‐wide association study identifies quantitative trait loci affecting hematological traits in an F2 intercross between Landrace and Korean native pigs

Changes affecting the status of health and robustness can bring about physiological alterations including hematological parameters in swine. To identify quantitative trait loci (QTL) associated with eight hematological traits (one leukocyte trait, six erythrocyte traits and one platelet trait), we conducted a genome‐wide association study using the PorcineSNP60K BeadChip in a resource population derived from an intercross between Landrace and Korean native pigs. A total of 36 740 SNPs from 816 F₂ progeny were analyzed for each blood‐related trait after filtering for quality control. Data were analyzed by the genome‐wide rapid association using mixed model and regression (GRAMMAR) approach. A total of 257 significant SNPs (P < 1.36 × 10^?6) on SSC3, 6, 8, 13 and 17 were identified for blood‐related traits in this study. Interestingly, the genomic region between 17.9 and 130 Mb on SSC8 was found to be significantly associated with red blood cell, mean corpuscular volume and mean corpuscular hemoglobin. Our results include the identification of five significant SNPs within five candidate genes (KIT, IL15, TXK, ARAP2 and ERG) for hematopoiesis. Further validation of these identified SNPs could give valuable information for understanding the variation of hematological traits in pigs.     

Meta‐analysis of genome‐wide association from genomic prediction models

Genome‐wide association (GWA) studies based on GBLUP models are a common practice in animal breeding. However, effect sizes of GWA tests are small, requiring larger sample sizes to enhance power of detection of rare variants. Because of difficulties in increasing sample size in animal populations, one alternative is to implement a meta‐analysis (MA), combining information and results from independent GWA studies. Although this methodology has been used widely in human genetics, implementation in animal breeding has been limited. Thus, we present methods to implement a MA of GWA, describing the proper approach to compute weights derived from multiple genomic evaluations based on animal‐centric GBLUP models. Application to real datasets shows that MA increases power of detection of associations in comparison with population‐level GWA, allowing for population structure and heterogeneity of variance components across populations to be accounted for. Another advantage of MA is that it does not require access to genotype data that is required for a joint analysis. Scripts related to the implementation of this approach, which consider the strength of association as well as the sign, are distributed and thus account for heterogeneity in association phase between QTL and SNPs. Thus, MA of GWA is an attractive alternative to summarizing results from multiple genomic studies, avoiding restrictions with genotype data sharing, definition of fixed effects and different scales of measurement of evaluated traits.     

A genome scan for quantitative trait loci affecting three ear traits in a White Duroc × Chinese Erhualian resource population

Chinese Erhualian pigs have larger and floppier ears compared with White Duroc pigs (small, half- or fully-pricked ears). To identify quantitative trait loci (QTL) for ear weight and area as well as erectness, a genome-wide scan with 194 microsatellites was performed in a White Duroc × Chinese Erhualian resource population (>1000 F₂ animals). Twenty-three genome-wide significant QTL and 12 suggestive QTL were identified. All QTL for ear erectness and size detected in two previous studies, bar two on SSC6 and 9, were confirmed here. The 1% genome-wide significant QTL at 70 cM on SSC5 and at 58 cM on SSC7 have profound and pleiotropic effects on the three ear traits, with Erhualian alleles increasing weight and area but decreasing erectness. Notably, the 95% confidence interval of the QTL for weight and area on SSC7 spanned only 3 cM. New QTL reaching 1% genome-wide significance were found on SSC8 (at 37 cM) for all three ear traits, on SSC4 and 16 for weight and area, and on SSCX for area. Unexpectedly, Erhualian alleles at these loci were associated with lighter and smaller or erect ear. Some new suggestive QTL were also found on other chromosome regions. Almost all the QTL for weight and area had essentially additive effects, while the QTL for erectness on SSC2, 5 and 7 showed not only additive effects but also partial dominance effects of Erhualian alleles. The two most significant QTL on SSC7 and SSC5 could be promising targets for fine mapping and identification of the causative mutations.     

A genome‐wide association study for equine recurrent airway obstruction in European Warmblood horses reveals a suggestive new quantitative trait locus on chromosome 13

Recurrent airway obstruction (RAO), also known as heaves, is an asthma‐like respiratory disease. Its development is strongly influenced by environmental risk factors such as sensitization and exposure to moldy hay, straw bedding and stabling indoors. A hereditary component has been documented in previous studies; however, so far no causative genetic variant that influences the risk of developing RAO has been identified. In this study, we revised an existing dataset and selected 384 horses for genotyping on the Affymetrix high‐density equine SNP array. We performed an allelic case–control genome‐wide association study, which revealed a suggestively significant association on equine chromosome 13 at 32 843 309 bp. This SNP is located in the protein‐coding gene TXNDC11, which is possibly involved in the folding process of the multiprotein complexes DUOX1 and DUOX2. In humans, these proteins are known to take part in regulating the production of H₂O₂ in the respiratory tract epithelium as well as in MUC5AC mucin expression. Therefore, TXNDC11 may be considered a functional candidate gene, and further research is needed to explore its potential role in RAO‐affected horses.     

Genome‐wide association analysis for quantitative trait loci influencing Warner–Bratzler shear force in five taurine cattle breeds

We performed a genome‐wide association study for Warner–Bratzler shear force (WBSF), a measure of meat tenderness, by genotyping 3360 animals from five breeds with 54 790 BovineSNP50 and 96 putative single‐nucleotide polymorphisms (SNPs) within μ‐calpain [HUGO nomenclature calpain 1, (mu/I) large subunit; CAPN1] and calpastatin (CAST). Within‐ and across‐breed analyses estimated SNP allele substitution effects (ASEs) by genomic best linear unbiased prediction (GBLUP) and variance components by restricted maximum likelihood under an animal model incorporating a genomic relationship matrix. GBLUP estimates of ASEs from the across‐breed analysis were moderately correlated (0.31–0.66) with those from the individual within‐breed analyses, indicating that prediction equations for molecular estimates of breeding value developed from across‐breed analyses should be effective for genomic selection within breeds. We identified 79 genomic regions associated with WBSF in at least three breeds, but only eight were detected in all five breeds, suggesting that the within‐breed analyses were underpowered, that different quantitative trait loci (QTL) underlie variation between breeds or that the BovineSNP50 SNP density is insufficient to detect common QTL among breeds. In the across‐breed analysis, CAPN1 was followed by CAST as the most strongly associated WBSF QTL genome‐wide, and associations with both were detected in all five breeds. We show that none of the four commercialized CAST and CAPN1 SNP diagnostics are causal for associations with WBSF, and we putatively fine‐map the CAPN1 causal mutation to a 4581‐bp region. We estimate that variation in CAST and CAPN1 explains 1.02 and 1.85% of the phenotypic variation in WBSF respectively.