Genome-wide footprints of pig domestication and selection revealed through massive parallel sequencing of pooled DNA

Population genetics reveals new introgression in the nucleus herd of min pigs

OBJECTIVE

Min pigs are a unique genetic resource among local pig breeds in China. They have more excellent characteristics in cold and stress resistance, good meat quality, and a high reproductive rate. However, the genetic structure and driving factors remain unclear in the nucleus herd. In this study, the genetic diversity of Min pigs was studied to reveal the formation mechanism of its unique genetic structure. We hope to protect and develop the genetic resources of Min pigs.

METHODS

We analyzed different types of genes to identify the genetic structure and gene introgression pattern of Min pigs. The nuclear DNA dataset includes information on 21 microsatellite loci and 6 Y-chromosome genes, and the mitochondrial D-loop gene is selected to represent maternal lineages. The above genes are all from the nucleus herd of Min pigs.

RESULTS

The results of genetic structure identification and analysis of potential exogenous gene introgression patterns indicate that the nucleus herd of Min pigs maintains a high level of genetic diversity (polymorphism information content = 0.713, expected heterozygosity = 0.662, observed heterozygosity = 0.612). Compared with other Asian pig breeds, the formation of Min pig breeds is more special. Gene introgression from European pig breeds to Min pigs has occurred, which is characterized by complete introgression of paternal genes and incomplete introgression of maternal genes.

CONCLUSION

Gene introgression caused by cross-breeding is not the main factor leading to the formation of the current genetic structure of Min pigs, but this process has increased the level of genetic diversity in the nucleus herd. Compared with the influence of gene introgression, our research suggest that artificial selection and environmental adaptive evolution make Min pigs form unique genetic characteristics.

Deciphering climate resilience in Indian cattle breeds by selection signature analyses

The signature of selection is a crucial concept in evolutionary biology that refers to the pattern of genetic variation which arises in a population due to natural selection. In the context of climate adaptation, the signature of selection can reveal the genetic basis of adaptive traits that enable organisms to survive and thrive in changing environmental conditions. Breeds living in diverse agroecological zones exhibit genetic "footprints" within their genomes that mirror the influence of climate-induced selective pressures, subsequently impacting phenotypic variance. It is assumed that the genomes of animals residing in these regions have been altered through selection for various climatic adaptations. These regions are known as signatures of selection and can be identified using various summary statistics. We examined genotypic data from eight different cattle breeds (Gir, Hariana, Kankrej, Nelore, Ongole, Red Sindhi, Sahiwal, and Tharparkar) that are adapted to diverse regional climates. To identify selection signature regions in this investigation, we used four intra-population statistics: Tajima's D, CLR, iHS, and ROH. In this study, we utilized Bovine 50 K chip data and four genome scan techniques to assess the genetic regions of positive selection for high-temperature adaptation. We have also performed a genome-wide investigation of genetic diversity, inbreeding, and effective population size in our target dataset. We identified potential regions for selection that are likely to be caused by adverse climatic conditions. We observed many adaptation genes in several potential selection signature areas. These include genes like HSPB2, HSPB3, HSP20, HSP90AB1, HSF4, HSPA1B, CLPB, GAP43, MITF, and MCHR1 which have been reported in the cattle populations that live in varied climatic regions. The findings demonstrated that genes involved in disease resistance and thermotolerance were subjected to intense selection. The findings have implications for marker-assisted breeding, understanding the genetic landscape of climate-induced adaptation, putting breeding and conservation programs into action.

Identification and analysis of circRNAs in the prefrontal cortices of wild boar and domestic pig

Domestication caused significant differences in morphology and behavior between wild and domestic pigs. However, the regulatory role of circRNA in this event is unclear. Here, we analyzed circRNA expression patterns in the prefrontal cortices of wild boar and domestic pigs to determine the potential role of circRNAs in domestication. We identified a total of 11,375 circRNAs and found that 349 and 354 circRNAs were up-regulated in wild boar and Rongchang pig, respectively. Functional enrichment analysis showed that host genes of significantly highly-expressed circRNAs in wild boar were significantly enriched in neural synapse-related categories and the categories of 'regulation of defense response (p = 0.028)' and 'neural retina development (p = 4.32 × 10-3)'. Host genes of significantly highly-expressed circRNAs in Rongchang pig were specifically involved in 'chordate embryonic development (p = 2.38 × 10-4)'. Additionally, we constructed circRNA-miRNA-mRNA regulatory axes in wild boar and Rongchang pig and found more regulatory axes in wild boar that potentially regulate synaptic activities. We identified multiple circRNAs that may be related to domesticated characteristics, such as ssc_circ_6179 (ssc_circ_6179-ssc-miR-9847-HRH3, related to aggression) and ssc_circ_3027 (ssc_circ_3027-ssc-miR-4334-5p-HCRTR1, related to attention). This study provides a resource for further investigation of the molecular basis of pig domestication.

Population Genetic Analysis of Six Chinese Indigenous Pig Meta-Populations Based on Geographically Isolated Regions

The diversification of indigenous pig breeds in China has resulted from multiple climate, topographic, and human cultural influences. The numerous indigenous pig breeds can be geographically divided into six meta-populations; however, their genetic relationships, contributions to genetic diversity, and genetic signatures remain unclear. Whole-genome SNP data for 613 indigenous pigs from the six Chinese meta-populations were obtained and analyzed. Population genetic analyses confirmed significant genetic differentiation and a moderate mixture among the Chinese indigenous pig meta-populations. The North China (NC) meta-population had the largest contribution to genetic and allelic diversity. Evidence from selective sweep signatures revealed that genes related to fat deposition and heat stress response (EPAS1, NFE2L2, VPS13A, SPRY1, PLA2G4A, and UBE3D) were potentially involved in adaptations to cold and heat. These findings from population genetic analyses provide a better understanding of indigenous pig characteristics in different environments and a theoretical basis for future work on the conservation and breeding of Chinese indigenous pigs.

Genomic signatures of artificial selection in fecundity of Pacific white shrimp, Penaeus vannamei

Penaeusvannamei is the most important economic shrimp in the world. Many selective breeding programs are carried out to improve its production and performance traits. Although significant differences in the reproductive ability of female P. vannamei under artificial breeding conditions have been reported, the genome-wide adaption of the reproductive ability of domesticated female P. vannamei is less investigated. In this study, whole-genome analysis was performed along with pooled DNA sequencing on two fecundity separated bulks, high fecundity bulk (HB), and low fecundity bulk (LB). Each bulk contained 30 individuals from 3 commercial populations. A sequencing depth of >30× was achieved for each bulk, leading to the identification of 625,181 and 629,748 single nucleotide polymorphisms (SNPs) in HB and LB, respectively. Fixation index (Fst) combined with p ratio allowed for the identification of 145 selective sweep regions, with a sequence length of 14.5 Mb, accounting for 0.59% of the genome. Among the 145 selective sweep regions, a total of 64,046 SNPs were identified, and further verification was performed by genotyping 50 candidate SNPs on 60 samples from the offspring of the three populations. Furthermore, 121 genes were screened from the sweep regions. GO annotation and KEGG enrichment analyses showed that partial genes were essential for fecundity regulation. This study provides important information for in-depth investigation of genomic characteristics for long-term selective breeding on the fecundity of female P. vannamei and will also be important for genome-assisted breeding of P. vannamei in the future.

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

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

