Genetic polymorphisms in ESR and FSHβ genes and their association with litter traits in Large White pigs

The estrogen receptor (ESR) gene and follicle-stimulating hormone β (FSHβ) gene are responsible for litter traits. The present study aimed to verify the polymorphisms of ESR and FSHβ and assess their effects on the litter traits in 201 Large White pigs. Four SNPs (g.C669T, g.A1296G, g.C1665T and g.A1755G) were found in ESR. The TT genotype at g.C1665T locus and AA genotype at g.A1755G locus could significantly increase the total litter size of the first litter of American Large White pigs (p < 0.05). Eight SNPs were found in exon 3 of FSHβ. The AA genotype at g.A511G locus, AA and AG genotypes at g.A617G locus, CC and CT genotypes at g.C630T locus, CT and TT genotypes at g.C652T locus, CT and TT genotypes at g.C735T locus, AA and AG genotypes at g.A746G, AA and AG genotypes at g.A921G and CT genotype at g.C678T could significantly increase the litter size of different strains of Large White pigs (p < 0.05). Our study revealed that the genetic variations of ESR and FSHβ were closely related to the litter trait of Large White pigs. Therefore, ESR and FSHβ genes could be used as molecular markers for the genetic selection of Large White pigs.

Functional genomics of reproduction in pigs: Are we there yet?

Reproductive failure is the main reason for culling females in swine herds and is both a financial and sustainability issue. Because reproductive traits are complex and lowly to moderately heritable, genomic selection within populations can achieve substantial genetic gain in reproductive efficiency. A better understanding of the physiological components affecting the expression of these traits will facilitate greater understanding of the genes affecting reproductive traits and is necessary to improve and optimize management strategies to maximize reproductive success of gilts and sows. Large-scale genotyping with single-nucleotide polymorphism (SNP) arrays are used for genome-wide association studies (GWAS) and have facilitated identification of positional candidate genes. Transcriptomic data can be used to weight SNP for GWAS and could lead to previously unidentified candidate genes. Resequencing and fine mapping of candidate genes are necessary to identify putative functional variants and some of these have been incorporated into new genotyping arrays. Sequence imputation and genotype by sequence are newer strategies that could reveal novel functional mutations. In this study, these approaches are discussed. Advantages and limitations are highlighted where additional research is needed.

Utilizing Variants Identified with Multiple Genome-Wide Association Study Methods Optimizes Genomic Selection for Growth Traits in Pigs

Improving the prediction accuracies of economically important traits in genomic selection (GS) is a main objective for researchers and breeders in the livestock industry. This study aims at utilizing potentially functional SNPs and QTLs identified with various genome-wide association study (GWAS) models in GS of pig growth traits. We used three well-established GWAS methods, including the mixed linear model, Bayesian model and meta-analysis, as well as 60K SNP-chip and whole genome sequence (WGS) data from 1734 Yorkshire and 1123 Landrace pigs to detect SNPs related to four growth traits: average daily gain, backfat thickness, body weight and birth weight. A total of 1485 significant loci and 24 candidate genes which are involved in skeletal muscle development, fatty deposition, lipid metabolism and insulin resistance were identified. Compared with using all SNP-chip data, GS with the pre-selected functional SNPs in the standard genomic best linear unbiased prediction (GBLUP), and a two-kernel based GBLUP model yielded average gains in accuracy by 4 to 46% (from 0.19 ± 0.07 to 0.56 ± 0.07) and 5 to 27% (from 0.16 ± 0.06 to 0.57 ± 0.05) for the four traits, respectively, suggesting that the prioritization of preselected functional markers in GS models had the potential to improve prediction accuracies for certain traits in livestock breeding.

Association of IL-4 and IL-4R Polymorphisms with Litter Size Traits in Pigs

Simple Summary

The IL-4 and IL-4R cytokine genes are responsible for immune response in the reproductive system and are related to embryonic implantation and fetal survival during pregnancy in females. However, to date, their effects on litter size traits in pigs have been not elucidated. Therefore, the present study was conducted to verify the porcine IL-4 and IL-4R polymorphisms and assess how they affect litter size traits in commercial pigs. The findings suggested that the porcine IL-4 g.134993898T > C and IL-4R c.1577A > T polymorphisms were associated with the litter size traits. Therefore, the porcine IL-4 and IL-4R genes may be potentially effective genetic markers to improve the litter size traits in pigs.

Abstract

The interleukin-4 (IL-4) and interleukin-4 receptor (IL-4R) are cytokines that are involved in the immune and reproductive systems. This study aimed to verify the polymorphisms in the porcine IL-4 and IL-4R genes and to assess their effects on litter size traits in commercial pigs. Single nucleotide polymorphisms (SNPs) in the porcine IL-4 and IL-4R genes were genotyped by the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. A non-coding SNP of IL-4 g.134993898T > C and a non-synonymous SNP of IL-4R c.1577A > T (amino acid change at position 526, Q526L) were found to be segregating in Landrace sows. The IL-4 g.134993898T > C polymorphism was significantly associated with the number of piglets weaned alive (NWA) trait. The IL-4R c.1577A > T polymorphism was significantly associated with the number born alive (NBA) and NWA traits. Moreover, the accumulation of favorable alleles of these two SNP markers revealed significant associations with the NBA, NWA, and mean weight of piglets at weaning (MWW) traits. These findings indicate that the porcine IL-4 and IL-4R genes may contribute to the reproductive traits of pigs and could be used as candidate genes to improve litter size traits in the pig breeding industry.

Potential functional variants in AHR signaling pathways are associated with age at puberty in swine

Puberty in female pigs is defined as age at first estrus and gilts that have an earlier age at puberty are more likely to have greater lifetime productivity. Because age at puberty is predictive for sow longevity and lifetime productivity, but not routinely measured in commercial herds, it would be beneficial to use genomic or marker-assisted selection to improve these traits. A GWAS at the US Meat Animal Research Center (USMARC) identified several loci associated with age at puberty in pigs. Candidate genes in these regions were scanned for potential functional variants using sequence information from the USMARC swine population founder animals and public databases. In total, 135 variants (SNP and insertion/deletions) in 39 genes were genotyped in 1284 phenotyped animals from a validation population sired by Landrace and Yorkshire industry semen using the Agena MassArray system. Twelve variants in eight genes were associated with age at puberty (P < 0.005) with estimated additive SNP effects ranging from 1.6 to 5.3 days. Nine of these variants were non-synonymous coding changes in AHR, CYP1A2, OR2M4, SDCCAG8, TBC1D1 and ZNF608, two variants were deletions of one and four codons in aryl hydrocarbon receptor, AHR, and the most significant SNP was near an acceptor splice site in the acetyl-CoA carboxylase alpha, ACACA. Several of the loci identified have a physiological and a genetic role in sexual maturation in humans and other animals and are involved in AHR-mediated pathways. Further functional validation of these variants could identify causative mutations that influence age at puberty in gilts and possibly sow lifetime productivity.

Genomics of Heat Tolerance in Reproductive Performance Investigated in Four Independent Maternal Lines of Pigs

