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

Genomic Prediction of Individual Inbreeding Levels for the Management of Genetic Diversity in Populations With Small Effective Size

In populations of small effective size (Ne), such as those in conservation programmes, companion animals or livestock species, inbreeding control is essential. Homozygosity-by-descent (HBD) segments provide relevant information in that context, as they allow accurate estimation of the inbreeding coefficient, provide locus-specific information and their length is informative about the "age" of inbreeding. Our objective was to evaluate tools for predicting HBD in future offspring based on parental genotypes, a problem equivalent to identifying segments identical-by-descent (IBD) among the four parental chromosomes. In total, we reviewed and evaluated 16 approaches using simulated and real data from populations with small Ne. The methods included model-based approaches as well as more computationally efficient rule-based approaches. The accuracy of the methods was then evaluated, including with low-density marker panels, genotyping-by-sequencing data and small groups of individuals, typical features of such populations. Two model-based approaches performed consistently well, while some rule-based approaches proved accurate for genome-wide predictions. The model-based approaches were particularly efficient when genomic information was sparse or degraded. Methods using phased data proved to be more accurate, while some approaches relying on unphased genotype data were sensitive to the assumed allele frequencies. In some settings, pedigree-based predictions ranked high for recent inbreeding levels. Finally, we showed that our evaluation is also informative about the accuracy of the methods for estimating relatedness and identifying IBD segments between pairs of present-day individuals. This study shows that future inbreeding can be accurately predicted, including at specific loci, but not all methods perform equally well.

Potential Adaptive Introgression From Dogs in Iberian Grey Wolves (Canis lupus)

Invading species along with increased anthropogenization may lead to hybridization events between wild species and closely related domesticates. As a consequence, wild species may carry introgressed alleles from domestic species, which is generally assumed to yield adverse effects in wild populations. The opposite evolutionary consequence, adaptive introgression, where introgressed genes are positively selected in the wild species, is possible but has rarely been documented. Grey wolves (Canis lupus) are widely distributed across the Holarctic and frequently coexist with their close relative, the domestic dog (C. familiaris). Despite ample opportunity, hybridization rarely occurs in most populations. Here we studied the geographically isolated grey wolves of the Iberian Peninsula, who have coexisted with a large population of loosely controlled dogs for thousands of years in a human-modified landscape. We assessed the extent and impact of dog introgression on the current Iberian grey wolf population by analysing 150 whole genomes of Iberian and other Eurasian grey wolves as well as dogs originating from across Europe and western Siberia. We identified almost no recent introgression and a small (< 5%) overall ancient dog ancestry. Using a combination of single scan statistics and ancestry enrichment estimates, we identified positive selection on six genes (DAPP1, NSMCE4A, MPPED2, PCDH9, MBTPS1, and CDH13) for which wild Iberian wolves carry alleles introgressed from dogs. The genes with introgressed and positively selected alleles include functions in immune response and brain functions, which may explain some of the unique behavioural phenotypes in Iberian wolves such as their reduced dispersal compared to other wolf populations.

Assessment of genetic diversity, population structure and wolf-dog hybridisation in the Eastern Romanian Carpathian wolf population

The Carpathian Mountains have been constantly inhabited by grey wolves and present one of the largest distribution areas in Europe, comprising between 2300 and 2700 individuals in Romania. To date, however, relatively little is known about the Romanian wolf population. We aimed to provide a first assessment of genetic diversity, population structure and wolf-dog hybridisation based on 444 mostly non-invasively collected samples in the Eastern Romanian Carpathians. Pack reconstruction and analysis of population genetic parameters were performed with mitochondrial DNA control-region sequencing and microsatellite genotyping. We found relatively high levels of genetic diversity, which is similar to values found in previous studies on Carpathian wolves from Poland and Slovakia, as well as to the long-lasting Dinaric-Balkan wolf population. We found no significant population structure in our study region, suggesting effective dispersal and admixture. Analysis of wolf-dog hybridisation using a Single Nucleotide Polymorphism panel optimised for hybrid detection revealed low rates of admixture between wolves and domestic dogs. Our results provide evidence for the existence of a genetically viable wolf population in the Romanian Carpathians. The genetic data obtained in this study may serve as valuable baseline information for the elaboration of monitoring standards and management plans for wolves in Romania.

