Search

The bobcat (Lynx rufus) is a medium-sized carnivore well adapted to various environments and an indicator species for landscape connectivity. It is one of the 4 species within the extant Lynx genus in the family Felidae. Because of its broad geographic distribution and central role in food webs, the bobcat is important for conservation. Here we present a high-quality de novo genome assembly of a male bobcat located in Mendocino County, CA, as part of the California Conservation Genomics Project (CCGP). The assembly was generated using the standard CCGP pipeline from a combination of Omni-C and HiFi technologies. The primary assembly comprises 76 scaffolds spanning 2.4 Gb, represented by a scaffold N50 of 142 Mb, a contig N50 of 66.2 Mb, and a BUSCO completeness score of 95.90%. The bobcat genome will be an important resource for the effective management and conservation of this species and comparative genomics exploration.

Species that hibernate generally live longer than would be expected based solely on their body size. Hibernation is characterized by long periods of metabolic suppression (torpor) interspersed by short periods of increased metabolism (arousal). The torpor–arousal cycles occur multiple times during hibernation, and it has been suggested that processes controlling the transition between torpor and arousal states cause ageing suppression. Metabolic rate is also a known correlate of longevity; we thus proposed the ‘hibernation–ageing hypothesis’ whereby ageing is suspended during hibernation. We tested this hypothesis in a well-studied population of yellow-bellied marmots (Marmota flaviventer), which spend 7–8 months per year hibernating. We used two approaches to estimate epigenetic age: the epigenetic clock and the epigenetic pacemaker. Variation in epigenetic age of 149 samples collected throughout the life of 73 females was modelled using generalized additive mixed models (GAMM), where season (cyclic cubic spline) and chronological age (cubic spline) were fixed effects. As expected, the GAMM using epigenetic ages calculated from the epigenetic pacemaker was better able to detect nonlinear patterns in epigenetic ageing over time. We observed a logarithmic curve of epigenetic age with time, where the epigenetic age increased at a higher rate until females reached sexual maturity (two years old). With respect to circannual patterns, the epigenetic age increased during the active season and essentially stalled during the hibernation period. Taken together, our results are consistent with the hibernation–ageing hypothesis and may explain the enhanced longevity in hibernators.

The genetic consequences of species-wide declines are rarely quantified because the timing and extent of the decline varies across the species' range. The sea otter (Enhydra lutris) is a unique model in this regard. Their dramatic decline from thousands to fewer than 100 individuals per population occurred range-wide and nearly simultaneously due to the 18th-19th century fur trade. Consequently, each sea otter population represents an independent natural experiment of recovery after extreme population decline. We designed sequence capture probes for 50 Mb of sea otter exonic and neutral genomic regions. We sequenced 107 sea otters from five populations that span the species range to high coverage (18-76×) and three historical Californian samples from ~1500 and ~200 years ago to low coverage (1.5-3.5×). We observe distinct population structure and find that sea otters in California are the last survivors of a divergent lineage isolated for thousands of years and therefore warrant special conservation concern. We detect signals of extreme population decline in every surviving sea otter population and use this demographic history to design forward-in-time simulations of coding sequence. Our simulations indicate that this decline could lower the fitness of recovering populations for generations. However, the simulations also demonstrate how historically low effective population sizes prior to the fur trade may have mitigated the effects of population decline on genetic health. Our comprehensive approach shows how demographic inference from genomic data, coupled with simulations, allows assessment of extinction risk and different models of recovery.

Long-term studies of mountain lions (Puma concolor) in Southern California have documented persistent small population sizes and the lowest genetic variation of any mountain lion population, except for the Federally endangered mountain lion subspecies, the Florida panther (Puma concolor coryi). There is overwhelming molecular evidence supporting inbreeding and low genetic diversity in these Southern California populations but there is a lack of phenotypical evidence of inbreeding depression. The primary goal of this study was to assess male mountain lions, in Southern California, for teratospermia (>60% abnormal sperm production), one of the first signs of inbreeding depression in mountain lions that are associated with decreased reproduction and population decline. From December 2019 to December 2020, we surveyed mountain lions during live captures, after mortality events, and in images collected from camera traps in the following populations: Santa Monica Mountains, Santa Susana Mountains, Santa Ana Mountains, and the Eastern Peninsular Range. Mountain lions were sampled for known physical abnormalities associated with inbreeding depression such as teratospermia, cryptorchidism, and distal tail kinks. For teratospermia, we extracted testes from five males post-mortem to assess sperm morphology. Epididymal sperm evaluations revealed all males were teratospermic. Across all samples, on average, 93% of observed spermatozoa were abnormal. We physically examined 32 mountain lions (males and females) for distal tail kinks, and we observed one individual affected. We examined 15 male mountain lions for cryptorchidism, and we observed one unilaterally cryptorchid male and one male with testes that differed significantly in size, likely reflecting asynchronous migration of the testes during puberty. Further, we identified three other animals in camera-trap images that had distal tail kinks, for a total of four. In conclusion, from December 2019 to December 2020, we identified nine individuals exhibiting physical signs of inbreeding depression. These reproductive and physical signs of inbreeding depression in Southern California mountain lions increase the urgency of conservation efforts in the region.