Improving swine climatic resilience through genomic selection has the potential to minimize welfare issues and increase the industry profitability. The main objective of this study was to investigate the genetic and genomic determinism of tolerance to heat stress in four independent purebred populations of swine. Three female reproductive traits were investigated: total number of piglets born (TNB), number of piglets born alive (NBA) and average birth weight (ABW). More than 80,000 phenotypic and 12,000 genotyped individuals were included in this study. Genomic random-regression models were fitted regressing the phenotypes of interest on a set of 95 environmental covariates extracted from public weather station records. The models yielded estimates of (genomic) reactions norms for individual pigs, as indicator of heat tolerance. Heat tolerance is a heritable trait, although the heritabilities are larger under comfortable than heat-stress conditions (larger than 0.05 vs. 0.02 for TNB; 0.10 vs. 0.05 for NBA; larger than 0.20 vs. 0.10 for ABW). TNB showed the lowest genetic correlation (-38%) between divergent climatic conditions, being the trait with the strongest impact of genotype by environment interaction, while NBA and ABW showed values slightly negative or equal to zero reporting a milder impact of the genotype by environment interaction. After estimating genetic parameters, a genome-wide association study was performed based on the single-step GBLUP method. Heat tolerance was observed to be a highly polygenic trait. Multiple and non-overlapping genomic regions were identified for each trait based on the genomic breeding values for reproductive performance under comfortable or heat stress conditions. Relevant regions were found on chromosomes (SSC) 1, 3, 5, 6, 9, 11, and 12, although there were important regions across all autosomal chromosomes. The genomic region located on SSC9 appears to be of particular interest since it was identified for two traits (TNB and NBA) and in two independent populations. Heat tolerance based on reproductive performance indicators is a heritable trait and genetic progress for heat tolerance can be achieved through genetic or genomic selection. Various genomic regions and candidate genes with important biological functions were identified, which will be of great value for future functional genomic studies.

Genetic and genomic characterization of vulva size traits in Yorkshire and Landrace gilts

BACKGROUND

Reproductive performance is critical for efficient swine production. Recent results indicated that vulva size (VS) may be predictive of reproductive performance in sows. Study objectives were to estimate genetic parameters, identify genomic regions associated, and estimate genomic prediction accuracies (GPA) for VS traits.

RESULTS

Heritability estimates of VS traits, vulva area (VA), height (VH), and width (VW) measurements, were moderately to highly heritable in Yorkshire, with 0.46 ± 0.10, 0.55 ± 0.10, 0.31 ± 0.09, respectively, whereas these estimates were low to moderate in Landrace, with 0.16 ± 0.09, 0.24 ± 0.11, and 0.08 ± 0.06, respectively. Genetic correlations within VS traits were very high for both breeds, with the lowest of 0.67 ± 0.29 for VH and VW for Landrace. Genome-wide association studies (GWAS) for Landrace, reveled genomic region associated with VS traits on Sus scrofa chromosome (SSC) 2 (154-157 Mb), 7 (107-110 Mb), 8 (4-6 Mb), and 10 (8-19 Mb). For Yorkshire, genomic regions on SSC 1 (87-91 and 282-287 Mb) and 5 (67 Mb) were identified. All regions explained at least 3.4% of the genetic variance. Accuracies of genomic prediction were moderate in Landrace, ranging from 0.30 (VH) to 0.61 (VA), and lower for Yorkshire, with 0.07 (VW) to 0.11 (VH). Between-breed and multi-breed genomic prediction accuracies were low.

CONCLUSIONS

Our findings suggest that VS traits are heritable in Landrace and Yorkshire gilts. Genomic analyses show that major QTL control these traits, and they differ between breed. Genomic information can be used to increase genetic gains for these traits in gilts. Additional research must be done to validate the GWAS and genomic prediction results reported in our study.

Porcine LIF gene polymorphisms and their association with litter size traits in four pig breeds

Leukemia inhibitory factor (LIF) is an important productivity-related gene in pigs. We found two polymorphisms — g.6646C>T and g.6988C>T — in exon 3 of LIF in pigs by using DNA sequencing and polymerase chain reaction-restriction fragment length polymorphism. Three genotypes were obtained and associated with litter size traits in Anqing Six-end-white (AQ), Wei (W), Wannan Black (WNB), and Large White (LW) pigs. At locus g.6646C>T, the g.6646C allele frequency variation was 0.6869 (AQ), 0.7473 (W), 1 (WNB), and 0.6852 (LW). In AQ pigs, sows with the TT genotype had higher total number of piglets born (TNB) and number of piglets born alive (NBA) in the first parity and multiparities (P < 0.01). In W and LW pigs, sows with the CC genotype had higher TNB and NBA in multiparities (P < 0.01). At locus g.6988C>T, the g.6988C allele frequency variation was 1 (AQ), 0.6154 (W), 1 (WNB), and 0.6667 (LW). The CC genotype significantly differed from CT or TT genotypes (P < 0.01) for TNB and NBA in W and LW pigs. Thus, LIF was shown to have a significant influence on litter size. Therefore, g.6646C>T and g.6988C>T loci of LIF could be potential marker-assisted selection tools for improving litter size in pig production.

A genomewide association study in divergently selected lines in rabbits reveals novel genomic regions associated with litter size traits

Uterine capacity (UC), defined as the total number of kits from unilaterally ovariectomized does at birth, has a high genetic correlation with litter size. The aim of our research was to identify genomic regions associated with litter size traits through a genomewide association study using rabbits from a divergent selection experiment for UC. A high-density SNP array (200K) was used to genotype 181 does from a control population, high and low UC lines. Traits included total number born (TNB), number born alive (NBA), number born dead, ovulation rate (OR), implanted embryos (IE) and embryo, foetal and prenatal survivals at second parity. We implemented the Bayes B method and the associations were tested by Bayes factors and the percentage of genomic variance (GV) explained by windows. Different genomic regions associated with TNB, NBA, IE and OR were found. These regions explained 7.36%, 1.27%, 15.87% and 3.95% of GV, respectively. Two consecutive windows on chromosome 17 were associated with TNB, NBA and IE. This genomic region accounted for 6.32% of GV of TNB. In this region, we found the BMP4, PTDGR, PTGER2, STYX and CDKN3 candidate genes which presented functional annotations linked to some reproductive processes. Our findings suggest that a genomic region on chromosome 17 has an important effect on litter size traits. However, further analyses are needed to validate this region in other maternal rabbit lines.

A Survey of Copy Number Variation in the Porcine Genome Detected From Whole-Genome Sequence

