A genome-wide scan for signatures of selection in Chinese indigenous and commercial pig breeds

Integrating spatial transcriptomics and single-nucleus RNA-seq revealed the specific inhibitory effects of TGF-β on intramuscular fat deposition

Intramuscular fat (IMF) is a complex adipose tissue within skeletal muscle, appearing specially tissue heterogeneous, and the factors influencing its formation remain unclear. In conditions such as diabetes, aging, and muscle wasting, IMF was deposited in abnormal locations in skeletal muscle, damaged the normal physiological functions of skeletal muscle. Here, we used Longissimus dorsi muscles from pigs with different IMF contents as samples and adopted a method combining spatial transcriptome (ST) and single-nucleus RNA-seq to identify the spatial heterogeneity of IMF. ST revealed that genes involved in TGF-β signaling pathways were specifically highly enriched in IMF. In lean pigs, IMF autocrine produces more TGF-β2, while in obese pigs, IMF received more endothelial-derived TGF-β1. In vitro experiments have proven that porcine endothelial cells in a simulated high-fat environment released more TGF-β1 than TGF-β2. Moreover, under obesity mice, the addition of TGF-β after muscle injury abolished IMF production and slowed muscle repair, whereas TGF-β inhibition accelerated muscle repair. Our findings demonstrate that the TGF-β pathway specifically regulates these processes, suggesting it as a potential therapeutic target for managing muscle atrophy in obese patients and enhancing muscle repair while reducing IMF deposition.

A meta-analysis of genetic and phenotypic diversity of European local pig breeds reveals genomic regions associated with breed differentiation for production traits

BACKGROUND

Intense selection of modern pig breeds has resulted in genetic improvement of production traits while the performance of local pig breeds has remained lower. As local pig breeds have been bred in extensive systems, they have adapted to specific environmental conditions, resulting in a rich genotypic and phenotypic diversity. This study is based on European local pig breeds that have been genetically characterized using DNA-pool sequencing data and phenotypically characterized using breed level phenotypes related to stature, fatness, growth, and reproductive performance traits. These data were analyzed using a dedicated approach to detect signatures of selection linked to phenotypic traits in order to uncover potential candidate genes that may underlie adaptation to specific environments.

RESULTS

Analysis of the genetic data of European pig breeds revealed four main axes of genetic variation represented by the Iberian and three modern breeds (i.e. Large White, Landrace, and Duroc). In addition, breeds clustered according to their geographical origin, for example French Gascon and Basque breeds, Italian Apulo Calabrese and Casertana breeds, Spanish Iberian, and Portuguese Alentejano breeds. Principal component analysis of the phenotypic data distinguished the larger and leaner breeds with better growth potential and reproductive performance from the smaller and fatter breeds with low growth and reproductive efficiency. Linking the signatures of selection with phenotype identified 16 significant genomic regions associated with stature, 24 with fatness, 2 with growth, and 192 with reproduction. Among them, several regions contained candidate genes with possible biological effects on stature, fatness, growth, and reproductive performance traits. For example, strong associations were found for stature in two regions containing, respectively, the ANXA4 and ANTXR1 genes, for fatness in a region containing the DNMT3A and POMC genes and for reproductive performance in a region containing the HSD17B7 gene.

CONCLUSIONS

In this study on European local pig breeds, we used a dedicated approach for detecting signatures of selection that were supported by phenotypic data at the breed level to identify potential candidate genes that may have adapted to different living environments and production systems.

Integrative Analysis of Transcriptomic and Lipidomic Profiles Reveals a Differential Subcutaneous Adipose Tissue Mechanism among Ningxiang Pig and Berkshires, and Their Offspring

Adipose tissue composition contributes greatly to the quality and nutritional value of meat. Transcriptomic and lipidomic techniques were used to investigate the molecular mechanisms of the differences in fat deposition in Ningxiang pigs, Berkshires and F1 offspring. Transcriptomic analysis identified 680, 592, and 380 DEGs in comparisons of Ningxiang pigs vs. Berkshires, Berkshires vs. F1 offspring, and Ningxiang pigs vs. F1 offspring. The lipidomic analysis screened 423, 252, and 50 SCLs in comparisons of Ningxiang pigs vs. Berkshires, Berkshires vs. F1 offspring, and Ningxiang pigs vs. F1 offspring. Lycine, serine, and the threonine metabolism pathway, fatty acid biosynthesis and metabolism-related pathways were significantly enriched in comparisons of Berkshires vs. Ningxiang pigs and Berkshires vs. F1 offspring. The DEGs (PHGDH, LOC110256000) and the SCLs (phosphatidylserines) may have a great impact on the glycine, serine, and the threonine metabolism pathway. Moreover, the DEGs (FASN, ACACA, CBR4, SCD, ELOV6, HACD2, CYP3A46, CYP2B22, GPX1, and GPX3) and the SCLs (palmitoleic acid, linoleic acid, arachidonic acid, and icosadienoic acid) play important roles in the fatty acid biosynthesis and metabolism of fatty acids. Thus, the difference in fat deposition among Ningxiang pig, Berkshires, and F1 offspring may be caused by differences in the expression patterns of key genes in multiple enriched KEGG pathways. This research revealed multiple lipids that are potentially available biological indicators and screened key genes that are potential targets for molecular design breeding. The research also explored the molecular mechanisms of the difference in fat deposition among Ningxiang pig, Berkshires, and F1 pigs, and provided an insight into selection for backfat thickness and the fat composition of adipose tissue for future breeding strategies.

Analysis of polygenic selection in purebred and crossbred pig genomes using generation proxy selection mapping

BACKGROUND

Artificial selection on quantitative traits using breeding values and selection indices in commercial livestock breeding populations causes changes in allele frequency over time at hundreds or thousands of causal loci and the surrounding genomic regions. In population genetics, this type of selection is called polygenic selection. Researchers and managers of pig breeding programs are motivated to understand the genetic basis of phenotypic diversity across genetic lines, breeds, and populations using selection mapping analyses. Here, we applied generation proxy selection mapping (GPSM), a genome-wide association analysis of single nucleotide polymorphism (SNP) genotypes (38,294-46,458 markers) of birth date, in four pig populations (15,457, 15,772, 16,595 and 8447 pigs per population) to identify loci responding to artificial selection over a period of five to ten years. Gene-drop simulation analyses were conducted to provide context for the GPSM results. Selected loci within and across each population of pigs were compared in the context of swine breeding objectives.

RESULTS

