A missense mutation in PPARD causes a major QTL effect on ear size in pigs

Two transects reveal remarkable variation in gene flow on opposite ends of a European toad hybrid zone

Speciation entails a reduction in gene flow between lineages. The rates at which genomic regions become isolated varies across space and time. Barrier markers are linked to putative genes involved in (processes of) reproductive isolation, and, when observed over two transects, indicate species-wide processes. In contrast, transect-specific putative barrier markers suggest local processes. We studied two widely separated transects along the 900 km hybrid zone between Bufo bufo and B. spinosus, in northern and southern France, for ~1200 RADseq markers. We used genomic and geographic cline analyses to identify barrier markers based on their restricted introgression, and found that some markers are transect-specific, while others are shared between transects. Twenty-six barrier markers were shared across both transects, of which some are clustered in the same chromosomal region, suggesting that their associated genes are involved in reduced gene flow across the entire hybrid zone. Transect-specific barrier markers were twice as numerous in the southern than in the northern transect, suggesting that the overall barrier effect is weaker in northern France. We hypothesize that this is consistent with a longer period of secondary contact in southern France. The smaller number of introgressed genes in the northern transect shows considerably more gene flow towards the southern (B. spinosus) than the northern species (B. bufo). We hypothesize that hybrid zone movement in northern France and hybrid zone stability in southern France explain this pattern. The Bufo hybrid zone provides an excellent opportunity to separate a general barrier effect from localized gene flow-reducing conditions.

Genome-wide association study reveals that the IBSP locus affects ear size in cattle

Ear size is a classical model for hot climate adaptation following the evolution, but the genetic basis of the traits associated with ear size remains to be elucidated. Here, we performed a genome-wide association study on 158 cattle to explain the genetic mechanism of ear size. One region on BTA6 between 36.79 and 38.80 Mb included 50 suggestive SNPs and 4 significant SNPs that were significantly associated with ear size. The most significant locus (P = 1.30 × 10-8) was a missense mutation (T250I) on the seventh exon of integrin-binding sialoprotein (IBSP), which had an allele substitution effect of 23.46 cm2 for ear size. Furthermore, this mutation will cause changes in the three-dimensional structure of the protein. To further identify genes underlying this typical feature, we performed a genome scan among nine cattle breeds with different ear sizes by using SweeD. Results suggested that IBSP was under positive selection among four breeds with relatively large ear sizes. The expression levels of IBSP in ear tissues of large- and small-ear cattle were significantly different. A haplotype diversity survey of this missense mutation in worldwide cattle breeds strongly implied that the origin of this missense mutation event was Bos taurus. These findings have important theoretical importance for the exploration of major genes associated with ear size and provide important molecular markers for the identification of cattle germplasm resources.

Identification of breed-specific genomic variants in Colombian Creole pig breeds by whole-genome sequencing

Dissecting genetic variation of local breeds is important for the success of conservation. In this research, we investigated the genomic variation of Colombian Creole (CR) pigs, with a focus on the breed-specific variants in the exonic region of 34 genes with reported effects on adaptive and economic traits. Seven individuals of each of the three CR breeds (CM, Casco de Mula; SP, San Pedreño; and ZU, Zungo) were whole-genome sequenced along with 7 Iberian (IB) pigs and 7 pigs of each of the four most used cosmopolitan (CP) breeds (Duroc, Landrace × Large White, and Pietrain). Molecular variability in CR (6,451,218 variants; from 3,919,242, in SP, to 4,648,069, in CM) was comparable to that in CP, but higher than in IB. For the investigated genes, SP pigs displayed less exonic variants (178) than ZU (254), CM (263), IB (200), and the individual CP genetic types (201 to 335). Sequence variation in these genes confirmed the resemblance of CR to IB and indicates that CR pigs, particularly ZU and CM, are not exempt from selective introgression of other breeds. A total of 50 exonic variants were identified as being potentially specific to CR, including a high-impact deletion in the intron between exons 15 and 16 of the leptin receptor gene, which was only found in CM and ZU. The identification of breed-specific variants in genes related to adaptive and economical traits can bolster the understanding of the role of gene-environment interactions on local adaptation and points the way for effective breeding and conservation of CR pigs.

Genetic structure and domestication footprints of the tusk, coat color, and ear morphology in East Chinese pigs

The domestication and artificial selection of wild boars have led to dramatic morphological and behavioral changes, especially in East Chinese (ECN) pigs. Here, we provide insights into the population structure and current genetic diversity of representative ECN pig breeds. We identify a 500-kb region containing six tooth development-relevant genes with almost completely different haplotypes between ECN pigs and Chinese wild boars or European domestic pigs. Notably, the c.195A>G missense mutation in exon 2 of AMBN may cause alterations in its protein structure associated with tusk degradation in ECN pigs. In addition, ESR1 may play an important role in the reproductive performance of ECN pigs. A major haplotype of the large lop ear-related MSRB3 gene and eight alleles in the deafness-related GRM7 gene may affect ear morphology and hearing in ECN pigs. Interestingly, we find that the two-end black (TEB) coat color in Jinhua pigs is most likely caused by EDNRB with genetic mechanisms different from other Chinese TEB pigs. This study identifies key loci that may be artificially selected in Chinese native pigs related to the tusk, coat color, and ear morphology, thus providing new insights into the genetic mechanisms of domesticated pigs.

Simultaneous testing of rule- and model-based approaches for runs of homozygosity detection opens up a window into genomic footprints of selection in pigs

BACKGROUND

Past selection events left footprints in the genome of domestic animals, which can be traced back by stretches of homozygous genotypes, designated as runs of homozygosity (ROHs). The analysis of common ROH regions within groups or populations displaying potential signatures of selection requires high-quality SNP data as well as carefully adjusted ROH-defining parameters. In this study, we used a simultaneous testing of rule- and model-based approaches to perform strategic ROH calling in genomic data from different pig populations to detect genomic regions under selection for specific phenotypes.

RESULTS

Our ROH analysis using a rule-based approach offered by PLINK, as well as a model-based approach run by RZooRoH demonstrated a high efficiency of both methods. It underlined the importance of providing a high-quality SNP set as input as well as adjusting parameters based on dataset and population for ROH calling. Particularly, ROHs ≤ 20 kb were called in a high frequency by both tools, but to some extent covered different gene sets in subsequent analysis of ROH regions common for investigated pig groups. Phenotype associated ROH analysis resulted in regions under potential selection characterizing heritage pig breeds, known to harbour a long-established breeding history. In particular, the selection focus on fitness-related traits was underlined by various ROHs harbouring disease resistance or tolerance-associated genes. Moreover, we identified potential selection signatures associated with ear morphology, which confirmed known candidate genes as well as uncovered a missense mutation in the ABCA6 gene potentially supporting ear cartilage formation.

