Symmetric expression of ohnologs encoding conserved antiviral responses in tetraploid common carp suggest absence of subgenome dominance after whole genome duplication.

Allopolyploids often experience subgenome dominance, with one subgenome showing higher levels of gene expression and greater gene retention. Here, we address the functionality of both subgenomes of allotetraploid common carp (Cyprinus carpio) by analysing a functional network of interferon-stimulated genes (ISGs) crucial in anti-viral immune defence. As an indicator of subgenome dominance we investigated retainment of a core set of ohnologous ISGs. To facilitate our functional genomic analysis a high quality genome was assembled (WagV4.0). Transcriptome data from an in vitro experiment mimicking a viral infection was used to infer ISG expression. Transcriptome analysis confirmed induction of 88 ISG ohnologs on both subgenomes. In both control and infected states, average expression of ISG ohnologs was comparable between the two subgenomes. Also, the highest expressing and most inducible gene copies of an ohnolog pair could be derived from either subgenome. We found no strong evidence of subgenome dominance for common carp.

A natural knockout of the MYO7A gene leads to pre-weaning mortality in pigs.

The pig breeding system provides a unique framework to study recessive defects and the consequence on the phenotype. We examined a commercial synthetic Duroc population for recessive defects and identified a haplotype on chromosome 9 significantly affecting pre-weaning mortality. To identify the causal variant underlying the mortality, we examined sequence data of four carrier animals and 21 non-carrier animals from the same population. The results yield a strong candidate causal stop-gained variant (NM_001099928.1:c.541C>T) affecting the MYO7A gene in complete linkage disequilibrium with the lethal haplotype. The variant leads to an impaired (p.Gln181*) MYO7A protein that truncates 2032 amino acids from the protein. We examined a litter from a carrier sow inseminated by a carrier boar. From the resulting piglets, two confirmed homozygous piglets suffered from severe balance difficulties and the inability to walk properly. The variant segregates at a carrier frequency of 8.2% in the evaluated population and will be gradually purged from the population, improving animal welfare. Finally, this 'natural knockout' will increase our understanding of the functioning of the MYO7A gene and provides a potential model for Usher syndrome in humans.

Genetic origin, admixture and population history of aurochs (Bos primigenius) and primitive European cattle.

This corrects the article DOI: 10.1038/hdy.2016.79.

Fine-mapping of BTA17 using imputed sequences for associations with de novo synthesized fatty acids in bovine milk.

The aim of this study was to fine-map a genomic region associated with milk fatty acids (FA) on Bos taurus autosome (BTA) 17. This genomic region has been discovered with 50,000 (50k) single nucleotide polymorphisms (SNP) imputed to 777,000 (777k) SNP. In this study, high-density genotypes were imputed to whole-genome sequences level to identify candidate gene(s) associated with milk FA composition on BTA17. Phenotypes and genotypes were available for 1,640 cows sampled in winter, and for 1,581 cows sampled in summer. Phenotypes consisted of gas chromatography measurements in winter and in summer milk samples of 6 individual FA and the indicator of de novo synthesis, C6:0-C14:0. Genotypes consisted of imputed 777k SNP, and 89 sequenced ancestors of the population of genotyped cows. In addition, 450 whole-genome sequences from the 1,000 Bull Genome Consortium were available. Using 495 Holstein-Friesian sequences as a reference population, the 777k SNP genotypes of the cows were imputed to sequence level. We then applied single-variant analyses with an animal model, and identified thousands of significant associations with C6:0, C8:0, C10:0, C12:0, C14:0, and C6:0-C14:0. For C8:0 in summer milk samples, the genomic region located between 29 and 34 Mbp on BTA17 revealed a total of 646 significant associations. The most significant associations [-log10(P-value) = 7.82] were 8 SNP in perfect linkage disequilibrium. After fitting one of these 8 SNP as a fixed effect in the model, and re-running the single-variant analyses, no further significant associations were found for any of the 6 FA or C6:0-C14:0. These findings suggest that one polymorphism underlying this QTL on BTA17 influences multiple de novo synthesized milk FA. Thirteen genes in the QTL region were identified and analyzed carefully. Six out of the 8 SNP that showed the strongest associations were located in the La ribonucleoprotein domain family, member 1B (LARP1B) gene, and we suggest LARP1B as a primary candidate gene. Another gene of interest for this QTL region might be PKL4. None of these suggested candidate genes have previously been associated with milk fat synthesis or milk FA composition.

