Search

The German Longhaired Pointer (GLP) breed is a versatile pointer dog breed. In the current study, we investigated the genetic diversity of these dogs based on SNP array data and compared it to 11 other pointer setter breeds. The results show that GLPs have a relatively low level of inbreeding among these pointer breeds. In addition, with the availability of pedigree information of the GLPs, we demonstrate that the correlation between pedigree-based inbreeding and genotype-based inbreeding coefficients was high (R = 0.89 and 0.85). By investigating population structure between these 12 pointer setter breeds we showed that GLP is a breed distinct from other pointers and shares common ancestry with a few other pointing breeds. Finally, we identified selection signatures in GLPs using the runs of homozygosity islands method. The most significant runs of homozygosity island was detected on chromosome 30 harboring the genes RYR3, FMN1, and GREM1. The RYR3 gene plays a role in skeletal muscle contraction while the FMN1 and GREM1 genes are involved in limb development. The selection on these three genes could have contributed to the excellent athletic performance of GLPs, which is an extremely important characteristic for this hunting dog.

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. Nasl. z nasl. zaslona. Opis vira z dne 3. 2. 02023. Skupno št. avtorjev: 21. Bibliografija: str. 142-143.

Thyroid carcinomas originating from follicular cells of the thyroid gland occur in both humans and dogs and they have highly similar histomorphologic patterns. In dogs, thyroid carcinomas have not been extensively investigated, especially concerning the familial origin of thyroid carcinomas. Here we report familial thyroid follicular cell carcinomas confirmed by histology in 54 Dutch origin German longhaired pointers. From the pedigree, 45 of 54 histopathologically confirmed cases are closely related to a pair of first-half cousins in the past, indicating a familial disease. In addition, genetics contributed more to the thyroid follicular cell carcinoma than other factors by an estimated heritability of 0.62 based on pedigree. The age of diagnosis ranged between 4.5 and 13.5 years, and 76% of cases were diagnosed before 10 years of age, implying an early onset of disease. We observed a significant higher pedigree-based inbreeding coefficient in the affected dogs (mean F 0.23) compared to unaffected dogs (mean F 0.14), suggesting the contribution of inbreeding to tumour development. The unique occurrence of familial thyroid follicular cell carcinoma in this dog population and the large number of affected dogs make this population an important model to identify the genetic basis of familial thyroid follicular cell carcinoma in this breed and may contribute to the research into pathogenesis, prevention and treatment in humans. This article is protected by copyright. All rights reserved.

Recent papers postulate that epifaunal organisms use artificial structures as stepping‐stones to spread to areas that are too distant to reach in a single generation. With thousands of artificial structures present in the North Sea, we test the hypothesis that these structures are connected by water currents and act as an interconnected reef. Population genetic structure of the Blue mussel, Mytilus edulis was expected to follow a pattern predicted by particle tracking models (PTM). Correlation between population genetic differentiation, based on microsatellite markers, and particle exchange was tested. Specimens of M. edulis were found at each location, although the PTM indicated that locations >85 km offshore were isolated from coastal sub‐populations. Fixation coefficient FST correlated with the number of arrivals in the PTM. However, the number of effective migrants per generation as inferred from coalescent simulations did not show a strong correlation with the arriving particles. Isolation by distance analysis showed no increase in isolation with increasing distance and we did not find clear structure among the populations. The marine stepping‐stone effect is obviously important for the distribution of M. edulis in the North Sea and it may influence ecologically comparable species in a similar way. In the absence of artificial shallow hard substrates, M. edulis would be unlikely to survive in offshore North Sea waters. Although we found an indication that FST was lower between connected locations, isolation by distance analysis showed no increase in isolation with increasing distance. Finally, we did not find clear structure among the populations.

