Genomic Analysis To Identify Signatures of Artificial Selection and Loci Associated with Important Economic Traits in Duroc Pigs

An intronic SNP affects skeletal muscle development by regulating the expression of TP63

Background

Porcine skeletal muscle development is pivotal for improving meat production. TP63, a transcription factor, regulates vital cellular processes, yet its role in skeletal muscle proliferation is unclear.

Methods

The effects of TP63 on skeletal muscle cell viability and proliferation were investigated using both mouse and porcine skeletal muscle myoblasts. Selective sweep analysis in Western pigs identified TP63 as a potential candidate gene for skeletal muscle development. The correlation between TP63 overexpression and cell proliferation was assessed using quantitative real-time PCR (RT-qPCR) and 5-ethynyl-2'-deoxyuridine (EDU).

Results

The study revealed a positive correlation between TP63 overexpression and skeletal muscle cell proliferation. Bioinformatics analysis predicted an interaction between MEF2A, another transcription factor, and the mutation site of TP63. Experimental validation through dual-luciferase assays confirmed that a candidate enhancer SNP could influence MEF2A binding, subsequently regulating TP63 expression and promoting skeletal muscle cell proliferation.

Conclusion

These findings offer experimental evidence for further exploration of skeletal muscle development mechanisms and the advancement of genetic breeding strategies aimed at improving meat production traits.

Genetic evaluation for production and body size traits using different animal models in purebred-Duroc pigs

Duroc pigs are popular crossbred terminal sires, and accurate assessment of genetic parameters in the population can help to rationalize breeding programmes. The principle aim of this study were to evaluate the genetic parameters of production (birth weight, BW; age at 115 kg, AGE; feed conversion ratio, FCR) and body size (body length, BL; body height, BH; front cannon circumference, FCC) traits of Duroc pigs. The second objective was to analyze the fit of different genetic assessment models. The variance components and correlations of BW (28,348 records), AGE (28,335 records), FCR (11,135 records), BL (31,544 records), BH (21,862 records), and FCC (14,684 records) traits were calculated by using DMU and AIREMLF90 from BLUPF90 package. In the common environment model, the heritability of BW, AGE, FCR, BL, BH, and FCC traits were 0.17 ± 0.014, 0.30 ± 0.019, 0.28 ± 0.024, 0.16 ± 0.013, 0.14 ± 0.017, and 0.081 ± 0.016, with common litter effect values of 0.25, 0.20, 0.18, 0.23, 0.19, and 0.16, respectively. According to the results of the Akaike information criterion (AIC) calculations, models with smaller AIC values have a better fit. We found that the common environment model with litter effects as random effects for estimating genetic parameters had a better fit. In this Model, the estimated genetic correlations between AGE with BW, FCR, BL, BH, and FCC traits were -0.28 (0.040), 0.76 (0.038), -0.71 (0.036), -0.44 (0.060), and -0.60 (0.073), respectively, with phenotypic correlations of -0.17, 0.52, -0.22, -0.13 and -0.24, respectively. In our analysis of genetic trends for six traits in the Duroc population from 2012 to 2021, we observed significant genetic trends for AGE, BL, and BH. Particularly noteworthy is the rapid decline in the genetic trend for AGE, indicating an enhancement in the pig's growth rate through selective breeding. Therefore, we believe that some challenging-to-select traits can benefit from the genetic correlations between traits. By selecting easily measurable traits, they can gain from synergistic selection effects, leading to genetic progress. Conducting population genetic parameter analysis can assist us in devising breeding strategies.

Landmarks in the history of selective sweeps