Genetic origin, admixture and population history of aurochs (Bos primigenius) and primitive European cattle.

The domestication of taurine cattle initiated ~10 000 years ago in the Near East from a wild aurochs (Bos primigenius) population followed by their dispersal through migration of agriculturalists to Europe. Although gene flow from wild aurochs still present at the time of this early dispersion is still debated, some of the extant primitive cattle populations are believed to possess the aurochs-like primitive features. In this study, we use genome-wide single nucleotide polymorphisms to assess relationship, admixture patterns and demographic history of an ancient aurochs sample and European cattle populations, several of which have primitive features and are suitable for extensive management. The principal component analysis, the model-based clustering and a distance-based network analysis support previous works suggesting different histories for north-western and southern European cattle. Population admixture analysis indicates a zebu gene flow in the Balkan and Italian Podolic cattle populations. Our analysis supports the previous report of gene flow between British and Irish primitive cattle populations and local aurochs. In addition, we show evidence of aurochs gene flow in the Iberian cattle populations indicating wide geographical distribution of the aurochs. Runs of homozygosity (ROH) reveal that demographic processes like genetic isolation and breed formation have contributed to genomic variations of European cattle populations. The ROH also indicate recent inbreeding in southern European cattle populations. We conclude that in addition to factors such as ancient human migrations, isolation by distance and cross-breeding, gene flow between domestic and wild-cattle populations also has shaped genomic composition of European cattle populations.

Accuracy of genomic prediction using imputed whole-genome sequence data in white layers.

There is an increasing interest in using whole-genome sequence data in genomic selection breeding programmes. Prediction of breeding values is expected to be more accurate when whole-genome sequence is used, because the causal mutations are assumed to be in the data. We performed genomic prediction for the number of eggs in white layers using imputed whole-genome resequence data including ~4.6 million SNPs. The prediction accuracies based on sequence data were compared with the accuracies from the 60 K SNP panel. Predictions were based on genomic best linear unbiased prediction (GBLUP) as well as a Bayesian variable selection model (BayesC). Moreover, the prediction accuracy from using different types of variants (synonymous, non-synonymous and non-coding SNPs) was evaluated. Genomic prediction using the 60 K SNP panel resulted in a prediction accuracy of 0.74 when GBLUP was applied. With sequence data, there was a small increase (~1%) in prediction accuracy over the 60 K genotypes. With both 60 K SNP panel and sequence data, GBLUP slightly outperformed BayesC in predicting the breeding values. Selection of SNPs more likely to affect the phenotype (i.e. non-synonymous SNPs) did not improve the accuracy of genomic prediction. The fact that sequence data were based on imputation from a small number of sequenced animals may have limited the potential to improve the prediction accuracy. A small reference population (n = 1004) and possible exclusion of many causal SNPs during quality control can be other possible reasons for limited benefit of sequence data. We expect, however, that the limited improvement is because the 60 K SNP panel was already sufficiently dense to accurately determine the relationships between animals in our data.

Genomic diversity and differentiation of a managed island wild boar population.

The evolution of island populations in natural systems is driven by local adaptation and genetic drift. However, evolutionary pathways may be altered by humans in several ways. The wild boar (WB) (Sus scrofa) is an iconic game species occurring in several islands, where it has been strongly managed since prehistoric times. We examined genomic diversity at 49 803 single-nucleotide polymorphisms in 99 Sardinian WBs and compared them with 196 wild specimens from mainland Europe and 105 domestic pigs (DP; 11 breeds). High levels of genetic variation were observed in Sardinia (80.9% of the total number of polymorphisms), which can be only in part associated to recent genetic introgression. Both Principal Component Analysis and Bayesian clustering approach revealed that the Sardinian WB population is highly differentiated from the other European populations (FST=0.126-0.138), and from DP (FST=0.169). Such evidences were mostly unaffected by an uneven sample size, although clustering results in reference populations changed when the number of individuals was standardized. Runs of homozygosity (ROHs) pattern and distribution in Sardinian WB are consistent with a past expansion following a bottleneck (small ROHs) and recent population substructuring (highly homozygous individuals). The observed effect of a non-random selection of Sardinian individuals on diversity, FST and ROH estimates, stressed the importance of sampling design in the study of structured or introgressed populations. Our results support the heterogeneity and distinctiveness of the Sardinian population and prompt further investigations on its origins and conservation status.

