Association of genetic variants and expression levels of porcine FABP4 and FABP5 genes.

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.

Using RNA-Seq SNP data to reveal potential causal mutations related to pig production traits and RNA editing.

RNA-Seq technology is widely used in quantitative gene expression studies and identification of non-annotated transcripts. However this technology also can be used for polymorphism detection and RNA editing in transcribed regions in an efficient and cost-effective way. This study used SNP data from an RNA-Seq assay to identify genes and mutations underlying production trait variations in an experimental pig population. The hypothalamic and hepatic transcriptomes of nine extreme animals for growth and fatness from an (Iberian × Landrace) × Landrace backcross were analyzed by RNA-Seq methodology, and SNP calling was conducted. More than 125 000 single nucleotide variants (SNVs) were identified in each tissue, and 78% were considered to be potential SNPs, those SNVs segregating in the context of this study. Potential informative SNPs were detected by considering those showing a homozygous or heterozygous genotype in one extreme group and the alternative genotype in the other group. In this way, 4396 and 1862 informative SNPs were detected in hypothalamus and liver respectively. Out of the 32 SNPs selected for validation, 25 (80%) were confirmed as actual SNPs. Association analyses for growth, fatness and premium cut yields with 19 selected SNPs were carried out, and four potential causal genes (RETSAT, COPA, RNMT and PALMD) were identified. Interestingly, new RNA editing modifications were detected and validated for the NR3C1:g.102797 (ss1985401074) and ACSM2B:g.13374 (ss1985401075) positions and for the COG3:g3.4525 (ss1985401087) modification previously identified across vertebrates, which could lead to phenotypic variation and should be further investigated.

Analysis of the porcine APOA2 gene expression in liver, polymorphism identification and association with fatty acid composition traits.

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.

Gene co-expression network analysis for porcine intramuscular fatty acid composition.

In pigs, meat quality depends markedly on the fatty acid (FA) content and composition of the intramuscular fat, which is partly determined by the gene expression in this tissue. The aim of this work was to identify the link between muscle gene expression and its FA composition. In an (Iberian × Duroc) × Duroc backcrossed pig population, we identified modules of co-expressed genes, and correlation analyses were performed for each of them versus the phenotypes, finding four relevant modules. Two of the modules were positively correlated with saturated FAs (SFAs) and monounsaturated FAs (MUFAs), while negatively correlated with polyunsaturated FAs (PUFAs) and the omega-6/omega-3 ratio. The gene-enrichment analysis showed that these modules had over-representation of pathways related with the biosynthesis of unsaturated FAs, the Peroxisome proliferator-activated receptor signalling pathway and FA elongation. The two other relevant modules were positively correlated with PUFA and the n-6/n-3 ratio, but negatively correlated with SFA and MUFA. In this case, they had an over-representation of pathways related with fatty and amino acid degradation, and with oxidative phosphorylation. Using a graphical Gaussian model, we inferred a network of connections between the genes within each module. The first module had 52 genes with 87 connections, and the most connected genes were ADIPOQ, which is related with FA oxidation, and ELOVL6 and FABP4, both involved in FA metabolism. The second module showed 196 genes connected by 263 edges, being FN1 and MAP3K11 the most connected genes. On the other hand, the third module had 161 genes connected by 251 edges and ATG13 was the top neighbouring gene, while the fourth module had 224 genes and 655 connections, and its most connected genes were related with mitochondrial pathways. Overall, this work successfully identified relevant muscle gene networks and modules linked with FA composition, providing further insights on how the physiology of the pigs influences FA composition.

Identification of differentially expressed genes in the oviduct of two rabbit lines divergently selected for uterine capacity using suppression subtractive hybridization.

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.

Transcriptional analysis of intramuscular fatty acid composition in the longissimus thoracis muscle of Iberian × Landrace back-crossed pigs.

