Mining the pig genome to investigate the domestication process

Exploring the mechanism of artificial selection signature in Chinese indigenous pigs by leveraging multiple bioinformatics database tools

BACKGROUND

Chinese indigenous pigs in Yunnan exhibit considerable phenotypic diversity, but their population structure and the biological interpretation of signatures of artificial selection require further investigation. To uncover population genetic diversity, migration events, and artificial selection signatures in Chinese domestic pigs, we sampled 111 Yunnan pigs from four breeds in Yunnan which is considered to be one of the centres of livestock domestication in China, and genotyped them using Illumina Porcine SNP60K BeadChip. We then leveraged multiple bioinformatics database tools to further investigate the signatures and associated complex traits.

RESULTS

Population structure and migration analyses showed that Diannanxiaoer pigs had different genetic backgrounds from other Yunnan pigs, and Gaoligongshan may undergone the migration events from Baoshan and Saba pigs. Intriguingly, we identified a possible common target of sharing artificial selection on a 265.09 kb region on chromosome 5 in Yunnan indigenous pigs, and the genes on this region were associated with cardiovascular and immune systems. We also detected several candidate genes correlated with dietary adaptation, body size (e.g., PASCIN1, GRM4, ITPR2), and reproductive performance. In addition, the breed-sharing gene MMP16 was identified to be a human-mediated gene. Multiple lines of evidence at the mammalian genome, transcriptome, and phenome levels further supported the evidence for the causality between MMP16 variants and the metabolic diseases, brain development, and cartilage tissues in Chinese pigs. Our results suggested that the suppression of MMP16 would directly lead to inactivity and insensitivity of neuronal activity and skeletal development in Chinese indigenous pigs.

CONCLUSION

In this study, the population genetic analyses and identification of artificial selection signatures of Yunnan indigenous pigs help to build an understanding of the effect of human-mediated selection mechanisms on phenotypic traits in Chinese indigenous pigs. Further studies are needed to fully characterize the process of human-mediated genes and biological mechanisms.

Identification and analysis of circRNAs in the prefrontal cortices of wild boar and domestic pig

Domestication caused significant differences in morphology and behavior between wild and domestic pigs. However, the regulatory role of circRNA in this event is unclear. Here, we analyzed circRNA expression patterns in the prefrontal cortices of wild boar and domestic pigs to determine the potential role of circRNAs in domestication. We identified a total of 11,375 circRNAs and found that 349 and 354 circRNAs were up-regulated in wild boar and Rongchang pig, respectively. Functional enrichment analysis showed that host genes of significantly highly-expressed circRNAs in wild boar were significantly enriched in neural synapse-related categories and the categories of 'regulation of defense response (p = 0.028)' and 'neural retina development (p = 4.32 × 10-3)'. Host genes of significantly highly-expressed circRNAs in Rongchang pig were specifically involved in 'chordate embryonic development (p = 2.38 × 10-4)'. Additionally, we constructed circRNA-miRNA-mRNA regulatory axes in wild boar and Rongchang pig and found more regulatory axes in wild boar that potentially regulate synaptic activities. We identified multiple circRNAs that may be related to domesticated characteristics, such as ssc_circ_6179 (ssc_circ_6179-ssc-miR-9847-HRH3, related to aggression) and ssc_circ_3027 (ssc_circ_3027-ssc-miR-4334-5p-HCRTR1, related to attention). This study provides a resource for further investigation of the molecular basis of pig domestication.

Anthropogenic hybridization and its influence on the adaptive potential of the Sardinian wild boar (Sus scrofa meridionalis)

The wild boar (Sus scrofa meridionalis) arrived in Sardinia with the first human settlers in the early Neolithic with the potential to hybridize with the domestic pig (S. s. domesticus) throughout its evolution on the island. In this paper, we investigated the possible microevolutionary effects of such introgressive hybridization on the present wild boar population, comparing Sardinian wild specimens with several commercial pig breeds and Sardinian local pigs, along with a putatively unadmixed wild boar population from Central Italy, all genotyped with a medium density SNP chip. We first aimed at identifying hybrids in the population using different approaches, then examined genomic regions enriched for domestic alleles in the hybrid group, and finally we applied two methods to find regions under positive selection to possibly highlight instances of domestic adaptive introgression into a wild population. We found three hybrids within the Sardinian sample (3.1% out of the whole dataset). We reported 11 significant windows under positive selection with a method that looks for overly differentiated loci in the target population, compared with other two populations. We also identified 82 genomic regions with signs of selection in the domestic pig but not in the wild boar, two of which overlapped with genomic regions enriched for domestic alleles in the hybrid pool. Genes in these regions can be linked with reproductive success. Given our results, domestic introgression does not seem to be pervasive in the Sardinian wild boar. Nevertheless, we suggest monitoring the possible spread of advantageous domestic alleles in the coming years.

