Increasing the impact of vertebrate scientific collections through 3D imaging: The openVertebrate (oVert) Thematic Collections Network.

The impact of preserved museum specimens is transforming and increasing by three-dimensional (3D) imaging that creates high-fidelity online digital specimens. Through examples from the openVertebrate (oVert) Thematic Collections Network, we describe how we created a digitization community dedicated to the shared vision of making 3D data of specimens available and the impact of these data on a broad audience of scientists, students, teachers, artists, and more. High-fidelity digital 3D models allow people from multiple communities to simultaneously access and use scientific specimens. Based on our multiyear, multi-institution project, we identify significant technological and social hurdles that remain for fully realizing the potential impact of digital 3D specimens.

A highly contiguous genome assembly for the Yellow Warbler (Setophaga petechia).

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.

An elevational shift facilitated the Mesoamerican diversification of Azure-hooded Jays (Cyanolyca cucullata) during the Great American Biotic Interchange.

The Great American Biotic Interchange (GABI) was a key biogeographic event in the history of the Americas. The rising of the Panamanian land bridge ended the isolation of South America and ushered in a period of dispersal, mass extinction, and new community assemblages, which sparked competition, adaptation, and speciation. Diversification across many bird groups, and the elevational zonation of others, ties back to events triggered by the GABI. But the exact timing of these events is still being revealed, with recent studies suggesting a much earlier time window for faunal exchange, perhaps as early as 20 million years ago (Mya). Using a time-calibrated phylogenetic tree, we show that the jay genus Cyanolyca is emblematic of bird dispersal trends, with an early, pre-land bridge dispersal from Mesoamerica to South America 6.3-7.3 Mya, followed by a back-colonization of C. cucullata to Mesoamerica 2.3-4.8 Mya, likely after the land bridge was complete. As Cyanolyca species came into contact in Mesoamerica, they avoided competition due to a prior shift to lower elevation in the ancestor of C. cucullata. This shift allowed C. cucullata to integrate itself into the Mesoamerican highland avifauna, which our time-calibrated phylogeny suggests was already populated by higher-elevation, congeneric dwarf-jays (C. argentigula, C. pumilo, C. mirabilis, and C. nanus). The outcome of these events and fortuitous elevational zonation was that C. cucullata could continue colonizing new highland areas farther north during the Pleistocene. Resultingly, four C. cucullata lineages became isolated in allopatric, highland regions from Panama to Mexico, diverging in genetics, morphology, plumage, and vocalizations. At least two of these lineages are best described as species (C. mitrata and C. cucullata). Continued study will further document the influence of the GABI and help clarify how dispersal and vicariance shaped modern-day species assemblages in the Americas.

De novo assembly of a chromosome-level reference genome for the California Scrub-Jay, Aphelocoma californica.

We announce the assembly of the first de novo reference genome for the California Scrub-Jay (Aphelocoma californica). The genus Aphelocoma comprises four currently recognized species including many locally adapted populations across Mesoamerica and North America. Intensive study of Aphelocoma has revealed novel insights into the evolutionary mechanisms driving diversification in natural systems. Additional insights into the evolutionary history of this group will require continued development of high-quality, publicly available genomic resources. We extracted high molecular weight genomic DNA from a female California Scrub-Jay from northern California and generated PacBio HiFi long-read data and Omni-C chromatin conformation capture data. We used these data to generate a de novo partially phased diploid genome assembly, consisting of two pseudo-haplotypes, and scaffolded them using inferred physical proximity information from the Omni-C data. The more complete pseudo-haplotype assembly (arbitrarily designated "Haplotype 1") is 1.35 Gb in total length, highly contiguous (contig N50 = 11.53 Mb), and highly complete (BUSCO completeness score = 97%), with comparable scaffold sizes to chromosome-level avian reference genomes (scaffold N50 = 66.14 Mb). Our California Scrub-Jay assembly is highly syntenic with the New Caledonian Crow reference genome despite ~10 million years of divergence, highlighting the temporal stability of the avian genome. This high-quality reference genome represents a leap forward in publicly available genomic resources for Aphelocoma, and the family Corvidae more broadly. Future work using Aphelocoma as a model for understanding the evolutionary forces generating and maintaining biodiversity across phylogenetic scales can now benefit from a highly contiguous, in-group reference genome.

