Introgression and the fate of domesticated genes in a wild mammal population

Selection against domestication alleles in introduced rabbit populations

Humans have moved domestic animals around the globe for thousands of years. These have occasionally established feral populations in nature, often with devastating ecological consequences. To understand how natural selection shapes re-adaptation into the wild, we investigated one of the most successful colonizers in history, the European rabbit. By sequencing the genomes of 297 rabbits across three continents, we show that introduced populations exhibit a mixed wild-domestic ancestry. We show that alleles that increased in frequency during domestication were preferentially selected against in novel natural environments. Interestingly, causative mutations for common domestication traits sometimes segregate at considerable frequencies if associated with less drastic phenotypes (for example, coat colour dilution), whereas mutations that are probably strongly maladaptive in nature are absent. Whereas natural selection largely targeted different genomic regions in each introduced population, some of the strongest signals of parallelism overlap genes associated with neuronal or brain function. This limited parallelism is probably explained by extensive standing genetic variation resulting from domestication together with the complex mixed ancestry of introduced populations. Our findings shed light on the selective and molecular mechanisms that enable domestic animals to re-adapt to the wild and provide important insights for the mitigation and management of invasive populations.

Maintenance of a narrow hybrid zone between native and introduced red foxes (Vulpes vulpes) despite conspecificity and high dispersal capabilities

Human-facilitated introductions of nonnative populations can lead to secondary contact between allopatric lineages, resulting in lineage homogenisation or the formation of stable hybrid zones maintained by reproductive barriers. We investigated patterns of gene flow between the native Sacramento Valley red fox (Vulpes vulpes patwin) and introduced conspecifics of captive-bred origin in California's Central Valley. Considering their recent divergence (20-70 kya), we hypothesised that any observed barriers to gene flow were primarily driven by pre-zygotic (e.g. behavioural differences) rather than post-zygotic (e.g. reduced hybrid fitness) barriers. We also explored whether nonnative genes could confer higher fitness in the human-dominated landscape resulting in selective introgression into the native population. Genetic analysis of red foxes (n = 682) at both mitochondrial (cytochrome b + D-loop) and nuclear (19,051 SNPs) loci revealed narrower cline widths than expected under a simulated model of unrestricted gene flow, consistent with the existence of reproductive barriers. We identified several loci with reduced introgression that were previously linked to behavioural divergence in captive-bred and domestic canids, supporting pre-zygotic, yet possibly hereditary, barriers as a mechanism driving the narrowness and stability of the hybrid zone. Several loci with elevated gene flow from the nonnative into the native population were linked to genes associated with domestication and adaptation to human-dominated landscapes. This study contributes to our understanding of hybridisation dynamics in vertebrates, particularly in the context of species introductions and landscape changes, underscoring the importance of considering how multiple mechanisms may be maintaining lineages at the species and subspecies level.

Introgression from local cultivars is a driver of agricultural adaptation in Argentinian weedy rice

Weedy rice, a pervasive and troublesome weed found across the globe, has often evolved through fertilization of rice cultivars with little importance of crop-weed gene flow. In Argentina, weedy rice has been reported as an important constraint since the early 1970s, and, in the last few years, strains with herbicide-resistance are suspected to evolve. Despite their importance, the origin and genetic composition of Argentinian weedy rice as well its adaptation to agricultural environments has not been explored so far. To study this, we conducted genotyping-by-sequencing on samples of Argentinian weedy and cultivated rice and compared them with published data from weedy, cultivated and wild rice accessions distributed worldwide. In addition, we conducted a phenotypic characterization for weedy-related traits, a herbicide resistance screening and genotyped accessions for known mutations in the acetolactate synthase (ALS) gene, which confers herbicide resistance. Our results revealed large phenotypic variability in Argentinian weedy rice. Most strains were resistant to ALS-inhibiting herbicides with a high frequency of the ALS mutation (A122T) present in Argentinian rice cultivars. Argentinian cultivars belonged to the three major genetic groups of rice: japonica, indica and aus while weeds were mostly aus or aus-indica admixed, resembling weedy rice strains from the Southern Cone region. Phylogenetic analysis supports a single origin for aus-like South American weeds, likely as seed contaminants from the United States, and then admixture with local indica cultivars. Our findings demonstrate that crop to weed introgression can facilitate rapid adaptation to agriculture environments.

Evaluating hybrid speciation and swamping in wild carnivores with a decision-tree approach

Hybridization is an important evolutionary force with a principal role in the origin of new species, known as hybrid speciation. However, ongoing hybridization can create hybrid swamping, in which parental genomes are completely lost. This can become a biodiversity threat if it involves species that have adapted to certain environmental conditions and occur nowhere else. Because conservation scientists commonly have a negative attitude toward hybrids, it is important to improve understanding of the influence of interspecific gene flow on the persistence of species. We reviewed the literature on species hybridization to build a list of all known cases in the order Carnivora. To examine the relative impact, we also noted level of introgression, whether fertile offspring were produced, and whether there was mention of negative or positive evolutionary effects (hybrid speciation and swamping). To evaluate the conservation implications of hybrids, we developed a decision-making tree with which to determine which actions should be taken to manage hybrid species. We found 53 hybrids involving 68 unique taxa, which is roughly 23% of all carnivore species. They mainly involved monophyletic (83%) and sympatric species (75%). For 2 species, the outcome of the assessment was to eliminate or restrict the hybrids: Ethiopian wolf (Canis simensis) and Scottish wildcat (Felis silvestris silvestris). Both species hybridize with their domestic conspecifics. For all other cases, we suggest hybrids be protected in the same manner as native species. We found no evidence of genomic extinction in Carnivora. To the contrary, some species appear to be of hybrid origin, such as the Asiatic black bear (Ursus thibetanus) and African golden wolf (Canis lupaster). Other positive outcomes of hybridization are novel genetic diversity, adaptation to extreme environments, and increased reproductive fitness. These outcomes are particularly valuable for counterbalancing genetic drift and enabling adaptive introgression in a human-dominated world.

