Search

Replication of porcine circovirus type 2 (PCV2), an important worldwide swine pathogen, has been demonstrated to be influenced by host genotype. Specifically, a missense DNA polymorphism (SYNGR2 p.Arg63Cys) within the SYNGR2 gene was demonstrated to contribute to variation in PCV2b viral load and subsequent immune response following infection. PCV2 is known to induce immunosuppression leading to an increase in susceptibility to subsequent infections with other viral pathogens such as porcine reproductive and respiratory syndrome virus (PRRSV). In order to assess the role of SYNGR2 p.Arg63Cys in co-infections, pigs homozygous for the favorable SYNGR2 p.63Cys (N = 30) and unfavorable SYNGR2 p.63Arg (N = 29) alleles were infected with PCV2b followed a week later by a challenge with PRRSV. A lower PCV2b viremia (P < 0.001) and PCV2-specific IgM antibodies (P < 0.005) were observed in SYNGR2 p.63Cys compared to SYNGR2 p.63Arg genotypes. No significant differences in PRRSV viremia and specific IgG antibodies were observed between SYNGR2 genotypes. Lung histology score, an indicator of disease severity, was lower in the pigs with SYNGR2 p.63Cys genotypes (P < 0.05). Variation in the lung histology scores within SYNGR2 genotypes suggests that additional factors, environmental and/or genetic, could be involved in disease severity. Lay Summary: Porcine circovirus type 2 (PCV2) is an important virus involved in the onset of a group of severe disease symptoms commonly known as porcine circovirus associated diseases (PCVAD). Vaccination options exist for PCV2, though the severity of PCVAD can be influenced by the presence of additional co-infecting pathogens, such as porcine reproductive and respiratory syndrome virus (PRRSV), for which vaccination is still a challenge. Host genetic resistance is a potential avenue for solving this problem. Previously, a genetic polymorphism in the SYNGR2 gene was found to be associated with PCV2b viremia and immune response. The aim of this stu

The BXD family has become one of the preeminent genetic reference populations to understand the genetic and environmental control of phenotypic variation. Here we evaluate the responses to different levels of fat in the diet using both chow diet (CD, 13-18% fat) and a high-fat diet (HFD, 45-60% fat). We studied cohorts of BXD strains, both inbred parents C57BL/6J and DBA/2J (commonly known as B6 and D2, respectively), as well as B6D2 and D2B6 reciprocal F1 hybrids. The comparative impact of genetic and dietary factors was analyzed by profiling a range of phenotypes, most prominently their cecum bacterial composition. The parents of the BXDs and F1 hybrids express limited differences in terms of weight and body fat gain on CD. In contrast, the strain differences on HFD are substantial for percent body fat, with DBA/2J accumulating 12.5% more fat than C57BL/6J (P < 0.0001). The F1 hybrids born to DBA/2J dams (D2B6F1) have 10.6% more body fat (P < 0.001) than those born to C57BL/6J dams. Sequence analysis of the cecum microbiota reveals important differences in bacterial composition among BXD family members with a substantial shift in composition caused by HFD. Relative to CD, the HFD induces a decline in diversity at the phylum level with a substantial increase in Firmicutes (+13.8%) and a reduction in Actinobacteria (-7.9%). In the majority of BXD strains, the HFD also increases cecal sIgA (P < 0.0001)-an important component of the adaptive immunity response against microbial pathogens. Host genetics modulates variation in cecum bacterial composition at the genus level in CD, with significant quantitative trait loci (QTLs) for Oscillibacter mapped to Chr 3 (18.7-19.2 Mb, LRS = 21.4) and for Bifidobacterium mapped to Chr 6 (89.21-89.37 Mb, LRS = 19.4). Introduction of HFD served as an environmental suppressor of these QTLs due to a reduction in the contribution of both genera (P < 0.001). Relations among liver metabolites and cecum bacterial composition were predomina

