Sexually selected differences in warbler plumage are related to a putative inversion on the Z chromosome.

Large structural variants in the genome, such as inversions, may play an important role in producing population structure and local adaptation to the environment through suppression of recombination. However, relatively few studies have linked inversions to phenotypic traits that are sexually selected and may play a role in reproductive isolation. Here, we found that geographic differences in the sexually selected plumage of a warbler, the common yellowthroat (Geothlypis trichas), are largely due to differences in the Z (sex) chromosome (males are ZZ), which contains at least one putative inversion spanning 40% (31/77 Mb) of its length. The inversions on the Z chromosome vary dramatically east and west of the Appalachian Mountains, which provides evidence of cryptic population structure within the range of the most widespread eastern subspecies (G. t. trichas). In an eastern (New York) and western (Wisconsin) population of this subspecies, female prefer different male ornaments; larger black facial masks are preferred in Wisconsin and larger yellow breasts are preferred in New York. The putative inversion also contains genes related to vision, which could influence mating preferences. Thus, structural variants on the Z chromosome are associated with geographic differences in male ornaments and female choice, which may provide a mechanism for maintaining different patterns of sexual selection in spite of gene flow between populations of the same subspecies.

Broad-scale seasonal climate tracking is a consequence, not a driver, of avian migratory connectivity.

Tracking climatic conditions throughout the year is often assumed to be an adaptive behaviour underlying seasonal migration patterns in animal populations. We investigate this hypothesis using genetic markers data to map migratory connectivity for 27 genetically distinct bird populations from 7 species. We found that the variation in seasonal climate tracking across our suite of populations at a continental scale is more likely a consequence, rather than a direct driver, of migratory connectivity, which is primarily shaped by energy efficiency-i.e., optimizing the balance between accessing available resources and movement costs. However, our results also suggest that regional-scale seasonal precipitation tracking affects population migration destinations, thus revealing a potential scale dependency of ecological processes driving migration. Our results have implications for the conservation of these migratory species under climate change, as populations tracking climate seasonally are potentially at higher risk if they adapt to a narrow range of climatic conditions.

Biogeographic factors contributing to the diversification of Euphoniinae (Aves, Passeriformes, Fringillidae): a phylogenetic and ancestral areas analysis.

Factors such as the Andean uplift, Isthmus of Panama, and climate changes have influenced bird diversity in the Neotropical region. Studying bird species that are widespread in Neotropical highlands and lowlands can help us understand the impact of these factors on taxa diversification. Our main objectives were to determine the biogeographic factors that contributed to the diversification of Euphoniinae and re-evaluate their phylogenetic relationships. The nextRAD and mitochondrial data were utilized to construct phylogenies. The ancestral distribution range was then estimated using a time-calibrated phylogeny, current species ranges, and neotropical regionalization. The phylogenies revealed two main Euphoniinae clades, Chlorophonia and Euphonia, similar to previous findings. Furthermore, each genus has distinctive subclades corresponding to morphology and geography. The biogeographic results suggest that the Andean uplift and the establishment of the western Amazon drove the vicariance of Chlorophonia and Euphonia during the Miocene. The Chlorophonia lineage originated in the Andes mountains and spread to Central America and the Mesoamerican highlands after the formation of the Isthmus of Panama. Meanwhile, the ancestral area of Euphonia was the Amazonas, from which it spread to trans-Andean areas during the Pliocene and Pleistocene due to the separation of the west lowlands from Amazonas due to the Northern Andean uplift. Chlorophonia and Euphonia species migrated to the Atlantic Forest during the Pleistocene through corridors from the East Andean Humid Forest and Amazonas. These two genera had Caribbean invasions with distinct geographic origins and ages. Finally, we suggested taxonomic changes in the genus Euphonia based on the study's phylogenetic, morphological, and biogeographic findings.

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.

Genomic analyses reveal poaching hotspots and illegal trade in pangolins from Africa to Asia.

The white-bellied pangolin (Phataginus tricuspis) is the world's most trafficked mammal and is at risk of extinction. Reducing the illegal wildlife trade requires an understanding of its origins. Using a genomic approach for tracing confiscations and analyzing 111 samples collected from known geographic localities in Africa and 643 seized scales from Asia between 2012 and 2018, we found that poaching pressures shifted over time from West to Central Africa. Recently, Cameroon's southern border has emerged as a site of intense poaching. Using data from seizures representing nearly 1 million African pangolins, we identified Nigeria as one important hub for trafficking, where scales are amassed and transshipped to markets in Asia. This origin-to-destination approach offers new opportunities to disrupt the illegal wildlife trade and to guide anti-trafficking measures.

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.

Genetic identification of avian samples recovered from solar energy installations.

