Genomics-informed conservation units reveal spatial variation in climate vulnerability in a migratory bird

Identifying genetic conservation units (CUs) in threatened species is critical for the preservation of adaptive capacity and evolutionary potential in the face of climate change. However, delineating CUs in highly mobile species remains a challenge due to high rates of gene flow and genetic signatures of isolation by distance. Even when CUs are delineated in highly mobile species, the CUs often lack key biological information about what populations have the most conservation need to guide management decisions. Here we implement a framework for CU identification in the Canada Warbler (Cardellina canadensis), a migratory bird species of conservation concern, and then integrate demographic modelling and genomic offset to guide conservation decisions. We find that patterns of whole genome genetic variation in this highly mobile species are primarily driven by putative adaptive variation. Identification of CUs across the breeding range revealed that Canada Warblers fall into two evolutionarily significant units (ESU), and three putative adaptive units (AUs) in the South, East, and Northwest. Quantification of genomic offset, a metric of genetic changes necessary to maintain current gene-environment relationships, revealed significant spatial variation in climate vulnerability, with the Northwestern AU being identified as the most vulnerable to future climate change. Alternatively, quantification of past population trends within each AU revealed the steepest population declines have occurred within the Eastern AU. Overall, we illustrate that genomics-informed CUs provide a strong foundation for identifying current and future regional threats that can be used to inform management strategies for a highly mobile species in a rapidly changing world.

Genetic and environmental drivers of migratory behavior in western burrowing owls and implications for conservation and management

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.

Genomic approaches to mitigating genetic diversity loss in declining populations

The accelerating pace of global biodiversity loss is exacerbated by habitat fragmentation and subsequent inbreeding in small populations. To address this problem, conservation practitioners often turn to assisted breeding programmes with the aim of enhancing genetic diversity in declining populations. Although genomic information is infrequently included in these efforts, it has the potential to significantly enhance the success of such programmes. In this study, we showcase the value of genomic approaches for increasing genetic diversity in assisted breeding efforts, specifically focusing on a highly inbred population of Western burrowing owls. To maximize genetic diversity in the resulting offspring, we begin by creating an optimal pairing decision tree based on sex, kinship and patterns of homozygosity across the genome. To evaluate the effectiveness of our strategy, we compare genetic diversity, brood size and nestling success rates between optimized and non-optimized pairs. Additionally, we leverage recently discovered correlations between telomere length and fitness across species to investigate whether genomic optimization could have long-term fitness benefits. Our results indicate that pairing individuals with contrasting patterns of homozygosity across the genome is an effective way to increase genetic diversity in offspring. Although short-term field-based metrics of success did not differ significantly between optimized and non-optimized pairs, offspring from optimized pairs had significantly longer telomeres, suggesting that genetic optimization can help reduce the risk of inbreeding depression. These findings underscore the importance of genomic tools for informing efforts to preserve the adaptive potential of small, inbred populations at risk of further decline.

Low-coverage whole genome sequencing for highly accurate population assignment: Mapping migratory connectivity in the American Redstart (Setophaga ruticilla)

Understanding the geographic linkages among populations across the annual cycle is an essential component for understanding the ecology and evolution of migratory species and for facilitating their effective conservation. While genetic markers have been widely applied to describe migratory connections, the rapid development of new sequencing methods, such as low-coverage whole genome sequencing (lcWGS), provides new opportunities for improved estimates of migratory connectivity. Here, we use lcWGS to identify fine-scale population structure in a widespread songbird, the American Redstart (Setophaga ruticilla), and accurately assign individuals to genetically distinct breeding populations. Assignment of individuals from the nonbreeding range reveals population-specific patterns of varying migratory connectivity. By combining migratory connectivity results with demographic analysis of population abundance and trends, we consider full annual cycle conservation strategies for preserving numbers of individuals and genetic diversity. Notably, we highlight the importance of the Northern Temperate-Greater Antilles migratory population as containing the largest proportion of individuals in the species. Finally, we highlight valuable considerations for other population assignment studies aimed at using lcWGS. Our results have broad implications for improving our understanding of the ecology and evolution of migratory species through conservation genomics approaches.

Where to draw the line? Expanding the delineation of conservation units to highly mobile taxa

