Genome-wide analysis reveals artificial selection on coat colour and reproductive traits in Chinese domestic pigs

DeepOCR: A multi-species deep-learning framework for accurate identification of open chromatin regions in livestock

A wealth of experimental evidence has suggested that open chromatin regions (OCRs) are involved in many critical biological activities, such as DNA replication, enhancer activity, and gene transcription. Accurately identifying OCRs in livestock species can provide critical insights into the distribution and characteristics of OCRs for disease treatment in livestock, thereby improving animal welfare. However, most current machine-learning methods for OCR prediction were originally designed for a limited number of model organisms, such as humans and some model organisms, and thus their performance on non-model organisms, specifically livestock, is often unsatisfactory. To bridge this gap, we propose DeepOCR, a lightweight depth-separable residual network model for predicting OCRs in livestock, including chicken, cattle, and sheep. DeepOCR integrates a single convolution layer and two improved residue structure blocks to extract and learn important features from the input DNA sequences. A fully connected layer was also employed to further process the extracted features and improve the robustness of the entire network. Our benchmarking experiments demonstrated superior prediction performance of DeepOCR compared to state-of-the-art approaches on testing datasets of the three species. The source code of DeepOCR is freely available for academic purposes at https://github.com/jasonzhao371/DeepOCR/. We anticipate DeepOCR servers as a practical and reliable computational tool for OCR-related studies in livestock species.

Whole-genome resequencing of Hu sheep identifies candidate genes associated with agronomic traits

The phenotypic diversity resulting from artificial or natural selection of sheep has made a significant contribution to human civilization. Hu sheep are a local sheep breed unique to China with high reproductive rates and rapid growth. Genomic selection signatures have been widely used to investigate the genetic mechanisms underlying phenotypic variation in livestock. Here, we conduct whole-genome sequencing of 207 Hu sheep and compare them with the wild ancestors of domestic sheep (Asiatic mouflon) to investigate the genetic characteristics and selection signatures of Hu sheep. Based on six signatures of selection approaches, we detect genomic regions containing genes related to reproduction (BMPR1B, BMP2, PGFS, CYP19, CAMK4, GGT5, and GNAQ), vision (ALDH1A2, SAG, and PDE6B), nervous system (NAV1), and immune response (GPR35, SH2B2, PIK3R3, and HRAS). Association analysis with a population of 1299 Hu sheep reveal those missense mutations in the GPR35 (GPR35 g.952651 A>G; GPR35 g.952496 C>T) and NAV1 (NAV1 g.84216190 C>T; NAV1 g.84227412 G>A) genes are significantly associated (P < 0.05) with immune and growth traits in Hu sheep, respectively. This research offers unique insights into the selection characteristics of Hu sheep and facilitates further genetic improvement and molecular investigations.

Copy number variants selected during pig domestication inferred from whole genome resequencing

Over extended periods of natural and artificial selection, China has developed numerous exceptional pig breeds. Deciphering the germplasm characteristics of these breeds is crucial for their preservation and utilization. While many studies have employed single nucleotide polymorphism (SNP) analysis to investigate the local pig germplasm characteristics, copy number variation (CNV), another significant type of genetic variation, has been less explored in understanding pig resources. In this study, we examined the CNVs of 18 Wanbei pigs (WBP) using whole genome resequencing data with an average depth of 12.61. We identified a total of 8,783 CNVs (~30.07 Mb, 1.20% of the pig genome) in WBP, including 8,427 deletions and 356 duplications. Utilizing fixation index (Fst), we determined that 164 CNVs were within the top 1% of the Fst value and defined as under selection. Functional enrichment analyses of the genes associated with these selected CNVs revealed genes linked to reproduction (SPATA6, CFAP43, CFTR, BPTF), growth and development (NR6A1, SMYD3, VIPR2), and immunity (PARD3, FYB2). This study enhances our understanding of the genomic characteristics of the Wanbei pig and offers a theoretical foundation for the future breeding of this breed.

Unveiling the Genetic Secrets of Chinese Indigenous Pigs from Guizhou Province: Diversity, Evolution and Candidate Genes Affecting Pig Coat Color

The local pig breeds in Guizhou possess exceptional meat quality, robust adaptability, and resilience to harsh feeding conditions, making them ideal for producing high-quality pork. With over 10 local pig breeds in the region, we focused on 7 specific breeds: Baixi pigs (BX), Congjiang Xiang pigs (CJX), Guanling pigs (GL), Jianhe White Xiang pigs (JHBX), Jiangkou Luobo pigs (JKLB), Kele pigs (KL), and Qiandong Hua pigs (QDH). Unfortunately, these breeds face threats such as introduced species and inbreeding, resulting in a decline in population size and numbers. To better protect and utilize these breeds, we employed genome-wide single-nucleotide polymorphism (SNP) markers to investigate the population structure, genetic diversity, and selection characteristics of 283 pigs across these seven breeds. Our findings revealed distinct ancestral sources between Chinese and Western pig breeds, as demonstrated by principal component analysis, adjacent tree analysis, and ADMIXTURE analysis. Notably, JHBX exhibited a distant genetic relationship from the other six local pig breeds in Guizhou province, showcasing unique genetic characteristics. While the genetic diversity of the six Chinese native pig populations, excluding JHBX, was generally moderate in Guizhou province, the JHBX population displayed low genetic diversity. Therefore, it is imperative to intensify selection efforts to prevent inbreeding decline in JHBX while further enhancing the protection measures for the other six pig populations. Additionally, we identified candidate genes influencing the size disparity among pigs in Guizhou province through signal selection. Our study outcomes serve as a reference for developing effective conservation and utilization plans for pig breeds in Guizhou province and deepen our understanding of the genetic mechanisms underlying pig body size.

Detection of Runs of Homozygosity and Identification of Candidate Genes in the Whole Genome of Tunchang Pigs

Tunchang pigs are an indigenous pig population in China known for their high tolerance to roughage, delicious meat, and fecundity. However, the number of Tunchang pigs has been declining due to the influence of commercial breeds and African swine fever, which could potentially lead to inbreeding. To assess the inbreeding level and the genetic basis of important traits in Tunchang pigs, our research investigated the patterns in "runs of homozygosity" (ROHs) using whole genome resequencing data from 32 Tunchang pigs. The study aimed to determine the length, number, coverage, and distribution model of ROHs in Tunchang pigs, as well as genomic regions with high ROH frequencies. The results of the study revealed that a total of 20,499,374 single-nucleotide polymorphisms (SNPs) and 1953 ROH fragments were recognized in 32 individuals. The ROH fragments in Tunchang pigs were predominantly short, ranging from 0.5 to 1 megabases (Mb) in length. Furthermore, the coverage of ROHs varied across chromosomes, with chromosome 3 having the highest coverage and chromosome 11 having the lowest coverage. The genetic diversity of Tunchang pigs was found to be relatively high based on the values of HE (expected heterozygosity), HO (observed heterozygosity), pi (nucleotide diversity), Ne (effective population size), and MAF (minor allele frequency). The average inbreeding coefficients of Tunchang pigs, as determined by three different methods (FHOM, FGRM, and FROH), were 0.019, 0.0138, and 0.0304, respectively. These values indicate that the level of inbreeding in Tunchang pigs is currently low. Additionally, the study identified a total of 13 ROH islands on all chromosomes, which in total contained 38,913 SNPs and 120 genes. These ROH islands included genes associated with economically important traits, including meat quality (GYS1, PHLPP1, SLC27A5, and CRTC1), growth and development (ANKS1A, TAF11, SPDEF, LHB, and PACSIN1), and environmental adaptation (SLC26A7). The findings of this research offer valuable perspectives on the present status of Tunchang pig resources and offer a reference for breeding conservation plans and the efficient utilization of Tunchang pigs in the future. By understanding the inbreeding level and genetic basis of important traits in Tunchang pigs, conservation efforts can be targeted towards maintaining genetic diversity and promoting the sustainable development of this indigenous pig population.

