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 inter-specific 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 Mb 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 RAD-seq 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.

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.

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.

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.

Reviving ghost alleles: Genetically admixed coyotes along the American Gulf Coast are critical for saving the endangered red wolf.

The last known red wolves were captured in southwestern Louisiana and eastern Texas in 1980 to establish a captive breeding population. Before their extirpation, gene flow with coyotes resulted in the persistence of endangered red wolf genetic variation in local coyote populations. We assessed genomic ancestry and morphology of coyotes in southwestern Louisiana. We detected that 38 to 62% of the coyote genomes contained red wolf ancestry acquired in the past 30 years and have an admixture profile similar to that of the canids captured before the extirpation of red wolves. We further documented a positive correlation between ancestry and weight. Our findings highlight the importance of hybrids and admixed genomes as a reservoir of endangered species ancestry for innovative conservation efforts. Together, this work presents an unprecedented system that conservation can leverage to enrich the recovery program of an endangered species.

The effects of age, sex, weight, and breed on canid methylomes.

Unlike genomes, which are static throughout the lifespan of an organism, DNA methylomes are dynamic. To study these dynamics, we developed quantitative models that measure the effect of multiple factors on DNA methylomes including, age, sex, weight, and genetics. We conducted our study in canids, which prove to be an ideal species to assess epigenetic moderators due to their extreme variability in size and well-characterized genetic structure. We collected buccal swabs from 217 canids (207 domestic dogs and 10 grey wolves) and used targeted bisulphite sequencing to measure methylomes. We also measured genotypes at over one thousand single nucleotide polymorphisms (SNPs). As expected, we found that DNA methylomes are strongly associated with age, enabling the construction of epigenetic clocks. However, we also identify novel associations between methylomes and sex, weight, and sterilization status, leading to accurate models that predict these factors. Methylomes are also affected by genetics, and we observe multiple associations between SNP loci and methylated CpGs. Finally, we show that several factors moderate the relationship between epigenetic ages and real ages, such as body weight, which increases epigenetic ageing. In conclusion, we demonstrate that the plasticity of DNA methylomes is impacted by myriad genetics and physiological factors, and that DNA methylation biomarkers are accurate predictors of age, sex and sterilization status.

Selection of both habitat and genes in specialized and endangered caribou.

Genetic mechanisms determining habitat selection and specialization of individuals within species have been hypothesized, but not tested at the appropriate individual level in nature. In this work, we analyzed habitat selection for 139 GPS-collared caribou belonging to three declining ecotypes sampled throughout Northwestern Canada. We used Resource Selection Functions (RSFs) comparing resources at used and available locations. We found that the three caribou ecotypes differed in their use of habitat suggesting specialization. On expected grounds, we also found differences in habitat selection between summer and winter, but also, originally, among the individuals within an ecotype. We next obtained Single Nucleotide Polymorphisms (SNPs) for the same caribou individuals, we detected those associated to habitat selection, and then identified genes linked to these SNPs. These genes had functions related in other organisms to habitat and dietary specializations, and climatic adaptations. We therefore suggest that individual variation in habitat selection was based on genotypic variation in the SNPs of individual caribou, indicating that genetic forces underlie habitat and diet selection in the species. We also suggest that the associations between habitat and genes that we detected may lead to lack of resilience in the species, thus contributing to caribou endangerment. Our work emphasizes that similar mechanisms may exist for other specialized, endangered species. This article is protected by copyright. All rights reserved.