This study aimed at identifying differential gene expression conditional on the fatty acid profile of the longissimus thoracis (Lt) muscle, a prime cut of economic relevance for fresh and cured pork production. A population of 110 Iberian (25%) × Landrace (75%) back-crossed pigs was used, because these two breeds exhibit extreme profiles of intramuscular saturated fatty acid, monounsaturated fatty acid (MUFA) and polyunsaturated fatty acid (PUFA) contents. Total RNA from Lt muscle was individually hybridized to GeneChip Porcine Genome arrays (Affymetrix). A principal component analysis was performed with data from the 110 animals to select 40 extreme animals based on the total fatty acid profile and the MUFA composition (MAP). Comparison of global transcription levels between extreme fatty acid profile pigs (n = 40) resulted in 219 differentially expressed probes (false discovery rate <0.10). Gene ontology, pathway and network analysis indicated that animals with higher percentages of PUFA exhibit a shift toward a more oxidative muscular metabolism state, with a raise in mitochondria function (PPARGC1A, ATF2), fatty acid uptake and oxidation (FABP5, MGLL). On the other hand, 87 probes were differentially expressed between MUFA composition groups (n = 40; false discovery rate <0.10). In particular, muscles rich in n-7 MUFA expressed higher levels of genes involved in lipid metabolism (GLUL, CRAT, PLA2G15) and lower levels of fatty acid elongation genes (ELOVL5). Moreover, the chromosomal position of FABP5, PAQR3, MGLL, PPARGC1A, GLUL and ELOVL5 co-localized with very relevant QTL for fat deposition and composition described in the same resource population. This study represents a complementary approach to identifying genes underlying these QTL effects.

Influence of the FecX(R) allele in heterozygous ewes on follicular population and outcomes of IVP and ET using LOPU-derived oocytes.

Ewes heterozygous for the FecX(R) allele (R+) in the bone morphogenetic protein 15 (BMP15) gene display increased ovulation rate and prolificacy. Besides this phenotypic advantage, the influence of the FecX(R) allele on follicle number and size, oocyte competence and in vitro production (IVP) remains undefined. With these aims, 8 R+ and 8 wild-type (++) ewes were subjected to 2 laparoscopic ovum pick-up (LOPU) trials (four sessions per trial; two with and two without FSH) and subsequent IVP and fresh embryo transfer. All follicles >3 mm were punctured (n = 1673). Genotype did not significantly affect the number of punctured follicles per ewe and session (10.4 and 10.2 in R+ and ++ untreated ewes, 17.4 and 14.3 in R+ and ++ FSH-treated ewes, respectively), but follicular diameter of R+ ewes was significantly reduced compared with ++ ewes (-0.2 mm in untreated and -0.8 mm in FSH-treated ewes; p < 0.01). R+ ewes showed higher recovery rate and increased numbers of total and suitable cumulus-oocyte complexes for in vitro maturation (IVM). Similar rates of day 8 blastocysts were observed in R+ (36.1%, 147/407) and ++ (32.6%, 100/307) ewes, but the final output of day 8 blastocysts per ewe and session was higher in R+ ewes (+0.75; p < 0.005), without differences in survival rate at birth of the transferred embryos (40.4%, 21/52 vs 36.4%, 16/44, respectively). In conclusion, a higher number of oocytes proven to be competent for in vitro development and embryo survival after transfer are recovered from R+ ewes, despite the lower mean size of their follicles at puncture.

Identification of genomic regions, genetic variants and gene networks regulating candidate genes for lipid metabolism in pig muscle.