Shallow genetic divergence and species delineations in the endemic Labeobarbus species flock of Lake Tana, Ethiopia.

To assess whether the species distinctions of Lake Tana's Labeobarbus spp. are supported by genetic information, microsatellite markers were used. A total of 376 Labeobarbus spp., belonging to 24 populations of 11 species from three regions of the lake (north, south and east), were sampled. Eight microsatellite markers were analysed. In general, differences between conspecific populations were smaller than differences between populations of different species. For six species, conspecific populations from different regions in the lake were consistently more similar than populations of other species from the same region. For four species this was not the case, while for one species two populations were similar, but different from the third population. River-spawning species appeared to be more distinct than presumed lake spawners. On the species level, there was a significant correlation between genetic and morphological differentiation, especially in morphological aspects associated with ecological functioning. This suggests that genetic differentiation arose together with adaptive radiation, although the overall genetic differentiation among the Lake Tana Labeobarbus spp. is small.

Systematic differences in the response of genetic variation to pedigree and genome-based selection methods.

Genomic selection (GS) is a DNA-based method of selecting for quantitative traits in animal and plant breeding, and offers a potentially superior alternative to traditional breeding methods that rely on pedigree and phenotype information. Using a 60 K SNP chip with markers spaced throughout the entire chicken genome, we compared the impact of GS and traditional BLUP (best linear unbiased prediction) selection methods applied side-by-side in three different lines of egg-laying chickens. Differences were demonstrated between methods, both at the level and genomic distribution of allele frequency changes. In all three lines, the average allele frequency changes were larger with GS, 0.056 0.064 and 0.066, compared with BLUP, 0.044, 0.045 and 0.036 for lines B1, B2 and W1, respectively. With BLUP, 35 selected regions (empirical P < 0.05) were identified across the three lines. With GS, 70 selected regions were identified. Empirical thresholds for local allele frequency changes were determined from gene dropping, and differed considerably between GS (0.167-0.198) and BLUP (0.105-0.126). Between lines, the genomic regions with large changes in allele frequencies showed limited overlap. Our results show that GS applies selection pressure much more locally than BLUP, resulting in larger allele frequency changes. With these results, novel insights into the nature of selection on quantitative traits have been gained and important questions regarding the long-term impact of GS are raised. The rapid changes to a part of the genetic architecture, while another part may not be selected, at least in the short term, require careful consideration, especially when selection occurs before phenotypes are observed.

Farm-by-farm analysis of microsatellite, mtDNA and SNP genotype data reveals inbreeding and crossbreeding as threats to the survival of a native Spanish pig breed.

The Chato Murciano (CM), a pig breed from the Murcia region in the southeastern region of Spain, is a good model for endangered livestock populations. The remaining populations are bred on approximately 15 small farms, and no herdbook exists. To assess the genetic threats to the integrity and survival of the CM breed, and to aid in designing a conservation program, three genetic marker systems - microsatellites, SNPs and mtDNA - were applied across the majority of the total breeding stock. In addition, mtDNA and SNPs were genotyped in breeds that likely contributed genetically to the current CM gene pool. The analyses revealed the levels of genetic diversity within the range of other European local breeds (H(e) = 0.53). However, when the eight farms that rear at least 10 CM pigs were independently analyzed, high levels of inbreeding were found in some. Despite the evidence for recent crossbreeding with commercial breeds on a few farms, the entire breeding stock remains readily identifiable as CM, facilitating the design of traceability assays. The genetic management of the breed is consistent with farm size, farm owner and presence of other pig breeds on the farm, demonstrating the highly ad hoc nature of current CM breeding. The results of genetic diversity and substructure of the entire breed, as well as admixture and crossbreeding obtained in the present study, provide a benchmark to develop future conservation strategies. Furthermore, this study demonstrates that identifying farm-based practices and farm-based breeding stocks can aid in the design of a sustainable breeding program for minority breeds.