Here, we aimed to identify and characterize genomic regions that differ between Groningen White Headed (GWH) breed and other cattle, and in particular to identify candidate genes associated with coat color and/or eye-protective phenotypes. Firstly, whole genome sequences of 170 animals from eight breeds were used to evaluate the genetic structure of the GWH in relation to other cattle breeds by carrying out principal components and model-based clustering analyses. Secondly, the candidate genomic regions were identified by integrating the findings from: a) a genome-wide association study using GWH, other white headed breeds (Hereford and Simmental), and breeds with a non-white headed phenotype (Dutch Friesian, Deep Red, Meuse-Rhine-Yssel, Dutch Belted, and Holstein Friesian); b) scans for specific signatures of selection in GWH cattle by comparison with four other Dutch traditional breeds (Dutch Friesian, Deep Red, Meuse-Rhine-Yssel and Dutch Belted) and the commercial Holstein Friesian; and c) detection of candidate genes identified via these approaches. The alignment of the filtered reads to the reference genome (ARS-UCD1.2) resulted in a mean depth of coverage of 8.7X. After variant calling, the lowest number of breed-specific variants was detected in Holstein Friesian (148,213), and the largest in Deep Red (558,909). By integrating the results, we identified five genomic regions under selection on BTA4 (70.2–71.3 Mb), BTA5 (10.0–19.7 Mb), BTA20 (10.0–19.9 and 20.0–22.7 Mb), and BTA25 (0.5–9.2 Mb). These regions contain positional and functional candidate genes associated with retinal degeneration (e.g., CWC27 and CLUAP1), ultraviolet protection (e.g., ERCC8), and pigmentation (e.g. PDE4D) which are probably associated with the GWH specific pigmentation and/or eye-protective phenotypes, e.g. Ambilateral Circumocular Pigmentation (ACOP). Our results will assist in characterizing the molecular basis of GWH phenotypes and the biological implications of its adaptation.

Familial thyroid cancer originating from follicular cells accounts for 5–15% of all the thyroid carcinoma cases in humans. Previously, we described thyroid follicular cell carcinomas in a large number of the Dutch German longhaired pointers (GLPs) with a likely autosomal recessive inheritance pattern. Here, we investigated the genetic causes of the disease using a combined approach of genome-wide association study and runs of homozygosity (ROH) analysis based on 170k SNP array genotype data and whole-genome sequences. A region 0–5 Mb on chromosome 17 was identified to be associated with the disease. Whole-genome sequencing revealed many mutations fitting the recessive inheritance pattern in this region including two deleterious mutations in the TPO gene, chr17:800788G&gt, A (686F&gt, V) and chr17:805276C&gt, T (845T&gt, M). These two SNP were subsequently genotyped in 186 GLPs (59 affected and 127 unaffected) and confirmed to be highly associated with the disease. The recessive genotypes had higher relative risks of 16.94 and 16.64 compared to homozygous genotypes for the reference alleles, respectively. This study provides novel insight into the genetic causes leading to the familial thyroid follicular cell carcinoma, and we were able to develop a genetic test to screen susceptible dogs.

Background China has the richest local chicken breeding resources in the world and is the world’s second largest producer of meat-type chickens. Development of a moderate-density SNP array for genetic analysis of chickens and breeding of meat-type chickens taking utility of those resources is urgently needed for conventional farms, breeding industry, and research areas. Results Eight representative local breeds or commercial broiler lines with 3 pools of 48 individuals within each breed/line were sequenced and supplied the major SNPs resource. There were 7.09 million - 9.41 million SNPs detected in each breed/line. After filtering using multiple criteria such as preferred incorporation of trait-related SNPs and uniformity of distribution across the genome, 52.18 K SNPs were selected in the final array. It consists of: (i) 19.22 K SNPs from the genomes of yellow-feathered, cyan-shank partridge and white-feathered chickens; (ii) 5.98 K SNPs related to economic traits from the Illumina 60 K SNP Bead Chip, which were found as significant associated SNPs with 15 traits in a Beijing-You crossed Cobb F2 resource population by genome-wide association study analysis; (iii) 7.63 K SNPs from 861 candidate genes of economic traits; (iv) the 0.94 K SNPs related to residual feed intake; and (v) 18.41 K from chicken SNPdb. The polymorphisms of 9 extra local breeds and 3 commercial lines were examined with this array, and 40 K - 47 K SNPs were polymorphic (with minor allele frequency > 0.05) in those breeds. The MDS result showed that those breeds can be clearly distinguished by this newly developed genotyping array. Conclusions We successfully developed a 55K genotyping array by using SNPs segregated from typical local breeds and commercial lines. Compared to the existing Affy 600 K and Illumina 60 K arrays, there were 21,41 K new SNPs included on our Affy 55K array. The results of the 55K genotyping data can therefore be imputed to high-density SNPs genotyping data. The array offers a wide range of potential applications such as genomic selection breeding, GWAS of interested traits, and investigation of diversity of different chicken breeds. Electronic supplementary material The online version of this article (10.1186/s12864-019-5736-8) contains supplementary material, which is available to authorized users.

