Coyotes can do 'puppy dog eyes' too: comparing interspecific variation in Canis facial expression muscles.

Facial expressions are critical for non-verbal communication. The Canis genus epitomizes the interplay between behaviour and morphology in the evolution of non-verbal communication. Recent work suggests that the levator anguli oculi medialis (LAOM) muscle is unique to dogs (Canis familiaris) within the Canis genus and evolved due to domestication. The LAOM raises the inner eyebrows, resulting in the 'puppy dog eyes' expression. Here, we test whether the LAOM is a derived trait in dogs by (i) examining the facial expression muscles of a closely related and ancestral wild Canis species, the coyote (C. latrans) and (ii) comparing our results with other Canis and canid taxa. We discover that coyotes have a well-developed LAOM like dogs, which differs from the modified/absent LAOM in grey wolves. Our findings challenge the hypothesis that the LAOM developed due to domestication. We suggest that the LAOM is a basal trait that was lost in grey wolves. Additionally, we find inter- and intraspecific variations in the size of the muscles of the outer ear, forehead, lips and rostrum, indicating potential adaptations related to sensory perception, communication and individual-level functional variations within canids. Together, this research expands our knowledge of facial expressions, their evolution and their role in communication.

Canine hyper-sociability structural variants associated with altered three-dimensional chromatin state.

Strong selection on complex traits can lead to skewed trait means and reduced trait variability in populations. An example of this phenomenon can be evidenced in allele frequency changes and skewed trait distributions driven by persistent human-directed selective pressures in domesticated species. Dog domestication is linked to several genomic variants; however, the functional impacts of these variants may not always be straightforward when found in non-coding regions of the genome. Four polymorphic transposable elements (TE) found within non-coding sites along a 5 Mb region on canine CFA6 have evolved due to directional selection associated with heightened human-directed hyper-sociability in domesticated dogs. We found that the polymorphic TE in intron 17 of the canine GTF2I gene, which was previously reported to be negatively correlated with canid human-directed hyper-sociability, is associated with altered chromatin looping and hence distinct cis-regulatory landscapes. We reported supporting evidence of an E2F1-DNA binding peak concordant with the altered loop and higher expression of GTF2I exon 18, indicative of alternative splicing. Globally, we discovered differences in pathways regulating the extra-cellular matrix with respect to TE copy number. Overall, we reported evidence suggesting an intriguing molecular convergence between the emergence of hypersocial behaviors in dogs and the same genes that, when hemizygous, produce human Williams Beuren Syndrome characterized by cranio-facial defects and heightened social behaviors. Our results additionally emphasize the often-overlooked potential role of chromatin architecture in social evolution.

Genetic Tracking of a Rabid Coyote (Canis latrans) Detected beyond a Rabies Enzootic Area in West Virginia, USA.

Wildlife translocation and cross-species transmission can impede control and elimination of emerging zoonotic diseases. Tracking the geographic origin of both host and virus (i.e., translocation versus local infection) may help determine the most effective response when high-risk cases of emerging pathogens are identified in wildlife. In May 2022, a coyote (Canis latrans) infected with the raccoon (Procyon lotor) rabies virus variant (RRV) was collected in Lewis County, West Virginia, USA, an area free from RRV. We applied host population genomics and RRV phylogenetic analyses to determine the most likely geographic origin of the rabid coyote. Coyote genomic analyses included animals from multiple eastern states bordering West Virginia, with the probable origin of the rabid coyote being the county of collection. The RRV phylogenetic analyses included cases detected from West Virginia and neighboring states, with most similar RRV sequences collected in a county 80 km to the northeast, within the oral rabies vaccination zone. The combined results suggest that the coyote was infected in an RRV management area and carried the RRV to Lewis County, a pattern consistent with coyote local movement ecology. Distant cross-species transmission and subsequent host movement presents a low risk for onward transmission in raccoon populations. This information helped with emergency response decision-making, thereby saving time and resources.

A genomic basis of vocal rhythm in birds.

