Quantitative trait loci for litter size and prenatal loss in a White Duroc × Chinese Erhualian resource population

Genomic evidence for human-mediated introgressive hybridization and selection in the developed breed

BACKGROUND

The pig (Sus Scrofa) is one of the oldest domesticated livestock species that has undergone extensive improvement through modern breeding. European breeds have advantages in lean meat development and highly-productive body type, whereas Asian breeds possess extraordinary fat deposition and reproductive performance. Consequently, Eurasian breeds have been extensively used to develop modern commercial breeds for fast-growing and high prolificacy. However, limited by the sequencing technology, the genome architecture of some nascent developed breeds and the human-mediated impact on their genomes are still unknown.

RESULTS

Through whole-genome analysis of 178 individuals from an Asian locally developed pig breed, Beijing Black pig, and its two ancestors from two different continents, we found the pervasive inconsistent gene trees and species trees across the genome of Beijing Black pig, which suggests its introgressive hybrid origin. Interestingly, we discovered that this developed breed has more genetic relationships with European pigs and an unexpected introgression from Asian pigs to this breed, which indicated that human-mediated introgression could form the porcine genome architecture in a completely different type compared to native introgression. We identified 554 genomic regions occupied 63.30 Mb with signals of introgression from the Asian ancestry to Beijing Black pig, and the genes in these regions enriched in pathways associated with meat quality, fertility, and disease-resistant. Additionally, a proportion of 7.77% of genomic regions were recognized as regions that have been under selection. Moreover, combined with the results of a genome-wide association study for meat quality traits in the 1537 Beijing Black pig population, two important candidate genes related to meat quality traits were identified. DNAJC6 is related to intramuscular fat content and fat deposition, and RUFY4 is related to meat pH and tenderness.

CONCLUSIONS

Our research provides insight for analyzing the origins of nascent developed breeds and genome-wide selection remaining in the developed breeds mediated by humans during modern breeding.

Identification of Signatures of Positive Selection That Have Shaped the Genomic Landscape of South African Pig Populations

South Africa boasts a diverse range of pig populations, encompassing intensively raised commercial breeds, as well as indigenous and village pigs reared under low-input production systems. The aim of this study was to investigate how natural and artificial selection have shaped the genomic landscape of South African pig populations sampled from different genetic backgrounds and production systems. For this purpose, the integrated haplotype score (iHS), as well as cross population extended haplotype homozygosity (XP-EHH) and Lewontin and Krakauer's extension of the Fst statistic based on haplotype information (HapFLK) were utilised. Our results revealed several population-specific signatures of selection associated with the different production systems. The importance of natural selection in village populations was highlighted, as the majority of genomic regions under selection were identified in these populations. Regions under natural and artificial selection causing the distinct genetic footprints of these populations also allow for the identification of genes and pathways that may influence production and adaptation. In the context of intensively raised commercial pig breeds (Large White, Kolbroek, and Windsnyer), the identified regions included quantitative loci (QTLs) associated with economically important traits. For example, meat and carcass QTLs were prevalent in all the populations, showing the potential of village and indigenous populations' ability to be managed and improved for such traits. Results of this study therefore increase our understanding of the intricate interplay between selection pressures, genomic adaptations, and desirable traits within South African pig populations.

Effect of NLR family pyrin domain containing 9 gene polymorphism on litter size in large white pigs

NLR family pyrin domain containing 9 (NLRP9) is a mammalian reproduction-related gene. In this study, we researched the associations between polymorphisms located in the coding sequence (CDS) of the NLRP9 gene, and both the total number of piglets born per litter (TNB) and the number of piglets born alive per litter (NBA) in Canada Large White pigs (CLW). We detected a single nucleotide polymorphism (SNP) within exon 3 (g.10910C > T). The allele frequencies at the NLRP9 locus were 0.474 for the C allele and 0.526 for the T allele. Three genotypes, CC, CT, and TT, occurred with frequencies of 0.216, 0.515, and 0.269, respectively. Sows with the CC genotype had the largest TNB and NBA, sows with TT had the smallest, and those with CT were in-between. This difference was statistically significant (p < 0.05). Furthermore, CC females grew faster than CT or TT females, and there was a significant relationship between NLRP9 polymorphism and the average daily gain (p < 0.05). Here, we provide the first evidence for a novel SNP in NLRP9 associated with litter size in CLW sows, which could be used as a genetic marker to improve litter size in pig breeding and production.