CONCLUSIONS

The results of this study highlight the strengths and unique features of rule- and model-based approaches as well as demonstrate their potential for ROH analysis in animal populations. We provide a workflow for ROH detection, evaluating the major steps from filtering for high-quality SNP sets to intersecting ROH regions. Formula-based estimations defining ROHs for rule-based method show its limits, particularly for efficient detection of smaller ROHs. Moreover, we emphasize the role of ROH detection for the identification of potential footprints of selection in pigs, displaying their breed-specific characteristics or favourable phenotypes.

Genome-Wide Assessment of Runs of Homozygosity and Estimates of Genomic Inbreeding in a Chinese Composite Pig Breed

The primary purpose of the current study was to assess the genetic diversity, runs of homozygosity (ROH) and ROH islands in a Chinese composite pig and explore hotspot regions for traces of selection. First, we estimated the length, number, and frequency of ROH in 262 Xidu black pigs using the Porcine SNP50 BeadChip and compared the estimates of inbreeding coefficients, which were calculated based on ROHs (F_ROH) and homozygosity (F_HOM). Our result shows that a total of 7,248 ROH exceeding 1Mb were detected in 262 pigs. In addition, Sus scrofa chromosome (SSC) 8 and SSC10, respectively, has the highest and lowest chromosome coverage by ROH. These results suggest that inbreeding estimation based on total ROH may be a useful method, especially for crossbreed or composite populations. We also calculated an inbreeding coefficient of 0.077 from the total ROH. Eight ROH islands were found in this study. These ROH islands harbored genes associated with fat deposition, muscular development, reproduction, ear shape, and adaptation, such as TRAF7, IGFBP7, XPO1, SLC26A8, PPARD, and OR1F1. These findings may help to understand the effects of environmental and artificial selection on the genome structure of composite pigs. Our results provide a basis for subsequent genomic selection (GS), and provides a reference for the hybrid utilization of other pig breeds.

Whole genome variants across 57 pig breeds enable comprehensive identification of genetic signatures that underlie breed features

Background

A large number of pig breeds are distributed around the world, their features and characteristics vary among breeds, and they are valuable resources. Understanding the underlying genetic mechanisms that explain across-breed variation can help breeders develop improved pig breeds.

Results

In this study, we performed GWAS using a standard mixed linear model with three types of genome variants (SNP, InDel, and CNV) that were identified from public, whole-genome, sequencing data sets. We used 469 pigs of 57 breeds, and we identified and analyzed approximately 19 million SNPs, 1.8 million InDels, and 18,016 CNVs. We defined six biological phenotypes by the characteristics of breed features to identify the associated genome variants and candidate genes, which included coat color, ear shape, gradient zone, body weight, body length, and body height. A total of 37 candidate genes was identified, which included 27 that were reported previously (e.g., PLAG1 for body weight), but the other 10 were newly detected candidate genes (e.g., ADAMTS9 for coat color).

Conclusion

Our study indicated that using GWAS across a modest number of breeds with high density genome variants provided efficient mapping of complex traits.

Supplementary Information

Supplementary information accompanies this paper at 10.1186/s40104-020-00520-8.

Ear type in sheep is associated with the MSRB3 locus

Ear morphology is an important determinant of sheep breeds. It includes different variable traits such as ear size and erectness, suggesting a polygenic architecture. Here, we performed a comprehensive genome-wide analysis to identify regions under selection for ear morphology in 515 sheep from 17 breeds fixed for diverse ear phenotypes using 34k SNP genotyping data. GWASs for two ear type traits, size and erectness, revealed a single genome-wide significant association on ovine chromosome 3. The derived marker alleles were enriched in sheep with large and/or floppy ears. The GWAS signal harboured the MSRB3 gene encoding methionine sulphoxide reductase B3, which has already been found to be associated with different ear types in other species. We attempted whole-genome resequencing to identify causal variant(s) within a 1 Mb interval around MSRB3. This experiment excluded major copy number variants in the interval, but failed to identify a compelling candidate causal variant. Fine-mapping suggested that the causal variant for large floppy ears most likely resides in a 175 kb interval downstream of the MSRB3 coding region.

Comparative analysis of Longissimus dorsi tissue from two sheep groups identifies differentially expressed genes related to growth, development and meat quality

From a genetic perspective, the advantages of crossbreeding in sheep are unclear. In the present study, a comparative transcriptomic analysis was performed using Longissimus dorsi tissues from two sheep groups in order to identify differentially expressed genes (DEGs) related to growth, development and meat quality. Compared to Small Tail Han sheep, a total of 874 DEGs were identified in the crossbred sheep. Among these DEGs, 30, 116 and 32 DEGs were related to growth, development and meat quality, respectively. Seven DEGs highlighted by functional analysis as playing crucial roles in growth, development and meat quality were validated by the gene-act-network and co-expression-network. The expression levels of DEG mRNAs and proteins were further confirmed using RT-qPCR and western blot analyses. The results were consistent with the comparative transcriptome data. The data from this transcriptomic analysis will help to understand genetic heterosis and molecular-assisted breeding in sheep.

Genome-wide detection of copy number variants in European autochthonous and commercial pig breeds by whole-genome sequencing of DNA pools identified breed-characterising copy number states

In this study, we identified copy number variants (CNVs) in 19 European autochthonous pig breeds and in two commercial breeds (Italian Large White and Italian Duroc) that represent important genetic resources for this species. The genome of 725 pigs was sequenced using a breed-specific DNA pooling approach (30-35 animals per pool) obtaining an average depth per pool of 42×. This approach maximised CNV discovery as well as the related copy number states characterising, on average, the analysed breeds. By mining more than 17.5 billion reads, we identified a total of 9592 CNVs (~683 CNVs per breed) and 3710 CNV regions (CNVRs; 1.15% of the reference pig genome), with an average of 77 CNVRs per breed that were considered as private. A few CNVRs were analysed in more detail, together with other information derived from sequencing data. For example, the CNVR encompassing the KIT gene was associated with coat colour phenotypes in the analysed breeds, confirming the role of the multiple copies in determining breed-specific coat colours. The CNVR covering the MSRB3 gene was associated with ear size in most breeds. The CNVRs affecting the ELOVL6 and ZNF622 genes were private features observed in the Lithuanian Indigenous Wattle and in the Turopolje pig breeds respectively. Overall, the genome variability unravelled here can explain part of the genetic diversity among breeds and might contribute to explain their origin, history and adaptation to a variety of production systems.