Vocal rhythm plays a fundamental role in sexual selection and species recognition in birds, but little is known of its genetic basis due to the confounding effect of vocal learning in model systems. Uncovering its genetic basis could facilitate identifying genes potentially important in speciation. Here we investigate the genomic underpinnings of rhythm in vocal non-learning Pogoniulus tinkerbirds using 135 individual whole genomes distributed across a southern African hybrid zone. We find rhythm speed is associated with two genes that are also known to affect human speech, Neurexin-1 and Coenzyme Q8A. Models leveraging ancestry reveal these candidate loci also impact rhythmic stability, a trait linked with motor performance which is an indicator of quality. Character displacement in rhythmic stability suggests possible reinforcement against hybridization, supported by evidence of asymmetric assortative mating in the species producing faster, more stable rhythms. Because rhythm is omnipresent in animal communication, candidate genes identified here may shape vocal rhythm across birds and other vertebrates.

The gut microbiomes of Channel Island foxes and island spotted skunks exhibit fine-scale differentiation across host species and island populations.

California's Channel Islands are home to two endemic mammalian carnivores: island foxes (Urocyon littoralis) and island spotted skunks (Spilogale gracilis amphiala). Although it is rare for two insular terrestrial carnivores to coexist, these known competitors persist on both Santa Cruz Island and Santa Rosa Island. We hypothesized that examination of their gut microbial communities would provide insight into the factors that enable this coexistence, as microbial symbionts often reflect host evolutionary history and contemporary ecology. Using rectal swabs collected from island foxes and island spotted skunks sampled across both islands, we generated 16S rRNA amplicon sequencing data to characterize their gut microbiomes. While island foxes and island spotted skunks both harbored the core mammalian microbiome, host species explained the largest proportion of variation in the dataset. We further identified intraspecific variation between island populations, with greater differentiation observed between more specialist island spotted skunk populations compared to more generalist island fox populations. This pattern may reflect differences in resource utilization following fine-scale niche differentiation. It may further reflect evolutionary differences regarding the timing of intraspecific separation. Considered together, this study contributes to the growing catalog of wildlife microbiome studies, with important implications for understanding how eco-evolutionary processes enable the coexistence of terrestrial carnivores-and their microbiomes-in island environments.

The value of hybrid genomes: Building two highly contiguous reference genome assemblies to advance Canis genomic studies.

Previous studies of canid population and evolutionary genetics have relied on high-quality domestic dog reference genomes that have been produced primarily for biomedical and trait mapping studies in dog breeds. However, the absence of highly contiguous genomes from other Canis species like the gray wolf and coyote, that represent additional distinct demographic histories, may bias inferences regarding interspecific genetic diversity and phylogenetic relationships. Here, we present single haplotype de novo genome assemblies for the gray wolf and coyote, generated by applying the trio-binning approach to long sequence reads generated from the genome of a female first-generation hybrid produced from a gray wolf and coyote mating. The assemblies were highly contiguous, with contig N50 sizes of 44.6 and 42.0 Mb for the wolf and coyote, respectively. Genome scaffolding and alignments between the two Canis assemblies and published dog reference genomes showed near complete collinearity, with one exception: a coyote-specific chromosome fission of chromosome 13 and fusion of the proximal portion of that chromosome with chromosome 8, retaining the Canis-typical haploid chromosome number of 2n = 78. We evaluated mapping quality for previous RADseq data from 334 canids and found nearly identical mapping quality and patterns among canid species and regional populations regardless of the genome used for alignment (dog, coyote, or gray wolf). These novel wolf and coyote genome reference assemblies will be important resources for proper and accurate inference of Canis demography, taxonomic evaluation, and conservation genetics.

Epigenetic changes to gene pathways linked to male fertility in ex situ black-footed ferrets.

