Large scale sequence-based screen for recessive variants allows for identification and monitoring of rare deleterious variants in pigs

Most deleterious variants are recessive and segregate at relatively low frequency. Therefore, high sample sizes are required to identify these variants. In this study we report a large-scale sequence based genome-wide association study (GWAS) in pigs, with a total of 120,000 Large White and 80,000 Synthetic breed animals imputed to sequence using a reference population of approximately 1,100 whole genome sequenced pigs. We imputed over 20 million variants with high accuracies (R2>0.9) even for low frequency variants (1-5% minor allele frequency). This sequence-based analysis revealed a total of 14 additive and 9 non-additive significant quantitative trait loci (QTLs) for growth rate and backfat thickness. With the non-additive (recessive) model, we identified a deleterious missense SNP in the CDHR2 gene reducing growth rate and backfat in homozygous Large White animals. For the Synthetic breed, we revealed a QTL on chromosome 15 with a frameshift variant in the OBSL1 gene. This QTL has a major impact on both growth rate and backfat, resembling human 3M-syndrome 2 which is related to the same gene. With the additive model, we confirmed known QTLs on chromosomes 1 and 5 for both breeds, including variants in the MC4R and CCND2 genes. On chromosome 1, we disentangled a complex QTL region with multiple variants affecting both traits, harboring 4 independent QTLs in the span of 5 Mb. Together we present a large scale sequence-based association study that provides a key resource to scan for novel variants at high resolution for breeding and to further reduce the frequency of deleterious alleles at an early stage in the breeding program.

Identification and characterisation of de novo germline structural variants in two commercial pig lines using trio-based whole genome sequencing

BACKGROUND

De novo mutations arising in the germline are a source of genetic variation and their discovery broadens our understanding of genetic disorders and evolutionary patterns. Although the number of de novo single nucleotide variants (dnSNVs) has been studied in a number of species, relatively little is known about the occurrence of de novo structural variants (dnSVs). In this study, we investigated 37 deeply sequenced pig trios from two commercial lines to identify dnSVs present in the offspring. The identified dnSVs were characterised by identifying their parent of origin, their functional annotations and characterizing sequence homology at the breakpoints.

RESULTS

We identified four swine germline dnSVs, all located in intronic regions of protein-coding genes. Our conservative, first estimate of the swine germline dnSV rate is 0.108 (95% CI 0.038-0.255) per generation (one dnSV per nine offspring), detected using short-read sequencing. Two detected dnSVs are clusters of mutations. Mutation cluster 1 contains a de novo duplication, a dnSNV and a de novo deletion. Mutation cluster 2 contains a de novo deletion and three de novo duplications, of which one is inverted. Mutation cluster 2 is 25 kb in size, whereas mutation cluster 1 (197 bp) and the other two individual dnSVs (64 and 573 bp) are smaller. Only mutation cluster 2 could be phased and is located on the paternal haplotype. Mutation cluster 2 originates from both micro-homology as well as non-homology mutation mechanisms, where mutation cluster 1 and the other two dnSVs are caused by mutation mechanisms lacking sequence homology. The 64 bp deletion and mutation cluster 1 were validated through PCR. Lastly, the 64 bp deletion and the 573 bp duplication were validated in sequenced offspring of probands with three generations of sequence data.

CONCLUSIONS

Our estimate of 0.108 dnSVs per generation in the swine germline is conservative, due to our small sample size and restricted possibilities of dnSV detection from short-read sequencing. The current study highlights the complexity of dnSVs and shows the potential of breeding programs for pigs and livestock species in general, to provide a suitable population structure for identification and characterisation of dnSVs.

Fine Mapping of a Major Backfat QTL Reveals a Causal Regulatory Variant Affecting the CCND2 Gene

Backfat is an important trait in pork production, and it has been included in the breeding objectives of genetic companies for decades. Although adipose tissue is a good energy storage, excessive fat results in reduced efficiency and economical losses. A large QTL for backfat thickness on chromosome 5 is still segregating in different commercial pig breeds. We fine mapped this QTL region using a genome-wide association analysis (GWAS) with 133,358 genotyped animals from five commercial populations (Landrace, Pietrain, Large White, Synthetic, and Duroc) imputed to the porcine 660K SNP chip. The lead SNP was located at 5:66103958 (G/A) within the third intron of the CCND2 gene, with the G allele associated with more backfat, while the A allele is associated with less backfat. We further phased the QTL region to discover a core haplotype of five SNPs associated with low backfat across three breeds. Linkage disequilibrium analysis using whole-genome sequence data revealed three candidate causal variants within intronic regions and downstream of the CCND2 gene, including the lead SNP. We evaluated the association of the lead SNP with the expression of the genes in the QTL region (including CCND2) in a large cohort of 100 crossbred samples, sequenced in four different tissues (lung, spleen, liver, muscle). Results show that the A allele increases the expression of CCND2 in an additive way in three out of four tissues. Our findings indicate that the causal variant for this QTL region is a regulatory variant within the third intron of the CCND2 gene affecting the expression of CCND2.

Using short read sequencing to characterise balanced reciprocal translocations in pigs

Background

A balanced constitutional reciprocal translocation (RT) is a mutual exchange of terminal segments of two non-homologous chromosomes without any loss or gain of DNA in germline cells. Carriers of balanced RTs are viable individuals with no apparent phenotypical consequences. These animals produce, however, unbalanced gametes and show therefore reduced fertility and offspring with congenital abnormalities. This cytogenetic abnormality is usually detected using chromosome staining techniques. The aim of this study was to test the possibilities of using paired end short read sequencing for detection of balanced RTs in boars and investigate their breakpoints and junctions.

Results

Balanced RTs were recovered in a blinded analysis, using structural variant calling software DELLY, in 6 of the 7 carriers with 30 fold short read paired end sequencing. In 15 non-carriers we did not detect any RTs. Reducing the coverage to 20 fold, 15 fold and 10 fold showed that at least 20 fold coverage is required to obtain good results. One RT was not detected using the blind screening, however, a highly likely RT was discovered after unblinding. This RT was located in a repetitive region, showing the limitations of short read sequence data. The detailed analysis of the breakpoints and junctions suggested three junctions showing microhomology, three junctions with blunt-end ligation, and three micro-insertions at the breakpoint junctions. The RTs detected also showed to disrupt genes.

Conclusions

We conclude that paired end short read sequence data can be used to detect and characterize balanced reciprocal translocations, if sequencing depth is at least 20 fold coverage. However, translocations in repetitive areas may require large fragments or even long read sequence data.

Detection of a Frameshift Deletion in the SPTBN4 Gene Leads to Prevention of Severe Myopathy and Postnatal Mortality in Pigs