Genomic diversity and signals of selection processes in wild and farm-reared red-legged partridges (Alectoris rufa)

The genetic dynamics of wild populations with releases of farm-reared reinforcements are very complex. These releases can endanger wild populations through genetic swamping or by displacing them. We assessed the genomic differences between wild and farm-reared red-legged partridges (Alectoris rufa) and described differential selection signals between both populations. We sequenced the whole genome of 30 wild and 30 farm-reared partridges. Both partridges had similar nucleotide diversity (π). Farm-reared partridges had a more negative Tajima's D and more and longer regions of extended haplotype homozygosity than wild partridges. We observed higher inbreeding coefficients (F_IS and F_ROH) in wild partridges. Selective sweeps (Rsb) were enriched with genes that contribute to the reproductive, skin and feather colouring, and behavioural differences between wild and farm-reared partridges. The analysis of genomic diversity should inform future decisions for the preservation of wild populations.

Pigs as Pets: Early Human Relations with the Sulawesi Warty Pig (Sus celebensis)

The Sulawesi warty pig (S. celebensis) is a wild and still-extant suid that is endemic to the Indonesian island of Sulawesi. It has long been theorised that S. celebensis was domesticated and/or deliberately introduced to other islands in Indonesia prior to the advent of the Neolithic farming transition in the region. Thus far, however, there has been no empirical support for this idea, nor have scientists critiqued the argument that S. celebensis was a pre-Neolithic domesticate in detail. Here, it is proposed that early foragers could have formed a relationship with S. celebensis that was similar in essence to the close association between Late Pleistocene foragers in Eurasia and the wild wolf ancestors of domestic dogs. That is, a longstanding practice of hunter-gatherers intensively socialising wild-caught S. celebensis piglets for adoption into human society as companion animals ('pets') may have altered the predator-prey dynamic, brought aspects of wild pig behaviour and reproduction under indirect human selection and control, and caused changes that differentiated human-associated pigs from their solely wild-living counterparts.

Genome-wide association study identifying genetic variants associated with carcass backfat thickness, lean percentage and fat percentage in a four-way crossbred pig population using SLAF-seq technology

Background

Carcass backfat thickness (BFT), carcass lean percentage (CLP) and carcass fat percentage (CFP) are important to the commercial pig industry. Nevertheless, the genetic architecture of BFT, CLP and CFP is still elusive. Here, we performed a genome-wide association study (GWAS) based on specific-locus amplified fragment sequencing (SLAF-seq) to analyze seven fatness-related traits, including five BFTs, CLP, and CFP on 223 four-way crossbred pigs.

Results

A total of 227, 921 highly consistent single nucleotide polymorphisms (SNPs) evenly distributed throughout the genome were used to perform GWAS. Using the mixed linear model (MLM), a total of 20 SNP loci significantly related to these traits were identified on ten Sus scrofa chromosomes (SSC), of which 10 SNPs were located in previously reported quantitative trait loci (QTL) regions. On SSC7, two SNPs (SSC7:29,503,670 and rs1112937671) for average backfat thickness (ABFT) exceeded 1% and 10% Bonferroni genome-wide significance levels, respectively. These two SNP loci were located within an intron region of the COL21A1 gene, which was a protein-coding gene that played an important role in the porcine backfat deposition by affecting extracellular matrix (ECM) remodeling. In addition, based on the other three significant SNPs on SSC7, five candidate genes, ZNF184, ZNF391, HMGA1, GRM4 and NUDT3 were proposed to influence BFT. On SSC9, two SNPs for backfat thickness at 6–7 ribs (67RBFT) and one SNP for CLP were in the same locus region (19 kb interval). These three SNPs were located in the PGM2L1 gene, which encoded a protein that played an indispensable role in glycogen metabolism, glycolysis and gluconeogenesis as a key enzyme. Finally, one significant SNP on SSC14 for CLP was located within the PLBD2 gene, which participated in the lipid catabolic process.