The GPSM identified 49 to 854 loci as under selection (Q-values less than 0.10) across 15 subsets of pigs based on combinations of populations. The number of significant associations increased when data were pooled across populations. In addition, several significant associations were identified in more than one population. These results indicate concurrent selection objectives, similar genetic architectures, and shared causal variants responding to selection across these pig populations. Negligible error rates (less than or equal to 0.02%) of false-positive associations were found when testing GPSM on gene-drop simulated genotypes, suggesting that GPSM distinguishes selection from random genetic drift in actual pig populations.

CONCLUSIONS

This work confirms the efficacy and the negligible error rates of the GPSM method in detecting selected loci in commercial pig populations. Our results suggest shared selection objectives and genetic architectures across swine populations. The identified polygenic selection highlights loci that are important to swine production.

Replicative Study in Performance-Related Genes of Brazilian Elite Soccer Players Highlights Genetic Differences from African Ancestry and Similarities between Professional and U20 Youth Athletes

Classically, genetic association studies have attempted to assess genetic polymorphisms related to human physiology and physical performance. However, the heterogeneity of some findings drives the research to replicate, validate, and confirmation as essential aspects for ensuring their applicability in sports sciences. Genetic distance matrix and molecular variance analyses may offer an alternative approach to comparing athletes' genomes with those from public databases. Thus, we performed a complete sequencing of 44 genomes from male Brazilian first-division soccer players under 20 years of age (U20_BFDSC). The performance-related SNP genotypes were obtained from players and from the "1000 Genomes" database (European, African, American, East Asian, and South Asian). Surprisingly, U20_BFDSC performance-related genotypes had significantly larger F_ST levels (p < 0.00001) than African populations, although studies using ancestry markers have shown an important similarity between Brazilian and African populations (12-24%). U20_BFDSC were genetically similar to professional athletes, showing the intense genetic selection pressure likely to occur before this maturation stage. Our study highlighted that performance-related genes might undergo selective pressure due to physical performance and environmental, cognitive, and sociocultural factors. This replicative study suggests that molecular variance and Wright's statistics can yield novel conclusions in exercise science.

Small Intestinal Digestive Functions and Feed Efficiency Differ in Different Pig Breeds

Small intestinal growth and health affect its digestion and absorption ability, while little information exists about the small intestinal morphology and function differences among the different pig breeds. Therefore, 90 healthy 35 days of age Taoyuan black (TB), Xiangcun black (XB), and Duroc (DR) pigs (30 pigs per breed) with similar body weight (BW) of the same breed were reared to 185 days of age to evaluate the potential relationship between feed efficiency and small intestinal morphology and function at 80, 125, and 185 days of age. The results show that the TB and XB pigs had lower initial and final BW, ADG, and ADFI and plasma CHO and LDL-C levels, whereas they had higher plasma LIP levels and jejunal trypsin, invertase, lactase, and maltase activities and higher DM, ADF, Tyr, Arg, and His digestibility at 80 days of age compared with the DR pigs. At 125 days of age, TB and XB pigs had lower apparent total tract digestibility and plasma CHO, HDL-C, LDL-C, and NH3 levels; XB pigs had lower DM and NDF digestibility, and TB pigs had higher jejunal lactase and maltase activities. At 185 days of age, TB and XB pigs had lower DM, EE, ADF, and GE digestibility, while having higher plasma ALT and UN levels; TB pigs had higher plasma AST level and jejunal chymase activity. Furthermore, the plasma free amino acid contents, small intestinal VH, and nutrient transporter expression levels differed at different ages. Therefore, the different pig breeds exhibited significantly different growth performance and small intestinal growth, mainly resulting from the differences in digestive enzymes and nutrient transporters in the small intestine.

Screening of candidate genes related to differences in growth and development between Chinese indigenous and Western pig breeds

Neijiang (NJ) and Yacha (YC) are two indigenous pig breeds in the Sichuan basin of China, displaying higher resistance to diseases, lower lean ratio, and slower growth rate than the commercial Western pig breed Yorkshire (YS). The molecular mechanisms underlying the differences in growth and development between these pig breeds are still unknown. In the present study, five pigs from NJ, YC, and YS breeds were subjected to the whole genome resequencing, and then the differential single-nucleotide polymorphisms (SNPs) were screened using a 10-kb window sliding in 1-kb step using the F_st method. Finally, 48,924, 48,543, and 46,228 nonsynonymous single-nucleotide polymorphism loci (nsSNPs) were identified between NJ and YS, NJ and YC, and YC and YS, which highly or moderately affected 2,490, 800, and 444 genes, respectively. Moreover, three nsSNPs were detected in the genes of acetyl-CoA acetyltransferase 1 (ACAT1) insulin-like growth factor 2 receptor (IGF2R), insulin-like growth factor 2 and mRNA-binding protein 3 (IGF2BP3), which potentially affected the transformation of acetyl-CoA to acetoacetyl-CoA and the normal functions of the insulin signaling pathways. Moreover, serous determinations revealed significantly lower acetyl-CoA content in YC than in YS, supporting that ACAT1 might be a reason explaining the differences in growth and development between YC and YS breeds. Contents of phosphatidylcholine (PC) and phosphatidic acid (PA) significantly differed between the pig breeds, suggesting that glycerophospholipid metabolism might be another reason for the differences between Chinese and Western pig breeds. Overall, these results might contribute basic information to understand the genetic differences determining the phenotypical traits in pigs.

A Poly(dA:dT) Tract in the IGF1 Gene Is a Genetic Marker for Growth Traits in Pigs

Insulin-like growth factor 1 (IGF1) is an important regulator of body growth, development, and metabolism. The poly(dA:dT) tract affects the accessibility of transcription factor binding sites to regulate transcription. Therefore, this study assessed the effects of two poly(dA:dT) tracts on the transcriptional activity of porcine IGF1. The luciferase assay results demonstrated that the poly(dA:dT) tract 2 (−264/−255) was a positive regulatory element for IGF1 gene expression, and the activities between the different lengths of the poly(dA:dT) tract 2 were significant (p<0.01). The transcription factor C/EBPα inhibited the transcription of IGF1 by binding to tract 2, and the expression levels between the lengths of tract 2 after C/EBPα binding were also statistically different (p<0.01). Only the alleles 10T and 11T were found in the tract 2 in commercial pig breeds, while the 9T, 10T, and 11T alleles were found in Chinese native pig breeds. The allele frequencies were in Hardy−Weinberg equilibrium in all pig breeds. The genotypes of tract 2 were significantly associated with the growth traits (days to 115 kg and average daily gain) (p<0.05) in commercial pig breeds. Based on these findings, it can be concluded that the tract 2 mutation could be applied as a candidate genetic marker for growth trait selection in pig breeding programs.

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.

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.