Genome-wide detection of selection signatures in Jianli pigs reveals novel cis-regulatory haplotype in EDNRB associated with two-end black coat color

BACKGROUND

Jianli pig, a renowned indigenous breed in China, has the characteristics of a two-end black (TEB) coat color, excellent meat quality, strong adaptability and increased prolificacy. However, there is limited information available regarding the genetic diversity, population structure and genomic regions under selection of Jianli pig. On the other hand, the genetic mechanism of TEB coat color has remained largely unknown.

RESULTS

In this study, the whole genome resequencing of 30 Jianli pigs within a context of 153 individuals representing 13 diverse breeds was performed. The population structure analysis revealed that Jianli pigs have close genetic relationships with the Tongcheng pig breed, their geographical neighbors. Three methods (observed heterozygosity, expected heterozygosity, and runs of homozygosity) implied a relatively high level of genetic diversity and, a low inbreeding coefficient in Jianli compared with other pigs. We used Fst and XP-EHH to detect the selection signatures in Jianli pigs compared with Asian wild boar. A total of 451 candidate genes influencing meat quality (CREBBP, ADCY9, EEPD1 and HDAC9), reproduction (ESR1 and FANCA), and coat color (EDNRB, MITF and MC1R), were detected by gene annotation analysis. Finally, to fine-map the genomic region for the two-end black (TEB) coat color phenotype in Jianli pigs, we performed three signature selection methods between the TEB coat color and no-TEB coat color pig breeds. The current study, further confirmed that the EDNRB gene is a candidate gene for TEB color phenotype found in Chinese pigs, including Jinhua pigs, and the haplotype harboring 25 SNPs in the EDNRB gene may promote the formation of TEB coat color. Further ATAC-seq and luciferase reporter assays of these regions suggest that the 25-SNPs region was a strong candidate causative mutation that regulates the TEB coat color phenotype by altering enhancer function.

CONCLUSION

Our results advanced the understanding of the genetic mechanism behind artificial selection, and provided further resources for the protection and breeding improvement of Jianli pigs.

Evaluation of the Genetic Diversity and Population Structure of Four Native Pig Populations in Gansu Province

Indigenous pig populations, including Bamei pigs (BM), Hezuo pigs (HZ), Huixian Qingni Black pigs (HX), and Minxian Black pigs (MX) in Gansu Province, live in a particular climate and a relatively closed geographical environment. These local pig breeds are characterized by excellent characteristics (e.g., cold tolerance, robust disease resistance, and superior meat quality). In the past few years, pig populations in Gansu Province have decreased significantly because of their poor lean meat percentage, high fat content, and slow growth rate. Maintaining the diversity of these four breeds can act as a source of new alleles to be incorporated into commercial breeds which are more susceptible to disease and less adaptable to changing conditions because of inbreeding. Genomic data analysis is adequate for determining the genetic diversity and livestock breeding population structure, even in local pig populations. However, the genetic diversity and population structure of the four native pig populations in Gansu Province are still unknown. Thus, we used "Zhongxin-I" porcine chip for the SNP detection of 102 individuals living on four pig conservation farms. A total of 57,466 SNPs were identified among the four pig breeds. The linkage disequilibrium (LD) plot showed that MX had the highest level of LD, followed by BM, HZ, and HX. The observed heterozygosity (Ho) in all four populations was higher than the expected heterozygosity (He). A principal component analysis (PCA) demonstrated that the four local pig populations were isolated. The identity displayed by the state matrix and G matrix heat map results indicated that small numbers of individuals among the four pig breeds had a high genetic distance and weak genetic relationships. The results of the population genetic structure of BM, HZ, HX, and MX pigs showed a slight genetic diversity loss. Our findings enabled us to better understand the genome characteristics of these four indigenous pig populations, which will provide novel insights for the future germplasm conservation and utilization of these indigenous pig populations.

Genome-wide population structure, homozygosity, and heterozygosity patterns of Nero Siciliano pig in the framework of Italian and cosmopolitan breeds

Analysis of genomic data is becoming more and more common for the effective management of livestock breeding programmes, even in the case of local populations. In this work, the genome-wide data of Nero Siciliano pig breed were compared to that of wild boar, Italian local and cosmopolitan breeds to investigate its genetic structure, and runs of homozygosity (ROH) and heterozygosity patterns. The Nero Siciliano has been reported to have the highest rate of genetic diversity among the Italian breeds, and a genetic variability comparable to that of the cosmopolitan breeds. Analyses of genomic structure and relationships underlined its proximity to wild boar, and an internal substructure probably linked to different family lines. The breed showed a low value of inbreeding estimated from ROH, and the highest diversity index among the Italian breeds, even if lower than that of the cosmopolitans. Four ROH islands in three chromosomes (SSC8, SSC11, and SSC14) and one heterozygosity-rich region (SSC1) were identified in Nero Siciliano, highlighting genomic regions related to productive QTL. Across breeds, SSC8 and SSC14 were the chromosomes with most ROH islands, with Mora Romagnola and wild boar showing the highest level of autozygosity. Chromosomes SSC2, SSC6, SSC8 and SSC13 showed the majority of runs of heterozygosity regions, mainly found in the cosmopolitan pig breeds, which reported several genes associated with health-related QTL. The outlined results can help to better identify the genomic profile of this local breed in order to plan matings, maintain adequate internal diversity and exploit the production system.

Genetic and Phenotypic Characteristics of Belted Pig Breeds: A Review

Belted pig breeds have unique, distinguishing phenotypic characteristics. This review summarises the current knowledge on pig breeds displaying a belted coat pattern. Belts of different widths and positions around the animal's trunk characterise specific pig breeds from all around the world. All the breeds included in the present paper have been searched through the FAO domestic animal diversity information system (DAD-IS), Every country was checked to identify all breeds described as having black or red piebald coat pattern variations. Advances in genomic technologies have made it possible to identify the specific genes and genetic markers associated with the belted phenotype and explore the genetic relationships between different local breeds. Thus, the origin, history, and production traits of these breeds, together with all the genomic information related to the mechanism of skin pigmentation, are discussed. By increasing our understanding of these breeds, we can appreciate the richness of our biological and cultural heritage and work to preserve the biodiversity of the world's animals.

Evaluation of the Genetic Diversity, Population Structure and Selection Signatures of Three Native Chinese Pig Populations

