Search

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.

Migration is driven by a combination of environmental and genetic factors, but many questions remain about those drivers. Potential interactions between genetic and environmental variants associated with different migratory phenotypes are rarely the focus of study. We pair low coverage whole genome resequencing with a de novo genome assembly to examine population structure, inbreeding, and the environmental factors associated with genetic differentiation between migratory and resident breeding phenotypes in a species of conservation concern, the western burrowing owl (Athene cunicularia hypugaea). Our analyses reveal a dichotomy in gene flow depending on whether the population is resident or migratory, with the former being genetically structured and the latter exhibiting no signs of structure. Among resident populations, we observed significantly higher genetic differentiation, significant isolation-by-distance, and significantly elevated inbreeding. Among migratory breeding groups, on the other hand, we observed lower genetic differentiation, no isolation-by-distance, and substantially lower inbreeding. Using genotype-environment association analysis, we find significant evidence for relationships between migratory phenotypes (i.e., migrant versus resident) and environmental variation associated with cold temperatures during the winter and barren, open habitats. In the regions of the genome most differentiated between migrants and residents, we find significant enrichment for genes associated with the metabolism of fats. This may be linked to the increased pressure on migrants to process and store fats more efficiently in preparation for and during migration. Our results provide a significant contribution toward understanding the evolution of migratory behavior and vital insight into ongoing conservation and management efforts for the western burrowing owl.

The extent of parallel genomic responses to similar selective pressures depends on a complex array of environmental, demographic, and evolutionary forces. Laboratory experiments with replicated selective pressures yield mixed outcomes under controlled conditions and our understanding of genomic parallelism in the wild is limited to a few well-established systems. Here, we examine genomic signals of selection in the eelgrass Zostera marina across temperature gradients in adjacent embayments. Although we find many genomic regions with signals of selection within each bay there is very little overlap in signals of selection at the SNP level, despite most polymorphisms being shared across bays. We do find overlap at the gene level, potentially suggesting multiple mutational pathways to the same phenotype. Using polygenic models we find that some sets of candidate SNPs are able to predict temperature across both bays, suggesting that small but parallel shifts in allele frequencies may be missed by independent genome scans. Together, these results highlight the continuous rather than binary nature of parallel evolution in polygenic traits and the complexity of evolutionary predictability.

Identifying areas of high evolutionary potential is a judicious strategy for developing conservation priorities in the face of environmental change. For wide-ranging species occupying heterogeneous environments, the evolutionary forces that shape distinct populations can vary spatially. Here, we investigate patterns of genomic variation and genotype-environment associations in the hermit thrush (Catharus guttatus), a North American songbird, at broad (across the breeding range) and narrow spatial scales (at a hybrid zone). We begin by building a genoscape or map of genetic variation across the breeding range and find five distinct genetic clusters within the species, with the greatest variation occurring in the western portion of the range. Genotype-environment association analyses indicate higher allelic turnover in the west than in the east, with measures of temperature surfacing as key predictors of putative adaptive genomic variation rangewide. Since broad patterns detected across a species' range represent the aggregate of many locally adapted populations, we investigate whether our broadscale analysis is consistent with a finer scale analysis. We find that top rangewide temperature-associated loci vary in their clinal patterns (e.g., steep clines vs. fixed allele frequencies) across a hybrid zone in British Columbia, suggesting that the environmental predictors and the associated candidate loci identified in the rangewide analysis are of variable importance in this particular region. However, two candidate loci exhibit strong concordance with the temperature gradient in British Columbia, suggesting a potential role for temperature-related barriers to gene flow and/or temperature-driven ecological selection in maintaining putative local adaptation. This study demonstrates how patterns identified at the broad (macrogeographic) scale can be validated by investigating genotype-environment correlations at the local (microgeographic) scale. Furthermore, our results highlight the importance of considering the spatial distribution of putative adaptive variation when assessing population-level sensitivity to climate change and other stressors.