Piglet mortality is a complex phenotype that depends on the environment, selection on piglet health, but also on the interaction between the piglet and sow. However, also monogenic recessive defects contribute to piglet mortality. Selective breeding has decreased overall piglet mortality by improving both mothering abilities and piglet viability. However, variants underlying recessive monogenic defects are usually not well captured within the breeding values, potentially drifting to higher frequency as a result of intense selection or genetic drift. This study describes the identification by whole-genome sequencing of a recessive 16-bp deletion in the SPTBN4 gene causing postnatal mortality in a pig breeding line. The deletion induces a frameshift and a premature stop codon, producing an impaired and truncated spectrin beta non-erythrocytic 4 protein (SPTBN4). Applying medium density single nucleotide polymorphism (SNP) data available for all breeding animals, a pregnant carrier sow sired by a carrier boar was identified. Of the resulting piglets, two confirmed homozygous piglets suffered from severe myopathy, hind-limb paralysis, and tremors. Histopathological examination showed dispersed degeneration and decrease of cross-striations in the dorsal and hind-limb muscle fibers of the affected piglets. Hence, the affected piglets are unable to walk or drink, usually resulting in death within a few hours after birth. This study demonstrates how growing genomic resources in pig breeding can be applied to identify rare syndromes in breeding populations, that are usually poorly documented and often are not even known to have a genetic basis. The study allows to prevent carrier-by-carrier matings, thereby gradually decreasing the frequency of the detrimental allele and avoiding the birth of affected piglets, improving animal welfare. Finally, these "natural knockouts" increase our understanding of gene function within the mammalian clade, and provide a potential model for human disease.

A QTL for Number of Teats Shows Breed Specific Effects on Number of Vertebrae in Pigs: Bridging the Gap Between Molecular and Quantitative Genetics

Modern breeding schemes for livestock species accumulate a large amount of genotype and phenotype data which can be used for genome-wide association studies (GWAS). Many chromosomal regions harboring effects on quantitative traits have been reported from these studies, but the underlying causative mutations remain mostly undetected. In this study, we combine large genotype and phenotype data available from a commercial pig breeding scheme for three different breeds (Duroc, Landrace, and Large White) to pinpoint functional variation for a region on porcine chromosome 7 affecting number of teats (NTE). Our results show that refining trait definition by counting number of vertebrae (NVE) and ribs (RIB) helps to reduce noise from other genetic variation and increases heritability from 0.28 up to 0.62 NVE and 0.78 RIB in Duroc. However, in Landrace, the effect of the same QTL on NTE mainly affects NVE and not RIB, which is reflected in reduced heritability for RIB (0.24) compared to NVE (0.59). Further, differences in allele frequencies and accuracy of rib counting influence genetic parameters. Correction for the top SNP does not detect any other QTL effect on NTE, NVE, or RIB in Landrace or Duroc. At the molecular level, haplotypes derived from 660K SNP data detects a core haplotype of seven SNPs in Duroc. Sequence analysis of 16 Duroc animals shows that two functional mutations of the Vertnin (VRTN) gene known to increase number of thoracic vertebrae (ribs) reside on this haplotype. In Landrace, the linkage disequilibrium (LD) extends over a region of more than 3 Mb also containing both VRTN mutations. Here, other modifying loci are expected to cause the breed-specific effect. Additional variants found on the wildtype haplotype surrounding the VRTN region in all sequenced Landrace animals point toward breed specific differences which are expected to be present also across the whole genome. This Landrace specific haplotype contains two missense mutations in the ABCD4 gene, one of which is expected to have a negative effect on the protein function. Together, the integration of largescale genotype, phenotype and sequence data shows exemplarily how population parameters are influenced by underlying variation at the molecular level.

Loss of function mutations in essential genes cause embryonic lethality in pigs

Lethal recessive alleles cause pre- or postnatal death in homozygous affected individuals, reducing fertility. Especially in small size domestic and wild populations, those alleles might be exposed by inbreeding, caused by matings between related parents that inherited the same recessive lethal allele from a common ancestor. In this study we report five relatively common (up to 13.4% carrier frequency) recessive lethal haplotypes in two commercial pig populations. The lethal haplotypes have a large effect on carrier-by-carrier matings, decreasing litter sizes by 15.1 to 21.6%. The causal mutations are of different type including two splice-site variants (affecting POLR1B and TADA2A genes), one frameshift (URB1), and one missense (PNKP) variant, resulting in a complete loss-of-function of these essential genes. The recessive lethal alleles affect up to 2.9% of the litters within a single population and are responsible for the death of 0.52% of the total population of embryos. Moreover, we provide compelling evidence that the identified embryonic lethal alleles contribute to the observed heterosis effect for fertility (i.e. larger litters in crossbred offspring). Together, this work marks specific recessive lethal variation describing its functional consequences at the molecular, phenotypic, and population level, providing a unique model to better understand fertility and heterosis in livestock.

Lethal recessives are mutations that cause early lethality in homozygous state that usually occur at very low frequency in wild and domestic populations. In livestock, however, those mutations might become more prevalent as a result of inbreeding. In this study, we report five such recessive lethal haplotypes that cause embryonic lethality in homozygous state in pigs. The causal mutations are of different type but all destroy the structure of essential genes involved in cellular housekeeping processes, essential for embryonic development. The lethal recessives have substantial impact on the population fitness affecting up to 3% of the population litters, causing the death of 0.52% of the total population of embryos. Moreover, these 'natural knockouts' can increase understanding of gene function within the mammalian clade. Together, our study will allow monitoring, and facilitate the purging and partial elimination of recessive lethal mutations in frequently used pig breeds.

Balancing selection on a recessive lethal deletion with pleiotropic effects on two neighboring genes in the porcine genome

Livestock populations can be used to study recessive defects caused by deleterious alleles. The frequency of deleterious alleles including recessive lethal alleles can stay at high or moderate frequency within a population, especially if recessive lethal alleles exhibit an advantage for favourable traits in heterozygotes. In this study, we report such a recessive lethal deletion of 212kb (del) within the BBS9 gene in a breeding population of pigs. The deletion produces a truncated BBS9 protein expected to cause a complete loss-of-function, and we find a reduction of approximately 20% on the total number of piglets born from carrier by carrier matings. Homozygous del/del animals die mid- to late-gestation, as observed from high increase in numbers of mummified piglets resulting from carrier-by-carrier crosses. The moderate 10.8% carrier frequency (5.4% allele frequency) in this pig population suggests an advantage on a favourable trait in heterozygotes. Indeed, heterozygous carriers exhibit increased growth rate, an important selection trait in pig breeding. Increased growth and appetite together with a lower birth weight for carriers of the BBS9 null allele in pigs is analogous to the phenotype described in human and mouse for (naturally occurring) BBS9 null-mutants. We show that fetal death, however, is induced by reduced expression of the downstream BMPER gene, an essential gene for normal foetal development. In conclusion, this study describes a lethal 212kb deletion with pleiotropic effects on two different genes, one resulting in fetal death in homozygous state (BMPER), and the other increasing growth (BBS9) in heterozygous state. We provide strong evidence for balancing selection resulting in an unexpected high frequency of a lethal allele in the population. This study shows that the large amounts of genomic and phenotypic data routinely generated in modern commercial breeding programs deliver a powerful tool to monitor and control lethal alleles much more efficiently.