Conclusions

A total of two regions on SSC7 and SSC9 and eight potential candidate genes were found for fatness-related traits in pigs. The results of this GWAS based on SLAF-seq will greatly advance our understanding of the genetic architecture of BFT, CLP, and CFP traits. These identified SNP loci and candidate genes might serve as a biological basis for improving the important fatness-related traits of pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08827-8.

The Boar War: Five Hot Factors Unleashing Boar Expansion and Related Emergency

The recent and ever-growing problem of boar (Sus scrofa forms including wild boar, hybrid and feral pig) expansion is a very complex issue in wildlife management. The damages caused to biodiversity and the economies are addressed in different ways by the various countries, but research is needed to shed light on the causal factors of this emergency before defining a useful collaborative management policy. In this review, we screened more than 280 references published between 1975–2022, identifying and dealing with five hot factors (climate change, human induced habitat modifications, predator regulation on the prey, hybridization with domestic forms, and transfaunation) that could account for the boar expansion and its niche invasion. We also discuss some issues arising from this boar emergency, such as epizootic and zoonotic diseases or the depression of biodiversity. Finally, we provide new insights for the research and the development of management policies.

Variability in porcine microRNA genes and its association with mRNA expression and lipid phenotypes

BACKGROUND

Mature microRNAs (miRNAs) play an important role in repressing the expression of a wide range of mRNAs. The presence of polymorphic sites in miRNA genes and their corresponding 3'UTR binding sites can disrupt canonical conserved miRNA-mRNA pairings, and thus modify gene expression patterns. However, to date such polymorphic sites in miRNA genes and their association with gene expression phenotypes and complex traits are poorly characterized in pigs.

RESULTS

By analyzing whole-genome sequences from 120 pigs and wild boars from Europe and Asia, we identified 285 single nucleotide polymorphisms (SNPs) that map to miRNA loci, and 109,724 SNPs that are located in predicted 7mer-m8 miRNA binding sites within porcine 3'UTR. In porcine miRNA genes, SNP density is reduced compared with their flanking non-miRNA regions. By sequencing the genomes of five Duroc boars, we identified 12 miRNA SNPs that were subsequently genotyped in their offspring (N = 345, Lipgen population). Association analyses of miRNA SNPs with 38 lipid-related traits and hepatic and muscle microarray expression phenotypes recorded in the Lipgen population were performed. The most relevant detected association was between the genotype of the rs319154814 (G/A) SNP located in the apical loop of the ssc-miR-326 hairpin precursor and PPP1CC mRNA levels in the liver (q-value = 0.058). This result was subsequently confirmed by qPCR (P-value = 0.027). The rs319154814 (G/A) genotype was also associated with several fatty acid composition traits.

CONCLUSIONS

Our findings show a reduced variability of porcine miRNA genes, which is consistent with strong purifying selection, particularly in the seed region that plays a critical role in miRNA binding. Although it is generally assumed that SNPs mapping to the seed region are those with the most pronounced consequences on mRNA expression, we show that a SNP mapping to the apical region of ssc-miR-326 is significantly associated with hepatic mRNA levels of the PPP1CC gene, one of its predicted targets. Although experimental confirmation of such an interaction is reported in humans but not in pigs, this result highlights the need to further investigate the functional effects of miRNA polymorphisms that are located outside the seed region on gene expression in pigs.

Characterization of a haplotype-reference panel for genotyping by low-pass sequencing in Swiss Large White pigs

BACKGROUND

The key-ancestor approach has been frequently applied to prioritize individuals for whole-genome sequencing based on their marginal genetic contribution to current populations. Using this approach, we selected 70 key ancestors from two lines of the Swiss Large White breed that have been selected divergently for fertility and fattening traits and sequenced their genomes with short paired-end reads.

RESULTS