Single-cell RNA sequencing and UPHLC-MS/MS targeted metabolomics offer new insights into the etiological basis for male cattle-yak sterility

The variety of species can be efficiently increased by interspecific hybridization. However, because the males in the hybrid progeny are usually sterile, this heterosis cannot be employed when other cattle and yaks are hybridized. While some system-level studies have sought to explore the etiological basis for male cattle-yak sterility, no systematic cellular analyses of this phenomenon have yet been performed. Here, single-cell RNA sequencing and UPHLC-MS/MS targeted metabolomics methods were used to study the differences in testicular tissue between 4-year-old male yak and 4-year-old male cattle-yak, providing new and comprehensive insights into the causes of male cattle-yak sterility. Cattle-yak testes samples detected 6 somatic cell types and one mixed germ cell type. Comparisons of these cell types revealed the more significant differences in Sertoli cells (SCs) and [Leydig cells and myoid cells (LCs_MCs)] between yak and cattle-yak samples compared to other somatic cell clusters. Even though the LCs and MCs from yaks and cattle-yaks were derived from the differentiation of the same progenitor cells, a high degree of overlap between LCs and MCs was observed in yak samples. Still, only a small overlap between LCs and MCs was observed in cattle-yak samples. Functional enrichment analyses revealed that genes down-regulated in cattle-yak SCs were primarily enriched in biological activity, whereas up-regulated genes in these cells were enriched for apoptotic activity. Furthermore, the genes of up-regulated in LCs_MCs of cattle-yak were significantly enriched in enzyme inhibitor and molecular function inhibitor activity. On the other hand, the genes of down-regulated in these cells were enriched for signal receptor binding, molecular function regulation, positive regulation of biological processes, and regulation of cell communication activity. The most significant annotated differences between yak and cattle-yak LCs_MCs were associated with cell-to-cell communication. While yak LCs_MCs regulated spermatogenic cells at spermatogonia, spermatocyte, and spermatid levels, no such relationships were found between cattle-yak LCs_MCs and germ cells. This may suggest that the somatic niche in male cattle-yak testes is a microenvironment that is ultimately not favorable for spermatogenesis.

Genomic prediction in the wild: A case study in Soay sheep

Genomic prediction, the technique whereby an individual's genetic component of their phenotype is estimated from its genome, has revolutionised animal and plant breeding and medical genetics. However, despite being first introduced nearly two decades ago, it has hardly been adopted by the evolutionary genetics community studying wild organisms. Here, genomic prediction is performed on eight traits in a wild population of Soay sheep. The population has been the focus of a >30 year evolutionary ecology study and there is already considerable understanding of the genetic architecture of the focal Mendelian and quantitative traits. We show that the accuracy of genomic prediction is high for all traits, but especially those with loci of large effect segregating. Five different methods are compared, and the two methods that can accommodate zero-effect and large-effect loci in the same model tend to perform best. If the accuracy of genomic prediction is similar in other wild populations, then there is a real opportunity for pedigree-free molecular quantitative genetics research to be enabled in many more wild populations; currently the literature is dominated by studies that have required decades of field data collection to generate sufficiently deep pedigrees. Finally, some of the potential applications of genomic prediction in wild populations are discussed.

Extensive genome introgression between domestic ferret and European polecat during population recovery in Great Britain

The European polecat (Mustela putorius) is a mammalian predator which occurs across much of Europe east to the Ural Mountains. In Great Britain, following years of persecution the range of the European polecat contracted and by the early 1900s was restricted to unmanaged forests of central Wales. The European polecat has recently undergone a population increase due to legal protection and its range now overlaps that of feral domestic ferrets (Mustela putorius furo). During this range expansion, European polecats hybridized with feral domestic ferrets producing viable offspring. Here, we carry out population-level whole-genome sequencing on 8 domestic ferrets, 19 British European polecats, and 15 European polecats from the European mainland. We used a range of population genomics methods to examine the data, including phylogenetics, phylogenetic graphs, model-based clustering, phylogenetic invariants, ABBA-BABA tests, topology weighting, and Fst. We found high degrees of genome introgression in British polecats outside their previous stronghold, even in those individuals phenotyped as “pure” polecats. These polecats ranged from presumed F1 hybrids (gamma = 0.53) to individuals that were much less introgressed (gamma = 0.2). We quantify this introgression and find introgressed genes containing Fst outliers associated with cognitive function and sight.

Genomic architecture of adaptive radiation and hybridization in Alpine whitefish

Adaptive radiations represent some of the most remarkable explosions of diversification across the tree of life. However, the constraints to rapid diversification and how they are sometimes overcome, particularly the relative roles of genetic architecture and hybridization, remain unclear. Here, we address these questions in the Alpine whitefish radiation, using a whole-genome dataset that includes multiple individuals of each of the 22 species belonging to six ecologically distinct ecomorph classes across several lake-systems. We reveal that repeated ecological and morphological diversification along a common environmental axis is associated with both genome-wide allele frequency shifts and a specific, larger effect, locus, associated with the gene edar. Additionally, we highlight the possible role of introgression between species from different lake-systems in facilitating the evolution and persistence of species with unique trait combinations and ecology. These results highlight the importance of both genome architecture and secondary contact with hybridization in fuelling adaptive radiation.