Discovery of selection-driven genetic differences of Duroc, Landrace, and Yorkshire pig breeds by EigenGWAS and Fst analyses

Pigs are one of the earliest domesticated animals and multiple breeds have been developed to meet the various demands of consumers. EigenGWAS is a novel strategy to identify candidate genes that underlying population genetic differences and to infer candidate regions under selection as well. In this study, EigenGWAS and F_st analyses were performed using the public re-sequencing data of three typical commercial pig breeds, Duroc, Landrace and Yorkshire. The intersection of genome-wide significant SNPs detected by EigenGWAS and top-ranked 1% SNPs of F_st results were treated as signals under selection. Using the data of all three breeds, 3062 signals under selection were detected and the nearby genomic regions within 300 kb upstream and downstream covered 6.54% of whole genome. Pairs of breeds were analysed along with the pathway analysis. The gene function enrichment results indicated that many candidate genes located in the genomic regions of the signals under selection were associated with biological processes related to growth, metabolism, reproduction, sensory perception, etc. Among the candidate genes, the FSHB, AHR, PTHLH, KDR and FST genes were reported to be associated with reproductive performance; the KIT, KITLG, MITF, MC1R and EDNRB genes were previously identified to affect coat colour; the RETREG1, TXNIP, BMP5, PPARD and RBP4 genes were reported to be associated with lipid metabolism and growth traits. The identified genetic differences across the three commercial breeds will advance understanding of the artificial selection history of pigs and the signals under selection will suggest potential uses in pig genomic breeding programmes.

pCADD: SNV prioritisation in Sus scrofa

BACKGROUND

In animal breeding, identification of causative genetic variants is of major importance and high economical value. Usually, the number of candidate variants exceeds the number of variants that can be validated. One way of prioritizing probable candidates is by evaluating their potential to have a deleterious effect, e.g. by predicting their consequence. Due to experimental difficulties to evaluate variants that do not cause an amino-acid substitution, other prioritization methods are needed. For human genomes, the prediction of deleterious genomic variants has taken a step forward with the introduction of the combined annotation dependent depletion (CADD) method. In theory, this approach can be applied to any species. Here, we present pCADD (p for pig), a model to score single nucleotide variants (SNVs) in pig genomes.

RESULTS

To evaluate whether pCADD captures sites with biological meaning, we used transcripts from miRNAs and introns, sequences from genes that are specific for a particular tissue, and the different sites of codons, to test how well pCADD scores differentiate between functional and non-functional elements. Furthermore, we conducted an assessment of examples of non-coding and coding SNVs, which are causal for changes in phenotypes. Our results show that pCADD scores discriminate between functional and non-functional sequences and prioritize functional SNVs, and that pCADD is able to score the different positions in a codon relative to their redundancy. Taken together, these results indicate that based on pCADD scores, regions with biological relevance can be identified and distinguished according to their rate of adaptation.

CONCLUSIONS

We present the ability of pCADD to prioritize SNVs in the pig genome with respect to their putative deleteriousness, in accordance to the biological significance of the region in which they are located. We created scores for all possible SNVs, coding and non-coding, for all autosomes and the X chromosome of the pig reference sequence Sscrofa11.1, proposing a toolbox to prioritize variants and evaluate sequences to highlight new sites of interest to explain biological functions that are relevant to animal breeding.

GWAS on Imputed Whole-Genome Resequencing From Genotyping-by-Sequencing Data for Farrowing Interval of Different Parities in Pigs

The whole-genome sequencing (WGS) data can potentially discover all genetic variants. Studies have shown the power of WGS for genome-wide association study (GWAS) lies in the ability to identify quantitative trait loci and nucleotides (QTNs). However, the resequencing of thousands of target individuals is expensive. Genotype imputation is a powerful approach for WGS and to identify causal mutations. This study aimed to evaluate the imputation accuracy from genotyping-by-sequencing (GBS) to WGS in two pig breeds using a resequencing reference population and to detect single-nucleotide polymorphisms (SNPs) and candidate genes for farrowing interval (FI) of different parities using the data before and after imputation for GWAS. Six hundred target pigs, 300 Landrace and 300 Large White pigs, were genotyped by GBS, and 60 reference pigs, 20 Landrace and 40 Large White pigs, were sequenced by whole-genome resequencing. Imputation for pigs was conducted using Beagle software. The average imputation accuracy (allelic R ²) from GBS to WGS was 0.42 for Landrace pigs and 0.45 for Large White pigs. For Landrace pigs (Large White pigs), 4,514,934 (5,533,290) SNPs had an accuracy >0.3, resulting an average accuracy of 0.73 (0.72), and 2,093,778 (2,468,645) SNPs had an accuracy >0.8, resulting an average accuracy of 0.94 (0.93). Association studies with data before and after imputation were performed for FI of different parities in two populations. Before imputation, 18 and 128 significant SNPs were detected for FI in Landrace and Large White pigs, respectively. After imputation, 125 and 27 significant SNPs were identified for dataset with an accuracy >0.3 and 0.8 in Large White pigs, and 113 and 18 SNPs were found among imputed sequence variants. Among these significant SNPs, six top SNPs were detected in both GBS data and imputed WGS data, namely, SSC2: 136127645, SSC5: 103426443, SSC6: 27811226, SSC10: 3609429, SSC14: 15199253, and SSC15: 150297519. Overall, many candidate genes could be involved in FI of different parities in pigs. Although imputation from GBS to WGS data resulted in a low imputation accuracy, association analyses with imputed WGS data were optimized to detect QTNs for complex trait. The obtained results provide new insight into genotype imputation, genetic architecture, and candidate genes for FI of different parities in Landrace and Large White pigs.

Genome-wide association and evolutionary analyses reveal the formation of swine facial wrinkles in Chinese Erhualian pigs

Wrinkles are uneven concave-convex folds, ridges or creases in skin. Facial wrinkles appear in head, typically increasing along with aging. However in several Chinese indigenous pigs, such as Erhualian pigs, rich facial wrinkles have been generated during the growth stages as one of their breed characteristics. To investigate the genetic basis underlying the development of swine facial wrinkles, we estimated the folding extent of facial wrinkles in a herd of Erhualian pigs (n=332), and then conducted genome-wide association studies and multi-trait meta-analysis for facial wrinkles using 60K porcine chips. We found that facial wrinkles had high heritability estimates of ~0.7 in Erhualian pigs. Notably, only one genome-wide significant QTL was detected at 34.8 Mb on porcine chromosome 7. The most significant SNP rs80983858 located at the 3255-bp downstream of candidate gene GRM4, and the G allele was of benefit to increase facial wrinkles. Evolutionary and selection analyses suggested that the haplotypes containing G allele were under artificial selection, which was consistent with local animal sacrificial custom praying for longevity. Our findings made important clues for further deciphering the molecular mechanism of swine facial wrinkles formation, and shed light on the research of skin wrinkle development in human or other mammals.