Thyroid carcinomas originating from follicular cells of the thyroid gland occur in both humans and dogs and they have highly similar histomorphologic patterns. In dogs, thyroid carcinomas have not been extensively investigated, especially concerning the familial origin of thyroid carcinomas. Here we report familial thyroid follicular cell carcinomas confirmed by histology in 54 Dutch origin German longhaired pointers. From the pedigree, 45 of 54 histopathologically confirmed cases are closely related to a pair of first-half cousins in the past, indicating a familial disease. In addition, genetics contributed more to the thyroid follicular cell carcinoma than other factors by an estimated heritability of 0.62 based on pedigree. The age of diagnosis ranged between 4.5 and 13.5 years, and 76% of cases were diagnosed before 10 years of age, implying an early onset of disease. We observed a significant higher pedigree-based inbreeding coefficient in the affected dogs (mean F 0.23) compared to unaffected dogs (mean F 0.14), suggesting the contribution of inbreeding to tumour development. The unique occurrence of familial thyroid follicular cell carcinoma in this dog population and the large number of affected dogs make this population an important model to identify the genetic basis of familial thyroid follicular cell carcinoma in this breed and may contribute to the research into pathogenesis, prevention and treatment in humans.

Structural variations (SVs) are an important source of phenotypic diversity in cattle. Here, 72 whole genome sequences representing taurine and zebu cattle were used to identify SVs. Applying multiple approaches, 16,738 SVs were identified. A comparison against the Database of Genomic Variants archives revealed that 1575 SVs were novel in our data. A novel duplication covering the entire GALNT15 gene, was observed only in N'Dama. A duplication, which was previously reported only in zebu and associated with navel length, was also observed in N'Dama. Investigation of a novel deletion located upstream of CAST13 gene and identified only in Italian cattle and zebu, revealed its introgressed origin in the former. Overall, our data highlights how the SVs distribution in cattle is also shaped by forces such as demographical differences and gene flow. The cattle SVs of this study and its meta-data can be visualized on an interactive genome browser at https://tinyurl.com/svCowArs

The growth of animals is a complex trait, in chicken resulting in a diverse variety of forms, caused by a heterogeneous genetic basis. Bantam chicken, known as an exquisite form of dwarfism, has been used for crossbreeding to create corresponding dwarf counterparts for native fowls in the Dutch populations. Here, we demonstrate the heterogeneity of the bantam trait in Dutch chickens and reveal the underlying genetic causes, using whole‐genome sequence data from matching pairs of bantam and normal‐sized breeds. During the bantam‐oriented crossbreeding, various bantam origins were used to introduce the bantam phenotype, and three major bantam sources were identified and clustered. The genome‐wide association studies revealed multiple genetic variants and genes associated with bantam phenotype, including HMGA2 and PRDM16, genes involved in body growth and stature. The comparison of associated variants among studies illustrated differences related to divergent bantam origins, suggesting a clear heterogeneity among bantam breeds. We show that in neo‐bantam breeds, the bantam‐related regions underwent a strong haplotype introgression from the bantam source, outcompeting haplotypes from the normal‐sized counterpart. The bantam heterogeneity is further confirmed by the presence of multiple haplotypes comprising associated alleles, which suggests the selection of the bantam phenotype is likely subject to a convergent direction across populations. Our study demonstrates that the diverse history of human‐mediated crossbreeding has contributed to the complexity and heterogeneity of the bantam phenotype.

