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.

Identification and characterisation of de novo germline structural variants in two commercial pig lines using trio-based whole genome sequencing.

BACKGROUND: De novo mutations arising in the germline are a source of genetic variation and their discovery broadens our understanding of genetic disorders and evolutionary patterns. Although the number of de novo single nucleotide variants (dnSNVs) has been studied in a number of species, relatively little is known about the occurrence of de novo structural variants (dnSVs). In this study, we investigated 37 deeply sequenced pig trios from two commercial lines to identify dnSVs present in the offspring. The identified dnSVs were characterised by identifying their parent of origin, their functional annotations and characterizing sequence homology at the breakpoints. RESULTS: We identified four swine germline dnSVs, all located in intronic regions of protein-coding genes. Our conservative, first estimate of the swine germline dnSV rate is 0.108 (95% CI 0.038-0.255) per generation (one dnSV per nine offspring), detected using short-read sequencing. Two detected dnSVs are clusters of mutations. Mutation cluster 1 contains a de novo duplication, a dnSNV and a de novo deletion. Mutation cluster 2 contains a de novo deletion and three de novo duplications, of which one is inverted. Mutation cluster 2 is 25 kb in size, whereas mutation cluster 1 (197 bp) and the other two individual dnSVs (64 and 573 bp) are smaller. Only mutation cluster 2 could be phased and is located on the paternal haplotype. Mutation cluster 2 originates from both micro-homology as well as non-homology mutation mechanisms, where mutation cluster 1 and the other two dnSVs are caused by mutation mechanisms lacking sequence homology. The 64 bp deletion and mutation cluster 1 were validated through PCR. Lastly, the 64 bp deletion and the 573 bp duplication were validated in sequenced offspring of probands with three generations of sequence data. CONCLUSIONS: Our estimate of 0.108 dnSVs per generation in the swine germline is conservative, due to our small sample size and restricted possibilities of dnSV detection from short-read sequencing. The current study highlights the complexity of dnSVs and shows the potential of breeding programs for pigs and livestock species in general, to provide a suitable population structure for identification and characterisation of dnSVs.

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.

A high-density SNP chip for genotyping great tit (Parus major) populations and its application to studying the genetic architecture of exploration behaviour.

High-density SNP microarrays ("SNP chips") are a rapid, accurate and efficient method for genotyping several hundred thousand polymorphisms in large numbers of individuals. While SNP chips are routinely used in human genetics and in animal and plant breeding, they are less widely used in evolutionary and ecological research. In this article, we describe the development and application of a high-density Affymetrix Axiom chip with around 500,000 SNPs, designed to perform genomics studies of great tit (Parus major) populations. We demonstrate that the per-SNP genotype error rate is well below 1% and that the chip can also be used to identify structural or copy number variation. The chip is used to explore the genetic architecture of exploration behaviour (EB), a personality trait that has been widely studied in great tits and other species. No SNPs reached genomewide significance, including at DRD4, a candidate gene. However, EB is heritable and appears to have a polygenic architecture. Researchers developing similar SNP chips may note: (i) SNPs previously typed on alternative platforms are more likely to be converted to working assays; (ii) detecting SNPs by more than one pipeline, and in independent data sets, ensures a high proportion of working assays; (iii) allele frequency ascertainment bias is minimized by performing SNP discovery in individuals from multiple populations; and (iv) samples with the lowest call rates tend to also have the greatest genotyping error rates.

Genome-wide population structure and admixture analysis reveals weak differentiation among Ugandan goat breeds.