Genome-wide single nucleotide polymorphism analysis reveals recent genetic introgression from domestic pigs into Northwest European wild boar populations.

Present-day genetic introgression from domestic pigs into European wild boar has been suggested in various studies. However, no hybrids have been identified beyond doubt mainly because available methods were unable to quantify the extent of introgression and rule out natural processes. Genetic introgression from domestic pigs may have far-reaching ecological consequences by altering traits like the reproduction rate or immunology of wild boar. In this study, we demonstrate a novel approach to investigate genetic introgression in a Northwest (NW) European wild boar data set using a genome-wide single nucleotide polymorphism (SNP) assay developed for domestic pigs. We quantified the extent of introgression using allele frequency spectrum analysis, in silico hybridization simulations and genome distribution patterns of introgressed SNPs. Levels of recent introgression in the study area were expected to be low, as pig farming practices are prevailingly intensive and indoors. However, evidence was found for geographically widespread presence of domestic pig SNPs in 10% of analysed wild boar. This was supported by the identification of two different pig mitochondrial DNA haplotypes in three of the identified hybrid wild boar, suggesting that introgression had occurred from multiple sources (pig breeds). In silico hybridization simulations showed that the level of introgression in the identified hybrid wild boar is equivalent to first-generation hybrids until fifth-generation backcrosses with wild boar. The distribution pattern of introgressed SNPs supported these assignments in four of nine hybrids. The other five hybrids are considered advanced-generation hybrids, resulting from interbreeding among hybrid individuals. Three of nine hybrids were genetically associated with a different wild boar population than the one in which they were sampled. This discrepancy suggests that genetic introgression has occurred through the escape or release of an already hybridized farmed wild boar stock. We conclude that genetic introgression from domestic pigs into NW European wild boar populations is more recent and more common than expected and that genome-wide SNP analysis is a promising tool to quantify recent hybridization in free-living populations.

Porcine colonization of the Americas: a 60k SNP story.

The pig, Sus scrofa, is a foreign species to the American continent. Although pigs originally introduced in the Americas should be related to those from the Iberian Peninsula and Canary islands, the phylogeny of current creole pigs that now populate the continent is likely to be very complex. Because of the extreme climates that America harbors, these populations also provide a unique example of a fast evolutionary phenomenon of adaptation. Here, we provide a genome wide study of these issues by genotyping, with a 60k SNP chip, 206 village pigs sampled across 14 countries and 183 pigs from outgroup breeds that are potential founders of the American populations, including wild boar, Iberian, international and Chinese breeds. Results show that American village pigs are primarily of European ancestry, although the observed genetic landscape is that of a complex conglomerate. There was no correlation between genetic and geographical distances, neither continent wide nor when analyzing specific areas. Most populations showed a clear admixed structure where the Iberian pig was not necessarily the main component, illustrating how international breeds, but also Chinese pigs, have contributed to extant genetic composition of American village pigs. We also observe that many genes related to the cardiovascular system show an increased differentiation between altiplano and genetically related pigs living near sea level.

Identification of high utility SNPs for population assignment and traceability purposes in the pig using high-throughput sequencing.

The objectives of this study were to develop breed-specific single nucleotide polymorphisms (SNPs) in five pig breeds sequenced with Illumina's Genome Analyzer and to investigate their usefulness for breed assignment purposes. DNA pools were prepared for Duroc, Landrace, Large White, Pietrain and Wild Boar. The total number of animals used for sequencing was 153. SNP discovery was performed by aligning the filtered reads against Build 7 of the pig genome. A total of 313,964 high confidence SNPs were identified and analysed for the presence of breed-specific SNPs (defined in this context as SNPs for which one of the alleles was detected in only one breed). There were 29,146 putative breed-specific SNPs identified, of which 4441 were included in the PorcineSNP60 beadchip. Upon re-examining the genotypes obtained using the beadchip, 193 SNPs were confirmed as being breed specific. These 193 SNPs were subsequently used to assign an additional 490 individuals from the same breeds, using the sequenced individuals as reference populations. In total, four breed assignment tests were performed. Results showed that for all methods tested 99% of the animals were correctly assigned, with an average probability of assignment of at least 99.2%, indicating the high utility of breed-specific markers for breed assignment and traceability. This study provides a blueprint for the way next-generation sequencing technologies can be used for the identification of breed-specific SNPs, as well as evidence that these SNPs may be a powerful tool for breed assignment and traceability of animal products to their breeds of origin.