Island ecosystems provide models to assess the impacts of isolation on population persistence. However, most studies of persistence have focused on a single species, without comparisons to other organisms they interact with in the ecosystem. The simple predator-prey system of moose and gray wolves on Isle Royale provides allows a direct contrast of genetic variation in a prey species with their natural predator. Wolves on Isle Royale exhibited signs of severe inbreeding depression, which nearly drove the population to extinction in 2019. In the relative absence of wolves, the moose population has thrived and exhibits no obvious signs of inbreeding depression despite being isolated for ∼120 years and having low genetic diversity. Here, we examine the genomic underpinnings of population persistence in the Isle Royale moose population. We document high levels of inbreeding in the population, roughly as high as the wolf population at the time of its decline. However, inbreeding in the moose population manifests in the form of intermediate-length runs of homozygosity indicative of gradual inbreeding, contrasting with the severe recent inbreeding observed in the wolf population. Using simulations, we demonstrate that this more gradual inbreeding in the moose population has resulted in an estimated 50% purging of the inbreeding load, helping to explain the continued persistence of the population. However, we also document notable increases in genetic load, which could eventually threaten population viability over the long term. Finally, we document low diversity in mainland North American moose populations due to a severe founder event occurring near the end of the Holocene. Overall, our results demonstrate a complex relationship between inbreeding, genetic diversity, and population viability that highlights the importance of maintaining isolated populations at moderate size to avert extinction from genetic factors.Significance statementIsolated wildlife populations face a high risk of extinction due in part to the deleterious consequences of inbreeding. Whether purifying natural selection can overcome these negative impacts by “purging” harmful recessive mutations is a topic of active debate. We characterized the extent of purging in an isolated moose population. Our results demonstrate signatures of gradual inbreeding in the population, ideal circumstances to facilitate purging. Using simulations, we demonstrate substantial potential for purging in the population, though we also show that fitness is reduced by small population size and inbreeding. Our findings provide insight into the mechanisms enabling persistence in isolated populations, with implications for conserving the growing number of isolated populations worldwide.

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.

The remarkable radiation of South American (SA) canids produced 10 extant species distributed across diverse habitats, including disparate forms such as the short-legged, hypercarnivorous bush dog and the long-legged, largely frugivorous maned wolf. Despite considerable research spanning nearly two centuries, many aspects of their evolutionary history remain unknown. Here, we analyzed 31 whole genomes encompassing all extant SA canid species to assess phylogenetic relationships, interspecific hybridization, historical demography, current genetic diversity, and the molecular bases of adaptations in the bush dog and maned wolf. We found that SA canids originated from a single ancestor that colonized South America 3.9 to 3.5 Mya, followed by diversification east of the Andes and then a single colonization event and radiation of Lycalopex species west of the Andes. We detected extensive historical gene flow between recently diverged lineages and observed distinct patterns of genomic diversity and demographic history in SA canids, likely induced by past climatic cycles compounded by human-induced population declines. Genome-wide scans of selection showed that disparate limb proportions in the bush dog and maned wolf may derive from mutations in genes regulating chondrocyte proliferation and enlargement. Further, frugivory in the maned wolf may have been enabled by variants in genes associated with energy intake from short-chain fatty acids. In contrast, unique genetic variants detected in the bush dog may underlie interdigital webbing and dental adaptations for hypercarnivory. Our analyses shed light on the evolution of a unique carnivoran radiation and how it was shaped by South American topography and climate change.

Climate change is leading to habitat shifts that threaten species persistence throughout California's unique ecosystems. Baseline biodiversity data would provide opportunities for habitats to be managed under short-term and long-term environmental change. Aiming to provide biodiversity data, the UC Conservation Genomics Consortium launched the California Environmental DNA (CALeDNA) program to be a citizen and community science biomonitoring initiative that uses environmental DNA (eDNA, DNA shed from organisms such as from fur, feces, spores, pollen or leaves). Now with results from 1,000 samples shared online, California biodiversity patterns are discoverable. Soil, sediment and water collected by researchers, undergraduates and the public reveal a new catalog of thousands of organisms that only slightly overlap with traditional survey bioinventories. The CALeDNA website lets users explore the taxonomic diversity in different ways, and researchers have created tools to help people new to eDNA to analyze community ecology patterns. Although eDNA results are not always precise, the program team is making progress to fit it into California's biodiversity management toolbox, such as for monitoring ecosystem recovery after invasive species removal or wildfire.

Human-driven habitat fragmentation and loss have led to a proliferation of small and isolated plant and animal populations with high risk of extinction. One of the main threats to extinction in these populations is inbreeding depression, which is primarily caused by recessive deleterious mutations becoming homozygous due to inbreeding. The typical approach for managing these populations is to maintain high genetic diversity, increasingly by translocating individuals from large populations to initiate a “genetic rescue.” However, the limitations of this approach have recently been highlighted by the demise of the gray wolf population on Isle Royale, which declined to the brink of extinction soon after the arrival of a migrant from the large mainland wolf population. Here, we use a novel population genetic simulation framework to investigate the role of genetic diversity, deleterious variation, and demographic history in mediating extinction risk due to inbreeding depression in small populations. We show that, under realistic models of dominance, large populations harbor high levels of recessive strongly deleterious variation due to these mutations being hidden from selection in the heterozygous state. As a result, when large populations contract, they experience a substantially elevated risk of extinction after these strongly deleterious mutations are exposed by inbreeding. Moreover, we demonstrate that, although genetic rescue is broadly effective as a means to reduce extinction risk, its effectiveness can be greatly increased by drawing migrants from small or moderate-sized source populations rather than large source populations due to smaller populations harboring lower levels of recessive strongly deleterious variation. Our findings challenge the traditional conservation paradigm that focuses on maximizing genetic diversity in small populations in favor of a view that emphasizes minimizing strongly deleterious variation. These insights have important implications for managing small and isolated populations in the increasingly fragmented landscape of the Anthropocene.

International audience; Aim. The majority of work done to gather information on the Earth's biodiversity has been carried out using in-situ data, with known issues related to epistemology (e.g., species determination and taxonomy), spatial uncertainty, logistics (time and costs), among others. An alternative way to gather information about spatial ecosystem variability is the use of satellite remote sensing. It works as a powerful tool for attaining rapid and standardized information. Several metrics used to calculate remotely sensed diversity of ecosystems are based on Shannon’s information theory, namely on the differences in relative abundance of pixel reflectances in a certain area. Additional metrics like the Rao’s quadratic entropy allow the use of spectral distance beside abundance, but they are point descriptors of diversity, that is they can account only for a part of the whole diversity continuum. The aim of this paper is thus to generalize the Rao’s quadratic entropy by proposing its parameterization for the first time. Innovation. The parametric Rao’s quadratic entropy, coded in R, (a) allows the representation of the whole continuum of potential diversity indices in one formula, and (b) starting from the Rao’s quadratic entropy, allows the explicit use of distances among pixel reflectance values, together with relative abundances. Main conclusions. The proposed unifying measure is an integration between abundance- and distance-based algorithms to map the continuum of diversity given a satellite image at any spatial scale. Being part of the rasterdiv R package, the proposed method is expected to ensure high robustness and reproducibility.