Using pedigree records, we estimated the effective population size of the dam and sire line to 72 and 44, respectively. In order to assess sequence variation in both lines, we sequenced the genomes of 70 boars at an average coverage of 16.69-fold. The boars explained 87.95 and 95.35% of the genetic diversity of the breeding populations of the dam and sire line, respectively. Reference-guided variant discovery using the GATK revealed 26,862,369 polymorphic sites. Principal component, admixture and fixation index (F_ST) analyses indicated considerable genetic differentiation between the lines. Genomic inbreeding quantified using runs of homozygosity was higher in the sire than dam line (0.28 vs 0.26). Using two complementary approaches, we detected 51 signatures of selection. However, only six signatures of selection overlapped between both lines. We used the sequenced haplotypes of the 70 key ancestors as a reference panel to call 22,618,811 genotypes in 175 pigs that had been sequenced at very low coverage (1.11-fold) using the GLIMPSE software. The genotype concordance, non-reference sensitivity and non-reference discrepancy between thus inferred and Illumina PorcineSNP60 BeadChip-called genotypes was 97.60, 98.73 and 3.24%, respectively. The low-pass sequencing-derived genomic relationship coefficients were highly correlated (r > 0.99) with those obtained from microarray genotyping.

CONCLUSIONS

We assessed genetic diversity within and between two lines of the Swiss Large White pig breed. Our analyses revealed considerable differentiation, even though the split into two populations occurred only few generations ago. The sequenced haplotypes of the key ancestor animals enabled us to implement genotyping by low-pass sequencing which offers an intriguing cost-effective approach to increase the variant density over current array-based genotyping by more than 350-fold.

Selective sweep analysis reveals extensive parallel selection traits between large white and Duroc pigs

In the process of pig genetic improvement, different commercial breeds have been bred for the same purpose, improving meat production. Most of the economic traits, such as growth and fertility, have been selected similarly despite the discrepant selection pressure, which is known as parallel selection. Here, 28 whole-genome sequencing data of Danish large white pigs with an approximately 25-fold depth each were generated, resulting in about 12 million high-quality SNPs for each individual. Combined with the sequencing data of 27 Duroc and 23 European wild boars, we investigated the parallel selection of Danish large white and Duroc pigs using two complementary methods, Fst and iHS. In total, 67 candidate regions were identified as the signatures of parallel selection. The genes in candidate regions of parallel selection were mainly associated with sensory perception, growth rate, and body size. Further functional annotation suggested that the striking consistency of the terms may be caused by the polygenetic basis of quantitative traits, and revealing the complex genetic basis of parallel selection. Besides, some unique terms were enriched in population-specific selection regions, such as the limb development-related terms enriched in Duroc-specific selection regions, suggesting unique selections of breed-specific selected traits. These results will help us better understand the parallel selection process of different breeds. Moreover, we identified several potential causal SNPs that may contribute to the pig genetic breeding process.

Evaluation of genetic richness among Vietnamese native pig breeds using microsatellite markers

We have elucidated genetic relationships of Vietnamese native pigs (VNP) using preliminarily collected samples by a single-nucleotide polymorphism (SNP) array. In order to confirm our previous results and compare with the results of a previous study using microsatellite (MS) markers, we aimed to characterize genetic diversity and population structure in wider varieties (24 breeds from 21 Provinces) of VNP across the country using 20 polymorphic MS markers recommended by ISAG/FAO (International Society for Animal Genetics/Food and Agriculture Organization) for diversity study. In this study, we collected 1,136 DNA samples of the VNPs and three exotic breeds. Our results revealed that the average number of alleles and allelic richness across the loci in VNPs were 10.0 and 7.6, which were higher than those of exotic breeds. Genomic components among VNPs were subjected to the sampling locations. Interestingly, Co Binh Thuan showed remarkable genetic feature compared to the other VNPs, because the habitation of Co Binh Thuan was relatively far from the other breeds. The results of this study provided useful information for exploitation, conservation, and development trends of the VNP breeds. More recently, African swine fever caused significant damage to most of the VNP populations. Therefore, our findings will help a reconstruction scheme of the VNP genetic resources.

Genomic diversity, linkage disequilibrium and selection signatures in European local pig breeds assessed with a high density SNP chip

Genetic characterization of local breeds is essential to preserve their genomic variability, to advance conservation policies and to contribute to their promotion and sustainability. Genomic diversity of twenty European local pig breeds and a small sample of Spanish wild pigs was assessed using high density SNP chips. A total of 992 DNA samples were analyzed with the GeneSeek Genomic Profiler (GGP) 70 K HD porcine genotyping chip. Genotype data was employed to compute genetic diversity, population differentiation and structure, genetic distances, linkage disequilibrium and effective population size. Our results point out several breeds, such as Turopolje, Apulo Calabrese, Casertana, Mora Romagnola and Lithuanian indigenous wattle, having the lowest genetic diversity, supported by low heterozygosity and very small effective population size, demonstrating the need of enhanced conservation strategies. Principal components analysis showed the clustering of the individuals of the same breed, with few breeds being clearly isolated from the rest. Several breeds were partially overlapped, suggesting genetic closeness, which was particularly marked in the case of Iberian and Alentejana breeds. Spanish wild boar was also narrowly related to other western populations, in agreement with recurrent admixture between wild and domestic animals. We also searched across the genome for loci under diversifying selection based on F_ST outlier tests. Candidate genes that may underlie differences in adaptation to specific environments and productive systems and phenotypic traits were detected in potentially selected genomic regions.

