Spatial genetic structure of European wild boar, with inferences on late-Pleistocene and Holocene demographic history.

European wildlife has been subjected to intensifying levels of anthropogenic impact throughout the Holocene, yet the main genetic partitioning of many species is thought to still reflect the late-Pleistocene glacial refugia. We analyzed 26,342 nuclear SNPs of 464 wild boar (Sus scrofa) across the European continent to infer demographic history and reassess the genetic consequences of natural and anthropogenic forces. We found that population fragmentation, inbreeding and recent hybridization with domestic pigs have caused the spatial genetic structure to be heterogeneous at the local scale. Underlying local anthropogenic signatures, we found a deep genetic structure in the form of an arch-shaped cline extending from the Dinaric Alps, via Southeastern Europe and the Baltic states, to Western Europe and, finally, to the genetically diverged Iberian peninsula. These findings indicate that, despite considerable anthropogenic influence, the deeper, natural continental structure is still intact. Regarding the glacial refugia, our findings show a weaker signal than generally assumed, but are nevertheless suggestive of two main recolonization routes, with important roles for Southern France and the Balkans. Our results highlight the importance of applying genomic resources and framing genetic results within a species' demographic history and geographic distribution for a better understanding of the complex mixture of underlying processes.

Corrigendum: A novel loss-of-function variant in transmembrane protein 263 (TMEM263) of autosomal dwarfism in chicken.

[This corrects the article DOI: 10.3389/fgene.2018.00193.].

Genomic consequences of a century of inbreeding and isolation in the Danish wild boar population.

Demographic events such as series of bottlenecks impact the genetic variation and adaptive potential of populations. European megafauna, such as wild boars (Sus scrofa), have experienced severe climatic and size fluctuations that have shaped their genetic variation. Habitat fragmentation and human-mediated translocations have further contributed to the complex demographic history of European wild boar. Danish wild boars represent an extreme case of a small and isolated population founded by four wild boars from Germany. Here, we explore the genetic composition of the Danish wild boar population in Klelund. We genotyped all 21 Danish wild boars that were recently transferred from the source population in Lille Vildmose into the Klelund Plantation to establish a novel wild boar population. We compared the Danish wild boars with high-density single-nucleotide polymorphism genotypes from a comprehensive reference set of 1263 wild and domesticated pigs, including 11 individuals from Ulm, one of two presumed founder locations in Germany. Our findings support the European wild background of the Danish population, and no traces of gene flow with wild or domesticated pigs were found. The narrow genetic origin of the Danish wild boars is illustrated by extremely long and frequent runs of homozygous stretches in their genomes, indicative of recent inbreeding. This study provides the first insights into one of the most inbred wild boar populations globally established a century ago from a narrow base of only four founders.

The Visayan Warty Pig (Sus cebifrons) Genome Provides Insight Into Chromosome Evolution and Sensory Adaptation in Pigs.

It is largely unknown how mammalian genomes evolve under rapid speciation and environmental adaptation. An excellent model for understanding fast evolution is provided by the genus Sus, which diverged relatively recently and lacks postzygotic isolation. Here, we present a high-quality reference genome of the Visayan warty pig, which is specialized to a tropical island environment. Comparing the genome sequences and chromatin contact maps of the Visayan warty pig (Sus cebifrons) and domestic pig (Sus scrofa), we characterized the dynamics of chromosomal structure evolution during Sus speciation, revealing the similar chromosome conformation as the potential biological mechanism of frequent postdivergence hybridization among Suidae. We further investigated the different signatures of adaptive selection and domestication in Visayan warty pig and domestic pig with specific emphasize on the evolution of olfactory and gustatory genes, elucidating higher olfactory diversity in Visayan warty pig and positive and relaxed evolution of bitter and fat taste receptors, respectively, in domestic pig. Our comprehensive evolutionary and comparative genome analyses provide insight into the dynamics of genomes and how these change over relative short evolutionary times, as well as how these genomic differences encode for differences in the phenotypes.