Copy number variations (CNVs) are gains and losses of large regions of genomic sequence between individuals of a species. Although CNVs have been associated with various phenotypic traits in humans and other species, the extent to which CNVs impact phenotypic variation remains unclear. In swine, as well as many other species, relatively little is understood about the frequency of CNV in the genome, sizes, locations, and other chromosomal properties. In this work, we identified and characterized CNV by utilizing whole-genome sequence from 240 members of an intensely phenotyped experimental swine herd at the U.S. Meat Animal Research Center (USMARC). These animals included all 24 of the purebred founding boars (12 Duroc and 12 Landrace), 48 of the founding Yorkshire-Landrace composite sows, 109 composite animals from generations 4 through 9, 29 composite animals from generation 15, and 30 purebred industry boars (15 Landrace and 15 Yorkshire) used as sires in generations 10 through 15. Using a combination of split reads, paired-end mapping, and read depth approaches, we identified a total of 3,538 copy number variable regions (CNVRs), including 1,820 novel CNVRs not reported in previous studies. The CNVRs covered 0.94% of the porcine genome and overlapped 1,401 genes. Gene ontology analysis identified that CNV-overlapped genes were enriched for functions related to organism development. Additionally, CNVRs overlapped with many known quantitative trait loci (QTL). In particular, analysis of QTL previously identified in the USMARC herd showed that CNVRs were most overlapped with reproductive traits, such as age of puberty and ovulation rate, and CNVRs were significantly enriched for reproductive QTL.

Effect of single-nucleotide polymorphisms on specific reproduction parameters in Hungarian Large White sows

The aim of this study was to reveal the effect of single-nucleotide polymorphisms (SNPs) on the total number of piglets born (TNB), the litter weight born alive (LWA), the number of piglets born dead (NBD), the average litter weight on the 21st day (M21D) and the interval between litters (IBL). Genotypes were determined on a high-density Illumina Porcine SNP 60K BeadChip. Data screening and data identification were performed by a multi-locus mixed-model. Statistical analyses were carried out to find associations between individual genotypes of 290 Hungarian Large White sows and the investigated reproduction parameters. According to the analysis outcome, three SNPs were identified to be associated with TNB. These loci are located on chromosomes 1, 6 and 13 (-log₁₀P = 6.0, 7.86 and 6.22, the frequencies of their minor alleles, MAF, were 0.298, 0.299 and 0.364, respectively). Two loci showed considerable association (-log₁₀P = 10.35 and 10.46) with LWA on chromosomes 5 and X, the MAF were 0.425 and 0.446, respectively. Seven loci were found to be associated with NBD. These loci are located on chromosomes 5, 6, 13, 14, 15, 16 and 18 (-log₁₀P = 10.95, 5.43, 8.29, 6.72, 6.81, 5.90, and 5.15, respectively). One locus showed association (-log₁₀P = 5.62) with M21D on chromosome 1 (the MAF was 0.461). Another locus was found to be associated with IBL on chromosome 8 (-log₁₀P = 7.56; the MAF was 0.438). The above-mentioned loci provide a straightforward possibility to assist selection by molecular tools and, consequently, to improve the competitiveness of the Hungarian Large White (HLW) breed.

Invited review: Piglet survival: benefits of the immunocompetence

Piglet mortality has a negative impact on animal welfare and public acceptance. Moreover, the number of weaned piglets per sow mainly determines the profitability of piglet production. Increased litter sizes are associated with lower birth weights and piglet survival. Decreased survival rates and performance of piglets make the control of diseases and infections within pig production even more crucial. Consequently, selection for immunocompetence becomes an important key aspect within modern breeding programmes. However, the phenotypic recording of immune traits is difficult and expensive to realize within farm routines. Even though immune traits show genetic variability, only few examples exist on their respective suitability within a breeding programme and their relationships to economically important production traits. The analysis of immune traits for an evaluation of immunocompetence to gain a generally improved immune response is promising. Generally, in-depth knowledge of the genetic background of the immune system is needed to gain helpful insights about its possible incorporation into breeding programmes. Possible physiological drawbacks for enhanced immunocompetence must be considered with regards to the allocation theory and possible trade-offs between the immune system and performance. This review aims to discuss the relationships between the immunocompetence of the pig, piglet survival as well as the potential of these traits to be included into a breeding strategy for improved robustness.

Whole-genome sequencing identifies potential candidate genes for reproductive traits in pigs

Reproductive performance is a complex quantitative trait, that is determined by multiple genes, regulatory pathways and environmental factors. A list of major genes with large effect have been detected, although multiple QTLs are identified. To identify candidate genes for pig prolificacy, whole genome variants from five high- and five low-prolificacy Yorkshire sows were collected using whole-genome resequencing. A total of 13,955,609 SNPs and 2,666,366 indels were detected across the genome. Common differential SNPs and indels were identified between the two groups of sows. Genes encoding components of the TGF-beta signaling pathway were enriched with the variations, including BMP5, BMP6, BMP7, ACVR1, INHBA, ZFYVE9, TGFBR2, DCN, ID4, BAMBI, and ACVR2A. Several differential variants within these genes related to reproductive traits were identified to be associated with litter size. A comparison of selective regions and published QTL data suggests that NEDD9, SLC39A11, SNCA, and UNC5D are candidate genes for reproduction traits.

Using genome wide association studies to identify common QTL regions in three different genetic backgrounds based on Iberian pig breed

One of the major limitation for the application of QTL results in pig breeding and QTN identification has been the limited number of QTL effects validated in different animal material. The aim of the current work was to validate QTL regions through joint and specific genome wide association and haplotype analyses for growth, fatness and premier cut weights in three different genetic backgrounds, backcrosses based on Iberian pigs, which has a major role in the analysis due to its high productive relevance. The results revealed nine common QTL regions, three segregating in all three backcrosses on SSC1, 0–3 Mb, for body weight, on SSC2, 3–9 Mb, for loin bone-in weight, and on SSC7, 3 Mb, for shoulder weight, and six segregating in two of the three backcrosses, on SSC2, SSC4, SSC6 and SSC10 for backfat thickness, shoulder and ham weights. Besides, 18 QTL regions were specifically identified in one of the three backcrosses, five identified only in BC_LD, seven in BC_DU and six in BC_PI. Beyond identifying and validating QTL, candidate genes and gene variants within the most interesting regions have been explored using functional annotation, gene expression data and SNP identification from RNA-Seq data. The results allowed us to propose a promising list of candidate mutations, those identified in PDE10A, DHCR7, MFN2 and CCNY genes located within the common QTL regions and those identified near ssc-mir-103-1 considered PANK3 regulators to be further analysed.

Genome-wide association study for reproductive traits in a Large White pig population

Using the PorcineSNP80 BeadChip, we performed a genome‐wide association study for seven reproductive traits, including total number born, number born alive, litter birth weight, average birth weight, gestation length, age at first service and age at first farrowing, in a population of 1207 Large White pigs. In total, we detected 12 genome‐wide significant and 41 suggestive significant SNPs associated with six reproductive traits. The proportion of phenotypic variance explained by all significant SNPs for each trait ranged from 4.46% (number born alive) to 11.49% (gestation length). Among them, 29 significant SNPs were located within known QTL regions for swine reproductive traits, such as corpus luteum number, stillborn number and litter size, of which one QTL region associated with litter size contained the ALGA0098819 SNP for total number born. Subsequently, we found that 376 functional genes contained or were near these significant SNPs. Of these, 14 genes—BHLHA15, OCM2, IL1B2, GCK, SMAD2, HABP2, PAQR5, GRB10, PRELID2, DMKN, GPI, GPIHBP1, ADCY2 and ACVR2B—were considered important candidates for swine reproductive traits based on their critical roles in embryonic development, energy metabolism and growth development. Our findings contribute to the understanding of the genetic mechanisms for reproductive traits and could have a positive effect on pig breeding programs.