A highly contiguous reference genome for the Steller's jay (Cyanocitta stelleri).

The Steller's jay is a familiar bird of western forests from Alaska south to Nicaragua. Here, we report a draft reference assembly for the species generated from PacBio HiFi long-read and Omni-C chromatin-proximity sequencing data as part of the California Conservation Genomics Project (CCGP). Sequenced reads were assembled into 352 scaffolds totaling 1.16 Gb in length. Assembly metrics indicate a highly contiguous and complete assembly with a contig N50 of 7.8 Mb, scaffold N50 of 25.8 Mb, and BUSCO completeness score of 97.2%. Repetitive elements span 16.6% of the genome including nearly 90% of the W chromosome. Compared with high-quality assemblies from other members of the family Corvidae, the Steller's jay genome contains a larger proportion of repetitive elements than 4 crow species (Corvus), but a lower proportion of repetitive elements than the California scrub-jay (Aphelocoma californica). This reference genome will serve as an essential resource for future studies on speciation, local adaptation, phylogeography, and conservation genetics in this species of significant biological interest.

Temporal stability of the hybrid zone between Calocitta magpie-jays revealed through comparison of museum specimens and iNaturalist photos.

Hybrid zones are natural experiments for the study of avian evolution. Hybrid zones can be dynamic, moving as species adjust to new climates and habitats, with unknown implications for species and speciation. There are relatively few studies that have comparable modern and historic sampling to assess change in hybrid zone location and width over time, and those studies have generally found mixed results, with many hybrid zones showing change over time, but others showing stability. The white-throated magpie-jay (Calocitta formosa) and black-throated magpie-jay (Calocitta colliei) occur along the western coast of Mexico and Central America. The two species differ markedly in throat color and tail length, and prior observation suggests a narrow hybrid zone in southern Jalisco where individuals have mixed throat color. This study aims to assess the existence and temporal stability of this putative hybrid zone by comparing throat color between georeferenced historical museum specimens and modern photos from iNaturalist with precise locality information. Our results confirm the existence of a narrow hybrid zone in Jalisco, with modern throat scores gradually increasing from the parental ends of the cline toward the cline center in a sigmoidal curve characteristic of hybrid zones. Our temporal comparison suggests that the hybrid zone has not shifted its position between historical (pre-1973) and modern (post-2005) time periods-a surprising result given the grand scale of habitat change to the western Mexican lowlands during this time. An anomalous pocket of white-throated individuals in the northern range of the black-throated magpie-jay hints at the possibility of prehistorical long-distance introduction. Future genomic data will help disentangle the evolutionary history of these lineages and better characterize how secondary contact is affecting both the DNA and the phenotype of these species.

Hybrid zone or hybrid lineage: a genomic reevaluation of Sibley's classic species conundrum in Pipilo towhees.

Hybrid zones can be studied by modeling clines of trait variation (e.g., morphology, genetics) over a linear transect. Yet, hybrid zones can also be spatially complex, can shift over time, and can even lead to the formation of hybrid lineages with the right combination of dispersal and vicariance. We reassessed Sibley's (1950) gradient between Collared Towhee (Pipilo ocai) and Spotted Towhee (Pipilo maculatus) in Central Mexico to test whether it conformed to a typical tension-zone cline model. By comparing historical and modern data, we found that cline centers for genetic and phenotypic traits have not shifted over the course of 70 years. This equilibrium suggests that secondary contact between these species, which originally diverged over 2 million years ago, likely dates to the Pleistocene. Given the amount of mtDNA divergence, parental ends of the cline have very low autosomal nuclear differentiation (FST = 0.12). Dramatic and coincident cline shifts in mtDNA and throat color suggest the possibility of sexual selection as a factor in differential introgression, while a contrasting cline shift in green back color hints at a role for natural selection. Supporting the idea of a continuum between clinal variation and hybrid lineage formation, the towhee gradient can be analyzed as one population under isolation-by-distance, as a two-population cline, and as three lineages experiencing divergence with gene flow. In the middle of the gradient, a hybrid lineage has become partly isolated, likely due to forested habitat shrinking and fragmenting as it moved upslope after the last glacial maximum and a stark environmental transition. This towhee system offers a window into the potential outcomes of hybridization across a dynamic landscape including the creation of novel genomic and phenotypic combinations and incipient hybrid lineages.