Understanding early predictors of sow fertility has the potential to improve genomic predictions. A custom SNP array (SowPro90 produced by Affymetrix) was developed to include genetic variants overlapping quantitative trait loci for age at puberty, one of the earliest indicators of sow fertility, as well as variants related to innate and adaptive immunity. The polymorphisms included in the custom genotyping array were identified using multiple genomic approaches including deep genomic and transcriptomic sequencing and genome-wide associations. Animals from research and commercial populations (n = 2,586) were genotyped for 103,476 SNPs included in SowPro90. To assess the quality of data generated, genotype concordance was evaluated between the SowPro90 and Porcine SNP60 BeadArray using a subset of common SNP (n = 44,708) and animals (n = 277). The mean genotype concordance rate per SNP was 98.4%. Differences in distribution of data quality were observed between the platforms indicating the need for platform specific thresholds for quality parameters. The optimal thresholds for SowPro90 (≥97% SNP and ≥93% sample call rate) were obtained by analyzing the data quality distribution and genotype concordance per SNP across platforms. At ≥97% SNP call rate, there were 42,151 SNPs (94.3%) retained with a mean genotype concordance of 98.6% across platforms. Similarly, ≥94% SNPs and ≥85% sample call rates were established as thresholds for Porcine SNP60 BeadArray. At ≥94% SNPs call rate, there were 41,043 SNPs (91.8%) retained with a mean genotype concordance of 98.6% across platforms. Final evaluation of SowPro90 array content (n = 103,476) at ≥97% SNPs and ≥93% sample call rates allowed retention of 89,040 SNPs (86%) for downstream analysis. The findings and strategy for quality control could be helpful in identifying consistent, high-quality genotypes for genomic evaluations, especially when integrating genotype data from different platforms.

Approximately 50% of sows are culled annually with more than one-third due to poor fertility. Our research demonstrated that age at puberty is an early pre-breeding indicator of reproductive longevity. Age at puberty can be measured early in life, has a moderate heritability, and is negatively correlated with lifetime number of parities. Detection of age at puberty is tedious and time consuming and is therefore not collected by the industry, which limits genetic progress. Genomic prediction is a viable approach to preselect gilts that will express puberty early and have superior reproductive longevity. The hypothesis that genetic variants explaining differences in age at puberty also explain differences in sow reproductive longevity was tested. Phenotypes, genotypes, and tissues from the UNL resource population (n > 1700) were used in genome-wide association analyses, genome, and RNA sequencing to uncover functional polymorphisms that could explain variation in puberty and reproductive longevity. A BeadArray including 56,424 SNP explained 25.2% of the phenotypic variation in age at puberty in a training set (n = 820). Evaluation of major windows and SNPs of subsequent batches of similar genetics (n = 412) showed that if all SNPs located in the major 1-Mb windows were tested, they explained a substantial amount of phenotypic variation (12.3 to 36.8%). Due to differences in linkage disequilibrium status, the most informative SNP from these windows explained a lower proportion of the variation (6.5 to 23.7%). To improve genomic predictive ability, the limited capability of BeadArray was enhanced by potential functional variants uncovered by genome sequencing of selected sires (n = 20; >20X). There were 11.2 mil. SNPs and 2.9 mil. indels discovered across sires and reference genomes. The role of gene expression differences in explaining phenotypic variation in age at puberty was investigated by RNA sequencing of the hypothalamic arcuate nucleus (ARC) in gilts (n = 37) w

Genetic variants associated with traits such as age at puberty and litter size could provide insight into the underlying genetic sources of variation impacting sow reproductive longevity and productivity. Genomewide characterization and gene expression profiling were used using gilts from the University of Nebraska–Lincoln swine resource population (n = 1,644) to identify genetic variants associated with age at puberty and litter size traits. From all reproductive traits studied, the largest fraction of phenotypic variation explained by the Porcine SNP60 BeadArray was for age at puberty (27.3%). In an evaluation data set, the predictive ability of all SNP from highranked 1-Mb windows (1 to 50%), based on genetic variance explained in training, was greater (12.3 to 36.8%) compared with the most informative SNP from these windows (6.5 to 23.7%). In the integrated data set (n = 1,644), the top 1% of the 1-Mb windows explained 6.7% of the genetic variation of age at puberty. One of the high-ranked windows detected (SSC2, 12–12.9 Mb) showed pleiotropic features, affecting both age at puberty and litter size traits. The RNA sequencing of the hypothalami arcuate nucleus uncovered 17 differentially expressed genes (adjusted P < 0.05) between gilts that became pubertal early (180 d of age). Twelve of the differentially expressed genes are upregulated in the late pubertal gilts. One of these genes is involved in energy homeostasis (FFAR2), a function in which the arcuate nucleus plays an important contribution, linking nutrition with reproductive development. Energy restriction during the gilt development period delayed age at puberty by 7 d but increased the probability of a sow to produce up to 3 parities (P < 0.05). Identification of pleotropic functional polymorphisms may improve accuracy of genomic prediction while facilitating a reduction in sow replacement rates and addressing welfare concerns.

