Department of Ecology and Evolutionary Biology

The Yellow Warbler (Setophaga petechia) is a small songbird in the wood-warbler family (Parulidae) that exhibits phenotypic and ecological differences across a widespread distribution and is important to California's riparian habitat conservation. Here, we present a high-quality de novo genome assembly of a vouchered female Yellow Warbler from southern California. Using HiFi long-read and Omni-C proximity sequencing technologies, we generated a 1.22 Gb assembly including 687 scaffolds with a contig N50 of 6.80 Mb, scaffold N50 of 21.18 Mb, and a BUSCO completeness score of 96.0%. This highly contiguous genome assembly provides an essential resource for understanding the history of gene flow, divergence, and local adaptation in Yellow Warblers and can inform conservation management of this charismatic bird species.

Adaptation to extreme environments often involves the evolution of dramatic physiological changes. To better understand how organisms evolve these complex phenotypic changes, the repeatability and predictability of evolution, and possible constraints on adapting to an extreme environment, it is important to understand how adaptive variation has evolved. Poeciliid fishes represent a particularly fruitful study system for investigations of adaptation to extreme environments due to their repeated colonization of toxic hydrogen sulfide-rich springs across multiple species within the clade. Previous investigations have highlighted changes in the physiology and gene expression in specific species that are thought to facilitate adaptation to hydrogen sulfide-rich springs. However, the presence of adaptive nucleotide variation in coding and regulatory regions and the degree to which convergent evolution has shaped the genomic regions underpinning sulfide tolerance across taxa are unknown. By sampling across seven independent lineages in which nonsulfidic lineages have colonized and adapted to sulfide springs, we reveal signatures of shared evolutionary rate shifts across the genome. We found evidence of genes, promoters, and putative enhancer regions associated with both increased and decreased convergent evolutionary rate shifts in hydrogen sulfide-adapted lineages. Our analysis highlights convergent evolutionary rate shifts in sulfidic lineages associated with the modulation of endogenous hydrogen sulfide production and hydrogen sulfide detoxification. We also found that regions with shifted evolutionary rates in sulfide spring fishes more often exhibited convergent shifts in either the coding region or the regulatory sequence of a given gene, rather than both.

Genetic diversity is a fundamental component of biodiversity. Examination of global patterns of genetic diversity can help highlight mechanisms underlying species diversity, though a recurring challenge has been that patterns may vary by molecular marker. Here, we compiled 6862 observations of genetic diversity from 492 species of marine fish and tested among hypotheses for diversity gradients: the founder effect hypothesis, the kinetic energy hypothesis, and the productivity-diversity hypothesis. We fit generalized linear mixed effect models (GLMMs) and explored the extent to which various macroecological drivers (latitude, longitude, temperature (SST), and chlorophyll-a concentration) explained variation in genetic diversity. We found that mitochondrial genetic diversity followed geographic gradients similar to those of species diversity, being highest near the Equator, particularly in the Coral Triangle, while nuclear genetic diversity did not follow clear geographic patterns. Despite these differences, all genetic diversity metrics were correlated with chlorophyll-a concentration, while mitochondrial diversity was also positively associated with SST. Our results provide support for the kinetic energy hypothesis, which predicts that elevated mutation rates at higher temperatures increase mitochondrial but not necessarily nuclear diversity, and the productivity-diversity hypothesis, which posits that resource-rich regions support larger populations with greater genetic diversity. Overall, these findings reveal how environmental variables can influence mutation rates and genetic drift in the ocean, caution against using mitochondrial macrogenetic patterns as proxies for whole-genome diversity, and aid in defining global gradients of genetic diversity.

Many animals and plants have species-typical annual cycles, but individuals vary in their timing of life-history events. Individual variation in fur replacement (moult) timing is poorly understood in mammals due to the challenge of repeated observations and longitudinal sampling. We examined factors that influence variation in moult duration and timing among elephant seals (Mirounga angustirostris). We quantified the onset and progression of fur loss in 1178 individuals. We found that an exceptionally rapid visible moult (7 days, the shortest of any mammals or birds), and a wide range of moult start dates (spanning 6-10× the event duration) facilitated high asynchrony across individuals (only 20% of individuals in the population moulting at the same time). Some of the variation was due to reproductive state, as reproductively mature females that skipped a breeding season moulted a week earlier than reproductive females. Moreover, individual variation in timing and duration within age-sex categories far outweighed (76-80%) variation among age-sex categories. Individuals arriving at the end of the moult season spent 50% less time on the beach, which allowed them to catch up in their annual cycles and reduce population-level variance during breeding. These findings underscore the importance of individual variation in annual cycles.

Reference genome assemblies have been created from multiple lineages within the Canidae family; however, despite its phylogenetic relevance as a basal genus within the clade, there is currently no reference genome for the gray fox (Urocyon cinereoargenteus). Here, we present a chromosome-level assembly for the gray fox (U. cinereoargenteus), which represents the most contiguous, non-domestic canid reference genome available to date, with 90% of the genome contained in just 34 scaffolds and a contig N50 and scaffold N50 of 59.4 and 72.9 Megabases, respectively. Repeat analyses identified an increased number of simple repeats relative to other canids. Based on mitochondrial DNA, our Vermont sample clusters with other gray fox samples from the northeastern United States and contains slightly lower levels of heterozygosity than gray foxes on the west coast of California. This new assembly lays the groundwork for future studies to describe past and present population dynamics, including the delineation of evolutionarily significant units of management relevance. Importantly, the phylogenetic position of Urocyon allows us to verify the loss of PRDM9 functionality in the basal canid lineage, confirming that pseudogenization occurred at least 10 million years ago.