Uganda has a large population of goats, predominantly from indigenous breeds reared in diverse production systems, whose existence is threatened by crossbreeding with exotic Boer goats. Knowledge about the genetic characteristics and relationships among these Ugandan goat breeds and the potential admixture with Boer goats is still limited. Using a medium-density single nucleotide polymorphism (SNP) panel, we assessed the genetic diversity, population structure and admixture in six goat breeds in Uganda: Boer, Karamojong, Kigezi, Mubende, Small East African and Sebei. All the animals had genotypes for about 46 105 SNPs after quality control. We found high proportions of polymorphic SNPs ranging from 0.885 (Kigezi) to 0.928 (Sebei). The overall mean observed (HO ) and expected (HE ) heterozygosity across breeds was 0.355 ± 0.147 and 0.384 ± 0.143 respectively. Principal components, genetic distances and admixture analyses revealed weak population sub-structuring among the breeds. Principal components separated Kigezi and weakly Small East African from other indigenous goats. Sebei and Karamojong were tightly entangled together, whereas Mubende occupied a more central position with high admixture from all other local breeds. The Boer breed showed a unique cluster from the Ugandan indigenous goat breeds. The results reflect common ancestry but also some level of geographical differentiation. admixture and f4 statistics revealed gene flow from Boer and varying levels of genetic admixture among the breeds. Generally, moderate to high levels of genetic variability were observed. Our findings provide useful insights into maintaining genetic diversity and designing appropriate breeding programs to exploit within-breed diversity and heterozygote advantage in crossbreeding schemes.

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

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

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.

After genome-wide association studies: Gene networks elucidating candidate genes divergences for number of teats across two pig populations.

Number of teats (NT) is an important trait affecting both piglet's welfare and the production level of pig farms. Biologically, embryonic mammary gland development requires the coordination of many signaling pathways necessary for the proper development of teats. Several QTL for NT have been identified; however, further analysis is still lacking. Therefore, gene networks derived from genomewide association study (GWAS) results can be used to examine shared pathways and functions of putative candidate genes. Besides, such analyses may also be helpful to understand the genetic diversity between populations for the same trait or traits. In this study, we identified significant SNP for Landrace-based (line C) and Large White-based (line D) dam lines. Besides, gene-transcription factor (TF) networks were constructed aiming to obtain the most likely candidate genes for NT in each line followed by a comparative analysis between both lines to access similarities or dissimilarities at the marker and gene level. We identified 24 and 19 significant SNP (Bayes factor ≥ 100) for lines C and D, respectively. Only 1 significant SNP overlapped both lines. Network analysis illustrated gene interactions consistent with known mammal's breast biology and captured known TF. We observed different sets of putative candidate genes for NT in each line evaluated that may have common effects on the phenotype. Based on these results, we demonstrated the importance of post-GWAS analyses increasing the biological understanding of relevant genes for a complex trait. Moreover, we believe that this genomic diversity across lines should be taken into account, considering breed-specific reference populations for genomic selection.

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.

Accuracy of genomic prediction using deregressed breeding values estimated from purebred and crossbred offspring phenotypes in pigs.

Genomic selection is applied to dairy cattle breeding to improve the genetic progress of purebred (PB) animals, whereas in pigs and poultry the target is a crossbred (CB) animal for which a different strategy appears to be needed. The source of information used to estimate the breeding values, i.e., using phenotypes of CB or PB animals, may affect the accuracy of prediction. The objective of our study was to assess the direct genomic value (DGV) accuracy of CB and PB pigs using different sources of phenotypic information. Data used were from 3 populations: 2,078 Dutch Landrace-based, 2,301 Large White-based, and 497 crossbreds from an F1 cross between the 2 lines. Two female reproduction traits were analyzed: gestation length (GLE) and total number of piglets born (TNB). Phenotypes used in the analyses originated from offspring of genotyped individuals. Phenotypes collected on CB and PB animals were analyzed as separate traits using a single-trait model. Breeding values were estimated separately for each trait in a pedigree BLUP analysis and subsequently deregressed. Deregressed EBV for each trait originating from different sources (CB or PB offspring) were used to study the accuracy of genomic prediction. Accuracy of prediction was computed as the correlation between DGV and the DEBV of the validation population. Accuracy of prediction within PB populations ranged from 0.43 to 0.62 across GLE and TNB. Accuracies to predict genetic merit of CB animals with one PB population in the training set ranged from 0.12 to 0.28, with the exception of using the CB offspring phenotype of the Dutch Landrace that resulted in an accuracy estimate around 0 for both traits. Accuracies to predict genetic merit of CB animals with both parental PB populations in the training set ranged from 0.17 to 0.30. We conclude that prediction within population and trait had good predictive ability regardless of the trait being the PB or CB performance, whereas using PB population(s) to predict genetic merit of CB animals had zero to moderate predictive ability. We observed that the DGV accuracy of CB animals when training on PB data was greater than or equal to training on CB data. However, when results are corrected for the different levels of reliabilities in the PB and CB training data, we showed that training on CB data does outperform PB data for the prediction of CB genetic merit, indicating that more CB animals should be phenotyped to increase the reliability and, consequently, accuracy of DGV for CB genetic merit.