Genome-Wide Association Study of Piglet Uniformity and Farrowing Interval

Piglet uniformity (PU) and farrowing interval (FI) are important reproductive traits related to production and economic profits in the pig industry. However, the genetic architecture of the longitudinal trends of reproductive traits still remains elusive. Herein, we performed a genome-wide association study (GWAS) to detect potential genetic variation and candidate genes underlying the phenotypic records at different parities for PU and FI in a population of 884 Large White pigs. In total, 12 significant SNPs were detected on SSC1, 3, 4, 9, and 14, which collectively explained 1–1.79% of the phenotypic variance for PU from parity 1 to 4, and 2.58–4.11% for FI at different stages. Of these, seven SNPs were located within 16 QTL regions related to swine reproductive traits. One QTL region was associated with birth body weight (related to PU) and contained the peak SNP MARC0040730, and another was associated with plasma FSH concentration (related to FI) and contained the SNP MARC0031325. Finally, some positional candidate genes for PU and FI were identified because of their roles in prenatal skeletal muscle development, fetal energy substrate, pre-implantation, and the expression of mammary gland epithelium. Identification of novel variants and candidate genes will greatly advance our understanding of the genetic mechanisms of PU and FI, and suggest a specific opportunity for improving marker assisted selection or genomic selection in pigs.

Genetic Architecture and Candidate Genes Identified for Follicle Number in Chicken

Follicular development has a major impact on reproductive performance. Most previous researchers focused on molecular mechanisms of follicular development. The genetic architecture underlying the number of follicle, however, has yet not to be thoroughly defined in chicken. Here we report a genome-wide association study for the genetic architecture determining the numbers of follicles in a large F₂ resource population. The results showed heritability were low to moderate (0.05-0.28) for number of pre-ovulatory follicles (POF), small yellow follicles (SYF) and atresia follicles (AF). The highly significant SNPs associated with SYF were mainly located on GGA17 and GGA28. Only four significant SNPs were identified for POF on GGA1. The variance partitioned across chromosomes and chromosome lengths had a linear relationship for SYF (R² = 0.58). The enriched genes created by the closest correspondent significant SNPs were found to be involved in biological pathways related to cell proliferation, cell cycle and cell survival. Two promising candidate genes, AMH and RGS3, were suggested to be prognostic biomarkers for SYF. In conclusion, this study offers the first evidence of genetic variance and positional candidate genes which influence the number of SYF in chicken. These identified informative SNPs may facilitate selection for an improved reproductive performance of laying hens.

A systematic survey to identify lethal recessive variation in highly managed pig populations

BACKGROUND

Lethal recessive variation can cause prenatal death of homozygous offspring. Although usually present at low-frequency in populations, the impact on individual fitness can be substantial. Until recently, the presence of recessive embryonic lethal variation could only be measured indirectly through reduced fertility. In this study, we estimate the presence of genetic loci associated with both early and late termination of development during gestation in pigs from the wealth of genome data routinely generated by a commercial breeding company.

RESULTS

We examined three commercial pig (Sus scrofa) populations for potentially deleterious genetic variation based on 80 K SNP-chip genotypes, and estimate the effects on reproductive traits. 24,000 pigs from three populations were analyzed for missing or depletion of homozygous haplotypes. We identified 145 haplotypes (ranging from 0.5-4 Mb in size) in the genome with complete absence or depletion of homozygous animals. Thirty-five haplotypes show a negative effect on at least one of the analysed reproductive traits (total number born, number of stillborn, and number of mummified piglets). One variant in particular appeared to result in relative late termination of development of fetuses, responsible for a significant fraction of observed stillborn piglets ('mummies'), as they die mid-gestation. Moreover, we identified the BMPER gene as a likely candidate underlying this phenomenon.

CONCLUSIONS

Our study shows that although lethal recessive variation is present, the frequency of these alleles is invariably low in these highly managed populations. Nevertheless, due to cumulative effects of deleterious variants, large numbers of affected offspring are produced. Furthermore, our study demonstrates the use of a large-scale commercial genetic experiment to systematically screen for 'natural knockouts' that can increase understanding of gene function.

Use of haplotypes to identify regions harbouring lethal recessive variants in pigs

Background

Lethal recessive genetic variants are maintained at relatively low frequencies in a population in the heterozygous state, but by definition are fatal and therefore unobserved in the homozygous state. Since haplotypes allow the tagging of rare and untyped genetic variants, they have potential for studying lethal recessive variants. In this study, we used a large commercial population to identify putative lethal recessive haplotypes that impact either the total number born (TNB) or the number born alive (NBA) as a proportion of the total number born (NBA/TNB). We also compared the use of haplotypes with a single nucleotide polymorphism (SNP)-by-SNP approach and examined the benefits of using additional haplotypes imputed from low-density genotype data for the detection of lethal recessive variants. Candidate haplotypes were identified using population-wide haplotype frequencies and within-family analyses. These candidate haplotypes were subsequently assessed for putative lethal recessive effects on TNB and NBA/TNB by comparing carrier-to-carrier matings with carrier-to-non-carrier matings.

Results

Using both medium-density and imputed low-density genotype data six regions were identified as containing putative lethal recessive haplotypes that had an effect on TNB. It is likely that these regions were related to at least four putative lethal recessive variants, each located on a different chromosome. Evidence for putative lethal recessive effects on TNB was found on chromosomes 1, 6, 10 and 14 using haplotypes. Using haplotypes from individuals genotyped only at medium-density or a SNP-by-SNP approach did not detect any lethal recessive effects. No lethal recessive haplotypes or SNPs were detected that had an effect on NBA/TNB.

Conclusions

We show that the use of haplotypes from combining medium-density and imputed low-density genotype data is superior for the identification of lethal recessive variants compared to both a SNP-by-SNP approach and to the use of only medium-density data. We developed a formal statistical framework that provided sufficient power to detect lethal recessive variants in species, which produce large full-sib families, while reducing false positive or type I errors. Applying this framework results in improvements in reproductive performance by purging lethal recessive alleles from a population in a timely and cost-effective manner.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-017-0332-3) contains supplementary material, which is available to authorized users.

Genomewide association of piglet responses to infection with one of two porcine reproductive and respiratory syndrome virus isolates

