BackgroundThe current availability of genotypes for very large numbers of single nucleotide polymorphisms (SNPs) is leading to more accurate estimates of inbreeding coefficients and more detailed approaches for detecting inbreeding depression. In the present study, genome-wide information was used to detect inbreeding depression for two reproductive traits (total number of piglets born and number of piglets born alive) in an ancient strain of Iberian pigs (the Guadyerbas strain) that is currently under serious danger of extinction.MethodsA total of 109 sows with phenotypic records were genotyped with the PorcineSNP60 BeadChip v1. Inbreeding depression was estimated using a bivariate animal model in which the inbreeding coefficient was included as a covariate. We used two different measures of genomic inbreeding to perform the analyses: inbreeding estimated on a SNP-by-SNP basis and inbreeding estimated from runs of homozygosity. We also performed the analyses using pedigree-based inbreeding.ResultsSignificant inbreeding depression was detected for both traits using all three measures of inbreeding. Genome-wide information allowed us to identify one region on chromosome 13 associated with inbreeding depression. This region spans from 27 to 54 Mb and overlaps with a previously detected quantitative trait locus and includes the inter-alpha-trypsin inhibitor gene cluster that is involved with embryo implantation.ConclusionsOur results highlight the value of high-density SNP genotyping for providing new insights on where genes causing inbreeding depression are located in the genome. Genomic measures of inbreeding obtained on a SNP-by-SNP basis or those based on the presence/absence of runs of homozygosity represent a suitable alternative to pedigree-based measures to detect inbreeding depression, and a useful tool for mapping studies. To our knowledge, this is the first study in domesticated animals using the SNP-by-SNP inbreeding coefficient to map specific regions within chromosomes associated with inbreeding depression.Electronic supplementary materialThe online version of this article (doi:10.1186/s12711-014-0081-5) contains supplementary material, which is available to authorized users.

BackgroundRecent studies in pigs have detected copy number variants (CNVs) using the Comparative Genomic Hybridization technique in arrays designed to cover specific porcine chromosomes. The goal of this study was to identify CNV regions (CNVRs) in swine species based on whole genome SNP genotyping chips.ResultsWe used predictions from three different programs (cnvPartition, PennCNV and GADA) to analyze data from the Porcine SNP60 BeadChip. A total of 49 CNVRs were identified in 55 animals from an Iberian x Landrace cross (IBMAP) according to three criteria: detected in at least two animals, contained three or more consecutive SNPs and recalled by at least two programs. Mendelian inheritance of CNVRs was confirmed in animals belonging to several generations of the IBMAP cross. Subsequently, a segregation analysis of these CNVRs was performed in 372 additional animals from the IBMAP cross and its distribution was studied in 133 unrelated pig samples from different geographical origins. Five out of seven analyzed CNVRs were validated by real time quantitative PCR, some of which coincide with well known examples of CNVs conserved across mammalian species.ConclusionsOur results illustrate the usefulness of Porcine SNP60 BeadChip to detect CNVRs and show that structural variants can not be neglected when studying the genetic variability in this species.

BackgroundNew advances in high-throughput technologies have allowed for the massive analysis of genomic data, providing new opportunities for the characterization of the transcriptome architectures. Recent studies in pigs have employed RNA-Seq to explore the transcriptome of different tissues in a reduced number of animals. The main goal of this study was the identification of differentially-expressed genes in the liver of Iberian x Landrace crossbred pigs showing extreme phenotypes for intramuscular fatty acid composition using RNA-Seq.ResultsThe liver transcriptomes of two female groups (H and L) with phenotypically extreme intramuscular fatty acid composition were sequenced using RNA-Seq. A total of 146 and 180 unannotated protein-coding genes were identified in intergenic regions for the L and H groups, respectively. In addition, a range of 5.8 to 7.3% of repetitive elements was found, with SINEs being the most abundant elements. The expression in liver of 186 (L) and 270 (H) lncRNAs was also detected. The higher reproducibility of the RNA-Seq data was validated by RT-qPCR and porcine expression microarrays, therefore showing a strong correlation between RT-qPCR and RNA-Seq data (ranking from 0.79 to 0.96), as well as between microarrays and RNA-Seq (r=0.72). A differential expression analysis between H and L animals identified 55 genes differentially-expressed between groups. Pathways analysis revealed that these genes belong to biological functions, canonical pathways and three gene networks related to lipid and fatty acid metabolism. In concordance with the phenotypic classification, the pathways analysis inferred that linolenic and arachidonic acids metabolism was altered between extreme individuals. In addition, a connection was observed among the top three networks, hence suggesting that these genes are interconnected and play an important role in lipid and fatty acid metabolism.ConclusionsIn the present study RNA-Seq was used as a tool to explore the liver transcriptome of pigs with extreme phenotypes for intramuscular fatty acid composition. The differential gene expression analysis showed potential gene networks which affect lipid and fatty acid metabolism. These results may help in the design of selection strategies to improve the sensorial and nutritional quality of pork meat.

In addition to its value in meat production, the pig is an interesting animal model for human digestive tract studies due to its physiological similarities. The aim of this study was to describe the microbiome composition, distribution and interaction along the Iberian pig intestinal tract and its role in whole-body energy homeostasis. The V3-V4 region of the 16S rRNA gene was amplified and sequenced from the microbiomes of five gut sections (duodenum, jejunum, ileum, and proximal and distal colon) in thirteen castrated male pigs. A total of 1,669 operational taxonomic units distributed in 179 genera were found among all samples. The two most abundant genera in the small intestine were Lactobacillus and Clostridium, while Prevotella was predominant in the colon. The colon samples were more similar among the pigs and richer in species than the small intestine samples were. In the small intestine, the metagenome prediction pointed to rapid internalization and conversion of the available simple carbohydrates for microbial proliferation and maintenance. In the colon, a competition among anaerobic bacteria for plant polysaccharide degradation to produce short chain fatty acids was found. This study confirms that the energy pathways of the gut microbiome differ along its sections and provides a description of the correlations between genera.

The leptin receptor gene (LEPR) is a candidate for traits related to growth and body composition, and is located on SSC6 in a region where fatness and meat composition quantitative trait loci (QTL) have previously been detected in several F2 experimental designs. The aims of this work were: (i) to fine map these QTL on a larger sample of animals and generations (F3 and backcross) of an Iberian x Landrace intercross and (ii) to examine the effects of LEPR alleles on body composition traits. Eleven single nucleotide polymorphisms (SNPs) were detected by sequencing LEPR coding regions in Iberian and Landrace pig samples. Three missense polymorphisms were genotyped by pyrosequencing in 33 F0, 70 F1, 418 F2, 86 F3 and 128 individuals coming from the backcross of four F2 males with 24 Landrace females. Thirteen microsatellites and one SNP were also genotyped. Traits analysed were: backfat thickness at different locations (BF(T)), intramuscular fat percentage (IMF(P)), eye muscle area (EM(A)), loin depth (LO(D)), weight of shoulder (SH(W)), weight of ribs (RIB(W)) and weight of belly bacon (BB(W)). Different statistical models were applied in order to evaluate the number and effects of QTL on chromosome 6 and the possible causality of the LEPR gene variants with respect to the QTL. The results support the presence of two QTL on SSC6. One, at position 60-100 cM, affects BF(T) and RIB(W). The other and more significant maps in a narrow region (130-132 cM) and affects BF(T), IMF(P), EM(A), LO(D), SH(W), RIB(W) and BB(W). Results also support the association between LEPR alleles and BF(T) traits. The possible functional implications of the analysed polymorphisms are considered.