Phenome-wide association is a novel reverse genetic strategy to analyze genome-to-phenome relations in human clinical cohorts. Here we test this approach using a large murine population segregating for approximately 5 million sequence variants, and we compare our results to those extracted from a matched analysis of gene variants in a large human cohort. For the mouse cohort, we amassed a deep and broad open-access phenome consisting of approximately 4,500 metabolic, physiological, pharmacological and behavioural traits, and more than 90 independent expression quantitative trait locus (QTL), transcriptome, proteome, metagenome and metabolome data sets--by far the largest coherent phenome for any experimental cohort (www.genenetwork.org). We tested downstream effects of subsets of variants and discovered several novel associations, including a missense mutation in fumarate hydratase that controls variation in the mitochondrial unfolded protein response in both mouse and Caenorhabditis elegans, and missense mutations in Col6a5 that underlies variation in bone mineral density in both mouse and human.

Tumor necrosis factor alpha (TNF-α) is a pro-inflammatory cytokine with a role in activating adaptive immunity to viral infections. By inhibiting the capacity of plasmacytoid dendritic cells to produce interferon-α and TNF-α, porcine circovirus 2 (PCV2) limits the maturation of myeloid dendritic cells and impairs their ability to recognize viral and bacterial antigens. Previously, we reported QTL for viremia and immune response in PCV2- infected pigs. In this study, we analyzed phenotypic and genetic relationships between TNFα protein levels, a potential indicator of predisposition to PCV2 co-infection, and PCV2 susceptibility. Following experimental challenge with PCV2b, TNF-α reached the peak at 21 days post-infection (dpi), at which time a difference was observed between pigs that expressed extreme variation in viremia and growth (P < 0.10). A genome-wide association study (n = 297) revealed that genotypes of 56 433 SNPs explained 73.9% of the variation in TNF-α at 21 dpi. Major SNPs were identified on SSC8, SSC10 and SSC14. Haplotypes based on SNPs from a SSC8 (9 Mb) 1-Mb window were associated with variation in TNF-α (P < 0.02), IgG (P = 0.05) and IgM (P < 0.13) levels at 21 dpi. Potential overlap of regulatory mechanisms was supported by the correlations between genomic prediction values of TNF-α and PCV2 antibodies (21 dpi, r > 0.22), viremia (14–21 dpi, P > 0.29) and viral load (r = 0.31, P < 0.0001). Characterization of the QTL regions uncovered genes that could influence variation in TNF-α levels as well as T- and B-cell development, which can affect disease susceptibility.

Genomic information could be used efficiently to improve traits that are expensive to measure, sex limited or expressed late in life. This study analyzed the phenotypic variation explained by major SNPs and windows for age at puberty in gilts, an indicator of reproductive longevity. A genome-wide association study using 56,424 SNPs explained 25.2% of the phenotypic variation in age at puberty in a training set (n = 820). All SNPs from the top 10% of 1-Mb windows explained 33.5% of the phenotypic variance compared to 47.1% explained by the most informative markers (n = 261). In an evaluation population, consisting of subsequent batches (n = 412), the predictive ability of all SNPs from the major 1-Mb windows was higher compared to the variance captured by the most informative SNP from each of these windows. The phenotypic variance explained in the evaluation population varied from 12.3% to 36.8% when all SNPs from major windows were used compared to 6.5–23.7% explained by most informative SNPs. The correlation between phenotype and genomic prediction values based on SNP effects estimated in the training population was marginal compared to their effects retrained in the evaluation population for all (0.46–0.81) or most informative SNPs (0.30–0.65) from major windows. An increase in genetic gain of 20.5% could be obtained if genomic selection included both sexes compared to females alone. The pleiotropic role of major genes such as AVPR1A could be exploited in selection of both age at puberty and reproductive longevity. Two supplemental files are attached (below).