Indigenous pig populations in Hainan Province live in tropical climate conditions and a relatively closed geographical environment, which has contributed to the formation of some excellent characteristics, such as heat tolerance, strong disease resistance and excellent meat quality. Over the past few decades, the number of these pig populations has decreased sharply, largely due to a decrease in growth rate and poor lean meat percentage. For effective conservation of these genetic resources (such as heat tolerance, meat quality and disease resistance), the whole-genome sequencing data of 78 individuals from 3 native Chinese pig populations, including Wuzhishan (WZS), Tunchang (TC) and Dingan (DA), were obtained using a 150 bp paired-end platform, and 25 individuals from two foreign breeds, including Landrace (LR) and Large White (LW), were downloaded from a public database. A total of 28,384,282 SNPs were identified, of which 27,134,233 SNPs were identified in native Chinese pig populations. Both genetic diversity statistics and linkage disequilibrium (LD) analysis indicated that indigenous pig populations displayed high genetic diversity. The result of population structure implied the uniqueness of each native Chinese pig population. The selection signatures were detected between indigenous pig populations and foreign breeds by using the population differentiation index (F_ST) method. A total of 359 candidate genes were identified, and some genes may affect characteristics such as immunity (IL-2, IL-21 and ZFYVE16), adaptability (APBA1), reproduction (FGF2, RNF17, ADAD1 and HIPK4), meat quality (ABCA1, ADIG, TLE4 and IRX5), and heat tolerance (VPS13A, HSPA4). Overall, the findings of this study will provide some valuable insights for the future breeding, conservation and utilization of these three Chinese indigenous pig populations.

Lymphoid enhancer binding factor 1 is associated with nose color in Yunong black pigs

Yunong black pig is an indigenous black pig breed being cultivated that has a pure black whole body. However, some individuals appear with a white spot on the nose. We performed case-control association studies and F_ST approaches in 76 animals with nose color records (26 white-nosed pigs vs. 50 black-nosed pigs) by Illumina Porcine SNP50 BeadChip data. In total, 76 SNPs, which included 2 genome-wide significant SNPs and 18 chromosome-wide suggestive SNPs, were identified by association study. The top-ranked 0.1% windows of F_ST results as signals under selection and 24 windows were selected. The lymphoid enhancer binding factor 1 was identified as candidate gene with strong signal in analyses of genome-wide association study and F_ST in black- and white-nosed pigs. Overall, our findings provide evidence that nose color is a heritable trait influenced by many loci. The results contribute to expand our understanding of pigmentation in pigs and provide SNP markers for skin color and related traits selection in Yunong black pigs. Additional research on the genetic link between nose pigmentation is needed.

Distinct traces of mixed ancestry in western commercial pig genomes following gene flow from Chinese indigenous breeds

Studying gene flow between different livestock breeds will benefit the discovery of genes related to production traits and provide insight into human historical breeding. Chinese pigs have played an indispensable role in the breeding of Western commercial pigs. However, the differences in the timing and volume of the contribution of pigs from different Chinese regions to Western pigs are not yet apparent. In this paper, we combine the whole-genome sequencing data of 592 pigs from different studies and illustrate patterns of gene flow from Chinese pigs into Western commercial pigs. We describe introgression patterns from four distinct Chinese indigenous groups into five Western commercial groups. There were considerable differences in the number and length of the putative introgressed segments from Chinese pig groups that contributed to Western commercial pig breeds. The contribution of pigs from different Chinese geographical locations to a given western commercial breed varied more than that from a specific Chinese pig group to different Western commercial breeds, implying admixture within Europe after introgression. Within different Western commercial lines from the same breed, the introgression patterns from a given Chinese pig group seemed highly conserved, suggesting that introgression of Chinese pigs into Western commercial pig breeds mainly occurred at an early stage of breed formation. Finally, based on analyses of introgression signals, allele frequencies, and selection footprints, we identified a ∼2.65 Mb Chinese-derived haplotype under selection in Duroc pigs (CHR14: 95.68-98.33 Mb). Functional and phenotypic studies demonstrate that this PRKG1 haplotype is related to backfat and loin depth in Duroc pigs. Overall, we demonstrate that the introgression history of domestic pigs is complex and that Western commercial pigs contain distinct traces of mixed ancestry, likely derived from various Chinese pig breeds.

Population Structure and Selection Signatures Underlying Domestication Inferred from Genome-Wide Copy Number Variations in Chinese Indigenous Pigs

Single nucleotide polymorphism was widely used to perform genetic and evolution research in pigs. However, little is known about the effect of copy number variation (CNV) on characteristics in pigs. This study performed a genome-wide comparison of CNVs between Wannan black pigs (WBP) and Asian wild boars (AWB), using whole genome resequencing data. By using Manta, we detected in total 28,720 CNVs that covered approximately 1.98% of the pig genome length. We identified 288 selected CNVs (top 1%) by performing Fst statistics. Functional enrichment analyses for genes located in selected CNVs were found to be muscle related (NDN, TMOD4, SFRP1, and SMYD3), reproduction related (GJA1, CYP26B1, WNT5A, SRD5A2, PTPN11, SPEF2, and CCNB1), residual feed intake (RFI) related (MAP3K5), and ear size related (WIF1). This study provides essential information on selected CNVs in Wannan black pigs for further research on the genetic basis of the complex phenotypic and provides essential information for direction in the protection and utilization of Wannan black pig.

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.

Genetic architecture and selection of Anhui autochthonous pig population revealed by whole genome resequencing

The genetic resources among pigs in Anhui Province are diverse, but their value and potential have yet to be discovered. To illustrate the genetic diversity and population structure of the Anhui pigs population, we resequenced the genome of 150 pigs from six representative Anhui pigs populations and analyzed this data together with the sequencing data from 40 Asian wild boars and commercial pigs. Our results showed that Anhui pigs were divided into two distinct types based on ancestral descent: Wannan Spotted pig (WSP) and Wannan Black pig (WBP) origins from the same ancestor and the other four populations origins from another ancestor. We also identified several potential selective sweep regions associated with domestication characteristics among Anhui pigs, including reproduction-associated genes (CABS1, INSL6, MAP3K12, IGF1R, INSR, LIMK2, PATZ1, MAPK1), lipid- and meat-related genes (SNX19, MSTN, MC5R, PRKG1, CREBBP, ADCY9), and ear size genes (MSRB3 and SOX5). Therefore, these findings expand the catalogue and how these genetic differences among pigs and this newly generated data will be a valuable resource for future genetic studies and for improving genome-assisted breeding of pigs and other domesticated animals.

Genetic diversity and population structure of Tongcheng pigs in China using whole-genome SNP chip

Tongcheng (TC) pigs, distinguished by their superior meat quality, are a Chinese indigenous pig breed. Recently, the genetic resources of TC pigs are under tremendous threat due to the introduction of cosmopolitan pig breeds and African swine fever disease. To promote their management and conservation, the present study assessed genetic diversity and population structure of TC pigs using single nucleotide polymorphism (SNP) markers. A total of 26, 999 SNPs were screened from 51, 315 SNPs in 68 TC pigs. The multi-dimensional scaling (MDS) analysis and neighbor-joining tree revealed that all 68 pigs were from a purebred population. The effective population size decreased over time, and it was 96 prior to generation 20. Both linkage disequilibrium (LD) and neutrality test indicated a low selection of TC pigs with average LD value of 0.15 ± 0.23. Genetic diversity results exhibited a minor allele frequency (MAF) of 0.23, observed heterozygosity (HO) of 0.32, expected heterozygosity (He) of 0.31, and nucleotide diversity (Pi) of 0.31. All these parameters indicated a remarkably high genetic diversity of TC pigs. Additionally, 184 runs of homozygosity (ROH) segments were detected from the whole genome of TC pigs with an average ROH length of 23.71Mb, ranging from 11.26Mb to 69.02 Mb. The highest ROH coverage was found on chromosome 1 (10.12%), while the lowest was on chromosome 18 (1.49%). The average inbreeding coefficients based on ROH (FROH) was 0.04%. Fourteen ROH islands containing 240 genes were detected on 9 different autosomes. Some of these 240 genes were overlapped with the genes related to biological processes such as immune function, reproduction, muscular development, and fat deposition, including FFAR2, FFAR4, MAPK8, NPY5R, KISS1, and these genes might be associated with such traits as meat quality and disease resistance in TC pigs. Taken together, population structure and genetic diversity results suggested that the TC pig represented a valuable genetic resource. However, TC pig breed conservation program remains to be further optimized to ensure adequate genetic diversity and avoid inbreeding depression. Our findings provide theoretical basis for formulating management and conservation strategies for TC pigs.

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.