Pig SOX9: Expression profiles of Sertoli cell (SCs) and a functional 18 bp indel affecting testis weight

Sex determining region Y-box 9 (SOX9), an important member of the SRY- type HMGbox (SOX) gene family, plays an important role in the regulation of mammalian reproduction, including sex differentiation during the embryonic development stage and spermatogenesis after birth. To explore the roles of polymorphism and expression of the SOX9 gene in the development of testes, we analyzed the indel of SOX9 in pigs and the corresponding expression level of the SOX9 gene in 7-day and 5-month-old porcine Sertoli cells. Results revealed that the DD haplotype of SOX9 gene as well as the ID genotype were significantly associated with larger testicular weight, while the II haplotype was closely related to the smaller testicular weight. More importantly, the SOX9 gene expression of ID genotyped group was significantly higher than that in II genotyped group. Our results first revealed that the indel polymorphism and expression of SOX9 were significantly associated with pig reproduction traits indicating the critical roles of SOX9 gene in testes development. The study provides a new clue for understanding the regulation of animal reproductive activities.

Molecular cloning of WIF1 and HMGA2 reveals ear-preferential expression while uncovering a missense mutation associated with porcine ear size in WIF1

Considerable diversity exists in porcine ear size, which is an important morphological feature of pig breeds. Previously, we localized four crucial candidate genes-high mobility group AT-hook 2 (HMGA2), LEM domain-containing 3 (LEMD3), methionine sulfoxide reductase B3 (MSRB3) and Wnt inhibitory factor 1 (WIF1)-on Sus Scrofa chromosome 5 affecting porcine ear size, then cloned LEMD3 and MSBR3. In this study, we performed rapid amplification of cDNA ends to obtain full-length cDNA sequences of 2338-bp WIF1 and 2998-bp HMGA2. Using quantitative real-time PCR, we revealed that WIF1 expression was highest in ear cartilage of 60-day-old pigs and that this is therefore a better candidate gene for ear size than HMGA2. We further screened coding sequence variants in both genes and identified only one missense mutation (WIF1:c.1167C>G) in a conserved epidermal growth factor-like domain from the mammalian WIF1 protein. The protein-altering mutation was significantly associated with ear size across the Large White × Minzhu hybrid and Beijing Black pig populations. When WIF1:c.1167C>G was included as fixed effect in the model to re-run a genome-wide association study in the Large White × Minzhu intercross population the P-value of the peak SNP on SSC5 from re-running the genome-wide association study dropped from 2.45E-12 to 7.33E-05. Taken together, the WIF1:c.1167C>G could be an important mutation associated with ear size. Our findings provide helpful information for further studies of the molecular mechanisms controlling porcine ear size.

Copy number variation in the MSRB3 gene enlarges porcine ear size through a mechanism involving miR-584-5p

Background

The size and type of ears are important conformation characteristics that distinguish pig breeds. A significant quantitative trait locus (QTL) for ear size has been identified on SSC5 (SSC for Sus scrofa chromosome) but the underlying causative gene and mutation remain unknown. Thus, our aim was to identify the gene responsible for enlarged ears in pig.

Results

First, we narrowed down the QTL region on SSC5 to a 137.85-kb interval that harbors only the methionine sulfoxide reductase B3 (MSRB3) gene. Then, we identified a 38.7-kb copy number variation (CNV) that affects the last two exons of MSRB3 and could be the candidate causative mutation for this QTL. This CNV showed complete concordance with genotype at the QTL of the founder animals in a white Duroc × Erhualian F₂ intercross and was found only in pigs from six Chinese indigenous breeds with large ears and from the Landrace breed with half-floppy ears. Moreover, it accounted for the significant association with ear size on SSC5 across the five pig populations tested. eQTL mapping revealed that this CNV was significantly associated with the expression of the microRNA (miRNA) miR-584-5p, which interacts with MSRB3, one of its target genes. In vivo and in vitro experiments confirmed that miR-584-5p inhibits the translation of MSRB3 mRNA. Taken together, these results led us to conclude that presence of the 38.7-kb CNV in the genome of some pig breeds affects ear size by altering the expression of miR-584-5p, which consequently hinders the expression of one of its target genes (e.g. MSRB3).

Conclusions

Our findings shed insight into the underlying mechanism of development of external ears in mammals and contribute to a better understanding of how the presence of CNV can regulate gene expression.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0442-6) contains supplementary material, which is available to authorized users.

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.

A genome-wide detection of selection signatures in conserved and commercial pig breeds maintained in Poland

BACKGROUND

Identification of selection signatures can provide a direct insight into the mechanism of artificial selection and allow further disclosure of the candidate genes related to the animals' phenotypic variation. Domestication and subsequent long-time selection have resulted in extensive phenotypic changes in domestic pigs, involving a number of traits, like behavior, body composition, disease resistance, reproduction and coat color. In this study, based on genotypes obtained from PorcineSNP60 Illumina assay we attempt to detect both diversifying and within-breed selection signatures in 530 pigs belonging to four breeds: Polish Landrace, Puławska, Złotnicka White and Złotnicka Spotted, of which the last three are a subject of conservative breeding and substantially represent the native populations.

RESULTS

A two largely complementary statistical methods were used for signatures detection, including: pairwise F_ST and relative extended haplotype homozygosity (REHH) test. Breed-specific diversifying selection signals included several genes involved in processes connected with fertility, growth and metabolism which are potentially responsible for different phenotypes of the studied breeds. The diversifying selection signals also comprised PPARD gene that was previously found to have a large effect on the shape of the external ear in pigs or two genes encoding neuropeptide Y receptors (Y2 and Y5) involved in fat deposition and stress response which are important features differentiating the studied breeds. REHH statistics allowed detecting several within-breed selection signatures overlapping with genes connected with a range of functions including, among others: metabolic pathways, immune system response or implantation and development of the embryo.

CONCLUSIONS

The study provides many potential candidate genes with implication for traits selected in the individual breeds and gives strong basis for further studies aiming at identification of sources of variation among the studied pig breeds.

Genome-wide association study reveals the locus responsible for microtia in Valle del Belice sheep breed