Reef-building coral species are experiencing an unprecedented decline owing to increasing frequency and intensity of marine heatwaves and associated bleaching-induced mortality. Closely related species from the Acropora hyacinthus species complex differ in heat tolerance and in their association with heat-tolerant symbionts. We used low-coverage full genome sequencing of 114 colonies monitored across the 2015 bleaching event in American Samoa to determine the genetic differences among four cryptic species (termed HA, HC, HD and HE) that have diverged in these species traits. Cryptic species differed strongly at thousands of single nucleotide polymorphisms across the genome which are enriched for amino acid changes in the bleaching-resistant species HE. In addition, HE also showed two particularly divergent regions with strong signals of differentiation. One approximately 220 kb locus, HES1, contained the majority of fixed differences in HE. A second locus, HES2, was fixed in HE but polymorphic in the other cryptic species. Surprisingly, non-HE individuals with HE-like haplotypes at HES2 were more likely to bleach. At both loci, HE showed particular sequence similarity to a congener, Acropora millepora. Overall, resilience to bleaching during the third global bleaching event was strongly structured by host cryptic species, buoyed by differences in symbiont associations between these species.

Habitat loss and alteration has driven many species into decline, often to the point of requiring protection and intervention to avert extinction. Genomic data provide the opportunity to inform conservation and recovery efforts with details about vital evolutionary processes with a resolution far beyond that of traditional genetic approaches. The tricolored blackbird (Agelaius tricolor) has suffered severe losses during the previous century largely due to anthropogenic impacts on their habitat. Using a dataset composed of a whole genome paired with reduced representation libraries (RAD-Seq) from samples collected across the species' range, we find evidence for panmixia using multiple methods, including PCA (no geographic clustering), admixture analyses (ADMIXTURE and TESS conclude K = 1), and comparisons of genetic differentiation (average FST = 0.029). Demographic modeling approaches recovered an ancient decline that had a strong impact on genetic diversity but did not detect any effect from the known recent decline. We also did not detect any evidence for selection, and hence adaptive variation, at any site, either geographic or genomic. These results indicate that species continues to have high vagility across its range despite population decline and habitat loss and should be managed as a single unit.

Recent research suggests that bird populations are declining at alarming rates across the United States. Over the last century, local, national, and international efforts to limit declines in bird populations have resulted in state and federal laws that now protect most of California’s 650 bird species. To comply with these protections, transportation infrastructure projects often face strict mitigation requirements, which can be expensive and cause delays. Understanding a project’s effects on specific bird populations can refine mitigation requirements and optimize infrastructure planning.The tools for identifying separate populations within species and understanding protected birds’ seasonal movement have been limited until recently. New genetic and genomic tools now provide a method for understanding population differentiation, which is vital to a wide array of conservation goals including estimating population declines, identifying potential for adaptation to stressors, measuring connectivity between populations, and estimating inbreeding. Because infrastructure projects can directly impact genetic diversity and connectivity, a toolkit to assess population structure and the distribution of genetic variation could aid in predicting and mitigating the impacts of such projects. As a test case, researchers at the University of California, Davis, sequenced entire genomes for 40 individual Anna’s hummingbirds (Calypte anna) from across California to identify breeding populations and develop a genetic toolkit to assign individuals to those populations. The presence of this species at bridge construction sites has resulted in construction delays in part because little information exists on the status of different populations within the species.This research brief summarizes findings and implications from the project.View the NCST Project Webpage

Recent studies have highlighted alarming rates of declines in bird populations across the country. The State of California is home to over 650 resident and migrant avian species. Legislation for protecting these species has existed for over a century now, yet tools for identifying populations and understanding seasonal movement remain limited. Recently, genetic and genomic tools have provided a method for understanding population structure, allowing for more informed delineation of management units. The goal of this project was to create a genetic toolkit for identifying breeding populations and assigning individuals to those populations. Ultimately, such tools could be used to assess population-level impacts when there are conflicts with birds at infrastructure construction sites. As a test case, we sequenced entire genomes for 40 individual Anna’s hummingbirds (Calypte anna) from across the state. Based on this initial data, we found low levels of differentiation between sampled locations, suggesting that C. anna in California are not subdivided into different population units. However, there was a weak signal of geography suggesting there may be localized genetic differences in a small proportion of the genome. Follow-up work will focus on a broader sampling across the state of California to clarify any possible population subdivision or geographical patterns of differentiation.View the NCST Project Webpage