Environmental variation can influence the reproductive success of species managed under human care and in the wild, yet the mechanisms underlying this phenomenon remain largely mysterious. Molecular mechanisms such as epigenetic modifiers are important in mediating the timing and progression of reproduction in humans and model organisms, but few studies have linked epigenetic variation to reproductive fitness in wildlife. Here, we investigated epigenetic variation in black-footed ferrets (Mustela nigripes), an endangered North American mammal reliant on ex situ management for survival and persistence in the wild. Despite similar levels of genetic diversity in human-managed and wild-born populations, individuals in ex situ facilities exhibit reproductive problems, such as poor sperm quality. Differences across these settings suggest that an environmentally driven decline in reproductive capacity may be occurring in this species. We examined the role of DNA methylation, one well-studied epigenetic modifier, in this emergent condition. We leveraged blood, testes, and semen samples from male black-footed ferrets bred in ex situ facilities and found tissue-type specificity in DNA methylation across the genome, although 1360 Gene Ontology terms associated with male average litter size shared functions across tissues. We then constructed gene networks of differentially methylated genomic sites associated with three different reproductive phenotypes to explore the putative biological impact of variation in DNA methylation. Sperm gene networks associated with average litter size and sperm count were functionally enriched for candidate genes involved in reproduction, development, and its regulation through transcriptional repression. We propose that DNA methylation plays an important role in regulating these reproductive phenotypes, thereby impacting the fertility of male ex situ individuals. Our results provide information into how DNA methylation may function in the alteration of reproductive pathways and phenotypes in artificial environments. These findings provide early insights to conservation hurdles faced in the protection of this rare species.

Transposons in the Williams-Beuren Syndrome Critical Region are Associated with Social Behavior in Assistance Dogs.

A strong signature of selection in the domestic dog genome is found in a five-megabase region of chromosome six in which four structural variants derived from transposons have previously been associated with human-oriented social behavior, such as attentional bias to social stimuli and social interest in strangers. To explore these genetic associations in more phenotypic detail-as well as their role in training success in a specialized assistance dog program-we genotyped 1001 assistance dogs from Canine Companions for Independence®, including both successful graduates and dogs released from the training program for behaviors incompatible with their working role. We collected phenotypes on each dog using puppy-raiser questionnaires, trainer questionnaires, and both cognitive and behavioral tests. Using Bayesian mixed models, we found strong associations (95% credibility intervals excluding zero) between genotypes and certain behavioral measures, including separation-related problems, aggression when challenged or corrected, and reactivity to other dogs. Furthermore, we found moderate differences in the genotypes of dogs who graduated versus those who did not; insertions in GTF2I showed the strongest association with training success (ß = 0.23, CI95% = - 0.04, 0.49), translating to an odds-ratio of 1.25 for one insertion. Our results provide insight into the role of each of these four transposons in canine sociability and may inform breeding and training practices for working dog organizations. Furthermore, the observed importance of the gene GTF2I supports the emerging consensus that variation in GTF2I genotypes and expression have important consequences for social behavior broadly.

Demographic history shapes North American gray wolf genomic diversity and informs species' conservation.

Effective population size estimates are critical information needed for evolutionary predictions and conservation decisions. This is particularly true for species with social factors that restrict access to breeding or experience repeated fluctuations in population size across generations. We investigated the genomic estimates of effective population size along with diversity, subdivision, and inbreeding from 162,109 minimally filtered and 81,595 statistically neutral and unlinked SNPs genotyped in 437 grey wolf samples from North America collected between 1986 and 2021. We found genetic structure across North America, represented by three distinct demographic histories of western, central, and eastern regions of the continent. Further, grey wolves in the northern Rocky Mountains have lower genomic diversity than wolves of the western Great Lakes and have declined over time. Effective population size estimates revealed the historical signatures of continental efforts of predator extermination, despite a quarter century of recovery efforts. We are the first to provide molecular estimates of effective population size across distinct grey wolf populations in North America, which ranged between Ne ~ 275 and 3050 since early 1980s. We provide data that inform managers regarding the status and importance of effective population size estimates for grey wolf conservation, which are on average 5.2-9.3% of census estimates for this species. We show that while grey wolves fall above minimum effective population sizes needed to avoid extinction due to inbreeding depression in the short term, they are below sizes predicted to be necessary to avoid long-term risk of extinction.