Contemporary selection on MHC genes in a free-living ruminant population

Genes within the major histocompatibility complex (MHC) are the most variable identified in vertebrates. Pathogen-mediated selection is believed to be the main force maintaining MHC diversity. However, relatively few studies have demonstrated contemporary selection on MHC genes. Here, we examine associations between MHC variation and several fitness measurements including total fitness and five fitness components, in 3400 wild Soay sheep (Ovis aries) monitored between 1989 and 2012. In terms of total fitness, measured as lifetime breeding success of all individuals born, we found haplotypes named C and D were associated with decreased and increased male total fitness respectively. In terms of fitness components, juvenile survival was associated with haplotype divergence while individual haplotypes (C, D and F) were associated with adult fitness components. Consistent with the increased male total fitness, the rarest haplotype D has increased in frequency throughout the study period more than expected under neutral expectations. Our results demonstrate contemporary natural selection is acting on MHC class II genes in Soay sheep and the mode of selection on specific fitness components can be different mode from selection on total fitness.

Consequences of Hybridization in Mammals: A Systematic Review

Hybridization, defined as breeding between two distinct taxonomic units, can have an important effect on the evolutionary patterns in cross-breeding taxa. Although interspecific hybridization has frequently been considered as a maladaptive process, which threatens species genetic integrity and survival via genetic swamping and outbreeding depression, in some cases hybridization can introduce novel adaptive variation and increase fitness. Most studies to date focused on documenting hybridization events and analyzing their causes, while relatively little is known about the consequences of hybridization and its impact on the parental species. To address this knowledge gap, we conducted a systematic review of studies on hybridization in mammals published in 2010-2021, and identified 115 relevant studies. Of 13 categories of hybridization consequences described in these studies, the most common negative consequence (21% of studies) was genetic swamping and the most common positive consequence (8%) was the gain of novel adaptive variation. The total frequency of negative consequences (49%) was higher than positive (13%) and neutral (38%) consequences. These frequencies are biased by the detection possibilities of microsatellite loci, the most common genetic markers used in the papers assessed. As negative outcomes are typically easier to demonstrate than positive ones (e.g., extinction vs hybrid speciation), they may be over-represented in publications. Transition towards genomic studies involving both neutral and adaptive variation will provide a better insight into the real impacts of hybridization.

Ancient hybridization patterns between bighorn and thinhorn sheep

Whole-genome sequencing has advanced the study of species evolution, including the detection of genealogical discordant events such as ancient hybridization and incomplete lineage sorting (ILS). The evolutionary history of bighorn (Ovis canadensis) and thinhorn (Ovis dalli) sheep present an ideal system to investigate evolutionary discordance due to their recent and rapid radiation and putative secondary contact between bighorn and thinhorn sheep subspecies, specifically the dark pelage Stone sheep (O. dalli stonei) and predominately white Dall sheep (O. dalli dalli), during the last ice age. Here, we used multiple genomes of bighorn and thinhorn sheep, together with snow (O. nivicola) and the domestic sheep (O. aries) as outgroups, to assess their phylogenomic history, potential introgression patterns and their adaptive consequences. Among the Pachyceriforms (snow, bighorn and thinhorn sheep) a consistent monophyletic species tree was retrieved; however, many genealogical discordance patterns were observed. Alternative phylogenies frequently placed Stone and bighorn as sister clades. This relationship occurred more often and was less divergent than that between Dall and bighorn. We also observed many blocks containing introgression signal between Stone and bighorn genomes in which coat colour genes were present. Introgression signals observed between Dall and bighorn were more random and less frequent, and therefore probably due to ILS or intermediary secondary contact. These results strongly suggest that Stone sheep originated from a complex series of events, characterized by multiple, ancient periods of secondary contact with bighorn sheep.

Introgression of ASIP and TYRP1 Alleles Explains Coat Color Variation in Valais Goats

The Valais Blackneck goat is a Swiss goat breed with a characteristic coat color phenotype. Before the revision of the breed standard in 1938, 4 different color varieties of Valais goats were known. Besides Blackneck animals resembling the modern breed standard, the brown and white Copperneck goat, the white Capra Sempione, and the greyish Grüenochte comprised the historic Valais goats. The brown pigmentation of Copperneck goats had previously been traced back to an introgression of a mutant TYRP1 allele from Toggenburg goats. In the present study, we identified additional introgression events of distinct ASIP alleles causing the remaining 2 rare coat color patterns within the Valais Blackneck goat breed. We identified the introgression of the AWt allele from Appenzell or Saanen goats in white Capra Sempione goats. Similarly, introgression of the Apc allele from Peacock goats resulted in the greyish Grüenochte phenotype. These results demonstrate past hybridization events between breeds that are separated today. A perfect genotype-phenotype association in 393 Valais goats supported the causality of the genotyped variants for the different coat color phenotypes. Our study gives insights into the introgression of functionally relevant copy number variant (CNV) alleles controlling pigmentation between goat breeds with strikingly different coat color patterns.

Genetic architecture and lifetime dynamics of inbreeding depression in a wild mammal

