Evaluation of the porcineACSL4gene as a candidate gene for meat quality traits in pigs

Full-Length Transcriptome and Gene Expression Analysis of Different Ovis aries Adipose Tissues Reveals Transcript Variants Involved in Lipid Biosynthesis

Sheep have historically been bred globally as a vital food source. To explore the transcriptome of adipose tissue and investigate key genes regulating adipose metabolism in sheep, adipose tissue samples were obtained from F1 Dorper × Hu sheep. High-throughput sequencing libraries for second- and third-generation sequencing were constructed using extracted total RNA. Functional annotation of differentially expressed genes and isoforms facilitated the identification of key regulatory genes and isoforms associated with sheep fat metabolism. SMRT-seq generated 919,259 high-accuracy cDNA sequences after filtering. Full-length sequences were corrected using RNA-seq sequences, and 699,680 high-quality full-length non-chimeric (FLNC) reads were obtained. Upon evaluating the ratio of total lengths based on FLNC sequencing, it was determined that 36,909 out of 56,316 multiple-exon isoforms met the criteria for full-length status. This indicates the identification of 330,375 full-length FLNC transcripts among the 370,114 multiple-exon FLNC transcripts. By comparing the reference genomes, 60,276 loci and 111,302 isoforms were identified. In addition, 43,423 new genes and 44,563 new isoforms were identified. The results identified 185 (3198), 394 (3592), and 83 (3286) differentially expressed genes (transcripts) between tail and subcutaneous, tail and visceral, and subcutaneous and visceral adipose tissues, respectively. Functional annotation and pathway analysis revealed the following observations. (1) Among the differentially expressed genes (DEGs) of TF and SF tissues, the downregulation of ACADL, ACSL6, and NC_056060.1.2536 was observed in SF, while FFAR4 exhibited upregulation. (2) Among the DEGs of TF and VF tissues, expressions of ACADL, ACSL6, COL1A1, COL1A2, and SCD were downregulated in VF, with upregulation of FFAR4. (3) Among SF and VF expressions of COL1A1, COL1A2, and NC_056060.1.2536 were downregulated in VF. Specific differentially expressed genes (ACADL, ACSL6, COL1A1, COL1A2, FFAR4, NC_056060.1.2536, and SCD) and transcripts (NC_056066.1.1866.16 and NC_056066.1.1866.22) were identified as relevant to fat metabolism. These results provide a dataset for further verification of the regulatory pathway associated with fat metabolism in sheep.

Acyl-CoA synthase ACSL4: an essential target in ferroptosis and fatty acid metabolism

ABSTRACT

Long-chain acyl-coenzyme A (CoA) synthase 4 (ACSL4) is an enzyme that esterifies CoA into specific polyunsaturated fatty acids, such as arachidonic acid and adrenic acid. Based on accumulated evidence, the ACSL4-catalyzed biosynthesis of arachidonoyl-CoA contributes to the execution of ferroptosis by triggering phospholipid peroxidation. Ferroptosis is a type of programmed cell death caused by iron-dependent peroxidation of lipids; ACSL4 and glutathione peroxidase 4 positively and negatively regulate ferroptosis, respectively. In addition, ACSL4 is an essential regulator of fatty acid (FA) metabolism. ACSL4 remodels the phospholipid composition of cell membranes, regulates steroidogenesis, and balances eicosanoid biosynthesis. In addition, ACSL4-mediated metabolic reprogramming and antitumor immunity have attracted much attention in cancer biology. Because it facilitates the cross-talk between ferroptosis and FA metabolism, ACSL4 is also a research hotspot in metabolic diseases and ischemia/reperfusion injuries. In this review, we focus on the structure, biological function, and unique role of ASCL4 in various human diseases. Finally, we propose that ACSL4 might be a potential therapeutic target.

SNP discovery and association study for growth, fatness and meat quality traits in Iberian crossbred pigs