CYP24A1 is associated with fetal mummification in pigs

Mummified piglets are among the leading causes of fertility loss and severely hamper reproductive performance in pigs. However, the contributions of genomic variation to the emergence of mummified piglets (MUM) have rarely been studied. This study aims to (1) elucidate the genetic architecture of MUM in sows of parity 1 - 3 using a single-step genome-wide association study (ssGWAS). The ssGWAS involved genotyping-by-sequencing of Large White and Landrace pig breeds. (2) Explore the biological role of the candidate genes at the cellular level. A total of 185 and 48 genome-wide significant SNPs are associated with MUM in Large White and Landrace pigs, explaining 0.01-36.52% genetic variance for different significant loci, respectively. All the significant SNPs are parity-specific, and the numerous, consecutive significant loci likely generated the nine significant peaks in different parities. Multiple candidate genes (including CYP24A1, FBXO30, and ARHGEF28) are associated with fetal congenital and maternal diseases. Collectively, CYP24A1 regulation contributes to steady-state levels of embryo development genes. CYP24A1 is involved in reproduction and, immune and gestational disorders. Thus, it is associated with known newborn death traits and MUM in Large White sows. Altogether, these results improve the current understanding of the genetic architecture of MUM and expand the knowledge on genetic variations for selecting against mummified piglets in pig breeding.

Retrotransposon Insertion Polymorphisms (RIPs) in Pig Reproductive Candidate Genes

Retrotransposons account for more than one-third of the pig reference genome. On account of the genome variability in different breeds, structural variation (SV) caused by retrotranspos-on-generated deletion or insertion (indel) may have a function in the genome. Litter size is one of the most important reproductive traits and significantly impacts profitability in terms of pig production. We used the method of bioinformatics, genetics, and molecular biology to make an analysis among different pig genomes. Predicted 100 SVs were annotated as retrotransposon indel in 20 genes related to reproductive performance. The PCR detection based on these predicted SVs revealed 20 RIPs in 20 genes, that most RIPs (12) were generated by SINE indel, and eight RIPs were generated by the ERV indel. We selected 12 RIPs to make the second round PCR detection in 24 individuals among nine pig breeds. The PCR detection results revealed that the RIP-A1CF-4 insertion in the breed of Bama, Large White, and Meishan only had the homozygous genotype but low to moderately polymorphisms were present in other breeds. We found that RIP-CWH43-9, RIP-IDO2-9, RIP-PRLR-6, RIP-VMP1-12, and RIP-OPN-1 had a rich polymorphism in the breed of Large White pigs. The statistical analysis revealed that RIP-CWH43-9 had a SINE insertion profitable to the reproductive traits of TNB and NBA but was significantly affected (p < 0.01) and (p < 0.05) in the reproductive traits of litter birthweight (LW) in Large White. On the other hand, the SINE insertion in IDO2-9 may be a disadvantage to the reproductive traits of LW, which was significantly affected (p < 0.05) in Large White. These two RIPs are significant in pig genome research and could be useful molecular markers in the breeding system.

Identification of Genomic Regions and Candidate Genes for Litter Traits in French Large White Pigs Using Genome-Wide Association Studies

The reproductive traits of sows are one of the important economic traits in pig production, and their performance directly affects the economic benefits of the entire pig industry. In this study, a total of 895 French Large White pigs were genotyped by GeneSeek Porcine 50K SNP Beadchip and four phenotypic traits of 1407 pigs were recorded, including total number born (TNB), number born alive (NBA), number healthy piglets (NHP) and litter weight born alive (LWB). To identify genomic regions and genes for these traits, we used two approaches: a single-locus genome-wide association study (GWAS) and a single-step GWAS (ssGWAS). Overall, a total of five SNPs and 36 genomic regions were identified by single-locus GWAS and ssGWAS, respectively. Notably, fourof all five significant SNPs were located in 10.72–11.06 Mb on chromosome 7, were also identified by ssGWAS. These regions explained the highest or second highest genetic variance in the TNB, NBA and NHP traits and harbor the protein coding gene ENSSSCG00000042180. In addition, several candidate genes associated with litter traits were identified, including JARID2, PDIA6, FLRT2 and DICER1. Overall, these novel results reflect the polygenic genetic architecture of the litter traits and provide a theoretical reference for the following implementation of molecular breeding.