The intramuscular fat content and fatty acid composition of porcine meat have a significant impact on its quality and nutritional value. This research aimed to investigate the expression of 45 genes involved in lipid metabolism in the longissimus dorsi muscle of three experimental pig backcrosses, with a 25% of Iberian background. To achieve this objective, we conducted an expression Genome-Wide Association Study (eGWAS) using gene expression levels in muscle measured by high-throughput real-time qPCR for 45 target genes and genotypes from the PorcineSNP60 BeadChip or Axiom Porcine Genotyping Array and 65 single nucleotide polymorphisms (SNPs) located in 20 genes genotyped by a custom-designed Taqman OpenArray in a cohort of 354 animals. The eGWAS analysis identified 301 eSNPs associated with 18 candidate genes (ANK2, APOE, ARNT, CIITA, CPT1A, EGF, ELOVL6, ELOVL7, FADS3, FASN, GPAT3, NR1D2, NR1H2, PLIN1, PPAP2A, RORA, RXRA and UCP3). Three cis-eQTL (expression quantitative trait loci) were identified for GPAT3, RXRA, and UCP3 genes, which indicates that a genetic polymorphism proximal to the same gene is affecting its expression. Furthermore, 24 trans-eQTLs were detected, and eight candidate regulatory genes were located in these genomic regions. Additionally, two trans-regulatory hotspots in Sus scrofa chromosomes 13 and 15 were identified. Moreover, a co-expression analysis performed on 89 candidate genes and the fatty acid composition revealed the regulatory role of four genes (FABP5, PPARG, SCD, and SREBF1). These genes modulate the levels of α-linolenic, arachidonic, and oleic acids, as well as regulating the expression of other candidate genes associated with lipid metabolism. The findings of this study offer novel insights into the functional regulatory mechanism of genes involved in lipid metabolism, thereby enhancing our understanding of this complex biological process.

Lack of Jun-N-terminal kinase 3 (JNK3) does not protect against neurodegeneration induced by 3-nitropropionic acid.

AIMS: 3-Nitropropionic acid (3-NP) is a toxin that replicates most of the clinical and pathophysiological symptoms of Huntington's disease, inducing neurodegeneration in the striatum due to the inhibition of mitochondrial succinate dehydrogenase. Different pathways have been implicated in the cell death induced by 3-NP in rodents. One of them is the Jun-N-terminal kinase (JNK) pathway, which may play a role in the neurodegenerative process in different diseases. Moreover, the lack of one isoform of JNK (JNK3) has been associated with neuroprotection in different experimental models of neurodegeneration. Therefore, in the present study the role of JNK3 in the experimental Huntington's model induced by 3-NP administration was evaluated. METHODS: 3-NP was intraperitoneally administered once a day for 3 days to wild-type and Jnk3-null mice. Coronal brain sections were used to determine cell death and astrogliosis in striatum. Western blots were performed to determine the involvement of different pathways in both wild-type and Jnk3-null mice. RESULTS: Although JNK activation was observed following 3-NP administration, the results indicate that the lack of JNK3 does not confer neuroprotection against 3-NP toxicity. Thus, other pathways must be involved in the neurodegeneration induced in this model. One of the possible pathways towards 3-NP-induced apoptosis could involve the calpains, as their activity was increased in wild-type and Jnk3-null mice. CONCLUSION: Although JNK3 is a key protein involved in cell death in different neurodegenerative diseases, the present study demonstrates that the lack of JNK3 does not confer neuroprotection against 3-NP-induced neuronal death.

Recombination rates across porcine autosomes inferred from high-density linkage maps.

Studies of the variation in recombination rate across the genome provide a better understanding of evolutionary genomics and are also an important step towards mapping and dissecting complex traits in domestic animals. With the recent completion of the porcine genome sequence and the availability of a high-density porcine single nucleotide polymorphism (SNP) array, it is now possible to construct a high-density porcine linkage map and estimate recombination rate across the genome. A total of 416 animals were genotyped with the Porcine SNP60BeadChip, and high-density chromosome linkage maps were constructed using CRI-MAP, assuming the physical order of the Sscrofa10 assembly. The total linkage map length was 2018.79 cM, using 658 meioses and 14,503 SNPs. The estimated average recombination rate across the porcine autosomes was 0.86 cM/Mb. However, a large variation in recombination rate was observed among chromosomes. The estimated average recombination rates (cM/Mb) per chromosome ranged from 0.48 in SSC1 to 1.48 in SSC10, displaying a significant negative correlation with the chromosome sizes. In addition, the analysis of the variation in the recombination rates taking 1-Mb sliding windows has allowed us to demonstrate the variation in recombination rates within chromosomes. In general, a larger recombination rate was observed in the extremes than in the centre of the chromosome. Finally, the ratio between female and male recombination rates was also inferred, obtaining a value of 1.38, with the heterogametic sex having the least recombination.