Genomic Analysis Reveals Specific Patterns of Homozygosity and Heterozygosity in Inbred Pigs

The inbred strain of miniature pig is an ideal model for biomedical research due to its high level of homozygosity. In this study, we investigated genetic diversity, relatedness, homozygosity, and heterozygosity using the Porcine SNP60K BeadChip in both inbred and non-inbred Wuzhishan pigs (WZSPs). Our results from multidimensional scaling, admixture, and phylogenetic analyses indicated that the inbred WZSP, with its unique genetic properties, can be utilized as a novel genetic resource for pig genome studies. Inbreeding depression and run of homozygosity (ROH) analyses revealed an average of 61 and 12 ROH regions in the inbred and non-inbred genomes of WZSPs, respectively. By investigating ROH number, length, and distribution across generations, we further briefly studied the impacts of recombination and demography on ROH in these WZSPs. Finally, we explored the SNPs with higher heterozygosity across generations and their potential functional implications in the inbred WZSP. We detected 56 SNPs showing constant heterozygosity with He = 1 across six generations in inbred pigs, while only one was found in the non-inbred population. Among these SNPs, we observed nine SNPs located in swine RefSeq genes, which were found to be involved in signaling and immune processes. Together, our findings indicate that the inbred-specific pattern of homozygosity and heterozygosity in inbred pigs can offer valuable insights for elucidating the mechanisms of inbreeding in farm animals.

Signatures of de-domestication in autochthonous pig breeds and of domestication in wild boar populations from MC1R and NR6A1 allele distribution

Autochthonous pig breeds are usually reared in extensive or semi-extensive production systems that might facilitate contact with wild boars and, thus, reciprocal genetic exchanges. In this study, we analysed variants in the melanocortin 1 receptor (MC1R) gene (which cause different coat colour phenotypes) and in the nuclear receptor subfamily 6 group A member 1 (NR6A1) gene (associated with increased vertebral number) in 712 pigs of 12 local pig breeds raised in Italy (Apulo-Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano and Sarda) and south-eastern European countries (Krškopolje from Slovenia, Black Slavonian and Turopolje from Croatia, Mangalitsa and Moravka from Serbia and East Balkan Swine from Bulgaria) and compared the data with the genetic variability at these loci investigated in 229 wild boars from populations spread in the same macro-geographic areas. None of the autochthonous pig breeds or wild boar populations were fixed for one allele at both loci. Domestic and wild-type alleles at these two genes were present in both domestic and wild populations. Findings of the distribution of MC1R alleles might be useful for tracing back the complex genetic history of autochthonous breeds. Altogether, these results indirectly demonstrate that bidirectional introgression of wild and domestic alleles is derived and affected by the human and naturally driven evolutionary forces that are shaping the Sus scrofa genome: autochthonous breeds are experiencing a sort of 'de-domestication' process, and wild resources are challenged by a 'domestication' drift. Both need to be further investigated and managed.

Genome-wide epigenetic landscape of pig lincRNAs and their evolution during porcine domestication

AIM

We aimed to identify previously unreported long intergenic noncoding RNAs (lincRNAs) in the porcine liver, an important metabolic tissue, and further illustrate the epigenomic landscapes and the evolution of lincRNAs.

MATERIALS & METHODS

We used porcine omics data and comprehensively analyzed and identified lincRNAs and their methylation, expression and evolutionary patterns during pig domestication.

RESULTS

LincRNAs exhibit highly methylated promoter and downstream regions, as well as lower expression levels and higher tissue specificity than protein-coding genes. We identified a batch of lincRNAs with selection signals that are associated with pig domestication, which are more highly expressed in the liver than in other tissues (19:10/8/6/3/2/1/1). Interestingly, the lincRNA linc-sscg1779 and its target gene C6, which is crucial in liver metabolism, are differentially expressed during pig domestication.