Genetic architecture underlying nascent speciation - The evolution of Eurasian pigs under domestication

Speciation is a process whereby the evolution of reproductive barriers leads to isolated species. Although many studies have addressed large-effect genetic footprints in the advanced stages of speciation, the genetics of reproductive isolation in nascent stage of speciation remains unclear. Here, we show that pig domestication offers an interesting model for studying the early stages of speciation in great details. Pig breeds have not evolved the large X-effect of hybrid incompatibility commonly observed between “good species.” Instead, deleterious epistatic interactions among multiple autosomal loci are common. These weak Dobzhansky–Muller incompatibilities confer partial hybrid inviability with sex biases in crosses between European and East Asian domestic pigs. The genomic incompatibility is enriched in pathways for angiogenesis, androgen receptor signaling and immunity, with an observation of many highly differentiated cis-regulatory variants. Our study suggests that partial hybrid inviability caused by pervasive but weak interactions among autosomal loci may be a hallmark of nascent speciation in mammals.

Integrated analysis of lncRNA, miRNA and mRNA reveals novel insights into the fertility regulation of large white sows

Background

Improving sow fertility is extremely important as it can lead to increased reproductive efficiency and thus profitability for swine producers. There are considerable differences in fertility rates among individual animals, but the underlying molecular mechanisms remain unclear. In this study, by using different types of RNA libraries, we investigated the complete transcriptome of ovarian tissue during the luteal (L) and follicular (F) phases of the estrous cycle in Large White pigs with high (H) and low (L) fecundity, and performed a comprehensive analysis of long noncoding RNAs (lncRNAs), mRNAs and micro RNAs (miRNAs) from 16 samples by combining RNA sequencing (RNA-seq) with bioinformatics.

Results

In total, 24,447 lncRNAs, 27,370 mRNAs, and 216 known miRNAs were identified in ovarian tissues. The genomic features of lncRNAs, such as length distribution and number of exons, were further analyzed. We selected a threshold of P < 0.05 and |log₂ (fold change)| ≥ 1 to obtain the differentially expressed lncRNAs, miRNAs and mRNAs by pairwise comparison (LH vs. LL, FH vs. FL). Bioinformatics analysis of these differentially expressed RNAs revealed multiple significantly enriched pathways (P < 0.05) that were closely involved in the reproductive process, such as ovarian steroidogenesis, lysosome, steroid biosynthesis, and the estrogen and GnRH signaling pathways. Moreover, bioinformatics screening of differentially expressed miRNAs that share common miRNA response elements (MREs) with lncRNAs and their downstream mRNA targets were performed. Finally, we constructed lncRNA–miRNA–mRNA regulation networks. The key genes in these networks were verified by Reverse Transcription Real-time Quantitative PCR (RT-qRCR), which were consistent with the results from RNA-Seq data.

Conclusions

These results provide further insights into the fertility of pigs andcan contribute to further experimental investigation of the functions of these genes.

SLA-11 mutations are associated with litter size traits in Large White and Chinese DIV pigs

Litter size is an important economic traits in pigs. SLA-11 gene is a member of SLA (swine leukocyte antigen) complex. In our previous study, the SLA-11 gene was differentially expressed in PMSG-hCG stimulated preovulatory ovarian follicles of Chinese Taihu and Large White sows. Here, we identified two mutations (c.754-132 T > C and c.1421 + 38 T > C) in SLA-11 gene and analyzed the associations of two SNPs with litter size traits in Large White (n = 263) and DIV (n = 117) sows. The results showed that in Large White pigs, SLA-11 c.754-132 CC sows produced 0.74 and 0.87 more pigs per litter for TNB and NBA of all parities than did TT sows (p < .05); In DIV pigs, SLA-11 c.754-132 CC sows produced 1.17 more pigs per litter for TNB of all parities than did TC sows (p < .05). In Large White pigs, SLA-11 c.1421 + 38 CC sows produced 0.9 more pigs per litter for TNB of all parities than did TT sows (p < .05), while in DIV pigs SLA-11 c.1421 + 38 CC sows produced 0.84 and 0.7 less pigs per litter for TNB and NBA of all parities than did TT sows (p < .05). Our research indicated that SLA-11 mutations were potential molecular markers for improving the litter size traits in pigs.

Genome-wide genetic structure and differentially selected regions among Landrace, Erhualian, and Meishan pigs using specific-locus amplified fragment sequencing