Patterns of reproduction and autozygosity distinguish the breeding from non-breeding gray wolves of Yellowstone National Park.

For species of management concern, accurate estimates of inbreeding and associated consequences on reproduction are crucial for predicting their future viability. However, few studies have partitioned this aspect of genetic viability with respect to reproduction in a group-living social mammal. We investigated the contributions of foundation stock lineages, putative fitness consequences of inbreeding, and genetic diversity of the breeding versus non-reproductive segment of the Yellowstone National Park gray wolf population. Our dataset spans 25 years and seven generations since reintroduction, encompassing 152 nuclear families and 329 litters. We found over 87% of the pedigree foundation genomes persisted and report influxes of allelic diversity from two translocated wolves from a divergent source in Montana. As expected for group-living species, mean kinship significantly increased over time but with minimal loss of observed heterozygosity. Strikingly, the reproductive portion of the population carried a significantly lower genome-wide inbreeding coefficients, autozygosity, and more rapid decay for linkage disequilibrium relative to the non-breeding population. Breeding wolves had significantly longer lifespans and lower inbreeding coefficients than non-breeding wolves. Our model revealed that the number of litters was negatively significantly associated with heterozygosity (R=-0.11). Our findings highlight genetic contributions to fitness, and the importance of the reproductively active individuals in a population to counteract loss of genetic variation in a wild, free-ranging social carnivore. It is crucial for managers to mitigate factors that significantly reduce effective population size and genetic connectivity, which supports the dispersion of genetic variation that aids in rapid evolutionary responses to environmental challenges.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture.

BACKGROUND: The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function. RESULTS: We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection. CONCLUSIONS: We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Region-specific reversal of epidermal planar polarity in the fancy rosette mouse.

The planar cell polarity (PCP) pathway collectively orients thousands of cells with respect to a body axis to direct cellular behaviors that are essential for embryonic morphogenesis. Hair follicles of the murine epidermis provide a striking readout of PCP activity in their uniform alignment along the entire skin surface. Here, we characterize, from the molecular to tissue-scale, PCP establishment in the rosette fancy mouse, a natural variant with posterior-specific whorls in its fur, to understand how epidermal polarity is coordinated across the tissue. We find that embryonic hair follicles of rosette mutants emerge with reversed orientations specifically in the posterior region, creating a mirror image of epidermal polarity. The rosette trait is associated with a missense mutation in the core PCP gene Fzd6 , which alters a consensus site for N-linked glycosylation and inhibits its membrane localization. Unexpectedly, this defect in Fzd6 trafficking, observed across the entire dorsal epidermis, does not interfere with the ability of other core PCP proteins to localize asymmetrically. Rather, the normally uniform axis of PCP asymmetry is disrupted and rotated in the posterior region such that polarity is reflected on either side of a transition zone. The result is a reversal of polarized cell movements that orient nascent follicles, specifically in the posterior of the embryo. Collectively, our multiscale analysis of epidermal polarity reveals PCP patterning can be regionally decoupled to produce the unique posterior whorls of the fancy rosette mouse. Summary: Region-specific rotation of the Planar Cell Polarity axis reverses posterior hair follicles in the fancy rosette mouse.

Region-specific reversal of epidermal planar polarity in the rosette fancy mouse.

The planar cell polarity (PCP) pathway collectively orients cells with respect to a body axis. Hair follicles of the murine epidermis provide a striking readout of PCP activity in their uniform alignment across the skin. Here, we characterize, from the molecular to tissue-scale, PCP establishment in the rosette fancy mouse, a natural variant with posterior-specific whorls in its fur, to understand how epidermal polarity is coordinated across the tissue. We find that rosette hair follicles emerge with reversed orientations specifically in the posterior region, creating a mirror image of epidermal polarity. The rosette trait is associated with a missense mutation in the core PCP gene Fzd6, which alters a consensus site for N-linked glycosylation, inhibiting its membrane localization. Unexpectedly, the Fzd6 trafficking defect does not block asymmetric localization of the other PCP proteins. Rather, the normally uniform axis of PCP asymmetry rotates where the PCP-directed cell movements that orient follicles are reversed, suggesting the PCP axis rotates 180°. Collectively, our multiscale analysis of epidermal polarity reveals PCP patterning can be regionally decoupled to produce posterior whorls in the rosette fancy mouse.