Identification of Genomic Signatures in Bullmastiff Dogs Using Composite Selection Signals Analysis of 23 Purebred Clades

Dog breeds represent canine sub-populations with distinctive phenotypic features and limited genetic diversity. We have established a resource to study breed-specific genetic diversity. Utilising genetic resources within our laboratory biobank, public domain genotype data and the phylogenetic framework of 23 breed clades, the primary objective for this study was to identify genomic regions that differentiate the Bullmastiff breed. Through application of a composite index analysis (CSS), genomic signatures were identified in Bullmastiffs when compared to the formative breeds, Mastiffs and Bulldogs, and to 22 other breed groups. Significant regions were identified on 15 chromosomes, with the most differentiated regions found on CFA1, CFA9, and CFA18. These regions may reflect genetic drift following establishment of the breed or the effects of selective breeding during development of the modern Bullmastiff. This was supported by analysis of genes from the identified genomic regions, including 458 genes from the multi-clade analysis, which revealed enriched pathways that may be related to characteristic traits and distinct morphology of the breed. The study demonstrates the utility of the CSS method in breed-specific genome analysis and advances our understanding of genetic diversity in Bullmastiff dogs.

Analysis of dog breed diversity using a composite selection index

During breed development, domestic dogs have undergone genetic bottlenecks and sustained selective pressures, as a result distinctive genomic diversity occurs to varying degrees within and between breed groups. This diversity can be identified using standard methods or combinations of these methods. This study explored the application of a combined selection index, composite selection signals (CSS), derived from multiple methods to an existing genotype dataset from three breed groups developed in distinct regions of Asia: Qinghai-Tibet plateau dogs (adapted to living at altitude), Xi dogs (with superior running ability) and Mountain hounds (used for hunting ability). The CSS analysis confirmed top ranked genomic regions on CFA10 and CFA21 in Qinghai-Tibet plateau dogs, CFA1 in Xi dogs and CFA5 in Mountain hounds. CSS analysis identified additional significant genomic regions in each group, defined by a total of 1,397, 1,475 and 1,675 significant SNPs in the Qinghai-Tibetan Plateau dogs, Xi dogs and Mountain hounds, respectively. Chitinase 3 Like 1 (CHI3L1) and Leucine Rich Repeat Containing G Protein-Coupled Receptor 6 (LGR6) genes were located in the top ranked region on CFA7 (0.02-1 Mb) in the Qinghai-Tibetan Plateau dogs. Both genes have been associated with hypoxia responses or altitude adaptation in humans. For the Xi dogs, the top ranked region on CFA25 contained the Transient Receptor Potential Cation Channel Subfamily C Member 4 (TRPC4) gene. This calcium channel is important for optimal muscle performance during exercise. The outstanding signals in the Mountain dogs were on CFA5 with 213 significant SNPs that spanned genes involved in cardiac development, sight and generation of biochemical energy. These findings support the use of the combined index approach for identifying novel regions of genome diversity in dogs. As with other methods, the results do not prove causal links between these regions and phenotypes, but they may assist in focusing future studies that seek to identify functional pathways that contribute to breed diversity.

Long-term Small Population Size, Deleterious Variation, and Altitude Adaptation in the Ethiopian Wolf, a Severely Endangered Canid

Ethiopian wolves, a canid species endemic to the Ethiopian Highlands, have been steadily declining in numbers for decades. Currently, out of 35 extant species, it is now one of the world's most endangered canids. Most conservation efforts have focused on preventing disease, monitoring movements and behavior, and assessing the geographic ranges of sub-populations. Here, we add an essential layer by determining the Ethiopian wolf's demographic and evolutionary history using high-coverage (∼40×) whole-genome sequencing from 10 Ethiopian wolves from the Bale Mountains. We observe exceptionally low diversity and enrichment of weakly deleterious variants in the Ethiopian wolves in comparison with two North American gray wolf populations and four dog breeds. These patterns are consequences of long-term small population size, rather than recent inbreeding. We infer the demographic history of the Ethiopian wolf and find it to be concordant with historic records and previous genetic analyses, suggesting Ethiopian wolves experienced a series of both ancient and recent bottlenecks, resulting in a census population size of fewer than 500 individuals and an estimated effective population size of approximately 100 individuals. Additionally, long-term small population size may have limited the accumulation of strongly deleterious recessive mutations. Finally, as the Ethiopian wolves have inhabited high-altitude areas for thousands of years, we searched for evidence of high-altitude adaptation, finding evidence of positive selection at a transcription factor in a hypoxia-response pathway [CREB-binding protein (CREBBP)]. Our findings are pertinent to continuing conservation efforts and understanding how demography influences the persistence of deleterious variation in small populations.