Biodiversity monitoring in conservation projects is essential to understand environmental status and recovery. However, traditional field surveys can be expensive, time-consuming, biased towards visual detection, and focused on measuring a limited set of taxa. Environmental DNA (eDNA) methods provide a new approach to biodiversity monitoring that has the potential to sample a taxonomically broader set of organisms with similar effort, but many of these approaches are still in the early stages of development and testing. Here, we use multilocus eDNA metabarcoding to understand how the removal of invasive red swamp crayfish impacts local biodiversity of a desert oasis ecosystem, as well as to detect crayfish both directly and indirectly. We tracked crayfish DNA signatures, microbial DNA associated with crayfish, and biodiversity of plant, fungal, animal, and bacterial communities through time. We were unsuccessful in detecting crayfish directly in either control tanks or oases using targeted metabarcoding primers for invertebrates and eukaryotes, similar to previous studies which have shown variable levels of success in detecting crayfish from environmental samples. However, we were successful in discerning a suite of 90 crayfish-associated taxa to serve as candidate bioindicators of invasive presence using 16S and Fungal ITS2 metabarcoding. Ranking these 90 taxa by their geographic distribution in eDNA surveys and by evidence of crayfish-associations in the literature, we support 9 taxa to be high-ranking, and suggest they be prioritized in future biomonitoring. Biodiversity analyses from five metabarcode loci including plants, animals, and both prokaryotic and eukaryotic microbes showed that communities differed but that species richness remained relatively similar between oases through time. Our results reveal that, while there are limitations of eDNA approaches to detect crayfish and other invasive species, microbial bioindicators offer a largely untapped biomonitoring opportunity for invasive species management, adding a valuable resource to a conservation manager’s toolkit.

In cases of severe wildlife population decline, a key question is whether recovery efforts will be impeded by genetic factors, such as inbreeding depression. Decades of excess mortality from gillnet fishing have driven Mexico’s vaquita porpoise ( Phocoena sinus ) to ~10 remaining individuals. We analyzed whole-genome sequences from 20 vaquitas and integrated genomic and demographic information into stochastic, individual-based simulations to quantify the species’ recovery potential. Our analysis suggests that the vaquita’s historical rarity has resulted in a low burden of segregating deleterious variation, reducing the risk of inbreeding depression. Similarly, genome-informed simulations suggest that the vaquita can recover if bycatch mortality is immediately halted. This study provides hope for vaquitas and other naturally rare endangered species and highlights the utility of genomics in predicting extinction risk.

Mojave Desert springs are fragile ecosystems, hosting endemic plants and animals, which are threatened by the increasing human demand for water and climate change. To develop management practices that will protect the groundwater‐dependent ecosystems at Mojave Desert springs, real‐time, low‐cost biodiversity monitoring, and assessments are required. Environmental DNA (eDNA) metabarcoding uses DNA shed from organisms (e.g., skin cells, feces, and pollen) that are present in water, air, soil, or sediment samples to assess community composition. This approach can increase the detection sensitivity for rare and elusive species, compared with expensive and time‐consuming conventional methods, which also require taxonomic expertise. This study tests the effectiveness of eDNA techniques in capturing the observed Mojave Desert spring biodiversity in the winter and spring of 2019 at four distinct, naturally occurring springs. We also test the utility of sample types (water vs. sediment) for capturing biodiversity. We found that each of the four Mojave Desert springs supports a unique biological community. Sediment samples contained the greatest biodiversity, but all sample types captured species observed in the field by humans or camera traps. We also found no statistical difference in species richness captured in winter and spring except for the Cytochrome Oxidase I marker, for which winter had greater biodiversity. This study supports the use of eDNA metabarcoding as an effective tool to mirror observation by human observers of ecological communities in desert springs. The study demonstrates the importance of appropriately timing eDNA field sampling, primer selection, and using field‐based surveys of wildlife and plants in addition to eDNA detection. This study also identified gaps in reference sequence databases for Mojave biodiversity and encourages collaboration of eDNA researchers with managers for effective conservation management plans.

Ecosystems globally are under threat from ongoing anthropogenic environmental change. Effective conservation management requires more thorough biodiversity surveys that can reveal system-level patterns and that can be applied rapidly across space and time. Using modern ecological models and community science, we integrate environmental DNA and Earth observations to produce a time snapshot of regional biodiversity patterns and provide multi-scalar community-level characterization. We collected 278 samples in spring 2017 from coastal, shrub, and lowland forest sites in California, a complex ecosystem and biodiversity hotspot. We recovered 16,118 taxonomic entries from eDNA analyses and compiled associated traditional observations and environmental data to assess how well they predicted alpha, beta, and zeta diversity. We found that local habitat classification was diagnostic of community composition and distinct communities and organisms in different kingdoms are predicted by different environmental variables. Nonetheless, gradient forest models of 915 families recovered by eDNA analysis and using BIOCLIM variables, Sentinel-2 satellite data, human impact, and topographical features as predictors, explained 35% of the variance in community turnover. Elevation, sand percentage, and photosynthetic activities (NDVI32) were the top predictors. In addition to this signal of environmental filtering, we found a positive relationship between environmentally predicted families and their numbers of biotic interactions, suggesting environmental change could have a disproportionate effect on community networks. Together, these analyses show that coupling eDNA with environmental predictors including remote sensing data has capacity to test proposed Essential Biodiversity Variables and create new landscape biodiversity baselines that span the tree of life.