A comparative investigation on H3K27ac enhancer activities in the brain and liver tissues between wild boars and domesticated pigs

Dramatic phenotypic differences between domestic pigs and wild boars (Sus scrofa) provide opportunities to investigate molecular mechanisms underlying the formation of complex traits, including morphology, physiology and behaviour. Most studies comparing domestic pigs and wild boars have focused on variations in DNA sequences and mRNA expression, but not on epigenetic changes. Here, we present a genome‐wide comparative study on H3K27ac enhancer activities and the corresponding mRNA profiling in the brain and liver tissues of adult Bama Xiang pigs (BMXs) and Chinese wild boars (CWBs). We identified a total of 1,29,487 potential regulatory elements, among which 11,241 H3K27ac peaks showed differential activity between CWBs and BMXs in at least one tissue. These peaks were overrepresented by binding motifs of FOXA1, JunB, ATF3 and BATF, and overlapped with differentially expressed genes that are involved in female mating behaviour, response to growth factors and hormones, and lipid metabolism. We also identified 4118 nonredundant super‐enhancers from ChIP‐Seq data on H3K27ac. Notably, we identified differentially active peaks located close to or within candidate genes, including TBX19, MSTN, AHR and P2RY1, which were identified in DNA sequence‐based population differentiation studies. This study generates a valuable dataset on H3K27ac profiles of the brain and liver from domestic pigs and wild boars, which helps gain insights into the changes in enhancer activities from wild boars to domestic pigs.

Machine-Learning Prospects for Detecting Selection Signatures Using Population Genomics Data

Natural selection has been given a lot of attention because it relates to the adaptation of populations to their environments, both biotic and abiotic. An allele is selected when it is favored by natural selection. Consequently, the favored allele increases in frequency in the population and neighboring linked variation diminishes, causing so-called selective sweeps. A high-throughput genomic sequence allows one to disentangle the evolutionary forces at play in populations. With the development of high-throughput genome sequencing technologies, it has become easier to detect these selective sweeps/selection signatures. Various methods can be used to detect selective sweeps, from simple implementations using summary statistics to complex statistical approaches. One of the important problems of these statistical models is the potential to provide inaccurate results when their assumptions are violated. The use of machine learning (ML) in population genetics has been introduced as an alternative method of detecting selection by treating the problem of detecting selection signatures as a classification problem. Since the availability of population genomics data is increasing, researchers may incorporate ML into these statistical models to infer signatures of selection with higher predictive accuracy and better resolution. This article describes how ML can be used to aid in detecting and studying natural selection patterns using population genomic data.

The de novo assembly of a European wild boar genome revealed unique patterns of chromosomal structural variations and segmental duplications

The rapid progress of sequencing technology has greatly facilitated the de novo genome assembly of pig breeds. However, the assembly of the wild boar genome is still lacking, hampering our understanding of chromosomal and genomic evolution during domestication from wild boars into domestic pigs. Here, we sequenced and de novo assembled a European wild boar genome (ASM2165605v1) using the long‐range information provided by 10× Linked‐Reads sequencing. We achieved a high‐quality assembly with contig N50 of 26.09 Mb. Additionally, 1.64% of the contigs (222) with lengths from 107.65 kb to 75.36 Mb covered 90.3% of the total genome size of ASM2165605v1 (~2.5 Gb). Mapping analysis revealed that the contigs can fill 24.73% (93/376) of the gaps present in the orthologous regions of the updated pig reference genome (Sscrofa11.1). We further improved the contigs into chromosome level with a reference‐assistant scaffolding method. Using the ‘assembly‐to‐assembly’ approach, we identified intra‐chromosomal large structural variations (SVs, length >1 kb) between ASM2165605v1 and Sscrofa11.1 assemblies. Interestingly, we found that the number of SV events on the X chromosome deviated significantly from the linear models fitting autosomes (R² > 0.64, p < 0.001). Specifically, deletions and insertions were deficient on the X chromosome by 66.14 and 58.41% respectively, whereas duplications and inversions were excessive on the X chromosome by 71.96 and 107.61% respectively. We further used the large segmental duplications (SDs, >1 kb) events as a proxy to understand the large‐scale inter‐chromosomal evolution, by resolving parental‐derived relationships for SD pairs. We revealed a significant excess of SD movements from the X chromosome to autosomes (p < 0.001), consistent with the expectation of meiotic sex chromosome inactivation. Enrichment analyses indicated that the genes within derived SD copies on autosomes were significantly related to biological processes involving nervous system, lipid biosynthesis and sperm motility (p < 0.01). Together, our analyses of the de novo assembly of ASM2165605v1 provides insight into the SVs between European wild boar and domestic pig, in addition to the ongoing process of meiotic sex chromosome inactivation in driving inter‐chromosomal interaction between the sex chromosome and autosomes.

Long-read assembly of the Chinese indigenous Ningxiang pig genome and identification of genetic variations in fat metabolism among different breeds

Advances in long-read sequencing technology and genome assembly provide an opportunity to improve the pig genome and reveal the full range of structural variations (SVs) between local Chinese and European pigs. To date, little is known about the genomes of some unique Chinese indigenous breeds, such as the Ningxiang pig. Here, we report the sequencing and assembly of a highly contiguous Ningxiang pig genome (NX) via an integration of PacBio single-molecule real-time sequencing, Illumina next-generation sequencing, BioNano optical mapping and Hi-C (chromosome conformation capture) approaches. The assembled genome comprises 2.44 Gb with a contig N50 of 26.1 Mb and 418 contigs in total. These contigs are organized into 121 scaffolds with a scaffold N50 of 139.0 Mb. More than 99.1% of the assembled sequence could be localized to 19 pseudochromosomes and is annotated with 20,914 protein-coding genes and 34.04% repetitive sequences. Comparisons between the NX and European Duroc assemblies revealed many SVs in genes involved in the immune system, nervous system, lipid metabolism and environmental adaptation. The genetic variants include 47 Chinese domestic pig-specific SVs and the associated 74 genes may contribute to the differences in domestic traits compared to European pigs. Moreover, single nucleotide polymorphisms (SNPs) identified from whole genome resequencing data of 73 Chinese pigs, representing 17 geographically isolated breeds, showed their specific genetic variations, population structure and evolutionary patterns. Finally, we explore transcriptional regulation in the first intron of the MYL4 gene, as the genomic SV (281-bp deletion) in Ningxiang pig promotes its subcutaneous fat compared to European pig breeds. This work identifies a set of Asian-specific SVs and SNPs, which will be important resources for modern pig breeding and genetic conservation.

PRE-1 Revealed Previous Unknown Introgression Events in Eurasian Boars during the Middle Pleistocene