A new haplotype-resolved turkey genome to enable turkey genetics and genomics research.

BACKGROUND: The domesticated turkey (Meleagris gallopavo) is a species of significant agricultural importance and is the second largest contributor, behind broiler chickens, to world poultry meat production. The previous genome is of draft quality and partly based on the chicken (Gallus gallus) genome. A high-quality reference genome of M. gallopavo is essential for turkey genomics and genetics research and the breeding industry. RESULTS: By adopting the trio-binning approach, we were able to assemble a high-quality chromosome-level F1 assembly and 2 parental haplotype assemblies, leveraging long-read technologies and genome-wide chromatin interaction data (Hi-C). From a total of 40 chromosomes (2n = 80), we captured 35 chromosomes in a single scaffold, showing much improved genome completeness and continuity compared to the old assembly build. The 3 assemblies are of higher quality than the previous draft quality assembly and comparable to the chicken assemblies (GRCg7) shown by the largest contig N50 (26.6 Mb) and comparable BUSCO gene set completeness scores (96-97%). Comparative analyses confirm a previously identified large inversion of around 19 Mbp on the Z chromosome not found in other Galliformes. Structural variation between the parent haplotypes was identified, which poses potential new target genes for breeding. CONCLUSIONS: We contribute a new high-quality turkey genome at the chromosome level, benefiting turkey genetics and other avian genomics research as well as the turkey breeding industry.

Genome-wide association analysis of adaptation to oxygen stress in Nile tilapia (Oreochromis niloticus).

BACKGROUND: Tilapia is one of the most abundant species in aquaculture. Hypoxia is known to depress growth rate, but the genetic mechanism by which this occurs is unknown. In this study, two groups consisting of 3140 fish that were raised in either aerated (normoxia) or non-aerated pond (nocturnal hypoxia). During grow out, fish were sampled five times to determine individual body weight (BW) gains. We applied a genome-wide association study to identify SNPs and genes associated with the hypoxic and normoxic environments in the 16th generation of a Genetically Improved Farmed Tilapia population. RESULTS: In the hypoxic environment, 36 SNPs associated with at least one of the five body weight measurements (BW1 till BW5), of which six, located between 19.48 Mb and 21.04 Mb on Linkage group (LG) 8, were significant for body weight in the early growth stage (BW1 to BW2). Further significant associations were found for BW in the later growth stage (BW3 to BW5), located on LG1 and LG8. Analysis of genes within the candidate genomic region suggested that MAPK and VEGF signalling were significantly involved in the later growth stage under the hypoxic environment. Well-known hypoxia-regulated genes such as igf1rb, rora, efna3 and aurk were also associated with growth in the later stage in the hypoxic environment. Conversely, 13 linkage groups containing 29 unique significant and suggestive SNPs were found across the whole growth period under the normoxic environment. A meta-analysis showed that 33 SNPs were significantly associated with BW across the two environments, indicating a shared effect independent of hypoxic or normoxic environment. Functional pathways were involved in nervous system development and organ growth in the early stage, and oocyte maturation in the later stage. CONCLUSIONS: There are clear genotype-growth associations in both normoxic and hypoxic environments, although genome architecture involved changed over the growing period, indicating a transition in metabolism along the way. The involvement of pathways important in hypoxia especially at the later growth stage indicates a genotype-by-environment interaction, in which MAPK and VEGF signalling are important components.

Accelerated discovery of functional genomic variation in pigs.