Microtia is a congenital deformity of the outer ear with phenotypes varying from a small auricle to total absence (anotia). The genetic basis is still poorly understood, and very few studies have been performed in sheep. Valle del Belice sheep is a breed showing microtia. The aim of this study was to identify the potential genomic regions involved in microtia in sheep. A total of 40 individuals, 20 with microtia and 20 normal, were genotyped with the Illumina OvineSNP50 BeadChip. The comparison among the results from a genome-wide association study, Fisher's exact test and F_ST analysis revealed a single strong association signal: rs419889303 on chromosome 1, located within intron 3 of the CLRN1 gene. Our study suggests for the first time that this novel candidate gene is responsible for microtia in sheep. Additional analysis based on the sequencing would help confirm our findings and allow for the proposal of a precise genetic basis for microtia in sheep.

Genomic and transcriptomic analyses reveal selection of genes for puberty in Bama Xiang pigs

The Bama Xiang pig (BMX) is a famous early-maturing Chinese indigenous breed with a two-end black coat. To uncover the genetic basis of the BMX phenotype, we conducted comparative genomic analyses between BMX and East Asian wild boars and Laiwu pigs, respectively. Genes under positive selection were enriched in pathways associated with gonadal hormone and melanin synthesis, consistent with the phenotypic changes observed during development in BMX pigs. We also performed differentially expressed gene analysis based on RNA-seq data from pituitary tissues of BMX and Large White pigs. The CTTNBP2NL, FRS2, KANK4, and KATNAL1 genes were under selection and exhibited expressional changes in the pituitary tissue, which may affect BMX pig puberty. Our study demonstrated the positive selection of early maturity in the development of BMX pigs and advances our knowledge on the role of regulatory elements in puberty evolution in pigs.

Genome-wide association study reveals novel genes for the ear size in sheep (Ovis aries)

Variations in ear size can be observed in livestock such as sheep; however, the genetic basis of variable ear size in sheep is still poorly understood. To investigate causative genes associated with ear size in sheep, a genome-wide association study was performed in 115 adult Duolang sheep with different-sized floppy ears using the Ovine Infinium HD BeadChip. We found 38 significant SNPs at the genome-wide or chromosome-wise 5% significance level after Bonferroni correction. The most significant association (P = 1.61 × 10^-6 ) was found at SNP rs402740419, located in the DCC gene, which plays a critical role in ear development. Also, we observed two additional significant SNPs, rs407891215 in PTPRD and rs407769095 in SOX5, both of which are functionally associated with ear developmental processes. Our results are useful for future sheep breeding and provide insights into the genetic basis of ear size development in sheep and other livestock.

Insertion/Deletion Within the KDM6A Gene Is Significantly Associated With Litter Size in Goat

A previous whole-genome association analysis identified lysine demethylase 6A (KDM6A), which encodes a type of histone demethylase, as a candidate gene associated to goat fecundity. KDM6A gene knockout mouse disrupts gametophyte development, suggesting that it has a critical role in reproduction. In this study, goat KDM6A mRNA expression profiles were determined, insertion/deletion (indel) variants in the gene identified, indel variants effect on KDM6A gene expression assessed, and their association with first-born litter size analyzed in 2326 healthy female Shaanbei white cashmere goats. KDM6A mRNA was expressed in all tissues tested (heart, liver, spleen, lung, kidney, muscle, brain, skin and testis); the expression levels in testes at different developmental stages [1-week-old (wk), 2, 3 wk, 1-month-old (mo), 1.5 and 2 mo] indicated a potential association with the mitosis-to-meiosis transition, implying that KDM6A may have an essential role in goat fertility. Meanwhile, two novel intronic indels of 16 bp and 5 bp were identified. Statistical analysis revealed that only the 16 bp indel was associated with first-born litter size (P < 0.01), and the average first-born litter size of individuals with an insertion/insertion genotype higher than that of those with the deletion/deletion genotype (P < 0.05). There was also a significant difference in genotype distributions of the 16 bp indel between mothers of single-lamb and multi-lamb litters in the studied goat population (P = 0.001). Consistently, the 16 bp indel also had a significant effect on KDM6A gene expression. Additionally, there was no significant linkage disequilibrium (LD) between these two indel loci, consistent with the association analysis results. Together, these findings suggest that the 16 bp indel in KDM6A may be useful for marker-assisted selection (MAS) of goats.

PigVar: a database of pig variations and positive selection signatures