Population structure of 3907 worldwide pigs and the introgression of Chinese indigenous pigs by European pigs

With the improvement in sequencing technology and the decrease in sequencing cost, increasing amounts of genomic data for pigs have been uploaded to public databases. However, no researchers have to date integrated all currently available data to uncover the global genetic status of pigs. Meanwhile, little is known about the introgression from European to Chinese pigs and its underlying influences. Therefore, we integrated the effective genotype data of 3907 pigs from 193 populations worldwide using population genetic analysis, gene flow analysis and a sharing-IBD study. These findings illustrate not only the population structure of 59 Chinese native breeds and others but also the amounts of gene flow and introgression that have occurred between Western and Chinese pigs. In addition, we demonstrate the presence of introgressed European haplotypes in Chinese indigenous breeds and identify relevant introgressed regions that contain genes associated with growth and feed efficiency. Moreover, we compare the introgression patterns of Western and Chinese pigs and further discuss possible explanations for why the level of introgression differs between Chinese pig breeds and Western modern breeds. Collectively, this study provides a fine global population structure analysis of pigs and presents evidence of European pigs being interbred with local breeds in China.

Detecting the differential genomic variants using cross-population phenotype-associated variant (XP-PAV) of the Landrace and Yorkshire pigs in Korea

Although there have been many genome-wide association studies (GWAS) and selective sweep analyses to understand pig genomic regions related to growth performance, these methods considered only the gene effect and selection signal, respectively. In this study, we suggest the cross-population phenotype associated variant (XP-PAV) analysis as a novel method to determine the genomic variants with different effects between the two populations. XP-PAV analysis could reveal the differential genetic variants between the two populations by considering the gene effect and selection signal simultaneously. In this study, we used daily weight gain (DWG) and back fat thickness (BF) as phenotypes and the Landrace and Yorkshire populations were used for XP-PAV analysis. The main aim was to reveal the differential selection by considering the gene effect between Landrace and Yorkshire pigs. In the gene ontology analysis of XP-PAV results, differential selective genes in DWG analysis were involved in the regulation of interleukin-2 production and cell cycle G2/M transition. The protein modification and glycerophospholipid biosynthetic processes were the most enriched terms in the BF analysis. Therefore, we could identify genetic differences for immune and several metabolic pathways between Landrace and Yorkshire breeds using the XP-PAV analysis. In this study, we expect that XP-PAV analysis will play a role in determining useful selective variants with gene effects and provide a new interpretation of the genetic differences between the two populations.

A Mix of Old British and Modern European Breeds: Genomic Prediction of Breed Composition of Smallholder Pigs in Uganda

Pig herds in Africa comprise genotypes ranging from local ecotypes to commercial breeds. Many animals are composites of these two types and the best levels of crossbreeding for particular production systems are largely unknown. These pigs are managed without structured breeding programs and inbreeding is potentially limiting. The objective of this study was to quantify ancestry contributions and inbreeding levels in a population of smallholder pigs in Uganda. The study was set in the districts of Hoima and Kamuli in Uganda and involved 422 pigs. Pig hair samples were taken from adult and growing pigs in the framework of a longitudinal study investigating productivity and profitability of smallholder pig production. The samples were genotyped using the porcine GeneSeek Genomic Profiler (GGP) 50K SNP Chip. The SNP data was analyzed to infer breed ancestry and autozygosity of the Uganda pigs. The results showed that exotic breeds (modern European and old British) contributed an average of 22.8% with a range of 2-50% while "local" blood contributed 69.2% (36.9-95.2%) to the ancestry of the pigs. Runs of homozygosity (ROH) greater than 2 megabase (Mb) quantified the average genomic inbreeding coefficient of the pigs as 0.043. The scarcity of long ROH indicated low recent inbreeding. We conclude that the genomic background of the pig population in the study is a mix of old British and modern pig ancestries. Best levels of admixture for smallholder pigs are yet to be determined, by linking genotypes and phenotypic records.

Tracing selection signatures in the pig genome gives evidence for selective pressures on a unique curly hair phenotype in Mangalitza

Selection for desirable traits and breed-specific phenotypes has left distinctive footprints in the genome of pigs. As representative of a breed with strong selective traces aiming for robustness, health and performance, the Mangalitza pig, a native curly-haired pig breed from Hungary, was investigated in this study. Whole genome sequencing and SNP chip genotyping was performed to detect runs of homozygosity (ROH) in Mangalitza and Mangalitza-crossbreeds. We identified breed specific ROH regions harboring genes associated with the development of the curly hair type and further characteristics of this breed. Further analysis of two matings of Mangalitza with straight-coated pig breeds confirmed an autosomal dominant inheritance of curly hair. Subsequent scanning of the genome for variant effects on this trait revealed two variants potentially affecting hair follicle development and differentiation. Validation in a large sample set as well as in imputed SNP data confirmed these variants to be Mangalitza-specific. Herein, we demonstrated how strong artificial selection has shaped the genome in Mangalitza pigs and left traces in the form of selection signatures. This knowledge on genomic variation promoting unique phenotypes like curly hair provides an important resource for futures studies unraveling genetic effects for special characteristics in livestock.

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

In the process of pig genetic improvement, different commercial breeds have been bred for the same purpose, improving meat production. Most of the economic traits, such as growth and fertility, have been selected similarly despite the discrepant selection pressure, which is known as parallel selection. Here, 28 whole-genome sequencing data of Danish large white pigs with an approximately 25-fold depth each were generated, resulting in about 12 million high-quality SNPs for each individual. Combined with the sequencing data of 27 Duroc and 23 European wild boars, we investigated the parallel selection of Danish large white and Duroc pigs using two complementary methods, Fst and iHS. In total, 67 candidate regions were identified as the signatures of parallel selection. The genes in candidate regions of parallel selection were mainly associated with sensory perception, growth rate, and body size. Further functional annotation suggested that the striking consistency of the terms may be caused by the polygenetic basis of quantitative traits, and revealing the complex genetic basis of parallel selection. Besides, some unique terms were enriched in population-specific selection regions, such as the limb development-related terms enriched in Duroc-specific selection regions, suggesting unique selections of breed-specific selected traits. These results will help us better understand the parallel selection process of different breeds. Moreover, we identified several potential causal SNPs that may contribute to the pig genetic breeding process.

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.

Selection signatures in tropical cattle are enriched for promoter and coding regions and reveal missense mutations in the damage response gene HELB

Background