Historical specimens and the limits of subspecies phylogenomics in the New World quails (Odontophoridae).

As phylogenomics focuses on comprehensive taxon sampling at the species and population/subspecies levels, incorporating genomic data from historical specimens has become increasingly common. While historical samples can fill critical gaps in our understanding of the evolutionary history of diverse groups, they also introduce additional sources of phylogenomic uncertainty, making it difficult to discern novel evolutionary relationships from artifacts caused by sample quality issues. These problems highlight the need for improved strategies to disentangle artifactual patterns from true biological signal as historical specimens become more prevalent in phylogenomic datasets. Here, we tested the limits of historical specimen-driven phylogenomics to resolve subspecies-level relationships within a highly polytypic family, the New World quails (Odontophoridae), using thousands of ultraconserved elements (UCEs). We found that relationships at and above the species-level were well-resolved and highly supported across all analyses, with the exception of discordant relationships within the two most polytypic genera which included many historical specimens. We examined the causes of discordance and found that inferring phylogenies from subsets of taxa resolved the disagreements, suggesting that analyzing subclades can help remove artifactual causes of discordance in datasets that include historical samples. At the subspecies-level, we found well-resolved geographic structure within the two most polytypic genera, including the most polytypic species in this family, Northern Bobwhites (Colinus virginianus), demonstrating that variable sites within UCEs are capable of resolving phylogenetic structure below the species level. Our results highlight the importance of complete taxonomic sampling for resolving relationships among polytypic species, often through the inclusion of historical specimens, and we propose an integrative strategy for understanding and addressing the uncertainty that historical samples sometimes introduce to phylogenetic analyses.

Combining Species Delimitation, Species Trees, and Tests for Gene Flow Clarifies Complex Speciation in Scrub-Jays.

Complex speciation, involving rapid divergence and multiple bouts of post-divergence gene flow, can obfuscate phylogenetic relationships and species limits. In North America, cases of complex speciation are common, due at least in part to the cyclical Pleistocene glacial history of the continent. Scrub-Jays in the genus Aphelocoma provide a useful case study in complex speciation because their range throughout North America is structured by phylogeographic barriers with multiple cases of secondary contact between divergent lineages. Here, we show that a comprehensive approach to genomic reconstruction of evolutionary history, i.e., synthesizing results from species delimitation, species tree reconstruction, demographic model testing, and tests for gene flow, is capable of clarifying evolutionary history despite complex speciation. We find concordant evidence across all statistical approaches for the distinctiveness of an endemic southern Mexico lineage (A. w. sumichrasti), culminating in support for the species status of this lineage under any commonly applied species concept. We also find novel genomic evidence for the species status of a Texas endemic lineage A. w. texana, for which equivocal species delimitation results were clarified by demographic modeling and spatially explicit models of gene flow. Finally, we find that complex signatures of both ancient and modern gene flow between the non-sister California Scrub-Jay (Aphelocoma californica) and Woodhouse's Scrub-Jay (Aphelocoma woodhouseii) result in discordant gene trees throughout the species' genomes despite clear support for their overall isolation and species status. In sum, we find that a multi-faceted approach to genomic analysis can increase our understanding of complex speciation histories, even in well-studied groups. Given the emerging recognition that complex speciation is relatively commonplace, the comprehensive framework that we demonstrate for interrogation of species limits and evolutionary history using genomic data can provide a necessary roadmap for disentangling the impacts of gene flow and incomplete lineage sorting (ILS) to better understand the systematics of other groups with similarly complex evolutionary histories. [Gene flow; genomics; introgression; phylogenomics; population genetics; species delimitation; species trees; systematics.].