Introgression events and population admixture occurred among Sus species across the Eurasian mainland in the Middle Pleistocene, which reflects the local adaption of different populations and contributes to evolutionary novelty. Previous findings on these population introgressions were largely based on extensive genome-wide single-nucleotide polymorphism information, ignoring structural variants (SVs) as an important alternative resource of genetic variations. Here, we profiled the genome-wide SVs and explored the formation of pattern-related SVs, indicating that PRE1-SS is a recently active subfamily that was strongly associated with introgression events in multiple Asian and European pig populations. As reflected by the three different combination haplotypes from two specific patterns and known phylogenetic relationships in Eurasian boars, we identified the Asian Northern wild pigs as having experienced introgression from European wild boars around 0.5–0.2 Ma and having received latitude-related selection. During further exploration of the influence of pattern-related SVs on gene functions, we found substantial sequence changes in 199 intron regions of 54 genes and 3 exon regions of 3 genes (HDX, TRO, and SMIM1), implying that the pattern-related SVs were highly related to positive selection and adaption of pigs. Our findings revealed novel introgression events in Eurasian wild boars, providing a timeline of population admixture and divergence across the Eurasian mainland in the Middle Pleistocene.

Pooled Sequencing Analysis of Geese (Anser cygnoides) Reveals Genomic Variations Associated With Feather Color

During the domestication of the goose a change in its feather color took place, however, the molecular mechanisms responsible for this change are not completely understood. Here, we performed whole-genome resequencing on three pooled samples of geese (feral and domestic geese), with two distinct feather colors, to identify genes that might regulate feather color. We identified around 8 million SNPs within each of the three pools and validated allele frequencies for a subset of these SNPs using PCR and Sanger sequencing. Several genomic regions with signatures of differential selection were found when we compared the gray and white feather color populations using the F_ST and Hp approaches. When we combined previous functional studies with our genomic analyses we identified 26 genes (KITLG, MITF, TYRO3, KIT, AP3B1, SMARCA2, ROR2, CSNK1G3, CCDC112, VAMP7, SLC16A2, LOC106047519, RLIM, KIAA2022, ST8SIA4, LOC106044163, TRPM6, TICAM2, LOC106038556, LOC106038575, LOC106038574, LOC106038594, LOC106038573, LOC106038604, LOC106047489, and LOC106047492) that potentially regulate feather color in geese. These results substantially expand the catalog of potential feather color regulators in geese and provide a basis for further studies on domestication and avian feather coloration.

Genomics of Adaptations in Ungulates

Ungulates are a group of hoofed animals that have long interacted with humans as essential sources of food, labor, clothing, and transportation. These consist of domesticated, feral, and wild species raised in a wide range of habitats and biomes. Given the diverse and extreme environments inhabited by ungulates, unique adaptive traits are fundamental for fitness. The documentation of genes that underlie their genomic signatures of selection is crucial in this regard. The increasing availability of advanced sequencing technologies has seen the rapid growth of ungulate genomic resources, which offers an exceptional opportunity to understand their adaptive evolution. Here, we summarize the current knowledge on evolutionary genetic signatures underlying the adaptations of ungulates to different habitats.

Identifying the unique characteristics of the Chinese indigenous pig breeds in the Yangtze River Delta region for precise conservation

Background

China is the country with the most abundant swine genetic resources in the world. Through thousands of years of domestication and natural selection, most of pigs in China have developed unique genetic characteristics. Finding the unique genetic characteristics and modules of each breed is an essential part of their precise conservation.

Results

In this study, we used the partial least squares method to identify the significant specific SNPs of 19 local Chinese pig breeds and 5 Western pig breeds. A total of 37,514 significant specific SNPs (p < 0.01) were obtained from these breeds, and the Chinese local pig breed with the most significant SNPs was Hongdenglong (HD), followed by Jiaxing black (JX), Huaibei (HB), Bihu (BH), small Meishan (SMS), Shengxian Hua (SH), Jiangquhai (JQ), Mi (MI), Chunan (CA), Chalu (CL), Jinhualiangtouwu (JHL), Fengjing (FJ), middle Meishan (MMS), Shanzhu (SZ), Pudong white (PD), Dongchuan (DC), Erhualian (EH), Shawutou (SW) and Lanxi Hua (LX) pig.

Furthermore, we identified the breeds with the most significant genes, GO terms, pathways, and networks using KOBAS and IPA and then ranked them separately. The results showed that the breeds with the highest number of interaction networks were Hongdenglong (12) and Huaibei (12) pigs. In contrast, the breeds with the lowest interaction networks were Shawutou (4) and Lanxi Hua pigs (3), indicating that Hongdenglong and Huaibei pigs might have the most significant genetic modules in their genome, whereas Shawutou and Lanxi Hua pigs may have the least unique characteristics. To some degree, the identified specific pathways and networks are related to the number of genes and SNPs linked to the specific breeds, but they do not appear to be the same. Most importantly, more significant modules were found to be related to the development and function of the digestive system, regulation of diseases, and metabolism of amino acids in the local Chinese pig breeds, whereas more significant modules were found to be related to the growth rate in the Western pig breeds.

Conclusion

Our results show that each breed has some relatively unique structural modules and functional characteristics. These modules allow us to better understand the genetic differences among local Chinese and Western pig breeds and therefore implement precise conservation methods. This study could provide a basis for formulating more effective strategies for managing and protecting these genetic resources in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07476-7.

Phenotypic and genomic relationships between vulva score categories and reproductive performance in first-parity sows

Background

One of the biggest challenges in the swine industry is to increase female reproductive efficiency. Recently, vulva score categories (VSC), assessed prior to puberty, has been proposed as an indicator trait of efficient reproductive performance in sows. The objective of this study was to validate the use of VSC as an indicator trait for reproductive performance, and to perform genetic and genomic analyses for VSC.

Methods

The phenotypic relationship of VSC, using a three-point scale: small (VSC-S), medium (VSC-M), and large (VSC-L), on reproductive performance was evaluated on three farms. VSC was measured at 15 weeks of age, for farms 1 and 2, and at 14 weeks of age for farm 3 on 3981 Yorkshire gilts, in which 1083 had genotypes (~ 50 K SNPs). Genetic parameters for VSC with reproductive traits were estimated using ssGBLUP. A Genome-wide association study (GWAS) for VSC was performed using BayesB.

Results

For the phenotypic analysis of VSC across datasets, differences in performance were identified there was a significant effect (P ≤ 0.05) for the interaction between Farm and VSC for total number dead (TND), and a trend (P < 0.10) for total number born (TNB). There were significant (P ≤ 0.05) pre-defined contrasts of VSC-S versus VSC-M + L on TNB, number born alive (NBA), TND, number of stillborn (NSB), and number of mummies (MUM). Heritability estimates for VSC as a categorical trait (VSCc) and a quantitative trait (VSCq) were 0.40 ± 0.02 and 0.83 ± 0.02, respectively, for across farm, 0.13 ± 0.07 and 0.20 ± 0.10, respectively, for Farm1, 0.07 ± 0.07 and 0.09 ± 0.09, respectively, for Farm2, and 0.20 ± 0.03 and 0.34 ± 0.05, respectively, for Farm3. For across farms, favorable genetic correlations estimates were found for TNB (0.28 ± 0.19) and NBA (0.26 ± 0.17). Within farms, moderate genetic correlations between VSC with reproductive traits were found for TNB (0.61 ± 0.47) and MUM (0.69 ± 0.47) for farm 1, for number of services until first farrow (NS; 0.69 ± 0.38) and unique service with successful first farrow (SFS; − 0.71 ± 0.38) for farm 3. Multiple genomic regions associated with VSC_c were identified. Of these, a QTL located on chromosome 3 at 33–34 Mb accounted for about 7.1% of the genetic variance for VSC_c and VSC_q. This region harbors the gene PRM1 that has been associated with early embryonic development in pigs.

Conclusions

The results support potential of VSC for improved reproductive efficiency on first-parity performance, but the results might depend on the interaction between environmental factors and VSC, as well as potentially additive genetics.