As typical Chinese indigenous pig breeds, Erhualian and Meishan have been widely used to produce new strain or breed in the world. However, the genetic basis of characteristics of these pig breeds is still limited. Moreover, considering cost and output of sequencing, it is necessary to further develop cost-effective method for pig genome screening. To contribute on this issue, we developed a SLAF-seq (specific-locus amplified fragment sequencing) method for pigs and applied it to analyze the genetic difference among Landrace, Erhualian, and Meishan pigs. A total of 453.75 million reads were produced by SLAF-seq. After quality-control, 165,670 SNPs (single nucleotide polymorphisms) were used in further analysis. The results showed that Landrace had distinct genetic relationship compared to Erhualian (F_ST = 0.5480) and Meishan (F_ST = 0.5800), respectively, while Erhualian and Meishan held the relatively close genetic relationship (F_ST = 0.2335). Furthermore, a genome-wide scanning revealed 268 differentially selected regions (DSRs) with 855 genes and 256 DSRs with 347 genes between Landrace and the two Chinese indigenous pig breeds and between Erhualian and Meishan, respectively. This study provides a new cost-effective method for pig genome study and might contribute to a better understanding on the formation mechanism of genetic difference among pigs with different geographical origins.

Genetic analysis of litter size in mice

We performed quantitative trait locus (QTL) mapping analysis for litter size (total number of pups born and/or number of pups born alive) in 255 backcross mice derived from C57BL/6J and RR/Sgn inbred mice. We identified one significant QTL on chromosome 7 and 4 suggestive QTLs on chromosomes 3, 5, 10 and 13. In addition, two suggestive QTLs were identified on chromosomes 1 and 4 for the number of stillbirth. These results suggested that both litter size and number of stillbirth were heritable traits, although they were controlled by distinct genes. The RR allele was associated with reduced litter size and increased stillbirth at all QTLs. Therefore, RR mothers were observed to have reduced prolificacy in this particular genetic cross.

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.

A genome-wide linkage analysis for reproductive traits in F2 Large White × Meishan cross gilts

Female reproductive performance traits in pigs have low heritabilities thus limiting improvement through traditional selective breeding programmes. However, there is substantial genetic variation found between pig breeds with the Chinese Meishan being one of the most prolific pig breeds known. In this study, three cohorts of Large White × Meishan F2 cross-bred pigs were analysed to identify quantitative trait loci (QTL) with effects on reproductive traits, including ovulation rate, teat number, litter size, total born alive and prenatal survival. A total of 307 individuals were genotyped for 174 genetic markers across the genome. The genome-wide analysis of the trait-recorded F2 gilts in their first parity/litter revealed one QTL for teat number significant at the genome level and a total of 12 QTL, which are significant at the chromosome-wide level, for: litter size (three QTL), total born alive (two QTL), ovulation rate (four QTL), prenatal survival (one QTL) and teat number (two QTL). Further support for eight of these QTL is provided by results from other studies. Four of these 12 QTL were mapped for the first time in this study: on SSC15 for ovulation rate and on SSC18 for teat number, ovulation rate and litter size.

A comprehensive survey of copy number variation in 18 diverse pig populations and identification of candidate copy number variable genes associated with complex traits

Background

Copy number variation (CNV) is a major source of structural variants and has been commonly identified in mammalian genome. It is associated with gene expression and may present a major genetic component of phenotypic diversity. Unlike many other mammalian genomes where CNVs have been well annotated, studies of porcine CNV in diverse breeds are still limited.

Result

Here we used Porcine SNP60 BeadChip and PennCNV algorithm to identify 1,315 putative CNVs belonging to 565 CNV regions (CNVRs) in 1,693 pigs from 18 diverse populations. Total 538 out of 683 CNVs identified in a White Duroc × Erhualian F₂ population fit Mendelian transmission and 6 out of 7 randomly selected CNVRs were confirmed by quantitative real time PCR. CNVRs were non-randomly distributed in the pig genome. Several CNV hotspots were found on pig chromosomes 6, 11, 13, 14 and 17. CNV numbers differ greatly among different pig populations. The Duroc pigs were identified to have the most number of CNVs per individual. Among 1,765 transcripts located within the CNVRs, 634 genes have been reported to be copy number variable genes in the human genome. By integrating analysis of QTL mapping, CNVRs and the description of phenotypes in knockout mice, we identified 7 copy number variable genes as candidate genes for phenotypes related to carcass length, backfat thickness, abdominal fat weight, length of scapular, intermuscle fat content of logissimus muscle, body weight at 240 day, glycolytic potential of logissimus muscle, mean corpuscular hemoglobin, mean corpuscular volume and humerus diameter.