Iberian pigs and its crosses are produced to obtain high-quality meat products. The objective of this work was to evaluate a wide panel of DNA markers, selected by biological and functional criteria, for association with traits related to muscle growth, fatness, meat quality and metabolism. We used 18 crossbred Iberian pigs with divergent postnatal growth patterns for whole genome sequencing and SNP discovery, with over 13 million variants being detected. We selected 1023 missense SNPs located on annotated genes and showing different allele frequencies between pigs with makerdly different growth patterns. We complemented this panel with 192 candidate SNPs obtained from literature mining and from muscle RNAseq data. The selected markers were genotyped in 480 Iberian × Duroc pigs from a commercial population, in which phenotypes were obtained, and an association study was performed for the 1005 successfully genotyped SNPs showing segregation. The results confirmed the effects of several known SNPs in candidate genes (such as LEPR, ACACA, FTO, LIPE or SCD on fatness, growth and fatty acid composition) and also disclosed interesting effects of new SNPs in less known genes such as LRIG3, DENND1B, SOWAHB, EPHX1 or NFE2L2 affecting body weight, average daily gain and adiposity at different ages, or KRT10, NLE1, KCNH2 or AHNAK affecting fatness and FA composition. The results provide a valuable basis for future implementation of marker-assisted selection strategies in swine and contribute to a better understanding of the genetic architecture of relevant traits.

MiR-23b Promotes Porcine Preadipocyte Differentiation via SESN3 and ACSL4

Genetic improvement of pork quality is one of the hot topics in pig germplasm innovation. Backfat thickness and intramuscular fat content are important indexes of meat quality. MiRNAs are becoming recognized as a crucial regulator of adipose development. Therefore, it is crucial to understand how miR-23b regulates fat metabolism at the molecular level. In the present study, Oil Red O staining, and Western blot were used to evaluate the effect of miR-23b on the differentiation of porcine preadipocytes. Dual-luciferase reporter gene assay, pulldown, and RIP were used to reveal the mechanism of miR-23b regulating cell differentiation. The findings demonstrated that miR-23b promotes the expression of adipogenic factors and increases the content of lipid droplets, thus promoting the differentiation of preadipocytes. Further research found that miR-23b can directly bind to the 3’UTR of SESN3 to regulate adipogenic differentiation. In addition, it was speculated that miR-23b controls cell differentiation by positively regulating the expression of ACSL4 in other ways. Here, we demonstrate that miR-23b promotes the differentiation of porcine preadipocytes by targeting SESN3 and promoting the expression of ACSL4. The present study is meaningful to the improvement of pork quality and the development of animal husbandry.

ACSL4 Directs Intramuscular Adipogenesis and Fatty Acid Composition in Pigs

Simple Summary

In the livestock industry, intramuscular fat content is an important indicator of the meat quality of domestic animals. The variations of the Acyl-CoA Synthetase Long-Chain Family Member 4 (ACSL4) gene locus are associated with intramuscular fat content in different pig populations, but the detailed molecular function of ACSL4 in pig intramuscular adipogenesis remains obscure. Our study reveals the function of ACSL4 in pig intramuscular adipogenesis and provides new clues for improving the palatability of meat and enhancing the nutritional value of pork for human health.

Abstract

The intramuscular fat is a major quality trait of meat, affecting sensory attributes such as flavor and texture. Several previous GWAS studies identified Acyl-CoA Synthetase Long Chain Family Member 4 (ACSL4) gene as the candidate gene to regulate intramuscular fat content in different pig populations, but the underlying molecular function of ACSL4 in adipogenesis within pig skeletal muscle is not fully investigated. In this study, we isolated porcine endogenous intramuscular adipocyte progenitors and performed ACSL4 loss- and gain-of-function experiments during adipogenic differentiation. Our data showed that ACSL4 is a positive regulator of adipogenesis in intramuscular fat cells isolated from pigs. More interestingly, the enhanced expression of ACSL4 in pig intramuscular adipocytes could increase the cellular content of monounsaturated and polyunsaturated fatty acids, such as gamma-L eicosapentaenoic acid (EPA), docosahexaenoic acid (DHA). The above results not only confirmed the function of ACSL4 in pig intramuscular adipogenesis and meat quality attributes, but also provided new clues for the improvement of the nutritional value of pork for human health.