We report a large deletion within the BBS9 gene that induces late fetal mortality in homozygous affected animals in a commercial pig population. This late fetal mortality causes the fetus to become encapsulated and desiccated during the remaining time of the pregnancy, a process called mummification. The unusually high carrier frequency for this lethal deletion (10.8%) likely results from its strong positive association with growth rate in heterozygous individuals, an important selection trait in the pig breeding industry. Interestingly, we show that the positive effect on growth is induced by a heterozygous loss-of-function of the BBS9 gene, associated with obesity in human and mouse. However, late fetal mortality is induced by insufficient expression of the BMPER gene located directly downstream of the deletion which affects its regulatory elements required for gene expression. Together, our study shows an unique example of allelic pleiotropy in which one allele (deletion) is responsible for both increased growth and late fetal mortality by affecting two different genes.

Weighted single-step GWAS and gene network analysis reveal new candidate genes for semen traits in pigs

Background

In recent years, there has been increased interest in the study of the molecular processes that affect semen traits. In this study, our aim was to identify quantitative trait loci (QTL) regions associated with four semen traits (motility, progressive motility, number of sperm cells per ejaculate and total morphological defects) in two commercial pig lines (L1: Large White type and L2: Landrace type). Since the number of animals with both phenotypes and genotypes was relatively small in our dataset, we conducted a weighted single-step genome-wide association study, which also allows unequal variances for single nucleotide polymorphisms. In addition, our aim was also to identify candidate genes within QTL regions that explained the highest proportions of genetic variance. Subsequently, we performed gene network analyses to investigate the biological processes shared by genes that were identified for the same semen traits across lines.

Results

We identified QTL regions that explained up to 10.8% of the genetic variance of the semen traits on 12 chromosomes in L1 and 11 chromosomes in L2. Sixteen QTL regions in L1 and six QTL regions in L2 were associated with two or more traits within the population. Candidate genes SCN8A, PTGS2, PLA2G4A, DNAI2, IQCG and LOC102167830 were identified in L1 and NME5, AZIN2, SPATA7, METTL3 and HPGDS in L2. No regions overlapped between these two lines. However, the gene network analysis for progressive motility revealed two genes in L1 (PLA2G4A and PTGS2) and one gene in L2 (HPGDS) that were involved in two biological processes i.e. eicosanoid biosynthesis and arachidonic acid metabolism. PTGS2 and HPGDS were also involved in the cyclooxygenase pathway.

Conclusions

We identified several QTL regions associated with semen traits in two pig lines, which confirms the assumption of a complex genetic determinism for these traits. A large part of the genetic variance of the semen traits under study was explained by different genes in the two evaluated lines. Nevertheless, the gene network analysis revealed candidate genes that are involved in shared biological pathways that occur in mammalian testes, in both lines.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0412-z) contains supplementary material, which is available to authorized users.

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.

Artificial selection on introduced Asian haplotypes shaped the genetic architecture in European commercial pigs

Early pig farmers in Europe imported Asian pigs to cross with their local breeds in order to improve traits of commercial interest. Current genomics techniques enabled genome-wide identification of these Asian introgressed haplotypes in modern European pig breeds. We propose that the Asian variants are still present because they affect phenotypes that were important for ancient traditional, as well as recent, commercial pig breeding. Genome-wide introgression levels were only weakly correlated with gene content and recombination frequency. However, regions with an excess or absence of Asian haplotypes (AS) contained genes that were previously identified as phenotypically important such as FASN, ME1, and KIT. Therefore, the Asian alleles are thought to have an effect on phenotypes that were historically under selection. We aimed to estimate the effect of AS in introgressed regions in Large White pigs on the traits of backfat (BF) and litter size. The majority of regions we tested that retained Asian deoxyribonucleic acid (DNA) showed significantly increased BF from the Asian alleles. Our results suggest that the introgression in Large White pigs has been strongly determined by the selective pressure acting upon the introgressed AS. We therefore conclude that human-driven hybridization and selection contributed to the genomic architecture of these commercial pigs.

Accounting for genetic architecture in single- and multipopulation genomic prediction using weights from genomewide association studies in pigs

We studied the effect of including GWAS results on the accuracy of single- and multipopulation genomic predictions. Phenotypes (backfat thickness) and genotypes of animals from two sire lines (SL1, n = 1146 and SL3, n = 1264) were used in the analyses. First, GWAS were conducted for each line and for a combined data set (both lines together) to estimate the genetic variance explained by each SNP. These estimates were used to build matrices of weights (D), which was incorporated into a GBLUP method. Single population evaluated with traditional GBLUP had accuracies of 0.30 for SL1 and 0.31 for SL3. When weights were employed in GBLUP, the accuracies for both lines increased (0.32 for SL1 and 0.34 for SL3). When a multipopulation reference set was used in GBLUP, the accuracies were higher (0.36 for SL1 and 0.32 for SL3) than in single-population prediction. In addition, putting together the multipopulation reference set and the weights from the combined GWAS provided even higher accuracies (0.37 for SL1, and 0.34 for SL3). The use of multipopulation predictions and weights estimated from a combined GWAS increased the accuracy of genomic predictions.

Linkage disequilibrium patterns and persistence of phase in purebred and crossbred pig (Sus scrofa) populations

Background

Genomic selection and genomic wide association studies are widely used methods that aim to exploit the linkage disequilibrium (LD) between markers and quantitative trait loci (QTL). Securing a sufficiently large set of genotypes and phenotypes can be a limiting factor that may be overcome by combining data from multiple breeds or using crossbred information. However, the estimated effect of a marker in one breed or a crossbred can only be useful for the selection of animals in another breed if there is a correspondence of the phase between the marker and the QTL across breeds. Using data of five pure pig (Sus scrofa) lines (SL1, SL2, SL3, DL1, DL2), one F₁ cross (DLF1) and two commercial finishing crosses (TER1 and TER2), the objectives of this study were: (i) to compare the equality of LD decay curves of different pig populations; and (ii) to evaluate the persistence of the LD phase across lines or final crosses.

Results

Almost all of the lines presented different extents of LD, except for the SL2 and DL3, both of which exhibited the same extent of LD. Similar levels of LD over large distances were found in crossbred and pure lines. The crossbred animals (DLF1, TER1 and TER2) presented a high persistence of phase with their parental lines, suggesting that the available porcine single nucleotide polymorphism (SNP) chip should be dense enough to include markers that have the same LD phase with QTL across crossbred and parental pure lines. The persistence of phase across pure lines varied considerably between the different line comparisons; however, correlations were above 0.8 for all line comparisons when marker distances were smaller than 50 kb.

Conclusions