Distinct domestication events, adaptation to different climatic zones, and divergent selection in productive traits have shaped the genomic differences between taurine and indicine cattle. In this study, we assessed the impact of artificial selection and environmental adaptation by comparing whole-genome sequences from European taurine and Asian indicine breeds and from African cattle. Next, we studied the impact of divergent selection by exploiting predicted and experimental functional annotation of the bovine genome.

Results

We identified selective sweeps in beef cattle taurine and indicine populations, including a 430-kb selective sweep on indicine cattle chromosome 5 that is located between 47,670,001 and 48,100,000 bp and spans five genes, i.e. HELB, IRAK3, ENSBTAG00000026993, GRIP1 and part of HMGA2. Regions under selection in indicine cattle display significant enrichment for promoters and coding genes. At the nucleotide level, sites that show a strong divergence in allele frequency between European taurine and Asian indicine are enriched for the same functional categories. We identified nine single nucleotide polymorphisms (SNPs) in coding regions that are fixed for different alleles between subspecies, eight of which were located within the DNA helicase B (HELB) gene. By mining information from the 1000 Bull Genomes Project, we found that HELB carries mutations that are specific to indicine cattle but also found in taurine cattle, which are known to have been subject to indicine introgression from breeds, such as N’Dama, Anatolian Red, Marchigiana, Chianina, and Piedmontese. Based on in-house genome sequences, we proved that mutations in HELB segregate independently of the copy number variation HMGA2-CNV, which is located in the same region.

Conclusions

Major genomic sequence differences between Bos taurus and Bos indicus are enriched for promoter and coding regions. We identified a 430-kb selective sweep in Asian indicine cattle located on chromosome 5, which carries SNPs that are fixed in indicine populations and located in the coding sequences of the HELB gene. HELB is involved in the response to DNA damage including exposure to ultra-violet light and is associated with reproductive traits and yearling weight in tropical cattle. Thus, HELB likely contributed to the adaptation of tropical cattle to their harsh environment.

Assessing breed integrity of Göttingen Minipigs

BACKGROUND

Göttingen Minipigs (GMP) is the smallest commercially available minipig breed under a controlled breeding scheme and is globally bred in five isolated colonies. The genetic isolation harbors the risk of stratification which might compromise the identity of the breed and its usability as an animal model for biomedical and human disease. We conducted whole genome re-sequencing of two DNA-pools per colony to assess genomic differentiation within and between colonies. We added publicly available samples from 13 various pig breeds and discovered overall about 32 M loci, ~ 16 M. thereof variable in GMPs. Individual samples were virtually pooled breed-wise. FST between virtual and DNA pools, a phylogenetic tree, principal component analysis (PCA) and evaluation of functional SNP classes were conducted. An F-test was performed to reveal significantly differentiated allele frequencies between colonies. Variation within a colony was quantified as expected heterozygosity.

RESULTS

Phylogeny and PCA showed that the GMP is easily discriminable from all other breads, but that there is also differentiation between the GMP colonies. Dependent on the contrast between GMP colonies, 4 to 8% of all loci had significantly different allele frequencies. Functional annotation revealed that functionally non-neutral loci are less prone to differentiation. Annotation of highly differentiated loci revealed a couple of deleterious mutations in genes with putative effects in the GMPs .

CONCLUSION

Differentiation and annotation results suggest that the underlying mechanisms are rather drift events than directed selection and limited to neutral genome regions. Animal exchange seems not yet necessary. The Relliehausen colony appears to be the genetically most unique GMP sub-population and could be a valuable resource if animal exchange is required to maintain uniformity of the GMP.

Changed Patterns of Genomic Variation Following Recent Domestication: Selection Sweeps in Farmed Atlantic Salmon

The introduction of wild Atlantic salmon into captivity, and their subsequent artificial selection for production traits, has caused phenotypic differences between domesticated fish and their wild counterparts. Identification of regions of the genome underling these changes offers the promise of characterizing the early biological consequences of domestication. In the current study, we sequenced a population of farmed European Atlantic salmon and compared the observed patterns of SNP variation to those found in conspecific wild populations. This identified 139 genomic regions that contained significantly elevated SNP homozygosity in farmed fish when compared to their wild counterparts. The most extreme was adjacent to versican, a gene involved in control of neural crest cell migration. To control for false positive signals, a second and independent dataset of farmed and wild European Atlantic salmon was assessed using the same methodology. A total of 81 outlier regions detected in the first dataset showed significantly reduced homozygosity within the second one, strongly suggesting the genomic regions identified are enriched for true selection sweeps. Examination of the associated genes identified a number previously characterized as targets of selection in other domestic species and that have roles in development, behavior and olfactory system. These include arcvf, sema6, errb4, id2-like, and 6n1-like genes. Finally, we searched for evidence of parallel sweeps using a farmed population of North American origin. This failed to detect a convincing overlap to the putative sweeps present in European populations, suggesting the factors that drive patterns of variation under domestication and early artificial selection were largely independent. This is the first analysis on domestication of aquaculture species exploiting whole-genome sequence data and resulted in the identification of sweeps common to multiple independent populations of farmed European Atlantic salmon.

Was the Cinta Senese Pig Already a Luxury Food in the Late Middle Ages? Ancient DNA and Archaeozoological Evidence from Central Italy

The Cinta senese is a pig breed, highly esteemed for its meat and derived products, characterized by a black coat with a typical white “belt” and documented by scant iconography, since the 13^th–14^th century in Italy. A piece of pottery showing a Cinta pig was found in the Graffignano castle (Northern Latium, Italy) dated 15th–16th centuries, spurring us to investigate the diet of the inhabitants. Ancient DNA analysis was carried out on 21 pig specimens on three nuclear SNPs: (1) g.43597545C>T, on the KIT gene, informative for the identification of the Cinta senese breed; (2) rs81460129, on an intergenic region in chr. 16, which discriminates between domestic pigs and wild boars, and; (3) a SNP on the ZFY/ZFX homologous genes, to determine the sex of the individuals. Our results indicate that the Cinta senese was present in Northern Latium in Late Medieval time, although it was not the only breed, and that pigs, including Cinta, interbred with wild boars, suggesting free-range breeding for all types of pigs. Moreover, the unexpected high proportion of young females may be considered as evidence for the wealth of the family inhabiting the castle.

Genomic diversity, linkage disequilibrium and selection signatures in European local pig breeds assessed with a high density SNP chip