An integral part of functional genomics studies is to assess the enrichment of specific biological terms in lists of genes found to be playing an important role in biological phenomena. Contrasting the observed frequency of annotated terms with those of the background is at the core of overrepresentation analysis (ORA). Gene Ontology (GO) is a means to consistently classify and annotate gene products and has become a mainstay in ORA. Alternatively, Medical Subject Headings (MeSH) offers a comprehensive life science vocabulary including additional categories that are not covered by GO. Although MeSH is applied predominantly in human and model organism research, its full potential in livestock genetics is yet to be explored. In this study, MeSH ORA was evaluated to discern biological properties of identified genes and contrast them with the results obtained from GO enrichment analysis. Three published datasets were employed for this purpose, representing a gene expression study in dairy cattle, the use of SNPs for genome-wide prediction in swine and the identification of genomic regions targeted by selection in horses. We found that several overrepresented MeSH annotations linked to these gene sets share similar concepts with those of GO terms. Moreover, MeSH yielded unique annotations, which are not directly provided by GO terms, suggesting that MeSH has the potential to refine and enrich the representation of biological knowledge. We demonstrated that MeSH can be regarded as another choice of annotation to draw biological inferences from genes identified via experimental analyses. When used in combination with GO terms, our results indicate that MeSH can enhance our functional interpretations for specific biological conditions or the genetic basis of complex traits in livestock species.

Background: Porcine circovirus 2 is the primary agent responsible for inducing a group of associated diseases known as Porcine Circovirus Associated Diseases (PCVAD), which can have detrimental effects on production efficiency as well as causing significant mortality. The objective of this study was to evaluate variation in viral replication, immune response and growth across pigs (n = 974) from different crossbred lines. The approach used in this study was experimental infection with a PCV2b strain of pigs at an average of 43 days of age. Results: The sequence of the PCV2b isolate used in the challenge was similar with a cluster of PCV2b isolates known to induce PCVAD and increased mortality rates. The swine leukocyte antigen class II (SLAII) profile of the population was diverse, with nine DQB1 haplotypes being present. Individual viremia and antibody profiles during challenge demonstrate variation in magnitude and time of viral surge and immune response. The correlations between PCV2 specific antibodies and average daily gain (ADG) were relatively low and varied between - 0.14 to 0.08 for IgM and −0.02 and 0.11 for IgG. In contrast, PCV2 viremia was an important driver of ADG decline following infection; a moderate negative correlation was observed between viral load and overall ADG (r = − 0.35, P < 0.001). The pigs with the lowest 10% level of viral load maintained a steady increase in weekly ADG (P < 0.0001) compared to the pigs that had the 10% greatest viral load (P < 0.55). In addition, the highly viremic group expressed higher IgM and IgG starting with d 14 and d 21 respectively, and higher tumor necrosis factor – alpha (TNF-α) at d 21 (P < 0.005), compared to low viremic group. Conclusions: Molecular sources of the observed differences in viremia and immune response could provide a better understanding of the host factors that influence the development of PCVAD and lead to improved knowledge of swine immunity.

Cattle are reared in diverse environments and collecting phenotypic body temperature (BT) measurements to characterize BT variation across diverse environments is difficult and expensive. To better understand the genetic basis of BT regulation, a genome-wide association study was conducted utilizing crossbred steers and heifers totaling 239 animals of unknown pedigree and breed fraction. During predicted extreme heat and cold stress events, hourly tympanic and vaginal BT devices were placed in steers and heifers, respect i v e l y. I n d i v i dua l s we r e genotyped with t h e BovineSNP50K_v2 assay and data analyzed using Bayesian models for area under the curve (AUC), a measure of BTover time, using hourly BT observations summed across 5-days (AUC summer 5-day (AUCS5D) and AUC winter 5-day (AUCW5D)). Posterior heritability estimates were moderate to high and were estimated to be 0.68 and 0.21 for AUCS5D and AUCW5D, respectively.Moderately positive correlations between direct genomic values for AUCS5D and AUCW5D (0.40) were found, although a small percentage of the top 5 % 1-Mb windows were in common. Different sets of genes were associated with BT during winter and summer, thus simultaneous selection for animals tolerant to both heat and cold appears possible.