One of the most enduring surprises about the genetic history of Late Pleistocene populations is that continuity is often disturbed by upheaval. In fact, studies that support population continuity are increasingly rare in humans, a variety of vertebrate taxa, and vascular plants (Hofreiter & Stewart 2009; Burbrink et al. 2016). Perhaps such continuity should not be expected as the Pleistocene is marked by episodes of climate change, glaciation and the invasions of humans into previously isolated areas. Although fossils are one of the primary sources for inferring population continuity, a problem with fossil material is that, even if similar morphological forms might exist in a place over time, they may not be from the same genetic lineage. There are now readily available methods to assess genetic continuity solely from DNA found in fossil material, provided the record is fairly continuous. In a From the Cover article in this issue of Molecular Ecology, Loog et al. (2020) apply some of these readily available methods to analyse mitochondrial genomes and model the demography of wolves over the last 50,000 years.

During the late Pleistocene, isolated lineages of hominins exchanged genes thus influencing genomic variation in humans in both the past and present. However, the dynamics of this genetic exchange and associated phenotypic consequences through time remain poorly understood. Gene exchange across divergent lineages can result in myriad outcomes arising from these dynamics and the environmental conditions under which it occurs. Here we draw from our collective research across various organisms, illustrating some of the ways in which gene exchange can structure genomic/phenotypic diversity within/among species. We present a range of examples relevant to questions about the evolution of hominins. These examples are not meant to be exhaustive, but rather illustrative of the diverse evolutionary causes/consequences of hybridization, highlighting potential drivers of human evolution in the context of hybridization including: influences on adaptive evolution, climate change, developmental systems, sex-differences in behavior, Haldane’s rule and the large X-effect, and transgressive phenotypic variation.

Within the Canidae, the African wild dog (Lycaon pictus) is the most specialized with regards to cursorial adaptations (specialized for running), having only four digits on their forefeet. In addition, this species is one of the few canids considered to be an obligate meat-eater, possessing a robust dentition for taking down large prey, and displays one of the most variable coat colorations amongst mammals. Here, we used comparative genomic analysis to investigate the evolutionary history and genetic basis for adaptations associated with cursoriality, hypercanivory, and coat color variation in African wild dogs. Genome-wide scans revealed unique amino acid deletions that suggest a mode of evolutionary digit loss through expanded apoptosis in the developing first digit. African wild dog-specific signals of positive selection also uncovered a putative mechanism of molar cusp modification through changes in genes associated with the sonic hedgehog (SHH) signaling pathway, required for spatial patterning of teeth, and three genes associated with pigmentation. Divergence time analyses suggest the suite of genomic changes we identified evolved ~1.7 Mya, coinciding with the diversification of large-bodied ungulates. Our results show that comparative genomics is a powerful tool for identifying the genetic basis of evolutionary changes in Canidae.

Psychological stress induced by exposure to predators has complex effects on the behaviour and physiology of prey species. This includes potential influences on gene expression mediated via stress-responsive physiological pathways such as the sympathetic nervous system and hypothalamus–pituitary–adrenal (HPA) axis. Laboratory studies have documented diverse transcriptional effects of predator-induced fear, but genomic responses to predator exposure in the wild remain poorly understood. Here, we used RNA-sequencing to investigate the leukocyte transcriptome response to chronic predator pressure in a well-studied population of wild yellow-bellied marmots, Marmota flaviventer. We assessed the genomic response to this stressor in three ways by (1) identifying differentially expressed individual genes across the genome, (2) assessing whether differentially expressed genes were statistically over-represented by functional categories and (3) testing for transcription factor activity that may mediate observed gene expression differences. We found 349 individual genes regulated in association with chronic predator presence, including transcripts known to regulate heat shock proteins, metabolism and DNA damage repair. Gene ontology analysis revealed that the majority of these differentially expressed genes were involved with the cellular response to stress, cellular metabolism and protein transport. Transcription factor analysis implicated glucocorticoid signalling in mediating these effects. Our work confirms that the physiological response to predator-induced stress is complex, initiating transcriptional activity in multiple processes and pathways. In addition to the canonical expectations that individuals exposed to predators mobilize HPA signalling and homeostasis pathways, we also detected activity in genes typically associated with human anxiety and cerebral function. This is the first study to demonstrate that leukocyte transcriptomes taken from animals living in a natural environment can reflect the complex ecology of fear.

Anticoagulant rodenticides (ARs) are used worldwide to control rodent populations. ARs bioaccumulate across trophic levels and threaten non-target wildlife. We investigated the prevalence of AR exposure in seven predator species in the rapidly developing Greater Cape Town region of South Africa – a mosaic of natural, urban, and agricultural areas within a global biodiversity hotspot. We focused sampling on caracals (Caracal caracal, n = 28) as part of a larger caracal ecology study, but also opportunistically sampled Cape Clawless otters (Aonyx capensis, n = 9), large-spotted genets (Genetta tigrina, n = 4), honey badger (Mellivora capensis, n = 1), water mongoose (Atilax paludinosus, n = 1), small gray mongoose (Galerella pulverulenta, n = 1), and Cape Eagle owl (Bubo capensis, n = 1). We tested livers from all species, and blood from ten caracals, for eight AR compounds to assess prevalence and amount of exposure for each compound. We used generalized linear models to test spatial, demographic, and seasonal risk factors for ten measures of AR exposure in caracals. We detected at least one of the four most toxic AR compounds in six species. Exposure was high for caracals (92%) and all species combined (81%). For caracals, proximity to vineyards was the most important AR exposure risk factor. Vineyards in Cape Town do not use ARs to protect their vines but do host commercial hospitality structures where ARs are used. Vineyards may thus link caracals that forage within vineyards to the rat poisons used in and around their commercial structures. Residue levels were unexpected in large-spotted genets and Cape Clawless otters, suggesting invertebrate vectors. ARs may present a cryptic threat to populations already vulnerable to increasing habitat loss, vehicle collisions, poachers and fire. Targeted mitigation should include a mix of environmentally responsible policies that reduce AR use, particularly in areas near wildlife habitat.