Biological mechanisms of growth performance and meat quality in porcine muscle tissue

Growth performance and meat quality are important traits for pig production. The aim of the present study was to investigate the molecular mechanisms underlying growth performance and meat quality, and to identify novel target molecules for predicting the growth performance and meat quality. The differentially expressed genes (DEGs) in Diannan small ears pigs (DSP) and Landrace pigs (LP) were assessed by RNA-sequencing analyzing technology. A total of 339 DEGs were obtained between DSP and LP. 146 DEGs were upregulated in LP compared with DSP and 193 DEGs were upregulated in DSP compared with LP. The DEGs were significantly enriched in 26 GO and 3 KEGG pathways. The protein-protein interaction (PPI) network with 201 nodes and 382 edges was constructed and 5 modules were extracted from the entire network. The identified upregulated expression of genes involved in glycolysis and myogenesis as well as extracellular matrix may be associated with fast body and muscle deposition rates in LP. Increased expression of genes involved in PPAR signaling pathway and fatty acid metabolism as well as oxidative phosphate processes could be related to the intramuscular fat deposition and meat quality in DSP. The present study may provide an improved understanding of the growth performance and meat quality.

Comparative Transcriptome Profile between Iberian Pig Varieties Provides New Insights into Their Distinct Fat Deposition and Fatty Acids Content

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.

Quercetin Improving Lipid Metabolism by Regulating Lipid Metabolism Pathway of Ileum Mucosa in Broilers

This study is aimed at evaluating the regulatory mechanism of quercetin on lipid metabolism in the ileum of broilers to better understand these pathways decreasing abdominal fat. 480 chickens were randomly divided into 4 groups (control, 0.02% quercetin, 0.04% quercetin, and 0.06% quercetin). Breast muscle, thigh muscle, and abdominal fat pad were removed and weighed at 42 d of age. Serum was obtained by centrifuging blood samples from the jugular vein (10 ml) to determine high-density lipoprotein (HDL), total cholesterol (TC), low-density lipoprotein (LDL), triglyceride (TG), leptin, and adiponectin using ELISA. About 5 g of the ileum was harvested and immediately frozen in liquid nitrogen for RNA-seq. Then, the confirmation of RNA-seq results by the Real-Time Quantitative PCR (RT-qPCR) method was evaluated using Pearson's correlation. Compared with control, abdominal fat percentage was significantly decreased with increasing quercetin supplementation, and the best result was obtained at 0.06% dietary quercetin supplementation (P < 0.01). Breast muscle percentage was significantly decreased at 0.02% quercetin (P < 0.01), and thigh muscle percentage tended to increase (P = 0.078). Meanwhile, 0.04% and 0.06% quercetin significantly decreased TG (P < 0.01), TC (P < 0.01), and LDL content (P < 0.05) in serum. Serum leptin and adiponectin contents were significantly increased by 0.04% and 0.06% dietary quercetin supplementation, compared with the control (P < 0.01). Analyses of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) database were used to identify differently expressed genes and lipid metabolism pathways. Quercetin decreased abdominal fat percentage through regulating fat digestion and absorption, glycerophospholipid metabolism, AMPK signaling pathway, fatty acid degradation, and cholesterol metabolism.

Evolution, expression profile, and regulatory characteristics of ACSL gene family in chicken (Gallus gallus)