Pigs are excellent large-animal models for medical research and a promising organ donor source for transplant patients. Next-generation sequencing technology has yielded a dramatic increase in the volume of genomic data for pigs. However, the limited amount of variation data provided by dbSNP, and non-congruent criteria used for calling variation, present considerable hindrances to the utility of this data. We used a uniform pipeline, based on GATK, to identify non-redundant, high-quality, whole-genome SNPs from 280 pigs and 6 outgroup species. A total of 64.6 million SNPs were identified in 280 pigs and 36.8 million in the outgroups. We then used LUMPY to identify a total of 7 236 813 structural variations (SVs) in 211 pigs. Positively selected loci were identified through five statistical tests of different evolutionary attributes of the SNPs. Combining the non-redundant variations and the evolutionary selective scores, we built the first pig-specific variation database, PigVar (http://www.ibiomedical.net/pigvar/), which is a web-based open-access resource. PigVar collects parameters of the variations including summary lists of the locations of the variations within protein-coding and long intergenic non-coding RNA (lincRNA) genes, whether the SNPs are synonymous or non-synonymous, their ancestral and derived states, geographic sampling locations, as well as breed information. The PigVar database will be kept operational and updated to facilitate medical research using the pig as model and agricultural research including pig breeding.

Database URL: http://www.ibiomedical.net/pigvar/

A combined genotype of three SNPs in the bovine gene is related to growth performance in Chinese cattle

Abstract. PPARD is involved in multiple biological processes, especially for those associated with energy metabolism. PPARD regulates lipid metabolism through up-regulate expression of genes associating with adipogenesis. This makes PPARD a significant candidate gene for production traits of livestock animals. Association studies between PPARD polymorphisms and production traits have been reported in pigs but are limited for other animals, including cattle. Here, we investigated the expression profile and polymorphism of bovine PPARD as well as their association with growth traits in Chinese cattle. Our results showed that the highest expression of PPARD was detected in kidney, following by adipose, which is consistent with its involvement in energy metabolism. Three SNPs of PPARD were detected and used to undergo selection pressure according the result of Hardy–Weinberg equilibrium analysis (P < 0.05). Moreover, all of these SNPs showed moderate diversity (0.25 < PIC < 0.5), indicating their relatively high selection potential. Association analysis suggested that individuals with the GAAGTT combined genotype of three SNPs detected showed optimal values in all of the growth traits analyzed. These results revealed that the GAAGTT combined genotype of three SNPs detected in the bovine PPARD gene was a significant potential genetic marker for marker-assisted selection in Chinese cattle. However, this should be further verified in larger populations before being applied to breeding.

PPARD is an Inhibitor of Cartilage Growth in External Ears

Peroxisome proliferator-activated receptor beta/delta (PPARD) is an important determinant of multiple biological processes. Our previous studies identified a missense mutation in the PPARD gene that significantly reduces its transcription activity, and consequently causes enlarged external ears in pigs. However, the mechanisms underlying the causality has remained largely unknown. Here, we show that PPARD retards the development of auricular cartilage by accelerating the apoptosis of cartilage stem/progenitor cells (CSPCs), the terminal differentiation of cartilage cells and the degradation of cartilage extracellular matrix in the auricle. At the transcription level, PPARD upregulates a set of genes that are associated with CSPCs apoptosis and chondrogenic differentiation, chondroblast differentiation and extracellular matrix degradation. ChIP-seq identified direct target genes of PPARD, including a well-documented gene for cartilage development: PPARG. We further show that compared to wild-type PPARD, the G32E mutant up-regulates the expression of PPARG and subsequently leads to the downregulation of critical genes that inhibit cartilage growth. These findings allow us to conclude that PPARD is an inhibitor of auricular cartilage growth in pigs. The causative mutation (G32E) in the PPARD gene attenuates the PPARD-mediated retardation of cartilage growth in the auricle, contributing to enlarged ears in pigs. The findings advance our understanding of the mechanisms underlying auricular development in mammals, and shed insight into the studies of innate pinna disorders and cartilage regeneration medicine in humans.

Imputation-Based Whole-Genome Sequence Association Study Rediscovered the Missing QTL for Lumbar Number in Sutai Pigs

Resequencing a number of individuals of various breeds as reference population and imputing the whole-genome sequences of individuals that were genotyped with medium-density chips to perform an association study is a very efficient strategy. Previously, we performed a genome-wide association study (GWAS) of lumbar number using 60K SNPs from the porcine Illumina chips in 418 Sutai pigs and did not detect any significant signals. Therefore, we imputed the whole-genome sequences of 418 Sutai individuals from 403 deeply resequenced reference individuals and performed association tests. We identified a quantitative trait locus (QTL) for lumbar number in SSC1 with a P value of 9.01E-18 that was close to the potential causative gene of NR6A1. The result of conditioning on the top SNP association test indicated that only one QTL was responsible for this trait in SSC1. The linkage disequilibrium (LD) drop test result for the condition of the reported potential causative mutation (c.575T > C missense mutation of NR6A1) indicated that this mutation was probably not the underlying mutation that affected lumbar number in our study. As the first trial of imputed whole-genome sequence GWAS in swine, this approach can be also powerful to investigate complex traits in pig like in human and cattle.

Transcription factor organic cation transporter 1 (OCT-1) affects the expression of porcine Klotho (KL) gene

Klotho (KL), originally discovered as an aging suppressor, is a membrane protein that shares sequence similarity with the β-glucosidase enzymes. Recent reports showed Klotho might play a role in adipocyte maturation and systemic glucose metabolism. However, little is known about the transcription factors involved in regulating the expression of porcine KL gene. Deletion fragment analysis identified KL-D2 (−418 bp to −3 bp) as the porcine KL core promoter. MARC0022311SNP (A or G) in KL intron 1 was detected in Landrace × DIV pigs using the Porcine SNP60 BeadChip. The pGL-D2-A and pGL-D2-G were constructed with KL-D2 and the intron fragment of different alleles and relative luciferase activity of pGL3-D2-G was significantly higher than that of pGL3-D2-A in the PK cells and ST cells. This was possibly the result of a change in KL binding ability with transcription factor organic cation transporter 1 (OCT-1), which was confirmed using electrophoretic mobility shift assays (EMSA) and chromatin immune-precipitation (ChIP). Moreover, OCT-1 regulated endogenous KL expression by RNA interference experiments. Our study indicates SNP MARC0022311 affects porcine KL expression by regulating its promoter activity via OCT-1.

A genome-wide association analysis for susceptibility of pigs to enterotoxigenic Escherichia coli F41

Enterotoxigenic Escherichia coli (ETEC) is a type of pathogenic bacteria that cause diarrhea in piglets through colonizing pig small intestine epithelial cells by their surface fimbriae. Different fimbriae type of ETEC including F4, F18, K99 and F41 have been isolated from diarrheal pigs. In this study, we performed a genome-wide association study to map the loci associated with the susceptibility of pigs to ETEC F41 using 39454 single nucleotide polymorphisms (SNPs) in 667 F2 pigs from a White Duroc×Erhualian F2 cross. The most significant SNP (ALGA0022658, P=5.59×10-13) located at 6.95 Mb on chromosome 4. ALGA0022658 was in high linkage disequilibrium (r 2>0.5) with surrounding SNPs that span a 1.21 Mb interval. Within this 1.21 Mb region, we investigated ZFAT as a positional candidate gene. We re-sequenced cDNA of ZFAT in four pigs with different susceptibility phenotypes, and identified seven coding variants. We genotyped these seven variants in 287 unrelated pigs from 15 diverse breeds that were measured with ETEC F41 susceptibility phenotype. Five variants showed nominal significant association (P<0.05) with ETEC F41 susceptibility phenotype in International commercial pigs. This study provided refined region associated with susceptibility of pigs to ETEC F41 than that reported previously. Further works are needed to uncover the underlying causal mutation(s).

SNP annotation-based whole genomic prediction and selection: an application to feed efficiency and its component traits in pigs

The study investigated genetic architecture and predictive ability using genomic annotation of residual feed intake (RFI) and its component traits (daily feed intake [DFI], ADG, and back fat [BF]). A total of 1,272 Duroc pigs had both genotypic and phenotypic records, and the records were split into a training (968 pigs) and a validation dataset (304 pigs) by assigning records as before and after January 1, 2012, respectively. SNP were annotated by 14 different classes using Ensembl variant effect prediction. Predictive accuracy and prediction bias were calculated using Bayesian Power LASSO, Bayesian A, B, and Cπ, and genomic BLUP (GBLUP) methods. Predictive accuracy ranged from 0.508 to 0.531, 0.506 to 0.532, 0.276 to 0.357, and 0.308 to 0.362 for DFI, RFI, ADG, and BF, respectively. BayesCπ100.1 increased accuracy slightly compared to the GBLUP model and other methods. The contribution per SNP to total genomic variance was similar among annotated classes across different traits. Predictive performance of SNP classes did not significantly differ from randomized SNP groups. Genomic prediction has accuracy comparable to observed phenotype, and use of genomic prediction can be cost effective by replacing feed intake measurement. Genomic annotation had less impact on predictive accuracy traits considered here but may be different for other traits. It is the first study to provide useful insights into biological classes of SNP driving the whole genomic prediction for complex traits in pigs.

Genetic architecture of fatty acid composition in the longissimus dorsi muscle revealed by genome-wide association studies on diverse pig populations

Background

Fatty acid composition in muscle is an important factor that affects the nutritive value and taste of pork. To investigate the genetic architecture of fatty acid composition of pork, we measured fatty acid contents in longissimus dorsi muscle of 1244 pigs from three divergent populations and conducted genome-wide association studies (GWAS) for fatty acid contents.

Results

We detected 26 genome-wide significant quantitative trait loci (QTL) on eight chromosomes (SSC for Sus scrofa) for eight fatty acids. These loci not only replicated previously reported QTL for C18:0 on SSC14 and C20:0 on SSC16, but also included several novel QTL such as those for C20:1 on SSC7, C14:0 on SSC9, and C14:0, C16:0 and C16:1 on SSC12. Furthermore, we performed a meta-analysis of GWAS on five populations, including the three populations that were investigated in this study and two additional populations that we had previously examined. This enhanced the strength of the associations detected between fatty acid composition and several marker loci, especially for those for C18:0 on SSC14 and C20:0 on SSC16. The genes ELOVL5, ELOVL6, ELOVL7, FASN, SCD and THRSP, which have functions that are directly relevant to fatty acid metabolism, are proximal to the top associated markers at six significant QTL.

Conclusions

The findings improve our understanding of the genetic architecture of fatty acid composition in pork and contribute to further fine-map and characterize genes that influence fatty acid composition.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0184-2) contains supplementary material, which is available to authorized users.

