Wild pedigrees: the way forward

Runs of homozygosity reveal contrasting histories of inbreeding across global lineages of the edible porcini mushroom, Boletus edulis

Inbreeding, the mating of individuals that are related through common ancestry, is of central importance in evolutionary and conservation biology due to its impacts on individual fitness and population dynamics. However, while advanced genomic approaches have revolutionised the study of inbreeding in animals, genomic studies of inbreeding are rare in plants and lacking in fungi. We investigated global patterns of inbreeding in the prized edible porcini mushroom Boletus edulis using 225 whole genomes from seven lineages distributed across the northern hemisphere. Genomic inbreeding was quantified using runs of homozygosity (ROHs). We found appreciable variation both among and within lineages, with some individuals having over 20% of their genomes in ROHs. Much of this variation could be explained by a combination of elevation and latitude, and to a lesser extent by predicted habitat suitability during the last glacial maximum. In line with this, the majority of ROHs were short, reflecting ancient common ancestry dating back approximately 200-1700 generations ago, while longer ROHs indicative of recent common ancestry (less than approximately 50 generations ago) were infrequent. Our study reveals the inbreeding legacy of major climatic events in a widely distributed forest mutualist, aligning with prevailing theories and empirical studies of the impacts of historical glaciation events on the dominant forest tree species of the northern hemisphere.

When birds of a feather flock together: Severe genomic erosion and the implications for genetic rescue in an endangered island passerine

The Seychelles magpie-robin's (SMR) five island populations exhibit some of the lowest recorded levels of genetic diversity among endangered birds, and high levels of inbreeding. These populations collapsed during the 20th century, and the species was listed as Critically Endangered in the IUCN Red List in 1994. An assisted translocation-for-recovery program initiated in the 1990s increased the number of mature individuals, resulting in its downlisting to Endangered in 2005. Here, we explore the temporal genomic erosion of the SMR based on a dataset of 201 re-sequenced whole genomes that span the past ~150 years. Our sample set includes individuals that predate the bottleneck by up to 100 years, as well as individuals from contemporary populations established during the species recovery program. Despite the SMR's recent demographic recovery, our data reveal a marked increase in both the genetic load and realized load in the extant populations when compared to the historical samples. Conservation management may have reduced the intensity of selection by increasing juvenile survival and relaxing intraspecific competition between individuals, resulting in the accumulation of loss-of-function mutations (i.e. severely deleterious variants) in the rapidly recovering population. In addition, we found a 3-fold decrease in genetic diversity between temporal samples. While the low genetic diversity in modern populations may limit the species' adaptability to future environmental changes, future conservation efforts (including IUCN assessments) may also need to assess the threats posed by their high genetic load. Our computer simulations highlight the value of translocations for genetic rescue and show how this could halt genomic erosion in threatened species such as the SMR.

Benchmarking kinship estimation tools for ancient genomes using pedigree simulations

There is growing interest in uncovering genetic kinship patterns in past societies using low-coverage palaeogenomes. Here, we benchmark four tools for kinship estimation with such data: lcMLkin, NgsRelate, KIN, and READ, which differ in their input, IBD estimation methods, and statistical approaches. We used pedigree and ancient genome sequence simulations to evaluate these tools when only a limited number (1 to 50 K, with minor allele frequency ≥0.01) of shared SNPs are available. The performance of all four tools was comparable using ≥20 K SNPs. We found that first-degree related pairs can be accurately classified even with 1 K SNPs, with 85% F1 scores using READ and 96% using NgsRelate or lcMLkin. Distinguishing third-degree relatives from unrelated pairs or second-degree relatives was also possible with high accuracy (F1 > 90%) with 5 K SNPs using NgsRelate and lcMLkin, while READ and KIN showed lower success (69 and 79% respectively). Meanwhile, noise in population allele frequencies and inbreeding (first-cousin mating) led to deviations in kinship coefficients, with different sensitivities across tools. We conclude that using multiple tools in parallel might be an effective approach to achieve robust estimates on ultra-low-coverage genomes.

Little evidence of inbreeding depression for birth mass, survival and growth in Antarctic fur seal pups

Inbreeding depression, the loss of offspring fitness due to consanguineous mating, is generally detrimental for individual performance and population viability. We investigated inbreeding effects in a declining population of Antarctic fur seals (Arctocephalus gazella) at Bird Island, South Georgia. Here, localised warming has reduced the availability of the seal's staple diet, Antarctic krill, leading to a temporal increase in the strength of selection against inbred offspring, which are increasingly failing to recruit into the adult breeding population. However, it remains unclear whether selection operates before or after nutritional independence at weaning. We therefore used microsatellite data from 885 pups and their mothers, and SNP array data from 98 mother-offspring pairs, to quantify the effects of individual and maternal inbreeding on three important neonatal fitness traits: birth mass, survival and growth. We did not find any clear or consistent effects of offspring or maternal inbreeding on any of these traits. This suggests that selection filters inbred individuals out of the population as juveniles during the time window between weaning and recruitment. Our study brings into focus a poorly understood life-history stage and emphasises the importance of understanding the ecology and threats facing juvenile pinnipeds.

Genetic architecture of inbreeding depression may explain its persistence in a population of wild red deer

Inbreeding depression is of major concern in declining populations, but relatively little is known about its genetic architecture in wild populations, such as the degree to which it is composed of large or small effect loci and their distribution throughout the genome. Here, we combine fitness and genomic data from a wild population of red deer to investigate the genomic distribution of inbreeding effects. Based on the runs of homozygosity (ROH)-based inbreeding coefficient, FROH, we use chromosome-specific inbreeding coefficients (FROHChr) to explore whether the effect of inbreeding varies between chromosomes. Under the assumption that within an individual the probability of being identical-by-descent is equal across all chromosomes, we used a multi-membership model to estimate the deviation of FROHChr from the average inbreeding effect. This novel approach ensures effect sizes are not overestimated whilst maximising the power of our available dataset of >3000 individuals genotyped on >35,000 autosomal SNPs. We find that most chromosomes confer a minor reduction in fitness-related traits, which when these effects are summed, results in the observed inbreeding depression in birth weight, survival and lifetime breeding success. However, no chromosomes had a significant detrimental effect compared to the overall effect of inbreeding, indicating no major effect loci. We conclude that in this population, inbreeding depression is likely the result of multiple mildly or moderately deleterious mutations spread across all chromosomes, which are difficult to detect with statistical confidence. Such mutations will be inefficiently purged, which may explain the persistence of inbreeding depression in this population.