This study showed that crossbred populations could be very useful as a reference for the selection of pure lines by means of the available SNP chip panel. Here, we also pinpoint pure lines that could be combined in a multiline training population. However, if multiline reference populations are used for genomic selection, the required density of SNP panels should be higher compared with a single breed reference population.

Electronic supplementary material

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

A genome-wide association study reveals a novel candidate gene for sperm motility in pigs

Sperm motility is one of the most widely used parameters in order to evaluate boar semen quality. However, this trait can only be measured after puberty. Thus, the use of genomic information appears as an appealing alternative to evaluate and improve selection for boar fertility traits earlier in life. With this study we aimed to identify SNPs with significant association with sperm motility in two different commercial pig populations and to identify possible candidate genes within the identified QTL regions. We performed a single-SNP genome-wide association study using genotyped animals from a Landrace-based (L1) and a Large White-based (L2) pig populations. For L1, a total of 602 animals genotyped for 42,551 SNPs were used in the association analysis. For L2, a total of 525 animals genotyped for 40,890 SNPs were available. After the association analysis, a false discovery rate q-value ≤0.05 was used as the threshold for significant association. No SNPs were significantly associated with sperm motility in L1, while six SNPs on Sus scrofa chromosome 1 (position 117.26-119.56Mb) were significant in L2. The mitochondrial methionyl-tRNA formyltransferase (MTFMT) gene, which affects translation efficiency of proteins in sperm cells, was identified as a putative candidate gene. The significant markers identified in this study may be useful to enhance the genetic improvement of sperm motility by selection of boars at an earlier age under a marker assisted selection strategy.

Asian low-androstenone haplotype on pig chromosome 6 does not unfavorably affect production and reproduction traits

European pigs that carry Asian haplotypes of a 1.94-Mbp region on pig chromosome 6 have lower levels of androstenone, one of the two main compounds causing boar taint. The objective of our study was to examine potential pleiotropic effects of the Asian low-androstenone haplotypes. A single nucleotide polymorphism marker, rs81308021, distinguishes the Asian from European haplotypes and was used to investigate possible associations of androstenone with production and reproduction traits. Eight traits were available from three European commercial breeds. For the two sow lines studied, a favorable effect on number of teats was detected for the low-androstenone haplotype. In one of these sow lines, a favorable effect on number of spermatozoa per ejaculation was detected for the low-androstenone haplotype. No unfavorable pleiotropic effects were found, which suggests that selection for low-androstenone haplotypes within the 1.94 Mbp would not unfavorably affect the other eight relevant traits.

A genome-wide association study reveals dominance effects on number of teats in pigs

Dominance has been suggested as one of the genetic mechanisms explaining heterosis. However, using traditional quantitative genetic methods it is difficult to obtain accurate estimates of dominance effects. With the availability of dense SNP (Single Nucleotide Polymorphism) panels, we now have new opportunities for the detection and use of dominance at individual loci. Thus, the aim of this study was to detect additive and dominance effects on number of teats (NT), specifically to investigate the importance of dominance in a Landrace-based population of pigs. In total, 1,550 animals, genotyped for 32,911 SNPs, were used in single SNP analysis. SNPs with a significant genetic effect were tested for their mode of gene action being additive, dominant or a combination. In total, 21 SNPs were associated with NT, located in three regions with additive (SSC6, 7 and 12) and one region with dominant effects (SSC4). Estimates of additive effects ranged from 0.24 to 0.29 teats. The dominance effect of the QTL located on SSC4 was negative (−0.26 teats). The additive variance of the four QTLs together explained 7.37% of the total phenotypic variance. The dominance variance of the four QTLs together explained 1.82% of the total phenotypic variance, which corresponds to one-fourth of the variance explained by additive effects. The results suggest that dominance effects play a relevant role in the genetic architecture of NT. The QTL region on SSC7 contains the most promising candidate gene: VRTN. This gene has been suggested to be related to the number of vertebrae, a trait correlated with NT.

Genomic analysis reveals selection for Asian genes in European pigs following human-mediated introgression

The independent domestication of local wild boar populations in Asia and Europe about 10,000 years ago led to distinct European and Asian pig breeds, each with very different phenotypic characteristics. During the Industrial Revolution, Chinese breeds were imported to Europe to improve commercial traits in European breeds. Here we demonstrate the presence of introgressed Asian haplotypes in European domestic pigs and selection signatures on some loci in these regions, using whole genome sequence data. The introgression signatures are widespread and the Asian haplotypes are rarely fixed. The Asian introgressed haplotypes are associated with regions harbouring genes involved in meat quality, development and fertility. We identify Asian-derived non-synonymous mutations in the AHR gene that associate with increased litter size in multiple European commercial lines. These findings demonstrate that increased fertility was an important breeding goal for early nineteenth century pig farmers, and that Asian variants of genes related to this trait were preferentially selected during the development of modern European pig breeds.

High-resolution association mapping of number of teats in pigs reveals regions controlling vertebral development

BACKGROUND

Selection pressure on the number of teats has been applied to be able to provide enough teats for the increase in litter size in pigs. Although many QTL were reported, they cover large chromosomal regions and the functional mutations and their underlying biological mechanisms have not yet been identified. To gain a better insight in the genetic architecture of the trait number of teats, we performed a genome-wide association study by genotyping 936 Large White pigs using the Illumina PorcineSNP60 Beadchip. The analysis is based on deregressed breeding values to account for the dense family structure and a Bayesian approach for estimation of the SNP effects.

RESULTS

The genome-wide association study resulted in 212 significant SNPs. In total, 39 QTL regions were defined including 170 SNPs on 13 Sus scrofa chromosomes (SSC) of which 5 regions on SSC7, 9, 10, 12 and 14 were highly significant. All significantly associated regions together explain 9.5% of the genetic variance where a QTL on SSC7 explains the most genetic variance (2.5%). For the five highly significant QTL regions, a search for candidate genes was performed. The most convincing candidate genes were VRTN and Prox2 on SSC7, MPP7, ARMC4, and MKX on SSC10, and vertebrae δ-EF1 on SSC12. All three QTL contain candidate genes which are known to be associated with vertebral development. In the new QTL regions on SSC9 and SSC14, no obvious candidate genes were identified.

CONCLUSIONS

Five major QTL were found at high resolution on SSC7, 9, 10, 12, and 14 of which the QTL on SSC9 and SSC14 are the first ones to be reported on these chromosomes. The significant SNPs found in this study could be used in selection to increase number of teats in pigs, so that the increasing number of live-born piglets can be nursed by the sow. This study points to common genetic mechanisms regulating number of vertebrae and number of teats.

On the relationship between an Asian haplotype on chromosome 6 that reduces androstenone levels in boars and the differential expression of SULT2A1 in the testis

Background

Androstenone is one of the major compounds responsible for boar taint, a pronounced urine-like odor produced when cooking boar meat. Several studies have identified quantitative trait loci (QTL) for androstenone level on Sus scrofa chromosome (SSC) 6. For one of the candidate genes in the region SULT2A1, a difference in expression levels in the testis has been shown at the protein and RNA level.

