Search

A better understanding of the environmental and host factors modulating gut microbiomes is a topic of greatest interest. Recent evidence suggests that genetic variation in the pig genome partially controls the composition of porcine gut microbiota.

Background The gut microbiota influences host performance playing a relevant role in homeostasis and function of the immune system. The aim of the present work was to identify microbial signatures linked to immunity traits and to characterize the contribution of host-genome and gut microbiota to the immunocompetence in healthy pigs. Results To achieve this goal, we undertook a combination of network, mixed model and microbial-wide association studies (MWAS) for 21 immunity traits and the relative abundance of gut bacterial communities in 389 pigs genotyped for 70K SNPs. The heritability (h2; proportion of phenotypic variance explained by the host genetics) and microbiability (m2; proportion of variance explained by the microbial composition) showed similar values for most of the analyzed immunity traits, except for both IgM and IgG in plasma that was dominated by the host genetics, and the haptoglobin in serum which was the trait with larger m2 (0.275) compared to h2 (0.138). Results from the MWAS suggested a polymicrobial nature of the immunocompetence in pigs and revealed associations between pigs gut microbiota composition and 15 of the analyzed traits. The lymphocytes phagocytic capacity (quantified as mean fluorescence) and the total number of monocytes in blood were the traits associated with the largest number of taxa (6 taxa). Among the associations identified by MWAS, 30% were confirmed by an information theory network approach. The strongest confirmed associations were between Fibrobacter and phagocytic capacity of lymphocytes (r = 0.37), followed by correlations between Streptococcus and the percentage of phagocytic lymphocytes (r = -0.34) and between Megasphaera and serum concentration of haptoglobin (r = 0.26). In the interaction network, Streptococcus and percentage of phagocytic lymphocytes were the keystone bacterial and immune-trait, respectively. Conclusions Overall, our findings reveal an important connection between gut microbiota composition and immunity traits in pigs, and highlight the need to consider both sources of information, host genome and microbial levels, to accurately characterize immunocompetence in pigs.

In recent years, the increase in awareness of antimicrobial resistance together with the societal demand of healthier meat products have driven attention to health-related traits in livestock production. Previous studies have reported medium to high heritabilities for these traits and described genomic regions associated with them. Despite its genetic component, health- and immunity-related traits are complex and its study by association analysis with genomic markers may be missing some information. To analyse multiple phenotypes and gene-by-gene interactions, systems biology approaches, such as the association weight matrix (AWM), allows combining genome wide association study results with network inference algorithms. The present study aimed to identify gene networks, key regulators and candidate genes associated to immunocompetence in pigs by integrating multiple health-related traits, enriched for innate immune phenotypes, using the AWM approach. The co-association network analysis unveiled a network comprised of 3,636 nodes (genes) and 451,407 edges (interactions), including a total of 246 regulators. From these, five genes (ARNT2,BRMS1L,MED12L,SUPT3HandTRIM25) were selected as key regulators as they were associated with the maximum number of genes with the minimum overlapping (1,827 genes in total). The five regulators were involved in pathways related to immunity such as lymphocyte differentiation and activation, platelet activation and degranulation, megakaryocyte differentiation, FcγR-mediated phagocytosis and response to nitric oxide, among others, but also in immunometabolism. Furthermore, we identified genes co-associated with the key regulators previously reported as candidate genes (e.g.,ANGPT1, CD4, CD36, DOCK1,PDE4B, PRKCE, PTPRCandSH2B3) for immunity traits in humans and pigs, but also new candidate ones (e.g.,ACSL3,CXADR,HBB,MMP12,PTPN6,WLS) that were not previously described. The co-association analysis revealed new regulators associated with health-related traits in pigs. This approach also identified gene-by-gene interactions and candidate genes involved in pathways related to cell fate and metabolic and immune functions. Our results shed new light in the regulatory mechanisms involved in pig immunity and reinforce the use of the pig as biomedical model.

[EN] Iberian pigs are meat quality models due to their high fat content, high intramuscular fat, and oleic fatty acid composition. These parameters present great variability and are differentiated among the lines that make up the Iberian pig population. However, there is little information on how the genetic expression influences quality across Iberian varieties. This study aimed to compare the muscle expression profile between two varieties of Iberian pig (Torbiscal and Retinto) and their reciprocal crosses, differentiated by fatness. Our results suggest that the Retinto variety, which has the greatest fat content amongst the studied Iberian varieties, showed a higher expression of genes related to adiposity. Likewise, a higher expression of genes related to lipolysis was found in the Torbiscal variety, described as having less fat content than Retinto. Further genetic variation analysis in these Iberian varieties showed relevant associations for SNP (Single Nucleotide Polymorphism), related to these differentially expressed genes, with the meat quality traits. Thus, our findings evidence that differences in the genetic architecture and expression of Iberian varieties might explain the variability in their fat content and composition and hence, their meat quality. The high deposition of intramuscular fat and the content of oleic fatty acid are characteristic of the Iberian pig. These two parameters present great variability and are differentiated amongst the varieties that make up the Iberian pig population. Although previous studies generated evidence for causal genes and polymorphisms associated to the adipogenic potential of the Iberian pig, there is little information about how genetic expression influences this trait's variability. The aim of this study was to analyses the expression profile between two varieties of Iberian pig (Torbiscal and Retinto) and their reciprocal crosses differentiated in their intramuscular fat (IMF) content and fatty acid (FA) composition in the Longissimus thoracis muscle using an RNA-seq approach. Our results corroborate that the Retinto variety is the fattiest amongst all studied varieties as its upregulated genes, such as FABP3 and FABP5, SLC27A1 and VEGFA among others, contribute to increasing adiposity. In its turn, Torbiscal pigs showed an upregulation of genes associated with the inhibition of fat deposition such as ADIPOQ and CPT1A. Further genetic variation analysis in these Iberian varieties showed relevant associations for SNP located within the differentially expressed genes with IMF and FA content. Thus, the differences found in the genetic architecture and the muscle transcriptome of these Iberian varieties might explain the variability in their fat content and composition and hence, their meat quality., This study was funded by the INIA grant RTA 2012-0054-C02-01* and the Centre for the Development of Industrial Technology grant IDI-20140447.