Genetic characterization of local breeds is essential to preserve their genomic variability, to advance conservation policies and to contribute to their promotion and sustainability. Genomic diversity of twenty European local pig breeds and a small sample of Spanish wild pigs was assessed using high density SNP chips. A total of 992 DNA samples were analyzed with the GeneSeek Genomic Profiler (GGP) 70 K HD porcine genotyping chip. Genotype data was employed to compute genetic diversity, population differentiation and structure, genetic distances, linkage disequilibrium and effective population size. Our results point out several breeds, such as Turopolje, Apulo Calabrese, Casertana, Mora Romagnola and Lithuanian indigenous wattle, having the lowest genetic diversity, supported by low heterozygosity and very small effective population size, demonstrating the need of enhanced conservation strategies. Principal components analysis showed the clustering of the individuals of the same breed, with few breeds being clearly isolated from the rest. Several breeds were partially overlapped, suggesting genetic closeness, which was particularly marked in the case of Iberian and Alentejana breeds. Spanish wild boar was also narrowly related to other western populations, in agreement with recurrent admixture between wild and domestic animals. We also searched across the genome for loci under diversifying selection based on F_ST outlier tests. Candidate genes that may underlie differences in adaptation to specific environments and productive systems and phenotypic traits were detected in potentially selected genomic regions.

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.

The first comprehensive description of the expression profile of genes involved in differential body growth and the immune system of the Jeju Native Pig and miniature pig

Sus scrofa provides a major source of animal protein for humans as well as being an excellent biomedical model. This study was carried out to understand, in detail, the genetic and functional variants of Jeju Native Pigs and miniature pigs through differential expression profiling of the genes controlling their immune response, growth performance, and meat quality. The Illumina HiSeq 2000 platform was used for generating 1.3 billion 90 bp paired-end reads, which were mapped to the S. scrofa genome using TopHat2. A total of 2481 and 2768 genes were differentially expressed with 8-log changes in muscle and liver samples, respectively. Five hundred forty-eight genes in muscle and 642 genes in liver samples had BLAST matches within the non-redundant database. GO process and pathway analyses showed enhanced biological processes related to the extracellular structural organization and skeletal muscle cell differentiation in muscle tissue, whereas the liver tissue shares functions related to the inflammatory response. Herein, we identify inflammatory regulatory genes in miniature pigs and growth response genes in Jeju Native Pigs, information which can provide a stronger base for the selection of breeding stock and facilitate further in vitro and in vivo studies for therapeutic purposes.

In silico approaches to discover the functional impact of non-synonymous single nucleotide polymorphisms in selective sweep regions of the Landrace genome

OBJECTIVE

The aim of this study was to discover the functional impact of non-synonymous single nucleotide polymorphisms (nsSNPs) that were found in selective sweep regions of the Landrace genome.

METHODS

Whole-genome re-sequencing data were obtained from 40 pigs, including 14 Landrace, 16 Yorkshire, and 10 wild boars, which were generated with the Illumina HiSeq 2000 platform. The nsSNPs in the selective sweep regions of the Landrace genome were identified, and the impacts of these variations on protein function were predicted to reveal their potential association with traits of the Landrace breed, such as reproductive capacity.

RESULTS

Total of 53,998 nsSNPs in the mapped regions of pigs were identified, and among them, 345 nsSNPs were found in the selective sweep regions of the Landrace genome which were reported previously. The genes featuring these nsSNPs fell into various functional categories, such as reproductive capacity or growth and development during the perinatal period. The impacts of amino acid sequence changes by nsSNPs on protein function were predicted using two in silico SNP prediction algorithms, i.e., sorting intolerant from tolerant and polymorphism phenotyping v2, to reveal their potential roles in biological processes that might be associated with the reproductive capacity of the Landrace breed.

CONCLUSION

The findings elucidated the domestication history of the Landrace breed and illustrated how Landrace domestication led to patterns of genetic variation related to superior reproductive capacity. Our novel findings will help understand the process of Landrace domestication at the genome level and provide SNPs that are informative for breeding.

Contrasting Patterns of Genomic Diversity Reveal Accelerated Genetic Drift but Reduced Directional Selection on X-Chromosome in Wild and Domestic Sheep Species

Analyses of genomic diversity along the X chromosome and of its correlation with autosomal diversity can facilitate understanding of evolutionary forces in shaping sex-linked genomic architecture. Strong selective sweeps and accelerated genetic drift on the X-chromosome have been inferred in primates and other model species, but no such insight has yet been gained in domestic animals compared with their wild relatives. Here, we analyzed X-chromosome variability in a large ovine data set, including a BeadChip array for 943 ewes from the world’s sheep populations and 110 whole genomes of wild and domestic sheep. Analyzing whole-genome sequences, we observed a substantially reduced X-to-autosome diversity ratio (∼0.6) compared with the value expected under a neutral model (0.75). In particular, one large X-linked segment (43.05–79.25 Mb) was found to show extremely low diversity, most likely due to a high density of coding genes, featuring highly conserved regions. In general, we observed higher nucleotide diversity on the autosomes, but a flat diversity gradient in X-linked segments, as a function of increasing distance from the nearest genes, leading to a decreased X: autosome (X/A) diversity ratio and contrasting to the positive correlation detected in primates and other model animals. Our evidence suggests that accelerated genetic drift but reduced directional selection on X chromosome, as well as sex-biased demographic events, explain low X-chromosome diversity in sheep species. The distinct patterns of X-linked and X/A diversity we observed between Middle Eastern and non-Middle Eastern sheep populations can be explained by multiple migrations, selection, and admixture during the domestic sheep’s recent postdomestication demographic expansion, coupled with natural selection for adaptation to new environments. In addition, we identify important novel genes involved in abnormal behavioral phenotypes, metabolism, and immunity, under selection on the sheep X-chromosome.

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.

A genome scan for selection signatures in Taihu pig breeds using next-generation sequencing

Taihu pig breeds are the most prolific breeds of swine in the world, and they also have superior economic traits, including high resistance to disease, superior meat quality, high resistance to crude feed and a docile temperament. The formation of these phenotypic characteristics is largely a result of long-term artificial or natural selection. Therefore, exploring selection signatures in the genomes of the Taihu pigs will help us to identify porcine genes related to productivity traits, disease and behaviour. In this study, we used both intra-population (Relative Extend Haplotype Homozygosity Test (REHH)) and inter-population (the Cross-Population Extend Haplotype Homozygosity Test (XPEHH); F-STATISTICS, F ST ) methods to detect genomic regions that might be under selection process in Taihu pig breeds. As a result, we found 282 (REHH) and 112 (XPEHH) selection signature candidate regions corresponding to 159.78 Mb (6.15%) and 62.29 Mb (2.40%) genomic regions, respectively. Further investigations of the selection candidate regions revealed that many genes under these genomic regions were related to reproductive traits (such as the TLR9 gene), coat colour (such as the KIT gene) and fat metabolism (such as the CPT1A and MAML3 genes). Furthermore, gene enrichment analyses showed that genes under the selection candidate regions were significantly over-represented in pathways related to diseases, such as autoimmune thyroid and asthma diseases. In conclusion, several candidate genes potentially under positive selection were involved in characteristics of Taihu pig. These results will further allow us to better understand the mechanisms of selection in pig breeding.

De novo assembly of mitochondrial genomes provides insights into genetic diversity and molecular evolution in wild boars and domestic pigs