Results

Haplotypes were predicted for the QTL region and their effects were estimated showing that haplotype 1 was consistently related with a lower level, and haplotype 2 with a higher level of androstenone. A recombinant haplotype allowed us to narrow down the QTL region from 3.75 Mbp to 1.94 Mbp. An RNA-seq analysis of the liver and testis revealed six genes that were differentially expressed between homozygotes of haplotypes 1 and 2. Genomic sequences of these differentially expressed genes were checked for variations within potential regulatory regions. We identified one variant located within a CpG island that could affect expression of SULT2A1 gene. An allele-specific expression analysis in the testis did not show differential expression between the alleles of SULT2A1 located on the different haplotypes in heterozygous animals. However a synonymous mutation C166T (SSC6: 49,117,861 bp in Sscrofa 10.2; C/T) was identified within the exon 2 of SULT2A1 for which the haplotype 2 only had the C allele which was higher expressed than the T allele, indicating haplotype-independent allelic-imbalanced expression between the two alleles. A phylogenetic analysis for the 1.94 Mbp region revealed that haplotype 1, associated with low androstenone level, originated from Asia.

Conclusions

Differential expression could be observed for six genes by RNA-seq analysis. No difference in the ratio of C:T expression of SULT2A1 for the haplotypes was found by the allele-specific expression analysis, however, a difference in expression between the C over T allele was found for a variation within SULT2A1, showing that the difference in androstenone levels between the haplotypes is not caused by the SNP in exon 2.

Improved estimation of inbreeding and kinship in pigs using optimized SNP panels

Background

Traditional breeding programs consider an average pairwise kinship between sibs. Based on pedigree information, the relationship matrix is used for genetic evaluations disregarding variation due to Mendelian sampling. Therefore, inbreeding and kinship coefficients are either over or underestimated resulting in reduction of accuracy of genetic evaluations and genetic progress. Single nucleotide polymorphism (SNPs) can be used to estimate pairwise kinship and individual inbreeding more accurately. The aim of this study was to optimize the selection of markers and determine the required number of SNPs for estimation of kinship and inbreeding.

Results

A total of 1,565 animals from three commercial pig populations were analyzed for 28,740 SNPs from the PorcineSNP60 Beadchip. Mean genomic inbreeding was higher than pedigree-based estimates in lines 2 and 3, but lower in line 1. As expected, a larger variation of genomic kinship estimates was observed for half and full sibs than for pedigree-based kinship reflecting Mendelian sampling. Genomic kinship between father-offspring pairs was lower (0.23) than the estimate based on pedigree (0.26). Bootstrap analyses using six reduced SNP panels (n = 500, 1000, 1500, 2000, 2500 and 3000) showed that 2,000 SNPs were able to reproduce the results very close to those obtained using the full set of unlinked markers (n = 7,984-10,235) with high correlations (inbreeding r > 0.82 and kinship r > 0.96) and low variation between different sets with the same number of SNPs.

Conclusions

Variation of kinship between sibs due to Mendelian sampling is better captured using genomic information than the pedigree-based method. Therefore, the reduced sets of SNPs could generate more accurate kinship coefficients between sibs than the pedigree-based method. Variation of genomic kinship of father-offspring pairs is recommended as a parameter to determine accuracy of the method rather than correlation with pedigree-based estimates. Inbreeding and kinship coefficients can be estimated with high accuracy using ≥2,000 unlinked SNPs within all three commercial pig lines evaluated. However, a larger number of SNPs might be necessary in other populations or across lines.

A de novo germline mutation in MYH7 causes a progressive dominant myopathy in pigs

Background

About 9% of the offspring of a clinically healthy Piétrain boar named ‘Campus’ showed a progressive postural tremor called Campus syndrome (CPS). Extensive backcross experiments suggested a dominant mode of inheritance, and the founder boar was believed to be a gonadal mosaic. A genome-scan mapped the disease-causing mutation to an 8 cM region of porcine chromosome 7 containing the MHY7 gene. Human distal myopathy type 1 (MPD1), a disease partially resembling CPS in pigs, has been associated with mutations in the MYH7 gene.

Results

The porcine MYH7 gene structure was predicted based on porcine reference genome sequence, porcine mRNA, and in comparison to the human ortholog. The gene structure was highly conserved with the exception of the first exon. Mutation analysis of a contiguous genomic interval of more than 22 kb spanning the complete MYH7 gene revealed an in-frame insertion within exon 30 of MYH7 (c.4320_4321insCCCGCC) which was perfectly associated with the disease phenotype and confirmed the dominant inheritance. The mutation is predicted to insert two amino acids (p.Ala1440_Ala1441insProAla) in a very highly conserved region of the myosin tail. The boar ‘Campus’ was shown to be a germline and somatic mosaic as assessed by the presence of the mutant allele in seven different organs.

Conclusion

This study illustrates the usefulness of recently established genomic resources in pigs. We have identified a spontaneous mutation in MYH7 as the causative mutation for CPS. This paper describes the first case of a disorder caused by a naturally occurring mutation in the MYH7 gene of a non-human mammalian species. Our study confirms the previous classification as a primary myopathy and provides a defined large animal model for human MPD1. We provide evidence that the CPS mutation occurred during the early development of the boar ‘Campus’. Therefore, this study provides an example of germline mosaicism with an asymptomatic founder.

The distal end of porcine chromosome 6p is involved in the regulation of skatole levels in boars

Background

Boar taint is an unpleasant condition of pork, mainly due to the accumulation of androstenone and skatole in male pigs at onset of puberty. This condition is the cause of considerable economic losses in the pig industry and the most common practice to control it is to castrate male piglets. Because of the economic and animal welfare concerns there is interest in developing genetic markers that could be used in selection schemes to decrease the incidence of boar taint. The Porcine 60 K SNP Beadchip was used to genotype 891 pigs from a composite Duroc sire line, for which skatole levels in fat had been collected.

Results

The genome-wide association study revealed that 16 SNPs (single nucleotide polymorphisms) located on the proximal region of chromosome 6 were significantly associated with skatole levels. These SNPs are grouped in three separate clusters located in the initial 6 Mb region of chromosome 6. The differences observed between the homozygote genotypes for SNPs in the three clusters were substantial, including a difference of 102.8 ng/g skatole in melted fat between the homozygotes for the ALGA0107039 marker. Single SNPs explain up to 22% of the phenotypic variance. No obvious candidate genes could be pinpointed in the region, which may be due to the need of further annotation of the pig genome.

Conclusions

This study demonstrated new SNP markers significantly associated with skatole levels in the distal region of chromosome 6p. These markers defined three independent clusters in the region, which contain a low number of protein-coding genes. The considerable differences observed between the homozygous genotypes for several SNPs may be used in future selection schemes to reduce skatole levels in pigs