The genotype-phenotype link is a major research topic in the life sciences but remains highly complex to disentangle. Part of the complexity arises from the number of genes contributing to the observed phenotype. Despite the vast increase of molecular data, pinpointing the causal variant underlying a phenotype of interest is still challenging. In this study, we present an approach to map causal variation and molecular pathways underlying important phenotypes in pigs. We prioritize variation by utilizing and integrating predicted variant impact scores (pCADD), functional genomic information, and associated phenotypes in other mammalian species. We demonstrate the efficacy of our approach by reporting known and novel causal variants, of which many affect non-coding sequences. Our approach allows the disentangling of the biology behind important phenotypes by accelerating the discovery of novel causal variants and molecular mechanisms affecting important phenotypes in pigs. This information on molecular mechanisms could be applicable in other mammalian species, including humans.

Genetic consequences of long-term small effective population size in the critically endangered pygmy hog.

Increasing human disturbance and climate change have a major impact on habitat integrity and size, with far-reaching consequences for wild fauna and flora. Specifically, population decline and habitat fragmentation result in small, isolated populations. To what extend different endangered species can cope with small population size is still largely unknown. Studies on the genomic landscape of these species can shed light on past demographic dynamics and current genetic load, thereby also providing guidance for conservation programs. The pygmy hog (Porcula salvania) is the smallest and rarest wild pig in the world, with current estimation of only a few hundred living in the wild. Here, we analyzed whole-genome sequencing data of six pygmy hogs, three from the wild and three from a captive population, along with 30 pigs representing six other Suidae. First, we show that the pygmy hog had a very small population size with low genetic diversity over the course of the past ~1 million years. One indication of historical small effective population size is the absence of mitochondrial variation in the six sequenced individuals. Second, we evaluated the impact of historical demography. Runs of homozygosity (ROH) analysis suggests that the pygmy hog population has gone through past but not recent inbreeding. Also, the long-term, extremely small population size may have led to the accumulation of harmful mutations suggesting that the accumulation of deleterious mutations is exceeding purifying selection in this species. Thus, care has to be taken in the conservation program to avoid or minimize the potential for further inbreeding depression, and guard against environmental changes in the future.

Parallel Genomic Changes Drive Repeated Evolution of Placentas in Live-Bearing Fish.

The evolutionary origin of complex organs challenges empirical study because most organs evolved hundreds of millions of years ago. The placenta of live-bearing fish in the family Poeciliidae represents a unique opportunity to study the evolutionary origin of complex organs, because in this family a placenta evolved at least nine times independently. It is currently unknown whether this repeated evolution is accompanied by similar, repeated, genomic changes in placental species. Here, we compare whole genomes of 26 poeciliid species representing six out of nine independent origins of placentation. Evolutionary rate analysis revealed that the evolution of the placenta coincides with convergent shifts in the evolutionary rate of 78 protein-coding genes, mainly observed in transporter- and vesicle-located genes. Furthermore, differences in sequence conservation showed that placental evolution coincided with similar changes in 76 noncoding regulatory elements, occurring primarily around genes that regulate development. The unexpected high occurrence of GATA simple repeats in the regulatory elements suggests an important function for GATA repeats in developmental gene regulation. The distinction in molecular evolution observed, with protein-coding parallel changes more often found in metabolic and structural pathways, compared with regulatory change more frequently found in developmental pathways, offers a compelling model for complex trait evolution in general: changing the regulation of otherwise highly conserved developmental genes may allow for the evolution of complex traits.

Prioritizing sequence variants in conserved non-coding elements in the chicken genome using chCADD.

The availability of genomes for many species has advanced our understanding of the non-protein-coding fraction of the genome. Comparative genomics has proven itself to be an invaluable approach for the systematic, genome-wide identification of conserved non-protein-coding elements (CNEs). However, for many non-mammalian model species, including chicken, our capability to interpret the functional importance of variants overlapping CNEs has been limited by current genomic annotations, which rely on a single information type (e.g. conservation). We here studied CNEs in chicken using a combination of population genomics and comparative genomics. To investigate the functional importance of variants found in CNEs we develop a ch(icken) Combined Annotation-Dependent Depletion (chCADD) model, a variant effect prediction tool first introduced for humans and later on for mouse and pig. We show that 73 Mb of the chicken genome has been conserved across more than 280 million years of vertebrate evolution. The vast majority of the conserved elements are in non-protein-coding regions, which display SNP densities and allele frequency distributions characteristic of genomic regions constrained by purifying selection. By annotating SNPs with the chCADD score we are able to pinpoint specific subregions of the CNEs to be of higher functional importance, as supported by SNPs found in these subregions are associated with known disease genes in humans, mice, and rats. Taken together, our findings indicate that CNEs harbor variants of functional significance that should be object of further investigation along with protein-coding mutations. We therefore anticipate chCADD to be of great use to the scientific community and breeding companies in future functional studies in chicken.