Asian low-androstenone haplotype on pig chromosome 6 does not unfavorably affect production and reproduction traits.

European pigs that carry Asian haplotypes of a 1.94-Mbp region on pig chromosome 6 have lower levels of androstenone, one of the two main compounds causing boar taint. The objective of our study was to examine potential pleiotropic effects of the Asian low-androstenone haplotypes. A single nucleotide polymorphism marker, rs81308021, distinguishes the Asian from European haplotypes and was used to investigate possible associations of androstenone with production and reproduction traits. Eight traits were available from three European commercial breeds. For the two sow lines studied, a favorable effect on number of teats was detected for the low-androstenone haplotype. In one of these sow lines, a favorable effect on number of spermatozoa per ejaculation was detected for the low-androstenone haplotype. No unfavorable pleiotropic effects were found, which suggests that selection for low-androstenone haplotypes within the 1.94 Mbp would not unfavorably affect the other eight relevant traits.

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.

Identification of species-specific novel transcripts in pig reproductive tissues using RNA-seq.

Although structural properties of the porcine reproductive system are shared by many placental mammals, some combination of these properties is unique to pigs. To explore whether genomic elements specific to pigs could potentially underlie this uniqueness, we made the first step to identify novel transcripts in two representative pig reproductive tissues by the technique of massively parallel sequencing. To automate the whole process, we built a computational pipeline, which can also be easily extended for similar studies in other species. In total, 5516 and 9061 novel transcripts were found, and 159 and 252 novel transcripts appear to be specific to pigs for the placenta and testis respectively. Furthermore, these novel transcripts were found to be enriched in quantitative trait loci (QTL) regions for reproduction traits in pigs. We validated eight of these novel transcripts by quantitative real-time PCR. With respect to their genomic organization and their functional relationship to reproduction, these transcripts need to be further validated and explored in various pig breeds to better comprehend the relevant aspects of pig physiology that contribute to reproductive performance.

Replicated high-density genetic maps of two great tit populations reveal fine-scale genomic departures from sex-equal recombination rates.

Linking variation in quantitative traits to variation in the genome is an important, but challenging task in the study of life-history evolution. Linkage maps provide a valuable tool for the unravelling of such trait-gene associations. Moreover, they give insight into recombination landscapes and between-species karyotype evolution. Here we used genotype data, generated from a 10k single-nucleotide polymorphism (SNP) chip, of over 2000 individuals to produce high-density linkage maps of the great tit (Parus major), a passerine bird that serves as a model species for ecological and evolutionary questions. We created independent maps from two distinct populations: a captive F2-cross from The Netherlands (NL) and a wild population from the United Kingdom (UK). The two maps contained 6554 SNPs in 32 linkage groups, spanning 2010 cM and 1917 cM for the NL and UK populations, respectively, and were similar in size and marker order. Subtle levels of heterochiasmy within and between chromosomes were remarkably consistent between the populations, suggesting that the local departures from sex-equal recombination rates have evolved. This key and surprising result would have been impossible to detect if only one population was mapped. A comparison with zebra finch Taeniopygia guttata, chicken Gallus gallus and the green anole lizard Anolis carolinensis genomes provided further insight into the evolution of avian karyotypes.

Genetic consequences of breaking migratory traditions in barnacle geese Branta leucopsis.