Half a century ago, a seminal article on the hitchhiking effect by Smith and Haigh inaugurated the concept of the selection signature. Selective sweeps are characterised by the rapid spread of an advantageous genetic variant through a population and hence play an important role in shaping evolution and research on genetic diversity. The process by which a beneficial allele arises and becomes fixed in a population, leading to a increase in the frequency of other linked alleles, is known as genetic hitchhiking or genetic draft. Kimura's neutral theory and hitchhiking theory are complementary, with Kimura's neutral evolution as the 'null model' and positive selection as the 'signal'. Both are widely accepted in evolution, especially with genomics enabling precise measurements. Significant advances in genomic technologies, such as next-generation sequencing, high-density SNP arrays and powerful bioinformatics tools, have made it possible to systematically investigate selection signatures in a variety of species. Although the history of selection signatures is relatively recent, progress has been made in the last two decades, owing to the increasing availability of large-scale genomic data and the development of computational methods. In this review, we embark on a journey through the history of research on selective sweeps, ranging from early theoretical work to recent empirical studies that utilise genomic data.

Identification of Selection Signatures and Loci Associated with Important Economic Traits in Yunan Black and Huainan Pigs

Henan Province is located in central China and rich in domestic pig populations; Huainan (HN) pigs are one of three Henan indigenous breeds with great performance, including early maturation, strong disease resistance and high meat quality. Yunan (YN) black pigs are a typical, newly cultivated breed, synthesized between HN pigs and American Duroc, and are subjected to selection for important traits, such as fast growth and excellent meat quality. However, the genomic differences, selection signatures and loci associated with important economic traits in YN black pigs and HN pigs are still not well understood. In this study, based on high-density SNP chip analysis of 159 samples covering commercial DLY (Duroc × Landrace × Large White) pigs, HN pigs and YN black pigs, we performed a comprehensive analysis of phylogenetic relationships and genetic diversity among the three breeds. Furthermore, we used composite likelihood ratio tests (CLR) and F-statistics (Fst) to identify specific signatures of selection associated with important economic traits and potential candidate genes. We found 147 selected regions (top 1%) harboring 90 genes based on genetic differentiation (Fst) in the YN-DLY group. In the HN-DLY group, 169 selected regions harbored 58 genes. In the YN-HN group, 179 selected regions harbored 77 genes. In addition, the QTLs database with the most overlapping regions was associated with triglyceride level, number of mummified pigs, hemoglobin and loin muscle depth for YN black pigs, litter size and intramuscular fat content for HN pigs, and humerus length, linolenic acid content and feed conversion ratio mainly in DLY pigs. Of note, overlapping 14 tissue-specific promoters’ annotation with the top Fst 1% selective regions systematically demonstrated the muscle-specific and hypothalamus-specific regulatory elements in YN black pigs. Taken together, these results contribute to an accurate knowledge of crossbreeding, thus benefitting the evaluation of production performance and improving the genome-assisted breeding of other important indigenous pig in the future.

Identification of Homozygous Regions With Adverse Effects on the Five Economic Traits of Duroc Pigs

Runs of homozygosity (ROH) are widely used to estimate genomic inbreeding, which is linked to inbreeding depression on phenotypes. However, the adverse effects of specific homozygous regions on phenotypic characteristics are rarely studied in livestock. In this study, the 50 K SNP data of 3,770 S21 Duroc (American origin) and 2,096 S22 Duroc (Canadian origin) pigs were used to investigate the harmful ROH regions on five economic traits. The results showed that the two Duroc lines had different numbers and distributions of unfavorable ROHs, which may be related to the different selection directions and intensities between the two lines. A total of 114 and 58 ROH segments were found with significant adverse effects on the economic traits of S21 and S22 pigs, respectively. Serval pleiotropic ROHs were detected to reduce two or multiple phenotypic performances in two Duroc populations. Candidate genes in these shared regions were mainly related to growth, fertility, immunity, and fat deposition. We also observed that some ROH genotypes may cause opposite effects on different traits. This study not only enhances our understanding of the adverse effects of ROH on phenotypes, but also indicates that ROH information could be incorporated into breeding programs to estimate and control the detrimental effects of homozygous regions.

Population Genetic Structure and Selection Signature Analysis of Beijing Black Pig