Inbreeding depression is ubiquitous, but we still know little about its genetic architecture and precise effects in wild populations. Here, we combine long-term life-history data with 417 K imputed SNP genotypes for 5952 wild Soay sheep to explore inbreeding depression on a key fitness component, annual survival. Inbreeding manifests in long runs of homozygosity (ROH), which make up nearly half of the genome in the most inbred individuals. The ROH landscape varies widely across the genome, with islands where up to 87% and deserts where only 4% of individuals have ROH. The fitness consequences of inbreeding are severe; a 10% increase in individual inbreeding FROH is associated with a 60% reduction in the odds of survival in lambs, though inbreeding depression decreases with age. Finally, a genome-wide association scan on ROH shows that many loci with small effects and five loci with larger effects contribute to inbreeding depression in survival.

Insight into the Genetic Population Structure of Wild Red Foxes in Poland Reveals Low Risk of Genetic Introgression from Escaped Farm Red Foxes

In this study we assessed the level of genetic introgression between red foxes bred on fur farms in Poland and the native wild population. We also evaluated the impact of a geographic barrier and isolation by distance on gene flow between two isolated subpopulations of the native red fox and their genetic differentiation. Nuclear and mitochondrial DNA was collected from a total of 308 individuals (200 farm and 108 wild red foxes) to study non-native allele flow from farm into wild red fox populations. Genetic structure analyses performed using 24 autosomal microsatellites showed two genetic clusters as being the most probable number of distinct populations. No strong admixture signals between farm and wild red foxes were detected, and significant genetic differentiation was identified between the two groups. This was also apparent from the mtDNA analysis. None of the concatenated haplotypes detected in farm foxes was found in wild animals. The consequence of this was that the haplotype network displayed two genetically distinct groups: farm foxes were completely separated from native ones. Neither the River Vistula nor isolation by distance had a significant impact on gene flow between the separated wild red fox subpopulations. The results of our research indicate a low probability of genetic introgression between farm and native red foxes, and no threat to the genetic integrity of this species.

Getting Back to Nature: Feralization in Animals and Plants

Formerly domesticated organisms and artificially selected genes often escape controlled cultivation, but their subsequent evolution is not well studied. In this review, we examine plant and animal feralization through an evolutionary lens, including how natural selection, artificial selection, and gene flow shape feral genomes, traits, and fitness. Available evidence shows that feralization is not a mere reversal of domestication. Instead, it is shaped by the varied and complex histories of feral populations, and by novel selection pressures. To stimulate further insight we outline several future directions. These include testing how 'domestication genes' act in wild settings, studying the brains and behaviors of feral animals, and comparative analyses of feral populations and taxa. This work offers feasible and exciting research opportunities with both theoretical and practical applications.

Introgression drives repeated evolution of winter coat color polymorphism in hares

Changing from summer-brown to winter-white pelage or plumage is a crucial adaptation to seasonal snow in more than 20 mammal and bird species. Many of these species maintain nonwhite winter morphs, locally adapted to less snowy conditions, which may have evolved independently. Mountain hares (Lepus timidus) from Fennoscandia were introduced into the Faroe Islands in 1855. While they were initially winter-white, within ∼65 y all Faroese hares became winter-gray, a morph that occurs in the source population at low frequency. The documented population history makes this a valuable model for understanding the genetic basis and evolution of the seasonal trait polymorphism. Through whole-genome scans of differentiation and single-nucleotide polymorphism (SNP) genotyping, we associated winter coat color polymorphism to the genomic region of the pigmentation gene Agouti, previously linked to introgression-driven winter coat color variation in the snowshoe hare (Lepus americanus). Lower Agouti expression in the skin of winter-gray individuals during the autumn molt suggests that regulatory changes may underlie the color polymorphism. Variation in the associated genomic region shows signatures of a selective sweep in the Faroese population, suggesting that positive selection drove the fixation of the variant after the introduction. Whole-genome analyses of several hare species revealed that the winter-gray variant originated through introgression from a noncolor changing species, in keeping with the history of ancient hybridization between the species. Our findings show the recurrent role of introgression in generating winter coat color variation by repeatedly recruiting the regulatory region of Agouti to modulate seasonal coat color change.

The genetic architecture of helminth-specific immune responses in a wild population of Soay sheep (Ovis aries)

Much of our knowledge of the drivers of immune variation, and how these responses vary over time, comes from humans, domesticated livestock or laboratory organisms. While the genetic basis of variation in immune responses have been investigated in these systems, there is a poor understanding of how genetic variation influences immunity in natural, untreated populations living in complex environments. Here, we examine the genetic architecture of variation in immune traits in the Soay sheep of St Kilda, an unmanaged population of sheep infected with strongyle gastrointestinal nematodes. We assayed IgA, IgE and IgG antibodies against the prevalent nematode Teladorsagia circumcincta in the blood plasma of > 3,000 sheep collected over 26 years. Antibody levels were significantly heritable (h² = 0.21 to 0.57) and highly stable over an individual’s lifespan. IgA levels were strongly associated with a region on chromosome 24 explaining 21.1% and 24.5% of heritable variation in lambs and adults, respectively. This region was adjacent to two candidate loci, Class II Major Histocompatibility Complex Transactivator (CIITA) and C-Type Lectin Domain Containing 16A (CLEC16A). Lamb IgA levels were also associated with the immunoglobulin heavy constant loci (IGH) complex, and adult IgE levels and lamb IgA and IgG levels were associated with the major histocompatibility complex (MHC). This study provides evidence of high heritability of a complex immunological trait under natural conditions and provides the first evidence from a genome-wide study that large effect genes located outside the MHC region exist for immune traits in the wild.