White-nosed coatis (Nasua narica) are widely distributed throughout North, Central, and South America, but the patterns of temporal and spatial diversification that have contributed to this distribution are unknown. In addition, the biogeographic history of procyonid species in the Americas remains contentious. Using sequences from three mitochondrial loci (Cytochrome b, NAHD5 and 16S rRNA; 2201 bp) and genotypes from 11 microsatellite loci, we analyzed genetic diversity to determine phylogeographic patterns, genetic structure, divergence times, and gene flow among Nasua narica populations throughout the majority of the species' range. We also estimated the ancestral geographic range of N. narica and other procyonid species. We found a high degree of genetic structure and divergence among populations that conform to five evolutionarily significant units. The most southerly distributed population (Panama) branched off much earlier (∼3.8 million years ago) than the northern populations (1.2 million years ago). Estimated gene flow among populations was low and mostly northwards and westwards. The phylogeographic patterns within N. narica are associated with geographic barriers and habitat shifts likely caused by Pliocene-Pleistocene climate oscillations. Significantly, our findings suggest the dispersal of N. narica was south-to-north beginning in the Pliocene, not in the opposite direction during the Pleistocene as suggested by the fossil record, and that the most recent common ancestor for coati species was most likely distributed in South or Central America six million years ago. Our study implies the possibility that the diversification of Nasua species, and other extant procyonid lineages, may have occurred in South America.

Species that hibernate live longer than would be expected based solely on their body size. Hibernation is characterized by long periods of metabolic suppression (torpor) interspersed by short periods of increased metabolism (arousal). The torpor-arousal cycles occur multiple times during hibernation, and it has been suggested that processes controlling the transition between torpor and arousal states cause aging suppression. Metabolic rate is also a known correlate of longevity, we thus proposed the ‘hibernation-aging hypothesis’ whereby aging is suspended during hibernation. We tested this hypothesis in a well-studied population of yellow-bellied marmots (Marmota flaviventer), which spend 7-8 months per year hibernating. We used two approaches to estimate epigenetic age: the epigenetic clock and the epigenetic pacemaker. Variation in epigenetic age of 149 samples collected throughout the life of 73 females was modeled using generalized additive mixed models (GAMM), where season (cyclic cubic spline) and chronological age (cubic spline) were fixed effects. As expected, the GAMM using epigenetic ages calculated from the epigenetic pacemaker was better able to detect nonlinear patterns in epigenetic age change over time. We observed a logarithmic curve of epigenetic age with time, where the epigenetic age increased at a higher rate until females reached sexual maturity (2-years old). With respect to circannual patterns, the epigenetic age increased during the summer and essentially stalled during the winter. Our enrichment analysis of age-related CpG sites revealed pathways related to development and cell differentiation, while the season-related CpGs enriched pathways related to central carbon metabolism, immune system, and circadian clock. Taken together, our results are consistent with the hibernation-aging hypothesis and may explain the enhanced longevity in hibernators.

Species that hibernate generally live longer than would be expected based solely on their body size. Hibernation is characterized by long periods of metabolic suppression (torpor) interspersed by short periods of increased metabolism (arousal). The torpor-arousal cycles occur multiple times during hibernation, and it has been suggested that processes controlling the transition between torpor and arousal states cause ageing suppression. Metabolic rate is also a known correlate of longevity; we thus proposed the 'hibernation-ageing hypothesis' whereby ageing is suspended during hibernation. We tested this hypothesis in a well-studied population of yellow-bellied marmots (Marmota flaviventer), which spend 7-8 months per year hibernating. We used two approaches to estimate epigenetic age: the epigenetic clock and the epigenetic pacemaker. Variation in epigenetic age of 149 samples collected throughout the life of 73 females was modelled using generalized additive mixed models (GAMM), where season (cyclic cubic spline) and chronological age (cubic spline) were fixed effects. As expected, the GAMM using epigenetic ages calculated from the epigenetic pacemaker was better able to detect nonlinear patterns in epigenetic ageing over time. We observed a logarithmic curve of epigenetic age with time, where the epigenetic age increased at a higher rate until females reached sexual maturity (two years old). With respect to circannual patterns, the epigenetic age increased during the active season and essentially stalled during the hibernation period. Taken together, our results are consistent with the hibernation-ageing hypothesis and may explain the enhanced longevity in hibernators.

Species that hibernate live longer than would be expected based solely on their body size. Hibernation is characterized by long periods of metabolic suppression (torpor) interspersed by short periods of increased metabolism (arousal). The torpor-arousal cycles occur multiple times during hibernation, and it has been suggested that processes controlling the transition between torpor and arousal states cause aging suppression. Metabolic rate is also a known correlate of longevity, we thus proposed the ‘hibernation-aging hypothesis’ whereby aging is suspended during hibernation. We tested this hypothesis in a well-studied population of yellow-bellied marmots (Marmota flaviventer), which spend 7-8 months per year hibernating. We used two approaches to estimate epigenetic age: the epigenetic clock and the epigenetic pacemaker. Variation in epigenetic age of 149 samples collected throughout the life of 73 females were modeled using generalized additive mixed models (GAMM), where season (cyclic cubic spline) and chronological age (cubic spline) were fixed effects. As expected, the GAMM using epigenetic ages calculated from the epigenetic pacemaker was better able to detect nonlinear patterns in epigenetic age change over time. We observed a logarithmic curve of epigenetic age with time, where the epigenetic age increased at a higher rate until females reached sexual maturity (2-years old). With respect to circannual patterns, the epigenetic age increased during the summer and essentially stalled during the winter. Our enrichment analysis of age-related CpG sites revealed pathways related to development and cell differentiation, while the season-related CpGs enriched pathways related to central carbon metabolism, immune system, and circadian clock. Taken together, our results are consistent with the hibernation-aging hypothesis and may explain the enhanced longevity in hibernators.