Conclusion

We revealed the distribution of the unprecedented number of 565 CNVRs in pig genome and investigated copy number variable genes as the possible candidate genes for phenotypic traits. These findings give novel insights into porcine CNVs and provide resources to facilitate the identification of trait-related CNVs.

Identification of the promoter region and genetic mutations of the porcine GALP gene

Galanin-like peptide (GALP) gene, encoding a member of the galanin family of neuropeptides involved in reproduction, was differentially expressed in PMSG-hCG stimulated pre-ovulatory ovarian follicles of Chinese Taihu and Large White sows in our previous study. In the present study, promoter region and genetic mutations of the porcine GALP gene were determined. A 1,322 bp contig in 5'-flanking region was predicted to contain 5 potential transcription promoters by Neural Network Promoter Prediction version 2.2. 5'-deletion expression in both CHO and hela cells showed that there were a negative regulatory element at -852 to -803 bp and a positive regulatory element at -1,318 to -1,269 bp. Comparative sequence analyses of Chinese Taihu and Large White GALP gene sequence revealed the c.*27C>G mutation in the 3'-UTR and the c.88-1225C>G mutation in intron 1, which can be detected by HhaI and AluI PCR-RFLP, respectively. The association analysis with litter size traits showed that at both loci CC and GG genotypes were different for NBA for all parities in DIV pigs (P < 0.05). However, two SNPs were not in significant linkage disequilibrium analyzed using SHEsis online software, and could be used in pig breeding individually.

A whole genome scan to detect quantitative trait loci for gestation length and sow maternal ability related traits in a White Duroc × Erhualian F2 resource population

Gestation length and maternal ability are important to improve the sow reproduction efficiency and their offspring survival. To map quantitative trait loci (QTL) for gestation length and maternal ability related traits including piglet survival rate and average body weight of piglets at weaning, more than 200 F2 sows from a White Duroc × Erhualian resource population were phenotyped. A genome-wide scan was performed with 194 microsatellite markers covering the whole pig genome. QTL analysis was carried out using a composite regression interval mapping method via QTL express. The results showed that total number of born piglets was significantly correlated with gestation length (r = -0.13, P < 0.05). Three QTL were detected on pig chromosome (SSC)2, 8 and 12 for gestation length. The QTL on SSC2 achieved the 5% genome-wide significant level and the QTL on SSC8 was consistent with previous reports. Four suggestive QTL were identified for maternal ability related traits including 1 QTL for survival rate of piglets at weaning on SSC8, 3 QTL for average body weight of piglet at weaning on SSC3, 11 and 13.

Whole-genome SNP association analysis of reproduction traits in the Finnish Landrace pig breed

Background

Good genetic progress for pig reproduction traits has been achieved using a quantitative genetics-based multi-trait BLUP evaluation system. At present, whole-genome single nucleotide polymorphisms (SNP) panels provide a new tool for pig selection. The purpose of this study was to identify SNP associated with reproduction traits in the Finnish Landrace pig breed using the Illumina PorcineSNP60 BeadChip.

Methods

Association of each SNP with different traits was tested with a weighted linear model, using SNP genotype as a covariate and animal as a random variable. Deregressed estimated breeding values of the progeny tested boars were used as the dependent variable and weights were based on their reliabilities. Statistical significance of the associations was based on Bonferroni-corrected P-values.

Results

Deregressed estimated breeding values were available for 328 genotyped boars. Of the 62 163 SNP in the chip, 57 868 SNP had a call rate > 0.9 and 7 632 SNP were monomorphic. Statistically significant results (P-value < 2.0E-06) were obtained for total number of piglets born in first and later parities and piglet mortality between birth and weaning in later parity, and suggestive associations (P-value < 4.0E-06) for piglet mortality between birth and weaning in first parity, number of stillborn piglets in later parity, first farrowing interval and second farrowing interval. Two of the statistically significant regions for total number of piglets born in first and later parities are located on chromosome 9 around 95 and 79 Mb. The estimated SNP effect in these regions was approximately one piglet between the two homozygote classes. By combining the two most significant SNP in these regions, favourable double homozygote animals are expected to have 1.3 piglets (P-value = 1.69E-08) more than unfavourable double homozygote animals. A region on chromosome 9 (66 Mb) was statistically significant for piglet mortality between birth and weaning in later parity (0.44 piglets between homozygotes, P-value = 6.94E-08).