Porcine methionine sulfoxide reductase B3: molecular cloning, tissue-specific expression profiles, and polymorphisms associated with ear size in Sus scrofa

Background

In Sus scrofa, methionine sulfoxide reductase B3 (MSRB3) is a crucial candidate gene for ear size, and an important conformational trait of pig breeds. However, challenges in MSRB3 cDNA amplification have prevented further identification of MSRB3 allelic variants influencing pig ear size.

Results

We cloned a full-length cDNA sequence of porcine MSRB3 by rapid-amplification of cDNA ends. The 3,765-bp gene contained a 5’-untranslated region (UTR) (190 bp), a coding region (552 bp), and a 3’-UTR (3,016 bp) and shared 84 %, 84 %, 87 %, 86 %, and 70 % sequence identities with human, orangutan, mouse, chicken, and zebrafish, respectively. The gene encoded a 183-amino acid protein, which shared 88 %, 91 %, 89 %, 86 %, and 67 % identities with human, orangutan, mouse, chicken, and zebrafish, respectively. Tissue expression analysis using qRT-PCR revealed that MSRB3 was expressed in the heart, liver, lung, kidney, spleen, ear, muscle, fat, lymph, skeletal, and hypothalamic tissues. Three single nucleotide polymorphisms (SNPs) were identified in MSRB3: c.-735C > T in the 5’ flanking region, c.2571 T > C in the 3’-UTR, and a synonymous mutation of c.484 T > C in the coding region. The SNPs c.-735C > T and c.2571 T > C were significantly associated with ear size in a Large White × Minzhu F2 population other than in Beijing Black pigs. Subsequently, at SNP c.-735C > T, the mRNA of MSRB3 was significantly higher expressed in ears of individuals with the TT genotype (Minzhu) than those with CC (Large White).

Conclusions

The porcine MSRB3 owned a 3,765-bp full-length cDNA sequence and was detected to express in ear tissue. Two SNPs of this gene were shown to be significantly associated with ear size in a Large White × Minzhu intercross population instead of Beijing Black pig population. What’s more, the individuals with higher mRNA expression of MSRB3 have larger ear sizes. These results provide useful information for further functional analyses of MSRB3 influencing ear size in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s40104-015-0060-x) contains supplementary material, which is available to authorized users.

The Functional SNPs in the 5' Regulatory Region of the Porcine PPARD Gene Have Significant Association with Fat Deposition Traits

Peroxisome proliferator-activated receptor delta (PPARD) is a key regulator of lipid metabolism, insulin sensitivity, cell proliferation and differentiation. In this study, we identified two Single Nucleotide Polymorphisms (SNPs, g.1015 A>G and g.1018 T>C) constituting four haplotypes (GT, GC, AC and AT) in the 5’ regulatory region of porcine PPARD gene. Functional analysis of the four haplotypes showed that the transcriptional activity of the PPARD promoter fragment carrying haplotype AC was significantly lower than that of the other haplotypes in 3T3-L1, C2C12 and PK-15 cells, and haplotype AC had the lowest binding capacities to the nuclear extracts. Transcription factor 7-like 2 (TCF7L2) enhanced the transcription activities of promoter fragments of PPARD gene carrying haplotypes GT, GC and AT in C2C12 and 3T3-L1 cells, and increased the protein expression of PPARD gene in C2C12 myoblasts. TCF7L2 differentially bound to the four haplotypes, and the binding capacity of TCF7L2 to haplotype AC was the lowest. There were significant associations between -655A/G and fat deposition traits in three pig populations including the Large White × Meishan F₂ pigs, France and American Large White pigs. Pigs with genotype GG had significantly higher expression of PPARD at both mRNA and protein level than those with genotype AG. These results strongly suggested that the SNPs in 5’ regulatory region of PPARD genes had significant impact on pig fat deposition traits.

The Genetic Architecture of Domestication in Animals

Domestication has been essential to the progress of human civilization, and the process itself has fascinated biologists for hundreds of years. Domestication has led to a series of remarkable changes in a variety of plants and animals, in what is termed the “domestication phenotype.” In domesticated animals, this general phenotype typically consists of similar changes in tameness, behavior, size/morphology, color, brain composition, and adrenal gland size. This domestication phenotype is seen in a range of different animals. However, the genetic basis of these associated changes is still puzzling. The genes for these different traits tend to be grouped together in clusters in the genome, though it is still not clear whether these clusters represent pleiotropic effects, or are in fact linked clusters. This review focuses on what is currently known about the genetic architecture of domesticated animal species, if genes of large effect (often referred to as major genes) are prevalent in driving the domestication phenotype, and whether pleiotropy can explain the loci underpinning these diverse traits being colocated.

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.

A genome scan for selection signatures in pigs