Joint effects of porcine leptin and leptin receptor polymorphisms on productivity and quality traits.

Leptin signalling plays a fundamental role in growth, fatness and body composition. The aim of this study was to investigate the porcine LEP gene sequence in an Iberian × Landrace experimental cross to identify polymorphisms associated with productivity and quality traits. Because of the documented effects on these traits of the LEPR c.1987C>T polymorphism, the LEP and LEPR c.1987C>T polymorphisms and their interactions have been jointly investigated. The LEP gene sequencing has allowed the identification of 39 polymorphisms, eight of which are novel. Three intronic SNPs, LEP g.1382C>T, LEP g.1387C>T and LEP g.1723A>G, have been genotyped, and association analyses have been carried out. Analyses of LEP g.1387C>T, fully linked to LEP g.1382C>T, have revealed additive effects on live and carcass weights and dominant effects on several backfat thickness measurements. Novel effects of both LEP and LEPR polymorphisms on fatty acid composition in subcutaneous fat have been detected, probably mediated through the effects on fatness. The results reported here suggest that the T alleles of both LEP g.1387C>T and LEPR c.1987C>T, which are fixed in the Iberian pigs, would lead to an increase in growth, fatness and saturated fatty acid content in fat, which could be explained by an increased feed intake.

Evaluation of the porcine ACSL4 gene as a candidate gene for meat quality traits in pigs.

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.

Evolutionary study of a potential selection target region in the pig.

Domestication, modern breeding and artificial selection have shaped dramatically the genomic variability of domestic animals. In livestock, the so-called FAT1 quantitative trait locus (QTL) in porcine chromosome 4 was the first QTL uncovered although, to date, its precise molecular nature has remained elusive. Here, we characterize the nucleotide variability of 13 fragments of ∼500 bp equally spaced in a 2 Mb region in the vicinity of the FAT1 region in a wide-diversity panel of 32 pigs. Asian and European animals, including local Mediterranean and international pig breeds, were sequenced. Patterns of genetic variability were very complex and varied largely across loci and populations; they did not reveal overall a clear signal of a selective sweep in any breed, although FABP4 fragment showed a significantly higher diversity. We used an approximate Bayesian computation approach to infer the evolutionary history of this SSC4 region. Notably, we found that European pig populations have a much lower effective size than their Asian counterparts: in the order of hundreds vs hundreds of thousands. We show also an important part of extant European variability is actually due to introgression of Asian germplasm into Europe. This study shows how a potential loss in diversity caused by bottlenecks and possible selective sweeps associated with domestication and artificial selection can be counterbalanced by migration, making it much more difficult the identification of selection footprints based on naive demographic assumptions. Given the small fragment analyzed here, it remains to be studied how these conclusions apply to the rest of the genome.

Partial short-read sequencing of a highly inbred Iberian pig and genomics inference thereof.

Despite dramatic reduction in sequencing costs with the advent of next generation sequencing technologies, obtaining a complete mammalian genome sequence at sufficient depth is still costly. An alternative is partial sequencing. Here, we have sequenced a reduced representation library of an Iberian sow from the Guadyerbas strain, a highly inbred strain that has been used in numerous QTL studies because of its extreme phenotypic characteristics. Using the Illumina Genome Analyzer II (San Diego, CA, USA), we resequenced ∼ 1% of the genome with average 4 × depth, identifying 68,778 polymorphisms. Of these, 55,457 were putative fixed differences with respect to the assembly, based on the genome of a Duroc pig, and 13,321 were heterozygous positions within Guadyerbas. Despite being highly inbred, the estimate of heterozygosity within Guadyerbas was ∼ 0.78 kb(-1) in autosomes, after correcting for low depth. Nucleotide variability was consistently higher at the telomeric regions than on the rest of the chromosome, likely a result of increased recombination rates. Further, variability was 50% lower in the X-chromosome than in autosomes, which may be explained by a recent bottleneck or by selection. We divided the whole genome in 500 kb windows and we analyzed overrepresented gene ontology terms in regions of low and high variability. Multi organism process, pigmentation and cell killing were overrepresented in high variability regions and metabolic process ontology, within low variability regions. Further, a genome wide Hudson-Kreitman-Aguadé test was carried out per window; overall, variability was in agreement with neutral expectations.