AimThe majority of work done to gather information on Earth diversity has been carried out by in-situ data, with known issues related to epistemology (e.g., species determination and taxonomy), spatial uncertainty, logistics (time and costs), among others. An alternative way to gather information about spatial ecosystem variability is the use of satellite remote sensing. It works as a powerful tool for attaining rapid and standardized information. Several metrics used to calculate remotely sensed diversity of ecosystems are based on Shannon’s Information Theory, namely on the differences in relative abundance of pixel reflectances in a certain area. Additional metrics like the Rao’s quadratic entropy allow the use of spectral distance beside abundance, but they are point descriptors of diversity, namely they can account only for a part of the whole diversity continuum. The aim of this paper is thus to generalize the Rao’s quadratic entropy by proposing its parameterization for the first time.InnovationThe parametric Rao’s quadratic entropy, coded in R, i) allows to represent the whole continuum of potential diversity indices in one formula, and ii) starting from the Rao’s quadratic entropy, allows to explicitly make use of distances among pixel reflectance values, together with relative abundances.Main conclusionsThe proposed unifying measure is an integration between abundance- and distance-based algorithms to map the continuum of diversity given a satellite image at any spatial scale.

In North American gray wolves, black coat color is dominantly inherited via a 3 base pair coding deletion in the canine beta defensin 3 (CBD103) gene. This 3 base pair deletion, called the KB allele, was introduced through hybridization with dogs and subsequently underwent a selective sweep that increased its frequency in wild wolves. Despite apparent positive selection, KBB wolves have lower fitness than wolves with the KyB genotype, even though the 2 genotypes show no observable differences in black coat color. Thus, the KB allele is thought to have pleiotropic effects on as-yet unknown phenotypes. Given the role of skin-expressed CBD103 in innate immunity, we hypothesized that the KB allele influences the keratinocyte gene expression response to TLR3 pathway stimulation and/or infection by canine distemper virus (CDV). To test this hypothesis, we developed a panel of primary epidermal keratinocyte cell cultures from 24 wild North American gray wolves of both Kyy and KyB genotypes. In addition, we generated an immortalized Kyy line and used CRISPR/Cas9 editing to produce a KyB line on the same genetic background. We assessed the transcriptome-wide responses of wolf keratinocytes to the TLR3 agonist polyinosinic:polycytidylic acid (polyI:C), and to live CDV. K locus genotype did not predict the transcriptional response to either challenge, suggesting that variation in the gene expression response does not explain pleiotropic effects of the KB allele on fitness. This study supports the feasibility of using cell culture methods to investigate the phenotypic effects of naturally occurring genetic variation in wild mammals.

In North American gray wolves, black coat color is dominantly inherited via a three base pair coding deletion in the canine beta defensin 3 (CBD103) gene. This three base pair deletion, called the KB allele, was introduced through hybridization with dogs and subsequently underwent a selective sweep that increased its frequency in wild wolves. Despite apparent positive selection, KBB wolves have significantly lower fitness than wolves with the KyB genotype, even though the two genotypes show no observable differences in black coat color. Thus, the KB allele is thought to have pleiotropic effects on as-yet unknown phenotypes. Given the role of skin-expressed CBD103 in innate immunity, we hypothesized that the KB allele influences the gene regulatory response to viral infection. To test this hypothesis, we developed a panel of primary epidermal keratinocyte cell cultures from 24 wild North American gray wolves (both Kyy and KyB genotypes) and generated immortalized Kyy and CBD103 knockout lines. We assessed the transcriptome-wide responses of wolf keratinocytes to polyinosinic:polycytidylic acid (polyI:C), which mimics infection by a double-stranded RNA virus, and to live canine distemper virus. Keratinocytes with the KyB genotype and with the Kyy genotype had similar gene regulatory responses to viral infection, suggesting that this response does not explain pleiotropic effects of the KB allele on fitness. This study supports the feasibility of using cell culture methods to investigate the phenotypic effects of naturally occurring genetic variation in wild mammals.

The most remarkable feature of our planet is the diversity of its life forms, ranging from viruses and nanobacteria to blue whales and giant sequoias to satanic leaf‐tailed geckos and leafy seadragons (look them up!). Life is found in essentially all environments on earth, and the number of species living on our planet is many times greater than we could have imagined a century ago. A well‐regarded estimate pegs the number of eukaryotic species on earth at 8.7 million (±1.3 million), of which fewer than 15% are currently described (Mora et al., 2011). The diversity of prokaryotes is less clear (and highly controversial), but an analysis of 1.6 billion 16S ribosomal RNA sequences estimated that 0.8–1.6 million prokaryotic operational taxonomic units exist globally (Louca et al., 2019). While we do not know how many species are currently extant, or have existed in the past, we do know that this biodiversity is valuable, providing food, fibre and medicine, furnishing ecosystem services such as water and air purification, nutrient cycling, pollination and carbon uptake, and contributing to technological innovations ranging from biotechnology to robotics to material science. Moreover, biodiversity underlies the cultural identity of human populations and is important to human health and well‐being. Geographically, species richness increases from the Polar Regions to the tropics in terrestrial and surface marine ecosystems. Thus, some countries, especially those in tropical and subtropical regions, are endowed with much greater biodiversity than others. Unfortunately, benefits arising from the access and utilization of this biodiversity have been unequally shared, with (paradoxically) biodiversity‐poor countries often accruing the lion's share of economic gains. There can be imbalances within countries as well, wherein some segments of the population obtain greater economic benefits from biodiversity and associated traditional knowledge than indigenous peoples. The “Nagoya Protocol on Access to Genetic Resources and the Fair and Equitable Sharing of Benefits Arising from their Utilization,” which came into force in 2014, is an international agreement designed to ensure that the benefits arising from biodiversity are shared equitably (https://www.cbd.int/abs/). However, few scientific journals require compliance with the Nagoya Protocol or the reporting of benefits from biodiversity research. In this editorial, we (the editors of Molecular Ecology and Molecular Ecology Resources) express support for the Nagoya Protocol and the principle of benefit sharing. We believe that scientific journals publishing research on biodiversity can play an important role in implementing the Nagoya Protocol and in reporting on benefits generated from such research. Below, we provide background on the Nagoya Protocol, discuss the kinds of benefits that may arise from biodiversity research, describe the rationale for reporting on these benefits and introduce changes to the journals’ Data Accessibility Statements to incorporate the requirements and goals of the Nagoya Protocol.