Comparative transcriptome analysis of winter yaks in plateau and plain

The yak is an important source for the people living and ecological environment in the Qinghai-Tibet Plateau. In every winter, many domestic yaks will lose bodyweight or dead under cold and food scarcity. Moving the plateau yaks to farm in the plain is a useful approach to reduce their environmental stress and gain more production. In this study, we measured growth, slaughter and beef quality traits every month and sequenced mRNA expression levels of muscles of two groups yaks living in plateau and plain respectively. We found there is significant difference (p-value <0.01) in the third (60 days), fourth (90 days), fifth (120 days) and sixth (150 days) weights between subpopulations in the plateau and plain. We identified 540 different expressed genes (DEGs) including 123 known genes and 417 unknown genes. Using the weighted correlation network analysis (WGCNA) to build a co-express network, the modules were strong relative to weight traits. The findings highlighted the underlying way and a relative network to yield a new view about gene expression between the yaks living plateau and plain.

Genome-wide analysis of genetic diversity and artificial selection in Large White pigs in Russia

Breeding practices adopted at different farms are aimed at maximizing the profitability of pig farming. In this work, we have analyzed the genetic diversity of Large White pigs in Russia. We compared genomes of historic and modern Large White Russian breeds using 271 pig samples. We have identified 120 candidate regions associated with the differentiation of modern and historic pigs and analyzed genomic differences between the modern farms. The identified genes were associated with height, fitness, conformation, reproductive performance, and meat quality.

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.

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.

The advancements, challenges, and future implications of the CRISPR/Cas9 system in swine research

Clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (CRISPR/Cas9) genome editing technology has dramatically influenced swine research by enabling the production of high-quality disease-resistant pig breeds, thus improving yields. In addition, CRISPR/Cas9 has been used extensively in pigs as one of the tools in biomedical research. In this review, we present the advancements of the CRISPR/Cas9 system in swine research, such as animal breeding, vaccine development, xenotransplantation, and disease modeling. We also highlight the current challenges and some potential applications of the CRISPR/Cas9 technologies.

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.

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.

Transcriptomic comparison of communally reared wild, domesticated and hybrid Atlantic salmon fry under stress and control conditions

BACKGROUND

Domestication is the process by which organisms become adapted to the human-controlled environment. Since the selection pressures that act upon cultured and natural populations differ, adaptations that favour life in the domesticated environment are unlikely to be advantageous in the wild. Elucidation of the differences between wild and domesticated Atlantic salmon may provide insights into some of the genomic changes occurring during domestication, and, help to predict the evolutionary consequences of farmed salmon escapees interbreeding with wild conspecifics. In this study the transcriptome of the offspring of wild and domesticated Atlantic salmon were compared using a common-garden experiment under standard hatchery conditions and in response to an applied crowding stressor.

RESULTS

Transcriptomic differences between wild and domesticated crosses were largely consistent between the control and stress conditions, and included down-regulation of environmental information processing, immune and nervous system pathways and up-regulation of genetic information processing, carbohydrate metabolism, lipid metabolism and digestive and endocrine system pathways in the domesticated fish relative to their wild counterparts, likely reflective of different selection pressures acting in wild and cultured populations. Many stress responsive functions were also shared between crosses and included down-regulation of cellular processes and genetic information processing and up-regulation of some metabolic pathways, lipid and energy in particular. The latter may be indicative of mobilization and reallocation of energy resources in response to stress. However, functional analysis indicated that a number of pathways behave differently between domesticated and wild salmon in response to stress. Reciprocal F1 hybrids permitted investigation of inheritance patterns that govern transcriptomic differences between these genetically divergent crosses. Additivity and maternal dominance accounted for approximately 42 and 25% of all differences under control conditions for both hybrids respectively. However, the inheritance of genes differentially expressed between crosses under stress was less consistent between reciprocal hybrids, potentially reflecting maternal environmental effects.

CONCLUSION

We conclude that there are transcriptomic differences between the domesticated and wild salmon strains studied here, reflecting the different selection pressures operating on them. Our results indicate that stress may affect certain biological functions differently in wild, domesticated and hybrid crosses and these should be further investigated.

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.

Lipid Deposition and Metabolism in Local and Modern Pig Breeds: A Review

Modern pig breeds, which have been genetically improved to achieve fast growth and a lean meat deposition, differ from local pig breeds with respect to fat deposition, fat specific metabolic characteristics and various other properties. The present review aimed to elucidate the mechanisms underlying the differences between fatty local and modern lean pig breeds in adipose tissue deposition and lipid metabolism, taking into consideration morphological, cellular, biochemical, transcriptomic and proteomic perspectives. Compared to modern breeds, local pig breeds accumulate larger amounts of fat, which generally contains more monounsaturated and saturated fatty acids; they exhibit a higher adipocyte size and higher activity of lipogenic enzymes. Studies using transcriptomic and proteomic approaches highlighted several processes like immune response, fatty-acid turn-over, oxidoreductase activity, mitochondrial function, etc. which differ between local and modern pig breeds.

Assessing selection signatures within and between selected lines of dual-purpose black and white and German Holstein cattle

The aim of this study was to detect selection signatures considering cows from the German Holstein (GH) and the local dual-purpose black and white (DSN) population, as well as from generated sub-populations. The 4654 GH and 261 DSN cows were genotyped with the BovineSNP50 Genotyping BeadChip. The geographical herd location was used as an environmental descriptor to create the East-DSN and West-DSN sub-populations. In addition, two further sub-populations of GH cows were generated, using the extreme values for solutions of residual effects of cows for the claw disorder dermatitis digitalis. These groups represented the most susceptible and most resistant cows. We used cross-population extended haplotype homozygosity methodology (XP-EHH) to identify the most recent selection signatures. Furthermore, we calculated Wright's fixation index (F_ST ). Chromosomal segments for the top 0.1 percentile of negative or positive XP-EHH scores were studied in detail. For gene annotations, we used the Ensembl database and we considered a window of 250 kbp downstream and upstream of each core SNP corresponding to peaks of XP-EHH. In addition, functional interactions among potential candidate genes were inferred via gene network analyses. The most outstanding XP-EHH score was on chromosome 12 (at 77.34 Mb) for DSN and on chromosome 20 (at 36.29-38.42 Mb) for GH. Selection signature locations harbored QTL for several economically important milk and meat quality traits, reflecting the different breeding goals for GH and DSN. The average F_ST value between GH and DSN was quite low (0.068), indicating shared founders. For group stratifications according to cow health, several identified potential candidate genes influence disease resistance, especially to dermatitis digitalis.