Polly: An R package for genotyping microsatellites and detecting highly polymorphic DNA markers from short-read data

Highly polymorphic markers, such as microsatellites, are invaluable for the study of natural populations. However, contemporary methods for genotyping highly polymorphic variants have serious drawbacks that impede their efficiency. We created Polly, an R package with C++ source code that uses Illumina short-read data to genotype microsatellites, detect highly polymorphic variants and identify clusters of highly polymorphic SNPs, indels and microsatellites. We tested Polly on short-read data from Xiphophorus birchmanni (Teleostei: Poeciliidae) and Arabidopsis thaliana, finding it to be efficient and accurate both for microsatellite genotyping and polymorphic marker detection. This program can be applied to any diploid population for which there exists short-read data and at least one scaffolded reference genome.

The estimation of additive genetic variance of body size in a wild passerine is sensitive to the method used to estimate relatedness among the individuals

Assessing additive genetic variance is a crucial step in predicting the evolutionary response of a target trait. However, the estimated genetic variance may be sensitive to the methodology used, e.g., the way relatedness is assessed among the individuals, especially in wild populations where social pedigrees can be inaccurate. To investigate this possibility, we investigated the additive genetic variance in tarsus length, a major proxy of skeletal body size in birds. The model species was the collared flycatcher (Ficedula albicollis), a socially monogamous but genetically polygamous migratory passerine. We used two relatedness matrices to estimate the genetic variance: (1) based solely on social links and (2) a genetic similarity matrix based on a large array of single-nucleotide polymorphisms (SNPs). Depending on the relatedness matrix considered, we found moderate to high additive genetic variance and heritability estimates for tarsus length. In particular, the heritability estimates were higher when obtained with the genetic similarity matrix instead of the social pedigree. Our results confirm the potential for this crucial trait to respond to selection and highlight methodological concerns when calculating additive genetic variance and heritability in phenotypic traits. We conclude that using a social pedigree instead of a genetic similarity matrix to estimate relatedness among individuals in a genetically polygamous wild population may significantly deflate the estimates of additive genetic variation.

Taking identity-by-descent analysis into the wild: Estimating realized relatedness in free-ranging macaques

Biological relatedness is a key consideration in studies of behavior, population structure, and trait evolution. Except for parent-offspring dyads, pedigrees capture relatedness imperfectly. The number and length of DNA segments that are identical-by-descent (IBD) yield the most precise estimates of relatedness. Here, we leverage novel methods for estimating locus-specific IBD from low coverage whole genome resequencing data to demonstrate the feasibility and value of resolving fine-scaled gradients of relatedness in free-living animals. Using primarily 4-6× coverage data from a rhesus macaque (Macaca mulatta) population with available long-term pedigree data, we show that we can call the number and length of IBD segments across the genome with high accuracy even at 0.5× coverage. The resulting estimates demonstrate substantial variation in genetic relatedness within kin classes, leading to overlapping distributions between kin classes. They identify cryptic genetic relatives that are not represented in the pedigree and reveal elevated recombination rates in females relative to males, which allows us to discriminate maternal and paternal kin using genotype data alone. Our findings represent a breakthrough in the ability to understand the predictors and consequences of genetic relatedness in natural populations, contributing to our understanding of a fundamental component of population structure in the wild.

Temporal Changes in Demography and Genetic Diversity of Oriental Storks at the Stage of Long-Term Captive Propagation and Reintroduction Initiation

Captive propagation and reintroduction are the major steps in the ex-situ conservation of locally extirpated endangered species in a historical region. In a species restoration project conducted in South Korea, we examined temporal changes in demographics and genetic diversity of oriental storks (Ciconia boyciana). Demographic and genetic data from 1996-2018 were analyzed for 80% of all captive and recently reintroduced individuals. Founder establishment and pair formation induced increases in population size and genetic diversity during the early stage of captive propagation. The degree of genetic diversity was found to become saturated and stable with long-term captive propagation. However, this might be a concern for future genetic diversity of both captive and reintroduced populations simultaneously due to the extraction of captive populations at the early stage of reintroduction. Our findings suggest that periodic evaluation of genetic diversity and selection for releasing individuals, using effective genetic markers, would assist in balancing the genetic diversity of the captive and reintroduced oriental storks at the early stage of reintroduction.

An empirical evaluation of the estimation of inbreeding depression from molecular markers under suboptimal conditions

Inbreeding depression (ID), the reduction in fitness due to inbreeding, is typically measured by the regression of the phenotypic values of individuals for a particular trait on their corresponding inbreeding coefficients (F). While genealogical records can provide these coefficients, they may be unavailable or incomplete, making molecular markers a useful alternative. The power to detect ID and its accuracy depend on the variation of F values of individuals, the sample sizes available, and the accuracy in the estimation of individual fitness traits and F values. In this study, we used Drosophila melanogaster to evaluate the effectiveness of molecular markers in estimating ID under suboptimal conditions. We generated two sets of 100 pairs of unrelated individuals from a large panmictic population and mated them for two generations to produce non-inbred and unrelated individuals (F = 0) and inbred individuals (full-sib progeny; F = 0.25). Using these expected genealogical F values, we calculated inbreeding depression for two fitness-related traits, pupae productivity and competitive fitness. We then sequenced the males from 17 non-inbred pairs and 17 inbred pairs to obtain their genomic inbreeding coefficients and estimate ID for the two traits. The scenario assumed was rather restrictive in terms of estimation of ID because: (1) the individuals belonged to the same generation of a large panmictic population, leading to low variation in individual F coefficients; (2) the sample sizes were small; and (3) the traits measured depended on both males and females while only males were sequenced. Despite the challenging conditions of our study, we found that molecular markers provided estimates of ID that were comparable to those obtained from simple pedigree estimations with larger sample sizes. The results therefore suggest that genomic measures of inbreeding are useful to provide estimates of inbreeding depression even under very challenging scenarios.

Kinship composition in mammals

Understanding the evolution of group-living and cooperation requires information on who animals live and cooperate with. Animals can live with kin, non-kin or both, and kinship structure can influence the benefits and costs of group-living and the evolution of within-group cooperation. One aspect of kinship structure is kinship composition, i.e. a group-level attribute of the presence of kin and/or non-kin dyads in groups. Despite its putative importance, the kinship composition of mammalian groups has yet to be characterized. Here, we use the published literature to build an initial kinship composition dataset in mammals, laying the groundwork for future work in the field. In roughly half of the 18 species in our sample, individuals lived solely with same-sex kin, and, in the other half, individuals lived with related and unrelated individuals of the same sex. These initial results suggest that it is not rare for social mammals to live with unrelated individuals of the same sex, highlighting the importance of considering indirect and direct fitness benefits as co-drivers of the evolution of sociality. We hope that our initial dataset and insights will spur the study of kinship structure and sociality towards new exciting avenues.