Cryptic diversity across the Trans-Mexican Volcanic Belt of Mexico in the montane bunchgrass lizard Sceloporus subniger/ (Squamata: Phrynosomatidae).

Sceloporus subniger Poglaygen Smith is a montane bunchgrass lizard distributed across pine-oak forests of central Mexico. Prompted by the discovery of a new population of this lizard in far western Mexico, and by recent studies suggesting S. subniger may be a composite of several distinct species, we examined in more detail the genetic structure of S. subniger. We generated a mitochondrial DNA (mtDNA) dataset from 81 specimens and an ultraconserved elements (UCE) dataset representing thousands of genomic regions from 12 specimens to specifically evaluate the genetic distinctiveness of populations from western Michoacán and adjacent Jalisco along with the newly discovered population in the Sierra de Mascota in western Jalisco. We also recorded morphological data from 47 museum specimens to compare to our genetic data. Results from our analyses of the genetic data, augmented by specimen measurements and scale counts, support the notion that S. subniger is indeed a composite of distinct species. Montane bunchgrass lizards from western Michoacán and adjacent Jalisco, and from the Sierra de Mascota in western Jalisco, each represent distinct new species, which we describe and name here.

A rapid and cost-effective pipeline for digitization of museum specimens with 3D photogrammetry.

Natural history collections are yielding more information as digitization brings specimen data to researchers, connects specimens across museums, and as new technologies allow for more large-scale data collection. Therefore, a key goal in specimen digitization is developing methods that both increase access and allow for the highest yield of phenomic data. 3D digitization is increasingly popular because it has the potential to meet both aspects of that key goal. However, current methods overlook or do not prioritize some of the most sought-after phenotypic traits, those involving the external appearance of specimens, especially color. Here, we introduce an efficient and cost-effective pipeline for 3D photogrammetry to capture the external appearance of natural history specimens and other museum objects. 3D photogrammetry aligns and compares sets of dozens, hundreds, or even thousands of photos to create 3D models. The hardware set-up requires little physical space and around $3,000 in initial investment, while the software pipeline requires $1,400/year for proprietary software subscriptions (with open-source alternatives). The creation of each 3D model takes 1-2 hours/specimen and much of the software pipeline is automated with minimal supervision required, including the onerous step of mesh processing. We showcase the method by creating 3D models for most of the type specimens in the Moore Laboratory of Zoology bird collection and show that digital bill measurements are comparable to hand-taken measurements. Color data, while not included as part of this pipeline, is easily extractable from the models and one of the most promising areas of data collection. Future advances can adapt the method for ultraviolet reflectance capture and increased efficiency and model quality. Combined with genomic data, phenomic data from 3D models including photogrammetry will open new doors to understanding organismal evolution.

More than skin and bones: Comparing extraction methods and alternative sources of DNA from avian museum specimens.

Next-generation sequencing has greatly expanded the utility and value of museum collections by revealing specimens as genomic resources. As the field of museum genomics grows, so does the need for extraction methods that maximize DNA yields. For avian museum specimens, the established method of extracting DNA from toe pads works well for most specimens. However, for some specimens, especially those of birds that are very small or very large, toe pads can be a poor source of DNA. In this study, we apply two DNA extraction methods (phenol-chloroform and silica column) to three different sources of DNA (toe pad, skin punch and bone) from 10 historical avian museum specimens. We show that a modified phenol-chloroform protocol yielded significantly more DNA than a silica column protocol (e.g., Qiagen DNeasy Blood & Tissue Kit) across all tissue types. However, extractions using the silica column protocol contained longer fragments on average than those using the phenol-chloroform protocol, probably as a result of loss of small fragments through the silica column. While toe pads yielded more DNA than skin punches and bone fragments, skin punches proved to be a reliable alternative source of DNA and might be especially appealing when toe pad extractions are impractical. Overall, we found that historical bird museum specimens contain substantial amounts of DNA for genomic studies under most extraction scenarios, but that a phenol-chloroform protocol consistently provides the high quantities of DNA required for most current genomic protocols.