The geographic distributions of some coastal marine species have appeared as cosmopolitan ever since they were first scientifically documented. In particular, for many benthic species that are associated with anthropogenic substrata, there is much speculation as to whether or not their broad distributions can be explained by natural mechanisms of dispersal. Here, we focused on two congeneric coastal crustaceans with cosmopolitan distributions—the tube-dwelling amphipodsJassa marmorataandJassa slatteryi. Both species are common elements of marine biofouling on nearly all kinds of artificial hard substrata in temperate to warm seas. We hypothesized that the two species’ modern occurrences across the oceans are the result of human shipping activities that started centuries ago. Mitochondrial DNA sequences of the CO1 fragment of specimens from distinct marine regions around the world were analysed, evaluating genetic structure and migration models and making inferences on putative native ranges of the twoJassaspecies. Populations of both species exhibited considerable genetic diversity with differing levels of geographic structure. For both species, at least two dominant haplotypes were shared among several geographic populations. Rapid demographic expansion and high migration rates between geographically distant regions support a scenario of ongoing dispersal all over the world. Our findings indicate that the likely former native range ofJ. marmoratais the Northwest Atlantic, whereas the likely former native range ofJ. slatteryiis the Northern Pacific region. As corroborated by the genetic connectivity between populations, shipping still appears to be the more successful vector of the two species’ dispersal when compared to natural mechanisms. Historical invasion events that likely started centuries ago, along with current ongoing dispersal, confirm these species’ identities as true “neocosmopolitans”.

Significance Archaeological evidence indicates that domestic pigs arrived in Europe, alongside farmers from the Near East ∼8,500 y ago, yet mitochondrial genomes of modern European pigs are derived from European wild boars. To address this conundrum, we obtained mitochondrial and nuclear data from modern and ancient Near Eastern and European pigs. Our analyses indicate that, aside from a coat color gene, most Near Eastern ancestry in the genomes of European domestic pigs disappeared over 3,000 y as a result of interbreeding with local wild boars. This implies that pigs were not domesticated independently in Europe, yet the first 2,500 y of human-mediated selection applied by Near Eastern Neolithic farmers played little role in the development of modern European pigs., Archaeological evidence indicates that pig domestication had begun by ∼10,500 y before the present (BP) in the Near East, and mitochondrial DNA (mtDNA) suggests that pigs arrived in Europe alongside farmers ∼8,500 y BP. A few thousand years after the introduction of Near Eastern pigs into Europe, however, their characteristic mtDNA signature disappeared and was replaced by haplotypes associated with European wild boars. This turnover could be accounted for by substantial gene flow from local European wild boars, although it is also possible that European wild boars were domesticated independently without any genetic contribution from the Near East. To test these hypotheses, we obtained mtDNA sequences from 2,099 modern and ancient pig samples and 63 nuclear ancient genomes from Near Eastern and European pigs. Our analyses revealed that European domestic pigs dating from 7,100 to 6,000 y BP possessed both Near Eastern and European nuclear ancestry, while later pigs possessed no more than 4% Near Eastern ancestry, indicating that gene flow from European wild boars resulted in a near-complete disappearance of Near East ancestry. In addition, we demonstrate that a variant at a locus encoding black coat color likely originated in the Near East and persisted in European pigs. Altogether, our results indicate that while pigs were not independently domesticated in Europe, the vast majority of human-mediated selection over the past 5,000 y focused on the genomic fraction derived from the European wild boars, and not on the fraction that was selected by early Neolithic farmers over the first 2,500 y of the domestication process.