A simulation-based evaluation of methods for estimating census population size of terrestrial game species from genetically-identified parent-offspring pairs

Estimates of wildlife population size are critical for conservation and management, but accurate estimates are difficult to obtain for many species. Several methods have recently been developed that estimate abundance using kinship relationships observed in genetic samples, particularly parent-offspring pairs. While these methods are similar to traditional Capture-Mark-Recapture, they do not need physical recapture, as individuals are considered recaptured if a sample contains one or more close relatives. This makes methods based on genetically-identified parent-offspring pairs particularly interesting for species for which releasing marked animals back into the population is not desirable or not possible (e.g., harvested fish or game species). However, while these methods have successfully been applied in commercially important fish species, in the absence of life-history data, they are making several assumptions unlikely to be met for harvested terrestrial species. They assume that a sample contains only one generation of parents and one generation of juveniles of the year, while more than two generations can coexist in the hunting bags of long-lived species, or that the sampling probability is the same for each individual, an assumption that is violated when fecundity and/or survival depend on sex or other individual traits. In order to assess the usefulness of kin-based methods to estimate population sizes of terrestrial game species, we simulated population pedigrees of two different species with contrasting demographic strategies (wild boar and red deer), applied four different methods and compared the accuracy and precision of their estimates. We also performed a sensitivity analysis, simulating population pedigrees with varying fecundity characteristics and various levels of harvesting to identify optimal conditions of applicability of each method. We showed that all these methods reached the required levels of accuracy and precision to be effective in wildlife management under simulated circumstances (i.e., for species within a given range of fecundity and for a given range of sampling intensity), while being robust to fecundity variation. Despite the potential usefulness of the methods for terrestrial game species, care is needed as several biases linked to hunting practices still need to be investigated (e.g., when hunting bags are biased toward a particular group of individuals).

Taking quantitative genomics into the wild

We organized this special issue to highlight new work and review recent advances at the cutting edge of 'wild quantitative genomics'. In this editorial, we will present some history of wild quantitative genetic and genomic studies, before discussing the main themes in the papers published in this special issue and highlighting the future outlook of this dynamic field.

Sex ratio and relatedness in the Griffon vulture (Gyps fulvus) population of Serbia

Background

Once a widespread species across the region of Southeast Europe, the Griffon vulture is now confined to small and isolated populations across the Balkan Peninsula. The population from Serbia represents its biggest and most viable population that can serve as an important reservoir of genetic diversity from which the birds can be used for the region's reintroduction programmes. The available genetic data for this valuable population are scarce and as a protected species that belongs to the highly endangered vulture group, it needs to be well described so that it can be properly managed and used as a restocking population. Considering the serious recent bottleneck event that the Griffon vulture population from Serbia experienced we estimated the overall relatedness among the birds from this population. Sex ratio, another important parameter that shows the vitality and strength of the population was evaluated as well.

Methods

During the annual monitoring that was performed in the period from 2013-2021, we collected blood samples from individual birds that were marked in the nests. In total, 169 samples were collected and each was used for molecular sexing while 58 presumably unrelated birds from different nests were used for inbreeding and relatedness analyses. The relatedness was estimated using both biparentally (10 microsatellite loci) and uniparentally (Cytb and D-loop I of mitochondrial DNA) inherited markers.

Results

The level of inbreeding was relatively high and on average it was 8.3% while the mean number of relatives for each bird was close to three. The sex ratio was close to 1:1 and for the analysed period of 9 years, it didn't demonstrate a statistically significant deviation from the expected ratio of 1:1, suggesting that this is a stable and healthy population. Our data suggest that, even though a relatively high level of inbreeding can be detected among the individual birds, the Griffon vulture population from Serbia can be used as a source population for restocking and reintroduction programmes in the region. These data combined with previously observed genetic differentiation between the populations from the Iberian and Balkan Peninsulas suggest that the introduction of foreign birds should be avoided and that local birds should be used instead.

Comparing genome-based estimates of relatedness for use in pedigree-based conservation management

Researchers have long debated which estimator of relatedness best captures the degree of relationship between two individuals. In the genomics era, this debate continues, with relatedness estimates being sensitive to the methods used to generate markers, marker quality, and levels of diversity in sampled individuals. Here, we compare six commonly used genome-based relatedness estimators (kinship genetic distance (KGD), Wang Maximum Likelihood (TrioML), Queller and Goodnight (R_xy ), Kinship INference for Genome-wide association studies (KING-robust), and Pairwise Relatedness (R_AB ), allele-sharing co-ancestry (AS)) across five species bred in captivity-including three birds and two mammals-with varying degrees of reliable pedigree data, using reduced-representation and whole genome resequencing data. Genome-based relatedness estimates varied widely across estimators, sequencing methods, and species, yet the most consistent results for known first order relationships were found using R_xy , R_AB , and AS. However, AS was found to be less consistently correlated with known pedigree relatedness than either R_xy or R_AB . Our combined results indicate there is not a single genome-based estimator that is ideal across different species and data types. To determine the most appropriate genome-based relatedness estimator for each new dataset, we recommend assessing the relative: (1) correlation of candidate estimators with known relationships in the pedigree and (2) precision of candidate estimators with known first-order relationships. These recommendations are broadly applicable to conservation breeding programs, particularly where genome-based estimates of relatedness can complement and complete poorly pedigreed populations. Given a growing interest in the application of wild pedigrees, our results are also applicable to in-situ wildlife management.

Age-specific, population-level pedigree of wild black bears provides insights into reproduction, paternity, and maternal effects on offspring apparent survival