Phylogenomic species delimitation in microendemic frogs of the Brazilian Atlantic Forest.

The advent of next-generation sequencing allows researchers to use large-scale datasets for species delimitation analyses, yet one can envision an inflection point where the added accuracy of including more loci does not offset the increased computational burden. One alternative to including all loci could be to prioritize the analysis of loci for which there is an expectation of high informativeness. Here, we explore the issue of species delimitation and locus selection with montane species from two anuran genera that have been isolated in sky islands across the southern Brazilian Atlantic Forest: Melanophryniscus (Bufonidae) and Brachycephalus (Brachycephalidae). To delimit species, we obtained genetic data using target enrichment of ultraconserved elements from 32 populations (13 for Melanophryniscus and 19 for Brachycephalus), and we were able to create datasets that included over 800 loci with no missing data. We ranked loci according to their number of parsimony-informative sites, and we performed species delimitation analyses using BPP with the most informative 10, 20, 40, 80, 160, 320, and 640 loci. We identified three types of phylogenetic node: nodes with either consistently high or low support regardless of the number of loci or their informativeness and nodes that were initially poorly supported where support became stronger as we included more data. When viewed across all sensitivity analyses, our results suggest that the current species richness in both genera is likely underestimated. In addition, our results show the effects of different sampling strategies on species delimitation using phylogenomic datasets.

Comparison of ultraconserved elements (UCEs) to microsatellite markers for the study of avian hybrid zones: a test in Aphelocoma jays.

OBJECTIVE: Hybrid zones are geographic regions where genetically distinct taxa interbreed, resulting in offspring of mixed ancestry. California Scrub-Jays (Aphelocoma californica) and Woodhouse's Scrub-Jays (A. woodhouseii) come into secondary contact and hybridize in western Nevada. Although previous work investigated divergence and gene flow between these species using a handful of microsatellite markers, the hybrid zone has not been studied using genome-scale markers, providing an opportunity to assess genome-wide introgression, test for a genetic basis for ecomorphological traits, and compare these estimates to those derived from microsatellites. RESULTS: Using variant sites flanking ultraconserved elements (UCEs), we performed population assignment and quantified hybrid ancestry for 16 individuals across the zone of secondary contact. Our study included 2468 SNPs distributed throughout the genome, allowing discrimination of genetic affinities of hybrid individuals that were similar to estimates from microsatellites. We show a relationship between bill and wing length and the genetic composition of individuals that was not found in prior work using microsatellites, suggesting a genetic basis for these traits. Our analyses demonstrate the utility of UCEs for the analysis of hybrid zones and provide a basis for future studies to identify the genomic architecture of speciation and phenotypic differences between these incipient species.

Dense Geographic and Genomic Sampling Reveals Paraphyly and a Cryptic Lineage in a Classic Sibling Species Complex.

Incomplete or geographically biased sampling poses significant problems for research in phylogeography, population genetics, phylogenetics, and species delimitation. Despite the power of using genome-wide genetic markers in systematics and related fields, approaches such as the multispecies coalescent remain unable to easily account for unsampled lineages. The Empidonax difficilis/Empidonax occidentalis complex of small tyrannid flycatchers (Aves: Tyrannidae) is a classic example of widely distributed species with limited phenotypic geographic variation that was broken into two largely cryptic (or "sibling") lineages following extensive study. Though the group is well-characterized north of the US Mexico border, the evolutionary distinctiveness and phylogenetic relationships of southern populations remain obscure. In this article, we use dense genomic and geographic sampling across the majority of the range of the E. difficilis/E. occidentalis complex to assess whether current taxonomy and species limits reflect underlying evolutionary patterns, or whether they are an artifact of historically biased or incomplete sampling. We find that additional samples from Mexico render the widely recognized species-level lineage E. occidentalis paraphyletic, though it retains support in the best-fit species delimitation model from clustering analyses. We further identify a highly divergent unrecognized lineage in a previously unsampled portion of the group's range, which a cline analysis suggests is more reproductively isolated than the currently recognized species E. difficilis and E. occidentalis. Our phylogeny supports a southern origin of these taxa. Our results highlight the pervasive impacts of biased geographic sampling, even in well-studied vertebrate groups like birds, and illustrate what is a common problem when attempting to define species in the face of recent divergence and reticulate evolution.