Up to date, the scarcity of publicly available complete mitochondrial sequences for European wild pigs hampers deeper understanding about the genetic changes following domestication. Here, we have assembled 26 de novo mtDNA sequences of European wild boars from next generation sequencing (NGS) data and downloaded 174 complete mtDNA sequences to assess the genetic relationship, nucleotide diversity, and selection. The Bayesian consensus tree reveals the clear divergence between the European and Asian clade and a very small portion (10 out of 200 samples) of maternal introgression. The overall nucleotides diversities of the mtDNA sequences have been reduced following domestication. Interestingly, the selection efficiencies in both European and Asian domestic pigs are reduced, probably caused by changes in both selection constraints and maternal population size following domestication. This study suggests that de novo assembled mitogenomes can be a great boon to uncover the genetic turnover following domestication. Further investigation is warranted to include more samples from the ever-increasing amounts of NGS data to help us to better understand the process of domestication.

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.

Signatures of Selection and Interspecies Introgression in the Genome of Chinese Domestic Pigs

Chinese domestic pigs have experienced strong artificial selection for thousands of years. However, the molecular mechanisms underlying the selection-causing phenotypic changes in Chinese domestic pigs are still largely unknown. Here we used whole-genome resequencing data of 54 pigs from 9 Chinese diverse breeds and 16 wild boars from 7 localities across China to identify genes that show evidence of positive selection in the process of domestication. A total of 14 candidate domestication regions were detected by selective sweep analyses of genetic differentiation and variability, and a set of genes in these candidate domestication regions were found to be related to metabolic process, development, reproduction, olfactory, behavior, and nervous system. The most promising candidate gene under selection - TBX19 - probably underlies the metabolic alteration and developmental traits, and may also associate with timidity of Chinese domestic pigs. Intriguingly, we found that the haplotype at TBX19 locus shared by nearly all Chinese domestic pigs was possibly introgressed from another Sus species. We also revealed the AHR gene associated with female reproduction is under strong positive selection. These results advance our understanding of the evolutionary history of Chinese domestic pigs and shed insights into identifying functionally important genes/mutations contributing to the phenotypic diversity in pigs.

Convergent and divergent genetic changes in the genome of Chinese and European pigs

Since 10,000 BC, continuous human selection has led to intense genetic and phenotypic changes in pig (Sus scrofa) domestication. Through whole genome analysis of 257 individuals, we demonstrated artificial unidirectional and bidirectional selection as the primary force to shape the convergent and divergent changes between Chinese domestic pigs (CHD) and European domestic pigs (EUD). We identified 31 genes in unidirectional selection regions that might be related to fundamental domestication requirements in pigs. And these genes belong predominantly to categories related to the nervous system, muscle development, and especially to metabolic diseases. In addition, 35 genes, representing different breeding preference, were found under bidirectional selection for the distinct leanness and reproduction traits between CHD and EUD. The convergent genetic changes, contributing physical and morphological adaption, represent the common concerns on pig domestication. And the divergent genetic changes reflect distinct breeding goals between Chinese and European pigs. Using ITPR3, AHR and NMU as examples, we explored and validated how the genetic variations contribute to the phenotype changes.

A Pathway-Centered Analysis of Pig Domestication and Breeding in Eurasia

Ascertaining the molecular and physiological basis of domestication and breeding is an active area of research. Due to the current wide distribution of its wild ancestor, the wild boar, the pig (Sus scrofa) is an excellent model to study these processes, which occurred independently in East Asia and Europe ca. 9000 yr ago. Analyzing genome variability patterns in terms of metabolic pathways is attractive since it considers the impact of interrelated functions of genes, in contrast to genome-wide scans that treat genes or genome windows in isolation. To that end, we studied 40 wild boars and 123 domestic pig genomes from Asia and Europe when metabolic pathway was the unit of analysis. We computed statistical significance for differentiation (Fst) and linkage disequilibrium (nSL) statistics at the pathway level. In terms of Fst, we found 21 and 12 pathways significantly differentiated at a q-value < 0.05 in Asia and Europe, respectively; five were shared across continents. In Asia, we found six significant pathways related to behavior, which involved essential neurotransmitters like dopamine and serotonin. Several significant pathways were interrelated and shared a variable percentage of genes. There were 12 genes present in >10 significant pathways (in terms of Fst), comprising genes involved in the transduction of a large number of signals, like phospholipase PCLB1, which is expressed in the brain, or ITPR3, which has an important role in taste transduction. In terms of nSL, significant pathways were mainly related to reproductive performance (ovarian steroidogenesis), a similarly important target trait during domestication and modern animal breeding. Different levels of recombination cannot explain these results, since we found no correlation between Fst and recombination rate. However, we did find an increased ratio of deleterious mutations in domestic vs. wild populations, suggesting a relaxed functional constraint associated with the domestication and breeding processes. Purifying selection was, nevertheless, stronger in significantly differentiated pathways than in random pathways, mainly in Europe. We conclude that pathway analysis facilitates the biological interpretation of genome-wide studies. Notably, in the case of pig, behavior played an important role, among other physiological and developmental processes.

The origin and evolution of fibromelanosis in domesticated chickens: Genomic comparison of Indonesian Cemani and Chinese Silkie breeds

Like Chinese Silkie, Indonesian Ayam Cemani exhibits fibromelanosis or dermal hyperpigmentation and possesses complex segmental duplications on chromosome 20 that involve the endothelin 3 gene, EDN3. A genomic region, DR1 of 127 kb, together with another region, DR2 of 171 kb, was duplicated by unequal crossing over, accompanied by inversion of one DR2. Quantitative PCR and copy number variation analyses on the Cemani genome sequence confirmed the duplication of EDN3. These genetic arrangements are identical in Cemani and Silkie, indicating a single origin of the genetic cause of Fm. The two DR1s harbor two distinct EDN3 haplotypes in a form of permanent heterozygosity, although they remain allelic in the ancestral Red Jungle Fowl population and some domesticated chicken breeds, with their allelic divergence time being as recent as 0.3 million years ago. In Cemani and Silkie breeds, artificial selection favoring the Fm phenotype has left an unambiguous record for selective sweep that extends in both directions from tandemly duplicated EDN3 loci. This highly homozygous tract is different in length between Cemani and Silkie, reflecting their distinct breeding histories. It is estimated that the Fm phenotype came into existence at least 6600-9100 years ago, prior to domestication of Cemani and Silkie, and that throughout domestication there has been intense artificial selection with strength s > 50% in each breed.

Characterization and management of long runs of homozygosity in parental nucleus lines and their associated crossbred progeny

BACKGROUND

In nucleus populations, regions of the genome that have a high frequency of runs of homozygosity (ROH) occur and are associated with a reduction in genetic diversity, as well as adverse effects on fitness. It is currently unclear whether, and to what extent, ROH stretches persist in the crossbred genome and how genomic management in the nucleus population might impact low diversity regions and its implications on the crossbred genome.