A machine learning approach for the identification of population-informative markers from high-throughput genotyping data: application to several pig breeds

Single nucleotide polymorphisms (SNPs) able to describe population differences can be used for important applications in livestock, including breed assignment of individual animals, authentication of mono-breed products and parentage verification among several other applications. To identify the most discriminating SNPs among thousands of markers in the available commercial SNP chip tools, several methods have been used. Random forest (RF) is a machine learning technique that has been proposed for this purpose. In this study, we used RF to analyse PorcineSNP60 BeadChip array genotyping data obtained from a total of 2737 pigs of 7 Italian pig breeds (3 cosmopolitan-derived breeds: Italian Large White, Italian Duroc and Italian Landrace, and 4 autochthonous breeds: Apulo-Calabrese, Casertana, Cinta Senese and Nero Siciliano) to identify breed informative and reduced SNP panels using the mean decrease in the Gini Index and the Mean Decrease in Accuracy parameters with stability evaluation. Other reduced informative SNP panels were obtained using Delta, Fixation index and principal component analysis statistics, and their performances were compared with those obtained using the RF-defined panels using the RF classification method and its derived Out Of Bag rates and correct prediction proportions. Therefore, the performances of a total of six reduced panels were evaluated. The correct assignment of the animals to its breed was close to 100% for all tested approaches. Porcine chromosome 8 harboured the largest number of selected SNPs across all panels. Many SNPs were included in genomic regions in which previous studies identified signatures of selection or genes (e.g. ESR1, KITL and LCORL) that could contribute to explain, at least in part, phenotypically or economically relevant traits that might differentiate cosmopolitan and autochthonous pig breeds. Random forest used as preselection statistics highlighted informative SNPs that were not the same as those identified by other methods. This might be due to specific features of this machine learning methodology. It will be interesting to explore if the adaptation of RF methods for the identification of selection signature regions could be able to describe population-specific features that are not captured by other approaches.

Genome-Wide Scan Identifies Selection Signatures in Chinese Wagyu Cattle Using a High-Density SNP Array

Selective breeding can lead to genetic diversity and diverse phenotypes in farm animals. Analysis of the genomic regions under selection can provide important insights into the genetic basis of complex traits. In this study, a high-density SNP array was used for analysis of genome selection signatures in Chinese Wagyu cattle. In total, we obtained 478,903 SNPs and 24,820 no-overlap regions for |iHS| (integrated haplotype score) estimations. Under the threshold of the top 1%, 239 regions were finally identified as candidate selected regions and 162 candidate genes were found based on the UMD3.1 genome assembly. These genes were reported to be associated with fatty acids, such as Bos taurus nitric oxide synthase 1 adaptor protein (NOS1AP), Bos taurus hydroxysteroid 17-beta dehydrogenase 7 (HSD17B7), Bos taurus WD repeat domain 7 (WDR7), Bos taurus ELOVL fatty acid elongase 2 (ELOVL2), Bos taurus calpain 1 (CAPN1), Bos taurus parkin RBR E3 ubiquitin protein ligase (PRKN, also known as PARK2), Bos taurus mitogen-activated protein kinase kinase 6 (MAP2K6), meat quality, including Bos taurus ADAM metallopeptidase domain 12 (ADAM12), Bos taurus 5'-aminolevulinate synthase 1 (ALAS1), Bos taurus small integral membrane protein 13 (SMIM13) and Bos taurus potassium two pore domain channel subfamily K member 2 (KCNK2), growth, and developmental traits, such as Bos taurus insulin like growth factor 2 receptor (IGF2R), Bos taurus RAR related orphan receptor A (RORA), Bos taurus fibroblast growth factor 14 (FGF14), Bos taurus paired box 6 (PAX6) and Bos taurus LIM homeobox 6 (LHX6). In addition, we identified several genes that are associated with body size and weight, including Bos taurus sorting nexin 29 (SNX29), Bos taurus zinc finger imprinted 2 (ZIM2), Bos taurus family with sequence similarity 110 member A (FAM110A), immune system, including Bos taurus toll like receptor 9 (TLR9), Bos taurus TAFA chemokine like family member 1 (TAFA1), Bos taurus glutathione peroxidase 8 (putative) (GPX8), Bos taurus interleukin 5 (IL5), Bos taurus PR domain containing 9 (PRDM9), Bos taurus glutamate ionotropic receptor kainate type subunit 2 (GRIK2) and feed intake efficiency, Bos taurus sodium voltage-gated channel alpha subunit 9 (SCN9A), Bos taurus relaxin family peptide/INSL5 receptor 4 (RXFP4), Bos taurus RNA polymerase II associated protein 3 (RPAP3). Moreover, four GO terms of biological regulation (GO:0009987, GO:0008152) and metabolic process (GO:0003824, GO:0005488) were found based on these genes. In addition, we found that 232 candidate regions (~18 Mb) overlapped with the Quantitative trait loci (QTL)regions extracted from cattle QTLdb. Our findings imply that many genes were selected for important traits in Chinese Wagyu cattle. Moreover, these results can contribute to the understanding of the genetic basis of the studied traits during the formation of this population.

Genome-Wide Signatures of Selection Detection in Three South China Indigenous Pigs

South China indigenous pigs are famous for their superior meat quality and crude feed tolerance. Saba and Baoshan pigs without saddleback were located in the high-altitude area of Yunnan Province, while Tunchang and Ding’an pigs with saddleback were located in the low-altitude area of Hainan Province. Although these pigs are different in appearance, the underlying genetic differences have not been investigated. In this study, based on the single-nucleotide polymorphism (SNP) genotypes of 124 samples, both the cross-population extended haplotype homozygosity (XP-EHH) and the fixation index (F_ST) statistic were used to identify potential signatures of selection in these pig breeds. We found nine potential signatures of selection detected simultaneously by two methods, annotated 22 genes in Hainan pigs, when Baoshan pigs were used as the reference group. In addition, eleven potential signatures of selection detected simultaneously by two methods, annotated 24 genes in Hainan pigs compared with Saba pigs. These candidate genes were most enriched in GO: 0048015~phosphatidylinositol-mediated signaling and ssc00604: Glycosphingolipid biosynthesis—ganglio series. These selection signatures were likely to overlap with quantitative trait loci associated with meat quality traits. Furthermore, one potential selection signature, which was associated with different coat color, was detected in Hainan pigs. These results contribute to a better understanding of the underlying genetic architecture of South China indigenous pigs.

Data Mining and Validation of AMPK Pathway as a Novel Candidate Role Affecting Intramuscular Fat Content in Pigs

Simple Summary