Understanding how immune responses vary in natural populations can give an insight into how infection affects the ability of hosts and parasites to survive and reproduce, and how this drives evolutionary and ecological dynamics. Yet, very little is known about how immune responses vary over an individual’s lifetime and how genes contribute to this variation under natural conditions. Our study investigates the genetic architecture of variation in three antibody types, IgA, IgE and IgG in a wild population of Soay sheep on the St Kilda archipelago in North-West Scotland. Using data collected over 26 years, we show that antibody levels have a heritable basis in lambs and adults and are stable over an individual’s lifetime. We also identify several genomic regions with large effects on immune responses. Our study offers the first insights into the genetic control of immunity in a wild population, which is essential to understand how immune profiles vary in challenging natural conditions and how natural selection maintains genetic variation in complex immune traits.

Genomic insights into the vulnerability of sympatric whitefish species flocks

The erosion of habitat heterogeneity can reduce species diversity directly but can also lead to the loss of distinctiveness of sympatric species through speciation reversal. We know little about changes in genomic differentiation during the early stages of these processes, which can be mediated by anthropogenic perturbation. Here, we analyse three sympatric whitefish species (Coregonus spp) sampled across two neighbouring and connected Swiss pre-alpine lakes, which have been differentially affected by anthropogenic eutrophication. Our data set comprises 16,173 loci genotyped across 138 whitefish using restriction-site associated DNA sequencing (RADseq). Our analysis suggests that in each of the two lakes, the population of a different, but ecologically similar, whitefish species declined following a recent period of eutrophication. Genomic signatures consistent with hybridization are more pronounced in the more severely impacted lake. Comparisons between sympatric pairs of whitefish species with contrasting ecology, where one is shallow benthic and the other one more profundal pelagic, reveal genomic differentiation that is largely correlated along the genome, while differentiation is uncorrelated between pairs of allopatric provenance with similar ecology. We identify four genomic loci that provide evidence of parallel divergent adaptation between the shallow benthic species and the two different more profundal species. Functional annotations available for two of those loci are consistent with divergent ecological adaptation. Our genomic analysis indicates the action of divergent natural selection between sympatric whitefish species in pre-alpine lakes and reveals the vulnerability of these species to anthropogenic alterations of the environment and associated adaptive landscape.

How to estimate kinship

The concept of kinship permeates many domains of fundamental and applied biology ranging from social evolution to conservation science to quantitative and human genetics. Until recently, pedigrees were the gold standard to infer kinship, but the advent of next‐generation sequencing and the availability of dense genetic markers in many species make it a good time to (re)evaluate the usefulness of genetic markers in this context. Using three published data sets where both pedigrees and markers are available, we evaluate two common and a new genetic estimator of kinship. We show discrepancies between pedigree values and marker estimates of kinship and explore via simulations the possible reasons for these. We find these discrepancies are attributable to two main sources: pedigree errors and heterogeneity in the origin of founders. We also show that our new marker‐based kinship estimator has very good statistical properties and behaviour and is particularly well suited for situations where the source population is of small size, as will often be the case in conservation biology, and where high levels of kinship are expected, as is typical in social evolution studies.

Widespread, long-term admixture between grey wolves and domestic dogs across Eurasia and its implications for the conservation status of hybrids

Hybridisation between a domesticated species and its wild ancestor is an important conservation problem, especially if it results in the introgression of domestic gene variants into wild species. Nevertheless, the legal status of hybrids remains unregulated, partially because of the limited understanding of the hybridisation process and its consequences. The occurrence of hybridisation between grey wolves and domestic dogs is well documented from different parts of the wolf geographic range, but little is known about the frequency of hybridisation events, their causes and the genetic impact on wolf populations. We analysed 61K SNPs spanning the canid genome in wolves from across Eurasia and North America and compared that data to similar data from dogs to identify signatures of admixture. The haplotype block analysis, which included 38 autosomes and the X chromosome, indicated the presence of individuals of mixed wolf-dog ancestry in most Eurasian wolf populations, but less admixture was present in North American populations. We found evidence for male-biased introgression of dog alleles into wolf populations, but also identified a first-generation hybrid resulting from mating between a female dog and a male wolf. We found small blocks of dog ancestry in the genomes of 62% Eurasian wolves studied and melanistic individuals with no signs of recent admixed ancestry, but with a dog-derived allele at a locus linked to melanism. Consequently, these results suggest that hybridisation has been occurring in different parts of Eurasia on multiple timescales and is not solely a recent phenomenon. Nevertheless, wolf populations have maintained genetic differentiation from dogs, suggesting that hybridisation at a low frequency does not diminish distinctiveness of the wolf gene pool. However, increased hybridisation frequency may be detrimental for wolf populations, stressing the need for genetic monitoring to assess the frequency and distribution of individuals resulting from recent admixture.

Low rates of hybridization between European wildcats and domestic cats in a human-dominated landscape

Hybridization between wild species and their domestic congeners is considered a major threat for wildlife conservation. Genetic integrity of the European wildcat, for instance, is a concern as they are outnumbered by domestic cats by several orders of magnitude throughout its range. We genotyped 1,071 individual wildcat samples obtained from hair traps and roadkills collected across the highly fragmented forests of western Central Europe, in Germany and Luxembourg, to assess domestic cat introgression in wildcats in human-dominated landscapes. Analyses using a panel of 75 autosomal SNPs suggested a low hybridization rate, with 3.5% of wildcat individuals being categorized as F1, F2, or backcrosses to either parental taxon. We report that results based on a set of SNPs were more consistent than on a set of 14 microsatellite markers, showed higher accuracy to detect hybrids and their class in simulation analyses, and were less affected by underlying population structure. Our results strongly suggest that very high hybridization rates previously reported for Central Europe may be partly due to inadequate choice of markers and/or sampling design. Our study documents that an adequately selected SNP panel for hybrid detection may be used as an alternative to commonly applied microsatellite markers, including studies relying on noninvasively collected samples. In addition, our finding of overall low hybridization rates in Central European wildcats provides an example of successful wildlife coexistence in human-dominated, fragmented landscapes.