Conservation units (CUs) are an essential tool for maximizing evolutionary potential and prioritizing areas across a species' range for protection when implementing conservation and management measures. However, current workflows for identifying CUs on the basis of neutral and adaptive genomic variation largely ignore information contained in patterns of isolation by distance (IBD), frequently the primary signal of population structure in highly mobile taxa, such as birds, bats, and marine organisms with pelagic larval stages. While individuals located on either end of a species' distribution may exhibit clear genetic, phenotypic, and ecological differences, IBD produces subtle changes in allele frequencies across space, making it difficult to draw clear boundaries for conservation purposes in the absence of discrete population structure. Here, we highlight potential pitfalls that arise when applying common methods for delineating CUs to continuously distributed organisms and review existing methods for detecting subtle breakpoints in patterns of IBD that can indicate barriers to gene flow in highly mobile taxa. In addition, we propose a new framework for identifying CUs in all organisms, including those characterized by continuous genomic differentiation, and suggest several possible ways to harness the information contained in patterns of IBD to guide conservation and management decisions.

The genetic basis of plumage coloration and elevation adaptation in a clade of recently diverged alpine and arctic songbirds

Trait genetic architecture plays an important role in the probability that variation in that trait leads to divergence and speciation. In some cases, speciation may be driven by the generation of novel phenotypes through the recombination of genes associated with traits that are important for local adaptation or sexual selection. Here, we investigate the genetic basis of three plumage color traits, and one ecological trait, breeding elevation, in a recent avian radiation, the North American rosy-finches (Leucosticte spp.). We identify unique genomic regions associated with each trait and highlight 11 candidate genes. Among these are well-characterized melanogenesis genes, including Mitf and Tyrp1, and previously reported hypoxia-related genes including Egln1. Additionally, we use mitochondrial data to date the divergence of rosy-finch clades which appear to have diverged within the past 250 ky. Given the low levels of genome-wide differentiation among rosy-finch taxa, and evidence for extensive introgression in North America, plumage coloration and adaptation to high elevations have likely played large roles in generating the observed patterns of lineage divergence. The relative independence of these candidate regions across the genome suggests that recombination might have led to multiple phenotypes, and subsequent rosy-finch speciation, over short periods of time.

Genotype‐environment associations across spatial scales reveal the importance of putative adaptive genetic variation in divergence

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.

Clock-linked genes underlie seasonal migratory timing in a diurnal raptor

Seasonal migration is a dynamic natural phenomenon that allows organisms to exploit favourable habitats across the annual cycle. While the morphological, physiological and behavioural changes associated with migratory behaviour are well characterized, the genetic basis of migration and its link to endogenous biological time-keeping pathways are poorly understood. Historically, genome-wide research has focused on genes of large effect, whereas many genes of small effect may work together to regulate complex traits like migratory behaviour. Here, we explicitly relax stringent outlier detection thresholds and, as a result, discover how multiple biological time-keeping genes are important to migratory timing in an iconic raptor species, the American kestrel (Falco sparverius). To validate the role of candidate loci in migratory timing, we genotyped kestrels captured across autumn migration and found significant associations between migratory timing and genetic variation in metabolic and light-input pathway genes that modulate biological clocks (top1, phlpp1, cpne4 and peak1). Further, we demonstrate that migrating individuals originated from a single panmictic source population, suggesting the existence of distinct early and late migratory genotypes (i.e. chronotypes). Overall, our results provide empirical support for the existence of a within-population-level polymorphism in genes underlying migratory timing in a diurnally migrating raptor.

A general theory of avian migratory connectivity

Birds exhibit a remarkable array of seasonal migrations. Despite much research describing migratory behaviour, the underlying forces driving how a species' breeding and wintering populations redistribute each year, that is, migratory connectivity, remain largely unknown. Here, we test the hypothesis that birds migrate in a way that minimises energy expenditure while considering intraspecific competition for energy acquisition, by developing a modelling framework that simulates an optimal redistribution of individuals between breeding and wintering areas. Using 25 species across the Americas, we find that the model accurately predicts empirical migration patterns, and thus offers a general explanation for migratory connectivity based on first ecological and energetic principles. Our model provides a strong basis for exploring additional processes underlying the ecology and evolution of migration, but also a framework for predicting how migration impacts local adaptation across seasons and how environmental change may affect population dynamics in migratory species.

An island-hopping bird reveals how founder events shape genome-wide divergence