Intramuscular fat (IMF) is increasingly being recognized as a key meat trait in the modern pork industry. The aims of this research were to identify potential signaling pathways associated with IMF content in the longissimus dorsi (LD) muscle of different pig breeds and investigate the gene expression levels in the screened signaling pathways. Our results indicated that the AMPK signaling pathway may be related to IMF deposition in the LD muscle of pigs. The results of qRT-PCR analysis showed that the expression of ten key hub genes (AMPK, ADIPOR1, ADIPOR2, LKB1, CAMKKβ, CPT1A, CPT1B, PGC-1α, CD36, and ACC1) differed between the LD muscle of Min and Large White pigs. The protein expression levels of AMPK, LKB1, CaMKK2, CPT1A, and ACC1 were similar to the genes expression patterns in the LD muscle of Large White pigs. The results of this study provide novel insights into the regulatory function of the AMPK signaling pathway in relation to the IMF content in the LD muscle of different pigs.

Abstract

Intramuscular fat (IMF) is an important economic trait for pork quality and a complex quantitative trait regulated by multiple genes. The objective of this work was to investigate the novel transcriptional effects of a multigene pathway on IMF deposition in the longissimus dorsi (LD) muscles of pigs. Potential signaling pathways were screened by mining data from three gene expression profiles in the Gene Expression Omnibus (GEO) database. We designed quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) arrays for the candidate signaling pathways to verify the results in the LD muscles of two pig breeds with different IMF contents (Large White and Min). Western blot analysis was used to detect the expression levels of several candidate proteins. Our results showed that the AMPK signaling pathway was screened via bioinformatics analysis. Ten key hub genes of this signaling pathway (AMPK, ADIPOR1, ADIPOR2, LKB1, CAMKKβ, CPT1A, CPT1B, PGC-1α, CD36, and ACC1) were differentially expressed between the Large White and Min pigs. Western blot analysis further confirmed that LKB1/CaMKK2-AMPK-ACC1-CPT1A axis dominates the activity of AMPK signaling pathway. Statistical analyses revealed that AMPK signaling pathway activity clearly varied among the two pig breeds. Based on these results, we concluded that the activation of the AMPK signaling pathway plays a positive role in reducing IMF deposition in pigs.

Population structure and genetic diversity of four Henan pig populations

To investigate the population structure and genetic diversity of Henan indigenous pig breeds, samples from a total of 78 pigs of 11 breeds were collected, including four pig populations from Henan Province, three Western commercial breeds, three Chinese native pig breeds from other provinces and one Asian wild boar. The genotyping datasets were obtained by genotyping-by-sequencing technology. We found a high degree of polymorphism and rapid linkage disequilibrium decay in Henan pigs. A neighbor-joining tree, principal component analysis and structure analysis revealed that the Huainan and Erhualian pigs were clustered together and that the Queshan black pigs were clearly grouped together but that the Nanyang and Yuxi pigs were extensively admixed with Western pigs. In addition, heterozygosity values might indicate that Henan indigenous pigs, especially the Queshan black and Huainan pigs, were subjected to little selection during domestication. The results presented here indicate that Henan pig breeds were admixed from Western breeds, especially Nanyang and Yuxi pigs. Therefore, establishment of purification and rejuvenation systems to implement conservation strategies is urgent. In addition, it is also necessary to accelerate genetic resources improvement and utilization using modern breeding technologies, such as genomic selection and genome-wide association studies.

Comparing genomic signatures of domestication in two Atlantic salmon (Salmo salar L.) populations with different geographical origins

Selective breeding and genetic improvement have left detectable signatures on the genomes of domestic species. The elucidation of such signatures is fundamental for detecting genomic regions of biological relevance to domestication and improving management practices. In aquaculture, domestication was carried out independently in different locations worldwide, which provides opportunities to study the parallel effects of domestication on the genome of individuals that have been selected for similar traits. In this study, we aimed to detect potential genomic signatures of domestication in two independent pairs of wild/domesticated Atlantic salmon populations of Canadian and Scottish origins, respectively. Putative genomic regions under divergent selection were investigated using a 200K SNP array by combining three different statistical methods based either on allele frequencies (LFMM, Bayescan) or haplotype differentiation (Rsb). We identified 337 and 270 SNPs potentially under divergent selection in wild and hatchery populations of Canadian and Scottish origins, respectively. We observed little overlap between results obtained from different statistical methods, highlighting the need to test complementary approaches for detecting a broad range of genomic footprints of selection. The vast majority of the outliers detected were population-specific but we found four candidate genes that were shared between the populations. We propose that these candidate genes may play a role in the parallel process of domestication. Overall, our results suggest that genetic drift may have override the effect of artificial selection and/or point toward a different genetic basis underlying the expression of similar traits in different domesticated strains. Finally, it is likely that domestication may predominantly target polygenic traits (e.g., growth) such that its genomic impact might be more difficult to detect with methods assuming selective sweeps.

Combined QTL and Genome Scan Analyses With the Help of 2b-RAD Identify Growth-Associated Genetic Markers in a New Fast-Growing Carp Strain

Common carp is one of the oldest and most popular cultured freshwater fish species both globally and in China. In a previous study, we used a carp strain with a long breeding tradition in China, named Huanghe, to create a new fast-growing strain by selection for fast growth for 6 years. The growth performance at 8 months of age has been improved by 20.84%. To achieve this, we combined the best linear unbiased prediction with marker-assisted selection techniques. Recent progress in genome-wide association studies and genomic selection in livestock breeding inspired common carp breeders to consider genome-based breeding approaches. In this study, we developed a 2b-RAD sequence assay as a means of investigating the quantitative trait loci in common carp. A total of 4,953,017,786 clean reads were generated for 250 specimens (average reads/specimen = 19,812,071) with BsaXI Restriction Enzyme. From these, 56,663 SNPs were identified, covering 50 chromosomes and 3,377 scaffolds. Principal component analysis indicated that selection and control groups are relatively clearly distinct. Top 1% of Fst values was selected as the threshold signature of artificial selection. Among the 244 identified loci, genes associated with sex-related factors and nutritional metabolism (especially fat metabolism) were annotated. Eighteen QTL were associated with growth parameters. Body length at 3 months of age and body weight (both at 3 and 8 months) were controlled by polygenic effects, but body size (length, depth, width) at 8 months of age was controlled mainly by several loci with major effects. Importantly, a single shared QTL (IGF2 gene) partially controlled the body length, depth, and width. By merging the above results, we concluded that mainly the genes related to neural pathways, sex and fatty acid metabolism contributed to the improved growth performance of the new Huanghe carp strain. These findings are one of the first investigations into the potential use of genomic selection in the breeding of common carp. Moreover, our results show that combining the Fst, QTL mapping and CRISPR–Cas9 methods can be an effective way to identify important novel candidate molecular markers in economic breeding programs.