Beijing Black pig is an excellent cultivated black pig breed in China, with desirable body shape, tender meat quality and robust disease resistance. To explore the level of admixture and selection signatures of Beijing Black pigs, a total number of 90 individuals covering nine pig breeds were used in our study, including Beijing Black pig, Large White, Landrace, Duroc, Lantang pig, Luchuan pig, Mashen pig, Huainan pig and Min pig. These animals were resequenced with 18.19 folds mapped read depth on average. Generally, we found that Beijing Black pig was genetically closer to commercial pig breeds by population genetic structure and genetic diversity analysis, and was also affected by Chinese domestic breeds Huainan pig and Min pig. These results are consistent with the cross-breeding history of Beijing Black pig. Selection signal detections were performed on three pig breeds, Beijing Black pig, Duroc and Large White, using three complementary methods (F_ST, θπ, and XP-EHH). In total, 1,167 significant selected regions and 392 candidate genes were identified. Functional annotations were enriched to pathways related to immune processes and meat and lipid metabolism. Finally, potential candidate genes, influencing meat quality (GPHA2, EHD1, HNF1A, C12orf43, GLTP, TRPV4, MVK, and MMAB), reproduction (PPP2R5B and MAP9), and disease resistance (OASL, ANKRD13A, and GIT2), were further detected by gene annotation analysis. Our results advanced the understanding of the genetic mechanism behind artificial selection of Beijing Black pigs, and provided theoretical basis for the subsequent breeding and genetic research of this breed.

Runs of Homozygosity Uncover Potential Functional-Altering Mutation Associated With Body Weight and Length in Two Duroc Pig Lines

Runs of homozygosity (ROH) are widely used to investigate genetic diversity, demographic history, and positive selection signatures of livestock. Commercial breeds provide excellent materials to reveal the landscape of ROH shaped during the intense selection process. Here, we used the GeneSeek Porcine 50K single-nucleotide polymorphism (SNP) Chip data of 3,770 American Duroc (AD) and 2,096 Canadian Duroc (CD) pigs to analyze the genome-wide ROH. First, we showed that AD had a moderate genetic differentiation with CD pigs, and AD had more abundant genetic diversity and significantly lower level of inbreeding than CD pigs. In addition, sows had larger levels of homozygosity than boars in AD pigs. These differences may be caused by differences in the selective intensity. Next, ROH hotspots revealed that many candidate genes are putatively under selection for growth, sperm, and muscle development in two lines. Population-specific ROHs inferred that AD pigs may have a special selection for female reproduction, while CD pigs may have a special selection for immunity. Moreover, in the overlapping ROH hotspots of two Duroc populations, we observed a missense mutation (rs81216249) located in the growth and fat deposition-related supergene (ARSB-DMGDH-BHMT) region. The derived allele of this variant originated from European pigs and was nearly fixed in Duroc pigs. Further selective sweep and association analyses indicated that this supergene was subjected to strong selection and probably contributed to the improvement of body weight and length in Duroc pigs. These findings will enhance our understanding of ROH patterns in different Duroc lines and provide promising trait-related genes and a functional-altering marker that can be used for genetic improvement of pigs.

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.

Detection of selection signatures for response to Aleutian mink disease virus infection in American mink

Aleutian disease (AD) is the most significant health issue for farmed American mink. The objective of this study was to identify the genomic regions subjected to selection for response to infection with Aleutian mink disease virus (AMDV) in American mink using genotyping by sequencing (GBS) data. A total of 225 black mink were inoculated with AMDV and genotyped using a GBS assay based on the sequencing of ApeKI-digested libraries. Five AD-characterized phenotypes were used to assign animals to pairwise groups. Signatures of selection were detected using integrated measurement of fixation index (F_ST) and nucleotide diversity (θπ), that were validated by haplotype-based (hap-FLK) test. The total of 99 putatively selected regions harbouring 63 genes were detected in different groups. The gene ontology revealed numerous genes related to immune response (e.g. TRAF3IP2, WDR7, SWAP70, CBFB, and GPR65), liver development (e.g. SULF2, SRSF5) and reproduction process (e.g. FBXO5, CatSperβ, CATSPER4, and IGF2R). The hapFLK test supported two strongly selected regions that contained five candidate genes related to immune response, virus–host interaction, reproduction and liver regeneration. This study provided the first map of putative selection signals of response to AMDV infection in American mink, bringing new insights into genomic regions controlling the AD phenotypes.

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.