A single nucleotide polymorphism set for paternal identification to reduce the costs of trait recording in commercial pig breeding

In animal breeding, recording of correct pedigrees is essential to achieve genetic progress. Markers on DNA are useful to verify the on-farm pedigree records (parental verification) but can also be used to assign parents retrospectively (parental identification). This approach could reduce the costs of recording for traits with low incidence, such as those related to diseases or mortality. In this study, SNP were used to assign the true sires of 368 purebred animals from a Duroc-based sire line and 140 crossbred offspring from a commercial pig population. Some of the sires were closely related. There were 3 full sibs and 17 half sibs among the true fathers and 4 full sibs and 35 half sibs among all putative fathers. To define the number of SNP necessary, 5 SNP panels (40, 60, 80, 100, and 120 SNP) were assembled from the Illumina PorcineSNP60 Beadchip (Illumina, San Diego, CA) based on minor allele frequency (>0.3), high genotyping call rate (≥90%), and equal spacing across the genome. For paternal identification considering only the 66 true sires in the data set, 60 SNP resulted in 100% correct assignment of the sire. By including additional putative sires (n = 304), 80 SNP were sufficient for 100% correct assignment of the sire. The following criteria were derived to identify the correct sire for the current data set: the logarithm of odds (LOD) score for assigning the correct sire was ≥5, the number of mismatches was ≤1, and the difference in the LOD score between the first and the second most likely sire was >5. If the correct sire was not present among all putative sires, the mean LOD for the most likely sire was close to zero or negative when using 100 SNP. More SNP would be needed for paternal identification if the number of putative sires increased and the degree of relatedness was greater than in the data set used here. The threshold for the number of mismatches can be adjusted according to the practical situation to account for the trade-off between false negatives and false positives. The latter can be avoided efficiently, ensuring that the correct father is being sampled. Nevertheless, a restriction on the number of putative sires is advisable to reduce the risk of assigning close relatives.

Pig genome sequence--analysis and publication strategy

Background

The pig genome is being sequenced and characterised under the auspices of the Swine Genome Sequencing Consortium. The sequencing strategy followed a hybrid approach combining hierarchical shotgun sequencing of BAC clones and whole genome shotgun sequencing.

Results

Assemblies of the BAC clone derived genome sequence have been annotated using the Pre-Ensembl and Ensembl automated pipelines and made accessible through the Pre-Ensembl/Ensembl browsers. The current annotated genome assembly (Sscrofa9) was released with Ensembl 56 in September 2009. A revised assembly (Sscrofa10) is under construction and will incorporate whole genome shotgun sequence (WGS) data providing > 30× genome coverage. The WGS sequence, most of which comprise short Illumina/Solexa reads, were generated from DNA from the same single Duroc sow as the source of the BAC library from which clones were preferentially selected for sequencing. In accordance with the Bermuda and Fort Lauderdale agreements and the more recent Toronto Statement the data have been released into public sequence repositories (Genbank/EMBL, NCBI/Ensembl trace repositories) in a timely manner and in advance of publication.

Conclusions

In this marker paper, the Swine Genome Sequencing Consortium (SGSC) sets outs its plans for analysis of the pig genome sequence, for the application and publication of the results.

A genome-wide association study on androstenone levels in pigs reveals a cluster of candidate genes on chromosome 6

BACKGROUND

In many countries, male piglets are castrated shortly after birth because a proportion of un-castrated male pigs produce meat with an unpleasant flavour and odour. Main compounds of boar taint are androstenone and skatole. The aim of this high-density genome-wide association study was to identify single nucleotide polymorphisms (SNPs) associated with androstenone levels in a commercial sire line of pigs. The identification of major genetic effects causing boar taint would accelerate the reduction of boar taint through breeding to finally eliminate the need for castration.

RESULTS

The Illumina Porcine 60K+SNP Beadchip was genotyped on 987 pigs divergent for androstenone concentration from a commercial Duroc-based sire line. The association analysis with 47,897 SNPs revealed that androstenone levels in fat tissue were significantly affected by 37 SNPs on pig chromosomes SSC1 and SSC6. Among them, the 5 most significant SNPs explained together 13.7% of the genetic variance in androstenone. On SSC6, a larger region of 10 Mb was shown to be associated with androstenone covering several candidate genes potentially involved in the synthesis and metabolism of androgens. Besides known candidate genes, such as cytochrome P450 A19 (CYP2A19), sulfotransferases SULT2A1, and SULT2B1, also new members of the cytochrome P450 CYP2 gene subfamilies and of the hydroxysteroid-dehydrogenases (HSD17B14) were found. In addition, the gene encoding the ss-chain of the luteinizing hormone (LHB) which induces steroid synthesis in the Leydig cells of the testis at onset of puberty maps to this area on SSC6. Interestingly, the gene encoding the alpha-chain of LH is also located in one of the highly significant areas on SSC1.

CONCLUSIONS

This study reveals several areas of the genome at high resolution responsible for variation of androstenone levels in intact boars. Major genetic factors on SSC1 and SSC6 showing moderate to large effects on androstenone concentration were identified in this commercial breeding line of pigs. Known and new candidate genes cluster especially on SSC6. For one of the most significant SNP variants, the difference in the proportion of animals surpassing the threshold of consumer acceptance between the two homozygous genotypes was as much as 15.6%.

Associations of functional candidate genes derived from gene-expression profiles of prenatal porcine muscle tissue with meat quality and muscle deposition

Ten genes (ANK1, bR10D1, CA3, EPOR, HMGA2, MYPN, NME1, PDGFRA, ERC1, TTN), whose candidacy for meat-quality and carcass traits arises from their differential expression in prenatal muscle development, were examined for association in 1700 performance-tested fattening pigs of commercial purebred and crossbred herds of Duroc, Pietrain, Pietrain x (Landrace x Large White), Duroc x (Landrace x Large White) as well as in an experimental F(2) population based on a reciprocal cross of Duroc and Pietrain. Comparative sequencing revealed polymorphic sites segregating across commercial breeds. Genetic mapping results corresponded to pre-existing assignments to porcine chromosomes or current human-porcine comparative maps. Nine of these genes showed association with meat-quality and carcass traits at a nominal P-value of < or = 0.05; PDGFRA revealed no association reaching the P < or = 0.05 threshold. In particular, HMGA2, CA3, EPOR, NME1 and TTN were associated with meat colour, pH and conductivity of loin 24 h postmortem; CA3 and MYPN exhibited association with ham weight and lean content (FOM) respectively at P-values of < 0.003 that correspond to false discovery rates of < 0.05. However, none of the genes showed significant associations for a particular trait across all populations. The study revealed statistical-genetic evidence for association of the functional candidate genes with traits related to meat quality and muscle deposition. The polymorphisms detected are not likely causal, but markers were identified that are in linkage disequilibrium with causal genetic variation within particular populations.

Congenital syndactyly in cattle: four novel mutations in the low density lipoprotein receptor-related protein 4 gene (LRP4)