Renewable energy production and development will drastically affect how we meet global energy demands, while simultaneously reducing the impact of climate change. Although the possible effects of renewable energy production (mainly from solar- and wind-energy facilities) on wildlife have been explored, knowledge gaps still exist, and collecting data from wildlife remains (when negative interactions occur) at energy installations can act as a first step regarding the study of species and communities interacting with facilities. In the case of avian species, samples can be collected relatively easily (as compared to other sampling methods), but may only be able to be identified when morphological characteristics are diagnostic for a species. Therefore, many samples that appear as partial remains, or "feather spots"-known to be of avian origin but not readily assignable to species via morphology-may remain unidentified, reducing the efficiency of sample collection and the accuracy of patterns observed. To obtain data from these samples and ensure their identification and inclusion in subsequent analyses, we applied, for the first time, a DNA barcoding approach that uses mitochondrial genetic data to identify unknown avian samples collected at solar facilities to species. We also verified and compared identifications obtained by our genetic method to traditional morphological identifications using a blind test, and discuss discrepancies observed. Our results suggest that this genetic tool can be used to verify, correct, and supplement identifications made in the field and can produce data that allow accurate comparisons of avian interactions across facilities, locations, or technology types. We recommend implementing this genetic approach to ensure that unknown samples collected are efficiently identified and contribute to a better understanding of wildlife impacts at renewable energy projects.

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.

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.

Strong habitat-specific phenotypic plasticity but no genome-wide differentiation across a rainforest gradient in an African butterfly.

Habitat-specific thermal responses are well documented in various organisms and likely determine the vulnerability of populations to climate change. However, the underlying roles of genetics and plasticity that shape such habitat-specific patterns are rarely investigated together. Here we examined the thermal plasticity of the butterfly Bicyclus dorothea originating from rainforest and ecotone habitats in Cameroon under common garden conditions. We also sampled wild-caught butterflies from forest and ecotone sites and used RADseq to explore genome-wide population differentiation. We found differences in the level of phenotypic plasticity across habitats. Specifically, ecotone populations exhibited greater sensitivity in wing eyespot features with variable development temperatures relative to rainforest populations. Known adaptive roles of wing eyespots in Bicyclus species suggest that this morphological plasticity is likely under divergent selection across environmental gradients. However, we found no distinct population structure of genome-wide variation between habitats, suggesting high level of ongoing gene flow between habitats is homogenizing most parts of the genome.

Widespread gene flow following range expansion in Anna's Hummingbird.

Anthropogenic changes have altered the historical distributions of many North American taxa. As environments shift, ecological and evolutionary processes can combine in complex ways to either stimulate or inhibit range expansion. Here, we examined the role of evolution in a rapid range expansion whose ecological context has been well-documented, Anna's Hummingbird (Calypte anna). Previous studies have suggested that the C. anna range expansion is the result of an ecological release facilitated by human-mediated environmental changes, where access to new food sources have allowed further filling of the abiotic niche. We examined the role of gene flow and adaptation during range expansion from their native California breeding range, north into Canada and east into New Mexico and Texas, USA. Using low coverage whole genome sequencing we found high genetic diversity, low divergence, and little evidence of selection on the northern and eastern expansion fronts. Additionally, there are no clear barriers to gene flow across the native and expanded range. The lack of selective signals between core and expanded ranges could reflect (i) an absence of novel selection pressure in the expanded range (supporting the ecological release hypothesis), (ii) swamping of adaptive variation due to high gene flow, or (iii) limitations of genome scans for detecting small shifts in allele frequencies across many loci. Nevertheless, our results provide an example where strong selection is not apparent during a rapid, contemporary range shift.

Winter connectivity and leapfrog migration in a migratory passerine.

Technological advances in migratory tracking tools have revealed a remarkable diversity in migratory patterns. One such pattern is leapfrog migration, where individuals that breed further north migrate to locations further south. Here, we analyzed migration patterns in the Painted Bunting (Passerina ciris) using a genetic-based approach. We started by mapping patterns of genetic variation across geographic space (called a genoscape) using 386 individuals from 25 populations across the breeding range. We then genotyped an additional 230 samples from 31 migration stopover locations and 178 samples from 16 wintering locations to map patterns of migratory connectivity. Our analyses of genetic variation across the breeding range show the existence of four genetically distinct groups within the species: Eastern, Southwestern, Louisiana, and Central groups. Subsequent assignment of migrating and wintering birds to genetic groups illustrated that birds from the Central group migrated during the fall via western Mexico or southern Texas, spent the winter from northeastern Mexico to Panama, and migrated north via the Gulf Coast of Mexico. While Louisiana birds overlapped with Central birds on their spring migratory routes along the Gulf Coast, we found that Louisiana birds had a more restricted wintering distribution in the Yucatan Peninsula and Central America. Further estimation of the straight-line distance from the predicted breeding location to the wintering location revealed that individuals sampled at lower winter latitudes traveled to greater distances (i.e., the predicted breeding area was further north; p > .001), confirming that these species exhibit a leapfrog migration pattern. Overall, these results demonstrate the utility of a genoscape-based approach for identifying range-wide patterns of migratory connectivity such as leapfrog migration with a high degree of clarity.

Evolutionary potential mitigates extinction risk under climate change in the endangered southwestern willow flycatcher.