Genome-wide analysis of expression QTL (eQTL) and allele-specific expression (ASE) in pig muscle identifies candidate genes for meat quality traits

BACKGROUND

Genetic analysis of gene expression level is a promising approach for characterizing candidate genes that are involved in complex economic traits such as meat quality. In the present study, we conducted expression quantitative trait loci (eQTL) and allele-specific expression (ASE) analyses based on RNA-sequencing (RNAseq) data from the longissimus muscle of 189 Duroc × Luchuan crossed pigs in order to identify some candidate genes for meat quality traits.

RESULTS

Using a genome-wide association study based on a mixed linear model, we identified 7192 cis-eQTL corresponding to 2098 cis-genes (p ≤ 1.33e-3, FDR ≤ 0.05) and 6400 trans-eQTL corresponding to 863 trans-genes (p ≤ 1.13e-6, FDR ≤ 0.05). ASE analysis using RNAseq SNPs identified 9815 significant ASE-SNPs in 2253 unique genes. Integrative analysis between the cis-eQTL and ASE target genes identified 540 common genes, including 33 genes with expression levels that were correlated with at least one meat quality trait. Among these 540 common genes, 63 have been reported previously as candidate genes for meat quality traits, such as PHKG1 (q-value = 1.67e-6 for the leading SNP in the cis-eQTL analysis), NUDT7 (q-value = 5.67e-13), FADS2 (q-value = 8.44e-5), and DGAT2 (q-value = 1.24e-3).

CONCLUSIONS

The present study confirmed several previously published candidate genes and identified some novel candidate genes for meat quality traits via eQTL and ASE analyses, which will be useful to prioritize candidate genes in further studies.

Detection of Selection Signatures Underlying Production and Adaptive Traits Based on Whole-Genome Sequencing of Six Donkey Populations

Simple Summary

After a long period of artificial selection, the donkey now presents a variety of body types and production performance values. In this experiment, we performed selective signal scanning on the second-generation resequencing data of six different varieties. The regions and candidate genes related to artificial selection were identified to provide reference for future breeding.

Abstract

Donkeys (Equus asinus) are an important farm animal. After long-term natural and artificial selection, donkeys now exhibit a variety of body sizes and production performance values. In this study, six donkey breeds, representing different regions and phenotypes, were used for second-generation resequencing. The sequencing results revealed more than seven million single nucleotide variants (SNVs), with an average of more than four million SNVs per species. We combined two methods, Z-transformed heterozygosity (ZHp) and unbiased estimates of pairwise fixation index (di) values, to analyze the signatures of selection. We mapped 11 selected regions and identified genes associated with coat color, body size, motion capacity, and high-altitude adaptation. These candidate genes included staining (ASIP and KITLG), body type (ACSL4, BCOR, CDKL5, LCOR, NCAPG, and TBX3), exercise (GABPA), and adaptation to low-oxygen environments (GLDC and HBB). We also analyzed the SNVs of the breed-specific genes for their potential functions and found that there are three varieties in the conserved regions with breed-specific mutation sites. Our results provide data to support the establishment of the donkey SNV chip and reference information for the utilization of the genetic resources of Chinese domestic donkeys.

Genome-Wide Detection of Selection Signatures in Duroc Revealed Candidate Genes Relating to Growth and Meat Quality