Red and blond Mangalitza pigs display a signature of divergent directional selection in the SLC45A2 gene

The Mangalitza lard-type pig breed is well known for its fat appearance and curly hair, and it is mainly distributed in Eastern Europe. Four main lines were created in the nineteenth century by artificial selection: Blond Mangalitza, Black Mangalitza, Swallow-Belly Mangalitza and Red Mangalitza. The Swallow-Belly line has a black coat combined with yellow-blond throat and underbelly. In the current work, we aimed to investigate if the colourations of Mangalitza pigs are genetically determined by one or a few loci whose frequencies have been modified by artificial selection. The results of selection scans, with HapFLK and BayeScan, and of a GWAS for coat colour highlighted the existence of one region on SSC16 (18-20 Mb) with potential effects on hair pigmentation (Red vs. Blond contrast). The analysis of the gene content of this region allowed us to detect the solute carrier family 45 member 2 (SLC45A2) locus as a candidate gene for this trait. The polymorphism of the SLC45A2 locus has been associated with reduced levels or the absence of melanin in several mammalian species. The genotyping of four missense polymorphisms evidenced that rs341599992:G > A and rs693695020:G > A SNPs are strongly but not fully associated with the red and blond coat colours of Mangalitza pigs, a result that was confirmed by performing a haplotype association test. The near fixation of alternative SLC45A2 genotypes in Red and Blond Mangalitza pigs provides a compelling example of the consequences of a divergent directional selection for coat colour in a domestic species.

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.

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.

Whole-genome sequencing reveals breed-differential CNVs between Tongcheng and Large White pigs

Large phenotypic differences have been observed between Tongcheng and Large White pigs. However, little is known about their genetic basis. This study performed a genome-wide comparison of CNVs between Tongcheng and Large White pigs using genome sequencing data. By combining the advantages of three different strategies (read depth, paired-end mapping and split read), we detected in total 18 687 CNVs that covered approximately 3.5% of the pig genome length for Tongcheng and Large White pigs. We identified 1864 breed-stratified CNVs (top 10%) by performing V_ST statistics. Functional enrichment analyses for genes located in breed-stratified CNVs were found to be involved in pigmentation, behavior, immune system and reproductive processes, which coincide with phenotypic differences between the two breeds. Using a systematic analysis of the genome and transcriptome data, we further identified four novel breed-differential CNVs on the functional genes (disease-resistant, DCUN1D2 and SPARCL1; lipid metabolism, PLEKHA2 and SLCO1A2). Subsequent PCR validation confirmed their accurate breakpoint positions in 33 Tongcheng pigs and 33 Large White pigs. This study provides essential information on differential CNVs for further research on the genetic basis of phenotypic differences between Tongcheng and Large White pigs.

An integration analysis based on genomic, transcriptomic and QTX information reveals credible candidate genes for fat-related traits in pigs

Meat quality improvement is of great interest to researchers in pig breeding and many researchers have identified abundant associated quantitative trait loci, genes and polymorphisms (QTXs) for fat-related traits. However, it is challenging to determine credible candidate genes from a mass of associations. The efficiency of identification of credible candidate genes in these QTXs is restricted by limited integration analyses of data from multiple omics. In this study, we constructed a 'candidate gene map' of fat-related traits in pigs based on published literature and the latest genome. In total, 6,861 QTXs, which covered 9,323 genes on the pig genome, were used. Combining the QTX hotspots and pathway analysis, we identified 180 candidate genes that may regulate the fat-related traits, and choose PNPLA2, PPARG, SREBF1, ACACA, PPARD and PPARA as credible candidate genes. In addition, we discussed the importance of incorporating transcriptome data and genomic data in causal gene identification, and the multi-omics information can effectively improve the credibility of identified candidate genes.

Whole genome sequence analysis reveals genetic structure and X-chromosome haplotype structure in indigenous Chinese pigs

Chinese indigenous pigs exhibit considerable phenotypic diversity, but their population structure and the genetic basis of agriculturally important traits need further exploration. Here, we sequenced the whole genomes of 24 individual pigs representing 22 breeds distributed throughout China. For comparison with European and commercial breeds (one pig per breed), we included seven published pig genomes with our new genomes for analyses. Our results showed that breeds grouped together based on morphological classifications are not necessarily more genetically similar to each other than to breeds from other groups. We found that genetic material from European pigs likely introgressed into five Chinese breeds. We have identified two new subpopulations of domestic pigs that encompass morphology-based criteria in China. The Southern Chinese subpopulation comprises the classical South Chinese Type and part of the Central China Type. In contrast, the Northern Chinese subpopulation comprises the North China Type, the Lower Yangtze River Basin Type, the Southwest Type, the Plateau Type, and the remainder of the Central China Type. Eight haplotypes and two recombination sites were identified within a conserved 40.09 Mb linkage-disequilibrium (LD) block on the X chromosome. Potential candidate genes (LEPR, FANCC, COL1A1, and PCCA) influencing body size were identified. Our findings provide insights into the phylogeny of Chinese indigenous pig breeds and benefit gene mining efforts to improve major economic traits.

Polymorphisms in 5' proximal regulating region of THRSP gene are associated with fat production in pigs

Chinese and imported pig breeds differ in fat production potential, which is associated with the polymorphisms in the 5' proximal regulating region (5'PRR) of thyroid hormone responsive gene (THRSP). In three Chinese breeds (Dingyuan, CDY; Wannanhua, CWH; and Jixi, CJX) and one introduced breed (Yorkshire, YKS), three variant sites were located at T/C-400, A/G-376, and G/A-98 in the 5'PRR. Chinese pig breeds had higher C-400 allele frequencies than YKS. The frequencies of A-376 in CDY and G-376 in CWH were about 0.8. G-98 allele frequencies in CWH and YKS were 0.8617 and 0.8149, respectively. TGG was the dominant haplotype in YKS, CGG in CWH and CJX, and CAA in CDY. According to haplotype frequency, four breeds were clustered into three types, which was consistent with the geographical distribution of the breeds. In CDY, the average backfat thickness (BFT) was the highest with the CC-400 genotype, followed by CT-400 and TT-400 genotypes. In YKS, the pigs with CC-400 or CT-400 genotypes had higher BFT and average daily weight gain, whereas those with CC-400 or TT-400 genotypes had larger lion-eye area. No significant difference was observed in carcass traits among different genotypes at the A/G-376 and G/A-98 loci. The mRNA abundance of THRSP expression for the CCAGAG genotype was significantly higher than that for CTAGAG or TTAGAG genotype. These results indicated that the polymorphisms and genotype distribution of THRSP were closely related to the potential for fat production in pig breeds, which were the result of adaptation to artificial selection and natural selection.

Genome-wide association studies for iris pigmentation and heterochromia patterns in Large White pigs

Eye colour genetics have been extensively studied in humans since the rediscovery of Mendel's laws. This trait was first interpreted using simplistic genetic models but soon it was realised that it is more complex. In this study, we analysed eye colour variability in a Large White pig population (n = 897) and report the results of GWASs based on several comparisons including pigs having four main eye colour categories (three with both pigmented eyes of different brown grades: pale, 17.9%; medium, 14.8%; and dark, 54.3%; another one with both eyes completely depigmented, 3.8%) and heterochromia patterns (heterochromia iridis - depigmented iris sectors in pigmented irises, 3.2%; heterochromia iridum - one whole eye iris of depigmented phenotype and the other eye with the iris completely pigmented, 5.9%). Pigs were genotyped with the Illumina PorcineSNP60 BeadChip and GEMMA was used for the association analyses. The results indicated that SLC45A2 (on chromosome 16, SSC16), EDNRB (SSC11) and KITLG (SSC5) affect the different grades of brown pigmentation of the eyes, the bilateral eye depigmentation defect and the heterochromia iridis defect recorded in this white pig population respectively. These genes are involved in several mechanisms affecting pigmentation. Significant associations for the eye depigmented patterns were also identified for SNPs on two SSC4 regions (including two candidate genes: NOTCH2 and PREX2) and on SSC6, SSC8 and SSC14 (including COL17A1 as candidate gene). This study provided useful information to understand eye pigmentation mechanisms, further valuing the pig as animal model to study complex phenotypes in humans.