Porcine reproductive and respiratory syndrome (PRRS) is a devastating disease in the swine industry. Identification of host genetic factors that enable selection for improved performance during PRRS virus (PRRSV) infection would reduce the impact of this disease on animal welfare and production efficiency. We conducted genomewide association study (GWAS) analyses of data from 13 trials of approximately 200 commercial crossbred nursery-age piglets that were experimentally infected with 1 of 2 type 2 isolates of PRRSV (NVSL 97-7985 [NVSL] and KS2006-72109 [KS06]). Phenotypes analyzed were viral load (VL) in blood during the first 21 d after infection (dpi) and weight gain (WG) from 0 to 42 dpi. We accounted for the previously identified QTL in the region on SSC4 in our models to increase power to identify additional regions. Many regions identified by single-SNP analyses were not identified using Bayes-B, but both analyses identified the same regions on SSC3 and SSC5 to be associated with VL in the KS06 trials and on SSC6 in the NVSL trials ( < 5 × 10); for WG, regions on SSC5 and SSC17 were associated in the NVSL trials ( < 3 × 10). No regions were identified with either method for WG in the KS06 trials. Except for the region on SSC4, which was associated with VL for both isolates (but only with WG for NVSL), identified regions did not overlap between the 2 PRRSV isolate data sets, despite high estimates of the genetic correlation between isolates for traits based on these data. We also identified genomic regions whose associations with VL or WG interacted with either PRRSV isolate or with genotype at the SSC4 QTL. Gene ontology (GO) annotation terms for genes located near moderately associated SNP ( < 0.003) were enriched for multiple immunologically (VL) and metabolism- (WG) related GO terms. The biological relevance of these regions suggests that, although it may increase the number of false positives, the use of single-SNP analyses and a relaxed threshold also increased the identification of true positives. In conclusion, although only the SSC4 QTL was associated with response to both PRRSV isolates, genes near associated SNP were enriched for the same GO terms across PRRSV isolates, suggesting that host responses to these 2 isolates are affected by the actions of many genes that function together in similar biological processes.

Novel effects of identified SNPs within the porcine Pregnancy-Associated Glycoprotein gene family (pPAGs) on the major reproductive traits in Hirschmann hybrid-line sows

This is the first study describing identification of SNPs within the multiple and polymorphic Pregnancy-Associated Glycoprotein gene family (PAGs) in the genome of the domestic pig (pPAGs). We identified pPAG-like (pPAG-L) genotypes in primiparous and multiparous farmed hybrid-line JSR Hirschmann (Hrn) sows (N=159), in which various novel associations with their phenotypes for the major reproductive traits have been discovered. Genomic DNA templates were isolated from the blood and different pPAG-L primers were used to amplify various regions by PCR. Electrophoretically-separated amplicons were selected, purified and sequenced. All identified SNPs were verified for possible pPAG2-L genotype associations with the major reproductive traits. In total, 196 SNPs were identified within the entire structure of the pPAG2-Ls, encompassing 9 exons and 8 (A-H) introns, resembling all aspartic proteinases. It was discovered that among all SNPs, one diplotype localized in exon 6 (657C>T/749G>C; pPAG2 ORF cDNA numbering; L34361) caused amino acid substitutions (Asp²²⁰→Asn and Ser²⁵⁰→Thr) in the polypeptide precursors and was associated with an increase in the number of live-born piglets (P≤0.05) in Hrn sows. In turn, co-localized SNP (504g>a; KF537535 numbering) in the intron F of the pPAG2-Ls, but only in the homozygotic genotype (gg), was associated with an increased number of live-born (P≤0.01) and weaned (P≤0.05) piglets in the Hrn sows. These results qualify the pPAG2-Ls as candidate genes of the main QTLs. The novel pPAG SNP profiles provide the basis for a diagnostic genotyping test required for early pre-selection of female/male piglets, presumably mainly useful in various breeding herds.

Genomewide association analysis for average birth interval and stillbirth in swine

Reproductive efficiency has a great impact on the economic success of pork production. Stillborn pigs and average birth interval contribute to the number of pigs born alive in a litter. To better understand the underlying genetics of these traits, a genomewide association study was undertaken. Samples of DNA were collected and tested using the Illumina Porcine SNP60 BeadChip from 798 females farrowing over a 4-yr period (all first parity). Birth intervals and piglet birth status (stillborn or alive) were determined by videotaping each farrowing event. A total of 41,148 SNP were tested using the Bayes C option of GenSel (version 4.61) and 1-Mb windows. These 1-Mb windows explained proportions of 0.017, 0.002, 0.032, 0.029, and 0.030 of the total variation, respectively, for litter average birth interval after deletion of the last piglet born, last birth interval in the litter, number of stillborn piglets ignoring the last piglet born, number of stillborns in the last birth position, and percent stillborn ignoring the last piglet. Significant 1-Mb nonoverlapping SNP windows were identified by using a conservative approach requiring 1-Mb windows to have a genetic variance ≥1.0% of genomic variance and these were considered to be QTL. Quantitative trait loci were located for number of stillborn piglets ignoring the last piglet born (1 QTL), number of stillborns in the last birth position (1 QTL), and percent stillborn ignoring the last piglet (3 QTL). In addition, 2, 13, 3, and 6 suggestive 1-Mb nonoverlapping SNP windows were identified for litter average birth interval after deletion of the last piglet born, number of stillborn piglets ignoring the last piglet born, number of stillborns in the last birth position, and percent stillborn ignoring the last piglet, respectively. Possible candidate genes affecting both birth interval and stillbirth included () and (). Possible genes affecting only birth interval included (), and (), and those affecting only stillbirth included (), LOC100518697 (a nostrin-like gene), and (). The QTL and the suggestive 1-Mb nonoverlapping SNP windows may lead to genetic markers for marker assisted selection, marker assisted management, or genomic selection applications in commercial pig populations.

A study of vertebra number in pigs confirms the association of vertnin and reveals additional QTL

Background

Formation of the vertebral column is a critical developmental stage in mammals. The strict control of this process has resulted in little variation in number of vertebrae across mammalian species and no variation within most mammalian species. The pig is quite unique as considerable variation exists in number of thoracic vertebrae as well as number of lumbar vertebrae. At least two genes have been identified that affect number of vertebrae in pigs yet considerable genetic variation still exists. Therefore, a genome-wide association (GWA) analysis was conducted to identify additional genomic regions that affect this trait.

Results

A total of 1883 animals were phenotyped for the number of ribs and thoracolumbar vertebrae as well as successfully genotyped with the Illumina Porcine SNP60 BeadChip. After data editing, 41,148 SNP markers were included in the GWA analysis. These animals were also phenotyped for kyphosis. Fifty-three 1 Mb windows each explained at least 1.0 % of the genomic variation for vertebrae counts while 16 regions were significant for kyphosis. Vertnin genotype significantly affected vertebral counts as well. The region with the largest effect for number of lumbar vertebrae and thoracolumbar vertebrae were located over the Hox B gene cluster and the largest association for thoracic vertebrae number was over the Hox A gene cluster. Genetic markers in significant regions accounted for approximately 50 % of the genomic variation. Less genomic variation for kyphosis was described by QTL regions and no region was associated with kyphosis and vertebra counts.

Conclusions

The importance of the Hox gene families in vertebral development was highlighted as significant associations were detected over the A, B and C families. Further evaluation of these regions and characterization of variants within these genes will expand our knowledge on vertebral development using natural genetic variants segregating in commercial swine.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0286-9) contains supplementary material, which is available to authorized users.

Stories and Challenges of Genome Wide Association Studies in Livestock - A Review