Wildlife pedigrees provide insights into ecological and evolutionary processes. DNA obtained from noninvasively collected hair is often used to determine individual identities for pedigrees and other genetic analyses. However, detection rates associated with some noninvasive DNA studies can be relatively low, and genetic data do not provide information on individual birth year. Supplementing hair DNA stations with video cameras should increase the individual detection rate, assuming accurate identification of individuals via video data. Video data can also provide birth year information for individuals captured as young of the year, which can enrich population-level pedigrees. We placed video cameras at hair stations and combined genetic and video data to reconstruct an age-specific, population-level pedigree of wild black bears during 2010-2020. Combining individual birth year with mother-offspring relatedness, we also estimated litter size, interlitter interval, primiparity, and fecundity. We used the Cormack-Jolly-Seber model in Program Mark to evaluate the effect of maternal identity on offspring apparent survival. We compared model rankings of apparent survival and parameter estimates based on combined genetic and video data with those based on only genetic data. We observed 42 mother-offspring relationships. Of these, 21 (50%) would not have been detected had we used hair DNA alone. Moreover, video data allowed for the cub and yearling age classes to be determined. Mean annual fecundity was 0.42 (95% CI: 0.27, 0.56). Maternal identity influenced offspring apparent survival, where offspring of one mother experienced significantly lower apparent survival (0.39; SE = 0.15) than that of offspring of four other mothers (0.89-1.00; SE = 0.00-0.06). We video-documented cub abandonment by the mother whose offspring experienced low apparent survival, indicating individual behaviors (e.g., maternal care) may scale up to affect population-level parameters (e.g., cub survival). Our findings provide insights into evolutionary processes and are broadly relevant to wildlife ecology and conservation.

Lack of inbreeding avoidance during mate selection in migratory monarch butterflies

Inbreeding is generally thought to have negative consequences for organismal health. However, despite the potential fitness effects, it remains surprisingly common among wild populations. In many cases, the complex factors that underlie mating dynamics make predicting whether individuals should or do avoid inbreeding quite challenging. One reason inbreeding may persist among species is that the likelihood of encountering relatives can be rare. Thus, even if inbreeding has severe consequences, selection to avoid mating with kin will be weak in species that are highly dispersed. Here we investigated if migratory monarch butterflies (Danaus plexippus), which are famous for their dispersal ability, actively avoid inbreeding. We found that neither female nor male monarchs choose mates based on relatedness. These results support the hypothesis that movement ecology can mask the deleterious effects of inbreeding and relax selection for active inbreeding avoidance behaviors. Overall, our data add to the growing list of studies showing that inbreeding avoidance is not the behavioral "default" for most species. We also highlight the implications that inbreeding may have on the declining populations of this iconic butterfly.

How much is enough? Sampling intensity influences estimates of reproductive variance in an introduced population

Conservation introductions to islands and fenced enclosures are increasing as in situ mitigations fail to keep pace with population declines. Few studies consider the potential loss of genetic diversity and increased inbreeding if released individuals breed disproportionately. As funding is limited and post-release monitoring expensive for conservation programs, understanding how sampling effort influences estimates of reproductive variance is useful. To investigate this relationship, we used a well-studied population of Tasmanian devils (Sarcophilus harrisii) introduced to Maria Island, Tasmania, Australia. Pedigree reconstruction based on molecular data revealed high variance in number of offspring per breeder and high proportions of unsuccessful individuals. Computational subsampling of 20%, 40%, 60%, and 80% of observed offspring resulted in inaccurate estimates of reproductive variance compared to the pedigree reconstructed with all sampled individuals. With decreased sampling effort, the proportion of inferred unsuccessful individuals was overestimated and the variance in number of offspring per breeder was underestimated. To accurately estimate reproductive variance, we recommend sampling as many individuals as logistically possible during the early stages of population establishment. Further, we recommend careful selection of colonizing individuals as they may be disproportionately represented in subsequent generations. Within the conservation management context, our results highlight important considerations for sample collection and post-release monitoring during population establishment.

Investigating population dynamics from parentage analysis in the highly endangered fan mussel Pinna nobilis

Understanding dispersal patterns is a major focus for conservation biology as it influences local survival and resilience in case of local disturbance, particularly for sessile species. Dispersal can be assessed through parentage analyses by estimating family structure and self-recruitment. This study documents the family structure of a pelagic spawner, Pinna nobilis, which is facing a major crisis that threatens its survival as most of its populations have been decimated by a parasite, Haplosporidium pinnae. In this context, we focused on a single population (Peyrefite, Banyuls-sur-mer, France) where 640 individuals were sampled in 2011, 2015, and 2018 and genotyped for 22 microsatellite markers. Genetic diversity was high and homogeneous among years, with mean allele numbers ranging between 13.6 and 14.8 and observed heterozygosities (H o) between 0.7121 and 0.7331. Low, but significant, genetic differentiations were found between 2011-2015 and 2015-2018. A parentage analysis described 11 clusters, including one prevailing, and revealed that 46.9% of individuals were involved in half-sib relationships, even between years, suggesting that source populations were recurrent year after year. There were few individuals resampled between years (30 in 2015 and 14 in 2018), indicating a rapid turnover. Considering the large number of half-sib relationships but the low number of relations per individual, we conclude that P. nobilis exhibit homogeneous reproductive success. Self-recruitment was not detected, making this population highly vulnerable as replenishment only relies on connectivity from neighboring populations. In the context of the pandemic caused by H. pinnae, these results will have to be considered when choosing a location to reintroduce individuals in potential future rescue plans.

Contributions of genetic and nongenetic sources to variation in cooperative behavior in a cooperative mammal

The evolution of cooperative behavior is a major area of research among evolutionary biologists and behavioral ecologists, yet there are few estimates of its heritability or its evolutionary potential, and long-term studies of identifiable individuals are required to disentangle genetic and nongenetic components of cooperative behavior. Here, we use long-term data on over 1800 individually recognizable wild meerkats (Suricata suricatta) collected over 30 years and a multigenerational genetic pedigree to partition phenotypic variation in three cooperative behaviors (babysitting, pup feeding, and sentinel behavior) into individual, additive genetic, and other sources, and to assess their repeatability and heritability. In addition to strong effects of sex, age, and dominance status, we found significant repeatability in individual contributions to all three types of cooperative behavior both within and across breeding seasons. Like most other studies of the heritability of social behavior, we found that the heritability of cooperative behavior was low. However, our analysis suggests that a substantial component of the repeatable individual differences in cooperative behavior that we observed was a consequence of additive genetic variation. Our results consequently indicate that cooperative behavior can respond to selection, and suggest scope for further exploration of the genetic basis of social behavior.

Relatedness within and between Agta residential groups