East Asian contributions to Dutch traditional and western commercial chickens inferred from mtDNA analysis.

Understanding the complex origin of domesticated populations is of vital importance for understanding, preserving and exploiting breed genetic diversity. Here, we aim to assess Asian contributions to European traditional breeds and western commercial chickens for mitochondrial genetic diversity. To this end, a 365-bp fragment of the chicken mtDNA D-loop region of 16 Dutch fancy breeds (113 individuals) was surveyed, comprising almost the entire breed diversity of The Netherlands. We also sequenced the same fragment for 160 commercial birds representing all important commercial types from multiple commercial companies that together represent more than 50% of the worldwide commercial value. We identified 20 different haplotypes. The haplotypes clustered into five clades. The commonest clade (E-clade) supposedly originates from the Indian subcontinent. In addition, both in commercial chicken and Dutch fancy breeds, many haplotypes were found with a clear East Asian origin. However, the erratic occurrence of many different East Asian mitochondrial clades indicates that there were many independent instances where breeders used imported exotic chickens for enhancing local breeds. Nucleotide diversity and haplotype diversity analyses showed the influence of the introgression of East Asian chicken on genetic diversity. All populations that had haplotypes of multiple origin displayed high inferred diversity, as opposed to most populations that had only a single mitochondrial haplotype signature. Most fancy breeds were found to have a much lower within-population diversity compared to broilers and layers, although this is not the case for mitochondrial estimates in fancy breeds that have multiple origin haplotypes.

Segregation distortion in chicken and the evolutionary consequences of female meiotic drive in birds.

As all four meiotic products give rise to sperm in males, female meiosis result in a single egg in most eukaryotes. Any genetic element with the potential to influence chromosome segregation, so that it is preferentially included in the egg, should therefore gain a transmission advantage; a process termed female meiotic drive. We are aware of two chromosomal components, centromeres and telomeres, which share the potential to influence chromosome movement during meioses and make the following predictions based on the presence of female meiotic drive: (1) centromere-binding proteins should experience rapid evolution as a result of a conflict between driving centromeres and the rest of the genome; and (2) segregation patterns should be skewed near centromeres and telomeres. To test these predictions, we first analyze the molecular evolution of seven centromere-binding proteins in nine divergent bird species. We find strong evidence for positive selection in two genes, lending support to the genomic conflict hypothesis. Then, to directly test for the presence of segregation distortion, we also investigate the transmission of approximately 9000 single-nucleotide polymorphisms in 197 chicken families. By simulating fair Mendelian meioses, we locate chromosomal regions with statistically significant transmission ratio distortion. One region is located near the centromere on chromosome 1 and a second region is located near the telomere on the p-arm of chromosome 1. Although these observations do not provide conclusive evidence in favour of the meiotic drive/genome conflict hypothesis, they do lend support to the hypothesis that centromeres and telomeres drive during female meioses in chicken.

Extent of linkage disequilibrium in chicken.

Many of the economically important traits in chicken are multifactorial and governed by multiple genes located at different quantitative trait loci (QTLs). The optimal marker density to identify these QTLs in linkage and association studies is largely determined by the extent of linkage disequilibrium (LD) around them. In this study, we investigated the extent of LD on two chromosomes in a white layer and two broiler chicken breeds. Pairwise levels of LD were calculated for 33 and 36 markers on chromosomes 10 and 28, respectively. We found that useful LD (i.e. an r(2) value higher than 0.3) in Nutreco chicken breed E5 (inbred) can extend to around 1 cM on chromosomes 10 and 28, although in a second region on chromosome 28 it extends to about 2.5 cM. The extent in breed Nutreco E3 (outbred) was very short in chromosome 10 (15 kb) but very much larger on chromosome 28, particularly in one region of depressed heterozygosity. The layer breed E2 (inbred) showed an extent of useful LD up to 4 cM on chromosome 10; the extent on chromosome 28 could not be assessed due to an erratic pattern of LD on that chromosome, although in one region LD appears to be in the order of 0.8 cM. This indicates that there may be very large differences in patterns of LD between different chicken breeds and different genomic regions.