When populations colonize new areas, both strong selection and strong drift can be experienced due to novel environments and small founding populations, respectively. Empirical studies have predominantly focused on the phenotype when assessing the role of selection, and limited neutral-loci when assessing founder-induced loss of diversity. Consequently, the extent to which processes interact to influence evolutionary trajectories is difficult to assess. Genomic-level approaches provide the opportunity to simultaneously consider these processes. Here, we examine the roles of selection and drift in shaping genomic diversity and divergence in historically documented sequential island colonizations by the silvereye (Zosterops lateralis). We provide the first empirical demonstration of the rapid appearance of highly diverged genomic regions following population founding, the position of which are highly idiosyncratic. As these regions rarely contained loci putatively under selection, it is most likely that these differences arise via the stochastic nature of the founding process. However, selection is required to explain rapid evolution of larger body size in insular silvereyes. Reconciling our genomic data with these phenotypic patterns suggests there may be many genomic routes to the island phenotype, which vary across populations. Finally, we show that accelerated divergence associated with multiple founding steps is the product of genome-wide rather than localized differences, and that diversity erodes due to loss of rare alleles. However, even multiple founder events do not result in divergence and diversity levels seen in evolutionary older subspecies, and therefore do not provide a shortcut to speciation as proposed by founder-effect speciation models.

A genoscape-network model for conservation prioritization in a migratory bird

Migratory animals are declining worldwide and coordinated conservation efforts are needed to reverse current trends. We devised a novel genoscape-network model that combines genetic analyses with species distribution modeling and demographic data to overcome challenges with conceptualizing alternative risk factors in migratory species across their full annual cycle. We applied our method to the long distance, Neotropical migratory bird, Wilson's Warbler (Cardellina pusilla). Despite a lack of data from some wintering locations, we demonstrated how the results can be used to help prioritize conservation of breeding and wintering areas. For example, we showed that when genetic, demographic, and network modeling results were considered together it became clear that conservation recommendations will differ depending on whether the goal is to preserve unique genetic lineages or the largest number of birds per unit area. More specifically, if preservation of genetic lineages is the goal, then limited resources should be focused on preserving habitat in the California Sierra, Basin Rockies, or Coastal California, where the 3 most vulnerable genetic lineages breed, or in western Mexico, where 2 of the 3 most vulnerable lineages overwinter. Alternatively, if preservation of the largest number of individuals per unit area is the goal, then limited conservation dollars should be placed in the Pacific Northwest or Central America, where densities are estimated to be the highest. Overall, our results demonstrated the utility of adopting a genetically based network model for integrating multiple types of data across vast geographic scales and better inform conservation decision-making for migratory animals.

Phylogenomic Data Reveal Widespread Introgression Across the Range of an Alpine and Arctic Specialist

Understanding how gene flow affects population divergence and speciation remains challenging. Differentiating one evolutionary process from another can be difficult because multiple processes can produce similar patterns, and more than one process can occur simultaneously. Although simple population models produce predictable results, how these processes balance in taxa with patchy distributions and complicated natural histories is less certain. These types of populations might be highly connected through migration (gene flow), but can experience stronger effects of genetic drift and inbreeding, or localized selection. Although different signals can be difficult to separate, the application of high-throughput sequence data can provide the resolution necessary to distinguish many of these processes. We present whole-genome sequence data for an avian species group with an alpine and arctic tundra distribution to examine the role that different population genetic processes have played in their evolutionary history. Rosy-finches inhabit high elevation mountaintop sky islands and high-latitude island and continental tundra. They exhibit extensive plumage variation coupled with low levels of genetic variation. Additionally, the number of species within the complex is debated, making them excellent for studying the forces involved in the process of diversification, as well as an important species group in which to investigate species boundaries. Total genomic variation suggests a broadly continuous pattern of allele frequency changes across the mainland taxa of this group in North America. However, phylogenomic analyses recover multiple distinct, well supported, groups that coincide with previously described morphological variation and current species-level taxonomy. Tests of introgression using D-statistics and approximate Bayesian computation reveal significant levels of introgression between multiple North American taxa. These results provide insight into the balance between divergent and homogenizing population genetic processes and highlight remaining challenges in interpreting conflict between different types of analytical approaches with whole-genome sequence data. [ABBA-BABA; approximate Bayesian computation; gene flow; phylogenomics; speciation; whole-genome sequencing.].

Rapid morphological divergence following a human-mediated introduction: the role of drift and directional selection