Traditional selection for sow reproductive longevity is ineffective due to low heritability and late expression of the trait. Incorporation of DNA markers into selection programs is potentially a more practical approach for improving sow lifetime productivity. Using a resource population of crossbred gilts, we explored pleiotropic sources of variation that influence age at puberty and reproductive longevity. Of the traits recorded before breeding, only age at puberty significantly affected the probability that females would produce a first parity litter. The genetic variance explained by 1-Mb windows of the sow genome, compared across traits, uncovered regions that influence both age at puberty and lifetime number of parities. Allelic variants of SNPs located on SSC5 (27–28 Mb), SSC8 (36–37 Mb) and SSC12 (1.2–2 Mb) exhibited additive effects and were associated with both early expression of puberty and a greater than average number of lifetime parities. Combined analysis of these SNPs showed that an increase in the number of favorable alleles had positive impact on reproductive longevity, increasing number of parities by up to 1.36. The region located on SSC5 harbors non-synonymous alleles in the arginine vasopressin receptor 1A (AVPR1A) gene, a G-protein-coupled receptor associated with social and reproductive behaviors in voles and humans and a candidate for the observed effects. This region is characterized by high levels of linkage disequilibrium in different lines and could be exploited in marker-assisted selection programs across populations to increase sow reproductive longevity. Includes supplementary materials (138 pp.)

The gastrointestinal tract harbors a complex and diverse microbiota that has an important role in host metabolism. Microbial diversity is influenced by a combination of environmental and host genetic factors and is associated with several polygenic diseases. In this study we combined next-generation sequencing, genetic mapping, and a set of physiological traits of the BXD mouse population to explore genetic factors that explain differences in gut microbiota and its impact on metabolic traits. Molecular profiling of the gut microbiota revealed important quantitative differences in microbial composition among BXD strains. These differences in gut microbial composition are influenced by host-genetics, which is complex and involves many loci. Linkage analysis defined Quantitative Trait Loci (QTLs) restricted to a particular taxon, branch or that influenced the variation of taxa across phyla. Gene expression within the gastrointestinal tract and sequence analysis of the parental genomes in the QTL regions uncovered candidate genes with potential to alter gut immunological profiles and impact the balance between gut microbial communities. A QTL region on Chr 4 that overlaps several interferon genes modulates the population of Bacteroides, and potentially Bacteroidetes and Firmicutes-the predominant BXD gut phyla. Irak4, a signaling molecule in the Toll-like receptor pathways is a candidate for the QTL on Chr15 that modulates Rikenellaceae, whereas Tgfb3, a cytokine modulating the barrier function of the intestine and tolerance to commensal bacteria, overlaps a QTL on Chr 12 that influence Prevotellaceae. Relationships between gut microflora, morphological and metabolic traits were uncovered, some potentially a result of common genetic sources of variation.

Stress is a major risk factor for numerous neuropsychiatric diseases. However, susceptibility to stress and the qualitative nature of stress effects on behavior differ markedly among individuals. This is partly because of the moderating influence of genetic factors. Inbred mouse strains provide a relatively stable and restricted range of genetic and environmental variability that is valuable for disentangling gene–stress interactions. Here, we screened a panel of inbred strains for anxiety- and depression-related phenotypes at baseline (trait) and after exposure to repeated restraint. Two strains, DBA/2J and C57BL/6J, differed in trait and restraint-induced anxiety-related behavior (dark/light exploration, elevated plus maze). Gene expression analysis of amygdala, medial prefrontal cortex, and hippocampus revealed divergent expression in DBA/2J and C57BL/6J both at baseline and after repeated restraint. Restraint produced strain-dependent expression alterations in various genes including glutamate receptors (e.g., Grin1, Grik1). To elucidate neuronal correlates of these strain differences, we performed ex vivo analysis of glutamate excitatory neurotransmission in amygdala principal neurons. Repeated restraint augmented amygdala excitatory postsynaptic signaling and altered metaplasticity (temporal summation of NMDA receptor currents) in DBA/2J but not C57BL/6J. Furthermore, we found that the C57BL/6J-like changes in anxiety-related behavior after restraint were absent in null mutants lacking the modulatory NMDA receptor subunit Grin2a, but not the AMPA receptor subunit Gria1. Grin2a null mutants exhibited significant (~30%) loss of dendritic spines on amygdala principal neurons under nonrestraint conditions. Collectively, our data support a model in which genetic variation in glutamatergic neuroplasticity in corticolimbic circuitry underlies phenotypic variation in responsivity to stress.

The DBA/2J mouse is not only the oldest inbred strain, but also one of the most widely used strains. DBA/2J exhibits many unique anatomical, physiological, and behavior traits. In addition, DBA/2J is one parent of the large BXD family of recombinant inbred strains [1]. The genome of the other parent of this BXD family— C57BL/6J—has been sequenced and serves as the mouse reference genome [2]. We sequenced the genome of DBA/2J using SOLiD and Illumina high throughput short read protocols to generate a comprehensive set of ~5 million sequence variants segregating in the BXD family that ultimately cause developmental, anatomical, functional and behavioral differences among these 80+ strains.