Correction to: Genome-wide SNP data unveils the globalization of domesticated pigs.

An amendment to this paper has been published and can be accessed via the original article.

Parallel Genetic Origin of Foot Feathering in Birds.

Understanding the genetic basis of similar phenotypes shared between lineages is a long-lasting research interest. Even though animal evolution offers many examples of parallelism, for many phenotypes little is known about the underlying genes and mutations. We here use a combination of whole-genome sequencing, expression analyses, and comparative genomics to study the parallel genetic origin of ptilopody (Pti) in chicken. Ptilopody (or foot feathering) is a polygenic trait that can be observed in domesticated and wild avian species and is characterized by the partial or complete development of feathers on the ankle and feet. In domesticated birds, ptilopody is easily selected to fixation, though extensive variation in the type and level of feather development is often observed. By means of a genome-wide association analysis, we identified two genomic regions associated with ptilopody. At one of the loci, we identified a 17-kb deletion affecting PITX1 expression, a gene known to encode a transcription regulator of hindlimb identity and development. Similarly to pigeon, at the second loci, we observed ectopic expression of TBX5, a gene involved in forelimb identity and a key determinant of foot feather development. We also observed that the trait evolved only once as foot-feathered birds share the same haplotype upstream TBX5. Our findings indicate that in chicken and pigeon ptilopody is determined by the same set of genes that affect similar molecular pathways. Our study confirms that ptilopody has evolved through parallel evolution in chicken and pigeon.

Mixed ancestry from wild and domestic lineages contributes to the rapid expansion of invasive feral swine.

Invasive alien species are a significant threat to both economic and ecological systems. Identifying the processes that give rise to invasive populations is essential for implementing effective control strategies. We conducted an ancestry analysis of invasive feral swine (Sus scrofa, Linnaeus, 1758), a highly destructive ungulate that is widely distributed throughout the contiguous United States, to describe introduction pathways, sources of newly emergent populations and processes contributing to an ongoing invasion. Comparisons of high-density single nucleotide polymorphism genotypes for 6,566 invasive feral swine to a comprehensive reference set of S. scrofa revealed that the vast majority of feral swine were of mixed ancestry, with dominant genetic associations to Western heritage breeds of domestic pig and European populations of wild boar. Further, the rapid expansion of invasive feral swine over the past 30 years was attributable to secondary introductions from established populations of admixed ancestry as opposed to direct introductions of domestic breeds or wild boar. Spatially widespread genetic associations of invasive feral swine to European wild boar deviated strongly from historical S. scrofa introduction pressure, which was largely restricted to domestic pigs with infrequent, localized wild boar releases. The deviation between historical introduction pressure and contemporary genetic ancestry suggests wild boar-hybridization may contribute to differential fitness in the environment and heightened invasive potential for individuals of admixed domestic pig-wild boar ancestry.

pCADD: SNV prioritisation in Sus scrofa.