Pangenome references address biases of reference genomes by storing a representative set of diverse haplotypes and their alignment, usually as a graph. Alternate alleles determined by variant callers can be used to construct pangenome graphs, but advances in long-read sequencing are leading to widely available, high-quality phased assemblies. Constructing a pangenome graph directly from assemblies, as opposed to variant calls, leverages the graph's ability to represent variation at different scales. Here we present the Minigraph-Cactus pangenome pipeline, which creates pangenomes directly from whole-genome alignments, and demonstrate its ability to scale to 90 human haplotypes from the Human Pangenome Reference Consortium. The method builds graphs containing all forms of genetic variation while allowing use of current mapping and genotyping tools. We measure the effect of the quality and completeness of reference genomes used for analysis within the pangenomes and show that using the CHM13 reference from the Telomere-to-Telomere Consortium improves the accuracy of our methods. We also demonstrate construction of a Drosophila melanogaster pangenome.

Ecological and evolutionary theories have proposed that species traits should be important in mediating species responses to contemporary climate change; yet, empirical evidence has so far provided mixed evidence for the role of behavioral, life history, or ecological characteristics in facilitating or hindering species range shifts. As such, the utility of trait-based approaches to predict species redistribution under climate change has been called into question. We develop the perspective, supported by evidence, that trait variation, if used carefully can have high potential utility, but that past analyses have in many cases failed to identify an explanatory value for traits by not fully embracing the complexity of species range shifts. First, we discuss the relevant theory linking species traits to range shift processes at the leading (expansion) and trailing (contraction) edges of species distributions and highlight the need to clarify the mechanistic basis of trait-based approaches. Second, we provide a brief overview of range shift-trait studies and identify new opportunities for trait integration that consider range-specific processes and intraspecific variability. Third, we explore the circumstances under which environmental and biotic context dependencies are likely to affect our ability to identify the contribution of species traits to range shift processes. Finally, we propose that revealing the role of traits in shaping species redistribution may likely require accounting for methodological variation arising from the range shift estimation process as well as addressing existing functional, geographical, and phylogenetic biases. We provide a series of considerations for more effectively integrating traits as well as extrinsic and methodological factors into species redistribution research. Together, these analytical approaches promise stronger mechanistic and predictive understanding that can help society mitigate and adapt to the effects of climate change on biodiversity.

Populations and species are threatened by human pressure, but their fate is variable. Some depleted populations, such as that of the northern elephant seal (Mirounga angustirostris), recover rapidly even when the surviving population was small. The northern elephant seal was hunted extensively and taken by collectors between the early 1800s and 1892, suffering an extreme population bottleneck as a consequence. Recovery was rapid and now there are over 200,000 individuals. We sequenced 260 modern and 8 historical northern elephant seal nuclear genomes to assess the impact of the population bottleneck on individual northern elephant seals and to better understand their recovery. Here we show that inbreeding, an increase in the frequency of alleles compromised by lost function, and allele frequency distortion, reduced the fitness of breeding males and females, as well as the performance of adult females on foraging migrations. We provide a detailed investigation of the impact of a severe bottleneck on fitness at the genomic level and report on the role of specific gene systems.

Range expansions, whether they are biological invasions or climate change-mediated range shifts, may have profound ecological and evolutionary consequences for plant-soil interactions. Range-expanding plants encounter soil biota with which they have a limited coevolutionary history, especially when introduced to a new continent. Past studies have found mixed results on whether plants experience positive or negative soil feedback interactions in their novel range, and these effects often change over time. One important theoretical explanation is that plants locally adapt to the soil pathogens and mutualists in their novel range. We tested this hypothesis in Dittrichia graveolens, an annual plant that is both expanding its European native range, initially coinciding with climate warming, and rapidly invading in California after human introduction. In parallel greenhouse experiments on both continents, we used plant genotypes and soils from five locations at the core and edge of each range to compare plant growth in soil inhabited by D. graveolens and nearby control microsites as a measure of plant-soil feedback. Plant-soil interactions were highly idiosyncratic across each range. On average, plant-soil feedbacks were more positive in the native range than in the exotic range. In line with the strongly heterogeneous pattern of soil responses along our biogeographic gradients, we found no evidence for evolutionary differentiation between plant genotypes from the core to edge of either range. Our results suggest that the evolution of plant-soil interactions during range expansion may be more strongly driven by local evolutionary dynamics varying across the range than by large-scale biogeographic shifts.

Substantial global attention is focused on how to reduce the risk of future pandemics. Reducing this risk requires investment in prevention, preparedness, and response. Although preparedness and response have received significant focus, prevention, especially the prevention of zoonotic spillover, remains largely absent from global conversations. This oversight is due in part to the lack of a clear definition of prevention and lack of guidance on how to achieve it. To address this gap, we elucidate the mechanisms linking environmental change and zoonotic spillover using spillover of viruses from bats as a case study. We identify ecological interventions that can disrupt these spillover mechanisms and propose policy frameworks for their implementation. Recognizing that pandemics originate in ecological systems, we advocate for integrating ecological approaches alongside biomedical approaches in a comprehensive and balanced pandemic prevention strategy.