Despite the enormous advances in genetics, links between phenotypes and genotypes have been made for only a few nonmodel organisms. However, such links can be essential to understand mechanisms of ecological speciation. The Costa Rican endemic Mangrove Warbler subspecies provides an excellent subject to study differentiation with gene flow, as it is distributed along a strong precipitation gradient on the Pacific coast with no strong geographic barriers to isolate populations. Mangrove Warbler populations could be subject to divergent selection driven by precipitation, which influences soil salinity levels, which in turn influences forest structure and food resources. We used single nucleotide polymorphisms (SNPs) and morphological traits to examine the balance between neutral genetic and phenotypic divergence to determine whether selection has acted on traits and genes with functions related to specific environmental variables. We present evidence showing: (a) associations between environmental variables and SNPs, identifying candidate genes related to bill morphology (BMP) and osmoregulation, (b) absence of population genetic structure in neutrally evolving markers, (c) divergence in bill size across the precipitation gradient, and (d) strong phenotypic differentiation (P ST) which largely exceeds neutral genetic differentiation (F ST) in bill size. Our results indicate an important role for salinity, forest structure, and resource availability in maintaining phenotypic divergence of Mangrove Warblers through natural selection. Our findings add to the growing body of literature identifying the processes involved in phenotypic differentiation along environmental gradients in the face of gene flow.

Population genomic surveys suggest that climate-associated genetic variation occurs widely across species, but whether it is sufficient to allow population persistence via evolutionary adaptation has seldom been quantified. To ask whether rapid adaptation in reef-building corals can keep pace with future ocean warming, we measured genetic variation at predicted warm-adapted loci and simulated future evolution and persistence in a high-latitude population of corals from Rarotonga, Cook Islands. Alleles associated with thermal tolerance were present but at low frequencies in this cooler, southerly locality. Simulations based on predicted ocean warming in Rarotonga showed rapid evolution of heat tolerance resulting in population persistence under mild warming scenarios consistent with low CO2 emission plans, RCP2.6 and RCP4.5. Under more severe scenarios, RCP6.0 and RCP8.5, adaptation was not rapid enough to prevent extinction. Population adaptation was faster for models based on smaller numbers of additive loci that determine thermal tolerance and for higher population growth rates. Finally, accelerated migration via transplantation of thermally tolerant individuals (1 to 5%/year) sped adaptation. These results show that cool-water corals can adapt to warmer oceans but only under mild scenarios resulting from international emissions controls. Incorporation of genomic data into models of species response to climate change offers a promising method for estimating future adaptive processes.

Ecological speciation with gene flow is widespread in nature [1], but it presents a conundrum: how are associations between traits under divergent natural selection and traits that contribute to assortative mating maintained? Theoretical models suggest that genetic mechanisms inhibiting free recombination between loci underlying these two types of traits (hereafter, "genetic coupling") can facilitate speciation [2-4]. Here, we perform a direct test for genetic coupling by mapping both divergent traits and female mate choice in a classic model of ecological speciation: sympatric benthic and limnetic threespine stickleback (Gasterosteus aculeatus). By measuring mate choice in F2 hybrid females, we allowed for recombination between loci underlying assortative mating and those under divergent ecological selection. In semi-natural mating arenas in which females had access to both benthic and limnetic males, we found that F2 females mated with males similar to themselves in body size and shape. In addition, we found two quantitative trait loci (QTLs) associated with female mate choice that also predicted female morphology along the benthic-limnetic trait axis. Furthermore, a polygenic genetic model that explains adaptation to contrasting benthic and limnetic feeding niches [5] also predicted F2 female mate choice. Together, these results provide empirical evidence that genetic coupling of assortative mating with traits under divergent ecological selection helps maintain species in the face of gene flow, despite a polygenic basis for adaptation to divergent environments.