Long chain acyl-CoA synthetases (ACSLs), which drive the conversion of long chain fatty acid into acyl-CoA, an ingredient of lipid synthesis, have been well-acknowledged to exert an indispensable role in many metabolic processes in mammals, especially lipid metabolism. However, in chicken, the evolutionary characteristics, expression profiles and regulatory mechanisms of ACSL gene family are rarely understood. Here, we analyzed the genomic synteny, gene structure, evolutionary event and functional domains of the ACSL gene family members using bioinformatics methods. The spatiotemporal expression profiles of ACSL gene family, and their regulatory mechanism were investigated via bioinformatics analysis incorporated with in vivo and in vitro estrogen-treated experiments. Our results indicated that ACSL2 gene was indeed evolutionarily lost in the genome of chicken. Chicken ACSLs shared an AMP-binding functional domain, as well as highly conversed ATP/AMP and FACS signature motifs, and were clustered into two clades, ACSL1/5/6 and ACSL3/4, based on high sequence similarity, similar gene features and conversed motifs. Chicken ACSLs showed differential tissue expression distributions, wherein the significantly decreased expression level of ACSL1 and the significantly increased expression level of ACSL5 were found, respectively, the expression levels of the other ACSL members remained unchanged in the liver of peak-laying hens versus pre-laying hens. Moreover, the transcription activity of ACSL1, ACSL3 and ACSL4 was silenced and ACSL6 was activated by estrogen, but no response to ACSL5. In conclusion, though having highly conversed functional domains, chicken ACSL gene family is organized into two separate groups, ACSL1/5/6 and ACSL3/4, and exhibits varying expression profiles and estrogen effects. These results not only pave the way for better understanding the specific functions of ACSL genes in avian lipid metabolism, but also provide a valuable evidence for gene family characteristics.

Diversity across major and candidate genes in European local pig breeds

The aim of this work was to analyse the distribution of causal and candidate mutations associated to relevant productive traits in twenty local European pig breeds. Also, the potential of the SNP panel employed for elucidating the genetic structure and relationships among breeds was evaluated. Most relevant genes and mutations associated with pig morphological, productive, meat quality, reproductive and disease resistance traits were prioritized and analyzed in a maximum of 47 blood samples from each of the breeds (Alentejana, Apulo-Calabrese, Basque, Bísara, Majorcan Black, Black Slavonian (Crna slavonska), Casertana, Cinta Senese, Gascon, Iberian, Krškopolje (Krškopoljski), Lithuanian indigenous wattle, Lithuanian White Old Type, Mora Romagnola, Moravka, Nero Siciliano, Sarda, Schwäbisch-Hällisches Schwein (Swabian Hall pig), Swallow-Bellied Mangalitsa and Turopolje). We successfully analyzed allelic variation in 39 polymorphisms, located in 33 candidate genes. Results provide relevant information regarding genetic diversity and segregation of SNPs associated to production and quality traits. Coat color and morphological trait-genes that show low level of segregation, and fixed SNPs may be useful for traceability. On the other hand, we detected SNPs which may be useful for association studies as well as breeding programs. For instance, we observed predominance of alleles that might be unfavorable for disease resistance and boar taint in most breeds and segregation of many alleles involved in meat quality, fatness and growth traits. Overall, these findings provide a detailed catalogue of segregating candidate SNPs in 20 European local pig breeds that may be useful for traceability purposes, for association studies and for breeding schemes. Population genetic analyses based on these candidate genes are able to uncover some clues regarding the hidden genetic substructure of these populations, as the extreme genetic closeness between Iberian and Alentejana breeds and an uneven admixture of the breeds studied. The results are in agreement with available knowledge regarding breed history and management, although largest panels of neutral markers should be employed to get a deeper understanding of the population's structure and relationships.

Transcriptomics Analysis on Excellent Meat Quality Traits of Skeletal Muscles of the Chinese Indigenous Min Pig Compared with the Large White Breed

The Min pig (Sus scrofa) is a well-known indigenous breed in China. One of its main advantages over European breeds is its high meat quality. Additionally, different cuts of pig also show some different traits of meat quality. To explore the underlying mechanism responsible for the differences of meat quality between different breeds or cuts, the longissimus dorsi muscle (LM) and the biceps femoris muscle (BF) from Min and Large White pigs were investigated using transcriptome analysis. The gene expression profiling identified 1371 differentially expressed genes (DEGs) between LM muscles from Min and Large White pigs, and 114 DEGs between LM and BF muscles from the same Min pigs. Gene Ontology (GO) enrichment of biological functions and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the gene products were mainly involved in the IRS1/Akt/FoxO1 signaling pathway, adenosine 5′-monophosphate-activated protein kinase (AMPK) cascade effects, lipid metabolism and amino acid metabolism pathway. Such pathways contributed to fatty acid metabolism, intramuscular fat deposition, and skeletal muscle growth in Min pig. These results give an insight into the mechanisms underlying the formation of skeletal muscle and provide candidate genes for improving meat quality. It will contribute to improving meat quality of pigs through molecular breeding.