Museum genomics reveals the speciation history of Dendrortyx wood-partridges in the Mesoamerican highlands.

Natural history collections are increasingly valued as genomic resources. Their specimens reflect the combined efforts of collectors and curators over hundreds of years. For many rare or endangered species, specimens are the only readily available source of DNA. We leveraged specimens from a historical collection to study the evolutionary history of wood-partridges in the genus Dendrortyx. The three Dendrortyx species are found in the highlands of central Mexico and Central America south to Costa Rica. One of these species is endangered, and in general, Dendrortyx are secretive and poorly represented in tissue collections. We extracted DNA from historical museum specimens and sequenced ultraconserved elements (UCEs) and mitochondrial DNA (mtDNA) to assess their phylogeny and divergence times. Phylogenies built from hundreds to thousands of nuclear markers were well resolved and largely congruent with an mtDNA phylogeny. The divergence times revealed an unusually old avian divergence across the Isthmus of Tehuantepec in the Pliocene around 3.6 million years ago. Combined with other recent studies, our results challenge the general pattern that highland bird divergences in Mesoamerica are relatively young and influenced by the Pleistocene glacial cycles compared to the older divergences of reptiles and plants, which are thought to overlap more with periods of mountain formation. We also found evidence for monophyletic genetic lineages in mountain ranges within the widespread D. macroura, which should be investigated further with integrative taxonomic methods. Our study demonstrates the power of museum genomics to provide insight into the evolutionary histories of groups where modern samples are lacking.

Cryptic diversity in the Mexican highlands: Thousands of UCE loci help illuminate phylogenetic relationships, species limits and divergence times of montane rattlesnakes (Viperidae: Crotalus).

With the continued adoption of genome-scale data in evolutionary biology comes the challenge of adequately harnessing the information to make accurate phylogenetic inferences. Coalescent-based methods of species tree inference have become common, and concatenation has been shown in simulation to perform well, particularly when levels of incomplete lineage sorting are low. However, simulation conditions are often overly simplistic, leaving empiricists with uncertainty regarding analytical tools. We use a large ultraconserved element data set (>3,000 loci) from rattlesnakes of the Crotalus triseriatus group to delimit lineages and estimate species trees using concatenation and several coalescent-based methods. Unpartitioned and partitioned maximum likelihood and Bayesian analysis of the concatenated matrix yield a topology identical to coalescent analysis of a subset of the data in bpp. ASTRAL analysis on a subset of the more variable loci also results in a tree consistent with concatenation and bpp, whereas the SVDquartets phylogeny differs at additional nodes. The size of the concatenated matrix has a strong effect on species tree inference using SVDquartets, warranting additional investigation on optimal data characteristics for this method. Species delimitation analyses suggest up to 16 unique lineages may be present within the C. triseriatus group, with divergences occurring during the Neogene and Quaternary. Network analyses suggest hybridization within the group is relatively rare. Altogether, our results reaffirm the Mexican highlands as a biodiversity hotspot and suggest that coalescent-based species tree inference on data subsets can provide a strongly supported species tree consistent with concatenation of all loci with a large amount of missing data.

Combining ultraconserved elements and mtDNA data to uncover lineage diversity in a Mexican highland frog (Sarcohyla; Hylidae).