Theoretical models relating to the evolution of human behaviour usually make assumptions about the kinship structure of social groups. Since humans were hunter-gatherers for most of our evolutionary history, data on the composition of contemporary hunter-gatherer groups has long been used to inform these models. Although several papers have taken a broad view of hunter-gatherer social organisation, it is also useful to explore data from single populations in more depth. Here, we describe patterns of relatedness among the Palanan Agta, hunter-gatherers from the northern Philippines. Across 271 adults, mean relatedness to adults across the population is r = 0.01 and to adult campmates is r = 0.074, estimates that are similar to those seen in other hunter-gatherers. We also report the distribution of kin across camps, relatedness and age differences between spouses, and the degree of shared reproductive interest between camp mates, a measure that incorporates affinal kinship. For both this this measure (s) and standard relatedness (r), we see no major age or sex differences in the relatedness of adults to their campmates, conditions that may reduce the potential for conflicts of interest within social groups.

Are immigrants outbred and unrelated? Testing standard assumptions in a wild metapopulation

Immigration into small recipient populations is expected to alleviate inbreeding and increase genetic variation, and hence facilitate population persistence through genetic and/or evolutionary rescue. Such expectations depend on three standard assumptions: that immigrants are outbred, unrelated to existing natives at arrival, and unrelated to each other. These assumptions are rarely explicitly verified, including in key field systems in evolutionary ecology. Yet, they could be violated due to non-random or repeated immigration from adjacent small populations. We combined molecular genetic marker data for 150-160 microsatellite loci with comprehensive pedigree data to test the three assumptions for a song sparrow (Melospiza melodia) population that is a model system for quantifying effects of inbreeding and immigration in the wild. Immigrants were less homozygous than existing natives on average, with mean homozygosity that closely resembled outbred natives. Immigrants can therefore be considered outbred on the focal population scale. Comparisons of homozygosity of real or hypothetical offspring of immigrant-native, native-native and immigrant-immigrant pairings implied that immigrants were typically unrelated to existing natives and to each other. Indeed, immigrants' offspring would be even less homozygous than outbred individuals on the focal population scale. The three standard assumptions of population genetic and evolutionary theory were consequently largely validated. Yet, our analyses revealed some deviations that should be accounted for in future analyses of heterosis and inbreeding depression, implying that the three assumptions should be verified in other systems to probe patterns of non-random or repeated dispersal and facilitate precise and unbiased estimation of key evolutionary parameters.

Genomic data of different resolutions reveal consistent inbreeding estimates but contrasting homozygosity landscapes for the threatened Aotearoa New Zealand hihi

Inbreeding can lead to a loss of heterozygosity in a population and when combined with genetic drift may reduce the adaptive potential of a species. However, there is uncertainty about whether resequencing data can provide accurate and consistent inbreeding estimates. Here, we performed an in-depth inbreeding analysis for hihi (Notiomystis cincta), an endemic and nationally vulnerable passerine bird of Aotearoa New Zealand. We first focused on subsampling variants from a reference genome male, and found that low-density data sets tend to miss runs of homozygosity (ROH) in some places and overestimate ROH length in others, resulting in contrasting homozygosity landscapes. Low-coverage resequencing and 50 K SNP array densities can yield comparable inbreeding results to high-coverage resequencing approaches, but the results for all data sets are highly dependent on the software settings employed. Second, we extended our analysis to 10 hihi where low-coverage whole genome resequencing, RAD-seq and SNP array genotypes are available. We inferred ROH and individual inbreeding to evaluate the relative effects of sequencing depth versus SNP density on estimating inbreeding coefficients and found that high rates of missingness downwardly bias both the number and length of ROH. In summary, when using genomic data to evaluate inbreeding, studies must consider that ROH estimates are heavily dependent on analysis parameters, data set density and individual sequencing depth.

Genomic and fitness consequences of inbreeding in an endangered carnivore

Reduced fitness through genetic drift and inbreeding is a major threat to small and isolated populations. Although previous studies have generally used genetically verified pedigrees to document effects of inbreeding and gene flow, these often fail to capture the whole inbreeding history of the species. By assembling a draft arctic fox (Vulpes lagopus) genome and resequencing complete genomes of 23 additional foxes born before and after a well-documented immigration event in Scandinavia, we here look into the genomic consequences of inbreeding and genetic rescue. We found a difference in genome-wide diversity, with 18% higher heterozygosity and 81% lower F_ROH in immigrant F1 compared to native individuals. However, more distant descendants of immigrants (F2, F3) did not show the same pattern. We also found that foxes with lower inbreeding had higher probability to survive their first year of life. Our results demonstrate the important link between genetic variation and fitness as well as the transient nature of genetic rescue. Moreover, our results have implications in conservation biology as they demonstrate that inbreeding depression can effectively be detected in the wild by a genomic approach.

Spatial familial networks to infer demographic structure of wild populations

In social species, reproductive success and rates of dispersal vary among individuals resulting in spatially structured populations. Network analyses of familial relationships may provide insights on how these parameters influence population-level demographic patterns. These methods, however, have rarely been applied to genetically derived pedigree data from wild populations.Here, we use parent-offspring relationships to construct familial networks from polygamous boreal woodland caribou (Rangifer tarandus caribou) in Saskatchewan, Canada, to inform recovery efforts. We collected samples from 933 individuals at 15 variable microsatellite loci along with caribou-specific primers for sex identification. Using network measures, we assess the contribution of individual caribou to the population with several centrality measures and then determine which measures are best suited to inform on the population demographic structure. We investigate the centrality of individuals from eighteen different local areas, along with the entire population.We found substantial differences in centrality of individuals in different local areas, that in turn contributed differently to the full network, highlighting the importance of analyzing networks at different scales. The full network revealed that boreal caribou in Saskatchewan form a complex, interconnected familial network, as the removal of edges with high betweenness did not result in distinct subgroups. Alpha, betweenness, and eccentricity centrality were the most informative measures to characterize the population demographic structure and for spatially identifying areas of highest fitness levels and family cohesion across the range. We found varied levels of dispersal, fitness, and cohesion in family groups. Synthesis and applications: Our results demonstrate the value of different network measures in assessing genetically derived familial networks. The spatial application of the familial networks identified individuals presenting different fitness levels, short- and long-distance dispersing ability across the range in support of population monitoring and recovery efforts.

Life stage-specific inbreeding depression in long-lived Pinaceae species depends on population connectivity