Response to Hansen Wheat et al.: Additional analysis further supports the early emergence of cooperative communication in dogs compared to wolves raised with more human exposure.

Here, we address Hansen Wheat et al.'s commentary in this journal in response to Salomons et al. Current Biology, 31(14), 3137-3144.E11, (2021). We conduct additional analyses in response to Hansen Wheat et al.'s two main questions. First, we examine the claim that it was the move to a human home environment which enabled the dog puppies to outperform the wolf puppies in gesture comprehension tasks. We show that the youngest dog puppies who had not yet been individually placed in raisers' homes were still highly skilled, and outperformed similar-aged wolf puppies who had higher levels of human interaction. Second, we address the claim that willingness to approach a stranger can explain the difference between dog and wolf pups' ability to succeed in gesture comprehension tasks. We explain the various controls in the original study that render this explanation insufficient, and demonstrate via model comparison that the covariance of species and temperament also make this parsing impossible. Overall, our additional analyses and considerations support the domestication hypothesis as laid out by Salomons et al. Current Biology, 31(14), 3137-3144.E11, (2021).

Ear mite infection restructures otic microbial networks in conservation-reliant Santa Catalina Island foxes (Urocyon littoralis catalinae).

Ceruminous gland tumours are highly prevalent in the ear canals of Santa Catalina Island foxes (Urocyon littoralis catalinae). Previous work suggests that tumours may result from a combination of ectoparasites, disruption of the host-associated microbiome, and host immunopathology. More specifically, ear mite infection has been associated with broad-scale microbial dysbiosis marked by secondary bacterial infection with the opportunistic pathogen Staphylococcus pseudintermedius. Together, ear mites and S. pseudintermedius probably sustain chronic inflammation and promote conditions suitable for tumour development. In the present study, we expanded upon this framework by constructing otic microbial community networks for mite-infected and uninfected foxes sampled in 2017-2019. Across sampling years, we observed consistent signatures of microbial dysbiosis in mite-infected ear canals, including reduced microbial diversity and shifted abundance towards S. pseudintermedius. Network analysis further revealed that mite infection disrupts overall community structure. In mite-infected networks, interaction strengths between taxa were generally weaker, and numerous subnetworks disappeared altogether. We also found that two strains of S. pseudintermedius connected to the main network, suggesting that multistrain biofilm formation may be occurring. In contrast, S. pseudintermedius is peripheral in the uninfected network, with its only connections including a second strain of S. pseudintermedius and the possible competitor Acinetobacter rhizosphaerae. Finally, the lineup of potential keystone taxa shifted across disease states. Fusobacteria spp., a carcinogenesis-promoting microbe, assumed a keystone role in the mite-infected community. Considered together, these findings provide insights into how mite infection may destabilize the microbiome and ultimately contribute to tumour development in this island endemic species.

Pleistocene parades of carnivores into North America.

The distribution and movement of species, broadly known as biogeography, is one of the fundamental subfields of ecology and evolutionary biology. However, significant mysteries remain about the processes that gave rise to the modern distribution of biodiversity across the globe. Over the last several decades, the genetic study of ancient and subfossil specimens has started to shed light on past migrations of some species, with a particular focus on humans and megafauna. In this issue of Molecular Ecology, Salis et al. (2021) use ancient mitogenomes and a new phylogeographic method to add an important new piece of evidence to the mystery of megafaunal migrations into North America during the Pleistocene. They found a striking synchronicity of brown bear (Ursus arctos) and lion (Panthera spp.) migrations across the Bering Land Bridge at several time points during the late Pleistocene, which highlights the lasting impact of sea level change on the prehistoric and modern dispersal of terrestrial carnivores across continents.

Genetic diversity and family groups detected in a coyote population with red wolf ancestry on Galveston Island, Texas.