Many ungulate populations have a complex history of isolation and translocation. Consequently, ungulate populations may have experienced substantial reductions in the level of overall gene flow, yet simultaneously have augmented levels of long-distance gene flow. To investigate the effect of this dual anthropogenic effect on the genetic landscape of ungulates, we genotyped 35K SNPs in 47 red deer (Cervus elaphus) of The Netherlands, including putative autochthonous relic populations as well as allochthonous populations established in private estates and rewilding areas. We applied FST and ordination analyses to determine the meta-population genetic structure and thereby the occurrence of hybridization. At population level, we investigated levels of inbreeding through individual-based diversity measures, including Runs of Homozygosity. We documented that both spatial genetic structure and within-population genetic variation differed markedly from patterns assumed from present-day abundance and distribution. Notwithstanding the small spatial scale, red deer populations formed distinct genetic clusters, and some had higher genetic similarity to distant than to nearby populations. Moreover, the putative autochthonous relic deer populations had much reduced levels of polymorphism and multi-locus heterozygosity, despite relatively large current population sizes. Accordingly, genomes of these deer contained a high proportion of long (>5Mb) Runs of Homozygosity. Whereas the observed high levels of inbreeding warrant defragmentation measures, the presence of adjacent autochthonous and allochthonous genetic stocks imply that facilitation of gene flow would cause genetic homogenization. Such distortions of the genetic landscape of ungulates creates management dilemmas that cannot be properly anticipated without baseline genetic monitoring.

Cattle domestication occurred at least twice independently and gave rise to the modern taurine and indicine cattle breeds. European cattle diversity is generally dominated by taurine cattle, although elevated levels of indicine ancestry have been recorded in several breeds from southern Europe. Here we use genome-wide high-density SNP genotyping data to investigate the taurine and indicine ancestry in southern European cattle, based on a dataset comprising 508 individuals from 23 cattle breeds of taurine, indicine and mixed ancestry, including three breeds from Central Italy known to exhibit the highest levels of indicine introgression among southern European breeds. Based on local genomic ancestry analyses, we reconstruct taurine and indicine ancestry genome-wide and along chromosomes. We scrutinise local genomic introgression signals and identify genomic regions that have introgressed from indicine into taurine cattle under positive selection, harbouring genes with functions related to body size and feed efficiency. These findings suggest that indicine-derived traits helped enhance Central Italian cattle through adaptive introgression. The identified genes could provide genomic targets for selection for improved cattle performance. Our findings elucidate the key role of adaptive introgression in shaping the phenotypic features of modern cattle, aided by cultural and livestock exchange among historic human societies.

Temperature and CO2 concentration during incubation have profound effects on broiler chick development, and numerous studies have identified significant effects on hatch heart weight (HW) as a result of differences in these parameters. Early life environment has also been shown to affect broiler performance later in life; it has thus been suggested that epigenetic mechanisms may mediate long-term physiological changes induced by environmental stimuli. DNA methylation is an epigenetic modification that can confer heritable changes in gene expression. Using reduced-representation bisulfite sequencing (RRBS), we assessed DNA methylation patterns in cardiac tissue of 84 broiler hatchlings incubated at two egg shell temperatures (EST; 37.8°C and 38.9°C) and three CO2 concentrations (0.1%, 0.4%, and 0.8%) from day 8 of incubation onward. We assessed differential methylation between EST treatments and identified 2,175 differentially methylated (DM) CpGs (1,121 hypermethylated, 1,054 hypomethylated at 38.9° vs. 37.8°) in 269 gene promoters and 949 intragenic regions. DM genes (DMGs) were associated with heart developmental processes, including cardiomyocyte proliferation and differentiation. We identified enriched binding motifs among DM loci, including those for transcription factors associated with cell proliferation and heart development among hypomethylated CpGs that suggest increased binding ability at higher EST. We identified 9,823 DM CpGs between at least two CO2 treatments, with the greatest difference observed between 0.8 and 0.1% CO2 that disproportionately impacted genes involved in cardiac muscle development and response to low oxygen levels. Using HW measurements from the same chicks, we performed an epigenome-wide association study (EWAS) for HW, and identified 23 significantly associated CpGs, nine of which were also DM between ESTs. We found corresponding differences in transcript abundance between ESTs in three DMGs (ABLIM2, PITX2, and THRSP). Hypomethylation of an exonic CpG in PITX2 at 38.9°C was associated with increased expression, and suggests increased cell proliferation in broiler hatchlings incubated at higher temperatures. Overall, these results identified numerous epigenetic associations between chick incubation factors and heart development that may manifest in long-term differences in animal performance.

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.

