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

High-throughput untargeted metabolomics reveals metabolites and metabolic pathways that differentiate two divergent pig breeds

Metabolomics can describe the molecular phenome and may contribute to dissecting the biological processes linked to economically relevant traits in livestock species. Comparative analyses of metabolomic profiles in purebred pigs can provide insights into the basic biological mechanisms that may explain differences in production performances. Following this concept, this study was designed to compare, on a large scale, the plasma metabolomic profiles of two Italian heavy pig breeds (Italian Duroc and Italian Large White) to indirectly evaluate the impact of their different genetic backgrounds on the breed metabolomes. We utilised a high-throughput untargeted metabolomics approach in a total of 962 pigs that allowed us to detect and relatively quantify 722 metabolites from various biological classes. The molecular data were analysed using a bioinformatics pipeline specifically designed for identifying differentially abundant metabolites between the two breeds in a robust and statistically significant manner, including the Boruta algorithm, which is a Random Forest wrapper, and sparse Partial Least Squares Discriminant Analysis (sPLS-DA) for feature selection. After thoroughly evaluating the impact of random components on missing value imputation, 100 discriminant metabolites were selected by Boruta and 17 discriminant metabolites (all included within the previous list) were identified with sPLS-DA. About half of the 100 discriminant metabolites had a higher concentration in one or the other breed (48 in Italian Large White pigs, with a prevalence of amino acids and peptides; 52 in Italian Duroc pigs, with a prevalence of lipids). These metabolites were from seven distinct super pathways and had an absolute mean value of percentage difference between the two breeds (|Δ|%) of 39.2 ± 32.4. Six of these metabolites had |Δ|%> 100. A general correlation network analysis based on Boruta-identified metabolites consisted of 31 singletons and 69 metabolites connected by 141 edges, with two large clusters (> 15 nodes), three medium clusters (3-6 nodes) and eight additional pairs, with most metabolites belonging to the same super pathway. The major cluster representing the lipids super-pathway included 24 metabolites, primarily sphingomyelins. Overall, this study identified metabolomic differences between Italian Duroc and Italian Large White pigs explained by the specific genetic background of the two breeds. These biomarkers can explain the biological differences between these two breeds and can have potential practical applications in pig breeding and husbandry.

The Evolution and Role of Molecular Tools in Measuring Diversity and Genomic Selection in Livestock Populations (Traditional and Up-to-Date Insights): A Comprehensive Exploration

Distinctive molecular approaches and tools, particularly high-throughput SNP genotyping, have been applied to determine and discover SNPs, potential genes of interest, indicators of evolutionary selection, genetic abnormalities, molecular indicators, and loci associated with quantitative traits (QTLs) in various livestock species. These methods have also been used to obtain whole-genome sequencing (WGS) data, enabling the implementation of genomic selection. Genomic selection allows for selection decisions based on genomic-estimated breeding values (GEBV). The estimation of GEBV relies on the calculation of SNP effects using prediction equations derived from a subset of individuals in the reference population who possess both SNP genotypes and phenotypes for target traits. Compared to traditional methods, modern genomic selection methods offer advantages for sex-limited traits, low heritability traits, late-measured traits, and the potential to increase genetic gain by reducing generation intervals. The current availability of high-density genotyping and next-generation sequencing data allow for genome-wide scans for selection. This investigation provides an overview of the essential role of advanced molecular tools in studying genetic diversity and implementing genomic selection. It also highlights the significance of adaptive selection in light of new high-throughput genomic technologies and the establishment of selective comparisons between different genomes. Moreover, this investigation presents candidate genes and QTLs associated with various traits in different livestock species, such as body conformation, meat production and quality, carcass characteristics and composition, milk yield and composition, fertility, fiber production and characteristics, and disease resistance.

Genetic basis of phenotypic convergence in pig terminal sires using pathway-based selection signature detection methods

The primary purpose of genetic improvement in lean pig breeds is to enhance production performance. Owing to their similar breeding directions, Duroc and Pietrain pigs are ideal models for investigating the phenotypic convergence underlying artificial selection. However, most important economic traits are controlled by a polygenic basis, so traditional strategies for detecting selection signatures may not fully reveal the genetic basis of complex traits. The pathway-based gene network analysis method utilizes each pathway as a unit, overcoming the limitations of traditional strategies for detecting selection signatures by revealing the selection of complex biological processes. Here, we utilized 13 122 398 high-quality SNPs from whole-genome sequencing data of 48 Pietrain pigs, 156 Duroc pigs and 36 European wild boars to detect selective signatures. After calculating FST and iHS scores, we integrated the pathway information and utilized the r/bioconductor graphite and signet packages to construct gene networks, identify subnets and uncover candidate genes underlying selection. Using the traditional strategy, a total of 47 genomic regions exhibiting parallel selection were identified. The enriched genes, including INO80, FZR1, LEPR and FAF1, may be associated with reproduction, fat deposition and skeletal development. Using the pathway-based selection signatures detection method, we identified two significant biological pathways and eight potential candidate genes underlying parallel selection, such as VTN, FN1 and ITGAV. This study presents a novel strategy for investigating the genetic basis of complex traits and elucidating the phenotypic convergence underlying artificial selection, by integrating traditional selection signature methods with pathway-based gene network analysis.