Human-modified canids in human-modified landscapes: The evolutionary consequences of hybridization for grey wolves and free-ranging domestic dogs

Introgressive hybridization between domestic animals and their wild relatives is an indirect form of human-induced evolution, altering gene pools and phenotypic traits of wild and domestic populations. Although this process is well documented in many taxa, its evolutionary consequences are poorly understood. In this study, we assess introgression patterns in admixed populations of Eurasian wolves and free-ranging domestic dogs (FRDs), identifying chromosomal regions with significantly overrepresented hybrid ancestry and assessing whether genes located within these regions show signatures of selection. Although the dog admixture proportion in West Eurasian wolves (2.7%) was greater than the wolf admixture proportion in FRDs (0.75%), the number and average length of chromosomal blocks showing significant overrepresentation of hybrid ancestry were smaller in wolves than FRDs. In wolves, 6% of genes located within these blocks showed signatures of positive selection compared to 23% in FRDs. We found that introgression from wolves may provide a considerable adaptive advantage to FRDs, counterbalancing some of the negative effects of domestication, which can include reduced genetic diversity and excessive tameness. In wolves, introgression from FRDs is mostly driven by drift, with a small number of positively selected genes associated with brain function and behaviour. The predominance of drift may be the consequence of small effective size of wolf populations, which reduces efficiency of selection for weakly advantageous or against weakly disadvantageous introgressed variants. Small wolf population sizes result largely from human-induced habitat loss and hunting, thus linking introgression rates to anthropogenic processes. Our results imply that maintenance of large population sizes should be an important element of wolf management strategies aimed at reducing introgression rates of dog-derived variants.

Unveiling the Fecal Microbiota in Two Captive Mexican Wolf (Canis lupus baileyi) Populations Receiving Different Type of Diets

Simple Summary

The Mexican wolf (Canis lupus baileyi) is an endangered canine. Both Mexico and the United States are currently collaborating to reproduce and reintroduce individuals to their original habitats. However, keeping these wolves in captivity represents a great commitment to meet their basic needs. Diet is a determining factor that is closely related to health and reproductive fitness. The type of diet that is fed to canines in captivity must provide the required nutrients for their development and welfare. The study of the fecal microbiota is a non-invasive way to establish the abundance and diversity of bacterial communities to determine if they are in a healthy condition. We analyzed data from two captive populations of Mexican wolves (i.e., northern and central Mexico) receiving different type of diets (Michilia population: mainly kibble vs. Ocotal population: mainly raw meat). The operational taxonomic units (OTUs) in Michilia resulted in 204 genera and 316 species, while in Ocotal there were 232 genera and 379 species. In the Michilia, dominance of bacteria that degrade carbohydrates was observed (related to kibble diet). In contrast, the Ocotal microbiota was dominated by protein-degrading bacteria (related to raw meat diet). The main outcomes generated in this study should help to enhance the welfare of the captive Mexican wolves to increase its numbers.

Abstract