Selection Signatures Underlying Dramatic Male Inflorescence Transformation During Modern Hybrid Maize Breeding

Inflorescence capacity plays a crucial role in reproductive fitness in plants, and in production of hybrid crops. Maize is a monoecious species bearing separate male and female flowers (tassel and ear, respectively). The switch from open-pollinated populations of maize to hybrid-based breeding schemes in the early 20th century was accompanied by a dramatic reduction in tassel size, and the trend has continued with modern breeding over the recent decades. The goal of this study was to identify selection signatures in genes that may underlie this dramatic transformation. Using a population of 942 diverse inbred maize accessions and a nested association mapping population comprising three 200-line biparental populations, we measured 15 tassel morphological characteristics by manual and image-based methods. Genome-wide association studies identified 242 single nucleotide polymorphisms significantly associated with measured traits. We compared 41 unselected lines from the Iowa Stiff Stalk Synthetic (BSSS) population to 21 highly selected lines developed by modern commercial breeding programs, and found that tassel size and weight were reduced significantly. We assayed genetic differences between the two groups using three selection statistics: cross population extended haplotype homozogysity, cross-population composite likelihood ratio, and fixation index. All three statistics show evidence of selection at genomic regions associated with tassel morphology relative to genome-wide null distributions. These results support the tremendous effect, both phenotypic and genotypic, that selection has had on maize male inflorescence morphology.

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

AIM

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

MATERIALS & METHODS

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

RESULTS

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

CONCLUSION

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genomic analysis reveals genes affecting distinct phenotypes among different Chinese and western pig breeds

The differences in artificial and natural selection have been some of the factors contributing to phenotypic diversity between Chinese and western pigs. Here, 830 individuals from western and Chinese pig breeds were genotyped using the reduced-representation genotyping method. First, we identified the selection signatures for different pig breeds. By comparing Chinese pigs and western pigs along the first principal component, the growth gene IGF1R; the immune genes IL1R1, IL1RL1, DUSP10, RAC3 and SWAP70; the meat quality-related gene SNORA50 and the olfactory gene OR1F1 were identified as candidate differentiated targets. Further, along a principal component separating Pudong White pigs from others, a potential causal gene for coat colour (EDNRB) was discovered. In addition, the divergent signatures evaluated by F_st within Chinese pig breeds found genes associated with the phenotypic features of coat colour, meat quality and feed efficiency among these indigenous pigs. Second, admixture and genomic introgression analysis were performed. Shan pigs have introgressed genes from Berkshire, Yorkshire and Hongdenglong pigs. The results of introgression mapping showed that this introgression conferred adaption to the local environment and coat colour of Chinese pigs and the superior productivity of western pigs.

Genome-Wide Signatures of Selection Reveal Genes Associated With Performance in American Quarter Horse Subpopulations

Selective breeding for athletic performance in various disciplines has resulted in population stratification within the American Quarter Horse (QH) breed. The goals of this study were to utilize high density genotype data to: (1) identify genomic regions undergoing positive selection within and among QH subpopulations; (2) investigate haplotype structure within each QH subpopulation; and (3) identify candidate genes within genomic regions of interest (ROI), as well as biological pathways, predicted to play a role in elite performance in each group. For that, 65K SNP genotyping data on 143 elite individuals from 6 QH subpopulations (cutting, halter, racing, reining, western pleasure, and working cow) were imputed to 2M SNPs. Signatures of selection were identified using F_ST-based (d_i ) and haplotype-based (hapFLK) analyses, accompanied by identification of local haplotype structure and sharing within subpopulations (hapQTL). Regions undergoing positive selection were identified on all 31 autosomes, and ROI on 2 chromosomes were identified by all 3 methods combined. Genes within each ROI were retrieved and used to identify pathways and genes that might contribute to performance in each subpopulation. These included, among others, candidate genes associated with skeletal muscle development, metabolism, and central nervous system development. This work improves our understanding of equine breed development, and provides breeders with a better understanding of how selective breeding impacts the performance of QH populations.

Genome-wide detection of selection signatures in Chinese indigenous Laiwu pigs revealed candidate genes regulating fat deposition in muscle

Background

Currently, genome-wide scans for positive selection signatures in commercial breed have been investigated. However, few studies have focused on selection footprints of indigenous breeds. Laiwu pig is an invaluable Chinese indigenous pig breed with extremely high proportion of intramuscular fat (IMF), and an excellent model to detect footprint as the result of natural and artificial selection for fat deposition in muscle.

Result

In this study, based on GeneSeek Genomic profiler Porcine HD data, three complementary methods, F_ST, iHS (integrated haplotype homozygosity score) and CLR (composite likelihood ratio), were implemented to detect selection signatures in the whole genome of Laiwu pigs. Totally, 175 candidate selected regions were obtained by at least two of the three methods, which covered 43.75 Mb genomic regions and corresponded to 1.79% of the genome sequence. Gene annotation of the selected regions revealed a list of functionally important genes for feed intake and fat deposition, reproduction, and immune response. Especially, in accordance to the phenotypic features of Laiwu pigs, among the candidate genes, we identified several genes, NPY1R, NPY5R, PIK3R1 and JAKMIP1, involved in the actions of two sets of neurons, which are central regulators in maintaining the balance between food intake and energy expenditure.

Conclusions

Our results identified a number of regions showing signatures of selection, as well as a list of functionally candidate genes with potential effect on phenotypic traits, especially fat deposition in muscle. Our findings provide insights into the mechanisms of artificial selection of fat deposition and further facilitate follow-up functional studies.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0622-y) contains supplementary material, which is available to authorized users.

Detection of Selection Signatures in Chinese Landrace and Yorkshire Pigs Based on Genotyping-by-Sequencing Data