BACKGROUND: Hybridization can be a conservation concern if genomic introgression leads to the loss of an endangered species' unique genome, or when hybrid offspring are sterile or less fit than their parental species. Yet hybridization can also be an adaptive management tool if rare populations are inbred and have reduced genetic variation, and there is the opportunity to enhance genetic variation through hybridization. The red wolf (Canis rufus) is a critically endangered wolf endemic to the eastern United States, where all extant red wolves are descended from 14 founders which has led to elevated levels of inbreeding over time. Red wolves were considered extirpated from the wild by 1980, but before they disappeared, they interbred with encroaching coyotes creating a genetically admixed population of canids along coastal Texas and Louisiana. In 2018, a genetic study identified individuals on Galveston Island, Texas with significant amounts of red wolf ancestry. We collected 203 fecal samples from Galveston for a more in-depth analysis of this population to identify the amount of red wolf ancestry present and potential mechanisms that support retention of red wolf ancestry on the landscape. RESULTS: We identified 24 individual coyotes from Galveston Island and 8 from mainland Texas with greater than 10% red wolf ancestry. Two of those individuals from mainland Texas had greater than 50% red wolf ancestry estimates. Additionally, this population had 5 private alleles that were absent in the North American reference canid populations used in this study, which included 107 southeastern coyotes, 19 captive red wolves, and 38 gray wolves, possibly representing lost red wolf genetic variation. We also identified several individuals on Galveston Island and the mainland of Texas that retained a unique red wolf mitochondrial haplotype present in the red wolf founding population. On Galveston Island, we identified a minimum of four family groups and found coyotes on the island to be highly related, but not genetically depauperate. We did not find clear associations between red wolf ancestry estimates and landscape features, such as open green space or developed areas. CONCLUSION: Our results confirm the presence of substantial red wolf ancestry persisting on Galveston Island and adjacent mainland Texas. This population has the potential to benefit future red wolf conservation efforts through novel reproductive techniques and possibly through de-introgression strategies, with the goals of recovering extinct red wolf genetic variation and reducing inbreeding within the species.

Disease outbreaks select for mate choice and coat color in wolves.

We know much about pathogen evolution and the emergence of new disease strains, but less about host resistance and how it is signaled to other individuals and subsequently maintained. The cline in frequency of black-coated wolves (Canis lupus) across North America is hypothesized to result from a relationship with canine distemper virus (CDV) outbreaks. We tested this hypothesis using cross-sectional data from wolf populations across North America that vary in the prevalence of CDV and the allele that makes coats black, longitudinal data from Yellowstone National Park, and modeling. We found that the frequency of CDV outbreaks generates fluctuating selection that results in heterozygote advantage that in turn affects the frequency of the black allele, optimal mating behavior, and black wolf cline across the continent.

Coyotes in New York City Carry Variable Genomic Dog Ancestry and Influence Their Interactions with Humans.

Coyotes are ubiquitous on the North American landscape as a result of their recent expansion across the continent. They have been documented in the heart of some of the most urbanized cities, such as Chicago, Los Angeles, and New York City. Here, we explored the genomic composition of 16 coyotes in the New York metropolitan area to investigate genomic demography and admixture for urban-dwelling canids in Queens County, New York. We identified moderate-to-high estimates of relatedness among coyotes living in Queens (r = 0.0-0.5) and adjacent neighborhoods, suggestive of a relatively small population. Although we found low background levels of domestic-dog ancestry across most coyotes in our sample (5%), we identified a male suspected to be a first-generation coyote-dog hybrid with 46% dog ancestry, as well as his two putative backcrossed offspring that carried approximately 25% dog ancestry. The male coyote-dog hybrid and one backcrossed offspring each carried two transposable element insertions that are associated with human-directed hypersociability in dogs and gray wolves. An additional, unrelated coyote with little dog ancestry also carried two of these insertions. These genetic patterns suggest that gene flow from domestic dogs may become an increasingly important consideration as coyotes continue to inhabit metropolitan regions.