Background

Isolated syndactyly in cattle, also known as mulefoot, is inherited as an autosomal recessive trait with variable penetrance in different cattle breeds. Recently, two independent mutations in the bovine LRP4 gene have been reported as the primary cause of syndactyly in the Holstein and Angus cattle breeds.

Results

We confirmed the previously described LRP4 exon 33 two nucleotide substitution in most of the affected Holstein calves and revealed additional evidence for allelic heterogeneity by the identification of four new LRP4 non-synonymous point mutations co-segregating in Holstein, German Simmental and Simmental-Charolais families.

Conclusion

We confirmed a significant role of LRP4 mutations in the pathogenesis of congenital syndactyly in cattle. The newly detected missense mutations in the LRP4 gene represent independent mutations affecting different conserved protein domains. However, the four newly described LRP4 mutations do still not explain all analyzed cases of syndactyly.

Variance component analysis of quantitative trait loci for pork carcass composition and meat quality on SSC4 and SSC11

In a previous study, QTL for carcass composition and meat quality were identified in a commercial finisher cross. The main objective of the current study was to confirm and fine map the QTL on SSC4 and SSC11 by genotyping an increased number of individuals and markers and to analyze the data using a combined linkage and linkage disequilibrium analysis method. A modified version of the method excludes linkage disequilibrium information from the analysis, enabling the comparison of results based on linkage information only or results based on combined linkage and linkage disequilibrium information. Nine additional paternal half-sib families were genotyped for 18 markers, resulting in a total of 1,855 animals genotyped for 15 and 13 markers on SSC4 and SSC11, respectively. The QTL affecting meat color on SSC4 was confirmed, whereas the QTL affecting LM weight could not be confirmed. The combined linkage and linkage disequilibrium analysis resulted in the identification of new significant effects for 14 traits on the 2 chromosomes. Heritabilities of the QTL effects ranged from 1.8 to 13.2%. The analysis contributed to a more accurate positioning of QTL and further characterized their phenotypic effect. However, results showed that even greater marker densities are required to take full advantage of linkage disequilibrium information and to identify haplotypes associated with favorable QTL alleles.

Identification of quantitative trait loci for carcass composition and pork quality traits in a commercial finishing cross

A QTL study for carcass composition and meat quality traits was conducted on finisher pigs of a cross between a synthetic Piétrain/Large White boar line and a commercial sow cross. The mapping population comprised 715 individuals evaluated for a total of 30 traits related to growth and fatness (4 traits), carcass composition (11 traits), and meat quality (15 traits). Offspring of 8 sires (n = 715) were used for linkage analysis and genotyped for 73 microsatellite markers covering 14 chromosomal regions representing approximately 50% of the pig genome. The regions examined were selected based on previous studies suggesting the presence of QTL affecting carcass composition or meat quality traits. Thirty-two QTL exceeding the 5% chromosome-wise significance level were identified. Among these, 5 QTL affecting 5 different traits were significant at the 1% chromosome-wise level. The greatest significance levels were found for a QTL affecting loin weight on SSC11 and a QTL with an effect on the Japanese color scale score of the loin on SSC4. About one-third of the identified QTL were in agreement with QTL previously reported. Results showed that QTL affecting carcass composition and meat quality traits segregated within commercial lines. Use of these results for marker-assisted selection offers opportunities for improving pork quality by within-line selection.

Genetic diversity within and between European pig breeds using microsatellite markers

An important prerequisite for a conservation programme is a comprehensive description of genetic diversity. The aim of this study was to use anonymous genetic markers to assess the between- and the within-population components of genetic diversity for European pig breeds at the scale of the whole continent using microsatellites. Fifty-eight European pig breeds and lines were analysed including local breeds, national varieties of international breeds and commercial lines. A sample of the Chinese Meishan breed was also included. Eleven additional breeds from a previous project were added for some analyses. Approximately 50 individuals per breed were genotyped for a maximum of 50 microsatellite loci. Substantial within-breed variability was observed, with the average expected heterozygosity and observed number of alleles per locus being 0.56 [range 0.43-0.68] and 4.5 respectively. Genotypic frequencies departed from Hardy-Weinberg expectations (P < 0.01) in 15 European populations, with an excess of homozygotes in 12 of them. The European breeds were on average genetically very distinct, with a Wright F(ST) index value of 0.21. The Neighbour-Joining tree drawn from the Reynolds distances among the breeds showed that the national varieties of major breeds and the commercial lines were mostly clustered around their breeds of reference (Duroc, Hampshire, Landrace, Large White and Piétrain). In contrast, local breeds, with the exception of the Iberian breeds, exhibited a star-like topology. The results are discussed in the light of various forces, which may have driven the recent evolution of European pig breeds. This study has consequences for the interpretation of biodiversity results and will be of importance for future conservation programmes.

Genetic diversity in European pigs utilizing amplified fragment length polymorphism markers

The use of DNA markers to evaluate genetic diversity is an important component of the management of animal genetic resources. The Food and Agriculture Organisation of the United Nations (FAO) has published a list of recommended microsatellite markers for such studies; however, other markers are potential alternatives. This paper describes results obtained with a set of amplified fragment length polymorphism (AFLP) markers as part of a genetic diversity study of European pig breeds that also utilized microsatellite markers. Data from 148 AFLP markers genotyped across samples from 58 European and one Chinese breed were analysed. The results were compared with previous analyses of data from 50 microsatellite markers genotyped on the same animals. The AFLP markers had an average within-breed heterozygosity of 0.124 but there was wide variation, with individual markers being monomorphic in 3-98% of the populations. The biallelic and dominant nature of AFLP markers creates a challenge for their use in genetic diversity studies as each individual marker contains limited information and AFLPs only provide indirect estimates of the allelic frequencies that are needed to estimate genetic distances. Nonetheless, AFLP marker-based characterization of genetic distances was consistent with expectations based on breed and regional distributions and produced a similar pattern to that obtained with microsatellites. Thus, data from AFLP markers can be combined with microsatellite data for measuring genetic diversity.

Genetic diversity analysis using lowly polymorphic dominant markers: the example of AFLP in pigs

DNA markers are commonly used for large-scale evaluation of genetic diversity in farm animals, as a component of the management of animal genetic resources. AFLP markers are useful for such studies as they can be generated relatively simply; however, challenges in analysis arise from their dominant scoring and the low level of polymorphism of some markers. This paper describes the results obtained with a set of AFLP markers in a study of 59 pig breeds. AFLP fingerprints were generated using four primer combinations (PC), yielding a total of 148 marker loci, and average harmonic mean of breed sample size was 37.3. The average proportion of monomorphic populations was 63% (range across loci: 3%-98%). The moment-based method of Hill and Weir (2004, Mol Ecol 13:895-908) was applied to estimate gene frequencies, gene diversity (F(ST)), and Reynolds genetic distances. A highly significant average F(ST) of 0.11 was estimated, together with highly significant PC effects on gene diversity. The variance of F(ST) across loci also significantly exceeded the variance expected under the hypothesis of AFLP neutrality, strongly suggesting the sensitivity of AFLP to selection or other forces. Moment estimates were compared to estimates derived from the square root estimation of gene frequency, as currently applied for dominant markers, and the biases incurred in the latter method were evaluated. The paper discusses the hypotheses underlying the moment estimations and various issues relating to the biallelic, dominant, and lowly polymorphic nature of this set of AFLP markers and to their use as compared to microsatellites for measuring genetic diversity.