Identifying signatures of selection can provide a straightforward insight into the mechanism of artificial selection and further uncover the causal genes related to the phenotypic variation. Based on Illumina Porcine60KSNP chip data, four complementary methods, Long-Range Haplotype (LRH), Tajima’s D, Cross Population Extend Haplotype Homozygosity Test (XPEHH) and F_ST, were implemented in this study to detect the selection signatures in the whole genome of one typical Chinese indigenous breed, Rongchang, one Chinese cultivated breed, Songliao, and two western breeds, Landrace and Yorkshire. False Discovery Rate (FDR) was implemented to control the false positive rates. In our study, a total of 159, 127, 179 and 159 candidate selection regions with average length of 0.80 Mb, 0.73 Mb, 0.78 Mb and 0.73 Mb were identified in Landrace, Rongchang, Songliao and Yorkshire, respectively, that span approximately 128.00 Mb, 92.38 Mb, 130.30 Mb and 115.40 Mb and account for approximately 3.74–5.33% of genome across all autosomes. The selection regions of 11.52 Mb shared by Landrace and Yorkshire were the longest when chosen pairs from the pool of the four breeds were examined. The overlaps between Yorkshire and Songliao, approximately 9.20 Mb, were greater than those of Yorkshire and Rongchang. Meanwhile, the overlaps between Landrace and Songliao were greater than those of Landrace and Rongchang but less than those of Songliao and Ronchang. Bioinformatics analysis showed that the genes/QTLs relevant to fertility, coat color, and ear morphology were found in candidate selection regions. Some genes, such as LEMD3, MC1R, KIT, TRHR etc. that were reported under selection, were confirmed in our study, and this analysis also demonstrated the diversity of breeds.

Genome-wide association analysis reveals genetic loci and candidate genes for meat quality traits in Chinese Laiwu pigs

Meat quality traits have economically significant impacts on the pig industry, and can be improved using molecular approaches in pig breeding. Since 1994 when the first genome-wide scan for quantitative trait loci (QTLs) in pig was reported, over the past two decades, numerous QTLs have been identified for meat quality traits by family based linkage analyses. However, little is known about the genetic variants for meat quality traits in Chinese purebred or outbred populations. To unveil it, we performed a genome-wide association study for 10 meat quality traits in Chinese purebred Laiwu pigs. In total, 75 significant SNPs (P < 1.01 × 10(-6)) and 33 suggestive SNPs (P < 2.03 × 10(-5)) were identified. On SSC12, a region between 56.22 and 61.49 Mb harbored a cluster of SNPs that were associated with meat color parameters (L*, lightness; a*, redness; b*, yellowness) and moisture content of longissimus muscle (LM) and semimembranosus muscle at the genome-wide significance level. A region on SSC4 also has pleiotropic effects on moisture content and drip loss of LM. In addition, this study revealed at least five novel QTLs and several candidate genes including 4-linked MYH genes (MYH1, MYH2, MYH3, and MYH13), MAL2, LPAR1, and PRKAG3 at four significant loci. Except for the SSC12 QTL, other QTLs are likely tissue-specific. These results provide new insights into the genetic basis of meat quality traits in Chinese Laiwu pigs and some significant SNPs reported here could be incorporated into the selection programs involving this breed.

Mining the pig genome to investigate the domestication process

Pig domestication began around 9000 YBP in the Fertile Crescent and Far East, involving marked morphological and genetic changes that occurred in a relatively short window of time. Identifying the alleles that drove the behavioural and physiological transformation of wild boars into pigs through artificial selection constitutes a formidable challenge that can only be faced from an interdisciplinary perspective. Indeed, although basic facts regarding the demography of pig domestication and dispersal have been uncovered, the biological substrate of these processes remains enigmatic. Considerable hope has been placed on new approaches, based on next-generation sequencing, which allow whole-genome variation to be analyzed at the population level. In this review, we provide an outline of the current knowledge on pig domestication by considering both archaeological and genetic data. Moreover, we discuss several potential scenarios of genome evolution under the complex mixture of demography and selection forces at play during domestication. Finally, we highlight several technical and methodological approaches that may represent significant advances in resolving the conundrum of livestock domestication.

A splice mutation in the PHKG1 gene causes high glycogen content and low meat quality in pig skeletal muscle

Glycolytic potential (GP) in skeletal muscle is economically important in the pig industry because of its effect on pork processing yield. We have previously mapped a major quantitative trait loci (QTL) for GP on chromosome 3 in a White Duroc × Erhualian F₂ intercross. We herein performed a systems genetic analysis to identify the causal variant underlying the phenotype QTL (pQTL). We first conducted genome-wide association analyses in the F₂ intercross and an F₁₉ Sutai pig population. The QTL was then refined to an 180-kb interval based on the 2-LOD drop method. We then performed expression QTL (eQTL) mapping using muscle transcriptome data from 497 F₂ animals. Within the QTL interval, only one gene (PHKG1) has a cis-eQTL that was colocolizated with pQTL peaked at the same SNP. The PHKG1 gene encodes a catalytic subunit of the phosphorylase kinase (PhK), which functions in the cascade activation of glycogen breakdown. Deep sequencing of PHKG1 revealed a point mutation (C>A) in a splice acceptor site of intron 9, resulting in a 32-bp deletion in the open reading frame and generating a premature stop codon. The aberrant transcript induces nonsense-mediated decay, leading to lower protein level and weaker enzymatic activity in affected animals. The mutation causes an increase of 43% in GP and a decrease of>20% in water-holding capacity of pork. These effects were consistent across the F₂ and Sutai populations, as well as Duroc × (Landrace × Yorkshire) hybrid pigs. The unfavorable allele exists predominantly in Duroc-derived pigs. The findings provide new insights into understanding risk factors affecting glucose metabolism, and would greatly contribute to the genetic improvement of meat quality in Duroc related pigs.

Glycogen storage diseases (GSD) are a group of inherited disorders characterized by storage of excess glycogen, which are mainly caused by the abnormality of a particular enzyme essential for releasing glucose from glycogen. GSD-like conditions have been described in a wide variety of species. Pigs are a valuable model for the study of human GSD. Moreover, pigs affected by GSD usually produce inferior pork with a lower ultimate pH (so-called “acid meat”) and less processing yield due to post-mortem degradation of the excess glycogen. So far, only one causal variant, PRKAG3 R225Q, has been identified for GSD in pigs. Here we reported a loss-of-function mutation in the PHKG1 gene that causes the deficiency of the glycogen breakdown, consequently leading to GSD and acid meat in Duroc-sired pigs. Eliminating the undesirable mutation from the breeding stock by a diagnostic DNA test will greatly reduce the incidence of GSD and significantly improve pork quality and productivity in the pig.