METHODS

We calculated a ROH statistic based on lengths of 5 (ROH5) or 10 (ROH10) Mb across the genome for genotyped Landrace (LA), Large White (LW) and Duroc (DU) dams. We simulated crossbred dam (LA × LW) and market [DU × (LA × LW)] animal genotypes based on observed parental genotypes and the ROH frequency was tabulated. We conducted a simulation using observed genotypes to determine the impact of minimizing parental relationships on multiple diversity metrics within nucleus herds, i.e. pedigree-(A), SNP-by-SNP relationship matrix or ROH relationship matrix. Genome-wide metrics included, pedigree inbreeding, heterozygosity and proportion of the genome in ROH of at least 5 Mb. Lastly, the genome was split into bins of increasing ROH5 frequency and, within each bin, heterozygosity, ROH5 and length (Mb) of ROH were evaluated.

RESULTS

We detected regions showing high frequencies of either ROH5 and/or ROH10 across both LW and LA on SSC1, SSC4, and SSC14, and across all breeds on SSC9. Long haplotypes were shared across parental breeds and thus, regions of ROH persisted in crossbred animals. Averaged across replicates and breeds, progeny had higher levels of heterozygosity (0.0056 ± 0.002%) and lower proportion of the genome in a ROH of at least 5 Mb (-0.015 ± 0.003%) than their parental genomes when genomic relationships were constrained, while pedigree relationships resulted in negligible differences at the genomic level. Across all breeds, only genomic data was able to target low diversity regions.

CONCLUSIONS

We show that long stretches of ROH present in the parents persist in crossbred animals. Furthermore, compared to using pedigree relationships, using genomic information to constrain parental relationships resulted in maintaining more genetic diversity and more effectively targeted low diversity regions.

MicroRNA in sperm from Duroc, Landrace and Yorkshire boars

Sperm contain microRNAs (miRNAs), which may have roles in epigenetic control. Regarding phylogenetic relationships among various swine breeds, Yorkshire and Landrace, are considered phenotypically and genetically very similar, but distinctly different from Duroc. The objective of the present study was to compare abundance of boar sperm miRNAs in these three breeds. Overall, 252 prioritized miRNAs were investigated using real-time PCR; relative expression of miRNAs in sperm was similar in Yorkshire and Landrace boars, but significantly different compared to Duroc. Seventeen miRNAs (hsa-miR-196a-5p, hsa-miR-514a-3p, hsa-miR-938, hsa-miR-372-3p, hsa-miR-558, hsa-miR-579-3p, hsa-miR-595, hsa-miR-648, hsa-miR-524-3p, hsa-miR-512-3p, hsa-miR-429, hsa-miR-639, hsa-miR-551a, hsa-miR-624-5p, hsa-miR-585-3p, hsa-miR-508-3p and hsa-miR-626) were down-regulated (P < 0.05; fold regulation ≤-2) in Yorkshire and Landrace sperm, compared to Duroc sperm. Furthermore, three miRNAs (hsa-miR-9-5p, hsa-miR-150-5p, and hsa-miR-99a-5p) were significantly up-regulated in Yorkshire and Landrace sperm compared to Duroc sperm, However, 240 miRNAs were not significantly different (within + 2 fold) between Yorkshire and Landrace sperm. We concluded that miRNAs in sperm were not significantly different between Yorkshire and Landrace boars, but there were significant differences between those two breeds and Duroc boars. Furthermore, integrated target genes for selected down-regulated miRNAs (identified via an in-silico method) appeared to participate in spermatogenesis and sperm functions.

Combined hybridization capture and shotgun sequencing for ancient DNA analysis of extinct wild and domestic dromedary camel

The performance of hybridization capture combined with next‐generation sequencing (NGS) has seen limited investigation with samples from hot and arid regions until now. We applied hybridization capture and shotgun sequencing to recover DNA sequences from bone specimens of ancient‐domestic dromedary (Camelus dromedarius) and its extinct ancestor, the wild dromedary from Jordan, Syria, Turkey and the Arabian Peninsula, respectively. Our results show that hybridization capture increased the percentage of mitochondrial DNA (mtDNA) recovery by an average 187‐fold and in some cases yielded virtually complete mitochondrial (mt) genomes at multifold coverage in a single capture experiment. Furthermore, we tested the effect of hybridization temperature and time by using a touchdown approach on a limited number of samples. We observed no significant difference in the number of unique dromedary mtDNA reads retrieved with the standard capture compared to the touchdown method. In total, we obtained 14 partial mitochondrial genomes from ancient‐domestic dromedaries with 17–95% length coverage and 1.27–47.1‐fold read depths for the covered regions. Using whole‐genome shotgun sequencing, we successfully recovered endogenous dromedary nuclear DNA (nuDNA) from domestic and wild dromedary specimens with 1–1.06‐fold read depths for covered regions. Our results highlight that despite recent methodological advances, obtaining ancient DNA (aDNA) from specimens recovered from hot, arid environments is still problematic. Hybridization protocols require specific optimization, and samples at the limit of DNA preservation need multiple replications of DNA extraction and hybridization capture as has been shown previously for Middle Pleistocene specimens.

Identification of signatures of selection for intramuscular fat and backfat thickness in two Duroc populations

Intramuscular fat (IMF) content is an important trait affecting the quality of pork. Two Duroc populations, one under positive selection for IMF and the other selected for decreased backfat but under stabilizing selection for IMF, were used to identify signatures of selection associated with IMF using 60,000 single-nucleotide polymorphism data. The effects of selection were analyzed between 2 lines or groups representing selected and control animals within each population using a discriminant analysis of principal components and Wright’s fixation index (FST). Moreover, extended haplotype homozygosity-based approaches were used to examine the changes in haplotype frequency due to recent selection. Each statistical method identified 10–20 selection signatures. A few haplotype-based signatures of selection agreed with results from a genome-wide association study (GWAS), while FST measures showed a better agreement with GWAS results. Agreement of marker-trait associations and signatures of selection was limited, and further examination will be necessary to understand the effect of selection on IMF and why some regions identified by GWAS did not appear to respond to the selection practiced. The genes in 21 consensus selection signatures were examined. Several genes with an effect on overall fatness were identified, but further research is needed to assess whether or not some of them could have a specific effect on IMF.

A decade of pig genome sequencing: a window on pig domestication and evolution

Insight into how genomes change and adapt due to selection addresses key questions in evolutionary biology and in domestication of animals and plants by humans. In that regard, the pig and its close relatives found in Africa and Eurasia represent an excellent group of species that enables studies of the effect of both natural and human-mediated selection on the genome. The recent completion of the draft genome sequence of a domestic pig and the development of next-generation sequencing technology during the past decade have created unprecedented possibilities to address these questions in great detail. In this paper, I review recent whole-genome sequencing studies in the pig and closely-related species that provide insight into the demography, admixture and selection of these species and, in particular, how domestication and subsequent selection of Sus scrofa have shaped the genomes of these animals.

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.