BACKGROUND: In animal breeding, identification of causative genetic variants is of major importance and high economical value. Usually, the number of candidate variants exceeds the number of variants that can be validated. One way of prioritizing probable candidates is by evaluating their potential to have a deleterious effect, e.g. by predicting their consequence. Due to experimental difficulties to evaluate variants that do not cause an amino-acid substitution, other prioritization methods are needed. For human genomes, the prediction of deleterious genomic variants has taken a step forward with the introduction of the combined annotation dependent depletion (CADD) method. In theory, this approach can be applied to any species. Here, we present pCADD (p for pig), a model to score single nucleotide variants (SNVs) in pig genomes. RESULTS: To evaluate whether pCADD captures sites with biological meaning, we used transcripts from miRNAs and introns, sequences from genes that are specific for a particular tissue, and the different sites of codons, to test how well pCADD scores differentiate between functional and non-functional elements. Furthermore, we conducted an assessment of examples of non-coding and coding SNVs, which are causal for changes in phenotypes. Our results show that pCADD scores discriminate between functional and non-functional sequences and prioritize functional SNVs, and that pCADD is able to score the different positions in a codon relative to their redundancy. Taken together, these results indicate that based on pCADD scores, regions with biological relevance can be identified and distinguished according to their rate of adaptation. CONCLUSIONS: We present the ability of pCADD to prioritize SNVs in the pig genome with respect to their putative deleteriousness, in accordance to the biological significance of the region in which they are located. We created scores for all possible SNVs, coding and non-coding, for all autosomes and the X chromosome of the pig reference sequence Sscrofa11.1, proposing a toolbox to prioritize variants and evaluate sequences to highlight new sites of interest to explain biological functions that are relevant to animal breeding.

The Genomes of the Livebearing Fish Species Poeciliopsis retropinna and Poeciliopsis turrubarensis Reflect Their Different Reproductive Strategies.

The evolution of a placenta is predicted to be accompanied by rapid evolution of genes involved in processes that regulate mother-offspring interactions during pregnancy, such as placenta formation, embryonic development, and nutrient transfer to offspring. However, these predictions have only been tested in mammalian species, where only a single instance of placenta evolution has occurred. In this light, the genus Poeciliopsis is a particularly interesting model for placenta evolution, because in this genus a placenta has evolved independently from the mammalian placenta. Here, we present and compare genome assemblies of two species of the livebearing fish genus Poeciliopsis (family Poeciliidae) that differ in their reproductive strategy: Poeciliopsis retropinna which has a well-developed complex placenta and P. turrubarensis which lacks a placenta. We applied different assembly strategies for each species: PacBio sequencing for P. retropinna (622-Mb assembly, scaffold N50 of 21.6 Mb) and 10× Genomics Chromium technology for P. turrubarensis (597-Mb assembly, scaffold N50 of 4.2 Mb). Using the high contiguity of these genome assemblies and near-completeness of gene annotations to our advantage, we searched for gene duplications and performed a genome-wide scan for genes evolving under positive selection. We find rapid evolution in major parts of several molecular pathways involved in parent-offspring interaction in P. retropinna, both in the form of gene duplications as well as positive selection. We conclude that the evolution of the placenta in the genus Poeciliopsis is accompanied by rapid evolution of genes involved in similar genomic pathways as found in mammals.

The type of bottleneck matters: Insights into the deleterious variation landscape of small managed populations.

Predictions about the consequences of a small population size on genetic and deleterious variation are fundamental to population genetics. As small populations are more affected by genetic drift, purifying selection acting against deleterious alleles is predicted to be less efficient, therefore increasing the risk of inbreeding depression. However, the extent to which small populations are subjected to genetic drift depends on the nature and time frame in which the bottleneck occurs. Domesticated species are an excellent model to investigate the consequences of population bottlenecks on genetic and deleterious variation in small populations. This is because their history is dominated by known bottlenecks associated with domestication, breed formation and intense selective breeding. Here, we use whole-genome sequencing data from 97 chickens representing 39 traditional fancy breeds to directly examine the consequences of two types of bottlenecks for deleterious variation: the severe domestication bottleneck and the recent population decline accompanying breed formation. We find that recently bottlenecked populations have a higher proportion of deleterious variants relative to populations that have been kept at small population sizes since domestication. We also observe that long tracts of homozygous genotypes (runs of homozygosity) are proportionally more enriched in deleterious variants than the rest of the genome. This enrichment is particularly evident in recently bottlenecked populations, suggesting that homozygosity of these variants is likely to occur due to genetic drift and recent inbreeding. Our results indicate that the timing and nature of population bottlenecks can substantially shape the deleterious variation landscape in small populations.