CONCLUSION

Although they may originate from noisy transcripts, lincRNAs may be subjected to artificial selection. This phenomenon implies the functional importance of lincRNAs in pig domestication.

The full mitochondrial genomes of Mangalica pig breeds and their possible origin

The mitogenomes of one animal of each of the three Mangalica breeds, Blonde, Red, and Swallow-belly were assembled from reads obtained by Next Generation Sequencing of the three genomes. Features of the mitogenomes were identical in the three breeds, apart from a second tRNA-Val gene on the L strand in Swallow-belly. Phylogenetic comparison of the three mitogenomes with 112 full mtDNA sequences clearly put Mangalicas into the European clade. Comparing the mitogenome of eight Mangalica animals revealed particular differences between them. The mitogenome of some Mangalicas was closely related to the Croatian Turopolje breed and this indicates either the common origin of their maternal lineages or admixture of some populations of the breeds. However, the origin of the mitogenome of certain purebred Mangalicas kept in the Hungarian Mangalica Gene Reserve still remains unknown.

Convergent and divergent genetic changes in the genome of Chinese and European pigs

Since 10,000 BC, continuous human selection has led to intense genetic and phenotypic changes in pig (Sus scrofa) domestication. Through whole genome analysis of 257 individuals, we demonstrated artificial unidirectional and bidirectional selection as the primary force to shape the convergent and divergent changes between Chinese domestic pigs (CHD) and European domestic pigs (EUD). We identified 31 genes in unidirectional selection regions that might be related to fundamental domestication requirements in pigs. And these genes belong predominantly to categories related to the nervous system, muscle development, and especially to metabolic diseases. In addition, 35 genes, representing different breeding preference, were found under bidirectional selection for the distinct leanness and reproduction traits between CHD and EUD. The convergent genetic changes, contributing physical and morphological adaption, represent the common concerns on pig domestication. And the divergent genetic changes reflect distinct breeding goals between Chinese and European pigs. Using ITPR3, AHR and NMU as examples, we explored and validated how the genetic variations contribute to the phenotype changes.

A Pathway-Centered Analysis of Pig Domestication and Breeding in Eurasia

Ascertaining the molecular and physiological basis of domestication and breeding is an active area of research. Due to the current wide distribution of its wild ancestor, the wild boar, the pig (Sus scrofa) is an excellent model to study these processes, which occurred independently in East Asia and Europe ca. 9000 yr ago. Analyzing genome variability patterns in terms of metabolic pathways is attractive since it considers the impact of interrelated functions of genes, in contrast to genome-wide scans that treat genes or genome windows in isolation. To that end, we studied 40 wild boars and 123 domestic pig genomes from Asia and Europe when metabolic pathway was the unit of analysis. We computed statistical significance for differentiation (Fst) and linkage disequilibrium (nSL) statistics at the pathway level. In terms of Fst, we found 21 and 12 pathways significantly differentiated at a q-value < 0.05 in Asia and Europe, respectively; five were shared across continents. In Asia, we found six significant pathways related to behavior, which involved essential neurotransmitters like dopamine and serotonin. Several significant pathways were interrelated and shared a variable percentage of genes. There were 12 genes present in >10 significant pathways (in terms of Fst), comprising genes involved in the transduction of a large number of signals, like phospholipase PCLB1, which is expressed in the brain, or ITPR3, which has an important role in taste transduction. In terms of nSL, significant pathways were mainly related to reproductive performance (ovarian steroidogenesis), a similarly important target trait during domestication and modern animal breeding. Different levels of recombination cannot explain these results, since we found no correlation between Fst and recombination rate. However, we did find an increased ratio of deleterious mutations in domestic vs. wild populations, suggesting a relaxed functional constraint associated with the domestication and breeding processes. Purifying selection was, nevertheless, stronger in significantly differentiated pathways than in random pathways, mainly in Europe. We conclude that pathway analysis facilitates the biological interpretation of genome-wide studies. Notably, in the case of pig, behavior played an important role, among other physiological and developmental processes.

A deep catalog of autosomal single nucleotide variation in the pig