The Mexican wolf (Canis lupus baileyi) was once distributed in southern United States and northern Mexico. It is an endangered subspecies detached from the gray wolf, and likely exemplifies one of the original migration waves of C. lupus into the new world. This is a canine whose individuals survive in specialized facilities, zoos, and museums as part of captive-breeding programs. In order to contribute to the improvement of the management of this species and favor its long-term conservation in Mexico, we aimed to evaluate the diversity and abundance of the fecal bacterial microbiota in two populations exposed to different types of diet: (1) Michilia (23° N, 104° W); kibble daily and raw meat sporadically, and (2) Ocotal (19° N, 99° W); raw meat daily and live animals periodically. Next generation sequencing (V3-V4 16S rRNA gene) by Illumina was implemented. The operational taxonomic units (OTUs) in Michilia resulted in 9 phyla, 19 classes, 34 orders, 61 families, 204 genera, and 316 species, while in Ocotal there were 12 phyla, 24 classes, 37 orders, 69 families, 232 genera, and 379 species. Higher estimated Chao1 richness, Shannon diversity, and core microbiota were observed in Ocotal. Differences (p < 0.05) between populations occurred according to the Bray–Curtis beta diversity index. In the Michilia, dominance of bacteria that degrade carbohydrates (Firmicutes, Lachnospiraceae, Blautia, Clostrodium, Eisenbergiella, Romboutsia, and Ruminococcus) was observed; they are abundant in kibble diets. In contrast, the Ocotal microbiota was dominated by protein-degrading bacteria (Fusobacteria, Fusobacteriaceae, and Fusobacteria), indicating a possible positive relation with a raw meat diet. The information generated in this study is fundamental to support the implementation of better management plans in the two populations considered here, as well as in different facilities of southern United States and Mexico, where this subspecies is kept in captivity for conservation purposes.

A sliver of the past: The decimation of the genetic diversity of the Mexican wolf

The endangered Mexican wolf (Canis lupus baileyi) is known to carry exceedingly low levels of genetic diversity. This could be (i) the result of long-term evolutionary patterns as they exist at the southernmost limit of the species distribution at a relatively reduced effective size, or (ii) due to rapid population decline caused by human persecution over the last century. If the former, purifying selection is expected to have minimized the impact of inbreeding. If the latter, rapid and recent declines in genetic diversity may have resulted in severe fitness consequences. To differentiate these hypotheses, we conducted comparative whole-genome analyses of five historical Mexican wolves (1907-1917) and 18 contemporary Mexican and grey wolves from North America and Eurasia. Based on whole-genome data, historical and modern Mexican wolves together form a discrete unit. Moreover, we found that modern Mexican wolves have reduced genetic diversity and increased inbreeding relative to the historical population, which was widespread across the southwestern United States and not restricted to Mexico as previously assumed. Finally, although Mexican wolves have evolved in sympatry with coyotes (C. latrans), we observed lower introgression between historical Mexican wolves and coyotes than with modern Mexican wolves, despite similarities in body size. Taken together, our data show that recent population declines probably caused the reduced level of genetic diversity, but not the observed differentiation of the Mexican wolves from other North American wolves.

Systematic Review on Local Ancestor Inference From a Mathematical and Algorithmic Perspective

Genotypic data provide deep insights into the population history and medical genetics. The local ancestry inference (LAI) (also termed local ancestry deconvolution) method uses the hidden Markov model (HMM) to solve the mathematical problem of ancestry reconstruction based on genomic data. HMM is combined with other statistical models and machine learning techniques for particular genetic tasks in a series of computer tools. In this article, we surveyed the mathematical structure, application characteristics, historical development, and benchmark analysis of the LAI method in detail, which will help researchers better understand and further develop LAI methods. Firstly, we extensively explore the mathematical structure of each model and its characteristic applications. Next, we use bibliometrics to show detailed model application fields and list articles to elaborate on the historical development. LAI publications had experienced a peak period during 2006-2016 and had kept on moving in the following years. The efficiency, accuracy, and stability of the existing models were evaluated by the benchmark. We find that phased data had higher accuracy in comparison with unphased data. We summarize these models with their distinct advantages and disadvantages. The Loter model uses dynamic programming to obtain a globally optimal solution with its parameter-free advantage. Aligned bases can be used directly in the Seqmix model if the genotype is hard to call. This research may help model developers to realize current challenges, develop more advanced models, and enable scholars to select appropriate models according to given populations and datasets.

The impact of identity by descent on fitness and disease in dogs