Molecular studies have uncovered significant diversity in the Mexican Highlands, leading to the description of many new endemic species. DNA approaches to this kind of species discovery have included both mitochondrial DNA (mtDNA) sequencing and multilocus genomic methods. While these marker types have often been pitted against one another, there are benefits to deploying them together, as linked mtDNA data can provide the bridge between uncovering lineages through rigorous multilocus genomic analysis and identifying lineages through comparison to existing mtDNA databases. Here, we apply one class of multilocus genomic marker, ultraconserved elements (UCEs), and linked mtDNA data to a species complex of frogs (Sarcohyla bistincta, Hylidae) found in the Mexican Highlands. We generated data from 1,891 UCEs, which contained 1,742 informative SNPs for S. bistincta and closely related species and captured mitochondrial genomes for most samples. Genetic analyses based on both whole loci and SNPs agree there are six to seven distinct lineages within what is currently described as S. bistincta. Phylogenies from UCEs and mtDNA mostly agreed in their topologies, and the few differences suggested a more complex evolutionary history of the mtDNA marker. Our study demonstrates that the Mexican Highlands still hold substantial undescribed diversity, making their conservation a particularly urgent goal. The Trans-Mexican Volcanic Range stands out as a significant geographic feature in Sarcohyla and may have acted as a dispersal corridor for S. bistincta to spread to the north. Combining multilocus genomic data with linked mtDNA data is a useful approach for identifying potential new species and associating them with already described taxa, which will be especially important in groups with undescribed subadult phenotypes and cryptic species.

Genome-wide signals of drift and local adaptation during rapid lineage divergence in a songbird.

The formation of independent evolutionary lineages involves neutral and selective factors, and understanding their relative roles in population divergence is a fundamental goal of speciation research. Correlations between allele frequencies and environmental variability can reveal the role of selection, yet the relative contribution of drift can be difficult to establish. Recently diversified taxa like the Oregon junco (Aves, Passerellidae, Junco hyemalis oreganus) of western North America provide ideal scenarios to apply genetic-environment association analyses (GEA) while controlling for population structure. Analysis of genome-wide SNP loci revealed marked genetic structure consisting of differentiated populations in isolated, dry southern mountain ranges, and less divergent, recently expanded populations in humid northern latitudes. We used correlations between genomic and environmental variance to test for three specific modes of evolutionary divergence: (a) drift in geographic isolation, (b) differentiation along continuous selective gradients and (c) isolation-by-adaptation. We found evidence of strong drift in southern mountains, but also signals of local adaptation driven by temperature, precipitation, elevation and vegetation, especially when controlling for population history. We identified numerous variants under selection scattered across the genome, suggesting that local adaptation can promote rapid differentiation when acting over multiple independent loci.

The gut of the finch: uniqueness of the gut microbiome of the Galápagos vampire finch.

BACKGROUND: Darwin's finches are a clade of 19 species of passerine birds native to the Galápagos Islands, whose biogeography, specialized beak morphologies, and dietary choices-ranging from seeds to blood-make them a classic example of adaptive radiation. While these iconic birds have been intensely studied, the composition of their gut microbiome and the factors influencing it, including host species, diet, and biogeography, has not yet been explored. RESULTS: We characterized the microbial community associated with 12 species of Darwin's finches using high-throughput 16S rRNA sequencing of fecal samples from 114 individuals across nine islands, including the unusual blood-feeding vampire finch (Geospiza septentrionalis) from Darwin and Wolf Islands. The phylum-level core gut microbiome for Darwin's finches included the Firmicutes, Gammaproteobacteria, and Actinobacteria, with members of the Bacteroidetes at conspicuously low abundance. The gut microbiome was surprisingly well conserved across the diversity of finch species, with one exception-the vampire finch-which harbored bacteria that were either absent or extremely rare in other finches, including Fusobacterium, Cetobacterium, Ureaplasma, Mucispirillum, Campylobacter, and various members of the Clostridia-bacteria known from the guts of carnivorous birds and reptiles. Complementary stable isotope analysis of feathers revealed exceptionally high δ15N isotope values in the vampire finch, resembling top marine predators. The Galápagos archipelago is also known for extreme wet and dry seasons, and we observed a significant seasonal shift in the gut microbial community of five additional finch species sampled during both seasons. CONCLUSIONS: This study demonstrates the overall conservatism of the finch gut microbiome over short (< 1 Ma) divergence timescales, except in the most extreme case of dietary specialization, and elevates the evolutionary importance of seasonal shifts in driving not only species adaptation, but also gut microbiome composition.