Cultural transmission of migratory traditions enables species to deal with their environment based on experiences from earlier generations. Also, it allows a more adequate and rapid response to rapidly changing environments. When individuals break with their migratory traditions, new population structures can emerge that may affect gene flow. Recently, the migratory traditions of the Barnacle Goose Branta leucopsis changed, and new populations differing in migratory distance emerged. Here, we investigate the population genetic structure of the Barnacle Goose to evaluate the consequences of altered migratory traditions. We used a set of 358 single nucleotide polymorphism (SNP) markers to genotype 418 individuals from breeding populations in Greenland, Spitsbergen, Russia, Sweden and the Netherlands, the latter two being newly emerged populations. We used discriminant analysis of principal components, FST , linkage disequilibrium and a comparison of geneflow models using migrate-n to show that there is significant population structure, but that relatively many pairs of SNPs are in linkage disequilibrium, suggesting recent admixture between these populations. Despite the assumed traditions of migration within populations, we also show that genetic exchange occurs between all populations. The newly established nonmigratory population in the Netherlands is characterized by high emigration into other populations, which suggests more exploratory behaviour, possibly as a result of shortened parental care. These results suggest that migratory traditions in populations are subject to change in geese and that such changes have population genetic consequences. We argue that the emergence of nonmigration probably resulted from developmental plasticity.

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.

SNP marker detection and genotyping in tilapia.

We have generated a unique resource consisting of nearly 175 000 short contig sequences and 3569 SNP markers from the widely cultured GIFT (Genetically Improved Farmed Tilapia) strain of Nile tilapia (Oreochromis niloticus). In total, 384 SNPs were selected to monitor the wider applicability of the SNPs by genotyping tilapia individuals from different strains and different geographical locations. In all strains and species tested (O. niloticus, O. aureus and O. mossambicus), the genotyping assay was working for a similar number of SNPs (288-305 SNPs). The actual number of polymorphic SNPs was, as expected, highest for individuals from the GIFT population (255 SNPs). In the individuals from an Egyptian strain and in individuals caught in the wild in the basin of the river Volta, 197 and 163 SNPs were polymorphic, respectively. A pairwise calculation of Nei's genetic distance allowed the discrimination of the individual strains and species based on the genotypes determined with the SNP set. We expect that this set will be widely applicable for use in tilapia aquaculture, e.g. for pedigree reconstruction. In addition, this set is currently used for assaying the genetic diversity of native Nile tilapia in areas where tilapia is, or will be, introduced in aquaculture projects. This allows the tracing of escapees from aquaculture and the monitoring of effects of introgression and hybridization.

Genetic correlation between heart ratio and body weight as a function of ascites frequency in broilers split up into sex and health status.

Ascites or pulmonary hypertension syndrome is a metabolic disorder in broilers. Male broilers have a higher BW and are therefore expected to be more prone to developing ascites than females. As genetic parameters might be affected by the ascites incidence, genetic parameters might differ between male and female broilers. The aims of this study were to estimate the heritability for the ratio of right ventricular weight to total ventricular weight (RATIO) and BW in male and female broilers, the genetic correlation between RATIO and BW separately for male and female broilers, and the genetic correlations between BW for ascitic and nonascitic broilers. Data were available from 7,856 broilers (3,819 males and 4,037 females). The broilers in the experiment were kept under a cold temperature regimen and increased CO(2) levels. In this study, we showed that the incidence of ascites is higher in male than in female broilers. Heritability estimates for BW at 7 wk of age were higher for male (0.22) than for female (0.17) broilers, and for RATIO heritability, estimates were higher for female (0.44) than for male (0.32) broilers. The genetic correlations between RATIO and BW measured at different ages changed from slightly positive at 2 wk of age to moderately negative at 7 wk. The change in genetic correlation was more extreme for male (from 0.01 to -0.62) than for female (from 0.13 to -0.24) broilers. The difference in ascites incidence between male and female broilers is the most likely reason for the difference in genetic correlations. The genetic correlation between BW traits measured in broilers with fluid in the abdomen and without fluid in the abdomen decreased from 0.91 at 2 wk to 0.69 at 7 wk. We conclude that under circumstances with ascites, data from male and female broilers should be analyzed separately.