Inbreeding depression (ID) is a fundamental selective pressure that shapes mating systems and population genetic structures in plants. Although it has been shown that ID varies over the life stages of shorter-lived plants, less is known about how the fitness effects of inbreeding vary across life stages in long-lived species. We conducted a literature survey in the Pinaceae, a tree family known to harbour some of the highest mutational loads ever reported. Using a meta-regression model, we investigated distributions of inbreeding depression over life stages, adjusting for effects of inbreeding levels and the genetic differentiation of populations within species. The final dataset contained 147 estimates of ID across life stages from 41 studies. 44 Fst estimates were collected from 40 peer-reviewed studies for the 18 species to aid genetic differentiation modelling. Partitioning species into fragmented and well-connected groups using Fst resulted in the best way (i.e. trade-off between high goodness-of-fit of the model to the data and reduced model complexity) to incorporate genetic connectivity in the meta-regression analysis. Inclusion of a life stage term and its interaction with the inbreeding coefficient (F) dramatically increased model precision. We observed that the correlation between ID and F was significant at the earliest life stage. Although partitioning of species populations into fragmented and well-connected groups explained little of the between-study heterogeneity, the inclusion of an interaction between life stage and population differentiation revealed that populations with fragmented distributions suffered lower inbreeding depression at early embryonic stages than species with well-connected populations. There was no evidence for increased ID in late life stages in well-connected populations, although ID tended to increase across life stages in the fragmented group. These findings suggest that life stage data should be included in inbreeding depression studies and that inbreeding needs to be managed over life stages in commercial populations of long-lived plants.

Pedigree-based assessment of recent population connectivity in a threatened rattlesnake

Managing endangered species in fragmented landscapes requires estimating dispersal rates between populations over contemporary timescales. Here, we developed a new method for quantifying recent dispersal using genetic pedigree data for close and distant kin. Specifically, we describe an approach that infers missing shared ancestors between pairs of kin in habitat patches across a fragmented landscape. We then applied a stepping-stone model to assign unsampled individuals in the pedigree to probable locations based on minimizing the number of movements required to produce the observed locations in sampled kin pairs. Finally, we used all pairs of reconstructed parent-offspring sets to estimate dispersal rates between habitat patches under a Bayesian model. Our approach measures connectivity over the timescale represented by the small number of generations contained within the pedigree and so is appropriate for estimating the impacts of recent habitat changes due to human activity. We used our method to estimate recent movement between newly discovered populations of threatened Eastern Massasauga rattlesnakes (Sistrurus catenatus) using data from 2996 RAD-based genetic loci. Our pedigree analyses found no evidence for contemporary connectivity between five genetic groups, but, as validation of our approach, showed high dispersal rates between sample sites within a single genetic cluster. We conclude that these five genetic clusters of Eastern Massasauga rattlesnakes have small numbers of resident snakes and are demographically isolated conservation units. More broadly, our methodology can be widely applied to determine contemporary connectivity rates, independent of bias from shared genetic similarity due to ancestry that impacts other approaches.

The thermal environment at fertilization mediates adaptive potential in the sea

Additive genetic variation for fitness at vulnerable life stages governs the adaptive potential of populations facing stressful conditions under climate change, and can depend on current conditions as well as those experienced by past stages or generations. For sexual populations, fertilization is the key stage that links one generation to the next, yet the effects of fertilization environment on the adaptive potential at the vulnerable stages that then unfold during development are rarely considered, despite climatic stress posing risks for gamete function and fertility in many taxa and external fertilizers especially. Here, we develop a simple fitness landscape model exploring the effects of environmental stress at fertilization and development on the adaptive potential in early life. We then test our model with a quantitative genetic breeding design exposing family groups of a marine external fertilizer, the tubeworm Galeolaria caespitosa, to a factorial manipulation of current and projected temperatures at fertilization and development. We find that adaptive potential in early life is substantially reduced, to the point of being no longer detectable, by genotype‐specific carryover effects of fertilization under projected warming. We interpret these results in light of our fitness landscape model, and argue that the thermal environment at fertilization deserves more attention than it currently receives when forecasting the adaptive potential of populations confronting climate change.

Sequenced-based paternity analysis to improve breeding and identify self-incompatibility loci in intermediate wheatgrass (Thinopyrum intermedium)

KEY MESSAGE

Paternity assignment and genome-wide association analyses for fertility were applied to a Thinopyrum intermedium breeding program. A lack of progeny between combinations of parents was associated with loci near self-incompatibility genes. In outcrossing species such as intermediate wheatgrass (IWG, Thinopyrum intermedium), polycrossing is often used to generate novel recombinants through each cycle of selection, but it cannot track pollen-parent pedigrees and it is unknown how self-incompatibility (SI) genes may limit the number of unique crosses obtained. This study investigated the potential of using next-generation sequencing to assign paternity and identify putative SI loci in IWG. Using a reference population of 380 individuals made from controlled crosses of 64 parents, paternity was assigned with 92% agreement using Cervus software. Using this approach, 80% of 4158 progeny (n = 3342) from a polycross of 89 parents were assigned paternity. Of the 89 pollen parents, 82 (92%) were represented with 1633 unique full-sib families representing 42% of all potential crosses. The number of progeny per successful pollen parent ranged from 1 to 123, with number of inflorescences per pollen parent significantly correlated to the number of progeny (r = 0.54, p < 0.001). Shannon's diversity index, assessing the total number and representation of families, was 7.33 compared to a theoretical maximum of 8.98. To test our hypothesis on the impact of SI genes, a genome-wide association study of the number of progeny observed from the 89 parents identified genetic effects related to non-random mating, including marker loci located near putative SI genes. Paternity testing of polycross progeny can impact future breeding gains by being incorporated in breeding programs to optimize polycross methodology, maintain genetic diversity, and reveal genetic architecture of mating patterns.

Leaf chemical defences and insect herbivory in oak: accounting for canopy position unravels marked genetic relatedness effects

BACKGROUND AND AIMS

Highly controlled experiments document that plant genetic diversity and relatedness can shape herbivore communities and patterns of herbivory. Evidence from the field is, however, scarce and inconsistent. We assessed whether a genetic signal underlying herbivory can be detected in oak woodlands when accounting for variation at smaller (within-tree) and larger (among-stand) scales.

METHODS

We tested relationships between tree genetic relatedness, leaf chemical defences and insect herbivory for different canopy layers in 240 trees from 15 pedunculate oak (Quercus robur) forest stands. We partitioned sources of variability in herbivory and defences among stands, individuals and branches.

KEY RESULTS

Leaf defences, insect herbivory and their relationship differed systematically between the upper and the lower tree canopy. When accounting for this canopy effect, the variation explained by tree genetic relatedness rose from 2.8 to 34.1 % for herbivory and from 7.1 to 13.8 % for leaf defences. The effect was driven by markedly stronger relationships in the upper canopy.