With the development of high-throughput genotyping techniques, selection signatures in the genome of domestic pigs have been extensively interrogated in the last decade. The Duroc, a major commercial pig breed famous for its fast growth rate and high lean ratio, has not been extensively studied focusing on footprints of intensively artificial selection in their genomes by a lot of re-sequencing data. The goal of this study was to investigate genomic regions under artificial selection and their contribution to the unique phenotypic traits of the Duroc using whole-genome resequencing data from 97 pigs. Three complementary methods (d_i, CLR, and iHH12) were implemented for selection signature detection. In Total, 464 significant candidate regions were identified, which covered 46.4 Mb of the pig genome. Within the identified regions, 709 genes were annotated, including 600 candidate protein-coding genes (486 functionally annotated genes) and 109 lncRNA genes. Genes undergoing selective pressure were significantly enriched in the insulin resistance signaling pathway, which may partly explain the difference between the Duroc and other breeds in terms of growth rate. The selection signatures identified in the Duroc population demonstrated positive pressures on a set of important genes with potential functions that are involved in many biological processes. The results provide new insights into the genetic mechanisms of fast growth rate and high lean mass, and further facilitate follow-up studies on functional genes that contribute to the Duroc's excellent phenotypic traits.

Whole-genome sequencing of European autochthonous and commercial pig breeds allows the detection of signatures of selection for adaptation of genetic resources to different breeding and production systems

Background

Natural and artificial directional selection in cosmopolitan and autochthonous pig breeds and wild boars have shaped their genomes and resulted in a reservoir of animal genetic diversity. Signatures of selection are the result of these selection events that have contributed to the adaptation of breeds to different environments and production systems. In this study, we analysed the genome variability of 19 European autochthonous pig breeds (Alentejana, Bísara, Majorcan Black, Basque, Gascon, Apulo-Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano, Sarda, Krškopolje pig, Black Slavonian, Turopolje, Moravka, Swallow-Bellied Mangalitsa, Schwäbisch-Hällisches Schwein, Lithuanian indigenous wattle and Lithuanian White old type) from nine countries, three European commercial breeds (Italian Large White, Italian Landrace and Italian Duroc), and European wild boars, by mining whole-genome sequencing data obtained by using a DNA-pool sequencing approach. Signatures of selection were identified by using a single-breed approach with two statistics [within-breed pooled heterozygosity (H_P) and fixation index (F_ST)] and group-based F_ST approaches, which compare groups of breeds defined according to external traits and use/specialization/type.

Results

We detected more than 22 million single nucleotide polymorphisms (SNPs) across the 23 compared populations and identified 359 chromosome regions showing signatures of selection. These regions harbour genes that are already known or new genes that are under selection and relevant for the domestication process in this species, and that affect several morphological and physiological traits (e.g. coat colours and patterns, body size, number of vertebrae and teats, ear size and conformation, reproductive traits, growth and fat deposition traits). Wild boar related signatures of selection were detected across all the genome of several autochthonous breeds, which suggests that crossbreeding (accidental or deliberate) occurred with wild boars.

Conclusions

Our findings provide a catalogue of genetic variants of many European pig populations and identify genome regions that can explain, at least in part, the phenotypic diversity of these genetic resources.

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.

Identifying Selection Signatures for Backfat Thickness in Yorkshire Pigs Highlights New Regions Affecting Fat Metabolism

Identifying the genetic basis of improvement in pigs contributes to our understanding of the role of artificial selection in shaping the genome. Here we employed the Cross Population Extended Haplotype Homozogysity (XPEHH) and the Wright's fixation index (F_ST) methods to detect trait-specific selection signatures by making phenotypic gradient differential population pairs, and then attempted to map functional genes of six backfat thickness traits in Yorkshire pigs. The results indicate that a total of 283 and 466 single nucleotide polymorphisms (SNPs) were identified as trait-specific selection signatures using F_ST and XPEHH, respectively. Functional annotation suggested that the genes overlapping with the trait-specific selection signatures such as OSBPL8, ASAH2, SMCO2, GBE1, and ABL1 are responsible for the phenotypes including fat metabolism, lean body mass and fat deposition, and transport in mouse. Overall, the study developed the methods of gene mapping on the basis of identification of selection signatures. The candidate genes putatively associated with backfat thickness traits can provide important references and fundamental information for future pig-breeding programs.