Nucleotide variability of the porcine SERPINA6 gene and the origin of a putative causal mutation associated with meat quality.

The serpin peptidase inhibitor, clade A, member 6 gene (SERPINA6), also known as corticosteroid-binding globulin or CBG, is involved in obesity and stress sensitivity. Previous studies have reported putative causal mutations within that gene in the porcine species. To characterize a hypothetical selective footprint, we have resequenced approximately 6 kb of coding and non-coding fragments in 20 pigs comprising domestic breeds and wild boars from Asia and Europe. Nucleotide variability was found to be far greater within Asian pig breeds than European breeds (π = 1% vs. 0.05%, respectively), which is consistent with pig evolutionary history. The putative causal amino acid substitution p.Gly307Arg (SNP c.919G>A) associated with meat quality (drip loss) was only detected in European domestic pig breeds, suggesting a very recent mutation that appeared after domestication in Europe. No support for positive selection was detected, as no reduction in levels of diversity surrounding the mutation was found in lean breeds with respect to wild boar.

A suitable duplex PCR for ovine embryo sex and genotype of PrnP gene determination for MOET-based selection programmes.

The objective of this study was to test the suitability of a duplex PCR assay for sex and scrapie resistance genotype determination in fresh embryos. Duplex PCR amplified a repetitive and specific fragment of Y chromosome, used for sex diagnosis, and a PrnP fragment. PrnP codons 134 and 156, and codon 171 were genotyped by restriction fragment length polymorphisms and allele-specific PCR, respectively, after re-amplification of PrnP fragment. The specificity of the method was first assessed by testing 359 blood samples from Rasa Aragonesa sheep breed (161 males and 198 females). No amplification failures and total agreement between genotypic and phenotypic sex were found. In the same way, PrnP genotype determination by duplex PCR assay was in agreement with the PrnP animal's genotype established by sequencing. Finally, 73 samples of 1-10 cells from compact morulae were aspirated through the zona pellucida and genotyped for sex and PrnP. The efficiency was 96% when three or more cells were sampled. These results confirm that the duplex PCR assay reported in this work can be used for rapid sex determination in ovine embryos, with a high efficiency and accuracy (96%) when three or more cells are sampled, allowing sexed fresh embryos of known PrnP genotype to be transferred in multiple ovulation and embryo transfer programmes.

The efficiency of in vitro ovine embryo production using an undefined or a defined maturation medium is determined by the source of the oocyte.

In vitro oocyte maturation can be influenced by oocyte source and maturation media composition. The aim of the present study was to compare the efficiency of a defined in vitro maturation medium (TCM199 supplemented with cysteamine and epidermal growth factor; Cys + EGF) with an undefined medium (TCM199 supplemented with follicle-stimulating hormone and follicular fluid; FSH + FF) for in vitro production (IVP) of ovine embryos, using oocytes obtained by laparoscopic ovum pick-up from FSH-stimulated [n=11; 158 cumulus-oocyte complexes (COCs)] and non-stimulated (n=16; 120 COCs) live ewes, as well as abattoir-derived oocytes (170 COCs). The produced blastocysts were vitrified and some of them were transferred to synchronized recipients. The best and the worst final yields of embryo IVP observed in this study were obtained using oocytes from FSH-stimulated ewes matured in FSH + FF (41.3%; 33/80) and in Cys + EGF (19.2%; 15/78) medium, respectively (p<0.01). No significant differences between both media were attained in the blastocyst development rate or in the final yield of embryo IVP using oocytes from non-stimulated ewes or abattoir-derived oocytes. The overall in vivo survival rate of the transferred vitrified blastocysts was 13.1% (8/61), without significant differences between oocyte sources or maturation media. In conclusion, under the experimental conditions of the present study, TCM199 supplemented with cysteamine and EGF is a convenient defined maturation medium for IVP of embryos from oocytes of live non-stimulated ewes or from oocytes of abattoir-derived ovaries. However, the best final yield of embryo IVP observed in this study was attained when oocytes came from FSH-stimulated donors and TCM199 was supplemented with FSH and follicular fluid.