Purpose—The size of neuronal populations is modulated by gene variants that influence cell production and survival, in turn influencing neuronal connectivity, function, and disease risk. The size of the dopaminergic amacrine (DA) cell population is a highly heritable trait exhibiting six-fold variation among inbred strains of mice, and is used here to identify genes that modulate the number of DA cells. Methods—The entire population was counted in retinal wholemounts from 37 genetically defined lines of mice, including six standard inbred strains, 25 recombinant inbred strains (AXB/BXA), reciprocal F1 hybrids, a chromosome (Chr) 7 consomic line, and three additional genetically modified lines. Results—We mapped much of this variation to a broad locus on Chr 7 (Dopaminergic amacrine cell number control, Chr 7). The Dacnc7 locus is flanked by two candidate genes known to modulate the number of other types of retinal neuron—the pro-apoptotic gene, Bax, and tyrosinase. The Tyr mutation was shown to modulate DA cell number modestly, although in the direction opposite that predicted. In contrast, Bax deficiency increased the population four-fold. Bax expression was significantly greater in the A/J strain relative to C57BL/6J, an effect that may be due to an SNP in a p53 consensus binding site known to modulate transcription. Finally, we note a strong candidate situated at the peak of the Dacnc7 locus, Lrrk1, a Parkinson’s disease gene exhibiting mis-sense mutations segregating within the AXB/BXA cross. Conclusions—Multiple polymorphic genes on Chr. 7 modulate the size of the population of DA cells.

Microarray analysis reveals distinctive signaling between the bed nucleus of the stria terminalis, nucleus accumbens, and dorsal striatum. Physiol Genomics 32: 283–298, 2008. First published October 2, 2006; doi:10.1152/physiolgenomics.00224.2006.—To identify distinct transcriptional patterns between the major subcortical dopamine targets commonly studied in addiction we studied differences in gene expression between the bed nucleus of the stria terminalis (BNST), nucleus accumbens (NAc), and dorsal striatum (dStr) using microarray analysis. We first tested for differences in expression of genes encoding transcripts for common neurotransmitter systems as well as calcium binding proteins routinely used in neuroanatomical delineation of brain regions. This a priori method revealed differential expression of corticotropin releasing hormone (Crh), the GABA transporter (Slc6a1), and prodynorphin (Pdyn) mRNAs as well as several others. Using a gene ontology tool, functional scoring analysis, and Ingenuity Pathway Analysis, we further identified several physiological pathways that were distinct among these brain regions. These two different analyses both identified calcium signaling, G-coupled protein receptor signaling, and adenylate cyclaserelated signaling as significantly different among the BNST, NAc, and dStr. These types of signaling pathways play important roles in, amongst other things, synaptic plasticity. Investigation of differential gene expression revealed several instances that may provide insight into reported differences in synaptic plasticity between these brain regions. The results support other studies suggesting that crucial pathways involved in neurotransmission are distinct among the BNST, NAc, and dStr and provide insight into the potential use of pharmacological agents that may target region-specific signaling pathways. Furthermore, these studies provide a framework for future mouse-mouse comparisons of transcriptional profiles after behavioral/pha

Suggestive QTL affecting raw firmness scores and average Instron force, tenderness, juiciness, and chewiness on cooked meat were mapped to pig chromosome 2 using a three-generation intercross between Berkshire and Yorkshire pigs. Based on its function and location, the calpastatin (CAST) gene was considered to be a good candidate for the observed effects. Several missense and silent mutations were identified in CAST and haplotypes covering most of the coding region were constructed and used for association analyses with meat quality traits. Results demonstrated that one CAST haplotype was significantly associated with lower Instron force and cooking loss and higher juiciness and, therefore, this haplotype is associated with higher eating quality. Some of the sequence variation identified may be associated with differences in phosphorylation of CAST by adenosine cyclic 3′, 5′-monophosphate- dependent protein kinase and may in turn explain the meat quality phenotypic differences. The beneficial haplotype was present in all the commercial breeds tested and may provide significant improvements for the pig industry and consumers because it can be used in marker-assisted selection to produce naturally tender and juicy pork without additional processing steps.