Undoubtedly livestock is one of the major contributors to the economy of any country. The economic value of livestock includes meat, dairy products, fiber, fertilizer etc. Understanding and identifying the associations of quantitative trait loci (QTL) with the economically important traits is believed to substantially benefit the livestock industry. The past two decades have seen a flurry of interest in mapping the QTL associated with traits of economic importance on the genome. With the availability of single nucleotide polymorphism chip of various densities it is possible to identify regions, QTL and genes on the genome that explain the association and its effect on the phenotype under consideration. Remarkable advancement has been seen in genome wide association studies (GWAS) since its inception till the present day. In this review we describe the progress and challenges of GWAS in various livestock species.

Genome-wide association study of body weight in Australian Merino sheep reveals an orthologous region on OAR6 to human and bovine genomic regions affecting height and weight

Background

Body weight (BW) is an important trait for meat production in sheep. Although over the past few years, numerous quantitative trait loci (QTL) have been detected for production traits in cattle, few QTL studies have been reported for sheep, with even fewer on meat production traits. Our objective was to perform a genome-wide association study (GWAS) with the medium-density Illumina Ovine SNP50 BeadChip to identify genomic regions and corresponding haplotypes associated with BW in Australian Merino sheep.

Methods

A total of 1781 Australian Merino sheep were genotyped using the medium-density Illumina Ovine SNP50 BeadChip. Among the 53 862 single nucleotide polymorphisms (SNPs) on this array, 48 640 were used to perform a GWAS using a linear mixed model approach. Genotypes were phased with hsphase; to estimate SNP haplotype effects, linkage disequilibrium blocks were identified in the detected QTL region.

Results

Thirty-nine SNPs were associated with BW at a Bonferroni-corrected genome-wide significance threshold of 1 %. One region on sheep (Ovis aries) chromosome 6 (OAR6) between 36.15 and 38.56 Mb, included 13 significant SNPs that were associated with BW; the most significant SNP was OAR6_41936490.1 (P = 2.37 × 10⁻¹⁶) at 37.69 Mb with an allele substitution effect of 2.12 kg, which corresponds to 0.248 phenotypic standard deviations for BW. The region that surrounds this association signal on OAR6 contains three genes: leucine aminopeptidase 3 (LAP3), which is involved in the processing of the oxytocin precursor; NCAPG non-SMC condensin I complex, subunit G (NCAPG), which is associated with foetal growth and carcass size in cattle; and ligand dependent nuclear receptor corepressor-like (LCORL), which is associated with height in humans and cattle.

Conclusions

The GWAS analysis detected 39 SNPs associated with BW in sheep and a major QTL region was identified on OAR6. In several other mammalian species, regions that are syntenic with this region have been found to be associated with body size traits, which may reflect that the underlying biological mechanisms share a common ancestry. These findings should facilitate the discovery of causative variants for BW and contribute to marker-assisted selection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0142-4) contains supplementary material, which is available to authorized users.

A Major Gene for Bovine Ovulation Rate

Half-sib daughters sired by a bull believed to be a carrier of a major gene for high ovulation rate were evaluated for ovulation rate and genotyped in an effort to both test the hypothesis of segregation of a major gene and to map the gene’s location. A total of 131 daughters were produced over four consecutive years at a University of Wisconsin-Madison research farm. All were evaluated for ovulation rate over an average of four estrous cycles using transrectal ultrasonography. The sire and all daughters were genotyped using a 3K SNP chip and the genotype and phenotype data were used in a linkage analysis. Subsequently, daughters recombinant within the QTL region and the sire were genotyped successively with 50K and 777K SNP chips to refine the location of the causative polymorphism. Positional candidate genes within the fine-mapped region were examined for polymorphism by Sanger sequencing of PCR amplicons encompassing coding and 5’ and 3’ flanking regions of the genes. Sire DNA was used as template in the PCR reactions. Strong evidence of a major gene for ovulation rate was observed (p<1x10^-28) with the gene localized to bovine chromosome 10. Fine-mapping subsequently reduced the location to a 1.2 Mb region between 13.6 and 14.8 Mb on chromosome 10. The location identified does not correspond to that for any previously identified major gene for ovulation rate. This region contains three candidate genes, SMAD3, SMAD6 and IQCH. While candidate gene screening failed to identify the causative polymorphism, three polymorphisms were identified that can be used as a haplotype to track inheritance of the high ovulation rate allele in descendants of the carrier sire.

Genome-wide association studies (GWAS) identify a QTL close to PRKAG3 affecting meat pH and colour in crossbred commercial pigs

Background

Improving meat quality is a high priority for the pork industry to satisfy consumers’ preferences. GWAS have become a state-of-the-art approach to genetically improve economically important traits. However, GWAS focused on pork quality are still relatively rare.

Results

Six genomic regions were shown to affect loin pH and Minolta colour a* and b* on both loin and ham through GWAS in 1943 crossbred commercial pigs. Five of them, located on Sus scrofa chromosome (SSC) 1, SSC5, SSC9, SSC16 and SSCX, were associated with meat colour. However, the most promising region was detected on SSC15 spanning 133–134 Mb which explained 3.51% - 17.06% of genetic variance for five measurements of pH and colour. Three SNPs (ASGA0070625, MARC0083357 and MARC0039273) in very strong LD were considered most likely to account for the effects in this region. ASGA0070625 is located in intron 2 of ZNF142, and the other two markers are close to PRKAG3, STK36, TTLL7 and CDK5R2. After fitting MARC0083357 (the closest SNP to PRKAG3) as a fixed factor, six SNPs still remained significant for at least one trait. Four of them are intragenic with ARPC2, TMBIM1, NRAMP1 and VIL1, while the remaining two are close to RUFY4 and CDK5R2. The gene network constructed demonstrated strong connections of these genes with two major hubs of PRKAG3 and UBC in the super-pathways of cell-to-cell signaling and interaction, cellular function and maintenance. All these pathways play important roles in maintaining the integral architecture and functionality of muscle cells facing the dramatic changes that occur after exsanguination, which is in agreement with the GWAS results found in this study.

Conclusions

There may be other markers and/or genes in this region besides PRKAG3 that have an important effect on pH and colour. The potential markers and their interactions with PRKAG3 require further investigation

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0192-1) contains supplementary material, which is available to authorized users.

Evaluation of quantitative trait loci affecting intramuscular fat and reproductive traits in pigs using marker-assisted introgression

We investigated the effects of previously identified quantitative trait loci (QTL) in an experimental backcross (BC) between Chinese Meishan pigs and commercial Duroc pigs. We performed marker-assisted introgression of two QTL for intramuscular fat (IMF) content (IMF population) and three QTL for reproductive traits (reproduction population) from a donor Meishan pig into a recipient Duroc pig. At the fourth BC generation of the IMF population and third BC generation of the reproduction population, carrier animals were selected for the production of animals homozygous for the QTL. Our previous studies have shown that the presence of a Meishan allele on the IMF QTL is associated with low IMF values, and the Meishan allele on the reproductive QTL is associated with large litters. In this study, the presence of a Duroc allele at the IMF QTL on SSC9 resulted in a 0.27% increase in IMF (additive effect = 0.27 ± 0.08), whereas the presence of a Meishan allele at the IMF QTL on SSC7 resulted in a 0.34% increase in IMF (additive effect = -0.34 ± 0.09). The presence of the Meishan allele at the IMF QTL on SSC7 thus had the opposite effect to our previous studies, that is, increased IMF. In the reproduction population, we observed no differences between the genotypes of the three QTL in regard to number of corpora lutea or litter size. Marker-assisted introgression at these QTL is thus unlikely to result in an associated increase in litter size. These results show that it is possible to introgress alleles from other breeds into a selection population using molecular markers; any unexpected results might be associated with the genetic background.