Genome-Wide Analysis Reveals Human-Mediated Introgression from Western Pigs to Indigenous Chinese Breeds

Genetic variations introduced via introgression from Western to Chinese pigs have contributed to the performance of Chinese breeds in traits such as growth rate and feed conversion efficiency. However, little is known about the underlying genomic changes that occurred during introgression and the types of traits affected by introgression. To address these questions, 525 animals were characterized using an SNP array to detect genomic regions that had been introgressed from European to indigenous Chinese breeds. The functions of genes located in introgressed regions were also investigated. Our data show that five out of six indigenous Chinese breeds show evidence of introgression from Western pigs, and eight introgressed genome regions are shared by five of the Chinese breeds. A region located on chr13: 12.8-13.1 M was affected by both introgression and artificial selection, and this region contains the glucose absorption related gene, OXSM, and the sensory related gene, NGLY. The results provide a foundation for understanding introgression from Western to indigenous Chinese pigs.

The fine-scale genetic structure and selection signals of Chinese indigenous pigs

Genome-wide SNP profiling has yielded insights into the genetic structure of China indigenous pigs, but has focused on a limited number of populations. Here, we present an analysis of population structure and signals of positive selection in 42 Chinese pig populations that represent the most extensive pig phenotypic diversity in China, using genotype data of 1.1 million SNPs on customized Beadchips. This unravels the fine-scale genetic diversity, phylogenic relationships, and population structure of these populations, which shows remarkably concordance between genetic clusters and geography with few exceptions. We also reveal the genetic contribution to North Chinese pigs from European modern pigs. Furthermore, we identify possible targets of selection in the Tibetan pig, including the well-characterized hypoxia gene (EPAS1) and several previously unrecognized candidates. Intriguingly, the selected haplotype in the EPAS1 gene is associated with higher hemoglobin contents in Tibetan pigs, which is different from the protective role of EPAS1 in the high-altitude adaptation in Tibetan dogs and their owners. Additionally, we present evidence for the causality between EDNRB variants and the two-end-black (TEB) coat color phenotype in all Chinese pig populations except the Jinhua pig. We hypothesize that distinct targets have been independently selected for the formation of the TEB phenotype in Chinese pigs of different geographic origins. This highlights the importance of characterizing population-specific genetic determinants for heritable phenotype in diverse pig populations.

Association analysis of mitochondrial DNA polymorphisms with oocyte number in pigs

In pigs, correlations between mitochondrial (mt) DNA polymorphisms and economic traits have been widely reported across and within swine breeds. In fecundity studies, the number of oocytes within ovaries was highly correlated with litter size. However, the effect of mitochondrial polymorphisms on porcine oocyte number remained unclear. In this study, 181 porcine ovaries were collected to analyse the relationship between oocyte number and mtDNA polymorphisms. There were considerable differences in oocyte numbers among different ovaries from commercial pig breeds, ranging from 2.7×105 to 1.3×106. Mitochondrial D-loop sequencing discovered 53 polymorphic sites. Association analysis revealed that 13 variations were associated with the number of oocytes (P<0.05). A C323T polymorphism showed the largest value between the C and T carriers, which differed at 105 oocytes (P<0.05). The 53 polymorphic sites generated 45 haplotypes, which clustered into two haplogroups, A and B. Haplogroup A had a higher number of oocytes than Haplogroup B (P<0.05), whereas Haplotype H6 in Haplogroup A had the highest number of oocytes (~7.5×105) of all haplotypes studied (P<0.05). The results of this study highlight a correlation between mtDNA polymorphisms and oocyte number, and suggest the potential application of mtDNA polymorphism analyses in pig selection and breeding practices.

Integrated Analysis of mRNA Expression, CpG Island Methylation, and Polymorphisms in the MITF Gene in Ducks (Anas platyrhynchos)

Microphthalmia-associated transcription factor (MITF) is a key regulator for the development and function of melanocytes in skin, eye, and plumage pigmentations. Thus, the MITF was selected as a candidate gene associated with plumage coloration in ducks. This study analyzed the mRNA expression, promoter methylation, and polymorphisms in the MITF gene in ducks with different plumage colors (Putian Black, Putian White, Liancheng White, and Longsheng Jade-green). No expression of the MITF melanin-specific isoform (MITF-M) was detected in white feather bulbs. By contrast, the mRNA expression levels of MITF-M were high in black feather bulbs. Bioinformatics analysis showed that two CpG islands were present in the promoter region of the MITF gene. The methylation level of the second CpG island was significantly lower in black feather bulbs than in white feather bulbs. However, the methylation level of the first CpG island was not different among the feather bulbs with various colors except Liancheng White feather bulbs. The methylation status of the whole CpG island significantly and negatively correlated with the mRNA expression of MITF-M (P < 0.05). Furthermore, four novel SNPs (single nucleotide polymorphisms) were identified in the 5′UTR, exon 4, intron 7, and intron 8 of the MITF gene. Allele T in g.39807T>G and allele G in g.40862G>A were the predominant alleles only found in Putian White, whereas the variant A allele in g.32813G>A exhibited a high allele frequency in Liancheng White. Collectively, these results contributed to the understanding of the function of the MITF gene in duck plumage coloration.

A genome-wide association study of reproduction traits in four pig populations with different genetic backgrounds

Objective

Genome-wide association study and two meta-analysis based on GWAS performed to explore the genetic mechanism underlying variation in pig number born alive (NBA) and total number born (TNB).

Methods

Single trait GWAS and two meta-analysis (single-trait meta analysis and multi-trait meta analysis) were used in our study for NBA and TNB on 3,121 Yorkshires from 4 populations, including three different American Yorkshire populations (n = 2,247) and one British Yorkshire populations (n = 874).

Results

The result of single trait GWAS showed that no significant associated single nucleotide polymorphisms (SNPs) were identified. Using single-trait meta analysis and multi-trait meta analysis within populations, 11 significant loci were identified associated with target traits. Spindlin 1, vascular endothelial growth factor A, forkhead box Q1, msh homeobox 1, and LHFPL tetraspan submily member 3 are five functionally plausible candidate genes for NBA and TNB. Compared to the single population GWAS, single-trait Meta analysis can improve the detection power to identify SNPs by integrating information of multiple populations. The multiple-trait analysis reduced the power to detect trait-specific loci but enhanced the power to identify the common loci across traits.

Conclusion

In total, our findings identified novel genes to be validated as candidates for NBA and TNB in pigs. Also, it enabled us to enlarge population size by including multiple populations with different genetic backgrounds and increase the power of GWAS by using meta analysis.

Genome-wide detection of genetic loci and candidate genes for teat number and body conformation traits at birth in Chinese Sushan pigs