Transitioning from in‐person to remote learning can present challenges for both the instructional team and the students. Here, we use our course “Biodiversity in the Age of Humans” to describe how we adapted tools and strategies designed for a flipped classroom to a remote learning format. Using anonymous survey data collected from students who attended the course either in‐person (2019) or remotely (2020), we quantify student expectations and experiences and compare these between years. We summarize our experience and provide ten “tips” or recommendations for a transition to remote learning, which we divide into three categories: (a) precourse instructor preparation; (b) outside of class use of online materials; and (c) during class student engagement. The survey results indicated no negative impact on student learning during the remote course compared to in‐person instruction. We found that communicating with students and assessing specific needs, such as access to technology, and being flexible with the structure of the course, simplified the transition to remote instruction. We also found that short, pre‐recorded videos that introduce subject materials were among the most valuable elements for student learning. We hope that instructors of undergraduate ecology and evolution courses can use these recommendations to help establish inclusive online learning communities that empower students to acquire conceptual knowledge and develop scientific inquiry and literacy skills., Transitioning from in‐person to remote learning can present challenges for both the instructional team and the students. Here, we use our course “Biodiversity in the Age of Humans” to describe how we adapted tools and strategies designed for a flipped classroom to remote learning. We provide ten “tips,” which we divide into three categories: (a) precourse instructor preparation; (b) outside of class use of online materials; and (c) during class student engagement. We hope that instructors of undergraduate ecology and evolution courses can use these recommendations to help establish inclusive online learning communities that empower students.

Large percentages of abnormal sperm, termed teratospermia, are associated with poor fertility in cats, many of which are threatened with extinction from their natural habitats. Even normal appearing spermatozoa from felids with teratospermia may have a compromised capacity for motility and fertilization indicating there are factors affecting the fertilization capacity of all sperm regardless of morphology. There was a comparative study conducted using the RNA-Seq approach to identify differentially expressed genes between morphologically normal and abnormal sperm from domestic cates with normospermia and teratospermia to elucidate genes and pathways associated with abnormal sperm function. Normal sperm from cats with teratospermia have a gene expression profile similar to abnormal sperm from males with teratospermia. There was also downregulation of cGMP pathways which may be associated with a lesser sperm motility in ejaculates from males with teratospermia. Kinase phosphorylation pathways also were downregulated in normal spermatozoa from ejaculates of males with teratospermia. Results indicate that analysis of sperm gene expression provides for a more precise assessment of sperm function in semen of cats with teratospermia and facilitates identification of molecular abnormalities that may lead to compromised fertilization capacity.

Dire wolves are considered to be one of the most common and widespread large carnivores in Pleistocene America1, yet relatively little is known about their evolution or extinction. Here, to reconstruct the evolutionary history of dire wolves, we sequenced five genomes from sub-fossil remains dating from 13,000 to more than 50,000 years ago. Our results indicate that although they were similar morphologically to the extant grey wolf, dire wolves were a highly divergent lineage that split from living canids around 5.7 million years ago. In contrast to numerous examples of hybridization across Canidae2,3, there is no evidence for gene flow between dire wolves and either North American grey wolves or coyotes. This suggests that dire wolves evolved in isolation from the Pleistocene ancestors of these species. Our results also support an early New World origin of dire wolves, while the ancestors of grey wolves, coyotes and dholes evolved in Eurasia and colonized North America only relatively recently. © 2021, The Author(s), under exclusive licence to Springer Nature Limited. Acknowledgements We thank the staff at the Carnegie Museum of Natural History, Cincinnati Museum Center, Danish Zoological Museum, Harrison Zoological Museum, Harvard Museum of Comparative Zoology, Idaho Museum of Natural History, Institute of Archaeology (Russian Academy of Sciences), Institute of Systematics and Animal Ecology (Russian Academy of Sciences), Institute of Zoology (Chinese Academy of Sciences), Instituto de Conservação da Natureza e das Florestas, Kansas Museum of Natural History, La Brea Tar Pits and Museum, Ludwig Maximilian University, McClung Museum, Museum of the Institute of Plant and Animal Ecology (Russian Academy of Sciences), Museum national d’Histoire naturelle, National Museums Scotland, Natural History Museum London, Naturalis Biodiversity Center, Naturhistorisches Museum Bern, Smithsonian National Museum of Natural History, Swedish Naturhistoriska Riksmuseet, SYLVATROP, US Bureau of Reclamation, University of California Museum of Paleontology, University of Texas at El Paso, University of Washington Burke Museum and the Zoological Institute (Russian Academy of Sciences; state assignment no. АААА-А19-119032590102-7) for access to specimens in their care; T. Barnosky, S. Bray, A. Farrell, R. Fischer, A. Harris, J. Harris, A. Henrici, P. Holroyd, R. MacPhee, T. Martin, A. Philpot, J. Saunders, J. Southon, G. Storrs, G. Takeuchi, X. Wang and C. Widga for assistance; and L. DeSantis for comments. A.M. used computational and storage services associated with the Hoffman2 Shared Cluster provided by UCLA Institute for Digital Research and Education’s Research Technology Group. DireGWC was sequenced using the Vincent J. Coates Genomics Sequencing Laboratory at UC Berkeley, supported by NIH S10 OD018174 Instrumentation Grant. We acknowledge the assistance of the Danish National High-Throughput Sequencing Centre, BGI-Europe, the Garvan Institute of Medical Research and the Australian Cancer Research Foundation (ACRF) Cancer Genomics Facility for assistance in Illumina and BGIseq500 data generation. A.R.P. was supported by a Marie Curie COFUND Junior Research Fellowship (Durham University). A.M. was supported by an NSF grant (award number: 1457106) and the QCB Collaboratory Postdoctoral Fellowship (UCLA). L.A.F.F., J.H., A.H.-B. and G.L. were supported by either European Research Council grant (ERC-2013-StG-337574-UNDEAD and ERC-2019-StG-853272-PALAEOFARM) and/or Natural Environmental Research Council grants (NE/K005243/1 and NE/K003259/1). K.S. was supported by a grant from Barrett, the Honors College at Arizona State University. A.T.O. was supported by the Strategic Initiative Funds, Office of the President, Arizona State University to the Institute of Human Origins DNA and Human Origins at Arizona State University project. L.A.F.F. was supported by a Junior Research Fellowship (Wolfson College, University of Oxford) and L.A.F.F. and A. Carmagnini were supported by the Wellcome Trust (210119/Z/18/Z). S.G. was supported by Carlsbergfondet grant CF14–0995 and Marie Skłodowska-Curie Actions grant 655732-WhereWolf. M.T.P.G. was supported by ERC Consolidator grant 681396-Extinction Genomics. B.S. and J.K. were supported by IMLS MG-30-17-0045-17 and NSF DEB-1754451. A.H.-B. was supported by the Leverhulme Trust (ECF-2017-315). A. Cooper, K.J.M. and H.H. were supported by the Australian Research Council. A.T.S. and G.G. were supported by Australian Government Research Training Program Scholarships. A.T.L. was supported by the Peter Buck Postdoctoral Fellowship from the Smithsonian Institution’s National Museum of Natural History. Y.V.K. was supported by the by State Assignment of the Sobolev Institute of Geology and Mineralogy.