Detection of a Frameshift Deletion in the SPTBN4 Gene Leads to Prevention of Severe Myopathy and Postnatal Mortality in Pigs.

Piglet mortality is a complex phenotype that depends on the environment, selection on piglet health, but also on the interaction between the piglet and sow. However, also monogenic recessive defects contribute to piglet mortality. Selective breeding has decreased overall piglet mortality by improving both mothering abilities and piglet viability. However, variants underlying recessive monogenic defects are usually not well captured within the breeding values, potentially drifting to higher frequency as a result of intense selection or genetic drift. This study describes the identification by whole-genome sequencing of a recessive 16-bp deletion in the SPTBN4 gene causing postnatal mortality in a pig breeding line. The deletion induces a frameshift and a premature stop codon, producing an impaired and truncated spectrin beta non-erythrocytic 4 protein (SPTBN4). Applying medium density single nucleotide polymorphism (SNP) data available for all breeding animals, a pregnant carrier sow sired by a carrier boar was identified. Of the resulting piglets, two confirmed homozygous piglets suffered from severe myopathy, hind-limb paralysis, and tremors. Histopathological examination showed dispersed degeneration and decrease of cross-striations in the dorsal and hind-limb muscle fibers of the affected piglets. Hence, the affected piglets are unable to walk or drink, usually resulting in death within a few hours after birth. This study demonstrates how growing genomic resources in pig breeding can be applied to identify rare syndromes in breeding populations, that are usually poorly documented and often are not even known to have a genetic basis. The study allows to prevent carrier-by-carrier matings, thereby gradually decreasing the frequency of the detrimental allele and avoiding the birth of affected piglets, improving animal welfare. Finally, these "natural knockouts" increase our understanding of gene function within the mammalian clade, and provide a potential model for human disease.

Author Correction: Hotspots of recent hybridization between pigs and wild boars in Europe.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Impact of merging commercial breeding lines on the genetic diversity of Landrace pigs.

BACKGROUND: The pig breeding industry has undergone a large number of mergers in the past decades. Various commercial lines were merged or discontinued, which is expected to reduce the genetic diversity of the pig species. The objective of the current study was to investigate the genetic diversity of different former Dutch Landrace breeding lines and quantify their relationship with the current Dutch Landrace breed that originated from these lines. RESULTS: Principal component analysis clearly divided the former Landrace lines into two main clusters, which are represented by Norwegian/Finnish Landrace lines and Dutch Landrace lines. Structure analysis revealed that each of the lines that are present in the Dutch Gene bank has a unique genetic identity. The current Dutch Landrace breed shows a high level of admixture and is closely related to the six former lines. The Dumeco N-line, which is conserved in the Dutch Gene bank, is poorly represented in the current Dutch Landrace. All seven lines (the six former and the current line) contribute almost equally to the genetic diversity of the Dutch Landrace breed. As expected, the current Dutch Landrace breed comprises only a small proportion of unique genetic diversity that was not present in the other lines. The genetic diversity level, as measured by Eding's core set method, was equal to 0.89 for the current Dutch Landrace breed, whereas total genetic diversity across the seven lines, measured by the same method, was equal to 0.99. CONCLUSIONS: The current Dutch Landrace breed shows a high level of admixture and is closely related to the six former Dutch Landrace lines. Merging of commercial Landrace lines has reduced the genetic diversity of the Landrace population in the Netherlands, although a large proportion of the original variation is maintained. Thus, our recommendation is to conserve breeding lines in a gene bank before they are merged.