Dogs and humans have coexisted together for thousands of years, but it was not until the Victorian Era that humans practiced selective breeding to produce the modern standards we see today. Strong artificial selection during the breed formation period has simplified the genetic architecture of complex traits and caused an enrichment of identity-by-descent (IBD) segments in the dog genome. This study demonstrates the value of IBD segments and utilizes them to infer the recent demography of canids, predict case-control status for complex traits, locate regions of the genome potentially linked to inbreeding depression, and to identify understudied breeds where there is potential to discover new disease-associated variants.

Domestic dogs have experienced population bottlenecks, recent inbreeding, and strong artificial selection. These processes have simplified the genetic architecture of complex traits, allowed deleterious variation to persist, and increased both identity-by-descent (IBD) segments and runs of homozygosity (ROH). As such, dogs provide an excellent model for examining how these evolutionary processes influence disease. We assembled a dataset containing 4,414 breed dogs, 327 village dogs, and 380 wolves genotyped at 117,288 markers and data for clinical and morphological phenotypes. Breed dogs have an enrichment of IBD and ROH, relative to both village dogs and wolves, and we use these patterns to show that breed dogs have experienced differing severities of bottlenecks in their recent past. We then found that ROH burden is associated with phenotypes in breed dogs, such as lymphoma. We next test the prediction that breeds with greater ROH have more disease alleles reported in the Online Mendelian Inheritance in Animals (OMIA). Surprisingly, the number of causal variants identified correlates with the popularity of that breed rather than the ROH or IBD burden, suggesting an ascertainment bias in OMIA. Lastly, we use the distribution of ROH across the genome to identify genes with depletions of ROH as potential hotspots for inbreeding depression and find multiple exons where ROH are never observed. Our results suggest that inbreeding has played a large role in shaping genetic and phenotypic variation in dogs and that future work on understudied breeds may reveal new disease-causing variation.

Global Phylogeographic and Admixture Patterns in Grey Wolves and Genetic Legacy of An Ancient Siberian Lineage

The evolutionary relationships between extinct and extant lineages provide important insight into species' response to environmental change. The grey wolf is among the few Holarctic large carnivores that survived the Late Pleistocene megafaunal extinctions, responding to that period's profound environmental changes with loss of distinct lineages and phylogeographic shifts, and undergoing domestication. We reconstructed global genome-wide phylogeographic patterns in modern wolves, including previously underrepresented Siberian wolves, and assessed their evolutionary relationships with a previously genotyped wolf from Taimyr, Siberia, dated at 35 Kya. The inferred phylogeographic structure was affected by admixture with dogs, coyotes and golden jackals, stressing the importance of accounting for this process in phylogeographic studies. The Taimyr lineage was distinct from modern Siberian wolves and constituted a sister lineage of modern Eurasian wolves and domestic dogs, with an ambiguous position relative to North American wolves. We detected gene flow from the Taimyr lineage to Arctic dog breeds, but population clustering methods indicated closer similarity of the Taimyr wolf to modern wolves than dogs, implying complex post-divergence relationships among these lineages. Our study shows that introgression from ecologically diverse con-specific and con-generic populations was common in wolves' evolutionary history, and could have facilitated their adaptation to environmental change.

Biological and Sociopolitical Sources of Uncertainty in Population Viability Analysis for Endangered Species Recovery Planning

Although population viability analysis (PVA) can be an important tool for strengthening endangered species recovery efforts, the extent to which such analyses remain embedded in the social process of recovery planning is often unrecognized. We analyzed two recovery plans for the Mexican wolf that were developed using similar data and methods but arrived at contrasting conclusions as to appropriate recovery goals or criteria. We found that approximately half of the contrast arose from uncertainty regarding biological data, with the remainder divided between policy-related decisions and mixed biological-policy factors. Contrasts arose from both differences in input parameter values and how parameter uncertainty informed the level of precaution embodied in resulting criteria. Policy-related uncertainty originated from contrasts in thresholds for acceptable risk and disagreement as to how to define endangered species recovery. Rather than turning to PVA to produce politically acceptable definitions of recovery that appear science-based, agencies should clarify the nexus between science and policy elements in their decision processes. The limitations we identify in endangered-species policy and how PVAs are conducted as part of recovery planning must be addressed if PVAs are to fulfill their potential to increase the odds of successful conservation outcomes.