The resolution of conservation genetic analyses has been limited until recently due to technological and computational challenges associated with genotyping multiple loci at once. In this review, we focus on how the development of high-throughput genotyping methods have enabled conservation genomics studies of wolves in North America. The gray wolf (Canis lupus) historically had a Holarctic distribution across widely varying environments, yet during the early twentieth century many populations declined due to direct persecution and other anthropogenic disturbances. First, we discuss genetic substructure and adaptive uniqueness among genetically and environmentally defined wolf ecotypes. Second, we focus on the new conservation implications revealed by studies having increased genomic resolution of the dynamics of reintroduced and re-established wolves, specifically Mexican and Pacific Northwest wolves. Mexican wolves, a distinct subspecies of North American wolf that inhabit a small area within the southwestern U.S. and Mexico, remain endangered despite decades since a reintroduction program began. How biologists and management agencies use scientific data to define the historical range of Mexican wolves will be critical to future reintroduction efforts. In the Pacific Northwest, admixture occurs between the distinct and declining coastal wolf ecotype and the more abundant reintroduced interior wolves. If coastal wolves obtain protection, then the Pacific Northwest wolves may also warrant protection. Therefore, more precise policies are needed for the management of admixed populations when one source is protected. We recommend that future conservation efforts should provide full protection for distinct ecotypes, support scientifically rigorous definitions of historical range to inform restoration, and enhance the legal status of admixed populations.

Admixture resulting from natural dispersal processes can potentially generate novel phenotypic variation that may facilitate persistence in changing environments or result in the loss of population-specific adaptations. Yet, under the US Endangered Species Act, policy is limited for management of individuals whose ancestry includes a protected taxon; therefore, they are generally not protected under the Act. This issue is exemplified by the recently re-established grey wolves of the Pacific Northwest states of Washington and Oregon, USA. This population was likely founded by two phenotypically and genetically distinct wolf ecotypes: Northern Rocky Mountain (NRM) forest and coastal rainforest. The latter is considered potentially threatened in southeast Alaska and thus the source of migrants may affect plans for their protection. To assess the genetic source of the re-established population, we sequenced a ~ 300 bp portion of the mitochondrial control region and ~ 5 Mbp of the nuclear genome. Genetic analysis revealed that the Washington wolves share ancestry with both wolf ecotypes, whereas the Oregon population shares ancestry with NRM forest wolves only. Using ecological niche modelling, we found that the Pacific Northwest states contain environments suitable for each ecotype, with wolf packs established in both environmental types. Continued migration from coastal rainforest and NRM forest source populations may increase the genetic diversity of the Pacific Northwest population. However, this admixed population challenges traditional management regimes given that admixture occurs between an adaptively distinct ecotype and a more abundant reintroduced interior form. Our results emphasize the need for a more precise US policy to address the general problem of admixture in the management of endangered species, subspecies, and distinct population segments.

Skin pigmentation and coat pigmentation are two of the best-studied examples of traits under natural selection given their quantifiable fitness interactions with the environment (e.g., camouflage) and signalling with other organisms (e.g., warning coloration). Previous morphological studies have found that skin pigmentation variation in the Virginia opossum (Didelphis virginiana) is associated with variation in precipitation and temperatures across its distribution range following Gloger's rule (lighter pigmentation in temperate environments). To investigate the molecular mechanism associated with skin pigmentation variation, we used RNA-Seq and quantified gene expression of wild opossums from tropical and temperate populations. Using differential expression analysis and a co-expression network approach, we found that expression variation in genes with melanocytic and immune functions is significantly associated with the degree of skin pigmentation variation and may be underlying this phenotypic difference. We also found evidence suggesting that the Wnt/β-catenin signalling pathway might be regulating the depigmentation observed in temperate populations. Based on our study results, we present several alternative hypotheses that may explain Gloger's rule pattern of skin pigmentation variation in opossum, including changes in pathogen diversity supporting a pathogen-resistant hypothesis, thermal stress associated with temperate environments, and pleiotropic and epistatic interactions between melanocytic and immune genes.