CONCLUSIONS

Our findings illustrate that considerable effects of the host plant genotype on levels of leaf chemical defences and associated insect herbivory can be detected in natural tree populations when within-individual variation is properly accounted for.

Power of a dual-use SNP panel for pedigree reconstruction and population assignment

The use of high-throughput, low-density sequencing approaches has dramatically increased in recent years in studies of eco-evolutionary processes in wild populations and domestication in commercial aquaculture. Most of these studies focus on identifying panels of SNP loci for a single downstream application, whereas there have been few studies examining the trade-offs for selecting panels of markers for use in multiple applications. Here, we detail the use of a bioinformatic workflow for the development of a dual-purpose SNP panel for parentage and population assignment, which included identifying putative SNP loci, filtering for the most informative loci for the two tasks, designing effective multiplex PCR primers, optimizing the SNP panel for performance, and performing quality control steps for downstream applications. We applied this workflow to two adjacent Alaskan Sockeye Salmon populations and identified a GTseq panel of 142 SNP loci for parentage and 35 SNP loci for population assignment. Only 50-75 panel loci were necessary for >95% accurate parentage, whereas population assignment success, with all 172 panel loci, ranged from 93.9% to 96.2%. Finally, we discuss the trade-offs and complexities of the decision-making process that drives SNP panel development, optimization, and testing.

Are shed hair genomes the most effective noninvasive resource for estimating relationships in the wild?

Knowledge of relationships in wild populations is critical for better understanding mating systems and inbreeding scenarios to inform conservation strategies for endangered species. To delineate pedigrees in wild populations, study genetic connectivity, study genotype-phenotype associations, trace individuals, or track wildlife trade, many identified individuals need to be genotyped at thousands of loci, mostly from noninvasive samples. This requires us to (a) identify the most common noninvasive sample available from identified individuals, (b) assess the ability to acquire genome-wide data from such samples, and (c) evaluate the quality of such genome-wide data, and its ability to reconstruct relationships between animals within a population.We followed identified individuals from a wild endangered tiger population and found that shed hair samples were the most common compared to scat samples, opportunistically found carcasses, and opportunistic invasive samples. We extracted DNA from these samples, prepared whole genome sequencing libraries, and sequenced genomes from these.Whole genome sequencing methods resulted in between 25%-98% of the genome sequenced for five such samples. Exploratory population genetic analyses revealed that these data were free of holistic biases and could recover expected population structure and relatedness. Mitochondrial genomes recovered matrilineages in accordance with long-term monitoring data. Even with just five samples, we were able to uncover the matrilineage for three individuals with unknown ancestry.In summary, we demonstrated that noninvasive shed hair samples yield adequate quality and quantity of DNA in conjunction with sensitive library preparation methods, and provide reliable data from hundreds of thousands of SNPs across the genome. This makes shed hair an ideal noninvasive resource for studying individual-based genetics of elusive endangered species in the wild.

Comparison of divergent breeding management strategies in two species of semi-captive eland in Senegal

Breeding management of small populations may have a critical influence on the development of population characteristics in terms of levels of genetic diversity and inbreeding. Two populations of antelope sister species - Critically Endangered Western Derby eland (Tauroragus derbianus derbianus) and non-native Least Concern Cape eland (Taurotragus oryx oryx) bred under different management strategies were studied in Senegal, Western Africa. The aims of the study were to compare the population genetic parameters of the two species and to test for the presence of interspecific hybrids. In total, blood and tissue samples from 76 Western Derby elands and 26 Cape elands were investigated, using 12 microsatellite markers. No hybrid individuals were detected in the sampled animals within the multispecies enclosure in Bandia Reserve, Senegal. The parameters of genetic polymorphism indicated much lower genetic diversity in Western Derby elands compared to Cape elands. On the other hand, the coefficient of inbreeding was low in both species. It is hypothesized that this could be a positive effect of strict population management of Western Derby elands, which, despite the loss of genetic diversity, minimizes inbreeding.

A comparison of pedigree, genetic and genomic estimates of relatedness for informing pairing decisions in two critically endangered birds: Implications for conservation breeding programmes worldwide

Conservation management strategies for many highly threatened species include conservation breeding to prevent extinction and enhance recovery. Pairing decisions for these conservation breeding programmes can be informed by pedigree data to minimize relatedness between individuals in an effort to avoid inbreeding, maximize diversity and maintain evolutionary potential. However, conservation breeding programmes struggle to use this approach when pedigrees are shallow or incomplete. While genetic data (i.e., microsatellites) can be used to estimate relatedness to inform pairing decisions, emerging evidence indicates this approach may lack precision in genetically depauperate species, and more effective estimates will likely be obtained from genomic data (i.e., thousands of genome-wide single nucleotide polymorphisms, or SNPs). Here, we compare relatedness estimates and subsequent pairing decisions using pedigrees, microsatellites and SNPs from whole-genome resequencing approaches in two critically endangered birds endemic to New Zealand: kakī/black stilt (Himantopus novaezelandiae) and kākāriki karaka/orange-fronted parakeet (Cyanoramphus malherbi). Our findings indicate that SNPs provide more precise estimates of relatedness than microsatellites when assessing empirical parent-offspring and full sibling relationships. Further, our results show that relatedness estimates and subsequent pairing recommendations using PMx are most similar between pedigree- and SNP-based approaches. These combined results indicate that in lieu of robust pedigrees, SNPs are an effective tool for informing pairing decisions, which has important implications for many poorly pedigreed conservation breeding programmes worldwide.

Conservation: The Costs of Inbreeding and of Being Inbred

Endangered species face a huge array of challenges, including the negative consequences of individuals having to breed with close genetic relatives. But just how costly is inbreeding in small populations? New research from an endangered bird species suggests that considering inbreeding could be crucial for conservation programmes.

In situ estimation of genetic variation of functional and ecological traits in Quercus petraea and Q.robur

BACKGROUND

Predicting the evolutionary potential of natural tree populations requires the estimation of heritability and genetic correlations among traits on which selection acts, as differences in evolutionary success between species may rely on differences for these genetic parameters. In situ estimates are expected to be more accurate than measures done under controlled conditions which do not reflect the natural environmental variance.

AIMS