The Pig Olfactory Brain: A Primer

Despite the fact that pigs are reputed to have excellent olfactory abilities, few studies have examined regions of the pig brain involved in the sense of smell. The present study provides an overview of the olfactory bulb, anterior olfactory nucleus, and piriform cortex of adult pigs using several approaches. Nissl, myelin, and Golgi stains were used to produce a general overview of the organization of the regions and confocal microscopy was employed to examine 1) projection neurons, 2) GABAergic local circuit neurons that express somatostatin, parvalbumin, vasoactive intestinal polypeptide, or calretinin, 3) neuromodulatory fibers (cholinergic and serotonergic), and 4) glia (astrocytes and microglia). The findings revealed that pig olfactory structures are quite large, highly organized and follow the general patterns observed in mammals.

Evidence for adaptation of porcine Toll-like receptors

Naturally endemic infectious diseases provide selective pressures for pig populations. Toll-like receptors (TLRs) represent the first line of immune defense against pathogens and are likely to play a crucial adaptive role for pig populations. This study was done to determine whether wild and domestic pig populations representing diverse global environments demonstrate local TLR adaptation. The genomic sequence encoding the ectodomain, responsible for interacting with pathogen ligands of bacterial (TLR1, TLR2 and TLR6) and viral (TLR3, TLR7 and TLR8) receptors, was obtained. Mitochondrial D-loop region sequences were obtained and a phylogenetic analysis using these sequences revealed a clear separation of animals into Asian (n = 27) and European (n = 40) clades. The TLR sequences were then analyzed for population-specific positive selection signatures within wild boars and domesticated pig populations derived from Asian and European clades. Using within-population and between-population tests for positive selection, a TLR2-derived variant 376A (126Thr), estimated to have arisen in 163,000 years ago with a frequency of 83.33 % within European wild boars, 98.00 % within domestic pig breeds of European origin, 40.00 % within Asian wild boars, and 11.36 % within Asian domestic pigs, was identified to be under positive selection in pigs of European origin. The variant is located within the N terminal domain of the TLR2 protein 3D crystal structure and could affect ligand binding. This study suggests the TLR2 gene contributing to responses to bacterial pathogens has been crucial in adaptation of pigs to pathogens.

Comparative Analysis of Muscle Transcriptome between Pig Genotypes Identifies Genes and Regulatory Mechanisms Associated to Growth, Fatness and Metabolism

Iberian ham production includes both purebred (IB) and Duroc-crossbred (IBxDU) Iberian pigs, which show important differences in meat quality and production traits, such as muscle growth and fatness. This experiment was conducted to investigate gene expression differences, transcriptional regulation and genetic polymorphisms that could be associated with the observed phenotypic differences between IB and IBxDU pigs. Nine IB and 10 IBxDU pigs were slaughtered at birth. Morphometric measures and blood samples were obtained and samples from Biceps femoris muscle were employed for compositional and transcriptome analysis by RNA-Seq technology. Phenotypic differences were evident at this early age, including greater body size and weight in IBxDU and greater Biceps femoris intramuscular fat and plasma cholesterol content in IB newborns. We detected 149 differentially expressed genes between IB and IBxDU neonates (p < 0.01 and Fold-Change > 1. 5). Several were related to adipose and muscle tissues development (DLK1, FGF21 or UBC). The functional interpretation of the transcriptomic differences revealed enrichment of functions and pathways related to lipid metabolism in IB and to cellular and muscle growth in IBxDU pigs. Protein catabolism, cholesterol biosynthesis and immune system were functions enriched in both genotypes. We identified transcription factors potentially affecting the observed gene expression differences. Some of them have known functions on adipogenesis (CEBPA, EGRs), lipid metabolism (PPARGC1B) and myogenesis (FOXOs, MEF2D, MYOD1), which suggest a key role in the meat quality differences existing between IB and IBxDU hams. We also identified several polymorphisms showing differential segregation between IB and IBxDU pigs. Among them, non-synonymous variants were detected in several transcription factors as PPARGC1B and TRIM63 genes, which could be associated to altered gene function. Taken together, these results provide information about candidate genes, metabolic pathways and genetic polymorphisms potentially involved in phenotypic differences between IB and IBxDU pigs associated to meat quality and production traits.

Exploring transcriptomic diversity in muscle revealed that cellular signaling pathways mainly differentiate five Western porcine breeds

Background

Among transcriptomic studies, those comparing species or populations can increase our understanding of the impact of the evolutionary forces on the differentiation of populations. A particular situation is the one of short evolution time with breeds of a domesticated species that underwent strong selective pressures. In this study, the gene expression diversity across five pig breeds has been explored in muscle. Samples came from: 24 Duroc, 33 Landrace, 41 Large White dam line, 10 Large White sire line and 39 Piétrain. From these animals, 147 muscle samples obtained at slaughter were analyzed using the porcine Agilent 44 K v1 microarray.

Results

A total of 12,358 genes were identified as expressed in muscle after normalization and 1,703 genes were declared differential for at least one breed (FDR < 0.001). The functional analysis highlighted that gene expression diversity is mainly linked to cellular signaling pathways such as the PI3K (phosphoinositide 3-kinase) pathway. The PI3K pathway is known to be involved in the control of development of the skeletal muscle mass by affecting extracellular matrix - receptor interactions, regulation of actin cytoskeleton pathways and some metabolic functions. This study also highlighted 228 spots (171 unique genes) that differentiate the breeds from each other. A common subgroup of 15 genes selected by three statistical methods was able to differentiate Duroc, Large White and Piétrain breeds.

Conclusions

This study on transcriptomic differentiation across Western pig breeds highlighted a global picture: mainly signaling pathways were affected. This result is consistent with the selection objective of increasing muscle mass. These transcriptional changes may indicate selection pressure or simply breed differences which may be driven by human selection. Further work aiming at comparing genetic and transcriptomic diversities would further increase our understanding of the consequences of human impact on livestock species.

Electronic supplementary material

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

Genome sequencing of the extinct Eurasian wild aurochs, Bos primigenius, illuminates the phylogeography and evolution of cattle

Background

Domestication of the now-extinct wild aurochs, Bos primigenius, gave rise to the two major domestic extant cattle taxa, B. taurus and B. indicus. While previous genetic studies have shed some light on the evolutionary relationships between European aurochs and modern cattle, important questions remain unanswered, including the phylogenetic status of aurochs, whether gene flow from aurochs into early domestic populations occurred, and which genomic regions were subject to selection processes during and after domestication. Here, we address these questions using whole-genome sequencing data generated from an approximately 6,750-year-old British aurochs bone and genome sequence data from 81 additional cattle plus genome-wide single nucleotide polymorphism data from a diverse panel of 1,225 modern animals.

Results