The domestic pigs have been undergone intense selection pressures for these development of interested traits following domestication and modern breeding. This has altered many traits in most of pig breeds, such as growth rate, body weight, fertility, and immunity. Thus, the objectives of this study were to (1) detect these selection signatures and identify the candidate genes that show evidences of recent artificial selection at the level of whole genome, (2) be beneficial to understand the relationship between genomic structure and phenotypic diversity, and (3) highlight the key roles of these candidate genes in growth and development in the two breeds. The data consisted of total raw number of 345570 single nucleotide polymorphisms (SNPs) in 1200 individuals from the Chinese Landrace pigs (L, n = 600) and Yorkshire pigs (Y, n = 600). Based on these SNPs data, two complementary methods, population differentiation (Fst) and composite likelihood ratio test (CLR), were carried out to detect the selection signatures in this study. A total of 540 potential selection regions (50 kb) which contained 111 candidate genes were detected for Landrace-Yorkshire pair (L-Y) by Fst. In addition, 73 and 125 candidate genes were found for Landrace pigs and Yorkshire pigs by CLR test based on 321 and 628 potential selection regions, respectively. Some candidate genes are associated with important traits and signaling pathways including the ACACA, MECR, COL11A1, GHR, IGF1R, IGF2R, IFNG, and MTOR gene. The ACACA and MECR gene are related to fatty acid biosynthesis. The COL11A1 gene is essential for the development of the normal differentiation. The GHR, IGF1R, and IGF2R gene are significant candidate genes which play major roles in the growth and development in animals. The IFNG gene is associated with some aspects of immune response. The MTOR gene regulates many signaling pathways and signaling transduction pathway.

Linkage disequilibrium and functional analysis of PRE1 insertion together with SNPs in the promoter region of IGFBP7 gene in different pig breeds

Polymorphisms in regions upstream of transcription initiation site may modify the transcriptional activity of target genes by changing promoter activity. This study aims to determine whether or not polymorphisms at porcine IGFBP7 promoter region affect gene expression. In this study, eight SNPs and one PRE1 insertion in this region were first confirmed. The PRE1 insertion was widespread in 20 Chinese indigenous breeds, but was not observed in three commercial breeds. A perfect linkage disequilibrium, consisting of six of those SNPs and a PRE1, was observed with two haplotypes (h1 and h2) in five pig breeds. The h1 haplotype had an overwhelming superiority distribution in Large White, Landrace, and Bama mini-pig; in turn, the h2 only existed in the PRE1 presence breeds. As the haplotypes and PRE1 were located at gene promoter regions, we further investigated the transfection of plasmids with three different fragments of IGFBP-7 promoter region (H1, H2, RF). The CMV promoter of the pEGFP-N1 was substituted by these three different fragments, respectively. Different transcriptional and translational activities of EGFP in PK-15 cells were observed in these three recombinant plasmids by quantitative real-time PCR and flow cytometric analysis. The results indicated that H1 had the higher transcriptional and translational activities of EGFP as compared to the H2 (P < 0.05, P < 0.05). As compared to the RF group, EGFP mRNA expression level was significantly higher in H1 groups (P < 0.05). The IGFBP-7 promoter polymorphisms detected in this study may be important functional variants and potential genetic markers for pig population genetic study.

Detection of genomic structural variations in Guizhou indigenous pigs and the comparison with other breeds

Genomic structural variation (SV) is noticed for the contribution to genetic diversity and phenotypic changes. Guizhou indigenous pig (GZP) has been raised for hundreds of years with many special characteristics. The present paper aimed to uncover the influence of SV on gene polymorphism and the genetic mechanisms of phenotypic traits for GZP. Eighteen GZPs were chosen for resequencing by Illumina sequencing platform. The confident SVs of GZP were called out by both programs of pindel and softSV simultaneously and compared with the SVs deduced from the genomic data of European pig (EUP) and the native pig outside of Guizhou, China (NPOG). A total of 39,166 SVs were detected and covered 27.37 Mb of pig genome. All of 76 SVs were confirmed in GZP pig population by PCR method. The SVs numbers in NPOG and GZP were about 1.8 to 1.9 times higher than that in EUP. And a SV hotspot was found out from the 20 Mb of chromosome X of GZP, which harbored 29 genes and focused on histone modification. More than half of SVs was positioned in the intergenic regions and about one third of SVs in the introns of genes. And we found that SVs tended to locate in genes produced multi-transcripts, in which a positive correlation was found out between the numbers of SV and the gene transcripts. It illustrated that the primary mode of SVs might function on the regulation of gene expression or the transcripts splicing process. A total of 1,628 protein-coding genes were disturbed by 1,956 SVs specific in GZP, in which 93 GZP-specific SV-related genes would lose their functions due to the SV interference and gathered in reproduction ability. Interestingly, the 1,628 protein-coding genes were mainly enriched in estrogen receptor binding, steroid hormone receptor binding, retinoic acid receptor binding, oxytocin signaling pathway, mTOR signaling pathway, axon guidance and cholinergic synapse pathways. It suggested that SV might be a reason for the strong adaptability and low fecundity of GZP, and 51 candidate genes would be useful for the configuration phenotype in Xiang pig breed.

A genome scan for selection signatures comparing farmed Atlantic salmon with two wild populations: Testing colocalization among outlier markers, candidate genes, and quantitative trait loci for production traits

Comparative genome scans can be used to identify chromosome regions, but not traits, that are putatively under selection. Identification of targeted traits may be more likely in recently domesticated populations under strong artificial selection for increased production. We used a North American Atlantic salmon 6K SNP dataset to locate genome regions of an aquaculture strain (Saint John River) that were highly diverged from that of its putative wild founder population (Tobique River). First, admixed individuals with partial European ancestry were detected using STRUCTURE and removed from the dataset. Outlier loci were then identified as those showing extreme differentiation between the aquaculture population and the founder population. All Arlequin methods identified an overlapping subset of 17 outlier loci, three of which were also identified by BayeScan. Many outlier loci were near candidate genes and some were near published quantitative trait loci (QTLs) for growth, appetite, maturity, or disease resistance. Parallel comparisons using a wild, nonfounder population (Stewiacke River) yielded only one overlapping outlier locus as well as a known maturity QTL. We conclude that genome scans comparing a recently domesticated strain with its wild founder population can facilitate identification of candidate genes for traits known to have been under strong artificial selection.

Genome-wide analysis reveals signatures of selection for important traits in domestic sheep from different ecoregions

Background

Throughout a long period of adaptation and selection, sheep have thrived in a diverse range of ecological environments. Mongolian sheep is the common ancestor of the Chinese short fat-tailed sheep. Migration to different ecoregions leads to changes in selection pressures and results in microevolution. Mongolian sheep and its subspecies differ in a number of important traits, especially reproductive traits. Genome-wide intraspecific variation is required to dissect the genetic basis of these traits.

Results

This research resequenced 3 short fat-tailed sheep breeds with a 43.2-fold coverage of the sheep genome. We report more than 17 million single nucleotide polymorphisms and 2.9 million indels and identify 143 genomic regions with reduced pooled heterozygosity or increased genetic distance to each other breed that represent likely targets for selection during the migration. These regions harbor genes related to developmental processes, cellular processes, multicellular organismal processes, biological regulation, metabolic processes, reproduction, localization, growth and various components of the stress responses. Furthermore, we examined the haplotype diversity of 3 genomic regions involved in reproduction and found significant differences in TSHR and PRL gene regions among 8 sheep breeds.