Expression genome-wide association study identifies key regulatory variants enriched with metabolic and immune functions in four porcine tissues

BACKGROUND

Integration of high throughput DNA genotyping and RNA-sequencing data enables the discovery of genomic regions that regulate gene expression, known as expression quantitative trait loci (eQTL). In pigs, efforts to date have been mainly focused on purebred lines for traits with commercial relevance as such growth and meat quality. However, little is known on genetic variants and mechanisms associated with the robustness of an animal, thus its overall health status. Here, the liver, lung, spleen, and muscle transcriptomes of 100 three-way crossbred female finishers were studied, with the aim of identifying novel eQTL regulatory regions and transcription factors (TFs) associated with regulation of porcine metabolism and health-related traits.

RESULTS

An expression genome-wide association study with 535,896 genotypes and the expression of 12,680 genes in liver, 13,310 genes in lung, 12,650 genes in spleen, and 12,595 genes in muscle resulted in 4,293, 10,630, 4,533, and 6,871 eQTL regions for each of these tissues, respectively. Although only a small fraction of the eQTLs were annotated as cis-eQTLs, these presented a higher number of polymorphisms per region and significantly stronger associations with their target gene compared to trans-eQTLs. Between 20 and 115 eQTL hotspots were identified across the four tissues. Interestingly, these were all enriched for immune-related biological processes. In spleen, two TFs were identified: ERF and ZNF45, with key roles in regulation of gene expression.

CONCLUSIONS

This study provides a comprehensive analysis with more than 26,000 eQTL regions identified that are now publicly available. The genomic regions and their variants were mostly associated with tissue-specific regulatory roles. However, some shared regions provide new insights into the complex regulation of genes and their interactions that are involved with important traits related to metabolism and immunity.

Genetic assessment and candidate genes identification for breed-specific characteristics of Qingyuan partridge chicken based on runs of homozygosity

BACKGROUND

Several core breeding and supporting lines of the Qingyuan partridge chicken, a representative local chicken breed in China, have been developed over 20 years. Consequently, its economic traits related to growth and reproduction have been significantly improved by breeding selection and commercial utilization, but some characteristic traits, such as partridge feathers, high meat quality and sufficient flavor, have always been retained. However, effective methods for genetic assessment and functional gene exploration of similar trait groups are lacking. The presence of identical haplotype fragments transmitted from parent to offspring results in runs of homozygosity (ROH), which offer an efficient solution. In this study, genomes of 134 Qingyuan partridge chickens representing two breeding populations and one preserved population were re-sequenced to evaluate the genetic diversity and explore functional genes by analyzing the diversity, distribution, and frequency of ROH.

RESULTS

The results showed a low level of genomic linkage and degree of inbreeding within both the bred and preserved populations, suggesting abundant genetic diversity and an adequate genetic potential of the Qingyuan partridge chicken. Throughout the long-term selection process, 21 genes, including GLI3, ANO5, BLVRA, EFNB2, SLC5A12, and SVIP, associated with breed-specific characteristics were accumulated within three ROH islands, whereas another 21 genes associated with growth traits including IRX1, IRX2, EGFR, TPK1, NOVA1, BDNF and so on were accumulated within five ROH islands.

CONCLUSIONS

These findings provide new insights into the genetic assessment and identification of genes with breed-specific and selective characteristics, offering a solid genetic basis for breeding and protection of Qingyuan partridge chickens.

Revelation of genetic diversity and genomic footprints of adaptation in Indian pig breeds