Nile tilapia (Oreochromis niloticus) is a globally significant aquaculture species rapidly gaining status as a farmed commodity. In West Africa, wild Nile tilapia genetic resources are abundant yet knowledge of fine-scale population structure and patterns of natural genetic variation are limited. Coinciding with this is a burgeoning growth in tilapia aquaculture in Ghana and other countries within the region underpinned by locally available genetic resources. Using 192 single nucleotide polymorphism (SNP) markers this study conducted a genetic survey of Nile tilapia throughout West Africa, sampling 23 wild populations across eight countries (Benin, Burkina Faso, Côte d’Ivoire, Ghana, Togo, Mali, Gambia and Senegal), representing the major catchments of the Volta, Niger, Senegal and Gambia River basins. A pattern of isolation-by-distance and significant spatial genetic structure was identified throughout West Africa (Global FST = 0.144), which largely corresponds to major river basins and, to a lesser extent, sub-basins. Two populations from the Gambia River (Kudang and Walekounda), one from the western Niger River (Lake Sélingué) and one from the upper Red Volta River (Kongoussi) showed markedly lower levels of diversity and high genetic differentiation compared to all other populations, suggesting genetically isolated populations occurring across the region. Genetic structure within the Volta Basin did not always follow the pattern expected for sub-river basins. This study identifies clear genetic structuring and differentiation amongst West African Nile tilapia populations, which concur with broad patterns found in previous studies. In addition, we provide new evidence for fine-scale genetic structuring within the Volta Basin and previously unidentified genetic differences of populations in Gambia. The 192 SNP marker suite used in this study is a useful tool for differentiating tilapia populations and we recommend incorporating this marker suite into future population screening of O. niloticus. Our results form the basis of a solid platform for future research on wild tilapia genetic resources in West Africa, and the identification of potentially valuable germplasm for use in ongoing breeding programs for aquaculture.

Summary Copy number variation (CNV) is a major type of structural genomic variation that is increasingly studied across different species for association with diseases and production traits. Established protocols for experimental detection and computational inference of CNVs from SNP array and next-generation sequencing data are available. We present the CNVRanger R/Bioconductor package which implements a comprehensive toolbox for structured downstream analysis of CNVs. This includes functionality for summarizing individual CNV calls across a population, assessing overlap with functional genomic regions, and genome-wide association analysis with gene expression and quantitative phenotypes. Availability and implementation http://bioconductor.org/packages/CNVRanger.

Offshore environments are increasingly invaded by man-made structures that form hard-substrate habitats for many marine species. Examples include oil and gas platforms, wind turbines and shipwrecks. One of the hypothesised effects is an increased genetic connectivity among natural populations due to new populations growing on man-made structures that may act as stepping stones. However, few data are available on genetic connectivity among organisms inhabiting artificial offshore structures. Here, we present a study on the common fouling amphipod Jassa herdmani from offshore structures in the southern North Sea. Partial mitochondrial DNA sequences (cytochrome-c-oxidase 1, N = 514) were obtained from artificial structures at 17 locations in the southern North Sea, including 13 shipwrecks, 2 wind turbines and 2 platforms. Samples from these locations were significantly differentiated, meaning that strong population structure exists for this species in the area. Levels of intraspecific variation were consistent with stable population sizes. No evidence was found for isolation by distance. Using coalescent simulations, the oldest population subdivision events were estimated to date back to the time the study area was flooded following the Last Glacial Maximum. We therefore tentatively conclude that J. herdmani may have colonised man-made structures from previously existing populations on the sea floor, and that the increase in offshore installations has not led to an overall increase in genetic connectivity for this species.

Chromosome inversions have clear effects on genome evolution and have been associated with speciation, adaptation, and the evolution of the sex chromosomes. In birds, these inversions may play an important role in hybridization of species and disassortative mating. We identified a large (≈64 Mb) inversion polymorphism in the great tit (Parus major) that encompasses almost 1,000 genes and more than 90% of Chromosome 1A. The inversion occurs at a low frequency in a set of over 2,300 genotyped great tits in the Netherlands with only 5% of the birds being heterozygous for the inversion. In an additional analysis of 29 resequenced birds from across Europe, we found two heterozygotes. The likely inversion breakpoints show considerable genomic complexity, including multiple copy number variable segments. We identified different haplotypes for the inversion, which differ in the degree of recombination in the center of the chromosome. Overall, this remarkable genetic variant is widespread among distinct great tit populations and future studies of the inversion haplotype, including how it affects the fitness of carriers, may help to understand the mechanisms that maintain it.

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.