Evidence for Selection-by-Environment but Not Genotype-by-Environment Interactions for Fitness-Related Traits in a Wild Mammal Population

How do environmental conditions influence selection and genetic variation in wild populations? There is widespread evidence for selection-by-environment interactions (S*E), but we reviewed studies of natural populations estimating the extent of genotype-by-environment interactions (G*E) in response to natural variation in environmental conditions and found that evidence for G*E appears to be rare within single populations in the wild. Studies estimating the simultaneous impact of environmental variation on both selection and genetic variation are especially scarce. Here, we used 24 years of data collected from a wild Soay sheep population to quantify how an important environmental variable, population density, impacts upon (1) selection through annual contribution to fitness and (2) expression of genetic variation, in six morphological and life history traits: body weight, hind leg length, parasite burden, horn length, horn growth, and testicular circumference. Our results supported the existence of S*E: selection was stronger in years of higher population density for all traits apart from horn growth, with directional selection being stronger under more adverse conditions. Quantitative genetic models revealed significant additive genetic variance for body weight, leg length, parasite burden, horn length, and testes size, but not for horn growth or our measure of annual fitness. However, random regression models found variation between individuals in their responses to the environment in only three traits, and did not support the presence of G*E for any trait. Our analyses of St Kilda Soay sheep data thus concurs with our cross-study review that, while natural environmental variation within a population can profoundly alter the strength of selection on phenotypic traits, there is less evidence for its effect on the expression of genetic variance in the wild.

Human-mediated secondary contact of two tortoise lineages results in sex-biased introgression

Human-mediated secondary contact of recently diverged taxa offers valuable opportunities for studying the evolutionary mechanisms involved in the establishment and maintenance of genetic boundaries between taxa. We used mitochondrial and microsatellite markers to examine a recently introduced population of the spur-thighed tortoise (Testudo graeca) of mixed origin in the Doñana National Park (SW Spain). The earliest records of tortoises in Doñana trace back to the 18th century, but several population reinforcements in the 20th century with animals from Morocco are well-documented. Consequently, different genetic lineages, which represent distinct subspecies, are thought to co-exist there. Our results confirmed the presence of distinct lineages by revealing that tortoises of the subspecies T. g. marokkensis were introduced into a local allochthonous T. g. graeca population. Unexpectedly, T. g. marokkensis haplotypes exclusively appeared in males, and admixture levels were statistically sex-biased toward males. The sex ratio of the population deviated from parity, with males being 2.36-fold more abundant than females. Our results indicated that population reinforcements had a strong effect on the genetic composition of this population and aggravated its sex ratio deviation. We predict that this sex-biased pattern of introgression is ephemeral and advocated to the near loss of T. g. marokkensis haplotypes.

Inbreeding depression by environment interactions in a free-living mammal population

Experimental studies often find that inbreeding depression is more severe in harsh environments, but the few studies of in situ wild populations available to date rarely find strong support for this effect. We investigated evidence for inbreeding depression by environment interactions in nine traits in the individually monitored Soay sheep population of St Kilda, using genomic inbreeding coefficients based on 37 037 single-nucleotide polymorphism loci, and population density as an axis of environmental variation. All traits showed variation with population density and all traits showed some evidence for depression because of either an individual's own inbreeding or maternal inbreeding. However, only six traits showed evidence for an interaction in the expected direction, and only two interactions were statistically significant. We identify three possible reasons why wild population studies may generally fail to find strong support for interactions between inbreeding depression and environmental variation compared with experimental studies. First, for species with biparental inbreeding only, the amount of observed inbreeding in natural populations is generally low compared with that used in experimental studies. Second, it is possible that experimental studies sometimes actually impose higher levels of stress than organisms experience in the wild. Third, some purging of the deleterious recessive alleles that underpin interaction effects may occur in the wild.

Genomic analysis reveals depression due to both individual and maternal inbreeding in a free-living mammal population

There is ample evidence for inbreeding depression manifested as a reduction in fitness or fitness-related traits in the focal individual. In many organisms, fitness is not only affected by genes carried by the individual, but also by genes carried by their parents, for example if receiving parental care. While maternal effects have been described in many systems, the extent to which inbreeding affects fitness directly through the focal individual, or indirectly through the inbreeding coefficients of its parents, has rarely been examined jointly. The Soay sheep study population is an excellent system in which to test for both effects, as lambs receive extended maternal care. Here, we tested for both maternal and individual inbreeding depression in three fitness-related traits (birthweight and weight and hindleg length at 4 months of age) and three fitness components (first-year survival, adult annual survival and annual breeding success), using either pedigree-derived inbreeding or genomic estimators calculated using ~37 000 SNP markers. We found evidence for inbreeding depression in 4-month hindleg and weight, first-year survival in males, and annual survival and breeding success in adults. Maternal inbreeding was found to depress both birthweight and 4-month weight. We detected more instances of significant inbreeding depression using genomic estimators than the pedigree, which is partly explained through the increased sample sizes available. In conclusion, our results highlight that cross-generational inbreeding effects warrant further exploration in species with parental care and that modern genomic tools can be used successfully instead of, or alongside, pedigrees in natural populations.

Conserved Genetic Architecture Underlying Individual Recombination Rate Variation in a Wild Population of Soay Sheep (Ovis aries)