Genomewide association and identification of candidate genes for ovulation rate in swine

Reproductive efficiency has a great impact on the economic success of pork production. Ovulation rate is an early component of reproduction efficiency and contributes to the number of pigs born in a litter. To better understand the underlying genetics of ovulation rate, a genomewide association study was undertaken. Samples of DNA were collected and tested using the Illumina Porcine SNP60 BeadChip from 1,180 females with ovulation measurements ranging from never farrowed to measurements taken after parity 2. A total of 41,848 SNP were tested using the Bayes C option of GenSel. After the Bayes C analysis, SNP were assigned to sliding windows of 5 consecutive SNP by chromosome-position order beginning with the first 5 SNP on SSC1 and ending with the last 5 SNP on SSCX. The 5-SNP windows were analyzed using the Predict option of GenSel. From the Predict analysis, putative QTL were selected having no overlap with other 5-SNP window groups, no overlap across chromosomes, and the highest genetic variation. These putative QTL were submitted to statistical testing using the bootstrap option of GenSel. Of the putative QTL tested, 80 were found to be statistically significant (P < 0.01). Ten QTL were found on SSC1, 12 on SSC2, 4 on SSC3, 8 on SSC4, 3 on SSC5, 3 on SSC6, 3 on SSC7, 4 on SSC8, 2 on SSC9, 4 on SSC10, 1 on SSC12, 4 on SSC13, 2 on SSC14, 4 on SSC15, 4 on SSC16, 6 on SSC17, 4 on SSC18, and 1 on SSCX. Sixteen QTL were found to be statistically significant at the P < 0.001 level. Six additional QTL were significant at the P = 0.001 level. These 22 QTL accounted for 71.10% of the total genetic variance. The most compelling candidate genes in these regions include Estrogen receptor 1, growth differentiation factor 9, and inhibin βA. These QTL, when combined with information on genes found in the same regions, should provide useful information that could be used for marker assisted selection, marker assisted management, or genomic selection applications in commercial pig populations.

Genetic analysis of reproductive traits and antibody response in a PRRS outbreak herd

Porcine reproductive and respiratory syndrome (PRRS) is the most economically significant disease impacting pig production in North America, Europe, and Asia, causing reproductive losses such as increased rates of stillbirth and mummified piglets. The objective of this study was to explore the genetic basis of host response to the PRRS virus (PRRSV) in a commercial multiplier sow herd before and after a PRRS outbreak, using antibody response and reproductive traits. Reproductive data comprising number born alive (NBA), number alive at 24 h (NA24), number stillborn (NSB), number born mummified (NBM), proportion born dead (PBD), number born dead (NBD), number weaned (NW), and number of mortalities through weaning (MW) of 5,227 litters from 1,967 purebred Landrace sows were used along with a pedigree comprising 2,995 pigs. The PRRS outbreak date was estimated from rolling averages of farrowing traits and was used to split the data into a pre-PRRS phase and a PRRS phase. All 641 sows in the herd during the outbreak were blood sampled 46 d after the estimated outbreak date and were tested for anti-PRRSV IgG using ELISA (sample-to-positive [S/P] ratio). Genetic parameters of traits were estimated separately for the pre-PRRS and PRRS phase data sets. Sows were genotyped using the PorcineSNP60 BeadChip, and genome-wide association studies (GWAS) were performed using method Bayes B. Heritability estimates for reproductive traits ranged from 0.01 (NBM) to 0.12 (NSB) and from 0.01 (MW) to 0.12 (NBD) for the pre-PRRS and PRRS phases, respectively. S/P ratio had heritability (0.45) and strong genetic correlations with most traits, ranging from -0.72 (NBM) to 0.73 (NBA). In the pre-PRRS phase, regions associated with NSB and PBD explained 1.6% and 3% of the genetic variance, respectively. In the PRRS phase, regions associated with NBD, NSB, and S/P ratio explained 0.8%, 11%, and 50.6% of the genetic variance, respectively. For S/P ratio, 2 regions on SSC 7 (SSC7) separated by 100 Mb explained 40% of the genetic variation, including a region encompassing the major histocompatibility complex, which explained 25% of the genetic variance. These results indicate a significant genomic component associated with PRRSV antibody response and NSB in this data set. Also, the high heritability and genetic correlation estimates for S/P ratio during the PRRS phase suggest that S/P ratio could be used as an indicator of the impact of PRRS on reproductive traits.

Genome-wide association analysis and genomic prediction of Mycobacterium avium subspecies paratuberculosis infection in US Jersey cattle

Paratuberculosis (Johne’s disease), an enteric disorder in ruminants caused by Mycobacterium avium subspecies paratuberculosis (MAP), causes economic losses in excess of $200 million annually to the US dairy industry. To identify genomic regions underlying susceptibility to MAP infection in Jersey cattle, a case-control genome-wide association study (GWAS) was performed. Blood and fecal samples were collected from ∼5,000 mature cows in 30 commercial Jersey herds from across the US. Discovery data consisted of 450 cases and 439 controls genotyped with the Illumina BovineSNP50 BeadChip. Cases were animals with positive ELISA and fecal culture (FC) results. Controls were animals negative to both ELISA and FC tests that matched cases on birth date and herd. Validation data consisted of 180 animals including 90 cases (positive to FC) and 90 controls (negative to ELISA and FC), selected from discovery herds and genotyped by Illumina BovineLD BeadChip (∼7K SNPs). Two analytical approaches were used: single-marker GWAS using the GRAMMAR-GC method and Bayesian variable selection (Bayes C) using GenSel software. GRAMMAR-GC identified one SNP on BTA7 at 68 megabases (Mb) surpassing a significance threshold of 5×10⁻⁵. ARS-BFGL-NGS-11887 on BTA23 (27.7 Mb) accounted for the highest percentage of genetic variance (3.3%) in the Bayes C analysis. SNPs identified in common by GRAMMAR-GC and Bayes C in both discovery and combined data were mapped to BTA23 (27, 29 and 44 Mb), 3 (100, 101, 106 and 107 Mb) and 17 (57 Mb). Correspondence between results of GRAMMAR-GC and Bayes C was high (70–80% of most significant SNPs in common). These SNPs could potentially be associated with causal variants underlying susceptibility to MAP infection in Jersey cattle. Predictive performance of the model developed by Bayes C for prediction of infection status of animals in validation set was low (55% probability of correct ranking of paired case and control samples).