Genomics for food safety and sustainable animal production

There is a growing concern in society about the safety of animal-derived food, the health and welfare of farm animals and the sustainability of current animal production systems. Along farm animal, breeding genomics may contribute to a solution for these concerns. The use of genomic analysis tools, to achieve genetic progress in typical out-bred populations of farm animals, seems to be more difficult compared to 'model' organisms or plants. However, identification of positional candidate genes may be accelerated by linkage disequilibrium (LD) mapping. Recording of sustainable traits requires a large financial and logistic input and the economic advantages for the market are not as clear as for traditional selection traits. Examples show that the major genes causing variability for similar traits in different species are rarely the same. Therefore, for breeding purposes genomic analysis of the species of interest is needed. The fundamental knowledge obtained on the genetic architecture of complex traits will open new perspectives for the use of DNA tests in selection schemes. For food safety and traceability, DNA-based techniques evolve for monitoring and early warning systems.

Genome scan reveals new coat color loci in exotic pig cross

The porcine genome was scanned to identify loci affecting coat color in an experimental cross between the Meishan breed and Dutch commercial lines. Linkage was studied in 1181 F(2) animals for 132 microsatellite markers and seven binary coat color scores: White, Black spotting, Speckle, Gray, Black, and specific color phenotypes for head and legs. The analyses were performed using interval mapping under various models. The study confirmed the existence of coat color loci on chromosome 8 and chromosome 6. One additional locus affecting White was detected on chromosome 5, possibly representing the porcine equivalent of the steel factor. Two new loci affecting Black were detected on chromosome 2. One of these showed exclusive maternal expression and mapped to a region where imprinted genes have been reported. The effect of the binary coding was tested by additional analyses excluding the white animals (>50% of F(2) animals). This showed that Black spotting was strongly influenced by the locus on chromosome 6 and the other color phenotypes were mainly influenced by the locus on chromosome 8. Epistatic effects were found between the loci on chromosomes 6 and 8 for Black spotting. For Black color, all combinations among chromosomes 2, 6, and 8 showed epistatic effects.

Detection and characterization of quantitative trait loci for meat quality traits in pigs

In an experimental cross between Meishan and Dutch Large White and Landrace lines, 785 F2 animals with carcass information and their parents were typed for molecular markers covering the entire porcine genome. Linkage was studied between these markers and eight meat quality traits. Quantitative trait locus analyses were performed using interval mapping by regression under two genetic models: 1) the line-cross approach, where the founder lines were assumed to be fixed for different QTL alleles and 2) a half-sib model where a unique allele substitution effect was fitted within each of the 38 half-sib families. The line-cross approach included tests for genomic imprinting and sex-specific QTL effects. In total, three genome-wide significant and 26 suggestive QTL were detected. The significant QTL on chromosomes 3, 4, and 13, affecting meat color, were only detected under the half-sib model. Failure of the line-cross approach to detect the meat color QTL suggests that the founder lines have similar allele frequencies for these QTL. This study provides information on new QTL affecting meat quality traits. It also shows the benefit of analyzing experimental data under different genetic and statistical models.

Improving the comparative map of SSC2p-q13 by sample sequencing of BAC clones

To improve the comparative map for pig chromosome 2 and increase the gene density on this chromosome, a porcine bacterial artificial chromosome (BAC) library was screened with 17 microsatellite markers and 18 genes previously assigned to pig chromosome 2. Fifty-one BAC clones located in the region of a maternally imprinted quantitative trait locus for backfat thickness (BFT) were identified. From these BACs 372 kb were sample sequenced. The average read length of a subclone was 442 basepair (bp). Contig assembly analysis showed that every bp was sequenced 1.28 times. Subsequently, sequences were compared with sequences in the nucleotide databases to identify homology with other mammalian sequences. Sequence identity was observed with sequences derived from 35 BACs. The average percentage identity with human sequences was 87.6%, with an average length of 143 bp. In total, sample sequencing of all BACs resulted in sequence identity with 29 human genes, 13 human expressed sequence tags (ESTs), 17 human genomic clones, one rat gene, one porcine gene and nine porcine ESTs. Eighteen genes located on human chromosome 11 and 19, and seven genes from other human locations, one rat gene and one porcine gene were assigned to pig chromosome 2 for the first time. The new genes were added to the radiation hybrid map at the same position as the locus from which the BAC that was sequenced was derived. In total 57 genes were placed on the radiation hybrid map of SSC2p-q13.

A whole-genome scan for quantitative trait loci affecting teat number in pigs

A whole-genome scan was conducted using 132 microsatellite markers to identify chromosomal regions that have an effect on teat number. For this purpose, an experimental cross between Chinese Meishan pigs and five commercial Dutch pig lines was used. Linkage analyses were performed using interval mapping by regression under line cross models including a test for imprinting effects. The whole-genome scan revealed highly significant evidence for three quantitative trait loci (QTL) affecting teat number, of which two were imprinted. Paternally expressed (i.e., maternally imprinted) QTL were found on chromosomes 2 and 12. A Mendelian expressed QTL was found on chromosome 10. The estimated additive effects showed that, for the QTL on chromosomes 10 and 12, the Meishan allele had a positive effect on teat number, but, for the QTL on chromosome 2, the Meishan allele had a negative effect on teat number. This study shows that imprinting may play an important role in the expression of teat number.

A high-resolution comparative RH map of porcine chromosome (SSC) 2

A high-resolution comparative map was constructed for porcine Chromosome (SSC) 2, where a QTL for back fat thickness (BFT) is located. A radiation hybrid (RH) map containing 33 genes and 25 microsatellite markers was constructed for this chromosome with a 3000-rad porcine RH panel. In total, 16 genes from human Chromosome (HSA) 11p, HSA19p, and HSA5q were newly assigned to SSC2. One linkage group was observed at LOD 3.0, and five linkage groups at LOD 4.0. Comparison of the porcine RH map with homologous human gene orders identified four conserved segments between SSC2 and HSA11, HSA19, and HSA5. Concerning HSA11, a rearrangement of gene order is observed. The segment HSA11p15.4-q13 is inverted on SSC2 when compared with the distal tip of SSC2p, which is homologous to HSA11p15.5. The boundaries of the conserved segments between human and pig were defined more precisely. This high-resolution comparative map will be a valuable tool for further fine mapping of the QTL area.