The aim of the current study was to estimate three genetic parameters (i.e. heritability, evolvability and genetic correlations) in a natural mixed oak stand composed of Quercus petraea and Quercus robur about 100 years old, for 58 traits of ecological and functional relevance (growth, reproduction, phenology, physiology, resilience, structure, morphology and defence).

METHODS

First we estimated genetic parameters directly in situ using realized genomic relatedness of adult trees and parentage relationships over two generations to estimate the traits additive variance. Secondly, we benefited from existing ex situ experiments (progeny tests and conservation collection) installed with the same populations, thus allowing comparisons of in situ heritability estimates with more traditional methods.

RESULTS

Heritability and evolvability estimates obtained with different methods varied substantially and showed large confidence intervals, however we found that in situ were less precise than ex situ estimates, and assessments over two generations (with deeper relatedness) improved estimates of heritability while large sampling sizes are needed for accurate estimations. At the biological level, heritability values varied moderately across different ecological and functional categories of traits, and genetic correlations among traits were conserved over the two species.

CONCLUSION

We identified limits for using realized genomic relatedness in natural stands to estimate the genetic variance, given the overall low variance of genetic relatedness and the rather low sampling sizes of currently used long term genetic plots in forestry. These limits can be overcome if larger sample sizes are considered, or if the approach is extended over the next generation.

Genetic Parameters for Growth and Kid Survival of Indigenous Goat under Smallholding System of Burundi

The goal of this study was to estimate genetic parameters for the growth, conformation, and survival of goat kids raised in smallholder farming systems in Burundi. To do this, measurements were taken on live weight, thoracic perimeter, length, and height at birth (n = 1538 animals), at 3 months (n = 1270 animals), at 6 months (n = 992 animals), at 9 months (n = 787 animals), and at 12 months (n = 705 animals). Kids were born between 2016 and 2019, from 645 dams and 106 bucks. Three bivariate animal models were used to estimate genetic parameters of body weight and conformation measurements as potential indicators of this weight. According to the measure, heritability was estimated between 15 and 17% and genetic correlations between 65 and 79%. An accelerated failure time animal model was used to estimate the heritability of survival for kids under one year, adjusted for birth weight. Goat survival was significantly prolonged by 0.64 days per kilogram of birth weight. The estimated heritability for this trait was 2%. Overall, these results suggest that a selection program could be implemented to improve animal growth, either directly on weight or indirectly on conformational traits. At the same time, efforts need to be made to improve rearing conditions to increase the survival of kids.

Heterozygosity-Fitness Correlations Reveal Inbreeding Depression in Neonatal Body Size in a Critically Endangered Rock Iguana

Inbreeding depression, though challenging to identify in nature, may play an important role in regulating the dynamics of small and isolated populations. Conversely, greater expression of genetic load can enhance opportunities for natural selection. Conditional expression concentrates these opportunities for selection and may lead to failure of detection. This study investigates the possibility for age-dependent expression of inbreeding depression in a critically endangered population of rock iguanas, Cyclura nubila caymanensis. We employ heterozygote-fitness correlations to examine the contributions of individual genetic factors to body size, a fitness-related trait. Nonsignificant reductions in homozygosity (up to 7%) were detected between neonates and individuals surviving past their first year, which may reflect natural absorption of inbreeding effects by this small, fecund population. The majority of variation in neonate body size was attributed to maternal or environmental effects (i.e., clutch identity and incubation length); however, heterozygosity across 22 microsatellite loci also contributed significantly and positively to model predictions. Conversely, effects of heterozygosity on fitness were not detectable when adults were examined, suggesting that inbreeding depression in body size may be age dependent in this taxon. Overall, these findings emphasize the importance of taking holistic, cross-generational approaches to genetic monitoring of endangered populations.

The role of selection and evolution in changing parturition date in a red deer population

Changing environmental conditions cause changes in the distributions of phenotypic traits in natural populations. However, determining the mechanisms responsible for these changes—and, in particular, the relative contributions of phenotypic plasticity versus evolutionary responses—is difficult. To our knowledge, no study has yet reported evidence that evolutionary change underlies the most widely reported phenotypic response to climate change: the advancement of breeding times. In a wild population of red deer, average parturition date has advanced by nearly 2 weeks in 4 decades. Here, we quantify the contribution of plastic, demographic, and genetic components to this change. In particular, we quantify the role of direct phenotypic plasticity in response to increasing temperatures and the role of changes in the population structure. Importantly, we show that adaptive evolution likely played a role in the shift towards earlier parturition dates. The observed rate of evolution was consistent with a response to selection and was less likely to be due to genetic drift. Our study provides a rare example of observed rates of genetic change being consistent with theoretical predictions, although the consistency would not have been detected with a solely phenotypic analysis. It also provides, to our knowledge, the first evidence of both evolution and phenotypic plasticity contributing to advances in phenology in a changing climate.

Adaptive genetic evolution and phenotypic plasticity both contribute to a two-week advancement of birth dates earlier in spring in a deer population subject to temperature warming over four decades.

Genome-Wide Analysis of SNPs Is Consistent with No Domestic Dog Ancestry in the Endangered Mexican Wolf (Canis lupus baileyi)

The Mexican gray wolf (Canis lupus baileyi) was historically distributed throughout the southwestern United States and northern Mexico. Extensive predator removal campaigns during the early 20th century, however, resulted in its eventual extirpation by the mid 1980s. At this time, the Mexican wolf existed only in 3 separate captive lineages (McBride, Ghost Ranch, and Aragón) descended from 3, 2, and 2 founders, respectively. These lineages were merged in 1995 to increase the available genetic variation, and Mexican wolves were reintroduced into Arizona and New Mexico in 1998. Despite the ongoing management of the Mexican wolf population, it has been suggested that a proportion of the Mexican wolf ancestry may be recently derived from hybridization with domestic dogs. In this study, we genotyped 87 Mexican wolves, including individuals from all 3 captive lineages and cross-lineage wolves, for more than 172000 single nucleotide polymorphisms. We identified levels of genetic variation consistent with the pedigree record and effects of genetic rescue. To identify the potential to detect hybridization with domestic dogs, we compared our Mexican wolf genotypes with those from studies of domestic dogs and other gray wolves. The proportion of Mexican wolf ancestry assigned to domestic dogs was only between 0.06% (SD 0.23%) and 7.8% (SD 1.0%) for global and local ancestry estimates, respectively; and was consistent with simulated levels of incomplete lineage sorting. Overall, our results suggested that Mexican wolves lack biologically significant ancestry with dogs and have useful implications for the conservation and management of this endangered wolf subspecies.