Genome wide screening of candidate genes for improving piglet birth weight using high and low estimated breeding value populations

Birth weight is an economically important trait in pig production because it directly impacts piglet growth and survival rate. In the present study, we performed a genome wide survey of candidate genes and pathways associated with individual birth weight (IBW) using the Illumina PorcineSNP60 BeadChip on 24 high (HEBV) and 24 low estimated breeding value (LEBV) animals. These animals were selected from a reference population of 522 individuals produced by three sires and six dam lines, which were crossbreds with multiple breeds. After quality-control, 43,257 SNPs (single nucleotide polymorphisms), including 42,243 autosomal SNPs and 1,014 SNPs on chromosome X, were used in the data analysis. A total of 27 differentially selected regions (DSRs), including 1 on Sus scrofa chromosome 1 (SSC1), 1 on SSC4, 2 on SSC5, 4 on SSC6, 2 on SSC7, 5 on SSC8, 3 on SSC9, 1 on SSC14, 3 on SSC18, and 5 on SSCX, were identified to show the genome wide separations between the HEBV and LEBV groups for IBW in piglets. A DSR with the most number of significant SNPs (including 7 top 0.1% and 31 top 5% SNPs) was located on SSC6, while another DSR with the largest genetic differences in F_ST was found on SSC18. These regions harbor known functionally important genes involved in growth and development, such as TNFRSF9 (tumor necrosis factor receptor superfamily member 9), CA6 (carbonic anhydrase VI) and MDFIC (MyoD family inhibitor domain containing). A DSR rich in imprinting genes appeared on SSC9, which included PEG10 (paternally expressed 10), SGCE (sarcoglycan, epsilon), PPP1R9A (protein phosphatase 1, regulatory subunit 9A) and ASB4 (ankyrin repeat and SOCS box containing 4). More importantly, our present study provided evidence to support six quantitative trait loci (QTL) regions for pig birth weight, six QTL regions for average birth weight (ABW) and three QTL regions for litter birth weight (LBW) reported previously by other groups. Furthermore, gene ontology analysis with 183 genes harbored in these 27 DSRs suggested that protein, metal, ion and ATP binding, viral process and innate immune response present important pathways for deciphering their roles in fetal growth or development. Overall, our study provides useful information on candidate genes and pathways for regulating birth weight in piglets, thus improving our understanding of the genetic mechanisms involved in porcine embryonic or fetal development.

Genome-wide association study of meat quality traits in a White Duroc×Erhualian F2 intercross and Chinese Sutai pigs

Thousands of QTLs for meat quality traits have been identified by linkage mapping studies, but most of them lack precise position or replication between populations, which hinder their application in pig breeding programs. To localize QTLs for meat quality traits to precise genomic regions, we performed a genome-wide association (GWA) study using the Illumina PorcineSNP60K Beadchip in two swine populations: 434 Sutai pigs and 933 F2 pigs from a White Duroc×Erhualian intercross. Meat quality traits, including pH, color, drip loss, moisture content, protein content and intramuscular fat content (IMF), marbling and firmness scores in the M. longissimus (LM) and M. semimembranosus (SM) muscles, were recorded on the two populations. In total, 127 chromosome-wide significant SNPs for these traits were identified. Among them, 11 SNPs reached genome-wise significance level, including 1 on SSC3 for pH, 1 on SSC3 and 3 on SSC15 for drip loss, 3 (unmapped) for color a*, and 2 for IMF each on SSC9 and SSCX. Except for 11 unmapped SNPs, 116 significant SNPs fell into 28 genomic regions of approximately 10 Mb or less. Most of these regions corresponded to previously reported QTL regions and spanned smaller intervals than before. The loci on SSC3 and SSC7 appeared to have pleiotropic effects on several related traits. Besides them, a few QTL signals were replicated between the two populations. Further, we identified thirteen new candidate genes for IMF, marbling and firmness, on the basis of their positions, functional annotations and reported expression patterns. The findings will contribute to further identification of the causal mutation underlying these QTLs and future marker-assisted selection in pigs.

A genome-wide association study to detect QTL for commercially important traits in Swiss Large White boars

The improvement of meat quality and production traits has high priority in the pork industry. Many of these traits show a low to moderate heritability and are difficult and expensive to measure. Their improvement by targeted breeding programs is challenging and requires knowledge of the genetic and molecular background. For this study we genotyped 192 artificial insemination boars of a commercial line derived from the Swiss Large White breed using the PorcineSNP60 BeadChip with 62,163 evenly spaced SNPs across the pig genome. We obtained 26 estimated breeding values (EBVs) for various traits including exterior, meat quality, reproduction, and production. The subsequent genome-wide association analysis allowed us to identify four QTL with suggestive significance for three of these traits (p-values ranging from 4.99×10⁻⁶ to 2.73×10⁻⁵). Single QTL for the EBVs pH one hour post mortem (pH1) and carcass length were on pig chromosome (SSC) 14 and SSC 2, respectively. Two QTL for the EBV rear view hind legs were on SSC 10 and SSC 16.

Heritability and Bayesian genome-wide association study of first service conception and pregnancy in Brangus heifers

Brangus [3/8 Brahman (Bos indicus) × 5/8 Angus (Bos taurus); n ≈ 800] heifers from 67 sires were used to estimate heritability and conduct a genome-wide association study (GWAS) for 2 binary fertility traits: first service conception (FSC) and heifer pregnancy (HPG). Genotypes were from 53,692 loci on the BovineSNP50 (Infinium Bead Chips, Illumina, San Diego, CA). Yearling heifers were estrous synchronized, bred by AI, and then exposed to natural service breeding. Reproductive ultrasound and DNA-based parentage testing were used to determine if the heifer conceived by AI or natural service, and code for FSC and HPG traits. Success rates for FSC and HPG were 53.3% and 78.0% ± 0.01%, and corresponding heritability estimates were 0.18 ± 0.07 and 0.10 ± 0.06, respectively. The models used in obtaining these heritability estimates and GWAS included fixed effects of year (i.e., 2005 to 2007), birth location, calving season, age of dam, and contemporary group. In GWAS, simultaneous associations of 1 Mb SNP windows with phenotype were undertaken with Bayes C analyses using GenSel software. The 1 Mb windows contained 21.3 ± 1.1 SNP. Analyses fitted a mixture model that treated SNP effects as random, with an assumed fraction pi = 0.9995 having no effect on phenotype. The windows that accounted for 1.0% of genetic variance were considered as QTL associated with FSC or HPG. Eighteen QTL existed on 15 chromosomes for the 2 traits. On average, each QTL accounted for 2.43% ± 0.2% of the genetic variance. Chromosome 8 harbored 2 QTL for FSC and 1 for HPG; however, these regions did not overlap. Chromosomes 3, 15, 16, 19, 24, 26, 27, 29, and X included QTL only for FSC, whereas chromosomes 2, 4, 10, 13, and 20 contained QTL only for HPG. The multitude of QTL detected for FSC and HPG in this GWAS involving Brangus heifers exemplifies the complex regulation of variation in heifer fertility traits of low heritability.