Conclusions

Our results provide useful genomic information for identifying genes or causal mutations associated with important economic traits in sheep and for understanding the genetic basis of adaptation to different ecological environments.

Electronic supplementary material

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

Comparing the transcriptomes of embryos from domesticated and wild Atlantic salmon (Salmo salar L.) stocks and examining factors that influence heritability of gene expression

BACKGROUND

Due to selective breeding, domesticated and wild Atlantic salmon are genetically diverged, which raises concerns about farmed escapees having the potential to alter the genetic composition of wild populations and thereby disrupting local adaptation. Documenting transcriptional differences between wild and domesticated stocks under controlled conditions is one way to explore the consequences of domestication and selection. We compared the transcriptomes of wild and domesticated Atlantic salmon embryos, by using a custom 44k oligonucleotide microarray to identify perturbed gene pathways between the two stocks, and to document the inheritance patterns of differentially-expressed genes by examining gene expression in their reciprocal hybrids.

RESULTS

Data from 24 array interrogations were analysed: four reciprocal cross types (W♀ × W♂, D♀ × W♂; W♀ × D♂, D♀ × D♂) × six biological replicates. A common set of 31,491 features on the microarrays passed quality control, of which about 62 % were assigned a KEGG Orthology number. A total of 6037 distinct genes were identified for gene-set enrichment/pathway analysis. The most highly enriched functional groups that were perturbed between the two stocks were cellular signalling and immune system, ribosome and RNA transport, and focal adhesion and gap junction pathways, relating to cell communication and cell adhesion molecules. Most transcripts that were differentially expressed between the stocks were governed by additive gene interaction (33 to 42 %). Maternal dominance and over-dominance were also prevalent modes of inheritance, with no convincing evidence for a stock effect.

CONCLUSIONS

Our data indicate that even at this relatively early developmental stage, transcriptional differences exist between the two stocks and affect pathways that are relevant to wild versus domesticated environments. Many of the identified differentially perturbed pathways are involved in organogenesis, which is expected to be an active process at the eyed egg stage. The dominant effects are more largely due to the maternal line than to the origin of the stock. This finding is particularly relevant in the context of potential introgression between farmed and wild fish, since female escapees tend to have a higher spawning success rate compared to males.

Structural Variant Detection by Large-scale Sequencing Reveals New Evolutionary Evidence on Breed Divergence between Chinese and European Pigs

In this study, we performed a genome-wide SV detection among the genomes of thirteen pigs from diverse Chinese and European originated breeds by next genetation sequencing, and constrcuted a single-nucleotide resolution map involving 56,930 putative SVs. We firstly identified a SV hotspot spanning 35 Mb region on the X chromosome specifically in the genomes of Chinese originated individuals. Further scrutinizing this region by large-scale sequencing data of extra 111 individuals, we obtained the confirmatory evidence on our initial finding. Moreover, thirty five SV-related genes within the hotspot region, being of importance for reproduction ability, rendered significant different evolution rates between Chinese and European originated breeds. The SV hotspot identified herein offers a novel evidence for assessing phylogenetic relationships, as well as likely explains the genetic difference of corresponding phenotypes and features, among Chinese and European pig breeds. Furthermore, we employed various SVs to infer genetic structure of individuls surveyed. We found SVs can clearly detect the difference of genetic background among individuals. This clues us that genome-wide SVs can capture majority of geneic variation and be applied into cladistic analyses. Characterizing whole genome SVs demonstrated that SVs are significantly enriched/depleted with various genomic features.

Signatures of positive selection in East African Shorthorn Zebu: A genome-wide single nucleotide polymorphism analysis

The small East African Shorthorn Zebu (EASZ) is the main indigenous cattle across East Africa. A recent genome wide SNP analysis revealed an ancient stable African taurine x Asian zebu admixture. Here, we assess the presence of candidate signatures of positive selection in their genome, with the aim to provide qualitative insights about the corresponding selective pressures. Four hundred and twenty-five EASZ and four reference populations (Holstein-Friesian, Jersey, N’Dama and Nellore) were analysed using 46,171 SNPs covering all autosomes and the X chromosome. Following F_ST and two extended haplotype homozygosity-based (iHS and Rsb) analyses 24 candidate genome regions within 14 autosomes and the X chromosome were revealed, in which 18 and 4 were previously identified in tropical-adapted and commercial breeds, respectively. These regions overlap with 340 bovine QTL. They include 409 annotated genes, in which 37 were considered as candidates. These genes are involved in various biological pathways (e.g. immunity, reproduction, development and heat tolerance). Our results support that different selection pressures (e.g. environmental constraints, human selection, genome admixture constrains) have shaped the genome of EASZ. We argue that these candidate regions represent genome landmarks to be maintained in breeding programs aiming to improve sustainable livestock productivity in the tropics.

Detection of selection signatures in Piemontese and Marchigiana cattle, two breeds with similar production aptitudes but different selection histories

BACKGROUND

Domestication and selection are processes that alter the pattern of within- and between-population genetic variability. They can be investigated at the genomic level by tracing the so-called selection signatures. Recently, sequence polymorphisms at the genome-wide level have been investigated in a wide range of animals. A common approach to detect selection signatures is to compare breeds that have been selected for different breeding goals (i.e. dairy and beef cattle). However, genetic variations in different breeds with similar production aptitudes and similar phenotypes can be related to differences in their selection history.

METHODS

In this study, we investigated selection signatures between two Italian beef cattle breeds, Piemontese and Marchigiana, using genotyping data that was obtained with the Illumina BovineSNP50 BeadChip. The comparison was based on the fixation index (Fst), combined with a locally weighted scatterplot smoothing (LOWESS) regression and a control chart approach. In addition, analyses of Fst were carried out to confirm candidate genes. In particular, data were processed using the varLD method, which compares the regional variation of linkage disequilibrium between populations.

RESULTS

Genome scans confirmed the presence of selective sweeps in the genomic regions that harbour candidate genes that are known to affect productive traits in cattle such as DGAT1, ABCG2, CAPN3, MSTN and FTO. In addition, several new putative candidate genes (for example ALAS1, ABCB8, ACADS and SOD1) were detected.

CONCLUSIONS

This study provided evidence on the different selection histories of two cattle breeds and the usefulness of genomic scans to detect selective sweeps even in cattle breeds that are bred for similar production aptitudes.

A genome scan for selection signatures in pigs

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