In the present study, the genetic diversity measures among four Indian domestic breeds of pig namely Agonda Goan, Ghurrah, Ghungroo, and Nicobari, of different agro-climatic regions of country were explored and compared with European commercial breeds, European wild boar and Chinese domestic breeds. The double digest restriction site-associated DNA sequencing (ddRADseq) data of Indian pigs (102) and Landrace (10 animals) were generated and whole genome sequencing data of exotic pigs (60 animals) from public data repository were used in the study. The principal component analysis (PCA), admixture analysis and phylogenetic analysis revealed that Indian breeds were closer in ancestry to Chinese breeds than European breeds. European breeds exhibited highest genetic diversity measures among all the considered breeds. Among Indian breeds, Agonda Goan and Ghurrah were found to be more genetically diverse than Nicobari and Ghungroo. The selection signature regions in Indian pigs were explored using iHS and XP-EHH, and during iHS analysis, it was observed that genes related to growth, reproduction, health, meat quality, sensory perception and behavior were found to be under selection pressure in Indian pig breeds. Strong selection signatures were recorded in 24.25-25.25 Mb region of SSC18, 123.25-124 Mb region of SSC15 and 118.75-119.5 Mb region of SSC2 in most of the Indian breeds upon pairwise comparison with European commercial breeds using XP-EHH. These regions were harboring some important genes such as EPHA4 for thermotolerance, TAS2R16, FEZF1, CADPS2 and PTPRZ1 for adaptability to scavenging system of rearing, TRIM36 and PGGT1B for disease resistance and CCDC112, PIAS1, FEM1B and ITGA11 for reproduction.

Genomic Scan for Runs of Homozygosity and Selective Signature Analysis to Identify Candidate Genes in Large White Pigs

Large White pigs are extensively utilized in China for their remarkable characteristics of rapid growth and the high proportion of lean meat. The economic traits of pigs, comprising reproductive and meat quality traits, play a vital role in swine production. In this study, 2295 individuals, representing three different genetic backgrounds Large White pig populations were used: 500 from the Canadian line, 295 from the Danish line, and 1500 from the American line. The GeneSeek 50K GGP porcine HD array was employed to genotype the three pig populations. Firstly, genomic selective signature regions were identified using the pairwise fixation index (FST) and locus-specific branch length (LSBL). By applying a top 1% threshold for both parameters, a total of 888 candidate selective windows were identified, harbouring 1571 genes. Secondly, the investigation of regions of homozygosity (ROH) was performed utilizing the PLINK software. In total, 25 genomic regions exhibiting a high frequency of ROHs were detected, leading to the identification of 1216 genes. Finally, the identified potential functional genes from candidate genomic regions were annotated, and several important candidate genes associated with reproductive traits (ADCYAP1, U2, U6, CETN1, Thoc1, Usp14, GREB1L, FGF12) and meat quality traits (MiR-133, PLEKHO1, LPIN2, SHANK2, FLVCR1, MYL4, SFRP1, miR-486, MYH3, STYX) were identified. The findings of this study provide valuable insights into the genetic basis of economic traits in Large White pigs and may have potential use in future pig breeding programs.

Genomic Association Analysis of Growth and Backfat Traits in Large White Pigs

The pig industry is significantly influenced by complex traits such as growth rate and fat deposition, which have substantial implications for economic returns. Over the years, remarkable genetic advancements have been achieved through intense artificial selection to enhance these traits in pigs. In this study, we aimed to investigate the genetic factors that contribute to growth efficiency and lean meat percentages in Large White pigs. Specifically, we focused on analyzing two key traits: age at 100 kg live weight (AGE100) and backfat thickness at 100 kg (BF100), in three distinct Large White pig populations-500 Canadian, 295 Danish, and 1500 American Large White pigs. By employing population genomic techniques, we observed significant population stratification among these pig populations. Utilizing imputed whole-genome sequencing data, we conducted single population genome-wide association studies (GWAS) as well as a combined meta-analysis across the three populations to identify genetic markers associated with the aforementioned traits. Our analyses highlighted several candidate genes, such as CNTN1-which has been linked to weight loss in mice and is potentially influential for AGE100-and MC4R, which is associated with obesity and appetite and may impact both traits. Additionally, we identified other genes-namely, PDZRN4, LIPM, and ANKRD22-which play a partial role in fat growth. Our findings provide valuable insights into the genetic basis of these important traits in Large White pigs, which may inform breeding strategies for improved production efficiency and meat quality.

Integrating genome-wide association studies and population genomics analysis reveals the genetic architecture of growth and backfat traits in pigs