Integration of liver gene co-expression networks and eGWAs analyses highlighted candidate regulators implicated in lipid metabolism in pigs

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.

Analysis of large versus small dogs reveals three genes on the canine X chromosome associated with body weight, muscling and back fat thickness

Domestic dog breeds display significant diversity in both body mass and skeletal size, resulting from intensive selective pressure during the formation and maintenance of modern breeds. While previous studies focused on the identification of alleles that contribute to small skeletal size, little is known about the underlying genetics controlling large size. We first performed a genome-wide association study (GWAS) using the Illumina Canine HD 170,000 single nucleotide polymorphism (SNP) array which compared 165 large-breed dogs from 19 breeds (defined as having a Standard Breed Weight (SBW) >41 kg [90 lb]) to 690 dogs from 69 small breeds (SBW ≤41 kg). We identified two loci on the canine X chromosome that were strongly associated with large body size at 82-84 megabases (Mb) and 101-104 Mb. Analyses of whole genome sequencing (WGS) data from 163 dogs revealed two indels in the Insulin Receptor Substrate 4 (IRS4) gene at 82.2 Mb and two additional mutations, one SNP and one deletion of a single codon, in Immunoglobulin Superfamily member 1 gene (IGSF1) at 102.3 Mb. IRS4 and IGSF1 are members of the GH/IGF1 and thyroid pathways whose roles include determination of body size. We also found one highly associated SNP in the 5'UTR of Acyl-CoA Synthetase Long-chain family member 4 (ACSL4) at 82.9 Mb, a gene which controls the traits of muscling and back fat thickness. We show by analysis of sequencing data from 26 wolves and 959 dogs representing 102 domestic dog breeds that skeletal size and body mass in large dog breeds are strongly associated with variants within IRS4, ACSL4 and IGSF1.

Transcriptomic Analysis Identifies Candidate Genes Related to Intramuscular Fat Deposition and Fatty Acid Composition in the Breast Muscle of Squabs (Columba)

Despite the fact that squab is consumed throughout the world because of its high nutritional value and appreciated sensory attributes, aspects related to its characterization, and in particular genetic issues, have rarely been studied. In this study, meat traits in terms of pH, water-holding capacity, intramuscular fat content, and fatty acid profile of the breast muscle of squabs from two meat pigeon breeds were determined. Breed-specific differences were detected in fat-related traits of intramuscular fat content and fatty acid composition. RNA-Sequencing was applied to compare the transcriptomes of muscle and liver tissues between squabs of two breeds to identify candidate genes associated with the differences in the capacity of fat deposition. A total of 27 differentially expressed genes assigned to pathways of lipid metabolism were identified, of which, six genes belonged to the peroxisome proliferator-activated receptor signaling pathway along with four other genes. Our results confirmed in part previous reports in livestock and provided also a number of genes which had not been related to fat deposition so far. These genes can serve as a basis for further investigations to screen markers closely associated with intramuscular fat content and fatty acid composition in squabs. The data from this study were deposited in the National Center for Biotechnology Information (NCBI)’s Sequence Read Archive under the accession numbers SRX1680021 and SRX1680022. This is the first transcriptome analysis of the muscle and liver tissue in Columba using next generation sequencing technology. Data provided here are of potential value to dissect functional genes influencing fat deposition in squabs.

Polymorphisms and expression analysis of SOX-6 in relation to porcine growth, carcass, and meat quality traits