Phylogenomic analyses place the aurochs as a distinct outgroup to the domestic B. taurus lineage, supporting the predominant Near Eastern origin of European cattle. Conversely, traditional British and Irish breeds share more genetic variants with this aurochs specimen than other European populations, supporting localized gene flow from aurochs into the ancestors of modern British and Irish cattle, perhaps through purposeful restocking by early herders in Britain. Finally, the functions of genes showing evidence for positive selection in B. taurus are enriched for neurobiology, growth, metabolism and immunobiology, suggesting that these biological processes have been important in the domestication of cattle.

Conclusions

This work provides important new information regarding the origins and functional evolution of modern cattle, revealing that the interface between early European domestic populations and wild aurochs was significantly more complex than previously thought.

Electronic supplementary material

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

A comprehensive mapping of the structure and gene organisation in the sheep MHC class I region

Background

The major histocompatibility complex (MHC) is a chromosomal region that regulates immune responsiveness in vertebrates. This region is one of the most important for disease resistance because it has been associated with resistance or susceptibility to a wide variety of diseases and because the MHC often accounts for more of the variance than other loci. Selective breeding for disease resistance is becoming increasingly common in livestock industries, and it is important to determine how this will influence MHC polymorphism and resistance to diseases that are not targeted for selection. However, in sheep the order and sequence of the protein coding genes is controversial. Yet this information is needed to determine precisely how the MHC influences resistance and susceptibility to disease.

Methods

CHORI bacterial artificial chromosomes (BACs) known to contain sequences from the sheep MHC class I region were sub-cloned, and the clones partially sequenced. The resulting sequences were analysed and re-assembled to identify gene content and organisation within each BAC. The low resolution MHC class I physical map was then compared to the cattle reference genome, the Chinese Merino sheep MHC map published by Gao, et al. (2010) and the recently available sheep reference genome.

Results

Immune related class I genes are clustered into 3 blocks; beta, kappa and a novel block not previously identified in other organisms. The revised map is more similar to Bovidae maps than the previous sheep maps and also includes several genes previously not annotated in the Chinese Merino BAC assembly and others not currently annotated in the sheep reference chromosome 20. In particular, the organisation of nonclassical MHC class I genes is similar to that present in the cattle MHC. Sequence analysis and prediction of amino acid sequences of MHC class I classical and nonclassical genes was performed and it was observed that the map contained one classical and eight nonclassical genes together with three possible pseudogenes.

Conclusions

The comprehensive physical map of the sheep MHC class I region enhances our understanding of the genetic architecture of the class I MHC region in sheep and will facilitate future studies of MHC function.

Electronic supplementary material

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

Detection of Selection Signatures on the X Chromosome in Three Sheep Breeds

Artificial selection has played a critical role in animal breeding. Detection of artificial selection footprints in genomic regions can provide insights for understanding the function of specific phenotypic traits and better guide animal breeding. To more fully understand the relationship between genomic composition and phenotypic diversity arising from breed development, a genome-wide scan was conducted using an OvineSNP50 BeadChip and integrated haplotype score and fixation index analyses to detect selection signatures on the X chromosome in three sheep breeds. We identified 49, 34, and 55 candidate selection regions with lengths of 27.49, 16.47, and 25.42 Mb in German Mutton, Dorper, and Sunit sheep, respectively. Bioinformatics analysis showed that some of the genes in these regions with selection signatures, such as BMP15, were relevant to reproduction. We also identified some selection regions harboring genes that had human orthologs, including BKT, CENPI, GUCY2F, MSN, PCDH11X, PLP1, VSIG4, PAK3, WAS, PCDH19, PDHA1, and SRPX2. The VSIG4 and PCDH11X genes are associated with the immune system and disease, PDHA1 is associated with biosynthetic related pathways, and PCDH19 is expressed in the nervous system and skin. These genes may be useful as candidate genes for molecular breeding.

Defining window-boundaries for genomic analyses using smoothing spline techniques

Background

High-density genomic data is often analyzed by combining information over windows of adjacent markers. Interpretation of data grouped in windows versus at individual locations may increase statistical power, simplify computation, reduce sampling noise, and reduce the total number of tests performed. However, use of adjacent marker information can result in over- or under-smoothing, undesirable window boundary specifications, or highly correlated test statistics. We introduce a method for defining windows based on statistically guided breakpoints in the data, as a foundation for the analysis of multiple adjacent data points. This method involves first fitting a cubic smoothing spline to the data and then identifying the inflection points of the fitted spline, which serve as the boundaries of adjacent windows. This technique does not require prior knowledge of linkage disequilibrium, and therefore can be applied to data collected from individual or pooled sequencing experiments. Moreover, in contrast to existing methods, an arbitrary choice of window size is not necessary, since these are determined empirically and allowed to vary along the genome.

Results

Simulations applying this method were performed to identify selection signatures from pooled sequencing F_ST data, for which allele frequencies were estimated from a pool of individuals. The relative ratio of true to false positives was twice that generated by existing techniques. A comparison of the approach to a previous study that involved pooled sequencing F_ST data from maize suggested that outlying windows were more clearly separated from their neighbors than when using a standard sliding window approach.

Conclusions

We have developed a novel technique to identify window boundaries for subsequent analysis protocols. When applied to selection studies based on F_ST data, this method provides a high discovery rate and minimizes false positives. The method is implemented in the R package GenWin, which is publicly available from CRAN.

A deep catalog of autosomal single nucleotide variation in the pig

A comprehensive catalog of variability in a given species is useful for many important purposes, e.g., designing high density arrays or pinpointing potential mutations of economic or physiological interest. Here we provide a genomewide, worldwide catalog of single nucleotide variants by simultaneously analyzing the shotgun sequence of 128 pigs and five suid outgroups. Despite the high SNP missing rate of some individuals (up to 88%), we retrieved over 48 million high quality variants. Of them, we were able to assess the ancestral allele of more than 39M biallelic SNPs. We found SNPs in 21,455 out of the 25,322 annotated genes in pig assembly 10.2. The annotation showed that more than 40% of the variants were novel variants, not present in dbSNP. Surprisingly, we found a large variability in transition / transversion rate along the genome, which is very well explained (R²=0.79) primarily by genome differences in in CpG content and recombination rate. The number of SNPs per window also varied but was less dependent of known factors such as gene density, missing rate or recombination (R²=0.48). When we divided the samples in four groups, Asian wild boar (ASWB), Asian domestics (ASDM), European wild boar (EUWB) and European domestics (EUDM), we found a marked correlation in allele frequencies between domestics and wild boars within Asia and within Europe, but not across continents, due to the large evolutive distance between pigs of both continents (~1.2 MYA). In general, the porcine species showed a small percentage of SNPs exclusive of each population group. EUWB and EUDM were predicted to harbor a larger fraction of potentially deleterious mutations, according to the SIFT algorithm, than Asian samples, perhaps a result of background selection being less effective due to a lower effective population size in Europe.

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.