Concern over rapid environmental shifts associated with climate change has led to a search for molecular markers of environmental tolerance. Climate-associated gene expression profiles exist for a number of systems, but have rarely been tied to fitness outcomes, especially in nonmodel organisms. We reciprocally transplanted corals between two backreef locations with more and less variable temperature regimes to disentangle effects of recent and native environment on survival and growth. Coral growth over 12 months was largely determined by local environment. Survival, however, was impacted by native environment; corals from the more variable environment had 22% higher survivorship. By contrast, corals native to the less variable environment had more variable survival. This might represent a "selective sieve" where poor survivors are filtered from the more stressful environment. We also find a potential fitness trade-off-corals with high survival under stressful conditions grew less in the more benign environment. Transcriptome samples taken a year before transplantation were used to examine gene expression patterns that predicted transplant survival and growth. Two separate clusters of coexpressed genes were predictive of survival in the two locations. Genes from these clusters are candidate biomarkers for predicting persistence of corals under future climate change scenarios.

In animals, introgression between species is often perceived as the breakdown of reproductive isolating mechanisms, but gene flow between incipient species can also represent a source for potentially beneficial alleles. Recently, genome-wide datasets have revealed clusters of differentiated loci ('genomic islands of divergence') that are thought to play a role in reproductive isolation and therefore have reduced gene flow. We use simulations to further examine the evolutionary forces that shape and maintain genomic islands of divergence between two subspecies of the migratory songbird, Swainson's thrush (Catharus ustulatus), which have come into secondary contact since the last glacial maximum. We find that, contrary to expectation, gene flow is high within islands and is highly asymmetric. In addition, patterns of nucleotide diversity at highly differentiated loci suggest selection was more frequent in a single ecotype. We propose a mechanism whereby beneficial alleles spread via selective sweeps following a post-glacial demographic expansion in one subspecies and move preferentially across the hybrid zone. We find no evidence that genomic islands are the result of divergent selection or reproductive isolation, rather our results suggest that differentiated loci both within and outside islands could provide opportunities for adaptive introgression across porous species boundaries.

Population response to environmental variation involves adaptation, acclimation, or both. For long-lived organisms, acclimation likely generates a faster response but is only effective if the rates and limits of acclimation match the dynamics of local environmental variation. In coral reef habitats, heat stress from extreme ocean warming can occur over several weeks, resulting in symbiont expulsion and widespread coral death. However, transcriptome regulation during short-term acclimation is not well understood. We examined acclimation during a 11-day experiment in the coral Acropora nana. We acclimated colonies to three regimes: ambient temperature (29 °C), increased stable temperature (31 °C), and variable temperature (29-33 °C), mimicking local heat stress conditions. Within 7-11 days, individuals acclimated to increased temperatures had higher tolerance to acute heat stress. Despite physiological changes, no gene expression changes occurred during acclimation before acute heat stress. However, we found strikingly different transcriptional responses to heat stress between acclimation treatments across 893 contigs. Across these contigs, corals acclimated to higher temperatures (31 °C or 29-33 °C) exhibited a muted stress response--the magnitude of expression change before and after heat stress was less than in 29 °C acclimated corals. Our results show that corals have a rapid phase of acclimation that substantially increases their heat resilience within 7 days and that alters their transcriptional response to heat stress. This is in addition to a previously observed longer term response, distinguishable by its shift in baseline expression, under nonstressful conditions. Such rapid acclimation may provide some protection for this species of coral against slow onset of warming ocean temperatures.

Reef-building coral species are experiencing an unprecedented decline owing to increasing frequency and intensity of marine heatwaves and associated bleaching-induced mortality. Closely related species from the Acropora hyacinthus species complex differ in heat tolerance and in their association with heat-tolerant symbionts. We used low-coverage full genome sequencing of 114 colonies monitored across the 2015 bleaching event in American Samoa to determine the genetic differences among four cryptic species (termed HA, HC, HD and HE) that have diverged in these species traits. Cryptic species differed strongly at thousands of single nucleotide polymorphisms across the genome which are enriched for amino acid changes in the bleaching-resistant species HE. In addition, HE also showed two particularly divergent regions with strong signals of differentiation. One approximately 220 kb locus, HES1, contained the majority of fixed differences in HE. A second locus, HES2, was fixed in HE but polymorphic in the other cryptic species. Surprisingly, non-HE individuals with HE-like haplotypes at HES2 were more likely to bleach. At both loci, HE showed particular sequence similarity to a congener, Acropora millepora . Overall, resilience to bleaching during the third global bleaching event was strongly structured by host cryptic species, buoyed by differences in symbiont associations between these species.