Conclusions

Three separate regions on chromosome 9 gave significant results for litter size and pig mortality. The frequencies of favourable alleles of the significant SNP are moderate in the Finnish Landrace population and these SNP are thus valuable candidates for possible marker-assisted selection.

Microsatellite mapping of quantitative trait loci affecting female reproductive tract characteristics in Meishan x Large White F(2) pigs

A QTL analysis of female reproductive data from a 3-generation experimental cross between Meishan and Large White pig breeds is presented. Six F(1) boars and 23 F(1) sows, progeny of 6 Large White boars and 6 Meishan sows, produced 502 F(2) gilts whose reproductive tract was collected after slaughter at 30 d of gestation. Five traits [i.e., the total weight of the reproductive tract, of the empty uterine horns, of the ovaries (WOV), and of the embryos], as well as the length of uterine horns (LUH), were measured and analyzed with and without adjustment for litter size. Animals were genotyped for a total of 137 markers covering the entire porcine genome. Analyses were carried out based on interval mapping methods, using a line-cross regression and a half-full sib maximum likelihood test. A total of 18 genome-wide significant (P < 0.05) QTL were detected on 9 different chromosomes (i.e., SSC 1, 5, 6, 7, 9, 12, 13, 18, and X). Five genome-wide significant QTL were detected for LUH, 4 for weight of the empty uterine horns and WOV, 2 for total weight of the reproductive tract, and 1 for weight of the embryos. Twenty-two additional suggestive QTL were also detected. The largest effects were obtained for LUH and WOV on SSC13 (9.2 and 7.0% of trait phenotypic variance, respectively). Meishan alleles had both positive (e.g., on SSC7) and negative effects (e.g., on SSC13) on the traits investigated. Moreover, the QTL were generally not fixed in founder breeds, and opposite effects were in some cases obtained in different families. Although reproductive tract characteristics had only a moderate correlation with reproductive performances, most of the major QTL detected in this study were previously reported as affecting female reproduction, generally with reduced significance levels. This study thus shows that focusing on traits with high heritability might help to detect loci involved in low heritability major traits for breeding.

A bi-dimensional genome scan for prolificacy traits in pigs shows the existence of multiple epistatic QTL

BACKGROUND

Prolificacy is the most important trait influencing the reproductive efficiency of pig production systems. The low heritability and sex-limited expression of prolificacy have hindered to some extent the improvement of this trait through artificial selection. Moreover, the relative contributions of additive, dominant and epistatic QTL to the genetic variance of pig prolificacy remain to be defined. In this work, we have undertaken this issue by performing one-dimensional and bi-dimensional genome scans for number of piglets born alive (NBA) and total number of piglets born (TNB) in a three generation Iberian by Meishan F(2) intercross.

RESULTS

The one-dimensional genome scan for NBA and TNB revealed the existence of two genome-wide highly significant QTL located on SSC13 (P < 0.001) and SSC17 (P < 0.01) with effects on both traits. This relative paucity of significant results contrasted very strongly with the wide array of highly significant epistatic QTL that emerged in the bi-dimensional genome-wide scan analysis. As much as 18 epistatic QTL were found for NBA (four at P < 0.01 and five at P < 0.05) and TNB (three at P < 0.01 and six at P < 0.05), respectively. These epistatic QTL were distributed in multiple genomic regions, which covered 13 of the 18 pig autosomes, and they had small individual effects that ranged between 3 to 4% of the phenotypic variance. Different patterns of interactions (a x a, a x d, d x a and d x d) were found amongst the epistatic QTL pairs identified in the current work.

CONCLUSIONS

The complex inheritance of prolificacy traits in pigs has been evidenced by identifying multiple additive (SSC13 and SSC17), dominant and epistatic QTL in an Iberian x Meishan F(2) intercross. Our results demonstrate that a significant fraction of the phenotypic variance of swine prolificacy traits can be attributed to first-order gene-by-gene interactions emphasizing that the phenotypic effects of alleles might be strongly modulated by the genetic background where they segregate.