Growth and fat deposition are complex traits, which can affect economical income in the pig industry. Due to the intensive artificial selection, a significant genetic improvement has been observed for growth and fat deposition in pigs. Here, we first investigated genomic-wide association studies (GWAS) and population genomics (e.g., selection signature) to explore the genetic basis of such complex traits in two Large White pig lines (n = 3,727) with the GeneSeek GGP Porcine HD array (n = 50,915 SNPs). Ten genetic variants were identified to be associated with growth and fatness traits in two Large White pig lines from different genetic backgrounds by performing both within-population GWAS and cross-population GWAS analyses. These ten significant loci represented eight candidate genes, i.e., NRG4, BATF3, IRS2, ANO1, ANO9, RNF152, KCNQ5, and EYA2. One of them, ANO1 gene was simultaneously identified for both two lines in BF100 trait. Compared to single-population GWAS, cross-population GWAS was less effective for identifying SNPs with population-specific effect, but more powerful for detecting SNPs with population-shared effects. We further detected genomic regions specifically selected in each of two populations, but did not observe a significant enrichment for the heritability of growth and backfat traits in such regions. In summary, the candidate genes will provide an insight into the understanding of the genetic architecture of growth-related traits and backfat thickness, and may have a potential use in the genomic breeding programs in pigs.

Assessing Genetic Diversity and Searching for Selection Signatures by Comparison between the Indigenous Livni and Duroc Breeds in Local Livestock of the Central Region of Russia

Indigenous pig breeds are mainly associated with the adaptive capacity that is necessary to respond adequately to climate change, food security, and livelihood needs, and natural resources conservation. Livni pigs are an indigenous fat-type breed farmed in a single farm in the Orel region and located in the Central European part of the Russian Federation. To determine the genomic regions and genes that are affected by artificial selection, we conducted the comparative study of two pig breeds with different breeding histories and breeding objectives, i.e., the native fat-type Livni and meat-type Duroc breeds using the Porcine GGP HD BeadChip, which contains ~80,000 SNPs. To check the Livni pigs for possible admixture, the Landrace and the Large White breeds were included into the study of genetic diversity as these breeds participated in the formation of the Livni pigs. We observed the highest level of genetic diversity in Livni pigs compared to commercial breeds (UHE = 0.409 vs. 0.319–0.359, p < 0.001; AR = 1.995 vs. 1.894–1.964, p < 0.001). A slight excess of heterozygotes was found in all of the breeds. We identified 291 candidate genes, which were localized within the regions under putative selection, including 22 and 228 genes, which were specific for Livni and Duroc breeds, respectively, and 41 genes common for both breeds. A detailed analysis of the molecular functions identified the genes, which were related to the formation of meat and fat traits, and adaptation to environmental stress, including extreme temperatures, which were different between breeds. Our research results are useful for conservation and sustainable breeding of Livni breed, which shows a high level of genetic diversity. This makes Livni one of the valuable national pig genetic resources.

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.

Genomic Diversity Profiling and Breed-Specific Evolutionary Signatures of Selection in Arunachali Yak

Arunachali yak, the only registered yak breed of India, is crucial for the economic sustainability of pastoralist Monpa community. This study intended to determine the genomic diversity and to identify signatures of selection in the breed. Previously available double digest restriction-site associated DNA (ddRAD) sequencing data of Arunachali yak animals was processed and 99,919 SNPs were considered for further analysis. The genomic diversity profiled based on nucleotide diversity, π (π = 0.041 in 200 bp windows), effective population size, Ne (Ne = 83) and Runs of homozygosity (ROH) (predominance of shorter length ROHs) was found to be optimum. Subsequently, 207 regions were identified to be under selective sweeps through de-correlated composite of multiple signals (DCMS) statistic which combined three individual test statistics viz. π, Tajima's D and |iHS| in non-overlapping 100 kb windows. Mapping of these regions revealed 611 protein-coding genes including KIT, KITLG, CDH12, FGG, FGA, FGB, PDGFRA, PEAR1, STXBP3, olfactory receptor genes (OR5K3, OR5H6 and OR1E1) and taste receptor genes (TAS2R1, TAS2R3 and TAS2R4). Functional annotation highlighted that biological processes like platelet aggregation and sensory perception were the most overrepresented and the associated regions could be considered as breed-specific signatures of selection in Arunachali yak. These findings point towards evolutionary role of natural selection in environmental adaptation of Arunachali yak population and provide useful insights for pursuing genome-wide association studies in future.

Assessment of Heterozygosity and Genome-Wide Analysis of Heterozygosity Regions in Two Duroc Pig Populations