mir-33a and mir-33b are co-transcribed with the SREBF2 and SREBF1 transcription factors, respectively. The main role of SREBF1 is the regulation of genes involved in fatty acid metabolism, while SREBF2 regulates genes participating in cholesterol biosynthesis and uptake. Our objective was to study the expression of both miR-33a and miR-33b, together with their host SREBF genes, in liver, adipose tissue and muscle to better understand the role of miR-33a/b in the lipid metabolism of pigs. In our study, the expression of miR-33a, miR-33b and SREBF2 in liver, adipose tissue, and muscle was studied in 42 BC1_LD (25% Iberian x 75% Landrace backcross) pigs by RT-qPCR. In addition, the expression of in-silico predicted target genes and fatty acid composition traits were correlated with the miR-33a/b expression. We observed different tissue expression patterns for both miRNAs. In adipose tissue and muscle a high correlation between miR-33a and miR-33b expression was found, whereas a lower correlation was observed in liver. The expression analysis of in-silico predicted target-lipid related genes showed negative correlations between miR-33b and CPT1A expression in liver. Conversely, positive correlations between miR-33a and PPARGC1A and USF1 gene expression in liver were observed. Lastly, positive and negative correlations between miR-33a/b expression and saturated fatty acid (SFA) and polyunsaturated fatty acid (PUFA) content, respectively, were identified. Overall, our results suggested that both miRNAs are differentially regulated and have distinct functions in liver, in contrast to muscle and adipose tissue. Furthermore, the correlations between miR-33a/b expression both with the expression of in-silico predicted target-lipid related genes and with fatty acid composition, opens new avenues to explore the role of miR33a/b in the regulation of lipid metabolism., This work was supported by the Spanish Ministerio de Economía y Competitividad (MINECO) with project AGL2017-82641-R (https://portal.mineco.gob.es/es-es/Paginas/default.aspx) and the Fondo Europeo de Desarrollo Regional (FEDER, http://ec.europa.eu/regional_policy/es/funding/erdf/). LCM was financially supported by an FPI grant from the AGL2014-56369-C2 project. MB was funded with a “Ramón y Cajal” contract (RYC-2013-12573) from the Spanish Ministerio de Economía y Competitividad. DCP was funded by a “Formació i Contractació de Personal Investigador Novell” (FI-DGR) Ph.D grant from the Generalitat de Catalunya (ECO/1788/2014). MP was funded by an FPI grant from the AGL2017-82641-R project. We acknowledge the support of the Spanish Ministerio de Economía y Competitividad for the “Severo Ochoa Programme for Centres of Excellence in R&D” 2016-2019 (SEV-2015-0533) to the Centre for Research in Agricultural Genomics and the CERCA Programme / Generalitat de Catalunya.

[Background]: Analyses of gut microbiome composition in livestock species have shown its potential to contribute to the regulation of complex phenotypes. However, little is known about the host genetic control over the gut microbial communities. In pigs, previous studies are based on classical “single-gene-single-trait” approaches and have evaluated the role of host genome controlling gut prokaryote and eukaryote communities separately., [Results]: In order to determine the ability of the host genome to control the diversity and composition of microbial communities in healthy pigs, we undertook genome-wide association studies (GWAS) for 39 microbial phenotypes that included 2 diversity indexes, and the relative abundance of 31 bacterial and six commensal protist genera in 390 pigs genotyped for 70 K SNPs. The GWAS results were processed through a 3-step analytical pipeline comprised of (1) association weight matrix; (2) regulatory impact factor; and (3) partial correlation and information theory. The inferred gene regulatory network comprised 3561 genes (within a 5 kb distance from a relevant SNP–P, [Conclusions]: Taken together, our results identified regulators, candidate genes, and mechanisms linked with microbiome modulation by the host. They further highlight the value of the proposed analytical pipeline to exploit pleiotropy and the crosstalk between bacteria and protists as significant contributors to host-microbiome interactions and identify genetic markers and candidate genes that can be incorporated in breeding program to improve host-performance and microbial traits., YRC was funded by Marie Skłodowska-Curie grant (P-Sphere) agreement No. 6655919 (EU). MB is recipient of a Ramon y Cajal post-doctoral fellowship (RYC-2013–12573) from the Spanish Ministry of Economy and Competitiveness. EMS was funded with a PhD fellowship FPU15/01733 awarded by the Spanish Ministry of Education and Culture. Part of the research presented in this publication was funded by Grants AGL2016–75432-R, AGL2017–88849-R awarded by the Spanish Ministry of Economy and Competitiveness. The authors belong to Consolidated Research Group TERRA (AGAUR, 2017 SGR 1719).

The inclusion of health-related traits, or functionally associated genetic markers, in pig breeding programs could contribute to produce more robust and disease resistant animals. The aim of the present work was to study the genetic determinism and genomic regions associated to global immunocompetence and health in a Duroc pig population. For this purpose, a set of 30 health-related traits covering immune (mainly innate), haematological, and stress parameters were measured in 432 healthy Duroc piglets aged 8 weeks. Moderate to high heritabilities were obtained for most traits and significant genetic correlations among them were observed. A genome wide association study pointed out 31 significantly associated SNPs at whole-genome level, located in six chromosomal regions on pig chromosomes SSC4, SSC6, SSC17 and SSCX, for IgG, γδ T-cells, C-reactive protein, lymphocytes phagocytic capacity, total number of lymphocytes, mean corpuscular volume and mean corpuscular haemoglobin. A total of 16 promising functionally-related candidate genes, including CRP, NFATC2, PRDX1, SLA, ST3GAL1, and VPS4A, have been proposed to explain the variation of immune and haematological traits. Our results enhance the knowledge of the genetic control of traits related with immunity and support the possibility of applying effective selection programs to improve immunocompetence in pigs.

Meat quality has an important genetic component and can be modified by the fatty acid (FA) composition and the amount of fat contained in adipose tissue and muscle. The present study aimed to find genomic regions associated with the FA composition in backfat and muscle (longissimus dorsi) in 439 pigs with three different genetic backgrounds but having the Iberian breed in common. Genome-wide association studies (GWAS) were performed between 38,424 single-nucleotide polymorphisms (SNPs) covering the pig genome and 60 phenotypic traits related to backfat and muscle FA composition. Nine significant associated regions were found in backfat on the Sus scrofa chromosomes (SSC): SSC1, SSC2, SSC4, SSC6, SSC8, SSC10, SSC12, and SSC16. For the intramuscular fat, six significant associated regions were identified on SSC4, SSC13, SSC14, and SSC17. A total of 52 candidate genes were proposed to explain the variation in backfat and muscle FA composition traits. GWAS were also reanalysed including SNPs on five candidate genes (ELOVL6, ELOVL7, FADS2, FASN, and SCD). Regions and molecular markers described in our study may be useful for meat quality selection of commercial pig breeds, although several polymorphisms were breed-specific, and further analysis would be needed to evaluate possible causal mutations., This work was supported by the Spanish Ministerio de Economía y Competitividad (MINECO) and the Fondo Europeo de Desarrollo Regional (FEDER) with Project References: AGL2014-56369-C2 and AGL2017-82641-R. D. Crespo-Piazuelo was funded by a “Formació i Contractació de Personal Investigador Novell” (FI-DGR) Ph.D grant from the Generalitat de Catalunya (ECO/1788/2014). L. Criado-Mesas was funded with a FPI grant from the AGL2014-56369-C2 project. M. Revilla was also funded by a FI-DGR (ECO/1639/2013). M. Ballester was financially supported by a “Ramón y Cajal” Contract (RYC-2013-12573) from the Spanish Ministerio de Economía y Competitividad. We acknowledge the support of the Spanish Ministerio de Economía y Competitividad for the “Severo Ochoa Programme for Centres of Excellence in R&D” 2016–2019 (SEV-2015-0533) to the Centre for Research in Agricultural Genomics and the CERCA Programme / Generalitat de Catalunya.

Background The pig gut microbiome harbors thousands of species of archaea, bacteria, viruses and eukaryotes such as protists and fungi. However, since the majority of published studies have been focused on prokaryotes, little is known about the diversity, host-genetic control, and contributions to host performance of the gut eukaryotic counterparts. Here we report the first study that aims at characterizing the diversity and composition of gut commensal eukaryotes in pigs, exploring their putative control by host genetics, and analyzing their association with piglets body weight. Results Fungi and protists from the faeces of 514 healthy Duroc pigs of two sexes and two different ages were characterized by 18S and ITS ribosomal RNA gene sequencing. The pig gut mycobiota was dominated by yeasts, with a high prevalence and abundance of Kazachstania spp. Regarding protists, representatives of four genera (Blastocystis, Neobalantidium, Tetratrichomonas and Trichomitus) were predominant in more than the 80% of the pigs. Heritabilities for the diversity and abundance of gut eukaryotic communities were estimated with the subset of 60d aged piglets (N = 390). The heritabilities of α-diversity and of the abundance of fungal and protists genera were low, ranging from 0.15 to 0.28. A genome wide association study reported genetic variants related to the fungal α-diversity and to the abundance of Blastocystis spp. Annotated candidate genes were mainly associated with immunity, gut homeostasis and metabolic processes. Additionally, we explored the association of gut commensal eukaryotes with piglet body weight. Our results pointed to a positive contribution of fungi from the Kazachstania genus, while protists displayed both positive (Blastocystis and Entamoeba) and negative (Trichomitus) associations with piglet body weight. Conclusions Our results point towards a minor and taxa specific genetic control over the diversity and composition of the pig gut eukaryotic communities. Moreover, we provide evidences of the associations between piglets’ body weight after weaning and members from the gut fungal and protist eukaryote community. Overall, this study highlights the relevance of considering, along with that of bacteria, the contribution of the gut eukaryote communities to better understand host-microbiome association and their role on pig performance, welfare and health.

Feed efficiency (FE) is one of the most economically and environmentally relevant traits in the animal production sector. The objective of this study was to gain knowledge about the genetic control of FE in rabbits. To this end, GWASs were conducted for individual growth under two feeding regimes (full feeding and restricted) and FE traits collected from cage groups, using 114 604 autosome SNPs segregating in 438 rabbits. Two different models were implemented: (1) an animal model with a linear regression on each SNP allele for growth trait; and (2) a two‐trait animal model, jointly fitting the performance trait and each SNP allele content, for FE traits. This last modeling strategy is a new tool applied to GWAS and allows information to be considered from non‐genotyped individuals whose contribution is relevant in the group average traits. A total of 189 SNPs in 17 chromosomal regions were declared to be significantly associated with any of the five analyzed traits at a chromosome‐wide level. In 12 of these regions, 20 candidate genes were proposed to explain the variation of the analyzed traits, including genes such as FTO, NDUFAF6 and CEBPA previously associated with growth and FE traits in monogastric species. Candidate genes associated with behavioral patterns were also identified. Overall, our results can be considered as the foundation for future functional research to unravel the actual causal mutations regulating growth and FE in rabbits., This study was financed by Feed‐a‐Gene, which has received funding from the European Union’s H2020 program under grant agreement no. 63353 and the Spanish projects RTA2014‐00015‐C2‐01, RTA2011‐00064‐00‐00 and RTI2018‐097610‐R‐I00. M. Ballester is financially supported by a Ramon y Cajal contract (RYC‐2013‐12573) from the Spanish Ministry of Economy and Competitiveness.

[Background]: Intramuscular fat (IMF) content and composition have a strong impact on the nutritional and organoleptic properties of porcine meat. The goal of the current work was to compare the patterns of gene expression and the genetic determinism of IMF traits in the porcine gluteus medius (GM) and longissimus dorsi (LD) muscles., [Results]: A comparative analysis of the mRNA expression profiles of the pig GM and LD muscles in 16 Duroc pigs with available microarray mRNA expression measurements revealed the existence of 106 differentially expressed probes (fold-change > 1.5 and q-value, [Conclusions]: We have detected significant differences in the mRNA expression patterns of the porcine LD and GM muscles, evidencing that the transcriptomic profile of the skeletal muscle tissue is affected by anatomical, metabolic and functional factors. A highly significant association with IMF composition on SSC14 was replicated in both muscles, highlighting the existence of a common genetic determinism, but we also observed the existence of a few associations whose magnitude and significance varied between LD and GM muscles., Part of the research presented in this publication was funded by grants AGL2013–48742-C2–1-R and AGL2013–48742-C2–2-R awarded by the Spanish Ministry of Economy and Competitiveness and grant 2014 SGR 1528 from the Agency for Management of University and Research Grants of the Generalitat de Catalunya. We also acknowledge the support of the Spanish Ministry of Economy and Competitiveness for the Center of Excellence Severo Ochoa 2016–2019 (SEV-2015-0533) grant awarded to the Center for Research in Agricultural Genomics. Tainã F Cardoso was funded with a fellowship from the CAPES Foundation-Coordination of Improvement of Higher Education, Ministry of Education (MEC) of the Federal Government of Brazil. Rayner González-Prendes and Emilio Mármol-Sánchez were funded with FPU Ph.D. grants from the Spanish Ministry of Education (FPU12/00860 and FPU15/01733). Maria Ballester is financially supported by a Ramon y Cajal contract (RYC-2013-12573) from the Spanish Ministry of Economy and Competitiveness. Thanks also to the CERCA Programme of the Generalitat de Catalunya for their support.

The gut microbiota has been evolving with its host along the time creating a symbiotic relationship. In this study, we assess the role of the host genome in the modulation of the microbiota composition in pigs. Gut microbiota compositions were estimated through sequencing the V3-V4 region of the 16S rRNA gene from rectal contents of 285 pigs. A total of 1,261 operational taxonomic units were obtained and grouped in 18 phyla and 101 genera. Firmicutes (45.36%) and Bacteroidetes (37.47%) were the two major phyla obtained, whereas at genus level Prevotella (7.03%) and Treponema (6.29%) were the most abundant. Pigs were also genotyped with a high-throughput method for 45,508 single nucleotide polymorphisms that covered the entire pig genome. Subsequently, genome-wide association studies were made among the genotypes of these pigs and their gut microbiota composition. A total of 52 single-nucleotide polymorphisms distributed in 17 regions along the pig genome were associated with the relative abundance of six genera; Akkermansia, CF231, Phascolarctobacterium, Prevotella, SMB53, and Streptococcus. Our results suggest 39 candidate genes that may be modulating the microbiota composition and manifest the association between host genome and gut microbiota in pigs., This work was funded by the Fondo Europeo de Desarrollo Regional (FEDER) and the Spanish Ministerio de Economía y Competitividad (MINECO AGL2014-56369-C2 and AGL2017-82641-R). D. Crespo-Piazuelo was funded by a “Formació i Contractació de Personal Investigador Novell” (FI-DGR) Ph.D grant from the Generalitat de Catalunya (ECO/1788/2014) and by the PiGutNet COST Action (www.pigutnet.eu) for a Short Term Scientific Mission at the GABI laboratory (INRA, France) under the supervision of J. Estellé. Contract of L. Migura-Garcia was supported by INIA and the European Social Fund. L. Criado-Mesas was funded with a FPI grant from the AGL2014-56369-C2 project. M. Revilla was also funded by a FI-DGR (ECO/1639/2013). M. Ballester was financially supported by a “Ramón y Cajal” contract (RYC-2013-12573) from the Spanish Ministry of Economy and Competitiveness. We acknowledge the support of the Spanish Ministry of Economy and Competitiveness for the “Severo Ochoa Programme for Centres of Excellence in R&D” 2016-2019 (SEV-2015-0533) grant awarded to the Centre for Research in Agricultural Genomics and the CERCA Programme/Generalitat de Catalunya.

The selection in commercial swine breeds for meat-production efficiency has been increasing among the past decades, reducing the intramuscular fat content, which has changed the sensorial and technological properties of pork. Through processes of natural adaptation and selective breeding, the accumulation of mutations has driven the genetic divergence between pig breeds. The most common and well-studied mutations are single-nucleotide polymorphisms (SNPs). However, insertions and deletions (indels) usually represents a fifth part of the detected mutations and should also be considered for animal breeding. In the present study, three different programs (Dindel, SAMtools mpileup, and GATK) were used to detect indels from Whole Genome Sequencing data of Iberian boars and Landrace sows. A total of 1,928,746 indels were found in common with the three programs. The VEP tool predicted that 1,289 indels may have a high impact on protein sequence and function. Ten indels inside genes related with lipid metabolism were genotyped in pigs from three different backcrosses with Iberian origin, obtaining different allelic frequencies on each backcross. Genome-Wide Association Studies performed in the Longissimus dorsi muscle found an association between an indel located in the C1q and TNF related 12 (C1QTNF12) gene and the amount of eicosadienoic acid (C20:2(n-6))., This work was funded by the Ministerio de Economía y Competitividad (MINECO) and the Fondo Europeo de Desarrollo Regional (FEDER) projects AGL2014-56369-C2-2-R and AGL2017-82641-R. DCP was funded by a “Formacio i Contractacio de Personal Investigador Novell”(FIDGR) Ph.D grant from the Generalitat de Catalunya (ECO/1788/2014). LCM was funded with a FPI grant from the AGL2014-56369-C2 project. MR was also funded by a FI-DGR (ECO/1639/2013). MB was financially supported by a “Ramo´n y Cajal” contract (RYC-2013-12573) from the Spanish Ministry of Economy and Competitiveness. We acknowledge the support of the Spanish Ministry of Economy and Competitiveness for the “Severo Ochoa Programme for Centres of Excellence in R&D” 2016-2019 (SEV-2015-0533) grant awarded to the Centre for Research in Agricultural Genomics and the CERCA Programme / Generalitat de Catalunya.

The aim of this work was to study the genetic basis of the backfat expression of lipid-related genes associated with meat quality traits in pigs. We performed a genome-wide association study with the backfat gene expression measured in 44 genes by qPCR and the PorcineSNP60 BeadChip genotypes in 115 Iberian x Landrace backcross animals. A total of 193 expression-associated SNPs located in 19 chromosomal regions were associated with expression levels of ACSM5, ELOVL6, FABP4, FADS2, and SLC27A4 genes. Three expression quantitative trail loci (eQTLs) corresponding to ACSM5, FABP4, and FADS2 were classified as cis-acting eQTLs, whereas the remaining 16 eQTLs have trans-regulatory effects. Remarkably, a SNP in the ACSM5 promoter region and a SNP in the 3′UTR region of FABP4 were the most associated polymorphisms with the ACSM5 and FABP4 expression levels, respectively. Moreover, relevant lipid-related genes mapped in the trans-eQTLs regions associated with the ACSM5, FABP4, FADS2, and SLC27A4 genes. Interestingly, a trans-eQTL hotspot on SSC13 regulating the gene expression of ELOVL6, ELOLV5, and SCD, three important genes implicated in the elongation and desaturation of fatty acids, was identified. These findings provide new data to further understand the functional regulatory mechanisms implicated in the variation of fatty acid composition in pigs., This work has been funded by the MINECO AGL2014-56369-C2-R and AGL2017-82641-R grants. MR was funded by a Formació i Contractació de Personal Investigador Novell (FI-DGR) PhD grant from Generalitat de Catalunya (ECO/1639/2013). MB was financially supported by a Ramon y Cajal contract (RYC-2013-12573) from the Spanish Ministry of Economy and Competitiveness, APO was funded by a Personal Investigador en Formación (PIF) PhD grant from the Universitat Autònoma de Barcelona (458-01-1/2011), DCP by a Formació i Contractació de Personal Investigador Novell (FI-DGR) PhD grant from Generalitat de Catalunya (ECO/1788/2014) and LCM was funded by a FPI PhD grant from the Spanish Ministry of Economy and Competitiveness (BES-2015-075403). Funded by grant INIA CPE03-010-C3. We acknowledge financial support from the “Severo Ochoa Programme for Centres of Excellence in R&D” 2016-2019 (SEV‐2015‐0533) and the CERCA Programme / Generalitat de Catalunya.

This study aims identifying candidate genes and pathways associated with feed efficiency (FE) in pigs. Liver and duodenum transcriptomes of 37 gilts showing high and low residual feed intake (RFI) were analysed by RNA-Seq. Gene expression data was explored through differential expression (DE) and weighted gene co-expression network analyses. DE analysis revealed 55 and 112 differentially regulated genes in liver and duodenum tissues, respectively. Clustering genes according to their connectivity resulted in 23 (liver) and 25 (duodenum) modules of genes with a co-expression pattern. Four modules, one in liver (with 444 co-expressed genes) and three in duodenum (gathering 37, 126 and 41 co-expressed genes), were significantly associated with FE indicators. Intra-module analyses revealed tissue-specific candidate genes; 12 of these genes were also identified as DE between individuals with high and low RFI. Pathways enriched by the list of genes showing DE and/or belonging to FE co-expressed modules included response to oxidative stress, inflammation, immune response, lipid metabolism and thermoregulation. Low overlapping between genes identified in duodenum and liver tissues was observed but heat shock proteins were associated to FE in both tissues. Our results suggest tissue-specific rather than common transcriptome regulatory processes associated with FE in pigs., This work was funded by the European Union Seventh Framework Programme (FP7/2007–2013) as part of the ECO-FCE project under grant agreement No 311794. Y. Ramayo-Caldas was funded by the European Union H2020 Research and Innovation programme under Marie Skłodowska-Curie grant (P-Sphere) agreement No 6655919. M. Ballester is recipient of a Ramón y Cajal post-doctoral fellowship (RYC-2013–12573) from the Spanish Ministry of Economy, Industry and Competitiveness (MINECO).

In the present study, liver co-expression networks and expression Genome Wide Association Study (eGWAS) were performed to identify DNA variants and molecular pathways implicated in the functional regulatory mechanisms of meat quality traits in pigs. With this purpose, the liver mRNA expression of 44 candidates genes related with lipid metabolism was analysed in 111 Iberian x Landrace backcross animals. The eGWAS identified 92 eSNPs located in seven chromosomal regions and associated with eight genes: CROT, CYP2U1, DGAT1, EGF, FABP1, FABP5, PLA2G12A, and PPARA. Remarkably, cis-eSNPs associated with FABP1 gene expression which may be determining the C18:2(n-6)/C18:3(n-3) ratio in backfat through the multiple interaction of DNA variants and genes were identified. Furthermore, a hotspot on SSC8 associated with the gene expression of eight genes was identified and the TBCK gene was pointed out as candidate gene regulating it. Our results also suggested that the PI3K-Akt-mTOR pathway plays an important role in the control of the analysed genes highlighting nuclear receptors as the NR3C1 or PPARA. Finally, sex-dimorphism associated with hepatic lipid metabolism was identified with over-representation of female-biased genes. These results increase our knowledge of the genetic architecture underlying fat composition traits., This work was funded by MINECO AGL2011-29821-C02 and AGL2014-56369-C2 grants. M. Ballester is financially supported by a Ramon y Cajal contract (RYC-2013-12573) from the Spanish Ministry of Economy and Competitiveness. J. Corominas was funded by a FPI PhD grant from the Spanish Ministry of Education (BES-2009-081223). M. Revilla was funded by a Formació i Contractació de Personal Investigador Novell (FI-DGR) PhD grant from Generalitat de Catalunya (ECO/1639/2013).

Copy number variations (CNVs) are important genetic variants complementary to SNPs, and can be considered as biomarkers for some economically important traits in domestic animals. In the present study, a genomic analysis of porcine CNVs based on next-generation sequencing data was carried out to identify CNVs segregating in an Iberian x Landrace backcross population and study their association with fatty acid composition and growth-related traits. A total of 1,279 CNVs, including duplications and deletions, were detected, ranging from 106 to 235 CNVs across samples, with an average of 183 CNVs per sample. Moreover, we detected 540 CNV regions (CNVRs) containing 245 genes. Functional annotation suggested that these genes possess a great variety of molecular functions and may play a role in production traits in commercial breeds. Some of the identified CNVRs contained relevant functional genes (e.g., CLCA4, CYP4X1, GPAT2, MOGAT2, PLA2G2A and PRKG1, among others). The variation in copy number of four of them (CLCA4, GPAT2, MOGAT2 and PRKG1) was validated in 150 BC1_LD (25% Iberian and 75% Landrace) animals by qPCR. Additionally, their contribution regarding backfat and intramuscular fatty acid composition and growth–related traits was analyzed. Statistically significant associations were obtained for CNVR112 (GPAT2) for the C18:2(n-6)/C18:3(n-3) ratio in backfat and carcass length, among others. Notably, GPATs are enzymes that catalyze the first step in the biosynthesis of both triglycerides and glycerophospholipids, suggesting that this CNVR may contribute to genetic variation in fatty acid composition and growth traits. These findings provide useful genomic information to facilitate the further identification of trait-related CNVRs affecting economically important traits in pigs., This work has been funded by the MICINN AGL2014-56369-C2 and MINECO AGL2011-29821-C02. MR was funded by a Formacio i Contractacio de Personal Investigador Catalunya (ECO/1639/2013). APO was funded by a Personal Investigador en Formacion (PIF) PhD grant from the Universitat Autònoma de Barcelona (458-01-1/2011), DCP by a Formacio i Contractacio de Personal Investigador Novell (FIDGR) PhD grant from Generalitat de Catalunya. (ECO/1788/2014) and EP was funded by National Council of Scientific and Technological Development of Brazil (CNPq 202243/2014-1). MB was financially supported by a Ramon y Cajal contract (RYC-2013-12573) from the Spanish Ministry of Economy and Competitiveness. Novell (FI-DGR) PhD grant from Generalitat de

The FABP4 and FABP5 genes, coding for fatty acid transport proteins, have long been studied as positional candidate genes for SSC4 QTL affecting fat deposition and composition traits in pigs. Polymorphisms in these genes, FABP4:g.2634_2635insC and FABP5:g.3000T>G, have previously been associated with fatness traits in an Iberian by Landrace cross (IBMAP). The aim of the present work was to evaluate the functional implication of these genetic variants. For this purpose, FABP4 and FABP5 mRNA expression levels in 114 BC1_LD animals (25% Iberian × 75% Landrace) were analyzed using real-time quantitative PCR in backfat and muscle. FABP4 gene expression in backfat, but not in muscle, was associated with FABP4:g.2634_2635insC. In contrast, FABP5:g.3000T>G was not associated with gene expression levels. An expression-based genome-wide association study highlighted the FABP4:g.2634_2635insC polymorphism as the polymorphism most associated with FABP4 gene expression in backfat. Furthermore, other genomic regions associated in trans with the mRNA expression of FABP4 in backfat and FABP5 in muscle were also identified. Finally, two putative transcription binding sites for PPARG and NR4A2 may be affected by the FABP4:g.2634_2635insC polymorphism, modifying FABP4 gene expression. Our results reinforce FABP4 as a candidate gene for fatness traits on SSC4., This work was funded by MINECO AGL2011-29821-C02 and MICINN AGL2014-56369-C2 grants. M. Ballester is financially supported by a Ramon y Cajal Contract (RYC-2013-12573) from the Spanish Ministry of Economy and Competitiveness. We acknowledge financial support from CERCA Programme/Generalitat de Catalunya and Severo Ochoa Programme for Centres of Excellence in R&D 2016–2019 (SEV-2015-0533). A. Puig-Oliveras was funded by a PIF Ph.D. grant from Universitat Autònoma de Barcelona (458-01-1/2011), M. Revilla was funded by a Formació i Contractació de Personal Investigador Novell (FI-DGR) Ph.D. Grant from Generalitat de Catalunya (ECO/1639/2013).

Ramírez, Óscar et al.-- Data archiving: Ancient, Iberian, Spanish wild boar and Guatemalan reads have been submitted to SRA (accession SRP044261), aligned mitochondrial fasta file, plink files with genotypic data have been deposited in Dryad (doi:10.5061/dryad.sd784)., Ancient DNA (aDNA) provides direct evidence of historical events that have modeled the genome of modern individuals. In livestock, resolving the differences between the effects of initial domestication and of subsequent modern breeding is not straight forward without aDNA data. Here, we have obtained shotgun genome sequence data from a sixteenth century pig from Northeastern Spain (Montsoriu castle), the ancient pig was obtained from an extremely well-preserved and diverse assemblage. In addition, we provide the sequence of three new modern genomes from an Iberian pig, Spanish wild boar and a Guatemalan Creole pig. Comparison with both mitochondrial and autosomal genome data shows that the ancient pig is closely related to extant Iberian pigs and to European wild boar. Although the ancient sample was clearly domestic, admixture with wild boar also occurred, according to the D-statistics. The close relationship between Iberian, European wild boar and the ancient pig confirms that Asian introgression in modern Iberian pigs has not existed or has been negligible. In contrast, the Guatemalan Creole pig clusters apart from the Iberian pig genome, likely due to introgression from international breeds., This study was funded by 2010ACOM00030 (AGAUR) to MS, FEDER and BFU2012-34157 grant (Spain) to CLF, Consolider CSD2007-00036 ‘Centre for Research in Agrigenomics’ and AGL2010-14822 grants (Spain) to MPE. OR is a postdoctoral Researcher from the JAEDOC program cofounded by ESF, EB is recipient of a FPI grant from ministry of Research (Spain), IO of a predoctoral fellowship from the Basque Government (DEUI, Spain), WBP is funded by COLCIENCIAS (Francisco Jose´ de Caldas fellowship 497/2009, Colombia).

Background Studies to identify markers associated with beef tenderness have focused on Warner–Bratzler shear force (WBSF) but the interplay between the genes associated with WBSF has not been explored. We used the association weight matrix (AWM), a systems biology approach, to identify a set of interacting genes that are co-associated with tenderness and other meat quality traits, and shared across the Charolaise, Limousine and Blonde d’Aquitaine beef cattle breeds. Results Genome-wide association studies were performed using ~500K single nucleotide polymorphisms (SNPs) and 17 phenotypes measured on more than 1000 animals for each breed. First, this multi-trait approach was applied separately for each breed across 17 phenotypes and second, between- and across-breed comparisons at the AWM and functional levels were performed. Genetic heterogeneity was observed, and most of the variants that were associated with WBSF segregated within rather than across breeds. We identified 206 common candidate genes associated with WBSF across the three breeds. SNPs in these common genes explained between 28 and 30 % of the phenotypic variance for WBSF. A reduced number of common SNPs mapping to the 206 common genes were identified, suggesting that different mutations may target the same genes in a breed-specific manner. Therefore, it is likely that, depending on allele frequencies and linkage disequilibrium patterns, a SNP that is identified for one breed may not be informative for another unrelated breed. Well-known candidate genes affecting beef tenderness were identified. In addition, some of the 206 common genes are located within previously reported quantitative trait loci for WBSF in several cattle breeds. Moreover, the multi-breed co-association analysis detected new candidate genes, regulators and metabolic pathways that are likely involved in the determination of meat tenderness and other meat quality traits in beef cattle. Conclusions Our results suggest that systems biology approaches that explore associations of correlated traits increase statistical power to identify candidate genes beyond the one-dimensional approach. Further studies on the 206 common genes, their pathways, regulators and interactions will expand our knowledge on the molecular basis of meat tenderness and could lead to the discovery of functional mutations useful for genomic selection in a multi-breed beef cattle context. Electronic supplementary material The online version of this article (doi:10.1186/s12711-016-0216-y) contains supplementary material, which is available to authorized users.

The aim of this work is to better understand the genetic mechanisms determining two complex traits affecting porcine meat quality: intramuscular fat (IMF) content and its fatty acid (FA) composition. With this purpose, expression Genome-Wide Association Study (eGWAS) of 45 lipid-related genes associated with meat quality traits in swine muscle (Longissimus dorsi) of 114 Iberian × Landrace backcross animals was performed. The eGWAS identified 241 SNPs associated with 11 genes: ACSM5, CROT, FABP3, FOS, HIF1AN, IGF2, MGLL, NCOA1, PIK3R1, PLA2G12A and PPARA. Three expression Quantitative Trait Loci (eQTLs) for IGF2, ACSM5 and MGLL were identified, showing cis-acting effects, whereas 16 eQTLs had trans regulatory effects. A polymorphism in the ACSM5 promoter region associated with its expression was identified. In addition, strong candidate genes regulating ACSM5, FOS, PPARA, PIK3R1, PLA2G12A and HIF1AN gene expression were also seen. Notably, the analysis highlighted the NR3C1 transcription factor as a strong candidate gene involved in the regulation of the 45 genes analysed. Finally, the IGF2, MGLL, MC2R, ARHGAP6 and NR3C1 genes were identified as potential regulators co-localizing within QTLs for fatness and growth traits in the IBMAP population. The results obtained increase our knowledge in the functional regulatory mechanisms involved in these complex traits., This work was funded by the Ministerio de Economía y Competitividad (MINECO AGL2011-29821-C02 and AGL2014-56369-C2). A. Puig-Oliveras was funded by a “Personal Investigador en Formació” (PIF) grant from Universitat Autònoma de Barcelona (458-01-1/2011). M. Revilla was funded by a Formació i Contractació de Personal Investigador Novell (FI-DGR) Ph.D grant from the Generalitat de Catalunya (ECO/1639/2013). M. Ballester is financially supported by a Ramon y Cajal contract (RYC-2013-12573) from the Spanish Ministry of Economy and Competitiveness.

APOA2 is a protein implicated in triglyceride, fatty acid and glucose metabolism. In pigs, the APOA2 gene is located on pig chromosome 4 (SSC4) in a QTL region affecting fatty acid composition, fatness and growth traits. In this study, we evaluated APOA2 as a candidate gene for meat quality traits in an Iberian × Landrace backcross population. The APOA2:c.131T>A polymorphism, located in exon 3 of APOA2 and determining a missense mutation, was associated with the percentage of hexadecenoic acid [C16:1(n–9)], linoleic acid [C18:2(n–6)], α-linolenic acid [C18:3(n–3)], dihomo-gamma-linolenic acid [C20:3(n–6)] and polyunsaturated fatty acids (PUFAs) in backfat. Furthermore, this SNP was associated with the global mRNA expression levels of APOA2 in liver and was used as a marker to determine allelic expression imbalance by pyrosequencing. We determined an overexpression of the T allele in heterozygous samples with a mean ratio of 2.8 (T/A), observing a high variability in the allelic expression among individuals. This result suggests that complex regulatory mechanisms, beyond a single polymorphism (e.g. epigenetic effects or multiple cis-acting polymorphisms), may be regulating APOA2 gene expression., This work was funded by MINECO AGL2011-29821-C02 and AGL2014-56369-C2 grants. M. Ballester is financially supported by a Ramon y Cajal contract (RYC-2013-12573) from the Spanish Ministry of Economy and Competitiveness. A. Puig-Oliveras was funded by a PIF PhD grant from Universitat Autònoma de Barcelona (458-01-1/2011), J. Corominas was funded by a FPI PhD grant from the Spanish Ministry of Education (BES-2009-081223), M. Revilla was funded by a Formació i Contractació de Personal Investigador Novell (FI-DGR) PhD grant from Generalitat de Catalunya (ECO/1639/2013) and J.A.P. Marchesi was funded by a CNPq scholarship from the Brazilian Science Without Borders Program.

[Background]: In previous studies on an Iberian x Landrace cross, we have provided evidence that supported the porcine ELOVL6 gene as the major causative gene of the QTL on pig chromosome 8 for palmitic and palmitoleic acid contents in muscle and backfat. The single nucleotide polymorphism (SNP) ELOVL6:c.-533C > T located in the promoter region of ELOVL6 was found to be highly associated with ELOVL6 expression and, accordingly, with the percentages of palmitic and palmitoleic acids in longissimus dorsi and adipose tissue. The main goal of the current work was to further study the role of ELOVL6 on these traits by analyzing the regulation of the expression of ELOVL6 and the implication of ELOVL6 polymorphisms on meat quality traits in pigs., [Results]: High-throughput sequencing of BAC clones that contain the porcine ELOVL6 gene coupled to RNAseq data re-analysis showed that two isoforms of this gene are expressed in liver and adipose tissue and that they differ in number of exons and 3’UTR length. Although several SNPs in the 3’UTR of ELOVL6 were associated with palmitic and palmitoleic acid contents, this association was lower than that previously observed with SNP ELOVL6:c.-533C > T. This SNP is in full linkage disequilibrium with SNP ELOVL6:c.-394G > A that was identified in the binding site for estrogen receptor alpha (ERα). Interestingly, the ELOVL6:c.-394G allele is associated with an increase in methylation levels of the ELOVL6 promoter and with a decrease of ELOVL6 expression. Therefore, ERα is clearly a good candidate to explain the regulation of ELOVL6 expression through dynamic epigenetic changes in the binding site of known regulators of ELOVL6 gene, such as SREBF1 and SP1., [Conclusions]: Our results strongly suggest the ELOVL6:c.-394G > A polymorphism as the causal mutation for the QTL on pig chromosome 8 that affects fatty acid composition in pigs., This work was funded by the Ministerio de Economía y Competitividad project AGL2011-29821-C02. J Corominas was funded by a FPI PhD grant from the Spanish Ministerio de Educación (BES-2009-081223), M Revilla was funded by a Formació i Contractació de Personal Investigador Novell (FI-DGR) PhD grant from Generalitat de Catalunya (ECO/1639/2013), A Puig-Oliveras was funded by a PIF PhD grant from the Universitat Autónoma de Barcelona (458-01-1/2011) and JAP Marchesi was funded by a CNPq scholarship from the Brazilian Science Without Borders Program.

β-lactoglobulin polymorphisms have been reported in the milk of different goat breeds, although no genetic variants affecting the protein have been characterized. In the present study, we amplified and sequenced the proximal promoter and the first six exons containing the entire coding region for the β-lactoglobulin gene in eleven goat breeds from Spain, France, Italy, Switzerland, Senegal and Asia to identify genetic variants. Fifteen polymorphisms were detected, nine in the promoter region and six in the exons of the β-lactoglobulin gene. All polymorphisms were single nucleotide substitutions with the exception of one deletion/insertion in the promoter region. The polymorphisms in the coding region did not produce any amino acid change. In addition, pyrosequencing technology was used to genotype four polymorphisms in the promoter region in 200 goats belonging to eleven breeds. Differences in allelic frequencies for these polymorphisms between breeds are described and a specific polymorphism for the Italian populations was identified. Finally, the analysis of association between these four promoter point mutations was investigated resulting in five haplotypes, GCGC being the most frequent haplotype in all breeds analysed.

Background Fatty acids (FA) play a critical role in energy homeostasis and metabolic diseases; in the context of livestock species, their profile also impacts on meat quality for healthy human consumption. Molecular pathways controlling lipid metabolism are highly interconnected and are not fully understood. Elucidating these molecular processes will aid technological development towards improvement of pork meat quality and increased knowledge of FA metabolism, underpinning metabolic diseases in humans. Results The results from genome-wide association studies (GWAS) across 15 phenotypes were subjected to an Association Weight Matrix (AWM) approach to predict a network of 1,096 genes related to intramuscular FA composition in pigs. To identify the key regulators of FA metabolism, we focused on the minimal set of transcription factors (TF) that the explored the majority of the network topology. Pathway and network analyses pointed towards a trio of TF as key regulators of FA metabolism: NCOA2, FHL2 and EP300. Promoter sequence analyses confirmed that these TF have binding sites for some well-know regulators of lipid and carbohydrate metabolism. For the first time in a non-model species, some of the co-associations observed at the genetic level were validated through co-expression at the transcriptomic level based on real-time PCR of 40 genes in adipose tissue, and a further 55 genes in liver. In particular, liver expression of NCOA2 and EP300 differed between pig breeds (Iberian and Landrace) extreme in terms of fat deposition. Highly clustered co-expression networks in both liver and adipose tissues were observed. EP300 and NCOA2 showed centrality parameters above average in the both networks. Over all genes, co-expression analyses confirmed 28.9% of the AWM predicted gene-gene interactions in liver and 33.0% in adipose tissue. The magnitude of this validation varied across genes, with up to 60.8% of the connections of NCOA2 in adipose tissue being validated via co-expression. Conclusions Our results recapitulate the known transcriptional regulation of FA metabolism, predict gene interactions that can be experimentally validated, and suggest that genetic variants mapped to EP300, FHL2, and NCOA2 modulate lipid metabolism and control energy homeostasis in pigs., This work was funded by MICINN project AGL2011-29821-C02 (Ministerio de Economía y Competitividad), and by the Innovation Consolider-Ingenio 2010 Program (CSD2007-00036, Centre for Research in Agrigenomics). We thank Dr. Nick Hudson for insightful suggestions and manuscript proofread. Y. Ramayo-Caldas was funded by a FPU PhD grant from the Spanish Ministerio de Educación (AP2008-01450).

[Background]: Fat content and fatty acid composition in swine are becoming increasingly studied because of their effect on sensory and nutritional quality of meat. A QTL (quantitative trait locus) for fatty acid composition in backfat was previously detected on porcine chromosome 8 (SSC8) in an Iberian x Landrace F2 intercross. More recently, a genome-wide association study detected the same genomic region for muscle fatty acid composition in an Iberian x Landrace backcross population. ELOVL6, a strong positional candidate gene for this QTL, contains a polymorphism in its promoter region (ELOVL6:c.-533C, [Results]: Two trait-associated SNP regions were identified at 93 Mb and 119 Mb on SSC8. The strongest statistical signals of both regions were observed for palmitoleic acid (C16:1(n-7)) content and C18:0/C16:0 and C18:1(n-7)/C16:1(n-7) elongation ratios. MAML3 and SETD7 are positional candidate genes in the 93 Mb region and two novel microsatellites in MAML3 and nine SNPs in SETD7 were identified. No significant association for the MAML3 microsatellite genotypes was detected. The SETD7:c.700G > T SNP, although statistically significant, was not the strongest signal in this region. In addition, the expression of MAML3 and SETD7 in liver and adipose tissue varied among animals, but no association was detected with the polymorphisms in these genes. In the 119 Mb region, the ELOVL6:c.-533C > T polymorphism showed a strong association with percentage of palmitic and palmitoleic fatty acids and elongation ratios in backfat., [Conclusions]: Our results suggest that the polymorphisms studied in MAML3 and SETD7 are not the causal mutations for the QTL in the 93 Mb region. However, the results for ELOVL6 support the hypothesis that the ELOVL6:c.-533C > T polymorphism has a pleiotropic effect on backfat and intramuscular fatty acid composition and that it has a role in the determination of the QTL in the 119 Mb region., This work was funded by MICINN AGL2008-04818-C03/GAN and MINECO AGL2011-29821-C02 and the Innovation Programme Consolider-Ingenio 2010 (CSD2007-00036). M. Revilla is a Master’s student of Animal Breeding and Biotechnology of Reproduction (Polytechnical University of Valencia and Autonomous University of Barcelona). Y. Ramayo-Caldas was funded by a FPU grant (AP2008-01450), J. Corominas by a FPI scholarship from the Ministry of Education (BES-2009-018223) and A. Puig-Oliveras by a PIF scholarship (458-01-1/2011).

[Background]: Besides having an impact on human health, the porcine muscle fatty acid profile determines meat quality and taste. The RNA-Seq technologies allowed us to explore the pig muscle transcriptome with an unprecedented detail. The aim of this study was to identify differentially-expressed genes between two groups of 6 sows belonging to an Iberian × Landrace backcross with extreme phenotypes according to FA profile., [Results]: We sequenced the muscle transcriptome acquiring 787.5 M of 75 bp paired-end reads. About 85.1% of reads were mapped to the reference genome. Of the total reads, 79.1% were located in exons, 6.0% in introns and 14.9% in intergenic regions, indicating expressed regions not annotated in the reference genome. We identified a 34.5% of the intergenic regions as interspersed repetitive regions. We predicted a total of 2,372 putative proteins. Pathway analysis with 131 differentially-expressed genes revealed that the most statistically-significant metabolic pathways were related with lipid metabolism. Moreover, 18 of the differentially-expressed genes were located in genomic regions associated with IMF composition in an independent GWAS study in the same genetic background. Thus, our results indicate that the lipid metabolism of FAs is differently modulated when the FA composition in muscle differs. For instance, a high content of PUFA may reduce FA and glucose uptake resulting in an inhibition of the lipogenesis. These results are consistent with previous studies of our group analysing the liver and the adipose tissue transcriptomes providing a view of each of the main organs involved in lipid metabolism., [Conclusions]: The results obtained in the muscle transcriptome analysis increase the knowledge of the gene regulation of IMF deposition, FA profile and meat quality, in terms of taste and nutritional value. Besides, our results may be important in terms of human health., This work has been funded by MICINN project AGL2011-29821-C02 (Ministerio de Economia y Competitividad) and by the Innovation ConsoliderIngenio 2010 Program (CSD2007-00036). A. Puig-Oliveras was funded by a Personal Investigador en Formacion (PIF) PhD grant from the Universitat Autonoma de Barcelona (458-01-1/2011). J. Corominas was funded by a Formacion de Personal Investigador (FPI) PhD grant from Spanish Ministerio de Educacion (BES-2009-018223), Y. Ramayo by a Formacion del Profesorado Universitario (FPU) PhD grant (AP2008-01450).

[Background]: Several QTLs have been identified for major economically relevant traits in livestock, such as growth and meat quality, revealing the complex genetic architecture of these traits. The use of network approaches considering the interactions of multiple molecules and traits provides useful insights into the molecular underpinnings of complex traits. Here, a network based methodology, named Association Weight Matrix, was applied to study gene interactions and pathways affecting pig conformation, growth and fatness traits., [Results]: The co-association network analysis underpinned three transcription factors, PPARγ, ELF1, and PRDM16 involved in mesoderm tissue differentiation. Fifty-four genes in the network belonged to growth-related ontologies and 46 of them were common with a similar study for growth in cattle supporting our results. The functional analysis uncovered the lipid metabolism and the corticotrophin and gonadotrophin release hormone pathways among the most important pathways influencing these traits. Our results suggest that the genes and pathways here identified are important determining either the total body weight of the animal and the fat content. For instance, a switch in the mesoderm tissue differentiation may determinate the age-related preferred pathways being in the puberty stage those related with the miogenic and osteogenic lineages; on the contrary, in the maturity stage cells may be more prone to the adipocyte fate. Hence, our results demonstrate that an integrative genomic co-association analysis is a powerful approach for identifying new connections and interactions among genes., [Conclusions]: This work provides insights about pathways and key regulators which may be important determining the animal growth, conformation and body proportions and fatness traits. Molecular information concerning genes and pathways here described may be crucial for the improvement of genetic breeding programs applied to pork meat production., This work has been funded by MICINN project AGL2011-29821-C02 (Ministerio de Economia y Competitividad), and by the Innovation Consolider-Ingenio 2010 Program (CSD2007-00036, Centre de Recerca en Agrigenomica). APO was funded by a Personal Investigador en Formacion (PIF) PhD grant from the Universitat Autonoma de Barcelona (458-01-1/2011). JC was funded by a Formacion de Personal Investigador (FPI) PhD grant from Spanish Ministerio de Educacion (BES-2009-018223), YRC by a Formacion del Profesorado Universitario (FPU) PhD grant from the Spanish Ministerio (AP2008-01450), and MR by a Formacio i Contractacio de Personal Investigador Novell (FI-DGR) PhD grant from Generalitat de Catalunya (ECO/1639/2013)

Suppressive subtractive hybridization libraries from oviduct at 62 h post-mating of two lines of rabbits divergently selected for uterine capacity were generated to identify differentially expressed genes. A total of 438 singletons and 126 contigs were obtained by cluster assembly and sequence alignment of 704 expressed sequence tags (ESTs), of which 54% showed homology to known proteins of the non-redundant NCBI databases. Differential screening by dot blot validated 71 ESTs, of which 47 showed similarity to known genes. Transcripts of genes were functionally annotated in the molecular function and the biological process gene ontology categories using the blast2go software and were assigned to reproductive developmental process, immune response, amino acid metabolism and degradation, response to stress and apoptosis terms. Finally, three interesting genes, PGR, HSD17B4 and ERO1L, were identified as overexpressed in the low line using RT-qPCR. Our study provides a list of candidate genes that can be useful to understanding the molecular mechanisms underlying the phenotypic differences observed in early embryo survival and development traits., This study was funded with project AGL2005-07624-C03.

[Background]: The ELOVL fatty acid elongase 6 (ELOVL6), the only elongase related to de novo lipogenesis, catalyzes the rate-limiting step in the elongation cycle by controlling the fatty acid balance in mammals. It is located on pig chromosome 8 (SSC8) in a region where a QTL affecting palmitic, and palmitoleic acid composition was previously detected, using an Iberian x Landrace intercross. The main goal of this work was to fine-map the QTL and to evaluate the ELOVL6 gene as a positional candidate gene affecting the percentages of palmitic and palmitoleic fatty acids in pigs., [Methodology and Principal Findings]: The combination of a haplotype-based approach and single-marker analysis allowed us to identify the main, associated interval for the QTL, in which the ELOVL6 gene was identified and selected as a positional candidate gene. A polymorphism in the promoter region of ELOVL6, ELOVL6:c.-533C>T, was highly associated with the percentage of palmitic and palmitoleic acids in muscle and backfat. Significant differences in ELOVL6 gene expression were observed in backfat when animals were classified by the ELOVL6:c.-533C>T genotype. Accordingly, animals carrying the allele associated with a decrease in ELOVL6 gene expression presented an increase in C16:0 and C16:1(n-7) fatty acid content and a decrease of elongation activity ratios in muscle and backfat. Furthermore, a SNP genome-wide association study with ELOVL6 relative expression levels in backfat showed the strongest effect on the SSC8 region in which the ELOVL6 gene is located. Finally, different potential genomic regions associated with ELOVL6 gene expression were also identified by GWAS in liver and muscle, suggesting a differential tissue regulation of the ELOVL6 gene., [Conclusions and Significance]: Our results suggest ELOVL6 as a potential causal gene for the QTL analyzed and, subsequently, for controlling the overall balance of fatty acid composition in pigs., This study has been funded by MICINN projects (AGL2008-04818-C03 and AGL2011-29821-CO2) and the Innovation Consolider-Ingenio 2012 Program (CSD2007-00036, Center for Research in Agrigenomics). J. Corominas was funded by a Formación de Personal Investigador (FPI) PhD grant from Spanish Ministerio de Educación (BES-2009-018223), Y. Ramayo by a Formación del Profesorado Universitario (FPU) PhD grant (AP2008-01450) and A. Puig was funded by a Personal Investigador en Formación (PIF) PhD grant from the Universitat Autónoma de Barcelona (458-01-1/2011).

At present, a wide range of molecular sex-typing protocols in wild and domestic animals are available. In pigs, most of these methods are based on PCR amplification of X–Y homologous genes followed by gel electrophoresis which is time-consuming and in some cases expensive. In this paper, we describe, for the first time, a SYBR green-based quantitative real-time PCR (qPCR) assay using an X-linked gene, the glycoprotein M6B, for genetic sexing of pigs. Taking into account the differences in the glycoprotein M6B gene copy number between genders, we determine the correct sex of 54 pig samples from either diaphragm or hair follicle from different breeds using the 2−ΔΔCT method for relative quantification. Our qPCR assay represents a quick, inexpensive, and reliable tool for sex determination in pigs. This new protocol could be easily adapted to other species in which the sex determination was required., This study has been funded by the Innovation Consolider-Ingenio 2010 Program (CSD2007-00036, Centre for Research in Agrigenomics). Y. Ramayo-Caldas was funded by a FPU PhD grant (AP2008-01450).

Long-chain acyl-CoA synthetase (ACSL) family members catalyse the formation of long-chain acyl-CoA from fatty acid, ATP and CoA, thus playing an important role in both de novo lipid synthesis and fatty acid catabolism. Previous studies in our group evaluated ACSL4 as a positional candidate gene for quantitative trait loci located on chromosome X in an Iberian × Landrace cross. A DQ144454:c.2645G>A SNP located in the 3′ untranslated region of the ACSL4 gene was associated with the percentages of oleic and monounsaturated fatty acids. The aim of the present work was to evaluate the functional implication of this genetic variant. An expression analysis was performed for 120 individuals with different genotypes for the DQ144454:c.2645G>A polymorphism using real-time quantitative PCR. Differences between genotypes were identified in liver, with the ACSL4 mRNA expression levels higher in animals with the G allele than in animals with the A allele. A SNP genome-wide association study with ACSL4 relative expression levels showed significant positions on chromosomes 6 and 12. Description of positional candidate genes for ACSL4 regulation on chromosomes 6 and 12 is provided., This study has been funded by MICINN project (AGL2008-04818-C03). J. Corominas was funded by a FPI PhD grant from the Spanish Ministerio de Educación (BES-2009-018223) and Y. Ramayo by a FPU PhD grant (AP2008-01450).

International audience; During the development of tissues, complex programs take place to reach terminally differentiated states with specific gene expression profiles. Epigenetic regulations such as histone modifications and chromatin condensation have been implicated in the short and long-term control of transcription. It has recently been shown that the 3D spatial organization of chromosomes in the nucleus also plays a role in genome function. Indeed, the eukaryotic interphase nucleus contains sub-domains that are characterized by their enrichment in specific factors such as RNA Polymerase II, splicing machineries or heterochromatin proteins which render portions of the genome differentially permissive to gene expression. The positioning of individual genes relative to these sub-domains is thought to participate in the control of gene expression as an epigenetic mechanism acting in the nuclear space. Here, we review what is known about the sub-nuclear organization of mammary epithelial cells in connection with gene expression and epigenetics. Throughout differentiation, global changes in nuclear architecture occur, notably with respect to heterochromatin distribution. The positions of mammary-specific genes relative to nuclear sub-compartments varies in response to hormonal stimulation. The contribution of tissue architecture to cell differentiation in the mammary gland is also seen at the level of nuclear organization, which is sensitive to microenvironmental stimuli such as extracellular matrix signaling. In addition, alterations in nuclear organization are concomitant with immortalization and carcinogenesis. Thus, the fate of cells appears to be controlled by complex pathways connecting external signal integration, gene expression, epigenetic modifications and chromatin organization in the nucleus.

Signal transducers and activators of transcription (STATs) are a family of transcription factors. STAT5A and 5B are two highly related proteins encoded by two distinct genes. Transgenic knockout mice studies have indicated the importance of STAT5 proteins for the regulation of both lactation and growth performance. Moreover, different studies determine the role of STAT5 proteins in the modulation of adipocyte function. In this study, we sequenced one fragment of STAT5B gene from animals of six breeds (Duroc, Iberian, Landrace, Large White, Pie´train and Meishan) to identify genetic variants. A G/A single nucleotide polymorphism in intron 14 creates a polymorphic PstI restriction site and was genotyped by polymerase chain reaction restriction fragment length polymorphism in the six breeds. Allele G was only present in Large White, Pie´train and Meishan populations, detecting only G allele in this last pig breed. The STAT5B gene was located by radiation hybrid mapping to porcine chromosome 12, within the confidence interval for the fatty acid composition quantitative trait loci, previously identified in an Iberian · Landrace cross.

The mitochondrial hypervariable region 1 (HVR1) sequence of 67 goats belonging to the Girgentana, Maltese and Derivata di Siria breeds was partially sequenced in order to present the first phylogenetic characterization of Sicilian goat breeds. These sequences were compared with published sequences of Indian and Pakistani domestic goats and wild goats. Mitochondrial lineage A was observed in most of the Sicilian goats. However, three Girgentana haplotypes were highly divergent from the Capra hircus clade, indicating that a new mtDNA lineage in domestic goats was found.

A recessive insertional mutation was identified in one line of transgenic mice for the caprine beta-lactoglobulin (betaLG) gene. High mortality after birth, a significant reduction in postnatal growth and adult body size, changes in the morphometric features of the head, and infertility are the most prominent phenotypic traits of the mutant animals. Molecular cloning and sequencing of the transgene insertion site showed that 22 copies of the betaLG transgene are inserted in an intronic region of the phospholipase C-beta1 (PLC-beta1) gene, which plays a pivotal role in modulating different cellular functions. As a result of the insertional mutation (PLC-beta1(betaLG) mutation), a hybrid messenger RNA (mRNA) between the mouse PLC-beta1 and the goat betaLG genes is transcribed. The tissue-specific pattern of expression of this hybrid mRNA in PLC-beta1(betaLG) homozygotes is equivalent to that of the endogenous PLC-beta1 mRNA in nontransgenic animals, which is reported for the first time in this species, but expression levels are significantly reduced. Although the hybrid PLCbeta1-betaLG mRNA contains all the essential information to produce a PLCbeta1 protein that could be activated, this protein was not detected by Western blot. The PLC-beta1(betaLG) mouse model described here represents a useful tool to investigate the role of the PLC-beta1 gene in the molecular mechanisms underlying growth and fertility.

Background New 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. Results The 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. Conclusions In 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., This work was funded by MICINN projects AGL2008-04818-C03/GAN and AGL2011-29821-C02 (Ministerio de Economía y Competitividad), and by the Innovation Consolider-Ingenio 2010 Program (CSD2007-00036, Centre for Research in Agrigenomics). Y. Ramayo-Caldas was funded by a FPU PhD grant from the Spanish Ministerio de Educación (AP2008-01450), J. Corominas was funded by a FPI PhD grant from the Spanish Ministerio de Educación (BES-2009-018223), A. Esteve-Codina is recipient of a FPI PhD fellowship from the Ministerio de Educación (BES-2008-005772), Spain.