Body conformation at birth and teat number are economically important traits in the pig industry, as these traits are usually explored to evaluate the growth and reproductive potential of piglets. To detect genetic loci and candidate genes for these traits, we performed a GWAS on 269 pigs from a recently developed Chinese breed (Sushan) using 38  128 informative SNPs on the Affymetrix Porcine SNP 55K Array. In total, we detected one genome-wide significant (P = 1.31e-6) SNP for teat number on chromosome X and 15 chromosome-wide significant SNPs for teat number, body weight, body length, chest circumference and cannon circumference at birth on chromosomes 1, 3, 4, 6, 7, 9, 10, 13, 14, 15, 17 and 18. The most significant SNP had an additive effect of 0.74 × total teat number, explaining 20% of phenotypic variance. Five significant SNPs resided in the previously reported quantitative trait loci for these traits and seven significant SNPs had a pleiotropic effect on multiple traits. Intriguingly, 12 of the genes nearest to the significant SNPs are functionally related to body conformation and teat number traits, including SPRED2, MKX, TMSB4X and ESR1. GO analysis revealed that candidate genes proximal to the significant SNPs were enriched in the G-protein coupled receptor and steroid hormone-mediated signaling pathway. Our findings shed light on the genetic basis of the measured traits and provide molecular markers especially for the genetic improvement of teat number in Sushan and related pigs.

Deep Genome Resequencing Reveals Artificial and Natural Selection for Visual Deterioration, Plateau Adaptability and High Prolificacy in Chinese Domestic Sheep

Sheep were one of the earliest domesticated animals. Both artificial and natural selection during domestication has resulted in remarkable changes in behavioral, physiological, and morphological phenotypes; however, the genetic mechanisms underpinning these changes remain unclear, particularly for indigenous Chinese sheep. In the present study, we performed pooled whole-genome resequencing of 338 sheep from five breeds representative of indigenous Chinese breeds and compared them to the wild ancestors of domestic sheep (Asian mouflon, Ovis orientalis) for detection of genome-wide selective sweeps. Comparative genomic analysis between domestic sheep and Asian mouflon showed that selected regions were enriched for genes involved in bone morphogenesis, growth regulation, and embryonic and neural development in domestic sheep. Moreover, we identified several vision-associated genes with funtional mutations, such as PDE6B (c.G2994C/p.A982P and c.C2284A/p.L762M mutations), PANK2, and FOXC1/GMSD in all five Chinese native breeds. Breed-specific selected regions were determined including genes such as CYP17 for hypoxia adaptability in Tibetan sheep and DNAJB5 for heat tolerance in Duolang sheep. Our findings provide insights into the genetic mechanisms underlying important phenotypic changes that have occurred during sheep domestication and subsequent selection.

Genomic structure of a crossbred Landrace pig population

Single nucleotide polymorphism (SNP) markers are used to study population structure and conservation genetics, which permits assessing similarities regarding the linkage disequilibrium and information about the relationship among individuals. To investigate the population genomic structure of 300 females and 25 males from a commercial maternal pig line we analyzed linkage disequilibrium extent, inbreeding coefficients using genomic and conventional pedigree data, and population stratification. The average linkage disequilibrium (r²) was 0.291 ± 0.312 for all adjacent SNPs, distancing less than 100 Kb (kilobase) between markers. The average inbreeding coefficients obtained from runs of homozygosity (ROH) and pedigree analyses were 0.119 and 0.0001, respectively. Low correlation was observed between the inbreeding coefficients possibly as a result of genetic recombination effect accounted for the ROH estimates or caused by pedigree identification errors. A large number of long ROHs might indicate recent inbreeding events in the studied population. A total of 36 homozygous segments were found in more than 30% of the population and these ROH harbor genes associated with reproductive traits. The population stratification analysis indicated that this population was possibly originated from two distinct populations, which is a result from crossings between the eastern and western breeds used in the formation of the line. Our findings provide support to understand the genetic structure of swine populations and may assist breeding companies to avoid a high level of inbreeding coefficients to maintain genetic diversity, showing the effectiveness of using genome-wide SNP information for quantifying inbreeding when the pedigree was incomplete or incorrect.

Introgression of Chinese haplotypes contributed to the improvement of Danish Duroc pigs

The distribution of Asian ancestry in the genome of Danish Duroc pigs was investigated using whole-genome sequencing data from European wild boars, Danish Duroc, Chinese Meishan and Bamaxiang pigs. Asian haplotypes deriving from Meishan and Bamaxiang occur widely across the genome. Signatures of selection on Asian haplotypes are common in the genome, but few of these haplotypes have been fixed. By defining 50-kb windows with more than 50% Chinese ancestry, which did not exhibit extreme genetic differentiation between Meishan and Bamaxiang as candidate regions, the enrichment of quantitative trait loci in candidate regions supports that Asian haplotypes under selection play an important role in contributing genetic variation underlying production, reproduction, meat and carcass, and exterior traits. Gene annotation of regions with the highest proportion of Chinese ancestry revealed genes of biological interest, such as NR6A1. Further haplotype clustering analysis suggested that a haplotype of Chinese origin around the NR6A1 gene was introduced to Europe and then underwent a selective sweep in European pigs. Besides, functional genes in candidate regions, such as AHR and PGRMC2, associated with fertility, and SAL1, associated with meat quality, were identified. Our results demonstrate the contribution of Asian haplotypes to the genomes of European pigs. Findings herein facilitate further genomic studies such as genomewide association study and genomic prediction by providing ancestry information of variants.

Identifying SNPs and candidate genes for three litter traits using single-step GWAS across six parities in Landrace and Large White pigs

Total number born (TNB), number born alive (NBA), and litter weight born alive (LWB) are critically important traits in pig production. The sow's parity is one of the major factors influencing litter traits. Because of monogenic or polygenic contributions and the presence of temporal gene effects in different sows' parities, it is difficult to clarify the biological and genetic background. To systematically explore the genetic mechanism of litter traits, we conducted 18 GWASs using single-step GWAS (ssGWAS) based on two breeds (908 Landrace and 1,130 Large White sow litter records) for each litter trait in different parities. A total of 300 Landrace and 300 Large White sows were genotyped by sequencing (GBS). ssGWAS was performed separately for each breed and each parity due to population stratification and temporal gene effect. In summary, we identified 80 (15 for Landrace and 65 for Large White), 227 (52 for Landrace, 175 for Large White), and 187 (34 for Landrace, 153 for Large White) single nucleotide polymorphisms (SNPs) affecting TNB, NBA, and LWB, respectively. Of them, we suggest that a total of 22 loci (SSC1: 125098202, SSC1: 117560058, SSC14: 147794697, SSC8: 84823302, SSC9: 143554876, and SSC9: 138766097 for Landrace; SSC1: 4023577, SSC1: 3859573, SSC1: 4891063, SSC16: 5197665, SSC10: 32050819, SSC13: 13552924, SSC13: 92819, SSC17: 3579607, SSC13: 196698221, SSC7: 30918403, SSC16: 46221484, SSC16: 46169204, SSC2: 41988642, SSC2: 44475457, SSC2: 42521875, and SSC7: 58411951 for Large White) are shared by TNB, NBA, and LWB. These results indicate the existence of gene temporal effect in each parity. Furthermore, our findings suggest four interesting candidate genes (FBXL7, ALDH1A2, LEPR, and DDX1) associated with litter traits in different parities that have a major effect on embryonic development progression. In conclusion, 22 crucial SNPs and four interesting candidate genes were identified for three litter traits across six parities. These findings advance our understanding of the genetic architecture of litter traits and confirm the presence of temporal gene effects in different parities. Importantly, functional validation studies for findings of particular interest are recommended in litter traits.