Heterozygosity can effectively reflect the diverse models of population structure and demographic history. However, the genomic distribution of heterozygotes and the correlation between regions of heterozygosity (runs of heterozygosity, ROHet) and phenotypes are largely understudied in livestock. The objective of this study was to identify ROHet in the Duroc pig genome, and investigate the relationships between ROHet and eight important economic traits. Here, we genotyped 3,770 American Duroc (S21) and 2,096 Canadian Duroc (S22) pigs using 50 K single nucleotide polymorphism array to analyze heterozygosity. A total of 145,010 and 84,396 ROHets were characterized for S21 and S22 populations, respectively. ROHet segments were mostly enriched in 1–2 Mb length classification (75.48% in S21 and 72.25% in S22). The average genome length covered by ROHet was 66.53 ± 12.20 Mb in S21 and 73.32 ± 13.77 Mb in S22 pigs. Additionally, we detected 20 and 13 ROHet islands in S21 and S22 pigs. Genes in these genomic regions were mainly involved in the biological processes of immunity and reproduction. Finally, the genome-wide ROHet-phenotypes association analysis revealed that 130 ROHets of S21 and 84 ROHets of S22 were significantly associated with eight economic traits. Among the candidate genes in the significant ROHet regions, 16 genes related to growth, metabolism, and meat quality were considered as candidate genes for important economic traits of pigs. This work preliminarily explores the effect of heterozygosity-rich regions in the pig genome on production performance and provides new insights for subsequent research on pig genetic improvement.

Trait-specific Selection Signature Detection Reveals Novel Loci of Meat Quality in Large White Pigs

In past decades, meat quality traits have been shaped by human-driven selection in the process of genetic improvement programs. Exploring the potential genetic basis of artificial selection and mapping functional candidate genes for economic traits are of great significance in genetic improvement of pigs. In this study, we focus on investigating the genetic basis of five meat quality traits, including intramuscular fat content (IMF), drip loss, water binding capacity, pH at 45 min (pH45min), and ultimate pH (pH24h). Through making phenotypic gradient differential population pairs, Wright’s fixation index (F_ST) and the cross-population extended haplotype homozogysity (XPEHH) were applied to detect selection signatures for these five traits. Finally, a total of 427 and 307 trait-specific selection signatures were revealed by F_ST and XPEHH, respectively. Further bioinformatics analysis indicates that some genes, such as USF1, NDUFS2, PIGM, IGSF8, CASQ1, and ACBD6, overlapping with the trait-specific selection signatures are responsible for the phenotypes including fat metabolism and muscle development. Among them, a series of promising trait-specific selection signatures that were detected in the high IMF subpopulation are located in the region of 93544042-95179724bp on SSC4, and the genes harboring in this region are all related to lipids and muscle development. Overall, these candidate genes of meat quality traits identified in this analysis may provide some fundamental information for further exploring the genetic basis of this complex trait.

Genomic Analysis Reveals Human-Mediated Introgression From European Commercial Pigs to Henan Indigenous Pigs

Introgression of genetic features from European pigs into Chinese pigs was reported possibly contributing to improvements in productivity traits, such as feed conversion efficiency and body size. However, the genomic differences from European pigs and the potential role of introgression in Henan indigenous pigs remains unclear. In this study, we found significant introgression from European pigs into the genome of Chinese indigenous pigs, especially in Henan indigenous pigs. The introgression in Henan indigenous pigs, particularly in the Nanyang black pig, was mainly derived from Duroc pigs. Most importantly, we found that the NR6A1, GPD2, and CSRNP3 genes were introgressed and reshaped by artificial selection, and these may have contributed to increases in pig body size and feed conversion efficiency. Our results suggest that human-mediated introgression and selection have reshaped the genome of Henan pigs and improved several of their desired traits. These findings contribute to our understanding of the history of Henan indigenous pigs and provide insights into the genetic mechanisms affecting economically important traits in pig populations.

Detection of selection signatures in farmed coho salmon (Oncorhynchus kisutch) using dense genome-wide information

Animal domestication and artificial selection give rise to gradual changes at the genomic level in populations. Subsequent footprints of selection, known as selection signatures or selective sweeps, have been traced in the genomes of many animal livestock species by exploiting variation in linkage disequilibrium patterns and/or reduction of genetic diversity. Domestication of most aquatic species is recent in comparison with land animals, and salmonids are one of the most important fish species in aquaculture. Coho salmon (Oncorhynchus kisutch), cultivated primarily in Chile, has been subjected to breeding programs to improve growth, disease resistance traits, and flesh color. This study aimed to identify selection signatures that may be involved in adaptation to culture conditions and traits of productive interest. To do so, individuals of two domestic populations cultured in Chile were genotyped with 200 thousand SNPs, and analyses were conducted using iHS, XP-EHH and CLR. Several signatures of selection on different chromosomal regions were detected across both populations. Some of the identified regions under selection contained genes such anapc2, alad, chp2 and myn, which have been previously associated with body weight in Atlantic salmon, or sec24d and robo1, which have been associated with resistance to Piscirickettsia salmonis in coho salmon. Findings in our study can contribute to an integrated genome-wide map of selection signatures, to help identify the genetic mechanisms of phenotypic diversity in coho salmon.