The aim of the study was to investigate single nucleotide polymorphisms (SNPs) and expression of SOX-6 to support its candidacy for growth, carcass, and meat quality traits in pigs. The first SNP, rs81358375, was associated with pH 45 min post mortem in loin (pH1L), the thickness of backfat and side fat, and carcass length in Pietrain (Pi) population, and related with backfat thickness and daily gain in Duroc × Pietrain F2 (DuPi) population. The other SNP, rs321666676, was associated with meat colour in Pi population. In DuPi population, the protein, not mRNA, level of SOX-6 in high pH1L pigs was significantly less abundant compared with low pH1L pigs, where microRNAs targeting SOX-6 were also differently regulated. This paper shows that SOX-6 could be a potential candidate gene for porcine growth, carcass, and meat quality traits based on genetic association and gene expression.

Analysis of porcine adipose tissue transcriptome reveals differences in de novo fatty acid synthesis in pigs with divergent muscle fatty acid composition

Background

In pigs, adipose tissue is one of the principal organs involved in the regulation of lipid metabolism. It is particularly involved in the overall fatty acid synthesis with consequences in other lipid-target organs such as muscles and the liver. With this in mind, we have used massive, parallel high-throughput sequencing technologies to characterize the porcine adipose tissue transcriptome architecture in six Iberian x Landrace crossbred pigs showing extreme phenotypes for intramuscular fatty acid composition (three per group).

Results

High-throughput RNA sequencing was used to generate a whole characterization of adipose tissue (backfat) transcriptome. A total of 4,130 putative unannotated protein-coding sequences were identified in the 20% of reads which mapped in intergenic regions. Furthermore, 36% of the unmapped reads were represented by interspersed repeats, SINEs being the most abundant elements. Differential expression analyses identified 396 candidate genes among divergent animals for intramuscular fatty acid composition. Sixty-two percent of these genes (247/396) presented higher expression in the group of pigs with higher content of intramuscular SFA and MUFA, while the remaining 149 showed higher expression in the group with higher content of PUFA. Pathway analysis related these genes to biological functions and canonical pathways controlling lipid and fatty acid metabolisms. In concordance with the phenotypic classification of animals, the major metabolic pathway differentially modulated between groups was de novo lipogenesis, the group with more PUFA being the one that showed lower expression of lipogenic genes.

Conclusions

These results will help in the identification of genetic variants at loci that affect fatty acid composition traits. The implications of these results range from the improvement of porcine meat quality traits to the application of the pig as an animal model of human metabolic diseases.

Dissecting structural and nucleotide genome-wide variation in inbred Iberian pigs

BACKGROUND

In contrast to international pig breeds, the Iberian breed has not been admixed with Asian germplasm. This makes it an important model to study both domestication and relevance of Asian genes in the pig. Besides, Iberian pigs exhibit high meat quality as well as appetite and propensity to obesity. Here we provide a genome wide analysis of nucleotide and structural diversity in a reduced representation library from a pool (n=9 sows) and shotgun genomic sequence from a single sow of the highly inbred Guadyerbas strain. In the pool, we applied newly developed tools to account for the peculiarities of these data.

RESULTS

A total of 254,106 SNPs in the pool (79.6 Mb covered) and 643,783 in the Guadyerbas sow (1.47 Gb covered) were called. The nucleotide diversity (1.31x10-3 per bp in autosomes) is very similar to that reported in wild boar. A much lower than expected diversity in the X chromosome was confirmed (1.79x10-4 per bp in the individual and 5.83x10-4 per bp in the pool). A strong (0.70) correlation between recombination and variability was observed, but not with gene density or GC content. Multicopy regions affected about 4% of annotated pig genes in their entirety, and 2% of the genes partially. Genes within the lowest variability windows comprised interferon genes and, in chromosome X, genes involved in behavior like HTR2C or MCEP2. A modified Hudson-Kreitman-Aguadé test for pools also indicated an accelerated evolution in genes involved in behavior, as well as in spermatogenesis and in lipid metabolism.

CONCLUSIONS

This work illustrates the strength of current sequencing technologies to picture a comprehensive landscape of variability in livestock species, and to pinpoint regions containing genes potentially under selection. Among those genes, we report genes involved in behavior, including feeding behavior, and lipid metabolism. The pig X chromosome is an outlier in terms of nucleotide diversity, which suggests selective constraints. Our data further confirm the importance of structural variation in the species, including Iberian pigs, and allowed us to identify new paralogs for known gene families.