A comprehensive catalog of variability in a given species is useful for many important purposes, e.g., designing high density arrays or pinpointing potential mutations of economic or physiological interest. Here we provide a genomewide, worldwide catalog of single nucleotide variants by simultaneously analyzing the shotgun sequence of 128 pigs and five suid outgroups. Despite the high SNP missing rate of some individuals (up to 88%), we retrieved over 48 million high quality variants. Of them, we were able to assess the ancestral allele of more than 39M biallelic SNPs. We found SNPs in 21,455 out of the 25,322 annotated genes in pig assembly 10.2. The annotation showed that more than 40% of the variants were novel variants, not present in dbSNP. Surprisingly, we found a large variability in transition / transversion rate along the genome, which is very well explained (R²=0.79) primarily by genome differences in in CpG content and recombination rate. The number of SNPs per window also varied but was less dependent of known factors such as gene density, missing rate or recombination (R²=0.48). When we divided the samples in four groups, Asian wild boar (ASWB), Asian domestics (ASDM), European wild boar (EUWB) and European domestics (EUDM), we found a marked correlation in allele frequencies between domestics and wild boars within Asia and within Europe, but not across continents, due to the large evolutive distance between pigs of both continents (~1.2 MYA). In general, the porcine species showed a small percentage of SNPs exclusive of each population group. EUWB and EUDM were predicted to harbor a larger fraction of potentially deleterious mutations, according to the SIFT algorithm, than Asian samples, perhaps a result of background selection being less effective due to a lower effective population size in Europe.

Genome-wide association analyses reveal significant loci and strong candidate genes for growth and fatness traits in two pig populations

BACKGROUND

Recently, genome-wide association studies (GWAS) have been reported on various pig traits. We performed a GWAS to analyze 22 traits related to growth and fatness on two pig populations: a White Duroc × Erhualian F2 intercross population and a Chinese Sutai half-sib population.

RESULTS

We identified 14 and 39 loci that displayed significant associations with growth and fatness traits at the genome-wide level and chromosome-wide level, respectively. The strongest association was between a 750 kb region on SSC7 (SSC for Sus scrofa) and backfat thickness at the first rib. This region had pleiotropic effects on both fatness and growth traits in F2 animals and contained a promising candidate gene HMGA1 (high mobility group AT-hook 1). Unexpectedly, population genetic analysis revealed that the allele at this locus that reduces fatness and increases growth is derived from Chinese indigenous pigs and segregates in multiple Chinese breeds. The second strongest association was between the region around 82.85 Mb on SSC4 and average backfat thickness. PLAG1 (pleiomorphic adenoma gene 1), a gene under strong selection in European domestic pigs, is proximal to the top SNP and stands out as a strong candidate gene. On SSC2, a locus that significantly affects fatness traits mapped to the region around the IGF2 (insulin-like growth factor 2) gene but its non-imprinting inheritance excluded IGF2 as a candidate gene. A significant locus was also detected within a recombination cold spot that spans more than 30 Mb on SSCX, which hampered the identification of plausible candidate genes. Notably, no genome-wide significant locus was shared by the two experimental populations; different loci were observed that had both constant and time-specific effects on growth traits at different stages, which illustrates the complex genetic architecture of these traits.

CONCLUSIONS

We confirm several previously reported QTL and provide a list of novel loci for porcine growth and fatness traits in two experimental populations with Chinese Taihu and Western pigs as common founders. We showed that distinct loci exist for these traits in the two populations and identified HMGA1 and PLAG1 as strong candidate genes on SSC7 and SSC4, respectively.

Genome-wide identification, classification and functional analyses of the bHLH transcription factor family in the pig, Sus scrofa

The basic helix-loop-helix (bHLH) transcription factors are one of the largest families of gene regulatory proteins and play crucial roles in genetic, developmental and physiological processes in eukaryotes. Here, we conducted a survey of the Sus scrofa genome and identified 109 putative bHLH transcription factor members belonging to super-groups A, B, C, D, E, and F, respectively, while four members were orphan genes. We identified 6 most significantly enriched KEGG pathways and 116 most significant GO annotation categories. Further comprehensive surveys in human genome and other 12 medical databases identified 72 significantly enriched biological pathways with these 113 pig bHLH transcription factors. From the functional protein association network analysis 93 hub proteins were identified and 55 hub proteins created a tight network or a functional module within their protein families. Especially, there were 20 hub proteins found highly connected in the functional interaction network. The present study deepens our understanding and provided insights into the evolution and functional aspects of animal bHLH proteins and should serve as a solid foundation for further for analyses of specific bHLH transcription factors in the pig and other mammals.