The complexity of global anthropogenic change makes forecasting species responses and planning effective conservation actions challenging. Additionally, important components of a species' adaptive capacity, such as evolutionary potential, are often not included in quantitative risk assessments due to lack of data. While genomic proxies for evolutionary potential in at-risk species are increasingly available, they have not yet been included in extinction risk assessments at a species-wide scale. In this study, we used an individual-based, spatially explicit, dynamic eco-evolutionary simulation model to evaluate the extinction risk of an endangered desert songbird, the southwestern willow flycatcher (Empidonax traillii extimus), in response to climate change. Using data from long-term demographic and habitat studies in conjunction with genome-wide ecological genomics research, we parameterized simulations that include 418 sites across the breeding range, genomic data from 225 individuals, and climate change forecasts spanning 3 generalized circulation models and 3 emissions scenarios. We evaluated how evolutionary potential, and the lack of it, impacted population trajectories in response to climate change. We then investigated the compounding impact of drought and warming temperatures on extinction risk through the mechanism of increased nest failure. Finally, we evaluated how rapid action to reverse greenhouse gas emissions would influence population responses and species extinction risk. Our results illustrate the value of incorporating evolutionary, demographic, and dispersal processes in a spatially explicit framework to more comprehensively evaluate the extinction risk of threatened and endangered species and conservation actions to promote their recovery.

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.

Genetic evidence for widespread population size expansion in North American boreal birds prior to the Last Glacial Maximum.

Pleistocene climate cycles are well documented to have shaped contemporary species distributions and genetic diversity. Northward range expansions in response to deglaciation following the Last Glacial Maximum (LGM; approximately 21 000 years ago) are surmised to have led to population size expansions in terrestrial taxa and changes in seasonal migratory behaviour. Recent findings, however, suggest that some northern temperate populations may have been more stable than expected through the LGM. We modelled the demographic history of 19 co-distributed boreal-breeding North American bird species from full mitochondrial gene sets and species-specific molecular rates. We used these demographic reconstructions to test how species with different migratory strategies were affected by glacial cycles. Our results suggest that effective population sizes increased in response to Pleistocene deglaciation earlier than the LGM, whereas genetic diversity was maintained throughout the LGM despite shifts in geographical range. We conclude that glacial cycles prior to the LGM have most strongly shaped contemporary genetic diversity in these species. We did not find a relationship between historic population dynamics and migratory strategy, contributing to growing evidence that major switches in migratory strategy during the LGM are unnecessary to explain contemporary migratory patterns.

Genetic and ecological drivers of molt in a migratory bird.

The ability of animals to sync the timing and location of molting (the replacement of hair, skin, exoskeletons or feathers) with peaks in resource availability has important implications for their ecology and evolution. In migratory birds, the timing and location of pre-migratory feather molting, a period when feathers are shed and replaced with newer, more aerodynamic feathers, can vary within and between species. While hypotheses to explain the evolution of intraspecific variation in the timing and location of molt have been proposed, little is known about the genetic basis of this trait or the specific environmental drivers that may result in natural selection for distinct molting phenotypes. Here we take advantage of intraspecific variation in the timing and location of molt in the iconic songbird, the Painted Bunting (Passerina ciris) to investigate the genetic and ecological drivers of distinct molting phenotypes. Specifically, we use genome-wide genetic sequencing in combination with stable isotope analysis to determine population genetic structure and molting phenotype across thirteen breeding sites. We then use genome-wide association analysis (GWAS) to identify a suite of genes associated with molting and pair this with gene-environment association analysis (GEA) to investigate potential environmental drivers of genetic variation in this trait. Associations between genetic variation in molt-linked genes and the environment are further tested via targeted SNP genotyping in 25 additional breeding populations across the range. Together, our integrative analysis suggests that molting is in part regulated by genes linked to feather development and structure (GLI2 and CSPG4) and that genetic variation in these genes is associated with seasonal variation in precipitation and aridity. Overall, this work provides important insights into the genetic basis and potential selective forces behind phenotypic variation in what is arguably one of the most important fitness-linked traits in a migratory bird.

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.

Estimating global genetic diversity loss.

A mathematical framework may help inform conservation efforts.

Limited domestic introgression in a final refuge of the wild pigeon.

Domesticated animals have been culturally and economically important throughout history. Many of their ancestral lineages are extinct or genetically endangered following hybridization with domesticated relatives. Consequently, they have been understudied compared to the ancestral lineages of domestic plants. The domestic pigeon Columba livia, which was pivotal in Darwin's studies, has maintained outsized cultural significance. Its role as a model organism spans the fields of behavior, genetics, and evolution. Domestic pigeons have hybridized with their progenitor, the Rock Dove, rendering the latter of dubious genetic status. Here, we use genomic and morphological data from the putative Rock Doves of the British Isles to identify relictual undomesticated populations. We reveal that Outer Hebridean Rock Doves have experienced minimal levels of introgression. Our results outline the contemporary status of these wild pigeons, highlighting the role of hybridization in the homogenization of genetic lineages.