International audience; Archaeological evidence indicates that pig domestication had begun by ∼10,500 y before the present (BP) in the Near East, and mitochondrial DNA (mtDNA) suggests that pigs arrived in Europe alongside farmers ∼8,500 y BP. A few thousand years after the introduction of Near Eastern pigs into Europe, however, their characteristic mtDNA signature disappeared and was replaced by haplotypes associated with European wild boars. This turnover could be accounted for by substantial gene flow from local Euro-pean wild boars, although it is also possible that European wild boars were domesticated independently without any genetic contribution from the Near East. To test these hypotheses, we obtained mtDNA sequences from 2,099 modern and ancient pig samples and 63 nuclear ancient genomes from Near Eastern and European pigs. Our analyses revealed that European domestic pigs dating from 7,100 to 6,000 y BP possessed both Near Eastern and European nuclear ancestry, while later pigs possessed no more than 4% Near Eastern ancestry, indicating that gene flow from European wild boars resulted in a near-complete disappearance of Near East ancestry. In addition, we demonstrate that a variant at a locus encoding black coat color likely originated in the Near East and persisted in European pigs. Altogether, our results indicate that while pigs were not independently domesticated in Europe, the vast majority of human-mediated selection over the past 5,000 y focused on the genomic fraction derived from the European wild boars, and not on the fraction that was selected by early Neolithic farmers over the first 2,500 y of the domestication process. domestication | evolution | gene flow | Neolithic

The divergence between indicine cattle (Bos indicus) and taurine cattle (Bos taurus) is estimated to have occurred approximately 250,000 years ago, but a small number of European cattle breeds still display shared ancestry with indicine cattle. Additionally, following the divergence of African and European taurine, the gene flow between African taurine and southern European cattle has also been proposed. However, the extent to which non‐European cattle ancestry is diffused across southern European cattle has not been investigated thoroughly. Also, in recent times, many local breeds have suffered severe reductions in effective population size. Therefore, in the present study, we investigated the pattern of genetic diversity in various European cattle based on single nucleotide polymorphisms (SNP) identified from whole‐genome sequencing data. Additionally, we also employed unlinked and phased SNP‐based approaches on high‐density SNP array data to characterize non‐European cattle ancestry in several southern European cattle breeds. Using heterozygosity‐based parameters, we concluded that, on average, nucleotide diversity is greater in southern European cattle than western European (British and commercial) cattle. However, an abundance of long runs of homozygosity (ROH) and the pattern of Linkage disequilibrium decay suggested recent bottlenecks in Maltese and Romagnola. High nucleotide diversity outside ROH indicated a highly diverse founder population for southern European and African taurine. We also show that Iberian cattle display shared ancestry with African cattle. Furthermore, we show that Podolica is an ancient cross‐bred between Indicine zebu and European taurine. Additionally, we also inferred similar ancestry profile of non‐European cattle ancestry in different Balkan and Italian cattle breeds which might be an indication of the common origin of indicine ancestry in these breeds. Finally, we discuss several plausible demographic scenarios which might account for the presence of non‐European cattle ancestry in these cattle breeds.