Rapid evolution of the MH class I locus results in different allelic compositions in recently diverged populations of Atlantic salmon.

We compared major histocompatibility class I allelic diversity in two currently reproductively isolated Atlantic salmon (Salmo salar) populations (Irish and Norwegian) with a common postglacial origin in order to test for among-population differences in allelic composition and patterns of recombination and point mutation. We also examined the evidence for adaptive molecular divergence at this locus by analyzing the rate of amino acid replacement in relation to a neutral expectation. Contrary to our prediction, and in contrast to the situation for other genetic markers, the two populations have almost nonoverlapping sets of major histocompatibility class I alleles. Although there is a strong signal of point mutation that predates population divergence, recent recombination, acting in similar, but not identical, ways in both populations appears to be a significant force in creating new alleles. Moreover, selection acting on peptide-binding residues seems to favor new recombinant alleles and is likely to be responsible for the rapid divergence between populations.

Patterns of variability at the major histocompatibility class II alpha locus in Atlantic salmon contrast with those at the class I locus.

In order to investigate the mechanisms creating and maintaining variability at the major histocompatibility (MH) class II alpha (DAA) locus we examined patterns of polymorphism in two isolated Atlantic salmon populations which share a common post-glacial origin. As expected from their common origin, but contrary to the observation at the MH class I locus, these populations shared the majority of DAA alleles: out of 17 sequences observed, 11 were common to both populations. Recombination seems to play a more important role in the origin of new alleles at the class II alpha locus than at the class I locus. A greater than expected proportion of sites inferred to be positively selected (potentially peptide binding residues, PBRs) were found to be involved in recombination events, suggesting a mechanism for increasing MH variability through an interaction between recombination and natural selection. Thus it appears that although selection and recombination are important mechanisms for the evolution of both class II alpha and class I loci in the Atlantic salmon, the pattern of variability differs markedly between these classes of MH loci.

Urinary metabolites and haemoglobin adducts as biomarkers of exposure to 1,3-butadiene: a basis for 1,3-butadiene cancer risk assessment.

Since 1,3-butadiene (BD) is a suspected human carcinogen, exposure to BD should be minimised and controlled. This study aimed at comparing the suitability of biomarkers for low levels of exposure to BD, and at exploration of the relative pathways of human metabolism of BD for comparison with experimental animals. Potentially sensitive biomarkers for BD are its urinary metabolites 1,2-dihydroxybutyl mercapturic acid (DHBMA, also referred to as MI) and 1- and 2-monohydroxy-3-butenyl mercapturic acid (MHBMA, also referred to as MII) and its haemoglobin (Hb) adducts 1- and 2-hydroxy-3-butenyl valine (MHBVal). In two field studies in BD-workers, airborne BD, MHBMA, DHBMA and MHBVal were determined. MHBMA proved more sensitive than DHBMA for monitoring recent exposures to BD and could measure 8-h time weighted average exposures as low as 0.13 ppm (0.29 mg/m(3)). The sensitivity of DHBMA was restricted by relatively high natural background levels in urine, of which the origin is currently unknown. MHBVal proved a sensitive method for monitoring cumulative exposures to BD at or above 0.35 ppm (0.77 mg/m(3)). Statistically significant relationships were found between either MHBMA or DHBMA and 8-h airborne BD levels, and between MHBVal adducts and average airborne BD levels over 60 days. The data showed a much higher rate of hydrolytic metabolism of BD in humans compared to animals, which was reflected in a much higher DHBMA/(MHBMA+DHBMA) ratio, and in much lower levels of MHBVal in humans, confirming in vitro results. Assuming a genotoxic mechanism, the data of this study coupled with our recent data on DNA and Hb binding in rodents, suggest that the cancer risk for humans from exposure to BD will be less than for the rat, and much less than for the mouse.