Meiotic recombination breaks down linkage disequilibrium (LD) and forms new haplotypes, meaning that it is an important driver of diversity in eukaryotic genomes. Understanding the causes of variation in recombination rate is important in interpreting and predicting evolutionary phenomena and in understanding the potential of a population to respond to selection. However, despite attention in model systems, there remains little data on how recombination rate varies at the individual level in natural populations. Here we used extensive pedigree and high-density SNP information in a wild population of Soay sheep (Ovis aries) to investigate the genetic architecture of individual autosomal recombination rates. Individual rates were high relative to other mammal systems and were higher in males than in females (autosomal map lengths of 3748 and 2860 cM, respectively). The heritability of autosomal recombination rate was low but significant in both sexes (h² = 0.16 and 0.12 in females and males, respectively). In females, 46.7% of the heritable variation was explained by a subtelomeric region on chromosome 6; a genome-wide association study showed the strongest associations at locus RNF212, with further associations observed at a nearby ∼374-kb region of complete LD containing three additional candidate loci, CPLX1, GAK, and PCGF3. A second region on chromosome 7 containing REC8 and RNF212B explained 26.2% of the heritable variation in recombination rate in both sexes. Comparative analyses with 40 other sheep breeds showed that haplotypes associated with recombination rates are both old and globally distributed. Both regions have been implicated in rate variation in mice, cattle, and humans, suggesting a common genetic architecture of recombination rate variation in mammals.

Asynchrony of senescence among phenotypic traits in a wild mammal population

The degree to which changes in lifespan are coupled to changes in senescence in different physiological systems and phenotypic traits is a central question in biogerontology. It is underpinned by deeper biological questions about whether or not senescence is a synchronised process, or whether levels of synchrony depend on species or environmental context. Understanding how natural selection shapes patterns of synchrony in senescence across physiological systems and phenotypic traits demands the longitudinal study of many phenotypes under natural conditions. Here, we examine the patterns of age-related variation in late adulthood in a wild population of Soay sheep (Ovis aries) that have been the subject of individual-based monitoring for thirty years. We examined twenty different phenotypic traits in both males and females, encompassing vital rates (survival and fecundity), maternal reproductive performance (offspring birth weight, birth date and survival), male rutting behaviour, home range measures, parasite burdens, and body mass. We initially quantified age-related variation in each trait having controlled for annual variation in the environment, among-individual variation and selective disappearance effects. We then standardised our age-specific trait means and tested whether age trajectories could be meaningfully grouped according to sex or the type of trait. Whilst most traits showed age-related declines in later life, we found striking levels of asynchrony both within and between the sexes. Of particular note, female fecundity and reproductive performance declined with age, but male annual reproductive success did not. We also discovered that whilst home range size and quality decline with age in females, home range size increases with age in males. Our findings highlight the complexity of phenotypic ageing under natural conditions and, along with emerging data from other wild populations and laboratory models, suggest that the long-standing hypothesis within evolutionary biology that fitness-related traits should senesce in a synchronous manner is seriously flawed.

Population structure and history of the Welsh sheep breeds determined by whole genome genotyping

Background

One of the most economically important areas within the Welsh agricultural sector is sheep farming, contributing around £230 million to the UK economy annually. Phenotypic selection over several centuries has generated a number of native sheep breeds, which are presumably adapted to the diverse and challenging landscape of Wales. Little is known about the history, genetic diversity and relationships of these breeds with other European breeds. We genotyped 353 individuals from 18 native Welsh sheep breeds using the Illumina OvineSNP50 array and characterised the genetic structure of these breeds. Our genotyping data were then combined with, and compared to, those from a set of 74 worldwide breeds, previously collected during the International Sheep Genome Consortium HapMap project.

Results

Model based clustering of the Welsh and European breeds indicated shared ancestry. This finding was supported by multidimensional scaling analysis (MDS), which revealed separation of the European, African and Asian breeds. As expected, the commercial Texel and Merino breeds appeared to have extensive co-ancestry with most European breeds. Consistently high levels of haplotype sharing were observed between native Welsh and other European breeds. The Welsh breeds did not, however, form a genetically homogeneous group, with pairwise F_ST between breeds averaging 0.107 and ranging between 0.020 and 0.201. Four subpopulations were identified within the 18 native breeds, with high homogeneity observed amongst the majority of mountain breeds. Recent effective population sizes estimated from linkage disequilibrium ranged from 88 to 825.

Conclusions

Welsh breeds are highly diverse with low to moderate effective population sizes and form at least four distinct genetic groups. Our data suggest common ancestry between the native Welsh and European breeds. These findings provide the basis for future genome-wide association studies and a first step towards developing genomics assisted breeding strategies in the UK.

Electronic supplementary material

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

Heterogeneity of genetic architecture of body size traits in a free-living population

Knowledge of the underlying genetic architecture of quantitative traits could aid in understanding how they evolve. In wild populations, it is still largely unknown whether complex traits are polygenic or influenced by few loci with major effect, due to often small sample sizes and low resolution of marker panels. Here, we examine the genetic architecture of five adult body size traits in a free-living population of Soay sheep on St Kilda using 37 037 polymorphic SNPs. Two traits (jaw and weight) show classical signs of a polygenic trait: the proportion of variance explained by a chromosome was proportional to its length, multiple chromosomes and genomic regions explained significant amounts of phenotypic variance, but no SNPs were associated with trait variance when using GWAS. In comparison, genetic variance for leg length traits (foreleg, hindleg and metacarpal) was disproportionately explained by two SNPs on chromosomes 16 (s23172.1) and 19 (s74894.1), which each explained >10% of the additive genetic variance. After controlling for environmental differences, females heterozygous for s74894.1 produced more lambs and recruits during their lifetime than females homozygous for the common allele conferring long legs. We also demonstrate that alleles conferring shorter legs have likely entered the population through a historic admixture event with the Dunface sheep. In summary, we show that different proxies for body size can have very different genetic architecture and that dense SNP helps in understanding both the mode of selection and the evolutionary history at loci underlying quantitative traits in natural populations.