Theory predicts that when populations are established by few individuals, random founder effects can facilitate rapid phenotypic divergence even in the absence of selective processes. However, empirical evidence from historically documented colonisations suggest that, in most cases, drift alone is not sufficient to explain the rate of morphological divergence. Here, using the human-mediated introduction of the silvereye (Zosterops lateralis) to French Polynesia, which represents a potentially extreme example of population founding, we reassess the potential for morphological shifts to arise via drift alone. Despite only 80 years of separation from their New Zealand ancestors, French Polynesian silvereyes displayed significant changes in body and bill size and shape, most of which could be accounted for by drift, without the need to invoke selection. However, signatures of selection at genes previously identified as candidates for bill size and body shape differences in a range of bird species, also suggests a role for selective processes in driving morphological shifts within this population. Twenty-four SNPs in our RAD-Seq dataset were also found to be strongly associated with phenotypic variation. Hence, even under population founding extremes, when it is difficult to reject drift as the sole mechanism based on rate tests of phenotypic shifts, the additional role of divergent natural selection in novel environments can be revealed at the level of the genome.

Genomic divergence across ecological gradients in the Central African rainforest songbird (Andropadus virens)

The little greenbul, a common rainforest passerine from sub-Saharan Africa, has been the subject of long-term evolutionary studies to understand the mechanisms leading to rainforest speciation. Previous research found morphological and behavioural divergence across rainforest-savannah transition zones (ecotones), and a pattern of divergence with gene flow suggesting divergent natural selection has contributed to adaptive divergence and ecotones could be important areas for rainforests speciation. Recent advances in genomics and environmental modelling make it possible to examine patterns of genetic divergence in a more comprehensive fashion. To assess the extent to which natural selection may drive patterns of differentiation, here we investigate patterns of genomic differentiation among populations across environmental gradients and regions. We find compelling evidence that individuals form discrete genetic clusters corresponding to distinctive environmental characteristics and habitat types. Pairwise F_ST between populations in different habitats is significantly higher than within habitats, and this differentiation is greater than what is expected from geographic distance alone. Moreover, we identified 140 SNPs that showed extreme differentiation among populations through a genomewide selection scan. These outliers were significantly enriched in exonic and coding regions, suggesting their functional importance. Environmental association analysis of SNP variation indicates that several environmental variables, including temperature and elevation, play important roles in driving the pattern of genomic diversification. Results lend important new genomic evidence for environmental gradients being important in population differentiation.

Are Antarctic minke whales unusually abundant because of 20th century whaling?

Severe declines in megafauna worldwide illuminate the role of top predators in ecosystem structure. In the Antarctic, the Krill Surplus Hypothesis posits that the killing of more than 2 million large whales led to competitive release for smaller krill-eating species like the Antarctic minke whale. If true, the current size of the Antarctic minke whale population may be unusually high as an indirect result of whaling. Here, we estimate the long-term population size of the Antarctic minke whale prior to whaling by sequencing 11 nuclear genetic markers from 52 modern samples purchased in Japanese meat markets. We use coalescent simulations to explore the potential influence of population substructure and find that even though our samples are drawn from a limited geographic area, our estimate reflects ocean-wide genetic diversity. Using Bayesian estimates of the mutation rate and coalescent-based analyses of genetic diversity across loci, we calculate the long-term population size of the Antarctic minke whale to be 670,000 individuals (95% confidence interval: 374,000-1,150,000). Our estimate of long-term abundance is similar to, or greater than, contemporary abundance estimates, suggesting that managing Antarctic ecosystems under the assumption that Antarctic minke whales are unusually abundant is not warranted.

Not as the crow flies: a historical explanation for circuitous migration in Swainson's thrush (Catharus ustulatus)

Many migratory songbirds follow circuitous migratory routes instead of taking the shortest path between overwintering and breeding areas. Here, we study the migration patterns in Swainson's thrush (Catharus ustulatus), a neartic-neotropical migrant songbird, using molecular genetic approaches. This species is presently separated into genetically distinct coastal and continental populations that diverged during the Late Pleistocene (as indicated by molecular dating), yet appear to have retained ancestral patterns of migration. Low nucleotide diversity, a star-like haplotype phylogeny and unimodal mismatch distributions all support the hypothesis that both the coastal and the continental populations have undergone recent demographic expansions. Nearctic-neotropical banding and genetic data show nearly complete segregation of migratory routes and of overwintering locations: coastal populations migrate along the Pacific Coast to overwintering sites in Central America and Mexico, whereas continental populations migrate along an eastern route to overwintering sites in Panama and South America. Nearctic-neotropical banding data also show that continental birds north, northwest and east of this migratory divide fly thousands of miles east before turning south. We conclude that circuitous migration in the Swainson's thrush is an artefact of a Late Pleistocene range expansion.