Background Sex-linked slow (SF) and fast (FF) feathering rates at hatch have been widely used in poultry breeding for autosexing at hatch. In chicken, the sex-linked K (SF) and k+ (FF) alleles are responsible for the feathering rate phenotype. Allele K is dominant and a partial duplication of the prolactin receptor gene has been identified as the causal mutation. Interestingly, some domesticated turkey lines exhibit similar slow- and fast-feathering phenotypes, but the underlying genetic components and causal mutation have never been investigated. In this study, our aim was to investigate the molecular basis of feathering rate at hatch in domestic turkey. Results We performed a sequence-based case–control association study and detected a genomic region on chromosome Z, which is statistically associated with rate of feathering at hatch in turkey. We identified a 5-bp frameshift deletion in the prolactin receptor (PRLR) gene that is responsible for slow feathering at hatch. All female cases (SF turkeys) were hemizygous for this deletion, while 188 controls (FF turkeys) were hemizygous or homozygous for the reference allele. This frameshift mutation introduces a premature stop codon and six novel amino acids (AA), which results in a truncated PRLR protein that lacks 98 C-terminal AA. Conclusions We present the causal mutation for feathering rate in turkey that causes a partial C-terminal loss of the prolactin receptor, and this truncated PRLR protein is strikingly similar to the protein encoded by the slow feathering K allele in chicken. Electronic supplementary material The online version of this article (10.1186/s12711-018-0380-3) contains supplementary material, which is available to authorized users.

Phylogenetic incongruence can be caused by analytical shortcomings or can be the result of biological processes, such as hybridization, incomplete lineage sorting and gene duplication. Differentiation between these causes of incongruence is essential to unravel complex speciation and diversification events. The phylogeny of the True Geese (tribe Anserini, Anatidae, Anseriformes) was, until now, contentious, i.e., the phylogenetic relationships and the timing of divergence between the different goose species could not be fully resolved. We sequenced nineteen goose genomes (representing seventeen species of which three subspecies of the Brent Goose, Branta bernicla) and used an exon-based phylogenomic approach (41,736 exons, representing 5887 genes) to unravel the evolutionary history of this bird group. We thereby provide general guidance on the combination of whole genome evolutionary analyses and analytical tools for such cases where previous attempts to resolve the phylogenetic history of several taxa could not be unravelled. Identical topologies were obtained using either a concatenation (based upon an alignment of 6,630,626 base pairs) or a coalescent-based consensus method. Two major lineages, corresponding to the genera Anser and Branta, were strongly supported. Within the Branta lineage, the White-cheeked Geese form a well-supported sub-lineage that is sister to the Red-breasted Goose (Branta ruficollis). In addition, two main clades of Anser species could be identified, the White Geese and the Grey Geese. The results from the consensus method suggest that the diversification of the genus Anser is heavily influenced by rapid speciation and by hybridization, which may explain the failure of previous studies to resolve the phylogenetic relationships within this genus. The majority of speciation events took place in the late Pliocene and early Pleistocene (between 4 and 2millionyears ago), conceivably driven by a global cooling trend that led to the establishment of a circumpolar tundra belt and the emergence of temperate grasslands. Our approach will be a fruitful strategy for resolving many other complex evolutionary histories at the level of genera, species, and subspecies. published

Aim We compared the power of different nuclear markers to investigate genetic structure of southern Balkan wild boar. We distinguished between historic events, such as isolation in different refugia during glacial periods, from recent demographic processes, such as naturally occurring expansions. Location Southern Balkans/Greece. Methods We sampled 555 wild boars from 20 different locations in southern Balkans/Greece. All individuals were analysed with 10 microsatellites and a subgroup of 91 with 49,508 single nucleotide polymorphisms (SNPs). Patterns of genetic structure and demographic processes were assessed with Bayesian clustering, linkage disequilibrium and past effective population size estimation analysis. Results Both microsatellite and SNP data analyses detected genetic structure caused by historic events and support the existence of three groups in the studied area. A hybrid zone between two of the groups was also detected. We also showed that genome-wide SNP data analysis can identify recent events in bottlenecked populations. Main conclusions We inferred the three groups diverged ~50,000–10,000 yr bp when populations contracted to different refugia. Our findings strengthened the evidence that the southern Balkan area was a glacial refugium including further local smaller refugia. Genome-wide genotyping inferred a recent population expansion that can mimic a ‘refugium within refugium’ scenario. It seems that microsatellite data tend to overestimate genetic structure when genetic drift happens in bottlenecked populations over a short distance. Therefore, genome-wide SNPs are more powerful at inferring phylogeography in natural populations, resolving inconsistencies from mitochondrial and microsatellite data sets.