A standardized method for quantifying unidirectional genetic introgression

Genetic introgression of domesticated to wild conspecifics is of great concern to the genetic integrity and viability of the wild populations. Therefore, we need tools that can be used for monitoring unidirectional gene flow from domesticated to wild populations. A challenge to quantitation of unidirectional gene flow is that both the donor and the recipient population may be genetically substructured and that the subpopulations are subjected to genetic drift and may exchange migrants between one another. We develop a standardized method for quantifying and monitoring domesticated to wild gene flow and demonstrate its usefulness to farm and wild Atlantic salmon as a model species. The challenge of having several wild and farm populations was circumvented by in silico generating one analytical center point for farm and wild salmon, respectively. Distributions for the probability that an individual is wild were generated from individual-based analyses of observed wild and farm genotypes using STRUCTURE. We show that estimates of proportions of the genome being of domesticated origin in a particular wild population can be obtained without having a historical reference sample for the same population. The main advantages of the method presented are the standardized way in which genetic processes within and between populations are taken into account, and the individual-based analyses giving estimates for each individual independent of other individuals. The method makes use of established software, and as long as genetic markers showing generic genetic differences between domesticated and wild populations are available, it can be applied to all species with unidirectional gene flow. Results from our method are easy to interpret and understand, and will serve as a powerful tool for management, especially because there is no need for a specific historical wild reference sample.

Linkage disequilibrium over short physical distances measured in sheep using a high-density SNP chip

The extent of linkage disequilibrium (LD) between genetic loci has implications for both association studies and the accuracy of genomic prediction. To characterise the persistence of LD in diverse sheep breeds, two SNP genotyping platforms were used. First, existing SNP genotypes from 63 breeds obtained using the ovine SNP50 BeadChip (49,034 loci) were used to estimate LD decay in populations with contrasting levels of genetic diversity. Given the paucity of marker pairs separated by short physical distances on the SNP50 BeadChip, genotyping was subsequently performed for four breeds using the recently developed ovine HD BeadChip that assays approximately 600,000 SNPs with an average genomic spacing of 5 kb. This facilitated a highly accurate estimate of LD over short genomic distances (<30 kb) and revealed LD varies considerably between sheep breeds. Further, sheep appear to contain generally lower levels of LD than do other domestic species, likely a reflection of aspects of their past population history.

Introgression from domestic goat generated variation at the major histocompatibility complex of Alpine ibex

The major histocompatibility complex (MHC) is a crucial component of the vertebrate immune system and shows extremely high levels of genetic polymorphism. The extraordinary genetic variation is thought to be ancient polymorphisms maintained by balancing selection. However, introgression from related species was recently proposed as an additional mechanism. Here we provide evidence for introgression at the MHC in Alpine ibex (Capra ibex ibex). At a usually very polymorphic MHC exon involved in pathogen recognition (DRB exon 2), Alpine ibex carried only two alleles. We found that one of these DRB alleles is identical to a DRB allele of domestic goats (Capra aegagrus hircus). We sequenced 2489 bp of the coding and non-coding regions of the DRB gene and found that Alpine ibex homozygous for the goat-type DRB exon 2 allele showed nearly identical sequences (99.8%) to a breed of domestic goats. Using Sanger and RAD sequencing, microsatellite and SNP chip data, we show that the chromosomal region containing the goat-type DRB allele has a signature of recent introgression in Alpine ibex. A region of approximately 750 kb including the DRB locus showed high rates of heterozygosity in individuals carrying one copy of the goat-type DRB allele. These individuals shared SNP alleles both with domestic goats and other Alpine ibex. In a survey of four Alpine ibex populations, we found that the region surrounding the DRB allele shows strong linkage disequilibria, strong sequence clustering and low diversity among haplotypes carrying the goat-type allele. Introgression at the MHC is likely adaptive and introgression critically increased MHC DRB diversity in the genetically impoverished Alpine ibex. Our finding contradicts the long-standing view that genetic variability at the MHC is solely a consequence of ancient trans-species polymorphism. Introgression is likely an underappreciated source of genetic diversity at the MHC and other loci under balancing selection.

The major histocompatibility complex (MHC), a crucial component of the defense against pathogens, contains the most polymorphic functional genes in vertebrate genomes. The extraordinary genetic variation is generally considered to be ancient. We investigated whether a previously neglected mechanism, introgression from related species, provides an additional source of MHC variation. We show that introgression from domestic goat dramatically increased genetic variation at the MHC of Alpine ibex, a species that had nearly gone extinct during the 18^th century, but has been restored to large numbers since. We show that Alpine ibex share one of only two alleles at a generally highly polymorphic MHC locus with domestic goats and that the chromosomal region containing the goat-type allele has a signature of recent introgression. Our finding contradicts the long-standing view that ancient trans-species polymorphism is the sole source of the extraordinary genetic variability at the MHC. Instead, we show that in Alpine ibex introgression generated genetic diversity at a MHC locus. Our study supports the view that loci favoring genetic polymorphism may be susceptible to adaptive introgression from related species and will encourage future research to identify unexpected signatures of introgression.