Building phenotype networks to improve QTL detection: a comparative analysis of fatty acid and fat traits in pigs.

Models in QTL mapping can be improved by considering all potential variables, i.e. we can use remaining traits other than the trait under study as potential predictors. QTL mapping is often conducted by correcting for a few fixed effects or covariates (e.g. sex, age), although many traits with potential causal relationships between them are recorded. In this work, we evaluate by simulation several procedures to identify optimum models in QTL scans: forward selection, undirected dependency graph and QTL-directed dependency graph (QDG). The latter, QDG, performed better in terms of power and false discovery rate and was applied to fatty acid (FA) composition and fat deposition traits in two pig F2 crosses from China and Spain. Compared with the typical QTL mapping, QDG approach revealed several new QTL. To the contrary, several FA QTL on chromosome 4 (e.g. Palmitic, C16:0; Stearic, C18:0) detected by typical mapping vanished after adjusting for phenotypic covariates in QDG mapping. This suggests that the QTL detected in typical mapping could be indirect. When a QTL is supported by both approaches, there is an increased confidence that the QTL have a primary effect on the corresponding trait. An example is a QTL for C16:1 on chromosome 8. In conclusion, mapping QTL based on causal phenotypic networks can increase power and help to make more biologically sound hypothesis on the genetic architecture of complex traits.

Understanding host-microbiota interactions in the commercial piglet around weaning.

Weaning is a critical period in the life of pigs with repercussions on their health and welfare and on the economy of the swine industry. This study aimed to assess the effect of the commercial early weaning on gut microbiota, intestinal gene expression and serum metabolomic response via an integrated-omic approach combining 16S rRNA gene sequencing, the OpenArray gene expression technology and 1H-NMR spectroscopy. Fourteen piglets from different litters were sampled for blood, jejunum tissue and caecal content two days before (- 2d), and three days after (+ 3d) weaning. A clearly differential ordination of caecal microbiota was observed. Higher abundances of Roseburia, Ruminococcus, Coprococcus, Dorea and Lachnospira genera in weaned piglets compared to prior to weaning showed the quick microbial changes of the piglets' gut microbiota. Downregulation of OCLN, CLDN4, MUC2, MUC13, SLC15A1 and SLC13A1 genes, also evidenced the negative impact of weaning on gut barrier and digestive functions. Metabolomic approach pinpointed significant decreases in choline, LDL, triglycerides, fatty acids, alanine and isoleucine and increases in 3-hydroxybutyrate after weaning. Moreover, the correlation between microbiota and metabolome datasets revealed the existence of metabolic clusters interrelated to different bacterial clusters. Our results demonstrate the impact of weaning stress on the piglet and give insights regarding the associations between gut microbiota and the animal gene activity and metabolic response.

Early socialization and environmental enrichment of lactating piglets affects the caecal microbiota and metabolomic response after weaning.

The aim of this study was to determine the possible impact of early socialization and an enriched neonatal environment to improve adaptation of piglets to weaning. We hypothesized that changes in the microbiota colonization process and in their metabolic response and intestinal functionality could help the animals face weaning stress. A total of 48 sows and their litters were allotted into a control (CTR) or an enriched treatment (ENR), in which piglets from two adjacent pens were combined and enriched with toys. The pattern of caecal microbial colonization, the jejunal gene expression, the serum metabolome and the intestinal physiology of the piglets were assessed before (-2 d) and after weaning (+ 3d). A differential ordination of caecal microbiota was observed after weaning. Serum metabolome suggested a reduced energetic metabolism in ENR animals, as evidenced by shifts in triglycerides and fatty acids, VLDL/LDL and creatine regions. The TLR2 gene showed to be downregulated in the jejunum of ENR pigs after weaning. The integration of gene expression, metabolome and microbiota datasets confirmed that differences between barren and enriched neonatal environments were evident only after weaning. Our results suggest that improvements in adaptation to weaning could be mediated by a better response to the post-weaning stress.